Explore the vast range of titles available.

The Max Planck Institute Grand Ensemble (MPI‐GE) is the largest ensemble of a single comprehensive climate model currently available, with 100 members for the historical simulations (1850–2005) and four forcing scenarios. It is currently the only large ensemble available that includes scenario representative concentration pathway (RCP) 2.6 and a 1% CO2 scenario. These advantages make MPI‐GE a powerful tool. We present an overview of MPI‐GE, its components, and detail the experiments completed. We demonstrate how to separate the forced response from internal variability in a large ensemble. This separation allows the quantification of both the forced signal under climate change and the internal variability to unprecedented precision. We then demonstrate multiple ways to evaluate MPI‐GE and put observations in the context of a large ensemble, including a novel approach for comparing model internal variability with estimated observed variability. Finally, we present four novel analyses, which can only be completed using a large ensemble. First, we address whether temperature and precipitation have a pathway dependence using the forcing scenarios. Second, the forced signal of the highly noisy atmospheric circulation is computed, and different drivers are identified to be important for the North Pacific and North Atlantic regions. Third, we use the ensemble dimension to investigate the time dependency of Atlantic Meridional Overturning Circulation variability changes under global warming. Last, sea level pressure is used as an example to demonstrate how MPI‐GE can be utilized to estimate the ensemble size needed for a given scientific problem and provide insights for future ensemble projects. Key Points The 100‐member MPI‐GE is currently the largest publicly available ensemble of a comprehensive climate model MPI‐GE currently has the most forcing scenarios of all large ensemble projects: RCP2.6, RCP4.5, RCP8.5, and 1% CO2 The power of MPI‐GE is to estimate the forced response and internal variability, including changing variability, to unprecedented precision

A review of the experimental protocol and motivation for DYAMOND, the first intercomparison project of global storm-resolving models, is presented. Nine models submitted simulation output for a 40-day (1 August–10 September 2016) intercomparison period. Eight of these employed a tiling of the sphere that was uniformly less than 5 km. By resolving the transient dynamics of convective storms in the tropics, global storm-resolving models remove the need to parameterize tropical deep convection, providing a fundamentally more sound representation of the climate system and a more natural link to commensurately high-resolution data from satellite-borne sensors. The models and some basic characteristics of their output are described in more detail, as is the availability and planned use of this output for future scientific study. Tropically and zonally averaged energy budgets, precipitable water distributions, and precipitation from the model ensemble are evaluated, as is their representation of tropical cyclones and the predictability of column water vapor, the latter being important for tropical weather.

ECHAM6, the sixth generation of the atmospheric general circulation model ECHAM, is described. Major changes with respect to its predecessor affect the representation of shortwave radiative transfer, the height of the model top. Minor changes have been made to model tuning and convective triggering. Several model configurations, differing in horizontal and vertical resolution, are compared. As horizontal resolution is increased beyond T63, the simulated climate improves but changes are incremental; major biases appear to be limited by the parameterization of small‐scale physical processes, such as clouds and convection. Higher vertical resolution in the middle atmosphere leads to a systematic reduction in temperature biases in the upper troposphere, and a better representation of the middle atmosphere and its modes of variability. ECHAM6 represents the present climate as well as, or better than, its predecessor. The most marked improvements are evident in the circulation of the extratropics. ECHAM6 continues to have a good representation of tropical variability. A number of biases, however, remain. These include a poor representation of low‐level clouds, systematic shifts in major precipitation features, biases in the partitioning of precipitation between land and sea (particularly in the tropics), and midlatitude jets that appear to be insufficiently poleward. The response of ECHAM6 to increasing concentrations of greenhouse gases is similar to that of ECHAM5. The equilibrium climate sensitivity of the mixed‐resolution (T63L95) configuration is between 2.9 and 3.4 K and is somewhat larger for the 47 level model. Cloud feedbacks and adjustments contribute positively to warming from increasing greenhouse gases. Key Points To describe ECHAM6, as it was configured for participation in CMIP5 To describe the climate of ECHAM6 To describe the climate sensitivity of ECHAM6

We discuss scientific features and computational performance of kilometre-scale global weather and climate simulations, considering the Icosahedral Non-hydrostatic (ICON) model and the Integrated Forecast System (IFS). Scalability measurements and a performance modelling approach are used to derive performance estimates for these models on upcoming exascale supercomputers. This is complemented by preliminary analyses of the model data that illustrate the importance of high-resolution models to gain improvements in the accuracy of convective processes, a better understanding of physics dynamics interactions and poorly resolved or parametrized processes, such as gravity waves, convection and boundary layer. This article is part of the theme issue ‘Multiscale modelling, simulation and computing: from the desktop to the exascale’.

