Search

Compared to traditional single-polarization radar, dual-polarization radar has a number of advantages for quantitative precipitation estimation because more information about the drop size distribution and hydrometeor type can be gleaned. In this paper, an improved dual-polarization rainfall methodology is proposed, which is driven by a region-based hydrometeor classification mechanism. The objective of this study is to incorporate the spatial coherence and self-aggregation of dual-polarization observables in hydrometeor classification and to produce robust rainfall estimates for operational applications. The S-band dual-polarization data collected from the NASA Polarimetric (NPOL) radar during the GPM Iowa Flood Studies (IFloodS) ground validation field campaign are used to demonstrate and evaluate the proposed rainfall algorithm. Results show that the improved rainfall method provides better performance than a few single- and dual-polarization algorithms in previous studies. This paper also investigates the impact of radar beam broadening on various rainfall algorithms. It is found that the radar-based rainfall products are less correlated with ground disdrometer measurements as the distance from the radar increases.

The CSU–CHILL radar is a dual-wavelength, dual-polarization weather radar system operating at S and X bands with coaxial beams. One of the capabilities of this radar system is the possibility of developing and/or validating algorithms across dual wavelengths and dual polarizations. This paper presents one such instance, showing how the rainfall field can be estimated either from the S- and X-band reflectivities or from the differential propagation phase at X band. To do so, the paper first presents a dual-wavelength attenuation correction method that uses the reflectivity measured at S band, as the constraint for the correction of the reflectivity measured at X band, and it describes how Mie scattering regions at X band may be detected from the retrieved path-integrated attenuation field. Then, the paper describes how the resulting specific attenuation field relates to rainfall and specific phase at X band, which can be obtained from dual-polarization data at a single wavelength as well, and shows examples. Finally, the paper looks at the relation between attenuation and the differential phase as a function of elevation angle for a few cases, which may be related to the drop size distribution and mean diameter, as well as temperature.

One pivotal factor affecting the accuracy of HF radar current measurements is the direction of arrival (DOA) estimation performance of the current signal. The beamforming technology or superresolution algorithm cannot always perform best in practical applications because of the phase errors existing in array channels. These phase errors, which cause uncertain estimation of DOA, lead to confused values in radial current maps. To solve this problem, this paper is focused on discussing the performances of two autocalibration methods using sea echoes for multifrequency high-frequency (MHF) radar current measurements. These two array calibration methods, based on maximum likelihood (ML) and multiple signal classification (MU), first seek single-DOA sea echoes and then gather them for array calibration using different cost functions. The ML and MU methods provide approximate mean phases, while the standard phase errors of the MU method are smaller. After array calibration using these two methods, the results show significant improvements in current retrievals. Comparisons between the MHF radar and ADCPs reveal that array calibration using the ML and MU methods also improves the estimation of radial currents clearly, with correlation coefficients over 0.93 and rms differences of 0.09–0.18 m s−1 at different operating frequencies and sampling locations. The performance of the bearing offset is also improved. Only small bearing offsets less than 10° exist in radial current measurements. Therefore, this paper demonstrates that array calibration is a crucial part for current measurements, especially for direction-finding HF radar.

In a two-part paper, radar rain-rate retrievals using specific attenuation A suggested by Ryzhkov et al. are thoroughly investigated. Continuous time series of overlapping measurements from two twin polarimetric X-band weather radars in Germany during the summers of 2011–13 are used to analyze various aspects of rain-rate retrieval, including miscalibration correction, mitigation of ground clutter contamination and partial beam blockage (PBB), sensitivity to precipitation characteristics, and the temperature assumptions of the R(A) technique. In this paper, the relations inherent to the R(A) method are used to estimate radar reflectivity Z from A and compare it to the measured Z in order to estimate PBB and calibration offsets for both radars. The fields of Z estimated from A for both radars are consistent, and the differences between Z(A) and measured Z are in good agreement with the ones calculated using either consistency relations between reflectivity at horizontal polarization ZH, differential reflectivity ZDR, and specific differential phase KDP in rain or a digital elevation model in the presence of PBB. In the analysis, the dependence of A on temperature appears to have minimal effects on the overall performance of the method. As expected, the difference between Z(A) and attenuation-corrected measured Z observations varies with rain type and exhibits a weak systematic dependency on rainfall intensity; thus, averaging over several rain events is required to obtain reliable estimates of the Z biases caused by radar miscalibration and PBB.

In a series of two papers, rain-rate retrievals based on specific attenuation A at radar X-band wavelength using the R(A) method presented by Ryzhkov et al. are thoroughly investigated. Continuous time series of overlapping measurements from two polarimetric X-band weather radars in Germany during the summers of 2011–13 are used to analyze various aspects of the method, like miscalibration correction, ground clutter contamination, partial beam blockage (PBB), sensitivity to precipitation characteristics, and sensitivity to temperature assumptions in the retrievals. In Part I of the series, the relations inherent to the R(A) method were used to calculate radar reflectivity Z from specific attenuation and it was compared with measured reflectivity to estimate PBB and calibration errors for both radars. In this paper, R(A) rain estimates are compared to R(Z) and R(KDP) retrievals using specific phase shift KDP. PBB and calibration corrections derived in Part I made the R(Z) rainfall estimates almost perfectly consistent. Accumulated over five summer months, rainfall maps showed strong effects of clutter contamination if R(KDP) is used and weaker impact on R(A). These effects could be reduced by processing the phase shift measurements with more resilience toward ground clutter contamination and by substituting problematic R(KDP) or R(A) estimates with R(Z). Hourly and daily accumulations from rain estimators are compared with rain gauge measurements; the results show that R(A) complemented by R(Z) in segments with low total differential phase shift correlates best with gauges and has the lowest bias and RMSE, followed by R(KDP) substituted with R(Z) at rain rates below 8 mm h−1.

