Reducing Errors in Velocity–Azimuth Display (VAD) Wind and Deformation Retrievals from Airborne Doppler Radars in Convective Environments

2020 ◽  
Vol 37 (12) ◽  
pp. 2251-2266
Author(s):  
Charles N. Helms ◽  
Matthew L. Walker McLinden ◽  
Gerald M. Heymsfield ◽  
Stephen R. Guimond

AbstractThe present study describes methods to reduce the uncertainty of velocity–azimuth display (VAD) wind and deformation retrievals from downward-pointing, conically scanning, airborne Doppler radars. These retrievals have important applications in data assimilation and real-time data processing. Several error sources for VAD retrievals are considered here, including violations to the underlying wind field assumptions, Doppler velocity noise, data gaps, temporal variability, and the spatial weighting function of the VAD retrieval. Specific to airborne VAD retrievals, we also consider errors produced due to the radar scans occurring while the instrument platform is in motion. While VAD retrievals are typically performed using data from a single antenna revolution, other strategies for selecting data can be used to reduce retrieval errors. Four such data selection strategies for airborne VAD retrievals are evaluated here with respect to their effects on the errors. These methods are evaluated using the second hurricane nature run numerical simulation, analytic wind fields, and observed Doppler radar radial velocities. The proposed methods are shown to reduce the median absolute error of the VAD wind retrievals, especially in the vicinity of deep convection embedded in stratiform precipitation. The median absolute error due to wind field assumption violations for the along-track and for the across-track wind is reduced from 0.36 to 0.08 m s−1 and from 0.35 to 0.24 m s−1, respectively. Although the study focuses on Doppler radars, the results are equally applicable to conically scanning Doppler lidars as well.

2014 ◽  
Vol 142 (2) ◽  
pp. 573-589 ◽  
Author(s):  
Wen-Chau Lee ◽  
Xiaowen Tang ◽  
Ben J.-D. Jou

Abstract The concept and mathematical framework of the distance velocity–azimuth display (DVAD) methodology is presented. DVAD uses rVd (Doppler velocity scaled by the distance from the radar to a gate, r) as the basis to display, interpret, and extract information from single Doppler radar observations. Both linear and nonlinear wind fields can be represented by the same Cartesian polynomial with different orders. DVAD is mathematically concise and superior to the velocity–azimuth display (VAD) in interpreting and deducing flow characteristics. The rVd pattern of a two-dimensional linear wind field is exclusively in the form of a bivariate quadratic equation representing conic sections (e.g., ellipse, parabola, and hyperbola) centered at the radar depending only on divergence and deformation. The presence of a constant background flow translates the conic sections to a different origin away from the radar. It is possible to graphically estimate the characteristics of a linear wind field from the conical sections without performing a VAD analysis. DVAD analysis can deduce quantitative flow characteristics by a least squares fitting and/or a derivative method, and is a natural way to account for nonlinearity. The rVd pattern behaves similar to a type of velocity potential in fluid mechanics where ∇(rVd) is a proxy of the true wind vector and is used to estimate the general flow pattern in the vicinity of the radar.


2015 ◽  
Vol 54 (7) ◽  
pp. 1538-1555
Author(s):  
Xiaowen Tang ◽  
Wen-Chau Lee ◽  
Yuan Wang

AbstractThe application of the distance velocity azimuth display (DVAD) method to the retrieval of vertical wind profiles from single-Doppler radar observations is presented in this study. It was shown that Doppler velocity observations at a constant altitude can be expressed as a single polynomial function for both linear and nonlinear wind fields in DVAD. Only a one-step least squares fitting of a polynomial function is required to obtain the vertical wind profile of a real wind field. The mathematic formulation of DVAD results in two advantages over the traditional nonlinear VAD method used for the nonlinear analysis of single-Doppler observations. First, the requirement of only one-step least squares fitting leads to robust performance when Doppler velocity observations are contaminated by unevenly distributed data noise and voids. Second, the degree of nonlinearity to properly represent a real wind field can be directly estimated in DVAD instead of being empirically determined in the traditional method. A proper nonlinear wind model for approximating the real wind field can be objectively derived using the DVAD method. The merits of DVAD as a quantitative single-Doppler analysis method were compared with the traditional method using both idealized and real datasets. Results show that the simplicity and robust performance of DVAD make it a good candidate for single-Doppler retrieval in operational use.


