Climate and vegetation controls on boreal zone energy exchange

Large, sustained soil CO2 efflux but rapid recovery of CH4 oxidation in post-harvest and post-fire stands in a mixedwood boreal forest

The net effect of forest disturbances, such as fires and harvesting, on soil greenhouse gas fluxes is determined by their impacts on both biological and physical factors, as well as the temporal dynamics of these effects post-disturbance. Although harvesting and fire may have distinct effects on soil carbon (C) dynamics, the temporal patterns in soil CO2 and CH4 fluxes and the potential differences between types of disturbances, remain poorly characterized in boreal forests. In this study, we measured soil CO2 and CH4 fluxes using a off-axis integrated cavity output spectroscopy system in snow-free seasons over two years in post-harvest and post-fire chronosequence sites within a mixedwood boreal forest in northwestern Ontario, Canada. Soil CO2 efflux showed a post-disturbance peak, with differing dynamics depending on the disturbance type: post-harvest stands exhibited a nearly tenfold increase (from ∼1 to ∼11 μmol CO2.m-2.s-1) from 1 to 9-10 years post-disturbance, followed by a steep decline; post-fire stands showed a more gradual increase, peaking at ∼6-7.2 μmol CO2.m-2.s-1 after ∼12-15 years. The youngest post-harvest stands were net sources of CH4,whereas post-fire stands were never net CH4 sources. In both disturbance types, the strength of the CH4 sink increased with stand age, approaching ∼2.4 nmol.m-2.s-1 by 15 years post-disturbance. Volumetric water content, bulk density, litter depth, and pH were significant predictors of CO2 fluxes; for CH4 fluxes, litter depth, pH, and the interaction of VWC and soil temperature were significant predictors in both disturbance types, with EC also showing a relationship in post-harvest stands. Our findings indicate that while soil CH4 oxidation rapidly recovers following disturbance, both post-harvest and post-fire stands show a multi-decade release of soil CO2 that is too large to be offset by C gains over this period.

Revealing the widespread potential of forests to increase low level cloud cover

Forests play a key role in humanity's current challenge to mitigate climate change thanks to their capacity to sequester carbon. Preserving and expanding forest cover is considered essential to enhance this carbon sink. However, changing the forest cover can further affect the climate system through biophysical effects. One such effect that is seldom studied is how afforestation can alter the cloud regime, which can potentially have repercussions on the hydrological cycle, the surface radiation budget and on planetary albedo itself. Here we provide a global scale assessment of this effect derived from satellite remote sensing observations. We show that for 67% of sampled areas across the world, afforestation would increase low level cloud cover, which should have a cooling effect on the planet. We further reveal a dependency of this effect on forest type, notably in Europe where needleleaf forests generate more clouds than broadleaf forests.

How tree species, tree size, and topographical location influenced tree transpiration in northern boreal forests during the historic 2018 drought

Trees in northern latitude ecosystems are projected to experience increasing drought stress as a result of rising air temperatures and changes in precipitation patterns in northern latitude ecosystems. However, most drought-related studies on high-latitude boreal forests (>50°N) have been conducted in North America, with few studies quantifying the response in European and Eurasian boreal forests. Here, we tested how daily whole-tree transpiration (Q, Liters day^-1 ) and Q normalized for mean daytime vapor pressure deficit (Q_DZ , Liters day^-1 kPa^-1 ) were affected by the historic 2018 drought in Europe. More specifically, we examined how tree species, size, and topographic position affected drought response in high-latitude mature boreal forest trees. We monitored 30 Pinus sylvestris (pine) and 30 Picea abies (spruce) trees distributed across a topographic gradient in northern Sweden. In general, pine showed a greater Q_DZ control compared to spruce during periods of severe drought (standardized precipitation-evapotranspiration index: SPEI < -1.5), suggesting that the latter are more sensitive to drought. Overall, Q_DZ reductions (using non-drought Q_DZ as reference) were less pronounced in larger trees during severe drought, but there was a species-specific pattern: Q_DZ reductions were greater in pine trees at high elevations and greater in spruce trees at lower elevations. Despite lower Q_DZ during severe drought, drought spells were interspersed with small precipitation events and overcast conditions, and Q_DZ returned to pre-drought conditions relatively quickly. This study highlights unique species-specific responses to drought, which are additionally driven by a codependent interaction among tree size, relative topographic position, and unique regional climate conditions.

