# KDP estimation

### From Lrose Wiki

LROSE procedure to estimate $\displaystyle{ K_{DP} }$ from $\displaystyle{ \phi_{DP} }$

## Contents

### Overview

$\displaystyle{ K_{DP} }$ is estimated from the range derivative of $\displaystyle{ \phi_{DP} }$, which accumulates along a radar beam as it passes through a precipitating medium. Regions of higher $\displaystyle{ K_{DP} }$ correspond to regions of increased raindrop size and concentration. $\displaystyle{ \phi_{DP} }$ must be pre-processed so that $\displaystyle{ K_{DP} }$ is useful.

The LROSE method for estimating $\displaystyle{ K_{DP} }$ is based on an updated version of the technique developed by Hubbert and Bringi (1995).

### Background and Objectives

As a radar beam propagates through a precipitating medium, the electromagnetic wave slows down resulting in a phase shift. In radar volumes where the hydrometeors are not spherical, the phase shifts of the horizontal and vertical pulses are different. $\displaystyle{ \phi_{DP} }$ measures the difference between these phase shifts. However, $\displaystyle{ \phi_{DP} }$ accumulates along the beam, which makes the field difficult to interpret. Calculating the range derivative of $\displaystyle{ \phi_{DP} }$ produces the $\displaystyle{ K_{DP} }$ field, which highlights areas of strong changes in $\displaystyle{ \phi_{DP} }$.

In order to ensure the KDP field is useful, pre-processing of the $\displaystyle{ \phi_{DP} }$ field is necessary. First, $\displaystyle{ \phi_{DP} }$ is a noisy field. We want the $\displaystyle{ K_{DP} }$ algorithm to identify regions of substantial change in $\displaystyle{ \phi_{DP} }$. Second, although phase shifts in Rayleigh-scattering regimes are dominated by changes in the wave propagation speed, Mie-scattering regimes induce a additional phase shift component related to backscatter differential phase (or phase shift on backscatter). That is, where nonspherical hydrometeors are large compared to the radar wavelength, constructive and destructive interference in the backscattered signal due to resonance can produce local phase shifts that aren't directly related to the liquid water content slowing down the radar beam. Refer to Figure 3 of Stepanian et al. (2016) for a good visualization of Mie-scattering effects. Before $\displaystyle{ K_{DP} }$ is calculated, these two issues need to be accounted for.

The primary objectives of the LROSE $\displaystyle{ K_{DP} }$ computation are:

• smooth the $\displaystyle{ \phi_{DP} }$ profile
• remove local increases in $\displaystyle{ \phi_{DP} }$ caused by backscatter differential phase

### Algorithm

(1) unfold $\displaystyle{ \phi_{DP} }$.

(2) apply a finite impulse response (FIR) filter to unfolded $\displaystyle{ \phi_{DP} }$, a specified number of times, to create 'smoothed $\displaystyle{ \phi_{DP} }$'. This step is similar to Hubbert and Bringi (1995), but uses fewer iterations.

(3) for increasing range, identify peaks in 'smoothed $\displaystyle{ \phi_{DP} }$' which are followed by a minimum. Remove these peaks, by trimming them down to the minimum that follows. We call the result the 'conditioned $\displaystyle{ \phi_{DP} }$'. This step removes the local increases in $\displaystyle{ \phi_{DP} }$ due to backscatter differential phase.

(4) apply FIR filter to 'conditioned $\displaystyle{ \phi_{DP} }$' a number of times, to provide further smoothing.

(5) Compute $\displaystyle{ K_{DP} }$ as the (slope of conditioned $\displaystyle{ \phi_{DP} }$ in range) / 2.0.