Search USGS⌕ Search

Geology topics

Carol A. Finn

Publications and source records attributed to Carol A. Finn.

At least 19 recordsLinked to original sources
USGS2025

Multi-scale geophysical imaging of a hydrothermal system in Yellowstone National Park, USA

Little is known about the local plumbing systems that fuel Yellowstone’s famous hot springs, geysers and mud pots. A multi-method, multi-scale geophysical investigation was carried out in the Obsidian Pool Thermal Area (OPTA) to: (i) delineate the lateral extent of the hydrothermal area and associated surface features; (ii) estimate the dimensions of the upflow zone and identify its main controlling structures; (iii) assess fluids circulation pathways from depth to surface. Ground and airborne geophysical data were acquired to connect local and regional scales, from shallow to large depths. Maps of surface electrical resistivity show a strong correlation with hydrothermal features. At in-termediate depths, electrical resistivity permits delineating the upper limit of the upflow zone, while Poisson’s ratio highlights differences in subsurface fluid content. Combining these results with surface observations and topographic information, we speculate that differential mixing of hydrothermal and fresh water could explain the wide diversity of features observed at OPTA. Low electrical resistivity observed at large depths also suggest that a vast upflow zone, controlled by rhyolite flows and conjugate faults, underlies the OPTA. We speculate that hydrothermal fluids rise along fractures and reach the surface in topographic lows to form hydrothermal features. Our results show that synoptic, multi-scale geophysical measurements provide a roadmap for understanding where and how geologic heterogeneity, topography, fluid-gas separation, and the mixing of thermal and meteoric waters conspire to produce the wide variety of Yellowstone’s renowned hydrothermal features.

Wyoming↗
USGS2025

State shifts in the deep Critical Zone drive landscape evolution in volcanic terrains

Understanding the near-surface environment where atmospheric and solid earth processes interact, often termed the “Critical Zone,” is important for assessing resources and building resilient societies. Here, we examine a volcanic landscape in the Oregon Cascade Range, an understudied Critical Zone setting that is host to major regional water resources, pervasive silicate weathering, and significant geohazards. We leverage a bedrock age chronosequence to show that the volcanic Critical Zone undergoes a structural shift, from depth extents of >1 km to meters, over timescales of ~1 My. We map an active groundwater volume comparable to major continental lakes, stored at the Cascade Range crest. This state shift makes volcanic landscape evolution a unique probe of deep coupling between Earth systems.

PNAS↗
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↗
USGS2024

Deep magmatic staging chambers for crustal layered mafic intrusions: An example from the Bushveld Complex of southern Africa

The deep mafic magmatic staging chambers of layered mafic intrusions have been conjectured but not imaged. Their existence has long been postulated from geochemical models which require multiple magma injections from staging chambers to account for their multi-scale igneous layering and variations in sources and degrees of crustal contamination. For the Bushveld Complex of southern Africa, the world’s largest layered mafic intrusion, seismic receiver functions identify a diffuse crust-mantle transition beneath the Complex that suggests high velocity lower crust and/or uppermost lithospheric mantle. Here we present 3D gravity modelling of the Bushveld Complex that includes dense material at the crust-mantle boundary, imaging for the first time, remnants of magma staging chambers. They underlie the whole Bushveld Complex and extend westwards to the Molopo Farms Complex in Botswana. Feeders to the Bushveld Complex coincide with intersections of major faults like the Thabazimbi-Murchison Lineament and Sugarbush Fault with the staging chamber. This identification of magmatic staging chambers beneath the Bushveld relates to similar geophysically imaged lower crustal features beneath the Duluth and Stillwater Complexes. Comparison of seismic and gravity data with geochemical models from other complexes aid in development of models for the magmatic architecture of layered mafic intrusions in general.

Precambrian Research↗
USGS2023

Electromagnetic and magnetic imaging of the Stillwater Complex, Montana, USA

Modelling and analysis of helicopter electromagnetic data result in resistivity and susceptibility models and derivatives of magnetic data that characterise shallow parts of the Stillwater Complex, critical for aiding exploration and expansion of globally scarce critical and battery mineral resources that include platinum group elements, nickel, copper and chromium. The magnetic susceptibly models derived from the electromagnetic data and the tilt derivative of the magnetic data image layering, mafic dikes, banded iron formation, and serpentinised peridotite. Known areas with contact-type mineralisation are generally characterised by low resistivities and susceptibilities where the volume of mineralised rock is large and/or the depth is shallow. We use iso-cluster and edge detection analysis of both resistivities and susceptibilities to identify potential mineralisation in poorly characterised regions as well as faults. Low resistivity layers beneath large landslides reflect water saturated porous slip surfaces which can interfere with drilling. This uncommon approach of tightly linking the resistivity and susceptibility models and magnetic anomaly data to rock property, surficial geologic, drill hole and soil geochemistry data to image the geology in the upper ∼100 m, aids identification of prospective mineralised regions as well landslides and faults that can impact mineral exploration and local hazards.

