Martian Dichotomy

The Martian dichotomy (also known as the global crustal dichotomy or hemispheric dichotomy of Mars) is the most prominent global-scale geomorphological, topographic, and geophysical feature of Mars, dividing the planet into two starkly contrasting hemispheres along an inclined great circle (~30° to the Martian equator): the ancient, heavily cratered, topographically elevated southern highlands (~60% of the planetary surface) and the younger, smooth, low-lying northern plains (Vastitas Borealis, ~40% of the surface). First recognized in global imagery returned by NASA's Mariner 9 orbiter in 1971–1972 and the Viking orbiters in 1976, and subsequently mapped in three dimensions by the Mars Global Surveyor Mars Orbiter Laser Altimeter (MOLA), Mars Reconnaissance Orbiter radar, and the InSight lander's broadband seismometer, the dichotomy manifests across three coupled planetary observables: a 4 to 8 km step in surface elevation, a ~25 to 30 km step in crustal thickness, and a sharp contrast in impact crater density and crustal remanent magnetism.

Dated stratigraphically to the Pre-Noachian epoch (~4.45 to 4.10 billion years ago, within the first 100 to 400 million years of Solar System history), the Martian dichotomy is the oldest preserved megastructure on any terrestrial planet. Within astrobiology and planetary hydrology, the topographic basin of the northern lowlands served as the terminal gravitational sink for all major Noachian valley networks and Hesperian catastrophic outflow channels, providing the basin geometry for a putative ancient northern ocean (Oceanus Borealis) and subsurface chaos terrains.

Precision laser altimetry from MOLA, spherical harmonic gravity fields from MRO and Mars Odyssey, and receiver-function seismic wave conversions recorded by NASA's InSight (Interior Exploration using Seismic Investigations, Geodesy and Heat Transport) mission at Elysium Planitia (2018–2022) have quantified the interior structure of the dichotomy:

1. Topographic Center-of-Mass / Center-of-Figure Offset (3.3 km)

Whereas the northern lowlands (Vastitas Borealis, Arcadia Planitia, Acidalia Planitia, and Utopia Planitia) lie at an average elevation of −4 to −5 km below the Martian equipotential datum (areoid), the southern highlands (Terra Cimmeria, Terra Sirenum, Noachis Terra, and Arabia Terra) stand at +1 to +4 km above datum. This global hemispherical asymmetry shifts Mars's geometric center of figure (COF) by 3.30 km southward along the polar and dichotomy axis relative to its gravitational center of mass (COM).

2. Bimodal Crustal Thickness Distribution (MOLA + Gravity + InSight Seismology)

By inverting Bouguer gravity anomalies anchored to the absolute crustal thickness measured beneath the InSight lander by P-to-S seismic receiver functions (which detected three intra-crustal discontinuities and the Mohorovičić crust–mantle boundary at 39 ± 8 km depth beneath the dichotomy transition zone of Elysium Planitia), geophysicists established that the Martian crust is strongly bimodal:

• Northern Lowlands Crust: Average thickness of H_north ≈ 30 to 38 km (thinning to 10–15 km beneath the Isidis and Utopia impact basins), covered by a 1-to-3-km veneer of Hesperian basaltic lavas and sedimentary deposits (the Vastitas Borealis Formation) that bury an underlying Pre-Noachian/Noachian basement loaded with radar- and gravity-detected "stealth" buried craters.
• Southern Highlands Crust: Average thickness of H_south ≈ 58 to 72 km, Airy-isostatically supporting the 4–6 km surface elevation contrast across a low-density basaltic/andesitic crust (ρ_crust ≈ 2,850–3,000 kg m⁻³ over a mantle of ρ_mantle ≈ 3,500 kg m⁻³).

3. Crustal Remanent Magnetism Contrast

Vector magnetometer measurements by Mars Global Surveyor and MAVEN revealed that the ancient southern highlands (particularly Terra Cimmeria and Terra Sirenum) possess intense crustal remanent magnetization (up to 1,500 nT at 100 km altitude, 10–20 times stronger than terrestrial continental anomalies), whereas the northern lowlands (and the giant Hellas, Argyre, and Utopia basins) are largely unmagnetized. This indicates either that the northern crust formed or was thermally demagnetized above the Curie temperature of magnetite (580 °C) and pyrrhotite (320 °C) after the Martian core dynamo shut down at ~4.1–4.0 Ga, or that the giant event that excavated the northern basin stripped away the magnetized crust before or after dynamo cessation.

The boundary separating the two hemispheres is sharp and steep along ~60% of its circumference (across Deuteronilus Mensae, Protonilus Mensae, Nilosyrtis Mensae, and Nepenthes Mensae) and obscured along the remaining ~40% by the younger Tharsis and Elysium volcanic rises (which host giant shield volcanoes, low-relief pyroclastic calderas or paterae, lobate flood-lava aprons or fluctūs, wrinkle-ridged plateaus or plana, and planetary lava tubes). Along the exposed escarpment separating the southern uplands from northern lowland basins (planitiae) and circumpolar sand seas (undae), Thomas A. Mutch and Robert P. Sharp (1973) defined the classic Martian fretted terrain:

• Flat-Topped Mesas (Mensae) and Knobby Outliers (Colles): Steep-walled, flat-topped remnant plateaus (mensae, 1 to 3 km high) and degraded clusters of conical hills (colles, such as Deuteronilus and Cydonia Colles) detached from the highland plateau by anastomosing flat-floored valleys widened by sapping, mass wasting, and periglacial ice creep—contrasting with the tectonic fracture networks (lineae), dark diapiric chaos spots (maculae), and ancient cratered albedo provinces (regiones) observed on outer Solar System ocean worlds.
• Lobate Debris Aprons (LDAs) and Lineated Valley Fill (LVF): Convex aprons extending 10 to 30 km outward from the base of dichotomy mesas. Subsurface sounding at 20 MHz by the SHARAD radar on MRO proved that these aprons are not dry talus piles, but debris-covered glaciers of nearly pure water ice (H₂O > 80%–90% by volume, permittivity ε′ ≈ 3.15) hundreds of meters thick, preserved beneath a protective meter-thick sublimation lag of rocky debris since prior high-obliquity glacial epochs.