2013 ◽  
Vol 30 (12) ◽  
pp. 2768-2787 ◽  
Author(s):  
Alessandro Battaglia ◽  
Simone Tanelli ◽  
Pavlos Kollias

Abstract Spaceborne Doppler radars have the potential to provide key missing observations of convective vertical air motions especially over the tropical oceans. Such measurements can improve understanding of the role of tropical convection in vertical energy transport and its interaction with the environment. Several millimeter wavelength Doppler radar concepts have been proposed since the 1990s. The Earth Clouds, Aerosols, and Radiation Explorer (EarthCARE) Cloud Profiling Radar (CPR) will be the first Dopplerized atmospheric radar in space but has not been optimized for Doppler measurements in deep convective clouds. The key challenge that constrains the CPR performance in convective clouds is the range–Doppler dilemma. Polarization diversity (PD) offers a solution to this constraint by decoupling the coherency (Doppler) requirement from the unambiguous range requirement. Careful modeling of the radar signal depolarization and its impact on radar receiver channel cross talk is needed to accurately assess the performance of the PD approach. The end-to-end simulator presented in this work allows reproduction of the signal sensed by a Doppler radar equipped with polarization diversity when overpassing realistic three-dimensional convective cells, with all relevant cross-talk sources accounted for. The notional study highlights that multiple scattering is the primary source of cross talk, highly detrimental for millimeter Doppler velocity accuracy. The ambitious scientific requirement of 1 m s−1 accuracy at 500-m integration for reflectivities above −15 dBZ are within reach for a W-band radar with a 2.5-m antenna with optimal values of the pulse-pair interval between 20 and 30 μs but only once multiple scattering and ghost-contaminated regions are screened out. The identification of such areas is key for Doppler accuracies and can be achieved by employing an interlaced pulse-pair mode that measures the cross and the copolar reflectivities. To mitigate the impact of attenuation and multiple scattering, the Ka band has been considered as either alternative or additional to the W band. However, a Ka system produces worse Doppler performances than a W-band system with the same 2.5-m antenna size. Furthermore, in deep convection it results in similar levels of multiple scattering and therefore it does not increase significantly the depth of penetration. In addition, the larger footprint causes stronger nonuniform beam-filling effects. One advantage of the Ka-band option is the larger Nyquist velocity that tends to reduce the Doppler accuracies. More significant benefits are derived from the Ka band when observing precipitation not as intense as the deep convection is considered here. This study demonstrates that polarization diversity indeed represents a very promising methodology capable of significantly reducing aliasing and Doppler moment estimate errors, two main error sources for Doppler velocity estimates in deep convective systems and a key step to achieving typical mission requirements for convection-oriented millimeter radar-based spaceborne missions.


2007 ◽  
Vol 135 (1) ◽  
pp. 93-117 ◽  
Author(s):  
John R. Stonitsch ◽  
Paul M. Markowski

