- About
- Editorial board
- Articles
- Special issues
- Highlight articles
- Manuscript tracking
- Subscribe to alerts
- Peer-review process
- For authors
- For editors and referees
- EGU publications
- Imprint
- Data protection

Journal cover
Journal topic
**Geoscientific Model Development**
An interactive open-access journal of the European Geosciences Union

Journal topic

- About
- Editorial board
- Articles
- Special issues
- Highlight articles
- Manuscript tracking
- Subscribe to alerts
- Peer-review process
- For authors
- For editors and referees
- EGU publications
- Imprint
- Data protection

- About
- Editorial board
- Articles
- Special issues
- Highlight articles
- Manuscript tracking
- Subscribe to alerts
- Peer-review process
- For authors
- For editors and referees
- EGU publications
- Imprint
- Data protection

**Model description paper**
22 Jun 2018

**Model description paper** | 22 Jun 2018

Observational operators for dual polarimetric radars in variational data assimilation systems (PolRad VAR v1.0)

^{1}Meteorological Research Institute, Japan Meteorological Agency, Tsukuba, 3050052, Japan^{2}Institute of Physics and Meteorology, University of Hohenheim, 70599 Stuttgart, Germany^{3}Japan Meteorological Agency, Tokyo, 1008122, Japan

^{1}Meteorological Research Institute, Japan Meteorological Agency, Tsukuba, 3050052, Japan^{2}Institute of Physics and Meteorology, University of Hohenheim, 70599 Stuttgart, Germany^{3}Japan Meteorological Agency, Tokyo, 1008122, Japan

**Correspondence**: Takuya Kawabata (tkawabat@mri-jma.go.jp)

**Correspondence**: Takuya Kawabata (tkawabat@mri-jma.go.jp)

Abstract

Back to toptop
We implemented two observational operators for dual polarimetric radars in two variational data assimilation systems: WRF Var, the Weather Research and Forecasting Model variational data assimilation system, and NHM-4DVAR, the nonhydrostatic variational data assimilation system for the Japan Meteorological Agency nonhydrostatic model. The operators consist of a space interpolator, two types of variable converters, and their linearized and transposed (adjoint) operators. The space interpolator takes account of the effects of radar-beam broadening in both the vertical and horizontal directions and climatological beam bending. The first variable converter emulates polarimetric parameters with model prognostic variables and includes attenuation effects, and the second one derives rainwater content from the observed polarimetric parameter (specific differential phase). We developed linearized and adjoint operators for the space interpolator and variable converters and then assessed whether the linearity of the linearized operators and the accuracy of the adjoint operators were good enough for implementation in variational systems. The results of a simple assimilation experiment showed good agreement between assimilation results and observations with respect to reflectivity and specific differential phase but not with respect to differential reflectivity.

Download & links

How to cite

Back to top
top
How to cite.

Kawabata, T., Schwitalla, T., Adachi, A., Bauer, H.-S., Wulfmeyer, V., Nagumo, N., and Yamauchi, H.: Observational operators for dual polarimetric radars in variational data assimilation systems (PolRad VAR v1.0), Geosci. Model Dev., 11, 2493–2501, https://doi.org/10.5194/gmd-11-2493-2018, 2018.

1 Introduction

Back to toptop
The Weather Research and Forecasting Model (WRF; Skamarock et al., 2008) is a widely used numerical weather model that was developed as a community model, and WRF Var, its data assimilation system (Barker et al., 2012), provides initial conditions for the model. NHM-4DVAR is a nonhydrostatic 4D-Var system for the Japan Meteorological Agency nonhydrostatic model (JMANHM; Saito, 2012) that functions at storm scale (Kawabata et al., 2007, 2014a). Many remote sensing data are available for NHM-4DVAR, such as the following: slant total delay, zenith total delay, and precipitable water vapor observed by global navigation satellite systems (GNSS; Kawabata et al., 2013); conventional radar data, including directly assimilated reflectivity data (Kawabata et al., 2011); and Doppler lidar data (Kawabata et al., 2014b). Because data assimilation associates observations with model fields, to make use of advanced observations, data assimilation methods need to be continuously developed and implemented into variational data assimilation systems.

Observations obtained by dual polarimetric radars are utilized by
operational systems at many meteorological and hydrological operation
centers (e.g., in the United States, France, Germany, and Japan) to improve
the accuracy of quantitative precipitation estimation (QPE). These radars
provide polarimetric parameters, including the horizontally polarized
reflectivity factor (*Z*_{H}), the vertically polarized reflectivity factor
(*Z*_{V}), differential reflectivity (*Z*_{DR}), and the specific differential
phase (*K*_{DP}). Many QPE methods that use these parameters have been
proposed (e.g., Jameson, 1991; Jameson and Caylor, 1994; Ryzhkov and
Zrnić, 1995; Anagnostou et al., 2008; Kim et al., 2010; Ryzhkov et al.,
2014; Adachi et al., 2015). Because QPE methods using dual polarimetric
radar parameters are expected to be better than methods using single
polarization radar data, we developed assimilation methods for dual
polarimetric radar observations for both WRF Var and NHM-4DVAR. The
objective of our study was thus to improve QPE, which was discussed in Bauer
et al. (2015) in the context of data assimilation with high resolution and a
rapid update cycle, and quantitative rainfall forecasts (QPF) through the
use of better analysis fields obtained by the assimilation of dual
polarimetric radar observations.

We chose an emulator (Zhang et al., 2001) and an estimator (Bringi and Chandrasekar, 2001) to use as forward operators after evaluating their accuracy (Kawabata et al., 2018). In addition, because both WRF Var and NHM-4DVAR consider only perturbations to rainwater in their tangent and adjoint models, our operators also deal only with rainwater and exclude ice particles. Although both the emulator of Jung et al. (2008a, b) and the estimator of Yokota et al. (2016) have been used previously as observational operators in ensemble Kalman filter data assimilation systems, to our knowledge, our study is their first implementation in variational assimilation systems. We refer to the current version of the operators as PolRad VAR v1.0.

The first author has mainly contributed to the WRF Var version of these operators developed over the rapid-update WRF 3D-Var system at the University of Hohenheim, Germany (see, e.g., Schwitalla et al., 2011; Schwitalla and Wulfmeyer, 2014; Bauer et al., 2015) and then to the version for NHM-4DVAR at the Meteorological Research Institute, Japan Meteorological Agency.

The scope of this paper is to provide the technical information on the observational operators and some evaluation results to help users understand the theoretical and practical aspects of the operators. The forward operators (space interpolator and variable converters) and their linearized (tangent linear) and transposed (adjoint) operators are described in Sect. 2. Section 3 describes the setup options of the observational operators, Sect. 4 presents verification and assimilation test results, and Sect. 5 is a summary.

