L波段探空仪湿度资料偏差订正及同化试验

L波段探空仪观测资料是基础资料之一,无论在天气预报还是在数值天气预报中都起着重要作用,其资料质量直接影响数值模式同化分析及降水预报准确性。通过对我国3种常用的L波段探空仪观测湿度的偏差特性比较,研发适合该仪器的偏差订正方案,并在GRAPES同化系统中加以试验应用。结果表明:L波段探空仪湿度观测资料与ECMWF再分析湿度场比较有偏干现象。多种偏差订正方案订正结果显示:湿度偏差值比订正前减小,特别是在500 hPa以上层次减小明显。在GRAPES分析同化系统中使用Vomel偏差订正方案,分析偏差减小5%;预报模式个例和连续试验中湿度观测订正后预报降水更接近实况,预报降水检验评分显著提高,故该订正方案在实际应用中表现出积极的正效果。     

中国L波段探空湿度观测资料的质量评估及偏差订正

L波段探空观测资料无论在天气预报还是数值预报中均为最基本和最重要的一类数据,而其湿度观测资料的质量对同化分析及降水预报有直接影响。通过用L波段探空湿度观测资料与不同类型的其他观测反演的湿度资料互校及与NCEP、GRAPES、EC等不同模式分析场为背景的湿度场比较,评估中国L波段探空湿度观测资料的质量状况,对探空湿度资料的质量有了新的认识,为更好地使用该资料提供依据。研究发现中国L波段探空湿度观测资料存在偏干的现象,特别是当背景场湿度大于60%时,观测湿度偏低更加明显。通过分析其偏差特征,找出了适合中国L波段探空湿度观测资料偏差特点的分段函数订正方法。个例试验表明,对探空湿度观测资料的偏差订正后,观测偏差明显减小,订正效果非常显著;模式降水强度预报能力有一定的提高。从连续试验检验的降水预报评分(TS)和预报偏差(Bias)看,中雨和暴雨的预报在探空湿度观测偏差订正后都表现出正效果。     

营口和鞍山城市气候变化对比分析及原因探讨

利用1951—2005年辽宁省沿海轻工业城市营口和重工业城市鞍山的气温、降水、云量、相对湿度及平均风速等气候资料,采用对比分析方法,分析了近55 a来2城市的各类气候变化特征及产生的原因。结果表明:近55 a来营口和鞍山平均气温呈递增趋势,降水、总云量、平均风速和相对湿度均呈递减趋势,只有低云量变化趋势不同,其中营口呈递增趋势,鞍山呈递减趋势,鞍山平均气温递增趋势及总云量、平均风速和相对湿度递减趋势均强于营口同类气候要素变化;营口和鞍山城市平均气温、相对湿度、低云量和风力气候要素特征趋势变率较大,降水和总云量趋势变率相对较小。     

A practical approach for the assimilation of cloudy infrared radiances and its evaluation using airs simulated observations

Abstract

A variational estimation procedure for the simultaneous retrieval of cloud parameters and thermodynamic profiles from infrared radiances is proposed. The method is based on a cloud emissivity model which accounts for the frequency dependence of cloud absorption and scattering and possible mixed phase situations. An effective cloud top height and emissivity are assumed. Monte Carlo experiments performed in a 1D‐var assimilation context using simulated Atmospheric Infrared Radiance Sounder (AIRS) observations from 100 channels demonstrate the substantial added value, in theory, of cloudy radiance assimilation as opposed to clear‐channel assimilation. Improved temperature and humidity retrievals are obtained for a broad layer above the cloud as well as below cloud level under partial cloud cover conditions. The impact is most pronounced in broken to overcast situations involving mid‐level clouds. In these situations, the effective cloud top height and emissivity are retrieved with estimated rms errors typically lower than 30 hPa and 3%, respectively. Expected relative errors on the retrieved effective particle size are of the order of 30–50%. The methodology is directly applicable to real hyperspectral infrared data upon inclusion, for local estimation, of the cloud parameters in the Canadian 4D‐var assimilation system.     

