Search

The low-lying and densely populated Southeast Asia (SEA) region is threatened by sea-level change. To better understand the mechanism of sea-level change in this region, the sea-level trends and variability in the SEA region are investigated over the historical period 1950–2014 and during 2015–2099 under Shared Socio-economic Pathway5-8.5 scenario forcing, based on the output of the high-resolution (~ 0.1° for ocean) global climate model MPI-ESM-ER. The results confirm that the SEA sterodynamic sea level and its components (i.e., thermosteric, halosteric, and manometric sea level) are rising with accelerations, which are superimposed on El Niño-Southern Oscillation (ENSO) driven variabilities. To understand changes of thermosteric and halosteric sea level, regional-mean ocean heat and freshwater contents are analysed based on physical processes of ocean transports and air-sea fluxes. We show that ENSO variability impacts the thermosteric sea level mainly through lateral ocean heat transports, while it impacts the halosteric sea level mainly through surface freshwater flux. In the projection, a decreased volume transport of the relatively cold water (with respect to the SEA average) from the Pacific warms the SEA region. At the same time, a freshening of the transported saline water (with respect to the SEA average) results in an increased ocean freshwater transport into the SEA region. Locally, the pathway of volume transport from the Pacific to the SEA region is shifting northward, which results in a weakened Indonesian throughflow and an enhanced South China Sea throughflow, both leading to changes of regional sea-level pattern. [ABSTRACT FROM AUTHOR]

The pre-industrial millennium is among the periods selected by the Paleoclimate Model Intercomparison Project (PMIP) for experiments contributing to the sixth phase of the Coupled Model Intercomparison Project (CMIP6) and the fourth phase of the PMIP (PMIP4). The past1000 transient simulations serve to investigate the response to (mainly) natural forcing under background conditions not too different from today, and to discriminate between forced and internally generated variability on interannual to centennial timescales. This paper describes the motivation and the experimental set-ups for the PMIP4-CMIP6 past1000 simulations, and discusses the forcing agents orbital, solar, volcanic, and land use/land cover changes, and variations in greenhouse gas concentrations. The past1000 simulations covering the pre-industrial millennium from 850 Common Era (CE) to 1849 CE have to be complemented by historical simulations (1850 to 2014 CE) following the CMIP6 protocol. The external forcings for the past1000 experiments have been adapted to provide a seamless transition across these time periods. Protocols for the past1000 simulations have been divided into three tiers. A default forcing data set has been defined for the Tier 1 (the CMIP6 past1000) experiment. However, the PMIP community has maintained the flexibility to conduct coordinated sensitivity experiments to explore uncertainty in forcing reconstructions as well as parameter uncertainty in dedicated Tier 2 simulations. Additional experiments (Tier 3) are defined to foster collaborative model experiments focusing on the early instrumental period and to extend the temporal range and the scope of the simulations. This paper outlines current and future research foci and common analyses for collaborative work between the PMIP and the observational communities (reconstructions, instrumental data). [ABSTRACT FROM AUTHOR]

The pre-industrial millennium is among the periods selected by the Paleoclimate Model Intercomparison Project (PMIP) for experiments contributing to the sixth phase of the Coupled Model Intercomparison Project (CMIP6) and the fourth phase of PMIP (PMIP4). The past1000 transient simulations serve to investigate the response to (mainly) natural forcing under background conditions not too different from today, and to discriminate between forced and internally generated variability on interannual to centennial time scales. This manuscript describes the motivation and the experimental set-ups for the PMIP4-CMIP6 past1000 simulations, and discusses the forcing agents: orbital, solar, volcanic, land-use/land-cover changes, and variations in greenhouse gas concentrations. The past1000 simulations covering the pre-industrial millennium from 850 Common Era (CE) to 1849 CE have to be complemented by historical simulations (1850 to 2014 CE) following the CMIP6 protocol. The external forcings for the past1000 experiments have been adapted to provide a seamless transition across these time periods. Protocols for the past1000 simulations have been divided into three tiers. A default forcing data set has been defined for the "tier-1" (the CMIP6 past1000) experiment. However, the PMIP community has maintained the flexibility to conduct coordinated sensitivity experiments to explore uncertainty in forcing reconstructions as well as parameter uncertainty in dedicated "tier-2" simulations. Additional experiments ("tier-3") are defined to foster collaborative model experiments focusing on the early instrumental period and to extend the temporal range and the scope of the simulations. This manuscript outlines current and future research foci and common analyses for collaborative work between the PMIP and the observational communities (reconstructions, instrumental data). [ABSTRACT FROM AUTHOR]

We use the coupled climate model MPI-ESM to show that for higher CO levels the El Niño-Southern Oscillation (ENSO) and the Pacific Decadal Oscillation (PDO) merge into a single mode of Pacific variability, regardless of present-day or Middle Miocene (~15 Ma) topographic boundary conditions. Hence, topographic differences-determining the landscape of past climates-play a smaller role for Pacific variability than previously thought. We attribute the single variability mode to resonance between these two oscillation patterns. In order to estimate the strength of the resonance we compute the spectral power of the ENSO and PDO time series and their coherence. We find that for both Middle Miocene and present-day topographic conditions, higher CO forcing leads to stronger resonance between ENSO and PDO. Our results show that (1) stronger CO forcing enhances Pacific variability resulting in stronger 'atmospheric bridge' and that (2) past climates are likely to exhibit Pacific variability corresponding either to ENSO, PDO, or our proposed single mode. [ABSTRACT FROM AUTHOR]

Pronounced positive decadal-scale temperature anomalies occurred in the Arctic region in the first half of the twentieth century, an episode known as the early twentieth century warming (ETCW). Analyzing a 3,000-year unperturbed climate simulation performed with the Max Planck Institute Earth System Model, we demonstrate that internal variability of the Northern Hemisphere climate system is sufficient to reproduce warm events matching the observed ETCW. We perform a superposed epoch analysis on simulated data and identify 26 Arctic warming episodes compatible with the ETCW. The simulated events reproduce, in their ensemble average, magnitude as well as spatial and temporal extent of the observed ETCW. In individual realizations, the ETCW-like events indicate that different patterns of internally generated decadal Arctic warming are possible, including pan-Arctic warming events. We investigate the dynamics that typically lead to the simulated warming events: positive oceanic heat transport anomalies that form in the North Atlantic initialize the warming events and trigger an ocean-ice-albedo feedback in the Barents Sea region. The consequent reduction in sea-ice extent leads to enhanced multi-year surface warming through strengthened ocean heat release to the atmosphere. The oceanic heat transport anomalies reduce to pre-event levels around the year of the maximum warming. However, the warming events typically lasts for another 5-7 years until the sea-ice extent recovers to pre-event conditions. [ABSTRACT FROM AUTHOR]

