Polarimetric X-band Research Articles

Attenuation correction for a high-resolution polarimetric X-band weather radar

Abstract. In 2007, IRCTR (Delft University of Technology) installed a new polarimetric X-band LFMCW radar (IDRA) at the meteorological observation site of Cabauw, The Netherlands. It provides plan position indicators (PPI) at a fixed elevation with a high range resolution of either 3 m or 30 m at a maximum observation range of 1.5 km and 15 km, respectively. IDRA aims to monitor precipitation events for the long-term analysis of the hydrological cycle. Due to the specifications of IDRA, the spatial and temporal variability of a large range of rainfall intensities (from drizzle to heavy convective rain) can be studied. Even though the usual observation range of IDRA is limited to 15 km, attenuation due to precipitation can be large enough to seriously affect the measurements. In this contribution we evaluate the application of a combined method to correct for the specific and the differential attenuation, and in the same vein estimate the parameters of the raindrop-size distribution. The estimated attenuations are compared to a phase constraint attenuation correction method.

Simultaneous Observations and Analysis of Severe Storms Using Polarimetric X-Band SAR and Ground-Based Weather Radar

Recent advances in synthetic aperture radar (SAR) technology have revived meteorological applications with this type of radar. SARs are designed for surface imaging, but now that several X-band multipolarization SAR satellites are in orbit, the attenuation and backscatter caused by precipitation can be better studied. The results presented here demonstrate some of the possibilities by analyzing observations from dual-polarization (HH, VV) TerraSAR-X (TSX) acquisitions over central Florida surrounding severe storms in August 2008. Simultaneous to the SAR acquisitions, WSR-88D ground weather radars in Melbourne and Tampa Bay, FL, collected reflectivity and radial velocity data; the observed strong precipitation cells from convective storms are colocated with severe attenuation in the corresponding SAR images. The observations from SAR measurements are explained quantitatively by converting ground radar reflectivity into spaceborne radar attenuation via a theoretical model. In addition, polarization analysis comparing the SAR image to two additional TSX acquisitions 11 days apart and without rain provides an indication of storm-induced propagation effects on X-band SAR. Specifically, the copolar ratio Z_dr and the copolar correlation differences exhibit behavior that is better explained by the precipitation impact versus surface changes. Multiple regions with varying ground cover, including urban, and storm characteristics are analyzed to highlight the complexity of meteorological research using SAR while revealing a potential use of the technology to investigate the storm structure.

Performance evaluation of high-resolution rainfall estimation by X-band dual-polarization radar for flash flood applications in mountainous basins

Different relations between surface rainfall rate, R , and high-resolution polarimetric X-band radar observations were evaluated using a dense network of rain gauge measurements over complex terrain in Central Italian Alps. The specific differential phase shift, K DP , rainfall algorithm ( R KDP ) although associated with low systematic error it exhibits low sensitivity to the spatial variability of rainfall as compared to the standard algorithm ( R STD ) that is based on the reflectivity-to-rainfall ( Z – R ) relationship. On the other hand, the dependence of the reflectivity measurement on the absolute radar calibration and the rain-path radar signal attenuation introduces significant systematic error on the R STD rainfall estimates. The study shows that adjusting the Z – R relationship for mean-field bias determined using the R KDP estimates as reference is the best technique for acquiring unbiased radar-rainfall estimates at fine space–time scales. Overall, the bias of the R KDP -adjusted Z – R estimator is shown to be lower than 10% for both storm cases, while the relative root-mean-square error is shown to range from 0.6 (convective storm) to 0.9 (stratiform storm). A vertical rainfall profile correction (VPR) technique is tested in this study for the stratiform storm case. The method is based on a newly developed VPR algorithm that uses the X-band polarimetric information to identify the properties of the melting layer and devices a precipitation profile that varies for each radar volume scan to correct the radar-rainfall estimates. Overall, when accounting for the VPR effect there is up to 70% reduction in the systematic error of the 3° elevation estimates, while the reduction in terms of relative root-mean-square error is limited to within 10%.

