Adjoint Assimilation in Marine and an Example of Application

This paper aims at a review of the work carried out to date on the adjoint assimilation of data in marine ecosystem models since 1995. The structure and feature of the adjoint assimilation in marine ecosystem models are also introduced. To illustrate the application of the adjoint technique and its merits, a 4-variable ecosystem model coupled with a 3-D physical model is established for the Bohai Sea and the Yellow Sea. The chlorophyll concentration data derived from the SeaWiFS ocean colour data are assimilated in the model with the technique. Some results are briefly presented.     

IMPACT OF DIABATIC PROCESSES IN AGCM ON 4-DIMENSIONAL VARIATIONAL DATA ASSIMILATION*

The impact of diabatic processes on 4-dimensional variational data assimilation (4D-Var) was studied using the 1995 version of NCEP's global spectral model with and without full physics.The adjoint was coded manually.A cost function measuring spectral errors of 6-hour forecasts to "observation" (the NCEP reanalysis data) was minimized using the L-BFGS (the limited memory quasi-Newton algorithm developed by Broyden,Fletcher,Goldfard and Shanno) for optimizing parameters and initial conditions.Minimization of the cost function constrained by an adiabatic version of the NCEP global model converged to a minimum with a significant amount of decrease in the value of the cost function.Minimization of the cost function using the diabatic model, however,failed after a few iterations due to discontinuities introduced by physical parameterizations.Examination of the convergence of the cost function in different spectral domains reveals that the large-scale flow is adjusted during the first 10 iterations,in which discontinuous diabatic parameterizations play very little role.The adjustment produced by the minimization gradually moves to relatively smaller scales between 10-20th iterations.During this transition period,discontinuities in the cost function produced by "on-off" switches in the physical parameterizations caused the cost function to stay in a shallow local minimum instead of continuously decreasing toward a deeper minimum. Next,a mixed 4D-Var scheme is tested in which large-scale flows are first adiabatically adjusted to a sufficient level,followed by a diabatic adjustment introduced after 10 to 20 iterations. The mixed 4D-Var produced a closer fit of analysis to observations,with 38% and 41% more decrease in the values of the cost function and the norm of gradient,respectively,than the standard diabatic 4D-Var,while the CPU time is reduced by 21%.The resulting optimal initial conditions improve the short-range forecast skills of 48-hour statistics.The detrimental effect of parameterization discontinuities on minimization was also reduced.     

On ensemble representation of the observation-error covariance in the Ensemble Kalman Filter

Evensen (2003) presents a modification of the Ensemble Kalman Filter (EnKF), in which the observation-error and background-error covariance matrices are both represented by ensembles, in contrast to the usual practice, where only the background error is so represented. It is shown that this modification can cause the ensemble to collapse to a single member, in the common situation where the number of observations is more than twice the number of ensemble members, and to be rank-deficient when the number of observations is greater than or equal to the ensemble size. It is also shown that some further modifications to the scheme, presented by Evensen as offering numerical efficiencies, can prevent this collapse. However, these latter modifications are shown in some simple numerical examples to require tuning to produce acceptable results, which are nevertheless inferior to those of the standard EnKF.Acknowledgements The author acknowledges useful discussions with Peter Steinle, and other participants at the EnKF workshop held in BMRC in November, 2003.     

LAPS分析场在一次强对流天气过程尺度分析中的应用

利用LAPS（local analysis and prediction system）同化系统融合多种观测资料,对2010年8月25日发生在上海的一次强对流天气过程进行中尺度分析。结果表明,这次强对流天气主要是由中、低空中尺度辐合系统直接触发形成的。强对流形成阶段,地面有分散辐合形成并逐渐加强,成为触发中尺度垂直环流的主要机制,垂直结构上出现低层辐合高层辐散的有利配置,风暴中心附近出现明显的上升气流区,中高层相对湿度显著增加。成熟阶段,强对流云体中心附近的对流层底层开始出现下沉气流,上升气流在其拖曳作用下明显倾斜。衰减阶段,下沉气流加强使中尺度环流动力结构和水汽供应受到破坏,垂直结构上转为底层辐散高层辐合。因此,与天气尺度分析相比,基于LAPS的中尺度分析能更深刻地揭示中小尺度系统的三维结构和时空演变特征。     

散射计和下投式探空仪资料对台风模拟的影响

