SEARCH · USGS Science

Results for “The Planetary Science Journal”

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

SourceResearch connectionFrom yearThrough year

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

1,665 records · Page 40Linked to original sources

USGS2025

Aeromagnetic and magnetotelluric imaging of west-central Idaho and the Stibnite-Yellow Pine mining district: A regional to district perspective

Aeromagnetic and magnetotelluric (MT) data are used to better understand the geology and mineral resources near the Stibnite-Yellow Pine mining district in central Idaho. The reduced-to-pole (RTP) transformation of regional-scale aeromagnetic data shows that allochthonous island-arc rocks west of the Salmon River suture are significantly more magnetic than the Laurentian continental rocks east of the suture and that the granitoids of the Idaho batholith have moderate to low magnetization in both early, metaluminous, and late, peraluminous phases. Application of tilt derivative to aeromagnetic data highlights major crustal-scale structures. The 5-km upward continued magnetic data indicate island-arc rocks have deep magnetic sources. The 110-km-long MT profile images resistivity structure to depths around 30 km. At shallow depths, resistivity corresponds to mapped geologic units, with moderate resistivities underlying volcanic and roof-pendant metasedimentary rocks and moderate to high resistivities occurring beneath the Idaho batholith. Crustal-scale moderate resistivities beneath the suture image the results of tectonomagmatic processes that accompanied suturing and translating allochthonous terranes. Low resistivity values beneath and fringing the batholith are derived from metasedimentary rocks that may have served as a melt source and reductant during melt generation and provided metals during later ore formation. In the Stibnite-Yellow Pine mining district, a high-resolution aeromagnetic compilation is shown to correlate with mapped lithologies and mineral deposit-related structures. The RTP transform distinguishes magnetic and nonmagnetic granitoid phases of the Idaho batholith. The tilt derivative highlights metasedimentary rocks, some of which are favorable ore hosts. The Meadow Creek fault hosts the Stibnite and Hangar Flats deposits and is imaged as a magnetic low due to hydrothermal alteration. Reconstructions of magnetic anomaly offsets and orebodies indicate around 3 km of post-95 Ma dextral separation, with some or all of the offset inferred to postdate the main Au mineralization episode (61–66 Ma).

Idaho↗

USGS2024

Parameter ESTimation with the Gauss–Levenberg–Marquardt algorithm: An intuitive guide

In this paper, we review the derivation of the Gauss–Levenberg–Marquardt (GLM) algorithm and its extension to ensemble parameter estimation. We explore the use of graphical methods to provide insights into how the algorithm works in practice and discuss the implications of both algorithm tuning parameters and objective function construction in performance. Some insights include understanding the control of both parameter trajectory and step size for GLM as a function of tuning parameters. Furthermore, for the iterative Ensemble Smoother (iES), we discuss the importance of noise on observations and show how iES can cope with non-unique outcomes based on objective function construction. These insights are valuable for modelers using PEST, PEST++, or similar parameter estimation tools.

Groundwater↗

USGS2024

Insights on gas hydrate formation and growth within an interbedded sand reservoir from well logging at the Qiongdongnan Basin, South China Sea

Although variable well log resolution and its control on saturation estimation has been studied, it has not been directly applied to a specific location to explore the nature of gas hydrate within a sand reservoir. We applied in-situ measurements of resistivities, neutron porosity, and gamma ray at two sites in the Qiongdongnan Basin, South China Sea (QDN-W05–2021 and QDN-W08–2021) to investigate the reservoir parameters of a hydrate-bearing sand reservoir. Our results show that gas hydrate is distributed in 5 zones with a total thickness of 10.7 m and an average saturation of 69% at the QDN-W05–2021 site, while they are distributed in 2 zones with a total thickness of 4.3 m and an average saturation of 49% at the QDN-W08–2021 site. We found that variances in saturations estimated from lateral-extra deep button (RX), phase shift (P40H-P40L), and attenuation (A40H-A40L) resistivities within the laterally mapped continuous sand body were affected by the nature of gas hydrate occurrences. Results indicate gas hydrate forms and accumulates at the center of the sand layer and tends to be less or not present toward the top and base. Integrated with seismic data, the in-situ measurements provide insights in the evolution of a mushroom-shaped, hydrate-gas reservoir system. In the system, free gas is likely horizontally transported from the top-center of the gas chimney to the surrounding areas in the early stage dominated by a warm-gas environment, whereas hydrate forms in the opposite pathway starting from the surrounding areas in the following stage with temperature reducing. Our study suggests that high-resolution in-situ measurements not only are a tool to identify the physical properties, but also can be used to help explain the physical process of hydrate growth and accumulation.

Marine Geology↗

USGS2025

The Long Island Sound and Watershed Metadata map application

