Planetary Contraction Research Articles

Lithospheric contraction concentric to Tharsis: 3D structural modeling of large thrust faults between Thaumasia highlands and Aonia Terra, Mars

Large thrust faults on Mars are caused by lithospheric planetary contraction. The geometry of these faults is linked with the mechanical behavior of the lithosphere. Tharsis, the largest volcano-tectonic province on Mars, controls the global tectonic pattern of the planet. Here, we present a study of five large thrust faults concentric to Tharsis, located between the Thaumasia Highlands and the Argyre impact basin. We applied a 3D structural modeling, using a combination of fault-parallel flow and trishear algorithms to estimate the geometry and kinematics of the faults at depth. The modeled faults show an upper planar part dipping 33° to 40°, rooting with a listric geometry into horizontal levels at 13–27 km depth, with fault slips of 801–3366 m. The general out-of-Tharsis vergence, the listric fault geometries and the deepening of the depth of faulting toward Thaumasia outline an incipient thrust wedge architecture. Assuming that the largest faults rooted at the Brittle-Ductile Transition, we calculate a heat flow at the time of faulting of 24–54 mW m−2. The obtained strength envelopes for dry and wet conditions show that all the strength of the lithosphere was located in the upper half of the crust.

A Morphometric Investigation of Large‐Scale Crustal Shortening on Mars

AbstractMars' surface exhibits abundant topographic expressions of large thrust fault‐related folds that have been attributed to global planetary contraction. Morphometric analyses of such structures provide insight into their growth history. With global THEMIS imagery and HRSC–MOLA topographic data, 49 thrusts with lengths between 35 and 544 km were mapped across Mars' surface. Assuming planar fault geometries with dips of 30°, the average maximum displacement‐length ratio (Dmax/L) of these structures is 6.1 × 10−3 ± 1.4 × 10−3, with smaller ratios observed for faults within the northern lowlands (2.9 × 10−3 ± 0.9 × 10−3) compared to the southern highlands (9.2 × 10−3 ± 1.9 × 10−3). However, these differences may be accounted for if mechanical layering in the northern lowland crust promotes either a shallowing of the fault dip angle relative to the southern highlands or the development of ramp‐flat geometries such that the topographic scarp height may under‐estimate the total fault displacement or a combination of these two scenarios together. Alternatively, these Dmax/L patterns may reflect hemispheric differences in the brittle‐ductile transition (BDT) depth; however, the observed pattern is stratigraphically inconsistent with the Martian crustal dichotomy, whereby the northern lowlands have thinner (or denser) crust and therefore presumably a deeper BDT than the southern highlands. Fault displacement‐length profiles are commonly asymmetric, with multiple local minima observed along their lengths. Spectral analysis of these profiles, using Fourier‐ and S‐Transforms, indicates power at a range of spatial frequencies, reflecting complex growth and linkage histories, with peak spectral frequency, or number of segments, being negatively correlated with the Dmax/L ratios.

Did the cessation of convection in Mercury's mantle allow for a dynamo supporting increase in heat loss from its core?

Mercury's large core generates a magnetic field and harbors a solidified inner component. These features constrain models of the planet's thermal history. The mantle provides the boundary condition on Mercury's core that determines heat loss. Recent studies suggest Mercury's mantle may have a higher thermal diffusivity than the silicate shells of other rocky bodies due to iron depletion. Considering the role of diffusivity, we model Mercury's mantle starting from a post magma-ocean state by calculating core-mantle boundary heat flow and the ratio of inner-core boundary radius to core-mantle boundary radius, fc, for periods comparable to the age of the solar system. Core-mantle boundary (CMB) heat flow is calculated using 3D mantle convection simulations for thermal diffusivities, κ, ranging from 1.0 - 3.0 ×10−6 m2/s and initial radiogenic uniform mantle heating rates, χ, of 0 - 40 pW/kg (that decay with a 3 Gyr half-life). Several scenarios can unfold for the range of parameters considered: these include cases featuring both the cessation of mantle convection and its continuation at 4.5 Gyr. We map a trend in κ - χ space and find that for some parameters present-day core heat flux and inner-core size estimates are satisfied for a Mercurian mantle that is cooling by conduction. The influence of initial sulfur fraction in the core was examined for a subset of cases. For a sulfur fraction of 0.10, fc falls between 0.1 and 0.55 which corresponds to planetary radius contractions below 7 km (since the Late Heavy Bombardment). However, fc exceeds 0.55 for a lower initial sulfur fraction and results in planetary contraction in excess of 7 km. In general, we find that the mean core heat flux reaches a temporal local minimum when the mantle transitions from a convective to a conductive regime and subsequently climbs before decreasing. The transition to conduction is delayed with increased mantle internal heating rate but the maximum mean heat flux from the core into the conducting and cooling mantle is always greater than the heat flux observed at the cessation of the stagnant-lid convection. We find that for an initially more vigorously convecting mantle the onset of inner-core growth is earlier while the basal heat flux is marginally reduced at present day. In general, we find that current-day conductive cooling of the Mercurian mantle can satisfy estimates on Mercury's core heat loss inferred from the strength of its magnetic field while also satisfying the limits on the present-day size of its inner core.

