Max Planck Institute Earth System Model Research Articles

Atmosphere–ocean–aerosol–chemistry–climate model SOCOLv4.0: description and evaluation

Abstract. This paper features the new atmosphere–ocean–aerosol–chemistry–climate model, SOlar Climate Ozone Links (SOCOL) v4.0, and its validation. The new model was built by interactively coupling the Max Planck Institute Earth System Model version 1.2 (MPI-ESM1.2) (T63, L47) with the chemistry (99 species) and size-resolving (40 bins) sulfate aerosol microphysics modules from the aerosol–chemistry–climate model, SOCOL-AERv2. We evaluate its performance against reanalysis products and observations of atmospheric circulation, temperature, and trace gas distribution, with a focus on stratospheric processes. We show that SOCOLv4.0 captures the low- and midlatitude stratospheric ozone well in terms of the climatological state, variability and evolution. The model provides an accurate representation of climate change, showing a global surface warming trend consistent with observations as well as realistic cooling in the stratosphere caused by greenhouse gas emissions, although, as in previous model versions, a too-fast residual circulation and exaggerated mixing in the surf zone are still present. The stratospheric sulfur budget for moderate volcanic activity is well represented by the model, albeit with slightly underestimated aerosol lifetime after major eruptions. The presence of the interactive ocean and a successful representation of recent climate and ozone layer trends make SOCOLv4.0 ideal for studies devoted to future ozone evolution and effects of greenhouse gases and ozone-destroying substances, as well as the evaluation of potential solar geoengineering measures through sulfur injections. Potential further model improvements could be to increase the vertical resolution, which is expected to allow better meridional transport in the stratosphere, as well as to update the photolysis calculation module and budget of mesospheric odd nitrogen. In summary, this paper demonstrates that SOCOLv4.0 is well suited for applications related to the stratospheric ozone and sulfate aerosol evolution, including its participation in ongoing and future model intercomparison projects.

Improving seasonal predictions of meteorological drought by conditioning on ENSO states

Useful hindcast skill of meteorological drought, assessed with the 3-month standardized precipitation index (SPI), has been so far limited to one lead month (time horizon of the prediction). Here, we quadruple that lead time by demonstrating useful skill up to lead month 4. To obtain useful hindcast skill of meteorological drought at these long lead times, we exploit well-known El Niño-Southern Oscillation (ENSO)–precipitation teleconnections through ENSO-state conditioning. We condition initialized seasonal SPI hindcasts, derived from the Max-Planck-Institute Earth System Model (MPI-ESM) over the period 1982–2013, on ENSO states by exploring significant agreements between two complementary analyses: hindcast skill ENSO–composites, and observed ENSO–precipitation correlations. Such conditioned hindcast skill of meteorological drought is in MPI-ESM significant and reliable for lead months 2 to 4 in equatorial South America and southern North America during these regions’ dry ENSO phases. When a region’s dry ENSO phase is present at the initialization in autumn (ASO), predictions of meteorological drought show useful hindcast skill for the upcoming winter (DJF) in the respective region. The area of this useful hindcast skill is further enlarged in both regions when the respective region’s dry ENSO phase is already present in the antecedent summer (conditioning on ENSO states in JJA). Active ENSO events constitute windows of opportunity for drought predictions that are insufficiently covered by typical predictability analyses. For these windows, we demonstrate predictive skill at unprecedented lead times with a single model whose output is not bias corrected. This contribution exemplifies the value of ENSO-state conditioning in identifying these windows of opportunity for regions that are arguably most affected by ENSO–precipitation teleconnections. During these regions’ dry ENSO phases, reliable predictive skill of meteorological drought is at long lead times particularly valuable and moves the frontier of meteorological drought predictions.

Incorporating the stable carbon isotope <sup>13</sup>C in the ocean biogeochemical component of the Max Planck Institute Earth System Model