The Long Island Sound and its watershed encompass an area of about 17,000 square miles and include the Connecticut, Housatonic, and Thames Rivers, which all drain to the sound. Dozens of organizations from government agencies, nonprofits, and Tribal Nations have developed projects and monitoring programs to analyze and protect the water resources of the watershed and sound. The abundance of data and lack of an existing searchable index require a centralized metadata repository to allow users to find water resources data more efficiently. The U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency and the Long Island Sound Study, has created an interactive map application to visualize and search for metadata information across organizations working to monitor and protect the Long Island Sound.

Connecticut, Massachusetts, New Hampshire, New Yor↗

USGS2026

Constraining source and path effects of large magnitude earthquakes using ground motion simulations

The purpose of this study is to use ground‐motion simulations to investigate ways in which source and path effects for large‐magnitude earthquakes can be represented in nonergodic ground‐motion models (GMMs). To achieve this, we designed a ground‐motion study in the San Francisco Bay Area that includes earthquakes with a broad range of magnitudes distributed uniformly on a fault plane, and sites covering a large range of rupture distances and azimuths. After running a large suite of kinematic simulations (magnitude 4–7), we then develop a nonergodic GMM with the simulated data. We find that trends in the within‐site residuals are affected significantly by the earthquake radiation pattern, rupture directivity, and slip patterns. Next, we modify an existing rupture directivity model to fit and remove the observed radiation pattern and rupture directivity from the residuals. We also minimize the contributions of slip patterns by averaging the within‐site residuals among multiple source realizations. Finally, after removing the source effects from the within‐site residuals, we compare the path effects computed with different magnitude groups using two approaches. The first approach only considers the small events that have the same shortest path to a site as the large events, whereas the second approach considers all small events on the fault plane. The results indicate that it is difficult to satisfactorily approximate the path effects of large events with those of small events using either approach, at least in the case of simulations.

California↗

USGS2025

Quality of groundwater used for domestic supply in the Gilroy-Hollister basin and surrounding areas, California, 2022

More than 2 million Californians rely on groundwater from domestic wells for drinking-water supply. This report summarizes a 2022 California Groundwater Ambient Monitoring and Assessment Priority Basin Project (GAMA-PBP) water-quality survey of 33 domestic and small-system drinking-water supply wells in the Gilroy-Hollister Valley groundwater basin and the surrounding areas, where more than 20,000 residents are estimated to utilize privately owned domestic wells. The study area includes the Llagas subbasin in the north, the North San Benito subbasin in the south, and the surrounding uplands. The study was focused on groundwater resources used for domestic drinking-water supply, which are mostly drawn from shallower parts of aquifer systems rather than those of groundwater resources used for public drinking-water supply in the same area. This assessment characterized the quality of ambient groundwater in the aquifer before filtration or treatment, rather than the quality of drinking water delivered to the tap. To provide context, the measured concentrations of constituents in groundwater were compared to Federal and California State regulatory and non-regulatory benchmarks for drinking-water quality. A grid-based method was used to estimate the areal proportions of groundwater resources used for domestic drinking wells that have water-quality constituents present at high concentrations (above the benchmark), moderate concentrations (between one-half of the benchmark and the benchmark for inorganic constituents, or between one-tenth of the benchmark and the benchmark for organic constituents), and low concentrations (less than one-half or one-tenth the benchmark for inorganic and organic constituents, respectively). This method provides statistically representative results at the study-area scale and permits comparisons to other GAMA-PBP study areas. In the study area, inorganic constituents in groundwater were greater than regulatory benchmarks (U.S. Environmental Protection Agency [EPA] or State of California maximum contaminant levels [MCLs]) for public drinking-water quality in 24 percent of domestic groundwater resources. The inorganic constituents present at concentrations greater than MCLs for drinking water were nitrate (as nitrogen), barium, chromium, and selenium. Total dissolved solids (TDS) or manganese were present at concentrations greater than the secondary maximum contaminant levels (SMCLs) that the State of California uses as aesthetic-based benchmarks in 48 percent of domestic groundwater resources. No volatile organic compounds or pesticide constituents were present at concentrations greater than regulatory benchmarks. Total coliform bacteria and enterococci were detected in 4 percent of domestic groundwater resources. Per- and polyfluoroalkyl substances (PFAS) were detected in 19 percent of domestic groundwater resources, and 10 percent had concentrations greater than recently enacted (April 2024) EPA MCLs. Physical and chemical factors from natural and anthropogenic sources that could affect the groundwater quality were evaluated using results from statistical testing of associations between constituent concentrations and potential explanatory variables. In this study, relevant physical factors include well construction characteristics, groundwater age, site proximity to groundwater recharge or discharge zones, and potential sources of contamination. Relevant chemical factors include the initial chemistry of the recharge water, the mineralogy of the aquifer sediments, and the subsequent shifts in chemistry as biologic and geologic reactions alter groundwater in the subsurface. Nitrate concentrations were correlated to agricultural land use, distance from the boundary of the Gilroy-Hollister Valley groundwater basin, and the proportion of modern (post-1950s) water captured by the well. Denitrification under anoxic redox conditions can mitigate some nitrate derived from fertilizer application. Total dissolved solids primarily were derived from water-rock interactions with soils and aquifer materials in the study area, but there were high concentrations where agricultural practices contributed additional TDS. Mineralogy of aquifer sediments and rocks also affect barium, selenium, boron, and chromium concentrations in the Gilroy-Hollister Valley groundwater basin. PFAS were positively correlated with urban land use and the proportion of modern water captured by the well.

