USGS Science⌕ Search

SEARCH · USGS Science

Results for “International Journal of Coal Geology”

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,660 records · Page 74Linked to original sources
USGS2025

Zircon as a pathfinder to REE mineralization

Carbonatites and alkaline silicate rocks are major primary sources of the rare earth elements (REE) and other critical metals, such as Nb. Despite the economic significance of these rocks, their formation and the processes of REE enrichment are poorly understood. Here, statistical analysis of a global dataset demonstrates that zircon geochemistry is a powerful recorder of REE metallogenesis and a potential pathfinder for REE deposits. Zircons from REE and Nb fertile intrusions lack Eu anomalies and have elevated Gd/Yb and Th/Yb, indicating they crystallised from magmas that originated from deep, oxidised and enriched mantle sources. Complexes with Nb enrichment have low U/Nb, reflecting an enriched mantle source, whereas high U/Nb in REE-only fertile intrusions suggest a subduction-metasomatised mantle source. Machine learning models demonstrate high accuracy in classifying zircon from barren and fertile deposits. Classification of detrital zircons shows that REE-enriched deposits correlate with supercontinent assembly, whereas Nb fertile complexes are associated with supercontinent breakup. This approach offers a new, mineral to global scale, petrologic and exploration tool that enhances understanding of REE metallogenesis.

Geochemical Perspectives Letters↗
USGS2025

Mapping eelgrass (Zostera marina) cover and biomass at Izembek Lagoon, Alaska, using in-situ field data and Sentinel-2 satellite imagery

The U.S. Geological Survey and the U.S. Fish and Wildlife Service have developed a three-tiered strategy for monitoring eelgrass ( Zostera marina ) beds at Izembek Lagoon, Alaska, that targets different spatial and temporal scales. The broadest-scale monitoring (tier-1) uses satellite imagery about every 5 years to delineate the spatial extent of eelgrass beds throughout the lagoon. This report describes the most recent (mid-2020s) tier-1 eelgrass monitoring at Izembek Lagoon. The monitoring effort began by canvasing all satellite imagery collected during summer, under clear daytime skies and at low-tide, since the last tier-1 effort in 2006. Two eelgrass maps of Izembek Lagoon were generated by first creating maps of spectrally unique classes from two Sentinel-2 satellite images collected on July 1, 2016, and August 14, 2020, then attributing those spectral classes with information about eelgrass conditions based on field data. Specifically, maps depicting various eelgrass metrics, such as percentage of cover and modeled biomass, were generated using summaries of the ground data that spatially intersected each spectral class. Comparisons of the 2016 and 2020 Sentinel-2 maps showing eelgrass distributional extent, as well as a 2006 Landsat map, indicated that areas where eelgrass presence may have declined during 2006–20 were most prevalent in the central part of Izembek Lagoon. More recently, during 2016-20, areas of possible biomass decline were more prevalent in the southern part of the lagoon. Monitoring eelgrass conditions at Izembek Lagoon with satellite imagery and concurrent ground data allows conditions to be compared over time, but the influences of tide levels, growing season phenology, and spatiotemporal co-registration accuracy should be considered when designing and interpreting change detection analyses.

Alaska↗
USGS2024

Case definitions for wildlife diseases

