46th Lunar and Planetary Science Conference (2015) 2891.pdf VOLCANIC SPREADING ON MARS: ROLE OF A BASAL DECOLLEMENT ON FAULTING AND MAGMA PROPAGATION. N. Le Corvec and P. J. McGovern, Lunar and Planetary Institute, USRA, 3600 Bay Area Blvd, Houston, TX, 77058. Introduction: Mars possesses the largest volcanic edifices in our solar system (e.g., Olympus Mons, Ascraeus Mons). These volcanoes have low slopes (<10°), a caldera at their summit, and flank terraces considered as low angle thrust faults. These structural characteristics have been related to two main processes [1, 2]: lithospheric flexure, and gravitational spreading. The extent of flank terraces varies from one volcano to another [3]. Such variation might represent differences in the geometry of the edifice, the thickness of the lithosphere, and the coefficient of friction and thus the amount of spreading. The goal is to investigate the influence of lithospheric flexure and graviational spreading on the type and extent of faulting evident at the surface of the volcano, and on the type of magmatic intrusion inside the volcanic edifice. Model: Axisymmetric finite element models using COMSOL Multiphysics [4] of volcanic loading of an elastic cone on top of an elastic lithosphere were created. For each model, a range of elastic lithosphere thicknesses Te were studied (Table 1). Finite Element Models Models Radius (km) 1 100 Initial Height (km) 5 2 370 3 Slope (°) Te (km) 2.86 20 - 30 - 40 22 3.40 75 - 100 - 125 160 12 4.29 20 - 50 - 80 4 230 14 3.48 50 - 75 - 100 5 170 18 6.04 60 - 100 - 140 Table 1: Geometric setup of the different finite elements models. In each model, the radius of the lithosphere was set to 2000 km to avoid boundary effects in the flexural response. The volcano and the lithosphere have a density of 2800 kg/m3, and share the same elastic parameters: Poisson ratio = 0.25 and Young’s Modulus E = 1011 Pa. In order to study volcanic spreading, we introduced a frictional contact between the volcanic edifice and the lithosphere. The contact was characterized by various friction coefficients (μ), from frictionless (μ = 0) to natural friction coefficients of basaltic rocks (μ=0.6) [5]. The left boundary of the model represents the symmetric axis, the upper boundary of the model is set free, rollers are applied to the right boundary of the model restricting horizontal displacement, and the bottom boundary is defined by an equilibrium balancing the displacement of the lithosphere using “Winkler” restoring forces [6] and counteracts the loading of the pre-stressed layered lithosphere [7]. Preliminary Results: We tested 5 different models, described in Table 1. The type of faulting at the surface of the volcano is mainly dependent on the ability of the volcano to spread (Fig. 1). In model 1, we observe that volcanoes with a frictionless base display circumferential thrust faults over most of their flanks. The faulting at the base of the flank is dependent on the thickness of the lithosphere (which controls the amplitude and wavelength of flexural response). Thin lithosphere will create radial normal faults at the base, whereas thicker lithospheres produce strike-slip and circumferential normal faults. Volcanoes with finite friction bases show circumferential thrust extent limited to their summit, the flank is mostly cut by radial thrust faults (Fig. 1) with their base cut by strike-slip faults only for thin lithospheres. Comparing our results with the minima elevation found on several Martian volcanoes [3], we plotted the vertical extent of circumferential thrusts against the ratio between the volcano initial height and the lithosphere thickness (Fig. 2). We observe that the largest vertical extent of circumferential thrusts along a volcano’s flank occur for the models with a frictionless base. These elevation minima are comparable with the ones found on Martian volcanoes (Fig. 2, shaded area) [3]. We observe that the vertical extent is however not a systematic function of flexural response. Thinner lithospheres produce smaller elevation minima for model 5, while it is the inverse for model 2. On the other hand, models 1, 3 and 4 show elevation minima for their median lithospheric thickness. Future Work: These results will be used to model natural example on Mars. The locations of circumferential thrusts could yield information sufficient to constrain the elastic thickness under the volcanic edifice. References: [1] Byrne P.K., et al. (2013) Geology, 41, 339-342. [2] McGovern P.J., et al. (2014) Geological Society, London, Special Publications, 401. [3] Byrne P.K., et al. (2009) Earth and Planetary Science Letters, 281, 1-13. [4] COMSOL Version 4.3b, 2013. [5] Schultz R.A. (1995) Rock Mechanics and Rock Engineering, 28, 1-15. [6] Watts A.B. (2001) Isostasy and Flexure of the Lithosphere, 458. [7] Galgana G.A., et al. (2011) J. Geophys. Res., 116, E03009. 46th Lunar and Planetary Science Conference (2015) 2891.pdf Figure 1: Plots of the differential stress and the type of faulting on the volcano’s flank as a function of distance from the volcano’s summit for model 1. Each curve represents a different coefficient of friction at the base of the volcano. The colors represent the type of faulting: cyan corresponds to circumferential thrust (CT); orange to radial thrust (RT), yellow to strike slip (SS), purple to radial normal (RN), and green to circumferential normal (CN). Figure 2: Plot of the elevation minima of circumferential thrust as a function of the ratio between the initial volcano height and the lithospheric thickness. The shaded area represents the range of elevation minima for Martian volcanoes (e.g., Arsia Mons, Ascraeus Mons, Olympus Mons) [3]. Inset: The surface of the lithosphere before flexure represents the basal reference surface (0%), the elevation of the volcano after deformation corresponds to 100%. Reworked from [3].
© Copyright 2024