Ensemble Forecasting of Tropical Cyclone Motion Using a Baroclinic Model

The purpose of this study is to investigate the effectiveness of two different ensemble forecasting （EF） techniques-the lagged-averaged forecast （LAF） and the breeding of growing modes （BGM）. In the BGM experiments, the vortex and the environment are perturbed separately （named BGMV and BGME）. Tropical cyclone （TC） motions in two difficult situations are studied： a large vortex interacting with its environment, and an apparent binary interaction. The former is Typhoon Yancy and the latter involves Typhoon Ed and super Typhoon Flo, all occurring during the Tropical Cyclone Motion Experiment TCM- 90. The model used is the baroclinic model of the University of New South Wales. The lateral boundary tendencies are computed from atmospheric analysis data. Only the relative skill of the ensemble forecast mean over the control run is used to evaluate the effectiveness of the EF methods, although the EF technique is also usecl to quantify forecast uncertainty in some studies. In the case of Yancy, the ensemble mean forecasts of each of the three methodologies are better than that of the control, with LAF being the best. The mean track of the LAF is close to the best track, and it predicts landfall over Taiwan. The improvements in LAF and the full BGM where both the environment and vortex are perturbed suggest the importance of combining the perturbation of the vortex and environment when the interaction between the two is appreciable. In the binary interaction case of Ed and Flo, the forecasts of Ed appear to be insensitive to perturbations of the environment and/or the vortex, which apparently results from erroneous forecasts by the model of the interaction between the subtropical ridge and Ed, as well as from the interaction between the two typhoons, thus reducing the effectiveness of the EF technique. This conclusion is reached through sensitivity experiments on the domain of the model and by adding or eliminating certain features in the model atmosphere. Nevertheless, the forecast tracks in some of the cases are improved over that of the control. On the other hand, the EF technique has little impact on the forecasts of Flo because the control forecast is already very close to the best track. The study provides a basis for the. future development of the EF technique. The limitations of this study are also addressed. For example, the above results are based on a small sample, and the study is actually a simulation, which is different than operational forecasting. Further tests of these EF techniques are proposed.

The Simulation of Five Tropical Cyclones by Sample Optimization of Ensemble Forecasting Based on the Observed Track and Intensity

The quality of ensemble forecasting is seriously affected by sample quality.In this study,the distributions of ensemble members based on the observed track and intensity of tropical cyclones(TCs)were optimized and their influence on the simulation results was analyzed.Simulated and observed tracks and intensities of TCs were compared and these two indicators were combined and weighted to score the sample.Samples with higher scores were retained and samples with lower scores were eliminated to improve the overall quality of the ensemble forecast.For each sample,the track score and intensity score were added as the final score of the sample with weight proportions of 10 to 0,9 to 1,8 to 2,7 to 3,6 to 4,5 to 5.These were named as“tr”,“91”,“82”,“73”,“64”,and“55”,respectively.The WRF model was used to simulate five tropical cyclones in the northwestern Pacific to test the ability of this scheme to improve the forecast track and intensity of these cyclones.The results show that the sample optimization effectively reduced the track and intensity error,“55”usually had better performance on the short-term intensity prediction,and“tr”had better performance in short-term track prediction.From the overall performance of the track and intensity simulation,“91”was the best and most stable among all sample optimization schemes.These results may provide some guidance for optimizing operational ensemble forecasting of TCs.

Dynamic Analogue Initialization for Ensemble Forecasting

This paper introduces a new approach for the initialization of ensemble numerical forecasting： Dynamic Analogue Initialization （DAI）. DAI assumes that the best model state trajectories for the past provide the initial conditions for the best forecasts in the future. As such, DAI performs the ensemble forecast using the best analogues from a full size ensemble. As a pilot study, the Lorenz63 and Lorenz96 models were used to test DAI＇s effectiveness independently. Results showed that DAI can improve the forecast significantly. Especially in lower-dimensional systems, DAI can reduce the forecast RMSE by ~50% compared to the Monte Carlo forecast （MC）. This improvement is because DAI is able to recognize the direction of the analysis error through the embedding process and therefore selects those good trajectories with reduced initial error. Meanwhile, a potential improvement of DAI is also proposed, and that is to find the optimal range of embedding time based on the error＇s growing speed.

