Storm Forecasting Research Articles

The Efficient Integration of Dust and Numerical Weather Prediction for Renewable Energy Applications

AbstractThe growing demand for renewable energy underscores the importance of accurate dust forecasting in regions with abundant wind and solar resources. However, leading real‐time global numerical weather prediction (NWP) models often lack dust modules due to computational constraints. Current “Near‐Real‐Time” dust forecasting services can only run after the completion of NWP, failing to meet the timeliness requirements for reporting power generation plans to the grids. This work proposes a global dust‐weather integrated (iDust) model development paradigm, efficiently incorporating dust modules into the dynamical core. Using about one‐eighth additional computing power, iDust extends global 12.5 km resolution NWP with dust prediction capabilities. iDust's forecasting abilities are evaluated against ECMWF CAMS forecast and NASA MERRA2 reanalysis, including verifications over China from March to May 2023 and three extreme dust events. Results show that iDust outperforms its counterparts in dust storm forecasting intensity and timing. Using iDust, global 12.5‐km 10‐day hourly dust storm forecast simulations initiated at 00UTC can produce results by 06UTC, enabling timely forecasting of severe dust storms with concentrations exceeding 1,000 μg/m3. This novel capability of iDust can meet the urgent forecasting needs of the renewable energy industry for extreme dust conditions, supporting the green energy transition.

Correspondence of variations of <i>AE</i> and <i>Apo</i> indices in 23–24 solar cycles

The auroral electrojet index AE is often used in forecasting models as a characteristic of a source of the disturbance propagation in the geosphere from the pole to middle and low latitudes. However, these data are no longer available digitally since January 2020. Instead of the AE−index, we suggest using the recently introduced 1 h Apo−index, given the close proximity of magnetometer networks for these indices at high latitudes and the availability of the Apo−index in real time. To this end their correlation is analyzed during 276 intense storms for 1995–2017. Storm profiles are constructed by method of superposed epoch with zero epoch time t0 = 0 taken at the threshold value of AE ≥ 1000 nT. A comparison is made of the storm profiles of AE(t), Apo(t), the interplanetary electric field E(t) and the solar wind speed Vsw(t) within 72 hours: 24 hours before the storm peak t0, and 48 hours after it. A good agreement is obtained between the sets AE(t) and Apo(t) with a correlation coefficient of 0.70. Comparison with the interplanetary parameters testifies on the correlation of AE(t) and Apo(t) with the electric field E(t) but absence of their coupling with the solar wind speed Vsw(t). A two−parametric formula is derived for dependence of the auroral electrojet index AE(t) on the interplanetary electric field E(t) and the geomagnetic Apo(t) index for the geomagnetic storm forecasting. In the absence of E(t) data, formulae for the dependence of AE(t) on Apo(t) is introduced for implementation in real time and the inverse dependence of Apo(t) on AE(t) for reconstruction of the 1 h Apo−index before 1995. Validation of the proposed models with data for 5 intense storms in 2018 has shown a close resemblance of the model with observation data of the AE−index with a high coefficient of determination R2 ranging from 0.62 to 0.81.

Assimilation of INSAT-3D Imager Retrieved Aerosol Optical Depth for Enhanced Dust Storm Prediction Using WRF-Chem Model

Abstract Frequent dust storms during the pre-monsoon season in northwest India significantly impact weather, air quality, and health, necessitating accurate predictions. This study demonstrates the first-time assimilation of aerosol optical depth (AOD) data from INSAT-3D into the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem), alongside previously assimilated AOD data from MODIS AQUA. Focusing on an extreme dust storm over the Indo-Gangetic Plain (IGP) from June 12-17, 2018, identical experiments were conducted with and without AOD assimilation. The Gridpoint Statistical Interpolation (GSI) three-dimensional variational (3DVAR) method was utilized for data assimilation, with simulations performed for 144 hours, adjusting for uncertainties in the model's background error covariance using different inflation factors. The model's performance was validated using the ECMWF Atmospheric Composition Reanalysis 4 (EAC4) product. The assimilated AOD analysis significantly improved bias, root mean square difference (RMSD), and correlation from [0.47, 0.71, 0.40] to [0.07, 0.38, 0.80]. The 24-hour AOD forecast improved by approximately 56% and 30% with INSAT-3D and MODIS AOD assimilation, respectively. AOD assimilation directly impacted cloud properties and radiative forcing, enhancing the 24-hour forecast of net shortwave downward flux (SWD) by approximately 51 W/m² and 24 W/m² for INSAT-3D and MODIS, respectively. These changes also affected vertical levels up to 5 km, modifying the rainwater mixing ratio (QRAIN) and surface rainfall forecasts. This study shows that integrating INSAT-3D AOD data provides substantial improvements in aerosol forecasts, crucial for making reliable and accurate forecasts of severe dust storms and air pollution episodes over the Indian region.

