Next Article in Journal
Error Characteristics and Scale Dependence of Current Satellite Precipitation Estimates Products in Hydrological Modeling
Next Article in Special Issue
Assimilating FY-4A Lightning and Radar Data for Improving Short-Term Forecasts of a High-Impact Convective Event with a Dual-Resolution Hybrid 3DEnVAR Method
Previous Article in Journal
Optical Classification of Lower Amazon Waters Based on In Situ Data and Sentinel-3 Ocean and Land Color Instrument Imagery
Previous Article in Special Issue
Climatology of Convective Storms in Estonia from Radar Data and Severe Convective Environments
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assimilation of Polarimetric Radar Data in Simulation of a Supercell Storm with a Variational Approach and the WRF Model

1
Hubei Key Laboratory for Heavy Rain Monitoring and Warning Research, Institute of Heavy Rain, China Meteorological Administration, Wuhan 430205, China
2
National Severe Storms Laboratory/NOAA, Norman, OK 73072, USA
3
Cooperative Institute for Mesoscale Meteorological Studies, The University of Oklahoma, Norman, OK 73072, USA
4
Key Laboratory of Mesoscale Severe Weather/Ministry of Education and School of Atmospheric Sciences, Nanjing University, Nanjing 210023, China
5
School of Meteorology, The University of Oklahoma, Norman, OK 73072, USA
*
Author to whom correspondence should be addressed.
Remote Sens. 2021, 13(16), 3060; https://doi.org/10.3390/rs13163060
Submission received: 17 June 2021 / Revised: 31 July 2021 / Accepted: 1 August 2021 / Published: 4 August 2021
(This article belongs to the Special Issue Atmospheric Radar for Severe Weather Research)

Abstract

:
Polarimetric radar data (PRD) have potential to be used in numerical weather prediction (NWP) models to improve convective-scale weather forecasts. However, thus far only a few studies have been undertaken in this research direction. To assimilate PRD in NWP models, a forward operator, also called a PRD simulator, is needed to establish the relation between model physics parameters and polarimetric radar variables. Such a forward operator needs to be accurate enough to make quantitative comparisons between radar observations and model output feasible, and to be computationally efficient so that these observations can be easily incorporated into a data assimilation (DA) scheme. To address this concern, a set of parameterized PRD simulators for the horizontal reflectivity, differential reflectivity, specific differential phase, and cross-correlation coefficient were developed. In this study, we have tested the performance of these new operators in a variational DA system. Firstly, the tangent linear and adjoint (TL/AD) models for these PRD simulators have been developed and checked for the validity. Then, both the forward operator and its adjoint model have been built into the three-dimensional variational (3DVAR) system. Finally, some preliminary DA experiments have been performed with an idealized supercell storm. It is found that the assimilation of PRD, including differential reflectivity and specific differential phase, in addition to radar radial velocity and horizontal reflectivity, can enhance the accuracy of both initial conditions for model hydrometer state variables and ensuing model forecasts. The usefulness of the cross-correlation coefficient is very limited in terms of improving convective-scale data analysis and NWP.

1. Introduction

The accuracy of the initial conditions of numerical weather prediction (NWP) models has a great impact on their forecast results, especially for high-resolution NWP. To build up the initial hydrometeor information which is the key to successful simulating and forecasting of convective storms, the inevitable process known as “spin-up” is needed for NWP models [1]. In order to reduce the negative influence of “spin-up” on the forecast skill of short-lived convective storms, Doppler weather radar observations, including the horizontal reflectivity (ZH) and radial velocity (Vr), are commonly used to improve the initial condition with advanced data assimilation (DA) methods, such as the three-dimensional variational (3DVAR) method [2,3,4,5,6,7], the four-dimensional variational (4DVAR) method [8,9,10,11,12], and the ensemble Kalman filter (EnKF) [13,14,15,16,17,18].
The detection of additional measurements representing the microphysical characteristics of precipitation systems, such as the differential reflectivity (ZDR), cross-correlation coefficient (ρhv), and differential phase (ϕHV), which is the range integral of the specific differential phase (KDP), benefits from polarimetric radar technology that has been widely concerned and developed rapidly after decades of effort. Nowadays, polarimetric radars are deployed and applied globally; in the U.S., for instance, the national radar network of WSR-88D has been upgraded with dual-polarization capability, and polarimetric radar data (PRD) are available nationally [19,20]. The PRD has been also successfully used in quantitative precipitation estimation (QPE) [21,22,23,24], hydrometeor classification (HC) [25,26,27], microphysics retrieval [28,29,30,31] and severe weather identification [32,33,34,35]. These mentioned advantages prompt us to further explore the possibility of assimilating PRD into NWP models to improve weather forecasts, especially for the prediction of convective-scale weather.
A forward observation operator, also known as the PRD simulator, is needed to establish the relation between model state variables and PRD for their assimilation in NWP models. Many polarimetric radar rain estimators have been developed and used for QPE [21,29,30]. These relations have also been extended for ice phases including snow and ice [36,37]. Because of the complex nature and nonlinear interaction in mixed phases, these aforementioned single species estimators are unable to fully establish the relation between the multiple species of model state variables and PRD for DA use.
To use PRD in the DA process of improving the accuracy of the initial conditions, or obtain the simulated PRD from model microphysics state variables, Pfeifer et al. [38] and Jung et al. [39] (hereafter J08) developed their own PRD simulators, respectively. Both studies include the contribution of ice and melting species to polarimetric parameters, but the former considers the attenuation effect of radar beams as well, while the latter emphasizes the important influence of melting ice particles. Then, Jung et al. [40] first directly assimilated PRD using the ensemble square-root Kalman filter (EnSRF) with a single-moment (SM) microphysical parameterization (MP) scheme, and then attained the improved analysis of the storm when polarimetric variables were assimilated. Applying rigorous scattering calculations via the T-matrix method, a more accurate and numerically integrated observation operator was developed by Jung et al. [41] (hereafter J10). Compared the simulation performance of J10 under SM and double-moment (DM) MP schemes, Jung et al. [42] found that some realistic polarimetric signatures reported in observational studies such as ZDR arc, midlevel ZDR, and ρhv ring [23,33] can only be reproduced by employing a DM MP scheme.
Li and Mecikalski [43,44] used empirical relations and retrieved the rain water mixing ratio from ZDR [45] and KDP [46], respectively, and then assimilated them in the Weather Research and Forecasting (WRF) model using the WRF 3DVAR system. The consequences demonstrate the benefits of utilizing polarimetric variables to improve the initial condition and even short-term forecasts for mesoscale convective systems (MCSs). Then, Li et al. [47] developed an ice-phase operator embedded in the WRF 3DVAR system for assimilation of KDP, which can also additionally adjust solid water (cloud ice and snow) through DA. Consequently, a positive impact of extra KDP assimilation on analysis fields of rainwater in the lower troposphere and snow in the mid- to upper troposphere for an MCS was found. Aimed at applying J08 under the variational DA framework, Wang and Liu [48] rebuilt an operator named RadZIceVar (a package containing a reflectivity forward operator with its associated TL/AD operators together) and indicated its successful implement in DA use and significant improvement on short-term (2-5 h) precipitation forecast. However, the direct assimilation was applicable for ZH only, and assimilation of polarimetric variables was not supported.
Meanwhile, using J08 with SM and DM MP schemes, Putnam et al. [16] noted several microphysics biases with the help of simulated PRD through qualitative and quantitative verification of ensemble probabilistic forecasts in a real case. Subsequently, by establishing look-up tables to increase the calculation efficiency of J10 during the DA procedure, Putnam et al. [17] first implemented the direct assimilation of ZDR through the EnKF method for a real supercell case in the literature. Through a series of observing system simulation experiments (OSSEs) of an idealized supercell storm, Zhu et al. [18] investigated the impact of assimilating ZDR within an EnKF framework by using J10 coupled with a DM MP scheme, and also discussed the potential influence of updating hydrometeor number concentrations for DA effect. Their conclusions show that the assimilation of ZDR can improve the accuracy of analyzed hydrometeor fields in terms of pattern and intensity, and that updating the number concentrations with mixing ratios is important for deciding whether the benefit of assimilating ZDR can be preserved. However, they are all limited to the use of partial PRD for assimilation, and significant errors in the magnitude of analyzed PRD do exist, indicating that the operational application of PRD assimilation remains challenging.
In addition, some polarimetric signatures found in convective storms provide the opportunity for indirect assimilation of PRD. Using the Advanced Regional Prediction System (ARPS) cloud analysis module [5,49,50], Carlin et al. [51] managed to assimilate PRD through the detection of ZDR column signatures [33,52], which often occur in the strong updraft region within intense thunderstorms, to modify the temperature and moisture of the ARPS model. Encouraging results were obtained with improved analyses (more coherent updraft) and forecasts (more realistic reflectivity structures with better skill scores) compared to the original cloud analysis scheme without polarimetric information included. There are limits to the widespread use of this approach, however, due to the absence of ZDR columns in relatively weak precipitation systems and its high dependence on empirical relationships between model state variables and PRD.
Recently, Zhang et al. [53] (hereafter Z21) developed a new set of simplified and parameterized PRD observation operators which can link NWP model state variables and PRD for DA use, and verified its validity and applicability by applying an ideal case and a real case respectively. They found that the parameterized operators yield results consistent with that of rigorous calculation in J10. However, Z21 increased computational efficiency significantly, resulting in costs less than one percent of the computing time of J10 to complete the same task, which makes it more suitable for DA applications. Thus, in this study, the applicability of Z21 is investigated under a variational DA framework, and the different DA effects of each polarimetric variable are studied by conducting various OSSEs.
The rest of this paper is organized as follows. The new parameterized PRD operators are briefly reviewed in Section 2 and the overview of the 3DVAR DA system is provided in the following section. In Section 4, the experimental design is given, and some preliminary experiment results are discussed in Section 5. A summary and conclusions for this study are presented in the last section.

