Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial variability”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Evolution of upwelling fronts off the coast of Portugal

We have analyzed the spatial and temporal variability of chlorophyll in the Atlantic Ocean off the Portugal coast using satellite imagery from 1979 to 1986 in conjunction with wind data. The spatial patterns in the imagery show the typical inshore-offshore and east-west gradients, while the temporal patterns show seasonal high winter values and low summer values. In the satellite imagery for the period 1979-1986, stable mesoscale and large scale patterns recurred despite interannual differences in absolute satellite-derived pigment concentrations.

Source record↗
NASA NTRS1989

Interannual variability in phytoplankton pigment distribution during the spring transition along the west coast of North America

A 5-year time series of coastal zone color scanner imagery (1980-1983, 1986) is used to examine changes in the large-scale pattern of chlorophyll pigment concentration coincident with the spring transition in winds and currents along the west coast of North America. The data show strong interannual variability in the timing and spatial patterns of pigment concentration at the time of the transition event. Interannual variability in the response of pigment concentration to the spring transition appears to be a function of spatial and temporal variability in vertical nutrient flux induced by wind mixing and/or the upwelling initiated at the time of the transition. Interannual differences in the mixing regime are illustrated with a one-dimensional mixing model.

Thomas, A. C.↗
DOE OSTI2023

A seamless approach for evaluating climate models across spatial scales

In regions of the world where topography varies significantly with distance, most global climate models (GCMs) have spatial resolutions that are too coarse to accurately simulate key meteorological variables that are influenced by topography, such as clouds, precipitation, and surface temperatures. One approach to tackle this challenge is to run climate models of sufficiently high resolution in those topographically complex regions such as the North American Regionally Refined Model (NARRM) subset of the Department of Energy’s (DOE) Energy Exascale Earth System Model version 2 (E3SM v2). Although high-resolution simulations are expected to provide unprecedented details of atmospheric processes, running models at such high resolutions remains computationally expensive compared to lower-resolution models such as the E3SM Low Resolution (LR). Moreover, because regionally refined and high-resolution GCMs are relatively new, there are a limited number of observational datasets and frameworks available for evaluating climate models with regionally varying spatial resolutions. As such, we developed a new framework to quantify the added value of high spatial resolution in simulating precipitation over the contiguous United States (CONUS). To determine its viability, we applied the framework to two model simulations and an observational dataset. We first remapped all the data into Hierarchical Equal-Area Iso-Latitude Pixelization (HEALPix) pixels. HEALPix offers several mathematical properties that enable seamless evaluation of climate models across different spatial resolutions including its equal-area and partitioning properties. The remapped HEALPix-based data are used to show how the spatial variability of both observed and simulated precipitation changes with resolution increases. This study provides valuable insights into the requirements for achieving accurate simulations of precipitation patterns over the CONUS. It highlights the importance of allocating sufficient computational resources to run climate models at higher temporal and spatial resolutions to capture spatial patterns effectively. Furthermore, the study demonstrates the effectiveness of the HEALPix framework in evaluating precipitation simulations across different spatial resolutions. This framework offers a viable approach for comparing observed and simulated data when dealing with datasets of varying spatial resolutions. By employing this framework, researchers can extend its usage to other climate variables, datasets, and disciplines that require comparing datasets with different spatial resolutions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Impacts of spatial heterogeneity of anthropogenic aerosol emissions in a regionally refined global aerosol–climate model

Abstract. Emissions of anthropogenic aerosol and their precursors are often prescribed in global aerosol models. Most of these emissions are spatially heterogeneous at model grid scales. When remapped from low-resolution data, the spatial heterogeneity in emissions can be lost, leading to large errors in the simulation. It can also cause the conservation problem if non-conservative remapping is used. The default anthropogenic emission treatment in the Energy Exascale Earth System Model (E3SM) is subject to both problems. In this study, we introduce a revised emission treatment for the E3SM Atmosphere Model (EAM) that ensures conservation of mass fluxes and preserves the original emission heterogeneity at the model-resolved grid scale. We assess the error estimates associated with the default emission treatment and the impact of improved heterogeneity and mass conservation in both globally uniform standard-resolution (∼ 165 km) and regionally refined high-resolution (∼ 42 km) simulations. The default treatment incurs significant errors near the surface, particularly over sharp emission gradient zones. Much larger errors are observed in high-resolution simulations. It substantially underestimates the aerosol burden, surface concentration, and aerosol sources over highly polluted regions, while it overestimates these quantities over less-polluted adjacent areas. Large errors can persist at higher elevation for daily mean estimates, which can affect aerosol extinction profiles and aerosol optical depth (AOD). We find that the revised treatment significantly improves the accuracy of the aerosol emissions from surface and elevated sources near sharp spatial gradient regions, with significant improvement in the spatial heterogeneity and variability of simulated surface concentration in high-resolution simulations. In the next-generation E3SM running at convection-permitting scales where the resolved spatial heterogeneity is significantly increased, the revised emission treatment is expected to better represent the aerosol emissions as well as their lifecycle and impacts on climate.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1987

