4DVAR data assimilation in the Intra-Americas Sea with the Regional Ocean Modeling System (ROMS)

We present the background, development, and preparation of a state-of-the-art 4D variational (4DVAR) data assimilation system in the Regional Ocean Modeling System (ROMS) with an application in the Intra-Americas Sea (IAS). This initial application with a coarse model shows the efficacy of the 4DVAR methodology for use within complex ocean environments, and serves as preparation for deploying an operational, real-time assimilation system onboard the Royal Caribbean Cruise Lines ship Explorer of the Seas. Assimilating satellite sea surface height and temperature observations with in situ data from the ship in 14 day cycles over 2 years from January 2005 through March 2007, reduces the observation-model misfit by over 75%. Using measures of the Loop Current dynamics, we show that the assimilated solution is consistent with observed statistics.     

4DVAR Data Assimilation with the Regional Ocean Modeling System (ROMS): Impact on the Water Mass Distributions in the Yellow Sea

In this study, we evaluate the performance of the recently developed incremental strong constraint 4-dimensional variational (4DVAR) data assimilation applied to the Yellow Sea (YS) using the Regional Ocean Modeling System (ROMS). Two assimilation experiments are compared: assimilating remote-sensed sea surface temperature (SST) and both the SST and in-situ profiles measured by shipboard CTD casts into a regional ocean modeling from January to December of 2011. By comparing the two assimilation experiments against a free-run without data assimilation, we investigate how the assimilation affects the hydrographic structures in the YS. Results indicate that the SST assimilation notably improves the model behavior at the surface when compared to the non-assimilative free-run. The SST assimilation also has an impact on the subsurface water structure in the eastern YS; however, the improvement is seasonally dependent, that is, the correction becomes more effective in winter than in summer. This is due to a strong stratification in summer that prevents the assimilation of SST from affecting the subsurface temperature. A significant improvement to the subsurface temperature is made when the in-situ profiles of temperature and salinity are assimilated, forming a tongue-shaped YS bottom cold water from the YS toward the southwestern seas of Jeju Island.     

基于ROMS模式的南海SST与SSH四维变分同化研究

卫星遥感观测获得了大量高分辨率的海面实时信息,包括海面温度（SST）和海面高度（SSH）等,同化进入数值模式可有效提升模拟精度。本文基于ROMS模式与四维变分同化方法（4DVAR）,使用AVHRR SST和AVISO SSH数据,开展了南海区域同化实验。为检验同化的效果,分别利用HYCOM再分析资料和Argo温盐实测数据分析了同化结果的海面高度、流场及温盐剖面的精度。对比结果表明,SST和SSH的同化能够改善ROMS的模拟结果:同化后海面高度场能够更为准确地捕捉海洋的中尺度特征,与HYCOM海面高度再分析资料相比,平均绝对偏差和均方根误差分别为0.054 m和0.066 m;与HYCOM 10 m层流场相比,东向与北向流速平均绝对偏差分别为0.12 m/s和0.11 m/s,相比未同化均提升约0.01 m/s;温盐同化结果与Argo温盐实测具有较高的一致性,温度和盐度平均绝对偏差为0.45℃、0.077,均方根误差为0.91℃、0.11,单个的温盐廓线对比说明,同化结果与HYCOM再分析资料精度相当。     

基于ROMS和4DVAR的沿轨与网格化SSH数据同化效果评价

Remote sensing products are significant in the data assimilation of an ocean model. Considering the resolution and space coverage of different remote sensing data, two types of sea surface height(SSH) product are employed in the assimilation, including the gridded products from AVISO and the original along-track observations used in the generation. To explore their impact on the assimilation results, an experiment focus on the South China Sea(SCS) is conducted based on the Regional Ocean Modeling System(ROMS) and the four-dimensional variational data assimilation(4 DVAR) technology. The comparison with EN4 data set and Argo profile indicates that, the along-track SSH assimilation result presents to be more accurate than the gridded SSH assimilation, because some noises may have been introduced in the merging process. Moreover, the mesoscale eddy detection capability of the assimilation results is analyzed by a vector geometry–based algorithm. It is verified that, the assimilation of the gridded SSH shows superiority in describing the eddy's characteristics, since the complete structure of the ocean surface has been reconstructed by the original data merging.     

The Regional Ocean Modeling System (ROMS) 4-dimensional variational data assimilation systems: Part I - System overview and formulation