Source contributions to two super dust storms over Northern China in March 2021 and the impact of soil moisture

Two extremely devastating super dust storms (SDS) hit Mongolia and Northern China in March 2021, causing many deaths and substantial economic damage. Accurate forecasting of dust storms is of great importance for avoiding or mitigating their effects. One of the most critical factors affecting dust emissions is soil moisture, but its value in desert exhibits significant uncertainty. In this study, model experiments were conducted to simulate dust emissions using four soil moisture datasets. The results were compared with observations to assess the effects of soil moisture on the dust emission strength. The Integrated Source Apportionment Method (ISAM) was used to track the dust sources and quantify the contribution from each source region to the dust load over the North China Plain (NCP), Korea peninsula, and western Japan. The results show large differences in the dust load depending on the soil moisture datasets used. The high soil moisture in the NCEP dataset results in substantial underestimation of the dust emission flux and PM10 concentration. Despite a minor overestimation of PM10 concentrations in many Northern China cities, the ERA5 dataset yields the best simulation performance. During the two SDS events, about 7.5 Mt dust was released from the deserts in Mongolia and 2.8 Mt from the deserts in China. Source apportionment indicates that the Mongolian Gobi Desert is the dominant source of PM10 in the NCP, Korea peninsula, and western Japan, accounting for 60 %–80 %, while Inner Mongolia contributed 10 %–20 %.

Winter Storm Severity Index in Alaska: Understanding the Usefulness for Impact-Based Winter Weather Severity Forecast Information

Abstract There is growing interest in impact-based decision support services to address complex decision-making, especially for winter storm forecasting. Understanding users’ needs for winter storm forecast information is necessary to make such impact-based winter forecasts relevant and useful to the diverse regions affected. A mixed-method social science research study investigated extending the winter storm severity index (WSSI) [operational for the contiguous United States (CONUS)] to Alaska, with consideration of the distinct needs of Alaskan stakeholders and the Alaskan climate. Data availability differences suggest the need for an Alaska-specific WSSI, calling for user feedback to inform the direction of product modifications. Focus groups and surveys in six regions of Alaska provided information on how the WSSI components, definitions, and categorization of impacts could align with stakeholder expectations and led to recommendations for the Weather Prediction Center to consider in developing the WSSI Alaska product. Overall, wind (strength and direction) and precipitation are key components to include. Air travel is a critical concern requiring wind and visibility information, while road travel is less emphasized (contrasting with CONUS needs). Special Weather Statements and Winter Storm Warnings are highly valued, and storm trajectory and transition (between precipitation types) information are the important contexts for decision-makers. Alaska is accustomed to and prepared for winter impacts but being able to understand how components (wind, snow, and ice) contribute to overall impact enhances the ability to respond and mitigate damage effectively. The WSSI adapted for Alaska can help address regional forecast needs, particularly valuable as the climate changes and typical winter conditions become more variable. Significance Statement Impact-based support services can assist decision-makers in prioritizing preparedness and mitigation actions related to winter storm events. The winter storm severity index adapted for specific considerations in Alaska (such as including wind and visibility components) can extend winter weather impact-based forecasting’s utility. Additionally, lessons learned from the process of adapting a national product to specific regional needs may inform best practices for gathering stakeholder input and feedback.

Understanding winter windstorm predictability over Europe

