SEARCH · USGS Science

Results for “North American Journal of Fisheries Management”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,662 records · Page 40Linked to original sources

USGS2025

Case study of deep learning image segmentation for the purposes of rapid 2D petrographic analysis in volcanic rocks

Automation using deep learning methods is a useful alternative to manual methods of petrographic segmentation, but often requires user familiarity with coding and/or algorithms. We examine the Dragonfly TM program's deep learning tools for application by users with a variety of skill levels as a method for petrographic image segmentation. An image processing methodology, bimodal image stacking, was created for low-input-data, high-efficacy training of models which can then be applied to varied samples. Using backscatter electron images we show that the resulting model segmentations agree with manual segmentation total and modal crystallinity values within 5%, and calculated plagioclase crystal size distribution (CSD) values within 2σ, despite limitations in discriminating mafic phases. Model creation and training takes <24 hours, 1–3 hours of which are supervised, and the resultant model can then be applied to new uncharacterized samples in <15 minutes per image. This allows for non-experts to create and utilize deep learning models to segment images of variable brightness and texture, at low user-time cost and resulting in size and shape data which are within uncertainty of manual segmentation. While some limitations are noted (for example, sieve-textured phases may need manual correction, and different minerals with similar BSE intensity may not be resolved as separate phases), this methodology can be utilized for general application of models to wide ranges of volcanic crystalline and bubble textures, and to create a library of models for rapid petrological analysis during volcanic eruptions.

Alaska↗

USGS1995

Methods for estimating the magnitude and frequency of peak discharges of rural, unregulated streams in Virginia

Methods are presented for estimating the peak discharges of rural, unregulated streams in Virginia. A Pearson Type III distribution is fitted to the logarithms of the unregulated annual peak-discharge records from 363 stream-gaging stations in Virginia to estimate the peak discharge at these stations for recurrence intervals of 2 to 500 years. Peak-discharge characteristics for 284 unregulated stations are divided into eight regions based on physiographic province, and regressed on basin characteristics, including drainage area, main channel length, main channel slope, mean basin elevation, percentage of forest cover, mean annual precipitation, and maximum rainfall intensity. Regression equations for each region are computed by use of the generalized least-squares method, which accounts for spatial and temporal correlation between nearby gaging stations. This regression technique weights the significance of each station to the regional equation based on the length of records collected at each cation, the correlation between annual peak discharges among the stations, and the standard deviation of the annual peak discharge for each station. Drainage area proved to be the only significant explanatory variable in four regions, while other regions have as many as three significant variables. Standard errors of the regression equations range from 30 to 80 percent. Alternate equations using drainage area only are provided for the five regions with more than one significant explanatory variable. Methods and sample computations are provided to estimate peak discharges at gaged and engaged sites in Virginia for recurrence intervals of 2, 5, 10, 25, 50, 100, 200, and 500 years, and to adjust the regression estimates for sites on gaged streams where nearby gaging-station records are available.

Virginia↗

USGS2025

Footprints of past mining in Alaska (USA) derived from high-resolution satellite imagery

Mapping the land area used for mining in the past is essential for guiding the remediation of affected landscapes and assessing the resource potential of related waste products. Despite significant recent progress delineating footprints of active and inactive mining globally, the known inventory of such mine lands remains incomplete. Here, I describe a new map dataset of footprints of land surface disturbance and waste at sites of past mining in Alaska (USA) based on visual interpretation of satellite imagery. This dataset maps 6–14 times the area of previous regional and global mine footprint maps in Alaska and is the first in the region to explicitly delineate mine waste landforms (e.g., tailings piles). The data are publicly available from the U.S. Geological Survey under a “no rights reserved” Creative Commons (CC0) license agreement.

Alaska↗

USGS2025

Enhanced hydrologic monitoring and characterization of groundwater drainage features

Groundwater drains to the land surface, generating the baseflow of streams, lakes, and wetlands. The hydrologic resilience of baseflow during prolonged dry periods and after disturbance can be assessed with evolving remote sensing analysis paired with localized monitoring of groundwater drainage features and creative model calibration strategies.

Nature Water↗

USGS2024

The dominance and growth of shallow groundwater resources in continuous permafrost environments

