We report on new observations of the orbital position of Phobos, and use them to obtain a new and improved estimate of the rate of secular acceleration in longitude due to tidal dissipation within Mars. Phobos is the inner-most natural satellite of Mars, and one of the few natural satellites in the solar system with orbital period shorter than the rotation period of its primary. As a result, any departure from a perfect elastic response by Mars in the tides raised on it by Phobos will cause a transfer of angular momentum from the orbit of Phobos to the spin of Mars. Since its discovery in 1877, Phobos has completed over 145,500 orbits, and has one of the best studied orbits in the solar system, with over 6000 earth-based astrometric observations, and over 300 spacecraft observations. As early as 1945, Sharpless noted that there is a secular acceleration in mean longitude, with rate (1.88 + 0.25) 10(exp -3) degrees per square year. In preparation for the 1989 Russian spacecraft mission to Phobos, considerable work was done compiling past observations, and refining the orbital model. All of the published estimates from that era are in good agreement. A typical solution (Jacobson et al., 1989) yields (1.249 + 0.018) 10(exp -3) degrees per square year. The MOLA instrument on MGS is a laser altimeter, and was designed to measure the topography of Mars. However, it has also been used to make observations of the position of Phobos. In 1998, a direct range measurement was made, which indicated that Phobos was slightly ahead of the predicted position. The MOLA detector views the surface of Mars in a narrow field of view, at 1064 nanometer wavelength, and can detect shadows cast by Phobos on the surface of Mars. We have found 15 such serendipitous shadow transit events over the interval from xx to xx, and all of them show Phobos to be ahead of schedule, and getting progressively farther ahead of the predicted position. In contrast, the cross-track positions are quite close to the predicted values. Assuming that the along-track discrepancy is small enough that we can linearize the corrections, we model the mean orbital longitude as a quadratic function of time, and solve for corrections to the constant, linear, and quadratic terms. The time span of the recent observations is insufficient to properly resolve this issue alone, but when the 127 years of prior observations are added, we find a solution which reduces misfit to the new observations considerably, and makes no significant change to the fit to earlier observations. Our estimate for the secular acceleration term is (1.367 + 0.006) degrees per square year. The corresponding rate of energy dissipation is 3.34 MW. From a geophysical perspective, a more interesting parameter than the secular acceleration itself is the tidal lag angle, or tidal quality factor Q, for Mars. Unfortunately, the limiting error source in that determination is remaining uncertainty in the tidal Love numbers at harmonic degrees 2,3, and even 4. Until those parameters are better constrained, improvement in the orbital model of Phobos will not provide corresponding benefits for understanding the interior of Mars.

Bills, Bruce↗

NASA NTRS2004

Modeling the Salar de Uyuni, Bolivia as an Equipotential Surface of Earth's Gravity Field

The salar de Uyuni is a massive dry salt lake that lies at the lowest point of an internal/drainage basin in the Bolivian Altiplano. Its topography is remarkable for its extraordinary flatness over almost a full degree of latitude and longitude. We surveyed a 54 x 45 km region of the salar with kinematic GPS in September, 2002 and found a topographic range of only 80 cm over the entire surveyed area. Furthermore, the survey revealed distinct surface features with several dominant wavelengths and orientations. Some of these appear to be aligned with orographic features that intersect the salar, leading us to conjecture that they are the surface expression of high-density mountains that have been buried by low-density basin sediments. Over the oceans, a similar correspondence between basin bathymetry and surface topography is exploited to map the seafloor using sea-surface satellite altimetry measurements, with the sea surface following geoid undulations due to the underwater mass distribution. On the salar, annual flooding creates a shallow lake whose surface also lies on a equipotential surface shaped by the distribution of underlying mass. The link to the actual salar surface is via the dissolution and redeposition of salt by the lake waters, which appears to push the system to an equilibrium of constant water depth and the coincidence of the shapes of the lake surface and bottom. To test our hypothesis about the origin of the surface features on the salar, we compare our GPS survey elevations with the equipotential surface generated from local gravity measurements in conjunction with gravity and potential values from the EGM96 global geopotential model. 50% of the variance of the GPS elevations can be explained by equipotential surface undulations from the EGM96 model alone, and an additional 40% is explained by the shorter-wavelength equipotential surface derived from local gravity. We examine the unexplained 10% of elevation variance from the standpoint of errors in the equipotential surface calculation and possible unmodelled surface processes.

Borsa, Adrian↗

NASA NTRS2004

Oceanic Tidal Mixing As a Contributor to Milankovitch-scale Climate Change

We propose that changes in the magnitude of oceanic tidal mixing on long time scales is an important, but previously unrecognized, contributor to global climate change. it is well known that Earth's orbital and rotational state changes significantly on 10(exp 4)-10(exp 5) year time scales, and that this influences the spatial and temporal pattern of incident radiation. It is widely supposed that climatic variations on these same time scales are, in large part, a response of the ocean-atmosphere-cryosphere system to this radiative forcing. Our proposal is that variations in the luni-solar tidal potential, induced by these same orbital and rotational variations, influences oceanic mixing and thus modulates meridional heat transport, by amounts which are competitive with the radiative forcing. There are some obvious differences between tidal potential and insolation. First is that the Sun and Moon both contribute to tides, whereas the radiation is entirely of solar origin. Second is that the Earth is transparent to gravity but opaque to radiation. Clipping associated with this opacity makes the radiation pattern temporal spectrum rather more complex than the tidal spectrum. A third point is that solar radiation directly delivers energy to Earth's surface whereas tidal mixing will only expedite lateral transport of heat in association with oceanic thermohaline circulation. The diurnal average insolation pattern is best parameterized via a Fourier series in time of year and Legendre polynomials in sine of latitude. Our present focus will be on the annual average terms. The Legendre degree n=0 term describes the global average insolation, and is nearly constant. The degree n=l term describes differences between northern and southern hemispheres, and the annual mean is zero. The degree n=2 term is the main contributor to the equator to pole variations, and varies with obliquity and orbital eccentricity, with the obliquity variation dominating. The lowest order decomposition of the tidal potential recognizes 3 constituents: semi-diurnal, diurnal, and long period, with solar and lunar contributions to each. Our present focus will be on long term variations in the mean square amplitude of the semi-diurnal constituent, with averaging over all the short period variations. For the solar tide that includes the day and year. For the lunar tide it includes the day, month, year, and the apsidal (8.85 year) and nodal (18.6 year) periods. We present calculations of the variations in radiative and tidal forcing for the past 3 million years. The two signals are quite similar. Both vary by approximately 1% of their respective mean values, are dominated by obliquity variations, and exhibit only secondary influence from orbital eccentricity.

Munk, Walter↗

NASA NTRS1993

An inversion of geoid and topography for mantle and crustal structure on Mars

Mars has the largest amplitude geoid anomalies and surface topography known on the terrestrial planets. A number of prior studies have analyzed Martian gravity anomalies and topography in terms of isostasy and flexure of the crust and lithosphere. Other studies have emphasized the role of mantle convection in producing gravity anomalies and topography in some regions of Mars. Geoid and topography observations for simultaneous estimates of density anomalies in the crust and mantle of Mars are inverted. In performing this study, a recent degree 50 spherical harmonic expansion of the Martian gravity field (GMM-l) and a corresponding resolution expansion of the USGS Mars topography model are used. However, our analysis is restricted to harmonic degrees up to L equals 25, which are better determined than the higher harmonics. This provides a half-wavelength horizontal resolution of 425 km.

Kiefer, Walter↗