Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sparse estimation”

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 55 records · Page 3

GES DISC Long-Term Data Analysis Services Using AIRS and CMS Methane Data as an Example

NASA Goddard Earth Sciences Data and Information Services Center (GES DISC) archives and distributes rich collections of data on atmospheric greenhouse gases from multiple satellite missions and model results. Among those greenhouse gases, atmospheric methane is a powerful greenhouse gas contributing ~0.5 (W/m^2) to total radiative forcing, and its concentration has increased by about 150% since 1750. Observations or estimates of methane emissions typically have sparse spatial and temporal coverage. The lack of comprehensive spatial and temporal coverage of methane source and sink observations has made analyzing atmospheric methane trends challenging. In this study the GES DISC aims to provide the community with the resources to better understand changes in atmospheric methane concentrations and the underlying causes. We will utilize methane datasets from Atmospheric Infrared Sounder (AIRS) retrieved methane concentration and three Carbon Monitoring System (CMS) methane emission datasets (in regions of North America, Canada, and Mexico) to compare AIRS methane growth with corresponding CMS regional methane emissions. Comparisons of AIRS methane growth rates and CMS methane emissions suggests wetland emissions may impact methane growth rate trends over North America. As the record for CMS methane data is extended, both datasets can be used in conjunction to better understand impacts on atmospheric methane trends. GES DISCs new anomaly tool can also be used on select datasets to further quantify trends in atmospheric greenhouse gases.

long-term data analysis↗

Application of Statistical Methods of Rain Rate Estimation to Data From The TRMM Precipitation Radar

The TRMM Precipitation Radar is well suited to statistical methods in that the measurements over any given region are sparsely sampled in time. Moreover, the instantaneous rain rate estimates are often of limited accuracy at high rain rates because of attenuation effects and at light rain rates because of receiver sensitivity. For the estimation of the time-averaged rain characteristics over an area both errors are relevant. By enlarging the space-time region over which the data are collected, the sampling error can be reduced. However. the bias and distortion of the estimated rain distribution generally will remain if estimates at the high and low rain rates are not corrected. In this paper we use the TRMM PR data to investigate the behavior of 2 statistical methods the purpose of which is to estimate the rain rate over large space-time domains. Examination of large-scale rain characteristics provides a useful starting point. The high correlation between the mean and standard deviation of rain rate implies that the conditional distribution of this quantity can be approximated by a one-parameter distribution. This property is used to explore the behavior of the area-time-integral (ATI) methods where fractional area above a threshold is related to the mean rain rate. In the usual application of the ATI method a correlation is established between these quantities. However, if a particular form of the rain rate distribution is assumed and if the ratio of the mean to standard deviation is known, then not only the mean but the full distribution can be extracted from a measurement of fractional area above a threshold. The second method is an extension of this idea where the distribution is estimated from data over a range of rain rates chosen in an intermediate range where the effects of attenuation and poor sensitivity can be neglected. The advantage of estimating the distribution itself rather than the mean value is that it yields the fraction of rain contributed by the light and heavy rain rates. This is useful in estimating the fraction of rainfall contributed by the rain rates that go undetected by the radar. The results at high rain rates provide a cross-check on the usual attenuation correction methods that are applied at the highest resolution of the instrument.

Meneghini, R.↗

An Indirect System Identification Technique for Stable Estimation of Continuous-Time Parameters of the Vestibulo-Ocular Reflex (VOR)

The vestibulo-ocular reflex (VOR) is a well-known dual mode bifurcating system that consists of slow and fast modes associated with nystagmus and saccade, respectively. Estimation of continuous-time parameters of nystagmus and saccade models are known to be sensitive to estimation methodology, noise and sampling rate. The stable and accurate estimation of these parameters are critical for accurate disease modelling, clinical diagnosis, robotic control strategies, mission planning for space exploration and pilot safety, etc. This paper presents a novel indirect system identification method for the estimation of continuous-time parameters of VOR employing standardised least-squares with dual sampling rates in a sparse structure. This approach permits the stable and simultaneous estimation of both nystagmus and saccade data. The efficacy of this approach is demonstrated via simulation of a continuous-time model of VOR with typical parameters found in clinical studies and in the presence of output additive noise.

Kukreja, Sunil L.↗

Regional and seasonal estimates of fractional storm coverage based on station precipitation observations

