skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Attention:

The NSF Public Access Repository (PAR) system and access will be unavailable from 10:00 PM ET on Friday, February 6 until 10:00 AM ET on Saturday, February 7 due to maintenance. We apologize for the inconvenience.


Title: Optimal model complexity for terrestrial carbon cycle prediction
Abstract. The terrestrial carbon cycle plays a critical role in modulating the interactions of climate with the Earth system, but different models often make vastly different predictions of its behavior. Efforts to reduce model uncertainty have commonly focused on model structure, namely by introducing additional processes and increasing structural complexity. However, the extent to which increased structural complexity can directly improve predictive skill is unclear. While adding processes may improve realism, the resulting models are often encumbered by a greater number of poorly determined or over-generalized parameters. To guide efficient model development, here we map the theoretical relationship between model complexity and predictive skill. To do so, we developed 16 structurally distinct carbon cycle models spanning an axis of complexity and incorporated them into a model–data fusion system. We calibrated each model at six globally distributed eddy covariance sites with long observation time series and under 42 data scenarios that resulted in different degrees of parameter uncertainty. For each combination of site, data scenario, and model, we then predicted net ecosystem exchange (NEE) and leaf area index (LAI) for validation against independent local site data. Though the maximum model complexity we evaluated is lower than most traditional terrestrial biosphere models, the complexity range we explored provides universal insight into the inter-relationship between structural uncertainty, parametric uncertainty, and model forecast skill. Specifically, increased complexity only improves forecast skill if parameters are adequately informed (e.g., when NEE observations are used for calibration). Otherwise, increased complexity can degrade skill and an intermediate-complexity model is optimal. This finding remains consistent regardless of whether NEE or LAI is predicted. Our COMPLexity EXperiment (COMPLEX) highlights the importance of robust observation-based parameterization for land surface modeling and suggests that data characterizing net carbon fluxes will be key to improving decadal predictions of high-dimensional terrestrial biosphere models.  more » « less
Award ID(s):
1942133
PAR ID:
10314190
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Biogeosciences
Volume:
18
Issue:
8
ISSN:
1726-4189
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract A significant warming effect on arctic tundra is greening. Although this increase in predominantly woody vegetation has been linked to increases in gross primary productivity, increasing temperatures also stimulate ecosystem respiration. We present a novel analysis from small-scale plot measurements showing that the shape of the temperature- and light-dependent sink-to-source threshold (where net ecosystem exchange (NEE) equals zero) differs between two tussock tundra ecosystems differing in leaf area index (LAI). At the higher LAI site, the threshold is exceeded (i.e the ecosystem becomes a source) at relatively higher temperatures under low light but at lower temperatures under high light. At the lower LAI site, the threshold is exceeded at relatively lower temperatures under low light but at higher temperatures under high light. We confirmed this response at a single site where LAI was experimentally increased. This suggests the carbon balance of the tundra may be sensitive to small increases in temperature under low light, but that this effect may be significantly offset by increases in LAI. Importantly, we found that this LAI effect is reversed under high light, and so in a warming tundra, greater vegetation cover could have a progressively negative effect on net carbon uptake. 
    more » « less
  2. Abstract Climate change is having significant impacts on Earth’s ecosystems and carbon budgets, and in the Arctic may drive a shift from an historic carbon sink to a source. Large uncertainties in terrestrial biosphere models (TBMs) used to forecast Arctic changes demonstrate the challenges of determining the timing and extent of this possible switch. This spread in model predictions can limit the ability of TBMs to guide management and policy decisions. One of the most influential sources of model uncertainty is model parameterization. Parameter uncertainty results in part from a mismatch between available data in databases and model needs. We identify that mismatch for three TBMs, DVM-DOS-TEM, SIPNET and ED2, and four databases with information on Arctic and boreal above- and belowground traits that may be applied to model parametrization. However, focusing solely on such data gaps can introduce biases towards simple models and ignores structural model uncertainty, another main source for model uncertainty. Therefore, we develop a causal loop diagram (CLD) of the Arctic and boreal ecosystem that includes unquantified, and thus unmodeled, processes. We map model parameters to processes in the CLD and assess parameter vulnerability via the internal network structure. One important substructure, feed forward loops (FFLs), describe processes that are linked both directly and indirectly. When the model parameters are data-informed, these indirect processes might be implicitly included in the model, but if not, they have the potential to introduce significant model uncertainty. We find that the parameters describing the impact of local temperature on microbial activity are associated with a particularly high number of FFLs but are not constrained well by existing data. By employing ecological models of varying complexity, databases, and network methods, we identify the key parameters responsible for limited model accuracy. They should be prioritized for future data sampling to reduce model uncertainty. 
    more » « less
  3. Improving our ability to understand and predict the dynamics of the terrestrial carbon cycle remains a pressing challenge despite a rapidly growing volume and diversity of Earth Observation data. State data assimilation represents a path forward via an iterative cycle of making process-based forecasts and then statistically reconciling these forecasts against numerous ground-based and remotely-sensed data constraints into a “reanalysis” data product that provides full spatiotemporal carbon budgets with robust uncertainty accounting. Here we report on an >100x expansion of the PEcAn+SIPNET reanalysis from 500 sites CONUS, 25 ensemble members, and 2 data constraints to 6400 sites across North America, 100 ensemble members, and 5 data constraints: GEDI and Landtrendr AGB, MODIS LAI, SoilGrids Soil C, and SMAP soil moisture. We also report on an ensemble-based machine learning (ML) downscaling to a 1km product that preserves spatial, temporal, and across-variable covariances and demonstrate the impacts of these covariances on uncertainty accounting (Fig. 1). Synergistically, we use the same ML models to assess what climate, vegetation, and soil variables explain the spatiotemporal variability in different C pools and fluxes. In addition, we review a wide range of ongoing validation activities, comparing the outputs of the reanalysis against withheld data from: Ameriflux and NEON NEE and LE; USFS Forest Inventory biomass, biomass increment, tree rings, soil C, and litter; and NEON soil C and soil respiration. Finally, we touch on ML analyses to diagnose and correct systematic biases and emulator-based recalibration efforts. 
    more » « less
  4. Globally, planted forests are rapidly replacing naturally regenerated stands but the implications for canopy structure, carbon (C) storage, and the linkages between the two are unclear. We investigated the successional dynamics, interlinkages and mechanistic relationships between wood net primary production (NPPw) and canopy structure in planted and naturally regenerated red pine (Pinus resinosa Sol. ex Aiton) stands spanning ≥ 45 years of development. We focused our canopy structural analysis on leaf area index (LAI) and a spatially integrative, terrestrial LiDAR-based complexity measure, canopy rugosity, which is positively correlated with NPPw in several naturally regenerated forests, but which has not been investigated in planted stands. We estimated stand NPPw using a dendrochronological approach and examined whether canopy rugosity relates to light absorption and light–use efficiency. We found that canopy rugosity increased similarly with age in planted and naturally regenerated stands, despite differences in other structural features including LAI and stem density. However, the relationship between canopy rugosity and NPPw was negative in planted and not significant in naturally regenerated stands, indicating structural complexity is not a globally positive driver of NPPw. Underlying the negative NPPw-canopy rugosity relationship in planted stands was a corresponding decline in light-use efficiency, which peaked in the youngest, densely stocked stand with high LAI and low structural complexity. Even with significant differences in the developmental trajectories of canopy structure, NPPw, and light use, planted and naturally regenerated stands stored similar amounts of C in wood over a 45-year period. We conclude that widespread increases in planted forests are likely to affect age-related patterns in canopy structure and NPPw, but planted and naturally regenerated forests may function as comparable long-term C sinks via different structural and mechanistic pathways. 
    more » « less
  5. Abstract Accelerated warming of the Arctic can affect the global climate system by thawing permafrost and exposing organic carbon in soils to decompose and release greenhouse gases into the atmosphere. We used a process-based biosphere model (DVM-DOS-TEM) designed to simulate biophysical and biogeochemical interactions between the soil, vegetation, and atmosphere. We varied soil and environmental parameters to assess the impact on cryohydrological and biogeochemical outputs in the model. We analyzed the responses of ecosystem carbon balances to permafrost thaw by running site-level simulations at two long-term tundra ecological monitoring sites in Alaska: Eight Mile Lake (EML) and Imnavait Creek Watershed (IMN), which are characterized by similar tussock tundra vegetation but differing soil drainage conditions and climate. Model outputs showed agreement with field observations at both sites for soil physical properties and ecosystem CO2fluxes. Model simulations of Net Ecosystem Exchange (NEE) showed an overestimation during the frozen season (higher CO2emissions) at EML with a mean NEE of 26.98 ± 4.83 gC/m2/month compared to observational mean of 22.01 ± 5.67 gC/m2/month, and during the fall months at IMN, with a modeled mean of 19.21 ± 7.49 gC/m2/month compared to observation mean of 11.9 ± 4.45 gC/m2/month. Our results underscore the importance of representing the impact of soil drainage conditions on the thawing of permafrost soils, particularly poorly drained soils, which will drive the magnitude of carbon released at sites across the high-latitude tundra. These findings can help improve predictions of net carbon releases from thawing permafrost, ultimately contributing to a better understanding of the impact of Arctic warming on the global climate system. 
    more » « less