Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear inverse modeling”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Inverse model based error detection in beamline optics

Optics tuning in transfer lines and LINACs can be challenging due to the fact that multiple combinations of machine settings can lead to the same diagnostic output. Moreover, the lack of a periodic solution can limit the ability to infer optics in the same way as rings from BPM signals. Model based approaches are often used to assist with the optics tuning in combination with optimization or parameter estimation. Here we have developed a novel approach using machine learning inverse models trained on a known configuration to detect variations in quadrupole settings without explicitly including them in the model. This paper shows a comparison of neural network models and linear models on both a simulation based study and experimental studies conducted at the AGS to RHIC transfer line at Brookhaven National Lab.

43 PARTICLE ACCELERATORS↗
NASA NTRS1999

Measuring Greenland Ice Mass Variation With Gravity Recovery and the Climate Experiment Gravity and GPS

The response of the Greenland ice sheet to climate change could significantly alter sea level. The ice sheet was much thicker at the last glacial maximum. To gain insight into the global change process and the future trend, it is important to evaluate the ice mass variation as a function of time and space. The Gravity Recovery and Climate Experiment (GRACE) mission to fly in 2001 for 5 years will measure gravity changes associated with the current ice variation and the solid earth's response to past variations. Our objective is to assess the separability of different change sources, accuracy and resolution in the mass variation determination by the new gravity data and possible Global Positioning System (GPS) bedrock uplift measurements. We use a reference parameter state that follows a dynamic ice model for current mass variation and a variant of the Tushingham and Peltier ICE-3G deglaciation model for historical deglaciation. The current linear trend is also assumed to have started 5 kyr ago. The Earth model is fixed as preliminary reference Earth model (PREM) with four viscoelastic layers. A discrete Bayesian inverse algorithm is developed employing an isotropic Gaussian a priori covariance function over the ice sheet and time. We use data noise predicted by the University of Texas and JPL for major GRACE error sources. A 2 mm/yr uplift uncertainty is assumed for GPS occupation time of 5 years. We then carry out covariance analysis and inverse simulation using GRACE geoid coefficients up to degree 180 in conjunction with a number of GPS uplift rates. Present-day ice mass variation and historical deglaciation are solved simultaneously over 146 grids of roughly 110 km x 110 km and with 6 time increments of 3 kyr each, along with a common starting epoch of the current trend. For present-day ice thickness change, the covariance analysis using GRACE geoid data alone results in a root mean square (RMS) posterior root variance of 2.6 cm/yr, with fairly large a priori uncertainties in the parameters and a Gaussian correlation length of 350 km. Simulated inverse can successfully recover most features in the reference present-day change. The RMS difference between them over the grids is 2.8 cm/yr. The RMS difference becomes 1.1 cm/yr when both are averaged with a half Gaussian wavelength of 150 km. With a fixed Earth model, GRACE alone can separate the geoid signals due to past and current load fairly well. Shown are the reference geoid signatures of direct and elastic effects of the current trend, the viscoelastic effect of the same trend starting from 5 kyr ago, the Post Glacial Rebound (PGR), and the predicted GRACE geoid error. The difference between the reference and inverse modeled total viscoelastic signatures is also shown. Although past and current ice mass variations are allowed the same spatial scale, their geoid signals have different spatial patterns. GPS data can contribute to the ice mass determination as well. Additional information is contained in the original.

Wu, Xiao-Ping↗
NASA NTRS2016

The Relaxation Matrix for Symmetric Tops with Inversion Symmetry. I. Effects of Line Coupling on Self-Broadened v (sub 1) and Pure Rotational Bands of NH3

The Robert-Bonamy formalism has been commonly used to calculate half-widths and shifts of spectral lines for decades. This formalism is based on several approximations. Among them, two have not been fully addressed: the isolated line approximation and the neglect of coupling between the translational and internal motions. Recently, we have shown that the isolated line approximation is not necessary in developing semi-classical line shape theories. Based on this progress, we have been able to develop a new formalism that enables not only to reduce uncertainties on calculated half-widths and shifts, but also to model line mixing effects on spectra starting from the knowledge of the intermolecular potential. In our previous studies, the new formalism had been applied to linear and asymmetric-top molecules. In the present study, the method has been extended to symmetric-top molecules with inversion symmetry. As expected, the inversion splitting induces a complete failure of the isolated line approximation. We have calculated the complex relaxation matrices of selfbroadened NH3. The half-widths and shifts in the ν1 and the pure rotational bands are reported in the present paper. When compared with measurements, the calculated half-widths match the experimental data very well, since the inapplicable isolated line approximation has been removed. With respect to the shifts, only qualitative results are obtained and discussed. Calculated off-diagonal elements of the relaxation matrix and a comparison with the observed line mixing effects are reported in the companion paper (Paper II).

