Steven Rutledge, Stephen Nesbitt, Robert Cifelli, and Timothy Lang Department of Atmospheric Science Colorado State University

Transcription

1 Report and Recommendations of the Global Precipitation Mission (GPM) Ground Validation (GV) Front Range Pilot Project Steven Rutledge, Stephen Nesbitt, Robert Cifelli, and Timothy Lang Department of Atmospheric Science Colorado State University Brooks Martner*, Sergey Matrosov*, and David Kingsmill* Environmental Technology Laboratories National Oceanic and Atmospheric Administration Kenneth Gage and Christopher Williams* Aeronomy Laboratory National Oceanic and Atmospheric Administration V. Bringi and V. Chandrasekar Department of Electrical and Computer Engineering Colorado State University Patrick Kennedy Colorado State University CHILL National Radar Facility Final Version 9 January 2005 Acknowledgements: This work was supported by NASA GPM-GV supplemental grant NNG04GF32A under the direction of Matthew Schwaller. *Also affiliated with the Cooperative Institute for Research in Environmental Studies University of Colorado

2 NASA GPM GV Front Range Pilot Project Final Report Page 2 TABLE OF CONTENTS Section 1: EXECUTIVE SUMARY AND PRELIMINARY RECOMMENDATIONS...3 Section 2: Introduction and Project Overview Experimental design...6 Section 3. Ground Validation Instrument Candidates Background Outline of Instrument Advantages and Limitations Specific Instrument Issues Addressed in the Pilot Project...17 Section 4. A Template for Uncertainty Analysis for GPM GV...19 Section 5. Addressing our GV Objectives S- and X- Band Polarimetric Radar Comparisons Profilers, Disdrometers, and Gauges for Ground Validation Polarimetric Retrievals of DSD Parameters by the 2-D Video Disdrometer (2DVD) Hydrometeor Identification by Polarimetric Radar...38 APPENDIX 1. Characteristics of the Scanning Radars, Profilers, and Disdrometers Used in the Pilot Project...41 APPENDIX 2: Comparisons of Algorithms for Calculating K dp and Corrections to Z eh and Z dr Measurements for the Effects of Attenuation at X-band...45 A2.1. Comparisons of Algorithms for Calculations of K dp...45 A2.2. Comparisons of Algorithms for Correcting Z eh and Z dr for Attenuation Effects...46 APPENDIX 3: Case study results...48 A June Case Study...48 A June 2004 Case Study...53 A July Case Study...59 REFERENCES...65

3 NASA GPM GV Front Range Pilot Project Final Report Page 3 Section 1: EXECUTIVE SUMARY AND PRELIMINARY RECOMMENDATIONS We are pleased to submit this report describing activities and findings from the Front Range GPM Pilot Project conducted during the summer of The purpose of the Pilot was to bring together a team of scientists and instruments to carry out preliminary design and evaluations that would aid planning for the GPM GV mid-latitude Supersite concept. The Pilot was a successful collaboration between scientists at Colorado State University (Departments of Atmospheric Science and Electrical and Computer Engineering), and NOAA s Aeronomy Laboratory and Environmental Technology Laboratories. Approximately 20 faculty, staff and graduate students participated in the project. Some key finding of the Pilot are: X-Band polarimetric radar is capable of estimating rain rates as low as ~2 mm hr -1 ; polarimetric-based rain estimates at S-band begin at 5-7 mm hr -1. At rain rates less than this, S-band radars must utilize the conventional Z-R estimator, which can be subject to significant error due to variability in the drop size distribution and calibration uncertainties. Under high rain rate conditions (rates of several tens of mm hr -1 ) X-Band signals are subject to severe attenuation losses. While X-Band attenuation procedures have been developed, complete signal extinction can occur, rendering precipitation estimation impossible. Therefore S-Band and X-Band form a complimentary combination that can apply the more accurate polarimetric rain estimators to a wider range of rainfall rates than either could handle alone. We also offer the recommendation that two vertically pointing profilers operating at different frequencies be collocated at the Supersites. One UHF profiler operating at 449 MHz would be used to estimate the vertical air motion and the other profiler operating near 2835 MHz (S-Band) would be used to estimate the Doppler fall velocity of the raindrops. Algorithms can be developed to estimate the air motion from the 449 MHz profiler spectra and used with the 2835 MHz profiler spectra to estimate the DSD from near the surface to just below the freezing level. 2-D Video and Joss-Waldvogel momentum disdrometers proved to be invaluable in characterizing calibration and algorithm uncertainty in all radar, profiler measurements made in the Pilot Project. Based on the data and results produced by the GPM Front Range Pilot Project, we preliminarily recommend the following instrumentation suite to accomplish the GPM GV task at the GPM continental Supersite 1 : Dual-wavelength (S- and X-Band) scanning dual-polarization Doppler radar. As discussed in section 5 and the Appendices, this system will be the core instrument for retrieving the 4-D evolution of precipitation rates (using polarimetric measurements at rain rates greater than about 2 mm hr -1 and reflectivity alone for lighter rain) and polarimetric hydrometeor identification methodology. The instrument will collect volume scans at 5-10 minute time resolution over the Supersite Multidimensional Observing Volume (MOV), providing the 3-D hydrometeor volume inputs required the Satellite Simulator Model (SSM). 1 We are assuming that the site will be deployed at the Department of Energy ARM-CART site near Lamont, Oklahoma. This is because facilities such as regular radiosonde launches, ceilometers, microwave water vapor path sensors, etc., present at the ARM-CART site will also be required for atmospheric state variable input into microphysical retrievals made by the Multidimensional Observing Volume (MOV) and as input into the the Satellite Simulator Model (SSM).

4 NASA GPM GV Front Range Pilot Project Final Report Page 4 To reduce uncertainty related to geometrically matching beams from two radars, it is proposed that the radar to be developed for GPM-GV have a single antenna for both wavelengths, which will be facilitated by a dual-wavelength dual-polarization feed horn currently in development by VertexRSI, Inc. While the beam width dimensions will most likely not be identical between the two radars (with the X-Band data having higher resolution than at S-Band), the pointing angles of the two frequencies will be matched such that direct comparisons (varying only in beam width) may be made between the S- and X- Band measurements. A system of this type will also provide a lower-cost alternative to a dual-antenna system or utilizing two radar systems for the radar site. Preliminary procurement cost estimate: $3-5 M. Dual-wavelength (S-Band and 449 MHz) vertically pointing Doppler radar system. The dual-frequency profiler system will provide the crucial microphysical measurement link between the scanning radars and the detailed, in situ microphysical measurements made on the ground. This system will provide measurements of radar reflectivity, vertical velocity spectra of hydrometeors (at S-Band and 449 MHz), and clear air velocity spectra (at 449 MHz) with 1 minute time resolution. This information can be used to retrieve radar reflectivity, fall speed, and DSD parameter profiles, which can be compared both to ground instrumentation measurements of the above parameters (at the lowest gate of the profilers) and to overlapping coverage with the dual-frequency radar measurements (at overlapping gates aloft). This capability will be crucial in providing a pathway to propagate uncertainty estimates between the ground sensors and the scanning polarimetric radar microphysical estimates. The profiler system will consist of 2 separate transmitter/antenna systems, due to the fact that the S-Band (449 MHz) systems are typically have an upward pointing dish (phased-array) antenna. These systems will be located side-by-side at some range from the radar about km away from the scanning polarimetric radar, calibrated with collocated surface disdrometer(s), and require relatively little maintenance compared with a scanning radar. Preliminary procurement cost estimate: $ K. Network of Surface disdrometers. This project has shown the wide applicability of the surface disdrometer as a rain rate and DSD measurement tool, as well as a calibration tool for profiling and scanning radars. Based on the scientific requirements of the Supersite, 2DVD disdrometers have demonstrated superiority relative to impact-type Joss-Waldvogel disdrometers. Outlined here are results and recommendations based on our findings, realizing that some combination in number of both types will likely be deployed due to budgetary issues. 2D video disdrometers (2DVDs). The 2DVD can provide direct in situ measurement of hydrometeor habit, number concentration, size, axis ratio, and fall speed at the surface. These measurements can be used to evaluate uncertainties in radar and profiler microphysical algorithm assumptions and retrievals. These measurements can also be used to directly simulate (and thus calibrate) profiler and scanning radar measurements of radar reflectivity Z e, as well as scanning polarimetric radar measurements of differential reflectivity Z dr and specific differential phase

5 NASA GPM GV Front Range Pilot Project Final Report Page 5 K dp. It is envisioned that several of these instruments will be placed at various ranges from the radar to evaluate uncertainties related to range effects in the scanning radar retrievals (which will be required to address uncertainties within the MOV), as well as at the profiler site as a validation and calibration standard for profiler measurements and retrievals. Calibration of the low-profile 2DVD used in this project has been demonstrated to be stable and straightforward, when necessary. Preliminary procurement cost estimate: 10 x $80 K each = $800 K Impact-type Joss-Waldvogel disdrometers (JWDs). This set of instrumentation can provide in situ measurements of the DSD and rain rate at the surface. The JWD incurs well-known DSD measurement errors due because of their impact-based technology, cannot measure hydrometeor shape, habit, or fall speed (these must be assumed), and are not able to be calibrated on site beyond gross comparisons with collocated gauges. Thus, these instruments provide less complete information than 2DVDs. However, they generally require less human attention in the field and can provide a lower cost alternative or supplement. Preliminary procurement cost estimate: 20 x $20 K each = $400 K Network of tipping bucket rain gauges. The Front Range Pilot Project had limited gauge resources available, such that their full utility was not demonstrated beyond point comparisons with disdrometer, profiler, and polarimetric radar rain rate estimates at the Platteville and BAO instrument sites. It is anticipated that a dense, broad rain gauge network will be required to evaluate the uncertainty in scanning polarimetric radar rainfall estimates over the MOV (10 x 10 km area). Assuming a triple redundant gauge network with 1-km resolution, this would require 300 tipping bucket gauges. However, the details of gauge network design are left to future work. Preliminary procurement cost estimate: 300 x $0.5 K each = $150 K.

6 NASA GPM GV Front Range Pilot Project Final Report Page 6 Section 2: Introduction and Project Overview The Pilot project s specific aims included the following: 1. Dual-wavelength radar DSD and rain rate estimate intercomparison, validation, and error characterization. Demonstration of dual-wavelength polarimetric radar network for creating rain estimates and documenting associated errors. The project used the S- Band CSU-CHILL radar at Greeley, CO and NOAA-ETL s X-Band radar at Erie, CO. The X-Band radar s improved phase sensitivity in light rain was evaluated against the S- Band radar s performance in such conditions. The S-Band s known insensitivity to attenuation in heavy rain was used to evaluate the X-Band radar s ability to correct for attenuation using its specific differential phase measurements, and thus diagnose heavy rain. 2. Profiler demonstration in the supersite concept. Selection of UHF profiler frequencies that will best complement S-Band profiler measurements at a midlatitude site and allow for the most accurate retrieval of drop size distribution characteristics as a further goal of the Pilot. A ancillary goal was to perform quantitative comparisons of drop size distribution (DSD) characteristics between the profilers and scanning radars in order to evaluate assumptions in the scanning radar retrieval technique (e.g., equilibrium drop shape relationship) as well as spatial variability of the DSD. 3. Rain rate and drop size distribution characterization in the context of supersite observations, rainfall regimes. The Pilot sought to demonstrate the complementary role played by rain gauges and surface disdrometers (both 2DVD and JWD types) in determining the error characteristics of multi-frequency profiler DSD estimates and dualfrequency radar DSD and rain estimates. 2.1 Experimental design Figure 2.1 shows the instrument locations for the Pilot. The CSU-CHILL S-Band polarimetric radar was located at its home base near Greeley, CO (marked CHILL in the figure). The NOAA- ETL X-Band polarimentric radar was deployed near Erie, CO (marked Erie-1). Both radars scanned narrow azimuth sectors over the continuously operating ground measurement sites, where the profilers, disdrometers, rain gauges, and surface meteorological stations were situated. The ground measurement sites were located at the NOAA-University of Colorado Platteville Atmospheric Observatory near Platteville, CO, and at the NOAA Boulder Atmospheric Observatory (BAO) near Erie, CO. Ground Measurement Sites Both the BAO and Platteville sites had profilers operating at 3 frequencies: one at S-Band to measure precipitation velocity spectra (and estimate DSD), and two in the UHF band: 915 and and 449 MHz. The UHF profilers were evaluated to determine their ability: 1) to resolve both the clear air and precipitation components of the radial velocity spectra in different precipitation

7 NASA GPM GV Front Range Pilot Project Final Report Page 7 environments (e.g., light vs. heavy rain); and 2) in concert with the S-Band profiler data, provide dual-wavelength DSD parameter retrievals with the smallest errors. Each site was instrumented with at least one JWD impact-type disdrometer, tipping-bucket rain gauge, and surface meteorological station. CSU s 2DVD was also installed at the BAO site until the end of June, at which time it was moved to Platteville. Radar scanning The X-Band radar nominally scanned 2 low-level PPIs and 2 RHIs covering the Platteville and BAO sites. Gate to gate resolution was set to either 150 m or meters, giving a maximum range of 38.8 or 28.8 km, respectively. The X-Band scan cycle took place over a period of one to two minutes. CHILL scanned a sector of approximately 40 degrees azimuth over the entire X- Band scan sector (and the instrumented sites). RHIs were also obtained over the instrumented sites. CHILL had two scanning modes depending on the rainfall regime: (1) a high-resolution 75 m range resolution low level scan mode, and (2) a standard resolution 150 m range resolution volume scan mode with a 6-8 minute scan repeat cycle (allowing storms vertical reflectivity and microphysical structure to be scanned by CHILL). Project timeline The overall data collection period for the Front Range Pilot Project was 15 May through June 30, All instruments operated during this time period with very little if any downtime. Data collection continued until 23 July on a target of opportunity basis. Several nice cases were obtained in this extension period. Table 2.1 shows the cases collected during the formal data collection period of the project, as well as those collected after June 30, 2004.