THE ENSEMBLE FORECASTING OF TROPICAL CYCLONE MOTION I:USING A PRIMITIVE EQUATION BAROTROPIC MODEL

Ensemble forecasting of tropical cyclone (TC) motion was studied using a primitive equation barotropic model by perturbing initial position and structure for 1979 1993 TC. The results show that TC initial position perturbation affects its track, but the ensemble mean is close to control forecast. Experiments was also performed by perturbing TC initial parameters which were used to generate TC initial field, and more improvement can be obtained by taking ensemble mean of selective member than selecting members randomly. The skill of 60 % 70 % of all cases is improved in selective ensemble mean. When the ambient steering current is weak, more improvement can be obtained over the control forecast.

Non-crossing Quantile Regression Neural Network as a Calibration Tool for Ensemble Weather Forecasts

Despite the maturity of ensemble numerical weather prediction(NWP),the resulting forecasts are still,more often than not,under-dispersed.As such,forecast calibration tools have become popular.Among those tools,quantile regression(QR)is highly competitive in terms of both flexibility and predictive performance.Nevertheless,a long-standing problem of QR is quantile crossing,which greatly limits the interpretability of QR-calibrated forecasts.On this point,this study proposes a non-crossing quantile regression neural network(NCQRNN),for calibrating ensemble NWP forecasts into a set of reliable quantile forecasts without crossing.The overarching design principle of NCQRNN is to add on top of the conventional QRNN structure another hidden layer,which imposes a non-decreasing mapping between the combined output from nodes of the last hidden layer to the nodes of the output layer,through a triangular weight matrix with positive entries.The empirical part of the work considers a solar irradiance case study,in which four years of ensemble irradiance forecasts at seven locations,issued by the European Centre for Medium-Range Weather Forecasts,are calibrated via NCQRNN,as well as via an eclectic mix of benchmarking models,ranging from the naïve climatology to the state-of-the-art deep-learning and other non-crossing models.Formal and stringent forecast verification suggests that the forecasts post-processed via NCQRNN attain the maximum sharpness subject to calibration,amongst all competitors.Furthermore,the proposed conception to resolve quantile crossing is remarkably simple yet general,and thus has broad applicability as it can be integrated with many shallow-and deep-learning-based neural networks.

Scientific Advances and Weather Services of the China Meteorological Administration’s National Forecasting Systems during the Beijing 2022 Winter Olympics

Since the Beijing 2022 Winter Olympics was the first Winter Olympics in history held in continental winter monsoon climate conditions across complex terrain areas,there is a deficiency of relevant research,operational techniques,and experience.This made providing meteorological services for this event particularly challenging.The China Meteorological Administration(CMA)Earth System Modeling and Prediction Centre,achieved breakthroughs in research on short-and medium-term deterministic and ensemble numerical predictions.Several key technologies crucial for precise winter weather services during the Winter Olympics were developed.A comprehensive framework,known as the Operational System for High-Precision Weather Forecasting for the Winter Olympics,was established.Some of these advancements represent the highest level of capabilities currently available in China.The meteorological service provided to the Beijing 2022 Games also exceeded previous Winter Olympic Games in both variety and quality.This included achievements such as the“100-meter level,minute level”downscaled spatiotemporal resolution and forecasts spanning 1 to 15 days.Around 30 new technologies and over 60 kinds of products that align with the requirements of the Winter Olympics Organizing Committee were developed,and many of these techniques have since been integrated into the CMA’s operational national forecasting systems.These accomplishments were facilitated by a dedicated weather forecasting and research initiative,in conjunction with the preexisting real-time operational forecasting systems of the CMA.This program represents one of the five subprograms of the WMO’s high-impact weather forecasting demonstration project(SMART2022),and continues to play an important role in their Regional Association(RA)II Research Development Project(Hangzhou RDP).Therefore,the research accomplishments and meteorological service experiences from this program will be carried forward into forthcoming highimpact weather forecasting activities.This article provides an overview and assessment of this program and the operational national forecasting systems.

Growth and Interactions of Multi-Source Perturbations in Convection-Allowing Ensemble Forecasts

