USGS ScienceSearch

SEARCH · USGS Science

Results for “Permafrost and Periglacial Processes”

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

Conceptual and numerical groundwater flow model of the Iowa River alluvial aquifer near Tama County, Iowa, 1980 through 2022

The Iowa River alluvial aquifer is an important source of water on the Meskwaki Settlement in Tama County, Iowa, which is land owned by the Sac & Fox Tribe of the Mississippi in Iowa (commonly known as the Meskwaki Nation). The U.S. Geological Survey constructed a groundwater flow model, including a conceptual and numerical model, of the Iowa River alluvial aquifer and underlying hydrogeologic units near the Meskwaki Settlement in Tama County, Iowa, for the period of January 1980–August 2022 to estimate the fraction of water pumped from the Iowa River alluvial aquifer by Meskwaki Settlement wells that is derived from streamflow depletion in the Iowa River and its tributaries. Streamflow depletion is a reduction in streamflow caused by groundwater pumping and includes the interception by groundwater production wells of water that otherwise would have been discharged to streams (called “captured groundwater discharge”) and induced infiltration of streamflow to the production wells. Calibrated model runs were performed with no simulated pumping and simulated pumping only at Meskwaki Settlement wells, and the change in simulated flow rates between the groundwater system and streams for the two model runs represents the amount of streamflow depletion in the Iowa River and tributary streams resulting from pumping at the Meskwaki Settlement wells. Streamflow depletion in the Iowa River and its tributaries as a percentage of simulated pumping at the Meskwaki Settlement wells was calculated by dividing this difference by the total simulated pumping rate for the Meskwaki Settlement wells. The model results demonstrate that the mean monthly streamflow depletion, including induced infiltration and captured discharge, in the Iowa River and its tributary streams as a percentage of mean monthly pumping at the Meskwaki Settlement wells was 97.4 percent and ranged from 65.4 to 112 percent. Of the total streamflow depletion, mean monthly induced recharge was 20.9 percent and ranged from 4.9 to 37.2 percent. Mean monthly captured discharge was 76.5 percent and ranged from 57.1 to 97.1 percent. These results indicate that most of the water pumped from the Meskwaki Settlement wells is the result of streamflow depletion, in the form of both induced infiltration and captured discharge.

Iowa

Mechanics and statistics of postseismic shaking

Analysis of two weeks of continuous post-seismic shaking after the 2019 M7.1 Ridgecrest, CA earthquake sequence using 4 nearby borehole seismometers reveals that continuous ground motions decay as Omori’s law in time and follow the Gutenberg-Richter distribution in logarithmic amplitude. The measured temporal decay in amplitudes agrees with predictions of the rate-and-state framework and indicates shaking amplitudes are proportional to the velocity of afterslip. Our ground motion-based statistical framework provides a basis to forecast shaking intensity in the minutes to hours after a large earthquake.

California

Near-real-time earthquake-induced fatality estimation using crowdsourced data and few-shot large-language models

When a damaging earthquake occurs, immediate information about casualties (e.g., fatalities and injuries) is critical for time-sensitive decision-making by emergency response and aid agencies in the first hours and days. Systems such as the Prompt Assessment of Global Earthquakes for Response (PAGER) by the U.S. Geological Survey (USGS) were developed to provide a forecast of such impacts within about 30 min of any significant earthquake globally. However, existing disaster-induced human loss estimation systems often rely on early casualty reports manually retrieved from global traditional media, which are labor-intensive, time-consuming, and have significant time latencies. Recent approaches use keyword matching and topic modeling to identify human casualty-relevant information from social media but tend to be error-prone when dealing with complex semantics in multi-lingual text data and parsing dynamically changing and conflicting human death and injury numbers shared by various unvetted sources in social media platforms. In this work, we introduce an end-to-end framework to significantly improve the timeliness and accuracy of global earthquake-induced human loss forecasting using multi-lingual, crowdsourced social media. Our framework integrates (i) a hierarchical casualty extraction model built upon large language models, prompt design, and few-shot learning to retrieve quantitative human loss claims from social media, (ii) a physical constraint-aware, dynamic-truth discovery model that discovers the truthful human loss from massive noisy and potentially conflicting human loss claims, and (iii) a Bayesian updating loss projection model that dynamically updates the final loss estimation using discovered truths. We test the framework in real-time on a series of global earthquake events in 2021 and 2022 and show that our framework effectively automates the retrieval of casualty information faster but with comparable accuracy to those now retrieved manually by the USGS. The code associated with this work is made available at: https://github.com/SusuXu-s-Lab/Hierarchical-Earthquake-Casualty-Information-Retrieval

