Geology Reports⌕ Search

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 18Linked to original sources
USGS1993

Hydraulic properties of the Madison aquifer system in the western Rapid City area, South Dakota

Available information on hydrogeology, data from borehole geophysical logs, and aquifer tests were used to determine the hydraulic properties of the Madison aquifer. From aquifer-test analysis, transmissivity and storage coefficient were determined for the Minnelusa and Madison aquifers, and vertical hydraulic conductivity (K v ') along with specific storage (Ss') for the Minnelusa confining bed. Borehole geophysical well logs were used to determine the thickness and location of the Minnelusa aquifer, the lower Minnelusa confining bed, and the Madison aquifer within the Madison Limestone. Porosity values determined from quantitative analysis of borehole geophysical well logs were used in analyzing the aquifer-test data. The average porosity at the two aquifer-test sites is about 10 percent in the Minnelusa aquifer, 5 percent in the lower Minnelusa confining bed, and 35 percent in the Madison aquifer. The first aquifer test, which was conducted at Rapid City production well #6, produced measured drawdown in the Minnelusa and Madison aquifers. Neuman and Witherspoon's method of determining the hydraulic properties of leaky two-aquifer systems was used to evaluate the aquifer-test data by assuming the fracture and solution-opening network is equivalent to a porous media. Analysis of the aquifer test for the Minnelusa aquifer yielded a transmissivity value of 12,000 feet squared per day and a storage coefficient of 3 x 10 -3 . The specific storage of the Minnelusa confining bed was 2 x 10 -7 per foot, and its vertical hydraulic conductivity was 0.3 foot per day. The transmissivity of the Madison aquifer at this site was 17,000 feet squared per day, and the storage coefficient was 2 x 10 -3 . The second aquifer test, which was conducted at Rapid City production well #5 (RC-5) produced measured drawdown only in the Madison aquifer. Hantush and Jacob's method of determining the hydraulic properties of leaky confined aquifers with no storage in the confining bed was used to evaluate the aquifer-test data by assuming the fracture and solution-opening network is equivalent to a porous media. The analysis of data from the RC-5 aquifer test showed that transmissivity was not equal in all directions. Hantush's method was used to determine the direction of radial anisotropy and magnitude of the major and minor axes of transmissivity. The major axis of transmissivity is at an angle of 42° east of north, and the transmissivity along this axis is about 56,000 feet squared per day. The minor axis of transmissivity is at an angle of 48° west of north, and the transmissivity along this axis is about 1,300 feet squared per day. The major axis of transmissivity intersects Cleghorn Springs, a large resurgent spring on the west edge of Rapid City. The shape of the potentiometric contours of the Madison aquifer near RC-5 agree with the orientation of the transmissivity ellipse. The average value of the storage coefficient from the isotropic analysis of the aquifer-test data was 3.5 x 10 -4 , and the average vertical hydraulic conductivity of the lower Minnelusa confining bed was 9.6 x 10 -3 foot per day.

South Dakota↗
USGS2025

A comprehensive geologic framework of the National Crustal Model for seismic hazard studies in the conterminous United States

A three-dimensional (3D) geologic framework has been developed for the conterminous United States (U.S.) as part of the U.S. Geological Survey National Crustal Model to enhance seismic hazard modeling. The geologic framework is created from geologic maps and multiple subsurface geologic unit boundaries including the base of the Miocene, Cenozoic, Phanerozoic, and the Mohorovičić discontinuity. Modifications are made to surficial geologic maps to remove discontinuities across state and country borders. The subsurface distribution of rock type and age is extrapolated from the surface, seeded with subsurface geologic information, and constrained by a map of basement geology. The framework provides the basis for estimates of subsurface seismic velocity and density that is needed to improve estimates of earthquake ground shaking and seismic hazard. The present framework greatly expands and updates a previously published 3D geologic framework of the western part of the U.S. that was itself a first-of-its-kind digital 3D portrayal of the nation.