8 NASA GPM GV Front Range Pilot Project Final Report Page 8 Table 2.1. GPM-FRPP 2004 Case List. Key is located at the top of the table.

9 NASA GPM GV Front Range Pilot Project Final Report Page 9 Figure 2.1. Map of the study areas showing the locations of experiment sites (bold plus signs) with instrument list, radar scan coverage areas (white lines). terrain (shaded), principal highways (grey lines), and county boundaries (thick opaque lines) also shown. KFTG is the Denver NEXRAD radar site.

10 NASA GPM GV Front Range Pilot Project Final Report Page 10 Section 3. Ground Validation Instrument Candidates 3.1. Background A preliminary conceptual plan for GPM s continental ground validation (GV) Supersite, shown in Figure 3.1, was developed prior to the Front Range Pilot Project, based on experience obtained during the Tropical Rainfall Measurement Mission (TRMM) and on requirements of the algorithm and data assimilation communities to define error statistics of the GV measurements. Basic GV requirements include obtaining accurate, high-resolution measurements of precipitation rate and its quantitative error statistics, as well as 3-D fields of hydrometeor type, independent of the satellites, within logistical and financial constraints. No one instrument can measure these quantities alone, thus a suite of complementary instrumentation is required. Furthermore, we recognize that one continental and one maritime GV Supersite will not adequately observe the full variety of precipitation characteristics that the satellite may observe. Hence, our vision is the two Supersites will constitute a ground-based core of instruments operating unattended around the clock, 365 days a year, for the several-year lifetime of the satellite program, providing a long-term dataset that is rigorously similar in form throughout. In addition to these routine data, these core Supersites will occasionally host expanded observation campaigns to target specific questions posed by the observational or algorithm community (i.e. to validate a specific algorithm assumption). Meanwhile, auxiliary or constellation GV sites may also operate in other precipitation regimes (e.g. snow, coastal rain, orographic precipitation), on a more limited basis in terms of instrumentation and duration, to address algorithm issues in weather that is inadequately observed at the core Supersites. In this report, we address some of the primary candidate instruments from the preliminary conceptual plan for the continental GV Supersite. The Front Range Pilot Project was specifically charged with investigating the GV utility of polarimetric scanning weather radars and Doppler profiling radars. These are big-ticket items with long lead times that justify an early analysis effort. To accomplish this work, the use of rain gauges and raindrop disdrometers was also regarded to be essential. Therefore, the scope of this report is limited to discussing these instruments, although the GV sites may also include others. We visualize a continental GV Supersite that covers overlapping measurement scales and capabilities with a dense network of rain gauges, several raindrop disdrometers, one or more Doppler dual-frequency profiling radars, and one dual-frequency polarimetric, scanning, precipitation radar. As background, basic information about the useful attributes and known problems associated with these instruments is outlined in this section. Conventional (non-polarimetric) radar is included in the outline for comparison purposes. This outline is intended only for use in preliminary guidance and is not comprehensive.

11 NASA GPM GV Front Range Pilot Project Final Report Page 11 Figure 3.1. A preliminary conceptual vision of the continental GV Supersite.

12 NASA GPM GV Front Range Pilot Project Final Report Page Outline of Instrument Advantages and Limitations The following tables outline the features of each instrument platform used in the Pilot Project, their capabilities, features that are attractive for their deployment at a GPM-GV site, as well as challenges and limitations that must be taken into consideration in their deployment. R = rainfall intensity, RA = rainfall accumulation, Z = radar reflectivity, DSD = drop size distribution Precipitation Gauges Physical Basis Measures time-resolved accumulation of precipitation mass in situ at the ground Attractive GV Attributes Direct measure of R and RA at a single point on the ground. Individual instrument is relatively inexpensive (~$500 for a very good tipping bucket gauge) Good temporal resolution (1 minute is typical) Instrument calibration is relatively simple Commercially available off-the-shelf Limitations & Difficulties Provides no area information, unless deployed in networks Provides no information on precipitation aloft Provides no information on DSD Large, dense networks for area coverage are expensive and maintenance is manpower-intensive Siting must avoid rain shadow of trees and buildings Wind-induced under-catch causes underestimates of R and RA Snowfall measurement requires more sophisticated gauges, shields, and siting Many new gauges that are more sophisticated than the traditional tipping bucket are available commercially, but have much shorter track records

13 NASA GPM GV Front Range Pilot Project Final Report Page 13 Raindrop Disdrometers Physical Basis Counts individual raindrops and determines their sizes in situ at the ground Attractive GV Attributes Provides measurements of DSD, from which many other rainfall parameters can be computed, including radar reflectivity, Z Direct measure of R and RA at a single point on the ground Good temporal resolution (1 minute is typical) 2DVD units also measure hydrometeor fall speeds and shapes A valuable tool for checking scanning radar and profiling radar estimates of rain parameters Commercially available off the shelf Limitations & Difficulties Provides no area information, unless deployed in networks Expensive (~$20K for momentum disdrometers; ~$80K for 2DVD disdrometers), thus networks of disdrometers must be relatively sparse compared to gauges Momentum disdrometers must assume drops are falling at terminal speed Momentum disdrometers have no feasible in-the-field calibration, although comparisons of RA with collocated gauges are useful indicators of calibration worthiness 2DVDs require more frequent human attention than momentum units

14 NASA GPM GV Front Range Pilot Project Final Report Page 14 Doppler Precipitation Profilers Physical Basis Reflectivity and spectrum of vertical speeds of hydrometeors aloft are remotely sensed by upward pointing Doppler radar and converted to vertical profiles of information on drop size distributions Attractive GV Attributes Provides profiles of precipitation conditions aloft at a single point Continuous, unattended operation Intermediate scale link between gauges and scanning radar sample volumes and the even larger sample volume of satellite radars and radiometers Provides information on DSDs aloft using Doppler techniques Provides continuous information about height of the melting layer Provides derived information about updraft/downdraft air motions aloft Commercially available off-the-shelf for most frequencies Valuable tool for checking scanning radar estimates of precipitation parameters Limitations & Difficulties Provides no area information, unless deployed in networks Expensive (~$100K for 2835-MHz; ~$200K for 915-MHz) DSD s can be derived confidently only if clear-air and hydrometeor Doppler spectra peaks can be clearly separated; this is often difficult and may require use of two frequencies (eg., 2835 and 449 MHz) Calibration generally requires a collocated raindrop disdrometer

15 NASA GPM GV Front Range Pilot Project Final Report Page 15 Conventional Scanning Weather Radars Physical Basis Reflectivity of hydrometeors is sensed remotely, mapped in 3D, and converted to estimates of precipitation intensity using empirical Z-R relations Attractive GV Attributes A single radar provides contiguous patterns of R and RA estimates over the entire satellite swatch width, typically with resolution ~1 km Provides 3-D information on precipitation across most of the satellite swath Operational networks are currently in place in several countries (eg. NEXRAD) Doppler data provides information about storm airflow patterns from which storm dynamical properties can be deduced Limitations & Difficulties Expensive (~$0.5-$1.0M for basic S-band system; totally autonomous systems are more expensive, NEXRAD units cost ~$10M) R and RA are estimated from Z, rather than directly measured Z-based estimates of R and RA are fraught with error sources associated with: o variable DSDs o improperly calibrated radar hardware o partial beam blockage (eg by terrain for low scans) o attenuation by rain (wavelengths < 10 cm) o inappropriate use of rain Z-R relations in regions of snow, melting snow, or hail Calibration is often difficult Provides no information on DSDs Difficult to distinguish rain from snow or hail Estimates of snowfall are highly inaccurate

16 NASA GPM GV Front Range Pilot Project Final Report Page 16 Dual-Polarimetric Scanning Radars Physical Basis Reflectivity, differential reflectivity, differential phase and other parameters returned from hydrometeors are measured remotely in two orthogonal polarization states, mapped in 3D, and converted to estimates of precipitation intensity, particle type, and DSD parameters using modeled relations involving particle shapes. Attractive GV Attributes A single radar provides contiguous patterns of R and RA estimates and several other precipitation parameters, such as median drop diameter, over the entire satellite swatch width, with good resolution Provides 3-D information on precipitation across much of the satellite swath Polarimetric capability allows reflectivity to be adjusted for attenuation effects at shorter wavelengths Doppler data provides information about storm airflow patterns from which storm dynamical properties can be deduced Provides information on DSDs using polarimetric techniques (complements DSD retrievals from profiling radars using Doppler techniques) Polarimetric measurables provide more accurate estimates of R and RA by avoiding many of the reflectivity-related rain estimation problems of conventional radar Provides information for identifying particle type (raindrops, snowflakes, hail, graupel, insects, etc.) Dual-wavelength retrievals of DSD information are possible (as will be done from the GPM-DPR) Polarimetric measurables are useful for checking hardware calibrations Limitations & Difficulties Expensive (~$3M for S-band, $1M for X-band) Not commercially available as off-the-shelf units Operational networks are not now available (but NEXRAD might be converted to polarimetric capability in the next decade) Two frequencies may be required for good polarimetric estimates of R and RA over the range of heavy to light rainfall rates Estimates of snow can benefit from the additional polarimetic information on particle shapes and from the use of two wavelengths, but these methods are still very experimental

17 NASA GPM GV Front Range Pilot Project Final Report Page Specific Instrument Issues Addressed in the Pilot Project Subsequent sections of the report focus in much greater detail on the particular radars and disdrometers used in the Pilot Project. See Appendix 1 for exact instrument characteristics of the instruments used in the Pilot Project. The emphasis of the investigations was to answer the following questions: 1. Does the addition of an X-band frequency to an S-band radar allow weaker rainfall rates to be accurately measured by polarimetric techniques? 2. What is the best frequency or combination of frequencies to be used by a Doppler profiling radar for retrieving raindrop drop size distributions (DSD) aloft? CSU and NOAA instruments were used in the Pilot Project because they were available for costeffectively testing these questions in Colorado. It is expected that the Pilot Project s findings will provide helpful guidance for designing new, even better, GPM radars, tailored for long-term dedicated use at the continental GV Supersite. Although conventional reflectivity-based radar estimates of rainfall rate have been used for half a century, the accuracies of the estimates leave much to be desired. As shown in the outline list of Section 3b, Z-R estimates of rainfall suffer from many factors, such as DSD variability and calibration uncertainties, that commonly degrade the results. New polarimetric radars avoid many of these reflectivity-based problems and offer a large number of new measurement parameters that are related to hydrometeor shapes (Zrnic and Ryzhkov 1996, Bringi and Chandrasekar 2001). These parameters are useful for identifying hydrometeor types (raindrops, snowflakes, hailstones, etc.) and for providing information about DSD parameters, such as median drop diameter, that cannot be obtained with conventional radar, as well as for improving estimates of rain intensity and accumulation. The ability to measure these parameters with polarimetric radar is crucial for fulfillment of the GV tasks of the Supersite. Most polarimetric radar research has been conducted at S-band (10-cm wavelength), such as with CSU s CHILL radar. Studies have shown S-band polarimetric radar using specific differential phase (K dp ) measurements can yield more accurate rainfall estimates than those available from reflectivity alone (e.g., Aydin et al. 1995, Ryzhkov and Zrnic 1996). However, this has been mainly limited to moderate to heavy rainfall situations, because K dp is too small to measure at S- band in lighter rain dominated by small, more nearly spherical drops. Importantly, the K dp sensitivity varies inversely with wavelength. Thus, shorter wavelength polarimetric radars, such as the NOAA/ETL X-band (3 cm), can measure the K dp parameter and use it to estimate rain rates in lighter rainfall relative to S-Band systems. The feasibility of using X-band in this manner has been demonstrated (Matrosov et el. 2002, Anagnostou et al. 2004). Therefore, a combination using S-band (for heavy-moderate rain) and X-band (for moderate-light rain) radar with matched beams is an attractive possibility. A primary goal of the Pilot Project was to test this possibility, using available (but not ideal) resources, and to determine the range of lighter rainfall rates to which X-band can extend polarimetric rain estimation applicability. A large fraction of the global rainfall occurs at light rain intensities. Thus, these light rainfall regimes are important and cannot be ignored in climate studies that GPM is designed to support.

18 NASA GPM GV Front Range Pilot Project Final Report Page 18 A second prime goal of the Pilot Project was to determine which frequency or combination of frequencies would be best for the Doppler profiling radar GV work. Doppler profilers have a demonstrated capability for retrieving DSD information aloft (Williams et al. 2000a). They can provide continuous vertical profiles of rain parameters, such as median drop diameter, on scales that bridge the gap between the tiny sampling volume of a rain gauge and the gigantic sample volumes of scanning radars and satellite instruments. A difficult aspect of the profiler DSD retrievals, however, involves separating peaks of the observed Doppler spectra that are caused by the terminal fall speed of the hydrometeors (Rayleigh scattering) from the peaks associated with the updraft/downdraft speeds of the air parcels (Bragg scattering) in which the precipitation particles are embedded (Williams et al. 2000b). The profiling radar measures the sum of these two contributions. However, higher frequency profilers are more sensitive to detecting the hydrometeors and lower frequencies are more sensitive to detecting the air motions (Gage et al. 1999). DSD retrieval techniques involving single and dual-frequency profilers have been developed to attack this problem. The Pilot Project used combinations of 915, 449, and MHz from NOAA/AL and NOAA/ETL to address this issue.

