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
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 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
  3. 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
  4. 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 multibeam 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.

     
    more » « less
  5. Abstract

    Tropical easterly waves (TEWs) are a recurrent mode of low‐latitude weather that are often convectively coupled and impact precipitation extremes. Previous work has examined the development of TEWs and their associated precipitation for individual seasons or regional domains, but no studies exist that document the importance of TEW precipitation globally. This study quantifies the precipitation associated with TEWs across the entire tropics using satellite (Integrated Multi‐satellitE Retrievals for the Global Precipitation Measurement [IMERG]) and reanalysis (Modern‐Era Retrospective analysis for Research and Applications, Version 2 [MERRA‐2]) data. Traditional space‐time filtering of precipitation reveals a mostly similar climatological power distribution for westward traveling, synoptic period disturbances corresponding to TEWs within all data sets. Using objective tracking, we find that areas with maximum TEW frequency such as the North Atlantic, Equatorial Pacific, and Indian Ocean have the highest accumulation of TEW‐associated precipitation. TEWs account for at most 30% of total annual precipitation in regions where they commonly occur and 1%–5% over much of the tropics. Vertically collocated storms, where the 850 and 700 hPa tracks correspond with each other, have higher conditional rain rates and indicate that waves with vertical development produce stronger and more organized convection. We find similar regional patterns using MERRA‐2 precipitation and latent heating, although the importance and contribution of TEWs to the background are reduced compared to IMERG. While the broad pattern of TEW associated precipitation in MERRA‐2 is like observations, the underestimation of rainfall contributions from TEWs, coupled with occasional false alarms in reanalysis data, suggests that MERRA‐2 does not capture organized convection within TEWs correctly.

     
    more » « less