Search NASA⌕ Search

SEARCH · Search NASA

Results for “temperature algorithm”

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 1,387 records · Page 77

The Influence of Radiosonde 'Age' on TRMM Field Campaign Soundings Humidity Correction

Hundreds of Vaisala sondes with a RS80-H Humicap thin-film capacitor humidity sensor were launched during the Tropical Rainfall Measuring Mission (TRMM) field campaigns in Large Scale Biosphere-Atmosphere held in Brazil (LBA) and in Kwajalein experiment (KWAJEX) held in the Republic of Marshall Islands. Using Six humidity error correction algorithms by Wang et al., these sondes were corrected for significant dry bias in the RS80-H data. It is further shown that sonde surface temperature error must be corrected for a better representation of the relative humidity. This error becomes prominent due to sensor arm-heating in the first 50-s data.

Roy, Biswadev↗

First Atmospheric Science Results from the Mars Exploration Rovers Mini-TES

Thermal infrared spectra of the martian atmosphere taken by the Miniature Thermal Emission Spectrometer (Mini-TES) were used to determine the atmospheric temperatures in the planetary boundary layer and the column-integrated optical depth of aerosols. Mini-TES observations show the diurnal variation of the martian boundary layer thermal structure, including a near-surface superadiabatic layer during the afternoon and an inversion layer at night. Upward-looking Mini-TES observations show warm and cool parcels of air moving through the Mini-TES field of view on a time scale of 30 seconds. The retrieved dust optical depth shows a downward trend at both sites.

Mars↗

NASA Tech Briefs, August 2008

Customizable Digital Receivers for Radar Two-Camera Acquisition and Tracking of a Flying Target Visual Data Analysis for Satellites A Data Type for Efficient Representation of Other Data Types Hand-Held Ultrasonic Instrument for Reading Matrix Symbols Broadband Microstrip-to-Coplanar Strip Double-Y Balun A Topographical Lidar System for Terrain-Relative Navigation Programmable Low-Voltage Circuit Breaker and Tester Electronic Switch Arrays for Managing Microbattery Arrays Topics covered include: Lower-Dark-Current, Higher-Blue-Response CMOS Imagers; Fabricating Large-Area Sheets of Single-Layer Graphene by CVD; Support for Diagnosis of Custom Computer Hardware; Providing Goal-Based Autonomy for Commanding a Spacecraft; Dynamic Method for Identifying Collected Sample Mass; Optimal Planning and Problem-Solving; Attitude-Control Algorithm for Minimizing Maneuver Execution Errors; Grants Document-Generation System; Heat-Storage Modules Containing LiNO3 3H2O and Graphite Foam; Precipitation-Strengthened, High-Temperature, High-Force Shape Memory Alloys; Improved Relief Valve Would Be Less Susceptible to Failure; Safety Modification of Cam-and-Groove Hose Coupling; Using Composite Materials in a Cryogenic Pump; Using Electronic Noses to Detect Tumors During Neurosurgery; Producing Newborn Synchronous Mammalian Cells; Smaller, Lower-Power Fast-Neutron Scintillation Detectors; Rotationally Vibrating Electric-Field Mill; Estimating Hardness from the USDC Tool-Bit Temperature Rise; Particle-Charge Spectrometer; Automated Production of Movies on a Cluster of Computers; FIDO-Class Development Rover; and Tone-Based Command of Deep Space Probes Using Ground Antennas.

Source record↗

Radiative Transfer and Retrievals in EOF Domain

The Infrared Atmospheric Sounding Interferometer (IASI) is a hyperspectral sensor with 8461 spectral channels and a nominal spectral resolution of 0.25 cm(sup -1). It is computationally intensive to perform radiative transfer calculations and inversions using all these channels. We will present a Principal Component-based Radiative Transfer Model (PCRTM) and a retrieval algorithm which perform all the necessary calculations in EOF domain. Since the EOFs are orthogonal to each other, only about 100 principal components are needed to represent the information content of the 8461 channels. The PCRTM provides the EOF coefficients and associated derivatives with respect to atmospheric and surface parameters needed by the inversion algorithm. The inversion algorithm is based on a non-linear Levenberg-Marquardt method with climatology covariance and a priori information as constraints. The retrieved parameters include atmospheric temperature, moisture and ozone profiles, cloud parameters, surface skin temperature, and surface emissivities. To make the retrieval system even more compact and stable. The atmospheric vertical profiles are compressed into the EOF space as well. The surface emissivities are also compressed into EOF space.

