The gravity of Mars is a natural phenomenon, due to the law of gravity, or gravitation, by which all things with mass around the planet Mars are brought towards it. It is weaker than Earth's gravity due to the planet's smaller mass. The average gravitational acceleration on Mars is 3.72076 ms^{−2} (about 38% of that of Earth) and it varies.^{[1]} In general, topographycontrolled isostasy drives the short wavelength freeair gravity anomalies.^{[2]} At the same time, convective flow and finite strength of the mantle lead to longwavelength planetaryscale freeair gravity anomalies over the entire planet.^{[3]}^{[4]} Variation in crustal thickness, magmatic and volcanic activities, impactinduced Mohouplift, seasonal variation of polar ice caps, atmospheric mass variation and variation of porosity of the crust could also correlate to the lateral variations.^{[5]}^{[6]}^{[7]}^{[8]}^{[9]} Over the years models consisting of an increasing but limited number of spherical harmonics have been produced. Maps produced have included freeair gravity anomaly, Bouguer gravity anomaly, and crustal thickness. In some areas of Mars there is a correlation between gravity anomalies and topography. Given the known topography, higher resolution gravity field can be inferred. Tidal deformation of Mars by the Sun or Phobos can be measured by its gravity. This reveals how stiff the interior is, and shows that the core is partially liquid. The study of surface gravity of Mars can therefore yield information about different features and provide beneficial information for future landing projects.
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,
where G is the universal constant of gravitation (commonly taken as G = 6.674 x 10^{−11} m^{3} kg^{−1} s^{−2}),^{[10]} M is the mass of Mars (most updated value: 6.41693 x 10^{23} kg),^{[11]} m is the mass of the satellite, r is the distance between Mars and the satellite, and is the angular velocity of the satellite, which is also equivalent to (T is the orbiting period of the satellite).
Therefore, , 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, nonspherical 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.
where are spherical coordinates of the test point.^{[12]} is longitude and is latitude. and are dimensionless harmonic coefficients of degree and order .^{[12]} is the Legendre polynomial of degree with and is the associated Legendre polynomial with . These are used to describe solutions of Laplace's equation.^{[12]} is the mean radius of the planet.^{[12]} The coefficient is sometimes written as .
The oldest technique in determining the gravity of Mars is through Earthbased observation. Later with the arrival of unmanned 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 Earthbased observations of the motions of Phobos and Deimos provide physical parameters including semimajor 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 oneway, twoway and threeway Doppler and range tracking applied.^{[16]} Oneway tracking means the data is transmitted in one way to the DSN from the spacecraft, while twoway and threeway involve transmitting signals from Earth to the spacecraft (uplink), and thereafter transponded coherently back to the Earth (downlink).^{[16]} The difference between twoway and threeway 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, nonconservative 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 (GMM3), 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 freeair 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 nonconservative forces applied to the spacecraft.^{[13]}
Gravity solutions  Authors  Year  Degree (m) and order (l) of the spherical harmonic solution
[Surface resolution (km)] 
Data source 