conterminous United States↗
USGS2025

Geologic map of the Buckner 7.5-Minute quadrangle, Louisa County, Virginia

The Buckner 7.5-minute quadrangle straddles three terrane boundaries in the Piedmont Physiographic Province in central Virginia: the Chopawamsic terrane, the Elk Hill Complex, and the Goochland terrane. In much of the map area, the Elk Hill Complex separates the Chopawamsic and Goochland terranes. Rocks of the Chopawamsic terrane include Ordovician metavolcanic and metasedimentary rocks of the Chopawamsic Formation, Ordovician to Silurian granodiorite sheet intrusions, and Paleozoic mafic intrusions. Silurian to Devonian rocks of the Quantico Formation, mostly garnet-mica schist, crop out in the northwesternmost part of the map area, and are in unconformable contact with rocks of the Chopawamsic Formation on the southeastern limb of the Quantico synclinorium. The main map unit in the Elk Hill Complex is Neoproterozoic mica gneiss, which is in pre-metamorphic fault contact with rocks of the Chopawamsic Formation to the west. The main map unit of the Goochland terrane is the Maidens Gneiss. Except for Jurassic diabase dikes, all rocks on the Buckner 7.5-minute quadrangle were metamorphosed to amphibolite facies during the Alleghanian orogeny and preserve multiple compositional and phyllosilicate penetrative foliations. Evidence of amphibolite-facies metamorphism during the Taconic orogeny is preserved in rocks of the Elk Hill Complex. The entire width of the Maidens Gneiss on the Buckner 7.5-minute quadrangle is within the Spotsylvania high-strain zone and amphibolite-facies mylonitic textures are pervasive. Quartz veins and Jurassic diabase dikes crosscut all older rocks of the quadrangle. Multiple levels of terrace deposits are present along and near the major streams of the quadrangle. The lower terrace deposits are likely remnants of former positions of the Little River on the landscape, whereas higher deposits may be remnants of former deposits of the Atlantic Coastal Plain that covered this portion of the Piedmont Province. A linear cluster of aftershocks from the magnitude 5.8 earthquake that occurred near Mineral, Virginia, in 2011 defines the Fredericks Hall fault, which is at depth on the Buckner quadrangle. Most of the aftershocks occurred in the core of the Elk Creek antiform and have no relation to faults mapped at the surface. Several abandoned crushed stone and building stone quarries, as well as a mica prospect, exist in the quadrangle.

Virginia↗
USGS2026

Utilization of multiple geochronology techniques to constrain the age of laterization and mineralization of the world-class Mount Weld rare earth element deposit, Western Australia

Pervasive chemical weathering on stable cratons may form thick regoliths and elemental enrichment, but constraining the age of regolith formation is challenging. In this study we utilize multiple geochronological techniques on different minerals from the world-class Mount Weld rare earth element (REE) deposit, formed by lateritic weathering of a carbonatite, to constrain the age of formation and provide insight into landscape evolution. The oldest dates, ca. 100 to 38 Ma, are from Lu-Hf dating of churchite [HREE(PO 4 )·2(H 2 O)], a heavy REE phosphate mineral. Growth bands on individual minerals show a younging outwards. 40 Ar/ 39 Ar geochronology of cryptomelane [K(Mn 4+ ,Mn 2+ )₈O₁₆] yielded dates from ca. 40 to 27 Ma. Similarly, (U-Th)/He geochronology of goethite [FeO(OH)] yielded dates ranging from ca. 45 to 19 Ma. Integrating results into regional constraints, suggests 1) churchite formed by mineral saturation in a karst-like setting below the water table from ca. 100 to 40 Ma, 2) with minor uplift and erosion, cryptomelane and goethite formed at or near the water table between ca. 45 and 19 Ma, 3) after ca. 15 to 10 Ma chemical weathering within the profile had ended. Other studies document that the region experienced minimal uplift and a wet, warm climate from ca. 100 Ma to 15 Ma. These conditions and the high carbonate content of the carbonatite promote extensive chemical weathering, a deep weathering profile, and the preservation of the weathered section. This study highlights the use of multiple geochronological techniques utilizing different minerals to provide insight into how laterites form and to constrain the timing and history of the formation of this important mineral deposit.