Trajectory models↗
DOE OSTI2024

Controlling cantilevered adaptive X-ray mirrors

Modeling the behavior of a prototype cantilevered X-ray adaptive mirror (held from one end) demonstrates its potential for use on high-performance X-ray beamlines. Similar adaptive mirrors are used on X-ray beamlines to compensate optical aberrations, control wavefronts and tune mirror focal distances at will. Controlled by 1D arrays of piezoceramic actuators, these glancing-incidence mirrors can provide nanometre-scale surface shape adjustment capabilities. However, significant engineering challenges remain for mounting them with low distortion and low environmental sensitivity. Finite-element analysis is used to predict the micron-scale full actuation surface shape from each channel and then linear modeling is applied to investigate the mirrors' ability to reach target profiles. Using either uniform or arbitrary spatial weighting, actuator voltages are optimized using a Moore–Penrose matrix inverse, or pseudoinverse, revealing a spatial dependence on the shape fitting with increasing fidelity farther from the mount.

47 OTHER INSTRUMENTATION↗
NASA NTRS2013

Analytical Algorithms to Quantify the Uncertainty in Remaining Useful Life Prediction

This paper investigates the use of analytical algorithms to quantify the uncertainty in the remaining useful life (RUL) estimate of components used in aerospace applications. The prediction of RUL is affected by several sources of uncertainty and it is important to systematically quantify their combined effect by computing the uncertainty in the RUL prediction in order to aid risk assessment, risk mitigation, and decisionmaking. While sampling-based algorithms have been conventionally used for quantifying the uncertainty in RUL, analytical algorithms are computationally cheaper and sometimes, are better suited for online decision-making. While exact analytical algorithms are available only for certain special cases (for e.g., linear models with Gaussian variables), effective approximations can be made using the the first-order second moment method (FOSM), the first-order reliability method (FORM), and the inverse first-order reliability method (Inverse FORM). These methods can be used not only to calculate the entire probability distribution of RUL but also to obtain probability bounds on RUL. This paper explains these three methods in detail and illustrates them using the state-space model of a lithium-ion battery.

Sankararaman, Shankar↗
NASA NTRS2020

The Fast Radio Burst Luminosity Function and Death Line in the Low-Twist Magnetar Model

We explore the burst energy distribution of fast radio bursts(FRBs)in the low-twist magnetar model ofWadiasingh & Timokhin(WT19). Motivated by the power-lawfluence distributions of FRB121102, we proposean elementary model for the FRB luminosity function of individual repeaters with an inversion protocol thatdirectly relates the power-law distribution index of magnetar short burstfluences to that for FRBs. The protocolindicates that the FRB energy scales virtually linearly with crust/field dislocation amplitude, if magnetar shortbursts prevail in the magnetoelastic regime. Charge starvation in the magnetosphere during bursts(required inWT19)for individual repeaters implies the predicted burstfluence distribution is narrow,3 decades for yieldingstrains and oscillation frequencies feasible in magnetar crusts. Requiring magnetic confinement and chargestarvation, we obtain a death line for FRBs, which segregates magnetars from the normal pulsar population,suggesting only the former will host recurrent FRBs. We convolve the burst energy distribution for individualmagnetars to define the distribution of luminosities in evolved magnetar populations. The broken power-lawluminosity function’s low-energy character depends on the population model, while the high-energy index tracesthat of individual repeaters. Independent of the evolved population, the broken power-law isotropic-equivalentenergy/luminosity function peaks at∼1037–1040erg with a low-energy cutoff at∼1037erg. Lastly, we considerthe localfluence distribution of FRBs andfind that it can constrain the subset of FRB-producing magnetarprogenitors. Our model suggests that improvements in sensitivity may reveal aflattening of the global FRBfluencedistribution and saturation in FRB rates.

Cosmology↗
NASA NTRS1993

