Search NASASearch

SEARCH · Search NASA

Results for “Gradient 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 37 records · Page 2

Estimation of terrain iso-gradients from a stochastic range data measurement matrix

The problem of estimating terrain iso-gradients for an autonomous roving vehicle is complicated by the measurement error inherent in the range data matrix. In this paper, the in-path and cross-path slopes at each data point in the multivariable range matrix are expressed in terms of parameters at that data point and those at its surrounding points. The sensitivity in the change of these slopes - due to the errors in measurements in range, elevation angle and azimuth angle are formulated. These sensitivities are used to determine the appropriate statistics of the gradient, which are then used to describe the terrain with a known probability of accuracy.

Shen, C. N.

Estimates of Terms in Ohm's Law During an Encounter with an Electron Diffusion Region

We present measurements from the Magnetospheric Multiscale (MMS) mission taken during a reconnection event on the dayside magnetopause which includes a passage through an electron diffusion region (EDR). The four MMS satellites were separated by about 10 km such that estimates of gradients and divergences allow a reasonable estimate of terms in the generalized Ohm's law, which is key to investigating the energy dissipation during reconnection. The strength and character of dissipation mechanisms determines how magnetic energy is released. We show that both electron pressure gradients and electron inertial effects are important, but not the only participants in reconnection near EDRs, since there are residuals of a few mVm (approximately 30-50%) of E+ U(sub e) x B (from the sum of these two terms) during the encounters. These results are compared to a simulation, which exhibits many of the observed features, but where relatively little residual is present.

Torbert, R. B.

Cosmic ray gradients in the heliosphere

Early estimates of the radial gradient made in the inner heliosphere are presented. In a discussion of discoveries in the outer heliosphere, consideration is given to low radial gradients and nonradial transport, the identification of different components, the large-scale organization of the heliospheric magnetic field, propagating modulation features, and latitude gradients. Current estimates of the spatial distribution of the cosmic ray intensity are presented as well.

Fillius, Walker

Mitral inertance in humans: critical factor in Doppler estimation of transvalvular pressure gradients

