Predicting lightning activity in Greece with the Weather Research and Forecasting (WRF) model

The evaluation and usefulness of lightning prediction for the Indian subcontinent are demonstrated. Implementation of the lightning parameterizations based on storm parameters, in the Weather Research and Forecasting (WRF) Model, with different microphysics schemes are carried out. With the availability of observed lightning measurements over Maharashtra from the lightning detection network (LDN), lightning cases have been identified during the pre-monsoon season of 2016–18. Lightning parameterization based on cloud top height defined by a reflectivity threshold factor of 20 dBZ is chosen. Initial analysis is carried out for 16 lightning events with four microphysical schemes for the usefulness in lightning prediction. Objective analysis is carried out and quantitative model performance (skill scores) is assessed based on observed data. The skills are evaluated for 10- and 50-km2 boxes from the 1-km domain. There is good POD of 0.86, 0.82, 0.85, and 0.84, and false alarm ratio (FAR) of 0.28, 0.25, 0.29, and 0.26 from WSM6, Thompson, Morrison, and WDM6, respectively. There is an overestimation in lightning flash with a spatial and temporal shift. The fractional skill score is evaluated as a function of spatial scale with neighborhoods from 25 to 250 km. These high skill scores and high degree of correlation between observations and model simulation gives us confidence to use the system for real-time operational forecast over India. The skill for 2019 and 2020 pre-monsoon are calculated to address the predictability of operational lightning prediction over India. Significance Statement A high-resolution model, namely, the Weather Research and Forecasting (WRF) Model, with multiple microphysics parameterization schemes and lightning parameterization is used here. The objective analysis is carried out for the lightning cases over India and the quantitative performance is assessed. The results highlight that there is fairly good probability of detection (POD) of 0.86, 0.82, 0.85, and 0.84 and false alarm ratio (FAR) of 0.28, 0.25, 0.29, and 0.26 from four different microphysical schemes (WSM6, Thompson, Morrison, and WDM6, respectively). These high skill scores and high degree of correlation between observations and model simulation gives us confidence to use the system for real-time operational forecast. The validation of lightning forecast system deployed over India for five pre-monsoon months in real time is carried out, which gives POD of 0.90, FAR of 0.64, hit rate of 0.57, and POFD of 0.50 for the whole Indian region.

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.

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.

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

<p>Superbolts are defined as lightning flashes that are a thousand times stronger than normal ones, and their occurrence is estimated to be less than 0.001% of total number of lightning on earth. The global distribution of these extremely powerful lightning flashes is remarkably different than that of regular lightning, which are concentrated in the well-known convective "chimneys" in tropical Africa, South-America and the maritime continent in South-East Asia. The physical mechanisms producing these powerful flashes remain unknown, and the puzzle is exacerbated by the fact that they are discovered mostly over oceans, in maritime winter storms.</p><p>The Mediterranean Sea is one of the most prolific regions where super-bolts occur, especially in the months November-January (Holzworth et al., 2019). We analyzed 8 years of lightning data obtained from the Israeli Lightning Detection Network (ILDN), defining a 200kA peak current threshold for superbolts. We mapped the spatial and temporal distribution of superbolts and their monthly frequency in winter season thunderstorms (DJF) in the eastern Mediterranean, and identified the meteorological and microphysical circumstances in such storms.</p><p>Our working hypothesis is that large amounts of desert dust aerosols, coming from the Sahara Desert, are ingested into maritime winter storms over the eastern Mediterranean. The large dust contributes to convective invigoration, enhanced freezing and efficient charge separation, implying that superbolts are more likely to occur in the presence of large dust. We will present the results of simulation conducted using the WRF-ELEC numerical model, and WRF with spectral bin microphysics coupled with Lynn et al.'s (2012) Dynamic Lightning Scheme (DLS) and the Lightning Potential Index (Yair et al., 2010; LPI), for selected case studies when an enhanced fraction of superbolts was observed.</p><p> </p><p>Holzworth, R. H., McCarthy, M. P., Brundell, J. B., Jacobson, A. R., and Rodger, C. J. (2019). Global distribution of superbolts. J. Geophys. Res. Atmos., 124., doi:10.1029/2019JD030975.</p><p>Lynn, B., Y. Yair, C. Price, G. Kelman and A. J. Clark (2012). Predicting cloud-to-ground and intracloud lightning in weather forecast models. Weather and Forecasting, 27, 1470-1488, doi:10.1175/WAF-D-11-00144.1.</p><p>Yair, Y., B. Lynn, C. Price, V. Kotroni, K. Lagouvardos, E. Morin, A. Mugnai, and M. d. C. Llasat (2010). Predicting the potential for lightning activity in Mediterranean storms based on the Weather Research and Forecasting (WRF) model dynamic and microphysical fields, J. Geophys. Res., 115, D04205, doi:10.1029/2008JD010868</p>

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.

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.

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.

<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>

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.

Medium frequency gravity wave characteristics obtained using Weather Research and Forecasting (WRF) model simulations

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.

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.