Linearization of Attitude-Control Error Dynamics

Direction cosines and quaternions are useful for representing rigid-body attitude because they exhibit no kinematic singularities. Each utilizes more variables than the minimum three degrees of freedom required to specify attitude. Therefore, application of a nonlinear inversion procedure to either formulation introduces singularities. Furthermore, in designing an attitude-control system, it is not appropriate to express attitude error as a difference of direction cosines (or quaternions). One should employ a measure of attitude error that not only is minimal but preserves orthogonal rotation properties as well. This note applies an inversion procedure to an appropriate measure of attitude error, so that the singularity occurs when the error reaches +/- 180 deg. This approach leads to the realization of a new model-follower attitude-control system that exhibits exact linear attitude-error dynamics.

Bach, Ralph↗
DOE OSTI2022

Effect of differently oriented interlayer phases on the radiation damage of Inconel-$\mathrm{Ni}$ multimetallic layered composite

Multimetallic layered composites (MMLCs) have shown an excellent potential for application under extreme environments, e.g., accident-tolerant fuel cladding, because of their low oxidation tendency and high corrosion resistance. Interfacial phases or complexions in nanocrystalline materials accelerate the annihilation of defects and enhance the radiation resistance of materials, making MMLCs with engineered interlayer phases compelling to deploy in extreme conditions. However, implementation of MMLCs in full capacity remained a challenge due to a lack of fundamental understanding of the underlying mechanisms governing the characteristics of the interface between the metallic layers. The precise role of interlayer phases in MMLCs and their interaction with defects, specifically under extreme conditions, is still unexplored. Pursuing atomistic simulations for various Inconel-Ni MMLCs model materials, we revealed accelerated defect mobility in interlayers with larger crystalline misorientation and the inverse relationship between the interface sink strength to the misorientation angle. Furthermore, we found a linear relation between interlayer misorientation angle with the density of radiation-induced defects and radiation enhanced displacements. So, our results indicate that radiation-induced material degradation is accelerated by the higher defect formation tendency of MMLCs with a high-angle interlayer interface.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Idealized Large-Eddy Simulations of Stratocumulus Advecting over Cold Water. Part I: Boundary Layer Decoupling

Abstract We explore the decoupling physics of a stratocumulus-topped boundary layer (STBL) moving over cooler water, a situation mimicking warm-air advection (WADV). We simulate an initially well-mixed STBL over a doubly periodic domain with the sea surface temperature decreasing linearly over time using the System for Atmospheric Modeling large-eddy model. Due to the surface cooling, the STBL becomes increasingly stably stratified, manifested as a near-surface temperature inversion topped by a well-mixed cloud-containing layer. Unlike the stably stratified STBL in cold-air advection (CADV) that is characterized by cumulus coupling, the stratocumulus deck in the WADV is unambiguously decoupled from the sea surface, manifested as weakly negative buoyancy flux throughout the subcloud layer. Without the influxes of buoyancy from the surface, the convective circulation in the well-mixed cloud-containing layer is driven by cloud-top radiative cooling. In such a regime, the downdrafts propel the circulation, in contrast to that in CADV regime for which the cumulus updrafts play a more determinant role. Such a contrast in convection regime explains the difference in many aspects of the STBLs including the entrainment rate, cloud homogeneity, vertical exchanges of heat and moisture, and lifetime of the stratocumulus deck, with the last being subject to a more thorough investigation in Part II. Finally, we investigate under what conditions a secondary stratus near the surface (or fog) can form in the WADV. We found that weaker subsidence favors the formation of fog whereas a more rapid surface cooling rate does not. Significant Statement The low-lying blanketlike clouds, called stratocumulus (Sc), reflect much incoming sunlight, substantially modulating Earth’s temperature. While much is known about how the Sc evolves when it moves over warmer water, few studies examine the opposite situation of Sc moving over colder water. We used a high-resolution numerical model to simulate such a case. When moving over cold water, the Sc becomes unambiguously decoupled from the water surface, distinctive from its warm counterpart in which the Sc interacts with the water surface via intermittent cauliflower-like clouds called cumulus clouds. Such decoupling influences many aspects of the Sc–sea surface system, which combine to alter the ability of the Sc to reflect sunlight, thereby influencing the climate. This work laid the foundation for future work that quantifies the contribution of such a decoupled Sc regime to Earth’s radiative budget and climate change.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2012