Liu, Xu↗

Performance of the GLAS Space Lidar Receiver Through Its Seven-Year Space Mission

NASA s Ice, Cloud, and land Elevation Satellite (ICESat) mission [1,2] carrying the Geoscience Laser Altimeter System (GLAS) Instrument, was launched on January 12, 2003. The three lasers on ICESat have made a total of 1.98 billion laser shot measurements of the Earth s surface and atmosphere during its 17 science data collection campaigns over its seven year operating lifetime. ICESat completed its science mission after the last laser stopped operating in October 2009. The spacecraft was de-orbited on August 30, 2010. The GLAS instrument carried 3 diode-pumped Q-switched Nd:YAG lasers, which emitted 6-nsec wide pulses at 1064 and 532 nm at a 40-Hz rate. There are three lidar receiver channels, a 1064 nm surface altimetry channel, a 1064 nm cloud backscattering lidar channel, and a 532 nm cloud and aerosol backscattering lidar channel. The altimetry and cloud backscatter channels used Si avalanche photodiode (APD) operated in analog mode as in the Mars Global Surveyor s Mars Orbital Laser Altimeter [3,4]. GLAS also utilized a number of new technologies and techniques for space lidar, including passively Q-switched diode-pumped Nd:YAG lasers, a 1-m diameter telescope, a temperature tuned etalon optical bandpass filter, Si APD single photon counting detectors, 1 Gsample/sec waveform digitizers, ultra stable clock oscillators, and digital signal processing and detection algorithms [5]. A global position system (GPS) receiver was used to provide the spacecraft position and epoch times. The ICESat mission provided a unique opportunity to monitor the lidar component performance in the space environment over a multi-year time period. We performed a number of engineering tests periodically to monitor the lidar receiver performance, including receiver sensitivity, timing precision, detector dark noise, etc. A series of engineering tests were also performed after the end of the science mission to evaluate the performance of the spare detector, oscillator, waveform digitizer, and GPS receiver. An experiment was conducted which pointed GLAS to Venus to test the receiver sensitivity to star light and to verify GLAS bore sight with respect to the spacecraft coordinate system. These tests provided unique data to assess the degradation and the rate of change of these key lidar components due to space radiation and aging. They also helped to validate new techniques to operate and calibrate future space lidars.

Sun, Xiaoli↗

Reflectivity, Emissivity, and Brightness from Multilayered Soil with Linearly Varying Permittivity at P- and L- bands

Passive sensing of root zone soil moisture is theoretically investigated using P and L band frequencies. The soil is considered to be a stratified medium. The reflectivity from the soil surface is found by considering the soil to be composed of a finely layered structure with layer widths that are small compared to a wavelength. A well-known transmission line algorithm is then used to compute the reflectivity at the surface. Emissivity is obtained as (1 – reflectivity). It assumed that the soil is at a constant temperature of 300o K when determining the brightness temperature. A soil moisture profile consisting of two homogeneous layers separated by a layer of linearly varying soil moisture is investigated. The study is made at an incident angle of 40O for horizontal polarization at 0.8 and 1.4 GHz. Two variations of the basic profile are considered: a transition from a layer of lesser moisture to a layer of greater moisture and the opposite. In both cases the P band signal shows greater penetration than the L band case. There is greater penetration in the transition to drier soils. In both cases the L band signals show sensitivity to only the upper surface of the soil as expected.

Ming Li↗

Investigating seasonal ENSO forecast amplitude calibration for GEOS-S2S-2

The GEOS-S2S-2 is a global coupled model and assimilation system, encompassing many aspects of the Earth climate system. This project looked at seasonal (nine-month) ensemble forecasts, over the period from 1982 through the present. We explored and validated a forecast amplitude correction technique used by the North American Multi-Model Ensemble (NMME), of which GEOS is a member, for Niño3.4 sea-surface temperature (SST) anomaly predictions. The method relies on deriving a set of correction factors based on the standard deviations of hindcast and observed SST anomalies. This algorithm was implemented at the Global Modeling and Assimilation Office (GMAO) and applied to the GEOS-S2S-2 Niño 3.4 hindcasts. In general, the correction reduced the forecast amplitude. This improved the quality of the forecasts when measured by the root mean square error (RMSE) of the ensemble mean forecast. RMSE decreased for most initialization months and forecast leads. The correction had a larger impact on seasons with strong El Niño-Southern Oscillation (ENSO) events than neutral seasons. The effects of the correction on the ensemble characteristics were also examined, and a variation of the algorithm was proposed that preserves the ensemble spread. The outcome of this project is a tool that can be used to calibrate GEOS-S2S-2 ENSO forecasts and potentially improve the precision of multi-model El Niño outlooks.

