基于ROMS模式利用集合最优插值同化沿轨海表面高度异常数据的研究

The ensemble optimal interpolation （EnOI） is applied to the regional ocean modeling system （ROMS） with the ability to assimilate the along-track sea level anomaly （TSLA）. This system is tested with an eddy-resolving system of the South China Sea （SCS）. Background errors are derived from a running seasonal ensemble to account for the seasonal variability within the SCS. A fifth-order localization function with a 250 km localization radius is chosen to reduce the negative effects of sampling errors. The data assimilation system is tested from January 2004 to December 2006. The results show that the root mean square deviation （RMSD） of the sea level anomaly decreased from 10.57 to 6.70 cm, which represents a 36.6% reduction of error. The data assimilation reduces error for temperature within the upper 800 m and for salinity within the upper 200 m, although error degrades slightly at deeper depths. Surface currents are in better agreement with trajectories of surface drifters after data assimilation. The variance of sea level improves significantly in terms of both the amplitude and position of the strong and weak variance regions after assimilating TSLA. Results with AGE error （AGE） perform better than no AGE error （NoAGE） when considering the improvements of the temperature and the salinity. Furthermore, reasons for the extremely strong variability in the northern SCS in high resolution models are investigated. The results demonstrate that the strong variability of sea level in the high resolution model is caused by an extremely strong Kuroshio intrusion. Therefore, it is demonstrated that it is necessary to assimilate the TSLA in order to better simulate the SCS with high resolution models.     

The effect of the wave-induced mixing on the upper ocean temperature in a climate model

The significant underestimation of sea surface temperature （SST） and the temperature in the upper ocean is one of common problems in present climate models. The influence of the wave-induced mixing on SST and the temperature in the upper ocean was examined based on a global climate model. The results from the model coupled with wave-induced mixing showed a significant improvement in the simulation of SST and the temperature in the upper ocean compared with those of the original model without wave effects. Although there has still a cold bias, the new simulation is much closer to the climatology, especially in the northern ocean and tropical ocean. This study indicates that some important physical processes in the accurate simulation of the ocean may be ignored in present climate models, and the wave-induced mixing is one of those factors. Thus, the wave-induced mixing （ or the effect of surface waves） should be incorporated properly into climate models in order to simulate or forecast the ocean, then climate system, more accurately.     

基于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 effects of sea surface temperature(SST) data assimilation in two regional ocean modeling systems were examined for the Yellow Sea(YS). The SST data from the Operational Sea Surface Temperature and Sea Ice Analysis(OSTIA) were assimilated. The National Marine Environmental Forecasting Center(NMEFC) modeling system uses the ensemble optimal interpolation method for ocean data assimilation and the Kunsan National University(KNU) modeling system uses the ensemble Kalman filter. Without data assimilation, the NMEFC modeling system was better in simulating the subsurface temperature while the KNU modeling system was better in simulating SST. The disparity between both modeling systems might be related to differences in calculating the surface heat flux, horizontal grid spacing, and atmospheric forcing data. The data assimilation reduced the root mean square error(RMSE) of the SST from 1.78°C(1.46°C) to 1.30°C(1.21°C) for the NMEFC(KNU) modeling system when the simulated temperature was compared to Optimum Interpolation Sea Surface Temperature(OISST) SST dataset. A comparison with the buoy SST data indicated a 41%(31%) decrease in the SST error for the NMEFC(KNU) modeling system by the data assimilation. In both data assimilative systems, the RMSE of the temperature was less than 1.5°C in the upper 20 m and approximately 3.1°C in the lower layer in October. In contrast, it was less than 1.0°C throughout the water column in February. This study suggests that assimilations of the observed temperature profiles are necessary in order to correct the lower layer temperature during the stratified season and an ocean modeling system with small grid spacing and optimal data assimilation method is preferable to ensure accurate predictions of the coastal ocean in the YS.     

引潮力对海洋环流模式的影响

The eight main tidal constituents have been implemented in the global ocean general circulation model with approximate 1° horizontal resolution.Compared with the observation data,the patterns of the tidal amplitudes and phases had been simulated fairly well.The responses of mean circulation,temperature and salinity are further investigated in the global sense.When implementing the tidal forcing,wind-driven circulations are reduced,especially those in coastal regions.It is also found that the upper cell transport of the Atlantic meridional overturning circulation(AMOC) reduces significantly,while its deep cell transport is slightly enhanced from 9×106m3/s to 10×106 m3/s.The changes of circulations are all related to the increase of a bottom friction and a vertical viscosity due to the tidal forcing.The temperature and salinity of the model are also significantly affected by the tidal forcing through the enhanced bottom friction,mixing and the changes in mean circulation.The largest changes occur in the coastal regions,where the water is cooled and freshened.In the open ocean,the changes are divided into three layers:cooled and freshened on the surface and below 3 000 m,and warmed and salted in the middle in the open ocean.In the upper two layers,the changes are mainly caused by the enhanced mixing,as warm and salty water sinks and cold and fresh water rises;whereas in the deep layer,the enhancement of the deep overturning circulation accounts for the cold and fresh changes in the deep ocean.     