Bistatic SAR experiment with the Ingara imaging radar

Synthetic aperture radar (SAR) operation in a bistatic configuration not only offers various advantages over its now well-established monostatic counterpart but also poses various challenges. As part of a research programme into the potential benefits and challenges of bistatic SAR, the Ingara fully polarimetric X-band airborne imaging radar system, developed and operated by the Australian Defence Science and Technology Organisation, was upgraded to conduct experimental SAR data collections in a bistatic geometry. Initial trials of the new bistatic SAR system were conducted between December 2007 and April 2008: these involved operation of the existing airborne radar in a fine-resolution (600 MHz bandwidth) circular spotlight-SAR mode, in conjunction with a newly developed fully polarimetric stationary ground-based bistatic receiver. These experimental trials produced a set of fully polarimetric simultaneously collected monostatic and bistatic SAR data, collected over a wide range of bistatic angles. Results from a preliminary analysis of the data have been encouraging: focussed fine-resolution imagery has been obtained, indicating the successful maintenance of synchronisation and phase stability between the independent airborne and ground-based systems. Furthermore, interferometric coherence has been demonstrated between single-pass simultaneously collected monostatic and bistatic images from the airborne and ground-based receivers, and between repeat-pass bistatic images from the ground-based receiver collected some 100 min apart. This study gives an overview of the Ingara bistatic SAR system, discusses the experiments and data processing and presents initial experimental results.

Using CASA IP1 to Diagnose Kinematic and Microphysical Interactions in a Convective Storm

Abstract Data from the Collaborative Adaptive Sensing of the Atmosphere (CASA) Integrated Project I (IP1) network of polarimetric X-band radars are used to observe a convective storm. A fuzzy logic hydrometeor identification algorithm is employed to study microphysical processes. Dual-Doppler techniques are used to analyze the 3D wind field. The scanning strategy, sensitivity, and low-level scanning focus of the radars are investigated for influencing bulk hydrometeor identification and dual-Doppler wind retrievals. Comparisons are made with the nearby S-band polarimetric Next Generation Weather Radar (NEXRAD) prototype radar (KOUN), for consistency. Lightning data are used as an independent indicator of storm evolution for comparison with radar observations. A new methodology for retrieving the vertical wind utilizing upward and variational integration techniques is employed and shown to illustrate trends in mean wind, with particularly good results at low levels. IP1 observations of a case on 10 June 2007 show the development of the updraft, subsequent graupel echo volume evolution, and a descending downdraft preceded by significant graupel in the midlevels, with updraft and graupel volumes leading the onset of lightning. Many of these trends are corroborated by KOUN. The high temporal resolution of three minutes and near-ground sampling provided by IP1 is integral to resolving up- and downdrafts, as well as hydrometeor evolution. IP1 coverage of the upper levels is diminished compared to KOUN, impacting the quality of the dual-Doppler derived vertical winds and ice echo volumes, although the low-level coverage helps to mitigate some errors. However, IP1 coverage of the low- to midlevels is demonstrated to be comparable or better than coverage by KOUN for this storm location.

Evaluating Polarimetric X-Band Radar Rainfall Estimators during HMT

Abstract Different relations between rainfall rate R and polarimetric X-band radar measurables were evaluated using the radar, disdrometer, and rain gauge measurements conducted during the 4-month-long field experiment. The specific differential phase shift KDP–based estimators generally show less scatter resulting from variability in raindrop size distributions than with the power-based relations. These estimators depend on model assumptions about the drop aspect ratios and are not applicable for lighter rainfalls. The polynomial approximation for the mean drop aspect ratio provides R–KDP relations that result overall in good agreement between the radar retrievals of rainfall accumulations and estimates from surface rain gauges. The accumulation data obtained from power estimators that use reflectivity Zeh and differential reflectivity ZDR measurements generally exhibit greater standard deviations with respect to the gauge measurements. Unlike the phase-based estimators, the power-based estimators have an advantage of being “point” measurements, thus providing continuous quantitative precipitation estimation (QPE) for the whole area of radar coverage. The uncertainty in the drop shape model can result in errors in the attenuation and differential attenuation correction procedures. These errors might provide biases of radar-derived QPE for the estimators that use power measurements. Overall, for all considered estimators, the radar-based total rainfall accumulations showed biases less than 10% (relative to gauges). The standard deviations of radar retrievals were about 23% for the mean Zeh–R relation, 17%–22% for the KDP-based estimators (depending on the drop shape model), and about 20%–32% for different Zeh–ZDR-based estimators. Comparing ZDR-based retrievals of mean mass raindrop size Dm (for Dm > 1 mm) with disdrometer-derived values reveals an about 20%–25% relative standard deviation between these two types of estimates.

