ALBERT

Description: The Earth System Model Evaluation Tool (ESMValTool) is a community diagnostics and performance metrics tool designed to improve comprehensive and routine evaluation of Earth system models (ESMs) participating in the Coupled Model Intercomparison Project (CMIP). It has undergone rapid development since the first release in 2016 and is now a well-tested tool that provides end-to-end provenance tracking to ensure reproducibility. It consists of (1) an easy-to-install, well-documented Python package providing the core functionalities (ESMValCore) that performs common preprocessing operations and (2) a diagnostic part that includes tailored diagnostics and performance metrics for specific scientific applications. Here we describe large-scale diagnostics of the second major release of the tool that supports the evaluation of ESMs participating in CMIP Phase 6 (CMIP6). ESMValTool v2.0 includes a large collection of diagnostics and performance metrics for atmospheric, oceanic, and terrestrial variables for the mean state, trends, and variability. ESMValTool v2.0 also successfully reproduces figures from the evaluation and projections chapters of the Intergovernmental Panel on Climate Change (IPCC) Fifth Assessment Report (AR5) and incorporates updates from targeted analysis packages, such as the NCAR Climate Variability Diagnostics Package for the evaluation of modes of variability, the Thermodynamic Diagnostic Tool (TheDiaTo) to evaluate the energetics of the climate system, as well as parts of AutoAssess that contains a mix of top–down performance metrics. The tool has been fully integrated into the Earth System Grid Federation (ESGF) infrastructure at the Deutsches Klimarechenzentrum (DKRZ) to provide evaluation results from CMIP6 model simulations shortly after the output is published to the CMIP archive. A result browser has been implemented that enables advanced monitoring of the evaluation results by a broad user community at much faster timescales than what was possible in CMIP5.

Description: Processes driving the production, transformation and transport of methane (CH4) in wetland ecosystems are highly complex. We present a simple calculation algorithm to separate open-water CH4 fluxes measured with automatic chambers into diffusion- and ebullition-derived components. This helps to reveal underlying dynamics, to identify potential environmental drivers and, thus, to calculate reliable CH4 emission estimates. The flux separation is based on identification of ebullition-related sudden concentration changes during single measurements. Therefore, a variable ebullition filter is applied, using the lower and upper quartile and the interquartile range (IQR). Automation of data processing is achieved by using an established R script, adjusted for the purpose of CH4 flux calculation. The algorithm was validated by performing a laboratory experiment and tested using flux measurement data (July to September 2013) from a former fen grassland site, which converted into a shallow lake as a result of rewetting. Ebullition and diffusion contributed equally (46 and 55 %) to total CH4 emissions, which is comparable to ratios given in the literature. Moreover, the separation algorithm revealed a concealed shift in the diurnal trend of diffusive fluxes throughout the measurement period. The water temperature gradient was identified as one of the major drivers of diffusive CH4 emissions, whereas no significant driver was found in the case of erratic CH4 ebullition events.

Description: The coupling of models is a commonly used approach when addressing the complex interactions between different components of earth systems. We demonstrate that this approach can result in a reduction of errors in wave forecasting, especially in dynamically complicated coastal ocean areas, such as the southern part of the North Sea – the German Bight. Here, we study the effects of coupling of an atmospheric model (COSMO) and a wind wave model (WAM), which is enabled by implementing wave-induced drag in the atmospheric model. The numerical simulations use a regional North Sea coupled wave–atmosphere model as well as a nested-grid high-resolution German Bight wave model. Using one atmospheric and two wind wave models simultaneously allows for study of the individual and combined effects of two-way coupling and grid resolution. This approach proved to be particularly important under severe storm conditions as the German Bight is a very shallow and dynamically complex coastal area exposed to storm floods. The two-way coupling leads to a reduction of both surface wind speeds and simulated wave heights. In this study, the sensitivity of atmospheric parameters, such as wind speed and atmospheric pressure, to the wave-induced drag, in particular under storm conditions, and the impact of two-way coupling on the wave model performance, is quantified. Comparisons between data from in situ and satellite altimeter observations indicate that two-way coupling improves the simulation of wind and wave parameters of the model and justify its implementation for both operational and climate simulations.

