Extracting multiple-scale observational information is critical for accurately reconstructing the structure of mesoscale circulation systems such as typhoon. The Space and Time Mesoscale Analysis System (STMAS) with multigrid data assimilation developed in Earth System Research Laboratory (ESRL) in National Oceanic and Atmospheric Administration (NOAA) has addressed this issue. Previous studies have shown the capability of STMAS to retrieve multiscale information in 2-dimensional Doppler radar radial velocity observations. This study explores the application of 3-dimensional (3D) Doppler radar radial velocities with STMAS for reconstructing a 3D typhoon structure. As for the first step, here, we use an idealized simulation framework. A two-scale simulated “typhoon” field is constructed and referred to as “truth,” from which randomly distributed conventional wind data and 3D Doppler radar radial wind data are generated. These data are used to reconstruct the synthetic 3D “typhoon” structure by the STMAS and the traditional 3D variational (3D-Var) analysis. The degree by which the “truth” 3D typhoon structure is recovered is an assessment of the impact of the data type or analysis scheme being evaluated. We also examine the effects of weak constraint and strong constraint on STMAS analyses. Results show that while the STMAS is superior to the traditional 3D-Var for reconstructing the 3D typhoon structure, the strong constraint STMAS can produce better analyses on both horizontal and vertical velocities.
National Basic Research Program2013CB430304National Natural Science Foundation of China4137601341376015413060064154104141506039National High-Tech R&D Program2013AA09A505National Programme on Global Change and Air-Sea InteractionGASI-01-01-12GASI-IPOVAI-041. Introduction
Doppler radar has long been a valuable observational tool in meteorology. Three-dimensional (3D) Doppler radar radial velocity data can provide an opportunity to estimate both horizontal and vertical velocities. Therefore, in recent years, Doppler radar data assimilation for short-term numerical weather forecasting or called nowcasting has become a focal point of research [1–4]. Lots of techniques have been developed to retrieve wind field from Doppler radar radial velocity observations [2–24].
In Doppler radar radial velocity data assimilation used in the above literatures, in a three-dimensional variational (3D-Var) framework, a background error covariance matrix is always needed to determine the spatial spreading of observational information. It is well known that an analysis field at different locations may have different correlation scales [25], which are difficult to be estimated. Unfortunately, the traditional 3D-Var always employs an empirical and static background error covariance matrix and therefore usually can only correct single-scale wavelength error. However, the errors in short wavelength scales cannot be sufficiently corrected until the long waves are corrected [25, 26].
To minimize the errors of long and short waves in turn, a sequential 3D-Var approach has been proposed by Xie et al. [25, 26], implemented by either a recursive filter [27] or a multigrid technique [28] at Global Systems Division (GSD) of Earth System Research Laboratory (ESRL) in National Oceanic and Atmospheric Administration (NOAA) for a Federal Aviation Agency (FAA) project joined by the research team from the Lincoln Laboratory in Massachusetts Institute of Technology (MIT). Since this system also uses the temporal observation information, it is called a Space and Time Mesoscale Analysis System (STMAS, thereafter; see Xie et al. [26]). The STMAS has been applied to assimilating 2-dimensional (2D) Doppler radar radial velocity data to improve the wind field analyses [29].
Here, we study the analysis of 3D Doppler radar radial velocities using the STMAS to reconstruct the 3D wind structure. As for the first step, this study is performed in a twin experiment framework. In the next section, we first briefly review the theory of the multigrid 3D-Var data assimilation scheme in the STMAS. Some important aspects of the STMAS techniques such as smoothing, constraint, and Doppler radar radial wind operators used in the cost function of the STMAS multigrid 3D-Var are described. Section 3 first introduces the twin experiment framework for 3D Doppler radar radial wind data assimilation with the STMAS and then gives the evaluation by comparing it to the traditional 3D-Var. Section 3 also examines the performance of the STMAS in weak and strong constraints for 3D Doppler radar radial velocity analysis. Conclusions and discussions are given in Section 4.
2. Smoothing, Constraint, and Radar Radial Wind Operators in STMAS
In this study, the STMAS implemented by the multigrid 3D-Var is applied to the analysis of 3D Doppler radar radial velocities. This method can extract long and short wavelength information in turn efficiently from observations and provide objective and accurate analysis. The basic idea of this multigrid implementation can be referred to Li et al. [28–30].
To assimilate 3D Doppler radar radial velocities, with the control variables being Xn=UnT,VnTT, where U and V represent zonal and meridional components of wind vector, the cost functional for the nth level grid is(1)Jn=Jbn+Jsn+Jcn+Jrnn=1,2,3,…,N,Jbn=12UnTUn+12VnTVn,Jsn=UnTSUUn+VnTSVVn,where the subscript b denotes the background term and s the smooth term, c the conventional observation data term, and r the radar radial wind observation data term. The smooth matrixes SU and SV in the smooth term are derived from the Laplacian of control variables U and V, respectively, at grid points. Let W represent the vertical component of wind vector. The details of the conventional observation data term and the 3D Doppler radar radial wind observation data term are as follows:(2)Jcn=12ΗnUn-YUTOU-1ΗnUn-YU+12ΗnVn-YVTOV-1ΗnVn-YV,(3)Jrn=12ΘsinΗnUn+ΘcosΗnVn+ΦΗnWn-YRTORn-1ΘsinΗnUn+ΘcosΗnVn+ΦΗnWn-YR,(4)Θsin=diagsinθ1cosφ1,sinθ2cosφ2,…,sinθMcosφM,(5)Θcos=diagcosθ1cosφ1,cosθ2cosφ2,…,cosθMcosφM,(6)Φ=diagsinφ1,sinφ2,…,sinφM,where M is the amount of radar radial wind observations, θ is the azimuth angle of the radar beam relative to north with positive clockwise, and φ is elevation angle of the radar beam. Of course, since radar scans at nonzero elevation angles, the fall speed of precipitation particles should be taken into account, and the algorithm of Sun and Crook [6] can be used to calculate terminal velocity. But for this study, we just neglect this terminal velocity, which does not lose its generality. The matrix ORn is an error covariance matrix for radar radial wind observation; its superscript -1 stands for the reverse matrix, and its subscript R represents the radar radial wind observation.
During the procedure of sequential multiscale analyses, the operators O, Θsin, Θcos, and Φ remain the same when the full observation dataset is used through all multigrid levels; therefore, the superscript n is omitted from these operators.
To make a strong constraint on these three components of wind vector, incompressible continuity equation ∂U/∂x+∂V/∂y+∂W/∂z=0 is employed and discretized to calculate vertical velocity from the other two horizontal components. The discretized incompressible continuity equation is as follows:(7)Wi,j,kn=Wi,j,k-1n+zkn-zk-1n2Ui+1,j,k-1n-Ui-1,j,k-1nxi+1n-xi-1n+Vi,j+1,k-1n-Vi,j-1,k-1nyj+1n-yj-1n+Ui+1,j,kn-Ui-1,j,knxi+1n-xi-1n+Vi,j+1,kn-Vi,j-1,knyj+1n-yj-1n.The adjoint codes are recursively developed for W represented by U and V.
The study domain covers a 500km×500km square region with 10 km thickness. The Doppler radar locates at the center (250 km, 250 km) of the study domain. A simulated typhoon field can be constructed by using the following function which consists of two subsection functions:(8)FL,D,r=1,0≤r≤D,exp-r-D4L4,r>D.This formula allows that the 2-order derivatives of this function exist. Let r=x-xc2+y-yc2 and Gx,y=A1FL1,D1,rx,y+A2FL2,D2,rx,y, where A1=-5 km2s−1, L1=70 km, D1=15 km, A2=-1 km2s−1, L2=20 km, and D2=2.5 km; then a stream function of a two-scale typhoon field can be constructed. The amplitude of one large-scale (L1/2=35 km) is scaled by A1 and the amplitude of one small-scale (L2/2=10 km) is scaled by A2. The locations of the maximum wind horizontal velocity of these two scales are 0.754L1+D1≈80 km for large-scale and 0.754L2+D2≈21 km for small-scale, respectively, from the typhoon center.
The typhoon center xc,yc is set at 300 km and 150 km. Then, the stream function and velocity potential function can be constructed as follows:(9)ψ=Gx,y1-exp-zD3,χ=-εGx,ycoszL3π,where ε=0.5, L3=10 km, and D3=1 km. The horizontal components can be expressed in terms of ψ and χ:(10)U=∂χ∂x-∂ψ∂y,V=∂χ∂y+∂ψ∂x.In this study, incompressibility is assumed. And the true vertical velocity field can be obtained by integrating the continuity equation ∂U/∂x+∂V/∂y+∂W/∂z=0. That is,(11)W=-∫0z∂U∂x+∂V∂ydz=-∫0z∂2χ∂x2+∂2χ∂y2dz=εL3π∂2G∂x2+∂2G∂y2sinzL3π,where the bottom boundary condition is Wx,y,0=0. The wind speed field in this simulated typhoon field contains two different scale information. The first one is about 35 km and the other is about 10 km. The radial wind, U component, and V component pattern of middle level of this simulated typhoon field and a section wind field across the center of this typhoon are shown in Figure 1. This typhoon pattern is located at the southeast part of the study domain, so only the southeast square part is shown for the detailed structure. This simulated typhoon wind field is referred to as the “truth” typhoon field in this twin experiment.
The synthetic “typhoon” field (serving as the “truth”) at middle vertical level (5 km). (a) Radial velocity, (b) U component, (c) V component, (d) wind vector, (e) W component, and (f) W distribution at an x-z section across the center of the synthetic “typhoon” with black arrows representing wind vector at this section (unit: m/s).
3.2. Observations
Then, Doppler radar radial velocity data are generated from the “truth” typhoon field with one-degree azimuth angle increment and 2500 m gate spacing and 2-degree elevation angle increment from 1 degree to 20 degrees by interpolating the “truth” velocity field to the radial velocity observations’ points and using the equation YR=ΘsinΗU+ΘcosΗV+ΦΗW, and the number of radial velocity data is 114480. Spatially coarse (1000), moderate (10000), and dense (100000) random distributed conventional observations are also generated, respectively, from the “truth” typhoon field by interpolating the “truth” velocity field to the conventional observations’ points.
3.3. Twin Experiments Setup
In the following, the above observational data are used to retrieve the simulated “typhoon” structure by the STMAS analysis method with weak or strong constraint and the traditional 3-dimensional variance (3D-Var) analysis with different correlation scales, respectively, and by comparing the analyzed results with the “truth,” performances of different analysis methods are discussed.
The error variances of radial velocity observations and conventional observations can be determined by the measurement error of instruments. But, here for simplicity, the same error variance is set for each kind of data. However, because the amount of radial velocity data is much larger than that of conventional data, a scaling scheme is used to balance the weights of these two types of observations. Thus, the conventional observation can have the same weight as that of radial velocity observation, which may comprise these two types of observations to get to a reasonable wind analysis.
The limited memory BFGS (Broyden-Fletcher-Goldfarb-Shanno) method [31] to solve the bound constrained optimization problem [32] is used as the minimization method in this study.
3.4. STMAS Results
Three level grids are employed ranging from about 31.25 km × 31.25 km × 2.5 km (n=1, and the number of grids is 17 × 17 × 5) to 7.8125 km × 7.8125 km × 0.625 km (n=3, and the number of grids is 65 × 65 × 17) with grid ratio being 0.5. The background is set to be zero for simplicity. The conventional data will be added gradually to investigate the impact of radar radial wind observations.
The STMAS analyses with the above 3D Doppler radar radial velocity data or conventional data which vary from coarse to dense are shown in Figures 2–4 and Table 1 gives the corresponding root mean square errors (RMSEs). Because the Doppler radar radial velocity data really provide some useful information of radial wind, the STMAS results by only using this type of data can make a good analysis on radial wind (see Figure 2(a)) and the two-scale information in radial wind has been captured by the STMAS. There is no information of tangential wind, and no statistical or empirical correlation information between U and V is used in this study. Therefore, only using radial velocity data, the STMAS analysis only can match the radial wind, and it cannot make a good analysis on U and V. This leads to large U and V RMSEs (15.7 m/s and 14.5 m/s, resp., see Table 1). With the 1000 coarse distributed conventional data only, the STMAS analysis only shows the large pattern of the typhoon wind (Figure 2(b)). With the combination of conventional data and radar radial wind data, the STMAS analysis significantly improves the typhoon structure (Figure 2(c)). The RMSE with the combined datasets is much smaller than the other one with an individual dataset (Table 1). Apparently, the radial velocity data provide additional information on radial direction to help the STMAS capture the detailed typhoon structure. If the conventional data change from the coarsest (totally 1000) to the densest (totally 100000), the detailed information on the typhoon is enhanced gradually (see Figures 2(b) and 2(c) corresponding to 1,000 conventional data and 1,000 conventional data combined with radar data, resp.; Figures 3(a) and 3(b) corresponding to 10,000 conventional data and 10,000 conventional data combined with radar data, resp.; Figures 4(a) and 4(b) corresponding to 100,000 conventional data and 100,000 conventional data combined with radar data, resp.). With the combination of radar radial velocity data and dense distributed conventional data, the STMAS can make almost perfect analysis for the typhoon structure (Figure 4(b)) and the RMSE is very small. And because dense distributed conventional data provide enough information, the RMSE improvements from the MG_DENSE to MG_RADAR_DENSE experiments are really negligible for the three wind components.
Root mean square errors (RMSE) of U, V, and W component in different experiments. MG_RADAR represents the STMAS results by only assimilating radial wind data. MG_COARSE represents the STMAS results by assimilating only coarse conventional data. MG_RADAR_COARSE represents the STMAS results by assimilating both radial wind data and coarse conventional data. MG_MODERATE and MG_RADAR_MODERATE are similar to MG_COARSE and MG_RADAR_COARSE but with the conventional data substituted by moderate conventional data. MG_DENSE and MG_RADAR_DENSE are similar to MG_COARSE and MG_RADAR_COARSE but with the conventional data substituted by dense conventional data. MG_RADAR_DENSE_WEAK is similar to MG_RADAR_DENSE but with weak constraint. T50_RADAR_DENSE represents the results of traditional 3D-Var with correlation scales being Lx=50 km, Ly=50 km, and Lz=5 km by assimilating both radial wind data and dense conventional data; T25_RADAR_DENSE is similar to T50_RADAR_DENSE but for Lx=25 km, Ly=25 km, and Lz=2.5 km; T12_RADAR_DENSE is similar to T50_RADAR_DENSE but for Lx=12.5 km, Ly=12.5 km, and Lz=1.25 km; T06_RADAR_DENSE is similar to T50_RADAR_DENSE but for Lx=6.25 km, Ly=6.25 km, and Lz=0.625 km.
Experiment
RMSE of U (m/s)
RMSE of V (m/s)
RMSE of W (m/s)
MG_RADAR
15.7032
14.5197
2.3663
MG_COARSE
9.2592
9.2195
1.4610
MG_RADAR_COARSE
7.9500
7.6029
1.3757
MG_MODERATE
3.2196
3.1966
0.8228
MG_RADAR_MODERATE
3.0649
2.6148
0.7930
MG_DENSE
1.0707
1.0849
0.5494
MG_RADAR_DENSE
1.0514
1.0135
0.5373
MG_RADAR_DENSE_WEAK
1.0747
1.0254
2.0840
T50_RADAR_DENSE
6.8944
7.1118
1.3832
T25_RADAR_DENSE
2.8425
2.7135
0.8172
T12_RADAR_DENSE
3.7725
3.6773
0.7787
T06_RADAR_DENSE
10.2998
10.2631
1.3669
STMAS analyses at middle vertical level. From left to right, these columns represent radial wind, U component, V component, and W component, respectively. (a) is the results by only assimilating radial wind data. (b) is the results by assimilating only coarse conventional data. (c) is the results by assimilating both radial wind data and coarse conventional data (unit: m/s).
(a) and (b) are similar to Figures 2(b) and 2(c) but with the conventional data substituted by moderate conventional data (unit: m/s).
(a) and (b) are similar to Figures 2(b) and 2(c) but with the conventional data substituted by dense conventional data. (c) is similar to (b) but with weak constraint (unit: m/s).
3.5. Strong Constraint versus Weak Constraint
To compare the performance of the STMAS with a strong constraint or weak constraint, the continuity equation is added as penalty term to make a weak constraint case (otherwise, the STMAS analysis is a strong constraint of the three components of wind vector). Then, the control variables become Xn=UnT,VnT,WnTT:(12)Jn=Jbn+Jsn+Jcn+Jrn+Jconstraintn,n=1,2,3,…,N.And background term is changed to (13)Jbn=12UnTUn+12VnTVn+12WnTWn.The smooth term is (14)Jsn=UnTSUUn+VnTSVVn+WnTSWWn.Here, SW is smooth matrix for W component of wind vector. The weak constraint term is (15)Jconstraintn=α∫∂Un∂x+∂Vn∂y+∂Wn∂z2dxdydz,where α is the penalty coefficient for weak constraint. The other terms, Jcn and Jrn, are the same as (2) and (3), respectively.
Assimilating the whole simulated Doppler radar radial velocity data and dense conventional data, the STMAS analyses with continuity equation as strong constraint or weak constraint are shown in Figures 4(b) and 4(c), respectively. And the RMSE is shown in Table 1. The coefficient of penalty term for weak constraint is set to 100.
As shown in Figure 4(b), the strong constraint forces these three components of wind vector to exactly satisfy the continuity equation. Then, while the horizontal components can be analyzed very well, the vertical component is reconstructed very well too. The weak constraint (Figure 4(c)) can only constrain the horizontal and vertical components of wind vector to satisfy the continuity equation to some degree. Whereas the order of vertical velocities is much smaller than that of horizontal ones and no direct observation of vertical velocities is available, the minimization of cost function primarily focuses on the horizontal parts. This causes little improvement on the vertical velocities even if the horizontal components can be well analyzed. Therefore, the vertical component RMSE score of MG_RADAR_DENSE_WEAK experiment is worse than that of MG_RADAR_DENSE experiment. When a large penalty coefficient α is used, although vertical velocity analysis can be improved a little bit suffering a lower accuracy of horizontal velocity, the accuracy of the analysis vertical velocities is still not high.
3.6. STMAS versus Traditional 3D-Var
For a traditional 3D-Var analysis, the cost function takes the form(16)J=Jb+Jc+Jr,Jb=12UTBU-1U+12VTBV-1V.For simplicity, here, we assume BU=BV=B, and no covariance between the two horizontal velocity components is considered. The background error covariance matrix B takes the following form [33]:(17)Bi,j=ahexp-Δxij2Lx2-Δyij2Ly2-Δzij2Lz2,where Lx, Ly, and Lz are characteristic length scales which reflect the extent of spatial correlation; xy and z are coordinates, and i and j are indexes of grid point; and ah is the first-guess error variance. Jc and Jr are similar to formulas (2) and (3), respectively, but with an omitted superscript n. Strong constraint is also imposed on U, V, and W component by using discretized continuity equation. To make a comparison with STMAS, the dense conventional data and Doppler radar radial velocity data are all used. The RMSE is shown in Table 1.
Based on Gaussian distribution, the traditional 3D-Var using correlation scale usually constructs the background error covariance matrix by an empirical correlation scale. Therefore, the traditional 3D-Var with a certain correlation scale only can analyze this kind of scale information. However, the “truth” typhoon field in this study contains two different spatial scales wind speed information (~35 km and ~10 km). The traditional 3D-Var with 50 km horizontal correlation scale can only capture the main pattern of this typhoon field (i.e., the long wave information) but lose the small-scale information and produce a smooth analysis (Figure 5(a)). The traditional 3D-Var with 25 km or 12.5 km horizontal correlation scale can capture some detailed information of the long wave as well as major short wave information, but detailed short wave information cannot be well analyzed (Figures 5(b) and 5(c)). On the contrary, the traditional 3D-Var with 5 km correlation scale can analyze short wave features but incorrectly treats long wave information. This leads to an erroneous analysis (Figure 5(d)). The two scales of this true typhoon field are 35 km and 10 km, respectively, and the large-scale (35 km) component covers most of study domain. In the experiment T25_RADAR_DENSE, the horizontal correlation scale is 25 km which is the closest to the large-scale (35 km) of the true typhoon field among all these traditional experiments; therefore, the horizontal wind velocity analysis is the best.
Similar to Figure 4(b) but with the results of traditional 3D-Var with correlation scales being Lx=50 km, Ly=50 km, and Lz=5 km (a); Lx=25 km, Ly=25 km, and Lz=2.5 km (b); Lx=12.5 km, Ly=12.5 km, and Lz=1.25 km (c); and Lx=6.25 km, Ly=6.25 km, and Lz=0.625 km (d) (unit: m/s).
From the vertical velocity distribution as well as wind vector shown in Figure 6, we can further verify the above argument. Near the typhoon center, there are two branches of vertical circulations on each side: one is narrow and the other is broad, corresponding to four upwelling zones in this section. Using the dense conventional data and Doppler radar radial velocity data, the STMAS can analyze these sets of vertical circulations and corresponding upwelling zones (Figure 6(a)). However, the traditional 3D-Var with horizontal correlation scales being 50 km, 25 km, or 12.5 km can only analyze two (Figure 6(b)) or three (Figures 6(c) and 6(d)) upwelling zones, since these kinds of correlation scales filter out the small-scale information and merge these upwelling zones. Although the traditional 3D-Var with 6.25 km horizontal correlation scale can distinguish the middle two upwelling zones, it incorrectly treats the large-scale information (Figure 6(e)).
Vertical velocity section (color shade) and wind vector in this section across the typhoon center. (a) STMAS; (b) tradition 3D-Var with Lx=50 km, Ly=50 km, and Lz=5 km; (c) tradition 3D-Var with Lx=25 km, Ly=25 km, and Lz=2.5 km; (d) tradition 3D-Var with Lx=12.5 km, Ly=12.5 km, and Lz=1.25 km; (e) tradition 3D-Var with Lx=6.25 km, Ly=6.25 km, and Lz=0.625 km (unit: m/s).
4. Conclusions and Discussions
Within an idealized simulation framework, the role of 3D Doppler radar radial velocity data for reconstructing 3D typhoon structures has been examined using the Space and Time Mesoscale Analysis System (STMAS). A two-scale simulated “typhoon” field is constructed and referred to as “truth,” from which randomly distributed conventional wind data and 3D Doppler radar radial wind data are generated. These data are used to reconstruct the synthetic 3D “typhoon” structure by the STMAS or the traditional 3D variational (3D-Var) analysis. The degree by which the “truth” 3D typhoon structure is recovered is an assessment of the impact of data type or analysis scheme being evaluated. The effects of weak or strong constraint on STMAS analysis have also been examined. We found that (1) the STMAS is superior to traditional 3D-Var for reconstructing the 3D typhoon structure, since the STMAS can retrieve multiscale information from observational network. (2) The radial velocity data provide additional useful information for the STMAS to reconstruct the detailed structure of 3D typhoon field. (3) Compared to a weak constraint, the strong constraint STMAS can produce better analyses on both horizontal and vertical velocities of the 3D typhoon structure.
This study gives us promising results. Challenges still remain when 3D radar radial velocity data are assimilated for the reconstruction and initialization of real typhoon structures in the future. First, given the fact that the real atmosphere is compressive, the wind vector may not satisfy the nondiffusivity continuity equation used in this study. Therefore, a full continuity equation should be used to make a more general strong constraint in future study. Second, the model error has not been taken into account in this study. The influence of model errors on typhoon reconstruction and initialization has to be addressed and how to deal with model errors could be an important research topic in the follow-up studies.
Competing Interests
The authors declare that there are no competing interests regarding the publication of this paper.
Acknowledgments
This research was jointly supported by grants of National Basic Research Program (2013CB430304), National Natural Science Foundation (41376013, 41376015, 41306006, 41541041, and 41506039), National High-Tech R&D Program (2013AA09A505), and National Programme on Global Change and Air-Sea Interaction (GASI-01-01-12 and GASI-IPOVAI-04) of China.
SunJ.CrookN. A.Real-time low-level wind and temperature analysis using single WSR-88D data200116111713210.1175/1520-0434(2001)016<0117:rtllwa>2.0.co;22-s2.0-0035245587WeygandtS. S.ShapiroA.DroegemeierK. K.Retrieval of model initial fields from single-Doppler observations of a supercell thunderstorm. Part I: single-Doppler velocity retrieval2002130343345310.1175/1520-0493(2002)130<0433:romiff>2.0.co;22-s2.0-0036509436WeygandtS. S.ShapiroA.DroegemeierK. K.Retrieval of model initial fields from single-Doppler observations of a supercell thunderstorm. Part II: thermodynamic retrieval and numerical prediction2002130345447610.1175/1520-0493(2002)13060;0454:romiff62;2.0.co;22-s2.0-0036509429LiouY.-C.ChiouJ.-L.ChenW.-H.YuH.-Y.Improving the model convective storm quantitative precipitation nowcasting by assimilating state variables retrieved from multiple-Doppler radar observations2014142114017403510.1175/mwr-d-13-00315.12-s2.0-84910020647SunJ.CrookN. A.Dynamical and microphysical retrieval from Doppler radar observations using a cloud model and its adjoint—part I: model development and simulated data experiments199754121642166110.1175/1520-0469(1997)054<1642:damrfd>2.0.co;22-s2.0-0031399793SunJ.CrookN. A.Dynamical and microphysical retrieval from Doppler radar observations using a cloud model and its adjoint. Part II: retrieval experiments of an observed Florida convective storm199855583585210.1175/1520-0469(1998)05560;0835:damrfd62;2.0.co;22-s2.0-0031806784GaoJ.XueM.ShapiroA.DroegemeierK. K.A variational method for the analysis of three-dimensional wind fields from two Doppler radars199912792128214210.1175/1520-0493(1999)127<2128:avmfta>2.0.co;22-s2.0-0033393125LindskogM.JärvinenH.MichelsonD. B.Development of radar wind data assimilation for the HIRLAM 3DVar200252GaoJ.-D.XueM.BrewsterK.DroegemeierK. K.A three-dimensional variational data analysis method with recursive filter for Doppler radars200421345746910.1175/1520-0426(2004)021<0457:atvdam>2.0.co;22-s2.0-1842735150XiaoQ.KuoY.-H.SunJ.LeeW.-C.LimE.GuoY.-R.BarkerD. M.Assimilation of Doppler radar observations with a regional 3DVAR system: impact of Doppler velocities on forecasts of a heavy rainfall case200544676878810.1175/jam2248.12-s2.0-24044451491GaoJ.XueM.LeeS.-Y.ShapiroA.XuQ.DroegemeierK. K.A three-dimensional variational single-Doppler velocity retrieval method with simple conservation equation constraint2006941–4112610.1007/s00703-005-0170-72-s2.0-33750517940HuM.XueM.BrewsterK.3DVAR and cloud analysis with WSR-88D level-II data for the prediction of the Fort Worth, Texas, tornadic thunderstorms. Part I: cloud analysis and its impact2006134267569810.1175/mwr3092.12-s2.0-33645323212HuM.XueM.GaoJ.BrewsterK.3DVAR and cloud analysis with WSR-88D level-II data for the prediction of the Fort Worth, Texas, tornadic thunderstorms—part II: impact of radial velocity analysis via 3DVAR2006134269972110.1175/mwr3093.12-s2.0-33645318980XiaoQ.SunJ.Multiple-radar data assimilation and short-range quantitative precipitation forecasting of a squall line observed during IHOP_20022007135103381340410.1175/mwr3471.12-s2.0-35748980004SugimotoS.CrookN. A.SunJ.XiaoQ.BarkerD.Assimilation of Doppler radar data with WRF 3DVAR: evaluation of its potential benefits to quantitative precipitation forecasting through Observing System Simulation Experiments200913740114029ZhangF. Q.WengY.SippelJ. A.MengZ.BishopC. H.Cloud-resolving hurricane initialization and prediction through assimilation of doppler radar observations with an ensemble Kalman filter200913772105212510.1175/2009mwr2645.12-s2.0-68249160487ZhaoK.XueM.Assimilation of coastal Doppler radar data with the ARPS 3DVAR and cloud analysis for the prediction of Hurricane Ike (2008)20093612L1280310.1029/2009GL0386582-s2.0-68749099447SchenkmanA. D.XueM.ShapiroA.BrewsterK.GaoJ.The analysis and prediction of the 8-9 May 2007 Oklahoma tornadic mesoscale convective system by assimilating WSR-88D and CASA radar data using 3DVAR2011139122424610.1175/2010MWR3336.12-s2.0-79951649710LiY. Z.WangX. G.XueM.Assimilation of radar radial velocity data with the wrf hybrid ensemble–3DVAR system for the prediction of hurricane ike (2008)2012140113507352410.1175/mwr-d-12-00043.12-s2.0-84867026914DuN.XueM.ZhaoK.MinJ.Impact of assimilating airborne Doppler radar velocity data using the ARPS 3DVAR on the analysis and prediction of Hurricane Ike (2008)201211717D1811310.1029/2012JD0176872-s2.0-84867009410AbhilashS.SahaiA. K.MohankumarK.GeorgeJ. P.DasS.Assimilation of Doppler weather radar radial velocity and reflectivity observations in WRF-3DVAR system for short-range forecasting of convective storms2012169112047207010.1007/s00024-012-0462-z2-s2.0-84868098928DongJ.XueM.Assimilation of radial velocity and reflectivity data from coastal WSR-88D radars using an ensemble Kalman filter for the analysis and forecast of landfalling hurricane Ike (2008)201313967146748710.1002/qj.19702-s2.0-84875381111ZhuL.WanQ.ShenX.MengZ.ZhangF.WengY.SippelJ.GaoY.ZhangY.YueJ.Prediction and predictability of high-impact Western Pacific landfalling Tropical Cyclone Vicente (2012) through convection-permitting ensemble assimilation of doppler radar velocity20161441214310.1175/MWR-D-14-00403.12-s2.0-84957837348ShenF.MinJ.XuD.Assimilation of radar radial velocity data with the WRF Hybrid ETKF-3DVAR system for the prediction of Hurricane Ike (2008)201616912713810.1016/j.atmosres.2015.09.0192-s2.0-84944930717XieY. F.KochS. E.McGinleyJ. A.AlbersS.WangN.A sequential variational analysis approach for mesoscale data assimilationProceedings of the 21st Conference on Weather Analysis and Forecasting/17th Conference on Numerical Weather Prediction, 15B.72005Washington, DC, USAAmerican Meteorological Societyhttp://ams.confex.com/ams/pdfpapers/93468.pdfXieY.KochS.McGinleyJ.AlbersS.BieringerP. E.WolfsonM.ChanM.A space-time multiscale analysis system: a sequential variational analysis approach201113941224124010.1175/2010mwr3338.12-s2.0-79955127661HeZ. J.XieY. F.LiW.LiD.HanG. J.LiuK. X.MaJ. R.Application of the sequential three-dimensional variational method to assimilating SST in a global ocean model20082561018103310.1175/2007jtecho540.12-s2.0-51749107003LiW.XieY.HeZ.HanG.LiuK.MaJ.LiD.Application of the multigrid data assimilation scheme to the China seas' temperature forecast200825112106211610.1175/2008JTECHO510.12-s2.0-65549083111LiW.XieY. F.DengS.-M.WangQ.Application of the multigrid method to the two-dimensional Doppler radar radial velocity data assimilation201027231933210.1175/2009jtecha1271.12-s2.0-77952479969LiW.XieY. F.HanG. J.A theoretical study on the multi-grid data assimilation scheme using a simple bilinear interpolation scheme20133238087FletcherR.19872ndChichester, UKJohn Wiley & SonsMR955799ByrdR. H.LuP.NocedalJ.ZhuC. Y.A limited memory algorithm for bound constrained optimization19951651190120810.1137/0916069MR1346301ZBL0836.65080DerberJ.RosatiA.A global oceanic data assimilation system19891991333134710.1175/1520-0485(1989)01960;1333:agodas62;2.0.co;2