skip to main content


Title: How Well Do Multisatellite Products Capture the Space–Time Dynamics of Precipitation? Part II: Building an Error Model through Spectral System Identification
Abstract

Satellite precipitation products, as all quantitative estimates, come with some inherent degree of uncertainty. To associate a quantitative value of the uncertainty to each individual estimate, error modeling is necessary. Most of the error models proposed so far compute the uncertainty as a function of precipitation intensity only, and only at one specific spatiotemporal scale. We propose a spectral error model that accounts for the neighboring space–time dynamics of precipitation into the uncertainty quantification. Systematic distortions of the precipitation signal and random errors are characterized distinctively in every frequency–wavenumber band in the Fourier domain, to accurately characterize error across scales. The systematic distortions are represented as a deterministic space–time linear filtering term. The random errors are represented as a nonstationary additive noise. The spectral error model is applied to the IMERG multisatellite precipitation product, and its parameters are estimated empirically through a system identification approach using the GV-MRMS gauge–radar measurements as reference (“truth”) over the eastern United States. The filtering term is found to be essentially low-pass (attenuating the fine-scale variability). While traditional error models attribute most of the error variance to random errors, it is found here that the systematic filtering term explains 48% of the error variance at the native resolution of IMERG. This fact confirms that, at high resolution, filtering effects in satellite precipitation products cannot be ignored, and that the error cannot be represented as a purely random additive or multiplicative term. An important consequence is that precipitation estimates derived from different sources shall not be expected to automatically have statistically independent errors.

Significance Statement

