Geyman, Emily C.; Wu, Ziman; Nadeau, Matthew D.; Edmonsond, Stacey; Turner, Andrew; Purkis, Sam J.; Howes, Bolton; Dyer, Blake; Ahm, Anne-Sofie C.; Yao, Nan; Deutsch, Curtis A.; Higgins, John A.; Stolper, Daniel A.; Maloof, Adam C.
Carbonate mud represents one of the most important geochemical archives for reconstructing ancient climatic, environmental, and evolutionary change from the rock record. Mud also represents a major sink in the global carbon cycle. Yet, there remains no consensus about how and where carbonate mud is formed. In this contribution, we present new geochemical data that bear on this problem, including stable isotope and minor and trace element data from carbonate sources in the modern Bahamas such as ooids, corals, foraminifera, and green algae.
Kim, Donghoon; Tracy, Sally J.; Smith, Raymond F.; Gleason, Arianna E.; Bolme, Cindy A.; Prakapenka, Vitali B.; Appel, Karen; Speziable, Sergio; Wicks, June K.; Berryman, Eleanor J.; Han, Sirus K.; Schoelmerich, Markus O.; Lee, Hae Ja; Nagler, Bob; Cunningham, Eric F.; Akin, Minta C.; Asimow, Paul D.; Eggert, Jon H.; Duffy, Thomas S.
The behavior of forsterite, Mg2SiO4, under dynamic compression is of fundamental importance for understanding its phase transformations and high-pressure behavior. Here, we have carried out an in situ X-ray diffraction study of laser-shocked poly- and single-crystal forsterite (a-, b-, and c- orientations) from 19 to 122 GPa using the Matter in Extreme Conditions end-station of the Linac Coherent Light Source. Under laser-based shock loading, forsterite does not transform to the high-pressure equilibrium assemblage of MgSiO3 bridgmanite and MgO periclase, as was suggested previously. Instead, we observe forsterite and forsterite III, a metastable polymorph of Mg2SiO4, coexisting in a mixed-phase region from 33 to 75 GPa for both polycrystalline and single-crystal samples. Densities inferred from X-ray diffraction data are consistent with earlier gas-gun shock data. At higher stress, the behavior observed is sample-dependent. Polycrystalline samples undergo amorphization above 79 GPa. For - and -oriented crystals, a mixture of crystalline and amorphous material is observed to 108 GPa, whereas the -oriented crystal adopts an unknown crystal structure at 122 GPa. The Q values of the first two sharp diffraction peaks of amorphous Mg2SiO4 show a similar trend with compression as those observed for MgSiO3 glass in both recent static and laser-compression experiments. Upon release to ambient pressure, all samples retain or revert to forsterite with evidence for amorphous material also present in some cases. This study demonstrates the utility of femtosecond free-electron laser X-ray sources for probing the time evolution of high-pressure silicates through the nanosecond-scale events of shock compression and release.
Berryman, Eleanor J.; Winey, J. M.; Gupta, Yogendra M.; Duffy, Thomas S.
Stishovite (rutile-type SiO2) is the archetype of dense silicates and may occur in post-garnet eclogitic rocks at lower-mantle conditions. Sound velocities in stishovite are fundamental to understanding its mechanical and thermodynamic behavior at high pressure and temperature. Here, we use plate-impact experiments combined with velocity interferometry to determine the stress, density, and longitudinal sound speed in stishovite formed during shock compression of fused silica at 44 GPa and above. The measured sound speeds range from 12.3(8) km/s at 43.8(8) GPa to 9.8(4) km/s at 72.7(11) GPa. The decrease observed at 64 GPa reacts a decrease in the shear modulus of stishovite, likely due to the onset of melting. By 72 GPa, the measured sound speed agrees with the theoretical bulk sound speed indicating loss of all shear stiffness due to complete melting. Our sound velocity results provide direct evidence for shock-induced melting, in agreement with previous pyrometry data.
The history of organismal evolution, seawater chemistry, and paleoclimate is recorded in layers of carbonate sedimentary rock. Meter-scale cyclic stacking patterns in these carbonates often are interpreted as representing sea level change. A reliable sedimentary proxy for eustasy would be profoundly useful for reconstructing paleoclimate, since sea level responds to changes in temperature and ice volume. However, the translation from water depth to carbonate layering has proven difficult, with recent surveys of modern shallow water platforms revealing little correlation between carbonate facies (i.e., grain size, sedimentary bed forms, ecology) and water depth. We train a convolutional neural network with satellite imagery and new field observations from a 3,000 km2 region northwest of Andros Island (Bahamas) to generate a facies map with 5 m resolution. Leveraging a newly-published bathymetry for the same region, we test the hypothesis that one can extract a signal of water depth change, not simply from individual facies, but from sequences of facies transitions analogous to vertically stacked carbonate strata. Our Hidden Markov Model (HMM) can distinguish relative sea level fall from random variability with ∼90% accuracy. Finally, since shallowing-upward patterns can result from local (autogenic) processes in addition to forced mechanisms such as eustasy, we search for statistical tools to diagnose the presence or absence of external forcings on relative sea level. With a new data-driven forward model that simulates how modern facies mosaics evolve to stack strata, we show how different sea level forcings generate characteristic patterns of cycle thicknesses in shallow carbonates, providing a new tool for quantitative reconstruction of ancient sea level conditions from the geologic record.
Amazonian deforestation causes systematic changes in regional dry season precipitation. Some of these changes at contemporary large scales (a few hundreds of kilometers) of deforestation have been associated with a ‘dynamical mesoscale circulation’, induced by the replacement of rough forest with smooth pasture. In terms of decadal averages, this dynamical mechanism yields increased precipitation in downwind regions and decreased precipitation in upwind regions of deforested areas. Daily, seasonal, and interannual variations in this phenomenon may exist, but have not yet been identified or explained. This study uses observations and numerical simulations to develop relationships between the dynamical mechanism and the local- and continental-scale atmospheric conditions across a range of time scales. It is found that the strength of the dynamical mechanism is primarily controlled by the regional-scale thermal and dynamical conditions of the boundary layer, and not by the continental- and global-scale atmospheric state. Lifting condensation level and wind speed within the boundary layer have large and positive correlations with the strength of the dynamical mechanism. The strength of these relationships depends on time scale and is strongest over the seasonal cycle. Overall, the dynamical mechanism is found to be strongest during times when the atmosphere is relatively stable. Hence, for contemporary large scales of deforestation this phenomenon is found to be the prevalent convective triggering mechanism during the dry and parts of transition seasons (especially during the dry-to-wet transition), significantly affecting the hydroclimate during this period.
Explosive volcanic eruptions have large climate impacts, and can serve as observable tests of the climatic response to radiative forcing. Using a high resolution climate model, we contrast the climate responses to Pinatubo, with symmetric forcing, and those to Santa Maria and Agung, which had meridionally asymmetric forcing. Although Pinatubo had larger global-mean forcing, asymmetric forcing strongly shifts the latitude of tropical rainfall features, leading to larger local precipitation/TC changes. For example, North Atlantic TC activity over is enhanced/reduced by SH-forcing (Agung)/NH-forcing (Santa Maria), but changes little in response to the Pinatubo forcing. Moreover, the transient climate sensitivity estimated from the response to Santa Maria is 20% larger than that from Pinatubo or Agung. This spread in climatic impacts of volcanoes needs to be considered when evaluating the role of volcanoes in global and regional climate, and serves to contextualize the well-observed response to Pinatubo.
Khanna, Jaya; Medvigy, David; Fueglistaler, Stephan; Walko, Robert
More than 20% Amazon rainforest has been cleared in the past three decades triggering important hydroclimatic changes. Small-scale (~few kilometers) deforestation in the 1980s has caused thermally-triggered atmospheric circulations that increase regional cloudiness and precipitation frequency. However, these circulations are predicted to diminish as deforestation increases. Here we use multi-decadal satellite records and numerical model simulations to show a regime shift in the regional hydroclimate accompanying increasing deforestation in Rondônia, Brazil. Compared to the 1980s, present-day deforested areas in downwind western Rondônia are found to be wetter than upwind eastern deforested areas during the local dry season. The resultant precipitation change in the two regions is approximately ±25% of the deforested area mean. Meso-resolution simulations robustly reproduce this transition when forced with increasing deforestation alone, showing a negligible role of large-scale climate variability. Furthermore, deforestation-induced surface roughness reduction is found to play an essential role in the present-day dry season hydroclimate. Our study illustrates the strong scale-sensitivity of the climatic response to Amazonian deforestation and suggests that deforestation is sufficiently advanced to have caused a shift from a thermally- to a dynamically-driven hydroclimatic regime.
Hogikyan, Allison; Resplandy, Laure; Yang, Wenchang; Fueglistaler, Stephan
Dataset constructed from GFDL-FLOR preindustrial control experiment run by Wenchang Yang (firstname.lastname@example.org) on Princeton University's tiger CPU. Processing by Allison Hogikyan (email@example.com) on Princeton University's tigress data processing node. June 2021.
The carbon isotopic (δ13C) composition of shallow-water carbonates often is interpreted to reflect the δ13C of the global ocean and is used as a proxy for changes in the global carbon cycle. However, local platform processes, in addition to meteoric and marine diagenesis, may decouple carbonate δ13C from that of the global ocean. To shed light on the extent to which changing sediment grain composition may produce δ13C shifts in the stratigraphic record, we present new δ13C measurements of benthic foraminifera, solitary corals, calcifying green algae, ooids, coated grains, and lime mud from the modern Great Bahama Bank (GBB). This survey of a modern carbonate environment reveals δ13C variability comparable to the largest δ13C excursions in the last two billion years of Earth history.
The prevalence of ooids in the stratigraphic record, and their association with shallow-water carbonate environments, make ooids an important paleoenvironmental indicator. Recent advances in the theoretical understanding of ooid morphology, along with empirical studies from Turks and Caicos, Great Salt Lake, and The Bahamas, have demonstrated that the morphology of ooids is indicative of depositional environment and hydraulic conditions. To apply this knowledge from modern environments to the stratigraphic record of Earth history, researchers measure the size and shape of lithified ooids on two-dimensional surfaces (i.e., thin sections or polished slabs), often assuming that random 2D slices intersect the nuclei and that the orientation of the ooids is known. Here we demonstrate that these assumptions rarely are true, resulting in errors of up to 35% on metrics like major axis length. We present a method for making 3D reconstructions by serial grinding and imaging, which enables accurate measurement of the morphology of individual ooids within an oolite, as well as the sorting and porosity of a sample. We also provide three case studies that use the morphology of ooids in oolites to extract environmental information. Each case study demonstrates that 2D measurements can be useful if the environmental signal is large relative to the error from 2D measurements. However, 3D measurements substantially improve the accuracy and precision of environmental interpretations. This study focuses on oolites, but errors from 2D measurements are not unique to oolites; this method can be used to extract accurate grain and porosity measurements from any lithified granular sample.
There has been considerable recent interest in the high-pressure behavior of silicon carbide, a potential major constituent of carbon-rich exoplanets. In this work, the atomic-level structure of SiC was determined through in situ X-ray diffraction under laser-driven ramp compression up to 1.5 TPa; stresses more than seven times greater than previous static and shock data. Here we show that the B1-type structure persists over this stress range and we have constrained its equation of state (EOS). Using this data we have determined the first experimentally based mass-radius curves for a hypothetical pure SiC planet. Interior structure models are constructed for planets consisting of a SiC-rich mantle and iron-rich core. Carbide planets are found to be ~10% less dense than corresponding terrestrial planets.
Vecchi, Gabriel A.; Landsea, Christopher; Zhang, Wei; Villarini, Gabriele; Knutson, Thomas
These are the data and scripts supporting the manuscript: Vecchi, Landsea, Zhang, Villarini and Knutson (2021): Changes in Atlantic Major Hurricane Frequency Since the Late-19th Century. Nature Communications.
This setup mimics ice lying above the drainage system. In the experiment, a fluid-filled blister is generated via liquid injection into the interface between a transparent elastic layer and a porous substrate. After injection of liquid, the fluid permeates from the blister through the porous substrate, the blister volume V(t) relaxes exponentially with time. Our lab experiments show that varying the permeability of the porous substrate k significantly impacts the relaxation timescale in the experiments.
The dataset contains the model file for the Global Adjoint Tomography Model 25 (GLAD-M25). The model file contains parameters defined on the spectral-element mesh and is recommend to be used in SPECFEM3D GLOBE for seismic wave simulation at the global scale.
In our study, we compare the three dimensional (3D) morphologic characteristics of Earth's first reef-building animals (archaeocyath sponges) with those of modern, photosynthetic corals. Within this repository are the 3D image data products for both groups of animals. The archaeocyath images were produced through serial grinding and imaging with the Grinding, Imaging, and Reconstruction Instrument at Princeton University. The images in this repository are the downsampled data products used in our study, and the full resolution (>2TB) image stacks are available upon request from the author. For the coral image data, the computed tomography (CT) images of all samples are included at full resolution. Also included in this repository are the manual and automated outline coordinates of the archaeocyath and coral branches, which can be directly used for morphological study.
The item included here is a collection of wave profiles collected and presented in the accompanying paper: Rucks, M. J., Winey, J. M., Toyoda, T., Gupta, Y. M., & Duffy, T. S. (in review). "Shock compression of fluorapatite to 120 GPa" Submitted to Journal of Geophysical Research: Planets.
This dataset includes individual CIF files with the refined structure of fluorapatite under compression to 61 GPa. The structures have been discussed in detail in the accompanying manuscript "Single-crystal X-ray diffraction of fluorapatite to 61 GPa"
This dataset contains all the model output used to generate the figures and data reported in the article "Climate, soil organic layer, and nitrogen jointly drive forest development after fire in the North American boreal zone". The data was generated during spring 2015 using the a modified version of the Ecosystem Demography model version 2, provided as a supplement accompanying the article. The data was generated using the computational resources supported by the PICSciE OIT High Performance Computing Center and Visualization Laboratory at Princeton University. The dataset contains a pdf Readme file which explains in detail how the data can be used. Users are recommended to go through this file before using the data.
This dataset contains all the data, model and MATLAB codes used to generate the figures and data reported in the article (DOI: 10.1002/2014JD022278). The data was generated during September 2013 and February 2014 using the Ocean-Land-Atmosphere Model also provided with this package. The data was generated using the computational resources supported by the PICSciE OIT High Performance Computing Center and Visualization Laboratory at Princeton University. The dataset contains a pdf Readme file which explains in detail how the data can be used. Users are recommended to go through this file before using the data.