TOPAZ4: an ocean-sea ice data assimilation system for the North Atlantic and Arctic

Abstract. We present a detailed description of TOPAZ4, the latest version of TOPAZ – a coupled ocean-sea ice data assimilation system for the North Atlantic Ocean and Arctic. It is the only operational, large-scale ocean data assimilation system that uses the ensemble Kalman filter. This means that TOPAZ features a time-evolving, state-dependent estimate of the state error covariance. Based on results from the pilot MyOcean reanalysis for 2003–2008, we demonstrate that TOPAZ4 produces a realistic estimate of the ocean circulation in the North Atlantic and the sea-ice variability in the Arctic. We find that the ensemble spread for temperature and sea-level remains fairly constant throughout the reanalysis demonstrating that the data assimilation system is robust to ensemble collapse. Moreover, the ensemble spread for ice concentration is well correlated with the actual errors. This indicates that the ensemble statistics provide reliable state-dependent error estimates – a feature that is unique to ensemble-based data assimilation systems. We demonstrate that the quality of the reanalysis changes when different sea surface temperature products are assimilated, or when in-situ profiles below the ice in the Arctic Ocean are assimilated. We find that data assimilation improves the match to independent observations compared to a free model. Improvements are particularly noticeable for ice thickness, salinity in the Arctic, and temperature in the Fram Strait, but not for transport estimates or underwater temperature. At the same time, the pilot reanalysis has revealed several flaws in the system that have degraded its performance. Finally, we show that a simple bias estimation scheme can effectively detect the seasonal or constant bias in temperature and sea-level.

[1]  D. Dumont,et al.  Modeling the Dynamics of the North Water Polynya Ice Bridge , 2009 .

[2]  V. Canuto,et al.  Ocean Turbulence. Part II: Vertical Diffusivities of Momentum, Heat, Salt, Mass, and Passive Scalars , 2002 .

[3]  Ron Kwok,et al.  Decline in Arctic sea ice thickness from submarine and ICESat records: 1958–2008 , 2009 .

[4]  Peter R. Oke,et al.  Representation Error of Oceanic Observations for Data Assimilation , 2008 .

[5]  R. Kwok,et al.  Evaluation of Arctic sea ice thickness simulated by Arctic Ocean Model Intercomparison Project models , 2012 .

[6]  S. Zalesak.  Fully multidimensional flux-corrected transport algorithms for fluids , 1979 .

[7]  Peter R. Oke,et al.  A Modeling Study of the Three-Dimensional Continental Shelf Circulation off Oregon. Part I: Model-Data Comparisons , 2002 .

[8]  Pierre Brasseur,et al.  A demonstration of ensemble-based assimilation methods with a layered OGCM from the perspective of operational ocean forecasting systems , 2003 .

[9]  T. Oki,et al.  Design of Total Runoff Integrating Pathways (TRIP)—A Global River Channel Network , 1998 .

[10]  Peter R. Oke,et al.  Ocean data assimilation: a case for ensemble optimal interpolation , 2010 .

[11]  Michael Steele,et al.  PHC: A Global Ocean Hydrography with a High-Quality Arctic Ocean , 2001 .

[12]  G. Halliwell,et al.  Evaluation of vertical coordinate and vertical mixing algorithms in the HYbrid-Coordinate Ocean Model (HYCOM) , 2004 .

[13]  R. Bleck,et al.  Global Ocean Prediction with the Hybrid Coordinate Ocean Model , 2009 .

[14]  Laurent Bertino,et al.  The TOPAZ monitoring and prediction system for the Atlantic and Arctic Oceans , 2008 .

[15]  P. L. Houtekamer,et al.  Ensemble Kalman filtering , 2005 .

[16]  J. Thepaut,et al.  The ERA‐Interim reanalysis: configuration and performance of the data assimilation system , 2011 .

[17]  M. Iskandarani,et al.  Comparison of advection schemes for high-order h–p finite element and finite volume methods , 2005 .

[18]  Françoise Ogor,et al.  ERS‐1/2 orbit improvement using TOPEX/POSEIDON: The 2 cm challenge , 1998 .

[19]  Jean-Marc Molines,et al.  Evaluation of high-resolution sea ice models on the basis of statistical and scaling properties of Arctic sea ice drift and deformation , 2009 .

[20]  S. Cohn,et al.  Ooce Note Series on Global Modeling and Data Assimilation Construction of Correlation Functions in Two and Three Dimensions and Convolution Covariance Functions , 2022 .

[21]  G. Evensen,et al.  Asynchronous data assimilation with the EnKF , 2010 .

[22]  C. Donlon,et al.  OSTIA : An operational, high resolution, real time, global sea surface temperature analysis system , 2007, OCEANS 2007 - Europe.

[23]  R. Kwok,et al.  Arctic ice-ocean simulation with optimized model parameters: Approach and assessment , 2011 .

[24]  Michael Karcher,et al.  A Synthesis of Exchanges Through the Main Oceanic Gateways to the Arctic Ocean , 2011 .

[25]  Timothy P. Boyer,et al.  World Ocean Atlas 2005 Volume 1: Temperature [+DVD] , 2006 .

[26]  Pavel Sakov,et al.  Relation between two common localisation methods for the EnKF , 2011 .

[27]  Peter R. Oke,et al.  Ocean Data Assimilation Systems for GODAE , 2009 .

[28]  Rainer Bleck,et al.  An oceanic general circulation model framed in hybrid isopycnic-Cartesian coordinates , 2002 .

