Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Geophysical Research”

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.

At least 883 records · Page 49Linked to original sources
USGS2012

Nitrate removal in deep sediments of a nitrogen-rich river network: A test of a conceptual model

Many estimates of nitrogen removal in streams and watersheds do not include or account for nitrate removal in deep sediments, particularly in gaining streams. We developed and tested a conceptual model for nitrate removal in deep sediments in a nitrogen-rich river network. The model predicts that oxic, nitrate-rich groundwater will become depleted in nitrate as groundwater upwelling through sediments encounters a zone that contains buried particulate organic carbon, which promotes redox conditions favorable for nitrate removal. We tested the model at eight sites in upwelling reaches of lotic ecosystems in the Waupaca River Watershed that varied by three orders of magnitude in groundwater nitrate concentration. We measured denitrification potential in sediment core sections to 30 cm and developed vertical nitrate profiles to a depth of about 1 m with peepers and piezometer nests. Denitrification potential was higher, on average, in shallower core sections. However, core sections deeper than 5 cm accounted for 70%, on average, of the depth-integrated denitrification potential. Denitrification potential increased linearly with groundwater nitrate concentration up to 2 mg NO 3 -N/L but the relationship broke down at higher concentrations (> 5 mg NO 3 -N/L), a pattern that suggests nitrate saturation. At most sites groundwater nitrate declined from high concentrations at depth to much lower concentrations prior to discharge into the surface water. The profiles suggested that nitrate removal occurred at sediment depths between 20 and 40 cm. Dissolved oxygen concentrations were much higher in deep sediments than in pore water at 5 cm sediment depth at most locations. The substantial denitrification potential in deep sediments coupled with the declines in nitrate and dissolved oxygen concentrations in upwelling groundwater suggest that our conceptual model for nitrate removal in deep sediments is applicable to this river network. Our results suggest that nitrate removal rates can be high in deep sediments of upwelling stream reaches, which may have implications for efforts to understand and quantify nitrogen transport and removal at larger scales.

Journal of Geophysical Research: Biogeosciences↗
USGS2012

Climatic limits on foliar growth during major droughts in the Southwestern U.S.A.

Pronounced droughts during the 1950s and 2000s in the Southwestern U.S.A. (SW) provide an opportunity to compare mesoscale ecosystem responses to anomalously dry conditions before and during the regional warming that started in the late 1970s. This year-round warming has produced fewer cool season freezes, losses in regional snowpack, an 8-10 day advance in spring onset, and hotter summers, all of which should affect vegetation differently across seasons and elevations. Here, we examine indices that represent climatic limits on foliar growth for both drought periods, and evaluate these indices for areas that experienced tree mortality during the 2000s drought. Relative to the 1950s drought, warmer conditions during the 2000s drought decreased the occurrence of temperatures too low for foliar growth at lower elevations in winter and higher elevations in summer. Higher vapor pressure deficits (VPDs) largely driven by warmer temperatures in the more recent drought were more limiting to foliar growth from spring through summer at lower and middle elevations. At many locations where tree mortality occurred during the 2000s drought, low-temperature constraints on foliar growth were extremely unlimiting, whereas VPD constraints were extremely limiting from early spring through late autumn. Our analysis shows that in physiographically complex regions like the SW, seasonality and elevational gradients are important for understanding vegetative responses to warming. It also suggests that continued warming will increase the degree to which VPD limits foliar growth during future droughts, and expand its reach to higher elevations and other seasons.

Journal of Geophysical Research G: Biogeosciences↗
USGS2012

Quantifying soil surface change in degraded drylands: shrub encroachment and effects of fire and vegetation removal in a desert grassland

Woody plant encroachment, a worldwide phenomenon, is a major driver of land degradation in desert grasslands. Woody plant encroachment by shrub functional types ultimately leads to the formation of a patchy landscape with fertile shrub patches interspaced with nutrient-depleted bare soil patches. This is considered to be an irreversible process of land and soil degradation. Recent studies have indicated that in the early stages of shrub encroachment, when there is sufficient herbaceous connectivity, fires (prescribed or natural) might provide some reversibility to the shrub encroachment process by negatively affecting shrub demography and homogenizing soil resources across patches within weeks to months after burning. A comprehensive understanding of longer term changes in microtopography and spatial patterning of soil properties following fire in shrub-encroached grasslands is desirable. Here, we investigate the changes in microtopography with LiDAR (light detection and ranging), vegetation recovery, and spatial pattering of soil properties in replicated burned, clipped, and control areas in a shrub-grass transition zone in the northern Chihuahuan Desert four years after prescribed fire or clipping. Results indicate a greater homogeneity in soil, microtopography, and vegetation patterning on burned relative to clipped and control treatments. Findings provide further evidence that disturbance by prescribed fire may allow for reversal of the shrub encroachment process, if the event occurs in the early stages of the vegetation shift. Improved understanding of longer-term effects of fire and associated changes in soil patterning can inform the use and role of fire in the context of changing disturbance regimes and climate.

