Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of the American Water Resources Association”

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

Use of precariously balanced rocks to constrain postglacial earthquake magnitudes in New England, United States

We use the age and fragility of 13 precariously balanced rocks (PBRs) in Massachusetts, New Hampshire, and southwestern Maine to evaluate the estimated probabilistic seismic hazard (PSH) from the U.S. Geological Survey’s National Seismic Hazard Model (NSHM), and also develop a map of maximum postglacial earthquake magnitudes in the region. The PBRs are erratic boulders left behind by the retreat of the Laurentide ice sheet between 13 and 18 ka, and published ice retreat maps are used to estimate their fragility ages. The fragilities and ages of the rocks appear to be consistent with the NSHM, in that the estimated earthquake magnitudes and resulting ground motions that would topple them are only expected on timescales longer than their fragility ages. This result, on the North American passive margin, contrasts with PBR studies at active plate boundaries, where the PBR constraints imply large reductions to estimated PSH. It may be that the PBRs are not old enough in the eastern United States to reveal the same PSH modeling issues as those revealed in the plate boundary settings. The study provides constraints on maximum magnitudes of postglacial earthquakes across much of Massachusetts, New Hampshire, and southwestern Maine, with the constraints limited to the ages of the PBRs.

Maine. Massachusetts, New Hampshire↗
USGS2024

Neotectonic mapping of Puerto Rico

Puerto Rico is part of the Puerto Rico-Virgin Islands microplate, along the Caribbean–North American plate boundary between the Puerto Rico trench subduction zone and the Muertos Trough incipient subduction zone. Despite recent seismicity and geodetically constrained deformation of ~3 mm/yr of left-lateral shear across the island, Quaternary fault locations remain largely uncertain. Preservation of recent faulting in the landscape is masked by distributed faulting, high weathering rates of the tropical climate, steep topography, frequent landsliding, and extensive agriculture and urbanization along coastlines and river valleys. We present remote neotectonic mapping of historical aerial imagery from the 1930s-60s and <1-m lidar topography, integrated with field observations, to create an updated active fault map of Puerto Rico. We focus on faults that offset younger geomorphic surfaces, with ages estimated by geologic mapping, OSL and radiocarbon dating, and morphologic interpretations. We present new evidence for Quaternary activity on seven faults, including the South Lajas, Salinas, Punta Montalva, Great Southern Puerto Rico, Cerro Goden, Parguera, and San Marcos faults. We find that active faulting occurs preferentially along the southern and western coasts and does not always spatially coincide with preexisting bedrock faults, possibly reflecting partial reactivation of older faults in the modern strain field.

Puerto Rico↗
USGS2024

Using surficial geologic maps, vegetation, and monitoring to address erosion impacts from grazing in Channel Islands National Park, California

Employing a map-unit classification scheme based on geomorphic process and age, the U.S. Geological Survey (USGS) mapped Quaternary surficial deposits of the five islands comprising Channel Islands National Park (CHIS), California, as no such maps previously existed. Mapping was motivated through an agreement with the National Park Service (NPS) to aid natural resource assessments, including post-grazing disturbance recovery and identification of mass wasting and tectonic hazards. The resulting detailed (1:12,000 scale) maps portray areas of upland erosional transport processes and alluvial, fluvial, eolian, beach, marine terrace, mass wasting, and mixed depositional processes at a scale commensurate with the process signal on the landscape. Detailed, GPS-located observations of sedimentology, pedogenic development, landscape position, and geomorphology constrain map unit identification. In addition, map-unit boundary delineation was determined in a GIS using a variety of high-resolution (sub-meter) aerial imagery and lidar-based DEMs complemented by derivative raster products. The GIS geodatabase schema used is GeMS compliant. Resolvable map-unit widths are as small as 10-15 meter. Locational accuracy of unit boundaries is estimated to be as low as 10 m or better. Detailed mapping is at a scale allowing determination of geomorphic transport processes operating on the landscape and conveys geologic information to land managers at a resolution commensurate with park resources. Map unit ages are constrained by relative soil chronostratigraphy, vertebrate and invertebrate fossils, and materials collected for dating control using radiogenic, luminescence, and cosmogenic radionuclide techniques. Associated topical studies examined connections between geology, surface water, and vegetation distribution and restoration.

California↗
USGS2026

Salinas Valley integrated hydrologic and reservoir operations models, Monterey and San Luis Obispo Counties, California