Abstract. The stable carbon isotopic composition (δ13C) is an important variable to study the ocean carbon cycle across different timescales. We include a new representation of the stable carbon isotope 13C into the HAMburg Ocean Carbon Cycle model (HAMOCC), the ocean biogeochemical component of the Max Planck Institute Earth System Model (MPI-ESM). 13C is explicitly resolved for all oceanic carbon pools considered. We account for fractionation during air–sea gas exchange and for biological fractionation ϵp associated with photosynthetic carbon fixation during phytoplankton growth. We examine two ϵp parameterisations of different complexity: ϵpPopp varies with surface dissolved CO2 concentration (Popp et al., 1989), while ϵpLaws additionally depends on local phytoplankton growth rates (Laws et al., 1995). When compared to observations of δ13C of dissolved inorganic carbon (DIC), both parameterisations yield similar performance. However, with regard to δ13C in particulate organic carbon (POC) ϵpPopp shows a considerably improved performance compared to ϵpLaws. This is because ϵpLaws produces too strong a preference for 12C, resulting in δ13CPOC that is too low in our model. The model also well reproduces the global oceanic anthropogenic CO2 sink and the oceanic 13C Suess effect, i.e. the intrusion and distribution of the isotopically light anthropogenic CO2 in the ocean. The satisfactory model performance of the present-day oceanic δ13C distribution using ϵpPopp and of the anthropogenic CO2 uptake allows us to further investigate the potential sources of uncertainty of the Eide et al. (2017a) approach for estimating the oceanic 13C Suess effect. Eide et al. (2017a) derived the first global oceanic 13C Suess effect estimate based on observations. They have noted a potential underestimation, but their approach does not provide any insight about the cause. By applying the Eide et al. (2017a) approach to the model data we are able to investigate in detail potential sources of underestimation of the 13C Suess effect. Based on our model we find underestimations of the 13C Suess effect at 200 m by 0.24 ‰ in the Indian Ocean, 0.21 ‰ in the North Pacific, 0.26 ‰ in the South Pacific, 0.1 ‰ in the North Atlantic and 0.14 ‰ in the South Atlantic. We attribute the major sources of underestimation to two assumptions in the Eide et al. (2017a) approach: the spatially uniform preformed component of δ13CDIC in year 1940 and the neglect of processes that are not directly linked to the oceanic uptake and transport of chlorofluorocarbon-12 (CFC-12) such as the decrease in δ13CPOC over the industrial period. The new 13C module in the ocean biogeochemical component of MPI-ESM shows satisfying performance. It is a useful tool to study the ocean carbon sink under the anthropogenic influences, and it will be applied to investigating variations of ocean carbon cycle in the past.

Detectability of Artificial Ocean Alkalinization and Stratospheric Aerosol Injection in MPI‐ESM

AbstractTo monitor the success of carbon dioxide removal (CDR) or solar radiation management (SRM) that offset anthropogenic climate change, the forced response to any external forcing is required to be detectable against internal variability. Thus far, only the detectability of SRM has been examined using both a stationary and nonstationary detection and attribution method. Here, the spatiotemporal detectability of the forced response to artificial ocean alkalinization (AOA) and stratospheric aerosol injection (SAI) as exemplary methods for CDR and SRM, respectively, is compared in Max Planck Institute Earth System Model (MPI‐ESM) experiments using regularized optimal fingerprinting and single‐model estimates of internal variability, while working under a stationary or nonstationary null hypothesis. Although both experiments are forced by emissions according to the Representative Concentration Pathway 8.5 (RCP8.5) and target the climate of the RCP4.5 scenario using AOA or SAI, detection timescales reflect the fundamentally different forcing agents. Moreover, detectability timescales are sensitive to the choice of null hypothesis. Globally, changes in the CO2 system in seawater are detected earlier than the response in temperature to AOA but later in the case of SAI. Locally, the detection time scales depend on the physical, chemical, and radiative impacts of CDR and SRM forcing on the climate system, as well as patterns of internal variability, which is highlighted for oceanic heat and carbon storage.

The climate change signal in the Mediterranean Sea in a regionally coupled atmosphere–ocean model

Abstract. We analyze the climate change signal in the Mediterranean Sea using the regionally coupled model REMO–OASIS–MPIOM (ROM; abbreviated from the regional atmosphere model, the OASIS3 coupler and the Max Planck Institute Ocean Model). The ROM oceanic component is global with regionally high horizontal resolution in the Mediterranean Sea so that the water exchanges with the adjacent North Atlantic and Black Sea are explicitly simulated. Simulations forced by ERA-Interim show an accurate representation of the present Mediterranean climate. Our analysis of the RCP8.5 (representative concentration pathway) scenario using the Max Planck Institute Earth System Model shows that the Mediterranean waters will be warmer and saltier throughout most of the basin by the end of this century. In the upper ocean layer, temperature is projected to have a mean increase of 2.7 ∘C, while the mean salinity will increase by 0.2 psu, presenting a decreasing trend in the western Mediterranean in contrast to the rest of the basin. The warming initially takes place at the surface and propagates gradually to deeper layers. Hydrographic changes have an impact on intermediate water characteristics, potentially affecting the Mediterranean thermohaline circulation in the future.

