全天候卫星微波观测资料变分同化研究进展

云和降水区微波观测包含大量与天气系统,特别是台风、暴雨等灾害性天气系统发生发展密切相关的大气信息,因此微波资料的全天候同化应用成为当前数值预报领域的热点研究问题。过去20年间,全球几大数值预报中心逐步开展了全天候同化技术的研究和业务应用,证实了全天候卫星微波观测资料能够改进模式中的质量、风场、湿度以及云和降水场的初始信息,从而改进数值预报模式的预报效果。通过梳理和评述全天候卫星微波观测资料同化方法,分析其中的关键技术问题和目前存在的困难和挑战,为未来在我国数值天气预报领域开展全天候同化研究提供依据。随着我国新一代数值天气预报模式的发展应用,加强我国全天候资料同化技术的研究将会在业务中发挥更大的科学效益和应用效益。     

论地震数值预报——关于我国地震预报研究发展战略的思考

随着中国工业化和社会现代化进程的加快 ,中国城市化进程必将进一步加速。如何减轻地震灾害的问题正变得日益严峻。尽管地震预测是一个举世公认的国际性科学难题 ,但在强化各种减轻地震灾害措施的同时 ,仍须大力推进地震预测研究。为此 ,需要打破长期徘徊在以地震前兆异常监测为基础的经验性预测局面 ,把注意力尽快转向研究以动力学为基础的数值预报。以GPS为代表的空间对地观测技术 ,岩石圈巨型高分辨率宽频带流动地震台阵观测技术以及数值模拟技术已经为地震数值预报研究提供了前所未有的技术基础。以地震数值预报为目标的GPS阵列地壳形变连续观测 ,高分辨率地壳上地幔结构探测 ,地壳动力学 ,地震孕育和破裂过程的理论、模拟试验和实际观测 ,数据同化和计算软件的开发应成为今后研究的重点。现在的问题是 ,需要积极借助数值天气预报的经验 ,强化多学科 ,多部门的组织协调 ,尽早在有条件的地区开展地震动力学数值预报的科学试验。地震数值预报研究必将极大地促进中国地震科学基础研究和地震预报的进一步发展。     

资料同化在空气质量预报中的应用

随着数值模式的不断完善和观测技术的不断提高，资料同化逐渐成为能够进一步提高数值预报水平的一种有效方法。20世纪70年代，资料同化开始引入空气质量预报领域，成为当前大气环境科学研究的一个新方向。简要介绍了资料同化的含义，较详细地介绍了Kalman滤波法、四维变分同化法、牛顿松弛法的基本思想和优缺点，重点阐述了国内外资料同化在空气质量预报中的研究应用情况，最后指出资料同化应用于空气质量预报时存在的问题和今后的研究方向。     

青藏高原气象学的研究进展和问题

分高原天气学、高原气候学、高原及邻近地区的大气环流、以及高原数值预报和模拟四方面简要回顾了新中国成立以来我国（也兼及国外）青藏高原气象学的主要进展，也提出今后研究中应注意的有关问题。     

Kalman滤波在气象数据同化中的发展与应用

气象学领域各种观测(特别是遥感遥测等非常规观测)数据的大量增多和数值天气预报模式的不断进步,推动气象数据同化技术不断发展。回顾了Kalman滤波在气象数据同化中的引入和几个发展阶段;介绍了Kalman滤波(尤其是简化Kalman滤波和总体Kalman滤波)在气象数据同化中的重要地位和应用进展。     

双校正模式下的大清河流域陆气耦合洪水预报研究

基于数值大气模式WRF、三维变分数据同化WRF-3DVar、河北雨洪模型以及实时校正模型ARMA,在北方半湿润半干旱地区的大清河流域构建了陆气耦合洪水预报系统,并利用2012、2013年发生的3场降雨洪水,对系统的降雨洪水预报结果进行分析。结果表明:雷达反射率与GTS数据的同时同化,可有效改善数值大气模式对中小尺度流域的降雨预报效果,从而降低系统的洪水预报误差,ARMA模型的应用,能够进一步提升系统的洪水预报精度,随着预见期的延长,系统的预报精度下降,但系统在6h预见期内仍表现出较好的应用效果。因此,在数据同化和实时校正的"双校正"模式下,陆气耦合洪水预报系统在延长洪水预报预见期的同时,具有较高的洪水预报精度,具有一定的应用前景。     

地面GPS探测大气的最新进展