Montana↗
USGS2022

The 180-km-long Meers-Willow Fault System in the Southern Oklahoma Aulacogen: A potential U.S. mid-continent seismic hazard

We integrate new high-resolution aeromagnetic data with seismic reflection data, well logs, satellite remote sensing, and field observations to provide a regional view of buried and exposed structures in the Southern Oklahoma Aulacogen and to assess their potential for future seismicity. Trends ranging from NW−SE to ∼E−W, peaking at 330° ± 4.5° and 280° ± 3°, dominate the magnetic lineaments of the Southern Oklahoma Aulacogen, reflecting basement contacts, dikes, and faults, including a previously unmapped ∼100-km-long basement fault, which is herein referred to as the Willow fault. The fault disrupts, truncates, and vertically offsets basement-related seismic reflectors and overlying Paleozoic strata up through the Permian reflectors. Surface deformation along the trend includes fault-parallel monoclinal folds, pervasive fractures, and fracture-hosted mud dikes in Permian evaporite units. These structures indicate a Permian or post-Permian reactivation of the fault. Along-strike, the Willow fault connects to the NW-trending, seismically active Meers Fault to comprise the ∼180-km-long Meers-Willow fault system, which potentially represents a major seismic hazard along the Southern Oklahoma Aulacogen. Fault slip potential analyses of the mapped potential fault traces show that seismic hazards are elevated where faults have steeper dips. Given some uncertainty in the regional stress state, we also show that hazards along the NW−SE to E−W trending faults vary considerably within the uncertainty range. We propose that the Meers-Willow fault system originated as a Cambrian aulacogen-scale, basement-rooted fault that was later reactivated as a left-lateral strike-slip fault (with ∼40 km displacement) during the late Paleozoic Ancestral Rocky Mountain orogeny, highlighting that lateral offset accommodated a major component of deformation during the orogen.

Oklahoma↗
USGS2022

Geophysical imaging of the Yellowstone hydrothermal plumbing system

The nature of Yellowstone National Park’s plumbing system linking deep thermal fluids to its legendary thermal features is virtually unknown. The prevailing concepts of Yellowstone hydrology and chemistry are that fluids reside in reservoirs with unknown geometries, flow laterally from distal sources and emerge at the edges of lava flows . Here we present a high-resolution synoptic view of pathways of the Yellowstone hydrothermal system derived from electrical resistivity and magnetic susceptibility models of airborne geophysical data . Groundwater and thermal fluids containing appreciable total dissolved solids significantly reduce resistivities of porous volcanic rocks and are differentiated by their resistivity signatures . Clay sequences mapped in thermal areas and boreholes typically form at depths of less than 1,000 metres over fault-controlled thermal fluid and/or gas conduits . We show that most thermal features are located above high-flux conduits along buried faults capped with clay that has low resistivity and low susceptibility. Shallow subhorizontal pathways feed groundwater into basins that mixes with thermal fluids from vertical conduits. These mixed fluids emerge at the surface, controlled by surficial permeability, and flow outwards along deeper brecciated layers. These outflows, continuing between the geyser basins, mix with local groundwater and thermal fluids to produce the observed geochemical signatures. Our high-fidelity images inform geochemical and groundwater models for hydrothermal systems worldwide.

Idaho, Wyoming↗
USGS2021

When Wyoming became Superior: Oblique convergence along the southern Trans-Hudson orogen

The Trans-Hudson orogen (THO) is one of the best-preserved Proterozoic orogens on Earth, largely unaffected by subsequent tectonism, yet its southern extent lies concealed beneath the North American Central Plains. A new 3D resistivity model over the southern orogen is developed and interpreted alongside borehole, potential field, and seismic reflection data. We present the first synoptic crustal view of the southern THO and a new tectonic model of the orogen. Our model reveals high-conductivity belts marking paleo-subduction zones while the orogen center consists of deeply exhumed relatively dense and mostly magnetic juvenile crust preserved between the deformed margins of the Wyoming and Superior cratons. Complex structure along the western margin suggests convergence began with oblique subduction and the northward transport of severed fragments of the Wyoming Province. High-conductivity belts are in places offset from upper-crustal geophysical boundaries, consistent with the thrusting of accreted arcs over the Archean margins during terminal closure.

