Search NASASearch

SEARCH · Search NASA

Results for “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 91 records · Page 5

An Analytical Calibration Approach for the Polarimetric Airborne C Band Radiometer

Passive microwave remote sensing is sensitive to the quantity and distribution of water in soil and vegetation. During summer 2000, the Microwave Geophysics Group a t the University of Michigan conducted the seventh Radiobrighness Energy Balance Experiment (REBEX-7) over a corn canopy in Michigan. Long time series of brightness temperatures, soil moisture and micrometeorology on the plot were taken. This paper addresses the calibration of the NASA GSFC polarimetric airborne C band microwave radiometer (ACMR) that participated in REBEX-7. These passive polarimeters are typically calibrated using an end-to-end approach based upon a standard artificial target or a well-known geophysical target. Analyzing the major internal functional subsystems offers a different perspective. The primary goal of this approach is to provide a transfer function that not only describes the system in its entire5 but also accounts for the contributions of each subsystem toward the final modified Stokes parameters. This approach does not assume that the radiometric system is linear as it does not take polarization isolation for granted, and it also serves as a realistic instrument simulator, a useful tool for future designs. The ACMR architecture can be partitioned into functional subsystems. The characteristics of each subsystem was extensively measured and the estimated parameters were imported into the overall dosed form system model. Inversion of the model yields a calibration for the modeled Stokes parameters with uncertainties of 0.2 K for the V and H polarizations and 2.4 K for the 3rd and 4th parameters. Application to the full Stokes parameters over a senescent cornfield is presented.

Pham, Hanh

An Atmospheric Science Observing System Simulation Experiment (OSSE) Environment

An atmospheric sounding mission starts with a wide range of concept designs involving measurement technologies, observing platforms, and observation scenarios. Observing system simulation experiment (OSSE) is a technical approach to evaluate the relative merits of mission and instrument concepts. At Jet Propulsion Laboratory (JPL), the OSSE team has developed an OSSE environment that allows atmospheric scientists to systematically explore a wide range of mission and instrument concepts and formulate a science traceability matrix with a quantitative science impact analysis. The OSSE environment virtually creates a multi-platform atmospheric sounding testbed (MAST) by integrating atmospheric phenomena models, forward modeling methods, and inverse modeling methods. The MAST performs OSSEs in four loosely coupled processes, observation scenario exploration, measurement quality exploration, measurement quality evaluation, and science impact analysis.

mission concepts

ISSM: Ice Sheet System Model

In order to have the capability to use satellite data from its own missions to inform future sea-level rise projections, JPL needed a full-fledged ice-sheet/iceshelf flow model, capable of modeling the mass balance of Antarctica and Greenland into the near future. ISSM was developed with such a goal in mind, as a massively parallelized, multi-purpose finite-element framework dedicated to ice-sheet modeling. ISSM features unstructured meshes (Tria in 2D, and Penta in 3D) along with corresponding finite elements for both types of meshes. Each finite element can carry out diagnostic, prognostic, transient, thermal 3D, surface, and bed slope simulations. Anisotropic meshing enables adaptation of meshes to a certain metric, and the 2D Shelfy-Stream, 3D Blatter/Pattyn, and 3D Full-Stokes formulations capture the bulk of the ice-flow physics. These elements can be coupled together, based on the Arlequin method, so that on a large scale model such as Antarctica, each type of finite element is used in the most efficient manner. For each finite element referenced above, ISSM implements an adjoint. This adjoint can be used to carry out model inversions of unknown model parameters, typically ice rheology and basal drag at the ice/bedrock interface, using a metric such as the observed InSAR surface velocity. This data assimilation capability is crucial to allow spinning up of ice flow models using available satellite data. ISSM relies on the PETSc library for its vectors, matrices, and solvers. This allows ISSM to run efficiently on any parallel platform, whether shared or distrib- ISSM: Ice Sheet System Model NASA's Jet Propulsion Laboratory, Pasadena, California uted. It can run on the largest clusters, and is fully scalable. This allows ISSM to tackle models the size of continents. ISSM is embedded into MATLAB and Python, both open scientific platforms. This improves its outreach within the science community. It is entirely written in C/C++, which gives it flexibility in its design, and the power/speed that C/C++ allows. ISSM is svn (subversion) hosted, on a JPL repository, to facilitate its development and maintenance. ISSM can also model propagation of rifts using contact mechanics and mesh splitting, and can interface to the Dakota software. To carry out sensitivity analysis, mesh partitioning algorithms are available, based on the Scotch, Chaco, and Metis partitioners that ensure equal area mesh partitions can be done, which are then usable for sampling and local reliability methods.

