Search

Swath altimetry is likely to revolutionize our ability to monitor and forecast ocean dynamics. To meet the requirements of the EU Copernicus Marine Service, a constellation of two wide-swath altimeters is envisioned for the long-term (post-2030) evolution of the Copernicus Sentinel 3 topography mission. A series of observing system simulation experiments (OSSEs) is carried out to quantify the expected performances. The OSSEs use a state-of-the-art high-resolution (1/12∘) global ocean data assimilation system similar to the one used operationally by the Copernicus Marine Service. Flying a constellation of two wide-swath altimeters will provide a major improvement of our capabilities to monitor and forecast the oceans. Compared to the present situation with three nadir altimeters flying simultaneously, the sea surface height (SSH) analysis and 7 d forecast error are globally reduced by about 50 % in the OSSEs. With two wide-swath altimeters, the quality of SSH 7 d forecasts is equivalent to the quality of SSH analysis errors from three nadir altimeters. Our understanding of ocean currents is also greatly improved (30 % improvements at the surface and 50 % at 300 m depth). The resolution capabilities will be drastically improved and will be closer to 100 km wavelength compared to about 250 km today. Flying a constellation of two wide-swath altimeters thus looks to be a very promising solution for the long-term evolution of the Sentinel 3 constellation and the Copernicus Marine Service.

Canada has the longest coastline in the world and includes diverse ocean environments, from the frozen waters of the Canadian Arctic Archipelago to the confluence region of Labrador and Gulf Stream waters on the east coast. There is a strong need for a pan-Canadian operational regional ocean prediction capacity covering all Canadian coastal areas in support of marine activities including emergency response, search and rescue, and safe navigation in ice-infested waters. Here we present the first pan-Canadian operational regional ocean analysis system developed as part of the Regional Ice Ocean Prediction System version 2 (RIOPSv2) running in operations at the Canadian Centre for Meteorological and Environmental Prediction (CCMEP). The RIOPSv2 domain extends from 26∘ N in the Atlantic Ocean through the Arctic Ocean to 44∘ N in the Pacific Ocean, with a model grid resolution that varies between 3 and 8 km. RIOPSv2 includes a multivariate data assimilation system based on a reduced-order extended Kalman filter together with a 3D-Var bias correction system for water mass properties. The analysis system assimilates satellite observations of sea level anomaly and sea surface temperature, as well as in situ temperature and salinity measurements. Background model error is specified in terms of seasonally varying model anomalies from a 10-year forced model integration, allowing inhomogeneous anisotropic multivariate error covariances. A novel online tidal harmonic analysis method is introduced that uses a sliding-window approach to reduce numerical costs and allow for the time-varying harmonic constants necessary in seasonally ice-infested waters. Compared to the Global Ice Ocean Prediction System (GIOPS) running at CCMEP, RIOPSv2 also includes a spatial filtering of model fields as part of the observation operator for sea surface temperature (SST). In addition to the tidal harmonic analysis, the observation operator for sea level anomaly (SLA) is also modified to remove the inverse barometer effect due to the application of atmospheric pressure forcing fields. RIOPSv2 is compared to GIOPS and shown to provide similar innovation statistics over a 3-year evaluation period. Specific improvements are found near the Gulf Stream for all model fields due to the higher model grid resolution, with smaller root mean squared (rms) innovations for RIOPSv2 of about 5 cm for SLA and 0.5 ∘C for SST. Verification against along-track satellite observations demonstrates the improved representation of mesoscale features in RIOPSv2 compared to GIOPS, with increased correlations of SLA (0.83 compared to 0.73) and reduced rms differences (12 cm compared to 14 cm). While the RIOPSv2 grid resolution is 3 times higher than GIOPS, the power spectral density of surface kinetic energy provides an indication that the effective resolution of RIOPSv2 is roughly double that of the global system (35 km compared to 66 km). Observations made as part of the Year of Polar Prediction (2017–2019) provide a rare glimpse at errors in Arctic water mass properties and show average salinity biases over the upper 500 m of 0.3–0.4 psu in the eastern Beaufort Sea in RIOPSv2.

