Geology ReportsSearch

SEARCH · Geology Reports

Results for “Elementa: Science of the Anthropocene”

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,687 records · Page 50Linked to original sources

Using the horizontal-to-vertical spectral ratio method to estimate thickness of the Barry Arm landslide, Prince William Sound, Alaska

Conducting detailed investigations of large landslides is difficult, especially in the subsurface, largely due to environmental factors such as steep slopes, difficult access, and numerous objective hazards. These factors have made it challenging to accurately estimate the depth to the failure surface of the Barry Arm landslide, a large (roughly 10 8 cubic meters), deep-seated bedrock landslide in Prince William Sound, Alaska, recognized in 2019. The landslide has exhibited accelerated movement in recent years and poses a potential tsunamigenic hazard if rapid failure occurs. Failure surface depth, equivalent to landslide thickness, is a necessary metric for landslide-volume calculations and associated tsunami wave models. In this report, we used seismic noise recorded by a seismometer located on the Barry Arm landslide in Alaska to calculate the horizontal-to-vertical spectral ratio (HVSR) to investigate the site fundamental frequency ( f 0 ) and depth of the failure surface. To ensure that observed peak frequencies in the spectral ratio were related to the underlying stratigraphy (and not caused by other noise sources like nearby glaciers, topographic resonance, weather, or human activities), we also calculated HVSRs using earthquake signals, HVSRs at other seismic stations within a 2.5-kilometer radius, and a standard spectral ratio between the landslide station and other sites. We observed multiple peaks in the landslide HVSR curves at 1.5 hertz (Hz), 4–5 Hz, and 7–11 Hz. The frequencies of these peaks were consistent at the landslide site through time and across methods and were dissimilar to those identified at other seismic stations in the area, making it unlikely the peaks were caused by local noise. Directional HVSRs calculated at 15-degree intervals showed amplification of the higher frequency peaks in the direction parallel to slip, indicating two-dimensional site effects. We used the distinct frequency peaks in the seismic record to develop a 4-layer conceptual model of the landslide wherein the top of the deepest layer represents the primary failure surface, or the boundary between damaged (mobile) and undamaged material. We inverted Rayleigh wave ellipticity curves within this 4-layer configuration with constraints on S-wave velocity and layer thickness based on analogous material properties identified in the literature. This was necessary absent any site-specific subsurface S-wave velocity data. The best-fitting models indicate a mean slope-normal depth to the failure surface of 188 (±9) meters (m), with additional stratigraphic boundaries at 4 and 20 m below ground surface, potentially representing layered motion. These results agree with and improve upon ranges estimated by previous studies and can support future modeling and assessment efforts at Barry Arm.

Alaska

The advantages of electric-grade heat storage in long closed-loop wells

Though long (e.g., tens of kilometers) closed-loop geothermal wells are under consideration for the extraction of subsurface heat, these wells might also serve as an efficient energy storage mechanism for electricity generation. Using a semi-analytic model solution, the potential for electric-grade heat storage as a function of ambient temperature (e.g., corresponds to depth of loop), well length, well diameter, and flow rate is evaluated. The following simplified cases for comparison of standard closed loop operation with energy storage operation were considered: [1] constant flow rate and constant temperature (90 °C) injection to represent the standard closed-loop base-case; and [2] constant flow rate and annual cycle sinusoidal temperature (90-150 °C) to represent seasonal (summer) charging while solar resources are peak. Initial temperatures for all scenarios considered herein are a uniform 175 °C. Electric-grade heat is assumed to be delivered whenever temperatures at the extraction point exceed 90 °C. For calculation purposes only, if temperature falls below 100 °C, it is assumed that heat delivered is sub-economic, so no electricity would be produced. For all scenarios, temperatures at the extraction well asymptotically approach the flow-weighted average injection temperature, but energy storage scenarios exhibited a damped time-varying signal that diminishes in magnitude with length of the loop. The asymptotic approach depends on initial temperatures in the rock and the heat extraction rate (a function of well diameter and flowrate). This analysis demonstrates that shorter closed loops can produce more electricity over time than longer closed-loops previously proposed for electricity production over typical engineering design lifetimes (e.g., 30 years). Although only a high-temperature scenario is considered herein, rock that is initially below boiling temperature would not host a standard closed-loop resource, but injection of hot water seasonally would asymptotically heat this low-temperature system to temperatures capable of electricity production. In other words, regardless of initial temperatures, closed loops could be used to store electricity with no critical minerals in the geothermal battery.