This study investigated the growth of forecast errors stemming from initial conditions(ICs),lateral boundary conditions(LBCs),and model(MO)perturbations,as well as their interactions,by conducting seven 36 h convectionallowing ensemble forecast(CAEF)experiments.Two cases,one with strong-forcing(SF)and the other with weak-forcing(WF),occurred over the Yangtze-Huai River basin(YHRB)in East China,were selected to examine the sources of uncertainties associated with perturbation growth under varying forcing backgrounds and the influence of these backgrounds on growth.The perturbations exhibited distinct characteristics in terms of temporal evolution,spatial propagation,and vertical distribution under different forcing backgrounds,indicating a dependence between perturbation growth and forcing background.A comparison of the perturbation growth in different precipitation areas revealed that IC and LBC perturbations were significantly influenced by the location of precipitation in the SF case,while MO perturbations were more responsive to convection triggering and dominated in the WF case.The vertical distribution of perturbations showed that the sources of uncertainties and the performance of perturbations varied between SF and WF cases,with LBC perturbations displaying notable case dependence.Furthermore,the interactions between perturbations were considered by exploring the added values of different source perturbations.For the SF case,the added values of IC,LBC,and MO perturbations were reflected in different forecast periods and different source uncertainties,suggesting that the combination of multi-source perturbations can yield positive interactions.In the WF case,MO perturbations provided a more accurate estimation of uncertainties downstream of the Dabie Mountain and need to be prioritized in the research on perturbation development.

Algorithm based on local breeding of growing modes for convection-allowing ensemble forecasting

We propose a method based on the local breeding of growing modes(LBGM) considering strong local weather characteristics for convection-allowing ensemble forecasting. The impact radius was introduced in the breeding of growing modes to develop the LBGM method. In the local breeding process, the ratio between the root mean square error(RMSE) of local space forecast at each grid point and that of the initial full-field forecast is computed to rescale perturbations. Preliminary evaluations of the method based on a nature run were performed in terms of three aspects: perturbation structure, spread,and the RMSE of the forecast. The experimental results confirm that the local adaptability of perturbation schemes improves after rescaling by the LBGM method. For perturbation physical variables and some near-surface meteorological elements, the LBGM method could increase the spread and reduce the RMSE of forecast,improving the performance of the ensemble forecast system.In addition, different from those existing methods of global orthogonalization approach, this new initial-condition perturbation method takes into full consideration the local characteristics of the convective-scale weather system, thus making convectionallowing ensemble forecast more accurate.

CHAR ACTERISTICS OF TROPICAL CYCLONE GENESIS FORECASTS AND UNDERDISPERSION IN HIGH-RESOLUTION ENSEMBLE FORECASTING WITH A STOCHASTIC KINETIC ENERGY BACKSCATTER SCHEME

This study evaluates ensemble forecasts with a stochastic kinetic energy backscatter scheme(SKEBS)to predict tropical cyclone(TC)genesis and also to characterize the related ensemble underdispersion.Several sets of ensemble forecasts are generated using an advanced research version of the Weather Research and Forecasting model at 5 km horizontal resolution to predict the genesis of Hurricane Ernesto(2006)and Typhoon Nuri(2008).Ensemble forecasts with SKEBS are compared against a control ensemble forecast with the WRF model using downscaled initial conditions derived from the NCEP Global Ensemble Forecasting System.It is found that ensemble forecasts with SKEBS are able to generate probabilistic forecasts for TC genesis and also capable of indicating the forecast uncertainties.Compared with the deterministic forecast that fails to predict the genesis of Typhoon Nuri,the ensemble forecast with SKEBS is able to produce the genesis forecast.However,the underdispersion of ensemble forecasts with SKEBS is also present in all cases in terms of the simulation period and over the whole model domain,TC environment,and inner core regions,although it is reduced near the TC inner core region.In addition,the initial perturbation–based ensemble forecasts shows slightly less underdispersion compared with the SKEBS ensembles.

The method of residual‑based bootstrap averaging of the forecast ensemble

This paper presents an optimization approach—residual-based bootstrap averaging(RBBA)—for different types of forecast ensembles.Unlike traditional residual-mean-square-error-based ensemble forecast averaging approaches,the RBBA method attempts to find optimal forecast weights in an ensemble and allows for their combi-nation into the most effective additive forecast.In the RBBA method,all the different types of forecasts obtain the optimal weights for ensemble residuals that are statisti-cally optimal in terms of the fitness function of the residuals.Empirical studies have been conducted to demonstrate why and how the RBBA method works.The experi-mental results based on the real-world time series of contemporary stock exchanges show that the RBBA method can produce ensemble forecasts with good generalization ability.

