A-TSCV – ESA project

Background

Measurements of the powerful, complex and highly variable ocean surface currents and surface waves are fundamental to our understanding of ocean circulation and air—sea interaction. This is because they influence the Earth system at time scales from wind waves, weather through to climate (GCOS, 2015). The motivation for better knowledge and understanding of ocean surface currents has its foundation in traditional maritime activities serving society at all levels including: shipping, maritime safety, marine operations, increasing maritime activities with sea-ice, fisheries, renewable energy, pollution events, environmental management, resource exploitation, ports and harbour operations, recreation, numerical weather prediction, ocean forecasting, and climate monitoring, amongst others. The actual velocity of a water parcel in contact with the atmosphere at any given location and time (i.e. the surface ocean current) is called the Total Surface Current Velocity (TSCV). TSCV is intimately related to waves (and thus wind stress) that through the wave dispersion relationship, reveal the fingerprint of TSCV.

A fundamental interface of the Earth system is between the ocean and the atmosphere through which all exchange of heat, momentum, gas and mass must be accomplished. On either side of this interface boundary layers are present, in which friction and mixing processes dominate where motion is not primarily horizontal, but also involves vertical exchange. In the ocean, this boundary layer typically has a dimension of the order of 10 m, characterised by a vertically homogeneous slab that moves horizontally (Price et al., 1986). However, the actively mixed layer is variable due to diurnal and seasonal stratification (Fig. 2.11), freshwater from rain and land run-off, and can be shallower than 1 m. The TSCV is characterised by the co-existence of different flow patterns and phenomena that have a physical expression on the sea surface. These expressions are embedded in the surface roughness of the ocean surface that can be measured.

The TSCV, denoted as vector U, has a markedly different impact on long-term drift behaviour compared to its usually defined components, the geostrophic current Ug, derived from altimetry, the mean wind-driven current UEk and the extremely important wave-induced Stokes drift, US (OOPC 2017). The long-term average velocity of particles at the ocean surface is well described by the sum of the three terms, U≈Ug+UEk+US. Each of the three terms plays an important and specific role. In the presence of waves, it is useful to separate the TSCV in a quasi-Eulerian current (as measured by a fixed current-meter located below the wave troughs) and the Stokes drift, the average motion due to waves alone (Jenkins, 1989). The vertical profile of the quasi-Eulerian current, its magnitude and surface angle, will depart from the simplified theory of Ekman (1905). In particular, the surface angle can range between 45 and 80°, depending on wind speed and stratification. As a result, the drift of floating material at the surface has very different patterns when one considers only geostrophic currents, with very weak convergence, Ekman currents, Stokes drift or the combination of the three (Fig. 1).

Figure 1. Numerical simulation of the distributions of surface particles (for example representing marine plastic debris) after 10 years, starting from a uniform distribution. The four panels show the effect of three components of the surface current combined, or their effect taken separately. (Adapted from Onink et al., 2019)