The area surrounding the Salinas Valley groundwater basin in Monterey and San Luis Obispo Counties of California is a highly productive agricultural area, contributes substantially to the local economy, and provides a substantial portion of vegetables and other agricultural commodities to the Nation. This region of California provides about half of the Nation’s lettuce, celery, broccoli, and spinach each year. Thus, this agricultural area provides substantial volumes of agricultural products not just for California but for the United States. Changes in population and increased agricultural development, which includes a shift toward more water-intensive crops, and climate variability, have put increasing demand on both surface-water and groundwater resources in the valley. This situation has resulted in water management challenges in the Salinas Valley that generally relate to the distribution of the water supply throughout the basin. Where and when the water is present in the surface and subsurface does not coincide with where and when the water is needed. Historically, to deal with the distribution issue, water has been used conjunctively in the valley. Conjunctive use is a water management strategy that coordinates surface-water and groundwater use to maximize water availability. Groundwater is used throughout the Salinas Valley to meet water demands when surface-water supplies are insufficient. The availability of surface water is constrained by climate. Precipitation and streamflow vary seasonally and year to year. Although there are two reservoirs in the Salinas Valley to capture and store water during wet periods, the only conveyance of reservoir water to coastal agricultural areas is the Salinas River. Increasing demand for groundwater and surface-water resources throughout the Salinas Valley has resulted in undesirable effects from unsustainable water use, such as surface-water depletion, groundwater-level declines, storage depletion in the principal aquifers, and seawater intrusion. To address these escalating issues, local communities, water management agencies, and groundwater sustainability agencies are evaluating how to sustainably manage both their surface-water and groundwater resources. To meet water demands and reduce the undesirable effects of unsustainable water use, continued conjunctive management of surface water and groundwater would ideally incorporate strategies to deal with increases in demand and climate variability. To evaluate the challenging water management issues in the Salinas Valley, the U.S. Geological Survey, Monterey County Water Resources Agency, and the Salinas Valley Basin Groundwater Sustainability Agency developed a comprehensive suite of models that represent the Salinas Valley hydrogeologic system called the Salinas Valley System Model. The geologic framework is known as the Salinas Valley Geologic Framework and was developed to characterize the subsurface using various topographic and geologic data sources, including information on hydrogeologic units, their surfaces and extents, geologic structures, lithology, and elevations from borehole data and cross sections, as well as details on faults and existing models. The surface-water model is called the Salinas Valley Watershed Model and simulates the Salinas River watershed. Monthly surface-water inflows into the integrated hydrologic model domain were simulated using the Salinas Valley Watershed Model. The historical model uses historical climate data, water and land use data, and reservoir releases to simulate agricultural operations, including landscape water demands, diversions, and reclaimed wastewater. The operational model adds an embedded reservoir operations framework to the simulation of the historical model that allows specified operational rules to simulate reservoir releases and changes in reservoir storage. The operational model assumes current reservoir operations and constant land use, which differs from historical conditions. Thus, the operational model is a hypothetical baseline model that can be used by local water managers to evaluate and quantify potential benefits of water supply projects. Together, the geologic framework, watershed, historical, and operational models form a tool that can be used to simulate irrigated agriculture and associated reservoir operations of the integrated hydrologic system of the Salinas Valley.

California↗
USGS2024

New methodology for assessing underground natural gas storage resources – Example from Michigan Basin, United States

Energy consumption in the United States (U.S.) and across the world is shifting away from traditional fossil fuels like coal and oil, and towards natural gas and renewable sources, including hydrogen. Because gas demand is typically greatest during cold seasons and renewable sources sometimes produce variable supplies, it is important to store energy for use when demand exceeds supply. Whereas batteries and tanks typically store energy above ground, geologic (underground) storage may be able to retain much greater quantities of energy over much longer durations (e.g., Matos et al. 2019; Buursink et al. 2023). Consequently, the U.S. Geological Survey (USGS) is developing a methodology to locate new geologic energy storage sites with an initial focus on depleted hydrocarbon reservoirs followed by estimates of associated pore space or gas storage capacity. Depleted hydrocarbon fields in Michigan Basin in the U.S. already host large quantities of injected natural gas (about 681 billion cubic feet of seasonal gas capacity in 42 gas storage facilities; U.S. Energy Information Administration 2024). Therefore, we used this basin as an example to assess potential new storage resources alongside existing storage capacity and carbon dioxide sequestration resources (e.g., Katz and Coats 1968; Haagsma et al. 2020). In Michigan Basin gas storage is taking place in reservoirs within several key hydrocarbon producing formations with different lithologies, including the informal Stray sandstone of the Michigan Formation and the pinnacle reefs of the Niagara Group (Gautier et al. 1995; Swezey et al. 2015). The assessment methodology introduced here consists of two steps relying on both well information and reservoir data.

