Multiple hypotheses/models have been put forward regarding the cooling history of the Earth. The search for life beyond Earth has brought these models into a new light as they connect to one of the two energy sources life can tap. The ability to discriminate between different Earth cooling models, and the utility of adopting such models to aid in the assessment of planetary habitability, has been hampered by a lack of uncertainty analysis. This motivates a layered uncertainty analysis for a range of thermal history models that have been applied to the Earth. The analysis evaluates coupled model input, initial condition, and structural uncertainty. Layered model uncertainty, together with data uncertainty and multiple working hypotheses (another form of uncertainty), means that results must be evaluated in a probabilistic sense even if the models are deterministic. For the Earth's cooling history uncertainty leads to ambiguity - multiple models, based on different hypotheses, can match data constraints. This has implications for using such models to forecast conditions for exoplanets that share Earth characteristics but are older than the Earth, i.e., it has implications for modeling the long-term life potential of terrestrial planets. Even for the most Earth-like planet we know of, the Earth itself, model uncertainty and ambiguity leads to large forecast spreads. Given that this comes from the planet with the most data constraints we should expect larger spreads for models of terrestrial planets in general. The layered uncertainty approach can be expanded by coupling planetary cooling models to climate models and propagating uncertainty between them to assess habitability from a probabilistic versus a binary view.

The search for an inhabited planet, other than our own, is a driver of planetary exploration in our solar system and beyond. Using information from our own planet to inform search strategies allows for a targeted search. It is, however, worth considering some span in the strategy and in a priori expectation. An inhabited Earth-like planet is one that would be similar to Earth in ways that extend beyond having biota. To facilitate analysis, we introduce a metric that extends from zero, for an inhabited planet that is like Earth in all other regards (i.e., zero differences), toward positive or negative values for planets that differ from Earth. The analysis shows how assessment of life potential in our galaxy changes more significantly if we find an inhabited planet that is less Earth-like (i.e., it quantifies how probability assessments improve with deviations from Earth-likeness). Discovering such planets could also provide a test of the strong form of the Gaia hypothesis - a test that has proved difficult using only the Earth as a laboratory. Lastly, we discuss how an Earth2.0 narrative, that has been presented to the public as a search strategy, comes with nostalgia-laden philosophical baggage that does not best serve exploration.

Thermal history models, that have been used to understand the geological history of Earth, are now being coupled to climate models to map conditions that allow planets to maintain surface water over geologic time - a criteria considered crucial for life. However, the lack of intrinsic uncertainty assessment has blurred guidelines for how thermal history models can be used toward this end. A model, as a representation of something real, is not expected to be complete. Unmodeled effects are assumed to be small enough that the model maintains utility for the issue(s) it was designed to address. The degree to which this holds depends on how unmodeled factors affect the certainty of model predictions. We quantify this intrinsic uncertainty for several parameterized thermal history models (a widely used subclass of planetary models). Single perturbation analysis is used 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 feedback (negative or positive) within a model. A perturbed physics analysis is then used to determine uncertainty shadows for model outputs. This provides probability distributions for model predictions and tests the structural stability of a model. That is, 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.

Changes that occur on our planet can be tracked back to one of two energy sources: the sun and the Earth's internal energy. The motion of tectonic plates, volcanism, mountain building and the reshaping of our planet's surface over geologic time depend on the Earth's internal energy. Tectonic activity is driven by internal energy and affects the rate at which energy is tapped, i.e., the cooling rate of our planet. Petrologic data indicate that cooling did not occur at a constant rate over geologic history. Interior cooling was mild until ~2.5 billion years ago and then increased (Figure 1). As the Earth cools, it cycles water between its rocky interior (crust and mantle) and its surface. Water affects the viscosity of mantle rock, which affects the pace of tectonics and, by association, Earth cooling. We present suites of thermal-tectonic history models, coupled to deep water cycling, to show that the petrologically constrained change in the Earth's cooling rate can be accounted for by variations in deep water cycling over geologic time. The change in cooling rate does not require a change in the global tectonic mode of the Earth. It can be accounted for by a change in the balance of water cycling between the Earth's interior and its surface envelopes. The nature and timing of that water cycling change can be correlated to a change in the nature of continental crust and an associated rise of atmospheric oxygen. The prediction that the rise of oxygen should then be correlated, in time, to the change in the Earth's cooling rate is consistent with data constraints.