Quasi-vertical Profiles Research Articles

An Automated Approach to Estimating Convective Boundary Layer Depth from Dual-Polarization WSR-88D Radar Observations

Abstract Convective boundary layer (CBL) depth can be estimated from dual-polarization WSR-88D radars using the product differential reflectivity ZDR because the CBL top is collocated with a local ZDR minimum produced by Bragg scatter at the interface of the CBL and the free troposphere. Quasi-vertical profiles (QVPs) of ZDR are produced for each radar volume scan and profiles from successive times are stitched together to form a time–height plot of ZDR from sunrise to sunset. QVPs of ZDR often show a low-ZDR channel that starts near the ground and rises during the morning and early afternoon, identifying the CBL top. Unfortunately, results show that this channel within the QVP can occasionally be misleading. This motivated creation of a new variable DVar, which combines ZDR with its azimuthal variance and is particularly helpful at identifying the CBL top during the morning hours. Two methods are developed to track the CBL top from QVPs of ZDR and DVar. Although each method has strengths and weaknesses, the best results are found when the two methods are combined using inverse variance weighting. The ability to detect CBL depth from routine WSR-88D radar scans rather than from rawinsondes or lidar instruments would vastly improve our understanding of CBL depth variations in the daytime by increasing the temporal and spatial frequencies of the observations. Significance Statement The daytime convective boundary layer (CBL) can increase in depth from a few hundred to a few thousand meters between sunrise and sunset and is strongly connected to temperature changes at Earth’s surface. Unfortunately, current observations of CBL depth primarily consist of measurements from twice daily rawinsonde launches at 97 locations across the United States. As a result, CBL depth observations lack fine spatial and temporal resolution and miss the daily cycle of CBL growth. This study seeks to fill the gaps in CBL depth observations by developing an automated method to estimate CBL depth from dual-polarization WSR-88D radar observations with a temporal resolution as fine as 5–10 min. These observations will greatly enhance our ability to observe and monitor CBL depth in real time.

Polarimetric radar observation of the melting layer during the pre-summer rainy season over South China

ABSTRACT The geometric and polarimetric characteristics of the melting layer (ML) during the pre-summer rainy season over South China (PRSC) are investigated from 389 S-band polarimetric radar volume scans. The existing ML detection algorithm is slightly improved, in which the ML is automatically detected by using vertical gradients of radar reflectivity factor for horizontal polarizations (Z H), and co-polar cross-correlation coefficient (ρ HV) from quasi-vertical profiles. And these profiles are obtained by azimuthal averaging of polarimetric radar variables at antenna elevation 19.5°. The comparison between the detected ML top and the 0°C isotherm heights from ERA5 demonstrates that the detection algorithm has high accuracy with an offset of 158 m. Due to larger snowflakes, the average ML thickness in PRSC is thicker (i.e. 150 m) than that in other climatic regions. The ML is characterized by large Z H, large differential reflectivity (Z DR) and small ρ HV, with average values of 32 dBZ, 0.91 dB and 0.94. Furthermore, the peak values of Z H and Z DR in the ML are obviously larger than that from rain regions and snow regions, while the opposite is true for ρ HV distribution.

Analysis of the Characteristics and Evolution Mechanisms of a Bow-Shaped Squall Line in East China Observed with Dual-Polarization Doppler Radars

To gain a deeper understanding of the formation and evolutionary mechanisms of a bow-shaped squall line (BSL) that occurred in East China on 10 May 2021, observations from S-band dual-polarization radars, a disdrometer and other instruments are used to investigate the characteristics and evolution of the kinematic, microphysical and radar echo structure within the squall line during its formative and mature stages. The results are as follows. The updraft induced by upper-level divergence and vertical thermal instability induced by the cold source at the middle and top of the troposphere provided environmental conditions suitable for the formation and strengthening of a squall line. The characteristics of the vertical vorticity at the leading edge of the squall line provided a good indication of its echo structure and evolutionary trend. The mechanism behind a new echo phenomenon—double high-differential reflectivity (ZDR) bands—observed in plan position indicator scans produced by the dual-polarization radar is investigated from the kinematic and microphysical structural perspectives. The evolutionary characteristics of the microphysical structure of the bulk of the squall line and its trailing stratiform cloud region are analyzed based on the quasi-vertical profiles retrieved from the S-band dual-polarization radar in Quzhou. Moreover, a conceptual model describing this type of BSL with a trailing region of stratiform rain in the warm sector is developed to provide technical support for the monitoring and early warning of BSLs.