Simulation of radar beam propagation is an important component of numerous radar applications in meteorology, including height assignment, quality control, and especially the so-called radar forward operator. Although beam propagation in the atmosphere depends on the refractive index and its vertical variation, which themselves depend on the actual state of the atmosphere, the most common method is to apply the 4/3 earth radius model, based on climatological standard conditions. Serious deviations from the climatological value can occur under so-called ducting conditions, where radar beams at low elevations can be trapped or propagate in a waveguide-like fashion, such that this model is unsuitable in this case. To account for the actual atmospheric conditions, sophisticated methods have been developed in literature. However, concerning the practical implementation of these methods, it was determined that the description in the literature is not always complete with respect to possible pitfalls for practical implementations.In this paper, a revised version of an existing method (one example for the above-mentioned “pitfall” statement) is introduced that exploits Snell’s law for spherically stratified media. From Snell’s law, the correct sign of the local elevation is a priori ambiguous, and the revised method explicitly applies (i) a total reflection criterion and (ii) another ad hoc criterion to solve the problem.Additionally, a new method, based on an ordinary differential equation with respect to range, is proposed in this paper that has no ambiguity.Sensitivity experiments are conducted to investigate the properties of these three methods. The results show that both the revised and new methods are robust under nonstandard conditions. But considering the need to catch an elevation sign ambiguity in the revised method (which cannot be excluded to fail in rare instances), the new method is regarded as more robust and unproblematic, for example, for applications in radar forward operators.

Several studies in the literature have developed approaches to diagnose hail storms from satellite-borne passive microwave imagery and build nearly global climatologies of hail. This paper uses spaceborne Ku-band radar measurements from the Global Precipitation Measurement (GPM) mission Dual-Frequency Precipitation Radar (DPR) to validate several passive microwave approaches. We assess the retrievals on the basis of how tightly they constrain the radar reflectivity at −20°C and how this measured radar reflectivity aloft varies geographically. The algorithm that combines minimum 19-GHz polarization corrected temperature (PCT) with a 37-GHz PCT depression normalized by tropopause height constrains the radar reflectivity most tightly and gives the least appearance of regional biases. A retrieval that is based on a 19-GHz PCT threshold of 261K also produces tightly clustered profiles of radar reflectivity, with little regional bias. An approach using regionally adjusted minimum 37-GHz PCT performs relatively well, but our results indicate it may overestimate hail in some subtropical and midlatitude regions. A threshold applied to the minimum 37-GHz PCT (≤230 K), without any scaling by region or probability of hail, overestimates hail in the tropics and underestimates beyond the tropics. For all retrieval approaches, storms identified as having hail tended to have radar reflectivity profiles that are consistent with general expectations for hailstorms (reflectivity > 50 dBZ below the 0°C level, and > 40 dBZ extending far above 0°C). Profiles from oceanic regions tended to have more rapidly decreasing reflectivity with height than profiles from other regions. Subtropical, high-latitude, and high-terrain land profiles had the slowest decreases of reflectivity with height.

Over the last two decades, the Central Weather Bureau of Taiwan and the U.S. National Severe Storms Laboratory have been involved in a research and development collaboration to improve the monitoring and prediction of river flooding, flash floods, debris flows, and severe storms for Taiwan. The collaboration resulted in the Quantitative Precipitation Estimation and Segregation Using Multiple Sensors (QPESUMS) system. The QPESUMS system integrates observations from multiple mixed-band weather radars, rain gauges, and numerical weather prediction model fields to produce high-resolution (1 km) and rapid-update (10 min) rainfall and severe storm monitoring and prediction products. The rainfall products are widely used by government agencies and emergency managers in Taiwan for flood and mudslide warnings as well as for water resource management. The 3D reflectivity mosaic and QPE products are also used in high-resolution radar data assimilation and for the verification of numerical weather prediction model forecasts. The system facilitated collaborations with academic communities for research and development of radar applications, including quantitative precipitation estimation and nowcasting. This paper provides an overview of the operational QPE capabilities in the Taiwan QPESUMS system.

