ECMWF全风速场集合预报结果的偏差订正与预报不一致性分析

采用卡尔曼滤波类型自适应误差订正法和滑动自适应权重法,对2012年夏季ECMWF 10 m全风速场集合预报结果进行偏差订正,对订正前后的预报结果进行评估,并通过Jumpiness指数对预报结果订正前后的预报不一致性特征进行分析。结果表明,卡尔曼滤波类型自适应误差订正法能有效降低集合预报的均方根误差,且当起报时刻为00时对中低纬度地区的订正效果更显著,当起报时刻为12时对中高纬度地区的订正效果更明显;卡尔曼滤波类型自适应误差订正法能有效改善Talagrand的U型或L型分布;由均方根误差分析结果知道,ECMWF 10 m全风速场集合预报本身存在较大的预报不一致性,经过卡尔曼滤波类型自适应误差订正后,集合预报的预报不一致性明显降低,偏差订正可有效改善集合预报的预报不一致性,且随着预报时效的延长,卡尔曼滤波误差法对预报不一致性的改善效果更加明显;从预报不一致性的发生次数特征来看,单点跳跃出现的次数最多,异号三点跳跃的次数最少;经过卡尔曼滤波类型自适应误差订正后,单点跳跃、异号两点跳跃、异号三点跳跃次数都有所下降。     

Prediction of the Caspian Sea level using ECMWF seasonal forecasts and reanalysis

The hydrological budget of the Caspian Sea (CS) is investigated using the European Centre for Medium-Range Weather Forecasts interim reanalysis (ERAi) and seasonal forecast (FCST) data with the aim of predicting the Caspian Sea Level (CSL) some months ahead. Precipitation and evaporation are used. After precipitation events over the Volga River, the discharge (Volga River discharge (VRD)) follows with delays, which are parameterized. The components of the water budget from ERAi and FCSTs are integrated to obtain time series of the CSL. Observations of the CSL and the VRD are used for comparison and tuning. The quality of ERAi data is sufficiently good to calculate the time variability of the CSL with a satisfactory accuracy. Already the storage of water within the Volga Basin allows forecasts of the CSL a few months ahead, and using the FCSTs of precipitation improves the CSL forecasts. The evaporation in the seasonal forecasts is deficient due to unrealistic sea surface temperatures over the CS. Impacts of different water budget terms on the CSL variability are shown by a variety of validation tools. The importance of precipitation anomalies over the catchment of the Volga River is confirmed, but also impacts from the two southern rivers (Sefidrud and Kura River) and the evaporation over the CS become obvious for some periods. When pushing the FCSTs beyond the limits of the seasonal FCSTs to 1 year, considerable forecast skill can still be found. Validating only FCSTs by the present approach, which show the same trend as one based on a statistical method, significantly enhances the skill scores.     

Progress towards better weather forecasts for city dwellers: from short range to climate change

Summary The current resolution of operational weather forecast model is not sufficient in general to explicitly resolve even the major cities of the World. As a consequence, urban areas have traditionally been neglected in such models. The introduction of tiled land surface models has enabled sub-gridscale landuse to be modelled, and hence has provided the opportunity to model cities within weather forecast models. However, to date there has been little effort made within the operational weather forecast community. At present there is only one operational centre that explicitly resolves urban areas. This centre includes a simple urban scheme within its mesoscale and global models, which has been shown to have a positive impact on the forecast. However, with the recent developments within urban meteorology there are now a variety of urban schemes, which vary in their complexity and parameter requirements, that would be suitable for operational weather forecast models. So it is likely that more operational models, and in particular mesoscale models, will include urban areas in the near future. With the majority of the World‘s population living in cities, the resilience of these cities to the impacts of climate change is also becoming of increasing interest. This means that urban areas will have to be included within climate change simulations, as well as weather forecast simulations, in the future. At present, only one climate change model has included a parametrisation for urban areas. However, this is likely to increase if work in this area grows rapidly.     

Multimodel based superensemble forecasts for short and medium range NWP over various regions of Africa

