ALBERT

Description: The European Reanalysis of Global Climate Observations 2 (ERA-CLIM2) is a European Union Seventh Framework Project started in January 2014 and due to be completed in December 2017. It aims to produce coupled reanalyses, which are physically consistent datasets describing the evolution of the global atmosphere, ocean, land surface, cryosphere, and the carbon cycle. ERA-CLIM2 has contributed to advancing the capacity for producing state-of-the-art climate reanalyses that extend back to the early twentieth century. ERA-CLIM2 has led to the generation of the first European ensemble of coupled ocean, sea ice, land, and atmosphere reanalyses of the twentieth century. The project has funded work to rescue and prepare observations and to advance the data-assimilation systems required to generate operational reanalyses, such as the ones planned by the European Union Copernicus Climate Change Service. This paper summarizes the main goals of the project, discusses some of its main areas of activities, and presents some of its key results.

Description: Ocean data assimilation systems encompass a wide range of scales that are difficult to control simultaneously using partial observation networks. All scales are not observable by all observation systems which is not easily taken into account in current ocean operational systems. The main reason for this difficulty is that the error covariance matrices are usually assumed to be local (e.g. using a localization algorithm in ensemble data assimilation systems), so that the large scale patterns are removed from the error statistics. To better exploit the observational information available for all scales in CMEMS assimilation systems, we investigate a new method to introduce scale separation in the assimilation scheme. The method is based on a spectral transformation of the assimilation problem and consists in carrying out the analysis with spectral localisation for the large scales and spatial localisation for the residual scales. The target is to improve the observational update of the large scale components of the signal by an explicit observational constraint applied directly on the large scales, and to restrict the use of spatial localisation to the small scale components of the signal. To evaluate our method, twin experiments are carried out with synthetic altimetry observations (simulating the JASON tracks), assimilated in a 1/4° model configuration of the North Atlantic and the Nordic Seas. Results show that the transformation to the spectral space and the spectral localization provides consistent ensemble estimates of the state of the system (in the spectral space,or after backward transformation to the spatial space). Combined with spatial localisation for the residual scales, the new scheme is able to provide a reliable ensemble update for all scales, with improved accuracy for the large scale; and the performance of the system can be checked explicitly and separately for all scales in the assimilation system.

Description: The authors investigate the variability of salinity in the Arctic Ocean and in the Nordic and Labrador Seas over recent years to see how the freshwater balance in the Arctic and the exchanges with the North Atlantic have been affected by the recent important sea ice melting, especially during the 2007 sea ice extent minimum. The Global Ocean Reanalysis and Simulations (GLORYS1) global ocean reanalysis based on a global coupled ocean–sea ice model with an average of 12-km grid resolution in the Arctic Ocean is used in this regard. Although no sea ice data and no data under sea ice are assimilated, simulation over the 2001–09 period is shown to represent fairly well the 2007 sea ice event and the different components accounting for the ocean and sea ice freshwater budget, compared to available observations. In the reanalysis, the 2007 sea ice minimum is due to an increase of the sea ice export through Fram Strait (25%) and an important sea ice melt in the Arctic (75%). Liquid freshwater is accumulated in the Beaufort gyre after 2002, in agreement with recent observations, and it is shown that this accumulation is due to both the sea ice melt and a spatial redistribution of the freshwater content in the Canadian Basin. In the Eurasian Basin, a very contrasting situation is found with an increase of the salinity. The effect of the sea ice melt is counterbalanced by an increase of the Atlantic inflow and a modification of the circulation north of Fram Strait after 2007. The authors suggest that a strong anomaly of the atmospheric conditions was responsible for this change of the circulation.

