Search for: All records

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. 

Abstract Measuring rainfall is complex, due to the high temporal and spatial variability of precipitation, especially in a changing climate, but it is of great importance for all the scientific and operational disciplines dealing with rainfall effects on the environment, human activities, and economy. Microwave (MW) telecommunication links carry information on rainfall rates along their path, through signal attenuation caused by raindrops, and can become measurements of opportunity, offering inexpensive chances to augment information without deploying additional infrastructures, at the cost of some smart processing. Processing satellite telecom signals bring some specific complexities related to the effects of rainfall boundaries, melting layer, and non-weather attenuations, but with the potential to provide worldwide precipitation data with high temporal and spatial samplings. These measurements have to be processed according to the probabilistic nature of the information they carry. An EnKF-based (Ensemble Kalman Filter) method has been developed to dynamically retrieve rainfall fields in gridded domains, which manages such probabilistic information and exploits the high sampling rate of measurements. The paper presents the EnKF method with some representative tests from synthetic 3D experiments. Ancillary data are assumed as from worldwide-available operational meteorological satellites and models, for advection, initial and boundary conditions, rain height. The method reproduces rainfall structures and quantities in a correct way, and also manages possible link outages. It results computationally viable also for operational implementation and applicable to different link observation geometries and characteristics. 

Abstract Accurate, physically based precipitation retrieval over global land surfaces is an important goal of the NASA/JAXA Global Precipitation Measurement Mission (GPM). This is a difficult problem for the passive microwave constellation, as the signal over radiometrically warm land surfaces in the microwave frequencies means that the measurements used are indirect and typically require inferring some type of relationship between an observed scattering signal and precipitation at the surface. GPM, with collocated radiometer and dual-frequency radar, is an excellent tool for tackling this problem and improving global retrievals. In the years following the launch of the GPM Core Observatory satellite, physically based passive microwave retrieval of precipitation over land continues to be challenging. Validation efforts suggest that the operational GPM passive microwave algorithm, the Goddard profiling algorithm (GPROF), tends to overestimate precipitation at the low (<5 mm h −1 ) end of the distribution over land. In this work, retrieval sensitivities to dynamic surface conditions are explored through enhancement of the algorithm with dynamic, retrieved information from a GPM-derived optimal estimation scheme. The retrieved parameters describing surface and background characteristics replace current static or ancillary GPROF information including emissivity, water vapor, and snow cover. Results show that adding this information decreases probability of false detection by 50% and, most importantly, the enhancements with retrieved parameters move the retrieval away from dependence on ancillary datasets and lead to improved physical consistency. 

Abstract This study focuses on the ability of the Global Precipitation Measurement (GPM) passive microwave sensors to detect and provide quantitative precipitation estimates (QPE) for extreme lake-effect snowfall events over the U.S. lower Great Lakes region. GPM Microwave Imager (GMI) high-frequency channels can clearly detect intense shallow convective snowfall events. However, GMI Goddard Profiling (GPROF) QPE retrievals produce inconsistent results when compared with the Multi-Radar Multi-Sensor (MRMS) ground-based radar reference dataset. While GPROF retrievals adequately capture intense snowfall rates and spatial patterns of one event, GPROF systematically underestimates intense snowfall rates in another event. Furthermore, GPROF produces abundant light snowfall rates that do not accord with MRMS observations. Ad hoc precipitation-rate thresholds are suggested to partially mitigate GPROF’s overproduction of light snowfall rates. The sensitivity and retrieval efficiency of GPROF to key parameters (2-m temperature, total precipitable water, and background surface type) used to constrain the GPROF a priori retrieval database are investigated. Results demonstrate that typical lake-effect snow environmental and surface conditions, especially coastal surfaces, are underpopulated in the database and adversely affect GPROF retrievals. For the two presented case studies, using a snow-cover a priori database in the locations originally deemed as coastline improves retrieval. This study suggests that it is particularly important to have more accurate GPROF surface classifications and better representativeness of the a priori databases to improve intense lake-effect snow detection and retrieval performance. 