Water is a limited resource in Arctic watersheds with continuous permafrost because freezing conditions in winter and the impermeability of permafrost limit storage and connectivity between surface water and deep groundwater. However, groundwater can still be an important source of surface water in such settings, feeding springs and large aufeis fields that are abundant in cold regions and generating runoff when precipitation is rare. Whether groundwater is sourced from suprapermafrost taliks or deeper regional aquifers will impact water availability as the Arctic continues to warm and thaw. Previous research is ambiguous about the role of deep groundwater, leading to uncertainty regarding Arctic water availability and changing water resources. We analyzed chemistry and residence times of spring, stream, and river waters in the continuous permafrost zone of Alaska, spanning the mountains to the coastal plain. Water chemistry and age tracers show that surface waters are predominately sourced from recent precipitation and have short (<50 y) subsurface residence times. Remote sensing indicates trends in the areal extent of aufeis over the last 37 y, and correlations between aufeis extent and previous year summer temperature. Together, these data indicate that surface waters in continuous permafrost regions may be impacted by short flow paths and shallow suprapermafrost aquifers that are highly sensitive to climatic and hydrologic change over annual timescales. Despite the lack of connection to regional aquifers, continued warming and permafrost thaw may promote deepening of the shallow subsurface aquifers and creation of shallow taliks, providing some resilience to Arctic freshwater ecosystems.

Proceedings of the National Academy of Sciences↗

USGS2024

On the provenance of field reports of the 1886 Charleston, South Carolina, earthquake: A seismo-historical whodunnit

Much of what is known about the effects of the 1886 Charleston, South Carolina, earthquake throughout the epicentral region can be attributed to meticulous field investigations by an individual with training in geology and engineering, Earle Sloan ( Clendenin, 1926 ). In a recent study, Bilham and Hough (2024) undertook a detailed analysis of the effects of the earthquake on railroads in the Charleston region, drawing heavily from Sloan’s reports. This exercise identified several inconsistencies in Sloan’s field reports, including understandable measurement imprecision, inferred data entry mistakes, and transcription errors. The study also begged the question, where was Sloan at the time of the mainshock and over the following week? And to what extent did he draw from secondhand information in compiling his reports? On this question Sloan’s reports were sometimes enigmatic, lending themselves to misinterpretation in contemporaneous as well as modern interpretations. Beyond the details that were germane for, and briefly summarized by, the studies of Bilham and Hough (2023 , 2024) , in this report we don our historical seismologist caps to chronicle Sloan’s activities following the earthquake. We summarize our inferences here for the benefit of future scholars who might attempt to retrace either Sloan’s footsteps or our own. This study also serves to highlight Sloan’s singular contributions to earthquake science, which were never published separately.

South Carolina↗

USGS2026

A novel drive-point multilevel system to investigate PFAS and other contaminants of global concern in the hyporheic zone of a wastewater effluent dominated stream

Contaminants found in treated wastewater discharged to streams, including pharmaceuticals and per- and polyfluoroalkyl substances (PFAS), are of global concern due to their deleterious effects on aquatic ecosystems and potential impacts to human health. Hyporheic zones have strong potential for contaminant attenuation. Assessing this potential requires collection of physical and biogeochemical data within the hyporheic zone. This study tested the applicability of a novel drive-point multilevel system (DP-MLS) for quantifying head profiles and characterizing contaminant concentrations in the hyporheic zone of a temperate region effluent dominated stream (EDS). DP-MLS, each with 4 ports, were installed in the stream bed at two sites, DS-1 and DS-2, 0.2 and 4.7 km downstream of the effluent outfall, respectively. Head profiles were measured and groundwater collected for analysis of pharmaceuticals and PFAS temporally over two years. The DP-MLS withstood rapid changes in stage, ice formation, and floating debris. Vertical hydraulic gradients (VHG) were generally upward but varied in magnitude indicating heterogeneity in hydraulic conductivity and variability in flow conditions. Upward VHG were also about 2X larger at DS-1 than at DS-2. Contaminant concentration profiles consistently showed penetration of pharmaceuticals and PFAS to 1 m below the bed at DS-2 while there was less penetration, lower groundwater concentrations, and more temporal variability in concentrations at DS-1. Integration of the physical and chemical data suggests weaker upwelling conditions at DS-2 are more easily reversed during periods of high stream stage, which could facilitate migration of wastewater contaminants into the bed. However, further studies incorporating other transport processes and reach scale dynamics are required to fully characterize these exchanges. Overall, this study demonstrates the efficacy of these novel DP-MLSs for characterization of the hyporheic zone and provides new insights into the occurrence, composition, and persistence of wastewater derived contaminants in the hyporheic zone of a well-studied EDS.

