4D tomographic reconstruction of the tropospheric wet refractivity using the concept of virtual reference station,case study: northwest of Iran

Iran enjoys a variety of climatological conditions. Moreover, numerical weather prediction (NWP) models are not assimilated with the meteorological data in Iran, the country suffering from poor spatial and temporal resolution of radiosonde measurements. These facts make modeling of troposphere impossible using the measurements and NWP. On the other hand, the global positioning system (GPS) has been emerged as a valuable tool for modeling and remote sensing of Earth’s atmosphere. This research is the first attempt to address the tropospheric wet refractivity modeling by GPS measurements in Iran. Changes of topography in the study area are taken into account. As a leading work, virtual reference stations (VRS) are used to fix the rank deficiency of the problem. The model space resolution matrix is used to achieve the optimum spatial resolution of the tomographic model and the optimum number of VRS stations. The accuracy of the developed model (KNTU1) is investigated by deploying radiosonde measurements.     

利用BERNESE 5.0解算地基GPS天顶湿延迟

地基GPS技术已被公认为观测大气水汽的最具潜力手段,而天顶湿延迟（ZWD）是地基GPS解算高精度水汽的关键量。瑞士伯尔尼大学天文研究所开发的BERNESE软件在解算天顶湿延迟方面独树一帜。以香港地基GPS连续运行参考站数据为解算实例,详细介绍了BERNESE软件解算ZWD的基本步骤和相关设置,并对有气象观测文件、数据跨天、跨周和定点解算ZWD情况的特殊设置进行了研究,结果表明BERNESE软件完全胜任ZWD解算工作。     

Impact of surface observations on simulation of rainfall over NCR Delhi using regional background error statistics in WRF-3DVAR model

In this study, efforts are made to improve the simulation of heavy rainfall events over National Capital Region (NCR) Delhi during 2010 summer monsoon, using additional observations from automatic weather stations (AWS). Two case studies have been carried out to simulate the relative humidity, wind speed and precipitation over NCR Delhi in 48-h model integrations; one from 00UTC, August 20, 2010, and the other from 00UTC, September 12, 2010. Several AWS installed over NCR Delhi in the recent past provide valuable surface observations, which are assimilated into state-of-the-art weather research and forecasting (WRF) model using the three-dimensional variational data assimilation (3DVAR). The quality of background error statistics (BES) is a key component in successful 3DVAR data assimilation in a mesoscale model. In this study, the domain-dependent regional background error statistics (RBS) are estimated using National Meteorological Center method in the months of August and September 2010 and then compared with the global background error statistics (GBS) in the WRF model. The model simulations are analyzed and validated against AWS and radiosonde observations to quantify the impact of RBS. The root mean square differences in the spatial distributions of precipitation, relative humidity and wind speed at the surface showed significant differences between both the global and regional BES. Similar differences are also observed in the vertical distributions along the latitudinal cross section at 28.5°N. Model-simulated fields are analyzed at five different surface stations and one upper air station located in NCR Delhi. It is found that in 24-h model simulation, the RBS significantly improves the model simulations in case of precipitation, relative humidity and wind speed as compared to GBS.     

利用ERA5资料进行桂林地区GNSS水汽反演精度分析

由于桂林地区地基GNSS站并未配置气象传感器，致使大量GNSS观测数据无法在大气水汽（PWV）监测中发挥作用.针对这一情况，本文将欧洲中期天气预报中心（ECMWF）最新发布的ERA5再分析资料中测站处的气压和温度气象数据加入到GNSS水汽反演中，并将反演结果与利用地面气象站反演的GNSS水汽做对比，以此评估ERA5在桂林地区反演GNSS水汽的精度和适用性.结果表明：1）以桂林地区2017年10个地面气象站的实测气压和温度数据为参考值，ERA5地表气压和温度的年均偏差分别为-0.35 hPa和0.86 K，年均均方根误差（RMSE）分别为0.65 hPa和1.66 K，该精度可用于GNSS水汽反演；2）以2017年6—7月GNSS利用地面气象站反演的PWV为参考值，ERA5反演的GNSS PWV的偏差和RMSE分别为0.17 mm和0.35 mm，且两者具较好的相关性和一致性.由此表明，ERA5地表温压产品可应用于桂林地区GNSS水汽反演，这些研究结果可为桂林地区的GNSS水汽反演及数据源的选用提供重要的参考依据.     

云南地基GPS观测大气可降水量变化特征