New Polarimetric Radar Algorithm for Melting-Layer Detection and Determination of Its Height

Abstract A novel polarimetric radar algorithm for melting-layer (ML) detection and determination of its height has been developed and tested for a large number of cold-season weather events. The algorithm uses radial profiles of the cross-correlation coefficient (ρhv or CC) at the lowest elevation angles (<5°–6°). The effects of beam broadening on the spatial distribution of CC have been taken into account via theoretical simulations of the radial profiles of CC assuming intrinsic vertical profiles of polarimetric radar variables within the ML with varying heights and depths of the ML. The model radial profiles of CC and their key parameters are stored in lookup tables and compared with the measured CC profiles. The matching of the model and measured CC radial profiles allows the algorithm to determine the “true” heights of the top and bottom of the ML, Ht and Hb, at distances up to 150 km from the radar. Integrating the CC information from all available antenna elevations makes it possible to produce accurate maps of Ht and Hb over large areas of radar coverage as opposed to the previous ML detection methods including the existing algorithm implemented on the U.S. network of the WSR-88Ds. The initial version of the algorithm has been implemented in C++ and tested for a multitude of cold-season weather events characterized by a low ML with different degrees of spatial nonuniformity including cases with sharp frontal boundaries and rain–snow transitions. The new ML detection algorithm (MLDA) exhibits robust performance, demonstrating spatial and temporal continuity, and showing general consistency of the ML designations matching those obtained from the regional model and the quasi-vertical profiles (QVP) methodology output.

Calibration of radar differential reflectivity using quasi-vertical profiles

Abstract. Accurate precipitation estimation with weather radars is essential for hydrological and meteorological applications. The differential reflectivity (ZDR) is a crucial weather radar measurement that helps to improve quantitative precipitation estimates using polarimetric weather radars. However, a system bias between the horizontal and vertical channels generated by the radar produces an offset in ZDR. Existing methods to calibrate ZDR measurements rely on the intrinsic values of the ZDR of natural targets (e.g. drizzle or dry snow) collected at high elevation angles (e.g. higher than 40∘ or even at 90∘), in which ZDR values close to 0 dB are expected. However, not all weather radar systems can scan at such high elevation angles or point the antenna vertically to collect precipitation measurements passing overhead. Therefore, there is a need to develop new methods to calibrate ZDR measurements using lower-elevation scans. In this work, we present and analyse a novel method for correcting and monitoring the ZDR offset using quasi-vertical profiles computed from scans collected at 9∘ elevations. The method is applied to radar data collected through 1 year of precipitation events by two operational C-band polarimetric weather radars in the UK. The proposed method shows a relative error of 0.1 dB when evaluated against the traditional approach based on ZDR measurements collected at 90∘ elevations. Additionally, the method is independently assessed using disdrometers located near the radar sites. The results showed a reasonable agreement between disdrometer-derived and radar-calibrated ZDR measurements.

Analysis of the microphysical properties of snowfall using scanning polarimetric and vertically pointing multi-frequency Doppler radars

Abstract. Radar dual-wavelength ratio (DWR) measurements from the Stony Brook Radar Observatory Ka-band scanning polarimetric radar (KASPR, 35 GHz), a W-band profiling radar (94 GHz), and a next-generation K-band (24 GHz) micro rain radar (MRRPro) were exploited for ice particle identification using triple-frequency approaches. The results indicated that two of the radar frequencies (K and Ka band) are not sufficiently separated; thus, the triple-frequency radar approaches had limited success. On the other hand, a joint analysis of DWR, mean Doppler velocity (MDV), and polarimetric radar variables indicated potential in identifying ice particle types and distinguishing among different ice growth processes and even in revealing additional microphysical details. We investigated all DWR pairs in conjunction with MDV from the KASPR profiling measurements and differential reflectivity (ZDR) and specific differential phase (KDP) from the KASPR quasi-vertical profiles. The DWR-versus-MDV diagrams coupled with the polarimetric observables exhibited distinct separations of particle populations attributed to different rime degrees and particle growth processes. In fallstreaks, the 35–94 GHz DWR pair increased with the magnitude of MDV corresponding to the scattering calculations for aggregates with lower degrees of riming. The DWR values further increased at lower altitudes while ZDR slightly decreased, indicating further aggregation. Particle populations with higher rime degrees had a similar increase in DWR but a 1–1.5 m s−1 larger magnitude of MDV and rapid decreases in KDP and ZDR. The analysis also depicted the early stage of riming where ZDR increased with the MDV magnitude collocated with small increases in DWR. This approach will improve quantitative estimations of snow amount and microphysical quantities such as rime mass fraction. The study suggests that triple-frequency measurements are not always necessary for in-depth ice microphysical studies and that dual-frequency polarimetric and Doppler measurements can successfully be used to gain insights into ice hydrometeor microphysics.

