Gravity of Mars

To understand the gravity of Mars, its gravitational field strength g and gravitational potential U are often measured. Simply, if Mars is assumed to be a static perfectly spherical body of radius R_M, provided that there is only one satellite revolving around Mars in a circular orbit and such gravitation interaction is the only force acting in the system, the equation would be

${\displaystyle {\frac {GMm}{r^{2}}}=mr\omega ^{2},}$ 

where G is the universal constant of gravitation (commonly taken as G = 6.674 × 10⁻¹¹ m³ kg⁻¹ s⁻²),^[10] M is the mass of Mars (most updated value: 6.41693 × 10²³ kg),^[11] m is the mass of the satellite, r is the distance between Mars and the satellite, and ${\displaystyle \omega }$  is the angular velocity of the satellite, which is also equivalent to ${\displaystyle {\frac {2\pi }{T}}}$  (T is the orbiting period of the satellite).

Therefore, ${\displaystyle g={\frac {GM}{R_{M}^{2}}}={\frac {r^{3}\omega ^{2}}{R_{M}^{2}}}={\frac {4r^{3}\pi ^{2}}{T^{2}R_{M}^{2}}}}$ , where R_M is the radius of Mars. With proper measurement, r, T, and R_M are obtainable parameters from Earth.

However, as Mars is a generic, non-spherical planetary body and influenced by complex geological processes, more accurately, the gravitational potential is described with spherical harmonic functions, following convention in geodesy; see Geopotential model.

${\displaystyle U(r,\lambda ,\psi )=-{\frac {GM}{r}}\left(1+\sum _{\ell =2}^{\ell =L}\left({\frac {R}{r}}\right)^{\ell }\left(C_{\ell 0}P_{\ell }^{0}(\sin \psi )+\sum _{m=1}^{+\ell }(C_{\ell m}\cos m\lambda +S_{\ell m}\sin m\lambda )P_{\ell }^{m}(\sin \psi )\right)\right),}$ ^[12]

where ${\displaystyle r,\psi ,\lambda }$  are spherical coordinates of the test point.^[12] ${\displaystyle \lambda }$  is longitude and ${\displaystyle \psi }$  is latitude. ${\displaystyle C_{\ell m}}$  and ${\displaystyle S_{\ell m}}$  are dimensionless harmonic coefficients of degree ${\displaystyle l}$  and order ${\displaystyle m}$ .^[12] ${\displaystyle P_{\ell }^{m}}$  is the Legendre polynomial of degree ${\displaystyle l}$  with ${\displaystyle m=0}$  and is the associated Legendre polynomial with ${\displaystyle m>0}$ . These are used to describe solutions of Laplace's equation.^[12] ${\displaystyle R}$  is the mean radius of the planet.^[12] The coefficient ${\displaystyle C_{\ell 0}}$  is sometimes written as ${\displaystyle J_{n}}$ .

1. The lower the degree ${\displaystyle \ell }$  and order ${\displaystyle m}$ , the longer wavelength of anomaly it represents. In turn, long-wavelength gravity anomaly is influenced by global geophysical structures.
2. The higher the degree ${\displaystyle \ell }$  and order ${\displaystyle m}$ , the shorter wavelength of anomaly it represents. For degree over 50, it has been shown that those variations have high correlation with the topography.^[13] Geophysical interpretation of surface features could further help deriving a more complete picture of the Martian gravity field, though misleading results could be produced.^[13]

The oldest technique in determining the gravity of Mars is through Earth-based observation. Later with the arrival of uncrewed spacecraft, subsequent gravity models were developed from radio tracking data.

Before the arrival of the Mariner 9 and Viking orbiter spacecraft at Mars, only an estimate of the Mars gravitational constant GM, i.e. the universal constant of gravitation times the mass of Mars, was available for deducing the properties of the Martian gravity field.^[14] GM could be obtained through observations of the motions of the natural satellites of Mars (Phobos and Deimos) and spacecraft flybys of Mars (Mariner 4 and Mariner 6).^[14]

