TOPEX/Poseidon高度计资料在全球海洋环流模式中的同化及分析

利用一个新的四维变分海洋资料同化系统LICOM-3DVM对TOPEX／Poseidon高度计资料进行了同化。该同化系统是在LASG／IAP气候海洋模式LICOM1．0的基础上建立起来的,所用的同化方法为三维变分映射资料同化方法3DVM。高度计观测资料是采取间接的方式进行同化,即先建立起二维海面高度距平场与三维温度场的统计关系,并由此通过观测的海面高度距平信息反演出“观测”的三维温度场,然后利用LICOM-3DVM四维变分同化系统将此反演的温度场同化到海洋模式中。作者设计了两组试验并对结果进行了比较分析,积分时间从1993年1月至2001年12月共9年时间。结果表明,由于上混合层相关系数较小,因此同化后海温没有改进;而在温跃层以及更深层次,同化后的海温均有很大程度的改善。从对赤道太平洋地区海温的气候态、季节变化和年际变化以及Nino3区的Nino指数的模拟情况来看,由于同化时将海面高度异常和海温异常之间的相关参数取为常数,没有考虑其季节和年际变化,因此,同化后对于赤道太平洋的年际变化没有改善。对于黑潮地区,由于模式的分辨率较低,同化之前没有很好地模拟出温度锋面,温度和盐度梯度都偏小,流速也偏弱;而同化后使得温度锋面和盐度梯度与WOA01更加吻合,流速增强。     

全球海洋资料同化系统ZFL_GODAS的研制和初步评估试验

本研究发展了一个全球海洋资料同化系统ZFL_GODAS。该系统是一个短期气候数值预测业务系统的子系统,为短期气候预测海气耦合模式提供全球海洋初始场。系统能够同化的观测资料包括卫星高度计资料、卫星海表温度(SST)资料,以及Argo、XBT、TAO等各种不同来源的现场温盐廓线资料。系统使用的海洋模式为中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室开发的气候系统海洋模式LICOM1.0,同化方案为集合最优插值(EnOI)方案。系统使用一个由海洋模式自由积分得到的静态样本来估计背景场误差协方差。这样的基于集合样本的背景场误差协方差具有多变量协变、各向异性的特征,且能反映海洋物理过程固有的空间尺度特征。针对EnOI同化程序的特点,开发了一套特色鲜明、负载均衡、高效的并行化同化程序。本文通过与不同类型观测资料的比较,对同化系统的性能进行了评估。通过比较海表温度和海面高度的年际变率,海表温度异常随时间的变化,SST、海面高度异常(SLA)以及次表层温盐预报产品的均方根误差,5年平均温度偏差廓线、平均盐度廓线、平均纬向流速廓线等发现:系统工作正常、同化效果较好;经过同化以后,各变量都更加接近观测,误差更小,与观测场的相关性更好,可以为短期气候预测系统提供较好的海洋初始场,也可以为物理海洋学的研究提供有效的再分析资料。     

The Impact of “Bad” Argo Profiles on Ocean Data Assimilation

