Formation Of Tharsis Research Articles

Recent faulting at the Claritas Rupes scarp on Mars

Endogenic processes have greatly affected the Martian surface, especially concentrating at several volcano-tectonic centers. The formation of Tharsis, a vast volcanic bulge, significantly influenced the western hemisphere of Mars. The associated loading stresses caused the formation of various sets of tectonic structures that might have remained active until today. However, surface evidence for very recent endogenic processes in Tharsis is sparse. Here, using the Context Camera (CTX) and High Resolution Imaging Science Experiment (HiRISE) images, we report the presence of fresh-appearing systems of local-scale scarps mostly developed at the southern edge of Tharsis, specifically at Claritas Rupes scarp. These scarps are spatially associated with depression centers situated at the base of Claritas Rupes, inside the Thaumasia Graben which is partially filled by volcanic deposits. The relationships between the studied scarps and present-day surficial processes and its deposits such as rockfalls indicate a young age of the scarps' development, namely, in the range of <1 Ma. The pristine topography (sharp and thin ridges), spatial distribution (association with depression centers), and regional geological context of Claritas Fossae lead us to interpret these scarps as surficial expressions of tectonic activity attributed to normal faulting. This could be related to Deep-seated Gravitational Slope Deformations (DGSDs) released by seismic activity related to ongoing subsidence of depression centers and/or reactivation of the listric normal Claritas Rupes fault. These observations imply that this region experienced long-lasting and multiple volcano-tectonic events and the formation of the youngest deformations could represent neotectonic activity, which has been active until recently, and possibly might be still ongoing.

The history of global strain and geodynamics on Mars

The tectonic record of Mars is dominated by compressional wrinkle ridges on volcanic surfaces, and these structures have been widely used as a record of tectonic and geodynamic evolution. This study analyzes the density of compressional tectonic structures and inferred strain as functions of time, using the lengths and heights of compressional ridges together with geologic estimates of surface age. Our analyses confirm an apparent peak in compressional strain in the late Noachian and early Hesperian, and comparatively lower values before and after. The lower tectonic strain in the early and middle Noachian relative to the early Hesperian reflects the incompleteness of the ancient compressional tectonic record of strain, as older surfaces should accumulate more strain than younger surfaces. This strain deficit necessitates other means of accommodating contractional strain in the ancient crust, including distributed strain and small-scale faulting. The abrupt decrease in the accumulated tectonic strain and strain rate after the early Hesperian reflects a rapid episode of global contraction followed by much lower rates, in conflict with models of steady-state mantle evolution. The decreasing strain rate may be explained by a changing mantle rheology due to volcanic outgassing. Alternatively, the strain history may be explained by the dominance of mantle plumes in the late Noachian and early Hesperian associated with the formation of Tharsis and the early Hesperian volcanic provinces, which could have led to a pulse of rapid contraction caused by disequilibrium cooling and volcanic outpourings. The tectonic record of strain on Mars – including the incomplete ancient record and evidence for strong departures from steady-state models – may have implications for the interpretation of the tectonic record and thermal evolution of other bodies as well.

Unraveling the geologic and tectonic history of the Memnonia-Sirenum region of Mars: Implications on the early formation of the Tharsis rise

The geologic and tectonic history of the Memnonia-Sirenum region of Mars has been reconstructed on the basis of geologic mapping and relative-age dating of rock units and structures using crater counts and superposition relations. Comparisons between the ages of rock unit formation and the orientations of the earliest structures associated with Tharsis demonstrate a strong interplay among tectonic, magmatic, and aqueous activity in the Memnonia-Sirenum region. Early Noachian geologic activity identified within this region includes the formation of large northeast-trending, extensional basins. Stratigraphic and tectonic evidence suggests that these basins formed concurrently with the early formation of Tharsis when the dynamo and associated magnetosphere were still in operation. In addition to the ancient basins, later Tharsis-influenced activity includes prominent tectonic fault systems associated with the formation of Tharsis and Arsia Mons. Detailed mapping of the faults indicates that structures associated with Memnonia and Sirenum Fossae resulted from the formation of Tharsis during the Middle to Late Noachian. Late-stage (Amazonian) activity also identified within this region includes valley networks and collapse depressions that are linked and/or occur among Tharsis-centered faults that cut features throughout this region, including Mangala Valles.