Partial beam blockage (PBB) correction is an indispensable step in weather radar data quality control and subsequent quantitative applications, such as precipitation estimation, especially in urban and/or complex terrain environments. This paper developed a novel PBB correction procedure based on the improved ZPHI method for attenuation correction and regional specific differential propagation phase (KDP)–reflectivity (ZH) relationship derived from in situ raindrop size distribution (DSD) measurements. The practical performance of this PBB correction technique was evaluated through comparing the spatial continuity of reflectivity measurements, the consistency between radar-measured and DSD-derived KDP and ZH relationships, as well as rainfall estimates based on R(ZH) and R(KDP). The results showed that through incorporating attenuation and PBB corrections (i) the spatial continuity of ZH measurements can effectively be enhanced; (ii) the distribution of radar-measured KDP versus ZH is more consistent with the DSD-derived KDP versus ZH; (iii) the measured ZH from a C-band radar in the PBB-affected area becomes more consistent with collocated S-band measurements, particularly in the rainstorm center area where ZH is larger than 30 dBZ; and (iv) rainfall estimates based on R(ZH) in the PBB-affected area are incrementally improved with better spatial continuity and the performance tends to be more comparable with R(KDP).

This paper examines nowcasts of precipitation from the High-Resolution Rapid Refresh (HRRRv2) model from the summer of 2017 along the Colorado Front Range. It was found that model nowcasts (2 h or less) of precipitation amount were less skillful than extrapolation of the KFTG WSR-88-D data at a spatial scale of 120 km. It was also found that local-scale (mesoscale) influences on rainfall intensity and amount have a much greater impact on rainfall intensity than large-scale (synoptic) influences. Thus, large-scale trends are not useful for modifying extrapolation nowcasts on the local scale. Errors in the HRRR nowcasts are attributed to an inability of the model and data assimilation to resolve convergence along outflow boundaries and other terrain-influenced mesogamma-scale flows that contribute to storm formation and evolution. While the HRRRv2 1-h nowcasts were strongly correlated with observed precipitation events, the nowcast precipitation amounts were in error by more than a factor of 2 about 50% of the time, with half of the cases being overestimates and half being underestimates. A large fraction of the HRRRv2 overestimates were associated with stratiform rain events. It is speculated that this was a result of misinterpretation of the radar bright band as more intense precipitation aloft by the data assimilation scheme. A large fraction of the HRRRv2 underestimates occurred when the data assimilation and model were unable to fully resolve the low-level convergence along small-scale, narrow boundaries that led to new storm initiation and/or storm growth.

This paper describes a technique for estimating the liquid water content (LWC) and a characteristic particle diameter in stratocumulus clouds using radar and lidar observations. The uncertainty in LWC estimate from radar and lidar measurements is significantly reduced once the characteristic particle diameter is known. The technique is independent of the drop size distribution. It is applicable for a broad range of W-band reflectivityZbetween −30 and 0 dBZand all values of lidar backscatterβobservations. No partitioning of cloud or drizzle is required on the basis of an arbitrary threshold ofZas in prior studies. A method for estimating droplet diameter and LWC was derived from the electromagnetic simulations of radar and lidar observations. In situ stratocumulus cloud and drizzle probe spectra were input to the electromagnetic simulation. The retrieved droplet diameter and LWC were validated using in situ measurements from the southeastern Pacific Ocean. The retrieval method was applied to radar and lidar measurements from the northeastern Pacific. Uncertainty in the retrieved droplet diameter and LWC that are due to the measurement errors in radar and lidar backscatter measurements are 7% and 14%, respectively. The retrieved LWC was validated using the concurrent G-band radiometer estimates of the liquid water path.

The Multi-Radar-Multi-Sensor (MRMS) system was transitioned into operations at the National Centers for Environmental Prediction in the fall of 2014. It provides high-quality and high-resolution severe weather and precipitation products for meteorology, hydrology, and aviation applications. Among processing modules, the radar data quality control (QC) plays a critical role in effectively identifying and removing various nonhydrometeor radar echoes for accurate quantitative precipitation estimation (QPE). Since its initial implementation in 2014, the radar QC has undergone continuous refinements and enhancements to ensure its robust performance across seasons and all regions in the continental United States and southern Canada. These updates include 1) improved melting-layer delineation, 2) clearance of wind farm contamination, 3) mitigation of corrupt data impacts due to hardware issues, 4) mitigation of sun spikes, and 5) mitigation of residual ground/lake/sea clutter due to sidelobe effects and anomalous propagation. This paper provides an overview of the MRMS radar data QC enhancements since 2014.

In a 2018 paper by Bukovčić et al., polarimetric bivariate power-law relations for estimating snowfall rate S and ice water content (IWC), and , were developed utilizing 2D video disdrometer snow measurements in Oklahoma. Herein, these disdrometer-based relations are generalized for the range of particle aspect ratios from 0.5 to 0.8 and the width of the canting angle distribution from 0° to 40° and are validated via analytical/theoretical derivations and simulations. In addition, a novel S(KDP, Zdr) polarimetric relation utilizing the ratio between specific differential phase KDP and differential reflectivity Zdr, , is derived. Both KDP and are proportionally affected by the ice particles’ aspect ratio and width of the canting angle distribution; therefore, the variables’ ratio tends to be almost invariant to the changes in these parameters. The S(KDP, Z) and S(KDP, Zdr) relations are applied to the polarimetric S-band WSR-88D data obtained from three geographical locations in Virginia, Oklahoma, and Colorado, and their performance is compared with estimations from the standard S(Z) relations and ground snow measurements. The polarimetric estimates of snow accumulations from the three cases exhibit smaller bias in comparison with the S(Z), indicating good potential for more reliable radar snow measurements.

