Water and carbon fluxes above European coniferous forests modelled with artificial neural networks

Artificial neural networks are used to select a minimal set of input variables to model water vapour and carbon exchange of coniferous forest ecosystems, independently of tree species and without detailed physiological information. Neural networks are used because of their power to fit highly non-linear relations between input and output-variables. Radiation, temperature, vapour pressure deficit and time of the day showed to be the dynamic input variables that determine ecosystem water fluxes. The same variables, together with projected leaf area index are needed for modelling CO₂-fluxes. The results for the individual sites show that the neural networks found mean water and carbon flux responses to the driving variables valid for all sites. The sensitivity analysis of the derived neural networks shows that the LAI-effect of the CO₂-flux model is overfitted because of the low variability of LAI. However, the predictions of CO₂-fluxes of sites not included in the calibration set indicate that the LAI-response of the network is reliable and that results can be used as a first estimate of the net ecosystem carbon exchange of the forest sites. Independent predictions of forest ecosystem vapour fluxes were equally satisfying as empirical models specifically calibrated for the individual sites. The results indicate that both short term water and carbon fluxes of European coniferous forests can be modelled without using detailed physiological and site specific information.     

Modelling radiative balance in a row-crop canopy: Row-soil surface net radiation partition

Row crops like vineyards undergo various and significant manipulations of training system and cultural practices, which strongly influence the quality of products. Variations of water vapour fluxes from the soil surface and the leaves in the row volume are closely linked to the ratio of energy available to each compartment. A physically realistic model of available energy partition between the rows and the soil surface is therefore a key factor towards optimization of such systems, and must be included in canopy models. A number of available models were not directly validated. The purpose of the study was therefore to design a model of net radiation partition and check it directly.The model of net radiation partition between rows (Rnv), considered as a whole, and intervening soil surface (Rns) of a row-crop canopy was developed from physically realistic yet simple assumptions:

-

global solar (short wave) radiation partition was calculated by a previously validated geometric model;

-

long-wave radiative fluxes between the soil surface, the rows and the atmosphere were calculated from the corresponding view factors, which only depended on canopy geometry;

-

atmospheric radiation was estimated by a simple empirical relation based on air temperature as the only input variable;

-

air temperature in the vicinity of leaves replaced leaf surface temperatures as a more convenient input variable, with little loss of information.

The input variables were incoming direct and diffuse solar radiation, soil surface mean temperature and air temperature near the leaves. The main parameters were soil and leaf albedos, row porosity and dimensions.A direct validation of the model was attempted by measuring net radiation above the canopy and at five positions above the soil surface in a vineyard of the Bordeaux area. The reliability of soil surface net radiation measurements was estimated by thorough error propagation analysis. When found significant, errors were corrected and finally soil surface net radiation data were corrected only for delay in direct downward solar radiation striking net radiometers, because canopy was discontinuous and height of net radiometers was not negligible compared to canopy height.In these conditions, model calculations were in agreement with measurements, although the model slightly underestimated Rns and therefore overestimated Rnv. As the mean error was about 20 W m⁻², and therefore compatible with instrument accuracy, the results were considered satisfactory.This available energy partition model is able to estimate radiative balance in various canopy systems and in various thermal environment conditions, leading to easier simulations of energy balance and water fluxes. It could therefore be a useful tool for optimizing row-crop canopies, taking fully into account any kind of present or future thermal environment.     

Mass conservation and atmospheric dynamics in the Regional Atmospheric Modeling System (RAMS)

This paper examines the spatio-temporal patterns of atmospheric carbon dioxide transport predicted by the Regional Atmospheric Modeling System (RAMS). Forty-eight hour simulations over northern New England incorporating a simple representation of the diurnal summertime surface carbon dioxide forcing arising from biological activity indicate that, in its native formulation, RAMS exhibits a significant degree of mass non-conservation. Domain-wide rates of non-physical mass gain and mass loss are as large as three percent per day which translates into approximately eleven parts per million per day for carbon dioxide — enough to rapidly dilute the signature of carbon dioxide fluxes arising from biological activity. Analysis shows that this is due to the approximation used by RAMS to compute the Exner function. Substitution of the exact, physically complete equation improves mass conservation by two orders of magnitude. In addition to greatly improving mass conservation, use of the complete Exner function equation has a substantial impact on the spatial pattern of carbon dioxide predicted by the model, yielding predictions differing from a conventional RAMS simulation by as much as forty parts per million. Such differences have important implications both for comparisons of modeled atmospheric carbon dioxide concentrations to observations and for carbon dioxide inversion studies, which use estimates of atmospheric transport of carbon dioxide in conjunction with measurements of atmospheric carbon dioxide concentrations to infer the spatio-temporal distribution of surface carbon dioxide fluxes. Furthermore, use of the complete Exner function equation affects the vertical velocity and water mixing ratio fields, causing significant changes in accumulated precipitation over the region.     

