Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tracking errors”

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,783 records · Page 99
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Improvements to Ozone Analyses Using Hyperspectral Sounders in the 9.6 μm Band

Previously, hyperspectral sounder brightness temperatures assimilated in the Goddard Earth Observing System Atmospheric Data Assimilation System (GEOS-ADAS) were limited to assimilating temperature and moisture. The ozone sensitive 9.6 m region is sensed by several hyperspectral sounders including AIRS (Atmospheric InfraRed Sounder), IASI (Infrared Atmospheric Sounding Interferometer), and CrIS (Cross-track Infrared Sounder). Direct assimilation of brightness temperatures in the 9.6 m region have been operational at ECMWF for several years (Dragani and McNally, 2013; Eresmaa et al., 2017). With this study, similar improvements using the GEOS-ADAS are presented. Channels were selected from available operational subsets evaluating information content and minimizing inter-channel correlation. Additionally, information such as channel selections made by other studies, and vertical sensitivities of ozone and temperature were considered in developing the study. The analyses produced show improvements verified against ozonesondes taken from SHADOZ (Southern Hemisphere Additional Ozonesondes), and WOUDC (World Ozone and Ultraviolet Data Center). Inclusion of inter-channel correlation error to the selected ozone channels are expected to improve ozone analyses further along with improve analyses overall. It is anticipated that inclusion of these ozone sensitive channels will be used to improve NASA GMAO products in the near future.

GEOS-ADAS↗
NASA NTRS1985

Computer simulation of TDRSS autotracking

The Tracking and Data Relay Satellite System (TDRSS) consists of synchronous-orbit satellites and a single ground station at White Sands, NM. Each TDRSS satellite has two 16-foot dish antennas with dual-frequency feeds for transmission and reception at S- and Ku-band. Autotracking is required on the Ku-band return link from user spacecraft. The TDRSS autotrack system was verified by computer simulation rather than by end-to-end hardware testing. One of the simulation models is described. Attention is given to an overview of system operation, a computer model description, and the simulation results. An example plot shows the sensitivity of system pointing error to user incidental amplitude modulation frequency.

Jones, S. C.↗
NASA NTRS1998

Assessment Study of the State of the Art in Adaptive Control and its Applications to Aircraft Control

Many papers relevant to reconfigurable flight control have appeared over the past fifteen years. In general these have consisted of theoretical issues, simulation experiments, and in some cases, actual flight tests. Results indicate that reconfiguration of flight controls is certainly feasible for a wide class of failures. However many of the proposed procedures although quite attractive, need further analytical and experimental studies for meaningful validation. Many procedures assume the availability of failure detection and identification logic that will supply adequately fast, the dynamics corresponding to the failed aircraft. This in general implies that the failure detection and fault identification logic must have access to all possible anticipated faults and the corresponding dynamical equations of motion. Unless some sort of explicit on line parameter identification is included, the computational demands could possibly be too excessive. This suggests the need for some form of adaptive control, either by itself as the prime procedure for control reconfiguration or in conjunction with the failure detection logic. If explicit or indirect adaptive control is used, then it is important that the identified models be such that the corresponding computed controls deliver adequate performance to the actual aircraft. Unknown changes in trim should be modelled, and parameter identification needs to be adequately insensitive to noise and at the same time capable of tracking abrupt changes. If however, both failure detection and system parameter identification turn out to be too time consuming in an emergency situation, then the concepts of direct adaptive control should be considered. If direct model reference adaptive control is to be used (on a linear model) with stability assurances, then a positive real or passivity condition needs to be satisfied for all possible configurations. This condition is often satisfied with a feedforward compensator around the plant. This compensator must be robustly designed such that the compensated plant satisfies the required positive real conditions over all expected parameter values. Furthermore, with the feedforward only around the plant, a nonzero (but bounded error) will exist in steady state between the plant and model outputs. This error can be removed by placing the compensator also in the reference model. Design of such a compensator should not be too difficult a problem since for flight control it is generally possible to feedback all the system states.

Kaufman, Howard↗
NASA NTRSDate not supplied

Cartesian Mesh Simulations for the Third AIAA Sonic Boom Prediction Workshop