2. Polarimetric Radar Observation Operator

2.1. Microphysics Models and Parametrization

SM and DM MP schemes are commonly used in NWP models [54,55,56,57,58]. In SM simulations of NWP models, water mixing ratios ( q ), which is directly relevant to water content ( W = ρ a q , where ρ a is the air density), is the only prognostic physics variable for hydrometeor physics. For DM MP schemes, the number concentration ( N t ) and the water mixing ratio ( q ) are both predicted variables. For a given two-parameter drop/particle size distribution (DSD/PSD) model and a specified hydrometeor x , both N t x and q x can be converted to mass-weighted diameter ( D m x = 4 ρ a q x π ρ x N t x 1 / 3 ) and water content ( W x = ρ x q x ); then, many other related quantities such as DSD/PSD parameter, rainfall rate and reflectivity are calculated by these two state variables. The detailed description of the derivation can be found in Z21.

2.2. Parameterized PRD Operators

Polarimetric radar variables (ZH, ZDR, KDP, ρhv) are expressed by the integrals of DSD/PSD weighted by the scattering amplitudes [59,60]. The scattering amplitudes are normally calculated by the T-matrix method, and the integrals are calculated numerically for accurate results. This was done in J10. However, the complex integral form and resulting expensive computational cost make it not convenient for variational DA use, which usually needs the operators to be efficient and easy calculation of the derivative.
Mahale et al. [61] represented polarimetric variables as simple functions of water mixing ratio q r and mass/volume-weighted diameter D m in the polynomial form for rain. In the case of ice phases and mixtures such as dry/melting snow, hail, and graupel, the calculation and parameterization are more complicated than those of rain induced by the increased variability in density during the melting stage and irregular shape, as well as the orientation of the particles. For simplification, a simple function of the percentage of melting ( γ x ) as following is used in Z21 to calculate the particle density ρ x for the melting hydrometeor species x:
ρ x = ρ d x 1 γ x 2 + ρ w γ x 2 ,
where γ x = q r / q r + q x (the same as J08), ρ d x stands for the dry/pure state of ice hydrometeor corresponding to the melting state, and ρ w is the density of water. Except for the mean axis ratio and standard deviation of the canting angles, Z21 makes some adjustments to allow a larger dynamic range of ZDR and ρhv. The shape and orientation of hydrometer particles follow the modeling and representation documented in J08. Similar to rain, for a given ice or mixed phase species x, including snow, hail, graupel, and their melting parts, PRD variables are calculated for a set of D m at a given γ x , and then parameterized as functions of D m as follows:
Z h x Z x α Z 0 γ x + α Z 1 γ x D m + α Z 2 γ x D m 2 + α Z 3 γ x D m 3 2 ,
Z dr x α d 0 γ x + α d 1 γ x D m + α d 2 γ x D m 2 ,
K DP x ρ a q w x ( α K 0 γ x + α K 1 γ x D m + α K 2 γ x D m 2 ) / ρ x ,
ρ hv x α ρ 0 γ x + α ρ 1 γ x D m + α ρ 2 γ x D m 2 ,
where the reflectivity factor is Z x = 11.25 × 10 3 ρ a q x π ρ x D m 3 ( D m , unit: mm). The fitting coefficients associated with γ x are provided in Table 1 and Table 2 of Z21, and more detailed descriptions for calculations, assumptions, and result analysis can also be found there.

3. The 3DVAR DA System

To test the assimilation performance of Z21 in a variational DA scheme, in this study, a 3DVAR DA system designed especially for radar data assimilation at the convective scale is employed. This 3DVAR system was initially developed at the Center for Analysis and Prediction of Storms (CAPS) and subsequently refined at National Severe Storms Laboratory (NSSL). The standard 3DVAR cost function can be written as follows [3]:
J x = 1 2 x x b T B 1 x x b + 1 2 H x y o T R 1 H x y o + J c x .
The cost function defined above contains three terms on the right-hand side. The first term is the background term, which defines the departure of the analysis vector x from the background vector x b weighted by the inverse of the NWP model background error covariance matrix B , where the analysis vector x contains the following variables: the three wind components ( u , v , and w ) and the mixing ratios of hydrometeors (rain water, snow, hail, graupel, cloud water, cloud ice). The background error covariance matrix is modeled by the product of a diagonal matrix of the standard deviation of the background error and a spatial recursive filter [3]. The standard deviations for the model variables are derived from the statistics of the Rapid Update Cycle (RUC) model’s 3 h forecasts over several years from 2001 to 2004 [62]. The second term is the observation term, which defines the departure of the analysis vector x from the observation vector y o , where H x is the forward operator which builds a “bridge” between the model state variables and observation data. In this study, the new PRD operators (Z21), which can use both the mixing ratio and the number concentration from a DM MP scheme, will be discussed later. In the second term, R is called the observation error covariance matrix which includes both instruments and representativeness errors. The observation errors’ standard deviation for PRD variables (Vr, ZH, ZDR, KDP, and ρhv) is specified empirically as 2 m/s, 3 dBZ, 0.5 dB, 0.5 °/km, and 0.1 respectively. Finally, as the third term on the right-side, J c x represents the dynamic or equation constraints. This term has proved to be very important, particularly for this convective-scale 3DVAR system. For example, the mass continuity equation imposed as a weak constraint is helpful in producing more suitable wind analysis [2,3].
Radar radial velocity Vr is regarded as a crucial variable in the DA procedure [3,4,10] because it provides important wind field information and dynamic characteristics of the precipitation event. In this 3DVAR DA system, the forward observation operator for radial velocity considering the effects of the Earth’s curvature is expressed as:
  V r = d h d r w + d s d r u sin φ + v cos φ ,
where h and s is radar beam height and distance, respectively, over and along the curving Earth’s surface, r is the radial slant range, and φ is the radar azimuth angle. The 4/3-effective Earth radius mode is assumed to follow for the propagation of radar beams [59]. The effect of hydrometeor fall velocity is neglected here, and will be considered in the follow-up real data case study.
For PRD, after the calculations for each species x available in specified MP scheme through Equations (2)–(5), the final variables for the pixel are obtained by summing the contributions from all species as follows:
Z H = 10 log [ Z h x ] ,
Z DR = 10 log Z h x Z h x / Z dr x ,
K DP = K DP x ,
ρ hv = Z h x Z d r 1 / 2 x ρ hv x Z h x Z dr 1 / 2 x .
For the 3DVAR method used here, considering the variations of observational spacing existing among different observations, multiple analysis passes are usually used to analyze different observation types with different filter scales. Indeed, applying multiple-pass with a recursive filter is proven superior to the common single-pass method theoretically [6,43]. Hence, in this study, each PRD variable is analyzed in a separate pass.
The new operators have been proved efficiently in PRD simulation and therefore particularly suitable for DA purposes (Z21), including both variational and ensemble DA. In such applications, the TL/AD operators, coupled with the forward operators, are required for the variational assimilation, which are not necessarily found in the EnKF [63]. As a consequence, complex operators such as those used in J10 are more applicable to be employed through the EnKF method for DA use. Despite this difficulty, due to the simple polynomial form of Z21, the programming of corresponding TL/AD operators is implemented smoothly according to the rules of adjoint code construction [64,65].
The AD operator is the transpose of the TL operator. Therefore, how well TL operator approximates the nonlinear forward operator determines how valid the gradient of the cost function is. Following [66], to verify the correctness of the TL operator, one can define a function R(a), given by
R a = f x 0 + a d x f x 0 a f d x ,
where f x and f x are forward, tangent linear operator respectively. The x 0 represents a random picked base state and d x represents a perturbation state. Here, R(a) is calculated over all points in the domain. According to the theory, if R(a) in Equation (12) is close to 1 for any small values of a (see Table 1), meaning that the TL operator approximates the nonlinear forward operator well. The test results of the TL operator are listed in Table 1, as a keeps decreasing from 1 to 1E-14, the R(a) values for all polarimetric variables (R_ZH, R_ZDR, R_KDP and R_ρhv) are approaching 1 uniformly and steadily in the case of double precision. However, when a goes down to an extra small value, the R(a) starts to increase, which is caused by the rounding error of the calculator floating-point arithmetic operation. The above outcomes enable us to conclude that the assimilation system with the Z21 integrated passes the validity test and that the subsequent experimental research can be carried out reasonably.

4. Experimental Design