Present‐Day Mars' Seismicity Predicted From 3‐D Thermal Evolution Models of Interior Dynamics

AbstractThe Interior Exploration using Seismic Investigations, Geodesy and Heat Transport mission, to be launched in 2018, will perform a comprehensive geophysical investigation of Mars in situ. The Seismic Experiment for Interior Structure package aims to detect global and regional seismic events and in turn offer constraints on core size, crustal thickness, and core, mantle, and crustal composition. In this study, we estimate the present‐day amount and distribution of seismicity using 3‐D numerical thermal evolution models of Mars, taking into account contributions from convective stresses as well as from stresses associated with cooling and planetary contraction. Defining the seismogenic lithosphere by an isotherm and assuming two end‐member cases of 573 K and the 1073 K, we determine the seismogenic lithosphere thickness. Assuming a seismic efficiency between 0.025 and 1, this thickness is used to estimate the total annual seismic moment budget, and our models show values between 5.7 × 1016 and 3.9 × 1019 Nm.

The interior configuration of planet Mercury constrained by moment of inertia and planetary contraction

This paper presents an analysis of present-day interior configuration models for Mercury considering cores of Fe-S or Fe-Si alloy, the latter possibly covered by a solid FeS layer, in light of the improved limit of planetary contraction of 7 km derived from MErcury Surface, Space ENvironment, GEochemistry, and Ranging observations of surface landforms. Density profiles, generated by a Monte Carlo approach, are constrained by Mercury's mass, polar moment of inertia (C), fraction of polar moment corresponding to its outer solid shell (Cm/C), and observed planetary contraction. Results show that the outer liquid core boundary is constrained to 1985–2090 km in radius, where large radii correspond to high Si and S core contents and high mantle densities or the presence of an FeS layer at the top of the outer core. The bulk core S and Si contents are within 2.8–8.9 wt % and above 8.5 wt %, respectively, where an increase of light element core content correlates positively with mantle density and core size. The size of the inner core is constrained by the observed planetary contraction to below 1454 or 1543 km in radius for bulk cores rich in S (near 8.9 wt %) or Si (near 25 wt %), respectively. For cores poor in light elements, inner cores up to 1690 km in radius remain consistent with the observed planetary contraction. Finally, we show that solid FeS at outer core conditions, previously argued to float on liquid Fe-S, may be denser than the residual liquid. This implies that a separate mechanism may be required to maintain an FeS layer at the suggested location.

Thermochemical evolution of Mercury's interior

A number of observations performed by the MESSENGER spacecraft can now be employed to better understand the evolution of Mercury's interior. Using recent constraints on interior structure, surface composition, volcanic and tectonic histories, we modeled the thermal and magmatic evolution of the planet. We ran a large set of Monte Carlo simulations based on one‐dimensional parametrized models, spanning a wide range of parameters. We complemented these simulations with selected calculations in 2‐D cylindrical and 3‐D spherical geometry, which confirmed the validity of the parametrized approach and allowed us to gain additional insight into the spatiotemporal evolution of mantle convection. Core radii of 1940 km, 2040 km, and 2140 km have been considered, and while in the first two cases several models satisfy the observational constraints, no admissible models were found for a radius of 2140 km. A typical thermal evolution scenario consists of an initial phase of mantle heating accompanied by planetary expansion and the production of a substantial amount of partial melt. The evolution subsequent to 2 Gyr is characterized by secular cooling that proceeds approximately at a constant rate and implies that planetary contraction should be ongoing today. Most of the models predict mantle convection to cease after 3–4 Gyr, indicating that Mercury may be no longer dynamically active. Finally, assuming the observed surface abundance of radiogenic elements to be representative for the entire crust, we determined bulk silicate concentrations of 35–62 ppb Th, 20–36 ppb U, and 290–515 ppm K, similar to those of other terrestrial planets.

The effect of composition on the evolution of giant and intermediate-mass planets

We model the evolution of planets with various masses and compositions. We investigate the effects of the composition and its depth dependence on the long-term evolution of the planets. The effects of opacity and stellar irradiation are also considered. It is shown that the change in radius due to various compositions can be significantly smaller than the change in radius caused by the opacity. Irradiation also affects the planetary contraction but is found to be less important than the opacity effects. We suggest that the mass-radius relationship used for characterization of observed extrasolar planets should be taken with great caution since different physical conditions can result in very different mass-radius relationships.