A hybrid origin of the Martian crustal dichotomy: Degree-1 convection antipodal to a giant impact

The Martian crustal dichotomy is the stark ∼5 km difference in surface elevation and ∼26 km difference in crustal thickness between the northern lowlands and southern highlands that originated within 100s of Myr of Mars' formation. The origin of the dichotomy has broad implications for the geodynamic history of Mars, but purely exogenic or endogenic theories so far cannot explain all of the large scale geophysical observations associated with dichotomy formation. A giant impact can produce the shape and slope of the dichotomy boundary, but struggles to explain Mars' remanent crustal magnetic signatures and the ultimate formation of Tharsis. Degree-1 mantle convection can relate the crustal dichotomy to the formation of Tharsis, but does not explain the elliptical dichotomy shape and must be initiated by a large pre-existing viscosity jump in the mantle. We propose a hybrid model of dichotomy formation in which a giant impact induces degree-1 convection with an upwelling antipodal to the impact site. In this scenario, a giant impact in the northern hemisphere excavates crust, creating an initial difference in crustal thickness and possibly composition between the two hemispheres. Over 10s to 100s of Myr, the dominant upwelling(s) would migrate to be under the thicker, insulating crust in the southern hemisphere, generating melt that further thickens the southern crust. We examine this process using 3-D mantle convection simulations, and find that a hemispherical difference in crustal thickness and composition caused by a giant impact can induce degree-1 convection with the upwelling(s) antipodal to the impact site in <100 Myr.

Timing of oceans on Mars from shoreline deformation.

Widespread evidence points to the existence of an ancient Martian ocean. Most compelling are the putative ancient shorelines in the northern plains. However, these shorelines fail to follow an equipotential surface, and this has been used to challenge the notion that they formed via an early ocean and hence to question the existence of such an ocean. The shorelines' deviation from a constant elevation can be explained by true polar wander occurring after the formation of Tharsis, a volcanic province that dominates the gravity and topography of Mars. However, surface loading from the oceans can drive polar wander only if Tharsis formed far from the equator, and most evidence indicates that Tharsis formed near the equator, meaning that there is no current explanation for the shorelines' deviation from an equipotential that is consistent with our geophysical understanding of Mars. Here we show that variations in shoreline topography can be explained by deformation caused by the emplacement of Tharsis. We find that the shorelines must have formed before and during the emplacement of Tharsis, instead of afterwards, as previously assumed. Our results imply that oceans on Mars formed early, concurrent with the valley networks, and point to a close relationship between the evolution of oceans on Mars and the initiation and decline of Tharsis volcanism, with broad implications for the geology, hydrological cycle and climate of early Mars.

The formation of the South Tharsis Ridge Belt: Basin and Range‐style extension on early Mars?

AbstractThe South Tharsis Ridge Belt (STRB) is located along the northeastern edge of Terra Sirenum and partially surrounds the southwestern part of Tharsis in an arc. It consists of 29 large ridges separated by distances 130 to 260 km, with average relief of 1.5 km above the surrounding plains. Because the STRB is among the oldest tectonic features associated both spatially and developmentally with Tharsis, it may provide key information on the early evolution of Tharsis and possibly pre‐Tharsis processes. Earlier studies concluded that the ridges formed through compressional tectonism. However, the shape, size, and separation of the ridges support the interpretation that the STRB resembles the extensional Basin and Range Province on Earth. Both regions are characterized by series of parallel mountain ranges separated by broad valleys. In this study, we evaluate both extensional and compressional hypotheses for the origin of the ridges using evidence from topography, deformed craters, crustal thickness models, and strain modeling. Though neither interpretation explains all aspects of the ridges, the topography of the ridges and crustal thinning associated with the western part of the ridge belt support an extensional origin. Strain models predict that Basin and Range‐style wide rifting would be expected for early Martian conditions. This extension may have been initiated by plume‐induced uplift in the early stages of Tharsis formation, but the large amount of extensional strain inferred in the western STRB must have been accommodated by compression elsewhere, possibly in the heavily deformed craters of western Terra Sirenum.