–  JP Gapcynski, RH Tolson and WH Michael Jr  1977  6^{[18]}  Tracking data of Mariner 9, Viking 1 and 2 spacecraft^{[18]} 
Geoide martien^{[19]}  G Balmino, B Moynot and N Vales  1982  18^{[19]}
[¬600 km] 
Tracking data of Mariner 9, Viking 1 and 2 spacecraft^{[19]} 
GMM1^{[20]}  DE Smith, FJ Lerch, RS Nerem, MT Zuber, GB Patel, SK Fricke and FG Lemoine  1993  50^{[20]}
[200–300 km] 
Tracking data of Mariner 9, Viking 1 and 2 spacecraft^{[20]} 
Mars50c^{[21]}  AS Konopliv, WL Sjogren  1995  50^{[21]}  Tracking data of Mariner 9, Viking 1 and 2 spacecraft^{[21]} 
GMM2B^{[14]}  FG Lemoine, DE Smith, DD Rowlands, MT Zuber, GA Neumann, DS Chinn, and DE Pavlis  2001  80^{[14]}  Tracking data of Mars Global Surveyor (MGS), and MOLAderived topography data ^{[14]} 
GGM1041C^{[22]}  FG Lemoine  2001  90^{[22]}  Tracking data of Mars Global Surveyor (MGS) and Mars Odyssey, and MOLAderived topography data^{[22]} 
MGS95J^{[23]}  AS Konopliv, CF Yoder, EM Standish, DN Yuan, WL Sjogren  2006  95^{[23]}
[~112 km] 
Tracking data of Mars Global Surveyor (MGS) and Mars Odyssey, and MOLAderived topography data ^{[23]} 
MGGM08A^{[7]}  JC Marty, G Balmino, J Duron, P Rosenblatt, S Le Maistre, A Rivoldini, V Dehant, T. Van Hoolst  2009  95^{[7]}
[~112 km] 
Tracking data of Mars Global Surveyor (MGS) and Mars Odyssey, and MOLAderived topography data^{[7]} 
MRO110B2^{[24]}  AS Konopliv, SW Asmar, WM Folkner, Ö Karatekin, DC Nunes, SE Smrekar, CF Yoder, MT Zuber  2011  110^{[24]}  Tracking data of Mars Global Surveyor (MGS), Mars Odyssey and Mars Reconnaissance Orbiter (MRO), and MOLAderived topography data^{[24]} 
MGM2011^{[1]}  C Hirt, SJ Claessens, M Kuhn, WE Featherstone  2012  [3 km (equator) – 125 km]^{[1]}  Gravity solution MRO110B2, and MOLAderived topography data^{[1]} 
GMM3^{[13]}  A Genova, S Goossens, FG Lemoine, E Mazarico, GA Neumann, DE Smith, MT Zuber  2016  120^{[13]}
[115 km] 
Mars Global Surveyor (MGS), Mars Odyssey and Mars Reconnaissance Orbiter (MRO)^{[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 deepprinted 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 Model1, or GMM1)^{[20]} then the subsequent models with higher completeness and degree and order up to 120 for the latest GMM3.^{[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 shortscale gravity model, utilizing the gravitytopography correlation in shortwavelength.^{[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 nonconservative 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 shortwavelength (locally varying) freeair gravity anomalies and topography. For regions with higher correlation, freeair 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 freeair gravity anomaly model typically has higher resolution for the southern hemisphere, as high as over 100 km.^{[13]}
Freeair 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 insitu analysis.
Since Bouguer gravity anomalies have strong link with depth of crustmantle 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 earlyNoachian buried surface.^{[5]}
Strong correlation between topography and shortwavelength freeair 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 (shortwavelength 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. small volcanoes Olympus Mons and the Tharsis Montes produce the smallest positive freeair 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^{3} 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 longwavelength gravity anomalies, are the lowdegree 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 lowdegree harmonics of the Earth's gravity field in order to understand its interior structure.^{[4]} It has been suggested that such longwavelength 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 longwavelength gravity anomalies, for example, the midAtlantic 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 sublimationcondensation 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,
^{[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 longwavelength 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 timevariable gravity due to such mass exchange,^{[24]}^{[8]}^{[31]}^{[32]} where is the degree while 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,
^{[32]}
^{[32]}
Data has indicated that the maximum mass variation of the southern polar cap is approximately 8.4 x 10^{15} kg,^{[13]} occurring near the autumnal equinox,^{[13]} while for that of the northern polar is approximately 6.2 x 10^{15} 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) x 10^{11} kg,^{[8]} while in South Pole it would decrease by (0.8 ± 0.6) x 10^{11} kg.^{[8]} In addition, the atmosphere would have decrease in term of the mass of carbon dioxide by (0.6 ± 0.6) x 10^{11} 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_{2} 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_{2} can tell quadrupolar deformation.^{[13]} Finding k_{2} is helpful in understanding the interior structure on Mars.^{[13]} The most updated k_{2} obtained by Genova's team is 0.1697 ± 0.0009.^{[13]} As if k_{2} 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 higherresolution 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,
^{[5]}
The residual Bouguer potential is contributed by the mantle.^{[5]} The undulation of the crustmantle 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 equalarea 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 latitudedependent.^{[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 preimpact 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 waterdeposited 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 selfgravitation. 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 subsolidus 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^{−3} 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^{−3},^{[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^{−3},^{[9]} 10% to 23% porosity could give a 200 kg m^{−3} 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 400km elevation orbit, could measure the range (distance) between the spacecraft and the ground surface through counting the roundtrip 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]}
As there is a large distance between Mars and Earth, immediate command to the lander is almost impossible and the landing relies highly on its autonomous system. It has been recognized that 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 manmade 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]}