为了评估同化QuikSCAT海平面气压场资料和Dropsonde资料对台风"Bilis"和"Kat-rina"模拟的影响,通过WRF-3DVAR变分同化系统将QuikSCAT海平面气压场资料和Dropsonde资料同化到WRF模式中,并以未同化任何资料和同化了QuikSCAT海平面气压场以及同化了Drop-sonde资料的数据为初始场,应用WRF模式进行模拟预报,并对得到的初始分析场和模拟预报得到的结果与台风实况报文进行了对比,得到的结论是:对台风初始时刻海平面气压场分布、初始台风中心定位、台风初始暖心结构的改善、台风路径、海平面最低气压场、海面最大风速和雷达回波的模拟方面,同化QuikSCAT海平面气压场的效果比同化Dropsonde资料明显,其中二者对海平面最低气压场的改进效果比海面最大风速的要好。但是两种资料的同化对初始相对湿度和台风强度细节变化方面都没有产生明显改善作用。同化QuikSCAT海平面气压场试验虽然仅仅同化了单层海平面气压场,但是在三维变分同化系统的动力约束作用下,也可以影响到中高层的各个要素场。     

A comparison of shelf observation platforms for assimilation in an eddy-resolving ocean model

An assessment of the likely benefits of assimilating in situ temperature (T) and salinity (S) observations from repeat glider transects and surface velocity observations from high-frequency radar arrays into an eddy-resolving ocean model is presented. The deployment of new shelf observation platforms around Australia is being undertaken through the Australian Integrated Marine Observing System program. In this study, various options for an observing system along the coast of New South Wales, Australia, are assessed for their benefits to an ocean forecast and reanalysis system. The forecast system considered here uses ensemble optimal interpolation (EnOI) for data assimilation. Using error estimates from the EnOI scheme, estimates of the theoretical analysis errors are calculated for different observing systems that include a range of remotely sensed and in situ observations. The results demonstrate that if HF radar observations are assimilated along with the standard components of the global ocean observing system, the analysis errors are likely to reduce by as much as 80% for velocity and 60% for T, S and sea-level in the vicinity of the observations. Owing to the relatively short along-shore decorrelation length-scales for T and S near the shelf, the glider observations are likely to provide the forecast system with a more modest gain.     

国家气象中心业务数值预报发展的回顾与展望

本文回顾了国家气象中心建立数值预报业务以来,国家级业务数值预报系统,主要是中、短期预报系统的发展历程,评估了目前业务数值预报的预报水平和基本特点,简述了我国业务数值预报发展的不同阶段取得的主要进展。通过对目前我国业务数值预报发展中存在的问题和不足的分析,讨论下一步的工作重点,并对未来的发展做初步展望。     

有云环境下MODIS亮温资料的变分同化Ⅱ——对暴雨预报的影响

卫星资料提供了大量关于云和雨的观测信息,在暴雨预报中可发挥巨大的作用,然而在数值模式资料同化中的应用水平仍然不高,特别是红外辐射资料的应用。由于有云环境下辐射传输过程的模拟难度很大,因此通常只同化晴空环境下的红外辐射资料。基于GRAPES-3DVAR(Global and Regional Assimilation and Prediction Enhanced System,全球/区域同化预报系统),根据RTTOV辐射传输模式(fast radiative transfer model for TOVS,快速辐射传输模式)的特点,增加云水含量、云冰水含量和云量作为同化系统控制变量,在改进辐射传输模式对红外资料模拟的同时,利用红外资料调整初始云参数和大气参数。针对2007年5月26日南海季风爆发后广东地区的一次暴雨过程,选取MODIS(Moderate Resolution Imaging Spectroradiometer,中分辨成像光谱仪)传感器水汽(第27)和云顶观测(第36)通道进行了同化试验,利用WRF(Weather Research and Forecasting Model,天气研究和预报模式)进行了数值模拟,结果表明同化MODIS资料,可以改进初始场水汽和温度分布,间接调整高空风场,调整趋势符合卫星观测,对短时降水预报有正面影响。     

Three-dimensional variational data assimilation of WindSat ocean surface winds for the genesis and forecasting of tropical storm Henri

With available high-resolution ocean surface wind vectors retrieved from the U.S. Naval Research Laboratorys WindSat on Coriolis, the impact of these data on genesis and forecasting of tropical storm Henri is examined using the non-hydrostatic, fifth-generation mesoscale model (MM5) of Pennsylvania State University-National Center for Atmospheric Research plus its newly released three-dimensional variational data assimilation (3DVAR) system. It is shown that the assimilation of the WindSat-retrieved ocean surface wind vectors in the 3DVAR system improves the model initialization fields by introducing a stronger vortex in the lower troposphere. As a result, the model reproduces the storm formation and track reasonably close to the observations. Compared to the experiment without the WindSat surface winds, the WindSat assimilation reduced an error between the model simulated track and observations of more than 80 km and also improved the storm intensity by nearly 2 hPa. It suggests that these data could provide early detection and prediction of tropical storms or hurricanes.