Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Natural Hazards and Earth Systems Sciences (NHESS)”

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

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

1,668 records · Page 74Linked to original sources
USGS2000

Monthly variability and possible sources of nitrate in ground water beneath mixed agricultural land use, Suwannee and Lafayette Counties, Florida

In an area of mixed agricultural land use in Suwannee and Lafayette Counties of northern Florida, water samples were collected monthly from 14 wells tapping the Upper Floridan aquifer during July 1998 through June 1999 to assess hydrologic and land-use factors affecting the variability in nitrate concentrations in ground water. Unusually high amounts of rainfall in September and October 1998 (43.5 centimeters total for both months) resulted in an increase in water levels in all wells in October 1998. This was followed by unusually low amounts of rainfall during November 1998 through May 1999, when rainfall was 40.7 centimeters below 30-year mean monthly values. The presence of karst features (sinkholes, springs, solution conduits) and the highly permeable sands that overlie the Upper Floridan aquifer provide for rapid movement of water containing elevated nitrate concentrations to the aquifer. Nitrate was the dominant form of nitrogen in ground water collected at all sites and nitrate concentrations ranged from less than 0.02 to 22 milligrams per liter (mg/L), as nitrogen. Water samples from most wells showed substantial monthly or seasonal fluctuations in nitrate concentrations. Generally, water samples from wells with nitrate concentrations higher than 10 mg/L showed the greatest amount of monthly fluctuation. For example, water samples from six of eight wells had monthly nitrate concentrations that varied by at least 5 mg/L during the study period. Water from most wells with lower nitrate concentrations (less than 6 mg/L) also showed large monthly fluctuations. For instance, nitrate concentrations in water from four sites showed monthly variations of more than 50 percent. Large fluctuations in nitrate concentrations likely result from seasonal agricultural practices (fertilizer application and animal waste spreading) at a particular site. For example, an increase in nitrate concentrations observed in water samples from seven sites in February or March 1999 most likely results from application of synthetic fertilizers during the late winter months. Lower nitrate concentrations were detected in water samples from five of eight wells sampled during high-flow conditions for the Suwannee River in March 1998 compared to low-flow conditions in November 1998. Evidence for reduction of nitrate due to denitrification reactions was observed at one site (AC-1), as indicated by elevated concentrations of nitrogen gas and a corresponding increase in nitrogen isotope (d15N-NO3) values with a decrease in nitrate concentrations. Denitrification is unlikely at other sites based on the presence of dissolved oxygen concentrations greater than 2 mg/L in ground water and no observed trend between nitrate concentrations and values d15N-NO3 values. Nitrate was the dominant nitrogen species in most monthly rainfall samples; however, ammonium concentrations were similar or greater than nitrate during November and December 1998. During February through May 1999, both nitrate and ammonium concentrations were substantially higher in monthly rainfall samples collected at the study area compared to mean monthly concentrations at the Bradford Forest site located east of the study area, which is part of the National Atmospheric Deposition Program/National Trends Network. Also, higher nitrogen deposition rates in the study area compared to those at Bradford Forest could indicate that substantial amounts of ammonia are volatilized from fertilizers and animal wastes, released to the atmosphere, and incorporated as nitrate and ammonium in rainfall deposited in the middle Suwannee River Basin. Ground-water samples from most sites had d15N-NO3 values that indicated a mixture of inorganic and organic sources of nitrogen, which corresponded to multiple land uses where both synthetic fertilizers and manure are used on fields near these sites. Distinct d15N-NO3 signatures, however, were observed at some sites. For example, water samples from areas of row-crop farming as the dominant land use had d15N-NO3 values less than 4 per mil, indicating an inorganic nitrogen source such as synthetic fertilizer. In contrast, d15N-NO3 values greater than 9 per mil were found in water samples from three sites where manure from dairy and or poultry operations is spread on fields throughout the year.

Florida↗
USGS2025

Origin and evolution of mafic volcanism associated with 3 m.y. of andesite production at the Goat Rocks volcanic cluster, southern Washington Cascade Range

