skip to main content

Title: Beyond the Pixel: Using Patterns and Multiscale Spatial Information to Improve the Retrieval of Precipitation from Spaceborne Passive Microwave Imagers
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.  more » « less
Award ID(s):
1839441 1839336
Author(s) / Creator(s):
Date Published:
Journal Name:
Journal of Atmospheric and Oceanic Technology
Page Range / eLocation ID:
1571 to 1591
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. We present a multi-sensor Bayesian passive microwave retrieval algorithm for flood inundation mapping at high spatial and temporal resolutions. The algorithm takes advantage of observations from multiple sensors in optical, short-infrared, and microwave bands, thereby allowing for detection and mapping of the sub-pixel fraction of inundated areas under almost all-sky conditions. The method relies on a nearest-neighbor search and a modern sparsity-promoting inversion method that make use of an a priori dataset in the form of two joint dictionaries. These dictionaries contain almost overlapping observations by the Special Sensor Microwave Imager and Sounder (SSMIS) on board the Defense Meteorological Satellite Program (DMSP) F17 satellite and the Moderate Resolution Imaging Spectroradiometer (MODIS) on board the Aqua and Terra satellites. Evaluation of the retrieval algorithm over the Mekong Delta shows that it is capable of capturing to a good degree the inundation diurnal variability due to localized convective precipitation. At longer timescales, the results demonstrate consistency with the ground-based water level observations, denoting that the method is properly capturing inundation seasonal patterns in response to regional monsoonal rain. The calculated Euclidean distance, rank-correlation, and also copula quantile analysis demonstrate a good agreement between the outputs of the algorithm and the observed water levels at monthly and daily timescales. The current inundation products are at a resolution of 12.5 km and taken twice per day, but a higher resolution (order of 5 km and every 3 h) can be achieved using the same algorithm with the dictionary populated by the Global Precipitation Mission (GPM) Microwave Imager (GMI) products. 
    more » « less
  3. Abstract

    Additive manufacturing, no longer reserved exclusively for prototyping components, can create parts with complex geometries and locally tailored properties. For example, multiple homogenous material sources can be used in different regions of a print or be mixed during printing to define properties locally. Additionally, heterogeneous composites provide an opportunity for another level of tuning properties through processing. For example, within particulate-filled polymer matrix composites before curing, the presence of an applied electric and/or magnetic fields can reorient filler particles and form hierarchical structures depending on the fields applied. Control of particle organization is important because effective material properties are highly dependent on the distribution of filler material within composites once cured. While previous work in homogenization and effective medium theories have determined properties based upon ideal analytic distributions of particle orientations and spatial location, this work expands upon these methods generating discrete distributions from quasi-Monte Carlo simulations of the electromagnetic processing event. Results of simulations provide predicted microarchitectures from which effective properties are determined via computational homogenization.

    These particle dynamics simulations account for dielectric and magnetic forces and torques in addition to hydrodynamic forces and hard particle separation. As such, the distributions generated are processing field dependent. The effective properties for a composite represented by this distribution are determined via computational homogenization using finite element analysis (FEA). This provides a path from constituents, through processing parameters to effective material properties. In this work, we use these simulations in conjunction with a multi-objective optimization scheme to resolve the relationships between processing conditions and effective properties, to inform field-assisted additive manufacturing processes.

    The constituent set providing the largest range of properties can be found using optimization techniques applied to the aforementioned simulation framework. This key information provides a recipe for tailoring properties for additive manufacturing design and production. For example, our simulation results show that stiffness for a 10% filler volume fraction can increase by 34% when aligned by an electric field as compared to a randomly distributed composite. The stiffness of this aligned sample is also 29% higher in the direction of the alignment than perpendicular to it, which only differs by 5% from the random case [1]. Understanding this behavior and accurately predicting composite properties is key to producing field processed composites and prints. Material property predictions compare favorably to effective medium theory and experimentation with trends in elastic and magnetic effective properties demonstrating the same anisotropic behavior as a result of applied field processing. This work will address the high computational expense of physics simulation based objective functions by using efficient algorithms and data structures. We will present an optimization framework using nested gradient searches for micro barium hexaferrite particles in a PDMS matrix, optimizing on composite magnetization to determine the volume fraction of filler that will provide the largest range of properties by varying the applied electric and magnetic fields.

    more » « less
  4. null (Ed.)
    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. 
    more » « less
  5. Abstract. Reconstructions of past temperature and precipitation are fundamental to modeling the Greenland Ice Sheet and assessing its sensitivity to climate. Paleoclimate information is sourced from proxy records and climate-model simulations; however, the former are spatially incomplete while the latter are sensitive to model dynamics and boundary conditions. Efforts to combine these sources of information to reconstruct spatial patterns of Greenland climate over glacial–interglacial cycles have been limited by assumptions of fixed spatial patterns and a restricted use of proxy data. We avoid these limitations by using paleoclimate data assimilation to create independent reconstructions of mean-annual temperature and precipitation for the last 20 000 years. Our method uses oxygen isotope ratios of ice and accumulation rates from long ice-core records and extends this information to all locations across Greenland using spatial relationships derived from a transient climate-model simulation. Standard evaluation metrics for this method show that our results capture climate at locations without ice-core records. Our results differ from previous work in the reconstructed spatial pattern of temperature change during abrupt climate transitions; this indicates a need for additional proxy data and additional transient climate-model simulations. We investigate the relationship between precipitation and temperature, finding that it is frequency dependent and spatially variable, suggesting that thermodynamic scaling methods commonly used in ice-sheet modeling are overly simplistic. Our results demonstrate that paleoclimate data assimilation is a useful tool for reconstructing the spatial and temporal patterns of past climate on timescales relevant to ice sheets. 
    more » « less