Elsevier

Icarus

Volume 225, Issue 1, July 2013, Pages 538-547
Icarus

Radial dike formation on Venus: Insights from models of uplift, flexure and magmatism

https://doi.org/10.1016/j.icarus.2013.04.020Get rights and content

Highlights

  • First numerical models integrating flexural uplift with magma reservoir rupture.

  • When shallow reservoirs fail they feed vertical radial dikes and surface eruptions.

  • Resulting edifice growth creates a shallow stress barrier inhibiting dike ascent.

  • This shallow stress barrier also helps redirect magma laterally via radial dikes.

  • Model results provide one mechanism for giant radial dike swarm formation on Venus.

Abstract

Giant radiating dike swarms on Venus are geomorphic manifestations of deep-seated volcano-tectonic processes. Through an initial series of experiments we assess if these features can be explained as a product of the relationship between magmatism and the flexural stress state caused by lithospheric uplift (due to a rising mantle diapir, plume or underplating processes) using gravitationally loaded, axisymmetric, elastic finite element models. Magma chambers situated in the upper (extensional) lithosphere fail at or near the crest, redirecting magma vertically into radially aligned dikes, but these are expected to feed surface eruptions; conditions promoting lateral propagation of dikes do not occur. Magma chambers located in the lower (compressional) lithosphere are predicted to generate horizontal intrusions (sills) upon midsection failure. When both uplift (from basal loading) and magmatic surface loading (from edifice growth) are introduced in sequence at discrete times, however, conditions that promote shallow lateral propagation of magma within radial dikes are created: (1) as magma from shallow reservoirs ascends and erupts, surface-loading volcanic edifices can form; (2) as the vertical loads increase, flexing the lithosphere downward, high compressional stresses produced within the edifice create a shallow stress barrier that makes it increasingly difficult for magma to complete the process of ascent and eruption; and then (3) stresses acting on dikes ascending from a reservoir toward the surface will thus increasingly favor fracture along their lateral margins, redirecting magma through radial dikes toward distal regions at shallow depth. Our results are consistent with observed patterns of laterally extensive dike systems at volcanic centers such as Mbokomu Mons and Becuma Mons. The models thus provide a robust mechanical explanation that for the first time links lithosphere flexure and reservoir pressurization to the ascent and lateral spreading of magma through radial dike systems on Venus, thereby providing key new spatial and temporal insights into an important magmatic process operating on multiple terrestrial planets.

Introduction

Giant radial lineament systems, thought to form primarily through lateral emplacement of dikes at shallow depth (e.g., Ernst and Baragar, 1992, Grosfils and Head, 1994a, Ernst et al., 1995), are among the most intriguing volcano-tectonic structures on Venus (Fig. 1). Initial reconnaissance mapping at C1-MIDR resolution (Grosfils and Head, 1994b) revealed ∼120 such features with an average radius of ∼325 km, although some examples extend more than 2000 km; however, active global mapping efforts using higher resolution FMIDR images are revealing ∼5-fold more examples (e.g., Ernst et al., 2003, Studd et al., 2010).

Although the presence of radiating dike systems on Venus is likely (inferred through the similarities of graben swarms to dike swarms on Earth), and a link to magma chambers is compelling in many instances (e.g., Grosfils and Head, 1994a, Ernst et al., 1995), we have only limited insight into the mechanical conditions that promote the formation of laterally extensive radial dike systems. More specifically, existing flexural models that predict stress patterns conducive to radial dike formation (through simple domical uplift) do not incorporate interactions with magma chambers (e.g. Stofan et al., 1991, Cyr and Melosh, 1993). On the other hand, simple models using pressurized magma chambers that predict radial dikes propagating laterally from the reservoir do not incorporate flexural effects (e.g., McKenzie et al., 1992, Parfitt et al., 1993, Koenig and Pollard, 1998, Gudmundsson, 2006). While the stress state beneath an edifice insufficient in size to create significant flexure can provide conditions in which radial dikes ascending vertically from a ruptured chamber transition to lateral propagation at shallow depth (Hurwitz et al., 2009), consistent with field observations at many small stratocones (e.g., Poland et al., 2008), this lateral redirection is essentially constrained to the sub-edifice region, and to date no model with appropriate boundary conditions has demonstrated circumstances that promote the growth of giant regional-scale radial dike swarms such as those observed on Venus.

In a previous study, Galgana et al. (2011a) presented an analysis of magma reservoir rupture caused by the interplay between reservoir inflation and flexural effects introduced by edifice loading at the surface. Here, expanding upon these results, we explore how basal uplift and flexure affect the rupture characteristics of inflating magma reservoirs. Guided by our results and those reported previously by Galgana et al. (2011a), we then show how time-dependent introduction of basal uplift followed by edifice loading creates stress conditions within the lithosphere conducive to giant radial dike swarm formation, with a transition from vertical to lateral propagation promoted by the presence of a stress trap at the edifice base.

Section snippets

Method

Following methodology described previously in Grosfils (2007) and Galgana et al. (2011a), we use COMSOL Multiphysics (COMSOL AB, 2008) to build axisymmetric finite element models of volcanic centers within the Venusian lithosphere (Fig. 2); most are purely elastic, but a few are viscoelastic. In either case, the models extend to 900 km in radius, a sufficient distance to avoid edge-related effects. The gravitationally loaded (lithostatic prestress) lithosphere is assumed to have a basaltic

