Simulating the effects of climate change on the carbon balance of North American high-latitude forests.

The large magnitude of predicted warming at high latitudes and the potential feedback of ecosystems to atmospheric CO2 concentrations make it important to quantify both warming and its effects on high-latitude carbon balance. We analysed long-term, daily surface meteorological records for 13 sites in Alaska and north-western Canada and an 82-y record of river ice breakup date for the Tanana River in interior Alaska. We found increases in winter and spring temperature extrema for all sites, with the greatest increases in spring minimum temperature, average 0.47 °C per 10 y, and a 0.7-day per 10 y advance in ice breakup on the Tanana River. We used the climate records to drive an ecosystem process model, BIOME_BGC, to simulate the effects of climate change on the carbon and water balances of boreal forest ecosystems. The growing season has lengthened by an average of 2.6 days per 10 y with an advance in average leaf onset date of 1.10 days per 10 y. This advance in the start of the active growing season correlates positively with progressively earlier ice breakup on the Tanana River in interior Alaska. The advance in the start of the growing season resulted in a 20% increase in net primary production for both aspen (Populus tremuloides) and white spruce (Picea glauca) stands. Aspen had a greater mean increase in maintenance respiration than spruce, whereas spruce had a greater mean increase in evapotranspiration. Average decomposition rates also increased for both species. Both net primary production and decomposition are enhanced in our simulations, suggesting that productive forest types may not experience a significant shift in net carbon flux as a result of climate warming.

Arctic and boreal ecosystems of western North America as components of the climate system.

Synthesis of results from several Arctic and boreal research programmes provides evidence for the strong role of high-latitude ecosystems in the climate system. Average surface air temperature has increased 0.3 °C per decade during the twentieth century in the western North American Arctic and boreal forest zones. Precipitation has also increased, but changes in soil moisture are uncertain. Disturbance rates have increased in the boreal forest; for example, there has been a doubling of the area burned in North America in the past 20 years. The disturbance regime in tundra may not have changed. Tundra has a 3-6-fold higher winter albedo than boreal forest, but summer albedo and energy partitioning differ more strongly among ecosystems within either tundra or boreal forest than between these two biomes. This indicates a need to improve our understanding of vegetation dynamics within, as well as between, biomes. If regional surface warming were to continue, changes in albedo and energy absorption would likely act as a positive feedback to regional warming due to earlier melting of snow and, over the long term, the northward movement of treeline. Surface drying and a change in dominance from mosses to vascular plants would also enhance sensible heat flux and regional warming in tundra. In the boreal forest of western North America, deciduous forests have twice the albedo of conifer forests in both winter and summer, 50-80% higher evapotranspiration, and therefore only 30-50% of the sensible heat flux of conifers in summer. Therefore, a warming-induced increase in fire frequency that increased the proportion of deciduous forests in the landscape, would act as a negative feedback to regional warming. Changes in thermokarst and the aerial extent of wetlands, lakes, and ponds would alter high-latitude methane flux. There is currently a wide discrepancy among estimates of the size and direction of CO2 flux between high-latitude ecosystems and the atmosphere. These discrepancies relate more strongly to the approach and assumptions for extrapolation than to inconsistencies in the underlying data. Inverse modelling from atmospheric CO2 concentrations suggests that high latitudes are neutral or net sinks for atmospheric CO2 , whereas field measurements suggest that high latitudes are neutral or a net CO2 source. Both approaches rely on assumptions that are difficult to verify. The most parsimonious explanation of the available data is that drying in tundra and disturbance in boreal forest enhance CO2 efflux. Nevertheless, many areas of both tundra and boreal forests remain net sinks due to regional variation in climate and local variation in topographically determined soil moisture. Improved understanding of the role of high-latitude ecosystems in the climate system requires a concerted research effort that focuses on geographical variation in the processes controlling land-atmosphere exchange, species composition, and ecosystem structure. Future studies must be conducted over a long enough time-period to detect and quantify ecosystem feedbacks.

Physiological control of water flux in conifers : A computer simulation model.

A water flux model with daily resolution is described which permits one to assess how changes in the rooting volume, amount of sapwood, leaf area and conductance properties interact to affect water uptake, internal storage, and transpiration.A root zone water compartment is defined for a particular tree on the basis of root depth, lateral extension and moisture holding characteristics of the soil. Water is taken up from different subcompartments of the root zone as a function of vertical position, soil water content, and water deficit within the sapwood. Excess water entering the root zone is channeled into runoff or seepage.The sapwood compartment of the model is restricted to the main stem of the tree and does not include sapwood in the branches or roots. The model assumes whatever water deficit is built up in the sapwood during the day will be replenished at night if the root zone water supply/capacity ratio exceeds 20%. A complex exponential equation describes the amount of water extractable from 20% to 0 capacity when no uptake is possible. The maximum change in volume of water in the sapwood of a large Douglas-fir is estimated to represent more than a 10 day supply for transpiration.Water loss through transpiration is predicted as a function of the mean daily absolute humidity deficit, leaf area, leaf conductance and daylength. Leaf conductance is controlled by predawn plant moisture stress which in turn is a function of the rooting zone water supply.The model incorporates two special constraints upon water uptake and transpiration. The first accounts for the effect of cold soil temperatures reducing the possible uptake by Douglas-fir to half at 2°C and to 0 at-2°C. The second represents a critical absolute humidity deficit sufficient to cause stomatal closure which results in leaf conductance being reduced to a minimum.The model is employed to compare trees of different sizes and those with different stomatal behavior. From this experience, it is suggested that future studies include, at a minimum, simultaneous measurements of: absolute humidity deficit, leaf area, sapwood volume and change in water content, predawn stress and leaf conductance.

