SEARCH · Engineering Papers

Results for “statistics of extremes”

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.

SourceResearch connectionFrom yearThrough year

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

At least 433 records · Page 24

All sources NASA NTRS DOE OSTI

NASA NTRS2021

A Quarter Century of Wind Spacecraft Discoveries

The Wind spacecraft, launched on November 1, 1994, is a critical element in NASA’s Heliophysics System Observatory (HSO)—a fleet of spacecraft created to understand the dynamics of the Sun-Earth system. The combination of its longevity (>25 years in service), its diverse complement of instrumentation, and high resolution and accurate measurements has led to it becoming the “standard candle” of solar wind measurements. Wind has over 55 selectable public data products with over ∼1,100 total data variables (including OMNI data products) on SPDF/CDAWeb alone. These data have led to paradigm shifting results in studies of statistical solar wind trends, magnetic reconnection, large-scale solar wind structures, kinetic physics, electromagnetic turbulence, the Van Allen radiation belts, coronal mass ejection topology, interplanetary and interstellar dust, the lunar wake, solar radio bursts, solar energetic particles, and extreme astrophysical phenomena such as gamma-ray bursts. This review introduces the mission and instrument suites then discusses examples of the contributions by Wind to these scientific topics that emphasize its importance to both the fields of heliophysics and astrophysics.

Lynn B. Wilson III↗

NASA NTRS1974

Electron capture in carbon dwarf supernovae

The rates of electron capture on heavier elements under the extreme conditions predicted for dwarf star supernovae have been computed, incorporating modifications that seem to be indicated by present experimental results. An estimate of the maximum possible value of such rates is also given. The distribution of nuclei in nuclear statistical equilibrium has been calculated for the range of expected supernovae conditions, including the effects of the temperature dependence of nuclear partition functions. These nuclide abundance distributions are then used to compute nuclear equilibrium thermodynamic properties. The effects of the electron capture on such equilibrium matter are discussed. In the context of the 'carbon detonation' supernova model, the dwarf central density required to ensure core collapse to a neutron star configuration is found to be slightly higher than that obtained by Bruenn (1972) with the electron capture rates of Hansen (1966).-

Mazurek, T. J.↗

NASA NTRS1964

Unmanned Spacecraft of the United States

In 1957 the first earth satellite ushered in the age of space flight. Since that historic event, space exploration has become a major national objective of both the United States and the Soviet Union. These two nations have attempted a total of well over 200 space flight missions. Other nations are also participating in various degrees in what will continue to grow as a cooperative world effort. In the years since 1957, man has successfully flown in earth orbit. He has initiated programs to land on the moon and return. He has made dramatic applications of earth satellites in meteorology, communications, navigation, and geodesy. A host of scientific satellites.continue to advance understanding of the earth's environment, the sun, and the stars. Automated spacecraft are being flown to the moon, deep into interplanetary space, and to the near planets, Mars and Venus. One of the most exciting technological aspects of space exploration has been the development of automated spacecraft. Most of the scientific exploration of space and the useful applications of space flight thus far have been made possible by automated spacecraft. Development of these spacecraft and their many complex subsystems is setting the pace today for many branches of science and technology. Guidance, computer, attitude control, power, telecommunication, instrumentation, and structural subsystems are being subjected to new standards of light weight, high efficiency, extreme accuracy, and unsurpassed reliability and quality. This publication reviews the automated spacecraft which have been developed and flown, or which are under active development in the United States by the National Aeronautics and Space Administration. From the facts and statistics contained herein, certain observations can be made and certain conclusions drawn.

Cortright, Edgar M.↗

NASA NTRS2010

Advances in Landslide Hazard Forecasting: Evaluation of Global and Regional Modeling Approach