International Journal of Disaster Risk Reduction

Monitoring marine eruptions

Introduction Submarine volcanoes produce much of the same seismicity and eruptive activity as subaerial volcanoes and can pose hazards to society. Although they can be monitored with similar techniques and methods as described in other chapters of this volume, their submerged location brings unique challenges. This chapter addresses these challenges and provides recommendations for monitoring volcanoes fully or partly in marine environments to meet the capabilities described in other chapters of this volume. The United States and its territories host dozens of submarine volcanoes with most (around 60) in the Commonwealth of the Northern Mariana Islands. Approximately 20 of the Northern Mariana Islands submarine volcanoes are known to be hydrothermally active, and 10 have confirmed eruptions since the 1950s (for example, Baker and others, 2008; Tepp and others, 2019a). Nine of those volcanoes were considered by the National Volcanic Threat Assessment (Ewert and others, 2018) to have a combination of eruptive type and summit depth that poses a higher risk of hazardous eruptions, although only one was listed as a moderate (level 3) threat. Other notable submarine volcanoes of interest to the United States that have historically erupted are Axial Seamount off the Washington State coast, Kamaʻehuakanaloa in Hawaiʻi, and Vailuluʻu seamount in American Samoa. All of these, however, have a low risk of hazards because of their depth (greater than 600 meters below sea level) and eruptive type and so are not included in the National Volcanic Threat Assessment. In addition to submarine volcanoes, the submerged flanks of island volcanoes can also be a source of hazardous submarine eruptions—for example, the 1877 eruption of Mauna Loa, Hawai‘i, in Kealakekua Bay (Wanless and others, 2006). The most notable submarine eruption in recent times was the 2022 eruption of Hunga Tonga–Hunga Haʻapai in Tonga, which was one of the largest eruptions on Earth in the past 100 years. It created a massive volcanic plume, unprecedented shock waves, and far-reaching tsunami (Lynett and others, 2022). Other recent submarine eruptions in the Pacific Ocean Basin have produced subaerial plumes that reached aircraft heights (Carey and others, 2014) and large pumice rafts that can affect marine traffic and harbors (for example, Jutzeler and others, 2014; Kornei, 2019). These examples illustrate the potential hazards of major submarine eruptions. Yet, submarine volcanoes are largely unmonitored, and many eruptions occur that are unnoticed or only identified hours or days afterward. Within U.S. territory, submarine volcanoes in the Northern Mariana Islands have been known to produce eruptive activity that can affect society. Reports from fishermen and other marine vessels in the Northern Mariana Islands have noted underwater explosions, sea-surface discoloration, and bubbling water, all of which are known to be signs of submarine volcanic activity. South Sarigan seamount, located about 160 kilometers (km) north of Saipan, erupted in 2010 from greater than 150 meters below the sea surface, resulting in a gas and ash plume that reached more than 11.9 km into the atmosphere (for example, Searcy, 2013; Embley and others, 2014), high enough to affect international air traffic. Precursory and co-eruptive seismicity was detected on the regional Northern Mariana Islands seismic network (Searcy, 2013) and on global monitoring instruments (Green and others, 2013). Monitoring of submarine volcanoes is best accomplished with marine-based instrumentation, which is also useful for monitoring small island volcanoes that may not have the land area necessary for comprehensive subaerial monitoring. The primary marine-based instrumentation used for submarine volcanoes includes ocean-bottom pressure sensors to assess sea-floor deformation, ocean-bottom seismometers (OBSs) to detect seismicity, and both moored and ocean-bottom hydrophones to detect submarine explosions. Other sensors offer important monitoring data, such as turbidity, temperature, and chemistry of hydrothermal emissions. Marine-based instruments are typically deployed in campaign-style networks with no real-time telemetry owing to cost considerations and technical limitations. However, when necessary, marine instruments can be operated in real time using cables to transmit data to land-based facilities; other technologies for this purpose are in use or in development, such as acoustic transmission from the instrument to a moored buoy (Matsumoto and others, 2016) and a winch-based system with a satellite antenna that is part of the instrument mooring (Matsumoto and others, 2019). Emerging technologies for marine-based monitoring may be considered as part of a future monitoring plan. These technologies include ocean gliders and floats with on-board hydrophones that have been used to record earthquakes and submarine eruptions (for example, Matsumoto and others, 2013; Sukhovich and others, 2015) and fiber-optic cables that have been used as strainmeters to detect earthquakes (for example, Marra and others, 2018; Lindsey and others, 2019). Land-based instruments and satellites can also provide some capability for monitoring submarine volcanoes, but they provide more limited observations than marine-based instrumentation.

