SEARCH · USGS Science

Results for “Bulletin of the Seismological Society of America”

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

SourceResearch connectionFrom yearThrough year

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

1,690 records · Page 22Linked to original sources

USGS2026

Mineral chemistry perspective on remobilization of stored magma at Kamakai'a Hills, Southwest Rift Zone of Kilauea, Island of Hawai'i, USA

Differentiated magmas stored in the rift zones of Kīlauea have received more attention in recent years following eruption of andesite during the early phase of 2018 lower East Rift Zone activity. Despite this growing interest, some of the most voluminous eruptions of differentiated rift zone magmas remain poorly studied. One such eruption, and the most voluminous exposed differentiated flow field at Kīlauea, is the Kamakaiʻa Hills. This eruption took place in the Southwest Rift Zone of Kīlauea, a region that is hypothesized to contain a long-lived rift zone reservoir. The Kamakaiʻa Hills flow field encompasses >250 × 10 6 m 3 of basaltic andesite and basalt compositions with a mineral assemblage of orthopyroxene + clinopyroxene + plagioclase during its early ʻaʻā phase and clinopyroxene + plagioclase + olivine during its late pāhoehoe phase. To better understand storage conditions and magma accumulation, this study focuses on major, minor, and trace elements from the mineral assemblage present within the early ʻaʻā and late pāhoehoe phases. The diversity of clinopyroxene and plagioclase compositions within the early ʻaʻā and late pāhoehoe phases, as well as diverse compositions of plagioclase and orthopyroxene within the early ʻaʻā phase, suggest multiple magma bodies and limited pre-eruption magma mixing within the broader Kamakaiʻa Hills reservoir. Oscillatory zoning patterns (particularly in clinopyroxene) imply processes such as recharge events, magma mixing or mingling, or convection within a differentially cooling, chemically stratified reservoir over protracted time intervals, whereas only limited resorbed mineral textures indicate incomplete mixing of heat and chemically distinct magmas during the dike intrusion that triggered the eruption. Mineral-mineral and mineral-melt thermobarometry indicate predominantly shallow (≤2.5 km depth) crustal storage conditions of the cooled, differentiated magma (∼1100 °C and cooler for the basaltic andesites) to hotter temperatures for the basalts (all >1100 °C). Despite the known large standard errors estimated for mineral-melt and mineral-mineral barometry (10s to >100 MPa), the calculated pressures and depths broadly correspond with earthquake swarm depths beneath the Kamakaiʻa Hills, and drill core and fluid inclusion barometry storage depths of differentiated magmas within the lower East Rift Zone. The Kamakaiʻa Hills differentiated magmas have H 2 O contents (∼0.5 wt%, using plagioclase-melt hygrometry) equivalent to typical Kīlauea basalts. Our data and interpretations demonstrate a complex, long-lived rift zone storage system that consisted of multiple magma bodies and was mobilized into eruption through intrusion of a hotter and more primitive summit-derived (uprift) magma.

Hawaii↗

USGS2025

Reconnaissance of potential alternate water supply sources for the City of Gary, West Virginia

