Path-integrated Attenuation Research Articles

Assessing rainfall radar errors with an inverse stochastic modelling framework

Abstract. Weather radar is a crucial tool for rainfall observation and forecasting, providing high-resolution estimates in both space and time. Despite this, radar rainfall estimates are subject to many error sources – including attenuation, ground clutter, beam blockage and drop-size distribution – with the true rainfall field unknown. A flexible stochastic model for simulating errors relating to the radar rainfall estimation process is implemented, inverting standard weather radar processing methods and imposing path-integrated attenuation effects, a stochastic drop-size-distribution field, and sampling and random errors. This can provide realistic weather radar images, of which we know the true rainfall field and the corrected “best-guess” rainfall field which would be obtained if they were observed in a real-world case. The structure of these errors is then investigated, with a focus on the frequency and behaviour of “rainfall shadows”. Half of the simulated weather radar images have at least 3 % of their significant rainfall rates shadowed, and 25 % have at least 45 km2 containing rainfall shadows, resulting in underestimation of the potential impacts of flooding. A model framework for investigating the behaviour of errors relating to the radar rainfall estimation process is demonstrated, with the flexible and efficient tool performing well in generating realistic weather radar images visually for a large range of event types.

A unified synergistic retrieval of clouds, aerosols, and precipitation from EarthCARE: the ACM-CAP product

Abstract. ACM-CAP provides a synergistic best-estimate retrieval of all clouds, aerosols, and precipitation detected by the atmospheric lidar (ATLID), cloud-profiling radar (CPR), and multi-spectral imager (MSI) aboard EarthCARE (Earth Cloud Aerosol and Radiation Explorer). While synergistic retrievals are now mature in many contexts, ACM-CAP is unique in that it provides a unified retrieval of all hydrometeors and aerosols. The Cloud, Aerosol, and Precipitation from mulTiple Instruments using a VAriational TEchnique (CAPTIVATE) algorithm allows for a robust accounting of observational and retrieval errors and the contributions of passive and integrated measurements and for enforcing physical relationships between components (e.g. the conservation of precipitating mass flux through the melting layer). We apply ACM-CAP to EarthCARE scenes simulated from numerical weather model forecasts and evaluate the retrievals against “true” quantities from the numerical model. The retrievals are well-constrained by observations from active and passive instruments and overall closely resemble the bulk quantities (e.g. cloud water content, precipitation mass flux, and aerosol extinction) and microphysical properties (e.g. cloud effective radius and median volume diameter) from the model fields. The retrieval performs best where the active instruments have a strong and unambiguous signal, namely in ice clouds and snow, which is observed by both ATLID and CPR, and in light to moderate rain, where the CPR signal is strong. In precipitation, CPR's Doppler capability permits enhanced retrievals of snow particle density and raindrop size. In complex and layered scenes where ATLID is obscured, we have shown that making a simple assumption about the presence and vertical distribution of liquid cloud in rain and mixed-phase clouds allows improved assimilation of MSI solar radiances. In combination with a constraint on the CPR path-integrated attenuation from the ocean surface, this leads to improved retrievals of both liquid cloud and rain in midlatitude stratiform precipitation. In the heaviest precipitation, both active instruments are attenuated and dominated by multiple scattering; in these situations, ACM-CAP provides a seamless retrieval of cloud and precipitation but is subject to a high degree of uncertainty. ACM-CAP's aerosol retrieval is performed in hydrometeor-free parts of the atmosphere and constrained by ATLID and MSI solar radiances. While the aerosol optical depth is well-constrained in the test scenes, there is a high degree of noise in the profiles of extinction. The use of numerical models to simulate test scenes has helped to showcase the capabilities of the ACM-CAP clouds, aerosols, and precipitation product ahead of the launch of EarthCARE.

Processing reflectivity and Doppler velocity from EarthCARE's cloud-profiling radar: the C-FMR, C-CD and C-APC products