Simulated climates using numerical atmospheric general circulation models (GCMs) have been shown to be highly sensitive to the fraction of GCM grid area assumed to be wetted during rain events. The model hydrologic cycle and land-surface water and energy balance are influenced by the parameter bar-kappa, which is the dimensionless fractional wetted area for GCM grids. Hourly precipitation records for over 1700 precipitation stations within the contiguous United States are used to obtain observation-based estimates of fractional wetting that exhibit regional and seasonal variations. The spatial parameter bar-kappa is estimated from the temporal raingauge data using conditional probability relations. Monthly bar-kappa values are estimated for rectangular grid areas over the contiguous United States as defined by the Goddard Institute for Space Studies 4 deg x 5 deg GCM. A bias in the estimates is evident due to the unavoidably sparse raingauge network density, which causes some storms to go undetected by the network. This bias is corrected by deriving the probability of a storm escaping detection by the network. A Monte Carlo simulation study is also conducted that consists of synthetically generated storm arrivals over an artificial grid area. It is used to confirm the bar-kappa estimation procedure and to test the nature of the bias and its correction. These monthly fractional wetting estimates, based on the analysis of station precipitation data, provide an observational basis for assigning the influential parameter bar-kappa in GCM land-surface hydrology parameterizations.

Gong, Gavin↗

Estimation of Faults in DC Electrical Power System

This paper demonstrates a novel optimization-based approach to estimating fault states in a DC power system. Potential faults changing the circuit topology are included along with faulty measurements. Our approach can be considered as a relaxation of the mixed estimation problem. We develop a linear model of the circuit and pose a convex problem for estimating the faults and other hidden states. A sparse fault vector solution is computed by using 11 regularization. The solution is computed reliably and efficiently, and gives accurate diagnostics on the faults. We demonstrate a real-time implementation of the approach for an instrumented electrical power system testbed, the ADAPT testbed at NASA ARC. The estimates are computed in milliseconds on a PC. The approach performs well despite unmodeled transients and other modeling uncertainties present in the system.

Gorinevsky, Dimitry↗

3D Shape Reconstruction of Small Bodies from Sparse Features

The autonomous approach of spacecraft to a small body (comet or asteroid) relies on using all available information at each phase of the approach. This paper presents new algorithms for global shape reconstructions from sparse tracked surface points. These methods leverage estimates from earlier phases, such as rotation pole, as well as a priori knowledge, such as a genus-0 body (i.e. without boundaries or topological holes). A mapping algorithm is proposed, which performs faithful reconstructions while enforcing genus-0 output through spherical parameterization. To estimate the shape of permanently shadowed regions of the body, a symmetry reconstruction method is added to the reconstruction algorithms. This method is shown to substantially increase the reconstruction accuracy but is subject to the symmetry of the body perpendicular to the rotation pole. The proposed mapping algorithm is compared to stateof- the-practice surface reconstruction algorithms, assessing their accuracy and ability to correctly generate genus-0 shape models for 2400 datasets and three small bodies. The proposed spherical parameterization algorithm performed consistently with the state-of-the-practice while being the only algorithm to always produce genus-0 shape models.

Nesnas, Issa↗

Underlying Fundamentals of Kalman Filtering for River Network Modeling

The grand challenge of producing hydrometeorological estimates every time and everywhere has motivated the fusion of sparse observations with dense numerical models, with a particular interest on discharge in river modeling. Ensemble methods are largely preferred as they enable the estimation of error properties, but at the expense of computational load and generally with underestimations. These imperfect stochastic estimates motivate the use of correction methods, that is, error localization and inflation, although the physical justifications for their optimality are limited. The purpose of this study is to use one of the simplest forms of data assimilation when applied to river modeling and reveal the underlying mechanisms impacting its performance. Our framework based on assimilating daily averaged in situ discharge measurements to correct daily averaged runoff was tested over a 4-yr case study of two rivers in Texas. Results show that under optimal conditions of inflation and localization, discharge simulations are consistently improved such that the mean values of Nash–Sutcliffe efficiency are enhanced from211.32 to 0.55 at observed gauges and from212.24 to21.10 at validation gauges. Yet, parameters controlling the inflation and the localization have a large impact on the performance. Further investigations of these sensitivities showed that optimal inflation occurs when compensating exactly for discrepancies in the magnitude of errors while optimal localization matches the distance traveled during one assimilation window. These results may be applicable to more advanced data assimilation methods as well as for larger applications motivated by upcoming river-observing satellite missions, such as NASA’s Surface Water and Ocean Topography mission.

Streamflow↗

Intercepted photosynthetically active radiation in wheat canopies estimated by spectral reflectance