A prototype global satellite-based landslide hazard algorithm has been developed to identify areas that exhibit a high potential for landslide activity by combining a calculation of landslide susceptibility with satellite-derived rainfall estimates. A recent evaluation of this algorithm framework found that while this tool represents an important first step in larger-scale landslide forecasting efforts, it requires several modifications before it can be fully realized as an operational tool. The evaluation finds that the landslide forecasting may be more feasible at a regional scale. This study draws upon a prior work's recommendations to develop a new approach for considering landslide susceptibility and forecasting at the regional scale. This case study uses a database of landslides triggered by Hurricane Mitch in 1998 over four countries in Central America: Guatemala, Honduras, EI Salvador and Nicaragua. A regional susceptibility map is calculated from satellite and surface datasets using a statistical methodology. The susceptibility map is tested with a regional rainfall intensity-duration triggering relationship and results are compared to global algorithm framework for the Hurricane Mitch event. The statistical results suggest that this regional investigation provides one plausible way to approach some of the data and resolution issues identified in the global assessment, providing more realistic landslide forecasts for this case study. Evaluation of landslide hazards for this extreme event helps to identify several potential improvements of the algorithm framework, but also highlights several remaining challenges for the algorithm assessment, transferability and performance accuracy. Evaluation challenges include representation errors from comparing susceptibility maps of different spatial resolutions, biases in event-based landslide inventory data, and limited nonlandslide event data for more comprehensive evaluation. Additional factors that may improve algorithm performance accuracy include incorporating additional triggering factors such as tectonic activity, anthropogenic impacts and soil moisture into the algorithm calculation. Despite these limitations, the methodology presented in this regional evaluation is both straightforward to calculate and easy to interpret, making results transferable between regions and allowing findings to be placed within an inter-comparison framework. The regional algorithm scenario represents an important step in advancing regional and global-scale landslide hazard assessment and forecasting.

Kirschbaum, Dalia B.↗

NASA NTRS2003

Vector-Ordering Filter Procedure for Data Reduction

The vector-ordering filter (VOF) technique involves a procedure for sampling a large population of data vectors to select a subset of data vectors that fully characterize the state space of the large population. The VOF technique enables a large reduction of the volume of data that must be handled in the automated monitoring system and method discussed in the two immediately preceding articles. In so doing, the VOF technique enables the development of data-driven mathematical models of a monitored asset from sets of data that would otherwise exceed the memory capacities of conventional engineering computers. Data-driven mathematical models have been shown to offer high fidelity for purposes of control and monitoring of assets. In practice, a collection of asset-operating observations is acquired with the intention that the collection contain observations characteristic of the full dynamic range of operation of the asset. Often, such a collection contains an extremely large number of observations, many of which are redundant. The VOF technique fills the need for a means to extract, from the original collection of observational data, a reduced data matrix that excludes redundant data while maintaining the full statistical character and dynamic range of the original data. The reduced data matrix can then be used as the input data for development of a mathematical model of the monitored asset, or as training data for a neural-network substitute for an explicit mathematical model of the asset. Alternatively, the reduced data matrix can, itself, be used directly as a mathematical model of the monitored asset, as is commonly done in multivariate state-estimation techniques. The original data are collected from the asset over a range of operating states and are put in matrix form. Each column vector in the original data matrix represents the signal values acquired at a particular operational state of the asset. Thus, the number of columns of the original data matrix equals the number of observed states and the number of rows in this matrix equals the number of signals acquired at each observation. In the VOF technique, one extracts the reduced data matrix from the original data matrix through the selection of a representative subset of the column (state) vectors.

Bickford, Randall L.↗

NASA NTRS2022

A Dynamic Landslide Hazard Monitoring Framework for the Lower Mekong Region

The Lower Mekong region is one of the most landslide-prone areas of the world. Despite the need for dynamic characterization of landslide hazard zones within the region, it is largely understudied for several reasons. Dynamic and integrated understanding of landslide processes requires landslide inventories across the region, which have not been available previously. Computational limitations also hamper regional landslide hazard assessment, including accessing and processing remotely sensed information. Finally, open-source software and modelling packages are required to address regional landslide hazard analysis. Leveraging an open-source data-driven global Landslide Hazard Assessment for Situational Awareness model framework, this study develops a region-specific dynamic landslide hazard system leveraging satellite-based Earth observation data to assess landslide hazards across the lower Mekong region. A set of landslide inventories were prepared from high-resolution optical imagery using advanced image-processing techniques. Several static and dynamic explanatory variables (i.e., rainfall, soil moisture, slope, relief, distance to roads, distance to faults, distance to rivers) were considered during the model development phase. An extreme gradient boosting decision tree model was trained for the monsoon period of 2015–2019 and the model was evaluated with independent inventory information for the 2020 monsoon period. The model performance demonstrated considerable skill using receiver operating characteristic curve statistics, with Area Under the Curve values exceeding 0.95. The model architecture was designed to use near-real-time data, and it can be implemented in a cloud computing environment (i.e., Google Cloud Platform) for the routine assessment of landslide hazards in the Lower Mekong region. This work was developed in collaboration with scientists at the Asian Disaster Preparedness Center as part of the NASA SERVIR Program’s Mekong hub. The goal of this work is to develop a suite of tools and services on accessible open-source platforms that support and enable stakeholder communities to better assess landslide hazard and exposure at local to regional scales for decision making and planning.

