Improving the precipitation forecast over the Eastern Mediterranean using a smoothed time‐lagged ensemble

[Departement_IRSTEA]Eaux [TR1_IRSTEA]ARCEAU

Since 1970, tropical cyclone (TC) track forecasts have improved steadily in the Atlantic basin. This improvement has been linked primarily to advances in numerical weather prediction (NWP) models. Concurrently, with few exceptions, the development and operational use of statistical track prediction schemes have experienced a relative decline. Statistical schemes provided the most accurate TC track forecasts until approximately the late 1980s. In this note, it is shown that increased reliance on the global NWP models does not always guarantee the best forecast. Here, Hurricane Ivan is used from the 2004 Atlantic TC season as a classical example, and reminder, of how strong climatological signals still can add substantial value to TC track forecasts, in the form of improved accuracy and increased timeliness at minimal computational cost. In an 8-day period in early September 2004, Hurricane Ivan was repeatedly, and incorrectly, forecast by 12 operational NWP models to move with a significant northward (poleward) component. It was found that the mean 24-h trajectory forecasts of a consensus of five commonly used NWP track prediction aids had a statistically significant right-of-track bias. Furthermore, the official track forecasts, which relied heavily on erroneous numerical guidance over this period, were also found to have significant poleward trajectory errors. At the same time, a climatology-based prediction technique, drawn entirely from the historical record of motion characteristics of TCs in geographical locations similar to Ivan, correctly and consistently indicated a more westward motion component, had a small directional spread, and was supported by a large number of archived cases. This climatological signal was in conflict with the deterministic NWP model output, and it is suggested that the large errors in the official track forecast for TC Ivan could have been reduced considerably by taking into greater account such a strong climatological signal. The potential impact of such an error reduction is a saving of lives and billions of dollars in both actual damage and unnecessary evacuations costs, for just this one hurricane. We also suggest that this simple strategy of examining the strength of the climatological signal be considered for all TCs to identify cases where the NWP and official forecasts differ significantly from strong, persistent climatological signals.

A comprehensive assessment of the forecast skill of various meteorological indices over the Indian region from different numerical weather prediction (NWP) models is lacking in the literature. In this study the performance of four NWP models, namely Global Ensemble Forecast System (GEFSv12), European Center for Medium Range Forecasting (ECMWF), Climate Forecast System (CFSv2) and Indian Institute of Tropical Meteorology (IITM) towards forecasting of precipitation, temperature and associated meteorological indices, is evaluated at short to medium timescales across the Indian region. Further, the effect of ocean atmospheric (OA) oscillations on the precipitation/temperature forecast skill from the different NWP models is also assessed. Results show that the NWP models are better in predicting the meteorological indices than the quantitative forecasts. The ECMWF model was found to be the best for PCP forecasting with CFSv2 performing poorly. For temperature the GEFSv12 model performance was the lowest, compared to rest of the models. The models show poor skill in forecasting monsoon season precipitation compared to non‐monsoon and the temperature forecasts from the NWP models are particularly poor for the northern basins. Skilful temperature forecasts are observed in the Northwestern, Indo‐Gangetic and Central basins for the CFSv2 and ECMWF models. The forecast skill of precipitation indices are higher in the northwestern, central and Indo‐Gangetic basins compared to the rest. The skill of precipitation indices, namely rainy days, extreme rainy days and consecutive wet days, is higher during the monsoon seasons while the prediction skill of consecutive dry days is higher during the non‐monsoon season. OA analysis revealed that the ENSO phases have dominant effect on the forecast skill of the precipitation only. The temperature and meteorological indices forecasts are not affected significantly by the OA phases. The outcomes of this study have implications towards irrigation scheduling and water resources management decision making in India.