Surface Energy Flux Estimation in Two Boreal Settings in Alaska Using a Thermal-Based Remote Sensing Model

Recent Arctic warming has led to changes in the hydrological cycle. Circum-Arctic and circumboreal ecosystems are showing evidence of “greening” and “browning” due to temperature warming leading to shrub encroachment, tree mortality and deciduousness. Increases in latent heat flux from increased evapotranspiration rates associated with deciduous-dominated ecosystems may be significant, because deciduous vegetation has extremely high-water use and water storage capacity compared to coniferous and herbaceous plant species. Thus, the impact of vegetation change in boreal ecosystems on regional surface energy balance is a significant knowledge gap that must be addressed to better understand observed trends in water use/availability and tree mortality. To this end, output from a two-source energy balance model (TSEB) with modifications for high latitude boreal ecosystems was evaluated using flux tower measurements and Terra/Aqua MODIS remote sensing data collected over the two largest boreal forest types in Alaska (birch and black spruce). Data under clear and overcast days and from leaf-out to senescence from 2012 to 2016 were used for validation. Using flux tower observations and local model inputs, modifications to the model formulation for soil heat flux, net radiation partitioning, and canopy transpiration were required for the boreal forest. These improvements resulted in a mean absolute percent difference of around 23% for turbulent daytime fluxes when surface temperature from the flux towers was used, similar to errors reported in other studies conducted in warmer climates. Results when surface temperature from Terra/Aqua MODIS estimates were used as model input suggested that these model improvements are pertinent for regional scale applications. Vegetation indices and LAI time-series from the Terra/Aqua MODIS products were confirmed to be appropriate for energy flux estimation in the boreal forest to describe vegetation properties (LAI and green fraction) when field observations are not available. Model improvements for boreal settings identified in this study will be implemented operationally over North America to map surface energy fluxes at regional scales using long time series of remote sensing estimates as part of NOAA’s GOES Evapotranspiration and Drought Information System.

Habitat patchiness, ecological connectivity and the uneven recovery of boreal stream ecosystems from an experimental drought

Ongoing climate change is increasing the occurrence and intensity of drought episodes worldwide, including in boreal regions not previously regarded as drought prone, and where the impacts of drought remain poorly understood. Ecological connectivity is one factor that might influence community structure and ecosystem functioning post-drought, by facilitating the recovery of sensitive species via dispersal at both local (e.g. a nearby habitat patch) and regional (from other systems within the same region) scales. In an outdoor mesocosm experiment, we investigated how impacts of drought on boreal stream ecosystems are altered by the spatial arrangement of local habitat patches within stream channels, and variation in ecological connectivity with a regional species pool. We measured basal ecosystem processes underlying carbon and nutrient cycling: (a) algal biomass accrual; (b) microbial respiration; and (c) decomposition of organic matter, and sampled communities of aquatic fungi and benthic invertebrates. An 8-day drought event had strong impacts on both community structure and ecosystem functioning, including algal accrual, leaf decomposition and microbial respiration, with many of these impacts persisting even after water levels had been restored for 3.5 weeks. Enhanced connectivity with the regional species pool and increased aggregation of habitat patches also affected multiple response variables, especially those associated with microbes, and in some cases reduced the effects of drought to a small extent. This indicates that spatial processes might play a role in the resilience of communities and ecosystem functioning, given enough time. These effects were however insufficient to facilitate significant recovery in algal growth before seasonal dieback began in autumn. The limited resilience of ecosystem functioning in our experiment suggests that even short-term droughts can have extended consequences for stream ecosystems in the world's vast boreal region, and especially on the ecosystem processes and services mediated by algal biofilms.

The Dynamics of Transpiration to Evapotranspiration Ratio under Wet and Dry Canopy Conditions in a Humid Boreal Forest

