USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Maps”

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,666 records · Page 64Linked to original sources
USGS1976

Uranium and trace elements in stream sediments as an exploration tool

More than 45 trace elements have been reported in anomalous amounts in the uranium ores. The specific suite of elements associated with any one uranium deposit varies according to deposit type and geologic province. The primary geochemical halo of uranium and associated trace elements in the host rock, together with secondary dispersion halo in soils and alluvium, offers a potential geochemical exploration target. Sediment from streams near low-grade uranium occurrences in arkosic sandstone of the Denver basin, Colorado, and sediment from streams in the igneous and metamorphic terrace near the Midnite mine of eastern Washington have been analyzed for acid-extractable U, Cu,Pb, Zn, Ni, Co, Fe, Mn, Mo, Cd, Cr, Ag, V, and Se, and for total Fe, Mn, Hg, As, and organic C. The -80 mesh fraction of the alluvium downstream from a small uraniferous limonite occurrence in the Denver basin shows anomalous concentrations of U, Pb, As, Ni, Mo, V, Zn, and Se. Anomalous concentrations of U persist farther downstream (as much as 1,000 m) than any of the trace elements. In another stream system southeast of Denver, above-background levels of U, Pb, As, Zn, V, Fe, and Mn occur near a uranium anomaly (defined by other methods) which has no surface mineralization. In the Midnite mine area, sediments of streams draining prospect areas showed anomalous concentrations of uranium (as much as 10 times background), but no accompanying anomalous trace-element concentrations, although analyses of Midnite mine ores show that a broad suite of trace elements are associated with uranium.

Open-File Report↗
USGS2024

Correlation analysis of groundwater and hydrologic data, Kaloko-Honokōhau National Historical Park, Hawai‘i

Designated in 1978, Kaloko-Honokōhau National Historical Park is located on the west coast of the Island of Hawaiʻi. The Kaloko-Honokōhau National Historical Park encompasses about 1,200 acres of coastal land and nearshore ecosystems, which include wetlands, anchialine pools (landlocked bodies of brackish water with hydrologic connections to the ocean), fishponds, a fishtrap, and coral reefs. These nearshore ecosystems are dependent on groundwater discharge with a freshwater component and provide habitat for threatened and endangered, endemic species, such as the orangeblack Hawaiian damselfly ( Megalagrion xanthomelas ) and the Hawaiian coot (ʻAlae keʻokeʻo, Fulica alai ). The populations of these native species, however, are threatened because of habitat loss related to urban development and environmental changes. Kaloko-Honokōhau National Historical Park is within the Keauhou aquifer system and the North Kona District, which experienced a 52 percent resident-population increase between 2000 and 2020 and a 41 percent visitor increase between 2008 and 2019. To support the current water demand associated with this growing population, groundwater is the primary source of freshwater used in the North Kona District, with about 15 million gallons of groundwater withdrawn from the Keauhou aquifer system per day since 2009. With anticipated development, future (2015–35) groundwater withdrawal from the Keauhou aquifer system is projected to be about 55 percent greater than recent (2012–14) withdrawal. Because Kaloko-Honokōhau National Historical Park is located within a coastal aquifer, natural and human-induced changes can affect the quality and quantity of groundwater, which can threaten groundwater-dependent ecosystems. To improve understanding of recent groundwater conditions, the U.S. Geological Survey, in cooperation with the National Park Service, undertook this study to document correlations between hydrologic time-series datasets from sites in and near Kaloko-Honokōhau National Historical Park using the nonparametric (distribution-free) Kendall’s tau statistical test. For the statistical analyses, dependent variables representing the groundwater system include groundwater level, the groundwater-level difference between pairs of sites, and specific conductance, and independent variables include datasets of sea level, rainfall, and groundwater withdrawal. About 34 percent of the 140 non-time-lagged Kendall’s tau statistical tests evaluated in this report are statistically significant ( p -value ≤ 0.050) with generally weak (0.1 ≤ tau ≤ 0.2) to moderate (0.2 ≤ tau ≤ 0.3) correlations. Groundwater levels measured at monitoring sites have the strongest correlation with the multivariate El Niño–Southern Oscillation index and withdrawal from production wells at the nearby Kohanaiki Private Club Community. Specific conductance is not consistently and significantly correlated with the independent hydrologic variables investigated in this report. Because the relations between hydrologic variables are commonly not instantaneous, a second set of correlations was evaluated after applying a range of time lags to the independent variable datasets. Relative to the non-time-lagged case (the set of correlations that did not use time-lagged independent variables), some of the time-lagged independent variables improved correlations with some of the dependent variables. For a particular independent variable, similar time lags were expected between the independent variable and dependent variable at all four monitoring sites. However, different time lags among the four sites sometimes produced the strongest correlations. This study identified several correlations that are statistically significant and hydrologically plausible, but the correlations could indicate that multiple concurrent factors are controlling the observed groundwater-system response, which might be better addressed using multivariate analyses. This study only investigates bivariate correlations, which may not explain all the variance in the data. The correlations analyzed in this report are limited by the quantity of available hydrologic data in the area near Kaloko-Honokōhau National Historical Park and are based on 14 years of time-series data, which were aggregated to a relatively coarse monthly temporal resolution that represents the minimum resolution common to all datasets.

