ALBERT

Description: We use 1D thermal history models and 3D numerical experiments to study the impact of dynamic thermal disequilibrium and large temporal variations of normal and shear stresses on the initiation of plate tectonics. Previous models that explored plate tectonics initiation from a steady state, single plate mode of convection concluded that normal stresses govern the initiation of plate tectonics, which based on our 1D model leads to plate yielding being more likely with increasing interior heat and planet mass for a depth-dependent Byerlee yield stress. Using 3D spherical shell mantle convection models in an episodic regime allows us to explore larger temporal stress variations than can be addressed by considering plate failure from a steady state stagnant lid configuration. The episodic models show that an increase in convective mantle shear stress at the lithospheric base initiates plate failure, which leads with our 1D model to plate yielding being less likely with increasing interior heat and planet mass. In this out-of-equilibrium and strongly time-dependent stress scenario, the onset of lithospheric overturn events cannot be explained by boundary layer thickening and normal stresses alone. Our results indicate that in order to understand the initiation of plate tectonics, one should consider the temporal variation of stresses and dynamic disequilibrium.

Description: Abstract Thermal history models, historically used to understand Earth's geologic history, are being coupled to climate models to map conditions that allow planets to maintain life. However, the lack of structural uncertainty assessment has blurred guidelines for how thermal history models can be used toward this end. Structural uncertainty is intrinsic to the modeling process. Model structure refers to the cause and effect relations that define a model and are assumed to adequately represent a particular real world system. Intrinsic/structural uncertainty is different from input and parameter uncertainties (which are often evaluated for thermal history models). A full uncertainty assessment requires that input/parametric and intrinsic/structural uncertainty be evaluated (one is not a substitute for the other). We quantify the intrinsic uncertainty for several parameterized thermal history models (a subclass of planetary models). We use single perturbation analysis to determine the reactance time of different models. This provides a metric for how long it takes low‐amplitude, unmodeled effects to decay or grow. Reactance time is shown to scale inversely with the strength of the dominant model feedback (negative or positive). A perturbed physics analysis is then used to determine uncertainty shadows for model outputs. This provides probability distributions for model predictions. It also tests the structural stability of a model (do model predictions remain qualitatively similar, and within assumed model limits, in the face of intrinsic uncertainty?). Once intrinsic uncertainty is accounted for, model outputs/predictions and comparisons to observational data should be treated in a probabilistic way.

Description: Consideration of the structure of dynamical equilibria in terrestrial planets using simplified descriptions of the relevant heat transport processes (rigid-lid convection, plate tectonics, heat-pipe volcanism) reveals that if the efficiency of plate-tectonic heat transport decreases at higher mantle temperature, then it cannot govern quasi-equilibrium dynamical evolution, and the system is always evolving away from the plate-tectonic regime. A planet on which plate tectonics is less efficient at higher temperature stays in heat-pipe mode longer, spends less time undergoing plate tectonics, and has a low and ever-decreasing Urey number during this phase. These conclusions are based solely on the structure of the equilibria in a system with less efficient plate tectonics in the past and are independent of the mechanisms leading to this behavior. Commonly used quasi-equilibrium approaches to planetary thermal evolution are likely not valid for planets in which heat transport becomes less efficient at higher temperature.

Description: The discovery of large terrestrial (~1 Earth mass (M e ) to 〈 10 M e ) extrasolar planets has prompted a debate as to the likelihood of plate tectonics on these planets. Canonical models assume classic basal heating scaling relationships remain valid for mixed heating systems with an appropriate internal temperature shift. Those scalings predict a rapid increase of convective velocities ( V rms ) with increasing Rayleigh numbers ( Ra ) and non-dimensional heating rates ( Q ). To test this we conduct a sweep of 3-D numerical parameter space for mixed heating convection in isoviscous spherical shells. Our results show that while V rms increases with increasing thermal Ra it does so at a slower rate than predicted by bottom heated scaling relationships. Further, the V rms decreases asymptotically with increasing Q . These results show that independent of specific rheologic assumptions (e.g., viscosity formulations, water effects, lithosphere yielding), the differing energetics of mixed and basally heated systems can explain the discrepancy between different modeling groups. High temperature, or young, planets with a large contribution from internal heating will operate in different scaling regimes compared to cooler temperature, or older, planets that may have a larger relative contribution from basal heating. Thus, differences in predictions as to the likelihood of plate tectonics on exoplanets may well result from different models being more appropriate to different times in the thermal evolution of a terrestrial planet (as opposed to different rheologic assumptions as has often been assumed).

