Search Results

You are looking at 1 - 10 of 17 items for

  • Author or Editor: Efi Foufoula-Georgiou x
  • Refine by Access: All Content x
Clear All Modify Search
Praveen Kumar and Efi Foufoula-Georgiou

Abstract

In this paper the authors study how zero intermittency in spatial rainfall, as described by the fraction of area covered by rainfall, changes with spatial scale of rainfall measurement or representation. A statistical measure of intermittency that describes the size distribution of “voids” (nonrainy areas imbedded inside rainy areas) as a function of scale is also introduced. Morphological algorithms are proposed for reconstructing rainfall intermittency at fine scales given the intermittency at coarser scales. These algorithms are envisioned to be useful in hydroclimatological studies where the rainfall spatial variability at the subgrid scale needs to be reconstructed from the results of synoptic- or mesoscale meteorological numerical models. The developed methodologies are demonstrated and tested using data from a severe springtime midlatitude squall line and a mild midlatitude winter storm monitored by a meteorological radar in Norman, Oklahoma.

Full access
Praveen Kumar and Efi Foufoula-Georgiou

Abstract

It has been observed that the finite-dimensional distribution functions of rainfall cannot obey simple scaling laws due to rainfall intermittency (mixed distribution with an atom at zero) and the probability of rainfall being an increasing function of area. Although rainfall fluctuations do not suffer these limitations, it is interesting to note that very few attempts have been made to study them in terms of their self-similarity characteristics. This is due to the lack of unambiguous definition of fluctuations in multidimensions. This paper shows that wavelet transforms offer a convenient and consistent method for the decomposition of inhomogeneous and anisotropic rainfall fields in two dimensions and that the components of this decomposition can be looked at as fluctuations of the rainfall field. It is also shown that under some mild assumptions, the component fields can be treated as homogeneous and thus are amenable to second-order analysis, which can provide useful insight into the nature of the process. The fact that wavelet transforms are a space-scale method also provides a convenient tool to study scaling characteristics of the process. Orthogonal wavelets are used, and these properties are investigated for a squall-line storm to study the presence of self-similarity.

Full access
Clément Guilloteau and Efi Foufoula-Georgiou

Abstract

The quantitative estimation of precipitation from orbiting passive microwave imagers has been performed for more than 30 years. The development of retrieval methods consists of establishing physical or statistical relationships between the brightness temperatures (TBs) measured at frequencies between 5 and 200 GHz and precipitation. Until now, these relationships have essentially been established at the “pixel” level, associating the average precipitation rate inside a predefined area (the pixel) to the collocated multispectral radiometric measurement. This approach considers each pixel as an independent realization of a process and ignores the fact that precipitation is a dynamic variable with rich multiscale spatial and temporal organization. Here we propose to look beyond the pixel values of the TBs and show that useful information for precipitation retrieval can be derived from the variations of the observed TBs in a spatial neighborhood around the pixel of interest. We also show that considering neighboring information allows us to better handle the complex observation geometry of conical-scanning microwave imagers, involving frequency-dependent beamwidths, overlapping fields of view, and large Earth incidence angles. Using spatial convolution filters, we compute “nonlocal” radiometric parameters sensitive to spatial patterns and scale-dependent structures of the TB fields, which are the “geometric signatures” of specific precipitation structures such as convective cells. We demonstrate that using nonlocal radiometric parameters to enrich the spectral information associated to each pixel allows for reduced retrieval uncertainty (reduction of 6%–11% of the mean absolute retrieval error) in a simple k-nearest neighbors retrieval scheme.

Open access
Deborah K. Nykanen, Efi Foufoula-Georgiou, and William M. Lapenta

Abstract

A coupled modeling framework is used in this study to investigate the effect of subgrid-scale rainfall variability on the spatial structure of the evolving storm and on other surface variables and water and energy fluxes. The Fifth-Generation Pennsylvania State University–National Center for Atmospheric Research Mesoscale Model coupled with the Biosphere–Atmosphere Transfer Scheme is combined with a dynamical/statistical scheme for statistically downscaling rainfall. Model simulations with and without including subgrid-scale rainfall variability are compared at the grid scale to quantify the propagation of small-scale rainfall heterogeneities through the nonlinear land–atmosphere system. It was found that including subgrid-scale rainfall variability (here on the order of 3 km) affects the spatial organization of the storm system itself, surface temperature, soil moisture, and sensible and latent heat fluxes. These effects were found to occur at spatial scales much larger than the scale at which rainfall variability was prescribed, illustrating the pronounced nonlinear spatial dynamics of the land–atmosphere system and its important role on hydrometeorological predictions.

Full access
Clément Guilloteau, Efi Foufoula-Georgiou, and Christian D. Kummerow

Abstract

The constellation of spaceborne passive microwave (MW) sensors, coordinated under the framework of the Precipitation Measurement Missions international agreement, continuously produces observations of clouds and precipitation all over the globe. The Goddard profiling algorithm (GPROF) is designed to infer the instantaneous surface precipitation rate from the measured MW radiances. The last version of the algorithm (GPROF-2014)—the product of more than 20 years of algorithmic development, validation, and improvement—is currently used to estimate precipitation rates from the microwave imager GMI on board the GPM core satellite. The previous version of the algorithm (GPROF-2010) was used with the microwave imager TMI on board TRMM. In this paper, TMI-GPROF-2010 estimates and GMI-GPROF-2014 estimates are compared with coincident active measurements from the Precipitation Radar on board TRMM and the Dual-Frequency Precipitation Radar on board GPM, considered as reference products. The objective is to assess the improvement of the GPM-era microwave estimates relative to the TRMM-era estimates and diagnose regions where continuous improvement is needed. The assessment is oriented toward estimating the “effective resolution” of the MW estimates, that is, the finest scale at which the retrieval is able to accurately reproduce the spatial variability of precipitation. A wavelet-based multiscale decomposition of the radar and passive microwave precipitation fields is used to formally define and assess the effective resolution. It is found that the GPM-era MW retrieval can resolve finer-scale spatial variability over oceans than the TRMM-era retrieval. Over land, significant challenges exist, and this analysis provides useful diagnostics and a benchmark against which future retrieval algorithm improvement can be assessed.

