Model Estimates of Non-Hydrostatic Stresses in the Martian Crust and Mantle: 1—Two-Level Model

Abstract

Regions of maximum shear and tension–compression stresses in the Martian interior have been revealed using two types of models: the elastic model and the model with an elastic lithosphere of varied thickness (150–500 km) positioned on a weak layer that has partially lost its elastic properties. The weakening is simulated by a ten-fold lower value of the shear modulus down to the core boundary. The numerical simulation applies Green’s functions (load number method) with the step of 1 × 1 grade along latitude and longitude down to a depth of 1000 km. The boundary condition is the expansion of the latest data on Martian topography and the gravitational field (model MRO120D) in spherical harmonics up to the degree and order of 90 in relation to the reference surface that is assumed an equilibrium spheroid. The considered two-level compensation model assumes nonequilibrium relief and density anomalies at the crust–mantle boundary to be the sources of the anomalous gravitational field. Calculations are performed for two test models of Martian internal structure with the crust mean thicknesses of 50 to 100 km and mean density of 2900 kg/m³. Considerable tangential and simultaneously compressive stresses occur under the Tharsis region. The main regions of high shear and simultaneously extentional stresses are located in the Hellas region crust and in the lithosphere of the following regions: Argyre Planitia, Mare Acidalium, Arcadia Planitia and Valles Marineris. The zone of high maximum shear and extentional stresses has been found at the base of the lithosphere under the Olympus volcano and that under the Elysium rise.