Forecast of AMD Quantity by a Series Tank Model in Three Stages: Case Studies in Two Closed Japanese Mines

There are about 100 sites of acid mine drainage (AMD) from abandoned/closed mines in Japan. For their sustainable treatment, future prediction of AMD quantity is crucial. In this study, AMD quantity was predicted for two closed mines in Japan based on a series tank model in three stages. The tank model parameters were determined from the relationship between the observed AMD quantity and the inflow of rainfall and snowmelt by using the Kalman filter and particle swarm optimization methods. The Automated Meteorological Data Acquisition System (AMeDAS) data of rainfall were corrected for elevation and by the statistical daily fluctuation model. The snowmelt was estimated from the AMeDAS data of rainfall, temperature, and sunshine duration by using mass and heat balance of snow. Fitting with one year of daily data was sufficient to obtain the AMD quantity model. Future AMD quantity was predicted by the constructed model using the forecast data of rainfall and temperature proposed by the Max Planck Institute–Earth System Model (MPI–ESM), based on the Intergovernmental Panel on Climate Change (IPCC) representative concentration pathway (RCP) 2.6 and RCP8.5 scenarios. The results showed that global warming causes an increase in the quantity and fluctuation of AMD, especially for large reservoirs and residence time of AMD. There is a concern that for mines with large AMD quantities, AMD treatment will be unstable due to future global warming.

Changes of Decadal SST Variations in the Subpolar North Atlantic under Strong CO2 Forcing as an Indicator for the Ocean Circulation’s Contribution to Atlantic Multidecadal Variability

AbstractThe question of whether ocean dynamics are relevant for basin-scale North Atlantic decadal temperature variability is the subject of ongoing discussions. Here, we analyze a set of simulations with a single climate model consisting of a 2000-yr preindustrial control experiment, a 100-member historical ensemble, and a 100-member ensemble forced with an incremental CO2 increase by 1% yr−1. Compared to previous approaches, our setup offers the following advantages: First, the large ensemble size allows us to robustly separate internally and externally forced variability and to robustly detect statistical links between different quantities. Second, the availability of different scenarios allows us to investigate the role of the background state for drivers of the variability. We find strong evidence that ocean dynamics, particularly ocean heat transport variations, form an important contribution to generate the Atlantic multidecadal variability (AMV) in the Max Planck Institute Earth System Model (MPI-ESM). Particularly the northwest North Atlantic is substantially affected by ocean circulation for the historical and preindustrial simulations. Anomalies of the Labrador Sea deep ocean density precede a change of the Atlantic meridional overturning circulation (AMOC) and heat advection to the region south of Greenland. Under strong CO2 forcing, the AMV–SST regression pattern shows crucial changes: SST variability in the northwestern part of the North Atlantic is strongly reduced, so that the AMV pattern in this scenario is dominated by the low-latitude branch. We found a connection to changes in the deep-water formation that cause a strong reduction of the mean AMOC and its variability. Consequently, ocean heat transport convergence becomes less important for the SST variability south of Greenland.

The Forced Response of the El Niño–Southern Oscillation–Indian Monsoon Teleconnection in Ensembles of Earth System Models