A Parameterized Inversion Model for Soil Moisture and Biomass from Polarimetric Backscattering Coefficients

A semi-empirical algorithm for the retrieval of soil moisture, root mean square (RMS) height and biomass from polarimetric SAR data is explained and analyzed in this paper. The algorithm is a simplification of the distorted Born model. It takes into account the physical scattering phenomenon and has three major components: volume, double-bounce and surface. This simplified model uses the three backscattering coefficients ( sigma HH, sigma HV and sigma vv) at low-frequency (P-band). The inversion process uses the Levenberg-Marquardt non-linear least-squares method to estimate the structural parameters. The estimation process is entirely explained in this paper, from initialization of the unknowns to retrievals. A sensitivity analysis is also done where the initial values in the inversion process are varying randomly. The results show that the inversion process is not really sensitive to initial values and a major part of the retrievals has a root-mean-square error lower than 5% for soil moisture, 24 Mg/ha for biomass and 0.49 cm for roughness, considering a soil moisture of 40%, roughness equal to 3cm and biomass varying from 0 to 500 Mg/ha with a mean of 161 Mg/ha

soil moisture↗
NASA NTRS2001

Our Sun V: A Bright Young Sun Consistent with Helioseismology and Warm Temperatures on Ancient Earth and Mars

The relatively warm temperatures required on early Earth and Mars have been difficult to account for with warming from greenhouse gases. A slightly more massive young Sun would be brighter than predicted by the standard solar model, simultaneously resolving this problem for both Earth and Mars. We computed high-precision solar models with seven initial masses, from Mi = 1.01 to 1.07 solar mass - the latter being the maximum permitted if the early Earth is not to lose its water via a moist greenhouse effect. The relatively modest early mass loss that is required remains consistent with observational limits on mass loss from young stars and with estimates of the past solar wind obtained from lunar rocks. We considered three types of mass loss rates: (1) a reasonable choice of a simple exponential decline, (2) an extreme step-function case that gives the maximum effect consistent with observations, and (3) the radical case of a linear decline which is inconsistent with the solar wind mass loss estimates from lunar rocks. Our computations demonstrated that mass loss leaves a fingerprint oil the Sun's internal structure large enough to be detectable with helioseismic observations. All of our mass-losing solar models were consistent with the helioseismic observations; in fact, our preferred mass-losing cases were in marginally better agreement with the helioseismology than the standard solar model was, although this difference was smaller than the effects of other uncertainties in the input physics and in the solar composition. Mass loss has only a relatively minor effect on the predicted lithium depletion; the major portion of the solar lithium depletion must still be due to rotational mixing. Thus the modest mass loss cases considered here cannot be ruled out by observed lithium depletions. For the three mass loss types considered, the preferred initial masses were 1.07 solar mass for the exponential case and 1.04 solar mass for the step-function and linear cases; all of these provided high enough solar fluxes at Mars 3.8 Gyr ago to be consistent with the existence of liquid water. For a more massive early Sun, the planets would have had to be closer to the young Sun in order to end up in their present orbits; the orbital radii of the planets would vary inversely with the solar mass. Both of these effects contribute to the fact that the early solar flux at the planets would have been considerably higher than that of the standard solar model at that time. In fact, the 1.07 solar mass exponential case has a flux at birth 5% higher than the present solar flux, while the radical 1.04 solar mass linear case has a nearly constant flux over the first 3 Gyr only about 10% lower than at present. The early solar evolution would be in the opposite direction in the H-R diagram to that of the standard Sun.

Sackmann, I.-Juliana↗
DOE OSTI2024

Simulations and models for the Richtmyer–Meshkov instability with broadband perturbations

