Assessment of the CMEMS global biogeochemical forecasting operational system, with assimilation of Ocean Colour data

Author(s):  
Julien Lamouroux ◽  
Alexandre Mignot ◽  
Coralie Perruche ◽  
Giovanni Ruggiero ◽  
Julien Paul ◽  
...  

<p>The operational production of data-assimilated biogeochemical state of the ocean is one of the challenging core projects of the Copernicus Marine Environment Monitoring Service (hereafter CMEMS). In that framework, Mercator Ocean International is in charge of improving the realism of its global 1⁄4° coupled physical-biogeochemical simulations, analyses and re‐analyses, and to develop an effective capacity to routinely estimate the biogeochemical state of the ocean, including, amongst others, the implementation of biogeochemical data assimilation. Primary objectives are to enhance the time representation of the seasonal cycle in the real time and reanalysis systems, and to provide a better control of the production in the equatorial regions.<br><br>In that framework, Mercator Ocean International has successfully updated its global biogeochemical analysis and forecasting system with an Ocean Color data assimilation capability. In this system, successfully commissioned in September 2019, the biogeochemical model (NEMO/PISCES) is offline coupled with the dynamical ocean (1/12° coarsened to 1/4° resolution) from the CMEMS global physical analysis and forecasting operational system (PSY4), at a daily frequency, and benefits from the assimilation of satellite (SSH-SST-SIC) and in situ physical data. Nevertheless, a biogeochemical climatological damping is activated in the biogeochemical model in order to mitigate the impact of some misconstrained processes (vertical velocities) from this physical data-assimilated forcing (under investigation). The dedicated assimilation of biogeochemical data relies on a simplified version of the SEEK filter, where the forecast error covariances are built from a fixed-basis - but seasonally variable - ensemble of anomalies computed from a multi-year numerical experiment (without biogeochemical data assimilation) with respect to a running mean. Regarding Ocean Colour observations, the system relies, as a first step, on the CMEMS Global Ocean surface chlorophyll concentration products, delivered in NRT.<br><br>The objective of this presentation is thus to provide (1) a short description of the implementation of the aforementioned data assimilation methodology in the forecasting system; (2) a synthesis of the assessment of this global biogeochemical forecasting system, by cross-comparing the assimilated solution with various datasets, both spatial (Ocean Colour) and in situ (BGC-Argo, GLODAP), and (3) a synthetic overview of the impact/benefit of the assimilation of the Ocean Colour data.</p>

2015 ◽  
Vol 12 (3) ◽  
pp. 1145-1186 ◽  
Author(s):  
V. Turpin ◽  
E. Remy ◽  
P. Y. Le Traon

Abstract. Observing System Experiments (OSEs) are carried out over a one-year period to quantify the impact of Argo observations on the Mercator-Ocean 1/4° global ocean analysis and forecasting system. The reference simulation assimilates sea surface temperature (SST), SSALTO/DUACS altimeter data and Argo and other in situ observations from the Coriolis data center. Two other simulations are carried out where all Argo and half of Argo data sets are withheld. Assimilating Argo observations has a significant impact on analyzed and forecast temperature and salinity fields at different depths. Without Argo data assimilation, large errors occur in analyzed fields as estimated from the differences when compared with in situ observations. For example, in the 0–300 m layer RMS differences between analyzed fields and observations reach 0.25 psu and 1.25 °C in the western boundary currents and 0.1 psu and 0.75 °C in the open ocean. The impact of the Argo data in reducing observation-model forecast error is also significant from the surface down to a depth of 2000 m. Differences between independent observations and forecast fields are thus reduced by 20 % in the upper layers and by up to 40 % at a depth of 2000 m when Argo data are assimilated. At depth, the most impacted regions in the global ocean are the Mediterranean outflow and the Labrador Sea. A significant degradation can be observed when only half of the data are assimilated. All Argo observations thus matter, even with a 1/4° model resolution. The main conclusion is that the performance of global data assimilation systems is heavily dependent on the availability of Argo data.