19 NASA GPM GV Front Range Pilot Project Final Report Page 19 Section 4. A Template for Uncertainty Analysis for GPM GV The ability to report quantitative uncertainty estimates in routine products generated at GPM Supersites will be a strict requirement of GV for GPM (cf., Yuter et al. 2002). This requirement will provide new challenges in understanding the error characteristics of the GV instrumentation as a function of regime (regime is defined here as a mode of precipitation type in a vertical column that defines a particular known range of microphysical and accompanying error characteristics within the spaceborne precipitation retrieval algorithms). Measurements made at the GV site must be able to not only quantify measurands and uncertainties related to observed vertical profiles of precipitation, but also be able to place those quantities within context of meteorological regime (i.e. weather and storm conditions which lead to similar storm types and thus similar error characteristics) such that error statistics determined at the GV site may be transferred to storm types, with presumably similar error characteristics, observed elsewhere by the satellites. Thus the focus of GV will not only revolve around determining measurands and uncertainties routinely required by the GPM algorithm developers, but also making the appropriate measurements to place those quantities in the context of regimes that are observed globally by the satellite. Consider a simple error model, varying as a function of regime (regimes represented by x). Recall we have defined the regime such that error characteristics are stationary within a predefined tolerance, defined by the algorithm developers. Note that the regime observed at a particular location varies as a function of time, so regimes in reality vary in both space and time. The example given here is a general form for estimating the true rain rate R(x) over a particular area using dual-frequency polarimetric radar, however it may be applied to other variables of interest. The GV ensemble measurements will provide their best estimate of rain rate R ˆ (x), plus a total error term (x). In mathematical form: R ˆ (x) = R(x) + (x). Error sources will vary based on the measurand of interest. The following sections will define, and discuss contributions to the total error by measurement error m (x), parameterization error p (x), and sampling error s (x) as such: (x) = m (x) + p (x) + s (x). Measurement error Within the GVS, measurement error will encompass both the instrumental and statistical uncertainties related to making a finite sample of measurements of the GV parameters. This will include calibration uncertainty (bias offsets) as well as random statistical fluctuations in the data (random perturbations). For the example of rain rate estimation with polarimetric radar data, Z- and Z dr -based rainfall estimators will be most affected by radar calibration uncertainty. This calibration uncertainty may be addressed with techniques ranging from polarimetric selfconsistency checks (Gorgucci et al. 1996, Carey and Rutledge 2000) to reflectivity, differential reflectivity, and specific differential phase comparisons with other instrumentation such as profilers, impact and video disdrometers, and rain gauges (Gage et al. 2000, Schuur et al. 2001). However, all measured input variables will be subject to statistical uncertainties. For example, statistical uncertainty in K dp estimated rain rates may be calculated as a function of the measured

20 NASA GPM GV Front Range Pilot Project Final Report Page 20 standard deviation of total differential phase ( dp ) with the following (Gorgucci et al 1999, Bringi and Chandrasekar 2001): ( m ) R = b ( 1/b dp ) 3 c, L [N (1/ N)] R where b, and c are the coefficients in the R-K dp relationship, L is the path length over which K dp is calculated, and N is the number of points within that path length. Variability in measurement error can be reduced through spatial averaging of the data, however, this cannot be done without introducing representativeness (sampling) uncertainty when comparing rain rates with rain rates measured by instruments of different sample volumes such as disdrometers or rain gauges (Tustison et al. 2001). Moreover, spatial averaging to reduce errors in K dp can introduce errors due to non-uniformity along the path length (Gorgucci et al. 2000). Parameterization Error Parameterization error revolves around uncertainties within the inversion process of converting measured quantities (i.e. radar reflectivity Z, differential reflectivity Z dr, specific differential phase K dp in the case of a polarimetric radar) into derived quantities (i.e. rain rate R, massweighted mean drop diameter D m, normalized intercept parameter N w, etc.). In other words, it is a quantitative assessment of uncertainties introduced by the model and assumptions selected to derive a particular quantity. For most quantities measured at the GVS, this uncertainty will be a strong function of the regime examined, as precipitation characteristics will directly impact the validity of the assumptions used. For example, DSD assumptions in the radar algorithm will be a strong function of regime, and measurements with known uncertainty can be placed in the context of these regimes. For the example of a polarimetric radar rainfall algorithm, parameterization error will arise when the natural variability of the parametric form of the drop size distribution or size-shape relationship differs from the assumed form. See Figure 4.1 for an example of how DSD measurements may be used to quantify rainfall algorithm parameterizations, in this case for drop axis ratio. Within the GVS, these error sources may be quantified through comparisons of the polarimetric radar assumptions and DSD measurements provided by profilers, impact and video disdrometers, or targeted aircraft observations. Representativeness error Tustison et al. (2001) and Tustison et al. (2003) are pioneering works quantifying representativeness (or sampling) uncertainties related to comparing precipitation measurements and estimates with differing sampling resolutions. As such, the Tustison et al. studies provide a template for defining sampling uncertainty at the GV site. Algorithms of this type must be generalized for the GVS to estimate this source of uncertainty in comparing measurements ranging in scale from the gauge and disdrometer sampling area scales (effective sampling volumes of O ~ m 3 ), to profiler and scanning radar sampling areas

21 NASA GPM GV Front Range Pilot Project Final Report Page 21 (O ~ m 3 ), to the satellite pixel scale (O ~ m 3 ). In addition, the impact of differences in temporal sampling among GVS instruments and satellite measurands must be quantified (Kummerow et al addresses this, albeit indirectly). The design of the GVS, specifically matching the beams of the dual-frequency system and the strategic placement of surface gauges, disdrometers, and profilers) must take into account the ability to quantify representativeness error both in a spatial and temporal sense. Synthesis Measurements and uncertainties from the various GVS instruments provide the basis for construction of a GVS error model. This model will be constructed using the methodology discussed above for the three error sources. Figure 4.2 shows a conceptual diagram, indicating how this would be addressed for a polarimetric radar rain rate algorithm. The suite of GVS and ancillary observations will be used for the following two tasks: (1) generate regime-dependent error statistics or lookup table of regime-observation covariance statistics by generating a catalog of regime-dependent error statistics that may be applied to radar measurements (i.e. generate regime-dependent (x) = m (x) + p (x) + s (x)), and (2) use the lookup table above and available observations to report the best estimate of the parameter of interest (i.e. R ˆ ± at the radar pixel) Task (1) provides an objective transfer standard that can be used to assess precipitation errors around the globe, outside the Supersite and constellation sites. Task (2) provides the information necessary for algorithm developers to quantitatively assess the assumptions in their radiative transfer models and determine where the largest amount of uncertainty remains.

22 NASA GPM GV Front Range Pilot Project Final Report Page 22 Figure 4.1. Figure 7.2 from Bringi and Chandrasekar (2001) showing experimental variability in axis ratio relationships relative to two functional fits to the data. The 2DVD, when deployed to the Supersite(s) will allow routine examination of this parameter, thus allowing characterization in uncertainty related to this assumption. Figure 4.2. Diagram summarizing error analysis procedure for radar surface rain estimates.

23 NASA GPM GV Front Range Pilot Project Final Report Page 23 Section 5. Addressing our GV Objectives 5.1. S- and X- Band Polarimetric Radar Comparisons Validation of the precipitation and microphysical estimates formulated by the GPM satelliteborne sensors is critically dependent upon obtaining an accurate characterization of the hydrometeor properties within the satellite s sampling volume. Mesoscale precipitation systems typically contain atmospheric volumes that amount to many thousands of cubic kilometers. Scanning meteorological radar systems afford the most practical method of collecting threedimensional observations of these precipitation systems in a timely manner (i.e., volume scan cycle times on the order of two to five minutes). Prior research has established that dual-polarization radar techniques provide improvements in the accuracy of radar estimates of the instantaneous rain rate over traditional reflectivity-only radar technology. Differential reflectivity (Z dr ), which is the ratio of the measured horizontal and vertical radar reflectivities, provides information on the drop size distribution within the radar volume in rain-only situations, thus providing an important constraint on rain rate and microphysical retrievals. In addition, the difference between the phases of the horizontally and vertically polarized return signals has been found to be a useful indicator of the existence of the oblate (larger diameter) raindrops that are typically present at higher rain rates, especially in mixed-phase environments where reflectivity (Z) and differential reflectivity (Z dr ) rainfall estimation methods fail. For a given drop diameter, the magnitude of this differential propagation phase shift ( dp ) is inversely related to the radar wavelength. Thus, the propagation phase shift (K dp ) observed by a short wavelength radar (i.e. 3 cm or X-Band) will be greater than that observed by a long wavelength radar (i.e. 11 cm or S-Band) when both systems are viewing the same precipitation path. While X-Band systems provide more sensitive K dp retrievals, they are also subject to attenuation in their Z and thus Z dr measurements in moderate to heavy rain, often becoming completely attenuated when propagating through paths of heavy rain. Thus, the S-Band system provides rain and microphysical retrievals in this situation where K dp sensitivity is not an issue. Furthermore, the estimation of K dp from differential propagation phase in heavy rain or mixed phase precipitation can prove difficult at X-band. A major goal of the Front Range Pilot Project was the exploration of the differences in polarimetric rainfall estimation and microphysical measurement capabilities at S and X-Bands (11 and 3.2 cm wavelengths). Detailed case study results are presented in the appendices of this report. A few key findings may be summarized as follows: 1.) X-Band radar is capable of detecting useable K dp magnitudes in rain rates as low as ~2 mm hr -1 ; useable S-Band K dp data cannot be obtained at these low rain rates. (21 June 2004 event; Appendix section ) The smallest reliabily detectable S-Band dp range gradient (K dp ) is ~0.1 o km -1. This corresponds to a rain rate of approximately 5-7 mm hr -1. The results in all three cases in Appendix 2 show that X-Band Kdp estimators are thus able to perform over a significantly larger area and extend polarimetric rain estimation capabilities lower into the rain rate distribution than those at S-Band.

24 NASA GPM GV Front Range Pilot Project Final Report Page 24 3.) Under high rain rate conditions (rates of several tens of mm hr -1 ) X-Band signals are subject to severe attenuation losses. While X-Band attenuation procedures have been developed, complete signal extinction can occur. (9 June 2004 case; Appendix section 2.1; also Appendix 1b). 4.) Both S- and X-Band multiparameter radar measurements (Z, Z dr, K dp ) may be used to estimate parameters of the rain DSD, such as parameters of the gamma DSD (D m, N w ; Ulbrich 1983). Appendix 2 sections b and c show comparisons of ground-based and profiler-based D m estimates with S- and X-Band estimates of D m collected during the GPM Pilot Project. 5.) Hydrometeor Identification may be employed reliably at both S- and X-Band in order to identify appropriate hydrometeor types within the radar scanning volume. See Section 5.4 for further discussion of the utility and implementation of hydrometeor ID for polarimetric radars for GPM. In summary, the ground-based scanning radar data collected during the Front Range Pilot Project explored the utility of X and S-Band polarimetric radar measurements under a wide variety of warm season precipitation regimes. The increased dp sensitivity at X-Band allows the improved accuracy of polarimetric rain rate estimation to be realized at rain rates as low as ~2 mm hr -1. At these low rain rates, the S-Band dp signal is too noisy to be useful; rain rates must reach levels of ~5 7 mm hr -1 to obtain stable dp range profiles at S-Band. However, to overcome the attenuation losses that impact X-Band measurements, S-Band observations are required to make rain rate estimates through long (several tens of kilometers), moderate-to-heavy precipitation-filled beam paths. In view of these results, to provide accurate rain rate estimations over the widest range of precipitation intensities, GPM satellite ground validation studies should make use of polarimetric radar data collected at both X and S-bands. Serious consideration should also be given to collocating the X and S-Band radars that cover the ground validation site, since this would ensure that the collection of dp data takes place with high time synchronization along a common beam path. The use of a single antenna site would also improve cost efficiency by permitting the two radars to share the expenses associated with site leasing, installation of commercial power and data communication lines, etc. A well-chosen radar site would provide high resolution, three-dimensional scan coverage over the other ground-based sensors (disdrometers, profilers, rain gauges, etc.) that are required to collect GPM validation data Profilers, Disdrometers, and Gauges for Ground Validation With NASA support the NOAA Aeronomy Laboratory (AL) has gained considerable experience using precipitation profilers to study the structure, evolution and variability of precipitating clouds. Precipitation measurements can be made with sufficient vertical resolution to categorize precipitation in deep and shallow convective systems and in stratiform systems. A recent focus of AL research with profilers has been to provide calibration and validation in support of satellite precipitation measurement missions such as the NASA Tropical Rainfall Measuring Mission