Analysis of the KK-distribution with medium grazing angle sea-clutter

Robust maritime surveillance with radar requires an accurate description of the backscatter from the sea. An estimated probability distribution of the backscatter is commonly used to determine the threshold for separating targets from clutter at a given false alarm rate. Data collected at medium to high grazing angles by the Defence Science Technology Organisation (DSTO) Ingara fully polarimetric X-band radar demonstrate that the commonly used K-distribution is not always adequate for modelling the probability distribution. This is especially the case for the horizontal polarisation and in regions of high backscatter where target detection can be a problem. An alternative proposed as a more accurate model in this region is known as the KK-distribution. The analysis presented in this study describes this model with the addition of multiple looks and a thermal noise component to produce greater accuracy in the mean and underlying shape. The threshold required to achieve a constant false alarm rate is then studied and compared with the K-distribution.

Transitions in the polarimetric radar scattering properties of the sea surface

The wide variability of sea surface conditions is reflected in a corresponding variability in sea clutter characteristics; moreover, these characteristics manifest a strong dependence on radar parameters such as frequency, waveform and polarisation. Ideally, radar systems should adapt their parameters to the prevailing sea conditions and viewing geometry so as to optimise performance. In practice this is seldom done, in large part because most radar systems have very limited capability to adapt, but also because the rules that might guide the adaptation process are not yet clearly established. Where an appreciable change in radar parameters or viewing geometry would yield only a slight change in clutter severity and a commensurately modest improvement in radar performance, there is little incentive to adapt. On the other hand, if there are sharp transitions in clutter properties as some radar parameter is varied, then exploitation of the transition may be highly advantageous. This study presents the results of a search for transitions in polarimetric X-band sea clutter, focusing on the dependence on scattering geometry. Experiments were conducted to study clutter variation with elevation angle and azimuth with respect to wind direction. The results demonstrate that transitions exist, that they can be surprisingly sharp, and that the consequences for target detectability may be significant.

Rainfall spatial variability observed by X-band weather radar and its implication for the accuracy of rainfall estimates

The main objective of this paper is to estimate the error in the rainfall derived from a polarimetric X-band radar, by comparison with the corresponding estimate of a rain gauge network. However the present analysis also considers the errors inherent to rain gauge, in particular instrumental and representativeness errors. A special emphasis is addressed to the spatial variability of the rainfall in order to appreciate the representativeness error of the rain gauge with respect to the 1 km square average, typical of the radar derived estimate. For this purpose the spatial correlation function of the rainfall is analyzed. The data set consists of 1-year radar data collected by the X-band polarimetric radar HYDRIX ®, located in Beauce region (80 km south of Paris). All data were processed in real time using the ZPHI ® algorithm. A dense 25 rain gauge network provided ground comparison data. The various sources of uncertainties (instrumental and representativeness) are then analyzed and quantified for each sensor.

Evaluation and Bias Removal of Multilook Effect on Entropy/Alpha/Anisotropy in Polarimetric SAR Decomposition