Simulation results are presented for all cases from the Third AIAA Sonic Boom Prediction Workshop. An inviscid, embedded-boundary Cartesian-mesh flow solver is used in conjunction with adjoint-based mesh adaptation to compute nearfield pressure signatures. Specialized techniques are applied to maximize accuracy and minimize cost on Cartesian meshes. The Richardson-based error estimate highlights regions of the signatures most sensitive to mesh refinement. Timing results and coarse, medium, and fine mesh sizes for nearfield cases demonstrate that the parallel decomposition approach is efficient in both computational time and wall-clock. Pressure signals are propagated to the ground using an augmented Burgers’ equation solver to predict boom carpets. Ground signatures and loudness metrics are presented for a standard atmosphere as well as more realistic atmospheric profiles, which affect overall noise levels and can significantly widen the boom carpet. Mesh convergence studies show that high sampling frequencies, around 500 kHz, are required for propagation, and the sampling frequency increases at large off-track angles with longer acoustic ray paths and propagation times. The numerical methods yield accurate results for predicting low sonic boom signatures while being among the least computationally expensive of the workshop.

ARMD↗
NASA NTRSDate not supplied

NASA’s Identified Risks of Adverse Outcomes Due to Inadequate Human Systems Integration Architecture in Human Spaceflight

The NASA Human System Risk Board (HSRB) has the overall responsibility for tracking the evolution of the top ~30 human system risks that it has identified to be associated with human spaceflight. As part of this process, the Board is charged with maintaining a consistent, integrated process to mitigate those risks, and developing evidence-based risk posture recommendations. One of the identified risks is due to inadequate human systems integration architecture (HSIA) and a driving factor of this risk is that given decreasing real-time ground support for execution of complex operations during future exploration missions, there is a possibility of adverse performance outcomes including that crew are unable to adequately respond to unanticipated critical malfunctions or detect safety critical procedural errors. The HSRB uses Directed Acyclic Graphs (DAGs) as a communication tool for describing how astronaut exposure to spaceflight hazards leads to meaningful mission-level health and performance outcomes and as the basis for understanding intermediate causal relationships between risk contributing factors and countermeasures that link hazards to outcomes. The HSIA risk DAG will be presented and described. Historically, critical malfunctions requiring Crew/MCC management occurred at a rate of 1.7 times per year for ISS averaged over the lifetime and 3-4 times per year in the burn in phase for the vehicle. These averages do not include EVA data, which greatly increases the incident rate. Prior experience from the Apollo program showed 10/11 crewed missions experienced significant anomalies where crew relied heavily on MCC expertise in real-time. These failure patterns are in line with those observed in other complex engineered systems (e.g., oil rigs, launch systems, commercial aviation, etc.) It is likely that general malfunction and error rates are > 10% for short duration missions (<30 days), based on past and current spaceflight operations data. Likelihood of adverse outcomes has the potential to increase as crew conduct work with new, complex systems and with less ground support. For Low Earth Orbit missions and Lunar missions less than 30 days, assuming minimal comm delays, disruptions and bandwidth limitations, malfunctions and errors can affect mission objectives and crew health but may be mitigated by ground support. For Lunar missions greater than 30 days and any potential Mars mission malfunctions and errors can have Loss of Crew and Loss of Mission consequences due to reduced ground support (communication delays, constraints and blackouts) for more complex operations, as well as reduced resupply and evacuation options.

Daniel M Buckland↗
NASA NTRS1972

Error analysis of earth physics satellite systems

Error analysis of distant-satellite-to-close-satellite range-rate, satellite-to-sea altimetry, and ground station to satellite range are made by simulations in which observational variances are assumed, observation equations are formed, and normal equations incremented. The final normal equation matrix is inverted to obtain standard deviations and correlation coefficients. The natural parameters solved for are the broad variations of the gravity field, represented by harmonic coefficients; local variations of gravity, represented by point masses; and the departure of the sea level from the geoid, represented by area means. A standard case of a low (263 km) polar close satellite, three equatorial geosynchronous satellites, and eight ground tracking stations is set up.

Kaula, W. M.↗
NASA NTRS2000

Geopotential Model Improvement Using POCM_4B Dynamic Ocean Topography Information: PGM2000A

