USGS Science⌕ Search

SEARCH · USGS Science

Results for “Data Series”

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.

887 records · Page 39Linked to original sources
USGS2025

Shallow geologic framework of the Mississippi Sound and the potential for sediment resources

The Mississippi Sound, an estuarine environment located between the mainland and barrier islands bordering the northern Gulf of America (formerly the Gulf of Mexico), serves as a vital ecosystem for the States of Mississippi and Alabama. Spanning approximately 100 kilometers from east to west and covering 1,400 square kilometers, the sound is home to marine industry and ports, and its shallow and brackish waters sustain a diverse array of marine life. Barrier islands along the southern edge of the sound separate the microtidal estuary from the Gulf of America. This protection from gulf wave action mediates current flow within the sound, resulting in predominantly fine-grained sediment deposition along the seafloor. This study, conducted by the U.S. Geological Survey in cooperation with the U.S. Army Corps of Engineers, provides insight on fluvial and tidal processes spanning the past 5,000 years. The report synthesizes existing research to provide a comprehensive overview of the sound geology, from Pleistocene origins to present-day morphology, and utilizes high-resolution single channel seismic profiles and sediment data to identify and map sedimentary deposits and morphologic features at and below the seafloor. Despite its ecological significance, the Mississippi Sound faces environmental challenges, including water-quality issues, habitat degradation, storm-induced erosion, and the ongoing threats of sea-level rise and environmental changes. This study uses the present-day understanding of the sound's geology to inform coastal management decisions, hazard assessment, and potential mineral resources.

Louisiana, Mississippi↗
USGS2024

Reference 1D seismic velocity models for volcano monitoring and imaging: Methods, models, and applications

Seismic velocity models of the crust are an integral part of earthquake monitoring systems at volcanoes. 1D models that vary only in depth are typically used for real‐time hypocenter determination and serve as critical reference models for detailed 3D imaging studies and geomechanical modeling. Such models are usually computed using seismic tomographic methods that rely on P ‐ and S ‐wave arrival‐time picks from numerous earthquakes recorded at receivers around the volcano. Traditional linearized tomographic methods that jointly invert for source locations, velocity structure, and station corrections depend critically on having reasonable starting values for the unknown parameters, are susceptible to local misfit minima and divergence, and often do not provide adequate uncertainty information. These issues are often exacerbated by sparse seismic networks, inadequate distributions of seismicity, and/or poor data quality common at volcanoes. In contrast, modern probabilistic global search methods avoid these issues only at the cost of increased computation time. In this article, we review both approaches and present example applications and comparisons at several volcanoes in the United States, including Mount Hood (Oregon), Mount St. Helens (Washington), the Island of Hawai’i, and Mount Cleveland (Alaska). We provide guidance on the proper usage of these methods as relevant to challenges specific to volcano monitoring and imaging. Finally, we survey‐published 1D P ‐wave velocity models from around the world and use them to derive a generic stratovolcano velocity model, which serves as a useful reference model for comparison and when local velocity information is sparse.

Seismological Research Letters↗
USGS2026

Escherichia coli monitoring and assessment in 2022 and 2023 after beach restoration at Lake St. Clair Metropark Beach, Macomb County, Michigan

Lake St. Clair Metropark Beach in Michigan has a history of closures because of elevated Escherichia coli ( E. coli ) concentrations in its recreational waters. To reduce closures, restoration projects were implemented in 2021 to deter waterfowl from congregating on the beach. In this study, the U.S. Geological Survey, in cooperation with the Michigan Department of the Environment, Great Lakes, and Energy and in collaboration with Huron-Clinton Metroparks and the Macomb County Health Department, monitored E. coli from 2022–23 in surface water, shallow groundwater, and sediment at Lake St. Clair Metropark Beach. Results were compared to data from a prerestoration (2018–19) study. A significant decrease in daily geometric mean E. coli concentrations in surface water was observed postrestoration, but the number of high concentration events increased. This resulted in more frequent beach closures postrestoration. Surface-sediment E. coli concentrations significantly decreased after restoration, and waterfowl populations generally decreased from 2021 to 2023, suggesting that the deterrence measures could be influencing E. coli concentrations in surface sediments and surface water. Groundwater E. coli concentrations were orders of magnitude higher than those in surface water and revealed no change correlated with restoration. Seepage measurements indicated that groundwater occasionally discharges into surface water, potentially providing a transport mechanism for E. coli to reach the lake. Continued monitoring and consideration of environmental factors could help to better understand the beach system.

Michigan↗
USGS2024

Preliminary observations of the April 5th, 2024, Mw4.8 New Jersey earthquake

