Description: We investigate the respective role of variations in subpolar deep water formation and Nordic Seas overflows for the decadal to multidecadal variability of the Atlantic meridional overturning circulation (AMOC). This is partly done by analysing long (order of 1000 years) control simulations with five coupled climate models. For all models, the maximum influence of variations in subpolar deep water formation is found at about 45° N, while the maximum influence of variations in Nordic Seas overflows is rather found at 55 to 60° N. Regarding the two overflow branches, the influence of variations in the Denmark Strait overflow is, for all models, substantially larger than that of variations in the overflow across the Iceland–Scotland Ridge. The latter might, however, be underestimated, as the models in general do not realistically simulate the flow path of the Iceland–Scotland overflow water south of the Iceland–Scotland Ridge. The influence of variations in subpolar deep water formation is, on multimodel average, larger than that of variations in the Denmark Strait overflow. This is true both at 45° N, where the maximum standard deviation of decadal to multidecadal AMOC variability is located for all but one model, and at the more classical latitude of 30° N. At 30° N, variations in subpolar deep water formation and Denmark Strait overflow explain, on multimodel average, about half and one-third respectively of the decadal to multidecadal AMOC variance. Apart from analysing multimodel control simulations, we have performed sensitivity experiments with one of the models, in which we suppress the variability of either subpolar deep water formation or Nordic Seas overflows. The sensitivity experiments indicate that variations in subpolar deep water formation and Nordic Seas overflows are not completely independent. We further conclude from these experiments that the decadal to multidecadal AMOC variability north of about 50° N is mainly related to variations in Nordic Seas overflows. At 45° N and south of this latitude, variations in both subpolar deep water formation and Nordic Seas overflows contribute to the AMOC variability, with neither of the processes being very dominant compared to the other.

Description: Pacific Ocean sea surface temperatures (SST) influence rainfall variability on multidecadal and interdecadal timescales in concert with the Pacific Decadal Oscillation (PDO) and Interdecadal Pacific Oscillation (IPO). Rainfall variations in locations such as Australia and North America are therefore linked to phase changes in the PDO. Furthermore, studies have suggested teleconnections exist between the western Indian Ocean and Pacific Decadal Variability (PDV), similar to those observed on interannual timescales related to the El Niño Southern Oscillation (ENSO). However, as instrumental records of rainfall are too short and sparse to confidently assess multidecadal climatic teleconnections, here we present four coral climate archives from Madagascar spanning up to the past 300 yr (1708–2008) to assess such decadal variability. Using spectral luminescence scanning to reconstruct past changes in river runoff, we identify significant multidecadal and interdecadal frequencies in the coral records, which before 1900 are coherent with Asian-based PDO reconstructions. This multidecadal relationship with the Asian-based PDO reconstructions points to an unidentified teleconnection mechanism that affects Madagascar rainfall/runoff, most likely triggered by multidecadal changes in North Pacific SST, influencing the Asian Monsoon circulation. In the 20th century we decouple human deforestation effects from rainfall-induced soil erosion by pairing luminescence with coral geochemistry. Positive PDO phases are associated with increased Indian Ocean temperatures and runoff/rainfall in eastern Madagascar, while precipitation in southern Africa and eastern Australia declines. Consequently, the negative PDO phase that started in 1998 may contribute to reduced rainfall over eastern Madagascar and increased precipitation in southern Africa and eastern Australia. We conclude that multidecadal rainfall variability in Madagascar and the western Indian Ocean needs to be taken into account when considering water resource management under a future warming climate.

Description: The last interglaciation (~130 to 116 ka) is a time period with a strong astronomically induced seasonal forcing of insolation compared to the present. Proxy records indicate a significantly different climate to that of the modern, in particular Arctic summer warming and higher eustatic sea level. Because the forcings are relatively well constrained, it provides an opportunity to test numerical models which are used for future climate prediction. In this paper we compile a set of climate model simulations of the early last interglaciation (130 to 125 ka), encompassing a range of model complexities. We compare the simulations to each other and to a recently published compilation of last interglacial temperature estimates. We show that the annual mean response of the models is rather small, with no clear signal in many regions. However, the seasonal response is more robust, and there is significant agreement amongst models as to the regions of warming vs cooling. However, the quantitative agreement of the model simulations with data is poor, with the models in general underestimating the magnitude of response seen in the proxies. Taking possible seasonal biases in the proxies into account improves the agreement, but only marginally. However, a lack of uncertainty estimates in the data does not allow us to draw firm conclusions. Instead, this paper points to several ways in which both modelling and data could be improved, to allow a more robust model–data comparison.

