Tuning without over-tuning: parametric uncertainty quantification for the NEMO ocean model

In this paper we discuss climate model tuning and present an iterative automatic tuning method from the statistical science literature. The method, which we refer to here as iterative refocussing (though also known as history matching), avoids many of the common pitfalls of automatic tuning procedures that are based on optimisation of a cost function; principally the over-tuning of a climate model due to using only partial observations. This avoidance comes by seeking to rule out parameter choices that we are confident could not reproduce the observations, rather than seeking the model that is closest to them (a procedure that risks over-tuning). We comment on the state of climate model tuning and illustrate our approach through 3 waves of iterative refocussing of the NEMO ORCA2 global ocean model run at 2° resolution. We show how at certain depths the anomalies of global mean temperature and salinity in a standard configuration of the model exceeds 10 standard deviations away from observations and show the extent to which this can be alleviated by iterative refocussing without compromising model performance spatially. We show how model improvements can be achieved by simultaneously perturbing multiple parameters, and illustrate the potential of using low resolution ensembles to tune NEMO ORCA configurations at higher resolutions

Surging of global surface temperature due to decadal legacy of ocean heat uptake

Global surface warming since 1850 consisted of a series of slowdowns (hiatus) followed by surges. Knowledge of a mechanism to explain how this occurs would aid development and testing of interannual to decadal climate forecasts. In this paper a global climate model is forced to adopt an ocean state corresponding to a hiatus (with negative Interdecadal Pacific Oscillation, IPO, and other surface features typical of a hiatus) by artificially increasing the background diffusivity for a decade before restoring it to its normal value and allowing the model to evolve freely. This causes the model to develop a decadal surge which overshoots equilibrium (resulting in a positive IPO state) leaving behind a modified, warmer climate for decades. Water mass transformation diagnostics indicate that the heat budget of the tropical Pacific is a balance between large opposite signed terms: surface heating/cooling due to air-sea heat flux is balanced by vertical mixing and ocean heat transport divergence. During the artificial hiatus, excess heat becomes trapped just above the thermocline and there is a weak vertical thermal gradient (due to the high artificial background mixing). When the hiatus is terminated, by returning the background diffusivity to normal, the thermal gradient strengthens to pre-hiatus values so that the mixing (diffusivity x thermal gradient) remains roughly constant. However, since the base layer just above the thermocline remains anomalously warm this implies a warming of the entire water column above the trapped heat which results in a surge followed by a prolonged period of elevated surface temperatures

Subpolar Atlantic Ocean mixed layer heat content variability is increasingly driven by an active ocean

Cold conditions in the upper layer of the subpolar North Atlantic Ocean, at a time of pervasive warming elsewhere, have provoked significant debate. Uncertainty arises both from potential causes (surface heat loss and ocean circulation changes) and characteristic timescales (interannual to multidecadal). Resolution of these uncertainties is important as cold conditions have been linked to recent European weather extremes and a decline in the Atlantic overturning circulation. Using observations, supported by high resolution climate model analysis, we show that a surprisingly active ocean regularly generates both cold and warm interannual anomalies in addition to those generated by surface heat exchange. Furthermore, we identify distinct sea surface temperature patterns that characterise whether the ocean or atmosphere has the strongest influence in a particular year. Applying these new insights to observations, we find an increasing role for the ocean in setting North Atlantic mixed layer heat content variability since 1960

Observed and projected changes in North Atlantic seasonal temperature reduction and their drivers

The autumn-winter seasonal temperature reduction (STR) of the surface North Atlantic Ocean is investigated with control and climate change simulations of a coupled model and an observation-based sea surface temperature (SST) data set. In the climate change simulation, an increase in the magnitude of the STR is found over much of the North Atlantic, and this change is particularly marked in sea-ice affected regions and the subpolar gyre. Similar results for the mid-high latitude North Atlantic are obtained in the observational analysis. In particular, both the observation and climate model based results show that the STR has increased in magnitude by up to 0.3°C per decade in the subpolar gyre over the period 1951–2020. Drivers for the stronger STR are explored with a focus on potential contributions from increases in either ocean heat loss or the sensitivity of SST to heat loss. Over a large part of the mid-high latitude North Atlantic surface heat loss is found to have weakened in recent decades and is therefore not responsible for the stronger STR (exceptions to this are the near-coastal areas where sea-ice loss is important). In contrast, analysis of daily sensible and latent heat flux data reveals that the sensitivity of SST to heat loss has increased indicating that this term has played a major role in the stronger STR. Areas of greater SST sensitivity (and greater STR) are associated with increased surface stratification brought about predominantly by warming of the northern ocean regions

