Earth's oxygenated history

We simulated Earth with oxygen concentrations between 1000 times less than the present atmospheric level (PAL), and up to 1.5 PAL.

Using WACCM6, we simulated the Earth with 1000 times less than the present atmospheric level (PAL) of oxygen, all the way up to 1.5 times PAL. We found lower ozone columns compared to previous 1D and 3D modelling, meaning that more ultraviolet radiation than previously anticipated could have reached the surface during the Proterozoic. See the following paper for more details: Cooke G. J., Marsh D. R., Walsh C., Black B. and Lamarque J.-F. 2022. A revised lower estimate of ozone columns during Earth’s oxygenated history. R. Soc. open sci.9211165211165. http://doi.org/10.1098/rsos.211165.

Toxic ozone (O3) in otherwise habitable exoplanets

Lethal surface O3 concentrations are possible on habitable zone exoplanets

O3 is a double-edged sword in Earth's atmosphere—vital for protecting life by filtering harmful ultraviolet (UV) radiation, yet a dangerous pollutant at ground level due to human activities. As a powerful oxidizer, ground-level O3 can be harmful (and in some cases lethal) to many organisms.

We investigated two exoplanets: TRAPPIST-1e and Proxima Centauri b, using the advanced Whole Atmosphere Community Climate Model version 6 (WACCM6). Through a series of twelve simulations for TRAPPIST-1 e and four for Proxima Centauri b, we explored various O2 concentrations. Our findings reveal that certain conditions lead to surface O3 levels far exceeding a 40 ppbv threshold which is deemed harmful to some lifeforms on Earth (e.g. vegetation). In fact, some scenarios reach up to 2200 ppbv surface miximg ratios. This would be fatal to most life on Earth, including humans.

However, not all scenarios are bleak. Some simulations show that surface O3 remains below harmful levels across significant areas, though some regions may still pose risks. This is an extra dimension that should be considered when evaluating the habitability of terrestrial exoplanets.

Tidally locked exoplanets: TRAPPIST-1e case study

Degenerate interpretations of O3 spectral features in exoplanet atmosphere observations due to stellar UV uncertainties: a 3D case study with TRAPPIST-1e

We simulated 10 different planetary scenarios for TRAPPIST-1e, a potentially habitable exoplanet: 5 cases with one incoming stellar spectrum, and 5 cases with another. One of these was from Peacock et al. (2019), and the other from Wilson et al. (2021). The latter has much lower ultraviolet (UV) emission than the former. UV emission is important for atmospheric chemistry and the habitability of the surface. Ultimately, the quiescent UV flux is unknown, and TRAPPIST-1 is an active star, so the incident flux will likely vary anyway. This is a general problem, one not isolated to the TRAPPIST-1 system. As we showed in our paper, it can have imporatnt consequences for observational ambiguities when the UV flux is not well-constrained.

More incident UV radiation allows increased photolysis of atmospheric constituents, including O2. O2 photolysis produces two O atoms. O3 is then made through the three body reaction O2 + O + M --> O3 + M, where M is any third body (usually O2 or N2 on Earth). Thus, the production rate of O3 depends not only on the density of O and O2, but on the density of the atmosphere as well. The final predicted amount of atmospheric O3 depends on the interplay between the three-dimensional nature of the planet (transport and no nightside photochemistry), the incoming UV radiation, and the O2 mixing ratio. With more UV, the peak O2 photolysis is closer to the ground. With less O2, UV can also penetrate further. The closer to the surface that O2 photolysis occurs, the greater the O3 production rate, provided that there is enough O2 and O to not limit the reaction rate.

This can mean some weird results occur for the predicted transmission spectra. The W21 pre-industrial (PI) case has undetectable O3 at visible wavelengths (0.6 µm), the lowest height feature at infrared wavelengths, but is actually the most detectable at UV wavelengths (0.25 µm). There are also some cases which overlap in the infrared: W21 1% PAL and P19 PI; W21 0.1% PAL and P19 1% PAL; and W21 10% PAL and P19 10% PAL.

Predicting future exoplanet observations

What would these exopanets look like through current and next generation telescopes?

Cooke et al. (2023), Variability due to climate and chemistry in observations of oxygenated Earth-analogue exoplanets has been published in MNRAS. We investigated the annual and seasonal variability in observations of oxygenated Earth-like exoplanets with next generation telescope concepts (LUVOIR and HabEx). Large variations in planetary brightness occur. Depending on the telescope used and its coronagraph, confirming annual and/or seasonal variability may be possible. We found that clouds, chemistry and surface albedo are all important when considering temporal variability. 3D chemistry climate models are crucial tools for predicting future observations and the detection possibilities that future missions may yield.

Code developments

Open-source python code in Jupyter Notebook on GitHub.

I have developed two Jupyter Notebooks which are both ongoing projects.

One is the Stellar Wind and Irradiance Module (SWIM). This code takes a spectra from the Mega-MUSCLES survey and scales that spectra to a particular exoplanet around the star you have chosen. It also rebins the spectra based on your model of choice. It is easy to use and the user only needs to clone the GitHub repository. All options are given in drop down menus.

The other is a Jupyter Notebook that uses drop down menus to select different types of climate plots from 3D climate data. The user can then plot different variables such as surface sea ice fraction, the atmospheric temperature, and the atmospheric ozone mixing ratio.