Decadal climate prediction (project GCEP)

Decadal prediction uses climate models forced by changing greenhouse gases, as in the International Panel for Climate Change, but unlike longer range predictions they also require initialization with observations of the current climate. In particular, the upper-ocean heat content and circulation have a critical influence. Decadal prediction is still in its infancy and there is an urgent need to understand the important processes that determine predictability on these timescales. We have taken the first Hadley Centre Decadal Prediction System (DePreSys) and implemented it on several NERC institute compute clusters in order to study a wider range of initial condition impacts on decadal forecasting, eventually including the state of the land and cryosphere. The eScience methods are used to manage submission and output from the many ensemble model runs required to assess predictive skill. Early results suggest initial condition skill may extend for several years, even over land areas, but this depends sensitively on the definition used to measure skill, and alternatives are presented. The Grid for Coupled Ensemble Prediction (GCEP) system will allow the UK academic community to contribute to international experiments being planned to explore decadal climate predictability.

A new perspective on warming of the global oceans

Changes in ocean circulation associated with internal climate variability have a major influence on upper ocean temperatures, particularly in regions such as the North Atlantic, which are relatively well-observed and therefore over-represented in the observational record. As a result, global estimates of upper ocean heat content can give misleading estimates of the roles of natural and anthropogenic factors in causing oceanic warming. We present a method to quantify ocean warming that filters out the natural internal variability from both observations and climate simulations and better isolates externally forced air-sea heat flux changes. We obtain a much clearer picture of the drivers of oceanic temperature changes, being able to detect the effects of both anthropogenic and volcanic influences simultaneously in the observed record. Our results show that climate models are capable of capturing in remarkable detail the externally forced component of ocean temperature evolution over the last five decades.

The impact of finer-resolution air-sea coupling on the intraseasonal oscillation of the Indian monsoon

The new HadKPP atmosphere–ocean coupled model is described and then used to determine the effects of sub-daily air–sea coupling and fine near-surface ocean vertical resolution on the representation of the Northern Hemisphere summer intra-seasonal oscillation. HadKPP comprises the Hadley Centre atmospheric model coupled to the K Profile Parameterization ocean-boundary-layer model. Four 30-member ensembles were performed that varied in oceanic vertical resolution between 1 m and 10 m and in coupling frequency between 3 h and 24 h. The 10 m, 24 h ensemble exhibited roughly 60% of the observed 30–50 day variability in sea-surface temperatures and rainfall and very weak northward propagation. Enhancing either only the vertical resolution or only the coupling frequency produced modest improvements in variability and only a standing intra-seasonal oscillation. Only the 1 m, 3 h configuration generated organized, northward-propagating convection similar to observations. Sub-daily surface forcing produced stronger upper-ocean temperature anomalies in quadrature with anomalous convection, which likely affected lower-atmospheric stability ahead of the convection, causing propagation. Well-resolved air–sea coupling did not improve the eastward propagation of the boreal summer intra-seasonal oscillation in this model. Upper-ocean vertical mixing and diurnal variability in coupled models must be improved to accurately resolve and simulate tropical sub-seasonal variability. In HadKPP, the mere presence of air–sea coupling was not sufficient to generate an intra-seasonal oscillation resembling observations.

The effect of Arabian Sea optical properties on SST biases and the South Asian summer monsoon in a coupled GCM