California↗

USGS2024

Conceptualization and simulation of groundwater flow and groundwater availability in the Boone and Roubidoux aquifers in northeastern Oklahoma, 1980–2017

Oklahoma Groundwater Law (Oklahoma Statute § 82-1020.5) requires that the Oklahoma Water Resources Board conduct hydrologic investigations to determine the maximum annual yield for the State’s groundwater basins. The Boone and Roubidoux aquifers (also known as the Springfield Plateau aquifer and Ozark aquifer, respectively) are bedrock aquifers that extend from northeastern Oklahoma into Kansas, Arkansas, and Missouri. At present (2024), the Oklahoma Water Resources Board has yet to legally issue orders for the final determination of maximum annual yields for the Boone and Roubidoux aquifers. To support determination of a maximum annual yield, the U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, developed a hydrogeologic framework, a conceptual groundwater-flow model, and a calibrated numerical groundwater-flow model for the Boone and Roubidoux aquifers. Three types of groundwater-availability scenarios were simulated by using the calibrated numerical model. These scenarios were used to (1) estimate equal-proportionate-share groundwater withdrawal rates (groundwater withdrawal applied equally over the aquifer), (2) quantify the potential effects of projected groundwater withdrawals on groundwater storage over a 50-year period, and (3) simulate the potential effects of a hypothetical 10-year drought. For the Boone aquifer, equal-proportionate-share groundwater withdrawal rates were 1.10, 0.98, and 0.96 acre-feet per acre per year for the 20-, 40-, and 50-year scenarios, respectively. For the Roubidoux aquifer, equal-proportionate-share groundwater withdrawal rates were 1.76, 1.34, and 1.25 acre-feet per acre per year for the 20-, 40-, and 50-year simulations, respectively. For the 50-year scenarios, stream seepage was minimally affected. Over the 10-year drought scenario, groundwater storage in the Boone and Roubidoux aquifers decreased by 660,451 acre-feet (6.7 percent) and 508,472 acre-feet (1.0 percent), respectively.

Oklahoma↗

USGS1996

Magnitude and frequency of low flows in the Suwannee River Water Management District, Florida

Low-flow frequency statistics for 20 gaging stations having at least 10 years of continuous record and 31 other stations having less than 10 years of continu ous record or a series of at least two low- flow measurements are presented for unregulated streams in the Suwannee River Water Management District in north-central Florida. Statistics for the 20 continuous-record stations included are the annual and monthly minimum consecutive-day average low- flow magnitudes for 1, 3, 7, 14, and 30 consecutive days for recurrence intervals of 2, 5, 10, 20, and, for some long-term stations, 50 years, based on records available through the 1994 climatic year.Only theannual statistics are given for the 31 other stations; these are for the 7- and 30-consecutive day periods only and for recurrence intervals of 2 and 10 years only. Annual low-flow frequency statistics range from zero for many small streams to 5,500 cubic feet per second for the annual 30- consecutive-day average flow with a recurrenceinterval of 2 years for the Suwannee River near Wilcox (station 02323500). Monthly low-flow frequency statistics range from zero for many small streams to 13,800 cubic feet per second for the minimum 30-consecutive-day average flow with a 2-year recurrence interval for the month of March for the same station. Generally, low-flow characteristics of streams in the Suwannee River Water Management District are controlled by climatic, topographic, and geologic fac tors. The carbonate Floridan aquifer system underlies, or is at the surface of, the entire District. The terrane's karstic nature results in many sinkholes and springs. In some places, springs may contribute greatly to low streamflow and the contributing areas of such springs may include areas outside the presumed surface drainage area of the springs. In other places, water may enter sinkholes within a drainage basin, then reappear in springs downstream from a gage. Many of the smaller streams in the District go dry or have no flow for several months in many years. In addition to the low-flow statistics, four synoptic low-flow measurement surveys were conducted on 161 sites during 1990, 1995, and 1996. The measurements were made to provide "snapshots" of flow conditions of streams throughout the Suwannee River Water Management District. Magnitudes of low flows during the 1990 series of measurements were in the range associated with minimum 7-consecutive-day 50-year recurrence interval to the minimum 7-consecutive-day 20-year recurrence interval, except in Taylor and Dixie Counties, where the magnitudes ranged from the minimum 7-consecutive-day 5-year flow level to the7-consecutive-day 2-year flow level. The magnitudes were all greater than the minimum 7- consecutive-day 2-year flow level during 1995 and 1996. Observations of no flow were recorded at many of the sites for all four series of measurements.