Abstract. Winter windstorms belong to the most damaging meteorological events in the extra-tropics. Their impact on society makes it essential to understand and improve seasonal forecasts of these extreme events. Skilful predictions on a seasonal timescale have been shown in previous studies by investigating hindcasts from various forecast centres. This study aims to explain storm forecast skill based on relevant dynamical factors. Therefore, a number of factors which are known to influence either windstorms directly or their synoptic relevant systems, mid-latitude cyclones, are investigated. These factors are analysed for their relation to windstorm forecast performance based on a reanalysis (ERA5) and the seasonal hindcast of the UK Met Office (Global Seasonal forecasting system version 5, GloSea5). Within GloSea5, relevant dynamical factors are (1) validated with respect to their physical connections to windstorms, (2) investigated with respect to the seasonal forecast skill of the factors themselves, and (3) assessed on the relevance and influence of their forecast performance to and on windstorm forecast skill. Although not all investigated factors reveal a clear and consistent influence on windstorm forecast skill over Europe, core factors like mean sea level pressure gradient, sea surface temperature, equivalent potential temperature and Eady growth rate show consistent results within these three steps: their physical connection is well represented in the model; these factors are skilfully predicted in storm-relevant regions, and, consequently, this skill leads to increased forecast skill of winter windstorms over Europe. This study thus explains existing forecast skill in winter windstorms but also indicates potential for further model developments to improve seasonal winter windstorm predictions.

Do AI models produce better weather forecasts than physics-based models? A quantitative evaluation case study of Storm Ciarán

There has been huge recent interest in the potential of making operational weather forecasts using machine learning techniques. As they become a part of the weather forecasting toolbox, there is a pressing need to understand how well current machine learning models can simulate high-impact weather events. We compare short to medium-range forecasts of Storm Ciarán, a European windstorm that caused sixteen deaths and extensive damage in Northern Europe, made by machine learning and numerical weather prediction models. The four machine learning models considered (FourCastNet, Pangu-Weather, GraphCast and FourCastNet-v2) produce forecasts that accurately capture the synoptic-scale structure of the cyclone including the position of the cloud head, shape of the warm sector and location of the warm conveyor belt jet, and the large-scale dynamical drivers important for the rapid storm development such as the position of the storm relative to the upper-level jet exit. However, their ability to resolve the more detailed structures important for issuing weather warnings is more mixed. All of the machine learning models underestimate the peak amplitude of winds associated with the storm, only some machine learning models resolve the warm core seclusion and none of the machine learning models capture the sharp bent-back warm frontal gradient. Our study shows there is a great deal about the performance and properties of machine learning weather forecasts that can be derived from case studies of high-impact weather events such as Storm Ciarán.

Improvement in Storm Surge Numerical Forecasting Based on Wave Buoys Data

The maximum wind speed radius of a strong typhoon making landfall is an important factor influencing the numerical forecasting of storm surges. A method for inverting the maximum wind speed radius of typhoons based on wave buoys data was designed to significantly reduce the error in 24 h storm surge forecasting in this paper, and an operation scheme was proposed to enhance the storm surge numerical forecasting system based on this method. Hangzhou Bay and the Yangtze River Estuary, which have been frequently impacted by typhoons over the past five years, were selected as the research area. Common schemes for the maximum wind speed radius were analyzed, and five ladder schemes (10, 15, 20, 25, and 30 km) were established for wave and storm surge numerical model verification of Typhoon Muifa in 2022. Based on a comparison of the wave hindcast results and wave buoys observation data, the wave hindcast result of the commonly used scheme (30 km) was significantly greater than that of the observation data, and the optimal scheme (15 km) closest to the observation data could be determined during the 48 h warning period. Moreover, it was difficult to identify the optimal scheme during the 48 h warning period based on the storm surge hindcast results. A 24 h storm surge numerical forecasting test was performed with the commonly used scheme (30 km) and the optimal scheme (15 km). The results showed that the root mean square error (RMSE) of the optimal scheme (15 km) was 34% lower than that of the commonly used scheme (30 km), while the maximum storm surge error was also reduced from 47.7% for the commonly adopted scheme (30 km) to 11.8% for the optimal scheme (15 km). The maximum storm surges under the optimal scheme (15 km) along Hangzhou Bay and the Yangtze River Estuary ranged from 1.9 to 2.2 m, which were closer to the observation data, and the maximum storm surge under the commonly used scheme (30 km) was 0.8~1.2 m greater than that under the optimal scheme (15 km).