The Regional Ocean Modeling System (ROMS) is one of the few community ocean general circulation models for which a 4-dimensional variational data assimilation (4D-Var) capability has been developed. The ROMS 4D-Var capability is unique in that three variants of 4D-Var are supported: a primal formulation of incremental strong constraint 4D-Var (I4D-Var), a dual formulation based on a physical-space statistical analysis system (4D-PSAS), and a dual formulation representer-based variant of 4D-Var (R4D-Var). In each case, ROMS is used in conjunction with available observations to identify a best estimate of the ocean circulation based on a set of a priori hypotheses about errors in the initial conditions, boundary conditions, surface forcing, and errors in the model in the case of 4D-PSAS and R4D-Var. In the primal formulation of I4D-Var the search for the best circulation estimate is performed in the full space of the model control vector, while for the dual formulations of 4D-PSAS and R4D-Var only the sub-space of linear functions of the model state vector spanned by the observations (i.e. the dual space) is searched. In oceanographic applications, the number of observations is typically much less than the dimension of the model control vector, so there are clear advantages to limiting the search to the space spanned by the observations. In the case of 4D-PSAS and R4D-Var, the strong constraint assumption (i.e. that the model is error free) can be relaxed leading to the so-called weak constraint formulation. This paper describes the three aforementioned variants of 4D-Var as they are implemented in ROMS. Critical components that are common to each approach are conjugate gradient descent, preconditioning, and error covariance models, which are also described. Finally, several powerful 4D-Var diagnostic tools are discussed, namely computation of posterior errors, eigenvector analysis of the posterior error covariance, observation impact, and observation sensitivity.     

Weak and strong constraint data assimilation in the inverse Regional Ocean Modeling System (ROMS): Development and application for a baroclinic coastal upwelling system

We describe the development and preliminary application of the inverse Regional Ocean Modeling System (ROMS), a four dimensional variational (4DVAR) data assimilation system for high-resolution basin-wide and coastal oceanic flows. Inverse ROMS makes use of the recently developed perturbation tangent linear (TL), representer tangent linear (RP) and adjoint (AD) models to implement an indirect representer-based generalized inverse modeling system. This modeling framework is modular. The TL, RP and AD models are used as stand-alone sub-models within the Inverse Ocean Modeling (IOM) system described in [Chua, B.S., Bennett, A.F., 2001. An inverse ocean modeling system. Ocean Modell. 3, 137–165.]. The system allows the assimilation of a wide range of observation types and uses an iterative algorithm to solve nonlinear assimilation problems. The assimilation is performed either under the perfect model assumption (strong constraint) or by also allowing for errors in the model dynamics (weak constraints). For the weak constraint case the TL and RP models are modified to include additional forcing terms on the right hand side of the model equations. These terms are needed to account for errors in the model dynamics.Inverse ROMS is tested in a realistic 3D baroclinic upwelling system with complex bottom topography, characterized by strong mesoscale eddy variability. We assimilate synthetic data for upper ocean (0–450 m) temperatures and currents over a period of 10 days using both a high resolution and a spatially and temporally aliased sampling array. During the assimilation period the flow field undergoes substantial changes from the initial state. This allows the inverse solution to extract the dynamically active information from the synthetic observations and improve the trajectory of the model state beyond the assimilation window. Both the strong and weak constraint assimilation experiments show forecast skill greater than persistence and climatology during the 10–20 days after the last observation is assimilated.Further investigation in the functional form of the model error covariance and in the use of the representer tangent linear model may lead to improvement in the forecast skill.     

The Regional Ocean Modeling System (ROMS) 4-dimensional variational data assimilation systems: Part II - Performance and application to the California Current System

The Regional Ocean Modeling System (ROMS) 4-dimensional variational (4D-Var) data assimilation systems have been systematically applied to the mesoscale circulation environment of the California Current to demonstrate the performance and practical utility of the various components of ROMS 4D-Var. In particular, we present a comparison of three approaches to 4D-Var, namely: the primal formulation of the incremental strong constraint approach; the dual formulation “physical-space statistical analysis system”; and the dual formulation indirect representer approach. In agreement with theoretical considerations all three approaches converge to the same ocean circulation estimate when using the same observations and prior information. However, the rate of convergence of the dual formulation was found to be inferior to that of the primal formulation. Other aspects of the 4D-Var performance that relate to the use of multiple outer-loops, preconditioning, and the weak constraint are also explored. A systematic evaluation of the impact of the various components of the 4D-Var control vector (i.e. the initial conditions, surface forcing and open boundary conditions) is also presented. It is shown that correcting for uncertainties in the model initial conditions exerts the largest influence on the ability of the model to fit the available observations. Various important diagnostics of 4D-Var are also examined, including estimates of the posterior error, the information content of the observation array, and innovation-based consistency checks on the prior error assumptions. Using these diagnostic tools, we find that more than 90% of the observations assimilated into the model provide redundant information. This is a symptom of the large percentage of satellite data that are used and to some extent the nature of the data processing employed. This is the second in a series of three papers describing the ROMS 4D-Var systems.     