Humid boreal forests are unique environments characterized by a cold climate, abundant precipitation, and high evapotranspiration. Transpiration ( E T ), as a component of evapotranspiration (E), behaves differently under wet and dry canopy conditions, yet very few studies have focused on the dynamics of transpiration to evapotranspiration ratio ( E T / E ) under transient canopy wetness states. This study presents field measurements of E T / E at the Montmorency Forest, Québec, Canada: a balsam fir boreal forest that receives ∼ 1600 mm of precipitation annually (continental subarctic climate; Köppen classification subtype Dfc). Half-hourly observations of E and E T were obtained over two growing seasons using eddy-covariance and sap flow (Granier’s constant thermal dissipation) methods, respectively, under wet and dry canopy conditions. A series of calibration experiments were performed for sap flow, resulting in species-specific calibration coefficients that increased estimates of sap flux density by 34 % ± 8 % , compared to Granier’s original coefficients. The uncertainties associated with the scaling of sap flow measurements to stand E T , especially circumferential and spatial variations, were also quantified. From 30 wetting–drying events recorded during the measurement period in summer 2018, variations in E T / E were analyzed under different stages of canopy wetness. A combination of low evaporative demand and the presence of water on the canopy from the rainfall led to small E T / E . During two growing seasons, the average E T / E ranged from 35 % ± 2 % to 47 % ± 3 % . The change in total precipitation was not the main driver of seasonal E T / E variation, therefore it is important to analyze the impact of rainfall at half-hourly intervals.

Spatial Heterogeneity of the Carbon Emission Effect Resulting from Urban Expansion among Three Coastal Agglomerations in China

Land-use change, particularly urban expansion, can greatly affect the carbon balance, both from the aspects of terrestrial ecosystems and anthropogenic carbon emissions. Coastal China is a typical region of rapid urban expansion, and obvious spatial heterogeneity exists from the north to south. However, the different urban change characteristics and the effect on carbon balance remain undetermined. By unifying the spatial-temporal resolution of carbon source and sink data, we effectively compared the carbon budgets of three coastal urban agglomerations in China. The results show that all of the three urban agglomerations have undergone an obvious urban expansion process, with the built-up area increasing from 1.03 × 104 km2 in 2000 to 3.06 × 104 km2 in 2013. For Beijing–Tianjin–Hebei (BTH), the built-up area gradually expanded. The built-up area in the Yangtze River Delta (YRD) gradually changed before 2007 but rapidly grew thereafter. The built-up expansion of the Pearl River Delta (PRD) passed through three growing stages and showed the largest mean patch size. Carbon emission spatial patterns in the three urban agglomerations are consistent with their economic development, from which the net ecosystem production (NEP) spatial patterns are very different. Compared to carbon emissions, NEP has a carbon sink effect and can absorb some carbon emissions, but the amounts were all much lower than the carbon emissions in the three urban agglomerations. The carbon sink effect in the Yangtze River Delta is the most obvious, with the Pearl River Delta following, and the lowest effect is in Beijing–Tianjin–Hebei. Finally, a scientific basis for policy-making is provided for viable CO2 emission mitigation policies.

Vertical air temperature gradients under the shade of two contrasting urban tree species during different types of summer days

Moderation of thermal energy balance through the canopies of urban trees is well known. However, a more functional and quantitative view of the heterogeneous urban environment and their influence on the below-canopy vertical air temperature gradients is largely missing. Throughout the summer 2016 we continuously measured air temperature at three different heights (at 1.5, 3 and 4.5m from the ground) under the canopies of two common but contrasting street tree species in respect of eco-physiology and morphology in Munich, Germany: Robinia pseudoacacia L. (ring porous) and Tilia cordata Mill. (diffuse porous). Along with air and surface temperature we also measured meteorological and edaphic variables and categorized summer time as cool, mild, summer and hot days. Global radiation, vapour pressure deficit and soil temperature increased as the days got warmer but precipitation, soil moisture and wind speed showed the reversed pattern. Overall, T. cordata trees with higher leaf area index and sap-wood area provided three times more transpiration than R. pseudoacacia. On an average air temperature gradient of outside to inside canopy dropped from 1.8°C to 1.3°C for T. cordata but from 1.5°C to only 0.5°C for R. pseudoacacia as the days got warmer. Vertical decline of air cooling effect was around 1°C from canopy to the near-ground (1.5m). Lower soil moisture but higher soil temperature suggested that cool air from the canopy mixed with a higher amount of sensible heat flux under the canopies of T. cordata compared to the R. pseudoacacia as the days got warmer. The study indicated a threshold for extreme hot days when grass surface evapotranspirational cooling will not be as effective and act like built surfaces rather deep shading from tree canopies will be important.