Polarimetric Observations and Simulations of Sublimating Snow: Implications for Nowcasting

AbstractSnow sublimating in dry air is a forecasting challenge and can delay the onset of surface snowfall and affect storm-total accumulations. Despite this, it remains comparatively less studied than other microphysical processes. Herein, the characteristics of sublimating snow and the potential for nowcasting snowfall reaching the surface are explored through the use of dual-polarization radar. Twelve cases featuring prolific sublimation were analyzed using range-defined quasi-vertical profiles (RDQVPs) and compared with environmental model analyses. Overall, reflectivity Z significantly decreases, differential reflectivity ZDR slightly decreases, and copolar-correlation coefficient ρhv remains nearly constant through the sublimation layer. Regions of enhanced specific differential phase Kdp were frequently observed in the sublimation layer and are believed to be polarimetric evidence of secondary ice production via sublimation. A 1D bin model was initialized using particle size distributions retrieved from the RDQVPs using numerous novel polarimetric snowretrieval relations for a wide range of forecast lead times, with the model environment evolving in response to sublimation. It was found that the model was largely able to predict the snowfall start time up to six hours in advance, with a 6-h median bias of just -18.5 minutes. A more detailed case study of the 08 December 2013 snowstorm in the Philadelphia region was also performed, demonstrating good correspondence with observations and examples of model fields (e.g., cooling rate) hypothetically available from such a tool. The proof-of-concept results herein demonstrate the potential benefits of incorporating spatially averaged radar data in conjunction with simple 1D models into the nowcasting process.

Detection of the melting level with polarimetric weather radar

Abstract. Accurate estimation of the melting level (ML) is essential in radar rainfall estimation to mitigate the bright band enhancement, classify hydrometeors, correct for rain attenuation and calibrate radar measurements. This paper presents a novel and robust ML-detection algorithm based on either vertical profiles (VPs) or quasi-vertical profiles (QVPs) built from operational polarimetric weather radar scans. The algorithm depends only on data collected by the radar itself, and it is based on the combination of several polarimetric radar measurements to generate an enhanced profile with strong gradients related to the melting layer. The algorithm is applied to 1 year of rainfall events that occurred over southeast England, and the results were validated using radiosonde data. After evaluating all possible combinations of polarimetric radar measurements, the algorithm achieves the best ML detection when combining VPs of ZH, ρHV and the gradient of the velocity (gradV), whereas, for QVPs, combining profiles of ZH, ρHV and ZDR produces the best results, regardless of the type of rain event. The root mean square error in the ML detection compared to radiosonde data is ∼200 m when using VPs and ∼250 m when using QVPs.

Hydrometeor classification of quasi-vertical profiles of polarimetric radar measurements using a top-down iterative hierarchical clustering method

Abstract. Correct, timely and meaningful interpretation of polarimetric weather radar observations requires an accurate understanding of hydrometeors and their associated microphysical processes along with well-developed techniques that automatize their recognition in both the spatial and temporal dimensions of the data. This study presents a novel technique for identifying different types of hydrometeors from quasi-vertical profiles (QVPs). In this new technique, the hydrometeor types are identified as clusters belonging to a hierarchical structure. The number of different hydrometeor types in the data is not predefined, and the method obtains the optimal number of clusters through a recursive process. The optimal clustering is then used to label the original data. Initial results using observations from the National Centre for Atmospheric Science (NCAS) X-band dual-polarization Doppler weather radar (NXPol) show that the technique provides stable and consistent results. Comparison with available airborne in situ measurements also indicates the value of this novel method for providing a physical delineation of radar observations. Although this demonstration uses NXPol data, the technique is generally applicable to similar multivariate data from other radar observations.