Description: We use a suite of 3-D numerical experiments to test and expand 2-D planar isoviscous scaling relationships of Moore [ 2008 ] for mixed heating convection in spherical geometry mantles. The internal temperature scaling of Moore [ 2008 ], when modified to account for spherical geometry, matches our experimental results to a high degree of fit. The heatflux through the boundary layers scale as a linear combination of internal (Q) and basal heating and the modified theory predictions match our experimental results. Our results indicate that boundary layer thickness and surface heat flux are not controlled by a local boundary layer stability condition (in agreement with the results of Moore [ 2008 ]), and are instead strongly influenced by boundary layer interactions. Subadiabtic mantle temperature gradients, in spherical 3D, are well described by a vertical velocity scaling based on discrete drips as opposed to a scaling based on coherent sinking sheets, which was found to describe 2D planar results. Root Mean Square (RMS) velocities are asymptotic for both low Q and high Q, with a region of rapid adjustment between asymptotes for moderate Q. RMS velocities are highest in the low Q asymptote, and decrease as internal heating is applied. The scaling laws derived by Moore [ 2008 ], and extended here, are robust and highlight the importance of differing boundary layer processes acting over variable Q and moderate Ra.

Description: SUMMARY Boundary layer theory is used to derive scaling relationships for plate stresses in a mantle convection system with a low-viscosity asthenosphere. The theory assumes a plate tectonic like mode of mantle convection with flow driven by an active upper boundary layer. The theory predicts that the confinement of horizontal mantle flow within a low-viscosity, sublithospheric channel can lead to an increase in plate stress compared to the case lacking a channel (even if the absolute viscosity of the sublithosphere mantle does not change between the two cases). The theory further predicts increasing shear stress with decreasing low-viscosity channel thickness. If the thickness of tectonic plates is determined dominantly by a dehydrated chemical lithosphere, then the plate normal stress is predicted to also increase with decreasing channel thickness. We use 3-D spherical shell simulations of mantle convection with temperature-, depth- and stress dependent rheology to test scaling trends. The simulations and theoretical scalings demonstrate that a low-viscosity layer (asthenosphere) can amplify convective stresses. If the level of convective stress plays a role in maintaining and/or reactivating plate boundaries, this suggests that a relatively thin low viscosity layer may help to maintain plate tectonics. The numerical simulations support this suggestion as they show that an increase in the thickness of a low viscosity channel can cause the system to transition from an active-lid mode of convection to a stagnant lid state. Collectively, the simulations and theoretical scalings lead to the conclusion that the role of the asthenosphere in maintaining plate tectonics does not come principally from a basal lubrication effect, associated with a low absolute asthenosphere viscosity, but, instead, from a mantle flow channelization effect, associated with a high viscosity contrast from the asthenosphere to the mantle below.

Description: We use simulations of mantle convection with surface yielding to show that multiple tectonic regimes are possible for equivalent system parameter values. Models with the same lithospheric strength parameters and the same vigor of convection can display different modes of tectonics. Within the region of multiple solutions, the evolutionary pathway of the system is the dominant factor that determines the tectonic mode (e.g., whether mantle convection operates in a plate tectonic like mode). The extent of the multiple regimes window is found to increase with the temperature-dependent viscosity contrast across the system. The implication for models that seek to predict the tectonic regimes of planets is that the temporal evolution of the planet needs to be taken into account. A further implication is that modeling studies can lead to different conclusions regarding the tectonic state of a planet, extra-solar planets in particular, despite the final model parameter values remaining equivalent.