Do the energy fluxes and surface conductance of boreal coniferous forests in Europe scale with leaf area?

Earth observing systems are now routinely used to infer leaf area index (LAI) given its significance in spatial aggregation of land surface fluxes. Whether LAI is an appropriate scaling parameter for daytime growing season energy budget, surface conductance (G_s ), water- and light-use efficiency and surface-atmosphere coupling of European boreal coniferous forests was explored using eddy-covariance (EC) energy and CO₂ fluxes. The observed scaling relations were then explained using a biophysical multilayer soil-vegetation-atmosphere transfer model as well as by a bulk G_s representation. The LAI variations significantly alter radiation regime, within-canopy microclimate, sink/source distributions of CO₂ , H₂ O and heat, and forest floor fluxes. The contribution of forest floor to ecosystem-scale energy exchange is shown to decrease asymptotically with increased LAI, as expected. Compared with other energy budget components, dry-canopy evapotranspiration (ET) was reasonably 'conservative' over the studied LAI range 0.5-7.0 m² m^-2 . Both ET and G_s experienced a minimum in the LAI range 1-2 m² m^-2 caused by opposing nonproportional response of stomatally controlled transpiration and 'free' forest floor evaporation to changes in canopy density. The young forests had strongest coupling with the atmosphere while stomatal control of energy partitioning was strongest in relatively sparse (LAI ~2 m² m^-2 ) pine stands growing on mineral soils. The data analysis and model results suggest that LAI may be an effective scaling parameter for net radiation and its partitioning but only in sparse stands (LAI <3 m² m^-2 ). This finding emphasizes the significance of stand-replacing disturbances on the controls of surface energy exchange. In denser forests, any LAI dependency varies with physiological traits such as light-saturated water-use efficiency. The results suggest that incorporating species traits and site conditions are necessary when LAI is used in upscaling energy exchanges of boreal coniferous forests.

Regional atmospheric cooling and wetting effect of permafrost thaw-induced boreal forest loss

In the sporadic permafrost zone of North America, thaw-induced boreal forest loss is leading to permafrost-free wetland expansion. These land cover changes alter landscape-scale surface properties with potentially large, however, still unknown impacts on regional climates. In this study, we combine nested eddy covariance flux tower measurements with satellite remote sensing to characterize the impacts of boreal forest loss on albedo, eco-physiological and aerodynamic surface properties, and turbulent energy fluxes of a lowland boreal forest region in the Northwest Territories, Canada. Planetary boundary layer modelling is used to estimate the potential forest loss impact on regional air temperature and atmospheric moisture. We show that thaw-induced conversion of forests to wetlands increases albedo: and bulk surface conductance for water vapour and decreases aerodynamic surface temperature. At the same time, heat transfer efficiency is reduced. These shifts in land surface properties increase latent at the expense of sensible heat fluxes, thus, drastically reducing Bowen ratios. Due to the lower albedo of forests and their masking effect of highly reflective snow, available energy is lower in wetlands, especially in late winter. Modelling results demonstrate that a conversion of a present-day boreal forest-wetland to a hypothetical homogeneous wetland landscape could induce a near-surface cooling effect on regional air temperatures of up to 3-4 °C in late winter and 1-2 °C in summer. An atmospheric wetting effect in summer is indicated by a maximum increase in water vapour mixing ratios of 2 mmol mol^-1 . At the same time, maximum boundary layer heights are reduced by about a third of the original height. In fall, simulated air temperature and atmospheric moisture between the two scenarios do not differ. Therefore, permafrost thaw-induced boreal forest loss may modify regional precipitation patterns and slow down regional warming trends.