Florida↗

USGS2024

New U-Pb geochronology and geochemistry of Paleozoic metaigneous rocks from western Yukon and eastern Alaska, cross-border synthesis, and implications for tectonic models

The tectonic evolution of and relation between the Yukon-Tanana terrane and the Lake George assemblage, as well as other associated tectonic assemblages in western Yukon and eastern Alaska, have been debated for decades. The Yukon-Tanana terrane is widely considered to be an allochthonous rifted fragment derived from the Laurentian continental margin, whereas the Lake George assemblage and associated assemblages are currently interpreted to be part of the parautochthonous continental margin of western North America (Laurentia). To address these topics, we present 40 new U-Pb zircon ages and 20 new whole-rock geochemical analyses. We incorporate these data into a new compilation of available geological mapping for a large area that straddles the Alaska-Yukon border, together with 34 previously published U-Pb age determinations and an extensive geochemical database of metaigneous rocks from Late Devonian to Early Mississippian and middle to late Permian assemblages in this area. Magmatism in the Lake George assemblage and related assemblages occurred in two pulses from about 371 to 360 and from about 358 to 347 million years ago (Ma); geochemical discrimination diagrams indicate a large crustal component, possibly indicative of arc magmatism, for felsic metaigneous rocks and a range of tectonic environments for mafic rocks. Magmatism in the Fortymile River and related assemblages, and parts of the Nasina assemblage—all parts of the Yukon-Tanana terrane—are mainly Early Mississippian and span a crystallization age range from about 361 to 343 Ma; geochemical discrimination diagrams for these rocks indicate primarily arc geochemical signatures for both mafic and felsic rocks. Middle to late Permian crystallization ages (about 261–253 Ma) are indicated for felsic metaigneous rocks in the Klondike assemblage and some of the felsic metaigneous rocks in the Nasina assemblage. Based on our mapping, we propose the existence of a possible unconformity between the Mississippian and Permian felsic metavolcanic rocks within the Nasina assemblage that is marked by sporadic occurrences of stretched-pebble conglomerate. Our combined database supports the well-established model of a magmatic arc comprising the Fortymile River and Finlayson assemblages of the rifted Yukon-Tanana terrane continental fragment on which a middle to late Permian arc (Klondike assemblage) was later built. The assemblages of the Yukon-Tanana terrane were subsequently intruded by Late Triassic to Early Jurassic granitoids, presumably during reaccretion of the Yukon-Tanana terrane to the continental margin. Permian and Late Triassic to Early Jurassic intrusions have not been mapped in the now structurally lower plate Lake George assemblage; their absence is one of the lines of evidence that have been used to support the parautochthonous, rather than allochthonous, origin of the Lake George assemblage and related assemblages. Our new data, together with previously published ranges of igneous crystallization ages and geochemical tectonic signatures of the Late Devonian to Early Mississippian magmatic rocks in the Lake George assemblage and associated assemblages and in the Fortymile River, Nasina, and correlated assemblages of the Yukon-Tanana terrane, indicate that the currently accepted interpretation of the Lake George assemblage and associated rocks being part of parauthochthonous North America is not the only possible interpretation of this tectonic entity. Approximately half of the dated intrusive rocks in the Lake George assemblage are contemporaneous with the metaigneous rocks of the Yukon-Tanana terrane arc (<361 Ma). We speculate that our approximately 361 Ma U-Pb age for quartz syenite in part of the North American continental margin in south-central Yukon defines the beginning of rifting of the Laurentian margin. Although the currently favored model of prolonged middle Paleozoic subduction and extension in both the Yukon-Tanana terrane and parautochthonous North America allows for simultaneous middle Paleozoic magmatism on both sides of the Slide Mountain Ocean, we now propose an alternative hypothesis in which the Lake George assemblage represents a deeper part of the rifted Yukon-Tanana terrane arc. If this is the case, the absence of Permian and Late Triassic to Early Jurassic arc rocks in the Lake George assemblage could be explained either by the arcs of these ages not being wide enough to have affected the Lake George assemblage or by tectonic displacement of these arc rocks away from the Lake George assemblage. Our approximately 259 Ma U-Pb zircon age and geochemical analyses of metarhyolite in the Seventymile terrane in Alaska, which comprises remnants of the back-arc basin that separated the Yukon-Tanana terrane from the Laurentian continental margin, confirm the presence of a late middle Permian volcanic arc component to the terrane. Our approximately 319 Ma U-Pb zircon age from the Chicken assemblage (as redefined in this study) in eastern Alaska, combined with previously reported fossil ages and a U-Pb zircon age from this assemblage, indicate that it is a Late Mississippian to Early Pennsylvanian arc assemblage. We propose several other relatively young, locally developed arc assemblages outboard of the ancient continental margin of Laurentia that may correlate with the Chicken assemblage, but we consider its origin to remain an enigma.