Abstract. The Earth Clouds, Aerosols and Radiation (EarthCARE) satellite mission is a joint effort by the European Space Agency (ESA) and the Japanese Aerospace Exploration Agency (JAXA). The EarthCARE mission features the first spaceborne 94 GHz cloud-profiling radar (CPR) with Doppler capability. The raw CPR observations and auxiliary information are used as input to three Level-2 (L2) algorithms: (1) C-APC: Antenna Pointing Characterization; (2) C-FMR: CPR feature mask and reflectivity; (3) C-CD: Corrected CPR Doppler Measurements. These algorithms apply quality control and corrections to the CPR primary measurements and derive important geophysical variables, such as hydrometeor locations, and best estimates of particle sedimentation fall velocities. The C-APC algorithm uses natural targets to introduce any corrections needed to the CPR raw Doppler velocities due to the CPR antenna pointing. The C-FMR product provides the feature mask based on only-reflectivity CPR measurements and quality-controlled radar-reflectivity profiles corrected for gaseous attenuation at 94 GHz. In addition, C-FMR provides best estimates of the path-integrated attenuation (PIA) and flags identifying the presence of multiple scattering in the CPR observations. Finally, the C-CD product provides the quality-controlled, bias-corrected mean Doppler velocity estimates (Doppler measurements corrected for antenna mispointing, non-uniform beam filling and velocity folding). In addition, the best estimate of the particle sedimentation velocity is estimated using a novel technique.

Understanding Errors in Cloud Liquid Water Path Retrievals Derived from CloudSat Path-Integrated Attenuation

Abstract An algorithm that derives the nonprecipitating cloud liquid water path Wcld from CloudSat using a surface reference technique (SRT) is presented. The uncertainty characteristics of the SRT are evaluated. It is demonstrated that an accurate analytical formulation for the pixel-scale precision can be derived. The average precision of the SRT is estimated to be 34 g m−2 at the individual pixel scale; however, precision systematically decreases from around 30 to 40 g m−2 as cloud fraction varies from 0% to 100%. The retrievals of clear-sky Wcld have a mean bias of 0.9 g m−2. Output from a large-eddy simulation coupled to a radar simulator shows that an additional bias of −8% may result from nonuniformity within the footprint of cloudy pixels. The retrieval yield for the SRT, measured relative to all warm clouds over ocean between 60°N and 60°S latitude is 43%. The SRT Wcld is compared with one estimate of Wcld from the Moderate Resolution Imaging Spectroradiometer (MODIS) using an adiabatic cloud profile and an effective radius derived from 3.7-μm reflectance. A strong correlation between the mean MODIS Wcld and SRT Wcld is found across diverse cloud regimes, but with biases in the mean Wcld that are cloud-regime dependent. Overall, the mean bias of the SRT relative to MODIS is −13.1 g m−2. Systematic underestimates of Wcld by the SRT resulting from nonuniform beamfilling cannot be ruled out as an explanation for the retrieval bias.

Correction of Path-Integrated Attenuation Estimates Considering the Soil Moisture Effect for the GPM Dual-Frequency Precipitation Radar

Abstract Spaceborne precipitation radars, including the Tropical Rainfall Measuring Mission’s Precipitation Radar (PR) and the Global Precipitation Measurement Mission’s Dual-Frequency Precipitation Radar (DPR), measure not only precipitation echoes but surface echoes as well, the latter of which are used to estimate the path-integrated attenuation (PIA) in the surface reference technique (SRT). In our previous study based on analyzing PR measurements, we found that attenuation-free surface backscattering cross sections (denoted by ) over land increased in the presence of precipitation. This behavior, called the soil moisture effect, causes an underestimate of the PIA by the SRT as the method does not explicitly consider this effect. In this study, measurements made by Ku-band Precipitation Radar (KuPR) and Ka-band Precipitation Radar (KaPR), which comprise the DPR, were analyzed to examine whether KuPR and KaPR exhibit similar dependencies on the soil moisture as does the PR. For both KuPR and KaPR, an increase in was observed for a large portion of the land area, except for forests and deserts. Results from the Hitschfeld–Bordan (HB) method suggest that increases with the surface precipitation rate for light precipitation events. Meanwhile, for heavy precipitation, owing to the degradation of the HB method, it is difficult to estimate quantitatively. Thus, a correction method for PIA that considers the soil moisture effect was developed and implemented into the DPR standard algorithm. With this correction, the surface precipitation rate estimates increased by approximately 18% for KuPR and 15% for the normal scan of KaPR over land.