25 NASA GPM GV Front Range Pilot Project Final Report Page 25 (TRMM). Observations obtained by AL for TRMM during TRMM Ground Validation field campaigns have provided important information on the vertical structure and temporal evolution of precipitating cloud systems (Gage et al., 2000, 2002). The TRMM profiler observations can be viewed on the AL web page ( Observations made during the field campaigns have been the subject of collaborative research with other TRMM researchers with an emphasis on the use of profilers to calibrate scanning radars used for TRMM ground validation research and the use of profilers to retrieve drop-size distributions and related precipitation parameters of interest to the TRMM Science Team. Validation of drop-size distributions used in algorithms is key to improving the retrieval of rainfall estimates from the TRMM satellite data. The profiler-based precipitation research described above also can be used to provide calibration of NEXRAD scanning radars as has recently been demonstrated for Melbourne, Florida. In related activities the Aeronomy Laboratory has teamed with the Environmental Technology Laboratory (ETL) in hydrometeorological studies in relation to the PACJET campaign on the California coast and on microphysical process studies utilizing profilers in the North American Monsoon Experiment (NAME). Results of several field campaigns indicate that profilers used together with disdrometers and rain gauges provide an effective independent means for calibration and validation of scanning radar estimates of precipitation. We briefly outline below the combined use of these instruments to provide reference precipitation parameters with error characteristics to monitor and continuously validate precipitation estimates from the scanning radar. The domain of interest is on the order of 10 km 10 km and the data collected as outlined here could be a primary component of real-time observations supporting any ground validation Supersite. Used with disdrometers and rain gauges profilers can provide an efficient means for tying ground-based observations from in situ sensors with the scanning radar observations at altitude. The profiler has an observing volume intermediate between the relatively small observing volume of rain gauges and disdrometers and the much larger observing volume of the scanning radar. It is problematic to utilize the ground-based in situ sensors to directly calibrate the scanning radar owing to the spatial separation of the surface measurements from the radar measurements at altitude and the large difference in observing volumes between the scanning radar and the in situ surface instruments. In the remainder of this section we focus on the use of the profiler as a transfer standard by comparing disdrometer and profiler measurements and profiler measurements with scanning radar observations. Here, we utilize recent studies by Gage et al. (2004) and Williams et al. (2005) to establish the precision of profiler and disdrometer measurements and to compare profiler and scanning radar observations. We also draw upon the GPM Front Range Pilot Study observations to illustrate the utility of combining these instruments in any ground validation effort. Gage et al. (2004) compared two collocated JWDs at Wallops Island to show that the precision of reflectivity determined from these instruments is about 1.5 dbz for one minute samples. These authors also show that the reflectivity difference from the collocated instruments is nearly normally distributed in dbz space and uncorrelated in time implying that mean reflectivities measured by these instruments can be obtained with a precision equal to the standard error of the

26 NASA GPM GV Front Range Pilot Project Final Report Page 26 mean. In a similar fashion Gage et al. (2004) and Williams et al. (2005) determine the precision of one minute profiler reflectivities from collocated instruments to be approximately 0.5 dbz. The precision of the profiler measurement is better than the precision of the disdrometer reflectivities owing to the averaging inherent in a larger sampling volume. The precision values quoted above are for stratiform precipitation and values for precipitation in highly convective conditions will be larger. The precision of measurement from the disdrometer and profiler instruments is good enough that the disdrometer can be used for absolute calibration of the profiler. The disdrometer also provides validation of the DSD parameters retrieved from the profiler. Figure reproduces time-height cross-sections of the precipitation event of June 2004 as seen at the BAO by the 2835 MHz profiler. This event was unusual since it included a period of very light rain that was confined to the lowest 1-2 km above the surface followed by a sequence of convective showers that continued through the remainder of the day. Thus, a great variety of conditions are represented in this event. The profiler contribution is immediately evident in its ability to clearly show the vertical structure of the precipitating clouds and their evolution throughout the day. Time series of reflectivity from the JWD and the profiler at the second range gate centered at 316 m AGL are compared in Figure with a scatter plot shown in Figure The comparisons in Figure highlight the statistical behavior of the reflectivity differences between the two instruments as a function of the JWD reflectivity. Note that the disdrometer appears to overestimate reflectivities compared to the profiler at low reflectivities where the observations in this event are most numerous. There is also a tendency for the disdrometer to underestimate reflectivities compared to the profiler at high reflectivities. This reflectivity dependent bias is inconsistent with the statistical bias associated with different size sampling volumes. We have noted these tendencies in other data sets and it now appears that they are associated with the disdrometer. These features need to be understood better in order to have complete confidence in the disdrometer for absolute calibration of the profiler. DSD Estimates from the S-band Profiler Compared with the JWD While continuous calibration and validation of other ground-based instruments is an important component of the Supersite observational program, the most important aspect is the use of the ground-based instruments to provide a reliable estimate of DSD parameters and their uncertainties. The GPM Front Range Pilot provided an opportunity to demonstrate the ability of the ground-based instruments to yield estimates of the DSD parameters within selected case studies. Below, we examine results from the June case study. There are several different forms of the drops size distribution N(D) in common use. For the purposes of this report we cite two of the more commonly used forms. The exponential distribution (Marshall and Palmer 1948) is N(D) = N 0 exp (- D) (5.2.1)

27 NASA GPM GV Front Range Pilot Project Final Report Page 27 where N 0 is the intercept parameter and is the slope parameter. The exponential distribution has 2 parameters and a Gamma DSD (Ulbrich, 1983) has an additional shape parameter µ and is expressed µ µ 4 + µ ( ) N( D) = N D D = N D 0 exp 0 exp D (5.2.2) Dm where D m is the mean mass-weighted diameter defined by 4 N( D) D dd D m = (5.2.3) 3 N( D) D dd for any DSD and by D m = µ + 4 / (5.2.4) for a Gamma DSD. When µ= 0, the Gamma distribution reduces to the two parameter exponential distribution. In recent years it is becoming more common for the DSD to be represented by a normalized DSD as given by Bringi and Chandrasekhar (2001) as N norm (D) = N(D)/N w (5.2.5) The normalizing factor N w includes contributions from the total liquid water content W and the mean drop size of the DSD (D m ). The normalization allows the shapes of different DSDs from different rain regimes to be compared relative to the total liquid water content and the mean drop size. The normalizing factor can be estimed from W (4.0) N w = (5.2.6) 4 wdm While the normalized DSD is expressed using 3 parameters, only 2 of these are directly related to physical processes. N w represents the amplitude scaling of the DSD and D m represents the mean drop size of the DSD. To estimate the width of the DSD in the liquid water content domain m,, the parameters D m and µ need to be combined using the second moment of the DSD 2 m = N w ( D D N w D m 2 3+ µ 4 + µ ) D exp D dd Dm 4 + µ exp D dd Dm 3+ µ (5.2.7) The June case provides an excellent opportunity to compare DSD parameters using the JWD and the 2835 MHz profiler. Figure shows time series of N w, D m and m retrieved from the second range gate of the S-band profiler centered 316 meters above the ground at the BAO in comparison with the values obtained from the JWD at the surface. While more work will be needed to completely define the error characteristics of these estimates, the agreement of these estimates is encouraging. A case study is presented next to support the recommendation that two vertically-pointing profilers operating at different frequencies be collocated at the Supersites. One profiler operating

28 NASA GPM GV Front Range Pilot Project Final Report Page 28 at 449 MHz would be used to estimate the vertical air motion and the other profiler operating at 2835 MHz would be used to estimate the Doppler fall velocity of the raindrops. Algorithms would be developed to estimate the air motion from the 449 MHz profiler spectra and used with the 2835 MHz profiler spectra to estimate the DSD from near the surface to just below the freezing level. Why use Different Frequency Profilers? Vertically-pointing Doppler radars used as precipitation profilers measure the Doppler velocity of hydrometeors at discrete velocity intervals at each radar range gate. Using the principle that raindrops of a particular size fall at a pre-determined fallspeed in still air, the power detected at each discrete velocity interval can be converted into the number of raindrops at that particular raindrop size. Thus, the raindrop size distribution (DSD) estimate from profiling radar observations can be viewed as a coordinate transformation from power and fallspeed to number of drops and diameter size. Two of the largest unknowns in DSDs retrieved from profiler observations are the vertical air motion and the amount of turbulence and wind shear within the radar pulse volume. The measured Doppler velocity at each discrete velocity interval is the combination of the raindrop terminal fallspeed and the air motion acting on the raindrops. Not knowing the air motion causes errors in the retrieved DSDs because the uncorrected spectra are shifted towards larger drops in downdrafts and shifted toward smaller drops in updrafts. Since the total reflectivity is conserved in this coordinate transformation, uncorrected spectra associated with raindrops in downdrafts will have fewer large drops and less rain amount than the correctly shifted spectra. Smaller drops and more rain are erroneously estimated for uncorrected spectra that are associated with updrafts. See Williams (2002) for more details. The influence of turbulence and wind shear within the radar pulse volume tends to spread the power associated with one size of raindrop into several velocity intervals. Thus, the Doppler velocity spectrum is broader than a spectrum that would be produced in non-turbulent and motionless air. To estimate the DSD, the observed Doppler velocity spectrum is modeled as a convolution of the unknown DSD fallspeed velocity spectrum and the broadening caused by the turbulence and wind shear effects. While the influence of spectral broadening will be addressed in modeling the DSD, the influence of the air motion needs to be resolved by direct measurements of the vertical air motion using profiling radars. There are several factors that must be considered when determining which profiler operating frequency should be used to estimate the vertical air motion during precipitation events. Radars with longer wavelengths (50, 449 and 915 MHz operating frequencies) can detect the vertical air motion through scattering off of refractivity turbulence (Bragg scattering). While these radars may detect precipitation, radars with shorter wavelengths (449, 915, 2835 MHz and 35 GHz operating frequencies) are much more sensitive to hydrometeors through scattering off of the distributed hard targets (Rayleigh scattering and Mie scattering at 10 GHz and above). The use of two profilers operating at two frequencies enables the air motion to be estimated with one profiler and the precipitation motion to be estimated with the other under nearly all meteorological conditions.

29 NASA GPM GV Front Range Pilot Project Final Report Page 29 The longer wavelength profilers should be used to estimate the vertical air motion since they are more sensitive to Bragg scattering. The downside of using longer wavelength radars includes larger antennas, the increased altitude of the lowest first usable range gate and the long time required between recorded Doppler velocity spectra for the multiple radar samples that are needed to improve the signal-to-noise ratio. The benefits of using higher frequency profilers to detect the precipitation motion include stable transmitters that enable better estimates of the reflectivity and shorter dwell times needed to acquire the spectra. The use of two profilers in the FRPP to retrieve estimates of DSD parameters is illustrated next using the example of 17 June Panels (a) - (c) in Figure show the Doppler velocity spectra profiles observed from the 449, 915, and 2835 MHz profilers at 1:31 UTC on 17 June These profiles are associated with the more intense portion of the rain event of 17 June with the freezing level near 2 km as indicated by the change in Doppler velocity from ~2 to ~8 m s -1 downward corresponding with the change in fallspeed of ice/snow melting into raindrops shown in Figure 5.2.5b. The black lines placed on top of the spectra profiles correspond to the altitude adjusted fallspeeds of individual raindrops with 1, 3, and 6 mm diameters. Panel (d) of Figure shows the Doppler velocity spectra from each profiler at the altitude of 750 meters. The air motion peak is resolved by the 449 and 915 MHz profilers, but is not resolved by the 2835 MHz profiler. This shows that both the 449 and 915 MHz profilers are sensitive to the Bragg scattering with the 449 MHz profiler being more sensitive than the 915 MHz profiler. All three profilers are sensitive to the Rayleigh scattering from the raindrops and show very good agreement in the precipitation portion of the Doppler velocity spectra. While the 449 MHz profiler can detect the air motion peak in the Doppler velocity spectrum, sophisticated algorithms are still needed to identify the air motion peak that is 3 orders of magnitude smaller than the peak in the precipitation signal. Also, while this example shows a relatively clean 449 MHz spectrum, many spectra have noise fluctuations that are comparable in size with the air motion peak. A multiple peak-picking routine was applied to the 449 MHz profiler spectra with the air motion estimate shown as black asterisks in Panel (a) of Figure Only the air motion estimates up to 2.2 km altitude are shown. Note that the air motion spectral peak is very hard to distinguish in the 915 MHz profiler observations and not present in the 2835 MHz profiler observations. From this analysis, it is recommended that two profilers operating at two different frequencies be collocated at the Supersites. One profiler operating at 449 MHz would be used to estimate the vertical air motion and the other profiler operating at 2835 MHz would be used to estimate the Doppler motion of the raindrops. Algorithms would be developed to estimate the air motion from the 449 MHz profiler spectra and used with the 2835 MHz profiler spectra to estimate the DSD from near the surface to just below the freezing level.

30 NASA GPM GV Front Range Pilot Project Final Report Page 30 Figure Time-altitude cross-sections of (a) reflectivity, (b) mean Doppler velocity, and (c) spectral width measured by the BAO S-band profiler for the June 2004 rain event.

31 NASA GPM GV Front Range Pilot Project Final Report Page 31 Figure Time series of JWD reflectivity (red) and S-band profiler reflectivity (blue) at 316 m above the ground for the June 2004 rain event. Figure Scatter plot of JWD reflectivity and S-band reflectivity at 316 m above the ground for the June 2004 rain event. Panel (a) shows the probability distribution function (PDF) of the reflectivity difference (S-band - JWD), Panel (b) shows the scatter plot of values, and Panel (c) shows the PDF of JWD reflectivity.

32 NASA GPM GV Front Range Pilot Project Final Report Page 32 Figure Rain drop size parameters estimated from the JWD (blue) and S-band profiler at 316 m above the ground (red). Panel (a) shows the Generalized Intercept Parameter, Nw, Panel (b) shows mean drop diameter, Dm, and Panel (c) shows the standard deviation of the Mass spectrum,.

33 NASA GPM GV Front Range Pilot Project Final Report Page 33 Figure Simultaneous spectra obtained from the three different frequency profilers at BAO. Panels (a), (b), and (c) shows the vertical profile of spectra collected from the 449 MHz, 915 MHz, and 2835 MHz profilers. Panel (d) shows the spectra for each profiler at the altitude of 750 meters above the ground. The air motion (Bragg scattering) portion of the spectrum occurs around 0 m s-1, and the raindrop motion (Rayleigh scattering) portion of the spectrum occurs at more downward velocities Polarimetric Retrievals of DSD Parameters by the 2-D Video Disdrometer (2DVD) This section describes the measurements provided by the 2DVD as deployed for the Pilot Project, as well as describes the methodology of comparing 2DVD measurements with those from scanning polarimetric radar. Measurements from the 2DVD can be used to directly measure the drop size distribution (DSD), rain rate, and particle fall speed through its optical technology (see Figure for a schematic, Figure for fall speeds measured during the Pilot Project). Here, the DSD is assumed to be of the normalized gamma form with parameters Nw, Dm and µ. Here, Nw is the normalized intercept parameter (in mm-1 m-3), Dm is the mass-weighted mean diameter (in mm) and µ is the

34 NASA GPM GV Front Range Pilot Project Final Report Page 34 shape of the DSD. The 1-min averaged N(D) is fitted to the normalized gamma form. Histograms of D o (median volume diameter which is close to D m ) and log 10 (N w ) are shown in Figures and for all the data gathered during the GPM-Pilot program. These histograms are mostly representative of stratiform rain. The average D o is 1.08 mm and average N w is around 2500 mm - 1 m -3 (which can be compared to the Marshall-Palmer value of 8000). The average µ was found to be close to 4. For Rayleigh scattering, the functional relationships Z h =f 1 (µ)n w D m 7, K dp =f 2 (µ)n w D m 5, and Z dr -1=f 3 (µ) D m 2 can be established (note that here Z dr is the ratio of reflectivites = Z h /Z v ). From these functional relations it is possible, via simulations, to construct algorithms for estimating N w, D m and (to a lesser accuracy) the shape parameter (µ) from radar measurements of Z h, Z dr and K dp (e.g., Gorgucci et al 2002). While such algorithms have been developed and applied at S-, C- and X- bands, the effects of attenuation at C- and X-bands implies that attenuation-correction procedures must be applied for Z h and Z dr prior to estimation of the DSD parameters. The 2DVD gamma DSD can also be used to simulate the radar observables such as Z h, Z dr and K dp and these can be compared with the corresponding radar measurements made by the NOAA/X-Band radar as illustrated in Figure The standard error bars for the 2DVD data are based on expected sampling errors (sensing area= 100 cm 2, time integration=1 min) while those for the radar are based on expected fluctuation errors. The radar data have been averaged over a m area centered over the location of the 2DVD. The bias in Z h was determined to be 2.7 db. The bias in Z dr was determined to be very small (<0.1 db). Note in Figure 5.3.5b that the minimum detectable K dp at X-band is close to 0.1 deg/km or around 2 mm/h. This figure demonstrates the advantage of X-band for measuring low rain rates. At higher rain rates, the X- band signal generally falls below noise level while this is not a problem at S-band. Thus, to optimally cover the full range of rain rates with polarimetric methods, a combined S/X-band radar system with dual-polarization at both frequencies is desired. Below 2 mm/h however, the standard Z-R method can be used. Having a network of 2DVD disdrometers (rather than the usual rain gages) will enable the radar Z h and Z dr data to be very accurately calibrated. At the same time it will enable the development and validation of DSD retrieval algorithms based on dual-frequency radar measurements of Z h, Z dr and K dp. The full power of dually-polarized radar techniques for rain measurement will be significantly enhanced by deploying a network of 2DVDs within the scanning range of a dual-frequency radar system.

35 NASA GPM GV Front Range Pilot Project Final Report Page 35 Figure 5.3.1: Schematic of 2DVD measurement principle. Figure 5.3.2: Terminal fallspeed ( Y-axis in m/s) versus D (X-axis in mm) 2D video measurements. Contours of log(n) are where N is the number of measurement pairs. Line fit to Gunn-Kinzer (1949) data. Filtering has been applied to eliminate mismatched drops.

36 NASA GPM GV Front Range Pilot Project Final Report Page 36 Figure 5.3.3: Histogram of median volume diameter D o from 2DVD of the DSD fitted with a normalized gamma form. Histogram based on all data acquired with the 2DVD during the GPM-Pilot project. Histogram dominated by stratiform rain type. Figure 5.3.4: As in figure except histogram of log(n w ) where N w is the Normalized intercept parameter. The Marshall-Palmer log(n w ) = 3.9 or N w =8000 mm -1 m -3.

37 NASA GPM GV Front Range Pilot Project Final Report Page 37 Figure 5.3.5: Intercomparison of Z h (top), K dp (middle) and Z dr (bottom) from XPOL radar and 1min averaged 2DVD measurements of DSD. Standard deviation bars for 2DVD data are based on expected sampling errors, while those from radar are based on expected signal fluctuation errors. Radar data averaged over a 700 X 700 m area centered on the 2DVD location.

38 NASA GPM GV Front Range Pilot Project Final Report Page Hydrometeor Identification by Polarimetric Radar A major benefit of radar polarimetry is the ability to identify the kinds of particles present in the beam sample volume, based on polarimetric measurables that are related to particle shapes. In addition to providing a foundation for studies of storm microphysical process, proper identification of particles is essential to prevent the use of inappropriate precipitation intensity estimators, such as applying rainfall Z-R relations to regions of snowfall. Particle identification information is not available from conventional radars. Hydrometeor identification using dual-polarization radar measurements Dual-polarization radar measurements of precipitation are sensitive to the hydrometeor properties such as shape, orientation, size, phase state, and fall behavior. Therefore, dual-polarization radar measurements of precipitation can be used effectively to identify hydrometeor types in precipitation (Bringi and Chandrasekar 2001). Fuzzy logic technique is well suited for classifying hydrometeors as argued by Liu and Chandrasekar (Liu and Chandrasekar 2000), because it is easier to implement than other techniques such as statistical decision method using the prior probability as well as the probability density function. The fuzzy logic classification also has the ability to identify hydrometeor types with overlapping and noise-contaminated measurements. The hydrometeor classification system uses typically five polarimetric radar measurements as input variables, namely reflectivity at horizontal polarization (Z h ), differential reflectivity (Z dr ), specific differential phase (K dp ), linear depolarization ratio (LDR), and correlation coefficient ( hv ) and Height (or Temperature) as corresponding environment factor. A fuzzy logic classification system has three principal components namely a) fuzzification b) inference and c) defuzzification. The block diagram of CSU fuzzy logic hydrometeor classification system is shown in Figure An example of hydrometeor classification with application to the data from the CSU-CHILL radar is shown in Figure The figure also shows the comparison of the classification result with in-situ airborne observations of the T-28 aircraft using instruments such as the two dimensional cloud particle measurement probe (2DC), High volume particle sampler (HVPS). The use of X-band polarimetric estimates to infer ice hydrometeor particle habits The NOAA/ETL X-band radar operates in the simultaneous transmission and simultaneous receiving (STSR) measurement mode which allows direct measurements of such important polarimetric parameters as differential phase shift, dp (and subsequent calculations of its range derivative, K dp ), differential reflectivity, Z dr, and the co-polar correlation coefficient, hv. These parameters can be used for hydrometer identification in a manner similar to the described in section The STSR measurement scheme does not allow measurements of another popular radar polarization parameter namely linear depolarization ratio (LDR) in the H-V polarization basis. Under this measurements scheme, however, it is possible to estimate circular

