Earth's climate, mantle, and core interact over geologic time scales. Climate influences whether plate tectonics can take place on a planet, with cool climates being favorable for plate tectonics because they enhance stresses in the lithosphere, suppress plate boundary annealing, and promote hydration and weakening of the lithosphere. Plate tectonics plays a vital role in the long-term carbon cycle, which helps to maintain a temperate climate. Plate tectonics provides long-term cooling of the core, which is vital for generating a magnetic field, and the magnetic field is capable of shielding atmospheric volatiles from the solar wind. Coupling between climate, mantle, and core can potentially explain the divergent evolution of Earth and Venus. As Venus lies too close to the sun for liquid water to exist, there is no long-term carbon cycle and thus an extremely hot climate. Therefore, plate tectonics cannot operate and a long-lived core dynamo cannot be sustained due to insufficient core cooling. On planets within the habitable zone where liquid water is possible, a wide range of evolutionary scenarios can take place depending on initial atmospheric composition, bulk volatile content, or the timing of when plate tectonics initiates, among other factors. Many of these evolutionary trajectories would render the planet uninhabitable. However, there is still significant uncertainty over the nature of the coupling between climate, mantle, and core. Future work is needed to constrain potential evolutionary scenarios and the likelihood of an Earth-like evolution.

The paleomagnetic record indicates the geodynamo has been active over much of Earth history with surprisingly little trend in paleointensity. Variability, however, is expected from models that predict a sharp increase in intensity following inner core nucleation (ICN) and implied by Neoproterozoic anomalies that hint at a highly variable field over several hundred million years. Here we demonstrate with a suite of numerical dynamos driven by a new thermal evolution model that the geodynamo could have transitioned from a multipolar to dipolar regime around 1.7Ga, then to a weak-field dynamo around 1.0Ga, and finally to a strong-field dipole following ICN around 650Ma that is maintained to the present day. The occurrence of a weak-field geodynamo in the Neoproterozoic may be consistent with the observed anomalous apparent polar wander paths and reversal behavior. Recovery to a dipolar geodynamo in the Phanerozoic could be a signature of inner core nucleation. Index terms: 1507, 1560, and 1521.

The origin of Earth's ancient magnetic field is an outstanding problem. It has recently been proposed that exsolution of MgO from the core may provide sufficient energy to drive an early geodynamo. Here we present new experiments on Mg partitioning between iron-rich liquids and silicate/oxide melts. Our results indicate that Mg partitioning depends strongly on the oxygen content in the iron-rich liquid, in contrast to previous findings that it depends only on temperature. Consequently, MgO exsolution during core cooling is drastically reduced and insufficient to drive an early geodynamo alone. Using the new experimental data, our thermal model predicts inner core nucleation at similar to 850 Ma and a nearly constant paleointensity.

The paleomagnetic record is central to our understanding of the history of the Earth. The orientation and intensity of magnetic minerals preserved in ancient rocks indicate the geodynamo has been alive since at least the Archean and possibly the Hadean. A paleomagnetic signature of the initial solidification of the inner core, arguably the singular most important event in core history, however, has remained elusive. In pursuit of this signature we investigate the assumption that the field is a geocentric axial dipole (GAD) over long time scales. We study a suite of numerical dynamo simulations from a paleomagnetic perspective to explore how long the field should be time-averaged to obtain stable paleomagnetic pole directions and intensities. We find that running averages over 20 - 40 kyr are needed to obtain stable paleomagnetic poles with alpha(95) <10 degrees, and over 40 - 120 kyr for alpha(95) <5 degrees, depending on the variability of the field. We find that models with higher heat flux and more frequent polarity reversals require longer time averages, and that obtaining stable intensities requires longer time averaging than obtaining stable directions. Running averages of local field intensity and inclination produce reliable estimates of the underlying dipole moment when reversal frequency is low. However, when heat flux and reversal frequency are increased we find that local observations tend to underestimate virtual dipole moment (VDM) by up to 50% and overestimate virtual axial dipolemoment (VADM) by up to 150%. A latitudinal dependence is found where VDM underestimates the true dipole moment more at low latitudes, while VADM overestimates the true axial dipole moment more at high latitudes. The cause for these observed intensity biases appears to be a contamination of the time averaged field by non-GAD terms, which grows with reversal frequency. We derive a scaling law connecting reversal frequency and site paleolatitude to paleointensity bias (ratio of observed to the true value). Finally we apply this adjustment to the PINT paleointensity record. These biases produce little change to the overall trend of a relatively flat but scattered intensity over the last 3.5 Ga. A more careful intensity adjustment applied during periods when the reversal frequency is known could reveal previously obscured features in the paleointensity record.

