Volcano dynamics vs tectonics on Mars: evidence from Pavonis Mons

https://doi.org/10.1016/j.jvolgeores.2020.107148Get rights and content

Highlights

  • Narrow grabens/pit chains penetrate up to the sub-crustal magma reservoir

  • Fractal clustering of features quantify the penetration depth to 20 and 100 km

  • Dyke-induced grabens on Pavonis Mons

  • Pavonis Mons experienced active rifting at early stages of formation followed by volcano growth

Abstract

Volcanic activity is widespread within the inner Solar system and it can be commonly observed on rocky planets. In this work, we analyse the structures of Pavonis Mons in the Tharsis volcanic province of Mars by performing structural mapping, azimuth, and topographic distribution of linear features on the flanks of Pavonis, such as grabens and pit chains. We tested whether their formation is to be ascribed to the volcano dynamics and magmatic activity or the tectonics. Through the length size distribution and fractal clustering analyses of the structural features, we found that large grabens are vertically confined in the upper mechanical layers of the brittle crust whereas pit chains penetrate the whole crust up to the magmatic source, indicating that they can be considered the main feeders of Pavonis Mons. We inverted the topography with dykes and faults models to test whether grabens at the surface are the expression of intrusions at depth and we suggest that thin dykes inducing normal faulting are the most likely mechanism. Furthermore, two azimuthal distribution of the grabens are identified: concentric grabens occur on the volcano summit while linear grabens at its base show NE-SW trend as the Tharsis Mons volcanos alignment. The occurrence of linear grabens suggests that Pavonis likely experienced a phase of active rifting with the formation of such structures, followed by a phase of volcano growth and concentric magma intrusions when volcano and magma chamber dynamics prevailed.

Introduction

Planetary grabens, identified as linear topographic lows bordered by normal faults, have been widely recognized on Mars and other rocky planets as well as on the Moon (Ernst et al., 2003; Watters et al., 2012; Klimczak, 2014). However, the origin of these structures is debated (Mège et al., 2003). On Earth linear grabens may form due to dyke intrusions (Gudmundsson and Loetveit, 2005; Gudmundsson et al., 2008; Koehn et al., 2019; Wyrick and Smart, 2009), normal faulting (Melosh and Williams, 1989; Koehn et al., 2019) or a combination of both -dyke-induced normal faulting above a dyke- (Gudmundsson and Loetveit, 2005) in extensional tectonic regimes (Mège et al., 2003; Acocella and Neri, 2009). Besides, concentric dykes can also originate from a pressurized magma chamber without any active tectonics (Acocella and Neri, 2009; Bistacchi et al., 2012) or, at the regional scale, they can be caused by interaction between magmatism and lithospheric flexure due to loading (Galgana et al., 2013; Grosfils et al., 2015). Therefore, both extensional tectonics and volcano dynamics can have played an important role in the formation of Mars grabens. The formation and evolution of the Tharsis volcanic province is still a matter of debate and several hypotheses have been suggested for the formation of such a huge volcanic province, including the presence of a plume on a stationary plate (Meyzen et al., 2014) or on a decoupled lithosphere-mantle with differential rotation (Zhong, 2008), a giant impact resulting in the formation of a large igneous province (Reese et al., 2004), episodic rollback of a subducting slab (Yin, 2012).

Here we explore whether volcano dynamics or tectonic processes control deformations at Pavonis Mons by analysing in detail i) the variation of graben azimuth distribution from its base to the summit, ii) the spatial clustering of grabens and pit chains in the volcano and iii) the volcano-tectonic structures on the volcano summit. We find that the azimuthal distribution of grabens shows two main patterns: concentric grabens occur on the volcano summit and linear grabens at its base. We interpret the linear grabens at the base of the volcano as evidence of active rifting followed by a phase of volcano growth and magma resurfacing when concentric grabens formed due to the volcano dynamics.

Section snippets

Tharsis

Tharsis is the largest volcanic province of Mars and among the largest of the Solar System (Byrne, 2020). It spans an area of over 20 million km2, corresponding to ~25% of the whole Martian surface, and it comprises twelve large volcanoes from tens to hundreds of kilometres wide and up to 22 km high (e.g. Carr, 1973, Carr, 1974; Plescia, 2004; Werner, 2009; Yin, 2012). In particular, some of the shield volcanoes developed on the utmost risen part of the Tharsis bulge, creating an elevated

Data set and methods

In this work we mapped all the structural features of Pavonis Mons related to tectonics and volcano dynamics, we then analysed the fractures length distribution and the spatial distribution of structures, and we attempted to explain the topography with models of dyke injections and fault slips.

Brittle structures