Seven potential sources of water, consisting of free-flowing discharge from abandoned coal mines at six locations and one abandoned flooded underground coal mine air shaft, were sampled for chemical analysis to assess the quality of the groundwater emanating from the seven mine sources. The six free-flowing mine discharge sources were also assessed for discharge by current-meter measurements on two separate occasions. The U.S. Geological Survey assessed these seven sources to provide information to the City of Gary, West Virginia (W. Va.), and the City of Gary’s consulting engineer with groundwater-quality and flow data to allow them to assess the seven sites as potential alternate sources of water for the City of Gary to augment its existing supply. For the six sites where discharge could be measured, discharge ranged from a minimum of 0.082 cubic feet per second (ft 3 /s) to a maximum of 3.685 ft 3 /s. Of the six sites measured, only two, Harmon Branch at Thorpe, W. Va. (USGS site 372201081303501) and the abandoned public-supply water wells near Havaco, W. Va. (USGS site 372358081344601), had discharge in excess of 1.00 ft 3 /s. Discharge from the abandoned public supply wells was 3.685 ft 3 /s on September 20, 2023, and 2.888 ft 3 /s on October 16, 2023, and discharge from Harmon Branch at Thorpe, W. Va., was 1.049 ft 3 /s on September 22, 2023, and 1.038 ft 3 /s on October 17, 2023. Discharge in the abandoned underground mine air shaft (USGS site 372224081340901) could not be assessed, but the air shaft drains an abandoned mine that likely contains water stored in approximately 1.7 square miles (mi 2 ) of abandoned underground coal mines in the Pocahontas No. 3 coal seam, and possibly an additional 0.9 mi 2 of leakage from the overlying Pocahontas No. 4 coal seam. Discharge for the six sites measured for the study was measured during a period between September 20 and October 18, 2023, and corresponded to the 12th to the 15th percentile of flow-duration statistics for the Tug Fork downstream of Elkhorn Creek at Welch, W. Va. streamgage (USGS site 03212750). Water-quality data for the seven sites sampled overall were acceptable with respect to drinking water standards. Of the 203 constituents analyzed, only a few failed to meet applicable U.S. Environmental Protection Agency (EPA) drinking water standards. Iron exceeded the 300 micrograms per liter (μg/L) secondary maximum contaminant level (SMCL) at only 1 of the 7 sites (14.3 percent) sampled. Iron concentrations ranged from a minimum of less than (<) 5.00 μg/L to a maximum of 724 μg/L with a median concentration of 7.62 μg/L. Manganese exceeded the 50.0 μg/L SMCL at 2 of the 7 sites (28.6 percent) sampled. Manganese concentrations ranged from a minimum of 1.93 μg/L to a maximum of 271 μg/L with a median concentration of 4.03 μg/L. No sites sampled exceeded the arsenic maximum contaminant level (MCL) of 10 μg/L. Arsenic concentrations ranged from a minimum of <0.100 μg/L to a maximum of 2.35 μg/L with a median arsenic concentration of 0.200 μg/L. None of the seven sites sampled for selenium for this study exceeded the EPA MCL of 50.0 μg/L. Selenium concentrations ranged from a minimum of <0.050 μg/L to a maximum of 5.26 μg/L with a median concentration of 3.21 μg/L. All seven sites were sampled for volatile organic compounds (VOCs), semivolatile organic compounds (SVOCs), and polychlorinated biphenyls (PCBs), but most had concentrations below the detection limit. Of the 10 PCB compounds analyzed for the seven sites sampled, none contained detectable concentrations of PCBs or Aroclor compounds. Of the 44 SVOCs analyzed at each of the seven sites sampled, only 1 SVOC, acenaphthene, was detected, at a concentration of 0.02 μg/L. Of the 96 VOCs analyzed, from each of the seven sites sampled, only two were found at detectable concentrations. Trichloromethane was detected only at 1 of the 7 (14.3 percent) sites sampled at a concentration of 0.027 μg/L, and benzene was detected at the same site and 3 additional sites (4 of the 7 sites or 57.1 percent of the sites sampled) at concentrations of 0.028, 0.029, 0.021, and 0.035 μg/L, but none exceeded the EPA MCL for benzene of 5.00 μg/L. Total coliform bacteria are ubiquitous in the environment, and their presence only suggests the potential for contamination by near-surface processes. Escherichia coli ( E. coli ) bacteria are derived from either human or animal fecal material and can be an indicator of potential contamination by pathogenic bacteria or viruses. Total coliform bacteria were detected at all 7 sites sampled at concentrations ranging from 17.5 to greater than (>) 2,420 most probable number per 100 mL (MPN/100 mL) of sample, with a median total coliform concentration of 1,553 MPN/100 mL. Escherichia coli bacteria were detected at 4 of the 7 sites sampled at concentrations ranging from <1 to 11.9 MPN/100 mL, with a median E. coli concentration of 5.1 MPN/100 mL.

West Virginia↗

USGS2025

The effects of line simplification on planform geometry

Data on maps should retain accuracy regardless of scale. Yet, as cartographic lines are generalized, there can be impacts on properties such as topology, density, and planform geometry. Here, we investigate the use of the Scale Specific Sinuosity (S3) metric (Stanislawski et al., 2023) to evaluate the effects of line simplification on planform geometry, which is the bends of streams in map view. We employ an open-source Python S3 workflow to characterize the geometry of five diverse stream channels in the United States. The original data are extracted from the U.S. Geological Survey National Hydrography Dataset 1:24,000-scale vector data (U.S. Geological Survey, 2000) (Table 1), and the simplification is done using the Visvalingam and Whyatt method (2017) with a simplification tolerance of 0.5, 1.0, 1.5, and 2.0 km. The S3 analysis is calculated at each level of simplification and S3 derivatives are generated. Derivatives include measures of sinuosity, fractal dimension, and the dominant bend wavelength. The findings show that the change in planform geometry is scale-dependent, though simplification will have little effect on straighter lines. The change becomes more apparent in complex lines as the degree of simplification aligns with the scale of the dominant bend geometries. These logical conclusions are evidence that the S3 is a useful metric for automated characterization of bend geometry regardless of line complexity.