During spring 2016 and spring 2017, a vertically pointing, S-band Frequency Modulated Continuous Wave radar (UMass FMCW) was deployed in northern Alabama under the auspices of the Verification of the Origins of Rotation in Tornadoes Experiment (VORTEX)-Southeast. In total, ~14 weeks of data were collected, in conditions ranging from quiescent clear skies to severe thunderstorms. The principal objective of these deployments was to characterize the boundary layer evolution near the VORTEX-Southeast domain. In this paper, we describe intermediate results in service of this objective. Specifically, we describe updates to the UMass FMCW system, document its deployments for VORTEX-Southeast, and apply four automated algorithms: 1) a dealiasing algorithm to the Doppler velocities, 2) a fuzzy logic scatterer classification scheme to separate precipitation from nonprecipitation observations, 3) a brightband/melting-layer identification algorithm for stratiform precipitation, and 4) an extended Kalman filter–based convective boundary layer depth (mixing height) measurement algorithm for nonprecipitation observations. Results from the latter two applications are qualitatively verified against retrieved soundings from a collocated thermodynamic profiling system.

Dealiasing is a common procedure in radar radial velocity quality control. Generally, there are two dealiasing steps: a continuity check and a reference check. In this paper, a modified version that uses azimuthal variance of radial velocity is introduced based on the integrating velocity–azimuth process (IVAP) method, referred to as the V-IVAP method. The new method can retrieve the averaged winds within a local area instead of averaged wind within a full range circle by the velocity–azimuth display (VAD) or the modified VAD method. The V-IVAP method is insensitive to the alias of the velocity, and provides a better way to produce reference velocities for a reference check. Instead of a continuity check, we use the IVAP method for a fine reference check because of its high-frequency filtering function. Then a dealiasing procedure with two steps of reference check is developed. The performance of the automatic dealiasing procedure is demonstrated by retrieving the wind field of a tornado. Using the dealiased radar velocities, the retrieved winds reveal a clear mesoscale vortex. A test based on radar network observations also has shown that the two-step dealiasing procedure based on V-IVAP and IVAP methods is reliable.

Kelvin-Helmholtz (KH) waves are a frequent source of turbulence in stratiform precipitation systems over mountainous terrain. KH waves introduce large eddies into otherwise laminar flow, with updrafts and downdrafts generating small-scale turbulence. When they occur in-cloud, such dynamics influence microphysical processes that impact precipitation growth and fallout. Part I of this paper used dual-Doppler, 2D wind and reflectivity measurements from an airborne cloud radar to demonstrate the occurrence of KH waves in stratiform orographic precipitation systems and identified four mechanisms for triggering KH waves. In Part II, we use similar observations to explore the effects of KH wave updrafts and turbulence on cloud microphysics. Measurements within KH wave updrafts reveal the production of liquid water in otherwise ice-dominated clouds, which can contribute to snow generation or enhancement via depositional and accretional growth. Fallstreaks beneath KH waves contain higher ice water content, composed of larger and more numerous ice particles, suggesting that KH waves and associated turbulence may also increase ice nucleation.A Large-Eddy Simulation (LES), designed to model the microphysical response to the KH wave eddies in mixed phase cloud, shows that depositional and accretional growth can be enhanced in KH waves, resulting in more precipitation when compared to a baseline simulation. While sublimation and evaporation occur in KH downdrafts, persistent supersaturation with respect to ice allows for net increase in ice mass. These modeling results and observations suggest that KH waves embedded in mixed-phase stratiform clouds may increase precipitation, although the quantitative impact remains uncertain.

Information about the characteristics of ice particles in clouds is necessary for improving our understanding of the states, processes, and subsequent modeling of clouds and precipitation for numerical weather prediction and climate analysis. Two NASA passive microwave radiometers, the satellite-borne Global Precipitation Measurement (GPM) Microwave Imager (GMI) and the aircraft-borne Conical Scanning Millimeter-Wave Imaging Radiometer (CoSMIR), measure vertically and horizontally polarized microwaves emitted by clouds (including precipitating particles) and Earth’s surface below. In this paper, GMI (or CoSMIR) data are analyzed with CloudSat (or aircraft-borne radar) data to find polarized difference (PD) signals not affected by the surface, thereby obtaining the information on ice particles. Statistical analysis of 4 years of GMI and CloudSat data, for the first time, reveals that optically thick clouds contribute positively to GMI PD at 166 GHz over all the latitudes and their positive magnitude of 166-GHz GMI PD varies little with latitude. This result suggests that horizontally oriented ice crystals in thick clouds are common from the tropics to high latitudes, which contrasts the result of Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) that horizontally oriented ice crystals are rare in optically thin ice clouds.