Geophysical Research Letters↗
USGS2021

Airborne geophysical imaging of weak zones on Iliamna Volcano, Alaska: Implications for slope stability

Water‐saturated, hydrothermally altered rocks reduce the strength of volcanic edifices and increase the potential for sector collapses and far‐traveled mass flows of unconsolidated debris. Iliamna Volcano is an andesitic stratovolcano located on the western side of the Cook Inlet, ∼225 km southwest of Anchorage and is a source of repeated avalanches. The widespread snow and ice cover on Iliamna Volcano make surface alteration difficult to identify. However, intense hydrothermal alteration significantly reduces both the electrical resistivity and magnetization of volcanic rock and can therefore be identified with airborne geophysical measurements. We use airborne electromagnetic and magnetic data to map snow and ice thickness and identify underlying alteration zones at Iliamna Volcano, Alaska. Resistivities were calculated to an average depth of >300 m, and a 3‐D susceptibility model extends from the surface to the base of the volcano, about 3,000 m below the summit. Geophysical models image low resistivity (<30 ohm‐m) and low susceptibilities near the summit of Iliamna and below its older vent complex, with the low susceptibilities indicating alteration up to ∼800 m in thickness. Thin conductors (∼50–100 m thick) on the edifice slopes coincide with recorded locations of repeated debris avalanches over the past ∼60 years and are attributed to saturated zones at high elevation. Three‐dimensional slope stability models based upon the geophysically constrained alteration distribution suggest the edifice of Iliamna is unstable and could lead to collapse scars ∼400 m deep near the current and former vent complexes.

Alaska↗
USGS2021

Geometry of the Bushveld Complex from 3D potential field modelling

A full three-dimensional (3D) potential field model of the central and southern Bushveld Complex reveals information about the Complex in areas obscured by younger geological cover. Previously, two-dimensional gravity models and a few magnetic models limited to certain sections of the Bushveld Complex have been used to propose geometries for the Rustenburg Layered Suite, especially in the western and eastern lobes. These models were often used to support different emplacement models. Although these models provided valuable information, two-and-a-half-dimensional (2.5D) potential field modelling is not well suited to modelling complex 3D geology. Also, in most cases, only the magnetic or gravity data were modelled, but jointly modelling both data sets better constrains the results, as was shown recently for a 3D model of the northern lobe. Joint 3D modelling of regional gravity and magnetic data combined with published crustal thickness models derived from broadband seismic tomography studies and constrained by density and susceptibility data, geologic mapping, boreholes and seismic reflection data were used to create a 3D model of the central and southeastern sections of the Bushveld Complex, as well as the southern part of the northern lobe. The model shows a complex geometry with thick continuous Rustenburg Layered Suite S in most of the western and southeastern lobes, but less continuous Rustenburg Layered Suite in the eastern lobe. Large domes or thick granites and granophyre in the latter interrupt the continuity of the Rustenburg Layered Suite and the western and eastern lobes are strictly speaking only partially connected in places. However, they are not separate intrusions, but one disconnected by pre-existing and synmagmatic updoming. Three possible feeders were modelled in the northern, western, and south-eastern lobes.

Precambrian Research↗
USGS2021

3-D Modeling of the Duluth Complex from geophysical data