Comparison of Nonlinear Local Lyapunov Vectors with Bred Vectors, Random Perturbations and Ensemble Transform Kalman Filter Strategies in a Barotropic Model

The breeding method has been widely used to generate ensemble perturbations in ensemble forecasting due to its simple concept and low computational cost. This method produces the fastest growing perturbation modes to catch the growing components in analysis errors. However, the bred vectors （BVs） are evolved on the same dynamical flow, which may increase the dependence of perturbations. In contrast, the nonlinear local Lyapunov vector （NLLV） scheme generates flow-dependent perturbations as in the breeding method, but regularly conducts the Gram-Schmidt reorthonormalization processes on the perturbations. The resulting NLLVs span the fast-growing perturbation subspace efficiently, and thus may grasp more com- ponents in analysis errors than the BVs. In this paper, the NLLVs are employed to generate initial ensemble perturbations in a barotropic quasi-geostrophic model. The performances of the ensemble forecasts of the NLLV method are systematically compared to those of the random pertur- bation （RP） technique, and the BV method, as well as its improved version--the ensemble transform Kalman filter （ETKF） method. The results demonstrate that the RP technique has the worst performance in ensemble forecasts, which indicates the importance of a flow-dependent initialization scheme. The ensemble perturbation subspaces of the NLLV and ETKF methods are preliminarily shown to catch similar components of analysis errors, which exceed that of the BVs. However, the NLLV scheme demonstrates slightly higher ensemble forecast skill than the ETKF scheme. In addition, the NLLV scheme involves a significantly simpler algorithm and less computation time than the ETKF method, and both demonstrate better ensemble forecast skill than the BV scheme.

Numerical simulation of an extreme haze pollution event over the North China Plain based on initial and boundary condition ensembles

The North China Plain often su ers heavy haze pollution events in the cold season due to the rapid industrial development and urbanization in recent decades.In the winter of 2015,the megacity cluster of Beijing Tianjin Hebei experienced a seven-day extreme haze pollution episode with peak PM2.5(particulate matter(PM)with an aerodynamic diameter≤2.5μm)concentration of 727μg m 3.Considering the in uence of meteorological conditions on pollu-tant evolution,the e ects of varying initial conditions and lateral boundary conditions(LBCs)of the WRF-Chem model on PM2.5 concentration variation were investigated through ensemble methods.A control run(CTRL)and three groups of ensemble experiments(INDE,BDDE,INBDDE)were carried out based on difierent initial conditions and LBCs derived from ERA5 reanalysis data and its 10 ensemble members.The CTRL run reproduced the meteorological conditions and the overall life cycle of the haze event reasonably well,but failed to capture the intense oscillation of the instantaneous PM2.5 concentration.However,the ensemble forecasting showed a considerable advantage to some extent.Compared with the CTRL run,the root-mean-square error(RMSE)of PM2.5 concentration decreased by 4.33%,6.91%,and 8.44%in INDE,BDDE and INBDDE,respectively,and the RMSE decreases of wind direction(5.19%,8.89%and 9.61%)were the dominant reason for the improvement of PM2.5 concentration in the three ensemble experiments.Based on this case,the ensemble scheme seems an e ective method to improve the prediction skill of wind direction and PM2.5 concentration by using the WRF-Chem model.

The Relationship between Deterministic and Ensemble Mean Forecast Errors Revealed by Global and Local Attractor Radii

It has been demonstrated that ensemble mean forecasts, in the context of the sample mean, have higher forecasting skill than deterministic(or single) forecasts. However, few studies have focused on quantifying the relationship between their forecast errors, especially in individual prediction cases. Clarification of the characteristics of deterministic and ensemble mean forecasts from the perspective of attractors of dynamical systems has also rarely been involved. In this paper, two attractor statistics—namely, the global and local attractor radii(GAR and LAR, respectively)—are applied to reveal the relationship between deterministic and ensemble mean forecast errors. The practical forecast experiments are implemented in a perfect model scenario with the Lorenz96 model as the numerical results for verification. The sample mean errors of deterministic and ensemble mean forecasts can be expressed by GAR and LAR, respectively, and their ratio is found to approach2^(1/2) with lead time. Meanwhile, the LAR can provide the expected ratio of the ensemble mean and deterministic forecast errors in individual cases.