2 Observational operators

Back to toptop
In variational data assimilation systems, a cost function is defined and then iteratively minimized until its gradient becomes zero. The cost function and its gradient are defined as

$$\begin{array}{ll}{\displaystyle}& {\displaystyle}J\left(\mathbf{x}\right)={\displaystyle \frac{\mathrm{1}}{\mathrm{2}}}{\left(\mathbf{x}-{\mathbf{x}}^{b}\right)}^{\mathrm{T}}{\mathbf{B}}^{-\mathrm{1}}\left(\mathbf{x}-{\mathbf{x}}^{b}\right)\\ \text{(1)}& {\displaystyle}& {\displaystyle}\phantom{\rule{1em}{0ex}}+{\displaystyle \frac{\mathrm{1}}{\mathrm{2}}}{\left(H\left(\mathbf{x}\right)-\mathbf{y}\right)}^{\mathrm{T}}{\mathbf{R}}^{-\mathrm{1}}\left(H\left(\mathbf{x}\right)-\mathbf{y}\right),\text{(2)}& {\displaystyle}& {\displaystyle}\mathrm{\nabla}J\left(\mathbf{x}\right)={\mathbf{B}}^{-\mathrm{1}}\left(\mathbf{x}-{\mathbf{x}}^{b}\right)+{\mathbf{H}}^{\mathrm{T}}{\mathbf{R}}^{-\mathrm{1}}\left(H\left(\mathbf{x}\right)-\mathbf{y}\right),\end{array}$$

where T denotes the transpose of a matrix; **x**, **x**^{b},
and **y** are model fields, first-guess fields, and observations,
respectively, and *H*(**x**), **H**, and **H**^{T}
represent the observational operators, their linearized operators (tangent
linear), and their transposed (adjoint) operators, respectively. The
observational operators work as variable converters (*H*_{v}) from
model fields **x** to observational values related to observations
**y** and as space interpolators (*H*_{s}) from model to
observational space as follows:

$$\begin{array}{}\text{(3)}& H\left(\mathbf{x}\right)={H}_{\mathrm{s}}{H}_{\mathrm{v}}\left(\mathbf{x}\right).\end{array}$$

We developed two types of variable converters, a single space interpolator, and their tangent linear and adjoint operators. Both WRF Var and NHM-4DVAR consider only perturbation to the mixing ratio of the rainwater and not to its number density in the tangent linear and adjoint models. However, in the tangent and adjoint operators described here (Sect. 2.2), the non-perturbed number density of rainwater is included. This variable is initialized to zero at the beginning or end of the operators, and this effect is directly considered in the cost functions of WRF Var and NHM-4DVAR, whereas its gradient is indirectly considered through perturbations of the mixing ratios of rainwater, water vapor, and other variables like temperature and pressure.

It is recommended that users of WRF Var run the system with CLOUD_CV (required) and the CV7 (optional) switches. The former adds the mixing ratios of rainwater to the default control variable set (Wang et al., 2013), and the latter replaces the control variables of stream function and velocity potential with momentum control variables to improve the performance of WRF simulations at high horizontal resolution (Sun et al., 2016). With these selections, the control variables in WRF Var are almost the same as those in NHM-4DVAR (Kawabata et al., 2011).

Among the many numerical precipitation scheme options (e.g., single-moment
scheme, large-scale condensation scheme) for WRF and JMANHM, we chose
double-moment schemes (WRF, Morrison et al., 2009; JMANHM, Hashimoto, 2008)
for our observational operators because such schemes predict both the number
density (*N*_{r}; m^{−3}) and the mixing ratio (*Q*_{r}; kg kg^{−1}) of
rainwater, whereas single-moment schemes predict only *Q*_{r}. Therefore, two
of three unknown parameters in the drop size distribution (DSD) function are
detected by the schemes. Following Morrison et al. (2009), the DSD function
is given by

$$\begin{array}{}\text{(4)}& N\left(D\right)={N}_{\mathrm{0}}{D}^{\mathit{\mu}}\mathrm{exp}\left(-\mathrm{\Lambda}D\right),\end{array}$$

where *D* (mm) is the raindrop diameter, *N*_{0} (mm^{−1} m^{−3}) is the
intercept parameter, *μ* is the shape parameter, and Λ
(mm^{−1}) is the slope parameter. Λ is given by

$$\begin{array}{}\text{(5)}& \mathrm{\Lambda}={\left({\displaystyle \frac{\mathit{\pi}{\mathit{\rho}}_{\mathrm{w}}{N}_{\mathrm{r}}}{{\mathrm{10}}^{\mathrm{3}}{\mathit{\rho}}_{\mathrm{a}}{Q}_{\mathrm{r}}}}\right)}^{{\scriptscriptstyle \frac{\mathrm{1}}{\mathrm{3}}}},\end{array}$$

where *ρ*_{w} is the density of water (997 kg m^{−3} in this study)
and *ρ*_{a} is air density (kg m^{−3}), a model diagnostic variable.
*ρ*_{a} and *N*_{0} are given by

$$\begin{array}{}\text{(6)}& {\displaystyle}& {\displaystyle}{\mathit{\rho}}_{\mathrm{a}}={\displaystyle \frac{p}{RT(\mathrm{1}+\mathrm{0.61}{q}_{\mathrm{v}})}},\text{(7)}& {\displaystyle}& {\displaystyle}{N}_{\mathrm{0}}={N}_{\mathrm{r}}\mathrm{\Lambda},\end{array}$$

where *p* is atmospheric pressure (Pa), *R* is the gas constant, *T* is temperature
(K), and *q*_{v} is the mixing ratio of water vapor (kg kg^{−1}).

In our study, the remaining unknown parameter *μ* is fixed at zero, and
*N*(*D*) is based on bulk sampling; the minimum and maximum values of *D* are set to
0.05 and 5 mm, respectively.

Because in the rainwater prognostic variables raindrops are assumed to be
spherical in both WRF and JMANHM, we introduce the axis ratio of a raindrop,
which is polynomial to *D* (Brandes et al., 2002, 2005), as follows:

$$\begin{array}{ll}{\displaystyle}& {\displaystyle}r=\mathrm{0.9951}+\mathrm{2.51}\times {\mathrm{10}}^{-\mathrm{2}}D-\mathrm{3.644}\times {\mathrm{10}}^{-\mathrm{2}}{D}^{\mathrm{2}}\\ \text{(8)}& {\displaystyle}& {\displaystyle}\phantom{\rule{1em}{0ex}}+\mathrm{5.303}\times {\mathrm{10}}^{-\mathrm{3}}{D}^{\mathrm{3}}-\mathrm{2.492}\times {\mathrm{10}}^{-\mathrm{4}}{D}^{\mathrm{4}}.\end{array}$$

