scholarly journals Towards a more objective evaluation of modelled land-carbon trends using atmospheric CO<sub>2</sub> and satellite-based vegetation activity observations

2013 ◽  
Vol 10 (6) ◽  
pp. 4189-4210 ◽  
Author(s):  
D. Dalmonech ◽  
S. Zaehle

Abstract. Terrestrial ecosystem models used for Earth system modelling show a significant divergence in future patterns of ecosystem processes, in particular the net land–atmosphere carbon exchanges, despite a seemingly common behaviour for the contemporary period. An in-depth evaluation of these models is hence of high importance to better understand the reasons for this disagreement. Here, we develop an extension for existing benchmarking systems by making use of the complementary information contained in the observational records of atmospheric CO2 and remotely sensed vegetation activity to provide a novel set of diagnostics of ecosystem responses to climate variability in the last 30 yr at different temporal and spatial scales. The selection of observational characteristics (traits) specifically considers the robustness of information given that the uncertainty of both data and evaluation methodology is largely unknown or difficult to quantify. Based on these considerations, we introduce a baseline benchmark – a minimum test that any model has to pass – to provide a more objective, quantitative evaluation framework. The benchmarking strategy can be used for any land surface model, either driven by observed meteorology or coupled to a climate model. We apply this framework to evaluate the offline version of the MPI Earth System Model's land surface scheme JSBACH. We demonstrate that the complementary use of atmospheric CO2 and satellite-based vegetation activity data allows pinpointing of specific model deficiencies that would not be possible by the sole use of atmospheric CO2 observations.

2012 ◽  
Vol 9 (11) ◽  
pp. 16087-16138 ◽  
Author(s):  
D. Dalmonech ◽  
S. Zaehle

Abstract. Terrestrial ecosystem models used for Earth system modelling show a significant divergence in future patterns of ecosystem processes, in particular carbon exchanges, despite a seemingly common behaviour for the contemporary period. An in-depth evaluation of these models is hence of high importance to achieve a better understanding of the reasons for this disagreement. Here, we develop an extension for existing benchmarking systems by making use of the complementary information contained in the observational records of atmospheric CO2 and remotely-sensed vegetation activity to provide a firm set of diagnostics of ecosystem responses to climate variability in the last 30 yr at different temporal and spatial scales. The selection of observational characteristics (traits) specifically considers the robustness of information given the uncertainties in both data and evaluation analysis. In addition, we provide a baseline benchmark, a minimum test that the model under consideration has to pass, to provide a more objective, quantitative evaluation framework. The benchmarking strategy can be used for any land surface model, either driven by observed meteorology or coupled to a climate model. We apply this framework to evaluate the offline version of the MPI-Earth system model's land surface scheme JSBACH. We demonstrate that the complementary use of atmospheric CO2 and satellite based vegetation activity data allows to pinpoint specific model failures that would not be possible by the sole use of atmospheric CO2 observations.


2014 ◽  
Vol 7 (6) ◽  
pp. 8649-8701 ◽  
Author(s):  
J. Ryder ◽  
J. Polcher ◽  
P. Peylin ◽  
C. Ottlé ◽  
Y. Chen ◽  
...  

Abstract. In Earth system modelling, a description of the energy budget of the vegetated surface layer is fundamental as it determines the meteorological conditions in the planetary boundary layer and as such contributes to the atmospheric conditions and its circulation. The energy budget in most Earth system models has long been based on a "big-leaf approach", with averaging schemes that represent in-canopy processes. Such models have difficulties in reproducing consistently the energy balance in field observations. We here outline a newly developed numerical model for energy budget simulation, as a component of the land surface model ORCHIDEE-CAN (Organising Carbon and Hydrology In Dynamic Ecosystems – CANopy). This new model implements techniques from single-site canopy models in a practical way. It includes representation of in-canopy transport, a multilayer longwave radiation budget, height-specific calculation of aerodynamic and stomatal conductance, and interaction with the bare soil flux within the canopy space. Significantly, it avoids iterations over the height of tha canopy and so maintains implicit coupling to the atmospheric model LMDz. As a first test, the model is evaluated against data from both an intensive measurement campaign and longer term eddy covariance measurements for the intensively studied Eucalyptus stand at Tumbarumba, Australia. The model performs well in replicating both diurnal and annual cycles of fluxes, as well as the gradients of sensible heat fluxes. However, the model overestimates sensible heat flux against an underestimate of the radiation budget. Improved performance is expected through the implementation of a more detailed calculation of stand albedo and a more up-to-date stomatal conductance calculation.