The dependence of future regional sea level changes on ocean model resolution is investigated based on Max Planck Institute Earth System Model (MPI-ESM) simulations with varying spatial resolution, ranging from low resolution (LR), high resolution (HR), to eddy-rich (ER) resolution. Each run was driven by the shared socioeconomic pathway (SSP) 5-8.5 (fossil-fueled development) forcing. For each run the dynamic sea level (DSL) changes are evaluated by comparing the time mean of the SSP5-8.5 climate change scenario for the years 2080–99 to the time mean of the historical simulation for the years 1995–2014. Respective results indicate that each run reproduces previously identified large-scale DSL change patterns. However, substantial sensitivity of the projected DSL changes can be found on a regional to local scale with respect to model resolution. In comparison to models with parameterized eddies (HR and LR), enhanced sea level changes are found in the North Atlantic subtropical region, the Kuroshio region, and the Arctic Ocean in the model version capturing mesoscale processes (ER). Smaller yet still significant sea level changes can be found in the Southern Ocean and the North Atlantic subpolar region. These sea level changes are associated with changes in the regional circulation. Our study suggests that low-resolution sea level projections should be interpreted with care in regions where major differences are revealed here, particularly in eddy active regions such as the Kuroshio, Antarctic Circumpolar Current, Gulf Stream, and East Australian Current. Significance Statement: Sea level change is expected to be more realistic when mesoscale processes are explicitly resolved in climate models. However, century-long simulations with eddy-resolving models are computationally expensive. Therefore, current sea level projections are based on climate models in which ocean eddies are parameterized. The representation of sea level by these models considerably differs from actual observations, particularly in the eddy-rich regions such as the Southern Ocean and the western boundary currents, implying erroneous ocean circulation that affects the sea level projections. Taking this into account, we review the sea level change pattern in a climate model with featuring an eddy-rich ocean model and compare the results to state-of-the-art coarser-resolution versions of the same model. We found substantial DSL differences in the global ocean between the different resolutions. Relatively small-scale ocean eddies can hence have profound large-scale effects on the projected sea level which may affect our understanding of future sea level change as well as the planning of future investments to adapt to climate change around the world. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Changes in meridional heat transports, carried either by the atmosphere (HTRA) or by the ocean (HTRO), have been proposed to explain the decadal to multidecadal climate variations in the Arctic. On the other hand, model simulations indicate that, at high northern latitudes, variations in HTRA and HTRO are strongly coupled and may even compensate each other. A multi-century control integration with the Max Planck Institute global atmosphere-ocean model is analyzed to investigate the relative role of the HTRO and HTRA variations in shaping the Arctic climate and the consequences of their possible compensation. In the simulation, ocean heat transport anomalies modulate sea ice cover and surface heat fluxes mainly in the Barents Sea/Kara Sea region and the atmosphere responds with a modified pressure field. In response to positive HTRO anomalies there are negative HTRA anomalies associated with an export of relatively warm air southward to Western Siberia and a reduced inflow of heat over Alaska and northern Canada. While the compensation mechanism is prominent in this model, its dominating role is not constant over long time scales. The presence or absence of the compensation is determined mainly by the atmospheric circulation in the Pacific sector of the Arctic where the two leading large-scale atmospheric circulation patterns determine the lateral fluxes with varying contributions. The degree of compensation also determines the heat available to modulate the large-scale Arctic climate. The combined effect of atmospheric and oceanic contributions has to be considered to explain decadal-scale warming or cooling trends. [ABSTRACT FROM AUTHOR]

This study aims at improving the forecast skill of climate predictions through the use of ocean synthesis data for initial conditions of a coupled climate model. For this purpose, the coupled model of the Max Planck Institute (MPI) for Meteorology, which consists of the atmosphere model ECHAM5 and the MPI Ocean Model (MPI-OM), is initialized with oceanic synthesis fields available from the German contribution to Estimating the Circulation and Climate of the Ocean (GECCO) project. The use of an anomaly coupling scheme during the initialization avoids the main problems with drift in the climate predictions. Thus, the coupled model is continuously forced to follow the density anomalies of the GECCO synthesis over the period 1952–2001. Hindcast experiments are initialized from this experiment at constant intervals. The results show predictive skill through the initialization up to the decadal time scale, particularly over the North Atlantic. Viewed over the time scales analyzed here (annual, 5-yr, and 10-yr mean), greater skill for the North Atlantic sea surface temperature (SST) is obtained in the hindcast experiments than in either a damped persistence or trend forecast. The Atlantic meridional overturning circulation hindcast closely follows that of the GECCO oceanic synthesis. Hindcasts of global-mean temperature do not obtain greater skill than either damped persistence or a trend forecast, owing to the SST errors in the GECCO synthesis, outside the North Atlantic. An ensemble of forecast experiments is subsequently performed over the period 2002–11. North Atlantic SST from the forecast experiment agrees well with observations until the year 2007, and it is higher than if simulated without the oceanic initialization (averaged over the forecast period). The results confirm that both the initial and the boundary conditions must be accounted for in decadal climate predictions. [ABSTRACT FROM AUTHOR]

The meridional overturning circulation (MOC) in the coupled ECHAM5/MPIOM exhibits variability at periods of near 30 years and near 60 years. The 30-year variability, referred to as interdecadal variability (IDV), exist in an ocean model driven by climatological atmospheric forcing, suggesting that it is maintained by ocean dynamics; the 60-year variability, the multidecadal variability (MDV), is only observed in the fully coupled model and therefore is interpreted as an atmosphere–ocean coupled mode. The coexistence of the 30-year IDV and the 60-year MDV provides a possible explanation for the widespread time scales observed in climate variables. Further analyses of the climatologically forced ocean model shows that, the IDV is related to the interplay between the horizontal temperature-dominated density gradients and the ocean circulation: temperature anomalies move along the cyclonic subpolar gyre leading to fluctuations in horizontal density gradients and the subsequent weakening and strengthening of the MOC. This result is consistent with that from less complex models, indicating the robustness of the IDV. We further show that, along the North Atlantic Current path, the sea surface temperature anomalies are determined by the slow LSW advection at the intermediate depth. [ABSTRACT FROM AUTHOR]

This paper analyzes regional sea level changes in a climate change simulation using the Max Planck Institute for Meteorology (MPI) coupled atmosphere–ocean general circulation model ECHAM5/MPI-OM. The climate change scenario builds on observed atmospheric greenhouse gas (GHG) concentrations from 1860 to 2000, followed by the International Panel on Climate Change (IPCC) A1B climate change scenario until 2100; from 2100 to 2199, GHG concentrations are fixed at the 2100 level. As compared with the unperturbed control climate, global sea level rises 0.26 m by 2100, and 0.56 m by 2199 through steric expansion; eustatic changes are not included in this simulation. The model’s sea level evolves substantially differently among ocean basins. Sea level rise is strongest in the Arctic Ocean, from enhanced freshwater input from precipitation and continental runoff, and weakest in the Southern Ocean, because of compensation of steric changes through dynamic sea surface height (SSH) adjustments. In the North Atlantic Ocean (NA), a complex tripole SSH pattern across the subtropical to subpolar gyre front evolves, which is consistent with a northward shift of the NA current. On interannual to decadal time scales, the SSH difference between Bermuda and the Labrador Sea correlates highly with the combined baroclinic gyre transport in the NA but only weakly with the meridional overturning circulation (MOC) and, thus, does not allow for estimates of the MOC on these time scales. Bottom pressure increases over shelf areas by up to 0.45 m (water column equivalent) and decreases over the Atlantic section in the Southern Ocean by up to 0.20 m. The separate evaluation of thermosteric and halosteric sea level changes shows that thermosteric anomalies are positive over most of the World Ocean. Because of increased atmospheric moisture transport from low to high latitudes, halosteric anomalies are negative in the subtropical NA and partly compensate thermosteric anomalies, but are positive in the Arctic Ocean and add to thermosteric anomalies. The vertical distribution of thermosteric and halosteric anomalies is highly nonuniform among ocean basins, reaching deeper than 3000 m in the Southern Ocean, down to 2200 m in the North Atlantic, and only to depths of 500 m in the Pacific Ocean by the end of the twenty-first century. [ABSTRACT FROM AUTHOR]

Analyses of a 500-yr control integration with the non-flux-adjusted coupled atmosphere–sea ice–ocean model ECHAM5/Max-Planck-Institute Ocean Model (MPI-OM) show pronounced multidecadal fluctuations of the Atlantic overturning circulation and the associated meridional heat transport. The period of the oscillations is about 70–80 yr. The low-frequency variability of the meridional overturning circulation (MOC) contributes substantially to sea surface temperature and sea ice fluctuations in the North Atlantic. The strength of the overturning circulation is related to the convective activity in the deep-water formation regions, most notably the Labrador Sea, and the time-varying control on the freshwater export from the Arctic to the convection sites modulates the overturning circulation. The variability is sustained by an interplay between the storage and release of freshwater from the central Arctic and circulation changes in the Nordic Seas that are caused by variations in the Atlantic heat and salt transport. The relatively high resolution in the deep-water formation region and the Arctic Ocean suggests that a better representation of convective and frontal processes not only leads to an improvement in the mean state but also introduces new mechanisms determining multidecadal variability in large-scale ocean circulation. [ABSTRACT FROM AUTHOR]

A densely spaced hydrographic survey of the northern Irminger Basin together with satellite-tracked near-surface drifters confirm the intense mesoscale variability within and above the Denmark Strait overflow. In particular, the drifters show distinct cyclonic vortices over the downslope edge of the outflow plume. Growing perturbations such as these can be attributed to the baroclinic instability of a density current. A primitive equation model with periodic boundaries is used to simulate the destabilization of an idealized dense filament on a continental slope that resembles the northeastern Irminger Basin. Unstable waves evolve rapidly if the initial temperature profile is perturbed with a sinusoidal anomaly that exceeds a certain cutoff wavelength. As the waves grow to large amplitudes isolated eddies of both signs develop. Anticyclones form initially within the dense filament and are rich in overflow water. In contrast, cyclones form initially with their center in the ambient water but wrap outflow water around their center, thus containing a mixture of both water types. The nonlinear advection of waters that were originally located within the front between both water masses contributes most significantly to the stronger intensification of the cyclones in comparison with anticyclones. The frontal waters carry positive relative vorticity into the center of the cyclone. The process bears therefore some resemblance to atmospheric frontal cyclogenesis. After saturation there is a bottom jet of overflow water that is confined by counterrotating eddies: anticyclones upslope and cyclones downslope of the overflow core. The parameter dependence of the maximum growth rate is studied, and the implications of eddy-induced mixing for the water mass modification is discussed. [ABSTRACT FROM AUTHOR]