Welcome to the first manual of “Case Definitions for Wildlife Diseases,” a “living” electronic publication. The plan is to add and update this manual’s case definitions periodically as warranted; thus, this manual will never be completed, and readers should download the latest versions of specific chapters (that is, definitions) when available. Constructive suggestions from readers are welcome and will help guide adjustments as this project progresses. The purpose of this manual is to provide case definitions for selected diseases of importance to wildlife in Canada and the United States. Case definitions provide standard sets of criteria for classifying the degree of certainty of a particular diagnosis and help improve surveillance data quality and comparability. Better data and standardization allow for improved data sharing, which increases geographic and species surveillance coverage and permits more robust analyses. The definitions included in this manual have been developed by veterinary pathologists, epidemiologists, and wildlife biologists primarily from the U.S. Geological Survey National Wildlife Health Center (NWHC) and Canadian Wildlife Health Cooperative (CWHC). Pathologists from each organization reviewed and finalized the definitions. Each case definition has been peer reviewed by two scientific experts before publication. This manual begins with the case definition template. This generic template includes four sections: “Individual, Place, and Time Criteria for Diagnosis and Testing,” “Field Criteria for Diagnosis,” “Laboratory Criteria for Diagnosis,” and “Epidemiological Linkage Criteria for Diagnosis” and can be used to guide development of new case definitions. Information in each section is then combined to provide an overall case classification. Disease diagnoses are classified as “Confirmed,” “Presumptive,” or “Suspected;” and evidence of a pathogen or toxin is classified as “Exposed” or “Present/Detected.” Each subsequent chapter is then a case definition for a specific disease of wildlife, and infectious and non-infectious diseases are included.

Techniques and Methods↗
USGS2025

Mapping potential sensitivity to hydrogeomorphic change in the UMRS riverscape

In 2020 the U.S. Geological Survey (USGS), as part of the Upper Mississippi River Restoration (UMRR) Program, began a new project to characterize potential hydrogeomorphic change associated with hydrogeomorphic units (HGUs) and their catenae (units linked by their association with sediment sources and flow origins). The goal of the project was to develop a geographic information system (GIS) database of HGUs for the Upper Mississippi River System (UMRS) available to both scientists and river managers working on UMRR studies and HREP planning and design studies. The characterization was based on a hydrogeomorphic change hierarchical classification developed previously for the UMRS. The products were generated with automated techniques in a GIS using systemic datasets. Landforms were mapped from the 2015 UMRS topobathymetric dataset with geomorphon (shape-based) tools tailored for the large riverscape, valley bottom environments in the UMRS. A clustering analysis was applied to the resulting landforms to identify HGUs associated with zones of perennial low flows, bankfull flows, and overbank floodplains. Catenae were assembled based on the proximity of the units to the main channel, tributary mouths, and side channels from previously published aquatic areas (USACE, 2018) coupled with least-cost flowpath linkages between potential sediment origins and planform change units developed by Rogala, Fitzpatrick, and Henderson (2020). These GIS-based analyses were successful at identifying a range of HGUs using an automated technique with available data across the entire riverscape, with emphasis on those that have the potential for hydrogeomorphic change. Most of the resulting features are depositional, as expected in a large river system. However, this is the first attempt of linking tributary inputs, side channel erosion and levee breaches with their depositional counterparts. The approach was successfully piloted in Pools 8 and 10 in the Upper Impounded Reach and Pool 14 in the Lower Impounded Reach, with next steps for application in reaches of the unimpounded section and Illinois River. This report emphasizes results from Pool 10, which was the focus of most of our attention during the pilot phase.

Illinois, Indiana, Iowa, Minnesota, Missouri, Sout↗
USGS2025

S/P amplitude ratios with Distributed Acoustic Sensing and application to earthquake focal mechanisms

Distributed acoustic sensing (DAS), which transforms a fiber optic cable into an array of high frequency strainmeters, has the potential to help us characterize earthquakes with a dense sampling of measurements. While earthquake focal mechanisms are frequently determined using P -wave polarities and S / P amplitude ratios with inertial seismometers, the dense sampling of DAS over potentially large portions of the focal sphere can aid our solutions. Here, we consider ∼200 regional earthquakes and thousands of S / P measurements on a DAS cable colocated with a network of inertial seismometers near Arcata, California. We demonstrate the S / P ratio measurements made on DAS are similar to those made on the inertial seismometers and can be used to constrain focal mechanism solutions.

California↗
USGS2024

Substrate Enhancement Pilot Project—Monitoring summary and evaluation, Kootenai River, Idaho, 2012–22