Improving Millimeter Radar Attenuation Corrections in High-Latitude Mixed-Phase Clouds via Radio Soundings and a Suite of Active and Passive Instruments

Supercooled liquid clouds are very frequent in high-latitude regions. In addition to their substantial effect on visible and infrared radiation, they affect the signal of millimeter radars by producing nonnegligible attenuation. Such attenuation must be properly corrected if the information of millimeter radars is used in quantitative retrievals for inferring ice microphysical properties. This study proposes a multisensor scheme for refining the vertical distribution of supercooled liquid water content (SLWC) compared to state-of-the-art methods that equipartition the liquid water path measured by microwave radiometer to all pixels identified as cloudy by the radars and warmer than −40 °C. Our methodology is applicable in high-latitude, mixed-phase environments based on the synergy between radar and lidar binary cloud phase masking, microwave radiometer, and radio sounding observations. The technique is demonstrated via data collected by the U.S. Department of Energy (DoE) Atmospheric Radiation Measurement (ARM) Program climate research facility at the North Slope of Alaska (NSA) and compared with the state-of-the-art methods. Path integrated attenuation (PIA) at W- and G-band frequencies (>95 GHz) is then assessed. Results indicate that the different in-cloud distributions of the liquid condensate lead to round-trip PIA discrepancies of cloudy volumes that range in [2, 5] dB at W- and G-band frequencies. These differences far exceed those encountered when changing some of the algorithm’s arbitrary assumptions and weighting functions.

Estimating Mean Surface Backscatter from GPM Surface Backscatter Observations

The radar on the Global Precipitation Measurement (GPM) mission observes precipitation at 13.6 GHz (Ku-band) and 35.6 GHz (Ka-band) and also receives echoes from the earth’s surface. Statistics of surface measurements for non-raining conditions are saved in a database for later use in estimating the precipitation path-integrated attenuation. Previous work by Meneghini and Jones (2011) showed that while averaging over larger latitude/longitude bins increase the number of samples, it can also increase sample variance due to spatial inhomogeneity in the data. As a result, Meneghini and Kim (2017) proposed a new, adaptive method of database construction, in which the number of measurements averaged depends on the spatial homogeneity. The purpose of this work is to re-visit previous, single-frequency results using dual-frequency data and optimal interpolation (kriging). Results include that (1) temporal inhomogeneity can create similar results as spatial, (2) Ka-band behavior is similar to Ku-band, (3) the Ku-/Ka-band difference has less spatial inhomogeneity than either band by itself, and (4) kriging and the adaptive method can reduce the sample variance. The author concludes that finer spatial and temporal resolution is necessary in constructing the database for single frequencies but less so for the Ku-/Ka-band difference. The adaptive approach reduces sample standard deviation with a relatively modest computational increase.

DATA OPTIMIZATION OF RAIN RADAR AND ITS COVERAGE EXPANSION USING RADAR NETWORK

This research developed network of two Santanu rain radars of the Center for Atmospheric Science and Technology LAPAN installed in Bandung (SANTANU 1) and Sumedang (SANTANU 2). The objective is to get larger coverage area which also comply data each other. The first step is data optimation to each radar detection results, particularly due to rain attenuation, and then validation with the transportable radar image. The shape analysis on a rain object detected in the slice of both radar images indicates 0.8672 of R 2 which show the width of both objectcs has a good similarity. Data optimation due to rain attenuation using attenuation correction for regional weather radar (KA) and path-integrated attenuation (PIA) indicate that the reflectivity improvement by KA is 10 times better than PIA. The object detected by Santanu rain radar is similar to the transportable radar for reflectivity larger than 30 dBZ. The next step is integration of both rain radar images to create mosaic image and then validation using the transportable radar data. The production of the mosaic image is tried out using three methods: image stitching, averaging, and Cressman scheme. The best result is presented by Cressman scheme which show more integrated slice area and maintain the features of the object. However, the processing time of this scheme is 17 minutes which is longer than image stitching method. The processing time depend on the number of grid data and the hardware memory. Validation on the mosaic result by using transportable radar data by implementing structural similarity index (SSIM) method yields 0,32 of average value. The improvement on SSIM value is obtained if speckle removal is applied by recording 59% of enhancement

