skip to main content

Title: Assessing dynamic vegetation model parameter uncertainty across Alaskan arctic tundra plant communities
As the Arctic region moves into uncharted territory under a warming climate, it is important to refine the terrestrial biosphere models (TBMs) that help us understand and predict change. One fundamental uncertainty in TBMs relates to model parameters, configuration variables internal to the model whose value can be estimated from data. We incorporate a version of the Terrestrial Ecosystem Model (TEM) developed for arctic ecosystems into the Predictive Ecosystem Analyzer (PEcAn) framework. PEcAn treats model parameters as probability distributions, estimates parameters based on a synthesis of available field data, and then quantifies both model sensitivity and uncertainty to a given parameter or suite of parameters. We examined how variation in 21 parameters in the equation for gross primary production influenced model sensitivity and uncertainty in terms of two carbon fluxes (net primary productivity and heterotrophic respiration) and two carbon (C) pools (vegetation C and soil C). We set up different parameterizations of TEM across a range of tundra types (tussock tundra, heath tundra, wet sedge tundra, and shrub tundra) in northern Alaska, along a latitudinal transect extending from the coastal plain near Utqiaġvik to the southern foothills of the Brooks Range, to the Seward Peninsula. TEM was most sensitive to more » parameters related to the temperature regulation of photosynthesis. Model uncertainty was mostly due to parameters related to leaf area, temperature regulation of photosynthesis, and the stomatal responses to ambient light conditions. Our analysis also showed that sensitivity and uncertainty to a given parameter varied spatially. At some sites, model sensitivity and uncertainty tended to be connected to a wider range of parameters, underlining the importance of assessing tundra community processes across environmental gradients or geographic locations. Generally, across sites, the flux of net primary productivity (NPP) and pool of vegetation C had about equal uncertainty, while heterotrophic respiration had higher uncertainty than the pool of soil C. Our study illustrates the complexity inherent in evaluating parameter uncertainty across highly heterogeneous arctic tundra plant communities. It also provides a framework for iteratively testing how newly collected field data related to key parameters may result in more effective forecasting of Arctic change. « less
Authors:
; ; ; ; ; ; ;
Award ID(s):
1636476 1936752
Publication Date:
NSF-PAR ID:
10313984
Journal Name:
Ecological Applications
ISSN:
1051-0761
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Global estimates of the land carbon sink are often based on simulations by terrestrial biosphere models (TBMs). The use of a large number of models that differ in their underlying hypotheses, structure and parameters is one way to assess the uncertainty in the historical land carbon sink. Here we show that the atmospheric forcing datasets used to drive these TBMs represent a significant source of uncertainty that is currently not systematically accounted for in land carbon cycle evaluations. We present results from three TBMs each forced with three different historical atmospheric forcing reconstructions over the period 1850–2015. We perform an analysis of variance to quantify the relative uncertainty in carbon fluxes arising from the models themselves, atmospheric forcing, and model-forcing interactions. We find that atmospheric forcing in this set of simulations plays a dominant role on uncertainties in global gross primary productivity (GPP) (75% of variability) and autotrophic respiration (90%), and a significant but reduced role on net primary productivity and heterotrophic respiration (30%). Atmospheric forcing is the dominant driver (52%) of variability for the net ecosystem exchange flux, defined as the difference between GPP and respiration (both autotrophic and heterotrophic respiration). In contrast, for wildfire-driven carbon emissions modelmore »uncertainties dominate and, as a result, model uncertainties dominate for net ecosystem productivity. At regional scales, the contribution of atmospheric forcing to uncertainty shows a very heterogeneous pattern and is smaller on average than at the global scale. We find that this difference in the relative importance of forcing uncertainty between global and regional scales is related to large differences in regional model flux estimates, which partially offset each other when integrated globally, while the flux differences driven by forcing are mainly consistent across the world and therefore add up to a larger fractional contribution to global uncertainty.

    « less
  2. Abstract

    An estimated 1700 Pg of carbon is frozen in the Arctic permafrost and the fate of this carbon is unclear because of the complex interaction of biophysical, ecological and biogeochemical processes that govern the Arctic carbon budget. Two key processes determining the region’s long-term carbon budget are: (a) carbon uptake through increased plant growth, and (b) carbon release through increased heterotrophic respiration (HR) due to warmer soils. Previous predictions for how these two opposing carbon fluxes may change in the future have varied greatly, indicating that improved understanding of these processes and their feedbacks is critical for advancing our predictive ability for the fate of Arctic peatlands. In this study, we implement and analyze a vertically-resolved model of peatland soil carbon into a cohort-based terrestrial biosphere model to improve our understanding of how on-going changes in climate are altering the Arctic carbon budget. A key feature of the formulation is that accumulation of peat within the soil column modifies its texture, hydraulic conductivity, and thermal conductivity, which, in turn influences resulting rates of HR within the soil column. Analysis of the model at three eddy covariance tower sites in the Alaskan tundra shows that the vertically-resolved soil column formulationmore »accurately captures the zero-curtain phenomenon, in which the temperature of soil layers remain at or near 0 °C during fall freezeback due to the release of latent heat, is critical to capturing observed patterns of wintertime respiration. We find that significant declines in net ecosystem productivity (NEP) occur starting in 2013 and that these declines are driven by increased HR arising from increased precipitation and warming. Sensitivity analyses indicate that the cumulative NEP over the decade responds strongly to the estimated soil carbon stock and more weakly to vegetation abundance at the beginning of the simulation.

    « less
  3. Abstract

    Arctic ecosystems are particularly vulnerable to climate change because of Arctic amplification. Here, we assessed the climatic impacts of low-end, 1.5 °C, and 2.0 °C global temperature increases above pre-industrial levels, on the warming of terrestrial ecosystems in northern high latitudes (NHL, above 60 °N including pan-Arctic tundra and boreal forests) under the framework of the Inter-Sectoral Impact Model Intercomparison Project phase 2b protocol. We analyzed the simulated changes of net primary productivity, vegetation biomass, and soil carbon stocks of eight ecosystem models that were forced by the projections of four global climate models and two atmospheric greenhouse gas pathways (RCP2.6 and RCP6.0). Our results showed that considerable impacts on ecosystem carbon budgets, particularly primary productivity and vegetation biomass, are very likely to occur in the NHL areas. The models agreed on increases in primary productivity and biomass accumulation, despite considerable inter-model and inter-scenario differences in the magnitudes of the responses. The inter-model variability highlighted the inadequacies of the present models, which fail to consider important components such as permafrost and wildfire. The simulated impacts were attributable primarily to the rapid temperature increases in the NHL and the greater sensitivity of northern vegetation to warming, which contrasted with themore »less pronounced responses of soil carbon stocks. The simulated increases of vegetation biomass by 30–60 Pg C in this century have implications for climate policy such as the Paris Agreement. Comparison between the results at two warming levels showed the effectiveness of emission reductions in ameliorating the impacts and revealed unavoidable impacts for which adaptation options are urgently needed in the NHL ecosystems.

    « less
  4. Abstract. Thaw and release of permafrost carbon (C) due to climate change is likely tooffset increased vegetation C uptake in northern high-latitude (NHL)terrestrial ecosystems. Models project that this permafrost C feedback mayact as a slow leak, in which case detection and attribution of the feedbackmay be difficult. The formation of talik, a subsurface layer of perenniallythawed soil, can accelerate permafrost degradation and soil respiration,ultimately shifting the C balance of permafrost-affected ecosystems fromlong-term C sinks to long-term C sources. It is imperative to understand andcharacterize mechanistic links between talik, permafrost thaw, andrespiration of deep soil C to detect and quantify the permafrost C feedback.Here, we use the Community Land Model (CLM) version 4.5, a permafrost andbiogeochemistry model, in comparison to long-term deep borehole data alongNorth American and Siberian transects, to investigate thaw-driven C sourcesin NHL (>55N) from 2000 to 2300. Widespread talik at depth isprojected across most of the NHL permafrost region(14million km2) by 2300, 6.2million km2 of which isprojected to become a long-term C source, emitting 10Pg C by 2100,50Pg C by 2200, and 120Pg C by 2300, with few signs ofslowing. Roughly half of the projected C source region is in predominantlywarm sub-Arctic permafrost following talik onset. This region emits only20Pg C by 2300, butmore »the CLM4.5 estimate may be biased low by notaccounting for deep C in yedoma. Accelerated decomposition of deep soilC following talik onset shifts the ecosystem C balance away from surfacedominant processes (photosynthesis and litter respiration), butsink-to-source transition dates are delayed by 20–200 years by highecosystem productivity, such that talik peaks early (2050s, although boreholedata suggest sooner) and C source transition peaks late(2150–2200). The remaining C source region in cold northern Arcticpermafrost, which shifts to a net source early (late 21st century), emits5 times more C (95Pg C) by 2300, and prior to talik formation dueto the high decomposition rates of shallow, young C in organic-rich soilscoupled with low productivity. Our results provide important clues signalingimminent talik onset and C source transition, including (1) late cold-season(January–February) soil warming at depth (2m),(2) increasing cold-season emissions (November–April), and (3) enhancedrespiration of deep, old C in warm permafrost and young, shallow C in organic-rich cold permafrost soils. Our results suggest a mosaic of processes thatgovern carbon source-to-sink transitions at high latitudes and emphasize theurgency of monitoring soil thermal profiles, organic C age and content, cold-season CO2 emissions, andatmospheric 14CO2 as key indicatorsof the permafrost C feedback.

    « less
  5. Abstract Wildfires are a major disturbance to forest carbon (C) balance through both immediate combustion emissions and post-fire ecosystem dynamics. Here we used a process-based biogeochemistry model, the Terrestrial Ecosystem Model (TEM), to simulate C budget in Alaska and Canada during 1986–2016, as impacted by fire disturbances. We extracted the data of difference Normalized Burn Ratio (dNBR) for fires from Landsat TM/ETM imagery and estimated the proportion of vegetation and soil C combustion. We observed that the region was a C source of 2.74 Pg C during the 31-year period. The observed C loss, 57.1 Tg C year −1 , was attributed to fire emissions, overwhelming the net ecosystem production (1.9 Tg C year −1 ) in the region. Our simulated direct emissions for Alaska and Canada are within the range of field measurements and other model estimates. As burn severity increased, combustion emission tended to switch from vegetation origin towards soil origin. When dNBR is below 300, fires increase soil temperature and decrease soil moisture and thus, enhance soil respiration. However, the post-fire soil respiration decreases for moderate or high burn severity. The proportion of post-fire soil emission in total emissions increased with burn severity. Net nitrogen mineralization graduallymore »recovered after fire, enhancing net primary production. Net ecosystem production recovered fast under higher burn severities. The impact of fire disturbance on the C balance of northern ecosystems and the associated uncertainties can be better characterized with long-term, prior-, during- and post-disturbance data across the geospatial spectrum. Our findings suggest that the regional source of carbon to the atmosphere will persist if the observed forest wildfire occurrence and severity continues into the future.« less