Entropy, alpha, and anisotropy (H/alpha/A) of the polarimetric target decomposition have been an effective and popular tool for polarimetric synthetic aperture radar (SAR) image analysis and for a geophysical parameter estimation. However, multilook processing can severely affect the values of these parameters. In this paper, a Monte Carlo simulation is used to evaluate and remove the bias generated by the multilook effect on these parameters for various media composed of grassland, forest, and urban returns. Due to insufficient averaging, entropy is underestimated, and anisotropy is overestimated. We also found that the bias in the alpha angle can be either underestimated or overestimated depending on scattering mechanisms. Based on simulation results, efficient bias removal procedures have been developed. In particular, the entropy bias can be precisely corrected, and the amount of correction is independent of the radar frequency and SAR systems. Data from L-band Advanced Land Observing Satellite/phased array type L-band SAR, German Aerospace Research Center (DLR)/enhanced SAR, Jet Propulsion Laboratory (JPL)/airborne SAR, and X-band polarimetric and interferometric SAR are used for demonstration in this paper.

Comparison of Two Raindrop Size Distribution Retrieval Algorithms for X-Band Dual Polarization Observations

Abstract In this study the authors evaluate two algorithms, the so-called beta (β) and constrained methods, proposed for retrieving the governing parameters of the “normalized” gamma drop size distribution (DSD) from dual-polarization radar measurements. The β method treats the drop axis ratio as a variable and computes drop shape and DSD parameters from radar reflectivity (ZH), differential reflectivity (ZDR), and specific differential phase shift (KDP). The constrained method assumes that the axis-ratio relation is fixed and computes DSD parameters from ZH, ZDR, and an empirical relation between the DSD slope and shape parameters. The two techniques are evaluated for polarimetric X-band radar observations by comparing retrieved DSD parameters with disdrometer observations and examining simulated radar parameters for consistency. Error effects on the β method and constrained method retrievals are analyzed. The β approach is found to be sensitive to errors in KDP and to be less consistent with observations. Large retrieved β values are found to be associated with large retrieved DSD shape parameters and small median drop diameters. The constrained method provides reasonable rain DSD retrievals that agree better with disdrometer observations.

Use of a Radar Simulator on the Output Fields from a Numerical Mesoscale Model to Analyze X-Band Rain Estimators

Abstract A full radar simulator, which works with the 3D output fields from a numerical mesoscale model, has been developed. This simulator uses a T-matrix code to calculate synthetic radar measurements, accounts for both backscattering and propagation effects, and includes polarimetric variables. The tool is modular to allow several options in the derivation of the synthetic radar variables. A measurement uncertainty can be taken into account on both the simulated reflectivities and the differential phase shift. A scheme can also be switched on to allow for the gate-to-gate variability of the rain drops size distribution or, also, their oblateness. This work was done in the framework of the installation in West Africa of a polarimetric X-band radar for the observation of tropical rain. Accordingly, the first objective pursued with this simulation setup is a detailed analysis of X-band polarimetric rain retrieval algorithms. Two retrieval schemes, a simple R–KDP formula and a profiler that uses both reflectivity and ϕDP, are tested. For that purpose the simulator is run on a model case study of an African squall line, then the two schemes are used to retrieve the rain rates from the synthetic radar variables and compare them to the original. The scores of the schemes are discussed and compared. The authors analyze the sensitivity of the results to the measurement uncertainty and also to several aspects of drop size distribution and drop shape variability.

Four component decomposition of polarimetric microwave scattering of rice using three dimensional structural data

In the present research, the authors examined the effects of the probability distribution of the target shape to the results of four component decomposition. The four component decomposition method decomposes polarimetric microwave data into four components, i.e. surface scattering, double bounce scattering, volume scattering and helix scattering. Prior to the present research, three dimensional structure of the rice was measured and modeled. The model was used to calculate the probability density funtion of the rice leaves at a viewpoint of arbitrary zenith and azimuth angles. With this probability density function, a covariance matrix of volume scattering for the rice was developed. This covariance matrix was embedded in the four component decomposition, and it was applied for the polarimetric X-band scattering data of rice, measured in an anechoic radio wave chamber in Niigata Univ. during September 25 to 26, 2007. It was found that the proposed four component decomposition produced more volume scattering and less surface scattering than the results produced by the standard four component analysis. As a result, it is concluded that the probability distribution of the target shape should be considered.

A Polarimetric Radar Approach to Identify Rain, Melting-Layer, and Snow Regions for Applying Corrections to Vertical Profiles of Reflectivity