Use of Surface Observations to Estimate Upper Air Humidity for the Objective Analysis of Relative Humidity over Indian Region

In the present study objective analyses of relative humidity (RH) at surface and at the levels of 850,700 and 500 hPa have been made using Gandin’s (1963) optimum interpolation scheme. As the horizontal resolution of the radiosonde stations is rather inadequate for upper air humidity analysis, a scheme has been developed, following Rasmussen (1982) to estimate the upper air RH from the surface observations like surface RH, present weather and cloud cover. The relative humidities at the levels 850, 700 and 500 hPa were related to the surface observations through three separate regression relations. The RH values at 850, 700 and 500 hPa levels were estimated from the surface RH, cloud coverage and present weather using the above regression relations and subsequently the objective analyses at 00 GMT for the period from 4 July to 8 July 1979, were made using these estimated data along with the observed radiosonde data. Objective analyses were also made for the same period using only the radiosonde data for comparison to study the impact of those estimated data. Root mean square errors were computed for all the five days by interpolating RH at the observing stations from the objectively analysed field and comparing them with the actually observed RH to examine how best the analyses (with and without estimated data) fitted the observations. Lastly they were compared with satellite cloud pictures. This study shows that the estimated upper air RH values have positive impact on the analysis of upper air RH and could be used over radiosonde date sparse region and even over oceanic regions.     

1961～2017年雅鲁藏布江河谷地区夏季气候暖干化趋势

基于1961～2017年青藏高原腹地雅鲁藏布江河谷地区4个站（拉萨、日喀则、泽当和江孜）夏季（6～8月）月平均气温、降水和相对湿度等观测资料,分析了该地区夏季气候年际和年代际演变特征,并探讨了气温、降水和相对湿度在年际和年代际时间尺度上的相互关系以及与总云量和地面水汽压的联系。结果表明:（1）1961～2017年该地区夏季气候出现了暖干化趋势。气温（相对湿度）显著升高（下降）,降水趋势变化不明显;本世纪初气温（相对湿度）均发生了显著的突变。（2）该地区夏季气候因子间在年际和年代际时间尺度上存在密切关系:气温与相对湿度和降水均存在明显的负相关,降水与相对湿度为正相关。（3）该地区夏季气候因子间的年际和年代际变化与同期总云量和地面水汽变化有关。1961～2017年总云量持续减少是气温显著升高的主要原因之一,气温的显著升高和降水变化不明显又造成了相对湿度的显著下降。     

Evaluation of cloud and precipitation parameterization using a single-column model: A TWP-ICE case study

Cloud and precipitation parameterization schemes are evaluated, and their sensitivity to the method and/or parameters used to determine cloud physical processes is examined using a singlecolumn version of the Unified Model (SCUM). In the experiment for TWP-ICE, cloud fraction is overestimated (underestimated) in the upper (lower) troposphere due to the wet (dry) bias. The precipitation rate is well simulated during the active monsoon period, but overestimated during the suppressed monsoon and clear skies periods. In the moist convection scheme, trigger condition and entrainment process affect the lower tropospheric humidity through the impact on convective occurrence frequency and intensity, respectively. Strengthening the trigger condition and using the adaptive entrainment method alleviate the low-level dry bias. In the microphysics scheme, more large-scale precipitation is produced with prognostic rain, due to rain sedimentation considering vertical velocity of rain drop, than with diagnostic rain. Less ice/snow deposition with the prognostic two-ice category results in lower ice water content and upper-level cloud fraction than with the diagnostic splitting method for the twoice category. In the cloud macrophysics scheme, the prognostic cloud fraction and cloud/ice water content scheme produces a larger cloud fraction and more cloud/ice water content than the diagnostic scheme, mainly due to detrainment from moist convection (cloud source) that surpasses the effect of convective heating and drying (cloud sink). This affects temperature by influencing the radiative, convective, and microphysical processes. The experiment with combined modifications in cloud and precipitation schemes shows that interaction between modified moist convection and cloud macrophysics schemes results in more alleviation of the cold bias not only at the lower levels but also at the upper levels.     