Nishan Kumar Biswas↗

DOE OSTI2021

Open-Loop Control of Adjustable Tuned Mass Dampers for Floating Wind Turbine Platforms: Preprint

Floating offshore wind turbines have several advantages over their land-based counterparts, including access to stronger, more consistent winds and proximity to population centers. However, increased costs from the platform have been a challenge inhibiting their widespread adoption. New damping technologies with tunable target frequencies and damping ratios promise a greater degree of control over platform motions, allowing hulls to be designed smaller and reducing the overall cost of the turbines. In this work, a control strategy is proposed for these adjustable tuned mass dampers (TMDs). A frequency-domain model of the system is developed, from which response amplitude operators (RAOs) for the platform's rigid-body modes can be obtained for different natural frequencies and damping ratios of the dampers. Using these RAOs, and JONSWAP spectrums of operational and extreme sea states, the performance of the various damper settings are compared by evaluating the standard deviation of a rigid-body motion of interest (platform pitch or heave). By finding the optimal damper setting for a range of representative design load cases (DLCs), a lookup function is made to return the damper target frequency given the peak period of the sea state. To determine the current sea state, we propose a simple estimator that uses wave measurements (e.g., from a measurement buoy) to determine the peak period and significant wave height. The significant wave height is determined using the statistics of the past 100 seconds of wave elevation information and the peak period is computed using a frequency locked loop. Initial results show that we can use these estimated values to control the TMD natural frequency with an open loop controller. The controlled system is tested in a range of environmental conditions determined by the International Electrotechnical Commission design load cases (DLCs), which include normal and extreme wind and wave models. In these tests, we compare the performance of 4 cases: (1) no TMD, (2) a constant TMD based on the worst case DLC, (3) a controlled TMD based on known wind and wave environments, and (4) a real-time controlled TMD using estimated wind and wave environments. The effect of hull-based TMD control is also compared to changes in traditional wind turbine control via blade pitch.

control↗

DOE OSTI2025

Temperature measurement of Quark-Gluon plasma at different stages

In a Quark-Gluon Plasma (QGP), the fundamental building blocks of matter, quarks and gluons, are under extreme conditions of temperature and density. A QGP could exist in the early stages of the Universe, and in various objects and events in the cosmos. The thermodynamic and hydrodynamic properties of the QGP are described by Quantum Chromodynamics (QCD) and can be studied in heavy-ion collisions. Despite being a key thermodynamic parameter, the QGP temperature is still poorly known. Thermal lepton pairs (e + e − and μ + μ − ) are ideal penetrating probes of the true temperature of the emitting source, since their invariant-mass spectra suffer neither from strong final-state interactions nor from blue-shift effects due to rapid expansion. Here we measure the QGP temperature using thermal e+e− production at the Relativistic Heavy Ion Collider (RHIC). The average temperature from the low-mass region (in-medium ρ 0 vector-meson dominant) is (2.01 ± 0.23) × 10 12 K, consistent with the chemical freeze-out temperature from statistical models and the phase transition temperature from Lattice QCD. The average temperature from the intermediate mass region (above the ρ0 mass, QGP dominant) is significantly higher at (3.25 ± 0.60) × 10 12 K. This work provides essential experimental thermodynamic measurements to map out the QCD phase diagram and understand the properties of matter under extreme conditions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Spatiotemporal Learning in Power Modules: Wavelet-Enhanced Forecasting of Thermomechanical Degradation