The Mesoproterozoic Duluth Complex in northeastern Minnesota is one of the major plutonic components of the Midcontinent Rift System and hosts a variety of copper-nickel sulfide and platinum-group element deposits. The Duluth Complex is composed of a series of individual mafic and felsic intrusions emplaced 1110-1098 Ma within Paleoproterozoic sedimentary rocks of the Animikie basin and volcanic flows of the Midcontinent Rift. Prior work has included 2-D modeling and qualitative geologic interpretations of gravity and magnetic data (e.g., Chandler, 1990; Chandler and Ferderer, 1989), much of which is still preliminary (V. Chandler, written commun., 2020). Three-dimensional modeling has been limited, with only one 3-D model created using Bouguer gravity data constrained by seismic-reflection interpretations as part of a PhD thesis (Allen, 1994). Given the complex geology of the area, 3-D modeling is useful for providing a complete picture of the variable densities, susceptibilities, and electrical resistivities throughout the Duluth Complex and associated volcanic rocks as well as their depth extent beneath sedimentary cover. Models of these geophysical properties at depth enable more accurate geologic mapping in the subsurface which can lead to an improved understanding of the formation history of the Duluth Complex. In this study, we use aeromagnetic data acquired between 1979-1991 (Chandler, 2007), Bouguer gravity data collected since 1950 (Chandler and Lively, 2019), and magnetotelluric data collected in 2019 to create new 2-D and 3-D geophysical models of the Duluth Complex constrained by seismic reflection, geologic, and rock property data. An inversion of the Bouguer gravity data for thickness of the Duluth Complex using constant densities of 3110 kg/m3 and 2670 kg/m3 for the Duluth Complex and surrounding crustal rocks, respectively, results in thicknesses ranging from ~3-28 km for the Duluth Complex and related intrusions and volcanic rocks (Figure 1A). A 3-D model of the magnetotelluric data reveals low resistivity anomalies at ~5-10 km depth below the northern margin of the Duluth Complex and below the Greenwood Lake intrusion (Figure 1B). We expect to encounter low resistivities at depth associated with the Paleoproterozoic Animikie basin, which makes up the floor of the Duluth Complex, and therefore interpret these anomalies as either the base of the complex or as fragments of Animikie sediments interfingered with igneous intrusive rocks. Finally, 3-D voxel models of density and susceptibility illuminate the subsurface distribution of rock properties below the Duluth Complex which, in combination with resistivity and thickness models, can be used to create a 3-D geologic map of this area.

Minnesota↗
USGS2020

Geological and thermal control of the hydrothermal system in northern Yellowstone Lake: Inferences from high resolution magnetic surveys

A multiscale magnetic survey of the northern basin of Yellowstone Lake was undertaken in 2016 as part of the Hydrothermal Dynamics of Yellowstone Lake Project (HD‐YLAKE)—a broad research effort to characterize the cause‐and‐effect relationships between geologic and environmental processes and hydrothermal activity on the lake floor. The magnetic survey includes lake surface, regional aeromagnetic, and near‐bottom autonomous underwater vehicle (AUV) data. The study reveals a strong contrast between the northeastern lake basin, characterized by a regional magnetic low punctuated by stronger local magnetic lows, many of which host hydrothermal vent activity, and the northwestern lake basin with higher‐amplitude magnetic anomalies and no obvious hydrothermal activity or punctuated magnetic lows. The boundary between these two regions is marked by a steep gradient in heat flow and magnetic values, likely reflecting a significant structure within the currently active ~20‐km‐long Eagle Bay‐Lake Hotel fault zone that may be related to the ~2.08‐Ma Huckleberry Ridge caldera rim. Modeling suggests that the broad northeastern magnetic low reflects both a shallower Curie isotherm and widespread hydrothermal activity that has demagnetized the rock. Along the western lake shoreline are sinuous‐shaped, high‐amplitude magnetic anomaly highs, interpreted as lava flow fronts of upper units of the West Thumb rhyolite. The AUV magnetic survey shows decreased magnetization at the periphery of the active Deep Hole hydrothermal vent. We postulate that lower magnetization in the outer zone results from enhanced hydrothermal alteration of rhyolite by hydrothermal condensates while the vapor‐dominated center of the vent is less altered.

Wyoming↗
USGS2020

Mapping the 3-D extent of the Stillwater Complex, Montana—Implications for potential platinum group element exploration and development

Geophysical models characterize the exposed and interpreted buried extent of the Stillwater Complex, critical for understanding the origin of the layered mafic intrusion and its associated high-grade platinum group element resources. The 3D models, constrained by gravity, magnetic, xenolith, seismic, borehole, and rock property data indicate that the likely maximum extent of the Stillwater Complex beneath Phanerozoic cover is ~10 times greater than its outcrop, ~2240 km 2 . The thickness values are poorly constrained but vary from ~7000 to 12,000 m, depending on crustal and mantle density variations and depths to the top of the lower crust and mantle. This thickness may include dense metasedimentary units of the basin into which the Stillwater Complex intruded. Using the modeled thickness results in a volume estimate of ~24,700 km 3 , albeit poorly constrained. New analyses of xenoliths from the Cretaceous Sliderock and Suzie Peak intrusions produce ages of 2706–2716 Ma, corresponding to the age of the Stillwater Complex, and 2813 Ma, corresponding to the age of Archean gneissic basement. Seismic reflectors in inferred Archean crystalline basement, possibly including the Stillwater Complex, dip ~25–30° north, with segments dipping as much as 70° north. Layered reflectors beneath the Phanerozoic sedimentary section and above the inferred Archean crystalline basement may represent metasedimentary units, perhaps a southern extension of the Mesoproterozoic Belt Basin. The potential field models and the seismic reflection data suggest that the Stillwater Complex was dipping northward prior to deposition of Cambrian strata, perhaps uplifted in the late Archean or Proterozoic as previously proposed, and that during Laramide times, faulting and intrusions highly disrupted the complex. Temperature measurements from boreholes help constrain the depths of feasible mining of the complex.

