Application of data assimilation to ocean and climate prediction

Ocean prediction systems are now able to analyse and predict temperature, salinity and velocity structures within the ocean by assimilating measurements of the ocean’s temperature and salinity into physically based ocean models. Data assimilation combines current estimates of state variables, such as temperature and salinity, from a computational model with measurements of the ocean and atmosphere in order to improve forecasts and reduce uncertainty in the forecast accuracy. Data assimilation generally works well with ocean models away from the equator but has been found to induce vigorous and unrealistic overturning circulations near the equator. A pressure correction method was developed at the University of Reading and the Met Office to control these circulations using ideas from control theory and an understanding of equatorial dynamics. The method has been used for the last 10 years in seasonal forecasting and ocean prediction systems at the Met Office and European Center for Medium-range Weather Forecasting (ECMWF). It has been an important element in recent re-analyses of the ocean heat uptake that mitigates climate change.

An assessment of air-sea heat fluxes from ocean and coupled reanalyses

Sixteen monthly air–sea heat flux products from global ocean/coupled reanalyses are compared over 1993–2009 as part of the Ocean Reanalysis Intercomparison Project (ORA-IP). Objectives include assessing the global heat closure, the consistency of temporal variability, comparison with other flux products, and documenting errors against in situ flux measurements at a number of OceanSITES moorings. The ensemble of 16 ORA-IP flux estimates has a global positive bias over 1993–2009 of 4.2 ± 1.1 W m−2. Residual heat gain (i.e., surface flux + assimilation increments) is reduced to a small positive imbalance (typically, +1–2 W m−2). This compensation between surface fluxes and assimilation increments is concentrated in the upper 100 m. Implied steady meridional heat transports also improve by including assimilation sources, except near the equator. The ensemble spread in surface heat fluxes is dominated by turbulent fluxes (>40 W m−2 over the western boundary currents). The mean seasonal cycle is highly consistent, with variability between products mostly <10 W m−2. The interannual variability has consistent signal-to-noise ratio (~2) throughout the equatorial Pacific, reflecting ENSO variability. Comparisons at tropical buoy sites (10°S–15°N) over 2007–2009 showed too little ocean heat gain (i.e., flux into the ocean) in ORA-IP (up to 1/3 smaller than buoy measurements) primarily due to latent heat flux errors in ORA-IP. Comparisons with the Stratus buoy (20°S, 85°W) over a longer period, 2001–2009, also show the ORA-IP ensemble has 16 W m−2 smaller net heat gain, nearly all of which is due to too much latent cooling caused by differences in surface winds imposed in ORA-IP.

ENSO and greenhouse warming

The El Niño/Southern Oscillation (ENSO) is the dominant climate phenomenon affecting extreme weather conditions worldwide. Its response to greenhouse warming has challenged scientists for decades, despite model agreement on projected changes in mean state. Recent studies have provided new insights into the elusive links between changes in ENSO and in the mean state of the Pacific climate. The projected slow-down in Walker circulation is expected to weaken equatorial Pacific Ocean currents, boosting the occurrences of eastward-propagating warm surface anomalies that characterize observed extreme El Niño events. Accelerated equatorial Pacific warming, particularly in the east, is expected to induce extreme rainfall in the eastern equatorial Pacific and extreme equatorward swings of the Pacific convergence zones, both of which are features of extreme El Niño. The frequency of extreme La Niña is also expected to increase in response to more extreme El Niños, an accelerated maritime continent warming and surface-intensified ocean warming. ENSO-related catastrophic weather events are thus likely to occur more frequently with unabated greenhouse-gas emissions. But model biases and recent observed strengthening of the Walker circulation highlight the need for further testing as new models, observations and insights become available.

Effect of surface restoring on subsurface variability in a climate model during 1949-2005

