scholarly journals Divergent predictions of carbon storage between two global land models: attribution of the causes through traceability analysis

2016 ◽  
Vol 7 (3) ◽  
pp. 649-658 ◽  
Author(s):  
Rashid Rafique ◽  
Jianyang Xia ◽  
Oleksandra Hararuk ◽  
Ghassem R. Asrar ◽  
Guoyong Leng ◽  
...  

Abstract. Representations of the terrestrial carbon cycle in land models are becoming increasingly complex. It is crucial to develop approaches for critical assessment of the complex model properties in order to understand key factors contributing to models' performance. In this study, we applied a traceability analysis which decomposes carbon cycle models into traceable components, for two global land models (CABLE and CLM-CASA′) to diagnose the causes of their differences in simulating ecosystem carbon storage capacity. Driven with similar forcing data, CLM-CASA′ predicted  ∼ 31 % larger carbon storage capacity than CABLE. Since ecosystem carbon storage capacity is a product of net primary productivity (NPP) and ecosystem residence time (τE), the predicted difference in the storage capacity between the two models results from differences in either NPP or τE or both. Our analysis showed that CLM-CASA′ simulated 37 % higher NPP than CABLE. On the other hand, τE, which was a function of the baseline carbon residence time (τ′E) and environmental effect on carbon residence time, was on average 11 years longer in CABLE than CLM-CASA′. This difference in τE was mainly caused by longer τ′E of woody biomass (23 vs. 14 years in CLM-CASA′), and higher proportion of NPP allocated to woody biomass (23 vs. 16 %). Differences in environmental effects on carbon residence times had smaller influences on differences in ecosystem carbon storage capacities compared to differences in NPP and τ′E. Overall, the traceability analysis showed that the major causes of different carbon storage estimations were found to be parameters setting related to carbon input and baseline carbon residence times between two models.

2015 ◽  
Vol 6 (2) ◽  
pp. 1579-1604
Author(s):  
R. Rafique ◽  
J. Xia ◽  
O. Hararuk ◽  
G. Asrar ◽  
Y. Wang ◽  
...  

