scholarly journals Technical Note: An efficient method for accelerating the spin-up process for process-based biogeochemistry models

2018 ◽  
Author(s):  
Yang Qu ◽  
Shamil Maksyutov ◽  
Qianlai Zhuang

Abstract. To better understand the role of terrestrial ecosystems in the global carbon cycle and their feedbacks to the global climate system, process-based biogeochemistry models need to be improved with respect to model parameterization and model structure. To achieve these improvements, the spin-up time for those differential equation-based models needs to be shortened. Here, an algorithm for a fast spin-up was developed and implemented in a biogeochemistry model, the Terrestrial Ecosystem Model (TEM). With the new spin-up algorithm, we showed that the model reached a steady state in less than 10 years of computing time, while the original method requires more than 200 years on average of model run. For the test sites with five different plant function types, the new method saves over 90 % of the original spin-up time in site-level simulations. In North America simulations, average spin-up time saving for all grid cells is 85 % for either daily or monthly version of TEM. The developed spin-up method shall greatly facilitate our future quantification of carbon dynamics at fine spatial and temporal scales.

2018 ◽  
Vol 15 (13) ◽  
pp. 3967-3973 ◽  
Author(s):  
Yang Qu ◽  
Shamil Maksyutov ◽  
Qianlai Zhuang

Abstract. To better understand the role of terrestrial ecosystems in the global carbon cycle and their feedbacks to the global climate system, process-based biogeochemistry models need to be improved with respect to model parameterization and model structure. To achieve these improvements, the spin-up time for those differential equation-based models needs to be shortened. Here, an algorithm for a fast spin-up was developed by finding the exact solution of a linearized system representing the cyclo-stationary state of a model and implemented in a biogeochemistry model, the Terrestrial Ecosystem Model (TEM). With the new spin-up algorithm, we showed that the model reached a steady state in less than 10 years of computing time, while the original method requires more than 200 years on average of model run. For the test sites with five different plant functional types, the new method saves over 90 % of the original spin-up time in site-level simulations. In North American simulations, average spin-up time savings for all grid cells is 85 % for either the daily or monthly version of TEM. The developed spin-up method shall be used for future quantification of carbon dynamics at fine spatial and temporal scales.


2011 ◽  
Vol 11 (2) ◽  
pp. 5379-5405 ◽  
Author(s):  
P. K. Patra ◽  
Y. Niwa ◽  
T. J. Schuck ◽  
C. A. M. Brenninkmeijer ◽  
T. Machida ◽  
...  

Abstract. Quantifying the fluxes of carbon dioxide (CO2) between the atmosphere and terrestrial ecosystems in all their diversity, across the continents, is important and urgent for implementing effective mitigating policies. Whereas much is known for Europe and North America for instance, in comparison, South Asia, with 1.6 billion inhabitants and considerable CO2 fluxes, remained terra incognita in this respect. We use regional measurements of atmospheric CO2 aboard a Lufthansa passenger aircraft between Frankfurt (Germany) and Chennai (India) at cruise altitude, in addition to the existing network sites for 2008, to estimate monthly fluxes for 64-regions using Bayesian inversion and transport model simulations. The applicability of the model's transport parameterization is confirmed using SF6, CH4 and N2O simulations for the CARIBIC datasets. The annual carbon flux obtained by including the aircraft data is twice as large as the fluxes simulated by a terrestrial ecosystem model that was applied to prescribe the fluxes used in the inversions. It is shown that South Asia sequestered carbon at a rate of 0.37±0.20 Pg C yr−1 (1Pg C = 1015 g of carbon in CO2) for the years 2007 and 2008. The seasonality and the strength of the calculated monthly fluxes are successfully validated using independent measurements of vertical CO2 profiles over Delhi and spatial variations at cruising altitude over Asia aboard Japan Airlines passenger aircraft.


2015 ◽  
Vol 19 (16) ◽  
pp. 1-21 ◽  
Author(s):  
Chang Liao ◽  
Qianlai Zhuang

Abstract Droughts dramatically affect plant production of global terrestrial ecosystems. To date, quantification of this impact remains a challenge because of the complex plant physiological and biochemical processes associated with drought. Here, this study incorporates a drought index into an existing process-based terrestrial ecosystem model to estimate the drought impact on global plant production for the period 2001–10. Global Moderate Resolution Imaging Spectroradiometer (MODIS) gross primary production (GPP) data products are used to constrain model parameters and verify the model algorithms. The verified model is then applied to evaluate the drought impact. The study indicates that droughts will reduce GPP by 9.8 g C m−2 month−1 during the study period. On average, drought reduces GPP by 10% globally. As a result, the global GPP decreased from 106.4 to 95.9 Pg C yr−1 while the global net primary production (NPP) decreased from 54.9 to 49.9 Pg C yr−1. This study revises the estimation of the global NPP and suggests that the future quantification of the global carbon budget of terrestrial ecosystems should take the drought impact into account.