Description: In the Kalman filter standard algorithm, the computational complexity of the observational update is proportional to the cube of the number y of observations (leading behavior for large y). In realistic atmospheric or oceanic applications, involving an increasing quantity of available observations, this often leads to a prohibitive cost and to the necessity of simplifying the problem by aggregating or dropping observations. If the filter error covariance matrices are in square root form, as in square root or ensemble Kalman filters, the standard algorithm can be transformed to be linear in y, providing that the observation error covariance matrix is diagonal. This is a significant drawback of this transformed algorithm and often leads to an assumption of uncorrelated observation errors for the sake of numerical efficiency. In this paper, it is shown that the linearity of the transformed algorithm in y can be preserved for other forms of the observation error covariance matrix. In particular, quite general correlation structures (with analytic asymptotic expressions) can be simulated simply by augmenting the observation vector with differences of the original observations, such as their discrete gradients. Errors in ocean altimetric observations are spatially correlated, as for instance orbit or atmospheric errors along the satellite track. Adequately parameterizing these correlations can directly improve the quality of observational updates and the accuracy of the associated error estimates. In this paper, the example of the North Brazil Current circulation is used to demonstrate the importance of this effect, which is especially significant in that region of moderate ratio between signal amplitude and observation noise, and to show that the efficient parameterization that is proposed for the observation error correlations is appropriate to take it into account. Adding explicit gradient observations also receives a physical justification. This parameterization is thus proved to be useful to ocean data assimilation systems that are based on square root or ensemble Kalman filters, as soon as the number of observations becomes penalizing, and if a sophisticated parameterization of the observation error correlations is required.

Description: In large-sized atmospheric or oceanic applications of square root or ensemble Kalman filters, it is often necessary to introduce the prior assumption that long-range correlations are negligible and force them to zero using a local parameterization, supplementing the ensemble or reduced-rank representation of the covariance. One classic algorithm to perform this operation consists of taking the Schur product of the ensemble covariance matrix by a local support correlation matrix. However, with this parameterization, the square root of the forecast error covariance matrix is no more directly available, so that any observational update algorithm requiring this square root must include an additional step to compute local square roots from the Schur product. This computation generates an additional numerical cost or produces high-rank square roots, which may deprive the observational update from its original efficiency. In this paper, it is shown how efficient local square root parameterizations can be obtained, for use with a specific square root formulation (i.e., eigenbasis algorithm) of the observational update. Comparisons with the classic algorithm are provided, mainly in terms of consistency, accuracy, and computational complexity. As an application, the resulting parameterization is used to estimate maps of dynamic topography characterizing a basin-scale ocean turbulent flow. Even with this moderate-sized system (a 2200-km-wide square basin with 100-km-wide mesoscale eddies), it is observed that more than 1000 ensemble members are necessary to faithfully represent the global correlation patterns, and that a local parameterization is needed to produce correct covariances with moderate-sized ensembles. Comparisons with the exact solution show that the use of local square roots is able to improve the accuracy of the updated ensemble mean and the consistency of the updated ensemble variance. With the eigenbasis algorithm, optimal adaptive estimates of scaling factors for the forecast and observation error covariance matrix can also be obtained locally at negligible additional numerical cost. Finally, a comparison of the overall computational cost illustrates the decisive advantage that efficient local square root parameterizations may have to deal simultaneously with a larger number of observations and avoid data thinning as much as possible.

Description: In Kalman filter applications, an adaptive parameterization of the error statistics is often necessary to avoid filter divergence, and prevent error estimates from becoming grossly inconsistent with the real error. With the classic formulation of the Kalman filter observational update, optimal estimates of general adaptive parameters can only be obtained at a numerical cost that is several times larger than the cost of the state observational update. In this paper, it is shown that there exists a few types of important parameters for which optimal estimates can be computed at a negligible numerical cost, as soon as the computation is performed using a transformed algorithm that works in the reduced control space defined by the square root or ensemble representation of the forecast error covariance matrix. The set of parameters that can be efficiently controlled includes scaling factors for the forecast error covariance matrix, scaling factors for the observation error covariance matrix, or even a scaling factor for the observation error correlation length scale. As an application, the resulting adaptive filter is used to estimate the time evolution of ocean mesoscale signals using observations of the ocean dynamic topography. To check the behavior of the adaptive mechanism, this is done in the context of idealized experiments, in which model error and observation error statistics are known. This ideal framework is particularly appropriate to explore the ill-conditioned situations (inadequate prior assumptions or uncontrollability of the parameters) in which adaptivity can be misleading. Overall, the experiments show that, if used correctly, the efficient optimal adaptive algorithm proposed in this paper introduces useful supplementary degrees of freedom in the estimation problem, and that the direct control of these statistical parameters by the observations increases the robustness of the error estimates and thus the optimality of the resulting Kalman filter.