Abstract. The current operational very short-term and short-term quantitative precipitation forecast (QPF) at the Meteorological Service of Catalonia (SMC) is made by three different methodologies: Advection of the radar reflectivity field (ADV), Identification, tracking and forecasting of convective structures (CST) and numerical weather prediction (NWP) models using observational data assimilation (radar, satellite, etc.). These precipitation forecasts have different characteristics, lead time and spatial resolutions. The objective of this study is to combine these methods in order to obtain a single and optimized QPF at each lead time. This combination (blending) of the radar forecast (ADV and CST) and precipitation forecast from NWP model is carried out by means of different methodologies according to the prediction horizon. Firstly, in order to take advantage of the rainfall location and intensity from radar observations, a phase correction technique is applied to the NWP output to derive an additional corrected forecast (MCO). To select the best precipitation estimation in the first and second hour (t+1 h and t+2 h), the information from radar advection (ADV) and the corrected outputs from the model (MCO) are mixed by using different weights, which vary dynamically, according to indexes that quantify the quality of these predictions. This procedure has the ability to integrate the skill of rainfall location and patterns that are given by the advection of radar reflectivity field with the capacity of generating new precipitation areas from the NWP models. From the third hour (t+3 h), as radar-based forecasting has generally low skills, only the quantitative precipitation forecast from model is used. This blending of different sources of prediction is verified for different types of episodes (convective, moderately convective and stratiform) to obtain a robust methodology for implementing it in an operational and dynamic way.

Forecasts from numerical weather prediction (NWP) models play a critical role in many sectors of the U.S. economy. Improvements to operational NWP model forecasts are generally assumed to provide significant economic savings through better decision-making. But is this true? Since 2014, several new versions of the High-Resolution Rapid Refresh (HRRR) model were released into operation within the National Weather Service. Practically, forecasts have an economic impact only if they lead to a different action than what would be taken under an alternative information set. And in many sectors, these decisions only need to be considered during certain weather conditions. We estimate the economic impacts of improvements made to the HRRR, using 12-h wind, precipitation, and temperature forecasts in several cases where they can have “economically meaningful” behavioral consequences. We examine three different components of the U.S. economy where such information matters: 1) better integration of wind energy resources into the electric grid, 2) increased worker output due to better precipitation forecasts that allow workers to arrive to their jobs on time, and 3) better decisions by agricultural producers in preparing for freezing conditions. These applications demonstrate some of the challenges in ascertaining the economic impacts of improved weather forecasts, including highlighting key assumptions that must be made to make the problem tractable. For these sectors, we demonstrate that there was a marked economic gain for the United States between HRRR versions 1 and 2 and a smaller, but still appreciable economic gain between versions 2 and 3.

Accuracy of humidity forecasts has been considered relatively unimportant to much of the operational numerical weather prediction (NWP) community. However, the U.S. Air Force is interested in accurate water vapor and cloud forecasts as end products. It is expected that the NWP community as a whole will become more involved in improving their humidity forecasts as they recognize the important role of accurate water vapor distributions in data assimilation, forecasts of temperature and precipitation, and climate change research. As a modeling community, we need to begin now to identify and rectify the systematic humidity forecast errors that are present in NWP models. This will allow us to take full advantage of the new types of remotely sensed water vapor and cloud measurements that are on the horizon. The research reported in this paper attempts to address this issue in a simple, straightforward manner, using the Phillips Laboratory Global Spectral Model (PL GSM). It was found that significant sy...

We construct the generalized inverse of a global numerical weather prediction (NWP) model, in order to prepare initial conditions for the model at time “t=0 hrs”. The inverse finds a weighted, least-squares best-fit to the dynamics for −24<t<0, to the previous initial condition att=−24, and to data att=−24,t=−18,t=−12 andt=0. That is, the inverse is a weak-constraint, four-dimensional variational assimilation scheme. The best-fit is found by solving the nonlinear Euler-Lagrange (EL) equations which determine the local extrema of a penalty functional. The latter is quadratic in the dynamical, initial and data residuals. The EL equations are solved using iterated representer expansions. The technique yields optimal conditioning of the very large minimization problem, which has ∼109 hydrodynamical and thermodynamical variables defined on a 4-dimensional, space-time grid.