Summary This study examines the predictability of weather over several regions in Africa using a multimodel superensemble technique developed at the Florida State University, which is an objective means of combining daily forecasts from multilevel global models. It is referred to as FSUSE and up to 7 different models are used to construct the superensemble. The benchmark reanalysis fields used are the precipitation data sets from CMORPH and all other global fields from ECMWF daily operational analysis. The FSUSE works by using multiple linear regression to derive weights from a comparison of each member model forecast to the benchmark analysis during a training period of the most recent 120 days, and these weights are passed to the forecast phase. This procedure removes the bias of each model and allows for an optimal linear combination of the individual model forecasts by taking account of the relative skill of each model to give a consensus forecast that is superior to the ensemble mean and all the members. Results show that bad models and poor analysis fields used during the training phase degrade the skill of the FSUSE. In the forecasts of rainfall events over all regions of Africa, the FSUSE root-mean-square (R M S) error, equitable threat skill score (E T S), and bias on the daily forecasts of rainfall were invariably superior to the best member model. The skills deteriorate as the forecast lead time in days increases, with the degradation being most significant beyond day 3. In all cases, the bias score of the FSUSE was approximately 1, while the anomaly correlation scores were to the order of 0.9. These scores indicate the robustness of the FSUSE forecasts. Over East Africa, the FSUSE forecasts were consistent with the spatial-temporal pattern of the Intertropical Convergence Zone (ITCZ), the main rain bearing synoptic mechanism across tropical Africa. Thus, in addition to superior forecasts, the use of FSUSE based data sets may provide a better understanding of the dynamical processes within the ITCZ over the region. These results could be further improved if the daily series of operational analysis had included gauge data and if the resolution were higher. It is hardly possible to get uniformly consistent and continuous daily observations over these diverse regions of Africa. However, given the availability of the satellite based estimates of daily rainfall, such as CMORPH and global analysis that are exchanged very fast nowadays, the FSUSE scheme for numerical weather predictions (N W P) provides useful medium range weather forecasts in real-time.     

基于ECMWF模式预报的台风降水地形订正方法

为了解决复杂地形条件下ECMWF模式预报的台风降水较实况显著小的问题,对Smith 1979年提出的地形降水方程进行改进,提出以饱和湿层高度作为方程积分上限,针对不同高度地形设定不同的降水效率;以无量纲湿弗劳德数大于1作为有、无地形降水的判据;利用ECWMF细网格预报场,通过迎风坡地形降水估算方程来订正模式预报的台风降水。用该地形降水订正方法对1617号台风“鲇鱼”的降水进行了订正预报。结果表明,虽然在一些小尺度地形区域会产生明显的空报,但是对于大尺度地形区域的强降水有显著的订正效果。对1513、1521和1614台风的订正结果进一步表明,该地形降水订正方法对改进台风极端降水预报效果显著。需要指出的是,采用的地形降水订正方法仅考虑了稳定条件下的地形降水,对于其他情形下的地形降水订正方法尚待进一步的研究。      

Asian summer monsoon prediction in ECMWF System 4 and NCEP CFSv2 retrospective seasonal forecasts

The seasonal prediction skill of the Asian summer monsoon is assessed using retrospective predictions (1982–2009) from the ECMWF System 4 (SYS4) and NCEP CFS version 2 (CFSv2) seasonal prediction systems. In both SYS4 and CFSv2, a cold bias of sea-surface temperature (SST) is found over the equatorial Pacific, North Atlantic, Indian Oceans and over a broad region in the Southern Hemisphere relative to observations. In contrast, a warm bias is found over the northern part of North Pacific and North Atlantic. Excessive precipitation is found along the ITCZ, equatorial Atlantic, equatorial Indian Ocean and the maritime continent. The southwest monsoon flow and the Somali Jet are stronger in SYS4, while the south-easterly trade winds over the tropical Indian Ocean, the Somali Jet and the subtropical northwestern Pacific high are weaker in CFSv2 relative to the reanalysis. In both systems, the prediction of SST, precipitation and low-level zonal wind has greatest skill in the tropical belt, especially over the central and eastern Pacific where the influence of El Nino-Southern Oscillation (ENSO) is dominant. Both modeling systems capture the global monsoon and the large-scale monsoon wind variability well, while at the same time performing poorly in simulating monsoon precipitation. The Asian monsoon prediction skill increases with the ENSO amplitude, although the models simulate an overly strong impact of ENSO on the monsoon. Overall, the monsoon predictive skill is lower than the ENSO skill in both modeling systems but both systems show greater predictive skill compared to persistence.     