The Richtmyer–Meshkov instability (RMI) is shock driven and affects many phenomena from inertial fusion to supernova explosions. The behavior of single-modes in the RMI has been studied extensively but less is known with the broadband perturbations that occur in applications. Here, we describe extensive numerical simulations and modeling of the RMI with broadband perturbations with an initial power spectrum of the form P ∝ km, where k is the wavenumber. The hydrodynamic condition is the same as in the θ-Group Collaboration [Thornber et al., “Late-time growth rate, mixing, and anisotropy in the multimode narrowband Richtmyer–Meshkov instability: The θ-group collaboration,” Phys. Fluids 29, 105107 (2017)] with a Mach ∼1.86 shock and gamma-law = 5/3 fluids with Atwood number A = −0.5. The bubble amplitude hB is found to grow in two stages. Initially, hB undergoes a phase-inversion and grows linearly in time (t) at a rate consistent with a linear theory. Asymptotically, hB grows as a power law ∼tθ when k|hB| > O(1), where k is the average wavenumber for the initial spectrum. The RMI behavior in simulations and modeling agree over a wide range of exponent m, spectral width, initial amplitude, and time. The agreement is quantified objectively using statistical analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Magnetic excitations of the hybrid multiferroic (ND 4 ) 2 FeCl 5 ·D 2 O

We report a comprehensive inelastic neutron scattering study of the hybrid molecule-based multiferroic compound ( ND 4 ) 2 FeCl 5 · D 2 O in the zero-field incommensurate cycloidal phase and the high-field quasicollinear phase. The spontaneous electric polarization changes its direction concurrently with the field-induced magnetic transition, from mostly aligned with the crystallographic a axis to the c axis. To account for such a change in polarization direction, the underlying multiferroic mechanism was proposed to switch from the spin-current model induced via the inverse Dzyaloshinskii-Moriya interaction to the p - d hybridization model. We perform a detailed analysis of the inelastic neutron data of ( ND 4 ) 2 FeCl 5 · D 2 O using linear spin-wave theory to quantify magnetic interaction strengths and investigate the possible impact of different multiferroic mechanisms on the magnetic couplings. Overall our result reveals that the spin dynamics of both multiferroic phases can be well described by a Heisenberg Hamiltonian with easy-plane anisotropy. We do not find notable differences between the optimal model parameters of the two phases. The hierarchy of exchange couplings and the balance among frustrated interactions remain the same between two phases, suggesting that magnetic interactions in ( ND 4 ) 2 FeCl 5 · D 2 O are much more robust than the electric polarization in response to delicate reorganizations of the electronic degrees of freedom in an applied magnetic field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Autonomous Inversion of In Situ Deformation Measurement Data for Injection-Induced Stress Change

Geologic carbon storage (GCS) is likely to play a key part of the global effort to dramatically reduce CO2 emissions and perhaps even reduce atmospheric CO2 concentrations through carbon negative operations. A critical part of effort to commercialize and widely deploy this technology is developing the capability to rapidly assimilate real-time monitoring data into a form that will enable site operators to make decisions to manage the safe and efficient operations. Two of the risks associate with GCS are the risk of inducing fractures in the sealing formations that can create leakage pathways and the risk of inducing earthquakes of sufficient magnitude to cause public concern, property damage, or safety risks. To properly manage these risks the site operator needs to know the initial state of stress, the change in stress induced by injection, and the relationship between operational parameters such as injection rate and pressure and the change in stress. Current methods of estimating the change in stress require choosing the type of constitutive model and the model parameters based on core, log, and geophysical data during the characterization phase, with little feedback from operational observations to validate or refine these choices. These characterization methods interrogate the geologic formations using length scales, loading rates or magnitudes that are quite different from those encountered by the actual storage system. It is shown that errors in the assumed constitutive response, even when informed by laboratory tests on core samples, are likely to be common, large, and underestimate the magnitude of stress change caused by injection. Recent advances in borehole-based strain instruments and borehole and surface-based tilt and displacement instruments have now enabled monitoring of the deformation of the storage system throughout its operational lifespan. This data can enable validation and refinement of the knowledge of the geomechanical properties and state of the system, but brings with it a challenge to transform the raw data into actionable knowledge. We demonstrate a method that uses automatic differentiation and a finite-element based geomechanical model perform a gradient-based deterministic inversion of geomechanical monitoring data. This approach allows autonomous integration of the instrument data without the need for time consuming manual interpretation and selection of updated model parameters. Furthermore, only isotropic linear elasticity is considered in this paper, the approach presented is very flexible as to what type of geomechanical constitutive response can be used. The approach is easily adaptable to nonlinear physics-based constitutive models to account for common rock behaviors such as creep and plasticity. The approach also enables training of machine learning-based constitutive models by allowing back propagation of errors through the finite element calculations. This enables strongly enforcing known physics, such as conservation of momentum and continuity, while allowing data-driven models to learn the truly unknown physics such as the constitutive or petrophysical responses.

