Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “potential inversion”

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.

1,109 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Integration of Soft Data Into Geostatistical Simulation of Categorical Variables

Uncertain or indirect “soft” data, such as geologic interpretation, driller’s logs, geophysical logs or imaging, offer potential constraints or “soft conditioning” to stochastic models of discrete categorical subsurface variables in hydrogeology such as hydrofacies. Previous bivariate geostatistical simulation algorithms have not fully addressed the impact of data uncertainty in formulation of the (co) kriging equations and the objective function in simulated annealing (or quenching). This paper introduces the geostatistical simulation code tsim-s, which accounts for categorical data uncertainty through a data “hardness” parameter. In generating geostatistical realizations with tsim-s, the uncertainty inherent to soft conditioning is factored into both 1) the data declustering and spatial correlation functions in cokriging and 2) the acceptance probability for change of category in simulated quenching. The degree or sensitivity to which soft data conditions a realization as a function of hardness can be quantified by mapping category probabilities derived from multiple realizations. In addition to point or borehole data, arrays of data (e.g., as derived from a depth-dependency function, probability map, or “prior realization”) can be used as soft conditioning. The tsim-s algorithm provides a theoretically sound and general framework for integrating datasets of variable location, resolution, and uncertainty into geostatistical simulation of categorical variables. A practical example shows how tsim-s is capable of generating a large-scale three-dimensional simulation including curvilinear features.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1992

Characterization of radiation-induced damage in high performance polymers by electron paramagnetic resonance imaging spectroscopy

The potential for long-term human activity beyond the Earth's protective magnetosphere is limited in part by the lack of detailed information on the effectiveness and performance of existing structural materials to shield the crew and spacecraft from highly penetrating space radiations. The two radiations of greatest concern are high energy protons emitted during solar flares and galactic cosmic rays which are energetic ions ranging from protons to highly oxidized iron. Although the interactions of such high-energy radiations with matter are not completely understood at this time, the effects of the incident radiation are clearly expected to include the formation of paramagnetic spin centers via ionization and bond-scission reactions in the molecular matrices of structural materials. Since this type of radiation damage is readily characterized by Electron Paramagnetic Resonance (EPR) spectroscopy, the NASA Langley Research Center EPR system was repaired and brought on-line during the 1991 ASEE term. A major goal of the 1992 ASEE term was to adapt the existing core of the LaRC EPR system to meet the requirements for EPR Imaging--a powerful new technique which provides detailed information on the internal structure of materials by mapping the spatial distribution of unpaired spin density in bulk media. Major impetus for this adaptation arises from the fact that information derived from EPRI complements other methods such as scanning electron microscopy which primarily characterize surface phenomena. The modification of the EPR system has been initiated by the construction of specially designed, counterwound Helmholtz coils which will be mounted on the main EPR electromagnet. The specifications of the coils have been set to achieve a static linear magnetic field gradient of 10 gauss/mm/amp along the principal (Z) axis of the Zeeman field. Construction is also in progress of a paramagnetic standard in which the spin distribution is known in all three dimensions. This sample will be used to assess the linearity of the magnetic field gradient and to ensure authentic image reconstruction. A second major task was to secure the computer capability to enable image reconstruction from projection data generated by the magnetic field gradients. To this end, commercially available and public domain software packages which perform inverse Fourier Transform and convoluted (filtered) back projection functions are being integrated into the existing EPR data processing system.

Suleman, Naushadalli K.↗
DOE OSTI2024

User Impressions and Gait Analysis of Exoskeleton Device Usage in Generalized Tank Farm Activities

Tank farm workers involved in nuclear cleanup activities perform physically demanding tasks, typically while wearing heavy personal protective equipment (PPE). Exoskeleton devices have the potential to bring considerable benefit to this industry but have not been thoroughly studied in the context of nuclear cleanup. In this paper, we examine the performance of exoskeletons during a series of tasks emulating jobs performed on tank farms while participants wore PPE commonly deployed by tank farm workers. The goal of this study was to evaluate the effects of commercially available lower-body exoskeletons on a user’s gait kinematics and user perceptions. Three participants each tested three lower-body exoskeletons in a 70-min protocol consisting of level treadmill walking, incline treadmill walking, weighted treadmill walking, a weight lifting session, and a hand tool dexterity task. Results were compared to a no exoskeleton baseline condition and evaluated as individual case studies. The three participants showed a wide spectrum of user preferences and adaptations toward the devices. Individual case studies revealed that some users quickly adapted to select devices for certain tasks while others remained hesitant to use the devices. Temporal effects on gait change and perception were also observed for select participants in device usage over the course of the device session. Here, device benefit varied between tasks, but no conclusive aggregate trends were observed across devices for all tasks. Evidence suggests that device benefits observed for specific tasks may have been overshadowed by the wide array of tasks used in the protocol.