2014 ◽  
Vol 28 (1) ◽  
pp. 272-291 ◽  
Author(s):  
Daniela Dalmonech ◽  
Sönke Zaehle ◽  
Gregor J. Schürmann ◽  
Victor Brovkin ◽  
Christian Reick ◽  
...  

Abstract The capacity of earth system models (ESMs) to make reliable projections of future atmospheric CO2 and climate is strongly dependent on the ability of the land surface model to adequately simulate the land carbon (C) cycle. Defining “adequate” performance of the land model requires an understanding of the contributions of climate model and land model errors to the land C cycle. Here, a benchmarking framework is applied based on significant, observed characteristics of the land C cycle for the contemporary period, for which sufficient evaluation data are available, to test the ability of the JSBACH land surface component of the Max Planck Institute Earth System Model (MPI-ESM) to simulate land C trends. Particular attention is given to the role of potential effects caused by climate biases, and therefore investigation is made of the results of model configurations in which JSBACH is interactively “coupled” to atmosphere and ocean components and of an “uncoupled” configuration, where JSBACH is driven by reconstructed meteorology. The ability of JSBACH to simulate the observed phase of phenology and seasonal C fluxes is not strongly affected by climate biases. Contrarily, noticeable differences in the simulated gross primary productivity and land C stocks emerge between coupled and uncoupled configurations, leading to significant differences in the decadal terrestrial C balance and its sensitivity to climate. These differences are strongly controlled by climate biases of the MPI-ESM, in particular those affecting soil moisture. To effectively characterize model performance, the potential effects of climate biases on the land C dynamics need to be considered during the development and calibration of land surface models.


2020 ◽  
Author(s):  
Dirk Barbi ◽  
Nadine Wieters ◽  
Paul Gierz ◽  
Fatemeh Chegini ◽  
Sara Khosravi ◽  
...  

Abstract. Earth system and climate modelling involves the simulation of processes on a wide range of scales and within and across various components of the Earth system. In practice, component models are often developed independently by different research groups and then combined using a dedicated coupling software. This procedure not only leads to a strongly growing number of available versions of model components and coupled setups but also to model- and system-dependent ways of obtaining and operating them. Therefore, implementing these Earth System Models (ESMs) can be challenging and extremely time-consuming, especially for less experienced modellers, or scientists aiming to use different ESMs as in the case of inter-comparison projects. To assist researchers and modellers by reducing avoidable complexity, we developed the ESM-Tools software, which provides a standard way for downloading, configuring, compiling, running and monitoring different models - coupled ESMs and stand-alone models alike - on a variety of High-Performance Computing (HPC) systems. (The ESM-Tools are equally applicable and helpful for stand-alone as for coupled models. In fact, the ESM-Tools are used as standard compile and runtime infrastructure for FESOM2, and currently also applied for ECHAM and ICON standalone simulations. As coupled ESMs are technically the more challenging tasks, we will focus on coupled setups, always implying that stand-alone models can benefit in the same way.) With the ESM-Tools, the user is only required to provide a short script consisting of only the experiment specific definitions, while the software executes all the phases of a simulation in the correct order. The software, which is well documented and easy to install and use, currently supports four ocean models, three atmosphere models, two biogeochemistry models, an ice sheet model, an isostatic adjustment model, a hydrology model and a land-surface model. ESM-Tools has been entirely re-coded in a high-level programming language (Python) and provides researchers with an even more user-friendly interface for Earth system modelling lately. The ESM-Tools were developed within the framework of the project Advanced Earth System Model Capacity, supported by the Helmholtz Association.


2019 ◽  
Vol 12 (8) ◽  
pp. 3503-3521 ◽  
Author(s):  
Simon P. K. Bowring ◽  
Ronny Lauerwald ◽  
Bertrand Guenet ◽  
Dan Zhu ◽  
Matthieu Guimberteau ◽  
...  

Abstract. Few Earth system models adequately represent the unique permafrost soil biogeochemistry and its respective processes; this significantly contributes to uncertainty in estimating their responses, and that of the planet at large, to warming. Likewise, the riverine component of what is known as the “boundless carbon cycle” is seldom recognised in Earth system modelling. The hydrological mobilisation of organic material from a ∼1330–1580 PgC carbon stock to the river network results in either sedimentary settling or atmospheric “evasion”, processes widely expected to increase with amplified Arctic climate warming. Here, the production, transport, and atmospheric release of dissolved organic carbon (DOC) from high-latitude permafrost soils into inland waters and the ocean are explicitly represented for the first time in the land surface component (ORCHIDEE) of a CMIP6 global climate model (Institut Pierre Simon Laplace – IPSL). The model, ORCHIDEE MICT-LEAK, which represents the merger of previously described ORCHIDEE versions MICT and LEAK, mechanistically represents (a) vegetation and soil physical processes for high-latitude snow, ice, and soil phenomena and (b) the cycling of DOC and CO2, including atmospheric evasion, along the terrestrial–aquatic continuum from soils through the river network to the coast at 0.5 to 2∘ resolution. This paper, the first in a two-part study, presents the rationale for including these processes in a high-latitude-specific land surface model, then describes the model with a focus on novel process implementations, followed by a summary of the model configuration and simulation protocol. The results of these simulation runs, conducted for the Lena River basin, are evaluated against observational data in the second part of this study.


