Impact of Ocean Mean Dynamic Topography on Satellite Data Assimilation

The response of an eddy-permitting ocean model to changes imposed by the use of different mean dynamic topographies (MDT) is analyzed in a multivariate assimilation context, allowing the evaluation of this impact, not only on the surface circulation, but also on the interior ocean representation. The assimilation scheme is a reduced-order sequential Kalman filter (SEEK). In a first set of experiments, high resolution sea surface temperature, along-track sea surface height and sea surface salinity from climatology are assimilated into a 1/3° resolution North and Tropical Atlantic version of the HYCOM model. In a second experiment, in situ profile data are assimilated in addition to the surface measurements.

The first set of experiments illustrates that important differences in the representation of the horizontal model circulation pattern are related to differences in the MDT used. The objective of assimilation is to improve the representation of the 3D ocean state. However, the imperfect representation of the mean dynamic topography appears to be an important limiting factor with regard to the degree of realism obtained in the simulated flow.

Vertical temperature and salinity profiles are key observations to drive a general circulation ocean model toward a more realistic state. The second set of experiments shows that assimilating them in addition to sea surface measurements is a far from trivial exercise. A specific difficulty is due to inconsistencies between the dynamic topography diagnosed from in situ observations and that diagnosed from sea surface height. These two fields obtained from different data sources do not contain exactly the same information. In order to overcome this difficulty, a strategy is proposed and validated.     

Improving the Coastal Mean Dynamic Topography by Geodetic Combination of Tide Gauge and Satellite Altimetry

Abstract

The ocean mean dynamic topography (MDT) is the surface representation of the ocean circulation. The MDT may be determined by the ocean approach, which involves temporal averaging of numerical ocean circulation model information, or by the geodetic approach, wherein the MDT is derived using the ellipsoidal height of the mean sea surface (MSS), or mean sea level (MSL) minus the geoid as the geoid. The ellipsoidal height of the MSS might be estimated either by satellite or coastal tide gauges by connecting the tide gauge datum to the Earth-centred reference frame. In this article we present a novel approach to improve the coastal MDT, where the solution is based on both satellite altimetry and tide gauge data using new set of 302 tide gauges with ellipsoidal heights through the SONEL network. The approach was evaluated for the Northeast Atlantic coast where a dense network of GNSS-surveyed tide gauges is available. The typical misfit between tide gauge and satellite or oceanographic MDT was found to be around 9?cm. This misfit was found to be mainly due to small scale geoid errors. Similarly, we found, that a single tide gauge places only weak constraints on the coastal dynamic topography.     

中国海及邻近海域卫星观测资料同化试验

利用1个基于POMgcs海洋模式和多重网格三维变分同化方法建立的中国海及邻近海域海面高与三维温盐流数值预报模型,通过一系列数值试验,研究了同化卫星测高和卫星遥感海面温度观测资料对该模型预报能力的影响。试验结果表明,同化卫星测高资料可明显改善海面高度与三维温度和盐度的分析预报效果,使1 200 m以上的温度预报误差减小0.16℃,并能有效提高对海洋中尺度现象的预报能力;同化卫星遥感海面温度对100 m以上的温度和盐度的预报效果有所改善,可使海面温度的预报误差减小10%。     

卫星高度计资料在三维海温和盐度数值预报中的应用

随着卫星遥感观测技术的发展,越来越多的卫星观测资料被应用于数值模式的同化研究中.基于国家海洋环境预报中心西北太平洋三维湿盐流预报系统,利用法国CLS中心的沿轨高度计资料的海表面高度异常的融合数据,结合基于三维变分的OVALS(ocean variational analysis system)同化系统,在垂向将海面高度...     

Improving Ocean State by Assimilating SARAL/AltiKa Derived Sea Level and Other Satellite-Derived Data in MITGCM