The two-year mean (1993-1994) Dynamic Ocean Topography (DOT) field implied by the POCM_4B circulation model was used to develop normal equations for DOT, in a surface spherical harmonic representation. These normal equations were combined with normal equations from satellite tracking data, surface gravity data, and altimeter data from TOPEX/Poseidon and ERS-1. Several least-squares combination solutions were developed in this fashion, by varying parameters such as the maximum degree of the estimated DOT and the relative weights of the different data. The solutions were evaluated in terms of orbit fit residuals, GPS/Leveling-derived undulations, and independent DOT information from in situ WOCE hydrographic data. An optimal solution was developed in this fashion which was originally presented at the 1998 EGS meeting in Nice, France. This model, designated here PGM2000A, maintains the orbit and land geoid modeling performance of EGM96, while improving its marine geoid modeling capability. In addition, PGM2000A's error spectrum is considerably more realistic than those of other contemporary gravitational models and agrees well with the error spectrum of EGM96. We will present the development and evaluation of PGM2000A, with particular emphasis on the weighting of the DOT information implied by POCM_4B. We will also present an inter-comparison of PGM2000A with the GRIM5-C1 and TEG-4 models. Directions for future work and problematic areas will be identified.

Pavlis, N. K.↗
NASA NTRS2012

Locating the LCROSS Impact Craters

The Lunar CRater Observations and Sensing Satellite (LCROSS) mission impacted a spent Centaur rocket stage into a permanently shadowed region near the lunar south pole. The Sheperding Spacecraft (SSC) separated approx. 9 hours before impact and performed a small braking maneuver in order to observe the Centaur impact plume, looking for evidence of water and other volatiles, before impacting itself. This paper describes the registration of imagery of the LCROSS impact region from the mid- and near-infrared cameras onboard the SSC, as well as from the Goldstone radar. We compare the Centaur impact features, positively identified in the first two, and with a consistent feature in the third, which are interpreted as a 20 m diameter crater surrounded by a 160 m diameter ejecta region. The images are registered to Lunar Reconnaisance Orbiter (LRO) topographical data which allows determination of the impact location. This location is compared with the impact location derived from ground-based tracking and propagation of the spacecraft's trajectory and with locations derived from two hybrid imagery/trajectory methods. The four methods give a weighted average Centaur impact location of −84.6796 deg, −48.7093 deg, with a 1σ uncertainty of 115 m along latitude, and 44 m along longitude, just 146 m from the target impact site. Meanwhile, the trajectory-derived SSC impact location is −84.719 deg, −49.61 deg, with a 1 alpha uncertainty of 3 m along the Earth vector and 75 m orthogonal to that, 766 m from the target location and 2.803 km south-west of the Centaur impact. We also detail the Centaur impact angle and SSC instrument pointing errors. Six high-level LCROSS mission requirements are shown to be met by wide margins. We hope that these results facilitate further analyses of the LCROSS experiment data and follow-up observations of the impact region

LCROSS↗
NASA NTRS2016

Isokinetic TWC Evaporator Probe: Development of the IKP2 and Performance Testing for the HAIC-HIWC Darwin 2014 and Cayenne 2015 Field Campaigns

A new Isokinetic Total Water Content Evaporator (IKP2) was downsized from a prototype instrument, specifically to make airborne measurements of hydrometeor total water content (TWC) in deep tropical convective clouds to assess the new ice crystal Appendix D icing envelope. The probe underwent numerous laboratory and wind tunnel investigations to ensure reliable operation under the difficult high altitude/speed/TWC conditions under which other TWC instruments have been known to either fail, or have unknown performance characteristics. The article tracks the testing and modifications of the IKP2 probe to ensure its readiness for three flight campaigns in 2014 and 2015. Comparisons are made between the IKP2 and the NASA Icing Research Tunnel reference values in liquid conditions, and to an exploratory technique estimating ice water content from a bulk ice capture cylinder method in glaciated conditions. These comparisons suggest that the initial target of 20% accuracy in TWC has been achieved and likely exceeded for tested TWC values in excess of about 0.5/cu gm. Uncertainties in the ice water content reference method have been identified. Complications are introduced in the necessary subtraction of an independently measured background water vapor concentration, errors of which are small at the colder flight temperatures, but increase rapidly with increasing temperature, and ultimately limit the practical use of the instrument in a tropical convective atmosphere to conditions colder than about 0 C. A companion article in this conference traces the accuracy of the components of the IKP2 to derive estimated system accuracy.

Meteorological Instruments↗
NASA NTRS2019

Gravitational Changes of the Earth's Free Oscillation From Earthquakes: Theory and Feasibility Study Using GRACE Inter‐Satellite Tracking

