skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


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
Award ID(s):
1905736
PAR ID:
10512042
Author(s) / Creator(s):
; ;
Publisher / Repository:
EGU
Date Published:
Journal Name:
Atmospheric Measurement Techniques
Volume:
17
Issue:
11
ISSN:
1867-8548
Page Range / eLocation ID:
3377 to 3399
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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. 
    more » « less
  2. 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. 
    more » « less
  3. 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) 
    more » « less
  4. 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) 
    more » « less
  5. Abstract The Echo Classification from COnvectivity (ECCO) algorithm identifies convective and stratiform types of radar echo in three dimensions. It is based on the calculation of reflectivity texture—a combination of the intensity and the heterogeneity of the radar echoes on each horizontal plane in a 3D Cartesian volume. Reflectivity texture is translated into convectivity, which is designed to be a quantitative measure of the convective nature of each 3D radar grid point. It ranges from 0 (100% stratiform) to 1 (100% convective). By thresholding convectivity, a more traditional qualitative categorization is obtained, which classifies radar echoes as convective, mixed, or stratiform. In contrast to previous algorithms, these echo-type classifications are provided on the full 3D grid of the reflectivity field. The vertically resolved classifications, in combination with temperature data, allow for subclassifications into shallow, mid-, deep, and elevated convective features, and low, mid-, and high stratiform regions—again in three dimensions. The algorithm was validated using datasets collected over the U.S. Great Plains during the PECAN field campaign. An analysis of lightning counts shows ∼90% of lightning occurring in regions classified as convective by ECCO. A statistical comparison of ECCO echo types with the well-established GPM radar precipitation-type categories show 84% (88%) of GPM stratiform (convective) echo being classified as stratiform (convective) or mixed by ECCO. ECCO was applied to radar grids for the continental United States, the United Arab Emirates, Australia, and Europe, illustrating its robustness and adaptability to different radar grid characteristics and climatic regions. 
    more » « less