ECMWF业务预报模式地面气温预报的不一致性特征研究

为探讨ECMWF业务预报模式(以下简称ECMWF)的地面气温预报不一致性问题,本文利用2015年12月1日—2016年11月30日业务预报中常用的地面气温预报数据,研究ECMWF地面气温预报产品在不同季节里的不一致性指数分布及变化特征。结果表明:各个季节不一致性指数有不同的特点,冬季不一致性指数最大,大值区主要分布在除华南和青藏高原外的大部分区域;而夏季不一致性指数最小,在青藏高原地区不一致性指数相对较大;春、秋两季不一致性指数大小均处于冬、夏季之间。此外,研究还发现冬季地面气温预报不一致性指数单日变化较大,而夏季较小。夏季不同起报时间的地面气温预报比较稳定。     

An application of chaotic theory in long range forecasts

Summary In this paper we combine chaotic theory with statistics to present three forecast models: a model of equal distance ind-dimensional phase space, a mode regression model of twelve units and a model of neighborhood model regression. Many experiments show that all these models can generate accurate forcasts.     

Application of multimodel ensemble techniques for real time district level rainfall forecasts in short range time scale over Indian region

During the summer monsoon (1 June to 30 September) 2007, real-time district level rainfall forecasts in short-range time scale were generated for Indian region applying multimodel ensemble technique. The pre-assigned grid point weights on the basis of correlation coefficients (CC) between the observed values and forecast values are determined for each constituent model at the resolution of 0.5° × 0.5° utilizing two seasons datasets (1 June to 30 September, 2005 and 2006), and the multimodel ensemble forecasts (day 1 and day 2 forecasts) are generated at the same resolution on a real-time basis. The ensemble forecast fields are then used to prepare forecasts for each district taking the average value of all grid points falling in a particular district. In this paper we examined the performance skill of the multimodel ensemble-based real-time district level short-range forecast of rainfall. It has clearly emerged from the results that the multimodel ensemble technique reported in this study is superior to each ensemble member. District wise performance of the ensemble rainfall forecast reveals that the technique, in general, is capable of providing reasonably good forecast skill over most districts of the country, particularly over the districts where the monsoon systems are dominant. Though the procedure shows appreciable skill to predict occurrence or non-occurrence of rainfall at the district level, it always underestimates rainfall amount, particularly in heavy rainfall events. Possible reasons of this failure may be due to model bias and poor data assimilation procedure.     

The role of horizontal resolution in simulating drivers of the global hydrological cycle

The role of atmospheric general circulation model (AGCM) horizontal resolution in representing the global energy budget and hydrological cycle is assessed, with the aim of improving the understanding of model uncertainties in simulating the hydrological cycle. We use two AGCMs from the UK Met Office Hadley Centre: HadGEM1-A at resolutions ranging from 270 to 60 km, and HadGEM3-A ranging from 135 to 25 km. The models exhibit a stable hydrological cycle, although too intense compared to reanalyses and observations. This over-intensity is explained by excess surface shortwave radiation, a common error in general circulation models (GCMs). This result is insensitive to resolution. However, as resolution is increased, precipitation decreases over the ocean and increases over the land. This is associated with an increase in atmospheric moisture transport from ocean to land, which changes the partitioning of moisture fluxes that contribute to precipitation over land from less local to more non-local moisture sources. The results start to converge at 60-km resolution, which underlines the excessive reliance of the mean hydrological cycle on physical parametrization (local unresolved processes) versus model dynamics (large-scale resolved processes) in coarser HadGEM1 and HadGEM3 GCMs. This finding may be valid for other GCMs, showing the necessity to analyze other chains of GCMs that may become available in the future with such a range of horizontal resolutions. Our finding supports the hypothesis that heterogeneity in model parametrization is one of the underlying causes of model disagreement in the Coupled Model Intercomparison Project (CMIP) exercises.     

