A process-based analysis of ocean heat uptake in an AOGCM with an eddy-permitting ocean component

About 90% of the anthropogenic increase in heat stored in the climate system is found the oceans. Therefore it is relevant to understand the details of ocean heat uptake. Here we present a detailed, process-based analysis of ocean heat uptake (OHU) processes in HiGEM1.2, an atmosphere-ocean general circulation model (AOGCM) with an eddy-permitting ocean component of 1/3 degree resolution. Similarly to various other models, HiGEM1.2 shows that the global heat budget is dominated by a downward advection of heat compensated by upward isopycnal diffusion. Only in the upper tropical ocean do we find the classical balance between downward diapycnal diffusion and upward advection of heat. The upward isopycnal diffusion of heat is located mostly in the Southern Ocean, which thus dominates the global heat budget. We compare the responses to a 4xCO2 forcing and an enhancement of the windstress forcing in the Southern Ocean. This highlights the importance of regional processes for the global ocean heat uptake. These are mainly surface fluxes and convection in the high latitudes, and advection in the Southern Ocean mid-latitudes. Changes in diffusion are less important. In line with the CMIP5 models, HiGEM1.2 shows a band of strong OHU in the mid-latitude Southern Ocean in the 4xCO2 run, which is mostly advective. By contrast, in the high-latitude Southern Ocean regions it is the suppression of convection that leads to OHU. In the enhanced windstress run, convection is strengthened at high Southern latitudes, leading to heat loss, while the magnitude of the OHU in the Southern mid-latitudes is very similar to the 4xCO2 results. Remarkably, there is only very small global OHU in the enhanced windstress run. The wind stress forcing just leads to a redistribution of heat. We relate the ocean changes at high southern latitudes to the effect of climate change on the Antarctic Circumpolar Current (ACC). It weakens in the 4xCO2 run and strengthens in the wind stress run. The weakening is due to a narrowing of the ACC, caused by an expansion of the Weddell Gyre, and a flattening of the isopycnals, which are explained by a combination of the wind stress forcing and increased precipitation

Stochastic Stability of Open-Ocean Deep Convection

Open-ocean deep convection is a highly variable and strongly nonlinear process that plays an essential role in the global ocean circulation. A new view of its stability is presented here, in which variability, as parameterised by stochastic forcing, is central. The use of an idealised deep convection box model allows analytical solutions and straightforward conceptual understand-ing, while retaining the main features of deep convection dynamics. In contrast to the generally abrupt stability changes in deterministic systems, measures of stochastic stability change smoothly in response to varying forcing parameters. These stochastic stability measures depend chie y on the residence times of the system in dierent regions of phase space, which need not contain a stable steady state in the deterministic sense. Deep convection can occur frequently even for parameter ranges in which it is deterministically unstable; this eect is denoted wandering unimodality. The stochastic stability concepts are readily applied to other components of the climate system. The results highlight the need to take climate variability into account when analysing the stability of a climate state.

Ocean heat uptake processes: a model intercomparison

We compare the quasi-equilibrium heat balances, as well as their responses to 4×CO2 perturbation, among three global climate models with the aim to identify and explain inter-model differences in ocean heat uptake (OHU) processes. We find that, in quasi-equilibrium, convective and mixed layer processes, as well as eddy-related processes, cause cooling of the subsurface ocean. The cooling is balanced by warming caused by advective and diapycnally diffusive processes. We also find that in the CO2-perturbed climates the largest contribution to OHU comes from changes in vertical mixing processes and the mean circulation, particularly in the extra-tropics, caused both by changes in wind forcing, and by changes in high-latitude buoyancy forcing. There is a substantial warming in the tropics, a significant part of which occurs because of changes in horizontal advection in extra-tropics. Diapycnal diffusion makes only a weak contribution to the OHU, mainly in the tropics, due to increased stratification. There are important qualitative differences in the contribution of eddy-induced advection and isopycnal diffusion to the OHU among the models. The former is related to the different values of the coefficients used in the corresponding scheme. The latter is related to the different tapering formulations of the isopycnal diffusion scheme. These differences affect the OHU in the deep ocean, which is substantial in two of the models, the dominant region of deep warming being the Southern Ocean. However, most of the OHU takes place above 2000 m, and the three models are quantitatively similar in their global OHU efficiency and its breakdown among processes and as a function of latitude