2011 ◽  
Vol 11 (9) ◽  
pp. 4163-4175 ◽  
Author(s):  
P. K. Patra ◽  
Y. Niwa ◽  
T. J. Schuck ◽  
C. A. M. Brenninkmeijer ◽  
T. Machida ◽  
...  

Abstract. Quantifying the fluxes of carbon dioxide (CO2) between the atmosphere and terrestrial ecosystems in all their diversity, across the continents, is important and urgent for implementing effective mitigating policies. Whereas much is known for Europe and North America for instance, in comparison, South Asia, with 1.6 billion inhabitants and considerable CO2 fluxes, remained terra incognita in this respect. We use regional measurements of atmospheric CO2 aboard a Lufthansa passenger aircraft between Frankfurt (Germany) and Chennai (India) at cruise altitude, in addition to the existing network sites for 2008, to estimate monthly fluxes for 64-regions using Bayesian inversion and transport model simulations. The applicability of the model's transport parameterization is confirmed using SF6, CH4 and N2O simulations for the CARIBIC datasets. The annual amplitude of carbon flux obtained by including the aircraft data is twice as large as the fluxes simulated by a terrestrial ecosystem model that was applied to prescribe the fluxes used in the inversions. It is shown that South Asia sequestered carbon at a rate of 0.37 ± 0.20 Pg C yr−1 (1 Pg C = 1015 g of carbon in CO2) for the years 2007 and 2008. The seasonality and the strength of the calculated monthly fluxes are successfully validated using independent measurements of vertical CO2 profiles over Delhi and spatial variations at cruising altitude over Asia aboard Japan Airlines passenger aircraft.


2021 ◽  
Vol 18 (23) ◽  
pp. 6245-6269
Author(s):  
Junrong Zha ◽  
Qianlai Zhuang

Abstract. Mosses are ubiquitous in northern terrestrial ecosystems, and play an important role in regional carbon, water and energy cycling. Current global land surface models that do not consider mosses may bias the quantification of regional carbon dynamics. Here we incorporate mosses as a new plant functional type into the process-based Terrestrial Ecosystem Model (TEM 5.0), to develop a new model (TEM_Moss). The new model explicitly quantifies the interactions between vascular plants and mosses and their competition for energy, water, and nutrients. Compared to the estimates using TEM 5.0, the new model estimates that the regional terrestrial soils currently store 132.7 Pg more C and will store 157.5 and 179.1 Pg more C under the RCP8.5 and RCP2.6 scenarios, respectively, by the end of the 21st century. Ensemble regional simulations forced with different parameters for the 21st century with TEM_Moss predict that the region will accumulate 161.1±142.1 Pg C under the RCP2.6 scenario and 186.7±166.1 Pg C under the RCP8.5 scenario over the century. Our study highlights the necessity of coupling moss into Earth system models to adequately quantify terrestrial carbon–climate feedbacks in the Arctic.


2021 ◽  
Vol 2 (1) ◽  
Author(s):  
Yuhao Feng ◽  
Haojie Su ◽  
Zhiyao Tang ◽  
Shaopeng Wang ◽  
Xia Zhao ◽  
...  

AbstractGlobal climate change likely alters the structure and function of vegetation and the stability of terrestrial ecosystems. It is therefore important to assess the factors controlling ecosystem resilience from local to global scales. Here we assess terrestrial vegetation resilience over the past 35 years using early warning indicators calculated from normalized difference vegetation index data. On a local scale we find that climate change reduced the resilience of ecosystems in 64.5% of the global terrestrial vegetated area. Temperature had a greater influence on vegetation resilience than precipitation, while climate mean state had a greater influence than climate variability. However, there is no evidence for decreased ecological resilience on larger scales. Instead, climate warming increased spatial asynchrony of vegetation which buffered the global-scale impacts on resilience. We suggest that the response of terrestrial ecosystem resilience to global climate change is scale-dependent and influenced by spatial asynchrony on the global scale.


2021 ◽  
Author(s):  
Franziska Lechleitner ◽  
Christopher C. Day ◽  
Oliver Kost ◽  
Micah Wilhelm ◽  
Negar Haghipour ◽  
...  