The National Aeronautics and Space Administration (NASA) has always included data reprocessing as a major component of every science mission. A final reprocessing is typically a part of mission closeout (known as phase F). The Tropical Rainfall Measuring Mission (TRMM) is currently in phase F, and NASA is preparing for the last reprocessing of all the TRMM precipitation data as part of the closeout. This reprocessing includes improvements in calibration of both the TRMM Microwave Imager (TMI) and the TRMM Precipitation Radar (PR). An initial step in the version 8 reprocessing is the improvement of geolocation. The PR calibration is being updated by the Japan Aerospace Exploration Agency (JAXA) using data collected as part of the calibration of the Dual-Frequency Precipitation Radar (DPR) on the Global Precipitation Measurement (GPM) Core Observatory. JAXA undertook a major effort to ensure TRMM PR and GPM Ku-band calibration is consistent.A major component of the TRMM version 8 reprocessing is to create consistent retrievals with the GPM version 05 (V05) retrievals. To this end, the TRMM version 8 reprocessing uses retrieval algorithms based on the GPM V05 algorithms. This approach ensures consistent retrievals from December 1997 (the beginning of TRMM) through the current ongoing GPM retrievals. An outcome of this reprocessing is the incorporation of TRMM data products into the GPM data suite. Incorporation also means that GPM file naming conventions and reprocessed TRMM data carry the V05 data product version. This paper describes the TRMM version 8 reprocessing, focusing on the improvements in TMI level 1 products.

The constellation of spaceborne passive microwave (MW) sensors, coordinated under the framework of the Precipitation Measurement Missions international agreement, continuously produces observations of clouds and precipitation all over the globe. The Goddard profiling algorithm (GPROF) is designed to infer the instantaneous surface precipitation rate from the measured MW radiances. The last version of the algorithm (GPROF-2014)—the product of more than 20 years of algorithmic development, validation, and improvement—is currently used to estimate precipitation rates from the microwave imager GMI on board the GPM core satellite. The previous version of the algorithm (GPROF-2010) was used with the microwave imager TMI on board TRMM. In this paper, TMI-GPROF-2010 estimates and GMI-GPROF-2014 estimates are compared with coincident active measurements from the Precipitation Radar on board TRMM and the Dual-Frequency Precipitation Radar on board GPM, considered as reference products. The objective is to assess the improvement of the GPM-era microwave estimates relative to the TRMM-era estimates and diagnose regions where continuous improvement is needed. The assessment is oriented toward estimating the “effective resolution” of the MW estimates, that is, the finest scale at which the retrieval is able to accurately reproduce the spatial variability of precipitation. A wavelet-based multiscale decomposition of the radar and passive microwave precipitation fields is used to formally define and assess the effective resolution. It is found that the GPM-era MW retrieval can resolve finer-scale spatial variability over oceans than the TRMM-era retrieval. Over land, significant challenges exist, and this analysis provides useful diagnostics and a benchmark against which future retrieval algorithm improvement can be assessed.

This paper presents a novel method for the automatic determination of the diagonal-loading level for robust adaptive beamforming on radar wind profilers. This method balances the degradation of the signal-to-interference ratio with that of the signal-to-noise ratio to maximize the detectability of the backscattered signals. Because radial wind velocities are usually estimated from the first moment of the spectrum of backscattered echoes, both the residual ground clutter and any increase in noise level degrade the detectability of atmospheric echoes. The proposed algorithm evaluates the power spectral density of the residual clutter and increased noise to determine the optimal diagonal-loading level by balancing these two factors. The results of numerical simulation show that, without the need to specify any user parameters, the proposed algorithm is stable and more effective at maximizing the signal-to-interference ratio than the conventional norm-constrained diagonal-loading approach. The stability and clutter suppression capability of the proposed algorithm are examined using data from the Program of the Antarctic Syowa Mesosphere–Stratosphere–Troposphere/Incoherent Scatter Radar.

The refractivity from clutter (RFC) technique has been proved to be an effective way to estimate atmospheric duct structure. An important issue for RFC is how to make the estimate more robust, especially in range-dependent ducting conditions. Traditionally, statistical inversion methods need a large number of forward propagation model runs to obtain an acceptable result. Especially when the parameter search space is multidimensional, these methods are prone to being trapped into local optimal solutions. Recently published results (Zhao and Huang) indicate that the adjoint parabolic equation (PE) method holds promise for real-time estimation of one-dimensional refractive index structure from radar sea clutter returns. This paper is aimed at extending the adjoint PE method to range-dependent evaporation duct cases, with a log-linear relationship describing duct structures. Numerical simulations are used to test the performance of this method and the results are compared with that retrieved using a genetic algorithm. Both noise-free and 3-dB additive Gaussian noise clutter simulations are considered, as well as linearly and nonlinearly varying duct height with range.

