Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher.
Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?
Some links on this page may take you to non-federal websites. Their policies may differ from this site.
-
Rapid sediment accumulation rates (SAR) in a fan delta situated on the rapidly uplifting footwall of the Taormina normal fault in NE Sicily preserves a rare record of earthquakes and base level change for a tightly coupled source to sink system. We use this sedimentary archive to reconstruct the kinematics and slip history of the fault and further an understanding of how tectonic forcing across various scales are encoded in stratigraphy. A revised luminescence-based age model indicates that ~82 m of the Pagliara fan-delta foreset facies was deposited in ~11 ka at a mean SAR of ~0.74 cm/yr during MIS 7. Syn-depositional terrestrial cosmogenic nuclide (TCN) determined paleoerosion rates of 0.91±0.12 mm/yr and 1.31 ±0.61 mm/yr are similar to published modern erosion rates for the Pagliara basin of 0.97 ±0.11 mm/yr. At the stratigraphic scale, a time series of magnetic susceptibility (c) sampled at 1 m intervals in the foresets displays four ~2,800 yr / 20 m-thick cycles of growing c, bounded by sharp decreases that do not coincide with changes in sediment texture. The c of the low-grade metamorphic bedrock in the source is 20-100 times weaker than the c of rubified soils mantling the hillslopes, which is comparable to the c of the delta sediments. We propose that large, bedrock-cored landslides quasi-periodically deliver weak c sediment to the delta that dilutes a c signal otherwise dominated by the stripping of soil-mantled hillslopes. We propose that centennial-scale recurrence interval earthquakes are most capable at triggering a basin-scale landslide only after channel incision has increased relief of hillslopes to the threshold condition, which requires millennia to achieve. At the landscape scale of delta geometry and location, the Pagliara delta accumulated in a hanging wall basin that has since been inverted. We reconstruct the history of base level fall for the delta from an inversion of fluvial topography and apportion that record to its rock uplift, delta deposition, and eustatic components. We show that footwall uplift has been unsteady over the past 600 ka ranging from -1 to 3 mm/yr. The integration of our stratigraphic- and landscape scale observations furthers our understanding of the natural hazards related to normal fault earthquakes and their impact on sediment dynamics in this steep, active tectonic setting.more » « lessFree, publicly-accessible full text available September 30, 2025
-
We explore how rock properties and channel morphology vary with rock type in Last Chance Canyon, Guadalupe Mountains, New Mexico, USA. The rocks here are composed of horizontally to near-horizontally interbedded carbonate and sandstone. This study focuses on first- and second-order channel sections, where the streams have a lower channel steepness index (ksn) upstream and transition to higher ksn values downstream. We hypothesize that differences in bed thickness and rock strength influence ksn values, both locally by influencing bulk bedrock strength and also nonlocally through the production of coarse sediment. We collected discontinuity intensity data (the length of bedding planes and fractures per unit area), Schmidt hammer rebound measurements, and measured the largest boulder at every 12.2 m elevation contour to test this hypothesis. Bedrock and boulder mineralogy were determined using a lab-based carbonate dissolution method. High-resolution orthomosaics and digital surface models (DSMs) were generated from drone and ground-based photogrammetry. The orthomosaics were used to map channel sections with exposed bedrock. The United States Geological Survey (USGS) 10 m digital elevation models (DEMs) were used to measure channel slope and hillslope relief. We find that discontinuity intensity is negatively correlated with Schmidt hammer rebound values in sandstone bedrock. Channel steepness tends to be higher where reaches are primarily incising through more thickly bedded carbonate bedrock and lower where more thinly bedded sandstone is exposed. Bedrock properties also influence channel morphology indirectly, through coarse sediment input from adjacent hillslopes. Thickly bedded rock layers on hillslopes erode to contribute larger colluvial sediment to adjacent channels, and these reaches have higher ksn values. Larger and more competent carbonate sediment armors both the carbonate and the more erodible sandstone and reduces steepness contrasts across rock types. We interpret that in the relatively steep, high-level ksn downstream channel sections, the slope is primarily controlled by the coarse alluvial cover. We further posit that the upstream low-level ksn reaches have a base level that is fixed by the steep downstream reaches, resulting in a stable configuration, where channel slopes have adjusted to lithologic differences and/or sediment armor.more » « less
-
Two major barriers hinder the holistic understanding of subsurface critical zone (CZ) evolution and its impacts: (a) an inability to measure, define, and share information and (b) a societal structure that inhibits inclusivity and creativity. In contrast to the aboveground portion of the CZ, which is visible and measurable, the bottom boundary is difficult to access and quantify. In the context of these barriers, we aim to expand the spatial reach of the CZ by highlighting existing and effective tools for research as well as the “human reach” of CZ science by expanding who performs such science and who it benefits. We do so by exploring the diversity of vocabularies and techniques used in relevant disciplines, defining terminology, and prioritizing research questions that can be addressed. Specifically, we explore geochemical, geomorphological, geophysical, and ecological measurements and modeling tools to estimate CZ base and thickness. We also outline the importance of and approaches to developing a diverse CZ workforce that looks like and harnesses the creativity of the society it serves, addressing historical legacies of exclusion. Looking forward, we suggest that to grow CZ science, we must broaden the physical spaces studied and their relationships with inhabitants, measure the “deep” CZ and make data accessible, and address the bottlenecks of scaling and data‐model integration. What is needed—and what we have tried to outline—are common and fundamental structures that can be applied anywhere and used by the diversity of researchers involved in investigating and recording CZ processes from a myriad of perspectives.more » « less
-
The Earth sciences lack representation of Latine scholars in academia and industry. Latine scholars encounter systematic barriers within academic environments globally, including discrimination and gender biases, lack of guidance due to cultural differences, funding limitations for non-citizens, access burden of visa applications, and language-barriers for non-native English speakers (Carrera et al., 2023; Valenzuela-Toro and Viglino, 2021).more » « less
-
Key Points Modeled ecosystem response to climate follows the “geo‐ecological law of distribution,” highlights the importance of ecohdyrologic refugia Woody Plant Encroachment is predicted as a three‐phase phenomenon: early establishment, rapid expansion, and woody plant equilibrium Regime shifts from grassland to shrubland are marked by vegetation cover thresholdsmore » « less
-
Abstract. Numerical simulation of the form and characteristics of Earth's surface provides insight into its evolution. Landlab is an open-source Python package that contains modularized elements of numerical models for Earth's surface, thus reducing time required for researchers to create new or reimplement existing models. Landlab contains a gridding engine which represents the model domain as a dual graph of structured quadrilaterals (e.g., raster) or irregular Voronoi polygon–Delaunay triangle mesh (e.g., regular hexagons, radially symmetric meshes, and fully irregular meshes). Landlab also contains components – modular implementations of single physical processes – and a suite of utilities that support numerical methods, input/output, and visualization. This contribution describes package development since version 1.0 and backward-compatibility-breaking changes that necessitate the new major release, version 2.0. Substantial changes include refactoring the grid, improving the component standard interface, dropping Python 2 support, and creating 31 new components – for a total of 58 components in the Landlab package. We describe reasons why many changes were made in order to provide insight for designers of future packages. We conclude by discussing lessons about the dynamics of scientific software development gained from the experience of using, developing, maintaining, and teaching with Landlab.more » « less