海洋资料同化对气候季节-年际预测技巧及初始场的影响试验

集合卡尔曼滤波(Ensemble Kalman filter, EnKF)是一种国内外广泛使用的海洋资料同化方案, 用集合成员的状态集合表征模式的背景误差协方差, 结合观测误差协方差, 计算卡尔曼增益矩阵, 有效地将观测信息添加到模式初始场中。由于季节、年际预测很大程度上受到初始场的影响, 因此资料同化可以提高模式的预测性能。本文在NUIST-CFS1.0预测系统逐日SST nudging的初始化方案上, 利用EnKF在每个月末将全场(full field)海表温度(sea surface temperature, SST)、温盐廓线(in-situ temperature and salinity profiles, T-S profiles)以及卫星观测海平面高度异常(sea level anomalies, SLA)观测资料同化到模式初始场中, 对比分析了无海洋资料同化以及加入同化后初始场的区别、加入海洋资料同化后模式提前1~24个月预测性能的差异以及对于厄尔尼诺-南方涛动(El Niño-southern oscillation, ENSO)预测技巧的影响。结果表明, 加入海洋资料同化能有效地改进初始场, 并且呈现随深度增加初始场改进越显著的特征。加入同化后, 对全球SST、次表层海水温度的平均预测技巧均有一定的提高, 也表现出随深度增加预测技巧改进越明显的特征。但加入海洋资料同化后, 模式对ENSO的预测技巧有所下降, 可能是由于模式误差的存在, 使得同化后的预测初始场从接近观测的状态又逐渐恢复到与模式动力相匹配的状态, 加剧了赤道太平洋冷舌偏西、中东部偏暖的气候平均态漂移。     

A multi-scale high-order recursive filter approach for the sea ice concentration analysis

With the development and deployment of observation systems in the ocean, more precise passive and active microwave data are becoming available for the weather forecasting and the climate monitoring. Due to the complicated variability of the sea ice concentration(SIC) in the marginal ice zone and the scarcity of high-precision sea ice data, how to use less data to accurately reconstruct the sea ice field has become an urgent problem to be solved. A reconstruction method for gridding observations ...     

Assimilation of surface currents into a regionalmodel over Qingdao coastal waters of China

Surface currents measured by high frequency （HF） radar arrays are assimilated into a regional ocean model over Qingdao coastal waters based on Kalman filter method. A series of numerical experiments are per- formed to evaluate the performance of the data assimilation schemes. In order to optimize the analysis pro- cedure in the traditional ensemble Kalman filter （ENKF）, a different analysis scheme called quasiensemble Kaman filter （QENKF） is proposed. The comparisons between the ENKF and the QENKF suggest that both them can improve the simulated error and the spatial structure. The estimations of the background error covariance （BEC） are also assessed by comparing three different methods： Monte Carlo method; Canadian quick covariance （CQC） method and data uncertainty engine （DUE） method. A significant reduction of the root-mean-square （RMS） errors between model results and the observations shows that the CQC method is able to better reproduce the error statistics for this coastal ocean model and the corresponding external forcing. In addition, the sensibility of the data assimilation system to the ensemble size is also analyzed by means of different scales of the ensemble size used in the experiments. It is found that given the balance of the computational cost and the forecasting accuracy, the ensemble size of 50 will be an appropriate choice in the Qingdao coastal waters.     

Assimilation of temperature and salinity using isotropic and anisotropic recursive filters in Tropic Pacific

A data assimilation scheme used in the updated Ocean three-dimensional Variational Assimilation System (OVALS),OVALS2,is described.Based on a recursive filter (RF) to estimate the background error covariance (BEC) over a predetermined scale,this new analysis system can be implemented with anisotropic and isotropic BECs.Similarities and differences of these two BEC schemes are briefly discussed and their impacts on the model simulation are also investigated.An idealized experiment demonstrates the ability of the updated analysis system to construct different BECs.Furthermore,a set of three years experiments is implemented by assimilating expendable bathythermograph (XBT) and ARGO data into a Tropical Pacific circulation model.The TAO and WOA01 data are used to validate the assimilation results.The results show that the model simulations are substantially improved by OVALS2.The inter-comparison of isotropic and anisotropic BEC shows that the corresponding temperature and salinity produced by the anisotropic BEC are almost as good as those obtained by the isotropic one.Moreover,the result of anisotropic RF is slightly closer to WOA01 and TAO than that of isotropic RF in some special area (e.g.the cold tongue area in the Tropic Pacific).     

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

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

