Engineering topics

Luthcke, S. B.

Publications and source records attributed to Luthcke, S. B..

At least 19 records

NASA NTRS2019

Improved Earth Oblateness Rate Reveals Increased Ice Sheet Losses and Mass-Driven Sea Level Rise

Satellite laser ranging (SLR) observations are routinely applied toward the estimation of dynamic oblateness, C(sub 20), which is the largest globally integrated component of Earth's time-variable gravity field. Since 2002, GRACE and GRACE Follow-On have revolutionized the recovery of higher spatial resolution features of global time-variable gravity, with SLR continuing to provide the most reliable estimates of C (sub 20).We quantify the effect of various SLR processing strategies on estimating C(sub 20) and demonstrate better signal recovery with the inclusion of GRACE-derived low-degree gravity information in the forward model. This improved SLR product modifies the Antarctic and Greenland Ice Sheet mass trends by -15.4 and -3.5 Gt/year, respectively, as compared to CSR TN11, and improves global mean sea level budget closure by modifying sea level rise by +0.08 mm/year. We recommend that this new C(sub 20) product be applied to RL06 GRACE data products for enhanced accuracy and scientific interpretation.

Loomis, B. D.↗

NASA NTRS2019

Regularization and Error Characterization of GRACE Mascons

We present a new global time-variable gravity mascon solution derived from Gravity Recovery and Climate Experiment (GRACE) Level 1B data. The new product from the NASA Goddard Space Flight Center (GSFC) results from a novel approach that combines an iterative solution strategy with geographical binning of inter-satellite range-acceleration residuals in the construction of time-dependent regularization matrices applied in the inversion of mascon parameters. This estimation strategy is intentionally conservative as it seeks to maximize the role of the GRACE measurements on the final solution while minimizing the influence of the regularization design process. We fully reprocess the Level 1B data in the presence of the final mascon solution to generate true post-fit inter-satellite residuals, which are utilized to confirm solution convergence and to validate the mascon noise uncertainties. We also present the mathematical case that regularized mascon solutions are biased, and that this bias, or leakage, must be combined with the estimated noise variance to accurately assess total mascon uncertainties. The estimated leakage errors are determined from the monthly resolution operators. We present a simple approach to compute the total uncertainty for both individual mascon and regional analysis of the GSFC mascon product, and validate the results in comparison with independent mascon solutions and calibrated Stokes uncertainties. Lastly, we present the new solution and uncertainties with global analyses of the mass trends and annual amplitudes, and compute updated trends for the global ocean, and the respective contributions of the Greenland Ice Sheet, Antarctic Ice Sheet, Gulf of Alaska, and terrestrial water storage. This analysis highlights the successful closure of the global mean sea level budget, that is, the sum of global ocean mass from the GSFC mascons and the steric component from Argo floats agrees well with the total determined from sea surface altimetry.

Loomis, B. D.↗

NASA NTRS2016

Mass Evolution of Mediterranean, Black, Red, and Caspian Seas from GRACE and Altimetry: Accuracy Assessment and Solution Calibration

We present new measurements of mass evolution for the Mediterranean, Black, Red, and Caspian Seas as determined by the NASA Goddard Space Flight Center (GSFC) GRACE time-variable global gravity mascon solutions. These new solutions are compared to sea surface altimetry measurements of sea level anomalies with steric corrections applied. To assess their accuracy, the GRACE and altimetry-derived solutions are applied to the set of forward models used by GSFC for processing the GRACE Level-1B datasets, with the resulting inter-satellite range acceleration residuals providing a useful metric for analyzing solution quality.

Loomis, B. D.↗

NASA NTRS2012

An Iterated Global Mascon Solution with Focus on Land Ice Mass Evolution

Land ice mass evolution is determined from a new GRACE global mascon solution. The solution is estimated directly from the reduction of the inter-satellite K-band range rate observations taking into account the full noise covariance, and formally iterating the solution. The new solution increases signal recovery while reducing the GRACE KBRR observation residuals. The mascons are estimated with 10-day and 1-arc-degree equal area sampling, applying anisotropic constraints for enhanced temporal and spatial resolution of the recovered land ice signal. The details of the solution are presented including error and resolution analysis. An Ensemble Empirical Mode Decomposition (EEMD) adaptive filter is applied to the mascon solution time series to compute timing of balance seasons and annual mass balances. The details and causes of the spatial and temporal variability of the land ice regions studied are discussed.