OHMIC DISSIPATION IN THE ATMOSPHERES OF HOT JUPITERS

Hot Jupiter atmospheres exhibit fast, weakly-ionized winds. The interaction of these winds with the planetary magnetic field generates drag on the winds and leads to ohmic dissipation of the induced electric currents. We study the magnitude of ohmic dissipation in representative, three-dimensional atmospheric circulation models of the hot Jupiter HD 209458b. We find that ohmic dissipation can reach or exceed 1% of the stellar insolation power in the deepest atmospheric layers, in models with and without dragged winds. Such power, dissipated in the deep atmosphere, appears sufficient to slow down planetary contraction and explain the typically inflated radii of hot Jupiters. This atmospheric scenario does not require a top insulating layer or radial currents that penetrate deep in the planetary interior. Circulation in the deepest atmospheric layers may actually be driven by spatially non-uniform ohmic dissipation. A consistent treatment of magnetic drag and ohmic dissipation is required to further elucidate the consequences of magnetic effects for the atmospheres and the contracting interiors of hot Jupiters.

East–west faults due to planetary contraction

Contraction, expansion and despinning have been common in the past evolution of Solar System bodies. These processes deform the lithosphere until it breaks along faults. Their characteristic tectonic patterns have thus been sought for on all planets and large satellites with an ancient surface. While the search for despinning tectonics has not been conclusive, there is good observational evidence on several bodies for the global faulting pattern associated with contraction or expansion, though the pattern is seldom isotropic as predicted. The cause of the non-random orientation of the faults has been attributed either to regional stresses or to the combined action of contraction/expansion with another deformation (despinning, tidal deformation, reorientation). Another cause of the mismatch may be the neglect of the lithospheric thinning at the equator or at the poles due either to latitudinal variation in solar insolation or to localized tidal dissipation. Using thin elastic shells with variable thickness, I show that the equatorial thinning of the lithosphere transforms the homogeneous and isotropic fault pattern caused by contraction/expansion into a pattern of faults striking east–west, preferably formed in the equatorial region. By contrast, lithospheric thickness variations only weakly affect the despinning faulting pattern consisting of equatorial strike-slip faults and polar normal faults. If contraction is added to despinning, the despinning pattern first shifts to thrust faults striking north–south and then to thrust faults striking east–west. If the lithosphere is thinner at the poles, the tectonic pattern caused by contraction/expansion consists of faults striking north/south. I start by predicting the main characteristics of the stress pattern with symmetry arguments. I further prove that the solutions for contraction and despinning are dual if the inverse elastic thickness is limited to harmonic degree two, making it easy to determine fault orientation for combined contraction and despinning. I give two methods for solving the equations of elasticity, one numerical and the other semi-analytical. The latter method yields explicit formulas for stresses as expansions in Legendre polynomials about the solution for constant shell thickness. Though I only discuss the cases of a lithosphere thinner at the equator or at the poles, the method is applicable for any latitudinal variation of the lithospheric thickness. On Iapetus, contraction or expansion on a lithosphere thinner at the equator explains the location and orientation of the equatorial ridge. On Mercury, the combination of contraction and despinning makes possible the existence of zonal provinces of thrust faults differing in orientation (north–south or east–west), which may be relevant to the orientation of lobate scarps.

The Interior Structure, Composition, and Evolution of Giant Planets

We discuss our current understanding of the interior structure and thermal evolution of giant planets. This includes the gas giants, such as Jupiter and Saturn, that are primarily composed of hydrogen and helium, as well as the "ice giants," such as Uranus and Neptune, which are primarily composed of elements heavier than H/He. The effect of different hydrogen equations of state (including new first-principles computations) on Jupiter's core mass and heavy element distribution is detailed. This variety of the hydrogen equations of state translate into an uncertainty in Jupiter's core mass of 18 M_Earth. For Uranus and Neptune we find deep envelope metallicities up to 0.95, perhaps indicating the existence of an eroded core, as also supported by their low luminosity. We discuss the results of simple cooling models of our solar system's planets, and show that more complex thermal evolution models may be necessary to understand their cooling history. We review how measurements of the masses and radii of the ~50 transiting extrasolar giant planets are changing our understanding of giant planets. In particular a fraction of these planets appear to be larger than can be accommodated by standard models of planetary contraction. We review the proposed explanations for the radii of these planets. We also discuss very young giant planets, which are being directly imaged with ground- and space-based telescopes.

Mercury: the enigmatic innermost planet

