Search

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 5 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 subsurface that is too shallow, 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, which are 12 %–42 % smaller than the latest borehole-based estimates (18.2 ZJ).

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.

In the Northern Hemisphere, the mid-6th century was one of the coldest periods of the last 2000 years, which was initiated by volcanic eruptions in 536 and 540 CE. Here, we study the effect of this volcanic double event on the climate and society in Scandinavia with a special focus on southern Norway. Using an ensemble of Max Planck Institute Earth system model transient simulations for 521–680 CE, temperature, precipitation, and atmospheric circulation patterns are analyzed. The simulated cooling magnitude is used as input for a growing degree day (GDD) model setup for three different study areas in southern Norway, representative of typical meteorological and landscape conditions. Pollen from bogs inside these study areas are analyzed at high resolution (1–3 cm sample intervals) to give insights into the validity of the GDD model setup with regard to the volcanic climate impact on the regional scale and to link the different data sets with the archeological records. We find that after the 536 and 540 CE double event, a maximum surface air cooling of up to 3.5 ∘C during the mean growing season is simulated regionally for southern Norway. With a 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. Archeological excavations show, however, a more complex pattern for the three areas with abandonment of farms and severe social impacts but also a continuation of occupation or a mix of those. Finally, we discuss the likely climatic and societal impacts of the 536 and 540 CE volcanic double event by synthesizing the new and available data sets for the whole Scandinavia.

State-of-the-art Earth system models typically employ grid spacings of O(100 km), which is too coarse to explicitly resolve main drivers of the flow of energy and matter across the Earth system. In this paper, we present the new ICON-Sapphire model configuration, which targets a representation of the components of the Earth system and their interactions with a grid spacing of 10 km and finer. Through the use of selected simulation examples, we demonstrate that ICON-Sapphire can (i) be run coupled globally on seasonal timescales with a grid spacing of 5 km, on monthly timescales with a grid spacing of 2.5 km, and on daily timescales with a grid spacing of 1.25 km; (ii) resolve large eddies in the atmosphere using hectometer grid spacings on limited-area domains in atmosphere-only simulations; (iii) resolve submesoscale ocean eddies by using a global uniform grid of 1.25 km or a telescoping grid with the finest grid spacing at 530 m, the latter coupled to a uniform atmosphere; and (iv) simulate biogeochemistry in an ocean-only simulation integrated for 4 years at 10 km. Comparison of basic features of the climate system to observations reveals no obvious pitfalls, even though some observed aspects remain difficult to capture. The throughput of the coupled 5 km global simulation is 126 simulated days per day employing 21 % of the latest machine of the German Climate Computing Center. Extrapolating from these results, multi-decadal global simulations including interactive carbon are now possible, and short global simulations resolving large eddies in the atmosphere and submesoscale eddies in the ocean are within reach.

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.

The climate of the Northern Hemisphere (NH) in the mid-6th century was one of the coldest during the last 2 millennia based on multiple paleo-proxies. While the onset of this cold period can be clearly connected to the volcanic eruptions in 536 and 540 Common Era (CE), the duration, extent, and magnitude of the cold period are uncertain. Proxy data are sparse for the first millennium, which compounds the uncertainties of the reconstructions. To better understand the mechanisms of the prolonged cooling, 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 to apply the external forcing as recommended in the Paleoclimate Modelling Intercomparison Project phase 4. After the four large eruptions in 536, 540, 574, and 626 CE, a significant mean surface climate response in the NH lasting up to 20 years is simulated. The 2 m air temperature shows a cooling over the Arctic in winter, corresponding to the increase in Arctic sea ice, mainly in the Labrador Sea and to the east of Greenland. The increase in sea-ice extent relates to a decrease in the northward ocean heat transport into the Arctic within the first 2 years after the eruptions and to an increase in the Atlantic meridional overturning circulation, which peaks 10 years after the eruptions. A decrease in the global ocean heat content is simulated after the eruptions that does not recover during the simulation period. These ocean–sea-ice interactions sustain the surface cooling, as the cooling lasts longer than is expected solely from the direct effects of the volcanic forcing, and are thus responsible for the multi-decadal surface cooling. In boreal summer, the main cooling occurs over the continents at midlatitudes. A dipole pattern develops with high sea level pressure and a decrease in both precipitation and evaporation poleward of 40∘ N. In addition, more pronounced cooling over land compared to ocean leads to an enhanced land–sea contrast. While our model ensemble simulations show a similar ∼20-year summer cooling over NH land after the eruptions as tree ring reconstructions, a volcanic-induced century-long cooling, as reconstructed from tree ring data, does not occur in our simulations.

