Skip to content

Latest commit

 

History

History
73 lines (51 loc) · 4.63 KB

README.md

File metadata and controls

73 lines (51 loc) · 4.63 KB
Raw

Solid Earth carbon degassing and sequestration since 1Ga

Tectono-thermodynamic deep carbon flux modelling workflows from Muller et al. (G3, in press) "Solid Earth carbon degassing and sequestration since 1 billion years ago" based on this plate model: Müller, R.D., Flament, N., Cannon, J., Tetley, M., Williams, S.E., Cao, X., Bodur, Ö.F., Zahirovic, S. and Merdith, A., 2022, A tectonic-rules based mantle reference frame since 1 billion years ago–implications for supercontinent cycles and plate-mantle system evolution, Solid Earth, 13, 1127.

Dependencies

Input files

Run the following notebooks and workflows to prepare the input files required to run these notebooks.

The notebooks will reference these grids using the following directory structure:

    # Change this: Directories to age grids and spreading rate grids
    grid_directory = "./Muller2022_InputGrids/"
    spreadrate_filename = grid_directory+"SpreadingRate/Muller2022_SPREADING_RATE_grid_{:.1f}Ma.nc"
    agegrid_filename = grid_directory+"SeafloorAge/Muller2022_SEAFLOOR_AGE_grid_{:.1f}Ma.nc"

This can be changed to suit the directory made for the input grids above, or the input grids can be saved to the following directory structure for ease of running the notebooks (i.e. they have been designed to run with the following directory structure):

  1. CarbonateSediment - carbonate sediment thickness grids from EarthByte's CarbonateSedimentThickness workflow.
  2. ContinentalMasks - passive margins, trenches and contoured continent polygon gpmls, as well as continent masks made from contoured/buffered continental polygons in EarthByte's continent contouring workflow. These use the Muller et al. (2022) ref frame. The current version is based on the plate model in webDAV: https://www.earthbyte.org/webdav/ftp/Data_Collections/Muller_etal_2022_SE/Muller_etal_2022_SE_Merdith21_PMAG-ref-frame-oceanic-crustal-agegrids_v1.2.1.zip

There is also ContinentalMasksNoArcs, which includes masks made from a manually-edited COB-continent combined file in /utils/COB_polygons_and_coastlines_combined_1000_0_Merdith_etal_no_arcs.gpml. We need this file to ignore island arcs when calculating continental platform degassing in notebook 4.

  1. CrustalCarbon - grids made using CO2_review_paper/Muller_etal_2022/utils/Min-Mean-Max-Crustal-Carbon.ipynb using the plate model in webDAV/
  2. SeafloorAge - grids made with /utils/A-SeafloorGrids.ipynb from the continent masks in /ContinentalMasks/. Considers a full spreading rate in units of mm/yr.
  3. SpreadingRate - as for SeafloorAge. Considers a full spreading rate in units of mm/yr.
  4. TotalSediment - total sediment thickness grids from EarthByte's predicting sediment thickness workflow.

Notebooks:

Once the input files above have been generated, run the following notebooks in order:

  1. 01-Sources-of-Carbon

  2. 02-Subducted-Carbon

  3. 03-Carbon-Degassing

  4. /H2O_review_paper/01-Sources-of-Water

  5. /H2O_review_paper/02-Subducted-Water

  6. /H2O_review_paper/06-Maps-of-water-storage.ipynb

  7. 04-Carbonate-Platform-Degassing

  8. 05-Atmospheric-Carbon

  9. 06-PlateTectonicTools

Some auxiliary notebooks for animations and plots are in /utils/:

  • 0A-Cumulative-Subducted-Carbon.ipynb
  • 0B-Miscellaneous-Plots.ipynb
  • 0C-Carbon-Panel-Plots.ipynb
  • 0D-Carbon-Panel-Plot-Videos.ipynb