我国逐日降水量格点化方法

国家气象信息中心(NMIC)和美国大气海洋局气候预测中心合作开发了"中国逐日格点降水量实时分析系统(V1.0)",并已在NMIC投入业务试运行。该系统基于我国2419个国家级地面气象站日降水量观测(08:00—08:00,北京时)数据,采用"基于气候背景场"的最优插值方法,实时生成空间分辨率为0.5°×0.5°的格点化日降水量资料。通过对汛期典型区域和单站降水过程的对比分析表明:该格点化产品的精度较高,能准确捕捉并再现每一次降水过程。误差分析表明:约91%的数据绝对误差小于1.0 mm/d。该产品在定量分析天气实况、检验天气气候模式精度、检验卫星产品精度等方面有应用前景。     

多普勒天气雷达地物回波特征及其识别方法改进

非气象因子会在雷达探测时对雷达资料造成污染,并导致雷达数据的质量问题,在雷达数据应用之前必须对被污染的距离库进行识别和处理。该文在现有基于模糊逻辑识别地物回波工作的基础上,发展适合于我国CINRAD/SA的地物回波识别方法,采用北京和天津雷达2005,2006年夏季部分时段体扫资料,同时对反射率因子和径向速度以及速度谱宽进行处理,得到不同回波的各种特征, 并对各种回波特征进行分析; 考虑到隶属函数的确定是地物识别准确率的关键, 运用CSI (critical success index)评判标准确定了模糊逻辑超折射地物回波识别的最佳线性梯形隶属函数;通过识别效果分析说明该方法在识别超折射地物回波中的作用。结果表明:运用改进后的模糊逻辑法可以更好地识别地物回波, 特别是那些超折射地物回波; 与原方法相比, 改进后的方法有效减少了对降水回波的误判。     

The microscale cooling effects of water sensitive urban design and irrigation in a suburban environment

There is still considerable uncertainty about precipitation at high elevation in mountain terrain due to the relatively few in situ measurements available and to the particular variability of the parameter. In this study, several spatialization techniques were tested, some for climatological time scale and others for daily fields, for precipitation over the western Alps for the period of 1990–2012. The study domain and period were chosen for the quality of available in situ observations and density of the network. First, a weather-type classification was established with a technique based on canonical correlation analysis combining large- and regional-scale data. The spatialization techniques applied for the climatological time scale were adapted from the Aurelhy method which uses elevation and principal components of the topography as predictors. The spatialization techniques applied to daily fields were based on kriging of daily rain gauges and used the climatological fields as predictors. This study aims to validate the advantage of using the climatology of the weather type of the day as predictor for daily fields over a monthly climatology. The climatology of the weather type of the day seems to demonstrate some small improvement.Finally, annual means over the period of 1990–2012 were produced using several methods, including some from accumulation of daily fields and others from the spatialization of in situ station means. Precipitation at high elevations and vertical climatological gradients were particularly scrutinized. Annual means based on sums of daily fields seem to have better performances.This paper only presents results for precipitation but temperature was also analysed.     

Summer heat waves over western Europe 1880–2003, their relationship to large-scale forcings and predictability

We investigate the large-scale forcing and teleconnections between atmospheric circulation (sea level pressure, SLP), sea surface temperatures (SSTs), precipitation and heat wave events over western Europe using a new dataset of 54 daily maximum temperature time series. Forty four of these time series have been homogenised at the daily timescale to ensure that the presence of inhomogeneities has been minimised. The daily data have been used to create a seasonal index of the number of heat waves. Using canonical correlation analysis (CCA), heat waves over western Europe are shown to be related to anomalous high pressure over Scandinavia and central western Europe. Other forcing factors such as Atlantic SSTs and European precipitation, the later as a proxy for soil moisture, a known factor in strengthening land–atmosphere feedback processes, are also important. The strength of the relationship between summer SLP anomalies and heat waves is improved (from 35%) to account for around 46% of its variability when summer Atlantic and Mediterranean SSTs and summer European precipitation anomalies are included as predictors. This indicates that these predictors are not completely collinear rather that they each have some contribution to accounting for summer heat wave variability. However, the simplicity and scale of the statistical analysis masks this complex interaction between variables. There is some useful predictive skill of summer heat waves using multiple lagged predictors. A CCA using preceding winter North Atlantic SSTs and preceding January to May Mediterranean total precipitation results in significant hindcast (1972–2003) Spearman rank correlation skill scores up to 0.55 with an average skill score over the domain equal to 0.28 ± 0.28. In agreement with previous studies focused on mean summer temperature, there appears to be some predictability of heat wave events on the decadal scale from the Atlantic Multidecadal Oscillation (AMO), although the long-term global mean temperature is also well related to western European heat waves. Combining these results with the observed positive trends in summer continental European SLP, North Atlantic SSTs and indications of a decline in European summer precipitation then possibly these long-term changes are also related to increased heat wave occurrence and it is important that the physical processes controlling these changes be more fully understood.     