Improving the Asian dust storm prediction using WRF-Chem through combinational optimization of physical parameterization schemes

This study aims to enhance the accuracy of the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem) in forecasting Asian dust storms (ADSs) by using the micro-Genetic Algorithm (μGA). We developed an optimization system---the WRF-Chem-μGA system---to seek the optimal combination of the planetary boundary layer (PBL) and land surface parameterization schemes, which are crucial for numerical forecast of dust storms. The optimization was conducted concerning meteorological and air quality variables, i.e., aerosol optical depth, PBL height, 2 m temperature, 2 m relative humidity, and 10 m wind speed, simultaneously for three ADS cases over the optimization domain, including South Korea. Among a total of 32 available combinations of physical parameterization scheme options (8 from PBL and 4 from land surface schemes), the optimized set through the WRF-Chem-μGA system consists of the Asymmetrical Convective Model version 2 (ACM2) for the PBL scheme and the Noah land surface model with Multiple Parameterization options (Noah-MP) for the land surface scheme. The optimized set showed an improvement ratio of up to 22.5 % in terms of the normalized RMSE for all meteorological and air quality variables, compared to various non-optimized sets of physical parameterization schemes for two additional ADS cases. The optimal set proposed in this study can be used comprehensively in numerical forecasts of various meteorological and air quality problems in the East Asian region, using the WRF-Chem model.

Stratospheric influence on the winter North Atlantic storm track in subseasonal reforecasts

Abstract. Extreme stratospheric polar vortex events, such as sudden stratospheric warmings (SSWs) or extremely strong polar vortex events, can have a significant impact on surface weather in winter. SSWs are most often associated with negative North Atlantic Oscillation (NAO) conditions, cold air outbreaks in the Arctic and a southward-shifted midlatitude storm track in the North Atlantic, while strong polar vortex events tend to be followed by a positive phase of the NAO, relatively warm conditions in the extratropics and a poleward-shifted storm track. Such changes in the storm track position and associated extratropical cyclone frequency over the North Atlantic and Europe can increase the risk of extreme windstorm, flooding or heavy snowfall over populated regions. Skillful predictions of the downward impact of stratospheric polar vortex extremes can therefore improve the predictability of extratropical winter storms on subseasonal timescales. However, there exists a strong inter-event variability in these downward impacts on the tropospheric storm track. Using ECMWF reanalysis data and reforecasts from the Subseasonal to Seasonal (S2S) Prediction Project database, we investigate the stratospheric influence on extratropical cyclones, identified with a cyclone detection algorithm. Following SSWs, there is an equatorward shift in cyclone frequency over the North Atlantic and Europe in reforecasts, and the opposite response is observed after strong polar vortex events, consistent with the response in reanalysis. However, although the response of cyclone frequency following SSWs with a canonical surface impact is typically captured well during weeks 1–4, less than 25 % of the reforecasts manage to capture the response following SSWs with a “non-canonical” impact. This suggests a possible overconfidence in the reforecasts with respect to reanalysis in predicting the canonical response after SSWs, although it only occurs in about two-thirds of the events. The cyclone forecasts following strong polar vortex events are generally more successful. Understanding the role of the stratosphere in subseasonal variability and predictability of storm tracks during winter can provide a key for reliable forecasts of midlatitude storms and their surface impacts.

Testing Stochastic and Perturbed Parameter Methods in an Experimental 1-km Warn-on-Forecast System Using NSSL’s Phased-Array Radar Observations