Abstract. Representations of the terrestrial carbon cycle in land models are becoming increasingly complex. It is crucial to develop approaches for critical assessment of the complex model properties in order to understand key factors contributing to models' performance. In this study, we applied a traceability analysis, which decomposes carbon cycle models into traceable components, to two global land models (CABLE and CLM-CASA') to diagnose the causes of their differences in simulating ecosystem carbon storage capacity. Driven with similar forcing data, the CLM-CASA' model predicted ~ 31 % larger carbon storage capacity than the CABLE model. Since ecosystem carbon storage capacity is a product of net primary productivity (NPP) and ecosystem residence time (τE), the predicted difference in the storage capacity between the two models results from differences in either NPP or τE or both. Our analysis showed that CLM-CASA' simulated 37 % higher NPP than CABLE due to higher rates of carboxylation (Vcmax) in CLM-CASA'. On the other hand, τE, which was a function the baseline carbon residence time (τ'E) and environmental effect on carbon residence time, was on average 11 years longer in CABLE than CLM-CASA'. The difference in τE was mainly found to be caused by longer τ'E in CABLE than CLM-CASA'. This difference in τE was mainly caused by longer τ'E of woody biomass (23 vs. 14 years in CLM-CASA') and higher proportion of NPP allocated to woody biomass (23 vs. 16 %). Differences in environmental effects on carbon residence times had smaller influences on differences in ecosystem carbon storage capacities compared to differences in NPP and τ'E. Overall; the traceability analysis is an effective method for identifying sources of variations between the two models.


2012 ◽  
Vol 5 (5) ◽  
pp. 1259-1271 ◽  
Author(s):  
J. Y. Xia ◽  
Y. Q. Luo ◽  
Y.-P. Wang ◽  
E. S. Weng ◽  
O. Hararuk

Abstract. The spin-up of land models to steady state of coupled carbon–nitrogen processes is computationally so costly that it becomes a bottleneck issue for global analysis. In this study, we introduced a semi-analytical solution (SAS) for the spin-up issue. SAS is fundamentally based on the analytic solution to a set of equations that describe carbon transfers within ecosystems over time. SAS is implemented by three steps: (1) having an initial spin-up with prior pool-size values until net primary productivity (NPP) reaches stabilization, (2) calculating quasi-steady-state pool sizes by letting fluxes of the equations equal zero, and (3) having a final spin-up to meet the criterion of steady state. Step 2 is enabled by averaged time-varying variables over one period of repeated driving forcings. SAS was applied to both site-level and global scale spin-up of the Australian Community Atmosphere Biosphere Land Exchange (CABLE) model. For the carbon-cycle-only simulations, SAS saved 95.7% and 92.4% of computational time for site-level and global spin-up, respectively, in comparison with the traditional method (a long-term iterative simulation to achieve the steady states of variables). For the carbon–nitrogen coupled simulations, SAS reduced computational cost by 84.5% and 86.6% for site-level and global spin-up, respectively. The estimated steady-state pool sizes represent the ecosystem carbon storage capacity, which was 12.1 kg C m−2 with the coupled carbon–nitrogen global model, 14.6% lower than that with the carbon-only model. The nitrogen down-regulation in modeled carbon storage is partly due to the 4.6% decrease in carbon influx (i.e., net primary productivity) and partly due to the 10.5% reduction in residence times. This steady-state analysis accelerated by the SAS method can facilitate comparative studies of structural differences in determining the ecosystem carbon storage capacity among biogeochemical models. Overall, the computational efficiency of SAS potentially permits many global analyses that are impossible with the traditional spin-up methods, such as ensemble analysis of land models against parameter variations.


2012 ◽  
Vol 5 (2) ◽  
pp. 803-836 ◽  
Author(s):  
J. Xia ◽  
Y. Luo ◽  
Y.-P. Wang ◽  
E. Weng ◽  
O. Hararuk

Abstract. The spin-up of land models to steady state of coupled carbon-nitrogen processes is computationally so costly that it becomes a~bottleneck issue for global analysis. In this study, we introduced a semi-analytical solution (SAS) for the spin-up issue. SAS is fundamentally based on the analytic solution to a set of equations that describe carbon transfers within ecosystems over time. SAS is implemented by three steps: (1) having an initial spin-up with prior pool-size values until net primary productivity (NPP) reaches steady state, (2) calculating quasi steady-state pool sizes by letting fluxes of the equations equal zero, and (3) having a final spin-up to meet the criterion of steady state. Step 2 is enabled by averaged time-varying variables over one period of repeated driving forcings. SAS was applied to both site-level and global scale spin-up of the Australian Community Atmosphere Biosphere Land Exchange (CABLE) model. For the carbon-cycle-only simulations, SAS saved 95.7% and 92.4% of computational time for site-level and global spin-up, respectively, in comparison with the traditional method. For the carbon-nitrogen-coupled simulations, SAS reduced computational cost by 84.5% and 86.6% for site-level and global spin-up, respectively. The estimated steady-state pool sizes represent the ecosystem carbon storage capacity, which was 12.1 kg C m−2 with the coupled carbon-nitrogen global model, 14.6% lower than that with the carbon-only model. The nitrogen down-regulation in modeled carbon storage is partly due to the 4.6% decrease in carbon influx (i.e., net primary productivity) and partly due to the 10.5% reduction in residence times. This steady-state analysis accelerated by the SAS method can facilitate comparative studies of structural differences in determining the ecosystem carbon storage capacity among biogeochemical models. Overall, the computational efficiency of SAS potentially permits many global analyses that are impossible with the traditional spin-up methods, such as ensemble analysis of land models against parameter variations.


2018 ◽  
Vol 31 (7) ◽  
pp. 2833-2851 ◽  
Author(s):  
Sha Zhou ◽  
Junyi Liang ◽  
Xingjie Lu ◽  
Qianyu Li ◽  
Lifen Jiang ◽  
...  

Terrestrial carbon cycle models have incorporated increasingly more processes as a means to achieve more-realistic representations of ecosystem carbon cycling. Despite this, there are large across-model variations in the simulation and projection of carbon cycling. Several model intercomparison projects (MIPs), for example, the fifth phase of the Coupled Model Intercomparison Project (CMIP5) (historical simulations), Trends in Net Land–Atmosphere Carbon Exchange (TRENDY), and Multiscale Synthesis and Terrestrial Model Intercomparison Project (MsTMIP), have sought to understand intermodel differences. In this study, the authors developed a suite of new techniques to conduct post-MIP analysis to gain insights into uncertainty sources across 25 models in the three MIPs. First, terrestrial carbon storage dynamics were characterized by a three-dimensional (3D) model output space with coordinates of carbon residence time, net primary productivity (NPP), and carbon storage potential. The latter represents the potential of an ecosystem to lose or gain carbon. This space can be used to measure how and why model output differs. Models with a nitrogen cycle generally exhibit lower annual NPP in comparison with other models, and mostly negative carbon storage potential. Second, a transient traceability framework was used to decompose any given carbon cycle model into traceable components and identify the sources of model differences. The carbon residence time (or NPP) was traced to baseline carbon residence time (or baseline NPP related to the maximum carbon input), environmental scalars, and climate forcing. Third, by applying a variance decomposition method, the authors show that the intermodel differences in carbon storage can be mainly attributed to the baseline carbon residence time and baseline NPP (>90% in the three MIPs). The three techniques developed in this study offer a novel approach to gain more insight from existing MIPs and can point out directions for future MIPs. Since this study is conducted at the global scale for an overview on intermodel differences, future studies should focus more on regional analysis to identify the sources of uncertainties and improve models at the specified mechanism level.


2004 ◽  
Vol 359 (1443) ◽  
pp. 463-476 ◽  
Author(s):  
Jeffrey Q. Chambers ◽  
Whendee L. Silver

Atmospheric changes that may affect physiological and biogeochemical processes in old–growth tropical forests include: (i) rising atmospheric CO 2 concentration; (ii) an increase in land surface temperature; (iii) changes in precipitation and ecosystem moisture status; and (iv) altered disturbance regimes. Elevated CO 2 is likely to directly influence numerous leaf–level physiological processes, but whether these changes are ultimately reflected in altered ecosystem carbon storage is unclear. The net primary productivity (NPP) response of old–growth tropical forests to elevated CO 2 is unknown, but unlikely to exceed the maximum experimentally measured 25% increase in NPP with a doubling of atmospheric CO 2 from pre–industrial levels. In addition, evolutionary constraints exhibited by tropical plants adapted to low CO 2 levels during most of the Late Pleistocene, may result in little response to increased carbon availability. To set a maximum potential response for a Central Amazon forest, using an individual–tree–based carbon cycling model, a modelling experiment was performed constituting a 25% increase in tree growth rate, linked to the known and expected increase in atmospheric CO 2 . Results demonstrated a maximum carbon sequestration rate of ca . 0.2 Mg C per hectare per year (ha −1 yr −1 , where 1 ha = 10 4 m 2 ), and a sequestration rate of only 0.05 Mg C ha −1 yr −1 for an interval centred on calendar years 1980–2020. This low rate results from slow growing trees and the long residence time of carbon in woody tissues. By contrast, changes in disturbance frequency, precipitation patterns and other environmental factors can cause marked and relatively rapid shifts in ecosystem carbon storage. It is our view that observed changes in tropical forest inventory plots over the past few decades is more probably being driven by changes in disturbance or other environmental factors, than by a response to elevated CO 2 . Whether these observed changes in tropical forests are the beginning of long–term permanent shifts or a transient response is uncertain and remains an important research priority.


2012 ◽  
Vol 117 (G3) ◽  
pp. n/a-n/a ◽  
Author(s):  
Ensheng Weng ◽  
Yiqi Luo ◽  
Weile Wang ◽  
Han Wang ◽  
Daniel J. Hayes ◽  
...  

2017 ◽  
Vol 122 (3) ◽  
pp. 528-545 ◽  
Author(s):  
T. L. Smallman ◽  
J.-F. Exbrayat ◽  
M. Mencuccini ◽  
A. A. Bloom ◽  
M. Williams

2010 ◽  
Vol 7 (3) ◽  
pp. 3735-3763 ◽  
Author(s):  
K. Fenn ◽  
Y. Malhi ◽  
M. Morecroft ◽  
C. Lloyd ◽  
M. Thomas

Abstract. There exist very few comprehensive descriptions of the productivity and carbon cycling of forest ecosystems. Here we present a description of the components of annual Net Primary Productivity (NPP), Gross Primary Productivity (GPP), autotrophic and heterotrophic respiration, and ecosystem respiration (RECO) for a temperate mixed deciduous woodland at Wytham Woods in southern Britain, calculated using "bottom-up" biometric and chamber measurements (leaf and wood production and soil and stem respiration). These are compared with estimates of these parameters from eddy-covariance measurements made at the same site. NPP was estimated as 7.0±0.8 Mg C ha−1 yr−1, and GPP as 20.3+1.0 Mg C ha−1 yr−1, a value which closely matched to eddy covariance-derived GPP value of 21.1 Mg C ha−1 yr−1. Annual RECO was calculated as 18.9±1.7 Mg C ha−1 yr−1, close to the eddy covariance value of 19.8 Mg C ha−1 yr−1; the seasonal cycle of biometric and eddy covariance RECO estimates also closely matched. The consistency between eddy covariance and biometric measurements substantially strengthens the confidence we attach to each as alternative indicators of site carbon dynamics, and permits an integrated perspective of the ecosystem carbon cycle. 37% of NPP was allocated below ground, and the ecosystem carbon use efficiency (CUE, = NPP/GPP) calculated to be 0.35±0.05, lower than reported for many temperate broadleaved sites.


Sign in / Sign up

Export Citation Format

Share Document