Hawaii↗
USGS2026

Temporal and spatial changes in seismic attenuation associated with inferred fluid migration in the 2016 central Apennines earthquake sequence

Prior work suggests that high‐frequency seismic attenuation acts as a highly sensitive proxy for crustal permeability and fluid mobility in fractured media. We test the hypothesis that the fault system responsible for the 2016–2017 Amatrice–Visso–Norcia–Capitignano sequence acted as an impermeable seal, compartmentalizing pressurized fluids until dynamic rupture triggered widespread fluid diffusion. By tracking across the sequence the spatiotemporal evolution of the S ‐wave anelastic attenuation parameter, we identify large, positive low‐frequency attenuation anomalies emerging within the hanging wall following the Amatrice mainshock and strictly preceding subsequent large ruptures. Conversely, we observe weaker, negative anomalies in the footwall, anticorrelated in time with those of the hanging wall, revealing a massive asymmetry in fluid redistribution and permeability evolution across the fault system. Furthermore, aftershock migration rates reveal distinct linear alignments in a distance‐reduced time space, allowing us to explicitly track and quantify episodes of lateral and upward fluid migration. These physically consistent patterns suggest that stress‐driven fluid diffusion directly weakens adjacent fault patches, dictating the spatiotemporal migration of seismicity. We conclude that near‐real‐time monitoring of seismic attenuation may help detect fluid redistribution in active fault systems and may provide useful information for time‐dependent seismic hazard assessment.

central Apennines↗
USGS2000

Pesticides in surface waters of the Upper Colorado River Basin, Colorado, 1996-98

Forty-four river, stream, and drain sites in the Upper Colorado River Basin in Colorado were sampled during 1996-98 to determine the occurrence and distribution of pesticides in the basin. In a fixed-station study, 57 surface-water samples were collected from October 1996 through January 1998 at four sites. Each site was sampled approximately monthly for up to a year, with more frequent sampling during the spring and summer growing season. In a synoptic study, surface-water samples were collected at 43 sites in the agricultural areas of the Grand Valley and the Uncompahgre River Valley in May 1998. Each site was sampled once. Pesticide concentrations generally were low and varied seasonally and across land uses. Thirty-five pesticides were detected at least once in the two studies, and 11 pesticides accounted for almost 82 percent of the detections. Herbicides were more commonly detected than insecticides, and only the herbicides alachlor and atrazine were detected in more than 50 percent of the samples. Carbofuran was the most commonly detected insecticide and was detected in 19 percent of the samples. Pesticide detections increased and were measured at higher concentrations in the summer months and at the agriculture sites. All pesticide concentrations were less than drinking-water standards, and most complied with human-health advisories and criteria for the protection of freshwater aquatic life. Pesticides in the Upper Colorado River Basin in Colorado were detected slightly less frequently and generally at lower concentrations than in 20 National Water-Quality Assessment Program study-unit basins that collected water-quality data from 1992 through 1996. Results from surface-water sampling conducted during 1996?98 in the Upper Colorado River Basin in Colorado indicate that beneficial uses of water were not being impaired by the presence of pesticides in surface waters in the basin.

