Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geodetics”

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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Manual for a workstation-based generic flight simulation program (LaRCsim), version 1.4

LaRCsim is a set of ANSI C routines that implement a full set of equations of motion for a rigid-body aircraft in atmospheric and low-earth orbital flight, suitable for pilot-in-the-loop simulations on a workstation-class computer. All six rigid-body degrees of freedom are modeled. The modules provided include calculations of the typical aircraft rigid-body simulation variables, earth geodesy, gravity and atmospheric models, and support several data recording options. Features/limitations of the current version include English units of measure, a 1962 atmosphere model in cubic spline function lookup form, ranging from sea level to 75,000 feet, rotating oblate spheroidal earth model, with aircraft C.G. coordinates in both geocentric and geodetic axes. Angular integrations are done using quaternion state variables Vehicle X-Z symmetry is assumed.

Jackson, E. Bruce↗
NASA NTRS1995

FORTRAN program for analyzing ground-based radar data: Usage and derivations, version 6.2

A postflight FORTRAN program called 'radar' reads and analyzes ground-based radar data. The output includes position, velocity, and acceleration parameters. Air data parameters are also provided if atmospheric characteristics are input. This program can read data from any radar in three formats. Geocentric Cartesian position can also be used as input, which may be from an inertial navigation or Global Positioning System. Options include spike removal, data filtering, and atmospheric refraction corrections. Atmospheric refraction can be corrected using the quick White Sands method or the gradient refraction method, which allows accurate analysis of very low elevation angle and long-range data. Refraction properties are extrapolated from surface conditions, or a measured profile may be input. Velocity is determined by differentiating position. Accelerations are determined by differentiating velocity. This paper describes the algorithms used, gives the operational details, and discusses the limitations and errors of the program. Appendices A through E contain the derivations for these algorithms. These derivations include an improvement in speed to the exact solution for geodetic altitude, an improved algorithm over earlier versions for determining scale height, a truncation algorithm for speeding up the gradient refraction method, and a refinement of the coefficients used in the White Sands method for Edwards AFB, California. Appendix G contains the nomenclature.

Haering, Edward A., Jr.↗
NASA NTRS1994

Detection of atmospheric pressure loading using very long baseline interferometry measurements

Loading of the Earth by the temporal redistribution of global atmospheric mass is likely to displace the positions of geodetic monuments by tens of millimeters both vertically and horizontally. Estimates of these displacements are determined by convolving National Meteorological Center (NMC) global values of atmospheric surface pressure with Farrell's elastic Green's functions. An analysis of the distances between radio telescopes determined by very long baseline interferometry (VLBI) between 1984 and 1992 reveals that in many of the cases studied there is a significant contribution to baseline length change due to atmospheric pressure loading. Our analysis covers intersite distances of between 1000 and 10,000 km and is restricted to those baselines measured more than 100 times. Accounting for the load effects (after first removing a best fit slope) reduces the weighted root-mean-square (WRMS) scatter of the baseline length residuals on 11 of the 22 baselines investigated. The slight degradation observed in the WRMS scatter on the remaining baselines is largely consistent with the expected statistical fluctuations when a small correction is applied to a data set having a much larger random noise. The results from all baselines are consistent with approximately 60% of the computed pressure contribution being present in the VLBI length determinations. Site dependent coefficients determined by fitting local pressure to the theoretical radial displacement are found to reproduce the deformation caused by the regional pressure to within 25% for most inland sites. The coefficients are less reliable at near coastal and island stations.

Vandam, T. M.↗
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

Strain accumulation in the Shumagin Islands: Results of initial GPS measurements