The atmosphere as a research window - Signal or noise

The interpretation of measurements by instruments with FOV that are large enough to contain pronounced radiometric variations is discussed, focusing on the retrieval of cloud properties from satellite radiance measurements. Examples are presented for three retrieval approaches. The use of a physical principle to establish a relationship between the measurement and the desired parameter is shown for cirrus-cloud optical depth and precipitation from passive microwave measurements. Examples of the use of inherent variability of the atmosphere and the statistical character of this variability are the spatial covariance method and the threshold method. Also, two examples using a hybrid of these two methods are presented. It is shown that the natural variability (noise) within the instrument FOV biases the relationship between the measurement and the parameter of interest which would be predicted by simple one-dimensional radiative-transfer theory. It is suggested that the noise of the measurement also contains signal, forming the basis for hybrid methods.

Stephens, Graeme L.↗
DOE OSTI2022

Increased Variability of Biomass Burning Emissions in CMIP6 Amplifies Hydrologic Cycle in the CESM2 Large Ensemble

Abstract Historical simulations performed for the Coupled Model Intercomparison Project Phase 6 used biomass burning emissions between 1997 and 2014 containing higher spatial and temporal variability compared to emission inventories specified for earlier years, and compared to emissions used in previous (e.g., CMIP5) simulation intercomparisons. Using the Community Earth System Model version 2 Large Ensemble, we show this increased biomass burning emissions variability leads to amplification of the hydrologic cycle poleward of 40°N. Notably, the high variability of biomass burning emissions leads to increased latent heat fluxes, column‐integrated precipitable water, and precipitation. Greater ocean heat uptake, weaker meridional energy transport from the tropics, greater atmospheric shortwave and longwave absorption, and lower relative humidity act to moderate this hydrologic cycle amplification. Our results suggest it is not only the secular changes (on multidecadal timescales) in biomass burning emissions that impact the hydrologic cycle, but also the shorter timescale variability in emissions.

Heyblom, Kyle B.↗
NASA NTRS1992

Plot-scale field experiment of surface hydrologic processes with EOS implications

Plot-scale hydrologic field studies were initiated at NASA Marshall Space Flight Center to a) investigate the spatial and temporal variability of surface and subsurface hydrologic processes, particularly as affected by vegetation, and b) develop experimental techniques and associated instrumentation methodology to study hydrologic processes at increasingly large spatial scales. About 150 instruments, most of which are remotely operated, have been installed at the field site to monitor ground atmospheric conditions, precipitation, interception, soil-water status, and energy flux. This paper describes the nature of the field experiment, instrumentation and sampling rationale, and presents preliminary findings.

Laymon, Charles A.↗
NASA NTRS1972

Dynamic response of a semi-infinite elastic cylinder containing an acoustic medium.

The solutions developed in this paper are based on the finite Hankel-Fourier-Laplace transform. The Love-Timoshenko shell equations are reduced through the application of the Fourier-Laplace transforms to the axial space and time variables. The equation of motion for the acoustic fluid contained within the shell are reduced through the application of the finite Hankel-Fourier-Laplace transform to the radial, axial spatial, and time variables. The boundary conditions at the fixed end of the semi-infinite shell are met through the application of loadings symmetric with respect to the origin and the application of a ring loading at the origin which is chosen so as to make the radial deflection there zero. Expressions are found for the transforms of the axial, tangential, and radial shell displacements, the axial, tangential, and radial fluid velocities, and the radiated fluid pressure. Numerical inversion of the Fourier-Laplace transform is accomplished in terms of a series of ultraspherical polynomials and Gauss-Hermite or Gauss-Laguerre quadrature.

Carpenter, R. F., III↗
NASA NTRS2001

Bio-Optical Measurements at Ocean Boundaries in Support of SIMBIOS