More than 3 m.y. of mafic volcanism near the Goat Rocks volcanic cluster in the southern Washington Cascade Range, USA, lends insight into the evolution of basalts and the subarc mantle at a long-lived, major arc volcanic locus. We contribute field observations, 40 Ar/ 39 Ar dates, paleomagnetic directions, and bulk rock and mineral compositions to characterize nine mafic units that erupted in association with the Goat Rocks volcanic cluster. The time frame of mafic volcanism, ca. 3.6 Ma to 60 ka, encompasses the lifespan of the central volcanic cluster (3.1 Ma to 115 ka), with a lull from ca. 2.7 Ma to 1.4 Ma. A climactic period of voluminous mafic activity and far-traveled lava flows, including construction of the Hogback Mountain shield volcano, coincided with voluminous andesite eruptions from the central volcanic cluster. The basaltic rocks in the Goat Rocks area are calc-alkaline to barely tholeiitic and have high field strength element depletion relative to large-ion lithophile elements characteristic of calc-alkaline basalts (CAB) of the Cascade volcanic arc. Unlike at neighboring andesitic volcanic centers (Mounts Adams, St. Helens, and Rainier), no other mafic end members such as high-aluminum olivine tholeiite (HAOT) or intraplate-type basalt (IPB) are present at or near the Goat Rocks volcanic cluster, although some of the calc-alkaline basalts in this study have IPB-like affinities. The Goat Rocks mafic units exhibit two main temporal trends in composition: (1) the most primitive basalts erupted earlier, compared to less primitive and more evolved compositions later, and (2) high field strength element concentrations are higher in the younger basalt units relative to the oldest two. In contrast to these temporal trends, the mafic units define two compositional groups that recur through time, a low-Sr and a high-Sr group, each with distinct trace element and Sr and Nd isotope ratios. Although radiogenic isotope ratios are generally aligned with High Cascades CAB and HAOT, some extend toward IPB of Mount Adams and Simcoe Mountains volcanic field. Olivine-dominated crystal fractionation at shallow pressure from a small range of parent magma compositions accounts for much of the variation among the basalts and basaltic andesites. A high-pressure fractionation model is plausible for only one of the youngest basalt units (basalt of Walupt Lake volcano). Mafic recharge and crustal assimilation accounts for the incompatible-element enriched composition of basaltic andesites erupted during construction of the largest andesitic centers, further supporting sustained basalt mass flux and thermal energy driving andesite genesis. We model the most primitive members of the Goat Rocks mafic units as partial melts of successively less depleted mantle in time. Variable degrees of fluxing with fluids and melts from subduction explain the distinction between high-Sr and low-Sr groups. We propose that mantle metasomatism by ancestral subduction and fluid-flux melting is heterogeneously distributed through the local subarc mantle and played a greater role in the genesis of the high-Sr basalt group. The limited range of primitive basalt types around the Goat Rocks volcanic cluster contrasts with the much greater diversity of basalts throughout the southern Washington to northern Oregon Cascade arc. On the other hand, the central volcanic cluster encompasses nearly the entire diversity observed at neighboring composite volcanoes. In the case of the Goat Rocks area at least, and perhaps attributable to the entire region, this means that the genesis of diverse intermediate magmas is independent from and does not require vastly different parental basalt compositions.

California, Oregon, Washington↗
USGS1996

Pesticides in surface water from three agricultural basins in south-central Georgia, 1993-97

Twenty-two of 43 pesticides analyzed were detected in 128 water samples collected from the Tucsawhatchee Creek, the Little River, and the Withlacoochee River. These streams drain agricultural basins in south-central Georgia and were sampled from March 1993 through June 1995. Herbicides were detected more frequently than insecticides. The most frequently detected herbicides were atrazine and metolachlor and the most frequently detected insecticide was carbaryl. Pesticide concentrations in the three streams were low and did not exceed U.S. Environmental Protection Agency drinking water standards. The maximum pesticide concentration was 2.6 mg/L (micrograms per liter) for propargite, a miticide detected in only one sample. The maximum concentrations of the remaining 21 pesticides were less than 0.25 mg/L. The median concentrations were equal to the method detection limit for all pesticides except atrazine (0.008 mg/L) and metolachlor (0.012 mg/L). The ratio of herbicide detections to nondetections was largest in the planting season, smaller in the harvest season and smallest in the fallow season for the three basins. The same pattern existed for the insecticide ratios in the Little River and the Withlacoochee River. Pairwise correlations between concentrations of atrazine and metolachlor, and four parameters (discharge, and concentrations of dissolved organic carbon, suspended organic carbon, and suspended sediment) were evaluated for each stream. The strongest correlations existed between metolachlor and mean daily discharge, and metolachlor and sediment in the Withlacoochee River. The only significant correlation for the Little River was between atrazine and suspended sediment.