Numerical modeling 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 are from the same position along the Earth's orbit. Most models use the classical 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 came to the following conclusions. (a) The largest cooling bias occurs in boreal autumn when the classical calendar is applied for the mid-Holocene and Last Interglacial, due to the fact that the vernal equinox is fixed on 21 March. (b) The sign of the temperature anomalies between the Last Interglacial and pre-industrial in boreal autumn can be reversed after the switch from the 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 classical seasonal cycle is applied. (d) Finally, month-length adjusted values for surface air temperature and precipitation are very similar to the day-length adjusted values, and therefore correcting the calendar based on the monthly model results can largely reduce the artificial bias. In addition, we examine the calendar effects in three transient simulations for 6–0 ka by AWI-ESM, MPI-ESM, and IPSL-CM. We find significant discrepancies between adjusted and unadjusted temperature values over 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 classical 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.

Human land use activities have resulted in large changes to the biogeochemical and biophysical properties of the Earth's surface, with consequences for climate and other ecosystem services. In the future, land use activities are likely to expand and/or intensify further to meet growing demands for food, fiber, and energy. As part of the World Climate Research Program Coupled Model Intercomparison Project (CMIP6), the international community has developed the next generation of advanced Earth system models (ESMs) to estimate the combined effects of human activities (e.g., land use and fossil fuel emissions) on the carbon–climate system. A new set of historical data based on the History of the Global Environment database (HYDE), and multiple alternative scenarios of the future (2015–2100) from Integrated Assessment Model (IAM) teams, is required as input for these models. With most ESM simulations for CMIP6 now completed, it is important to document the land use patterns used by those simulations. Here we present results from the Land-Use Harmonization 2 (LUH2) project, which smoothly connects updated historical reconstructions of land use with eight new future projections in the format required for ESMs. The harmonization strategy estimates the fractional land use patterns, underlying land use transitions, key agricultural management information, and resulting secondary lands annually, while minimizing the differences between the end of the historical reconstruction and IAM initial conditions and preserving changes depicted by the IAMs in the future. The new approach builds on a similar effort from CMIP5 and is now provided at higher resolution (0.25∘×0.25∘) over a longer time domain (850–2100, with extensions to 2300) with more detail (including multiple crop and pasture types and associated management practices) using more input datasets (including Landsat remote sensing data) and updated algorithms (wood harvest and shifting cultivation); it is assessed via a new diagnostic package. The new LUH2 products contain > 50 times the information content of the datasets used in CMIP5 and are designed to enable new and improved estimates of the combined effects of land use on the global carbon–climate system.

The Flux-Anomaly-Forced Model Intercomparison Project (FAFMIP) aims to investigate the spread in simulations of sea-level and ocean climate change in response to CO2 forcing by atmosphere–ocean general circulation models (AOGCMs). It is particularly motivated by the uncertainties in projections of ocean heat uptake, global-mean sea-level rise due to thermal expansion and the geographical patterns of sea-level change due to ocean density and circulation change. FAFMIP has three tier-1 experiments, in which prescribed surface flux perturbations of momentum, heat and freshwater respectively are applied to the ocean in separate AOGCM simulations. All other conditions are as in the pre-industrial control. The prescribed fields are typical of pattern and magnitude of changes in these fluxes projected by AOGCMs for doubled CO2 concentration. Five groups have tested the experimental design with existing AOGCMs. Their results show diversity in the pattern and magnitude of changes, with some common qualitative features. Heat and water flux perturbation cause the dipole in sea-level change in the North Atlantic, while momentum and heat flux perturbation cause the gradient across the Antarctic Circumpolar Current. The Atlantic meridional overturning circulation (AMOC) declines in response to the heat flux perturbation, and there is a strong positive feedback on this effect due to the consequent cooling of sea-surface temperature in the North Atlantic, which enhances the local heat input to the ocean. The momentum and water flux perturbations do not substantially affect the AMOC. Heat is taken up largely as a passive tracer in the Southern Ocean, which is the region of greatest heat input, while the weakening of the AMOC causes redistribution of heat towards lower latitudes. Future analysis of these and other phenomena with the wider range of CMIP6 FAFMIP AOGCMs will benefit from new diagnostics of temperature and salinity tendencies, which will enable investigation of the model spread in behaviour in terms of physical processes as formulated in the models.