Western Australia↗
USGS2024

The geometry of fault reactivation and uplift along the central part of the Maacama fault zone, northern California Coast Ranges (USA)

Fault reactivation of bedrock structures in active fault zones influences stress state and earthquake rupture phenomena through the introduction of weak slip surfaces that impact fault zone geometry and width. Yet, geometric relationships between modern faults and older reactivated faults are difficult to quantify in rocks that have experienced multiple deformation episodes. We used new geologic mapping, geomorphic tools, and structural modeling to quantify rock uplift and subsurface fault geometry of the central part of the Maacama Fault Zone near Ukiah, California, USA, and the surrounding area. Results suggest that the northern Mayacamas Mountains are in a tectonically driven disequilibrium, with differential rock uplift focused on the western side of the range. Steeply east-dipping fault surfaces and splays characterize the geometry of the Maacama Fault Zone. We mapped two newly identified faults to the east of the main Maacama Fault, the Cow Mountain–Mill Creek Fault, and Willow Creek Fault, which align with a moderately east-dipping cluster of microseismicity between 4–10 km depth beneath the Mayacamas Mountains. Static stress modeling on the Maacama Fault Zone and newly identified faults to the east quantify slip tendency values of 0.5–0.4, which suggests that the faults are moderately to poorly suited for slip in the modern stress field and may be weak. We infer that modern uplift is driven by oblique reverse, up-to-the-east, dip-slip motion on the reactivated Cenozoic Cow Mountain–Mill Creek and Willow Creek Faults as material is advected through a restraining bend on the Maacama Fault. This study shows that reactivated bedrock faults increase the fault zone width and introduce fault surfaces that contribute a component of vertical deformation and uplift in major strike-slip fault zones. Deformation is accommodated on an interconnected network of new and reactivated faults that delineate a complex seismic hazard.

California↗
USGS2025

Application of Hydrologic Simulation Program—FORTRAN (HSPF) as part of an integrated hydrologic model for the Salinas Valley, California

The U.S. Geological Survey (USGS), in cooperation with the Monterey County Water Resources Agency, conducted studies to help evaluate the surface-water and groundwater resources of the Salinas Valley study area, consisting of the entire Salinas River watershed and several smaller, adjacent coastal watersheds draining into Monterey Bay. The Salinas Valley study area is a highly productive agricultural region that depends on the coordinated use of surface water and groundwater to meet demand for irrigation and public water supply. To continue to meet these demands, a better understanding of the historical water balance and the effects of water-resource development on the long-term sustainability of water resources in the Salinas Valley study area is needed.

California↗
USGS2025

Revisiting an enigma on California's north coast: The Mw6.5 Fickle Hill earthquake of 21 December 1954

Many earthquakes occur along the North Coast of California in the vicinity of the Mendocino Triple Junction (MTJ), where the Pacific, Gorda, and North American (NA) plates meet, and on the adjacent plate boundaries. The MTJ marks the nexus of the Mendocino and San Andreas faults with the Cascadia subduction zone (CSZ). Historically, most large earthquakes around the MTJ have been within the offshore Gorda plate and its subducted portion beneath the NA plate. North of the MTJ, active faults mapped in the NA plate are part of the CSZ fold‐and‐thrust belt. Although some events have been detected in the NA plate, no large historic events have been associated with mapped surface faults. The 21 December 1954 M w 6.5 earthquake in Humboldt County is one possible exception. Using published data from catalogs and articles, unpublished data from Berkeley’s archives, and S ‐ P times interpreted from two U.S. Coast and Geodetic Survey (USCGS) accelerometers, we determine a probability cloud for the earthquake’s hypocenter using NonLinLoc. The highest probability location lies beneath Fickle Hill just east of the city of Arcata, California, at 40.87° N, 124.03° W, and ∼11 km depth. Using P ‐wave polarities from Berkeley stations and the digitized waveforms from the accelerometers, we find that the focal mechanism most consistent with the data indicates thrust movement with strike, dip, and rake of 350°, 10°, and 90°, respectively, at a depth of 14 km. Given the depth uncertainties of both this event and the megathrust, this implies that the earthquake most likely took place on the subduction interface rather than on the mapped faults in the Mad River fault zone that trend 322° and dip to the northeast. The revisited intensity in the epicentral region also supports a location beneath Fickle Hill to the east of the city of Arcata, California.