Chihuahuan Desert↗
USGS2018

Three-dimensional modeling of fine sediment transport by waves and currents in a shallow estuary

A suspended sediment transport model is implemented in the unstructured‐grid SUNTANS model and applied to study fine‐grained sediment transport in South San Francisco Bay. The model enables calculation of suspension of bottom sediment based on combined forcing of tidal currents and wind waves. We show that accurate results can be obtained by employing two‐size classes which are representative of microflocs and macroflocs in the Bay. A key finding of the paper is that the critical calibration parameter is the ratio of the erosion of the microflocs to macroflocs from the bed. Different values of this erosion ratio are needed on the shallow shoals and deeper channels because of the different nature of the sediment dynamics in these regions. Application of a spatially variable erosion ratio and critical shear stress for erosion is shown to accurately reproduce observed suspended sediment concentration at four‐field sites located along a cross‐channel transect. The results reveal a stark contrast between the behavior of the suspended sediment concentration on the shoals and in the deep channel. Waves are shown to resuspend sediments on the shoals, although tidal and wind‐generated currents are needed to mix the thin wave‐driven suspensions into the water column. The contribution to the suspended sediment concentration in the channel by transport from the shoals is similar in magnitude to that due to local resuspension. However, the local contribution is in phase with strong bottom currents which resuspend the sediments, while the contribution from the shoals peaks during low‐water slack tide.

San Francisco Bay↗
USGS2018

The thermophysical properties of the Bagnold Dunes, Mars: Ground truthing orbital data

We compare the thermophysical properties and particle sizes derived from the Mars Science Laboratory rover's Ground Temperature Sensor of the Bagnold dunes, specifically Namib dune, to those derived orbitally from Thermal Emission Imaging System, ultimately linking these measurements to ground truth particle sizes determined from Mars Hand Lens Imager images. In general, we find that all three datasets report consistent particle sizes for the Bagnold dunes (~110–350 μm and are within measurement and model uncertainties), indicating that particle sizes of homogeneous materials inferred from temperature measurements and thermophysical models are reliable. Furthermore, we examine the effects of two physical characteristics that could influence the modeled thermal inertia and particle sizes, including (1) fine‐scale (centimeter to meter scale) ripples and (2) thin layering of indurated/armored materials. To first order, we find that small‐scale ripples and thin (approximately centimeter scale) layers do not significantly affect the determination of bulk thermal inertia from orbital thermal data using a single nighttime temperature. Modeling of a layer of coarse or indurated material reveals that a thin layer (< ~5 mm; similar to what was observed by the Curiosity rover) would not significantly change the observed thermal properties of the surface and would be dominated by the properties of the underlying material. Thermal inertia and particle sizes of relatively homogeneous materials derived from nighttime orbital data should be considered as reliable, as long as there are no significant subpixel anisothermality effects (e.g., lateral mixing of multiple thermophysically distinct materials).

Journal of Geophysical Research: Planets↗
USGS2018

Consequences of abrading bed load on vertical and lateral bedrock erosion in a curved experimental channel

In this study, we conducted multiple physical experiments to estimate the efficacy and spatial pattern of erosion by abrading sediment moving through a simple U‐shaped channel bend with erodible bed and banks. The experiments showed that in the bend, lateral abrasion followed a monotonically increasing linear relationship with sediment feed rate. However, vertical incision had a more complex relation with the sediment feed rate, with an initial increase in abrasion as the feed rate increased followed by a decrease in abrasion of the bed as cover effects became dominant at higher feed rates. Bank erosion was large in places where the width and the lateral slope of the point bar were relatively large. On the other hand, in places where the width of the point bar was smaller, the bedrock bed was eroded primarily along the boundary of the point bar, resulting in a bedrock bench near the outer bank.

Shimanto River↗
USGS2018

A model ensemble for projecting multi‐decadal coastal cliff retreat during the 21st century