Iowa↗

USGS2026

The eruptive behavior of distributed volcanism forming low shield edifices—A case study of Sentinel-Arlington volcanic field, U.S.A.

Distributed volcanic fields are present in various tectonic settings worldwide, and their characteristics reflect differing influences from magmatic and tectonic processes. In the southwestern United States alone, there are 37 Quaternary distributed volcanic fields. After the primary period of extensional tectonics in the southern Basin and Range 15–5 million years ago, the Sentinel-Arlington volcanic field developed in southwestern Arizona between 4 and 1 million years ago. The Sentinel-Arlington volcanic field consists primarily of low relief shield volcanoes, a type of distributed volcanism with poorly understood temporal evolution. The Sentinel-Arlington volcanic field is less than 200 kilometers (km) from the Colorado Plateau, Gulf of California, and southern San Andreas Fault system. This work identifies and examines controls on the emplacement of the Sentinel-Arlington volcanic field by documenting shallow and surficial structures as well as eruption characteristics and style through time. The Sentinel-Arlington volcanic field consists of 21 volcanoes with a total of 33 vents over an area of about 770 square kilometers (km 2 ). The prominence of low relief shield volcanoes may be explained by ascent of basaltic magmas through thin Basin and Range crust, without much crustal contamination, and low viscosities common to mafic compositions. Typical eruption characteristics involve the construction of low relief shield volcanoes followed by Strombolian fissure eruptions at the summits or near-summit medial areas that produce scoria lapilli, which may weld to form agglomerate. The total lifetime erupted volume of about 4.3 cubic kilometers (km 3 ) represents an average eruptive flux of approximately 2x10 -3 km 3 per thousand years (k.y.). This erupted volume is low relative to Neogene basaltic intraplate distributed volcanic fields worldwide, which typically range from 10 -3 to 1 km 3 k.y. -1 . Sentinel-Arlington volcanic field eruptions were likely triggered by intermittent rejuvenation of transient magmatic zones that exist in thinned crust. Instantaneous flux from point sources feeding the lava flows is estimated to be on the order of 10 -1 to 10 cubic meters per second.

Arizona↗

USGS1996

Techniques to estimate generalized skew coefficients of annual peak streamflow for natural basins in Texas

This report presents two techniques to estimate generalized skew coefficients used for log-Pearson Type III peak-streamflow frequency analysis of natural basins in Texas. A natural basin has less than 10 percent impervious cover, and less than 10 percent of its drainage area is controlled by reservoirs. The estimation of generalized skew coefficients is based on annual peak and historical peak streamflow for all U.S. Geological Survey streamflow-gaging stations having at least 20 years of annual peak-streamflow record from natural basins in Texas. Station skew coefficients calculated for each of 255 Texas stations were used to estimate generalized skew coefficients for Texas. One technique to estimate generalized skew coefficients involved the use of regression equations developed for each of eight regions in Texas, and the other involved development of a statewide map of generalized skew coefficients. The weighted mean of the weighted mean standard errors of the regression equations for the eight regions is 0.36 log 10 skew units, and the weighted mean standard error of the map is 0.35 log 10 skew units. The technique based on the map is preferred for estimating generalized skew coefficients because of its smooth transition from one region of the State to another.

Texas↗

USGS2024

A global assessment of SAOCOM-1 L-band stripmap data for InSAR characterization of volcanic, tectonic, cryospheric, and anthropogenic deformation