Assimilation of satellite-derived surface datasets has been explored in the study. Three types of surface data, namely sea level anomaly, sea surface temperature and sea surface salinity, have been used in various data assimilation experiments. The emphasis has been on the extra benefit arising out of the additional sea level assimilation and hence there are two parallel runs, in one of which sea level assimilation has been withheld. The model used is a state-of-the art ocean general circulation model (OGCM) and the assimilation method is the widely used singular evolutive extended Kalman filter (SEEK). Evaluation of the assimilation skill has been carried out by comparing the simulated depth of the 20°C isotherm with the same quantity measured by buoys and Argo floats. Simulated subsurface temperature and salinity profiles have also been compared with the same profiles measured by Argo floats. Finally, surface currents in the assimilation runs have been compared with currents measured by several off-equatorial buoys. Addition of sea level has been found to substantially improve the quality of simulation. An important feature that has been effectively simulated by the addition of sea level in the assimilation scheme is the near-surface temperature inversion (2-3°C) in the northern Bay of Bengal.     

Upper Ocean Four-Dimensional Variational Data Assimilation in the Arabian Sea and Bay of Bengal

A regional ocean circulation model with four-dimensional variational data assimilation scheme is configured to study the ocean state of the Indian Ocean region (65°E–95°E; 5°N–20°N) covering the Arabian Sea (AS) and Bay of Bengal (BoB). The state estimation setup uses 10 km horizontal resolution and 5 m vertical resolution in the upper ocean. The in-situ temperature and salinity, satellite-derived observations of sea surface height, and blended (in-situ and satellite-derived) observations of sea surface temperature alongwith their associated uncertainties are used for data assimilation with the regionally configured ocean model. The ocean state estimation is carried out for 61 days (1 June to 31 July 2013). The assimilated fields are closer to observations compared to other global state estimates. The mixed layer depth (MLD) of the region shows deepening during the period of assimilation with AS showing higher MLD compared to the BoB. An empirical forecast equation is derived for the prediction of MLD using the air–sea forcing variables as predictors. The surface and sub-surface (50 m) heat and salt budget tendencies of the region are also investigated. It is found that at the sub-surface, only the advection and diffusion temperature and salt tendencies are important.     

气候模式中海洋数据同化对热带降水偏差的影响

本文采用海洋卫星观测海表温度(SST)和海面高度异常(SLA)数据,对国家海洋局第一海洋研究所地球系统模式FIO-ESM(First Institute of Oceanography Earth System Model version 1.0)中海洋模式分量进行了集合调整卡尔曼滤波(EAKF)同化,对比分析了大气环流、湿度和云量对海洋数据同化的响应,探讨了海洋同化对热带降水模拟偏差的影响。结果表明:海洋数据同化能有效改善海表温度和上层海洋热含量的模拟,30°S~30°N纬度带内年平均SST的绝均差降低60%。同化后大气模式模拟的赤道两侧信风得到明显改善,上升气流在赤道以北热带地区增强而在赤道以南热带地区减弱,热带降水模拟的动力结构更为合理,水汽和云量分布也更切合实际。热带年平均降水的空间分布和强度在同化后均得到改善,赤道以南的纬向年平均降水峰值显著降低,降水偏差明显减小,同化后30°S~30°N纬度带内年平均降水绝均差降低35%。     

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.     

Impact of Aquarius Sea-Surface Salinity Assimilation in Improving the Ocean Analysis Over Indian Ocean

The distribution of ocean salinity controls the density field and thereby plays a major role in influencing the ocean dynamics. It has been a challenging task to understand the variability of salinity structure in the regions of large fresh water discharge and high precipitation such as Bay of Bengal (BoB). Recent advancement in satellite technology has made possible the measurement of sea surface salinity (SSS). Aquarius is the satellite which measured the global SSS for the period 2011 to 2015. In the present study, we assimilated Aquarius SSS in the Global Ocean Data Assimilation System based on 3DVAR technique. The assimilation of Aquarius SSS resulted in reduced biases in salinity not only at the surface, but also in the vertical distribution of salinity and better captured the temporal variations of salinity structure in sensitive regions, such as the Bay of Bengal. In addition, the assimilation of SSS showed marginal improvement in ocean thermal structure over data sparse regions of Indian Ocean. It is also shown that the assimilation of Aquarius SSS has improved the stratification in the upper Ocean which is the key factor in the observed improvement in ocean analysis.     

The impact of ocean data assimilation on seasonal predictions based on the National Climate Center climate system model