2012 ◽  
Vol 9 (2) ◽  
pp. 687-744 ◽  
Author(s):  
D. A. Ford ◽  
K. P. Edwards ◽  
D. Lea ◽  
R. M. Barciela ◽  
M. J. Martin ◽  
...  

Abstract. As part of the GlobColour project, daily chlorophyll-a observations, derived using remotely sensed ocean colour data from the MERIS, MODIS and SeaWiFS sensors, are produced. The ability of these products to be assimilated into a pre-operational global coupled physical-biogeochemical model has been tested, on both a hindcast and near-real-time basis, and the impact on the system assessed. The assimilation was found to immediately and significantly improve the bias, root mean square error and correlation of modelled surface chlorophyll concentration compared to the GlobColour observations, an improvement which was sustained throughout the year and in every ocean basin. Errors against independent in situ chlorophyll observations were also reduced, both at and beneath the ocean surface. However the model fit to in situ observations was not consistently better than that of climatology, due to errors in the underlying model. The assimilation scheme used is multivariate, updating all biogeochemical model state variables at all depths. Consistent changes were found in the other model variables, with reduced errors against in situ observations of nitrate and pCO2, and evidence of improved representation of zooplankton concentration. Annual mean surface fields of nutrients, alkalinity and carbon variables remained of similar quality compared to climatology. The near-real-time GlobColour products were found to be sufficiently reliable for operational purposes, and of benefit to both operational-style systems and reanalyses.


2014 ◽  
Vol 11 (4) ◽  
pp. 5399-5441 ◽  
Author(s):  
L. Visinelli ◽  
S. Masina ◽  
M. Vichi ◽  
A. Storto

Abstract. Prognostic simulations of ocean carbon distribution are largely dependent on an adequate representation of physical dynamics. In this work we show that the assimilation of temperature and salinity in a coupled ocean-biogeochemical model significantly improves the reconstruction of the carbonate system variables over the last two decades. For this purpose, we use the NEMO ocean global circulation model, coupled to the Biogeochemical Flux Model (BFM) in the global PELAGOS configuration. The assimilation of temperature and salinity is included into the coupled ocean-biogeochemical model by using a variational assimilation method. The use of ocean physics data assimilation improves the simulation of alkalinity and dissolved organic carbon against the control run as assessed by comparing with independent time series and gridded datasets. At the global scale, the effects of the assimilation of physical variables in the simulation of pCO2 improves the seasonal cycle in all basins, getting closer to the SOCAT estimates. Biases in the partial pressure of CO2 with respect to data that are evident in the control run are reduced once the physical data assimilation is used. The root mean squared errors in the pCO2 are reduced by up to 30% depending on the ocean basin considered. In addition, we quantify the relative contribution of biological carbon uptake on surface pCO2 by performing another simulation in which biology is neglected in the assimilated run.


2016 ◽  
Vol 46 (3) ◽  
pp. 713-729 ◽  
Author(s):  
Weiwei Fu

AbstractA three-dimensional variational data assimilation (3DVAR) method is implemented in a coupled physical–biogeochemical (CPB) model in the Baltic Sea. This study carries out a 10-yr assimilation experiment with satellite sea surface temperature (SST) and observed in situ temperature (T) and salinity (S) profiles. The impact of the assimilation is assessed with the focus on how the biogeochemical model responds to the improved hydrodynamics. The assimilation of temperature and salinity data yields considerable improvements in the physical model. On a basin scale, the mean bias of SST, T, S, and mixed layer depth (MLD) is decreased by 0.18°C (57%), 0.31°C (49%), 0.34 psu (43%), and 1.8 m (43%), respectively. More importantly, the biogeochemical simulation is improved in response to the physical data assimilation. Compared with in situ observations, the mean biases of chlorophyll a (Chl), dissolved inorganic nitrogen (DIN) and phosphorus (DIP) are decreased by 0.09 mg m−3 (15.5%), 0.19 mmol m−3 (9%), and 0.15 mmol m−3 (23%). Physical data assimilation also improves the simulated variability of Chl, DIN, and DIP and their correlations with observation. Compared with satellite observations, the mean bias of surface chlorophyll is reduced by 0.10–0.32 mg m−3 especially in the Skagerrak–Kattegat area and Bornholm basin. The decrease of total Chl change is caused by different mechanisms for winter and summer. While the deepened mixed layer acts as a dilution factor in winter, strengthened stratification agrees well with the decrease of chlorophyll in summer. In the vertical, relatively large changes of DIN and DIP occur below 60 m, which corresponds to the mean permanent halocline depth (~60–80 m) of the Baltic Sea.