Larour, Eric

Research In Nonlinear Flight Control for Tiltrotor Aircraft Operating in the Terminal Area

The research during the first year of the effort focused on the implementation of the recently developed combination of neural net work adaptive control and feedback linearization. At the core of this research is the comprehensive simulation code Generic Tiltrotor Simulator (GTRS) of the XV-15 tilt rotor aircraft. For this research the GTRS code has been ported to a Fortran environment for use on PC. The emphasis of the research is on terminal area approach procedures, including conversion from aircraft to helicopter configuration. This report focuses on the longitudinal control which is the more challenging case for augmentation. Therefore, an attitude command attitude hold (ACAH) control augmentation is considered which is typically used for the pitch channel during approach procedures. To evaluate the performance of the neural network adaptive control architecture it was necessary to develop a set of low order pilot models capable of performing such tasks as, follow desired altitude profiles, follow desired speed profiles, operate on both sides of powercurve, convert, including flaps as well as mastangle changes, operate with different stability and control augmentation system (SCAS) modes. The pilot models are divided in two sets, one for the backside of the powercurve and one for the frontside. These two sets are linearly blended with speed. The mastangle is also scheduled with speed. Different aspects of the proposed architecture for the neural network (NNW) augmented model inversion were also demonstrated. The demonstration involved implementation of a NNW architecture using linearized models from GTRS, including rotor states, to represent the XV-15 at various operating points. The dynamics used for the model inversion were based on the XV-15 operating at 30 Kts, with residualized rotor dynamics, and not including cross coupling between translational and rotational states. The neural network demonstrated ACAH control under various circumstances. Future efforts will include the implementation into the Fortran environment of GTRS, including pilot modeling and NNW augmentation for the lateral channels. These efforts should lead to the development of architectures that will provide for fully automated approach, using similar strategies.

Calise, A. J.

A remanent and induced magnetization model of Magsat vector anomalies over the west African craton

Scalar and vector Magsat anomalies over the west African craton are analyzed by forward and inverse models. A forward model of the Man shield is based on Liberia. Induced magnetization contrasts due to sporadic iron-formations and to regional metamorphic rocks, and a contrast in remanent magnetization within the lower crust are included. This combination reproduces the location, magnitude and adopted local zero level of anomalies in the initial Magsat maps. An inverse model of the Reguibat shield estimates the magnetization contrast of its lithosphere, and when magnetism is restricted to shallower than 75 km both shields can be represented by a susceptibility contrast of +0.02. A residual anomaly between the shields involves a relative deficiency of induced magnetization along with other causes.

Toft, P. B.

Objective Operational Utilization of Satellite Microwave Scatterometer Observations of Tropical Cyclones