Lithospheric stresses due to uplift

Models of buoyant subsurface loading produce a flexural stress state with high σD in the upper and lower parts of the lithosphere, separated by a σD minimum in the neutral plane (Fig. 3). The upper lithosphere is predominantly under extension, while the lower lithosphere is under compression. The principal stress axis patterns show σE = σθ (out of plane) and σC = σz (vertical) throughout the upper lithosphere, predicting the formation of radially aligned structures in the upper lithosphere, either

Predicted lithosphere intrusions and magma ascent

Our models have important implications for the formation of radial dike swarms from magma chamber pressurization. In the upper portions of an uplifted lithosphere, a pressurized chamber fails at the crest area (Fig. 10), and predicted intrusion orientations are radial (Fig. 3). Thus, initial magma transport from the top of an upper lithosphere chamber will occur via vertically propagating, radially-oriented dikes (Fig. 9), a result similar to those derived from recent half-space models (

Conclusions

We present volcanologically motivated mechanical models that for the first time (to our knowledge) can explain the genesis of giant radiating dike systems found on the Venusian lithosphere, though we caution that additional explanations must continue to be sought as only about half of the radial systems on Venus are affiliated with elevated central topography (Grosfils and Head, 1994b). Our finite element models show that lithosphere fracture orientations and magma reservoir failure modes are

Acknowledgments

We thank Oded Aharonson (editor), Lionel Wilson and an anonymous reviewer for their thoughtful reviews and constructive comments that significantly improved the manuscript. This work is supported by NASA Grant NNX08AL77G from the Planetary Geology and Geophysics program and NASA CAN grant to the Lunar and Planetary Institute, USRA, with additional student involvement funded by the Summer Undergraduate Research Program of Pomona College. This is LPI Contribution No. 1728.

References (37)

  • COMSOL AB, 2008. COMSOL Multiphysics Reference Manual: Version 3.5,...
  • J. Dieterich

    Growth and persistence of Hawaiian volcanic rift zones

    J. Geophys. Res.

    (1988)
  • R. Ernst et al.

    Evidence from magnetic fabric for the flow pattern of magma in the Mackenzie giant radiating dyke swarm

    Nature

    (1992)
  • G. Galgana et al.

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

    J. Geophys. Res.

    (2011)
  • Galgana, G., McGovern, P., Grosfils, E., 2011b. The development of giant radiating dike swarms on Venus from coupled...
  • R. Grimm et al.

    The crust of Venus

  • E. Grosfils et al.

    Emplacement of a radiating dike swarm in western Vinmara Planitia, Venus: Interpretation of the regional stress field orientation and subsurface magmatic configuration

    Earth, Moon, Planets

    (1994)
  • E. Grosfils et al.

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

    Geophys. Res. Lett.

    (1994)
  • Cited by (20)

    • Magma propagation and emplacement within the central intrusive complex of Summer Coon stratovolcano, Colorado

      2021, Journal of Volcanology and Geothermal Research
      Citation Excerpt :

      The directions of the principal stresses control dike orientation, with dikes propagating perpendicular to the least compressive stress (σ3) and parallel to the maximum compressive stress (σ1; e.g., Nakamura, 1977; Rubin and Pollard, 1987). Normal stress gradients or unfavorable principal stress orientations may result in the arrest of ascending dikes near the base of an edifice and potentially a transition to sub-horizontal propagation (Pinel and Jaupart, 2000; Pinel and Jaupart, 2004; Poland et al., 2008; Hurwitz et al., 2009; Galgana et al., 2013). However, it remains unclear to what extent the overpressure (combined excess magma pressure and magma buoyancy minus the normal stress; Gudmundsson, 2002, 2012) within a dike, or host rock heterogeneities within the CIC, contributes to these mechanisms.

    • Volcano dynamics vs tectonics on Mars: evidence from Pavonis Mons

      2021, Journal of Volcanology and Geothermal Research
      Citation Excerpt :

      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.

    • Magma ascent at floor-fractured craters diagnoses the lithospheric stress state on the Moon

      2020, Earth and Planetary Science Letters
      Citation Excerpt :

      Lessons learned from Earth applied volcanology have been transferred to planetology with success. In particular, similar giant dike swarms observed on Venus have been interpreted either in terms of local plumbing system or large scale rising diapir effects (Ernst et al., 1995; Galgana et al., 2013). More recent studies of magma transport on Earth have relied not only on the magmatic intrusion shape and orientation but also on vent location to gain quantitative information on the stress state, otherwise difficult to constrain, within volcanic edifices (Maccaferri et al., 2017) or rifting areas (Maccaferri et al., 2014).

    • Hot climate inhibits volcanism on Venus: Constraints from rock deformation experiments and argon isotope geochemistry

      2017, Physics of the Earth and Planetary Interiors
      Citation Excerpt :

      Therefore, if magma transport through the lower to middle Venusian crust is dominated by diapirism, then a lower fraction of crust-situated melt can reach the surface and erupt, relative to Earth. Occasionally however, a sill may form that is large enough to generate enough uplift to initiate faulting in the brittle crust, creating a set of concentric vertical faults (see Galgana et al., 2013). If the magma reaches these faults (or vice versa) then melts can propagate upwards, forming ring-dikes or arachnoids (Head et al., 1992; Donahue and Russell, 1997; Basilevsky and Head, 2003; Wilson, 2009).

    View all citing articles on Scopus
    View full text