Radar observations are derived from measurements of the scattering of
electromagnetic waves by raindrops. The first converter is based on fitting
functions that relate equivolume diameters *D* to scattering amplitude (Zhang
et al., 2001). The backscattering amplitudes are represented by a power-law
function as follows:

$$\begin{array}{}\text{(9)}& \left|{S}_{\mathrm{h},\mathrm{v}}\left(D\right)\right|={\mathit{\alpha}}_{\mathrm{h},\mathrm{v}}{D}^{{\mathit{\beta}}_{\mathrm{h},\mathrm{v}}},\end{array}$$

where the coefficients *α*_{h,v} and *β*_{h,v} are
determined by fitting *D* to the backscattering amplitudes |*S*_{h,v}| calculated by the
*T*-matrix method (Mishchenko et al., 1996). The difference between the
horizontal and vertical forward-scattering amplitudes is defined as

$$\begin{array}{}\text{(10)}& \mathit{\text{Re}}\left({f}_{\mathrm{h}}\left(D\right)-{f}_{\mathrm{v}}\left(D\right)\right)={\mathit{\alpha}}_{\mathrm{k}}{D}^{{\mathit{\beta}}_{\mathrm{k}}},\end{array}$$

where *f*_{h}(D) and *f*_{v}(D) represent the
horizontal and vertical forward-scattering amplitudes, and *α*_{k} and
*β*_{k} are determined by the fitting. Zhang et al. (2001) proposed
fitting functions for S-band radars, and Kawabata et al. (2018) derived new
fitting parameters for C-band radars. Following Zhang et al. (2001),
horizontal (H) and vertical (V) reflectivity factors are

$$\begin{array}{}\text{(11)}& {\displaystyle}{Z}_{\mathrm{H},\mathrm{V}}={\displaystyle \frac{\mathrm{4}{\mathit{\lambda}}^{\mathrm{4}}}{{\mathit{\pi}}^{\mathrm{4}}{\left|{K}_{\mathrm{w}}\right|}^{\mathrm{2}}}}\left({\mathit{\alpha}}_{\mathrm{h},\mathrm{v}}^{\mathrm{2}}{N}_{\mathrm{0}}{\mathrm{\Lambda}}^{-\left(\mathrm{2}{\mathit{\beta}}_{\mathrm{h},\mathrm{v}}+\mathrm{1}\right)}\mathrm{\Gamma}\left(\mathrm{2}{\mathit{\beta}}_{\mathrm{h},\mathrm{v}}+\mathrm{1}\right)\right),\end{array}$$

where *λ* (m) is the radar wavelength, *K*_{w} is a constant, defined
as ${K}_{\mathrm{w}}=\left(\mathit{\epsilon}-\mathrm{1}\right)/\left(\mathit{\epsilon}+\mathrm{2}\right)$, where
*ε* is the complex dielectric constant of water estimated as a
function of wavelength and temperature (Sadiku, 1985), and Γ
represents the Gamma function. The horizontal reflectivity *Z*_{H} is
converted to conventional reflectivity *Z*_{h} (dBZ) by

$$\begin{array}{}\text{(12)}& {Z}_{\mathrm{h}}=\mathrm{10}{\mathrm{log}}_{\mathrm{10}}\left({Z}_{\mathrm{H}}\right),\end{array}$$

and *Z*_{DR} (dB) is defined as

$$\begin{array}{}\text{(13)}& {\displaystyle}{Z}_{\mathrm{DR}}=\mathrm{10}{\mathrm{log}}_{\mathrm{10}}\left({Z}_{\mathrm{H}}/{Z}_{\mathrm{v}}\right)={Z}_{\mathrm{h}}-{Z}_{\mathrm{v}}.\end{array}$$

*K*_{DP} (^{∘} km^{−1}) is defined as

$$\begin{array}{}\text{(14)}& {\displaystyle}{K}_{\mathrm{DP}}={\displaystyle \frac{\mathrm{180}\mathit{\lambda}}{\mathit{\pi}}}{N}_{\mathrm{0}}{\mathit{\alpha}}_{\mathrm{k}}{\mathrm{\Lambda}}^{-\left({\mathit{\beta}}_{\mathrm{k}}+\mathrm{1}\right)}\mathrm{\Gamma}\left({\mathit{\beta}}_{\mathrm{k}}+\mathrm{1}\right).\end{array}$$

The attenuation effects are calculated as follows:

$$\begin{array}{}\text{(15)}& {\displaystyle}& {\displaystyle}{Z}_{\mathrm{h}}^{\mathrm{att}}\left(x\right)={Z}_{\mathrm{h}}\left(x\right)-\mathrm{2}\underset{\mathrm{0}}{\overset{x}{\int}}{A}_{\mathrm{H}}\left(s\right)\mathrm{d}s,\text{(16)}& {\displaystyle}& {\displaystyle}{Z}_{\mathrm{DR}}^{\mathrm{att}}\left(x\right)={Z}_{\mathrm{DR}}\left(x\right)-\mathrm{2}\underset{\mathrm{0}}{\overset{x}{\int}}{A}_{\mathrm{DP}}\left(s\right)\mathrm{d}s,\end{array}$$

where ${Z}_{\mathrm{h}}^{\mathrm{att}}$ and ${Z}_{\mathrm{DR}}^{\mathrm{att}}$ represent
attenuated *Z*_{h} and *Z*_{DR}, respectively. *A*_{H} and *A*_{DP} are the
specific attenuation (dB km^{−1}) and the specific differential
attenuation (dB km^{−1}), respectively, defined as

$$\begin{array}{}\text{(17)}& {\displaystyle}& {\displaystyle}{A}_{\mathrm{H}}={\mathit{\alpha}}_{\mathrm{H}}{K}_{\mathrm{DP}}^{{\mathit{\beta}}_{\mathrm{H}}},\text{(18)}& {\displaystyle}& {\displaystyle}{A}_{\mathrm{DP}}={\mathit{\alpha}}_{\mathrm{d}}{K}_{\mathrm{DP}}^{{\mathit{\beta}}_{\mathrm{d}}}.\end{array}$$

The values of the coefficients *α*_{h}, *α*_{v},
*α*_{k}, *α*_{H}, and *α*_{d} and *β*_{h}, *β*_{v}, *β*_{k}, *β*_{H}, and
*β*_{d} for C-band in these equations are listed in Table 1.
Hereafter, this converter is called FIT.

FIT is also applicable for X- and S-bands by replacing the coefficients. Although we already prepared the coefficients for all bands in the source codes, the users should carefully investigate their validity.