利用2007年云南地基GPS站点观测资料,分析GPS反演的大气可降水量(PWV)变化特征,并用探空、实际降水量资料和GPS反演结果进行比较。结果表明:GPS/PWV能反映云南降水的季节变化特征,海拔较低的测站普遍比同期海拔较高的测站测得的GPS/PWV值高;GPS/PWV值与探空得到的大气水汽总量随时间演变趋势基本一致,其相关系数均达0.89;GPS/PWV变化周期和实际降水发生的周期基本相同,降水大多为GPS/PWV值连续增加达到峰值(或从峰值开始下降)后开始;GPS/PWV上升幅度较大或位于高位可作为连续性强降水过程出现的预报指标,但使用GPS/PWV峰值作预报指标时,还应考虑季节因素。     

利用地基北斗站反演大气水汽总量的精度检验

采用上海市气象局建立的北斗气象站的2014年观测数据和我国自主研发的精密导航数据处理软件PANDA (position and navigation data analysist) 实现了基于北斗数据的大气水汽总量 (precipitable water vapor,PWV) 反演,并将利用北斗卫星信号解算的大气水汽总量 (W_BD) 结果与目前较为成熟的GPS卫星反演结果 (W_GPS) 和无线电探空反演结果 (W_Radio) 进行对比,研究表明:反演的W_BD与W_GPS的均方根误差均低于3.5 mm,反演的W_BD与W_Radio的均方根误差为3.6 mm,两种对比方式的相关系数均在0.95以上,反演方法以及地理位置的差异对于反演结果有一定影响;反演的W_BD能够很好地反映出大气中水汽的变化特征,对于气象短时临近预报、气候分析有指示作用。     

数值模式预报延迟量与GPS测量的比较研究

利用长江三角洲地区GPS观测网中4个站点的资料，研究了2002年梅雨期静力延迟和湿延迟的时、空变化，以上海宝山站为例用探空资料对GPS资料进行了验证，并将MM5模式24h的预报结果与GPS测量资料进行比较，以检验MM5模式对静力延迟和湿延迟的预报能力，在此基础上提出了订正模式预报静力延迟误差的方法。研究结果表明：GPS是一种与探空观测有可比性的大气静力延迟和湿延迟测量手段。GPS站点测得的静力延迟和湿延迟呈现出不同特征，静力延迟有明显的周期变化，湿延迟的变化与天气过程相关联。模式初始场对静力延迟有很好的描述能力，而对湿延迟却存在低估现象。MM5模式24h滚动预报基本能反映出静力延迟的日变化趋势，而模式在预报静力延迟的逐时次变化方面还存在一定的偏差。模式对湿延迟变化有一定的预报能力，但也存在低估现象。模式对总延迟的预报总体偏低，偏差主要集中在-3cm附近。可以通过订正地面气压的方法减小模式对静力延迟的预报误差。     

GPS/PWV资料同化在强降水过程中的定量作用评估

基于WRF模式及其三维变分同化系统,分别采用江苏GPS观测网的大气可降水量资料及区域内探空和地面气象站资料,对2011年8月25日江苏苏南地区的一次区域性暴雨和2008年苏皖地区的一次局地特大暴雨过程进行了同化试验,用以比较分析GPS/PWV、探空和地面观测资料同化对强降水预报的定量作用。结果表明:探空和地面资料的同化通过对模式动力和热力场的影响,在强降水中心附近形成了强烈的辐合上升运动和热力不稳定条件,直接改进了强降水中心分布结构和强度特征,对数值模拟的成功与否起着决定性作用。而GPS/PWV在探空和地面资料同化的基础上,将使水汽条件得到增强且更有组织性,无论对降水中心强度还是位置都具有较为显著的改进作用。     

IGRA全球站点探空数据集介绍及其资料质量状况分析

美国国家气候数据中心（NCDC）研制的全球站点无线电探空资料数据集（The Integrated Global Radiosonde Archive,简称IGRA）包含全球1538个探空站点从1938年至今的多层次、多要素逐日探空资料。该数据是目前全球范围内时空密度最高、资料最为完整的一套高质量探空站点资料数据集。对该数据集概况及探空资料的质量控制过程进行了概述,并对该资料存在的质量问题进行分析。结果表明：该数据集虽然进行了质量控制,但是该数据集中部分资料的质量仍然存在问题,如该套资料中1973年之前的中国区域探空温度资料存在较大的问题,与国家气象信息中心整理的中国区域探空资料相比,两套资料中各站在1973年以前的不一致率都达到了70 %以上。另外,该套资料中风速资料在1975-1978年以及1981年错误较多,这5 a中出现的错误资料占所有错误资料的50.74 %。而这5 a中错误资料的分布也比较集中,风速出错站点主要分布在美国南部和阿拉斯加、非洲中部和南部以及俄罗斯。全球大部分探空台站通过全球通讯系统进行传输资料集中在1980年以后,因此本数据集中1980年以后的探空资料质量和完整性相对较高。结论和分析为用户使用该数据集提供了很好的参考依据。     