Tank farm↗
NASA NTRS2020

The AERONET Version 3 Aerosol Retrieval Algorithm, Associated Uncertainties and Comparisons to Version 2

The Aerosol Robotic Network (AERONET) Version 3 (V3) aerosol retrieval algorithm is described, which is based on the Version 2 (V2) algorithm with numerous updates. Comparisons of V3 aerosol retrievals to those of V2 are presented, along with a new approach to estimate uncertainties in many of the retrieved aerosol parameters. Changes in the V3 aerosol retrieval algorithm include (1) a new polarized radiative transfer code (RTC), which replaced the scalar RTC of V2, (2) detailed characterization of gas absorption by adding NO2 and H2O to specify total gas absorption in the atmospheric column, specification of vertical profiles of all the atmospheric species, (3) new bidirectional reflectance distribution function (BRDF) parameters for land sites adopted from the MODIS BRDF/Albedo product, (4) a new version of the extraterrestrial solar flux spectrum, and (5) a new temperature correction procedure of both direct Sun and sky radiance measurements. The potential effect of each change in V3 on single scattering albedo (SSA) retrievals was analyzed. The operational almucantar retrievals of V2 versus V3 were compared for four AERONET sites: GSFC, Mezaira, Mongu, and Kanpur. Analysis showed very good agreement in retrieved parameters of the size distributions. Comparisons of SSA retrievals for dust aerosols (Mezaira) showed a good agreement in 440 nm SSA, while for longer wavelengths V3 SSAs are systematically higher than those of V2, with the largest mean difference at 675 nm due to cumulative effects of both extraterrestrial solar flux and BRDF changes. For non-dust aerosols, the largest SSA deviation is at 675 nm due to differences in extraterrestrial solar flux spectrums used in each version. Further, the SSA 675 nm mean differences are very different for weakly (GSFC) and strongly (Mongu) absorbing aerosols, which is explained by the lower sensitivity to a bias in aerosol scattering optical depth by less absorbing aerosols. A new hybrid (HYB) sky radiance measurement scan is introduced and discussed. The HYB combines features of scans in two different planes to maximize the range of scattering angles and achieve scan symmetry, thereby allowing for cloud screening and spatial averaging, which is an advantage over the principal plane scan that lacks robust symmetry. We show that due to an extended range of scattering angles, HYB SSA retrievals for dust aerosols exhibit smaller variability with solar zenith angles (SZAs) than those of almucantar (ALM), which allows extension of HYB SSA retrievals to SZAs less than 50∘ to as small as 25∘. The comparison of SSA retrievals from closely time-matched HYB and ALM scans in the 50 to 75∘ SZA range showed good agreement with the differences below ∼0.005. We also present an approach to estimate retrieval uncertainties which utilizes the variability in retrieved parameters generated by perturbing both measurements and auxiliary input parameters as a proxy for retrieval uncertainty. The perturbations in measurements and auxiliary inputs are assumed as estimated biases in aerosol optical depth (AOD), radiometric calibration of sky radiances combined with solar spectral irradiance, and surface reflectance. For each set of Level 2 Sun/sky radiometer observations, 27 inputs corresponding to 27 combinations of biases were produced and separately inverted to generate the following statistics of the inversion results: average, standard deviation, minimum and maximum values. From these statistics, standard deviation (labeled U27) is used as a proxy for estimated uncertainty, and a lookup table (LUT) approach was implemented to reduce the computational time. The U27 climatological LUT was generated from the entire AERONET almucantar (1993–2018) and hybrid (2014–2018) scan databases by binning U27s in AOD (440 nm), Angström exponent (AE, 440–870 nm), and SSA (440, 675, 870, 1020 nm). Using this LUT approach, the uncertainty estimates U27 for each individual V3 Level 2 retrieval can be obtained by interpolation using the corresponding measured and inverted combination of AOD, AE, and SSA.

