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

Wind-driven oscillations in meridional overturning circulations near the equator. Part II: idealized simulations

Large-amplitude [±100 Sv (1 Sv ≡ 106 m3 s−1)], high-frequency oscillations in the Pacific Ocean’s meridional overturning circulation within 10° of the equator have been found in integrations of the NEMO ocean general circulation model. Part I of this paper showed that these oscillations are dominated by two bands of frequencies with periods close to 4 and 10 days and that they are driven by the winds within about 10° of the equator. This part shows that the oscillations can be well simulated by small-amplitude, wind-driven motions on a horizontally uniform, stably stratified state of rest. Its main novelty is that, by focusing on the zonally integrated linearized equations, it presents solutions for the motions in a basin with sloping side boundaries. The solutions are found using vertical normal modes and equatorial meridional modes representing Yanai and inertia–gravity waves. Simulations of 16-day-long segments of the time series for the Pacific of each of the first three meridional and vertical modes (nine modes in all) capture between 85% and 95% of the variance of matching time series segments diagnosed from the NEMO integrations. The best agreement is obtained by driving the solutions with the full wind forcing and the full pressure forces on the bathymetry. Similar results are obtained for the corresponding modes in the Atlantic and Indian Oceans. Slower variations in the same meridional and vertical modes of the MOC are also shown to be well simulated by a quasi-stationary solution driven by zonal wind and pressure forces

The North Atlantic subpolar circulation in an eddy-resolving global ocean model

The subpolar North Atlantic represents a key region for global climate, but most numerical models still have well-described limitations in correctly simulating the local circulation patterns. Here, we present the analysis of a 30-year run with a global eddy-resolving (1/12°) version of the NEMO ocean model. Compared to the 1° and 1/4° equivalent versions, this simulation more realistically represents the shape of the Subpolar Gyre, the position of the North Atlantic Current, and the Gulf Stream separation. Other key improvements are found in the representation of boundary currents, multi-year variability of temperature and depth of winter mixing in the Labrador Sea, and the transport of overflows at the Greenland–Scotland Ridge. However, the salinity, stratification and mean depth of winter mixing in the Labrador Sea, and the density and depth of overflow water south of the sill, still present challenges to the model. This simulation also provides further insight into the spatio-temporal development of the warming event observed in the Subpolar Gyre in the mid 1990s, which appears to coincide with a phase of increased eddy activity in the southernmost part of the gyre. This may have provided a gateway through which heat would have propagated into the gyre's interior

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

Sustained intensification of the Aleutian Low induces weak tropical Pacific sea surface warming

It has been proposed that externally forced trends in the Aleutian Low can induce a basin-wide Pacific sea surface temperature (SST) response that projects onto the pattern of the Pacific Decadal Oscillation (PDO). To investigate this hypothesis, we apply local atmospheric nudging in an intermediate-complexity climate model to isolate the effects of an intensified winter Aleutian Low sustained over several decades. An intensification of the Aleutian Low produces a basin-wide SST response with a similar pattern to the model's internally generated PDO. The amplitude of the SST response in the North Pacific is comparable to the PDO, but in the tropics and southern subtropics the anomalies induced by the imposed Aleutian Low anomaly are a factor of 3 weaker than for the internally generated PDO. The tropical Pacific warming peaks in boreal spring, though anomalies persist year-round. A heat budget analysis shows the northern subtropical Pacific SST response is predominantly driven by anomalous surface turbulent heat fluxes in boreal winter, while in the equatorial Pacific the response is mainly due to meridional heat advection in boreal spring. The propagation of anomalies from the extratropics to the tropics can be explained by the seasonal footprinting mechanism, involving the wind–evaporation–SST feedback. The results show that low-frequency variability and trends in the Aleutian Low could contribute to basin-wide anomalous Pacific SST, but the magnitude of the effect in the tropical Pacific, even for the extreme Aleutian Low forcing applied here, is small. Therefore, external forcing of the Aleutian Low is unlikely to account for observed decadal SST trends in the tropical Pacific in the late 20th and early 21st centuries