A decomposition analysis of CO2 emissions: evidence from Malaysia’s tourism industry

The most important question raised from issues of environmental degradation is how economic activities bring about changes that will result in pollution. In the pursuit of tourism economy, contrary to popular interest, the travel and tourism (T&T) industry may cause environmental damages through the emissions of carbon dioxide (CO₂) from energy consumption in areas such as transportation and delivery of amenities. Given this major concern, this paper attempts to investigate the linkage between tourism and CO₂ emissions in Malaysia between 1981 and 2011. In particular, this study fills the knowledge gap by taking a closer look at the impact of international tourist arrivals on CO₂ emissions by sector – electricity and heat generation and transport. Results from the bound test method suggest that there exists a long-run relationship among the variables under consideration when CO₂ emissions become the dependent variable. The original result is similarly robust to alternatives, which are CO₂ emissions from sectors of electricity and heat generation and transport. Furthermore, the vector error correction model causality analysis indicates a causal relationship between tourism and CO₂ emissions by transport and electricity and heat generation. Subsequently, several tourism-related policies are drawn from these findings.     

Acknowledging the spatial heterogeneity in modelling/reconstructing carbon dioxide exchange in a northern aapa mire

Chamber method is commonly used to measure the CO₂ exchange from plant communities. Due to low time resolution, actual measurements reflect only momentary CO₂ exchange rates. Therefore, a common way to derive seasonal or annual estimates is to establish models describing the response of CO₂ exchange to environmental variables, and then to reconstruct the CO₂ exchange over the desired time period. There are several alternative ways to obtain the CO₂ balance for the entire mire: models can be parameterized by individual sample plots, plant communities or the entire site. Similarly, the CO₂ balance can be reconstructed by plots, plant communities or the entire site. We tested how the choice of the modelling and reconstruction approach influences the CO₂ exchange estimates for the entire mire and for individual sample plots and plant communities. We measured the CO₂ exchange in a spatially heterogeneous sedge-dominated northern aapa mire for two growing seasons. We observed high spatial variation in CO₂ balance between the plant communities. We noticed that when the CO₂ balances of individual sample plots or plant communities are of interest, using a model appropriate for the entire site may result in biased estimates. In worst case the different modelling approaches may turn the CO₂ balance of an individual sample plot from positive to negative. Further, while using the whole ecosystem approach in modelling, the superior ability of chamber method in acknowledging spatial variation is lost. While the modelled growing season CO₂ balance of the mire ranged from 232 to 625 g CO₂ m⁻² depending on the chosen modelling and reconstruction approach, the average estimates still remained within the uncertainty range of one another. Acknowledgement of the spatial variation in plant community level makes the areal estimate more robust to varying weather conditions. Further, the reliability of estimates is improved by explicit formulation of the choices behind the modelling and reconstruction units reflecting the study objectives.     

The performance of RAMS in representing the convective boundary layer structure in a very steep valley

Data from a comprehensive field study in the Riviera Valley of Southern Switzerland are used to investigate convective boundary layer structure in a steep valley and to evaluate wind and temperature fields, convective boundary layer height, and surface sensible heat fluxes as predicted by the mesoscale model RAMS. Current parameterizations of surface and boundary layer processes in RAMS, as well as in other mesoscale models, are based on scaling laws strictly valid only for flat topography and uniform land cover. Model evaluation is required to investigate whether this limits the applicability of RAMS in steep, inhomogeneous terrain. One clear-sky day with light synoptic winds is selected from the field study. Observed temperature structure across and along the valley is nearly homogeneous while wind structure is complex with a wind speed maximum on one side of the valley. Upvalley flows are not purely thermally driven and mechanical effects near the valley entrance also affect the wind structure. RAMS captured many of the observed boundary layer characteristics within the steep valley. The wind field, temperature structure, and convective boundary layer height in the valley are qualitatively simulated by RAMS, but the horizontal temperature structure across and along the valley is less homogeneous in the model than in the observations. The model reproduced the observed net radiation, except around sunset and sunrise when RAMS does not take into account the shadows cast by the surrounding topography. The observed sensible heat fluxes fall within the range of simulated values at grid points surrounding the measurement sites. Some of the scatter between observed and simulated turbulent sensible heat fluxes are due to sub-grid scale effects related to local topography.     