Faults, grabens and pit chains characterize Pavonis Mons volcano as well as its surroundings (Scott et al., 1998). Extensional features appear to be concentric to the shape of Pavonis and concentrated at the slope change at the base of the volcano, often exhibiting very well-developed concentric sets of km-wide grabens and normal faults that can be up to 200 km-long. On the north-western and south-eastern sides of Pavonis, there are sets of grabens up to 5 km wide, more than 100 km long and a

Volcano summit zone

Pavonis Mons and its surroundings are characterized by faulting, grabens, pit chains, lava tube, lava flows as well as fan-shaped deposits at the toe of the edifice, and small volcanic fields (Shean et al., 2005; Bleacher et al., 2009). Also, the summit of Pavonis Mons shows a very complex structural setting consisting of two-nested calderas (Scott et al., 1998), with the smaller younger caldera that grew in the south-western part of the larger and older caldera (late Amazonian age; Robbins et

Spatial Analysis

The spatial analysis of grabens and pit-chains in Pavonis Mons has been performed in terms of length size distribution, self-similar clustering and azimuth distribution. Two main distributions can be used when analysing fault population statistics: size distribution and spatial distribution. The former focuses on the geometric properties of features such as fault length or displacement, while the latter analyses the properties of the whole population by relating each fault to the others, such

Topographic modelling

We analyzed the topography to determine what type of magmatic and tectonic sources may have played a role in constructing the current topography. A high-resolution DEM was extracted on a sector of Pavonis and we selected for the modelling an area with a graben over 10 km long and 500 m deep on the flank of Pavonis Mons (Fig. 11). To match the observed topography, we tested models of dyke intrusion, faulting and dyke-induced faulting above a dyke (Okada, 1992), assuming the conventional elastic

Discussion

The azimuth distribution of grabens (Fig. 3) shows different patterns depending on the elevation. Grabens developed above the 9 km of elevation, i.e. at an elevation greater than that of the caldera floor, show a concentric pattern (whose azimuths are homogeneously distributed) consistent with inflation and/or deflation dynamics of the volcano. On the other hand, grabens at an elevation lower than 9 km have a distribution with a well-defined NE-SW trend, parallel to the alignment of the Tharsis

Conclusions

In this study, we analyzed grabens and pit chains associated with Pavonis Mons in the Tharsis volcanic province on Mars. From the analysis of their topographic distribution, azimuth, size distribution and clustering we were able to discern whether they are the result of tectonics related to Tharsis province activity or internal volcano dynamics. The analyses show that faults related to large grabens are confined in a mechanical layering in the upper layers of the brittle crust, whereas deeper

Data availability

All the planetary data used in this paper are freely available at NASA PDS Geosciences Node at https://pds.nasa.gov and high-level processed products (MOLA MEGDR) are available at the USGS Astrogeology Science Center https://www.usgs.gov/centers/astrogeology-science-center. Raw CTX image processing, mosaicking and custom DTMs were generated through USGS Astrogeology ISIS3 software (https://isis.astrogeology.usgs.gov) and Ames Stereo Pipeline plugin //ti.arc.nasa.gov/tech/asr/groups/intelligent-robotics/ngt/stereo/

CRediT authorship contribution statement

Riccardo Pozzobon: Conceptualization, Investigation, Methodology, Formal analysis. Diana Orlandi: Investigation, Visualization. Carolina Pagli: Formal analysis, Investigation. Francesco Mazzarini: Conceptualization, Supervision.

Declaration of Competing Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Acknowledgements

C.P. acknowledges support by the University of Pisa grant PRA_2018_19, Italy.

References (94)

  • A. Gudmundsson

    Fracture dimensions, displacement and fluid transport

    J. Struct. Geol.

    (2000)
  • A. Gudmundsson et al.

    Dyke emplacement in a layered and faulted rift zone

    J. Volcanol. Geotherm. Res.

    (2005)
  • A. Gudmundsson et al.

    Effects of internal structure and local stresses on fracture propagation, deflection and arrest in fault zones

    J. Struct. Geol.

    (2010)
  • A. Gudmundsson et al.

    Length-displacement scaling and fault growth

    Tectonophysics

    (2013)
  • D. Koehn et al.

    Modelling of extension and dyking-induced collapse faults and fissures in rifts

    J. Struct. Geol.

    (2019)
  • F. Mazzarini et al.

    Vein development during folding in the upper brittle crust: The case of tourmaline-rich veins of eastern Elba Island, northern Tyrrhenian Sea, Italy

    J. Struct. Geol.

    (2011)
  • D. Mège et al.

    A plume tectonics model for the Tharsis province, Mars

    Planet. Space Sci.

    (1996)
  • J.B. Plescia

    Geology of the Uranius group volcanic constructs: Uranius Patera, Ceraunius Tholus, and Uranius Tholus

    Icarus

    (2000)
  • S.J. Robbins et al.

    The volcanic history of Mars: High-resolution crater-based studies of the calderas of 20 volcanoes

    Icarus

    (2011)
  • J. Ruiz et al.

    The thermal evolution of Mars as constrained by paleo-heat flows

    Icarus

    (2011)
  • K.L. Tanaka et al.

    The digital global geologic map of Mars: Chronostratigraphic ages, topographic and crater morphologic characteristics, and updated resurfacing history

    Planet. Space Sci.

    (2014)
  • J.J. Walsh et al.

    Fractal analysis of fracture pattern using the standard box-counting technique: Valid and invalid methodologies

    J. Struct. Geol.

    (1993)
  • S.C. Werner

    The global martian volcanic evolutionary history

    Icarus

    (2009)
  • D.Y. Wyrick et al.

    Dike-induced deformation and Martian graben systems

    J. Volcanol. Geotherm. Res.

    (2009)
  • T.D. Barnie et al.

    Does the lava lake of Erta ’Ale volcano respond to regional magmatic and tectonic events? An investigation using Earth Observation data

    Geol. Soc. Spec. Publ.

    (2015)
  • L. Becerril et al.

    Depth of origin of magma in eruptions

    Sci. Rep.

    (2013)
  • A. Benedicto et al.

    Layer thickness and shape of faults

    Geophys. Res. Lett.

    (2003)
  • E. Bonnet et al.

    Scaling of fracture systems in geological media

    Rev. Geophys.

    (2001)
  • S. Bouley et al.

    Late Tharsis formation and implications for early Mars

    Nature

    (2016)
  • O.P. Bour et al.

    A statistical scaling model for fracture network geometry, with validation on a multiscale mapping of a joint network (Hornelen Basin, Norway)

    J. Geophys. Res.

    (2002)
  • Buchan, K., L., Ernst, R., 2019. Dyke Swarms of the World: A Modern Perspective, Dyke Swarms of the World: A Modern...
  • P.K. Byrne

    A comparison of inner Solar System volcanism

    Nat. Astron.

    (2020)
  • P.K. Byrne et al.

    A volcanotectonic survey of Ascraeus Mons, Mars. J. Geophys. Res.

    (2012)
  • P.K. Byrne et al.

    A sagging-spreading continuum of large volcano structure. Geology

    (2013)
  • M.H. Carr

    Volcanism on Mars

    J. Geophys. Res.

    (1973)
  • M.H. Carr

    Tectonism and volcanism of the Tharsis region of Mars

    J. Geophys. Res.

    (1974)
  • P. Cervelli et al.

    Estimating source parameters from deformation data, with an application to the March 1997 earthquake swarm off the Izu Peninsula, Japan

    J. Geophys. Res. Solid Earth

    (2001)
  • A. Clauset et al.

    Power-law distributions in empirical data

    SIAM Rev.

    (2009)
  • G. Currenti et al.

    FEM-based inversion for heterogeneous fault mechanisms: application at Etna volcano by DInSAR data

    Geophys. J. Int.

    (2010)
  • C. Darcel et al.

    Connectivity properties of two-dimensional fracture networks with stochastic fractal correlation

    Water Resour. Res.

    (2003)
  • C.S. Edwards et al.

    Mosaicking of global planetary image datasets: 1. Techniques and data processing for Thermal Emission Imaging System (THEMIS) multi-spectral data

    J. Geophys. Res.

    (2011)
  • D.A. Ferrill et al.

    Dilational fault slip and pit chain formation on Mars

    GSA Today

    (2004)
  • G.A. Galgana et al.

    Evolution of large Venusian volcanoes: insights from coupled models of lithospheric flexure and magma reservoir pressurization

    J. Geophys. Res.

    (2011)
  • E.B. Grosfils et al.

    The global distribution of giant radiating dike swarms on Venus: Implications for the global stress state

    Geophys. Res. Lett.

    (1994)
  • E.B. Grosfils et al.

    Elastic models of magma reservoir mechanics: A key tool for investigating planetary volcanism

    Geol. Soc. Spec. Publ.

    (2015)
  • A. Gudmundsson

    Fluid pressure and stress drop in fault zones

    Geophys. Res. Lett.

    (1999)
  • A. Gudmundsson et al.

    2013. Relations between the scaling exponents, entropies, and energies of fracture networks. Bull.Géol

    France

    (2013)
  • Cited by (4)

    • Passive magmatism on Earth and Earth-like planets

      2022, Geosystems and Geoenvironment
    View full text