Deformation in the Shumagin seismic gap has been monitored with repeated trilateration (EDM) in the 1980-1987 interval and with the Global Positioning System (GPS) in the 1987-1991 interval. The geodetic network extends for 100-km across the Shumagin Islands to the Alaska Peninsula. Results from the GPS surveys are consistent with those previously reported for the EDM surveys: we failed to detect significant strain accumulation in the N30 deg W direction of plate convergence. Using the method of simultaneous reduction for position and strain rates, we found the average rate of extension in the direction of plate convergence to be -25 +/- 25 nanostrain/yr (nstrain/yr) during the 1987-1991 interval of GPS surveys compared with -20 +/- 15 nstrain/yr during the 1981-1987 interval of complete EDM surveys. We found a marginally significant -26 +/- 12 nstrain/yr extension rate in the 1981-1991 interval covered by the combined EDM and GPS surveys. Strain rates are higher, but not significantly so, in the part of the network closest to the trench. Spatial variation in the deformation is observed in the 1980-1991 average station velocities, where three of the four stations closest to the trench have an arcward velocity of a few mm/yr. The observed strain rates are an order of magnitude lower than the -200 nstrain/yr rate predicted by dislocation models.

Larson, Kristine M.↗
NASA NTRS1994

Site velocities before and after the Loma Prieta and Gulf of Alaska earthquakes determined from VLBI

We use geodetic data from Very Long Baseline Interferometry (VLBI) to determine the pre- and postseismic velocities of two sites. We then place limits on variations in interseismic strain buildup. The 1987 and 1988 Gulf of Alaska earthquakes (each Ms = 7.6) broke the Pacific plate interior. During the earthquakes the Cape Yakataga site moved 78 mm toward southwest. During the 1989 Loma Prieta earthquake (Ms = 7.1) the Fort Ord site moved 48 mm toward north. Baselines (a) from Fairbanks to Cape Yakataga and (b) from Mojave to Fort Ord change at nearly the same rate before and after the earthquakes. Postseismic transients, which we determine from differences between post- and preseismic rates, are minor: at Cape Yakataga the transient is 3 +/- 4 mm in a postseismic interval of 23 months, and at Fort Ord the transient is 6 +/- 5 mm in 21 months. The slip beneath the Loma Prieta rupture needed to generate the Fort Ord transient is 0.22 +/- 0.19 m, one-tenth the coseismic slip (2 m). We analyze elastic lithosphere-viscous asthenosphere models to determine that the characteristic time describing exponential decay in deep fault slip is longer than 6 years. The VLBI measurements are consistent with uniform interseismic strain buildup. They disagree with fast postseismic rates caused by an asthenosphere with very low viscosity.

Argus, Donald F.↗
NASA NTRS1994

Vertical ocean loading amplitudes from VLBI measurements

Vertical ocean loading amplitudes are determined by analysis of IRIS geodetic Very Long Baseline Interferometry (VLBI) data. The 4 nearly diurnal (K(sub 1), P(sub 1), O(sub 1), Q(sub 1)) and 4 nearly semidiurnal (K(sub 2), S(sub 2), M(sub 2), and N(sub 2)) component amplitudes can be inferred from the data with accuracies of 1-2 mm. Uncertainties of total displacements can approach 1 cm. Empirically determined total displacements are considerably larger than the values calculated from two geophysical models. The Scherneck model is found to give a better representation of VLBI delay data than the model of Pagiatakis by about 3 mm in residuals at 6 sites. Empirical estimation of the ocean loading amplitudes reduces Chi squared by 3067 for the 96 additional degrees of freedom in a fit to 273,000 IRIS VLBI observations, and reduces the RMS residuals by 3 mm relative to a fit using the fixed Scherneck model. Vertical ocean loading amplitudes can thus be inferred from VLBI data at a level which improves the overall model, but detailed assessment of individual tidal components is presently obscured by incomplete modeling at the tidal frequencies.

Sovers, O. J.↗
NASA NTRS1994

Evaluation of very long baseline interferometry atmospheric modeling improvements