The interception of photosynthetically active radiation (PAR) was evaluated relative to greenness and normalized difference (MSS 7-5/7+5) for five planting dates of wheat for 1978-79 and 1979-80 in Phoenix. Intercepted PAR was calculated from a model driven by leaf area index and stage of growth. Linear relationships were found between greenness and normalized difference with a separate model representing growth and senescence of the crop. Normalized difference was a significantly better model and would be easier to apply than the empirically derived greenness parameter. For the leaf area growth portion of the season the model between PAR interception and normalized difference was the same over years, however, for the leaf senescence the models showed more variability due to the lack of data on measured interception in sparse canopies. Normalized difference could be used to estimate PAR interception directly for crop growth models.

Hatfield, J. L.↗

Influence of Persistent Wind Scour on the Surface Mass Balance of Antarctica

Accurate quantification of surface snow accumulation over Antarctica is a key constraint for estimates of the Antarctic mass balance, as well as climatic interpretations of ice-core records. Over Antarctica, near-surface winds accelerate down relatively steep surface slopes, eroding and sublimating the snow. This wind scour results in numerous localized regions (< or = 200 sq km) with reduced surface accumulation. Estimates of Antarctic surface mass balance rely on sparse point measurements or coarse atmospheric models that do not capture these local processes, and overestimate the net mass input in wind-scour zones. Here we combine airborne radar observations of unconformable stratigraphic layers with lidar-derived surface roughness measurements to identify extensive wind-scour zones over Dome A, in the interior of East Antarctica. The scour zones are persistent because they are controlled by bedrock topography. On the basis of our Dome A observations, we develop an empirical model to predict wind-scour zones across the Antarctic continent and find that these zones are predominantly located in East Antarctica. We estimate that approx. 2.7-6.6% of the surface area of Antarctica has persistent negative net accumulation due to wind scour, which suggests that, across the continent, the snow mass input is overestimated by 11-36.5 Gt /yr in present surface-mass-balance calculations.

Antarctica↗

Utilization of satellite-derived cloud cover to improve the estimation of latent heat release in a winter extratropical cyclone

A simple scheme for utilizing satellite cloud pattern information to improve the estimation of parameterized stable and convective latent heat release is proposed and tested for an extratropical cyclone occurring over the United States during January 9-11, 1975. Calculations are done over the standard North American rawinsonde data region and over a sparse data region simulated by ignoring one-third of the original rawinsonde stations. Results indicate 15 percent and 30 percent improvements in latent heating estimates derived from the total and sparse data sets, respectively. A larger improvement was found using the sparse data set because the parameterized latent heat release estimates obtained using these data were poorer than those obtained using the total data set.

Lin, S.-J.↗

Anthropogenic Nitrogen Inputs and Impacts on Oceanic N2O Fluxes in the Northern Indian Ocean: the Need for an Integrated Observation and Modelling Approach

Anthropogenically-derived nitrogen input to the northern Indian Ocean has increased significantly in recent decades, based on both observational and model derived estimates This external nutrient source is supplied by atmospheric deposition and riverine fluxes, and has the potential to affect the vulnerable biogeochemical systems of the Arabian Sea and Bay of Bengal, influencing productivity and oceanic production of the greenhouse-gas nitrous-oxide (N2O). We summarize current estimates of this external nitrogen source to the northern Indian Ocean from observations and models, highlight implications for regional marine N2O emissions using model-based analyses, and make recommendations for measurement and model needs to improve current estimates and future predictions of this impact. Current observationally-derived estimates of deposition and riverine nitrogen inputs are limited by sparse measurements and uncertainties on accurate characterization of nitrogen species composition. Ocean model assessments of the impact of external nitrogen sources on regional marine N2O production in the northern Indian Ocean estimate potentially significant changes but also have large associated uncertainties. We recommend an integrated program of basin-wide measurements combined with high-resolution modeling and more detailed characterization of nitrogen-cycle process to address these uncertainties and improve current estimates and predictions.

Parvadha Suntharalingam↗

Canopy reflectance modeling in a tropical wooded grassland

