Lightning super-bolts in Eastern Mediterranean winter thunderstorms

Lightning is a naturally occurring phenomenon with significant socioeconomic and environmental impacts, underlining the need for the application of advanced methods to predict it. This study aims to evaluate the performance of the Weather Research and Forecasting (WRF) model in the prediction of the lightning potential index (LPI), the lightning threat index (LTI), and the product of convective available potential energy (CAPE) and precipitation (CAPE × P). Based on the intensity of convection and the number of lightning flashes, we simulated seven thundercloud events that occurred in Tehran between 2008 and 2013. The simulated LPI, LTI, and CAPE × P values are compared both qualitatively and quantitatively against ground‐based lighting data obtained from the world wide lightning location network (WWLLN). LPI, LTI, and CAPE × P predict the location of lightning with relatively good accuracy, although LPI outperforms LTI and CAPE × P. We also compared the simulated area‐averaged LPI, LTI, and CAPE × P against the number of lightning flashes from the WWLLN data, based on which LPI shows a better performance. Overall, LPI can be used as an effective index for the prediction of lightning. As LPI is based on cloud microphysics and LTI and CAPE × P are thermodynamic‐based methods, a better performance of LPI implies that a method based on the microphysical approach can better predict lightning. The performance of LPI, LTI, and CAPE × P is also sensitive to the horizontal resolution of model simulations, with a considerable improvement in simulations with higher horizontal resolutions.

Rainfall dynamics play a vital role in tropical peatland by providing sufficient water to keep peat moist throughout the year. Therefore, information of rainfall data either historical or forecasting data has risen in recent decades especially for an alert system of fire. Here the Weather and Research Forecasting (WRF) model may act as a tool to provide forecasting weather data. This study aims to do parameterization on WRF parameters for peatland in Sumatra, and to perform bias correction on the WRF’s rainfall output with observed data. We performed stepwise calibration to choose the best five physical schemes of WRF for use in the study area. The output WRF’s rainfall was bias corrected by spatially observed rainfall data for 2019 at day resolution. Our results showed the following schemes namely (i) Eta scheme for cloud microphysical parameters; (ii) GD scheme for cumulus cloud parameters, (iii) MYJ scheme for planetary boundary layer parameters; (iv) RRTM for longwave radiation; and (v) New Goddard schemes for shortwave radiation are best combination for being used to predict rainfall in maritime continent. The spatially interpolated observed rainfall with the Inverse Distance Weighting (IDW) was outperformed for calibration process of WRF’s rainfall as shown by statistical indicators used in this study. Further, the findings have contributed to advance knowledge of rainfall forecasting in maritime continent, particularly in providing data to support the development of fire danger rating system for Indonesian peatland.

Abstract. Italy is a territory characterized by complex topography with the Apennines mountain range crossing the entire peninsula and its highest peaks in central Italy. Using the latter as our area of interest and the snow seasons 2018/19, 2019/20 and 2020/21, the goal of this study is to investigate the ability of a simple single-layer and a more sophisticated multi-layer snow cover numerical model to reproduce the observed snow height, snow water equivalent and snow extent in the central Apennines, using for both models the same forecast weather data as meteorological forcing. We here consider two well-known ground surface and soil models: (i) Noah LSM, an Eulerian model which simulates the snowpack as a bulk single layer, and (ii) Alpine3D, a multi-layer Lagrangian model which simulates the snowpack stratification. We adopt the Weather Research and Forecasting (WRF) model to produce the meteorological data to drive both Noah LSM and Alpine3D at a regional scale with a spatial resolution of 3 km. While Noah LSM is already online-coupled with the WRF model, we develop here a dedicated offline coupling between WRF and Alpine3D. We validate the WRF simulations of surface meteorological variables in central Italy using a dense network of automatic weather stations, obtaining correlation coefficients higher than 0.68, except for wind speed, which suffered from the model underestimation of the real elevation. The performances of both WRF–Noah and WRF–Alpine3D are evaluated by comparing simulated and measured snow height, snow height variation and snow water equivalent, provided by a quality-controlled network of automatic and manual snow stations located in the central Apennines. We find that WRF–Alpine3D can predict better than WRF–Noah the snow height and the snow water equivalent, showing a correlation coefficient with the observations of 0.9 for the former and 0.7 for the latter. Both models show similar performances in reproducing the observed daily snow height variation; nevertheless WRF–Noah is slightly better at predicting large positive variations, while WRF–Alpine3D can slightly better simulate large negative variations. Finally we investigate the abilities of the models in simulating the snow cover area fraction, and we show that WRF–Noah and WRF–Alpine3D have almost equal skills, with both models overestimating it. The equal skills are also confirmed by Jaccard and the average symmetric surface distance indices.

