Impact of assimilating a merged sea-ice thickness from CryoSat-2 and SMOS in the Arctic reanalysis

Accurately forecasting the sea-ice thickness (SIT) in the Arctic is a major challenge. The new SIT product (referred to as CS2SMOS) merges measurements from the CryoSat-2 and SMOS satellites on a weekly basis during the winter. The impact of assimilating CS2SMOS data is tested for the TOPAZ4 system – the Arctic component of the Copernicus Marine Environment Monitoring Services (CMEMS). TOPAZ4 currently assimilates a large set of ocean and sea-ice observations with the Deterministic Ensemble Kalman Filter (DEnKF). Two parallel reanalyses are conducted without (Official run) and with (Test run) assimilation of CS2SMOS data from 19 March 2014 to 31 March 2015. Since only mapping errors were provided in the CS2SMOS observation, an arbitrary term was added to compensate for the missing errors, but was found a posteriori too large. The SIT bias (too thin) is reduced from 16 to 5&thinsp;cm and the standard errors decrease from 53 to 38&thinsp;cm (by 28&thinsp;%) when compared to the assimilated SIT. When compared to independent SIT observations, the error reduction is 24&thinsp;% against the ice mass balance (IMB) buoy 2013F and by 12.5&thinsp;% against SIT data from the IceBridge campaigns. The improvement of sea-ice volume persists through the summer months in the absence of CS2SMOS data. Comparisons to sea-ice drift from the satellites show that dynamical adjustments reduce the drift errors around the North Pole by about 8&thinsp;%–9&thinsp;% in December 2014 and February 2015. Finally, using the degrees of freedom for signal (DFS), we find that CS2SMOS makes the prime source of information in the central Arctic and in the Kara Sea. We therefore recommend the assimilation of C2SMOS for Arctic reanalyses in order to improve the ice thickness and the ice drift.</p

An eddy resolving tidal-driven model of the South China Sea assimilating along-track SLA data using the EnOI

The upper ocean circulation in the South China Sea (SCS) is driven by the Asian monsoon, the Kuroshio intrusion through the Luzon Strait, strong tidal currents, and a complex topography. Here, we demonstrate the benefit of assimilating along-track altimeter data into a nested configuration of the HYbrid Coordinate Ocean Model that includes tides. Including tides in models is important because they interact with the main circulation. However, assimilation of altimetry data into a model including tides is challenging because tides and mesoscale features contribute to the elevation of ocean surface at different time scales and require different corrections. To address this issue, tides are filtered out of the model output and only the mesoscale variability is corrected with a computationally cheap data assimilation method: the Ensemble Optimal Interpolation (EnOI). This method uses a running selection of members to handle the seasonal variability and assimilates the track data asynchronously. The data assimilative system is tested for the period 1994–1995, during which time a large number of validation data are available. Data assimilation reduces the Root Mean Square Error of Sea Level Anomalies from 9.3 to 6.9 cm and improves the representation of the mesoscale features. With respect to the vertical temperature profiles, the data assimilation scheme reduces the errors quantitatively with an improvement at intermediate depth and deterioration at deeper depth. The comparison to surface drifters shows an improvement of surface current by approximately −9% in the Northern SCS and east of Vietnam. Results are improved compared to an assimilative system that does not include tides and a system that does not consider asynchronous assimilation

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

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

Scaling Observation Error for Optimal Assimilation of CCI SST Data into a Regional HYCOM EnOI System

South Africa currently possesses no operational ocean forecasting system for the purpose of predicting ocean state variables including temperature,salinity and velocity. Substantial initial efforts towards this goal have been made and resulted in a system using a regional Hybrid Coordinate Ocean Model (HYCOM) along with the Ensemble Optimal Interpolation (EnOI)assimilation scheme. Assimilating only sea surface temperature (SST) observations from the Operational Sea Surface Temperature and Sea Ice Analysis (OSTIA) product into the system resulted in a degraded forecast. Aiming to address this, Climate Change Initiative (CCI) SSTs are assimilated into the system in an effort to improve the forecast skill. Observation errors in the assimilated product are used in the EnOI to determine whether more confidence should be placed in the model or observations in producing the analysis, but overconfidence in observations can shock the model and result in failure. To tweak the impact of the assimilation, a scaling factor is applied in the assimilation code. A scaling factor of 25 was found to produce a favourable result with lowest mean root mean square error (RMSE;1.098C) between the model and observations over time. Postulating the error to be overconfident, a floor value is introduced in order to set a minimum value for the observation error thereby reducing confidence in the observations. These experiments fared less favourably with a floor value of 0.5 and a scaling factor of 15 producing the best mean RMSE (1.118C)

Optimising assimilation of sea ice concentration in an Earth system model with a multicategory sea ice model

A data assimilation method capable of constraining the sea ice of an Earth system model in a dynamically consistent manner has the potential to enhance the accuracy of climate reconstructions and predictions. Finding such a method is challenging because the sea ice dynamics is highly non-linear, and sea ice variables are strongly non-Gaussian distributed and tightly coupled to the rest of the Earth system - particularly thermodynamically with the ocean. We investigate key practical implementations for assimilating sea ice concentration - the predominant source of observations in polar regions - with the Norwegian Climate Prediction Model that combines the Norwegian Earth System Model with the Ensemble Kalman Filter. The performances of the different configurations are investigated by conducting 10-year reanalyses in a perfect model framework. First, we find that with a flow-dependent assimilation method, strongly coupled ocean-sea ice assimilation outperforms weakly coupled (sea ice only) assimilation. An attempt to prescribe the covariance between the ocean temperature and the sea ice concentration performed poorly. Extending the ocean updates below the mixed layer is slightly beneficial for the Arctic hydrography. Second, we find that solving the analysis for the multicategory instead of the aggregated ice state variables greatly reduces the errors in the ice state. Updating the ice volumes induces a weak drift in the bias for the thick ice category that relates to the postprocessing of unphysical thicknesses. Preserving the ice thicknesses for each category during the assimilation mitigates the drift without degrading the performance. The robustness and reliability of the optimal setting is demonstrated for a 20-year reanalysis. The error of sea ice concentration reduces by 50% (65%), sea ice thickness by 25% (35%), sea surface temperature by 33% (23%) and sea surface salinity by 11% (25%) in the Arctic (Antarctic) compared to a reference run without assimilation

Supermodeling Improving Predictions with an Ensemble of Interacting Models

The modeling of weather and climate has been a success story. The skill of forecasts continues to improve and model biases continue to decrease. Combining the output of multiple models has further improved forecast skill and reduced biases. But are we exploiting the full capacity of state-of-the-art models in making forecasts and projections? Supermodeling is a recent step forward in the multimodel ensemble approach. Instead of combining model output after the simulations are completed, in a supermodel individual models exchange state information as they run, influencing each other's behavior. By learning the optimal parameters that determine how models influence each other based on past observations, model errors are reduced at an early stage before they propagate into larger scales and affect other regions and variables. The models synchronize on a common solution that through learning remains closer to the observed evolution. Effectively a new dynamical system has been created, a supermodel, that optimally combines the strengths of the constituent models. The supermodel approach has the potential to rapidly improve current state-of-the-art weather forecasts and climate predictions. In this paper we introduce supermodeling, demonstrate its potential in examples of various complexity, and discuss learning strategies. We conclude with a discussion of remaining challenges for a successful application of supermodeling in the context of state-of-the-art models. The supermodeling approach is not limited to the modeling of weather and climate, but can be applied to improve the prediction capabilities of any complex system, for which a set of different models exists

Synoptic-scale analysis of mechanisms driving surface chlorophyll dynamics in the North Atlantic

Several hypotheses have been proposed for the onset of the spring phytoplankton bloom in the North Atlantic. Our main objective is to examine which bottom-up processes can best predict the annual increase in surface phytoplankton concentration in the North Atlantic by applying novel phenology algorithms to ocean colour data. We construct indicator fields and time series which, in various combinations, provide models consistent with the principle dynamics previously proposed. Using a multimodel inference approach, we investigate the evidence supporting these models and how it varies in space. We show that, in terms of bottom-up processes alone, there is a dominant physical mechanism, namely mixed-layer shoaling, that best predicts the interannual variation in the initial increase in surface chlorophyll across large sectors of the North Atlantic. We further show that different regions are governed by different physical phenomena and that wind-driven mixing is a common component, with either heat flux or light as triggers. We believe these findings to be relevant to the ongoing discussion on North Atlantic bloom onset

Propagation of Thermohaline Anomalies and Their Predictive Potential along the Atlantic Water Pathway

We assess to what extent seven state-of-the-art dynamical prediction systems can retrospectively predict winter sea surface temperature (SST) in the subpolar North Atlantic and the Nordic seas in the period 1970-2005. We focus on the region where warm water flows poleward (i.e., the Atlantic water pathway to the Arctic) and on interannual-to-decadal time scales. Observational studies demonstrate predictability several years in advance in this region, but we find that SST skill is low with significant skill only at a lead time of 1-2 years. To better understand why the prediction systems have predictive skill or lack thereof, we assess the skill of the systems to reproduce a spatiotemporal SST pattern based on observations. The physical mechanism underlying this pattern is a propagation of oceanic anomalies from low to high latitudes along the major currents, the North Atlantic Current and the Norwegian Atlantic Current. We find that the prediction systems have difficulties in reproducing this pattern. To identify whether the misrepresentation is due to incorrect model physics, we assess the respective uninitialized historical simulations. These simulations also tend to misrepresent the spatiotemporal SST pattern, indicating that the physical mechanism is not properly simulated. However, the representation of the pattern is slightly degraded in the predictions compared to historical runs, which could be a result of initialization shocks and forecast drift effects. Ways to enhance predictions could include improved initialization and better simulation of poleward circulation of anomalies. This might require model resolutions in which flow over complex bathymetry and the physics of mesoscale ocean eddies and their interactions with the atmosphere are resolved. SIGNIFICANCE STATEMENT: In this study, we find that dynamical prediction systems and their respective climate models struggle to realistically represent ocean surface temperature variability in the eastern subpolar North Atlantic and Nordic seas on interannual-to-decadal time scales. In previous studies, ocean advection is proposed as a key mechanism in propagating temperature anomalies along the Atlantic water pathway toward the Arctic Ocean. Our analysis suggests that the predicted temperature anomalies are not properly circulated to the north; this is a result of model errors that seems to be exacerbated by the effect of initialization shocks and forecast drift. Better climate predictions in the study region will thus require improving the initialization step, as well as enhancing process representation in the climate models

Seasonal-to-decadal predictions with the ensemble Kalman filter and the Norwegian Earth System Model: a twin experiment

Here, we firstly demonstrate the potential of an advanced flow dependent data assimilation method for performing seasonal-to-decadal prediction and secondly, reassess the use of sea surface temperature (SST) for initialisation of these forecasts. We use the Norwegian Climate Prediction Model (NorCPM), which is based on the Norwegian Earth System Model (NorESM) and uses the deterministic ensemble Kalman filter to assimilate observations. NorESM is a fully coupled system based on the Community Earth System Model version 1, which includes an ocean, an atmosphere, a sea ice and a land model. A numerically efficient coarse resolution version of NorESM is used. We employ a twin experiment methodology to provide an upper estimate of predictability in our model framework (i.e. without considering model bias) of NorCPM that assimilates synthetic monthly SST data (EnKF-SST). The accuracy of EnKF-SST is compared to an unconstrained ensemble run (FREE) and ensemble predictions made with near perfect (i.e. microscopic SST perturbation) initial conditions (PERFECT). We perform 10 cycles, each consisting of a 10-yr assimilation phase, followed by a 10-yr prediction. The results indicate that EnKF-SST improves sea level, ice concentration, 2 m atmospheric temperature, precipitation and 3-D hydrography compared to FREE. Improvements for the hydrography are largest near the surface and are retained for longer periods at depth. Benefits in salinity are retained for longer periods compared to temperature. Near-surface improvements are largest in the tropics, while improvements at intermediate depths are found in regions of large-scale currents, regions of deep convection, and at the Mediterranean Sea outflow. However, the benefits are often small compared to PERFECT, in particular, at depth suggesting that more observations should be assimilated in addition to SST. The EnKF-SST system is also tested for standard ocean circulation indices and demonstrates decadal predictability for Atlantic overturning and sub-polar gyre circulations, and heat content in the Nordic Seas. The system beats persistence forecast and shows skill for heat content in the Nordic Seas that is close to PERFECT

Advances in reconstructing the AMOC using sea surface observations of salinity

The Atlantic meridional overturning circulation (AMOC) is one of the main drivers of climate variability at decadal and longer time scales. As there are no direct multi-decadal observations of this key circulation, the reconstruction of past AMOC variations is essential. This work presents a step forward in reconstructing the AMOC using climate models and time-varying surface nudging of salinity and temperature data, for which independent multi-decadal observed series are available. A number of nudging protocols are explored in a perfect model framework to best reproduce the AMOC variability accommodating to the characteristics of SST and SSS available products. As reference SST products with sufficient space and time coverage are available, we here choose to focus on the limitations associated to SSS products with the goal of providing protocols using independent salinity products. We consider a global gridded dataset and, additionally, a coarser SSS dataset restricted to the Atlantic and with a quite low spatial resolution (order of 10 degrees vs. 2 for the model grid). We show how, using the latter, we can improve the efficiency of the nudging on the AMOC reconstruction by adding a high-resolution annual cycle to the coarse resolution SSS product as well as a spatial downscaling to account for SSS gradient. The final protocol retained for the coarse SSS data is able to reconstruct a 100-year long AMOC period (average of 10.18 Sv and a standard deviation of 1.39 Sv), with a correlation of 0.76 to the target and a RMSE of 0.99 Sv. These values can be respectively compared to 0.85 and 0.75 Sv when using the global salinity surface observations. This work provides a first step towards understanding the limitations and prospects of historical AMOC reconstructions using different sea surface salinity datasets for the surface nudging