综述了利用地面ＧＰＳ站探测大气的技术，包括基本原理，方法，应用，最新进展和结果的介绍以及存在于该技术的误差源的讨论。从地面ＧＰＳ站的观测和表面气象参数可得出天顶湿延迟的估计，再由它和可降水汽量之间的转换关系就可得出ＰＷＶ估计。这种方法求得的ＰＷＶ时间序列不仅能提高天气预报的准确度而且有助于气候变化的研究和数值天气预报模型的改进，这就是新兴学科ＧＰＳ气象学形成的基础。     

云分析预报方法研究进展

云作为地球大气系统的重要组成部分,不仅影响着气候变化和天气系统的发展演变,还与航空活动密切相关,一直以来是空军和民航部门非常关注的气象要素之一。在云探测、资料同化和反演方法发展的基础上,从实际业务保障和数值模式发展需求出发,综述国内外云分析、预报方法和云分析预报系统开发的研究成果,分析各类方法的优势和不足,明确国内外研究的主要差距,并探讨国内未来研究的方向。云分析方法中,探空对云廓线识别较好,卫星可见光和红外资料在云顶信息反演方面优势明显,多普勒雷达能够获取对流层中层和底层的云信息,而毫米波雷达能够很好地反映云三维结构信息,发展潜力巨大。云预报方法中,传统的统计和诊断方法发展较为成熟,而考虑了大气温湿和云微物理状况的大气辐射传输模式正演模拟云顶亮温的方法是未来的发展趋势。加强云探测技术,综合利用云分析预报方法,借鉴国外先进云分析预报系统的设计理念,积极开发我国自主的云分析预报系统,推动天气预报、航空气象保障和数值预报模式的发展将会是我国云研究的重要方面。     

大气廓线物理反演的最优化方法进展

从最优化数学理论角度对大气廓线物理反演以及卫星辐射率资料直接同化中的最优化算法进行了回顾。分析了各种方法的优点和缺点、联系和差别。总结了卫星大气遥感反演问题的求解思路。对大气廓线反演研究中几种主要的目标函数和寻优策略进行了分析,着重分析了目前作为各数值预报中心和卫星数据处理中心业务数值产品核心算法的牛顿非线性迭代法的不足之处,并对其改进途径进行了探讨。引入了Levenberg-Marquardt方法及信赖域方法用于大气廓线反演,使反演算法的收敛性质得到改善。     

中国地壳运动观测网络在地球科学研究中的应用前景

简要地介绍了中国地壳运动观测网络的概况和目前的进展,论述了地壳运动观测网络在地球科学研究中的重要作用和应用前景。其应用潜力主要体现在:①在中国地壳运动监测研究中发挥重要作用,成为中国地震预测预报研究的重要基础;②为地球动力学研究提供重要的依据;③建立维持ITRF地球参考框架、研究地球自转、极移和章动及变化;④精化、加密全国大地网和大地水准面;⑤在气象学中的应用;⑥研究电离层电子浓度及其变化规律;⑦提供精密近实时的GPS轨道参数;⑧为广域差分定位奠定基础。     

Abel inversion for deriving refractivity profile from down-looking GPS radio occultation: simulation analysis

Down-looking (DL) Global Positioning System (GPS) radio occultation can produce an estimate of the atmospheric refractivity profile. The main observations are the bending angle as a function of the impact parameter. DL provides both negative as well as positive elevation angle measurements. Abel inversion can be operated on a profile of partial bending angle found by subtracting the positive elevation measurement from the negative one with the same impact parameter. Abel inversion requires the spherical symmetrical assumption. Basically, partial bending calculation removes the ionospheric bending and hence it is possible to use a single frequency GPS receiver. The current paper introduces a simulation data for the case of a receiver on mountain top. The simulation uses model refractivity from MSISE-90 model as well as radiosonde data. Random noises are added to the bending angle profile before inversion. The result shows that it is possible to produce accurate vertical refractivity profile below the receiver altitude. The calculation of the water vapor profile is also made using temperature profile information from the MSISE-90 model as well as radiosonde. The errors in the retrieved vapor profile are always less than 0.1?mbar.     

基于地基GPS与MODIS遥感影像的映射函数与大气可降水汽分析

为了提高短期和中小尺度灾害性天气可降水汽预报的准确性和实时性,利用西安地基GPS数据和MODIS遥感影像,采用NMF、GMF和VMF1映射函数设计了4种获取天顶干延迟初始值的方案,探讨了映射函数对计算天顶总延迟的影响。结果表明：在卫星高度角为10°时,VMF1映射函数反演效果和精度最佳,GMF映射函数次之,NMF映射函数最弱,而在卫星高度角为15°时,上述3种映射函数效果相当。最后,结合综合水汽含量与可降水汽的关系,通过计算地基GPS和MODIS遥感影像综合水汽含量,得到了地基GPS和MODIS遥感影像可降水汽（PPWV-GPS,PPWV-MODIS）,并拟合了二者之间的线性关系,结果为PPWV-MODIS=1.421 7PPWV-GPS-2.143,相关系数为0.952 1。     

