Engineering topics

Davis, J. L.

Publications and source records attributed to Davis, J. L..

At least 19 records

NASA NTRS2013

Assessment of the Mars Science Laboratory Entry, Descent, and Landing Simulation

On August 5, 2012, the Mars Science Laboratory rover, Curiosity, successfully landed inside Gale Crater. This landing was only the seventh successful landing and fourth rover to be delivered to Mars. Weighing nearly one metric ton, Curiosity is the largest and most complex rover ever sent to investigate another planet. Safely landing such a large payload required an innovative Entry, Descent, and Landing system, which included the first guided entry at Mars, the largest supersonic parachute ever flown at Mars, and a novel and untested Sky Crane landing system. A complete, end-to-end, six degree-of-freedom, multi-body computer simulation of the Mars Science Laboratory Entry, Descent, and Landing sequence was developed at the NASA Langley Research Center. In-flight data gathered during the successful landing is compared to pre-flight statistical distributions, predicted by the simulation. These comparisons provide insight into both the accuracy of the simulation and the overall performance of the vehicle.

Way, David W.↗

NASA NTRS2013

MARS Science Laboratory Post-Landing Location Estimation Using Post2 Trajectory Simulation

The Mars Science Laboratory (MSL) Curiosity rover landed safely on Mars August 5th, 2012 at 10:32 PDT, Earth Received Time. Immediately following touchdown confirmation, best estimates of position were calculated to assist in determining official MSL locations during entry, descent and landing (EDL). Additionally, estimated balance mass impact locations were provided and used to assess how predicted locations compared to actual locations. For MSL, the Program to Optimize Simulated Trajectories II (POST2) was the primary trajectory simulation tool used to predict and assess EDL performance from cruise stage separation through rover touchdown and descent stage impact. This POST2 simulation was used during MSL operations for EDL trajectory analyses in support of maneuver decisions and imaging MSL during EDL. This paper presents the simulation methodology used and results of pre/post-landing MSL location estimates and associated imagery from Mars Reconnaissance Orbiter s (MRO) High Resolution Imaging Science Experiment (HiRISE) camera. To generate these estimates, the MSL POST2 simulation nominal and Monte Carlo data, flight telemetry from onboard navigation, relay orbiter positions from MRO and Mars Odyssey and HiRISE generated digital elevation models (DEM) were utilized. A comparison of predicted rover and balance mass location estimations against actual locations are also presented.

Davis, J. L.↗

NASA NTRS2007

Global Distortion of GPS Networks Associated with Satellite Antenna Model Errors

Recent studies of the GPS satellite phase center offsets (PCOs) suggest that these have been in error by approx.1 m. Previous studies had shown that PCO errors are absorbed mainly by parameters representing satellite clock and the radial components of site position. On the basis of the assumption that the radial errors are equal, PCO errors will therefore introduce an error in network scale. However, PCO errors also introduce distortions, or apparent deformations, within the network, primarily in the radial (vertical) component of site position that cannot be corrected via a Helmert transformation. Using numerical simulations to quantify the effects of PC0 errors, we found that these PCO errors lead to a vertical network distortion of 6-12 mm per meter of PCO error. The network distortion depends on the minimum elevation angle used in the analysis of the GPS phase observables, becoming larger as the minimum elevation angle increases. The steady evolution of the GPS constellation as new satellites are launched, age, and are decommissioned, leads to the effects of PCO errors varying with time that introduce an apparent global-scale rate change. We demonstrate here that current estimates for PCO errors result in a geographically variable error in the vertical rate at the 1-2 mm/yr level, which will impact high-precision crustal deformation studies.

Cardellach, E.↗

NASA NTRS2007

GRACE Gravity Data Constrain Ancient Ice Geometries and Continental Dynamics over Laurentia

The free-air gravity trend over Canada, derived from the Gravity Recovery and Climate Experiment (GRACE) satellite mission, robustly isolates the gravity signal associated with glacial isostatic adjustment (GIA) from the longer-time scale mantle convection process. This trend proves that the ancient Laurentian ice complex was composed of two large domes to the west and east of Hudson Bay, in accord with one of two classes of earlier reconstructions. Moreover, GIA models that reconcile the peak rates contribute approximately 25 to approximately 45% to the observed static gravity field, which represents an important boundary condition on the buoyancy of the continental tectosphere.

Tamisiea, M. E.↗

NASA NTRS2006

Accuracy of High-Rate GPS for Seismology