Colorado↗
USGS2024

Characteristics of the fault damage zone From high-resolution seismic imaging along the Palos Verdes Fault, California

The distribution and intensity of fault damage zones provides insight into fault activity and its relationship to fluid flow in the crust. Presently, measures of the in-situ distribution of fault damage remain limited and along-strike studies are rare. This study focuses on an offshore section Palos Verdes Fault damage zone that spans 28 km, near Los Angeles, California. To investigate the previously unresolved shallow (∼400 m below the seafloor) fault damage zone we use densely spaced (∼500 m line separation) newly collected sparker multichannel seismic lines and sub-bottom profiles. The combination of high-resolution acquisition methods and specialized seismic processing workflows provide improved imaging of shallow faulting. We apply a multi-trace similarity technique to identify discontinuities in the seismic data that may be attributed to faults and fractures. This fault detection approach reveals diverse fault damage patterns on adjacent seismic profiles. However, a discernible damage zone pattern emerges by stacking multiple damage detection profiles along strike. We find that peak damage identified in this way corresponds to the active main fault strand, confirmed in this study, and thus the technique may be useful for identifying active fault strands elsewhere. Additionally, we observe that the variable width of the damage zone along strike is controlled by fault obliquity. Furthermore, our observations reveal a correlation between fault damage and seafloor fluid seeps visible in the water column, suggesting that damage plays a role in controlling fluid flow around the fault.

California↗
USGS2025

Near-fault amplification and ground motion variability during the 2019 Ridgecrest, California sequence

We estimate ground-motion variability near the 2019 M 7.1 Ridgecrest earthquake sequence. Accurate seismic hazard estimation requires understanding ground-motion spatial correlations, yet many studies lack the dense station coverage needed to resolve small-scale variability. The 2019 M 7.1 Ridgecrest earthquake sequence presents a unique opportunity to examine ground motions and their spatial correlations at a range of interstation distances. The permanent seismic network was augmented with hundreds of temporary stations including several fault-crossing nodal arrays. We compute the event ( δE i ) and within-event ( δW ij ) residuals from the observed peak ground velocity and peak ground acceleration data to isolate potential sources of ground-motion variability. We then compare δW ij between station pairs that record an event to understand the semivariance of the ground motion versus interstation distance. By fitting an exponential model to the semivariances, we determine a correlation range of 25 km for the Ridgecrest region. Although the exponential model fits the broad-scale increase of semivariance with interstation distance, we also observe smaller-scale trends. We find that ground motions are less correlated for station pairs that are near or across faults that ruptured during the 2019 Ridgecrest sequence. We also find large, positive median δW ij with relative values 2–3 times larger than nearby stations for individual stations’ near-fault traces. Near-fault amplification and greater ground-motion variability can delineate fault zones and may locally increase the seismic hazard.

California↗
USGS2026

Core microbiomes as a potential fingerprinting method of Western USA dust sources

Introduction: Changing frequency and intensity of dust emissions impacts ecosystems and human health. Dust carries microbes, nutrients, heavy metals, and other materials that may change environmental biogeochemistry at deposition sites. Identifying dust sources provides key information on where and when mitigation strategies should be employed. However, commonly used geochemical or isotopic tracers are often not capable of distinguishing between geographic regions. Methods: We explored whether soil bacterial communities may provide distinct fingerprints of dust sources in the western United States. We identified bacterial core communities of dust from ten locations monitored by the National Wind Erosion Research Network (NWERN) with varied land use (cropland, rangeland, and playa), and compared communities to location, soil, and regional characteristics. Samples were collected monthly from Modified Wilson and Cooke (MWAC) samplers, composited by season (spring, summer, and fall), and analyzed using 16S rRNA sequencing. Results: We found distinct bacterial core communities that reflected dust source characteristics. In order of importance, precipitation levels ( p = 0.0001), location ( p = 0.0001), soil texture ( p = 0.0001), seasonality ( p = 0.0001), and elevation (p = 0.0002) were correlated with bacterial community composition. Discussion: Distinct bacterial core communities were associated with site characteristics such as biocrusts, playas, and military base proximity. Our results suggest that the use of core microbiomes may offer a fingerprinting method to identify dust source regions.