AbstractWe study the teleconnection between El Niño–Southern Oscillation (ENSO) and the Indian summer monsoon (IM) in large ensemble simulations, the Max Planck Institute Earth System Model (MPI-ESM), and the Community Earth System Model (CESM1). We characterize ENSO by the June–August Niño-3 box-average SST and the IM by the June–September average precipitation over India, and define their teleconnection in a changing climate as an ensemble-wise correlation. To test robustness, we also consider somewhat different variables that can characterize ENSO and the IM. We utilize ensembles converged to the system’s snapshot attractor for analyzing possible changes in the teleconnection. Our main finding is that the teleconnection strength is typically increasing on the long term in view of appropriately revised ensemble-wise indices. Indices involving a more western part of the Pacific reveal, furthermore, a short-term but rather strong increase in strength followed by some decrease at the turn of the century. Using the station-based Southern Oscillation index (SOI) as opposed to area-based indices leads to the identification of somewhat more erratic trends, but the turn-of-the-century “bump” is well detectable with it. All this is in contrast, if not in contradiction, to the discussion in the literature of a weakening teleconnection in the late twentieth century. We show here that this discrepancy can be due to any of three reasons: 1) ensemble-wise and temporal correlation coefficients used in the literature are different quantities; 2) the temporal moving correlation has a high statistical variability but possibly also persistence; or 3) MPI-ESM does not represent the Earth system faithfully.

Improvement in the decadal prediction skill of the North Atlantic extratropical winter circulation through increased model resolution

Abstract. In this study the latest version of the MiKlip decadal hindcast system is analyzed, and the effect of an increased horizontal and vertical resolution on the prediction skill of the extratropical winter circulation is assessed. Four different metrics – the storm track, blocking, cyclone and windstorm frequencies – are analyzed in the North Atlantic and European region. The model bias and the deterministic decadal hindcast skill are evaluated in ensembles of five members in a lower-resolution version (LR, atm: T63L47, ocean: 1.5∘ L40) and a higher-resolution version (HR, atm: T127L95, ocean: 0.4∘ L40) of the MiKlip system based on the Max Planck Institute Earth System model (MPI-ESM). The skill is assessed for the lead winters 2–5 in terms of the anomaly correlation of the quantities' winter averages using initializations between 1978 and 2012. The deterministic predictions are considered skillful if the anomaly correlation is positive and statistically significant. While the LR version shows common shortcomings of lower-resolution climate models, e.g., a storm track that is too zonal and southward displaced as well as a negative bias of blocking frequencies over the eastern North Atlantic and Europe, the HR version counteracts these biases. Cyclones, i.e., their frequencies and characteristics like strength and lifetime, are particularly better represented in HR. As a result, a chain of significantly improved decadal prediction skill between all four metrics is found with the increase in the spatial resolution. While the skill of the storm track is significantly improved primarily over the main source region of synoptic activity – the North Atlantic Current – the other extratropical quantities experience a significant improvement primarily downstream thereof, i.e., in regions where the synoptic systems typically intensify. Thus, the skill of the cyclone frequencies is significantly improved over the central North Atlantic and northern Europe, the skill of the blocking frequencies is significantly improved over the Mediterranean, Scandinavia and eastern Europe, and the skill of the windstorms is significantly improved over Newfoundland and central Europe. Not only is the skill improved with the increase in resolution, but the HR system itself also exhibits significant skill over large areas of the North Atlantic and European sector for all four circulation metrics. These results are particularly promising regarding the high socioeconomic impact of European winter windstorms and blocking situations.

Time Series Analysis of CMIP5 Model and Observed Sea Surface Temperature Anomaly Along Indian Coastal Zones

Joseph, D.; Liya, V.B.; Rojith, G.; Zacharia, P.U., and Grinson, G., 2019. Time series analysis of CMIP5 model and observed sea surface temperature anomaly along Indian coastal zones. In: Jithendran, K.P.; Saraswathy, R.; Balasubramanian, C.P.; Kumaraguru Vasagam, K.P.; Jayasankar, V.; Raghavan, R.; Alavandi, S.V., and Vijayan, K.K. (eds.), BRAQCON 2019: World Brackishwater Aquaculture Conference. Journal of Coastal Research, Special Issue No. 86, pp. 239–247. Coconut Creek (Florida), ISSN 0749-0208.Analysis of the time series Sea Surface Temperature (SST) variations is a key element in understanding the climate change impacts on the phenology, trophodynamics, distribution, and catch of commercial marine fish species. As SST projections are mostly model dependent, the discrepancy in the model and observed values needs to be elucidated so as to derive accurate interpretations and conclusions. In this paper, the decadal and seasonal variations of SST anomaly over the four coastal zones of India were analysed. The selected period of the study is from 1968 to 2017 and the data obtained from the International Comprehensive Ocean Atmosphere Data Sets (ICOADS) and Max Planck Institute Earth System Model (MPI-ESM) was used. The Coupled Model Intercomparison Project 5 (CMIP5) model data were used for comparison with the observed value. The linear trend of observed and modelled values was inferred using the least square method for four seasons. The seasonal variation of SST anomaly of observed data in four coastal zones of India reveals that Northeast zone exhibits the least trend of warming in all seasons, whereas Northwest zone shows the highest trend of warming. SST anomaly of the first decade in all zones exhibited negative values in all seasons, while last decade shows positive value, indicating warming trend. In other decades, no uniform SST trend for all seasons was observed. However, the Northeast zone is exceptional with negative anomaly during all decades. The warming trend was observed in all coasts for both model and observed values. The model and observed SST anomaly follows an almost similar trend, but with noticeable differences in values among both. Owing to the differences in the model and observed values, it could be emphasised that error corrections needs to be applied in futuristic SST projections and related studies of Indian fisheries.

