Showing papers by "Frank G. Lemoine published in 2014"

Accuracy assessment of global barotropic ocean tide models

Abstract: The accuracy of state-of-the-art global barotropic tide models is assessed using bottom pressure data, coastal tide gauges, satellite altimetry, various geodetic data on Antarctic ice shelves, and independent tracked satellite orbit perturbations. Tide models under review include empirical, purely hydrodynamic (“forward”), and assimilative dynamical, i.e., constrained by observations. Ten dominant tidal constituents in the diurnal, semidiurnal, and quarter-diurnal bands are considered. Since the last major model comparison project in 1997, models have improved markedly, especially in shallow-water regions and also in the deep ocean. The root-sum-square differences between tide observations and the best models for eight major constituents are approximately 0.9, 5.0, and 6.5 cm for pelagic, shelf, and coastal conditions, respectively. Large intermodel discrepancies occur in high latitudes, but testing in those regions is impeded by the paucity of high-quality in situ tide records. Long-wavelength components of models tested by analyzing satellite laser ranging measurements suggest that several models are comparably accurate for use in precise orbit determination, but analyses of GRACE intersatellite ranging data show that all models are still imperfect on basin and subbasin scales, especially near Antarctica. For the M2 constituent, errors in purely hydrodynamic models are now almost comparable to the 1980-era Schwiderski empirical solution, indicating marked advancement in dynamical modeling. Assessing model accuracy using tidal currents remains problematic owing to uncertainties in in situ current meter estimates and the inability to isolate the barotropic mode. Velocity tests against both acoustic tomography and current meters do confirm that assimilative models perform better than purely hydrodynamic models.

Lunar interior properties from the GRAIL mission

Abstract: The Gravity Recovery and Interior Laboratory (GRAIL) mission has sampled lunar gravity with unprecedented accuracy and resolution. The lunar GM, the product of the gravitational constant G and the mass M, is very well determined. However, uncertainties in the mass and mean density, 3345.56 ± 0.40 kg/m3, are limited by the accuracy of G. Values of the spherical harmonic degree-2 gravity coefficients J2 and C22, as well as the Love number k2 describing lunar degree-2 elastic response to tidal forces, come from two independent analyses of the 3 month GRAIL Primary Mission data at the Jet Propulsion Laboratory and the Goddard Space Flight Center. The two k2 determinations, with uncertainties of ~1%, differ by 1%; the average value is 0.02416 ± 0.00022 at a 1 month period with reference radius R = 1738 km. Lunar laser ranging (LLR) data analysis determines (C − A)/B and (B − A)/C, where A < B < C are the principal moments of inertia; the flattening of the fluid outer core; the dissipation at its solid boundaries; and the monthly tidal dissipation Q = 37.5 ± 4. The moment of inertia computation combines the GRAIL-determined J2 and C22 with LLR-derived (C − A)/B and (B − A)/C. The normalized mean moment of inertia of the solid Moon is Is/MR2 = 0.392728 ± 0.000012. Matching the density, moment, and Love number, calculated models have a fluid outer core with radius of 200–380 km, a solid inner core with radius of 0–280 km and mass fraction of 0–1%, and a deep mantle zone of low seismic shear velocity. The mass fraction of the combined inner and outer core is ≤1.5%.

GRGM900C: A degree 900 lunar gravity model from GRAIL primary and extended mission data

Abstract: We have derived a gravity field solution in spherical harmonics to degree and order 900, GRGM900C, from the tracking data of the Gravity Recovery and Interior Laboratory (GRAIL) Primary (1 March to 29 May 2012) and Extended Missions (30 August to 14 December 2012). A power law constraint of 3.6 ×10−4/l2 was applied only for degree l greater than 600. The model produces global correlations of gravity, and gravity predicted from lunar topography of ≥ 0.98 through degree 638. The model's degree strength varies from a minimum of 575–675 over the central nearside and farside to 900 over the polar regions. The model fits the Extended Mission Ka-Band Range Rate data through 17 November 2012 at 0.13 μm/s RMS, whereas the last month of Ka-Band Range-Rate data obtained from altitudes of 2–10 km fit at 0.98 μm/s RMS, indicating that there is still signal inherent in the tracking data beyond degree 900.

The gravity field, orientation, and ephemeris of Mercury from MESSENGER observations after three years in orbit

Abstract: We have analyzed 3 years of radio tracking data from the MESSENGER spacecraft in orbit around Mercury and determined the gravity field, planetary orientation, and ephemeris of the innermost planet. With improvements in spatial coverage, force modeling, and data weighting, we refined an earlier global gravity field both in quality and resolution, and we present here a spherical harmonic solution to degree and order 50. In this field, termed HgM005, uncertainties in low-degree coefficients are reduced by an order of magnitude relative to earlier global fields, and we obtained a preliminary value of the tidal Love number k2 of 0.451 ± 0.014. We also estimated Mercury's pole position, and we obtained an obliquity value of 2.06 ± 0.16 arcmin, in good agreement with analysis of Earth-based radar observations. From our updated rotation period (58.646146 ± 0.000011 days) and Mercury ephemeris, we verified experimentally the planet's 3:2 spin-orbit resonance to greater accuracy than previously possible. We present a detailed analysis of the HgM005 covariance matrix, and we describe some near-circular frozen orbits around Mercury that could be advantageous for future exploration.

External Evaluation of the Terrestrial Reference Frame: Report of the Task Force of the IAG Sub-commission 1.2