39 NASA GPM GV Front Range Pilot Project Final Report Page 39 depolarization ratio (CDR) if the phase shifts between horizontal and vertical polarizations of the radar transmitter and receiver are accounted for (Matrosov 2004). Both LDR and CDR depend on hydrometeor shapes, but unlike LDR, CDR does not depend on hydrometer orientation in the radar polarization plane. This is an obvious advantage of CDR over LDR for the purpose of hydrometeor identification at closer ranges when propagation effects are not very severe. Another advantage of CDR is that it is larger than LDR which extends the applicability of depolarization measurements. ETL has extensively used CDR (and other closely associated depolarization ratios) for the purpose of identifying ice hydrometers. The methods can distinguish between prolately shaped hydrometeors, such as columns or needle ice crystals and oblately shaped hydrometeors, such as plate or dendrite crystals. It can also obtain quantitative measures of their shapes (aspect ratios). The primary scanning mode used for hydrometer identification is RHI scans that sweep from low elevation angles through the zenith. The depolarization patterns observed on these scans are used to distinguish between hydrometeor types. CDR values at an elevation of around 40 degrees can be used to estimate their aspect ratios. Most of the ETL studies have been conducted in cold non-precipitating and weakly precipitating clouds with the use of a Ka-band ( = 8.7 mm) cloud radar (Matrosov et al. 2001, Reinking et al. 2002), which has a direct measurement of depolarization ratio, but similar methods can be applied at X-band. For example, with the development of CDR capability within the STSR measurement mode (Matrosov 2004), it became possible to apply some of the Ka-band approaches for ice hydrometeor type identification and shape estimations to the NOAA/ETL X- band measurements. As an illustration, Figure shows X-band CDR estimates in the snow region of a winter storm above the melting level. The DR elevation angle dependence at 5.5 km near the echo top shows a distinct pattern of decreasing depolarization with increasing the elevation angle. Such a pattern is characteristic of single, pristine dendrites (Matrosov et al. 2002). At 3.5 km, DR values remain low and do not change much with the radar elevation angle. This pattern suggests highly aggregated and/or rimed snowflakes. This illustrates, for example, the utility of the X-band polarimetric measurements for characterizing the evolution of snowflakes as they fall. Figure Block diagram of CSU fuzzy hydrometeor classification system.

40 NASA GPM GV Front Range Pilot Project Final Report Page 40 Figure (a) Vertical structure of radar measurements (Z h, Z dr ) and the hydrometeor classification result corresponding to the case of June 22, Black line in Z h field: T28 aircraft track along time and altitude (o: start time, +: end time), Dotted line in Z dr field: Melting layer detected, (b) 2DC image data, and (c) HVPS image data PACJET-03, February 16, 2003, 00:01 UTC, azim=6 o a) red: height=5.5 km b) blue: height=3.5 km -14 CDR (s) (db) a) -22 b) radar elevation angle (deg) d/d/wallops/cdr/sldr/rhidr Figure CDR elevation angle patterns for pristine dendrites (a) and aggregated/rimed snowflakes (b).

41 NASA GPM GV Front Range Pilot Project Final Report Page 41 APPENDIX 1. Characteristics of the Scanning Radars, Profilers, and Disdrometers Used in the Pilot Project Basic operating characteristics are listed in this section for the scanning and profiling radars and the disdrometers used in the Front Range Pilot Project. These may serve as useful comparison benchmarks when designing or selecting similar instruments for dedicated use at the GPM GV sites. A number of variations of these instruments, such as other commercial disdrometers, are also available from different sources, but it is not within the scope of this report to address these. CHILL S-Band Radar General Description: NSF-funded, transportable, S-Band polarimetric research radar facility Frequency: GHz ( = 11.1 cm) Transmitters: Two upgraded FPS-18 klystron systems which separately drive the antenna s horizontal and vertical polarization ports. Typical Peak Transmit Power: 800 kw Antenna: 8.5-m diameter center-fed parabolic, 42 db net gain Beam width: 1.0 deg., circular Pulse Length: selectable from 0.33 to 1.0 µs (resolution = 50 to 150 m) PRF: Nominally ~1 khz, selectable from 0.8 to 1.2 khz Maximum Range: km Scans: PPI, RHI, sector, fixed-pointing; vertically-pointed (90 o ) elevation available Scan Rates: up to 20 deg/sec; ~6 deg/sec for typical alternating polarization PPI scans Receivers / data system: Two Lassen Research DRX digital receivers Sensitivity: approx. 0 dbz at 60 km range (receiver noise power level ~ -112 dbm) Polarizations: alternating H,V; simultaneous H and V (including adjustable phase lag between H and V for circular / elliptical net polarizations) More Detailed Information: Brunkow et al (2000); web: NOAA/ETL X-Band Radar Major Capabilities: transportable, scanning, Doppler, dual-polarization Primary Uses: measurements of precipitation, boundary layer airflow, and ocean surface features Designer/Manufacturer: NOAA/ETL Frequency: 9.34 GHz ( = 3.2 cm) Peak Transmit Power: ~20 kw Antenna: 3.1-m diameter parabolic, 44 db gain Beam width: 0.9 deg., circular Pulse Length: selectable from 0.05 to 1.0 µs (resolution = 7.5 to 150 m) PRF: selectable (2000 typical); double pulse method used to extend Nyquist Maximum Useful Range: ~50 km is typical Scans: PPI, RHI, sector, fixed-beam, all with elevations through zenith & below horizon; Scan Rates: up to 30 deg/sec; ~6 deg/sec for typical uses Beam Rate: selectable, 8 beams/sec typical Sensitivity: approx. 0 dbz at 25 km range Polarizations: simultaneous equal HV (STSR); pulse-to-pulse HV switching or circular are also available with pre-project retrofitting Doppler processing: pulse pairs (typically), and time series (limited use) Data system: VME-based with DSP & SPARC (ETL-built)

