An official website of the United States government
Data for: 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).
more »
« less
- Award ID(s):
- 1925990
- PAR ID:
- 10560802
- Publisher / Repository:
- Dryad
- Date Published:
- Subject(s) / Keyword(s):
- FOS: Earth and related environmental sciences FOS: Earth and related environmental sciences Silicate weathering ophiolite Miocene
- Format(s):
- Medium: X Size: 961006094 bytes
- Size(s):
- 961006094 bytes
- Right(s):
- Creative Commons Zero v1.0 Universal
- Sponsoring Org:
- National Science Foundation
More Like this
-
The ~2,000-km-long Central Range of New Guinea is a hotspot of modern carbon sequestration due to the chemical weathering of igneous rocks with steep topography in the warm wet tropics. These high mountains formed in a collision between the Australian plate and ophiolite-bearing volcanic arc terranes, but poor resolution of the uplift and exhumation history has precluded assessments of the impact on global climate change. Here, we develop a palinspastic reconstruction of the Central Range orogen with existing surface geological constraints and seismic data to generate time–temperature paths and estimate volumes of eroded material. New (U-Th)/He thermochronology data reveal rapid uplift and regional denudation between 10 and 6 Mya. Erosion fluxes from the palinspastic reconstruction, calibrated for time with the thermochronological data, were used as input to a coupled global climate and weathering model. This model estimates 0.6 to 1.2 °C of cooling associated with the Late Miocene rise of New Guinea due to increased silicate weathering alone, and this CO2sink continues to the present. Our data and modeling experiments support the hypothesis that tropical arc-continent collision and the rise of New Guinea contributed to Neogene cooling due to increased silicate weathering.more » « less
-
GEOCLIM model as it was used for the experiments presented in Maffre, P., Swanson-Hysell, N.L., and Godd\'eris, Y. (2021), Limited carbon cycle response to increased sulfide weathering due to oxygen feedback, ]Geophysical Research Letters, doi:10.1029/2021GL094589. GEOCLIM is cluster of models, more or less adaptable, computing geochemical cycles of several species (C, O...) at geological timescale. The core of the model is an ocean-atmosphere chemistry model (advection-reaction) COMBINE. Ocean an atmosphere are discretized in 10 reservoirs (boxes). This core is closely associated to an early diagenesis module computing the "output" fluxes (burial of elements in marine sediments) for each box. It is also associated to an continental weathering module computing the "input" fluxes. This weathering module is spatially-resolved (using a geographic mesh grid), its resolution is adaptable, and several options exist for the silicate weathering part. This triplet is indirectly coupled to a climate model (GCM). Climate simulations must be run before using GEOCLIM, for a range of CO2 levels. Any climate model can be used, as long as it computes surface air temperature and continental runoff. Oceanic temperature from the climate model can be used, or parameterized if not available. Climate fields are then interpolated on the "CO2 dimension", at the current atmospheric CO2 computed by GEOCLIM, and used to compute continental weathering, and oceanic boxes temperature. The resolution (ie, the spatial grid) of the continental weathering module must be the same than the GCM. How to run the model: After downloading the present repository, type ./make_test testname (testname being one of "ERA5", "GFDL", "CESM", "paleo" and "ascii"). This command will compile and execute a short GEOCLIM run and compare the output to a reference template. This allows to verify that the compilation and execution of the model are performed without error, and yield the same results than reference runs. If not, the command should tell what type of error was encountered (see also section Frequent issues at the end of this file). You could try the 5 tests to make sure everything works as excepted. If the tests are conclusive, follow those step to create your run: Compile the code with build_GEOCLIM (specifying the model set of components and resolution). Try ./build_GEOCLIM --help for more information. This command create an executable file in 'executable/' Configure the run by editing the files 'config/IO_CONDITIONS' and 'config/cond_p20.dat' (name of the run, initial condition, forcing fields, solver parameters...). This post-compilation configuration must be consistent with the specified set of components and resolution. Run the model with executable/exec_name (name of the executable created at the compilation step). Alternatively, and if you are running the model on a cluster, you can use the files in 'job/' to submit the run as a batch process. 'job/submit_chain.sh' is a script designed for submitting a series of runs, each new one starting from the end of the previous one. See also section Multiple runs and job submission in Run GEOCLIM.more » « less
-
Soil piping (concentrated leak erosion) is a major contributor to soil erosion in many parts of the world, and collapse of eroded pipes can result in the formation of gullies and sinkholes or trigger slope instability. Despite these significant impacts, there is little understanding of factors controlling pipe collapse, and how water within the pipe influences moisture levels within a slope. In this study, physical models were employed on unsaturated model slopes with pre-formed macropores to investigate how soil properties, pipe characteristics, and hydraulic conditions govern internal erosion processes and slope stability. Experiments simulated shallow field conditions (0.45 m overburden) using 4 mm and 12 mm pipes to establish preferential flow paths, while varying model parameters including initial compaction moisture content and density, pipe condition (absent, closed, or open), slope angle, and model width. Volumetric water content sensors monitored moisture evolution, while cameras captured slope responses to subsurface flow. Results demonstrate that initial compaction conditions (water content and density), pipe size, hydraulic connectivity, and pipe condition control internal erosion processes and slope stability.more » « less
-
Mathelier, Anthony (Ed.)Abstract Motivation Methods to model dynamic changes in gene expression at a genome-wide level are not currently sufficient for large (temporally rich or single-cell) datasets. Variational autoencoders offer means to characterize large datasets and have been used effectively to characterize features of single-cell datasets. Here, we extend these methods for use with gene expression time series data. Results We present RVAgene: a recurrent variational autoencoder to model gene expression dynamics. RVAgene learns to accurately and efficiently reconstruct temporal gene profiles. It also learns a low dimensional representation of the data via a recurrent encoder network that can be used for biological feature discovery, and from which we can generate new gene expression data by sampling the latent space. We test RVAgene on simulated and real biological datasets, including embryonic stem cell differentiation and kidney injury response dynamics. In all cases, RVAgene accurately reconstructed complex gene expression temporal profiles. Via cross validation, we show that a low-error latent space representation can be learnt using only a fraction of the data. Through clustering and gene ontology term enrichment analysis on the latent space, we demonstrate the potential of RVAgene for unsupervised discovery. In particular, RVAgene identifies new programs of shared gene regulation of Lox family genes in response to kidney injury. Availability and implementation All datasets analyzed in this manuscript are publicly available and have been published previously. RVAgene is available in Python, at GitHub: https://github.com/maclean-lab/RVAgene; Zenodo archive: http://doi.org/10.5281/zenodo.4271097. Supplementary information Supplementary data are available at Bioinformatics online.more » « less