Conference Paper↗

USGS2026

Colored shaded-relief bathymetric and acoustic-backscatter maps of Jenkinson Lake with orthomosaic of the Sly Park Creek and Hazel Creek area, California

The Caldor Fire was ignited on August 14, 2021, and burned almost 222,000 acres (898 square kilometers) in forested terrain of the central and western Sierra Nevada, California. During the subsequent two months, the fire burned nearly all of Sly Park Creek watershed in El Dorado County. The El Dorado Irrigation District manages the water supply for the area using storage in Jenkinson Lake, a 1.6-kilometer- (1.0-mile-) wide and 3.6-kilometer- (2.2-mile-) long reservoir, located south of the town of Pollock Pines. Several weeks after the fire, the U.S. Geological Survey began investigations into post-fire landscape responses, including sediment yield, by measuring new sediment deposition in Jenkinson Lake. This study focused on the collection and processing of bathymetric and acoustic-backscatter data, as well as onshore aerial imagery in and around Jenkinson Lake, to support wildfire science after the Caldor Fire. A colored shaded-relief bathymetric map (sheet 1) and an acoustic backscatter map (sheet 2) show the lake floor morphology and backscatter intensities.

California↗

USGS2025

Apatite geo-thermochronology and geochemistry constrain Oligocene-Miocene growth and geodynamics of the northeastern Tibetan Plateau

Understanding the geodynamics of plateau evolution requires examining the spatial and temporal aspects of mountain building in the northeastern Tibetan Plateau, which are still under debate. Here we integrate apatite geo-thermochronological and geochemical data from the Oligocene-Miocene succession of the Xunhua Basin to elucidate the evolution of the regional topography. The results suggest sediment provenance changes, at ca. 28, 20 and 12 Ma, indicating topographic growth of the West Qinling, Laji Shan, and Jishi Shan source areas, respectively. Our findings reveal Oligocene-Miocene stepwise deformation and middle Miocene stress reorganization within the northeastern Tibetan Plateau. We attribute this to Oligocene lithospheric foundering beneath the south-central Tibetan Plateau and the middle Miocene initiation of sinistral strike-slip faults due to the underthrusting of India and North China. This study highlights the roles of lithospheric removal and continental underthrusting in plateau growth, and the value of multi-proxy apatite analyses in provenance studies.

Northeastern Tibetan Plateau↗

USGS2025

New constraints on location and timing of the Great Lakes tectonic zone, central Upper Peninsula, Michigan, USA

The Great Lakes tectonic zone (GLTZ) forms the boundary between the Wawa–Abitibi and Minnesota River Valley subprovinces within the Archean Superior Province. The GLTZ is concealed for all of its 1100 km length, except for a segment in the central Upper Peninsula of Michigan. There, it is exposed as a northwest-striking mylonite zone along a 11 km segment, extending to the onlap of Paleozoic rocks to the east. Farther east, its location has been unknown. Here, we use aeromagnetic and gravity data to develop interpretations of the expression of the GLTZ and to define its extent under cover. Aeromagnetic gradients over the mylonite zone are interpreted to be produced by structurally juxtaposed rocks with varying magnetizations. Gravity data show a regional gradient along the GLTZ, produced by the juxtaposition of a dense greenstone belt on the north against lower-density gneisses and granites on the south. The GLTZ is interpreted to extend ∼55 km under cover to the east. The GLTZ is terminated on the east by the buried eastern arm of the ca. 1100 Ma Midcontinent Rift. An undeformed granitic dike that cuts the mylonitic foliation produces a U–Pb apatite age of 2523 ± 33 Ma, implying no major post-Archean shearing occurred, and is at odds with previous interpretations of major Proterozoic reactivation. A granite intrusion in the Minnesota River Valley subprovince produces a Pb–Pb zircon age of 2606.9 ± 3.6/7.4 Ma. This suggests that magmatism related to the Sacred Heart orogeny, previously known in Minnesota, extended to Michigan.

Michigan↗

USGS2026

Regional characterization of coal resources in the Gulf Coast Coal Region—In support of coal as a potential source for critical minerals and rare earth elements