We built a device for translating a GPS antenna on a positioning table to simulate the ground motions caused by an earthquake. The earthquake simulator is accurate to better than 0.1 mm in position, and provides the "ground truth" displacements for assessing the technique of high-rate GPS. We found that the root-mean-square error of the 1-Hz GPS position estimates over the 15-min duration of the simulated seismic event was 2.5 mm, with approximately 96% of the observations in error by less than 5 mm, and is independent of GPS antenna motion. The error spectrum of the GPS estimates is approximately flicker noise, with a 50% decorrelation time for the position error of approx.1.6 s. We that, for the particular event simulated, the spectrum of dependent error in the GPS measurements. surface deformations exceeds the GPS error spectrum within a finite band. More studies are required to determine whether a generally optimal bandwidth exists for a target group of seismic events.

Elosegui, P.↗

NASA NTRS2005

Constraining Hydrological and Cryospheric Mass Flux in Southeastern Alaska using Space-Based Gravity Measurements

Watersheds draining into the Gulf of Alaska (GoA) experience large seasonal and inter-annual variations of water in the form of rain, snow, and ice, but accurate constraints on these variations have been difficult to obtain. Over larger geographic regions, water variations can be inferred directly from the Gravity Recovery and Climate Experiment (GRACE) data. However, because GoA variations occur over such a small region, the inferred average value of water flux increases as the applied smoothing of the GRACE data decreases. We use this observed scaling together with scaling results obtained from forward models to infer a seasonal amplitude of 115 plus or minus 20 cubic kilometers of water and an average contribution to sea level rise over the two years of data of 0.31 plus or minus 0.09 millimeters per year. These results suggest that accelerated melting that began in the late 1990s, as inferred from altimetry, continues unabated.

Tamisiea, M. E.↗

NASA NTRS2005

Regional Deformation Studies with GRACE and GPS

GRACE data indicate large seasonal variations in gravity that have been shown to be to be related to climate-driven fluxes of surface water. Seasonal redistribution of surface mass deforms the Earth, and our previous study using GRACE data demonstrate that annual radial deformations of +/-13 mm in the region of Amazon River Basin were observed by both GRACE and ten GPS sites in the region. For the GRACE determinations, we estimate in a least-squares solution for each Stokes coefficient parameters that represent the amplitudes of the annual variation. We then filter these parameters based on a statistical test that uses the scatter of the postfit residuals. We demonstrate by comparison to the GPS amplitudes that this method is more accurate, for this region, than Gaussian smoothing. Our model for the temporal behavior of the gravity coefficients includes a rate term, and although the time series are noisy, the glacial isostatic adjustment signal over Hudson s Bay can be observed. .

Davis, J. L.↗

NASA NTRS2004

Climate-Driven Deformation of the Solid Earth from GRACE and GPS

GRACE data indicate large seasonal variations in gravity that are assumed to be related to climate-driven fluxes of space water. Seasonal redistribution of surface mass should deform the Earth, and our calculations using GRACE data suggest vertical deformations of approx. 13 mm in the region of greatest flux, the Amazon River Basin. To test the GRACE gravity-hydrology connection, we analyzed GPS data acquired from sites in this region. After accounting for degree 1 variations not observable with GRACE, we find that annual deformation measured with GPS correlates highly with predictions calculated from GRACE measurements. These results confirm the variations in surface water sensed by GRACE, which are significantly larger than those predicted by some hydrology models. The results also demonstrate that GRACE can be an important tool for monitoring deformation of the Earth, and suggest that combined analysis of GRACE and GPS may be a useful approach for estimation of geocenter variations.

Davis, J. L.↗

NASA NTRS1994

Determination of tidal h Love number parameters in the diurnal band using an extensive VLBI data set

We use over a decade of geodetic Very Long Baseline Interferometry (VLBI) data to estimate parameters in a resonance expansion of the frequency dependence of the tidal h(sub 2) Love number within the diurnal band. The resonance is associated with the retrograde free core nutation (RFCN). We obtain a value for the real part of the resonance strength of (-0.27 +/- 0.03) x 10(exp -3); a value of -0.19 x 10(exp -3) is predicted theoretically. Uncertainties in the VLBI estimates of the body tide radial displacement amplitudes are approximately 0.5 mm (1.1 mm for the K1 frequency), but they do not yield sufficiently small Love number uncertainties for placing useful constraints on the frequency of the RFCN, given the much smaller uncertainties obtained from independent analyses using nutation or gravimetric data. We also consider the imaginary part of the tidal h(sub 2) Love number. The estimated imaginary part of the resonance strength is (0.00 +/- 0.02) x 10(exp -3). The estimated imaginary part of the nonresonant component of the Love number implies a phase angle in the diurnal tidal response of the Earth of 0.7 deg +/- 0.5 deg (lag).