42 NASA GPM GV Front Range Pilot Project Final Report Page 42 Data rates: ~100 MB/hr for pulse pairs Operators Required: 1 person typically on site, but the radar can run unattended for a few days using programmed sequences of scans Platform: 16-m flatbed trailer, or 3 sea containers Calibration: trihedral corner reflector (Martner et al. 2003), then adjustments based on comparisons with disdrometer; Z dr values checked in drizzle More Detailed Information: Martner et al. (2001); NOAA/AL S-Band Vertically Pointing Precipitation Radar Profiler Major Capabilities: low-power, vertically-directed, continuous unattended operation, transportable, Doppler, high vertical resolution Primary Uses: continuous measurements of vertical structure of precipitating clouds, radar reflectivity, precipitation fall velocities, reference reflectivity estimates of Drop-size distributions Designer/Manufacturer: NOAA/AL (Available commercially from Vaisala) Frequency: GHz ( = 10.6 cm) Peak Transmit Power: ~380 W Antenna: 3-m shrouded dish Beam width: 3.2 deg., Pulse Length: selectable from 30 m to 500 m PRF: ~ 10,000 Hz Dwell Time: selectable, typically sec/mode Doppler processing: Full Doppler spectra (typically, 256 points) of coherently averaged return signals Data rates: moments: 0.5 MB/hr; Spectra 60 MB/hr Data Products: Doppler spectra and moments produced every sec; spectral moments are for reflectivity, Doppler velocity and spectral width Operators Required: The profiler is designed to run unattended and can be checked remotely with occasional visits to retrieve data etc. Calibration: reflectivity calibration based on collocated Joss-Waldvogel disdrometer More Detailed Information: Ecklund et al. (1999); Gage et al. (2002) NOAA/FSL UHF (449 MHz) Wind/Temperature Radar Profiler Major Capabilities: Doppler beam swinging radar profiler, five beams (zenith and Oblique beams at 75 degree elevation) continuous unattended operation, RASS capability Primary Uses: continuous real-time measurements of vertical profiles of horizontal wind and temperature, Doppler spectra recorded, part of NOAA s wind profiler network Designer/Manufacturer: NOAA/FSL (all commercially available components) Frequency: 449 MHz ( = 66.7 cm) Peak Transmit Power: ~ 5 kw Antenna: 10-m COCO array Beam width: 4.7 deg., Pulse Length: selectable 250, 500, and 1000 meters PRF: 2 modes - ~ 7,750 and 13,500 Dwell Time: 2 modes 36 and 24 sec. at each antenna position Doppler processing: Full Doppler spectra (typically, 256 points) of coherently averaged return signals Data rates: wind/temperature profiles twice/hour, moments and spectra twice/minute Total data ~ 120 Mb/day Data Products: Wind and temperature profiles every half-hour (typical), Doppler spectra

43 NASA GPM GV Front Range Pilot Project Final Report Page 43 and moments twice per minute Operators Required: The profiler is designed to run unattended Calibration: Doppler calibration with signal generator, theoretical reflectivity calibration only More Detailed Information: Weber et. al. (1990) NOAA/ETL UHF (915 MHz) Wind/Temperature Radar Profiler Major Capabilities: Doppler beam swinging radar profiler, three beams (zenith, NE and SE beams at 75 degree elevation), continuous unattended operation, RASS capability Primary Uses: continuous real-time measurements of vertical profiles of horizontal wind and temperature, Doppler spectra recorded Designer/Manufacturer: NOAA/ETL Frequency: 915 MHz ( = cm) Peak Transmit Power: ~ 5.3 kw Antenna: 10-m offset parabola with offset feeds Beam width: 2.3 deg., Pulse Length: selectable 75, 150, and300 meters PRF: ~ 16,000 Hz typically (selectable 1,000 to 30,000) Dwell Time: selectable, typically sec/mode for each beam position Doppler processing: Full Doppler spectra (typically, 256 points) of coherently averaged return signals Data rates: wind/ temperature profiles twice/hour (typical), moments and spectra twice/minute (typical), Total data ~ 400Mb/day Data Products: Wind and temperature profiles every half-hour (typical), Doppler spectra and moments every 30 sec. (typical) Operators Required: The profiler is designed to run unattended with occasional visits to retrieve data etc. Calibration: Doppler calibration with signal generator, theoretical reflectivity calibration only More Detailed Information: Strauch et. al. (1984) 2D-Video Disdrometer Common Name: 2DVD Manufacturer: Joanneum Research., Graz, Austria Capabilities: counts and determines sizes, fall speeds, shapes, and orientations of hydrometeors, including raindrops and snowflakes; parameters such as rainfall rate, radar reflectivity factor, and liquid water content are computed from the measured raindrop size distribution. Physical Basis: particles cast shadows onto optical sensor arrays as they fall through two vertically separated light beam sheets. Drop Diameter Size Range: ~ mm Diameter Resolution: ~0.2 mm Sensor Area: 100 cm 2 Integration Time: selectable from 15 sec to 24 h; 60 sec is typical Operation: unattended Field Calibration: steel balls of known diameter are dropped through sensor More information: Kruger and Krajewski (2002),

44 NASA GPM GV Front Range Pilot Project Final Report Page 44 Momentum Disdrometer Common Name: Joss-Waldvogel disdrometer Manufacturer: Distromet, Ltd., Basel, Switzerland, since 1968 Capabilities: counts and determines sizes of raindrops; parameters such as rainfall rate, radar reflectivity factor, and liquid water content are derived from the measured drop size distribution. Physical Basis: senses the impact momentum of each drop and converts to diameter by assuming drops are falling at their terminal velocities Drop Diameter Size Range: ~ mm Number of Size Bins: 127 (but routinely only 20 are output) Sensor Area: 50 cm 2 Temporal Resolution: selectable, 60 sec is typical Operation: unattended More information: Joss and Waldvogel (1967),

45 NASA GPM GV Front Range Pilot Project Final Report Page 45 APPENDIX 2: Comparisons of Algorithms for Calculating K dp and Corrections to Z eh and Z dr Measurements for the Effects of Attenuation at X-band A2.1. Comparisons of Algorithms for Calculations of K dp Specific differential phase shift on propagation, K dp, is a very important polarimetric radar parameter. K dp -based rainfall estimates generally show less variability to the details of DSDs compared to the estimates from Z-R relations (Zrnic and Ryzhkov 1996). One attractive attribute of X-band radars is that K dp values are approximately scaled with frequency so at X-band frequencies they are significantly higher than those at longer wavelengths (e.g., S-band or C- band). This allows polarimetric estimates of rain rate and accumulations to be obtained for lighter rainfall. As it was shown in section 5.1, the NOAA/ETL X-band radar provided reliable K dp estimates in the light-to-moderate rains with reflectivities in the 28 to 35 dbz range where CHILL S-band phase measurements were too noisy to calculate meaningful K dp values. Thus the use of X-band frequencies in radar polarimetry allows extending the applicability of the phasebased rainfall retrievals from moderate-to-heavy rains (where they are available with S-band measurements) to lighter rains potentially improving the accuracy of retrievals of such rains. K dp, however, is not directly measured by the polarimetric radars but rather it is calculated as a range derivative of total measured differential phase shift dp. Since dp measurements are somewhat noisy, different algorithms were suggested to filter differential phase shift data. The filtering approaches typically exclude dp points that correspond to non-meteorological echoes, spurious hardware responses and, to some extent, backscatter phase shifts. K dp are then calculated as range derivate of filtered/smoothened dp data at some range interval. CSU has a lot of experience in calculating K dp at S-band from CHILL data. CSU algorithms have been used for many years and proved to be robust. X-band polarimetry has a significantly shorter weather application history than S-band polarimetry, thus it was important to compare NOAA ETL K dp calculation algorithm (Matrosov et al. 2002) which is being used in the standard processing of X- band radar data with the CSU algorithm. For the sake of comparisons, a subset of X-band radar data from the Pilot Project was chosen as a test case. The CSU and ETL algorithms were both applied to this subset and corresponding results were compared. Figure A2.1 shows scatter-plots of K dp data calculated using these algorithms. The CSU algorithm uses an iterative procedure with telescoping sliding window range interval, which changes according to the reflectivity values (Bringi and Chandrasekar 2001). The ETL algorithm uses a fixed value of the sliding window range interval (5 km in this example). K dp in this algorithm is calculated as a slope of the best linear fit to good dp data contained in this interval using the least squares technique. It can be seen from Figure A2.1 that, in spite of differing approaches, K dp values calculated using ETL and CSU algorithms agree closely. The relative bias between ETL and CSU K dp data is about 10% and the relative standard deviation is typically within 30%. These values are comparable to ones due to statistical uncertainty of K dp data itself, thus the agreement between the two different algorithms is considered to be quite satisfactory.

46 NASA GPM GV Front Range Pilot Project Final Report Page 46 A2.2. Comparisons of Algorithms for Correcting Z eh and Z dr for Attenuation Effects Due to attenuation of radar signals in rain and a higher rate of attenuation for horizontally polarized signals compared to vertically polarized signals, measured values of Z eh and Z dr are biased low. There is a relatively strong relation between differential propagation phase and attenuation/differential attenuation, so the correction of measured values of Z eh and Z dr can be introduced based on dp data (Bringi and Chandrasekar 2001). At S-band frequencies, the effects of attenuation and differential attenuation are generally small and, in many practical situations, can be neglected. At X-band, however, these effects (especially at longer ranges) are substantial. For a long time, before the introduction of polarimetry, they had precluded a wider use of X- band radars for quantitative measurements of rain. Strict validation of the correction schemes is a complicated task. However, comparisons between correction results obtained with different algorithms can provide a certain degree of confidence in these algorithms if the results are consistent which each other. The same subset of the X-band radar data described above, was also used to compare the different attenuation/differential attenuation algorithms. The CSU method uses the reflectivity correction algorithms based on ZPHI approach (Testud et al. 2000). In contrast, the ETL method uses a relatively simple correction to the measured reflectivity, which scales it with the measured differential phase, and assumes the mean drop shape aspect ratio from (Brandes et al. 2004) applies. Figure A2.2 depicts scatter plots of horizontal polarization X-band reflectivities corrected using these two algorithms. It can be seen that the results of the corrections have a narrow distribution along the 1-to-1 line, with the majority of the differences not exceeding about 1 db, and very few of them exceeding 2 db. There is a slight (about 0.5 db) bias with CSU data being larger for Z eh >40 dbz. This bias can be explained by slightly different drop shape models used by ETL and CSU. Given that a typical uncertainty of reflectivity measurements is about 1 db, the agreement between the two methods for correcting reflectivity data should be considered quite satisfactory. Figure A2.3 shows comparisons of corrected Z dr values calculated with the CSU and ETL algorithms. Both algorithms scale differential attenuation corrections linearly with measured differential phase shift, although CSU uses an additional constraint, which assumes that Z dr at the far reaches of good data in each beam is equal to the mean Z dr that corresponds to the observed value of corrected reflectivity there. Due to this constraint, the coefficient between Z dr correction and differential phase in the CSU correction scheme changes from beam-to-beam, while in the ETL algorithm this coefficient is constant for the entire data set. Figure A2.3 shows the algorithm differences, on average, do not result in significant differences of the corrected Z dr values. The mean bias between CSU and ETL data is close to 0 db for Z dr < 2 db, and is slightly negative (of the order db, which is close to the uncertainty of differential reflectivity measurements) for larger Z dr values. Again, the agreement is quite satisfactory. Overall, the results presented in this appendix indicate consistency in the procedures. The more mature S-band algorithms showed good agreement with the ones more recently developed independently for X-band. The agreement reinforces our confidence that these basic measurement and correction schemes are working properly.

47 NASA GPM GV Front Range Pilot Project Final Report Page 47 Figure A2.1 Scatter-plot of Kdp (db/km) values calculated using the NOAA/ETL and CSU algorithms (left frame) and the corresponding difference (right frame). Colors represent density of points in log scale. Figure A2.2. Scatter-plot of corrected reflectivities (dbz) calculated using the NOAA ETL and CSU algorithms (left) and the corresponding difference (right). Colors represent density of points in log scale. Figure A2.3. Same as Figure A2 but for differential reflectivity (db).

48 NASA GPM GV Front Range Pilot Project Final Report Page 48 APPENDIX 3: Case study results A June Case Study An example of the basic pros and cons of rainfall estimations using X-Band radar: The 9 June 2004 GPM Front Range Pilot Project case. The use of polarimetric technology has been found to improve the accuracy of radar-based rainfall estimation as compared to the results obtained using only the intensity (reflectivity) of the return signal. This improvement is due the ability of dual polarization radars to remotely sense critical physical characteristics of the illuminated precipitation particles such their mean axis ratio, degree of preferential orientation, etc. Useful polarimetric measurements in precipitation can be made at both 11 cm (S-Band) and 3 cm (X-Band) radar wavelengths. Each wavelength, however, presents somewhat different strengths and weaknesses in the remote measurement of the full spectrum of rainfall rates. The exploration of these wavelength-dependent differences was a major goal of the GPM Front Range Pilot Project. At shorter radar wavelengths such as 3 cm, signal attenuation presents a major restriction to data collection in heavy rainfall conditions. Propagation through heavy rainfall (several 10s of mm per hour rain rates) removes significant amounts of energy from a short wavelength radar pulse. The path integrated attenuation losses can become so great that no detectable signal is returned to the radar. In contrast, at the longer 11 cm wavelength, signal attenuation losses are usually negligible. An example of this attenuation difference is provided in the data collected during the afternoon hours of 9 June Several strong thunderstorms moved across the Pilot Project network on this day. At 2228 UTC one such storm was located approximately 80 km south southwest of the CSU-CHILL radar. A Plan Position Indicator (PPI) plot of the CSU-CHILL reflectivity field at the 0.5 o elevation is shown in Figure A3.1.1a. The intense (> ~45 dbz) portion of the echo is organized into an asymmetrical Y shape centered near X= -22, Y= -65 km. The corresponding X-Band depiction of the storm as seen by the NOAA radar is shown in Figure A3.1.1b. Since the storm is closer to the NOAA radar, the 1.8 o PPI sweep is shown. (The rings depict the beam heights in km above the radar). The two PPI scans shown in Figure A3.1.1 were made within 20 seconds of each other. An attenuation correction based on the accumulated propagation differential phase ( dp ) has been applied to the X-Band reflectivity data. ( dp is the phase difference between the horizontally and vertically polarized return signals at a common range location. The development of a significant lag in the H phase return relative to that of the V phase indicates that the beam has propagated through many oblate (large) raindrops, implying that appreciable signal attenuation has occurred). The general Y shaped echo pattern seen in the CHILL data is also present in the X-Band plot. However, attenuation effects have weakened portions of the right hand branch of the Y. Signal attenuation is complete at NOAA radar ranges beyond the point marked by an X (X=14, Y= -28 in Figure A3.1.1b). This same geographic location is also marked with an X in the unattenuated CSU-CHILL reflectivity plot

