Geology ReportsSearch

SEARCH · Geology Reports

Results for “Ecology of Freshwater Fish”

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

Cave records reveal recent origin of North America’s deepest canyon

We explore how and when Hells Canyon, North America’s deepest river gorge (~2,400 m deep), formed, addressing these fundamental questions first posed by W. Lindgren [ The Gold Belt of the Blue Mountains of Oregon (1901)]. Existing hypotheses about the canyon’s formation and timing of incision remain speculative due to a lack of direct constraints and geomorphic analysis in the canyon. Herein, we combine cosmogenic nuclide dating of cave-bound river deposits, river profile analysis, and numerical modeling to provide the first direct age constraints and systematic analysis of incision processes at work in Hells Canyon. Our study reveals a significant drainage capture triggered rapid incision at ~2.1 ± 1.0 Ma, establishing the Snake River’s modern route into the Columbia River system. The increased drainage area and subsequent increase in stream power resulted in the rapid incision of Hells Canyon and the formation of tributary knickpoints (KPs) that decrease in elevation away from the capture location. Cosmogenic dating of cave deposits indicates incision rates increased from ~0.01 to ~0.16 mm y −1 . Numerical modeling of the stream capture supports these observations, demonstrating how abrupt drainage area increase drives rapid river incision. Our findings from Hells Canyon provide a well-constrained example of how drainage capture can dramatically shape the evolution of a major river gorge.

Idaho, Oregon, Washington

Geophysical architecture and geochronology of the Neoarchean Mentor anorthosite intrusive complex, northwestern Minnesota: Largest anorthosite complex of the Superior Province?

The Neoarchean Mentor anorthosite intrusive complex (MAIC) lies within the Wawa subprovince in northwestern Minnesota, in a region where the Wawa, Quetico, and Wabigoon subprovinces are juxtaposed in close proximity. Archean rocks are entirely concealed, and interpretations are developed from aeromagnetic, gravity, and borehole samples. The MAIC includes both anorthosite (dense, weakly magnetized) and oxide-rich gabbro (strongly magnetized) lithologies. Anorthosite is coarse-grained to megacrystic, intensely altered, and locally brecciated. Pervasive epidote alteration enhanced the density of the anorthosite via introduction of higher-density mineral assemblages, explaining why the MAIC produces a significant gravity high. Oxide-rich gabbro forms a border phase of the MAIC and has potential for vanadium, chromium, and titanium mineralization, and produces a strong aeromagnetic high. The MAIC is interpreted to extend over an area of 640 km 2 , making it the largest known anorthosite complex of the Superior Province, as measured by preserved areal extent. Modeling indicates the MAIC extends more than 6 km into the subsurface. A new Pb–Pb zircon age of 2737.2 ± 4.5 Ma is interpretated as the crystallization age of anorthosite within the MAIC, showing that the MAIC formed well before the ca. 2690 Ma Shebandowanian orogeny, and raising new questions about correlations with other parts of the Wawa subprovince. A low-density batholith, here informally called the Fertile batholith, is interpreted to intrude the southern part of the MAIC. A new Pb–Pb zircon age of 2701.1 ± 6 Ma is interpreted as the magmatic age of the Fertile batholith.

Minnesota

Geologic framework and hydrostratigraphy of the Edwards and Trinity aquifers within Hays County, Texas

During 2023–24, the U.S. Geological Survey, in cooperation with the Edwards Aquifer Authority, revised a previous publication of the geologic framework and hydrostratigraphy of the Edwards and Trinity aquifers that was completed during 2018 within Hays County, Texas. The purpose of this report is to present the updated geologic framework and hydrostratigraphy of the rocks containing the Edwards and Trinity aquifers in Hays County from field observations of the surficial expressions of the rocks. The report includes a detailed 1:24,000-scale hydrostratigraphic map with names and descriptions of the geologic framework and hydrostratigraphic units (HSUs) in the study area. The study includes updates to the interpretation of the Kainer Formation of the Edwards Group with the addition of a burrowed unit between the basal nodular and dolomitic members. Hydrostratigraphy was also updated with the addition of the Seco Pass HSU for the burrowed member of the Kainer Formation. The study also includes updates to the interpretation of the hydrostratigraphy of the Trinity aquifer with the addition of the cavernous HSU at the top of the upper zone of the Trinity aquifer and the Herff Falls HSU between the Bulverde and Rust HSUs of the middle zone of the Trinity aquifer. This updated report provides additional information about a complex aquifer system. The complexity in the aquifer system results from a combination of the original depositional history, bioturbation, development of primary and secondary porosity, postdepositional diagenesis, fracturing, and faulting.