Isoneutral control of effective diapycnal mixing in numerical ocean models with neutral rotated diffusion tensors

It is well known that there are infinite number of ways of constructing a globally-defined density variable for the ocean, with each possible density variable having a priori its own distinct diapycnal diffusivity. Because no globally-defined density variable can be exactly neutral, numerical ocean models tend to use rotated diffusion tensors mixing separately in the directions parallel and perpendicular to the local neutral vector at rates defined by the isoneutral and dianeutral mixing coefficients respectively. To constrain these mixing coefficients from observations, one widely used tool are inverse methods based on Walin-type water masses analyses. Such methods, however, can only constrain the diapycnal diffusivity of the globally defined density variable $\gamma$ —such as $\sigma_2$ —that underlies the inverse method. To use such a method to constrain the dianeutral mixing coefficient therefore requires understanding the relations between the different diapycnal diffusivities. However, this is complicated by the fact that the effective diapycnal diffusivity experienced by is necessarily partly controlled by isoneutral diffusion owing to the unavoidable misalignment between iso- surfaces and the neutral directions. Here, this effect is quantified by evaluating the effective diapycnal diffusion coefficient pertaining to five widely used density variables: Jackett and McDougall (1997) $\gamma^n$, Lorenz reference state density $\rho_{ref}$ of Saenz et al. (2015), and three potential density variables $\sigma_0$, $\sigma_2$ and $\sigma_4$. Computations are based on the World Ocean Circulation Experiment climatology, assuming either a uniform value for the isoneutral mixing coefficient or spatially varying values inferred from an inverse calculation. Isopycnal mixing 15 contributions to the effective diapycnal mixing yield values consistently larger than 10^(-3) m^2/s in the deep ocean for all density variables, with $\gamma^n$ suffering the least from the isoneutral control of effective diapycnal mixing, and $\sigma_0$ the most. These high values are due to spatially localised large values of non-neutrality, mostly in the deep Southern Ocean. Removing only 5% of these high values on each density surface reduces the effective diapycnal diffusivities to less than 10^(-4) m^2/s. The main implication of this work is to highlight the conceptual and practical difficulties of relating the diapycnal mixing diffusivities inferred from global budgets or inverse methods relying on Walin-like water mass analyses to locally defined dianeutral diffusivities. Doing so requires the ability to separate the relative contribution of isoneutral mixing from the effective diapycnal mixing. Because it corresponds to a special case of Walin-type water mass analysis, the determination of spurious diapycnal mixing based on monitoring the evolution of the Lorenz reference state may also be affected by the above issues when using a realistic nonlinear equation of state. The present results thus suggest that part of previously published spurious diapycnal mixing estimates could be due to isoneutral mixing contamination

A traceable physical calibration of the vertical advection-diffusion equation for modelling ocean heat uptake

The classic vertical advection-diffusion (VAD) balance is a central concept in studying the ocean heat budget, in particular in simple climate models (SCMs). Here we present a new framework to calibrate the parameters of the VAD equation to the vertical ocean heat balance of two fully-coupled climate models that is traceable to the models’ circulation as well as to vertical mixing and diffusion processes. Based on temperature diagnostics, we derive an effective vertical velocity w∗ and turbulent diffusivity k∗ for each individual physical process. In steady-state, we find that the residual vertical velocity and diffusivity change sign in mid-depth, highlighting the different regional contributions of isopycnal and diapycnal diffusion in balancing the models’ residual advection and vertical mixing. We quantify the impacts of the time-evolution of the effective quantities under a transient 1%CO2 simulation and make the link to the parameters of currently employed SCMs

The influence of eddy parameterizations on the transport of the Antarctic Circumpolar Current in coupled climate models