Detecting internal defects in power electronics packages is critical for their performance and reliability, especially under extreme operating conditions, as these defects can lead to catastrophic failure if not properly addressed. Confocal scanning acoustic microscopy (C-SAM) plays a key role in the nondestructive evaluation of bond layer degradation within a power electronics package by detecting defects such as delamination, voids, and cracks. However, accurately quantifying and predicting these defects from C-SAM images remains a significant challenge due to the low noise-to-signal ratio, which typically arises from both imaging process and bond patterns itself. In this paper, we explore machine learning strategies for processing C-SAM images and providing predictive models of defect growth. We use C-SAM images of sintered copper and sintered silver samples, which are obtained under accelerated thermal experiments, as the representative dataset for our study. We investigate the effect of Fourier transforms and wavelet transforms on these datasets to remove high-frequency noise and address noise across multiple scales with histogram equalization to enhance the contrast and improve the visibility of defects. As a result, defect boundaries can be clearly distinguished, enabling more accurate tracking of their growth over time. We then employ different time-series forecasting algorithms on the denoised images to formulate an image-based lifetime prediction model. Statistical models and deep-learning techniques are trained on images obtained in the early stages of thermal shock, and defect growth in the later stages is predicted. Our work serves as a preliminary attempt to improve the accuracy of lifetime prediction models of power electronics packages, which is critical under extreme operating environments.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2020

Exploring topography-based methods for downscaling subgrid precipitation for use in Earth System Models

Topography exerts major control on land surface processes. To improve representation of topographic impacts on land surface processes, a new topography-based subgrid structure has been introduced to the Energy Exascale Earth System Model to represent the subgrid heterogeneity of surface elevation. To take advantage of the new subgrid structure for improving land surface modeling, this study explores four variations of topography-based methods for downscaling grid precipitation to the corresponding subgrids. In the first three methods, the deviation of the subgrid precipitation from the grid’s is equal to the grid precipitation multiplied by the ratio of the elevation difference between the subgrid and grid mean to a specified elevation equals to the grid elevation, the difference between the maximum and minimum subgrid elevation, and the maximum subgrid elevation, respectively. The second method limits the ratio to 0.5 to avoid extreme values on mountains and the third method accounts for the slope effect. The fourth method is similar to the third method except that the Froude Number is used to limit the numerator in a blocking regime of the ambient flow. The downscaled precipitation is evaluated using the PRISM precipitation data over the U.S. using statistical metrics. Results show that by accounting for topographic slope besides elevation, the third and fourth methods show clear advantages over the first and second methods. Furthermore, introducing the Froude Number in the fourth method improves downscaling skill and shows consistent advantages over the third method in areas with larger subgrid heterogeneity across different grids sizes.

Tesfa, Teklu K.↗

NASA NTRS1993

Percolation technique for galaxy clustering

We study percolation in mass and galaxy distributions obtained in 3D simulations of the CDM, C + HDM, and the power law (n = -1) models in the Omega = 1 universe. Percolation statistics is used here as a quantitative measure of the degree to which a mass or galaxy distribution is of a filamentary or cellular type. The very fast code used calculates the statistics of clusters along with the direct detection of percolation. We found that the two parameters mu(infinity), characterizing the size of the largest cluster, and mu-squared, characterizing the weighted mean size of all clusters excluding the largest one, are extremely useful for evaluating the percolation threshold. An advantage of using these parameters is their low sensitivity to boundary effects. We show that both the CDM and the C + HDM models are extremely filamentary both in mass and galaxy distribution. The percolation thresholds for the mass distributions are determined.

Klypin, Anatoly↗

NASA NTRS1997

Subacute diabetic proximal neuropathy

OBJECTIVE: To evaluate the clinical, electrophysiologic, autonomic, and neuropathologic characteristics and the natural history of subacute diabetic proximal neuropathy and its response to immunotherapy. MATERIAL AND METHODS: For the 12-year period from 1983 to 1995, we conducted a retrospective review of medical records of Mayo Clinic patients with diabetes who had subacute onset and progression of proximal weakness. The responses of treated versus untreated patients were compared statistically. RESULTS: During the designated study period, 44 patients with subacute diabetic proximal neuropathy were encountered. Most patients were middle-aged or elderly, and no sex preponderance was noted. The proximal muscle weakness often was associated with reduced or absent lower extremity reflexes. Associated weight loss was a common finding. Frequently, patients had some evidence of demyelination on nerve conduction studies, but it invariably was accompanied by concomitant axonal degeneration. The cerebrospinal fluid protein concentration was usually increased. Diffuse and substantial autonomic failure was generally present. In most cases, a sural nerve biopsy specimen suggested demyelination, although evidence of an inflammatory infiltrate was less common. Of 12 patients who received treatment (with prednisone, intravenous immune globulin, or plasma exchange), 9 had improvement of their conditions, but 17 of 29 untreated patients (59%) with follow-up also eventually had improvement, albeit at a much slower rate. Improvement was usually incomplete. CONCLUSION: We suggest that the entity of subacute diabetic proximal neuropathy is an extensive and severe variant of bilateral lumbosacral radiculoplexopathy, with some features suggestive of an immune-mediated cause. It differs from chronic inflammatory demyelinating polyradiculoneuropathy in that most cases have a more restricted distribution and seem to be monophasic and self-limiting. The efficacy of immunotherapy is unproved, but such intervention may be considered in the severe and progressive cases or ones associated with severe neuropathic pain.