The early 19th century was the coldest period over the past 500 years, when strong tropical volcanic events and a solar minimum coincided. The 1809 unidentified eruption and the 1815 Tambora eruption happened consecutively during the Dalton minimum of solar irradiance; however, the relative role of the two forcing (volcano and solar) agents is still unclear. In this study, we examine the responses from a set of early 19th century simulations with combined and separated volcanic and solar forcing agents, as suggested in the protocol for the past1000 experiment of the Paleoclimate Modelling Intercomparison Project – Phase 4 (PMIP4). From 20-member ensemble simulations with the Max Planck Institute Earth system model (MPI-ESM1.2-LR), we find that the volcano- and solar-induced surface cooling is additive in the global mean/large scale, regardless of combining or separating the forcing agents. The two solar reconstructions (SATIRE (Spectral and Total Irradiance REconstruction-Millennia model) and PMOD (Physikalisch-Meteorologisches Observatorium Davos)) contribute to a cooling before and after 1815 of ∼0.05 and ∼0.15 K monthly average near-surface air cooling, respectively, indicating a limited solar contribution to the early 19th century cold period. The volcanic events provide the main cooling contributions, inducing a surface cooling that peaks at ∼0.82 K for the 1809 event and ∼1.35 K for Tambora. After the Tambora eruption, the temperature in most regions increases toward climatology largely within 5 years, along with the reduction of volcanic forcing. In the northern extratropical oceans, the temperature increases slowly at a constant rate until 1830, which is related to the reduction of seasonality and the concurrent changes in Arctic sea-ice extent. The albedo feedback of Arctic sea ice is found to be the main contributor to the Arctic amplification of the cooling signal. Several non-additive responses to solar and volcanic forcing happen on regional scales. In the atmosphere, the stratospheric polar vortex tends to strengthen when combining both volcano and solar forcing, even though the two forcing agents separately induce opposite-sign changes in stratospheric temperatures and zonal winds. In the ocean, when combining the two forcings, additional surface cold water propagates to the northern extratropics from the additional solar cooling in the tropics, which results in regional cooling along the propagation. Overall, this study not only quantifies the surface responses from combinations of the volcano and solar forcing, but also highlights the components that cannot be simply added from the responses of the individual forcing agents, indicating that a relatively small forcing agent (such as solar in early 19th century) can impact the response from the large forcing (such as the 1815 Tambora eruption) when considering regional climates. [ABSTRACT FROM AUTHOR]

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]

Current global warming results in an uptake of heat by the Earth system, which is distributed among the different components of the climate system. However, current‐generation climate models deliver heat inventory and partitioning estimates of Earth system components that differ from recent observations. Here we investigate the global heat distribution under warming by using fully‐coupled CMIP6 Earth system model experiments, including a version of the MPI‐ESM with a deep land model component, accommodating the required space for more realistic terrestrial heat storage. The results show that sufficiently deep land models exert increased subsurface land heat uptake, leading to a heat uptake partitioning among the Earth system components that is closer to observational estimates. The results are relevant for the understanding of Earth's heat partitioning and highlight the importance of the land heat sink in the Earth heat inventory. Plain Language Summary: Global warming is associated with heat accumulation in the Earth system due to the intensification of the greenhouse effect. The available heat is distributed unevenly throughout the climate subsystems: the ocean, land, atmosphere, and cryosphere. Overall, the current generation of climate models captures this partitioning well but, on average, shows an overestimation of the ocean heat uptake and an underestimation of the land heat uptake. Previous studies have shown that the lack of heat input into the land comes from structural limitations in the land model components used. In this study, we account for these shortcomings, which greatly improve the land heat uptake in simulations of future climate scenarios. We find that, as a result, the fraction of simulated heat taken up by the ocean is reduced. This leads to a heat distribution among the climate subsystems that is closer to observational estimates. Our results highlight that land heat uptake is relevant for the Earth system heat distribution and that future research should consider modeling approaches including a more realistic land heat sink. Key Points: Considering a sufficient land depth in current‐generation Earth system models facilitates realistic land heat uptake estimates under warmingIncreased land heat uptake leads to an adjusted global heat distribution and a reduced overestimation of simulated ocean heat uptakeThe land heat sink is relevant for the projection of the global energy inventory and its partitioning under global warming [ABSTRACT FROM AUTHOR]

The early 19th century was the coldest period over the past 500 years, when strong tropical volcanic events and a solar minimum coincided. The 1809 unidentified eruption and the 1815 Tambora eruption happened consecutively during the Dalton minimum of solar irradiance; however, the relative role of the two forcing (volcano and solar) agents is still unclear. In this study, we examine the effects from combinations of one volcanic with two different solar forcing reconstructions (SATIRE and PMOD) suggested in the protocol for the past1000 experiment of the Paleoclimate Modelling Intercomparison Project - Phase 4 (PMIP4) by simulating the early 19th century climate. From 20-member ensemble simulations with the Max Planck Institute Earth System Model (MPI-ESM1.2-LR), we find that the volcano- and solar-induced surface cooling is in general additive, regardless of combining or separating the forcing agents. The two solar reconstructions (SATIRE and PMOD) contribute on average ~0.05 K/month and ~0.15 K/month surface air cooling, respectively, indicating a limited solar contribution to the early 19th century cold period. The volcanic events provide the main cooling contributions, inducing a surface cooling peak of ~0.82 K for the 1809 event and ~1.35 K for Tambora. After the Tambora eruption, the cooling in most regions reduces largely within 5 years when a global cooling of ~0.34 K is reached, along with the reduction of volcanic forcing. In the northern extratropical oceans, the cooling reduces only slowly with a constant rate until 1830, which is related to the reduction of seasonality and the increased Arctic sea-ice extent. Also, the albedo feedback of Arctic sea ice is found to be the main contributor to the Arctic amplification of the cooling signal. Several non-additive responses to solar and volcanic forcing happen on regional scales. In the atmosphere, the polar vortex tends to strengthen when combining both volcano and solar forcing, even though the two forcing agents separately induce opposite responses. In the ocean, when combining the two forcings, additional surface cold water propagates to the northern extra-tropics from the additional solar cooling in the tropics, which results in regional cooling along the propagation. Overall, this study not only quantifies the surface responses from combinations of volcano and solar forcing, but also highlights the components that cannot be simply added from the responses of the individual forcing agents, indicating that a relatively small forcing agent (such as solar in early 19th century) can impact the response from the large forcing (such as the Tambora eruption) when considering regional climates. [ABSTRACT FROM AUTHOR]

The literature offers conflicting findings about which regions contribute most to increases in the global radiative feedback after a forcing increase. This paper explains the disagreement by discriminating between two common definitions of the local feedback, which use either local temperature or global temperature as their basis. Although the two definitions of feedback have been previously compared in aquaplanet models with slab oceans, here the definitions are compared for the first time in an atmosphere-ocean general circulation model (MPI-ESM1.2) integrated over four doublings of atmospheric CO2 concentrations. Large differences between the definitions can be seen in all feedbacks, but especially in the temperature and water vapor feedbacks. Differences of up to 10 W m-2 K-1 over the Southern Ocean can be explained by the pattern of surface warming, which weights the local feedbacks and reduces their contribution to the global mean. This finding is, however, dependent on the resolution of analysis, because the local-temperature definition is mathematically inconsistent across spatial scales. Furthermore, attempts to estimate the effect of "pattern weighting" by separating local feedbacks and warming patterns at the gridcell level fail, because the radiative change in key tropical regions is also determined by tropospheric stability via the global circulation. These findings indicate that studies of regional feedback change are more sensitive to methodological choices than previously thought, and that the tropics most likely dominate regional contributions to global radiative feedback change on decadal to centennial time scales. [ABSTRACT FROM AUTHOR]