This study has demonstrated that high-resolution scatterometer measurements in tropical cyclones and other high-marine surface wind regimes may be retrieved accurately for wind speeds up to about 35 mls (1-hour average at 10 m) when the scatterometer data are processed through a revised geophysical model function, and a spatial adaptive algorithm is applied which utilizes the fact that wind direction is so tightly constrained in tile inner core of severe marine storms that wind direction may be prescribed from conventional data. This potential is demonstrated through case studies with NSCAT data in a severe West Pacific Typhoon (Violet, 1996) and an intense North Atlantic hurricane (Lili, 1996). However, operational scatterometer winds from NSCAT and QuickScat in hurricanes and severe winter storms are biased low in winds above 25 m/s. We have developed an inverse model to specify the entire surface wind field about a tropical cyclone from operational QuickScat scatterometer measurements within 150 nm of a storm center with the restriction that only wind speeds up to 20 m/s are used until improved model function are introduced. The inverse model is used to specify the wind field over the entire life-cycle of Hurricane Floyd (1999) for use to drive an ocean wave model. The wind field compares very favorably with wind fields developed from the copious aircraft flight level winds obtained in this storm.

Cardone, Vincent J.

Modeling and Inverse Controller Design for an Unmanned Aerial Vehicle Based on the Self-Organizing Map

The next generation of aircraft will have dynamics that vary considerably over the operating regime. A single controller will have difficulty to meet the design specifications. In this paper, a SOM-based local linear modeling scheme of an unmanned aerial vehicle (UAV) is developed to design a set of inverse controllers. The SOM selects the operating regime depending only on the embedded output space information and avoids normalization of the input data. Each local linear model is associated with a linear controller, which is easy to design. Switching of the controllers is done synchronously with the active local linear model that tracks the different operating conditions. The proposed multiple modeling and control strategy has been successfully tested in a simulator that models the LoFLYTE UAV.

Cho, Jeongho

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.

Thermal radiation from Titan's atmosphere

The temperature-inversion model of Titan's atmosphere is modified to include the thermal emission from a plausible amount of C2H2. Infrared photometry of Titan from 8 to 14 microns agrees with model calculations for 0.5 cm-atm C2H6 and 1.0 cm-atm C2H2 and an optically thin dust layer in a temperature-inversion region at 160 K. The random band approximation used to calculate C2H6 opacity has been verified by laboratory transmission measurements. Measurements of the C2H2 emission at 13.7 microns may provide an indirect check on the presence of H2 on Titan. Existing observations of Titan appear to be fully consistent with the inversion model if the surface radius equals 2700 km, the atmospheric radius equals 2900 km, and the surface temperature equals 78 K. The energy balance depends sensitively on these three parameters.

Caldwell, J.

Optimization as a Tool for Consistency Maintenance in Multi-Resolution Simulation

The need for new approaches to the consistent simulation of related phenomena at multiple levels of resolution is great. While many fields of application would benefit from a complete and approachable solution to this problem, such solutions have proven extremely difficult. We present a multi-resolution simulation methodology that uses numerical optimization as a tool for maintaining external consistency between models of the same phenomena operating at different levels of temporal and/or spatial resolution. Our approach follows from previous work in the disparate fields of inverse modeling and spacetime constraint-based animation. As a case study, our methodology is applied to two environmental models of forest canopy processes that make overlapping predictions under unique sets of operating assumptions, and which execute at different temporal resolutions. Experimental results are presented and future directions are addressed.

TEMPORAL RESOLUTION

Top-down Estimate of Dust Emissions Through Integration of MODIS and MISR Aerosol Retrievals With the Geos-chem Adjoint Model

Predicting the influences of dust on atmospheric composition, climate, and human health requires accurate knowledge of dust emissions, but large uncertainties persist in quantifying mineral sources. This study presents a new method for combined use of satellite-measured radiances and inverse modeling to spatially constrain the amount and location of dust emissions. The technique is illustrated with a case study in May 2008; the dust emissions in Taklimakan and Gobi deserts are spatially optimized using the GEOSChem chemical transport model and its adjoint constrained by aerosol optical depth (AOD) that are derived over the downwind dark-surface region in China from MODIS (Moderate Resolution Imaging Spectroradiometer) reflectance with the aerosol single scattering properties consistent with GEOS-chem. The adjoint inverse modeling yields an overall 51% decrease in prior dust emissions estimated by GEOS-Chem over the Taklimakan-Gobi area, with more significant reductions south of the Gobi Desert. The model simulation with optimized dust emissions shows much better agreement with independent observations from MISR (Multi-angle Imaging SpectroRadiometer) AOD and MODIS Deep Blue AOD over the dust source region and surface PM10 concentrations. The technique of this study can be applied to global multi-sensor remote sensing data for constraining dust emissions at various temporal and spatial scales, and hence improving the quantification of dust effects on climate, air quality, and human health.