Montana↗
USGS2018

On the feasibility of real-time mapping of the geoelectric field across North America

A review is given of the present feasibility for accurately mapping geoelectric fields across North America in near-realtime by modeling geomagnetic monitoring and magnetotelluric survey data. Should this capability be successfully developed, it could inform utility companies of magnetic-storm interference on electric-power-grid systems. That real-time mapping of geoelectric fields is a challenge is reflective of (1) the spatiotemporal complexity of geomagnetic variation, especially during magnetic storms, (2) the sparse distribution of ground-based geomagnetic monitoring stations that report data in realtime, (3) the spatial complexity of three-dimensional solid-Earth impedance, and (4) the geographically incomplete state of continental-scale magnetotelluric surveys.

Open-File Report↗
USGS2018

Combining multiphase groundwater flow and slope stability models to assess stratovolcano flank collapse in the Cascade Range

Hydrothermal alteration can create low‐permeability zones, potentially resulting in elevated pore‐fluid pressures, within a volcanic edifice. Strength reduction by rock alteration and high pore‐fluid pressures have been suggested as a mechanism for edifice flank instability. Here we combine numerical models of multiphase heat transport and groundwater flow with a slope‐stability code that incorporates three‐dimensional distributions of strength and pore‐water pressure to address the following questions: (1) What permeability distributions and contrasts produce elevated pore‐fluid pressures in a stratovolcano? (2) What are the effects of these elevated pressures on flank stability? (3) Finally, what are the effects of magma intrusion on potential flank failure in an edifice? Simulation results show that under a range of plausible parameters, water tables in a stratovolcano can be elevated or perched. These elevated water tables result in universally lower stability (lower factor of safety) compared with equivalent dry edifices, indicating a higher likelihood of flank collapse. Low‐permeability (<1 × 10 −17 m 2 ) layers such as altered pyroclastic deposits or breccias can result in locally saturated regions (perched water) and lower factors of safety near the ground surface but may actually reduce liquid water saturation and pore pressures in the core of the edifice and thus may favor small, shallow collapses over larger, deeper collapses. Magma intrusion into the base of the edifice increases pore‐fluid pressures and decreases the factor of safety. However, the shear strength of edifice rocks also exerts a significant control on stability, so both mechanical properties and pore‐fluid pressures are important for stability assessments.

Washington↗
USGS2018

Three-dimensional geophysical mapping of shallow water saturated altered rocks at Mount Baker, Washington: Implications for slope stability

Water-saturated hydrothermal alteration reduces the strength of volcanic edifices, increasing the potential for catastrophic sector collapses that can lead to far traveled and destructive debris flows. Intense hydrothermal alteration significantly lowers the resistivity and magnetization of volcanic rock and therefore hydrothermally altered rocks can be identified with helicopter electromagnetic and magnetic measurements. Geophysical models constrained by rock properties and geologic mapping show that intensely altered rock is restricted to two small (500 m diameter), >150 m thick regions around Sherman Crater and Dorr Fumarole Field at Mount Baker, Washington. This distribution of alteration contrasts with much thicker and widespread alteration encompassing the summits of Mounts Adams and Rainier prior to the 5600 year old Osceola collapse, which is most likely due to extreme erosion and the limited duration of summit magmatism at Mount Baker. In addition, the models suggest that the upper ~300 m of rock contains water which could help to lubricate potential debris flows. Slope stability modeling incorporating the geophysically modeled distribution of alteration and water indicates that the most likely and largest (~0.1 km 3 ) collapses are from the east side of Sherman Crater. Alteration at Dorr Fumarole Field raises the collapse hazard there, but not significantly because of its lower slope angles. Geochemistry and analogs from other volcanoes suggest a model for the edifice hydrothermal system.

Washington↗