Colorado, Nevada, New Mexico, North Dakota, Oklaho↗
USGS2025

Widespread anhydrite saturation in Laramide-age arc magmas of southwestern USA

Anhydrite is considered a rare mineral phase in magmas, with only ∼33 documented occurrences worldwide. However, anhydrite readily decomposes in the near-surface environment, making it difficult to recognize its former presence in rocks collected at or near Earth’s surface. In such samples, only small anhydrite inclusions fully shielded within other minerals can have survived. During a recent field trip to the southwestern USA, we sampled 17 Laramide-age (ca. 40−80 Ma) magma systems, most of which are associated with porphyry copper deposits. A systematic search for anhydrite inclusions preserved within apatite, amphibole, plagioclase, and quartz phenocrysts in ∼100 rock samples by optical microscopy and Raman spectroscopy revealed that each of these 17 magma systems was at least temporarily anhydrite-saturated. Also, most previously identified magmatic anhydrite-bearing intrusions are associated with porphyry copper deposits, and both intrusive and volcanic rocks containing magmatic anhydrite show high Sr/Y ratios. These observations suggest that anhydrite saturation and porphyry copper formation are linked via magma fractionation at high pressure. Compared to average arc magmas, anhydrite-bearing magmas are unusually oxidized and sulfur-rich and seem to also be unusually water-rich. Hence, our preferred interpretation is that magma generation and/or fractionation at high pressure promotes the formation of superhydrous and oxidized magmas, which in turn promotes high sulfur contents and ultimately the precipitation of anhydrite. The high mineralization potential of these magmas does not need to result from their high sulfur content but could be caused by other properties of high-pressure magmas.

Arizona, New Mexico↗
USGS2026

Rapid characterization of the 2025 Mw 8.8 Kamchatka, Russia earthquake

The 29 July 2025 M w 8.8 Kamchatka, Russia, earthquake was the sixth largest instrumentally recorded earthquake. This event was seismically well observed at regional and teleseismic distances, but publicly available near‐source data were sparse at the time of the event, presenting unique challenges for rapid source and impact characterization. The U.S. Geological Survey (USGS) National Earthquake Information Center provides global real‐time monitoring for earthquakes, including rapid response information products that estimate source characteristics, shaking, and the resulting impacts. We describe the USGS rapid response earthquake information products following the Kamchatka event and discuss their implications for ongoing hazards in the region. We describe potential improvements to our response workflows motivated by this event, including more rapid constraints on source geometries and the automated selection of fault geometries for finite‐fault inversions. The rapid response products together support the interpretation of a unilateral southwestward rupture with significant slip on the southwestern end of the rupture extent. The M w 8.8–9.0 event in 1952, which ruptured a comparable extent of the Kuril–Kamchatka subduction interface, has many similarities to the 2025 rupture. This illustrates that slip deficits may remain following great earthquakes and highlights the usefulness of comparative studies between historic and modern events.

Kuril-Kamchatka subduction zone↗
USGS1976

Availability of ground water in the area surrounding the Trident submarine construction facility, Kitsap County, Washington

