Silicate Weathering Model
The package of the silicate weathering model codes include data and code to calculate the modern continental silicate weathering flux and select the best parameters of the model.
Reference
Zuo, Haoyue & Liu, Yonggang & Li, Gaojun & Xu, Zhifang & Zhao, Liang & Guo, Zhengtang & Hu, Yongyun. (2023). A Revised Model of Global Silicate Weathering Considering the Influence of Vegetation Cover on Erosion Rate. 10.5194/gmd-2023-199.
Package Download
https://doi.org/10.5281/zenodo.8423769
About the codes:
This package includes two python files and several data used in the python code
Python file:
1.flux_for_different_runoff_temp_slope_v3_sta.py
This is the file for the first part described in Experiments that can test several different climate and slope datasets to the model
--In this file, we use two different observation dataset to evaluate the model and "sta" variable is the one control the choice
--For the erosion rate revised case, we use "B" and "Be" variable to control and "B" is the erosion rate from TSS and "Be" is the erosion rate from both TSS and Be isotope.
2.new_test_parameter_update_Ec_LAI_global_Etotal.py
This is the file that we redo the parameter selection process and we use 72 cores in this python code
--It will output the .csv result of all the tests which include the basinal river flux and total silicate weathering flux.
--In this code, we don't exclude the possibility that the Ca2++Mg2+ concentration of sediment is higher than metamorphic but when we plot will exclude it.
Some data include: amazon_basin_results.csv, amazon_basins_360_720_global.nc, G1999_basins_360_720.nc, land_area_360_720_PD_match_clim_slope_lith.nc, modern_runoff_v1_interp.nc, modern_runoff_v2_interp.nc, modern_slope, modern_temp_interp.nc
are from Park et al. 2020.
Lithology file: glim_wgs84_0point5deg is from https://doi.pangaea.de/10.1594/PANGAEA.788537