Rainfall-induced attenuation correction for two operational dual-polarization C-band radars in the Netherlands

AbstractThe Royal Netherlands Meteorological Institute (KNMI) operates two operational dual-polarization C-band weather radars providing 2-D radar rainfall products. Attenuation can result in severe underestimation of rainfall amounts, particularly in convective situations that are known to have high impact on society. In order to improve the radar-based precipitation estimates, two attenuation correction methods are evaluated and compared: 1) Modified Kraemer (MK) method, i.e. Hitschfeld-Bordan where parameters of the power-law Zh–kh relation are adjusted such that reflectivities in the entire dataset do not exceed 59 dBZh and attenuation correction is limited to 10 dB; 2) using two-way path-integrated attenuation computed from the dual-polarization moment specific differential phase Kdp (Kdp method). In both cases the open-source Python library wradlib is employed for the actual attenuation correction. A radar voxel only contributes to the computed path-integrated attenuation if its height is below the forecasted freezing level height from the numerical weather prediction model HARMONIE-AROME. The methods are effective in improving hourly and daily quantitative precipitation estimation (QPE), where the Kdp method performs best. The verification against rain gauge data shows that the underestimation diminishes from 55% to 37% for hourly rainfall for the Kdp method when the gauge indicates more than 10 mm of rain in that hour. The improvement for the MK method is less pronounced, with a resulting underestimation of 40%. The stability of the MK method holds a promise for application to data archives from single-polarization radars.

Cloud Assumption of Precipitation Retrieval Algorithms for the Dual-Frequency Precipitation Radar

AbstractAn assumption related to clouds is one of uncertain factors in precipitation retrievals by the Dual-Frequency Precipitation Radar (DPR) on board the Global Precipitation Measurement (GPM) Core Observatory. While an attenuation due to cloud ice is negligibly small for Ku and Ka bands, attenuation by cloud liquid water is larger in the Ka band and estimating precipitation intensity with high accuracy from Ka-band observations can require developing a method to estimate the attenuation due to cloud liquid water content (CLWC). This paper describes a CLWC database used in the DPR level-2 algorithm for the GPM V06A product. In the algorithm, the CLWC value is assumed using the database with inputs of precipitation-related variables, temperature, and geolocation information. A calculation of the database was made using the 3.5-km-mesh global atmospheric simulation derived from the Nonhydrostatic Icosahedral Atmospheric Model (NICAM) global cloud-system-resolving model. Impacts of current CLWC assumptions for surface precipitation estimates were evaluated by comparisons of precipitation retrieval results between default values and 0 mg m−3 of the CLWC. The impacts were quantified by the normalized mean absolute difference (NMAD) and the NMAD values showed 2.3% for the Ku, 9.9% for the Ka, and 6.5% for the dual-frequency algorithms in global averages, while they were larger in the tropics than in high latitudes. Effects of the precipitation estimates from the CLWC assumption were examined further in terms of retrieval processes affected by the CLWC assumption. This study emphasizes the CLWC assumption provided more effects on the precipitation estimates through estimating path-integrated attenuation due to rain.

An Improvement to Precipitation Inversion Model Using Oblique Earth–Space Link Based on the Melting Layer Attenuation