A new method for driving a One-Dimensional Stratiform Cold (1DSC) cloud model with Weather Research and Forecasting (WRF) model outputs was developed by conducting numerical experiments for a typical large-scale stratiform rainfall event that took place on 4–5 July 2004 in Changchun, China. Sensitivity test results suggested that, with hydrometeor profiles extracted from the WRF outputs as the initial input, and with continuous updating of soundings and vertical velocities (including downdraft) derived from the WRF model, the new WRF-driven 1DSC modeling system (WRF-1DSC) was able to successfully reproduce both the generation and dissipation processes of the precipitation event. The simulated rainfall intensity showed a time-lag behind that observed, which could have been caused by simulation errors of soundings, vertical velocities and hydrometeor profiles in the WRF output. Taking into consideration the simulated and observed movement path of the precipitation system, a nearby grid point was found to possess more accurate environmental fields in terms of their similarity to those observed in Changchun Station. Using profiles from this nearby grid point, WRF-1DSC was able to reproduce a realistic precipitation pattern. This study demonstrates that 1D cloud-seeding models do indeed have the potential to predict realistic precipitation patterns when properly driven by accurate atmospheric profiles derived from a regional shortrange forecasting system. This opens a novel and important approach to developing an ensemble-based rain enhancement prediction and operation system under a probabilistic framework concept.

The fire weather of south-east Australia from 1985 to 2009 has been simulated using the Weather Research and Forecasting (WRF) model. The US National Oceanic and Atmospheric Administration Centers for Environmental Prediction and National Center for Atmospheric Research reanalysis supplied the lateral boundary conditions and initial conditions. The model simulated climate and the reanalysis were evaluated against station-based observations of the McArthur Forest Fire Danger Index (FFDI) using probability density function skill scores, annual cumulative FFDI and days per year with FFDI above 50. WRF simulated the main features of the FFDI distribution and its spatial variation, with an overall positive bias. Errors in average FFDI were caused mostly by errors in the ability of WRF to simulate relative humidity. In contrast, errors in extreme FFDI values were driven mainly by WRF errors in wind speed simulation. However, in both cases the quality of the observed data is difficult to ascertain. WRF run with 50-km grid spacing did not consistently improve upon the reanalysis statistics. Decreasing the grid spacing to 10 km led to fire weather that was generally closer to observations than the reanalysis across the full range of evaluation metrics used here. This suggests it is a very useful tool for modelling fire weather over the entire landscape of south-east Australia.

Assimilating GNSS PWV and radiosonde meteorological profiles to improve the PWV and rainfall forecasting performance from the Weather Research and Forecasting (WRF) model over the South China

Coupling of the Weather Research and Forecasting Model with AERMOD for pollutant dispersion modeling. A case study for PM10 dispersion over Pune, India

This research evaluates the quality of high-impact rainfall forecasts across Vietnam and its sub-climate regions. The 3-day rainfall forecast products evaluated include the European Centre for Medium-Range Weather Forecasts (ECMWF) High-Resolution Integrated Forecasting System (IFS) and its downscaled outputs using the Weather Research and Forecasting (WRF) model with the Advanced Research WRF core (WRF-ARW): direct downscaling and downscaling with data assimilation. A full 5-year validation period from 2019 to 2025 was processed. The validation focused on basic rainfall thresholds and also considered the distribution of skill scores for intense events and extreme events. The validations revealed systematic errors (bias) in the models at low rainfall thresholds. The forecast skill was the lowest for northern regions, while the central regions exhibited the highest. For regions strongly affected by terrain, high-resolution downscaling with local observation data assimilation is necessary to improve the detectability of extreme events.