Description: A new Earth system model, the Flexible Ocean and Climate Infrastructure (FOCI), is introduced. A first version of FOCI consists of a global high-top atmosphere (ECHAM6.3) and an ocean model (NEMO3.6) as well as sea ice (LIM2) and land surface model components (JSBACH), which are coupled through the OASIS3-MCT software package. FOCI includes a number of optional modules which can be activated depending on the scientific question of interest. In the atmosphere, interactive stratospheric chemistry can be used (ECHAM6-HAMMOZ) to study, for example, the effects of the ozone hole on the climate system. In the ocean, a biogeochemistry model (MOPS) is available to study the global carbon cycle. A unique feature of FOCI is the ability to explicitly resolve mesoscale ocean eddies in specific regions. This is realized in the ocean through nesting; first examples for the Agulhas Current and the Gulf Stream systems are described here. FOCI therefore bridges the gap between coarse-resolution climate models and global high-resolution weather prediction and ocean-only models. It allows to study the evolution of the climate system on regional and seasonal to (multi-) decadal scales. The development of FOCI resulted from a combination of the long-standing expertise in ocean and climate modeling in several research units and divisions at GEOMAR. FOCI will thus be used to complement and interpret long-term observations in the Atlantic, enhance the process understanding of the role of mesoscale oceanic eddies for large-scale oceanic and atmospheric circulation patterns, study feedback mechanisms with stratospheric processes, estimate future ocean acidification, improve the simulation of the Atlantic Meridional Overturning Circulation changes and their influence on climate, ocean chemistry and biology. In this paper we present both the scientific vision for the development of FOCI as well as some technical details. This includes a first validation of the different model components using several configurations of FOCI. Results show that the model in its basic configuration runs stably under pre-industrial control as well as under historical forcing, and produces a mean climate and variability which compares well with observations, reanalysis products and other climate models. The nested configurations reduce some long-standing biases in climate models and are an important step forward to include the atmospheric response in multi-decadal eddy-rich configurations.

Description: We examine the impact of horizontal resolution and model time step on the climate of the OpenIFS version 43r3 atmospheric general circulation model. A series of simulations for the period 1979–2019 are conducted with various horizontal resolutions (i.e. ∼100, ∼50, and ∼25 km) while maintaining the same time step (i.e. 15 min) and using different time steps (i.e. 60, 30, and 15 min) at 100 km horizontal resolution. We find that the surface zonal wind bias is significantly reduced over certain regions such as the Southern Ocean and the Northern Hemisphere mid-latitudes and in tropical and subtropical regions at a high horizontal resolution (i.e. ∼25 km). Similar improvement is evident too when using a coarse-resolution model (∼100 km) with a smaller time step (i.e. 30 and 15 min). We also find improvements in Rossby wave amplitude and phase speed, as well as in weather regime patterns, when a smaller time step or higher horizontal resolution is used. The improvement in the wind bias when using the shorter time step is mostly due to an increase in shallow and mid-level convection that enhances vertical mixing in the lower troposphere. The enhanced mixing allows frictional effects to influence a deeper layer and reduces wind and wind speed throughout the troposphere. However, precipitation biases generally increase with higher horizontal resolutions or smaller time steps, whereas the surface air temperature bias exhibits a small improvement over North America and the eastern Eurasian continent. We argue that the bias improvement in the highest-horizontal-resolution (i.e. ∼25 km) configuration benefits from a combination of both the enhanced horizontal resolution and the shorter time step. In summary, we demonstrate that, by reducing the time step in the coarse-resolution (∼100 km) OpenIFS model, one can alleviate some climate biases at a lower cost than by increasing the horizontal resolution.