Mitrovica, J. X.↗

NASA NTRS1994

A spectral formalism for computing three-dimensional deformations due to surface loads. 1: Theory

We outline a complete spectral formalism for computing high spatial resolution three-dimensional deformations arising from the surface mass loading of a spherically symmetric planet. The main advantages of the formalism are that all surface mass loads are always described using a consistent mathematical representation and that calculations of deformation fields for various spatial resolutions can be performed by simpley altering the spherical harmonic degree truncation level of the procedure. The latter may be important when incorporating improved observational constraints on a particular surface mass load, when considering potential errors in the computed field associated with mass loading having a spatial scale unresolved by the observational constraints, or when treating a number of global surface mass loads constrained with different spatial resolutions. The advantages do not extend to traditional 'Green's function' approaches which involve surface element discretizations of the global mass loads. Another advantage of the spectral formalism, over the Green's function approach, is that a posteriori analyses of the computed deformation fields are easily performed. In developing the spectral formalism, we consider specific cases where the Earth's mantle is assumed to respond as an elastic, slightly anelastic, or linear viscoelastic medium. In the case of an elastic or slightly anelastic mantle rheology the spectral response equations incorporate frequency dependent Love numbers. The formalism can therefore be used, for example, to compute the potentially resonant deformational response associated with the free core nutation and Chandler wobble eigenfunctions. For completeness, the spectral response equations include both body forces, as arise from the gravitational attraction of the Sun and the Moon, and surface mass loads. In either case, and for both elastic and anelastic mantle rheologies, we outline a pseudo-spectral technique for computing the ocean adjustment associated with the total gravitational perturbation induced by the external forcing. Three-dimensional deformations computed using the usual Love number approach are generally referenced to an origin at the center of mass of the undeformed planet. We derive a spectral technique for transforming the results to an origin located at the center of mass of the deformed planet.

Mitrovica, J. X.↗

NASA NTRS1994

A spectral formalism for computing three-dimensional deformations due to surface loads. 2: Present-day glacial isostatic adjustment

Using a spherically symmetric, self-gravitating, linear viscoelastic Earth model, we predict present-day three-dimensional surface deformation rates and baseline evolutions arising as a consequence of the late Pleistocene glacial cycles. In general, we use realistic models for the space-time geometry of the final late Pleistocene deglaciation event and incorporate a gravitationally self-consistent ocean meltwater redistribution. The predictions of horizontal velocity presented differ significantly, in both their amplitude and their spatial variation, from those presented in earlier analysis of others which adopted simplified models of both the late Pleistocene ice history and the Earth rheology. An important characteristic of our predicted velocity fields is that the melting of the Laurentide ice sheet over Canada is capable of contributing appreciably to the adjustment in Europe. The sensitivity of the predictions to variations in mantle rheology is investigated by considering a number of different Earth models, and by computing appropriate Frechet kernels. These calculations suggest that the sensitivity of the deformations to the Earth's rheology is significant and strongly dependent on the location of the site relative to the ancient ice sheet. The effects on the predictions of three-dimensional deformation rates of altering the ice history or adopting approximate models for the ocean meltwater redistribution have also been considered and found to be important (the former especially so). Finally, for a suite of Earth models we provide predictions of the velocity of a number of baselines in North America and Europe. We find that, in general, both radial and tangential motions contribute significantly to baseline length changes, and that these contributions are a strong function of the Earth model. We have, furthermore, found a set of Earth models which, together with the ICE-3G deglaciation chronology, produce predictions of baseline length changes that are consistent with very long baseline interferometry measurements of baselines within Europe.

Mitrovica, J. X.↗

NASA NTRS1993

Constraining proposed combinations of ice history and Earth rheology using VLBI determined baseline length rates in North America

We predict the present-day rates of change of the lengths of 19 North American baselines due to the glacial isostatic adjustment process. Contrary to previously published research, we find that the three dimensional motion of each of the sites defining a baseline, rather than only the radial motions of these sites, needs to be considered to obtain an accurate estimate of the rate of change of the baseline length. Predictions are generated using a suite of Earth models and late Pleistocene ice histories, these include specific combinations of the two which have been proposed in the literature as satisfying a variety of rebound related geophysical observations from the North American region. A number of these published models are shown to predict rates which differ significantly from the VLBI observations.

Mitrovica, J. X.↗

NASA NTRS1993

Constraining proposed combinations of ice history and earth rheology using VLBI determined baseline length rates in North America