Michigan↗
USGS2024

Status of water-level altitudes and long-term and short-term water-level changes in the Chicot and Evangeline (undifferentiated) and Jasper aquifers, greater Houston area, Texas, 2023

Since the early 1900s, groundwater withdrawn from the primary aquifers that compose the Gulf Coast aquifer system—the Chicot, Evangeline, and Jasper aquifers—has been an important source of water in the greater Houston area, Texas. This report, prepared by the U.S. Geological Survey in cooperation with the Harris-Galveston Subsidence District, City of Houston, Fort Bend Subsidence District, Lone Star Groundwater Conservation District, and Brazoria County Groundwater Conservation District, is one in an annual series of reports depicting the status of water-level altitudes and water-level changes in these aquifers in the greater Houston area. In this report, the Chicot and Evangeline aquifers are treated as a single aquifer for the purposes of providing annual assessments of regional-scale water-level altitudes and water-level changes over time. In 2023, shaded depictions of water-level altitudes for the Chicot and Evangeline aquifers (undifferentiated) ranged from about 286 feet (ft) below the North American Vertical Datum of 1988 (NAVD 88) to about 169 ft above NAVD 88. The largest decline in water-level altitudes indicated by the 1977–2023 long-term water-level-change map was in south-central Montgomery County southeast of The Woodlands. In comparison, the 1990–2023 long-term water-level-change map depicts the largest declines in water-level altitudes in localized areas at or near certain wells in parts of northwestern Harris County and south-central Montgomery County. The largest rise in water-level altitudes for 1977–2023 is depicted in a relatively large area in southeastern Harris County, whereas the largest rise in water-level altitudes for 1990–2023 is depicted in a relatively large area in central Harris County. The 5-year short-term water-level-change map depicts the largest declines at three wells in northern Fort Bend County, one well in western Harris County, and three wells in south-central Montgomery County and the largest rise at one well in central Harris County. The 1-year short-term water-level-change map depicts the largest declines at one well in northern Fort Bend County and two wells in southwestern Harris County and the largest rises at one well in northern Brazoria County and one well in south-central Montgomery County. In 2023, shaded depictions of water-level altitudes for the Jasper aquifer ranged from about 242 ft below NAVD 88 to about 218 ft above NAVD 88. The 2000–23 long-term water-level-change map depicts water-level declines throughout the study area where water-level-measurement data from the aquifer were collected, with the largest declines in north-central Harris County and south-central Montgomery County south of The Woodlands. The 5-year short-term water-level-change map depicts the largest declines at two wells in central Montgomery County near Conroe and two wells in south-central Montgomery County southeast of The Woodlands and the largest rise at one well in western Montgomery County. The 1-year short-term water-level-change map depicts the largest declines at four wells in south-central Montgomery County southeast of The Woodlands and one well in central Montgomery County near Conroe and the largest rises at two wells in western Montgomery County.

Texas↗
USGS2024

Geologic map and structure sections along the southern part of the Bartlett Springs Fault Zone and adjacent areas from Cache Creek to Lake Berryessa, northern Coast Ranges, California