On 5 April 2024, 10:23 a.m. local time, a moment magnitude 4.8 earthquake struck Tewksbury Township, New Jersey, about 65 km west of New York City. Millions of people from Virginia to Maine and beyond felt the ground shaking, resulting in the largest number (>180,000) of U.S. Geological Survey (USGS) “Did You Feel It?” reports of any earthquake. A team deployed by the Geotechnical Extreme Events Reconnaissance Association and the National Institute of Standards and Technology documented structural and nonstructural damage, including substantial damage to a historic masonry building in Lebanon, New Jersey. The USGS National Earthquake Information Center reported a focal depth of about 5 km, consistent with a lack of signal in Interferometric Synthetic Aperture Radar data. The focal mechanism solution is strike slip with a substantial thrust component. Neither mechanism’s nodal plane is parallel to the primary northeast trend of geologic discontinuities and mapped faults in the region, including the Ramapo fault. However, many of the relocated aftershocks, for which locations were augmented by temporary seismic deployments, form a cluster that parallels the general northeast trend of the faults. The aftershocks lie near the Tewksbury fault, north of the Ramapo fault.

New Jersey↗
USGS2026

Decadal-scale characteristics of natural and anthropogenic dune morphology along North Carolina Barrier Islands (SE Atlantic Coast)

Dunes offer protection to coastal communities from storms, yet quantitative research characterizing human influences on dune systems across decades and 10–100 km spatial scales constitutes a known data gap. This study analyzed a multidecadal archive of dune crest elevations and positions to map mesoscale spatiotemporal trends in foredune morphology across the barrier island coastline of North Carolina, SE Atlantic coast. Results are summarized within four categories of coastline representing a range of coastal management intensities, integrating beach width and dune lateral accommodation space to assess the comparative influence of differing intensities of anthropogenic intervention on dune morphology. These categories account for the presence of structures, roads, and alongshore connectivity to anthropogenic sediment sources (beach nourishment). These categories show distinct dune and beach system characteristics. Dunes are largest in the areas with the longest legacy of coastal management. As categories decrease in development intensity, dune and beach system characteristics generally depart from the large, static dunes and narrow beaches consistent with hold-the-line management. Where small or transgressive dunes exist on heavily developed coasts, they coincide with relatively strong hydrodynamic forcing or sediment-limiting geology. There is evidence of indirect coastal management influence on undeveloped barrier islands with alongshore connectivity to developed barriers. While this analysis was conducted at the mesoscale, the results provide researchers with the framework to make insights at smaller scales than those demonstrated in this study.

North Carolina↗
USGS2024

Global geologic map of Europa

Discovered by Galileo Galilei more than 400 years ago and imaged in detail by the Voyager 2 Galileo spacecraft, Jupiter’s icy moon Europa has been a source of intrigue. A range of science investigations indicate that it contains the key ingredients for habitability, notably energy, chemistry, and liquid water. Europa’s surface is geologically complex and, based on the dearth of impact craters, interpreted to be as young as ~60 Ma. The array of geologic features that characterize the surface include extensive ridged plains, regions of broad disruption termed chaotic terrain, long, quasi-linear ridges that span thousands of kilometers, and bands as much as 60 kilometers wide and that extend 100s of kilometers. These features, along with other geophysical measurements, indicate the presence of a global briny liquid water ocean beneath the ice shell. It was not until the arrival of the Galileo spacecraft in 1995 that the true nature and level of complexity of the surface was revealed. Although image data returned by Galileo provided insight into the structure of a variety of regions, the entire satellite has yet to be observed at a regional scale (less than 250 meters per pixel) and the detailed geologic nature of much of its surface remains a mystery. Establishing the global context of the distribution and timing of Europan geologic units forms a basis to understand regional and local scale processes, serves as a tool for the planning of future missions, and most of all is essential to gaining insight into the potential habitability of this icy world.

Scientific Investigations Map↗
USGS2026

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

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

California↗
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↗
USGS2025

Linking stream-reach nitrogen loads and groundwater “reachsheds” to inform wastewater-nitrogen management actions, Cape Cod, Massachusetts

Study Region Cape Cod, Massachusetts, U.S.A. Study Focus Anthropogenic nitrogen (N) is a key factor in degrading groundwater and surface-water quality, particularly in coastal New England where onsite wastewater systems are prevalent. This study evaluated whether direct N-load measurements in streams on Cape Cod, Massachusetts, coupled with flow-path information from groundwater-flow models, can effectively identify potential land areas where nitrogen mitigation could substantially reduce loads to receiving waters. Nitrogen fluxes were measured along 63 stream reaches during winter and summer and paired with simulated groundwater recharge areas to identify and rank potential areas for reduction of nonpoint-source N inputs. New Hydrologic Insights for the Region Reach-scale nitrate-N loads ranged from −39.1–1182 kg-N/yr per 100 m of stream, indicating spatially variable groundwater inputs across seasons. “Reachsheds” — areas contributing groundwater recharge to specific stream reaches — were delineated using a regional groundwater-flow model. Strong correlations were found between observed N loads and land-use characteristics, especially the number of septic systems and total N inputs from the sum of considered sources. Observed N loads were moderately correlated with recharge area size and wastewater flow estimates. Correlating reach-specific groundwater N loads with land use and parcel-scale nitrogen-yield data identified reachsheds with the highest potential for N load reduction. This approach enables targeted implementation of restoration efforts to optimize nutrient management and support regional load reduction.