To assess changes in substrate conditions and the efficacy of artificially placed substrates at select sites on the Kootenai River near Bonners Ferry, Idaho, the U.S. Geological Survey, in cooperation with the Kootenai Tribe of Idaho, completed repeat bathymetric, velocimetric, and underwater videography surveys. Collectively, three project sites throughout the Kootenai River make up the Substrate Enhancement Pilot Project (SEPP), an effort intended to improve spawning and egg incubation viability at locations identified to be aquatic habitat limited for the endangered Kootenai River white sturgeon (Acipenser transmontanus). Following the placement of coarse substrates at each site, bathymetric, velocimetric, and underwater videography data were collected from 2012 to 2022 to assess the role of deposition and erosion on maintaining suitable white sturgeon spawning and incubation substrates. Minimal erosion and deposition occurred at all Substrate Enhancement Pilot Project sites, according to interannual and intra-annual net volumetric changes between bathymetric surveys. Depending on the timing of bathymetric surveys relative to the annual peak streamflow conditions, isolated locations of deposition or erosion were observed at each site and the potential for deposition or erosion was supported by measured mean depth-averaged velocities. This study concluded that variability of deposition and scour were common at each site throughout the monitoring period and may be attributed to fluctuations in streamflow. Repeat bathymetric, underwater videography, and velocity mapping surveys were used to verify the interstitial spaces and surfaces of substrates at each SEPP site remained free of fine sediments for intervals longer than a year but were susceptible to deposition between high streamflow events.

Idaho↗
USGS2025

Do bighorn sheep use desert tortoise burrow spoil piles as mineral licks in southern California?

In a previous study, bighorn sheep ( Ovis canadensis ) were periodically photographed by trail cameras when they visited desert tortoise ( Gopherus agassizii ) burrows in the San Bernardino Mountains near Palm Springs, California, USA. The authors suggested that bighorn sheep may utilize the excavated spoil material from burrows as a mineral lick. To test that supposition, we collected soil samples from the spoil piles and undisturbed surface soil up to 1 m away at two desert tortoise study sites with bighorn sheep at the aforementioned site and another nearby site. We hypothesized that sodium concentrations would be greater in spoil than in nearby surface soil samples. Common major cations were quantified in soil samples using inductively coupled plasma-mass spectrometry following microwave-assisted digestion to solubilize the fraction of elements that could become biologically or environmentally available. Mean sodium concentrations were greater in burrow spoil samples than in undisturbed soil samples at both sites, likely due to aqueous leaching of sodium from surface to subsurface soil horizons. Although we did not collect evidence of geophagy by bighorn sheep, the behavior is well known. However, trail camera photographs of bighorn sheep with their head next to the entrance of tortoise burrows suggests that they may periodically utilize sodium and other elements brought to the surface by desert tortoises that are otherwise unavailable to bighorn sheep on the surface.

California↗
USGS2024

Ungulate migrations of the Western United States, volume 4

Broadly distributed across the Western United States, ungulates (hooved mammals) play an important role in ecosystem function by affecting vegetation communities and forming the prey base for large carnivores. Additionally, ungulates provide economic benefits to regional communities through tourism and hunting and hold cultural significance for many Tribal communities. Many ungulates migrate seasonally between distinct summer and winter ranges to take advantage of spatially and temporally variable food sources and avoid threats such as predators and deep snow. Increasingly, these migrations are threatened by the growing human footprint and associated subdivisions, energy development, and increased traffic volume. Efforts to study ungulate populations and conserve their migrations received support in recent years from the U.S. Department of the Interior Secretarial Order No. 3362, which provided Federal support for enhancing habitat quality for ungulates across the Western States. In response to Secretarial Order No. 3362, the U.S. Geological Survey (USGS) established the Corridor Mapping Team, a collaboration among USGS and participating State and Federal wildlife management agencies and numerous Tribal Nations. Together, the Corridor Mapping Team maps ungulate migrations throughout the Western United States in the USGS “Ungulate Migrations of the Western United States” report series. This report (volume 4) details migrations and seasonal ranges from 31 new herds throughout nine Western States. Additionally, this report includes updates to two herds published in previous reports. Including this report, the report series has provided the mapped migrations and seasonal ranges of 182 unique herds and has provided a map-based inventory of the documented ungulate migrations across the Western United States for biologists, managers, policy makers, and conservation practitioners. This report also discusses how the mapping efforts associated with the Corridor Mapping Team can be used to guide management and policy regarding renewable energy development and ungulate disease, specifically chronic wasting disease, in the Western United States.

