An ensemble of eddy-permitting global ocean reanalyses from the MyOcean project

A set of four eddy-permitting global ocean reanalyses produced in the framework of the MyOcean project have been compared over the altimetry period 1993–2011. The main differences among the reanalyses used here come from the data assimilation scheme implemented to control the ocean state by inserting reprocessed observations of sea surface temperature (SST), in situ temperature and salinity profiles, sea level anomaly and sea-ice concentration. A first objective of this work includes assessing the interannual variability and trends for a series of parameters, usually considered in the community as essential ocean variables: SST, sea surface salinity, temperature and salinity averaged over meaningful layers of the water column, sea level, transports across pre-defined sections, and sea ice parameters. The eddy-permitting nature of the global reanalyses allows also to estimate eddy kinetic energy. The results show that in general there is a good consistency between the different reanalyses. An intercomparison against experiments without data assimilation was done during the MyOcean project and we conclude that data assimilation is crucial for correctly simulating some quantities such as regional trends of sea level as well as the eddy kinetic energy. A second objective is to show that the ensemble mean of reanalyses can be evaluated as one single system regarding its reliability in reproducing the climate signals, where both variability and uncertainties are assessed through the ensemble spread and signal-to-noise ratio. The main advantage of having access to several reanalyses differing in the way data assimilation is performed is that it becomes possible to assess part of the total uncertainty. Given the fact that we use very similar ocean models and atmospheric forcing, we can conclude that the spread of the ensemble of reanalyses is mainly representative of our ability to gauge uncertainty in the assimilation methods. This uncertainty changes a lot from one ocean parameter to another, especially in global indices. However, despite several caveats in the design of the multi-system ensemble, the main conclusion from this study is that an eddy-permitting multi-system ensemble approach has become mature and our results provide a first step towards a systematic comparison of eddy-permitting global ocean reanalyses aimed at providing robust conclusions on the recent evolution of the oceanic state.

Annular modes and apparent eddy feedbacks in the Southern Hemisphere

Lagged correlation analysis is often used to infer intraseasonal dynamical effects but is known to be affected by non-stationarity. We highlight a pronounced quasi-two-year peak in the anomalous zonal wind and eddy momentum flux convergence power spectra in the Southern Hemisphere, which is prima facie evidence for non-stationarity. We then investigate the consequences of this non-stationarity for the Southern Annular Mode and for eddy momentum flux convergence. We argue that positive lagged correlations previously attributed to the existence of an eddy feedback are more plausibly attributed to non-stationary interannual variability external to any potential feedback process in the mid-latitude troposphere. The findings have implications for the diagnosis of feedbacks in both models and re-analysis data as well as for understanding the mechanisms underlying variations in the zonal wind.

Inter-comparison of stratospheric mean-meridional circulation and eddy mixing among six reanalysis data sets

The stratospheric mean-meridional circulation (MMC) and eddy mixing are compared among six meteorological reanalysis data sets: NCEP-NCAR, NCEP-CFSR, ERA-40, ERA-Interim, JRA-25, and JRA-55 for the period 1979–2012. The reanalysis data sets produced using advanced systems (i.e., NCEP-CFSR, ERA-Interim, and JRA-55) generally reveal a weaker MMC in the Northern Hemisphere (NH) compared with those produced using older systems (i.e., NCEP/NCAR, ERA-40, and JRA-25). The mean mixing strength differs largely among the data products. In the NH lower stratosphere, the contribution of planetary-scale mixing is larger in the new data sets than in the old data sets, whereas that of small-scale mixing is weaker in the new data sets. Conventional data assimilation techniques introduce analysis increments without maintaining physical balance, which may have caused an overly strong MMC and spurious small-scale eddies in the old data sets. At the NH mid-latitudes, only ERA-Interim reveals a weakening MMC trend in the deep branch of the Brewer–Dobson circulation (BDC). The relative importance of the eddy mixing compared with the mean-meridional transport in the subtropical lower stratosphere shows increasing trends in ERA-Interim and JRA-55; this together with the weakened MMC in the deep branch may imply an increasing age-of-air (AoA) in the NH middle stratosphere in ERA-Interim. Overall, discrepancies between the different variables and trends therein as derived from the different reanalyses are still relatively large, suggesting that more investments in these products are needed in order to obtain a consolidated picture of observed changes in the BDC and the mechanisms that drive them.