SAOCOM-1 is an L-band (23.5 cm) synthetic aperture radar (SAR) constellation made up of two satellites launched in 2018 and 2020 by Comisión Nacional de Actividades Espaciales (CONAE, Argentina). In this contribution, we present a global summary of interferometric SAR (InSAR) observations of ground deformation with SAOCOM-1 stripmap data for tracking volcanic, tectonic, glacier, and anthropogenic deformation. These examples include: 1) episodes of unrest at volcanoes in the Aleutian Islands, Southern Andes, and Italy, with line-of-sight (LOS) deformation from 4 cm/yr in InSAR time series to ~70 cm in interferograms; 2) dike intrusions in Hawai’i; 3) earthquakes in the Andean fold and thrust belt and the East Anatolian fault; 4) ice flow of the Southern Patagonia icefield; and 5) subsidence due to lithium brine extraction in the Salar de Atacama basin (northern Chile). Comparisons between SAOCOM-1, ALOS-2 SM3, Sentinel-1, and TerraSAR-X/ TanDEM-X/PAZ (TSX/TDX/PAZ) mean velocities from InSAR time series show a 1:1 ± 3% correlation in the LOS velocity, which highlights the high accuracy of SAOCOM-1 data. The minimum deformation that we measured in individual interferograms is 4 ± 0.6 cm. One limitation of SAOCOM-1 is the lack of a global acquisition program, which reduces its global and broader applications. Considering the repeat periods, background observation program, and lack of a controlled orbital tube, the best suited targets for SAOCOM-1 InSAR are two. First, volcanoes that deform with secular rates located in vegetated regions in mid- and high-latitudes, and/or that undergo transient episodes of fast deformation in which C-band coherence is lost quickly. Second, glaciers where coherence can be sustained during the repeat period of eight days.

IEEE Transactions on Geoscience and Remote Sensing↗

USGS2025

Predictive regressive models of recent marsh sediment thickness improve the quantification of coastal marsh sediment budgets

Coastal marsh wetlands experience variations in vertical gains and losses through time, which have allowed them to infill relict topography and record variations in drivers. The stratigraphic unit associated with the development of the marsh also reflects the long-term importance of key ecosystem services supplied by the marsh environment, including carbon storage and storm mitigation. Mapping these coastal wetland sediments and the marsh unit thickness is challenging as traditional coastal geophysical tools are not easily deployable (acoustic methods) or are unreliable in saline-soil environments (e.g., ground-penetrating radar), leaving core-based methods the most viable mapping method. In the present study, we utilized prior information on the geologic architecture of the region to select spatial and physical metrics that likely persisted throughout evolution of the marsh during the late Holocene. We then assessed the individual and collective power of these metrics to predict marsh thickness observed from cores. Employing regressive predictive models powered by these data, we improve the quantification of marsh thickness for a coastal fringing marsh within the Grand Bay estuary in Mississippi and Alabama (USA). The information gained from this approach yields improved estimates of the carbon stocks in this environment. Additionally, the stored sediment masses reflect the past, and potential future, persistence of the Grand Bay marsh under historical and present marsh-estuarine sediment exchange fluxes. Such improvements to both the sediment budget of recent marsh stratigraphic units and the spatial extent provide new resources for comparison with large-scale landscape models, the latter of which may be used, when validated, to predict future change and ecosystem transformations.

Alabama. Mississippi↗

USGS2024

Use of a numerical groundwater-flow model and projected climate scenarios to simulate the effects of future climate conditions on base flow for reach 1 of the Washita River alluvial aquifer and Foss Reservoir storage, western Oklahoma

To better understand the relation between climate variability and future groundwater resources in reach 1 of the Washita River alluvial aquifer and Foss Reservoir in western Oklahoma, the U.S. Geological Survey, in cooperation with the Bureau of Reclamation, used a previously published numerical groundwater-flow model and climate-model data to investigate changes in base flow and reservoir storage by evaluating three scenarios. The three projected climate scenarios were (1) a central-tendency scenario, (2) a warmer/drier scenario, and (3) a less-warm/wetter scenario. To estimate future base flow and groundwater availability in western Oklahoma, specifically in reach 1 of the Washita River alluvial aquifer, downscaled climate-model data from 231 Coupled Model Intercomparison Project phase 5 (CMIP5) projections coupled with a previously published numerical groundwater-flow model were used to compare the effects of different climate scenarios on the aquifer. Changes in base flow and groundwater-level elevations during a 30-year baseline scenario (1985–2014) and the three 30-year projected climate scenarios (2050–79) under central-tendency, warmer/drier, and less-warm/wetter climatic conditions were assessed by using the calibrated model. In the simulations, the amount of base flow and reservoir storage declined in the central-tendency and warmer/drier scenarios compared to the amount of base flow and reservoir storage under historical climatic conditions (baseline scenario). Mean annual change in reservoir storage decreased from the baseline scenario the most in the warmer/drier scenario, followed by the central-tendency scenario, but increased in the less-warm/wetter scenario compared to the baseline scenario. At the end of the simulation period (2079), the largest magnitude differences in groundwater-level elevations in all three projected climate scenarios relative to the baseline scenario occurred upstream from Foss Reservoir. Results from incorporating downscaled climate projections into localized numerical groundwater-flow models can highlight potential future changes in and implications for groundwater resources and availability.