Long term Earth-based observations of the motions of Phobos and Deimos provide physical parameters including semi-major axis, eccentricity, inclination angle to the Laplacian plane etc.,^[15] which allow calculation of the ratio of solar mass to the mass of Mars, moment of inertia and coefficient of the gravitational potential of Mars, and give initial estimates of the gravity field of Mars.^[15]

Precise tracking of spacecraft is of prime importance for accurate gravity modeling, as gravity models are developed from observing tiny perturbation of spacecraft, i.e. small variation in velocity and altitude. The tracking is done basically by the antennae of the Deep Space Network (DSN), with one-way, two-way and three-way Doppler and range tracking applied.^[16] One-way tracking means the data is transmitted in one way to the DSN from the spacecraft, while two-way and three-way involve transmitting signals from Earth to the spacecraft (uplink), and thereafter transponded coherently back to the Earth (downlink).^[16] The difference between two-way and three-way tracking is, the former one has the same signal transmitter and receiver on Earth, while the latter one has the transmitter and receiver at different locations on Earth.^[16] The use of these three types of tracking data enhances the coverage and quality of the data, as one could fill in the data gap of another.^[16]

Doppler tracking is a common technique in tracking the spacecraft, utilizing radial velocity method, which involves detection of Doppler shifts.^[13] As the spacecraft moves away from us along line of sight, there would be redshift of signal, while for the reverse, there would be blueshift of signal. Such technique has also been applied for observation of the motion of exoplanets.^[17] While for the range tracking, it is done through measurement of round trip propagation time of the signal.^[13] Combination of Doppler shift and range observation promotes higher tracking accuracy of the spacecraft.

The tracking data would then be converted to develop global gravity models using the spherical harmonic equation displayed above. However, further elimination of the effects due to affect of solid tide, various relativistic effects due to the Sun, Jupiter and Saturn, non-conservative forces (e.g. angular momentum desaturations (AMD), atmospheric drag and solar radiation pressure) have to be done,^[13] otherwise, considerable errors result.

The latest gravity model for Mars is the Goddard Mars Model 3 (GMM-3), produced in 2016, with spherical harmonics solution up to degree and order 120.^[13] This model is developed from 16 years of radio tracking data from Mars Global Surveyor (MGS), Mars Odyssey and Mars Reconnaissance Orbiter (MRO), as well as the MOLA topography model and provides a global resolution of 115 km.^[13] A separate free-air gravity anomaly map, Bouguer gravity anomaly map and a map of crustal thickness were produced along with this model.^[13] Compared with MRO110C and other previous models, major improvement of the estimation of the gravity field comes from more careful modeling of the non-conservative forces applied to the spacecraft.^[13]

The techniques in tracking the spacecraft and geophysical interpretation of surface features can affect the resolution of the strength of gravity field. The better technique favors spherical harmonic solutions to higher degrees and orders. Independent analysis on Mariner 9 and Viking Orbiter tracking data yielded a degree and order of 6 spherical harmonic solution.,^[18] Further combination of the two data sets, along with correlation of anomalies with volcanic features (positive anomaly) and deep-printed depression (negative anomaly) assisted by image data allows a degree and order of 18 spherical harmonic solution produced.^[19] Further use of spatial a priori constraint method, which had taken the topography into account in solving the Kaula power law constraint, had favored model of up to degree 50 spherical harmonic solution in global resolution (Goddard Mars Model-1, or GMM-1)^[20] then the subsequent models with higher completeness and degree and order up to 120 for the latest GMM-3.^[13]

Therefore, gravity models nowadays are not directly produced through transfer of the measured gravity data to any spatial information system because there is difficulty in producing model with sufficiently high resolution. Topography data obtained from the MOLA instrument aboard the Mars Global Surveyor thus becomes a useful tool in producing a more detailed short-scale gravity model, utilizing the gravity-topography correlation in short-wavelength.^[13] However, not all regions on Mars show such correlation, notably the northern lowland and the poles.^[13] Misleading results could be easily produced, which could lead to wrong geophysics interpretation.^[13]