GRACE satellites have detected regional‐scale preseismic, coseismic, and postseismic gravity changes associated with great earthquakes during the GRACE era (2002‐2017). Earthquakes also excite global‐scale transient gravity changes associated with free oscillations that may be discerned for a few days. In this study, we examine such global gravity changes due to Earth's free oscillations and quantify how they affect GRACE measurements. We employ the normal mode formalism to synthesize the global gravity changes after the 2004 Sumatra earthquake and simulate the (gravitational) free oscillation signals manifested in the GRACE K‐band ranging (KBR) measurements. Using the Kaula orbit perturbation theory, we show how GRACE inter‐satellite distances are perturbed through a complex coupling of eigenfrequencies of the normal modes with the Earth's rotation rate and the GRACE satellites' orbital frequency. It is found that a few gravest normal modes can generate range‐rate perturbations as large as 0.2 μm/s, which are comparable to actual errors of GRACE KBR ranging and accelerometer instruments. Wavelet time‐frequency analysis of the GRACE KBR residual data in December 2004 reveals the existence of a significant transient signal after the 2004 Sumatra earthquake. This transient signal is characterized by a frequency of ~0.022 mHz that could be potentially associated with the largest excitation due to the “football” mode of the Earth's free oscillation. However, the results are also affected by low‐frequency noise of the GRACE accelerometers. Improved space‐borne gravitational instrumentation may open new opportunities to study the Earth's interior and earthquakes independently from global seismological analysis.

Ghobadi-Far, Khosro↗
NASA NTRS1982

Starprobe - Coronal plasma turbulence effects on tracking and telemetry

The telemetry and tracking problems expected for Starprobe, a solar flyby with a perihelion of 4 solar radii, are discussed. The Starprobe is intended to obtain measurements of the solar quadrupole moment within an accuracy of 1/100,000,000, allowing more accurate modelling of the solar interior. Obstacles are perceived for Doppler velocity measurements due to RF signal amplitude and phase scintillations created by the solar corona. Model calculations are presented for the propagation effects in the maximum disturbance area, and an open loop two-way Doppler system is described. The system comprises X-band and S-band uplinks and downlinks to eliminate the plasma-induced phase scintillations, means of correcting unfolding errors caused by S-band amplitude scintillations, and a new method of estimating the Doppler frequency.

Armstrong, J. W.↗
NASA NTRS1999

Orbit Determination for Mars Global Surveyor During Mapping

The Mars Global Surveyor (MGS) spacecraft reached a low-altitude circular orbit on February 4, 1999, after the termination of the second phase of aerobraking. The MGS spacecraft carries the Mars Orbiter Laser Altimeter (MOLA) whose primary goal is to derive a global, geodetically referenced 0.2 deg x 0.2 deg topographic grid of Mars with a vertical accuracy of better than 30 meters. During the interim science orbits in the' Hiatus mission phase (October - November 1997), and the Science Phasing Orbits (March - April, 1998, and June - July 1998) 208 passes of altimeter data were collected by the MOLA instrument. On March 1, 1999 the first ten orbits of MOLA altimeter data from the near-circular orbit were successfully returned from MGS by the Deep Space Network (DSN). Data will be collected from MOLA throughout the Mapping phase of the MCS mission, or for at least one Mars year (687 days). Whereas the interim orbits of Hiatus and SPO were highly eccentric, and altimeter data were only collected near periapsis when the spacecraft was below 785 km, the Mapping orbit of MGS is near circular, and altimeter data will be collected continuously at a rate of 10 Hz. The proper analysis of the altimeter data requires that the orbit of the MGS spacecraft be known to an accuracy comparable to that of the quality of the altimeter data. The altimeter has an ultimate precision of 30 cm on mostly flat surfaces, so ideally the orbits of the MGS spacecraft should be known to this level. This is a stringent requirement, and more realistic goals of orbit error for MGS are ten to thirty meters. In this paper we will discuss the force and measurement modelling required to achieve this objective. Issues in force modelling include the proper modelling of the gravity field of Mars, and the modelling of non-conservatives forces, including the development of a 'macro-model', in a similar fashion to TOPEX/POSEIDON and TDRSS. During Cruise and Aerobraking, the high gain antenna (HGA) was stowed on the +X face of the spacecraft. On March 29, 1999 the HGA will be deployed on a meter long boom which will remain Earth-pointing while the instrument panel (including the MOLA instrument) remains pointed at nadir. The tracking data must be corrected for the regular motion of the high gain antenna with respect to the center of mass, and the success of the MGS determination during Mapping will depend on correctly accounting for this offset in the measurement model.