Ensemble Forecast: A New Approach to Uncertainty and Predictability

Ensemble techniques have been used to generate daily numerical weather forecasts since the 1990s in numerical centers around the world due to the increase in computation ability. One of the main purposes of numerical ensemble forecasts is to try to assimilate the initial uncertainty （initial error） and the forecast uncertainty （forecast error） by applying either the initial perturbation method or the multi-model/multiphysics method. In fact, the mean of an ensemble forecast offers a better forecast than a deterministic （or control） forecast after a short lead time （3-5 days） for global modelling applications. There is about a 1-2-day improvement in the forecast skill when using an ensemble mean instead of a single forecast for longer lead-time. The skillful forecast （65% and above of an anomaly correlation） could be extended to 8 days （or longer） by present-day ensemble forecast systems. Furthermore, ensemble forecasts can deliver a probabilistic forecast to the users, which is based on the probability density function （PDF） instead of a single-value forecast from a traditional deterministic system. It has long been recognized that the ensemble forecast not only improves our weather forecast predictability but also offers a remarkable forecast for the future uncertainty, such as the relative measure of predictability （RMOP） and probabilistic quantitative precipitation forecast （PQPF）. Not surprisingly, the success of the ensemble forecast and its wide application greatly increase the confidence of model developers and research communities.

A Regional Ensemble Forecast System for Stratiform Precipitation Events in the Northern China Region.Part Ⅱ:Seasonal Evaluation for Summer 2010

In this study, the Institute of Atmospheric Physics, Chinese Academy of Sciences - regional ensemble forecast system （IAP-REFS） described in Part I was further validated through a 65-day experiment using the summer season of 2010. The verification results show that IAP-REFS is skillful for quantitative precipitation forecasts （QPF） and probabilistic QPF, but it has a systematic bias in forecasting near-surface variables. Applying a 7-day running mean bias correction to the forecasts of near-surface variables remarkably improved the reliability of the forecasts. In this study, the perturbation extraction and inflation method （proposed with the single case study in Part I） was further applied to the full season with different inflation factors. This method increased the ensemble spread and improved the accuracy of forecasts of precipitation and near-surface variables. The seasonal mean profiles of the IAP-REFS ensemble indicate good spread among ensemble members and some model biases at certain vertical levels.

A Regional Ensemble Forecast System for Stratiform Precipitation Events in Northern China.Part I:A Case Study

A single-model, short-range, ensemble forecasting system （Institute of Atmospheric Physics, Regional Ensemble Forecast System, IAP REFS） with 15-km grid spacing, configured with multiple initial conditions, multiple lateral boundary conditions, and multiple physics parameterizations with 11 ensemble members, was developed using the Weather and Research Forecasting Model Advanced Research modeling system for prediction of stratiform precipitation events in northern China. This is the first part of a broader research project to develop a novel cloud-seeding operational system in a probabilistic framework. The ensemble perturbations were extracted from selected members of the National Center for Environmental Prediction Global Ensemble Forecasting System （NCEP GEFS） forecasts, and an inflation factor of two was applied to compensate for the lack of spread in the GEFS forecasts over the research region. Experiments on an actual stratiform precipitation case that occurred on 5-7 June 2009 in northern China were conducted to validate the ensemble system. The IAP REFS system had reasonably good performance in predicting the observed stratiform precipitation system. The perturbation inflation enlarged the ensemble spread and alleviated the underdispersion caused by parent forecasts. Centering the extracted perturbations on higher-resolution NCEP Global Forecast System forecasts resulted in less ensemble mean root-mean-square error and better accuracy in probabilistic quantitative precipitation forecasts （PQPF）. However, the perturbation inflation and recentering had less effect on near-surface-level variables compared to the mid-level variables, and its influence on PQPF resolution was limited as well.

AN ENSEMBLE FORECAST EXPERIMENT OF A LANDING TYPHOON