The Regional Ocean Modeling System (ROMS) 4-dimensional variational data assimilation systems: Part III - Observation impact and observation sensitivity in the California Current System

The critical role played by observations during ocean data assimilation was explored when the Regional Ocean Modeling System (ROMS) 4-dimensional variational (4D-Var) data assimilation system was applied sequentially to the California Current circulation. The adjoint of the 4D-Var gain matrix was used to quantify the impact of individual observations and observation platforms on different aspects of the 4D-Var circulation estimates during both analysis and subsequent forecast cycles. In this study we focus on the alongshore and cross-shore transport of the California Current System associated with wind-induced coastal upwelling along the central California coast. The majority of the observations available during any given analysis cycle are from satellite platforms in the form of SST and SSH, and on average these data exert the largest controlling influence on the analysis increments and forecast skill of coastal transport. However, subsurface in situ observations from Argo floats, CTDs, XBTs and tagged marine mammals often have a considerable impact on analyses and forecasts of coastal transport, even though these observations represent a relatively small fraction of the available data at any particular time.During 4D-Var the observations are used to correct for uncertainties in the model control variables, namely the initial conditions, surface forcing, and open boundary conditions. It is found that correcting for uncertainties in both the initial conditions and surface forcing has the largest impact on the analysis increments in alongshore transport, while the cross-shore transport is controlled mainly by the surface forcing. The memory of the circulation associated with the control variable increments was also explored in relation to 7 day forecasts of the coastal circulation. Despite the importance of correcting for surface forcing uncertainties during analysis cycles, the coastal transport during forecast cycles initialized from the analyses has less memory of the surface forcing corrections, and is controlled primarily by the analysis initial conditions.Using the adjoint of the entire 4D-Var system we have also explored the sensitivity of the coastal transport to changes in the observations and the observation array. A single integration of the adjoint of 4D-Var can be used to predict the change that occurs when observations from different platforms are omitted from the 4D-Var analysis. Thus observing system experiments can be performed for each data assimilation cycle at a fraction of the computational cost that would be required to repeat the 4D-Var analyses when observations are withheld. This is the third part of a three part series describing the ROMS 4D-Var systems.     

Cycling the Representer Method for 4D-variational data assimilation with the Navy Coastal Ocean Model

The Cycling Representer Method, which is a technique for solving 4D-variational data assimilation problems, has been demonstrated to improve the assimilation accuracy with simpler nonlinear models. In this paper, the Cycling Representer Method will be used to assimilate an array of ADCP velocity observations with the Navy Coastal Ocean Model (NCOM). Experiments are performed in a high-resolution Mississippi Bight domain for the entire month of June, 2004 and demonstrate the usefulness of this assimilation technique in a realistic application.The Representer Method is solved by minimizing a cost function containing the weighted squared errors of velocity measurements, initial conditions, boundary conditions, and model dynamics. NCOM, however, is a highly nonlinear model and in order to converge towards the global minimum of this cost function, NCOM is linearized about a background state using tangent linearization. The stability of this tangent linearized model (TLM) is a very sensitive function of the background state, the level of nonlinearity of the model, open boundary conditions, and the complexity of the bathymetry and flow field. For the Mississippi Bight domain, the TLM is stable for only about a day. Due to this short TLM stability time period, the Representer Method is cycled by splitting the time period of the assimilation problem into short intervals. The interval time period needs to be such that it is short enough for the TLM to be stable, but long enough to minimize the loss of information due to reducing the temporal correlation of the dynamics and data. For each new cycle, a background is created as a nonlinear forecast from the previous cycle’s assimilated solution. This background, along with the data that falls within this new cycle, is then used to calculate a new assimilated solution. The experiments presented in this paper demonstrate the improvement of the assimilated solution as the time window of the cycles is reduced to 1 day. The 1-day cycling, however, was only optimal for the first half of the experiment. This is because there was a strong wind event near the middle of June that significantly reduced the stability of the 1-day cycling and caused substantial errors in the assimilation. Therefore, the 12-h cycling worked best for the second half of the experiment. This paper also demonstrates that the forecast skill is improved as the assimilation system progresses through the cycles.     