California, Idaho, Nevada, Oregon, Washington↗
USGS2026

Development of a hydrologic flow model of the Zuñi Salt Lake and surrounding areas, west-central New Mexico

The terminal Zuñi Salt Lake is in a maar in west-​central New Mexico and contains hypersaline water that has long been used by Native Americans for religious purposes and the collection of salt. The U.S. Geological Survey (USGS), in cooperation with the Zuni Tribe of the Zuni Reservation, New Mexico, and Bureau of Reclamation, developed the Zuñi Salt Lake hydrologic flow model to simulate the steady-state conditions of the lake and surrounding groundwater-flow system. The model is a tool that can be used to analyze the potential hydrologic responses to resource development in the region under different water-use scenarios, which can support negotiations that could help protect future use of the lake. The USGS modular hydrologic model, MODFLOW 6, was used to simulate the hydrologic system of the lake and surrounding areas. To define the layering of the model, a three-dimensional hydrogeologic framework was constructed, which included seven informal hydrogeologic units and the locations of faults that truncate and offset the units. Calibration of the model was performed using the iterative ensemble smoother included in a parameter estimation software suite (PEST++). The iterative ensemble smoother approach resulted in a posterior parameter ensemble of model realizations with parameter and simulated values that show similar success in matching observations. Groundwater flux to the Zuñi Salt Lake simulated by using the posterior parameter distribution ranged from 26 to 730 acre-feet per year with a mean, median, and standard deviation of 530, 570, and 170 acre-feet per year, respectively. The relative contribution of groundwater discharge to the lake was simulated using a particle-tracking model for MODFLOW (MODPATH) to be from, in order of decreasing contribution, the Mesaverde, Cenozoic, Mancos, and Dakota hydrogeologic units. Simulations of future conditions using the ensemble of models can provide an empirical estimate of forecast uncertainty without substantial additional analysis.

New Mexico↗
USGS2024

Slip rate for the Rose Canyon fault through San Diego, California, based on analysis of GPS data: Evidence for a potential Rose Canyon–San Miguel-Vallecitos fault connection?

The Rose Canyon fault is the southern extension of the larger Newport–Inglewood–Rose Canyon fault system, which represents a major structural boundary in the Inner Continental Borderland (ICB) offshore of southern California. Ten to fifteen percent of total plate boundary motion in southern California is thought to be accommodated by the faults of the ICB, but the exact distribution of slip is uncertain. With an onshore segment, the Rose Canyon fault offers an opportunity to measure the slip rate using traditional geodetic methods. In this study, we use Global Positioning System (GPS) surface velocities from a combined campaign and continuous GPS network to constrain elastic models of the Rose Canyon fault. We then compare the observed surface velocities with proposed conceptual models of regional fault connections that facilitate the transfer of slip into the Rose Canyon fault to assess how well the observations are explained by the models. The results of elastic half‐space models suggest that the Rose Canyon fault may be slipping toward the higher end of geologic estimates, with the preferred model indicating a slip rate of 2.4 ± 0.5 mm/yr. Although limited in terms of near‐fault benchmarks, we find an improved model fit using an asymmetrical elastic half‐space model and a higher slip rate, suggesting a potential rheological contrast across the Rose Canyon fault, similar to observations from the northern Newport–Inglewood fault segments. Observed GPS surface velocities, background seismicity, and gravity anomalies south of San Diego Bay point toward a more easterly trace for the Rose Canyon fault, suggesting a possible connection with the San Miguel–Vallecitos fault system. Such a connection could increase the potential rupture lengths of future earthquakes and have important consequences for regional seismic hazards.