Data assimilation (DA) is an emerging topic in palaeoclimatology and one of the key challenges in this field. Assimilating proxy-based continental mean temperature reconstructions into the MPI-ESM model showed a lack of information propagation to small spatial scales . Here, we investigate whether this lack of regional skill is due to the methodology or to errors in the assimilated reconstructions. Error separation is fundamental, as it can lead to improvements in DA methods. We address the question by performing a new set of simulations, using two different sets of target data; the proxy-based PAGES 2K reconstructions (DA-P scheme), and the HadCRUT3v instrumental observations (DA-I scheme). Again, we employ ensemble-member selection DA using the MPI-ESM model, and assimilate Northern Hemisphere (NH) continental mean temperatures; the simulated period is 1850–1949 AD. Both DA schemes follow the large-scale target and observed climate variations well, but the assimilation of instrumental data improves the performance. This improvement cannot be seen for Asia, where the limited instrumental coverage leads to errors in the target data and low skill for the DA-I scheme. No skill on small spatial scales is found for either of the two DA schemes, demonstrating that errors in the assimilated data are not the main reason for the unrealistic representation of the regional temperature variability in Europe and the NH. It can thus be concluded that assimilating continental mean temperatures is not ideal for providing skill on small spatial scales.

Different ensemble-based data assimilation (DA) approaches for palaeoclimate reconstructions have been recently undertaken, but no systematic comparison among them has been attempted. We compare an off-line and an on-line ensemble-based method, with the testing period being the 17th century, which led into the Maunder Minimum. We use a low-resolution version of Max Planck Institute for Meteorology Earth System Model (MPI-ESM) to assimilate the Past Global Changes (PAGES) 2k continental temperature reconstructions. In the off-line approach, the ensemble for the entire simulation period is generated first and then the ensemble is used in combination with the empirical information to produce the analysis. In contrast, in the on-line approach, the ensembles are generated sequentially for sub-periods based on the analysis of previous sub-periods. Both schemes perform better than the simulations without DA. The on-line method would be expected to perform better if the assimilation led to states of the slow components of the climate system that are close to reality and the system had sufficient memory to propagate this information forward in time. In our comparison, which is based on analysing correlations and differences between the analysis and the proxy-based reconstructions, we find similar skill for both methods on the continental and hemispheric scales. This indicates either a lack of control of the slow components in our setup or a lack of skill in the information propagation on decadal timescales. Additional experiments are however needed to check whether the conclusions reached in this particular setup are valid in other cases. Although the performance of the two schemes is similar and the on-line method is more difficult to implement, the temporal consistency of the analysis in the on-line method makes it in general preferable.

The response of the global climate-carbon cycle system to an extremely large Northern Hemisphere mid-latitude volcanic eruption is investigated using ensemble integrations with the comprehensive Earth System Model MPI-ESM. The model includes dynamical compartments of the atmosphere and ocean and interactive modules of the terrestrial biosphere as well as ocean biogeochemistry. The MPI-ESM was forced with anomalies of aerosol optical depth and effective radius of aerosol particles corresponding to a super eruption of the Yellowstone volcanic system. The model experiment consists of an ensemble of fifteen model integrations that are started at different pre-ENSO states of a control experiment and run for 200 years after the volcanic eruption. The climate response to the volcanic eruption is a maximum global monthly mean surface air temperature cooling of 3.8 K for the ensemble mean and from 3.3 K to 4.3 K for individual ensemble members. Atmospheric pCO2 decreases by a maximum of 5 ppm for the ensemble mean and by 3 ppm to 7 ppm for individual ensemble members approximately 6 years after the eruption. The atmospheric carbon content only very slowly returns to near pre-eruption level at year 200 after the eruption. The ocean takes up carbon shortly after the eruption in response to the cooling, changed wind fields and ice cover. This physics-driven uptake is weakly counteracted by a reduction of the biological export production mainly in the tropical Pacific. The land vegetation pool shows a decrease by 4 GtC due to reduced short-wave radiation that has not been present in a smaller scale eruption. The gain of the soil carbon pool determines the amplitude of the CO2 perturbation and the long-term behaviour of the overall system: an initial gain caused by reduced soil respiration is followed by a rather slow return towards pre-eruption levels. During this phase, the ocean compensates partly for the reduced atmospheric carbon content in response to the land's gain. In summary, we find that the volcanic eruption has long-lasting effects on the carbon cycle: After 200 years, the ocean and the land carbon pools are still different from the pre-eruption state by 3 GtC and 4 GtC, respectively, and the land carbon pools (vegetation and soil) show some long-lasting local anomalies that are only partly visible in the global signal.

