卫星遥感海表温度资料和高度计资料的变分同化

利用国家气候中心正在发展的第二代全球海洋资料同化系统(BCC_GODAS2.0),针对多变量同化的协调性问题,发展了一种基于三维变分框架(3DVAR)下的高度计和海表温度(SST)相互约束的同化方法。该方法使海面高度和SST资料在同一个动力约束关系下进行同化。在一般方法中,海面高度和SST观测项是代价函数中2个独立的观测项,海面高度项引入动力高度计算公式,海表温度项用统计关系进行垂向投影。在代价函数的实际求解的计算过程中,虽然其总体积分效应受海面高度观测的约束,但整个水柱中各层温盐分析变量的调整是无序的。针对这个问题,文章提出一种新的同化方案。该方案将SST的观测项并入海面高度观测项中,海面高度的一部分,确切说是上层海洋部分,由SST决定,因此至少在SST的统计关系能影响到深度的上层海洋,在代价函数的求解过程中,温盐的调整是受较强的统计关系约束的,而这种统计关系的有效性已经在很多SST的同化试验中被其他学者广泛应用并证明。利用该方法,对1993—1997年的AVHRR卫星遥感海表温度资料进行变分同化试验,用TAO、OISST和SODA数据集进行检验证明,通过对卫星遥感资料的同化能够有效改进对海洋温度和盐度的估...     

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

集合卡尔曼滤波(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的预测技巧有所下降, 可能是由于模式误差的存在, 使得同化后的预测初始场从接近观测的状态又逐渐恢复到与模式动力相匹配的状态, 加剧了赤道太平洋冷舌偏西、中东部偏暖的气候平均态漂移。     

黄海、渤海TOPEX/Poseidon高度计资料潮汐伴随同化

首先将大约10a的TOPEX/Poseidon(T/P)高度计资料沿星下轨迹点做潮汐调和分析,提取得到各分潮的调和常数,利用伴随同化方法,同化到二维非线性潮汐数值模式中,模拟了黄海、渤海区域M₂,S₂,O₁,K₁等4个潮汐分潮,并根据计算结果给出了各分潮的同潮图.将计算值与观测值的进行偏差统计,结果表明计算值与验潮站资料符合良好.研究过程中做了两类试验:一类试验是针对不同的参数进行优化,一类试验是针对不同的资料进行同化.第一类试验表明:将开边界条件和底摩擦系数同时作为模型优化的控制参数,其结果明显优于单独优化开边界条件;第二类试验表明:同时同化高度计资料与验潮站资料,比单独同化其中任一种资料,对模式计算结果都有较好的改进.研究结果表明,采用伴随同化方法,利用T/P高度计资料和验潮站资料作为同化数据能有效改进模拟结果,用来反演黄海、渤海的潮波系统是可行的.     

卫星高度计波高数据同化对西北太平洋海浪数值预报的影响评估

通过数值试验验证卫星高度计波高数据同化对西北太平洋3 d海浪预报的改进效果。驱动海浪模式的强迫场采用国家海洋环境预报中心基于MM5模式预报的风场,波高数据同化使用的观测数据是Jason-1卫星高度计有效波高。用最优插值数据同化方法获得海浪有效波高的最优估计并重构相应的海浪方向谱,以此为初始场进行为期3 d的数值预报试验。与没有同化的预报进行了比较和分析,结果表明卫星高度计海浪数据同化对0~72 h预报有不同程度的明显改善,改进程度随预报时效的增加而减少。     

北太平洋副热带潜沉率及其变化中海面风的作用

采用中国科学院大气物理研究所的高分辨率逐日风场驱动的全球海洋模式(LICOM1·0)对北太平洋海域的潜沉(Subduction)过程和副热带模态水形成区潜沉率的年际变化进行了数值模拟,并将模拟结果与同化的海洋模式资料(SODA)进行了比较。研究结果表明,该高分辨率的海洋模式对北太平洋的绝大部分海域晚冬混合层底水质点的运动方向和路径的数值模拟结果较好,模式模拟的副热带环流比SODA资料中的副热带环流流速强;模式模拟的混合层深度比SODA资料中的混合层深度深,更接近观测;模式中副热带海域的潜沉率大于SODA资料中的副热带海域的潜沉率。模式结果表明,副热带环流和副极地环流交界处是潜沉过程发生的最主要区,该区气候平均的潜沉率超过100m/a,最大为150m/a,海面风变异引起的海洋平流的年际和年代际变化,是该区潜沉率发生年际和年代际变化的主要原因;在太平洋副热带东部模态水形成区,气候平均的潜沉率超过50m/a,在该区潜沉率的年际变化中,局地风应力旋度决定的Ekman抽吸要比海洋平流效应更加重要。     

卫星高度计波高资料的同化试验分析

结合经验的台风风场模型和NCEP再分析资料构造了南海台风风场,并以此作为一个第三代海浪模式的输入模拟了1999年约克(York)台风经过南海时所产生的台风浪场.模拟结果显示,尽管文中构造的南海台风风场比较接近真实的风场,但模拟出的台风浪场相对于TOPEX/Poseidon高度计的观测波高仍有一定差异,因此进而采用了一种简单的最优插值同化方法,开展了波高资料的初步同化试验.试验结果显示,资料同化显著地改进了海浪模式的(后)预报精度,而这种改进在涌浪区更加显著,同化有效维持时间大约是43h.     

最优化插值同化方法在预报南海台风浪中的应用

由经验的holland台风模型和NCEP再分析风场资料相结合构造出南海台风风场，结果较好地符合了TOPEX／Poseidon(T／P)卫星高度计观测的风速分布。以此作为第三代海浪模式的输入风场，模拟了1999年约克(York)台风经过南海海域的台风浪，并利用T／P卫星高度计观测的有效波高资料对模式进行同化。结果显示，同化影响半径取为2000km效果较好，同化影响时间是35h，同化改善了模式预报的精度。     

COSMIC掩星资料同化对台风“天兔”预报影响的试验

台风数值预报的关键之一在于给定准确的模式初始场。本文使用气象、电离层和气候星座观测系统(COSMIC)的卫星资料,通过资料同化技术改善数值天气预报模式的初始场,进而评估数据同化技术对台风(或热带气旋)预报的改善效果。使用我国自主研发的数值天气预报模式,对1319号台风“天兔”进行模拟试验,对温度场、湿度场、涡度散度场及假相当位温场进行诊断分析。试验结果表明,与同化常规探空资料相比,COSMIC资料的同化可减小模式预报台风路径误差,台风附近降水落区更为集中。研究评估了掩星资料同化对模式的改进作用,并能为台风预报技术和防灾减灾提供参考。     

基于适应性观测误差和温盐偏差控制的EnKF温盐廓线资料同化试验

通过开展2008年夏季南海北部开放航次CTD的温盐廓线数据资料同化试验,本文采取了观测误差适应的方法来防止EnKF滤波发散问题;同时,从背景误差协方差和温盐模式偏差关系入手,在同化中引入温盐控制来减小模式偏差对同化结果的影响。对于改进的同化方案进行了试验验证,并用卫星高度计观测数据,OSCAR流速数据,走航ADCP数据作为独立观测数据检验。结果证明新的EnKF同化策略能够有效地减小温盐均方根误差。同时整个同化系统能有效地改善高度场和流场的模拟。     

MASNUM-WAM海浪模式集合Kalman滤波同化研究——Ⅱ.集合样本对同化效果的影响

背景误差相关结构的确定是影响海浪同化效果的关键因素之一。集合Kalman滤波是一种较为成熟的同化方法,其可以对背景误差进行实时更新和动态估计,现已广泛应用于海洋和大气领域的研究。本文基于MASNUM-WAM海浪模式,分别采用静态样本集合Kalman滤波和EAKF方法,针对2014年全球海域开展海浪数据同化实验,同化资料为Jason-2卫星高度计数据,利用Saral卫星高度计资料对同化实验结果进行检验。结果表明,两组同化方案均有效提高了海浪模式的模拟水平,EAKF方案在风场变化较大的西风带区域表现显著优于静态样本集合Kalman滤波方案,但总体上两者相差不大。综合考虑计算成本和同化效果,静态样本集合Kalman滤波方案更适用于海浪业务化预报。     

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).     