Abstract Dual-Doppler observations acquired by a network of mobile radars deployed in the Oklahoma panhandle on 3 June 2002 are used to document the kinematic structure and evolution of a front. The data were collected during the International H2O Project on a mission to study the initiation of deep convection. Synchronized scanning allowed for the synthesis of three-dimensional wind fields for nearly 5.5 h of the 1557–0000 UTC period. The front initially moved southward as a cold front, stalled, and later retreated northward as a warm front. Deep convection failed to be initiated along the front. In situ thermodynamic measurements obtained by a mobile mesonet were used to document changes in the density gradient at the surface. This paper examines the relationships among the changes in baroclinity, the thermally direct frontal circulation, updraft intensity, alongfront updraft variability, and the intensity of vortices along the front. Increases in the front-normal density gradient tended to be associated with increases in the thermally direct frontal circulation, as expected. Increases in the front-normal density gradient were also associated with an increase in the tilt of the frontal updraft as well as an increase in the contiguity of the updraft along the front, termed the “slabularity.” During periods when the front-normal density gradient and associated thermally direct frontal circulation were weak, the kinematic fields were dominated by boundary layer convection and the slabularity of the front was reduced. Intensification of the front-normal density gradient was accompanied by an increase in the horizontal wind shear and the intensity of vortices that were observed along the front. The vortices modulated the vertical velocity field along the front and therefore the slabularity, too. Thus, although the slabularity was a strong function of the strength of the thermally direct frontal circulation, the slabularity appeared to be modified by vortices in complex ways. Possible implications of the observations for convection initiation are also discussed, particularly with respect to updraft tilt and slabularity.


2015 ◽  
Vol 54 (8) ◽  
pp. 1792-1808 ◽  
Author(s):  
Lin Tian ◽  
Gerald M. Heymsfield ◽  
Anthony C. Didlake ◽  
Stephen Guimond ◽  
Lihua Li

AbstractThe velocity–azimuth display (VAD) analysis technique established for ground-based scanning radar is applied to the NASA High-Altitude Imaging Wind and Rain Airborne Profiler (HIWRAP). The VAD technique provides a mean vertical profile of the horizontal winds for each complete conical scan of the HIWRAP radar. One advantage of this technique is that it has shown great value for data assimilation and for operational forecasts. Another advantage is that it is computationally inexpensive, which makes it suitable for real-time retrievals. The VAD analysis has been applied to the HIWRAP data collected during NASA’s Genesis and Rapid Intensification Processes (GRIP) mission. The traditional dual-Doppler analysis for deriving wind fields in the nadir plane is also presented and is compared with the VAD analysis. The results show that the along-track winds from the VAD technique and dual-Doppler analysis agree in general. The VAD horizontal winds capture the mean vortex structure of two tropical cyclones, and they are in general agreement with winds from nearby dropsondes. Several assumptions are made for the VAD technique. These assumptions include a stationary platform for each HIWRAP scan and constant vertical velocity of the hydrometeors along each complete scan. As a result, the VAD technique can produce appreciable errors in regions of deep convection such as the eyewall, whereas in stratiform regions the retrieval errors are minimal. Despite these errors, the VAD technique can still adequately capture the larger-scale structure of the hurricane vortex given a sufficient number of flight passes over the storm.


2008 ◽  
Vol 25 (11) ◽  
pp. 1939-1954 ◽  
Author(s):  
Michel Chong ◽  
Nabil Lamrani ◽  
Martin Hagen

Abstract The problem of sidelobe contamination of bistatic apparent Doppler velocity measurements involved in a bistatic Doppler radar network is examined. So far in the context of 3D wind field analysis, by combining a traditional Doppler radar with one or more bistatic receivers, identification and hence removal of regions of high degrees of contamination were necessarily crucial steps to obtaining reliable wind fields. This study proposes an alternative solution to the forced rejection of bistatic Doppler data suspected to be contaminated by sidelobe echoes, on the basis of restoring the nonmeasured “actual” (i.e., noncontaminated) bistatic Doppler velocity from both monostatic radar and bistatic receiver measurements. The correction method is based on a modeled expression of the observed bistatic apparent Doppler velocity defined as the reflectivity-weighted average of actual Doppler velocity of particles within individual volume samples, including the antenna gain pattern of both transmitting and receiving radars. The searched actual Doppler velocity is a solution of an underdetermined inverse problem that can be handled as a constrained linear inversion problem, through a variational least squares analysis method. The performances of the proposed method are analyzed, using simulated radar observations involving one remote receiver. An example of application to experimental data collected by the Deutsches Zentrum für Luft und Raumfahrt (DLR) bistatic Doppler radar network within a moderate precipitation system observed on 8 May 2000 in Germany is also presented. Pseudo-Doppler observations of a tropical squall-line system are used to quantify the effective improvement of the correction method on the bistatic Doppler velocity and hence the retrieved 3D wind field. Statistics of the differences are presented between observed and idealized (sidelobe free) velocity structures on the one hand, and corrected and idealized velocity structures on the other hand. Clearly shown is the very low level of the corrected minus idealized differences (mean and standard deviation) against the significantly high level of the observed minus idealized differences. As previously observed, maximum correction occurs in regions of potentially high gradients of reflectivity. It is also found that regions of low observed minus idealized differences remain unchanged after correction, which means that the sidelobe-correction method only acts on needed regions and does not introduce any artificial modification.


