Search

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]

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]

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]

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]

In the Northern Hemisphere, the mid-6th century was one of the coldest periods of the last 2000 years, as indicated by both proxy records and Earth System Model (ESM) simulations. This cold period was initiated by volcanic eruptions in 536 CE and 540 CE. Evidence from historical sources, archaeological findings and proxy records suggests that the extent and severity of this volcanic induced cooling was spatially heterogeneous, and that the effect on society resulted in adaptation and resilience at some locations, whereas social crisis has been indicated at others. Here, we study the effect of the volcanic double event in 536 CE and 540 CE on the climate and society in Scandinavia with a special focus on Southern Norway. Using an ensemble of Max Planck Institute ESM transient simulations for 521-680 CE, the temperature, precipitation and atmospheric circulation patterns are studied. The simulated cooling magnitude is then used as input for the growing degree day (GDD) model set-up for Southern Norway. This GDD model indicates the possible effects on agriculture for three different study areas in Southern Norway, representative of typical meteorological and landscape conditions. Pollen from bogs and archaeological records inside the study area are then analysed at high resolution (1-3 cm sample intervals) to give insights into the validity of the GDD model set-up with regard to the volcanic climate impact on the regional scale, and to link the different types of data sets. After the 536/540 CE double event, a maximum surface air cooling of up to 3.5 °C during the mean growing season is simulated regionally in Southern Norway. With a worst-case scenario cooling of 3 °C, the GDD model indicates crop failures were likely in our northernmost and western study areas, while crops were more likely to mature in the southeastern study area. These results are in agreement with the pollen records from the respective areas. During the sixth century, excavations show an abandonment of farms, severe social impact but also a continuation of occupation or a mix of those. In addition, archaeological findings from one of the excavation sites suggest wetter conditions for the mid sixth century in Scandinavia, as simulated by individual ensemble members. Finally, we discuss the likely climate and societal impacts of the 536/540 CE volcanic double event by synthesising the new and available data sets for the whole Scandinavia. [ABSTRACT FROM AUTHOR]

Numerical modelling enables a comprehensive understanding not only of the Earth's system today, but also of the past. To date, a significant amount of time and effort has been devoted to paleoclimate modeling and analysis, which involves the latest and most advanced Paleoclimate Modelling Intercomparison Project phase 4 (PMIP4). The definition of seasonality, which is influenced by slow variations in the Earth's orbital parameters, plays a key role in determining the calculated seasonal cycle of the climate. In contrast to the classical calendar used today, where the lengths of the months and seasons are fixed, the angular calendar calculates the lengths of the months and seasons according to a fixed number of degrees along the Earth's orbit. When comparing simulation results for different time intervals, it is essential to account for the angular calendar to ensure that the data for comparison is from the same position along the Earth's orbit. Most models use the classical "fixed-length" calendar, which can lead to strong distortions of the monthly and seasonal values, especially for the climate of the past. Here, by analyzing daily outputs from multiple PMIP4 model simulations, we examine calendar effects on surface air temperature and precipitation under mid-Holocene, last interglacial, and pre-industrial climate conditions. We conclude that: (a) The largest cooling bias occurs in autumn when the classical calendar is applied for the mid-Holocene and last interglacial. (b) The sign of the temperature anomalies between the Last Interglacial and pre-industrial in boreal autumn can be reversed after the switch from classical to angular calendar, particularly over the Northern Hemisphere continents. (c) Precipitation over West Africa is overestimated in boreal summer and underestimated in boreal autumn when the "fixed-length" seasonal cycle is applied. (d) Finally, correcting the calendar based on the monthly model results can reduce the biases to a large extent, but not completely eliminate them. In addition, we examine the calendar effects in 3 transient simulations for 6-0 ka by AWI-ESM, MPI-ESM, and IPSL. We find significant discrepancies between adjusted and unadjusted temperature values over ice-free continents for both hemispheres in boreal autumn. While for other seasons the deviations are relatively small. A drying bias can be found in the summer monsoon precipitation in Africa (in the "fixed-length" calendar), whereby the magnitude of bias becomes smaller over time. Overall, our study underlines the importance of the application of calendar transformation in the analysis of climate simulations. Neglecting the calendar effects could lead to a profound artificial distortion of the calculated seasonal cycle of surface air temperature and precipitation. One important fact to be noted here is that the discrepancy in seasonality under different calendars is an analysis bias and is highly depends on the choice of the reference position/date (usually the vernal equinox, which is set to 31th March) on the Earth's ellipse around the sun. Different model groups may apply different reference dates, so ensuring a consistent reference date and seasonal definition is key when we compare results across multiple models. [ABSTRACT FROM AUTHOR]

The climate in the Northern Hemisphere (NH) of the mid-6th century was one of the coldest during the last two millennia. The onset of this cold period is attributed to the volcanic double eruption event in 536 and 540 Common Era (CE) based on multiple paleo-proxies. Recently, there has been a debate about how long lasting and cold this volcanic induced cold period actually was. To better understand this, we analyze new transient simulations over the Common Era and enhance the representation of mid 6th to 7th century climate by additional ensemble simulations covering 520-680 CE. We use the Max Planck Institute Earth System Model and apply external forcing as recommended in the Paleo Model Intercomparison Project, Phase 4. After the four large eruptions in 536, 540, 574, and 626 CE, a significant surface climate response up to 20 years is simulated. The Northern Hemisphere 2 m air temperature, and precipitation decreases up to 2 K, and 0.2 mm day-1, respectively, and sea ice area increases up to 1.5 ·1012 m². The global ocean heat content decreases drastically by 1.5·1023 Jm-1,which is significant for 30-40 years, and does not totally recover during the entire study period. The surface maps reveal atmospheric circulation changes with a hemispheric dipole pattern and land see contrast in the first two years after the eruptions. Poleward of ~45°N higher sea level pressure and a decrease in hydrological variables occur, accompanied by a land see contrast, with decreased values over land and an increase in values over the ocean, which is especially pronounced for evaporation during boreal summer. During boreal winter, a positive North Atlantic Oscillation develops in the first year after (three out of the four) large eruptions. Analysing underlying mechanisms in the North Atlantic reveals that complex interaction between sea-ice expansion, changes in barotropic streamfunction and meridional overturning circulation leads to a reduction in the ocean heat transport, which then further enhances sea ice expansion impacting NH surface climate up to 20 years. Temperature records reconstructed from tree-rings in the NH agree well with the model simulations and show a similar ~ 20 year cooling after the eruptions. A century of surface cooling starting in the mid-6th century, as shown from local tree-ring records from the Alps and Altai, does not occur in our volcanic climate model simulations, nor in the NH tree-ring compilation. [ABSTRACT FROM AUTHOR]