Supplementary Data & Software for "Balancing the marine sulfur cycle", Johnson & Adkins (submitted)

These files are the supplementary data and software for the manuscript "Balancing the marine sulfur cycle" by Johnson & Adkins (submitted).

The MATLAB workspace file "DeepSeaSO4_ClusterAnalysisData.mat" contains organized data structures with DSDP, ODP, and IODP site metadata, porosity data, pore water sulfate concentration data, methane concentration data, total organic carbon wt% data, and calcium carbonate wt% data + their corresponding depths. These data served as the basis for the cluster analysis performed by the included MATLAB script "SO4_conc_cluster_analysis.m" in the study.

The "DatasetS1.xlsx" spreadsheet file contains all compiled pyrite sulfur isotopic data, pyrite abundance data, and total sulfur abundance data included within the literature compilation discussed in the manuscript. Pyrite sulfur isotopic data references are listed and are cited in the main text. DOI URLs for the total sulfur abundance data from the PANGAEA repository are also included in the spreadsheet. These data have additionally been saved into the MATLAB workspace files "Pyrite_d34S_Comp.mat" and "Pyrite Abundance_Comp.mat" for import into MATLAB by other scripts.

The MATLAB workspace file "POROSITY.mat" includes vectors of latitude (in degrees), longitude (in degrees), water depth (in mbsl), sediment depth (in mbsf), porosity (in vol%), and corresponding PANGAEA DOI URLs along with the data tables from which this information was extracted. This file is called by "PorosityPlotter.m" to extract initial porosity values, compile coring methods, and make plots of initial porosity as a function of water depth for each coring method.

The MATLAB script "hypsometry.m" loads the global topographic data included in "topography.mat", calculates associated surface areas across different ocean depth intervals, and determines average porosity and sedimentary pyrite parameters across those intervals in tandem with "S_Iso_Mass.m". "S_Iso_Mass.m" uses the averages to calculate the steady-state sulfur isotopic composition of seawater.

The script "phi_waterdepth_analysis.m" loads the cluster analysis data and uses the corresponding porosity data to calculate linear regressions between initial porosity and water depth for use by other scripts; the workspace files "phi_waterdepth_fit.mat" and "phi_waterdepth_trimmed_fit.mat" are saved versions of the regressions used in the manuscript and called by "hypsometry.m" for calculating averages. 

The additional included files "PyriteSAbundanceCompilationSources.pdf" and "PorosityCompilationSources.pdf" files include full citations to the data included in the pyrite S abundance and porosity compilations (respectively).