Arizona, California, Colorado, New Mexico, Nevada,↗
USGS2026

Site response models based on geometric parameters for southern California sedimentary basins

Site response in sedimentary basins is influenced by complex three-dimensional (3D) features, including trapping of seismic waves, focusing of seismic energy and basin resonance. Current ground motion models (GMMs) incorporate basin effects using one-dimensional parameters like V S30 and shear wave velocity isosurface depths, which are limited in capturing lateral and 3D effects. To address these limitations, we develop seismic site response models based on novel parameters that represent multi-dimensional properties of the Los Angeles Basin (LAB) geometry and shear wave velocity. We define a basin shape for the LAB using depth to subsurface geologic interfaces associated with the oldest sedimentary deposits (depth to a particular shear wave velocity horizon, i.e., 1.5 km/s - z 1.5 ) and the depth to the crystalline basement ( z cb ) which are determined using geologic cross sections and community seismic velocity model profiles. We explore a suite of geometric descriptors computed for the LAB and southern California, from which three parameters with the greatest predictive potential are selected and evaluated using empirical ground motion residual analyses in combination with the Boore et al. GMM. The results demonstrate that the zonal heterogeneity index ( ), standard deviation of the absolute difference between z 1.5 and z cb ( ) and standard deviation of z cb ( ) each provide a reduction in site-to-site variability ( ϕ S2S ) of empirical GMMs. The reduction in ϕ S2S is period-dependent, with average decreases of 3%, 26% and 6% for , , and , respectively. Although these reductions are modest from an engineering application perspective, they are statistically significant, underscoring the inherent difficulty in fully characterising complex basin effects. Collectively, these findings indicate that the inclusion of basin-specific geometric parameters yields measurable, albeit incremental, improvements in site response prediction and establishes a framework for the progressive refinement of seismic hazard characterisation within sedimentary basins.

California↗
USGS2025

Three-dimensional seismic velocity model for the Cascadia Subduction Zone with shallow soils and topography, version 1.7

The U.S. Geological Survey’s seismic velocity model for the Cascadia Subduction Zone provides P- and S-wave velocity ( V P and V S , respectively) information from 40.2° to 50.0° N. latitude and −129.0° to −121.0° W. longitude, and is used to support a variety of research topics, including three-dimensional (3D) earthquake simulations and seismic hazard assessment in the Pacific Northwest. This report describes an update to the previous version (v) 1.6 of the 3D seismic velocity model for the Cascadia Subduction Zone. This new model (herein referred to as v1.7) contains more detailed near-surface structure for improved earthquake ground motion modeling. Updated features include the addition of a new shallow soil velocity model in the top few hundred meters and the option of adding user-specified topography. Although v1.6 of the Cascadia seismic velocity model has a minimum V S of 600 meters per second (m/s), the new model (v1.7) has a minimum V S of approximately 40 m/s. Overall, this update will allow for more accurate ground motion estimates from 3D simulations of scenario earthquakes in the Cascadia Subduction Zone region.

British Columbia, California, Oregon, Washington↗
USGS2026

Modeled groundwater and surface-water interactions surrounding Mobile Bay, Alabama, 2008–15

The U.S. Geological Survey, in cooperation with the Gulf Coast Ecosystem Restoration Council, has used MODFLOW 6 to develop a groundwater-flow model to simulate groundwater and surface-water interactions in the Mobile Bay, Alabama, area. The model results indicated that, on average, groundwater discharge near the coastline is equal to 2.5 percent of the surface water that flows into the bay. The model was also used to determine how changes in recharge, sea level, and groundwater pumping affect groundwater levels and discharge rates. The results indicate that more groundwater discharge occurred in the winter and spring when recharge was higher, sea level was lower, and groundwater pumping was lower than during the summer and fall. Additionally, the amount of emergent groundwater was closely related to sea level; when sea level was higher, there was more area with emergent groundwater. Furthermore, the depth of nonemergent groundwater was related to trends in recharge and pumping. During periods of increased recharge and reduced pumping, a significant portion of the model area exhibited a depth to the water table of less than 1 meter below the land surface. Conversely, during periods of decreased recharge and heightened pumping, much of the model area showed a depth to the water table ranging from 1 to 5 meters.