Sea cliff retreat rates are expected to accelerate with rising sea levels during the 21 st century. Here we develop an approach for a multi‐model ensemble that efficiently projects time‐averaged sea cliff retreat over multi‐decadal time scales and large (>50 km) spatial scales. The ensemble consists of five simple 1‐D models adapted from the literature that relate sea cliff retreat to wave impacts, sea level rise (SLR), historical cliff behavior, and cross‐shore profile geometry. Ensemble predictions are based on Monte Carlo simulations of each individual model, which account for the uncertainty of model parameters. The consensus of the individual models also weights uncertainty, such that uncertainty is greater when predictions from different models do not agree. A calibrated, but unvalidated, ensemble was applied to the 475 km‐long coastline of Southern California (USA), with 4 SLR scenarios of 0.5, 0.93, 1.5, and 2 m by 2100. Results suggest that future retreat rates could increase relative to mean historical rates by more than two‐fold for the higher SLR scenarios, causing an average total land loss of 19 – 41 m by 2100. However, model uncertainty ranges from +/‐ 5 – 15 m, reflecting the inherent difficulties of projecting cliff retreat over multiple decades. To enhance ensemble performance, future work could include weighting each model by its skill in matching observations in different morphological settings

Journal of Geophysical Research F: Earth Surface↗
USGS2018

Enhancement of a parsimonious water balance model to simulate surface hydrology in a glacierized watershed

The U.S. Geological Survey monthly water balance model (MWBM) was enhanced with the capability to simulate glaciers in order to make it more suitable for simulating cold region hydrology. The new model, MWBMglacier, is demonstrated in the heavily glacierized and ecologically important Copper River watershed in Southcentral Alaska. Simulated water budget components compared well to satellite‐based observations and ground measurements of streamflow, evapotranspiration, snow extent, and total water storage, with differences ranging from 0.2% to 7% of the precipitation flux. Nash Sutcliffe efficiency for simulated and observed streamflow was greater than 0.8 for six of eight stream gages. Snow extent matched satellite‐based observations with Nash Sutcliffe efficiency values of greater than 0.89 in the four Copper River ecoregions represented. During the simulation period 1949 to 2009, glacier ice melt contributed 25% of total runoff, ranging from 12% to 45% in different tributaries, and glacierized area was reduced by 6%. Statistically significant ( p < 0.05) decreasing and increasing trends in annual glacier mass balance occurred during the multidecade cool and warm phases of the Pacific Decadal Oscillation, respectively, reinforcing the link between climate perturbations and glacier mass balance change. The simulations of glaciers and total runoff for a large, remote region of Alaska provide useful data to evaluate hydrologic, cryospheric, ecologic, and climatic trends. MWBM glacier is a valuable tool to understand when, and to what extent, streamflow may increase or decrease as glaciers respond to a changing climate.

Journal of Geophysical Research F: Earth Surface↗
USGS2018

A global empirical model for near real-time assessment of seismically induced landslides

Earthquake-triggered landslides are a significant hazard in seismically active regions, but our ability to assess the hazard they pose in near real-time is limited. In this study, we present a new globally applicable model for seismically induced landslides based on the most comprehensive global dataset available; we use 23 landslide inventories that span a range of earthquake magnitudes and climatic and tectonic settings. We use logistic regression to relate the presence and distribution of earthquake-triggered landslides with spatially distributed estimates of ground shaking, topographic slope, lithology, land-cover type, and a topographic index designed to estimate variability in soil wetness to provide an empirical model of landslide distribution. We tested over 100 combinations of independent predictor variables to find the best-fitting model, using a diverse set of statistical tests. Blind validation tests show the model accurately estimates the distribution of available landslide inventories. The results indicate that the model is reliable and stable, with high “balanced accuracy” (correctly vs. incorrectly classified pixels) for the majority of test events. A cross validation analysis shows high balanced accuracy for a majority of events as well. By combining near-real time estimates of ground shaking with globally available landslide susceptibility data, this model provides a tool to estimate the distribution of co-seismic landslide hazard within minutes of the occurrence of any earthquake worldwide for which a USGS ShakeMap is available.

Journal of Geophysical Research: Earth Surface↗
USGS2018

Long-term evolution of sand transport through a river network: Relative influences of a dam versus natural changes in grain size from sand waves

