Search

Abstract We evaluate the new icosahedral nonhydrostatic atmospheric (ICON‐A) general circulation model of the Max Planck Institute for Meteorology that is flexible to be run at grid spacings from a few tens of meters to hundreds of kilometers. A simulation with ICON‐A at a low resolution (160 km) is compared to a not‐tuned fourfold higher‐resolution simulation (40 km). Simulations using the last release of the ECHAM climate model (ECHAM6.3) are also presented at two different resolutions. The ICON‐A simulations provide a compelling representation of the climate and its variability. The climate of the low‐resolution ICON‐A is even slightly better than that of ECHAM6.3. Improvements are obtained in aspects that are sensitive to the representation of orography, including the representation of cloud fields over eastern‐boundary currents, the latitudinal distribution of cloud top heights, and the spatial distribution of convection over the Indian Ocean and the Maritime Continent. Precipitation over land is enhanced, in particular at high‐resolution ICON‐A. The response of precipitation to El Niño sea surface temperature variability is close to observations, particularly over the eastern Indian Ocean. Some parameterization changes lead to improvements, for example, with respect to rain intensities and the representation of equatorial waves, but also imply a warmer troposphere, which we suggest leads to an unrealistic poleward mass shift. Many biases familiar to ECHAM6.3 are also evident in ICON‐A, namely, a too zonal SPCZ, an inadequate representation of north hemispheric blocking, and a relatively poor representation of tropical intraseasonal variability.

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.

A long-standing task in climate research has been to distinguish between anthropogenic climate change and natural climate variability. A prerequisite for fulfilling this task is the understanding of the relative roles of external drivers and internal variability of climate and the carbon cycle. Here, we present the first ensemble simulations over the last 1200 years with a comprehensive Earth system model including a fully interactive carbon cycle. Applying up-to-date reconstructions of external forcing including the recent low-amplitude estimates of solar variations, the ensemble simulations reproduce temperature evolutions consistent with the range of reconstructions. The 20th-century warming trend stands out against all pre-industrial trends within the ensemble. Volcanic eruptions are necessary to explain variations in pre-industrial climate such as the Little Ice Age; yet only the strongest, repeated eruptions lead to cooling trends that differ significantly from the internal variability across all ensemble members. The simulated atmospheric CO2 concentrations exhibit a stable carbon cycle over the pre-industrial era with multi-centennial variations somewhat smaller than in the observational records. Early land-cover changes have modulated atmospheric CO2 concentrations only slightly. We provide a model-based quantification of the sensitivity (termed γ) of the global carbon cycle to temperature for a variety of climate and forcing conditions. We diagnose a distinct dependence of γ on the forcing strength and time-scales involved, thus providing a possible explanation for the systematic difference in the observational estimates for different segments of the last millennium.

We apply different aerosol and aerosol precursor emission scenarios reflecting possible future control strategies for air pollution in the ECHAM5-HAM model, and simulate the resulting effect on the Earth's radiation budget. We use two opposing future mitigation strategies for the year 2030: one in which emission reduction legislation decided in countries throughout the world are effectively implemented (current legislation; CLE 2030) and one in which all technical options for emission reductions are being implemented independent of their cost (maximum feasible reduction; MFR 2030). We consider the direct, semi-direct and indirect radiative effects of aerosols. The total anthropogenic aerosol radiative forcing defined as the difference in the top-of-the-atmosphere radiation between 2000 and pre-industrial times amounts to −2.00 W/m2. In the future this negative global annual mean aerosol radiative forcing will only slightly change (+0.02 W/m2) under the "current legislation" scenario. Regionally, the effects are much larger: e.g. over Eastern Europe radiative forcing would increase by +1.50 W/m2 because of successful aerosol reduction policies, whereas over South Asia it would decrease by −1.10 W/m2 because of further growth of emissions. A "maximum feasible reduction" of aerosols and their precursors would lead to an increase of the global annual mean aerosol radiative forcing by +1.13 W/m2. Hence, in the latter case, the present day negative anthropogenic aerosol forcing could be more than halved by 2030 because of aerosol reduction policies and climate change thereafter will be to a larger extent be controlled by greenhouse gas emissions. We combined these two opposing future mitigation strategies for a number of experiments focusing on different sectors and regions. In addition, we performed sensitivity studies to estimate the importance of future changes in oxidant concentrations and the importance of the aerosol microphysical coupling within the range of expected future changes. For changes in oxidant concentrations caused by future air pollution mitigation, we do not find a significant effect for the global annual mean radiative aerosol forcing. In the extreme case of only abating SO2 or carbonaceous emissions to a maximum feasible extent, we find deviations from additivity for the radiative forcing over anthropogenic source regions up to 10% compared to an experiment abating both at the same time.