Abstract. Flash floods are among the most devastating and lethal natural hazards. In 2018, three flash-flood episodes resulted in 46 casualties in the deserts of Israel and Jordan alone. This paper presents the hydrometeorological analysis and forecasting of a substantial storm (25–27 April 2018) that hit an arid desert basin (Zin, ∼1400 km2, southern Israel) claiming 12 human lives. This paper aims to (a) spatially assess the severity of the storm, (b) quantify the timescale of the hydrological response, and (c) evaluate the available operational precipitation forecasting. Return periods of the storm's maximal rain intensities were derived locally at 1 km2 resolution using weather radar data and a novel statistical methodology. A high-resolution grid-based hydrological model was used to study the intra-basin flash-flood magnitudes which were consistent with direct information from witnesses. The model was further used to examine the hydrological response to different forecast scenarios. A small portion of the basin (1 %–20 %) experienced extreme precipitation intensities (75- to 100-year return period), resulting in a local hydrological response of a high magnitude (10- to 50-year return period). Hillslope runoff, initiated minutes after the intense rainfall occurred, reached the streams and resulted in peak discharge within tens of minutes. Available deterministic operational precipitation forecasts poorly predicted the hydrological response in the studied basins (tens to hundreds of square kilometers) mostly due to location inaccuracy. There was no gain from assimilating radar estimates in the numerical weather prediction model. Therefore, we suggest using deterministic forecasts with caution as it might lead to fatal decision making. To cope with such errors, a novel cost-effective methodology is applied by spatially shifting the forecasted precipitation fields. In this way, flash-flood occurrences were captured in most of the subbasins, resulting in few false alarms.

Numerical weather prediction (NWP) models produce a quantitative precipitation forecast (QPF), which is vital for a wide range of applications, especially for accurate flash flood forecasting. The under- and over-estimation of forecast uncertainty pose operational risks and often encourage overly conservative decisions to be made. Since NWP models are subject to many uncertainties, the QPFs need to be post-processed. The NWP biases should be corrected prior to their use as a reliable data source in hydrological models. In recent years, several post-processing techniques have been proposed. However, there is a lack of research on post-processing the real-time forecast of NWP models considering bias lead-time dependency for short- to medium-range forecasts. The main objective of this study is to use the total least squares (TLS) method and the lead-time dependent bias correction method—known as dynamic weighting (DW)—to post-process forecast real-time data. The findings show improved bias scores, a decrease in the normalized error and an improvement in the scatter index (SI). A comparison between the real-time precipitation and flood forecast relative bias error shows that applying the TLS and DW methods reduced the biases of real-time forecast precipitation. The results for real-time flood forecasts for the events of 2002, 2007 and 2011 show error reductions and accuracy improvements of 78.58%, 81.26% and 62.33%, respectively.

Artificial intelligence and numerical weather prediction models: A technical survey

Two methods for assimilating radar reflectivity into the COSMO numerical weather prediction (NWP) model were compared to precipitation forecasts. The first method assimilated observed radar reflectivity, and the second one assimilated observed and extrapolated radar reflectivity. The assimilation technique was based on the correction of the model's water vapour mixing ratio. The extrapolation was performed by the COTREC method and was 1 hour long. The model's horizontal resolution was 2.8 km.The comparison of methods was based on verification of the observed and forecast hourly precipitation. The comparison was performed for the 1st, 2nd and 3rd hours of each forecast. On the whole, 45 forecasts from nine days of convective precipitation were evaluated for each hour. The evaluation included subjective verification and the following objective skill scores: Fractions Skill Scores, SAL and a measure based on a categorical‐probabilistic approach.The results confirmed that assimilation complemented by the extrapolated data improves the accuracy of precipitation forecasts. The improvement was obvious in a majority of the single forecasts studied, and it is confirmed by all evaluation techniques. COSMO forecasts that used the extrapolation showed reasonable competence in forecasting for the first and the second hours. Copyright © 2011 Royal Meteorological Society

Improving the Turbulence Coupling between High Resolution Numerical Weather Prediction Models and Lagrangian Particle Dispersion Models