As discussed previously, PRD shows potential for enhancing storm simulations in terms of intensity and structure and ensuing short-range forecasts [16,17,18,40,42,43,44,51]. However the forward operators used thus far are either too simple [43,44] or too complicated [39,41]. The recently proposed parameterized PRD operators in Z21 can be used in different microphysics schemes and are relatively efficient.
In this study, a fully compressible and non-hydrostatic Advanced Research Weather Research and Forecasting (WRF-ARW) model (version 3.7.1) is used in a three-dimensional space for a supercell storm simulation. The horizontal resolution is 1 km with 121 grid points in the east-west and north-south directions. Simultaneously, 51 stretched levels up to 20 km (50 hPa, approximately) above ground level (AGL) are chosen. Open boundary conditions for the lateral and Rayleigh damping for the top boundary are used in this idealized case. The whole integration of the WRF-ARW model for simulation is 3 hours.
A modified sounding from a classic supercell event that occurred on 20 May 1977 in Del City, Oklahoma is used to provide the environmental field of the storm. A thermal bubble with a potential temperature perturbation of 3 K is added in the truth field to initiate convection [14,67,68,69]. This warm bubble is centered at the location of (x = 90 km, y = 40 km, z = 1.5 km) and has a 10 km horizontal radius and 1.5 km vertical radius inside the model domain. The Milbrandt and Yau scheme [55], a DM MP scheme consisting of the predicted mixing ratios of six-class hydrometeors, such as cloud water ( q c ), cloud ice ( q i ), rain water ( q r ), snow ( q s ), hail ( q h ), and graupel ( q g ), and their corresponding number concentrations ( N c , N i , N r , N s , N h , and N g , respectively), is applied in description of microphysical processes in the idealized case. The time step for the integration of the prognostic equations is set to be 3 seconds. The standard 1.5-order TKE closure scheme is chosen for the turbulence parameterization, and other configurations are the same as that of Ge et al. [67].
During the 3 h truth simulation, the cloud forms around 10 min, rain water appears at 15 min, ice hydrometeors are generated at 20 min, and a single convective storm develops in the first 30 min (not shown). The storm reaches its mature stage at 45 min, then starts to split, evolves and slightly weakens. At 2 hours into the model integration, the right-splitting cell tends to dominate, as indicated by a clear hook echo and strong updraft.
Figure 1 shows the cycled DA and forecast scheme employed in OSSEs. Based on the truth simulation with WRF, all polarimetric radar observations (Vr, ZH, ZDR, KDP, and ρhv) can be simulated and used in the iterations of 3DVAR analysis to update the model state variables ( u a , v a , w a , and q x a ). After the analysis step, a 5 min forecast procedure is performed and the result is the background field ( u b , v b , w b , and q x b ) as the input of the DA analysis step in the next cycle. Considering the mixing ratios of precipitation are too small at the early stage of storm initiation, the idealized data are generated after the 30 min of model integration. Two assumed S-band radars are located in the southeast and southwest corner, respectively, and are performed on the VCP11 scan mode of WSR-88D radar (14 tilts from 0.5° to 19.5°). Idealized PRD is derived from the truth simulation output using the forward operators in Equations (7)–(11).
All OSSEs presented in this study are listed in Table 2. The control run assimilating Vr and ZH data, labeled as ExpVrZh, is considered as a reference for comparison purposes. Sensitivity experiments ExpVrZhZdr, ExpVrZhKdp, and ExpVrZhRhv are the same as the control run but with the additional assimilation of ZDR, KDP, and ρhv, respectively. ExpVrZhPol is a special run with all PRD assimilated. These sensitivity experiments are performed to help understand the impact of PRD on the analysis of hydrometer variables. Therein, the number concentrations of hydrometeors will not be updated during the DA cycles on account of the complexity introduced by their wide dynamic range (from 0 to above 1E12 m-3) and the problem of a significant increase of the degrees of freedom. The explanation of all abbreviations used in this paper can be found in Appendix A.

5. Results of 3DVAR Analysis

5.1. The Root Mean Square Error Analysis

The Root Mean Square Errors (RMSEs) are commonly used as a criterion for the accuracy judgment of DA [14,18,40,68,69]. To be consistent with the previous requirement of ZH assimilation, the RMSEs of analysis and forecast displayed in Figure 2 are only calculated over precipitation grids with ZH exceeding 10 dBZ.
At a first glance of Figure 2, for all model state variables and simulated radar observations, the RMSEs of analysis and forecast have a generally decreasing tendency over time for each OSSE. Since number concentrations are not updated with mixing ratios during the DA process, there may be a potential imbalance between them, leading to fluctuations in the early DA stages. Compared to ExpVrZh (black line), ExpVrZhZdr (red line) produces slightly better analyses and forecasts for most variables, where more pronounced improvements occur in the last few cycles, especially in ice hydrometeor species q i (Figure 2e), q h , and q g (Figure 2h). Since KDP is most sensitive to the presence of liquid water, the RMSE of q r (Figure 2f) of ExpVrZhKdp (blue line) has the most obvious difference. The better convergence through the whole DA cycle suggests that the satisfied rain water field much closer to the truth can be analyzed with additional KDP assimilation. In fact, as with ExpVrZhZdr, it also plays a positive role in obtaining a smaller RMSE in the last few cycles for q i , q h , and q g . For ExpVrZhRhv (cyan line), the relatively distinct variations only occur in the first or last few cycles, but as a whole, its RMSEs almost follow that of ExpVrZh with no significant discrepancies, indicating the limited effect of ρhv assimilation.
In terms of observed radar variables, the RMSEs of ExpVrZhZdr converge better and faster, indicating its simulated PRD is superior to that of ExpVrZh, especially for ZDR (Figure 2j) and ZH (Figure 2i). The above conclusions agree well with those of Jung et al. [40], except for benefits shown earlier in DA cycles here. This is largely the same as the result of Zhu et al. [18]. Apparently, ExpVrZhKdp produces consistently better KDP (Figure 2k) in both analysis and forecast cycles, as expected. Moreover, there are some subtle advancements of ZH in the last few cycles, but no contributions to the simulation of ZDR and ρhv (Figure 2l). As confirmed by the variation of the RMSEs, the assimilation of ρhv cannot bring benefits in observation space as well. Although both ZH and ρhv can converge reasonably, the influence of ρhv assimilation is too small to be noticed, which may be caused by the relatively small variation range of itself (on the order of 0.01). Noted that the reverse changing problem (the RMSE decreases in the forecast cycle, but increases in the corresponding analysis cycle) occurring in ZDR and KDP can be resolved by assimilating the corresponding PRD variable.
In the aspect of RMSEs, different abilities and effects in improving the accuracy of the model state variables and simulation of PRD have been presented for each polarimetric variable. Superior and more balanced results are expected to be achieved by jointly assimilating all PRD. As shown in purple line in Figure 2, although some fluctuations arise in w (Figure 2c), q s (Figure 2g), q h and q g at certain cycles, the best convergent trends are obtained in ExpVrZhPol for almost all kinds of hydrometeors, especially in later DA cycles. For simulated PRD, ExpVrZhPol also acquires better RMSE convergence compared to other experiments, especially for ZDR and ρhv. Actually, with respect to most variables, the sharpest downward variation range occurring in the first cycle of ExpVrZhPol indirectly reveals the greatest potential capability of revising the model analysis field.

5.2. Evaluation of PRD Assimilation

The impacts of PRD assimilation on the horizontal distribution of intensity and polarimetric characteristics of the simulated storm are further analyzed in Figure 3, Figure 4, Figure 5 and Figure 6. Compared with the true ZH (Figure 3a), all OSSEs represent the hook echo signature and the strong ZH core at the leading edge of the storm. Among them, ExpVrZh (Figure 3b) and ExpVrZhRhv (Figure 3e) have very similar results with obviously high ZH (greater than 15 dBZ overestimated) in the rear stratiform precipitation area (hereafter SPA). However, the ZH intensity described by ExpVrZhZdr (Figure 3c) is superior to them, both in the front convective precipitation area (CPA) and the rear SPA. However, the weak echo (<20 dBZ) on the northern edge is over-adjusted to deviate from the truth. On the contrary, although the weak echo is well analyzed, a much smaller ZH exists in the SPA of ExpVrZhKdp (Figure 3d). In terms of intensity and structure, the overall result of ExpVrZhPol (Figure 3f) is much closer to the truth, particularly in the echo edge and rear SPA, except for a certain overestimation in the CPA. There is a small ZDR patch (<1 dB) in the southeast corner of the main storm in ExpVrZh (Figure 4b), and also the same in ExpVrZhRhv (Figure 4e). Beyond that, the ZDR in the SPA is significantly smaller than the truth (Figure 4a), with the largest difference over 1.5 dB. ExpVrZhKdp (Figure 4d) has even worse simulation with much smaller ZDR. However, the intensity of ZDR in ExpVrZhZdr (Figure 4c) is much enhanced in the fore-mentioned two underestimated areas, particularly in the SPA. However, to some extent, overestimation persists in the weak echo at the edge of storm. ExpVrZhPol (Figure 4f) equally gives the best analysis among all OSSEs, except for the underestimation that occurred in the south-east of the storm. Although the dominant distribution of KDP is well exhibited in ExpVrZh (Figure 5b), ExpVrZhZdr (Figure 5c) and ExpVrZhRhv (Figure 5e) without KDP assimilation, the central strength of the leading edge and the hook echo are clearly weaker. The assimilation of KDP is beneficial that the simulated KDP shown in Figure 5d is perfectly matched with the truth (Figure 5a) in both intensity and structure. Undoubtedly, as a result of additional assimilation of KDP, ExpVrZhPol (Figure 5f) also has a similar good performance to ExpVrZhKdp. Based on the simulated result of ExpVrZh (Figure 6b) which is pretty close to the true ρhv (Figure 6a), ExpVrZhRhv (Figure 6e) mainly brings improvement in details, whereas ExpVrZhZdr (Figure 6c) and ExpVrZhKdp (Figure 6d) make further overestimation of ρhv in certain areas. Affected by ZDR and KDP, ExpVrZhPol (Figure 6f) does not exhibit better ρhv distribution, but instead has an underestimated area (<0.9).
To further evaluate the simulation of the vertical structure of the storm in different DA schemes, the vertical slices of analyses through the intensity core and hook echo region (black dashed line in Figure 3a) are shown in Figure 7, Figure 8, Figure 9 and Figure 10. For ZH, the whole structure and extended height of the strong echo are both close to the truth (Figure 7a), with the exceptions of the ring-shape core (>55 dBZ) in ExpVrZh (Figure 7b), ExpVrZhZdr (Figure 7c), and ExpVrZhRhv (Figure 7e). Additionally, ZH in these OSSEs are overestimated markedly in the middle level of the forward flank (5 km horizontally, the same hereafter), which is reduced properly in ExpVrZhKdp (Figure 7d). Additionally, the overestimation in the low level of the rear SPA (30~40 km) has been alleviated by assimilating ZDR, such that the intensity is much closer to the truth and the echo goes from ungrounded to grounded. ExpVrZhRhv shows no enhancement compared with ExpVrZh, or even worse results (overestimation in some areas, Figure 7e vs. Figure 7a and Figure 7b). More than that, the ahead-developing convective cell (~6 km) and inner strong echo of the main storm are both better described in ExpVrZhKdp and ExpVrZhPol (Figure 7f). In comparison with the truth (Figure 8a), the signature of the ZDR column (~5 km) related to the updraft within the storm is more evidently discerned in ExpVrZdr (Figure 8c) and ExpVrKdp (Figure 8d), whereas the sagged high ZDR area (~25 km) attributed to the existence of a small number of large raindrops is well retained in ExpVrZh (Figure 8b), ExpVrZdr, and ExpVrRhv (Figure 8e). Modifications of amplitude and shape of ZDR can be noticed in the tail (~35 km) in ExpVrZhZdr and ExpVrZhKdp. With respect to above features, the best analysis result in terms of intensity and structure is given in ExpVrZhPol (Figure 8f). ExpVrZhKdp reveals an extremely high degree of similarity with the truth (Figure 9d vs. Figure 9a), including the KDP core associated with high liquid water content (LWC) within the storm, an overshooting structure linked to the strong updraft rushing through the melting layer, and even the small KDP area hanging ahead. The effect of KDP assimilation is maintained in ExpVrZhPol (Figure 9f), showing almost identical distribution characteristics to that of ExpVrZhKdp. Comparatively, other OSSEs (Figure 9b,c,e) nearly fail to reproduce these mentioned signatures. The melting layer associated with the melting hydrometeors is precisely described by a ρhv drop in these OSSEs, but with a wider melting band and a few underestimated areas above. Other than a small positive contribution from the ρhv assimilation (Figure 10e) in the low-value area located in the middle level of the forward edge (~5 km) of the storm, there are no notable improvements in the distribution of analyzed ρhv for different OSSEs.