The double-moment cloud microphysics scheme from ECHAM4 that predicts both the mass mixing ratios and number concentrations of cloud droplets and ice crystals has been coupled to the size-resolved aerosol scheme ECHAM5-HAM. ECHAM5-HAM predicts the aerosol mass, number concentrations and mixing state. The simulated liquid, ice and total water content and the cloud droplet and ice crystal number concentrations as a function of temperature in stratiform mixed-phase clouds between 0 and −35° C agree much better with aircraft observations in the ECHAM5 simulations. ECHAM5 performs better because more realistic aerosol concentrations are available for cloud droplet nucleation and because the Bergeron-Findeisen process is parameterized as being more efficient. The total anthropogenic aerosol effect includes the direct, semi-direct and indirect effects and is defined as the difference in the top-of-the-atmosphere net radiation between present-day and pre-industrial times. It amounts to −1.9 W m−2 in ECHAM5, when a relative humidity dependent cloud cover scheme and aerosol emissions representative for the years 1750 and 2000 from the AeroCom emission inventory are used. The contribution of the cloud albedo effect amounts to −0.7 W m−2. The total anthropogenic aerosol effect is larger when either a statistical cloud cover scheme or a different aerosol emission inventory are employed because the cloud lifetime effect increases.

The evolution of the global aerosol system from 1860 to 2100 is investigated through a transient atmosphere-ocean General Circulation Model climate simulation with interactively coupled atmospheric aerosol and oceanic biogeochemistry modules. The microphysical aerosol module HAM incorporates the major global aerosol cycles with prognostic treatment of their composition, size distribution, and mixing state. Based on an SRES A1B emission scenario, the global mean sulfate burden is projected to peak in 2020 while black carbon and particulate organic matter show a lagged peak around 2070. From present day to future conditions the anthropogenic aerosol burden shifts generally from the northern high-latitudes to the developing low-latitude source regions with impacts on regional climate. Atmospheric residence- and aging-times show significant alterations under varying climatic and pollution conditions. Concurrently, the aerosol mixing state changes with an increasing aerosol mass fraction residing in the internally mixed accumulation mode. The associated increase in black carbon causes a more than threefold increase of its co-single scattering albedo from 1860 to 2100. Mid-visible aerosol optical depth increases from pre-industrial times, predominantly from the aerosol fine fraction, peaks at 0.26 around the sulfate peak in 2020 and maintains a high level thereafter, due to the continuing increase in carbonaceous aerosols. The global mean anthropogenic top of the atmosphere clear-sky short-wave direct aerosol radiative perturbation intensifies to −1.1 W m−2 around 2020 and weakens after 2050 to −0.6 W m−2, owing to an increase in atmospheric absorption. The demonstrated modifications in the aerosol residence- and aging-times, the microphysical state, and radiative properties challenge simplistic approaches to estimate the aerosol radiative effects from emission projections.

This paper describes a set of sensitivity experiments to several formulations of orography. Three sets are considered: a "Standard" orography consisting of an envelope orography produced originally for the ECMWF model, a"Navy" orography directly from the US Navy data and a "Scripps" orography based on the data set originally compiled several years ago at Scripps. The last two are mean orographies which do not use the envelope enhancement. A new filtering technique for handling the problem of Gibbs oscillations in spectral models has been used to produce the "Navy" and "Scripps" orographies, resulting in smoother fields than the "Standard" orography. The sensitivity experiments show that orography is still an important factor in controlling the model performance even in this class of models that use a semi-lagrangian formulation for water vapour, that in principle should be less sensitive to Gibbs oscillations than the Eulerian formulation. The largest impact can be seen in the stationary waves (asymmetric part of the geopotential at 500 mb) where the differences in total height and spatial pattern generate up to 60 m differences, and in the surface fields where the Gibbs removal procedure is successful in alleviating the appearance of unrealistic oscillations over the ocean. These results indicate that Gibbs oscillations also need to be treated in this class of models. The best overall result is obtained using the "Navy" data set, that achieves a good compromise between amplitude of the stationary waves and smoothness of the surface fields.