Main

Before the emergence of a recognizable form of plate tectonics on Earth circa 3 billion years ago (Ga), the Earth was resurfaced by earlier forms of volcanotectonics1. Despite the importance for understanding planetary evolution and the environments where life originated on this planet, there is little consensus on the style and process of crustal resurfacing in the Eoarchean and Paleoarchean (3.2–4.0 Ga)2. Understanding Archaean crustal evolution relies heavily on models and geochemical proxies within rare, highly metamorphosed and metasomatized terranes because most of our planet’s ancient crust has been lost to erosion and younger crustal recycling3.

Geological exploration of other rocky planets provides clues to early crustal evolution and volcanotectonic processes because other bodies exemplify comparable geosystems with variable mass, gravity and composition4. Observed through the lenses of different gravity fields, bulk planet compositions and heat flow, it is possible to test models for crustal resurfacing and so discover the steps that led to plate tectonics and other forms of crustal recycling5. The planet Mars represents a particularly valuable puzzle piece in this regard. Approximately 70% of the Martian surface is >3 billion years old, and 45% is older than 3.6 billion years6. Though the crust is heavily impact cratered, the ancient geologic record remains largely intact and therefore provides an invaluable window into early solar system geological conditions7 and clues to early crustal evolution.

One essential element in the pathway toward plate tectonics is the existence of continents composed of felsic rock8, but it remains unclear how the earliest tronjemite–tonalite–granodiorite (TTG) or granite–monzonite–syenite (GMS) plutons and related dacite–rhyolite volcanic systems formed on Earth9. It is possible that plume-dominated volcanism in an Icelandic-type setting resulted in rhyolitic volcanism through fractional crystallization driven by melting of weathered, hydrated volcanic materials buried in a high geothermal gradient10. On the Archaean Earth, lithospheric delamination resulted in crustal recycling of altered, hydrated crust11,12,13,14, and a similar process has been hypothesized for ancient Mars15,16,17.

Mars has long been considered a basaltic world with little igneous diversity, but discoveries in the past decade of enclaves of evolved plutonic and volcanic rocks18,19,20,21,22,23,24 hint that this paradigm requires reframing. This work explores the geology and volcanology of the Eridania basins and surrounding terrains in Terra Cimmeria and Terra Sirenum (180° E, 30° S). This region contains the strongest signature of remnant crustal magnetism on the planet25, a high potassium anomaly26, large contrasts in crustal thickness27 and volcanically and structurally complex Hesperian–Noachian (>3 Ga) crust6,28,29,30,31 within a unique Martian terrane or crustal block32. In addition, the region contains evidence for an ancient Martian sea approximately three times the size of the modern Caspian Sea on Earth with associated widespread alteration and hydration of the crust33. Though this part of Terra Cimmeria and Terra Sirenum is recognized for having some volcanic vents and structures6,28,31,34,35, the provenance of olivine-rich flood lavas and the widespread (1.8 × 106 km2), voluminous (106 km3) ‘Electris’ suite of airfall deposits28,36 remains enigmatic37. This work describes a diverse suite of volcanic structures and associations of those structures with tectonics and crustal properties. The diverse volcanism is associated with felsic volcanic compositions which are unlike any other suite of recognized deposits or volcanic region on Mars.

Volcanic morphology and morphometrics

The Eridania region of Mars contains at least 63 (identified here) and possibly hundreds of volcanoes, many of which were previously unrecognized (Fig. 1, Supplementary Table 1 and Supplementary Fig. 1). We classified the structures into four categories based on their morphologies and morphometric characteristics: (1) volcanic domes, (2) stratovolcanoes, (3) pyroclastic shields and (4) caldera complexes (Fig. 2). Volcanic domes are characterized by their smaller size (average diameter, d = 14 km), low height (average height, h = 0.8 km), viscous flow features, lobate plan view, sharp edges or escarpments (Supplementary Fig. 1) and high average height/diameter ratio (h/d) of 0.056 (Fig. 3). Stratovolcanoes have some similar features to domes, such as steep escarpments and some lobate flow features elevated on the structure, but overall have a shallow, conical morphology and often contain summit or flank craters, sometimes structurally linked through faults and fractures. They are typically 40–80 km in diameter and are the tallest of the structures described here, with heights reaching 2–3 km above the surrounding terrain with moderate h/d (0.036) (Fig. 3). Pyroclastic shields are large features approximately 100 km in diameter with low heights (1.1 km) and low corresponding h/d (.011). They are quasi-circular in plan view and contain a summit or near-summit crater. A striking and important characteristic of the pyroclastic shields and stratovolcanoes is that they contain dense valley networks eroded into their flanks (Extended Data Fig. 1), which we propose is a reflection of the susceptibility to erosion of their pyroclastic rock. Caldera complexes are unique because they are destructional volcanic features, not constructional ones38, hence their negative average h/d value (−0.015). They form through explosive volcanism and surface collapse but build no central edifice. A critical aspect of their morphologies is the occurrence of multiple structurally linked topographic depressions lacking the characteristic features of impact craters such as raised rims, ejecta and central peaks38. Some of the plains-style caldera complexes in the Terra Cimmeria and Terra Sirenum terrains have as many as 14 structurally linked depressions.