5.3. Evaluation of Hydrometeor Analysis

Figure 11, Figure 12, Figure 13 and Figure 14 show the mixing ratios of cloud ice q i , rain water q r , hail q h , and graupel q g in the truth and analysis fields of OSSEs, respectively, through identical vertical slices as in Figure 7. These OSSEs with PRD assimilation partially (mainly in the right part) and slightly (within 0.1 g/kg) strengthen q i more than that in ExpVrZh (Figure 11b), which is still much weaker compared to the truth (Figure 11a). Besides the basic distribution of q r , well depicted in OSSEs, the outstretched structure corresponding to the convection initiation ahead (~5 km) is better analyzed in ExpVrZhZdr (Figure 12c) and ExpVrZhKdp (Figure 12d). The latter better portrays the rear (>30 km) relatively small q r (0.1~0.2 g/kg), while the former is slightly over-adjusted (>0.2 g/kg). The overshooting pattern (~10 km) resulting from the strong updraft within the storm and suggesting the presence of super-cooled water is clearly displayed in ExpVrZhKdp as well, completely consistent with the truth (Figure 12a). However, the more obvious q h (Figure 13c) core (>1.8 g/kg) and enhanced q g (Figure 14c) in the middle- (>1.5 g/kg) and upper level (>4.0 g/kg) can be seen in ExpVrZhZdr. As expected, these mentioned benefits induced by ZDR or KDP assimilation are also retained in ExpVrZhPol. In addition, there are many other improvements such as the preferable descriptions of relatively small q i (<0.2 g/kg), q h (<0.6 g/kg), and q g (<4.5 g/kg) in the upper left corner f, Figure 13f and Figure 14f), the superior q r analysis much closer to the truth (Figure 12f vs. Figure 12a), and so on. However, apart from minor differences, all analyzed mixing ratios in ExpVrZhRhv are largely the same as the corresponding ones in ExpVrZh (Figure 11e, Figure 12e, Figure 13e and Figure 14e vs. Figure 11b, Figure 12b, Figure 13b and Figure 14b, respectively). Additionally, analyses of cloud water q c and snow q s yield similar results as other hydrometeor variables, that is, q c is similar to q r , and q s is similar to other ice species ( q i , q g , and q h ), which are not detailed here.
In general, the above results further demonstrate the potential positive effects of ZDR assimilation on the analysis of solid hydrometeor particles, more benefits of KDP assimilation in improving the simulation of rain water rather than ice hydrometeors, and the limited usefulness of ρhv in terms of enhancing model hydrometeor analysis. Attributing to the combined positive impacts from PRD assimilations, more comprehensive hydrometeor analysis results in both intensity and pattern can be obtained. The conclusions are also in accordance with the findings of RMSE analysis.
Different PRD variables have different information contents and error effects corresponding to each hydrometeor. ZDR increases as the raindrop size increases or ice particles melt. Therefore, assimilation of ZDR yields improved analysis of all hydrometeors. KDP depends on both the raindrop size and the number concentration, hence, its assimilation improves the analysis of q r . Since ρhv reduction occurs only in the melting/mixing regions which have a very small percentage of the storm and it has a small dynamic range with relatively large errors, the improvement with assimilation of ρhv is minimal.

5.4. Evaluation of Forecast

To investigate the improvement effect of modified initial conditions from different DA schemes on model predictions, the 1 h forecasts beginning at the 120 min of a storm are made at 15 min intervals for these OSSEs. In addition, the mean RMSEs of forecasts are calculated as seen in Figure 2 and displayed in Figure 15. Unsurprisingly, the increasing trend of RMSE with time is exhibited in all model state variables and simulated PRD, but the variation characteristics for different OSSEs are still existing diversely.
Taking all parameters into account, ExpVrZhPol (purple line) generally outperforms other OSSEs. Although the differences in the horizontal wind of u and v (Figure 15a,b) are almost indistinguishable, it can be seen upon careful inspection that ExpVrZhPol has a relatively low RMSE on the whole. Furthermore, except for the better performance of ExpVrZhZdr (red line) in the intermediate stage (30-45 min) of prediction, ExpVrZhPol also has the best overall prediction effect on vertical velocity w (Figure 15c). An analogous situation appears in the perturbation potential temperature θ (Figure 15d) and perturbation pressure p (Figure 15e) as well. For hydrometeors, the remarkable advancements resulting from combined assimilation of PRD can be found in cloud water q c (Figure 15g), rain water q r (Figure 15i), snow q s (Figure 15j), and graupel q g (Figure 15l) almost throughout the entire prediction period, and mainly in both ends of the period for cloud ice q i (Figure 15h) and hail q h (Figure 15k). While taken individually, with respect to RMSEs, ExpVrZhZdr has a good performance in ice hydrometeors ( q i , q h , and q g ) forecasts, suggesting the potential advantages of ZDR assimilation in increasing the prediction accuracy of ice species. Further, due to the high sensitivity of KDP to liquid water, ExpVrZhKdp (blue line) illustrates more superiorities in the forecasting of water vapor q v (Figure 15f), q c , and q r . Comparatively speaking, ExpVrZhRhv (cyan line) has fully comparable behaviors with ExpVrZh (black line), demonstrating the limited benefits of ρhv assimilation.
In observation space, ExpVrZhPol also manifests superior overall performance, thoroughly in ZH (Figure 15m), and mainly in the late forecast stage (after 30 min) of ZDR (Figure 15n), KDP (Figure 15o) and ρhv (Figure 15p), which is of great significance for quantitative applications of PRD (precipitation estimation, hydrometeor classification, and so on). By contrast, ExpVrZhKdp is the second choice, having a suboptimal effect comprehensively. It is noteworthy that the RMSEs of ExpVrZhRhv nearly coincide with that of ExpVrZh for all PRD variables, further indicating the minor positive impact of ρhv assimilation on the forecasts.

6. Summary and Conclusions

In this paper, based on the new simplified and parameterized polarimetric radar observation operators developed by Z21, the corresponding TL/AD operators necessarily required by the variational DA procedure are constructed. The gradient check of the cost function in the minimization step of analysis cycles proves the operators work reasonably. In fact, as mentioned previously, the updating of number concentrations which affect the DA effect [18] is not considered here and is worthy of further research.
To evaluate the influences of assimilating the differential reflectivity ZDR, specific differential phase KDP, cross-correlation coefficient ρhv, as well as horizontal reflectivity ZH and radial velocity Vr, cycled 3DVAR DA and forecast schemes are performed. The truth is generated with a 3 h simulation of an idealized supercell storm using the WRF model. The PRD is constructed from the idealized truth simulations via Z21 forward operators. Then, these observations are assimilated into the WRF model at a 5 min interval with the initial DA beginning at 30 min of model integration.
Using the experiment of assimilating both ZH and Vr as the benchmark run, OSSEs with additional assimilation of ZDR, KDP, and ρhv, respectively, are performed to examine the impact of each polarimetric variable on DA results. Additionally, an extra OSSE which assimilates all PRD is conducted to see whether a more balanced result can be achieved. The results suggest that the ZDR assimilation adjusts all model state variables to varying degrees, and helps reduce RMSEs at almost each DA cycle. Moreover, the assimilation process makes some polarimetric signatures, such as hook echo, ZDR column and melting layer, more obvious, and improves the distributions of hydrometeors and observations much closer to the truth as well. High sensitivity to liquid water of KDP is believed to be responsible for the more improved analysis of rainwater, and the intensity and pattern of PRD are also enhanced with additional assimilation of KDP, especially for ZH and KDP itself. The small discrepancy between OSSEs with and without ρhv assimilation indicates its limited ability in adjusting model state variables and improving the analysis accuracy. In contrast, due to the accumulation of favorable DA effects of different polarimetric variables, assimilating all PRD obtains more balanced analysis results, as expected. The above polarimetric signatures are also well presented in observation space, and many other improvements have been made, including the optimization of rain water analysis which is perfectly matched with the reference truth, and better descriptions of some small value areas in mixing ratios of cloud water, snow, and graupel. These different assimilation results are attributed to the different information contents and error effects from PRD variables. It is noted that these characteristics of assimilation PRD can change depending on the type of storm and model microphysics as well as observation errors, which needs further study. A series of prediction experiments, which initiate at the 120 min of different OSSEs and run over a 1 h period at a 15 min interval, are conducted to further study the impact of improved initial conditions on ensuing forecasts. The RMSEs of forecasts are calculated and the results show that the ZDR assimilation can bring obvious benefits for forecasts of ice hydrometeor species including cloud ice, hail, and graupel and even for the vertical wind field and perturbation pressure in the late stage of prediction. Assimilating KDP is found to improve the prediction of water vapor, cloud water, and rainwater, which is consistent with its performance in the liquid water analysis. There is no visible promoting effect on the prediction of most variables, suggesting that ρhv assimilation is not very helpful for the improvement of NWP. Similar to analysis results, the OSSE assimilating all PRD shows a more comprehensive performance in terms of RMSEs, having the most accurate (smallest RMSE) results of horizontal wind field, perturbation potential temperature, cloud water, rain water, snow, hail, and graupel in the end of forecast period. In addition, the forecasts in observation space are also better, with higher accuracy in the second half of period.
Although these findings illustrate the potential usefulness of PRD assimilation in improving convective-scale NWP, additional important research questions remain. A subsequent study will apply the PRD assimilation to real data cases to assess its benefits to real-world forecasts. Furthermore, the feasibility and capability of this new parameterized polarimetric radar operators under different MP schemes require investigation as well.