ENSO↗

Seasonal Surface Spectral Emissivity Derived from VIIRS Data

Surface emissivity is essential for many remote-sensing applications including the retrieval of surface skin temperature from satellite-based infrared measurements, the determination of cloud detection thresholds, and the estimation of the surface longwave radiation emission, an important component of the energy budget of the surface-atmosphere interface. The CERES (Clouds and the Earth’s Radiant Energy System) Project is measuring broadband shortwave and longwave radiances and deriving cloud properties from the MODIS on Terra and Aqua and from the VIIRS on NOAA-19 and NOAA-20 orbiters to produce combined global radiation and cloud property data sets. Zhou et al. (IEEE Trans. Geosci. Remote Sens., 49, 2011) used Infrared Atmospheric Sounding Interferometer (IASI) data to create a high spectral resolution surface emissivity atlas for remote sensing and modeling applications. The IASI measures spectral radiances between 3.62 and 15.5 µm. The VIIRS I4 channel width is from 3.55 to 3.93 µm, while MODIS Band 20 is from 3.66 to 3.84 µm. Comparisons of top-of-atmosphere (TOA) radiance calculations with MODIS and VIIRS observations for these bands suggest that the IASI emissivity atlas near 3.7 µm may not be suitable for CERES cloud retrievals. In this paper, the IASI emissivities for the VIIRS and MODIS bands centered near 11µm are used to derive surface skin temperature from nighttime MODIS/VIIRS data. The Goddard Earth Observing System for Instrument Teams (GEOS-IT) numerical weather analyses provide temperature and water vapor profiles fused to correct the observed radiances for atmospheric absorption and emission. Global seasonal emissivity maps are then derived for the VIIRS and MODIS 3.7µm bands that are consistent with the derived skin temperatures and the observed TOA radiances. These seasonal climatology maps will be validated and used in CERES Edition 5 and other CERES-related cloud retrieval algorithms to provide improved clear-sky radiances and derived cloud properties.

Surface Emissivity↗

Thermal Actuator Identification and Control for Thermomechanical Real-Time Cyber–Physical Testing

Thermomechanical cyber–physical testing enables two-way thermal coupling between a numerical and an experimental subsystem. The interactions between the numerical model and the physical specimen occur through transfer systems, which enforce interface conditions. Thus, efficient control methodologies are necessary to achieve the desired interface interaction through thermal actuators with minimal error. This study introduces a novel thermal transfer system that imposes distributed cooling (or heating) thermal loads on a physical subsystem. First, the thermal actuator is identified considering switching-mode continuous dynamics for heating and cooling conditions. A switching-mode estimation algorithm is adopted to estimate the operating thermal cycle of the actuator in real-time. A control system is developed to experimentally impose the desired temperature and reduce tracking error (i.e., the error between the desired and actual temperature) under different thermal cycles. The identification and control of the thermal transfer system are then validated through a set of experiments considering different temperature rates of change. The developed control system is found to effectively minimize tracking errors in real-time cyber–physical experiments.

Herta Montoya↗

Efficient GPU-Accelerated MultiSource Global Fit Pipeline for LISA Data Analysis

The large-scale analysis task of deciphering gravitational-wave signals in the LISA data stream will be difficult, requiring a large amount of computational resources and extensive development of computational methods. Its high dimensionality, multiple model types, and complicated noise profile require a global fit to all parameters and input models simultaneously. In this work, we detail our global fit algorithm, called “Erebor,” designed to accomplish this challenging task. It is capable of analyzing current state-of-the-art datasets and then growing into the future as more pieces of the pipeline are completed and added. We describe our pipeline strategy, the algorithmic setup, and the results from our analysis of the LDC2A Sangria dataset, which contains massive black hole binaries, compact galactic binaries, and a parametrized noise spectrum whose parameters are unknown to the user. The Erebor algorithm includes three unique and very useful contributions: GPU acceleration for enhanced computational efficiency; ensemble Markov Chain Monte Carlo (MCMC) sampling with multiple MCMC walkers per temperature for better mixing and parallelized sample creation; and special online updates to reversible-jump (or transdimensional) sampling distributions to ensure sampler mixing and accurate initial estimates for detectable sources in the data.We recover posterior distributions for all 15 (6) of the injected massive black hole binaries (MBHB) in the LDC2A training (hidden) dataset. We catalog ∼12000 galactic binaries (∼8000 as high confidence detections) for both the training and hidden datasets. All of the sources and their posterior distributions are provided in publicly available catalogs.