Late Tharsis formation and implications for early Mars.

The Tharsis region is the largest volcanic complex on Mars and in the Solar System. Young lava flows cover its surface (from the Amazonian period, less than 3 billion years ago) but its growth started during the Noachian era (more than 3.7 billion years ago). Its position has induced a reorientation of the planet with respect to its spin axis (true polar wander, TPW), which is responsible for the present equatorial position of the volcanic province. It has been suggested that the Tharsis load on the lithosphere influenced the orientation of the Noachian/Early Hesperian (more than 3.5 billion years ago) valley networks and therefore that most of the topography of Tharsis was completed before fluvial incision. Here we calculate the rotational figure of Mars (that is, its equilibrium shape) and its surface topography before Tharsis formed, when the spin axis of the planet was controlled by the difference in elevation between the northern and southern hemispheres (hemispheric dichotomy). We show that the observed directions of valley networks are also consistent with topographic gradients in this configuration and thus do not require the presence of the Tharsis load. Furthermore, the distribution of the valleys along a small circle tilted with respect to the equator is found to correspond to a southern-hemisphere latitudinal band in the pre-TPW geographical frame. Preferential accumulation of ice or water in a south tropical band is predicted by climate model simulations of early Mars applied to the pre-TPW topography. A late growth of Tharsis, contemporaneous with valley incision, has several implications for the early geological history of Mars, including the existence of glacial environments near the locations of the pre-TPW poles of rotation, and a possible link between volcanic outgassing from Tharsis and the stability of liquid water at the surface of Mars.

Alignments of volcanic features in the southern hemisphere of Mars produced by migrating mantle plumes

Mars shows alignments of volcanic landforms in its southern hemisphere, starting from the equatorial regions and converging towards the South Pole, and visible at global scale. These composite alignments of volcanoes, calderas, shields, vents, heads of valley networks and massifs between the equatorial regions and the southern polar region define twelve different lines, fitted by rhumb lines (loxodromes), that I propose to be the traces of mantle plumes. The morphology of the volcanic centres changes along some of the alignments suggesting different processes of magma emplacement and eruptive style. The diameters of the volcanic centres and of the volcanic provinces are largest at Tharsis and Elysium, directly proportional to the number of alignments starting from them. A minor presence of unaligned volcanic features is observed on the northern lowlands and on the highlands outside the 12 major alignments. The heads of channels commonly interpreted as fluvial valleys are aligned with the other volcanic centres; unaltered olivine is present along their bed-floors, raising severe doubts as to their aqueous origin. Several hypotheses have tried to explain the formation of Tharsis with the migration of a single mantle plume under the Martian lithosphere, but the discovery of twelve alignments, six starting from Tharsis, favours the hypothesis of several mantle plumes as predicted by the model of the Southern Polar Giant Impact (SPGI) and provides a new view on the formation of the volcanic provinces of Mars.

Martian crustal dichotomy and Tharsis formation by partial melting coupled to early plume migration