Mechanisms for late 20th and early 21st century decadal AMOC variability

Recent studies using data from the OSNAP observational campaign and from numerical ocean models suggest that the Iceland Basin and the Irminger Sea may be more significant for formation of upper North Atlantic Deep Water than the Labrador Sea. Here, we present a set of hindcast integrations of a global 1/4° NEMO simulation from 1958 until nearly the present day, forced with three standard forcing data sets. We use the surface-forced stream function, estimated from surface buoyancy fluxes, along with the overturning stream function, similarly defined in potential density space, to investigate the causal link between surface forcing and decadal variability in the strength of the Atlantic meridional overturning circulation (AMOC). We use the stream functions to demonstrate that watermasses in the simulations are transformed to higher densities as they propagate around the subpolar gyre from their formation locations in the north-east Atlantic and the Irminger Sea, consistent with the picture emerging from observations. The surface heat loss from the Irminger Sea is confirmed to be the dominant mechanism for decadal AMOC variability, with the heat loss anomaly from the Labrador Sea having about half the magnitude. A scalar metric based on the surface-forced stream function, accumulated in time, is found to be a good predictor of changes in the overturning strength. The AMOC variability is shown to be related to that of the North Atlantic Oscillation (NAO), primarily through the surface heat flux, itself dominated by the air-sea temperature difference, but also with some local feedback from the SST to the surface fluxes

Chaotic variability of the Atlantic meridional overturning circulation at subannual time scales

This study describes the intra- to interannual variability of the Atlantic meridional overturning circulation (AMOC) and the relative dynamical contributions to the total variability in an eddy-resolving 1/128 resolution ocean model. Based on a 53-yr-long hindcast and two 4-yr-long ensembles, we assess the total AMOC variability as well as the variability arising from small differences in the ocean initial state that rapidly imprints on the mesoscale eddy fields and subsequently on large-scale features. This initial-condition-dependent variability will henceforth be referred to as “chaotic” variability. We find that intra-annual AMOC fluctuations are mainly driven by the atmospheric forcing, with the chaotic variability fraction never exceeding 26% of the total variance in the whole meridional Atlantic domain. To understand the nature of the chaotic variability we decompose the AMOC (into its Ekman, geostrophic, barotropic, and residual components). The barotropic and geostrophic AMOC contributions exhibit strong, partly compensating fluctuations, which are linked to chaotic spatial variations of currents over topography. In the North Atlantic, the largest chaotic divergence of ensemble members is found around 248, 388, and 648N. At 26.58N, where the AMOC is monitored by the RAPID– MOCHA array, the chaotic fraction of the AMOC variability is 10%. This fraction is slightly overestimated with the reconstruction methodology as used in the observations (∼15%). This higher fraction of chaotic variability is due to the barotropic contribution not being completely captured by the monitoring system. We look at the strong AMOC decline observed in 2009/10 and find that the ensemble spread (our measure for chaotic variability) was not particularly large during this event

Interactions between the stratospheric polar vortex and Atlantic circulation on seasonal to multi-decadal timescales

Variations in the strength of the Northern Hemisphere winter polar stratospheric vortex can influence surface variability in the Atlantic sector. Disruptions of the vortex, known as sudden stratospheric warmings (SSWs), are associated with an equatorward shift and deceleration of the North Atlantic jet stream, negative phases of the North Atlantic Oscillation, and cold snaps over Eurasia and North America. Despite clear influences at the surface on sub-seasonal timescales, how stratospheric vortex variability interacts with ocean circulation on decadal to multi-decadal timescales is less well understood. In this study, we use a 1000 year preindustrial control simulation of the UK Earth System Model to study such interactions, using a wavelet analysis technique to examine non-stationary periodic signals in the vortex and ocean. We find that intervals which exhibit persistent anomalous vortex behaviour lead to oscillatory responses in the Atlantic Meridional Overturning Circulation (AMOC). The origin of these responses appears to be highly non-stationary, with spectral power in vortex variability at periods of 30 and 50 years. In contrast, AMOC variations on longer timescales (near 90-year periods) are found to lead to a vortex response through a pathway involving the equatorial Pacific and quasi-biennial oscillation. Using the relationship between persistent vortex behaviour and the AMOC response established in the model, we use regression analysis to estimate the potential contribution of the 8-year SSW hiatus interval in the 1990s to the recent negative trend in AMOC observations. The result suggests that approximately 30 % of the trend may have been caused by the SSW hiatus