Massachusetts↗
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↗
USGS2026

Black rails in Baja California

Summarizing records from 1905 to 2025, we document the presence and distribution of the California black rail ( Laterallus jamaicensis coturniculus ) in northwestern México, particularly Baja California and the Colorado River Delta. This subspecies has a restricted distribution in the western U.S. and northwestern México. We found few records over the past century: sparse detections in the early 1900s and an increase in reported records after 2000. Over 120 years, we compiled 35 documented records of California black rail occurrence in Baja California with a minimum reported count of 56 rails. Our call-broadcast survey methods reconfirmed previously occupied sites and documented the species at a previously unreported site. Surveys in the Baja California portion of the Colorado River Delta yielded 20 detections across nine survey events between 2000 and 2022. In April 2003, standardized surveys detected 12 California black rails: 7 in Bahía de San Quintín, 3 in Arroyo San Telmo, and 2 in Arroyo El Rosario. Our recent survey work (2022–2025) recorded maximum monthly counts of 31 detections at Bahía de San Quintín and 15 detections at Estero de Punta Banda, the latter representing a previously undocumented location for the subspecies. These surveys also yielded the first photographic and acoustic documentation of California black rail in México. The contribution of spatial and temporal data provided by citizen science platforms increased our understanding of the subspecies’ distribution. Enforceable protection of this subspecies’ habitat is needed to safeguard its persistence.

Baja California↗
USGS2024

Tungsten skarn quantitative mineral resource assessment and gold, rare earth elements, graphite, and uranium qualitative assessments of the Kuldjuktau and Auminzatau Ranges, in the central Kyzylkum region, Uzbekistan

A new quantitative mineral resource assessment for tungsten skarn was conducted for the Auminzatau and Kuldjuktau mountain ranges in Central Uzbekistan, along with qualitative assessments of orogenic gold, rare earth elements (REEs), amorphous graphite, and uranium. By integrating a variety of geological, geochemical, geophysical, and remote sensing data sets, estimates of undiscovered tungsten skarn deposits in permissive tracts are combined with grade and tonnage distributions of known deposits to generate probabilistic estimates of undiscovered resources. Undiscovered deposits in Auminzatau are estimated to contain median resources of 98 thousand metric tons (kt) of WO 3 with a 70 percent (%) probability of at least 28 kt and a 10% probability of at least 468 kt, of which 16 kt to 293 kt may be economic to extract. In Kuldjuktau, the undiscovered deposits are estimated to contain median resources of 27 kt of WO 3 with a 60% probability of at least 12 kt and a 10% probability of at least 208 kt, of which 5 kt to 132 kt may be economic to extract. Our results suggest that the Auminzatau–Kuldjuktau Mountains area is highly prospective for additional discovery of significant Au and U resources and has low prospectivity for discovery of significant REE and graphite resources.

Minerals↗
USGS2025

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

Flood-inundation maps for Río de la Plata in and near Comerío, Puerto Rico, 2025

Digital flood-inundation maps for a 3.1-mile reach of Río de la Plata in and near Comerío, Puerto Rico, were created by the U.S. Geological Survey (USGS). Water-surface profiles were computed for the stream reach by using a one-dimensional steady-state step-backwater model. The model was calibrated to the current (2025) stage-streamflow relation (rating curve 11.0) for the USGS streamgage 50043800, Río de la Plata at Comerío, Puerto Rico. The resulting hydraulic model was then used to compute 16 water-surface profiles for water levels (flood stages) ranging from 10.00 to 40.00 feet at the streamgage and ranging from “action stage” to above “major flood stage” as reported by the National Weather Service. The 40.00-foot stage was selected because it exceeds the peak stage of 34.86 ft recorded during Hurricane Maria at the USGS streamgage 50043800, Río de Plata at Comerío, Puerto Rico. The simulated water-surface profiles were then used in combination with a digital elevation model derived from light detection and ranging data to map the inundated areas associated with each flood profile. The flood-inundation maps and the supporting hydraulic model produced by this study can be used by emergency managers and local officials to assess flood mitigation strategies and to define flood hazard areas to help protect life and property, to coordinate flood response activities such as evacuations and road closures, and to aid post-flood recovery efforts.