Improved estimates of water cycle change from ocean salinity: the key role of ocean warming

Changes in the global water cycle critically impact environmental, agricultural, and energy systems relied upon by humanity (Jiménez Cisneros et al 2014 Climate Change 2014: Impacts, Adaptation, and Vulnerability (Cambridge: Cambridge University Press)). Understanding recent water cycle change is essential in constraining future projections. Warming-induced water cycle change is expected to amplify the pattern of sea surface salinity (Durack et al 2012 Science 336 455–8). A puzzle has, however, emerged. The surface salinity pattern has amplified by 5%–8% since the 1950s (Durack et al 2012 Science 336 455–8, Skliris et al 2014 Clim. Dyn. 43 709–36) while the water cycle is thought to have amplified at close to half that rate (Durack et al 2012 Science 336 455–8, Skliris et al 2016 Sci. Rep. 6 752). This discrepancy is also replicated in climate projections of the 21st century (Durack et al 2012 Science 336 455–8). Using targeted numerical ocean model experiments we find that, while surface water fluxes due to water cycle change and ice mass loss amplify the surface salinity pattern, ocean warming exerts a substantial influence. Warming increases near-surface stratification, inhibiting the decay of existing salinity contrasts and further amplifying surface salinity patterns. Observed ocean warming can explain approximately half of observed surface salinity pattern changes from 1957–2016 with ice mass loss playing a minor role. Water cycle change of 3.6% ± 2.1% per degree Celsius of surface air temperature change is sufficient to explain the remaining observed salinity pattern change

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

Re-emergence of North Atlantic subsurface ocean temperature anomalies in a seasonal forecast system

A high-resolution coupled ocean atmosphere model is used to study the effects of seasonal re-emergence of North Atlantic subsurface ocean temperature anomalies on northern hemisphere winter climate. A 50-member control ensemble is integrated from 1 September 2007 to 28 February 2008 and compared with a parallel ensemble with perturbed ocean initial conditions. The perturbation consists of a density-compensated subsurface Atlantic temperature anomaly corresponding to the observed subsurface temperature anomaly for September 2010. The experiment is repeated for two atmosphere horizontal resolutions (~ 60 km and ~ 25 km) in order to determine whether the sensitivity of the atmosphere to re-emerging temperature anomalies is dependent on resolution. A wide range of re-emergence behavior is found within the perturbed ensembles. While the observations seem to indicate that most of the re-emergence is occurring in November, most members of the ensemble show re-emergence occurring later in the winter. However, when re-emergence does occur it is preceded by an atmospheric pressure pattern that induces a strong flow of cold, dry air over the mid-latitude Atlantic, and enhances oceanic latent heat loss. In response to re-emergence (negative SST anomalies), there is reduced latent heat loss, less atmospheric convection, a reduction in eddy kinetic energy and positive low-level pressure anomalies downstream. Within the framework of a seasonal forecast system the results highlight the atmospheric conditions required for re-emergence to take place and the physical processes that may lead to a significant effect on the winter atmospheric circulation

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

The Atlantic Meridional Overturning Circulation in High-Resolution Models

The Atlantic meridional overturning circulation (AMOC) represents the zonally integrated stream function of meridional volume transport in the Atlantic Basin. The AMOC plays an important role in transporting heat meridionally in the climate system. Observations suggest a heat transport by the AMOC of 1.3 PW at 26°N—a latitude which is close to where the Atlantic northward heat transport is thought to reach its maximum. This shapes the climate of the North Atlantic region as we know it today. In recent years there has been significant progress both in our ability to observe the AMOC in nature and to simulate it in numerical models. Most previous modeling investigations of the AMOC and its impact on climate have relied on models with horizontal resolution that does not resolve ocean mesoscale eddies and the dynamics of the Gulf Stream/North Atlantic Current system. As a result of recent increases in computing power, models are now being run that are able to represent mesoscale ocean dynamics and the circulation features that rely on them. The aim of this review is to describe new insights into the AMOC provided by high‐resolution models. Furthermore, we will describe how high‐resolution model simulations can help resolve outstanding challenges in our understanding of the AMOC