The U.S. Weather Surveillance Radar-1988 Doppler (WSR-88D) network has provided meteorologists and hydrologists with quantitative precipitation observations at an unprecedented high spatial–temporal resolution since its deployment in the mid-1990s. Since each single radar can only cover a maximum range of 460 km, a mosaic of multiple-radar observations is needed to generate any national-scale products. The Multi-Radar Multi-Sensor (MRMS) system utilizes a physically based two-dimensional mosaicking algorithm of the WSR-88D data to generate seamless national quantitative precipitation estimation (QPE) products. For areas covered by multiple radars, the mosaicking scheme first determines if precipitation is present by checking the lowest-altitude observation. If the lowest observed radar data indicate no precipitation, then the mosaicked value is set to no precipitation. Otherwise, a weighted mean of multiple-radar observations is taken as the mosaicked value. The weighting function is based on multiple factors, including the distance from the radar and the height of the observation with respect to the melting layer. The mosaic algorithm uses the physically lowest radar observations with no/little blockage while maintaining a spatial continuity in the mosaicked field. The performance of the MRMS seamless radar mosaic algorithm was examined for various precipitation events of different characteristics. The results of these case evaluations are presented in this paper.

This paper presents an adaptive scan strategy concept for phased array weather radars (PAWR) with the objective of increasing the scan speed and capturing features of the storm system while maintaining the measurement accuracy. The adaptive scan strategy is developed based on the space–time variability of the storm under observation. Quickly evolving regions are scanned more often and the spatial sampling resolution is matched to the spatial scale. A model that includes the interaction between space and time is used to extract spatial and temporal scales of the medium and to define scanning regions. The temporal scale constrains the radar revisit time, while the measurement accuracy controls the radar’s dwell time. These conditions are employed in a task scheduler that works on a ray-by-ray basis and is designed to balance task priority and radar resources. The scheduler algorithm also includes an optimization procedure for minimizing radar scan time. The model and the scan strategy are demonstrated using simulation data. The results show that the proposed scan strategy can reduce the scan time significantly without compromising data quality.

This paper describes a method of processing marine radar signals for the purpose of generating phase-resolved surface elevation maps as well as statistical measures of ocean surface wave fields. The method is well suited to the processing of data collected by marine radars because it allows for the incorporation of effects dependent on the radar look direction relative to the propagation direction of ocean waves. Applications to Doppler radar and backscattered power measurements are described, and example results are presented using simulated radar data.

The newly developed High-Performance Instrumented Airborne Platform for Environmental Research (HIAPER) Cloud Radar (HCR) is an airborne, W-band, dual-polarization, Doppler research radar that fits within an underwing pod on the National Center for Atmospheric Research Gulfstream-V HIAPER aircraft. On 2 February 2015, the HCR was flown on its maiden research voyage over a cyclone along the Northeast coast of the United States. Six straight flight legs were flown over 6 h between the northern tip of Delaware Bay and Bangor, Maine, crossing the rain–snow line, and passing directly over Boston, Massachusetts, which received over 16 in. of snow during the event. The HCR, which recorded reflectivity, radial velocity, spectral width, and linear depolarization ratio with a 0.7° beam, was pointed at nadir from a flight altitude of 12,800 m (42,000 ft). The along-track resolution ranged between 20 and 200 m, depending on range, at aircraft speeds varying between 200 and 275 m s−1. The range resolution was 19.2 m. Remarkably detailed finescale structures were found throughout the storm system, including cloud-top generating cells, upright elevated convection, layers of turbulence, vertical velocity perturbations across the melting level, gravity waves, boundary layer circulations, and other complex features. Vertical velocities in these features ranged from 1 to 5 m s−1, and many features were on scales of 5 km or less. The purpose of this paper is to introduce the HCR and highlight the remarkable finescale structures revealed within this Northeast U.S. cyclone by the HCR.

This paper presents a novel, simple method to correct reflectivity measurements of weather radars that operate in attenuation-influenced frequency bands using observations from less attenuated radar systems. In recent years radar systems operating in the X-band frequency range have been developed to provide precipitation fields for areas of special interest in high temporal (≤1 min) and spatial (≤250 m) resolution in complement to nationwide radar networks. However, X-band radars are highly influenced by attenuation. C- and S-band radars typically have coarser resolution (250 m–1 km and 5 min) but are less affected by attenuation.Correcting for attenuation effects in simple (non-Doppler) single-polarized X-band radars remains challenging and is often dependent on restriction parameters, for example, those derived from mountain returns. Therefore, these algorithms are applicable only in limited areas. The method proposed here uses measurements from C-band radars and hence can be applied in all regions covered by nationwide C- (or S-) band radar networks. First, a single scan of X-band radar measurements is used exemplary to identify advantages and disadvantages of the novel algorithm compared to a standard single radar algorithm. The performance of the correction algorithms in different types of precipitation is examined in nine case studies. The proposed method provides very promising results for each type of precipitation. Additionally, it is evaluated in a 5-month comparison with Micro Rain Radar (MRR) observations. The bias between uncorrected X-band radar and MRR data is nearly eliminated by the attenuation correction algorithm, and the RMSE is reduced by 20% while the correlation of ~0.9 between both systems remains nearly constant.

