Search

Abstract The non‐local model of mixing based on internal wave breaking, IDEMIX, is implemented as an enhancement of a turbulent kinetic energy closure model in three non‐eddy resolving general circulation ocean models that differ in the discretization and choice of computational grids. In IDEMIX internal wave energy is generated by an energy flux resulting from near‐inertial waves induced by wind forcing at the surface, and at the bottom, by an energy flux that parameterizes the transfer of energy between baroclinic and barotropic tides. In all model simulations with IDEMIX, the mixing work is increased compared to the reference solutions without IDEMIX, reaching values in better agreement with finestructure observations. Furthermore, the horizontal structure of the mixing work is more realistic as a consequence of the heterogeneous forcing functions. All models with IDEMIX simulate deeper thermocline depths related to stronger shallow overturning cells in the Indo‐Pacific. In the North Atlantic, deeper mixed layers in simulations with IDEMIX are associated with an increased Atlantic overturning circulation and an increase of northward heat transports toward more realistic values. The response of the deep Indo‐Pacific overturning circulation and the weak bottom cell of the Atlantic to the inclusion of IDEMIX is incoherent between the models, suggesting that additional unidentified processes and numerical mixing may confound the analysis. Applying different tidal forcing functions leads to simulation differences that are small compared to differences between the different models or between simulations with IDEMIX and without IDEMIX.

The non‐local model of mixing based on internal wave breaking, IDEMIX, is implemented as an enhancement of a turbulent kinetic energy closure model in three non‐eddy resolving general circulation ocean models that differ in the discretization and choice of computational grids. In IDEMIX internal wave energy is generated by an energy flux resulting from near‐inertial waves induced by wind forcing at the surface, and at the bottom, by an energy flux that parameterizes the transfer of energy between baroclinic and barotropic tides. In all model simulations with IDEMIX, the mixing work is increased compared to the reference solutions without IDEMIX, reaching values in better agreement with finestructure observations. Furthermore, the horizontal structure of the mixing work is more realistic as a consequence of the heterogeneous forcing functions. All models with IDEMIX simulate deeper thermocline depths related to stronger shallow overturning cells in the Indo‐Pacific. In the North Atlantic, deeper mixed layers in simulations with IDEMIX are associated with an increased Atlantic overturning circulation and an increase of northward heat transports toward more realistic values. The response of the deep Indo‐Pacific overturning circulation and the weak bottom cell of the Atlantic to the inclusion of IDEMIX is incoherent between the models, suggesting that additional unidentified processes and numerical mixing may confound the analysis. Applying different tidal forcing functions leads to simulation differences that are small compared to differences between the different models or between simulations with IDEMIX and without IDEMIX. Plain Language Summary: Waves in the ocean interior play a fundamental role for ocean dynamics since they can carry energy over long distances and, once they break, lead to turbulent mixing. This turbulent mixing can cause dense water masses to rise from the deep ocean with a direct impact on large‐scale currents. The wave dynamics occur on spatial scales that cannot be resolved in global ocean or climate models. To account for these processes, we apply the new parameterization IDEMIX that describes internal wave generation, propagation, and mixing. Using three different ocean models with and without IDEMIX ensures that we can identify model‐specific effects of the parameterization and discriminate them from those independent of the model. We find that the simulated mixing patterns agree better with observations once IDEMIX is applied. Large‐scale currents and the vertical temperature distribution are substantially affected by the internal wave parameterization. Whether this leads to an improved agreement with observed currents and water mass properties depends on the specific model and on numerical effects. In most cases, simulations with IDEMIX are not very sensitive to details of how the internal wave model is driven by tidal energy input. Key Points: the IDEMIX closure for a consistent representation of internal wave‐induced mixing is evaluated in three state‐of‐the‐art ocean modelsonly in simulations with IDEMIX can the models reproduce the magnitude and spatial variability of the observed mixing workmost changes with IDEMIX can be attributed to stronger mixing, but some effects are confounded by other processes and numerical mixing [ABSTRACT FROM AUTHOR]