1960~2000年黄河流域太阳总辐射气候变化规律研究

运用数据集群技术，在黄河流域建立了用日照百分率拟合太阳总辐射的不同时空尺度估算式，对比分析了不同数据集群下总辐射的拟合精度。利用逆距离加权插值法，将获得的黄河流域及其周边35个日射站1~12月总辐射拟合的经验系数进行空间内插，获得了黄河流域1~12月总辐射拟合经验系数的空间分布。结合黄河流域及其周边164个常规气象站日照百分率观测资料，对黄河流域1960~2000年太阳总辐射进行了计算，分析了其气候变化的时空分布规律。结果表明:黄河流域近年来太阳总辐射呈下降趋势，在季节上主要表现在夏季和冬季。     

Evaluating regional atmospheric water vapour estimates derived from GPS and short‐range forecasts of the Canadian Global Environmental Multiscale model in southern Alberta

Abstract

Integrated atmospheric moisture has been derived from a network of Global Positioning System (GPS) receivers established in southern Alberta. GPS receivers and post‐processing techniques provide the ability to estimate integrated precipitable water vapour (PWV) at temporal and spatial scales not usually available using conventional observational techniques and without costly expendables. GPS‐derived PWV was evaluated during the Alberta GPS Atmospheric Moisture Evaluation (A‐GAME) using nearby radiosonde observations from the Airdrie, Olds‐Didsbury and Sundre airports during field campaigns in the summers of 2003 and 2004. For the 2004 A‐GAME period, the regional (15 km) Global Environmental Multiscale model (GEM)‐modelled PWV was compared to the GPS derived PWV using a distance weighting approach. GEM model performance was assessed with regards to prognosis time (from 0 to 9 hours), grid cell elevation, location and the presence of storms in the study region. The results show that there is good agreement between radiosonde‐derived PWV and PWV derived from nearby GPS sites with correlations (r²) ranging from 0.76 to 0.84; the GPS‐derived PWV showed a small dry bias averaging 0.6 mm. When compared to GPS‐derived PWV, GEM model performance was found to be favourable out to the hour‐3 prognosis with an overall correlation (r²) of 0.63. Performance decreased with increasing prognosis time and as a result of the presence of storm activity in the study region but did not decrease with increasing grid cell elevation.     

地面有效辐射气候学模型评估和参数优化

基于中国19个辐射站1993-2012年的地面辐射平衡资料和气象资料,分析评估了布朗特法、彭曼法、别尔良德法、FAO24法、FAO56-PM法、邓根云法和童宏良法7种参数化方案计算中国地面有效辐射的适用性;并以均方根误差最小为目标函数,利用步长加速法和多元回归法迭代求解最优参数,建立适合于中国的最优参数化逐日有效辐射估算方法。结果表明:参与评估的7种方案都不同程度低估了中国的有效辐射;从全中国总体误差水平看,童宏良法的平均绝对百分比误差和均方根误差小于其他6种方案,分别为27.0%和24.5 W/m2,估算效果较好;其次是彭曼法和邓根云法;FAO56-PM法精度较低,不适用于中国的有效辐射估算。针对单站来说,邓根云法在东部平原地区的精度最高,童宏良法由于考虑了海拔高度的订正,适用于西部高原地区。相关分析表明水汽压是影响有效辐射估算误差的最关键因素,因此根据水汽压的地理分布规律,分东部区和西部区建立分区方案。基于观测资料建立的全中国方案和分区方案的均方根误差分别为20.8和21.4 W/m2,精度均高于已有参与评估的7种方案;而且在绝大多数站点,分区方案的误差小于全中国方案,所以划分东部区和西部区进行有效辐射模型参数化很有必要。同时发现,分区方案在西部区明显优于邓根云法,在东部区明显优于童宏良法,因此推荐其作为中国有效辐射的计算方法。      

Long-term variations of surface air temperature at major industrial cities of India