Luthcke, S. B.↗

NASA NTRS2012

The Effect of Geocenter Motion on Jason-2 Orbits and the Mean Sea Level

We compute a series of Jason-2 GPS and SLR/DORIS-based orbits using ITRF2005 and the std0905 standards (Lemoine et al. 2010). Our GPS and SLR/DORIS orbit data sets span a period of 2 years from cycle 3 (July 2008) to cycle 74 (July 2010). We extract the Jason-2 orbit frame translational parameters per cycle by the means of a Helmert transformation between a set of reference orbits and a set of test orbits. We compare the annual terms of these time-series to the annual terms of two different geocenter motion models where biases and trends have been removed. Subsequently, we include the annual terms of the modeled geocenter motion as a degree-1 loading displacement correction to the GPS and SLR/DORIS tracking network of the POD process. Although the annual geocenter motion correction would reflect a stationary signal in time, under ideal conditions, the whole geocenter motion is a non-stationary process that includes secular trends. Our results suggest that our GSFC Jason-2 GPS-based orbits are closely tied to the center of mass (CM) of the Earth consistent with our current force modeling, whereas GSFC's SLR/DORIS-based orbits are tied to the origin of ITRF2005, which is the center of figure (CF) for sub-secular scales. We quantify the GPS and SLR/DORIS orbit centering and how this impacts the orbit radial error over the globe, which is assimilated into mean sea level (MSL) error, from the omission of the annual term of the geocenter correction. We find that for the SLR/DORIS std0905 orbits, currently used by the oceanographic community, only the negligence of the annual term of the geocenter motion correction results in a 4.67 plus or minus 3.40 mm error in the Z-component of the orbit frame which creates 1.06 plus or minus 2.66 mm of systematic error in the MSL estimates, mainly due to the uneven distribution of the oceans between the North and South hemisphere.

Melachroinos, S. A.↗

NASA NTRS2011

A New Global Mascon Solution Tuned for High-Latitude Ice Studies

A new global mascon solution has been developed with I-arc-degree spatial and IO-day temporal sampling. The global mas cons are estimated from the reduction of nearly 8 years of GRACE K-band range-rate data. Temporal and anisotropic spatial constraints have been applied for land, ocean and ice regions. The solution construction and tuning is focused towards the Greenland and Antarctic ice sheets (GIS and AIS) as well as the Gulf of Alaska mountain glaciers (GoA). Details of the solution development will be discussed, including the mascon parameter definitions, constraints, and the tuning of the constraint damping factor. Results will be presented, exploring the spatial and temporal variability of the ice sheets and GoA regions. A detailed error analysis will be discussed, including solution dependence on iteration, damping factor, forward modeling, and multitechnique comparisons. We also investigate the fundamental resolution of the solution and the spatial correlation of ice sheet inter-annual change. Finally, we discuss future improvements, including specific constraint application for the rest of the major land ice regions and improvements in solution regularization.

Luthcke, S. B.↗

NASA NTRS2011

Estimation and Validation of Oceanic Mass Circulation from the GRACE Mission

Since the launch of the Gravity Recovery And Climate Experiment (GRACE) in March 2002, the Earth's surface mass variations have been monitored with unprecedented accuracy and resolution. Compared to the classical spherical harmonic solutions, global high-resolution mascon solutions allows the retrieval of mass variations with higher spatial and temporal sampling (2 degrees and 10 days). We present here the validation of the GRACE global mascon solutions by comparing mass estimates to a set of about 100 ocean bottom pressure (OSP) records, and show that the forward modelling of continental hydrology prior to the inversion of the K-band range rate data allows better estimates of ocean mass variations. We also validate our GRACE results to OSP variations modelled by different state-of-the-art ocean general circulation models, including ECCO (Estimating the Circulation and Climate of the Ocean) and operational and reanalysis from the MERCATOR project.

Boy, J.-P.↗

NASA NTRS2010

Improving Global Mass Flux Solutions from Gravity Recovery and Climate Experiment (GRACE) Through Forward Modeling and Continuous Time Correlation