背景误差协方差矩阵不同求逆方案在高度计资料同化试验中的应用比较

针对海表高度计资料的同化,考查了背景误差协方差矩阵的不同求逆方案对同化效果的影响。所使用的求逆方案包括避免求逆的经验正交函数方案(EOF/EOF_var)、递归滤波方案(RF/RF_var)以及采用初等变换法直接求逆的方案(Inv)。基于上述方案开展了热带太平洋地区2002年1-7月的TOPEX/Poseidon高度计资料同化试验,并利用SODA再分析资料和TAO观测资料评估了各方案对温度场的同化效果,主要得到如下结果:与SODA相比较,Inv方案对模式温度场改进甚微,其余四种方案在100~300 m深度之间对温度场改进较多,在其它深度范围内则改进较少;与TAO观测相比较,EOF_var、RF_var方案对模式温度场改进最多,EOF和RF方案次之,Inv方案则对温度场改进甚少。     

基于集合同化方法的南海北部最优观测实验

根据集合同化理论及最优观测方法,应用基于ROMS模式的南海北部区域模拟结果作为集合样本,对该区域的最优观测廓线以及最优观测航次进行了评估和选择.最优观测方法是一种不需要直接同化观测,而通过最小化分析场误差方差理论值来选择最优观测点的方法.最优观测方法根据卡曼滤波理论中协方差矩阵的更新方程,使用一组模式或观测得到的历史样...     

中国近海现场海洋观测系统设计评估

中国科学院正在发展一个在中国近海(包括黄海、东海和南海)现场海洋观测系统。观测系统包括3个沿岸观测站点、4个近海离岸浮标和由观测船只按固定航线做的船舶观测断面。观测站点、浮标和断面的位置已经预先确定,这个计划在2008-2011实施。利用基于卡尔曼理论的样本集合方法对这样一个能够监测大尺度的季节和年季变率的观测系统设计进行了评估。根据卡尔曼滤波理论,用集合样本的方法能够给出经过同化这个观测系统位置的观测资料后能够减少多少分析误差和分析场的不确定性。用2个来自不同模式、不同分辨率的模式的结果作为集合样本来计算静态的背景误差协方差,这2套样本分别是来自分辨率是0.5°×0.5°的模式同化结果和高分辨0.125°×0.125°的模式结果。由这2个不同资料得到的结果是一致的。发现来自3个近岸和4个离岸浮标得到的观测能够有效地减少SST在渤海、黄海、东海和南海中部的分析误差。然而在越南东部和台湾东部海域,分析误差减少的百分比相对要小。最后,给出了中国近海最优的观测位置序列设计。     

Evaluating the contribution of satellite measurements to the re-construction of three-dimensional ocean temperature fields in combination with Argo profiles

Assimilation systems absorb both satellite measurements and Argo observations. This assimilation is essential to diagnose and evaluate the contribution from each type of data to the reconstructed analysis, allowing for better configuration of assimilation parameters. To achieve this, two comparative reconstruction schemes were designed under the optimal interpolation framework. Using a static scheme, an in situ-only field of ocean temperature was derived by correcting climatology with only Argo ...     

Coupling a two-way nested primitive equation model and a statistical SST predictor of the Ligurian Sea via data assimilation

A primitive equation model and a statistical predictor are coupled by data assimilation in order to combine the strength of both approaches. In this work, the system of two-way nested models centred in the Ligurian Sea and the satellite-based ocean forecasting (SOFT) system predicting the sea surface temperature (SST) are used. The data assimilation scheme is a simplified reduced order Kalman filter based on a constant error space. The assimilation of predicted SST improves the forecast of the hydrodynamic model compared to the forecast obtained by assimilating past SST observations used by the statistical predictor. This study shows that the SST of the SOFT predictor can be used to correct atmospheric heat fluxes. Traditionally this is done by relaxing the model SST towards the climatological SST. Therefore, the assimilation of SOFT SST and climatological SST are also compared.     