The later modifications of gravity model include taking other non-conservative forces acting on spacecraft into account, including atmospheric drag, solar radiation pressure, Mars reflected solar radiation pressure, Mars thermal emission, and spacecraft thrusting which despins or desaturates the angular moment wheels.^[14] In addition, Martian precession and third body attraction due to the Sun, Moon and planets, which could affect the spacecraft orbit, as well as relavistic effects on the measurements should also be corrected.^[7] These factors could lead to offset of the true gravity field. Accurate modeling is thus required to eliminate the offset. Such work is still ongoing.

Many researchers have outlined the correlation between short-wavelength (locally varying) free-air gravity anomalies and topography. For regions with higher correlation, free-air gravity anomalies could be expanded to higher degree strength through geophysical interpretation of surface features,^[13] so that the gravity map could offer higher resolution. It has been found that the southern highland has high gravity/topography correlation but not for the northern lowland.^[13] Therefore, the resolution of free-air gravity anomaly model typically has higher resolution for the southern hemisphere, as high as over 100 km.^[13]

Free-air gravity anomalies are relatively easier to measure than the Bouguer anomalies as long as topography data is available because it does not need to eliminate the gravitational effect due to the effect of mass surplus or deficit of the terrain after the gravity is reduced to sea level. However, to interpret the crustal structure, further elimination of such gravitational effect is necessary so that the reduced gravity would only be the result of the core, mantle and crust below datum.^[5] The product after elimination is the Bouguer anomalies. However, density of the material in building up the terrain would be the most important constraint in the calculation, which may vary laterally on the planet and is affected by porosity and geochemistry of the rock.^[5][9] Relevant information could be obtained from Martian meteorites and in-situ analysis.

Since Bouguer gravity anomalies have strong links with depth of crust-mantle boundary, one with positive Bouguer anomalies may mean that it has a thinner crust composed of lower density material and is influenced more strongly by the denser mantle, and vice versa. However, it could also be contributed by the difference in density of the erupted volcanic load and sedimentary load, as well as subsurface intrusion and removal of material.^[5][6][25] Many of these anomalies are associated with either geological or topographic features.^[5] Few exception includes the 63°E, 71°N anomaly,^[5] which may represent an extensive buried structure as large as over 600 km, predated the early-Noachian buried surface.^[5]

Strong correlation between topography and short-wavelength free-air gravity anomalies has been shown for both study of the gravity field of the Earth and the Moon,^[2] and it can be explained by the wide occurrence of isostasy.^[2][26] High correlation is expected for degree over 50 (short-wavelength anomaly) on Mars.^[13] And it could be as high as 0.9 for degrees between 70 and 85.^[13] Such correlation could be explained by flexural compensation of topographic loads.^[2][26] It is noted that older regions on Mars are isostatically compensated when the younger region are usually only partially compensated.^[13]

Different volcanic constructs could behave differently in terms of gravity anomalies. Volcanoes Olympus Mons and the Tharsis Montes produce the smallest positive free-air gravity anomalies in the solar system.^[5] Alba Patera, also a volcanic rise, north of the Tharsis Montes, however, produces negative Bouguer anomaly, though its extension is similar to that of Olympus Mons.^[5] And for the Elysium Mons, its center is found to have slight increase in Bouguer anomalies in an overall broad negative anomaly context in the Elysium rise.^[5]

The knowledge of anomaly of volcanoes, along with density of the volcanic material, would be useful in determining the lithospheric composition and crustal evolution of different volcanic edifices.^[27] It has been suggested that the extruded lava could range from andesite (low density) to basaltic (high density) and the composition could change during the construction of the volcanic shield, which contributes to the anomaly.^[27] Another scenario is it is possible for high density material intruded beneath the volcano.^[27][6] Such setting has already been observed over the famous Syrtis major, which has been inferred to have an extinct magma chamber with 3300 kg m³ underlying the volcano, evident from positive Bouguer anomaly.^[6]