Solar influences on dynamical coupling between the stratosphere and troposphere

We use a simplified atmospheric general circulation model (AGCM) to investigate the response of the lower atmosphere to thermal perturbations in the lower stratosphere. The results show that generic heating of the lower stratosphere tends to weaken the sub-tropical jets and the tropospheric mean meridional circulations. The positions of the jets, and the extent of the Hadley cells, respond to the distribution of the stratospheric heating, with low latitude heating displacing them poleward, and uniform heating displacing them equatorward. The patterns of response to the low latitude heating are similar to those found to be associated with solar variability in previous observational data analysis, and to the effects of varying solar UV radiation in sophisticated AGCMs. In order to investigate the chain of causality involved in converting the stratospheric thermal forcing to a tropospheric climate signal we conduct an experiment which uses an ensemble of model spin-ups to analyse the time development of the response to an applied stratospheric perturbation. We find that the initial effect of the change in static stability at the tropopause is to reduce the eddy momentum flux convergence in this region. This is followed by a vertical transfer of the momentum forcing anomaly by an anomalous mean circulation to the surface, where it is partly balanced by surface stress anomalies. The unbalanced part drives the evolution of the vertically integrated zonal flow. We conclude that solar heating of the stratosphere may produce changes in the circulation of the troposphere even without any direct forcing below the tropopause. We suggest that the impact of the stratospheric changes on wave propagation is key to the mechanisms involved.

Langmuir turbulence and deeply penetrating jets in an unstratified mixed layer

The influence of surface waves and an applied wind stress is studied in an ensemble of large eddy simulations to investigate the nature of deeply penetrating jets into an unstratified mixed layer. The influence of a steady monochromatic surface wave propagating parallel to the wind direction is parameterized using the wave-filtered Craik-Leibovich equations. Tracer trajectories and instantaneous downwelling velocities reveal classic counterrotating Langmuir rolls. The associated downwelling jets penetrate to depths in excess of the wave's Stokes depth scale, δs. Qualitative evidence suggests the depth of the jets is controlled by the Ekman depth scale. Analysis of turbulent kinetic energy (tke) budgets reveals a dynamical distinction between Langmuir turbulence and shear-driven turbulence. In the former, tke production is dominated by Stokes shear and a vertical flux term transports tke to a depth where it is dissipated. In the latter, tke production is from the mean shear and is locally balanced by dissipation. We define the turbulent Langmuir number Lat = (v*/Us)0.5 (v* is the ocean's friction velocity and Us is the surface Stokes drift velocity) and a turbulent anisotropy coefficient Rt = /( + ). The transition between shear-driven and Langmuir turbulence is investigated by varying external wave parameters δs and Lat and by diagnosing Rt and the Eulerian mean and Stokes shears. When either Lat or δs are sufficiently small the Stokes shear dominates the mean shear and the flow is preconditioned to Langmuir turbulence and the associated deeply penetrating jets.

The Storm-Track Response to Idealized SST Perturbations in an Aquaplanet GCM

The tropospheric response to midlatitude SST anomalies has been investigated through a series of aquaplanet simulations using a high-resolution version of the Hadley Centre atmosphere model (HadAM3) under perpetual equinox conditions. Model integrations show that increases in the midlatitude SST gradient generally lead to stronger storm tracks that are shifted slightly poleward, consistent with changes in the lower-tropospheric baroclinicity. The large-scale atmospheric response is, however, highly sensitive to the position of the SST gradient anomaly relative to that of the subtropical jet in the unperturbed atmosphere. In particular, when SST gradients are increased very close to the subtropical jet, then the Hadley cell and subtropical jet is strengthened while the storm track and eddy-driven jet are shifted equatorward. Conversely, if the subtropical SST gradients are reduced and the midlatitude gradients increased, then the storm track shows a strong poleward shift and a well-separated eddy-driven jet is produced. The sign of the SST anomaly is shown to play a secondary role in determining the overall tropospheric response. These findings are used to provide a new and consistent interpretation of some previous GCM studies concerning the atmospheric response to midlatitude SST anomalies.