Scientific Investigations Report

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

Geologic map of the Guinevere Planitia quadrangle (V–30), Venus

This 1:5,000,000-scale geologic map of the Guinevere Planitia quadrangle divides the region into 15 geologic material units, defined using Magellan synthetic aperture radar (12.6-centimeter-wavelength radar system; 75 meters per pixel) datasets and including upland terrain units (2.4 percent of the surface area), plains materials units (59 percent), flow materials associated with named and unnamed eruptive centers (37.2 percent), small volcanic edifices, and impact crater materials (1.4 percent). Upland terrain units consist of tessera and lineated upland material, plains materials consist of Guinevere regional plains and Guinevere lineated and mottled plains, and flow materials consist of lobate flow material and plains-forming flow material. Specific lobate flows associated with Atanua Mons, Tuli Mons, Var Mons, and Uilata Fluctus are mapped separately. Other mapped units are impact crater material and small volcanic edifice. In addition to geologic units, we mapped linear features that show patterns of deformation or flow across the quadrangle. These consist of faults, wrinkle ridges, broad arches, channels, troughs, and flow direction indicators. The map region also contains several small volcanic features: shields, depressions, and craters. These, in combination with the plains, large volcanoes, and coronae, show the pervasive influence of volcanism across Venusian lowlands. The rims of nine identified impact features are delineated; large bright and dark haloes, which in some cases are associated with individual impact craters, are mapped as surficial mantling deposits. We documented spatial relationships using the stratigraphic and cross-cutting relationships of the quadrangle’s geologic units and features to provide a synthesis of the region’s geologic history. The upland terrain of the quadrangle indicates intense tectonic deformation and uplift. It is exposed as embayed remnants, typically within the plains, and represents the oldest geologic materials locally and across the region. Guinevere plains and the plains-forming flow unit appear to be assemblages of volcanic flows from multiple sources, including distinct coronae and corona-like structures. The temporal evolution of Guinevere lineated and mottled plains was likely protracted, with continued formation of small volcanic edifices over a long period. The morphologic and radar brightness characteristics of volcanoes in the region indicate their growth may have involved (1) multiple large-scale eruptive centers with recognizable spatial and temporal sequences, (2) extensive lava flow fields with a multitude of flows producing complex, overlapping patterns, and (3) numerous small volcanic edifices, including shields, domes, and cones. Although geologic patterns common to other regions of Venus are evident in the Guinevere Planitia quadrangle, local relative age relationships are inconsistent or unclear, preventing robust stratigraphic correlation. The mapping results do, however, indicate complicated local sequences of volcanic and tectonic activity.

