Radar Reflectivity–Based Model Initialization Using Specified Latent Heating (Radar-LHI) within a Diabatic Digital Filter or Pre-Forecast Integration

The High-Resolution Rapid Refresh (HRRR) is a convection-allowing implementation of the Advanced Research version of the Weather Research and Forecasting (WRF-ARW) Model with hourly data assimilation that covers the conterminous United States and Alaska and runs in real time at the NOAA/National Centers for Environmental Prediction (NCEP). Implemented operationally at NOAA/NCEP in 2014, the HRRR features 3-km horizontal grid spacing and frequent forecasts (hourly for CONUS and 3-hourly for Alaska). HRRR initialization is designed for optimal short-range forecast skill with a particular focus on the evolution of precipitating systems. Key components of the initialization are radar-reflectivity data assimilation, hybrid ensemble-variational assimilation of conventional weather observations, and a cloud analysis to initialize stratiform cloud layers. From this initial state, HRRR forecasts are produced out to 18 h every hour, and out to 48 h every 6 h, with boundary conditions provided by the Rapid Refresh system. Between 2014 and 2020, HRRR development was focused on reducing model bias errors and improving forecast realism and accuracy. Improved representation of the planetary boundary layer, subgrid-scale clouds, and land surface contributed extensively to overall HRRR improvements. The final version of the HRRR (HRRRv4), implemented in late 2020, also features hybrid data assimilation using flow-dependent covariances from a 3-km, 36-member ensemble (“HRRRDAS”) with explicit convective storms. HRRRv4 also includes prediction of wildfire smoke plumes. The HRRR provides a baseline capability for evaluating NOAA’s next-generation Rapid Refresh Forecast System, now under development. Significance Statement NOAA’s operational hourly updating, convection-allowing model, the High-Resolution Rapid Refresh (HRRR), is a key tool for short-range weather forecasting and situational awareness. Improvements in assimilation of weather observations, as well as in physics parameterizations, have led to improvements in simulated radar reflectivity and quantitative precipitation forecasts since the initial implementation of HRRR in September 2014. Other targeted development has focused on improved representation of the diurnal cycle of the planetary boundary layer, resulting in improved near-surface temperature and humidity forecasts. Additional physics and data assimilation changes have led to improved treatment of the development and erosion of low-level clouds, including subgrid-scale clouds. The final version of HRRR features storm-scale ensemble data assimilation and explicit prediction of wildfire smoke plumes.

In this study, we evaluate the ability of the High-Resolution Rapid Refresh (HRRR) model to forecast cloud characteristics through comparison of observed and simulated satellite brightness temperatures (BTs) and radar reflectivity during different weather phenomena in December 2021: the Mayfield, Kentucky, tornado on 11 December, a heavy snow event in Minnesota from 10 to 11 December, and the Midwest derecho on 15 December. This is done to illustrate the importance of examining model accuracy across a range of weather phenomena. Observation and forecast objects were created using the Method for Object-Based Diagnostic Evaluation (MODE). HRRR accurately depicted the spatial displacements between observation cloud (defined using BTs) and radar reflectivity objects, namely, the centers of cloud objects are to the east of the radar objects for the tornado and derecho events, and generally west of the radar objects for the snow event. However, HRRR had higher (less intense) simulated BTs and higher (more intense) radar reflectivity than the observations for the tornado event. Simulated radar reflectivity is higher and BTs are lower than the observations during the middle of the snow event. Also, simulated radar reflectivity is higher and BTs are lower than the observations during the derecho event. Of the three weather events, the HRRR forecasts are most accurate for the snow event, based on the object-based threat score, followed by the derecho and tornado events. The tornado event has lower accuracy because matches between paired simulated and observation objects are worse than for the snow event, with less similarity in size forecast objects and greater distance between paired object centers. Significance Statement The purpose of this study is to assess the accuracy of forecast cloud and radar objects, defined using simulated satellite brightness temperatures and radar reflectivity, from the High-Resolution Rapid Refresh (HRRR) model. This assessment was conducted for a tornado, snow, and derecho event from December 2021. Results from these three events indicate that the HRRR model accurately represents the observed displacement between the center of cloud and radar objects for the tornado and derecho events, and is the most accurate overall for the snow event.