Ocean Science ◽  
2016 ◽  
Vol 12 (1) ◽  
pp. 257-274 ◽  
Author(s):  
V. Turpin ◽  
E. Remy ◽  
P. Y. Le Traon

Abstract. Observing system experiments (OSEs) are carried out over a 1-year period to quantify the impact of Argo observations on the Mercator Ocean 0.25° global ocean analysis and forecasting system. The reference simulation assimilates sea surface temperature (SST), SSALTO/DUACS (Segment Sol multi-missions dALTimetrie, d'orbitographie et de localisation précise/Data unification and Altimeter combination system) altimeter data and Argo and other in situ observations from the Coriolis data center. Two other simulations are carried out where all Argo and half of the Argo data are withheld. Assimilating Argo observations has a significant impact on analyzed and forecast temperature and salinity fields at different depths. Without Argo data assimilation, large errors occur in analyzed fields as estimated from the differences when compared with in situ observations. For example, in the 0–300 m layer RMS (root mean square) differences between analyzed fields and observations reach 0.25 psu and 1.25 °C in the western boundary currents and 0.1 psu and 0.75 °C in the open ocean. The impact of the Argo data in reducing observation–model forecast differences is also significant from the surface down to a depth of 2000 m. Differences between in situ observations and forecast fields are thus reduced by 20 % in the upper layers and by up to 40 % at a depth of 2000 m when Argo data are assimilated. At depth, the most impacted regions in the global ocean are the Mediterranean outflow, the Gulf Stream region and the Labrador Sea. A significant degradation can be observed when only half of the data are assimilated. Therefore, Argo observations matter to constrain the model solution, even for an eddy-permitting model configuration. The impact of the Argo floats' data assimilation on other model variables is briefly assessed: the improvement of the fit to Argo profiles do not lead globally to unphysical corrections on the sea surface temperature and sea surface height. The main conclusion is that the performance of the Mercator Ocean 0.25° global data assimilation system is heavily dependent on the availability of Argo data.


Ocean Science ◽  
2012 ◽  
Vol 8 (5) ◽  
pp. 751-771 ◽  
Author(s):  
D. A. Ford ◽  
K. P. Edwards ◽  
D. Lea ◽  
R. M. Barciela ◽  
M. J. Martin ◽  
...  

Abstract. As part of the GlobColour project, daily chlorophyll a observations, derived using remotely sensed ocean colour data from the MERIS, MODIS and SeaWiFS sensors, are produced. The ability of these products to be assimilated into a pre-operational global coupled physical-biogeochemical model has been tested, on both a hindcast and near-real-time basis, and the impact on the system assessed. The assimilation was found to immediately and considerably improve the bias, root mean square error and correlation of modelled surface chlorophyll concentration compared to the GlobColour observations, an improvement which was sustained throughout the year and in every ocean basin. Errors against independent in situ chlorophyll observations were also reduced, both at and beneath the ocean surface. However, the model fit to in situ observations was not consistently better than that of climatology, due to errors in the underlying model. The assimilation scheme used is multivariate, updating all biogeochemical model state variables at all depths. The other variables were not degraded by the assimilation, with annual mean surface fields of nutrients, alkalinity and carbon variables remaining of similar quality compared to climatology. There was evidence of improved representation of zooplankton concentration, and reduced errors were seen against in situ observations of nitrate and pCO2, but too few observations were available to conclude about global model skill. The near-real-time GlobColour products were found to be sufficiently reliable for operational purposes, and of benefit to both operational-style systems and reanalyses.


Author(s):  
Nemesio Rodriguez-Fernandez ◽  
Patricia de Rosnay ◽  
Clement Albergel ◽  
Philippe Richaume ◽  
Filipe Aires ◽  
...  

The assimilation of Soil Moisture and Ocean Salinity (SMOS) data into the ECMWF (European Centre for Medium Range Weather Forecasts) H-TESSEL (Hydrology revised - Tiled ECMWF Scheme for Surface Exchanges over Land) model is presented. SMOS soil moisture (SM) estimates have been produced specifically by training a neural network with SMOS brightness temperatures as input and H-TESSEL model SM simulations as reference. This can help the assimilation of SMOS information in several ways: (1) the neural network soil moisture (NNSM) data have a similar climatology to the model, (2) no global bias is present with respect to the model even if regional differences can exist. Experiments performing joint data assimilation (DA) of NNSM, 2 metre air temperature and relative humidity or NNSM-only DA are discussed. The resulting SM was evaluated against a large number of in situ measurements of SM obtaining similar results to those of the model with no assimilation, even if significant differences were found from site to site. In addition, atmospheric forecasts initialized with H-TESSEL runs (without DA) or with the analysed SM were compared to measure of the impact of the satellite information. Although, NNSM DA has an overall neutral impact in the forecast in the Tropics, a significant positive impact was found in other areas and periods, especially in regions with limited in situ information. The joint NNSM, T2m and RH2m DA improves the forecast for all the seasons in the Southern Hemisphere. The impact is mostly due to T2m and RH2m, but SMOS NN DA alone also improves the forecast in July- September. In the Northern Hemisphere, the joint NNSM, T2m and RH2m DA improves the forecast in April-September, while NNSM alone has a significant positive effect in July-September. Furthermore, forecasting skill maps show that SMOS NNSM improves the forecast in North America and in Northern Asia for up to 72 hours lead time.


2015 ◽  
Vol 19 (12) ◽  
pp. 4831-4844 ◽  
Author(s):  
C. Draper ◽  
R. Reichle

Abstract. A 9 year record of Advanced Microwave Scanning Radiometer – Earth Observing System (AMSR-E) soil moisture retrievals are assimilated into the Catchment land surface model at four locations in the US. The assimilation is evaluated using the unbiased mean square error (ubMSE) relative to watershed-scale in situ observations, with the ubMSE separated into contributions from the subseasonal (SMshort), mean seasonal (SMseas), and inter-annual (SMlong) soil moisture dynamics. For near-surface soil moisture, the average ubMSE for Catchment without assimilation was (1.8 × 10−3 m3 m−3)2, of which 19 % was in SMlong, 26 % in SMseas, and 55 % in SMshort. The AMSR-E assimilation significantly reduced the total ubMSE at every site, with an average reduction of 33 %. Of this ubMSE reduction, 37 % occurred in SMlong, 24 % in SMseas, and 38 % in SMshort. For root-zone soil moisture, in situ observations were available at one site only, and the near-surface and root-zone results were very similar at this site. These results suggest that, in addition to the well-reported improvements in SMshort, assimilating a sufficiently long soil moisture data record can also improve the model representation of important long-term events, such as droughts. The improved agreement between the modeled and in situ SMseas is harder to interpret, given that mean seasonal cycle errors are systematic, and systematic errors are not typically targeted by (bias-blind) data assimilation. Finally, the use of 1-year subsets of the AMSR-E and Catchment soil moisture for estimating the observation-bias correction (rescaling) parameters is investigated. It is concluded that when only 1 year of data are available, the associated uncertainty in the rescaling parameters should not greatly reduce the average benefit gained from data assimilation, although locally and in extreme years there is a risk of increased errors.


2018 ◽  
Author(s):  
Benoît Tranchant ◽  
Elisabeth Remy ◽  
Eric Greiner ◽  
Olivier Legalloudec

Abstract. Monitoring Sea Surface Salinity (SSS) is important for understanding and forecasting the ocean circulation. It is even crucial in the context of the acceleration of the water cycle. Until recently, SSS was one of the less observed essential ocean variables. Only sparse in situ observations, most often closer to 5 meters deep than the surface, were available to estimate the SSS. The recent satellite missions of ESA's SMOS, NASA's Aquarius, and now SMAP have made possible for the first time to measure SSS from space. The SSS drivers can be quite different than the temperature ones. The model SSS can suffer from significant errors coming not only from the ocean dynamical model but also the atmospheric precipitation and evaporation as well as ice melting and river runoff. Satellite SSS can bring a valuable additional constraint to control the model salinity. In the framework of the SMOS Nino 2015 ESA project (https://www.godae-oceanview.org/projects/smos-nino15/), the impact of satellite SSS data assimilation is assessed with the Met Office and Mercator Ocean global ocean analysis and forecasting systems with a focus on the Tropical Pacific region. This article presents the analysis of an Observing System Experiment (OSE) conducted with the 1/4° resolution Mercator Ocean analysis and forecasting system. SSS data assimilation constrains the model SSS to be closer to the observations in a coherent way with the other data sets already routinely assimilated in an operational context. Globally, the SMOS SSS assimilation has a positive impact in salinity over the top 30 meters. Comparisons to independent data sets show a small but positive impact. The sea surface height (SSH) has also been impacted by implying a reinforcement of TIWs during the El-Niño 2015/16 event. Finally, this study helped us to progress in the understanding of the biases and errors that can degrade the SMOS SSS performance.


2009 ◽  
Vol 6 (11) ◽  
pp. 2333-2353 ◽  
Author(s):  
M. Vichi ◽  
S. Masina

Abstract. Global Ocean Biogeochemistry General Circulation Models are useful tools to study biogeochemical processes at global and large scales under current climate and future scenario conditions. The credibility of future estimates is however dependent on the model skill in capturing the observed multi-annual variability of firstly the mean bulk biogeochemical properties, and secondly the rates at which organic matter is processed within the food web. For this double purpose, the results of a multi-annual simulation of the global ocean biogeochemical model PELAGOS have been objectively compared with multi-variate observations from the last 20 years of the 20th century, both considering bulk variables and carbon production/consumption rates. Simulated net primary production (NPP) is comparable with satellite-derived estimates at the global scale and when compared with an independent data-set of in situ observations in the equatorial Pacific. The usage of objective skill indicators allowed us to demonstrate the importance of comparing like with like when considering carbon transformation processes. NPP scores improve substantially when in situ data are compared with modeled NPP which takes into account the excretion of freshly-produced dissolved organic carbon (DOC). It is thus recommended that DOC measurements be performed during in situ NPP measurements to quantify the actual production of organic carbon in the surface ocean. The chlorophyll bias in the Southern Ocean that affects this model as well as several others is linked to the inadequate representation of the mixed layer seasonal cycle in the region. A sensitivity experiment confirms that the artificial increase of mixed layer depths towards the observed values substantially reduces the bias. Our assessment results qualify the model for studies of carbon transformation in the surface ocean and metabolic balances. Within the limits of the model assumption and known biases, PELAGOS indicates a net heterotrophic balance especially in the more oligotrophic regions of the Atlantic during the boreal winter period. However, at the annual time scale and over the global ocean, the model suggests that the surface ocean is close to a weakly positive autotrophic balance in accordance with recent experimental findings and geochemical considerations.


Sign in / Sign up

Export Citation Format

Share Document