Title: Dual adaptive differential threshold method for automated detection of faint and strong echo features in radar observations of winter storms
Abstract. Radar observations of winter storms often exhibit locally enhanced linear features in reflectivity, sometimes labeled as snow bands. We have developed a new, objective method for detecting locally enhanced echo features in radar data from winter storms. In comparison to convective cells in warm season precipitation, these features are usually less distinct from the background echo and often have more fuzzy or feathered edges. This technique identifies both prominent, strong features and more subtle, faint features. A key difference from previous radar reflectivity feature detection algorithms is the combined use of two adaptive differential thresholds, one that decreases with increasing background values and one that increases with increasing background values. The algorithm detects features within a snow rate field rather than reflectivity and incorporates an underestimate and overestimate of feature areas to account for uncertainties in the detection. We demonstrate the technique on several examples from the US National Weather Service operational radar network. The feature detection algorithm is highly customizable and can be tuned for a variety of data sets and applications. more »« less
Tomkins, Laura M.; Yuter, Sandra E.; Miller, Matthew A.; Allen, Luke R.
(, Atmospheric Measurement Techniques)
Abstract. In winter storms, enhanced radar reflectivity is often associated with heavy snow. However, some higher reflectivities are the result of mixed precipitation including melting snow. The correlation coefficient (a dual-polarization radar variable) can identify regions of mixed precipitation, but this information is usually presented separately from reflectivity. Especially under time pressure, radar data users can mistake regions of mixed precipitation for heavy snow because of the high cognitive load associated with comparing data in two fields while simultaneously attempting to discount a portion of the high reflectivity values. We developed an image muting method for regional radar maps that visually de-emphasizes the high reflectivity values associated with mixed precipitation. These image muted depictions of winter storm precipitation structures are useful for analyzing regions of heavy snow and monitoring real-time weather conditions.
Tobin, Dana M.; Kumjian, Matthew R.; Oue, Mariko; Kollias, Pavlos
(, Journal of the Atmospheric Sciences)
Abstract The discovery of a polarimetric radar signature indicative of hydrometeor refreezing has shown promise in its utility to identify periods of ice pellet production. Uniquely characterized well below the melting layer by locally enhanced values of differential reflectivity ( Z DR ) within a layer of decreasing radar reflectivity factor at horizontal polarization ( Z H ), the signature has been documented in cases where hydrometeors were completely melted prior to refreezing. However, polarimetric radar features associated with the refreezing of partially melted hydrometeors have not been examined as rigorously in either an observational or microphysical modeling framework. Here, polarimetric radar data—including vertically pointing Doppler spectral data from the Ka-band Scanning Polarimetric Radar (KASPR)—are analyzed for an ice pellets and rain mixture event where the ice pellets formed via the refreezing of partially melted hydrometeors. Observations show that no such distinct localized Z DR enhancement is present, and that values instead decrease directly beneath enhanced values associated with melting. A simplified, explicit bin microphysical model is then developed to simulate the refreezing of partially melted hydrometeors, and coupled to a polarimetric radar forward operator to examine the impacts of such refreezing on simulated radar variables. Simulated vertical profiles of polarimetric radar variables and Doppler spectra have similar features to observations, and confirm that a Z DR enhancement is not produced. This suggests the possibility of two distinct polarimetric features of hydrometeor refreezing: ones associated with refreezing of completely melted hydrometeors, and those associated with refreezing of partially melted hydrometeors. Significance Statement There exist two pathways for the formation of ice pellets: refreezing of fully melted hydrometeors, and refreezing of partially melted hydrometeors. A polarimetric radar signature indicative of fully melted hydrometeor refreezing has been extensively documented in the past, yet no study has documented the refreezing of partially melted hydrometeors. Here, observations and idealized modeling simulations are presented to show different polarimetric radar features associated with partially melted hydrometeor refreezing. The distinction in polarimetric features may be beneficial to identifying layers of supercooled liquid drops within transitional winter storms.
This dataset contains over 14,000 hours of regional radar mosaics over the northeast US from 600+ winter storm days between 1996-2023. Winter storm days are defined when at least 2 out of 15 surface stations in the northeast US (see attached map) produced at least 1 inch of snow over the 24 hour period. Sequences of these mosaics aid in analyzing the precipitation area and the structures within winter storms. Radar reflectivity data is combined from the first, lowest (0.5 degree) elevation angle from 12 NEXRAD WSR-88D radars in the northeast US (see attached). The scans occur every 5-10 minutes from each radar depending on the radar scan settings. The time label of the regional map is based on the scan time central radar, KOKX (Upton, NY). Scans from other radars in the region are used for that time as long as they are within 8 minutes of the KOKX scan. The polar radar data from each radar is interpolated to a regional 1202 km x 1202 km Cartesian grid with 2 km grid spacing covering 35.73-46.8 degN and 66.36-81.85 degW. Where the radar domains overlap, we take the highest reflectivity value. For dates after dual-polarization integration (2012 onwards), files contain the correlation coefficient (RHO_HV) field and a binary field that can be used to “image mute” the reflectivity which reduces the visual prominence of melting and mixed precipitation commonly mistaken for heavy snow. Image muting is applied where radar reflectivity is ≥ 20 dBZ and RHO_HV is ≤ 0.97. This product is different from other widely used radar mosaics such as the MRMS produced by NOAA since it does not interpolate to a constant altitude and thus preserves the finer scale details in the reflectivity field. Because the data used to create these mosaics are not interpolated to a constant altitude, the altitude varies over the region (altitudes of radar scan used at each grid point are provided as a field for each data file). This data set is specifically designed to analyze fine-scale structures in winter storms. Part 1 contains files pre-dual polarization integration (1996-2012)Part 2 contains files post-dual polarization integration (2012-2023)
This dataset contains over 14,000 hours of regional radar mosaics over the northeast US from 600+ winter storm days between 1996-2023. Winter storm days are defined when at least 2 out of 15 surface stations in the northeast US (see attached map) produced at least 1 inch of snow over the 24 hour period. Sequences of these mosaics aid in analyzing the precipitation area and the structures within winter storms. Radar reflectivity data is combined from the first, lowest (0.5 degree) elevation angle from 12 NEXRAD WSR-88D radars in the northeast US (see attached). The scans occur every 5-10 minutes from each radar depending on the radar scan settings. The time label of the regional map is based on the scan time central radar, KOKX (Upton, NY). Scans from other radars in the region are used for that time as long as they are within 8 minutes of the KOKX scan. The polar radar data from each radar is interpolated to a regional 1202 km x 1202 km Cartesian grid with 2 km grid spacing covering 35.73-46.8 degN and 66.36-81.85 degW. Where the radar domains overlap, we take the highest reflectivity value. For dates after dual-polarization integration (2012 onwards), files contain the correlation coefficient (RHO_HV) field and a binary field that can be used to “image mute” the reflectivity which reduces the visual prominence of melting and mixed precipitation commonly mistaken for heavy snow. Image muting is applied where radar reflectivity is ≥ 20 dBZ and RHO_HV is ≤ 0.97. This product is different from other widely used radar mosaics such as the MRMS produced by NOAA since it does not interpolate to a constant altitude and thus preserves the finer scale details in the reflectivity field. Because the data used to create these mosaics are not interpolated to a constant altitude, the altitude varies over the region (altitudes of radar scan used at each grid point are provided as a field for each data file). This data set is specifically designed to analyze fine-scale structures in winter storms. Part 1 contains files pre-dual polarization integration (1996-2012)Part 2 contains files post-dual polarization integration (2012-2023)
Majewski, Adam; French, Jeffrey R
(, Journal of the Atmospheric Sciences)
Abstract Shear and buoyancy gradients, often observed in midlatitude baroclinic and orographic winter storms, produce discrete layers of turbulence. These turbulent layers modify the distribution of supercooled liquid water (SLW), whose presence enables hydrometeors to grow to precipitation sizes faster than through vapor-ice deposition alone. Both the Wegener–Bergeron–Findeisen process and riming require SLW—heterogeneously distributed in response to the dynamic forcings at all superposed scales. The University of Wyoming W-band cloud radar Doppler spectrum width characterizes the air motion turbulence intensity in mixed-phase layer clouds after avoiding fall speed dominated regions through comparison to coarser radar turbulence metrics. Embedded layers of turbulent air motion are compared to quiescent cloud regions in either/both the upwind and downwind directions. Median radar reflectivity profiles characterize the vertical growth of hydrometeors in the vicinity of identified layers, and differences in these vertical reflectivity gradients comparing turbulent to nonturbulent regions quantify enhanced hydrometeor growth over the layer. Over the entirety of the Seeded and Natural Orographic Wintertime Clouds—the Idaho Experiment, this parameter demonstrates a statistically significant increase, −13.6 dBZekm−1(from −1.7 to −24.5, 95% computed confidence), in radar reflectivity echo power with distance downward for embedded turbulent layers compared to quiescent cloud nearby. The increased vertical particle growth rate for turbulent layers appears to result from spatially heterogeneous phase partitioning, increased SLW mass and extent, and enhanced collision/collection rates in these layers. These first two conditions are examined individually where turbulent layers or fall streaks are sampled in situ, while the latter agrees with modeling results but can only be inferred herein. Significance StatementTurbulent mixing in mixed-phase clouds is understood to enhance cloud hydrometeor growth. This study quantifies these effects on observed airborne W-band radar reflectivity over an entire field campaign targeting midlatitude winter storms and calls into question whether this linkage is diagnosed properly if at all in forecast and bulk microphysical models with coarse (greater than 500 m) grid spacing or vertical resolution.
Tomkins, Laura M, Yuter, Sandra E, and Miller, Matthew A. Dual adaptive differential threshold method for automated detection of faint and strong echo features in radar observations of winter storms. Retrieved from https://par.nsf.gov/biblio/10512042. Atmospheric Measurement Techniques 17.11 Web. doi:10.5194/amt-17-3377-2024.
Tomkins, Laura M, Yuter, Sandra E, & Miller, Matthew A. Dual adaptive differential threshold method for automated detection of faint and strong echo features in radar observations of winter storms. Atmospheric Measurement Techniques, 17 (11). Retrieved from https://par.nsf.gov/biblio/10512042. https://doi.org/10.5194/amt-17-3377-2024
Tomkins, Laura M, Yuter, Sandra E, and Miller, Matthew A.
"Dual adaptive differential threshold method for automated detection of faint and strong echo features in radar observations of winter storms". Atmospheric Measurement Techniques 17 (11). Country unknown/Code not available: EGU. https://doi.org/10.5194/amt-17-3377-2024.https://par.nsf.gov/biblio/10512042.
@article{osti_10512042,
place = {Country unknown/Code not available},
title = {Dual adaptive differential threshold method for automated detection of faint and strong echo features in radar observations of winter storms},
url = {https://par.nsf.gov/biblio/10512042},
DOI = {10.5194/amt-17-3377-2024},
abstractNote = {Abstract. Radar observations of winter storms often exhibit locally enhanced linear features in reflectivity, sometimes labeled as snow bands. We have developed a new, objective method for detecting locally enhanced echo features in radar data from winter storms. In comparison to convective cells in warm season precipitation, these features are usually less distinct from the background echo and often have more fuzzy or feathered edges. This technique identifies both prominent, strong features and more subtle, faint features. A key difference from previous radar reflectivity feature detection algorithms is the combined use of two adaptive differential thresholds, one that decreases with increasing background values and one that increases with increasing background values. The algorithm detects features within a snow rate field rather than reflectivity and incorporates an underestimate and overestimate of feature areas to account for uncertainties in the detection. We demonstrate the technique on several examples from the US National Weather Service operational radar network. The feature detection algorithm is highly customizable and can be tuned for a variety of data sets and applications.},
journal = {Atmospheric Measurement Techniques},
volume = {17},
number = {11},
publisher = {EGU},
author = {Tomkins, Laura M and Yuter, Sandra E and Miller, Matthew A},
}
Warning: Leaving National Science Foundation Website
You are now leaving the National Science Foundation website to go to a non-government website.
Website:
NSF takes no responsibility for and exercises no control over the views expressed or the accuracy of
the information contained on this site. Also be aware that NSF's privacy policy does not apply to this site.