Soil Temperature Influences on Root Resistance of Pinus contorta Seedlings.

The influence of low temperature in the root zone on water uptake in lodgepole pine (Pinus contorta Dougl. ex Loud.) was studied under laboratory conditions. To remove soil hydraulic influences, two-year-old seedlings were transferred to solution cultures and maintained in temperature controlled water baths. Short term measurements of leaf conductance, leaf water potential and tritiated water movement were taken at root temperatures from 22 C down to 0 C. Root resistance was calculated to be 67% of total plant resistance at 7 C and 93% at 0 C. In addition an Arrhenius break was found in a plant resistance versus temperature plot, suggesting a significant change with temperature in the membrane pathway in the root water uptake system.

Simulating the effects of climatic variation on stem carbon accumulation of a ponderosa pine stand: comparison with annual growth increment data.

Simulation models of ecosystem processes may be necessary to separate the long-term effects of climate change on forest productivity from the effects of year-to-year variations in climate. The objective of this study was to compare simulated annual stem growth with measured annual stem growth from 1930 to 1982 for a uniform stand of ponderosa pine (Pinus ponderosa Dougl.) in Montana, USA. The model, FOREST-BGC, was used to simulate growth assuming leaf area index (LAI) was either constant or increasing. The measured stem annual growth increased exponentially over time; the differences between the simulated and measured stem carbon accumulations were not large. Growth trends were removed from both the measured and simulated annual increments of stem carbon to enhance the year-to-year variations in growth resulting from climate. The detrended increments from the increasing LAI simulation fit the detrended increments of the stand data over time with an R(2) of 0.47; the R(2) increased to 0.65 when the previous year's simulated detrended increment was included with the current year's simulated increment to account for autocorrelation. Stepwise multiple linear regression of the detrended increments of the stand data versus monthly meteorological variables had an R(2) of 0.37, and the R(2) increased to 0.47 when the previous year's meteorological data were included to account for autocorrelation. Thus, FOREST-BGC was more sensitive to the effects of year-to-year climate variation on annual stem growth than were multiple linear regression models.

Simulating effects of fire on northern Rocky Mountain landscapes with the ecological process model FIRE-BGC.

A mechanistic, biogeochemical succession model, FIRE-BGC, was used to investigate the role of fire on long-term landscape dynamics in northern Rocky Mountain coniferous forests of Glacier National Park, Montana, USA. FIRE-BGC is an individual-tree model-created by merging the gap-phase process-based model FIRESUM with the mechanistic ecosystem biogeochemical model FOREST-BGC-that has mixed spatial and temporal resolution in its simulation architecture. Ecological processes that act at a landscape level, such as fire and seed dispersal, are simulated annually from stand and topographic information. Stand-level processes, such as tree establishment, growth and mortality, organic matter accumulation and decomposition, and undergrowth plant dynamics are simulated both daily and annually. Tree growth is mechanistically modeled based on the ecosystem process approach of FOREST-BGC where carbon is fixed daily by forest canopy photosynthesis at the stand level. Carbon allocated to the tree stem at the end of the year generates the corresponding diameter and height growth. The model also explicitly simulates fire behavior and effects on landscape characteristics. We simulated the effects of fire on ecosystem characteristics of net primary productivity, evapotranspiration, standing crop biomass, nitrogen cycling and leaf area index over 200 years for the 50,000-ha McDonald Drainage in Glacier National Park. Results show increases in net primary productivity and available nitrogen when fires are included in the simulation. Standing crop biomass and evapotranspiration decrease under a fire regime. Shade-intolerant species dominate the landscape when fires are excluded. Model tree increment predictions compared well with field data.

Regional assessment of boreal forest productivity using an ecological process model and remote sensing parameter maps.

An ecological process model (BIOME-BGC) was used to assess boreal forest regional net primary production (NPP) and response to short-term, year-to-year weather fluctuations based on spatially explicit, land cover and biomass maps derived by radar remote sensing, as well as soil, terrain and daily weather information. Simulations were conducted at a 30-m spatial resolution, over a 1205 km(2) portion of the BOREAS Southern Study Area of central Saskatchewan, Canada, over a 3-year period (1994-1996). Simulations of NPP for the study region were spatially and temporally complex, averaging 2.2 (+/- 0.6), 1.8 (+/- 0.5) and 1.7 (+/- 0.5) Mg C ha(-1) year(-1) for 1994, 1995 and 1996, respectively. Spatial variability of NPP was strongly controlled by the amount of aboveground biomass, particularly photosynthetic leaf area, whereas biophysical differences between broadleaf deciduous and evergreen coniferous vegetation were of secondary importance. Simulations of NPP were strongly sensitive to year-to-year variations in seasonal weather patterns, which influenced the timing of spring thaw and deciduous bud-burst. Reductions in annual NPP of approximately 17 and 22% for 1995 and 1996, respectively, were attributed to 3- and 5-week delays in spring thaw relative to 1994. Boreal forest stands with greater proportions of deciduous vegetation were more sensitive to the timing of spring thaw than evergreen coniferous stands. Similar relationships were found by comparing simulated snow depth records with 10-year records of aboveground NPP measurements obtained from biomass harvest plots within the BOREAS region. These results highlight the importance of sub-grid scale land cover complexity in controlling boreal forest regional productivity, the dynamic response of the biome to short-term interannual climate variations, and the potential implications of climate change and other large-scale disturbances.