Based on the power-law relation between rain intensity and corresponding specific attenuation, precipitation inversion by oblique earth–space link (OEL) has become a supplement to conventional method. However, the attenuation caused by melting layer has a significant effect on its performance. According to the microphysics, thermodynamics, and scattering properties of melting hydrometeors, a melting layer model is introduced to improve OEL inversion approach. Meanwhile, the observations by a microrain radar (MRR) located at Nanjing Reference Climatological Station (NRCS) are used and plenty of numerical simulations are also performed, which aim to validate the melting layer model. The choice of raindrop size distribution (DSD) determines the simulated values of radar reflectivity factor (RRF) and path-integrated attenuation (PIA). The gamma distribution gives the results of simulated RRF that are consistent with MRR observations. Considering the nonspherical shape of melting hydrometeors, its effects on attenuation calculation are discussed based on <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$T$ </tex-math></inline-formula> -matrix, which indicates the error due to hydrometeor shape is negligible below Ka-band. The values of PIA derived from our melting layer model show good agreements with MRR observations. In the simulation experiment, we found the average multiplier of melting layer attenuation is a function of rain intensity and frequency. The relation is established by using the nonlinear least-squares method. According to the multiplier, a novel model is proposed to estimate the melting layer attenuation. Compared with other existing attenuation models, it shows a similar performance in different frequency bands. Moreover, we introduce an equivalent multiplier into OEL rainfall inversion model, which demonstrates the obvious enhancement of effectiveness and robustness of improved model in the validation experiment.

What Can We Learn from the CloudSat Radiometric Mode Observations of Snowfall over the Ice-Free Ocean?

The quantification of global snowfall by the current observing system remains challenging, with the CloudSat 94 GHz Cloud Profiling Radar (CPR) providing the current state-of-the-art snow climatology, especially at high latitudes. This work explores the potential of the novel Level-2 CloudSat 94 GHz Brightness Temperature Product (2B-TB94), developed in recent years by processing the noise floor data contained in the 1B-CPR product; the focus of the study is on the characterization of snow systems over the ice-free ocean, which has well constrained emissivity and backscattering properties. When used in combination with the path integrated attenuation (PIA), the radiometric mode can provide crucial information on the presence/amount of supercooled layers and on the contribution of the ice to the total attenuation. Radiative transfer simulations show that the location of the supercooled layers and the snow density are important factors affecting the warming caused by supercooled emission and the cooling induced by ice scattering. Over the ice-free ocean, the inclusion of the 2B-TB94 observations to the standard CPR observables (reflectivity profile and PIA) is recommended, should more sophisticated attenuation corrections be implemented in the snow CloudSat product to mitigate its well-known underestimation at large snowfall rates. Similar approaches will also be applicable to the upcoming EarthCARE mission. The findings of this paper are relevant for the design of future missions targeting precipitation in the polar regions.

Estimating total attenuation using Rayleigh targets at cloud top: applications in multilayer and mixed-phase clouds observed by ground-based multifrequency radars

Abstract. At millimeter wavelengths, attenuation by hydrometeors, such as liquid droplets or large snowflakes, is generally not negligible. When using multifrequency ground-based radar measurements, it is common practice to use the Rayleigh targets at cloud top as a reference in order to derive attenuation-corrected reflectivities and meaningful dual-frequency ratios (DFRs). By capitalizing on this idea, this study describes a new quality-controlled approach that aims at identifying regions of cloud where particle growth is negligible. The core of the method is the identification of a “Rayleigh plateau”, i.e., a large enough region near cloud top where the vertical gradient of DFR remains small. By analyzing co-located Ka–W band radar and microwave radiometer (MWR) observations taken at two European sites under various meteorological conditions, it is shown how the resulting estimates of differential path-integrated attenuation (ΔPIA) can be used to characterize hydrometeor properties. When the ΔPIA is predominantly produced by cloud liquid droplets, this technique alone can provide accurate estimates of the liquid water path. When combined with MWR observations, this methodology paves the way towards profiling the cloud liquid water, quality-flagging the MWR retrieval for rain and drizzle contamination, and/or estimating the snow differential attenuation.

Mind the gap – Part 2: Improving quantitative estimates of cloud and rain water path in oceanic warm rain using spaceborne radars

