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

A window on the deep ocean: the special value of ocean bottom pressure for monitoring the large-scale, deep-ocean circulation

We show how, by focusing on bottom pressure measurements particularly on the global continental slope, it is possible to avoid the “fog” of mesoscale variability which dominates most observables in the deep ocean. This makes it possible to monitor those aspects of the ocean circulation which are most important for global scale ocean variability and climate. We therefore argue that such measurements should be considered an important future component of the Global Ocean Observing System, to complement the present open-ocean and coastal elements. Our conclusions are founded on both theoretical arguments, and diagnostics from a fine-resolution ocean model that has realistic amplitudes and spectra of mesoscale variability. These show that boundary pressure variations are coherent over along-slope distances of tens of thousands of kilometres, for several vertical modes. We illustrate the value of this in the model Atlantic, by determining the time for boundary and equatorial waves to complete a circuit of the northern basin (115 and 205 days for the first and second vertical modes), showing how the boundary features compare with basin-scale theoretical models, and demonstrating the ability to monitor the meridional overturning circulation using these boundary measurements. Finally, we discuss applicability to the real ocean and make recommendations on how to make such measurements without contamination from instrumental drift

The low-resolution version of HadGEM3 GC3.1: development and evaluation for global climate

A new climate model, HadGEM3 N96ORCA1, is presented that is part of the GC3.1 configuration of HadGEM3. N96ORCA1 has a horizontal resolution of ~135 km in the atmosphere and 1° in the ocean and requires an order of magnitude less computing power than its medium-resolution counterpart, N216ORCA025, while retaining a high degree of performance traceability. Scientific performance is compared both to observations and the N216ORCA025 model. N96ORCA1 reproduces observed climate mean and variability almost as well as N216ORCA025. Patterns of biases are similar across the two models. In the north-west Atlantic, N96ORCA1 shows a cold surface bias of up to 6K, typical of ocean models of this resolution. The strength of the Atlantic meridional overturning circulation (16 to 17 Sv) matches observations. In the Southern Ocean, a warm surface bias (up to 2K) is smaller than in N216ORCA025 and linked to improved ocean circulation. Model El Niño/Southern Oscillation and Atlantic Multidecadal Variability are close to observations. Both the cold bias in the Northern hemisphere (N96ORCA1) and the warm bias in the Southern hemisphere (N216ORCA025) develop in the first few decades of the simulations. As in many comparable climate models, simulated interhemispheric gradients of top-of-atmosphere radiation are larger than observations suggest, with contributions from both hemispheres. HadGEM3 GC3.1 N96ORCA1 constitutes the physical core of the UK Earth System Model (UKESM1) and will be used extensively in the Coupled Model Intercomparison Project 6 (CMIP6), both as part of UKESM1 and as a stand-alone coupled climate model

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

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

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