LISA↗

Initial exploration of a novel transient arrest system involving fuel heating

A preliminary analysis on a novel accident response system to diminish the severity of super- critical transients was conducted. The novel accident response system, called the instant shock arrest system, involves using electricity to heat the nuclear fuel at the onset of a large accidental reactivity insertion. This system is specifically designed for reactors with metallic fuel, such that the fuel is capable of conducting electricity, and being resistively heated. A reactor dynamics model of the advanced test reactor was created using the point kinetics equations and a linear reactivity feedback model to simulate how the system would effect the maximum fuel temperatures experienced during the transient. Transients with the instant shock arrest system were compared to those without it. It was found that the instant shock arrest system initially heated the fuel more than the unaffected transient but the negative reactivity inserted from such heating was enough to lower the maximum fuel temperature experienced during the transient. After simulating six different accident scenarios with reactivity insertions ranging from 0.5 to 1.3 dollar, it was found that an optimal system response could reduce peak fuel temperatures during the transient by 3.5% to 5%. Furthermore, discussion was given on how the optimal system response could be obtained using relatively simple numerical optimization algorithms due to the smoothness of the optimization problem. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

INITIAL EXPLORATION OF A NOVEL TRANSIENT ARREST SYSTEM INVOLVING FUEL HEATING (Presentation)

A preliminary analysis on a novel accident response system to diminish the severity of supercritical transients was conducted. The novel accident response system, called the instant shock arrest system, involves using electricity to heat the nuclear fuel at the onset of a large accidental reactivity insertion. This system is specifically designed for reactors with metallic fuel, such that the fuel is capable of conducting electricity, and being resistively heated. A reactor dynamics model of the advanced test reactor was created using the point kinetics equations and a linear reactivity feedback model to simulate how the system would effect the maximum fuel temperatures experienced during the transient. Transients with the instant shock arrest system were compare to those without it. It was found that the instant shock arrest system initially heated the fuel more than the unaffected transient but the negative reactivity inserted from such heating was enough to lower the maximum fuel temperature experienced during the transient. After simulating six different accident scenarios with reactivity insertions ranging from 0.5 \$ to 1.3 \$, it was found that the an optimal system response could reduce peak fuel temperatures during the transient by 3.5% to 5%. Furthermore, discussion was given on how the optimal system response could be obtained using relatively simple numerical optimization algorithms due to the smoothness of the optimization problem.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Analysis of Lobe Power Calculator and Indication System with Physics and Cycle Based Models

The Advanced Test Reactor (ATR) at INL measures reactor power through two methods, thermal and Nitrogen-16 (N-16) activity. Water power calculator (WPC) is a thermal power system that measures flow and temperature to determine the thermal quadrant powers. The N-16 system utilizes a beta chamber detector that outputs reactivity levels to calculate lobe power through an algorithm called lobe power calculation and indication system (LPCIS). The LPCIS utilizes the N-16 system and the WPC system to determine reactor core power levels. The WPC provides accurate calculations of quadrant and total reactor thermal power. With the use of WPC measurements, thermal-to-N-16 (T2N) power ratios are produced to determine if the two indication systems agree on core power. Relative magnitude equations are used to utilize N-16 coefficients and multipliers to improve the indications of the LPCIS. This is crucial for maintaining safety limits. Currently, the LPCIS system uses linear equations and matrices to calculate lobe power through multipliers and coefficients. Advancements in technology and system upgrades have assisted system engineers at INL with the objective to reach a more dynamic system. In return the system demonstrates an increase in the accuracy of power reading while maintaining safety margins. The new proposed coefficient and multiple method implements a cycle specific, and a physics-based model intended for changing coefficients during operation. This calibration experiment focuses on power splitting, outer shim and neck shim, as well as fuel burning into the RDAS weighting factors. Results show that the Physics learning method yields a smaller error margin inside of the desired power range for the data set 166-A. This proved true for both constrained and unconstrained testing. This is most likely due to the physics data fitting approach that resulted in favorable coefficients and multipliers for 166-A. Continuing to improve the physics-based model will help improve the power accuracy of the LPCIS system.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Assessment of Calibration Approaches for the Stress Relaxation Test