Reconstruction of eddies by assimilating satellite altimeter data into Princeton Ocean Model

An optimal interpolation assimilation model for satellite altimetry data is developed based on Princeton Ocean Model (POM), which is applied in a quasi-global domain, by the method of isotropic correlation between sea level anomaly (SLA) and sea temperature anomaly. The perfor- mance of this assimilation model is validated by the modeled results of SLA and the current patterns. Comparisons between modeling and satellite data show that both the magnitudes and distribution patterns of the simulated SLA are improved by assimilation. The most signiˉcant improvement is that meso-scale systems, e.g., eddies, are well reconstructed. The evolution of an eddy located in the northwest Paciˉc Ocean is traced by using the assimilation model. Model results show that during three months the eddy migrated southwestward for about 6 degrees before merging into the Kuroshio. The three dimensional structure of this eddy on 12 August 2001 is further analyzed. The strength of this warm, cyclonic eddy decreases with the increase of depth. The eddy shows di?erent horizontal patterns at di?erent layers, and the SLA and temperature ˉelds agree with each other well. This study suggests that this kind of data assimilation is economic and reliable for eddy reconstruction, and can be used as a promising technique in further studies of ocean eddies as well as other ˉne circulation structures.     

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.     

Three-dimensional structure of an observed cyclonic mesoscale eddy in the Northwest Pacific and its assimilation experiment

Mesoscale eddies play an important role in modulating the ocean circulation. Many previous studies on the threedimensional structure of mesoscale eddies were mainly based on composite analysis, and there are few targeted observations for individual eddies. A cyclonic eddy surveyed during an oceanographic cruise in the Northwest Pacific Ocean is investigated in this study. The three-dimensional structure of this cyclonic eddy is revealed by observations and simulated by the four-dimensional variational data assimilation(4 DVAR) system combined with the Regional Ocean Modeling System. The observation and assimilation results together present the characteristics of the cyclonic eddy. The cold eddy has an obvious dual-core structure of temperature anomaly.One core is at 50–150 m and another is at 300–550 m, which both have the average temperature anomaly of approximately-3.5°C. The salinity anomaly core is between 250 m and 500 m, which is approximately-0.3. The horizontal velocity structure is axis-asymmetric and it is enhanced on the eastern side of the cold eddy. In the assimilation experiment, sea level anomaly, sea surface temperature, and in situ measurements are assimilated into the system, and the results of assimilation are close to the observations. Based on the high-resolution assimilation output results, the study also diagnoses the vertical velocity in the mesoscale eddy, which reaches the maximum of approximately 10 m/d. The larger vertical velocity is found to be distributed in the range of 0.5 to 1 time of the normalized radius of the eddy. The validation of the simulation result shows that the 4 DVAR method is effective to reconstruct the three-dimensional structure of mesoscale eddy and the research is an application to study the mesoscale eddy in the Northwest Pacific by combining observation and assimilation methods.     

Upper ocean high resolution regional modeling of the Arabian Sea and Bay of Bengal

In this paper, effort is made to demonstrate the quality of high-resolution regional ocean circulation model in realistically simulating the circulation and variability properties of the northern Indian Ocean(10°S–25°N,45°–100°E) covering the Arabian Sea(AS) and Bay of Bengal(BoB). The model run using the open boundary conditions is carried out at 10 km horizontal resolution and highest vertical resolution of 2 m in the upper ocean.The surface and sub-surface structure of hydrographic variables(temperature and salinity) and currents is compared against the observations during 1998–2014(17 years). In particular, the seasonal variability of the sea surface temperature, sea surface salinity, and surface currents over the model domain is studied. The highresolution model's ability in correct estimation of the spatio-temporal mixed layer depth(MLD) variability of the AS and BoB is also shown. The lowest MLD values are observed during spring(March-April-May) and highest during winter(December-January-February) seasons. The maximum MLD in the AS(BoB) during December to February reaches 150 m (67 m). On the other hand, the minimum MLD in these regions during March-April-May becomes as low as 11–12 m. The influence of wind stress, net heat flux and freshwater flux on the seasonal variability of the MLD is discussed. The physical processes controlling the seasonal cycle of sea surface temperature are investigated by carrying out mixed layer heat budget analysis. It is found that air-sea fluxes play a dominant role in the seasonal evolution of sea surface temperature of the northern Indian Ocean and the contribution of horizontal advection, vertical entrainment and diffusion processes is small. The upper ocean zonal and meridional volume transport across different sections in the AS and BoB is also computed. The seasonal variability of the transports is studied in the context of monsoonal currents.     

