Notebooks for computing cumulative subducted mantle carbon based on the Muller et al. (2025) plate model.
The following Python libraries are required to run this workflow:
The required input files for this workflow are:
- All files in
\utils\, including the filesfft_gaussian_filter.py,C_serp_array_individual.npz,serpentinite_spreading_grids.pyandwater_lithosphere_storage.py. - Input grids, which should follow the structure:
MantleCarbon/
│
├── 01-Sources-of-Carbon.ipynb
├── 02-Subducted-Carbon.ipynb
├── 03-Cumulative-Subducted-Carbon.ipynb
│
├── utils/
│ └── ...
│
└── Grids/
│
└── InputGrids/
├── CarbonateSediment/
├── ContinentalMasks/
├── CrustalCarbon/
├── SeafloorAge/
├── SpreadingRate/
└── TotalSediment/
The input grids should be prepared as follows, depending on the plate reconstruction model required, see EarthByte Group's list of available global and regional plate motion models:
CarbonateSedimentfrom the GitHub workflow CarbonateSedimentThicknessContinentalMasks,SeafloorAgeandSpreadingRatemay be found in the EarthByte Group's list of available global and regional plate motion models, subject to availability. Otherwise, the grids can be produced fromgplately'sSeafloorGridmodule - see the sample notebook here. Otherwise, grids may also be available for the selected plate model elsewhere, like Zenodo or webDAV.CrustalCarbonfrom\utils\, which itself requires theSeafloorAgeandSpreadingRategrids.TotalSedimentfrom the GitHub workflow. More details here.
- 01-Sources-of-Carbon.ipynb - Produce grids for total contained carbon for each reservoir.
- 02-Subducted-Carbon.ipynb - Interpolate the grids from (1) at trenches to produce subducted carbon per reservoir.
- 03-Cumulative-Subducted-Carbon.ipynb - Smooth the grids from (2) and accumulate them per timesteps to develop cumulative subducted carbon.