The combination of inorganic semiconductors with organic thin films promises new strategies for the realization of complex hybrid optoelectronic devices. Oxidative chemical vapor deposition (oCVD) of conductive polymers offers a flexible and scalable path towards high-quality three-dimensional inorganic/organic optoelectronic structures. Here, hole-conductive poly(3,4-ethylenedioxythiophene) (PEDOT) grown by oxidative chemical vapor deposition is used to fabricate transparent and conformal wrap-around p-type contacts on three-dimensional microLEDs with large aspect ratios, a yet unsolved challenge in three-dimensional gallium nitride technology. The electrical characteristics of two-dimensional reference structures confirm the quasi-metallic state of the polymer, show high rectification ratios, and exhibit excellent thermal and temporal stability. We analyze the electroluminescence from a three-dimensional hybrid microrod/polymer LED array and demonstrate its improved optical properties compared with a purely inorganic microrod LED. The findings highlight a way towards the fabrication of hybrid three-dimensional optoelectronics on the sub-micron scale.

Reconstructions of the global mean annual temperature evolution during the Holocene yield conflicting results. One temperature reconstruction shows global cooling during the late Holocene. The other reconstruction reveals global warming. Here we show that both a global warming mode and a cooling mode emerge when performing a spatio-temporal analysis of annual temperature variability during the Holocene using data from a transient climate model simulation. The warming mode is most pronounced in the tropics. The simulated cooling mode is determined by changes in the seasonal cycle of Arctic sea-ice that are forced by orbital variations and volcanic eruptions. The warming mode dominates in the mid-Holocene, whereas the cooling mode takes over in the late Holocene. The weighted sum of the two modes yields the simulated global temperature trend evolution. Our findings have strong implications for the interpretation of proxy data and the selection of proxy locations to compute global mean temperatures.

A new release of the Max Planck Institute for Meteorology Earth System Model version 1.2 (MPI-ESM1.2) is presented. The development focused on correcting errors in and improving the physical processes representation, as well as improving the computational performance, versatility, and overall user friendliness. In addition to new radiation and aerosol parameterizations of the atmosphere, several relatively large, but partly compensating, coding errors in the model's cloud, convection, and turbulence parameterizations were corrected. The representation of land processes was refined by introducing a multilayer soil hydrology scheme, extending the land biogeochemistry to include the nitrogen cycle, replacing the soil and litter decomposition model and improving the representation of wildfires. The ocean biogeochemistry now represents cyanobacteria prognostically in order to capture the response of nitrogen fixation to changing climate conditions and further includes improved detritus settling and numerous other refinements. As something new, in addition to limiting drift and minimizing certain biases, the instrumental record warming was explicitly taken into account during the tuning process. To this end, a very high climate sensitivity of around 7 K caused by low-level clouds in the tropics as found in an intermediate model version was addressed, as it was not deemed possible to match observed warming otherwise. As a result, the model has a climate sensitivity to a doubling of CO 2 over preindustrial conditions of 2.77 K, maintaining the previously identified highly nonlinear global mean response to increasing CO 2 forcing, which nonetheless can be represented by a simple two-layer model., (©2019. The Authors.)

There is growing evidence that singing can have a positive effect on language learning, but few studies have explored its benefit for children who have recently migrated to a new country. In the present study, recently migrated children ( N = 35) received three 40-min sessions where all students learnt the lyrics of two songs designed to simulate language learning through alternating teaching modalities (singing and speaking). Children improved their language knowledge significantly including on tasks targeting the transfer of grammatical skills, an area largely neglected in previous studies. This improvement was sustainable over the retention interval. However, the two teaching modalities did not show differential effects on cued recall of song lyrics indicating that singing and speaking are equally effective when used in combination with one another. Taken together, the data suggest that singing may be useful as an additional teaching strategy, irrespective of initial language proficiency, warranting more research on songs as a supplement for grammar instruction.

The climate of the North Atlantic region exhibits fluctuations on decadal timescales that have large societal consequences. Prominent examples include hurricane activity in the Atlantic, and surface-temperature and rainfall variations over North America, Europe and northern Africa. Although these multidecadal variations are potentially predictable if the current state of the ocean is known, the lack of subsurface ocean observations that constrain this state has been a limiting factor for realizing the full skill potential of such predictions. Here we apply a simple approach-that uses only sea surface temperature (SST) observations-to partly overcome this difficulty and perform retrospective decadal predictions with a climate model. Skill is improved significantly relative to predictions made with incomplete knowledge of the ocean state, particularly in the North Atlantic and tropical Pacific oceans. Thus these results point towards the possibility of routine decadal climate predictions. Using this method, and by considering both internal natural climate variations and projected future anthropogenic forcing, we make the following forecast: over the next decade, the current Atlantic meridional overturning circulation will weaken to its long-term mean; moreover, North Atlantic SST and European and North American surface temperatures will cool slightly, whereas tropical Pacific SST will remain almost unchanged. Our results suggest that global surface temperature may not increase over the next decade, as natural climate variations in the North Atlantic and tropical Pacific temporarily offset the projected anthropogenic warming.