Wintertime variability of Mediterranean precipitation and its links with large-scale circulation anomalies

A standard principal component analysis has been performed over the Mediterranean and over the larger European region on monthly precipitation anomalies for the winters between 1979 and 1995. The main centres of action of the associated EOFs are very similar for the two regions and the two sets of PCs are highly correlated with each other. Focusing on the Mediterranean region, the same analysis has been performed using 500?hPa geopotential height monthly anomalies taken from the operational NCEP analysis. Comparing the two sets of PCs associated with upper-air and surface data, a strong correlation has been found suggesting the presence of a two-way link between regional precipitation patterns and large-scale circulation anomalies. For both fields, the largest fraction of variance is explained by the North Atlantic Oscillation, while smaller but still substantial fractions are explained by other known patterns of large-scale variability such as the Eastern Atlantic pattern and the Euro-Atlantic blocking. No detectable connection has been found between Mediterranean precipitation patterns and El Niño SST anomalies during winter. With respect to temporal variability, significant trends have been found over most European areas during the winters considered. The associated pattern is characterised by a substantial increase of precipitation over western Scandinavia and a general decrease over southern Europe. This result is confirmed by analysing data from stations located in northern Italy.     

Numerical Simulation of Intense Precipitation Events South of the Alps: Sensitivity to Initial Conditions and Horizontal Resolution

Summary  In this paper we describe the results of several numerical experiments performed with the limited area model LAMBO, based on a 1989 version of the NCEP (National Center for Environmental Prediction) ETA model, operational at ARPA-SMR since 1993. The experiments have been designed to assess the impact of different horizontal resolutions and initial conditions on the quality and detail of the forecast, especially as regards the precipitation field in the case of severe flood events. For initial conditions we developed a mesoscale data assimilation scheme, based on the nudging technique. The scheme makes use of upper air and surface meteorological observations to modify ECMWF (European Centre for Medium Range Weather Forecast) operational analyses, used as first-guess fields, in order to better describe smaller scales features, mainly in the lower troposphere. Three flood cases in the Alpine and Mediterranean regions have been simulated with LAMBO, using a horizontal grid spacing of 15 and 5 km and starting either from ECMWF initialised analysis or from the result of our mesoscale analysis procedure. The results show that increasing the resolution generally improves the forecast, bringing the precipitation peaks in the flooded areas close to the observed values without producing many spurious precipitation patterns. The use of mesoscale analysis produces a more realistic representation of precipitation patterns giving a further improvement to the forecast of precipitation. Furthermore, when simulations are started from mesoscale analysis, some model-simulated thermodynamic indices show greater vertical instability just in the regions where strongest precipitation occurred. Received March 2, 1999/Revised May 30, 1999     

Spatial modelling of air temperature and precipitation for Andorra (Pyrenees) from daily circulation patterns

