2020ApJ...893..114S

Multiple hypotheses/models have been put forward regarding Earth's cooling history. Searching for life beyond Earth has brought these models into a new light as they connect to an energy source that life can tap. Discriminating between different cooling models and adapting them to aid in the assessment of planetary habitability has been hampered by a lack of uncertainty quantification. Here, we provide an uncertainty quantification that accounts for a range of interconnected model uncertainties. This involved calculating over a million individual model evolutions to determine uncertainty metrics. Accounting for uncertainties means that model results must be evaluated in a probabilistic sense, even though the underlying models are deterministic. The uncertainty analysis was used to quantify the degree to which different models can satisfy observational constraints on the Earth's cooling. For the Earth's cooling history, uncertainty leads to ambiguity-multiple models, based on different hypotheses, can match observations. This has implications for using such models to forecast conditions for exoplanets that share Earth characteristics but are older than the Earth, i.e., ambiguity has implications for modeling the long-term life potential of terrestrial planets. Even for the most earthlike planet we know of, the Earth itself, model uncertainty and ambiguity leads to large forecast spreads. Given that Earth has the best data constraints, we should expect larger spreads for models of terrestrial planets, in general. The uncertainty analysis provided here can be expanded by coupling planetary cooling models to climate models and propagating uncertainty between them to assess habitability from a probabilistic view.