The equatorial Pacific is a major component of global biogeochemical cycles, due to upwelling that occurs from the coast of South America to beyond 180 deg. This upwelling has significant implications for global CO2 fluxes, as well as primary and secondary production. In addition, this region of the world's oceans represents a large oceanic province over which validation data for Sea-Viewing Wide Field-of-view Sensor (SeaWiFS) are necessary. This project consists of a mooring program and supporting cruise-based measurements aimed at quantifying the spectrum of biological and chemical variability in the equatorial Pacific and obtaining validation data for SeaWiFS. The project has the following general objectives: (1) to understand the relationships between physical forcing, primary production, nutrient supply and the exchange of carbon dioxide between ocean and atmosphere in the equatorial Pacific; (2) to describe the biological and chemical responses to climate and ocean variability; (3) to describe the spatial, seasonal and inter-annual variability in near surface plant pigments, primary production, carbon dioxide and nutrient distributions; and (4) to obtain near real-time bio-optical measurements for validation of SeaWiFS and subsequent ocean color sensors.

Chavez, Francisco P.↗
NASA NTRS2016

Uncertainty Assessment of the NASA Earth Exchange Global Daily Downscaled Climate Projections (NEX-GDDP) Dataset

The NASA Earth Exchange Global Daily Downscaled Projections (NEX-GDDP) dataset is comprised of downscaled climate projections that are derived from 21 General Circulation Model (GCM) runs conducted under the Coupled Model Intercomparison Project Phase 5 (CMIP5) and across two of the four greenhouse gas emissions scenarios (RCP4.5 and RCP8.5). Each of the climate projections includes daily maximum temperature, minimum temperature, and precipitation for the periods from 1950 through 2100 and the spatial resolution is 0.25 degrees (approximately 25 km x 25 km). The GDDP dataset has received warm welcome from the science community in conducting studies of climate change impacts at local to regional scales, but a comprehensive evaluation of its uncertainties is still missing. In this study, we apply the Perfect Model Experiment framework (Dixon et al. 2016) to quantify the key sources of uncertainties from the observational baseline dataset, the downscaling algorithm, and some intrinsic assumptions (e.g., the stationary assumption) inherent to the statistical downscaling techniques. We developed a set of metrics to evaluate downscaling errors resulted from bias-correction ("quantile-mapping"), spatial disaggregation, as well as the temporal-spatial non-stationarity of climate variability. Our results highlight the spatial disaggregation (or interpolation) errors, which dominate the overall uncertainties of the GDDP dataset, especially over heterogeneous and complex terrains (e.g., mountains and coastal area). In comparison, the temporal errors in the GDDP dataset tend to be more constrained. Our results also indicate that the downscaled daily precipitation also has relatively larger uncertainties than the temperature fields, reflecting the rather stochastic nature of precipitation in space. Therefore, our results provide insights in improving statistical downscaling algorithms and products in the future.

climate projection↗
NASA NTRS2016

Uncertainty Assessment of the NASA Earth Exchange Global Daily Downscaled Climate Projections (NEX-GDDP) Dataset

The NASA Earth Exchange Global Daily Downscaled Projections (NEX-GDDP) dataset is comprised of downscaled climate projections that are derived from 21 General Circulation Model (GCM) runs conducted under the Coupled Model Intercomparison Project Phase 5 (CMIP5) and across two of the four greenhouse gas emissions scenarios (RCP4.5 and RCP8.5). Each of the climate projections includes daily maximum temperature, minimum temperature, and precipitation for the periods from 1950 through 2100 and the spatial resolution is 0.25 degrees (approximately 25 km by 25 km). The GDDP dataset has received warm welcome from the science community in conducting studies of climate change impacts at local to regional scales, but a comprehensive evaluation of its uncertainties is still missing. In this study, we apply the Perfect Model Experiment framework (Dixon et al. 2016) to quantify the key sources of uncertainties from the observational baseline dataset, the downscaling algorithm, and some intrinsic assumptions (e.g., the stationary assumption) inherent to the statistical downscaling techniques. We developed a set of metrics to evaluate downscaling errors resulted from bias-correction ("quantile-mapping"), spatial disaggregation, as well as the temporal-spatial non-stationarity of climate variability. Our results highlight the spatial disaggregation (or interpolation) errors, which dominate the overall uncertainties of the GDDP dataset, especially over heterogeneous and complex terrains (e.g., mountains and coastal area). In comparison, the temporal errors in the GDDP dataset tend to be more constrained. Our results also indicate that the downscaled daily precipitation also has relatively larger uncertainties than the temperature fields, reflecting the rather stochastic nature of precipitation in space. Therefore, our results provide insights in improving statistical downscaling algorithms and products in the future.