Introduction Located in the Coast Ranges of northern California, the Bartlett Springs Fault Zone is the easternmost fault in the San Andreas Fault system in northern California. The fault is a right-lateral, strike-slip structure considered capable of producing an earthquake of moment magnitude 7. The purpose of this mapping is to better characterize the geology and earthquake hazards associated with the southern part of the Bartlett Springs Fault Zone and to help identify any evidence of active uplift on the faults bounding the Coast Ranges. Although the area immediately surrounding the Bartlett Springs Fault Zone is sparsely populated, its southern segment presents a potential seismic hazard to northern California communities as far away as the San Francisco Bay region and Sacramento. There are also nearby water resources, mineral resources, and public lands used for public recreation. The Coast Ranges of northern California are a series of northwest-southeast-oriented mountain ranges and valleys located north of the San Francisco Bay region, between the Pacific Ocean to the west and the Sacramento Valley to the east. The region has rugged terrain, high mountain peaks that reach more than 2,400 meters above sea level, isolated and narrow valley bottoms on which most human settlements are located, and large drainage systems that tend to follow the northwest-southeast-oriented topographic grain. The physiographic character of the region is shaped by its bedrock geology, deformational history, and active faulting. The basement rocks of the northern Coast Ranges consist of the Franciscan Complex and the Great Valley complex, the latter of which consists of two informal units, the Coast Range ophiolite and the Great Valley sequence. The Franciscan Complex and the Great Valley complex are in structural contact along the Coast Range Fault, a regional-scale structure and fundamental crustal boundary. The Franciscan Complex and the Great Valley complex are superposed by active, northwest-southeast-striking strike-slip faults that are associated with seismicity swarms. These active strike-slip faults can produce moderate to large earthquakes that have moment magnitudes of 7–8. In places, these active structures bound large ranges and valleys, suggesting that much of the modern topographic expression is the result of active deformation processes. This report contains new 1:24,000-scale geologic mapping along the southern part of the Bartlett Springs Fault Zone between Clear Lake and Lake Berryessa. The map area spans 738 square kilometers in northern Napa County, southern Lake County, and parts of Yolo and Colusa Counties. The south and east borders of the map are 90 kilometers north of San Francisco and 70 kilometers west of Sacramento, respectively. The map area is within the Knoxville mining district, which has a history of mercury and gold mining dating back to the mid-19th century. The two main towns in the region, Lower Lake and Clearlake, California, are west-northwest of the map area. Approximately 71,000 people live in the cities and rural communities located within a 40-kilometer radius of the center of the map area. The bedrock geology, cross sections, and structural data presented herein are critical for evaluating the long-term evolution of the Bartlett Springs Fault Zone. This work will supplement studies on local seismic hazards, liquefaction potential, landslide hazards, earthquake geology, natural resources, groundwater resources, engineering geology, and tectonic history by providing the background information for site-specific investigations on these subjects.

California↗
USGS2026

Cotton farming pesticides affect ileal microbiota activity expressions of virulome but not resistome or metabolic pathways in a sedentary wild passerine

The increased use of agrochemicals to enhance crop production has had detrimental environmental effects including implication in the sharp decline of North American farmland-breeding birds. Here, using a combination of deep shotgun metatranscriptomics and pesticide exposure data, we sought to assess whether exposure to cotton ( Gossypium spp.) production had a differential effect on ileum multi-kingdom microbial activity, metabolism, anti-microbial resistance, and virulence factors of sedentary northern mockingbirds ( Mimus polyglottos ) sampled from two cotton-producing areas (16 birds in total) and one uncultivated area (7 birds) in Texas, USA. Both Shannon Index values (Adj. r 2 = 0.174, F (1,21) = 5.633, p = 0.027) and a Mantel test (Spearman ρ = 0.184, p = 0.013) supported a relationship between metabolically active microbiota Bray–Curtis dissimilarities and differences in pesticide mixtures among study areas. Virulence factor richness (Adj. r 2 = 0.182, F (1,21) = 5.890, p = 0.024), Shannon Index (Adj. r 2 = 0.231, F (1,21) = 7.612, p = 0.012), and load (sum of virulence factor abundances; Adj. r 2 = 0.160, F (1,21) = 5.194, p = 0.033) were related to total pesticide load (total quantity of pesticides). We found no pesticide effects on expression of either antimicrobial resistance genes or metabolic pathways.

Texas↗
USGS1990

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

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

Tennessee↗
USGS2025

Status of water-level altitudes and long-term and short-term water-level changes in the Chicot and Evangeline (undifferentiated) and Jasper aquifers, greater Houston area, Texas, 2024

