Bragg Scatter Detection by the WSR-88D. Part I: Algorithm Development
Total Page:16
File Type:pdf, Size:1020Kb
VOLUME 34 JOURNAL OF ATMOSPHERIC AND OCEANIC TECHNOLOGY MARCH 2017 Bragg Scatter Detection by the WSR-88D. Part I: Algorithm Development a,b c a LINDSEY M. RICHARDSON, JEFFREY G. CUNNINGHAM, W. DAVID ZITTEL, a a,c d,e ROBERT R. LEE, RICHARD L. ICE, VALERY M. MELNIKOV, f,g f,h NICOLE P. HOBAN, AND JOSHUA G. GEBAUER a Radar Operations Center, National Weather Service, Norman, Oklahoma b Centuria Corporation, Norman, Oklahoma c 557th Weather Wing, Offutt Air Force Base, Nebraska d Cooperative Institute for Mesoscale Meteorological Studies, University of Oklahoma, Norman, Oklahoma e National Severe Storms Laboratory, Norman, Oklahoma f Research Experiences for Undergraduates, National Weather Center, Norman, Oklahoma g North Carolina State University, Raleigh, North Carolina h University of Oklahoma, Norman, Oklahoma (Manuscript received 25 January 2016, in final form 6 December 2016) ABSTRACT Studies have shown that echo returns from clear-air Bragg scatter (CABS) can be used to detect the height of the convective boundary layer and to assess the systematic differential reflectivity (ZDR) bias for a radar site. However, these studies did not use data from operational Weather Surveillance Radar-1988 Doppler (WSR-88D) or data from a large variety of sites. A new algorithm to automatically detect CABS from any operational WSR-88D with dual-polarization capability while excluding contamination from precipitation, biota, and ground clutter is presented here. Visual confirmation and tests related to the sounding parameters’ relative humidity slope, refractivity gradient, and gradient Richardson number are used to assess the algo- rithm. Results show that automated detection of CABS in operational WSR-88D data gives useful ZDR bias information while omitting the majority of contaminated cases. Such an algorithm holds potential for radar calibration efforts and Bragg scatter studies in general. 1. Introduction Bragg scatter is caused by turbulent inhomogeneities with sizes around one-half of a transmitted radar wave- a. Bragg scatter background length (e.g., Cowley 1995; Hardy and Katz 1969; Knight Since the completion of the dual-polarization upgrade and Miller 1993; Doviak and Zrnic´ 2006, chapter 11). to the Weather Surveillance Radar-1988 Doppler (WSR- Bragg scatter has been observed as a layer in clear air and 88D), users have been able to make more detailed developing clouds, and it is mostly associated with re- assessments of the atmosphere. One phenomenon of in- fractivity gradients (Atlas 1959; Ottersten 1969). Studies terest is clear-air Bragg scatter (CABS). Specifically, it has have found that refractivity gradients are related to gra- been shown that CABS can be used to assist with radar dients of moisture in dynamically unstable regions calibration and detection of the height of the convective (Ottersten 1969; Hardy and Katz 1969; Hardy and boundary layer (Melnikov et al. 2011, 2013a; Cunningham Ottersten 1969). In maritime environments, moisture et al. 2013; Zittel et al. 2014). This paper describes the proved to be a more important factor than temperature, development and testing of an automated algorithm for and the primary generation mechanisms consisted of detecting CABS on operational WSR-88Ds. turbulent mixing and detrainment/entrainment of cloudy air (Knight and Miller 1993, 1998; Cohn 1994; Gage et al. 1999; Davison et al. 2013a,b). Corresponding author e-mail: Lindsey M. Richardson, lindsey.m. Melnikov et al. (2011, 2013a) used Bragg scatter to [email protected] detect the height of the turbulent convective boundary DOI: 10.1175/JTECH-D-16-0030.1 For information regarding reuse of this content and general copyright information, consult the AMS Copyright Policy (www.ametsoc.org/ PUBSReuseLicenses). 465 466 JOURNAL OF ATMOSPHERIC AND OCEANIC TECHNOLOGY VOLUME 34 layer in a continental environment. Information about (e.g., Melnikov et al. 2011)wereusedtoformulatea the dual-polarization field characteristics allowed them specific setup for capturing CABS data. Reflectivity (Z), to easily distinguish biota and precipitation from CABS. velocity (V), spectrum width (W), correlation coefficient Their results showed that Bragg scatter corresponds (rHV or CC), and signal-to-noise ratio (SNR) are used to with strong vertical gradients of humidity. Both Davison identify a potential CABS layer and filter out contami- et al. and Melnikov et al. suggest that soundings could be nants (described further in sections 2d and 2e). These used to indicate layers conducive to producing Bragg fields are collectively referred to as base data.3 scattering (discussed in section 3b). There are six main steps for identifying a potential CABS layer and estimating the Z bias from the re- b. Practical application of Bragg scatter DR turns on a WSR-88D: The turbulent eddies associated with CABS are ran- 1) Find radar data within certain spatial limits over a domly oriented and thus have an intrinsic differential certain time domain. reflectivity (Z ) near 0 dB on an unbiased system DR 2) Create a Z histogram of data within the spatial and (Melnikov et al. 2011, 2013a). External targets with known temporal limits. intrinsic Z values,suchasCABS,canbeusedtoin- DR 3) Create a separate Z histogram of range gates that dependently verify the systematic bias of Z (Z bias) DR DR DR pass base data filters using the same spatial limits of a radar site. Thus, CABS without contamination can applied over the specified time domain. be a potential estimator of the Z bias. The Z bias DR DR 4) Use statistical filters to assess Z histogram data for within the WSR-88D can be introduced via an engineering- DR statistical validity and potential contamination from derived internal parameter known as ZDR .This Offset non-Bragg sources. ZDR is applied automatically to the measured Z Offset DR 5) Use a precipitation filter to further reduce the likeli- field (Cunningham et al. 2013; Melnikov et al. 2013b). If the hood of contamination. ZDR fails to capture correctly some aspect of hard- Offset 6) Calculate the mode of the Z histogram if it passes ware bias, the result is a bias in Z . Thus, an error in DR DR all previous filters. ZDROffset translates to a ZDR bias. CABS returns with ZDR estimates not near 0 dB reflect a bias in the radar. For this study, specific range limits, elevation angles, and Results from initial development and testing of an volume coverage patterns (VCPs) mitigate clutter and automated algorithm to collect ZDR bias estimates from some precipitation contamination (described further in CABS in radar data are presented here. Datasets in this section 2b). Values of Z from passing range gates within study span from October 2013 through September 2014 the spatial limits are collected into a histogram over a from over 130 WSR-88D radar sites across the United certain time window for further statistical testing. We used States.1 Archived operational Level II2 data were the 1700–1900 UTC time window for our initial testing, processed offline in a MATLAB environment from though CABS is not limited to this time domain (section sites across the contiguous United States (CONUS) 2c). A ZDR histogram is created using the same VCP, and outside-CONUS (OCONUS) sites, such as Alaska, spatial, and temporal limits, but only range gates that pass Hawaii, and Puerto Rico. Section 2 describes an algorithm thebasedatafilters(section 2d) are included in the dis- to detect CABS with the operational dual-polarization tribution. A test for the sample size and spread of the ZDR WSR-88Ds with a focus on ZDR calibration aspects. histogram is used to check the statistical validity of the Section 3 covers results from visual confirmation and distribution and potential contamination (section 2e). compares the algorithm output to sounding tests fol- Next, the separate Z histogram for precipitation contam- lowed by a summary and a discussion in section 4. ination (section 2f) is checked. Finally, once these filters arepassed,themodefromtheZDR histogram is calculated as an assessment of Z bias from CABS-like returns. 2. Bragg scatter detection algorithm DR b. Volume coverage pattern and range limits a. Algorithm overview CABS layers are expected to lie at the top of a Preliminary visual confirmation and knowledge of boundary layer (convective or marine) a majority of the the dual-polarization characteristics of Bragg scatter 1 3 Radar data were not always available every day for each radar The differential phase (FDP) is also a dual-polarization base site due to radar downtime/data feed errors. variable, but the FDP characteristics of CABS is similar to the 2 See Crum et al. (1993) for the distinction between the various characteristics of light precipitation. Thus, we opted to not use FDP levels of radar data available. for this study. MARCH 2017 R I C H A R D S O N E T A L . 467 time. Several studies show that the convective at 2.58. Additionally, many sites use VCP 31 exclu- boundary layer is generally no higher than 3 km sively in light precipitation (especially snow) cases. (;10 000 ft) above ground level (AGL; Kaimal et al. c. Time considerations 1982; vanZanten et al. 1999; Stensrud 2007, chapter 5; Heinselman et al. 2009). To cover the heights where For initial widespread testing, a time window of 1700– CABS has been observed, a range limit of 10 # R # 1900 UTC was chosen to analyze the October 2013– 80 km was selected. Given typical WSR-88D scanning September 2014 data. This time frame corresponds well angles, CABS is less likely to be detected beyond with heating and convective boundary layer mixing (in 80 km (;43 n mi) and still fill the radar beam4 in op- the central plains of the United States). However, it erational WSR-88D volume scans. Additionally, data could be too early for some western areas and too late within 10 km (;5 n mi) are susceptible to ground for some eastern areas.