Characteristics of the Bright Band Based on Quasi-Vertical Profiles of Polarimetric Observations from an S-Band Weather Radar Network

Bright band (BB) characteristics obtained via dual-polarization weather radars elucidate thermodynamic and microphysical processes within precipitation systems. This study identified BB using morphological features from quasi-vertical profiles (QVPs) of polarimetric observations, and their geometric, thermodynamic, and polarimetric characteristics were statistically examined using nine operational S-band weather radars in South Korea. For comparable analysis among weather radars in the network, the calibration biases in reflectivity (ZH) and differential reflectivity (ZDR) were corrected based on self-consistency. The cross-correlation coefficient (ρHV) bias in the weak echo regions was corrected using the signal-to-noise ratio (SNR). First, we analyzed the heights of BBPEAK derived from the ZH as a function of season and compared the heights of BBPEAK derived from the ZH, ZDR, and ρHV. The heights of BBPEAK were highest in the summer season when the surface temperature was high. However, they showed distinct differences depending on the location (e.g., latitude) within the radar network, even in the same season. The height where the size of melting particles was at a maximum (BBPEAK from the ZH) was above that where the oblateness of these particles maximized (BBPEAK from ZDR). The height at which the inhomogeneity of hydometeors was at maximum (BBPEAK from the ρHV) was also below that of BBPEAK from the ZH. Second, BB thickness and relative position of BBPEAK were investigated to characterize the geometric structure of the BBs. The BB thickness increased as the ZH at BBBOTTOM increased, which indicated that large snowflakes melt more slowly than small snowflakes. The geometrical structure of the BBs was asymmetric, since the melting particles spent more time forming the thin shell of meltwater around them, and they rapidly collapsed to form a raindrop at the final stage of melting. Third, the heights of BBTOP, BBPEAK, and BBBOTTOM were compared with the zero-isotherm heights. The dry-temperature zero-isotherm heights were between BBTOP and BBBOTTOM, while the wet-bulb temperature zero-isotherm heights were close to the height of BBPEAK. Finally, we examined the polarimetric observations to understand the involved microphysical processes. The correlation among ZH at BBTOP, BBPEAK, and BBBOTTOM was high (>0.94), and the ZDR at BBBOTTOM was high when the BB’s intensity was strong. This proved that the size and concentration of snowflakes above the BB influence the size and concentration of raindrops below the BB. There was no depression in the ρHV for a weak BB. Finally, the mean profile of the ZH and ZDR depended on the ZH at BBBOTTOM. In conclusion, the growth process of snowflakes above the BB controls polarimetric observations of BB.

A Polarimetric Radar Analysis of Ice Microphysical Processes in Melting Layers of Winter Storms Using S-Band Quasi-Vertical Profiles

AbstractQuasi-vertical profiles (QVPs) obtained from a database of U.S. WSR-88D data are used to document polarimetric characteristics of the melting layer (ML) in cold-season storms with high vertical resolution and accuracy. A polarimetric technique to define the top and bottom of the ML is first introduced. Using the QVPs, statistical relationships are developed to gain insight into the evolution of microphysical processes above, within, and below the ML, leading to a statistical polarimetric model of the ML that reveals characteristics that reflectivity data alone are not able to provide, particularly in regions of weak reflectivity factor at horizontal polarization ZH. QVP ML statistics are examined for two regimes in the ML data: ZH ≥ 20 dBZ and ZH < 20 dBZ. Regions of ZH ≥ 20 dBZ indicate locations of MLs collocated with enhanced differential reflectivity ZDR and reduced copolar correlation coefficient ρhv, while for ZH < 20 dBZ a well-defined ML is difficult to discern using ZH alone. Evidence of large ZDR up to 4 dB, backscatter differential phase δ up to 8°, and low ρhv down to 0.80 associated with lower ZH (from −10 to 20 dBZ) in the ML is observed when pristine, nonaggregated ice falls through it. Positive correlation is documented between maximum specific differential phase KDP and maximum ZH in the ML; these are the first QVP observations of KDP in MLs documented at S band. Negative correlation occurs between minimum ρhv in the ML and ML depth and between minimum ρhv in the ML and the corresponding enhancement of ZH (ΔZH = ZHmax − ZHrain).