Description: Black shale sediments from the Barremian to Aptian South Atlantic document intense and widespread burial of marine organic carbon during the initial stages of seafloor spreading between Africa and South America. The enhanced sequestration of atmospheric CO2 makes these young ocean basins potential drivers of the Early Cretaceous carbon cycle and climate perturbations. The opening of marine gateways between initially restricted basins and related circulation and ventilation changes are a commonly invoked explanation for the transient formation and disappearance of these regional carbon sinks. However, large uncertainties in paleogeographic reconstructions limit the interpretation of available paleoceanographic data and prevent any robust model-based quantifications of the proposed circulation and carbon burial changes. Here, we present a new approach to assess the principal controls on the Early Cretaceous South Atlantic and Southern Ocean circulation changes under full consideration of the uncertainties in available boundary conditions. Specifically, we use a large ensemble of 36 climate model experiments to simulate the Barremian to Albian progressive opening of the Falkland Plateau and Georgia Basin gateways with different configurations of the proto-Drake Passage, the Walvis Ridge, and atmospheric CO2 concentrations. The experiments are designed to complement available geochemical data across the regions and to test circulation scenarios derived from them. All simulations show increased evaporation and intermediate water formation at subtropical latitudes that drive a meridional overturning circulation whose vertical extent is determined by the sill depth of the Falkland Plateau. Densest water masses formed in the southern Angola Basin and potentially reached the deep Cape Basin as Walvis Ridge Overflow Water. Paleogeographic uncertainties are as important as the lack of precise knowledge of atmospheric CO2 levels for the simulated temperature and salinity spread in large parts of the South Atlantic. Overall temperature uncertainties are up to 15 °C and increase significantly with water depth. The ensemble approach reveals temporal changes in the relative importance of geographic and radiative forcings for the simulated oceanographic conditions and, importantly, nonlinear interactions between them. Progressive northward opening of the highly restricted Angola Basin increased the sensitivity of local overturning and upper ocean stratification to atmospheric CO2 concentrations due to large-scale changes in the hydrological cycle, while the chosen proto-Drake Passage depth is critical for the ocean dynamics and CO2 response in the southern South Atlantic. Finally, the simulated processes are integrated into a recent carbon burial framework to document the principal control of the regional gateway evolution on the progressive shift from the prevailing saline and oxygen-depleted subtropical water masses to the dominance of ventilated high-latitude deep waters.

Description: Events of extreme precipitation pose a hazard to many parts of Europe but are typically not well represented in climate models. Here, we evaluate daily extreme precipitation over Europe during 1982–2019 in observations (GPCC), reanalysis (ERA5) and a set of atmosphere-only simulations at low- (100 km), medium- (50 km) and high- (25 km) horizontal resolution with identical vertical resolutions using OpenIFS (version 43r3). We find that both OpenIFS simulations and reanalysis underestimate the rates of extreme precipitation compared to observations. The biases are largest for the lowest resolution (100 km) and decrease with increasing horizontal resolution (50 and 25 km) simulations in all seasons. The sensitivity to horizontal resolution is particularly high in mountain regions (such as the Alps, Scandinavia, Iberian Peninsula), likely linked to the sensitivity of vertical velocity to the representation of topography. The sensitivity of precipitation to model resolution increases dramatically with increasing percentiles, with modest biases in the 70th–80thpercentile range and large biases above the 99th percentile range. We also find that precipitation above the 99th percentile mostly consists of large-scale precipitation (~80 %) in winter, while in summer it is mostly large-scale precipitation in Northern Europe (~70 %) and convective precipitation in Southern Europe (~70 %). Compared to ERA5, the OpenIFS overestimates large-scale precipitation extremes in winter, but underestimates in summer. The discrepancy between OpenIFS and ERA5 decreases with increasing horizontal resolutions. We also examine the sensitivity of extreme precipitation to model time step and find that the convective contribution to extreme precipitation is more sensitive to the model time step than the horizontal resolution. This is likely due to the sensitivity of convective activity to model time step. On the other hand, the large-scale contribution to extreme precipitation is more sensitive to horizontal resolution than the model time step, which may be due to sharper fronts and steeper topography at higher horizontal resolution.