Alaska↗

USGS2026

Bedrock geologic map of the Eagle Lake quadrangle, Essex County, New York

The bedrock geology of the 7.5-minute Eagle Lake quadrangle, Essex County, New York, consists of deformed and metamorphosed Mesoproterozoic gneisses of the Adirondack Highlands unconformably overlain by weakly deformed lower Paleozoic sedimentary rocks of the Champlain Valley. The Mesoproterozoic rocks occur on the eastern edge of the Adirondack Highlands and represent an extension of the Grenville Province of Laurentia. Granulite facies Mesoproterozoic paragneiss, marble, and amphibolite hosted the emplacement of an anorthosite-mangerite-charnockite-granite (AMCG) suite, now exposed mostly as orthogneiss, at approximately 1.18–1.15 giga-annum (Ga, billion years before present). The earliest of four phases of deformation (D1) predated AMCG magmatism and is characterized by gneissosity, rarely preserved F1 isoclinal folds, and migmatite in the paragneiss host rocks. A sample of hornblende quartz syenite from the AMCG suite, collected from an abandoned railroad cut on Old Furnace Road, yielded a U-Pb zircon age of 1,149±10 million years before present. D2 deformation produced a composite penetrative gneissosity, migmatite, and isoclinal F2 folds. Towards the end of D2, felsic magmatism (including the regionally extensive Lyon Mountain Granite Gneiss, abbreviated “LMG”) spread by penetrative migration as semiconcordant alkali feldspar granite sheets subparallel to S2 into the previously deformed lithologies. The LMG crystallized at approximately 1.15 to 1.14 Ga and displays synkinematic F2 folds thus constraining the time of D2 deformation. Exhumation of the Marcy anorthosite began during D3 along a mylonitic extensional detachment, as a type of core complex. Protracted D3 produced F3 folds exhibited in regional domes and basins, such as the Hammondville antiform, reactivation of the S2 foliation, partial melting, metamorphism, metasomatism, iron ore remobilization, and intrusion of magnetite-bearing pegmatite both as layer-parallel sills and crosscutting dikes. D4 created NE- and NW-trending boudinage, local high-grade ductile shear zones, and crosscutting granitic pegmatite dikes. Kilometer (km)-scale lineaments readily observed in lidar data are Ediacaran mafic dikes and Phanerozoic brittle faults. Lower Paleozoic rocks are part of the Early Cambrian to Late Ordovician great American carbonate bank on the ancient margin of Laurentia. The Potsdam Sandstone preserves the Cambrian stratigraphy in outliers above the Great Unconformity. The Paleozoic rocks are weakly folded and block faulted. Parts of the quadrangle are covered by undifferentiated glacial deposits, but much of the quadrangle contains only a variably thick, veneer of unmapped glacial till over significant areas of exposed bedrock. The map also shows waste rock piles and locations of historical mining operations. This study was undertaken to improve our understanding of the bedrock geology in the Adirondack Highlands, establish a modern framework for 1:24,000-scale bedrock geologic mapping in the Adirondack Mountains, and provide a modern context for historical mines. This Scientific Investigations Map of the Eagle Lake 7.5-minute quadrangle consists of a map sheet, an explanatory pamphlet, and a geographic information system database that includes bedrock geologic units, faults, outcrops, and structural geologic information. The map sheet includes a bedrock geologic map, a correlation of map units, a description of map units, an explanation of map symbols, and two cross sections. The explanatory pamphlet includes a discussion of the geology.

New York↗

USGS2025

Base-flow sampling to enhance understanding of the groundwater flow component of nitrogen loading in small watersheds draining into Long Island Sound