Different depressions also behave differently in Bouguer anomaly. Giant impact basins like Argyre, Isidis, Hellas and Utopia basins also exhibit very strong positive Bouguer anomalies in circular manner.^[5] These basins have been debated for their impact crater origin. If they are, the positive anomalies may be due to uplift of Moho, crustal thinning and modification events by sedimentary and volcanic surface loads after impacting.^[5][25]

But at the same time there are also some large basins that are not associated with such positive Bouguer anomaly, for example, Daedalia, northern Tharsis and Elysium, which are believed to be underlain by the northern lowland plain.^[5]

In addition, certain portions of Coprates, Eos Chasma and Kasei Valles are also found to have positive Bouguer anomalies,^[5] though they are topographic depressions. This may suggest that these depressions are underlain by shallow dense intrusion body.^[5]

Global gravity anomalies, also termed as long-wavelength gravity anomalies, are the low-degree harmonics of the gravity field,^[4] which cannot be attributed to local isostasy, but rather finite strength of the mantle and density differences in the convection current.^[13][3][4] For Mars, the largest component of Bouguer anomaly is the degree one harmonic, which represents the mass deficit in the southern hemisphere and excess in the northern hemisphere.^[5] The second largest component corresponds to the planet flattening and Tharsis bulge.^[5]

Early study of the geoid in the 1950s and 1960s has focused on the low-degree harmonics of the Earth's gravity field in order to understand its interior structure.^[4] It has been suggested that such long-wavelength anomalies on Earth could be contributed by the sources located in deep mantle and not in the crust, for example, caused by the density differences in driving the convection current,^[4][28] which has been evolving with time. The correlation between certain topography anomalies and long-wavelength gravity anomalies, for example, the mid-Atlantic ridge and Carlsberg ridge, which are topography high and gravity high on the ocean floor, thus became the argument for the convection current idea on Earth in the 1970s,^[29][30] though such correlations are weak in the global picture.

Another possible explanation for the global scale anomalies is the finite strength of the mantle (in contrast to zero stress), which makes the gravity deviated from hydrostatic equilibrium.^[3] For this theory, because of the finite strength, flow may not exist for most of the region that are understressed.^[3] And the variations of density of the deep mantle could be the result of chemical inhomogeneities associated with continent separations,^[3] and scars left on Earth after the torn away of the moon.^[3] These are the cases suggested to work when slow flow is allowed to happen under certain circumstances.^[3] However, it has been argued that the theory may not be physically feasible.^[4]

The sublimation-condensation cycle of carbon dioxide on Mars between the atmosphere and cryosphere (polar ice cap) operates seasonally.^[8] This cycle contributes as almost the only variable accounting for changes in gravity field on Mars.^[8] The measured gravitational potential of Mars from orbiters could be generalized as the equation below,

${\displaystyle V({\text{Mars}})=V({\text{solid planet}})+V({\text{seasonal caps}}+{\text{atmosphere}})}$ ^[8]

In turn, when there is more mass in the seasonal caps due to the more condensation of carbon dioxide from the atmosphere, the mass of the atmosphere would drop. They have inverse relationship with each other. And the change in mass has direct effect towards the measured gravitational potential.

The seasonal mass exchange between the northern polar cap and southern polar cap exhibits long-wavelength gravity variation with time.^[8][13] Long years of continuous observation has found that the determination of even zonal, normalized gravity coefficient C_{l=2, m=0}, and odd zonal, normalized gravity coefficient C_{l=3, m=0} are crucial for outlining the time-variable gravity due to such mass exchange,^{[24][8][31][32]} where ${\displaystyle l}$  is the degree while ${\displaystyle m}$  is the order. More commonly, they are represented in form of C_lm in research papers.

If we regard the two poles as two distinct point masses, then, their masses are defined as,

${\displaystyle M_{NP}={\frac {C_{20}+C_{30}}{2}}\,M_{\text{Mars}}}$ ^[32]

${\displaystyle M_{SP}={\frac {C_{20}-C_{30}}{2}}\,M_{\text{Mars}}}$ ^[32]

Data has indicated that the maximum mass variation of the southern polar cap is approximately 8.4 × 10¹⁵ kg,^[13] occurring near the autumnal equinox,^[13] while for that of the northern polar is approximately 6.2 × 10¹⁵ kg,^[13] occurring in between the winter solstice and spring equinox.^[13]