Based on the daily sea level pressure (SLP) circulation catalogue obtained by Esteban, Martin-Vide and Mases, Int J Climatol 26:1501–1515, (2006) for Western Europe, high-resolution maps of daily maximum and minimum temperature, mean daily precipitation and daily precipitation probability have been obtained for Andorra (Pyrenees). The 20 daily-circulation patterns cover the period 1960–2001 and were generated using new approaches based on the rotated principal component analysis (PCA) and clustering technique. The final maps of Andorra associated with each circulation pattern have been constructed using altitude, latitude, continentality and solar radiation as multiple regression predictors (Ninyerola, Pons and Roure, Int J Climatol 20:1823–1841, 2000). The daily temperature and rainfall series used from Andorran, French and Catalan/Spanish weather stations have been checked for data quality. The results confirm the complexity of the spatial distribution of meteorological phenomena over mountainous areas such as in Andorra, and show the importance of the Mediterranean and Atlantic influence upon the climate of this country of the Pyrenees. On the other hand, different tests have been made showing that the classification results could improve the resulting interpolated climate maps by the use of the circulation-pattern frequencies.     

A strategy for merging objective estimates of global daily precipitation from gauge observations,satellite estimates,and numerical predictions

This paper describes a strategy for merging daily precipitation information from gauge observations, satellite estimates(SEs), and numerical predictions at the global scale. The strategy is designed to remove systemic bias and random error from each individual daily precipitation source to produce a better gridded global daily precipitation product through three steps.First, a cumulative distribution function matching procedure is performed to remove systemic bias over gauge-located land areas. Then, the overall biases in SEs and model predictions(MPs) over ocean areas are corrected using a rescaled strategy based on monthly precipitation. Third, an optimal interpolation(OI)–based merging scheme(referred as the HL-OI scheme)is used to combine unbiased gauge observations, SEs, and MPs to reduce random error from each source and to produce a gauge—satellite–model merged daily precipitation analysis, called BMEP-d(Beijing Climate Center Merged Estimation of Precipitation with daily resolution), with complete global coverage. The BMEP-d data from a four-year period(2011–14) demonstrate the ability of the merging strategy to provide global daily precipitation of substantially improved quality.Benefiting from the advantages of the HL-OI scheme for quantitative error estimates, the better source data can obtain more weights during the merging processes. The BMEP-d data exhibit higher consistency with satellite and gauge source data at middle and low latitudes, and with model source data at high latitudes. Overall, independent validations against GPCP-1DD(GPCP one-degree daily) show that the consistencies between BMEP-d and GPCP-1DD are higher than those of each source dataset in terms of spatial pattern, temporal variability, probability distribution, and statistical precipitation events.     

Circulation patterns, daily precipitation in Portugal and implications for climate change simulated by the second Hadley Centre GCM

Based on principal component analysis (PCA) and a k-means clustering algorithm, daily mean sea level pressure (MSLP) fields over the northeastern Atlantic and Western Europe, simulated by the Hadley Centre's second generation coupled ocean-atmosphere GCM (HADCM2) control run (HADCM2CON), are validated by comparison with the observed daily MSLP fields. It is clear that HADCM2 reproduces daily MSLP fields and its seasonal variability over the region very well, despite suffering from some deficiencies, such as the systematic displacement of the atmospheric centres of action. Four daily circulation patterns, previously identified from the observed daily MSLP fields over the area and well related to daily precipitation in Portugal, were also well classified from the daily MSLP fields simulated by HADCM2. The model can also simulate rather successfully the relationships between the four daily circulation patterns and daily precipitation in southern Portugal. However, compared with observations, daily precipitation intensities simulated by the model are too weak in southern Portugal. Nevertheless, HADCM2 represents a considerable improvement relative to the UKTR experiment. The results described here imply that it is doubtful whether regional precipitation scenarios provided by HADCM2 can be directly applied in impact studies and that a downscaling technique, based on daily circulation patterns, might be successful in reproducing local and regional precipitation characteristics. Moreover, the four circulation patterns can also be clearly identified in the two perturbed experiments, one under greenhouse gases forcing only (HADCM2GHG) and the other under additional forcing of sulphate aerosol (HADCM2SUL), although changes in the frequencies of occurrence of certain circulation patterns are found. Nevertheless, the observed links between regional precipitation in southern Portugal and large-scale atmospheric circulation seem likely to hold in the model's perturbed climate. It is therefore credible to use those links to downscale large-scale atmospheric circulation from GCM simulations to obtain future precipitation scenarios in southern Portugal. Received: 21 August 1998 / Accepted: 28 May 1999     