Non-NASA Center↗

NASA NTRS2013

Prediction of Seasonal Climate-induced Variations in Global Food Production

Consumers, including the poor in many countries, are increasingly dependent on food imports and are therefore exposed to variations in yields, production, and export prices in the major food-producing regions of the world. National governments and commercial entities are paying increased attention to the cropping forecasts of major food-exporting countries as well as to their own domestic food production. Given the increased volatility of food markets and the rising incidence of climatic extremes affecting food production, food price spikes may increase in prevalence in future years. Here we present a global assessment of the reliability of crop failure hindcasts for major crops at two lead times derived by linking ensemble seasonal climatic forecasts with statistical crop models. We assessed the reliability of hindcasts (i.e., retrospective forecasts for the past) of crop yield loss relative to the previous year for two lead times. Pre-season yield predictions employ climatic forecasts and have lead times of approximately 3 to 5 months for providing information regarding variations in yields for the coming cropping season. Within-season yield predictions use climatic forecasts with lead times of 1 to 3 months. Pre-season predictions can be of value to national governments and commercial concerns, complemented by subsequent updates from within-season predictions. The latter incorporate information on the most recent climatic data for the upcoming period of reproductive growth. In addition to such predictions, hindcasts using observations from satellites were performed to demonstrate the upper limit of the reliability of crop forecasting.

variations↗

NASA NTRS2015

Measuring Skew in Average Surface Roughness as a Function of Surface Preparation

Characterizing surface roughness is important for predicting optical performance. Better measurement of surface roughness reduces polishing time, saves money and allows the science requirements to be better defined. This study characterized statistics of average surface roughness as a function of polishing time. Average surface roughness was measured at 81 locations using a Zygo white light interferometer at regular intervals during the polishing process. Each data set was fit to a normal and Largest Extreme Value (LEV) distribution; then tested for goodness of fit. We show that the skew in the average data changes as a function of polishing time.

Stahl, Mark↗

DOE OSTI2023

An upper bound for extreme temperatures over midlatitude land

Heatwaves damage societies worldwide and are intensifying with global warming. Several mechanistic drivers of heatwaves, such as atmospheric blocking and soil moisture-atmosphere feedback, are well-known for their ability to raise surface air temperature. However, what limits the maximum surface air temperature in heatwaves remains unclear; this became evident during recent Northern Hemisphere heatwaves which achieved temperatures far beyond the upper tail of the observed statistical distribution. Here, we present evidence for the hypothesis that convective instability limits annual maximum surface air temperatures (TXx) over midlatitude land. We provide a theory for the corresponding upper bound of midlatitude temperatures, which accurately describes the observed relationship between temperatures at the surface and in the midtroposphere. We show that known heatwave drivers shift the position of the atmospheric state in the phase space described by the theory, changing its proximity to the upper bound. This theory suggests that the upper bound for midlatitude TXx should increase 1.9 times as fast as 500-hPa temperatures at the time and location of TXx occurrences. Using empirical 500-hPa warming, we project that the upper bound of TXx over Northern Hemisphere midlatitude land (40°N to 65°N) will increase about twice as fast as global mean surface air temperature, and TXx will increase faster than this bound over regions that dry on the hottest days.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Self-Organizing Map-Based Resilience Quantification and Resilient Control of Distribution Systems Under Extreme Events

