Data for: The rise of New Guinea and the fall of Neogene global temperatures
Data files
Sep 01, 2023 version files 960.84 MB
-
erosion_projection.py
-
erosion_transect_map.png
-
gdss_output_NG_slope_0.5-0.0.nc
-
gdss_output_NG_slope_10.0-9.6.nc
-
gdss_output_NG_slope_10.7-10.0.nc
-
gdss_output_NG_slope_11.4-10.7.nc
-
gdss_output_NG_slope_12.0-11.4.nc
-
gdss_output_NG_slope_13.2-12.0.nc
-
gdss_output_NG_slope_14.4-13.2.nc
-
gdss_output_NG_slope_15.6-14.4.nc
-
gdss_output_NG_slope_3.0-0.5.nc
-
gdss_output_NG_slope_4.0-3.0.nc
-
gdss_output_NG_slope_4.8-4.0.nc
-
gdss_output_NG_slope_6.1-4.8.nc
-
gdss_output_NG_slope_6.8-6.1.nc
-
gdss_output_NG_slope_7.3-6.8.nc
-
gdss_output_NG_slope_7.9-7.3.nc
-
gdss_output_NG_slope_8.1-7.9.nc
-
gdss_output_NG_slope_8.6-8.1.nc
-
gdss_output_NG_slope_9.0-8.6.nc
-
gdss_output_NG_slope_9.6-9.0.nc
-
geoclim.py
-
make_slope_inputs.py
-
README.md
-
Total_erosion.csv
-
transect_erosion.pdf
-
Transect_location.dbf
-
Transect_location.prj
-
Transect_location.shp
-
Transect_location.shx
Abstract
This dataset contains the output files of numerical simulations of the model GEOCLIM-DynSoil-steady-state, which were conducted for the study "The rise of New Guinea and the fall of Neogene global temperatures". Those outputs are modelled fields of erosion rate, silicate weathering flux, regolith thickness and primary phase depletion, on the island of New Guinea, at 30-minute resolution, and for 19 time slices from 15 Ma to present-day. This dataset also contains additional Python scripts for generating inputs of the model (namely, slope field) and for projecting the 2D New Guinea erosion field on the 1D transect highlighted in the abovementioned study.
README: Data for the study "The rise of New Guinea and the fall of Neogene global temperatures"
This dataset contains the output files of numerical simulations of the model GEOCLIM-DynSoil-steady-state,
that were conducted for the study "The rise of New Guinea and the fall of Neogene global temperatures".
Those outputs are modeled fields of erosion rate, silicate weathering flux, regolith thickness and primary phases depletion,
on the island of New Guinea, at 30 minutes resolution, and for 19 time slices from 15 Ma to present-day.
This dataset also contains additional Python scripts for generating inputs of the model (namely, slope field)
and for projecting the 2D New Guinea erosion field on the 1D transect highlighted in the above-mentioned study.
files description
- The 19 files "gdss_output_NG_slope_X.X-X.X.nc" are output files from the model GEOCLIM-DynSoil-steady-state. "X.X-X.X" indicates the time slice targeted by the simulation (e.g., "14.4-13.2" is "14.4 Ma -- 13.2 Ma). Metadata in each file (netCDF format) describes the fields computed by the model and outputted in the files. Instruction to reproduce those outputs are given in the corresponding branch of the model's Github repository (https://doi.org/10.5281/zenodo.8245945)
- "Total_erosion.csv" is a 3-column table indicating, for each time interval targeted by the modeling study, the beginning and ending of the time interval, and the total eroded material during that interval, as computed by the model MOVE2D, on the New Guinea 1D transect (see "Transect_locations" files) Because the transect is 1D, the "total eroded material" during a time interval is a 2D value, and its units is km^2.
- "make_slope_inputs.py" is the Python script that generates the slope field for each time interval, accordingly to the erosion rate during that interval. It uses, as input, a modern (0.5-0Ma) slope field, that can be found in the model's Github repository https://doi.org/10.5281/zenodo.8245945.
- "erosion_projection.py": script to project the erosion field computed by GEOCLIM-DynSoil-steady-state on the 1D transect, and generate the 2 figures (already present in the current repertory):
- "erosion_transect_map.png": map the New Guinea erosion with the transect location.
- "transect_erosion.pdf": plot of the erosion rate along the transect.
- "geoclim.py": Python script storing the erosion function taken from the GEOCLIM model. This function may be used by "erosion_projection.py", as an alternative way to recompute erosion "from scratch".
- "Transect_location.*": shapefile indicating the location of 1D New Guinea transect (line).
External references
- Github repository of the model code: https://github.com/piermafrost/GEOCLIM-dynsoil-steady-state, release "v2-NG". The release associated with the study "The rise of New Guinea..." is archived on Zenodo: https://doi.org/10.5281/zenodo.8245945
- Associated article: "Peter E. Martin, Francis A. Macdonald, Nadine McQuarrie, Rebecca M. Flowers, and Pierre J.Y. Maffre: The rise of New Guinea and the fall of Neogene global temperatures. Proceedings of the National Academy of Sciences, (accepted in August 2023)"
Usage notes
Dataset structure
- "GDSS\_outputs/": the raw outputs from GEOCLIM-DynSoil-steady-state simulations, 1 file per time slice. NetCDF format (including self-describing metadata).
- "slope\_construction/": scripts generating the slope field (input of GEOCLIM-DynSoil-steady-state model) for each time slice.
- "transect\_erosion\_projection/": scripts to project the erosion field computed by GEOCLIM-DynSoil-steady-state on the 1D transect, and generate figures.
External references
- Github repository of the model code: https://github.com/piermafrost/GEOCLIM-dynsoil-steady-state, release "v2-NG"
The release associated with the study "The rise of New Guinea..." is archived on Zenodo: https://doi.org/10.5281/zenodo.8245945 - Associated article: "Peter E. Martin, Francis A. Macdonald, Nadine McQuarrie, Rebecca M. Flowers, and Pierre J.Y. Maffre: The rise of New Guinea and the fall of Neogene global temperatures. Proceedings of the National Academy of Sciences, (accepted in August 2023)"