Alabama, Mississippi↗
USGS2026

Regreening, restoring, and reconnecting a southwestern wetland ecosystem – the Zeedyk wetland

Alluvial wetland ecosystems are vital as biodiversity hotspots but are increasingly threatened by anthropogenic stressors and drought. These pressures are especially acute in arid and semi-arid regions, where eco-hydrologic connectivity is fragile and recovery is slow. This study quantifies the efficacy of nature-based solutions, particularly the ‘Zeedyk approach,’ which employs low-tech Natural Infrastructure in Dryland Streams (NIDS)—including rock detention structures—to slow surface water, raise groundwater tables, and restore wetland function at a spring-fed wetland in Cebolla Canyon, New Mexico, U.S.A. Our results depict a Restoration Feedback Loop that captures stages of change from a healthy wetland in 1935, altered by 20th-century agriculture and grazing, to the re-establishment of the historical flow regime by 2024 documented through an 89-year archive of aerial imagery (1935–2024). By the end of our study period, the Spring-Fed Wetland had expanded by roughly 229% of the original 1935 area, to 4.13 ha. Using 40 years of satellite data, we assess changes in vegetation and hydrology with remote sensing indices. Spatial and temporal analyses reveal significant increases in vegetation greenness and wetness, particularly in an Expanded Wetland subregion, which exhibited ∼3.5x higher wetness and ∼1.5x higher greenness trends compared to adjacent areas. Monthly metrics highlight seasonal variability, with increases in greenness linked to monsoonal rainfall and lateral water redistribution, indicating that restoration impacts extend beyond the primary wetland. This study demonstrates the utility of cloud-based platforms like Google Earth Engine and USGS EarthExplorer for long-term monitoring of wetland restoration, while quantifying the efficacy of the ‘Zeedyk approach’ and demonstrating its potential as a scalable method to restore and conserve wetland meadows in other arid and semi-arid landscapes.

New Mexico↗
USGS2026

Lunar analog study using portable gamma-ray neutron detector: Radiochemical mapping of silicic-to-basaltic volcanic terrains in the San Francisco Volcanic Field, Arizona