Recent variability in West African monsoon rainfall (WAMR) has been shown to be influenced by multiple ocean–atmosphere modes, including the El Niño Southern Oscillation, Atlantic Multidecadal Oscillation and the Interdecadal Pacific Oscillation. How these modes will change in response to long term forcing is less well understood. Here we use four transient simulations driven by changes in orbital forcing and greenhouse gas concentrations over the past 6000 years to examine the relationship between West African monsoon rainfall multiscale variability and changes in the modes associated with this variability. All four models show a near linear decline in monsoon rainfall over the past 6000 years in response to the gradual weakening of the interhemispheric gradient in sea surface temperatures. The only indices that show a long-term trend are those associated with the strengthening of the El Niño Southern Oscillation from the mid-Holocene onwards. At the interannual-to-decadal timescale, WAMR variability is largely influenced by Pacific–Atlantic – Mediterranean Sea teleconnections in all simulations; the exact configurations are model sensitive. The WAMR interannual-to-decadal variability depicts marked multi-centennial oscillations, with La Niña/negative Pacific Decadal Oscillation and a weakening and/or poleward shift of subtropical high-pressure systems over the Atlantic favoring wet WAMR anomalies. The WAMR interannual-to-decadal variability also depicts an overall decreasing trend throughout the Holocene that is consistent among the simulations. This decreasing trend relates to changes in the North Atlantic and Gulf of Guinea Sea Surface Temperature variability. [ABSTRACT FROM AUTHOR]

The North Atlantic Ocean gyres play a crucial role in distributing heat, salinity and nutrients throughout the ocean basin. An important question is: how could these gyre circulations evolve under global warming? Using the Max Planck Institute Grand Ensemble (MPI-GE) simulations, we investigate the mechanisms behind the changes in the mean state of the gyre circulations in the North Atlantic under increased greenhouse gas forcing. We find two distinct phases in the changes of the gyre circulation. In the first phase of 2 K warming in the global mean surface temperature (GMST) from the pre-industrial values, the sub-polar gyre (SPG) intensifies. The intensification starts from the eastern side of the SPG and is controlled by the density difference between the east and the centre of the SPG. In contrast, in the next phase of further warming beyond 2 K, the SPG keeps constant in strength except the southern part. Moreover, a northward shift of the zonal winds leads to a northward shift of the sub-tropical gyre (STG). The latter causes warmer conditions in the inter-gyre region and a weakening of the southern part of the SPG. These changes in the mean state of the gyre circulations have associated signatures in the basin wide ocean heat transports and on the North Atlantic Ocean surface temperatures. [ABSTRACT FROM AUTHOR]

The anthropogenically-intensified greenhouse effect has caused a radiative imbalance at the top of the atmosphere during the industrial period. This, in turn, has led to an energy surplus in various components of the Earth system, with the ocean storing the largest part. The land contribution ranks second with the latest observational estimates based on borehole temperature profiles, which quantify the terrestrial energy surplus to be 6 % in the last five decades, whereas studies based on state-of-the-art climate models scale it down to 2 %. This underestimation stems from land surface models (LSMs) having a too shallow subsurface, which severely constrains the land heat uptake simulated by Earth System Models (ESMs). A forced simulation of the last 2000 years with the Max Planck Institute ESM (MPI-ESM) using a deep LSM captures 4 times more heat than the standard shallow MPI-ESM simulations in the historical period, well above the estimates provided by other ESMs. However, deepening the LSM does not remarkably affect the simulated surface temperature. It is shown that the heat stored during the historical period by an ESM using a deep LSM component can be accurately estimated by considering the surface temperatures simulated by the ESM using a shallow LSM and propagating them with a standalone forward model. This result is used to derive estimates of land heat uptake using all available observational datasets, reanalysis products, and state-of-the-art ESM experiments. This approach yields values of 10.5-16.0 ZJ for 1971-2018, slightly smaller than the latest borehole-based estimates (18.2 ZJ). [ABSTRACT FROM AUTHOR]

Small‐to‐moderate volcanic eruptions can lead to significant surface cooling when they occur clustered, as observed in recent decades. In this study, based on new high‐resolution ice‐core data from Greenland, we produce a new volcanic forcing data set that includes several small‐to‐moderate eruptions not included in prior reconstructions and investigate their climate impacts of the early 19th century through ensemble simulations with the Max Planck Institute Earth System Model. We find that clustered small‐to‐moderate eruptions produce significant additional global surface cooling (∼0.07 K) during the period 1812–1820, superposing with the cooling by large eruptions in 1809 (unidentified location) and 1815 (Tambora). This additional cooling helps explain the reconstructed long‐lasting cooling after the large eruptions, but simulated regional impacts cannot be confirmed with reconstructions due to a low signal‐to‐noise ratio. This study highlights the importance of small‐to‐moderate eruptions for climate simulations as their impacts can be comparable with that of solar irradiance changes. Plain Language Summary: Volcanic eruptions can influence global climate through the emission of sulfuric acids shielding Earth from incoming solar radiation. Previous volcanic reconstructions based on ice‐cores from the polar regions, however, only considered very strong volcanic eruptions. In this study, based on new ice‐core measurements from Greenland, we reconstruct for the first time volcanic sulfur emissions from small to medium‐sized eruptions and investigate their impact on climate in the early 19th century through experiments with the Max Planck Institute Earth System Model (MPI‐ESM1.2‐LR). We find that clustering of small to medium‐sized eruptions can cause significant global surface cooling (∼0.07 K), which during the 1812–1820 period amplified the cooling caused by the two known large eruptions of the period (1809 unidentified and 1815 Tambora). This additional surface cooling from small eruptions helps explain the long‐lasting cooling after the two strong eruptions generally found in the reconstruction, but the simulated regional impacts cannot be fully confirmed with reconstructions that are too noisy. This study highlights the importance of including small‐to‐moderate eruptions for climate model simulations as their impacts are comparable with that of solar irradiance forcing. Key Points: A new ice‐core based reconstruction of volcanic sulfate in the atmosphere (1733–1895) includes small‐to‐moderate eruptionsSmall‐to‐moderate eruptions can induce significant surface cooling and help explain the long‐lasting cooling in the early 19th centuryRegional cooling from small‐to‐moderate eruptions may be influenced by the circulation changes from the 1815 Tambora for over 10 years [ABSTRACT FROM AUTHOR]

The Northern Hemisphere mid‐latitudes will be exposed to hydroclimatic risk in next coming decades because the subtropical expansion. However, it is not clear when the anthropogenic signal will emerge from the internal climate variability. For this purpose, we investigate the time of emergence (ToE) of the hemispheric and regional shift of northern subtropical margins in the Max Planck Institute Grand Ensemble. For several indicators, the ToE of the poleward shift of Northern subtropical margin will not occur by the end of the 21st century, neither at regional nor at hemispheric scale. The exceptions are the Mediterranean/Middle East and, to a lesser degree, Western Pacific, where the ToE would occur earlier. According to our results, given the fundamental role played by internal variability, trends of Northern Hemisphere subtropical poleward shift that have been identified over last decades in reanalyses cannot be considered as robust signals of anthropogenic climate change. Plain Language Summary: The projected poleward shift of subtropical margins will affect densely populated areas of the Northern Hemisphere that have limited or scarce water resources. Therefore, it is fundamental for societies in these areas to know when the effects of climate change on the regional environment will become evident, because this information can make the difference for effective adaptation and mitigation strategies. The position of the subtropical margins is intimately related to the zonal‐mean Hadley Circulation, but is also affected by regional processes related to the distribution of continents, oceans and to their characteristics. In our study we show that the poleward shift of Northern subtropics has strong regional connotations and only in few areas, for example, the Mediterranean/Middle East and, to a lesser extent the Western Pacific, the anthropogenic climate change signal will emerge from the natural climate variability before the end of the 21st century. Key Points: Recent trends of subtropical margins are not a robust signal of anthropogenic climate change in the Northern HemisphereThe internal climate variability prevails over the forced response until 2100 on Northern subtropical margin poleward shift at hemispheric scaleThe poleward shift of Northern subtropics will affect mainly Mediterranean/Middle East and Western Pacific sectors [ABSTRACT FROM AUTHOR]

