Simulations from: Dark land surfaces allow for refugia that could support photosynthetic life on the surface of Snowball Earth
Data files
Jun 10, 2026 version files 106.10 GB
-
010co2_020alb.tar.gz
4.79 GB
-
010co2_030alb.tar.gz
4.74 GB
-
010co2_040alb.tar.gz
4.71 GB
-
025co2_030alb.tar.gz
4.77 GB
-
050co2_020alb.tar.gz
4.84 GB
-
050co2_025alb.tar.gz
4.81 GB
-
050co2_030alb.tar.gz
4.79 GB
-
050co2_035alb.tar.gz
4.77 GB
-
050co2_040alb.tar.gz
4.76 GB
-
100co2_030alb.tar.gz
4.81 GB
-
200co2_020alb.tar.gz
4.88 GB
-
200co2_030alb.tar.gz
4.83 GB
-
200co2_040alb.tar.gz
4.80 GB
-
moderncam4b_050p20p40a_port_base_output.tar
35.70 GB
-
port_experiments_output.tar
7.72 GB
-
README.md
9.62 KB
-
snowball_experiment_rundirs.tar
34.81 MB
-
snowball_port_rundirs.tar
7.21 MB
-
snowball_slim_albedo_inputfiles.tar
336.19 MB
-
variable_table.csv
12.98 KB
Abstract
Photosynthetic eukaryotic algae survived the Neoproterozoic Snowball Earth events, indicating that liquid-water refugia existed somewhere on the surface. We examine the potential for refugia at the coldest time of a snowball event, before CO2 had risen and with high-albedo ice on the frozen ocean, before it became darkened by dust deposition. We use the Community Earth System Model to simulate a “modern” Snowball Earth (i.e., with continents in their current configuration), in which the ocean surface has frozen to the equator as “sea glaciers”, hundreds of meters thick, flowing like ice shelves. Despite global mean surface temperatures below -60°C, some areas of the land surface reach above-freezing temperatures because they are darker than the ice-covered ocean. With low CO2 (10 ppm) and land-surface albedo characteristic of bright sand-deserts (0.4), 0.1 percent of the land surface could host liquid water seasonally; this increases to 12 percent for darker land of albedo characteristic of polar deserts (0.2). Given a water source to these locations, such as a narrow bay intruding from the ocean like the modern Red Sea, or summer precipitation or melt of nearby glaciers, we argue that photosynthetic life would have the potential to survive on or near warm, bare land. The abundance of warm land and consequently potential refugia increases more strongly in response to reducing the land albedo than to higher CO2, for the same global radiative forcing due to the concentration of forcing over land.
https://doi.org/10.5061/dryad.5hqbzkhcf
In this repository, we provide input data and output from CESM from a series of model experiments testing the habitability of Snowball Earth.
Description of experiments
Experiments associated with this project sample a range of bare land surface albedo and atmospheric CO2 concentrations as described in the paper. The output data is contained within the tar file for each experiment as described in this table:
| CO2 | Visible albedo/NIR albedo | Broadband Solar Albedo | Output Files tarball name |
|---|---|---|---|
| 10ppm | 0.1/0.3 | 0.2 | 010co2_020alb.tar.gz |
| 10ppm | 0.2/0.4 | 0.3 | 010co2_030alb.tar.gz |
| 10ppm | 0.3/0.5 | 0.4 | 010co2_040alb.tar.gz |
| 25ppm | 0.2/0.4 | 0.3 | 025co2_030alb.tar.gz |
| 50ppm | 0.1/0.3 | 0.2 | 050co2_020alb.tar.gz |
| 50ppm | 0.15/0.35 | 0.25 | 050co2_025alb.tar.gz |
| 50ppm | 0.2/0.4 | 0.3 | 050co2_030alb.tar.gz |
| 50ppm | 0.25/0.45 | 0.35 | 050co2_035alb.tar.gz |
| 50ppm | 0.3/0.5 | 0.40 | 050co2_040alb.tar.gz |
| 100ppm | 0.2/0.4 | 0.30 | 100co2_030alb.tar.gz |
| 200ppm | 0.1/0.3 | 0.2 | 200co2_020alb.tar.gz |
| 200ppm | 0.2/0.4 | 0.3 | 200co2_030alb.tar.gz |
| 200ppm | 0.3/0.5 | 0.4 | 200co2_040alb.tar.gz |
Additional experiments to quantify the radiative forcing of changing albedo and CO2 were completed with a configuration of the model that runs only the radiation code called "port". Output from the following simulations is included:
| co2 | Visible/NIR albedo | experiment name |
|---|---|---|
| base case: 50ppm | base case: 0.2/0.4 | moderncam4b_050p20p40a_port_base |
| 50ppm | 0.2/0.4 | moderncam4b_050p20p40a_port_p20p40alb_case |
| 10ppm | 0.2/0.4 | moderncam4b_050p20p40a_port_0.2xco2_case |
| 100ppm | 0.2/0.4 | moderncam4b_050p20p40a_port_2.0xco2_case |
| 200ppm | 0.2/0.4 | moderncam4b_050p20p40a_port_4.0xco2_case |
| 50ppm | 0.1/0.3 | moderncam4b_050p20p40a_port_p10p30alb_case |
| 50ppm | 0.15/0.35 | moderncam4b_050p20p40a_port_p15p35alb_case |
| 50ppm | 0.25/0.45 | moderncam4b_050p20p40a_port_p25p45alb_case |
| 50ppm | 0.3/0.5 | moderncam4b_050p20p40a_port_p30p50alb_case |
The base case is in the tarball:
moderncam4b_050p20p40a_port_base_output.tar
The experiment output for all experiments is in the tarball:
port_experiments_output.tar
Description of the data and file structure
This dataset is organized into zipped tarballs for each simulation as well as tarballs of the case directories and input files used to generate these outputs.
Atmosphere data is provided for all experiments, and land and ice data is provided for a select set of experiments (warmest and coldest) to contextualize the ice and land. Solar and land input data files are also provided.
Variables
Output from the model simulations is stored in NetCDF files. The model output for the full model experiments is given in a time series of monthly means for each variable. For the port radiation runs, there is one file for each month that contains all variables.
The following variables are analyzed in the Shum et al. publication
| Variable Name | Model Domain | Variable Long Name | Unit |
|---|---|---|---|
| flat_ai | ice | latent heat flux | W/m2 |
| FLNR | port | Net longwave flux at tropopause | W/m2 |
| FSNR | port | Net solar flux at tropopause | W/m2 |
| FSUTOA | atm | Upwelling solar flux at top of atmosphere | W/m2 |
| LANDFRAC | atm | Fraction of sfc area covered by land | fraction |
| LHFLX | atm | Surface latent heat flux | W/m2 |
| MML_snow | lnd | MML total amount of snow in snow bucket (hydrology) | kg/m2 |
| PRECC | atm | Convective precipitation rate (liq+ice) | m/s |
| PRECL | atm | Large-scale (stable) precipitation rate (liq+ice) | m/s |
| PS | atm | Surface pressure | Pa |
| QREFHT | atm | Reference height humidity | kg/kg |
| RELHUM | atm | Relative humidity | percent |
| SHFLX | atm | Surface sensible heat flux | W/m2 |
| SNOWHLND | atm | Water equivalent snow depth | m |
| SOLIN | atm | Solar insolation | W/m2 |
| T | atm | Temperature | K |
| TREFHT | atm | Reference height temperature | K |
| TS | atm | Surface temperature (radiative) | K |
| TSMX | atm | Maximum surface temperature over output period | K |
| U10 | atm | 10m wind speed | m/s |
| UU | atm | Zonal velocity squared | m2/s2 |
| vsnon | ice | snow depth on ice categories | m |
| VV | atm | Meridional velocity squared | m2/s2 |
The relevant metadata for the variables is contained in the netCDF files and can also be found in the included file
varibles_table.csv
Further information about CESM variables can be found in the CESM output variable directory.
Sharing/Access information
Data for this project is open-access and available from this repository.
We acknowledge high-performance computing support from Cheyenne (https://doi.org/10.5065/D6RX99HX) provided by NCAR’s Computational and Information Systems Laboratory, sponsored by the National Science Foundation.
Code/Software
Analysis code to create the figures and calculations in the manuscript associated with this dataset is archived on Zenodo (https://doi.org/10.5281/zenodo.18247372).
NetCDF files can be opened and viewed using the xarray python package.
Source code for CESM, the model used in this study, is archived and publicly accessible online (https://doi.org/10.5281/zenodo.3895306). We thank all scientists, software engineers, and administration who contributed to CESM2's development. The CESM user guide describes how to run CESM and how different inputs (e.g., namelist settings) are used when running the model. The CESM tutorial is likely helpful for anyone looking to run CESM.
Input files that are custom to the experiments in this paper are included in this tarball:
snowball_slim_albedo_inputfiles.tar
These input files are NetCDF files that specify the surface albedo and are needed to reproduce the model simulations.
We additionally include the case directories for all experiments, which include additional configuration details needed to repeat these experiments in standard CESM formats in two tarballs, one for the full model experiments, and another for the port experiments:
snowball_experiment_rundirs.tar
snowball_port_rundirs.tar
These configuration files are likely only relevant to people already familiar with running the Community Earth System Model (CESM). Within each tarball, there are directories for each simulation (i.e., moderncam4b_200p30p50a). The upper-level files in each experiment directory contain .xml files that specify run configurations. Model namelist files also specify configurations (i.e., user_nl_cam for the atmospheric model subcomponent). The SourceMods directory contains Fortran source files which document changes to the default Community Earth System Model (CESM) code. All the files in each experiment subdirectory can be opened in a text editor.