Lemoine, F. G.↗
NASA NTRS2017

Monitoring Land Surface Albedo and Vegetation Dynamics Using High Spatial and Temporal Resolution Synthetic Time Series from Landsat and the MODIS BRDF/NBAR/Albedo Product

Seasonal vegetation phenology can significantly alter surface albedo which in turn affects the global energy balance and the albedo warmingcooling feedbacks that impact climate change. To monitor and quantify the surface dynamics of heterogeneous landscapes, high temporal and spatial resolution synthetic time series of albedo and the enhanced vegetation index (EVI) were generated from the 500-meter Moderate Resolution Imaging Spectroradiometer (MODIS) operational Collection V006 daily BRDF (Bidirectional Reflectance Distribution Function) / NBAR (Nadir BRDF-Adjusted Reflectance) / albedo products and 30-meter Landsat 5 albedo and near-nadir reflectance data through the use of the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM). The traditional Landsat Albedo (Shuai et al., 2011) makes use of the MODIS BRDFAlbedo products (MCD43) by assigning appropriate BRDFs from coincident MODIS products to each Landsat image to generate a 30-meter Landsat albedo product for that acquisition date. The available cloud free Landsat 5 albedos (due to clouds, generated every 16 days at best) were used in conjunction with the daily MODIS albedos to determine the appropriate 30-meter albedos for the intervening daily time steps in this study. These enhanced daily 30-meter spatial resolution synthetic time series were then used to track albedo and vegetation phenology dynamics over three Ameriflux tower sites (Harvard Forest in 2007, Santa Rita in 2011 and Walker Branch in 2005). These Ameriflux sites were chosen as they are all quite nearby new towers coming on line for the National Ecological Observatory Network (NEON), and thus represent locations which will be served by spatially paired albedo measures in the near future. The availability of data from the NEON towers will greatly expand the sources of tower albedometer data available for evaluation of satellite products. At these three Ameriflux tower sites the synthetic time series of broadband shortwave albedos were evaluated using the tower albedo measurements with a Root Mean Square Error (RMSE) less than 0.013 and a bias within the range of 0.006. These synthetic time series provide much greater spatial detail than the 500 meter gridded MODIS data, especially over more heterogeneous surfaces, which improves the efforts to characterize and monitor the spatial variation across species and communities. The mean of the difference between maximum and minimum synthetic time series of albedo within the MODIS pixels over a subset of satellite data of Harvard Forest (16 kilometers by 14 kilometers) was as high as 0.2 during the snow-covered period and reduced to around 0.1 during the snow-free period. Similarly, we have used STARFM to also couple MODIS Nadir BRDF-Adjusted Reflectances (NBAR) values with Landsat 5 reflectances to generate daily synthetic times series of NBAR and thus Enhanced Vegetation Index (NBAR-EVI) at a 30-meter resolution. While normally STARFM is used with directional reflectances, the use of the view angle corrected daily MODIS NBAR values will provide more consistent time series. These synthetic times series of EVI are shown to capture seasonal vegetation dynamics with finer spatial and temporal details, especially over heterogeneous land surfaces.

Vegetation Index↗
NASA NTRS1992

The TOPSAR interferometric radar topographic mapping instrument

The NASA DC-8 AIRSAR instrument was augmented with a pair of C-band antennas displaced across track to form an interferometer sensitive to topographic variations of the Earth's surface. The antennas were developed by the Italian consortium Co.Ri.S.T.A., under contract to the Italian Space Agency (ASI), while the AIRSAR instrument and modifications to it supporting TOPSAR were sponsored by NASA. A new data processor was developed at JPL for producing the topographic maps, and a second processor was developed at Co.Ri.S.T.A. All the results presented below were processed at JPL. During the 1991 DC-8 flight campaign, data were acquired over several sites in the United States and Europe, and topographic maps were produced from several of these flight lines. Analysis of the results indicate that statistical errors are in the 2-3 m range for flat terrain and in the 4-5 m range for mountainous areas.

Zebker, Howard A.↗
NASA NTRS2000

Calibrating Atmospheric Delay for the Cassini Gravitational Wave Experiment