MODIS aerosol retrievals

Irrigation Requirement Estimation using MODIS Vegetation Indices and Inverse Biophysical Modeling; A Case Study for Oran, Algeria

Human demand for food influences the water cycle through diversion and extraction of fresh water needed to support agriculture. Future population growth and economic development alone will substantially increase water demand and much of it for agricultural uses. For many semi-arid lands, socio-economic shifts are likely to exacerbate changes in climate as a driver of future water supply and demand. For these areas in particular, where the balance between water supply and demand is fragile, variations in regional climate can have potentially predictable effect on agricultural production. Satellite data and biophysically-based models provide a powerful method to quantify the interactions between local climate, plant growth and water resource requirements. In irrigated agricultural lands, satellite observations indicate high vegetation density while the precipitation amount indicates otherwise. This inconsistency between the observed precipitation and the observed canopy leaf density triggers the possibility that the observed high leaf density is due to an alternate source of water, irrigation. We explore an inverse process approach using observations from the Moderate Resolution Imaging Spectroradiometer (MODIS), climatological data, and the NASA's Simple Biosphere model, SiB2, to quantitatively assess water demand in a semi-arid agricultural land by constraining the carbon and water cycles modeled under both equilibrium (balance between vegetation and prevailing local climate) and nonequilibrium (water added through irrigation) conditions. We postulate that the degree to which irrigated lands vary from equilibrium conditions is related to the amount of irrigation water used. We added water using two distribution methods: The first method adds water on top of the canopy and is a proxy for the traditional spray irrigation. The second method allows water to be applied directly into the soil layer and serves as proxy for drip irrigation. Our approach indicates that over the study site, for the month of July, spray irrigation resulted in an irrigation amount of about 1.4 mm per occurrence with an average frequency of occurrence of 24.6 hours. The simulated total monthly irrigation for July was 34.85 mm. In contrast, the drip irrigation resulted in less frequent irrigation events with an average water requirement about 57% less than that simulated during the spray irrigation case. The efficiency of the drip irrigation method rests on its reduction of the canopy interception loss compared to the spray irrigation method. When compared to a country-wide average estimate of irrigation water use, our numbers are quite low. We would have to revise the reported country level estimates downward to 17% or less

Bounoua, L.

Diffusion kinetics of 3 He in pyroxene and plagioclase and applications to cosmogenic exposure dating and paleothermometry in mafic rocks

In this study, we investigate the diffusivity of cosmogenic 3 He in a variety of plagioclase and pyroxene compositions, and its application to paleothermometry and exposure dating in these minerals, through stepwise degassing experiments. While cosmogenic 3 He has been utilized for exposure dating in pyroxene for decades due to its retentivity, plagioclase, often found along with pyroxene in mafic rocks, is generally less retentive of cosmogenic noble gas. However, the diffusivity of 3 He in either plagioclase or pyroxene has not yet been measured quantitatively. A challenge in measuring diffusion kinetics by step-degassing experiments in poorly retentive minerals is the fact that significant amounts of He can be lost prior to the experiment. To address this issue, we apply a forward “multiple diffusion domain” (MDD) inversion model that includes model predictions of initial gas loss during irradiation and storage of the samples to account for this observation and add constraints to the diffusion parameters. We find that 3 He diffusivity in plagioclase appears to be highly variable. This variability can be explained by the MDD inversion models' inability to constrain the diffusion parameters when significant gas has been lost during irradiation and/or prolonged storage prior to experiment analysis, resulting in an overestimation of 3 He retentivity. Plagioclase samples that were kept frozen after irradiation, to limit this initial gas loss, yielded the most reliable estimate of diffusion kinetics. We find that 3 He in plagioclase is diffusively lost at Earth's surface temperatures on a timescale of 100 years and is therefore unsuitable for surface temperature paleothermometry. On the contrary, we find cosmogenic 3 He in pyroxene to be retentive at Earth's surface temperatures on a 1-million-year timescale.

