Data assimilation for initialization of seasonal forecasts

This article reviews the requirements for a data assimilation system from the perspective of initializing seasonal forecasts. It provides a historical perspective of the developments in ocean data assimilation and ocean observing systems. It also discusses the differences between state estimation and initialization, and presents a brief assessment of different initialization strategies. The value of assimilating ocean data to estimate the ocean state and to initialize seasonal forecasts is demonstrated. However, it is also shown that the assumption of unbiased models in conventional data assimilation methods is not suitable for the production of long temporal records of ocean initial states. This is due to the combined effect of model-forcing error and the changing nature of the observing system. Bias correction algorithms are therefore important in the estimation of long records of ocean states. In the equatorial ocean, the delicate balance between the mass and the velocity fields should be preserved in order to maintain realistic circulations. The most common approach for initializing seasonal forecasts is the so-called full uncoupled initialization, which basically consists of producing an ocean reanalysis by assimilating ocean observations into an ocean model driven by atmospheric fluxes. Alternative approaches are the so-called anomaly initialization, which only attempts to initialize the anomalous state without any attempt of correcting mean; the latter is usually conducted in coupled mode, but coupled and anomaly initialization are not synonymous, and there are approaches where the initialization of the full state is done in coupled mode. The relative value of the approaches is system dependent, but as a long-term strategy the full initialization in coupled mode is more promising

Historical reconstruction of the Atlantic Meridional Overturning Circulation from the ECMWF operational ocean reanalysis

A reconstruction of the Atlantic Meridional Overturning Circulation (MOC) for the period 1959–2006 has been derived from the ECMWF operational ocean reanalysis. The reconstruction shows a wide range of time-variability, including a downward trend. At 26N, both the MOC intensity and changes in its vertical structure are in good agreement with previous estimates based on trans-Atlantic surveys. At 50N, the MOC and strength of the subpolar gyre are correlated at interannual time scales, but show opposite secular trends. Heat transport variability is highly correlated with the MOC but shows a smaller trend due to the warming of the upper ocean, which partially compensates for the weakening of the circulation. Results from sensitivity experiments show that although the time-varying upper boundary forcing provides useful MOC information, the sequential assimilation of ocean data further improves the MOC estimation by increasing both the mean and the time variability

A Drift-Free Decadal Climate Prediction System for the Community Earth System Model

Performance of a newly developed decadal climate prediction system is examined using the low-resolution Community Earth System Model (CESM). To identify key sources of predictability and determine the role of upper and deeper ocean data assimilation, we first conduct a series of perfect model experiments. These experiments reveal the importance of upper ocean temperature and salinity assimilation in reducing sea surface temperature biases. However, to reduce biases in the sea surface height, data assimilation below 300 m in the ocean is necessary, in particular for high-latitude regions. The perfect model experiments clearly emphasize the key role of combined three-dimensional ocean temperature and salinity assimilation in reproducing mean state and model trajectories. Applying this knowledge to the realistic decadal climate prediction system, we conducted an ensemble of ocean assimilation simulations with the fully coupled CESM covering the period 1960–2014. In this system, we assimilate three-dimensional ocean temperature and salinity data into the ocean component of CESM. Instead of assimilating direct observations, we assimilate temperature and salinity anomalies obtained from the ECMWF Ocean Reanalysis version 4 (ORA-S4). Anomalies are calculated relative to the sum of the ORA-S4 climatology and an estimate of the externally forced signal. As a result of applying the balanced ocean conditions to the model, our hindcasts show only very little drift and initialization shocks. This new prediction system exhibits multiyear predictive skills for decadal climate variations of the Atlantic meridional overturning circulation (AMOC) and North Pacific decadal variability

Predicting El Niño in 2014 and 2015.

Early in 2014 several forecast systems were suggesting a strong 1997/98-like El Niño event for the following northern hemisphere winter 2014/15. However the eventual outcome was a modest warming. In contrast, winter 2015/16 saw one of the strongest El Niño events on record. Here we assess the ability of two operational seasonal prediction systems to forecast these events, using the forecast ensembles to try to understand the reasons underlying the very different development and outcomes for these two years. We test three hypotheses. First we find that the continuation of neutral ENSO conditions in 2014 is associated with the maintenance of the observed cold southeast Pacific sea surface temperature anomaly; secondly that, in our forecasts at least, warm west equatorial Pacific sea surface temperature anomalies do not appear to hinder El Niño development; and finally that stronger westerly wind burst activity in 2015 compared to 2014 is a key difference between the two years. Interestingly, in these years at least, this interannual variability in wind burst activity is predictable. ECMWF System 4 tends to produce more westerly wind bursts than Met Office GloSea5 and this likely contributes to the larger SST anomalies predicted in this model in both years

Derive observable ocean climate indicators from seasonal forecast

Derive the user-relevant indicators defined in Milestone MS7 from the ensemble of ECMWF and CMCC seasonal forecasts systems contributing to C3

Skill assessment of ECV/EOV from seasonal forecast

Assess the seasonal forecast skill of selected ocean variables - SST, OHC300m, and SSH - from the ensemble of ECMWF and CMCC seasonal forecasts systems contributing to C3

Comparative Analysis of Upper Ocean Heat Content Variability from Ensemble Operational Ocean Analyses

Upper ocean heat content (HC) is one of the key indicators of climate variability on many time-scales extending from seasonal to interannual to long-term climate trends. For example, HC in the tropical Pacific provides information on thermocline anomalies that is critical for the longlead forecast skill of ENSO. Since HC variability is also associated with SST variability, a better understanding and monitoring of HC variability can help us understand and forecast SST variability associated with ENSO and other modes such as Indian Ocean Dipole (IOD), Pacific Decadal Oscillation (PDO), Tropical Atlantic Variability (TAV) and Atlantic Multidecadal Oscillation (AMO). An accurate ocean initialization of HC anomalies in coupled climate models could also contribute to skill in decadal climate prediction. Errors, and/or uncertainties, in the estimation of HC variability can be affected by many factors including uncertainties in surface forcings, ocean model biases, and deficiencies in data assimilation schemes. Changes in observing systems can also leave an imprint on the estimated variability. The availability of multiple operational ocean analyses (ORA) that are routinely produced by operational and research centers around the world provides an opportunity to assess uncertainties in HC analyses, to help identify gaps in observing systems as they impact the quality of ORAs and therefore climate model forecasts. A comparison of ORAs also gives an opportunity to identify deficiencies in data assimilation schemes, and can be used as a basis for development of real-time multi-model ensemble HC monitoring products. The OceanObs09 Conference called for an intercomparison of ORAs and use of ORAs for global ocean monitoring. As a follow up, we intercompared HC variations from ten ORAs -- two objective analyses based on in-situ data only and eight model analyses based on ocean data assimilation systems. The mean, annual cycle, interannual variability and longterm trend of HC have been analyze

Skills of the user-relevant ocean indicators

This document demonstrates the capability of seasonal forecasting systems to predict observable and user-relevant ocean climate indicators

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

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