general circulation model (GCM)↗
NASA NTRS2008

Atmospheric Parameter Climatologies from AIRS: Monitoring Short-, and Longer-Term Climate Variabilities and 'Trends'

The AIRS instrument is currently the best space-based tool to simultaneously monitor the vertical distribution of key climatically important atmospheric parameters as well as surface properties, and has provided high quality data for more than 5 years. AIRS analysis results produced at the GODDARD/DAAC, based on Versions 4 & 5 of the AIRS retrieval algorithm, are currently available for public use. Here, first we present an assessment of interrelationships of anomalies (proxies of climate variability based on 5 full years, since Sept. 2002) of various climate parameters at different spatial scales. We also present AIRS-retrievals-based global, regional and 1x1 degree grid-scale "trend"-analyses of important atmospheric parameters for this 5-year period. Note that here "trend" simply means the linear fit to the anomaly (relative the mean seasonal cycle) time series of various parameters at the above-mentioned spatial scales, and we present these to illustrate the usefulness of continuing AIRS-based climate observations. Preliminary validation efforts, in terms of intercomparisons of interannual variabilities with other available satellite data analysis results, will also be addressed. For example, we show that the outgoing longwave radiation (OLR) interannual spatial variabilities from the available state-of-the-art CERES measurements and from the AIRS computations are in remarkably good agreement. Version 6 of the AIRS retrieval scheme (currently under development) promises to further improve bias agreements for the absolute values by implementing a more accurate radiative transfer model for the OLR computations and by improving surface emissivity retrievals.

Molnar, Gyula↗
NASA NTRS1994

Spatial patterns in the length of the sea ice season in the Southern Ocean, 1979-1986

The length of the sea ice season summarizes in one number the ice coverage conditions for an individual location for an entire year. It becomes a particularly valuable variable when mapped spatially over a large area and examined for regional and interannual differences, as is done here for the Southern Ocean over the years 1979-1986, using the satellite passive microwave data of the Nimbus 7 scanning multichannel microwave radiometer. Three prominent geographic anomalies in ice season lengths occur consistently in each year of the data set, countering the general tendency toward shorter ice seasons from south to north: (1) in the Weddell Sea the tendency is toward shorter ice seasons from southwest to northeast, reflective of the cyclonic ice/atmosphere/ocean circulations in the Weddell Sea region. (2) Directly north of the Ross Ice Shelf anomalously short ice seasons occur, lasting only 245-270 days, in contrast to the perennial ice coverage at comparable latitudes in the southern Bellingshausen and Amundsen Seas and in the western Weddell Sea. The short ice season off the Ross Ice Shelf reflects the consistently early opening of the ice cover each spring, under the influence of upwelling along the continental slope and shelf and atmospheric forcing from winds blowing off the Antarctic continent. (3) In the southern Amundsen Sea, anomalously short ice seasons occur adjacent to the coast, owing to the frequent existence of coastal polynyas off the many small ice shelves bordering the sea. Least squares trends in the ice season lengths over the 1979-1986 period are highly coherent spatially, with overall trends toward shorter ice seasons in the northern Weddell and Bellingshausen seas and toward longer ice seasons in the Ross Sea, around much of East Antarctica, and in a portion of the south central Weddell Sea.

Parkinson, Claire L.↗
NASA NTRS2006

Diurnal Variability of Vertical Structure from a TRMM Passive Microwave "Virtual Radar" Retrieval

Robust description of the diurnal cycle from TRMM observations is complicated by the limitations of Low Earth Orbit (LEO) sampling; from a 'climatological' perspective, sufficient sampling must exist to control for both spatial and seasonal variability, before tackling an additional diurnal component (e.g., with 8 additional 3-hourly or 24 1-hourly bins). For documentation of vertical structure, the narrow sample swath of the TRMM Precipitation Radar limits the resolution of any of these components. A neural-network based 'virtual radar" retrieval has been trained and internally validated, using multifrequency / multipolarization passive microwave(TM1) brightness temperatures and textures parameters and lightning (LIS) observations, as inputs, and PR volumetric reflectivity as targets (outputs). By training the algorithms (essentially highly multivariate, nonlinear regressions) on a very large sample of high-quality co-located data from the center of the TRMM swath, 3D radar reflectivity and derived parameters (VIL, IWC, Echo Tops, etc.) can be retrieved across the entire TMI swath, good to 8-9% over the dynamic range of parameters. As a step in the retrieval (and as an output of the process), each TMI multifrequency pixel (at 85 GHz resolution) is classified into one of the 25 archetypal radar profile vertical structure "types", previously identified using cluster analysis. The dynamic range of retrieved vertical structure appears to have higher fidelity than the current (Version 6) experimental GPROF hydrometeor vertical structure retrievals. This is attributable to correct representation of the prior probabilities of vertical structure variability in the neural network training data, unlike the GPROF cloud-resolving model training dataset used in the V6 algorithms. The LIS lightning inputs are supplementary inputs, and a separate offline neural network has been trained to impute (predict) LIS lightning from passive-microwave-only data. The virtual radar retrieval is thus, in principle, extensible to Aqua/AMSR-E and NPOESS/CMIS passive microwave instruments. The virtual radar approach yields a threefold increase in effective sampling from the mission, albeit of lower-quality "retrieved" data, reducing the variance of local estimates by one third (or the standard deviation by-0.57). In this talk, the variance reduction is leveraged to more finely resolve global diurnal variability in both space and time (local hour).