Abstract: Ideally, the origin of the Terrestrial Reference Frame (TRF) is defined as the center of mass of the whole Earth system, the time evolution of its orientation is such that no global net rotation of the whole Earth’s surface is possible and the TRF scale is specified through the adoption of some physical constants and time-scale. These parameters need to be accurately determined since their choice has an influence on many Earth’s science applications. The aim of the task force “External evaluation of the Terrestrial Reference Frame” is to review all the applications for which the TRF accuracy is of fundamental importance. As the TRF choice has an influence on the interpretation of the results in these specific applications, we investigate if some evaluation procedures could be established. We classified the methods that allow evaluation of the TRF using ground, geodetic data or models that have not been used in the TRF construction, based on their expected contributions. Some of these methods have been applied to the latest International Terrestrial Reference System realizations and the results are presented here. Although further analysis will be necessary to deliver a more precise error budget, our findings demonstrate that the most recent realizations of the ITRS are more accurate than the previous in terms of origin and scale rate definition. The current level of ITRF2008 accuracy is likely to be at the level of 0.5 mm/year along each origin component and better than 0.3 mm/year in the scale rate according to the most recent studies.

The International DORIS Service (IDS) - Recent Developments in Preparation for ITRF2013

Abstract: The International DORIS Service (IDS) was created in 2003 under the umbrella of the International Association of Geodesy (IAG) to foster scientific research related to the French DORIS tracking system and to deliver scientific products, mostly related to the International Earth rotation and Reference systems Service (IERS). We first present some general background related to the DORIS system (current and planned satellites, current tracking network and expected evolution) and to the general IDS organization (from Data Centers, Analysis Centers and Combination Center). Then, we discuss some of the steps recently taken to prepare the IDS submission to ITRF2013 (combined weekly time series based on individual solutions from several Analysis Centers). In particular, recent results obtained from the Analysis Centers and the Combination Center show that improvements can still be made when updating physical models of some DORIS satellites, such as Envisat, Cryosat-2 or Jason-2. The DORIS contribution to ITRF2013 should also benefit from the larger number of ground observations collected by the last generation of DGXX receivers (first instrument being onboard Jason-2 satellite). In particular for polar motion, sub-milliarcsecond accuracy seems now to be achievable. Weekly station positioning internal consistency also seems to be improved with a larger DORIS constellation.

Estimated SLR station position and network frame sensitivity to time-varying gravity

Abstract: This paper evaluates the sensitivity of ITRF2008-based satellite laser ranging (SLR) station positions estimated weekly using LAGEOS-1/2 data from 1993 to 2012 to non-tidal time-varying gravity (TVG). Two primary methods for modeling TVG from degree-2 are employed. The operational approach applies an annual GRACE-derived field, and IERS recommended linear rates for five coefficients. The experimental approach uses low-order/degree $$4\times 4$$ coefficients estimated weekly from SLR and DORIS processing of up to 11 satellites (tvg4x4). This study shows that the LAGEOS-1/2 orbits and the weekly station solutions are sensitive to more detailed modeling of TVG than prescribed in the current IERS standards. Over 1993–2012 tvg4x4 improves SLR residuals by 18 % and shows 10 % RMS improvement in station stability. Tests suggest that the improved stability of the tvg4x4 POD solution frame may help clarify geophysical signals present in the estimated station position time series. The signals include linear and seasonal station motion, and motion of the TRF origin, particularly in Z. The effect on both POD and the station solutions becomes increasingly evident starting in 2006. Over 2008–2012, the tvg4x4 series improves SLR residuals by 29 %. Use of the GRGS RL02 $$50\times 50$$ series shows similar improvement in POD. Using tvg4x4, secular changes in the TRF origin Z component double over the last decade and although not conclusive, it is consistent with increased geocenter rate expected due to continental ice melt. The test results indicate that accurate modeling of TVG is necessary for improvement of station position estimation using SLR data.

High-resolution local gravity model of the south pole of the Moon from GRAIL extended mission data

Abstract: We estimated a high-resolution local gravity field model over the south pole of the Moon using data from the Gravity Recovery and Interior Laboratory's extended mission. Our solution consists of adjustments with respect to a global model expressed in spherical harmonics. The adjustments are expressed as gridded gravity anomalies with a resolution of 1/6° by 1/6° (equivalent to that of a degree and order 1080 model in spherical harmonics), covering a cap over the south pole with a radius of 40°. The gravity anomalies have been estimated from a short-arc analysis using only Ka-band range-rate (KBRR) data over the area of interest. We apply a neighbor-smoothing constraint to our solution. Our local model removes striping present in the global model; it reduces the misfit to the KBRR data and improves correlations with topography to higher degrees than current global models. Key Points We present a high-resolution gravity model of the south pole of the Moon Improved correlations with topography to higher degrees than global models Improved fits to the data and reduced striping that is present in global models

The effect of seasonal and long-period geopotential variations on the GPS orbits

Abstract: We examine the impact of using seasonal and long-period time-variable gravity field (TVG) models on GPS orbit determination, through simulations from 1994 to 2012. The models of time-variable gravity that we test include the GRGS release RL02 GRACE-derived 10-day gravity field models up to degree and order 20 (grgs20x20), a 4 A— 4 series of weekly coefficients using GGM03S as a base derived from SLR and DORIS tracking to 11 satellites (tvg4x4), and a harmonic fit to the above 4 A— 4 SLR---DORIS time series (goco2s_fit2). These detailed models are compared to GPS orbit simulations using a reference model (stdtvg) based on the International Earth Rotation Service (IERS) and International GNSS Service (IGS) repro1 standards. We find that the new TVG modeling produces significant along, cross-track orbit differences as well as annual, semi-annual, draconitic and long-period effects in the Helmert translation parameters (Tx, Ty, Tz) of the GPS orbits with magnitudes of several mm. We show that the simplistic TVG modeling approach used by all of the IGS Analysis Centers, which is based on the models provided by the IERS standards, becomes progressively less adequate following 2006 when compared to the seasonal and long-period TVG models.