DOI: 10.5061/dryad.p2ngf1vx3
title: Data for: The rise of New Guinea and the fall of Neogene global temperatures

description: 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.
description: 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"  <br>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)"
description: # 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)"