2012 ◽  
Vol 12 (15) ◽  
pp. 6679-6698 ◽  
Author(s):  
K. Schmidt ◽  
M. Hagen ◽  
H. Höller ◽  
E. Richard ◽  
H. Volkert

Abstract. The three-hour life-cycle of the isolated thunderstorm on 15 July 2007 during the Convective and Orographically-induced Precipitation Study (COPS) is documented in detail, with a special emphasis on the rapid development and mature phases. Remote sensing techniques as 5-min rapid scans from geostationary satellites, combined velocity retrievals from up to four Doppler-radars, the polarimetric determination of hydrometeors and spatio-temporal occurrences of lightning strokes are employed to arrive at a quantification of the physical parameters of this, during the COPS period, singular event. Inner cloud flow fields are available from radar multiple Doppler analyses at four consecutive times separated by 15 min-intervals. They contain horizontal winds of around 15 m s−1 and updrafts exceeding 5 m s−1, the latter collocated with lightning strokes. Reflectivity and polarimetric data indicate the existence of hail at the 2 km level around 14:40. Furthermore, polarimetric and Doppler radar variables indicate intense hydrometeor variability and turbulence corresponding to an enhanced variance of the retrieved 3-D wind fields. Profiles of flow and hydrometeor statistics over the entire cloud volume provide reference data for high-resolution numerical weather prediction runs in research mode. The study embarks from two movie-loops of geostationary satellite imagery (as Supplement), which provide an intuitive distinction of six phases making up the entire life-cycle of the thunderstorm. It concludes with a triple-image loop, juxtaposing a close-up of the cloud motion as seen by Meteosat, simulated brightness temperature (as a proxy for clouds seen by the infrared satellite channel), and a perspective view on the model generated system of cloud cells. The simulation suggests that several updrafts fed from a low level convergence line eventually removed the convective inhibition and set deep convection in motion. A shear line in the radial velocity relative to the Feldberg radar site shows good agreement beween observation and simulation, whereas the onset location of deep convection exhibits a horizontal discrepancy of 15 km. A quantitative schematic of the isolated thunderstorm synthesizes all retrieved characteristics.


2018 ◽  
Vol 35 (8) ◽  
pp. 1649-1663 ◽  
Author(s):  
Yu-Chieng Liou ◽  
Howard B. Bluestein ◽  
Michael M. French ◽  
Zachary B. Wienhoff

AbstractA three-dimensional data assimilation (3DVar) least squares–type single-Doppler velocity retrieval (SDVR) algorithm is utilized to retrieve the wind field of a tornadic supercell using data collected by a mobile, phased-array, Doppler radar [Mobile Weather Radar (MWR) 05XP] with very high temporal resolution (6 s). It is found that the cyclonic circulation in the hook-echo region can be successfully recovered by the SDVR algorithm. The quality of the SDVR analyses is evaluated by dual-Doppler syntheses using data collected by two mobile Doppler radars [Doppler on Wheels 6 and 7 (DOW6 and DOW7, respectively)]. A comparison between the SDVR analyses and dual-Doppler syntheses confirms the conclusion reached by an earlier theoretical analysis that because of the temporally discrete nature of the radar data, the wind speed retrieved by single-Doppler radar is always underestimated, and this underestimate occurs more significantly for the azimuthal (crossbeam) wind component than for the radial (along beam) component. However, the underestimate can be mitigated by increasing the radar data temporal resolution. When the radar data are collected at a sufficiently high rate, the azimuthal wind component may be overestimated. Even with data from a rapid scan, phased-array, Doppler radar, our study indicates that it is still necessary to calculate the SDVR in an optimal moving frame of reference. Finally, the SDVR algorithm’s robustness is demonstrated. Even with a temporal resolution (2 min) much lower than that of the phased-array radar, the cyclonic flow structure in the hook-echo region can still be retrieved through SDVR using data observed by DOW6 or DOW7, although a difference in the retrieved fields does exist. A further analysis indicates that this difference is caused by the location of the radars.