Abstract Understanding the physical drivers of seasonal hydroclimatic variability and improving predictive skill remains a challenge with important socioeconomic and environmental implications for many regions around the world. Physics-based deterministic models show limited ability to predict precipitation as the lead time increases, due to imperfect representation of physical processes and incomplete knowledge of initial conditions. Similarly, statistical methods drawing upon established climate teleconnections have low prediction skill due to the complex nature of the climate system. Recently, promising data-driven approaches have been proposed, but they often suffer from overparameterization and overfitting due to the short observational record, and they often do not account for spatiotemporal dependencies among covariates (i.e., predictors such as sea surface temperatures). This study addresses these challenges via a predictive model based on a graph-guided regularizer that simultaneously promotes similarity of predictive weights for highly correlated covariates and enforces sparsity in the covariate domain. This approach both decreases the effective dimensionality of the problem and identifies the most predictive features without specifying them a priori. We use large ensemble simulations from a climate model to construct this regularizer, reducing the structural uncertainty in the estimation. We apply the learned model to predict winter precipitation in the southwestern United States using sea surface temperatures over the entire Pacific basin, and demonstrate its superiority compared to other regularization approaches and statistical models informed by known teleconnections. Our results highlight the potential to combine optimally the space–time structure of predictor variables learned from climate models with new graph-based regularizers to improve seasonal prediction. 

Abstract Spectral PCA (sPCA), in contrast to classical PCA, offers the advantage of identifying organized spatiotemporal patterns within specific frequency bands and extracting dynamical modes. However, the unavoidable trade-off between frequency resolution and robustness of the PCs leads to high sensitivity to noise and overfitting, which limits the interpretation of the sPCA results. We propose herein a simple nonparametric implementation of sPCA using the continuous analytic Morlet wavelet as a robust estimator of the cross-spectral matrices with good frequency resolution. To improve the interpretability of the results, especially when several modes of similar amplitude exist within the same frequency band, we propose a rotation of the complex-valued eigenvectors to optimize their spatial regularity (smoothness). The developed method, called rotated spectral PCA (rsPCA), is tested on synthetic data simulating propagating waves and shows impressive performance even with high levels of noise in the data. Applied to global historical geopotential height (GPH) and sea surface temperature (SST) daily time series, the method accurately captures patterns of atmospheric Rossby waves at high frequencies (3–60-day periods) in both GPH and SST and El Niño–Southern Oscillation (ENSO) at low frequencies (2–7-yr periodicity) in SST. At high frequencies the rsPCA successfully unmixes the identified waves, revealing spatially coherent patterns with robust propagation dynamics. 

Abstract Precipitation estimation based on passive microwave (MW) observations from low-Earth-orbiting satellites is one of the essential variables for understanding the global climate. However, almost all validation studies for such precipitation estimation have focused only on the surface precipitation rate. This study investigates the vertical precipitation profiles estimated by two passive MW-based retrieval algorithms, i.e., the emissivity principal components (EPC) algorithm and the Goddard profiling algorithm (GPROF). The passive MW-based condensed water content profiles estimated from the Global Precipitation Measurement Microwave Imager (GMI) are validated using the GMI + Dual-Frequency Precipitation Radar combined algorithm as the reference product. It is shown that the EPC generally underestimates the magnitude of the condensed water content profiles, described by the mean condensed water content, by about 20%–50% in the middle-to-high latitudes, while GPROF overestimates it by about 20%–50% in the middle-to-high latitudes and more than 50% in the tropics. Part of the EPC magnitude biases is associated with the representation of the precipitation type (i.e., convective and stratiform) in the retrieval algorithm. This suggests that a separate technique for precipitation type identification would aid in mitigating these biases. In contrast to the magnitude of the profile, the profile shapes are relatively well represented by these two passive MW-based retrievals. The joint analysis between the estimation performances of the vertical profiles and surface precipitation rate shows that the physically reasonable connections between the surface precipitation rate and the associated vertical profiles are achieved to some extent by the passive MW-based algorithms. 