<p>Water vapor (WV) is the most variable greenhouse gas in the atmosphere, which acts as a key feature in climate change studies and plays a crucial role in global warming. Its spatiotemporal distribution is necessary for understanding the hydrological cycle, and consequently can be used as an input factor in climatological studies at global, regional, and local scales. Integrated water vapor (IWV), which is defined as the amount of vertically integrated water vapor, can also augment atmospheric modeling at local and regional scales because it is frequently used in energy budget and evapotranspiration assessments. Currently, there are numerous existing atmospheric models which are able to estimate IWV amount, nevertheless, they fail to obtain extremely accurate results compared with in-situ measurements such as radiosondes. Here, we present a new methodology for improving Weather Research and Forecast (WRF) model predictions accuracy, by using data assimilation technique, which combines estimated 2D IWV regional maps, derived from GPS tropospheric path delays, along with the WRF numerical model output to generate an optimal approximation of the evolving sate of the system. This is done as opposed to pervious works, which assimilated single point measurements, either from radiosondes or GPS zenith wet delay (ZTD) estimation, demonstrating some extant of improvement in the WRF prediction accuracy compare to the standalone WRF numerical runs. Using the suggested technique, our results shows a decrease of up to 30% in the root mean square difference relative to the radiosonde data for WRF predictions assimilated with 2D GPS-IWV regional maps compare to the standalone WRF numerical runs.</p>

The Weather Research and Forecasting (WRF) Model is used to dynamically downscale ERA-Interim global reanalysis data to test its performance as a regional climate model (RCM) for the Great Lakes region (GLR). Four cumulus parameterizations and three spectral nudging techniques applied to moisture are evaluated based on 2-m temperature and precipitation accumulation in the Great Lakes drainage basin (GLDB). Results are compared to a control simulation without spectral nudging, and additional analysis is presented showing the contribution of each nudged variable to temperature, moisture, and precipitation. All but one of the RCM test simulations have a dry precipitation bias in the warm months, and the only simulation with a wet bias also has the least precipitation error. It is found that the inclusion of spectral nudging of temperature dramatically improves a cold-season cold bias, and while the nudging of moisture improves simulated annual and diurnal temperature ranges, its impact on precipitation is complicated. Significance Statement Global climate models are vital to understanding our changing climate. While many include a coarse representation of the Great Lakes, they lack the resolution to represent effects like lake effect precipitation, lake breeze, and surface air temperature modification. Therefore, using a regional climate model to downscale global data is imperative to correctly simulate the land–lake–atmosphere feedbacks that contribute to regional climate. Modeling precipitation is particularly important because it plays a direct role in the Great Lakes’ water cycle. The purpose of this study is to identify the configuration of the Weather Research and Forecasting Model that best simulates precipitation and temperature in the Great Lakes region by testing cumulus parameterizations and methods of nudging the regional model toward the global model.

Due to climate change and associated longer and more frequent droughts, the risk of forest fires increases. To address this, the Institute of Meteorology and Water Management implemented a system for forecasting fire weather in Poland. The Fire Weather Index (FWI) system, developed in Canada, has been adapted to work with meteorological fields derived from the high-resolution (2.5 km) Weather Research and Forecasting (WRF) model. Forecasts are made with 24- and 48-h lead times. The purpose of this work is to present the validation of the implemented system. First, the results of the WRF model were validated using in situ observations from ~70 synoptic stations. Second, we used the correlation method and Eastaugh’s percentile analysis to assess the quality of the FWI index. The data covered the 2019 fire season and were analysed for the whole forest area in Poland. Based on the presented results, it can be concluded that the FWI index (calculated based on the WRF model) has a very high predictive ability of fire risk. However, the results vary by region, distance from human habitats, and size of fire.

The New England region of the northeastern United States has a land use history characterized by forest clearing for agriculture and other uses during European colonization and subsequent reforestation following widespread farm abandonment. Despite these broad changes, the potential influence on local and regional climate has received relatively little attention. This study investigated wintertime (December through March) climate impacts of reforestation in New England using a high-resolution (4 km) multiphysics ensemble of the Weather Research and Forecasting Model. In general, the conversion from mid-1800s cropland/grassland to forest led to warming, but results were sensitive to physics parameterizations. The 2-m maximum temperature (T2max) was most sensitive to choice of land surface model, 2-m minimum temperature (T2min) was sensitive to radiation scheme, and all ensemble members simulated precipitation poorly. Reforestation experiments suggest that conversion of mid-1800s cropland/grassland to present-day forest warmed T2max +0.5 to +3 K, with weaker warming during a warm, dry winter compared to a cold, snowy winter. Warmer T2max over forests was primarily the result of increased absorbed shortwave radiation and increased sensible heat flux compared to cropland/grassland. At night, T2min warmed +0.2 to +1.5 K where deciduous broadleaf forest replaced cropland/grassland, a result of decreased ground heat flux. By contrast, T2min of evergreen needleleaf forest cooled –0.5 to –2.1 K, primarily owing to increased ground heat flux and decreased sensible heat flux.

