In terrestrial high‐latitude regions, observations indicate recent changes in snow cover, permafrost, and soil freeze–thaw transitions due to climate change. These modifications may result in ...temporal shifts in the growing season and the associated rates of terrestrial productivity. Changes in productivity will influence the ability of these ecosystems to sequester atmospheric CO2. We use the terrestrial ecosystem model (TEM), which simulates the soil thermal regime, in addition to terrestrial carbon (C), nitrogen and water dynamics, to explore these issues over the years 1960–2100 in extratropical regions (30–90°N). Our model simulations show decreases in snow cover and permafrost stability from 1960 to 2100. Decreases in snow cover agree well with National Oceanic and Atmospheric Administration satellite observations collected between the years 1972 and 2000, with Pearson rank correlation coefficients between 0.58 and 0.65. Model analyses also indicate a trend towards an earlier thaw date of frozen soils and the onset of the growing season in the spring by approximately 2–4 days from 1988 to 2000. Between 1988 and 2000, satellite records yield a slightly stronger trend in thaw and the onset of the growing season, averaging between 5 and 8 days earlier. In both, the TEM simulations and satellite records, trends in day of freeze in the autumn are weaker, such that overall increases in growing season length are due primarily to earlier thaw. Although regions with the longest snow cover duration displayed the greatest increase in growing season length, these regions maintained smaller increases in productivity and heterotrophic respiration than those regions with shorter duration of snow cover and less of an increase in growing season length. Concurrent with increases in growing season length, we found a reduction in soil C and increases in vegetation C, with greatest losses of soil C occurring in those areas with more vegetation, but simulations also suggest that this trend could reverse in the future. Our results reveal noteworthy changes in snow, permafrost, growing season length, productivity, and net C uptake, indicating that prediction of terrestrial C dynamics from one decade to the next will require that large‐scale models adequately take into account the corresponding changes in soil thermal regimes.
Wildfire is a common occurrence in ecosystems of northern high latitudes, and changes in the fire regime of this region have consequences for carbon feedbacks to the climate system. To improve our ...understanding of how wildfire influences carbon dynamics of this region, we used the process‐based Terrestrial Ecosystem Model to simulate fire emissions and changes in carbon storage north of 45°N from the start of spatially explicit historically recorded fire records in the twentieth century through 2002, and evaluated the role of fire in the carbon dynamics of the region within the context of ecosystem responses to changes in atmospheric CO2 concentration and climate. Our analysis indicates that fire plays an important role in interannual and decadal scale variation of source/sink relationships of northern terrestrial ecosystems and also suggests that atmospheric CO2 may be important to consider in addition to changes in climate and fire disturbance. There are substantial uncertainties in the effects of fire on carbon storage in our simulations. These uncertainties are associated with sparse fire data for northern Eurasia, uncertainty in estimating carbon consumption, and difficulty in verifying assumptions about the representation of fires that occurred prior to the start of the historical fire record. To improve the ability to better predict how fire will influence carbon storage of this region in the future, new analyses of the retrospective role of fire in the carbon dynamics of northern high latitudes should address these uncertainties.
We carried out a field experiment to evaluate the effect of moisture regime on microbial biomass and activity in birch litter in the Alaskan taiga. Litter bags were placed in one of three treatments: ...continuously moist (0.5 cm water d
−1), cycling (0.5 cm water weekly), and `natural', which experienced two natural dry–wet cycles of 2 weeks dry followed by rain. The experiment lasted for 1 month. Each week we collected litter bags and analyzed microbial respiration and biomass C and N. In the last two cycles we analyzed bacterial substrate use on Biolog GN plates. There were strong overall correlations between biomass, respiration and litter moisture content. However, the different treatments had significantly different rates of respiration, biomass and respiratory quotient (
qCO
2) that could not be explained by moisture content directly. The natural treatment had lower respiration rates and biomass than the wet or cycling samples, indicating that the 2-week droughts in the natural treatment reduced microbial populations and activity to a greater degree than did shorter droughts. Episodic drying and rewetting considerably decreased the number of Biolog substrates used. This experiment showed that the size and functioning of the litter microbial community was strongly affected by its stress history.
Four tundra and taiga soils were experimentally subjected to three freeze-thaw cycles (5 days each at −5°C and +5°C). After each thaw, there was an initial pulse (<24 h) in microbial respiration. The ...total amount of C respired in each thaw period was largest during the first cycle and decreased in successive cycles. Three cycles caused a net increase in total respired C relative to the +5°C control in wet meadow tundra soil, a reduction in birch soil, and had no net effect in either alder-poplar or tussock tundra soil. These different patterns apparently resulted from differences in the quality of the soil organic matter and the relative activity of the microbial biomass. Net N mineralization was generally enhanced relative to the +5°C controls in the first cycle, but was inhibited in the third cycle, similar to what was observed with C. Over multiple freeze-thaw cycles, the initial response of C and N mineralization appear to be driven by release from the microbial biomass, while over the longer-term, the response is driven by the reduction in attack on soil organic matter resulting from a reduced microbial population.
ABSTRACT
There is substantial evidence that soil thermal dynamics are changing in terrestrial ecosystems of the Northern Hemisphere and that these dynamics have implications for the exchange of ...carbon between terrestrial ecosystems and the atmosphere. To date, large‐scale biogeochemical models have been slow to incorporate the effects of soil thermal dynamics on processes that affect carbon exchange with the atmosphere. In this study we incorporated a soil thermal module (STM), appropriate to both permafrost and non‐permafrost soils, into a large‐scale ecosystem model, version 5.0 of the Terrestrial Ecosystem Model (TEM). We then compared observed regional and seasonal patterns of atmospheric CO2 to simulations of carbon dynamics for terrestrial ecosystems north of 30°N between TEM 5.0 and an earlier version of TEM (version 4.2) that lacked a STM. The timing of the draw‐down of atmospheric CO2 at the start of the growing season and the degree of draw‐down during the growing season were substantially improved by the consideration of soil thermal dynamics. Both versions of TEM indicate that climate variability and change promoted the loss of carbon from temperate ecosystems during the first half of the 20th century, and promoted carbon storage during the second half of the century. The results of the simulations by TEM suggest that land‐use change in temperate latitudes (30–60°N) plays a stronger role than climate change in driving trends for increased uptake of carbon in extratropical terrestrial ecosystems (30–90°N) during recent decades. In the 1980s the TEM 5.0 simulation estimated that extratropical terrestrial ecosystems stored 0.55 Pg C yr−1, with 0.24 Pg C yr−1 in North America and 0.31 Pg C yr−1 in northern Eurasia. From 1990 through 1995 the model simulated that these ecosystems stored 0.90 Pg C yr−1, with 0.27 Pg C yr−1 stored in North America and 0.63 Pg C yr−1 stored in northern Eurasia. Thus, in comparison to the 1980s, simulated net carbon storage in the 1990s was enhanced by an additional 0.35 Pg C yr−1 in extratropical terrestrial ecosystems, with most of the additional storage in northern Eurasia. The carbon storage simulated by TEM 5.0 in the 1980s and 1990s was lower than estimates based on other methodologies, including estimates by atmospheric inversion models and remote sensing and inventory analyses. This suggests that other issues besides the role of soil thermal dynamics may be responsible, in part, for the temporal and spatial dynamics of carbon storage of extratropical terrestrial ecosystems. In conclusion, the consideration of soil thermal dynamics and terrestrial cryospheric processes in modeling the global carbon cycle has helped to reduce biases in the simulation of the seasonality of carbon dynamics of extratropical terrestrial ecosystems. This progress should lead to an enhanced ability to clarify the role of other issues that influence carbon dynamics in terrestrial regions that experience seasonal freezing and thawing of soil.
Nine ecosystem process models were used to predict CO2 and water vapor exchanges by a 150‐year‐old black spruce forest in central Canada during 1994–1996 to evaluate and improve the models. Three ...models had hourly time steps, five had daily time steps, and one had monthly time steps. Model input included site ecosystem characteristics and meteorology. Model predictions were compared to eddy covariance (EC) measurements of whole‐ecosystem CO2 exchange and evapotranspiration, to chamber measurements of nighttime moss‐surface CO2 release, and to ground‐based estimates of annual gross primary production, net primary production, net ecosystem production (NEP), plant respiration, and decomposition. Model‐model differences were apparent for all variables. Model‐measurement agreement was good in some cases but poor in others. Modeled annual NEP ranged from −11 g C m−2 (weak CO2 source) to 85 g C m−2 (moderate CO2 sink). The models generally predicted greater annual CO2 sink activity than measured by EC, a discrepancy consistent with the fact that model parameterizations represented the more productive fraction of the EC tower “footprint.” At hourly to monthly timescales, predictions bracketed EC measurements so median predictions were similar to measurements, but there were quantitatively important model‐measurement discrepancies found for all models at subannual timescales. For these models and input data, hourly time steps (and greater complexity) compared to daily time steps tended to improve model‐measurement agreement for daily scale CO2 exchange and evapotranspiration (as judged by root‐mean‐squared error). Model time step and complexity played only small roles in monthly to annual predictions.
Summary
Historical and projected climate trends for high latitudes show substantial temporal and spatial variability. To identify uncertainties in simulating carbon (C) dynamics for pan‐Arctic ...tundra, we compare the historical and projected responses of tundra C storage from 1921 to 2100 between simulations by the Terrestrial Ecosystem Model (TEM) for the pan‐Arctic and the Kuparuk River Basin, which was the focus of an integrated study of C dynamics from 1994 to 1996. In the historical period from 1921 to 1994, the responses of net primary production (NPP) and heterotrophic respiration (RH) simulated for the Kuparuk River Basin and the pan‐Arctic are correlated with the same factors; NPP is positively correlated with net nitrogen mineralization (NMIN) and RH is negatively correlated with mean annual soil moisture. In comparison to the historical period, the spatially aggregated responses of NPP and RH for the Kuparuk River Basin and the pan‐Arctic in our simulations for the projected period have different sensitivities to temperature, soil moisture and NMIN. In addition to being sensitive to soil moisture during the projected period, RH is also sensitive to temperature and there is a significant correlation between RH and NMIN. We interpret the increases in NPP during the projected period as being driven primarily by increases in NMIN, and that the correlation between NPP and temperature in the projected period is a result primarily of the causal linkage between temperature, RH, and NMIN. Although similar factors appear to be controlling simulated regional‐and biome‐scale C dynamics, simulated C dynamics at the two scales differ in magnitude with higher increases in C storage simulated for the Kuparuk River Basin than for the pan‐Arctic at the end of the historical period and throughout the projected period. Also, the results of the simulations indicate that responses of C storage show different climate sensitivities at regional and pan‐Arctic spatial scales and that these sensitivities change across the temporal scope of the simulations. The results of the TEM simulations indicate that the scaling of C dynamics to a region of arctic tundra may not represent C dynamics of pan‐Arctic tundra because of the limited spatial variation in climate and vegetation within a region relative to the pan‐Arctic. For reducing uncertainties, our analyses highlight the importance of incorporating the understanding gained from process‐level studies of C dynamics in a region of arctic tundra into process‐based models that simulate C dynamics in a spatially explicit fashion across the spatial domain of pan‐Arctic tundra. Also, efforts to improve gridded datasets of
historical climate for the pan‐Arctic would advance the ability to assess the responses of C dynamics for pan‐Arctic tundra in a more realistic fashion. A major challenge will be to incorporate topographic controls over soil moisture in assessing the response of C storage for pan‐Arctic tundra.
Summary
We are developing a process‐based modelling approach to investigate how carbon (C) storage of tundra across the entire Arctic will respond to projected climate change. To implement the ...approach, the processes that are least understood, and thus have the most uncertainty, need to be identified and studied. In this paper, we identified a key uncertainty by comparing the responses of C storage in tussock tundra at one site between the simulations of two models – one a global‐scale ecosystem model (Terrestrial Ecosystem Model, TEM) and one a plot‐scale ecosystem model (General Ecosystem Model, GEM). The simulations spanned the historical period (1921–94) and the projected period (1995–2100). In the historical period, the model simulations of net primary production (NPP) differed in their sensitivity to variability in climate. However, the long‐term changes in C storage were similar in both simulations, because the dynamics of heterotrophic respiration (RH) were similar in both models. In contrast, the responses of C storage in the two model simulations diverged during the projected period. In the GEM simulation for this period, increases in RH tracked increases in NPP, whereas in the TEM simulation increases in RH lagged increases in NPP. We were able to make the long‐term C dynamics of the two simulations agree by parameterizing TEM to the fast soil C pools of GEM. We concluded that the differences between the long‐term C dynamics of the two simulations lay in modelling the role of the recalcitrant soil C. These differences, which reflect an incomplete understanding of soil processes, lead to quite different projections of the response of pan‐Arctic C storage to global change. For example, the reference parameterization of TEM resulted in an estimate of cumulative C storage of 2032 g C m−2 for moist tundra north of 50°N, which was substantially higher than the 463 g C m−2 estimated for a parameterization of fast soil C dynamics. This uncertainty in the depiction of the role of recalcitrant soil C in long‐term ecosystem C dynamics resulted from our incomplete understanding of controls over C and N transformations in Arctic soils. Mechanistic studies of these issues are needed to improve our ability to model the response of Arctic ecosystems to global change.
The role of carbon (C) and nitrogen (N) interactions on sequestration of atmospheric CO2 in black spruce ecosystems across North America was evaluated with the Terrestrial Ecosystem Model (TEM) by ...applying parameterizations of the model in which C–N dynamics were either coupled or uncoupled. First, the performance of the parameterizations, which were developed for the dynamics of black spruce ecosystems at the Bonanza Creek Long-Term Ecological Research site in Alaska, were evaluated by simulating C dynamics at eddy correlation tower sites in the Boreal Ecosystem Atmosphere Study (BOREAS) for black spruce ecosystems in the northern study area (northern site) and the southern study area (southern site) with local climate data. We compared simulated monthly growing season (May to September) estimates of gross primary production (GPP), total ecosystem respiration (RESP), and net ecosystem production (NEP) from 1994 to 1997 to available field-based estimates at both sites. At the northern site, monthly growing season estimates of GPP and RESP for the coupled and uncoupled simulations were highly correlated with the field-based estimates (coupled: R2= 0.77, 0.88 for GPP and RESP; uncoupled: R2 = 0.67, 0.92 for GPP and RESP). Although the simulated seasonal pattern of NEP generally matched the field-based data, the correlations between field-based and simulated monthly growing season NEP were lower (R2 = 0.40, 0.00 for coupled and uncoupled simulations, respectively) in comparison to the correlations between field-based and simulated GPP and RESP. The annual NEP simulated by the coupled parameterization fell within the uncertainty of field-based estimates in two of three years. On the other hand, annual NEP simulated by the uncoupled parameterization only fell within the field-based uncertainty in one of three years. At the southern site, simulated NEP generally matched field-based NEP estimates, and the correlation between monthly growing season field-based and simulated NEP (R2 = 0.36, 0.20 for coupled and uncoupled simulations, respectively) was similar to the correlations at the northern site. To evaluate the role of N dynamics in C balance of black spruce ecosystems across North America, we simulated historical and projected C dynamics from 1900 to 2100 with a global-based climatology at 0.5° resolution (latitude × longitude) with both the coupled and uncoupled parameterizations of TEM. From analyses at the northern site, several consistent patterns emerge. There was greater inter-annual variability in net primary production (NPP) simulated by the uncoupled parameterization as compared to the coupled parameterization, which led to substantial differences in inter-annual variability in NEP between the parameterizations. The divergence between NPP and heterotrophic respiration was greater in the uncoupled simulation, resulting in more C sequestration during the projected period. These responses were the result of fundamentally different responses of the coupled and uncoupled parameterizations to changes in CO2 and climate. Across North American black spruce ecosystems, the range of simulated decadal changes in C storage was substantially greater for the uncoupled parameterization than for the coupled parameterization. Analysis of the spatial variability in decadal responses of C dynamics revealed that C fluxes simulated by the coupled and uncoupled parameterizations have different sensitivities to climate and that the climate sensitivities of the fluxes change over the temporal scope of the simulations. The results of this study suggest that uncertainties can be reduced through (1) factorial studies focused on elucidating the role of C and N interactions in the response of mature black spruce ecosystems to manipulations of atmospheric CO2 and climate, (2) establishment of a network of continuous, long-term measurements of C dynamics across the range of mature black spruce ecosystems in North America, and (3) ancillary measurements in the network to elucidate the role of C and N interactions in exchange of CO2 with the atmosphere.