Excessive nitrogen discharge is a major concern for the Long Island Sound. Programs have been implemented to reduce point sources of nitrogen to the sound, but little is known about the nonpoint sources. This study aims to better understand the current groundwater contributions of nitrogen from nonpoint sources in the Long Island Sound watershed. During the spring and summer of 2022, the U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency, collected water-quality samples to analyze nutrients (nitrogen and phosphorus), chloride, and bromide at 45 stations in the Long Island Sound watershed in Connecticut, New York, and Rhode Island. The stations were in small drainage watersheds (5 to 30 square kilometers) in the southern part of the Long Island Sound watershed. During two separate synoptic sampling events, water-quality samples and instantaneous streamflow measurements were collected under base-flow conditions (where the streamflow is dominated by groundwater inputs rather than overland flow or runoff flow). One sampling event was in the nongrowing season (April 24–25, 2022), and the other was in the growing season (June 30–July 1, 2022). To calculate instantaneous nitrogen loads and yields, streamflow was measured at the time of sample collection. Nitrogen concentrations, loads, and yields varied among sampling stations and by season. Total filtered nitrogen concentrations were generally lower in the nongrowing season (from less than 0.14 to 1.9 milligrams per liter) than in the growing season (from less than 0.23 to 3.0 milligrams per liter). Nitrate plus nitrite concentrations showed little variation between the nongrowing and growing seasons. Unfiltered ammonia plus organic nitrogen concentrations were generally lower in the nongrowing season (from less than 0.07 to 0.83 milligram per liter) than in the growing season (from 0.11 to 0.98 milligram per liter). In contrast, total filtered and unfiltered nitrogen loads and yields were higher in the nongrowing season than during the growing season, likely because streamflows were higher during the nongrowing season. Total unfiltered nitrogen yields during the nongrowing season ranged from less than 0.15 to 5.0 kilograms per square kilometer per day. Total unfiltered nitrogen yields during the growing season ranged from less than 0.12 to 2.5 kilograms per square kilometer per day. Total filtered nitrogen yields during the nongrowing season ranged from less than 0.13 to 5.2 kilograms per square kilometer per day. Total filtered nitrogen yields during the growing season ranged from less than 0.06 to 2.5 kilograms per square kilometer per day.

Connecticut, New York, Rhode Island↗

USGS2026

Regional models for postfire debris-flow likelihood and rainfall thresholds across the western United States

The U.S. Geological Survey (USGS) uses an empirical model developed with logistic regression (the ‘M1’ model) to rapidly assess debris-flow likelihood and to identify quantitative rainfall thresholds for debris flows after wildfire in the western United States. The M1 model was calibrated to a debris-flow inventory from southern California (United States) and has been applied throughout the western United States. Limited spatial coverage in the calibration dataset has motivated evaluation of M1 model accuracy outside the calibration region (e.g., the Sierra Nevada or the eastern Cascade Range, United States). Previous test cases showed that M1 overpredicts debris-flow likelihood and underpredicts rainfall thresholds for some locations (e.g., Arizona, northern California, Colorado, New Mexico, United States). We sought to improve the regional applicability of a debris-flow likelihood model by expanding the debris-flow inventory used for calibration, testing multiple potential models and generating an updated model framework. The updated inventory includes 3788 observations from 67 burned areas paired with short duration rainfall ratios. The updated model framework consists of a modified model structure and sets of coefficients calibrated separately to the entire updated inventory and to subsets of the inventory that intersect three Environmental Protection Agency (EPA) Level 2 ecoregions (Mediterranean California, Upper Gila Mountains and Western Cordillera). Comparisons of predictions from the updated models with observed rainfall and debris-flow activity show that the updated models outperform the M1 model by ~15%–60% and improve the uniformity of predictive performance across the western United States. The updated models also reduce false positive rates relative to M1 and generate rainfall thresholds that are better aligned with relative differences in regional climatology and debris-flow activity.

Arizona, California, Colorado, Idaho, Montana, Nev↗

USGS2025

Review and synthesis of the applications of machine learning to coalbed methane recovery

Over the last 30 years, a substantial literature has evolved on the use of machine learning (ML) to assess, predict, and improve the efficiency of coalbed methane (CBM) recovery. In the United States, the production of CBM declined as shale gas production matured, but CBM continues to be an important energy resource in other parts of the world. ML applications that have the potential to improve CBM reservoir management and production forecasts, and to increase exploration and operational efficiency, are still of significant interest. The integration of geostatistical techniques into the CBM ML applications has been largely absent but represents an opportunity for improvement. The literature demonstrates the widespread interest in, and applicability of, ML algorithms applied to CBM problems, and that they continue to result in improvements in predictive performance. However, (1) much of the research is more academic than operational, (2) many results are based on simulations, or small or proprietary datasets, (3) ML performance information can be inconsistent and sometimes entirely omitted, (4) most methodologies are unique to the specific CBM situation and likely not generalizable, (5) no standard data repositories are available to directly compare the performance of competing algorithms, and (6) the spatial component is often omitted. Finally, relatively new ML protocols involving causality analysis and reinforced learning, as well as hybrid workflows combining both supervised and unsupervised learning, are anticipated to dominate the future investigations. Integration of geostatistical and geospatial analysis with ML should enhance performance.

Book chapter↗

USGS2025

Beavers in the Tualatin River Basin, northwestern Oregon