Advanced manufacturing technologies are enabling the next generation of superalloys for extreme temperatures, pressures, and longevity. There is a need for Accelerated Creep Test (ACT) methods to reduce the time needed to implement new creep resistant superalloys. Without the development of ACTs, the qualification of new superalloys using conventional creep tests can take 11+ years corresponding to the 100,000-hour service life of long-lived turbomachinery components. A refined stress relaxation test (SRT) is introduced to succeed real-time conventional creep tests (CCTs). The stress relaxation test generates extensive stress versus creep-strain-rate data within a single specimen. Since the specimen is held below the elastic limit and relatively little creep deformation has occurred, several SRTs at different temperatures can be performed using a single specimen. To obtain data, specimens are rapidly loaded in force-control to just below the elastic limit. The specimen is switched to displacement-control and held at a fixed displacement while stress relaxation occurs. This generates a stress versus time curve that must be converted into stress versus creep-strain-rate. The stress versus creep-strain-rate curves are utilized to calculate creep activation energy and produce a creep-strain-rate master curve. The objective of this study is to determine the best calibration approach for SRT: regression of stress, regression of creep strain, and the finite difference of creep strain. It is recommended that the differentiation be applied; however, finite difference applied to the raw data produces inconsistent calculations of the creep-strain-rate. In this study, a MATLAB algorithm is written to evaluate all three calibration options where, stress and creep deformation are regressed into natural logarithm functions. The resulting creep-strain-rate curves are compared to determine the most reliable and consistent method. Analyzing and evaluating the results the regression of stress and regression or creep strain calibration methods prove to provide accurate results. These methods will help replace CCTs by conducting less expensive ACTs, while increasing the overall efficiency of the development of new materials.

36 MATERIALS SCIENCE↗

Supercooled Liquid Water Detection Capabilities from Ka-Band Doppler Profiling Radars: Moment-Based Algorithm Formulation and Assessment

The occurrence of supercooled liquid water in mixed-phase cloud (MPC) affects their cloud microphysical and radiative properties. The prevalence of MPCs in the mid- and high latitudes translates these effects to significant contributions to Earth’s radiative balance and hydrological cycle. The current study develops and assesses a radar-only, moment-based phase partition technique for the demarcation of supercooled liquid water volumes in arctic, MPC conditions. The study utilizes observations from the Ka band profiling radar, the collocated high spectral resolution lidar, and ambient temperature profiles from radio sounding deployments following a statistical analysis of 5.5 years of data (January 2014–May 2019) from the Atmospheric Radiation Measurement observatory at the North Slope of Alaska. The ice/liquid phase partition occurs via a per-pixel, neighborhood-dependent algorithm based on the premise that the partitioning can be deduced by examining the mean values of locally sampled probability distributions of radar-based observables and then compare those against the means of climatologically derived, per-phase probability distributions. Analyzed radar observables include linear depolarization ratio (LDR), spectral width, and vertical gradients of reflectivity factor and radial velocity corrected for vertical air motion. Results highlight that the optimal supercooled liquid water detection skill levels are realized for the radar variable combination of spectral width and reflectivity vertical gradient, suggesting that radar-based polarimetry, in the absence of full LDR spectra, is not as critical as Doppler capabilities. The cloud phase masking technique is proven particularly reliable when applied to cloud tops with an Equitable Threat Score (ETS) of 65%; the detection of embedded supercooled layers remains much more uncertain (ETS = 27%).

54 ENVIRONMENTAL SCIENCES↗

Convection equation modeling: A non-iterative direct matrix solution algorithm for use with SINDA