We describe Earth's mass flux from April 2003 through November 2008 by deriving a time series of mas cons on a global 2deg x 2deg equal-area grid at 10 day intervals. We estimate the mass flux directly from K band range rate (KBRR) data provided by the Gravity Recovery and Climate Experiment (GRACE) mission. Using regularized least squares, we take into account the underlying process dynamics through continuous space and time-correlated constraints. In addition, we place the mascon approach in the context of other filtering techniques, showing its equivalence to anisotropic, nonsymmetric filtering, least squares collocation, and Kalman smoothing. We produce mascon time series from KBRR data that have and have not been corrected (forward modeled) for hydrological processes and fmd that the former produce superior results in oceanic areas by minimizing signal leakage from strong sources on land. By exploiting the structure of the spatiotemporal constraints, we are able to use a much more efficient (in storage and computation) inversion algorithm based upon the conjugate gradient method. This allows us to apply continuous rather than piecewise continuous time-correlated constraints, which we show via global maps and comparisons with ocean-bottom pressure gauges, to produce time series with reduced random variance and full systematic signal. Finally, we present a preferred global model, a hybrid whose oceanic portions are derived using forward modeling of hydrology but whose land portions are not, and thus represent a pure GRACE-derived signal.

Sabaka, T. J.↗

NASA NTRS2010

Global Mass Flux Solutions from GRACE: A Comparison of Parameter Estimation Strategies - Mass Concentrations Versus Stokes Coefficients

The differences between mass concentration (mas con) parameters and standard Stokes coefficient parameters in the recovery of gravity infonnation from gravity recovery and climate experiment (GRACE) intersatellite K-band range rate data are investigated. First, mascons are decomposed into their Stokes coefficient representations to gauge the range of solutions available using each of the two types of parameters. Next, a direct comparison is made between two time series of unconstrained gravity solutions, one based on a set of global equal area mascon parameters (equivalent to 4deg x 4deg at the equator), and the other based on standard Stokes coefficients with each time series using the same fundamental processing of the GRACE tracking data. It is shown that in unconstrained solutions, the type of gravity parameter being estimated does not qualitatively affect the estimated gravity field. It is also shown that many of the differences in mass flux derivations from GRACE gravity solutions arise from the type of smoothing being used and that the type of smoothing that can be embedded in mas con solutions has distinct advantages over postsolution smoothing. Finally, a 1 year time series based on global 2deg equal area mascons estimated every 10 days is presented.

Rowlands, D. D.↗

NASA NTRS2009

A Simulation Study of Multi-Beam Altimetry for Lunar Reconnaissance Orbiter and Other Planetary Missions

The combined use of altimetry, Earth-based Doppler and Earth-based range measurements in the lunar reconnaissance orbiter (LRO) mission (Chin et al. in Space Sci Rev 129:391-419, 2007) has been examined in a simulation study. It is found that in the initial phases of the mission orbit and altimeter geolocation accuracies should be better than 10m in the radial component and 60m overall. It is demonstrated that LRO's precise 1-way laser range measurement from Earth-based stations (Smith et al. in Proceedings of the 15th International Laser Ranging Workshop, Canberra, Australia, October 15-20, 2006) will be useful for gravity recovery. The advantages of multiple laser beams are demonstrated for altimeter calibration, orbit determination and gravity recovery in general planetary settings as well as for LRO.

Rowlands, D. D.↗

NASA NTRS2006

13 Years of TOPEX/POSEIDON Precision Orbit Determination and the 10-fold Improvement in Expected Orbit Accuracy

Launched in the summer of 1992, TOPEX/POSEIDON (T/P) was a joint mission between NASA and the Centre National d Etudes Spatiales (CNES), the French Space Agency, to make precise radar altimeter measurements of the ocean surface. After the remarkably successful 13-years of mapping the ocean surface T/P lost its ability to maneuver and was de-commissioned January 2006. T/P revolutionized the study of the Earth s oceans by vastly exceeding pre-launch estimates of surface height accuracy recoverable from radar altimeter measurements. The precision orbit lies at the heart of the altimeter measurement providing the reference frame from which the radar altimeter measurements are made. The expected quality of orbit knowledge had limited the measurement accuracy expectations of past altimeter missions, and still remains a major component in the error budget of all altimeter missions. This paper describes critical improvements made to the T/P orbit time series over the 13-years of precise orbit determination (POD) provided by the GSFC Space Geodesy Laboratory. The POD improvements from the pre-launch T/P expectation of radial orbit accuracy and Mission requirement of 13-cm to an expected accuracy of about 1.5-cm with today s latest orbits will be discussed. The latest orbits with 1.5 cm RMS radial accuracy represent a significant improvement to the 2.0-cm accuracy orbits currently available on the T/P Geophysical Data Record (GDR) altimeter product.