Temporal and spatial nonuniformity in supplies of water and sand in a river network leads to sand transport that is in local disequilibrium with the upstream sand supply. In such river networks, sand is transported downstream as elongating waves in which coupled changes in grain size and transport occur. Depending on the magnitude of each sand‐supplying event and the interval between such events, changes in bed‐sand grain size associated with sand‐wave passage may more strongly regulate sand transport than do changes in water discharge. When sand transport is controlled more by episodic resupply of sand than by discharge, upstream dam construction may exacerbate or mitigate sand‐transport disequilibria, thus leading to complicated and difficult‐to‐predict patterns of deposition and erosion. We analyzed all historical sediment‐transport data and embarked on a 4‐year program of continuous sediment‐transport measurements to describe disequilibrium sand transport in a river network. Results indicate that sand transport in long river segments can evolve over ≥50‐year timescales following rare large sand‐supplying events. These natural changes in sand transport in distal downstream river segments can be larger than those caused by an upstream dam. Because there is no way to know a priori whether sand transport in a river has changed in response to changes in the upstream sand supply, contemporary continuous measurements of sand transport are required for accurate sand loads and budgeting. Analysis of only historical sediment‐transport measurements, as is common in the literature, may lead to incorrect conclusions with respect to current or future sediment‐transport conditions.

Journal of Geophysical Research: Earth Surface↗
USGS2018

Evolving environmental and geometric controls on Columbia Glacier’s continued retreat

Geometry strongly controls the dynamic behavior of marine‐terminating (tidewater) glaciers, significantly influencing advance and retreat cycles independent of climate. Yet the recent, nearly ubiquitous retreat of tidewater glaciers suggests that changes in atmospheric and oceanic forcing may also drive dynamic change. To isolate the influence of geometry on tidewater glacier dynamics, we analyzed detailed observational time series from 2012 to 2016 for two tidewater glaciers with shared dynamic histories and environmental forcing: Columbia Glacier and its former tributary (Post Glacier) in southcentral Alaska. We find that although terminus retreat has driven decadal‐scale changes in dynamics of the Columbia‐Post system, environmental factors contribute to short‐term (i.e., seasonal) dynamic variability. In particular, analysis of force balance time series indicates that observed variations in speed result from seasonal changes to the subglacial hydrologic system and associated changes in basal drag. Variations in terminus position only drive noticeable speed change when the terminus retreats from regions of relatively high basal drag. In agreement with long‐term analyses of Columbia Glacier, we find that terminus geometry can perturb the timing of seasonal ice flow patterns. Specifically, our data support the idea that retreat of a glacier terminus into deeper water is accompanied by a shift in the primary control on frontal ablation. Although our analysis focuses on two Alaskan glaciers, our data suggest that changes in the relative importance of surface meltwater and buoyancy effects on submarine melting and/or calving may manifest as a shift in terminus change seasonality and offer a mechanism to identify frontal ablation controls.

Alaska↗
USGS2018

Atmospheric and surface climate associated with 1986–2013 wildfires in North America

We analyze climate simulations conducted with the RegCM3 regional climate model on 50‐ and 15‐km model grids to diagnose the dependence of wildfire incidence and area burned variations on monthly climate long‐term means and anomalies over North America for the period 1986–2013. We created a new wildfire database by merging the Fire Program Analysis Fire‐Occurrence Database, the National Interagency Fire Center Fire History Data, and the Canadian National Fire Database. The database includes 2,083,865 daily fire starts that burned a total of 1.25 × 10 8 ha in North America. We derive long‐term climatologies, standardized gamma indices, and composite climate anomalies of atmospheric circulation (500‐hPa height and wind) and various surface fields (e.g., solar radiation, soil moisture, vapor pressure deficit, and latent and sensible heat fluxes) to illustrate the climatology of burned area. The immediate and lagged monthly atmospheric circulation and surface climate anomalies differentiate high‐ and low‐fire years and the role of El Niño–Southern Oscillation in wildfire occurrence. Our approach demonstrates the association of the seasonal cycles of wildfire and climate and the strong role of climatic variability in modulating the seasonal cycle as a control of wildfire on monthly time scales.

Journal of Geophysical Research: Biogeosciences↗
USGS2018

Fire, flood, and drought: Extreme climate events alter flow paths and stream chemistry