Conference Paper

Biogeochemical and community ecology responses to the wetting of non-perennial streams

Transitions between dry and wet hydrologic states are the defining characteristic of non-perennial rivers and streams, which constitute the majority of the global river network. Although past work has focused on stream drying characteristics, there has been less focus on how hydrology, ecology and biogeochemistry respond and interact during stream wetting. Wetting mechanisms are highly variable and can range from dramatic floods and debris flows to gradual saturation by upwelling groundwater. This variation in wetting affects ecological and biogeochemical functions, including nutrient processing, sediment transport and the assembly of biotic communities. Here we synthesize evidence describing the hydrological mechanisms underpinning different types of wetting regimes, the associated biogeochemical and organismal responses, and the potential scientific and management implications for downstream ecosystems. This combined multidisciplinary understanding of wetting dynamics in non-perennial streams will be key to predicting and managing for the effects of climate change on non-perennial ecosystems.

Nature Water

Biodegradation of chlorinated ethenes at a karst site in middle Tennessee

This report presents results of field and laboratory investigations examining the biodegradation of chlorinated ethenes in a karst aquifer contaminated with trichloroethylene (TCE). The study site, located in Middle Tennessee, was selected because of the presence of TCE degradation byproducts in the karst aquifer and available site hydrologic and chlorinated-ethene information. Additional chemical, biological, and hydrologic data were gathered to evaluate whether the occurrence of TCE degradation byproducts in the karst aquifer was the result of biodegradation within the aquifer or simply transport into the aquifer. Geochemical analysis established that sulfate-reducing conditions, essential for reductive dechlorination of chlorinated solvents, existed in parts of the contaminated karst aquifer. Other areas of the aquifer fluctuated between anaerobic and aerobic conditions and contained compounds associated with cometabolism, such as ethane, methane, ammonia, and dissolved oxygen. A large, diverse bacteria population inhabits the contaminated aquifer. Bacteria known to biodegrade TCE and other chlorinated solvents, such as sulfate-reducers, methanotrophs, and ammonia-oxidizers, were identified from karst-aquifer water using the RNA-hybridization technique. Results from microcosms using raw karst-aquifer water found that aerobic cometabolism and anaerobic reductive-dechlorination degradation processes were possible when appropriate conditions were established in the microcosms. These chemical and biological results provide circumstantial evidence that several biodegradation processes are active in the aquifer. Additional site hydrologic information was developed to determine if appropriate conditions persist long enough in the karst aquifer for these biodegradation processes to be significant. Continuous monitoring devices placed in four wells during the spring of 1998 indicated that pH, specific conductance, dissolved oxygen, and oxidation-reduction potentials changed very little in areas isolated from active ground-water flow paths. These stable areas in the karst aquifer had geochemical conditions and bacteria conducive to reductive dechlorination of chlorinated ethenes. Other areas of the karst aquifer were associated with active ground-water flow paths and fluctuated between anaerobic and aerobic conditions in response to rain events. Associated with this dynamic environment were bacteria and geochemical conditions conducive to cometabolism. In summary, multiple lines of evidence developed from chemical, biological, and hydrologic data demonstrate that a variety of biodegradation processes are active in this karst aquifer.

Tennessee

Archean to Mesozoic–Cenozoic seismic crustal structure: Implications for geological and biological evolution