Abstract This work documents the ICON‐Earth System Model (ICON‐ESM V1.0), the first coupled model based on the ICON (ICOsahedral Non‐hydrostatic) framework with its unstructured, icosahedral grid concept. The ICON‐A atmosphere uses a nonhydrostatic dynamical core and the ocean model ICON‐O builds on the same ICON infrastructure, but applies the Boussinesq and hydrostatic approximation and includes a sea‐ice model. The ICON‐Land module provides a new framework for the modeling of land processes and the terrestrial carbon cycle. The oceanic carbon cycle and biogeochemistry are represented by the Hamburg Ocean Carbon Cycle module. We describe the tuning and spin‐up of a base‐line version at a resolution typical for models participating in the Coupled Model Intercomparison Project (CMIP). The performance of ICON‐ESM is assessed by means of a set of standard CMIP6 simulations. Achievements are well‐balanced top‐of‐atmosphere radiation, stable key climate quantities in the control simulation, and a good representation of the historical surface temperature evolution. The model has overall biases, which are comparable to those of other CMIP models, but ICON‐ESM performs less well than its predecessor, the Max Planck Institute Earth System Model. Problematic biases are diagnosed in ICON‐ESM in the vertical cloud distribution and the mean zonal wind field. In the ocean, sub‐surface temperature and salinity biases are of concern as is a too strong seasonal cycle of the sea‐ice cover in both hemispheres. ICON‐ESM V1.0 serves as a basis for further developments that will take advantage of ICON‐specific properties such as spatially varying resolution, and configurations at very high resolution.

Abstract This paper investigates the causation for the decline of the Atlantic Meridional Overturning Circulation (AMOC) from approximately 17 Sv to about 9 Sv, when the atmospheric resolution of the Max Planck Institute‐Earth System Model is enhanced from ∼1° to ∼0.5°. The results show that the slowdown of the AMOC is caused by the cessation of deep convection. In most modeling studies, this is thought to be controlled by buoyancy fluxes in the convective regions, for example, by surface freshwater flux that is introduced locally or via enormous input from glacier or iceberg melts. While we find that freshwater is still the key to the reduction of AMOC seen in the higher‐resolution run, the freshening of the North Atlantic does not need to be directly caused by local freshwater fluxes. Instead, it can be caused indirectly through winds via a reduced wind‐driven gyre circulation and salinity transport associated to this circulation, as seen in the higher‐resolution run.

Abstract The MPI‐ESM1.2 is the latest version of the Max Planck Institute Earth System Model and is the baseline for the Coupled Model Intercomparison Project Phase 6 and current seasonal and decadal climate predictions. This paper evaluates a coupled higher‐resolution version (MPI‐ESM1.2‐HR) in comparison with its lower‐resolved version (MPI‐ESM1.2‐LR). We focus on basic oceanic and atmospheric mean states and selected modes of variability, the El Niño/Southern Oscillation and the North Atlantic Oscillation. The increase in atmospheric resolution in MPI‐ESM1.2‐HR reduces the biases of upper‐level zonal wind and atmospheric jet stream position in the northern extratropics. This results in a decrease of the storm track bias over the northern North Atlantic, for both winter and summer season. The blocking frequency over the European region is improved in summer, and North Atlantic Oscillation and related storm track variations improve in winter. Stable Atlantic meridional overturning circulations are found with magnitudes of ~16 Sv for MPI‐ESM1.2‐HR and ~20 Sv for MPI‐ESM1.2‐LR at 26°N. A strong sea surface temperature bias of ~5°C along with a too zonal North Atlantic current is present in both versions. The sea surface temperature bias in the eastern tropical Atlantic is reduced by ~1°C due to higher‐resolved orography in MPI‐ESM‐HR, and the region of the cold‐tongue bias is reduced in the tropical Pacific. MPI‐ESM1.2‐HR has a well‐balanced radiation budget and its climate sensitivity is explicitly tuned to 3 K. Although the obtained reductions in long‐standing biases are modest, the improvements in atmospheric dynamics make this model well suited for prediction and impact studies.

During a development stage global climate models have their properties adjusted or tuned in various ways to best match the known state of the Earth's climate system. These desired properties are observables, such as the radiation balance at the top of the atmosphere, the global mean temperature, sea ice, clouds and wind fields. The tuning is typically performed by adjusting uncertain, or even non-observable, parameters related to processes not explicitly represented at the model grid resolution. The practice of climate model tuning has seen an increasing level of attention because key model properties, such as climate sensitivity, have been shown to depend on frequently used tuning parameters. Here we provide insights into how climate model tuning is practically done in the case of closing the radiation balance and adjusting the global mean temperature for the Max Planck Institute Earth System Model (MPI-ESM). We demonstrate that considerable ambiguity exists in the choice of parameters, and present and compare three alternatively tuned, yet plausible configurations of the climate model. The impacts of parameter tuning on climate sensitivity was less than anticipated.