Description: The performance of improved versions of the four earth system models (ESMs) CNRM, EC-Earth, HadGEM, and MPI-ESM is assessed in comparison to their predecessor versions used in Phase 5 of the Coupled Model Intercomparison Project. The earth System Model Evaluation Tool (ESMValTool) is applied to evaluate selected climate phenomena in the models against observations. This is the first systematic application of the ESMValTool to assess and document the progress made during an extensive model development and improvement project. This study focuses on the South Asian (SAM) and West African (WAM) monsoons, the coupled equatorial climate, and Southern Ocean clouds and radiation, which are known to exhibit systematic biases in present-day ESMs. The analysis shows that the tropical precipitation in three out of four models is significantly improved. Two of three updated coupled models show an improved representation of tropical sea surface temperatures with one coupled model not exhibiting a double Inter-Tropical Convergence Zone. Simulated cloud amounts and cloud-radiation interactions are improved over the Southern Ocean. Improvements are also seen in the simulation of the SAM and WAM, although systematic biases remain in regional details and the timing of monsoon rainfall. Analysis of simulations with EC-Earth at different horizontal resolutions from T159 up to T1279 shows that the synoptic-scale variability of precipitation over the SAM and WAM regions improves with higher model resolution. The results suggest the reasonably good agreement of modeled and observed mean WAM and SAM rainfall in lower resolution models may be a result of unrealistic intensity distributions.

Description: Drainage and cultivation of fen peatlands creates complex small-scale mosaics of soils with extremely variable soil organic carbon (SOC) stocks and groundwater-level (GWL). To date, it remains unclear if such sites are sources or sinks for greenhouse gases like CO2 and CH4, especially if used for cropland. As individual control factors like GWL fail to account for this complexity, holistic approaches combining gas fluxes with the underlying processes are required to understand the carbon (C) gas exchange of drained fens. It can be assumed that the stocks of SOC and N located above the variable GWL – defined as dynamic C and N stocks – play a key role in the regulation of plant- and microbially mediated C gas fluxes of these soils. To test this assumption, the present study analysed the C gas exchange (gross primary production – GPP, ecosystem respiration – Reco, net ecosystem exchange – NEE, CH4) of maize using manual chambers for four years. The study sites were located near Paulinenaue, Germany. Here we selected three soils, which represent the full gradient in pedogenesis, GWL and SOC stocks (0–1 m) of the fen peatland: (a) Haplic Arenosol (AR; 8 kg C m−2); (b) Mollic Gleysol (GL; 38 kg C m−2); and (c) Hemic Histosol (HS; 87 kg C m−2). Daily GWL data was used to calculate dynamic SOC (SOCdyn) and N (Ndyn) stocks. Average annual NEE differed considerably among sites, ranging from 47 ± 30 g C m−2 a−1 at AR to −305 ± 123 g C m−2 a−1 at GL and −127 ± 212 g C m−2 a−1 at HS. While static SOC and N stocks showed no significant effect on C fluxes, SOCdyn and Ndyn and their interaction with GWL strongly influenced the C gas exchange, particularly NEE and the GPP:Reco ratio. Moreover, based on nonlinear regression analysis, 86% of NEE variability was explained by GWL and SOCdyn. The observed high relevance of dynamic SOC and N stocks in the aerobic zone for plant and soil gas exchange likely originates from the effects of GWL-dependent N availability on C formation and transformation processes in the plant-soil system, which promote CO2 input via GPP more than CO2 emission via Reco. The process-oriented approach of dynamic C and N stocks is a promising, potentially generalizable method for system-oriented investigations of the C gas exchange of groundwater-influenced soils and could be expanded to other nutrients and soil characteristics. However, in order to assess the climate impact of arable sites on drained peatlands, it is always necessary to consider the entire range of groundwater-influenced mineral and organic soils and their respective areal extent within the soil landscape.

Description: We examined parameter optimisation in the JSBACH (Kaminski et al., 2013; Knorr and Kattge, 2005; Reick et al., 2013) ecosystem model, applied to two boreal forest sites (Hyytiälä and Sodankylä) in Finland. We identified and tested key parameters in soil hydrology and forest water and carbon-exchange-related formulations, and optimised them using the adaptive Metropolis (AM) algorithm for Hyytiälä with a 5-year calibration period (2000–2004) followed by a 4-year validation period (2005–2008). Sodankylä acted as an independent validation site, where optimisations were not made. The tuning provided estimates for full distribution of possible parameters, along with information about correlation, sensitivity and identifiability. Some parameters were correlated with each other due to a phenomenological connection between carbon uptake and water stress or other connections due to the set-up of the model formulations. The latter holds especially for vegetation phenology parameters. The least identifiable parameters include phenology parameters, parameters connecting relative humidity and soil dryness, and the field capacity of the skin reservoir. These soil parameters were masked by the large contribution from vegetation transpiration. In addition to leaf area index and the maximum carboxylation rate, the most effective parameters adjusting the gross primary production (GPP) and evapotranspiration (ET) fluxes in seasonal tuning were related to soil wilting point, drainage and moisture stress imposed on vegetation. For daily and half-hourly tunings the most important parameters were the ratio of leaf internal CO2 concentration to external CO2 and the parameter connecting relative humidity and soil dryness. Effectively the seasonal tuning transferred water from soil moisture into ET, and daily and half-hourly tunings reversed this process. The seasonal tuning improved the month-to-month development of GPP and ET, and produced the most stable estimates of water use efficiency. When compared to the seasonal tuning, the daily tuning is worse on the seasonal scale. However, daily parametrisation reproduced the observations for average diurnal cycle best, except for the GPP for Sodankylä validation period, where half-hourly tuned parameters were better. In general, the daily tuning provided the largest reduction in model–data mismatch. The models response to drought was unaffected by our parametrisations and further studies are needed into enhancing the dry response in JSBACH.