Characteristics of Langmuir Turbulence in the Ocean Mixed Layer

This study uses large-eddy simulation (LES) to investigate the characteristics of Langmuir turbulence through the turbulent kinetic energy (TKE) budget. Based on an analysis of the TKE budget a velocity scale for Langmuir turbulence is proposed. The velocity scale depends on both the friction velocity and the surface Stokes drift associated with the wave field. The scaling leads to unique profiles of nondimensional dissipation rate and velocity component variances when the Stokes drift of the wave field is sufficiently large compared to the surface friction velocity. The existence of such a scaling shows that Langmuir turbulence can be considered as a turbulence regime in its own right, rather than a modification of shear-driven turbulence. Comparisons are made between the LES results and observations, but the lack of information concerning the wave field means these are mainly restricted to comparing profile shapes. The shapes of the LES profiles are consistent with observed profiles. The dissipation length scale for Langmuir turbulence is found to be similar to the dissipation length scale in the shear-driven boundary layer. Beyond this it is not possible to test the proposed scaling directly using available data. Entrainment at the base of the mixed layer is shown to be significantly enhanced over that due to normal shear turbulence.

Surface wave processes in air-sea interaction

We review briefly recent progress on understanding the role of surface waves on the marine atmospheric boundary layer and the ocean mixed layer and give a global perspective on these processes by analysing ERA-40 data. Ocean surface waves interact with the marine atmospheric boundary layer in two broad regimes: (i) the conventional wind-driven wave regime, when fast winds blow over slower moving waves, and (ii) a wave-driven wind regime when long wavelength swell propagates under low winds, and generates a wave-driven jet in the lower part of the marine boundary layer. Analysis of ERA-40 data indicates that the wave-driven wind regime is as prevalent as the conventional wind-driven regime. Ocean surface waves also change profoundly mixing in the ocean mixed layer through generation of Langmuir circulation. Results from large-eddy simulation are used here to develop a scaling for the resulting Langmuir turbulence, which is a necessary step in developing a parametrization of the process. ERA-40 data is then used to show that the Langmuir regime is the predominant regime over much of the global ocean, providing a compelling motivation for parameterising this process in ocean general circulation models.

U.K. HiGEM: The New U.K. High-Resolution Global Environment Model—Model Description and Basic Evaluation

This article describes the development and evaluation of the U.K.’s new High-Resolution Global Environmental Model (HiGEM), which is based on the latest climate configuration of the Met Office Unified Model, known as the Hadley Centre Global Environmental Model, version 1 (HadGEM1). In HiGEM, the horizontal resolution has been increased to 0.83° latitude × 1.25° longitude for the atmosphere, and 1/3° × 1/3° globally for the ocean. Multidecadal integrations of HiGEM, and the lower-resolution HadGEM, are used to explore the impact of resolution on the fidelity of climate simulations. Generally, SST errors are reduced in HiGEM. Cold SST errors associated with the path of the North Atlantic drift improve, and warm SST errors are reduced in upwelling stratocumulus regions where the simulation of low-level cloud is better at higher resolution. The ocean model in HiGEM allows ocean eddies to be partially resolved, which dramatically improves the representation of sea surface height variability. In the Southern Ocean, most of the heat transports in HiGEM is achieved by resolved eddy motions, which replaces the parameterized eddy heat transport in the lower-resolution model. HiGEM is also able to more realistically simulate small-scale features in the wind stress curl around islands and oceanic SST fronts, which may have implications for oceanic upwelling and ocean biology. Higher resolution in both the atmosphere and the ocean allows coupling to occur on small spatial scales. In particular, the small-scale interaction recently seen in satellite imagery between the atmosphere and tropical instability waves in the tropical Pacific Ocean is realistically captured in HiGEM. Tropical instability waves play a role in improving the simulation of the mean state of the tropical Pacific, which has important implications for climate variability. In particular, all aspects of the simulation of ENSO (spatial patterns, the time scales at which ENSO occurs, and global teleconnections) are much improved in HiGEM.

Estimating soil water fluxes from soil water records obtained using dielectric sensors