A recently proposed model links the formation and early evolution of the Tharsis volcanic province on Mars to the preexisting hemispheric dichotomy (Zhong, 2009). A key aspect of this model is the assumption of a deep lithospheric root below the thicker crust of the southern highlands. We implemented a parameterization of partial melting into the 3‐D spherical shell mantle convection code CitcomS in order to investigate whether the required lithospheric thickness variation can be generated self‐consistently by partial melting when stiffening of the melt residue due to devolatilization is considered. The rate of melt production strongly depends on the mantle temperature, and additional strong coupling between the flow and partial melting is introduced through the stiffening effect on the melt residue. We find that it is possible to generate a lithospheric keel by partial melting above a single upwelling that excites a relative rotation between the one‐plate lithosphere and the mantle below while producing the amount of melt distributed in a broad region constrained to one hemisphere that is necessary to form the crustal dichotomy. This scenario thus offers an internal mechanism for the Martian dichotomy formation and validates the hypothesis of Zhong (2009).

CO 2 solubility in Martian basalts and Martian atmospheric evolution

To understand possible volcanogenic fluxes of CO 2 to the Martian atmosphere, we investigated experimentally carbonate solubility in a synthetic melt based on the Adirondack-class Humphrey basalt at 1–2.5 GPa and 1400–1625 °C. Starting materials included both oxidized and reduced compositions, allowing a test of the effect of iron oxidation state on CO 2 solubility. CO 2 contents in experimental glasses were determined using Fourier transform infrared spectroscopy (FTIR) and Fe 3+/Fe T was measured by Mössbauer spectroscopy. The CO 2 contents of glasses show no dependence on Fe 3+/Fe T and range from 0.34 to 2.12 wt.%. For Humphrey basalt, analysis of glasses with gravimetrically-determined CO 2 contents allowed calibration of an integrated molar absorptivity of 81,500 ± 1500 L mol −1 cm −2 for the integrated area under the carbonate doublet at 1430 and 1520 cm −1. The experimentally determined CO 2 solubilities allow calibration of the thermodynamic parameters governing dissolution of CO 2 vapor as carbonate in silicate melt, K II , ( Stolper and Holloway, 1988) as follows: ln K II 0 = - 15.42 ± 0.20 , Δ V 0 = 20.85 ± 0.91 cm 3 mol −1, and Δ H 0 = −17.96 ± 10.2 kJ mol −1. This relation, combined with the known thermodynamics of graphite oxidation, facilitates calculation of the CO 2 dissolved in magmas derived from graphite-saturated Martian basalt source regions as a function of P, T, and f O 2 . For the source region for Humphrey, constrained by phase equilibria to be near 1350 °C and 1.2 GPa, the resulting CO 2 contents are 51 ppm at the iron–wüstite buffer (IW), and 510 ppm at one order of magnitude above IW (IW + 1). However, solubilities are expected to be greater for depolymerized partial melts similar to primitive shergottite Yamato 980459 (Y 980459). This, combined with hotter source temperatures (1540 °C and 1.2 GPa) could allow hot plume-like magmas similar to Y 980459 to dissolve 240 ppm CO 2 at IW and 0.24 wt.% of CO 2 at IW + 1. For expected magmatic fluxes over the last 4.5 Ga of Martian history, magmas similar to Humphrey would only produce 0.03 and 0.26 bars from sources at IW and IW + 1, respectively. On the other hand, more primitive magmas like Y 980459 could plausibly produce 0.12 and 1.2 bars at IW and IW + 1, respectively. Thus, if typical Martian volcanic activity was reduced and the melting conditions cool, then degassing of CO 2 to the atmosphere may not be sufficient to create greenhouse conditions required by observations of liquid surface water. However, if a significant fraction of Martian magmas derive from hot and primitive sources, as may have been true during the formation of Tharsis in the late Noachian, that are also slightly oxidized (IW + 1.2), then significant contribution of volcanogenic CO 2 to an early Martian greenhouse is plausible.

Mars without the equilibrium rotational figure, Tharsis, and the remnant rotational figure