The transport of the Antarctic Circumpolar Current (ACC) varies strongly across the coupled GCMs (general circulation models) used for the IPCC AR4. This note shows that a large fraction of this across-model variance can be explained by relating it to the parameterization of eddy-induced transports. In the majority of models this parameterization is based on the study by Gent and McWilliams (1990). The main parameter is the quasi-Stokes diffusivity kappa (often referred to less accurately as ’’thickness diffusion’’). The ACC transport and the meridional density gradient both correlate strongly with kappa across those models where kappa is a prescribed constant. In contrast, there is no correlation with the isopycnal diffusivity jiso across the models. The sensitivity of the ACC transport to kappa is larger than to the zonal wind stress maximum. Experiments with the fast GCM FAMOUS show that changing kappa directly affects the ACC transport by changing the density structure throughout the water column. Our results suggest that this limits the role of the wind stress magnitude in setting the ACC transport in FAMOUS. The sensitivities of the ACC and the meridional density gradient are very similar across the AR4 GCMs (for those models where kappa is a prescribed constant) and among the FAMOUS experiments. The strong sensitivity of the ACC transport to kappa needs careful assessment in climate models

BGC-val: a model- and grid-independent Python toolkit to evaluate marine biogeochemical models

The biogeochemical evaluation toolkit, BGC-val, is a model- and grid-independent Python toolkit that has been built to evaluate marine biogeochemical models using a simple interface. Here, we present the ideas that motivated the development of the BGC-val software framework, introduce the code structure, and show some applications of the toolkit using model results from the Fifth Climate Model Intercomparison Project (CMIP5). A brief outline of how to access and install the repository is presented in Appendix A, but the specific details on how to use the toolkit are kept in the code repository. The key ideas that directed the toolkit design were model and grid independence, front-loading analysis functions and regional masking, interruptibility, and ease of use. We present each of these goals, why they were important, and what we did to address them. We also present an outline of the code structure of the toolkit illustrated with example plots produced by the toolkit. After describing BGC-val, we use the toolkit to investigate the performance of the marine physical and biogeochemical quantities of the CMIP5 models and highlight some predictions about the future state of the marine ecosystem under a business-as-usual CO2 concentration scenario (RCP8.5)

Earth system music: music generated from the United Kingdom Earth System Model (UKESM1)

Scientific data are almost always represented graphically in figures or in videos. With the ever-growing interest from the general public in understanding climate sciences, it is becoming increasingly important that scientists present this information in ways that are both accessible and engaging to non-experts. In this pilot study, we use time series data from the first United Kingdom Earth System Model (UKESM1) to create six procedurally generated musical pieces. Each of these pieces presents a unique aspect of the ocean component of the UKESM1, either in terms of a scientific principle or a practical aspect of modelling. In addition, each piece is arranged using a different musical progression, style and tempo. These pieces were created in the Musical Instrument Digital Interface (MIDI) format and then performed by a digital piano synthesiser. An associated video showing the time development of the data in time with the music was also created. The music and video were published on the lead author's YouTube channel. A brief description of the methodology was also posted alongside the video. We also discuss the limitations of this pilot study and describe several approaches to extend and expand upon this work

Preindustrial control simulations with HadGEM3-GC3.1 for CMIP6

Pre‐industrial control simulations with the HadGEM3‐GC3.1 climate model are presented at two resolutions. These are N216ORCA025, which has a horizontal resolution of 60km in the atmosphere and 0.25° in the ocean, and N96ORCA1, which has a horizontal resolution of 130km in the atmosphere and 1° in the ocean. The aim of this study is to document the climate variability in these simulations, make comparisons against present‐day observations (albeit under different forcing), and discuss differences arising due to resolution. In terms of interannual variability in the leading modes of climate variability the two resolutions behave generally very similarly. Notable differences are in the westward extent of El‐Niño and the pattern of Atlantic multidecadal variability, in which N216ORCA025 compares more favourably to observations, and in the Antarctic Circumpolar Current, which is far too weak in N216ORCA025. In the North Atlantic region, N216ORCA025 has a stronger and deeper AMOC, which compares well against observations, and reduced biases in temperature and salinity in the North Atlantic subpolar gyre (NA SPG). These simulations are being provided to the sixth Coupled Model Intercomparison Project (CMIP6) and provide a baseline against which further forced experiments may be assessed