Abstract. The intrinsic small spatial scales and low-reflectivity structure of oceanic warm precipitating clouds suggest that millimeter spaceborne radars are best suited to providing quantitative estimates of cloud and rain liquid water paths (LWPs). This assertion is based on their smaller horizontal footprint; high sensitivities; and a wide dynamic range of path-integrated attenuations associated with warm-rain cells across the millimeter wavelength spectrum, with diverse spectral responses to rain and cloud partitioning. State-of-the-art single-frequency radar profiling algorithms of warm rain seem to be inadequate because of their dependence on uncertain assumptions about the rain–cloud partitioning and because of the rain microphysics. Here, high-resolution cloud-resolving model simulations for the Rain in Cumulus over the Ocean field study and a spaceborne forward radar simulator are exploited to assess the potential of existing and future spaceborne radar systems for quantitative warm-rain microphysical retrievals. Specifically, the detrimental effects of nonuniform beam filling on estimates of path-integrated attenuation (PIA), the added value of brightness temperature (TB) derived adopting radiometric radar modes, and the performances of multifrequency PIA and/or TB combinations when retrieving liquid water paths partitioned into cloud (c-LWPs) and rain (r-LWPs) are assessed. Results show that (1) Ka- and W-band TB values add useful constraints and are effective at lower LWPs than the same-frequency PIAs; (2) matched-beam combined TB values and PIAs from single-frequency or multifrequency radars can significantly narrow down uncertainties in retrieved cloud and rain liquid water paths; and (3) the configuration including PIAs, TB values and near-surface reflectivities for the Ka-band–W-band pairs in our synthetic retrieval can achieve an RMSE of better than 30 % for c-LWPs and r-LWPs exceeding 100 g m−2.

Preliminary investigation of the relationship between differential phase shift and path-integrated attenuation at the X band frequency in an Alpine environment

Abstract. The RadAlp experiment aims at developing advanced methods for rainfall and snowfall estimation using weather radar remote sensing techniques in high mountain regions for improved water resource assessment and hydrological risk mitigation. A unique observation system has been deployed since 2016 in the Grenoble region of France. It is composed of an X-band radar operated by Météo-France on top of the Moucherotte mountain (1901 m above sea level; hereinafter MOUC radar). In the Grenoble valley (220 m above sea level; hereinafter a.s.l.), we operate a research X-band radar called XPORT and in situ sensors (weather station, rain gauge and disdrometer). In this paper we present a methodology for studying the relationship between the differential phase shift due to propagation in precipitation (Φdp) and path-integrated attenuation (PIA) at X band. This relationship is critical for quantitative precipitation estimation (QPE) based on polarimetry due to severe attenuation effects in rain at the considered frequency. Furthermore, this relationship is still poorly documented in the melting layer (ML) due to the complexity of the hydrometeors' distributions in terms of size, shape and density. The available observation system offers promising features to improve this understanding and to subsequently better process the radar observations in the ML. We use the mountain reference technique (MRT) for direct PIA estimations associated with the decrease in returns from mountain targets during precipitation events. The polarimetric PIA estimations are based on the regularization of the profiles of the total differential phase shift (Ψdp) from which the profiles of the specific differential phase shift on propagation (Kdp) are derived. This is followed by the application of relationships between the specific attenuation (k) and the specific differential phase shift. Such k–Kdp relationships are estimated for rain by using drop size distribution (DSD) measurements available at ground level. Two sets of precipitation events are considered in this preliminary study, namely (i) nine convective cases with high rain rates which allow us to study the ϕdp–PIA relationship in rain, and (ii) a stratiform case with moderate rain rates, for which the melting layer (ML) rose up from about 1000 up to 2500 m a.s.l., where we were able to perform a horizontal scanning of the ML with the MOUC radar and a detailed analysis of the ϕdp–PIA relationship in the various layers of the ML. A common methodology was developed for the two configurations with some specific parameterizations. The various sources of error affecting the two PIA estimators are discussed, namely the stability of the dry weather mountain reference targets, radome attenuation, noise of the total differential phase shift profiles, contamination due to the differential phase shift on backscatter and relevance of the k–Kdp relationship derived from DSD measurements, etc. In the end, the rain case study indicates that the relationship between MRT-derived PIAs and polarimetry-derived PIAs presents an overall coherence but quite a considerable dispersion (explained variance of 0.77). Interestingly, the nonlinear k–Kdp relationship derived from independent DSD measurements yields almost unbiased PIA estimates. For the stratiform case, clear signatures of the MRT-derived PIAs, the corresponding ϕdp value and their ratio are evidenced within the ML. In particular, the averaged PIA∕ϕdp ratio, a proxy for the slope of a linear k–Kdp relationship in the ML, peaks at the level of the copolar correlation coefficient (ρhv) peak, just below the reflectivity peak, with a value of about 0.42 dB per degree. Its value in rain below the ML is 0.33 dB per degree, which is in rather good agreement with the slope of the linear k–Kdp relationship derived from DSD measurements at ground level. The PIA∕ϕdp ratio remains quite high in the upper part of the ML, between 0.32 and 0.38 dB per degree, before tending towards 0 above the ML.