Burghardt, Jeffrey A.↗
NASA NTRS1984

Optimization and mechanisms of mistuning in cascades

In the present inverse design procedure for the optimum mistuning of a high bypass ratio shroudless fan that is modeled as a cascade of blades (each with a single torsional degree-of-freedom), linearized supersonic aerodynamic theory is used to compute the unsteady aerodynamic forces in the 'influence coefficient' form at a typical blade section. The mistuning pattern is then numerically optimized in order to achieve a specified increase in the aeroelastic stability margin with a minimum amount of mistuning. If the blades are self-damped, an optimized mistuning pattern can be found that achieves a given stability margin for a much lower level of mistuning than required for the alternate mistuning pattern, which requires only two blade frequencies and is relatively insensitive to implementation errors.

Crawley, E. F.↗
NASA NTRS1995

In situ phytoplankton absorption, fluorescence emission, and particulate backscattering spectra determined from reflectance

An inverse model was developed to extract the absortion and scattering (elastic and inelastic) properties of oceanic constituents from surface spectral reflectance measurements. In particular, phytoplankton spectral absorption coefficients, solar-stimulated chlorophyll a fluorescence spectra, and particle backscattering spectra were modeled. The model was tested on 35 reflectance spectra obtained from irradiance measurements in optically diverse ocean waters (0.07 to 25.35 mg/cu m range in surface chlorophyll a concentrations). The universality of the model was demonstrated by the accurate estimation of the spectral phytoplankton absorption coefficents over a range of 3 orders of magnitude (rho = 0.94 at 500 nm). Under most oceanic conditions (chlorophyll a less than 3 mg/cu m) the percent difference between measured and modeled phytoplankton absorption coefficents was less than 35%. Spectral variations in measured phytoplankton absorption spectra were well predicted by the inverse model. Modeled volume fluorescence was weakly correlated with measured chl a; fluorescence quantum yield varied from 0.008 to 0.09 as a function of environment and incident irradiance. Modeled particle backscattering coefficients were linearly related to total particle cross section over a twentyfold range in backscattering coefficents (rho = 0.996, n = 12).

Roesler, Collin S.↗
NASA NTRS1998

A Multipole Expansion Method for Analyzing Lightning Field Changes

Changes in the surface electric field are frequently used to infer the locations and magnitudes of lightning-caused changes in thundercloud charge distributions. The traditional procedure is to assume that the charges that are effectively deposited by the flash can be modeled either as a single point charge (the Q-model) or a point dipole (the P-model). The Q-model has 4 unknown parameters and provides a good description of many cloud-to-ground (CG) flashes. The P-model has 6 unknown parameters and describes many intracloud (IC) discharges. In this paper, we introduce a new analysis method that assumes that the change in the cloud charge can be described by a truncated multipole expansion, i.e., there are both monopole and dipole terms in the unknown source distribution, and both terms are applied simultaneously. This method can be used to analyze CG flashes that are accompanied by large changes in the cloud dipole moment and complex IC discharges. If there is enough information content in the measurements, the model can also be generalized to include quadrupole and higher order terms. The parameters of the charge moments are determined using a 3-dimensional grid search in combination with a linear inversion, and because of this, local minima in the error function and the associated solution ambiguities are avoided. The multipole method has been tested on computer simulated sources and on natural lightning at the NASA Kennedy Space Center and USAF Eastern Range.

Koshak, William J.↗
NASA NTRS1990

Sedimentary Facies Analysis Using AVIRIS Data: A Geophysical Inverse Problem

AVIRIS data can be used to quantitatively analyze and map sedimentary lithofacies. The observed radiance spectra can be reduced to 'apparent reflectance' spectra by topographic and reflectance characterization of several field sites within the image. These apparent reflectance spectra correspond to the true reflectance at each pixel, multiplied by an unknown illumination factor (ranging in value from zero to one). The spatial abundance patterns of spectrally defined lithofacies and the unknown illumination factors can be simultaneously derived using constrained linear spectral unmixing methods. Estimates of the minimum uncertainty in the final results (due to noise, instrument resolutions, degree of illumination and mixing systematics) can be made by forward and inverse modeling. Specific facies studies in the Rattlesnake Hills region of Wyoming illustrate the successful application of these methods.

Boardmann, Joe W.↗