49 NASA GPM GV Front Range Pilot Project Final Report Page 49 (X=-19, Y=-65.5 in Figure A3.1.1a). A general inspection of the PPI plots from the two radars shows that severe attenuation has resulted in a complete loss of echo detection by the X-Band radar in the southeastern portions of the storm. While the 3 cm radar signal suffers badly from attenuation effects, its dp response to raindrops is more sensitive than that of the 11 cm wavelength radar. This sensitivity is due to the fact that the dp magnitude is related to the ratio between the raindrop diameter and the radar wavelength. Thus, neglecting Mie scattering complications, the differential propagation phase shift along a common beam path will increase approximately 3.7 times as fast (11 / 3) at X-Band as it will at S-Band. The rate of differential propagation phase shift per kilometer of beam path range (specific differential propagation phase, K dp ) is a polarimetric radar quantity of particular interest. It has been established that a rain rate estimator that is relatively unaffected by typical variations in the raindrop size distribution is given by a power law of the form: R(mm hr -1 )=a(k dp ) b, where a and b are constants. Also, unlike estimators based on signal power, phase based rain rate estimators are fairly insensitive to partial beam blockage and are completely independent of radar system calibration errors. Thus, the ability to develop a useable K dp signal at lower rain rates allows X-Band radars to extend the accuracy improvements of phase-based polarimetric radar estimators into lighter rain areas. An example of this capability is shown in Figure A These data are from the same PPI sweeps shown in Figure A3.1.1; the plot area has been shifted towards the northwest to center on a region of lower reflectivities. The CSU-CHILL S-Band reflectivity pattern is depicted in Figure A3.1.2a. The overlaid jagged line indicates the point along each radial of X-Band data at which K dp first reaches or exceeds 0.1 km -1 in a run of at least five contiguous range gates. This is taken to indicate the initial range from the NOAA radar at which a useable K dp signal begins. In general, valid X-Band K dp measurements start to become available in areas where the S-Band reflectivity exceeds ~20-25 dbz. Figure A3.1.2b is a comparison of the K dp fields at the two radar wavelengths. The jagged line marking the onset of useable X-Band K dp data has been repeated from Figure A3.1.1a. It is apparent from Figures A3.1.2a and b that valid S-Band K dp information does not become available until reflectivity levels reach ~35 40 dbz closer to the echo core, while X- Band K dp data are available throughout much of the lower reflectivity areas. The accuracy of radar-based precipitation estimates can be meaningfully improved through the use of dual polarization techniques. The precipitation regimes in which these improvements can be realized vary with the radar wavelength. Shorter wavelength (X-Band) radars can detect the differential propagation phase signal at lower rain rates. While corrections for attenuation can be calculated and applied in real time, complete 3 cm wavelength signal loss can still occur when moderate or greater rainfall rates exist along sufficiently long segments of the beam path. Attenuation effects are minimal at longer wavelengths, permitting S-Band polarimetric radars to continue to make accurate precipitation rate estimates in the presence of heavy rain and hail.

50 NASA GPM GV Front Range Pilot Project Final Report Page 50 One method of exploiting the joint benefits of X and S-Band polarimetric observations for GPM ground validation purposes might be realized by designing a single, dual frequency antenna system. Such a matched-beam system would guarantee that observations made at the two frequencies were collected at identical antenna pointing angles and times. The use of a common beam path would also greatly improve the range correlation of the X-Band and S-Band K dp fields. Additionally, improvements to the X-Band sensitivity beyond the level achieved by the NOAA D radar system during the Front Range Pilot Project would be desirable to combat the precipitation-induced signal attenuation losses documented in this appendix.

51 NASA GPM GV Front Range Pilot Project Final Report Page 51 Figure A3.1.1a: CSU-CHILL 0.5 o elevation angle PPI scan at 2228:30 UTC on 9 June Plotted data are reflectivity levels (dbz) coded according to the scale on the right edge of the figure. Curved solid lines are the heights of the center of the beam in km AGL referenced to the radar elevation. Horizontal axes are labeled in km distances from the radar. The marked reference location is a common geographic point that also appears in Figure A3.1.1b. Figure A3.1.1b: NOAA-X-Band 1.8 o elevation angle PPI scan at 2228:50 UTC on 9 June Reflectivity data color scaling, beam height rings, and horizontal axes labeling as in Figure A3.1.1a. Hand-drawn solid red line marks the distant range edge of the echo pattern that is improperly represented due to the effects of signal attenuation at X-Band.

52 NASA GPM GV Front Range Pilot Project Final Report Page 52 Figure A3.1.2a: CSU-CHILL 0.5 o elevation angle PPI scan at 2228:30 UTC on 9 June Color-filled data are reflectivity levels (dbz) shaded according to the color scale along the upper edge of the plot. Irregular solid line connects the points at which the K dp values from the NOAA-D radar first reliably reach or exceed 0.1 deg km -1 in the NOAA-D PPI scan shown in Fig1b. (See text for the details). The NOAA-X-Band radar is located just northwest of this figure s plot domain. Figure A3.1.2b: CHILL 0.5 o elevation angle PPI scan at 2228:30 UTC on 9 June Color-filled data are specific differential phase (K dp ; deg km -1 ) values shaded according to the scale along the upper edge of the plot. As in Figure A3.1.1a, the irregular solid line marks the start of valid X-Band K dp data coverage.

53 NASA GPM GV Front Range Pilot Project Final Report Page 53 A June 2004 Case Study Widespread precipitation developed along the Front Range during the late morning and afternoon hours of 21 June in response to a cold frontal passage and subsequent low-level upslope flow (Figure A.3.2.1). Light-to-moderate rain observed around the BAO ground validation site between 1800 and 2000 UTC provided a good opportunity to test X-band polarimetric approaches to rainfall estimation for rain rates less than about 5 mm/h where the use of S-band polarimetry is limited due to noisy and unreliable measurements of differential reflectivity and differential phase shift. X-band differential phase data are about factor of 3 larger than those at S-band and are expected to improve rainfall retrieval results in lighter rain rate situations. Figure A shows comparisons of reflectivities calculated from the drop size distributions (DSDs) measured by the ETL JWD and CSU 2DVD with reflectivities observed by the CHILL radar (elevation 0.5 deg) and the NOAA ETL X-band radar (elevation 1.8 deg) in their corresponding pixels above the BAO site. X-band data were corrected for gaseous and rain attenuation. The S-band data were also corrected for gaseous absorption; however, no hydrometeor correction was applied due to the relative insensitivity of S-band radar to the relatively small hydrometeors observed on 21 June (Figure A.3.2.1). Rainfall rates from the two ground-based disdrometers are also shown. The JWD clock was about 2 min ahead of UTC, which explains some time offset in the corresponding JWD data. Given the differences in sampling geometry between the two radars 2 and the ground instruments, the overall agreement between different reflectivity measurements shown in Figure A is rather good. The average bias of X-band reflectivities for this time period was just about -1 db (which is accounted for in Figure A.3.2.1). Figure A shows time series of specific differential phase shift (K dp ) from CHILL and X- band phase measurements above BAO 3. While S-band K dp values from CHILL are generally very noisy, X-band K dp values exhibit significantly less noise with maximum values of around deg km -1 for rainfall maxima observed at about and 19.2 decimal UTC. Moreover, in contrast to the S-band estimates relatively few X-band K dp data points are negative, and most of these negative points correspond to Z eh < 25 dbz and rainfall rates less than about mm h -1. Figure A shows the accumulation time series as calculated from different radar estimators and ground based instruments deployed at BAO. Two different collocated high resolution (0.01 ) tipping bucket type rain gauges measured 3.88 and 3.66 mm of total accumulation for the 2-hour period between 1800 and 2000 UTC, respectively. It is reasonable to assume that a ground truth value of the total rain accumulation for this event can be 2 The center of the X-band radar beam is approximately 0.1 km above the BAO site while the center of the CHILL radar beam is about 0.3 km above the BAO ground instrumentation. Moreover, the width of the CHILL radar beam over BAO is approximately a factor of 10 larger than the corresponding X-band beam width (1 km vs 0.1 km). 3 X-band K DP is calculated using a least squares technique over a 4.95 km path length while S-band K DP is calculated using a digital filter method (Hubbert and Bringi 1995) over a 3.1 km path length. Sensitivity tests comparing the two K DP methodologies showed similar results for the 21 June event.

54 NASA GPM GV Front Range Pilot Project Final Report Page 54 approximated as a mean of these two gauge values (i.e., 3.77 mm). Total accumulation values as calculated from the DSD spectra recorded by the JWD and 2DVD were 4.04 and 3.51 mm, thus providing +7% and -7% biases with respect to the gauge mean value, correspondingly. Three (two) different rainfall rate estimators were used to calculate instantaneous rates and accumulations using X-band (S-band) radar data. One estimator for each radar used the generic NEXRAD Z-R relation (Z H =300 R 1.4 ). The other two X-band estimators utilized the mean X- band K dp -R relation (R=15K dp 0.85 ) and the polarimetric estimator which uses combined measurements of Z dr, K dp and Z H to adjust for changes in the drop shape parameter as described by Matrosov et al. (2002). For very low values of K dp (<0.06 deg km -1 ) and/or reflectivity (< 26.5 dbz) polarimetric estimates at X-band become unreliable, and both X-band polarimetric estimators use the NEXRAD Z-R relation value of rainfall rate. The second estimator used for S-band radar data is the blended CSU-CHILL algorithm (Petersen et al. 1999; Carey and Rutledge, 2000; Cifelli et al. 2002), which uses polarimetric information (Z H, Z dr, K dp ) to distinguish among different hydrometeor types and utilizes an optimization procedure to calculate rainfall rates at each pixel within the radar volume. The blended algorithm consists of 5 different rainfall estimators: 1. R(K dp ) = / 0 * 50.7*(K dp ) [0.85] (method 1) 2. R(K dp,z dr ) = 90.8*(K dp ) [0.93] *10 [-1.69*0.1* Z dr ] (method 2), where Z dr is in db units 3. R(Z H,Z dr ) = 6.7*10 [-3] *Z H [0.927] *10 [ * Z dr ] (method 3), where Z H is in linear units and Z dr is in db 4. R(Z H ) = (Z H /300) [1/1.4] (method 4), where Z H is in linear units 5. R(Z HRAIN ) = (Z HRAIN /300) [1/1.4] (method 5), where Z HRAIN is in linear units In (5) above, Z HRAIN is calculated using a difference reflectivity method (Golestani et al. 1989) to determine the ice fraction in the radar volume (Carey and Rutledge, 2000; Cifelli et al. 2002). The logic of the CSU-CHILL blended algorithm used to pick the optimal rainfall estimator for a given location is described in Figure A The Z-R relation-based estimator provided practically identical values of total accumulation (i.e., 3.40 mm) for both CHILL and X-band radar data (Figure A.3.2.4). It results in a -10% bias relative to the gauge data. The blended CSU-CHILL algorithm significantly underestimated the total accumulation due a combination of misidentification of ice hydrometeors 4 and noisy K dp values. During the UTC time period, the blended algorithm often identified large (spurious) ice fractions in the CHILL sampling volume (not shown). As indicated in Figure A.3.2.5, the combination of large ice fraction and noisy K dp requires the blended algorithm to use R(Z HRAIN ) method 5. As the ice fraction approaches 1.0, Z HRAIN approaches zero, by definition with a corresponding reduction in rainfall rate. It should be noted that the difference reflectivity and rain line methodology was originally developed for convective events with relatively large drops and corresponding large Z H (Golestani et al. 1989). The methodology may not be appropriate for 21 June-type situations of wide spread stratiform rain, with relatively weak radar echoes and small drops (described below). 4 The blended algorithm utilizes the difference reflectivity vs. Z H relation (i.e., rain-line fit) to determine the ice fraction.

55 NASA GPM GV Front Range Pilot Project Final Report Page 55 The K dp based estimator and the combined polarimetric estimator at X-band provided 3.59 mm and 4.13 mm of total accumulation resulting in -4.8% and +9% biases, respectively. It should be noted that similar small biases of the polarimetric X-band estimators were also recorded in other field experiments (e.g. Matrosov et al. 2002) while the Z-R - based algorithms often provide more significant biases than -10% which was observed in this event. Figure A shows rainfall rates calculated from ETL JWD DSDs and X-band K dp -based estimates of R at the times when the above mentioned conditions for the use of polarimetric information were satisfied. It can be seen that X-band differential phase information provides useful information when rainfall rates are greater then about 2 mm h -1, and K dp -based values of R are in good agreement with the disdrometer data. Overall, the 21 June event illustrates the value of X-band radar polarimetric measurements in extending rainfall accumulation estimates to lighter rains (2 mm h -1 <R< 5 mm h -1 ). ), where similar S-band polarimetry methods fail to due to differential phase measurements that are too weak and noisy. In this case, the polarimetric estimates from either radar were not significantly better than those based entirely on reflectivity. Generally, however, reflectivity-based estimates have significantly poorer accuracies due to many degrading factors that are less problematic for polarimetric radar methods (Section 3). Although the S-band K dp -based rainfall estimates were not useable in this case, S-band did nevertheless provide good polarimetric-based (Z dr ) information on characteristic drop sizes in this light rainfall. Figure A shows comparisons of mass-weighted drop diameters, D m, calculated from the ETL JWD DSD spectra with D m values retrieved from X- and S-band radar measurements of differential reflectivity, Z dr. While there is a general agreement between the JWD and radar estimates (within mm) and similar trends in drop diameters are seen, a few of X-band estimates might be biased low. Additional analyses are needed to understand the cause of the X-band Z dr low bias.