L波段与59-701探空系统相对湿度对比分析

通过统计全国探空系统早期换型时获取的59个高空站L波段雷达-GTS1型电子探空仪系统相对于59-701探空系统在1000 hPa到200 hPa之间各规定等压面相对湿度的差值,分析探空系统换型对于相对湿度一致性的影响。结果表明:GTS1型探空仪测得的相对湿度明显比59型探空仪低,且差值随高度增加而增大,近地层二者相对湿度差值低于5%,200 hPa高度二者差值达到20%以上;冬季两套探空系统相对湿度的差值明显大于夏季,且差值随高度分布情况不尽相同。分析表明:相对湿度差值的变化除了与探空仪施放过程中外界温度变化相关,还与湿度变化趋势和变化幅度密切相关,湿度传感器存在湿滞回线和滞后现象。分析还发现,太原无线电一厂与上海长望气象科技有限公司制造的59型探空仪的性能差异不明显,且GTS1型探空仪与59型探空仪的两种湿度传感器在白天受太阳辐射的残余影响差异也不明显。     

两种探空仪观测湿度垂直分布及其应用比较

对2010年8月在云南腾冲利用芬兰Vaisala RS80和低温霜点仪 (Cryogenic Frostpoint Hygrometer,CFH) 两种探空仪测量大气湿度的垂直分布进行对比分析,同时比较它们白天和夜间测量误差的差别,并对国产GTS1,RS80和CFH共3种探空仪测量水汽总量与地基GPS遥测结果进行比较。结果表明:RS80湿度测值在整个对流层比CFH测值偏干 (23.7±18.5)%;因太阳辐射白天RS80偏干较夜间更明显,比夜间偏干 (13.5±14.8)%。而在对流层上层向平流层过渡区域内RS80湿度数据基本无效。CFH在低温、低湿环境下对湿度能有效测量,但在湿度较高的对流层低层测值偏高,导致比较中CFH水汽总量平均比GPS遥测的水汽总量偏高 (4.3±2.0) mm (样本数为11),而RS80,GTS1与GPS的水汽总量差别分别是 (0.2±1.4) mm (样本数为12), (-0.2±2.2) mm (样本数为43)。地基GPS遥测的水汽总量对对流层上层至平流层的水汽变化不敏感。由于RS80测量相对湿度在高空偏低,通过RS80相对湿度测值来确定中、高云结果是偏低的,特别是对6000 m以上的高云判别上,RS80相对湿度的探测几乎很难甄别到云的存在。     

The Impact of Satellite Radiance Data Assimilation within a Frequently Updated Regional Forecast System Using a GSI-based Ensemble Kalman Filter

A regional ensemble Kalman filter (EnKF) data assimilation (DA) and forecast system was recently established based on the Gridpoint Statistical Interpolation (GSI) analysis system. The EnKF DA system was tested with continuous threehourly updated cycles followed by 18-h deterministic forecasts from every three-hourly ensemble mean analysis. Initial tests showed negative to neutral impacts of assimilating satellite radiance data due to the improper bias correction procedure. In this study, two bias correction schemes within the established EnKF DA system are investigated and the impact of assimilating additional polar-orbiting satellite radiance is also investigated. Two group experiments are conducted. The purpose of the first group is to evaluate the bias correction procedure. Two online bias correction methods based on GSI 3DVar and EnKF algorithms are used to assimilate AMSU-A radiance data. Results show that both variational and EnKF-based bias correction procedures effectively reduce the observation and background radiance differences, achieving positive impacts on forecasts. With proper bias correction, we assimilate full radiance observations including AMSU-A, AMSU-B, AIRS, HIRS3/4, and MHS in the second group. The relative percentage improvements(RPIs) for all forecast variables compared to those without radiance data assimilation are mostly positive, with the RPI of upper-air relative humidity being the largest. Additionally, precipitation forecasts on a downscaled 13-km grid from 40-km EnKF analyses are also improved by radiance assimilation for almost all forecast hours.     

A preliminary investigation of winter air pollution potential at fort simpson,northwest territories