Aliaksandr Siniuk↗
NASA NTRS1997

Effects of Fuel Preheat on Soot Formation in Microgravity Laminar Diffusion Flames

Nonbuoyant flames offer themselves as an attractive and promising platform to gain a better understanding of soot mechanisms. The effects of buoyancy can be eliminated temporarily in drop towers which sustain brief intervals of reduced gravity-typically lower than 10(exp -3)g- extending up to several seconds at a time. Microgravity facilities have been employed to show that nonbuoyant flames are longer, wider and sootier than their normal-gravity counterparts. Sunderland et al. recently verified the existence of smoke point in laminar nonbuoyant flames. As reported, microgravity flames operating above their smoke point displayed a blunt tip and much broader soot-containing regimes in comparison to their buoyant counterparts. Mortazavi et al. established that residence times in microgravity laminar jet diffusion flames with Re=0(100) tend to be proportional to burner diameter and inversely proportional to burner exit velocity. This offers the capability to alter residence times in nonbuoyant laminar jet diffusion flames when varying the burner exit diameters and velocities. Megaridis et al. presented a quantitative definition of the soot-field structure within laminar microgravity jet diffusion flames which operated well above their smoke point. The experimental methodology involved a full-field laser-light extinction technique and jet diffusion flames of nitrogen-diluted (50% vol.) acetylene fuel burning in quiescent air at atmospheric pressure. The work was conducted at the 2.2s drop tower of the NASA Lewis Research Center (NASA-LeRC). Parallel work on 1-g flames was also presented in (6) to facilitate comparisons on the effect of gravity on the soot fields. As reported, the soot spatial distributions in 0-g flames did not change in a detectable manner after 1s within a typical 2.2s experiment. During that period, the soot field was shown to sustain a pronounced annular structure throughout the luminous nonbuoyant-flame zone. The maximum soot volume fraction measured at 0-g was nearly a factor of two higher than that at 1-g, thus confirming the enhanced sooting tendency of nonbuoyant flames. Greenberg and Ku presented a similar study and reported trends that matched those of for the 50% (vol.) nitrogen-diluted acetylene fuel. Furthermore, they examined pure acetylene flames and reported similar trends with respect to the influence of gravity on maximum soot volume fractions and flame cross-section-averaged soot loadings. Both studies clearly demonstrated the improved spatial resolution of microgravity flames compared to their normal-gravity counterparts. The current study evaluates the influence of moderate fuel preheat on soot formation within 0-g laminar gas jet diffusion flames. While fuel temperature variations have little influence on residence times in 1-g, they have a much more significant effect in 0-g. The primary objective of this program is to quantify this effect and its consequences on sooting by comparing soot volume fraction distributions under preheated and unpreheated-fuel conditions. Furthermore, the current work aims at expanding the limited soot database available for nonbuoyant flames. Soot fields in such flames can be used to perform additional tests of recently developed soot sub-models which have the potential to become powerful predicting tools in combustion design.

Konsur, Bogdan↗
NASA NTRSDate not supplied

Advancing Satellite-Constrained Modeled Air-Sea CO 2 Fluxes With a Focus on the Strength of the Southern Ocean Carbon Sink