Bergelin, Marie [Berkeley Geochronology Center (BG

Low-Cost Evaluation of EO-1 Hyperion and ALI for Detection and Biophysical Characterization of Forest Logging in Amazonia (NCC5-481)

Major uncertainties exist regarding the rate and intensity of logging in tropical forests worldwide: these uncertainties severely limit economic, ecological, and biogeochemical analyses of these regions. Recent sawmill surveys in the Amazon region of Brazil show that the area logged is nearly equal to total area deforested annually, but conversion of survey data to forest area, forest structural damage, and biomass estimates requires multiple assumptions about logging practices. Remote sensing could provide an independent means to monitor logging activity and to estimate the biophysical consequences of this land use. Previous studies have demonstrated that the detection of logging in Amazon forests is difficult and no studies have developed either the quantitative physical basis or remote sensing approaches needed to estimate the effects of various logging regimes on forest structure. A major reason for these limitations has been a lack of sufficient, well-calibrated optical satellite data, which in turn, has impeded the development and use of physically-based, quantitative approaches for detection and structural characterization of forest logging regimes. We propose to use data from the EO-1 Hyperion imaging spectrometer to greatly increase our ability to estimate the presence and structural attributes of selective logging in the Amazon Basin. Our approach is based on four "biogeophysical indicators" not yet derived simultaneously from any satellite sensor: 1) green canopy leaf area index; 2) degree of shadowing; 3) presence of exposed soil and; 4) non-photosynthetic vegetation material. Airborne, field and modeling studies have shown that the optical reflectance continuum (400-2500 nm) contains sufficient information to derive estimates of each of these indicators. Our ongoing studies in the eastern Amazon basin also suggest that these four indicators are sensitive to logging intensity. Satellite-based estimates of these indicators should provide a means to quantify both the presence and degree of structural disturbance caused by various logging regimes. Our quantitative assessment of Hyperion hyperspectral and ALI multi-spectral data for the detection and structural characterization of selective logging in Amazonia will benefit from data collected through an ongoing project run by the Tropical Forest Foundation, within which we have developed a study of the canopy and landscape biophysics of conventional and reduced-impact logging. We will add to our base of forest structural information in concert with an EO-1 overpass. Using a photon transport model inversion technique that accounts for non-linear mixing of the four biogeophysical indicators, we will estimate these parameters across a gradient of selective logging intensity provided by conventional and reduced impact logging sites. We will also compare our physical ly-based approach to both conventional (e.g., NDVI) and novel (e.g., SWIR-channel) vegetation indices as well as to linear mixture modeling methods. We will cross-compare these approaches using Hyperion and ALI imagers to determine the strengths and limitations of these two sensors for applications of forest biophysics. This effort will yield the first physical ly-based, quantitative analysis of the detection and intensity of selective logging in Amazonia, comparing hyperspectral and improved multi-spectral approaches as well as inverse modeling, linear mixture modeling, and vegetation index techniques.

Asner, Gregory P.

Polar Motion Constraints on Models of the Fortnightly Tide