Comerio↗
USGS2024

The dominance and growth of shallow groundwater resources in continuous permafrost environments

Water is a limited resource in Arctic watersheds with continuous permafrost because freezing conditions in winter and the impermeability of permafrost limit storage and connectivity between surface water and deep groundwater. However, groundwater can still be an important source of surface water in such settings, feeding springs and large aufeis fields that are abundant in cold regions and generating runoff when precipitation is rare. Whether groundwater is sourced from suprapermafrost taliks or deeper regional aquifers will impact water availability as the Arctic continues to warm and thaw. Previous research is ambiguous about the role of deep groundwater, leading to uncertainty regarding Arctic water availability and changing water resources. We analyzed chemistry and residence times of spring, stream, and river waters in the continuous permafrost zone of Alaska, spanning the mountains to the coastal plain. Water chemistry and age tracers show that surface waters are predominately sourced from recent precipitation and have short (<50 y) subsurface residence times. Remote sensing indicates trends in the areal extent of aufeis over the last 37 y, and correlations between aufeis extent and previous year summer temperature. Together, these data indicate that surface waters in continuous permafrost regions may be impacted by short flow paths and shallow suprapermafrost aquifers that are highly sensitive to climatic and hydrologic change over annual timescales. Despite the lack of connection to regional aquifers, continued warming and permafrost thaw may promote deepening of the shallow subsurface aquifers and creation of shallow taliks, providing some resilience to Arctic freshwater ecosystems.

Proceedings of the National Academy of Sciences↗
USGS2025

The Mammoth magnetic anomaly, Pinal County, Arizona

A high-resolution Earth Mapping Resources Initiative airborne geophysical survey was flown in the southwest North American porphyry copper province to improve bedrock geologic maps and to identify areas that have unrecognized critical mineral resource potential. During the review of the aeromagnetic data, a distinctly monopolar-shaped, negative magnetic anomaly was observed at a flight elevation of 200 m above the ground with a maximum amplitude of –9500 nT. We have named this the Mammoth magnetic anomaly (MMA) because it is centered 12 km northeast of the town of Mammoth, Arizona, USA. The total field anomaly (TFA) contour of –500 nT enclosing the MMA defines an elongate shape measuring 2.5 km long by 1 km wide that trends northwest–southeast. Given the striking nature of this negative, monopolar-shaped magnetic anomaly, we conducted a ground campaign in May 2025 to determine its authenticity and potential relationship to critical mineral endowment. The MMA was confirmed on the ground with a TFA approaching –46,000 nT. Total magnetic intensity (TMI) observations routinely fell below the 18,000 nT operating floor of an industry-standard cesium-vapor total field magnetometer, and extremely low TMI measurements were corroborated along coincident traverse lines using two high dynamic range, but lower sensitivity, smartphone vector magnetometers. The lowest TMI values recorded by both smartphone magnetometers were 1000 nT and confirmed with multiple adjacent and crossing lines. Field observations suggest that this magnetic feature is caused by strong remanent magnetization within fine-grained magnetite hosted within locally altered Pinal Schist.

Arizona↗
USGS2023

Eruptive history within the vicinity of Al Madīnah in northern Harrat Rahat, Kingdom of Saudi Arabia

The northernmost part of the Harrat Rahat volcanic field contains early Pleistocene to Holocene mafic eruptive products within the vicinity of the city of Al Madīnah, Kingdom of Saudi Arabia. A detailed geologic investigation into the eruptive history of a 570 square kilometer (km 2 ) area covering Al Madīnah and the surrounding area has yielded 33 mapped Quaternary volcanic units consisting of lava flows, scoria cones, and shield volcanoes. These eruptive products consist of continental, intraplate alkalic and minor transitional basalts, hawaiites, and a single mugearite that were emplaced from at least 1,014±14 thousand years ago (ka) to a single Holocene eruption in 1256 C.E. Lava flows are generally 10 to 15 kilometers (km) long (but can reach 23 km long), 1 to 3 km wide, and at least 10 meters thick. Most of the mapped units erupted episodically between 400 and 340 ka and 180 and 100 ka. Despite small individual volumes (less than 1 cubic kilometers dense rock equivalent), each unit represents eruption of a distinct magma batch that was strongly influenced by clinopyroxene, olivine, and plagioclase fractionation. Some of these units are interpreted to have undergone magma mixing pre- and (or) syneruptively. Integrating eruption ages, geochemistry, and paleomagnetic data yields evidence that some eruptions were temporally and (or) spatially clustered. Aligned scoria cones and elongate vent edifices were constructed atop fissure vent systems that reflect the local stress field, which controls dike ascent through the middle and upper crust.

Professional Paper↗