Texas

Sedimentological and geochemical characterization of lacustrine deposits of the Babouri-Figuil basin, northern Cameroon: Implications for source rocks distribution and petroleum exploration

The West and Central African Rift System (WCARS) refers to the series of Cretaceous rift basins where commercial hydrocarbon accumulations have been discovered. Some of the WCARS frontier basins are currently being investigated to increase our understanding of these basins in light of new commercial discoveries. The present study was performed in the Babouri-Figuil Basin (BFB), which is genetically related to the WCARS and constitutes an area of interest in terms of petroleum prospecting, where the distribution of petroleum source rocks and potential targets for petroleum exploration across the entire basin is poorly understood. For the current study, an integrated facies analysis along with organic and inorganic geochemical techniques were applied to the basin's Cretaceous deposits with the aim of reconstructing the paleodepositional environment, assessing factors that triggered the input of organic matter, and providing a spatial overview of the organic matter accumulation in the basin based on outcrop samples. An alluvial fan-lacustrine-braided river system is inferred from the facies analysis of the stratigraphic sequence consisting of conglomerate, sandstone, siltstone, limestone, marlstone, and claystone. Bulk analysis of organic matter reveals that black shale and massive claystone are the main prospective petroleum source rocks in the basin. Inorganic geochemical analyses reveal the influence of anoxic conditions, moderate to high primary productivity, and low terrigenous inputs in organic matter enrichment. The formations rich in organic matter are predominantly concentrated in the western and eastern parts of the basin which may represent areas with depressions, characterized by high accommodation space. In terms of the regional context of the WCARS rift basins, typical hydrocarbon exploration in the BFB may target basal-conglomerate, sandstone beds situated directly above or/and interbedded with the Lower Cretaceous source rocks, and the Upper Cretaceous sandstone beds. Basement rocks (granite, granodiorite, and gneisses) and oil shale deposits may represent potential unconventional hydrocarbon exploration. The current integrated study provides an insight that should guide future hydrocarbon exploration campaigns in the basin.

Babouri-Figuil basin

3D viscoelastic models of slip-deficit rate along the Cascadia subduction zone

Interseismic deformation in the Pacific Northwest is constrained by the horizontal crustal velocity field derived from the Global Positioning System (GPS) in addition to vertical rates derived from GPS, leveling, and tide gauge measurements. Such measurements were folded in to deformation models of fault slip rates as part of the 2023 National Seismic Hazard Model (NSHM) update. Here I build upon one of the contributing models, the viscoelastic earthquake-cycle model of Pollitz [2022]. This model permits inclusion of effects of time-dependent viscoelastic relaxation within earthquake cycles (i.e., ‘ghost transients’) and laterally variable elastic and/or ductile material properties. I lever-age these capabilities to incorporate the Cascadia megathrust into Western U.S.-wide deformation models in which crustal fault slip rates are estimated simultaneously with slip deficit rates along the interplate boundary between the descending Juan de Fuca plate and North American plate. This effort includes construction of a margin-wide model of viscoelastic structure founded on the Slab 2.0 model and probes different models of the ductile properties of the surrounding oceanic asthenosphere, continental lower crust, and mantle asthenosphere. This results in new estimates of the distribution of slip deficit rate along the ∼ 1000 km long margin, highlights the importance of correcting for glacial-isostatic adjustment effects, and permits assessment of sensitivity of results to assumed ductile properties.

California, Oregon, Washington

Exploring the uncertainty of machine learning models and geostatistical mapping of rare earth element potential in Indiana coals, USA

