Solar radiation incident on Mars and the outer planets: Latitudinal,seasonal, and atmospheric effects

Calculations of the daily solar radiation incident at the tops of the atmospheres of Mars and the outer planets and its variability with latitude and season are presented in a series of figures and tables similar to those for Earth in The Smithsonian Meteorological Tables. The changes in the latitudinal and seasonal distributions of daily surface insolation during the great Martian dust storm of 1971 (when Martian atmospheric optical depth increased from about τ = 0.1 to 2.0 were significant and dramatically illustrate the effect of atmospheric aerosols on surface insolation; i.e., the mean annual daily insolation at the poles decreased by more than a factor of 100 as τ increased from 0.1 to 2.0.     

Lidar atmospheric measurements on Mars and Earth

The LIDAR instrument operating from the surface of Mars on the Phoenix Mission measured vertical profiles of atmospheric dust and water ice clouds at temperatures around −65 °C. An equivalent lidar system was utilized for measurements in the atmosphere of Earth where dust and cloud conditions are similar to Mars. Coordinated aircraft in situ sampling provided a verification of lidar measurement and analysis methods and also insight for interpretation of lidar derived optical parameters in terms of the dust and cloud microphysical properties. It was found that the vertical distribution of airborne dust above the Australian desert is quite similar to what is observed in the planetary boundary layer above Mars. Comparison with the in situ sampling is used to demonstrate how the lidar derived optical extinction coefficient is related to the dust particle size distribution. The lidar measurement placed a constraint on the model size distribution that has been used for Mars. Airborne lidar measurements were also conducted to study cirrus clouds that form in the Earth’s atmosphere at a similar temperature and humidity as the clouds observed with the lidar on Mars. Comparison with the in situ sampling provides a method to derive the cloud ice water content (IWC) from the Mars lidar measurements.     

Seasonal buffering of atmospheric pressure on Mars

An isothermal reservoir of carbon dioxide in gaseous contact with the Martian atmosphere would reduce the amplitude and advance the phase of global atmospheric pressure fluctuations caused by seasonal growth and decline of polar CO₂ frost caps. Adsorbed carbon dioxide in the upper ～10 m of Martian regolith is sufficient to buffer the present atmosphere on a seasonal basis. Available observations and related polar cap models do not confirm or refute the operation of such a mechanism. Implications for the amplitude and phase of seasonal pressure fluctuations are subject to direct test by the upcoming Viking mission to Mars.     

Assimilation of thermal emission spectrometer atmospheric data during the Mars Global Surveyor aerobraking period

The Thermal Emission Spectrometer aboard the Mars Global Surveyor spacecraft has produced an extensive atmospheric data set, beginning during aerobraking and continuing throughout the extended scientific mapping phase. Temperature profiles for the atmosphere below about 40 km, surface temperatures and total dust and water ice opacities, can be retrieved from infrared spectra in nadir viewing mode. This paper describes assimilation of nadir retrievals from the spacecraft aerobraking period, L_S=190°–260°, northern hemisphere autumn to winter, into a Mars general circulation model. The assimilation scheme is able to combine information from temperature and dust optical depth retrievals, making use of a model forecast containing information from the assimilation of earlier observations, to obtain a global, time-dependent analysis. Given sufficient temperature retrievals, the assimilation procedure indicates errors in the a priori dust distribution assumptions even when lacking dust observations; in this case there are relatively cold regions above the poles compared to a model which assumes a horizontally-uniform dust distribution. One major reason for using assimilation techniques is in order to investigate the transient wave behavior on Mars. Whilst the data from the 2-h spacecraft mapping orbit phase is much more suitable for assimilation, even the longer (45–24 h) period aerobraking orbit data contain useful information about the three-dimensional synoptic-scale martian circulation which the assimilation procedure can reconstruct in a consistent way. Assimilations from the period of the Noachis regional dust storm demonstrate that the combined assimilation of temperature and dust retrievals has a beneficial impact on the atmospheric analysis.     

Numerical simulation of atmospheric bore waves on Mars

The Viking Orbiters imaged early morning, long, linear wave clouds along the flanks of the Tharsis volcanoes during late northern spring and early summer. These clouds are believed to be a product of either an atmospheric bore wave or a hydraulic jump generated by nightly katabatic winds. The Mars Regional Atmospheric Modeling System was used to study the interaction of the katabatic flows with the surrounding atmosphere to determine what mechanism is responsible for the clouds. Simulations at Ls=90°, 100°, 142°, 180°, 270°, and 358° were conducted focusing on the eastern flank of Olympus Mons. Model results compare well with Viking observations and closely approximate theoretical treatments of atmospheric bores. Strong downslope flows are simulated during the night, with a bore wave forming on and behind a well-defined katabatic front. The observed seasonality of the clouds was reproduced in the simulations; the bore was deeper and faster during northern summer and weakest during the winter. When the bore was strong, it was undular in form, and generated vertically propagating gravity waves in the atmosphere above. During the winter, the atmospheric structure was such that any gravity waves generated damped with height. Less atmospheric water vapor abundance during northern winter, as compared to the summer, is also a factor in the seasonality of the wave clouds. This study concludes that bore waves are the most likely mechanism for the generation of the observed linear wave clouds.     

Global patterns in cloud forms on Mars

Mariner 9 images and all Viking orbiter images through July 1979 were searched for cloud forms. A computer-accessible catalog was assembled, consisting of a classification of cloud type (lee wav wave, for example) and properties (directionality, wavelength, for example). Lee wave directionality shows a pattern and seasonal variation at high latitudes which is consistent with predictions of theoretical modeling. Fog and haze occurrence shows no obvious correlation with water abundance or any other simple causal factor. Lee waves are rare at equatorial latitudes. Plumes (probably dust) occur preferentially at locations where strong boundary layer convection is expected.     

Global and regional ridge patterns on Mars

Over 16,000 wrinkle ridges on Mars have been classified, mapped, and digitized to provide a large computer-accessible data base for analyzing regions subjected to possible compressive stresses. The survey has revealed major compressive structures that occur well beyond the Tharsis-dominated hemisphere. The large variety of affected geologic terrains indicate stresses not simply localized in the intercrater plains. One major area of inferred compression occurs in the southern cratered highlands near longitude 180°W where major ridges and scarps extend over 3000 km. The occurrence and orientation of many ridges are locally controlled by ancient impact basins. The Chryse basin in particular has an important effect on ridges in northern Lunae Planum. The removal of all basin-concentric ridges reveals, however, a complex global pattern. Although such patterns may yet be controlled by heavily degraded impact basins, major regional trends also emerge that appear to require broader scale global stresses. Most ridges in the Western Hemisphere are shown to be orthogonal to three centers corresponding to Tharsis (1°N, 122°W), Syria (12°S, 104°W), and Chryse (19°N, 47°W). Ridges not included in these three sets are generally more random and highly localized. Most, but by no means all, ridges in the Eastern Hemisphere are controlled by the Hellas and Isidis basins. A simple global grid is not yet identifiable and may not be recognizable owing to the large number of regional patterns.     

Non-uniform seasonal defrosting of subpolar dune field on Mars

Images of the sub-polar regions of Mars show existence of various seasonal-depending features, including mosaics of dark patches called Dalmatian spots. We interpret formation of dark patches on dune fields by means of the non-uniform sublimation of winter CO₂ ice deposit.     

The atmospheric temperatures over Olympus Mons on Mars: An atmospheric hot ring

We study the thermal fields over Olympus Mons separating seasons (northern spring and summer against southern spring and summer) and local time observations (day side versus night side). Temperature vertical profiles retrieved from Planetary Fourier Spectrometer on board Mars Express (PFS-MEX) data have been used. In many orbits (running north to south along a meridian) passing over the top of the volcano there is evidence of a hot air on top of the volcano, of two enhancement of the air temperature both north and south of it and in between a collar of air that is colder than nearby at low altitudes, and warmer than nearby at high altitudes. Mapping together many orbits passing over or near the volcano we find that the hot air has the tendency to form an hot ring around it. This hot structure occurs mostly between LT = 10.00 and 15.00 and during the northern summer. Distance of the hot structure from the top of the volcano is about 600 km (10° of latitude). The hot atmospheric region is 300-420 km (5-7°) wide. Hot ring temperature contrasts of about 40 K occur at 2 km above the surface and decrease to 20 K at 5 km and to 10 K at 10 km. The atmospheric circulation over an area of 40° × 40° (latitudes and longitudes) is affected by the topography and the presence of Olympus Mons (−133°W, 18°N). We discuss also the thermal stability of the atmosphere over the selected area using the potential temperatures. The temperature field over the top of the volcano shows unstable atmosphere within 10 km from the surface. Finally, we interpret the hot temperatures around volcano as an adiabatic compression of down-welling branch coming from over the top of volcano.Different air temperature profiles are observed in the same seasons during the night, or in different seasons. In northern spring-summer during the night the isothermal contours do not show the presence of the volcano until we reach close to the surface very much, where a thermal inversion is observed. The surface temperature shows higher values (by 10 K) in correspondence of the scarp (an abrupt altimetry variation of roughly 5 km) on south (6°N) and north (30°N) sides of volcano. During the southern spring-summer, on the contrary the isothermal curves run parallel to the surface even on top the volcano, just like the GCM have predicted.     

Albedo control of seasonal South Polar cap recession on Mars

Over the last few decades, General Circulation Models (GCM) have been used to simulate the current martian climate. The calibration of these GCMs with the current seasonal cycle is a crucial step in understanding the climate history of Mars. One of the main climatic signals currently used to validate GCMs is the annual atmospheric pressure cycle. It is difficult to use changes in seasonal deposits on the surface of Mars to calibrate the GCMs given the spectral ambiguities between CO₂ and H₂O ice in the visible range. With the OMEGA imaging spectrometer covering the near infra-red range, it is now possible to monitor both types of ice at a spatial resolution of about 1 km. At global scale, we determine the change with time of the Seasonal South Polar Cap (SSPC) crocus line, defining the edge of CO₂ deposits. This crocus line is not symmetric around the geographic South Pole. At local scale, we introduce the snowdrop distance, describing the local structure of the SSPC edge. Crocus line and snowdrop distance changes can now be used to calibrate GCMs. The albedo of the seasonal deposits is usually assumed to be a uniform and constant parameter of the GCMs. In this study, albedo is found to be the main parameter controlling the SSPC recession at both global and local scale. Using a defrost mass balance model (referred to as D-frost) that incorporates the effect of shadowing induced by topography, we show that the global SSPC asymmetry in the crocus line is controlled by albedo variations. At local scale, we show that the snowdrop distance is correlated with the albedo variability. Further GCM improvements should take into account these two results. We propose several possibilities for the origin of the asymmetric albedo control. The next step will be to identify and model the physical processes that create the albedo differences.     

Hydrogen peroxide on Mars: evidence for spatial and seasonal variations

Hydrogen peroxide (H₂O₂) has been suggested as a possible oxidizer of the martian surface. Photochemical models predict a mean column density in the range of 10¹⁵-10¹⁶ cm⁻². However, a stringent upper limit of the H₂O₂ abundance on Mars (9×10¹⁴ cm⁻²) was derived in February 2001 from ground-based infrared spectroscopy, at a time corresponding to a maximum water vapor abundance in the northern summer (30 pr. μm, Ls=112°). Here we report the detection of H₂O₂ on Mars in June 2003, and its mapping over the martian disk using the same technique, during the southern spring (Ls=206°) when the global water vapor abundance was ∼10 pr. μm. The spatial distribution of H₂O₂ shows a maximum in the morning around the sub-solar latitude. The mean H₂O₂ column density (6×10¹⁵ cm⁻²) is significantly greater than our previous upper limit, pointing to seasonal variations. Our new result is globally consistent with the predictions of photochemical models, and also with submillimeter ground-based measurements obtained in September 2003 (Ls=254°), averaged over the martian disk (Clancy et al., 2004, Icarus 168, 116-121).     

Mars: North-polar atmospheric warming during dust storms

Measurements of winter polar atmospheric temperatures on Mars are presented for the two 1977 dust storms observed by Viking. These data show a strong time and spatial dependence of the polar warming and provide important boundary conditions for dynamical models of the phenomenon.     

Mars atmospheric water vapor abundance: 1996-1997

Measurements of martian atmospheric water vapor made throughout L_s = 18.0°-146.4° (October 3, 1996-July 12, 1997) show changes in Mars humidity on hourly, daily, and seasonal time scales. Because our observing program during the 1996-1997 Mars apparition did not include concomitant measurement of nearby CO₂ bands, high northern latitude data were corrected for dust and aerosol extinction assuming an optical depth of 0.8, consistent with ground-based and HST imaging of northern dust storms. All other measurements with airmass greater than 3.5 were corrected using a total optical depth of 0.5. Three dominant results from this data set are as follows: (1) pre- and post-opposition measurements made with the slit crossing many hours of local time on Mars’ Earth-facing disk show a distinct diurnal pattern with highest abundances around and slightly after noon with low abundances in the late afternoon, (2) measurements of water vapor over the Mars Pathfinder landing site (Carl Sagan Memorial Station) on July 12, 1997, found 21 ppt μm in the spatial sector centered near 19° latitude, 36° longitude while abundances around the site varied from as low as 6 to as high as 28 ppt μm, and (3) water vapor abundance is patchy on hourly and daily time scales but follows the usual seasonal trends.     

Mars high resolution gravity fields from MRO, Mars seasonal gravity, and other dynamical parameters

With 2 years of tracking data collection from the MRO spacecraft, there is noticeable improvement in the high frequency portion of the spherical harmonic Mars gravity field. The new JPL Mars gravity fields, MRO110B and MRO110B2, show resolution near degree 90. Additional years of MGS and Mars Odyssey tracking data result in improvement for the seasonal gravity changes which compares well to global circulation models and Odyssey neutron data and Mars rotation and precession (). Once atmospheric dust is accounted for in the spacecraft solar pressure model, solutions for Mars solar tide are consistent between data sets and show slightly larger values (k₂ = 0.164 ± 0.009, after correction for atmospheric tide) compared to previous results, further constraining core models. An additional 4 years of Mars range data improves the Mars ephemeris, determines 21 asteroid masses and bounds solar mass loss (dGM_Sun/dt < 1.6 × 10⁻¹³ GM_Sun year⁻¹).     

Retrieval algorithm for atmospheric dust properties from Mars Global Surveyor Thermal Emission Spectrometer data during global dust storm 2001A

We present a new parameter retrieval algorithm for Mars Global Surveyor Thermal Emission Spectrometer data. The algorithm uses Newtonian first-order sensitivity functions of the infrared spectrum in response to variations in physical parameters to fit a model spectrum to the data at 499, 1099, and 1301 cm⁻¹. The algorithm iteratively fits the model spectrum to data to simultaneously retrieve dust extinction optical depth, effective radius, and surface temperature. There are several sources of uncertainty in the results. The assumed dust vertical distribution can introduce errors in retrieved optical depth of a few tens of percent. The assumed dust optical constants can introduce errors in both optical depth and effective radius, although the systematic nature of these errors will not affect retrieval of trends in these parameters. The algorithm does not include the spectral signature of water ice, and hence data needs to be filtered against this parameter before the algorithm is applied. The algorithm also needs sufficient dust spectral signature, and hence surface-to-atmosphere temperature contrast, to successfully retrieve the parameters. After the application of data filters the algorithm is both relatively accurate and very fast, successfully retrieving parameters, as well as meaningful parameter variability and trends from tens of thousands of individual spectra on a global scale (Elteto, A., Toon, O.B. [2010]. Icarus, this issue). Our results for optical depth compare well with TES archive values when corrected by the single scattering albedo. Our results are on average 1–4 K higher in surface temperatures from the TES archive values, with greater differences at higher optical depths. Our retrieval of dust effective radii compare well with the retrievals of Wolff and Clancy (Wolff, M.J., Clancy, R.T. [2003]. J. Geophys. Res. 108 (E9), 5097) for the corresponding data selections from the same orbits.     

Mars Global Surveyor Thermal Emission Spectrometer (TES) observations of variations in atmospheric dust optical depth over cold surfaces

The Mars Global Surveyor Thermal Emission Spectrometer (TES) instrument has returned over 200 million thermal infrared spectra of Mars taken between March 1999 and August 2004. This represents one of the most complete records of spatial and temporal changes of the martian atmosphere ever recorded by an orbiting spacecraft. Previous reports of the standard TES retrieval of aerosol optical depth have been limited to those observations taken over surfaces with temperatures above 210 K, limiting the spatiotemporal coverage of Polar Regions with TES. Here, we present an extension to the standard TES retrieval that better models the effects of cold surfaces below 200 K. This modification allows aerosol optical depth to be retrieved from TES spectra over a greater spatiotemporal range than was previously possible, specifically in Polar Regions. This new algorithm is applied to the Polar Regions to show the seasonal variability in dust and ice optical depth for the complete temporal range of the TES database (Mars Year 24, Ls=104°, 1 March 1999 to Mars Year 24, Ls=82°, 31 August 2004).     

Global atmospheric monitoring with SCIAMACHY

SCIAMACHY (SCanning Imaging Absorption spectrometer for Atmospheric CHartographY) is a space based spectrometer designed to measure sunlight transmitted, reflected and scattered by the Earth atmosphere or surface. It is a contribution to the Envisat-1 satellite to be launched in late 1999.SCIAMACHY measurements will provide amounts and distribution of 0₃, BrO, OCl0, ClO, S0₂, H₂CO, N0₂, CO, CO₂, CH₄, H₂O, N₂0, pressure, temperature, aerosol, radiation, cloud cover and cloud top height from atmospheric measurements in nadir, limb and occultation geometry.By the combination of the near simultaneous limb and nadir observations SCIAMACHY is one of a limited number of instruments which is able to detect tropospheric column amounts of 0₃, N0₂, CO, CH₄, H₂O, N₂0, S0₂, H₂CO, and BrO down to the planetary boundary layer under cloud free conditions.     

The atmospheric circulation and dust activity in different orbital epochs on Mars

A general circulation model is used to evaluate changes to the circulation and dust transport in the martian atmosphere for a range of past orbital conditions. A dust transport scheme, including parameterized dust lifting, is incorporated within the model to enable passive or radiatively active dust transport. The focus is on changes which relate to surface features, as these may potentially be verified by observations. Obliquity variations have the largest impact, as they affect the latitudinal distribution of solar heating. At low obliquities permanent CO₂ ice caps form at both poles, lowering mean surface pressures. At higher obliquities, solar insolation peaks at higher summer latitudes near solstice, producing a stronger, broader meridional circulation and a larger seasonal CO₂ ice cap in winter. Near-surface winds associated with the main meridional circulation intensify and extend polewards, with changes in cap edge position also affecting the flow. Hence the model predicts significant changes in surface wind directions as well as magnitudes. Dust lifting by wind stress increases with obliquity as the meridional circulation and associated near-surface winds strengthen. If active dust transport is used, then lifting rates increase further in response to the larger atmospheric dust opacities (hence circulation) produced. Dust lifting by dust devils increases more gradually with obliquity, having a weaker link to the meridional circulation. The primary effect of varying eccentricity is to change the impact of varying the areocentric longitude of perihelion, l, which determines when the solar forcing is strongest. The atmospheric circulation is stronger when l aligns with solstice rather than equinox, and there is also a bias from the martian topography, resulting in the strongest circulations when perihelion is at northern winter solstice. Net dust accumulation depends on both lifting and deposition. Dust which has been well mixed within the atmosphere is deposited preferentially over high topography. For wind stress lifting, the combination produces peak net removal within western boundary currents and southern midlatitude bands, and net accumulation concentrated in Arabia and Tharsis. In active dust transport experiments, dust is also scoured from northern midlatitudes during winter, further confining peak accumulation to equatorial regions. As obliquity increases, polar accumulation rates increase for wind stress lifting and are largest for high eccentricities when perihelion occurs during northern winter. For dust devil lifting, polar accumulation rates increase (though less rapidly) with obliquity above o=25°, but increase with decreasing obliquity below this, thus polar dust accumulation at low obliquities may be increasingly due to dust lifted by dust devils. For all cases discussed, the pole receiving most dust shifts from north to south as obliquity is increased.     

Mars: Photodesorption from mineral surfaces and its effects on atmospheric stability

A mechanism has been proposed for uv-accelerated desorption from Fe²⁺ sites on mineral surfaces that satisfies kinetic constraints determined in the laboratory by Huguenin. The process is an integral step of the photochemical weathering mechanism for producing dust on Mars, and it now appears that it may play primary roles in stabilizing CO₂ against dissociation by sunlight and in controlling the oxidation state of the atmosphere. We propose that adsorption occurs at octahedrally coordinated Fe²⁺ surface sites to form seven-coordinate transition-state complexes. These complexes acquire 16–18 kcal mole^?1 of ligand field stabilization energy. During illumination (λ ≤ 0.35 μm), electrons are photoemitted from the surfaced Fe²⁺, temporarily oxidizing them to Fe³⁺. Fe³⁺ has no ligand field stabilization energy, and the complexes lose 16–18 kcal mole^?1 of stabilization energy. This is a large fraction of the 19- to 28-kcal mole^?1 activation energy for dissociating the complexes, and desorption should proceed spontaneously. The gases that were observed to undergo adsorption-photodesorption include O₂, CO₂, CO, H₂O, N₂, and Ar. Photodesorption can drive several catalytic reactions, one of which is the oxidation of CO to CO₂. The rate of this reaction should be limited by the supply of CO and O₂ to the surface to ～2 × 10¹² cm^?2 sec^?1 (column photodissociation rate of CO₂). By including this surface reaction in models of Martian atmospheric CO₂ chemistry, CO₂ can be stabilized against photodissociation with eddy diffusion coefficients of only 3 × 10⁵?1 × 10⁷ cm² sec^?1 below 40 km, raising to ～ 10⁹ cm² sec^?1 at 140 km. Odd hydrogen is not needed to catalyze the oxidation of CO below 40 km, and odd hydrogen mixing ratios need only to be $\text{f}{}_{\mathrm{<mtext>H</mtext>}}\text{? 10}{}^{\mathrm{?10}}$ to depress ozone concentrations below the observed upper limit in equatorial regions. Another catalytic reaction that should be driven by photodesorption on Mars is $\text{20}\text{H}{}^{\mathrm{?}}{}_{\mathrm{<mtext>(</mtext><mtext>ads</mtext><mtext>)</mtext>}}\text{→}\text{H}{}_2\text{O}\text{+}\text{1}\text{2}\text{O}{}_{\mathrm{2(g)}}\text{+ 2e}{}^{\mathrm{?}}{}_{\mathrm{<mtext>crystal</mtext>}}$. This is an important source of atmospheric O₂, amounting to 7 × 10¹³?2 × 10¹⁷ O₂ molecules cm^?2 yr^?1, and it could have a significant effect on atmospheric oxidation state.     

The effect of the internal structure of Mars on its seasonal loading deformations

Mars is continuously subjected to surface loading induced by seasonal mass changes in the atmosphere and ice caps due to the CO₂ sublimation and condensation process. It results in surface deformations and in time variations of gravity. Large wavelength annual and semi-annual variations of gravity (particularly zonal coefficients ΔJn) have been determined using present day geodetic satellite measurements. However loading deformations have been poorly studied for a planet like Mars. In this paper, we compute these deformations and their effect on spacecraft orbiting around Mars. Loading deformations of terrestrial planet are typically investigated assuming a spherical planet, radially symmetric. The mean radial structure of Mars is not well known. In particular the radius of the liquid or solid core remains not precisely determined. One may then wonder what is the effect of these uncertainties on loading deformations. Moreover, Mars presents a strong topography and probably large lateral variations of crustal thickness (relative to the Earth). The paper answer the questions of what is the effect of such lateral heterogeneities on surface deformations, and is the classical way to calculate loading deformation well adapted for a planet like Mars. In order to answer these questions we have investigated theoretically loading deformations of Mars-like planets. We first investigated classical load Love numbers. We show that for degrees inferior to 10, the load Love numbers mainly depend on the radius of the core and on its state, and that for degree greater than 10, they depend on the mean radius of mantle-crust interface. Using a General Circulation Model (GCM) of atmosphere and ice caps dynamics we show that loading vertical displacements have a 4-5 cm magnitude and present a North-South pattern with periodic transitions. Finally we investigated the effect of lateral variations of the crustal thickness on these loading deformations. We show that thickness heterogeneities perturb the deformations and the time variation of gravity at about 0.5%. However this perturbation on ΔJn is only about 1‰ due to main direct attraction of surface fluid layers. We conclude that lateral variations of crustal thickness are today negligible. However, observation of load Love numbers would bring information on the radial internal structure of the planet, particularly on the core radius. ΔJn study would permit to infer the load Love number , particularly for degree 2 and 3, knowing surface fluid layer dynamics. However load Love numbers are quite small (about 0.05), and despite the present good agreement between GCM and ΔJn observations, will only be estimated in the near future when a slightly better precision in observation and modeling will make it possible to infer these numbers. The investigation of load Love number , which are larger than numbers, would be particularly interesting. It would permit to study degree 1 contribution of atmosphere and ice caps dynamics, which is the most important component of surface fluid dynamics on Mars. Surface displacement measurements would be necessary on a few places near the pole regions, which may be possible in the future, with a project involving precise positioning of a lander on the surface of Mars.