In long term speaking, it has been found that the mass of ice stored in North Pole would increase by (1.4 ± 0.5) × 10¹¹ kg,^[8] while in South Pole it would decrease by (0.8 ± 0.6) × 10¹¹ kg.^[8] In addition, the atmosphere would have decrease in term of the mass of carbon dioxide by (0.6 ± 0.6) × 10¹¹ kg in long term as well.^[8] Due to existence of uncertainties, it is unclear whether migration of material from the South Pole to the North Pole is ongoing, though such a possibility cannot be ruled out.^[8]

The two major tidal forces acting on Mars are the solar tide and Phobos tide.^[13] Love number k₂ is an important proportional dimensionless constant relating the tidal field acting to the body with the multipolar moment resulting from the mass distribution of the body. Usually k₂ can tell quadrupolar deformation.^[13] Finding k₂ is helpful in understanding the interior structure on Mars.^[13] The most updated k₂ obtained by Genova's team is 0.1697 ± 0.0009.^[13] As if k₂ is smaller than 0.10 a solid core would be indicated, this tells that at least the outer core is liquid on Mars,^[31] and the predicted core radius is 1520–1840 km.^[31]

However, current radio tracking data from MGS, ODY and MRO does not allow the effect of phase lag on the tides to be detected because it is too weak and needs more precise measurement on the perturbation of spacecraft in the future.^[13]

No direct measurement of crustal thickness on Mars is currently available. Geochemical implications from SNC meteorites and orthopyroxenite meteorite ALH84001 suggested that mean crustal thickness of Mars is 100–250 km.^[33] Viscous relaxation analysis suggested that the maximum thickness is 50–100 km. Such thickness is critical in maintaining hemispheric crustal variations and preventing channel flow.^[34] Combination studies on geophysics and geochemistry suggested that average crustal thickness could be down to 50 ± 12 km.^[35]

Measurement of gravity field by different orbiters allows higher-resolution global Bouguer potential model to be produced.^[5] With local shallow density anomalies and effect of core flattening eliminated,^[5] the residual Bouguer potential is produced, as indicated by the following equation:

${\displaystyle U_{RB}=U_{B}-U_{\text{core}}-U_{\text{local}}}$ ^[5]

The residual Bouguer potential is contributed by the mantle.^[5] The undulation of the crust-mantle boundary, or the Moho surface, with mass of terrain corrected, should have resulted in varying residual anomaly.^[5] In turn, if undulating boundary is observed, there should be changes in crustal thickness.

Global study of residual Bouguer anomaly data indicates that crustal thickness of Mars varies from 5.8 km to 102 km.^[5] Two major peaks at 32 km and 58 km are identified from an equal-area histogram of crustal thickness.^[5] These two peaks are linked to the crustal dichotomy of Mars.^[5] Almost all the crust thicker than 60 km are contributed by the southern highland, with generally uniform thickness.^[5] And the northern lowland in general has thinner crust. The crustal thickness of the Arabia Terra region and northern hemisphere are found to be latitude-dependent.^[5] The more southward towards the Sinai Planum and Lunae Planum, the more thickened the crust is.^[5]

Among all regions, the Thaumasia and Claritis contain the thickest portion of crust on Mars that account for the histogram > 70 km.^[5] The Hellas and Argyre basins are observed to have crust thinner than 30 km,^[5] which are the exceptionally thin area in the southern hemisphere.^[5] Isidis and Utopia are also observed to have significant crustal thinning,^[5] with the center of Isidis basins believed to have the thinnest crust on Mars.^[5]

Crustal thinning is believed to have taken place underneath almost all the major impact craters.^[5] Crustal excavation, modification through emplacement of volcanic material and crustal flow taking place in the weak lithosphere are the possible causes.^[5] With the pre-impact crust excavated, gravitational restoration would take place through central mantle uplift, so that the mass deficit of cavity could be compensated by the mass of the uplifted denser material.^[5]