Estimates of the near-fortnightly Mf ocean tide from Topex/Poseidon satellite altimetry and from numerical solutions to the shallow water equations agree reasonably well, at least in their basin-scale features. For example, both show that the Pacific Ocean tide lags the Atlantic tide by roughly 30 degrees. There are hints of finer scale agreements in the elevation fields, but noise levels are high. In contrast, estimates of Mf currents are only weakly constrained by the TP data, because high-wavenumber Rossby waves (with intense currents) are associated with relatively small perturbations in surface elevation. As a result, a wide range of Mf current fields are consistent with both the TP data and the hydrodynamic equations within a priori plausible misfit bounds. We find that a useful constraint on the Mf currents is provided by independent estimates of the Earth's polar motion. At the Mf period polar motion shows a weak signal (both prograde and retrograde) which must be almost entirely caused by the ocean tide. We have estimated this signal from the SPACE2000 time series, after applying a broad-band correction for atmospheric angular momentum. Although the polar motion estimates have relatively large uncertainties, they are sufficiently precise to fix optimum data weights in a global ocean inverse model of Mf. These weights control the tradeoff between fitting a prior hydrodynamic model of Mf and fitting the relatively noisy T/P measurements of Mf. The predicted polar motion from the final inverse model agrees remarkably well with the Mf polar motion observations. The preferred model is also consistent with noise levels suggested by island gauges, and it is marginally consistent with differences observed by subsetting the altimetry (to the small extent that this is possible). In turn, this new model of the Mf ocean tide allows the ocean component to be removed from Mf estimates of length of day, thus yielding estimates of complex Love numbers less contaminated by oceanic effects than has hitherto been possible.

Ray, Richard D.

An explicit canopy BRDF model and inversion

Based on a rigorous canopy radiative transfer equation, the multiple scattering radiance is approximated by the asymptotic theory, and the single scattering radiance calculation, which requires an numerical intergration due to considering the hotspot effect, is simplified. A new formulation is presented to obtain more exact angular dependence of the sky radiance distribution. The unscattered solar radiance and single scattering radiance are calculated exactly, and the multiple scattering is approximated by the delta two-stream atmospheric radiative transfer model. The numerical algorithms prove that the parametric canopy model is very accurate, especially when the viewing angles are smaller than 55 deg. The Powell algorithm is used to retrieve biospheric parameters from the ground measured multiangle observations.

Liang, Shunlin

Characteristics of Total Column CO2 Retrievals from the OCO Missions: Biases, Information Content and Implications for Flux Inversions

The Orbiting Carbon Observatory-2 and Orbiting Carbon Observatory-3, launched in 2015 and 2019, respectively, are intended to collect and deliver high-resolution observations of CO2 with unprecedented space and time coverage. Observations of CO2 from these remote-sensing missions (also known as XCO2, or column-based average, dry air mole fraction of CO2) are then used by the global carbon cycle community to answer a wide range of science questions, from the distribution and quantification of global and regional CO2 source-sink patterns to quantification of anthropogenic sources at urban scales. Even though we have had the OCO-2 mission flying for a few years now, the retrieval algorithms are continuously evolving and improving to deliver XCO2 retrievals with very high precision and high accuracy (or low biases). In this presentation, we will discuss a simple yet effective quantitative framework that has been developed by the OCO-2 flux team to evaluate the information content of these XCO2 retrievals as soon as they are released, i.e., with lower latency than full-scale flux inversions. This framework serves as a precursor to advanced inverse modeling frameworks and is intended to provide an early but accurate assessment of the signal present in the satellite retrievals, the robustness of that signal, and the ability of these retrievals to resolve patterns in CO2 surface fluxes that cannot be resolved by our current network of surface sites. Specific results will tackle a tiered set of questions that are being addressed using this framework: (a) what are the distribution of retrievals in the different modes of operation and how do they vary in space and time? (b) what is the information that is being given to the inverse modeling frameworks from the space-based data, information above and beyond what is provided by the in-situ data? and (c) how do these factors influence our choices for doing flux inversions with the satellite retrievals? While the primary focus of the results will be on application of this technique to mature OCO-2 retrievals, we will show early results for a couple of months of OCO-3 retrievals. For the time-period that the retrievals from the two missions overlap, we will highlight how this framework allows us to effortlessly put the information from OCO-3 and OCO-2 on an equal footing, thus enabling easy comparison between the two pioneering missions.

Chatterjee, Abhishek