The second converter (hereafter KD) converts observed *K*_{DP} to rainwater
content (*Q*_{rain}) according to the following relation:

$$\begin{array}{}\text{(19)}& {\displaystyle}{Q}_{\mathrm{rain}}=\mathrm{3.565}{\left({\displaystyle \frac{{K}_{\mathrm{DP}}}{f}}\right)}^{\mathrm{0.77}},\end{array}$$

where *f* (GHz) is the radar frequency and the power-law coefficients are from
Bringi and Chandrasekar (2001). *Q*_{rain} in the model is defined as
*Q*_{rain}=*Q*_{r}*ρ*_{a} (kg m^{−3}). Note that Eq. (19) is applicable
not only for C-band but also X- and S-bands by putting their frequencies in
*f*. Equations (4)–(19) follow Kawabata et al. (2018), and we put the equations with
different order in this paper for reader convenience to
understand the flow of implementations of the forward, tangent linear, and
adjoint codes.

Because only *p*, *T*, and *q*_{v} are perturbed in WRF Var and NHM-4DVAR, the
linearized form of Eq. (6) is

$$\begin{array}{ll}{\displaystyle}& {\displaystyle}\mathrm{\Delta}{\mathit{\rho}}_{\mathrm{a}}={\displaystyle \frac{\mathrm{\Delta}p}{RT\left(\mathrm{1}+\mathrm{0.61}{q}_{\mathrm{v}}\right)}}-{\displaystyle \frac{p\mathrm{\Delta}T}{R{T}^{\mathrm{2}}\left(\mathrm{1}+\mathrm{0.61}{q}_{\mathrm{v}}\right)}}\\ \text{(20)}& {\displaystyle}& {\displaystyle}\phantom{\rule{1em}{0ex}}-{\displaystyle \frac{\mathrm{0.61}p\mathrm{\Delta}{q}_{\mathrm{v}}}{RT{\left(\mathrm{1}+\mathrm{0.61}{q}_{\mathrm{v}}\right)}^{\mathrm{2}}}},\end{array}$$

and the perturbations of Λ and *N*_{0} are given as

$$\begin{array}{}\text{(21)}& {\displaystyle}& {\displaystyle}\mathrm{\Delta}\mathrm{\Lambda}={\displaystyle \frac{\mathrm{1}}{\mathrm{3}}}\mathrm{\Lambda}\left(-\mathrm{\Delta}{Q}_{\mathrm{r}}{Q}_{\mathrm{r}}^{-\mathrm{1}}-\mathrm{\Delta}{\mathit{\rho}}_{\mathrm{a}}{\mathit{\rho}}_{\mathrm{a}}^{-\mathrm{1}}\right),\text{(22)}& {\displaystyle}& {\displaystyle}\mathrm{\Delta}{N}_{\mathrm{0}}={N}_{\mathrm{r}}\mathrm{\Delta}\mathrm{\Lambda},\end{array}$$

where Δ*Q*_{r} and Δ*N*_{0} are perturbations of the
mixing ratio and number density of rainwater, respectively. Note that the
perturbation of *N*_{r} is not considered in the adjoint model (see
Sect. 2). Thus, the perturbations of *Z*_{H,V}, *Z*_{DR}, and
*K*_{DP} are represented as follows.

$$\begin{array}{ll}{\displaystyle}& {\displaystyle}\mathrm{\Delta}{Z}_{\mathrm{H},\mathrm{V}}={\displaystyle \frac{\mathrm{4}{\mathit{\lambda}}^{\mathrm{4}}}{{\mathit{\pi}}^{\mathrm{4}}{\left|{K}_{\mathrm{w}}\right|}^{\mathrm{2}}}}\left({\mathit{\alpha}}_{\mathrm{h},\mathrm{v}}^{\mathrm{2}}\mathrm{\Gamma}\left(\mathrm{2}{\mathit{\beta}}_{\mathrm{h},\mathrm{v}}+\mathrm{1}\right)\left(\mathrm{\Delta}{N}_{\mathrm{0}}{\mathrm{\Lambda}}^{-\left(\mathrm{2}{\mathit{\beta}}_{\mathrm{h},\mathrm{v}}+\mathrm{1}\right)}\right.\right.\\ \text{(23)}& {\displaystyle}& {\displaystyle}\phantom{\rule{1em}{0ex}}\left.\left.-\left(\mathrm{2}{\mathit{\beta}}_{\mathrm{h},\mathrm{v}}+\mathrm{1}\right)\mathrm{\Delta}\mathrm{\Lambda}{N}_{\mathrm{0}}{\mathrm{\Lambda}}^{-\left(\mathrm{2}{\mathit{\beta}}_{\mathrm{h},\mathrm{v}}+\mathrm{2}\right)}\right)\right)\text{(24)}& {\displaystyle}& {\displaystyle}\mathrm{\Delta}{Z}_{\mathrm{DR}}=\mathrm{\Delta}{Z}_{\mathrm{h}}-\mathrm{\Delta}{Z}_{\mathrm{v}}{\displaystyle}& {\displaystyle}\mathrm{\Delta}{K}_{\mathrm{DP}}={\displaystyle \frac{\mathrm{180}\mathit{\lambda}}{\mathit{\pi}}}{\mathit{\alpha}}_{\mathrm{k}}\mathrm{\Gamma}\left({\mathit{\beta}}_{\mathrm{k}}+\mathrm{1}\right)\left(\mathrm{\Delta}{N}_{\mathrm{0}}{\mathrm{\Lambda}}^{-\left({\mathit{\beta}}_{\mathrm{k}}+\mathrm{1}\right)}\right.\\ \text{(25)}& {\displaystyle}& {\displaystyle}\phantom{\rule{1em}{0ex}}\left.-\left({\mathit{\beta}}_{\mathrm{k}}+\mathrm{1}\right)\mathrm{\Delta}\mathrm{\Lambda}{N}_{\mathrm{0}}{\mathrm{\Lambda}}^{-\left({\mathit{\beta}}_{\mathrm{k}}+\mathrm{2}\right)}\right)\end{array}$$

Finally, the perturbations of *A*_{H} and *A*_{DP} are

$$\begin{array}{}\text{(26)}& {\displaystyle}& {\displaystyle}\mathrm{\Delta}{A}_{\mathrm{H}}={\mathit{\alpha}}_{\mathrm{H}}{\mathit{\beta}}_{\mathrm{H}}\mathrm{\Delta}{K}_{\mathrm{DP}}{{K}_{\mathrm{DP}}}^{{\mathit{\beta}}_{\mathrm{H}}-\mathrm{1}},\text{(27)}& {\displaystyle}& {\displaystyle}\mathrm{\Delta}{A}_{\mathrm{DP}}={\mathit{\alpha}}_{\mathrm{d}}{\mathit{\beta}}_{\mathrm{d}}\mathrm{\Delta}{K}_{\mathrm{DP}}{{K}_{\mathrm{DP}}}^{{\mathit{\beta}}_{\mathrm{d}}-\mathrm{1}}.\end{array}$$