Abstract: The MPI‐ESM1.2 is the latest version of the Max Planck Institute Earth System Model and is the baseline for the Coupled Model Intercomparison Project Phase 6 and current seasonal and decadal climate predictions. This paper evaluates a coupled higher‐resolution version (MPI‐ESM1.2‐HR) in comparison with its lower‐resolved version (MPI‐ESM1.2‐LR). We focus on basic oceanic and atmospheric mean states and selected modes of variability, the El Niño/Southern Oscillation and the North Atlantic Oscillation. The increase in atmospheric resolution in MPI‐ESM1.2‐HR reduces the biases of upper‐level zonal wind and atmospheric jet stream position in the northern extratropics. This results in a decrease of the storm track bias over the northern North Atlantic, for both winter and summer season. The blocking frequency over the European region is improved in summer, and North Atlantic Oscillation and related storm track variations improve in winter. Stable Atlantic meridional overturning circulations are found with magnitudes of ~16 Sv for MPI‐ESM1.2‐HR and ~20 Sv for MPI‐ESM1.2‐LR at 26°N. A strong sea surface temperature bias of ~5°C along with a too zonal North Atlantic current is present in both versions. The sea surface temperature bias in the eastern tropical Atlantic is reduced by ~1°C due to higher‐resolved orography in MPI‐ESM‐HR, and the region of the cold‐tongue bias is reduced in the tropical Pacific. MPI‐ESM1.2‐HR has a well‐balanced radiation budget and its climate sensitivity is explicitly tuned to 3 K. Although the obtained reductions in long‐standing biases are modest, the improvements in atmospheric dynamics make this model well suited for prediction and impact studies. [ABSTRACT FROM AUTHOR]

Abstract 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.

We describe the evolution of Arctic sea ice as modeled by the Max Planck Institute for Meteorology's Earth System Model (MPI‐ESM). The modeled spatial distribution and interannual variability of the sea‐ice cover agree well with satellite observations and are improved relative to the model's predecessor ECHAM5/MPIOM. An evaluation of modeled sea‐ice coverage based on sea‐ice area gives, however, conflicting results compared to an evaluation based on sea‐ice extent and is additionally hindered by uncertainties in the observational record. Simulated trends in sea‐ice coverage for the satellite period range from more strongly negative than observed to positive. The observed evolution of Arctic sea ice is incompatible with modeled internal variability and probably caused by external forcing. Simulated drift patterns agree well with observations, but simulated drift speed is generally too high. Simulated sea‐ice volume agrees well with volume estimates of the PIOMAS reanalysis for the past few years. However, a preceding Arctic wide decrease in sea‐ice volume starts much earlier in MPI‐ESM than in PIOMAS. Analyzing this behavior in MPI‐ESM's ocean model MPIOM, we find that the modeled volume trend depends crucially on the specific choice of atmospheric reanalysis forcing, which casts some doubt on the reliability of estimates of volume trends. In our CMIP5 scenario simulations, we find a substantial delay in sea‐ice response to increasing CO2 concentration; a seasonally ice‐free Arctic can result for a CO2 concentration of around 500 ppm. Simulated winter sea‐ice coverage drops rapidly to near ice‐free conditions once the mean Arctic winter temperature exceeds −5°C.

MPI‐ESM is a new version of the global Earth system model developed at the Max Planck Institute for Meteorology. This paper describes the ocean state and circulation as well as basic aspects of variability in simulations contributing to the fifth phase of the Coupled Model Intercomparison Project (CMIP5). The performance of the ocean/sea‐ice model MPIOM, coupled to a new version of the atmosphere model ECHAM6 and modules for land surface and ocean biogeochemistry, is assessed for two model versions with different grid resolution in the ocean. The low‐resolution configuration has a nominal resolution of 1.5°, whereas the higher resolution version features a quasiuniform, eddy‐permitting global resolution of 0.4°. The paper focuses on important oceanic features, such as surface temperature and salinity, water mass distribution, large‐scale circulation, and heat and freshwater transports. In general, these integral quantities are simulated well in comparison with observational estimates, and improvements in comparison with the predecessor system are documented; for example, for tropical variability and sea ice representation. Introducing an eddy‐permitting grid configuration in the ocean leads to improvements, in particular, in the representation of interior water mass properties in the Atlantic and in the representation of important ocean currents, such as the Agulhas and Equatorial current systems. In general, however, there are more similarities than differences between the two grid configurations, and several shortcomings, known from earlier versions of the coupled model, prevail.