California↗
USGS1996

Hydrogeologic investigation and simulation of ground-water flow in the Upper Floridan Aquifer of north-central Florida and southwestern Georgia and delineation of contributing areas for selected city of Tallahassee, Florida, water-supply wells

A 4-year investigation of the Upper Floridan aquifer and ground-water flow system in Leon County, Florida, and surrounding counties of north-central Florida and southwestern Georgia began in 1990. The purpose of the investigation was to describe the ground-water flow system and to delineate the contributing areas to selected City of Tallahassee, Florida, water-supply wells. The investigation was prompted by the detection of low levels of tetrachloroethylene in ground-water samples collected from several of the city's water-supply wells. Hydrologic data and previous studies indicate that; ground-water flow within the Upper Floridan aquifer can be considered steady-state; the Upper Floridan aquifer is a single water-bearing unit; recharge is from precipitation; and that discharge occurs as spring flow, leakage to rivers, leakage to the Gulf of Mexico, and pumpage. Measured transmissivities of the aquifer ranged from 1,300 ft 2 /d (feet squared per day) to 1,300,000 ft2/d. Steady-state ground-water flow in the Upper Floridan aquifer was simulated using a three-dimensional ground- water flow model. Transmissivities ranging from less than 5,000 ft2/d to greater than 11,000,000 ft 2 /d were required to calibrate to observed conditions. Recharge rates used in the model ranged from 18.0 inches per year in areas where the aquifer was unconfined to less than 2 inches per year in broad areas where the aquifer was confined. Contributing areas to five Tallahassee water-supply wells were simulated by particle- tracking techniques. Particles were seeded in model cells containing pumping wells then tracked backwards in time toward recharge areas. The contributing area for each well was simulated twice, once assuming a porosity of 25 percent and once assuming a porosity of 5 percent. A porosity of 25 percent is considered a reasonable average value for the Upper Floridan aquifer; the 5 percent porosity simulated the movement of ground-water through only solution-enhanced bedding plains and fractures. The contributing areas were generally elliptical in shape, reflecting the influence of the sloping potentiometric surface. The contributing areas delineated for a 5 percent porosity were always much larger than those determined using a 25 percent porosity. The lowest average ground-water velocity computed within a contributing area, using a 25 percent porosity, was 1.0 ft/d (foot per day) and the highest velocity was 1.6 ft/d. The lowest average ground-water velocity, determined using a 5 percent porosity, was 2.4 ft/d and the highest was 7.4 ft/d. The contributing areas for each of the five wells was also determined analytically and compared to the model-derived areas. The upgradient width of the simulated contributing areas were larger than the upgradient width of the analytically determined contributing areas for four of the five wells. The model could more accurately delineate contributing areas because of the ability to simulate wells as partially penetrating and by incorporating complex, three-dimensional aquifer characteristics, which the analytical method could not.

Florida, Georgia↗
USGS2024

A seismic nodal deployment to understand magmatic structure in the vicinity of the Pahala earthquake swarm