Calibrated high-temporal-resolution precipitation forecasts are desirable for a range of applications, for example, flood prediction in fast-rising rivers. However, high-temporal-resolution precipitation observations may not be available to support the establishment of calibration methods, particularly in regions with low population density or in developing countries. We present a new method to produce calibrated hourly precipitation ensemble forecasts from daily observations. Precipitation forecasts are taken from a high-resolution convective-scale numerical weather prediction (NWP) model run at the hourly time step. We conduct three experiments to develop the new calibration method: (i) calibrate daily precipitation totals and disaggregate daily forecasts to hourly; (ii) generate pseudohourly observations from daily precipitation observations, and use these to calibrate hourly precipitation forecasts; and (iii) combine aspects of (i) and (ii). In all experiments, we use the existing Bayesian joint probability model to calibrate the forecasts and the well-known Schaake shuffle technique to instill realistic spatial and temporal correlations in the ensembles. As hourly observations are not available, we use hourly patterns from the NWP as the template for the Schaake shuffle. The daily member matching method (DMM), method (iii), produces the best-performing ensemble precipitation forecasts over a range of metrics for forecast accuracy, bias, and reliability. The DMM method performs very similarly to the ideal case where hourly observations are available to calibrate forecasts. Overall, valuable spatial and temporal information from the forecast can be extracted for calibration with daily data, with a slight trade-off between forecast bias and reliability.

Improving lead time for forecasting floods is important to minimize property damage and ensure the safety of the public and emergency services during flood events. Numerical weather prediction (NWP) models are important components of flood forecasting systems and have been vital in extending forecasting lead time under complex weather and terrain conditions. However, NWP forecasts still have significant uncertainty associated with the precipitation fields that are the main inputs of the hydrologic models and thus the resulting flood forecasts. An issue often overlooked is the importance of correctly representing variability over a range of different temporal scales. To address this gap, here a new wavelet-based method for postprocessing NWP precipitation forecasts is proposed. First, precipitation forecasts are decomposed into the frequency domain using a wavelet transform, providing estimates of the amplitudes and phases of the time series at different frequencies. Quantile mapping is then used to correct bias in the amplitudes of each frequency. Randomized phases are used to generate an ensemble of realizations of the precipitation forecasts. The postprocessed precipitation forecasts are reconstructed by taking the inverse of adjusted time-frequency decompositions with the corrected amplitudes and randomized phases. The proposed method was used to postprocess NWP precipitation forecasts in the Sydney region of Australia. There is a significant improvement in postprocessed precipitation forecasts across multiple time scales in terms of bias and temporal and spatial correlation structures. The postprocessed precipitation fields can be used for the modeling of fully distributed hydrologic systems, improving runoff stimulation, flood depth estimation, and flood early warning. Significance Statement A new method accounting for the timing and spatial errors of NWP precipitation forecasts is proposed, and it can improve the skill of forecasts across multiple time scales, especially at short lead times. The proposed method provides a practical and effective way to correct these errors by incorporating spatiotemporal neighborhood information through the frequency domain using sophisticated wavelet transforms. With systematic timing and spatial errors removed, precipitation forecasts will be more skillful, and hydrological modeling using the postprocessed forecasts can provide higher accuracy of streamflow estimation.

Accurate wind speed forecasting is essential for daily activities and socio‐economic production. This study employs an explainable machine learning (ML) model to correct wind speed forecasts from a numerical weather prediction (NWP) model over weather stations across lead times of 1–24 hr in Zhejiang Province, China. Eleven meteorological variables derived from the NWP model are used as predictors. The results demonstrate that ML model substantially improves the original NWP forecasts, reducing the root mean square error of approximately 20%–40% at most stations and about 10% at certain coastal stations. The corrected forecasts also consistently achieve higher Kling‐Gupta Efficiency scores and Forecast Accuracy across all lead times. To enhance model interpretability, we apply Split Count and SHapley Additive exPlanation analyses. The results indicate that the friction velocity (UST) is the most influential predictor and reveal a threshold beyond which the ML model's correction performance improves markedly. The ML model performance slightly declines in coastal regions, likely due to the inconsistent representation of UST between the land and water surface in the NWP model. Monthly analysis further shows superior model performance in November and December, when UST values are higher. This study underscores the value of explainable ML models in improving NWP wind forecasts and provides insights for correcting meteorological variables in NWP models.