2016 ◽  
Vol 9 (1) ◽  
pp. 223-245 ◽  
Author(s):  
J. Ryder ◽  
J. Polcher ◽  
P. Peylin ◽  
C. Ottlé ◽  
Y. Chen ◽  
...  

Abstract. In Earth system modelling, a description of the energy budget of the vegetated surface layer is fundamental as it determines the meteorological conditions in the planetary boundary layer and as such contributes to the atmospheric conditions and its circulation. The energy budget in most Earth system models has been based on a big-leaf approach, with averaging schemes that represent in-canopy processes. Furthermore, to be stable, that is to say, over large time steps and without large iterations, a surface layer model should be capable of implicit coupling to the atmospheric model. Surface models with large time steps, however, have difficulties in reproducing consistently the energy balance in field observations. Here we outline a newly developed numerical model for energy budget simulation, as a component of the land surface model ORCHIDEE-CAN (Organising Carbon and Hydrology In Dynamic Ecosystems – CANopy). This new model implements techniques from single-site canopy models in a practical way. It includes representation of in-canopy transport, a multi-layer long-wave radiation budget, height-specific calculation of aerodynamic and stomatal conductance, and interaction with the bare-soil flux within the canopy space. Significantly, it avoids iterations over the height of the canopy and so maintains implicit coupling to the atmospheric model LMDz (Laboratoire de Météorologie Dynamique Zoomed model). As a first test, the model is evaluated against data from both an intensive measurement campaign and longer-term eddy-covariance measurements for the intensively studied Eucalyptus stand at Tumbarumba, Australia. The model performs well in replicating both diurnal and annual cycles of energy and water fluxes, as well as the vertical gradients of temperature and of sensible heat fluxes.


2021 ◽  
Vol 14 (6) ◽  
pp. 4051-4067
Author(s):  
Dirk Barbi ◽  
Nadine Wieters ◽  
Paul Gierz ◽  
Miguel Andrés-Martínez ◽  
Deniz Ural ◽  
...  

Abstract. Earth system and climate modelling involves the simulation of processes on a wide range of scales and within and across various compartments of the Earth system. In practice, component models are often developed independently by different research groups, adapted by others to their special interests and then combined using a dedicated coupling software. This procedure not only leads to a strongly growing number of available versions of model components and coupled setups but also to model- and high-performance computing (HPC)-system-dependent ways of obtaining, configuring, building and operating them. Therefore, implementing these Earth system models (ESMs) can be challenging and extremely time consuming, especially for less experienced modellers or scientists aiming to use different ESMs as in the case of intercomparison projects. To assist researchers and modellers by reducing avoidable complexity, we developed the ESM-Tools software, which provides a standard way for downloading, configuring, compiling, running and monitoring different models on a variety of HPC systems. It should be noted that ESM-Tools is not a coupling software itself but a workflow and infrastructure management tool to provide access to increase usability of already existing components and coupled setups. As coupled ESMs are technically the more challenging tasks, we will focus on coupled setups, always implying that stand-alone models can benefit in the same way. With ESM-Tools, the user is only required to provide a short script consisting of only the experiment-specific definitions, while the software executes all the phases of a simulation in the correct order. The software, which is well documented and easy to install and use, currently supports four ocean models, three atmosphere models, two biogeochemistry models, an ice sheet model, an isostatic adjustment model, a hydrology model and a land-surface model. Compared to previous versions, ESM-Tools has lately been entirely recoded in a high-level programming language (Python) and provides researchers with an even more user-friendly interface for Earth system modelling. ESM-Tools was developed within the framework of the Advanced Earth System Model Capacity project, supported by the Helmholtz Association.


2016 ◽  
Author(s):  
G. J. Schürmann ◽  
T. Kaminski ◽  
C. Köstler ◽  
N. Carvalhais ◽  
M. Voßbeck ◽  
...  