This study examines the effect of seasonally varying chlorophyll on the climate of the Arabian Sea and South Asian monsoon. The effect of such seasonality on the radiative properties of the upper ocean is often a missing process in coupled general circulation models and its large amplitude in the region makes it a pertinent choice for study to determine any impact on systematic biases in the mean and seasonality of the Arabian Sea. In this study we examine the effects of incorporating a seasonal cycle in chlorophyll due to phytoplankton blooms in the UK Met Office coupled atmosphere-ocean GCM HadCM3. This is achieved by performing experiments in which the optical properties of water in the Arabian Sea - a key signal of the semi-annual cycle of phytoplankton blooms in the region - are calculated from a chlorophyll climatology derived from Sea-viewing Wide Field-of-View Sensor (SeaWiFS) data. The SeaWiFS chlorophyll is prescribed in annual mean and seasonally-varying experiments. In response to the chlorophyll bloom in late spring, biases in mixed layer depth are reduced by up to 50% and the surface is warmed, leading to increases in monsoon rainfall during the onset period. However when the monsoons are fully established in boreal winter and summer and there are strong surface winds and a deep mixed layer, biases in the mixed layer depth are reduced but the surface undergoes cooling. The seasonality of the response of SST to chlorophyll is found to depend on the relative depth of the mixed layer to that of the anomalous penetration depth of solar fluxes. Thus the inclusion of the effects of chlorophyll on radiative properties of the upper ocean acts to reduce biases in mixed layer depth and increase seasonality in SST.

Wind-driven mixing below the oceanic mixed layer

This study describes the turbulent processes in the upper ocean boundary layer forced by a constant surface stress in the absence of the Coriolis force using large-eddy simulation. The boundary layer that develops has a two-layer structure, a well-mixed layer above a stratified shear layer. The depth of the mixed layer is approximately constant, whereas the depth of the shear layer increases with time. The turbulent momentum flux varies approximately linearly from the surface to the base of the shear layer. There is a maximum in the production of turbulence through shear at the base of the mixed layer. The magnitude of the shear production increases with time. The increase is mainly a result of the increase in the turbulent momentum flux at the base of the mixed layer due to the increase in the depth of the boundary layer. The length scale for the shear turbulence is the boundary layer depth. A simple scaling is proposed for the magnitude of the shear production that depends on the surface forcing and the average mixed layer current. The scaling can be interpreted in terms of the divergence of a mean kinetic energy flux. A simple bulk model of the boundary layer is developed to obtain equations describing the variation of the mixed layer and boundary layer depths with time. The model shows that the rate at which the boundary layer deepens does not depend on the stratification of the thermocline. The bulk model shows that the variation in the mixed layer depth is small as long as the surface buoyancy flux is small.

On the distortion of turbulence by a progressive surface wave