In this paper, the operational Global Precipitation Measurement (GPM) mission combined radar–radiometer algorithm is thoroughly described. The operational combined algorithm is designed to reduce uncertainties in GPM Core Observatory precipitation estimates by effectively integrating complementary information from the GPM Dual-Frequency Precipitation Radar (DPR) and the GPM Microwave Imager (GMI) into an optimal, physically consistent precipitation product. Although similar in many respects to previously developed combined algorithms, the GPM combined algorithm has several unique features that are specifically designed to meet the GPM objectives of deriving, based on GPM Core Observatory information, accurate and physically consistent precipitation estimates from multiple spaceborne instruments, and ancillary environmental data from reanalyses. The algorithm features an optimal estimation framework based on a statistical formulation of the Gauss–Newton method, a parameterization for the nonuniform distribution of precipitation within the radar fields of view, a methodology to detect and account for multiple scattering in Ka-band DPR observations, and a statistical deconvolution technique that allows for an efficient sequential incorporation of radiometer information into DPR precipitation retrievals.

The Multi-Radar Multi-Sensor (MRMS) system, which was developed at the National Severe Storms Laboratory and the University of Oklahoma, was made operational in 2014 at the National Centers for Environmental Prediction. The MRMS system consists of the Warning Decision Support System–Integrated Information suite of severe weather and aviation products, and the quantitative precipitation estimation products created by the National Mosaic and Multi-sensor Quantitative Precipitation Estimation system. Products created by the MRMS system are at a spatial resolution of approximately 1 km, with 33 vertical levels, updating every 2 min over the conterminous United States and southern Canada. This paper describes initial operating capabilities for the severe weather and aviation products that include a three-dimensional mosaic of reflectivity; guidance for hail, tornado, and lightning hazards; and nowcasts of storm location, height, and intensity.

By fitting a parametric vortex model directly to aliased radar radial velocities scanned from a hurricane, the maximum tangential velocity and its radial distance from the hurricane vortex center can be estimated by the recently developed alias-robust vortex analysis. This vortex analysis can be refined to produce a suitable reference radial velocity field on each tilt of a radar scan for the reference check in the first main step of dealiasing. This paper presents the techniques developed for the refinements and shows how and to what extent the refined vortex analysis can improve the reference check and thus enhance the dealiased data coverage, especially over severely aliased data areas around the vortex core of a hurricane or typhoon. In addition, stringent threshold conditions are used in the reference check and the subsequent continuity check to ensure the accepted data are free of alias or almost so. The robustness and improved performance of the method are exemplified by the results tested with severely aliased radial velocities scanned by operational WSR-88D radars from hurricanes in the United States and by operational China New Generation Weather Radar (CINRAD) base data format SA radars from typhoons in China.

The Global Precipitation Measurement (GPM) Dual-Frequency Precipitation Radar (DPR) algorithms consist of modules. This paper describes version 4 (V4) of GPM DPR level 2 (L2) classification (CSF) modules, which consist of two single-frequency (SF) modules—that is, Ku-only and Ka-only modules—and a dual-frequency (DF) module. Each CSF module detects bright band (BB) and classifies rain into three major types, that is, stratiform, convective, and other. The Ku-only and Ka-only CSF modules use algorithms that are similar to the Tropical Rainfall Measuring Mission (TRMM) rain type classification algorithm 2A23. The DF CSF module uses a new method called the measured dual-frequency ratio (DFRm) method for the rain type classification and the detection of BB. It is shown that the Ku-only CSF module and the DF CSF module produce almost indistinguishable rain type counts in a statistical sense. It is also shown that the DFRm method in the DF CSF module improves the detection of BB.

To detect wake vortices in all weather conditions, lidar and radar sensors are complementary. It is important to determine accurately the limitations of the radar detection in clear air and to determine the parameters influencing the detection. This paper is the first study to present the simulation and analysis of the radar signatures of wake vortices on the basis of state-of-the-art 3D large-eddy simulations. The setup of the large-eddy simulations is described and the evolution of the wake vortices is illustrated for the three simulated cases. A modified version of the model developed by Muschinski et al. is used here for the calculation of the radar-backscattered signal and the second moment of the Doppler spectrum that is used to estimate the wake vortex circulation. The simulator uses two different radar configurations for the retrieval of the circulation, which is compared with the circulation obtained directly from the large-eddy simulations. The retrieval of the circulation is fairly accurate if the level of backscattered signal is high enough above the detection threshold of the radar.

Cumulus parameterizations in general circulation models (GCMs) frequently apply mass-flux schemes in their description of tropical convection. Mass flux constitutes the product of the fractional area covered by cumulus clouds in a model grid box and the vertical velocity within the cumulus clouds. The cumulus area fraction profiles can be derived from precipitating radar reflectivity volumes. However, the vertical velocities are difficult to observe, making the evaluation of mass-flux schemes difficult. In this paper, the authors develop and evaluate a parameterization of vertical velocity in convective (cumulus) clouds using only radar reflectivities collected by a C-band polarimetric research radar (CPOL), operating at Darwin, Australia. The parameterization is trained using vertical velocity retrievals from a dual-frequency wind profiler pair located within the field of view of CPOL. The parametric model uses two inputs derived from CPOL reflectivities: the 0-dBZ echo-top height (0-dBZ ETH) and a height-weighted column reflectivity index (ZHWT). The 0-dBZ ETH determines the shape of the vertical velocity profile, while ZHWT determines its strength. The evaluation of these parameterized vertical velocities using (i) the training dataset, (ii) an independent wind-profiler-based dataset, and (iii) 1 month of dual-Doppler vertical velocity retrievals indicates that the statistical representation of vertical velocity is reasonably accurate up to the 75th percentile. However, the parametric model underestimates the extreme velocities. The method allows for the derivation of cumulus mass flux and its variability on current GCM scales based only on reflectivities from precipitating radar, which could be valuable to modelers.