Author Contributions

Conceptualization, M.D. and J.G.; methodology, J.G., G.Z. and Y.W.; validation, M.D., J.G. and G.Z.; formal analysis, M.D. and J.G.; investigation, M.D.; resources, J.G., P.L.H. and C.C.; data curation, M.D.; writing—original draft preparation, M.D.; writing—review and editing, J.G., G.Z. and P.L.H.; visualization, M.D.; supervision, J.G., P.L.H. and C.C.; project administration, J.G.; funding acquisition, M.D. and C.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was jointly supported by the Key Program for International S&T Cooperation Projects of China, grant number 2016YFE0109400; the National Natural Science Foundation of China (NNSFC), grant number 42005121; and the Science and Technology Development Foundation of Hubei Meteorological Bureau, grant number 2019Y03.

Data Availability Statement

Not available.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

All abbreviations and acronyms used in this paper are summarized briefly in Table A1.
Table A1. Explanation of abbreviations and acronyms used in the paper.
Table A1. Explanation of abbreviations and acronyms used in the paper.
Abbreviation/AcronymExplanation
PRDPolarimetric radar data
NWPNumerical weather prediction
DAData assimilation
TLTangent linear model
ADAdjoint model
3DVARThree-dimensional variational system/method
4DVARFour-dimensional variational system/method
EnKFEnsemble Kalman filter
EnSRFEnsemble square-root Kalman filter
VrRadial velocity
ZHHorizontal reflectivity
ZDRDifferential reflectivity
ϕHVDifferential phase
KDPSpecific differential phase
ρhvCross-correlation coefficient
QPEQuantitative precipitation estimation
HCHydrometeor classification
MPMicrophysical parameterization scheme of model
SMSingle-moment scheme
DMDouble-moment scheme
WRFWeather research and forecasting model
WRF-ARWAdvanced research weather research and forecasting model
MCS(s)Mesoscale convective system(s)
OSSE(s)Observing system simulation experiment(s)
DSDDrop size distribution
PSDParticle size distribution
q Mixing ratio
N t Number concentration
W Water content
ρ a Air density
ρ w Water density
D m Mass-weighted diameter
γ x Percentage of melting
APRSAdvanced Regional Prediction System
CAPSCenter for Analysis and Prediction of Storms
NSSLNational Severe Storms Laboratory
J08Represents the article of Jung et al., (2008)
J10Represents the article of Jung et al., (2010)
Z21Represents the article of Zhang et al., (2021)
x b Background vector in the cost function
x Analysis vector in the cost function
y o Observation vector in the cost function
B Background error covariance matrix of model
R Observation error covariance matrix of model
H x Forward operator
J c x Constraints in the cost function
RUCRapid Update Cycle
AGLAbove ground level
u Horizontal wind in u-direction
v Horizontal wind in v-direction
w Vertical velocity
θ Perturbation potential temperature
q c Mixing ratio of cloud water
q i Mixing ratio of cloud ice
q r Mixing ratio of rain water
q s Mixing ratio of snow
q h Mixing ratio of hail
q g Mixing ratio of graupel
N c Number concentration of cloud water
N i Number concentration of cloud ice
N r Number concentration of rain water
N s Number concentration of snow
N h Number concentration of hail
N g Number concentration of graupel
VCPVolume coverage pattern
WSR-88DWeather surveillance radar—1988 Doppler
RMSE(s)Root mean square error(s)
ExpVrZhExperiment on assimilation of Vr and ZH
ExpVrZhZdrExperiment on assimilation of Vr, ZH, and ZDR
ExpVrZhKdpExperiment on assimilation of Vr, ZH, and KDP
ExpVrZhRhvExperiment on assimilation of Vr, ZH, and ρhv
ExpVrZhPolExperiment on assimilation of Vr, ZH, ZDR, KDP, and ρhv
CPAConvective precipitation area
SPAStratiform precipitation area
LWCLiquid water content