Abstract The success of the National Severe Storms Laboratory’s (NSSL) experimental Warn-on-Forecast System (WoFS) to provide useful probabilistic guidance of severe and hazardous weather is mostly due to the frequent assimilation of observations, especially radar observations. Phased-array radar (PAR) technology, which is a potential candidate to replace the current U.S. operational radar network, would allow for even more rapid assimilation of radar observations by providing full-volumetric scans of the atmosphere every ∼1 min. Based on previous studies, more frequent PAR data assimilation can lead to improved forecasts, but it can also lead to ensemble underdispersion and suboptimal observation assimilation. The use of stochastic and perturbed parameter methods to increase ensemble spread is a potential solution to this problem. In this study, four stochastic and perturbed parameter methods are assessed using a 1-km-scale version of the WoFS and include the stochastic kinetic energy backscatter (SKEB) scheme, the physically based stochastic perturbation (PSP) scheme, a fixed perturbed parameters (FPP) method, and a novel surface-model scheme blending (SMSB) method. Using NSSL PAR observations from the 9 May 2016 tornado outbreak, experiments are conducted to assess the impact of the methods individually, in different combinations, and with different cycling intervals. The results from these experiments reveal the potential benefits of stochastic and perturbed parameter methods for future versions of the WoFS. Stochastic and perturbed parameter methods can lead to more skillful forecasts during periods of storm development. Moreover, a combination of multiple methods can result in more skillful forecasts than using a single method. Significance Statement Phased-array radar technology allows for more frequent assimilation of radar observations into ensemble forecast systems like the experimental Warn-on-Forecast System. However, more frequent radar data assimilation can eventually cause issues for prediction systems due to the lack of ensemble spread. Thus, the purpose of this study is to explore the use of four stochastic and perturbed parameter methods in a next-generation Warn-on-Forecast System to generate ensemble spread and help prevent the issues from frequent radar data assimilation. Results from this study indicate the stochastic and perturbed parameter methods can improve forecasts of storms, especially during storm development.

Kombinace indexů stability pro předpovědi bouřek

Large amount of both historical and modern indices calculated from radio-souding data was developed for the storm prediction. In this work we processed 75 storm indices and other parameters, which could potentially be connected with the storm development. First part of the article is devoted to the effectivity of particular indices with special emphasis on the comparison of various modifications of the convective available potential energy (CAPE). Afterwards there were examined possibilities of the use of various combinations of indices and other parameters, and benefits of these combinations for the thunderstorm forecast. The forecast of the afternoon storms based on midday soundings is analysed in detail, as for this time point, a combination of indices has apparently the biggest benefit.

A Case Study Investigating the Low Summertime CAPE Behavior in the Global Forecast System

Abstract Convective available potential energy (CAPE) is an important index for storm forecasting. Recent versions (v15.2 and v16) of the Global Forecast System (GFS) predict lower values of CAPE during summertime in the continental United States than analysis and observation. We conducted an evaluation of the GFS in simulating summertime CAPE using an example from the Unified Forecast System Case Study collection to investigate the factors that lead to the low CAPE bias in GFS. Specifically, we investigated the surface energy budget, soil properties, and near-surface and upper-level meteorological fields. Results show that the GFS simulates smaller surface latent heat flux and larger surface sensible heat flux than the observations. This can be attributed to the slightly drier-than-observed soil moisture in the GFS that comes from an offline global land data assimilation system. The lower simulated CAPE in GFS v16 is related to the early drop of surface net radiation with excessive boundary layer cloud after midday when compared with GFS v15.2. A moisture-budget analysis indicates that errors in the large-scale advection of water vapor does not contribute to the dry bias in the GFS at low levels. Common Community Physics Package single-column model (SCM) experiments suggest that with realistic initial vertical profiles, SCM simulations generate a larger CAPE than runs with GFS IC. SCM runs with an active LSM tend to produce smaller CAPE than that with prescribed surface fluxes. Note that the findings are only applicable to this case study. Including more warm-season cases would enhance the generalizability of our findings. Significance Statement Convective available potential energy (CAPE) is one of the key parameters for severe weather analysis. The low bias of CAPE is identified by forecasters as one of the key issues for the NOAA operational global numerical weather prediction model, Global Forecast System (GFS). Our case study shows that the lower CAPE in GFS is related to the drier atmosphere than observed within the lowest 1 km. Further investigations suggest that it is related to the drier atmosphere that already exists in the initial conditions, which are produced by the Global Data Assimilation System, in which an earlier 6-h GFS forecast is combined with current observations. It is also attributed to the slightly lower simulated soil moisture than observed. The lower CAPE in GFS v16 when compared with GFS v15.2 in the case analyzed here is related to excessive boundary layer cloud formation beginning at midday that leads to a drop of net radiation reaching the surface and thus less latent heat feeding back to the low-level atmosphere.