Since the early 1900s, groundwater withdrawn from the primary aquifers that compose the Gulf Coast aquifer system—the Chicot, Evangeline, and Jasper aquifers—has been an important source of water in the greater Houston area, Texas. This report, prepared by the U.S. Geological Survey in cooperation with the Harris-Galveston Subsidence District, City of Houston, Fort Bend Subsidence District, Lone Star Groundwater Conservation District, and Brazoria County Groundwater Conservation District, is one in an annual series of reports depicting the status of water-level altitudes and water-level changes in these aquifers in the greater Houston area. In this report, the Chicot and Evangeline aquifers are treated as a single aquifer for the purposes of providing annual assessments of regional-scale water-level altitudes and water-level changes over time. In 2024, shaded depictions of estimated water-level altitudes for the Chicot and Evangeline aquifers (undifferentiated) ranged from about 301 feet (ft) below the North American Vertical Datum of 1988 (NAVD 88) to about 184 ft above NAVD 88. The largest decline in water-level altitudes depicted by the 1977–2024 long-term water-level-change map was in south-central Montgomery County. In comparison, the 1990–2024 long-term water-level-change map depicts the largest declines in water-level altitudes in an area northwest of The Woodlands and in an area of northern Waller County. The largest rise in water-level altitudes for 1977–2024 is depicted in an area of east-central Harris County, whereas the largest rise in water-level altitudes for 1990–2024 is depicted in an area of central Harris County. The 5-year short-term water-level-change map depicts the largest declines in several parts of the study area, but these declines are concentrated primarily in northern Fort Bend County, southwestern Harris County, and south-central Montgomery County. The largest rise for 2019–24 is depicted at a well in northern Fort Bend County. The 1-year short-term water-level-change map depicts the largest declines at a well in northern Fort Bend County and a well in west-central Harris County. The largest rise for 2023–24 is depicted at a well in east-central Fort Bend County. In 2024, shaded depictions of estimated water-level altitudes for the Jasper aquifer ranged from about 255 ft below NAVD 88 to about 321 ft above NAVD 88. The 2000–24 long-term water-level-change map depicts the largest water-level decline in an area of central San Jacinto County; the largest rise is depicted in an area of central Grimes County. The 5-year short-term water-level-change map depicts the largest declines across parts of central and southern Montgomery County and at one well in north-central Harris County. The largest rise for 2019–24 is depicted at a well centered on the Montgomery-Grimes County line. The 1-year short-term water-level-change map depicts the largest declines at two wells in south-central Montgomery County and one well in northwestern Montgomery County on the west side of Lake Conroe. The largest rises during 2023–24 are depicted at one well in northwestern Montgomery County and one well in south-central Montgomery County.

Texas↗
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

Fingerprinting magmatic REE deposit sources with zircon petrochronology

Carbonatites and associated alkaline silicate rocks are of considerable economic interest due to their enrichments in rare earth elements. The petrogenesis and source(s) of these complexes, however, are poorly understood. Models propose either mantle plume-derived carbon-rich melts or a mantle source enriched by subduction-related metasomatism. We use zircon trace elements to fingerprint carbonatite-alkaline complex sources, focusing on the economically significant 1.4 Ga Mountain Pass intrusive suite (MPIS). Autocrystic zircon from MPIS alkaline silicate rocks are enriched in Th and U relative to Nb suggesting a subduction influenced, oxidized source region. Lower Sc/Yb, higher Ti concentrations, and the absence of Eu anomalies in autocrystic zircon suggest derivation from a less hydrous, hotter, and deeper mantle source relative to arc-related 1.8–1.6 Ga inherited zircon. These data are best explained by syn- to post-collisional Mesoproterozoic reactivation of a subduction-metasomatized mantle source. The source and tectonic setting of the MPIS contrasts with Mesoproterozoic Gifford Creek and Bayan Obo carbonatites which exhibit plume-like compositions and are associated with rifts, suggesting different sources and tectonic settings for these economically significant deposits. We demonstrate that zircon petrochronology is a robust method for distinguishing carbonatite sources and can inform more targeted exploration strategies for critical mineral resources.

Conference Paper↗
USGS2024

Earthquake relocations delineate discrete a fault network and deformation corridor throughout Southeast Alaska and Southwest Yukon

Deformation in southeastern Alaska and southwest Yukon is governed by the subduction and translation of the Pacific-Yakutat plates relative to the North American plate in the St. Elias region. Despite notable historical seismicity and major regional faults, studies of the region between the Fairweather and Denali faults are complicated by glacial coverage and the remote setting. In the last decade, significant improvements have been made to the density of regional broadband seismometer networks. We relocate more than 5,000 earthquakes between 2010 and 2021 in the region of southeastern Alaska and southwestern Yukon utilizing these improved seismic networks. With reductions in catalog uncertainty, particularly in depth, we quantify the thickness of the seismogenic layer in the crust throughout the region and locate seismicity on a shallow network of upper-crustal faults. Relocated earthquakes, combined with an updated focal-mechanism catalog, permit estimating and classifying motion of active faults. This includes mapping the Totschunda-Fairweather “Connector” fault, which plays an important role in explaining regional deformation, and identifying new faults like the Kathleen Lake fault. We draw similarities between our seismic observations and simplified conceptual models of regional tectonics, which describe a dominant transpressional regime and localized slip partitioning. Our results support a hypothesis where current deformation is taking place on a well-defined and evolved network of shallow faults in the corridor between the Totschunda-Fairweather “Connector” and Denali faults.