This report provides general information on the water resources--with emphasis on ground-water occurrence and availability--in that part of Kitsap County (referred to as Trident Impact Area) that would be most affected by the development of the Trident submarine-construction facility at Bangor. Ground water occurs principally in a lower, confined aquifer near and below sea level and an upper, unconfined aquifer above sea level; both aquifers are extensive beneath the study area and are composed of sand and gravel layers separated by poorly permeable silt, clay, and till layers. Locally throughout the area efficiently constructed wells tapping each of the two aquifers are capable of yielding 100 to 1,300 or more gallons per minute. The estimated 1970 water use in the study area averaged about 13 million gal/day (gallons per day); of this amount about 9 million gal/day came from surface-water sources--from a large reservoir outside the study area--and about 4 million gal/day came from ground water pumped from the two aquifers. Anticipated water use soon will be about 18 to 21 million gal/day; virtually all the additional quantity required (about 5 to 8 million gal/day) above present use must come from ground-water sources. Preliminary evaluation of the aquifers suggests that an additional 1.5 million gal/day can be developed from the upper aquifer and 7 million gal/day from the lower aquifer. Existing wells tapping the lower aquifer might yield additional water and increase the total yield in the area by 3.5 million gal/day, and new wells drilled in selected areas could produce an additional 3.5 million gal/day from this aquifer. However, additional, large-scale ground-water withdrawal from the lower aquifer could induce saltwater intrusion into wells situated in coastal areas. Recommended future studies include additional data collection, monitoring of chloride concentrations in water from large-yield coastal wells, confirmation of well yields and capacities, and test-well drilling and analysis.

Washington↗
USGS2025

Landsliding follows signatures of wildfire history and vegetative regrowth in a steep coastal shrubland

Five years after the deadly and destructive 9 January 2018 Montecito debris flows (Santa Barbara County, California, USA), an atmospheric river storm on 9 January 2023 triggered widespread landsliding that affected many of the same drainages in the Santa Ynez Mountains. Using high-resolution aerial imagery, we identified >10,000 landslides over an ∼160 km 2 area. Most of the landslides were shallow (<1 m in depth) translational debris slides that initiated on steep (∼40°), south-facing hillslopes, with the highest incidence of landsliding in a sandstone-dominated bedrock unit. The landslides mobilized into debris flows and delivered substantial quantities of sediment downstream, which contributed to costly infrastructure impairments. We detected order-of-magnitude differences in landslide density across the study area that could not be attributed to variations in geomorphology (topographic aspect and slope), geology (bedrock type), or hydrology (seasonal antecedent rainfall, peak hourly storm rainfall intensity, total storm rainfall), which are usually considered relevant factors for shallow landsliding. Rather, we found that vegetation regrowth following wildfire was likely a relevant factor associated with the highly variable landslide densities. Hillslopes with less and different types of vegetation regrowth after fire appear to have been more susceptible to shallow landslides. We identify a possible vegetation control on postfire landsliding, which highlights an opportunity for hypothesis testing using more advanced techniques to track the evolution of vegetation cover and vegetation type in steep shrubland environments following wildfire.

California↗
USGS2025

Site response and wave propagation effects in the eastern United States

Fourier amplitude spectra from regional earthquakes in the eastern United States are used in a parametric inversion for source, path, and site effects. Five earthquakes are selected for analysis during the installation of the United States National Seismic Network (US), Earthscope’s USArray Transportable Array (TA), and other temporary arrays to maximize station coverage. A global search algorithm is used to solve for site response from 0.1 to 15 Hz, corner frequency, geometrical spreading ( r - γ ), and frequency dependent anelastic attenuation in the form Q(f) = Q o f α . Tradeoff between moment and geometric spreading is handled by fixing the moment. The tradeoff between corner frequency and Q(f) is solved by selecting the value of corner frequency that minimizes an objective function defined over all stations. Values of site response and attenuation parameters show a strong spatial correlation with the physiographic provinces of the eastern United States. Site response for the Atlantic Coastal Plain is consistent with previous work using spectral ratios relative to a reference site, defined by strong resonance peaks correlated with the thickness of sediments. Site response for the other physiographic provinces is markedly different from the coastal plain, with a lack of distinct resonance peaks and a broad moderate high at frequences from 0.1 to 0.5 Hz consistent with the hard-rock geology of the regions. Like site response, Q(f) has a strong correlation with physiographic province, showing lower values on the coastal plain and higher values inland. Geometric spreading exponent, γ, decreases with increasing hypocenter distance from just above 1 at a few tens of kilometers to 0.9 at 500 km. The limited range in geometric spreading values is attributed to starting the Fourier transform window at the S ‐wave arrival for all distances and averaging over multiple wave types.