Polarimetric Radar Variables in the Layers of Melting and Dendritic Growth at X Band—Implications for a Nowcasting Strategy in Stratiform Rain

AbstractTime series of quasi–vertical profiles (QVPs) from 52 stratiform precipitation events observed with the polarimetric X-band radar in Bonn, Germany (BoXPol), between 2013 and 2016 have been statistically analyzed to infer microphysical processes shaping the dendritic-growth-layer (DGL) and melting-layer (ML) signatures including surface rainfall. Specific differential phase KDP in the ML shows an average correlation of 0.65 with surface rainfall for these cases. Radar reflectivity decreases below the ML by about 2 dB on average while differential reflectivity ZDR is hardly affected, which suggests rain evaporation as the dominating effect. Estimated ice water content or snow water equivalent precipitation rate S in the DGL is correlated with surface rain rates with lead times of 30 min and longer, which opens a pathway for radar-based nowcasting of stratiform precipitation tendencies. Trajectories of snow generated aloft down to the surface are constructed from wind profiles derived both from the nearest radiosounding and radar-based velocity azimuth displays (VAD) to narrow down the location at which the DGL-generated snow reaches the surface as rain. The lagged correlation analysis between KDP in the DGL and reflectivity ZH at that location demonstrates the superiority of the VAD information. Correlation coefficients up to 0.80 with lead times up to 120 min provide a proof of concept for future nowcasting applications that are based on DGL monitoring. Statistical relations found in our QVP dataset provide a database for estimating intrinsic polarimetric variables from the usual azimuth and elevation scans within and in the vicinity of the ML.

Sensitivity of Forward-Modeled Bending Angles to Vertical Interpolation of Refractivity for Radio Occultation Data Assimilation

Assimilation of radio occultation (RO) observations into numerical weather prediction (NWP) models has improved forecasts, where RO is typically one of the top five observational systems contributing to forecast accuracy. By measuring the phase delay of radio waves traversing Earth’s atmosphere between global positioning system (GPS) and low-Earth orbiting satellites, RO obtains quasi-vertical profiles of bending angles (BA) of the radio waves’ trajectories. BA are the RO observation most often assimilated into NWP models, but since they are not computed or analyzed in the models, they must be computed from model data using a forward model. First, model refractivity N is computed from variables specified on the model’s vertical levels, then using the Abel integral, BA are computed from N. The forward model requires vertical differentiation of N, and accurate differentiation requires vertical interpolation of N between model levels. The interpolation results in errors, which then propagate through the forward model to produce BA errors. In this study, we investigate the sensitivity of forward-modeled BA to five different methods of vertical interpolation of N between model levels to determine a method that minimizes interpolation errors. We use RO-observed N to isolate the interpolation errors and determine an accurate method that can be applied to any NWP model. Of the five methods investigated, the log-spline interpolation reduces N and BA errors the most, regardless of the vertical resolution of the model grid.

Convective Boundary Layer Depth Estimation from S-Band Dual-Polarization Radar

AbstractThis study investigates Bragg scatter signatures in dual-polarization radar observations, which are defined by low differential reflectivity values, as a proxy for convective boundary layer (CBL) depth. Using data from the WSR-88D in Twin Lakes, Oklahoma (KTLX), local minima in quasi-vertical profiles of are found to provide a reasonable estimate of CBL depth when compared with depth estimates from upper-air soundings from Norman, Oklahoma (KOUN), during 2014. The 243 Bragg scatter and upper-air sounding CBL depth estimates have a correlation of 0.90 and an RMSE of 254 m. Using Bragg scatter as a proxy for CBL depth was expanded to other seasons and locations—performing well in Wilmington, Ohio; Fairbanks, Alaska; Tucson, Arizona; Minneapolis, Minnesota; Albany, New York; Portland, Oregon; and Tampa, Florida—showing its potential usefulness in monitoring CBL depth throughout the year in a variety of geographic locations and meteorological conditions.

A Polarimetric Analysis of Ice Microphysical Processes in Snow, Using Quasi-Vertical Profiles