Lemoine, F. G.↗

NASA NTRS2005

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.↗

NASA NTRS2004

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.↗

NASA NTRS2003

An Evaluation of Recent Gravity Models wrt. Altimeter Satellite Missions

With the launch of CHAMP and GRACE, we have entered a new phase in the history of satellite geodesy. For the first time, geopotential models are now available based almost exclusively on satellite-satellite tracking either with GPS in the case of the CHAMP-based geopotential models, or co-orbital intersatellite ultra-precise ranging in the case of GRACE. Different groups have analyzed these data, and produced a series of geopotential models (e.g., EIGENlS, EIGEN2, GGM0lS, GGMOlC) that incorporate the new data. We will compare the performance of these "newer" geopotential models with the standard models now used for computations, (e.g., JGM-3, BGM-96, PGS7727, and GRIMS-C1) for TOPEX, JASON, Geosat-Follow-On (GFO), and Envisat using standard metrics such as SLR RMS of fit, altimeter crossovers, and orbit overlaps. Where covariances are available we can evaluate the predicted geographically correlated orbit error. These predicted results can be compared with the Earth-fixed differences between dynamic and reduced-dynamic orbits to test the predictive accuracy of the covariances, as well as to calibrate the error of the solutions.

Lemoine, Frank G.↗

NASA NTRS2002

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.↗

NASA NTRS2002

Recent Results from CHAMP Tracking and Accelerometer Data Analysis

The CHAMP mission's unique combination of sensors and orbit configuration will enable unprecedented improvements in modeling and understanding the Earth's static gravity field and its temporal variations. CHAMP is the first of two missions (GRACE to be launched in the early part of 02') that combine a new generation of Global Positioning System (GPS) receivers, a high precision three-axis accelerometer, and star cameras for the precision attitude determination. In order to isolate the gravity signal for science investigations, it is necessary to perform a detailed reduction and analysis of the GPS and Satellite Laser Ranging (SLR) tracking data in conjunction with the accelerometer and attitude data. Precision orbit determination based on the GPS and SLR tracking data will isolate the orbit perturbations, while the accelerometer data will be used to distinguish the non-gravitational forces from those due to the geopotential (static, and time varying). In preparation for the CHAMP and GRACE missions, extensive modifications have been made to NASA/GSFC's GEODYN orbit determination software to enable the simultaneous reduction of spacecraft tracking (e.g. GPS and SLR), three-axis accelerometer and precise attitude data. Several weeks of CHAMP tracking and accelerometer data have been analyzed and the results will be presented. Precision orbit determination analysis based on tracking data alone in addition to results based on the simultaneous reduction of tracking and accelerometer data will be discussed. Results from a calibration of the accelerometer will be presented along with the results from various orbit determination strategies.

Luthcke, S. B.↗

NASA NTRS2002

CHAMP Tracking and Accelerometer Data Analysis Results

The CHAMP (Challenging Minisatellite Payload) mission's unique combination of sensors and orbit configuration will enable unprecedented improvements in modeling and understanding the Earth's static gravity field and its temporal variations. CHAMP is the first of two missions (GRACE (Gravity Recovery and Climate Experiment) to be launched in the later part of '01) that combine a new generation of GPS (Global Positioning System) receivers, a high precision three axis accelerometer, and star cameras for the precision attitude determination. In order to isolate the gravity signal for science investigations, it is necessary to perform a detailed reduction and analysis of the GPS and SLR tracking data in conjunction with the accelerometer and attitude data. Precision orbit determination based on the GPS and SLR (Satellite Laser Ranging) tracking data will isolate the orbit perturbations, while the accelerometer data will be used to distinguish the surface forces from those due to the geopotential (static, and time varying). In preparation for the CHAMP and GRACE missions, extensive modifications have been made to NASA/GSFC's GEODYN orbit determination software to enable the simultaneous reduction of spacecraft tracking (e.g. GPS and SLR), three axis accelerometer and precise attitude data. Several weeks of CHAMP tracking and accelerometer data have been analyzed and the results will be presented. Precision orbit determination analysis based on tracking data alone in addition to results based on the simultaneous reduction of tracking and accelerometer data will be discussed. Results from a calibration of the accelerometer will be presented along with the results from various orbit determination strategies. Gravity field modeling status and plans will be discussed.

Lemoine, Frank G.↗

NASA NTRS2002

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.↗