eastern United States↗
USGS2026

Development of liquefaction-targeted design basis in U.S. seismic provisions

Seismic design criteria for new buildings in the United States have historically focused on life safety performance objectives through minimizing the potential for structural collapse. Development of design criteria to meet this performance objective has evolved over time, leading to the current, risk-targeted maximum considered earthquake (MCE R ) design basis. Corresponding MCE R seismic design loads incorporate the full range of probabilistic ground motion hazard information available from U.S. Geological Survey National Seismic Hazard Models, as well as a representation of the uncertainty in structural collapse capacity. Despite these considerable advances in structural design criteria, guidelines for geotechnical hazards such as liquefaction-induced damage, which continue to be based on uniform ground-shaking hazard (i.e., MCE ground motions), are not well-connected with risk-based, collapse prevention performance objectives. This study, undertaken as part of Building Seismic Safety Council efforts to update recommended seismic provisions through the National Earthquake Hazards Reduction Program, presents probabilistic, liquefaction-targeted design objectives and corresponding ground motion parameters, with a focus on improving consistency in first-level screening criteria for liquefaction hazard assessment. This study identifies a potential maximum acceptable annualized probability of liquefaction triggering ( P L,max ) of 2.5% in 50 years, above which liquefaction consequences and potential mitigation measures would be considered. The P L,max -based design basis is calibrated to effective probabilistic liquefaction hazard levels obtained in practice using current MCE peak ground accelerations, and could be implemented via a new liquefaction-targeted maximum considered earthquake (MCE L ) peak ground acceleration ( PGA L ). The proposed PGA L is derived from probabilistic liquefaction hazard curves and includes uncertainties inherent to liquefaction hazard modeling. The potential design effects of the proposed PGA L are demonstrated via comparison with current MCE PGA levels, effects on first-level liquefaction screening analyses, and improvements in the consistency of P L,max levels across the United States, and for a wide range of site conditions.

conterminous United States↗
USGS2008

Geology and geochemistry of Jurassic plutonic rocks, Baboquivari Mountains, south-central Arizona

Among the plutons of the Jurassic magmatic arc segment in southern Arizona, northern Sonora, and southern California, two distinctive rock types dominate, occurring together in numerous mountain ranges: porphyritic, titanite-bearing hornblende-biotite granodiorite; and biotite leucogranite. These characteristic rocks are particularly well preserved and well exposed around Kitt Peak, in the Baboquivari Mountains, southern Arizona. The Middle to early Late Jurassic Kitt Peak Plutonic Suite (KPPS) comprises three units: Aguirre Peak Quartz Diorite (APQD;<170 Ma, >=165 Ma, U-Pb), chiefly hornblende mesodiorite; Kitt Peak Granodiorite (KPGD; 165 Ma), hornblende-biotite granodiorite and monzogranite; and Pavo Kug Granite (PKG; 159 Ma), equigranular biotite leucogranite. This range also hosts a fourth, regionally unusual, Jurassic plutonic unit: the Baboquivari Peak Perthite Granite, related to post-arc, Late Jurassic crustal extension. Binary variation diagrams, incompatible element patterns, and REE spectra indicate APQD and KPGD are consanguineous. APQD basaltic mesodiorite evidently approximates the parental magma for the APQD–KPGD series. More evolved rocks of this series evidently formed by combinations of fractional crystallization and mixing with or assimilation of crustal material. PKG is less directly related to APQD and KPGD than these two are to one another.

Arizona, Sonora↗
USGS2026