Accurate estimation of the soil water balance (SWB) is important for a number of applications (e.g. environmental, meteorological, agronomical and hydrological). The objective of this study was to develop and test techniques for the estimation of soil water fluxes and SWB components (particularly infiltration, evaporation and drainage below the root zone) from soil water records. The work presented here is based on profile soil moisture data measured using dielectric methods, at 30-min resolution, at an experimental site with different vegetation covers (barley, sunflower and bare soil). Estimates of infiltration were derived by assuming that observed gains in the soil profile water content during rainfall were due to infiltration. Inaccuracies related to diurnal fluctuations present in the dielectric-based soil water records are resolved by filtering the data with adequate threshold values. Inconsistencies caused by the redistribution of water after rain events were corrected by allowing for a redistribution period before computing water gains. Estimates of evaporation and drainage were derived from water losses above and below the deepest zero flux plane (ZFP), respectively. The evaporation estimates for the sunflower field were compared to evaporation data obtained with an eddy covariance (EC) system located elsewhere in the field. The EC estimate of total evaporation for the growing season was about 25% larger than that derived from the soil water records. This was consistent with differences in crop growth (based on direct measurements of biomass, and field mapping of vegetation using laser altimetry) between the EC footprint and the area of the field used for soil moisture monitoring. Copyright (c) 2007 John Wiley & Sons, Ltd.

Do High-Frequency Eddies Contribute to Low-Frequency Teleconnection Tendencies?*

An isentropic potential vorticity (PV) budget analysis is employed to examine the role of synoptic transients, advection, and nonconservative processes as forcings for the evolution of the low-frequency PV anomalies locally and those associated with the North Atlantic Oscillation (NAO) and the Pacific–North American (PNA) pattern. Specifically, the rate of change of the low-frequency PV is expressed as a sum of tendencies due to divergence of eddy transport, advection by the low-frequency flow (hereafter referred to as advection), and the residual nonconservative processes. The balance between the variances and covariances of these terms is illustrated using a novel vector representation. It is shown that for most locations, as well as for the PNA pattern, the PV variability is dominantly driven by advection. The eddy forcing explains a small amount of the tendency variance. For the NAO, the role of synoptic eddy fluxes is found to be stronger, explaining on average 15% of the NAO tendency variance. Previous studies have not assessed quantitively how the various forcings balance the tendency. Thus, such studies may have overestimated the role of eddy fluxes for the evolution of teleconnections by examining, for example, composites and regressions that indicate maintenance, rather than evolution driven by the eddies. The authors confirm this contrasting view by showing that during persistent blocking (negative NAO) episodes the eddy driving is relatively stronger.

The role of eddies in driving the tropospheric response to stratospheric heating perturbations

A simplified general circulation model has been used to investigate the chain of causality whereby changes in tropospheric circulation and temperature are produced in response to stratospheric heating perturbations. Spinup ensemble experiments have been performed to examine the evolution of the tropospheric circulation in response to such perturbations. The primary aim of these experiments is to investigate the possible mechanisms whereby a tropospheric response to changing solar activity over the 11-yr solar cycle could be produced in response to heating of the equatorial lower stratosphere. This study therefore focuses on a stratospheric heating perturbation in which the heating is largest in the tropics. For comparison, experiments are also performed in which the stratosphere is heated uniformly at all latitudes and in which it is heated preferentially in the polar region. Thus, the mechanisms discussed have a wider relevance for the impact of stratospheric perturbations on the troposphere. The results demonstrate the importance of changing eddy momentum fluxes in driving the tropospheric response. This is confirmed by the lack of a similar response in a zonally symmetric model with fixed eddy forcing. Furthermore, it is apparent that feedback between the tropospheric eddy fluxes and tropospheric circulation changes is required to produce the full model response. The quasigeostrophic index of refraction is used to diagnose the cause of the changes in eddy behavior. It is demonstrated that the latitudinal extent of stratospheric heating is important in determining the direction of displacement of the tropospheric jet and storm track.

Formation and circulation of the cold intermediate layer in the Gulf of Saint Lawrence

