Browse
Abstract
A key strategy in obtaining complete global coverage of high-resolution precipitation is to combine observations from multiple fields, such as the intermittent passive microwave observations, precipitation propagated in time using motion vectors, and geosynchronous infrared observations. These separate precipitation fields can be combined through weighted averaging, which produces estimates that are generally superior to the individual parent fields. However, the process of averaging changes the distribution of the precipitation values, leading to an increase in precipitating area and a decrease in the values of high precipitation rates, a phenomenon observed in IMERG. To mitigate this issue, we introduce a new scheme called SHARPEN (Scheme for Histogram Adjustment with Ranked Precipitation Estimates in the Neighborhood), which recovers the distribution of the averaged precipitation field based on the idea of quantile mapping applied to the local environment. When implemented in IMERG, precipitation estimates from SHARPEN exhibit a distribution that resembles that of the original instantaneous observations, with matching precipitating area and peak precipitation rates. Case studies demonstrate its improved ability in bridging between the parent precipitation fields. Evaluation against ground observations reveals a distinct improvement in precipitation detection skill, but also a slightly reduced correlation likely because of a sharper precipitation field. The increased computational demand of SHARPEN can be mitigated by striding over multiple grid boxes, which has only marginal impacts on the accuracy of the estimates. SHARPEN can be applied to any precipitation algorithm that produces an average from multiple input precipitation fields and is being considered for implementation in IMERG V07.
Abstract
A key strategy in obtaining complete global coverage of high-resolution precipitation is to combine observations from multiple fields, such as the intermittent passive microwave observations, precipitation propagated in time using motion vectors, and geosynchronous infrared observations. These separate precipitation fields can be combined through weighted averaging, which produces estimates that are generally superior to the individual parent fields. However, the process of averaging changes the distribution of the precipitation values, leading to an increase in precipitating area and a decrease in the values of high precipitation rates, a phenomenon observed in IMERG. To mitigate this issue, we introduce a new scheme called SHARPEN (Scheme for Histogram Adjustment with Ranked Precipitation Estimates in the Neighborhood), which recovers the distribution of the averaged precipitation field based on the idea of quantile mapping applied to the local environment. When implemented in IMERG, precipitation estimates from SHARPEN exhibit a distribution that resembles that of the original instantaneous observations, with matching precipitating area and peak precipitation rates. Case studies demonstrate its improved ability in bridging between the parent precipitation fields. Evaluation against ground observations reveals a distinct improvement in precipitation detection skill, but also a slightly reduced correlation likely because of a sharper precipitation field. The increased computational demand of SHARPEN can be mitigated by striding over multiple grid boxes, which has only marginal impacts on the accuracy of the estimates. SHARPEN can be applied to any precipitation algorithm that produces an average from multiple input precipitation fields and is being considered for implementation in IMERG V07.
Abstract
Applications like drought monitoring and forecasting can profit from the global and near-real-time availability of satellite-based precipitation estimates once their related uncertainties and challenges are identified and treated. To this end, this study evaluates the IMERG V06B Late Run precipitation product from the Global Precipitation Measurement mission (GPM), a multisatellite product that combines space-based radar, passive microwave (PMW), and infrared (IR) data into gridded precipitation estimates. The evaluation is performed on the spatiotemporal resolution of IMERG (0.1° × 0.1°, 30 min) over the Netherlands over a 5-yr period. A gauge-adjusted radar precipitation product from the Royal Netherlands Meteorological Institute (KNMI) is used as reference, against which IMERG shows a large positive bias. To find the origin of this systematic overestimation, the data are divided into seasons, rainfall intensity ranges, echo top height (ETH) ranges, and categories based on the relative contributions of IR, morphing, and PMW data to the IMERG estimates. Furthermore, the specific radiometer is identified for each PMW-based estimate. IMERG’s detection performance improves with higher ETH and rainfall intensity, but the associated error and relative bias increase as well. Severe overestimation occurs during low-intensity rainfall events and is especially linked to PMW observations. All individual PMW instruments show the same pattern: overestimation of low-intensity events and underestimation of high-intensity events. IMERG misses a large fraction of shallow rainfall events, which is amplified when IR data are included. Space-based retrieval of shallow and low-intensity precipitation events should improve before IMERG can become accurate over the middle and high latitudes.
Abstract
Applications like drought monitoring and forecasting can profit from the global and near-real-time availability of satellite-based precipitation estimates once their related uncertainties and challenges are identified and treated. To this end, this study evaluates the IMERG V06B Late Run precipitation product from the Global Precipitation Measurement mission (GPM), a multisatellite product that combines space-based radar, passive microwave (PMW), and infrared (IR) data into gridded precipitation estimates. The evaluation is performed on the spatiotemporal resolution of IMERG (0.1° × 0.1°, 30 min) over the Netherlands over a 5-yr period. A gauge-adjusted radar precipitation product from the Royal Netherlands Meteorological Institute (KNMI) is used as reference, against which IMERG shows a large positive bias. To find the origin of this systematic overestimation, the data are divided into seasons, rainfall intensity ranges, echo top height (ETH) ranges, and categories based on the relative contributions of IR, morphing, and PMW data to the IMERG estimates. Furthermore, the specific radiometer is identified for each PMW-based estimate. IMERG’s detection performance improves with higher ETH and rainfall intensity, but the associated error and relative bias increase as well. Severe overestimation occurs during low-intensity rainfall events and is especially linked to PMW observations. All individual PMW instruments show the same pattern: overestimation of low-intensity events and underestimation of high-intensity events. IMERG misses a large fraction of shallow rainfall events, which is amplified when IR data are included. Space-based retrieval of shallow and low-intensity precipitation events should improve before IMERG can become accurate over the middle and high latitudes.
Abstract
This paper describes a new Passive Microwave Empirical Cold Surface Classification Algorithm (PESCA) developed for snow-cover detection and characterization by using passive microwave satellite measurements. The main goal of PESCA is to support the retrieval of falling snow, since several studies have highlighted the influence of snow-cover radiative properties on the falling-snow passive microwave signature. The developed method is based on the exploitation of the lower-frequency channels (<90 GHz), common to most microwave radiometers. The method applied to the conically scanning Global Precipitation Measurement (GPM) Microwave Imager (GMI) and the cross-track-scanning Advanced Technology Microwave Sounder (ATMS) is described in this paper. PESCA is based on a decision tree developed using an empirical method and verified using the AutoSnow product built from satellite measurements. The algorithm performance appears to be robust both for sensors in dry conditions (total precipitable water < 10 mm) and for mean surface elevation < 2500 m, independent of the cloud cover. The algorithm shows very good performance for cold temperatures (2-m temperature below 270 K) with a rapid decrease of the detection capabilities between 270 and 280 K, where 280 K is assumed as the maximum temperature limit for PESCA (overall detection statistics: probability of detection is 0.98 for ATMS and 0.92 for GMI, false alarm ratio is 0.01 for ATMS and 0.08 for GMI, and Heidke skill score is 0.72 for ATMS and 0.69 for GMI). Some inconsistencies found between the snow categories identified with the two radiometers are related to their different viewing geometries, spatial resolution, and temporal sampling. The spectral signatures of the different snow classes also appear to be different at high frequency (>90 GHz), indicating potential impact for snowfall retrieval. This method can be applied to other conically scanning and cross-track-scanning radiometers, including the future operational EUMETSAT Polar System Second Generation (EPS-SG) mission microwave radiometers.
Abstract
This paper describes a new Passive Microwave Empirical Cold Surface Classification Algorithm (PESCA) developed for snow-cover detection and characterization by using passive microwave satellite measurements. The main goal of PESCA is to support the retrieval of falling snow, since several studies have highlighted the influence of snow-cover radiative properties on the falling-snow passive microwave signature. The developed method is based on the exploitation of the lower-frequency channels (<90 GHz), common to most microwave radiometers. The method applied to the conically scanning Global Precipitation Measurement (GPM) Microwave Imager (GMI) and the cross-track-scanning Advanced Technology Microwave Sounder (ATMS) is described in this paper. PESCA is based on a decision tree developed using an empirical method and verified using the AutoSnow product built from satellite measurements. The algorithm performance appears to be robust both for sensors in dry conditions (total precipitable water < 10 mm) and for mean surface elevation < 2500 m, independent of the cloud cover. The algorithm shows very good performance for cold temperatures (2-m temperature below 270 K) with a rapid decrease of the detection capabilities between 270 and 280 K, where 280 K is assumed as the maximum temperature limit for PESCA (overall detection statistics: probability of detection is 0.98 for ATMS and 0.92 for GMI, false alarm ratio is 0.01 for ATMS and 0.08 for GMI, and Heidke skill score is 0.72 for ATMS and 0.69 for GMI). Some inconsistencies found between the snow categories identified with the two radiometers are related to their different viewing geometries, spatial resolution, and temporal sampling. The spectral signatures of the different snow classes also appear to be different at high frequency (>90 GHz), indicating potential impact for snowfall retrieval. This method can be applied to other conically scanning and cross-track-scanning radiometers, including the future operational EUMETSAT Polar System Second Generation (EPS-SG) mission microwave radiometers.
Abstract
Surface snowfall rate estimates from the Global Precipitation Measurement (GPM) mission’s Core Observatory sensors and the CloudSat radar are compared to those from the Multi-Radar Multi-Sensor (MRMS) radar composite product over the continental United States during the period from November 2014 to September 2020. The analysis includes the Dual-Frequency Precipitation Radar (DPR) retrieval and its single-frequency counterparts, the GPM Combined Radar Radiometer Algorithm (CORRA), the CloudSat Snow Profile product (2C-SNOW-PROFILE), and two passive microwave retrievals, i.e., the Goddard Profiling algorithm (GPROF) and the Snow Retrieval Algorithm for GMI (SLALOM). The 2C-SNOW retrieval has the highest Heidke skill score (HSS) for detecting snowfall among the products analyzed. SLALOM ranks second; it outperforms GPROF and the other GPM algorithms, all detecting only 30% of the snow events. Since SLALOM is trained with 2C-SNOW, it suggests that the optimal use of the information content in the GMI observations critically depends on the precipitation training dataset. All the retrievals underestimate snowfall rates by a factor of 2 compared to MRMS. Large discrepancies (RMSE of 0.7–1.5 mm h−1) between spaceborne and ground-based snowfall rate estimates are attributed to the complexity of the ice scattering properties and to the limitations of the remote sensing systems: the DPR instrument has low sensitivity, while the radiometric measurements are affected by the confounding effects of the background surface emissivity and of the emission of supercooled liquid droplet layers.
Abstract
Surface snowfall rate estimates from the Global Precipitation Measurement (GPM) mission’s Core Observatory sensors and the CloudSat radar are compared to those from the Multi-Radar Multi-Sensor (MRMS) radar composite product over the continental United States during the period from November 2014 to September 2020. The analysis includes the Dual-Frequency Precipitation Radar (DPR) retrieval and its single-frequency counterparts, the GPM Combined Radar Radiometer Algorithm (CORRA), the CloudSat Snow Profile product (2C-SNOW-PROFILE), and two passive microwave retrievals, i.e., the Goddard Profiling algorithm (GPROF) and the Snow Retrieval Algorithm for GMI (SLALOM). The 2C-SNOW retrieval has the highest Heidke skill score (HSS) for detecting snowfall among the products analyzed. SLALOM ranks second; it outperforms GPROF and the other GPM algorithms, all detecting only 30% of the snow events. Since SLALOM is trained with 2C-SNOW, it suggests that the optimal use of the information content in the GMI observations critically depends on the precipitation training dataset. All the retrievals underestimate snowfall rates by a factor of 2 compared to MRMS. Large discrepancies (RMSE of 0.7–1.5 mm h−1) between spaceborne and ground-based snowfall rate estimates are attributed to the complexity of the ice scattering properties and to the limitations of the remote sensing systems: the DPR instrument has low sensitivity, while the radiometric measurements are affected by the confounding effects of the background surface emissivity and of the emission of supercooled liquid droplet layers.
Abstract
Rainfall retrieval algorithms for passive microwave radiometers often exploit the brightness temperature depression due to ice scattering at high-frequency channels (≥85 GHz) over land. This study presents an alternate method to estimate the daily rainfall amount using the emissivity temporal variation (i.e., Δe) under rain-free conditions at low-frequency channels (19, 24, and 37 GHz). Emissivity is derived from 10 passive microwave radiometers, including the Global Precipitation Measurement (GPM) Microwave Imager (GMI), the Advanced Microwave Scanning Radiometer 2 (AMSR2), three Special Sensor Microwave Imager/Sounders (SSMIS), the Advanced Technology Microwave Sounder (ATMS), and four Advanced Microwave Sounding Units-A (AMSU-A). Four different satellite combination schemes are used to derive the Δe for daily rainfall estimates. They are all 10 satellites, 5 imagers, 6 satellites with very different equator crossing times, and GMI only. Results show that Δe from all 10 satellites has the best performance with a correlation of 0.60 and RMSE of 6.52 mm, compared with the Integrated Multisatellite Retrievals for GPM (IMERG) Final run product. The 6-satellites scheme has comparable performance with the all-10-satellites scheme. The 5-imagers scheme performs noticeably worse with a correlation of 0.49 and RMSE of 7.28 mm, while the GMI-only scheme performs the worst with a correlation of 0.25 and RMSE of 11.36 mm. The inferior performance from the 5-imagers and GMI-only schemes can be explained by the much longer revisit time, which cannot accurately capture the emissivity temporal variation.
Abstract
Rainfall retrieval algorithms for passive microwave radiometers often exploit the brightness temperature depression due to ice scattering at high-frequency channels (≥85 GHz) over land. This study presents an alternate method to estimate the daily rainfall amount using the emissivity temporal variation (i.e., Δe) under rain-free conditions at low-frequency channels (19, 24, and 37 GHz). Emissivity is derived from 10 passive microwave radiometers, including the Global Precipitation Measurement (GPM) Microwave Imager (GMI), the Advanced Microwave Scanning Radiometer 2 (AMSR2), three Special Sensor Microwave Imager/Sounders (SSMIS), the Advanced Technology Microwave Sounder (ATMS), and four Advanced Microwave Sounding Units-A (AMSU-A). Four different satellite combination schemes are used to derive the Δe for daily rainfall estimates. They are all 10 satellites, 5 imagers, 6 satellites with very different equator crossing times, and GMI only. Results show that Δe from all 10 satellites has the best performance with a correlation of 0.60 and RMSE of 6.52 mm, compared with the Integrated Multisatellite Retrievals for GPM (IMERG) Final run product. The 6-satellites scheme has comparable performance with the all-10-satellites scheme. The 5-imagers scheme performs noticeably worse with a correlation of 0.49 and RMSE of 7.28 mm, while the GMI-only scheme performs the worst with a correlation of 0.25 and RMSE of 11.36 mm. The inferior performance from the 5-imagers and GMI-only schemes can be explained by the much longer revisit time, which cannot accurately capture the emissivity temporal variation.
Abstract
High-resolution datasets offer the potential to improve our understanding of spatial and temporal precipitation patterns and storm structures. The goal of this study is to evaluate the similarities and differences of object-based storm characteristics as observed using space- or land-based sensors. The Method of Object-based Diagnostic Evaluation (MODE) Time Domain (MTD) is used to identify and track storm objects in two high-resolution merged datasets: the Integrated Multisatellite Retrievals for Global Precipitation Measurement (IMERG) final product V06B and gauge-corrected ground-radar-based Multi-Radar Multi-Sensor (MRMS) quantitative precipitation estimations. Characteristics associated with landfalling hurricanes were also examined as a separate category of storm. The results reveal that IMERG and MRMS agree reasonably well across many object-based storm characteristics. However, there are some discrepancies that are statistically significant. MRMS storms are more concentrated, with smaller areas and higher peak intensities, which implies higher flash flood risks associated with the storms. On the other hand, IMERG storms can travel longer distances with a higher volume of precipitation, which implies higher risk of riverine flooding. Agreement between the datasets is higher for faster-moving hurricanes in terms of the averaged intensity. Finally, MRMS indicates a higher average precipitation intensity during the hurricane’s lifetime. However, in non-hurricanes, the opposite result was observed. This is likely related to MRMS having higher resolution; monitoring the hurricanes from many viewing angles, leading to different signal saturation properties compared to IMERG; and/or the dominance of droplet aggregation effects over evaporation effects at lower altitudes.
Abstract
High-resolution datasets offer the potential to improve our understanding of spatial and temporal precipitation patterns and storm structures. The goal of this study is to evaluate the similarities and differences of object-based storm characteristics as observed using space- or land-based sensors. The Method of Object-based Diagnostic Evaluation (MODE) Time Domain (MTD) is used to identify and track storm objects in two high-resolution merged datasets: the Integrated Multisatellite Retrievals for Global Precipitation Measurement (IMERG) final product V06B and gauge-corrected ground-radar-based Multi-Radar Multi-Sensor (MRMS) quantitative precipitation estimations. Characteristics associated with landfalling hurricanes were also examined as a separate category of storm. The results reveal that IMERG and MRMS agree reasonably well across many object-based storm characteristics. However, there are some discrepancies that are statistically significant. MRMS storms are more concentrated, with smaller areas and higher peak intensities, which implies higher flash flood risks associated with the storms. On the other hand, IMERG storms can travel longer distances with a higher volume of precipitation, which implies higher risk of riverine flooding. Agreement between the datasets is higher for faster-moving hurricanes in terms of the averaged intensity. Finally, MRMS indicates a higher average precipitation intensity during the hurricane’s lifetime. However, in non-hurricanes, the opposite result was observed. This is likely related to MRMS having higher resolution; monitoring the hurricanes from many viewing angles, leading to different signal saturation properties compared to IMERG; and/or the dominance of droplet aggregation effects over evaporation effects at lower altitudes.
Abstract
This study evaluates a machine learning–based precipitation ensemble technique (MLPET) over three mountainous tropical regions. The technique, based on quantile regression forests, integrates global satellite precipitation datasets from CMORPH, PERSIANN, GSMaP (V6), and 3B42 (V7) and an atmospheric reanalysis precipitation product (EI_GPCC) with daily soil moisture, specific humidity, and terrain elevation datasets. The complex terrain study areas include the Peruvian and Colombian Andes in South America and the Blue Nile in East Africa. Evaluation is performed at a daily time scale and 0.25° spatial resolution based on 13 years (2000–12) of reference rainfall data derived from dense in situ rain gauge networks. The technique is evaluated using K-fold, separately in each region, and leave-one-region-out validation experiments. Comparison of MLPET with the individual satellite and reanalysis precipitation datasets used for the blending and the recent Multi-Source Weighted-Ensemble Precipitation (MSWEP) global precipitation product exhibited improved systematic and random error statistics for all regions. In addition, it is shown that observations are encapsulated well within the ensemble envelope generated by the blending technique.
Abstract
This study evaluates a machine learning–based precipitation ensemble technique (MLPET) over three mountainous tropical regions. The technique, based on quantile regression forests, integrates global satellite precipitation datasets from CMORPH, PERSIANN, GSMaP (V6), and 3B42 (V7) and an atmospheric reanalysis precipitation product (EI_GPCC) with daily soil moisture, specific humidity, and terrain elevation datasets. The complex terrain study areas include the Peruvian and Colombian Andes in South America and the Blue Nile in East Africa. Evaluation is performed at a daily time scale and 0.25° spatial resolution based on 13 years (2000–12) of reference rainfall data derived from dense in situ rain gauge networks. The technique is evaluated using K-fold, separately in each region, and leave-one-region-out validation experiments. Comparison of MLPET with the individual satellite and reanalysis precipitation datasets used for the blending and the recent Multi-Source Weighted-Ensemble Precipitation (MSWEP) global precipitation product exhibited improved systematic and random error statistics for all regions. In addition, it is shown that observations are encapsulated well within the ensemble envelope generated by the blending technique.
Abstract
The study evaluated a numerical weather model (WRF)-based satellite precipitation adjustment technique with 81 heavy precipitation events that occurred in three tropical mountainous regions (Colombia, Peru, and Taiwan). The technique was applied on two widely used near-real-time global satellite precipitation products—the National Oceanic and Atmospheric Administration (NOAA) Climate Prediction Center morphing technique (CMORPH) and the Global Satellite Mapping of Precipitation project (GSMaP)—for each precipitation event. The WRF-adjusted satellite products along with the near-real-time and gauge-adjusted satellite products as well as the WRF simulation were evaluated by independent gauge networks at daily scale and event total scale. Results show that the near-real-time precipitation products exhibited severe underestimation relative to the gauge observations over the three tropical mountainous regions. The underestimation tended to be larger for higher rainfall accumulations. The WRF-based satellite adjustment provided considerable improvements to the near-real-time CMORPH and GSMaP products. Moreover, error metrics show that WRF-adjusted satellite products outperformed the gauge-adjusted counterparts for most of the events. The effectiveness of WRF-based satellite adjustment varied with events of different physical processes. Thus, the technique applied on satellite precipitation estimates of these events may exhibit inconsistencies in the bias correction.
Abstract
The study evaluated a numerical weather model (WRF)-based satellite precipitation adjustment technique with 81 heavy precipitation events that occurred in three tropical mountainous regions (Colombia, Peru, and Taiwan). The technique was applied on two widely used near-real-time global satellite precipitation products—the National Oceanic and Atmospheric Administration (NOAA) Climate Prediction Center morphing technique (CMORPH) and the Global Satellite Mapping of Precipitation project (GSMaP)—for each precipitation event. The WRF-adjusted satellite products along with the near-real-time and gauge-adjusted satellite products as well as the WRF simulation were evaluated by independent gauge networks at daily scale and event total scale. Results show that the near-real-time precipitation products exhibited severe underestimation relative to the gauge observations over the three tropical mountainous regions. The underestimation tended to be larger for higher rainfall accumulations. The WRF-based satellite adjustment provided considerable improvements to the near-real-time CMORPH and GSMaP products. Moreover, error metrics show that WRF-adjusted satellite products outperformed the gauge-adjusted counterparts for most of the events. The effectiveness of WRF-based satellite adjustment varied with events of different physical processes. Thus, the technique applied on satellite precipitation estimates of these events may exhibit inconsistencies in the bias correction.
Abstract
Monitoring changes of precipitation phase from space is important for understanding the mass balance of Earth’s cryosphere in a changing climate. This paper examines a Bayesian nearest neighbor approach for prognostic detection of precipitation and its phase using passive microwave observations from the Global Precipitation Measurement (GPM) satellite. The method uses the weighted Euclidean distance metric to search through an a priori database populated with coincident GPM radiometer and radar observations as well as ancillary snow-cover data. The algorithm performance is evaluated using data from GPM official precipitation products, ground-based radars, and high-fidelity simulations from the Weather Research and Forecasting Model. Using the presented approach, we demonstrate that the hit probability of terrestrial precipitation detection can reach to 0.80, while the probability of false alarm remains below 0.11. The algorithm demonstrates higher skill in detecting snowfall than rainfall, on average by 10%. In particular, the probability of precipitation detection and its solid phase increases by 11% and 8%, over dry snow cover, when compared to other surface types. The main reason is found to be related to the ability of the algorithm in capturing the signal of increased liquid water content in snowy clouds over radiometrically cold snow-covered surfaces.
Abstract
Monitoring changes of precipitation phase from space is important for understanding the mass balance of Earth’s cryosphere in a changing climate. This paper examines a Bayesian nearest neighbor approach for prognostic detection of precipitation and its phase using passive microwave observations from the Global Precipitation Measurement (GPM) satellite. The method uses the weighted Euclidean distance metric to search through an a priori database populated with coincident GPM radiometer and radar observations as well as ancillary snow-cover data. The algorithm performance is evaluated using data from GPM official precipitation products, ground-based radars, and high-fidelity simulations from the Weather Research and Forecasting Model. Using the presented approach, we demonstrate that the hit probability of terrestrial precipitation detection can reach to 0.80, while the probability of false alarm remains below 0.11. The algorithm demonstrates higher skill in detecting snowfall than rainfall, on average by 10%. In particular, the probability of precipitation detection and its solid phase increases by 11% and 8%, over dry snow cover, when compared to other surface types. The main reason is found to be related to the ability of the algorithm in capturing the signal of increased liquid water content in snowy clouds over radiometrically cold snow-covered surfaces.
Abstract
Hydrological applications rely on the availability and quality of precipitation products, especially model- and satellite-based products for use in areas without ground measurements. It is known that the quality of model- and satellite-based precipitation products is complementary: model-based products exhibit high quality during cold seasons while satellite-based products are better during warm seasons. To explore the complementary behavior of the quality of the precipitation products, this study uses 2-m air temperature as auxiliary information to evaluate high-resolution (0.1°/hourly) precipitation estimates from the Weather Research and Forecasting (WRF) Model and from the version 5 Integrated Multisatellite Retrievals for GPM (IMERG) algorithm (i.e., early and final runs). The products are evaluated relative to the reference NCEP Stage IV precipitation estimates over the central United States during August 2015–July 2017. Results show that the IMERG final-run estimates are nearly unbiased, while the IMERG early-run and the WRF estimates are positively biased. The WRF estimates exhibit high correlations with the reference data when the temperature falls below 280 K. The IMERG estimates, both early and final runs, do so when the temperature exceeds 280 K. Moreover, the complementary behavior of the WRF and the IMERG products conditioned on air temperature does not vary with either season or location.
Abstract
Hydrological applications rely on the availability and quality of precipitation products, especially model- and satellite-based products for use in areas without ground measurements. It is known that the quality of model- and satellite-based precipitation products is complementary: model-based products exhibit high quality during cold seasons while satellite-based products are better during warm seasons. To explore the complementary behavior of the quality of the precipitation products, this study uses 2-m air temperature as auxiliary information to evaluate high-resolution (0.1°/hourly) precipitation estimates from the Weather Research and Forecasting (WRF) Model and from the version 5 Integrated Multisatellite Retrievals for GPM (IMERG) algorithm (i.e., early and final runs). The products are evaluated relative to the reference NCEP Stage IV precipitation estimates over the central United States during August 2015–July 2017. Results show that the IMERG final-run estimates are nearly unbiased, while the IMERG early-run and the WRF estimates are positively biased. The WRF estimates exhibit high correlations with the reference data when the temperature falls below 280 K. The IMERG estimates, both early and final runs, do so when the temperature exceeds 280 K. Moreover, the complementary behavior of the WRF and the IMERG products conditioned on air temperature does not vary with either season or location.