Explaining how Mars acquired a hemispherical crustal thickness step during the Pre-Noachian eon (>4.1 Ga) centers on two major geophysical models (alongside hybrid impact-triggered mantle plume models):

1. The Exogenic Model: The Giant Borealis Basin Oblique Megaimpact

First proposed by Don E. Wilhelms and Steven W. Squyres (1984) for a single impact (and by Frey and Schultz for multiple overlapping basins), and revitalized in 2008 through three-dimensional hydrocode simulations and elliptical boundary reconstructions by Jeffrey C. Andrews-Hanna, Maria T. Zuber, Francis Nimmo, and Margarita M. Marinova:

• Once the younger volcanic load of the Tharsis rise is mathematically subtracted using gravity-topography admittance, the dichotomy boundary traces a single, continuous elliptical basin measuring 8,500 km by 10,600 km, centered at 67° N, 208° E (the Borealis Basin—the largest impact structure in the Solar System, four times wider than the Moon's South Pole–Aitken basin or Mars's Hellas basin).
• Three-dimensional smoothed-particle hydrodynamics (SPH) simulations demonstrate that an oblique impact (impact angle 30° to 60°) by a protoplanetary Differentiated bolide 1,600 to 2,700 km in diameter traveling at 6 to 10 km s⁻¹ (kinetic energy ~10²⁹ to 10³⁰ J) excavates the northern crust to a depth of 20–30 km and blasts ejecta preferentially onto the southern hemisphere without melting the entire global mantle into a deep magma ocean (which would have erased the crustal thickness step by viscous relaxation).

2. The Endogenic Model: Degree-1 Spherical Harmonic Mantle Convection

Proposed by Donald J. Wise, Matthew Golombek, and Phinney (1979) and modeled dynamically by Shijie Zhong and Maria T. Zuber (2001) and J. H. Roberts and Zhong (2006):

• In a stagnant-lid planet whose mantle viscosity increases with depth across the olivine-to-wadsleyite/ringwoodite phase transition (or above a dense basal magma ocean layer), early mantle convection naturally organizes into a spherical harmonic degree-1 planform—consisting of a single giant cylindrical mantle upwelling plume in one hemisphere and a single antipodal downwelling in the other.
• Under the upwelling-in-the-south model, decompression melting above a superplume thickened the southern highland crust by magmatic underplating and intrusion; conversely, under the upwelling-in-the-north / downwelling-in-the-south model, horizontal mantle traction dragged primordial buoyant crust southward to pile up the thick southern highlands (analogous to downwelling models for tessera plateaus on Venus).

Regardless of its Pre-Noachian tectonic or impact origin, the Martian dichotomy controlled the entire subsequent hydrological and astrobiological evolution of Mars:

• Global Hydraulic Gradient and Clay Synthesis: Throughout the late Noachian and early Hesperian (3.9–3.6 Ga), precipitation and groundwater migrating down the global south-to-north topographic gradient carved branching valley networks across the dichotomy escarpment, depositing phyllosilicate clays (Fe/Mg smectites) and carbonate-bearing delta fans into dichotomy-margin craters such as Jezero Crater (explored by the Perseverance rover), Gale Crater (Curiosity rover), and Mawrth Vallis.
• Catastrophic Outflow Channels and Putative Paleo-Shorelines: During the Hesperian epoch, pressurized subsurface aquifers ruptured across chaos terrains (such as Arsinoes and Hydraotes Chaos), unleashing catastrophic floods down Kasei, Ares, Tiu, and Maja Valles directly into the northern Borealis basin. Timothy J. Parker et al. (1989, 1993) and J. Taylor Perron et al. (2007) identified two concentric circum-basin benches along the dichotomy boundary—the Arabia shoreline and the Deuteronilus shoreline (which becomes equipotential once true polar wander induced by Tharsis loading is removed)—interpreted as the paleo-shorelines of a primordial northern ocean (Oceanus Borealis) that may have covered a third of Mars before freezing into subsurface cryospheric reservoirs and modern high-latitude ice sheets, leaving only transient seasonal slope processes such as recurring slope lineae active at the present-day surface.

1. Andrews-Hanna, J. C., Zuber, M. T., & Banerdt, W. B. (2008). "The Borealis basin and the origin of the martian crustal dichotomy." Nature, 453(7199), 1212–1215. https://doi.org/10.1038/nature07011
2. Marinova, M. M., Aharonson, O., & Asphaug, E. (2008). "Mega-impact formation of the Mars hemispheric dichotomy." Nature, 453(7199), 1216–1219. https://doi.org/10.1038/nature07070
3. Zhong, S., & Zuber, M. T. (2001). "Degree-1 mantle convection and the crustal dichotomy on Mars." Earth and Planetary Science Letters, 189(1–2), 75–84. https://doi.org/10.1016/S0012-821X(01)00345-400345-4)
4. Watters, T. R., McGovern, P. J., & Irwin, R. P., III (2007). "Hemispheres Apart: The Crustal Dichotomy on Mars." Annual Review of Earth and Planetary Sciences, 35(1), 621–652. https://doi.org/10.1146/annurev.earth.35.031306.140220
5. Knapmeyer-Endrun, B., et al. (2021). "Thickness and structure of the martian crust from InSight seismic data." Science, 373(6553), 438–443. https://doi.org/10.1126/science.abf8966