Challenge and Motivation: The ocean plays a critical role in mitigating climate change by removing approximately a quarter of annual anthropogenic CO 2 emissions from the atmosphere. Model-based estimates point to the Southern Ocean as a key marine region, responsible for approximately 40 % of the anthropogenic carbon uptake by the global ocean. However, the contemporary strength of the Southern Ocean carbon sink has recently come into question. On the one hand, airborne-based observations of atmospheric CO 2 gradients indicate that the Southern Ocean represents a strong net sink of atmospheric CO 2 , consistent in magnitude with atmospheric inversion estimates and surface-ocean partial pressure of CO 2 (pCO 2 )-based products. On the other hand, estimates of pCO 2 based on in situ pH measurements taken by biogeochemical profiling floats yield strong wintertime outgassing fluxes that greatly reduce the Southern Ocean’s annually integrated CO 2 uptake. This uncertainty in the strength of the Southern Ocean air-sea CO 2 flux and its role in the global carbon cycle hinders our ability to constrain global carbon fluxes, one of the major goals of NASA’s Carbon Monitoring System (CMS). Opportunity: The NASA Ocean Biogeochemical Model (NOBM) produces near-global pCO 2 and air-sea CO 2 flux estimates that are currently included into the NASA’s Goddard Earth Observing System (GEOS) models in support of the CMS effort to monitor global carbon fluxes. The NOBM assimilates ocean color data to improve the representation of biogeochemical fluxes and overcome spatial and temporal gaps in the space-based retrievals. Here, we propose to advance the satellite-constrained flux estimates by investigating the uncertainties in the Southern Ocean air-sea CO 2 flux produced by the NOBM, and assess the value that remote sensing ocean color data can have in providing improved estimates of carbon fluxes in the ocean. Our proposed work includes the delivery of refined in situ float-based carbon fluxes to serve as a constraint on the model-based estimates. Taking advantage of the model’s integration of satellite ocean color data to represent multiple phytoplankton groups, we propose to deliver maps of biogenic carbon export specific to each modeled phytoplankton type and investigate the role of ecological plankton complexity in regulating marine carbon uptake and export. Goals: (a) Delivery of seasonally-adjusted float-based Southern Ocean air-sea CO 2 fluxes: We will produce updated and improved float-based air-sea CO 2 fluxes that will serve as a bias-reduced float-based constraint to evaluate our model-based estimates of the NOBM. (b) Investigation of uncertainties in Southern Ocean air-sea CO 2 flux from the NOBM: Modeled air-sea carbon fluxes will be evaluated against the updated float product as well as ship- and airborne-based data to identify uncertainties and potential model deficiencies. (c) Delivery of model-based carbon export partitioning by phytoplankton functional types (PFTs): We will produce depth-resolved maps of particulate organic export production integrated for all phytoplankton groups and allocated to each individual PFT in the model. The expected significance of this goal is to quantify the role that the functional-oriented diversity in phytoplankton groups represented in the NOBM plays in regulating air-sea CO 2 fluxes in the Southern Ocean.

Lionel A. Arteaga↗
NASA NTRS2023

Propulsion-Airframe Integration for Conceptual Redesign of a Low-Boom Supersonic Transport

A low-boom supersonic transport was generated in a previous low-boom multidisciplinary optimization (MDO) study that used computational fluid dynamics (CFD) off-body pressure to compute the undertrack sonic boom ground signature and calibrated low-fidelity aerodynamics analyses to compute the mission performance metrics. This low-boom aircraft, referred to as the Mach 1.7 40-PAX concept, can carry 40 passengers for a low-boom overland mission with cruise Mach 1.7 for the airport pairs over the continental US and has the potential to achieve an undertrack sonic boom ground noise level below 70 perceived level of decibels at the start of overland cruise (SOC). The engine for the Mach 1.7 40-PAX concept was designed using the Numerical Propulsion System Simulation (NPSS) for the mission analysis and modeled as a flow-through nacelle for CFD-based sonic boom analysis. To understand how the engine plume affects the undertrack ground signature, the Mach 1.7 40-PAX concept is redesigned in this paper after replacing the flow-through nacelles by CFD engines for sonic boom analysis using aeropropulsive CFD simulation. A process is developed for approximation of the NPSS engine at SOC by a CFD engine for aeropropulsive CFD simulation. The generated CFD engine has the identical nozzle boundary conditions and approximately the same mass flow and thrust as those of the NPSS engine at SOC. Then, the outer mold line of the configuration with the CFD engines is redesigned to approximately restore the low-boom characteristics of the Mach 1.7 40-PAX concept at SOC. Finally, the CFD simulation data for the redesigned concept with the CFD engines is used to calibrate the low-fidelity aerodynamic analyses for mission analysis of this concept. The cyclic dependency of the involved disciplinary analyses is resolved using an iteration method for a consistent coupling of the mission analysis, NPSS engine analysis, and low-boom redesign using the aeropropulsive CFD simulation. This low-boom redesign study is used as an example to demonstrate how the propulsion-airframe integration could be implemented for conceptual design of supersonic transports that satisfy both the low-boom and mission performance requirements.