We use a revised partitioning of the planet figure into equilibrium and nonequilibrium contributions that takes into account the presence of an elastic lithosphere to study the Martian gravity field and shape. The equilibrium contribution is associated with the present rotational figure, and the nonequilibrium contribution is dominated by Tharsis and a remnant rotational figure supported by the elastic lithosphere that traces the paleopole location prior to the formation of Tharsis. We calculate the probability density functions for Tharsis' size and location, the paleopole location, and the global average thickness of the elastic lithosphere at the time Tharsis was emplaced. Given the observed degree‐3 spherical harmonic gravity coefficients, the expected Tharsis center location is 258.6 ± 4.2°E, 9.8 ± 0.9°N, where the uncertainties represent the 90% confidence interval. Given this Tharsis center location and the observed degree‐2 spherical harmonic gravity coefficients, the expected paleopole location prior to the emplacement of Tharsis is 259.5 ± 49.5°E, 71.1−14.4+17.5°N, and the expected elastic lithospheric thickness at the time of loading is 58−32+34km. Our estimated paleopole colatitude implies 18.9−17.5+14.4°of true polar wander (TPW) driven by the emplacement of Tharsis, in disagreement with previous studies that invoke large TPW. The remnant rotational figure is visible in both the nonequilibrium degree‐2 geoid (areoid) without Tharsis and the nonequilibrium degree‐2 topography without Tharsis. The remnant rotational figure is also visible in the total nonequilibrium geoid without Tharsis, but it is not visible in the total nonequilibrium topography without Tharsis due to the strong signal of the north‐south dichotomy. Shorter wavelength geological features become significantly more visible in the geoid with the removal of the long wavelength contributions of the equilibrium rotational figure, Tharsis, and the remnant rotational figure. Removal of the equilibrium rotational figure and Tharsis from the topography reveals a better defined north‐south dichotomy boundary.

Sulfuric acid aerosols in the atmospheres of the terrestrial planets

Clouds and hazes composed of sulfuric acid are observed to exist or postulated to have once existed on each of the terrestrial planets with atmospheres in our solar system. Venus today maintains a global cover of clouds composed of a sulfuric acid/water solution that extends in altitude from roughly 50 km to roughly 80 km. Terrestrial polar stratospheric clouds (PSCs) form on stratospheric sulfuric acid aerosols, and both PSCs and stratospheric aerosols play a critical role in the formation of the ozone hole. Stratospheric aerosols can modify the climate when they are enhanced following volcanic eruptions, and are a current focus for geoengineering studies. Rain is made more acidic by sulfuric acid originating from sulfur dioxide generated by industry on Earth. Analysis of the sulfur content of Martian rocks has led to the hypothesis that an early Martian atmosphere, rich in SO 2 and H 2O, could support a sulfur-infused hydrological cycle. Here we consider the plausibility of frozen sulfuric acid in the upper clouds of Venus, which could lead to lightning generation, with implications for observations by the European Space Agency's Venus Express and the Japan Aerospace Exploration Agency's Venus Climate Orbiter (also known as Akatsuki). We also present simulations of a sulfur-rich early Martian atmosphere. We find that about 40 cm/yr of precipitation having a pH of about 2.0 could fall in an early Martian atmosphere, assuming a surface temperature of 273 K, and SO 2 generation rates consistent with the formation of Tharsis. This modeled acid rain is a powerful sink for SO 2, quickly removing it and preventing it from having a significant greenhouse effect.

The formation of Tharsis on Mars: What the line‐of‐sight gravity is telling us