In summer-fall 2022, 80 three-component SmartSolo IGU-BD3C-5 nodal seismometers were deployed surrounding the Pāhala seismic swarm on the Island of Hawaiʻi, with the goal of improving seismicity catalogs, and seismic velocity images of the crust and upper mantle in this region. The Pāhala swarm, located south of Mauna Loa and Kīlauea, has been the site of a multiyear sustained swarm of seismicity at depths of ~ 25 – 40 km, with order of magnitude increases in rate in 2015, and then again in 2019. This seismicity is possibly related to the input of magma from the mantle plume below, which may then be subsequently transported to volcanic edifices. However, these processes remain enigmatic, in part due to a lack of precise earthquake locations and seismic velocity models in this region. Here we provide an overview of the deployment, an assessment of the quality of the collected data, and discuss the viability of the dataset for local earthquake relocation, tomography, and teleseismic receiver functions. Through comparisons with proximal permanent broadband and short period instruments, we find that the nodes produce high quality data, particularly at periods shorter than 5 s, although we find, document, and correct discrepancies with the gain and polarities of the instruments. We successfully record signals from teleseismic earthquakes, even at periods longer than 5 s (the corner of the flat response of the nodes). We also record local earthquakes, including details related to source characteristics. This indicates that the data is likely to prove useful for investigations using both local and teleseismic earthquake signals to better understand the connections between the deep and shallow magmatic systems of Hawaiʻi. While this deployment provides a snapshot in time, its success may provide a useful benchmark for future studies as the volcanic systems of Hawaiʻi continue to evolve in the future.

Hawaii↗
USGS2024

Surficial geology and Quaternary fault map of the Las Vegas Valley, Clark County, Nevada

The Surficial Geology and Quaternary Fault Map of the Las Vegas Valley, Clark County, Nevada is a 1:50,000-scale compilation of published 1:24,000-scale geologic maps integrated with new field and desktop mapping. This geologic map compilation and GIS database are part of a broader study on the Quaternary faults in the Las Vegas Valley, the results of which are summarized in this report. The map compilation utilized a variety of digital base maps to evaluate the characteristics of surficial deposits and Quaternary fault scarps including lidar data in undeveloped parts of the valley and orthorectified historical aerial photos in the urbanized center of the valley. The map distinguishes twenty surficial map units including alluvial-fan and ephemeral-wash deposits, groundwater discharge deposits of the Las Vegas Formation, and alluvium deposits of Las Vegas Wash. Quaternary faults include the Quaternary Las Vegas Valley fault system in the center of the valley and the Frenchman Mountain fault system along the eastern margin. All Quaternary fault traces in the valley were evaluated, and the previously published fault mapping was modified as needed. The new mapping and uniform correlation of surficial deposits displaced by Quaternary faults yields a consistent valley-wide characterization of the recency of fault activity. The map and GIS database include a Quaternary geochronology compilation of 119 previously published surficial deposit ages of various types, 35 new luminescence ages, and 3 new radiocarbon ages. The new age data further bracket the ages of the Las Vegas basin Quaternary stratigraphy and provide new constraints on the timing of Quaternary fault activity.

Nevada↗
USGS2025

The role of alkali bicarbonate-sulfate brines in the genesis of carbonatite REE resources at the Bear Lodge Alkaline Complex, Wyoming

Rare-earth element (REE) resources in the Bear Lodge Alkaline Complex, Wyoming, are hosted in a variably leached carbonatite dike swarm spatially related to bodies of diatreme breccia. This study examines fluid inclusions in carbonatite dikes, peripheral fluorite breccias, and smoky quartz veins to reconstruct the physiochemical conditions of REE mineralization. Results reveal a multi-stage evolution of carbonatitic fluids: 1) high-temperature (330–432°C), CO 2 -rich magmatic fluids, 2) REE-enriched alkali bicarbonate-sulfate brines, and 3) low-temperature meteoric (117–182°C) or diluted magmatic fluids. Multiphase inclusions (stage-2) contain burbankite, nahcolite, and alkali sulfate daughter crystals, linking the brines to early burbankite [(Na,Ca) 3 (Sr,Ba,REE) 3 (CO 3 ) 5 ] mineralization in carbonatite. REE mobility has been linked to the stability of aqueous complexes with ligands such as Cl-, SO 4 2- , and HCO 3 - . However, light REEs were not initially transported significant distances, rather they crystallized as burbankite from brines or brine-melts associated with the carbonatite magma. Crystallization of burbankite removed Na from the system, leaving residual fluids with high K/Na ratios and promoting potassic fenitization throughout the complex. Fractionation of light and heavy REEs and involvement of low temperature (<300°C), high K/Na bicarbonate brines may have facilitated transport and formation of peripheral occurrences that are more enriched in heavy REEs at Bear Lodge.