A vivid scientific debate exists on the nature of the Atlantic Multidecadal Variability (AMV) as an intrinsic rather than predominantly forced climatic phenomenon, and on the role of ocean circulation. Here, we use a multi-millennial unperturbed control simulation and a Holocene simulation with slow-varying greenhouse gas and orbital forcing performed with the low-resolution version of the Max Planck Institute Earth System Model to illustrate thermohaline conditions associated with twelve events of strong AMV that are comparable, in the surface anomalies, to observations in their amplitudes (~ 0.3 °C) and periods (~ 80 years). The events are associated with recurrent yet spatially diverse same-sign anomalous sea-surface temperature and salinity fields that are substantially symmetric in the warm-to-cold and following cold-to-warm transitions and only partly superpose with the long-term spatial AMV pattern. Subpolar cold-fresh anomalies develop in the deep layers during the peak cold phase of strong AMV events, often in association with subtropical warm-salty anomalies yielding a meridional dipole pattern. The Atlantic meridional overturning circulation (AMOC) robustly weakens during the warm-to-cold transition of a strong AMV event and recovers thereafter, with surface salinity anomalies being potential precursors of such overturning changes. A Holocene simulation with the same model including volcanic forcing can disrupt the intrinsic AMV–AMOC connection as post-eruption periods often feature an AMOC strengthening forced by the volcanically induced surface cooling. Overall, our results support the AMV as a potential intrinsic feature of climate, whose episodic strong anomalous events can display different shades of spatial patterns and timings for the warm-to-cold and subsequent cold-to-warm transitions. Attribution of historical AMV fluctuations thus requires full consideration of the associated surface and subsurface thermohaline conditions and assessing the AMOC–AMV relation. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

We develop a data assimilation scheme with the Icosahedral Non-hydrostatic Earth System Model (ICON-ESM) for operational decadal and seasonal climate predictions at the German weather service. For this purpose, we implement an Ensemble Kalman Filter to the ocean component as a first step towards a weakly coupled data assimilation. We performed an assimilation experiment over the period 1960–2014. This ocean-only assimilation experiment serves to initialize 10-year long retrospective predictions (hindcasts) started each year on 1 November. On multi-annual time scales, we find predictability of sea surface temperature and salinity as well as oceanic heat and salt contents especially in the North Atlantic. The mean Atlantic Meridional Overturning Circulation is realistic and the variability is stable during the assimilation. On seasonal time scales, we find high predictive skill in the tropics with highest values in variables related to the El Niño/Southern Oscillation phenomenon. In the Arctic, the hindcasts correctly represent the decreasing sea ice trend in winter and, to a lesser degree, also in summer, although sea ice concentration is generally much too low in both hemispheres in summer. However, compared to other prediction systems, prediction skill is relatively low in regions apart from the tropical Pacific due to the missing atmospheric assimilation. Further improvements of the simulated mean state of ICON-ESM, e.g. through fine-tuning of the sea ice and the oceanic circulation in the Southern Ocean, are expected to improve the predictive skill. In general, we demonstrate that our data assimilation method is successfully initializing the oceanic component of the climate system. [ABSTRACT FROM AUTHOR]

The Atlantic Meridional Overturning Circulation (AMOC) is a key feature of the North Atlantic with global ocean impacts. The AMOC's response to past changes in forcings during the Holocene provides important context for the coming centuries. Here, we investigate AMOC trends using an emerging set of transient simulations using multiple global climate models for the past 6,000 years. Although some models show changes, no consistent trend in overall AMOC strength during the mid‐to‐late Holocene emerges from the ensemble. We interpret this result to suggest no overall change in AMOC, which fits with our assessment of available proxy reconstructions. The decadal variability of the AMOC does not change in ensemble during the mid‐ and late‐Holocene. There are interesting AMOC changes seen in the early Holocene, but their nature depends a lot on which inputs are used to drive the experiment. Plain Language Summary: The Atlantic Meridional Overturning Circulation (AMOC) is a deep ocean circulation system that is both important for climate and vulnerable to climate changes. Here we use a set of multiple climate models to look at how the AMOC responded to changes in climate drivers over the past few thousand years. The changes are only small in all of the models, and do not agree in their direction. The AMOC naturally varies on decadal timescales, but we do not see any strong trends in its variability either. We consider these simulations to indicate that the overall AMOC has not changed over the past 6,000 years, which fits with recent data reconstructions. Key Points: A multi‐model ensemble of Holocene transient simulations by general circulation models has been assembledAlthough some models show changes, no consistent trend in overall Atlantic Meridional Overturning Circulation (AMOC) strength during the mid‐to‐late Holocene emerges from the ensembleWe interpret this result to suggest no overall change in AMOC, which fits with our assessment of available proxy reconstructions [ABSTRACT FROM AUTHOR]

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 as of yet. We use the Max Planck Institute Earth System Model (MPI-ESM) and find that - in contrast to the heating effect which was reported in previous studies - cyanobacteria have a considerable cooling effect on tropical climatological sea surface temperature (SST) in the order of 0.5 K. This cooling is caused by local shading of subtropical subsurface water that is upwelled at the equator and in eastern boundary upwelling systems. Implications for the climate system include an expansion of the Hadley cells and a westward shift of the Walker 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 ~20%. This study emphasizes the sensitivity of the tropical climate system to light absorption by the specific phytoplankton group of cyanobacteria due to its regulative effect on tropical SST. Generally, including the phytoplankton-dependent light attenuation instead of a globally uniform attenuation depth improves some of the major model temperature biases, indicating the relevance of taking into account this bio-physical feedback in climate models. [ABSTRACT FROM AUTHOR]

The Atlantic meridional overturning circulation (AMOC) is an important part of our climate system. The AMOC is predicted to weaken under climate change; however, theories suggest that it may have a tipping point beyond which recovery is difficult, hence showing quasi-irreversibility (hysteresis). Although hysteresis has been seen in simple models, it has been difficult to demonstrate in comprehensive global climate models. Here, we outline a set of experiments designed to explore AMOC hysteresis and sensitivity to additional freshwater input as part of the North Atlantic Hosing Model Intercomparison Project (NAHosMIP). These experiments include adding additional freshwater (hosing) for a fixed length of time to examine the rate and mechanisms of AMOC weakening and whether the AMOC subsequently recovers once hosing stops. Initial results are shown from eight climate models participating in the Sixth Coupled Model Intercomparison Project (CMIP6). The AMOC weakens in all models as a result of the freshening, but once the freshening ceases, the AMOC recovers in half of the models, and in the other half it stays in a weakened state. The difference in model behaviour cannot be explained by the ocean model resolution or type nor by details of subgrid-scale parameterisations. Likewise, it cannot be explained by previously proposed properties of the mean climate state such as the strength of the salinity advection feedback. Instead, the AMOC recovery is determined by the climate state reached when hosing stops, with those experiments where the AMOC is weakest not experiencing a recovery. [ABSTRACT FROM AUTHOR]

The effect of anthropogenic climate change in the ocean is challenging to project because atmosphere-ocean general circulation models (AOGCMs) respond differently to forcing. This study focuses on changes in the Atlantic Meridional Overturning Circulation (AMOC), ocean heat content (Δ OHC), and the spatial pattern of ocean dynamic sea level ( Δ ζ ). We analyse experiments following the FAFMIP protocol, in which AOGCMs are forced at the ocean surface with standardised heat, freshwater and momentum flux perturbations, typical of those produced by doubling CO 2 . Using two new heat-flux-forced experiments, we find that the AMOC weakening is mainly caused by and linearly related to the North Atlantic heat flux perturbation, and further weakened by a positive coupled heat flux feedback. The quantitative relationships are model-dependent, but few models show significant AMOC change due to freshwater or momentum forcing, or to heat flux forcing outside the North Atlantic. AMOC decline causes warming at the South Atlantic-Southern Ocean interface. It does not strongly affect the global-mean vertical distribution of Δ OHC, which is dominated by the Southern Ocean. AMOC decline strongly affects Δ ζ in the North Atlantic, with smaller effects in the Southern Ocean and North Pacific. The ensemble-mean Δ ζ and Δ OHC patterns are mostly attributable to the heat added by the flux perturbation, with smaller effects from ocean heat and salinity redistribution. The ensemble spread, on the other hand, is largely due to redistribution, with pronounced disagreement among the AOGCMs. [ABSTRACT FROM AUTHOR]