Low-boom supersonic transport↗
DOE OSTI2023

Encoding nonlinear and unsteady aerodynamics of limit cycle oscillations using nonlinear sparse Bayesian learning

This article investigates the applicability of a recently proposed, nonlinear sparse Bayesian learning (NSBL) algorithm to identify and estimate the complex aerodynamics of limit cycle oscillations. NSBL provides a semi-analytical framework for determining the data-optimal sparse model nested within a (potentially) over-parameterized model. This is particularly relevant to nonlinear dynamical systems where modelling approaches involve the use of physics-based and data-driven components. In such cases, the data-driven components, where analytical descriptions of the physical processes are not readily available, are often prone to overfitting, meaning that the empirical aspects of these models will often involve the calibration of an unnecessarily large number of parameters. While an overparameterized model may fit the observed data well, such models may be inadequate for making predictions in regimes that are different from those wherein the data were recorded. In view of this, it is desirable to not only calibrate the model parameters, but also identify the optimal compromise between data fit and model complexity. In this article, we exhibit the optimal model discovery for an aeroelastic system wherein the structural dynamics are well-known and described by a differential equation model, coupled with a semi-empirical aerodynamic model for laminar separation flutter, resulting in low-amplitude limit cycle oscillations (LCO). To illustrate the performance of the algorithm, in this article, we use synthetic data and demonstrate the ability of the algorithm to correctly rediscover the optimal model and model parameters, given a known data-generating model. The synthetic data are generated from a forward simulation of a known differential equation model with parameters selected so as to mimic the dynamics observed in wind-tunnel experiments. Subsequently, we demonstrate the performance of the algorithm for model selection using noisy LCO data from wind tunnel experiments. As there is no ground truth available for the experimental data case, we provide a comparison between NSBL and Bayesian model selection to validate the results, and demonstrate the use of NSBL as an efficient alternative to traditional methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Big Sky Regional Carbon Sequestration Partnership (Phase III Final Scientific/Technical Report)

The Big Sky Carbon Sequestration Partnership (BSCSP) pursued a Phase III demonstration project at Kevin Dome in north central Montana. Kevin Dome covers approximately 700 square miles and is a naturally occurring CO 2 reservoir that is flanked by oil and gas fields. The carbon dioxide (CO 2 ) is in the upper Devonian Duperow (carbonate) formation and does not reach the spill point of the dome; therefore, the dome has potential as a CO 2 sequestration reservoir, a CO 2 supply, or as both if anthropogenic sources and enhanced oil recovery (EOR) operations are associated with the dome. Kevin Dome could potentially act as a buffer to continue accepting anthropogenic CO 2 when EOR flooding operations are interrupted or completed. The project objective was to produce one million tonnes of CO 2 from the gas cap of Kevin Dome, pipe it laterally, inject, and re-store it in the brine leg of the same formation to test the hub / buffer storage concept. This was to be accomplished by drilling up to five production wells, building a short pipeline and compression facilities, and drilling an injection well and several monitoring wells. BSCSP commenced outreach and site characterization activities including acquiring baseline data for near-surface insurance monitoring, acquiring 3-dimensional, 9-component surface seismic over the project area, drilling two test wells (one in the production area and one in the injection area), coring key intervals, and performing comprehensive logging. Well tests of those wells revealed two barriers to the project. The production (Danielson 33-17) well was perforated in multiple zones but failed to produce any significant CO 2 . This was despite being drilled in the near vicinity of a historic well that had produced 3700 MCF per day in a drill stem test. Modeling indicated that this was likely due to a phase change during production causing a temperature drop resulting in hydrate and/or water ice formation that clogged the formation. Tests of the injection zone (Wallewein 22-1) well indicated total dissolved solids (TDS) slightly below the EPA required 10,000 parts per million (ppm) for a Class VI underground injection control permit. While the project was initiated before Class VI rules were promulgated, and this was an experimental project (seemingly qualified for a Class V permit), the Environmental Protection Agency (EPA) indicated that injection would require a Class VI permit. The low salinity result was unexpected as contours plotted based on regional formation water quality data indicated an expected TDS above 20,000 ppm, and wells between the recharge zone and the Wallewein well tested above 10,000 ppm. Faced with the inability to obtain an injection permit, the demonstration project could not proceed. However, the project had generated valuable samples and data on a large natural analog including 32 sq. mi. of 3-D, 9-C seismic, 430 ft. of carbonate core covering seven different depositional environments taken from areas with, and without the presence of CO 2 , 30 ft. of core of two caprocks, a tight carbonate and an anhydrite, a full set of modern logs on both wells, and well tests. DOE decided to re-scope the project around completing studies utilizing this data. This report covers both the initial scope and the re-scope (Task / Section numbers preceded with an R). While the report covers a wide range of project activities, highlights of this work include: Development of a geostatic model using neural nets to match well logs to facies and using multi-waveform seismic to inform reservoir heterogeneity; Unique mechanical testing of permeability – stress relationship in two caprock materials; Development of full waveform inversion to generate a high resolution velocity model; Model development for dual permeability (fracture and matrix) systems to better account for matrix-matrix interactions; Joint seismic wave inversion (including the first quadr-joint inversion) exhibiting better imaging of a challenging reservoir zone in stiff rock; Core flow and core flood results on a reactive carbonate; and Innovative laboratory measurements of seismic response of fractured core as a function of fluid fill.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2000