Sequential forecasting of the surface and subsurface conditions in the Japan Sea

This study estimates a realistic change of the Japan Sea by assimilating satellite measurements into an eddy-resolving circulation model. Suboptimal but feasible assimilation schemes of approximate filtering and nudging play essential roles in the system. The sequential update of error covariance significantly outperforms the asymptotic covariance in the sequential assimilation due to the irregular sampling patterns from multiple altimeter satellites. The best estimates show an average rms difference of only 1.2°C from the radiometer data, and also explain about half of the sea level variance measured by the altimeter observation. The subsurface conditions associated with the mesoscale variabilities are also improved, especially in the Tsushima Warm Current region. It is demonstrated that the forecast limit strongly depends on variable, depth, and location.     

中国近海以及邻近海域海洋再分析产品(CORA)中季节内变化的评估

A regional reanalysis product—China Ocean Reanalysis(CORA)—has been developed for the China's seas and the adjacent areas. In this study, the intraseasonal variabilities(ISVs) in CORA are assessed by comparing with observations and two other reanalysis products(ECCO2 and SODA). CORA shows a better performance in capturing the intraseasonal sea surface temperatures(SSTs) and the intraseasonal sea surface heights(SSHs) than ECCO2 and SODA do, probably due to its high resolution, stronger response to the intraseasonal forcing in the atmosphere(especially the Madden-Julian Oscillation), and more available regional data for assimilation. But at the subsurface, the ISVs in CORA are likely to be weaker than reality, which is probably attributed to rare observational data for assimilation and weak diapycnal eddy diffusivity in the CORA model. According to the comparison results, CORA is a good choice for the study related to variabilities at the surface, but cares have to be taken for the study focusing on the subsurface processes.     

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.     

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

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