The adjoint operators are represented by the transposed form of Eqs. (20)–(27), which is (tangent linear)^{T}.
As an example, the adjoint of Eq. (27) is

$$\begin{array}{}\text{(28)}& \mathrm{\Delta}{K}_{\mathrm{DP}}=\mathrm{\Delta}{K}_{\mathrm{DP}}+{\mathit{\alpha}}_{\mathrm{d}}{\mathit{\beta}}_{\mathrm{d}}{{K}_{\mathrm{DP}}}^{{\mathit{\beta}}_{\mathrm{d}}-\mathrm{1}}\mathrm{\Delta}{A}_{\mathrm{DP}}.\end{array}$$

Because *K*_{DP} in Eq. (19) is an observed value, it is not necessary to
linearize the equation. However, the equation that relates *Q*_{rain} to
*Q*_{r} (Sect. 2.1.2) needs to be linearized as follows:

$$\begin{array}{}\text{(29)}& \mathrm{\Delta}{Q}_{\mathrm{rain}}=\mathrm{\Delta}{Q}_{\mathrm{r}}{\mathit{\rho}}_{\mathrm{a}}+{Q}_{\mathrm{r}}\mathrm{\Delta}{\mathit{\rho}}_{\mathrm{a}}.\end{array}$$

The transposed form of this equation is used for the adjoint model (see Sect. 2.2.1).

Space interpolators in data assimilation systems map the model space to the
observational space according to the representativeness of the observations.
In the case of radar data, the effect of beam broadening stands for the
representativeness, typically for a beam width of approximately
1.0^{∘}. The broadening is characterized by a Gaussian distribution
orthogonal to the direction of the radar beam. Most previous studies (e.g.,
Seko et al., 2004; Wattrelot et al., 2014), except Zeng et al. (2016),
consider only vertical beam broadening because numerical models have
horizontal grid spacings of several kilometers, whereas they have vertical
grid spacings in the lower troposphere of less than 1 km. However,
data assimilation systems must have sub-kilometer horizontal grid spacings
as well (e.g., Kawabata et al., 2014a; Miyoshi et al., 2016) so that the
space interpolators can take account of horizontal beam broadening. In
addition, several phased-array radars recently deployed in Japan have
different beam widths in the vertical and horizontal directions. Our
operator thus considers beam broadening in both the vertical and horizontal
directions.

In addition, it is important for the space interpolator to include
beam-bending effects, which depend on atmospheric conditions. In this study,
the bending is determined by considering the climatological vertical
gradient of the refractive index of the atmosphere in accordance with the
effective Earth radius model (Doviak and Zrnić, 1993) following Haase
and Crewell (2000), who showed statistically that the climatological
refractive index is close to the actual refractive index at elevation angles
higher than 1^{∘} instead of by considering the actual atmospheric
conditions, although Zeng et al. (2014) developed an excellent radar
simulator that considers the actual refractivity of the atmosphere.

Remote sensing observations usually have higher spatial resolutions than model grid spacings. To avoid correlations of the observational errors in such high-resolution data, it is necessary either to thin the data or to use “super observations”. In this study, we chose the super observation method, in which observations are averaged over each model grid cell. Super observation methods also have the advantage that they remove undesirable fluctuations associated with sub-grid-scale phenomena, the assimilation of which makes the numerical model unnecessarily noisy (e.g., Seko et al., 2004; Zhang et al., 2009).

First, we calculated the path of the center of the radar beam in the model domain, including its elevation, azimuth, and bending angles (Fig. 1a). Once sufficient data are included within a model grid cell, they are averaged and mapped onto an interpolation point along the radar beam (IP in Fig. 1). This value at this point is a “super observation”, and it is compared with the modeled value, which was interpolated by using Gaussian weights (Fig. 1b). Moreover, we also developed the tangent and adjoint codes of the space interpolator.

3 Setup options

Back to toptop
The operators are controlled by the namelist (“namelist.polradar”) as
follows.

&name_obs o_dir='/home/usr/datadir',
o_stn(1)='OFT', o_stn(2)='TUR', icnv=0/

Here, “o_dir” is the directory for the input observational
data, “o_stn” indicates the station names of radar sites,
where “max_stn”, the number of names, is set in
“da_setup_obs_structures_polradar.inc” in WRF Var and in
“obs_dual_pol.f90:” in NHM-4DVAR, and “icnv”
is a switch for the selection of the observational operator, where “0” and
“5” mean FIT and KD, respectively.

In addition, a file that defines for each radar area where the beam is
blocked by topography, named “beam_block_rate_$*{*radar_site*}*.dat”, must be supplied by the user. This file is made
by another program and should be prepared before the assimilation.

4 Results

Back to toptop
In this section, we examine the linearity of only the FIT variable converter; it is not necessary to examine the linearity of the KD converter because of the intrinsic linearity of Eq. (19). We evaluated the linearity of FIT by performing a Taylor expansion. If the original equation is given as

$$\begin{array}{}\text{(30)}& \mathbf{y}=H\left(\mathbf{x}\right),\end{array}$$

then the linearized equation is defined as

$$\begin{array}{}\text{(31)}& \mathit{\delta}\mathbf{y}=\mathbf{H}\mathit{\delta}\mathbf{x}.\end{array}$$

If the linear equation is derived with no errors, the following Taylor expansion of Eq. (31),

$$\begin{array}{}\text{(32)}& {\displaystyle \frac{\mathrm{|}H\left(\mathbf{x}+\mathit{\alpha}\mathit{\delta}\mathbf{x}\right)-H\left(\mathbf{x}\right)\mathrm{|}}{\left|\mathit{\alpha}\right|\mathrm{|}\mathbf{H}\mathit{\delta}\mathbf{x}\mathrm{|}}}=\mathrm{1}+O\left(\mathit{\alpha}\right),\end{array}$$

should be accurate within the rounding error of the computer. The results
for *Z*_{H}, *Z*_{V}, and *K*_{DP} in Eqs. (11) and (14) are 1.00 when *α* is 10^{−7} to 10^{−15}.

Regarding the adjoint operator, we evaluated the following equation:

$$\begin{array}{}\text{(33)}& {\left(\mathbf{H}\mathit{\delta}\mathbf{x}\right)}^{\mathrm{T}}\left(\mathbf{H}\mathit{\delta}\mathbf{x}\right)=\mathit{\delta}{\mathbf{x}}^{\mathrm{T}}\left[{\mathbf{H}}^{\mathrm{T}}\left(\mathbf{H}\mathit{\delta}\mathbf{x}\right)\right],\end{array}$$

where the left-hand side of Eq. (33) is calculated using the tangent linear
operator, and on the right-hand side, the output variables of the tangent
linear operator are input into the adjoint operator. This equation must be
accurate within the rounding error. In FIT, the difference between the left-
and right-hand sides was −8.215650382 × 10^{−15}, which we
consider accurate enough.

We conducted two simple data assimilation tests. Observational errors of
*Z*_{h}, *Z*_{DR}, *K*_{DP}, and *Q*_{rain}, which were determined after the
statistical examination (Kawabata et al., 2018), were 15.0 dBZ, which is the
same as in Kawabata et al. (2011) at 2.0 dB, 4.0^{∘} km^{−1}, and 4 g m^{−3}, respectively. These errors are homogeneous in space, which means
observational error covariances are diagonal.

The first one was done using NHM-4DVAR with actual radar data from the C-band
dual polarimetric radar at the Meteorological Research Institute in Tsukuba,
Japan (Yamauchi et al., 2012; Adachi et al., 2013). In this experiment, both
radial velocity data and the polarimetric parameters of *Z*_{H},
*Z*_{DR}, and *K*_{DP} were assimilated in FIT, and radial
velocity and *Q*_{rain} derived from *K*_{DP} were assimilated
in KD. The assimilation window was from 21:00 to 21:05 UTC on 23 June 2014, a
day on which intense hail fell in Tokyo, Japan. The horizontal resolution of
NHM-4DVAR was 2 km and the length of the assimilation window was 5 min, 11
PPI data from 0.5 to 4.8^{∘} elevations with the azimuth resolution of
0.7^{∘} and the range resolution of 150 m were assimilated. PPI
data were assimilated at the exact observation time as far as the time interval of
NHM-4DVAR (10 s in this case) permits. The background errors were described
in Kawabata et al. (2007, 2011). Analysis (KD and FIT) and observational
(OBS) fields of *Z*_{h}, *Z*_{DR}, and*K*_{DP} are
shown in Fig. 2, which displays the whole assimilation domain. Although there
was no rain region in the first-guess field (FG; Fig. 2d), *Z*_{h} in
KD was comparable to that in OBS from the standpoint of rainfall distribution
and intensity, but *Z*_{h} in FIT covered a much smaller area than it
did in OBS. This smaller coverage may be due to nonlinearity in FIT. In KD,
we can see quite small values of *K*_{DP} (Fig. 2f), but good
agreement with OBS in its horizontal distribution, while *Z*_{h} looks
better than *K*_{DP}. *K*_{DP} values were smaller in both KD
and FIT than in OBS. This result is similar to that of a statistical analysis
performed by Kawabata et al. (2018). In contrast, *Z*_{DR} values in
KD and FIT were larger than OBS over large areas. This result implies that
the calculation of the axis ratio of raindrops (Eq. 8) may need modification
because in the FG field, *Z*_{DR} values and coverage were already too
large in comparison with those of OBS.

The second one was done using WRF 3D-Var with actual radar data from the DWD
radar network (Helmert et al., 2014) for the same case with “Case 1”
described in Kawabata et al. (2018). The horizontal resolution of WRF 3D-Var
was 2 km, and polarimetric parameters and rainwater content in single PPI
data by Offenthal radar were assimilated (see Kawabata et al., 2018, for
detailed information on the observation). The background errors were
calculated with ensemble simulations by WRF initialized by ECMWF analysis
using the “gen_be” tool compiled in WRFDA. Observational
errors were the same as the first case. From the increments of
polarimetric parameters (Fig. 3), although quite small impacts are seen,
similar patterns are recognized in both methods and larger impacts of
*Z*_{h} and *Z*_{DR} were produced in FIT and KD, respectively.

In both cases, the radial velocity data were assimilated with the same method as in Sun and Crook (1997).

5 Summary

Back to toptop
We implemented two variable converters for polarimetric radars in the WRF variational data assimilation system (WRF Var) and the JMANHM data assimilation system (NHM-4DVAR). FIT simulates polarimetric parameters using a double-moment cloud microphysics scheme, and KD estimates rainwater contents with the observed specific differential phase. Both FIT and KD are applicable not only for C-band but also X- and S-bands. The advantage of FIT over KD is that it includes theoretically precise formulations for both the mixing ratio and number density of rainwater, as well as attenuation effects, whereas KD has advantages due to its linear formulation and small computational cost.

These operators work in conjunction with an advanced space interpolator, which considers (1) beam broadening in three dimensions, (2) different beam widths in the vertical and horizontal directions, and (3) the climatological beam-bending effect. The interpolator also simulated attenuation effects.

Tangent and adjoint operators of the two variable converters and the space
interpolator were developed and implemented along with the forward
operators. In a simple data assimilation experiment, we succeeded in
assimilating actual polarimetric observations and obtained reasonable
results with both the FIT and KD operators, except for *Z*_{DR}. However, our
results show a need for further improvements of the *K*_{DP} and *Z*_{DR}
estimates. It would be possible to overcome the weaknesses of the *Z*_{h}
distributions in FIT and FG through assimilation–forecast cycles and/or by
adding other types of observation data, such as conventional observations,
Doppler (water vapor) lidar data, and water vapor data observed by GNSS.
Furthermore, it is necessary to improve quality controls (QCs) for
polarimetric parameters, although the same QCs were applied as described in
Kawabata et al. (2018) and the impact of axis ratio (Eq. 8) and observational
errors on assimilations will be investigated, and it is necessary to
estimate more appropriate observational errors (e.g., Wulfmeyer et al.,
2016). These challenges would improve QPE and QPF with the current forms of
the operators.

Code and data availability

Back to toptop
Code and data availability.