Ships, sailing in favorable wind or obstructed by wind, will operate with different output power of the engines, and the exhaust emissions will be different even though the ships are sailing at the same ground speed. In this study, the influence of wind was taken into consideration; the ship emission inventory (0.025° × 0.025°) in China of a full year (2014) was reassessed. A speed modification model was employed to figure out the actual output speed of ships by integrating AIS data and the hourly wind. The Weather Research and Forecast (WRF) model was applied to predict the hourly real‐time wind field. The spatial and temporal changes in emissions between the results calculated by the proposed method and our previous study were presented. Overall, when considering the influence of wind, the total ship emissions for the year would increase. In this study, the total estimated emissions of SO2, NOx, PM10, PM2.5, HC, and CO in the area were 1.286 × 106, 2.583 × 106, 2.135× 105, 1.967 × 105, 1.522 × 105, and 3.053 × 105 t(metric tons) in 2014, respectively. Under the influence of wind, the proportion of the regions’ emissions to the total was close to that of the previous study. For SO2 and NOx, emissions presented significant monthly variations. On a monthly timescale, the difference in emissions was more obvious between the results considering and not considering the wind, relative to that on a yearly basis. This study adjusted the method of ship emission estimation, which modified the ship emission inventory over an hourly timescale.

This study is intended to simulate and better understand mean and extreme precipitation over the Loess Plateau (LP) in China using the Weather Research and Forecasting (WRF) model. We performed a long-term (1980–2011) simulation with the WRF model at a 10 km horizontal resolution, forced by the ERA-Interim reanalysis. A series of sensitivity experiments were conducted to investigate the influence of resolution and physical parameterization schemes on the simulation ability of the WRF model over the LP. The region has a dominant semi-arid climate and is adversely affected by extreme precipitation. Results show that the WRF model produced better simulations at a 10 km resolution for the LP with complex terrain than it did at coarser resolutions. WRF simulations for precipitation over the LP were most sensitive to cumulus schemes, moderately sensitive to planetary boundary layer schemes, and less sensitive to microphysics schemes. The WRF model not only adequately captured the spatial distribution of precipitation in the LP but also reproduced its magnitude and variability at different time scales. Moreover, the WRF model captured 73% of the observed extreme precipitation events with an intensity of more than 12 mm day− 1, which cause substantial soil erosion in the LP. Although it was unable to capture 27% of the observed events, the WRF model reasonably reproduced ERA-Interim extreme precipitation and its dynamic processes. These analyses strongly suggest that the biased extreme precipitation generated with the WRF model most likely resulted from inaccurate ERA-Interim forcing data. This study provides a better understanding of extreme precipitation simulations as well as insight into predictions of soil erosion in the LP.

Atmospheric wet delay caused by Precipitable Water Vapor (PWV) significantly impacts the performance of many geodetic surveying systems such as Global Navigation Satellite System (GNSS). In this study, we use wet delay corrections forecast by the Weather Research and Forecasting (WRF) model to enhance GNSS Precise Point Positioning (PPP) during two observation periods with two different weather conditions, that is, period 1: March 01 to 14, 2020 (average PWV: 23.5 kg/m2) and period 2: June 02 to 15, 2020 (flooding weather with average PWV: 55.6 kg/m2), over the South China. PWV data from 277 to 263 GNSS stations are assimilated into WRF model to enhance the WRF water vapor forecasting capability for period 1 and period 2, respectively. Wet delay corrections from two different WRF configurations, that is, WRF no data assimilation and WRF with assimilation of GNSS PWV, are used to augment the PPP. Totally, eight WRF‐enhanced PPP schemes are tested. The results show that WRF‐enhanced PPP schemes generally have a better positioning performance in the up component than traditional PPP. After using WRF wet delay corrections, for static mode, the vertical positioning accuracy is improved by 14.6% and 33.7% for period 1 and period 2, respectively. The corresponding convergence time are reduced by 41.8% and 25.0% for period 1 and period 2, respectively. For kinematic mode, the positioning accuracy improvements in the up component reach 13.8% and 19.0% for period 1 and period 2, respectively. The kinematic PPP convergence time is reduced by up to 8.2% for period 1.