Strong wind events can cause severe economic loss and societal impacts. Peak winds over Wyoming and Colorado occur during the wintertime months, and in the first part of this two-part series, it was shown that the High-Resolution Rapid Refresh (HRRR) model displays large negative biases with respect to strong wind events. In this part of the study, we address two questions: 1) does increasing the horizontal resolution improve the representation of strong wind events over this region, and 2) are the biases in HRRR-forecasted winds related to the selected planetary boundary layer (PBL), surface layer (SL), and/or land surface model (LSM) parameterizations? We conduct Weather Research and Forecasting (WRF) Model simulations to address these two main questions. Increasing the horizontal resolution leads to a small improvement in the simulation of strong wind speeds over the complex terrain of Wyoming and Colorado. In general, changes in the PBL, SL, and LSM parameterizations show much larger changes in simulated wind speeds compared to increasing the model resolution alone. Specifically, changing from the Mellor–Yamada–Nakanishi–Niino scheme to the Mellor–Yamada–Janjić PBL and SL schemes results in nearly no change in the r2 values, but there is a decrease in the magnitude of the strong wind speed bias (from −12.52 to −10.16 m s−1). We attribute these differences to differences in the diagnosis of surface wind speeds and mixing in the boundary layer. Further analysis is conducted to determine the value of 1-km forecasts of strong winds compared with wind speed diagnostics commonly used by the National Weather Service. Significance Statement Motivated by prior studies showing low skill in the prediction of strong winds using state-of-the-art weather forecast models, in this study, we aim to investigate two questions: 1) does increasing the horizontal resolution improve the prediction of strong wind events over the complex terrain of Wyoming and Colorado, and 2) are the biases in High-Resolution Rapid Refresh (HRRR) forecasted winds related to the selected planetary boundary layer, surface layer, and/or land surface model parameterizations? We find that increasing the horizontal resolution provides a slight improvement in the prediction of strong winds. Further, considerable improvement in the prediction of strong winds is found for varying boundary layer, surface layer, and land surface parameterizations.

Traditional pixel-versus-pixel forecast evaluation scores such as the critical success index (CSI) provide a simple way to compare the performances of different forecasts; however, they offer little information on how to improve a particular forecast. This paper strives to demonstrate what additional information an object-based forecast evaluation tool such as the Method for Object-Based Diagnostic Evaluation (MODE) can provide in terms of assessing numerical weather prediction models’ convective storm forecasts. Forecast storm attributes evaluated by MODE in this paper include storm size, intensity, orientation, aspect ratio, complexity, and number of storms. Three weeks of the High Resolution Rapid Refresh (HRRR) model’s precipitation forecasts during the summer of 2010 over the eastern two-thirds of the contiguous United States were evaluated as an example to demonstrate the methodology. It is found that the HRRR model was able to forecast convective storm characteristics rather well either as a function of time of day or as a function of storm size, although significant bias does exist, especially in terms of storm number and storm size. Another interesting finding is that the model’s ability of forecasting new storm initiation varies substantially by regions, probably as a result of its different skills in forecasting convection driven by different forcing mechanisms (i.e., diurnal heating vs synoptic-scale frontal systems).

Latent heating (LH) is an important quantity in both weather forecasting and climate analysis, being the essential factor driving convective systems. Yet, inferring LH rates from our current observing systems is challenging at best. For climate studies, LH has been retrieved from the Precipitation Radar on the Tropical Rainfall Measuring Mission (TRMM) using model simulations in the look-up table (LUT) that relates instantaneous radar profiles to corresponding heating profiles. These radars, first on TRMM and then Global Precipitation Measurement Mission (GPM), provide a continuous record of LH. However, temporal resolution is too coarse to have a significant impacts on forecast models. In operational forecast models such as High-Resolution Rapid Refresh, convection is initiated from LH derived from ground based radar. Despite the high spatial and temporal resolution of ground-based radars, one disadvantage of using these sources is that its data are only available over well observed land areas. This study develops a method to derive LH from the Geostationary Operational Environmental Satellite-16 (GOES-16) in near-real time. Even though the visible and infrared channels on the Advanced Baseline Imager (ABI) provide mostly cloud top information, rapid changes in cloud top visible and infrared properties, when formulated as a LUT similar to those used by the TRMM and GPM radars, can equally be used to derive LH profiles for convective regions based on model simulations with a convective classification scheme and channel 14 (11.2 μm) brightness temperatures. Convective regions detected by GOES-16 are assigned LH from the LUT, and they are compared with LH from the Next Generation Weather Radar (NEXRAD) and one of the Dual-frequency Precipitation Radar (DPR) products, the Goddard Convective-Stratiform Heating (CSH). LH obtained from GOES-16 show similar magnitude with NEXRAD and CSH, and vertical distribution of LH is also very similar with CSH. One month analysis of total LH from convective clouds from GOES-16 and NEXRAD shows good correlation between the two products. Finally LH profiles from GOES-16 and NEXRAD are applied to WRF simulations for convective initiation and their results are compared to investigate their impacts in precipitation forecasts. Results show that LH from GOES-16 have similar impacts as NEXRAD, and improves the forecast significantly.