Characterization of Electronic Materials HgZnSe and HgZnTe Using Innovative and Conventional Techniques

HgZnSe and HgZnTe are electronic materials of interest for potential IR detector and focal plane array applications due to their improved strength and compositional stability over HgCdTe, but they are difficult to grow on Earth and to fully characterize. Conventional contact methods of characterization, such as Hall and van der Paw, although adequate for many situations are typically labor intensive and not entirely suitable where only very small samples are available. To adequately characterize and compare properties of electronic materials grown in low earth orbit with those grown on Earth, innovative techniques are needed that complement existing methods. This paper describes the implementation and test results of a unique non-contact method of characterizing uniformity, mobility, and carrier concentration together with results from conventional methods applied to HgZnSe and HgZnTe. The innovative method has advantages over conventional contact methods since it circumvents problems of possible contamination from alloying electrical contacts to a sample and also has the capability to map a sample. Non- destructive mapping, the determination of the carrier concentration and mobility at each place on a sample, provides a means to quantitatively compare, at high spatial resolution, effects of microgravity on electronic properties and uniformity of electronic materials grown in low-Earth orbit with Earth grown materials. The mapping technique described here uses a 1mm diameter polarized beam of radiation to probe the sample. Activation of a magnetic field, in which the sample is placed, causes the plane of polarization of the probe beam to rotate. This Faraday rotation is a function of the free carrier concentration and the band parameters of the material. Maps of carrier concentration, mobility, and transmission generated from measurements of the Faraday rotation angles over the temperature range from 300K to 77K will be presented. New information on band parameters, obtained by combining results from conventional Hall measurements of the free carrier concentration with Faraday rotation measurements, will also be presented. One example of how this type of information was derived is illustrated in the following figure which shows Faraday rotation vs wavelength modeled for Hg(l-x)ZnxSe at a temperature of 300K and x=0.07. The plasma contribution, total Faraday rotation, and interband contribution to the Faraday rotation, are designated in the Figure as del(p), FR tot, and del(i) respectively. Experimentally measured values of FR tot, each indicated by + , agree acceptably well with the model at the probe wavelength of 10.6 microns. The model shows that at the probe wavelength, practically all the rotation is due to the plasma component, which can be expressed as delta(sub p)= 2pi(e(sup 3))NBL/c(sup 2)nm*(sup 2) omega(sup 2). In this equation, delta(sub p) is the rotation angle due to the free carrier plasma, N is the free carrier concentration, B the magnetic field strength, L the thickness of the sample, n the index of refraction, omega the probe radiation frequency, c the speed of light, e the electron charge, and m* the effective mass. A measurement of N by conventional techniques, combined with a measurement of the Faraday rotation angle allows m* to be accurately determined since it is an inverse square function.

Tanton, George↗
NASA NTRS1998

A Conceptual Approach to Assimilating Remote Sensing Data to Improve Soil Moisture Profile Estimates in a Surface Flux/Hydrology Model: Overview - Part 1