The dynamic causal links between CO2 emissions from transport,real GDP,energy use and international tourism

This study examines the dynamic causality relationship between international tourism and carbon dioxide (CO₂) emissions from transport, real gross domestic product and energy use. The vector error correction model and Granger causality test approach have been used to investigate these relationships for the top ten international tourism destinations spanning the period 1995–2013. Results reveal a unidirectional causality running from CO₂ emissions to economic growth without feedback; a bidirectional causality between economic growth and energy use; a bidirectional causality between international tourism and economic growth; and a bidirectional causality between international tourism and energy use. They also suggest that energy use and international tourism both contribute to the decrease of emissions level coming from transport sector, while economic growth leads to the increase of CO₂ emissions. This study can be used in policy recommendations by encouraging countries to use clean energy and to stimulate tourism sector for combating global warming.     

Modeling the carbon dynamics of the La Grande hydroelectric complex in northern Quebec

A mechanistic semi-empirical carbon cycle model of the La Grande reservoir complex in northern Quebec, Canada was conceived in order to investigate the climate impact of such a large alteration of the continental water cycle. The model includes inputs from the drainage basin, organic matter release from flooded soils, CO₂ emissions across the water-atmosphere interface and sedimentation. Most input data stems from previous research by our group on those ecosystems. The model includes the seven reservoirs of the La Grande complex and was run for periods of 50 and 100 years. Terrigeneous dissolved, particulate and suspended soil carbon fluxes and concentrations were computed. Over 100 years, 31.3 × 10¹² g C are released from flooded soils, equivalent to 28-29% of inputs from the drainage basin. 40-74% of dissolved organic carbon is mineralized. CO₂ fluxes over 100 years are 50.5-79.8 × 10¹² g C, 46.4-67.9 × 10¹² g C more than in the absence of reservoirs. The increase in mineralization of organic matter and in CO₂ emissions is a result of the increase in cumulated water residence time due to the creation of the reservoirs. Changes in other carbon sinks and sources likely offset a part of this additional carbon flux to the atmosphere. In the first years following flooding of the reservoir, organic carbon release from flooded soils exceeds CO₂ emissions, implying the downstream export of large quantities of eroded soil organic carbon. After this initial period, CO₂ emissions are fuelled by organic carbon originating from the drainage basin.     

Moderate seawater acidification does not elicit long-term metabolic depression in the blue mussel Mytilus edulis

Marine organisms are exposed to increasingly acidic oceans, as a result of equilibration of surface ocean water with rising atmospheric CO₂ concentrations. In this study, we examined the physiological response of Mytilus edulis from the Baltic Sea, grown for 2 months at 4 seawater pCO₂ levels (39, 113, 243 and 405 Pa/385, 1,120, 2,400 and 4,000 μatm). Shell and somatic growth, calcification, oxygen consumption and \textNH₄ ⁺ {\text{NH}}_{4}^{ + } excretion rates were measured in order to test the hypothesis whether exposure to elevated seawater pCO₂ is causally related to metabolic depression. During the experimental period, mussel shell mass and shell-free dry mass (SFDM) increased at least by a factor of two and three, respectively. However, shell length and shell mass growth decreased linearly with increasing pCO₂ by 6–20 and 10–34%, while SFDM growth was not significantly affected by hypercapnia. We observed a parabolic change in routine metabolic rates with increasing pCO₂ and the highest rates (+60%) at 243 Pa. \textNH₄ ⁺ {\text{NH}}_{4}^{ + } excretion rose linearly with increasing pCO₂. Decreased O:N ratios at the highest seawater pCO₂ indicate enhanced protein metabolism which may contribute to intracellular pH regulation. We suggest that reduced shell growth under severe acidification is not caused by (global) metabolic depression but is potentially due to synergistic effects of increased cellular energy demand and nitrogen loss.     

Long-run dynamics of renewable energy consumption on carbon emissions and economic growth in the European union

