Search NASASearch

Engineering topics

Cox, C. M.

Publications and source records attributed to Cox, C. M..

Excitation of Earth Rotation Variations "Observed" by Time-Variable Gravity

Time variable gravity measurements have been made over the past two decades using the space geodetic technique of satellite laser ranging, and more recently by the GRACE satellite mission with improved spatial resolutions. The degree-2 harmonic components of the time-variable gravity contain important information about the Earth s length-of-day and polar motion excitation functions, in a way independent to the traditional "direct" Earth rotation measurements made by, for example, the very-long-baseline interferometry and GPS. In particular, the (degree=2, order= 1) components give the mass term of the polar motion excitation; the (2,O) component, under certain mass conservation conditions, gives the mass term of the length-of-day excitation. Combining these with yet another independent source of angular momentum estimation calculated from global geophysical fluid models (for example the atmospheric angular momentum, in both mass and motion terms), in principle can lead to new insights into the dynamics, particularly the role or the lack thereof of the cores, in the excitation processes of the Earth rotation variations.

Chao, Ben F.

The Use of Mascons to Resolve Time-Variable Gravity from GRACE

We have analyzed GRACE Level 1B data to resolve time-variable gravity using a local mascon approach. The spherical harmonic solutions released to date resolve the signal from surface hydrology over land areas at spatial scales of 750 to 1000 km over one month intervals [Wahr et al., 2004; Tapley et al., 2004]. In our local approach, we solve explicitly for the mass of water in surface blocks using only the KBRR data collected as GRACE overflies the region of interest. The local representation of gravity minimizes leakage of errors from other areas due to aliasing or mismodelling. In this paper, we report on the analysis of GRACE data from January 2003 through August 2004 over three regions: the Amazon, the Indian subcontinent, and the continental United States. We solve for mass change at 10-day intervals using 4 deg x 4 deg blocks. We give an overview of our latest results, and we present the results of error analyses, and comparisons to both hydrology models and in-situ data.

Lemoine, Frank G.

Unique Approaches to Analysis of Time-Variable Gravity from GRACE

We have developed an innovative analysis strategy for analysis of GRACE data. We have developed a capability to recover local/regional gravity changes using non-global functional representations (Le. surface anomalies vs. global spherical harmonics) h m the GRACE data. Our approach can take regularly or irregularly shaped regions, populate them with surface anomaly blocks of suitable area and solve for the resulting mass flux with respect to a mean field. The surface mass or gravity anomalies benefit from the application of spatial and temporal constraints to add stability to the solution. In this paper we discuss the analysis of four months of GRACE Level 1B data (accelerometry, intersatellite data, attitude information and precise orbits) from July to October 2003, recently released to the GRACE Science Team. We compare and contrast this local approach to gravity recovery, with the more conventional approach using global spherical harmonics. We review simulations of this technique which allow us to pinpoint optimum strategies for applications of this local gravity recovery approach.

Lemoine, Frank G.

Time-Variable Gravity Signal due to Extratropic Pacific Water Mass Redistribution

Using the satellite-laser-ranging (SLR) data, Cox and Chao [2002] reported the detection of a large post-1998 anomaly (in the form of a positive jump) in the time series of Earth s lowest-degree gravity harmonic 52, or the dynamic oblateness. Among several groups now examining the mass redistribution in the global geophysical fluids in search of the cause(s), we report here a temporally coinciding anomalies found in the extratropic north + south Pacific basins. Clearly seen in the leading EOFPC mode for extratropic Pacific, these anomalies occurred in sea-surface height, sea-surface temperature, and temperature- and salinity-depth profiles. We based our analysis on two different data sources: TOPEX/Poseidon altimetry, and the ECCO ocean general circulation model output assimilating T/P data. The magnitude of these changes, when converted to equivalent J2 change, appears to be a few times too small to explain the observed J2 directly. These findings, and the fact that the anomalies occurred following the strong 1997-98 El Nino, suggest strong geophysical connection of the interannual-to-decadal variation of 52 with the Pacific Decadal Oscillation (PDO) and the ultimate global-change processes that cause PDO. More work is underway, and additional independent data sources are examined, paying close attention to the fact that the J2 anomaly has been reversing back to normal since 2001. These include: (1) cryospheric contributions (melting of glaciers and ice sheets); (2) land hydrological contributions; (3) polar sea influences ( e g , via deep flow); (4) fluid flow in Earth's core; (5) time-variable gravity signals from SLR in higher harmonic degree/order, including J3,J4, (2,1), and (2,2) coefficients, considering their lower signal-to-noise ratios; (6) Earth rotation data in terms of length-of-day and polar motion.