Recent studies have found cold biases in a fraction of Argo profiles (hereinafter referred to as bad Array for Real-time Geostrophic Oceanography (Argo) profiles) due to the pressure drifts during 2003 and 2006. These bad Argo profiles have had an important impact on in situ observation-based global ocean heat content estimates. This study investigated the impact of bad Argo profiles on ocean data assimilation results that were based on observations from diverse ocean observation systems, such as in situ profiles (e.g., Argo, expendable bathythermograph (XBT), and Tropical Atmosphere Ocean (TAO), remote-sensing sea surface temperature products and satellite altimetry between 2004 and 2006. Results from this work show that the upper ocean heat content analysis is vulnerable to bad Argo profiles and demonstrate a cooling trend in the studied period despite the multiple independent data types that were assimilated. When the bad Argo profiles were excluded from the assimilation, the decreased heat content disappeared and a warming occurred. Combination of satellite altimetry and mass variation data from gravity satellite demonstrated an increase, which agrees well with the increased heat content. Additionally, when an additional Argo profile quality control procedure was utilized that simply removed the profiles that presented static unstable water columns, the results were very similar to those obtained when the bad Argo profiles were excluded from the assimilation. This indicates that an ocean data assimilation that uses multiple data sources with improved quality control could be less vulnerable to a major observation system failure, such as a bad Argo event.     

Data assimilation into a numerical equatorial ocean model. II. Assimilation experiments

A sequence of numerical experiments is conducted using a linear, semi-spectral equatorial ocean model and an advanced data assimilation scheme. The numerical model is based on decomposition of the oceanic fields into Kelvin and Rossby waves belonging to the baroclinic modes of a stratified equatorial ocean. The assimilation procedure finds that solution to the model equations that best fits, in the generalized least-squares sense, all observations made within some specified space-time interval. All experiments are of the ‘identical twin’ type; synthetic data are generated by sampling the observable fields produced by a control run of the model, then the data are assimilated using the same model. The sequence of numerical experiments serves two purposes; to demonstrate the performance of the assimilation procedure in the context of a fully three-dimensional, time-varying equatorial ocean model; and to examine the utility of specified data sets, in particular, observations of sea level, in estimating the state of the equatorial ocean. The results indicate that the assimilation procedure works very well when sufficient data are provided. However, sea-level data alone are not sufficient and must be supplemented with subsurface observations if more than a few baroclinic modes are allowed in the model ocean. The required amount of supplementary subsurface data (in the form of density profiles in these experiments) can be reduced by imposing smoothness contraints on the recovered model solution.     

Reducing Drift and Bias of a Global Ocean Model by Frequency-Dependent Nudging

Frequency-dependent nudging is applied to a coarse resolution (nominal 1°) global ocean model to suppress its drift and bias, and the impact of the nudging on the skill of the model is assessed. The nudging is applied to temperature and salinity in frequency bands centred on 0 and 1 cycles per year. As expected, the nudging significantly reduces the biases in the long-term mean and annual cycle of temperature, salinity, and sea level. By comparing the simulated (i) sea surface temperature with operational analyses based on observations, (ii) vertical profiles of temperature and salinity with observations made by Argo floats, and (iii) sea level with altimeter observations, it is shown that the skill of the model in simulating variability about the annual cycle is also improved. The potential benefit of applying frequency-dependent nudging to the ocean component of a coupled atmosphere–ocean model is discussed.     

Assimilation of Feng-Yun-3B Satellite Microwave Humidity Sounder Data over Land

The ECMWF has been assimilating Feng-Yun-3B(FY-3B) satellite microwave humidity sounder(MWHS) data over ocean in an operational forecasting system since 24 September 2014. It is more difficult, however, to assimilate microwave observations over land and sea ice than over the open ocean due to higher uncertainties in land surface temperature, surface emissivity and less effective cloud screening. We compare approaches in which the emissivity is retrieved dynamically from MWHS channel 1 [150 GHz(vertical polarization)] with the use of an evolving emissivity atlas from 89 GHz observations from the MWHS onboard NOAA and EUMETSAT satellites. The assimilation of the additional data over land improves the fit of short-range forecasts to other observations, notably ATMS(Advanced Technology Microwave Sounder) humidity channels, and the forecast impacts are mainly neutral to slightly positive over the first five days. The forecast impacts are better in boreal summer and the Southern Hemisphere. These results suggest that the techniques tested allow for effective assimilation of MWHS/FY-3B data over land.     

GTS的温盐资料在BCC_GODAS中的同化结果分析

分析了从GTS(全球无线通讯系统)获得的2002-2007年海洋温盐观测资料在国家气候中心第2代全球海洋资料同化系统(BCC_GODAS 2.0)中的同化结果。与SODA(简易海洋同化数据)资料的比较表明:GTS中的海洋温盐资料同化对模式温盐场的改进之处主要表现在混合层暖区的范围和中心强度、温跃层中温度槽脊的深度、温跃层附近的温度梯度以及盐度高、低值区的范围和中心强度等方面,同化后全球温盐场的均方根误差得到一定程度的降低。挑选位于不同海区的单点温盐廓线与ARGO(地转海洋学实时观测阵)观测作了进一步比较,结果表明:大多数情况下,同化后温盐廓线的均方根误差得到明显降低,模拟的温盐垂向分布特征也更为准确。与TAO(热带大气海洋观测网)资料的比较也同样表明:同化后的温盐场特征会得到一定程度改善。     

A Regional Ocean Reanalysis System for Coastal Waters of China and Adjacent Seas

A regional ocean reanalysis system for the coastal waters of China and adjacent seas has been developed by the National Marine Data and Information Service(NMDIS).It produces a dataset package called CORA (China ocean reanalysis).The regional ocean model used is based on the Princeton Ocean Model with a generalized coordinate system(POMgcs).The model is parallelized by NMDIS with the addition of the wave breaking and tidal mixing processes into model parameterizations.Data assimilation is a sequential three-dimensional variational(3D-Var) scheme implemented within a multigrid framework.Observations include satellite remote sensing sea surface temperature(SST),altimetry sea level anomaly(SLA),and temperature/salinity profiles.The reanalysis fields of sea surface height,temperature,salinity,and currents begin with January 1986 and are currently updated every year. Error statistics and error distributions of temperature,salinity and currents are presented as a primary evaluation of the reanalysis fields using sea level data from tidal gauges,temperature profiles,as well as the trajectories of Argo floats.Some case studies offer the opportunity to verify the evolution of certain local circulations.These evaluations show that the reanalysis data produced provide a good representation of the ocean processes and phenomena in the coastal waters of China and adjacent seas.     

An ocean data assimilation system in the Indian Ocean and west Pacific Ocean

The development and application of a regional ocean data assimilation system are among the aims of the Global Ocean Data Assimilation Experiment. The ocean data assimilation system in the regions including the Indian and West Pacific oceans is an endeavor motivated by this goal. In this study, we describe the system in detail. Moreover, the reanalysis in the joint area of Asia, the Indian Ocean, and the western Pacific Ocean(hereafter AIPOcean) constructed using multi-year model integration with data assimilation is used to test the performance of this system. The ocean model is an eddy-resolving,hybrid coordinate ocean model. Various types of observations including in-situ temperature and salinity profiles(mechanical bathythermograph, expendable bathythermograph, Array for Real-time Geostrophic Oceanography, Tropical Atmosphere Ocean Array, conductivity–temperature–depth, station data), remotely-sensed sea surface temperature, and altimetry sea level anomalies, are assimilated into the reanalysis via the ensemble optimal interpolation method. An ensemble of model states sampled from a long-term integration is allowed to change with season, rather than remaining stationary. The estimated background error covariance matrix may reasonably reflect the seasonality and anisotropy. We evaluate the performance of AIPOcean during the period 1993–2006 by comparisons with independent observations, and some reanalysis products. We show that AIPOcean reduces the errors of subsurface temperature and salinity, and reproduces mesoscale eddies. In contrast to ECCO and SODA products, AIPOcean captures the interannual variability and linear trend of sea level anomalies very well. AIPOcean also shows a good consistency with tide gauges.     

GNSS反演资料在GRAPES_Meso三维变分中的应用

为了进一步提高GRAPES_Meso的分析和预报效果,该文在GRAPES_Meso三维变分同化系统中建立了同化GNSS/RO反演的大气资料的观测算子,实现了对GNSS/RO反演的大气资料的同化应用,并通过2013年7月1个月的同化和预报试验分析了GNSS/RO反演大气资料对GRAPES_Meso模式系统分析和预报的影响。结果表明:增加了GNSS/RO反演大气资料的同化后,GRAPES_Meso位势高度场的分析误差明显减小,平均分析误差减小约8%,预报误差略有减小,平均预报误差减小约1%;湿度场的分析误差和预报误差变化不明显,常规观测资料稀少的青藏高原地区的降水预报技巧有所提高,小雨到大雨的ETS (equitable threat score) 评分提高约0.01,对全国及其他分区的降水预报技巧总体上有正效果。     

FY-3A微波资料在“莫拉克”台风预报中的同化试验

我国新一代极轨气象卫星FY-3A于2008年5月26日发射成功,携带的微波垂直探测仪与NOAA系列卫星的ATOVS性能相似。为研究微波垂直探测仪资料在台风数值预报中的作用,实现我国FY-3A卫星的微波探测资料的直接同化,达到改进台风预报的目的,利用FY-3A微波探测资料,以WRF-3DVar系统为基础,针对2009年第8号台风"莫拉克"路径和强度预报,开展数值预报直接同化技术研究。试验结果表明,直接同化FY-3A微波资料对数值模式初始场改进要优于仅仅同化常规观测资料,对缺乏观测资料的海洋上改进尤为明显,模式初始场更加合理地反映海上台风环流形势以及温湿条件,海上台风的模式初始位置也得到了校正;经过FY-3A微波资料三维变分直接同化后,区域中尺度模式对台风路径预报效果具有积极的改善作用。     

A New Global Four-Dimensional Variational Ocean Data Assimilation System and Its Application

A four-dimensional variational data assimilation （4DVar） system of the LASG/IAP Climate Ocean Model, version 1.0 （LICOM1.0）, named LICOM-3DVM, has been developed using the three-dimensional variational data assimilation of mapped observation （3DVM）, a 4DVar method newly proposed in the past two years. Two experiments with 12-year model integrations were designed to validate it. One is the assimilation run, called ASSM, which incorporated the analyzed weekly sea surface temperature （SST） fields from Reynolds and Smith （OISST） between 1990 and 2001 once a week by the LICOM-3DVM. The other is the control run without any assimilation, named CTL. ASSM shows that the simulated temperatures of the upper ocean （above 50 meters）, especially the SST of equatorial Pacific, coincide with the Tropic Atmosphere Ocean （TAO） mooring data, the World Ocean Atlas 2001 （WOA01） data and the Met Office Hadley Centre＇s sea ice and sea surface temperature （HadISST） data. It decreased the cold bias existing in CTL in the eastern Pacific and produced a Nifio index that agrees with observation well. The validation results suggest that the LICOM-3DVM is able to effectively adjust the model results of the ocean temperature, although it＇s hard to correct the subsurface results and it even makes them worse in some areas due to the incorporation of only surface data. Future development of the LICOM-3DVM is to include subsurface in situ observations and satellite observations to further improve model simulations.     

CYGNSS海表风场观测数据验证及其可能的应用

海表面风场可以用于获取许多大气和海洋现象的信号,高质量、高时空分辨率的海表面风场数据产品将有利于海洋-大气动力过程的研究.本文使用全球热带系泊浮标阵列计划（Global Tropical Moored Array Programs）的锚定浮标风场数据和西沙通量塔气象观测资料验证了Cyclone Global Navigation Satellite System （CYGNSS）的35°N~35°S海面遥感风场观测数据.结果表明,CYGNSS海表面风场与实测资料存在着2.17 m/s左右的平均均方根误差（RMSD）,它可能源于观测数据和卫星遥感资料的观测误差,以及两者在空间和时间上未严格匹配而引起的代表性误差.另外,CYGNSS海表面风速的时间演变与实测资料非常一致,展现了CYGNSS在研究海洋-大气能量和动量交换过程方面的潜在应用价值.本文使用Madden-Julian Oscillation （MJO）和赤道东部印度洋上升流事件作为两个个例,说明了CYGNSS海表面风场资料的潜在应用价值.     

风云三号微波观测资料的海雾同化模拟

数值模式边界层物理过程和初值场条件的欠缺是海雾模拟准确率偏低的主要原因。本文为改进模式初始场,开展针对海雾模拟的卫星观测资料同化试验,将质量控制和偏差订正后的FY-3A卫星微波湿度计（MWHS）和微波温度计（MWTS）的优选通道数据,经3DVar（Three-dimensional variational data assimilation）进入WRF模式以试验其对黄、渤海海雾模拟的影响。通过分析静止气象卫星检测到的海雾区模拟大气温、湿场同化分析增量,发现代表环境场条件的海雾类型及模式对其模拟能力的差异,显著影响了同化效果,表现为同化对模式模拟能力较强的平流冷型海雾改进明显,对模拟效果不甚理想的非典型混合过程中的暖型海雾阶段则基本没有改进效果。为寻找原因,对包括海雾区低层大气模拟场逆温结构在内的温湿度场与邻近探空观测进行了对比,分析了随时间演变的海雾格点温、湿场同化分析增量,发现冷型海雾区格点同化分析增量能弥补观测—模拟差异,使气温调减,相对湿度调增,同时水汽和液态水也出现负相关的变化,边界层相关热力动力场同化分析增量在垂直方向也有配合迹象,相比而言,主体是暖型海雾的非典型过程则未见此类现象和其他的有益调整迹象。     

利用ARGO资料改进海洋资料同化和海洋模式中的物理过程

国际ARGO计划的实施每年将可提供多达 10万个剖面 (0～ 2 0 0 0m水深 )的海水温度和盐度资料 ,这些资料的获取无疑将会大大促进海洋和大气科学的发展 ,使人们加深对海洋过程的了解 ,揭示海 气相互作用的机理 ,为长期天气预报和短期气候预测提供模式初始场 ,提高长期天气预报和短期气候预测的能力。如何利用这些资料开展研究工作以及在实际业务中应用这些资料是目前大气和海洋科学界的一个前沿课题。本研究将ARGO浮标资料引入了国家气候中心的NCC GODAS同化系统 ,结果分析表明 ,同化ARGO资料后所得到的海温场在三大洋中不仅在温度数值的大小 ,而且在分布形式方面都与观测场具有较好的一致性 ,可以很好地反映出观测到的冬季和夏季海温的分布形式以及海温的季节变化特征和异常特征。本研究还应用最新的ARGO海洋观测资料 ,通过建立新的热带西太平洋次表层海温参数化方案 ,改进了Zebiak Cane(1987)海洋模式 (ZC模式 ) ,克服了ZC模式几乎没有模拟赤道西太平洋表层和次表层海温变化能力的缺陷。在ZC模式中引用新的次表层海温参数化方案后 ,在赤道西太平洋不仅次表层海温的模拟得到了改善 ,对海面温度异常的模拟也有了较大的改进 ,不仅模拟出了赤道西太平洋表层和次表层海温异常的年际变化特征 ,也模拟出了与观     

Sea-Ice Concentration Multivariate Assimilation for the Canadian East Coast in a Coupled Sea Ice–Ocean Model

The present study focuses on the impact of ocean state (i.e., salinity and temperature) updates on the sea-ice analysis and short-term forecast in an assimilative sea ice–ocean coupled system. A relatively simple sea-ice assimilation scheme was applied to the sea ice–ocean coupled North Atlantic Nucleus for European Modelling of the Ocean (NEMO) system with a focus on the Canadian East Coast. In this assimilation scheme the ocean state was updated directly based on the correlations between the model's sea-ice concentration and the upper ocean salinity and temperature. These correlations were based on a limited time ensemble generated by applying random perturbations to the atmospheric forcing fields. High deviations in the sea-ice conditions were found along the ice edge, implying that the sea-ice edge position is sensitive to small atmospheric forcing variations. Assimilation runs with and without ocean state updates (i.e., sea-ice concentration nudging) were conducted and compared for the winter of 2002. Both continuous and intermittent assimilation schemes were examined. In a continuous sea-ice assimilation experiment, the ocean direct update is unnecessary. When the sea-ice updates are introduced intermittently the ocean state has to be altered to accommodate them, or they will be rapidly diminished by the model's dynamics. The correlations between sea-ice concentration and ocean salinity and temperature based on the first 15 days of January were used for corrections during the entire winter season when, in addition to thermodynamic processes, dynamic processes are responsible for, and even dominate, sea-ice evolution on the Labrador and Newfoundland shelves. This was an adequate choice as was demonstrated by the results of the study which showed that the experiments with ocean state adjustments generated more accurate short-term sea-ice forecasts.     

Argo大洋观测资料的同化及其在短期气候预测和海洋分析中的应用

国际Argo(Array for Real-time Geostrophic Oceanography)计划的实施,提供了前所未有的全球深海大洋0～2000 m水深范围内的海水温度和盐度观测资料,在大气和海洋科研业务中应用这一全新的资料,是深入认识大气和海洋变异、提高我国气候预测、海洋监测分析和预报能力的一个关键所在.通过开发非线性温—盐协调同化方案和利用同化高度计资料来调整模式的温度和盐度场,建立了可同化包括Argo等多种海洋观测资料的全球海洋资料变分同化系统,提高了对全球海洋的监测分析能力.实现了海洋资料同化系统与全球海气耦合模式的耦合,显著提高了短期气候预测水平.利用Argo资料改进了海洋动力模式中的物理过程参数化方案,有效提高了海洋模式对真实大洋的模拟能力和对厄尔尼诺/拉尼娜的预测能力.开发了利用Argo浮标漂流轨迹推算全球海洋表层和中层流的方法,提高了推算的全球表层流、中层流资料质量,有效弥补了洋流观测的匮乏.     

4D-variational assimilation over the Bay of Bengal using the Adjoint of a Variable Resolution Global Ocean Model

Summary The adjoint technique has been widely used over the last two decades in applications involving the dynamics of the Atmosphere and Ocean. This present study applies this technique in the assimilation of oceanic parameters such as temperature and net-surface heat flux, in a variable resolution Ocean General Circulation Model (OGCM). The National Center for Environmental Prediction (NCEP) 1999 July–August global data (sea-surface temperature and derived surface heat flux) and vertical temperature profiles collected from ORV Sagar Kanya and INS Sagardhwani during the Bay of Bengal Monsoon Experiment (BOBMEX 1999) are used for assimilation over the Bay of Bengal. The net-surface heat flux and temperature obtained after assimilation show good agreement with their respective seasonal patterns over the whole Indian Ocean, though there are some differences in the numerical values of heat fluxes. The vertical profiles of temperature after assimilation over the Bay of Bengal show strong resemblance with the corresponding ship observations, emphasizing the importance of the adjoint technique in ocean data assimilation.     

基于自动站资料的WRF-EnSRF陆面同化系统的效果检验及应用

基于集合平方根滤波方法(En SRF)同化方法和NOAH陆面模式的WRF-En SRF陆面同化系统,同化了江苏省70个自动站资料进行试验,研究加入不同的同化资料(地表温度、10 cm土壤温度、20 cm土壤温度)及初始扰动强度的大小对陆面数据同化系统性能的影响,以及对不同区域(降水大值区和降水小值区)的分析场进行效果对比,并且检验了同化系统在一次典型的梅雨锋暴雨的同化效果,证明了这个系统的有效性和可行性。对于资料选取试验,比较全场平均的同化时刻分析场模拟观测相对真实观测的均方根误差可以得到:同化地表温度资料并且初始扰动强度1 K的时候同化效果最理想。对于选定的降水大值区和降水小值区来讲,降水大值区的土壤温度和土壤湿度分析场更加接近于真实场。运用于一次梅雨锋暴雨的同化实验,对于最后一个同化时次的分析场作为背景场做集合预报,最终证明预报结果是有效的。土壤温度、土壤湿度、地表温度和近地面风场的预报结果都较用NCEP再分析资料直接做预报作为控制试验的结果有不同程度的改进。这说明该系统应用于实际同化中的性能较为良好,可以应用于实际土壤湿度与温度的预报。     

Impacts of XBT, TAO, Altimetry and ARGO Observations on the Tropical Pacific Ocean Data Assimilation

This study aims at assessing the relative impacts of four major components of the tropical Pacific Ocean observing system on assimilation of temperature and salinity fields. Observations were collected over a period between January 2001 through June 2003 including temperature data from the expendable bathythermographs (XBT), thermistor data from the Tropical Ocean Global Atmosphere Tropical Atmosphere-Ocean (TOGA-TAO) mooring array, sea level anomalies from the Topex/Poseidon and Jason-1 altimetry (T/P-J),and temperature and salinity profiles from the Array for Real-time Geostrophic Oceanography (ARGO) floats.An efficient three-dimensional variational analysis-based method was introduced to assimilate the above data into the tropical-Pacific circulation model. To evaluate the impact of the individual component of the observing system, four observation system experiments were carried out. The experiment that assimilated all four components of the observing system was taken as the reference. The other three experiments were implemented by withholding one of the four components. Results show that the spatial distribution of the data influences its relative contribution. XBT observations produce the most distinguished effects on temperature analyses in the off-equatorial region due to the large amount of measurements and high quality.Similarly, the impact of TAO is dominant in the equatorial region due to the focus of the spatial distribution.The Topex/Poseidon-Jason-1 can be highly complementary where the XBT and TAO observations are sparse.The contribution of XBT or TAO on the assimilated salinity is made by the model dynamics because no salinity observations from them are assimilated. Therefore, T/P-J, as a main source for providing salinity data, has been shown to have greater impacts than either XBT or TAO on the salinity analysis. Although ARGO includes the subsurface observations, the relatively smaller number of observation makes it have the smallest contribution to the assimilation system.