We determine the improvement in baseline length precision and accuracy using new atmospheric delay mapping functions and MTT by analyzing the NASA Crustal Dynamics Project research and development (R&D) experiments and the International Radio Interferometric Surveying (IRIS) A experiments. These mapping functions reduce baseline length scatter by about 20% below that using the CfA2.2 dry and Chao wet mapping functions. With the newer mapping functions, average station vertical scatter inferred from observed length precision (given by length repeatabilites) is 11.4 mm for the 1987-1990 monthly R&D series of experiments and 5.6 mm for the 3-week-long extended research and development experiment (ERDE) series. The inferred monthly R&D station vertical scatter is reduced by 2 mm or by 7 mm is a root-sum-square (rss) sense. Length repeatabilities are optimum when observations below a 7-8 deg elevation cutoff are removed from the geodetic solution. Analyses of IRIS-A data from 1984 through 1991 and the monthly R&D experiments both yielded a nonatmospheric unmodeled station vertical error or about 8 mm. In addition, analysis of the IRIS-A exeriments revealed systematic effects in the evolution of some baseline length measurements. The length rate of change has an apparent acceleration, and the length evolution has a quasi-annual signature. We show that the origin of these effects is unlikely to be related to atmospheric modeling errors. Rates of change of the transatlantic Westford-Wettzell and Richmond-Wettzell baseline lengths calculated from 1988 through 1991 agree with the NUVEL-1 plate motion model (Argus and Gordon, 1991) to within 1 mm/yr. Short-term (less than 90 days) variations of IRIS-A baseline length measurements contribute more than 90% of the observed scatter about a best fit line, and this short-term scatter has large variations on an annual time scale.

Macmillan, D. S.↗
NASA NTRS1994

Gravity model improvement using the DORIS tracking system on the SPOT 2 satellite

A high-precision radiometric satellite tracking system, Doppler Orbitography and Radio-positioning Integrated by Satellite system (DORIS), has recently been developed by the French space agency, Centre National d'Etudes Spatiales (CNES). DORIS was designed to provide tracking support for missions such as the joint United States/French TOPEX/Poseidon. As part of the flight testing process, a DORIS package was flown on the French SPOT 2 satellite. A substantial quantity of geodetic quality tracking data was obtained on SPOT 2 from an extensive international DORIS tracking network. These data were analyzed to assess their accuracy and to evaluate the gravitational modeling enhancements provided by these data in combination with the Goddard Earth Model-T3 (GEM-T3) gravitational model. These observations have noise levels of 0.4 to 0.5 mm/s, with few residual systematic effects. Although the SPOT 2 satellite experiences high atmospheric drag forces, the precision and global coverage of the DORIS tracking data have enabled more extensive orbit parameterization to mitigate these effects. As a result, the SPOT 2 orbital errors have been reduced to an estimated radial accuracy in the 10-20 cm RMS range. The addition of these data, which encompass many regions heretofore lacking in precision satellite tracking, has significantly improved GEM-T3 and allowed greatly improved orbit accuracies for Sun-synchronous satellites like SPOT 2 (such as ERS 1 and EOS). Comparison of the ensuing gravity model with other contemporary fields (GRIM-4C2, TEG2B, and OSU91A) provides a means to assess the current state of knowledge of the Earth's gravity field. Thus, the DORIS experiment on SPOT 2 has provided a strong basis for evaluating this new orbit tracking technology and has demonstrated the important contribution of the DORIS network to the success of the TOPEX/Poseidon mission.

Nerem, R. S.↗
NASA NTRS1994

Modeling of periodic great earthquakes on the San Andreas fault: Effects of nonlinear crustal rheology