Satellite precipitation products are nowadays widely used for climate and environmental research, water management, risk analysis, and decision support at the local, regional, and global scales. For all these applications, knowledge about the accuracy of the products is critical for their usability. However, products are not systematically provided with a quantitative measure of the uncertainty associated with each individual estimate. Various parametric error models have been proposed for uncertainty quantification, mostly assuming that the uncertainty is only a function of the precipitation intensity at the pixel and time of interest. By projecting satellite precipitation fields and their retrieval errors into the Fourier frequency–wavenumber domain, we show that we can explicitly take into account the neighboring space–time multiscale dynamics of precipitation and compute a scale-dependent uncertainty.

 
more » « less
Award ID(s):
1839441 1839336
NSF-PAR ID:
10370417
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
American Meteorological Society
Date Published:
Journal Name:
Journal of Hydrometeorology
Volume:
23
Issue:
9
ISSN:
1525-755X
Page Range / eLocation ID:
p. 1383-1399
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract Many existing models that predict landslide hazards utilize ground-based sources of precipitation data. In locations where ground-based precipitation observations are limited (i.e., a vast majority of the globe), or for landslide hazard models that assess regional or global domains, satellite multisensor precipitation products offer a promising near-real-time alternative to ground-based data. NASA’s global Landslide Hazard Assessment for Situational Awareness (LHASA) model uses the Integrated Multisatellite Retrievals for Global Precipitation Measurement (IMERG) product to issue hazard “nowcasts” in near–real time for areas that are currently at risk for landsliding. Satellite-based precipitation estimates, however, can contain considerable systematic bias and random error, especially over mountainous terrain and during extreme rainfall events. This study combines a precipitation error modeling framework with a probabilistic adaptation of LHASA. Compared with the routine version of LHASA, this probabilistic version correctly predicts more of the observed landslides in the study region with fewer false alarms by high hazard nowcasts. This study demonstrates that improvements in landslide hazard prediction can be achieved regardless of whether the IMERG error model is trained using abundant ground-based precipitation observations or using far fewer and more scattered observations, suggesting that the approach is viable in data-limited regions. Results emphasize the importance of accounting for both random error and systematic satellite precipitation bias. The approach provides an example of how environmental prediction models can incorporate satellite precipitation uncertainty. Other applications such as flood and drought monitoring and forecasting could likely benefit from consideration of precipitation uncertainty. 
    more » « less
  2. Abstract

    Understanding the nature and origin of errors in satellite precipitation products is important for applications and product improvement. Here we propose a new error decomposition scheme incorporating precipitation event (continuous rainy periods) information to characterize satellite errors. Under this framework, the errors are attributed to the inaccuracies in event occurrence, timing (event start/end time), and intensity. The Integrated MultisatellitE Retrieval for Global Precipitation Measurement (IMERG) is used as our test product to apply the method over CONUS. The above‐listed factors contribute approximately 30%, 20%, and 50% to the total bias, respectively. Significant asymmetry exists in the temporal distribution of biases throughout events: early event endings cause threefold more precipitation amount bias than late event beginnings, while early event beginnings cause fourfold more bias than late event endings. Dominant contributors vary across seasons and regions. The proposed error decomposition provides insight into sources of error for improved retrievals.

     
    more » « less
  3. Abstract Line intensity mapping (LIM) proposes to efficiently observe distant faint galaxies and map the matter density field at high redshift.Building upon the formalism in a companion paper,we first highlight the degeneracies between cosmology and astrophysics in LIM.We discuss what can be constrained from measurements of the mean intensity and redshift-space power spectra.With a sufficient spectral resolution, the large-scale redshift-space distortions of the 2-halo term can be measured, helping to break the degeneracy between bias and mean intensity.With a higher spectral resolution, measuring the small-scale redshift-space distortions disentangles the 1-halo and shot noise terms.Cross-correlations with external galaxy catalogs or lensing surveys further break degeneracies.We derive requirements for experiments similar to SPHEREx, HETDEX, CDIM, COMAP and CONCERTO.We then revisit the question of the optimality of the LIM observables, compared to galaxy detection, for astrophysics and cosmology.We use a matched filter to compute the luminosity detection threshold for individual sources.We show that LIM contains information about galaxies too faint to detect, in the high-noise or high-confusion regimes.We quantify the sparsity and clustering bias of the detected sources and compare them to LIM, showing in which cases LIM is a better tracer of the matter density.We extend previous work by answering these questions as a function of Fourier scale, including for the first time the effect of cosmic variance, pixel-to-pixel correlations, luminosity-dependent clustering bias and redshift-space distortions. 
    more » « less
  4. Abstract As more global satellite-derived precipitation products become available, it is imperative to evaluate them more carefully for providing guidance as to how well precipitation space-time features are captured for use in hydrologic modeling, climate studies and other applications. Here we propose a space-time Fourier spectral analysis and define a suite of metrics which evaluate the spatial organization of storm systems, the propagation speed and direction of precipitation features, and the space-time scales at which a satellite product reproduces the variability of a reference “ground-truth” product (“effective resolution”). We demonstrate how the methodology relates to our physical intuition using the case study of a storm system with rich space-time structure. We then evaluate five high-resolution multi-satellite products (CMORPH, GSMaP, IMERG-early, IMERG-final and PERSIANN-CCS) over a period of two years over the southeastern US. All five satellite products show generally consistent space-time power spectral density when compared to a reference ground gauge-radar dataset (GV-MRMS), revealing agreement in terms of average morphology and dynamics of precipitation systems. However, a deficit of spectral power at wavelengths shorter than 200 km and periods shorter than 4 h reveals that all satellite products are excessively “smooth”. The products also show low levels of spectral coherence with the gauge-radar reference at these fine scales, revealing discrepancies in capturing the location and timing of precipitation features. From the space-time spectral coherence, the IMERG-final product shows superior ability in resolving the space-time dynamics of precipitation down to 200 km and 4 h scales compared to the other products. 
    more » « less
  5. null (Ed.)
    ABSTRACT The launch of NOAA’s latest generation of geostationary satellites known as the Geostationary Operational Environmental Satellite (GOES)-R Series has opened new opportunities in quantifying precipitation rates. Recent efforts have strived to utilize these data to improve space-based precipitation retrievals. The overall objective of the present work is to carry out a detailed error budget analysis of the improved Self-Calibrating Multivariate Precipitation Retrieval (SCaMPR) algorithm for GOES-R and the passive microwave (MW) combined (MWCOMB) precipitation dataset used to calibrate it with an aim to provide insights regarding strengths and weaknesses of these products. This study systematically analyzes the errors across different climate regions and also as a function of different precipitation types over the conterminous United States. The reference precipitation dataset is Ground-Validation Multi-Radar Multi-Sensor (GV-MRMS). Overall, MWCOMB reveals smaller errors as compared to SCaMPR. However, the analysis indicated that that the major portion of error in SCaMPR is propagated from the MWCOMB calibration data. The major challenge starts with poor detection from MWCOMB, which propagates in SCaMPR. In particular, MWCOMB misses 90% of cool stratiform precipitation and the overall detection score is around 40%. The ability of the algorithms to quantify precipitation amounts for the Warm Stratiform, Cool Stratiform, and Tropical/Stratiform Mix categories is poor compared to the Convective and Tropical/Convective Mix categories with additional challenges in complex terrain regions. Further analysis showed strong similarities in systematic and random error models with both products. This suggests that the potential of high-resolution GOES-R observations remains underutilized in SCaMPR due to the errors from the calibrator MWCOMB. 
    more » « less