Local temperature is one of the major climatic elements to record the changes in the atmospheric environment brought about by industrialization and urbanization. Long-term variations of seasonal and annual surface air temperature at six major industrial cities of India (Calcutta, Bombay, Madras, Bangalore, Pune and Delhi) have been studied, using data for the past 86 to 112 yr. Comparative analysis of the temperature data of six nonindustrial stations has also been done. The long-term change in the temperature has been evaluated by linear trend. Calcutta, Bombay and Bangalore have shown significant warning trend, while Delhi has shown a cooling trend. The trends at Madras and Pune are not significant. The nonindustrial stations did not show significant trends. In general, there was either a cooling tendency or cessation of warming, after the late 1950s at most of the industrial cities.     

云南地基GPS水汽解算方案及精度检验

在分析GAMIT水汽解算方案特点的基础上,利用云南6站地基GPS水汽探测资料,进行了不同解算方案计算结果的分析,并与GPS探空资料的PWV值进行比较,结果显示:不同解算方案对云南地基GPS水汽反演结果有显著的影响;在进行本地化后,GAMIT软件包对云南GPS水汽反演的精度有显著提高;不同季节的云南地基GPS反演水汽值与高精度探空水汽值的数值和变化趋势极为一致,二者的均方根差小于2 mm,说明该解算方案的地基GPS水汽反演结果可用.     

GNSS/PWV与风云四号A星GIIRS水汽廓线融合应用研究

中国新一代地球静止气象卫星风云四号A星（FY-4A）搭载的干涉式大气垂直探测仪（Geostationary Interferometric Infrared Sounder, GIIRS）以红外高光谱干涉分光方式探测三维大气温湿结构,取得了在静止轨道上探测大气的突破性进展。地基全球导航卫星系统（Global Navigation Satellite System,GNSS）是一种连续监测大气可降水量（Precipitable Water Vapor,PWV）的有效手段,基于2018年6—8月中国地基GNSS站监测的PWV和FY-4A/GIIRS水汽廓线的业务产品以及常规无线电探空资料,开展GNSS/PWV与FY-4A/GIIRS水汽廓线快速融合应用,以提高卫星资料反演大气水汽廓线的精度。结果表明:与常规无线电探空相比,FY-4A/GIIRS水汽廓线产品在大气低层均方根误差（Root Mean Square Error,RMSE）为4.5 g/kg,700 hPa为2.4 g/kg,500 hPa以上因水汽含量较低RSME小于1.5 g/kg。GNSS/PWV与FY-4A/GIIRS水汽廓线融合后,FY-4A/GIIRS水汽廓线误差整层RMSE减小20%,从近地层到600 hPa RMSE平均减小20%—25%,尤其是850—700 hPa改善最明显,极大改善了卫星水汽反演资料的可用性。对一次多系统影响的暴雨天气过程应用分析表明,GNSS/PWV和FY-4A/GIIRS融合产品可获得高时、空密度的大气水汽廓线,对强降水的临近预报有非常重要的支撑作用。      

Simulation and Regionalization of Daily Global Solar Radiation: A Case Study in Quebec,Canada

Global solar radiation (GSR) is essential for agricultural and plant growth modelling, air and water heating analyses, and solar electric power systems. However, GSR gauging stations are scarce compared with stations for monitoring common meteorological variables such as air temperature and relative humidity. In this study, one power function, three linear regression, and three non-linear models based on an artificial neural network (ANN) are developed to extend short records of daily GSR for meteorological stations where predictors (i.e., temperature and/or relative humidity) are available. The seven models are then applied to 19 meteorological stations located across the province of Quebec (Canada). On average, the root-mean-square errors (RMSEs) for ANN-based models are 0.33–0.54?MJ?m^?2?d^?1 smaller than those for the power function and linear regression models for the same input variables, indicating that the non-linear ANN-based models are more efficient in simulating daily GSR. Regionalization potential of the seven models is also evaluated for ungauged stations where predictors are available. The power function and the three linear regression models are tested by interpolating spatially correlated at-site coefficients using universal kriging or by applying a leave-one-out calibration procedure for spatially uncorrelated at-site coefficients. Regional ANN-based models are also developed by training the model based on the leave-one-out procedure. The RMSEs for regional ANN models are 0.08–0.46?MJ?m^?2?d^?1 smaller than for other models using the same input conditions. However, the regional ANN-based models are more sensitive to new station input values compared with the other models. Maps of interpolated coefficients and regional equations of the power function and the linear regression models are provided for direct application to the study area.     

中国区域性大气加权平均温度线性模型精度评估