Georgia↗
USGS2024

Hypothetical CO2 leakage into, and hydrological plume management within, an underground source of drinking water at a proposed CO2 storage facility, Kemper County, Mississippi, USA

A large Geologic Carbon Sequestration (GCS) hub has been proposed in Kemper County, Mississippi. The target injection interval consists of numerous Cretaceous-aged deep saline aquifers overlain by a competent and extensive regional sealing layer. Above the seal, the deepest Underground Source of Drinking Water (USDW) at the site is the Eutaw aquifer of the Eutaw Group and McShan Formation, undifferentiated. To assess potential risks of leakage from the deep sequestration reservoir, a model of a portion of the Cretaceous Eutaw Group was constructed in this study. Simulations tested various permeabilities, hypothetical leakage rates, and plume mitigation strategies utilizing existing wells. Results suggest that, under the influence of regional groundwater flow fields, leaking CO 2 would effectively bypass the existing wells, and to influence this migration would require very large water extraction rates. Therefore, to ensure plume detection, monitoring for leakage at the injection wells themselves is very important.

Mississippi↗
USGS2024

Analyzing spatial distributions and alignments of pitted cone features in Utopia Planitia on Mars

Martian geomorphology and surface features provide links to understanding past geologic processes such as fluid movement, local and regional tectonics, and feature formation mechanisms. Pitted cones are common features in the northern plains basins of Mars. They have been proposed to have formed from upwelling volatile-rich fluids, such as magma or water-sediment slurries. In this study, we map the spatial distributions of pitted cone features across the Utopia Planitia (UP) basin. Using the Average Nearest Neighbor technique, we find that pitted cone features appear to be clustered across the basin, occurring in a narrow band around the circumferential rim of UP. Additionally, we find that pitted cone features also appear to be aligned in chains that are sub-parallel with the basin margin. Parallel bands of cones generally follow elevation contours of the UP basin, which suggests elevation or a correlated factor plays a major role in pitted cone and cone chain formation. We propose that pitted cones and cone chains may be related to vertical fractures formed around the UP basin rim from subsidence of infilling basin material.

Icarus↗
USGS2025

Constraining large magnitude event source and path effects 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 events can be represented in non-ergodic GMMs. While we initially developed computation techniques using CyberShake simulations, the range of magnitudes and source-site combinations is not adequate to replicate what is observed empirically. We therefore designed a new ground motion simulation study, which includes earthquakes with a large 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 simulations (M4-M7), we then develop a non-ergodic GMM with the simulation data. We find that the within-site residuals are dominated by the 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 with the large events, while the second approach considers all small events on the fault plane. The results indicate that the path effects of large events cannot be satisfactorily approximated with that of small events using either approach.

Conference Paper↗
USGS2026

The effect of temperature, flow, density and disease on smallmouth bass (Micropterus dolomieu) early life growth rate

Individual growth rates of fish are influenced by gradually changing environmental conditions and rapid-onset events. However, there is a paucity of information on drivers of spatiotemporal variability in growth rates of fishes across large spatial extents. Using a 36-year dataset of smallmouth bass ( Micropterus dolomieu ) length-at-age data across 10 river reaches in Pennsylvania, USA ( n = 54,068 individuals), we estimated annual early life growth rates (mm year −1 ) and assessed the effects of relative abundance, summer temperature, summer streamflow rate and a disease outbreak on growth rate. We found that temperature had a positive effect on growth rate, varied spatially, and had a stronger effect in diseased reaches. There was also evidence for an impact of disease on early life growth through density-dependent mechanisms, with growth rates increasing from ∼85 to 101 mm·year −1 following disease-related decreases in abundance. This study adds to our understanding of the factors that shape early life growth, including rapid-onset events like disease which can have immediate and lasting effects on growth rate trajectory.

Pennsylvania↗
USGS1976

Transport and dispersion of fluorescent tracer particles for the dune-bed condition, Atrisco Feeder Canal near Bernalillo, New Mexico