华南5-6月降水量场统计特征与站网关系

由华南区15、42个测站构成了相互独立的低、高密度两种站网。用1958-2000年5-6月逐日降水资料，分析了两种站网华南前汛期降水量场特征的异同。统计表明：1）两种站网给出的华南区域平均季、月降水量无显著差异，候、日降水量有显著差异；2）两种站网给出的逐日气候及其异常场均存在显著差异。     

基于雷达资料4DVar的低层热动力反演系统及其在北京奥运期间的初步应用分析

在变分多普勒雷达分析系统(VDRAS)中,通过对雷达资料的质量控制和预处理、云尺度模式的暖雨参数化方案、中尺度初猜场的插值分析方法、循环同化的冷启动和热启动中尺度背景场的计算方案等进行改进,实现了VDRAS对低层动力和热力场的分析反演及其在北京奥运期间的实时应用.改进后的VDRAS利用四维变分(4DVar)同化技术和一...     

Impact of the lateral boundary conditions resolution on dynamical downscaling of precipitation in mediterranean spain

Conclusions on the General Circulation Models (GCMs) horizontal and temporal optimum resolution for dynamical downscaling of rainfall in Mediterranean Spain are derived based on the statistical analysis of mesoscale simulations of past events. These events correspond to the 165 heavy rainfall days during 1984–1993, which are simulated with the HIRLAM mesoscale model. The model is nested within the European Centre for Medium-Range Weather Forecasts atmospheric grid analyses. We represent the spectrum of GCMs resolutions currently applied in climate change research by using varying horizontal and temporal resolutions of these analyses. Three sets of simulations are designed using input data with 1°, 2° and 3° horizontal resolutions (available at 6 h intervals), and three additional sets are designed using 1° horizontal resolution with less frequent boundary conditions updated every 12, 24 and 48 h. The quality of the daily rainfall forecasts is verified against rain-gauge observations using correlation and root mean square error analysis as well as Relative Operating Characteristic curves. Spatial distribution of average precipitation fields are also computed and verified against observations. For the whole Mediterranean Spain, model skill is not appreciably improved when using enhanced spatial input data, suggesting that there is no clear benefit in using high resolution data from General Circulation Model for the regional downscaling of precipitation under the conditions tested. However, significant differences are found in verification scores when boundary conditions are interpolated less frequently than 12 h apart. The analysis is particularized for six major rain bearing flow regimes that affect the region, and differences in model performance are found among the flow types, with slightly better forecasts for Atlantic and cold front passage flows. A remarkable spatial variability in forecast quality is found in the domain, with an overall tendency for higher Relative Operating Characteristic scores in the west and north of the region and over highlands, where the two previous flow regimes are quite influential. The findings of this study could be of help for dynamical downscaling design applied to future precipitation scenarios in the region, as well as to better establish confidence intervals on its results.     

逐时云迹风资料同化对暴雨预报的模拟试验

文中基于不同的云迹风同化方案,用GRAPES模式对2005年7月11-12日长江中下游一次暴雨强降水过程进行了云迹风资料同化试验及数值模拟,通过对比分析不同方案所得的分析场及预报场的差异,研究逐时云迹风资料三维变分同化对分析场及暴雨预报的影响.首先,根据连续性原理及双通道各层次云迹风资料的误差分析,分3个步骤对7月11日00:00-12:00 UTC共12个时次的双通道云迹风资料进行了初步的质量控制;然后,把经过质量控制的云迹风资料放入基于GRAPES 3D-VAR三维变分同化方案开发的逐小时循环同化系统中进行同化,将得出的分析场与单一时次未经质量控制的云迹风资料同化得出的分析场进行对比,探讨了逐时云迹风资料同化对数值预报分析场的影响;最后,把同化后的分析场作为初始场,用 GRAPES模式对 2005年7月11-12日长江流域暴雨过程做24 h降水预报试验,分析两个同化方案所模拟得到的预报场的差异.结果表明,经资料的筛选、同经纬度单点通道的选择及资料的稀疏化3个步骤控制后,各层次云迹风资料的误差有明显减小;加入经质量控制的逐时云迹风资料其三维变分同化可以提高分析场中风压场及水汽场的质量;而且在暴雨预报试验中可以相对更准确地预报暴雨落区及雨强.     