Based on radio‐sonde data collected during nine consecutive days in February 1974, an examination was made of the daytime variation of the winter mixing depth and air pollution potential in the Fort Simpson area of the Mackenzie Valley. Under anti‐cyclonic weather conditions, mixing depths were generally low (< 100 m) or non‐existent. The generation of a mixed layer was found to be primarily associated with extensive low‐ to mid‐level cloud cover and precipitation. During this period of investigation, ventilation coefficients were extremely low. This suggests that the area has a great risk of high pollution potential in winter.     

The sublimation of falling snow over the Mackenzie River Basin

The sublimation of falling snow may be an important component of the atmospheric water budget of the Mackenzie River Basin and many parts of the Arctic. To investigate this issue, a simple sublimation model is used along with surface precipitation observations and sonde data obtained during the autumn 1994 Beaufort and Arctic Storms Experiment (BASE). Model results are then compared with actual precipitation measurements at Inuvik and Tuktoyaktuk, sites in Northern Canada, to approximate mass loss due to sublimation. The sublimation results are found to vary in concert with cloud base height, precipitation intensity aloft and the nature of the precipitation. Atmospheric conditions are furthermore examined over a wide range of the Arctic, especially the Mackenzie River Basin, to assess to what degree the results can be generalized. The presence of a relatively dry near-surface layer, a favourable environment for sublimation, is a key feature of most sites during the early autumn storm period. Estimates of sublimational mass losses are found over Inuvik and Tuktoyaktuk using sonde derived cloud base heights and temperature and humidity profiles. Sublimation losses for such sites are found to be of the order of 40–60%, which shows that sublimation is indeed a significant process over the Mackenzie Basin and needs to be well handled in climate models. However, increasing the vertical resolution of the sublimation model to that of climate scales can dramatically affect predicted sublimation amounts; how to properly account for sublimation then remains a difficult task.     

Intercomparison of Humidity and Temperature Sensors: GTS1, Vaisala RS80, and CFH

GTS1 digital radiosonde, developed by the Shanghai Changwang Meteorological Science and Technology Company in 1998, is now widely used in operational radiosonde stations in China. A preliminary comparison of simultaneous humidity measurements by the GTS1 radiosonde, the Vaisala RS80 radiosonde, and the Cryogenic Frostpoint Hygrometer (CFH), launched at Kunming in August 2009, reveals a large dry bias produced by the GTS1 humidity sensor. The average relative dry bias is in the order of 10% below 500 hPa, increasing rapidly to 30% above 500 hPa, and up to 55% at 310 hPa. A much larger dry bias is observed in the daytime, and this daytime effect increases with altitude. The GTS1 radiosonde fails to respond to humidity changes in the upper troposphere, and sometimes even in the middle troposphere. The failure of GTS1 in the middle and upper troposphere will result in significant artificial humidity shifts in radiosonde climate records at stations in China where a transition from mechanical to digital radiosondes has occurred. A comparison of simultaneous temperature observations by the GTS1 radiosonde and the Vaisala RS80 radiosonde suggests that these two radiosondes provide highly reproducible temperature measurements in the troposphere, but produce opposite biases for daytime and nighttime measurements in the stratosphere. In the stratosphere, the GTS1 shows a warm bias (<0.5 K) in the daytime and a relatively large cool bias (-0.2 K to -1.6 K) at nighttime.     

L波段电子探空仪相对湿度误差研究及其应用

为了提高高空大气探测数据准确度, 我国从2002年1月开始推广使用L波段雷达-电子探空仪探测系统, 用携带碳湿敏元件的数字式电子探空仪代替用肠膜测湿元件的59型探空仪进行相对湿度探测。但大量的实测探空资料表明:相对湿度探空曲线仍然存在较大误差。利用能测到-30℃低温的高精度湿度校准设备在-30~30℃试验温度范围内对碳湿敏元件进行大量静态测试, 在进一步了解碳湿敏电阻校准线随温度变化的特征基础上, 结合实际探空资料, 修正工厂的相对湿度订正原理和公式, 从而可以提高碳湿敏电阻在高湿端的测量准确度, 提高判断云层垂直位置的准确性, 进而提高L波段探空仪温度测量的精度。     