Description: Soil moisture is an essential climate variable (ECV) of major importance for land–atmosphere interactions and global hydrology. An appropriate representation of soil moisture dynamics in global climate models is therefore important. Recently, a first multidecadal, observation-based soil moisture dataset has become available that provides information on soil moisture dynamics from satellite observations (ECVSM, essential climate variable soil moisture). The present study investigates the potential and limitations of this new dataset for several applications in climate model evaluation. We compare soil moisture data from satellite observations, reanalysis and simulations from a state-of-the-art land surface model and analyze relationships between soil moisture and precipitation anomalies in the different dataset. Other potential applications like model parameter optimization or model initialization are not investigated in the present study. In a detailed regional study, we show that ECVSM is capable to capture well the interannual and intraannual soil moisture and precipitation dynamics in the Sahelian region. Current deficits of the new dataset are critically discussed and summarized at the end of the paper to provide guidance for an appropriate usage of the ECVSM dataset for climate studies.

Description: Due to biases in the output of climate models, a bias correction is often needed to make the output suitable for use in hydrological simulations. In most cases only the temperature and precipitation values are bias corrected. However, often there are also biases in other variables such as radiation, humidity and wind speed. In this study we tested to what extent it is also needed to bias correct these variables. Responses to radiation, humidity and wind estimates from two climate models for four large-scale hydrological models are analysed. For the period 1971–2000 these hydrological simulations are compared to simulations using meteorological data based on observations and reanalysis; i.e. the baseline simulation. In both forcing datasets originating from climate models precipitation and temperature are bias corrected to the baseline forcing dataset. Hence, it is only effects of radiation, humidity and wind estimates that are tested here. The direct use of climate model outputs result in substantial different evapotranspiration and runoff estimates, when compared to the baseline simulations. A simple bias correction method is implemented and tested by rerunning the hydrological models using bias corrected radiation, humidity and wind values. The results indicate that bias correction can successfully be used to match the baseline simulations. Finally, historical (1971–2000) and future (2071–2100) model simulations resulting from using bias corrected forcings are compared to the results using non-bias corrected forcings. The relative changes in simulated evapotranspiration and runoff are relatively similar for the bias corrected and non bias corrected hydrological projections, although the absolute evapotranspiration and runoff numbers are often very different. The simulated relative and absolute differences when using bias corrected and non bias corrected climate model radiation, humidity and wind values are, however, smaller than literature reported differences resulting from using bias corrected and non bias corrected climate model precipitation and temperature values.

Description: It is well known that output from climate models cannot be used to force hydrological simulations without some form of preprocessing to remove the existing biases. In principle, statistical bias correction methodologies act on model output so the statistical properties of the corrected data match those of the observations. However, the improvements to the statistical properties of the data are limited to the specific timescale of the fluctuations that are considered. For example, a statistical bias correction methodology for mean daily temperature values might be detrimental to monthly statistics. Also, in applying bias corrections derived from present day to scenario simulations, an assumption is made on the stationarity of the bias over the largest timescales. First, we point out several conditions that have to be fulfilled by model data to make the application of a statistical bias correction meaningful. We then examine the effects of mixing fluctuations on different timescales and suggest an alternative statistical methodology, referred to here as a cascade bias correction method, that eliminates, or greatly reduces, the negative effects.

Description: Many studies have investigated potential climate change impacts on regional hydrology; less attention has been given to the components of uncertainty that affect these scenarios. This study quantifies uncertainties resulting from (i) General Circulation Models (GCMs), (ii) Regional Climate Models (RCMs), (iii) bias-correction of RCMs, and (iv) hydrological model parameterization using a multi-model framework. This consists of three GCMs, three RCMs, three bias-correction techniques, and sets of hydrological model parameters. The study is performed for the Lech watershed (~ 1000 km2), located in the Northern Limestone Alps, Austria. Bias-corrected climate data are used to drive the hydrological model HQsim to simulate runoff under present (1971–2000) and future (2070–2099) climate conditions. Hydrological model parameter uncertainty is assessed by Monte Carlo sampling. The model chain is found to perform well under present climate conditions. However, hydrological projections are associated with high uncertainty, mainly due to the choice of GCM and RCM. Uncertainty due to bias-correction is found to have greatest influence on projections of extreme river flows, and the choice of method(s) is an important consideration in snowmelt systems. Overall, hydrological model parameterization is least important. The study also demonstrates how an improved understanding of the physical processes governing future river flows can help focus attention on the scientifically tractable elements of the uncertainty.