Probing the deep critical zone beneath the Luquillo Experimental Forest, Puerto Rico

Heather L. Buss; Susan L. Brantley; Fred Scatena; Katya Bazilevskaya; Alex E. Blum; Marjorie S. Schulz; Rafael Jimenez; Arthur F. White; G. Rother; D. Cole

doi:10.1002/esp.3409

USGS · 70040792

Probing the deep critical zone beneath the Luquillo Experimental Forest, Puerto Rico

Heather L. Buss·Susan L. Brantley·Fred Scatena·Katya Bazilevskaya·Alex E. Blum·Marjorie S. Schulz·Rafael Jimenez·Arthur F. White·G. Rother·D. Cole

View original source ↗Export metadata ↓

Abstract

Recent work has suggested that weathering processes occurring in the subsurface produce the majority of silicate weathering products discharged to the world's oceans, thereby exerting a primary control on global temperature via the well-known positive feedback between silicate weathering and CO 2 . In addition, chemical and physical weathering processes deep within the critical zone create aquifers and control groundwater chemistry, watershed geometry and regolith formation rates. Despite this, most weathering studies are restricted to the shallow critical zone (e.g. soils, outcrops). Here we investigate the chemical weathering, fracturing and geomorphology of the deep critical zone in the Bisley watershed in the Luquillo Critical Zone Observatory, Puerto Rico, from two boreholes drilled to 37.2 and 27.0 m depth, from which continuous core samples were taken. Corestones exposed aboveground were also sampled. Weathered rinds developed on exposed corestones and along fracture surfaces on subsurface rocks slough off of exposed corestones once rinds attain a thickness up to ~1 cm, preventing the corestones from rounding due to diffusion limitation. Such corestones at the land surface are assumed to be what remains after exhumation of similar, fractured bedrock pieces that were observed in the drilled cores between thick layers of regolith. Some of these subsurface corestones are massive and others are highly fractured, whereas aboveground corestones are generally massive with little to no apparent fracturing. Subsurface corestones are larger and less fractured in the borehole drilled on a road where it crosses a ridge compared with the borehole drilled where the road crosses the stream channel. Both borehole profiles indicate that the weathering zone extends to well below the stream channel in this upland catchment; hence weathering depth is not controlled by the stream level within the catchment and not all of the water in the watershed is discharged to the stream

Explore related subjects

Earth Surface Processes and Landforms Article

Source-reported bounding extent: 18.3093884924° to 18.3233015696° latitude; -65.7478094101° to -65.7264590263° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

Heather L. Buss, Susan L. Brantley, Fred Scatena, Katya Bazilevskaya, Alex E. Blum, Marjorie S. Schulz, Rafael Jimenez, Arthur F. White, G. Rother, D. Cole. 2013-04-30. Probing the deep critical zone beneath the Luquillo Experimental Forest, Puerto Rico. https://doi.org/10.1002/esp.3409

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

USGS2026

Deep learning models for mapping surficial geology in selected physiographic regions of New York

Surficial geologic mapping is required for decisions including infrastructure and water-resource development and management. Deep learning, a type of machine learning that uses training data to self-learn and perform tasks, is explored as a tool to help increase efficiency of the labour- and time-intensive mapping process. Deep learning models were trained, and their potential to aid in mapping surficial geology was explored for two physiographic regions in New York: the high-relief Allegheny Plateau and the low-relief Erie-Ontario Lowlands. Key to the development of deep learning models is the availability of highly detailed surficial geologic maps in each of these regions from which the models can learn. Through exploring different groupings of surficial deposits and their corresponding spatial data, deep learning models were iteratively developed to reproduce published mapping for greater than 79% of the training areas, and with similar accuracy in test areas within the same physiographic region. These models were trained using only two inputs: high-resolution lidar data and previously published surficial geologic maps. This straightforward approach is intended to make the methods and models created easily reproducible using widely accessible datasets. This research describes the strengths and limitations of lidar-derived surficial geologic models and how broadly grouping surficial types, characteristics of different physiographic regions, and testing models in physiographic regions outside of their training areas affect model creation and performance. In addition, the models created could be used as a tool for mapping surficial geology in similar physiographic regions, in addition to establishing a framework for creating similar models elsewhere.

New York↗

USGS2026

Regional models for postfire debris-flow likelihood and rainfall thresholds across the western United States

The U.S. Geological Survey (USGS) uses an empirical model developed with logistic regression (the ‘M1’ model) to rapidly assess debris-flow likelihood and to identify quantitative rainfall thresholds for debris flows after wildfire in the western United States. The M1 model was calibrated to a debris-flow inventory from southern California (United States) and has been applied throughout the western United States. Limited spatial coverage in the calibration dataset has motivated evaluation of M1 model accuracy outside the calibration region (e.g., the Sierra Nevada or the eastern Cascade Range, United States). Previous test cases showed that M1 overpredicts debris-flow likelihood and underpredicts rainfall thresholds for some locations (e.g., Arizona, northern California, Colorado, New Mexico, United States). We sought to improve the regional applicability of a debris-flow likelihood model by expanding the debris-flow inventory used for calibration, testing multiple potential models and generating an updated model framework. The updated inventory includes 3788 observations from 67 burned areas paired with short duration rainfall ratios. The updated model framework consists of a modified model structure and sets of coefficients calibrated separately to the entire updated inventory and to subsets of the inventory that intersect three Environmental Protection Agency (EPA) Level 2 ecoregions (Mediterranean California, Upper Gila Mountains and Western Cordillera). Comparisons of predictions from the updated models with observed rainfall and debris-flow activity show that the updated models outperform the M1 model by ~15%–60% and improve the uniformity of predictive performance across the western United States. The updated models also reduce false positive rates relative to M1 and generate rainfall thresholds that are better aligned with relative differences in regional climatology and debris-flow activity.

Arizona, California, Colorado, Idaho, Montana, Nev↗

USGS2026

Numerical modeling of Late Pleistocene to Holocene earthquake-induced progressive rock slope damage in the paraglacial Serpentine valley, Prince William Sound, Alaska

Landslides in deglaciating fjords pose a potential tsunami threat to nearby communities; however, processes contributing to long-term progressive rock damage and landslide conditioning remain poorly constrained in paraglacial settings. Here, we analyse the role of earthquake-induced rock mass damage over late Pleistocene to Holocene time scales, as a conditioning factor for modern landslides, using distinct element numerical modelling to assess spatial and temporal patterns of fracture propagation influenced by varying glacier thickness. Conceptualised numerical models were parameterised by in situ rock mass, glacial and topographic conditions in Serpentine valley, located in Prince William Sound, Alaska, where several large landslides are actively developing along the western valley wall. Results show that although rigid glacier buttressing reduces co-seismic rock mass damage, it does not suppress it completely, and damage occurs both above and below the glacier surface elevation. We demonstrate that topography, and especially steep slopes with topographic convexity, as well as preexisting damage of joint networks and faults inherited from tectonic and exhumation induced stresses, exert primary control on the location of new co-seismic damage. Simulations representing a simplified deglaciation sequence over the past 25 ka, with a series of 10 evenly spaced earthquakes, generated rock mass damage patterns that qualitatively match in situ landslide structural and kinematic observations at one instability in Serpentine valley. Our conceptual study helps clarify the role of repeated seismicity over glacial timescales as a long-term conditioning process for paraglacial rock slope failure, highlighting spatial and temporal patterns of progressive damage accumulation, with outcomes relevant for modern landslide hazard assessment.

Alaska↗