We predict the present-day rates of change of the lengths of 19 North American baselines due to the glacial isostatic adjustment process. Contrary to previously published research, we find that the three-dimensional motion of each of the sites defining a baseline, rather than only the radial motions of these sites, needs to be considered to obtain an accurate estimate of the rate of change of the baseline length. Predictions are generated using a suite of Earth models and late Pleistocene ice histories; these include specific combinations of the two which have been proposed in the literature as satisfying a variety of rebound related geophysical observations from the North American region. A number of these published models are shown to predict rates which differ significantly from the Very Long Base Interferometry (VLBI) observations.

Mitrovica, J. X.↗

NASA NTRS1991

Geodesy by radio interferometry - Water vapor radiometry for estimation of the wet delay

An important source of error in VLBI estimates of baseline length is unmodeled variations of the refractivity of the neutral atmosphere along the propagation path of the radio signals. This paper presents and discusses the method of using data from a water vapor radiomete (WVR) to correct for the propagation delay caused by atmospheric water vapor, the major cause of these variations. Data from different WVRs are compared with estimated propagation delays obtained by Kalman filtering of the VLBI data themselves. The consequences of using either WVR data or Kalman filtering to correct for atmospheric propagation delay at the Onsala VLBI site are investigated by studying the repeatability of estimated baseline lengths from Onsala to several other sites. The repeatability obtained for baseline length estimates shows that the methods of water vapor radiometry and Kalman filtering offer comparable accuracies when applied to VLBI observations obtained in the climate of the Swedish west coast. For the most frequently measured baseline in this study, the use of WVR data yielded a 13 percent smaller weighted-root-mean-square (WRMS) scatter of the baseline length estimates compared to the use of a Kalman filter. It is also clear that the 'best' minimum elevationi angle for VLBI observations depends on the accuracy of the determinations of the total propagation delay to be used, since the error in this delay increases with increasing air mass.

Elgered, G.↗

NASA NTRS1991

Effects of atmospheric modeling errors on determinations of baseline vectors from very long baseline interferometry

An investigation is conducted in order to expose errors in the formula for the elevation angle dependence of the hydrostatic atmospheric propagation delay, called the mapping function. A series of special VLBI experiments that involve a large fraction of group delay data from very low elevation angles are performed. Twenty-two experiments, each 26 hours in duration, were performed using VLBI antennas at Goldstone, California and at Westford, Massachusetts. Elevation angle cutoff tests were used to examine the contribution of these errors to the estimates of the vertical coordinate of site position. This contribution is determined to be about 19 mm. It is also found that, for the spacing of 1-2 months between experiments, the mapping function errors do not exhibit a coherent annual signature, but appear to be random in the long term.

Davis, J. L.↗

NASA NTRS1991

Peat analyses in the Hudson Bay Lowlands using ground penetrating radar

The use of ground penetrating radar (GPR) as a means to determine peak thickness and estimate peat volume in the Hudson Bay Lowlands of Canada is examined. Ground-based and airborne GPR data were acquired so as to extrapolate measurements to larger scales. While the ground-based measurements did an excellent job in determining peat depth, the airborne techniques did a fair job a low altitudes and demonstrated great promise with additional system engineering changes.

Pelletier, R. E.↗

NASA NTRS1988

VLBI geodesy - 2 parts-per-billion precision in length determinations for transcontinental baselines

VLBI was to make twenty-two independent measurements, between September 1984 and December 1986, of the length of the 3900-km baseline between the Mojave site in California and the Haystack/Westford site in Massachusetts. These experiments differ from the typical geodetic VLBI experiments in that a large fraction of observations is obtained at elevation angles between 4 and 10 deg. Data from these low elevation angles allow the vertical coordinate of site position, and hence the baseline length, to be estimated with greater precision. For the sixteen experiments processed thus far, the weighted root-mean-square scatter of the estimates of the baseline length is 8 mm.

Davis, J. L.↗

NASA NTRS1988

Methods of correction for the 'wet' atmosphere in estimating baseline lengths from VLBI

The error in VLBI estimates of baseline length caused by unmodelled variations in the propagation path through the atmosphere is greater for longer baselines. Series of estimates of baseline lengths obtained using different methods to correct for the propagation delay caused by atmospheric water vapor are presented and discussed. The main methods are use of data from: a water-vapor radiometer (WVR) and Kalman-filtering of the VLBI data themselves to estimate the propagation delay. Since the longest timespan of WVR data associated with geodetic VLBI experiments was obtained at the Onsala Space Observatory in Sweden, results are presented for the following three baselines: (1) Onsala-Wettzell, FRG (920 km), (2) Onsala-Haystack/Westford, MA (5600 km), and (3) Onsala-Owens Valley (7914 km).

Elgered, G.↗