New Analysis of Mercury Laser Altimeter Crossovers to Improve Geodetic Constraints by MESSENGER
No abstract available
Bertone, S.↗
Engineering topics
Sabaka, T. J.
Publications and source records attributed to Sabaka, T. J..
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.↗
Antarctica, Greenland and Gulf of Alaska Land-ice Evolution from an Iterated GRACE Global Mascon Solution
We have determined the ice mass evolution of the Antarctica and Greenland ice sheets (AIS and GIS) and Gulf of Alaska (GOA) glaciers from a new GRACE global solution of equal-area surface mass concentration parcels (mascons) in equivalent height of water. The mascons were estimated directly from the reduction of the inter-satellite K-band range-rate (KBRR) 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 were estimated with 10 day and 1 arc degree equal-area sampling, applying anisotropic constraints. An ensemble empirical mode decomposition adaptive filter was applied to the mascon time series to compute annual mass balances. The details and causes of the spatial and temporal variability of the land-ice regions studied are discussed. The estimated mass trend over the total GIS, AIS and GOA glaciers for the time period 1 December 2003 to 1 December 2010 is -380 plus or minus 31 Gt a(exp -1), equivalent to -1.05 plus or minus 0.09 mma(exp -1) sea-level rise. Over the same time period we estimate the mass acceleration to be -41 plus or minus 27 Gt a(exp -2), equivalent to a 0.11 plus or minus 0.08 mm a(exp -2) rate of change in sea level. The trends and accelerations are dependent on significant seasonal and annual balance anomalies.
MASCON↗
High Degree and Order Gravity Fields of the Moon Derived from GRAIL Data
The Gravity Recovery and Interior Laboratory (GRAIL) spacecraft conducted the mapping of the gravity field of the Moon from March 1, 2012 to May 29, 2012. The twin spacecraft acquired highly precise K Band range-rate (KBRR) intersatellite ranging data and Deep Space Network (DSN) data during this prime mission phase from altitudes of 15 to 75 km above the lunar surface over three lunar months. We have processed these data using the NASA GSFC GEODYN orbit determination and geodetic parameter estimation program, and we have determined gravity fields up to degree and order 420 in spherical harmonics. The new gravity solutions show improved correlations with LOLA-derived topography to high degree and order and resolve many lunar features in the geopotential with a resolution of less than 30 km, including for example the central peak of the crater Tycho. We discuss the methodology used for the processing of the GRAIL data, the quality of the orbit determination on the GRAIL satellites and the derivation of the solutions, and their evaluation with independent data, including Lunar Prospector. We show that with these new GRAIL gravity solutions, we can now fit the low altitude, extended mission Lunar Prospector tracking data better than with any previous gravity model that included the LP data.
Lemoine, F. G.↗
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.↗
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.↗
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.↗
Magnetic Spectra of Earth and Mars
The spectral method for distinguishing crustal from core-source magnetic fields has been re-examined, modified, and applied to both a comprehensive geomagnetic field model and an altitude normalized magnetic map of Mars. Additional information is contained in the original extended abstract.
Voorhies, C. V.↗
On Magnetic Spectra of Earth and Mars
The spectral method for distinguishing crustal from core-source magnetic fields is reexamined, modified, and applied to both a comprehensive geomagnetic field model and an altitude normalized magnetic map of Mars. The observational spectra are fairly fitted by theoretical forms expected from certain elementary classes of magnetic sources. For Earth we find fields from a core of radius 3512 +/- 64 km, in accord with the seismologic core radius of 3480 km, and a crust represented by a shell of random dipolar sources at radius 6367 +/- 14 km, near the planetary mean radius of 6371.2 km. For Mars we find no sign of a core-source field, only a field from a crust represented in same way, but at radius 3344 +/- 10 km, about 46 km below the planetary mean radius of 3389.5 km, and with sources about 9.6 +/- 3.2 times stronger.
Voorhies, C. V.↗
An initial analysis of the data from the Polar Orbiting Geophysical (POGS) Satellite
The Polar Orbiting Geophysical Satellite (POGS) was launched in 1990 to measure the geomagnetic field. POGS data from selected magnetically quiet days was chosen, quality checked and deleted where thought to be erroneous. A time and position correction was applied. The resulting data was fit to a degree 13 spherical harmonic model. Evaluation of the quality of the data indicates that it is sufficient for definition of the low degree (approximately less than 8) portion of the geomagnetic field. Further correction of the data time and position may improve this quality.
Langel, R. A.↗
Geomagnetic model investigations for 1980 - 1989: A model for strategic defense initiative particle beam experiments and a study in the effects of data types and observatory bias solutions
Two suites of geomagnetic field models were generated at the request of Los Alamos National Lab. concerning Strategic Defense Initiative (SDI) research. The first is a progression of five models incorporating MAGSAT data and data from a sequence of batches as a priori information. The batch sequence is: post 1979.5 observatory data, post 1980 land survey and selected aeromagnetic and marine survey data, a special White Sands (NM) area survey by Project Magnet with some additional post 1980 marine survey data, and finally DE-2 satellite data. These models are of 13th deg and order in their main field terms, and deg and order 10 in their first derivative temporal terms. The second suite consists of four models based solely upon post 1983.5 observatory and survey data. They are of deg and order 10 in main field and 8 in a first deg Taylor series. A comprehensive error analysis was applied to both series, which accounted for error sources such as the truncated core and crustal fields, and the neglected Sq and low deg crustal fields. Comparison of the power spectrum of the MGST (10/81) model with those of this series show good agreement.
Langel, Robert A.↗
Processing of DMSP magnetic data: Handbook of programs, tapes, and datasets
The DMSP F-7 satellite was an operational Air Force meteorological satellite which carried a magnetometer for geophysical measurements. The magnetometer was located within the body of the spacecraft in the presence of large spacecraft fields. In addition to stray magnetic fields, the data have inherent position and time inaccuracies. Algorithms were developed to identify and remove time varying magnetic field noise from the data. These algorithms are embodied in an automated procedure which fits a smooth curve through the data and then identifies outliers and which filters the predominant Fourier component of noise from the data. Techniques developed for Magsat were then modified and used to attempt determination of the spacecraft fields, of any rotation between the magnetometer axes and the spacecraft axes, and of any scale changes within the magnetometer itself. Software setup and usage are documented and program listings are included in the Appendix. The initial and resulting data are archived on magnetic cartridge and the formats are documented.
Langel, R. A.↗
Processing of DMSP magnetic data and its use in geomagnetic field modeling
The DMSP F-7 satellite is an operational Air Force meteorological satellite which carried a magnetometer for geophysical measurements. The magnetometer was located within the body of the spacecraft in the presence of large spacecraft fields. In addition to stray magnetic fields, the data have inherent position and time inaccuracies. Algorithms were developed to identify and remove time varying magnetic field noise from the data. Techniques developed for Magsat were then modified and used to attempt determination of the spacecraft fields, of any rotation between the magnetometer axes and the spacecraft axes, and of any scale changes within the magnetometer itself. The corrected data were then used to attempt to model the geomagnetic field. This was done in combination with data from Magsat, from the standard magnetic observatories, from aeromagnetic and other survey data, and from DE-2 spacecraft field data. Future DMSP missions can be upgraded in terms of geomagnetic measurements by upgrading the time and position information furnished with the data, placing the magnetometer at the end of the boom, upgrading the attitude determination at the magnetometer, and increasing the accuracy of the magnetometer.
Ridgway, J. R.↗
Uncertainty estimates in geomagnetic field modeling
This paper presents an extension of the conventional uncertainty analysis which characterizes the sources of uncertainty in the coefficients of the geomagnetic field models. The new formalism accounts for the systematic errors introduced by the omission of such sources as the presence of crustal fields, the external fields, and the field from the truncated terms. The usefulness of this formalism depends on two critical conditions. The first of these is the knowledge of the statistical properties of the fields whose parameters are not solved in the analysis, i.e., the crustal field and the external field. The second critical point in the practical use of the method is the approximation used for the weight matrix.
Langel, R. A.↗
Definitive IGRF models for 1945, 1950, 1955, and 1960
Definitive internal geomagnetic field models are obtained for the epochs 1945, 1950, 1955, and 1960, with the solution including local biases for the fixed observatories. In order to derive a priori models for the epochs, the International Geomagnetic Reference Field (12/83) model (based on Magsat data) was projected backward in time using a spline fit to annual secular-variation models. An uncertainty model was developed which included estimates of the effects of crustal and core fields not represented by the model. Differences between biases from year to year are found to be within the predicted bounds.
Langel, R. A.↗