Using ground radar overlaps to verify the retrieval of calibration bias estimates from spaceborne platforms

Abstract. Many institutions struggle to tap into the potential of their large archives of radar reflectivity: these data are often affected by miscalibration, yet the bias is typically unknown and temporally volatile. Still, relative calibration techniques can be used to correct the measurements a posteriori. For that purpose, the usage of spaceborne reflectivity observations from the Tropical Rainfall Measuring Mission (TRMM) and Global Precipitation Measurement (GPM) platforms has become increasingly popular: the calibration bias of a ground radar (GR) is estimated from its average reflectivity difference to the spaceborne radar (SR). Recently, Crisologo et al. (2018) introduced a formal procedure to enhance the reliability of such estimates: each match between SR and GR observations is assigned a quality index, and the calibration bias is inferred as a quality-weighted average of the differences between SR and GR. The relevance of quality was exemplified for the Subic S-band radar in the Philippines, which is greatly affected by partial beam blockage. The present study extends the concept of quality-weighted averaging by accounting for path-integrated attenuation (PIA) in addition to beam blockage. This extension becomes vital for radars that operate at the C or X band. Correspondingly, the study setup includes a C-band radar that substantially overlaps with the S-band radar. Based on the extended quality-weighting approach, we retrieve, for each of the two ground radars, a time series of calibration bias estimates from suitable SR overpasses. As a result of applying these estimates to correct the ground radar observations, the consistency between the ground radars in the region of overlap increased substantially. Furthermore, we investigated if the bias estimates can be interpolated in time, so that ground radar observations can be corrected even in the absence of prompt SR overpasses. We found that a moving average approach was most suitable for that purpose, although limited by the absence of explicit records of radar maintenance operations.

A Method for Assessing Relative Skill in Retrieving Cloud and Precipitation Properties in Next-Generation Cloud Radar and Radiometer Orbiting Observatories

AbstractA Bayesian Markov chain Monte Carlo (MCMC) algorithm is utilized to compare the skill of an A-Train-like observing system with a cloud, convection, and precipitation (CCP) observing system like that contemplated for the 2020s by the 2017 National Academy of Sciences Decadal Survey. The main objective is to demonstrate a framework for observational trade space studies. This initial work focuses on weakly precipitating warm shallow cumulus constructed from in situ data. Radiative computations are based on Mie theory with spherical assumptions. Simulated measurements in the CCP configuration consist of W- and Ka-band radar reflectivity and path-integrated attenuation, 31 and 94 GHz brightness temperatures (Tb), and visible and near-infrared reflectances. The collection of measurements in the CloudSat configuration is identical, but includes a single 94 GHz radar frequency, and the uncertainty in the 94 GHz microwave brightness temperature is increased to mimic the CloudSat Tb product. The experiments demonstrate that it remains a challenge to diagnose cloud properties in the presence of light rain because of the tendency of microwave remote sensing to respond to the higher moments of the hydrometeor populations. Rain properties are significantly better constrained than cloud properties, even in the optimal CCP configuration. The addition of Ka-band measurements places substantial constraints on the precipitation rain effective radius and rain rates. The Tb offers important information regarding the column-integrated condensate mass, the measurement accuracy of which appears more likely to affect the retrievals of clouds with low liquid water path. The constraints provided by reflectances are largely restricted to regions near the cloud top, particularly in the raining cases.