Studying lunar silicic volcanism provides key insights into the Moon’s crustal evolution, magmatic processes, and volcanic history. The Lunar-VISE (Lunar Vulkan Imaging and Spectroscopy Explorer) mission will investigate the Gruithuisen domes, a unique lunar region hypothesized to have formed through silicic volcanism. Using instruments on a Firefly Aerospace lander and a Honeybee Robotics rover, Lunar-VISE will analyze mineralogy, geochemistry, and surface properties to determine the origin and evolution of the domes, with a gamma-ray and neutron spectrometer (LV-GRNS) among its payload instruments. In preparation for this mission, we conducted preliminary fieldwork using a handheld gamma-ray neutron detector with NaI(Tl) and Cs LiYCl :Ce (CLYC) scintillators. This study focuses on various rhyolitic and basaltic volcanic centers in the San Francisco Volcanic Field (SFVF) near Flagstaff, Arizona—specifically Sugarloaf Peak, Bonito Lava Flow, and Robinson Mountain, a region containing several well-characterized lunar analog sites. The SFVF was selected for this study due to its broad compositional diversity spanning basaltic to rhyolitic compositions, its well-preserved volcanic morphologies, and its extensive use in previous NASA field campaigns and astronaut training exercises, making it an ideal terrestrial laboratory for testing planetary exploration techniques. We measured natural radioactivity, specifically from potassium (K), thorium (Th), and uranium (U) and other elements in their decay chains, which serve as diagnostic tracers of magmatic differentiation and crustal evolution processes, to assess geochemical variability across compositionally diverse terrains. The detector’s sensitivity was assessed across varying concentration levels. Regionally averaged concentrations of radioisotopes were determined by gamma-ray spectroscopy at selected sites. Our measurements reveal significant compositional variations between sites, with Sugarloaf Peak (rhyolitic, silicic dome) exhibiting the highest average radioisotope concentrations (K: 3.29 ± 0.24 wt%, U: 14.35 ± 1.81 ppm, Th: 27.14 ± 2.43 ppm), while Bonito Lava Flow (K: 1.13 ± 0.08 wt%, U: 3.86 ± 0.52 ppm, Th: 7.52 ± 1.02 ppm), and Robinson Mountain (K: 1.46 ± 0.16 wt%, U: 5.59 ± 1.06 ppm, Th: 11.75 ± 1.98 ppm) show lower concentrations aligned with basaltic compositions. Gamma-ray fluxes were elevated by a factor of approximately three to four at Sugarloaf Peak relative to nearby basaltic terrains, consistent with expected geochemical differentiation patterns. Furthermore, at meter scales, proximity to geological features significantly affects measurements. Th concentrations adjacent to a cliff face at Sugarloaf’s base were 28% higher than values measured 3-5 meters away from the same feature, demonstrating localized compositional heterogeneity. Analysis of station-by-station measurements reveals that K, U, and Th concentrations show an elevation trend at Robinson Mountain and generally higher concentrations near Sugarloaf Peak’s summit, suggesting progressive magmatic differentiation and/or the presence of more evolved lithologies at higher elevations. By mapping these elements on Earth using GRNS, we aim to optimize measurement techniques for Lunar-VISE and similar rover-borne missions by demonstrating the applicability of portable GRNS instruments for planetary surface exploration in an analog environment.

Arizona↗
USGS2025

Earthquake stress-drop values delineate spatial variations in maximum shear stress in the Japanese forearc lithosphere

Earthquake stress drop (Δσ) may increase with depth and stress in the brittle lithosphere. However, the range of uncertainty in Δσ and the lack of constraints on absolute stress make it difficult to establish whether they are correlated. Here, we investigate Δσ dependence on depth and maximum shear stress ( τ max ) based on ~11 years of seismicity in the northeastern Japanese forearc following the 2011 Tohoku-Oki megathrust earthquake. We interpret Δσ estimates computed using both individual spectra and spectral-ratio methods and find that Δσ exhibits a clear depth dependence within the seismically active upper ~60 km of the forearc lithosphere ( ~ 0.8 MPa per 10 km). We further compare Δσ values with quantitative τ max estimates from finite-element models of force balance. We find that median Δσ values increase with τ max in the brittle forearc lithosphere and that earthquake stress release is proportional to τ max . The dependence of Δσ on τ max explains the apparent depth dependence of Δσ and suggests that average Δσ values provide a relative measure of the stress at failure. In the northeastern Japanese forearc, Δσ values remained roughly constant in the decade following the Tohoku-Oki earthquake, suggesting negligible changes in failure stress in the forearc since the mainshock.

Communications Earth and Environment↗
USGS2025

Airborne geophysical efforts for critical mineral systems mapping in the southern Midcontinent, USA

The increasing global demand for critical minerals to support energy and technological advancement has accelerated exploration and research efforts for these essential resources. Since 2019, the United States Geological Survey (USGS) Earth Mapping Resources Initiative (EMRI) has worked to modernize geologic mapping of the Nation to better understand its critical mineral resources. To further this initiative, the USGS has flown a series of high-resolution airborne magnetic and radiometric surveys over large areas of the southern Midcontinent. The surveys cover known critical mineral deposits and areas with the potential to host additional critical minerals based on the presence of one or more overlapping mineral systems. One aspect of EMRI emphasizes close collaboration between the USGS and the Association of American State Geologists, as well as other government and industry partners to leverage geophysical, geological, and geochemical expertise on both regional and local scales. The EMRI high-resolution airborne survey data provide new insights into the geophysical framework of the southern Midcontinent and its critical mineral endowment. Additionally, discoveries made from the data have directed new studies for critical mineral exploration.