Seasonal responses of terrestrial ecosystem water-use efficiency to climate change

Ecosystem water-use efficiency (EWUE) is an indicator of carbon-water interactions and is defined as the ratio of carbon assimilation (GPP) to evapotranspiration (ET). Previous research suggests an increasing long-term trend in annual EWUE over many regions and is largely attributed to the physiological effects of rising CO2 . The seasonal trends in EWUE, however, have not yet been analyzed. In this study, we investigate seasonal EWUE trends and responses to various drivers during 1982-2008. The seasonal cycle for two variants of EWUE, water-use efficiency (WUE, GPP/ET), and transpiration-based WUE (WUEt , the ratio of GPP and transpiration), is analyzed from 0.5° gridded fields from four process-based models and satellite-based products, as well as a network of 63 local flux tower observations. WUE derived from flux tower observations shows moderate seasonal variation for most latitude bands, which is in agreement with satellite-based products. In contrast, the seasonal EWUE trends are not well captured by the same satellite-based products. Trend analysis, based on process-model factorial simulations separating effects of climate, CO2 , and nitrogen deposition (NDEP), further suggests that the seasonal EWUE trends are mainly associated with seasonal trends of climate, whereas CO2 and NDEP do not show obvious seasonal difference in EWUE trends. About 66% grid cells show positive annual WUE trends, mainly over mid- and high northern latitudes. In these regions, spring climate change has amplified the effect of CO2 in increasing WUE by more than 0.005 gC m(-2)  mm(-1)  yr(-1) for 41% pixels. Multiple regression analysis further shows that the increase in springtime WUE in the northern hemisphere is the result of GPP increasing faster than ET because of the higher temperature sensitivity of GPP relative to ET. The partitioning of annual EWUE to seasonal components provides new insight into the relative sensitivities of GPP and ET to climate, CO2, and NDEP.

Nutritive quality and protein production from grain legumes in a boreal climate

BACKGROUND

Boreal cropping systems are heavily focused on the production of small-grain cereals; to improve their resilience to climate change and to achieve food and feed security, diversification is needed. This study investigated the potential of faba bean, narrow-leafed lupin and lentil as protein crops in southern Finland, where faba bean is traditional but the other two are novel.

RESULTS

Early cultivars of narrow-leafed lupin and lentil matured adequately. Protein concentration in faba bean was, at 32%, higher than the world average of 29%, while those of narrow-leafed lupin and lentil were close to their world averages. Protein yields decreased in the order faba bean > narrow-leafed lupin > lentil. Lipid content of faba bean and lentil was about 1.2% and that of narrow-leafed lupin about 5.5%, and fatty acid composition was largely oleic and linoleic in all three species.

CONCLUSION

Both lentil and narrow-leafed lupin can be added to the range of feed and food crops produced at high latitudes in Europe. While faba bean produces the greatest protein yield and lysine concentration, the higher sulfur amino acid concentration in lupin, its oil content and its adaptation to acid, sandy soils not suitable for faba bean make it an attractive alternative.

Metrics for biogeophysical climate forcings from land use and land cover changes and their inclusion in life cycle assessment: a critical review

The regulation by vegetation of heat, momentum, and moisture exchanges between the land surface and the atmosphere is a major component in Earth's climate system. By altering surface biogeophysics, anthropogenic land use activities often perturb these exchanges and thereby directly affect climate. Although long recognized scientifically as being important, biogeophysical climate forcings from land use and land cover changes (LULCC) are rarely included in life cycle assessment (LCA). Here, I review climate metrics for characterizing biogeophysical climate forcings from LULCC, focusing mostly on those that do not require coupled land-atmosphere climate models to compute. I discuss their merits, highlight their pros and cons in terms of their compatibility with the LCA framework, outline near-term practical guidelines and solutions for their integration, and point to areas of longer term research needs in both the climate science and LCA research communities.