多源资料重构三维海温场的先验误差估计和诊断分析

针对遥感资料与实况观测融合重构三维海温场的问题，改进了最优插值中的先验统计信息估计，并研究了背景场和观测项的各自影响机制。在最优插值中，将背景场分别取为平均气候场（静态方案）和遥感合成场（动态方案）；在最优融合前，用后验诊断方法迭代优化了背景误差和观测误差协方差；融合完成后，用观测空间诊断误差和模式格点指标对比分析两种方案。主要结论:（1）动态方案相对于静态方案分析误差减小的绝对量值由二者的观测误差之差决定、减小的相对程度由二者的背景误差之比决定；（2）背景场的调整作用在中高纬占优势、在赤道海域与观测项调整作用相当，是由协方差相关尺度的准纬向分布特征决定的；（3）动态分析场的中尺度信号谱能量整体比静态分析场高1~3个量级，但在热带海域相当；（4）遥感资料通过分辨实况观测所无法分辨出来的海温中尺度特征，达到减小误差、提高有效分辨率的效果。     

黄海海雾短时临近预报中云水路径的EnKF同化研究

在海雾的短时临近预报中,初始场的水汽凝结状态扮演着重要角色。为了改进初始场的云水含量,本文提出直接同化雾体云水信息的思路。针对2011年5月一次大范围的黄海海雾,借助EnKF (Ensemble Kalman Filter)方法,尝试进行了极轨卫星反演云水路径数据的同化试验。结果表明:(1)通过利用EnKF将云水混合比增加到背景场和分析场的控制变量中,构建云水观测数据与背景场之间的关系,实现云水路径数据的直接同化是可行的;(2)同化云水路径可显著改善海面气温与湿度状态,大幅提高海雾预报效果;(3)EnKF能够基于集合体动态统计流依赖的背景误差协方差是其取得良好同化效果的主要原因。值得指出的是,受集合样本误差的影响,需要特别关注云水含量与风之间的相关关系。     

各向异性背景场误差协方差构建及在“凡亚比”台风的应用

利用相临过去时段预报结果中同一时刻不同时效的模式预报场差异,计算预报误差协方差,并基于集合-变分混合同化系统将其与静态背景场误差协方差结合,从而在同化系统中构建了具有各向异性和一定流依赖特征的背景场误差协方差。单点观测理想试验显示本方案改善了静态模型化背景场误差协方差的各向同性和流依赖性问题。“凡亚比”台风的一系列同化及模拟试验表明,从台风路径、强度等方面本文方案的效果都要优于三维变分法。本文方案在不需要集合预报,计算量与三维变分法相当的情况下,给同化系统引入了各向异性、一定流依赖特征的背景误差协方差,因此本方案适于在计算资源较为紧缺情况下,对时效要求较高的预报业务中应用。     

斜压罗斯贝变形半径优化的误差相关尺度及其对最优插值效果的改进

利用最优插值方法进行三维温盐场重构时,如何正确估算背景场误差协方差至关重要。针对背景场误差协方差主要取决于误差相关尺度的问题,本文提出了用中国东部海域及其附近海域的最新的高分辨率气候态数据进行斜压罗斯贝变形半径优化其误差相关尺度计算的研究方案和技术途径;对比分析了均一化相关尺度方案和法国ISAS系统尺度方案,讨论了变形半径对最优插值的影响。结果表明:均一化相关尺度方案的均方根误差小于ISAS方案,但温度场过于平滑,难以刻画一些重要的物理现象;相比而言,本文提出的基于变形半径的相关尺度方案在取2倍变形半径时,不仅均方根误差在各水平层较小,且温度场能够更好地刻画四国海盆海域涡旋及黑潮影响的温度场三维结构。由于实际海洋中各层物理过程的尺度存在差异,实际应用时各层的最优尺度设置也应有所不同。     

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.     

利用海洋温度剖面与海表盐度反演盐度剖面方法研究