We analyze the cycle of great earthquakes along the San Andreas fault with a finite element numerical model of deformation in a crust with a nonlinear viscoelastic rheology. The viscous component of deformation has an effective viscosity that depends exponentially on the inverse absolute temperature and nonlinearity on the shear stress; the elastic deformation is linear. Crustal thickness and temperature are constrained by seismic and heat flow data for California. The models are for anti plane strain in a 25-km-thick crustal layer having a very long, vertical strike-slip fault; the crustal block extends 250 km to either side of the fault. During the earthquake cycle that lasts 160 years, a constant plate velocity v(sub p)/2 = 17.5 mm yr is applied to the base of the crust and to the vertical end of the crustal block 250 km away from the fault. The upper half of the fault is locked during the interseismic period, while its lower half slips at the constant plate velocity. The locked part of the fault is moved abruptly 2.8 m every 160 years to simulate great earthquakes. The results are sensitive to crustal rheology. Models with quartzite-like rheology display profound transient stages in the velocity, displacement, and stress fields. The predicted transient zone extends about 3-4 times the crustal thickness on each side of the fault, significantly wider than the zone of deformation in elastic models. Models with diabase-like rheology behave similarly to elastic models and exhibit no transient stages. The model predictions are compared with geodetic observations of fault-parallel velocities in northern and central California and local rates of shear strain along the San Andreas fault. The observations are best fit by models which are 10-100 times less viscous than a quartzite-like rheology. Since the lower crust in California is composed of intermediate to mafic rocks, the present result suggests that the in situ viscosity of the crustal rock is orders of magnitude less the rock viscosity determined in the laboratory.

Reches, Ze'ev↗
NASA NTRS1994

The analysis of star catalogs. 1: An intercomparison among the FK3, the FK4, and the FK5

Progress in the construction of the fundamental catalogs allows one to ask the following questions: 'Now that we have the FK5 with which to measure by, how good was the FK4? The FK3? How realistic were the error estimates they gave for their equatorial coordinates and their proper motions? Have their residual systematic differences continued to decrease in amplitude?' and so forth. Similar questions could be asked of the GC and the N30, and a companion paper will similarly analyze them. Finally, the identical series of questions could be asked of the AGK2, AGK3, AGK3U sequence. In this paper we develop a general method to objectively investigate, post facto, the true quality of star catalog error estimates. We present a complete mathematical formulation of the questions above and delineate the calculations needed to address them. Our initial numerical applications are restricted to a straightforward examination of the mean errors of the coordinates and proper motions in the FK3 (using the FK4 and the FK5 for comparison) and in the FK4 (using the FK5 for comparison). However, we can simply and consistently interpret the results of these computations -- for the FK catalogs fail to follow the predictions of mathematical statistics -- only if the FK3 and the FK4 contain higher angular frequency systematic errors. These errors, at the mean epochs of place of the older catalogs, are larger than the random errors in the FK5 system. We expect that our new ability to discern these features means that we can also prevent them during compilation. Finally, note that the concepts behind our procedures are completely general and can be applied to any set of compiled data -- astrometric, photometric, geodetic, and so forth.

Bucciarelli, B.↗
NASA NTRS1994

Topography and tectonics of the central New Madrid seismic zone: Results of numerical experiements using a three-dimensional boundary element program

We present results of a series of numerical experiments designed to test hypothetical mechanisms that derive deformation in the New Madrid seismic zone. Experiments are constrained by subtle topography and the distribution of seismicity in the region. We use a new boundary element algorithm that permits calcuation of the three-dimensional deformation field. Surface displacement fields are calculated for the New Madrid zone under both far-field (plate tectonics scale) and locally derived driving strains. Results demonstrate that surface displacement fields cannot distinguish between either a far-field simple or pure shear strain field or one that involves a deep shear zone beneath the upper crustal faults. Thus, neither geomorphic nor geodetic studies alone are expected to reveal the ultimate driving mechanism behind the present-day deformation. We have also tested hypotheses about strain accommodation within the New Madrid contractional step-over by including linking faults, two southwest dipping and one vertical, recently inferred from microearthquake data. Only those models with step-over faults are able to predict the observed topography. Surface displacement fields for long-term, relaxed deformation predict the distribution of uplift and subsidence in the contractional step-over remarkably well. Generation of these displacement fields appear to require slip on both the two northeast trending vertical faults and the two dipping faults in the step-over region, with very minor displacements occurring during the interseismic period when the northeast trending vertical faults are locked. These models suggest that the gently dippling central step-over fault is a reverse fault and that the steeper fault, extending to the southeast of the step-over, acts as a normal fault over the long term.