southern Midcontinent↗
USGS2025

Assessing flood water infiltration and storage in a restored floodplain

In urban areas, floodplain restoration is gaining prominence as a strategy for restoring the natural functions of floodplain ecosystems and reducing flood risk. This has spurred research into potential interactions between floodwaters, the hyporheic zone, and the floodplain aquifer. An urban restored stream in Wisconsin, USA, was used as a case study to examine four methods to estimate floodplain infiltration and storage during overbank floods. We characterised flood-related infiltration over a 4-year period from 2018 through 2021 by simultaneously and continuously measuring groundwater levels and vertical temperature profiles with stream water levels linked to high-resolution flood inundation maps. High-resolution topographic data helped to quantify surface floodplain storage and the unsaturated soil volume relative to flood stage. Infiltration estimates from the simple methods align well with those from the more complex methods; however, the complex methods provide additional insights about the factors influencing infiltration. Results from all methods indicate that the volume of water that vertically infiltrates during floods is likely small relative to the total volume of the flood, with 0.08%–0.52% of flood water infiltrating into the floodplain, on average. Spatially variable vertical hydraulic gradients, driven by flood depth, groundwater level, and permeability, imply heterogeneous patterns of infiltration across the floodplain. Gradients favourable for infiltration typically occurred during the onset of flooding but, over the study period, were mostly (98% of the time) favourable for groundwater discharge to the channel (non-flood periods). These findings highlight the importance of considering surface-groundwater dynamics, floodplain soils, and unsaturated floodplain volume in defining the benefits of floodplain infiltration for flood attenuation.

Wisconsin↗
USGS2024

Inset groundwater-flow models for the Cache and Grand Prairie Critical Groundwater Areas, northeastern Arkansas

The water resources in the Mississippi alluvial plain, located in parts of Missouri, Kentucky, Tennessee, Mississippi, Louisiana, and Arkansas, supports a multibillion-dollar agricultural industry that relies heavily on pumping of groundwater for irrigation of crops and aquaculture. The primary source of groundwater for agricultural-related pumping is the Mississippi River Valley alluvial aquifer, which has declined in storage for decades; secondary groundwater sources include the middle Claiborne aquifer and Wilcox aquifer system. Two areas in northeastern Arkansas that lie within the Mississippi alluvial plain, part of the Cache and Grand Prairie regions, have been designated as Critical Groundwater Areas owing to decades of groundwater declines that resulted from past and current water use. The multidisciplinary Mississippi Alluvial Plain project, led by the U.S. Geological Survey, and funded by their Water Availability and Use Science Program, included objectives to develop numerical groundwater models in focus regions, including the part of the Cache and Grand Prairie regions of northeastern Arkansas. Two inset models were developed using the child model capabilities of MODFLOW 6, the U.S. Geological Survey’s Modular Hydrologic Model simulation software. Both models, called the Cache model and Grand Prairie model, simulated the groundwater system and surface-water/groundwater interactions for the Mississippi River Valley alluvial aquifer and underlying Tertiary-age aquifers and confining units to the Midway confining unit. Each model was spatially discretized into 500-meter x 500-meter orthogonal cells on a grid with 5-meter constant-thickness vertical layers that represented the Mississippi River Valley alluvial aquifer and increasing thickness layers for the aquifers and confining units below the alluvial aquifer. The Cache and Grand Prairie models were calibrated with the PEST++ iterative ensemble smoother Version 5 and employed high dimensional parameterization schemes of 13,740 and 30,436 parameters, respectively. The Cache mean absolute residual for groundwater-level observations within each model domain for the priority well was 1.58 meters. Grand Prairie mean absolute residuals for the alluvial aquifer and middle Claiborne aquifer groundwater-level observations were 2.71 and 10.78 meters, respectively. The groundwater budgets for the Cache and Grand Prairie models were characterized by substantial outflows to irrigation wells, which constituted about 52 and 54 percent of all outflows, with the primary source of water to those wells being releases from unconfined aquifer storage.

Arkansas↗