Numerical modeling of Late Pleistocene to Holocene earthquake-induced progressive rock slope damage in the paraglacial Serpentine valley, Prince William Sound, Alaska

Landslides in deglaciating fjords pose a potential tsunami threat to nearby communities; however, processes contributing to long-term progressive rock damage and landslide conditioning remain poorly constrained in paraglacial settings. Here, we analyse the role of earthquake-induced rock mass damage over late Pleistocene to Holocene time scales, as a conditioning factor for modern landslides, using distinct element numerical modelling to assess spatial and temporal patterns of fracture propagation influenced by varying glacier thickness. Conceptualised numerical models were parameterised by in situ rock mass, glacial and topographic conditions in Serpentine valley, located in Prince William Sound, Alaska, where several large landslides are actively developing along the western valley wall. Results show that although rigid glacier buttressing reduces co-seismic rock mass damage, it does not suppress it completely, and damage occurs both above and below the glacier surface elevation. We demonstrate that topography, and especially steep slopes with topographic convexity, as well as preexisting damage of joint networks and faults inherited from tectonic and exhumation induced stresses, exert primary control on the location of new co-seismic damage. Simulations representing a simplified deglaciation sequence over the past 25 ka, with a series of 10 evenly spaced earthquakes, generated rock mass damage patterns that qualitatively match in situ landslide structural and kinematic observations at one instability in Serpentine valley. Our conceptual study helps clarify the role of repeated seismicity over glacial timescales as a long-term conditioning process for paraglacial rock slope failure, highlighting spatial and temporal patterns of progressive damage accumulation, with outcomes relevant for modern landslide hazard assessment.

Alaska↗
USGS2026

Scenarios to assess the future water availability in the Mississippi River Valley Alluvial Aquifer for the Cache River and Grand Prairie Regions of Arkansas

The U.S. Geological Survey, as part of the Arkansas Groundwater Initiative, developed forecast scenarios using previously calibrated MODFLOW 6 groundwater models that focused on the Cache and Grand Prairie Critical Groundwater Areas to assess the impact of future climate and water management strategies on the Mississippi River Valley alluvial aquifer. A Soil Water Balance model was used to forecast recharge and irrigation water use. The forecast scenario period was from January 1, 2019, through December 31, 2055, with monthly stress periods. Twenty scenarios were simulated and included seven alternate climate forecasts, five 13 general groundwater pumping reduction scenarios (round 1), and groundwater pumping reduction scenarios by crop type and for the Bayou Meto Water Management Project and Grand Prairie Area Demonstration Project (round 2). Declines in saturated thickness within the Cache Critical Groundwater Area were larger for 18 of the 20 scenarios as compared to outside of the Critical Groundwater Area. The largest average increase in saturated thickness inside the Critical Groundwater Area was 6.4 m which occurred for the round 1, 50 percent reduction scenario. Automatic reductions in groundwater pumping by MODFLOW 6 in the Cache simulation ranged from 0.02 to 13.1 percent of total groundwater pumping. For the Grand Prairie model domain, the average change in saturated thickness of the Mississippi River Valley alluvial aquifer inside the Critical Groundwater Area for the forecast period ranged between -6.6 to 1.7 m. The average saturated thickness of the Mississippi River Valley alluvial aquifer inside the Grand Prairie Critical Groundwater Area declined for 16 of the 20 scenarios. The average reduction in requested groundwater pumping for all scenarios inside the Grand Prairie Critical Groundwater Area was 25.1 percent, and the largest reduction was 46.5 percent.

ESS Open Archive↗
USGS2026

Determining Volcanic Risk in Auckland (DEVORA) Research Programme—A transdisciplinary approach to address the challenge of distributed volcanism in an urban environment