We use >4500 measurements of crustal structure to investigate the seismic structure of continental crust, Archean to Mesozoic–Cenozoic. The mean crustal thickness of continents, including their margins, is 36.5 km. We find that Archean, Paleoproterozoic, and Mesoproterozoic crust have similar mean crustal thickness (41 km), seismic velocities, Vp/Vs ratio, and density. Crusts of these ages span 3 Ga, from 4.0 to 1.0 Ga, and their similarity of physical properties (thickness, density, and seismic velocities) suggests that the process of crustal formation may have been similar during this time period, which covers 66% of Earth's history. Neoproterozoic crustal properties differ significantly from Mesoproterozoic and older crusts. The mean thickness of Neoproterozoic crust is 32 km, some 9 km thinner than the mean thickness of Archean, Paleoproterozoic, and Mesoproterozoic crust. A lithospheric root with a thickness of 150–200 km underlies Archean, Paleoproterozoic, and Mesoproterozoic crusts, and this root resists lithospheric rifting and crustal extension and thinning. Based on these observations, particularly the lithospheric thickness, we conclude that Archean, Paleoproterozoic, and Mesoproterozoic lithosphere are unique and together form the stable cratonic nuclei, defined as the thick (150–250 km), long-lived cores of continents. Higher mantle temperatures during the Archean, Paleoproterozoic, and Mesoproterozoic may have played a key role in the formation of the thick lithospheric roots. The second key finding is that the similarities of mean seismic properties indicate that the process of crustal formation operating in the Archean continued during the Paleoproterozoic and Mesoproterozoic. The thin (mean value < 135 km) lithospheric root beneath Neoproterozoic and younger crust may be related to the steady decrease in mantle temperature through time. Neoproterozoic and Paleozoic crust have similar physical properties, and these eras are characterized by pronounced biodiversification, including the renowned Garden of Ediacara, the Cambrian Explosion, and the Great Ordovician Biodiversification. Mesozoic–Cenozoic crust is the most diverse and reflects the current tectonic and magmatic processes of crustal formation.

Book chapter

Deep lithospheric controls on surface deformation and seismicity around the East Anatolian Fault Zone and A3 Triple Junction

The East Anatolian Fault Zone (EAFZ) is a plate-bounding strike-slip fault capable of hosting large earthquakes, as demonstrated by the extremely damaging February 2023 M w 7.8 and M w 7.7 mainshocks of the Kahramanmaraş earthquake sequence. Deformation related to this boundary, part of the Anatolia-Arabia-Africa (A 3 ) Triple Junction, is diffuse, as was shown by part of this earthquake sequence occurring on a northern splay of the EAFZ (the Sürgü-Çardak Fault Zone; SCFZ). Controls on surface deformation are commonly linked to stress in the brittle upper crust, but the complex deformation and seismicity patterns in this region may also reflect deeper processes, such as variations in the location and extent of the strong Arabian Plate lithospheric mantle. Seismic tomography indicates that the Arabian Plate underthrusts Anatolia as far north as the SCFZ and extends as far west as the central Adana Basin, coincident with a zone of relatively deep (>30 km) strike-slip seismogenesis that has produced M w >6 earthquakes. By investigating the relationship between deformation since the inception of the EAFZ (ca. 5 Ma), seismic structure, and seismicity, we infer that the SCFZ will become the future SE boundary of the Anatolian Plate as part of the evolving A 3 Triple Junction.

GSA Today

Sensitivity of benthic biota and toxicity of cadmium, cobalt, copper, nickel, lead, and zinc mixtures in Washington, United States, and British Columbia, Canada

Relative sensitivities and responses of juvenile white sturgeon ( Acipenser transmontanus ), Hyalella azteca , two families of mayfly (Ephemerellidae, Heptageniidae), one family of caddisfly (Brachycentridae), and a natural community of benthic macroinvertebrates (BMI) to multiple metals are predicted using previously collected laboratory and field samples and a metal mixture model. Biological responses in single metal exposures are used to parameterize toxicity functions, which include accumulations of hydrogen and selected metals on biological receptors, intrinsic potencies of hydrogen and metals, sensitivities of organisms, and times of exposure. The model then is used to predict responses in multiple metal laboratory exposures and field-collected porewater. The following sensitivity sequence in porewater was determined based on endpoints of survival or total abundance: juvenile white sturgeon greater than (>) Ephemerellidae family > Hyalella azteca > Heptageniidae family about equal to (~) benthic macroinvertebrate community > Brachycentridae family. The fraction of porewater samples that are predicted to have adverse impacts on benthic biota (20-percent or greater negative response) depends on organism sensitivities and metal toxicities, and ranges from 44 to 48 percent for juvenile white sturgeon, 23 to 26 percent for the Ephemerellidae family, 16 to 22 percent for Hyalella azteca , 5 to 8 percent for the Heptageniidae family and BMI community, and 0 percent for the caddisfly family. The most toxic porewater in the upper Columbia River Basin (UCR) is at the backwater bar site at Deadmans Eddy and China Bend. The model also indicates that the element responsible for the most toxic conditions in UCR porewater is copper for all organisms, except Hyalella azteca and the metal-insensitive Brachycentridae family. Copper and lead result in the most toxic conditions for Hyalella azteca . This approach and results can aid in assessing metal toxicity and its potential risk to aquatic biota in ecosystems impacted by historical mining activities.