Assessing climate change impacts on water resources in the Benue River Basin, Northern Cameroon

The Benue River Basin (BRB) is a major tributary of the Niger River Basin (NRB) and the second-largest river in Cameroon. It serves many water resource functions including irrigation, hydroelectricity production and navigation. Previous research has indicated that recent climate change (CC) has had significant impacts on local and regional hydrological regimes of this watershed. In this study, CC scenarios were integrated with a hydrological model to evaluate the influence of CC on water resources in the BRB. Historical and projected scenarios of dynamically downscaled temperature and precipitation from the REMO regional climate model (RCM) forced by the boundary conditions data of the Europe-wide Consortium Earth System Model (EC-ESM) and the Max Planck Institute-Earth System Model (MPI-ESM) general circulation models (GCMs) were used. The historical runs of the REMO simulations were first evaluated after which downscaled temperature and precipitation data were used as input for the HBV-Light hydrological model to simulate water balance components. The mean climate and hydrological variables for the historical (1981–2005) and the two future periods (2041–2065 and 2071–2095) were compared to assess the potential impact of CC on water resources in the middle and late twenty-first century under three greenhouse gases (GHGs) concentration scenarios, the Representative Concentration Pathways (RCPs) 2.6, 4.5 and 8.5. Our results show that (a) the HBV-Light hydrological model could effectively simulate the streamflow change in the BRB; (b) annual precipitation will decrease between 1 and 10% while both annual temperature and potential evapotranspiration (PET) will increase between 8–18 and 6–30%, respectively, under both scenarios, models and future periods; c) the combination of reduced precipitation and increase of PET results in a significant decrease in streamflow in the BRB (up to 51%) and this will move the basin to a much drier environmental state. Therefore, CC adaptation strategies and future development planning in this region must consider these important decreases of discharge.

Different response of surface temperature and air temperature to deforestation in climate models

Abstract. When quantifying temperature changes induced by deforestation (e.g., cooling in high latitudes, warming in low latitudes), satellite data, in situ observations, and climate models differ concerning the height at which the temperature is typically measured/simulated. In this study the effects of deforestation on surface temperature, near-surface air temperature, and lower atmospheric temperature are compared by analyzing the biogeophysical temperature effects of large-scale deforestation in the Max Planck Institute Earth System Model (MPI-ESM) separately for local effects (which are only apparent at the location of deforestation) and nonlocal effects (which are also apparent elsewhere). While the nonlocal effects (cooling in most regions) influence the temperature of the surface and lowest atmospheric layer equally, the local effects (warming in the tropics but a cooling in the higher latitudes) mainly affect the temperature of the surface. In agreement with observation-based studies, the local effects on surface and near-surface air temperature respond differently in the MPI-ESM, both concerning the magnitude of local temperature changes and the latitude at which the local deforestation effects turn from a cooling to a warming (at 45–55∘ N for surface temperature and around 35∘ N for near-surface air temperature). Subsequently, our single-model results are compared to model data from multiple climate models from the Climate Model Intercomparison Project (CMIP5). This inter-model comparison shows that in the northern midlatitudes, both concerning the summer warming and winter cooling, near-surface air temperature is affected by the local effects only about half as strongly as surface temperature. This study shows that the choice of temperature variable has a considerable effect on the observed and simulated temperature change. Studies about the biogeophysical effects of deforestation must carefully choose which temperature to consider.

