Skip to main content
Dryad logo

Varied oxygen simulations with WACCM6 (Proterozoic to pre-industrial atmosphere)


Cooke, Gregory et al. (2021), Varied oxygen simulations with WACCM6 (Proterozoic to pre-industrial atmosphere), Dryad, Dataset,


The history of molecular oxygen (O2) in Earth’s atmosphere is still debated; however, geological evidence supports at least two major episodes where O2 increased by an order of magnitude or more: the Great Oxidation Event (GOE) and the Neoproterozoic Oxidation Event. Oconcentrations have likely fluctuated (between 10−3 and 1.5 times the present atmospheric level) since the GOE ∼ 2.4 Gyr ago, resulting in a time-varying ozone (O3) layer. Using a three-dimensional (3D) chemistry climate model, we simulate changes in O3 in Earth’s atmosphere since the GOE and consider the implications for surface habitability, and glaciation during the Mesoproterozoic. We find lower O3 columns (reduced by up to 4.68 times for a given O2 level) compared to previous work; hence, higher fluxes of biologically harmful UV radiation would have reached the surface. Reduced O3 leads to enhanced tropospheric production of the hydroxyl radical (OH) which then substantially reduces the lifetime of methane (CH4). We show that a CH4 supported greenhouse effect during the Mesoproterozoic is highly unlikely. The reduced O3 columns we simulate have important implications for astrobiological and terrestrial habitability, demonstrating the relevance of 3D chemistry-climate simulations when assessing paleoclimates and the habitability of faraway worlds.


The data was originally produced using the Earth System Model WACCM6. The corresponding paper that describes WACCM can be found at the following DOI: WACCM6 is a model configuration of the Community Earth System Model version 2 (CESM2). We used CESM2.1.3, which can be downloaded from the following URL: 

The atmospheric data from each simulation was output in terms of monthly means (h0 files; h for history), and 5 day averages (h2 files). For both h0 and h2 files, the last 4 years of each simulation was averaged using netCDF operators. ncra for time average (h0 and h2 files), and ncwa -alon for zonal mean (h2 files). 

The data has been processed using the Python programming language. Two .py files are included named and The python files may not be well-commented in all places.

To produce the solar file which represent the Sun 2 billion years ago (2 Gyr), we used an existing solar model by Claire et al. 2012 ( The model can be downloaded here:

The datasets Kasting_Catling_Fig10_data.csv and ROCKE3D_Fig3_O3_col.csv are from the following papers: and, respectively. The data was recovered using WebPlotDigitizer (

Usage Notes

To use the data:

The data can be downloaded and processed using Python or a programming language that can handle netCDF data. The python (.py) files included were used by the authors for producing figures and analysis for the corresponding research article to this dataset. To use the python files, you need to correctly set up the paths to the correct directories and files in and h0 files should be in path + 'Ozone_CH4_paper/h0/', h2 files in path + 'Ozone_CH4_paper/h2/', and solar files in path + 'Ozone_CH4_paper/Solar_files/'. The path can be any directory you have permissions for. For example, on a Linux system, using /localhome/user/ is how the path was set up for the python file when used by the authors. To run the python code, run the file. If you have not used python before, downloading python via Anaconda and installing xarray should allow you to see the code (see here:, and here: Figures produced will be deposited into: path + 'Ozone_CH4_paper/ 

The h0 and h2 files correspond to the following simulations, which are described in Table 1 of the associated manuscript:

150% PAL:;  
50% PAL:;
10% PAL:;  
5% PAL:;
1% PAL:;
CH4 em1:
CH4 em0.1:
YS 4xCO2:
0.5% PAL:;
0.1% PAL:;

In addition, there are four files for calculating the three major loss channels of methane (see Fig. 8 in the associated article at
10% PAL:
1% PAL:
0.1% PAL:

Note that b.e21.BWma1850 corresponds to a WACM6 (CESM2.1) pre-industrial (year 1850) simulation. f19_g17 describes the grid setup. 96 x 144 latitude x longitude, corresponding to 1.875 degrees x 2.5 degrees. The model years of the dataset are given by 0014-0017, for years 14, 15, 16, and 17. 'zm' means a zonal mean was applied to the file.

A file is provided which contains gaussian weights for the area of grid cells (gaussian weights vary with latitude):

Solar files are as follows:

Present day Sun:
Sun's spectrum 2 Gyr ago (2.0 Ga):


Science and Technology Facilities Council, Award: ST/T506230/1

National Center for Atmospheric Research

Science and Technology Facilities Council, Award: MR/T040726/1

Science and Technology Facilities Council, Award: ST/T000287/1

Science and Technology Facilities Council, Award: ST/R000549/1