Geometric/optical canopy reflectance modeling and spatial/spectral pattern recognition are used to study the form and structure of savanna in West Africa. An invertible plant canopy reflectance model is tested for its ability to estimate the amount of woody vegetation cover in areas of sparsely wooded grassland from remotely sensed data. Dry woodlands and wooded grasslands, commonly referred to as savannas, are important ecologically and economically in Africa, and cover approximately forty percent of the continent by some estimates. The Sahelian and Sudanian savanna make up the important and sensitive transition zone between the tropical forests and the arid Saharan region. The depletion of woody cover, used for fodder and fuel in these regions, has become a very severe problem for the people living there. LANDSAT Thematic Mapper (TM) data is used to stratify woodland and wooded grassland into areas of relatively homogeneous canopy cover, and then by applying an invertible forest canopy reflectance model to estimate directly the height and spacing of the trees in the stands. Since height and spacing are proportional to biomass in some cases, a successful application of the segmentation/modeling techniques will allow direct estimation of woody biomass, as well as cover density, over significant areas of these valuable and sensitive ecosystems. Sahelian savanna sites in the Gourma area of Mali being used by the NASA/GIMMS project (Global Inventory Modeling and Monitoring System, at Goddard Space Flight Center), in conjunction with CIPEA/Mali (Centre International pour l'Elevage en Afrique) will be used for testing the canopy model. The model will also be tested in a Sudanian zone crop/woodland area in the Region of Segou, Mali.

Simonett, D.↗

Algorithms for Autonomous GS Orbit Determination and Formation Flying

This final report for our study of autonomous Global Positioning System (GPS) satellite orbit determination comprises two sections. The first is the Ph.D. dissertation written by Michael C. Moreau entitled, "GPS Receiver Architecture for Autonomous Navigation in High Earth Orbits." Dr. Moreau's work was conducted under both this project and a NASA GSRP. His dissertation describes the key design features of a receiver specifically designed for autonomous operation in high earth orbits (HEO). He focused on the implementation and testing of these features for the GSFC PiVoT receiver. The second part is a memo describing a robust method for autonomous initialization of the orbit estimate given very little a priori information and sparse measurements. This is a key piece missing in the design of receivers for HEO.

Moreau, Michael C.↗

The OMPS Limb Profiler Instrument: Two-Dimensional Retrieval Algorithm

The upcoming Ozone Mapper and Profiler Suite (OMPS), which will be launched on the NPOESS Preparatory Project (NPP) platform in early 2011, will continue monitoring the global distribution of the Earth's middle atmosphere ozone and aerosol. OMPS is composed of three instruments, namely the Total Column Mapper (heritage: TOMS, OMI), the Nadir Profiler (heritage: SBUV) and the Limb Profiler (heritage: SOLSE/LORE, OSIRIS, SCIAMACHY, SAGE III). The ultimate goal of the mission is to better understand and quantify the rate of stratospheric ozone recovery. The focus of the paper will be on the Limb Profiler (LP) instrument. The LP instrument will measure the Earth's limb radiance (which is due to the scattering of solar photons by air molecules, aerosol and Earth surface) in the ultra-violet (UV), visible and near infrared, from 285 to 1000 nm. The LP simultaneously images the whole vertical extent of the Earth's limb through three vertical slits, each covering a vertical tangent height range of 100 km and each horizontally spaced by 250 km in the cross-track direction. Measurements are made every 19 seconds along the orbit track, which corresponds to a distance of about 150km. Several data analysis tools are presently being constructed and tested to retrieve ozone and aerosol vertical distribution from limb radiance measurements. The primary NASA algorithm is based on earlier algorithms developed for the SOLSE/LORE and SAGE III limb scatter missions. All the existing retrieval algorithms rely on a spherical symmetry assumption for the atmosphere structure. While this assumption is reasonable in most of the stratosphere, it is no longer valid in regions of prime scientific interest, such as polar vortex and UTLS regions. The paper will describe a two-dimensional retrieval algorithm whereby the ozone distribution is simultaneously retrieved vertically and horizontally for a whole orbit. The retrieval code relies on (1) a forward 2D Radiative Transfer code (to model limb radiances within a non-uniform atmosphere and evaluate 2D analytical partial derivatives) and (2) an optimal estimator inversion routine. The algorithm uses the typically sparse nature of the kernel matrices as well as fast matrix inversion techniques to allow for fast inversion of limb data with efficient memory management (as was done for MIPAS data processing). While the method has so far only been developed in the context of Single Scatter, the paper will show how the CPU intensive Multiple Scatter modeling can be implemented using parallel CPU processing. Initial results will be presented in terms of retrieved ozone profiles and code performance.

Rault, Didier F.↗

Development of the Ames Global Hyperspectral Synthetic Dataset

This study develops the surface BRDF (bidirectional reflectance distribution function) product of the Ames Global Hyperspectral Synthetic Dataset (AGHSD), based on the corresponding MODIS products, to support the NASA Surface Biology and Geology mission development. A main challenge in deriving a hyperspectral dataset from the multi-band satellite products is how to identify a succinct yet robust algorithm that allow us to infer BRDF at unobserved wavelengths based on the few observed bands. Using the theories of radiative transfer in vegetation canopies, we arrive at a simple equation that accurately approximates hyperspectral surface BRDF as the weighted sum of components from the soil and the vegetation. Each of the components is modeled by the product of the spectrally-dependent optical properties of a surface element (the spectra of the soil surface reflectance, the leaf single albedo, or the canopy scattering coefficient) and a spectrally-independent bidirectional scattering function. The optical properties of the soil and the vegetation can be obtained from existing spectral libraries or model simulations. The bidirectional scattering functions are represented by the Ross-Thick-Li-Sparse BRDF model, where the linear coefficients are estimated with regression analysis from the multi-band MODIS data. We validate the algorithm with simulations by Monte Carlo Ray Tracing model experiments, and the results are highly consistent with the theoretic derivation. We apply the algorithm to generate the AGHSD BRDF product at 1km and 8-day resolutions for the year of 2019. The results are biogeochemically and physically coherent and consistent, and thus serve the goal to support the science and application development of the SBG community.

Hyperspectral↗

Collaborative Pose Estimation of An Unknown Target Using Multiple Spacecraft

A reliable method for pose estimation of an unknown and uncooperative space target using monocular vision remains an open problem. Vision-based pose determination can be challenging in case of unfavorable illumination, time-varying conditions due to rotational motion and relative orbit, and scale ambiguity resolution. To address these challenges, we propose a novel collaborative pose determination algorithm called Multi- Spacecraft Simultaneous Estimation of Pose and Shape algorithm or M-SEPS.Within M-SEPS, a team of chaser spacecraft, each equipped with a monocular camera, exchange information over a local network to jointly estimate the relative kinematic state of the target and its sparse shape landmarks. In this approach, each spacecraft processes its own images and observes particular target landmarks in parallel and in a distributed fashion. Then, the local network is exploited by the spacecraft to share their consensus proposals and aggregate them to achieve the joint estimate. We validate our algorithm using simulations of relative orbits and observations, captured by each chaser spacecraft. To the best of the authors’ knowledge, this is the first cooperative, vision-based algorithm for estimating the pose and shape of a space object for an arbitrary number of spacecraft.

Chung, Soon-Jo↗

Demonstration of precise estimation of polar motion parameters with the global positioning system: Initial results

Data from the Global Positioning System (GPS) were used to determine precise polar motion estimates. Conservatively calculated formal errors of the GPS least squares solution are approx. 10 cm. The GPS estimates agree with independently determined polar motion values from very long baseline interferometry (VLBI) at the 5 cm level. The data were obtained from a partial constellation of GPS satellites and from a sparse worldwide distribution of ground stations. The accuracy of the GPS estimates should continue to improve as more satellites and ground receivers become operational, and eventually a near real time GPS capability should be available. Because the GPS data are obtained and processed independently from the large radio antennas at the Deep Space Network (DSN), GPS estimation could provide very precise measurements of Earth orientation for calibration of deep space tracking data and could significantly relieve the ever growing burden on the DSN radio telescopes to provide Earth platform calibrations.

Lichten, S. M.↗

Estimating the Noise Floor of Sonic Boom Metrics Across the USA

NASA is building the X-59 Quiet Supersonic Technology aircraft to produce low noise sonic booms for a series of community noise surveys across the USA. Survey participants will rate their perception of the low booms from supersonic X-59 flyovers. Several noise metrics are proposed to quantify the noise dose: A-, B-, D-, E-weighted Sound Exposure Level, Stevens Perceived Level, and Indoor Sonic Boom Annoyance Predictor. Sparse measurements across the survey area will be used to estimate community noise exposure. The level of these low booms may be comparable to the ambient noise level in some locations, leading to uncertainty in noise exposure estimations. This uncertainty may necessitate increased reliance on sonic boom propagation predictions for exposure estimation. Low-boom signal to ambient noise ratio is one way to quantify uncertainty in measured sonic boom levels. An empirical relationship between A-weighted ambient level and sonic boom metric levels is used in conjunction with the National Park Service’s L50 SPL map to estimate sonic boom metric ambient levels across the USA. The estimate of ambient sonic boom metric levels will aid in X-59 test planning and execution.

sonic boom↗