Based on the Global Regional Assimilation and Prediction System-Tropical Cyclone Model(GRAPES-TCM),an ensemble forecast experiment was performed,in which Typhoon Wipha during the period immediately prior to landfall was selected for the study and the breeding of growing mode(BGM) method was used to perturb the initial conditions of the vortex field and the environment field.The results of the experiment indicate that each member had a different initial status in BGM processing and they show a reasonable spread among members along with the forecast phase.Changes in the large-scale field,thermodynamic structure,and spread among members took place when Wipha made landfall.The steering effect of the large-scale field and the interaction between the thermodynamics and the dynamics resulted in different tracks of the members.Meanwhile,the forecast uncertainty increased.In summary,the ensemble mean did not perform as well as the control forecast,but the cluster mean provided some useful information,and performed better than the control in some instances.The position error was 34 km for 24 h forecast,153 km for 48 h forecast,and 191 km for 66 h forecast.The strike probability chart qualitatively described the forecast uncertainty.

A STUDY ON THE ENSEMBLE FORECAST REAL-TIME CORRECTION METHOD

Using real-time correction technology for typhoons, this paper discusses real-time correction for forecasting the track of four typhoons during 2009 and 2010 in Japan, Beijing, Guangzhou, and Shanghai. It was determined that the short-time forecast effect was better than the original objective mode. By selecting four types of integration schemes after multiple mode path integration for those four objective modes, the forecast effect of the multi-mode path integration is better, on average, than any single model. Moreover, multi-mode ensemble forecasting has obvious advantages during the initial 36 h.

Evaluation of TIGGE Ensemble Forecasts of Precipitation in Distinct Climate Regions in Iran

The application of numerical weather prediction （NWP） products is increasing dramatically. Existing reports indicate that ensemble predictions have better skill than deterministic forecasts. In this study, numerical ensemble precipitation forecasts in the TIGGE database were evaluated using deterministic, dichotomous （yes/no）, and probabilistic techniques over Iran for the period 2008-16. Thirteen rain gauges spread over eight homogeneous precipitation regimes were selected for evaluation. The Inverse Distance Weighting and Kriging methods were adopted for interpolation of the prediction values, downscaled to the stations at lead times of one to three days. To enhance the forecast quality, NWP values were post-processed via Bayesian Model Averaging. The results showed that ECMWF had better scores than other products. However, products of all centers underestimated precipitation in high precipitation regions while overestimating precipitation in other regions. This points to a systematic bias in forecasts and demands application of bias correction techniques. Based on dichotomous evaluation, NCEP did better at most stations, although all centers overpredicted the number of precipitation events. Compared to those of ECMWF and NCER UKMO yielded higher scores in mountainous regions, but performed poorly at other selected stations. Furthermore, the evaluations showed that all centers had better skill in wet than in dry seasons. The quality of post-processed predictions was better than those of the raw predictions. In conclusion, the accuracy of the NWP predictions made by the selected centers could be classified as medium over Iran, while post-processing of predictions is recommended to improve the quality.

Evaluation of WRF-based Convection-Permitting Multi-Physics Ensemble Forecasts over China for an Extreme Rainfall Event on 21 July 2012 in Beijing

On 21 July 2012,an extreme rainfall event that recorded a maximum rainfall amount over 24 hours of 460 mm,occurred in Beijing,China. Most operational models failed to predict such an extreme amount. In this study,a convective-permitting ensemble forecast system（CEFS）,at 4-km grid spacing,covering the entire mainland of China,is applied to this extreme rainfall case. CEFS consists of 22 members and uses multiple physics parameterizations. For the event,the predicted maximum is 415 mm d^-1 in the probability-matched ensemble mean. The predicted high-probability heavy rain region is located in southwest Beijing,as was observed. Ensemble-based verification scores are then investigated. For a small verification domain covering Beijing and its surrounding areas,the precipitation rank histogram of CEFS is much flatter than that of a reference global ensemble. CEFS has a lower（higher） Brier score and a higher resolution than the global ensemble for precipitation,indicating more reliable probabilistic forecasting by CEFS. Additionally,forecasts of different ensemble members are compared and discussed. Most of the extreme rainfall comes from convection in the warm sector east of an approaching cold front. A few members of CEFS successfully reproduce such precipitation,and orographic lift of highly moist low-level flows with a significantly southeasterly component is suggested to have played important roles in producing the initial convection. Comparisons between good and bad forecast members indicate a strong sensitivity of the extreme rainfall to the mesoscale environmental conditions,and,to less of an extent,the model physics.