An ensemble optimal interpolation (EnOI) data assimilation method is applied in the BCC_CSM1.1 to investigate the impact of ocean data assimilations on seasonal forecasts in an idealized twin experiment framework. Pseudo-observations of sea surface temperature (SST), sea surface height (SSH), sea surface salinity (SSS), temperature and salinity (T/S) profiles were first generated in a free model run. Then, a series of sensitivity tests initialized with predefined bias were conducted for a one-year period; this involved a free run (CTR) and seven assimilation runs. These tests allowed us to check the analysis field accuracy against the “truth”. As expected, data assimilation improved all investigated quantities; the joint assimilation of all variables gave more improved results than assimilating them separately. One-year predictions initialized from the seven runs and CTR were then conducted and compared. The forecasts initialized from joint assimilation of surface data produced comparable SST root mean square errors to that from assimilation of T/S profiles, but the assimilation of T/S profiles is crucial to reduce subsurface deficiencies. The ocean surface currents in the tropics were better predicted when initial conditions produced by assimilating T/S profiles, while surface data assimilation became more important at higher latitudes, particularly near the western boundary currents. The predictions of ocean heat content and mixed layer depth are significantly improved initialized from the joint assimilation of all the variables. Finally, a central Pacific El Ni?o was well predicted from the joint assimilation of surface data, indicating the importance of joint assimilation of SST, SSH, and SSS for ENSO predictions.     

OSTIA数据在中国近海业务化环流模型中的同化应用

The prediction of sea surface temperature(SST) is an essential task for an operational ocean circulation model. A sea surface heat flux, an initial temperature field, and boundary conditions directly affect the accuracy of a SST simulation. Here two quick and convenient data assimilation methods are employed to improve the SST simulation in the domain of the Bohai Sea, the Yellow Sea and the East China Sea(BYECS). One is based on a surface net heat flux correction, named as Qcorrection(QC), which nudges the flux correction to the model equation; the other is ensemble optimal interpolation(En OI), which optimizes the model initial field. Based on such two methods, the SST data obtained from the operational SST and sea ice analysis(OSTIA) system are assimilated into an operational circulation model for the coastal seas of China. The results of the simulated SST based on four experiments, in 2011, have been analyzed. By comparing with the OSTIA SST, the domain averaged root mean square error(RMSE) of the four experiments is 1.74, 1.16, 1.30 and 0.91°C, respectively; the improvements of assimilation experiments Exps 2, 3 and 4 are about 33.3%, 25.3%, and 47.7%, respectively.Although both two methods are effective in assimilating the SST, the En OI shows more advantages than the QC,and the best result is achieved when the two methods are combined. Comparing with the observational data from coastal buoy stations, show that assimilating the high-resolution satellite SST products can effectively improve the SST prediction skill in coastal regions.     

Assimilation of sea surface temperature in the Northwest Pacific Ocean and its marginal seas using the ensemble Kalman filter

Satellite-borne sea surface temperature (SST) data were assimilated with the ensemble Kalman filter (EnKF) in a Northwest Pacific Ocean circulation model to examine the effect of data assimilation. The model domain included the northwestern part of the Pacific Ocean and its marginal seas, such as the Yellow Sea and East/Japan Sea. The performance of the data assimilation was evaluated by comparing the simulated ocean state with that observed. Spatially averaged root-mean-squared errors in the SST and sea surface height (SSH) decreased by 0.44 °C and 4 cm, respectively, by the assimilation. The results of the numerical experiments substantiated the effectiveness of the SST assimilation via the EnKF for all marginal seas, as well as the Kuroshio region. The benefit of the data assimilation depended on the characteristics of each marginal sea. The variation of the SST in the East/Japan Sea and the Kuroshio extension (KE) region were improved 34% and those in the Yellow Sea 12.5%. The variation of the SSH was improved approximately 36% in the KE region. This large improvement was achieved in the deep-water regions because assimilation of SST data corrected the separation point of the western boundary currents, such as the Kuroshio and the East Korea Warm Current, and the associated horizontal surface currents. The SST assimilation via the EnKF also improved the subsurface temperature profiles. The effectiveness of SST assimilation was seasonally dependent, with the improvement being relatively larger in winter than in summer, which was related to the seasonal variation of the vertical mixing and stratification in the ocean surface layer.     

南海北部海洋四维变分同化试验研究