This paper examines long-run and short-run dynamics of renewable energy consumption on carbon dioxide (CO₂) emissions and economic growth in the European Union. This study employs cointegration tests, Granger causality tests and vector error correction estimates to examine the direction of Granger causality, the long-run dynamics of economic growth and energy variables on carbon emissions. This study analyses time series data from the World Development Indicators over the period from1961 to 2012. The results of this study support a link between renewable energy consumption, economic growth, industrialization, exports and CO₂ emissions in the long-run and short-run. The results support that the sign of the long-run dynamics from the endogenous variables to the CO₂ emissions variable is negative and significant, which implies that the energy and environmental policies of the European Union aimed at curbing CO₂ emissions must have been effective in the long-term. Furthermore, renewable energy consumption and exports have significant negative impact on CO₂ emissions in the short-run. However, industrialization and economic growth have positive impact on CO₂ emissions in the short-run. The results suggest that both economic growth and industrialization must have been achieved at the cost of harming the environment. The finding suggests that the increasing consumption of renewable energy tends to play an important role in curbing carbon emissions in the region.     

Estimating marine biogeochemical rates of the carbonate pH system—A Kalman filter tested

Oxygen (O₂), nitrate (NO₃), dissolved inorganic carbon (DIC) or pCO₂, and pH or total alkalinity (TA), are useful indices of marine chemical, physical and biological processes operating on varying time-scales. Although these properties are increasingly being monitored at high frequency, they have not been extensively used for studying ecosystem dynamics. We test whether we can estimate time-evolving biogeochemical rates (e.g. primary production, respiration, calcification and carbonate dissolution, and nitrification) from synthetic high frequency time-series of O₂, NO₃, DIC, pCO₂, TA or pH. More specifically, a Kalman filter has been implemented in a very simplified biogeochemical model describing the dynamics of O₂, NO₃, DIC and TA and linking the concentration data to biogeochemical fluxes. Different sets of concentration data are assimilated and biogeochemical rates are estimated. The frequency of assimilation required to get acceptable results is investigated and is compared with the frequency of sampling in the field or in controlled experimental settings.Smoothing of the data to remove data noise before assimilation improves the estimation of the biogeochemical rates. The best estimated rates are obtained when assimilating O₂, NO₃ and TA although the assimilation of DIC instead of TA also gives satisfactory results. In case pH or pCO₂ is assimilated rather than DIC or TA, the linearization of the (now nonlinear) observation equation introduces perturbations and the Kalman filter behaves suboptimal. We conclude that, given the resolution of data required, the tool has potential to estimate biogeochemical rates of the carbonate system under controlled settings.     

Assessment of canopy stomatal conductance models using flux measurements

Stomatal conductance (g) is a key parameter in controlling energy and water exchanges between canopy and the atmosphere. Stomatal conductance models proposed by Ball, Woodrow and Berry (BWB) and Leuning have been increasingly used in land surface schemes. In a recent study, a new diagnostic index was developed by Wang et al. to examine the response of g to humidity and new models were proposed to resolve problems identified in the BWB and Leuning models. This approach is theoretically sound, but relies on canopy latent heat and CO₂ fluxes and environmental variables at the leaf surface which are not available at most eddy correlation (EC) observation sites. In this study, we tested the diagnostic index by empirically correcting EC measurements to canopy-level fluxes and by replacing leaf surface variables by their corresponding ambient air variables, and re-examined the stomatal conductance models of BWB, Leuning, and Wang et al. We found that the impact of the above modifications on the evaluation of g–humidity relationships is very small. This study provides a practical approach to investigate the stomatal response to humidity using routine EC measurements.     

Introducing effects of temperature and CO2 elevation on tree growth into a statistical growth and yield model

Impacts of elevated temperature and CO₂ on tree growth were introduced into a statistical growth and yield model for Finnish conditions based on corresponding predictions obtained from a physiological growth model. This one-way link between models was made by means of species-specific transfer functions describing the increase in stem volume growth of trees as a function of elevated temperature and CO₂, stand density and the tree's competition status in a stand of Scots pine (Pinus sylvestris), silver birch (Betula pendula) and Norway spruce (Picea abies). This method allows the inner dynamics of the statistical model to be followed when the impacts of temperature and CO₂ elevation on tree growth are introduced into the calculation of volume growth and further allocated between diameter and height growth. In this way compatibility with previous predictions of tree growth by means of statistical models and related model systems under current climatic conditions could be retained.The performance of the statistical model with species-specific transfer functions was evaluated by comparing its predictions with corresponding predictions given by a physiological model under conditions of elevated temperature and CO₂. These calculations revealed that the growth response of individual trees to elevated temperature and CO₂ can be introduced into the statistical model from a physiological growth model with an outcome that results in fairly satisfactory growth responses at the stand level as well.     