Extreme climate events—such as hurricanes, droughts, extreme precipitation, and wildfires—have the potential to alter watershed processes and stream response. Yet due to the destructive and hazardous nature and unpredictability of such events, capturing their hydrochemical signal is challenging. A 5‐year postwildfire study of stream chemistry in the Fourmile Creek watershed, Colorado Front Range, USA, focused on high‐frequency storm sampling. During the study, the watershed was impacted by three additional extreme climate events—drought and two periods of extreme rainfall totals. These events altered concentration‐discharge relationships in ways that elucidate how hydrologic flow paths and source material availability affect stream water chemistry. Reduced infiltration after wildfire led to overland flow during thunderstorms, which conveyed ash and soil into streams. This resulted in elevated stream concentrations of constituents elevated in ash—Ca, K, Mg, alkalinity, and dissolved organic carbon—along with sediment and nitrate. Subsurface flow paths were bypassed, leading to low concentrations of Na and SiO 2 , which are bedrock derived and not elevated in ash. During drought conditions, when stream discharge was <20% of average, concentrations of sediment, dissolved organic carbon, and Ca fell below average concentrations, but SiO 2 did not. Extreme rainfall totals saturated the subsurface and led to prolonged elevated stream discharge. Concentration‐discharge relationships for bedrock‐derived constituents, such as Ca and SiO 2 , were altered in that time period, while those for dissolved organic carbon were not. Previous disturbances, including historical mining, also affect stream chemistry, and water‐quality impairment can be exacerbated by extreme climate events.

Colorado↗
USGS2018

Flooding alters plant-mediated carbon cycling independently of elevated atmospheric CO2 concentrations

Plant‐mediated processes determine carbon (C) cycling and storage in many ecosystems; how plant‐associated processes may be altered by climate‐induced changes in environmental drivers is therefore an essential question for understanding global C cycling. In this study, we hypothesize that environmental alterations associated with near‐term climate change can exert strong control on plant‐associated ecosystem C cycling and that investigations along an extended hydrologic gradient may give mechanistic insight into C cycling. We utilize a mesocosm approach to investigate the response of plant, soil, and gaseous C cycling to changing hydrologic regimes and elevated atmospheric carbon dioxide (CO 2) concentrations expected by 2100 in a coastal salt marsh in Louisiana, USA. Although elevated CO 2 had no significant effects on C cycling, we demonstrate that greater average flooding depth stimulated C exchange, with higher rates of labile C decomposition, plant CO 2 assimilation, and soil C respiration. Greater average flooding depth also significantly decreased the soil C pool and marginally increased the aboveground biomass C pool, leading to net losses in total C stocks. Further, flooding depths along an extended hydrologic gradient garnered insight into decomposition mechanisms that was not apparent from other data. In C‐4 dominated salt marshes, sea‐level rise will likely overwhelm effects of elevated CO 2 with climate change. Deeper flooding associated with sea‐level rise may decrease long‐term soil C pools and quicken C exchange between soil and atmosphere, thereby threatening net C storage in salt marsh habitats. Manipulative studies will be indispensable for understanding biogeochemical cycling under future conditions.

Journal of Geophysical Research: Biogeosciences↗
USGS2018

Impact of pore fluid chemistry on fine-grained sediment fabric and compressibility

Fines, defined here as grains or particles, less than 75 μm in diameter, exist nearly ubiquitously in natural sediment, even those classified as coarse. Macroscopic sediment properties, such as compressibility, which relates applied effective stress to the resulting sediment deformation, depend on the fabric of fines. Unlike coarse grains, fines have sizes and masses small enough to be more strongly influenced by electrical interparticle forces than by gravity. These electrical forces acting through pore fluids are influenced by pore fluid chemistry changes. Macroscopic property dependence on pore fluid chemistry must be accounted for in sediment studies involving subsurface flow and sediment stability analyses, as well as in engineered flow situations such as groundwater pollutant remediation, hydrocarbon migration, or other energy resource extraction applications. This study demonstrates how the liquid limit‐based electrical sensitivity index can be used to predict sediment compressibility changes due to pore fluid chemistry changes. Laboratory tests of electrical sensitivity, sedimentation, and compressibility illustrate mechanisms linking microscale and macroscale processes for selected pure, end‐member fines. A specific application considered here is methane extraction via depressurization of gas hydrate‐bearing sediment, which causes a dramatic pore water salinity drop concurrent with sediment being compressed by the imposed effective stress increase.

Journal of Geophysical Research↗
USGS2018

Interseismic ground deformation and fault slip rates in the greater San Francisco Bay Area from two decades of space geodetic data