Oklahoma↗

USGS2026

Subaerially exposed Iceberg Lake sediments: An exceptional record of historical subaqueous earthquake disturbance at the eastern edge of the Alaskan-Aleutian subduction zone

Paleoseismic records are limited in the Yakutat Terrane (eastern edge of the Alaskan-Aleutian subduction zone) due to the extensive ice cover that hinders traditional methods such as trenching of the faults, but lacustrine sediments offer an alternative archive. We investigated lakebed sediments exposed after recent outburst floods (1999 CE) at Iceberg Lake, a glacier-dammed basin whose stratigraphy was revealed after the lake drained and partially eroded. We logged outcrops across the basin and sampled their sediments. Between annually laminated background deposits, we identified sediment gravity flow beds and in-situ soft-sediment deformation structures (convolutions, sand blows, and fractures) interpreted to be earthquake-induced. Our age model links some of the uppermost depositional and deformational events to the 1958, 1964 and 1979 CE earthquakes. These results demonstrate that Iceberg Lake was a sensitive recorder of seismic shaking and its sediments hold strong potential for producing a paleoseismic record for the northern Yakutat Terrane.

Alaska↗

USGS2024

Prediction of regional broadband strong ground motions using a teleseismic source model of the 18 April 2014 Mw 7.3 Papanoa, Mexico, earthquake

To estimate predicted ground motion from a teleseismic slip model, we use a low‐ and high‐frequency hybrid method to simulate the regional, strong ground motions observed following the 18 April 2014 moment magnitude ( ⁠M w ⁠ ) 7.3 Papanoa, Mexico, earthquake. To generate the regional ground motion at low frequencies (<1 Hz), a teleseismically derived, finite‐fault, kinematic model is used to define the earthquake source, taking into account slip‐model variations identified with a parameter sampling approach that considers possible errors in the fault geometry, the hypocenter depth, and the rupture velocity. A 3D crustal model is used to calculate the low‐frequency ground motions using a finite‐element calculation that includes topography and considers variations in the source model to estimate the uncertainty in the calculations. High frequencies (>1 Hz) are added using a 1D full‐wave propagation code that estimates uncertainties by considering multiple random distributions of slip with different spatial correlation lengths. The synthetic, broadband (0.05–10.0 Hz) ground motions are obtained by combining the low‐ and high‐frequency portions match filtered at 1 Hz. These synthetic ground motions are compared with the regional observations using velocity records, peak ground acceleration, and medians of the orientation‐independent response spectra of the horizontal components (RotD50) calculated at periods of 0.2, 0.3, 0.5, 1.0, 2.0, 3.0, 5.0, 7.5, and 10.0 s. The results indicate that ground motions estimated at these periods using our hybrid approach based primarily on a teleseismically derived source model are comparable to the values observed for the 2014 Papanoa earthquake at regional distances. The approach could be used to estimate strong‐motion spectral levels expected for regions with limited local and regional recordings and could also fill in magnitude or distance gaps in ground‐motion prediction relations utilized in the assessment of seismic hazard.

Papanoa↗

USGS2025

The systematics of stable hydrogen (δ2H) and oxygen (δ18O) isotopes and tritium (3H) in the hydrothermal system of the Yellowstone Plateau volcanic field, USA

To improve our understanding of hydrothermal activity on the Yellowstone Plateau volcanic field, we collected and analyzed a large data set of δ 2 H, δ 18 O, and the 3 H concentrations of circum-neutral and alkaline waters. We find that (a) hot springs are fed by recharge throughout the volcanic plateau, likely focused through fractured, permeable tuff units. Previous work had stressed the need for light δ 2 H water recharge restricted to the northern part of the plateau or recharge during past cold periods. However, new data from the Y-7 drill hole suggests that recharge is not restricted to a certain area or a cold period. (b) δ 18 O values of thermal waters in the geyser basins are shifted from the global meteoric water line by temperature-dependent water-rock reactions with higher subsurface temperatures resulting in a greater shift. (c) Large temporal variations in the isotopic composition of meteoric water recharge and small temporal variability in the isotopic composition of hot spring discharge implies that the volume of groundwater in, and around the Yellowstone caldera is substantially larger than the volume of annual water recharge. (d) Hot springs discharged through different rhyolitic units correlate with identifiable differences in δ 2 H and δ 18 O compositions, 3 H concentrations, and water chemistry that imply equilibration at different temperatures and travel along different flow paths. (e) Based on measured 3 H concentrations, we calculate that hot spring waters in the central part of the geyser basins mostly contain <2% post-1950 meteoric water, whereas waters discharged at the basin margins contain larger fractions of post-1950s meteoric water.