References

  1. Sun, J.Z.; Xue, M.; Wilson, J.W.; Zawadzki, I.; Ballard, S.P.; Onvlee-Hooimeyer, J.; Joe, P.; Barker, D.M.; Li, P.-W.; Golding, B.; et al. Use of NWP for nowcasting convective precipitation: Recent progress and challenges. Bull. Am. Meteorol. Soc. 2014, 95, 409–426. [Google Scholar] [CrossRef] [Green Version]
  2. Gao, J.D.; Xue, M.; Shapiro, A.; Droegemeier K., K. A variational analysis for the retrieval of three-dimensional mesoscale wind fields from two Doppler radars. Mon. Weather Rev. 1999, 127, 2128–2142. [Google Scholar] [CrossRef]
  3. Gao, J.D.; Xue, M.; Brewster, K.; Droegemeier, K.K. A three-dimensional data analysis method with recursive filter for Doppler radars. J. Atmos. Ocean. Technol. 2004, 21, 457–469. [Google Scholar] [CrossRef]
  4. Gao, J.D.; Stensrud, D.J. Assimilation of reflectivity data in a convective-scale, cycled 3DVAR framework with hydrometeor classification. J. Atmos. Sci. 2012, 69, 1054–1065. [Google Scholar] [CrossRef]
  5. Hu, M.; Xue, M.; Brewster, K. 3DVAR and cloud analysis with WSR-88D level-II data for the prediction of Fort Worth tornadic thunderstorms. Part I: Cloud analysis and its impact. Mon. Weather Rev. 2006, 134, 675–698. [Google Scholar] [CrossRef]
  6. Xie, Y.; Koch, S.E.; McGinley, J.A.; Albers, S.; Bieringer, P.; Wolfson, M.; Chan, M. A space and time multiscale analysis system: A sequential variational analysis approach. Mon. Weather Rev. 2011, 139, 1224–1240. [Google Scholar] [CrossRef]
  7. Lai, A.W.; Gao, J.D.; Koch, S.E.; Wang, Y.H.; Pan, S.J.; Fierro, A.O.; Cui, C.G.; Min, J.Z. Assimilation of radar radial velocity, reflectivity, and pseudo–water vapor for convective-scale NWP in a variational framework. Mon. Weather Rev. 2019, 147, 2877–2900. [Google Scholar] [CrossRef]
  8. Sun, J.Z.; Crook, N.A. Dynamical and microphysical retrieval from Doppler radar observations using a cloud model and its adjoint. Part I: Model development and simulated data experiments. J. Atmos. Sci. 1997, 54, 1642–1661. [Google Scholar] [CrossRef]
  9. Sun, J.Z.; Crook, N.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 storm. J. Atmos. Sci. 1998, 55, 835–852. [Google Scholar] [CrossRef]
  10. Sun, J.Z. Convective-scale assimilation of radar data: Progress and challenges. Q. J. R. Meteor. Soc. 2005, 131, 3439–3463. [Google Scholar] [CrossRef]
  11. Wang, H.L.; Sun, J.Z.; Zhang, X.; Huang, X.-Y.; Auligne, T. Radar data assimilation with WRF 4D-Var. Part I: System development and preliminary testing. Mon. Weather Rev. 2013, 141, 2224. [Google Scholar] [CrossRef]
  12. Sun, J.Z.; Wang, H.L. Radar data assimilation with WRF 4D-Var. Part II: Comparison with 3D-Var for a squall line over the U.S. Great Plains. Mon. Weather Rev. 2013, 141, 2245–2264. [Google Scholar] [CrossRef]
  13. Zhang, F.; Snyder, C.; Sun, J.Z. Impacts of initial estimate and observations on the convective-scale data assimilation with an ensemble Kalman filter. Mon. Weather Rev. 2004, 132, 1238–1253. [Google Scholar] [CrossRef]
  14. Tong, M.; Xue, M. Ensemble Kalman filter assimilation of Doppler radar data with a compressible nonhydrostatic model: OSS experiments. Mon. Weather Rev. 2005, 133, 1789–1807. [Google Scholar] [CrossRef] [Green Version]
  15. Thompson, T.E.; Wicker, L.; Wang, X.; Potvin, C. A comparison between the Local Ensemble Transform Kalman Filter and the Ensemble Square Root Filter for the assimilation of radar data in convective-scale models. Q. J. R. Meteor. Soc. 2014, 141, 1163–1176. [Google Scholar] [CrossRef]
  16. Putnam, B.J.; Xue, M.; Jung, Y.; Snook, N.; Zhang, G.F. Ensemble probabilistic prediction of a mesoscale convective system and associated polarimetric radar variables using single-moment and double-moment microphysics schemes and EnKF radar data assimilation. Mon. Weather Rev. 2017, 145, 2257–2279. [Google Scholar] [CrossRef]
  17. Putnam, B.J.; Xue, M.; Jung, Y.; Snook, N.; Zhang, G.F. Ensemble Kalman filter assimilation of polarimetric radar observations for the 20 May 2013 Oklahoma tornadic supercell case. Mon. Weather Rev. 2019, 147, 2511–2533. [Google Scholar] [CrossRef]
  18. Zhu, K.F.; Xue, M.; Ouyang, K.; Jung, Y. Assimilating polarimetric radar data with an ensemble Kalman filter: OSSEs with a tornadic supercell storm simulated with a double-moment microphysics scheme. Q. J. R. Meteor. Soc. 2020, 146, 1880–1900. [Google Scholar] [CrossRef]
  19. Kumjian, M.R. 2013: Principles and applications of dualpolarization weather radar. Part I: Description of the polarimetric radar variables. J. Oper. Meteor. 2013, 1, 226–242. [Google Scholar] [CrossRef]
  20. Zhao, K.; Huang, H.; Wang, M.J.; Lee, W.-C.; Chen, G.; Wen, L.; Wen, J.; Zhang, G.F.; Xue, M.; Yang, Z.W.; et al. Recent progress in dual-polarization radar research and applications in China. Adv. Atmos. Sci. 2019, 36, 961–974. [Google Scholar] [CrossRef]
  21. Sachidananda, M.; Zrnic, D.S. Rain rate estimates from differential polarimetric measurements. J. Atmos. Ocean. Technol. 1987, 4, 588–598. [Google Scholar] [CrossRef] [Green Version]
  22. Gorgucci, E.; Scarchilli, G.; Chandrasekar, V.; Bringi, V.N. Rainfall estimation from polarimetric radar measurements: Composite algorithms immune to variability in raindrop shape–size relation. J. Atmos. Ocean. Technol. 2001, 18, 1773–1786. [Google Scholar] [CrossRef]
  23. Ryzhkov, A.V.; Giangrande, S.E.; Schuur, T.J. Rainfall estimation with a polarimetric prototype of WSR-88D. J. Appl. Meteor. 2005, 44, 502–515. [Google Scholar] [CrossRef] [Green Version]
  24. Huang, H.; Zhao, K.; Zhang, G.F.; Lin, Q.; Wen, L.; Chen, G.; Yang, Z.W.; Wang, M.J.; Hu, D.M. Quantitative precipitation estimation with operational polarimetric radar measurements in southern China: A differential phase-based variational approach. J. Atmos. Ocean. Technol. 2018, 35, 1253–1271. [Google Scholar] [CrossRef]
  25. Straka, J.M.; Zrnic, D.S.; Ryzhkov, A.V. Bulk hydrometeor classification and quantification using polarimetric radar data: Synthesis of relations. J. Appl. Meteor. 2000, 39, 1341–1372. [Google Scholar] [CrossRef]
  26. Park, H.; Ryzhkov, A.V.; Zrnic, D.S.; Kim, K. The hydrometeor classification algorithm for the polarimetric WSR-88D: Description and application to an MCS. Weather Forecast. 2009, 24, 730–748. [Google Scholar] [CrossRef]
  27. Yang, J.; Zhao, K.; Zhang, G.F.; Chen, G.; Huang, H.; Chen, H.N. A Bayesian hydrometeor classification algorithm for C-band polarimetric radar. Remote Sens. 2019, 11, 1884. [Google Scholar] [CrossRef] [Green Version]
  28. Vivekanandan, J.; Ellis, S.M.; Oye, R.; Zrnic, D.S.; Ryzhkov, A.V.; Straka, J. Cloud microphysics retrieval using S-band dual-polarization radar measurements. Bull. Am. Meteor. Soc. 1999, 80, 381–388. [Google Scholar] [CrossRef]
  29. Brandes, E.A.; Zhang, G.F.; Vivekanandan, J. Comparison of polarimetric radar drop size distribution retrieval algorithms. J. Atmos. Ocean. Technol. 2004, 21, 584–598. [Google Scholar] [CrossRef] [Green Version]
  30. Zhang, G.F.; Vivekanandan, J.; Brandes, E.A. A method for estimating rain rate and drop size distribution from polarimetric radar measurements. IEEE Trans. Geosci. Remote Sens. 2001, 39, 830–841. [Google Scholar] [CrossRef] [Green Version]
  31. Cao, Q.; Zhang, G.F.; Xue, M. A variational approach for retrieving raindrop size distribution from polarimetric radar measurements in the presence of attenuation. J. Appl. Meteor. Climatol. 2013, 52, 169–185. [Google Scholar] [CrossRef] [Green Version]
  32. Kenney, P.C.; Detwiler, A.G. A case study of the origin of hail in a multicell thunderstorm using in situ aircraft and polarimetric radar data. J. Appl. Meteor. Climatol. 2003, 42, 1679–1690. [Google Scholar] [CrossRef] [Green Version]
  33. Kumjian, M.R.; Ryzhkov, A.V. Polarimetric signatures in supercell thunderstorms. J. Appl. Meteor. Climatol. 2008, 47, 1940–1961. [Google Scholar] [CrossRef]
  34. Kumjian, M.R.; Ryzhkov, A.V.; Melnikov, V.M.; Schuur, T.J. Rapid-scan super-resolution observations of a cyclic supercell with a dual-polarimetric WSR-88D. Mon. Weather Rev. 2010, 138, 3762–3786. [Google Scholar] [CrossRef]
  35. Van Den Broeke, M.S. Polarimetric radar metrics related to tornado life cycles and intensity in supercell storms. Mon. Weather Rev. 2017, 145, 3671–3686. [Google Scholar] [CrossRef]
  36. Vivekanandan, J.; Bringi, V.N.; Hagen, M.; Meischner, P. Polarimetric radar studies of atmospheric ice particles. IEEE Trans. Geosci. Remote Sens. 1994, 32, 1–10. [Google Scholar] [CrossRef] [Green Version]
  37. Ryzhkov, A.V.; Zrnic, D.S.; Gordon, B.A. Polarimetric method for ice water content determination. J. Appl. Meteor. 1998, 37, 125–134. [Google Scholar] [CrossRef]
  38. Pfeifer, M.; Craig, G.C.; Hagen, M.; Keil, C. A polarimetric radar forward operator for model evaluation. J. Appl. Meteor. Climatol. 2008, 47, 3202–3220. [Google Scholar] [CrossRef] [Green Version]
  39. Jung, Y.; Zhang, G.F.; Xue, M. Assimilation of simulated polarimetric radar data for a convective storm using the ensemble Kalman filter. Part I: Observation operators for reflectivity and polarimetric variables. Mon. Weather Rev. 2008, 136, 2228–2245. [Google Scholar] [CrossRef]
  40. Jung, Y.; Xue, M.; Zhang, G.F.; Straka, J.M. Assimilation of simulated polarimetric radar data for convective storm using the ensemble Kalman filter. Part II: Impact of polarimetric data on storm analysis. Mon. Weather Rev. 2008, 136, 2246–2260. [Google Scholar] [CrossRef]
  41. Jung, Y.; Xue, M.; Zhang, G.F. Simultaneous estimation of microphysical parameters and the atmospheric state using simulated polarimetric radar data and an ensemble Kalman filter in the presence of an observation operator error. Mon. Weather Rev. 2010, 138, 539–562. [Google Scholar] [CrossRef] [Green Version]
  42. Jung, Y.; Xue, M.; Zhang, G.F. Simulations of polarimetric radar signatures of a supercell storm using a two-moment bulk microphysics scheme. J. Appl. Meteor. Climatol. 2010, 49, 146–163. [Google Scholar] [CrossRef]
  43. Li, X.L.; Mecikalski, J.R. Assimilation of the dual-polarization Doppler radar data for a convective storm with a warm-rain radar forward operator. J. Geophys. Res. Atmos. 2010, 115, D16208. [Google Scholar] [CrossRef] [Green Version]
  44. Li, X.L.; Mecikalski, J.R. Impact of the dual-polarization Doppler radar data on two convective storms with a warm-rain radar forward operator. Mon. Weather Rev. 2012, 140, 2147–2167. [Google Scholar] [CrossRef]
  45. Ulbrich, C.W.; Atlas, D. Assessment of the contribution of different polarization to improved rainfall measurements. Radio Sci. 1984, 19, 49–57. [Google Scholar] [CrossRef]
  46. Bringi, V.N.; Chandrasekar, V. Polarimetric Doppler Weather Radar: Principles and Applications; Cambridge University Press: Cambridge, UK, 2001; pp. 1–636. [Google Scholar]
  47. Li, X.L.; Mecikalski, J.R.; Posselt, D. An ice-phase microphysics forward model and preliminary results of polarimetric radar data assimilation. Mon. Weather Rev. 2017, 145, 683–708. [Google Scholar] [CrossRef]
  48. Wang, S.Z.; Liu, Z.Q. A radar reflectivity operator with ice-phase hydrometeors for variational data assimilation (RadZIceVar v1.0) and its evaluation with real radar data. Geosci. Model Dev. 2019, 12, 4031–4051. [Google Scholar] [CrossRef] [Green Version]
  49. Xue, M.; Droegemeier, K.K.; Wong, V.; Shapiro, A.; Brewster, K.; Carr, F.; Weber, D.; Liu, Y.; Wang, D. The Advanced Regional Prediction System (ARPS)—A multi-scale nonhydrostatic atmospheric simulation and prediction tool. Part II: Model physics and applications. Meteor. Atmos. Phys. 2001, 76, 143–165. [Google Scholar] [CrossRef]
  50. Xue, M.; Wang, D.H.; Gao, J.D.; Brewster, K.; Droegemeier, K.K. The Advanced Regional Prediction System (ARPS), storm-scale numerical weather prediction and data assimilation. Meteor. Atmos. Phys. 2003, 82, 139–170. [Google Scholar] [CrossRef]
  51. Carlin, J.T.; Gao, J.D.; Snyder, J.C.; Ryzhkov, A.V. Assimilation of ZDR columns for improving the spinup and forecast of convective storms in storm-scale models: Proof-of-concept experiments. Mon. Weather Rev. 2017, 145, 5033–5057. [Google Scholar] [CrossRef]
  52. Kumjian, M.R.; Khain, A.P.; Benmoshe, N.; Ilotoviz, E.; Ryzhkov, A.V.; Phillips, V.T. The anatomy and physics of ZDR columns: Investigating a polarimetric radar signature with a spectral bin microphysical model. J. Appl. Meteor. Climatol. 2014, 53, 1820–1843. [Google Scholar] [CrossRef]
  53. Zhang, G.F.; Gao, J.D.; Du, M.Y. Parameterized forward operators for simulation and assimilation of polarimetric radar data with numerical weather predictions. Adv. Atmos. Sci. 2021, 38, 737–754. [Google Scholar] [CrossRef]
  54. Lin, Y.L.; Farley, R.D.; Orville, H.D. Bulk parameterization of the snow field in a cloud model. J. Climate Appl. Meteor. 1983, 22, 1065–1092. [Google Scholar] [CrossRef] [Green Version]
  55. Milbrandt, J.A.; Yau, M.K. A multimoment bulk microphysics parameterization. Part I: Analysis of the role of the spectral shape parameter. J. Atmos. Sci. 2005, 62, 3051–3064. [Google Scholar] [CrossRef] [Green Version]
  56. Morrison, H.; Pinto, J.O. Mesoscale modeling of springtime arctic mixed-phase clouds using a new two-moment bulk microphysics scheme. J. Atmos. Sci. 2005, 62, 3683–3704. [Google Scholar] [CrossRef]
  57. Thompson, G.; Field, P.R.; Rasmussen, R.M.; Hall, W.R. Explicit forecasts of winter precipitation using an improved bulk microphysics scheme. Part II: Implementation of a new snow parameterization. Mon. Weather Rev. 2008, 136, 5095–5115. [Google Scholar] [CrossRef]
  58. Mansell, E.R.; Ziegler, C.L.; Bruning, E.C. Simulated electrification of a small thunderstorm with two–moment bulk microphysics. J. Atmos. Sci. 2010, 67, 171–194. [Google Scholar] [CrossRef]
  59. Doviak, R.J.; Zrnic, D.S. Doppler Radar and Weather Observations, 2nd ed.; Academic Press: New York, NY, USA, 1993; pp. 1–562. [Google Scholar]
  60. Zhang, G.F. Weather Radar Polarimetry; CRC Press: Boca Raton, FL, USA, 2016; pp. 1–304. [Google Scholar]
  61. Mahale, V.N.; Zhang, G.F.; Xue, M.; Gao, J.D.; Reeves, H.D. Variational retrieval of rain microphysics and related parameters from polarimetric radar data with a parameterized operator. J. Atmos. Ocean. Technol. 2019, 36, 2483–2500. [Google Scholar] [CrossRef]
  62. Benjamin, S.G.; Devenyi, D.; Weygandt, S.S.; Brundage, K.J.; Brown, J.M.; Grell, G.A.; Kim, D.; Schwartz, B.E.; Smirnova, T.G.; Smith, T.L. An hourly assimilation-forecast cycle: The RUC. Mon. Weather Rev. 2004, 132, 495–518. [Google Scholar] [CrossRef]
  63. Evensen, G. The ensemble Kalman filter: Theoretical formulation and practical implementation. Ocean. Technol. 2003, 53, 343–367. [Google Scholar] [CrossRef]
  64. Errico, R.M. What is an adjoint model? Bull. Am. Meteor. Soc. 1997, 78, 2577–2591. [Google Scholar] [CrossRef]
  65. Giering, R.; Kaminski, T. Recipes for adjoint code construction. ACM Trans. Math. Software 1998, 24, 437–474. [Google Scholar] [CrossRef]
  66. Kalnay, E. Atmospheric Modeling, Data Assimilation, and Predictability; Cambridge University Press: Cambridge, UK, 2003; pp. 1–341. [Google Scholar]
  67. Ge, G.Q.; Gao, J.D.; Xue, M. Impacts of assimilating measurements of different state variables with a simulated supercell storm and three-dimensional variational method. Mon. Weather Rev. 2013, 141, 2759–2777. [Google Scholar] [CrossRef] [Green Version]
  68. Gao, J.D.; Stensrud, D.J. Some observing system simulation experiments with a hybrid 3DEnVAR system for storm-scale radar data assimilation. Mon. Weather Rev. 2014, 142, 3326–3346. [Google Scholar] [CrossRef]
  69. Gao, J.D.; Fu, C.H.; Stensrud, D.J.; Kain, J.S. OSSEs for an ensemble 3DVAR data assimilation system with radar observations of convective storms. J. Atmos. Sci. 2016, 73, 2403–2426. [Google Scholar] [CrossRef]