Description: Super-continental insulation refers to an increase in mantle temperature below a supercontinent due to the heat transfer inefficiency of thick, stagnant continental lithosphere relative to thinner, subducting oceanic lithosphere. We use thermal network theory, numerical simulations, and laboratory experiments to provide tighter physical insight into this process. We isolate two end-member dynamic regimes. In the thermally well mixed regime the insulating effect of continental lithosphere can not cause a localized increase in mantle temperature due to the efficiency of lateral mixing in the mantle. In this regime the potential temperature of the entire mantle is higher than it would be without continents, the magnitude depending on the relative thickness of continental and oceanic lithosphere (i.e., the insulating effects of continental lithosphere are communicated to the entire mantle). Thermal mixing can be short circuited if subduction zones surround a supercontinent or if the convective flow pattern of the mantle becomes spatially fixed relative to a stationary supercontinent. This causes a transition to the thermal isolation regime: The potential temperature increases below a supercontinent whereas the potential temperature below oceanic domains drops such that the average temperature of the whole mantle remains constant. Transition into this regime would thus involve an increase in the suboceanic viscosity, due to local cooling, and consequently a decrease in the rate of oceanic lithosphere overturn. Transition out of this regime can involve the unleashing of flow driven by a large lateral temperature gradient, which will enhance global convective motions. Our analysis highlights that transitions between the two states, in either direction, will effect not only the mantle below a supercontinent but also the mantle below oceanic regions. This provides a larger set of predictions that can be compared to the geologic record to help determine if a hypothesized super-continental thermal effect did or did not occur on our planet.

Description: Tectonic plate motions reflect dynamical contributions from subduction processes (i.e., classical “slab-pull” forces) and lateral pressure gradients within the asthenosphere (“asthenosphere-drive” forces), which are distinct from gravity forces exerted by elevated mid-ocean ridges (i.e., classical “ridge-push” forces). Here we use scaling analysis to show that the extent to which asthenosphere-drive contributes to plate motions depends on the lateral dimension of plates and on the relative viscosities and thicknesses of the lithosphere and asthenosphere. Whereas slab-pull forces always govern the motions of plates with a lateral extent greater than the mantle depth, asthenosphere-drive forces can be relatively more important for smaller (shorter wavelength) plates, large relative asthenosphere viscosities or large asthenosphere thicknesses. Published plate velocities, tomographic images and age-binned mean shear wave velocity anomaly data allow us to estimate the relative contributions of slab-pull and asthenosphere-drive forces for the motions of the Atlantic and Pacific plates. Whereas the Pacific plate is driven largely by slab pull, the Atlantic plate is predicted to be strongly driven by basal forces related to viscous coupling to strong asthenospheric flow, consistent with recent observations related to the stress state of North America. In addition, compared to the East Pacific Rise (EPR), the relatively large lateral pressure gradient near the Mid-Atlantic Ridge (MAR) is expected to produce significantly steeper dynamic topography. Thus, the relative importance of this plate-driving force may partly explain why the flanking topography at the EPR is smoother than at the MAR. Our analysis also indicates that this plate-driving force was more significant, and heat loss less efficient, in Earth's hotter past compared with its cooler present state. This type of trend is consistent with thermal history modeling results which require less efficient heat transfer in Earth's past.

Description: Abstract The current goals of the astrobiology community are focused on developing a framework for the detection of biosignatures, or evidence thereof, on objects inside and outside of our solar system. A fundamental aspect of understanding the limits of habitable environments (surface liquid water) and detectable signatures thereof is the study of where the boundaries of such environments can occur. Such studies provide the basis for understanding how a once inhabitable planet might come to be uninhabitable. The archetype of such a planet is arguably Earth's sibling planet, Venus. Given the need to define the conditions that can rule out bio‐related signatures of exoplanets, Venus provides a unique opportunity to explore the processes that led to a completely uninhabitable environment by our current definition of the term. Here we review the current state of knowledge regarding Venus, particularly in the context of remote‐sensing techniques that are being or will be employed in the search for and characterization of exoplanets. We discuss candidate Venus analogs identified by the Kepler and TESS exoplanet missions and provide an update to exoplanet demographics that can be placed in the potential runaway greenhouse regime where Venus analogs are thought to reside. We list several major outstanding questions regarding the Venus environment and the relevance of those questions to understanding the atmospheres and interior structure of exoplanets. Finally, we outline the path toward a deeper analysis of our sibling planet and the synergy to exoplanetary science.