scholarly journals ESM-Tools version 5.0: a modular infrastructure for stand-alone and coupled Earth system modelling (ESM)

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.

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.


2021 ◽  
Author(s):  
Dirk Barbi ◽  
Miguel Andrés-Martínez ◽  
Deniz Ural ◽  
Luisa Cristini ◽  
Paul Gierz ◽  
...  

<p>During the last two decades, modern societies have gradually understood the urge to tackle the climate change challenge, and consequently, a growing number of national and international initiatives have been launched with the aim of better understanding the Earth System. In this context, Earth System Modelling (ESM) has rapidly expanded, leading to a large number of research groups targeting the many components of the system at different scales and with different levels of interactions between components. This has led to the development of increasing number of models, couplings, versions tuned to address different scales or scenarios, and model-specific compilation and operating procedures. This operational complexity makes the implementation of multiple models excessively time consuming especially for less experienced modellers.</p><p>ESM-Tools is an open-source modular software written in Python, aimed to overcome many of the difficulties associated to the operation of ESMs. ESM-Tools allows for downloading, compiling and running a wide range of ESM models and coupled setups in the most important HPC facilities available in Germany. It currently supports multiple models for ocean, atmosphere, biochemistry, ice sheet, isostatic adjustment, hydrology, and land-surface, and six ocean-atmosphere and two ice-sheet-ocean-atmosphere coupled setups, through two couplers (included modularly through ESM-Interface). The tools are coded in Python while all the component and coupling information is contained in easy-to-read YAML files. The front-end user is required to provide only a short script written in YAML format, containing the experiment specific definitions. This user-friendly interface makes ESM-Tools a convenient software for training and educational purposes. Simultaneously, its modularity and the separation between the component-specific information and tool scripts facilitates the implementation and maintenance of new components, couplings and versions. ESM-Tools team of scientific programmers provides also user support, workshops and detailed documentation. The ESM-Tools were developed within the framework of the project Advance Earth System Model Capacity, supported by Helmholtz Association and has become one of the main pillars of the German infrastructure for Climate Modelling.</p>


2018 ◽  
Author(s):  
Trung Nguyen-Quang ◽  
Jan Polcher ◽  
Agnès Ducharne ◽  
Thomas Arsouze ◽  
Xudong Zhou ◽  
...  

Abstract. This study presents a revised river routing scheme (RRS) for the Organising Carbon and Hydrology in Dynamic Ecosystems (ORCHIDEE) land surface model. The revision is carried out to benefit from the high resolution topography provided the Hydrological data and maps based on SHuttle Elevation Derivatives at multiple Scales (HydroSHEDS), processed to a resolution of approximately 1 kilometer. The RRS scheme of the ORCHIDEE uses a unit-to-unit routing concept which allows to preserve as much of the hydrological information of the HydroSHEDS as the user requires. The evaluation focuses on 12 rivers of contrasted size and climate which contribute freshwater to the Mediterranean Sea. First, the numerical aspect of the new RRS is investigated, to identify the practical configuration offering the best trade-off between computational cost and simulation quality for ensuing validations. Second, the performance of the revised scheme is evaluated against observations at both monthly and daily timescales. The new RRS captures satisfactorily the seasonal variability of river discharges, although important biases come from the water budget simulated by the ORCHIDEE model. The results highlight that realistic streamflow simulations require accurate precipitation forcing data and a precise river catchment description over a wide range of scales, as permitted by the new RRS. Detailed analyses at the daily timescale show promising performances of this high resolution RRS for replicating river flow variation at various frequencies. Eventually, this RRS is well adapted for further developments in the ORCHIDEE land surface model to assess anthropogenic impacts on river processes (e.g. damming for irrigation operation).


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.


2010 ◽  
Vol 2 (2) ◽  
Author(s):  
Diandong Ren

AbstractBased on a 2-layer land surface model, a rather general variational data assimilation framework for estimating model state variables is developed. The method minimizes the error of surface soil temperature predictions subject to constraints imposed by the prediction model. Retrieval experiments for soil prognostic variables are performed and the results verified against model simulated data as well as real observations for the Oklahoma Atmospheric Surface layer Instrumentation System (OASIS). The optimization scheme is robust with respect to a wide range of initial guess errors in surface soil temperature (as large as 30 K) and deep soil moisture (within the range between wilting point and saturation). When assimilating OASIS data, the scheme can reduce the initial guess error by more than 90%, while for Observing Simulation System Experiments (OSSEs), the initial guess error is usually reduced by over four orders of magnitude.Using synthetic data, the robustness of the retrieval scheme as related to information content of the data and the physical meaning of the adjoint variables and their use in sensitivity studies are investigated. Through sensitivity analysis, it is confirmed that the vegetation coverage and growth condition determine whether or not the optimally estimated initial soil moisture condition leads to an optimal estimation of the surface fluxes. This reconciles two recent studies.With the real data experiments, it is shown that observations during the daytime period are the most effective for the retrieval. Longer assimilation windows result in more accurate initial condition retrieval, underlining the importance of information quantity, especially for schemes assimilating noisy observations.