The North Atlantic subpolar gyre (SPG) plays a crucial role in determining the regional ocean surface temperature (SST), which has profound implications on the surrounding continental and coastal climate. Here, we analyze the Max Planck Institute-Grand Ensemble global warming experiments and show that the SPG can evolve in two distinct phases under continuous global warming. In the first phase, as the global mean surface temperature approaches 2-K warming, the eastern SPG intensifies in combination with a weakening Atlantic meridional overturning circulation (AMOC), accompanied by a cooling of subpolar North Atlantic SST, known as the warming hole. The associated oceanic fingerprint matches with the observations over the last 15 years, where an intensification and cooling of the eastern SPG is related to salinity reduction at the eastern side of the SPG. However, for further warming beyond 2 K, in spite of a continuous decline in the AMOC, a northward shift of the mean zonal wind extends the subtropical gyre northward with an associated disruption of the eastern SPG intensification, resulting in the cessation of the warming hole. Therefore, a shift from the initially dominating oceanic drivers to the atmospheric driver results into a two-phase evolution of the North Atlantic Ocean SPG circulation and the associated SST under continuous global warming. [ABSTRACT FROM AUTHOR]

The Atlantic meridional overturning circulation (AMOC) is an important part of our climate system. The AMOC is predicted to weaken under climate change, however there are theories that it may have a tipping point beyond which recovery is difficult, hence showing quasi-irreversibility (hysteresis). Although hysteresis has been seen in simple models, it has been difficult to demonstrate in comprehensive global climate models. Here we outline a set of experiments designed to explore AMOC hysteresis and sensitivity to additional freshwater input as part of the North Atlantic hosing model intercomparison project (NAHosMIP). These experiments include adding additional freshwater (hosing) for a fixed length of time to examine the rate and mechanisms or AMOC weakening, and whether the AMOC subsequently recovers once hosing stops. Initial results are shown from eight climate models participating in the Sixth Coupled Model Intercomparison Project CMIP6). The AMOC weakens in all models from the freshening, but once the freshening ceases, the AMOC recovers in half of the models, and in the other half it stays in a weakened state. The difference in model behaviour cannot be explained by the ocean model resolution or type, or by details of subgridscale parameterizations. Nor can it be explained by previously proposed properties of the mean climate state such as the strength of the salinity advection feedback. Instead the AMOC recovery is determined by the climate state reached when hosing stops, with those experiments where the AMOC is weakest not experiencing a recovery. [ABSTRACT FROM AUTHOR]

Mechanistic causes for sea level (SL) change patterns are analyzed as they emerge from the Coupled Model Intercomparison Project Phase 6 (CMIP6) endorsed Flux-Anomaly-Forced Model Intercomparison Project (FAFMIP) coupled climate experiments imposing individual forcing anomalies in wind stress, heatflux and freshwater flux to the Max-Planck-Institute Earth System Model (MPI-ESM). It appears that the heat flux perturbations have the largest effect on the sea level. In contrast, the direct impact of momentum and freshwater flux anomalies on SL anomalies appear to be limited to some region e.g. the Southern Ocean, Arctic Ocean and to some extent the North Pacific and North Atlantic Ocean. We find that thermosteric changes dominate the total SL change over large parts of the global ocean, except north of 60 °N where halosteric changes prevail. An analysis of added and redistributed components of heat and freshwater further suggests that the added component dominates the thermosteric SL and the redistributed component dominates the halosteric SL. Due to feedback processes a superposition of all forcing components together leads to the simulated sea level changes in each individual experiment. As a result, large surface heat flux anomalies over the Atlantic lead to wind stress change outside of the Atlantic through teleconnections, which in turn appear to be the primary driving agent for changes of sea level outside of the Atlantic in all three experiments. The associated wind driven Sverdrup stream function implicates that outside of the Atlantic most of the feedback can be explained by changes in the Sverdrup circulation. [ABSTRACT FROM AUTHOR]

An on-line, ensemble-based data assimilation (DA) method is performed to reconstruct the climate for 1750-1850 AD, and the performance is evaluated on large and small spatial scales. We use a low-resolution version of the Max Planck Institute for Meteorology MPI-ESM model and assimilate the PAGES 2K continental mean temperature reconstructions for the Northern Hemisphere (NH). The ensembles are generated sequentially for sub-periods based on the analysis of previous sub-periods. The assimilation has good skill for large-scale temperatures, but there is no agreement between the DA analysis and proxy-based reconstructions for small-scale temperature patterns within Europe or with reconstructions for the North Atlantic Oscillation (NAO) index. To explain the lack of added value in small spatial scales, a maximum covariance analysis (MCA) of links between NH temperature and sea level pressure is performed based on a control simulation with MPI-ESM. For annual values, winter and spring the Northern Annular Mode (NAM) is the pattern that is most closely linked to the NH continental temperatures, while for summer and autumn it is a wave-like pattern. This link is reproduced in the DA for winter, spring and annual means, providing potential for constraining the NAM/NAO phase and in turn regional temperature variability. It is shown that the lack of actual small-scale skill is likely due to the fact that the link might be too weak, as the NH continental mean temperatures are not the best predictors for large-scale circulation anomalies, or that the PAGES 2K temperatures include noise. Both factors can lead to circulation anomalies in the DA analysis that are substantially different from reality, leading to unrealistic representation of small-scale temperature variability. Moreover, we show that even if the true amplitudes of the leading MCA circulation patterns were known, there is still a large amount of unexplained local temperature variance. Based on these results, we argue that assimilating temperature reconstructions with a higher spatial resolution might improve the DA performance. [ABSTRACT FROM AUTHOR]

Responses of summer [June-August (JJA)] temperature and precipitation to large volcanic eruptions are analyzed using the millennial simulations of the earth system model developed at the Max Planck Institute for Meteorology. The model was driven by up-to-date reconstructions of external forcing, including natural forcing (solar and volcanic) and anthropogenic forcing (land-cover change and greenhouse gases). Cooling anomalies after large volcanic eruptions are seen on a nearly global scale. The cooling in the Northern Hemisphere (NH) is stronger than in the Southern Hemisphere (SH), and cooling is stronger over the continents than over the oceans. The precipitation decreases in the tropical and subtropical regions in the first summer after large volcanic eruptions. The cooling, with amplitudes of up to −0.6°C, is also seen over eastern China. East Asia is dominated by northerly wind anomalies, and the corresponding summer rainfall exhibits a coherent reduction over the entirety of eastern China. The tropospheric mean temperature anomalies indicate that there is coherent cooling over East Asia and the tropical ocean after large volcanic eruptions. The cooling over the middle-to-high latitudes of East Asia is stronger than over the tropical ocean. This temperature anomaly pattern suggests a reduced land-sea thermal contrast and favors a weaker East Asian summer monsoon (EASM) circulation. Analysis of the radiative fluxes at the top of the atmosphere (TOA) suggests that the reduction in shortwave radiation after large volcanic eruptions is nearly twice as large as the reduction in emitted longwave radiation, a net loss of radiative energy that cools the surface and lower troposphere. [ABSTRACT FROM AUTHOR]

In projections of twenty-first century climate, Arctic sea ice declines and at the same time exhibits strong interannual anomalies. Here, we investigate the potential to predict these strong sea-ice anomalies under a perfect-model assumption, using the Max-Planck-Institute Earth System Model in the same setup as in the Coupled Model Intercomparison Project Phase 5 (CMIP5). We study two cases of strong negative sea-ice anomalies: a 5-year-long anomaly for present-day conditions, and a 10-year-long anomaly for conditions projected for the middle of the twenty-first century. We treat these anomalies in the CMIP5 projections as the truth, and use exactly the same model configuration for predictions of this synthetic truth. We start ensemble predictions at different times during the anomalies, considering lagged-perfect and sea-ice-assimilated initial conditions. We find that the onset and amplitude of the interannual anomalies are not predictable. However, the further deepening of the anomaly can be predicted for typically 1 year lead time if predictions start after the onset but before the maximal amplitude of the anomaly. The magnitude of an extremely low summer sea-ice minimum is hard to predict: the skill of the prediction ensemble is not better than a damped-persistence forecast for lead times of more than a few months, and is not better than a climatology forecast for lead times of two or more years. Predictions of the present-day anomaly are more skillful than predictions of the mid-century anomaly. Predictions using sea-ice-assimilated initial conditions are competitive with those using lagged-perfect initial conditions for lead times of a year or less, but yield degraded skill for longer lead times. The results presented here suggest that there is limited prospect of predicting the large interannual sea-ice anomalies expected to occur throughout the twenty-first century. [ABSTRACT FROM AUTHOR]