Wyoming↗
USGS2025

Hydrogeologic framework and conceptual model of the Red River alluvial aquifer east of Lake Texoma, southeastern Oklahoma, 1980–2022

The 1973 Oklahoma Groundwater Law (Oklahoma Statutes §82-1020.5) requires that the Oklahoma Water Resources Board conduct hydrologic investigations of the State’s groundwater basins to support a determination of the maximum annual yield for each groundwater basin. At present (2025), the Oklahoma Water Resources Board has not established a maximum annual yield for the Red River alluvial aquifer east of Lake Texoma. To support the evaluation and determination of a maximum annual yield, a hydrogeologic framework and conceptual groundwater-flow model were developed to assess groundwater availability in the Red River alluvial aquifer east of Lake Texoma. The scope of this hydrologic investigation is the alluvium and terrace containing the Red River alluvial aquifer in Oklahoma between Lake Texoma, the Texas State line, and the Arkansas State line, an extent referred to in this report as “the eastern part of the Red River alluvial aquifer.” Parts of the alluvium and terrace extent in Arkansas and Texas are included in some analyses to address hydrologic influences from outside the aquifer’s boundaries in Oklahoma. The eastern part of the Red River alluvial aquifer in southeastern Oklahoma consists of approximately 401,280 acres of Quaternary alluvium and terrace deposits associated with the Red River and its major tributaries. Mean annual recharge to the aquifer for the 1980–2022 study period was estimated to be 8.62 inches per year, or 17.98 percent of the mean annual precipitation over the same period (47.94 inches). This mean annual recharge rate is equivalent to an inflow of approximately 288,250 acre-feet per year for the eastern part of the Red River alluvial aquifer. Recharge estimated using the Soil-Water-Balance code accounts for 98.7 percent of the conceptual-model inflows to the eastern part of the Red River alluvial aquifer. Saturated-zone evapotranspiration accounts for 11.9 percent and net streambed seepage accounts for 87.4 percent of the outflows in the conceptual model.

Arkansas, Oklahoma, Texas↗
USGS2025

Methodology and technical input for the 2025 U.S. List of Critical Minerals—Assessing the potential effects of mineral commodity supply chain disruptions on the U.S. economy

The Secretary of the Interior, acting through the Director of the U.S. Geological Survey, is tasked by section 7002 (“Mineral Security”) of title VII (“Critical Minerals”) of the Energy Act of 2020 (Public Law 116–260, December 27, 2020, 116th Congress) with reviewing and revising the methodology used to evaluate mineral commodity supply risk and the U.S. List of Critical Minerals (LCM) no less than every 3 years. Following two previous LCM assessments, this analysis represents the latest technical input for evaluating each mineral commodity’s supply risk and determining their recommended status on the LCM. We evaluated mineral commodity supply risk using two criteria: (1) an economic effects assessment that quantified the potential effects of various trade disruption scenarios on the U.S. economy, and (2) an examination of whether the mineral commodity’s U.S. supply chain relied on a sole domestic producer that represented a single point of failure. For the first criterion, postdisruption equilibrium quantities and prices for each mineral commodity were calculated based on their price elasticities of supply and demand and the availability of excess production capacity for each yearlong foreign trade disruption scenario. Subsequently, a nonlinear optimization routine was used with detailed economic input-output tables to estimate the potential economic effects on the U.S. economy of over 1,200 scenarios for 84 mineral commodities. After accounting for the probability of each scenario’s occurrence, the overall results are presented in terms of changes in U.S. gross domestic product (GDP) by individual industry and the economy overall. The results, which ranged from a net decrease in U.S. GDP of nearly $4.5 billion to a net increase of $33 million, largely reflect U.S. import dependency and world production concentration. Using the Jenks natural breaks optimization method, a statistical classification technique, we categorized the mineral commodities into several classes based on this overall risk quantification. Mineral commodities with annualized probability-weighted net decreases in U.S. GDP greater than $2 million were recommended for inclusion on the LCM. If a mineral commodity did not meet the threshold for inclusion on the LCM under the first criterion, its domestic supply chain was examined under the second criterion, which recommended a mineral commodity for inclusion on the LCM if there was only a single domestic producer. Ultimately, the two criteria resulted in the recommendation of the addition of six mineral commodities (in descending risk order, potash, silicon, copper, silver, rhenium, and lead) to and the removal of two mineral commodities (arsenic and tellurium) from the LCM. By using an economic effects assessment, the results of this analysis provide a prioritization that can also be compared directly against other risk analyses and the cost of various risk mitigation strategies.