The hourly updating High-Resolution Rapid Refresh (HRRR) model is evaluated with regard to its ability to predict the areal extent of cold-season precipitation and accurately depict the timing and location of regions of snow, rain, and mixed-phase precipitation on the ground. Validation of the HRRR forecasts is performed using observations collected by the Automated Surface Observing System (ASOS) stations across the eastern two-thirds of the United States during the 2010–11 cold season. The results show that the HRRR is able to reliably forecast precipitation extent during the cold season. In particular, the location and areal extent of both snow and rain are very well predicted. Depiction of rain-to-snow transitions and freezing rain is reasonably good; however, the associated evaluation scores are significantly lower than for either snow or rain. The analyses suggest the skill in accurately depicting precipitation extent and phase (i.e., rain, snow, and mixed phase) depends on the size and organization of a weather system. Typically, larger synoptically forced weather systems are better predicted than smaller weather systems, including the associated rain-to-snow transition or freezing-rain areas. Offsets in space or time (i.e., causing misses and false alarms) have a larger effect on the model performance for smaller weather systems.

Tropical cyclone tornadoes (TCTORs) are a hazard to life and property during landfalling tropical cyclones (TCs). The threat is often spread over a wide area within the TC envelope and must be continually evaluated as the TC moves inland and dissipates. To anticipate the risk of TCTORs, forecasters may use high-resolution, rapidly updating model analyses and short-range forecasts such as the Rapid Refresh (RAP) and High-Resolution Rapid Refresh (HRRR), and an ingredients-based approach similar to that used for forecasting continental midlatitude tornadoes. Though RAP and HRRR errors have been identified in typical midlatitude convective environments, this study evaluates the performance of the RAP and the HRRR within the TC envelope, with particular attention given to sounding-derived parameters previously identified as useful for TCTOR forecasting. A sample of 1730 observed upper-air soundings is sourced from 13 TCs that made landfall along the U.S. coastline between 2017 and 2019. The observed soundings are paired with their corresponding model gridpoint soundings from the RAP analysis, RAP 12-h forecast, and HRRR 12-h forecast. Model errors are calculated for both the raw sounding variables of temperature, dewpoint, and wind speed, as well as for the selected sounding-derived parameters. Results show a moist bias that worsens with height across all model runs. There are also statistically significant underpredictions in stability-related parameters such as convective available potential energy (CAPE) and kinematic parameters such as vertical wind shear.

Effects of latent heat in various cloud microphysics processes on autumn rainstorms with different intensities on Hainan Island, China

Abstract. The ability to detect convective regions and to add latent heating to drive convection is one of the most important additions to short-term forecast models such as National Oceanic and Atmospheric Administration's (NOAA's) High-Resolution Rapid Refresh (HRRR) model. Since radars are most directly related to precipitation and are available in high temporal resolution, their data are often used for both detecting convection and estimating latent heating. However, radar data are limited to land areas, largely in developed nations, and early convection is not detectable from radars until drops become large enough to produce significant echoes. Visible and infrared sensors on a geostationary satellite can provide data that are more sensitive to small droplets, but they also have shortcomings: their information is almost exclusively from the cloud top. Relatively new geostationary satellites, Geostationary Operational Environmental Satellite-16 and Satellite-17 (GOES-16 and GOES-17), along with Himawari-8, can make up for this lack of vertical information through the use of very high spatial and temporal resolutions, allowing better observation of bubbling features on convective cloud tops. This study develops two algorithms to detect convection at vertically growing clouds and mature convective clouds using 1 min GOES-16 Advanced Baseline Imager (ABI) data. Two case studies are used to explain the two methods, followed by results applied to 1 month of data over the contiguous United States. Vertically growing clouds in early stages are detected using decreases in brightness temperatures over 10 min. For mature convective clouds which no longer show much of a decrease in brightness temperature, the lumpy texture from rapid development can be observed using 1 min high spatial resolution reflectance data. The detection skills of the two methods are validated against Multi-Radar/Multi-Sensor System (MRMS), a ground-based radar product. With the contingency table, results applying both methods to 1-month data show a relatively low false alarm rate of 14.4 % but missed 54.7 % of convective clouds detected by the radar product. These convective clouds were missed largely due to less lumpy texture, which is mostly caused by optically thick cloud shields above.