British Columbia, Washington

Integration of rupture directivity models for the US National Seismic Hazard Model

Several rupture directivity models (DMs) have been developed in recent years to describe the near-source spatial variations in ground motion amplitudes related to propagation of rupture along the fault. We recently organized an effort towards incorporating these directivity effects into the USGS National Seismic Hazard Model (NSHM), by first evaluating the community's work and potential methods to implement directivity adjustments into probabilistic seismic hazard analysis (PSHA). Guided by this evaluation and comparison among the considered DMs, we selected an approach that can be readily implemented into the USGS hazard software, that provides an azimuthally varying adjustment to the median ground motion and its aleatory variability. This method allows assessment of the impact on hazard levels and provides a platform to test the DM amplification predictions using a generalized coordinate system, necessary for consistent calculation of source-to-site distance terms for complex ruptures. We give examples of the directivity-related impact on hazard, progressing from a simple, hypothetical rupture, to more complex fault systems, composed of multiple rupture segments and sources. The directivity adjustments were constrained to strike-slip faulting, where DMs have good agreement. We find that rupture directivity adjustments using a simple median and aleatory adjustment approach can impact hazard both from a site perspective and on a regional scale, increasing shaking off the end of the fault trace up to 30--40\% and potentially reducing it for sites along strike. Statewide hazard maps of California show that the change in shaking along major faults can be a factor to consider for assessing long-period (>ls) near-source effects within the USGS NSHM going forward, reaching up to 10--20\%. Finally, we suggest consideration of minimum parameter ranges and baseline requirements as future DMs are developed to minimize single approach adaptations, to enable more consistent application within both ground motion and hazard studies.

Earthquake Spectra

Wangyanite, PdNi8S8, a new Pd end-member mineral of the pentlandite group from the J-M reef, Stillwater Complex, Montana, USA