为了研究四维变分同化方法在南海北部海洋数值预报中的适用性,使用海洋区域模式(ROMS),建立了南海北部海洋资料四维变分同化系统,进行了温盐廓线和海面温度数据同化试验,初步对比分析了三种四维变分实现方法的同化效果。研究结果表明,四维变分同化方法具有较好的同化效果,其中,增量强约束方法(I4DVar)具有较好的稳定性,其稳定性高于4DPSAS和R4DVar。本文研究成果为建立南海业务化海洋四维变分同化及预报系统奠定技术基础。     

The impact of mean dynamic topography on a sea-level anomaly assimilation in the South China Sea based on an eddy-resolving model

The sea-level anomaly (SLA) from a satellite altimeter has a high accuracy and can be used to improve ocean state estimation by assimilation techniques. However, the lack of an accurate mean dynamic topography (MDT) is still a bothersome issue in an ocean data assimilation. The previous studies showed that the errors in MDT have significant impacts on assimilation results, especially on the time-mean components of ocean states and on the time variant parts of states via nonlinear ocean dynamics. The temporal-spatial differences of three MDTs and their impacts on the SLA analysis are focused on in the South China Sea (SCS). The theoretical analysis shows that even for linear models, the errors in MDT have impacts on the SLA analysis using a sequential data assimilation scheme. Assimilation experiments, based on EnOI scheme and HYCOM, with three MDTs from July 2003 to June 2004 also show that the SLA assimilation is very sensitive to the choice of different MDTs in the SCS with obvious differences between the experimental results and observations in the centre of the SCS and in the vicinity of the Philippine Islands. A new MDT for assimilation of SLA data in the SCS was proposed. The results from the assimilation experiment with this new MDT show a marked reduction (increase) in the RMSEs (correlation coefficient) between the experimental and observed SLA. Furthermore, the subsurface temperature field is also improved with this new MDT in the SCS.     

State Estimation of the North Pacific Ocean by a Four-Dimensional Variational Data Assimilation Experiment

A four-dimensional variational data assimilation system has been applied to an experiment to describe the dynamic state of the North Pacific Ocean. A synthesis of available observational records and a sophisticated ocean general circulation model produces a dynamically consistent dataset, which, in contrast to the nudging approach, provides realistic features of the seasonally-varying ocean circulation with no artificial sources/sinks for temperature and salinity fields. This new dataset enables us to estimate heat and water mass transports in addition to the qualification of water mass formation and movement processes. A sensitivity experiment on our assimilation system reveals that the origin of the North Pacific Intermediate Water can be traced back to the Sea of Okhotsk and the Bering Sea in the subarctic region and to the subtropical Kuroshio region further south. These results demonstrate that our data assimilation system is a very powerful tool for the identification and characterization of ocean variabilities and for our understanding of the dynamic state of ocean circulation. This revised version was published online in July 2006 with corrections to the Cover Date.     

A global coupled sea ice–ocean model

A new sea ice model, GELATO, was developed at Centre National de Recherches Météorologiques (CNRM) and coupled with OPA global ocean model. The sea ice model includes elastic–viscous–plastic rheology, redistribution of ice floes of different thicknesses, and it also takes into account leads, snow cover and snow ice formation. Climatologies of atmospheric surface parameters are used to perform a 20-year global ocean–sea ice simulation, in order to compute surface heat fluxes from diagnosed sea ice or ocean surface temperature. A surface salinity restoring term is applied only to ocean grid cells with no sea ice to avoid significant surface salinity drifts, but no correction of sea surface temperature is introduced. In the Arctic the use of an ocean model substantially improves the representation of sea ice, and particularly of the ice edge in all seasons, as advection of heat and salt can be more accurately accounted for than in the case of, for example, a sea ice–ocean mixed layer model. In contrast, in the Antarctic, a region where ocean convective processes bear a much stronger influence in shaping sea ice characteristics, a better representation of convection and probably of sea ice (for example, of frazil sea ice, brine rejection) would be needed to improve the simulation of the annual cycle of the sea ice cover. The effect of the inclusion of several ice categories in the sea ice model is assessed by running a sensitivity experiment in which only one category of sea ice is considered, along with leads. In the Arctic, such an experiment clearly shows that a multicategory sea ice model better captures the position of the sea ice edge and yields much more realistic sea ice concentrations in most of the region, which is in agreement with results from Bitz et al. [J. Geophys. Res. 106 (C2) (2001) 2441–2463].     