Chao, B. F.

Gravity Models from CHAMP and other Satellite Data

CHAMP spacecraft is the first of a series of new spacecraft missions that are revolutionizing our ability to model the Earth's geopotential. We report on the analysis of over 100 days of CHAMP data in 2001 and 2002, merged with tracking data of other satellites such as Jason, Topex, GFO, Starlette, Stella, Spot-2, as well as satellite altimetry. We find that the CHAMP-only component of these solutions is a significant improvement over pre-CHAMP satellite only models with respect to the high degree information expressed by the geopotential model coefficients. For example, the variance of the differences with altimeter-derived anomalies through degree 70 is 2.80 mGal(sup 2) for the CHAMP-only solution based on 87 days of data vs. 10.19 mGal(sup 2) for EGM96S. Nonetheless, in order to model properly the various resonances to which different satellites are sensitive, we must include other satellite data. We evaluate the performance of these new CHAMP derived solutions with EGM96 and the EIGEN series of solutions. We review carefully the performance of these models for altimetric satellites.

Lemoime, Frank G.

An Assessment of Gravity Recovery with CHAMP Data

The CHAMP mission, launched in July 2000, is the first in the series of mapping missions for the Earth's geopotential scheduled for the first decade of the new millenium. Its unique contributions compared to all the previous generation of satellites whose data have been included in Earth geopotential models are the precision global tracking with GPS data, and the availability of precision accelerometry data to model the nonconservative forces. Over the past year we have implemented extensive modifications to our GEODYN orbit determination processing code and ancillary data preprocessors to process the GPS and accelerometry data from missions such as CHAMP and GRACE. We report on the analysis of up to 60 days of CHAMP data and how these data contribute to Earth geopotential solutions where the base model is a derivative of EGM96. Preliminary results with only 12.5 days of data processed clearly show the ability of the CHAMP data to improve the modeling of the zonals (1=10 to 40), the m-dailies, the primary resonance terms, and the sectoral harmonics. We will detail the results of our calibrations of the CHAMP accelerometry and assess the quality of test solutions that include these CHAMP data.

Lemoine, F. G.

Gravity Fields from CHAMP Mission Data

The CHAMP mission, launched in July 2000, is the first in a series of missions that will revolutionize our ability to model the Earth s geopotential. The CHAMP spacecraft is equipped for precision tracking by the Global Positioning System (GPS) and Satellite Laser Ranging (SLR) along with a precision accelerometer to provide measurements of the surface forces. Preliminary satellite-only geopotential solutions with only 30 days of CHAMP data are, by some criteria, as strong as solutions made from tracking data collected over the previous 30 years of the space age. Compared to EGM96, CHAMP makes notable contributions in regions where the terrestrial data (surface gravimetry and altimetry) were weak, for example in the polar regions, in the Amazon and the Himalayas. The CHAMP data allow us to separate the geoid from the dynamic ocean topography (DOT) up to at least degree 25 rather than just under degree 20 as in EGM96. We report on satellite-only and combination models that incorporate up to 100 days of CHAMP data as well as other satellite data. We report on our updated processing of the CHAMP tracking and accelerometer data and evaluate the performance of the geopotential models using a variety of tests.

Lemoine, Frank G.

Time-Variable Gravity Signal Due to Extratropic Pacific Water Mass Redistribution

Cox and Chao [2002] reported the detection of a large anomaly in the form of a positive "jump" in the time series of Earth's lowest-degree gravity harmonic J2, or the dynamic oblateness, during 1998. This prompted us to examine the mass redistribution in the global oceans. We report here a seesaw of the sea-surface height (SSH) in the extratropic north + south Pacific basins -- the leading (nonseasonal) EOF/PC mode in SSH derived from the 10-year TOPEX/Poseidon altimetry data in the extratropic Pacific region. The mode underwent a step-like jump with time evolution that match remarkably well with the observed J2 anomaly. However, the magnitude is several times too small to explain the observed J2, even if assuming the SSH jump was all mass-induced (as opposed to any steric effect which causes no time-variable gravity signal). If one accepts the notion that this extratropic Pacific seesaw is part of the geophysical process that produced the observed 1998 J2 anomaly, then this finding suggests strong geophysical connection of the interannual-to-decadal variation of J2 with the Pacific Decadal Oscillation (PDO), as the time series of the above EOF/PC mode is actually a formally defined PDO Index series.

Chao, Benjamin F.

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.

The Development of the Joint NASA GSFC and the National Imagery and Mapping Agency (NIMA) Geopotential Model EGM96