Line‐of‐sight (LOS) spacecraft acceleration profiles from the Radio Science Experiment and topography from the Mars Orbiter Laser Altimeter (MOLA) instrument of the Mars Global Surveyor (MGS) are analyzed to estimate the effective elastic thickness (Te) for various regions of Tharsis. We identify a buried basin flanking the Thaumasia Highlands at the southeastern margin of Tharsis. Assuming that this basin results from lithospheric flexure from surface loading by the Thaumasia Highlands, we fit LOS profiles across the feature with a thin‐shell, elastic flexure model and find the mountain belt to reflect a value of Te ∼ 20 km consistent with a Noachian formation age. We also determine admittances from LOS profiles for five regions across Tharsis and fit them with theoretical admittances calculated using the flexural model. Crater density, surface density, and predominant surface age are found to vary systematically across Tharsis while Te does not. The highest surface density and lowest Te values are obtained for the western portion of Tharsis where crater densities are lowest. Our results imply the majority of the topographic rise was emplaced within the Noachian irrespective of the surface ages. Topographic loading and resurfacing (i.e., volcanic activity) persisted into the Amazonian while becoming increasingly confined to the western margin where the youngest surface ages are found and the eruptive style transitioned from effusive volcanism to shield‐forming volcanism as Te increased.

Strike‐slip faults on Mars: Observations and implications for global tectonics and geodynamics

The tectonic evolution of Mars has been driven primarily by the interaction of Tharsis‐induced loading stresses with a uniform contractional stress field, leading to global assemblages of graben and wrinkle ridges. Until recently, strike‐slip faults have appeared to be largely absent from the tectonic record. We here present evidence for a new set of Noachian to Early Hesperian strike‐slip faults southwest of Tharsis. The best‐preserved fault in the newly mapped set has a length of ∼200 km and measured lateral offset of 5–9 km. We use the locations of both previously mapped and newly identified strike‐slip faults to investigate the tectonic and geodynamic evolution of Mars. We model the tectonic evolution of the planet in response to the evolving balance between the loading and contractional stresses, finding widespread regions in which strike‐slip faulting was favored. The observed tectonic history is consistent with a scenario in which loading‐induced stresses peaked early in Mars history, followed by the growth of contractional stresses, leading to a shift in the preferred locus of strike‐slip faulting from regions southwest of Tharsis in the Noachian to Early Hesperian, to northwest of Tharsis in the Early Amazonian. The contractional strain history inferred from the ages and locations of the strike‐slip faults is consistent with the strains calculated from models of the thermal evolution of Mars in which the gradual secular cooling of the planet over time was augmented by plume‐induced contraction during Tharsis formation in the Noachian.

Equilibrium rotational stability and figure of Mars

Studies extending over three decades have concluded that the current orientation of the martian rotation pole is unstable. Specifically, the gravitational figure of the planet, after correction for a hydrostatic form, has been interpreted to indicate that the rotation pole should move easily between the present position and a site on the current equator, 90° from the location of the massive Tharsis volcanic province. We demonstrate, using general physical arguments supported by a fluid Love number analysis, that the so-called non-hydrostatic theory is an inaccurate framework for analyzing the rotational stability of planets, such as Mars, that are characterized by long-term elastic strength within the lithosphere. In this case, the appropriate correction to the gravitational figure is the equilibrium rotating form achieved when the elastic lithospheric shell (of some thickness LT) is accounted for. Moreover, the current rotation vector of Mars is shown to be stable when the correct non-equilibrium theory is adopted using values consistent with recent, independent estimates of LT. Finally, we compare observational constraints on the figure of Mars with non-equilibrium predictions based on a large suite of possible Tharsis-driven true polar wander (TPW) scenarios. We conclude, in contrast to recent comparisons of this type based on a non-hydrostatic theory, that the reorientation of the pole associated with the development of Tharsis was likely less than 15° and that the thickness of the elastic lithosphere at the time of Tharsis formation was at least ∼ 50 km . Larger Tharsis-driven TPW is possible if the present-day gravitational form of the planet at degree 2 has significant contributions from non-Tharsis loads; in this case, the most plausible source would be internal heterogeneities linked to convection.

Evolution of the Tharsis region of Mars: insights from magnetic field observations