56 NASA GPM GV Front Range Pilot Project Final Report Page 56 Figure A CHILL radar reflectivity PPI at 18:16:52 UTC on 21 June Color shading represents radar reflectivity (dbz) as indicated by the scale in the figure. The CHILL radar is located at the zero range ring. Approximate locations of the BAO, PLT, and X-band (X-POL) facilities, as well as BRF and TCA rain gauge locations are indicated. 21-June-2004, BAO comparisons (X-band, az=157 o, el=1.8 o, RC h =79.7 db, RC v =79.4 db) Reflectivity, Z eh (dbz) X-band reflectivity corrected for attenuation (BAO pixel at 6.45 km) + 1dB (to remove bias) Z e (Rayleigh spheres) from ETL JWD spectra corrected for dead time and altitude CHILL S-band measured reflectivities (BAO pixel, ael=0.5 o, range=54.5 km) Z e calculated from 2dVD disdrometer at BAO rainfall rate from ETL JWD rainfall rate from 2DVD disdrometer Rainfall rate (mm h -1 ) Time (UTC) Figure A Time series of reflectivitivites and rain rates observed over BAO during the event of 21 June 2004.

57 NASA GPM GV Front Range Pilot Project Final Report Page June-2004, BAO comparisons from X-band (LMS standard processing =4.95 km) from CHILL (digital filter method =3.1 km) K DP (deg/km) time (UTC) (X-band, az=157 o, el=1.8 o ) Figure A Time series of S-band and X-band K dp values over BAO during the event of 21 June June-2004, BAO comparisons from K DP -based estimator (R=c 15 K 0.85 DP ) (X-band data) from combined polarimetric estimator R=c 1.1K 0.5 DP Z eh Z dr from NEXRAD Z-R relation (Z eh =300R 1.4 ) (X-band data) from ETL JWD (corrected for altutude and dead-time) from ET6 black gauge #3 (0.01 resolution) from ET7 standard gauge (0.01 resolution) 4 from NEXRAD Z-R relation (Z eh =300R 1.4 ) (CHILL data) 4 from CHILL blended algorithm from 2DVD data accumulation (mm) (X-band, az=157 o, el=1.8 o ) (X-band polarimetric estimators are used when K DP >0.06 deg/km and Z eh >25 dbz) time (UTC) Figure A Rainfall accumulations from different radar estimators and ground-based instruments, 21 June 2004.

58 NASA GPM GV Front Range Pilot Project Final Report Page 58 Figure A Flow diagram of the CSU-CHILL blended algorithm. Rain estimation method is indicated by number corresponding to equations in the text. 21-June-2004, BAO comparisons 8 7 from K DP -based estimator (R=c 15 K DP 0.85 ) (X-band data) from ETL JWD (corrected for altutude and dead-time) 6 rainfall rate (m h -1 ) (X-band, az=157 o, el=1.8 o ) (X-band polarimetric estimators are used when K DP >0.06 deg/km and Z eh >25 dbz) time (UTC) Figure A Rainfall rates from JWD and X-band K dp -based estimators, 21 June 2004.

59 NASA GPM GV Front Range Pilot Project Final Report Page June-2004, BAO comparisons 3 D m (from X-band Z DR measurements) D m (from ETL JWD (corrected for altutude and dead-time) D m (from CHILL Z DR measurements) mass-weighted drop diamter (mm) 2 1 (X-band, az=157 o, el=1.8 o ) time (UTC) Figure A Mass-weighted drop diameters as estimated from JWD DSDs and radar Z dr data, 21 June A July Case Study The 15 July case over the Platteville site offered an opportunity to compare both radars at similar ranges (30.4 km from CHILL, 27.8 km from X-Band) over the NOAA-AL profilers, 2 momentum JW disdrometers, and the 2D Video Disdrometer. Vertical Structure and DSD characteristics The July 2004 case brought two differing storm cells to the Platteville Atmospheric Observatory in terms of their vertical structure and rainfall intensity. Figure A3.3.1, panel (a) shows a time-height plot of S-Band reflectivity measured by the NOAA-AL profiler for the time period encompassing the two cells of interest. A newly developed convective cell passed directly over the ground site, beginning around 2230 UTC 15 July (Julian Day ). Note reflectivities above 40 dbz above the melting level (4.5 km MSL in this case), indicating the presence of convective precipitation in the first portion of this cell. Also note the tilt of the reflectivity profiles in the lowest 1 km AGL likely due to the advection of precipitation particles by the ambient wind. The second precipitation event was longer lasting and had stratiform rain characteristics as indicated by the presence of a radar bright band throughout most of the event. Figure A3.3.1, panel (b) shows a time-height plot of S-Band/449 MHz retrieval of massweighted mean diameter D m, retrieved using a technique similar to Schafer et al. (2002) described in section 3. Note that the time axes are identical in the two plots; however, the height scale has been expanded in panel (b) to examine the rain region in detail. In the first convective cell, there exist high mass-weighted particle sizes (many regions exceeding 2.5 mm) within the high-reflectivity convective region. Within the more stratiform event later on, modal D m values range within smaller values, between 1 to 2 mm.

60 NASA GPM GV Front Range Pilot Project Final Report Page 60 To compare both surface based and scanning radar estimates of D m, Figure A3.3.2 shows time series of (a) profiler and polarimetric radar estimated- and (b) profiler and altitude-corrected surface disdrometer estimated-d m for the time period shown in Figure A For the radarprofiler comparison in panel (a), profiler D m estimates were averaged over the vertical extent of the half-power beamwidth of the main lobe of the CHILL 1.5 and X-Band 1.8 beams, respectively. The height ranges of the scanning radar beams used to vertically average the profiler data are indicated by the white boxes in panel (b) of Figure A Mass-weighted mean diameter is calculated from the radar measured differential reflectivity Z dr using the following equations: D m =1.619Z dr at S-Band (Bringi and Chandrasekar 2001), and D m =1.6Z dr 0.49 at X-Band (Matrosov 2004, personal communication). From the time series in Figure A3.3.2a, it is apparent that the profiler sees overall higher peak D m values within the convective cell than either radar. The X-Band estimates higher D m than CHILL for this cell, and the X-Band D m values are overall in better agreement with the profilers during this first cell than CHILL. In the later stratiform precipitation, again the X-Band D m values are higher than CHILL s, however CHILL s values seem in better agreement with the profiler measurements made at its scan elevation than for X-Band. Figure A3.3.2b shows a D m profiler-surface disdrometer time series comparison for these two cells. Profiler D m retrievals are taken from the lowest bin in the data, at 186 m AGL, while the AL- and CU-Joss-Waldvogel Disdrometers (AL-JWD and CU-JWD) and the CSU 2DVD are collocated at the surface. Good agreement overall is present between the profilers given the sample volume differences between the profiler and the disdrometer, with the profiler retrievals having somewhat more temporal variability. Figure A3.3.3 shows the statistical variability of the comparisons and time period shown in Figure A3.3.2 for the (a) profiler-polarimetric radar and (b) profiler-surface disdrometer D m retrievals. The figure uses a box and whiskers presentation, where the central line and vertical extent and center line of the colored box represent the median and quartile variability in the parameter, while the whiskers denote the range 1.5 times the difference between the upper quartile and the lower quartile. Individual outliers that fall outside the whiskers are also plotted. In Figure A3.3.3a, it is apparent that for this case, CHILL s retrieved D m values are closer to those estimated by the profiler retrieval within its sample volume than those measured by X- Band. The distribution of X-Band retrievals is offset by about 0.5 mm above the profiler retrievals. In Figure A3.3.2b, profiler-disdrometer comparisons show median D m values are consistent to within 0.2 mm of each other. However, there is significantly more measurement variability associated with the profiler measurements relative to the JWDs, and especially the 2DVD disdrometer. However, this comparison is one demonstration of making quantitative comparisons of DSD estimates and errors as one moves from ground in situ instrumentation, to the vertically pointing profiler and scanning polarimetric radar remote sensing estimates. Developing methodology to propagate measurement, parameterization, and sampling errors from one instrument to another will be a subject of future work in this vein.

61 NASA GPM GV Front Range Pilot Project Final Report Page 61 Polarimetric variables and Rain Accumulations Figure A3.3.4a shows time series of measured equivalent reflectivities from CHILL and X-Band radars from the 1.5 and 1.8 elevation scans, respectively, as well as calculated equivalent altitude-corrected reflectivities from the surface disdrometers. There is good agreement for most time periods within a few db, and there is no significant offset in reflectivity for this case between all these instruments and the S-Band profiler reflectivities shown in Figure A3.3.1 (not shown). However, during the case the X-Band suffered from second trip contamination from a cell just beyond the 28 km range of the X-Band scan in the direction of Platteville. Second-trip reflectivity contamination is obvious in the circled regions of Figure A3.3.4a. In Figure A3.3.4b, a time series of CHILL and X-Band measured Z dr and least-squares fit (4.95 km window) K dp is plotted. Second trip-contaminated points are circled in Figure A3.3.4b. Excluding these regions, it is apparent that X-Band s Z dr values are, on average, about 0.3 db higher than CHILL Z dr s. A difference of about +0.3 db in X-Band Z dr would explain the relatively high D m values reported by X-Band in this case in Figures A3.3.2 and A K dp values within the first cell in Figure A3.3.4b differ by a factor of ~2 between X-Band and CHILL, while for the second cell CHILL K dp is useless because there is insufficient liquid water path in this low reflectivity region, while a few points of useable X-Band K dp appear late in this event. Figure A3.3.5 shows estimated rain accumulations from the various instruments for the 15 July case. For the first cell, CHILL and X-Band Z-based rain estimates (using Z=300 R 1.4 ) overestimate surface rain by factors of roughly 2.5 and 4, respectively. The major reasons for this disagreement are inappropriateness of the Z-R relationship for this cell and the tilting and vertical gradient of rainfall rate within the cell producing significant echo above, but not at the surface (Figure A3.3.1a). For this cell, the CHILL polarimetric blended algorithm (described in Appendix 3.2) only reduces the Z-R estimate by 0.3 mm. An X-Band blended algorithm was applied in this case as well. The algorithm logic is identical to the CHILL blended algorithm, except (1) polarimetric Z e thresholds of 26.5 dbz and K dp threshold of 0.06 km -1 were used instead of 38 dbz and 0.3 km -1 in the CHILL algorithm, and (2) the X-Band polarimetric rain rate equations are the following: R(K dp ) =15K dp 0.88 R(K dp,z dr ) = 28.6K dp ( 1.37)Z dr. R(Z h,z dr ) = ( )Z 1.07 h ( 5.97)Z dr In this first cell, the X-Band Blended Algorithm rain accumulations fall within the variability of the surface rain accumulations. The rain estimation method for this case is shown in Figure A3.3.5 at both S- and X-Band. For this first cell, both radars use R-K dp within the heaviest rain areas. Given the observed factor-of-2 difference in K dp between the two wavelengths (a factor of 3 is predicted by theory), the high R-K dp rain rates estimated from CHILL relative to the other instrumentation, it is hypothesized that there is some error in K dp measured at S-Band in the core of this cell. This may be because of the small size of the cell introducing error in the calculation of K dp, or water-coated ice particles might be leading to anomalously high K dp within the core of this cell within the S-Band sample volume (e.g. Vivekanandan et al. 1990). In the second cell, both Z-based and blended algorithms are biased high because of second-trip contamination in reflectivity (Z) and specific differential phase measurements (K dp ). Arrows in

62 NASA GPM GV Front Range Pilot Project Final Report Page 62 the figure highlight these areas. When these points are removed by discarding points where the X-Band reflectivities are more than 10 dbz higher than the CHILL reflectivities in the Z-R rain accumulations, and also removing points with anomalously high K dp values given Z (points with K dp > Z were found to be anomalous), then much more reasonable agreement is produced between the X-Band and other rain estimates, especially for the blended polarimetric estimator (values with second trip. This example points out that the ability to perform self-consistency checks between the two frequencies and among polarimetric estimators in removing artifacts (such as second trip, clutter, non-meteorological targets) will be a great advantage in GPM versus using traditional single-frequency Doppler radars. Within this second cell, blended radar estimates fall close to or within the variability of the surface gauge and disdrometer measurements. The blended algorithms appropriately decrease high Z-R estimates from both radars through inclusion of primarily Z-Z dr methodology (see Figure A3.3.6). K dp is also used within the X-Band measurements, however reflectivities in this cell are too weak to utilize K dp at S-Band (see Figure A3.3.4a). Figure A (a) Time series of the vertical profile of measured equivalent reflectivity [dbz e ] from the NOAA-AL S-Band profiler. (b) Time-series of the vertical profile of S-Band/449 MHz retrieved mass-weighted mean diameter D m. The sample volume of the scanning radars is denoted by the white rectangles. Figure A (a) Time series of D m estimates from the polarimetric scanning radars and average profiler-retrieved values within the scanning radars sample volume. (b) Time series of Dm estimates from the surface disdrometers and profiler-retrieved values at 186 m AGL.

63 NASA GPM GV Front Range Pilot Project Final Report Page 63 Figure A Box-and-whisker plot of retrieved and measured D m values for the time period shown in Figure A See expanded description in the text. Figure A Time series of measured (a) equivalent radar reflectivity factor (dbz e ) for the polarimetric scanning radars and disdrometers and (b) differential reflectivity (Z dr ) and specific differential phase (K dp ) from the scanning polarimetric radars. Areas of second-trip contamination in the X-Band measurements are circled.

64 NASA GPM GV Front Range Pilot Project Final Report Page 64 Figure A Time series of rain accumulation (mm) from the S- and X-Band radars Z-R and blended polarimetric techniques, surface rain gauges, and surface disdrometers. Second-trip-corrected X-Band rain rates appear for the second cell. See the text for more details. Figure A Rain-estimation method used in the CHILL and X-Band blended algorithm.