Gomberg, Joan↗
NASA NTRS1995

Relative motions of the Australian, Pacific and Antarctic plates estimated by the Global Positioning System

Global Positioning System (GPS) measurements spanning approximately 3 years have been used to determine velocities for 7 sites on the Australian, Pacific and Antarctic plates. The site velocities agree with both plate model predictions and other space geodetic techniques. We find no evidence for internal deformation of the interior of the Australian plate. Wellington, New Zealand, located in the Australian-Pacific plate boundary zone, moves 20 +/- 5 mm/yr west-southwest relative to the Australian plate. Its velocity lies midway between the predicted velocities of the two plates. Relative Euler vectors for the Australia-Antarctica and Pacific-Antarctica plates agree within one standard deviation with the NUVEL-1A predictions.

Larson, Kristine M.↗
NASA NTRS1994

Determination of the extragalactic-planetary frame tie from joint analysis of radio interferometric and lunar laser ranging measurements

Very Long Baseline Interferometry (VLBI) observations of extragalactic radio sources provide the basis for defining an accurate non-rotating reference frame in terms of angular positions of the sources. Measurements of the distance from the Earth to the Moon and to the inner planets provide the basis for defining an inertial planetary ephemeris reference frame. The relative orientation, or frame tie, between these two reference frames is of interest for combining Earth orientation measurements, for comparing Earth orientation results with theories referred to the mean equator and equinox, and for determining the positions of the planets with respect to the extragalactic reference frame. This work presents an indirect determination of the extragalactic-planetary frame tie from a combined reduction of VLBI and Lunar Laser Ranging (LLR) observations. For this determination, data acquired by LLR tracking stations since 1969 have been analyzed and combined with 14 years of VLBI data acquired by NASA's Deep Space Network since 1978. The frame tie derived from this joint analysis, with an accuracy of 0.003 sec, is the most accurate determination obtained so far. This result, combined with a determination of the mean ecliptic (defined in the rotating sense), shows that the mean equinox of epoch J2000 is offset from the x-axis of the extragalactic frame adopted by the International Earth Rotation Service for astrometric and geodetic applications by 0.078 sec +/- 0.010 sec along the y-direction and y 0.019 sec +/- 0.001 sec. along the z-direction.

Folkner, W. M.↗
NASA NTRS1994

Atmospheric pressure loading parameters from very long baseline interferometry observations

Atmospheric mass loading produces a primarily vertical displacement of the Earth's crust. This displacement is correlated with surface pressure and is large enough to be detected by very long baseline interferometry (VLBI) measurements. Using the measured surface pressure at VLBI stations, we have estimated the atmospheric loading term for each station location directly from VLBI data acquired from 1979 to 1992. Our estimates of the vertical sensitivity to change in pressure range from 0 to -0.6 mm/mbar depending on the station. These estimates agree with inverted barometer model calculations (Manabe et al., 1991; vanDam and Herring, 1994) of the vertical displacement sensitivity computed by convolving actual pressure distributions with loading Green's functions. The pressure sensitivity tends to be smaller for stations near the coast, which is consistent with the inverted barometer hypothesis. Applying this estimated pressure loading correction in standard VLBI geodetic analysis improves the repeatability of estimated lengths of 25 out of 37 baselines that were measured at least 50 times. In a root-sum-square (rss) sense, the improvement generally increases with baseline length at a rate of about 0.3 to 0.6 ppb depending on whether the baseline stations are close to the coast. For the 5998-km baseline from Westford, Massachusetts, to Wettzell, Germany, the rss improvement is about 3.6 mm out of 11.0 mm. The average rss reduction of the vertical scatter for inland stations ranges from 2.7 to 5.4 mm.

Macmillan, D. S.↗
NASA NTRS1995