Latent heat exchange in the boreal and arctic biomes

In this study latent heat flux (λE) measurements made at 65 boreal and arctic eddy-covariance (EC) sites were analyses by using the Penman-Monteith equation. Sites were stratified into nine different ecosystem types: harvested and burnt forest areas, pine forests, spruce or fir forests, Douglas-fir forests, broadleaf deciduous forests, larch forests, wetlands, tundra and natural grasslands. The Penman-Monteith equation was calibrated with variable surface resistances against half-hourly eddy-covariance data and clear differences between ecosystem types were observed. Based on the modeled behavior of surface and aerodynamic resistances, surface resistance tightly control λE in most mature forests, while it had less importance in ecosystems having shorter vegetation like young or recently harvested forests, grasslands, wetlands and tundra. The parameters of the Penman-Monteith equation were clearly different for winter and summer conditions, indicating that phenological effects on surface resistance are important. We also compared the simulated λE of different ecosystem types under meteorological conditions at one site. Values of λE varied between 15% and 38% of the net radiation in the simulations with mean ecosystem parameters. In general, the simulations suggest that λE is higher from forested ecosystems than from grasslands, wetlands or tundra-type ecosystems. Forests showed usually a tighter stomatal control of λE as indicated by a pronounced sensitivity of surface resistance to atmospheric vapor pressure deficit. Nevertheless, the surface resistance of forests was lower than for open vegetation types including wetlands. Tundra and wetlands had higher surface resistances, which were less sensitive to vapor pressure deficits. The results indicate that the variation in surface resistance within and between different vegetation types might play a significant role in energy exchange between terrestrial ecosystems and atmosphere. These results suggest the need to take into account vegetation type and phenology in energy exchange modeling.

Climate-driven uncertainties in modeling terrestrial energy and water fluxes: a site-level to global-scale analysis

We used a land surface model constrained using data from flux tower sites, to analyze the biases in ecosystem energy and water fluxes arising due to the use of meteorological reanalysis datasets. Following site-level model calibration encompassing major vegetation types from the tropics to the northern high-latitudes, we repeated the site and global simulations using two reanalysis datasets: the NCEP/NCAR and the CRUNCEP. In comparison with the model simulations using observed meteorology from sites, the reanalysis-driven simulations produced several systematic biases in net radiation (Rn ), latent heat (LE), and sensible heat (H) fluxes. These include: (i) persistently positive tropical/subtropical biases in Rn using the NCEP/NCAR, and gradually transitioning to negative Rn biases in the higher latitudes; (ii) large positive H biases in the tropics/subtropics using the NCEP/NCAR; (iii) negative LE biases using the NCEP/NCAR above 40°N; (iv) high tropical LE using the CRUNCEP in comparison with observationally derived global estimates; and (v) flux-partitioning biases from canopy and ground components. Across vegetation types, we investigated the role of the meteorological drivers (shortwave and longwave radiation, atmospheric humidity, temperature, precipitation) and their seasonal biases in controlling these reanalysis-driven uncertainties. At the global scale, our site-level analysis explains several model-data differences in the LE and H fluxes when compared with observationally derived global estimates of these fluxes. Using our results, we discuss the implications of site-level model calibration on subsequent regional/global applications to study energy and hydrological processes. The flux-partitioning biases presented in this study have potential implications on the couplings among terrestrial carbon, energy, and water fluxes, and for the calibration of land-atmosphere parameterizations that are dependent on LE/H partitioning.

Interannual variation of the Bowen ratio in a subtropical coniferous plantation in southeast China, 2003-2012