Alaska, Yukon↗
USGS2026

Groundwater budget for the Mountain Home area, southern Idaho, 2022–23

The U.S. Geological Survey, with funding from the Idaho Department of Water Resources, developed a groundwater budget for the Mountain Home area in southern Idaho for irrigation year 2023 (November 1, 2022–October 31, 2023). This study focused on the water balance across the Cinder Cone Butte Critical Groundwater Area (CGWA), Mountain Home Groundwater Management Area (GWMA), and the rest of the study area (RoSA), compiling data from various sources, including precipitation records, groundwater level measurements, metered groundwater pumpage data, surface water diversions and evapotranspiration (ET) estimates derived from remote sensing satellite imagery, and ground-based reference data. Key inflow components included recharge from applied surface water irrigation (which incorporates incidental recharge from irrigation practices and conveyance losses), estimated tributary streamflow, and estimated mountain block recharge. The key outflow components were groundwater pumpage for irrigation, municipal, industrial, and domestic uses, and ET. Recharge from applied irrigation and mountain block recharge were the largest inflows, and groundwater pumpage for irrigation was the largest outflow. The CGWA had a positive groundwater budget residual of 2,170 acre-feet (acre-ft), which contrasts with observed long-term groundwater level declines and historical trends of storage depletion. This positive residual is likely associated with unquantified outflows, including lateral groundwater flow out of the subregion, or other complexities, such as overestimated tributary contributions relative to the actual recharge for the 2023 water budget. The GWMA exhibited a positive residual of 56,563 acre-ft, primarily owing to recharge from applied surface water irrigation and areal recharge during a wetter-than-average year, which allowed irrigation entities to deliver more water from in-basin and out-of-basin reservoirs. The RoSA showed a large positive residual of 124,933 acre-ft. The interpretation of these positive residuals must account for significant uncertainties, including estimations of areal recharge, tributary streamflow (particularly losses and diversions), ET, the volume of surface water loss to the Snake River, lateral groundwater flows between subregions and across study area boundaries, and the unquantified groundwater discharge to the Snake River. These uncertainties, in combination with the complex hydrogeologic controls on water movement and limitations of remotely sensed data, directly affect the accuracy of water availability assessments. Future data collection efforts would help reduce these uncertainties and support water resource management decisions in the Mountain Home area. Key efforts could include installing additional streamflow gaging stations (particularly to quantify tributary losses and gains and surface water losses to the Snake River), improving groundwater pumpage metering, and validating remotely sensed ET data with ground-based measurements. Furthermore, to better quantify unrepresented or highly uncertain fluxes, focused investigations on groundwater discharge to the Snake River, lateral groundwater flows between subregions and across study area boundaries, and a more robust determination of the actual influence and volume of mountain block recharge would help refine future water availability assessments for the Mountain Home area.

Idaho↗
USGS1998

Assessment of the hydraulic connection between ground water and the Peace River, west-central Florida