Initializing the ocean for decadal predictability studies is a challenge, as it requires reconstructing the little observed subsurface trajectory of ocean variability. In this study we explore to what extent surface nudging using well-observed sea surface temperature (SST) can reconstruct the deeper ocean variations for the 1949–2005 period. An ensemble made with a nudged version of the IPSLCM5A model and compared to ocean reanalyses and reconstructed datasets. The SST is restored to observations using a physically-based relaxation coefficient, in contrast to earlier studies, which use a much larger value. The assessment is restricted to the regions where the ocean reanalyses agree, i.e. in the upper 500 m of the ocean, although this can be latitude and basin dependent. Significant reconstruction of the subsurface is achieved in specific regions, namely region of subduction in the subtropical Atlantic, below the thermocline in the equatorial Pacific and, in some cases, in the North Atlantic deep convection regions. Beyond the mean correlations, ocean integrals are used to explore the time evolution of the correlation over 20-year windows. Classical fixed depth heat content diagnostics do not exhibit any significant reconstruction between the different existing observation-based references and can therefore not be used to assess global average time-varying correlations in the nudged simulations. Using the physically based average temperature above an isotherm (14 °C) alleviates this issue in the tropics and subtropics and shows significant reconstruction of these quantities in the nudged simulations for several decades. This skill is attributed to the wind stress reconstruction in the tropics, as already demonstrated in a perfect model study using the same model. Thus, we also show here the robustness of this result in an historical and observational context.

The impact of westerly wind bursts and ocean initial state on the development, and diversity of El Niño events

Westerly wind bursts (WWBs) that occur in the western tropical Pacific are believed to play an important role in the development of El Niño events. Here, following the study of Lengaigne et al. (Clim Dyn 23(6):601–620, 2004), we conduct numerical simulations in which we reexamine the response of the climate system to an observed wind burst added to a coupled general circulation model. Two sets of twin ensemble experiments are conducted (each set has control and perturbed experiments). In the first set, the initial ocean heat content of the system is higher than the model climatology (recharged), while in the second set it is nearly normal (neutral). For the recharged state, in the absence of WWBs, a moderate El Niño with a maximum warming in the central Pacific (CP) develops in about a year. In contrast, for the neutral state, there develops a weak La Niña. However, when the WWB is imposed, the situation dramatically changes: the recharged state slides into an El Niño with a maximum warming in the eastern Pacific, while the neutral set produces a weak CP El Niño instead of previous La Niña conditions. The different response of the system to the exact same perturbations is controlled by the initial state of the ocean and the subsequent ocean–atmosphere interactions involving the interplay between the eastward shift of the warm pool and the warming of the eastern equatorial Pacific. Consequently, the observed diversity of El Niño, including the occurrence of extreme events, may depend on stochastic atmospheric processes, modulating El Niño properties within a broad continuum.

Increasing frequency of extreme El Niño events due to greenhouse warming

El Niño events are a prominent feature of climate variability with global climatic impacts. The 1997/98 episode, often referred to as ‘the climate event of the twentieth century’1, 2, and the 1982/83 extreme El Niño3, featured a pronounced eastward extension of the west Pacific warm pool and development of atmospheric convection, and hence a huge rainfall increase, in the usually cold and dry equatorial eastern Pacific. Such a massive reorganization of atmospheric convection, which we define as an extreme El Niño, severely disrupted global weather patterns, affecting ecosystems4, 5, agriculture6, tropical cyclones, drought, bushfires, floods and other extreme weather events worldwide3, 7, 8, 9. Potential future changes in such extreme El Niño occurrences could have profound socio-economic consequences. Here we present climate modelling evidence for a doubling in the occurrences in the future in response to greenhouse warming. We estimate the change by aggregating results from climate models in the Coupled Model Intercomparison Project phases 3 (CMIP3; ref. 10) and 5 (CMIP5; ref. 11) multi-model databases, and a perturbed physics ensemble12. The increased frequency arises from a projected surface warming over the eastern equatorial Pacific that occurs faster than in the surrounding ocean waters13, 14, facilitating more occurrences of atmospheric convection in the eastern equatorial region.

Interdecadal changes on the seasonal prediction of the western North Pacific summer climate around the late 1970s and early 1990s