Predictability of Multiyear Trends of the Pacific Decadal Oscillation in an MPI‐ESM Hindcast Ensemble

AbstractWhile hindcast skill for the Pacific Decadal Oscillation (PDO) has so far been limited to a few years, we present hindcast skill for PDO trends up to 10 years ahead. Our analysis is based on an initialized hindcast ensemble with the global Max Planck Institute Earth System Model (MPI‐ESM). As in previous studies, we find hindcast skill limited to a few years, when we first construct a lead‐year time series, from which we second calculate the PDO. We find similar hindcast skill when we first calculate the PDO for each start year and second construct a lead‐year time series. However, we find hindcast skill considerably increased, when we first calculate the PDO for each start year, second estimate multiyear trends, and third construct a lead‐year time series. Our results suggest hindcast skill for the low‐frequency variability of the PDO, which holds important implications for predictability analyses of other modes of long‐term climate variability.

Light absorption by marine cyanobacteria affects tropical climate mean state and variability

Abstract. Observations indicate that positively buoyant marine cyanobacteria, which are abundant throughout the tropical and subtropical ocean, have a strong local heating effect due to light absorption at the ocean surface. How these local changes in radiative heating affect the climate system on the large scale is unclear. We use the Max Planck Institute Earth System Model (MPI-ESM), include light absorption by cyanobacteria, and find a considerable cooling effect on tropical sea surface temperature (SST) in the order of 0.5 K on a climatological timescale. This cooling is caused by local shading of subtropical subsurface water by cyanobacteria that is upwelled at the Equator and in eastern boundary upwelling systems. Implications for the climate system include a westward shift of the Walker circulation and a weakening of the Hadley circulation. The amplitude of the seasonal cycle of SST is increased in large parts of the tropical ocean by up to 25 %, and the tropical Pacific interannual variability is enhanced by approx. 20 %. This study emphasizes the sensitivity of the tropical climate system to light absorption by cyanobacteria due to its regulative effect on tropical SST. Generally, including phytoplankton-dependent light attenuation instead of a globally uniform attenuation depth improves some of the major model temperature biases, indicating the relevance of taking this biophysical feedback into account in climate models.

Interactive ocean bathymetry and coastlines for simulating the last deglaciation with the Max Planck Institute Earth System Model (MPI-ESM-v1.2)

Abstract. As ice sheets grow or decay, the net flux of freshwater into the ocean changes and the bedrock adjusts due to isostatic adjustments, leading to variations in the bottom topography and the oceanic boundaries. This process was particularly intense during the last deglaciation due to the high rates of ice-sheet melting. It is, therefore, necessary to consider transient ocean bathymetry and coastlines when attempting to simulate the last deglaciation with Earth system models (ESMs). However, in most standard ESMs the land-sea mask is fixed throughout simulations because the generation of a new ocean model bathymetry implies several levels of manual corrections, a procedure that is hardly doable very often for long runs. This is one of the main technical problems towards simulating a complete glacial cycle with general circulation models. For the first time, we present a tool allowing for an automatic computation of bathymetry and land-sea mask changes in the Max Planck Institute Earth System Model (MPI-ESM). The algorithms developed in this paper can easily be adapted to any free-surface ocean model that uses the Arakawa-C grid in the horizontal and z-grid in the vertical including partial bottom cells. The strategy applied is described in detail and the algorithms are tested in a long-term simulation demonstrating the reliable behaviour. Our approach guarantees the conservation of mass and tracers at global and regional scales; that is, changes in a single grid point are only propagated regionally. The procedures presented here are an important contribution to the development of a fully coupled ice sheet–solid Earth–climate model system with time-varying topography and will allow for transient simulations of the last deglaciation considering interactive bathymetry and land-sea mask.

The climate of a retrograde rotating Earth