Scientific Investigations Map

The 3D National Topography Model Call for Action—Part 1. The 3D Hydrography Program

The U.S. Geological Survey is initiating the 3D Hydrography Program (3DHP), the first systematic remapping of the Nation’s surface waters since the original 1:24,000-scale topographic mapping program was active from 1947 to 1992. Building on decades of experience maintaining the National Hydrography Dataset (NHD), the Watershed Boundary Dataset (WBD), and the NHDPlus High Resolution (NHDPlus HR), the 3DHP will completely refresh the Nation’s hydrography data and improve discovery and sharing of water-related data. The design of the 3DHP is based on the results of a study that estimated that the fully implemented program would have the potential to provide more than $1 billion in benefits to Federal, State, Tribal, Territorial, and local governments and to private and nonprofit organizations every year, in addition to myriad societal benefits. The 3DHP would directly support better decision making regarding water resources by providing more accurate, complete, and integrated information than is currently available. The 3DHP datasets will include a three-dimensional (3D) hydrography network generated from and integrated with elevation data from the 3D Elevation Program (3DEP) to better represent stream gradients and channel conditions, along with waterbodies, hydrologic units, hydrologically enhanced elevation and other surfaces, and more consistent and accurate attributes. The 3DHP datasets will inherit key attributes of the NHD, WBD, and NHDPlus HR, and they also will include new attributes and links to other data such as the U.S. Fish and Wildlife Service National Wetlands Inventory, groundwater data, and engineered hydrologic systems such as stormwater networks. The 3DHP will be designed to provide a set of open and interoperable web-based tools, maps, and data catalogs, creating a robust system for users to reference their information about water; the system elements are collectively referred to as the “infostructure.” The 3DHP and the infostructure can provide a foundational geospatial underpinning for the Internet of Water, a community-based effort to modernize tools and technologies to share water data. As proposed, the 3DHP would begin providing products and services to the public in 2024.

Circular

“Leaky weirs” capture alluvial deposition and enhance seasonal mountain-front recharge in dryland streams

“Leaky weirs” are rock structures installed in dryland streams, which are anchored into exposed bedrock, loosely cemented, and designed to allow water to slowly pass through. They are being tested at a ranch in southeastern Arizona, USA, to restore and conserve the historic range and desert wetlands. Data are collected to assess how leaky weirs impact surface water, subsurface water, and groundwater recharge—including stream discharge, timing, and depth of infiltration, and groundwater elevations. Three adjacent watersheds, two with outlets just below leaky weirs and one with leaky weirs farther upstream, were instrumented with water-level loggers, wildlife cameras, and crest stage instruments with temperature sensors in the soil. As most groundwater recharge is assumed to be focused along the mountain fronts in this region, mountain-block recharge is also evaluated to differentiate between the two using isotope analyses. Finally, a single, late-season flood event is scrutinized to consider the leaky weir effect on all monitored components in the water budget. Results indicated groundwater flow is primarily from the mountains to the east via older, regional mountain-block recharge. However, the development of shallow alluvial aquifers is supported by the leaky weirs, that slow flows, capture permeable sediments, and allow infiltration, thus enhancing mountain-front recharge. In turn, these new pockets of water help support the restoration of historic wetlands. Sediment accumulates where leaky weirs are installed, reducing flashy peak flows, and resulting in a series of infiltration ponds along the channel that support vegetation during growing seasons and recharge the shallow aquifer during non-growing seasons.

Arizona

Metallurgical coal—Deposits, production, resources, market dynamics, and supply chain risks