There is increasing interest in extracting critical minerals (CM), including rare earth elements (REE), from coals in the United States to address the overreliance on imported REE. The Gulf Coast and the Fort Union regions are the two major lignite-bearing regions within the United States. Recent REE and CM studies of the lignite in these regions have indicated that the coals may be a viable source material for REE and CM extraction. To evaluate in-place coal as a potential source of REE and CM, the coal resources need to be quantified. This study presents the results of a regional analysis of Gulf Coast Coal Region lignite and bituminous coal resources that might be available as potential sources of REE and CM. The resource analysis used kriging methods to develop isopleth maps of cumulative coal thickness throughout the region using data from drill holes and measured mine and outcrop sections for a total of 31,181 data points. The estimated total remaining coal resource in the Gulf Coast Coal Region is about 83 billion metric tons in the upper 90 meters (approximately 300 feet) of the subsurface. The current estimate of remaining coal resources includes approximately 1.88 billion metric tons of coal that have been extracted in the region since 1983. Texas accounted for 39 percent (32 billion metric tons) of the remaining total resource, followed by Mississippi (24 percent, 20 billion metric tons ), Louisiana (14 percent, 12 billion metric tons), Tennessee (10 percent, 8.5 billion metric tons), and Arkansas (6 percent, 5.1 billion metric tons). Alabama, Georgia, Kentucky, and Missouri each accounted for less than 5 percent of the total resource; Georgia had the smallest resource estimated at 7 million metric tons. Here we report the first known statewide lignite resource estimates for Georgia, Kentucky (820 million metric tons), and Missouri (1,800 million metric tons). A comparison of the results of this study with those of previous Gulf Coast Coal Region and Fort Union Coal Region resource studies is difficult because each study used different data sources, assessment methodologies, overburden depths, and qualifying coal thicknesses.

Alabama, Arkansas, Florida, Georgia, Illinois, Lou↗

USGS2024

Chloride concentrations in groundwater from the western part of the Southern Hills regional aquifer system, Louisiana, 2021–22

Groundwater is heavily used for public supply and industrial uses in the Baton Rouge, Louisiana, area. Lowered water levels resulting from groundwater withdrawals have induced the movement of saltwater towards wells in East Baton Rouge and West Baton Rouge Parishes. Saltwater intrusion has the potential to affect water supply infrastructure, reduce water availability for some uses, and increase treatment costs. To document current conditions, samples were collected from 161 wells screened in 10 aquifers of the Southern Hills regional aquifer system during November 2021 through February 2022. The results were compared with historical data to identify where chloride concentrations are increasing, which could indicate that saltwater intrusion is occurring. Saltwater intrusion, to varying degrees and areal extents, was observed in most of the 10 aquifers. The limited availability of monitoring wells near or within some of the known saltwater plume areas restricts tracking of the movement or delineation of the plumes’ current extents.

Louisiana↗

USGS2024

Simulated effects of projected 2014–40 withdrawals on groundwater flow and water levels in the New Jersey Coastal Plain

Abstract Groundwater flow between 2014 through 2040 was simulated in the New Jersey Coastal Plain based on three withdrawal scenarios. Two of the scenarios were based on projected population trends and the assumption of water conservation; the nominal water-loss scenario projected a status quo in the efficiency of water loss in the delivery systems whereas the optimal water-loss scenario projected a better water-loss efficiency resulting in less withdrawals. The third scenario assumes that all wells will withdraw water at their full allocation level which is generally much more than reported withdrawals in 2013 or projected under the other two scenarios. Maps and summaries of heads and drawdowns are presented for nine confined aquifers. All the aquifers have areas with heads below sea level by 2040. Of the three scenarios, the drawdowns are most extreme in the full allocation scenarios; there are large areas of head decline greater than 20 feet in 5 of the 9 confined aquifers. The exceptions are the Vincentown aquifer, despite some areas of large drawdown in the vicinity of wells, and the three Potomac-Raritan-Magothy (PRM) aquifers where withdrawals are regulated by Critical Area restrictions. The nominal and optimal water-loss scenarios have some areas of head declines; most are less than 15 feet. The simulation of these scenarios shows some extensive areas of head recovery as well—especially in the aquifers that are regulated by the Critical Area restrictions. Budgets of inflow and outflow components were calculated for 44 hydrologic budget areas (HBAs). The budget analysis shows that the water movement is complex and varies based on the aquifer geometry and location of pumping wells. Flow components between the unconfined and confined parts of the system were summarized by HUC11 (hydrologic unit code 11) basins.

New Jersey↗

USGS2025