A tracer technique in which mineral particles were coated with fluorescent dyes was used to study the rates of transport and dispersion of sediment particles`of various diameters and specific gravities for a dune-bed condition in an alluvial channel. The experiment was conducted in the Atrisco Feeder Canal near Bernalillo, N. Mex., between May 1 and July 14, 1967. A continuous point source of tracers, approximated by injections at 5- or 10-minute intervals, was maintained for 7 days so that the steady-dilution procedure could be used to calculate the transport rate of bed material. After termination of the injection process, the spatial-integration procedure was used to follow the movement of the tracers downstream and to calculate the transport rate. Samples of the bed material in transport and the accompanying tracers moving along the surface of the dune bed were obtained periodically throughout the study with the "dustpan" sampler especially designed for fluorescent tracer studies. In addition, the spatial distributions of the tracers in the dune bed were determined three times during the study by core sampling.

New Mexico↗
USGS2025

Reply to, “Comment on ‘The 1886 Charleston, South Carolina, earthquake: Relic railroad offset reveals rupture,’ by Roger Bilham and Susan E. Hough”

We welcome this opportunity to respond to Pratt et al. (2024) (hereinafter P24). Bilham and Hough (2023) proposed a “first-cut” elastic deformation model for the 1886 earthquake, a quantitative source model constrained by identified coseismic constraints. A key observation was the measurement of a lateral offset of a railroad line south of Summerville, leading to a model with predominately dextral slip and minor convergence, from which we concluded that active faulting had raised the Penholoway Marine Terrace >6 m since ∼770 ka. P24 questioned these constraints and proposed an alternative rupture model with predominantly reverse slip. This alternative model is neither consistent with coseismic constraints nor with other geophysical data. In a revised model presented here, we recognize that uplift of the Penholoway Terrace is confined to the eastern edge of the terrace, which we conclude results from active folding and tectonic transpression centered on the dextral fault that offset the railroad in 1886.

South Carolina↗
USGS2024

Will there be water? Climate change, housing needs, and future water demand in California

Climate change in California is expected to alter future water availability, impacting water supplies needed to support future housing growth and agriculture demand. In groundwater-dependent regions like California's Central Coast, new land-use related water demand and decreasing recharge is already stressing depleted groundwater basins. We developed a spatially explicit state-and-transition simulation model that integrates climate, land-use change, water demand, and groundwater gain-loss to examine the impact of future climate and land use change on groundwater balance and water demand in five counties along the Central Coast from 2010 to 2060. The model incorporated downscaled groundwater recharge projections based on a Warm/Wet and a Hot/Dry climate future from a spatially explicit hydrological process-based model. Two urbanization projections from a parcel-based, regional urban growth model representing 1) recent historical and 2) state-mandated housing growth projections were used as alternative spatial targets for future urban growth. Agricultural projections were based on recent historical trends from remote sensing data. Annual projected changes in groundwater balance were calculated as the difference between land-use related water demand, based on historical estimates, and climate-driven recharge plus agriculture return flows. Results indicate that future changes in climate-driven groundwater recharge, coupled with cumulative increases in agricultural water demand, result in overall declines in future groundwater balance, with a Hot/Dry future resulting in cumulative groundwater decline in all but Santa Cruz County. Cumulative declines by 2060 are especially prominent in San Luis Obispo (−2.9 to −5.1 Bm 3 ) and Monterey counties (−6.5 to −8.7 Bm 3 ), despite limited changes in agricultural water demand over the model period. These two counties show declining groundwater reserves in a Warm/Wet future as well, while San Benito and Santa Barbara County barely reach equilibrium. These results suggest future groundwater supplies may not be able to keep pace with regional demand and declining climate-driven recharge, resulting in a potential reduction in water security in the region. However, our county-scale projections showed new housing and associated water demand does not conflict with California's groundwater sustainability goals. Rather, future climate coupled with increasing agricultural groundwater demand may reduce water security in some counties, potentially limiting available groundwater supplies for new housing.

California↗
USGS2025

Gravity and magnetic surveys of the Skaergaard intrusion, East Greenland

Aeromagnetic and gravity surveys of the Skaergaard intrusion in East Greenland were carried out in July–August 1971 as part of a grant to the University of Oregon Center for Volcanology to refine the models of crystallization and differentiation of the intrusion, specifically to test whether the intrusion is underlain by dense rocks of a reservoir 20 kilometers (km) thick (referred to as a “hidden zone”). The Skaergaard intrusion is a source of platinum group elements that are critical mineral resources for many technologies, and because no new data have been collected these legacy datasets remain a valuable asset. The total-intensity aeromagnetic survey was flown in early July 1971 with a proton precession magnetometer at a constant barometric altitude of 1.5 km (5,000 feet) with a nominal line spacing of 1 km. Two gravimeters were used to acquire 168 stations of which 86 were at known altitudes (mainly sea level) and 82 had altitudes measured by altimetry in late July–August 1971. Finally, a north-south ground vertical-intensity magnetic traverse was completed across the intrusion together with collection of oriented hand specimens. The hand specimens were measured for remnant magnetization and density, along with density measurements of more specimens collected by expedition geologists for other purposes. The intrusion is composed of layered gabbro with extensive crystal fractionation that is dense and strongly reversely polarized. After terrain correction and standard Bouguer gravity reduction, the gravity anomaly dataset was corrected for all rock above sea level using the density measurements of the various zones of the intrusion and the topographic and geologic maps (variable density Bouguer gravity reduction). A large regional gradient in the gravity anomaly data was removed using orthogonal polynomial fitting to the gridded data. The zonal volumes of rock below sea level were calculated from the dipping polygonal layer gravity model of the intrusion below sea level and combined with elliptic cross–section cylinders for the various zones above sea level to approximate the original zonal volumes of the intrusion. The residual gravity anomaly of 18–20 milligals (mGal) was only about half of the expected anomaly if a large hidden zone proposed from petrologic considerations were present, and both two-dimensional and three-dimensional models imply that the exposed series of intrusion zones explain the gravity anomaly by their down-dip extension below sea level together with a small hidden-zone volume. A three-dimensional model of the exposed rocks and their down-dip extension below sea level also can account for the aeromagnetic anomaly with little or no requirement for hidden-zone rock. The middle and upper zone units of the intrusion contain the most magnetite and account for most of the aeromagnetic anomaly.

Skaergaard intrusion↗
USGS2024

Static and dynamic strain in the 1886 Charleston, South Carolina, earthquake

During the 1886 Mw 7.3 Charleston, South Carolina, earthquake, three railroads emanating from the city were exposed to severe shaking. Expansion joints in segmented railroad tracks are designed to allow railroad infrastructure to withstand a few parts in 10,000 of thermoelastic strain. We show that, in 1886, transient contractions exceeding this limiting value buckled rails, and transient extensions pulled rails apart. Calculated values for dynamic strain in the meizoseismal region are in reasonable agreement with those anticipated from the relation between strain and moment magnitude proposed by Barbour et al. (2021) and exceed estimated tectonic strain released by the earthquake by an order of magnitude. Almost all of the documented disturbances of railroad lines, including evidence for shortening of the rails, can thus be ascribed to the effects of dynamic strain changes, not static strain. Little or no damage to railroads was reported outside the estimated 10 −4 dynamic strain contour. The correspondence between 10 −3 and 2×10 −4 contours of dynamic strain and Mercalli intensity 9 and 8, anticipated from the dependence of each quantity on peak ground velocity, suggests it may be possible to use railroad damage to quantitatively estimate shaking intensity. At one location, near Rantowles, ≈20 km west of Charleston, a photograph of buckled track taken one day after the earthquake has been cited as evidence for shallow dextral slip and has long focused a search for a causal fault in this region. Photogrammetric analysis reveals that the buckle was caused by transient contraction of <10 cm with no dextral offset. Our results further weaken the evidence for faulting in the swamps and forests south of the Ashley River in 1886, hitherto motivated by the photograph and limited macroseismic evidence for high‐intensity shaking.

South Carolina↗
USGS2024

Improving crop-specific groundwater use estimation in the Mississippi Alluvial Plain: Implications for integrated remote sensing and machine learning approaches in data-scarce regions

Study region The Mississippi Alluvial Plain (MAP) in the United States (US). Study focus Understanding local-scale groundwater use, a critical component of the water budget, is necessary for implementing sustainable water management practices. The MAP is one of the most productive agricultural regions in the US and extracts more than 11 km 3 /year for irrigation activities. Consequently, groundwater-level declines in the MAP region pose a substantial challenge to water sustainability, and hence, we need reliable groundwater pumping monitoring solutions to manage this resource appropriately. New hydrological insights for the region We incorporate remote sensing datasets and machine learning to improve an existing lookup table-based model of groundwater use previously developed by the U.S. Geological Survey (USGS). Here, we employ Distributed Random Forests, an ensemble machine learning algorithm to predict annual and monthly groundwater use (2014–2020) throughout this region at 1-km resolution, using pumping data from existing flowmeters in the Mississippi Delta. Our model compares favorably with the existing USGS model, with higher R 2 (0.51 compared to 0.42 in the previous model), and lower root mean square error (RMSE) and mean absolute error (MAE)— 0.14 m and 0.09 m, respectively in our model, compared to 0.15 m and 0.1 m in the previous model. Therefore, this work advances our ability to predict groundwater use in regions with scarce or limited in-situ groundwater withdrawal data availability.

Journal of Hydrology Regional Studies↗
USGS1990

Effects of storm-water runoff on local ground-water quality, Clarksville, Tennessee

Storm-related water-quality data were collected at a drainage-well site and at a spring site in Clarksville, Tennessee, to define the effects of storm-water runoff on the quality of ground water in the area. A dye-trace test verified the direct hydraulic connection between the drainage well and Mobley Spring. Samples of storm run off and spring flow were collected at these sites for nine storms during the period February to October 1988. Water samples were collected also from Mobley Spring and two other springs and two observation wells in the area during dry-weather conditions to assess the general quality of ground water in an urban karst terrain. Evaluation of the effect of storm-water runoff on the quality of local ground water is complicated by the presence of other sources of contaminants in the area Concentrations and load for most major constituents were much smaller in storm-water runoff at the drainage well than in the discharge of Mobley Spring, indicating that much of the chemical constituent load discharged from the spring comes from sources other than the drainage well. However, for some of the minor constituents associated with roadway runoff (arsenic, copper, lead, organic carbon, and oil and grease), the drainage well contributed relatively large amounts of these constituents to local ground water during storms. The close correlation between concentrations of total organic carbon and concentrations of most trace metals at the drainage-well and Mobley Spring sites indicates that these constituents are transported together. Many trace metals were flushed early during each runoff event. Mean storm loads for copper, lead, zinc, and four nutrient species (total nitrogen, ammonia nitrogen, total phosphorus, and orthophosphorus) in storm-water runoff at the drainage-well site were lower than mean storm load predicted from an existing regression model. The overprediction by the model may be a result of the small size of the drainage area relative to the range of drainage areas used in the development of the models, or to the below-normal amounts of rainfall during the period of sampling for this investigation. Loads& in storm-water runoff for 22 constituents were extrapolated from sampled storms to total loads for the period February to October 1988. Calculated loads for trace metals for the period ranged from 0.030pound.s for cadmium to 12pound.s for strontium. Loads of the primary nutrients ranged from 0.97pounds for nitrite as nitrogen to 34pounds of organic nitrogen. Storm-water quality at the drainage-well and Mobley Spring sites was compared to background water quality of the local aquifer; as characterized by dry-weather samples from three springs and two observation wells in the Clarksville area. Concentrations of total-recoverable cadmium, chromium, copper, lead, and nickel were higher in many stormwater samples from both the drainage-well and Mobley Spring sites than in samples from any other site. In addition, concentrations of total organic carbon, methylene blue active substances, and total-recoverable oil and grease were generally higher in storm-water samples from the drainage-well site than in any ground-water sample. Densities of fecal coliform and fecal streptococcus bacteria and concentrations of total recoverable iron, manganese, and methylene blue active substances in storm samples from the drainage-well site exceeded the maximum contaminant levels listed in Tennessee’s drinking-water standards (1988) by as much as 2,500 and 5,500 colonies per 100 milliliters, and 2.7, 0.29, and 0.05 milligrams per liter, respectively. Densities of fecal coliform and fecal streptococcus bacteria and concentrations of total-recoverable iron, manganese, and lead in storm samples from Mobley Spring exceeded the maximum contaminant levels by as much as 500 and 4,500 colonies per 100 milliliters, and 18.7,0.65, and 0.02 milligrams per liter, respectively. For iron, manganese, and bacteria, these undesirable levels are not necessarily attributable to storm-water recharge, because concentrations of these constituents also exceeded drinking-water standards in one or more of the dry-weather samples from selected springs and observation wells in the area.

Tennessee↗
USGS2025

Relating systematic molecular and textural properties of graptolite pyrolyzed via gold tube hydrous pyrolysis: Implications for thermal proxies in lower Paleozoic marine shales

A series of gold tube pyrolysis experiments (72 h, 300–550 °C, 50 MPa) conducted on a graptolite-rich lower Paleozoic marine shale generated pyrolysis residues for a comprehensive evaluation of the molecular and structural variability of three types of graptolite periderm. Organic petrology, Raman spectroscopy, and field emission scanning electron microscopy (FE-SEM) with energy dispersive spectroscopy (EDS) were combined to evaluate the thermal evolution process. The three types of graptolite periderm, namely granular, non-granular, and nodular graptolite, were analyzed by Raman spectroscopy wherein point measurements were obtained after the maceral was identified and the location verified by organic petrology. Distinct thermal evolution pathways among non-granular, granular, and nodular graptolite periderms were recorded. The evolution patterns of the Raman parameters, particularly D1 and G bands, highlight the differences in geochemical composition of the graptolite periderm types and the alteration of molecular structure with increasing thermal maturity. Raman parameters D1 (position of the D1 peak), G-FWHM (full width at half maximum of the G peak), and ratios D1-FWHM/G-FWHM (full width at half maximum of the D1 peak ratioed to G-FWHM) and A D1 /A G (ratio of D1 and G peak intensities) showed effectiveness in assessing thermal maturity. Bireflectance with increasing gold tube pyrolysis temperature followed a hierarchy: non-granular > granular > nodular, reflecting different molecular alignment intensities. Qualitative FE-SEM evaluation showed that fine-grained mineral inclusions (primarily Fe-sulfide as determined via EDS) were associated with the graptolite populations, with granular graptolite containing greater amounts of coarser-grained (e.g., ∼300–1400 nm) mineral inclusions relative to non-granular and nodular graptolite, which contain finer-grained (e.g., ∼100–200 nm) inclusions difficult to resolve with optical microscopy. These findings are investigated to highlight the mechanisms that drive organic matter evolution within graptolite during thermal maturation, as well as to explore some of the limitations of using spectroscopic parameters as thermal maturity proxies.

International Journal of Coal Geology↗
USGS2025

REDPy: A Python tool for automated repeating earthquake detection and visualization

Detecting and cataloging seismic events are among the most fundamental tasks in seismology. Many standardized tools for these tasks exist, including the open‐source package repeating earthquake detector in Python (REDPy). REDPy generates an organized catalog of seismic events from continuous waveform data, in which events are automatically separated into groups (“families”) by their waveform similarity through cross‐correlation. REDPy also automatically generates various outputs that allow a user to visualize important trends in the catalog, which may be used in real time or in retrospective analyses to allow rapid identification of interesting features. The code was designed for near‐real‐time volcano monitoring but is applicable across a broad range of use cases in seismology and seismoacoustics. In this article, the utility and performance of REDPy are demonstrated on two highly seismogenic volcanic eruption sequences: the onset of the dome‐building eruption of Mount St. Helens, Washington, from 2004 to 2005, and the entirety of the summit caldera collapse sequence of Kīlauea, Hawai‘i, in 2018. This article is meant to be a companion to the documentation of the code; in addition to detailing the basic required inputs, script functionality, and resulting outputs, the reasonings behind several important design decisions are also discussed.

Seismological Research Letters↗
USGS2025

Divide and conquer: Separating the two probabilities in seismic phase picking

There are two fundamental probabilities in the seismic phase picking process – the probability of the existence of a seismic phase (detection probability) and the probability of correctly identifying the phase arrival time (timing probability). The nearly ubiquitous approach in developing deep learning phase picking models is to use a kernel, such as a truncated Gaussian, to mask the labeled phase arrival time, and train a segmentation model. Once a model is trained, the times of the peaks in the output are taken as phase arrival times (picks) and the height of the peaks are taken as “probability” of the picks. Here, we show that this “probability” represents neither the detection nor the timing probabilty because this approach forces the output to follow the shape of the kernel. We introduce an approach using two models to estimate these two distinct probabilities. We use a binary classifier with a calibrated confidence to address the detection probability and a multi-class classifier to obtain a probability mass function to address the timing probability. This new approach makes the deep learning-based phase picking process more interpretable and gives us options to logically control seismic monitoring workflows.

Geophysical Journal International↗
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↗