Plain Language Summary Metallurgical coal (met coal; consumed to produce coke for steelmaking) must meet specific chemical and physical specifications. In 2023, the conterminous United States produced 66 million short tons (mst) of met coal, consumed 15.85 mst domestically, exported 51.1 mst, and imported 0.7 mst. Most met coal was produced in the Appalachian Basin, but there are also resources that meet the specifications for met coal in the Western United States.

conterminous United States

Subduction zone earthquake catalog separation tool: Implementation in the USGS 2025 Puerto Rico and U.S. Virgin Islands National Seismic Hazard Model

The U.S. Geological Survey (USGS) periodically releases updates to National Seismic Hazard Model (NSHM) for the United States and its territories leveraging current scientific knowledge and methodologies to guide public policy, building codes, and risk assessments regarding potential ground shaking due to earthquakes that may result in infrastructure damage. In subduction zones, there is a need to separate the earthquake catalog into tectonic regimes to create specific seismicity models for which the most appropriate ground‐motion models are then applied. Here, we describe newly developed methods and software, called CatSep, that classifies subduction zone events into three primary tectonic regimes: crustal, interface, and intraslab. This method incorporates information about the location of the earthquake relative to the subducting slab, the depth of the Mohorovičić discontinuity, and the earthquake’s moment tensor. Applying this method is a first step in the NSHM workflow for regions covering U.S. subduction zones. Results using this subduction zone earthquake catalog separation tool for the 2025 Puerto Rico and U.S. Virgin Islands NSHM earthquake catalog are presented and analyzed.

Puerto Rico, U.S. Virgin Islands

Interplay between tectonics and submarine mass transport deposits in Cortes Basin: New high-resolution geophysics in the Outer California Borderland

The Outer California Borderland (OCB) is an active transform plate boundary offshore Southern California, where the relationship between faulting and submarine mass transport deposits (MTDs) remains poorly understood. Onshore paleoseismic data provide high-resolution earthquake records, whereas marine geophysical data capture longer-term histories. Offshore fault systems pose hazards to infrastructure and dense coastal populations, particularly when linked to submarine landslides. We present new high-resolution geophysical data set (cruise SR2303), including bathymetric and CHIRP sub-bottom data integrated with legacy seismic reflection data and chronostratigraphic constraints from ODP Site 1012 to examine Quaternary MTD recurrence and tectonic controls in the Cortes Basin, OCB. Bathymetry shows deformational features, including slide scarps and previously unmapped fault segments with evidence of Holocene activity. CHIRP profiles reveal 10 stacked MTDs in the East Cortes Basin and 8 in the West Cortes Basin, spanning ∼752 ka with an average recurrence of ∼83.6 ± 1 ka. Acoustic imaging shows 7 MTD intervals coinciding with fault offset increments and fault growth suggesting earthquake-triggered mass wasting. A strong association between MTD occurrences and sea-level extremes also supports glacio-eustatic contribution to slope failure. Stratigraphic correlations suggest quasi-synchronous MTDs across the eastern and western areas, likely triggered by larger eathquakes in the Quaternary. Although the identified MTDs occur relatively far from the Southern California coast, they still pose a potential tsunamigenic hazard requiring further assessment. Moreover, if linked to earthquakes along major strike-slip faults, for example, the Ferrelo fault, the MTDs may provide valuable proxies to constrain rupture scenarios and fault connectivity within the understudied OCB.

California

Effect of land cover type on 3D deformation recovery from synthetically deformed high-resolution satellite optical imagery