<p>Terrestrial ecosystems are intimately linked with the global climate system, but their response to ongoing and future anthropogenic climate change remains poorly understood. Reconstructing the response of terrestrial ecosystem processes over past periods of rapid and substantial climate change can serve as a tool to better constrain the sensitivity in the ecosystem-climate response.</p><p>In this talk, we will present a new reconstruction of soil respiration in the temperate region of Western Europe based on speleothem carbon isotopes (δ<sup>13</sup>C). Soil respiration remains poorly constrained over past climatic transitions, but is critical for understanding the global carbon cycle and its response to ongoing anthropogenic warming. Our study builds upon two decades of speleothem research in Western Europe, which has shown clear correlation between δ<sup>13</sup>C and regional temperature reconstructions during the last glacial and the deglaciation, with exceptional regional coherency in timing, amplitude, and absolute δ<sup>13</sup>C variation. By combining innovative multi-proxy geochemical analysis (δ<sup>13</sup>C, Ca isotopes, and radiocarbon) on three speleothems from Northern Spain, and quantitative forward modelling of processes in soil, karst, and cave, we show how deglacial variability in speleothem δ<sup>13</sup>C is best explained by increasing soil respiration. Our study is the first to quantify and remove the effects of prior calcite precipitation (PCP, using Ca isotopes) and bedrock dissolution (open vs closed system, using the radiocarbon reservoir effect) from the speleothem δ<sup>13</sup>C signal to derive changes in respired δ<sup>13</sup>C over time. Our approach allows us to estimate the temperature sensitivity of soil respiration (Q<sub>10</sub>), which is higher than current measurements, suggesting that part of the speleothem signal may be related to a change in the composition of the soil respired δ<sup>13</sup>C. This is likely related to changing substrate through increasing contribution from vegetation biomass with the onset of the Holocene.</p><p>These results highlight the exciting possibilities speleothems offer as a coupled archive for quantitative proxy-based reconstructions of climate and ecosystem conditions.</p>


2017 ◽  
Author(s):  
Joe R. Melton ◽  
Reinel Sospedra-Alfonso ◽  
Kelly E. McCusker

Abstract. We investigate the application of clustering algorithms to represent sub-grid scale variability in soil texture for use in a global-scale terrestrial ecosystem model. Our model, the coupled Canadian Land Surface Scheme – Canadian Terrestrial Ecosystem Model (CLASS-CTEM), is typically implemented at a coarse spatial resolution (ca. 2.8° × 2.8°) due to its use as the land surface component of the Canadian Earth System Model (CanESM). CLASS-CTEM can, however, be run with tiling of the land surface as a means to represent sub-grid heterogeneity. We first determined that the model was sensitive to tiling of the soil textures via an idealized test case before attempting to cluster soil textures globally. To cluster a high-resolution soil texture dataset onto our coarse model grid, we use two linked algorithms (OPTICS (Ankerst et al., 1999; Daszykowski et al., 2002) and Sander et al. (2003)) to provide tiles of representative soil textures for use as CLASS-CTEM inputs. The clustering process results in, on average, about three tiles per CLASS-CTEM grid cell with most cells having four or less tiles. Results from CLASS-CTEM simulations conducted with the tiled inputs (Cluster) versus those using a simple grid-mean soil texture (Gridmean) show CLASS-CTEM, at least on a global scale, is relatively insensitive to the tiled soil textures, however differences can be large in arid or peatland regions. The Cluster simulation has generally lower soil moisture and lower overall vegetation productivity than the Gridmean simulation except in arid regions where plant productivity increases. In these dry regions, the influence of the tiling is stronger due to the general state of vegetation moisture stress which allows a single tile, whose soil texture retains more plant available water, to yield much higher productivity. Although the use of clustering analysis appears promising as a means to represent sub-grid heterogeneity, soil textures appear to be reasonably represented for global scale simulations using a simple grid-mean value.


2019 ◽  
Author(s):  
Tea Thum ◽  
Silvia Caldararu ◽  
Jan Engel ◽  
Melanie Kern ◽  
Marleen Pallandt ◽  
...  

Abstract. The dynamics of terrestrial ecosystems are shaped by the coupled cycles of carbon, nitrogen and phosphorus, and strongly depend on the availability of water and energy. These interactions shape future terrestrial biosphere responses to global change. Many process-based models of the terrestrial biosphere have been gradually extended from considering carbon-water interactions to also including nitrogen, and later, phosphorus dynamics. This evolutionary model development has hindered full integration of these biogeochemical cycles and the feedbacks amongst them. Here we present a new terrestrial ecosystem model QUINCY (QUantifying Interactions between terrestrial Nutrient CYcles and the climate system), which is formulated around a consistent representation of element cycling in terrestrial ecosystems. This new model includes i) a representation of plant growth which separates source (e.g. photosynthesis) and sink (growth rate of individual tissues, constrained by nutrients, temperature, and water availability) processes; ii) the acclimation of many ecophysiological processes to meteorological conditions and/or nutrient availabilities; iii) an explicit representation of vertical soil processes to separate litter and soil organic matter dynamics; iv) a range of new diagnostics (leaf chlorophyll content; 13C, 14C, and 15N isotope tracers) to allow for a more in-depth model evaluation. We present the model structure and provide an assessment of its performance against a range of observations from global-scale ecosystem monitoring networks. We demonstrate that the framework is capable of consistently simulating ecosystem dynamics across a large gradient in climate and soil conditions, as well as across different plant functional types. To aid this understanding we provide an assessment of the model's sensitivity to its parameterisation and the associated uncertainty.


Sign in / Sign up

Export Citation Format

Share Document