The Determining Volcanic Risk in Auckland (DEVORA) Research Programme was launched in 2008 to address the challenges associated with monogenetic volcanism in an urban setting and to enhance volcanic risk management in Tāmaki Makaurau Auckland in Aotearoa New Zealand. It is a multi-agency, increasingly transdisciplinary (defined here as research that transcends traditional disciplinary boundaries by integrating diverse types of knowledge, perspectives, and methods from academic and non-academic participants to create novel solutions to complex problems), and collaborative research program jointly led by Waipapa Taumata Rau University of Auckland and Earth Sciences New Zealand (ESNZ; formerly GNS Science), with core funding from Natural Hazards Commission Toka Tū Ake (NHC; formerly the Earthquake Commission, EQC) and Te Kaunihera o Tāmaki Makaurau Auckland Council (AC). The primary research focus of DEVORA is to investigate the geologic history, volcanic hazards, and risk posed by the basaltic intraplate Auckland Volcanic Field. Disruption from ash fall and gas from other Aotearoa New Zealand volcanoes is also considered. DEVORA’s work to explore exposure and vulnerability in Tāmaki Makaurau Auckland is also useful for assessing risks from other non-volcanic natural hazards, such as seismic and tsunami hazards. The greater Tāmaki Makaurau Auckland region has an ethnically and socio-economically diverse population of approximately 1.7 million, representing about one-third of the Aotearoa New Zealand population, and hosts critical infrastructure of national significance. The size and nature of the populace, consequential economic base, and important infrastructure within Tāmaki Makaurau Auckland mean that the effects of a volcanic eruption would be felt nationally, including through the disruption of air travel to Aotearoa New Zealand. The hazards from such an eruption could potentially affect hundreds of thousands of people, businesses, and lifelines (critical infrastructure). A considerable challenge for emergency and risk managers is the monogenetic nature of the volcanic field. It is not known where or when the next eruption will occur, how much warning we may get before an eruption, nor how an eruption and its effects might unfold. In this contribution, we highlight the concept and collaborative intent of the DEVORA Programme and show how it has evolved over the 16 years since its inception. We describe how DEVORA has unified more than 100 researchers (including more than 50 graduate students) and numerous stakeholders to address key issues facing Tāmaki Makaurau Auckland and describe how research findings are being implemented into policy and communicated to stakeholder agencies and the public. We also illustrate the broader influence of the DEVORA Programme and provide some learnings that might benefit others embarking on similar integrated projects, especially those focused on distributed volcanism in and near populated areas.

Auckland↗
USGS1998

Technique for simulating peak-flow hydrographs in Maryland

The efficient design and management of many bridges, culverts, embankments, and flood-protection structures may require the estimation of time-of-inundation and (or) storage of floodwater relating to such structures. These estimates can be made on the basis of information derived from the peak-flow hydrograph. Average peak-flow hydrographs corresponding to a peak discharge of specific recurrence interval can be simulated for drainage basins having drainage areas less than 500 square miles in Maryland, using a direct technique of known accuracy. The technique uses dimensionless hydrographs in conjunction with estimates of basin lagtime and instantaneous peak flow. Ordinary least-squares regression analysis was used to develop an equation for estimating basin lagtime in Maryland. Drainage area, main channel slope, forest cover, and impervious area were determined to be the significant explanatory variables necessary to estimate average basin lagtime at the 95-percent confidence interval. Qualitative variables included in the equation adequately correct for geographic bias across the State. The average standard error of prediction associated with the equation is approximated as plus or minus (+/-) 37.6 percent. Volume correction factors may be applied to the basin lagtime on the basis of a comparison between actual and estimated hydrograph volumes prior to hydrograph simulation. Three dimensionless hydrographs were developed and tested using data collected during 278 significant rainfall-runoff events at 81 stream-gaging stations distributed throughout Maryland and Delaware. The data represent a range of drainage area sizes and basin conditions. The technique was verified by applying it to the simulation of 20 peak-flow events and comparing actual and simulated hydrograph widths at 50 and 75 percent of the observed peak-flow levels. The events chosen are considered extreme in that the average recurrence interval of the selected peak flows is 130 years. The average standard errors of prediction were +/- 61 and +/- 56 percent at the 50 and 75 percent of peak-flow hydrograph widths, respectively.

Maryland↗