Abstract. To enhance understanding of Earth's climate, numerical experiments are performed contrasting a retrograde and prograde rotating Earth using the Max Planck Institute Earth system model. The experiments show that the sense of rotation has relatively little impact on the globally and zonally averaged energy budgets but leads to large shifts in continental climates, patterns of precipitation, and regions of deep water formation. Changes in the zonal asymmetries of the continental climates are expected given ideas developed more than a hundred years ago. Unexpected was, however, the switch in the character of the European–African climate with that of the Americas, with a drying of the former and a greening of the latter. Also unexpected was a shift in the storm track activity from the oceans to the land in the Northern Hemisphere. The different patterns of storms and changes in the direction of the trades influence fresh water transport, which may underpin the change of the role of the North Atlantic and the Pacific in terms of deep water formation, overturning and northward oceanic heat transport. These changes greatly influence northern hemispheric climate and atmospheric heat transport by eddies in ways that appear energetically consistent with a southward shift of the zonally and annually averaged tropical rain bands. Differences between the zonally averaged energy budget and the rain band shifts leave the door open, however, for an important role for stationary eddies in determining the position of tropical rains. Changes in ocean biogeochemistry largely follow shifts in ocean circulation, but the emergence of a “super” oxygen minimum zone in the Indian Ocean is not expected. The upwelling of phosphate-enriched and nitrate-depleted water provokes a dominance of cyanobacteria over bulk phytoplankton over vast areas – a phenomenon not observed in the prograde model. What would the climate of Earth look like if it would rotate in the reversed (retrograde) direction? Which of the characteristic climate patterns in the ocean, atmosphere, or land that are observed in a present-day climate are the result of the direction of Earth's rotation? Is, for example, the structure of the oceanic meridional overturning circulation (MOC) a consequence of the interplay of basin location and rotation direction? In experiments with the Max Planck Institute Earth system model (MPI-ESM), we investigate the effects of a retrograde rotation in all aspects of the climate system. The expected consequences of a retrograde rotation are reversals of the zonal wind and ocean circulation patterns. These changes are associated with major shifts in the temperature and precipitation patterns. For example, the temperature gradient between Europe and eastern Siberia is reversed, and the Sahara greens, while large parts of the Americas become deserts. Interestingly, the Intertropical Convergence Zone (ITCZ) shifts southward and the modeled double ITCZ in the Pacific changes to a single ITCZ, a result of zonal asymmetries in the structure of the tropical circulation. One of the most prominent non-trivial effects of a retrograde rotation is a collapse of the Atlantic MOC, while a strong overturning cell emerges in the Pacific. This clearly shows that the position of the MOC is not controlled by the sizes of the basins or by mountain chains splitting the continents in unequal runoff basins but by the location of the basins relative to the dominant wind directions. As a consequence of the changes in the ocean circulation, a “super” oxygen minimum zone develops in the Indian Ocean leading to upwelling of phosphate-enriched and nitrate-depleted water. These conditions provoke a dominance of cyanobacteria over bulk phytoplankton over vast areas, a phenomenon not observed in the prograde model.

Closing the energy balance using a canopy heat capacity and storage concept – a physically based approach for the land component JSBACHv3.11

Abstract. Land surface–atmosphere interaction is one of the most important characteristic for understanding the terrestrial climate system, as it determines the exchange fluxes of energy and water between the land and the overlying air mass. In several current climate models, it is common practice to use an unphysical approach to close the surface energy balance within the uppermost soil layer with finite thickness and heat capacity. In this study, a different approach is investigated by means of a physically based estimation of the canopy heat storage (SkIn+). Therefore, as a first step, results of an offline simulation of the land component JSBACH of the Max Planck Institute Earth system model (MPI-ESM) – constrained with atmospheric observations – are compared to energy fluxes and water fluxes derived from eddy covariance measurements observed at the CASES-99 field experiment in Kansas, where shallow vegetation prevails. This comparison of energy and evapotranspiration fluxes with observations at the site-level provides an assessment of the model's capacity to correctly reproduce the diurnal cycle. Following this, a global coupled land–atmosphere experiment is performed using an AMIP (Atmospheric Model Intercomparison Project) type simulation over 30 years to evaluate the regional impact of the SkIn+ scheme on a longer timescale, in particular, with respect to the effect of the canopy heat storage. The results of the offline experiment show that SkIn+ leads to a warming during the day and to a cooling at night relative to the old reference scheme, thereby improving the performance in the representation of the modeled surface fluxes on diurnal timescales. In particular: nocturnal heat releases unrealistically destroying the stable boundary layer disappear and phase errors are removed. On the global scale, for regions with no or low vegetation and a pronounced diurnal cycle, the nocturnal cooling prevails due to the fact that stable conditions at night maintain the delayed response in temperature, whereas the daytime turbulent exchange amplifies it. For the tropics and boreal forests as well as high latitudes, the scheme tends to warm the system.