Abstract Global gridded precipitation products have proven essential for many applications ranging from hydrological modeling and climate model validation to natural hazard risk assessment. They provide a global picture of how precipitation varies across time and space, specifically in regions where ground-based observations are scarce. While the application of global precipitation products has become widespread, there is limited knowledge on how well these products represent the magnitude and frequency of extreme precipitation—the key features in triggering flood hazards. Here, five global precipitation datasets (MSWEP, CFSR, CPC, PERSIANN-CDR, and WFDEI) are compared to each other and to surface observations. The spatial variability of relatively high precipitation events (tail heaviness) and the resulting discrepancy among datasets in the predicted precipitation return levels were evaluated for the time period 1979–2017. The analysis shows that 1) these products do not provide a consistent representation of the behavior of extremes as quantified by the tail heaviness, 2) there is strong spatial variability in the tail index, 3) the spatial patterns of the tail heaviness generally match the Köppen–Geiger climate classification, and 4) the predicted return levels for 100 and 1000 years differ significantly among the gridded products. More generally, our findings reveal shortcomings of global precipitation products in representing extremes and highlight that there is no single global product that performs best for all regions and climates. 

Abstract The Global Precipitation Measurement (GPM) constellation of spaceborne sensors provides a variety of direct and indirect measurements of precipitation processes. Such observations can be employed to derive spatially and temporally consistent gridded precipitation estimates either via data-driven retrieval algorithms or by assimilation into physically based numerical weather models. We compare the data-driven Integrated Multisatellite Retrievals for GPM (IMERG) and the assimilation-enabled NASA-Unified Weather Research and Forecasting (NU-WRF) model against Stage IV reference precipitation for four major extreme rainfall events in the southeastern United States using an object-based analysis framework that decomposes gridded precipitation fields into storm objects. As an alternative to conventional “grid-by-grid analysis,” the object-based approach provides a promising way to diagnose spatial properties of storms, trace them through space and time, and connect their accuracy to storm types and input data sources. The evolution of two tropical cyclones are generally captured by IMERG and NU-WRF, while the less organized spatial patterns of two mesoscale convective systems pose challenges for both. NU-WRF rain rates are generally more accurate, while IMERG better captures storm location and shape. Both show higher skill in detecting large, intense storms compared to smaller, weaker storms. IMERG’s accuracy depends on the input microwave and infrared data sources; NU-WRF does not appear to exhibit this dependence. Findings highlight that an object-oriented view can provide deeper insights into satellite precipitation performance and that the satellite precipitation community should further explore the potential for “hybrid” data-driven and physics-driven estimates in order to make optimal usage of satellite observations. 

Abstract This study presents the Precipitation Estimation from Remotely Sensed Information Using Artificial Neural Networks–Dynamic Infrared Rain Rate (PDIR-Now) near-real-time precipitation dataset. This dataset provides hourly, quasi-global, infrared-based precipitation estimates at 0.04° × 0.04° spatial resolution with a short latency (15–60 min). It is intended to supersede the PERSIANN–Cloud Classification System (PERSIANN-CCS) dataset previously produced as the near-real-time product of the PERSIANN family. We first provide a brief description of the algorithm’s fundamentals and the input data used for deriving precipitation estimates. Second, we provide an extensive evaluation of the PDIR-Now dataset over annual, monthly, daily, and subdaily scales. Last, the article presents information on the dissemination of the dataset through the Center for Hydrometeorology and Remote Sensing (CHRS) web-based interfaces. The evaluation, conducted over the period 2017–18, demonstrates the utility of PDIR-Now and its improvement over PERSIANN-CCS at all temporal scales. Specifically, PDIR-Now improves the estimation of rain/no-rain days as demonstrated by a critical success index (CSI) of 0.53 compared to 0.47 of PERSIANN-CCS. In addition, PDIR-Now improves the estimation of seasonal and diurnal cycles of precipitation as well as regional precipitation patterns erroneously estimated by PERSIANN-CCS. Finally, an evaluation is carried out to examine the performance of PDIR-Now in capturing two extreme events, Hurricane Harvey and a cluster of summer thunderstorms that occurred over the Netherlands, where it is shown that PDIR-Now adequately represents spatial precipitation patterns as well as subdaily precipitation rates with a correlation coefficient (CORR) of 0.64 for Hurricane Harvey and 0.76 for the Netherlands thunderstorms.