Weak constraint data assimilation for tides in the Arctic Ocean

Estimation of the open-boundary inputs solving a weak constraint variational formulation for an Arctic tide model is considered as an ill-posed problem in the sense that the solution is very sensitive to the data noise and to grid size. Mathematically, spatial discretization of a cost function to be minimized and penalization of normal flow through the open boundary act as regularization of the problem. An heuristic choosing rule for the regularization parameter is applied to assess a suitable spatial resolution and the weight referred to the open boundary penalty. It is shown that these provide a better fit of the solution to a control data set compared with a finer grid, the value of the energy flux through the open boundary being in agreement with other model estimates. The M₂ solution obtained is much closer to the control data than other modern solutions while the accuracy of the simulated K₁ constituent is within the same error level. The tidal maps for these waves exhibit certain distinctions in comparison with other charts.     

Impact on oceanic dynamics from assimilation of satellite surface height anomaly data into the Hybrid Coordinate Ocean Model Ocean Model (HYCOM) over the Atlantic Ocean

Observed along-track data of sea surface height anomalies (SSHAs) over the Atlantic Ocean from the Jason-1 and Jason-2 satellites were assimilated into the Hybrid-Coordinate Ocean Model (HYCOM) with the Ensemble Optimal Interpolation scheme (EnOI). The impact of assimilation of SSHA with focus on oceanic dynamics was investigated. Time series of analyzed and forecasted values were compared with a model free run with the same forcing but without assimilation. In addition, the results were compared with an independent run, the so-called HYCOM + NCODA analysis from the US Navy. The study shows that the assimilation technique with some modifications allowed substantial improvement in the 24 h ocean prediction by reducing the forecast errors in comparison with the free run. It is also shown that the analyzed sea surface fields contain mesoscale and synoptic variability, which are poorly seen in the free run.     

Detecting massive green algae (Ulva prolifera) blooms in the Yellow Sea and East China Sea using Geostationary Ocean Color Imager (GOCI) data

The historically massive bloom of the green macroalgae Ulva prolifera reported in June?CAugust 2008 around the Qingdao, Yellow Sea, East China Sea and Japan coasts has recurred in a similar season and region. On June 13, 2011, around Qingdao, China, the world??s first Geostationary Ocean Color Imager (GOCI) detected an enormous bloom of floating green algae, which originated from the nearshore Subei Bank, China. The large floating green algae patches were observed along and across the Yellow Sea and in the East China Sea during 2011 summer by various oceanic cruises. To detect the massive macroalgae blooms from space, we analyzed their spectral characteristics from in situ optical measurements and satellite-derived green algae spectra. An ??Index of floating Green Algae for GOCI?? (IGAG) was developed from the multiple spectral band ratios using three wavelengths (555, 660, 745 nm), which the spectral response of green algae reflected at 555, 745, and 865 nm and absorbed at 660 and 680 nm. The results were compared with those obtained by the normalized difference vegetation index (NDVI), enhanced vegetation index (EVI), and Korea Ocean Satellite Center (KOSC) approaches. An advantage of the IGAG method was that muted or subtle signals of floating green algae were enhanced and separated from surrounding complex water signals. Although maps of floating green algae derived by the other approaches delineated dense green algae, they were less sensitive to subtle (less dense) features and in cases of nearby cloudy or complex water conditions. The floating green algae maps from IGAG provided a more robust estimate of wide floating green algae blooms than those derived using NDVI, EVI, or KOSC approaches. The IGAG approach should be useful for tracing and monitoring changes in green algae blooms on regional and global scales.     

俄罗斯北极海(喀拉海)的生物地球化学