The planet Mercury, a difficult object for study by astronomical observation and spacecraft exploration alike, poses multiple challenges to our general understanding of the inner planets. Mercury’s anomalously high uncompressed density implies a metal fraction of 60% or more by mass, an extreme outcome of planetary formational processes common to the inner solar system. Whether that outcome was the result of chemical gradients in the early solar nebula or removal by impact or vaporization of most of the silicate shell from a differentiated protoplanet can potentially be distinguished on the basis of the chemical composition of the present crust. Our understanding of the geological evolution of Mercury and how it fits within the known histories of the other terrestrial planets is restricted by the limited coverage and resolution of imaging by the only spacecraft to have visited the planet. The role of volcanism in Mercury’s geological history remains uncertain, and the dominant tectonic structures are lobate scarps interpreted as recording an extended episode of planetary contraction, issues that require global imaging to be fully examined. That Mercury has retained a global magnetic field when larger terrestrial planets have not stretches the limits of standard hydromagnetic dynamo theory and has led to proposals for a fossil field or for exotic dynamo scenarios. Hypotheses for field generation can be distinguished on the basis of the geometry of Mercury’s internal field, and the existence and size of a fluid outer core on Mercury can be ascertained from measurements of the planet’s spin axis orientation and gravity field and the amplitude of Mercury’s forced librations. The nature of Mercury’s polar deposits, suggested to consist of volatile material cold-trapped on the permanently shadowed floors of high-latitude impact craters, can be tested by remote sensing of the composition of Mercury’s surface and polar atmosphere. The extremely dynamic exosphere, which includes a number of species derived from Mercury’s surface, offers a novel laboratory for exploring the nature of the complex and changing interactions among the solar wind, a small magnetosphere, and a solid planet. Recent ground-based astronomical measurements and several new theoretical developments set the stage for the in-depth exploration of Mercury by two spacecraft missions within the coming decade.

The tectonic and volcanic history of mercury as inferred from studies of scarps, ridges, troughs, and other lineaments

This paper consists of separate descriptive and interpretive sections which address the issues of tectonic and volcanic modification of the Mercurian surface. In section 1 a nongenetic classification system based on planimetric form is outlined and applied to scarps, ridges, troughs, and other lineaments on Mercury. Nine landform classes are defined, and type examples of each are presented and discussed. The resulting set of observational constraints provides the basis for the following evolutionary chronology for Mercury developed in section 2. After formation of a solid lithosphere but before the end of heavy meteoritic bombardment, tidal spindown and planetary contraction owing to cooling combined to produce a global pattern of linear fractures and arcuate thrust faults. Throughout this period, intercrater plains formed by extrusion from linear vents or by local lithospheric melting. Near the end of heavy meteoritic bombardment the Caloris basin‐forming event drastically altered the morphology of much of the Mercurian surface. Soon after basin excavation, fluid material covered its floor to form smooth plains. Planimetrically irregular ridges formed inside Caloris in response to compressive stresses induced by regional subsidence. Gravitational readjustment to emplacement of plains material outside Caloris was responsible for scarp and ridge formation there. Uplift of plains inside Caloris produced extensional fractures which transect the previously established pattern of ridges. At roughly the same time, intense modification of preexisting landforms occurred in a region roughly antipodal to Caloris Basin to produce a distinctive terrain unit characterized by scattered hills and linear, open‐ended troughs. By this time, volcanic plains emplacement was largely restricted inside large, fresh craters. Irregular scarps and ridges developed in these plains in response to local stresses, perhaps induced by magma redistribution at depth. The period following the most recent smooth plains formation on Mercury has been characterized by relatively light cratering and only mild landform modification.

Some aspects of core formation in Mercury

Evidence for a large metallic core in Mercury is all indirect; the internal magnetic field may imply a convective dynamo; the surface geology is suggestive of large-scale differentiation; and thermal history calculations based on cosmochemical models for Mercury predict core formation. The presence of a core will not be confirmed by an accurate measure of J 2, probably two orders of magnitude larger than the hydrostatic value. Core infall on Mercury would be accompanied by an increase in planetary radius of 17 km, an increase in the mean internal temperature of 700°K, and substantial melting of the mantle. If Mercury's core differentiated from an originally homogeneous planet, such an event must predate most of the present surface. Subsequent to core formation, Mercury's radius has decreased by about 2km due primarily to cooling of the lithosphere, in agreement with the estimate by Strom and others of the amount of planetary contraction based on photogeologic measurements of length, dip and throw on the global system of lobate scarps. A convective dynamo mechanism for Mercury's magnetic field is in apparent conflict with cosmochemical models that do not predict a substantial source of heat, most probably radiogenic, in Mercury's core. Without such a heat source, the core would solidify within about 1 b.y. after core infall, producing an unacceptably large contraction in Mercury's radius.