The Atlantic multidecadal variability (AMV) modulates the North Atlantic surface ocean variability and affects decadal climates over the globe; its underlying mechanisms remain, however, under debate. In this study, we use a multi‐model ensemble of transient past‐millennium (850–1849) and unperturbed preindustrial control simulations contributing to the paleoclimate modeling intercomparison project—phase 4 (PMIP4) to decompose the AMV signal into the internal AMV and the external signal. The internal component of AMV exhibits no robust behavior across simulations during periods of major forcing such as strong volcanic eruptions, whereas the external forced temperature responds to volcanic eruptions with an immediate radiative cooling followed, in some simulations, by a sequence of damped multidecadal oscillations. The internal component tightly relates with the Atlantic meridional overturning circulation (AMOC) and dominates the fluctuations of AMV; whereas the external signal has limited impacts on AMOC and explains ∼25% of the AMV variance over the past millennium. Plain Language Summary: The decadal‐scale changes of the sea surface temperature over the North Atlantic is referred to as the Atlantic multidecadal variability (AMV), which modulates the North Atlantic surface ocean variability and affects decadal climates over the globe; however, how different forcing processes can affect AMV remains debated. In this study, a multi‐model ensemble participating in the paleoclimate modeling intercomparison project—phase 4 (PMIP4) of preindustrial millennium simulation (850–1849 CE) is used to decompose the AMV signal into the internal AMV and the external signal and to compare it with the preindustrial control simulation not including external forcing (such as volcanic eruptions, solar irradiance, and anthropogenic greenhouse gas). The internal component of AMV exhibits no robust signature across simulations during periods of major forcing such as strong volcanic eruptions, whereas the external forced temperature responds to volcanic eruptions with an immediate radiative cooling followed, in some simulations, by a sequence of damped multidecadal oscillations. The internal component of AMV tightly relates with the Atlantic subsurface ocean circulation called Atlantic meridional overturning circulation (AMOC) and dominates the variations of AMV, while the external signal has limited impacts on AMOC and explains ∼25% of the AMV variance over past millennium. Key Points: The internal and external contributions to Atlantic multidecadal variability (AMV) over the past millennium can be separated by a pattern‐based calculation of AMVThe internal AMV has a robust tripole pattern and dominates the relationship with meridional overturning circulationThe North Atlantic response to post‐eruption global cooling has a damped behavior and is identified as a more uniform spatial pattern [ABSTRACT FROM AUTHOR]

This paper investigates the impact of different ocean initialization strategies on the forecast skill of decadal prediction experiments performed with the ECHAM5/Max Planck Institute Ocean Model (MPI-OM) coupled model. The ocean initializations assimilate three-dimensional temperature and salinity anomalies from two different ocean state estimates, the ocean reanalysis of the German contribution to Estimating the Circulation and Climate of the Ocean (GECCO) and an ensemble of MPI-OM ocean experiments forced with the NCEP-NCAR atmospheric reanalysis. The results show that North Atlantic and Mediterranean sea surface temperature (SST) variations can be skillfully predicted up to a decade ahead and with greater skill than by both uninitialized simulations and persistence forecasts. The regional distribution of SST predictive skill is similar in both initialization approaches; however, higher skill is found for the NCEP hindcasts than for the GECCO hindcasts when a combination of predictive skill measures is used. Skillful predictions of surface air temperature are obtained over northwestern Europe, northern Africa, and central-eastern Asia. The North Atlantic subpolar gyre region stands out as the region with the highest predictive skill beyond the warming trend, in both SST and upper-ocean heat-content predictions. Here the NCEP hindcasts deliver the best results due to a more accurate initialization of the observed variability. The dominant mechanism for North Atlantic climate predictability is of dynamical origin and can be attributed to the initialization of the Atlantic meridional overturning circulation, thus explaining the reoccurrence of high predictive skill within the second pentad of the hindcasts experiments. The results herein demonstrate that ocean experiments forced with the observed history of the atmospheric state constitute a simple but successful alternative strategy for the initialization of skillful climate predictions over the next decade. [ABSTRACT FROM AUTHOR]

The current study deals with a parameterization of diapycnal diffusivity in an ocean model. The parameterization estimates the diapycnal diffusivity depending on the location of tidal-related energy dissipation over rough topography. The scheme requires a bottom roughness map that can be chosen depending on the scales of topographic features. Here, we implement the parameterization on an ocean general circulation model, and we examine the sensitivity of the modeled circulations to different spatial scales of the modeled bottom roughness. We compare three simulations that include the tidal mixing scheme using bottom roughness calculated at three different ranges of spatial scales, with the largest scale varying up to 200 km. Three main results are discussed. First, the dependence of the topographic spectra with depth, characterized by an increase in spectral energy over short length scales in the deep ocean, influences the vertical profile of the diffusivity. Second, the changes in diffusivities lead to different equilibrium solutions in the Atlantic meridional overturning circulation and bottom circulation. In particular, the lower cell of the Atlantic overturning and the bottom water transport in the Pacific Ocean are stronger for stronger diffusivities at the corresponding basins and depths, and the strongest when using the small-scale roughness map. Third, a comparison of the density fields of the three simulations with the density field of World Ocean Atlas dataset, from which the models are initialized, shows that among the simulations with three different roughness maps, the one using small-scale bottom roughness map has the smallest density bias. [ABSTRACT FROM AUTHOR]

The decadal-centennial variations of East Asian summer monsoon (EASM) and the associated rainfall change during the past millennium are simulated using the earth system model developed at the Max Planck Institute for Meteorology. The model was driven by up-to-date reconstructions of external forcing including the recent low-amplitude estimates of solar variations. Analysis of the simulations indicates that the EASM is generally strong during the Medieval Warm Period (MWP; A.D. 1000-1100) and weak during the Little Ice Age (LIA; A.D. 1600-1700). The monsoon rainband exhibits a meridional tripolar pattern during both epochs. Excessive (deficient) precipitation is found over northern China (35°-42°N, 100°-120°E) but deficient (excessive) precipitation is seen along the Yangtze River valley (27°-34°N, 100°-120°E) during the MWP (LIA). Both similarities and disparities of the rainfall pattern between the model results herein and the proxy data have been compared, and reconstructions from Chinese historical documents and some geological evidence support the results. The changes of the EASM circulation including the subtropical westerly jet stream in the upper troposphere and the western Pacific subtropical high (WPSH) in the middle and lower troposphere are consistent with the meridional shift of the monsoon rain belt during both epochs. The meridional monsoon circulation changes are accompanied with anomalous southerly (northerly) winds between 20° and 50°N during the MWP (LIA). The land-sea thermal contrast change caused by the effective radiative forcing leads to the MWP and LIA monsoon changes. The 'warmer land-colder ocean' anomaly pattern during the MWP favors a stronger monsoon, while the 'colder land-warmer ocean' anomaly pattern during the LIA favors a weaker monsoon. [ABSTRACT FROM AUTHOR]

Attempts to predict changes in Atlantic Meridional Overturning Circulation (AMOC) have yielded little success to date. Here, we demonstrate predictability for monthly mean AMOC strength at 26.5°N for up to 4 years in advance. This AMOC predictive skill arises predominantly from the basin-wide upper-mid-ocean geostrophic transport, which in turn can be predicted because we have skill in predicting the upper-ocean zonal density difference. Ensemble forecasts initialized between January 2008 and January 2011 indicate a stable AMOC at 26.5°N unfit at least 2014, despite a brief wind-induced weakening in 2010. Because AMOC influences many aspects of climate, our results establish AMOC as an important potential carrier of climate predictability. [ABSTRACT FROM AUTHOR]

Abstract: This study explores the effects of the Changjiang (also called the Yangtze River) river discharge (CRD) on the density stratifications and associated sea surface temperature (SST) changes using a global ocean general circulation model with regional focus on the Yellow and East China Seas (YECS). It is found that CRD increases the SST in summer through a barrier layer (BL) formation that tends to enhance stratification at the mixed layer base, and thus reduces both vertical mixing and entrainment. This process is effective, particularly in August, after the CRD reaches its maximum in July. The SST difference between the composites of flood and drought years confirms that the surface warming is related to surface freshening by the CRD. This result suggests that the BL induced by the CRD is an important contributor to the surface heat budget in the YECS. [Copyright &y& Elsevier]

