Browse
Abstract
Using passive microwave brightness temperatures Tb from the Global Precipitation Measurement (GPM) Microwave Imager (GMI) and hydrometeor identification (HID) data from dual-polarization ground radars, empirical lookup tables are developed for a multifrequency estimation of the likelihood a precipitation column includes certain hydrometeor types, as a function of Tb . Eight years of collocated Tb and HID data from the GPM Validation Network are used for development and testing of the GMI-based HID retrieval, with 2015–20 used for training and 2021–22 used for testing the GMI-based HID retrieval. The occurrence of profiles with hail and graupel are both slightly underpredicted by the lookup tables, but the percentage of profiles predicted is highly correlated with the percentage observed (0.98 correlation coefficient for hail and 0.99 for graupel). By having snow appear before rain in the hierarchy, the sample size for rain, without ice aloft, is fairly small, and the percentage of rain profiles is less than snow for all Tb .
Abstract
Using passive microwave brightness temperatures Tb from the Global Precipitation Measurement (GPM) Microwave Imager (GMI) and hydrometeor identification (HID) data from dual-polarization ground radars, empirical lookup tables are developed for a multifrequency estimation of the likelihood a precipitation column includes certain hydrometeor types, as a function of Tb . Eight years of collocated Tb and HID data from the GPM Validation Network are used for development and testing of the GMI-based HID retrieval, with 2015–20 used for training and 2021–22 used for testing the GMI-based HID retrieval. The occurrence of profiles with hail and graupel are both slightly underpredicted by the lookup tables, but the percentage of profiles predicted is highly correlated with the percentage observed (0.98 correlation coefficient for hail and 0.99 for graupel). By having snow appear before rain in the hierarchy, the sample size for rain, without ice aloft, is fairly small, and the percentage of rain profiles is less than snow for all Tb .
Abstract
The Dual-Frequency Precipitation Radar (DPR), which consists of a Ku-band precipitation radar (KuPR) and a Ka-band precipitation radar (KaPR) on board the GPM Core Observatory, cannot observe precipitation at low altitudes near the ground contaminated by surface clutter. This near-surface region is called the blind zone. DPR estimates the clutter-free bottom (CFB), which is the lowest altitude not included in the blind zone, and estimates precipitation at altitudes higher than the CFB. High CFBs, which are common over mountainous areas, represent obstacles to detection of shallow precipitation and estimation of low-level enhanced precipitation. We compared KuPR data with rain gauge data from Da-Tun Mountain of northern Taiwan acquired from March 2014 to February 2020. A total of 12 cases were identified in which the KuPR missed some rainfall with intensity of >10 mm h−1 that was observed by rain gauges. Comparison of KuPR profile and ground-based radar profile revealed that shallow precipitation in the KuPR blind zone was missed because the CFB was estimated to be higher than the lower bound of the range free from surface echoes. In the original operational algorithm, CFB was estimated using only the received power data of the KuPR. In this study, the CFB was identified by the sharp increase in the difference between the received powers of the KuPR and the KaPR at altitude affected by surface clutter. By lowering the CFB, the KuPR succeeded in detection and estimation of shallow precipitation.
Significance Statement
The Dual-Frequency Precipitation Radar (DPR) on board the GPM Core Observatory cannot capture precipitation in the low-altitude region near the ground contaminated by surface clutter. This region is called the blind zone. The DPR estimates the clutter-free bottom (CFB), which is the lower bound of the range free from surface echoes, and uses data higher than CFB. DPR consists of a Ku-band precipitation radar (KuPR) and a Ka-band precipitation radar (KaPR). KuPR missed some shallow precipitation more than 10 mm h−1 in the blind zone over Da-Tun Mountain of northern Taiwan because of misjudged CFB estimation. Using both the KuPR and the KaPR, we improved the CFB estimation algorithm, which lowered the CFB, narrowed the blind zone, and improved the capability to detect shallow precipitation.
Abstract
The Dual-Frequency Precipitation Radar (DPR), which consists of a Ku-band precipitation radar (KuPR) and a Ka-band precipitation radar (KaPR) on board the GPM Core Observatory, cannot observe precipitation at low altitudes near the ground contaminated by surface clutter. This near-surface region is called the blind zone. DPR estimates the clutter-free bottom (CFB), which is the lowest altitude not included in the blind zone, and estimates precipitation at altitudes higher than the CFB. High CFBs, which are common over mountainous areas, represent obstacles to detection of shallow precipitation and estimation of low-level enhanced precipitation. We compared KuPR data with rain gauge data from Da-Tun Mountain of northern Taiwan acquired from March 2014 to February 2020. A total of 12 cases were identified in which the KuPR missed some rainfall with intensity of >10 mm h−1 that was observed by rain gauges. Comparison of KuPR profile and ground-based radar profile revealed that shallow precipitation in the KuPR blind zone was missed because the CFB was estimated to be higher than the lower bound of the range free from surface echoes. In the original operational algorithm, CFB was estimated using only the received power data of the KuPR. In this study, the CFB was identified by the sharp increase in the difference between the received powers of the KuPR and the KaPR at altitude affected by surface clutter. By lowering the CFB, the KuPR succeeded in detection and estimation of shallow precipitation.
Significance Statement
The Dual-Frequency Precipitation Radar (DPR) on board the GPM Core Observatory cannot capture precipitation in the low-altitude region near the ground contaminated by surface clutter. This region is called the blind zone. The DPR estimates the clutter-free bottom (CFB), which is the lower bound of the range free from surface echoes, and uses data higher than CFB. DPR consists of a Ku-band precipitation radar (KuPR) and a Ka-band precipitation radar (KaPR). KuPR missed some shallow precipitation more than 10 mm h−1 in the blind zone over Da-Tun Mountain of northern Taiwan because of misjudged CFB estimation. Using both the KuPR and the KaPR, we improved the CFB estimation algorithm, which lowered the CFB, narrowed the blind zone, and improved the capability to detect shallow precipitation.
Abstract
Several decades of continuous improvements in satellite precipitation algorithms have resulted in fairly accurate level-2 precipitation products for local-scale applications. Numerous studies have been carried out to quantify random and systematic errors at individual validation sites and regional networks. Understanding uncertainties at larger scales, however, has remained a challenge. Temporal changes in precipitation regional biases, regime morphology, sampling, and observation-vector information content, all play important roles in defining the accuracy of satellite rainfall retrievals. This study considers these contributors to offer a quantitative estimate of uncertainty in recently produced global precipitation climate data record. Generated from intercalibrated observations collected by a constellation of passive microwave (PMW) radiometers over the course of 30 years, this data record relies on Global Precipitation Measurement (GPM) mission enterprise PMW precipitation retrieval to offer a long-term global monthly precipitation estimates with corresponding uncertainty at 5° scales. To address changes in the information content across different constellation members the study develops synthetic datasets from GPM Microwave Imager (GMI) sensor, while sampling- and morphology-related uncertainties are quantified using GPM’s dual-frequency precipitation radar (DPR). Special attention is given to separating precipitation into self-similar states that appear to be consistent across environmental conditions. Results show that the variability of bias patterns can be explained by the relative occurrence of different precipitation states across the regions and used to calculate product’s uncertainty. It is found that at 5° spatial scale monthly mean precipitation uncertainties in tropics can exceed 10%.
Abstract
Several decades of continuous improvements in satellite precipitation algorithms have resulted in fairly accurate level-2 precipitation products for local-scale applications. Numerous studies have been carried out to quantify random and systematic errors at individual validation sites and regional networks. Understanding uncertainties at larger scales, however, has remained a challenge. Temporal changes in precipitation regional biases, regime morphology, sampling, and observation-vector information content, all play important roles in defining the accuracy of satellite rainfall retrievals. This study considers these contributors to offer a quantitative estimate of uncertainty in recently produced global precipitation climate data record. Generated from intercalibrated observations collected by a constellation of passive microwave (PMW) radiometers over the course of 30 years, this data record relies on Global Precipitation Measurement (GPM) mission enterprise PMW precipitation retrieval to offer a long-term global monthly precipitation estimates with corresponding uncertainty at 5° scales. To address changes in the information content across different constellation members the study develops synthetic datasets from GPM Microwave Imager (GMI) sensor, while sampling- and morphology-related uncertainties are quantified using GPM’s dual-frequency precipitation radar (DPR). Special attention is given to separating precipitation into self-similar states that appear to be consistent across environmental conditions. Results show that the variability of bias patterns can be explained by the relative occurrence of different precipitation states across the regions and used to calculate product’s uncertainty. It is found that at 5° spatial scale monthly mean precipitation uncertainties in tropics can exceed 10%.
Abstract
The accuracy and uncertainty of radar echo-top heights estimated by ground-based radars remain largely unknown despite their critical importance for applications ranging from aviation weather forecasting to severe weather diagnosis. Because the vantage point of space is more suited than that of ground-based radars for the estimation of echo-top heights, the use of spaceborne radar observations is explored as an external reference for cross comparison. An investigation has been carried out across the conterminous United States by comparing the NOAA/National Severe Storms Laboratory Multi-Radar Multi-Sensor (MRMS) system with the space-based radar on board the NASA–JAXA Global Precipitation Measurement satellite platform. No major bias was assessed between the two products. An annual cycle of differences is found, driven by an underestimation of the stratiform cloud echo-top heights and an overestimation of the convective ones. The investigation of the systematic biases for different radar volume coverage patterns (VCP) shows that scanning strategies with fewer tilts and greater voids as VCP 21/121/221 contribute to overestimations observed for high MRMS tops. For VCP 12/212, the automated volume scan evaluation and termination (AVSET) function increases the radar cone of silence, causing overestimations when the echo top lies above the highest elevation scan. However, it seems that for low echo tops the shorter refresh rates contribute to mitigate underestimations, especially in stratiform cases.
Abstract
The accuracy and uncertainty of radar echo-top heights estimated by ground-based radars remain largely unknown despite their critical importance for applications ranging from aviation weather forecasting to severe weather diagnosis. Because the vantage point of space is more suited than that of ground-based radars for the estimation of echo-top heights, the use of spaceborne radar observations is explored as an external reference for cross comparison. An investigation has been carried out across the conterminous United States by comparing the NOAA/National Severe Storms Laboratory Multi-Radar Multi-Sensor (MRMS) system with the space-based radar on board the NASA–JAXA Global Precipitation Measurement satellite platform. No major bias was assessed between the two products. An annual cycle of differences is found, driven by an underestimation of the stratiform cloud echo-top heights and an overestimation of the convective ones. The investigation of the systematic biases for different radar volume coverage patterns (VCP) shows that scanning strategies with fewer tilts and greater voids as VCP 21/121/221 contribute to overestimations observed for high MRMS tops. For VCP 12/212, the automated volume scan evaluation and termination (AVSET) function increases the radar cone of silence, causing overestimations when the echo top lies above the highest elevation scan. However, it seems that for low echo tops the shorter refresh rates contribute to mitigate underestimations, especially in stratiform cases.
Abstract
Parameters of the normalized gamma particle size distribution (PSD) have been retrieved from the Precipitation Image Package (PIP) snowfall observations collected during the International Collaborative Experiment–PyeongChang Olympic and Paralympic winter games (ICE-POP 2018). Two of the gamma PSD parameters, the mass-weighted particle diameter D mass and the normalized intercept parameter NW , have median values of 1.15–1.31 mm and 2.84–3.04 log(mm−1 m−3), respectively. This range arises from the choice of the relationship between the maximum versus equivalent diameter, D mx–D eq, and the relationship between the Reynolds and Best numbers, Re–X. Normalization of snow water equivalent rate (SWER) and ice water content W by NW reduces the range in NW , resulting in well-fitted power-law relationships between SWER/NW and D mass and between W/NW and D mass. The bulk descriptors of snowfall are calculated from PIP observations and from the gamma PSD with values of the shape parameter μ ranging from −2 to 10. NASA’s Global Precipitation Measurement (GPM) mission, which adopted the normalized gamma PSD, assumes μ = 2 and 3 in its two separate algorithms. The mean fractional bias (MFB) of the snowfall parameters changes with μ, where the functional dependence on μ depends on the specific snowfall parameter of interest. The MFB of the total concentration was underestimated by 0.23–0.34 when μ = 2 and by 0.29–0.40 when μ = 3, whereas the MFB of SWER had a much narrower range (from −0.03 to 0.04) for the same μ values.
Abstract
Parameters of the normalized gamma particle size distribution (PSD) have been retrieved from the Precipitation Image Package (PIP) snowfall observations collected during the International Collaborative Experiment–PyeongChang Olympic and Paralympic winter games (ICE-POP 2018). Two of the gamma PSD parameters, the mass-weighted particle diameter D mass and the normalized intercept parameter NW , have median values of 1.15–1.31 mm and 2.84–3.04 log(mm−1 m−3), respectively. This range arises from the choice of the relationship between the maximum versus equivalent diameter, D mx–D eq, and the relationship between the Reynolds and Best numbers, Re–X. Normalization of snow water equivalent rate (SWER) and ice water content W by NW reduces the range in NW , resulting in well-fitted power-law relationships between SWER/NW and D mass and between W/NW and D mass. The bulk descriptors of snowfall are calculated from PIP observations and from the gamma PSD with values of the shape parameter μ ranging from −2 to 10. NASA’s Global Precipitation Measurement (GPM) mission, which adopted the normalized gamma PSD, assumes μ = 2 and 3 in its two separate algorithms. The mean fractional bias (MFB) of the snowfall parameters changes with μ, where the functional dependence on μ depends on the specific snowfall parameter of interest. The MFB of the total concentration was underestimated by 0.23–0.34 when μ = 2 and by 0.29–0.40 when μ = 3, whereas the MFB of SWER had a much narrower range (from −0.03 to 0.04) for the same μ values.
Abstract
We use dual-polarization C-band data collected in the Southern Ocean to examine the properties of snow observed during a voyage in the austral summer of 2018. Using existing forward modeling formalisms based on an assumption of Rayleigh scattering by soft spheroids, an optimal estimation algorithm is implemented to infer snow properties from horizontally polarized radar reflectivity, the differential radar reflectivity, and the specific differential phase. From the dual-polarization observables, we estimate ice water content qi, the mass-mean particle size Dm, and the exponent of the mass–dimensional relationship bm that, with several assumptions, allow for evaluation of snow bulk density, and snow number concentration. Upon evaluating the uncertainties associated with measurement and forward model errors, we determine that the algorithm can retrieve qi, Dm, and bm within single-pixel uncertainties conservatively estimated in the range 120%, 60%, and 40%, respectively. Applying the algorithm to open-cellular convection in the Southern Ocean, we find evidence for secondary ice formation processes within multicellular complexes. In stratiform precipitation systems we find snow properties and infer processes that are distinctly different from the shallow convective systems with evidence for riming and aggregation being common. We also find that embedded convection within the frontal system produces precipitation properties consistent with graupel. Examining 5 weeks of data, we show that snow in open-cellular cumulus has higher overall bulk density than snow in stratiform precipitation systems with implications for interpreting measurements from space-based active remote sensors.
Abstract
We use dual-polarization C-band data collected in the Southern Ocean to examine the properties of snow observed during a voyage in the austral summer of 2018. Using existing forward modeling formalisms based on an assumption of Rayleigh scattering by soft spheroids, an optimal estimation algorithm is implemented to infer snow properties from horizontally polarized radar reflectivity, the differential radar reflectivity, and the specific differential phase. From the dual-polarization observables, we estimate ice water content qi, the mass-mean particle size Dm, and the exponent of the mass–dimensional relationship bm that, with several assumptions, allow for evaluation of snow bulk density, and snow number concentration. Upon evaluating the uncertainties associated with measurement and forward model errors, we determine that the algorithm can retrieve qi, Dm, and bm within single-pixel uncertainties conservatively estimated in the range 120%, 60%, and 40%, respectively. Applying the algorithm to open-cellular convection in the Southern Ocean, we find evidence for secondary ice formation processes within multicellular complexes. In stratiform precipitation systems we find snow properties and infer processes that are distinctly different from the shallow convective systems with evidence for riming and aggregation being common. We also find that embedded convection within the frontal system produces precipitation properties consistent with graupel. Examining 5 weeks of data, we show that snow in open-cellular cumulus has higher overall bulk density than snow in stratiform precipitation systems with implications for interpreting measurements from space-based active remote sensors.
Abstract
With the launch of the Global Precipitation Measurement Dual-Frequency Precipitation Radar (GPM-DPR) in 2014, renewed interest in retrievals of snowfall in the atmospheric column has occurred. The current operational GPM-DPR retrieval largely underestimates surface snowfall accumulation. Here, a neural network (NN) trained on data that are synthetically derived from state-of-the-art ice particle scattering models and measured in situ particle size distributions (PSDs) is used to retrieve two parameters of the PSD: liquid equivalent mass-weighted mean diameter
Abstract
With the launch of the Global Precipitation Measurement Dual-Frequency Precipitation Radar (GPM-DPR) in 2014, renewed interest in retrievals of snowfall in the atmospheric column has occurred. The current operational GPM-DPR retrieval largely underestimates surface snowfall accumulation. Here, a neural network (NN) trained on data that are synthetically derived from state-of-the-art ice particle scattering models and measured in situ particle size distributions (PSDs) is used to retrieve two parameters of the PSD: liquid equivalent mass-weighted mean diameter
Abstract
Two power-law relations linking equivalent radar reflectivity factor Z e and snowfall rate S are derived for a K-band Micro Rain Radar (MRR) and for a W-band cloud radar. For the development of these Z e –S relationships, a dataset of calculated and measured variables is used. Surface-based video-disdrometer measurements were collected during snowfall events over five winters at the high-latitude site in Hyytiälä, Finland. The data from 2014 to 2018 include particle size distributions (PSD) and their fall velocities, from which snowflake masses were derived. The K- and W-band Z e values are computed using these surface-based observations and snowflake scattering properties as provided by T-matrix and single-particle scattering tables, respectively. The uncertainty analysis shows that the K-band snowfall-rate estimation is significantly improved by including the intercept parameter N 0 of the PSD calculated from concurrent disdrometer measurements. If N 0 is used to adjust the prefactor of the Z e –S relationship, the RMSE of the snowfall-rate estimate can be reduced from 0.37 to around 0.11 mm h−1. For W-band radar, a Z e –S relationship with constant parameters for all available snow events shows a similar uncertainty when compared with the method that includes the PSD intercept parameter. To demonstrate the performance of the proposed Z e –S relationships, they are applied to measurements of the MRR and the W-band microwave radar for Arctic clouds at the Arctic research base operated by the German Alfred Wegener Institute Helmholtz Centre for Polar and Marine Research (AWI) and the French Polar Institute Paul Emile Victor (IPEV) (AWIPEV) in Ny-Ålesund, Svalbard, Norway. The resulting snowfall-rate estimates show good agreement with in situ snowfall observations while other Z e –S relationships from literature reveal larger differences.
Abstract
Two power-law relations linking equivalent radar reflectivity factor Z e and snowfall rate S are derived for a K-band Micro Rain Radar (MRR) and for a W-band cloud radar. For the development of these Z e –S relationships, a dataset of calculated and measured variables is used. Surface-based video-disdrometer measurements were collected during snowfall events over five winters at the high-latitude site in Hyytiälä, Finland. The data from 2014 to 2018 include particle size distributions (PSD) and their fall velocities, from which snowflake masses were derived. The K- and W-band Z e values are computed using these surface-based observations and snowflake scattering properties as provided by T-matrix and single-particle scattering tables, respectively. The uncertainty analysis shows that the K-band snowfall-rate estimation is significantly improved by including the intercept parameter N 0 of the PSD calculated from concurrent disdrometer measurements. If N 0 is used to adjust the prefactor of the Z e –S relationship, the RMSE of the snowfall-rate estimate can be reduced from 0.37 to around 0.11 mm h−1. For W-band radar, a Z e –S relationship with constant parameters for all available snow events shows a similar uncertainty when compared with the method that includes the PSD intercept parameter. To demonstrate the performance of the proposed Z e –S relationships, they are applied to measurements of the MRR and the W-band microwave radar for Arctic clouds at the Arctic research base operated by the German Alfred Wegener Institute Helmholtz Centre for Polar and Marine Research (AWI) and the French Polar Institute Paul Emile Victor (IPEV) (AWIPEV) in Ny-Ålesund, Svalbard, Norway. The resulting snowfall-rate estimates show good agreement with in situ snowfall observations while other Z e –S relationships from literature reveal larger differences.
Abstract
Few statistical downscaling applications have provided gridded products that can provide downscaled values for a no-gauge area as is done by dynamical downscaling. In this study, a gridded statistical downscaling scheme is presented to downscale summer precipitation to a dense grid that covers North China. The main innovation of this scheme is interpolating the parameters of single-station models to this dense grid and assigning optimal predictor values according to an interpolated predictand–predictor distance function. This method can produce spatial dependence (spatial autocorrelation) and transmit the spatial heterogeneity of predictor values from the large-scale predictors to the downscaled outputs. Such gridded output at no-gauge stations shows performances comparable to that at the gauged stations. The area mean precipitation of the downscaled results is comparable to other products. The main value of the downscaling scheme is that it can obtain reasonable outputs for no-gauge stations.
Abstract
Few statistical downscaling applications have provided gridded products that can provide downscaled values for a no-gauge area as is done by dynamical downscaling. In this study, a gridded statistical downscaling scheme is presented to downscale summer precipitation to a dense grid that covers North China. The main innovation of this scheme is interpolating the parameters of single-station models to this dense grid and assigning optimal predictor values according to an interpolated predictand–predictor distance function. This method can produce spatial dependence (spatial autocorrelation) and transmit the spatial heterogeneity of predictor values from the large-scale predictors to the downscaled outputs. Such gridded output at no-gauge stations shows performances comparable to that at the gauged stations. The area mean precipitation of the downscaled results is comparable to other products. The main value of the downscaling scheme is that it can obtain reasonable outputs for no-gauge stations.
Abstract
Large hail is a primary contributor to damages and loss around the world, in both agriculture and infrastructure. The sensitivity of passive microwave radiometer measurements to scattering by hail led to the development of proxies for severe hail, most of which use brightness temperature thresholds from 37-GHz and higher-frequency microwave channels on board weather satellites in low-Earth orbit. Using 16+ years of data from the Tropical Rainfall Measuring Mission (TRMM; 36°S–36°N), we pair TRMM brightness temperature–derived precipitation features with surface hail reports in the United States to train a hail retrieval on passive microwave data from the 10-, 19-, 37-, and 85-GHz channels based on probability curves fit to the microwave data. We then apply this hail retrieval to features in the Global Precipitation Measurement (GPM) domain (from 69°S to 69°N) to develop a nearly global passive microwave–based climatology of hail. The extended domain of the GPM satellite into higher latitudes requires filtering out features that we believe are over icy and snowy surface regimes. We also normalize brightness temperature depression by tropopause height in an effort to account for differences in storm depth between the tropics and higher latitudes. Our results show the highest hail frequencies in the region of northern Argentina through Paraguay, Uruguay, and southern Brazil; the central United States; and a swath of Africa just south of the Sahel. Smaller hot spots include Pakistan, eastern India, and Bangladesh. A notable difference between these results and many prior satellite-based studies is that central Africa, while still active in our climatology, does not rival the aforementioned regions in retrieved hailstorm frequency.
Abstract
Large hail is a primary contributor to damages and loss around the world, in both agriculture and infrastructure. The sensitivity of passive microwave radiometer measurements to scattering by hail led to the development of proxies for severe hail, most of which use brightness temperature thresholds from 37-GHz and higher-frequency microwave channels on board weather satellites in low-Earth orbit. Using 16+ years of data from the Tropical Rainfall Measuring Mission (TRMM; 36°S–36°N), we pair TRMM brightness temperature–derived precipitation features with surface hail reports in the United States to train a hail retrieval on passive microwave data from the 10-, 19-, 37-, and 85-GHz channels based on probability curves fit to the microwave data. We then apply this hail retrieval to features in the Global Precipitation Measurement (GPM) domain (from 69°S to 69°N) to develop a nearly global passive microwave–based climatology of hail. The extended domain of the GPM satellite into higher latitudes requires filtering out features that we believe are over icy and snowy surface regimes. We also normalize brightness temperature depression by tropopause height in an effort to account for differences in storm depth between the tropics and higher latitudes. Our results show the highest hail frequencies in the region of northern Argentina through Paraguay, Uruguay, and southern Brazil; the central United States; and a swath of Africa just south of the Sahel. Smaller hot spots include Pakistan, eastern India, and Bangladesh. A notable difference between these results and many prior satellite-based studies is that central Africa, while still active in our climatology, does not rival the aforementioned regions in retrieved hailstorm frequency.