为解决海洋中大量观测数据只含有温度剖面而缺乏盐度观测的问题, 基于历史观测的温盐剖面资料, 考虑到盐度卫星数据的发展, 采用回归分析方法, 在孟加拉湾建立了盐度与温度、经纬度、表层盐度的关系, 并对不同反演方法的反演结果进行检验评估。结果发现, 在不引入海表盐度(sea surface salinity, SSS)时, 最佳反演模型是温度、温度的二次项与经纬度确定的回归模型, 而SSS的引入则可以进一步优化反演结果。将反演结果与观测结果进行对比, 显示用反演的盐度剖面计算的比容海面高度误差超过2cm, 而引入SSS后的误差低于1.5cm。SSS的引入能够较为真实地反映海洋盐度场的垂直结构和内部变化特征, 既能够捕捉到对上混合层有重要影响的SSS信号, 又能够反映盐度在跃层上的季节内变化以及盐度障碍层的季节变化。水团分析显示, 与气候态相比, 盐度反演结果可以更好地表征海洋上层水团的变化特征。     

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.     

黄、东海海洋对于台风过程的响应

强天气过程下海洋的响应是海洋环境预报的重要内容,它的研究对于防灾减灾、远洋运输、水产养殖等有积极的意义。本文中作者运用一个改进后具有模拟和预报能力的河口海洋模式ECOM-si,引入全强迫条件,对一个真实的强台风下的海洋响应进行了研究。模拟结果表明,海洋对台风过程有强烈的响应,强台风引起SST出现大幅降低,最大达5℃,其中大风抽吸和大风夹卷影响最大;大风引起的平流输送在实际情况中对海水温度、盐度的水平分布有重要影响:台风诱导海水上翻,会使得海洋的混合层明显加深,最大达20—30m;海洋在台风作用下,在上层海洋产生明显的气旋式流场,海面产生的明显下陷可达30cm,台风中心、气旋式流场中心和海面下陷中心三个位并不重合。同时在台风登陆位置附近产生风暴潮,最大增水可达0.8—1.0m,但沿海各地最大增水时间不相同。     

Formation and transport of intermediate water masses in a model of the Pacific Ocean

A basin-wide ocean general circulation model of the Pacific Ocean was used to investigate how the interior restoration in the Okhotsk Sea and the isopycnal diffusion affect the circulation and intermediate water masses. Four numerical experiments were conducted, including a run with the same isopycnal and thickness diffusivity of 1.0×103 m2/s, a run employing the interior restoration of temperature and salinity in the Okhotsk Sea with a time scale of 3 months, a run that is the same as the first run except for the enhanced isopycnal mixing, and a final run with the combination of the restoration in the Okhotsk Sea and large isopycnal diffusivity. Simulated results show that the intermediate water masses reproduced in the first run are relatively weak. An increase in isopycnal diffusivity can improve the simulation of both Antarctic and North Pacific intermediate waters, mainly increasing the transport in the interior ocean, but inhibiting the outflow from the Okhotsk Sea. The interior restoration generates the reverse current from the observation in the Okhotsk Sea, whereas the simulation of the temperature and salinity is improved in the high latitude region of the Northern Hemisphere because of the reasonable source of the North Pacific Intermediate Water. A comparison of vertical profiles of temperature and salinity along 50°N between the simulation and observations demonstrates that the vertical mixing in the source region of intermediate water masses is very important.     

An effectivemethod for improving the accuracy of Argo objective analysis

Based on the optimal interpolation objective analysis of the Argo data, improvements are made to the em- pirical formula of a background error covariance matrix widely used in data assimilation and objective anal- ysis systems. Specifically, an estimation of correlation scales that can improve effectively the accuracy of Ar- go objective analysis has been developed. This method can automatically adapt to the gradient change of a variable and is referred to as ＂gradient-dependent correlation scale method＂. Its effect on the Argo objective analysis is verified theoretically with Gaussian pulse and spectrum analysis. The results of one-dimensional simulation experiment show that the gradient-dependent correlation scales can improve the adaptability of the objective analysis system, making it possible for the analysis scheme to fully absorb the shortwave information of observation in areas with larger oceanographic gradients. The new scheme is applied to the Argo data obiective analysis system in the Pacific Ocean. The results are obviously improved.