The limits of detection for earthquake surface deformation in the spatial domain have improved with advances in remote sensing imagery data availability, resolution, and analysis. Sub-pixel correlation and digital elevation model (DEM) differencing from sub-meter, earthquake-spanning satellite optical imagery has enhanced surface rupture mapping and deformation measurements. However, knowledge of measurement accuracy and uncertainty is limited. To address this, we construct orthophotos and digital elevation models (DEMs) from repeat high resolution (∼0.5 m) satellite optical imagery along two sections of the Garlock fault, California with clear fault geomorphology and differing land cover. We deform later sets of DEMs and images with synthetic earthquakes containing both diffuse and discrete horizontal and vertical displacements. Sub-pixel image correlation and DEM differencing demonstrate how vegetation degrades recovered displacement accuracy. In barren land cover, horizontal displacements are detectable to an expected ∼1/10th-pixel size. With shrubs, trees, and grass, detectable displacements increase to >1/2-pixel size, and filtering results by correlation score and using elevation values as input rather than image values improves accuracy. Vertical displacement detection thresholds remain lower in vegetation, at >1-pixel size. Higher slope angles degrade displacement recovery, worsened by vegetation. Diminishing seasonal separation improves accuracy over vegetated regions, though not to the level achieved in barren environments. These results will inform research and operational efforts on the utility of high resolution satellite optical imagery for detecting deformation in varied land cover. Furthermore, they reveal where alternative measurements, such as from LiDAR or radar interferometry, are required to mitigate the effects of vegetation and capture fine-scale crustal deformation.

California

Groundwater flowpath characteristics drive variability in per- and polyfluoroalkyl substances (PFAS) loading across a stream-wetland system

Groundwater-dependent ecosystems in areas with industrial land use are at risk of exposure to a PFAS chemicals. We investigated one such system with several known PFAS source areas, where high and low permeability sediments (glacial) coupled with groundwater-lake and groundwater/surface-water interactions created complex ‘source to seep’ dynamics. Using heat-tracing and chemical methods, numerous preferential groundwater discharge zones were identified and sampled across the upper Quashnet River stream-wetland system in Mashpee, MA, USA, downgradient of Joint Base Cape Cod (JBCC). Surface-water and groundwater samples were analyzed for 40 PFAS compounds between March and October 2022. Samples were collected from groundwater seeps identified as preferential discharge points (PDPs), wells upgradient of the stream-wetland system, contributing flow-through kettle lakes, and along Quashnet River surface-waters. PFAS from sampled waters contained perfluorinated carboxylates (PFPea, PFHxA, PFNA), perfluorinated sulfonates (PFBS, PFPeS, PFHxS, PFOS), fluorotelomer sulfonates (6:2, 8:2 FtS), and perfluoroalkyl sulfonamides (PFOSA). Samples from PDPs and wells had measured PFAS concentrations ranging from non-detect to 4677 ng/L ng/L (mean = 418 ng/L, std. = 709 ng/L), and a range of deuterium excess values (3.2 to 15.9 per mil) indicative of varying degrees of groundwater-lake interaction prior to groundwater flowpath emergence at PDPs. Correlations ( p < 0.01) between deuterium excess, %PFAS precursors, and terminal PFAS compounds highlighted potential precursor transformations associated with lake-groundwater exchange along flowpaths sourcing PDPs. However, some seepages had higher total PFAS concentrations (>1000 ng/L) than upgradient kettle lakes despite showing lake (evaporative) isotopic signatures, indicating the potential for groundwater flowpath convergence at wetland discharge zones and the influence of lakebed PFAS precursor reactions. Results from these synoptic surveys address gaps in the existing PFAS literature by demonstrating the importance of subsurface fate and transport on PFAS compound concentrations and mass loading in preferential groundwater discharge zones.

Massachusetts

Petrogenesis and mineralization potential of the Bradley Peak komatiitic basalts, Wyoming Province

The ca. 2.72 Ga Bradley Peak greenstone terrane in the Wyoming Province contains spinifex- and cumulate-textured komatiitic volcanic rocks that may host Ni-Cu-PGE sulfide mineralization. New whole rock geochemistry classifies these rocks as Al-undepleted komatiitic basalts derived from a parental melt with ~17 wt. % MgO. Isotopic data (Sm-Nd, Re-Os) and REE profiles suggest a plume source and possible assimilation of juvenile footwall volcanic rocks. Calculated sulfur concentrations at sulfide saturation and low-detection PGE geochemistry indicate these magmas did not saturate a sulfide liquid resulting in sulfide mineralization. Although the Bradley Peak komatiitic basalt flows are unlikely to host magmatic sulfide mineralization, they are likely the source of Au for regional orogenic gold deposits.