The NASA Goddard Space Flight Center (GSFC), the National Imagery and Mapping Agency (NIMA), and The Ohio State University (OSU) have collaborated to develop an improved spherical harmonic model of the Earth's gravitational potential to degree 360. The new model, Earth Gravitational Model 1996 (EGM96), incorporates improved surface gravity data, altimeter-derived gravity anomalies from ERS-1 and from the GEOSAT Geodetic Mission (GM), extensive satellite tracking data-including new data from Satellite Laser Ranging (SLR), the Global Postioning System (GPS), NASA's Tracking and Data Relay Satellite System (TDRSS), the French DORIS system, and the US Navy TRANET Doppler tracking system-as well as direct altimeter ranges from TOPEX/POSEIDON (T/P), ERS-1, and GEOSAT. The final solution blends a low-degree combination model to degree 70, a block-diagonal solution from degree 71 to 359, and a quadrature solution at degree 360. The model was used to compute geoid undulations accurate to better than one meter (with the exception of areas void of dense and accurate surface gravity data) and realize WGS84 as a true three-dimensional reference system. Additional results from the EGM96 solution include models of the dynamic ocean topography to degree 20 from T/P and ERS-1 together, and GEOSAT separately, and improved orbit determination for Earth-orbiting satellites.

Lemoine, F. G.

Modelling and Estimation of Dynamic Ocean Topography Within Global Geopotential Solutions

In this poster, we will present a report on the ongoing investigation "Improved Gravitational and Dynamic Height Models Through the use of Oceanographic Data." we have used a two year mean (1993-1994) of the Dynamic ocean Topography (DOT) field implied by the Semtner and Chervin POCM_4B model, and developed normal equations in surface spherical harmonics to degree and order 30. These normal equations, were combined with normal equations derived from TOPEX and ERS-1 altimeter data over the same time period. Combination solutions (based on satellite tracking data, altimeter data, surface gravity data and OCM data) were were developed Test solutions were obtained estimating the DOT field to 2Ox2O and 3ox3O. These solutions were tested with independent DOT values computed over 38 WOCE hydrographic sections, which contained a total of 3072 stations and represented 216000 km of travelled lines. The weighted standard deviation of the differences between the DOT obtained from the hydrographic data and and the field estimated from the joint combination solutions was computed for each of the test models, the weighted standard deviation for the baseline combination solution excluding the POCM4_B data was 9.7 cm for a DOT solution to 30x30, where introducing the POCM4_B data into the combination model reduced the standard deviation to 9.2 cm, indicating the the introduction of oceanographic information benefits the solution for the dynamic ocean topography. We will discuss the weighting schemes applied and the method of solution. Another aspect of our investigation involves testing alternate parametric representations of the dynamic height field. We looked at the alternative representations in terms of the Proudman functions (PF), and compared these to the use of spherical harmonics (SH) to represent the dynamic ocean topography, using once again the 1993 and 1994 output of the POCM4_B OCM as the reference model. A significant advantage of PF's over SH's is that the former require no "fill-in" values over areas where the DOT is undefined. We show that the PF and SH results using equal number of parameters agree quite well in overall content. In addition, PFs appear better suited for representation of high frequency signals close to the signal boundary. The development of the PF solutions for the DOT will be described and the statistics of the comparisons will be presented.

Lemoine, Frank G.

Production and quality assurance automation in the Goddard Space Flight Center Flight Dynamics Facility