A comparison of two global ocean-ice coupled models with different horizontal resolutions

A global eddy-permitting ocean-ice coupled model with a horizontal resolution of 0.25 by 0.25 is established on the basis of Modular Ocean Model version 4 (MOM4) and Sea Ice Simulator (SIS). Simulation results are compared with those of an intermediate resolution ocean-ice coupled model with a horizontal resolution of about 1 by 1 . The results show that the simulated ocean temperature, ocean current and sea ice concentration from the eddy-permitting model are better than those from the intermediate resolution model. However, both the two models have the common problem of ocean general circulation models (OGCMs) that the majority of the simulated summer sea surface temperature (SST) is too warm while the majority of the simulated subsurface summer temperature is too cold. Further numerical experiments show that this problem can be alleviated by incorporating the non-breaking surface wave-induced vertical mixing into the vertical mixing scheme for both eddy-permitting and intermediate resolution models.     

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%。     

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.     

The seasonal variation of undercurrent and temperature in the equatorial Pacific jointly derived from buoy measurement and assimilation analysis

Based on the TOGA-TAO buoy chain observed data in the equatorial Pacific and the assimilation analysis results from SODA(simple ocean data assimilation analysis), the role of the meridional cells in the subsurface of the tropical Pacific was discussed. It was found that, the seasonal varying direction of EUC (the quatorial Undercurrent)in the Peacific is westwards beginning from the eastern equatorial Pacific in the boreal spring. The meridional cell south of the equator plays important role on this seasonal change of EUC.On the other hand, although the varying direction is westwards,the seasonal variation of temperature in the same region gets its minimum values in the boreal autumn beginning from the eastern equatorial Pacific.The meridional cell north of the equator is most responsible for the seasonal temperature variation in the eastern equatorial Pacific while the meridional cell south of the equator mainly controls the seasonal temperature change in the central Pacific. It is probably true that the asymmetry by the equator is an important factor influencing the seasonal cycle of EUC and temperature in the tropical Pacific.     

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.     

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.     

南海及邻近海域大气-海浪-海洋耦合精细化数值预报系统研制

本文主要介绍了南海及邻近海域大气-海浪-海洋耦合精细化数值预报系统的研制概况。预报区域为99°E~135°E，15°S~45°N，包括渤海、黄海、东海和南海及其周边海域。为了给耦合预报模式提供较准确的预报初始场，在预报开始之前，分别进行了海浪模式和海洋模式的前24小时同化后报模拟。海浪模式和海洋模式都采用了集合调整Kalman滤波同化方法，海浪模式同化了Jason-2有效波高数据；海洋模式同化了SST数据、MADT数据和ARGO剖面数据。为了改进海洋温度和盐度的模拟，我们在海洋模式的垂向混合方案中引入波致混合和内波致混合的作用。预报系统的运行主要包括两个阶段，首先海浪模式和海洋模式进行了2014年1月至2015年10月底的同化后报模拟，强迫场源自欧洲气象中心的六小时的再分析数据产品。然后耦合预报系统将同化后报模拟的结果作为初始场进行了14个月的耦合预报。预报产品包括大气产品（气温、风速风向、气压等）、海浪产品（有效波高和波向等）、海流产品（温度、盐度和海流等）。一系列观测资料的检验比较表明该大气-海浪-海洋耦合精细化数值预报系统的预报结果较为可靠，可以为南海及周边海洋资源开发和安全保障提供数据和信息产品服务。     

Validating and assessing the sensitivity of the climate model with an ocean general circulation model developed at the Institute of Atmospheric Physics,Russian Academy of Sciences

A new version of the Institute of Atmospheric Physics, Russian Academy of Sciences (IAP RAS), climate model (CM) has been developed using an ocean general circulation model instead of the statistical-dynamical ocean model applied in the previous version. The spatial resolution of the new ocean model is 3° in latitude and 5° in longitude, with 25 unevenly spaced vertical levels. In the previous version of the oceanic model, as in the atmospheric model, the horizontal resolution was 4.5° in latitude and 6° in longitude, with four vertical levels (the upper quasi-homogeneous layer, seasonal thermocline, abyssal ocean, and bottom friction layer). There is no correction for the heat and momentum fluxes between the atmosphere and ocean in the new version of the IAP RAS CM. Numerical experiments with the IAP RAS CM have been performed under current initial and boundary conditions, as well as with an increasing concentration of atmospheric carbon dioxide. The main simulated atmospheric and oceanic fields agree quite well with observational data. The new version’s equilibrium temperature sensitivity to atmospheric CO₂ doubling was found to be 2.9 K. This value lies in the mid-range of estimates (2–4.5 K) obtained from simulations with state-of-the-art models of different complexities.