The pressure-velocity relationship across the normal mitral valve is approximated by the Bernoulli equation DeltaP = 1/2 rhoDeltav(2) + M. dv/dt, where DeltaP is the atrioventricular pressure difference, rho is blood density, v is transmitral flow velocity, and M is mitral inertance. Although M is indispensable in assessing transvalvular pressure differences from transmitral flow, this term is poorly understood. We measured intraoperative high-fidelity left atrial and ventricular pressures and simultaneous transmitral flow velocities by using transesophageal echocardiography in 100 beats (8 patients). We computed mean mitral inertance (M) by M = integral((DeltaP)-(1/2 x rho v(2))dt/integral(dv/dt)dt and we assessed the effect of the inertial term on the transmitral pressure-flow relation. ranged from 1.03 to 5.96 g/cm(2) (mean = 3.82 +/- 1.22 g/cm(2)). DeltaP calculated from the simplified Bernoulli equation (DeltaP = 1/2. rhov(2)) lagged behind (44 +/- 11 ms) and underestimated the actual peak pressures (2.3 +/- 1.1 mmHg). correlated with left ventricular systolic pressure (r = -0.68, P < 0.0001) and transmitral pressure gradients (r = 0.65, P < 0.0001). Because mitral inertance causes the velocity to lag significantly behind the actual pressure gradient, it needs to be considered when assessing diastolic filling and the pressure difference across normal mitral valves.

NASA Discipline Cardiopulmonary

Revisiting the Least-squares Procedure for Gradient Reconstruction on Unstructured Meshes

The accuracy of the least-squares technique for gradient reconstruction on unstructured meshes is examined. While least-squares techniques produce accurate results on arbitrary isotropic unstructured meshes, serious difficulties exist for highly stretched meshes in the presence of surface curvature. In these situations, gradients are typically under-estimated by up to an order of magnitude. For vertex-based discretizations on triangular and quadrilateral meshes, and cell-centered discretizations on quadrilateral meshes, accuracy can be recovered using an inverse distance weighting in the least-squares construction. For cell-centered discretizations on triangles, both the unweighted and weighted least-squares constructions fail to provide suitable gradient estimates for highly stretched curved meshes. Good overall flow solution accuracy can be retained in spite of poor gradient estimates, due to the presence of flow alignment in exactly the same regions where the poor gradient accuracy is observed. However, the use of entropy fixes has the potential for generating large but subtle discretization errors.

Mavriplis, Dimitri J.

Numerical optimization in Hilbert space using inexact function and gradient evaluations

Trust region algorithms provide a robust iterative technique for solving non-convex unstrained optimization problems, but in many instances it is prohibitively expensive to compute high accuracy function and gradient values for the method. Of particular interest are inverse and parameter estimation problems, since function and gradient evaluations involve numerically solving large systems of differential equations. A global convergence theory is presented for trust region algorithms in which neither function nor gradient values are known exactly. The theory is formulated in a Hilbert space setting so that it can be applied to variational problems as well as the finite dimensional problems normally seen in trust region literature. The conditions concerning allowable error are remarkably relaxed: relative errors in the gradient error condition is automatically satisfied if the error is orthogonal to the gradient approximation. A technique for estimating gradient error and improving the approximation is also presented.

Carter, Richard G.

Large Volcanic Rises on Venus

Large volcanic rises on Venus have been interpreted as hotspots, or the surface manifestation of mantle upwelling, on the basis of their broad topographic rises, abundant volcanism, and large positive gravity anomalies. Hotspots offer an important opportunity to study the behavior of the lithosphere in response to mantle forces. In addition to the four previously known hotspots, Atla, Bell, Beta, and western Eistla Regiones, five new probable hotspots, Dione, central Eistla, eastern Eistla, Imdr, and Themis, have been identified in the Magellan radar, gravity and topography data. These nine regions exhibit a wider range of volcano-tectonic characteristics than previously recognized for venusian hotspots, and have been classified as rift-dominated (Atla, Beta), coronae-dominated (central and eastern Eistla, Themis), or volcano-dominated (Bell, Dione, western Eistla, Imdr). The apparent depths of compensation for these regions ranges from 65 to 260 km. New estimates of the elastic thickness, using the 90 deg and order spherical harmonic field, are 15-40 km at Bell Regio, and 25 km at western Eistla Regio. Phillips et al. find a value of 30 km at Atla Regio. Numerous models of lithospheric and mantle behavior have been proposed to interpret the gravity and topography signature of the hotspots, with most studies focusing on Atla or Beta Regiones. Convective models with Earth-like parameters result in estimates of the thickness of the thermal lithosphere of approximately 100 km. Models of stagnant lid convection or thermal thinning infer the thickness of the thermal lithosphere to be 300 km or more. Without additional constraints, any of the model fits are equally valid. The thinner thermal lithosphere estimates are most consistent with the volcanic and tectonic characteristics of the hotspots. Estimates of the thermal gradient based on estimates of the elastic thickness also support a relatively thin lithosphere (Phillips et al.). The advantage of larger estimates of the thermal lithospheric thickness is that they provide an explanation for the apparently modest levels of geologic activity on Venus over the last half billion years.

Smrekar, Suzanne E.

Plasmapause diffusion

The Bohm diffusion coefficient and observed electrostatic wave scattering are used as the bases of estimates of the smoothing effect that diffusion may have on steep plasmapause density gradients. The estimate for diffusion resulting from scattering by observed electrostatic waves is found to be much lower than that of the perpendicular Bohm diffusion coefficient for characteristic plasma temperatures and magnetic fields. This diffusion rate estimate may be too small, however, if the wave amplitudes are significantly higher for steep plasmapauses. The effects are therefore negligible for most considerations of macroscopic plasmapause dynamics, but may be significant in limiting drift wave instabilities and similar phenomena driven by the steepness of the plasmapause density gradient.

Horwitz, J. L.

Multi-Pass Sequential Mini-Batch Stochastic Gradient Descent Algorithms for Noise Covariance Estimation in Adaptive Kalman Filtering

Estimation of unknown noise covariances in a Kalman filter is a problem of significant practical interest in a wide array of applications. Although this problem has a long history, reliable algorithms for their estimation were scant, and necessary and sufficient conditions for identifiability of the covariances were in dispute until recently. Necessary and sufficient conditions for covariance estimation and a batch estimation algorithm were presented in our previous study. This paper presents stochastic gradient descent algorithms for noise covariance estimation in adaptive Kalman filters that are an order of magnitude faster than the batch method for similar or better root mean square error. More significantly, these algorithms are applicable to non-stationary systems where the noise covariances can occasionally jump up or down by an unknown magnitude. The computational efficiency of the new algorithms stems from adaptive thresholds for convergence, recursive fading memory estimation of the sample cross-correlations of the innovations, and accelerated stochastic gradient descent algorithms. The comparative evaluation of the proposed methods on a number of test cases demonstrates their computational efficiency and accuracy.

Adaptive Kalman filtering

Effect of land use change on the carbon cycle in Amazon soils

The overall goal of this study was to provide a quantitative understanding of the cycling of carbon in the soils associated with deep-rooting Amazon forests. In particular, we wished to apply the understanding gained by answering two questions: (1) what changes will accompany the major land use change in this region, the conversion of forest to pasture? and (2) what is the role of carbon stored deeper than one meter in depth in these soils? To construct carbon budgets for pasture and forest soils we combined the following: measurements of carbon stocks in above-ground vegetation, root biomass, detritus, and soil organic matter; rates of carbon inputs to soil and detrital layers using litterfall collection and sequential coring to estimate fine root turnover; C-14 analyses of fractionated SOM and soil CO2 to estimate residence times; C-13 analyses to estimate C inputs to pasture soils from C-4 grasses; soil pCO2, volumetric water content, and radon gradients to estimate CO2 production as a function of soil depth; soil respiration to estimate total C outputs; and a model of soil C dynamics that defines SOM fractions cycling on annual, decadal, and millennial time scales.

Trumbore, Susan E.

Prospects for UT1 Measurements from VLBI Intensive Sessions

Very Long Baseline Interferometry (VLBI) Intensives are one-hour single baseline sessions to provide Universal Time (UT1) in near real-time up to a delay of three days if a site is not e-transferring the observational data. Due to the importance of UT1 estimates for the prediction of Earth orientation parameters, as well as any kind of navigation on Earth or in space, there is not only the need to improve the timeliness of the results but also their accuracy. We identify the asymmetry of the tropospheric delays as the major error source, and we provide two strategies to improve the results, in particular of those Intensives which include the station Tsukuba in Japan with its large tropospheric variation. We find an improvement when (1) using ray-traced delays from a numerical weather model, and (2) when estimating tropospheric gradients within the analysis of Intensive sessions. The improvement is shown in terms of reduction of rms of length-of-day estimates w.r.t. those derived from Global Positioning System observations

Boehm, Johannes

Retrieval Assimilation and Modeling of Atmospheric Water Vapor from Ground- and Space-Based GPS Networks: Investigation of the Global and Regional Hydrological Cycles

Uncertainty over the response of the atmospheric hydrological cycle (particularly the distribution of water vapor and cloudiness) to anthropogenic forcing is a primary source of doubt in current estimates of global climate sensitivity, which raises severe difficulties in evaluating its likely societal impact. Fortunately, a variety of advanced techniques and sensors are beginning to shed new light on the atmospheric hydrological cycle. One of the most promising makes use of the sensitivity of the Global Positioning System (GPS) to the thermodynamic state, and in particular the water vapor content, of the atmosphere through which the radio signals propagate. Our strategy to derive the maximum benefit for hydrological studies from the rapidly increasing GPS data stream will proceed in three stages: (1) systematically analyze and archive quality-controlled retrievals using state-of-the-art techniques; (2) employ both currently available and innovative assimilation procedures to incorporate these determinations into advanced regional and global atmospheric models and assess their effects; and (3) apply the results to investigate selected scientific issues of relevance to regional and global hydrological studies. An archive of GPS-based estimation of total zenith delay (TZD) data and water vapor where applicable has been established with expanded automated quality control. The accuracy of the GPS estimates is being monitored; the investigation of systematic errors is ongoing using comparisons with water vapor radiometers. Meteorological packages have been implemented. The accuracy and utilization of the TZD estimates has been improved by implementing a troposphere gradient model. GPS-based gradients have been validated as real atmospheric moisture gradients, establishing a link between the estimated gradients and the passage of weather fronts. We have developed a generalized ray tracing inversion scheme that can be used to analyze occultation data acquired from space- or land-based receivers. The National Center for Atmospheric Research mesoscale model (version MM5) has been adapted for Southern California, and assimilation studies are underway. Additional information is contained in the original.

Dickey, Jean O.

Biomass Estimation from Simulated GEDI, ICESat-2 and NISAR Across Environmental Gradients in Sonoma County, California

Estimates of the magnitude and distribution of aboveground carbon in Earth’s forests remain uncertain, yet knowledge of forest carbon content at a global scale is critical for forest management in support of climate mitigation. In light of this knowledge gap, several upcoming spaceborne missions aim to map forest aboveground biomass, and many new biomass products are expected from these datasets. As these new missions host different technologies, each with relative strengths and weaknesses for biomass retrieval, as well as different spatial resolutions, consistently comparing or combining biomass estimates from these new datasets will be challenging. This paper presents a demonstration of an inter-comparison of biomass estimates from simulations of three NASA missions (GEDI, ICESat-2 and NISAR) over Sonoma county in California, USA. We use a high resolution, locally calibrated airborne lidar map as our reference dataset, and emphasize the importance of considering uncertainties in both reference maps and spaceborne estimates when conducting biomass product validation. GEDI and ICESat-2 were simulated from airborne lidar point clouds, while UAVSAR’s L-band backscatter was used as a proxy for NISAR. To estimate biomass for the lidar missions we used GEDI’s footprint-level biomass algorithms, and also adapted these for application to ICESat-2. For UAVSAR, we developed a locally trained biomass model, calibrated against the ALS reference map. Each mission simulation was evaluated in comparison to the local reference map at its native product resolution (25 m, 100m transect, and 1 ha) yielding RMSEs of 57%, 75%, and 89% for GEDI, NISAR, and ICESat-2 respectively. RMSE values increased for GEDI’s power beam during simulated daytime conditions (64%), coverage beam during nighttime conditions (72%), and coverage beam daytime conditions (87%). We also test the application of GEDI’s biomass modeling framework for estimation of biomass from ICESat-2, and fine that ICESat-2 yields reasonable biomass estimates, particularly in relatively short, open canopies. Results suggest that while all three missions will produce datasets useful for biomass mapping, tall, dense canopies such as those found in Sonoma County present the greatest challenges for all three missions, while steep slopes also prove challenging for single-date SAR based biomass retrieval. Our methods provide guidance for the inter-comparison and validation of spaceborne biomass estimates through the use of airborne lidar reference maps, and could be repeated with on-orbit estimates in any area with high quality field plot and ALS data. These methods allow for regional interpretations and filtering of multi-mission biomass estimates toward improved wall-to-wall biomass maps through data fusion.

GEDI

Estimation of diastolic intraventricular pressure gradients by Doppler M-mode echocardiography

Previous studies have shown that small intraventricular pressure gradients (IVPG) are important for efficient filling of the left ventricle (LV) and as a sensitive marker for ischemia. Unfortunately, there has previously been no way of measuring these noninvasively, severely limiting their research and clinical utility. Color Doppler M-mode (CMM) echocardiography provides a spatiotemporal velocity distribution along the inflow tract throughout diastole, which we hypothesized would allow direct estimation of IVPG by using the Euler equation. Digital CMM images, obtained simultaneously with intracardiac pressure waveforms in six dogs, were processed by numerical differentiation for the Euler equation, then integrated to estimate IVPG and the total (left atrial to left ventricular apex) pressure drop. CMM-derived estimates agreed well with invasive measurements (IVPG: y = 0.87x + 0.22, r = 0.96, P < 0.001, standard error of the estimate = 0.35 mmHg). Quantitative processing of CMM data allows accurate estimation of IVPG and tracking of changes induced by beta-adrenergic stimulation. This novel approach provides unique information on LV filling dynamics in an entirely noninvasive way that has previously not been available for assessment of diastolic filling and function.

Non-NASA Center

Atmospheric gradients from very long baseline interferometry observations

Azimuthal asymmetries in the atmospheric refractive index can lead to errors in estimated vertical and horizontal station coordinates. Daily average gradient effects can be as large as 50 mm of delay at a 7 deg elevation. To model gradients, the constrained estimation of gradient paramters was added to the standard VLBI solution procedure. Here the analysis of two sets of data is summarized: the set of all geodetic VLBI experiments from 1990-1993 and a series of 12 state-of-the-art R&D experiments run on consecutive days in January 1994. In both cases, when the gradient parameters are estimated, the overall fit of the geodetic solution is improved at greater than the 99% confidence level. Repeatabilities of baseline lengths ranging up to 11,000 km are improved by 1 to 8 mm in a root-sum-square sense. This varies from about 20% to 40% of the total baseline length scatter without gradient modeling for the 1990-1993 series and 40% to 50% for the January series. Gradients estimated independently for each day as a piecewise linear function are mostly continuous from day to day within their formal uncertainties.

Macmillan, D. S.

Parameter estimation for terrain modeling from gradient data

A method is developed for modeling terrain surfaces for use on an unmanned Martian roving vehicle. The modeling procedure employs a two-step process which uses gradient as well as height data in order to improve the accuracy of the model's gradient. Least square approximation is used in order to stochastically determine the parameters which describe the modeled surface. A complete error analysis of the modeling procedure is included which determines the effect of instrumental measurement errors on the model's accuracy. Computer simulation is used as a means of testing the entire modeling process which includes the acquisition of data points, the two-step modeling process and the error analysis. Finally, to illustrate the procedure, a numerical example is included.

Dangelo, K. R.