Identifying predictability and the corresponding sources for the western North Pacific (WNP) summer climate in the case of non-stationary teleconnections during recent decades benefits for further improvements of long-range prediction on the WNP and East Asian summers. In the past few decades, pronounced increases on the summer sea surface temperature (SST) and associated interannual variability are observed over the tropical Indian Ocean and eastern Pacific around the late 1970s and over the Maritime Continent and western–central Pacific around the early 1990s. These increases are associated with significant enhancements of the interannual variability for the lower-tropospheric wind over the WNP. In this study, we further assess interdecadal changes on the seasonal prediction of the WNP summer anomalies, using May-start retrospective forecasts from the ENSEMBLES multi-model project in the period 1960–2005. It is found that prediction of the WNP summer anomalies exhibits an interdecadal shift with higher prediction skills since the late 1970s, particularly after the early 1990s. Improvements of the prediction skills for SSTs after the late 1970s are mainly found around tropical Indian Ocean and the WNP. The better prediction of the WNP after the late 1970s may arise mainly from the improvement of the SST prediction around the tropical eastern Indian Ocean. The close teleconnections between the tropical eastern Indian Ocean and WNP summer variability work both in the model predictions and observations. After the early 1990s, on the other hand, the improvements are detected mainly around the South China Sea and Philippines for the lower-tropospheric zonal wind and precipitation anomalies, associating with a better description of the SST anomalies around the Maritime Continent. A dipole SST pattern over the Maritime Continent and the central equatorial Pacific Ocean is closely related to the WNP summer anomalies after the early 1990s. This teleconnection mode is quite predictable, which is realistically reproduced by the models, presenting more predictable signals to the WNP summer climate after the early 1990s.

Understanding advances in the simulation of the Madden-Julian Oscillation in a numerical weather prediction model.

The Madden-Julian Oscillation (MJO) is the dominant mode of intraseasonal variability in the Trop- ics. It can be characterised as a planetary-scale coupling between the atmospheric circulation and organised deep convection that propagates east through the equatorial Indo-Pacific region. The MJO interacts with weather and climate systems on a near-global scale and is a crucial source of predictability for weather forecasts on medium to seasonal timescales. Despite its global signifi- cance, accurately representing the MJO in numerical weather prediction (NWP) and climate models remains a challenge. This thesis focuses on the representation of the MJO in the Integrated Forecasting System (IFS) at the European Centre for Medium-Range Weather Forecasting (ECMWF), a state-of-the-art NWP model. Recent modifications to the model physics in Cycle 32r3 (Cy32r3) of the IFS led to ad- vances in the simulation of the MJO; for the first time the observed amplitude of the MJO was maintained throughout the integration period. A set of hindcast experiments, which differ only in their formulation of convection, have been performed between May 2008 and April 2009 to asses the sensitivity of MJO simulation in the IFS to the Cy32r3 convective parameterization. Unique to this thesis is the attribution of the advances in MJO simulation in Cy32r3 to the mod- ified convective parameterization, specifically, the relative-humidity-dependent formulation for or- ganised deep entrainment. Increasing the sensitivity of the deep convection scheme to environmen- tal moisture is shown to modify the relationship between precipitation and moisture in the model. Through dry-air entrainment, convective plumes ascending in low-humidity environments terminate lower in the atmosphere. As a result, there is an increase in the occurrence of cumulus congestus, which acts to moisten the mid-troposphere. Due to the modified precipitation-moisture relationship more moisture is able to build up which effectively preconditions the tropical atmosphere for the transition to deep convection. Results from this thesis suggest that a tropospheric moisture control on convection is key to simulating the interaction between the physics and large-scale circulation associated with the MJO.

Eddies and variability in the Mozambique Channel

Between 1995 and 2000, on average 4 eddies per year are observed from satellite altimetry to propagate southward through the Mozambique Channel, into the upstream Agulhas region. Further south, these eddies have been found to control the timing and frequenc yof Agulhas ring shedding. Within the Mozambique Channel, anomalous SSH amplitudes rise to 30 cm ; in agreement with in situ measured velocities. Comparison of an observed velocit ysection with GCM model results shows that the Mozambique Channel eddies in these models are too surface intensified. Also, the number of eddies formed in the models is in disagreement with our observational analysis. Moored current meter measurements observing the passage of three eddies in 2000 are extended to a 5-year time series b yreferencing the anomalous surface currents estimated from altimeter data to a s ynoptic LADCP velocit y measurement. The results show intermittent edd ypassage at the mooring location. A statistical analysis of SSH observations in different parts of the Mozambique Channel shows a southward decrease of the dominant frequency of the variability, going from 7 per year in the extension of the South Equatorial Current north of Madagascar to 4 per year south of Madagascar. The observations suggest that frequency reduction is related to the Rossb ywaves coming in from the east