The impact of forthcoming wide-swath altimetry missions on the ocean analysis and forecasting system was investigated by means of OSSEs (observing system simulation experiments). These experiments were performed with a regional data assimilation system, implemented in the Iberian–Biscay–Ireland (IBI) region, at 1∕12° resolution using simulated observations derived from a fully eddy-resolving free simulation at 1∕36° resolution over the same region. The objective of the experiments was to assess the ability of different satellite constellations to constrain the ocean analyses and forecasts, considering both along-track altimeters and future wide-swath missions; consequently, the capability of the data assimilation techniques used in the Mercator Ocean operational system to effectively combine the different kinds of measurements was also investigated. These assessments were carried out as part of a European Space Agency (ESA) study on the potential role of wide-swath altimetry in future versions of the European Union Copernicus programme. The impact of future wide-swath altimetry data is evident for investigating the reliability of sea level values in OSSEs. The most significant results were obtained when looking at the sensitivity of the system to wide-swath instrumental error: considering a constellation of three nadir and two accurate (small instrumental error) wide-swath altimeters, the error in ocean analysis was reduced by up to 50 % compared to conventional altimeters. Investigating the impact of the repetitivity of the future measurements, the results showed that two wide-swath missions had a major impact on sea-level forecasting – increasing the accuracy over the entire time window of the 5-day forecasts – compared with a single wide-swath instrument. A spectral analysis underlined that the contributions of wide-swath altimetry data observed in ocean analyses and forecast statistics were mainly due to the more accurate resolution, compared with along-track data, of ocean variability at spatial scales smaller than 100 km. Considering the ocean currents, the results confirmed that the information provided by wide-swath measurements at the surface is propagated down the water column and has a considerable impact (30 %) on ocean currents (up to a depth of 300 m), compared with the present constellation of altimeters. The ocean analysis and forecasting systems used here are those currently used by the Copernicus Marine Environment and Monitoring Service (CMEMS) to provide operational services and ocean reanalysis. The results obtained in the OSSEs considering along-track altimeters were consistent with those derived from real data (observing system experiments, OSEs). OSSEs can also be used to assess the potential of new observing systems, and in this study the results showed that future constellations of altimeters will have a major impact on constraining the CMEMS ocean analysis and forecasting systems and their applications.

Since 19 October 2016, and in the framework of Copernicus Marine Environment Monitoring Service (CMEMS), Mercator Ocean has delivered real-time daily services (weekly analyses and daily 10-day forecasts) with a new global 1∕12° high-resolution (eddy-resolving) monitoring and forecasting system. The model component is the NEMO platform driven at the surface by the IFS ECMWF atmospheric analyses and forecasts. Observations are assimilated by means of a reduced-order Kalman filter with a three-dimensional multivariate modal decomposition of the background error. Along-track altimeter data, satellite sea surface temperature, sea ice concentration, and in situ temperature and salinity vertical profiles are jointly assimilated to estimate the initial conditions for numerical ocean forecasting. A 3D-VAR scheme provides a correction for the slowly evolving large-scale biases in temperature and salinity.This paper describes the recent updates applied to the system and discusses the importance of fine tuning an ocean monitoring and forecasting system. It details more particularly the impact of the initialization, the correction of precipitation, the assimilation of climatological temperature and salinity in the deep ocean, the construction of the background error covariance and the adaptive tuning of observation error on increasing the realism of the analysis and forecasts.The scientific assessment of the ocean estimations are illustrated with diagnostics over some particular years, assorted with time series over the time period 2007–2016. The overall impact of the integration of all updates on the product quality is also discussed, highlighting a gain in performance and reliability of the current global monitoring and forecasting system compared to its previous version.

Since December 2010, the MyOcean global analysis and forecasting system has consisted of the Mercator Océan NEMO global 1/4° configuration with a 1/12° nested model over the Atlantic and the Mediterranean. The open boundary data for the nested configuration come from the global 1/4° configuration at 20° S and 80° N. The data are assimilated by means of a reduced-order Kalman filter with a 3-D multivariate modal decomposition of the forecast error. It includes an adaptive-error estimate and a localization algorithm. A 3-D-Var scheme provides a correction for the slowly evolving large-scale biases in temperature and salinity. Altimeter data, satellite sea surface temperature and in situ temperature and salinity vertical profiles are jointly assimilated to estimate the initial conditions for numerical ocean forecasting. In addition to the quality control performed by data producers, the system carries out a proper quality control on temperature and salinity vertical profiles in order to minimise the risk of erroneous observed profiles being assimilated in the model. This paper describes the recent systems used by Mercator Océan and the validation procedure applied to current MyOcean systems as well as systems under development. The paper shows how refinements or adjustments to the system during the validation procedure affect its quality. Additionally, we show that quality checks (in situ, drifters) and data sources (satellite sea surface temperature) have as great an impact as the system design (model physics and assimilation parameters). The results of the scientific assessment are illustrated with diagnostics over the year 2010 mainly, assorted with time series over the 2007–2011 period. The validation procedure demonstrates the accuracy of MyOcean global products, whose quality is stable over time. All monitoring systems are close to altimetric observations with a forecast RMS difference of 7 cm. The update of the mean dynamic topography corrects local biases in the Indonesian Throughflow and in the western tropical Pacific. This improves also the subsurface currents at the Equator. The global systems give an accurate description of water masses almost everywhere. Between 0 and 500 m, departures from in situ observations rarely exceed 1 °C and 0.2 psu. The assimilation of an improved sea surface temperature product aims to better represent the sea ice concentration and the sea ice edge. The systems under development are still suffering from a drift which can only be detected by means of a 5-yr hindcast, preventing us from upgrading them in real time. This emphasizes the need to pursue research while building future systems for MyOcean2 forecasting.