A simulation to study the feasibility of improving the temporal resolution of LAGEOS geodynamic solutions by using a sequential process noise filter

A key drawback to estimating geodetic and geodynamic parameters over time based on satellite laser ranging (SLR) observations is the inability to accurately model all the forces acting on the satellite. Errors associated with the observations and the measurement model can detract from the estimates as well. These 'model errors' corrupt the solutions obtained from the satellite orbit determination process. Dynamical models for satellite motion utilize known geophysical parameters to mathematically detail the forces acting on the satellite. However, these parameters, while estimated as constants, vary over time. These temporal variations must be accounted for in some fashion to maintain meaningful solutions. The primary goal of this study is to analyze the feasibility of using a sequential process noise filter for estimating geodynamic parameters over time from the Laser Geodynamics Satellite (LAGEOS) SLR data. This evaluation is achieved by first simulating a sequence of realistic LAGEOS laser ranging observations. These observations are generated using models with known temporal variations in several geodynamic parameters (along track drag and the J(sub 2), J(sub 3), J(sub 4), and J(sub 5) geopotential coefficients). A standard (non-stochastic) filter and a stochastic process noise filter are then utilized to estimate the model parameters from the simulated observations. The standard non-stochastic filter estimates these parameters as constants over consecutive fixed time intervals. Thus, the resulting solutions contain constant estimates of parameters that vary in time which limits the temporal resolution and accuracy of the solution. The stochastic process noise filter estimates these parameters as correlated process noise variables. As a result, the stochastic process noise filter has the potential to estimate the temporal variations more accurately since the constraint of estimating the parameters as constants is eliminated. A comparison of the temporal resolution of solutions obtained from standard sequential filtering methods and process noise sequential filtering methods shows that the accuracy is significantly improved using process noise. The results show that the positional accuracy of the orbit is improved as well. The temporal resolution of the resulting solutions are detailed, and conclusions drawn about the results. Benefits and drawbacks of using process noise filtering in this type of scenario are also identified.

Hartman, Brian Davis↗
NASA NTRS1995

Dynamics of the solid earth

The work performed under this contract provides for the improvement of modeling of troposphere propagation delay in very long base interferometry (VLBI) data, defining of required accuracy of ancillary data type(s) to measure tropospheric behavior and to develop a clear correspondence between tropospheric behavior and the quality of geodetic VLBI measurements.

Knight, Curtis A.↗
NASA NTRS1995

An a priori model for the reduction of nutation observations: KSV(1994.3) nutation series

We discuss the formulation of a new nutation series to be used in the reduction of modern space geodetic data. The motivation for developing such a series is to develop a nutation series that has smaller short period errors than the IAU 1980 nutation series and to provide a series that can be used with techniques such as the Global Positioning System (GPS) that have sensitivity to nutations but can directly separate the effects of nutations from errors in the dynamical force models that effect the satellite orbits. A modern nutation series should allow the errors in the force models for GPS to be better understood. The series is constructed by convolving the Kinoshita and Souchay rigid Earth nutation series with an Earth response function whose parameters are partly based on geophysical models of the Earth and partly estimated from a long series (1979-1993) of very long baseline interferometry (VLBI) estimates of nutation angles. Secular rates of change of the nutation angles to represent corrections to the precession constant and a secular change of the obliquity of the ecliptic are included in the theory. Time dependent amplitudes of the Free Core Nutation (FCN) that is most likely excited by variations in atmospheric pressure are included when the geophysical parameters are estimated. The complex components of the prograde annual nutation are estimated simultaneously with the geophysical parameters because of the large contribution to the nutation from the S(sub 1) atmospheric tide. The weighted root mean square (WRMS) scatter of the nutation angle estimates about this new model are 0.32 mas and the largest correction to the series when the amplitudes of the ten largest nutations are estimated is 0.18 +/- 0.03 mas for the in phase component of the prograde 18. 6 year nutation.

Herring, T. A.↗