据1995—2003年SIRRO计划的研究成果,喀拉海是研究河流—海洋体系相互作用过程独一无二的地区。巨大的西伯利亚叶尼塞河和鄂毕河注入这个浅海。1995—2003年科学研究船“鲍利斯.彼得罗夫”号对喀拉海进行了国际性考察。提供的工作成果总结了俄罗斯科学院地球化学与分析化学研究所完成的研究成果。对整个海区200多个测站研究了沉积物的水化学参数,有机碳和碳酸盐碳的含量及同位素成分,水中浮游生物和悬浮物质以及烃类和溶解CO2。在大西洋水进入喀拉海区δ13C有机的变化范围为-22‰~-24‰,而在东北海区叶尼塞和鄂毕河口区则为-27‰~-30‰…     

Spectral wave data assimilation for the prediction of waves in the North Sea

Spectral observations from pitch-and-roll buoys have been assimilated in a North Sea wave model, in order to study their impact on the wave analysis and forecast. The assimilation is based on Optimal Interpolation (OI) of a limited number of characteristic spectral parameters. In a case study, the propagation of the corrections through the model domain is followed, and it is clarified for which wave conditions the data assimilation has the largest influence on the forecast: this is especially the case for swell waves with long travel times between the assimilation site and the location where validation is carried out. A 1-year test has been carried out in which an analysis and subsequent forecast were produced four times a day. From a statistical analysis of the results a modest but systematic improvement of the 12-h forecast is found. When only swell cases are selected, the impact is more pronounced. It is argued that for shelf seas like the North Sea, more progress is to be expected from extension of the ‘conventional' observations network (buoys and wave radars) than from satellite measurements.     

An ocean data assimilation system and reanalysis of the World Ocean hydrophysical fields

A new version of the ocean data assimilation system (ODAS) developed at the Hydrometcentre of Russia is presented. The assimilation is performed following the sequential scheme analysis–forecast–analysis. The main components of the ODAS are procedures for operational observation data processing, a variational analysis scheme, and an ocean general circulation model used to estimate the first guess fields involved in the analysis. In situ observations of temperature and salinity in the upper 1400-m ocean layer obtained from various observational platforms are used as input data. In the new ODAS version, the horizontal resolution of the assimilating model and of the output products is increased, the previous 2D-Var analysis scheme is replaced by a more general 3D-Var scheme, and a more flexible incremental analysis updating procedure is introduced to correct the model calculations. A reanalysis of the main World Ocean hydrophysical fields over the 2005–2015 period has been performed using the updated ODAS. The reanalysis results are compared with data from independent sources.     

El Ni(n)o phenomenon in simple ocean data assimilation data

To study how the air and sea interact with each other during El Nino/La Nina onsets, extended associate pattern analysis (EAPA) is adopted with the simple ocean data assimilation (SODA) data. The results show that as El Nino/La Nina' s parents their behaviors are quite different, there does not exist a relatively independent tropical atmosphere but does exist a relatively independent tropical Pacific Ocean because the air is heated from the bottom surface instead of the top surface and of much stronger baroclinic instability than the sea and has a very large inter-tropical convergence zone covering the most tropical Pacific Ocean. The idea that it is the wester burst and wind convergence, coming from middle latitudes directly that produce the seawater eastward movement and meridional convergence in the upper levels and result in the typical El Nino sea surface temperature warm signal is confirmed again.     

Asynchronous data assimilation with the EnKF

This study revisits the problem of assimilation of asynchronous observations, or four-dimensional data assimilation, with the ensemble Kalman filter (EnKF). We show that for a system with perfect model and linear dynamics the ensemble Kalman smoother (EnKS) provides a simple and efficient solution for the problem: one just needs to use the ensemble observations (that is, the forecast observations for each ensemble member) from the time of observation during the update, for each assimilated observation. This recipe can be used for assimilating both past and future data; in the context of assimilating generic asynchronous observations we refer to it as the asynchronous EnKF. The asynchronous EnKF is essentially equivalent to the four-dimensional variational data assimilation (4D-Var). It requires only one forward integration of the system to obtain and store the data necessary for the analysis, and therefore is feasible for large-scale applications. Unlike 4D-Var, the asynchronous EnKF requires no tangent linear or adjoint model.     

Benthic foraminifera communities of the Andaman Sea (Indian Ocean)

Eleven communities of benthic foraminifera were distinguished on the basis of the results of a qualitative analysis of the structure the population at 94 stations at sea depths down to 4180 m in the Andaman Sea. The habitats of the communities are formed under the influence of the environmental factors, which are determined by the latitudinal and bathymetric zonality of the basin, the water masses, the currents, the upwellings, the coastal runoff, and the contents of calcium carbonate and oxygen in the near-bottom waters and organic carbon in the sediments.