Spatial trends in abundance of long-finned pilot whales, white-beaked dolphins and harbour porpoises in West Greenland

Data from an aerial line transect survey conducted off West Greenland during August–September 2007 were used to estimate the abundance of long-finned pilot whales (Globicephala melas), white-beaked dolphins (Lagenorhynchus albirostris) and harbour porpoises (Phocoena phocoena). The abundance of each species was estimated using mark-recapture distance sampling techniques to correct for perception bias, and correction factors for time spent at the surface were applied. The fully corrected abundance estimates were 8,133 long-finned pilot whales, 11,984 white-beaked dolphins and 33,271 harbour porpoises. Based on density surface modelling methods, a count model with a generalised additive model formulation was used to relate abundance to spatial variables. Response curves indicated that the preferred habitats were deep offshore areas in Midwest Greenland for pilot whales, deep water over steep seabed slopes in South Greenland for white-beaked dolphins and relatively shallow inshore waters in Midwest–South Greenland for harbour porpoises. The abundance estimates and spatial trends for the three species are the first obtained from Greenland.     

Experimental ecology of the temperate scleractinian coral Astrangia danae I. Partition of respiration,photosynthesis and calcification between host and symbionts

Calcification, photosynthesis and respiration of the scleractinian coral Astrangia danae were calculated from the changes in total alkalinity, pH, calculated total CO₂, and oxygen concentration produced by colonies incubated in glass jars. A correction for changes in ammonia, nitrate and nitrite was taken into account and the method evaluated. The fluxes of oxygen and CO₂ were highly correlated (r=0.99). The statistical error of alkalinity determinations was less than 10% of the changes observed in the slowest calcifying samples. Metabolism of polyparium alone was estimated by difference after removal of tissue and reincubation of bare corallum. Zooxanthellae concentration in the polyps was obtained from cell counts made on homogenates of polyp tissue. The calculated photosynthetic rate of the zooxanthellae in vivo was 25 mol O₂ (10⁸ cell)^-1 h^-1 at a light intensity of 120 Ein m^-2 s^-1. In corals having 0.5x10⁹ zooxanthellae/dm² of colony area up to 8% of the total photosynthesis was attributed to the corallum microcosm. Polyp respiration, photosynthesis, and CaCO₃ uptake rates were all much higher than rates previously reported from A. danae, apparently because in these experiments the organisms were better fed. This increased photosynthesis in turn enhanced calcification still further. The symbiosis therefore appears to provide a growth advantage even to fed corals, under the conditions of these experiments.     

Photosynthetic carbon acquisition in the red alga Gracilaria conferta

In this continuing study on photosynthesis of the marine red alga Gracilaria conferta, it was found that ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) in crude extracts had a K _m (CO₂) of 85 M. Since seawater contains only ca. 10 M CO₂, it appears that this alga must possess a CO₂ concetrating system in order to supply sufficient CO₂ to the vicinity of the enzyme. Because this species is a C₃ plant (and thus lacks the C₄ system for concentrating CO₂), but can utilize HCO₃ ^- as an exogenous carbon source, we examined whether HCO₃ ^- uptake could be the initial step of such a CO₂ concetrating system. The surface pH of G. conferta thalli was 9.4 during photosynthesis. At this pH, estimated maximal uncatalyzed HCO₃ ^- dehydration (CO₂ formation) within the unstirred layer was too slow to account for measured phostosynthetic rates, even in the presence of an external carbonic anhydrase inhibitor. This observation, and the marked pH increase in the unstirred layer following the onset of light, suggests that a HCO₃ ^- transport system (probably coupled to transmembrane H⁺/OH^- fluxes) operates at the plasmalemma level. The involvement of surface-bound carbonic anhydrase in such a system remains, however, obscure. The apparent need of marine macroalgae such as G. conferta for CO₂ concentrating mechanisms is discussed with regard to their low affinity of Rubisco to CO₂ and the low rate of CO₂ supply in water. The close similarity between rates of Rubisco carboxylation and measured photosynthesis further suggests that the carboxylase activity, rather than inorganic carbon transport and intercoversion events, could be an internal limiting factor for photosynthetic rates of G. conferta.     

Fernmeßverfahren