2018 ◽  
Vol 75 (1) ◽  
pp. 143-162 ◽  
Author(s):  
Udai Shimada ◽  
Masahiro Sawada ◽  
Hiroyuki Yamada

A ground-based Doppler radar observed the rapid intensification (RI) of Typhoon Goni (2015) for 24 h immediately after it completed an eyewall replacement cycle. Goni’s RI processes were examined by using radar reflectivity and wind fields retrieved by the ground-based velocity track display (GBVTD) technique. The maximum wind at 2-km altitude increased by 30 m s−1 during the first 6 h of RI, and it further increased by 20 m s−1 during the subsequent 12 h. Around the onset of RI, relatively strong outflow (>2 m s−1) was present both inside and outside the radius of maximum wind (RMW) above the boundary layer (BL), suggesting the existence of supergradient flow in and just above the BL. Despite this outflow, angular momentum increased inside the RMW. The low-level RMW contracted rapidly from 50 to 33 km, causing the RMW to slope greatly outward with height. The radius of maximum reflectivity was a few kilometers inside the RMW. A budget analysis of absolute angular momentum showed that the outflow contributed to the contraction of the tangential wind field. During RI, eyewall convection was enhanced, and a well-defined eye appeared. The low-level outflow changed into inflow immediately outside the RMW. Then the tangential wind field and high inertial stability region expanded radially outward, followed by the formation of an outer reflectivity maximum at twice the RMW. The contraction speed of the low-level RMW slowed down.


2008 ◽  
Vol 136 (3) ◽  
pp. 995-1012 ◽  
Author(s):  
Ben Jong-Dao Jou ◽  
Wen-Chau Lee ◽  
Su-Ping Liu ◽  
Yu-Cheng Kao

Abstract The primary circulation of atmospheric vortices, such as tropical cyclones and tornadoes, can be estimated from single-Doppler radar observations using the ground-based velocity track display (GBVTD) algorithm. The GBVTD algorithm has limitations in the following four areas: 1) distortion in the retrieved asymmetric wind fields, 2) a limited analysis domain, 3) the inability to resolve the cross-beam component of the mean wind, and 4) the inability to separate the asymmetric tangential and radial winds. This paper presents the generalized velocity track display (GVTD) algorithm, which eliminates the first two limitations inherent in the GBVTD technique and demonstrates the possibility of subjectively estimating the mean wind vector when its signature is visible beyond the influence of the vortex circulation. In this new paradigm, the GVTD algorithm fits the atmospheric vortex circulation to a new variable VdD/RT in a linear azimuth angle (θ′), rather than the Doppler velocity Vd in a nonlinear angle (ψ), which is used in GBVTD. Key vortex kinematic structures (e.g., mean wind, axisymmetric tangential wind, etc.) in the VdD/RT space simplify the interpretation of the radar signature and eliminate the geometric distortion inherent in the Vd display. This is a significant improvement in diagnosing vortex structures in both operations and research. The advantages of using VdD/RT are illustrated using analytical atmospheric vortices, and the properties are compared with GBVTD. The characteristics of the VdD/RT display of Typhoon Gladys (1994) can be approximated by a constant mean wind plus an axisymmetric vortex.


Sign in / Sign up

Export Citation Format

Share Document