The key to the improvement of the quality of ensemble forecasts assessing the inherent flow uncertainties is the choice of the initial ensemble perturbations. To generate such perturbations, the breeding of growing modes approach has been used for the past two decades. Here, the fastest-growing error modes of the initial model state are estimated. However, the resulting bred vectors (BVs) mainly point in the phase space direction of the leading Lyapunov vector and therefore favor one direction of growing errors. To overcome this characteristic and obtain growing modes pointing to Lyapunov vectors different from the leading one, an orthogonalization implemented as a singular value decomposition based on the similarity between the BVs is applied. This transformation is similar to that used in the ensemble transform technique currently in operational use at NCEP but with certain differences in the metric used and in the implementation. In this study, results of this approach using BVs generated in the Ensemble Forecasting System (EFS) based on the global numerical weather prediction model GME of the German Meteorological Service are presented. The gain in forecast performance achieved with the orthogonalized BV initialization is shown by using different probabilistic forecast scores evaluating ensemble reliability, variance, and resolution. For a 3-month period in summer 2007, the results are compared to forecasts generated with simple BV initializations of the same ensemble prediction system as well as operational ensemble forecasts from ECMWF and NCEP. The orthogonalization vastly improves the GME–EFS scores and makes them competitive with the two other centers.

A new release of the Max Planck Institute for Meteorology Earth System Model version 1.2 (MPI‐ESM1.2) is presented. The development focused on correcting errors in and improving the physical processes representation, as well as improving the computational performance, versatility, and overall user friendliness. In addition to new radiation and aerosol parameterizations of the atmosphere, several relatively large, but partly compensating, coding errors in the model's cloud, convection, and turbulence parameterizations were corrected. The representation of land processes was refined by introducing a multilayer soil hydrology scheme, extending the land biogeochemistry to include the nitrogen cycle, replacing the soil and litter decomposition model and improving the representation of wildfires. The ocean biogeochemistry now represents cyanobacteria prognostically in order to capture the response of nitrogen fixation to changing climate conditions and further includes improved detritus settling and numerous other refinements. As something new, in addition to limiting drift and minimizing certain biases, the instrumental record warming was explicitly taken into account during the tuning process. To this end, a very high climate sensitivity of around 7 K caused by low‐level clouds in the tropics as found in an intermediate model version was addressed, as it was not deemed possible to match observed warming otherwise. As a result, the model has a climate sensitivity to a doubling of CO2 over preindustrial conditions of 2.77 K, maintaining the previously identified highly nonlinear global mean response to increasing CO2 forcing, which nonetheless can be represented by a simple two‐layer model. Key Points An updated version of the Max Planck Institute for Meteorology Earth System Model (MPI‐ESM1.2) is presented The model includes both code corrections and parameterization improvements Despite this, the model maintains an equilibrium climate sensitivity, which rises with warming

The double Inter-Tropical Convergence Zone (double-ITCZ) bias has been a persistent tropical precipitation bias over many climate model generations. This motivates investigations of whether increasing resolution and discarding parameterizations improves the representation of the large-scale atmospheric circulation and ITCZ. In this work, we study the double-ITCZ bias in an ICON XPP resolution hierarchy spanning from parameterized to explicitly described deep convection within a consistent framework. We demonstrate that the double ITCZ persists across horizontal resolutions from 160 to 5 km in specified sea-surface temperature simulations, independent of deep-convective and non-orographic gravity wave parameterizations. Changes in the treatment of near-surface wind speed within the turbulence parameterization can reduce the bias. However, we highlight that a key driver of the double-ITCZ bias seems to lie in the insufficient moisture transport from the subtropics to the inner tropics. The resulting dry bias in tropical near-surface moisture reduces deep convection over the Warm Pool, leading to a weakened Walker circulation. These biases ultimately culminate in the double-ITCZ feature. Increasing the near-surface wind speed limiter improves tropical near-surface moisture but exacerbates the bias in the moisture source, increasing the inner-tropical contribution at the expense of the subtropics. This degrades the representation of the global circulation, energy balance, and teleconnections. Additionally, we show that parameter adjustments at low resolution are informative of the response to the same parameter adjustments at high resolution. Our findings showcase the benefits of models supporting a range of resolutions and underline the importance of continuing the development of non-discardable parameterizations.