AbstractThis study implements a new quasi-vertical profile (QVP) methodology to investigate the microphysical evolution and significance of intriguing winter polarimetric signatures and their statistical correlations. QVPs of transitional stratiform and pure snow precipitation are analyzed using WSR-88D S-band data, alongside their corresponding environmental thermodynamic High-Resolution Rapid Refresh model analyses. QVPs of KDP and ZDR are implemented to demonstrate their value in interpreting elevated ice processes. Several fascinating and repetitive signatures are observed in the QVPs for differential reflectivity ZDR and specific differential phase KDP, in the dendritic growth layer (DGL), and at the tops of clouds. The most striking feature is maximum ZDR (up to 6 dB) in the DGL occurring near the −10-dBZ ZH contour within low KDP and during shallower and warmer cloud tops. Conversely, maximum KDP (up to 0.3° km−1) in the DGL occurs within low ZDR and during taller and colder cloud tops. Essentially, ZDR and KDP in the DGL are anticorrelated and strongly depend on cloud-top temperature. Analyses also show correlations indicating larger ZDR within lower ZH in the DGL and larger KDP within greater ZH in the DGL. The high-ZDR regions are likely dominated by growth of a mixture of highly oblate dendrites and/or hexagonal plates, or prolate needles. Regions of high KDP are expected to be overwhelmed with snow aggregates and crystals with irregular or nearly spherical shapes, seeded at cloud tops. Furthermore, QVP indications of hexagonal plate crystals within the DGL are verified using in situ microphysical measurements, demonstrating the reliability of QVPs in evaluating ice microphysics in upper regions of winter clouds.

Polarimetric Radar Observations of Precipitation Type and Rate from the 2–3 March 2014 Winter Storm in Oklahoma and Arkansas

AbstractA powerful winter storm affected the south-central United States in early March 2014, accompanied by elevated convective cells with hail and high rates of sleet, freezing rain, and snow. During portions of the event the thermal profile exhibited a shallow surface cold layer and warm, unstable air aloft. Precipitation falling into the cold layer refroze into ice pellets and was accompanied by a polarimetric refreezing signature and numerous crowdsourced surface ice pellet reports. Quasi-vertical profiles of the polarimetric variables indicated an enhanced reflectivity factor ZHH below the melting layer bright band and enhanced low-level differential reflectivity ZDR values coincident with surface ice pellet reports. Freezing rain rate was highest in areas with high ZHH and specific differential phase KDP values at low levels. High snow rates were most closely associated with 1- and 1.5-km ZHH values, though KDP and ZDR also appeared to show some ability to distinguish high snow rate. Numerous elevated convective cells contained rotating updrafts that appeared to contribute to storm longevity and intensity. Most contained well-defined ZDR maxima or columns and relatively high base-scan ZDR values. Several contained polarimetric signatures consistent with heavy mixed-phase precipitation and hail; social media reports indicated that large hail was produced by some of the storms.

Quasi-Vertical Profiles—A New Way to Look at Polarimetric Radar Data

AbstractA novel methodology is introduced for processing and presenting polarimetric data collected by weather surveillance radars. It involves azimuthal averaging of radar reflectivity Z, differential reflectivity ZDR, cross-correlation coefficient ρhv, and differential phase ΦDP at high antenna elevation, and presenting resulting quasi-vertical profiles (QVPs) in a height-versus-time format. Multiple examples of QVPs retrieved from the data collected by S-, C-, and X-band dual-polarization radars at elevations ranging from 6.4° to 28° illustrate advantages of the QVP technique. The benefits include an ability to examine the temporal evolution of microphysical processes governing precipitation production and to compare polarimetric data obtained from the scanning surveillance weather radars with observations made by vertically looking remote sensors, such as wind profilers, lidars, radiometers, cloud radars, and radars operating on spaceborne and airborne platforms. Continuous monitoring of the melting layer and the layer of dendritic growth with high vertical resolution, and the possible opportunity to discriminate between the processes of snow aggregation and riming, constitute other potential benefits of the suggested methodology.

Resist profile control of image reversal process in contact lithography

Image reversal with AZ 5214 resist is characterized by contact lithography, using a statistical experimental approach. Overhanging and quasi vertical profiles which are of particular interest in GaAs technology (lift-off and anisotropic etching processing) are obtained simply by varying the UV3 exposure dose and submicron features are achieved with a low dispersion (3 σ = 0.12 μm) using a thick resist.