The hydraulic connection between the Peace River and the underlying aquifers along the length of the Peace River from Bartow to Arcadia was assessed to evaluate flow exchanges between these hydrologic systems. Methods included an evaluation of hydrologic and geologic records and seismic-reflection profiles, seepage investigations, and thermal infrared imagery interpretation. Along the upper Peace River, a progressive long-term decline in streamflow has occurred since 1931 due to a lowering of the potentiometric surface of the Upper Floridan aquifer by as much as 60 feet because of intensive ground-water withdrawals for phosphate mining and agriculture. Another effect from lowering the potentiometric surface has been the cessation of flow at several springs located near and within the Peace River channel, including Kissengen Spring, that once averaged a flow of about 19 million gallons a day. The lowering of ground-water head resulted in flow reversals at locations where streamflow enters sinkholes along the streambed and floodplain. Hydrogeologic conditions along the Peace River vary from Bartow to Arcadia. Three distinctive hydrogeologic areas along the Peace River were delineated: (1) the upper Peace River near Bartow, where ground-water recharge occurs; (2) the middle Peace River near Bowling Green, where reversals of hydraulic gradients occur; and (3) the lower Peace River near Arcadia, where ground-water discharge occurs. Seismic-reflection data were used to identify geologic features that could serve as potential conduits for surface-water and ground-water exchange. Depending on the hydrologic regime, this exchange could be recharge of surface water into the aquifer system or discharge of ground water into the stream channel. Geologic features that would provide pathways for water movement were identified in the seismic record; they varied from buried irregular surfaces to large-scale subsidence flexures and vertical fractures or enlarged solution conduits. Generally, the upper Peace River is characterized by a shallow, buried irregular top of rock, numerous observed sinkholes, and subsidence depressions. The downward head gradient provides potential for the Peace River to lose water to the ground-water system. Along the middle Peace River area, head gradients alternate between downward and upward, creating both recharging and discharging ground-water conditions. Seismic records show that buried, laterally continuous reflectors in the lower Peace River pinch out in the middle Peace River streambed. Small springs have been observed along the streambed where these units pinch out. This area corresponds to the region where highest ground-water seepage volumes were measured during this study. Further south, along the lower Peace River, upward head gradients provide conditions for ground-water discharge into the Peace River. Generally, confinement between the surficial aquifer and the confined ground-water systems in this area is better than to the north. However, localized avenues for surface-water and ground-water interactions may exist along discontinuities observed in seismic reflectors associated with large-scale flexures or subsidence features. Ground-water seepage gains or losses along the Peace River were quantified by making three seepage runs during periods of: (1) low base flow, (2) high base flow, and (3) high flow. Low and high base-flow seepage runs were performed along a 74-mile length of the Peace River, between Bartow and Nocatee. Maximum losses of 17.3 cubic feet per second (11.2 million gallons per day) were measured along a 3.2-mile reach of the upper Peace River. The high-flow seepage run was conducted to quantify losses in the Peace River channel and floodplain between Bartow and Fort Meade. Seepage losses calculated during high-flow along a 7.2-mile reach of the Peace River, from the Clear Springs Mine bridge to the Mobil Mine bridge, were approximately 10 percent of the river flow, or 118 cubic feet per second. Calculated seepages along the Peace River in Hardee and De Soto Counties were inconclusive, because most seepages were within the range of discharge measurement error. Two continuous aerial thermal infrared imagery surveys were conducted to locate sites of ground-water discharge along the Peace River. Although temperature and hydrologic conditions were ideal to observe spring flow using thermal infrared imaging techniques, no sources of ground-water discharge were identified using this method. Diffuse ground-water seepage may, however, provide significant ground-water discharge.

Florida↗
USGS2025

Assessing legacy nitrogen in groundwater using numerical models of the Long Island aquifer system, New York

Nitrogen transported along groundwater flow paths in coastal aquifers can contribute substantially to nitrogen loading into surface water receptors, particularly in hydrologic systems dominated by groundwater discharge. Nitrogen entrained in the aquifer is a function of land use and associated nitrogen sources at the time of groundwater recharge, which may differ considerably from present-day sources. Legacy nitrogen can result in substantial discrepancies between observed present-day nitrogen loading to surface water receptors and loading estimated from present-day sources. Additionally, legacy nitrogen can continue to discharge into surface waters after nitrogen mitigation actions have been undertaken. Here, we use a numerical modeling framework to compare three methods of estimating time-varying historical nitrogen loads to four water bodies (receptors) on eastern Long Island, New York. The methods span a range of data requirements and process complexity, from instantaneous receptor loads calculated from steady-state groundwater contributing areas, to transient loads estimated by explicitly simulating legacy groundwater nitrogen transport over a century with large changes in nitrogen sources and hydrologic conditions. The effects of legacy nitrogen on estimated receptor loads varied temporally and spatially within the study area. Depending on antecedent nitrogen inputs and hydrologic conditions, historical annual nitrogen loads estimated from transient simulations accounting for legacy nitrogen can be quite similar (<10% difference) or substantially different (±100%) from those estimated from simpler instantaneous methods. Continued input of present-day nitrogen sources using methods that account for legacy nitrogen results in asymptotic increases in receptor nitrogen loads over time, indicating that simulated present-day receptor nitrogen loads are not in equilibrium with present-day inputs. For these receptors in disequilibrium, models simulating transient groundwater nitrogen transport could be used to account for legacy nitrogen lag times to help resource managers evaluate the potential effectiveness of proposed nitrogen mitigation actions.

EarthArXiv↗
USGS2025

Combining scanning electron microscopy, X-ray diffraction, and X-ray fluorescence to characterize shear zones at the Pogo gold deposit, Alaska