2015 ◽  
Vol 16 (4) ◽  
pp. 1502-1520 ◽  
Author(s):  
Elizabeth A. Clark ◽  
Justin Sheffield ◽  
Michelle T. H. van Vliet ◽  
Bart Nijssen ◽  
Dennis P. Lettenmaier

Abstract A common term in the continental and oceanic components of the global water cycle is freshwater discharge to the oceans. Many estimates of the annual average global discharge have been made over the past 100 yr with a surprisingly wide range. As more observations have become available and continental-scale land surface model simulations of runoff have improved, these past estimates are cast in a somewhat different light. In this paper, a combination of observations from 839 river gauging stations near the outlets of large river basins is used in combination with simulated runoff fields from two implementations of the Variable Infiltration Capacity land surface model to estimate continental runoff into the world’s oceans from 1950 to 2008. The gauges used account for ~58% of continental areas draining to the ocean worldwide, excluding Greenland and Antarctica. This study estimates that flows to the world’s oceans globally are 44 200 (±2660) km3 yr−1 (9% from Africa, 37% from Eurasia, 30% from South America, 16% from North America, and 8% from Australia–Oceania). These estimates are generally higher than previous estimates, with the largest differences in South America and Australia–Oceania. Given that roughly 42% of ocean-draining continental areas are ungauged, it is not surprising that estimates are sensitive to the land surface and hydrologic model (LSM) used, even with a correction applied to adjust for model bias. The results show that more and better in situ streamflow measurements would be most useful in reducing uncertainties, in particular in the southern tip of South America, the islands of Oceania, and central Africa.


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

<p>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‐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‐metre temperature, while in this study, the wider hydrological impact is evaluated and documented.</p><p>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.</p>


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.


2018 ◽  
Vol 11 (12) ◽  
pp. 4965-4985 ◽  
Author(s):  
Trung Nguyen-Quang ◽  
Jan Polcher ◽  
Agnès Ducharne ◽  
Thomas Arsouze ◽  
Xudong Zhou ◽  
...  

Abstract. The river routing scheme (RRS) in the Organising Carbon and Hydrology in Dynamic Ecosystems (ORCHIDEE) land surface model is a valuable tool for closing the water cycle in a coupled environment and for validating the model performance. This study presents a revision of the RRS of the ORCHIDEE model that aims to benefit from the high-resolution topography provided by the Hydrological data and maps based on SHuttle Elevation Derivatives at multiple Scales (HydroSHEDS), which is processed to a resolution of approximately 1 km. Adapting a new algorithm to construct river networks, the new RRS in ORCHIDEE allows for the preservation of as much of the hydrological information from HydroSHEDS as the user requires. The evaluation focuses on 12 rivers of contrasting size and climate which contribute freshwater to the Mediterranean Sea. First, the numerical aspect of the new RRS is investigated, in order to identify the practical configuration offering the best trade-off between computational cost and simulation quality for ensuing validations. Second, the performance of the new scheme is evaluated against observations at both monthly and daily timescales. The new RRS satisfactorily captures the seasonal variability of river discharge, although important biases stem from the water budget simulated by the ORCHIDEE model. The results highlight that realistic streamflow simulations require accurate precipitation forcing data and a precise river catchment description over a wide range of scales, as permitted by the new RRS. Detailed analyses at the daily timescale show the promising performance of this high-resolution RRS with respect to replicating river flow variation at various frequencies. Furthermore, this RRS may also eventually be well adapted for further developments in the ORCHIDEE land surface model to assess anthropogenic impacts on river processes (e.g. damming for irrigation operation).


2019 ◽  
Author(s):  
Ting Sun ◽  
Sue Grimmond

Abstract. Accurate and agile modelling of the climate of cities is essential for urban climate services. The Surface Urban Energy and Water balance Scheme (SUEWS) is a state-of-the-art, widely used, urban land surface model (ULSM) which simulates urban-atmospheric interactions by quantifying the energy, water and mass fluxes. Using SUEWS as the computation kernel, SuPy (SUEWS in Python), stands on the Python-based data stack to streamline the pre-processing, computation and post-processing that are involved in the common modelling-centred urban climate studies. This paper documents the development of SuPy, which includes the SUEWS interface modification, F2PY (Fortran to Python) configuration and Python frontend implementation. In addition, the deployment of SuPy via PyPI (Python Package Index) is introduced along with the automated workflow for cross-platform compilation. This makes SuPy available for all mainstream operating systems (Windows, Linux, and macOS). Furthermore, three online tutorials in Jupyter notebooks are provided to users of different levels to become familiar with SuPy urban climate modelling. The SuPy package represents a significant enhancement that supports existing and new model applications, reproducibility, and enhanced functionality.


Sign in / Sign up

Export Citation Format

Share Document