Wangyanite (IMA2024-008a), ideally PdNi 8 S 8 , is a Pd end-member mineral of the pentlandite group that was discovered in the J-M reef of the Stillwater Complex, Montana, USA. Wangyanite occurs as anhedral-subhedral granular crystals 200–400 µm in size, associated with isoferroplatinum, braggite, pentlandite, and chalcopyrite interstitial to plagioclase grains within anorthosite. Wangyanite exhibits a yellowish brown color with a black streak and a metallic luster. It is brittle with uneven fractures, and has a calculated density of 5.14 g/cm 3 . The mineral does not show discernible pleochroism, bireflectance, or anisotropy. It has an average composition of 9.95 wt.% Pd, 31.95 wt.% Ni, 25.02 wt.% Fe, 0.57 wt.% Co, 31.74 wt.% S, totaling 99.23 wt.%. The empirical formula, based on eight sulfur atoms per formula unit, is (Pd 0.76 Co 0.08 ) Σ0.84 (Ni 4.39 Fe 3.60 ) Σ7.99 S 8 . Wangyanite has a cubic cell with a space group of Fm -3 m (#225), having lattice parameters of a = 10.1167(12) Å, V = 1035.4(4) Å 3 , and Z = 4. Its crystal structure has been solved by single-crystal three-dimensional electron diffraction study. The strongest X-ray diffraction lines of wangyanite are claculated at [ d in Å ( I %)( hkl )]: 5.841(14.03)(111), 3.050(100)(311), 1.947(29.16)(115,333), 1.264(11.66)(800), 3.577(8.79)(220), 2.920(20.82)(222), and 2.321(9.34)(331). Wangyanite shares the same crystal structure as pentlandite, but the octahedrally coordinated site is mainly occupied by Pd in wangyanite. Based on the textural features and previous experimental Pd-Fe-Ni-S phase system, wangyanite could form by peritectic reaction between braggite, pentlandite and sulfide liquid. These mineral associations are stable in a Ni-Pd-rich sulfide melt system at about 550 °C or even lower temperature. Therefore, wangyanite can potentially serve as an indicator of the presence of Pd-rich residual melts. The mineral is named in honor of Prof. Christina Yan Wang, a well-known researcher on platinum-group element (PGE) occurrences and enrichment mechanisms in mafic-ultramafic intrusions, notably those deposits related to the Emeishan large igneous province in China.

Montana

Exploratory water well, St. George Island, Florida

A 1,026-foot well was drilled into the Floridan aquifer on St. George Island in search of a supply of water for residential development in 1971. This is the only well on the island on record to have been drilled below 200 feet for water. The well was abandoned; water flowing from the well was saline. It had a specific conductance of 7,500 micromhos per centimeter at 25° Celsius. The specific conductances of samples taken at drilling depths 217-1,026 feet ranged from 3,140 to 39,600, with the maximum value at 470 ft. A composite geophysical and lithologic log of the well shows spontaneous potential, single-point resistance, gamma-ray, neutron, gamma-gamma, and caliper logs. Well-site notes, a written description of cuttings samples, a graph showing the relation of specific conductance of water samples to the depth of the well at time of sampling, and four chemical analyses of the water are included in the report.

Florida

High-resolution magnetic survey using an unoccupied aerial vehicle to constrain buried lava flow geometry, volume, and eruptive history of Little Cones, Crater Flat, Nevada

Magnetic surveys are an important tool used to augment geologic mapping in distributed volcanic fields. Using magnetic anomalies, it is possible to model the geometry of shallowly buried volcanic features, such as conduits, sills, and lava flows. This subsurface mapping is important for understanding eruption dynamics and emplacement of lava flows, and it sometimes reveals buried volcanoes no longer visible at the surface. These data are critical to better interpret the numbers, styles, and magnitudes of eruptions in distributed volcanic fields and their associated volcanic hazards. New advances in unoccupied aerial vehicles (UAVs) offer an attractive middle range of resolution and aerial coverage between ground-based magnetic surveys and aeromagnetic surveys. Here, we present the results of a UAV fluxgate magnetic survey of the Little Cones, Nevada, scoria cones, which have been the target of previous ground and aeromagnetic surveys. The magnetic anomalies at Little Cones are of interest because the surrounding alluvium conceals lava flows that erupted from Little Cones, making it very difficult to understand the volume and morphology of lava flows from geologic mapping alone. Nonlinear inversion of UAV-collected magnetic data were used to model the thickness and morphology of buried Little Cones’ lava flows with higher precision than achieved previously. The sequence of events and calculated flow characteristics are then interpreted. The total volume of Little Cones, including concealed lava flows, is approximately 0.016 cubic kilometer, and the initial sheet flow erupted in less than 24 hours. The findings presented herein demonstrate that UAV-based magnetic surveys are a reliable method of data collection and an efficient alternative to other survey methods, facilitating development of a three-dimensional perspective of distributed volcanic fields.

Nevada

The role of bedrock circulation depth and porosity in mountain streamflow response to prolonged drought

Quantitative understanding is lacking on how the depth of active groundwater circulation in bedrock affects mountain streamflow response to a multi‐year drought. We use an integrated hydrological model to explore the sensitivity of a variety of streamflow metrics to bedrock circulation depth and porosity under a plausible extreme drought scenario lasting up to 5 years. Endmember depth versus hydraulic conductivity relationships and porosity values for fractured crystalline rock are simulated. With drought, a deeper circulation system with higher drainable porosity more effectively buffers minimum flow and significantly limits perennial stream loss in comparison to a shallow circulation system. Streamflow buffering is accomplished through extensive groundwater storage loss. However, deeper circulation systems experience prolonged recovery from drought in comparison to storage‐limited shallow systems. Research highlights the importance of characterizing the deeper bedrock hydrogeology in mountainous watersheds to better understand and predict drought impacts on stream ecosystem health and water resource sustainability.

Colorado

Polystyrene microplastics alter the accumulation and elimination dynamics of silver nanoparticles in Daphnia magna

Rationale Microplastics (MPs) can interact with engineered nanomaterials and alter their environmental fate and bioavailability. However, their influence on the bioaccumulation dynamics of silver nanoparticles (AgNPs) in aquatic filter feeders under environmentally relevant conditions remains poorly understood. We hypothesized that polystyrene microplastics (PS-MPs) alter the uptake, elimination and overall bioaccumulation dynamics of AgNPs in Daphnia magna . Methodology Adult D. magna were exposed to isotopically labeled citrate-coated 109 AgNPs in the presence and absence of PS-MPs. Waterborne uptake, dietary uptake, assimilation efficiency, food ingestion and elimination were quantified experimentally and incorporated into a biodynamic model to predict steady-state silver (Ag) concentrations under environmentally relevant exposure scenarios. Results PS-MPs increased waterborne Ag accumulation by approximately eightfold compared with AgNPs exposure alone. However, Ag elimination was substantially faster, with 98% of accumulated Ag eliminated after five days of depuration compared with 66% without PS-MPs. PS-MPs also reduced food ingestion rates (IRs), while assimilation efficiency remained largely unchanged. Biodynamic modeling predicted that steady-state Ag accumulation was approximately fivefold greater in the presence of PS-MPs, with waterborne exposure becoming the dominant accumulation pathway. Discussion These findings suggest that PS-MPs alter bioaccumulation dynamics of AgNPs in filter-feeding organisms by enhancing organism-associated Ag during waterborne exposure while accelerating Ag elimination, likely through particle-associated transport and gut egestion. Collectively, these results underscore the complexity of Ag bioaccumulation processes when MPs and AgNPs interact and highlight the importance of studies conducted under environmentally relevant conditions.

Environmental Chemistry

Surface-wave relocation and characterization of the October 2023 tsunamigenic seismic unrest near Sofugan volcano, Izu Islands, Japan

A moderate-magnitude earthquake swarm occurred in the remote Izu Islands region of Japan between October 1 and 8, 2023. The swarm included 151 shallow earthquakes cataloged by the U.S. Geological Survey, which notably included a roughly 2.5-hr episode of 15 successive magnitude (M) < 5.5 earthquakes. Origin times were coincident with regionally recorded tsunami waves, but tsunamigenesis for moderate-magnitude earthquakes is uncommon, indicating that volcanic activity generated the ocean displacements. Leveraging a surface-wave relative relocation approach, we estimate precise epicentroid locations for the remote swarm. Final epicentroids and caldera analogs indicate a three-stage model to explain swarm activity: (a) caldera pressurization due to magma intrusion, (b) depressurization via dike propagation away from the caldera, and (c) eruption corresponding with caldera reactivation either by collapse or additional intrusion.

Sofugan volcano, Izu Islands

Complex hydrology and variability of nitrogen sources in a karst watershed

Streams draining karst areas with rapid groundwater transit times may respond relatively quickly to nitrogen reduction strategies, but the complex hydrologic network of interconnected sinkholes and springs is challenging for determining the placement and effectiveness of management practices. This study aims to inform nitrogen reduction strategies in a representative agricultural karst setting of the Chesapeake Bay watershed (Fishing Creek watershed, Pennsylvania) with known elevated nitrate contamination and a previous documented groundwater residence time of less than a decade. During baseflow conditions, streamflow did not increase with drainage area. Headwaters and the main stem lost substantial flow to sinkholes until eventually discharging along large springs downstream. Seasonal hydrologic conditions shift the flow and nitrogen load spatially among losing and gaining stream sections. A compilation of nitrogen source inputs with the geochemistry and the pattern of enrichment of δ 15 N and δ 18 O suggest that the nitrogen in streams and springs during baseflow represents a mixture of manure, fertilizer, and wastewater sources with low potential for denitrification. The pH and calcite saturation index increased along generalized flow paths from headwaters to springs and indicate shorter groundwater residence times in baseflow during the spring versus summer. Given the substantial investment in management practices, fixed monitoring sites could incorporate synoptic water sampling to properly monitor long-term progress and help inform management actions in karst watersheds. Although karst watersheds have the potential to respond to nitrogen reduction strategies due to shorter groundwater residence times, high nitrogen inputs, effectiveness of conservation practices, and release of legacy nutrients within the karst cavities could confound progress of water quality goals.

Pennsylvania

Stream habitat characteristics at selected sites in the Georgia-Florida coastal plain

Habitat characterization is part of a multidisciplinary approach to water-quality assessment implemented by the National Water-Quality Assessment Program. Habitat data were collected in the Georgia-Florida Coastal Plain study unit at 24 sites during 1993-95. Data were collected for habitat characteristics at three spatial scales: basin, segment, and reach. Basin data include physiography, land resource provinces, and land use, providing a description of the environmental setting at each site. Segment data include length, gradient, and sinuosity. A Kendall correlation analysis performed on segment characteristics and the log-of-basin area showed a correlation between segment gradient and the log-of-basin area and a correlation between sinuosity and segment length. Reach data consist of field-collected measurements of both instream and riparian habitats. Sand and detritus were the most common channel-bed substrates among the sampled sites. Measurements of channel width, water depth, and bank width and height were used to create cross-sectional profiles of each sampled area. Elevations of selected durations plotted on cross sections illustrated the percentage of time that the banks were inundated at each site. Sites were divided into two groups based on duration of bank inundation (less than or equal to 1 percent and greater than 1 percent). Bank woody vegetation was also sampled and a clustering algorithm known as Two-Way INdicator SPecies ANalysis (TWINSPAN) was used to analyze these data. TWINSPAN divided the sites into two groups based on their vegetation composition. A statistical comparison of the two types of site groups (duration of bank inundation and vegetation) was performed. The significant association between these groups was consistent with the hypothesis that inundation frequency affected riparian vegetation.

Florida, Georgia

Don’t Let Negatives Hold You Back: Accounting for Underlying Physics and Natural Distributions of Hydrothermal Systems When Selecting Negative Training Sites Leads to Better Machine Learning Predictions

Selecting negative training sites is an important challenge to resolve when utilizing machine learning (ML) for predicting hydrothermal resource favorability because ideal models would discriminate between hydrothermal systems (positives) and all types of locations without hydrothermal systems (negatives). The Nevada Machine Learning project (NVML) fit an artificial neural network to identify areas favorable for hydrothermal systems by selecting 62 negative sites where the research team had confidence that no hydrothermal resource exists. Herein, we compare the implications of the expert selection of negatives (i.e., the NVML strategy) with a random sample strategy, where it is assumed that areas outside the favorable structural ellipses defined by NVML are negative. Because hydrothermal systems are sparse, it is highly probable that, in the absence of a favorable geological structure, hydrothermal favorability is low. We compare three training strategies: 1) the positive and negative labeled examples from NVML; 2) the positive examples from NVML with randomly selected negatives in equal frequency as NVML; and 3) the positive examples from NVML with randomly selected negatives reflecting the expected natural distribution of hydrothermal systems relative to the total area. We apply these training strategies to the NVML feature data (input data) using two ML algorithms (XGBoost and logistic regression) to create six favorability maps for hydrothermal resources. When accounting for the expected natural distribution of hydrothermal systems, we find that XGBoost performs better than the NVML neural network and its negatives. Model validation was less reliable using F1 scores, a common performance metric, than comparing probability estimates at known positives, likely because of the extreme natural class imbalance and the lack of negatively labeled sites. This work demonstrates that expert selection of negatives for training in NVML likely imparted modeling bias. Accounting for the sparsity of hydrothermal systems and all the types of locations without hydrothermal systems allows us to create better models for predicting hydrothermal resource favorability.