Midlatitude storms transport warm and moist air poleward and upward, releasing latent heat. Latent heating is thus organized by thecirculation but then modifies temperature gradients and winds, constituting a nonlinear feedback. We define the latent heating feedbackas the effects that arise from latent heating being coupled with the circulation. Because of its nonlinearity, the climatic effects of thisfeedback are difficult to isolate and remain poorly understood.By decoupling latent heating from the circulation in an atmospheric general circulation model, we show that the latent heating feedbackenhances storm track eddy diffusivity, modifying eddy heat fluxes beyond changes in mean baroclinicity. Simultaneously, tracked stormsoccur at lower latitudes, intensify more, and propagate further poleward, while the subtropical jet strengthens as coupled latent heatingpreserves lower latitude baroclinicity. The feedback response supports the idea that diabatic effects cause the “too zonal, tooequatorward” storm track biases in climate models.Finally, we extend the analysis to climate change experiments where we isolate the contribution from the latent heating feedback onstorm intensity and eddy kinetic energy as the world warms. The feedback is most important in summer where it accounts for most of thechanges in eddy kinetic energy. In winter, the feedback is constrained. Isolating the latent heatingfeedback helps to quantify how storminess changes as the atmosphere warms, which climate models currently struggle with.

Relatively little is known about how topography affects convective storms. The first step toward understanding these effects is to investigate how topography affects storm environments. Unfortunately, the effects of topography on convective environments are not easily observed directly. Instead, it is necessary to resort to using output from the High-Resolution Rapid Refresh (HRRR). The HRRR’s 3-km grid spacing can resolve some of the larger-scale topographic effects. Popular convective storm forecasting parameters obtained from the HRRR are averaged on convective days from February to September 2013–15. It is surmised that most of the day-to-day variability attributable to synoptic- and mesoscale meteorological influences is removed by averaging; the remaining horizontal heterogeneity in parameters related to instability and vertical wind shear is due to the hemispheric-scale meridional temperature and pressure gradient, and likely also topographic influences, especially where recurring longitudinal variations in instability, wind shear, etc. are found. Anomalies are sensitive to the ambient low-level wind direction (i.e., whether winds are locally blowing upslope or downslope), especially for parameters that depend on the low-level vertical shear. The statistical significance of local maxima and minima is demonstrated by comparing the amplitudes of the anomalies to bootstrapped estimates of the standard errors.

Rainfall is a fundamental process within the Earth's hydrological cycle because it represents a principal forcing term in surface water budgets, while its energetics corollary, latent heating, is the principal source of atmospheric diabatic heating well into the middle latitudes. Latent heat production itself is a consequence of phase changes between the vapor, liquid, and frozen states of water. The properties of the vertical distribution of latent heat release modulate large-scale meridional and zonal circulations within the Tropics, as well as modify the energetic efficiencies of midlatitude weather systems. This paper highlights the retrieval of latent heating from satellite measurements generated by the Tropical Rainfall Measuring Mission (TRMM) satellite observatory, which was launched in November 1997 as a joint American–Japanese space endeavor. Since then, TRMM measurements have been providing credible four-dimensional accounts of rainfall over the global Tropics and subtropics, information that can be used to estimate the space–time structure of latent heating across the Earth's low latitudes. A set of algorithm methodologies for estimating latent heating based on precipitation-rate profile retrievals obtained from TRMM measurements has been under continuous development since the advent of the mission s research program. These algorithms are briefly described, followed by a discussion of the latent heating products that they generate. The paper then provides an overview of how TRMM-derived latent heating information is currently being used in conjunction with global weather and climate models, concluding with remarks intended to stimulate further research on latent heating retrieval from satellites.