Corrections to ocean surface forcing in the California Current System using 4D variational data assimilation

The option for surface forcing correction, recently developed in the 4D-variational (4DVAR) data assimilation systems of the Regional Ocean Model System (ROMS), is presented. Assimilation of remotely-sensed (satellite sea surface height anomaly and sea surface temperature) and in situ (from mechanical and expendable bathythermographs, Argo floats and CTD profiles) oceanic observations has been applied in a realistic, high resolution configuration of the California Current System (CCS) to sequentially correct model initial conditions and surface forcing, using the Incremental Strong constraint version of ROMS-4DVAR (ROMS-IS4DVAR). Results from both twin and real data experiments are presented where it is demonstrated that ROMS-IS4DVAR always reduces the difference between the model and the observations that are assimilated. However, without corrections to the surface forcing, the assimilation of surface data can degrade the temperature structure at depth. When using surface forcing adjustment in ROMS-IS4DVAR the system does not degrade the temperature structure at depth, because differences between the model and surface observations can be reduced through corrections to surface forcing rather than to temperature at depth. However, corrections to surface forcing can generate abnormal spatial and temporal variability in the structure of the wind stress or surface heat flux fields if not properly constrained. This behavior can be partially controlled via the choice of decorrelation length scales that are assumed for the forcing errors. Abnormal forcing corrections may also arise due to the effects of model error which are not accounted for in IS4DVAR. In particular, data assimilation tends to weaken the alongshore wind stress in an attempt to reduce the rate of coastal upwelling, which seems to be too strong due to other sources of error. However, corrections to wind stress and surface heat flux improve systematically the ocean state analyses. Trends in the correction of surface heat fluxes indicate that, given the ocean model used and its potential limitations, the heat flux data from the Coupled Ocean–Atmosphere Mesoscale Prediction System (COAMPS) used to impose surface conditions in the model are generally too low except in spring-summer, in the upwelling region, where they are too high. Comparisons with independent data provide confidence in the resulting forecast ocean circulation on timescales ～14 days, with less than 1.5 °C, 0.3 psu, and 9 cm RMS error in temperature, salinity and sea surface height anomaly, respectively, compared to observations.     

基于四维变分同化方法的南海中尺度涡后报实验

海洋中尺度涡在本质上是属于满足准地转平衡的大尺度运动,因此理论上,其在短时间内的运动将主要受到准地转平衡关系的约束,而外部强迫场的影响在短期内不会明显改变其运动特征。基于上述思想,我们提出了一种基于四维变分同化初始场的中尺度涡旋预报方案。为了检验该方案的可行性,本文使用区域海洋模式(regional ocean modeling system, ROMS)和其内建的增量强约束四维变分同化(incremental strong constraint four dimensional variational, I4D-Var)模块,建立了一个南海海洋同化模拟系统。首先,通过I4D-Var方法将AVISO卫星高度计资料同化到海洋数值模拟中,获得了理想的中尺度涡同化模拟结果。同化、模式模拟和观测三者的中尺度涡统计结果表明,该同化系统模拟的南海中尺度涡的路径、半径、海表高度异常和振幅等特征信息与AVISO(Archiving ValidationandInterpolationofSatelliteOceanographicData)观测结果高度吻合,同时在深度上的分析表明,涡旋对应的温度、盐度和密度均得到有效的调整。然后,将该同化系统的模拟结果做为初始场,对某一特定时段的南海中尺度涡进行了后报模拟和结果的定量化分析。通过比较后报模拟与观测资料中对应涡旋的海表面高度异常(sea surface height anomalies, SSHA)相关系数、涡心差距和半径绝对误差,证明该方案的中尺度涡后报时效至少可达10 d以上。后报实验结果验证了该中尺度涡预报方案的可行性,从而为中尺度涡的预报提供一定的理论基础和可行性方案。     

模块化海洋数据同化系统(MODAS)研究

模块化海洋数据同化系统(Modular Ocean Data Assimilation System,MODAS)是美国海军用来确定全球海洋三维温度和盐度场的主要工具。MODAS通过同化卫星遥感测得的海面温度和海面高度,产生一种动态气候态,能够更接近地预报出海洋的真实状况。文中介绍了MODAS的应用价值和构建原理,对中国海军数字化海洋战场环境的建设具有重要的借鉴作用。