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

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 more » 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.

« less
 ;  ;  ;  ;  
Award ID(s):
1839441 1839336
Publication Date:
Journal Name:
Journal of Hydrometeorology
Page Range or eLocation-ID:
p. 1383-1399
American Meteorological Society
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 importancemore »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.« less
  2. 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 questionsmore »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.« less
  3. 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 coherencemore »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.« less
  4. 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,more »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.« less
  5. Abstract

    Models of bathymetry derived from satellite radar altimetry are essential for modeling many marine processes. They are affected by uncertainties which require quantification. We propose an uncertainty model that assumes errors are caused by the lack of high‐wavenumber content within the altimetry data. The model is then applied to a tsunami hazard assessment. We build a bathymetry uncertainty model for northern Chile. Statistical properties of the altimetry‐predicted bathymetry error are obtained using multibeam data. We find that a Von Karman correlation function and a Laplacian marginal distribution can be used to define an uncertainty model based on a random field. We also propose a method for generating synthetic bathymetry samples conditional to shipboard measurements. The method is further extended to account for interpolation uncertainties, when bathymetry data resolution is finer than10 km. We illustrate the usefulness of the method by quantifying the bathymetry‐induced uncertainty of a tsunami hazard estimate. We demonstrate that tsunami leading wave predictions at middle/near field tide gauges and buoys are insensitive to bathymetry uncertainties in Chile. This result implies that tsunami early warning approaches can take full advantage of altimetry‐predicted bathymetry in numerical simulations. Finally, we evaluate the feasibility of modeling uncertainties in regions without multibeammore »data by assessing the bathymetry error statistics of 15 globally distributed regions. We find that a general Von Karman correlation and a Laplacian marginal distribution can serve as a first‐order approximation. The standard deviation of the uncertainty random field model varies regionally and is estimated from a proposed scaling law.

    « less