Since PolRad VAR v1.0 for NHM-4DVAR belongs to the Meteorological Research Institute of the Japan Meteorological Agency and is not publicly available, any researchers interested in the code are encouraged to contact the corresponding author and sign a contract for license to get the code. PolRad VAR v1.0 for WRF Var is currently being implemented into the community version of WRF Var and will be accessible at the WRF repository (http://www2.mmm.ucar.edu/wrf/users/downloads.html, last access: 19 June 2018) in the near future. Any researchers interested in the current form of the code can get it from the corresponding author via e-mail.

Competing interests

Back to toptop
Competing interests.

The authors declare that they have no conflict of interest.

Acknowledgements

Back to toptop
Acknowledgements.

We are grateful to the WRF development and support team. This study was
partly supported by the Catchments As Organized System (CAOS) project funded
by the German Research Foundation (FOR 1598), the Japanese Ministry of
Education, Culture, Sports, Science and Technology “Advancement of
meteorological and global environmental predictions utilizing observational
Big Data” project, the Japan Science and Technology Agency (CREST)
“Innovating Big Data Assimilation technology for revolutionizing
very-short-range severe weather prediction” project (JPMJCR1312), and
Grants-in-Aid for Scientific Research no. 16H04054 “Study of optimum
perturbation methods for ensemble data assimilation” and no. 17H02962
“Study on uncertainties in convection initiations and developments using a
particle filter” from the Japan Society for the Promotion of Science.
We also express our sincere
thanks to the reviewers, Juanzhen Sun and Soichiro Sugimoto, and the editor, Josef Koller, for their valuable contributions to this paper.

Edited by: Josef Koller

Reviewed by: Juanzhen Sun and Soichiro Sugimoto

References

Back to toptop
Adachi, A., Kobayashi, T., Yamauchi, H., and Onogi, S.: Detection of potentially hazardous convective clouds with a dual-polarized C-band radar, Atmos. Meas. Tech., 6, 2741–2760, https://doi.org/10.5194/amt-6-2741-2013, 2013.

Adachi, A., Kobayashi, T., and Yamauchi, H.: Estimation of raindrop size distribution and rainfall rate from polarimetric radar measurements at attenuating frequency based on the self-consistency principle, J. Meteorol. Soc. Jpn., 93, 359–388, 2015.

Anagnostou, M. N., Anagnostou, E. N., Vivekanandan, J., and Ogden, F. L.: Comparison of two raindrop size distribution retrieval algorithms for X-band dual polarization observations, J. Hydrometeorol., 9, 589–600, 2008.

Barker, D., Huang, X.-Y., Liu, Z., Auligné, T., Zhang, X., Rugg, S., Ajjaji, R., Bourgeois, A., Bray, J., Chen, Y., Demirtas, M., Guo, Y.-R., Henderson, T., Huang, W., Lin, H.-C., Michalakes, J., Rizvi, S., and Zhang, X.: The Weather Research and Forecasting Model's Community Variational/Ensemble Data Assimilation System: WRFDA, B. Am. Meteorol. Soc., 93, 831–843, 2012.

Bauer, H.-S., Schwitalla, T., Wulfmeyer, V., Bakhshaii, A., Ehret, U., Neuper, M., and Caumont, O.: Quantitative precipitation estimation based on high-resolution numerical weather prediction and data assimilation with WRF – a performance test, Tellus A, 67, 25047, https://doi.org/10.3402/tellusa.v67.25047, 2015.

Brandes, E. A., Zhang, G., and Vivekanandan, J.: Experiments in rainfall estimation with a polarimetric radar in a subtropical environment, J. Appl. Meteorol., 41, 674–685, 2002.

Brandes, E. A., Zhang, G., and Vivekanandan, J.: Corrigendum, J. Appl. Meteorol., 44, 186–186, 2005.

Bringi, V. N. and Chandrasekar, V.: Polarimetric Doppler weather radar-principle and application, Cambridge University Press, Cambridge, UK, 636 pp., 2001.

Doviak, R. J. and Zrnić, D. S.: Doppler radar and weather observations, Second edition, Dover publications, New York, USA, 562 pp., 1993.

Haase, G. and Crewell, S.: Simulation of radar reflectivities using a mesoscale weather forecast model, Water Resour. Res., 36, 2221–2231, 2000.

Hashimoto, A.: A double-moment cloud microphysical parameterization, Supplemental report at Numerical Prediction Division at Meteorological Agency, 54, 81–92, 2008 (in Japanese).

Helmert, K., Tracksdorf, P., Steinert, J., Werner, M., Frech, M., Rathmann, N., Hengstebeck, T., Mott, M., Schumann, S., and Mammen, T.: DWDs new radar network and postprocessing algorithm chain, Abstract of 8th European Conference on Radar in Meteorology and Hydrology, 1–5 September 2014, Garmisch-Partenkirchen, Germany, 6 pp., available at: http://www.pa.op.dlr.de/erad2014/programme/ExtendedAbstracts/237_Helmert.pdf (last access: 19 June 2018), 2014.

Jameson, A. R.: A comparison of microwave techniques for measuring rainfall, J. Appl. Meteorol., 30, 32–54, 1991.

Jameson, A. R. and Caylor, I. J.: A new approach to estimating rainwater content by radar using propagation differential phase shift, J. Appl. Meteorol. Clim., 11, 311–322, 1994.

Jung, Y., Zhang, G., and 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., 136, 2228–2245, 2008a.

Jung, Y., Xue, M., Zhang, G., and Straka, J. M.: Assimilation of simulated polarimetric radar data for a convective storm using the ensemble Kalman filter. Part II: impact of polarimetric data on storm analysis, Mon. Weather Rev., 136, 2246–2260, 2008b.

Kawabata, T., Seko, H., Saito, K., Kuroda, T., Tamiya, K., Tsuyuki, T., Honda, Y., and Wakazuki, Y.: An Assimilation and Forecasting Experiment of the Nerima Heavy Rainfall with a Cloud-Resolving Nonhydrostatic 4-Dimensional Variational Data Assimilation System, J. Meteorol. Soc. Jpn., 85, 255–276, https://doi.org/10.2151/jmsj.85.255, 2007.

Kawabata, T., Kuroda, T., Seko, H., and Saito, K.: A cloud-resolving 4D-Var assimilation experiment for a local heavy rainfall event in the Tokyo metropolitan area, Mon. Weather Rev., 139, 1911–1931, 2011.

Kawabata, T., Shoji, Y., Seko, H., and Saito, K.: A numerical study on a mesoscale convective system over a subtropical island with 4D-Var assimilation of GPS slant total delays, J. Meteorol. Soc. Jpn., 91, 705–721, 2013.

Kawabata, T., Ito, K., and Saito, K.: Recent Progress of the NHM-4DVAR towards a Super-High Resolution Data Assimilation, SOLA, 10, 145–149, 2014a.

Kawabata, T., Iwai, H., Seko, H., Shoji, Y., Saito, K., Ishii, S., and Mizutani, K.: Cloud-Resolving 4D-Var Assimilation of Doppler Wind Lidar Data on a Meso-Gamma-Scale Convective System, Mon. Weather Rev., 142, 4484–449, 2014b.

Kawabata, T., Bauer, H.-S., Schwitalla, T., Wulfmeyer, V., and Adachi, A.: Evaluation of forward operators for polarimetric radars aiming for data assimilation, J. Meteorol. Soc. Jpn., 96A, 157–174, 2018.

Kim, D. S., Maki, M., and Lee, D. I.: Retrieval of three-dimensional raindrop size distribution using X-band polarimetric radar data, J. Atmos. Ocean. Tech., 27, 1265–1285, 2010.

Mishchenko, M. I., Travis, L. D., and Mackowski, D. W.: T-matrix computations of light scattering by nonspherical particles: A review, J. Quant. Spectrosc. Ra., 55, 535–575, 1996.

Miyoshi, T., Kunii, M., Ruiz, J., Lien, G.-Y., Satoh, S., Ushio, T., Bessho, K., Seko, H., Tomita, H., and Ishikawa, Y.: “Big Data Assimilation” Revolutionizing Severe Weather Prediction, B. Am. Meteorol. Soc., 97, 1347–1354, 2016.

Morrison, H., Thompson, G., and Tatarskii, V.: Impact of cloud microphysics on the development of trailing stratiform precipitation in a simulated squall line: Comparison of one- and two-moment schemes, Mon. Weather Rev., 137, 991–1007, 2009.

Ryzhkov, A., Diederich, M., Zhang, P., and Simmer, C.: Potential utilization of specific attenuation for rainfall estimation, mitigation of partial beam blockage, and radar networking., J. Atmos. Ocean. Tech., 31, 599–619, 2014.

Ryzhkov, A. V. and Zrnić, D. S.: Comparison of dual-polarization radar estimators of rain, J. Atmos. Ocean. Tech., 12, 249–256, 1995.

Sadiku, M. N. O.: Refractive index of snow at microwave frequencies, Appl. Optics, 24, 572–575, 1985.

Saito, K.: The JMA nonhydrostatic model and its application to operation and research. Atmospheric Model Applications, edited by: Yucel, I., Intech, London, UK, 85–110, 2012.

Schwitalla, T. and Wulfmeyer, V.: Radar data assimilation experiments using the IPM WRF Rapid Update Cycle, Meteorol. Z., 23, 79–102, https://doi.org/10.1127/0941-2948/2014/0513, 2014.

Schwitalla, T., Bauer, H.-S., Wulfmeyer, V., and Aoshima, F.: High-resolution simu-lation over central Europe: Assimilation experiments with WRF 3DVAR during COPS IOP9c, Q. J. Roy. Meteor. Soc., 137, 156–175, https://doi.org/10.1002/qj.721, 2011.

Seko, H., Kawabata, T., Tsuyuki, T., Nakamura, H., Koizumi, K., and Iwabuchi, T.: Impacts of GPS-derived Water Vapor and Radial Wind Measured by Doppler Radar on Numerical Prediction of Precipitation, J. Meteorol. Soc. Jpn., 82, 473–489, 2004.

Skamarock, W., C., Klemp, J. B., Dudhia, J., Gill, D. O., Barker, D. M., Duda, M. G., Huang, X.-Y., Wang, W., and Powers, J. G.: A description of the advanced research WRF version 3, NCAR Tech. Note, National Center for Atmospheric Research, Colorado, USA, 2008.

Sun, J. and 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., 54, 1642–1661, 1997.

Sun, J., Wang, H., Tong, W., Zhang, Y., Lin, C.-Y., and Xu, D.: Comparison of the Impacts of Momentum Control Variables on High-Resolution Variational Data Assimilation and Precipitation Forecasting, Mon. Weather Rev., 144, 149–169, 2016.

Wang, H., Sun, J., Fan, S., and Huang, X.-Y.: Indirect Assimilation of Radar Reflectivity with WRF 3D-Var and Its Impact on Prediction of Four Summertime Convective Events, J. Appl. Meteorol. Clim., 52, 889–902, 2013.

Wattrelot, E., Caumont, O., and Mahfouf, J.-F.: Operational Implementation of the 1D+3D-Var Assimilation Method of Radar Reflectivity Data in the AROME Model, Mon. Weather Rev., 142, 1852–1873, 2014.

Wulfmeyer, V., Muppa, S. K., Behrendt, A., Hammann, E., Späth, F., Sorbjan, Z., Turner, D. D., and Hardesty, R. M.: Determination of Convective Boundary Layer Entrainment Fluxes, Dissipation Rates, and the Molecular Destruction of Variances: Theoretical Description and a Strategy for Its Confirmation with a Novel Lidar System Synergy, J. Atmos. Sci., 73, 667–692, https://doi.org/10.1175/JAS-D-14-0392.1, 2016.

Yamauchi, H., Adachi, A., Suzuki, O., and Kobayashi, T.: Precipitation estimate of a heavy rain event using a C-band solid-state polarimetric radar, 7th European Conf. on Radar in Meteorology and Hydrology, 24–29 June 2012, Toulouse, France, available at: http://www.meteo.fr/cic/meetings/2012/ERAD/extended_abs/SP_412_ext_abs.pdf (last access: 19 June 2018), 2012.

Yokota, S., Seko, H., Kunii, M., Yamauchi, H., and Niino, H.: The tornadic supercell on the Kanto Plain on 6 May 2012: Polarimetric radar and surface data assimilation with EnKF and ensemble-based sensitivity analysis, Mon. Weather Rev., 144, 3133–3157, 2016.

Zeng, Y., Blahak, U., Neuper, M., and Jerger, D.: Radar Beam Tracing Methods Based on Atmospheric Refractive Index, J. Atmos. Ocean. Tech., 31, 2650–2670, 2014.

Zeng, Y., Blahak, U., and Epperlein, D.: An Efficient Modular Volume Scanning Radar Forward Operator for NWP-models: Description and coupling to the COSMO-model, Q. J. Roy. Meteor. Soc., 142, 3234–3256, 2016.

Zhang, F., Weng, Y., Sippel, J. A., Meng, Z., and Bishop, C. H.: Cloud-resolving hurricane initialization and prediction through assimilation of Doppler radar observations with an ensemble Kalman filter, Mon. Weather Rev., 137, 2105–2125, 2009.

Zhang, G., Vivekanandan, J., and Brandes, E.: A method for estimating rain rate and drop size distribution from polarimetric radar measurements, IEEE T. Geosci. Remote, 39, 830–841, 2001.

Short summary

We implemented two observational operators for dual polarimetric radars in two variational data assimilation systems: WRF Var and NHM-4DVAR. The operators consist of a space interpolator and two types of variable converters. The first variable converter emulates polarimetric parameters with model prognostic variables, and the second derives rainwater content from the observed polarimetric parameter. The system worked properly in verification and assimilation tests.

We implemented two observational operators for dual polarimetric radars in two variational data...

Geoscientific Model Development

An interactive open-access journal of the European Geosciences Union