Delayed-mode reprocessing of in situ sea level data for the Copernicus Marine Service

Abstract. The number of tide gauges providing coastal sea level data has significantly increased in recent decades. They help in the issue of coastal hazard warnings, in the forecasting (indirectly through models) of storm surges and tsunamis, and in operational oceanography applications. These data are automatically quality controlled in near-real time in the Copernicus Marine Service. A new initiative seeks to provide delayed-mode reprocessed data for the Copernicus Marine Service by developing a new product and upgrading the software used in its automated quality control. Several new modules, such as buddy checking or the detection of attenuated data, are implemented. The new product was launched in November of 2022. The entire reprocessing is discussed in detail. An example of the information that can be extracted from the delayed-mode reprocessed product is also given.

Forecast of Intense Geomagnetic Storms in February-March 2023, Based on Measurements of Galactic Cosmic Ray Intensity Fluctuations

Forecast of Intense Geomagnetic Storms in February-March 2023, Based on Measurements of Galactic Cosmic Ray Intensity Fluctuations

Regional and Seasonal Biases in Convection-Allowing Model Forecasts of Near-Surface Temperature and Moisture

Abstract This study investigates regional, seasonal biases in convection-allowing model forecasts of near-surface temperature and dewpoint in areas of particular importance to forecasts of severe local storms. One method compares model forecasts with objective analyses of observed conditions in the inflow sectors of reported tornadoes. A second method captures a broader sample of environments, comparing model forecasts with surface observations under certain warm-sector criteria. Both methods reveal a cold bias across all models tested in Southeast U.S. cool-season warm sectors. This is an operationally important bias given the thermodynamic sensitivity of instability-limited severe weather that is common in the Southeast cool season. There is not a clear bias across models in the Great Plains warm season, but instead more varied behavior with differing model physics. Significance Statement The severity of thunderstorms and the types of hazards they produce depend in part on the low-level temperature and moisture in the near-storm environment. It is important for numerical forecast models to accurately represent these fields in forecasts of severe weather events. We show that the most widely used short-term, high-resolution forecast models have a consistent cold bias of about 1 K (up to 2 K in certain cases) in storm environments in the southeastern U.S. cool season. Human forecasters must recognize and adjust for this bias, and future model development should aim to improve it.

Investigating the Storm Surge and Flooding in Shenzhen City, China

Tropical cyclones affecting Shenzhen city have shown a remarkable tendency to increase in both intensity and quantity, highlighting the urgency of accurate forecasts of storm surges and flooding for effective planning and mitigation. Utilizing satellite and field observations together with the advanced high-resolution baroclinic wave–current model (SCHISM), a comprehensive investigation aimed at storm surge and flooding in Shenzhen was conducted. Statistical work of historical tropical cyclones revealed that Shenzhen was most vulnerable to cyclones propagating from the southeast toward the northwest and passing Shenzhen down the Pearl River Estuary. Thus, a representative, i.e., super typhoon Hato (2017), was selected for further study. Validations of numerical results suggested satisfactory model performance in mapping the wave, tide, and surge processes. Remarkable differences in spatiotemporal distribution and intensity of storm surge and flooding were found along the Shenzhen coast, which was dominated by the propagation of far-field surge and tidal waves, cooperation between wind direction and coastline orientation, estuary morphology, and the land terrain. Intervention of wave–current interaction improved the simulation of the surge and flooding and triggered an earlier occurrence time of the maximum surge in specific areas. The Pearl River discharge significantly elevated the sea level height inside the estuary and contributed to a more severe surge. Given the extremely complicated river networks and huge freshwater flux of Pearl River and the increasing trend of concurrent heavy precipitation of tropical cyclones, future investigations on compound flooding were suggested.

Classifying and bounding geomagnetic storms based on the SYM-H and ASY-H indices