Estimation of dynamic geologic CO2 storage resources in the Illinois Basin, including effects of brine extraction, anisotropy, and hydrogeologic heterogeneity

Since the vast majority of carbon dioxide (CO 2 ) storage resources in the United States are in deep saline aquifers, optimizing the use of these saline storage resources could be crucial for efficient development of geologic CO 2 storage (GCS) resources and basin- or larger-scale deployment of GCS in the country. Maximum CO 2 injection rates can be enhanced by extracting brine from the CO 2 storage unit. However, disposal of the extracted brine is both a technological and economic challenge. The lowest-cost option would likely be reinjection of the extracted brine into another formation above or below the CO 2 storage unit. Therefore, it is important to estimate brine injectivity as it will constrain the potential to increase CO 2 injectivity at an injection site that has access to multiple geologic storage units where either CO 2 or brine can be injected. Using a simulation-optimization framework, coupled with a non-isothermal, multiphase CO 2 -water-salt equation-of-state module, we developed a computationally efficient method for evaluating optimization of simultaneous CO 2 injection, brine extraction, and brine (re)injection at hypothetical injection sites deployed across a geologic basin. The Illinois basin is ideal for testing our methodology because it contains multiple geologic storage units with seals in between them to isolate injection of CO 2 in one unit from interfering with the injection of either brine or CO 2 in another unit above or below it. In addition, we investigated the relative effects of variation in key geologic parameters as well as two reservoir structures (hydrogeologic heterogeneity/anisotropy and homogeneity/isotropy) on CO 2 injectivities and enhancement of CO 2 injectivity through extracting brine. Results suggest that permeability, depth, and especially thickness of the storage unit could be the most influential parameters determining CO 2 injectivity. They also suggest that only injecting CO 2 into the storage unit with the greatest injectivity, enhancing that unit’s injectivity by extracting brine, and disposing of the produced brine in other suitable units could maximize total CO 2 injectivity in limited regions of the basin. At the majority of simulated injection sites, however, we found that injecting CO 2 into all of the accessible and suitable storage units was more likely to maximize the CO 2 storage resource.

Illinois, Kentucky, Indiana↗

USGS2024

Earthquake effects surveyed during the nineteenth century as ecological features of Chinookan tidelands

Lasting effects of a Cascadia earthquake in 1700 were documented during surveys of Chinookan tidelands near the mouth of the Columbia River between 1805 and 1868. The effects resemble estuarine consequences, near Anchorage, of the 1964 Alaska earthquake: fatal drowning of subsided meadows and forests by post-earthquake tides, rebirth of marshes and forests through post-earthquake sedimentation, and uplift. Chinookan remains of killed forests were recorded by James Graham Cooper, John J. Lowell, and Cleveland Rockwell. Cooper, attached to a railroad survey and the Smithsonian Institution, wrote of redcedar stumps and trunks standing dead in tidal marshes of Shoalwater (now Willapa) Bay. Two such snags served as bearing trees for Lowell as he platted a Shoalwater Bay township under contract with the General Land Office. Rockwell, of the US Coast Survey, flecked landward edges of tidal flats west of Astoria with symbols that evoke remains of a bygone spruce forest. The Lewis and Clark Expedition, while in that area in 1805–1806, mapped and puzzled over tideland vegetation that post-1700 succession helps explain.

Alaska, British Columbia, California, Oregon, Wash↗

USGS2025

The 3D Elevation Program—Supporting the economy of Massachusetts