Abstract. We describe the Max Planck Institute Carbon Cycle Data Assimilation System (MPI-CCDAS) built around the tangent-linear version of the land surface scheme of the MPI-Earth System Model v1 (JSBACH). The simulated terrestrial biosphere processes (phenology and carbon balance) were constrained by observations of the fraction of photosynthetically active radiation (TIP-FAPAR product) and by observations of atmospheric CO2 at a global set of monitoring stations for the years 2005–2009. The system successfully, and computationally efficiently, improved average foliar area and northern extra-tropical seasonality of foliar area when constrained by TIP-FAPAR. Global net and gross carbon fluxes were improved when constrained by atmospheric CO2, although the system tended to underestimate tropical productivity. Assimilating both data streams jointly allowed the MPI-CCDAS to match both observations (TIP-FAPAR and atmospheric CO2) equally well as the single data stream assimilation cases, therefore overall increasing the appropriateness of the resultant parameter values and biosphere dynamics. Our study thus highlights the role of the TIP-FAPAR product in stabilising the underdetermined atmospheric inversion problem and demonstrates the value of multiple-data stream assimilation for the simulation of terrestrial biosphere dynamics. The constraint on regional gross and net CO2 flux patterns is limited through the parametrisation of the biosphere model. We expect improvement on that aspect through a refined initialisation strategy and inclusion of further biosphere observations as constraints.


2021 ◽  
Author(s):  
Gabriele Arduini ◽  
Ervin Zsoter ◽  
Hannah Cloke ◽  
Elisabeth Stephens ◽  
Christel Prudhomme

&lt;p&gt;Snow processes, with the water stored in the snowpack and released as snowmelt, are very important components of the water balance, in particular in high latitude and mountain regions. The evolution of the snow cover and the timing of the snow melt can have major impact on river discharge. Land surface models are used in Earth System models to compute exchanges of water, energy and momentum between the atmosphere and the surface underneath, and also to compute other components of the hydrological cycle. In order to improve the snow representation, a new multi-layer snow scheme is under development in the HTESSEL land surface model of the European Centre for Medium&amp;#8208;Range Weather Forecasts (ECMWF) Integrated Forecasting System (IFS), to replace the current single-layer snow scheme used in HTESSEL. The new scheme has already been shown to improve snow and 2&amp;#8208;metre temperature, while in this study, the wider hydrological impact is evaluated and documented.&lt;/p&gt;&lt;p&gt;The analysis is done in the reanalysis context by comparing two ERA5-forced offline HTESSEL experiments. The runoff output of HTESSEL is coupled to the CaMa-Flood hydrodynamic model in order to derive river discharge. The analysis is done globally for the period between 1980-2018. The evaluation was carried out using over 1000 discharge observation time-series with varying catchment size. The hydrological response of the multi-layer snow scheme is generally positive, but in some areas the improvement is not clear and can even be negative with deteriorated signal in river discharge. Further investigation is needed to understand the complex hydrological impact of the new snow scheme, making sure it contributes to an improved description of all hydrological components of the Earth System.&lt;/p&gt;


2017 ◽  
Author(s):  
Michail Alvanos ◽  
Theodoros Christoudias

Abstract. This paper presents an application of GPU accelerators in Earth system modelling. We focus on atmospheric chemical kinetics, one of the most computationally intensive tasks in climate-chemistry model simulations. We developed a software package that automatically generates CUDA kernels to numerically integrate atmospheric chemical kinetics in the global climate model ECHAM/MESSy Atmospheric Chemistry (EMAC), used to study climate change and air quality scenarios. A source-to-source compiler outputs a CUDA compatible kernel, by parsing the FORTRAN code generated by the Kinetic Pre-Processor (KPP) general analysis tool. All Rosenbrock methods that are available in the KPP numerical library are supported. Performance evaluation, using Fermi and Pascal CUDA-enabled GPU accelerators shows achieved speedups of 4.5× and 22.4× respectively of the kernel execution time. A node-to-node real-world production performance comparison shows a 1.75× speed-up over the non-accelerated application using the KPP 3-stage Rosenbrock solver. We provide a detailed description of the code optimizations used to improve the performance including memory optimizations, control code simplification, and reduction of idle time. The accuracy and correctness of the accelerated implementation are evaluated by comparing to the CPU-only version of the application. The relative difference is found to be less than 0.00005 % when comparing the output of the accelerated kernel the CPU-only code, within the target level of relative accuracy (relative error tolerance) of 0.1 %. The approach followed, including the computational workload division and the developed GPU solver code can potentially be used as the basis for hardware acceleration of numerous geoscientific models that rely on KPP for atmospheric chemical kinetics applications.


Sign in / Sign up

Export Citation Format

Share Document