Prediction skill of the Madden and Julian Oscillation in dynamical extended range forecasts

The Madden and Julian Oscillation (MJO) is the most prominent mode of intraseasonal variations in the tropical region. It plays an important role in climate variability and has a significant influence on medium-to-extended ranges weather forecasting in the tropics. This study examines the forecast skill of the oscillation in a set of recent dynamical extended range forecasts (DERF) experiments performed by the National Centers for Environmental Prediction (NCEP). The present DERF experiments were done with the reanalysis version of the medium range forecast (MRF) model and include 50-day forecasts, initialized once-a-day (0Z) with reanalyses fields, for the period between 1 January, 1985, and 31 December, 1989. The MRF model shows large mean errors in representing intraseasonal variations of the large-scale circulation, especially over the equatorial eastern Pacific Ocean. A diagnostic analysis has considered the different phases of the MJO and the associated forecast skill of the MRF model. Anomaly correlations on the order of 0.3 to 0.4 indicate that skillful forecasts extend out to 5 to 7 days lead-time. Furthermore, the results show a slight increase in the forecast skill for periods when convective anomalies associated with the MJO are intense. By removing the mean errors, the analysis shows systematic errors in the representation of the MJO with weaker than observed upper level zonal circulations. The examination of the climate run of the MRF model shows the existence of an intraseasonal oscillation, although less intense (50–70%) and with faster (nearly twice as fast) eastward propagation than the observed MJO. The results indicate that the MRF model likely has difficulty maintaining the MJO, which impacts its forecast. A discussion of future work to improve the representation of the MJO in dynamical models and assess its prediction is presented. Received: 28 December 1998 / Accepted: 27 September 1999     

江淮流域夏季持续性强降水的低频特征分析

基于1986-2015年湖南逐日降水资料、同期美国气象环境预报中心（NCEP）和美国国家大气研究中心（NCAR）再分析资料,通过分析强低频振荡年的汛期强降水特征和低频环流场演变对强降水的影响,建立了湖南省汛期延伸期强降水过程预报指数。结果表明：（1）汛期33%的强降水过程均发生在具有显著30~60 d低频振荡的年份中,且大多位于低频降水峰值阶段。（2）通过对强低频振荡年进行合成发现,在活跃位相,南亚高压偏强偏东,副热带高压偏西偏强,这种环流配置导致中国南方大部分地区的高层环流为辐散,底层环流为辐合,有利于降水的产生。在中断位相,南亚高压呈东西带状分布且其位置偏西、强度偏弱,副热带高压偏东偏弱,使得向湖南地区输送水汽的西南气流减弱,进入降水中断期。（3）基于低频散度场不同位相的变化特征,选取了与低频降水相关的两个关键区,从而建立延伸期预报指数,该指数对低频降水显著年的强降水回报准确率能够达到73%。（4）前期4月黑潮的海温异常（SSTA）可作为湖南省强低频振荡年的预测指标。

     

Effects of heterogeneous vegetation on the surface hydrological cycle

Using the three-layer variable infiltration capacity (VIC-3L) hydrological model and the successive interpolation approach (SIA) of climate factors, the authors studied the effect of different land cover types on the surface hydrological cycle. Daily climate data from 1992 to 2001 and remotely-sensed leaf area index (LAI) are used in the model. The model is applied to the Baohe River basin, a subbasin of the Yangtze River basin, China, with an area of 2500 km2. The vegetation cover types in the Baohe River basin consist mostly of the mixed forest type (85%). Comparison of the modeled results with the observed discharge data suggests that: (1) Daily discharges over the period of 1992–2001 simulated with inputs of remotely-sensed land cover data and LAI data can generally produce observed discharge variations, and the modeled annual total discharge agrees with observations with a mean difference of 1.4%. The use of remote sensing images also makes the modeled spatial distributions of evapotranspiration physically meaningful. (2) The relative computing error (RCE) of the annual average discharge is ?24.8% when the homogeneous broadleaf deciduous forestry cover is assumed for the watershed. The error is 21.8% when a homogeneous cropland cover is assumed and ?14.32% when an REDC (Resource and Environment Database of China) land cover map is used. The error is reduced to 1.4% when a remotely-sensed land cover at 1000-m resolution is used.     