Impact of assimilation of ATOVS temperature and humidity and SSM/I total precipitable water on the simulation of a monsoon depression

The present study explored the effect of assimilation of Advanced TIROS Vertical Sounder (ATOVS) temperature and humidity profiles and Spectral sensor microwave imager (SSM/I) total precipitable water (TPW) on the simulation of a monsoon depression which formed over the Arabian Sea during September 2005 using the Weather Research and Forecast model. The three-dimensional variational (3DVAR) data assimilation technique has been employed for the purpose of assimilation of satellite observations. Statistical scores like “equitable threat score,” “bias score,” “forecast impact,” and “improvement parameter” have been used to examine the impact of the above-mentioned satellite observations on the numerical simulation of a monsoon depression. The diagnostics of this study include verification of the vertical structure of depression, in terms of temperature anomaly profiles and relative vorticity profiles with observations/analysis. Additional diagnostics of the study include the analysis of the heat budget and moisture budget. Such budget studies have been performed to provide information on the role of cumulus convection associated with the depression. The results of this study show direct and good evidence of the impact of the assimilation of the satellite observations using 3DVAR on the dynamical and thermodynamical features of a monsoon depression along with the effect of inclusion of satellite observation on the spatial pattern of the simulated precipitation associated with the depression. The “forecast impact” parameter calculated for the wind speed provides good evidence of the positive impact of the assimilation of ATOVS temperature and humidity profiles and SSM/I TPW on the model simulation, with the assimilation of the ATOVS profiles showing better impact in terms of a more positive value of the “forecast impact” parameter. The results of the study also indicate the improvement of the forecast skill in terms of “equitable threat score” and “bias score” due to the assimilation of satellite observation.     

高精度GPS观测中的负荷效应

讨论了高精度GPS观测中的负荷效应问题。根据地表质量负荷理论的Green函数方法，计算了地球表面质量负荷（包括大气、非潮汐海洋和陆地水）对中国部分IGS台站垂直位移的影响。其中，大气负荷对各台站的影响比较一致，变化幅度约为20 mm，非潮汐海洋对沿海台站的负荷影响显著，其变化幅度达到10 mm，水负荷对各台站的影响差异较大，最大的变化幅度约20 mm；总的负荷效应达到厘米量级，与GPS观测结果的比较说明在高精度的GPS观测的时间序列中，存在明显的季节性变化，而该变化与表面负荷有着直接的关系。特别是对于房山和乌鲁木齐2个台站，GPS观测的垂直变化中的季节项基本上就是来自于地球表面的质量负荷；而昆明和拉萨这2个台站尽管位于板块运动活跃区域，垂直季节性变化也主要源于地表负荷。因此扣除这些负荷影响的GPS观测资料将更有利于研究地壳的运动特征。     

A global positioning measurement system for regional geodesy in the caribbean

Low cost, portable receivers using signals from satellites of the Global Positioning System (GPS) will enable precision geodetic observations to be made on a large scale. A number of important geophysical questions relating to plate-motion kinematics and dynamics can be addressed with this measurement capability. We describe a plan to design and validate a GPS-based geodetic system, and to demonstrate its capability in California, Mexico and the Caribbean region. The Caribbean program is a prototype for a number of regional geodetic networks to be globally distributed.

In 1985, efforts will be concentrated on understanding and minimizing error sources. Two dominant sources of error are uncertainties in the orbit ephemeris of the GPS satellites, and uncertainties in the correction for signal delay due to variable tropospheric water vapor. Orbit ephemeris uncertainties can be minimized by performing simultaneous satellite observations with GPS receivers at known (fiducial) points. Water vapor corrections can be made by performing simultaneous line-of-sight measurements of integrated water vapor content with ground-based water vapor radiometers. Specific experiments to validate both concepts are outlined.

Caribbean measurements will begin in late 1985 or early 1986. Key areas of measurement are the northern strike-slip boundary, and the western convergent boundary. Specific measurement plans in both regions are described.     

A new ensemble-based data assimilation algorithm to improve track prediction of tropical cyclones