This study employs a multi-method analytical approach to characterize the mineralogical, geochemical, and textural properties of fault rocks from the Pogo gold mine in the Yukon-Tanana Upland, central Alaska. Specifically, we examine cataclasites, to document the structural and geochemical evolution of shear zones and their associations with gold mineralization. To investigate the shear zone, we integrate portable X-ray fluorescence (pXRF), scanning electron microscopy-based automated mineralogy (SEM-AM), X-ray diffraction (XRD), and high-resolution micro-X-ray fluorescence (micro-XRF) mapping. These methods collectively provide insights into bulk and trace element chemistry, mineralogical composition, and deformation-related textures across multiple scales. Handheld pXRF enables rapid geochemical screening, guiding SEM-AM and XRD analyses to ensure consistent mineralogical interpretation. X-ray diffraction identifies and quantifies crystalline phases, while SEM-AM produces high-resolution mineral maps, revealing mineral abundances, grain-scale textures, and gold associations. Micro-XRF mapping further refines our understanding by showing visual trace element distributions at sub-millimetre resolution. By integrating these techniques, we improve our understanding of the nature and geochemistry of Pogo shear zones, their role in gold mineralization, and support metallurgical processing strategies. This approach enhances exploration models and resource characterization for structurally complex gold deposits.

Alaska↗
USGS2024

Spatial distribution of API gravity and gas/oil ratios for petroleum accumulations in Upper Cretaceous strata of the San Miguel, Olmos, and Escondido Formations of the south Texas Maverick Basin—Implications for petroleum migration and charge history

The Maverick Basin of south Texas is currently undergoing active exploration and production of gas and oil from tight sandstone reservoirs. The most productive tight sandstones in the basin are in the Upper Cretaceous San Miguel, Olmos, and Escondido Formations. These units are second only to the Eagle Ford Shale in terms of cumulative production volumes. The structural history of the Maverick Basin, from rifting to subsidence to exhumation, has had a profound effect on the characteristics of these reservoirs and the petroleum resources contained therein. This U.S. Geological Survey review of the production history of these strata reflects a recent shift from conventional production to horizontal drilling (unconventional) that exploits low permeability reservoirs in previously overlooked areas of existing oil and gas fields in southern Texas, typically outside of established field boundaries. To investigate the physical properties of the Maverick Basin hydrocarbon accumulations, this case study compiled American Petroleum Institute (API) gravity measurements and calculated cumulative gas/oil ratios (GOR) for thousands of producing wells from the San Miguel, Olmos, and Escondido Formations. Maps were generated from the compiled well production data to show the spatial heterogeneity of API gravity and GOR values for the three formations within the Maverick Basin and immediately outside the basin to the northeast. Within the Maverick Basin, the spatial patterns of API gravity values indicate lighter oils downdip towards the southern basin edge. GOR values indicative of wet and dry gases within the basin are seen interspersed, with values that correspond to black and heavy oils. Differences in the spatial patterns of the petroleum properties within the Maverick Basin are interpreted as effects of Eocene basin inversion caused by Laramide orogenic deformation, and the resulting reservoir exhumation of basin strata. East of the Maverick Basin, spatial distributions of API gravity and GOR values show progressively heavier oils updip to the northwest, grading to dry gases downdip to the southeast, which correlates to the oil and gas windows of the underlying Eagle Ford Shale. Correlation of API gravity and GOR values from the San Miguel, Olmos, and Escondido Formations with thermal maturity data from the Eagle Ford Shale suggests that the Eagle Ford Shale is the petroleum source, and that petroleum migration was approximately vertical for areas to the east of the Maverick Basin. The discontinuity of API gravity and GOR properties within the Maverick Basin implies a complex petroleum charge history, possibly involving the remigration of petroleum and the addition of petroleum from other source intervals in Mexico, to the southwest. Depressurization of exhumed, overpressured reservoirs of the San Miguel, Olmos, and Escondido Formations can explain the intermittent occurrence of gas production throughout the southern Maverick Basin by exsolution of gas from formation brines and the resulting dry gas flushing of hydrocarbon-charged reservoirs. The introduction of dry gas through flushing can, in turn, explain why the patterns of API gravity and GOR values are so dissimilar in the Maverick Basin. This process has implications for possible future production of unconventional resources from undiscovered tight-gas reservoirs in strata of the San Miguel, Olmos, and Escondido Formations, and a different approach to petroleum exploration may be needed in the Maverick Basin relative to exploration techniques applied in other basins within the northern Gulf of Mexico.

Texas↗