Fig. 1: A topographic map of the Eridania region of Mars.
figure 1

MOLA/HRSC elevation data are draped onto hillshade data. Blue tones are used to show the extent and depth of a large inland sea. Volcanic structures described in this study are classified by morphology and morphometrics.

Fig. 2: Morphology of four classes of volcanoes in the Eridania region.
figure 2

a,b,c,d, Type examples of a volcanic dome (a), stratovolcano (b), pyroclastic shield (c) and caldera complex (d) are shown with MOLA/HRSC elevation data draped over Mars Context Imager data. e, Topographic profiles of each feature represent SW–NE profile lines at each site. Scale bars, 5 km.

Fig. 3: Morphometric data for volcanic structures.
figure 3

Data for volcanic domes, stratovolcanoes and pyroclastic shields in the Eridania region described in this study are compared to similar measurements for structures on the Earth, Moon and Mars. Diameter and height/diameter ratios are shown for rhyolitic–dacitic domes, stratovolcanoes and shield volcanoes on Earth, large shields on Mars and large shields on the Moon.

Morphometrics illustrate important trends regarding Martian volcanology and also shine a light on critical differences between volcanism on the Earth and Mars. The smallest features described here are volcanic domes, and the largest-diameter features are shield volcanoes; stratovolcanoes have similar heights to shields but with higher h/d (Fig. 3). The Martian volcanic structures discovered here follow the same trend as comparable classes of volcanoes on Earth, but they are shifted towards larger diameters, and therefore the Martian features have lower h/d values (Fig. 3). This can be explained in the context of Mars’s lower-gravity environment (38% of Earth), resulting not only in more explosive volcanism but dispersal of pyroclasts further from the vents, resulting in lower profile features39. One key difference on Mars is the occurrence of very large shield volcanoes, comparable in diameter to those on the Moon, but much taller. Those features formed largely through sustained hot-spot volcanism over billions of years, into the Amazonian (<3 Ga) (ref. 40). The older, smaller, more subtle volcanic features described here could be present throughout the ancient crust of Mars41 and could provide a fundamental connection between source and sink for many highland deposits currently lacking provenance42.

Composition of volcanic materials

Thermal infrared remote sensing of the Eridania region demonstrates a diversity of compositions among the volcanic materials. Decorrelation-stretched Thermal Emission Imaging System (THEMIS) images (Fig. 4 and Supplementary Fig. 2) show mafic to ultramafic surfaces as purple-blue and more felsic surfaces as yellow24 (Fig. 4a), linking qualitative compositional variation to morphological context43 (Extended Data Fig. 2). Surface composition is quantitatively constrained by fitting a spline curve to the atmospherically corrected surface emissivity spectra (Fig. 4b), providing a central wavelength of the Si–O stretching absorption in silicates44, which is correlated with SiO2per cent (Fig. 4c). THEMIS spline fits suggest that the felsic volcaniclastic surfaces and volcanic domes have relatively high silica content (Fig. 4b,c and Supplementary Fig. 4). By comparing the THEMIS spline fits to those of rocks measured in the laboratory and adjusted to similar spectral sampling, it is possible to quantitatively constrain the SiO2 content of the Martian volcanic materials. If the Martian rocks are alkaline volcanics, they are likely 64–69% SiO2 and if they are glass-rich, they could be 71–75% SiO2. Given their high potassium content and morphology, it is unlikely they are glass-poor subalkaline rocks, but if they are in this category, the silica content is estimated at 58–59% SiO2 (Fig. 4c). It bears noting that basaltic sand on the felsic surfaces could make the silica content appear lower than it is in the primary rock. High silica content in the same surfaces is supported by high values of the Thermal Emission Spectrometer (TES) silica index, which independently measures Si–O bending absorptions in the region (Extended Data Fig. 3). TES and THEMIS data both suggest little or no quartz and high feldspar/pyroxene ratios. The high-silica volcanic deposits generally correspond to light-toned, buff-coloured, layered or massively bedded, surface draping material in high-resolution images (Fig. 4d,e), but THEMIS data of some central uplifts of impact craters shows that silica-rich material exists to at least 5 km depth (Supplementary Fig. 2).