A Prototype Quantitative Precipitation Estimation Algorithm for Operational S-Band Polarimetric Radar Utilizing Specific Attenuation and Specific Differential Phase. Part I: Algorithm Description

Abstract A prototype quantitative precipitation estimate (QPE) algorithm that utilizes specific attenuation A and specific differential phase KDP was developed for inclusion into the Multi-Radar Multi-Sensor (MRMS) system and the Weather Surveillance Radar-1988 Doppler (WSR-88D) network. Special attention is given to the optimization of the factor α used for computation of a path-integrated attenuation from a total span of differential phase along the propagation path in rain. It is suggested to estimate α from a slope of the ZDR dependence on Z in rain. The use of real-time adjusted α allows us to capture the variations of the drop size distributions, and therefore improve the QPE accuracy. It is demonstrated that the factor α is generally higher for tropical rain type compared to continental rain. Since the R(A) approach is only valid for pure rainfall, the R(KDP) relation is suggested as a complement in areas of hail contamination. The paper contains a description of the basic version of the R(A) and R(KDP) algorithm and recommendations for its further optimization.

Relating GPM Radar Reflectivity Profile Characteristics to Path-Integrated Attenuation

The Global Precipitation Measurement (GPM) mission was launched in February 2014; its dual-frequency precipitation radar (DPR) operates at both Ku- and Ka-band s. Attenuation in precipitation is typically not negligible, especially at Ka-band. Hence, attenuation correction is an important part of the GPM DPR retrieval algorithm. The operational algorithm uses a path-integrated attenuation (PIA) obtained by comparing the measured surface return with that expected either from a nearby, nonprecipitating area or from the same area, acquired at a previous, nonprecipitating time. This surface reference technique has worked well in most situations but can result in erroneously low estimates of the path attenuation in situations with nonuniform filling of the radar beam, especially at Ka-band due to its larger attenuation. This paper explores the existence of relationships between the Ka-band PIA and the characteristics of the measured reflectivity profiles. The author finds that PIA is, indeed, related to reflectivity profiles, with strongest correlation between the PIA and the measured rainfall dual-frequency reflectivity ratio just above the surface. This relationship could be used as an estimator in the cases with severe nonuniform beam filling.

Relationships between Extreme Rain Rates and Convective Intensities from the Perspectives of TRMM and WSR-88D Radars

AbstractTRMM PR 2A25, version 7 (V7), retrievals of reflectivity Z and rainfall rate R are compared with WSR-88D dual-polarimetric S-band radar data for 28 radars over the southeastern United States after matching their horizontal resolution and sampling. TRMM Ku-band measurements are converted to S-band approximations to more directly compare reflectivity estimates. Rain rates are approximated from WSR-88D data using the CSU–hydrometeor identification rainfall optimization (HIDRO) algorithm. Tropics-wide TRMM retrievals confirm previous findings of a low overlap fraction between extreme convective intensity, as approximated by the maximum 40-dBZ height, and extreme near-surface rain rates. WSR-88D data also confirm this low overlap but show that it is likely higher than TRMM PR retrievals indicate. For maximum 40-dBZ echo heights that extend above the freezing level, mean WSR-88D reflectivities at low levels are approximately 2 dB higher than TRMM PR reflectivities. Higher WSR-88D-retrieved rain rates for a given low-level reflectivity combine with these higher low-level reflectivities for a given maximum 40-dBZ height to produce rain rates that are approximately double those retrieved by the TRMM PR for maximum 40-dBZ heights that extend above the freezing level. TRMM PR path-integrated attenuation, and WSR-88D specific differential phase, differential reflectivity, and hail fraction indicate that the TRMM PR 2A25 V7 algorithm is possibly misidentifying low–midlevel hail and/or graupel as greater attenuating liquid, or vice versa. This misidentification, coupled with underestimation of path-integrated attenuation caused by nonuniform beamfilling and higher rain rates produced by specific differential phase (KDP)–R than Z–R relationships, results in low-biased 2A25 V7 rain rates in intense convection.