The combination of remote sensing methods like Doppler lidar and FTIR allows remote determination of mass fluxes of gases. Doppler lidar measures the three-dimensional wind vector in the vicinity of diffuse sources, or the velocity of air in a chimney plume if an industrial complex is monitored. FTIR is a multicomponent remote sensing method for gas concentrations. The Fourier transformation of an interferogram of a Michelson interferometer within a FTIR system converts the recorded intensity (function of optical path length) to a spectral signal (function of wavenumber). Both information, velocity and concentration, give the mass fluxes of the tracer (gas). A first test was performed at Munich-Nord power station with FTIR and cw-Doppler lidar. Fluxes of CO₂, CO, NO, and HCl were determined. The results are in good agreement with the fluxes measured byin situ instruments of the power station. The method can be used to control industrial complexes from an outside observation site.     

Effects of pCO2 on physiology and skeletal mineralogy in a tidal pool coralline alga Corallina elongata

Marine organisms inhabiting environments where pCO₂/pH varies naturally are suggested to be relatively resilient to future ocean acidification. To test this hypothesis, the effect of elevated pCO₂ was investigated in the articulated coralline red alga Corallina elongata from an intertidal rock pool on the north coast of Brittany (France), where pCO₂ naturally varied daily between 70 and 1000 μatm. Metabolism was measured on algae in the laboratory after they had been grown for 3 weeks at pCO₂ concentrations of 380, 550, 750 and 1000 μatm. Net and gross primary production, respiration and calcification rates were assessed by measurements of oxygen and total alkalinity fluxes using incubation chambers in the light and dark. Calcite mol % Mg/Ca (mMg/Ca) was analysed in the tips, branches and basal parts of the fronds, as well as in new skeletal structures produced by the algae in the different pCO₂ treatments. Respiration, gross primary production and calcification in light and dark were not significantly affected by increased pCO₂. Algae grown under elevated pCO₂ (550, 750 and 1000 μatm) formed fewer new structures and produced calcite with a lower mMg/Ca ratio relative to those grown under 380 μatm. This study supports the assumption that C. elongata from a tidal pool, where pCO₂ fluctuates over diel and seasonal cycles, is relatively robust to elevated pCO₂ compared to other recently investigated coralline algae.     

Responses of marine benthic microalgae to elevated CO2

Increasing anthropogenic CO₂ emissions to the atmosphere are causing a rise in pCO₂ concentrations in the ocean surface and lowering pH. To predict the effects of these changes, we need to improve our understanding of the responses of marine primary producers since these drive biogeochemical cycles and profoundly affect the structure and function of benthic habitats. The effects of increasing CO₂ levels on the colonisation of artificial substrata by microalgal assemblages (periphyton) were examined across a CO₂ gradient off the volcanic island of Vulcano (NE Sicily). We show that periphyton communities altered significantly as CO₂ concentrations increased. CO₂ enrichment caused significant increases in chlorophyll a concentrations and in diatom abundance although we did not detect any changes in cyanobacteria. SEM analysis revealed major shifts in diatom assemblage composition as CO₂ levels increased. The responses of benthic microalgae to rising anthropogenic CO₂ emissions are likely to have significant ecological ramifications for coastal systems.     

Effects of boreal forest management practices on the climate impact of CO2 emissions from bioenergy

In Life Cycle Assessment (LCA), carbon dioxide (CO₂) emissions from biomass combustion are traditionally assumed climate neutral if the bioenergy system is CO₂ flux neutral, i.e. the quantity of CO₂ released approximately equals the amount of CO₂ sequestered in biomass. This convention is a plausible assumption for fast growing biomass species, but is inappropriate for slower growing biomass, like forests. In this case, the climate impact from biomass combustion can be potentially underestimated if CO₂ emissions are ignored, or overestimated, if biogenic CO₂ is considered equal to anthropogenic CO₂. The estimation of the effective climate impact should take into account how the CO₂ fluxes are distributed over time: the emission of CO₂ from bioenergy approximately occurs at a single point in time, while the absorption by the new trees is spread over several decades. Our research target is to include this dynamic time dimension in unit-based impact analysis, using a boreal forest stand as case study. The boreal forest growth is modelled with an appropriate function, and is investigated under different forestry regimes (affecting the growth rate and the year of harvest). Specific atmospheric decay functions for biomass-derived CO₂ are then elaborated for selected combinations of forest management options. The contribution to global warming is finally quantified using the GWP_bio index as climate metric. Results estimates the effects of these practices on the characterization factor used for the global warming potential of CO₂ from bioenergy, and point out the key role played by the selected time horizon.