Growing interest in beaver-assisted restoration in the Tualatin River Basin of northwestern Oregon motivated a series of studies by the U.S. Geological Survey to assess the capacity of the stream network to support beaver dams and to evaluate the effects of beaver dams and ponds on urban streams. This multichapter volume describes the data collection from 2016–17 and the findings of these studies, which were done in partnership with Clean Water Services. Chapter A documents the locations of beaver dams in the Tualatin River Basin and how many beaver dams the stream network could support with existing and improved riparian vegetation. Beaver dam capacity was estimated by modifying existing tools to account for the low gradient of many streams in the Tualatin River Basin. Chapter B describes the effects of beaver dams and ponds on hydrologic and hydraulic responses of storm flows. Hydrologic and hydraulic responses for two urban stream reaches were compared with and without beaver dams and ponds and for a range of streamflow conditions using two-dimensional hydraulic models. Chapter C characterizes the effects of beaver dams and ponds on the transport and deposition of suspended sediment. Continuous turbidity, discrete suspended-sediment samples, and streamflow measurements collected during storms and base-flow periods were used to assess: (1) suspended-sediment loads upstream and downstream from two beaver-affected reaches, and (2) seasonal and longitudinal turbidity patterns. Chapter D describes the effects of beaver dams and ponds on longitudinal, spatial, and seasonal water-quality patterns. Continuous and synoptic water-quality data were collected along urban stream reaches, and net ecosystem production was calculated for two beaver-affected reaches. The findings of these studies illustrate that the effects of beaver dams and ponds on hydrology, hydraulics, suspended-sediment transport and deposition, and water quality are dependent on the characteristics of a stream reach (for example, channel gradient, groundwater exchange, and riparian vegetation) and the characteristics of beaver dams and ponds along that reach. This information can be used to consider the implications of beaver-assisted restoration in the Tualatin River Basin and the effects of beaver dams and ponds in urban streams.

Oregon↗

USGS2025

Hydrogeologic framework of the Mountain Home area, southern Idaho

In the arid western Snake River Plain around the City of Mountain Home, Idaho, declining groundwater levels concern agricultural, municipal, and other water users who rely on groundwater for sustenance because surface-water resources are limited. The U.S. Geological Survey developed this hydrogeologic framework to provide an updated characterization of groundwater resources in the western Snake River Plain around the City of Mountain Home. The hydrogeologic framework comprises: (1) a conceptual description of hydrogeologic units, (2) a three-dimensional hydrogeologic model and borehole database, (3) a map of groundwater levels and change, and (4) a discussion of groundwater occurrence and movement within the study area. Hydrogeologic units were defined based on existing literature and the borehole database compiled for this study; the five hydrogeologic units are granite, rhyolite, basalt, fine-grained sediments, and coarse-grained sediments. Each unit can bear water, but the main regional aquifer in the study area occurs in the basalt and fine-grained sediment units with depth to water ranging from 150 to 765 feet. A perched groundwater zone near the City of Mountain Home is primarily hosted in basalt and used domestically with most depths to water ranging from 30 to 100 feet. Interflow zones, scoria, and vertical fractures create heterogeneity within the basalt hydrogeologic unit that exerts strong control on groundwater movement, creating horizontal perching conditions and zones of enhanced vertical conductivity that facilitate downward groundwater percolation. In the fine- and coarse-grained sediments and rhyolite units, inferred faults both impede and enhance groundwater movement. The borehole database was constructed by digitizing 540 well-driller reports and was used to build a three-dimensional hydrogeologic framework model which reasonably represents the spatial distribution of hydrogeologic units in the study area. Generally, fine-grained sediments underlie much of the study area, with basalt concentrated in the central and western study area and rhyolite and granite in the uplands to the north. Groundwater levels were measured in 180 wells in March and November 2023; these data were used to develop water-table contour maps and describe groundwater-level change over an irrigation season. Groundwater generally flows south-southwest to the Snake River and groundwater levels declined across most of the study area (from 0.03 to 22.01 feet) between spring and autumn 2023, which is consistent with long-term declines in the Cinder Cone Butte Critical Groundwater Area and Mountain Home Groundwater Management Area. Groundwater levels rose (0.6 to 15.44 feet) over the irrigation season in most wells in the perched groundwater zone near the City of Mountain Home and near the Snake River, indicating the importance of surface-water recharge to groundwater in areas where surface water irrigation occurs. In aggregate, this hydrogeologic framework provides an updated characterization of and new insights into groundwater resources in the study area to help inform water resources management.

Idaho↗

USGS2024

Flood-inundation maps for the Cuyahoga River at Jaite, Ohio, 2024

Digital flood-inundation maps for a nearly 6-mile reach of the Cuyahoga River at Jaite, Ohio, were created by the U.S. Geological Survey (USGS) in cooperation with the Northeast Ohio Regional Sewer District Board of Trustees. The maps depict estimates of the extent and depth of flooding corresponding to selected water levels (stages) at USGS streamgage 04206425 on the Cuyahoga River at Jaite, Ohio. Water-surface profiles were computed for the stream reach by using a one-dimensional steady-state step-backwater model. The hydraulic model was calibrated to the current USGS streamgage data and then used to compute 15 water-surface profiles for flood stages at 1-foot intervals referenced to the streamgage datum and ranging from 6 to 20 feet, which correspond to below “action stage” to “major flood stage” as reported by the National Weather Service. The simulated water-surface profiles were then used with a geographic information system digital elevation model derived from light detection and ranging data to delineate the areas flooded at each stage. These maps, along with current stage data from the USGS streamgage and forecasted high-flow stages from the National Weather Service, can provide emergency management personnel and residents with information that is critical for flood response activities such as evacuations and road closures, as well as for postflood recovery efforts.