Fuzzy rule-based downscaling of precipitation

Summary A fuzzy rule-based methodology for downscaling local hydrological variables from large-scale atmospheric circulation is presented. The method is used to estimate the frequency distribution of daily precipitation conditioned on daily geopotential fields. The task is accomplished in two steps. First, the exceedence probabilities corresponding to selected precipitation thresholds are estimated by fuzzy rules defined between geopotential fields (premises) and exceedence events (response). Then a continuous probability distribution is constructed from the discrete exceedence probabilities and the observed behaviour of precipitation. The methodology is applied to precipitation measured at Essen, a location in the Ruhr catchment, Germany. Ten years of precipitation data (1970–1979) were used for training and another ten years (1980–1989) for validation. The 700 hPa geopotential fields are used to characterise large-scale circulation. The application example demonstrates that this direct downscaling method is able to capture the relationship between premises and the response; namely both the estimated exceedence probabilities and the frequency distribution reproduce the empirical data observed in the validation period.     

多普勒天气雷达1 h降水产品的质量评估

为评估多普勒天气雷达1 h降水量(OHP)产品在降水预报中的质量, 从南通新一代多普勒天气雷达的体扫资料及其覆盖范围内太湖流域34个雨量站的同期逐小时降水资料中选取了19个降水个例。以空报率、漏报率、准确预报临界指数、平均误差、平均绝对误差和均方根误差等, 为评估降水预报质量的指标。根据雨量站实测日降水量(P)和1 h降水量(P₁)大小分级, 对不同情况下的OHP进行质量评估。结果表明:(1)OHP的空报率和漏报率较高, 导致整体的准确预报临界指数偏低。(2)与不同的P和P₁量级相对应, 不同的OHP预报质量有明显差异。除P<10 mm时和P₁<1 mm时稍有高估外, 其他量级都存在低估, 且随着P的增加, 空报率和漏报率都减小。(3)对2007年9月18日的降水个例分析表明:OHP与实测的落区基本一致, 但降水强度上有偏差。不同站点的雨强偏差不同, 距南通雷达站100 km左右的OHP误差较小。     

多普勒天气雷达资料同化对江淮暴雨模拟的影响

利用GSI同化系统 (Gridpoint Statistical Interpolation System) 对我国多普勒天气雷达资料进行直接循环同化分析,并采用WRF-ARW 3.5.1中尺度模式对2013年我国夏季江淮流域典型暴雨过程进行模拟试验。结果表明:经过质量控制的雷达径向风、反射率因子资料经GSI同化系统同化后,可形成合理的分析增量。仅同化径向风,模拟的风场与实况更接近,模拟的降雨落区与观测雨带位置更加接近。仅同化反射率因子,对水平风场的直接调整比较小,通过水凝物含量调整,对水平和垂直风场进行调整,对降水的落区影响较小,主要影响模拟的降水强度。同时同化两种资料,能更好地反映风场特征,并改善强降水的落区和强度的模拟。模拟改善最明显是在积分12~36 h时段内,该时段有效雷达资料量较多,表明同化雷达资料对暴雨模拟确实具有正效果。     

Radar observations of snow formation in a warm pre‐frontal snowband

Abstract

Radar reflectivity measurements and sounding data were analyzed to investigate snowfall production in a long‐lasting snowband that formed in advance of a warm surface front moving across Alberta. The sounding data indicated that the band could have been forced by slantwise overturning during the release of moist symmetric instability combined with frontogenesis. The stability analysis presented here is novel in that it includes ice phase thermodynamics, neglected in previous studies of slantwise convection.