[29]  Ryan D. Torn,et al.  Performance Characteristics of a Pseudo-Operational Ensemble Kalman Filter , 2008 .

[30]  Geir Evensen,et al.  Coordinate Transformation on a Sphere Using Conformal Mapping , 1999 .

[31]  G. Evensen,et al.  Analysis Scheme in the Ensemble Kalman Filter , 1998 .

[32]  Ola M. Johannessen,et al.  Norwegian Remote Sensing Experiment: Evaluation of the Nimbus 7 scanning multichannel microwave radiometer for sea ice research , 1983 .

[33]  Peter A. Rochford,et al.  Efficient and Accurate Bulk Parameterizations of Air-Sea Fluxes for Use in General Circulation Models , 2000 .

[34]  Geir Evensen,et al.  Assimilation of ice concentration in a coupled ice–ocean model, using the Ensemble Kalman filter , 2003 .

[35]  Craig M. Lee,et al.  Closing the loop – Approaches to monitoring the state of the Arctic Mediterranean during the International Polar Year 2007–2008 , 2011 .

[36]  S. Levitus,et al.  Warming of the World Ocean , 2000 .

[37]  Geir Evensen,et al.  The Ensemble Kalman Filter: theoretical formulation and practical implementation , 2003 .

[38]  G. Bigg,et al.  Impact of flow through the Canadian Archipelago and Bering Strait on the North Atlantic and Arctic circulation: An ocean modelling study , 2002 .

[39]  Peter R. Oke,et al.  The Bluelink ocean data assimilation system (BODAS) , 2008 .

[40]  Timothy P. Boyer,et al.  Warming of the world ocean, 1955–2003 , 2005 .

[41]  William S. Curran,et al.  A/I: a synthesis , 1982, ACM-SE 20.

[42]  S. Gorshkov,et al.  World ocean atlas , 1976 .

[43]  P. Oke,et al.  Implications of the Form of the Ensemble Transformation in the Ensemble Square Root Filters , 2008 .

[44]  G. Evensen.  Sequential data assimilation with a nonlinear quasi‐geostrophic model using Monte Carlo methods to forecast error statistics , 1994 .

[45]  P. Jones,et al.  The Twentieth Century Reanalysis Project , 2009 .

[46]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[47]  Remy Baraille,et al.  The HYCOM (HYbrid Coordinate Ocean Model) data assimilative system , 2007 .

[48]  James Cummings,et al.  A comparison of sequential assimilation schemes for ocean prediction with the HYbrid Coordinate Ocean Model (HYCOM): Twin experiments with static forecast error covariances , 2011 .

[49]  Erik Andersson,et al.  Influence‐matrix diagnostic of a data assimilation system , 2004 .

[50]  Jens Schröter,et al.  Data assimilation for marine monitoring and prediction: The MERCATOR operational assimilation systems and the MERSEA developments , 2005 .

[51]  Massimo Bonavita,et al.  The ensemble Kalman filter in an operational regional NWP system: preliminary results with real observations , 2008 .

[52]  Peter R. Oke,et al.  A deterministic formulation of the ensemble Kalman filter : an alternative to ensemble square root filters , 2008 .

[53]  W. Teague,et al.  A Comparison Between the Generalized Digital Environmental Model and Levitus climatologies , 1990 .

[54]  J. Morison,et al.  Internal wave dissipation under sea ice , 1985 .

[55]  W. Budgell,et al.  Mass and heat transports in the NE Barents Sea: Observations and models , 2009 .

[56]  Thomas M. Smith,et al.  Improved Global Sea Surface Temperature Analyses Using Optimum Interpolation , 1994 .

[57]  J. Dorandeu,et al.  Effects of Global Mean Atmospheric Pressure Variations on Mean Sea Level Changes from TOPEX/Poseidon , 1999 .

[58]  H. Loeng,et al.  The seasonal cycle in the Atlantic transport to the Barents Sea during the years 1997–2001 , 2004 .

[59]  M. Maqueda,et al.  Sensitivity of a global sea ice model to the treatment of ice thermodynamics and dynamics , 1997 .

[60]  J. Johannessen,et al.  Steric height variability in the Nordic Seas , 2007 .

[61]  Ron Kwok,et al.  Improved modeling of the Arctic halocline with a subgrid‐scale brine rejection parameterization , 2009 .

[62]  W. Large,et al.  Open Ocean Momentum Flux Measurements in Moderate to Strong Winds , 1981 .

[63]  K. A. Orvik,et al.  Volume and Heat Transports to the Arctic Ocean Via the Norwegian and Barents Seas , 2008 .

[64]  E. Hunke,et al.  An Elastic–Viscous–Plastic Model for Sea Ice Dynamics , 1996 .

[65]  Istvan Szunyogh,et al.  Efficient data assimilation for spatiotemporal chaos: A local ensemble transform Kalman filter , 2005, physics/0511236.

[66]  Craig H. Bishop,et al.  Adaptive sampling with the ensemble transform Kalman filter , 2001 .

[67]  E. Fahrbach,et al.  Variation of Measured Heat Flow Through the Fram Strait Between 1997 and 2006 , 2008 .

[68]  R. Kwok.  Outflow of Arctic Ocean Sea Ice into the Greenland and Barents Seas: 1979-2007 , 2009 .

[69]  Timothy P. Boyer,et al.  World ocean atlas 2013. Volume 2, Salinity , 2002 .

[70]  J. Whitaker,et al.  The Twentieth Century Reanalysis Project 3 , 2011 .