Knowledge of the amount of water in the soil is of great importance to many earth science disciplines. Soil moisture is a key variable in controlling the exchange of water and energy between the land surface and the atmosphere. Thus, soil moisture information is valuable in a wide range of applications including weather and climate, runoff potential and flood control, early warning of droughts, irrigation, crop yield forecasting, soil erosion, reservoir management, geotechnical engineering, and water quality. Despite the importance of soil moisture information, widespread and continuous measurements of soil moisture are not possible today. Although many earth surface conditions can be measured from satellites, we still cannot adequately measure soil moisture from space. Research in soil moisture remote sensing began in the mid 1970s shortly after the surge in satellite development. Recent advances in remote sensing have shown that soil moisture can be measured, at least qualitatively, by several methods. Quantitative measurements of moisture in the soil surface layer have been most successful using both passive and active microwave remote sensing, although complications arise from surface roughness and vegetation type and density. Early attempts to measure soil moisture from space-borne microwave instruments were hindered by what is now considered sub-optimal wavelengths (shorter than 5 cm) and the coarse spatial resolution of the measurements. L-band frequencies between 1 and 3 GHz (10-30 cm) have been deemed optimal for detection of soil moisture in the upper few centimeters of soil. The Electronically Steered Thinned Array Radiometer (ESTAR), an aircraft-based instrument operating a 1,4 GHz, has shown great promise for soil moisture determination. Initiatives are underway to develop a similar instrument for space. Existing space-borne synthetic aperture radars (SARS) operating at C- and L-band have also shown some potential to detect surface wetness. The advantage of radar is its much higher resolution than passive microwave systems, but it is currently hampered by surface roughness effects and the lack of a good algorithm based on a single frequency and single polarization. In addition, its repeat frequency is generally low (about 40 days). In the meantime, two new radiometers offer some hope for remote sensing of soil moisture from space. The Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI), launched in November 1997, possesses a 10.65 GHz channel and the Advanced Microwave Scanning Radiometer (AMSR) on both the ADEOS-11 and Earth Observing System AM-1 platforms to be launched in 1999 possesses a 6.9 GHz channel. Aside from issues about interference from vegetation, the coarse resolution of these data will provide considerable challenges pertaining to their application. The resolution of TMI is about 45 km and that of AMSR is about 70 km. These resolutions are grossly inconsistent with the scale of soil moisture processes and the spatial variability of factors that control soil moisture. Scale disparities such as these are forcing us to rethink how we assimilate data of various scales in hydrologic models. Of particular interest is how to assimilate soil moisture data by reconciling the scale disparity between what we can expect from present and future remote sensing measurements of soil moisture and modeling soil moisture processes. It is because of this disparity between the resolution of space-based sensors and the scale of data needed for capturing the spatial variability of soil moisture and related properties that remote sensing of soil moisture has not met with more widespread success. Within a single footprint of current sensors at the wavelengths optimal for this application, in most cases there is enormous heterogeneity in soil moisture created by differences in landcover, soils and topography, as well as variability in antecedent precipitation. It is difficult to interpret the meaning of 'mean' soil moisture under such conditions and even more difficult to apply such a value. Because of the non-linear relationships between near-surface soil moisture and other variables of interest, such as surface energy fluxes and runoff, mean soil moisture has little applicability at such large scales. It is for these reasons that the use of remote sensing in conjunction with a hydrologic model appears to be of benefit in capturing the complete spatial and temporal structure of soil moisture. This paper is Part I of a four-part series describing a method for intermittently assimilating remotely-sensed soil moisture information to improve performance of a distributed land surface hydrology model. The method, summarized in section II, involves the following components, each of which is detailed in the indicated section of the paper or subsequent papers in this series: Forward radiative transfer model methods (section II and Part IV); Use of a Kalman filter to assimilate remotely-sensed soil moisture estimates with the model profile (section II and Part IV); Application of a soil hydrology model to capture the continuous evolution of the soil moisture profile within and below the root zone (section III); Statistical aggregation techniques (section IV and Part II); Disaggregation techniques using a neural network approach (section IV and Part III); and Maximum likelihood and Bayesian algorithms for inversely solving for the soil moisture profile in the upper few cm (Part IV).

Crosson, William L.↗