Due to climate change, extreme weather events are occurring more frequently and with increasing impact. This trend poses a significant challenge for distribution system operators (DSO) to ensure that there is uninterrupted power supply to critical loads in their networks. To embed resilience into DSO's decision-making, resilience needs to be first quantified and then integrated into the system-level optimization. Therefore, this paper first develops a novel self-organizing map (SOM) based method (called SomRes) to quantify the time-varying resilience index of a system that can leverage the powerful classification property of SOMs and removes some of the disadvantages of subjective weight assignment methods. Using SomRes, a resilient resource allocation and operational dispatch algorithm is further developed to enhance system resilience against extreme events by considering the SomRes resilience index directly as the feedback. Here, the proposed resilience quantification approach is benchmarked with a state-of-the-art approach and the efficacy of the proposed resilient dispatch algorithm is demonstrated through several deterministic and statistical case studies on the IEEE 123-bus distribution system. Simulation studies show that the proposed SomRes quantification method is an appropriate indicator of system resilience, and the resilient resource allocation and dispatch strategy can significantly reduce critical load shedding under varying event propagation scenarios.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2026

Optimizing enzymes for plastic upcycling using machine learning design and high throughput experiments

Plastic use is ubiquitous in the modern world, and polyethylene terephthalate (PET) is one of the most abundantly produced plastics (and the most highly produced polyester), with ~65 million metric tons manufactured annually. To the consumer, PET is likely most recognizable as the plastic used to make beverage bottles. Like many plastics, traditional mechanical or chemical means of PET deconstruction and upcycling are costly and inefficient. Because of these challenges, recycled plastic is generally of lower quality and is more expensive to produce than virgin plastic derived from petroleum. Ultimately, this results in most plastic ending up as waste. We view plastic waste as an underutilized resource which, with the development of more efficient and high-quality recycling processes, could (1) generate significant economic value while (2) decreasing petroleum usage and greenhouse gas emissions, as well as (3) minimizing its negative environmental and health impacts. Biocatalytic recycling, or biomanufacturing the basic building blocks of new plastic from plastic waste, is a promising approach to plastic reuse that complements existing recycling technologies. Recently, biological enzymes capable of breaking down PET have garnered significant attention as an attractive means of dealing with the plastic problem. These enzymes are currently undergoing pilot studies for implementation in industrial-scale enzyme-based recycling. However, there are significant limitations to current enzymes, including the need to perform costly pre-processing of the plastic waste before the enzymes are able to work. Further optimization of these enzymes is necessary to make these technologies competitive, and ultimately incentivise industry-wide adoption of this biology-based green recycling technology. n this work we demonstrate a means to design and generate performant biological enzymes, capable of efficiently deconstructing plastic waste. Specifically, we applied recent advances in artificial intelligence, machine learning, and statistical analysis to design new versions and discover natural enzymes capable of breaking down PET. We focused on optimizing key properties that are important for industrial-scale enzymatic recycling such as pH and thermotolerance. Normal testing of enzymatic plastic-deconstruction is extremely labor intensive and so through this work we also developed a robotic-assisted experimental pipeline capable of characterizing thousands of candidate enzymes. The results of this iterative, AI-guided, multi-discipline approach have led to increases in enzymatic breakdown of over 150X over starting enzymes. This work supports the rapidly developing and transformative field of biocatalytic solutions to environmental problems beyond the discovery and predictive understanding of enzymes for polymer recycling, and has wide implications for tackling numerous energy problems such as carbon capture and fixation (e.g., engineering carbon monoxide dehydrogenase and the rubisco-pathway), biomining (e.g., design of lanthanide-binding proteins) and biomanufacturing (e.g., lignin-deconstruction enzymes).

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

NASA NTRS2003

Stochastic Models of Tropical Rain-Rate Fields

Because of the extreme variability of rain rate in space and time and the difficulties with remote sensing methods of measuring rain rates, accurate determination of rainfall over large areas and time periods has long been a problem for hydrologists, meteorogists, and climatologists. A number of statistical models of rain have been developed in order to investigate the impact of rain variability on satellite remote sensing methods, validation of satellite rain products, and generation of rain maps with accompanying error estimates. These models may be useful in examining 'sub-grid scale' issues in representing precipitation in numerical mdoels. A stochastic model will first be described which can generate time-dependent high-resolution spatial rain fields with space and time correlations similar to those seen in rain data, as well as representing the presence of areas with zero rain rate and log-normally distributed rain rates where there is rain. A simpler model derived from this, formulated in the spectral domain, seems to imply fractal-like rain statistics at small scales when fit to rain data.

Bell, Thomas L.↗