气候模式分辨率对气溶胶气候效应模拟结果的影响

气候模式分辨率作为影响模式模拟结果的重要因素,其对气溶胶与云相互作用的影响尚未全面认识。利用公共大气模型CAM5.3在3种分辨率（2°、1°、0.5°）下,分别采用2000年和1850年气溶胶排放情景进行试验,检验提高分辨率是否能改进气候模式的模拟能力,分析不同分辨率下气溶胶气候效应的异同,探索模式分辨率对气溶胶气候效应数值模拟结果的影响。通过观测资料与模式结果对比发现,提高分辨率可以明显改进模式对总云量、云短波辐射强迫的模拟能力,0.5°分辨率下模拟结果与观测更接近,其他变量并无明显改善。在不同分辨率下,全球平均的气溶胶气候效应较为一致,总云量、云水路径均增加,云短波和长波辐射强迫均加强,而云顶的云滴有效半径和降水均减小,地面气温降低。不同分辨率下,气溶胶增加引起的气溶胶光学厚度、云水路径、地面温度、云短波和长波辐射强迫变化的纬向平均分布相似但大小存在差异;而降水和云量变化的纬向分布与大小均存在较大差异,在区域尺度上还存在较大的不确定性。全球平均而言, 0.5°分辨率下气溶胶的间接辐射强迫相比1°分辨率下的结果降低了2.5%,相比2°分辨率下的结果降低了6.4%。提高模式分辨率可以部分改进模式模拟能力,同时,气溶胶的间接效应随着模式分辨率的提高而减弱。但气溶胶引起的云量、降水的变化在不同分辨率下差异较大,存在较大的不确定性。      

三维变分和物理初始化方法相结合同化多普勒雷达资料的试验研究

以美国新近研发的天气研究预报模式(WRF)配置的三维变分(3D-Var)同化系统WRF 3D-Var为平台,结合物理初始化方法(Physical Initialization,简称PI)来同化多普勒雷达径向风和回波强度观测资料.其基本做法是首先用物理初始化方法由雷达回波资料估计出比湿、云水混合比和垂直速度,然后用估计的比湿和云水混合比对模式的相应变量进行调整,最后再将估计出的垂直速度作为一种新的观测类型添加到现有的WRF 3D-Var目标函数中,同时以WRF 3D-Var提供的方法直接同化径向风.针对2002年6月19日的一次强对流性降水过程和2003年7月5日的一次梅雨锋暴雨过程进行了一组同化多普勒雷达径向风和回波资料的试验研究.同化结果表明:分析变量的增量场和观测的雷达回波有很好的对应关系.在雷达回波区,有正的比湿增量、云水含量增量和垂直速度增量,并且水平风增量在此辐合;在没有雷达回波的地方有负的垂直速度增量.预报结果表明,调整云水含量对降水预报改善不明显,调整比湿对降水预报改进明显,直接用物理初始化估计出的垂直速度替代模式的初始垂直速度,对降水预报改进不明显,但以新的方案同化雷达资料能有效地缩短模式的起转时间(spin-up time),明显改进短时降水预报.     

探空和飞机观测资料联合同化对台风“苏迪罗”（2015）数值模拟的影响研究

基于中尺度数值天气预报模式WRF和WRFDA三维变分同化系统,针对2015年台风“苏迪罗”二次登陆过程,利用NCEP ADP提供的无线电探空、飞机报观测资料,同时同化两种资料,并进行36 h模拟预报,从台风移动路径、强度、降水及模式初始场改进等方面分析了同化模拟效果。结果表明,同时同化两种资料能够有效改善台风移动路径、中心附近最低气压的模拟。对浙闽降水关键区24 h降水的TS评分结果表明,该同化方案对台风降水的预报评分有一定程度的改进,尤其是中雨和特大暴雨的改进最明显。对于模式初始场的温度、相对湿度、纬向风、经向风在各高度上的均方根误差,除湿度外均有不同程度的减小。并且同化对于台风结构的模拟也有所调整。      