The detailed spatial variations of strain accumulation and creep on major faults in the northern San Francisco Bay Area (North Bay), which are important for seismic potential and evaluation of natural hazards, remain poorly understood. Here we combine interferometric synthetic aperture radar data from the ERS‐1/2 and Envisat satellites between 1992 and 2010 with continuous and campaign GPS data to obtain a high spatial and temporal coverage of ground deformation of the North Bay. The SAR data from both ascending and descending orbits are combined to separate horizontal and vertical components of the deformation. We jointly invert the horizontal component of the mean velocities derived from these data to infer the deep strike‐slip rates on major locked faults. We use the estimated deep rates to simulate the long‐wavelength deformation due to interseismic elastic strain accumulation along these locked faults. After removing the long‐wavelength signal from the InSAR horizontal mean velocity field, we estimate fault‐parallel surface creep rates of up to 2 mm/year along the central section of the Rodgers Creek fault and surface creep rates ranging between 2 and 4 mm/year along the Concord fault. No surface creep is geodetically resolved along the West Napa and Green Valley fault zones. We identified characteristically repeating earthquakes on the Rodgers Creek fault, the West Napa fault, the Green Valley fault, and the Concord fault. Nontectonic deformation in the Geysers geothermal field and in Late Cenozoic basins (Rohnert Park and Sonoma basins) are also observed, likely due to hydrological and sediment‐compaction processes, respectively.

Journal of Geophysical Research B: Solid Earth↗
USGS2019

Thermal, deformation, and degassing remote sensing time-series (A.D. 2000-2017) at the 47 most active volcanoes in Latin America: Implications for volcanic systems

Volcanoes are hazardous to local and global populations, but only a fraction are continuously monitored by ground-based sensors. For example, in Latin America, more than 60% of Holocene volcanoes are unmonitored, meaning long-term multi-parameter datasets of volcanic activity are rare and sparse. We use satellite observations of degassing, thermal anomalies, and surface deformation spanning 17 years at 47 of the most active volcanoes in Latin America, and compare these datasets to ground-based observations archived by the Global Volcanism Program (GVP). This first comparison of multi-satellite time-series on a regional scale provides information regarding volcanic behavior during, non-, pre-, syn- and post-eruptive periods. For example, at Copahue volcano, deviations from background activity in all three types of satellite measurements were manifested months to years in advance of renewed eruptive activity in 2012. By quantifying the amount of degassing, thermal output, and deformation measured at each of these volcanoes, we test the classification of these volcanoes as open or closed volcanic systems. We find that ~28% of the volcanoes do not fall into either classification and the rest show elements of both, demonstrating a dynamic range of behavior that can change over time. Finally, we recommend how volcano monitoring could be improved through better coordination of available satellite-based capabilities and new instruments.

Journal of Geophysical Research↗
USGS2019

Heat and mass transport in a vapor-dominated hydrothermal area in Yellowstone National Park, USA: Inferences from magnetic, electrical, electromagnetic, subsurface temperature and diffuse CO2 flux measurements

Vapor‐dominated hydrothermal systems are characterized by localized and elevated heat and gas flux. In these systems, steam and gas ascend from a boiling water reservoir, steam condenses beneath a low‐permeability cap layer, and liquid water descends, driven by gravity (“heat pipe” model). We combine magnetic, electromagnetic, and geoelectrical methods and CO 2 flux and subsurface temperature measurements in the Solfatara Plateau Thermal Area in the Yellowstone Caldera to address several fundamental questions: (1) What are the structural and/or lithological controls on heat and mass transport in vapor‐dominated areas? (2) What is the geometry and size of convecting multiphase thermal plumes? (3) Are thermal plumes associated with subsurface rock alteration and demagnetization? Magnetic and electromagnetic data inversions suggest an asymmetric 50‐ to 100‐m thick basin of glacial deposits with the thickest part adjacent to the margin of a rhyolite flow. The 3‐D electrical conductivity model in the glacial basin reveals a narrow vertical conductor interpreted as a focused multiphase plume, which coincides at the ground surface with the heat and CO 2 flux maxima. The magnetic data suggest that destruction of magnetic minerals due to rock alteration associated with the hydrothermal plume occurs mainly near the ground surface. We propose a model where the buoyant multiphase plume forms in response to decompression, boiling, and phase separation of pressurized thermal groundwater that discharges from the brecciated base of a rhyolite flow into the basin of glacial deposits. Results from multiphase groundwater flow and heat transport numerical simulations corroborate the first‐order characteristics of this model.

Wyoming↗