Description: Since October 19, 2016, and in the framework of Copernicus Marine Environment Monitoring Service (CMEMS), Mercator Ocean delivers in 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 3D multivariate modal decomposition of the forecast 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 of 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 forecast error covariance and the adaptive tuning of observations 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 products quality is also discussed, highlighting a gain in performance and reliability of the current global monitoring and forecasting system compared to its previous version.

Description: We generate ocean biogeochemical model ensembles including several kinds of stochastic parameterizations. The NEMO stochastic modules are complemented by integrating a subroutine to calculate variable anisotropic spatial scales, which are of particular importance in high-resolution coastal configurations. The domain covers the Bay of Biscay at 1/36° resolution, as a case study for open-ocean and coastal shelf dynamics. At first, we identify uncertainties from assumptions subject to erroneous atmospheric forcing, ocean model improper parameterizations and ecosystem state uncertainties. The error regimes are found to be mainly driven by the wind forcing, with the rest of the perturbed tendencies locally augmenting the ensemble spread. Biogeochemical uncertainties arise from inborn ecosystem model errors and from errors in the physical state. Model errors in physics are found to have larger impact on chlorophyll spread than those of the ecosystem. In a second step, the ensembles undergo verification with respect to observations, focusing on upper-ocean properties. We investigate the statistical consistency of prior model errors and observation estimates, in view of joint uncertainty vicinities, associated with both sources of information. OSTIA-SST L4 distribution appears to be compatible with ensembles perturbing physics, since vicinities overlap, enabling data assimilation. The most consistent configuration for SLA along-track L3 data is in the Abyssal plain, where the spread is increased due to mesoscale eddy decorrelation. The largest statistical SLA biases are observed in coastal regions, sometimes to the point that vicinities become disjoint. Missing error processes in relation to SLA hint at the presence of high-frequency error sources currently unaccounted for, potentially leading to ill-posed assimilation problems. Ecosystem model-data samples with respect to Ocean Colour L4 appear to be compatible with each other only at times, with data assimilation being marginally well-posed. In a third step, we illustrate the potential influence of those uncertainties on data assimilation impact exercise, by means of multivariate representers and EnKF-type incremental analysis for a few members. Corrections on physical properties are associated with large-scale biases between model and data, with diverse characteristics in the open-ocean and the shelves. The increments are often characteristic of the underlying mesoscale features, chlorophyll included due to the vertical velocity field. Small scale local corrections are visible over the shelves. Chlorophyll information seems to have a very measurable potential impact on physical variables.

Description: Ocean data assimilation systems encompass a wide range of scales that are difficult to control simultaneously using partial observation networks. All scales are not observable by all observation systems, which is not easily taken into account in current ocean operational systems. The main reason for this difficulty is that the error covariance matrices are usually assumed to be local (e.g. using a localisation algorithm in ensemble data assimilation systems), so that the large-scale patterns are removed from the error statistics. To better exploit the observational information available for all scales in the assimilation systems of the Copernicus Marine Environment Monitoring Service, we investigate a new method to introduce scale separation in the assimilation scheme. The method is based on a spectral transformation of the assimilation problem and consists in carrying out the analysis with spectral localisation for the large scales and spatial localisation for the residual scales. The target is to improve the observational update of the large-scale components of the signal by an explicit observational constraint applied directly on the large scales and to restrict the use of spatial localisation to the small-scale components of the signal. To evaluate our method, twin experiments are carried out with synthetic altimetry observations (simulating the Jason tracks), assimilated in a 1/4∘ model configuration of the North Atlantic and the Nordic Seas. Results show that the transformation to the spectral domain and the spectral localisation provides consistent ensemble estimates of the state of the system (in the spectral domain or after backward transformation to the spatial domain). Combined with spatial localisation for the residual scales, the new scheme is able to provide a reliable ensemble update for all scales, with improved accuracy for the large scale; and the performance of the system can be checked explicitly and separately for all scales in the assimilation system.