Abstract This article describes polarimetric X-band radar-based quantitative precipitation estimations (QPE) under conditions of low freezing levels when, even at the lowest possible elevation angles, radar resolution volumes at longer ranges are in melting-layer or snow regions while it rains at the ground. A specifically adjusted vertical-profile-of-reflectivity (VPR) approach is introduced. The mean VPR is constructed based on the range–height indicator scans, and the effects of smoothing of brightband (BB) features with range are accounted for. A principal feature of the suggested QPE approach is the determination of the reflectivity BB boundaries and freezing-level heights on a beam-by-beam basis using the copolar correlation coefficient ρhv, which is routinely available from the X-band radar measurements. It is shown that this coefficient provides a robust discrimination among the regions of rain, melting hydrometeors, and snow. The freezing-level estimates made using ρhv were within 100–200 m from the corresponding estimates of the 0° isotherm heights from radiosonde soundings. The suggested VPR approach with the polarimetric determination of the reflectivity BB boundaries was used for QPE during the wintertime deployment of the NOAA X-band radar as part of the 2006 Hydrometeorological Test Bed (HMT-06) field experiment in the California Sierra Nevada foothills. It is shown that this approach noticeably improves radar-rainfall accumulation estimates. The use of the HMT-06 mean X-band reflectivity–rain-rate (Zeh–R) relation resulted in an approximately 65% relative standard deviation of radar estimates from the surface rain gauges if no VPR adjustments were made. Applying the VPR approach with polarimetric detection of the melting layer resulted in reduction of the corresponding relative standard deviation by about a factor of 2.

An interferometric coherence optimization method in radar polarimetry for high-resolution imagery

This paper investigates to what extent a new interferometric coherence optimization in radar polarimetry allows the separation of point scatterers located in the same resolution cell according to their interferometric phases. An interferometric coherence definition called the single-mechanism coherence is introduced, and the corresponding optimization method is briefly discussed. This method was first validated theoretically when no volume decorrelation occurs. Then, it has been applied to simple target measurements acquired in an anechoic chamber, and to an X-band polarimetric and interferometric synthetic aperture radar image containing man-made targets. In both cases, the single-mechanism coherence optimization enables to resolve the interferometric phases of several scattering centers inside the same resolution cell.

The Utility of X-Band Polarimetric Radar for Quantitative Estimates of Rainfall Parameters

Abstract The utility of X-band polarimetric radar for quantitative retrievals of rainfall parameters is analyzed using observations collected along the U.S. west coast near the mouth of the Russian River during the Hydrometeorological Testbed project conducted by NOAA’s Environmental Technology and National Severe Storms Laboratories in December 2003 through March 2004. It is demonstrated that the rain attenuation effects in measurements of reflectivity (Ze) and differential attenuation effects in measurements of differential reflectivity (ZDR) can be efficiently corrected in near–real time using differential phase shift data. A scheme for correcting gaseous attenuation effects that are important at longer ranges is introduced. The use of polarimetric rainfall estimators that utilize specific differential phase and differential reflectivity data often provides results that are superior to estimators that use fixed reflectivity-based relations, even if these relations were derived from the ensemble of drop size distributions collected in a given geographical region. Comparisons of polarimetrically derived rainfall accumulations with data from the high-resolution rain gauges located along the coast indicated deviation between radar and gauge estimates of about 25%. The ZDR measurements corrected for differential attenuation were also used to retrieve median raindrop sizes, D0. Because of uncertainties in differential reflectivity measurements, these retrievals are typically performed only for D0 > 0.75 mm. The D0 estimates from an impact disdrometer located at 25 km from the radar were in good agreement with the radar retrievals. The experience of operating the transportable polarimetric X-band radar in the coastal area that does not have good coverage by the National Weather Service radar network showed the value of such radar in filling the gaps in the network coverage. The NOAA X-band radar was effective in covering an area up to 40–50 km in radius offshore adjacent to a region that is prone to flooding during wintertime landfalling Pacific storms.