State-of-the-art Earth system models typically employ grid spacings of O(100 km), which is too coarse to explicitly resolve main drivers of the flow of energy and matter across the Earth system. In this paper, we present the new ICON-Sapphire model configuration, which targets a representation of the components of the Earth system and their interactions with a grid spacing of 10 km and finer. Through the use of selected simulation examples, we demonstrate that ICON-Sapphire can (i) be run coupled globally on seasonal timescales with a grid spacing of 5 km, on monthly timescales with a grid spacing of 2.5 km, and on daily timescales with a grid spacing of 1.25 km; (ii) resolve large eddies in the atmosphere using hectometer grid spacings on limited-area domains in atmosphere-only simulations; (iii) resolve submesoscale ocean eddies by using a global uniform grid of 1.25 km or a telescoping grid with the finest grid spacing at 530 m, the latter coupled to a uniform atmosphere; and (iv) simulate biogeochemistry in an ocean-only simulation integrated for 4 years at 10 km. Comparison of basic features of the climate system to observations reveals no obvious pitfalls, even though some observed aspects remain difficult to capture. The throughput of the coupled 5 km global simulation is 126 simulated days per day employing 21 % of the latest machine of the German Climate Computing Center. Extrapolating from these results, multi-decadal global simulations including interactive carbon are now possible, and short global simulations resolving large eddies in the atmosphere and submesoscale eddies in the ocean are within reach.

Classical numerical models for the global atmosphere, as used for numerical weather forecasting or climate research, have been developed for conventional central processing unit (CPU) architectures. This hinders the employment of such models on current top-performing supercomputers, which achieve their computing power with hybrid architectures, mostly using graphics processing units (GPUs). Thus also scientific applications of such models are restricted to the lesser computer power of CPUs. Here we present the development of a GPU-enabled version of the ICON atmosphere model (ICON-A), motivated by a research project on the quasi-biennial oscillation (QBO), a global-scale wind oscillation in the equatorial stratosphere that depends on a broad spectrum of atmospheric waves, which originates from tropical deep convection. Resolving the relevant scales, from a few kilometers to the size of the globe, is a formidable computational problem, which can only be realized now on top-performing supercomputers. This motivated porting ICON-A, in the specific configuration needed for the research project, in a first step to the GPU architecture of the Piz Daint computer at the Swiss National Supercomputing Centre and in a second step to the JUWELS Booster computer at the Forschungszentrum Jülich. On Piz Daint, the ported code achieves a single-node GPU vs. CPU speedup factor of 6.4 and allows for global experiments at a horizontal resolution of 5 km on 1024 computing nodes with 1 GPU per node with a turnover of 48 simulated days per day. On JUWELS Booster, the more modern hardware in combination with an upgraded code base allows for simulations at the same resolution on 128 computing nodes with 4 GPUs per node and a turnover of 133 simulated days per day. Additionally, the code still remains functional on CPUs, as is demonstrated by additional experiments on the Levante compute system at the German Climate Computing Center. While the application shows good weak scaling over the tested 16-fold increase in grid size and node count, making also higher resolved global simulations possible, the strong scaling on GPUs is relatively poor, which limits the options to increase turnover with more nodes. Initial experiments demonstrate that the ICON-A model can simulate downward-propagating QBO jets, which are driven by wave–mean flow interaction.

The new Max‐Planck‐Institute Earth System Model (MPI‐ESM) is used in the Coupled Model Intercomparison Project phase 5 (CMIP5) in a series of climate change experiments for either idealized CO2‐only forcing or forcings based on observations and the Representative Concentration Pathway (RCP) scenarios. The paper gives an overview of the model configurations, experiments related forcings, and initialization procedures and presents results for the simulated changes in climate and carbon cycle. It is found that the climate feedback depends on the global warming and possibly the forcing history. The global warming from climatological 1850 conditions to 2080–2100 ranges from 1.5°C under the RCP2.6 scenario to 4.4°C under the RCP8.5 scenario. Over this range, the patterns of temperature and precipitation change are nearly independent of the global warming. The model shows a tendency to reduce the ocean heat uptake efficiency toward a warmer climate, and hence acceleration in warming in the later years. The precipitation sensitivity can be as high as 2.5% K−1 if the CO2 concentration is constant, or as small as 1.6% K−1, if the CO2 concentration is increasing. The oceanic uptake of anthropogenic carbon increases over time in all scenarios, being smallest in the experiment forced by RCP2.6 and largest in that for RCP8.5. The land also serves as a net carbon sink in all scenarios, predominantly in boreal regions. The strong tropical carbon sources found in the RCP2.6 and RCP8.5 experiments are almost absent in the RCP4.5 experiment, which can be explained by reforestation in the RCP4.5 scenario. Key Points The climate feedback in MPI‐ESM is non‐linear and depends on the forcing history Ocean heat uptake is reduced in a warmer climate. Patterns of temperature and precipitation changes are robust for RCP26/45/85.