[ 1] The local heat content and formation rate of the cold intermediate layer (CIL) in the Gulf of Saint Lawrence are examined using a combination of new in situ wintertime observations and a three-dimensional numerical model. The field observations consist of five moorings located throughout the gulf over the period of November 2002 to June 2003. The observations demonstrate a substantially deeper surface mixed layer in the central and northeast gulf than in regions downstream of the buoyant surface outflow from the Saint Lawrence Estuary. The mixed-layer depth in the estuary remains shallow (< 60 m) throughout winter, with the arrival of a layer of near-freezing waters between 40 and 100 m depth in April. An eddy-permitting ice-ocean model with realistic forcing is used to hindcast the period of observation. The model simulates well the seasonal evolution of mixed-layer depth and CIL heat content. Although the greatest heat losses occur in the northeast, the most significant change in CIL heat content over winter occurs in the Anticosti Trough. The observed renewal of CIL in the estuary in spring is captured by the model. The simulation highlights the role of the northwest gulf, and in particular, the separation of the Gaspe Current, in controlling the exchange of CIL between the estuary and the gulf. In order to isolate the effects of inflow through the Strait of Belle Isle on the CIL heat content, we examine a sensitivity experiment in which the strait is closed. This simulation shows that the inflow has a less important effect on the CIL than was suggested by previous studies.

The cumulus-capped boundary layer. I: Modelling transports in the cloud layer

Scalar-flux budgets have been obtained from large-eddy simulations (LESs) of the cumulus-capped boundary layer. Parametrizations of the terms in the budgets are discussed, and two parametrizations for the transport term in the cloud layer are proposed. It is shown that these lead to two models for scalar transports by shallow cumulus convection. One is equivalent to the subsidence detrainment form of convective tendencies obtained from mass-flux parametrizations of cumulus convection. The second is a flux-gradient relationship that is similar in form to the non-local parametrizations of turbulent transports in the dry-convective boundary layer. Using the fluxes of liquid-water potential temperature and total water content from the LES, it is shown that both models are reasonable diagnostic relations between fluxes and the vertical gradients of the mean fields. The LESs used in this study are for steady-state convection and it is possible to treat the fluxes of conserved thermodynamic variables as independent, and ignore the effects of condensation. It is argued that a parametrization of cumulus transports in a model of the cumulus-capped boundary layer should also include an explicit representation of condensation. A simple parametrization of the liquid-water flux in terms of conserved variables is also derived.

The physical properties of the atmosphere in the new Hadley Centre Global Environmental Model (HadGEM1). Part II: Aspects of variability and regional climate

The performance of the atmospheric component of the new Hadley Centre Global Environmental Model (HadGEM1) is assessed in terms of its ability to represent a selection of key aspects of variability in the Tropics and extratropics. These include midlatitude storm tracks and blocking activity, synoptic variability over Europe, and the North Atlantic Oscillation together with tropical convection, the Madden-Julian oscillation, and the Asian summer monsoon. Comparisons with the previous model, the Third Hadley Centre Coupled Ocean-Atmosphere GCM (HadCM3), demonstrate that there has been a considerable increase in the transient eddy kinetic energy (EKE), bringing HadGEM1 into closer agreement with current reanalyses. This increase in EKE results from the increased horizontal resolution and, in combination with the improved physical parameterizations, leads to improvements in the representation of Northern Hemisphere storm tracks and blocking. The simulation of synoptic weather regimes over Europe is also greatly improved compared to HadCM3, again due to both increased resolution and other model developments. The variability of convection in the equatorial region is generally stronger and closer to observations than in HadCM3. There is, however, still limited convective variance coincident with several of the observed equatorial wave modes. Simulation of the Madden-Julian oscillation is improved in HadGEM1: both the activity and interannual variability are increased and the eastward propagation, although slower than observed, is much better simulated. While some aspects of the climatology of the Asian summer monsoon are improved in HadGEM1, the upper-level winds are too weak and the simulation of precipitation deteriorates. The dominant modes of monsoon interannual variability are similar in the two models, although in HadCM3 this is linked to SST forcing, while in HadGEM1 internal variability dominates. Overall, analysis of the phenomena considered here indicates that HadGEM1 performs well and, in many important respects, improves upon HadCM3. Together with the improved representation of the mean climate, this improvement in the simulation of atmospheric variability suggests that HadGEM1 provides a sound basis for future studies of climate and climate change.