Figure 1. Illustration of DA and forecast cycles used for the idealized case study. Observations represent the assimilated PRD including radar radial velocity, reflectivity, and all polarimetric variables simulated from the truth simulation output of the WRF run. The observations are assimilated by the 3DVAR system based on the background field of previous 5 min forecast result.
Figure 1. Illustration of DA and forecast cycles used for the idealized case study. Observations represent the assimilated PRD including radar radial velocity, reflectivity, and all polarimetric variables simulated from the truth simulation output of the WRF run. The observations are assimilated by the 3DVAR system based on the background field of previous 5 min forecast result.
Remotesensing 13 03060 g001
Figure 2. The mean RMSEs of the forecasts and analyses throughout the 5 min 3DVAR DA cycles for ExpVrZh (black line), ExpVrZhZdr (red line), ExpVrZhKdp (blue line), ExpVrZhRhv (cyan line), and ExpVrZhPol (purple line), respectively, averaged over grid points where the observed (simulated) ZH is greater than 10 dBZ. (a) u , (b) v , (c) w , mixing ratios of (d) cloud water q c , (e) cloud ice q i , (f) rain water q r , (g) snow q s , (h) hail q h , and graupel q g , and the same calculating procedure as truth simulations for polarimetric radar variables (i) ZH, (j) ZDR, (k) KDP, and (l) ρhv via Z21 forward operators.
Figure 2. The mean RMSEs of the forecasts and analyses throughout the 5 min 3DVAR DA cycles for ExpVrZh (black line), ExpVrZhZdr (red line), ExpVrZhKdp (blue line), ExpVrZhRhv (cyan line), and ExpVrZhPol (purple line), respectively, averaged over grid points where the observed (simulated) ZH is greater than 10 dBZ. (a) u , (b) v , (c) w , mixing ratios of (d) cloud water q c , (e) cloud ice q i , (f) rain water q r , (g) snow q s , (h) hail q h , and graupel q g , and the same calculating procedure as truth simulations for polarimetric radar variables (i) ZH, (j) ZDR, (k) KDP, and (l) ρhv via Z21 forward operators.
Remotesensing 13 03060 g002
Figure 3. The reflectivity ZH (unit: dBZ) at 1 km altitude at 80 min for (a) truth, and analysis fields of (b) ExpVrZh, (c) ExpVrZhZdr, (d) ExpVrZhKdp, (e) ExpVrZhRhv, and (f) ExpVrZhPol, respectively. The two assumed radars are located at the bottom-left and top-right corners as shown in Figure 3a.
Figure 3. The reflectivity ZH (unit: dBZ) at 1 km altitude at 80 min for (a) truth, and analysis fields of (b) ExpVrZh, (c) ExpVrZhZdr, (d) ExpVrZhKdp, (e) ExpVrZhRhv, and (f) ExpVrZhPol, respectively. The two assumed radars are located at the bottom-left and top-right corners as shown in Figure 3a.
Remotesensing 13 03060 g003
Figure 4. The same as Figure 3, but for the differential reflectivity ZDR (unit: dB) of (a) truth, and analysis fields of (b) ExpVrZh, (c) ExpVrZhZdr, (d) ExpVrZhKdp, (e) ExpVrZhRhv, and (f) ExpVrZhPol, respectively.
Figure 4. The same as Figure 3, but for the differential reflectivity ZDR (unit: dB) of (a) truth, and analysis fields of (b) ExpVrZh, (c) ExpVrZhZdr, (d) ExpVrZhKdp, (e) ExpVrZhRhv, and (f) ExpVrZhPol, respectively.
Remotesensing 13 03060 g004
Figure 5. The same as Figure 3, but for the specific differential phase KDP (unit: °/km) of (a) truth, and analysis fields of (b) ExpVrZh, (c) ExpVrZhZdr, (d) ExpVrZhKdp, (e) ExpVrZhRhv, and (f) ExpVrZhPol, respectively.
Figure 5. The same as Figure 3, but for the specific differential phase KDP (unit: °/km) of (a) truth, and analysis fields of (b) ExpVrZh, (c) ExpVrZhZdr, (d) ExpVrZhKdp, (e) ExpVrZhRhv, and (f) ExpVrZhPol, respectively.
Remotesensing 13 03060 g005
Figure 6. The same as Figure 3, but for the cross-correlation coefficient ρhv (dimensionless) of (a) truth, and analysis fields of (b) ExpVrZh, (c) ExpVrZhZdr, (d) ExpVrZhKdp, (e) ExpVrZhRhv, and (f) ExpVrZhPol, respectively.
Figure 6. The same as Figure 3, but for the cross-correlation coefficient ρhv (dimensionless) of (a) truth, and analysis fields of (b) ExpVrZh, (c) ExpVrZhZdr, (d) ExpVrZhKdp, (e) ExpVrZhRhv, and (f) ExpVrZhPol, respectively.
Remotesensing 13 03060 g006
Figure 7. Vertical cross-sections of reflectivity ZH along the black dashed line in Figure 3a for (a) truth, and analysis fields of (b) ExpVrZh, (c) ExpVrZhZdr, (d) ExpVrZhKdp, (e) ExpVrZhRhv, and (f) ExpVrZhPol, respectively.
Figure 7. Vertical cross-sections of reflectivity ZH along the black dashed line in Figure 3a for (a) truth, and analysis fields of (b) ExpVrZh, (c) ExpVrZhZdr, (d) ExpVrZhKdp, (e) ExpVrZhRhv, and (f) ExpVrZhPol, respectively.
Remotesensing 13 03060 g007
Figure 8. The same as Figure 7, but for the differential reflectivity ZDR of (a) truth, and analysis fields of (b) ExpVrZh, (c) ExpVrZhZdr, (d) ExpVrZhKdp, (e) ExpVrZhRhv, and (f) ExpVrZhPol, respectively.
Figure 8. The same as Figure 7, but for the differential reflectivity ZDR of (a) truth, and analysis fields of (b) ExpVrZh, (c) ExpVrZhZdr, (d) ExpVrZhKdp, (e) ExpVrZhRhv, and (f) ExpVrZhPol, respectively.
Remotesensing 13 03060 g008
Figure 9. The same as Figure 7, but for the specific differential phase KDP of (a) truth, and analysis fields of (b) ExpVrZh, (c) ExpVrZhZdr, (d) ExpVrZhKdp, (e) ExpVrZhRhv, and (f) ExpVrZhPol, respectively.
Figure 9. The same as Figure 7, but for the specific differential phase KDP of (a) truth, and analysis fields of (b) ExpVrZh, (c) ExpVrZhZdr, (d) ExpVrZhKdp, (e) ExpVrZhRhv, and (f) ExpVrZhPol, respectively.
Remotesensing 13 03060 g009
Figure 10. The same as Figure 7, but for the cross-correlation coefficient ρhv of (a) truth, and analysis fields of (b) ExpVrZh, (c) ExpVrZhZdr, (d) ExpVrZhKdp, (e) ExpVrZhRhv, and (f) ExpVrZhPol, respectively.
Figure 10. The same as Figure 7, but for the cross-correlation coefficient ρhv of (a) truth, and analysis fields of (b) ExpVrZh, (c) ExpVrZhZdr, (d) ExpVrZhKdp, (e) ExpVrZhRhv, and (f) ExpVrZhPol, respectively.
Remotesensing 13 03060 g010
Figure 11. The same as Figure 7, but for the mixing ratio of cloud ice q i (unit: g/kg) of (a) truth, and analysis fields of (b) ExpVrZh, (c) ExpVrZhZdr, (d) ExpVrZhKdp, (e) ExpVrZhRhv, and (f) ExpVrZhPol, respectively.
Figure 11. The same as Figure 7, but for the mixing ratio of cloud ice q i (unit: g/kg) of (a) truth, and analysis fields of (b) ExpVrZh, (c) ExpVrZhZdr, (d) ExpVrZhKdp, (e) ExpVrZhRhv, and (f) ExpVrZhPol, respectively.
Remotesensing 13 03060 g011
Figure 12. The same as Figure 7, but for the mixing ratio of rain water q r (unit: g/kg) of (a) truth, and analysis fields of (b) ExpVrZh, (c) ExpVrZhZdr, (d) ExpVrZhKdp, (e) ExpVrZhRhv, and (f) ExpVrZhPol, respectively.
Figure 12. The same as Figure 7, but for the mixing ratio of rain water q r (unit: g/kg) of (a) truth, and analysis fields of (b) ExpVrZh, (c) ExpVrZhZdr, (d) ExpVrZhKdp, (e) ExpVrZhRhv, and (f) ExpVrZhPol, respectively.
Remotesensing 13 03060 g012
Figure 13. The same as Figure 7, but for the mixing ratio of hail q h (unit: g/kg) of (a) truth, and analysis fields of (b) ExpVrZh, (c) ExpVrZhZdr, (d) ExpVrZhKdp, (e) ExpVrZhRhv, and (f) ExpVrZhPol, respectively.
Figure 13. The same as Figure 7, but for the mixing ratio of hail q h (unit: g/kg) of (a) truth, and analysis fields of (b) ExpVrZh, (c) ExpVrZhZdr, (d) ExpVrZhKdp, (e) ExpVrZhRhv, and (f) ExpVrZhPol, respectively.
Remotesensing 13 03060 g013
Figure 14. The same as Figure 7, but for the mixing ratio of graupel q g (unit: g/kg) of (a) truth, and analysis fields of (b) ExpVrZh, (c) ExpVrZhZdr, (d) ExpVrZhKdp, (e) ExpVrZhRhv, and (f) ExpVrZhPol, respectively.
Figure 14. The same as Figure 7, but for the mixing ratio of graupel q g (unit: g/kg) of (a) truth, and analysis fields of (b) ExpVrZh, (c) ExpVrZhZdr, (d) ExpVrZhKdp, (e) ExpVrZhRhv, and (f) ExpVrZhPol, respectively.
Remotesensing 13 03060 g014
Figure 15. The mean RMSEs of 1 h forecasts at a 15 min intervals initiated at 120 min for ExpVrZh (black line), ExpVrZhZdr (red line), ExpVrZhKdp (blue line), ExpVrZhRhv (cyan line), and ExpVrZhPol (purple line), respectively, averaged over grid points with observed (simulated) ZH exceeding 10 dBZ. (a) u , (b) v , (c) w , (d) perturbation potential temperature θ , (e) perturbation pressure p , mixing ratios of (f) water vapor q v , (g) cloud water q c , (h) cloud ice q i , (i) rain water q r , (j) snow q s , (k) hail q h , and (l) graupel q g , and simulated PRD in observation space: (m) ZH, (n) ZDR, (o) KDP, and (p) ρhv.
Figure 15. The mean RMSEs of 1 h forecasts at a 15 min intervals initiated at 120 min for ExpVrZh (black line), ExpVrZhZdr (red line), ExpVrZhKdp (blue line), ExpVrZhRhv (cyan line), and ExpVrZhPol (purple line), respectively, averaged over grid points with observed (simulated) ZH exceeding 10 dBZ. (a) u , (b) v , (c) w , (d) perturbation potential temperature θ , (e) perturbation pressure p , mixing ratios of (f) water vapor q v , (g) cloud water q c , (h) cloud ice q i , (i) rain water q r , (j) snow q s , (k) hail q h , and (l) graupel q g , and simulated PRD in observation space: (m) ZH, (n) ZDR, (o) KDP, and (p) ρhv.
Remotesensing 13 03060 g015
Table 1. Results of test for Z21 TL/AD operators
Table 1. Results of test for Z21 TL/AD operators
αR_ZHR_ZDRR_KDPR_ρhv
1E00.6990451650731240.4830749819578291.025288191328010.490697018550403
1E-10.9548866430476960.8845977844118370.9952784603349510.931147252389009
1E-20.9952278049580500.9876140134760370.9995213338035390.992331436494260
1E-30.9995199839647570.9987521968253160.9999520658284010.999224559666156
1E-40.9999519702304810.9998751269655560.9999952059042660.999922370046164
1E-50.9999951969133380.9999875131400140.9999995205778620.999992241080020
1E-60.9999995202215900.9999987724737370.9999999520193590.999999175462863
1E-70.9999999597772961.000000060871910.9999999976009950.999999828029847
1E-81.000000010300941.000002923085631.000000023462921.00001425319475
1E-90.9999993366523491.000023505296651.000000314409530.999973038437884
1E-101.000014493745671.000038178161930.9999996678615020.999114397669863
1E-111.000199747108501.000468314212641.000015831562150.999457853977071
1E-121.000368159256520.9951418631173891.000048158963440.996023290904985
1E-131.010472888137901.089409997595691.002149440047410.686912614417231
1E-141.010472888137901.386887266311851.018313140693343.43456307208615
Table 2. List of experiments.
Table 2. List of experiments.
ExperimentsObservationsDescription
ExpVrZhVr+ZHVr and ZH assimilated
ExpVrZhZdrVr+ZH+ZDRAs ExpVrZh with additional ZDR assimilated
ExpVrZhKdpVr+ZH+KDPAs ExpVrZh with additional KDP assimilated
ExpVrZhRhvVr+ZH+ρhvAs ExpVrZh with additional ρhv assimilated
ExpVrZhPolVr+ZH+ZDR+KDP+ρhvAs ExpVrZh with all PRD assimilated
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Du, M.; Gao, J.; Zhang, G.; Wang, Y.; Heiselman, P.L.; Cui, C. Assimilation of Polarimetric Radar Data in Simulation of a Supercell Storm with a Variational Approach and the WRF Model. Remote Sens. 2021, 13, 3060. https://doi.org/10.3390/rs13163060

AMA Style

Du M, Gao J, Zhang G, Wang Y, Heiselman PL, Cui C. Assimilation of Polarimetric Radar Data in Simulation of a Supercell Storm with a Variational Approach and the WRF Model. Remote Sensing. 2021; 13(16):3060. https://doi.org/10.3390/rs13163060

Chicago/Turabian Style

Du, Muyun, Jidong Gao, Guifu Zhang, Yunheng Wang, Pamela L. Heiselman, and Chunguang Cui. 2021. "Assimilation of Polarimetric Radar Data in Simulation of a Supercell Storm with a Variational Approach and the WRF Model" Remote Sensing 13, no. 16: 3060. https://doi.org/10.3390/rs13163060

APA Style

Du, M., Gao, J., Zhang, G., Wang, Y., Heiselman, P. L., & Cui, C. (2021). Assimilation of Polarimetric Radar Data in Simulation of a Supercell Storm with a Variational Approach and the WRF Model. Remote Sensing, 13(16), 3060. https://doi.org/10.3390/rs13163060

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop