ALBERT

Description: 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.

Description: We consider the estimation of the grazing preferences parameters of zooplankton in ocean ecosystem models with ensemble-based Kalman filters. These parameters are introduced to model the relative diet composition of zooplankton that consists of phytoplankton, small size-classes of zooplankton and detritus. They are positive values and their sum is equal to one. However, the sum-to-one constraint cannot be guaranteed by ensemble-based Kalman filters when parameters are bounded. Therefore, a reformulation of the parameterization is proposed. We investigate two types of variable transformations for the estimation of positive sum-to-one constrained parameters that lead to the estimation of a new set of parameters with normal or bounded distributions. These transformations are illustrated and discussed with twin experiments performed with the 1-D coupled model GOTM-NORWECOM with Gaussian anamorphosis extensions of the deterministic ensemble Kalman filter (DEnKF).

Description: We review the procedures and challenges that must be considered when using geoid data derived from the Gravity and steady-state Ocean Circulation Explorer (GOCE) mission in order to constrain the circulation and water mass representation in an ocean general circulation model. It covers the combination of the geoid information with time-mean sea level information derived from satellite altimeter data, to construct a mean dynamic topography (MDT), and considers how this complements the time-varying sea level anomaly, also available from the satellite altimeter. We particularly consider the compatibility of these different fields in their spatial scale content, their temporal representation, and in their error covariances. These considerations are very important when the resulting data are to be used to estimate ocean circulation and its corresponding errors. We describe the further steps needed for assimilating the resulting dynamic topography information into an ocean circulation model using three different operational forecasting and data assimilation systems. We look at methods used for assimilating altimeter anomaly data in the absence of a suitable geoid, and then discuss different approaches which have been tried for assimilating the additional geoid information. We review the problems that have been encountered and the lessons learned in order the help future users. Finally we present some results from the use of GRACE geoid information in the operational oceanography community and discuss the future potential gains that may be obtained from a new GOCE geoid.

Description: Data assimilation methods often use an ensemble to represent the background error covariance. Two approaches are commonly used; a simple one with a static ensemble, or a more advanced one with a dynamic ensemble. The latter is often non-practical due to its high computational requirements. Some recent studies suggested using a hybrid covariance, which is a linear combination of the covariances represented by a static and a dynamic ensemble. Here, the use of the hybrid covariance is first extensively tested with a quasi-geostrophic model and with different analysis schemes, namely the Ensemble Kalman Filter (EnKF) and the Ensemble Square Root Filter (ESRF). The hybrid covariance ESRF (ESRF-OI) is more accurate and more stable than the hybrid covariance EnKF (EnKF-OI), but the overall conclusions are similar regardless of the analysis scheme used. The benefits of using the hybrid covariance are large compared to both the static and the dynamic methods with a small dynamic ensemble. The benefits over the dynamic methods become negligible, but remain, for large dynamic ensembles. The optimal value of the hybrid blending coefficient appears to decrease exponentially with the size of the dynamic ensemble. Finally, we consider a realistic application with the assimilation of altimetry data in a hybrid coordinate ocean model (HYCOM) for the Gulf of Mexico, during the shedding of Eddy Yankee (2006). A 10-member EnKF-OI is compared to a 10-member EnKF and a static method called the Ensemble Optimal Interpolation (EnOI). While 10 members seem insufficient for running the EnKF, the 10-member EnKF-OI reduces the forecast error compared to the EnOI, and improves the positions of the fronts.

Description: We consider the application of the Ensemble Kalman Filter (EnKF) to a coupled ocean ecosystem model (HYCOM-NORWECOM). Such models, especially the ecosystem models, are characterized by strongly non-linear interactions active in ocean blooms and present important limitations for the use of data assimilation methods based on linear statistical analysis. Besides the non-linearity of the model, one is confronted with physical/biological limitations, the analysis state having to be consistent with the model, especially with the constraints of positiveness of some variables. Furthermore the non-Gaussian distributions of the biogeochemical variables break an important assumption of the linear analysis, leading to a loss of optimality of the filter. We present an extension of the EnKF dealing with these limitations by introducing a non-linear change of variables (anamorphosis function) in order to execute the analysis step in a Gaussian space. We present also the initial results of the application of this non-Gaussian extension of the EnKF to the assimilation of simulated chlorophyll surface concentration data in a North Atlantic configuration of the HYCOM NORWECOM coupled model.

Description: A reanalysis of the North Atlantic spring bloom in 2007 was produced using the real-time analysis from the TOPAZ North Atlantic and Arctic forecasting system. The TOPAZ system uses a hybrid coordinate general circulation ocean model and assimilates physical observations: sea surface anomalies, sea surface temperatures, and sea-ice concentrations using the Ensemble Kalman Filter. This ocean model was coupled to an ecosystem model, NORWECOM (Norwegian Ecological Model System), and the TOPAZ-NORWECOM coupled model was run throughout the spring and summer of 2007. The ecosystem model was run online, restarting from analyzed physical fields (result after data assimilation) every 7 days. Biological variables were not assimilated in the model. The main purpose of the study was to investigate the impact of physical data assimilation on the ecosystem model. This was determined by comparing the results to those from a model without assimilation of physical data. The regions of focus are the North Atlantic and the Arctic Ocean. Assimilation of physical variables does not affect the results from the ecosystem model significantly. The differences between the weekly mean values of chlorophyll are normally within 5–10% during the summer months, and the maximum difference of ~20% occurs in the Arctic, also during summer. Special attention was paid to the nutrient input from the North Atlantic to the Nordic Seas and the impact of ice-assimilation on the ecosystem. The ice-assimilation increased the phytoplankton concentration: because there was less ice in the assimilation run, this increased both the mixing of nutrients during winter and the area where production could occur during summer. The forecast was also compared to remotely sensed chlorophyll, climatological nutrients, and in-situ data. The results show that the model reproduces a realistic annual cycle, but the chlorophyll concentrations tend to be between 0.1 and 1.0 mg chla/m3 too low during winter and spring and 1–2 mg chla/m3 too high during summer. Surface nutrients on the other hand are generally lower than the climatology throughout the year.