Wyoming↗

USGS2010

Divisions of geologic time—Major chronostratigraphic and geochronologic units

Effective communication in the geosciences requires consistent uses of stratigraphic nomenclature, especially divisions of geologic time. A geologic time scale is composed of standard stratigraphic divisions based on rock sequences and is calibrated in years. Over the years, the development of new dating methods and the refinement of previous methods have stimulated revisions to geologic time scales. Advances in stratigraphy and geochronology require that any time scale be periodically updated. Therefore, Divisions of Geologic Time, which shows the major chronostratigraphic (position) and geochronologic (time) units, is intended to be a dynamic resource that will be modified to include accepted changes of unit names and boundary age estimates. This fact sheet is a modification of USGS Fact Sheet 2007-3015 by the U.S. Geological Survey Geologic Names Committee.

Fact Sheet↗

USGS2026

Constraining source and path effects of large magnitude earthquakes using ground motion simulations

The purpose of this study is to use ground‐motion simulations to investigate ways in which source and path effects for large‐magnitude earthquakes can be represented in nonergodic ground‐motion models (GMMs). To achieve this, we designed a ground‐motion study in the San Francisco Bay Area that includes earthquakes with a broad range of magnitudes distributed uniformly on a fault plane, and sites covering a large range of rupture distances and azimuths. After running a large suite of kinematic simulations (magnitude 4–7), we then develop a nonergodic GMM with the simulated data. We find that trends in the within‐site residuals are affected significantly by the earthquake radiation pattern, rupture directivity, and slip patterns. Next, we modify an existing rupture directivity model to fit and remove the observed radiation pattern and rupture directivity from the residuals. We also minimize the contributions of slip patterns by averaging the within‐site residuals among multiple source realizations. Finally, after removing the source effects from the within‐site residuals, we compare the path effects computed with different magnitude groups using two approaches. The first approach only considers the small events that have the same shortest path to a site as the large events, whereas the second approach considers all small events on the fault plane. The results indicate that it is difficult to satisfactorily approximate the path effects of large events with those of small events using either approach, at least in the case of simulations.

California↗

USGS2025

The impact of burial diagenesis on soil-formed minerals in paleosols using stable isotopes of phyllosilicates and carbonate clumped isotopes

To understand the effects of burial diagenesis on the stable isotope geochemistry of soil-formed clay and carbonate minerals in paleosols, samples were collected from seven cores, spanning middle- to upper-Pennsylvanian strata of the Illinois Basin, with varied maximum burial depths of 1–3 km. Mixed-layer illite-smectite and kaolinite mixtures give δ 2 H and δ 18 O values of −83 ‰ to −36 ‰ and 11.9 ‰ to 21.1 ‰ (VSMOW), respectively. After carbonates were screened petrographically for diagenetic textures using transmitted light and cathodoluminescence, measured clumped isotope Δ 47 values range from 0.504 to 0.563 ‰ (I-CDES). Resulting mineral formation temperatures for phyllosilicate mineral mixtures are 28 to 66 °C (mean = 47 °C), whereas T(Δ 47 ) estimates for calcites are 36 to 61 °C (mean = 45 °C). Calculated δ 18 O water values from which phyllosilicate minerals and calcites precipitated under isotopic equilibrium ranges from −7.1 to −1.2 ‰ and − 1.4 to +4.9 ‰, respectively. Closed and open-system phyllosilicate-fluid exchange modeling indicates that phyllosilicate alteration occurred in the presence of a low temperature brine or meteoric water and is interpreted to occur in a layer-by-layer illitization transformation. Due to the lack of diagenetic textures and positively correlated T(Δ 47 ) and δ 18 O water , calcites are interpreted to have undergone solid-state bond reordering. Despite low to moderate temperatures (<125 °C) and varying depths of shallow burial (1–3 km), solid-state transformation of phyllosilicates and calcites indicates paleosols had prolonged exposure to burial conditions which has implications for the use of paleosol minerals for paleoenvironmental reconstructions.

Illinois, Indiana, Kentucky↗