Geomagnetic storms, which are significant disturbances in Earth’s magnetic field, pose a natural hazard with potentially harmful effects on technological and electrical systems. Various indices such as Kp, Dst, and the newer SYM-H and ASY-H are used to measure storm intensity. However, there is no consensus among researchers regarding storm classification based on intensity and asymmetry for the SYM-H and ASY-H indices, as existing classifications mainly rely on the Dst and Kp indices. This study proposes a classification method based on the cumulative distribution function of the SYM-H and ASY-H indices, applying industry-wide percentiles to determine storm classes. Using percentiles and a superposed epoch analysis, storms are objectively classified based on intensity and occurrence probability. This classification approach has applications in storm forecasting and risk assessment.

Using Object-Based Verification to Assess Improvements in Forecasts of Convective Storms between Operational HRRR Versions 3 and 4

Abstract The object-based verification procedure described in a recent paper by Duda and Turner was expanded herein to compare forecasts of composite reflectivity and 6-h precipitation objects between the two most recent operational versions of the High-Resolution Rapid Refresh (HRRR) model, versions 3 and 4, over an expanded set of warm season cases in 2019 and 2020. In addition to analyzing all objects, a reduced set of forecast–observation object pairs was constructed by taking the best forecast match to a given observation object for the purposes of bias-reduction and unequivocal object comparison. Despite the apparent signal of improved scalar metrics such as the object-based threat score in HRRRv4 compared to HRRRv3, no statistically significant differences were found between the models. Nonetheless, many object attribute comparisons revealed indications of improved forecast performance in HRRRv4 compared to HRRRv3. For example, HRRRv4 had a reduced overforecasting bias for medium- and large-sized reflectivity objects, and all objects during the afternoon. HRRRv4 also better replicated the distribution of object complexity and aspect ratio. Results for 6-h precipitation also suggested superior performance in HRRRv4 over HRRRv3. However, HRRRv4 was worse with centroid displacement errors and more severely overforecast objects with a high maximum precipitation amount. Overall, this exercise revealed multiple forecast deficiencies in the HRRR, which enables developers to direct development efforts on detailed and specific endeavors to improve model forecasts. Significance Statement This work builds upon the authors’ prior work in assessing model forecast quality using an alternative verification method—object-based verification. In this paper we verified two versions of the same model (one an upgrade from the other) that were making forecasts covering the same time window, using the object-based verification method. We found that the updated model was not statistically significantly better, although there were indications it performed better in certain aspects such as capturing the change in the number of storms during the daytime. We were able to identify specific problem areas in the models, which helps us direct model developers in their efforts to further improve the model.

An Investigation on Joint Data Assimilation of a Radar Network and Ground-Based Profiling Platforms for Forecasting Convective Storms

Abstract A summer convective precipitation case, occurring in eastern China on 16–17 July 2020, is selected to investigate the impact of joint assimilation of ground-based profiling platforms and weather radars on forecasting convective storms using observational system simulation experiments (OSSEs). The simulated profiling platforms include the Doppler wind lidar (DWL), a wind profiler (WP), and a microwave radiometer (MWR). Results show that joint assimilation of WP and radar data produces a better analysis of convective dynamical structure than joint assimilation of DWL and radar data, since WP detects deeper layer winds. Joint assimilation of MWR and radar data enables rapid adjustment of temperature and humidity and thus, avoids the potential errors introduced by the latent heat term of the radar diabatic initialization in the early stage. Profiling observations in a horizontal spacing of 80 km provide fewer benefits for convective forecasting, while reducing the spacing to 40 km can dramatically improve model analysis and forecasts. Joint assimilation of multiple profiling observations in a 20-km horizontal spacing with radar data exhibits a beneficial synergistic effect and mitigates “the ramp-down issue” during the forecast stage. Assimilating profiling observations with an update interval less than 30 min does not have as pronounced an effect on convective forecasts as horizontal spacing. Furthermore, assimilating profiling observations at a 20-km horizontal spacing can obtain accurate mesoscale background environment and forecast storms with an ability comparable to radar data assimilation. This work emphasizes the need to consider implementing a joint mesoscale detection system that incorporates weather radars and profiling observations for leveraging convective storm forecasting.