skip to main content


Title: A Polarimetric Radar Operator and Application for Convective Storm Simulation
In this study, a polarimetric radar forward model operator was developed for the Weather Research and Forecasting (WRF) model that was based on a scattering algorithm using the T-matrix methodology. Three microphysics schemes—Thompson, Morrison 2-moment, and Milbrandt-Yau 2-moment—were supported in the operator. This radar forward operator used the microphysics, thermodynamic, and wind fields from WRF model forecasts to compute horizontal reflectivity, radial velocity, and polarimetric variables including differential reflectivity (ZDR) and specific differential phase (KDP) for S-band radar. A case study with severe convective storms was used to examine the accuracy of the radar operator. Output from the radar operator was compared to real radar observations from the Weather Surveillance Radar–1988 Doppler (WSR-88D) radar. The results showed that the radar forward operator generated realistic polarimetric signatures. The distribution of polarimetric variables agreed well with the hydrometer properties produced by different microphysics schemes. Similar to the observed polarimetric signatures, radar operator output showed ZDR and KDP columns from low-to-mid troposphere, reflecting the large amount of rain within strong updrafts. The Thompson scheme produced a better simulation for the hail storm with a ZDR hole to indicate the existence of graupel in the low troposphere.  more » « less
Award ID(s):
1746119
NSF-PAR ID:
10389641
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Atmosphere
Volume:
13
Issue:
5
ISSN:
2073-4433
Page Range / eLocation ID:
645
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. In this study, a polarimetric radar forward model operator was developed for the Weather Research and Forecasting (WRF) model that was based on a scattering algorithm using the T-matrix methodology. Three microphysics schemes—Thompson, Morrison 2-moment, and Milbrandt-Yau 2-moment—were supported in the operator. This radar forward operator used the microphysics, thermodynamic, and wind fields from WRF model forecasts to compute horizontal reflectivity, radial velocity, and polarimetric variables including differential reflectivity (ZDR) and specific differential phase (KDP) for S-band radar. A case study with severe convective storms was used to examine the accuracy of the radar operator. Output from the radar operator was compared to real radar observations from the Weather Surveillance Radar–1988 Doppler (WSR-88D) radar. The results showed that the radar forward operator generated realistic polarimetric signatures. The distribution of polarimetric variables agreed well with the hydrometer properties produced by different microphysics schemes. Similar to the observed polarimetric signatures, radar operator output showed ZDR and KDP columns from low-to-mid troposphere, reflecting the large amount of rain within strong updrafts. The Thompson scheme produced a better simulation for the hail storm with a ZDR hole to indicate the existence of graupel in the low troposphere. 
    more » « less
  2. null (Ed.)
    A multi-radar analysis of the 20 May 2013 Moore, Oklahoma, U.S. supercell is presented using three Weather Surveillance Radars 1988 Doppler (WSR-88Ds) and PX-1000, a rapid-scan, polarimetric, X-band radar, with a focus on the period between 1930 and 2008 UTC, encompassing supercell maturation through rapid tornado intensification. Owing to the 20-s temporal resolution of PX-1000, a detailed radar analysis of the hook echo is performed on (1) the microphysical characteristics through a hydrometeor classification algorithm (HCA)—inter-compared between X- and S-band for performance evaluation—including a hail and debris class and (2) kinematic properties of the low-level mesocyclone (LLM) assessed through ΔVr analyses. Four transient intensifications in ΔVr prior to tornadogenesis are documented and found to be associated with two prevalent internal rear-flank downdraft (RFD) momentum surges, the latter surge coincident with tornadogenesis. The momentum surges are marked by a rapidly advancing reflectivity (ZH) gradient traversing around the LLM, descending reflectivity cores (DRCs), a drop in differential reflectivity (ZDR) due to the advection of smaller drops into the hook echo, a decrease in correlation coefficient (ρhv), and the detection of debris from the HCA. Additionally, volumetric analyses of ZDR and specific differential phase (KDP) signatures show general diffusivity of the ZDR arc even after tornadogenesis in contrast with explosive deepening of the KDP foot downshear of the updraft. Similarly, while the vertical extent of the ZDR and KDP columns decrease leading up to tornadogenesis, the phasing of these signatures are offset after tornadogenesis, with the ZDR column deepening the lagging of KDP. 
    more » « less
  3. Abstract

    A variational retrieval of rain microphysics from polarimetric radar data (PRD) has been developed through the use of S-band parameterized polarimetric observation operators. Polarimetric observations allow for the optimal retrieval of cloud and precipitation microphysics for weather quantification and data assimilation for convective-scale numerical weather prediction (NWP) by linking PRD to physical parameters. Rain polarimetric observation operators for reflectivity ZH, differential reflectivity ZDR, and specific differential phase KDP were derived for S-band PRD using T-matrix scattering amplitudes. These observation operators link the PRD to the physical parameters of water content W and mass-/volume-weighted diameter Dm for rain, which can be used to calculate other microphysical information. The S-band observation operators were tested using a 1D variational retrieval that uses the (nonlinear) Gauss–Newton method to iteratively minimize the cost function to find an optimal estimate of Dm and W separately for each azimuth of radar data, which can be applied to a plan position indicator (PPI) radar scan (i.e., a single elevation). Experiments on two-dimensional video disdrometer (2DVD) data demonstrated the advantages of including ΦDP observations and using the nonlinear solution rather than the (linear) optimal interpolation (OI) solution. PRD collected by the Norman, Oklahoma (KOUN) WSR-88D on 15 June 2011 were used to successfully test the retrieval method on radar data. The successful variational retrieval from the 2DVD and the radar data demonstrate the utility of the proposed method.

     
    more » « less
  4. Abstract. Radar dual-wavelength ratio (DWR) measurements from the Stony Brook RadarObservatory Ka-band scanning polarimetric radar (KASPR, 35 GHz), a W-bandprofiling radar (94 GHz), and a next-generation K-band (24 GHz) micro rainradar (MRRPro) were exploited for ice particle identification using triple-frequency approaches. The results indicated that two of the radarfrequencies (K and Ka band) are not sufficiently separated; thus, thetriple-frequency radar approaches had limited success. On the other hand, ajoint analysis of DWR, mean Doppler velocity (MDV), andpolarimetric radar variables indicated potential in identifying ice particletypes and distinguishing among different ice growth processes and even inrevealing additional microphysical details. We investigated all DWR pairs in conjunction with MDV from the KASPRprofiling measurements and differential reflectivity (ZDR) and specificdifferential phase (KDP) from the KASPR quasi-vertical profiles. TheDWR-versus-MDV diagrams coupled with the polarimetric observables exhibiteddistinct separations of particle populations attributed to different rimedegrees and particle growth processes. In fallstreaks, the 35–94 GHz DWRpair increased with the magnitude of MDV corresponding to the scatteringcalculations for aggregates with lower degrees of riming. The DWR valuesfurther increased at lower altitudes while ZDR slightly decreased,indicating further aggregation. Particle populations with higher rimedegrees had a similar increase in DWR but a 1–1.5 m s−1 largermagnitude of MDV and rapid decreases in KDP and ZDR. The analysisalso depicted the early stage of riming where ZDR increased with theMDV magnitude collocated with small increases in DWR. This approach willimprove quantitative estimations of snow amount and microphysical quantitiessuch as rime mass fraction. The study suggests that triple-frequencymeasurements are not always necessary for in-depth ice microphysical studiesand that dual-frequency polarimetric and Doppler measurements cansuccessfully be used to gain insights into ice hydrometeor microphysics. 
    more » « less
  5. This study evaluates moist physics in the Weather Research and Forecasting (WRF) Model using observations collected during the Olympic Mountains Experiment (OLYMPEX) field campaign by the Global Precipitation Measurement (GPM) satellite, including data from the GPM Microwave Imager (GMI) and Dual-Frequency Precipitation Radar (DPR) instruments. Even though WRF using Thompson et al. microphysics was able to realistically simulate water vapor concentrations approaching the barrier, there was underprediction of cloud water content and rain rates offshore and over western slopes of terrain. We showed that underprediction of rain rate occurred when cloud water was underpredicted, establishing a connection between cloud water and rain-rate deficits. Evaluations of vertical hydrometeor mixing ratio profiles indicated that WRF produced too little cloud water and rainwater content, particularly below 2.5 km, with excessive snow above this altitude. Simulated mixing ratio profiles were less influenced by coastal proximity or midlatitude storm sector than were GMI profiles. Evaluations of different synoptic storm sectors suggested that postfrontal storm sectors were simulated most realistically, while warm sectors had the largest errors. DPR observations confirm the underprediction of rain rates noted by GMI, with no dependence on whether rain occurs over land or water. Finally, WRF underpredicted radar reflectivity below 2 km and overpredicted above 2 km, consistent with GMI vertical mixing ratio profiles.

     
    more » « less