The sun is known to be a good target for weather radar calibration. In this paper high-resolution raster scans of the sun at high elevations will be used to derive the antenna pattern of weather radar, without being affected by beam propagation effects and reflections close to the earth’s surface. It is shown that this pattern matches well to pattern measurements using a point source. Hence, a good estimation of the real antenna pattern can be derived using the sun.Furthermore, formulas to extract undistorted antenna patterns from the sun, even at high-elevation angles, are derived. The signal processing required to achieve high sensitivity for the antenna pattern measurements will be described. Important parts of the antenna pattern—for example, sidelobes—become visible when using long integration times.The polarimetric receiver channel cross-correlation coefficient is proposed as a figure of merit of the cross-polar isolation of the antenna and hence the cross-polar pattern. The results are also compared to point source measurements. This illustrates how an unpolarized signal source like the sun can be used to derive polarimetric variables.

This study is the third part of a paper series investigating the polarimetric radar properties of melting hail and application of those properties for operational polarimetric hail detection and determination of its size. The results of theoretical simulations in Part I were used to develop a hail size discrimination algorithm (HSDA) described in Part II. The HSDA uses radar reflectivity Z, differential reflectivity ZDR, and cross-correlation coefficient ρhv along with melting-level height within a fuzzy-logic scheme to distinguish among three hail size classes: small hail (with diameter D < 2.5 cm), large hail (2.5 < D < 5.0 cm), and giant hail (D > 5.0 cm). The HSDA validation is performed using radar data collected by numerous WSR-88D sites and more than 3000 surface hail reports obtained from the Severe Hazards Analysis and Verification Experiment (SHAVE). The original HSDA version was modified in the process of validation, and the modified algorithm demonstrates probability of detection of 0.594, false-alarm ratio of 0.136, and resulting critical success index (CSI) equal to 0.543. The HSDA outperformed the current operational single-polarization hail detection algorithm, which only provides a single hail size estimate per storm and is characterized by CSI equal to 0.324. It is shown that HSDA is particularly sensitive to the quality of ZDR measurements, which might be affected by possible radar miscalibration and anomalously high differential attenuation.

Detection of meteorological radar signals is often carried out using power averaging with noise subtraction either in the time domain or the spectral domain. This paper considers the relative signal processing gain of these two methods, showing a clear advantage for spectral-domain processing when normalized spectral width is less than ~0.1. A simple expression for the optimal discrete Fourier transform (DFT) length to maximize signal processing gain is presented that depends only on the normalized spectral width and the time-domain weighting function. The relative signal processing gain between noncoherent power averaging and spectral processing is found to depend on a variety of parameters, including the radar wavelength, spectral width, available observation time, and the false alarm rate. Expressions presented for the probability of detection for noncoherent and spectral-based processing also depend on these same parameters. Results of this analysis show that DFT-based processing can provide a substantial advantage in signal processing gain and probability of detection, especially when the normalized spectral width is small and when a large number of samples are available. Noncoherent power estimation can provide superior probability of detection when the normalized spectral width is greater than ~0.1, especially when the desired false alarm rate exceeds 10%.

Rapid advancements of computer technologies in recent years made the real-time transferring and integration of high-volume, multisource data at a centralized location a possibility. The Multi-Radar Multi-Sensor (MRMS) system recently implemented at the National Centers for Environmental Prediction demonstrates such capabilities by integrating about 180 operational weather radars from the conterminous United States and Canada into a seamless national 3D radar mosaic with very high spatial (1 km) and temporal (2 min) resolution. The radar data can be integrated with high-resolution numerical weather prediction model data, satellite data, and lightning and rain gauge observations to generate a suite of severe weather and quantitative precipitation estimation (QPE) products. This paper provides an overview of the initial operating capabilities of MRMS QPE products.

A raindrop size distribution (DSD) retrieval method for a weather radar network consisting of several X-band dual-polarization radars is proposed. An iterative maximum likelihood (ML) estimator for DSD retrieval in a single radar was developed in the authors’ previous work, and the proposed algorithm in this paper extends the single-radar retrieval to radar-networked retrieval, where ML solutions in each single-radar node are integrated based on a Bayesian scheme in order to reduce estimation errors and to enhance accuracy. Statistical evaluations of the proposed algorithm were carried out using numerical simulations. The results with eight radar nodes showed that the bias and standard errors are −0.05 and 0.09 in log(Nw); and Nw (mm−1 m−3) and 0.04 and 0.09 in D0 (mm) in an environment with fluctuations in dual-polarization radar measurements (normal distributions with standard deviations of 0.8 dBZ, 0.2 dB, and 1.5° in ZHm, ZDRm, and ΦDPm, respectively). Further error analyses indicated that the estimation accuracy depended on the number of radar nodes, the ranges of varying μ, the raindrop axis ratio model, and the system bias errors in dual-polarization radar measurements.