Impact of SST diurnal cycle on ENSO asymmetry

The dominant mode of inter-annual variability in the tropical Pacific is El Niño-Southern Oscillation (ENSO). ENSO is not symmetric in the sense that El Niño is generally stronger than La Niña. However, many CMIP5 models, including the Max Planck Institute Earth System Model (MPI-ESM), produce an almost symmetric ENSO. This paper shows that, when resolving the intra-daily air-sea interactions by coupling the atmospheric and oceanic model components once per hour, the simulated ENSO asymmetry is improved. The improvement is closely related to the simulated diurnal cycle of sea surface tempreature (SST). In the central tropical Pacific, the simulated diurnal range of SST is about 0.2 °C, up to 10% of SST anomalies of the simulated ENSO events. During El Niño events, the diurnal cycle of SST anomalies enhances the atmospheric moist instability, whereby triggering more convection in the central tropical Pacific. During La Niña and normal years, however, the mean convection is not changed by the included diurnal cycle of SST anomalies. As a result, the anomalies of the trades, which are directly related to the convection, are stronger during El Niño years than that during La Niña years, making El Niño to be stronger than La Niña via Bjerknes feedback. These results obtained with a low resolution MPI-ESM are further confirmed by simulation with the same model at higher spatial resolutions, suggesting that the role of the intra-daily air-sea interactions for the ENSO asymmetry is independent of model resolutions.

Plant functional diversity affects climate–vegetation interaction

Abstract. We present how variations in plant functional diversity affect climate–vegetation interaction towards the end of the African Humid Period (AHP) in coupled land–atmosphere simulations using the Max Planck Institute Earth system model (MPI-ESM). In experiments with AHP boundary conditions, the extent of the “green” Sahara varies considerably with changes in plant functional diversity. Differences in vegetation cover extent and plant functional type (PFT) composition translate into significantly different land surface parameters, water cycling, and surface energy budgets. These changes have not only regional consequences but considerably alter large-scale atmospheric circulation patterns and the position of the tropical rain belt. Towards the end of the AHP, simulations with the standard PFT set in MPI-ESM depict a gradual decrease of precipitation and vegetation cover over time, while simulations with modified PFT composition show either a sharp decline of both variables or an even slower retreat. Thus, not the quantitative but the qualitative PFT composition determines climate–vegetation interaction and the climate–vegetation system response to external forcing. The sensitivity of simulated system states to changes in PFT composition raises the question how realistically Earth system models can actually represent climate–vegetation interaction, considering the poor representation of plant diversity in the current generation of land surface models.

The role of local sea surface temperature pattern changes in shaping climate change in the North Atlantic sector

Beside its global effects, climate change is manifested in many regionally pronounced features mainly resulting from changes in the oceanic and atmospheric circulation. Here we investigate the influence of the North Atlantic SST on shaping the winter-time response to global warming. Our results are based on a long-term climate projection with the Max Planck Institute Earth System Model (MPI-ESM) to investigate the influence of North Atlantic sea surface temperature pattern changes on shaping the atmospheric climate change signal. In sensitivity experiments with the model’s atmospheric component we decompose the response into components controlled by the local SST structure and components controlled by global/remote changes. MPI-ESM simulates a global warming response in SST similar to other climate models: there is a warming minimum—or ”warming hole”—in the subpolar North Atlantic, and the sharp SST gradients associated with the Gulf Stream and the North Atlantic Current shift northward by a few a degrees. Over the warming hole, global warming causes a relatively weak increase in rainfall. Beyond this, our experiments show more localized effects, likely resulting from future SST gradient changes in the North Atlantic. This includes a significant precipitation decrease to the south of the Gulf Stream despite increased underlying SSTs. Since this region is characterised by a strong band of precipitation in the current climate, this is contrary to the usual case that wet regions become wetter and dry regions become drier in a warmer climate. A moisture budget analysis identifies a complex interplay of various processes in the region of modified SST gradients: reduced surface winds cause a decrease in evaporation; and thermodynamic, modified atmospheric eddy transports, and coastal processes cause a change in the moisture convergence. The changes in the the North Atlantic storm track are mainly controlled by the non-regional changes in the forcing. The impact of the local SST pattern changes on regions outside the North Atlantic is small in our setup.