Introduction Massachusetts extends from the mountains of the Appalachian system in the west of the State to the sandy beaches and rocky shorelines of the Atlantic coast in the east. Inland topographic data support a wide range of important activities, including geologic mapping, transportation planning, forest and wildlife management, quantifying ecological services, water supply protection, commonwealth-wide infrastructure planning, local site planning, and flood-plain management. Nearshore bathymetry can be used to support coastal portions of the Commonwealth by addressing the combined threats of ocean warming, strong storm surge, and rising sea levels. The maintenance and (or) expansion of Massachusetts ports (for instance, Boston, New Bedford) and Cape Cod sediment management depends upon the accurate mapping of bathymetry and the frequent influx of sediment and redeposition. Critical applications that address the broad range of requirements depend on light detection and ranging (lidar) data that provide a highly detailed three-dimensional (3D) model of the Earth’s surface and aboveground features. The 3D Elevation Program (3DEP) is managed by the U.S. Geological Survey (USGS) in partnership with Federal, State, Tribal, U.S. territorial, and local agencies to acquire consistent lidar coverage at quality level 2 or better to meet the many needs of the Nation and Massachusetts. The status of available and in-progress 3DEP baseline lidar data in Massachusetts is shown in figure 1. 3DEP baseline lidar data include quality level 2 or better, 1-meter or better digital elevation models, and lidar point clouds, and must meet the Lidar Base Specification version 1.2 ( https://www.usgs.gov/3dep/lidarspec ) or newer requirements. The National Enhanced Elevation Assessment identified user requirements and conservatively estimated that availability of lidar data would result in at least $1.23 million in new benefits annually to Massachusetts. The top 10 Massachusetts business uses for 3D elevation data, which are based on the estimated annual conservative benefits of 3DEP, are shown in table 2.

Massachusetts↗

USGS2025

Comparative properties of saponitic fault gouge and serpentinite muds cored from mud volcanoes of the Mariana subduction zone

We obtained 12 core samples for physical and chemical characterization from three serpentinite mud volcanoes (Yinazao, Asùt Tesoru, and Fantangisña) located on the forearc of the Mariana subduction system, that were drilled during International Ocean Discovery Program Expedition 366. Two samples from the Fantangisña mud volcano are interpreted to be clay-rich fault gouges derived from the subduction channel. Their bulk compositions are intermediate between the serpentinites and oceanic basalts. The oceanic crustal materials in the gouges have been thoroughly metasomatized and the serpentinites extensively altered to the trioctahedral, Mg-rich smectite clays saponite and corrensite. The only relict phases in clasts of crustal rock are accessory Ti- and P-bearing minerals. The two fault gouge samples have lower frictional strengths ( μ < 0.2) than the serpentinites ( μ = 0.2–0.4), and their measured permeabilities are also somewhat lower. Their physical and compositional properties correspond to saponitic gouges from other faults that juxtapose serpentinite against crustal rocks, in particular gouges from the two creeping traces of the San Andreas Fault recovered in the core from the San Andreas Fault Observatory at Depth. The décollement beneath Fantangisña mud volcano is thus expected to be very weak and likely characterized by stable slip.

Geophysics, Geochemistry, Geosystems↗

USGS1999

The mineral economy of Brazil--Economia mineral do Brasil

This study depicts the Brazilian government structure, mineral legislation and investment policy, taxation, foreign investment policies, environmental laws and regulations, and conditions in which the mineral industry operates. The report underlines Brazil's large and diversified mineral endowment. A total of 37 mineral commodities, or groups of closely related commodities, is discussed. An overview of the geologic setting of the major mineral deposits is presented. This report is presented in English and Portuguese in pdf format.

Data Series↗

USGS2025

Gravity and magnetic surveys of the Skaergaard intrusion, East Greenland

Aeromagnetic and gravity surveys of the Skaergaard intrusion in East Greenland were carried out in July–August 1971 as part of a grant to the University of Oregon Center for Volcanology to refine the models of crystallization and differentiation of the intrusion, specifically to test whether the intrusion is underlain by dense rocks of a reservoir 20 kilometers (km) thick (referred to as a “hidden zone”). The Skaergaard intrusion is a source of platinum group elements that are critical mineral resources for many technologies, and because no new data have been collected these legacy datasets remain a valuable asset. The total-intensity aeromagnetic survey was flown in early July 1971 with a proton precession magnetometer at a constant barometric altitude of 1.5 km (5,000 feet) with a nominal line spacing of 1 km. Two gravimeters were used to acquire 168 stations of which 86 were at known altitudes (mainly sea level) and 82 had altitudes measured by altimetry in late July–August 1971. Finally, a north-south ground vertical-intensity magnetic traverse was completed across the intrusion together with collection of oriented hand specimens. The hand specimens were measured for remnant magnetization and density, along with density measurements of more specimens collected by expedition geologists for other purposes. The intrusion is composed of layered gabbro with extensive crystal fractionation that is dense and strongly reversely polarized. After terrain correction and standard Bouguer gravity reduction, the gravity anomaly dataset was corrected for all rock above sea level using the density measurements of the various zones of the intrusion and the topographic and geologic maps (variable density Bouguer gravity reduction). A large regional gradient in the gravity anomaly data was removed using orthogonal polynomial fitting to the gridded data. The zonal volumes of rock below sea level were calculated from the dipping polygonal layer gravity model of the intrusion below sea level and combined with elliptic cross–section cylinders for the various zones above sea level to approximate the original zonal volumes of the intrusion. The residual gravity anomaly of 18–20 milligals (mGal) was only about half of the expected anomaly if a large hidden zone proposed from petrologic considerations were present, and both two-dimensional and three-dimensional models imply that the exposed series of intrusion zones explain the gravity anomaly by their down-dip extension below sea level together with a small hidden-zone volume. A three-dimensional model of the exposed rocks and their down-dip extension below sea level also can account for the aeromagnetic anomaly with little or no requirement for hidden-zone rock. The middle and upper zone units of the intrusion contain the most magnetite and account for most of the aeromagnetic anomaly.

Skaergaard intrusion↗

USGS2025

The geology of Canadian potash: A critical mineral for feeding the world

Potash, potassium-bearing water-soluble salt, is the primary global economic source of potassium. Potash is recognized as a critical mineral in Canada as it is the largest source of potassium used in fertilizers. It is essential for global agricultural productivity and food security. Canada is the world’s largest potash exporter with vast deposits in the widely mined Prairie Evaporite of Saskatchewan, which formed in the epicontinental Elk Point Basin during the Middle Devonian. Potash is also found in the Windsor Group of Atlantic Canada where it formed in a series of tectonically active basins during the Mississippian that have undergone substantial post-depositional subsurface alteration and deformation. Potash deposits were mined in New Brunswick up until 2016. Both deposits are salt giants, recording times in the geologic record of extensive and long-lasting evaporite genesis under arid conditions in restricted seas. This paper reviews the geological and economic significance of Canadian potash, including (1) the genesis of each deposit, (2) diagenetic, erosional, and tectonic modification, and (3) exploration and mining in each basin. Underdeveloped regions, possible undiscovered resources, environmental considerations, and the importance of sustainable practices in light of climate change and socioeconomic risks are also addressed.

Facets↗

USGS2024

Global variability of the composition and temperature at the 410-km discontinuity from receiver function analysis of dense arrays

Seismic boundaries caused by phase transitions between olivine polymorphs in Earth's mantle provide thermal and compositional markers that inform mantle dynamics. Seismic studies of the mantle transition zone often use either global averaging with sparse arrays or regional sampling from a single dense array. The intermediate approach of this study utilizes many densely spaced seismic arrays distributed around the globe. We systematically compute teleseismic P-to-S receiver functions for each seismic array and invert for the 1-D seismic velocity structure of the mantle transition zone beneath each array to facilitate a comparison between densely sampled regions. We stack 3,600 receiver functions on average at 67 arrays in total. The stack is used in a probabilistic inversion to estimate the mantle transition zone interface depths and velocities beneath each array. We focus on the 410-km discontinuity (410) because it is a prominent seismic interface that is clearly linked to a single mineral phase transition between olivine and wadsleyite. The depths and velocity contrasts of the 410 are mapped to temperatures and compositions using mineral physics constraints. The depth of the 410 ranges from ∼405–440 km, which is consistent with a ∼360 K temperature range in a dry mantle and a ∼260 K temperature range in a wet mantle (2 wt. % water). The Vs contrast across the 410 ranges from ∼2.5–8 %, which is consistent with ∼20–70 vol. % olivine composition in a dry mantle and ∼25–80 vol. % in a wet mantle. The bulk composition of the upper mantle near the 410-km discontinuity is typically considered to be well-mixed because there is no thermodynamic impediment to convection at the olivine to wadsleyite phase transition. However, the wide range of inferred olivine content from our study suggests that there are large lateral variations in the bulk composition of the upper mantle near the 410-km discontinuity.

Earth and Planetary Science Letters↗

USGS2024

Water-quality characteristics of the Red River of the North and tributaries in the Fargo-Moorhead metropolitan area, North Dakota, 2019–22

The Flood Risk Management Project was initiated in 2008 in the Fargo-Moorhead metropolitan area to reduce flood risk, flood damages, and flood protection costs in the Fargo-Moorhead metropolitan area. In cooperation with the U.S. Army Corps of Engineers, the U.S. Geological Survey initiated a water-quality monitoring study to describe the water-quality characteristics of the Red River of the North and its tributaries in the Fargo-Moorhead metropolitan area during the preconstruction period of the Flood Risk Management Project from October 1, 2019, to October 1, 2022. The monitoring study included the collection of discrete and continuous water-quality data and streamflow monitoring at selected sites that integrated and enhanced existing monitoring programs within the study area. Discrete samples collected at 10 sites in the Fargo-Moorhead metropolitan area were analyzed for major ions, trace elements, nutrients, suspended sediment, pesticides, and fecal indicator bacteria. In general, major ion concentrations were higher at sites on the tributaries (Wild Rice, Sheyenne, and Maple Rivers) compared to sites on the Red River of the North. In general, bicarbonate, calcium, magnesium, and sulfate represented most of the dissolved ions measured in samples collected at the 10 sites. Calcium, chloride, fluoride, potassium, silica, and sodium were also measured in samples, but they represented a smaller portion of the total dissolved ions. Sulfate was the most dominant dissolved ion that had the highest concentrations among the major ions measured in samples. A total of 18 trace elements were analyzed in discrete samples. Several of the trace elements had concentrations below the laboratory reporting level in all of the samples, including antimony, beryllium, cadmium, chromium, silver, and thallium. Sites on the Wild Rice River generally had the highest concentrations of arsenic, barium, boron, manganese, and nickel compared to the other sites. Nutrients analyzed in discrete samples included filtered and unfiltered concentrations of ammonia, nitrate plus nitrite, phosphorus, and organic carbon. The median filtered ammonia concentration at most sites was less than the laboratory reporting level of 0.03 milligram per liter as nitrogen except for the Sheyenne River at Harwood, North Dakota (U.S. Geological Survey [USGS] station 05060400), and Red River of the North near Georgetown, Minnesota (USGS station 05062130). The lowest median unfiltered nitrate plus nitrite concentration was measured at sites on the Red River of the North upstream from the Fargo-Moorhead metropolitan area and the highest median was at sites on the Red River of the North downstream from the Fargo-Moorhead metropolitan area compared to all other sites. The increase in nitrate plus nitrite concentrations could reflect the effect of the wastewater-treatment plant discharge that enters the Red River of the North upstream from the site located downstream from the Fargo-Moorhead metropolitan area and from urban runoff. Phosphorus (unfiltered) concentrations were generally higher at sites on the Maple and Sheyenne Rivers compared to the other sites and were higher at sites on the Red River of the North downstream from the Fargo-Moorhead metropolitan area compared to sites upstream on the Red River of the North. Suspended-sediment concentrations were generally highest at sites in the Sheyenne River and lowest in the upstream Red River of the North sites. Suspended-sediment concentration was highly variable in samples collected at the 10 sites, mostly influenced by the occurrence of snowmelt and rainfall-runoff events. The Sheyenne River near Kindred, N. Dak. (USGS station 05059000) had the largest range in sediment concentrations in samples collected at the 10 sites. For all sites other than the Sheyenne River near Kindred, N. Dak., 95 percent or more of the suspended sediment had particle diameter sizes less than 0.0625 millimeter in 50 percent of the samples (median). Of the 102 pesticides and pesticide degradates analyzed, 45 constituents had no detectable concentrations in any of the 17 samples collected at five sites. The remaining 57 pesticides had at least one detection in the samples collected at the five sites. The sites on the Wild Rice River (near Abercrombie, N. Dak., USGS station 05053000, and near St. Benedict, N. Dak., USGS station 05053500) and Sheyenne River near Kindred, N. Dak., had fewer pesticide detections compared to the Maple River below Mapleton, N.Dak. (USGS station 05060100) and the Red River of the North at Fargo, N. Dak (USGS station 05054000) and near Georgetown, Minn. Patterns in annual loads generally followed the same pattern as streamflow at the 10 sites for water years 2020–22. A water year is the 12-month period from October 1 to September 30 and is designated by the calendar year in which it ends. The greatest loads for all constituents were delivered at the two downstream sites on the Red River of the North; sites that also had the highest annual streamflows among the sites and the greatest loads were delivered in water year 2020 when the highest streamflows occurred at the sites. Likewise, the least loads for most constituents were at the Maple River and were least in 2021 compared to the other years because of low-streamflow conditions. Water-quality measurements continuously recorded at the Red River of the North at Hickson, N. Dak. (USGS station 05051522); Red River of the North at Fargo, N. Dak.; and Red River of the North near Georgetown, Minn. included water temperature, specific conductance, dissolved oxygen, pH, and turbidity. Specific conductance values were similar for the Red River of the North near Hickson, N. Dak., and Red River of the North at Fargo, N. Dak., when compared to the Red River of the North near Georgetown, Minn. that had higher values than the other two sites. Dissolved oxygen concentrations and pH were similar among the three sites on the Red River. The patterns in turbidity were mostly related to streamflow conditions and were similar among the three sites on the Red River of the North.

Minnesota, North Dakota↗