The active roles of sensible heat supply from the Kuroshio/Kuroshio Extension in the rapid development of an extratropical cyclone, which occurred in the middle of January 2013, were examined by using a regional cloud-resolving model. In this study, a control experiment and three sensitivity experiments without sensible and latent heat fluxes from the warm currents were conducted. When the cyclone intensified, sensible heat fluxes from these currents become prominent around the cold conveyor belt (CCB) in the control run. Comparisons among the four runs revealed that the sensible heat supply facilitates deepening of the cyclone’s central pressure, CCB development, and enhanced latent heating over the bent-back front. The sensible heat supply enhances convectively unstable conditions within the atmospheric boundary layer along the CCB. The increased convective instability is released by the forced ascent associated with frontogenesis around the bent-back front, eventually promoting updraft and resultant latent heating. Additionally, the sensible heating leads to an increase in the water vapor content of the saturated air related to the CCB through an increase in the saturation mixing ratio. This increased water vapor content reinforces the moisture flux convergence at the bent-back front, contributing to the activation of latent heating. Previous research has proposed a positive feedback process between the CCB and latent heating over the bent-back front in terms of moisture supply from warm currents. Considering the above two effects of the sensible heat supply, this study revises the positive feedback process.

<p>Heavy precipitation forecast remains one of the biggest problems in numerical weather prediction. Modern remote sensing systems allow tracking of rapidly developing convective processes and provide additional data for numerical weather models practically in real time. Assimilation of Doppler weather radar data also allows to specify the position and intensity of convective processes in atmospheric numerical models.</p><p>The primary objective of this study is to evaluate the impact of Doppler  radar reflectivity and velocity assimilation in the WRF-ARW mesoscale model for the territory of Belarus in different seasons of the year. Specifically, we focus on the short-range numerical forecasting of mesoscale convective systems passage over the territory of Belarus in 2017-2019 with assimilated radar data.</p><p>Proceeding with weather radar observations available for our cases, we first perform the necessary processing of the raw radar data to eliminate noise, reflections and other kinds of clutter. For identification of non-meteorological noise fuzzy echo classification was used. Then we use the WRF-DA (3D-Var) system to assimilate the processed radar observations from 3 Belarusian Doppler weather radar in the WRF model. Assimilating both radar reflectivity and radial velocity data in the model we aim to better represent not only the distribution of clouds and their moisture content, but also the detailed dynamical aspects of convective circulation. Finally, we analyze WRF modelling output obtained with assimilated radar data and compare it with available meteorological observations and with other model runs (including control runs with no data assimilation or with assimilation of conventional weather stations data only), paying special attention to the accuracy of precipitation forecast 12 hours in advance.</p>

Monitoring the fuel moisture content (FMC) of 10 h dead vegetation is crucial for managing and mitigating the impact of wildland fires. The combination of in situ FMC observations, numerical weather prediction (NWP) models, and satellite retrievals has facilitated the development of machine learning (ML) models to estimate 10 h dead FMC retrievals over the contiguous US (CONUS). In this study, ML models were trained using variables from the National Water Model, the High-Resolution Rapid Refresh (HRRR) NWP model, and static surface properties, along with surface reflectances and land surface temperature (LST) retrievals from the Visible Infrared Imaging Radiometer Suite (VIIRS) instrument on the Suomi-NPP satellite system. Extensive hyper-parameter optimization resulted in skillful FMC models compared to a daily climatography RMSE (+44%) and an hourly climatography RMSE (+24%). Notably, VIIRS retrievals played a significant role as predictors for estimating 10 h dead FMC, demonstrating their importance as a group due to their high band correlation. Conversely, individual predictors within the HRRR group exhibited relatively high importance according to explainability techniques. Removing both HRRR and VIIRS retrievals as model inputs led to a significant decline in performance, particularly with worse RMSE values when excluding VIIRS retrievals. The importance of the VIIRS predictor group reinforces the dynamic relationship between 10 h dead fuel, the atmosphere, and soil moisture. These findings underscore the significance of selecting appropriate data sources when utilizing ML models for FMC prediction. VIIRS retrievals, in combination with selected HRRR variables, emerge as critical components in achieving skillful FMC estimates.