1. Introduction
The intensity of Bragg backscatter from refractive index perturbations, at scales half the centimetric and metric wavelengths of atmospheric radars, is measured to estimate the refractive index structure parameter
Monitoring of the CBL is very important for forecasting the timing and likelihood of storm initiation. Heinselman et al. (2009) and Elmore et al. (2012) show that if the reflectivity field obtained with the Weather Surveillance Radar-1988 Doppler (WSR-88D) in “clear air” exhibits an elevated maximum, its height correlates well with the top of the CBL. Our observations show that this type of reflectivity field is only one of several types that can be detected from Bragg scatterers, which we describe herein. To observe the “fine” structure of
In absence of echoes from atmospheric biota, radar wind profilers estimate
In the next section, we discuss types of scatter from and above CBL observed with KOUN running the enhanced detection capability mode.
2. Types of structures from Bragg scatter
a. Two layers of Bragg scatter
Various types of reflectivity fields in clear air are apparent from observations with the KOUN dual-polarization radar. An example of elevated layers is from 1627 UTC 19 February 2008 (Fig. 1a), where three layers of echo are seen.
The lowest layer is from the ground to around 500 m above ground level (AGL) and stretches horizontally across the entire cross section from the radar location to 50 km. The 1200 UTC Norman, Oklahoma, sounding captured a strong stable layer from the ground to 950 hPa with weak northerly winds (Fig. 2a). This stable layer is associated with a shallow frontal boundary that moved into southern Oklahoma overnight. Above the shallow stable layer is a nearly isothermal layer that extends to 900 hPa. By the time of the 1627 UTC radar observations, the 2-m temperatures approached 9°C, suggesting a convective boundary layer depth of approximately 500 m from the sounding. This depth is estimated by assuming that a dry adiabatic layer, characteristic of CBLs, stretches from the surface upward until it connects with the temperature profile from the 1200 UTC sounding. The potential temperature of this layer is calculated using the observed 2-m temperature and surface pressure, and it is assumed that the sounding structure above the CBL does not evolve with time. An examination of Rapid Update Cycle (RUC) analyses (Benjamin et al. 2004) from 1700 UTC also indicates a boundary layer depth of near 500 m. These boundary layer depth estimates agree well with the depth of the lowest echo layer from the KOUN radar (Fig. 1a). This result suggests that the dual-polarization radar observations can provide information on the depth of the convective boundary layer.
Two other layered echoes are seen between 1.5 and 2.5 km and between 7.5 and 9.5 km AGL (Fig. 1a). Such layers appear differently in different azimuthal directions. Visible and infrared satellite imagery (not shown) indicates the presence of thin upper-level clouds across central Oklahoma. Observers from Tinker Air Force Base (AFB) and Oklahoma City at 1650 UTC reported cloud bases at 7.6 and 8.5 km AGL, respectively. These reported cloud bases agree well with the layer between 7.5 and 9.5 km in Fig. 1a, indicating the layer is produced by particulate scattering (ZDR of about 1 dB) from nonprecipitating clouds. The potential of the WSR-88D to map cloud structures is shown by Melnikov et al. (2011b). However, no clouds are reported at lower levels, whereas the radar shows a layer near 2 km AGL. The radar operator [virtual machine monitor (VMM)] visually confirmed the lack of low-level clouds at this time during the data collection period. Furthermore, the corresponding 1700 UTC RUC analysis has relative humidity (RH) values below 50% between 2 and 3 km AGL, suggesting that sufficient moisture is not present for cloud formation. Closest in time rawinsonde profiles from 1200 UTC (Fig. 2a) exhibit a strong gradient of relative humidity (thus, a strong gradient in potential refractive index), wind speed, and direction near the height of the layer. The top of this echoing layer is just below the stable temperature stratification seen in Fig. 2a. Thus, we hypothesize that this layered radar echo is produced by intense refractive index perturbations associated with shear-induced turbulent mixing within a strong vertical gradient of potential refractive index. The height of maximal reflectivity obtained from a nearby National Oceanic and Atmospheric Administration (NOAA) wind profiler (Fig. 1b) also is in accord with the height of maximum reflectivity observed with KOUN, providing further support for the hypothesis that shear-induced turbulent mixing produced the layer of larger
b. Convective mixing above a layer of Bragg scatter
Another example of the potential use of
More curious is the layer of enhanced scatter from 1.5 to 4 km AGL (Fig. 3). Although conditions at the radar site indicated clear skies directly overhead, broken cloud cover was reported at Oklahoma City and Tinker AFB, with cloud bases at 1.5 and 2.0 km, respectively, although cloudy conditions were observed at KOUN an hour earlier. The horizontal distance between KOUN and Oklahoma City is roughly 32 km, and the radar cross section in Fig. 3a extends from KOUN toward Oklahoma City. The 0000 UTC sounding from Norman indicates a cloud mixing layer from the top of the boundary layer to around 3 km AGL (Fig. 2b). A visible satellite image (not presented here) shows cumulus congestus clouds in north-central Oklahoma just to the north of KOUN. We submit that developing clouds are responsible for the power increase of echoes observed with the Purcell profiler for the heights between 1.5 and 4 km. The presence of developing clouds suggests that the larger values of
c. A layer of Bragg scatter with more complex inner structure
The final example of Bragg scatter is from the period 1–2 March 2008 (Fig. 4, profiler data are not available for these days). One can see a weakly reflecting layer up to 2 km AGL that stretches horizontally across the entire cross section, with locally stronger scatter at heights of about 1.7 km (Fig. 4a). This stronger reflectivity layer coincides with the height of strong gradients of relative humidity (Fig. 4c) capped with a temperature inversion, indicating the top of the convective boundary layer. Thus, the boundary layer depth estimated from the sounding again agrees well with the depth of the 1.7-km-deep reflectivity layer. However, sunset on 1 March 2008 occurred at 2326 UTC, in between the times of the two cross sections shown in Fig. 4. While the intermittent reflections at heights of about 1.7 km decrease gradually between the two observation times, separated by 53 min, the reflectivity layer depth remains nearly constant. This suggests that the radar is observing a residual layer of active turbulence above the developing stable surface layer as modeled in a large-eddy simulation by Sorbjan (1997). This layer was visible by the radar for a few hours after sunset, suggesting that the radar can provide information on the time scale over which turbulence dissipates in the convective afternoon to the stable evening boundary layer transition.
3. Conclusions
Radar observations from KOUN, running the enhanced detection capability mode, exhibit various types of Bragg scatter: reflectivity layers corresponding to the convective boundary layer (the lowest layers in Figs. 1a, 3a,b, and the layer in Figs. 4a,b), turbulent layers within strong gradients in relative humidity (Fig. 1a), nonprecipitating cloud (Fig. 1a), strong convective plumes (Fig. 3) with developing weakly reflecting clouds above them, and a layer of weak reflectivity associated with decaying boundary layer turbulence (Fig. 4). Refractive index perturbations, associated with the mixing of humidity gradients at the top of the CBL, scatter a sufficient amount of radar signals to be detected by the WSR-88Ds. Thus, one of meteorological applications of mapping of Bragg scatter is monitoring of the depth of the CBL. It was demonstrated in section 2 that enhanced radar detection capability allows observations of not only layers of Bragg scatter but also developing cumulus congestus clouds prior to the onset of thunderstorms.
Our analysis of Bragg scatter structures shows that the tops of the lowest echo layers correlate well with estimates of the CBL depth obtained from atmospheric rawinsonde soundings. Radar data also indicate local inhomogeneities in the CBL. Radar observations of CBL depth could provide an important constraint on the changes in water vapor, pollutants, and turbulence within the boundary layer. Current model predictions of CBL depth often differ from observations by a factor of 2 (Bright and Mullen 2002; Stensrud and Weiss 2002), suggesting that routine observations of CBL depth would provide new information that could be used advantageously in data assimilation systems.
We collected radar data and presented our results in the RHI format because this format shows the vertical profiles naturally, whereas the operational WSR-88Ds collect data in conical scans [plan position indicator (PPI) format)]. However, RHIs can easily be obtained from a dense collection of PPIs. Because the monitoring of the boundary layer is typically performed in a prestorm environment, there is sufficient time to make a dense PPI scanning. We have collected data in the short-pulse mode (the range resolution of 250 m) to have a more meaningful estimate of
Acknowledgments
The authors thank three anonymous reviewers for their constructive and helpful comments, which improved the presentation. Funding for this study was provided by NOAA's Office of Oceanic and Atmospheric Research under NOAA–University of Oklahoma Cooperative Agreement NA17RJ1227 U.S. Department of Commerce.
REFERENCES
Benjamin, S. G., and Coauthors, 2004: An hourly assimilation–forecast cycle: The RUC. Mon. Wea. Rev., 132, 495–518.
Bright, D. R., and Mullen S. L. , 2002: Short-range ensemble forecasts of precipitation during the Southwest monsoon. Wea. Forecasting, 17, 1080–1100.
Doviak, R. J., and Berger M. J. , 1980: Turbulence and waves in the optically clear planetary boundary layer. Radio Sci., 15, 297–317.
Doviak, R. J., and Zrnić D. S. , 2006: Doppler Radar and Weather Observations. 2nd ed. Dover Publications, 562 pp.
Elmore, K. L., Heinselman P. L. , and Stensrud D. J. , 2012: Using WSR-88D data and insolation estimates to determine convective boundary layer depth. J. Atmos. Oceanic Technol., 29, 581–588.
Fairall, C. W., 1991: The humidity and temperature sensitivity of clear-air radars in the convective boundary layer. J. Appl. Meteor., 30, 1064–1074.
Heinselman, P. L., Stensrud D. J. , Hluchan R. M. , Spencer P. L. , Burke P. C. , and Elmore K. L. , 2009: Radar reflectivity–based estimates of mixed layer depth. J. Atmos. Oceanic Technol., 26, 229–239.
Holleman, I., van Gasteren H. , and Bouten W. , 2008: Quality assessment of weather radar wind profiles during bird migration. J. Atmos. Oceanic Technol., 25, 2188–2198.
Lang, T. J., Rutledge S. A. , and Smith J. L. , 2004: Observations of quasi-symmetric echo patterns in clear air with the CSU–CHILL polarimetric radar. J. Atmos. Oceanic Technol., 21, 1182–1189.
Lehmann, V., 2012: Optimal Gabor-frame-expansion-based intermittent-clutter-filtering method for radar wind profiler. J. Atmos. Oceanic Technol., 29, 141–158.
Melnikov, V. M., Doviak R. J. , Zrnić D. S. , and Stensrud D. J. , 2011a: Mapping Bragg scatter with a polarimetric WSR-88D. J. Atmos. Oceanic Technol., 28, 1273–1285.
Melnikov, V. M., Zrnić D. S. , Doviak R. J. , Chilson P. B. , Mechem D. B. , and Kogan Y. L. , 2011b: Prospects of the WSR-88D radar for cloud studies. J. Appl. Meteor. Climatol., 50, 859–872.
Sorbjan, Z., 1997: Decay of convective turbulence revisited. Bound.-Layer Meteor., 82, 501–515.
Stensrud, D. J., and Weiss S. J. , 2002: Mesoscale model ensemble forecast of the 3 May 1999 tornado outbreak. Wea. Forecasting, 17, 526–543.
Tatarskii, V. I., 1971: The Effects of Turbulent Atmosphere on Wave Propagation. National Technical Information Service, 472 pp.
Wilczak, J. M., and Coauthors, 1995: Contamination of wind profiler data by migrating birds: Characteristics of corrupted data and potential solutions. J. Atmos. Oceanic Technol., 12, 449–467.
Wilson, J. W., Weckwerth T. M. , Vivekanandan J. , Wakimoto R. M. , and Russell R. W. , 1994: Boundary layer clear-air radar echoes: Origin of echoes and accuracy of deriving winds. J. Atmos. Oceanic Technol., 11, 1184–1206.
Wyngaard, J. C., and LeMone M. A. , 1980: Behavior of the refractive index structure parameter in the entraining convective boundary layer. J. Atmos. Sci., 37, 1573–1585.
Zrnić, D. S., and Ryzhkov A. V. , 1998: Observations of insects and birds with a polarimetric radar. IEEE Trans. Geosci. Remote Sens., 36, 661–668.