Radar reflectivity fields were analyzed to determine the total snow content and the mass outflow rate as factors of time. The peak value of total snow content was 17 kilotons per km of snowband, and the peak mass outflow rate was 10 tons s^‐1 km^‐1. The snowfall rate averaged across the cloud base was about 0.8 cm h^‐1, and the average snow content remained close to 0.2 g m^‐1. The characteristic time (defined as the ratio of total snow content over mass outflow rate) was about 30 minutes, which is approximately the time needed for the growth of snowflakes by aggregation in the observed temperature range. The precipitation efficiency of the snowband, defined as the ratio of snow mass outflow to water vapour inflow was estimated to be 14%. The precipitation production values observed in the Alberta snowband are compared with previous estimates reported for frontal rainbands and Alberta thunderstorms.     

Estimation of Indian summer monsoon rainfall using Kalpana-1 VHRR data and its validation using rain gauge and GPCP data

In the present study, an attempt has been made to estimate and validate the daily and monthly rainfall during the Indian summer monsoon seasons of 2008 and 2009 using INSAT (Indian National Satellite System) Multispectral Rainfall Algorithm (IMSRA) technique utilizing Kalpana-1 very high resolution radiometer (VHRR) measurements. In contrary to infrared (IR), microwave (MW) rain rates are based on measurements that sense precipitation in clouds and do not rely merely on cloud top temperature. Geostationary satellites provide broad coverage and frequent refresh measurements but microwave measurements are accurate but sparse. IMSRA technique is the combination of the infrared and microwave measurements which make use of the best features of both IR- and MW-based rainfall estimates. The development of this algorithm included two major steps: (a) classification of rain-bearing clouds using proper cloud classification scheme utilizing Kalpana-1 IR and water vapor (WV) brightness temperatures (Tb) and (b) collocation of Kalpana-1 IR brightness temperature with Tropical Rainfall Measuring Mission (TRMM)-Precipitation Radar (PR) surface rain rate and establishment of a regression relation between them. In this paper, the capability of IMSRA as an operational algorithm has been tested for the two monsoon seasons 2008 and 2009. For this, IMSRA has been used to estimate daily and monthly rainfall and has been intercompared on daily and monthly scales with TRMM Multisatellite Precipitation Analysis (TMPA)-3B42 V6 product and Global Precipitation Climatology Project (GPCP) rain product during these two monsoon years. The daily and monthly IMSRA rainfall has also been validated against ground-based observations from Automatic Weather Station (AWS) Rain Gauge and Buoy data. The algorithm proved to be in good correlation with AWS data over land up to 0.70 for daily rain estimates except orographic regions like North-East and South-West India and 0.72 for monthly rain estimates. The validation with Buoys gives the reasonable correlation of 0.49 for daily rain estimates and 0.66 for monthly rain estimates over Tropical Indian Ocean.     

近57a南疆帕米尔地区春季降水季节内差异及环流异常分析

基于1961—2017年帕米尔地区3站的日降水资料及NECP/NCAR再分析资料,利用气候诊断分析及多元统计学方法,研究了近57 a来帕米尔地区春季降水特征及季节内差异,并讨论了降水偏多年大气环流异常特征。研究表明:(1)近57 a来,帕米尔地区季节尺度上和月尺度上降水增加明显,旬尺度上有不显著的增加趋势。21世纪初帕米尔地区进入新的多雨期,降水正距平年份强度和频率均有显著增加。(2)500 hPa高度场上3月东欧沿岸脊发展、东欧—西西伯利亚槽加深,4—5月欧洲沿岸槽加深、乌拉尔山脊发展、巴尔喀什湖槽加深是帕米尔地区降水异常偏多的关键系统及指标。(3)降水异常年的高空急流较常年强度更强、位置更偏东,在旬尺度上是一个逐渐减弱西退的过程;其低层风场距平分为西南风距平辐合型、西南风和东北风距平辐合型和气旋式环流距平辐合型,均有利于低层辐合。(4)降水异常偏多年水汽输送以偏西路径为主,同时还有偏东、西南和偏南路径,水汽辐合强度较常年更强,更有利于降水的产生。