The global energy cycle of stationary and transient atmospheric waves: Results from ECMWF analyses

Summary The role of stationary (monthly mean) and transient (departure from monthly mean) waves within the atmospheric energy cycle is examined using global analyses from the European Centre for Medium Range Weather Forecasts (ECMWF) for the period 1980–1987. Only January and July averages are considered.It is confirmed that planetary stationary waves are basically baroclinic. Their contribution to the globally averaged energy cycle of the atmosphere is comparable to that of the transient waves. In January they contribute about 40% to the baroclinic conversion (CA) from zonal mean to eddy available potential energy. Local values for the northern hemisphere even show a predominant role of the stationary wave conversions over those originating from transient waves. Part of the available potential energy of stationary waves (A _SE) is converted to kinetic energy by warm air rising and cold air sinking. Nonlinear energy conversion, which can be interpreted as destruction of stationary temperature waves by transients, is the second sink forA _SE. The order of magnitude of these two processes is similar.Barotropic nonlinear conversions, though negligible in the global average, reveal large conversion rates between the mean positions of the polar and the subtropical jets. Their orientation is suggestive of a tendency to increase stationary wave kinetic energyK _SE at its local minimum between the jets at the expense of the synoptic scale transients.While all terms of the energy cycle related to stationary waves reveal a predominance of the planetary scale (zonal wave numbers 1–3) transient waves are governed by synoptic scale waves (zonal wave numbers 4–9) only with respect to the baroclinic and barotropic conversions: a significant amount of transient wave energy (50% for the global average ofA _TE) is due to planetary scale waves.With 15 Figures     

Emission scenario dependencies in climate change assessments of the hydrological cycle

Anthropogenic global warming will lead to changes in the global hydrological cycle. The uncertainty in precipitation sensitivity per 1 K of global warming across coupled atmosphere-ocean general circulation models (AOGCMs) has been actively examined. On the other hand, the uncertainty in precipitation sensitivity in different emission scenarios of greenhouse gases (GHGs) and aerosols has received little attention. Here we show a robust emission-scenario dependency (ESD); smaller global precipitation sensitivities occur in higher GHG and aerosol emission scenarios. Although previous studies have applied this ESD to the multi-AOGCM mean, our surprising finding is that current AOGCMs all have the common ESD in the same direction. Different aerosol emissions lead to this ESD. The implications of the ESD of precipitation sensitivity extend far beyond climate analyses. As we show, the ESD potentially propagates into considerable biases in impact assessments of the hydrological cycle via a widely used technique, so-called pattern scaling. Since pattern scaling is essential to conducting parallel analyses across climate, impact, adaptation and mitigation scenarios in the next report from the Intergovernmental Panel on Climate Change, more attention should be paid to the ESD of precipitation sensitivity.     

基于TIGGE资料的地面气温延伸期多模式集成预报