The Cassini spacecraft together with one of the stations of the Deep Space Network (DSN) have been instrumented to carry out extremely precise Doppler tracking that will be used to search for the direct evidence of Gravitational radiation passing through our solar system. The two-way communications link between the ground antenna and the spacecraft constitute an "antenna" for gravitational waves that would perturb the phase of the RF signal in the link. The experiments will be carried out during the long cruise phase of the spacecraft on its journey to Saturn and will be sensitive to gravitational wave perturbations larger than the noise level fluctuations of 3 x 10(exp -15) as measured in the Allan Standard Deviation domain. The use of simultaneous, coherent X- and Ka-Band up and down signals will reduce the errors associated with charged particle fluctuations in the interplanetary medium and Earth's ionosphere to a negligible level. The primary fluctuations in the phase of the signals both to and from the spacecraft are expected to be caused by fluctuations in water vapor in the Earth's atmosphere We have designed and are testing a new atmospheric calibration system (with duplicate components) to sense line-of-sight water vapor and its physical temperature with a goal of calibrating 95% or more of tropospheric path delay fluctuations during the Cassini Gravitational Wave Experiment (GWE). The critical component of the calibration system consists of a newly designed water vapor radiometer having a I degree sensing beamwidth and 0.01 K brightness temperature stability over hour time scales. Auxiliary instrumentation includes a microwave temperature profiler to retrieve the vertical distribution of the vapor physical temperature, and surface meteorology. A detailed error budget has been developed to account for all of the possible sources of error during calibration of the GWE and will be discussed. Two identical calibration systems have been constructed in order to provide capability as well as backup during the actual GWE experiment. We will report on side-by-side testing of the vapor and temperature sensing components of this calibrations system as well as a comparison with a short baseline radio interferometric measurement at our Goldstone complex.

Resch, G. M.↗
NASA NTRS2016

Isokinetic TWC Evaporator Probe: Development of the IKP2 and Performance Testing for the HAIC-HIWC Darwin 2014 and Cayenne-2015 Field Campaigns

A new Isokinetic Total Water Content Evaporator (IKP2) was downsized from a prototype instrument, specifically to make airborne measurements of hydrometeor total water content (TWC) in deep tropical convective clouds to assess the new ice crystal Appendix D icing envelope. The probe underwent numerous laboratory and wind tunnel investigations to ensure reliable operation under the difficult high altitude/speed/TWC conditions under which other TWC instruments have been known to either fail, or have unknown performance characteristics. The article tracks the testing and modifications of the IKP2 probe to ensure its readiness for three flight campaigns in 2014 and 2015. Comparisons are made between the IKP2 and the NASA Icing Research Tunnel reference values in liquid conditions, and to an exploratory technique estimating ice water content from a bulk ice capture cylinder method in glaciated conditions. These comparisons suggest that the initial target of 20 percent accuracy in TWC has been achieved and likely exceeded for tested TWC values in excess of about 0.5 gm (exp -3). Uncertainties in the ice water content reference method have been identified. Complications are introduced in the necessary subtraction of an independently measured background water vapour concentration, errors of which are small at the colder flight temperatures, but increase rapidly with increasing temperature, and ultimately limit the practical use of the instrument in a tropical convective atmosphere to conditions colder than about 0 degrees C. A companion article in this conference traces the accuracy of the components of the IKP2 to derive estimated system accuracy.

Meteorological Instruments↗
NASA NTRS1990

The application of Encke's method to long arc orbit determination solutions

The Laser Geodynamics Satellite (LAGEOS) was launched on May 4, 1976 to provide geophysical measurements by means of laser ranging techniques. To date, over twelve years of laser range measurements have been collected from various tracking stations located around the world. Laser range measurements to LAGEOS have contributed to studies of earth rotation, plate tectonics, global baseline, and the gravity field as well as many other areas. Some of these studies are based upon the determination of a single, continuous orbit for LAGEOS for time spans on the order of several years. Current studies at the University of Texas Center for Space Research include the precision orbit determination of LAGEOS for arc lengths of up to 12.8 years which represents over 31,000 orbital revolutions. These long arc studies have led to the implementation of Encke's method to improve the convergence of the batch filter while reducing numerical integration errors. While the technique has been successfully applied to arc lengths of up to 12.8 years, the results presented focus on the solution of a six-year arc.

Lundberg, J. B.↗
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↗