Geothermal Resources Council Transactions

Suitability of ground-motion models for seismic hazard assessment in Puerto Rico and the U.S. Virgin Islands

We perform linear mixed-effects analyses with a ground-motion dataset to evaluate how well ground-motion models (GMMs) fit active crustal, subduction interface, and subduction intraslab earthquakes in Puerto Rico and the US Virgin Islands. Most of the GMMs reproduce the spatial variation in peak ground motions with earthquake magnitude and rupture distance but predict ground motions 0.3–1.0 natural log units (35%–270%) greater than observations. Two GMMs developed for Puerto Rico that are based on ground-motion records from mostly magnitude 4–5 earthquakes do not perform as well. We attribute the period-dependent overprediction in observed ground motions to differences between the observed site response and the linear site response in the GMMs. Consequently, we developed region-specific GMMs by adjusting the period-dependent linear site response coefficients and period-dependent constant coefficients to remove most of the bias between predicted ground motions and observations. For the analysis, we compile ground-motion records and process waveforms to build a dataset with 10,127 records at 72 stations from 849 magnitude 4.0–6.4 earthquakes between 1 January 2006 and 31 March 2024. The earthquakes include active crustal, subduction interface, and subduction intraslab events. We evaluate the GMMs using the time-averaged shear wave speed in the top 30 m ( ), which we compile from site surveys and proxy values computed from horizontal to vertical spectral ratios. Site terms exhibit strong consistency across GMMs and crustal and subduction earthquakes, indicating that the linear mixed-effects analysis successfully isolates the effects of local site response. The event terms show little spatial correlation and more substantial variability than in other regions, which we attribute to uncertainties in the earthquake magnitudes. This analysis guides the selection of the GMMs for the 2025 update of the National Seismic Hazard Model for Puerto Rico and the US Virgin Islands.

Puerto Rico, U.S. Virgin Islands