Mars Global Surveyor (MGS) observations of crustal magnetic fields over Tharsis provide new constraints on models for the thermal and magmatic evolution of this region. We analyze the distribution of magnetic field anomalies over Tharsis surface units of Noachian, Hesperian and Amazonian age. These data suggest that early Noachian crust underlies the Tharsis province, and formed contemporaneously with the existence of a martian dynamo. This crust either pre-dates the formation of Tharsis, or formed during the earlier phases of Tharsis volcanism. The preservation of strong magnetic field anomalies over some of the earliest Noachian and topographically high units, together with the observation of magnetic field anomalies over Hesperian- and Amazonian-age surface units, indicate that a large fraction of the magnetized crust has remained cool (below the blocking temperature of the magnetic carrier) throughout the construction of Tharsis. Moreover, the distributions of magnetic anomaly amplitudes over Noachian, Hesperian, and Amazonian surface units suggest that the youngest units overlie sites of prolonged intrusion and have undergone a greater extent of thermal demagnetization. The absence of magnetic anomalies around the Tharsis Montes and Olympus Mons argues for strong, localized heating, as would be expected at volcanic centers. We show that end-member models for progressive thermal demagnetization of a Noachian magnetized crustal layer are consistent with the anomaly amplitude distributions. We integrate the magnetic field observations with constraints from tectonics, gravity, and topography, and present a revised scenario for the evolution of the Tharsis region.

Lineament Analysis and Geophysical Modelling of the Alba Patera Region on Mars

Global data sets of images, topography and gravity are available for Mars from several orbiter missions. At the eve of new global data from Mars Global Surveyor (MGS), the capabilities of 3D geophysical modelling based on areal topography and gravity data combined with geologic-tectonic image interpretation is demonstrated here. A unique structure is chosen for the model calculations: the Alba Patera volcanic complex at the northern border of the Tharsis rise. Five groups of graben are discriminated: Ceraunius Fossae, Catenae, Tantalus Fossae (radial group) radial to the Tharsis rise, mainly associated to the formation of Tharsis, and Alba and Tantalus Fossae (circular group), younger than the other graben and circular around Alba Patera. Combining 3D elastic flexure of the lithosphere due to a 3D topographic surface load with 3D gravity models results in a rather thick lithosphere (150–200 km) and thick crust (60–100 km). In another model estimate it has been assumed that the circular grabens are induced by the stresses from the surface load of Alba Patera. In a first order calculation the surface stresses under a point load have been determined resulting in a good correlation of the stress maximum with the location of the circular grabens for a 50-km thick lithosphere. This is in accordance with earlier results from this method, but in contradiction with the thick lithosphere derived from flexure-gravity models. One possibility for this contradiction may be that the different models represent two evolutionary points of Alba Patera. (1) The correlation of stresses with the circular grabens may represent an older stage of evolution with a thinner lithosphere. (2) The flexure-gravity models represent a younger to present stage with a thick lithosphere. The results of the lithosphere thicknesses are compared with an admittance calculation and different thermal evolution models which determine comparable thicknesses (≈150 km). More detailed models including 3D stress models should wait for new data sets from MGS. The results from the lineament analysis and geophysical modelling are summarized in an evolution model for Alba Patera.

Shape of the northern hemisphere of Mars from the Mars Orbiter Laser Altimeter (MOLA)

Eighteen profiles of ∼N‐S‐trending topography from the Mars Orbiter Laser Altimeter (MOLA) are used to analyze the shape of Mars' northern hemisphere. MOLA observations show smaller northern hemisphere flattening than previously thought. The hypsometric distribution is narrowly peaked with &gt;20% of the surface lying within 200 m of the mean elevation. Low elevation correlates with low surface roughness, but the elevation and roughness may reflect different mechanisms. Bouguer gravity indicates less variability in crustal thickness and/or lateral density structure than previously expected. The 3.1‐km offset between centers of mass and figure along the polar axis results in a pole‐to‐equator slope at all longitudes. The N‐S slope distribution also shows a subtle longitude‐dependent variation that may represent the antipodal effect of the formation of Tharsis.

Mars without Tharsis