In a metal, as in Earth's core, the thermal and electrical conductivities are assumed to be correlated. In a planetary dynamo this implies a contradiction: that both electrical conductivity, which makes it easier to induce current and magnetic field, and conductive heat transport, which hinders thermal convection, should increase simultaneously. Here we show that this contradiction implies that the magnetic induction rate peaks at a particular value of electrical and thermal conductivity and derive the low- and high-conductivity limits for thermal dynamo action. A dynamo regime diagram is derived as a function of electrical conductivity and temperature for Earth's core that identifies four distinct dynamo regimes: no dynamo, thermal dynamo, compositional dynamo, and thermocompositional dynamo. Estimates for the temperature-dependent electrical conductivity of the core imply that the geodynamo may have come close to its high-conductivity "no dynamo" limit prior to inner core nucleation, consistent with recent paleomagnetic observations.

Proxima Centauri b provides an unprecedented opportunity to understand the evolution and nature of terrestrial planets orbiting M dwarfs. Although Proxima Cen b orbits within its star's habitable zone, multiple plausible evolutionary paths could have generated different environments that may or may not be habitable. Here, we use 1-D coupled climate-photochemical models to generate self-consistent atmospheres for several evolutionary scenarios, including high-O-2, high-CO2, and more Earth-like atmospheres, with both oxic and anoxic compositions. We show that these modeled environments can be habitable or uninhabitable at Proxima Cen b's position in the habitable zone. We use radiative transfer models to generate synthetic spectra and thermal phase curves for these simulated environments, and use instrument models to explore our ability to discriminate between possible planetary states. These results are applicable not only to Proxima Cen b but to other terrestrial planets orbiting M dwarfs. Thermal phase curves may provide the first constraint on the existence of an atmosphere. We find that James Webb Space Telescope (JWST) observations longward of 10 mu m could characterize atmospheric heat transport and molecular composition. Detection of ocean glint is unlikely with JWST but may be within the reach of larger-aperture telescopes. Direct imaging spectra may detect O-4 absorption, which is diagnostic of massive water loss and O-2 retention, rather than a photosynthetic biosphere. Similarly, strong CO2 and CO bands at wavelengths shortward of 2.5 mu m would indicate a CO2-dominated atmosphere. If the planet is habitable and volatile-rich, direct imaging will be the best means of detecting habitability. Earth-like planets with microbial biospheres may be identified by the presence of CH4-which has a longer atmospheric lifetime under Proxima Centauri's incident UV-and either photosynthetically produced O-2 or a hydrocarbon haze layer.

We describe a software package called VPLanet that simulates fundamental aspects of planetary system evolution over Gyr timescales, with a focus on investigating habitable worlds. In this initial release, eleven physics modules are included that model internal, atmospheric, rotational, orbital, stellar, and galactic processes. Many of these modules can be coupled to simultaneously simulate the evolution of terrestrial planets, gaseous planets, and stars. The code is validated by reproducing a selection of observations and past results. VPLanet is written in C and designed so that the user can choose the physics modules to apply to an individual object at runtime without recompiling, i.e., a single executable can simulate the diverse phenomena that are relevant to a wide range of planetary and stellar systems. This feature is enabled by matrices and vectors of function pointers that are dynamically allocated and populated based on user input. The speed and modularity of VPLanet enables large parameter sweeps and the versatility to add/remove physical phenomena to assess their importance. VPLanet is publicly available from a repository that contains extensive documentation, numerous examples, Python scripts for plotting and data management, and infrastructure for community input and future development.