A rapid-distortion model is developed to investigate the interaction of weak turbulence with a monochromatic irrotational surface water wave. The model is applicable when the orbital velocity of the wave is larger than the turbulence intensity, and when the slope of the wave is sufficiently high that the straining of the turbulence by the wave dominates over the straining of the turbulence by itself. The turbulence suffers two distortions. Firstly, vorticity in the turbulence is modulated by the wave orbital motions, which leads to the streamwise Reynolds stress attaining maxima at the wave crests and minima at the wave troughs; the Reynolds stress normal to the free surface develops minima at the wave crests and maxima at the troughs. Secondly, over several wave cycles the Stokes drift associated with the wave tilts vertical vorticity into the horizontal direction, subsequently stretching it into elongated streamwise vortices, which come to dominate the flow. These results are shown to be strikingly different from turbulence distorted by a mean shear flow, when `streaky structures' of high and low streamwise velocity fluctuations develop. It is shown that, in the case of distortion by a mean shear flow, the tendency for the mean shear to produce streamwise vortices by distortion of the turbulent vorticity is largely cancelled by a distortion of the mean vorticity by the turbulent fluctuations. This latter process is absent in distortion by Stokes drift, since there is then no mean vorticity. The components of the Reynolds stress and the integral length scales computed from turbulence distorted by Stokes drift show the same behaviour as in the simulations of Langmuir turbulence reported by McWilliams, Sullivan & Moeng (1997). Hence we suggest that turbulent vorticity in the upper ocean, such as produced by breaking waves, may help to provide the initial seeds for Langmuir circulations, thereby complementing the shear-flow instability mechanism developed by Craik & Leibovich (1976). The tilting of the vertical vorticity into the horizontal by the Stokes drift tends also to produce a shear stress that does work against the mean straining associated with the wave orbital motions. The turbulent kinetic energy then increases at the expense of energy in the wave. Hence the wave decays. An expression for the wave attenuation rate is obtained by scaling the equation for the wave energy, and is found to be broadly consistent with available laboratory data.

Have aerosols caused the observed Atlantic multidecadal variability?

Identifying the prime drivers of the twentieth-century multidecadal variability in the Atlantic Ocean is crucial for predicting how the Atlantic will evolve in the coming decades and the resulting broad impacts on weather and precipitation patterns around the globe. Recently, Booth et al. showed that the Hadley Centre Global Environmental Model, version 2, Earth system configuration (HadGEM2-ES) closely reproduces the observed multidecadal variations of area-averaged North Atlantic sea surface temperature in the twentieth century. The multidecadal variations simulated in HadGEM2-ES are primarily driven by aerosol indirect effects that modify net surface shortwave radiation. On the basis of these results, Booth et al. concluded that aerosols are a prime driver of twentieth-century North Atlantic climate variability. However, here it is shown that there are major discrepancies between the HadGEM2-ES simulations and observations in the North Atlantic upper-ocean heat content, in the spatial pattern of multidecadal SST changes within and outside the North Atlantic, and in the subpolar North Atlantic sea surface salinity. These discrepancies may be strongly influenced by, and indeed in large part caused by, aerosol effects. It is also shown that the aerosol effects simulated in HadGEM2-ES cannot account for the observed anticorrelation between detrended multidecadal surface and subsurface temperature variations in the tropical North Atlantic. These discrepancies cast considerable doubt on the claim that aerosol forcing drives the bulk of this multidecadal variability.

The role of oceans and sea ice in abrupt transitions between multiple climate states

The coupled climate dynamics underlying large, rapid, and potentially irreversible changes in ice cover are studied. A global atmosphere–ocean–sea ice general circulation model with idealized aquaplanet geometry is forced by gradual multi-millennial variations in solar luminosity. The model traverses a hysteresis loop between warm ice-free conditions and cold glacial conditions in response to ±5 W m−2 variations in global, annual-mean insolation. Comparison of several model configurations confirms the importance of polar ocean processes in setting the sensitivity and time scales of the transitions. A “sawtooth” character is found with faster warming and slower cooling, reflecting the opposing effects of surface heating and cooling on upper-ocean buoyancy and, thus, effective heat capacity. The transition from a glacial to warm, equable climate occurs in about 200 years. In contrast to the “freshwater hosing” scenario, transitions are driven by radiative forcing and sea ice feedbacks. The ocean circulation, and notably the meridional overturning circulation (MOC), does not drive the climate change. The MOC (and associated heat transport) collapses poleward of the advancing ice edge, but this is a purely passive response to cooling and ice expansion. The MOC does, however, play a key role in setting the time scales of the transition and contributes to the asymmetry between warming and cooling.

The GODAE high resolution sea surface temperature pilot project

Sea surface temperature (SST) measurements are required by operational ocean and atmospheric forecasting systems to constrain modeled upper ocean circulation and thermal structure. The Global Ocean Data Assimilation Experiment (GODAE) High Resolution SST Pilot Project (GHRSST-PP) was initiated to address these needs by coordinating the provision of accurate, high-resolution, SST products for the global domain. The pilot project is now complete, but activities continue within the Group for High Resolution SST (GHRSST). The pilot project focused on harmonizing diverse satellite and in situ data streams that were indexed, processed, quality controlled, analyzed, and documented within a Regional/Global Task Sharing (R/GTS) framework implemented in an internationally distributed manner. Data with meaningful error estimates developed within GHRSST are provided by services within R/GTS. Currently, several terabytes of data are processed at international centers daily, creating more than 25 gigabytes of product. Ensemble SST analyses together with anomaly SST outputs are generated each day, providing confidence in SST analyses via diagnostic outputs. Diagnostic data sets are generated and Web interfaces are provided to monitor the quality of observation and analysis products. GHRSST research and development projects continue to tackle problems of instrument calibration, algorithm development, diurnal variability, skin temperature deviation, and validation/verification of GHRSST products. GHRSST also works closely with applications and users, providing a forum for discussion and feedback between SST users and producers on a regular basis. All data within the GHRSST R/GTS framework are freely available. This paper reviews the progress of GHRSST-PP, highlighting achievements that have been fundamental to the success of the pilot project.

Impact of brine-induced stratification on the glacial carbon cycle

During the cold period of the Last Glacial Maximum (LGM, about 21 000 years ago) atmospheric CO2 was around 190 ppm, much lower than the pre-industrial concentration of 280 ppm. The causes of this substantial drop remain partially unresolved, despite intense research. Understanding the origin of reduced atmospheric CO2 during glacial times is crucial to comprehend the evolution of the different carbon reservoirs within the Earth system (atmosphere, terrestrial biosphere and ocean). In this context, the ocean is believed to play a major role as it can store large amounts of carbon, especially in the abyss, which is a carbon reservoir that is thought to have expanded during glacial times. To create this larger reservoir, one possible mechanism is to produce very dense glacial waters, thereby stratifying the deep ocean and reducing the carbon exchange between the deep and upper ocean. The existence of such very dense waters has been inferred in the LGM deep Atlantic from sediment pore water salinity and δ18O inferred temperature. Based on these observations, we study the impact of a brine mechanism on the glacial carbon cycle. This mechanism relies on the formation and rapid sinking of brines, very salty water released during sea ice formation, which brings salty dense water down to the bottom of the ocean. It provides two major features: a direct link from the surface to the deep ocean along with an efficient way of setting a strong stratification. We show with the CLIMBER-2 carbon-climate model that such a brine mechanism can account for a significant decrease in atmospheric CO2 and contribute to the glacial-interglacial change. This mechanism can be amplified by low vertical diffusion resulting from the brine-induced stratification. The modeled glacial distribution of oceanic δ13C as well as the deep ocean salinity are substantially improved and better agree with reconstructions from sediment cores, suggesting that such a mechanism could have played an important role during glacial times.

A continuum model of melt pond evolution on Arctic sea ice

[1] During the Northern Hemisphere summer, absorbed solar radiation melts snow and the upper surface of Arctic sea ice to generate meltwater that accumulates in ponds. The melt ponds reduce the albedo of the sea ice cover during the melting season, with a significant impact on the heat and mass budget of the sea ice and the upper ocean. We have developed a model, designed to be suitable for inclusion into a global circulation model (GCM), which simulates the formation and evolution of the melt pond cover. In order to be compatible with existing GCM sea ice models, our melt pond model builds upon the existing theory of the evolution of the sea ice thickness distribution. Since this theory does not describe the topography of the ice cover, which is crucial to determining the location, extent, and depth of individual ponds, we have needed to introduce some assumptions. We describe our model, present calculations and a sensitivity analysis, and discuss our results.

On the relationship between decadal buoyancy anomalies and variability of the Atlantic meridional overturning circulation

Owing to the role of the Atlantic meridional overturning circulation (AMOC) in ocean heat transport, AMOC variability is thought to play a role in climate variability on a wide range of time scales. This paper focuses on the potential role of the AMOC in climate variability on decadal time scales. Coupled and ocean-only general circulation models run in idealized geometries are utilized to study the relationships between decadal AMOC and buoyancy variability and determine whether the AMOC plays an active role in setting sea surface temperature on decadal time scales.DecadalAMOC variability is related to changes in the buoyancy field along the western boundary according to the thermal wind relation. Buoyancy anomalies originate in the upper ocean of the subpolar gyre and travel westward as baroclinic Rossby waves. When the buoyancy anomalies strike the western boundary, they are advected southward by the deep western boundary current, leading to latitudinally coherent AMOC variability. The AMOC is observed to respond passively to decadal buoyancy anomalies: although variability of the AMOC leads to meridional ocean heat transport anomalies, these transports are not responsible for creating the buoyancy anomalies in the subpolar gyre that drive AMOC variability.

On the observability of turbulent transport rates by Argo: supporting evidence from an inversion experiment

Although estimation of turbulent transport parameters using inverse methods is not new, there is little evaluation of the method in the literature. Here, it is shown that extended observation of the broad scale hydrography by Argo provides a path to improved estimates of regional turbulent transport rates. Results from a 20 year ocean state estimate produced with the ECCO v4 non-linear inverse modeling framework provide supporting evidence. Turbulent transport parameter maps are estimated under the constraints of fitting the extensive collection of Argo profiles collected through 2011. The adjusted parameters dramatically reduce misfits to in situ profiles as compared with earlier ECCO solutions. They also yield a clear reduction in the model drift away from observations over multi-century long simulations, both for assimilated variables (temperature and salinity) and independent variables (bio-geochemical tracers). Despite the minimal constraints imposed specifically on the estimated parameters, their geography is physically plausible and exhibits close connections with the upper ocean ocean stratification as observed by Argo. The estimated parameter adjustments furthermore have first order impacts on upper-ocean stratification and mixed layer depths over 20 years. These results identify the constraint of fitting Argo profiles as an effective observational basis for regional turbulent transport rates. Uncertainties and further improvements of the method are discussed.

Towards predictive understanding of regional climate change

Regional information on climate change is urgently needed but often deemed unreliable. To achieve credible regional climate projections, it is essential to understand underlying physical processes, reduce model biases and evaluate their impact on projections, and adequately account for internal variability. In the tropics, where atmospheric internal variability is small compared with the forced change, advancing our understanding of the coupling between long-term changes in upper-ocean temperature and the atmospheric circulation will help most to narrow the uncertainty. In the extratropics, relatively large internal variability introduces substantial uncertainty, while exacerbating risks associated with extreme events. Large ensemble simulations are essential to estimate the probabilistic distribution of climate change on regional scales. Regional models inherit atmospheric circulation uncertainty from global models and do not automatically solve the problem of regional climate change. We conclude that the current priority is to understand and reduce uncertainties on scales greater than 100 km to aid assessments at finer scales.

Atmospheric and oceanic contributions to irreducible forecast uncertainty of Arctic surface climate

Uncertainty of Arctic seasonal to interannual predictions arising from model errors and initial state uncertainty has been widely discussed in the literature, whereas the irreducible forecast uncertainty (IFU) arising from the chaoticity of the climate system has received less attention. However, IFU provides important insights into the mechanisms through which predictability is lost, and hence can inform prioritization of model development and observations deployment. Here, we characterize how internal oceanic and surface atmospheric heat fluxes contribute to IFU of Arctic sea ice and upper ocean heat content in an Earth system model by analyzing a set of idealized ensemble prediction experiments. We find that atmospheric and oceanic heat flux are often equally important for driving unpredictable Arctic-wide changes in sea ice and surface water temperatures, and hence contribute equally to IFU. Atmospheric surface heat flux tends to dominate Arctic-wide changes for lead times of up to a year, whereas oceanic heat flux tends to dominate regionally and on interannual time scales. There is in general a strong negative covariance between surface heat flux and ocean vertical heat flux at depth, and anomalies of lateral ocean heat transport are wind-driven, which suggests that the unpredictable oceanic heat flux variability is mainly forced by the atmosphere. These results are qualitatively robust across different initial states, but substantial variations in the amplitude of IFU exist. We conclude that both atmospheric variability and the initial state of the upper ocean are key ingredients for predictions of Arctic surface climate on seasonal to interannual time scales.