Research article 18 May 2021
Research article  18 May 2021
Comparison of singleDoppler and multipleDoppler wind retrievals in Hurricane Matthew (2016)
 Department of Atmospheric Science, Colorado State University, Fort Collins, Colorado, USA
 Department of Atmospheric Science, Colorado State University, Fort Collins, Colorado, USA
Correspondence: TingYu Cha (tingyu@colostate.edu)
Hide author detailsCorrespondence: TingYu Cha (tingyu@colostate.edu)
Hurricane Matthew (2016) was observed by the groundbased polarimetric Next Generation Weather Radar (NEXRAD) in Miami (KAMX) and the National Oceanic and Atmospheric Administration WP3D (NOAA P3) airborne tail Doppler radar near the coast of the southeastern United States for several hours, providing a novel opportunity to evaluate and compare single and multipleDoppler wind retrieval techniques for tropical cyclone flows. The generalized velocity track display (GVTD) technique can retrieve a subset of the wind field from a single groundbased Doppler radar under the assumption of nearly axisymmetric rotational wind, but it has been shown to have errors from the aliasing of unresolved wind components. An improved technique that mitigates errors due to storm motion is derived in this study, although some spatial aliasing remains due to limited information content from the singleDoppler measurements. A splinebased variational wind retrieval technique called SAMURAI can retrieve the full threedimensional wind field from airborne radar fore–aft pseudodualDoppler scanning, but it has been shown to have errors due to temporal aliasing from the nonsimultaneous Doppler measurements. A comparison between the two techniques shows that the axisymmetric tangential winds are generally comparable between the two techniques, and the improved GVTD technique improves the accuracy of the retrieval. Fourier decomposition of asymmetric kinematic and convective structure shows more discrepancies due to spatial and temporal aliasing in the retrievals. The strengths and weaknesses of each technique for studying tropical cyclone structure are discussed and suggest that complementary information can be retrieved from both single and dualDoppler retrievals. Future improvements to the asymmetric flow assumptions in singleDoppler analysis and steadystate assumptions in pseudodualDoppler analysis are required to reconcile differences in retrieved tropical cyclone structure.
Doppler radar can provide highresolution wind measurements within tropical cyclones (TCs), but the measurement is limited to the projection of the wind along the radial direction of the radar beam. Therefore, wind retrieval techniques are required in order to identify the convective and kinematic structure of TCs from either single or multipleDoppler observations. While multipleDoppler retrievals are generally superior for deriving threedimensional winds, measurements from two or more radars are not generally available and are often not simultaneous. In addition to the presence of an airborne Doppler radar with fore–aft capability or multiple radars with sufficient range and geometry around the TC, a steadystate assumption during the Doppler radar observation period is required to synthesize the wind fields into one snapshot in time. The steadystate assumption is less severe for singleDoppler wind retrievals, but more assumptions about the unresolved components of the flow are required. Previous studies have shown the intercomparison of dualDoppler wind fields from two orthogonal flight legs and a groundbased tworadar network (Jorgensen et al., 1983; Hildebrand and Mueller, 1985). Several other studies have investigated both single and multiDoppler techniques for retrieving TC wind fields (Lee et al., 1994; Crum et al., 1998; Reasor et al., 2000; Lee et al., 1999; Jou et al., 2008; Bell et al., 2012), but the strengths and weaknesses of different techniques have not been compared and addressed fully. In this study, groundbased singleDoppler and airborne dualDoppler observations simultaneously sampling Hurricane Matthew (2016) are analyzed to provide the first comprehensive comparison between groundbased singleDoppler and airborne multiDoppler wind retrieval techniques in a TC.
One of the primary dualDoppler platforms for TC studies is the airborne National Oceanic and Atmospheric Administration WP3D (NOAA P3) tail Doppler radar (TDR), which can obtain kinematic structure for storms well away from the US coast. The TDR has been used since the early 1980s for airborne radar data collection in hurricane reconnaissance and research missions. In early usage, the wind field was reconstructed using the pseudodualDoppler analysis of two flight legs that were perpendicular to each other known as the “L” pattern (Marks and Houze, 1984, 1987). The L pattern takes a period of 0.5 to 1 h to complete, such that some slowly evolving wind asymmetries can be deduced (Marks et al., 1992), but more rapidly evolving structure cannot be retrieved by this technique.
In Lee et al. (1994), the velocity track display (VTD) technique was proposed to retrieve the TC kinematic structure from a single airborne Doppler radar. The VTD technique takes advantage of the Doppler signatures of a vortex with a dipole pattern of approaching and receding velocities. The linear least squares method is utilized to fit Fourier basis functions onto the Doppler velocity, such that a subset of the wind field can be retrieved for each radius and altitude based on the Fourier coefficients. VTD retrievals from Hurricane Gloria (1985) were compared with pseudodualDoppler analysis constructed from two orthogonal flight legs under the assumption that the circulation was in a steady state over a 1–2 h period. Lee et al. (1994) concluded that spatial and temporal aliasing in the pseudodualDoppler analysis over long periods tended to create artificial higher wave numbers and reduce the wind maxima compared with the singleDoppler retrievals. The study was one of the first to quantify the accuracy of lowwavenumber TC asymmetries from multipleDoppler analysis due to the evolving weather.
The fore–aft scanning technique (FAST) was proposed by Hildebrand et al. (1986) in which the antenna alternately scans forward and then aft of the flight track within a few seconds, which could produce more accurate local wind estimations. This scanning methodology can gather data faster, mitigating some of the impact from weather evolution, and no longer requires the execution an Lshaped flight pattern (Gamache et al., 1995). The FAST approach produces horizontal winds with high accuracy, and it can yield vertical motion through the mass continuity equation or a multibeam technique with multiple aircraft (Jorgensen et al., 1996). The FAST scanning strategy has been used almost exclusively in recent years and has provided significant advances in our understanding of TC structure and dynamics (Lee et al., 2003; Houze, 2010).
Although dualDoppler observations can be used to assess snapshots of highresolution kinematic and convective structure, airborne reconnaissance and research missions are rare events in most of the countries impacted by TCs. The threedimensional airflow structure can also be retrieved from the dualDoppler observations when the system is detected by two groundbased radars. General sources of error in the intercomparison of groundbased and airborne dualDoppler observations include instrument effects, algorithm effects, and sampling effects (Hildebrand and Mueller, 1985). Instrument effects include the effects of attenuation, signaltonoise ratio, and the number of radar samples that could be caused by the radar processor design or measurement technique. These effects are likely to be most influential with marginal signaltonoise, but random velocity errors of up to 1 m s^{−1} are possible with many radar designs, including airborne radars (Hildebrand et al., 1994). Algorithm effects include the effects of the interpolation to the Cartesian grid, multiDoppler geometry, the solution method and its associated assumptions, and the derivation of the vertical velocities. Sampling effects include the effects of data spacing and density, geometry of flight tracks, temporal changes in the storm, advection, and data collection period. One of the longlasting problems is the length of time required for each flight leg with airborne Doppler radar (Ray and Stephenson, 1990). The temporal effects can degrade the analysis if the data collection takes too long. Jorgensen et al. (1983) quantitatively compared the wind fields of homogeneous precipitation derived from the two pseudoorthogonal flight legs and two groundbased dualDoppler observations. Their measurements showed agreement in the horizontal wind fields but small discrepancies in the vertical velocities of about 0.5–1 m s^{−1} for the airborne system and about 0.2 m s^{−1} for the groundbased system. The discrepancy was attributed to uncertainties in the pointing angle of the airborne system and a long data collection period.
Groundbased dualDoppler radar observations of TCs are usually limited to the observation of storms that happen to develop or move within the domain covered by the radars and extensive radar baselines (Jou et al., 1996). The range limitation of groundbased radar observations can restrict the operational exploitation of dualDoppler measurements due to the large spacing between coastal radars and limited dualDoppler lobes. As such, singleDoppler retrieval techniques are often required to estimate the vortex structure.
The groundbased velocity track display (GBVTD) was developed by Lee et al. (1999) as an extension of VTD for stationary radar scanning geometry. The GBVTD technique provides a new way to examine axisymmetric and asymmetric structures of a TC near landfall from single groundbased Doppler radar, and it has been successfully utilized in several studies (Lee et al., 2000; Lee and Bell, 2007; Zhao et al., 2008, 2012; Shimada et al., 2018). One limitation is that the radial distance between the radar and the storm center has to be large enough to sample the tangential component of the vortex circulation in order to minimize the geometric distortion. Additionally, the GBVTD technique cannot fully separate the asymmetric components of the tangential and radial wind or the mean environmental wind due to spatial aliasing. Extensions to the GBVTD technique have been developed to better resolve asymmetries when multiple radars are available (extended groundbased velocity track display, EGBVTD; Liou et al., 2006) or to resolve the mean wind when sufficient scatterers are near the radar (modified groundbased velocity track display, MGBVTD; Chen et al., 2013).
The generalized velocity track display (GVTD) (Jou et al., 2008) is a technique that improves upon GBVTD by introducing an aspect ratio calculated by multiplying the distance of each gate (D) by measured Doppler velocity (V_{d}) and then scaling by the distance between the radar and the TC center (R_{T}). Key vortex kinematic structures displayed in the ${V}_{\mathrm{d}}D/{R}_{\mathrm{T}}$ space simplify the interpretation of the radar signature and eliminate the geometric distortion inherited in the V_{d} space (Jou et al., 1996). GVTD expands ${V}_{\mathrm{d}}D/{R}_{\mathrm{T}}$ into Fourier coefficients in a linear coordinate (θ^{′}) rather than expanding V_{d} in a nonlinear coordinate (ψ^{′}) in GBVTD. The geometry and symbols of GVTD are displayed in Fig. 1. The retrieved wind field from GVTD is no longer limited by the analysis domain due to the required approximation of cos α in GBVTD (Eq. 5 in Lee et al., 1999), and the retrieved asymmetric structures are without distortion. The percentage errors of the retrieved wave number 2 and 3 asymmetries are negligible (<1 %) in general, which agrees well with the analytical solutions. Also, the GVTD formulation can be applied to the extensions of the velocity track display (VTD) techniques (e.g., GBVTDsimplex, Lee and Marks, 2000) to improve their performance. Therefore, it not only expands the capability of using groundbased Doppler radar data in TC forecasts but provides researchers with an opportunity to examine TC kinematic and some derived dynamic variables in detail (such as vertical velocity, angular momentum, and vertical vorticity).
Jou et al. (2008) tested the GVTD technique with idealized vortices and confirmed that the modification by the aspect ratio is beneficial to retrieve the vortex kinematic structure. Although the GVTD technique largely improves and extends the capability of GBVTD, several limitations of the singleDoppler measurement remain. In particular, the equations are still underdetermined and require a closure assumption in order to retrieve the asymmetric wind field.
The primary motivation of this study is to compare the strengths and weaknesses of the GVTD technique and an airborne dualDoppler windsynthesis analysis in a real case, which has not been done in previous studies. Hurricane Matthew (2016) was observed by the Next Generation Weather Radar (NEXRAD) radar in Miami (KAMX), concurrently with the National Oceanic and Atmospheric Administration WP3D (NOAA P3, hereafter P3) airborne tail Doppler radar when it approached the southeastern US. The KAMX radar has a larger data coverage, but the P3 has better spatial resolution due to its closer range. As airborne and groundbased radar collected data simultaneously, the case provides a unique opportunity to evaluate the wind retrievals. The datasets and methodology are described in Sect. 2. The improved algorithm of the GVTD is formulated in Sect. 3. In Sect. 4, the improved GVTD algorithm is applied to the NEXRAD data from Hurricane Matthew (2016) and compared with the retrieved winds from the dualDoppler analysis. A summary of our results and conclusions are presented in Sect. 5.
Hurricane Matthew (2016) was the first category 5 hurricane in the Atlantic Basin since 2007, and it caused widespread damage across its destructive path. When Matthew moved parallel to the east coast of Florida, it was observed simultaneously by the KAMX singleDoppler radar at a 5 min interval and the P3 TDR from 19:00 UTC on 6 October to 00:00 UTC on 7 October with four flight passes through the center during Matthew's eyewall replacement process.
The flight track of the P3 and the detecting range of the KAMX radar are displayed in Fig. 2a. All radar sweep files were initially processed using an automated quality control (QC) script using National Center for Atmospheric Research (NCAR) SoloII software (Bell et al., 2013) and then manually edited to unfold the Doppler velocity aliasing and remove the discontinuities and noise echoes. A coordinate transform and interpolation were applied to the KAMX radar fields from the original plan position indicators (PPIs) to the constantaltitude plan position indicators (CAPPIs) in Cartesian coordinates using Radx2Grid in the Lidar Radar Open Software Environment (LROSE) software (Bell, 2019). The gridded domain is 400 km×400 km with a horizontal grid spacing of 1 km and vertical grid spacing of 0.5 km. While the vertical resolution may be a bit fine, we focus on the horizontal structure which is appropriately resolved for the given sampling. The gridded data were further analyzed using the Vortex Objective Radar Tracking and Circulation (VORTRAC) software in LROSE to interpolate onto a cylindrical coordinate and obtain the kinematic structure by the improved GVTD algorithm formulated in Sect. 3.1.
The P3 was involved in the reconnaissance mission from 19:00 UTC on 6 October to 00:00 UTC on 7 October. The P3 was equipped with the TDR, which scanned in FAST mode in order to obtain pseudodualDoppler measurements. The TDR documented the intensification and weakening stages of Matthew's eyewall replacement cycle (ERC). The P3 flew four radial passes through the center of the TC, with each pass being 30 to 60 min apart. The time windows of the four passes are listed in Table 1, and the location of each pass is shown in Fig. 2a. The storm center and storm motion were both estimated from the Hurricane Research Division (HRD) aircraftderived dynamic center (Willoughby and Chelmow, 1982). The analysis track for dualDoppler analyses was linearly interpolated from each dynamic center using the derived storm motion. The dualDoppler analysis was synthesized with each of the P3 radial passes at 1 km horizontal spline nodal spacing and 0.5 km vertical nodal spacing using SAMURAI software (Bell et al., 2012) in LROSE, with a 4Δ× Gaussian filter in the horizontal and 2Δ× filter in the vertical applied. SAMURAI is a threedimensional variational data assimilation tool that uses a finite element basis to estimate the most likely state of the atmosphere given a set of observations. The nodal spacing of the finite elements should be smaller than the data spacing in order to accurately represent a spline function that can depict the spatial scales resolved by a given data sampling (e.g., Koch et al., 1983; Ooyama, 1987, 2002). For the P3 TDR in 2016, the data spacing is limited in the alongtrack direction to ∼1.4 km due to the rotation rate of the radar. With the chosen spline nodal spacing and Gaussian filter length the minimum resolved scale is ∼4 km in the horizontal, or approximately 2.85 times the alongtrack data spacing. Largerscale features such as low azimuthal wave number structures are well resolved by the analysis. Additional algorithm effects of SAMURAI have been tested in Bell et al. (2012) where the analysis has high fidelity to observations with low noise, with linear correlations of 0.99 and linear slope and bias values near one and zero, respectively. The analysis was initially done on a Cartesian coordinate and then interpolated onto a cylindrical coordinate with azimuthal resolution of 1^{∘} and radial resolution of 1 km. The four passes were analyzed in detail to examine the changes in kinematic structure with high spatial resolution over 4 h, with a particular focus on the first pass in this study.
The Doppler velocity in Jou et al. (2008) is decomposed into tangential, radial, and mean wind components, where storm motion is an implicit element in the mean wind component. The mean wind (V_{M}) is the horizontal average of the environmental flow at each altitude following the procedure proposed by Marks et al. (1992), which can be used to calculate the vertical wind shear. While the divergence of the environmental flow may not be zero as the environmental flow is not constant in the horizontal direction (Chan and Gray, 1982; Chan, 1984), we assume a horizontally homogeneous mean wind in the following derivation. This assumption is different from the original derivation in Jou et al. (2008) where the mean wind is a function of radius and height. Storm motion (U_{S}, V_{S}) is defined here as the deep layer motion vector over the whole vortex that varies only in time and does not vary with height or radius. The remaining terms in the GVTD formulation, namely the tangential (V_{T}) and radial (V_{R}) winds, are functions of radius, height, and azimuth to the storm center. Despite the fact that the environmental wind is an important factor to determine storm motion, the storm motion and environmental wind are not the same component due to vertical wind shear.
The storm motion and mean wind component were originally combined together in GVTD and affected the retrievals of the mean wind as well as the axisymmetric (wave number 0) tangential and radial winds. Notably, they do not have an influence on retrieving the phase and magnitude of the GVTD asymmetric components of tangential wind. However, when the storm moves fast or there is a high deviation in direction between the mean wind and storm motion, they can produce errors. To resolve this error, we rederive the GVTD technique and separate the storm motion from the mean wind component.
3.1 Mathematical formulation
Following the symbols and geometry utilized in the GBVTD (Lee et al., 1999) and GVTD (Jou et al., 2008) techniques, the addition of the storm motion to the geometry is illustrated in Fig. 1. We start with the horizontal projection of the Doppler velocity:
where ϕ is the elevation angle of the radar, θ_{d} is the mathematical angle of the radar measured from the east, θ_{M} is the direction of mean wind, and ψ is the angle composed by the measured radar beam to the radar and the measured radar beam to the storm center. Note that ${\widehat{V}}_{\mathrm{d}}$ neglects the contribution from the terminal velocity (v_{t}) and vertical velocity (w) (Eq. 2 in Lee et al., 1999). The contribution from w and v_{t} is small if the elevation angle of the radar beam is low ($<\mathrm{1}{}^{\circ}$). In this study, the storm motion (U_{S}, V_{S}) is added into the equation as an independent, known variable that projects onto the Doppler velocity, such that
where $\mathit{\psi}=\mathit{\theta}{\mathit{\theta}}_{\mathrm{d}}$. Rearranging Eq. (2), we obtain
The radar angle θ_{d} can be denoted as
Plugging Eq. (4) into Eq. (3) and approximating ${\widehat{V}}_{\mathrm{d}}/\mathrm{cos}\mathit{\varphi}$ with V_{d} (only valid when the elevation angle is low),
Rearranging Eq. (5) and letting $\mathit{\theta}{}^{\prime}=\mathit{\theta}{\mathit{\theta}}_{\mathrm{T}}$, we obtain
Decomposing ${V}_{\mathrm{d}}D/{R}_{\mathrm{T}}$, V_{T}, and V_{R} into Fourier components in the θ^{′} coordinates,
where A_{n} (V_{T}C_{n} and V_{R}C_{n}) and B_{n} (V_{T}S_{n} and V_{R}S_{n}) are the amplitude of the azimuthal wave number n cosine and sine components, as defined in Lee et al. (1999) and Jou et al. (2008), respectively. Substituting Eqs. (7), (8), and (9) into Eq. (6), we obtain the following expressions for the relation between the Fourier coefficients and wind components:
The Fourier coefficients can be rearranged to obtain each wind component of the vortex as follows:
Plugging Eq. (16) into Eq. (15),
Equations (15)–(19) correspond to Eqs. (16)–(20) in Jou et al. (2008) with the additional terms of storm motion in Eqs. (15)–(17). Equation (20) is an updated version of Eq. (15) to minimize the unknown terms after plugging in the V_{M}cos (θ_{T}−θ_{M}). One caveat of the V_{R}C_{0} updated form is that the axisymmetric radial wind cannot be derived when R = R_{T} because of the singular point. The derivation shows that the storm motion is aliased on the components of mean wind as well as wave number 0 tangential and wave number 0 radial wind. As the storm motion can be accurately estimated over successive radar volumes, the above equations can yield more accurate estimation of these wind components. However, the separation of storm motion and mean wind does not solve the underdetermined problem that the number of unknown variables is greater than the number of equations. We apply the same closure assumption as that used in Lee et al. (1999) and Jou et al. (2008): the asymmetric component of radial wind is much smaller than the asymmetric component of tangential wind, so the terms associated with radial wind asymmetries can be ignored. This closure assumption may not be applicable within the boundary layer or outflow layer where the radial wind asymmetries can be substantial. Future research is required to improve this closure assumption and retrieve the asymmetric radial wind.
With the above equations, the alongbeam component of the mean flow, axisymmetric (n=0) tangential wind, axisymmetric radial wind, and asymmetric tangential winds (n=1–2) can be retrieved by performing a linear least squares fit on ${V}_{\mathrm{d}}D/{R}_{\mathrm{T}}$ in a TCcentered cylindrical coordinate. All data within 1 km radiuswide annulus are included in the linear least squares fit. To deal with missing data in observational radar data and reduce the influence of outliers (Matejka and Srivastava, 1991), the truncation of the Fourier series follows Lee et al. (2000) (Table 2), which is consistent with the restriction of maximum allowable gap size in Lorsolo and Aksoy (2012). Lorsolo and Aksoy (2012) have shown that the maximum allowable gap size varies with the number of gaps and the noise. If more gaps are present in the signal, the maximum allowable gap size is greater than originally suggested in Lee et al. (2000). We allowed the maximum wave number up to wave number 2 in the retrieved tangential wind in this study to reduce retrieval errors.
3.2 GVTDsimplex center finding
The GVTD algorithm can be highly sensitive to the center location. Jou et al. (2008) showed that the uncertainty of the center location cannot exceed 5 km, or about 20 % of the radius of the maximum wind (RMW), in order to have a reasonable wind retrieval. There are several ways to identify TC center, such as the geometric center (Griffin et al., 1992), wind center (Wood and Brown, 1992), dynamic center (Willoughby and Chelmow, 1982), and vorticity center (Marks et al., 1992). The centers identified by different methods are not necessarily colocated, and the range of uncertainties is a few kilometers or more. As both vorticity centers estimated from the GVTD technique and dynamic centers derived from the aircraft reconnaissance were available, the comparison of the different centers is required in order to have a better result of wind retrieval.
The GBVTDsimplex algorithm is a method to identify TC vorticity centers using singleDoppler radar data developed by Lee and Marks (2000). The simplex center is found by maximizing the mean tangential wind within an axisymmetric TC with three operations on a simplex: reflection, contraction, and expansion (Lee and Marks, 2000; Harasti et al., 2004). The GBVTDsimplex algorithm reduces the uncertainties in estimating TC position and improves the quality of the GBVTDretrieved TC circulation. The deviation of the true centers to the GBVTDsimplex center in an idealized TC is approximately 340 m. In this study, as the GVTD technique has better wind field estimation, we conducted the GVTDsimplex method to estimate the centers following the GBVTDsimplex algorithm. By maximizing GVTDretrieved mean tangential wind, using the GVTD technique to estimate the TC vorticity center has higher accuracy than the GBVTD technique.
The GVTDsimplex method is performed as follows:

Doppler velocities on a CAPPI are interpolated onto a cylindrical grid with 1 km radial spacing centered at a given TC center.

The TC center possessing the highest mean tangential wind within an axisymmetric TC is found.

Three operations are used on a simplex to search for a new maximum or minimum in the field around the simplex: reflection, contraction, and expansion. We put the dynamic center from HRD as the first guess. The operation process would start from this point and find the circulation center.
To compare the performance of GVTDsimplex and GBVTDsimplex centers, the dynamic centers from HRD (Willoughby and Chelmow, 1982) are treated as the reference center because they are consistently centered at the geometric center of the storm over our analysis period compared to the simplex centers. The centers are interpolated from a few dynamic centers with a series of spline curves every 2 min, so the centers are connected into a continuous track. Figure 2b shows that the centers retrieved by the GVTDsimplex method are not consistent with that of the retrievals from the GBVTDsimplex method, which is due to the more accurate estimation of axisymmetric tangential wind. Although the GVTDsimplex centers follow more closely to the dynamic centers over the 35 h of observation, the location of the GVTDsimplex center is still variable, similar to GBVTDsimplex center (Lee and Marks, 2000) and not fully consistent with different radar retrievals (such as KAMX vs. KMLB, and KMLB vs. KJAX). As the dynamic centers are qualitatively and quantitatively better than the simplex centers in terms of consistency and are independent of assumptions in both the airborne dualDoppler and singleDoppler retrievals, our study utilized the dynamic centers to perform the GVTD technique and cylindrical decomposition of the airborne wind field.
4.1 Wave number 0 tangential wind retrieval
The dualDoppler analyses from the four aircraft passes into Hurricane Matthew are optimal for evaluating the performance of the GVTD technique because we can obtain all coefficients (Eqs. 16 to 20) from the Fourier decomposition of V_{T} and V_{R}, known storm motion, and mean wind components, which ensures the comparability of the wind field. To evaluate the improved algorithm, the wind fields from the dualDoppler analyses were resampled into Doppler velocity that would be observed by the KAMX radar. The subsequent analyses use the observations at the altitude of 4 km, so the groundbased 0.5^{∘} radar elevation beam can detect the TC inner core. The wind fields were then retrieved from the resampled velocities using the original and improved GVTD algorithms (Table 1). Figure 3a and b show the Doppler velocity observed by the KAMX radar and the Doppler velocity projected from the dualDoppler analysis using the 18:55–19:40 UTC aircraft data as an example. Figure 3c and d show the reflectivity field from the KAMX radar observation and dualDoppler analysis, respectively. The singleDoppler observations have much more missing data in the eye and the moat due to the reduced sensitivity to weak echoes at longer range from the radar. Additionally, the maximum Doppler velocity of single Doppler observations is not colocated with the maximum Doppler velocity of the resampled dualDoppler analysis. This discrepancy is believed to be due to temporal aliasing from the extended sampling period of the aircraft pass compared with the groundbased scanning. The issue of temporal aliasing of airborne Doppler analysis will be discussed in more detail in Sect. 4.2. Despite the aliasing, the dualDoppler analysis provides a relatively complete, consistent wind field and can reasonably be assumed as the “truth” for the purposes of algorithm evaluation.
Figure 4 displays the results from the improved GVTD algorithm compared to the original GVTD algorithm and the “true” wave number 0 (axisymmetric) tangential wind derived from the four dualDoppler analyses. The green dashed line (hereafter optimal solution) shows the improved GVTD retrieved axisymmetric tangential wind (V_{T}C_{0}) when all the terms are known. This optimal retrieval from the resampled ${V}_{\mathrm{d}}D/{R}_{\mathrm{T}}$ from the dualDoppler analysis using Eq. (17) has the least deviation from the truth (black line) in general. When the storm motion is removed from the V_{T}C_{0} retrieval (blue dashed line, hereafter optimal but without storm motion), the axisymmetric wind in the inner eyewall (between 10 and 30 km) does not differ significantly from the truth. As the radius increases, the storm motion term becomes more important to the retrieval of wave number 0 tangential wind. Between the radius of 55 to 65 km, the impact of neglecting the storm motion term can be up to 3–4 m s^{−1}. In the optimal retrieval, all terms in Eq. (17) are known, but in practice, the retrieved Fourier coefficients are underdetermined and a closure assumption is required to retrieve V_{T}C_{0}, typically by neglecting the crossbeam mean wind and asymmetric radial wind. The orange line (hereafter original GVTD) in Fig. 4 shows the result of the original GVTD algorithm invoking this closure assumption where only the B_{1} and B_{3} coefficients are used to retrieve V_{T}C_{0}. The red line (hereafter improved GVTD) shows the improvement to the GVTD algorithm by adding in the storm motion terms, while still invoking the necessary closure assumption. The deviations from the truth caused by the closure assumption occur at different radii depending on the aircraft pass, suggesting that the neglected terms have different magnitudes in different parts of the storm. For example, in pass 1 the deviations are most pronounced from 55 to 70 km radius, whereas in pass 4 the deviations are most pronounced from 15 to 25 km radius.
Rootmeansquare (RMS) differences in the above V_{T}C_{0} solutions and integrated perturbation pressure deficit (optimal, optimal but without storm motion, and original GVTD and improved GVTD algorithms) averaged over the four passes are shown in Table 3. For the RMS difference in V_{T}C_{0}, the optimal solution has the least deviation from the truth. Including the storm motion terms decreases the RMS differences by about 0.8 m s^{−1} in both optimal versus optimal without storm motion and improved GVTD versus original GVTD. Interestingly, the RMS difference in the improved GVTD algorithm has a similar magnitude to the optimal solution but without storm motion, suggesting that the influence of storm motion herein is similar to neglecting terms in Eq. (17).
The perturbation pressure deficit is integrated from r=10 to 70 km using the gradient wind balance equation (Lee et al., 2000). The integrated perturbation pressure deficit retrievals from different methods agree well with the truth (∼1 hPa RMS in general). Similar to the results of the RMS difference in V_{T}C_{0}, including the storm motion terms decreases the RMS differences by about 0.2 hPa. The RMS difference in the improved GVTD algorithm from the singleDoppler retrieval has the least deviation from the truth, suggesting that the perturbation pressure deficit derived from the singleDoppler observations has high fidelity.
The V_{T}C_{0} values retrieved by the original (orange dots) and improved (red dots) GVTD algorithm from the KAMX observations alone (Fig. 3a) are also shown in Fig. 4. The axisymmetric tangential winds retrieved from the KAMX radar alone are in generally good agreement with the winds from the dualDoppler analyses. As in the dualDoppler retrieval, the storm motion term has a relatively small impact on the retrieval in the inner eyewall because $R/{R}_{\mathrm{T}}$ is small. The deviation between the retrievals becomes larger beyond the radius of 40 km, but the improved GVTD algorithm is generally more consistent with the truth compared with the original GVTD algorithm. The RMS differences show a reduction of 0.4 m s^{−1} between the original and improved GVTD algorithms compared with the truth (Table 3).
A nonparametric Wilcoxon signedrank test is conducted to test the null hypothesis that two paired sets of the RMS differences derived from the original and improved GVTD algorithms are drawn from the same distribution. The RMS difference between the original and improved GVTD algorithms is statistically significant with a p value <0.001 using both the projected dualDoppler winds and singleDoppler velocities, indicating that we can reject the null hypothesis at the 1 % significance level (99 % confidence). The statistics suggest that the RMS differences distribution of wave number 0 tangential wind retrieved from the original GVTD algorithm are likely to be larger than those from the improved GVTD method. While it is a relatively small reduction in the RMS difference in the current case, the statistically significant difference in this algorithm error contributes to an overall reduction in the total error from instrument, algorithm, and sampling contributions. The comparison demonstrates that the inclusion of the storm motion has better consistency of the wave number 0 tangential wind intensity with the dualDoppler reference solution and that the improved GVTD technique can provide insightful information about TC kinematic structure from groundbased singleDoppler radar data.
4.2 Asymmetric wind retrievals
Both Lee et al. (1999) and Jou et al. (2008) have tested the GBVTD and GVTD techniques with analytic datasets, respectively, and confirmed the accuracy of these two methods. Murillo et al. (2011) further assessed the ability of GBVTD to retrieve lowwavenumber wind structure by comparing simultaneous singleDoppler and dualDoppler retrievals in Hurricane Danny (1999) using groundbased radar observations. However, to the authors' knowledge, no previous studies have conducted a detailed comparison of the Fourier decomposition of the wind field between the groundbased singleDoppler and airborne dualDoppler wind retrievals. We have already shown that wave number 0 tangential wind from the dualDoppler analyses can be accurately retrieved by the GVTDimproved algorithm in Sect. 4.1. With the confidence of the performance of the GVTD technique, the next step is to examine the lowwavenumber structure and applicability of a steadystate assumption of the dualDoppler analysis.
Figure 5a displays ${V}_{\mathrm{d}}D/{R}_{\mathrm{T}}$ on a ring at the RMW (18 km) derived from the dualDoppler and singleDoppler analyses. For clarity, we use the term “harmonics” to refer to Fourier decomposition in ${V}_{\mathrm{d}}D/{R}_{\mathrm{T}}$ around the ring and “wave numbers” to refer to the components of tangential wind which are constructed from combinations of the observed Doppler harmonics. ${V}_{\mathrm{d}}D/{R}_{\mathrm{T}}$ constructed from both analyses are dominated by a single harmonic pattern with a very similar magnitude. Nevertheless, a small difference in the ${V}_{\mathrm{d}}D/{R}_{\mathrm{T}}$ pattern between the two analyses is evident. The discrepancies of the pattern and magnitude suggest that higher harmonics of the ${V}_{\mathrm{d}}D/{R}_{\mathrm{T}}$ coefficients are different.
Figure 5b shows the residuals after the subtraction of the harmonics 0 and 1 components from the total ${V}_{\mathrm{d}}D/{R}_{\mathrm{T}}$, with the harmonic 2 components of ${V}_{\mathrm{d}}D/{R}_{\mathrm{T}}$ highlighted. The residuals include the harmonic 2 and higher components of ${V}_{\mathrm{d}}D/{R}_{\mathrm{T}}$. The residuals from the singleDoppler analysis (light red line) are dominated by a wave number 2 component, whereas the residuals from the dualDoppler analysis (light green line) are without a clear pattern. The peak value of the residual from the singleDoppler analysis is 8 m s^{−1}, but the peak value of the residual from the dualDoppler analysis is about 3 m s^{−1}. To quantitatively compare the two retrievals, Table 4 shows the retrieved magnitude of harmonics 0–3 around the RMW. The magnitude of harmonics 0, 1, and 3 are similar at the RMW for the two analyses with around a 1 m s^{−1} deviation. The similarity in harmonic 1 is primarily responsible for the similarity in the wave number 0 tangential wind (Fig. 4). Harmonic 2 shows a much larger deviation of 6 m s^{−1}, resulting in different wave number 1 tangential wind retrievals.
Figure 6 shows the harmonics of ${V}_{\mathrm{d}}D/{R}_{\mathrm{T}}$ as a function of radius for the resampled dualDoppler wind field (green dashed line) and singleDoppler observations (red dots). Note that the ordinate on each panel is different due to the varying magnitudes of each coefficient. The deviations of the A_{0}, A_{1}, and B_{1} coefficients between the two analyses within the eyewall region (15–25 km) are less than 2 m s^{−1}. The deviations in the A_{0}, A_{1}, and B_{1} coefficients are larger from the eye to the inner edge of the eyewall (10–15 km) and outside of the eyewall region (beyond 25 km), but they generally have a similar magnitude. On the other hand, larger discrepancies are apparent with A_{2} and B_{2}. A_{2} is similar between 10 and 20 km but has a large 7 m s^{−1} difference outside 20 km. B_{2} is different at nearly all radii, with up to 7 m s^{−1} differences. The different amplitudes of A_{2} and B_{2} indicate that the retrieved wave number 1 tangential wind phase and magnitude are inconsistent (Eq. 19). Discrepancies are also evident in the A_{3} and B_{3} coefficients, indicating inconsistencies in the retrieved wave number 2 tangential wind phase and magnitude as well.
We hypothesize that the discrepancies in retrieved wave number 1 and 2 tangential winds are due to the steadystate assumption in the dualDoppler wind synthesis into one snapshot. Two different types of sampling errors in effect arise from the steadystate assumption: the first is due to the time lag between fore and aft beams, and the second is due to the length of time used to composite the multiDoppler into a single snapshot. While both produce some errors, the latter is more consequential when considering that the temporal evolution of the phenomena is faster than the period of data collection, resulting in evolution over the flight pass. The “local” wind may be correct, but the overall structure is distorted by collapsing to a single time. For example, the propagation velocity of a wave number 2 vortex Rossby wave (VRW) is half of the symmetric tangential wind velocity (Lamb, 1932; Guinn and Schubert, 1993; Kuo et al., 1999). A propagating wave number 2 asymmetry could then alias onto other wave numbers, contributing to a discrepancy in wave number 1 tangential wind.
To test the hypothesis, the phases of maximum wave number 1 and wave number 2 tangential winds retrieved from the 5 min singleDoppler observations are examined in Fig. 7 for the temporal evolution during the first flight pass from 19:07 to 19:40 UTC. The amplitude and phase (Eqs. 18 and 19) of wave number 1 and 2 tangential winds are denoted in polar coordinates by the radius and azimuth, respectively. The wave number 1 tangential wind (Fig. 7a) generally remained unchanged throughout the first pass with a magnitude between 8 and 12 m s^{−1} and phase to the east to northeast (same as the wave number 1 reflectivity, not shown here). Environmental vertical wind shear derived from the Statistical Hurricane Intensity Prediction Scheme dataset (SHIPS) points to the northeast direction with a magnitude of 7 m s^{−1}, suggesting that the wave number 1 distribution is forced by the vertical wind shear to be consistently in the downshearright quadrant.
The wave number 2 tangential wind (Fig. 7b) propagated cyclonically during the flight pass (same as the wave number 2 reflectivity, not shown here) with a magnitude up to 7 m s^{−1}. The propagation of the wave number 2 tangential wind is estimated to be 285^{∘} from 19:07 to 19:40 UTC, which is 35 m s^{−1}, or 63 % of ${V}_{{\mathrm{T}}_{\mathrm{max}}}$. The propagation of the wave number 2 tangential wind is roughly consistent with linear VRW theory (Kuo et al., 1999; Cha et al., 2020).
According to the singleDoppler analysis, the wave number 2 VRW propagated approximately 43^{∘} every 5 min on average. The wave number 2 maximum amplitudes were located in the northeastern and southwestern quadrants as the P3 approached the inner eyewall from the southwest (Fig. 2), such that the TDR would see the maximum component in the southwestern eyewall. The VRW rotated cyclonically as the P3 crossed the 36 km diameter eye at 120 m s^{−1}, such that by the time the P3 exited the eyewall the radar would begin to see the wave number 2 wind minimum on the northeastern side. The southwestern maximum and northeastern minimum would then appear as a wave number 1 component under a steadystate assumption. The analysis supports our hypothesis that the propagation of wave number 2 tangential wind is aliased onto the steady wave number 1 component, resulting in a reduced amplitude and a phase shift in A_{2} and B_{2} in the dualDoppler analysis. As the axisymmetric radial wind is influenced by the harmonics 2 and 3 as well as the wave number 2 radial wind component (Eq. 20), we cannot fully validate the axisymmetric radial wind retrieval with the current dataset. Lee et al. (2006) shows that the Lamb solution of the wave number 2 radial wind has a comparable magnitude to the wave number 2 tangential wind but with a phase shift, so the wave number 0 radial wind retrieval is uncertain when VRWs are present. An evaluation of the accuracy of the axisymmetric radial wind retrieval is not included in this study.
The above analysis suggests that both the dualDoppler and singleDoppler analysis have strengths and weaknesses. The retrieved wave number 0 component of tangential wind has been shown to be comparable between an airborne dualDoppler wind synthesis and singleDoppler retrieval due to the fact that a steadystate assumption for one flight pass is usually valid to apply to the axisymmetric circulation. The axisymmetric kinematic structure of Hurricane Matthew revealed in the dualDoppler analysis is consistent with the singleDoppler retrievals across multiple aircraft passes (Cha, 2018). However, the steadystate assumption may not be applicable with more rapidly evolving features, such as convective bursts, rotation of mesovortices, and propagation of a wave number 2 or higher VRWs. The evolution during the aircraft pass could potentially impact the asymmetric reflectivity and wind structure retrieval from a dual or multiDoppler analysis due to temporal aliasing. On the other hand, dualDoppler analysis can retrieve highquality and detailed threedimensional structure, whereas a groundbased singleDoppler radar retrieval is often limited by spatial aliasing and the observing distance from the radar to TC. Nevertheless, the analyses presented herein demonstrate that complementary information can be retrieved from both single and dualDoppler retrievals.
4.3 An idealized experiment with a propagating wave number 2 asymmetry
To further test our hypothesis, we performed an idealized experiment to sample a rotating wave number 2 asymmetry by an aircraft with a realistic scanning strategy, similar to that first presented by Bell et al. (2007). The idealized framework allows for minimization of instrument and algorithm errors in order to isolate the sampling error due to the steadystate assumption. The simulated aircraft penetrates the TC eyewall at a 120 m s^{−1} airspeed on a south to north flight track, which is a typical straightline flight track during NOAA P3 operational reconnaissance. The radar has “perfect” FAST scanning with no instrument error and ideal geometry from the straight flight track, such that a steadystate wind field can be retrieved with minimal error (not shown). The TC wind field is a symmetric Rankine vortex with a 20 km RMW and V_{max} of 50 m s^{−1}, representative of a strong mature TC, with an added propagating wave number 2 linear Rankine edge wave with an epsilon of 3 km corresponding to an amplitude of 7.1 m s^{−1} as in Lee et al. (2006). The initial phase of the asymmetry is oriented from east to west, and it rotates around the RMW cyclonically at $\mathrm{1}/\mathrm{2}$ V_{max}. The simulated aircraft takes 22 min to cross the 160 km distance through the storm, during which time the wave number 2 asymmetry rotates to have a north–south orientation. Therefore, the timemean phase of the asymmetry is 45^{∘}, oriented from southwest to northeast over the course of the sampling period.
Figure 8a shows the derived tangential wind for a wave number 2 rotation during the flight pass. The overall magnitude of the wind field is close to the prescribed timemean structure but is distorted, suggesting temporal aliasing from the extended sampling period of the aircraft pass due to the propagation of the wave number 2 asymmetry. Figure 8b–d show the retrieved wave number 1, 2, and 3 tangential wind fields from the Fourier decomposition. Although only a wave number 2 component is prescribed in the experiment, the retrieved wave number 1 and 3 components are up to 1.5–2 m s^{−1}. The orientation of the flow is close to the “true” timemean structure shown in black contours, and the Fourier wave number 2 component has the generally correct phase. We note that the amplitude of the wave number 2 component correctly changes sign at the RMW due to the construction of a pure vorticity wave and is consistent with the prescribed flow and phase. The exact asymmetry retrieved depends to a large extent on the specifics of the flight track. However, we note that over the 22 min period, a groundbased singleDoppler retrieval would see the storm about four times with different phase orientations of the asymmetry (not shown). While idealized, this experiment provides further evidence that the steadystate assumption to synthesize the data into one snapshot with the presence of evolving features could result in temporal aliasing and result in potential misinterpretation of the asymmetric flow field.
Tropical cyclone (TC) wind retrieval techniques from Doppler radar observations are salient to identify the convective and kinematic structure and evolution of TCs. The current study is believed to be the first detailed comparison of the retrieved wind field between groundbased singleDoppler and airborne dualDoppler wind techniques. Hurricane Matthew (2016) was observed by the polarimetric Next Generation Weather Radar (NEXRAD) in Miami (KAMX) and National Oceanic and Atmospheric Administration WP3D (NOAA P3) airborne radar near the coast of Florida for several hours, providing a novel opportunity to evaluate and compare single and dualDoppler wind retrieval techniques for TC flows.
Jou et al. (2008) has shown that the GVTD technique improves the capability of the GBVTD and can retrieve a subset of the wind field from a single groundbased Doppler radar under the assumption of nearly axisymmetric rotational wind. Here, we present an improved technique that mitigates errors due to storm motion, which yields a more accurate estimation of mean wind as well as wave number 0 tangential and radial winds. The aliasing of unresolved wind components from the singleDoppler measurements remains, so another closure assumption for the GVTD technique is needed to improve the accuracy of the retrieval in the future. In addition to the wind improvements, the GVTDsimplex centers more closely follow the dynamic centers compared with GBVTDsimplex centers due to the more accurate estimation of axisymmetric tangential wind. However, the GVTDsimplex centers are still variable and not fully consistent between retrievals from different radars (such as KAMX vs. KMLB, and KMLB vs. KJAX), suggesting further improvement is needed in the technique. Nevertheless, the improvement of TC wind estimation from singleDoppler radar presented herein can be useful in deducing storm structure for research purposes and assimilating realtime TC intensity and structure for forecasts (Zhao et al., 2012).
The strengths and weaknesses of different platform observing capabilities and retrieval techniques for TC wind fields are discussed in this study. The full threedimensional wind field from airborne radar fore–aft pseudodualDoppler scanning can be retrieved by a splinebased technique called SAMURAI, but it is shown to have errors due to temporal aliasing from the nonsimultaneous Doppler measurements when the temporal evolution of the phenomena is faster than the period of data collection. SingleDoppler wind retrieval by the GVTD technique can be obtained at a 5 min interval, but spatial aliasing remains due to the limited measurements from only the radial direction of the radar beam. A comparison between the two techniques shows that the retrieved axisymmetric (wave number 0) component of tangential winds is generally comparable between the two techniques, and the improved GVTD technique improves the accuracy of the retrieval. Fourier decomposition of asymmetric kinematic and convective structure shows more discrepancies between the two techniques due to spatial and temporal aliasing in the retrievals. The propagation of wave number 2 tangential wind is found to be aliased onto the steady wave number 1 component retrieved by the dualDoppler analysis, causing a reduced amplitude and phase shift of wave number 1 tangential wind. The steadystate assumption for lowwavenumber structure retrievals from pseudodualDoppler analysis may not be applicable with the presence of rapidly evolving features, and the temporal evolution during a flight pass should be considered. On the other hand, the steadystate assumption is less severe for singleDoppler wind retrievals, and rapidly evolving phenomena, such as vortex Rossby waves, have been documented by several studies using singleDoppler observations (Corbosiero et al., 2006; Cha et al., 2020). The analyses presented herein demonstrate that complementary information can be retrieved from both single and dualDoppler retrievals. Future work on the timedependent analysis of asymmetric structure in dualDoppler retrieval and asymmetric flow closure assumptions in singleDoppler retrieval are required to reconcile differences in retrieved TC structure.
The underlying Julia code used for the analysis in this study is available from the corresponding author upon request. The Vortex Objective Radar Tracking And Circulation (VORTRAC) and SAMURAI software are part of the Lidar Radar Open Software Environment (LROSE) and are available in Bell (2019, https://doi.org/10.5281/zenodo.2532758). The data used in this paper are available through Zenodo: https://doi.org/10.5281/zenodo.4427194 (Cha, 2020).
The study was designed by both authors. TYC carried out the analysis and made the figures under the supervision of MMB. TYC wrote the paper with contributions from MMB.
The authors declare that they have no conflict of interest.
This research was supported by the National Science Foundation (grant no. OAC1661663). We would like to thank the NOAA Aircraft Operations Center and the Hurricane Research Division of the Atlantic Oceanographic and Meteorological Laboratory for collecting the airborne tail Doppler radar data used for this study, and we are grateful to the National Weather Service for the groundbased radar data. We thank Naufal Razin and the three anonymous reviewers for their constructive and helpful comments that improved the quality of the paper.
This research has been supported by the National Science Foundation, Office of Advanced Cyberinfrastructure (grant no. OAC1661663).
This paper was edited by S. Joseph Munchak and reviewed by three anonymous referees.
Bell, M. M.: nsflrose/lroseblaze: lroseblaze20190105 (Version lroseblaze20190105), Zenodo [data set], https://doi.org/10.5281/zenodo.2532758, 2019. a, b
Bell, M. M., Lee, W.C., Houze, R. A., Smull, B. F., and Marks, F. D.: Circular Flight Tracks for Deducing Tropical Cyclone Innercore Structure Using Airborne Doppler Radar, 33rd Conference on Radar Meteorology, Cairns, Australia, [poster], 6–10 August 2007. a
Bell, M. M., Montgomery, M. T., and Emanuel, K. A.: Air–sea enthalpy and momentum exchange at major hurricane wind speeds observed during CBLAST, J. Atmos. Sci., 69, 3197–3222, https://doi.org/10.1175/JASD110276.1, 2012. a, b, c
Bell, M. M., Lee, W.C., Wolff, C. A., and Cai, H.: A Solobased automated quality control algorithm for airborne tail Doppler radar data, J. Appl. Meteorol. Clim., 52, 2509–2528, https://doi.org/10.1175/JAMCD120283.1, 2013. a
Cha, T.Y.: Eyewall Replacement Cycle of Hurricane Matthew (2016) observed by Doppler radars, MS thesis, Colorado State University, 2018. a
Cha, T.Y.: Data for amt2020240, Zenodo, https://doi.org/10.5281/zenodo.4427194, 2020. a
Cha, T.Y., Bell, M. M., Lee, W.C., and DesRosiers, A. J.: Polygonal Eyewall Asymmetries During the Rapid Intensification of Hurricane Michael (2018), Geophys. Res. Lett., 47, e2020GL087919, https://doi.org/10.1029/2020GL087919, 2020. a, b
Chan, J. C.: An observational study of the physical processes responsible for tropical cyclone motion, J. Atmos. Sci., 41, 1036–1048, https://doi.org/10.1175/15200469(1984)041<1036:AOSOTP>2.0.CO;2, 1984. a
Chan, J. C. and Gray, W. M.: Tropical cyclone movement and surrounding flow relationships, Mon. Weather Rev., 110, 1354–1374, https://doi.org/10.1175/15200493(1982)110<1354:TCMASF>2.0.CO;2, 1982. a
Chen, X., Zhao, K., Lee, W.C., JongDao Jou, B., Xue, M., and Harasti, P. R.: The Improvement to the Environmental Wind and Tropical Cyclone Circulation Retrievals with the Modified GBVTD (MGBVTD) Technique, J. Appl. Meteorol. Clim., 52, 2493–2508, https://doi.org/10.1175/JAMCD13031.1, 2013. a
Corbosiero, K. L., Molinari, J., Aiyyer, A. R., and Black, M. L.: The structure and evolution of Hurricane Elena (1985). Part II: Convective asymmetries and evidence for vortex Rossby waves, Mon. Weather Rev., 134, 3073–3091, https://doi.org/10.1175/MWR3250.1, 2006. a
Crum, T. D., Saffle, R. E., and Wilson, J. W.: An Update on the NEXRAD Program and Future WSR88D Support to Operations, Weather Forecast., 13, 253–262, https://doi.org/10.1175/15200434(1998)013<0253:AUOTNP>2.0.CO;2, 1998. a
Gamache, J. F., Marks, F. D., and Roux, F.: Comparison of Three Airborne Doppler Sampling Techniques with Airborne In Situ Wind Observations in Hurricane Gustav (1990), J. Atmos. Ocean. Tech., 12, 171–181, https://doi.org/10.1175/15200426(1995)012<0171:COTADS>2.0.CO;2, 1995. a
Griffin, J. S., Burpee, R. W., Marks, F. D., and Franklin, J. L.: Realtime airborne analysis of aircraft data supporting operational hurricane forecasting, Weather Forecast., 7, 480–490, https://doi.org/10.1175/15200434(1992)007<0480:RTAAOA>2.0.CO;2, 1992. a
Guinn, T. A. and Schubert, W. H.: Hurricane Spiral Bands, J. Atmos. Sci., 50, 3380–3403, https://doi.org/10.1175/15200469(1993)050<3380:HSB>2.0.CO;2, 1993. a
Harasti, P. R., McAdie, C. J., Dodge, P. P., Lee, W.C., Tuttle, J., Murillo, S. T., and Marks, F. D.: RealTime Implementation of SingleDoppler Radar Analysis Methods for Tropical Cyclones: Algorithm Improvements and Use with WSR88D Display Data, Weather Forecast., 19, 219–239, https://doi.org/10.1175/15200434(2004)019<0219:RIOSRA>2.0.CO;2, 2004. a
Hildebrand, P. H. and Mueller, C. K.: Evaluation of Meteorological Airborne Doppler Radar. Part I: DualDoppler Analyses of Air Motions, J. Atmos. Ocean. Tech., 2, 362–380, https://doi.org/10.1175/15200426(1985)002<0362:EOMADR>2.0.CO;2, 1985. a, b
Hildebrand, P. H., Walther, C., and Frush, C.: The NCAR Electra Doppler radar, in: Part I: Evaluation of scientific needs, Preprints, 23rd Conf. on Radar Meteorology, Snowmass, CO, 22–26 September 1986, Amer. Meteor. Soc, 147–150, 1986. a
Hildebrand, P. H., Walther, C. A., Frush, C. L., Testud, J., and Baudin, F.: The ELDORA/ASTRAIA airborne Doppler weather radar: goals, design, and first field tests, Proceedings of the IEEE, 82, 1873–1890, https://doi.org/10.1109/5.338076, 1994. a
Houze Jr., R. A.: Clouds in Tropical Cyclones, Mon. Weather Rev., 138, 293–344, https://doi.org/10.1175/2009MWR2989.1, 2010. a
Jorgensen, D., Matejka, T., and DuGranrut, J.: Multibeam techniques for deriving wind fields from airborne Doppler radars, Meteorol. Atmos. Phys., 59, 83–104, https://doi.org/10.1007/BF01032002, 1996. a
Jorgensen, D. P., Hildebrand, P. H., and Frush, C. L.: Feasibility Test of an Airborne PulseDoppler Meteorological Radar, J. Appl. Meteorol. Clim., 22, 744–757, https://doi.org/10.1175/15200450(1983)022<0744:FTOAAP>2.0.CO;2, 1983. a, b
Jou, B. J.D., Deng, S.M., and Chang, B.: Determination of typhoon center and radius of maximum wind by using Doppler radar, Atmos. Sci., 24, 1–24, 1996. a, b
Jou, B. J.D., Lee, W.C., Liu, S.P., and Kao, Y.C.: Generalized VTD retrieval of atmospheric vortex kinematic structure. Part I: Formulation and error analysis, Mon. Weather Rev., 136, 995–1012, https://doi.org/10.1175/2007MWR2116.1, 2008. a, b, c, d, e, f, g, h, i, j, k, l, m
Koch, S. E., desJardins, M., and Kocin, P. J.: An Interactive Barnes Objective Map Analysis Scheme for Use with Satellite and Conventional Data, J. Appl. Meteorol. Clim., 22, 1487–1503, https://doi.org/10.1175/15200450(1983)022<1487:AIBOMA>2.0.CO;2, 1983. a
Kuo, H.C., Williams, R., and Chen, J.H.: A possible mechanism for the eye rotation of Typhoon Herb, J. Atmos. Sci., 56, 1659–1673, https://doi.org/10.1175/15200469(1999)056<1659:APMFTE>2.0.CO;2, 1999. a, b
Lamb, H.: Hydrodynamics, 6th edn., Dover, 732 pp., 1932. a
Lee, W.C. and Bell, M. M.: Rapid intensification, eyewall contraction, and breakdown of Hurricane Charley (2004) near landfall, Geophys. Res. Lett., 34, L02802, https://doi.org/10.1029/2006GL027889, 2007. a
Lee, W.C. and Marks, F. D.: Tropical cyclone kinematic structure retrieved from singleDoppler radar observations. Part II: The GBVTDsimplex center finding algorithm, Mon. Weather Rev., 128, 1925–1936, https://doi.org/10.1175/15200493(2000)128<1925:TCKSRF>2.0.CO;2, 2000. a, b, c, d
Lee, W.C., Marks, F. D., and Carbone, R. E.: Velocity track display – A technique to extract realtime tropical cyclone circulations using a single airborne Doppler radar, J. Atmos. Ocean. Tech., 11, 337–356, https://doi.org/10.1175/15200426(1994)011<0337:VTDTTE>2.0.CO;2, 1994. a, b, c
Lee, W.C., Jou, B. J.D., Chang, P.L., and Deng, S.M.: Tropical cyclone kinematic structure retrieved from singleDoppler radar observations. Part I: Interpretation of Doppler velocity patterns and the GBVTD technique, Mon. Weather Rev., 127, 2419–2439, https://doi.org/10.1175/15200493(1999)127<2419:TCKSRF>2.0.CO;2, 1999. a, b, c, d, e, f, g, h
Lee, W.C., Jou, B. J.D., Chang, P.L., and Marks, F. D.: Tropical cyclone kinematic structure retrieved from singleDoppler radar observations. Part III: Evolution and structures of Typhoon Alex (1987), Mon. Weather Rev., 128, 3982–4001, https://doi.org/10.1175/15200493(2000)129<3982:TCKSRF>2.0.CO;2, 2000. a, b, c, d
Lee, W.C., Marks, F. D., and Walther, C.: Airborne Doppler radar data analysis workshop, B. Am. Meteorol. Soc., 84, 1063–1075, 2003. a
Lee, W.C., Harasti, P. R., Bell, M. M., Jou, B. J.D., and Chang, M.H.: Doppler velocity signatures of idealized elliptical vortices, TAO: Terr. Atmos. Ocean. Sci., 17, 429–446, 2006. a, b
Liou, Y.C., Wang, T.C. C., Lee, W.C., and Chang, Y.J.: The retrieval of asymmetric tropical cyclone structures using Doppler radar simulations and observations with the extended GBVTD technique, Mon. Weather Rev., 134, 1140–1160, https://doi.org/10.1175/MWR3107.1, 2006. a
Lorsolo, S. and Aksoy, A.: Wavenumber analysis of azimuthally distributed data: Assessing maximum allowable gap size, Mon. Weather Rev., 140, 1945–1956, https://doi.org/10.1175/MWRD1100219.1, 2012. a, b
Marks, F. D. and Houze, R. A.: Airborne Doppler Radar Observations in Hurricane Debby, B. Am. Meteorol. Soc., 65, 569–582, https://doi.org/10.1175/15200477(1984)065<0569:ADROIH>2.0.CO;2, 1984. a
Marks, F. D. and Houze, R. A.: Inner Core Structure of Hurricane Alicia from Airborne Doppler Radar Observations, J. Atmos. Sci., 44, 1296–1317, https://doi.org/10.1175/15200469(1987)044<1296:ICSOHA>2.0.CO;2, 1987. a
Marks, F. D., Houze, R. A., and Gamache, J. F.: Dualaircraft investigation of the inner core of Hurricane Norbert. Part I: Kinematic structure, J. Atmos. Sci., 49, 919–942, https://doi.org/10.1175/15200469(1992)049<0919:DAIOTI>2.0.CO;2, 1992. a, b, c
Matejka, T. and Srivastava, R. C.: An Improved Version of the Extended VelocityAzimuth Display Analysis of SingleDoppler Radar Data, J. Atmos. Ocean. Tech., 8, 453–466, https://doi.org/10.1175/15200426(1991)008<0453:AIVOTE>2.0.CO;2, 1991. a
Murillo, S. T., Lee, W.C., Bell, M. M., Barnes, G. M., Marks, F. D., and Dodge, P. P.: Intercomparison of GroundBased Velocity Track Display (GBVTD)Retrieved Circulation Centers and Structures of Hurricane Danny (1997) from Two Coastal WSR88Ds, Mon. Weather Rev., 139, 153–174, https://doi.org/10.1175/2010MWR3036.1, 2011. a
Ooyama, K. V.: ScaleControlled Objective Analysis, Mon. Weather Rev., 115, 2479–2506, https://doi.org/10.1175/15200493(1987)115<2479:SCOA>2.0.CO;2, 1987. a
Ooyama, K. V.: The CubicSpline Transform Method: Basic Definitions and Tests in a 1D Single Domain, Mon. Weather Rev., 130, 2392 – 2415, https://doi.org/10.1175/15200493(2002)130<2392:TCSTMB>2.0.CO;2, 2002. a
Ray, P. S. and Stephenson, M.: Assessment of the Geometric and Temporal Errors Associated with Airborne Doppler Radar Measurements of a Convective Storm, J. Atmos. Ocean. Tech., 7, 206–217, https://doi.org/10.1175/15200426(1990)007<0206:AOTGAT>2.0.CO;2, 1990. a
Reasor, P. D., Montgomery, M. T., Marks, F. D., and Gamache, J. F.: Lowwavenumber structure and evolution of the hurricane inner core observed by airborne dualDoppler radar, Mon. Weather Rev., 128, 1653–1680, https://doi.org/10.1175/15200493(2000)128<1653:LWSAEO>2.0.CO;2, 2000. a
Shimada, U., Sawada, M., and Yamada, H.: Doppler Radar Analysis of the Rapid Intensification of Typhoon Goni (2015) after Eyewall Replacement, J. Atmos. Sci., 75, 143–162, https://doi.org/10.1175/JASD170042.1, 2018. a
Willoughby, H. E. and Chelmow, M. B.: Objective determination of hurricane tracks from aircraft observations, Mon. Weather Rev., 110, 1298–1305, https://doi.org/10.1175/15200493(1982)110<1298:ODOHTF>2.0.CO;2, 1982. a, b, c
Wood, V. T. and Brown, R. A.: Effects of radar proximity on singleDoppler velocity signatures of axisymmetric rotation and divergence, Mon. Weather Rev., 120, 2798–2807, https://doi.org/10.1175/15200493(1992)120<2798:EORPOS>2.0.CO;2, 1992. a
Zhao, K., Lee, W.C., and Jou, B. J.D.: Single Doppler radar observation of the concentric eyewall in Typhoon Saomai, 2006, near landfall, Geophys. Res. Lett., 35, L07807, https://doi.org/10.1029/2007GL032773, 2008. a
Zhao, K., Xue, M., and Lee, W.C.: Assimilation of GBVTDretrieved winds from singleDoppler radar for shortterm forecasting of super typhoon Saomai (0608) at landfall, Q. J. Roy. Meteor. Soc., 138, 1055–1071, https://doi.org/10.1002/qj.975, 2012. a, b
 Abstract
 Introduction
 Datasets and methodology
 The GVTD technique improvement
 Wind retrieval comparison between singleDoppler and airborne dualDoppler analyses
 Conclusions
 Code and data availability
 Author contributions
 Competing interests
 Acknowledgements
 Financial support
 Review statement
 References
 Abstract
 Introduction
 Datasets and methodology
 The GVTD technique improvement
 Wind retrieval comparison between singleDoppler and airborne dualDoppler analyses
 Conclusions
 Code and data availability
 Author contributions
 Competing interests
 Acknowledgements
 Financial support
 Review statement
 References