FORTE 2.0: a fast, parallel and flexible coupled climate model

FORTE 2.0 is an intermediate-resolution coupled atmosphere–ocean general circulation model (AOGCM) consisting of the Intermediate General Circulation Model 4 (IGCM4), a T42 spectral atmosphere with 35σ layers, coupled to Modular Ocean Model – Array (MOMA), a 2∘ × 2∘ ocean with 15 z-layer depth levels. Sea ice is represented by a simple flux barrier. Both the atmosphere and ocean components are coded in Fortran. It is capable of producing a stable climate for long integrations without the need for flux adjustments. One flexibility afforded by the IGCM4 atmosphere is the ability to configure the atmosphere with either 35σ layers (troposphere and stratosphere) or 20σ layers (troposphere only). This enables experimental designs for exploring the roles of the troposphere and stratosphere, and the faster integration of the 20σ layer configuration enables longer duration studies on modest hardware. A description of FORTE 2.0 is given, followed by the analysis of two 2000-year control integrations, one using the 35σ configuration of IGCM4 and one using the 20σ configuration

Drivers of exceptionally cold North Atlantic Ocean temperatures and their link to the 2015 European heat wave

The North Atlantic and Europe experienced two extreme climate events in 2015: exceptionally cold ocean surface temperatures and a summer heat wave ranked in the top ten over the past 65 years. Here, we show that the cold ocean temperatures were the most extreme in the modern record over much of the mid-high latitude North-East Atlantic. Further, by considering surface heat loss, ocean heat content and wind driven upwelling we explain for the first time the genesis of this cold ocean anomaly. We find that it is primarily due to extreme ocean heat loss driven by atmospheric circulation changes in the preceding two winters combined with the re-emergence of cold ocean water masses. Furthermore, we reveal that a similar cold Atlantic anomaly was also present prior to the most extreme European heat waves since the 1980s indicating that it is a common factor in the development of these events. For the specific case of 2015, we show that the ocean anomaly is linked to a stationary position of the Jet Stream that favours the development of high surface temperatures over Central Europe during the heat wave. Our study calls for an urgent assessment of the impact of ocean drivers on major European summer temperature extremes in order to provide better advance warning measures of these high societal impact events

Labrador Sea subsurface density as a precursor of multidecadal variability in the North Atlantic: a multi-model study

The subpolar North Atlantic (SPNA) is a region with prominent decadal variability that has experienced remarkable warming and cooling trends in the last few decades. These observed trends have been preceded by slow-paced increases and decreases in the Labrador Sea density (LSD), which are thought to be a precursor of large-scale ocean circulation changes. This article analyses the interrelationships between the LSD and the wider North Atlantic across an ensemble of coupled climate model simulations. In particular, it analyses the link between subsurface density and the deep boundary density, the Atlantic Meridional Overturning Circulation (AMOC), the subpolar gyre (SPG) circulation, and the upper-ocean temperature in the eastern SPNA. All simulations exhibit considerable multidecadal variability in the LSD and the ocean circulation indices, which are found to be interrelated. LSD is strongly linked to the strength of the subpolar AMOC and gyre circulation, and it is also linked to the subtropical AMOC, although the strength of this relationship is model-dependent and affected by the inclusion of the Ekman component. The connectivity of LSD with the subtropics is found to be sensitive to different model features, including the mean density stratification in the Labrador Sea, the strength and depth of the AMOC, and the depth at which the LSD propagates southward along the western boundary. Several of these quantities can also be computed from observations, and comparison with these observation-based quantities suggests that models representing a weaker link to the subtropical AMOC might be more realistic