Open-File Report↗
USGS2024

Pesticide occurrence in shallow groundwater in three regions of agricultural land use: Baldwin County, the Wiregrass region, and the Tennessee River valley region of Alabama, 2009–20

As part of a cooperative investigation between the U.S. Geological Survey and the Alabama Department of Agriculture and Industries, a network of 22 groundwater wells were sampled from 2014 through 2020 for about 230 pesticide and pesticide degradate compounds. Wells were located in three regions of intensive agricultural land use in Alabama: Baldwin County, the Wiregrass region, and the Tennessee River valley region. Metolachlor sulfonic acid, a degradate of the herbicide metolachlor, was the most frequently detected compound, occurring in about 70 percent of the samples. Three other compounds, metolachlor, atrazine, and 2-chloro-4-isopropylamino-6-amino-s-triazine, were also detected in more than half of the samples. Metolachlor and its degradates accounted for 33 of the 50 greatest compound concentrations study-wide, including the maximum pesticide concentration across all compounds (62,500 nanograms per liter). The frequency and magnitude of detections of many specific pesticide compounds varied among the three regions, but all detected pesticide concentrations were well below the U.S. Environmental Protection Agency maximum contaminant levels and applicable human health benchmarks. Sample results were combined with results of previous (2009–13) sampling to provide a continuous time-series of data for 2009–20. More than half of the 289 pesticide compounds analyzed during 2009–20 were not detected in any samples. Only four compounds were detected at great enough frequency throughout the 10 sampling years to evaluate patterns of change through time. Metolachlor and its degradate, metolachlor sulfonic acid, were frequently detected in all regions. Atrazine and its degradate, 2-chloro-4-isopropylamino-6-amino-s-triazine, were also detected in wells from all regions, but the variability and magnitude of concentrations were greatest in the Tennessee River valley region. No apparent temporal pattern in concentrations was found.

Alabama↗
USGS2024

Summary of annual repeat magnetotelluric surveys of the Geysers geothermal field

As part of a multi-year project, funded by the California Energy Commission, annual magnetotelluric (MT) surveys have been collected at The Geysers geothermal field in northern California with the goal of measuring temporal changes within the steam field. The repeat surveys started in 2021 and repeated a 2017 survey (Peacock et al., 2020) with further extension to the southern part of the geothermal field. Temporal variations in the MT transfer functions are observed to be spatially coherent and compartmentalized. Mapping residual phase tensor ellipses demonstrates the direction of maximum change is often aligned with existing fracture orientations. Three dimensional inversion of the MT data, using the inversion results from previous years as the starting model, indicates that the steam reservoir has generally become more resistive over time (~10%), suggesting more steam in the field. A few pockets within the steam field have become more electrically conductive over time and are collocated with injection wells, suggesting either more fluid content in those zones, less steam, or more saline fluids.

California↗