Ohio↗

USGS2025

Geochemical evidence for the origin of late Quaternary loess, Seward Peninsula, Alaska

Loess is the most widespread surficial deposit in the state of Alaska. Although loess of last glacial age is common in mid-continental North America, records of last glacial loess in Alaska have been elusive. Here we report a record of last glacial loess on the Seward Peninsula, along with studies of a transect of loess localities also thought to date to the last glacial period. Highly diverse hypotheses on the origin of loess on the Seward Peninsula include local volcanic rocks and mountain ranges to the north and south, as well as offshore deposits. We present particle size and mineralogical and geochemical data that indicate that the most likely sources for loess in this part of eastern Beringia are glaciogenic silt particles of the Kobuk River and Noatak River, both of which drain the Brooks Range to the northeast. Northeasterly paleowinds were therefore responsible for loess on the Seward Peninsula. Such winds probably occurred during late summer or early autumn, when river discharge was at a minimum but temperatures were still above freezing. Northeasterly winds during the last glacial period are in good agreement with a recent paleoclimate model.

Alaska↗

USGS2000

Regional evaluation of evapotranspiration in the Everglades

Nine sites in the Florida Everglades were selected and instrumented for collection of data necessary for evapotranspiration-determination using the Bowen-ratio energy-budget method. The sites were selected to represent the sawgrass or cattail marshes, wet prairie, and open-water areas that constitute most of the natural Everglades system. At each site, measurements necessary for evapotranspiration (ET) calculation and modeling were automatically made and stored on-site at 15- or 30-minute intervals. Data collected included air temperature and humidity at two heights, wind speed and direction, incoming solar radiation, net solar radiation, water level and temperature, soil moisture content, soil temperature, soil heat flux, and rainfall. Data summarized in this report were collected from January 1996 through December 1997, and the development of site-specific and regional models of ET for this period is described. Latent heat flux is the energy flux density equivalent of the ET rate. Modified Priestley-Taylor models of latent heat flux as a function of selected independent variables were developed at each site. These models were used to fill in periods of missing latent heat flux measurement, and to develop regional models of the entire Everglades region. The regional models may be used to estimate ET in wet prairie, sawgrass or cattail marsh, and open-water portions of the natural Everglades system. The models are not applicable to forested areas or to the brackish areas adjacent to Florida Bay. Two types of regional models were developed. One type of model uses measurements of available energy at a site, together with incoming solar energy and water depth, to estimate hourly ET. This available-energy model requires site data for net radiation, water heat storage, and soil heat flux, as well as data for incoming solar radiation and water depth. The other type of model requires only incoming solar energy, air temperature, and water depth data to provide estimates of hourly ET. The second model thus uses data that are more readily available than the data required for the available-energy model. Computed ET mean annual totals for all nine sites for the 1996-97 period ranged from 42.4 inches per year at a site where the water level is below land surface for several months each year to 57.4 inches per year at an open-water site with no emergent vegetation. Although the density of photosynthetically-active plant leaves has been shown to relate directly to ET in some studies, it does not appear to relate directly to ET in the Everglades, based on comparison of annual ET data with leaf-area index, defined as the Normalized Difference Vegetation Index (NDVI), data from satellite imagery. NDVI and ET appear to be inversely related in the Everglades. The greatest ET rates occurred at open-water sites where the NDVI data indicated the lowest leaf-area index. Among the remaining vegetated sites, there is no clear relation between ET and NDVI, though the highest ET rate corresponded to the lowest NDVI and one of the lowest ET rates corresponded to the highest NDVI value. The variation in ET follows a seasonal pattern, with lowest monthly ET totals occurring in December through February, and highest ET occurring in May through August. The greatest range in monthly ET among all nine sites for the 2-year period occurred at site 3: from 1.81 inches in December 1997 to 6.84 inches in July 1996. A study to compare the Bowen-ratio/energy balance method of ET measurement with the eddy-correlation method was done at one site from June 22, 1998, through September 28, 1998. This comparison indicated that both methods gave comparable values of the Bowen ratio, but there was a considerable difference in available energy measured by the two methods. The mean of all 30-minute measured turbulent heat fluxes from the eddy-correlation apparatus for June 22 through September 29, 1998, was 137.4 watts per square meter, and the mean of the corresponding measured energy was 163.6 watts per square meter, or about 20 percent greater. The disagreement in mean energy fluxes measured by the two methods is problematical and is not fully understood. Although the difference seems to be related to friction velocity, and is practically non-existent at values of friction velocity greater than 0.3 meter per second, the "correctness" of either method cannot be determined with the data available.

Florida↗