The significant power in the Martian gravitational field due to the Tharsis rise may mask or modify gravitational signatures that contain important information on Martian geophysical processes. In order to isolate gravity signals in regions where the field is significantly affected by Tharsis as well to investigate characteristics of the global field, we present a numerical technique to “remove” Tharsis from the Martian gravitational field. Our analysis will show that to first order Tharsis can be represented as a sixth‐degree spherical harmonic zonal gravitational feature in a reference frame in which the axis of symmetry runs through the center of the province. We produced the “Mars without Tharsis” (MWT) field by subtracting the gravitational signature of Tharsis from the full field and rotating the spherical harmonics back to the original coordinate system, in which the z‐axis is coincident with the planetary rotation axis. Our study yields limits on the moment of inertia factor of 0.361&lt;C/MR2&lt;0.366 for the range of allowable non‐hydrostatic contributions to the flattening. The analysis indicates that the formation of Tharsis may have resulted in re‐orientation of the spin axis, though if the planet retained rotational fossil bulge due to a thick elastic lithosphere then changes in the global figure could have been inhibited. The zonal and tesseral components of the degree 2 field in the GMM‐I and rotated configurations indicates that Tharsis is ∼5 times more axisymmetric about its “pole” than Mars as a whole is about the rotation axis. With the removal of Tharsis' long wavelength gravitational power, regional structures become better defined in a spatial sense while maintaining their local dynamic range. The Olympus Mons volcanic shield displays clear evidence of a gravitational flexural moat. Our method for removing Tharsis may also be applied to analyze other global‐scale geophysical features that may be considered axisymmetric to first order, such as mantle plumes, major volcanic constructs, or large impact basins.

Stress models for Tharsis formation, Mars

A critical survey is presented of most stress models proposed for the formation of the tectonic structures in the Tharsis volcano-tectonic province on Mars and provides new constraints for further models. First papers, in the 1970s, attempted to relate the Tharsis formation to asthenospheric movements and lithosphere loading by magma bodies. These processes were then quantified in terms of stress trajectory and magnitude models in elastic lithosphere (e.g. Banerdt et al., J. Geophys. Res. 87(B12), 9723–9733, 1982). Stresses generated by dynamic lithosphere uplift were rapidly dismissed because of the poor agreement between the stress trajectories provided by the elastic models and the structural observations. The preferred stress models involved lithosphere loading, inducing isostatic compensation, and then lithosphere flexure. Some incomsistency with structural interpretation of Viking imagery has been found. In the early 1990s, an attempt to solve this problem resulted in a model involving the existence of a Tharsis-centred brittle crustal cap, deteched from the strong mantle by a weak crustal layer (Tanaka et al., J. Geophys. Res. 96(E1), 15617–15633, 1991). Such a configuration should produce loading stresses akin to those predicted by some combination of the two loading modes. This model has not been quantified yet, however it is expected to reconcile stress trajectories and most structural patterns. Nevertheless, some inconsistencies with observed structures are also expected to remain. Parallel to this approach focused on loading mechanisms, the idea that volcanism and tectonic structures could be related to mantle circulation began to be considered again through numerical convection experiments, whose results have however not been clearly correlated with surface observations. Structural clues to early Tharsis dynamic uplift are reported. These structures have little to do with those predicted by elastic stress modelling of dynamic lithosphere uplift. They denote the existence of unsteady stress trajectories responsible for surface deformations that cannot be readily predicted by elastic models. These structures illustrate that improving current stress models for Tharsis formation shall come from deeper consideration of rock failure criterion and load growth in the lithosphere (e.g. Schultz and Zuber, J. Geophys. Res. 99(E7), 14691–14702, 1994). Improvements should also arise from better understanding rheological layering in the lithosphere and its evolution with time, and from consideration of stress associated to magma emplacement in the crust, which may have produced many tectonic structures before loading stress resulting from magma freezing became significant (Mège and Masson, Planet. Space Sci. 44, 1499–1546, 1996a).