Giant impact basins Utopia, Hellas, Argyre and Isidis are some of the most prominent examples.^[5] Utopia, an impact basin located in northern lowland, is filled by light and water-deposited sedimentary material and has slightly thickened crust at the center.^[5] This is potentially due to large resurfacing process in the northern lowland.^[5] While for Hellas, Argyre and Isidis basins, they have great Moho uplifted relief and exhibit annuli of diffuse thickened crust beyond the crustal rim.^[5]

But on the contrary, almost all the Martian basins with diameter of 275 km < D < 1000 km are associated with low amplitude surface and low amplitude Moho relief.^[25] Many are even found to have negative free air gravity anomaly, though evidence has shown that all of them should have experienced gravity high (positive free air gravity anomaly).^[25] These have been suggested not caused by erosion and burial alone, as the adding of material into the basin would in fact increase the gravity strength rather than decrease it.^[25] Thus viscous relaxation should have been taking place.^[25] High heat flux and high water content in the early Martian crust favored viscous relaxation.^[25] These two factors have made the crust more ductile. The basin topography of the craters would be subjected to greater stress due to self-gravitation. Such stress would drive crustal flow and therefore decay of relief. The giant impact basins are the exceptions that have not experienced viscous relaxation, as crustal thinning has made the crust too thin to sustain sub-solidus crustal flow.^[5][25]

The most recent crustal density model RM1 developed in 2017 gives the bulk crustal density to be 2582 ± 209 kg m⁻³ for Mars,^[9] which represents a global average value.^[9] Lateral variation of the crustal density should exist.^[9] For example, over the volcanic complexes, local density is expected to be as high as 3231 ± 95 kg m⁻³,^[9] which matched the meteorite data and previous estimations. In addition, the density of the northern hemisphere is in general higher than that of the southern hemisphere,^[9] which may imply that the latter is more porous than the former.

To achieve the bulk value, porosity could play an important role. If the mineral grain density is chosen to be 3100 kg m⁻³,^[9] 10% to 23% porosity could give a 200 kg m⁻³ drop in the bulk density.^[9] If the pore spaces are filled with water or ice, bulk density decrease is also expected.^[9] A further drop in bulk density could be explained by increasing density with depth,^[9] with the surface layer more porous than the deeper Mars, and the increase of density with depth also has geographical variation.^[9]

The areoid is a planetary geoid that represents the gravitational and rotational equipotential figure of Mars, analogous to the concept of geoid ("sea level") on Earth.^[5][36][37] This has been set as the reference frame for developing the MOLA Mission Experiment Gridded Data Records (MEGDRs),^[5][36] which is a global topography model. The topography model is important in mapping the geomorphological features and understanding different kinds of processes on Mars.

To derive the areoid, two parts of works are required. First, as gravity data is essential for identifying the position of the center of mass of the planet,^[36] which is largely affected by the distribution of the mass of the interior, radio tracking data of spacecraft is necessary.^[36] This was largely done by the Mars Global Surveyor (MGS).^[5][36] Then, the MOLA 2 instrument aboard the MGS, which operates at 400-km elevation orbit, could measure the range (distance) between the spacecraft and the ground surface through counting the round-trip time of flight of the pulse from the instrument.^[36] Combination of these two works allows the areoid as well as the MEGDRs to be constructed. Based on the above, the areoid has taken the radius as the mean radius of the planet at the equator as 3396 km.^[5][36]

There is a large distance between Mars and Earth, so immediate command to a lander is almost impossible. A landing relies highly on its autonomous system. To avoid failure, precise understanding of the gravity field of Mars is essential for the landing projects, so that offsetting factors and uncertainties of gravitational effects could be minimized, allowing for a smooth landing progress.^[38][39] The first ever man-made object landing on Mars, the Mars 2 lander, crashed for an unknown reason. Since the surface environment of Mars is complex, composed of laterally varying morphological patterns, in order to avoid rock hazard the landing progress should be further assisted by employment of LIDAR on site in determining the exact landing position and other protective measures.^[38][39]

• Geology of Mars
• Mars carbonate catastrophe