青藏高原探空大气水汽偏差及订正方法研究

水汽是大气的主要成分和降水的主要物质来源.青藏高原大气水汽分布对区域天气和气候有很大影响,为了探讨探空观测的大气水汽总量（R）资料的可靠性,本文以地基GPS遥感的大气水汽总量（G）为参照标准,对拉萨（1999~2010年）和那曲（2003年）的R进行对比分析和偏差（R-G）订正.结果表明:近10多年拉萨站R比G明显偏小,偏小程度随使用不同的探空仪而异.GZZ-2型机械探空仪和GTS-1型电子探空仪多年平均的PW偏差分别为-8.8%和-3.9%,随机误差分别为17.6%和13.6%.近10多年PW偏差变化呈减少趋势,这与探空仪性能改进有关.分析发现,青藏高原PW偏差具有明显季节变化和日变化特征,夏季比冬季明显,1200 UTC比0000 UTC明显.拉萨站GZZ-2型和GTS-1型探空仪在1200 UTC多年平均的PW偏差分别为-15.8%和-7.3%,在0000 UTC分别为-1.6%和-0.4%.那曲站GZZ-2型探空仪在1200 UTC和0000 UTC的PW偏差分别为-12.4%和-0.3%.分析还表明,太阳辐射加热与气温的日变化和季节变化是造成高原PW偏差日变化和季节变化的重要原因.据此,提出了高原PW偏差的订正方法,并以拉萨和那曲站为例进行PW偏差订正,订正后的PW系统偏差显著减少,随机误差也相应得到了改善.     

43a来我国城市气候和太阳辐射的变化特征

利用1961-2003年气候观测资料和辐射资料,综合分析了我国10个城市43a的云量、日照百分率、相对湿度、气温、降水以及到达地面的太阳辐射等要素的变化趋势。结果表明：（1）43a来10个城市中,多数城市的总云量减少而低云量增加;（2）各城市的日照百分率呈减少趋势;除乌鲁木齐外,其他城市相对湿度呈减小趋势;各地气温都有所升高;降水量变化的地区性差别较大而总体上变化幅度不明显;（3）除昆明外,各地太阳总辐射呈下降趋势,且总辐射减少主要是由直接辐射减少引起的。由此揭示,43a来我国城市正经历着以变暖变干为主的气候变化过程,人类活动特别是工业化和城市化进程,对气候系统产生了重要影响。     

1979—2012年中国探空相对湿度资料的非均一性检验与订正

利用加拿大环境部气候研究中心研发的PMTred非均一性检验方法,以ERA interim资料作为参考序列,应用中国区域各探空台站详细的元数据信息为主要断点判断依据,对1979—2012年我国125个探空台站各标准等压面月平均探空相对湿度资料进行了非均一性检验和订正。并结合详细的元数据信息分析了造成我国探空相对湿度序列非均一问题形成的主要原因及对资料的影响程度。结果表明：仪器换型、探测系统资料整理计算方式变化以及辐射订正方法改变和探测系统升级等是造成这一时段中国区域月平均探空相对湿度资料不均一的主要原因。其中,2002年以后的L波段雷达 电子探空仪换型造成了相对湿度资料非常明显的不连续问题,这主要是由于早期59 701型探空仪穿云挂水,造成探测到的相对湿度资料明显偏湿。两个观测时次相对湿度序列不均一的台站数和断点数随高度的增加而增加。各标准等压面上月平均相对湿度序列不均一的探空台站平均订正幅度也随着高度的升高而增大,并且订正量为负值的比例在整套订正资料中所占的比例较高,说明我国的探空月平均相对湿度原始观测资料有明显高估的问题。1979—2012年以来,从全国的情况来看,订正前全国850~300 hPa 5个标准等压面全部为相对湿度降低趋势。但是订正后,850~300 hPa这5个标准等压面在1979—2012年相对湿度的变化趋势均不显著。