Analysis of interdecadal variation of tropical Pacific thermocline based on assimilated data

1IntroductionAfterNittaandYamada(1989)andTren-berth(1990)haveinitiallydescribedthe1977climaticvariationovertheNorthPacificwithstep-likeshiftofwintermeansealevelpressurefield,thePacificinterdecadaloscillation(PDO)hasextensivelybeenstudied.TheworkaboutthecauseofPDOfromtheviewoftheNorthPacificincludesthermodynamicexplanationoflow-frequencysea-surfacetemperature(SST)affectedbyatmosphericrandomforcing(FrankignoulandHasselmann,1977),themeridionalchangeofthermaltransportbysubtropicalgyreandhen…     

A 4D-variational approach applied to an eddy-permitting North Atlantic configuration: Synthetic and real data assimilation of altimeter observations

The increasing number of oceanic observations calls for the use of synthetic methods to provide consistent analyses of the oceanic variability that will support a better understanding of the underlying mechanisms. In this study, a 1/3° eddy-permitting model of the North Atlantic (from 20°S to 70°N) is combined with a 4D-variational method to estimate the oceanic state from altimeter observations. This resolution allows a better extraction of the physical content of altimeter data since the model spatial scales are more consistent with the data than coarser assimilation exercises because of a lower error in model representativity. Several strategies for the assimilation window are tested through twin experiments carried out under the following conditions: different window lengths and either a quasi-static (also known as progressive) variational assimilation with progressive extension of the window, or a simpler direct method without prior assimilation. From our set of experiments, the most efficient strategy is the use of both a simple direct assimilation method and a 90-day window. The assimilation of synthetic altimeter data constrains the model-temperature, -salinity and -velocity fields mainly over the first 1300 m where the error is the largest. Improvements occur not only in quiescent regions, but also in more energetic meso-scale regimes. Despite the existence of model- and surface forcing-errors as well as large errors in the first guess, the assimilation of real altimeter data proves to be consistent with our twin experiments. Indeed, the analyses show a better detachment of the Gulf Stream, weaker regional biases and more accurate positions for meso-scale structures. Independent hydrographic data (Argo floats and CTD cruises) and transports estimates along the OVIDE 2002 cruise show an improvement of the analysed oceanic state with respect to the assimilation-free case though water mass properties are still incorrectly represented. After assimilation, the North Atlantic heat transport in the model is in good agreement with independent estimates based on hydrographic data.     

Data assimilation of partitioned HF radar wave data into Wavewatch III

In this study the assimilation of HF radar data into a high resolution, coastal Wavewatch III model is investigated. An optimal interpolation scheme is used to assimilate the data and the design of a background error covariance matrix which reflects the local conditions and difficulties associated with a coastal domain is discussed. Two assimilation schemes are trialled; a scheme which assimilates mean parameters from the HF radar data and a scheme which assimilates partitioned spectral HF radar data. This study demonstrates the feasibility of assimilating partitioned wave data into a coastal domain. The results show that the assimilation schemes provide satisfactory improvements to significant wave heights but more mixed results for mean periods. The best improvements are seen during a stormy period with turning winds. During this period the model is deficient at capturing the change in wave directions and the peak in the waveheights, while the high sea state ensures good quality HF radar data for assimilation. The study also suggests that there are both physical and practical advantages to assimilating partitioned wave data compared to assimilating mean parameters for the whole spectrum.     

A comparison between 3DVAR and EnOI techniques for satellite altimetry data assimilation

Two conceptually different assimilation schemes, three dimensional variational (3DVAR) assimilation and Ensemble Optimum Interpolation (EnOI) are compared in the context of satellite altimetric data assimilation. Similarities and differences of the two schemes are briefly discussed and their impacts on the model simulation are investigated.With a tropical Pacific ocean model, two assimilation experiments of sea level anomaly (SLA) data from TOPEX/Poseidon are performed for 5 years from 1997 to 2001. Annual mean states of temperature and salinity fields are compared with analysis data and some independent observations. It is found that EnOI generally produces moderate improvements on both temperature and salinity fields, while changes induced by 3DVAR assimilation are strong and vary remarkably in different areas. For instance, 3DVAR tends to excessively modify the temperature field along the thermocline depth and even deteriorate the simulation, but it is more effective than EnOI below the thermocline depth. However, for the salinity field 3DVAR outperforms EnOI nearly for almost the whole layer. As the difference relative to the WOA01 analysis is compared, it is apparently reduced to below 0.3 psu in most areas in the 3DVAR experiment. On the other hand, the pattern of difference in the EnOI experiment resembles that of the simulation and the magnitude is only diminished to some extent. One advantage of EnOI is that it yields more consistent improvements even in areas where there are large model errors. It is more reliable than 3DVAR in such a sense. It is also revealed that the T–S relation plays a very important role in altimetric data assimilation. Further, the distinct performance of the two schemes can be partly accounted for by their inherent assumptions and settings.