Fig. 4: Infrared and image data showing felsic volcanic materials in the Eridania region.
figure 4

a, A THEMIS decorrelation-stretched image (I08708002) shows felsic materials as yellow and mafic material as purple surrounding a stratovolcano (labelled ‘SV’) at 170.3° E, 35° S. b,c, Surface emissivity spectra of the felsic (yellow) and mafic (purple) surfaces are compared to global average TES surface type 1 (ST1) and surface type 2 (ST2). Spline fits to the THEMIS data shown in b correspond to 9.36–9.5 μm for the yellow felsic surfaces and 10.74 μm for the purple mafic, and when compared to the plot shown in c of the spline fits for laboratory rocks and glasses, these values suggest high-silica compositions (yellow zone in c) as described in the main text. d,e, HiRISE data show that the felsic material is buff- to pink-coloured massive (d) to finely layered surface deposits (e), which we interpret as volcanic ash (images ESP_025301_1440 and ESP_023745_1465, respectively). Locations of HiRISE infrared-red-blue/green images are indicated with arrows in a.

Near-infrared reflectance spectra from the Compact Reconnaissance Imaging Spectrometer for Mars (CRISM) of subaerially deposited volcanic deposits are generally spectrally unremarkable, with the exception of their evidence for bound water that could occur as zeolites or clay minerals (Supplementary Fig. 3) formed through chemical weathering. By contrast, CRISM spectra of ancient volcanic deposits in the deep basin below interpreted sea level contain abundant and diverse Fe- and Mg-rich clay minerals, carbonates, and other alteration phases33. Mafic volcanics deposited subaerially occur as flood lavas and show strong absorptions related to olivine and pyroxene. CRISM is sensitive to some hydrous igneous minerals such as mica and amphibole, but neither is detected.

Volcanotectonics and crustal recycling

The morphologically and compositionally diverse volcanism observed in Terra Cimmeria and Terra Sirenum points to an episode of magmatic evolution and possible crustal recycling. The volcanoes occur within a zone of potassium-enrichment (Fig. 5a), but the footprint of the Gamma Ray Spectrometer (GRS) data is large (500 km) and it is not possible to trace the K to host rocks from those data. While it would be useful to know the sodium abundance in the region, GRS is not able to measure this element. The K likely occurs within feldspars and volcanic glass, both of which appear relatively abundant in the area from TES and THEMIS data.

Fig. 5: Volcanic structures in the context of regional crustal properties.
figure 5

a, MOLA elevation data show that the Eridania basin is composed of a series of 1- to 3-km-deep subbasins. Volcanoes identified in this study are marked with black squares. b, The volcanoes occur in a zone of high crustal potassium. c,d, The basins (c) have 10 to 20 km thinner crusts. Wrinkle ridges marked by black lines (c,d) occur in NE- and NW-trending lineaments, which also mark zones (d) of disruption of the remnant crustal magnetism. Lat., latitude; Br, radial magnetic component.

The volcanic structures occur within and adjacent to the Eridania basin, which is composed of a WNW-trending suite of five or more large (200 km diameter), deep (1–3 km) subbasins45 (Fig. 5b). It is likely that volcanism occurred throughout the basins, but the morphologies of the volcanic vents are very different in the basin floor because they occurred beneath 400–1,200 m of water in the Eridania sea, resulting in 20–50 bars of pressure33. The range of volcanic styles observed proximal to the basins occurred because these volcanoes erupted subaerially (Extended Data Fig. 1).

Based on gravity measurements, the crustal thickness in the basins is 10–20 km thinner than in the immediately adjacent highlands (Fig. 5c). Strong E–W-trending remnant magnetic stripes thought to originate from magmatic intrusions at depths of tens of kilometres46,47are disrupted in the terrain where the basins occur (Fig. 5d), suggesting a connection between deep subsurface processes and surface expression. Though the basins could have an impact origin, they are unlike Martian impact basins because they have much lower circularity (average of 0.83) compared to Martian impact craters such as Gusev Crater, Jezero Crater, Gale Crater or Isidis basin (average of 0.98; Supplementary Fig. 5.); they have concave-up floors rather than flat floors and lack any typical features of impact such as raised rims, ejecta, terraced walls or central peaks.