This paper proposes a new ensemble-based algorithm that assimilates the vertical rain structure retrieved from microwave radiometer and radar measurements in a regional weather forecast model, by employing a Bayesian framework. The goal of the study is to evaluate the capability of the proposed technique to improve track prediction of tropical cyclones that originate in the North Indian Ocean. For this purpose, the tropical cyclone Jal has been analyzed by the community mesoscale weather model, weather research and forecasting (WRF). The ensembles of prognostic variables such as perturbation potential temperature (θ, k), perturbation geopotential (?, m²/s²), meridional (U) and zonal velocities (V) and water vapor mixing ratio (q _v , kg/kg) are generated by the empirical orthogonal function technique. An over pass of the tropical rainfall-measuring mission (TRMM) satellite occurred on 06th NOV 0730 UTC over the system, and the observations from the radiometer and radar on board the satellite(1B11 data products) are inverted using a combined in-home radiometer-radar retrieval technique to estimate the vertical rain structure, namely the cloud liquid water, cloud ice, precipitation water and precipitation ice. Each ensemble is input as a possible set of initial conditions to the WRF model from 00 UTC which was marched in time till 06th NOV 0730 UTC. The above-mentioned hydrometeors from the cloud water and rain water mixing ratios are then estimated for all the ensembles. The Bayesian filter framework technique is then used to determine the conditional probabilities of all the candidates in the ensemble by comparing the retrieved hydrometeors through measured TRMM radiances with the model simulated hydrometeors. Based on the posterior probability density function, the initial conditions at 06 00 UTC are then corrected using a linear weighted average of initial ensembles for the all prognostic variables. With these weighted average initial conditions, the WRF model has been run up to 08th Nov 06 UTC and the predictions are then compared with observations and the control run. An ensemble independence study was conducted on the basis of which, an optimum of 25 ensembles is arrived at. With the optimum ensemble size, the sensitivity of prognostic variables was also analyzed. The model simulated track when compared with that obtained with the corrected set of initial conditions gives better results than the control run. The algorithm can improve track prediction up to 35 % for a 24 h forecast and up to 12 % for a 54 h forecast.     

Study About Ocean Eddy Effect on Strong Convection in Local Atmosphere over the Kuroshio Extension Region

The Kuroshio Extension （KE） is the key area where the water heats the atmosphere in the northwestern Pacific Ocean in winter. Previous studies show that the active eddies in the KE area can affect sea surface temperature and thus sea surface winds. The present study reviewed the progress about the influences of the eddies on local atmosphere in recent years. Analysis and comparison were made especially for the achievements from shipboard sounding data, satellite observations and numerical experiments. Based on the geostrophic adaptation theory involved in atmospheric anomalies induced by the eddies, the following new scientific deductions were suggested: Air pressure adjustment mechanism dominated in the atmospheric response to eddies under the conditions of weaker wind speed over the eddies. The influence of eddies was often limited in the atmospheric boundary layer. On the other hand, vertical mixing mechanism played a major role in the response of the atmosphere to warm (cold) eddies when air moved faster over the eddies. Surface wind speed increased (decreased) over the warm (cold) water. Significant wind convergence took place downwind the warm water, and large amount of water vapor was transported also downwind from the warm water surface. The positive feedback between water vapor condensation and rising air forced by the surface convergence provided necessary conditions for the development of strong convection in atmosphere. These deductions will be conducive to further depicting the impact of oceanic eddies on the atmosphere quantitatively.     

空间GPS无线电掩星反演大气参数方法及其应用

给出了空间无线电掩星反演大气参数的原理及其误差估计方法,介绍了国内外在该领域的研究进展。并针对该技术的一些特点,指出了若干尚需解决的问题。     

呼伦贝尔市2016年3月31日-4月2日暴雪过程天气学特征研究

利用常规观测资料及NCEP 1.0°×1.0°再分析资料,对2016年3月31日至4月2日发生于呼伦贝尔市的一次暴雪天气过程的环流形势、成因机制进行分析,以期总结出此次暴雪的特殊之处,为今后的暴雪天气预报提供可参考的经验。从天气学角度详细地分析此次暴雪过程的高低空影响系统的生消演变,特别对产生极端降水的水汽输送条件、水汽辐合情况、局地水汽聚集以及垂直运动情况进行详尽的分析。研究结果表明：此次暴雪天气过程是由短波扰动叠加低空锋区,导致斜压有效位能释放转变为扰动动能,并促使扰动发展加深为大型涡动而产生;水汽通量分布形态在一定程度上决定了水汽输送效果,“均匀狭长”的分布形态具有更高的水汽输送效率;涡度差动平流与温度平流表明动力因子与热力因子在本次过程中对垂直运动均有显著贡献,而水平散度作为直接反映质量汇集、流失的参量可以更为直观地反映垂直运动的强弱与分布;整层水汽通量散度积分作为与降水强度直接相关的物理量,对于降水量级的反映异常精准,配合中低空气流强度与方向,可对降水落区与时段进行精确判断。