Wyoming

Exploring the uncertainty of machine learning models and geostatistical mapping of rare earth element potential in Indiana coals, USA

Rare earth elements and yttrium (REEs) have a wide range of applications in high- and low-carbon technologies. The strategic significance of REEs has grown due to their expanding applications in manufacturing industries and the constrained availability of these essential resources. This research explores the applicability of machine learning models and their uncertainty for assessing the REE potential in coal beds using various coal parameters as inputs. The work focuses on developing a predictive model based on geological variables, excluding considerations related to potential shifts in the commodities market. The Indiana Coal Quality Database was used as the data source. The promising and unpromising indicators derived from the outlook coefficient of samples from the database were used as the REE potential indicator for machine learning classification models. The filter-based approach with bootstrap was used to evaluate the importance of the coal parameters and their prediction uncertainties. Four machine learning methods (linear discriminant analysis (LDA), random forest (RF), support vector machine (SVM), and artificial neural networks (ANN), a data balancing and augmentation approach (Synthetic Minority Over-sampling Technique), and bootstrap resampling techniques were used for building the models and evaluating their prediction capabilities under uncertainty. It was determined that the SVM bootstrap model with ten-times balanced and augmented data provided superior results compared with other models. Finally, stochastic spatial maps of the REE potential within the coal basin were generated using sequential indicator simulation. The spatial maps of the REE potential showed that a 29% area of the Indiana section of the Illinois coal basin has economic potential of REEs, with 90% confidence.

Indiana

Estimating the hypothetical endowment of critical minerals and other commodities in porphyry copper mine waste in the Four Corners states, USA

Society is fundamentally dependent upon commodities that are used in end-use products for the aerospace, defense, energy, telecommunication, and transportation sectors, resulting in centuries of mining to supply these commodities and materials. Waste from these mining operations can remain on the landscape indefinitely, but there is a lack of national understanding of the distribution and scale of such waste features. The renewable energy transition will continue to increase demand for critical minerals and will result in increasing volumes of mine waste on the Earth’s surface. Reprocessing mine waste can reduce environmental risks and recover needed commodities to match growing demand for societal growth. Therefore, understanding the approximate abundance of commodities that may be available for recovery within mine waste features can be an important piece of domestic critical mineral supply.

Arizona, Colorado, New Mexico, Utah

Potentiometric surface maps and groundwater-level hydrographs for confined aquifers of the New Jersey Coastal Plain, 2018

The U.S. Geological Survey, in cooperation with the New Jersey Department of Environmental Protection (NJDEP), prepared potentiometric surface maps for 10 confined aquifers of the New Jersey Coastal Plain physiographic province based on water-level measurements collected during late 2018 and early 2019 from 951 wells in New Jersey and parts of Pennsylvania and Delaware. Maps were prepared for the confined Cohansey aquifer, Rio Grande water-bearing zone, Atlantic City 800-foot sand, Piney Point aquifer, Vincentown aquifer, Wenonah-Mount Laurel aquifer, Englishtown aquifer system, and the upper, middle, and lower aquifers of the Potomac-Raritan-Magothy aquifer system. Potentiometric surface maps indicate regional cones of depression in the following aquifers and the counties in which they are centered: Atlantic City 800-foot sand in Atlantic County, the Piney Point aquifer in Cumberland County, the Wenonah-Mount Laurel aquifer and Englishtown aquifer system in Monmouth and Ocean Counties, the Wenonah-Mount Laurel aquifer in Camden and Gloucester County, the upper aquifer of the Potomac-Raritan-Magothy aquifer system in Ocean County, and the upper, middle, and lower aquifers of the Potomac-Raritan-Magothy aquifer system in Camden County. Cones of depression with smaller areal extents were in the confined Cohansey aquifer, the Rio Grande water-bearing zone, the Atlantic City 800-foot sand centered in Cape May County, the Piney Point aquifer centered in Ocean County, the Wenonah-Mount Laurel aquifer in Salem and Burlington Counties, the Englishtown aquifer system in Camden County, the upper aquifer of the Potomac-Raritan-Magothy aquifer system in Monmouth County, and the middle aquifer of the Potomac-Raritan-Magothy aquifer system in Monmouth, Ocean, and Salem Counties. No cone of depression was interpreted in the Vincentown aquifer. Long-term hydrographs are presented for 75 wells spanning each of the 10 confined aquifers, and contain a mix of discrete water-level measurements and daily mean water levels based on continuously recorded 15-minute data. Changes of water levels during 2014–19, as indicated by the hydrographs, were compared with those of previous periods to assess any departures from historical data. During 2014–19, water levels were stable and fluctuated within similar ranges as previous periods in the following aquifers and locations: all wells in the confined Cohansey aquifer, the Rio Grande water-bearing zone, the Vincentown aquifer, the Englishtown aquifer system, the Piney Point aquifer wells in Burlington and Ocean Counties, six of eight wells in the Wenonah-Mount Laurel aquifer, all wells in the upper and lower aquifers of the Potomac-Raritan-Magothy aquifer system outside NJDEP Critical Areas, and all wells in the middle aquifer of the Potomac-Raritan-Magothy aquifer system except those within NJDEP Critical Area II. Increasing water levels in 2014–19, ongoing since historical periods, were indicated in the following aquifers and locations: Atlantic County wells in the Piney Point aquifer, all wells in the upper and lower aquifers of the Potomac-Raritan-Magothy aquifer system outside NJDEP Critical Areas, and all wells in the middle aquifer of the Potomac-Raritan-Magothy aquifer system within NJDEP Critical Area II. Water levels in the Atlantic City 800-foot sand also increased during 2014–19 in wells in Atlantic County and northern Cape May County closer to the center of the cone of depression in that aquifer, which is a response unique to this period and absent from previous periods. During 2013–19, continued decreasing water levels, ongoing since previous periods, were indicated by hydrographs of Atlantic City 800-foot sand wells in southern Cape May County, Piney Point aquifer wells in Cumberland County where the regional cone of depression is located, and two wells in the Wenonah-Mount Laurel aquifer—070478, which in 2014–19 departed from previous periods, and 330020, which continued a gradual decrease throughout its period of record.

New Jersey

Paleomagnetism of the Harrat Rahat volcanic field, Kingdom of Saudi Arabia—Geologic unit correlations and geomagnetic cryptochron identifications

Paleomagnetic rock samples were collected from 173 drill sites in the Quaternary alkali basaltic volcanic field of northern Harrat Rahat, Kingdom of Saudi Arabia. Laboratory measurements on these samples established that lava flows and vent complexes—identified and mapped from field characteristics, rock types, and compositions as products of single or temporally close eruptions—typically record single, or very similar, directions of remanent magnetization. Correlations defined through geologic mapping, spatial association, geochemistry, geochronology, and identical mean remanent directions indicate at least 16 brief episodes of temporally clustered eruptions. These episodes had durations of a few centuries or less. Anomalous remanent magnetic directions were found for at least 13 mapped lavas of northern Harrat Rahat, which demonstrate that they were acquired during brief geomagnetic cryptochrons during the Brunhes Normal Polarity Chron. These uncharacteristic directions enhance the opportunity to identify common eruptive episodes, and to better understand and evaluate assessments of eruption ages based on 40 Ar/ 39 Ar geochronology. Combining paleomagnetic and regional archaeomagnetic results for the youngest eruptions allows us to evaluate their historical age assignments and, in one case, refute a previously assigned provisional age.

Harrat Rahat volcanic field