Rare earth elements and yttrium (REEs) have a wide range of applications in high- and low-carbon technologies. The strategic significance of REEs has grown due to their expanding applications in manufacturing industries and the constrained availability of these essential resources. This research explores the applicability of machine learning models and their uncertainty for assessing the REE potential in coal beds using various coal parameters as inputs. The work focuses on developing a predictive model based on geological variables, excluding considerations related to potential shifts in the commodities market. The Indiana Coal Quality Database was used as the data source. The promising and unpromising indicators derived from the outlook coefficient of samples from the database were used as the REE potential indicator for machine learning classification models. The filter-based approach with bootstrap was used to evaluate the importance of the coal parameters and their prediction uncertainties. Four machine learning methods (linear discriminant analysis (LDA), random forest (RF), support vector machine (SVM), and artificial neural networks (ANN), a data balancing and augmentation approach (Synthetic Minority Over-sampling Technique), and bootstrap resampling techniques were used for building the models and evaluating their prediction capabilities under uncertainty. It was determined that the SVM bootstrap model with ten-times balanced and augmented data provided superior results compared with other models. Finally, stochastic spatial maps of the REE potential within the coal basin were generated using sequential indicator simulation. The spatial maps of the REE potential showed that a 29% area of the Indiana section of the Illinois coal basin has economic potential of REEs, with 90% confidence.

Indiana

Estimation of the accessible and useful resource base for electric-grade enhanced geothermal systems (EGS) resources of the Great Basin, USA

Scientists with the U.S. Geological Survey (USGS) recently completed a provisional assessment of the electric-grade geothermal resources associated with the low-permeability geologic formations of the Great Basin, USA, where resources are assumed to be accessible using enhanced geothermal systems (EGS) technologies (i.e., the engineering of sufficient permeability to facilitate efficient heat extraction). This assessment required estimation of the accessible resource base (electric-grade heat [>90ºC] at depths where drilling and stimulation are deemed achievable using current technology) and useful resource (heat that can be extracted from the accessible region). Electric-grade heat can be estimated from existing temperature models. The accessible resource base can be estimated as the electric-grade heat that exists at depths shallower than 6 km based on the limitations of current drilling and stimulation technologies, along with evidence for sustained natural fracture conductivity at depth. The useful part of the accessible heat can be estimated as the product of three efficiencies and factors: the heat extraction efficiency, the viable geology factor, and the reservoir spacing efficiency. The accessible and useful parts of the resource can be estimated in units of heat, or in units of electric power using an electrical conversion efficiency, which is a function of resource temperature. We also estimate the ranges for each of the efficiencies and describe the motivations behind the choice of best estimates used for the recent assessment. An analytic solution is provided for the useful resource above any depth (in units of electric power), where efficiency estimation assumes nearly steady heat extraction rates that cool reservoirs to 90ºC over 30 years of power generation.

Great Basin

Effect of deicing chemicals on the hydrologic environment in Massachusetts; evaluation of surface-water data collection

The objective of this study is to develop predictive relationships that can be used to describe the impact of highway deicing salts on the hydrologic environment. A method, presented in an earlier report, for estimating yearly mean chloride concentrations from estimated or actual runoff and salt-application data was applied to the period 1971 through1974 and the results compared with yearly mean chloride concentrations computed from records of specific conductance. Estimates were found to be from minus 63 percent to plus 56 percent in error when the error was expressed as the difference between estimated and computed values as a percentage of the computed values. Chloride concentrations used in this study are computed from specific conductance/chloride concentration relationships and records of specific conductance. Preliminary results from graphic analyses and least-squares regression analysis show that relationships between measured values of specific conductance and chloride concentration have correlation coefficients ranging from 0.30 to 0.97. Analysis of streamflow for all major dissolved constituents is recommended with the purpose of attempting to describe the variations in the specific conductance/chloride concentration relationships.

Massachusetts

Hydrogeologic mapping and three-dimensional geologic modeling of glacial deposits in a multicounty area of southeastern Michigan, northeastern Indiana, and northwestern Ohio

The glacial deposits underlying southeastern Michigan, northeastern Indiana, and northwestern Ohio are a substantial source of water to communities, agriculture, and industry in the region. Previous efforts to characterize aquifer materials in the area cited a need for additional information about the underlying hydrogeologic characteristics and related groundwater availability as well as improved mapping of the extent and properties of the glacial deposits. Recent U.S. Geological Survey multi-State compilations of water-well drilling records have greatly increased access to high-resolution geologic data, particularly in glacial depositional environments. This study by the U.S. Geological Survey, in cooperation with the Ohio Environmental Protection Agency, uses processed data from the State-managed collections of well records to characterize the glacial deposits in the study area using two methods. The first method creates two-dimensional maps of basic hydrogeologic information commonly required for assessments of groundwater availability, including (1) total thickness of glacial deposits, (2) total thickness of coarse-grained deposits, (3) specific-capacity-based transmissivity and hydraulic conductivity, and (4) texture-based estimated equivalent horizontal and vertical hydraulic conductivity and transmissivity. The second method builds a hydrogeologic framework of the complex glacial aquifer through construction of a volumetric geologic model by using three-dimensional kriging. Results of the volumetric model indicate that aquifer materials are primarily concentrated in the western parts of the study area near the Indiana-Ohio border. Coarse-grained sediments are also present as surficial deposits in the north of the study area where intermixing glacial advances created complex distributions of unconsolidated deposits. Two-dimensional maps of hydrogeologic properties support the volumetric model, showing thicknesses of coarse-grained deposits that reach up to 250 feet in the western sections of the study area and progressively thin to near absence in the east. Visualization of the aquifer materials with a volumetric model generally shows a highly discontinuous distribution of coarse- and fine-grained materials, with no clearly defined boundaries to delineate the extent of the aquifer. Comparisons of cross sections derived from the volumetric model with existing published maps support previous near-surface hydrogeologic interpretations while filling gaps where data are sparse, particularly in deeper parts of the aquifer. Both the two-dimensional maps and the volumetric model provide data that can directly inform assessments of groundwater availability, in addition to having future applications to studies of groundwater flow and transport.

Indiana, Michigan, Ohio

Seismic response of rock towers at the Trona Pinnacles (U.S.A.) to the 2019 Ridgecrest earthquake sequence: Theory, observations, and models

We analyze the seismic response of a class of fragile geologic features (FGFs), referred to as rock towers (RTs) at the Trona Pinnacles, a group of RTs in southern California that suffered strong shaking during the 2019 Ridgecrest earthquake sequence. FGFs, including RTs, may provide maximum constraints on past earthquake shaking intensity, and thereby support probabilistic seismic hazard assessments (PSHAs). In a rare case study, we explore the hypothesis that RT structural integrity is time dependent, as damage accumulates progressively. We develop finite‐element method (FEM) models of the RTs using photogrammetric shape models. We validate the models by comparing numerical simulations of their response to broadband ground shaking with low‐intensity seismic recordings obtained at the Pinnacles. Results of our simulations are in good agreement with the seismic recordings of actual earthquake aftershocks. We next use the results of the FEM models to analyze the response and evolution of RTs. Our analyses elucidate the influence of geometry over their seismic response, providing a rationale that may explain the rarity of slender RTs at Trona: high‐aspect‐ratio structures that respond in bending develop detrimental tensile stresses that crack the rock, whereas low‐aspect‐ratio ones’ response also includes shearing, which does not compromise material integrity as much as tension. Field measurements with a rebound hammer support this finding, suggesting that the material around the base of slender rocks has been weakened relative to other parts of the RT. We also study how to define simplified mechanical models (“archetypes”) to predict the natural frequencies of RTs. Results from our work illuminate the fundamental mechanisms of seismic response and progressive failure of RTs, and open new avenues of research to potentially incorporate these geologic features as long‐return period constraints on PSHA, in ways analogous to those of the widely used precariously balanced rocks.

California

A history of cryohydrogeology modeling and recent advancements through the integration of solute transport

Groundwater flow systems and permafrost are interrelated because permafrost thaw enhances permeability, while groundwater flow can advect heat and accelerate permafrost thaw (McKenzie et al. 2021). Given amplified climate change in cold regions, there is renewed interest in ‘cryohydrogeology’, the study of groundwater in cold regions. Many data-driven studies have shown that permafrost thaw is leading to activated aquifers and increased baseflow across the pan-Arctic region (e.g. Walvoord and Striegl 2007, Evans et al. 2020). Empirical evidence of a subsurface ‘replumbing’ (Walvoord and Kurylyk 2016) in permafrost regions raises questions about the fate of sequestered contaminants in the North (Langer et al. 2023). We will discuss the history of and emerging opportunities in cryohydrogeological modeling, with a focus on recent contaminant transport modeling.

Conference Paper

The potential impact of three-dimensional distributed slip models derived from real-time GNSS data on the performance of the ShakeAlert earthquake early warning system for slab interface earthquakes

The ShakeAlert® earthquake early warning (EEW) system is designed to warn users of imminent strong ground motion with sufficient time to take protective actions. ShakeAlert currently uses three algorithms to characterize the earthquake source. One estimates the location and magnitude using the first few seconds of the P wave and, while fast, tends to underestimate magnitude for M w 7.0+ earthquakes. A second estimates the location, orientation, length, and corresponding magnitude of a line source using observed peak ground acceleration and contributes primarily to M w 5.5+ earthquakes. The third infers earthquake magnitude from peak ground displacement measured using Global Navigation Satellite System (GNSS) data and offers nonsaturating magnitudes for M w 7.0+ earthquakes. Other EEW algorithms exist that infer temporally evolving spatially variable slip on a 3D fault surface from real‐time GNSS data, information that might enable more accurate and timely alerts in the event of large‐magnitude subduction interface earthquakes. Here, we evaluate the potential contribution of one such algorithm, BEFORES ( Minson et al. , 2014 ), to improve ShakeAlert performance through a simulated real‐time implementation of Bayesian evidence‐based fault orientation and real‐time earthquake slip (BEFORES) and other ShakeAlert algorithms using data for eight M w 7.6+ earthquakes. The test results demonstrate that BEFORES can produce well‐constrained and accurate magnitude estimates as soon as or sooner than other EEW algorithms, in turn enabling it to increase the amount of warning time users receive in many cases. However, with a modified Mercalli intensity (MMI) threshold of 3.5, which is commonly used for issuing alerts, BEFORES would tend to alert large geographic regions that did not feel strong shaking (MMI 6+). This effect can be mitigated using a higher alert threshold of MMI 4.5 without negative impact on the amount of warning time obtainable with BEFORES.

Bulletin of the Seismological Society of America

Surface-wave relocation and characterization of the 2024–2025 dyking episode along the Fentale–Dofen segment of the Ethiopian rift

From late-December 2024 to mid-March 2025, a 50-km-long dyke intrusion triggered over 300 earthquakes (magnitude 4 to 5.9) between Fentale and Dofen volcanoes along the Northern Main Ethiopian Rift. Dyke intrusions periodically occur along the Fentale–Dofen magmatic segment and are an expression of ongoing rift extension. Preliminary analyses using interferometric synthetic aperture radar revealed extensive ground deformation (up to 60 cm), which closely matched the temporal and spatial evolution of surface manifestations and earthquake locations from global catalogues. While global catalogues are critical for real-time monitoring, the precision of locations in remote and or sparsely instrumented regions can be low. In this investigation, we present surface-wave relocation results of the dyking episode that began near Fentale volcano in December 2024. We estimate relative locations using differential traveltimes measured from regional-to-teleseismic distance surface-wave observations of earthquakes reported by the U.S. Geological Survey. Relative relocations reduce the initial region of diffuse seismicity to a 50-km-long narrow band bounding the strike of surface manifestations and the zone of maximum surface deformation. We demonstrate the precision of surface-wave relocations over incremental time periods, capturing the progression of dyking from seismic onset through seismic migration and caldera subsidence. Results showcase the utility of surface-wave relocations in the characterization of dyking episodes and provide complementary insights into the current understanding of the Fentale–Dofen volcanic plumbing system.

Geophysical Journal International

Monazite and xenotime U-Pb geochronology and thermometry of the Blue Ridge and Inner Piedmont of North Carolina: Implications for the thermal-metamorphic evolution of the southern Appalachian metamorphic “core”

The southern Appalachian orogen preserves a complex distribution of metamorphism and deformation varying in timing, magnitude, and spatial extent. These complexities give rise to disparate interpretations for southern Appalachian tectonic evolution, which complicates the testing and interpretation of tectonic models in this system. New monazite (Mnz) and xenotime (Xtm) laser ablation split stream (LASS) analyses alongside Mnz-Xtm thermometry in the orogenic core in the eastern Blue Ridge (EBR), western Inner Piedmont (WIP), and Cat Square terranes (CST) of North Carolina yield new constraints that define distinct pro- and retrograde metamorphic events. The EBR preserves two prograde thermal events: the Taconic (∼470-440 Ma, >660°C) and Neoacadian (∼380-340 Ma, 600–700°C), separated by a period of cooling (exhumation?) and followed by garnet breakdown from 339 to 329 Ma. Evidence of pervasive Neoacadian ductile deformation in the EBR is largely limited to the Brevard fault zone (BFZ), indicating that a major rheological gradient existed across the BFZ during the Neoacadian and early Alleghanian. Southeast of the BFZ, in the WIP and CST, monazite data define a protracted Neoacadian evolution from early mineral growth at ∼405 Ma at ∼450–600°C to >700°C at ∼360 Ma, followed by early Alleghanian retrograde metamorphism and deformation (<345 Ma, 350–500°C). These constraints, together with previously reported thermobarometric data, define a P-T-t evolution for the WIP and CST consistent with Neoacadian crustal flow, while the coeval presence of a thermal-rheological boundary along the BFZ further supports a model of Neoacadian crustal “escape” flow within the orogen.

North Carolina

Origin of the Pd/Pt ratio of the J-M Reef, Stillwater Complex, Montana, USA

The J-M reef of the Stillwater Complex is characterized by a high Pd/Pt ratio (mean ~3.8 with a standard error of 0.03) with a homogeneous geospatial distribution at the deposit scale. In this contribution, we demonstrate that the Pd/Pt ratio of the reef is the product of equilibration of an immiscible sulfide liquid with a silicate melt rich in Pd relative to Pt. Despite the high tenors of the J-M reef sulfides (avg 2,700 ppm Pt and 770 ppm Pt), numerical modeling shows that the parental melts did not have extraordinary Pd and Pt concentrations. Instead, the initial composition of a plausible parental silicate melt can have Pd and Pt contents well within the expected range of a normal, mantle-derived partial melt (i.e., ~10–20 ppb for both Pd and Pt with Pd/Pt of ~1). The relative differences in the partitioning behavior of Pt and Pd between sulfide liquid and silicate melt are unlikely to produce a consistent Pd/Pt ratio across a wide range of silicate melt to sulfide liquid mass ratios (i.e., R factors). Instead, the pre-emplacement fractionation of Pt alloy from S-undersaturated silicate magma accounts for the homogeneous and high Pd/Pt ratio of the J-M reef. We show that batch equilibration of sulfide liquid with silicate melt can produce the high Pd/Pt ratios of the reef if the partition coefficients between sulfide liquid and silicate melt for Pd and Pt are extremely high (>10 6 ). In an alternative model, Pd enrichment could be achieved by sulfide upgrading in resident footwall mush even if the partition coefficients between sulfide liquid and silicate melt are relatively small (between 10 4 and 10 6 ) because the instantaneous mass ratio of silicate melt to sulfide liquid is small (R ≈ 100–700), so the partitioning behavior of Pt and Pd has little impact on the composition of sulfide liquid.

Montana

State of science, gap analysis, and prioritization for southeastern United States water-quality impacts from coastal storms—Fiscal year 2023 program report to the Water Resources Mission Area from the Water Availability Impacts of Extreme Events Program—Hurricanes

Tropical cyclones (coastal storm events that include tropical depressions, tropical storms, and hurricanes) cause landscape-scale disturbances that can lead to impaired water quality and thus reduce water availability for use. Stakeholders and scientists at local and national scales have illustrated a need for understanding these risks to water quality. A regional and comprehensive understanding of the impacts of tropical storms and hurricanes on surface-water and groundwater quality—and thus water availability—is lacking for potentially impacted coastal and inland areas. As the U.S. Geological Survey considers development of tools to predict the extent to which water-quality impacts of hurricanes affect water availability, an assessment of the state of the science of hurricane impacts is needed, including a gap analysis and prioritization of data and science needs. This assessment focuses on the southeastern coastal States.

southeastern United States

Analysis of minimum 7-day discharges and estimation of minimum 7-day, 2-year discharges for streamflow-gaging stations in the Brazos River basin, Texas

The minimum 7-day, 2-year discharge (7Q2) was estimated for each of 97 streamflow-gaging stations in the Brazos River Basin, Texas, with at least 10 years of unregulated or 10 years of regulated discharge data using the Pearson Type III distribution fitted by L-moment statistics of the annual minimum 7-day discharges. Forty-six stations have at least 10 years of unregulated discharge data only, 26 stations have at least 10 years of regulated discharge data only, and 25 stations have at least 10 years of unregulated and 10 years of regulated discharge data. Statistical tests were applied to the annual minimum 7-day discharges to determine if the unregulated and regulated data have trends; and to determine if the data before regulation are significantly different from the data after regulation. Results of the Mann-Kendall test indicate a significant trend at the 95-percent confidence level for 8 of the 46 unregulated-only stations, for 6 of the 26 regulated-only stations, and for 5 stations before regulation and for 10 stations after regulation of the 25 stations with both unregulated and regulated data. Results of the Wilcoxon rank-sum test indicate a significant difference between the unregulated and regulated annual minimum 7-day discharges at the 95-percent confidence level for 16 of the 25 stations with both types of data. The 7Q2 is zero for 31 of the 46 unregulated-only stations and for 9 of the 26 regulated-only stations. The 7Q2 also is zero for 9 stations before regulation and for 3 stations after regulation of the 25 stations with both unregulated and regulated data. The 7Q2 ranges from 33 to 631 cubic feet per second for the eight regulated-only stations on the mainstem Brazos River. For the seven stations on the mainstem Brazos River with at least 10 years of unregulated and 10 years of regulated discharge data, the 7Q2 ranges from 0 to 568 cubic feet per second before regulation and from 0.30 to 670 cubic feet per second after regulation. The 7Q2 increased during regulation for 19 of the 25 stations with both unregulated and regulated data. The effect of regulation for most stations generally was an increase in the magnitude of the 7Q2.

Texas

Are the horizontal-to-vertical spectral ratios of earthquakes and microtremors the same?

We consider the similarities and differences between earthquake and microtremor horizontal‐to‐vertical spectral ratios (eHVSR and mHVSR, respectively) using a dataset of 161 sites in southern California. Quantitative comparisons are made in terms of the eHVSR and mHVSR lognormal median curves, as well as the frequencies and amplitudes associated with the fundamental‐ and higher‐mode resonances where present. The results show only 58% of the eHVSR–mHVSR pairs agree in terms of their median curve and only 25% of the eHVSR–mHVSR pairs agree in terms of shared resonances, which increases to 68% if flat HVSRs are considered equivalent. Furthermore, while the shared resonances match very well in terms of frequency (root mean square error, RMSE, <0.11 Hz), the amplitudes of those resonances do not agree (RMSE >1.6). These findings demonstrate that while eHVSR and mHVSR agree at some sites, they are not equivalent at all sites. To investigate if the agreement between eHVSR and mHVSR could be related to features of the microtremor data, earthquake recordings, and/or the site conditions, three machine learning (ML) models at varying levels of interpretability are presented. The ML models—which include multivariate logistic regression, gradient‐boosted trees, and support vector machines—show only partial success at using site‐specific data to predict whether eHVSR and mHVSR will likely agree in terms of their median curve (accuracy of 78%) and number of resonances (accuracy of 84%). Therefore, we conclude that while eHVSR and mHVSR can be quite similar in terms of resonant frequencies at some sites, they are not identical at all sites. Furthermore, preliminary evidence shows that the agreement of eHVSR and mHVSR can be predicted a priori given features of the microtremor measurements, earthquake recordings, and site conditions, although a larger dataset will be necessary for developing a robust predictive model.

California