The sensitivity of the climate-biogeochemistry system to volcanic eruptions is investigated using the comprehensive Earth System Model developed at the Max Planck Institute for Meteorology. The model includes an interactive carbon cycle with modules for terrestrial biosphere as well as ocean biogeochemistry. The volcanic forcing is based on a recent reconstruction for the last 1200 yr. An ensemble of five simulations is performed and the averaged response of the system is analysed in particular for the largest eruption of the last millennium in the year 1258. After this eruption, the global annual mean temperature drops by 1 K and recovers slowly during 10 yr. Atmospheric CO concentration declines during 4 yr after the eruption by ca. 2 ppmv to its minimum value and then starts to increase towards the pre-eruption level. This CO decrease is explained mainly by reduced heterotrophic respiration on land in response to the surface cooling, which leads to increased carbon storage in soils, mostly in tropical and subtropical regions. The ocean acts as a weak carbon sink, which is primarily due to temperature-induced solubility. This sink saturates 2 yr after the eruption, earlier than the land uptake. [ABSTRACT FROM AUTHOR]

Sea ice variability in the Barents Sea and its impact on climate are analyzed using a 465-year control integration of a global coupled atmosphere–ocean–sea ice model. Sensitivity simulations are performed to investigate the response to an isolated sea ice anomaly in the Barents Sea. The interannual variability of sea ice volume in the Barents Sea is mainly determined by variations in sea ice import into Barents Sea from the Central Arctic. This import is primarily driven by the local wind field. Horizontal oceanic heat transport into the Barents Sea is of minor importance for interannual sea ice variations but is important on longer time scales. Events with strong positive sea ice anomalies in the Barents Sea are due to accumulation of sea ice by enhanced sea ice imports and related NAO-like pressure conditions in the years before the event. Sea ice volume and concentration stay above normal in the Barents Sea for about 2 years after an event. This strongly increases the albedo and reduces the ocean heat release to the atmosphere. Consequently, air temperature is much colder than usual in the Barents Sea and surrounding areas. Precipitation is decreased and sea level pressure in the Barents Sea is anomalously high. The large-scale atmospheric response is limited with the main impact being a reduced pressure over Scandinavia in the year after a large ice volume occurs in the Barents Sea. Furthermore, high sea ice volume in the Barents Sea leads to increased sea ice melting and hence reduced surface salinity. Generally, the climate response is smallest in summer and largest in winter and spring. [ABSTRACT FROM AUTHOR]

The ocean meridional overturning circulation (MOC) plays a central role for the climate in the Atlantic realm. Since scenarios for future climate change indicate a significant reduction of the MOC under global warming, an assessment of variations and trends of the real MOC is important. Using observations at ocean weather ship (OWS) stations and along oceanic sections, we examine the hydrographic information that can be used to determine MOC trends via its signature in water mass properties obtained from model simulations with the climate model ECHAM5/MPI-OM. We show that temperature trends at mid-latitudes provide useful indirect measure of large-scale changes of deep circulation: A mid-depth warming is related to MOC weakening and a cooling to MOC strengthening. Based on our model experiments, we argue that a continuation of measurements at key OWS sites may contribute to a timely detection of a possible future MOC slowdown and to separate the signal from interannual-to-multidecadal MOC variability. The simulations suggest that the subsurface hydrographic information related to MOC has a lower variability than the MOC trend measured directly. Based on our model and the available long-term hydrographic data, we estimate non-significant MOC trends for the last 80 years. For the twenty-first century, however, the model simulations predict a significant MOC decline and accompanied mid-depth warming trend. [ABSTRACT FROM AUTHOR]

For the first time, we compare the effects of four different ocean vertical mixing schemes on the mean state of the ocean and atmosphere in the Max Planck Institute Earth System Model (MPI-ESM1.2). These four schemes are namely the default (PP) scheme, the K -profile parameterization (KPP) from the Community Vertical Mixing (CVMix) library, a recently implemented scheme based on turbulent kinetic energy (TKE), and a recently developed prognostic scheme for internal wave dissipation, energy, and mixing (IDEMIX) to replace the often assumed constant background diffusivity in the ocean interior. In this study, the IDEMIX scheme is combined with the TKE scheme (collectively called the TKE + IDEMIX scheme) to provide an energetically more consistent framework for mixing, as it does not rely on the unwanted effect of creating spurious energy for mixing. Energetic consistency can have implications on the climate. Therefore, we focus on the effects of TKE + IDEMIX on the climate mean state and compare them with the first three schemes that are commonly used in other models but are not energetically consistent. We find warmer sea surface temperatures (SSTs) in the North Atlantic and Nordic Seas using KPP or TKE(+ IDEMIX), which is related to 10 % higher overflows that cause a stronger and deeper upper cell of the Atlantic meridional overturning circulation (AMOC) and thereby an enhanced northward heat transport and higher inflow of warm and saline water from the Indian Ocean into the South Atlantic. Saltier subpolar North Atlantic and Nordic Seas lead to increased deep convection and thus to the increased overflows. Due to the warmer SSTs, the extratropics of the Northern Hemisphere become warmer with TKE(+ IDEMIX), weakening the meridional gradient and thus the jet stream. With KPP, the tropics and the Southern Hemisphere also become warmer without weakening the jet stream. Using an energetically more consistent scheme (TKE + IDEMIX) produces a more heterogeneous and realistic pattern of vertical eddy diffusivity, with lower diffusivities in deep and flat-bottom basins and elevated turbulence over rough topography. IDEMIX improves in particular the diffusivity in the Arctic Ocean and reduces the warm bias in the Atlantic Water layer. We conclude that although shortcomings due to model resolution determine the global-scale bias pattern, the choice of the vertical mixing scheme may play an important role for regional biases. [ABSTRACT FROM AUTHOR]

Land-atmosphere interactions are often interpreted as local effects, whereby the soil state drives local atmospheric conditions and feedbacks originate. However, nonlocal mechanisms can significantly modulate land-atmosphere exchanges and coupling. We make use of GCMs at different resolutions (low ~1° and high ~0.25°) to separate the two contributions to coupling: better represented local processes versus the influence of improved large-scale circulation. We use a two-legged metric, complemented by a process-based assessment of four CMIP6 GCMs. Our results show that weakening, strengthening, and relocation of coupling hot spots occur at high resolution globally. The northward expansion of the Sahel hot spot, driven by nonlocal mechanisms, is the most notable change. The African easterly jet's horizontal wind shear is enhanced in JJA due to better resolved orography at high resolution. This effect, combined with enhanced easterly moisture flux, favors the development of African easterly waves over the Sahel. More precipitation and soil moisture recharge produce strengthening of the coupling, where evapotranspiration remains controlled by soil moisture, and weakening where evapotranspiration depends on atmospheric demand. In SON, the atmospheric influence is weaker, but soil memory helps to maintain the coupling between soil moisture and evapotranspiration and the relocation of the hot spot at high resolution. The multimodel agreement provides robust evidence that atmospheric dynamics determines the onset of land-atmosphere interactions, while the soil state modulates their duration. Comparison of precipitation, soil moisture, and evapotranspiration against satellite data reveals that the enhanced moistening at high resolution significantly reduces model biases, supporting the realism of the hot-spot relocation. [ABSTRACT FROM AUTHOR]

Sea levels of different atmosphere–ocean general circulation models (AOGCMs) respond to climate change forcing in different ways, representing a crucial uncertainty in climate change research. We isolate the role of the ocean dynamics in setting the spatial pattern of dynamic sea-level (ζ) change by forcing several AOGCMs with prescribed identical heat, momentum (wind) and freshwater flux perturbations. This method produces a ζ projection spread comparable in magnitude to the spread that results from greenhouse gas forcing, indicating that the differences in ocean model formulation are the cause, rather than diversity in surface flux change. The heat flux change drives most of the global pattern of ζ change, while the momentum and water flux changes cause locally confined features. North Atlantic heat uptake causes large temperature and salinity driven density changes, altering local ocean transport and ζ. The spread between AOGCMs here is caused largely by differences in their regional transport adjustment, which redistributes heat that was already in the ocean prior to perturbation. The geographic details of the ζ change in the North Atlantic are diverse across models, but the underlying dynamic change is similar. In contrast, the heat absorbed by the Southern Ocean does not strongly alter the vertically coherent circulation. The Arctic ζ change is dissimilar across models, owing to differences in passive heat uptake and circulation change. Only the Arctic is strongly affected by nonlinear interactions between the three air-sea flux changes, and these are model specific. [ABSTRACT FROM AUTHOR]