以欧洲中期天气预报中心的再分析资料ERA5为参考数据,评估由探空数据建立的中国区域88个单站大气加权平均温度（T_m）与地表气温（T_s）线性关系模型的精度.各站T_m-T_s线性模型计算的T_m（计算值）与ERA5 气压层数据积分所得的T_m（参考值）间偏差均方根值（RMSE）为1.8~5.5 K.不同站模型计算值与参考值间存在-1.22~4.54 K 的系统性偏差,且绝大多数测站（82个站）系统性偏差为正值,即模型计算值总体上大于参考值.补偿各站系统性偏差后,模型计算值与参考值间RMSE降为1.5~3.5 K.与使用中国区域统一模型相比,使用单站模型平均能提高0.6 K的T_m计算精度,尤其在中国西部、西北和内蒙区域,精度提高可达1~3.9 K.对所有测站模型计算值和参考值间偏差时序进行分析,发现超过半数测站的偏差存在明显季节性变化.     

3种方法对中国探空温度资料均一性检验对比分析

利用美国国家气候数据研究中心开发的Pairwise均一性检验订正方法,加拿大环境研究中心开发的PMFT方法和PMT方法,结合20世纪再分析资料和详细的元数据信息对1951-2008年中国123个探空台站的温度资料进行均一性检验和订正。结果表明：3种方法对中国探空温度资料订正幅度和订正前后的趋势变化存在差异。产生差异的主要原因,一方面是由于Pairwise方法设计的问题,导致对中国大部分台站同时改变辐射误差订正方法和系统升级造成的断点检测能力较弱;另一方面是由于PMFT方法未使用参考序列,不能去除气候变化本身的变化趋势,导致PMFT方法在断点判断上存在漏判的现象。而PMT方法结合20世纪再分析资料的客观判断方法较适合中国探空资料的均一性检验和订正,其订正结果表明,利用该方法得到的中国探空温度数据的统计学特征在低层小于全球尺度探空温度的订正结果,而在高层基本和全球探空数据的订正尺度相当。两个间断点和0.0-0.2℃的订正量所占的比例最大。从全国的来看,对流层为增温趋势,这种增温趋势随高度的增加逐渐减弱,至对流层顶100 hPa转为弱的降温趋势。     

An Operational Precipitable Water Vapor Retrieval Algorithm for Fengyun-2F/VISSR Using a Modified Three-Band Physical Split-Window Method

The Visible and Infrared Spin-Scan Radiometer(VISSR) onboard the Fengyun-2(FY-2) satellite can provide valuable thermal infrared observations to help create a precipitable water vapor(PWV) product with high spatial and temporal resolutions. The current FY-2/VISSR PWV product in operation is produced by using a traditional two-band physical split-window(PSW) method, which produces low quality results under dry atmospheric conditions. Based on the sensitivity characteristics of FY-2 F/VISSR water vapor channel and two split-window channels to atmospheric water vapor, this study developed a new, robust operational PWV retrieval algorithm for FY-2 F to improve the operational precision of the current PWV product. The algorithm uses a modified three-band PSW method, which adds a scale for the water vapor channel in the improved three-band PSW method. Integrated PWV products from the radiosonde data in 2016 are used here to validate the precision of the PWV retrieved by the modified three-band and traditional two-band PSW methods. The mean bias, root mean square error(RMSE), and correlation coefficient of the PWV retrieved by the modified three-band PSW method are 0.28 mm, 4.53 mm, and 0.969, respectively. The accuracy is much better than the PWV retrieved by the two-band method, whose mean bias, RMSE, and correlation coefficient are 12.67 mm, 29.35 mm, and 0.23. Especially, in mid-or high-latitude regions, the RMSE of the PWV is improved from 10 to 2 mm by changing the inversion in the two-band method to the modified three-band PSW method. Furthermore, the modified three-band PSW results show a better consistency with the radiosonde PWV at any zonal belt and season than the two-band PSW results. This new algorithm could significantly improve the quality of the current FY-2 F/VISSR PWV product, especially at sites where the actual PWV are lower than 15 mm.     

区域GPS网对流层延迟直接推算可降水量研究

针对武汉地区GPS气象网资料进行GPS对流层延迟直接推算可降水量的研究。推导了对流层延迟直接推算可降水量的模型,并对模型结果进行检验。从武汉东湖站GPS对流层延迟与无线电探空可降水量的比较中可知,两者具有很好的相关性,相关系数达到0.931;由该站对流层延迟转换的可降水量与无线电探空可降水量的比较可得,均方根为4.45 mm,相关系数为0.905,对流层延迟转换的可降水量与GPS可降水量的均方根为2.23 mm,相关系数为0.988。说明在没有气象数据的地区,对流层延迟直接推算的可降水量可以作为气象短期预报的参考依据。