The interannual variation of the Bowen ratio, through its effect on the warming extent of available energy to the ecosystem land surface air, heavily influences the ecosystem microclimate and affects the hydrological cycle at both regional and global scales. Although the precipitation amount in southeast China is not expected to change greatly as a result of climate change, the precipitation frequency may be altered in the future. We explored the interannual variation of the Bowen ratio and its affecting mechanisms based on eddy covariance measurements in a subtropical plantation in southeast China during 2003-2012. The results indicated that the annual mean Bowen ratio was 0.35 ± 0.06, with a range of 0.29-0.45. The Bowen ratio during the dry season (July-October) positively correlated with the annual Bowen ratio (R(2) = 0.85, p<0.001). The effective precipitation frequency during the dry season, through its positive effect on shallow soil water content, indirectly and negatively affected the annual Bowen ratio. Between 2003 and 2012, the annual Bowen ratio exhibited a marginally significant decreasing trend (p = 0.061), meanwhile the effective precipitation frequency and shallow soil water content during the dry season increased significantly (p<0.001). The annual Bowen ratio may decrease further if the effective precipitation frequency and shallow soil water content during the dry season follow similar trends in the future. The warming effect of available energy to the surface air of our studied plantation may decline with the decreasing annual Bowen ratio.

Climate change implications of shifting forest management strategy in a boreal forest ecosystem of Norway

Empirical models alongside remotely sensed and station measured meteorological observations are employed to investigate both the local and global direct climate change impacts of alternative forest management strategies within a boreal ecosystem of eastern Norway. Stand-level analysis is firstly executed to attribute differences in daily, seasonal, and annual mean surface temperatures to differences in surface intrinsic biophysical properties across conifer, deciduous, and clear-cut sites. Relative to a conifer site, a slight local cooling of −0.13 °C at a deciduous site and −0.25 °C at a clear-cut site were observed over a 6-year period, which were mostly attributed to a higher albedo throughout the year. When monthly mean albedo trajectories over the entire managed forest landscape were taken into consideration, we found that strategies promoting natural regeneration of coniferous sites with native deciduous species led to substantial global direct climate cooling benefits relative to those maintaining current silviculture regimes – despite predicted long-term regional warming feedbacks and a reduced albedo in spring and autumn months. The magnitude and duration of the cooling benefit depended largely on whether management strategies jointly promoted an enhanced material supply over business-as-usual levels. Expressed in terms of an equivalent CO2 emission pulse at the start of the simulation, the net climate response at the end of the 21st century spanned −8 to −159 Tg-CO2-eq., depending on whether near-term harvest levels increased or followed current trends, respectively. This magnitude equates to approximately −20 to −300% of Norway's annual domestic (production) emission impact. Our analysis supports the assertion that a carbon-only focus in the design and implementation of forest management policy in boreal and other climatically similar regions can be counterproductive – and at best – suboptimal if boreal forests are to be used as a tool to mitigate global warming.

Simulating effects of fire disturbance and climate change on boreal forest productivity and evapotranspiration

We used a terrestrial ecosystem process model, BIOME-BGC, to investigate historical climate change and fire disturbance effects on regional carbon and water budgets within a 357,500 km(2) portion of the Canadian boreal forest. Historical patterns of increasing atmospheric CO2, climate change, and regional fire activity were used as model drivers to evaluate the relative effects of these impacts to spatial patterns and temporal trends in forest net primary production (NPP) and evapotranspiration (ET). Historical trends of increasing atmospheric CO2 resulted in overall 13% and 5% increases in annual NPP and ET from 1994 to 1996, respectively. NPP was found to be relatively sensitive to changes in air temperature (T(a)), while ET was more sensitive to precipitation (P) change within the ranges of observed climate variability (e.g., +/-2 degrees C for T(a) and +/-20% for P). In addition, the potential effect of climate change related warming on NPP is exacerbated or offset depending on whether these changes are accompanied by respective decreases or increases in precipitation. Historical fire activity generally resulted in reductions of both NPP and ET, which consumed an average of approximately 6% of annual NPP from 1959 to 1996. Areas currently occupied by dry conifer forests were found to be subject to more frequent fire activity, which consumed approximately 8% of annual NPP. The results of this study show that the North American boreal ecosystem is sensitive to historical patterns of increasing atmospheric CO2, climate change and regional fire activity. The relative impacts of these disturbances on NPP and ET interact in complex ways and are spatially variable depending on regional land cover and climate gradients.