Boccippio, Dennis J.↗
DOE OSTI2022

A physical method for downscaling land surface temperatures using surface energy balance theory

Fine-resolution land surface temperature (LST) derived from thermal infrared remote sensing images is a good indicator of surface water status and plays an essential role in the exchange of energy and water between land and atmosphere. A physical surface energy balance (SEB)-based LST downscaling method (DTsEB) is developed to downscale coarse remotely sensed thermal infrared LST products with fine-resolution visible and near-infrared data. Here, the DTsEB method is advantageous for its ability to mechanically interrelate surface variables contributing to the spatial variation of LST, to quantitatively weigh the contributions of each related variable within a physical framework, and to efficaciously avoid the subjective selection of scaling factors and the establishment of statistical regression relationships. The applicability of the DTsEB method was tested by downscaling 12 scenes of 990 m Moderate Resolution Imaging Spectroradiometer (MODIS) and aggregated Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) LST products to 90 m resolution at six overpass times between 2005 and 2015 over three 9.9 km by 9.9 km cropland (mixed by grass, tree, and built-up land) study areas. Three typical LST downscaling methods, namely the widely applied TsHARP, the later developed least median square regression downscaling (LMS) and the geographically weighted regression (GWR), were introduced for intercomparison. The results showed that the DTsEB method could more effectively reconstruct the subpixel spatial variations in LST within the coarse-resolution pixels and achieve a better downscaling accuracy than the TsHARP, LMS and GWR methods. The DTsEB method yielded, on average, root mean square errors (RMSEs) of 2.01 K and 1.42 K when applied to the MODIS datasets and aggregated ASTER datasets, respectively, which were lower than those obtained with the TsHARP method, with average RMSEs of 2.41 K and 1.71 K, the LMS method, with average RMSEs of 2.35 K and 1.63 K, and the GWR method, with average RMSEs of 2.38 K and 1.64 K, respectively. The contributions of the related surface variables to the subpixel spatial variation in the LST varied both spatially and temporally and were different from each other. In summary, the DTsEB method was demonstrated to outperform the TsHARP, LMS, and GWR methods and could be used as a good alternative for downscaling LST products from coarse to fine resolution with high robustness and accuracy.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Long-Duration High-Resolution Large-Eddy Simulations of Hurricane Laura (2020) for Infrastructure Design

This dataset provides long-duration large-eddy simulations (LES) of Hurricane Laura (2020), capturing the turbulent boundary-layer wind field during the storm's passage at horizontal grid spacings of 55.55 m and 10 m. Hurricane Laura made landfall as a Category 4 storm in southwestern Louisiana at 06:00 UTC on 27 August 2020, producing widespread damage to coastal infrastructure. The simulations characterize the spatial and temporal variability of the turbulent wind at heights spanning the lowest 300 m of the boundary layer, throughout the approach and passage of the storm's eyewall over a fixed location. To resolve the spatial heterogeneity of wind and turbulence conditions across the storm, seven independent refined LES runs are performed at radial locations spanning the storm relative to its direction of motion. The refinement positions are defined on a normalized radial coordinate r ̂=r/R, where R≈21 km is Laura's radius of maximum wind, and lie between r ̂=-1.5 and r ̂=+1.5 in ∆r ̂=0.5 increments along a line perpendicular to the storm's track at the midpoint of the simulation. Negative r ̂ corresponds to locations south-southwest of the storm center and positive r ̂ to locations north-northeast. The dataset is segmented into seven radial refinement groups.

17 WIND ENERGY↗
DOE OSTI2025