The determination of the boundary conditions for a component-level analysis, applying discrete finite element and finite difference modeling techniques often requires an analysis of complex coupled phenomenon that cannot be described algebraically. For example, an analysis of the temperature field of a coldplate surface with an integral fluid loop requires a solution to the parabolic heat equation and also requires the boundary conditions that describe the local fluid temperature. However, the local fluid temperature is described by a convection equation that can only be solved with the knowledge of the locally-coupled coldplate temperatures. Generally speaking, it is not computationally efficient, and sometimes, not even possible to perform a direct, coupled phenomenon analysis of the component-level and boundary condition models within a single analysis code. An alternative is to perform a disjoint analysis, but transmit the necessary information between models during the simulation to provide an indirect coupling. For this approach to be effective, the component-level model retains full detail while the boundary condition model is simplified to provide a fast, first-order prediction of the phenomenon in question. Specifically for the present study, the coldplate structure is analyzed with a discrete, numerical model (SINDA) while the fluid loop convection equation is analyzed with a discrete, analytical model (direct matrix solution). This indirect coupling allows a satisfactory prediction of the boundary condition, while not subjugating the overall computational efficiency of the component-level analysis. In the present study a discussion of the complete analysis of the derivation and direct matrix solution algorithm of the convection equation is presented. Discretization is analyzed and discussed to extend of solution accuracy, stability and computation speed. Case studies considering a pulsed and harmonic inlet disturbance to the fluid loop are analyzed to assist in the discussion of numerical dissipation and accuracy. In addition, the issues of code melding or integration with standard class solvers such as SINDA are discussed to advise the user of the potential problems to be encountered.

Schrage, Dean S.↗

Evaluation of Spaceborne L-band Radiometer Measurements for Terrestrial Freeze/Thaw Retrievals in Canada

The landscape freeze/thaw (FT) state has an important impact on the surface energy balance, carbon fluxes, and hydrologic processes; the timing of spring melt is linked to active layer dynamics in permafrost areas. L-band (1.4 GHz) microwave emission could allow the monitoring of surface state dynamics due to its sensitivity to the pronounced permittivity difference between frozen and thawed soil. The aim of this paper is to evaluate the performance of both Aquarius and Soil Moisture and Ocean Salinity (SMOS) L-band passive microwave measurements using a polarization ratio-based algorithm for landscape FT monitoring. Weekly L-band satellite observations are compared with a large set of reference data at 48 sites across Canada spanning three environments: tundra, boreal forest, and prairies. The reference data include in situ measurements of soil temperature (Tsoil) and air temperature (Tair), and Moderate Resolution Imaging Spectroradiometer (MODIS) land surface temperature (LST) and snow cover area (SCA) products. Results show generally good agreement between Lband FT detection and the surface state estimated from four reference datasets. The best apparent accuracies for all seasons are obtained using Tair as the reference. Aquarius radiometer 2 (incidence angle of 39.6) data gives the best accuracies (90.8), while for SMOS the best results (87.8 of accuracy) are obtained at higher incidence angles (55- 60). The FT algorithm identifies both freeze onset and end with a delay of about one week in tundra and two weeks in forest and prairies, when compared to Tair. The analysis shows a stronger FT signal at tundra sites due to the typically clean transitions between consistently frozen and thawed conditions (and vice versa) and the absence of surface vegetation. Results in the prairies were poorer because of the influence of vegetation growth in summer (which decreases the polarization ratio) and the high frequency of ephemeral thaw events during winter. Freeze onset and end maps created from the same algorithm applied to SMOS and Aquarius measurements characterize similar FT patterns over Canada. This study shows the potential of using L-band spaceborne observations for FT monitoring, but underlines some limitations due to ice crusts in the snowpack, liquid water content in snow cover during the spring freeze to thaw transition, and vegetation growth.

soil freeze/thaw↗

Application of thermal model for pan evaporation to the hydrology of a defined medium, the sponge

A technique is presented which estimates pan evaporation from the commonly observed values of daily maximum and minimum air temperatures. These two variables are transformed to saturation vapor pressure equivalents which are used in a simple linear regression model. The model provides reasonably accurate estimates of pan evaporation rates over a large geographic area. The derived evaporation algorithm is combined with precipitation to obtain a simple moisture variable. A hypothetical medium with a capacity of 8 inches of water is initialized at 4 inches. The medium behaves like a sponge: it absorbs all incident precipitation, with runoff or drainage occurring only after it is saturated. Water is lost from this simple system through evaporation just as from a Class A pan, but at a rate proportional to its degree of saturation. The contents of the sponge is a moisture index calculated from only the maximum and minium temperatures and precipitation.

Trenchard, M. H.↗