The Flight Dynamics Facility (FDF) at the NASA Goddard Space Flight Center (GSFC) generates numerous products for NASA-supported spacecraft, including the Tracking and Data Relay Satellites (TDRS's), the Hubble Space Telescope (HST), the Extreme Ultraviolet Explorer (EUVE), and the space shuttle. These products include orbit determination data, acquisition data, event scheduling data, and attitude data. In most cases, product generation involves repetitive execution of many programs. The increasing number of missions supported by the FDF has necessitated the use of automated systems to schedule, execute, and quality assure these products. This automation allows the delivery of accurate products in a timely and cost-efficient manner. To be effective, these systems must automate as many repetitive operations as possible and must be flexible enough to meet changing support requirements. The FDF Orbit Determination Task (ODT) has implemented several systems that automate product generation and quality assurance (QA). These systems include the Orbit Production Automation System (OPAS), the New Enhanced Operations Log (NEOLOG), and the Quality Assurance Automation Software (QA Tool). Implementation of these systems has resulted in a significant reduction in required manpower, elimination of shift work and most weekend support, and improved support quality, while incurring minimal development cost. This paper will present an overview of the concepts used and experiences gained from the implementation of these automation systems.

Chapman, K. B.

Improved solution accuracy for TDRSS-based TOPEX/Poseidon orbit determination

Orbit determination results are obtained by the Goddard Space Flight Center (GSFC) Flight Dynamics Division (FDD) using a batch-least-squares estimator available in the Goddard Trajectory Determination System (GTDS) and an extended Kalman filter estimation system to process Tracking and Data Relay Satellite (TDRS) System (TDRSS) measurements. GTDS is the operational orbit determination system used by the FDD in support of the Ocean Topography Experiment (TOPEX)/Poseidon spacecraft navigation and health and safety operations. The extended Kalman filter was implemented in an orbit determination analysis prototype system, closely related to the Real-Time Orbit Determination System/Enhanced (RTOD/E) system. In addition, the Precision Orbit Determination (POD) team within the GSFC Space Geodesy Branch generated an independent set of high-accuracy trajectories to support the TOPEX/Poseidon scientific data. These latter solutions use the geodynamics (GEODYN) orbit determination system with laser ranging and Doppler Orbitography and Radiopositioning integrated by satellite (DORIS) tracking measurements. The TOPEX/Poseidon trajectories were estimated for November 7 through November 11, 1992, the timeframe under study. Independent assessments were made of the consistencies of solutions produced by the batch and sequential methods. The batch-least-squares solutions were assessed based on the solution residuals, while the sequential solutions were assessed based on primarily the estimated covariances. The batch-least-squares and sequential orbit solutions were compared with the definitive POD orbit solutions. The solution differences were generally less than 2 meters for the batch-least-squares and less than 13 meters for the sequential estimation solutions. After the sequential estimation solutions were processed with a smoother algorithm, position differences with POD orbit solutions of less than 7 meters were obtained. The differences among the POD, GTDS, and filter/smoother solutions can be traced to differences in modeling and tracking data types, which are being analyzed in detail.

Doll, C. E.

Evaluation of semiempirical atmospheric density models for orbit determination applications

This paper presents the results of an investigation of the orbit determination performance of the Jacchia-Roberts (JR), mass spectrometer incoherent scatter 1986 (MSIS-86), and drag temperature model (DTM) atmospheric density models. Evaluation of the models was performed to assess the modeling of the total atmospheric density. This study was made generic by using six spacecraft and selecting time periods of study representative of all portions of the 11-year cycle. Performance of the models was measured for multiple spacecraft, representing a selection of orbit geometries from near-equatorial to polar inclinations and altitudes from 400 kilometers to 900 kilometers. The orbit geometries represent typical low earth-orbiting spacecraft supported by the Goddard Space Flight Center (GSFC) Flight Dynamics Division (FDD). The best available modeling and orbit determination techniques using the Goddard Trajectory Determination System (GTDS) were employed to minimize the effects of modeling errors. The latest geopotential model available during the analysis, the Goddard earth model-T3 (GEM-T3), was employed to minimize geopotential model error effects on the drag estimation. Improved-accuracy techniques identified for TOPEX/Poseidon orbit determination analysis were used to improve the Tracking and Data Relay Satellite System (TDRSS)-based orbit determination used for most of the spacecraft chosen for this analysis. This paper shows that during periods of relatively quiet solar flux and geomagnetic activity near the solar minimum, the choice of atmospheric density model used for orbit determination is relatively inconsequential. During typical solar flux conditions near the solar maximum, the differences between the JR, DTM, and MSIS-86 models begin to become apparent. Time periods of extreme solar activity, those in which the daily and 81-day mean solar flux are high and change rapidly, result in significant differences between the models. During periods of high geomagnetic activity, the standard JR model was outperformed by DTM. Modification of the JR model to use a geomagnetic heating delay of 3 hours, as used in DTM, instead of the 6.7-hour delay produced results comparable to or better than the DTM performance, reducing definitive orbit solution ephermeris overlap differences by 30 to 50 percent. The reduction in the overlap differences would be useful for mitigating the impact of geomagnetic storms on orbit prediction.

Cox, C. M.

Analysis of the Cosmic Background Explorer (COBE) helium venting through effects on the orbit

An analysis is described which shows that a helium venting thrust level that is 60 percent higher than the predicted 2.3 millinewtons may be the cause of the unusually high premission-analysis and on-orbit support orbital decay rates and drag correction factors found for the Cosmic Background Explorer satellite. Depletion of the helium allowed a determination of the average mass flow rate and specific impulse. The average mass flow rate of 3.3 mg/sec confirmed the thrust levels indicated by the venting thrust analysis. Liquid helium is used for cooling instruments aboard the satellite.

Cox, C. M.

Helicopter noise

Strategies for the reduction of noise at and around heliports are discussed. Heliport planning and operations and piloting techniques are addressed.

Cox, C. M.