Description: Data assimilation methods often use an ensemble to represent the background error covariance. Two approaches are commonly used; a simple one with a static ensemble, or a more advanced one with a dynamic ensemble. The latter is often non-practical due to its high computational requirements. Some recent studies suggested using a hybrid covariance, which is a linear combination of the covariances represented by a static and a dynamic ensemble. Here, the use of the hybrid covariance is first extensively tested with a quasi-geostrophic model and with different analysis schemes, namely the Ensemble Kalman Filter (EnKF) and the Ensemble Square Root Filter (ESRF). The hybrid covariance ESRF (ESRF-OI) is more accurate and more stable than the hybrid covariance EnKF (EnKF-OI), but the overall conclusions are similar regardless of the analysis scheme used. The benefits of using the hybrid covariance are large compared to both the static and the dynamic methods with a small dynamic ensemble. The benefits over the dynamic methods become negligible, but remain, for large dynamic ensembles. The optimal value of the hybrid blending coefficient appears to decrease exponentially with the size of the dynamic ensemble. Finally, we consider a realistic application with the assimilation of altimetry data in a hybrid coordinate ocean model (HYCOM) for the Gulf of Mexico, during the shedding of Eddy Yankee (2006). A 10-member EnKF-OI is compared to a 10-member EnKF and a static method called the Ensemble Optimal Interpolation (EnOI). While 10 members seem insufficient for running the EnKF, the 10-member EnKF-OI reduces the forecast error compared to the EnOI, and improves the positions of the fronts.

Description: A 4th order advection scheme is applied in a nested eddy-resolving Hybrid Coordinate Ocean Model (HYCOM) of the greater Agulhas Current system for the purpose of testing advanced numerics as a means for improving the model simulation for eventual operational implementation. Model validation techniques comparing sea surface height variations, sea level skewness and variogram analyses to satellite altimetry measurements quantify that generally the 4th order advection scheme improves the realism of the model simulation. The most striking improvement over the standard 2nd order momentum advection scheme, is that the southern Agulhas Current is simulated as a well-defined meandering current, rather than a train of successive eddies. A better vertical structure and stronger poleward transports in the Agulhas Current core contribute toward a better southwestward penetration of the current, and its temperature field, implying a stronger Indo-Atlantic inter-ocean exchange. It is found that the transport, and hence this exchange, is sensitive to the occurrences of mesoscale features originating upstream in the Mozambique Channel and southern East Madagascar Current, and that the improved HYCOM simulation is well suited for further studies of these inter-actions.

Description: 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 are important because they interact with the main circulation. However, assimilation of altimetry data into a model including tides is challenging because tides and mesoscales 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). The 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 cm to 6.9 cm and improves the representation of the mesoscale features. With respect to the vertical temperature profiles, the data assimilation scheme improves the results at intermediate depth, but a slight degradation of the results at the surface is noted. The comparison to surface drifters shows an improvement of surface current by approximately −8.3%, with largest improvements in the Northern SCS and east of Vietnam.

Description: A reanalysis of the North Atlantic spring bloom in 2007 was produced using the real-time analyses from the TOPAZ (Towards an Operational Prediction system for the North Atlantic European coastal Zones) North Atlantic and Arctic forecasting system. The TOPAZ system uses a hybrid coordinate general circulation ocean model and assimilates physical observations: sea surface anomalies, sea surface temperatures, and sea-ice concentrations using the Ensemble Kalman Filter. This ocean model was coupled to an ecosystem model, NORWECOM (Norwegian Ecological Model System), and the TOPAZ-NORWECOM coupled model was run throughout the spring and summer of 2007. The ecosystem model was run online, restarting from analyzed physical fields (result after data assimilation) every 7 days. Biological variables were not assimilated in the model. The forecast was compared to remotely sensed chlorophyll and in-situ data. The impact of physical data assimilation on the ecosystem model was determined by comparing the results to those from a model without assimilation of physical data. The regions of focus are the North Atlantic and the Arctic Ocean. The results show that the model reproduces a realistic annual cycle, but the chlorophyll concentrations tend to be too low during winter and spring and too high during summer. Surface nutrients on the other hand are generally too low throughout the year. Assimilation of physical variables does not affect the results from the ecosystem model significantly. The differences between the weekly mean values of chlorophyll are normally within 5–10% during the summer months, and the maximum difference of ~20% occurs in the Arctic, also during summer. Special attention was paid to the nutrient input from the North Atlantic to the Nordic Seas and the impact of ice-assimilation on the ecosystem. The ice-assimilation increased the phytoplankton concentration: because there was less ice in the assimilation run, this increased both the mixing of nutrients during winter and the area where production could occur during summer.