基于TIGGE资料中心提供的CMC、ECMWF、UKMO及NCEP四个集合预报中心2008年7月1日-9月30日北半球中纬度地区地面气温10 ～ 15 d延伸期集合预报产品,首先采用Tala-grand分布及离散度—误差关系评估了单个预报系统的预报性能,然后分别利用多模式集成平均(Ensemble Mean,EMN)、消除偏差集成平均(Bias-Removed Ensemble Mean,BREM)及多模式超级集合(Multi-model Superensemble,SUP)对地面气温进行多模式集成预报试验.由于逐日的延伸期预报准确率相对较低,因此人们更关注延伸期预报对天气过程的预报准确率.对各个集合预报系统的逐日预报资料以及逐日“观测”资料做滑动平均,并对处理后的资料进行多模式集成,最后对超级集合预报的训练期长度进行调试,以获得最佳训练期长度.结果表明,四个集合预报系统的离散度相对于均方根误差都偏小,ECMWF预报效果最好,NCEP次之,UKMO预报效果最差.EMN、BREM及SUP三种多模式集成方法的预报效果均优于单个系统且SUP对预报效果的改善最明显.滑动平均后,预报误差进一步降低,且滑动步长越长,误差越小.对于SUP的训练期,逐日预报和3d滑动平均10～12 d预报最佳训练期长度为75 d;13 ～ 15 d预报最佳训练期长度为35 d;5 d及7d滑动平均其训练期长度在各个时效均以35 d为宜.     

Forecast skill of the tropical intraseasonal oscillation in the NCEP GFS dynamical extended range forecasts

This study examines the forecast performance of tropical intraseasonal oscillation (ISO) in recent dynamical extended range forecast (DERF) experiments conducted with the National Centers for Environmental Prediction (NCEP) Global Forecasting System (GFS) model. The present study extends earlier work by comparing prediction skill of the northern winter ISO (Madden-Julian Oscillation) between the current and earlier experiments. Prediction skill for the northern summer ISO is also investigated. Since the boreal summer ISO exhibits northward propagation as well as eastward propagation along the equator, forecast skill for both components is computed. For the 5-year period from 1 January, 1998 through 31 December, 2002, 30-day forecasts were made once a day. Compared to the previous DERF experiment, the current model has shown some improvements in forecasting the ISO during winter season so that the skillful forecasts (anomaly correlation>0.6) for upper-level zonal wind anomaly extend from the previous shorter-than 5 days out to 7 days lead-time. A similar level of skill is seen for both northward and eastward propagation components during the summer season as in the winter case. Results also show that forecasts from extreme initial states are more skillful than those from null phases for both seasons, extending the skillful range by 3–6 days. For strong ISO convection phases, the GFS model performs better during the summer season than during the winter season. In summer forecasts, large-scale circulation and convection anomalies exhibit northward propagation during the peak phase. In contrast, the GFS model still has difficulties in sustaining ISO variability during the northern winter as in the previous DERF run. That is, the forecast does not maintain the observed eastward propagating signals associated with large-scale circulation; rather the forecast anomalies appear to be stationary at their initial location and decay with time. The NCEP Coupled Forecast System produces daily operational forecasts and its predication skill of the MJO will be reported in the future.     

Response of hydrological cycle to recent climate changes in the Tibetan Plateau

The Tibetan Plateau (TP) surfaces have been experiencing an overall rapid warming and wetting while wind speed and solar radiation have been declining in the last three decades. This study investigated how climate changes influenced the hydrological cycle on the TP during 1984??2006. To facilitate the analysis, a land surface model was used to simulate surface water budget at all CMA (China Meteorological Administration) stations on the TP. The simulated results were first validated against observed ground temperature and observation-derived heat flux on the western TP and observed discharge trends on the eastern TP. The response of evaporation and runoff to the climate changes was then analyzed. Major finding are as follows. (1) Surface water balance has been changed in recent decades. Observed precipitation shows insignificant increasing trends in central TP and decreasing trends along the TP periphery while evaporation shows overall increasing trends, leading to decreased discharge at major TP water resource areas (semi-humid and humid zones in the eastern and southern TP). (2) At the annual scale, evaporation is water-limited in dry areas and energy-limited (radiation and air temperature) in wet areas; these constraints can be interpreted by the Budyko-curve. Evaporation in autumns and winters was strongly controlled by soil water storage in summers, weakening the dependence of evaporation on precipitation at seasonal scales. (3) There is a complementary effect between the simulated actual evaporation and potential evaporation, but this complementary relationship may deviate from Bouchet??s hypothesis when vapor pressure deficit (or air temperature) is too low, which suppresses the power of vapor transfer.