Hydrologic Model Data for the East Fork Poplar Creek Watershed Simulated with the Advanced Terrestrial Simulator (ATS): Streamflow and Network Expansion–Contraction Dynamics

This dataset supports hydrologic modeling and stream network expansion–contraction analysis for the East Fork Poplar Creek (EFPC) Watershed in Tennessee. It includes a Jupyter notebook for model setup, model configuration files, simulation outputs, and derived products used to evaluate model performance and investigate stream dynamics under varying hydrologic conditions. The dataset was generated using the Watershed Workflow Python package and the Advanced Terrestrial Simulator (ATS), enabling integrated surface–subsurface hydrologic simulations using a stream-aligned mesh. Outputs include high-resolution time series of streamflow, active network length, water table depth, and related hydrologic variables. Also included are spatially explicit stream persistency indices and classifications of reaches as perennial or non-perennial. These data facilitate reproducibility and support further research on stream intermittency and variability in network extent.The model data archive is organized in following directories:1) model_setup_inputsContains the Watershed Workflow Jupyter notebooks (accessed through any open source code editor), selected input datasets, and resulting ATS input files, including XML files (access through any open source code editor), computational mesh (.exo files can be viewed using Paraview), and meteorological forcing files (.h5 files can be accessed through h5py python package and HDFView open source software). 2) model_outputsIncludes ATS simulation outputs relevant to this study. Time series of spatially integrated or averaged variables (e.g., streamflow, water table depth) are provided as CSV files. Select spatial fields (e.g., ponded depth and water table depth) are saved as pickled Python objects to reduce file size, and can be accessed through pickle package in Python. Key geometry objects from Watershed Workflow—such as the surface mesh and river tree—are also included to support analysis of streamflow persistency and expansion–contraction dynamics. These files can also be accessed through Watershed Workflow Python package.3) model_evaluationProvides observed streamflow time series and field survey-based flow regime classifications used to evaluate model performance. Jupyter notebooks for processing ATS outputs and comparing model predictions with observations to build confidence in the model prior to scientific analysis are also included.4) Q_L_relationshipsContains workflows for generating time series of discharge, active network length, and related hydrologic variables used in the stream network expansion–contraction analysis. Includes routines for delineating baseflow-dominated periods. For each catchment, notebooks and processed data (as pickled DataFrames accessed through Pandas Python package) are provided. 5) figure_scriptsProvides the Jupyter notebooks used to generate the figures presented in the paper.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2002

A Stochastic Model of Space-Time Variability of Tropical Rainfall: I. Statistics of Spatial Averages

Global maps of rainfall are of great importance in connection with modeling of the earth s climate. Comparison between the maps of rainfall predicted by computer-generated climate models with observation provides a sensitive test for these models. To make such a comparison, one typically needs the total precipitation amount over a large area, which could be hundreds of kilometers in size over extended periods of time of order days or months. This presents a difficult problem since rain varies greatly from place to place as well as in time. Remote sensing methods using ground radar or satellites detect rain over a large area by essentially taking a series of snapshots at infrequent intervals and indirectly deriving the average rain intensity within a collection of pixels , usually several kilometers in size. They measure area average of rain at a particular instant. Rain gauges, on the other hand, record rain accumulation continuously in time but only over a very small area tens of centimeters across, say, the size of a dinner plate. They measure only a time average at a single location. In making use of either method one needs to fill in the gaps in the observation - either the gaps in the area covered or the gaps in time of observation. This involves using statistical models to obtain information about the rain that is missed from what is actually detected. This paper investigates such a statistical model and validates it with rain data collected over the tropical Western Pacific from ship borne radars during TOGA COARE (Tropical Oceans Global Atmosphere Coupled Ocean-Atmosphere Response Experiment). The model incorporates a number of commonly observed features of rain. While rain varies rapidly with location and time, the variability diminishes when averaged over larger areas or longer periods of time. Moreover, rain is patchy in nature - at any instant on the average only a certain fraction of the observed pixels contain rain. The fraction of area covered by rain decreases, as the size of a pixel becomes smaller. This means that within what looks like a patch of rainy area in a coarse resolution view with larger pixel size, one finds clusters of rainy and dry patches when viewed on a finer scale. The model makes definite predictions about how these and other related statistics depend on the pixel size. These predictions were found to agree well with data. In a subsequent second part of the work we plan to test the model with rain gauge data collected during the TRMM (Tropical Rainfall Measuring Mission) ground validation campaign.

Kundu, Prasun K.↗