The Stokes drift has a direct contribution to the TSCV, but also a very strong influence on the other component by influencing mixing (Ardhuin and Jenkins 2006, Belcher et al., 2012). It is the ageostrophic motions (i.e. those not in balance with the pressure gradients), in this case UEk+US, that are important for surface convergence, and thus the aggregation of floating material such as plastic litter, but also define cross-shelf exchanges, including the export of carbon from land to the deep ocean. Owing to the properties of the mixing layer, the TSCV may be representative of a layer that varies in thickness between 0.10 m and tens of meters depending on stratification, and it has been very difficult to measure routinely. Stratification of the ocean refers to its density separation into layers owing to the dependence of the density on temperature, salinity and pressure. Stratification can lead to effective de-coupling of the upper-ocean water from that at depth (e.g. during strong diurnal heating or shallow halocline
conditions leading to a ‘slippery’ shallow layer where water at the surface can flow freely over that at depth (see Kudryavtsev and Soloviev, 1990). Only the measurements in the top few metres of the ocean, such as provided by un-drogued and shallow-drogued drifters (e.g. Lumpkin et al., 2017) that follow the surface-water motion (Novelli et al., 2017) or coastal HF-radars, can be representative of the surface mixed layer in most conditions, providing an estimate of the TSCV.

Alternative measurement techniques, at a fixed position or with drifting instruments drogued below the depth where the wave-induced Stokes drift US becomes negligible, measure the Eulerian mean flow UE, i.e. the total flow U-US (Jenkins, 1989, Niiler et al., 1995). By measuring the sea-level gradients, altimeter measurements miss the difference between the geostrophic current that often dominates below the mixed layer: the difference between the geostrophic current that often dominates below the mixed layer and the TSCV is 25 cm s-1 in the daytime.

The TSCV can be formally defined as the Lagrangian mean velocity at the sea surface (Andrews and McIntyre, 1978): the TSCV is the actual velocity of water parcels that are right at the ocean surface in contact with the atmosphere. The depth of the surface layer representative of the TSCV in the upper ocean depends on its vertical stratification. It includes a multiscale continuum of variability across all space and time scales governed by quasi-random forcing perturbations and nonlinear interactions at a given time and location. Motion at the ocean surface is the result of a superposition of different forces acting on the same ocean surface and those connected to the deeper ocean layers. Only in recent years have numerical ocean models started to include the necessary large range of space and time scales to represent these and the reader is referred to https://bit.ly/2XsokAe for one example realisation (D. Menemenlis personal communication, see Torres et al., 2018 for analyses). Thus, the processes resolved today in state-of-the-art models are only part of the full range of motions that contribute to the upper-ocean velocity field (e.g. Fox-Kemper et al., 2019).

Today, circulation models are capable of simulating ocean currents in response to atmospheric wind, tides, and buoyancy forcing. A question remains: How well do such simulations represent reality?

Their full validity in dynamical regimes readily accessible through potential future satellites missions (e.g. EE-9 candidate mission SKIM, EE-10 candidate mission Harmony), has never been tested except at a few locations (e.g. Scott et al., 2010) or using the sparse global Surface Velocity Program drifter array (Elipot et al., 2016). A synthesis of the Copernicus Marine Environmental Monitoring System (CMEMS) model performance on surface current was made by Rémy et al. (2019). For the mean error, the modelled zonal current is typically within 5 cm s-1 of the SVP drifter data, but up to 40 cm s-1 along the equator for the western half of the Pacific. Random errors typically lead to 60 km distances of three-day drift end-points, corresponding to a constant 23 cm s-1 error over 3 days.

Providing high-resolution maps of the surface current will reveal details that will help test theories and models, involving especially the wind-driven currents and wave-driven drift and mixing. On the other hand, model results will also be essential for interpreting such measurements that are being sampled irregularly for some of the relevant time scales, with ‘entangled’ contributions of rapid physical processes in snapshots measured at each satellite pass. Artificial intelligence or machine learning can be used to identify patterns and statistical descriptions embedded in datasets with applications to Earth-system science that will certainly emerge in unforeseen areas. Essential aspects here involve: response of the upper ocean to wind forcing, the vertical structure of resulting flow fields, the dissipation of wind energy, resulting mixing and transports and their convergences.

Of particular importance is the revolution offered by fully coupled ocean–atmosphere model systems from climate to weather forecasting time scales (e.g. Williams et al., 2019). It is the TSCV and wave spectrum that are most relevant to such systems because wind–wave–current interaction governs, to first order, all exchanges of heat, gas, mass and momentum between the atmosphere and the ocean. A particular new aspect of modelling involve Stokes drifts and the impact on momentum and matter transports (Belcher et al., 2012). While global NOP traditionally were not coupled to surface-wave models, great efforts have been made recently towards a consistent ocean–wave–atmosphere coupling (Harris 2018, Lemarié et al., 2015, 2019, Beljaars et al., 2017), also including sea ice (Boutin et al., 2019). Respective parameterisations show success with respect to improving tracer fields and transport characteristics. However, only very limited tests could be done owing to the lack of relevant observations. For the first time in the history of oceanography and ocean modelling, missions such as SKIM/Harmony could supply these observations and thus will allow for further improvement of existing parameterisation. This is expected to lead to a new modelling framework for ocean and coupled climate models.

Models are usually used alone or in combination with observations to study ocean or environmental processes and mechanisms as well as to predict changes. However, they can also be improved by combining them with the observations in a formal sense through data assimilation. Today, significant incremental progress has been made in ocean forecasting (e.g. Schiller et al., 2015) with the inclusion of tides, improved surface forcing/surface fields and waves/current interactions and many models stand ready to assimilate wave-current data such as that from SKIM/Harmony (e.g. Tonani et al., 2015, Chassignet and Sandery, 2013, Villas Bôas et al., 2019). Satellite data are being assimilated on a regular basis, including sea surface temperature, ocean colour and, particularly, altimetry – providing a dynamical pressure boundary condition. SKIM/Harmony data would be entirely complementary by providing a surface layer horizontal velocity field. Impacts have to be determined, though, as this type of data were never assimilated before. To assimilate Stokes drift data requires further developments before the full potential can be inferred. We envision that along with the next generation of ocean models containing active Stokes drift modules, also the assimilation capabilities of these models will be extended SKIM/Harmony data will be an essential driver and benefit to such an evolution.

Given the acute responsibility of teams working in extremely challenging situations such as oil, chemical spill or harmful algal bloom response, marine search and rescue activities (where TSCV estimates are a critical input for success) access to regular repeat-coverage TSCV is essential. Figure 2 shows the progress envisaged 30 years ago in climate models for integrating a dedicated ocean-wave model (Hasselmann, 1990), which is now a reality for some of the climate models used in the IPCC AR6.

Figure 2. Left climate model integration over the past 30 years to (right) the development of operational integrated ocean modelling system (CMEMS). The steady increase in the resolution of models requires a strong dynamic coupling of ocean and atmosphere, via the TSCV for the momentum and mechanical energy equations. (After K. Hasselmann, University of Hamburg, and CMEMS)

The same type of effort, after combining biogeochemistry and sea ice with ocean circulation, is taking place in the Copernicus Marine Environment Monitoring Service. In that context, given the higher resolution (CMEMS will operate a 3-km global ocean model by 2025), a proper representation of surface layer processes is much more critical (e.g. Renault et al., 2017, 2019). Ocean—atmosphere coupling, including TSCV, is one of the major research and development efforts both for CMEMS and ECMWF because it is critical for both the oceanic and atmospheric dynamics. Data assimilation for coupled systems is one particular area of intense research, and TSCV data, when available, would be a key variable to control the state of the coupled model. New measurements of TSCV SKIM/Harmony will certainly accelerate and contribute to these endeavours in a significant capacity.

Satellite TSCV measurements are not yet available to the ocean modelling community since no satellite system has been implemented. In the next few years the NASA/CNES SWOT mission will be launched and the EE-9 SKIM, EE-10 Harmony and NASA WaCM concepts are being studied in detail. However, The technical challenges of assimilating satellite TSCV and the challenging task of assessing the quality of satellite TSCV measurements using limited in situ reference data have not been developed to a satisfactory level. Furthermore, assessing the impact of the assimilation using objective tools and reporting requires significant effort.

Objectives

The GODAE OceanView Science Team (GOVST) convened the Observing System Evaluation Task Team (OS-Eval TT) to evaluate the impact of different measurement systems. This is achieved by designing, implementing and reporting the results of specific Observing System Experiments (OSE). Each OSE produces an Observation Impact Statement Report (OISR) which describes in detail the end-to-end process of the OSE and the results obtained. OISR were proposed to the OS-Eval TT by ESA at the 3rd GOVST meeting held at ESA HQ, Paris November 2011. Through dedicated OSE experiments, GOVST is able to inform Space Agencies on the utility and performance of satellite data products and formulate specific requirements for ocean observations on the basis of improved understanding of data utility.

This A-TSCV study is focused on the design, implementation and reporting of an Observing System Evaluation of satellite Doppler-derived waves and TSCV using synthetic waves and TSCV measurements (e.g. from the SKIM L2 SKIMulator system). The primary output will be a GOV-ST Observation Impact Statement Report focused on satellite waves and TSCV, journal publications and a workshop dedicated to the findings and approach taken by the study team. In addition, it is anticipated that the study will include the development of a better understanding of suitable observation operators required by the assimilation system(s) and a more complete understanding of how to utilise Doppler/Along-Track Interferometry (ATI) derived waves and TSCV in ocean modelling systems.

Against this background, the aim of the A-TSCV study is:

The Objectives for the A-TSCV study are:

Team

The consortium is composed of four companies: Met Office, Mercator Océan International, OceanDataLab and OceanNext. The Met Office, the prime contractor and sole interface to ESA, will be responsible for WP1, WP4 and WP6. Mercator Océan International (sub-contractor) will be responsible for WP2 and WP5. OceanDataLab (sub-contractor) will be responsible for WP3. OceanNext (sub-contractor) will also be involved in WP3. In addition to the main consortium, we also include two sets of scientific advisors from Ifremer and Thales Alenia Space who will not be funded in the project.

Met Office

Dr. Matt Martin (M) is the project manager. He is a Met Office Science Fellow and the manager of the Marine Data Assimilation group within the Ocean Forecasting R&D team at the Met Office. He is co-chair of the Data Assimilation Task Team in OceanPredict.

Dr. Andrew Saulter (M) will contribute to WPs 2, 3 and 4. He is the manager of the Met Office Surge, Waves and Met Ocean Projects team. He has developed the Met Office’s global and regional wave and storm surge forecasting systems for many years and is a member of the CMEMS waves working group.

Dr. Jennifer Waters (F) is assigned WP leader for WP4. She is a senior scientist in the Marine Data Assimilation Group within the Ocean Forecasting R&D team at the Met Office. She has developed the Met Office’s global ocean data assimilation capabilities over a number of years, prior to which she completed a PhD in the assimilation of wave data from HF radar.

Dr. Robert King (M) is to work on WP4.2. He is a senior scientist in the Marine Data Assimilation group within the Ocean Forecasting R&D team at the Met Office. He has been involved in various international projects focussed on assessing the impact of observations in operational ocean-forecasting systems, including E-AIMS, ERA-Clim2, AtlantOS and ESA SMOS-NINO15.

Kirsten Wilmer-Becker (F) is assigned WP leader for WP1.2. She is the OceanPredict Programme Office Coordinator, responsible for coordinating OceanPredict communication and activities, organising OceanPredict events, developing and maintaining the OceanPredict website and managing the OceanPredict budget.

Mercator Océan international

Dr. Stéphane Law Shune (M) is working in the “Forecasting and analysis systems evolution” team of Mercator Ocean on several projects related to waves. Stéphane has a B.Sc. in Applied Physics, a M.Sc. in Meteorology/Oceanography and a thesis in Physical Oceanography.  His area of expertise covers (1) surface layer dynamics and ocean-wave interactions in numerical models, (2) implementation and validation of wave reanalysis and (3) surface current assessment for Lagrangian dynamics.

Dr. Jean-Michel Lellouche (M) is at the head of the “Forecasting and analysis systems evolution” team, within the Operational Oceanography department. He is responsible of the preparation of the future real time forecasting systems. He has a perfect knowledge of the systems configuration and often run short OSEs to test the operational systems.

Dr. Elisabeth Rémy (F) is at the head of the “Observations for ocean analysis and forecasting systems” team within the Department. She has a long experience in ocean modelling, data assimilation and more recently observation impact studies. She coordinates activities on the benefit of future observation deployment for the Mercator Océan operational ocean analysis and forecasts. She is currently involved in the EuroSea H2020 project, SWOT science team and GODAE OSE/OSSE task team.

CapGemini DEMS

Dr Isabelle Mirouze (F) is a consultant for Mercator Océan International and will mainly contribute to WP4. She has a long experience in ocean data assimilation on different systems. She had taken part in various projects working on error covariance estimate and modelling, observation operators, OSSEs, and atmosphere-ocean coupled data assimilation.

OceanDataLab

Dr. Lucile Gaultier (F) is a Research Engineer, working on improving the reconstruction of ocean surface current using various remote sensing data as well as on the use of the synergy of data and Lagrangian diagnostics to improve our understanding of the ocean surface dynamics. She is also involved in SKIM proposal to simulate skim-like data and evaluate the ocean surface current products performance.

Dr. Fabrice Collard (M) is president of OceanDataLab IFREMER spinoff at CERSAT, working on ocean remote sensing multi-sensor synergy methods and tools with a special focus on ocean surface wind, wave and currents. He has a long experience in ocean surface current measurements from his post-poc on HF radars to the development of SAR Doppler techniques (ENVISAT and now Sentinel-1) but also drifters development/deployments and analysis (DRIFT4SKIM).  He participated in several projects related to upper ocean circulation such as ESA Globcurrent project, member of ESA SKIM Mission advisory group. ESA radial velocity assessment from Sentinel-1 Doppler and platform Gyros, EU TOPVOYS for ship routing optimisation.

Datlas

Dr. Clément Ubelmann (M) is a research scientist now at Datlas. He holds a PhD in Oceanography with a strong expertise on data assimilation. His main interests are the inversion of Ocean surface satellite data to reconstruct mesoscale and sub-mesoscale fields, involving data assimilation or data-driven techniques. He is strongly involved in the SWOT mission to develop new algorithms for calibration and mapping the future high-resolution altimetry. Firstly as a postdoc and research scientist at JPL (2009-2011) and after at CLS (2011-2019). He is going to join Ocean Next (march 2020). During the recent years, he also focused on surface current reconstruction including high-frequency ageostrophic processes. He contributed to the SKIM-phase A , designing total surface current reconstruction algorithms combining Doppler and Altimetry.

  Schedule and milestones        Work Logic        Work packages        Deliverables  

Schedule and milestones

The proposed schedule and main milestones of the project are given below. The project will run for two years, divided into six phases:

• Task 1: Outreach, Communication and Promotion
• Task 2: OSE design
• Task 3: Data preparation and collection
• Task 4: OSE implementation
• Task 5: OSE output processing, validation and analysis
• Task 6: International workshop and final reporting

The A-TSCV web portal will be implemented during the first 3 months of the project. The A-TSCV database will be ready at T0+12. WPs 2, 3, 4 and 5 are mainly sequential with some overlap allowing the next task to start before the end of the current task. Task 6 will start at T0+18 and mainly consists in the organisation of the international workshop.

  Back top top 

Work logic

  Back top top 

Work package descriptions

WP1:   Project management & communication (Met Office)

• Project management (Met Office)
• Communication (Met Office)

WP2:  OSE design (MOI)

• TSCV data assimilation design (Met Office)
• OSE design (MOI)

WP3:  Data collection (OceanDataLab)

• SKIM data generation (OceanDataLab)
• Standard observations generation (OceanDataLab)

WP4:  OSE implementation (Met Office)

• Met Office OSE implementation (Met Office)
• MOI OSE implementation (MOI)

WP5:  Data processing and analysis (MOi)

• Writing papers (MOI)
• OISR (Met Office)

WP6:  International workshop and final reporting (Met Office)

• Final reporting (Met Office)
• International workshop (MOi)

  Back top top 

Deliverables

Web-stories (WS)

The project activities are described through web-stories which are published at least twice per year with input from project partners. The web-stories are  linked to a collection of high-quality graphics produced by the project.

WS-1: Developing assimilation of TSCV data

Idealised observation experiments and estimation of velocity forecast error covariances
The ocean total surface current velocity (TSCV) is the Lagrangian mean velocity at the    instantaneous sea surface (Marié et al., 2020). Accurate forecasting of the ocean TSCV is important for applications such as search and rescue, tracking marine plastic and for coupled ocean/atmosphere/sea-ice/wave forecasting…..  read more or download as pdf.

WS-2: Simulating an altimetry observing system from a global high resolution model

A recent experiment to test the impact of a future spaceborne Doppler surface current observing system in combination with an altimetry observing system.
High-resolution Ocean-General-Circulation-Models (OGCMs) provide realistic scenes of ocean circulation which are extremely useful to build, test and evaluate the impact of an observing system for data assimilation. This is the common technique known as Observing System Simulation Experiment (OSSE) where a non-assimilated model (a “free run”) is considered as a ground-truth used to generate realistic – but synthetic – observations….  read more or download as pdf.

WS-3: Assimilation of TSCV data in a global ocean forecasting system

An experiment assessing the impact of assimilating synthetic TSCV data in the Met Office global ocean forecasting system. The study demonstrates that assimilation of TSCV data has the potential to significantly improve our prediction of the ocean state.
Observing System Simulation Experiments (OSSEs) are used to assess the impact of assimilating potential Total Surface Current Velocities (TSCV) in the Met Office global ocean forecasting system (FOAM; Blockley et al., 2014). Synthetic observations are generated from a Nature Run (NR) which is a high resolution model run. We generate synthetic observations for all standard observation types …..  read more or download as pdf.

WS-4: TSCV observations assimilated in an ocean global configuration

Experiments have been run to assess the impact of assimilating TSCV data in the Mercator Océan International forecasting system. Preliminary results show a nice improvement for the surface and sub-surface currents.

The impact of assimilating TSCV data in the Mercator Océan International forecasting system (SAM) is assessed through Observing System Simulation Experiments (OSSEs). Two experiments are run with the same settings and configuration, assimilating the same synthetic observation networks simulated from a high-resolution free run (Nature Run; NR). Only one of the experiments (ASM) assimilates synthetic TSCV observations. Its comparison to the Control experiment (CTL) allows us to highlight the impact of such data…..  read more or download as pdf.

High-resolution graphics

Graphics produced during the project and linked to the web-stories can be downloaded from this page.

Graphics from WS-1: “Idealised observation experiments and estimation of velocity forecast error covariances”. 

Click on image to view enlargement and/or download the file.

Surface increments for speed, temperature, salinity and SSH. From MOI (top) and FOAM (bottom).

Zonally averaged horizontal forecast error correlation length scales for unbalanced surface U and V. These are estimated by fitting a Gaussian function with two correlation scales to the NMC error covariance data.  Black and blue lines are length scales in the x-direction and y-direction, respectively. Dashed and solid lines are the long and short scale, respectively.

U vertical forecast error correlations with the surface. Plot (a) shows the global mean correlations plotted against a normalising depth. For the green, blue and black line the normalising quantity is the global mean MldRho, MldZ and Ekman depth respectively. The horizontal red line shows where the normalised depth is 1 and the vertical red line is the value of a Gaussian function when the depth variable equals the correlation length scale. The shaded region shows the standard deviation of the error correlations. Plot (b) shows a latitudinal section of the zonal mean vertical correlations with the surface. The green, blue and black lines are the zonal mean MldRho, MldZ and Ekman depth, respectively.

Graphics from WS-2: “Simulating an altimetry observing system from a global high resolution model”. 

Click on image to view enlargement and/or download the file.

Scheme of the Observing System Simulation experiment

Snapshots of the 6h difference SSH in the LLC-4320 simulation for theSSH (upper panel), SSH minus the estimated tides (middle panel) and SSH minus the estimated tides and DAC (bottom panel)?

Graphics from WS-3: “Assimilation of TSCV data in a global ocean forecasting system”. 

Click on image to view enlargement and/or download the file.

Global RMSE (solid line) and Mean Error (dashed line) for T, S, U and V profiles and SSH time series calculate between the 1^st of February and 31^st of July 2009. Blue is the control experiment, green is the A-TSCV experiment.

Difference in the surface zonal velocity RMSE for the control and A-TSCV experiments calculated for July 2009. Blue shows regions where the A-TSCV experiment had a lower RMSE than the control while red shows regions where the A-TSCV experiment has a higher RMSE.

Global surface zonal velocity RMSE (mean) for the control and A-TSCV experiments as a function of forecast lead time. The statistics were calculated by running 7 day forecasts every 7 days between February and December 2009 to produce a total of 47 forecasts. Blue is the control experiment, green is the A-TSCV experiment.

A-TSCV Twitter page

Tweets by a_tscv

  Reference Documents        Reports        Presentations        Publications       Data 

Reference documents

Project management area (accessible only by project members)

A-TSCV Workshop –  13 June 2023

The A-TSCV project is organising a workshop on the results of the project activities. Please visit the workshop pages here: https://oceanpredict.org/events/a-tscv-workshop/

This workshop aims at sharing experience in assimilating ocean surface velocity data. Today, HF radar data are assimilated in operational systems. Future space missions should observe total surface currents. Lesson learnt from past experience and design experiments to assess the impact of future observations will be discussed as well as assimilation methods well suited for those observations. Current status and challenges for efficiently assimilating surface ocean current observations from diverse observing platforms will be identified, in research and operational context. Observation gaps and requirements for operational systems will be discussed.

General topics to be addressed in this workshop include:

• Assimilation of direct velocity observations from HF radar, drifting platforms, ADCP, future space missions, SAR measurements, …
• Assimilation methods for surface currents, including Lagrangian and Eulerian approaches,
• Surface current estimation via data assimilation from coastal to global scale.