The margins of the irregularly or rhombohedrally shaped basins are generally defined by structures including large escarpments and wrinkle ridges or thrust faults (Fig. 5c; Extended Data Fig. 4)6,29, The amount of offset of the thrusts can be difficult to determine, but minimum values reach up to 15 km calculated where circular impact or volcanic craters have been overthrust by the hanging wall (Extended Data Fig. 5). Strikingly, three large subbasins, Atlantis Chaos, Caralis Chaos and Simois Colles, are separated only by the presence of a 140-km-long WNW-trending fold complex (Fig. 6). The antiformal structure bears some characteristics of a diapir or salt dome, but no salts are detected in the THEMIS data. A three-dimensional perspective (Fig. 6b) shows steeply dipping layers flanking the structure. High-resolution images of the small impact craters on fold axes show that surface exposures of layering are not deflected by topography, indicating that the layers have near vertical dips (Extended Data Fig. 6).

Fig. 6: Folded rocks in a diapir-like structure.
figure 6

a, MOLA/HRSC elevation data draped onto THEMIS daytime infrared data illustrate that three subbasins (Atlantis, Caralis and Simois) are separated only by a large antiformal structure. Black arrows indicate dip direction of layers. b, An HRSC colour image draped over topography shows steeply dipping layers in the north flank of the antiform. Perspective is looking SE from the north side of the structure. Arrows point to the traces of steeply dipping layers.

The suite of volcanic, structural, geochemical, geophysical and topographic observations within and around the Eridania basins can be understood in the context of vertical tectonics. Before plate tectonics took hold on the Archaean Earth circa 2.5 billion years ago, crustal evolution on this planet was governed in part by vertical tectonics driven by eclogitization and delamination of the lithosphere12,13. Observational tests for convective overturn of the crust due to gravitational instabilities from Archaean terranes include eruption of thick volcanic deposits, crustal downwarping and the formation of basin-marginal diapiric-like folds13. The incorporation of hydrated, hydrothermally altered and chemically weathered volcanics into the deep crust and upper mantle via catastrophic, episodic overturn through vertical tectonics was likely an important step in evolution of felsic rocks and ultimately continental material on Earth11,12, and a similar process might have operated on early Mars, as evidenced in the Eridania region. The observation of evolved igneous rocks, structurally defined basins and folded volcanic rocks and sediments in diapiric structures are all consistent with a model of vertical recycling of the ancient crust of Mars.

The wide range of volcanic styles evident in the Eridania region suggests that similar structures, including domes, stratovolcanoes, small pyroclastic shields and calderas38, could be much more widespread in the ancient Martian crust than is currently recognized41. Evidence for diverse volcanic compositions in highland terranes also challenges the long-held paradigm that Mars is purely a basaltic world and without crustal recycling. Further, Eridania contains a wide range of likely subaqueous and coastal volcanoes analogous to the hypothesized geology of the Eoarchean or Mesoarchean Earth. The chemically and physically diverse settings present in this region of Mars provide a close analogue to both seafloor and land-based hydrothermal origin of life scenarios proposed for the Earth48,49.

Methods

Morphometric analyses

Volcanic features were identified using THEMIS daytime infrared data and Mars Context Imager images, including the global mosaic created by the Murray Lab. Where available, High Resolution Imaging Science Experiment (HiRISE) visible images were also used to describe and interpret volcanic morphologies. Visible and infrared images were combined with elevation data obtained at 200 metres per pixel resolution from the Mars Orbiter Laser Altimeter (MOLA) and the High Resolution Stereo Camera (HRSC). Elevation and slope data were instrumental in the identification of new features and essential for quantifying diameter and height characteristics.

Compositional mapping

The wavelength of the emissivity minimum between 8 and 12 µm corresponding to Si–O stretching was used to estimate silica content. To more precisely estimate the wavelength position of the surface emissivity minimum in multispectral data, we calculated a cubic spline interpolation to atmospherically corrected THEMIS emissivity spectra (bands 3–9). The wavelength position of the spline-fit minimum was then mapped on a pixel-by-pixel basis. We also applied cubic spline fits to laboratory thermal emission spectra (degraded using THEMIS filter response functions) from felsic, anorthositic, mafic and ultramafic rocks (Supplementary Information). Additional spline-fit minimum values were calculated for a suite of volcanic rocks, ranging from basalts to dacites.