Open access
Clément Guilloteau, Efi Foufoula-Georgiou, Christian D. Kummerow, and Veljko Petković

Abstract

The scattering of microwaves at frequencies between 50 and 200 GHz by ice particles in the atmosphere is an essential element in the retrieval of instantaneous surface precipitation from spaceborne passive radiometers. This paper explores how the variable distribution of solid and liquid hydrometeors in the atmospheric column over land surfaces affects the brightness temperature (TB) measured by GMI at 89 GHz through the analysis of Dual-Frequency Precipitation Radar (DPR) reflectivity profiles along the 89-GHz beam. The objective is to refine the statistical relations between observed TBs and surface precipitation over land and to define their limits. As GMI is scanning with a 53° Earth incident angle, the observed atmospheric volume is actually not a vertical column, which may lead to very heterogeneous and seemingly inconsistent distributions of the hydrometeors inside the beam. It is found that the 89-GHz TB is mostly sensitive to the presence of ice hydrometeors several kilometers above the 0°C isotherm, up to 10 km above the 0°C isotherm for the deepest convective systems, but is a modest predictor of the surface precipitation rate. To perform a precise mapping of atmospheric ice, the altitude of the individual ice clusters must be known. Indeed, if variations in the altitude of ice are not accounted for, then the high incident angle of GMI causes a horizontal shift (parallax shift) between the estimated position of the ice clusters and their actual position. We show here that the altitude of ice clusters can be derived from the 89-GHz TB itself, allowing for correction of the parallax shift.

Open access
Ardeshir M. Ebtehaj, Rafael L. Bras, and Efi Foufoula-Georgiou

Abstract

Using satellite measurements in microwave bands to retrieve precipitation over land requires proper discrimination of the weak rainfall signals from strong and highly variable background Earth surface emissions. Traditionally, land retrieval methods rely on a weak signal of rainfall scattering on high-frequency channels and make use of empirical thresholding and regression-based techniques. Because of the increased surface signal interference, retrievals over radiometrically complex land surfaces—snow-covered lands, deserts, and coastal areas—are particularly challenging for this class of retrieval techniques. This paper evaluates the results by the recently proposed Shrunken Locally Linear Embedding Algorithm for Retrieval of Precipitation (ShARP) using data from the Tropical Rainfall Measuring Mission (TRMM) satellite. The study focuses on a radiometrically complex region, partly covering the Tibetan highlands, Himalayas, and Ganges–Brahmaputra–Meghna River basins, which is unique in terms of its diverse land surface radiation regime and precipitation type, within the TRMM domain. Promising results are presented using ShARP over snow-covered land surfaces and in the vicinity of coastlines, in comparison with the land rainfall retrievals of the standard TRMM 2A12, version 7, product. The results show that ShARP can significantly reduce the rainfall overestimation due to the background snow contamination and markedly improve detection and retrieval of rainfall in the vicinity of coastlines. During the calendar year 2013, compared to TRMM 2A25, it is demonstrated that over the study domain the root-mean-square difference can be reduced up to 38% annually, while the improvement can reach up to 70% during the cold months of the year.

Full access
Clément Guilloteau, Antonios Mamalakis, Lawrence Vulis, Phong V. V. Le, Tryphon T. Georgiou, and Efi Foufoula-Georgiou

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.

Open access
Jamie L. Smedsmo, Efi Foufoula-Georgiou, Venugopal Vuruputur, Fanyou Kong, and Kelvin Droegemeier

Abstract

An understanding of the vertical structure of clouds is important for remote sensing of precipitation from space and for the parameterization of cloud microphysics in numerical weather prediction (NWP) models. The representation of cloud hydrometeor profiles in high-resolution NWP models has direct applications in inversion-type precipitation retrieval algorithms [e.g., the Goddard profiling (GPROF) algorithm used for retrieval of precipitation from passive microwave sensors] and in quantitative precipitation forecasting. This study seeks to understand how the vertical structure of hydrometeors (liquid and frozen water droplets in a cloud) produced by high-resolution NWP models with explicit microphysics, henceforth referred to as cloud-resolving models (CRMs), compares to observations. Although direct observations of 3D hydrometeor fields are not available, comparisons of modeled and observed radar echoes can provide some insight into the vertical structure of hydrometeors and, in turn, into the ability of CRMs to produce precipitation structures that are consistent with observations. Significant differences are revealed between the vertical structure of observed and modeled clouds of a severe midlatitude storm over Texas for which the surface precipitation is reasonably well captured. Possible reasons for this discrepancy are presented, and the need for future research is highlighted.

Full access
Clement Guilloteau, Efi Foufoula-Georgiou, Pierre Kirstetter, Jackson Tan, and George J. Huffman

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 that 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 multisatellite products (CMORPH, GSMaP, IMERG-Early, IMERG-Final, and PERSIANN-CCS) over a period of 2 years over the southeastern United States. 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.

Open access