SEARCH · Geology Reports

Results for “Interfaces”

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.

At least 883 records · Page 49Linked to original sources

USGS2018

Regeneration of Salicaceae riparian forests in the Northern Hemisphere: A new framework and management tool

Human activities on floodplains have severely disrupted the regeneration of foundation riparian shrub and tree species of the Salicaceae family ( Populus and Salix spp.) throughout the Northern Hemisphere. Restoration ecologists initially tackled this problem from a terrestrial perspective that emphasized planting. More recently, floodplain restoration activities have embraced an aquatic perspective, inspired by the expanding practice of managing river flows to improve river health (environmental flows) . However, riparian Salicaceae species occupy floodplain and riparian areas, which lie at the interface of both terrestrial and aquatic ecosystems along watercourses . Thus, their regeneration depends on a complex interaction of hydrologic and geomorphic processes that have shaped key life-cycle requirements for seedling establishment . Ultimately, restoration needs to integrate these concepts to succeed. However, while regeneration of Salicaceae is now reasonably well-understood, the literature reporting restoration actions on Salicaceae regeneration is sparse, and a specific theoretical framework is still missing. Here, we have reviewed 105 peer-reviewed published experiences in restoration of Salicaceae forests, including 91 projects in 10 world regions, to construct a decision tree to inform restoration planning through explicit links between the well-studied biophysical requirements of Salicaceae regeneration and 17 specific restoration actions, the most popular being planting (in 55% of the projects), land contouring (30%), removal of competing vegetation (30%), site selection (26%), and irrigation (24%). We also identified research gaps related to Salicaceae forest restoration and discuss alternative, innovative and feasible approaches that incorporate the human component.

Journal of Environmental Management↗

USGS2011

Polymethylene-interrupted fatty acids: Biomarkers for native and exotic mussels in the Laurentian Great Lakes

Freshwater organisms synthesize a wide variety of fatty acids (FAs); however, the ability to synthesize and/or subsequently modify a particular FA is not universal, making it possible to use certain FAs as biomarkers. Herein we document the occurrence of unusual FAs (polymethylene-interrupted fatty acids; PMI-FAs) in select freshwater organisms in the Laurentian Great Lakes. We did not detect PMI-FAs in: (a) natural seston from Lake Erie and Hamilton Harbor (Lake Ontario), (b) various species of laboratory-cultured algae including a green alga ( Scenedesmus obliquus ), two cyanobacteria ( Aphanizomenon flos-aquae and Synechococystis sp.), two diatoms ( Asterionella formosa , Diatoma elongatum ) and a chrysophyte ( Dinobryon cylindricum ) or, (c) zooplankton ( Daphnia spp., calanoid or cyclopoid copepods) from Lake Ontario, suggesting that PMI-FAs are not substantively incorporated into consumers at the phytoplankton–zooplankton interface. However, these unusual FAs comprised 4-6% of total fatty acids (on a dry tissue weight basis) of native fat mucket ( Lampsilis siliquoidea ) and plain pocketbook ( L. cardium ) mussels and in invasive zebra ( Dreissena polymorpha ) and quagga ( D. bugensis ) mussels. We were able to clearly partition Great Lakes' mussels into three separate groups (zebra, quagga, and native mussels) based solely on their PMI-FA profiles. We also provide evidence for the trophic transfer of PMI-FAs from mussels to various fishes in Lakes Ontario and Michigan, further underlining the potential usefulness of PMI-FAs for tracking the dietary contribution of mollusks in food web and contaminant-fate studies.

Michigan, Ontario, Prince Edward Island, Wisconsin↗

USGS2013

Great Lakes rivermouth ecosystems: scientific synthesis and management implications

At the interface of the Great Lakes and their tributary rivers lies the rivermouths, a class of aquatic ecosystem where lake and lotic processes mix and distinct features emerge. Many rivermouths are the focal point of both human interaction with the Great Lakes and human impacts to the lakes; many cities, ports, and beaches are located in rivermouth ecosystems, and these human pressures often degrade key ecological functions that rivermouths provide. Despite their ecological uniqueness and apparent economic importance, there has been relatively little research on these ecosystems as a class relative to studies on upstream rivers or the open-lake waters. Here we present a synthesis of current knowledge about ecosystem structure and function in Great Lakes rivermouths based on studies in both Laurentian rivermouths, coastal wetlands, and marine estuarine systems. A conceptual model is presented that establishes a common semantic framework for discussing the characteristic spatial features of rivermouths. This model then is used to conceptually link ecosystem structure and function to ecological services provided by rivermouths. This synthesis helps identify the critical gaps in understanding rivermouth ecology. Specifically, additional information is needed on how rivermouths collectively influence the Great Lakes ecosystem, how human alterations influence rivermouth functions, and how ecosystem services provided by rivermouths can be managed to benefit the surrounding socioeconomic networks.

Great Lakes↗

USGS2015

Application of the FluEgg model to predict transport of Asian carp eggs in the Saint Joseph River (Great Lakes tributary)

The Fluvial Egg Drift Simulator (FluEgg) is a three-dimensional Lagrangian model that simulates the movement and development of Asian carp eggs until hatching based on the physical characteristics of the flow field and the physical and biological characteristics of the eggs. This tool provides information concerning egg development and spawning habitat suitability including: egg plume location, egg vertical and travel time distribution, and egg-hatching risk. A case study of the simulation of Asian carp eggs in the Lower Saint Joseph River, a tributary of Lake Michigan, is presented. The river hydrodynamic input for FluEgg was generated in two ways — using hydroacoustic data and using HEC-RAS model data. The HEC-RAS model hydrodynamic input data were used to simulate 52 scenarios covering a broad range of flows and water temperatures with the eggs at risk of hatching ranging from 0 to 93% depending on river conditions. FluEgg simulations depict the highest percentage of eggs at risk of hatching occurs at the lowest discharge and at peak water temperatures. Analysis of these scenarios illustrates how the interactive relation among river length, hydrodynamics, and water temperature influence egg transport and hatching risk. An improved version of FluEgg, which more realistically simulates dispersion and egg development, is presented. Also presented is a graphical user interface that facilitates the use of FluEgg and provides a set of post-processing analysis tools to support management decision-making regarding the prevention and control of Asian carp reproduction in rivers with or without Asian carp populations.

Michigan↗

USGS2026

Geospatial data identify potentially restorable coastal wetland habitats in the Laurentian Great Lakes

Identifying areas with high restoration potential can guide investment and support management of coastal wetland habitats in the Laurentian Great Lakes. Restoration is a priority in the region since extensive degradation and alteration from anthropogenic activities reduced coastal wetlands to only a small fraction of historical extent. Our study evaluated the likelihood of coastal wetland restoration success using multi-criteria geospatial modeling in six restoration assessments focused on the U.S. portion of the Great Lakes coastline. We incorporated parameters into our model that positively influenced restoration potential, such as hydroperiod, wetland soil types, surface water flow networks, and conservation and recreation status of the land. Alternatively, factors that negatively influenced restoration potential were the proportion of impervious surfaces and other important human land uses of permeable soils. The resulting composite restoration index and online interfaces provide valuable decision-support tools for ecosystem managers. The modeled amount of restorable coastal wetland habitat and the characteristics of areas identified as restorable varied among the six restoration assessments, with the highest restoration index values occurring predominantly in low-elevation areas with frequent inundation that were near flowlines hydrologically connected to the Great Lakes. Conversely, areas with a greater proportion of impervious surfaces and anthropogenic land uses (e.g., golf courses) received low restoration index values. Our analysis highlights the heterogeneity of coastal wetland restoration potential throughout the Great Lakes, as well as the need for careful consideration of natural and anthropogenic land use factors to execute successful coastal wetland restoration.

Great Lakes↗

USGS2026

The Great Lakes aquatic tissue analysis repository (GLATAR)

Empiricists and modellers use information on energy density, proximate composition, stable isotopes, fatty acids, thiamine, and bioaccumulative tracers (e.g., PCBs and mercury) to understand the state and inter-relationships of aquatic food webs. Data exist in many published and unpublished sources, but are not consolidated in an easily accessible database that would serve as a vital resource to i) provide basic estimates of these diet-derived measures of body composition, ii) understand sources of variation in the underlying data, iii) facilitate exploration and development of data proxies, and iv) assist in study design. We designed GLATAR (Great Lakes Aquatic Tissue Analysis Repository, glatar.org ) to address this need. GLATAR is an open-access, searchable database linked to a web-based toolbox to visualise and generate user-defined summaries on diet-derived ecological metrics, with a focus on taxa of importance to the Great Lakes. GLATAR currently contains over 50,000 records on energy density, chemical tracers, and proximate body composition from 67 species of fish and 72 invertebrate taxa. We hope this user-friendly interface will entice others to upload their published and unpublished data to the repository, enriching the breadth of data accessible to researchers. In this way, GLATAR will become a ‘living’ and interactive resource for empiricists and modellers working in freshwater ecosystems as they make critical decisions related to growth, production, and consumption across a diverse group of economically and ecologically important aquatic species.

Journal of Great Lakes Research↗

USGS2014

Alternative waste residue materials for passive in situ prevention of sulfide-mine tailings oxidation: A field evaluation

Novel solutions for sulfide-mine tailings remediation were evaluated in field-scale experiments on a former tailings repository in northern Sweden. Uncovered sulfide-tailings were compared to sewage-sludge biosolid amended tailings over 2 years. An application of a 0.2 m single-layer sewage-sludge amendment was unsuccessful at preventing oxygen ingress to underlying tailings. It merely slowed the sulfide-oxidation rate by 20%. In addition, sludge-derived metals (Cu, Ni, Fe, and Zn) migrated and precipitated at the tailings-to-sludge interface. By using an additional 0.6 m thick fly-ash sealing layer underlying the sewage sludge layer, a solution to mitigate oxygen transport to the underlying tailings and minimize sulfide-oxidation was found. The fly-ash acted as a hardened physical barrier that prevented oxygen diffusion and provided a trap for sludge-borne metals. Nevertheless, the biosolid application hampered the application, despite the advances in the effectiveness of the fly-ash layer, as sludge-borne nitrate leached through the cover system into the underlying tailings, oxidizing pyrite. This created a 0.3 m deep oxidized zone in 6-years. This study highlights that using sewage sludge in unconventional cover systems is not always a practical solution for the remediation of sulfide-bearing mine tailings to mitigate against sulfide weathering and acid rock drainage formation.

Skellefteï¿½;Stockholm↗

USGS2008

A seepage meter designed for use in flowing water

Seepage meters provide one of the most direct means to measure exchange of water across the sediment-water interface, but they generally have been unsuitable for use in fluvial settings. Although the seepage bag can be placed inside a rigid container to minimize velocity head concerns, the seepage cylinder installed in the sediment bed projects into and disrupts the flow field, altering both the local-scale fluid exchange as well as measurement of that exchange. A low-profile seepage meter designed for use in moving water was tested in a seepage meter flux tank where both current velocity and seepage velocity could be controlled. The conical seepage cylinder protrudes only slightly above the sediment bed and is connected via tubing to a seepage bag or flowmeter positioned inside a rigid shelter that is located nearby where current velocity is much slower. Laboratory and field tests indicate that the net effect of the small protrusion of the seepage cylinder into the surface water flow field is inconsequentially small for surface water currents up to 65 cm s-1. Current velocity affects the variability of seepage measurements; seepage standard deviation increased from ???2 to ???6 cm d-1 as current velocity increased from 9 to 65 cm s-1. Substantial bias can result if the shelter is not placed to minimize hydraulic gradient between the bag and the seepage cylinder.

Journal of Hydrology↗

USGS2011

Subsurface transport of orthophosphate in five agricultural watersheds, USA

Concentrations of dissolved orthophosphate (ortho P) in the unsaturated zone, groundwater, tile drains, and groundwater/stream water interfaces were assessed in five agricultural watersheds to determine the potential for subsurface transport. Concentrations of iron oxides were measured in the aquifer material and adsorption of ortho P on oxide surfaces was assessed by geochemical modeling. Attenuation of ortho P in these aquifers was attributed primarily to sorption onto iron oxides, and in one location onto clay minerals. Only one location showed a clear indication of phosphorus transport to a stream from groundwater discharge, although groundwater did contribute to the stream load elsewhere. Subsurface ortho P movement at a site in California resulted in a plume down gradient from orchards, which was attenuated by a 200 m thick riparian zone with natural vegetation. Iron oxides had an effect on phosphorus movement and concentrations at all locations, and groundwater chemistry, especially pH, exerted a major control on the amount of phosphorus adsorbed. Groundwater pH at a site in Maryland was below 5 and that resulted in complete sequestration of phosphorus and no movement toward the stream. Geochemical modeling indicated that as the surfaces approached saturation, groundwater concentrations of ortho P rise rapidly.

Journal of Hydrology↗

USGS2014

Quantity and quality of groundwater discharge in a hypersaline lake environment

Geophysical and geochemical surveys were conducted to understand groundwater discharge to Great Salt Lake (GSL) and assess the potential significance of groundwater discharge as a source of selenium (Se). Continuous resistivity profiling (CRP) focusing below the sediment/water interface and fiber-optic distributed temperature sensing (FO-DTS) surveys were conducted along the south shore of GSL. FO-DTS surveys identified persistent cold-water temperature anomalies at 10 separate locations. Seepage measurements were conducted at 17 sites (mean seepage rate = 0.8 cm/day). High resistivity anomalies identified by the CRP survey were likely a mirabilite (Na 2 SO 4 ·10H 2 O) salt layer acting as a semi-confining layer for the shallow groundwater below the south shore of the lake. Positive seepage rates measured along the near-shore areas of GSL indicate that a ∼1-m thick oolitic sand overlying the mirabilite layer is likely acting as a shallow, unconfined aquifer. Using the average seepage rate of 0.8 cm/day over an area of 1.6 km 2 , an annual Se mass loading to GSL of 23.5 kg was estimated. Determination of R/Ra values (calculated 3 He/ 4 He ratio over the present-day atmospheric 3 He/ 4 He ratio) <1 and tritium activities of 1.2–2.0 tritium units in groundwater within and below the mirabilite layer indicates a convergence of regional and local groundwater flow paths discharging into GSL. Groundwater within and below the mirabilite layer obtains its high sulfate salinity from the dissolution of mirabilite. The δ 34 S and δ 18 O isotopic values in samples of dissolved sulfate from the shallow groundwater below the mirabilite are almost identical to the isotopic signature of the mirabilite core material. The saturation index calculated for groundwater samples using PHREEQC indicates the water is at equilibrium with mirabilite. Water samples collected from GSL immediately off shore contained Se concentrations that were 3–4 times higher than other sampling sites >25 km offshore from the study site and may be originating from less saline groundwater seeps mixing with the more saline water from GSL. Additional evidence for mixing with near shore seeps is found in the δD and δ 18 O isotopic values and Br:Cl ratios. Geochemical modeling for a water sample collected in the vicinity of the study area indicates that under chemically reducing conditions, arsenic- (As) bearing minerals could dissolve while Se-bearing minerals will likely precipitate out of solution, possibly explaining why the shallow groundwater below and within the mirabilite salt layer contains low concentrations of Se (0.9–2.3 μg/L).

Utah↗

USGS2014

Practical limitations on the use of diurnal temperature signals to quantify groundwater upwelling

Groundwater upwelling to streams creates unique habitat by influencing stream water quality and temperature; upwelling zones also serve as vectors for contamination when groundwater is degraded. Temperature time series data acquired along vertical profiles in the streambed have been applied to simple analytical models to determine rates of vertical fluid flux. These models are based on the downward propagation characteristics (amplitude attenuation and phase-lag) of the surface diurnal signal. Despite the popularity of these models, there are few published characterizations of moderate-to-strong upwelling. We attribute this limitation to the thermodynamics of upwelling, under which the downward conductive signal transport from the streambed interface occurs opposite the upward advective fluid flux. Governing equations describing the advection–diffusion of heat within the streambed predict that under upwelling conditions, signal amplitude attenuation will increase, but, counterintuitively, phase-lag will decrease. Therefore the extinction (measurable) depth of the diurnal signal is very shallow, but phase lag is also short, yielding low signal to noise ratio and poor model sensitivity. Conversely, amplitude attenuation over similar sensor spacing is strong, yielding greater potential model sensitivity. Here we present streambed thermal time series over a range of moderate to strong upwelling sites in the Quashnet River, Cape Cod, Massachusetts. The predicted inverse relationship between phase-lag and rate of upwelling was observed in the field data over a range of conditions, but the observed phase-lags were consistently shorter than predicted. Analytical solutions for fluid flux based on signal amplitude attenuation return results consistent with numerical models and physical seepage meters, but the phase-lag analytical model results are generally unreasonable. Through numerical modeling we explore reasons why phase-lag may have been over-predicted by the analytical models, and develop guiding relations of diurnal temperature signal extinction depth based on stream diurnal signal amplitude, upwelling magnitude, and streambed thermal properties that will be useful in designing future experiments.

Journal of Hydrology↗

USGS2019

Multi-scale preferential flow processes in an urban streambed under variable hydraulic conditions

Spatially preferential flow processes occur at nested scales at the sediment-water interface (SWI), due in part to sediment heterogeneities, which may be enhanced in flashy urban streams with heavy road sand influence. However, several factors, including the flow-rate dependence of preferential hyporheic flow and discrete groundwater discharge zones are commonly overlooked in reach-scale models of groundwater/surface water exchange. Using a series of controlled-head tracer-injection experiments coupled with cm-scale geophysics within the highly reactive upper 30 cm of the hyporheic zone of an urban stream, we quantified the flow dependence of local less-mobile porosity volume, mass-transfer rate coefficient, and the resulting local residence time in the less-mobile pore space at three controlled downward fluid fluxes (0.8, 2, and 3 m/d). Experiments were performed in two adjacent streambed locations, representing different sediment bulk vertical permeability. Less-mobile porosity parameters were generally substantial and similar between the two streambed locations; though a more competent, thin, organic layer at ∼15 cm depth in one location strongly impacted tracer loading, flushing dynamics, and local residence times. Increased downward flux led to (1) a decrease in less-mobile porosity residence time in all experiments, and (2) an increase in less-mobile porosity fraction for most experiments. Additionally, at the larger stream reach-scale, surface electrodes for electrical resistivity measurement were installed along 22 m of the wetted stream channel. These surface electrode measurements were collected during a natural storm flow event, which revealed widespread, short-term, flushing (e.g. <3 h) of the hyporheic zone with stream water, followed by longer-term (e.g. >60 h) flushing of the SWI with riparian zone groundwater. Flow dependence of preferential hyporheic zone flowpaths, like in the controlled tracer experiments, was also observed in these reach-scale electrical resistivity tomography measurements. Our findings reveal that the spatial and temporal dependence of preferential flow processes create highly dynamic SWI conditions that will affect the physical and coupled biogeochemical functions of the SWI in urbanized, sand-impacted streams.

Journal of Hydrology↗

USGS2023

Modeling the dynamic penetration depth of post-1950s water in unconfined aquifers using environmental tracers: Central Valley, California

The penetration depth of post-1950s recharge (D-1950) in aquifers is a marker that is frequently used to identify groundwater that is susceptible to anthropogenic contamination. Here, we compute D-1950 values at wells, interpolate them in space, and project them across time to map the moving front of modern recharge in four dimensions in the Central Valley aquifer system, California, USA. Tracers of groundwater age (tritium, carbon-14, noble gases, sulfur hexafluoride, and chlorofluorocarbons) were collected at 650 wells spatially distributed throughout the Central Valley and were fit to a lumped-parameter model that assumes a logarithmic age-depth profile in the aquifer. For samples where tritium was present (>0.3 tritium units), the model was used to predict D-1950 at wells screened above or across the modern-premodern interface (n = 484). Wells with samples where tritium was absent (≤0.3 tritium units) were used to define the depth beyond which groundwater is completely premodern (n = 166). Predicted D-1950 values were below the depth of screen bottoms for wells where groundwater is completely modern, and above the depth of screen tops for wells where groundwater is completely premodern. The interpolated surface of D-1950 is dynamic, less prone to extreme values, and produces maps with lower interpolation errors due to a higher spatial density of wells than maps based on the depth of premodern groundwater. Between 2005 and 2025, D-1950 is expected to deepen by 11 and 12 m in the northern and southern parts of the Central Valley, respectively. Areas where D-1950 increases rapidly are likely to see increases in nitrate and other anthropogenic contaminants associated with the downward moving front of modern water.

California↗

USGS2023

Incorporating temperature into seepage loss estimates for a large unlined irrigation canal

Quantifying seepage losses from unlined irrigation canals is necessary to improve water use and conservation. The use of heat as a tracer is widely used in quantifying seepage rates across the sediment–water interface. In this study, field observations and two-dimensional numerical models were used to simulate seepage losses during the 2018 and 2019 irrigation season in the Truckee Canal system. Nineteen transects were instrumented with temperature probes and stage recording devices for inverse modeling to derive seepage flux and volumetric losses over the 39 km length of canal. The numerical models for each transect were calibrated and validated using the two-year dataset. Soil zones and observation data were used in each numerical model to help guide calibration of vertical and lateral heat and fluid fluxes. Model simulations were used to derive multivariable regression equations that consider stage, temperature, and hydraulic gradient. The results demonstrate the value of long-term datasets that illustrate the seasonality of groundwater levels, siltation, stage, and temperature on seepage rates. Seepage rates estimated by the numerical models range from 0.16 to 4.6 m 3 /d m −1 . Total annual volumetric losses estimated for 2018 and 2019 were 1.6 × 10 -2 to 1.2 × 10 -2 km 3 , respectively. The seepage losses estimated by this study account for 32 % to 41 % of the inflow volumes. Regression models were able to reproduce seepage time-series simulated by the numerical models reasonably well. In arid environments, water diverted into irrigation canals may be influenced by seasonal variations in temperature sufficient to influence the water accounting of conveyed surface flows.

Nevada↗

USGS2024

Indications of preferential groundwater seepage feeding northern peatland pools

Groundwater seepage from underlying permeable glacial sedimentary structures, such as eskers, has been hypothesized to directly feed pools in northern peat bogs. These hypotheses directly contradict classical peat bog models for ombrogenous systems, wherein meteoric water is the sole water input to these systems. Variations in the underlying mineral sediment in contact with the peat imply that unrecognized hydrogeologic connectivity may exist with pools in northern peat bogs, particularly where high permeability materials are in contact with the peat. Seepage dynamics originating from these structural variations were investigated using a suite of thermal and hydrogeophysical methods deployed around pools in a peat bog of northeastern Maine, USA. Thermal characterization methods mapped anomalies that were confirmed as matrix seepage or preferential flow pathways (PFPs). Geochemical methods were employed at identified thermal anomalies to confirm upwelling of solute-rich groundwater. Conduits around pools were associated with surficial terminations of suspected peat pipes, based on the inference of pathways extending down into the peat, that focus flow through PFPs in the peat matrix. Discharge also occurred through the peat matrix adjacent to suspected pipe structures and matrix seepage rates were quantified using analysis of diurnal temperature signals recorded at multiple depths. Seepage rates, with a maximum of nearly 0.4 m/d, were measured at localized points around pools. Periods of synchronized temperatures paired with highly muted diurnal temperature signals, recorded in diurnal temperature with depth data, were interpreted qualitatively as activation of strong upward discharge rates through suspected peat pipes. These time periods correlated strongly with local precipitation events around the peatland. Ground-penetrating radar surveys revealed discontinuities in the low permeability glacio-marine clay at the mineral sediment-peat interface, interpreted to be regional glacial esker deposits, which were located beneath and around pools. Heat tracing, specific conductance contrasts, seepage rates, and trace metal concentrations all imply groundwater seepage originating from underlying permeable glacial esker deposits and directly sourcing pools. Preferential groundwater inputs into northern peat bogs may play a key role in developing and maintaining pool systems, with enhanced solute transport impacting peatland ecology, water resources, and carbon cycling.

Maine↗

USGS2025

Impact of gas/liquid phase change of CO2 during injection for sequestration

CO 2 sequestration in deep saline formations is an effective and important process to control the rapid rise in CO 2 emissions. The process of injecting CO 2 requires reliable predictions of the stress in the formation and the fluid pressure distributions – particularly since monitoring of the CO 2 migration is difficult – to mitigate leakage, prevent induced seismicity, and analyze wellbore stability. A key aspect of CO 2 is the gas–liquid phase transition at the temperatures and pressures of relevance to leakage and sequestration, which has been recognized as being critical for accurate predictions but has been challenging to model without ad hoc empiricisms. This paper presents a robust multiphase thermodynamics-based poromechanics model to capture the complex phase transition behavior of CO 2 and predict the stress and pressure distribution under super- and sub- critical conditions during the injection process. A finite element implementation of the model is applied to analyze the behavior of a multiphase porous system with CO 2 as it displaces the fluid brine phase. We find that if CO 2 undergoes a phase transition in the geologic reservoir, the spatial variation of the density is significantly affected, and the migration mobility of CO 2 decreases in the reservoir. A key feature of our approach is that we do not a priori assume the location of the CO 2 gas/liquid interface – or even if it occurs at all – but rather, this is a prediction of the model, along with the spatial variation of the phase of CO 2 and the change of the saturation profile due to the phase change.

Journal of the Mechanics and Physics of Solids↗

USGS2004

Acoustic measurements of the 1999 basaltic eruption of Shishaldin volcano, Alaska 1. Origin of Strombolian activity

The 1999 basaltic eruption of Shishaldin volcano (Alaska, USA) displayed both classical Strombolian activity and an explosive Subplinian plume. Strombolian activity at Shishaldin occurred in two major phases following the Subplinian activity. In this paper, we use acoustic measurements to interpret the Strombolian activity. Acoustic measurements of the two Strombolian phases show a series of explosions that are modeled by the vibration of a large overpressurised cylindrical bubble at the top of the magma column. Results show that the bubble does not burst at its maximum radius, as expected if the liquid film is stretched beyond its elasticity. But bursting occurs after one cycle of vibration, as a consequence of an instability of the air-magma interface close to the bubble minimum radius. During each Strombolian period, estimates of bubble length and overpressure are calculated. Using an alternate method based on acoustic power, we estimate gas velocity to be 30-60 m/s, in very good agreement with synthetic waveforms. Although there is some variation within these parameters, bubble length and overpressure for the first Strombolian phase are found to be ??? 82 ?? 11 m and 0.083 MPa. For the second Strombolian phase, bubble length and overpressure are estimated at 24 ?? 12 m and 0.15 MPa for the first 17 h after which bubble overpressure shows a constant increase, reaching a peak of 1.4 MPa, just prior to the end of the second Strombolian phase. This peak suggests that, at the time, the magma in the conduit may contain a relatively large concentration of small bubbles. Maximum total gas volume and gas fluxes at the surface are estimated to be 3.3 ?? 107 and 2.9 ?? 103 m3/s for the first phase and 1.0 ?? 108 and 2.2 ?? 103 m3/s for the second phase. This gives a mass flux of 1.2 ?? 103 and 8.7 ?? 102 kg/s, respectively, for the first and the second Strombolian phases. ?? 2004 Elsevier B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗

USGS2006

Geodetic observations and modeling of magmatic inflation at the Three Sisters volcanic center, central Oregon Cascade Range, USA

Tumescence at the Three Sisters volcanic center began sometime between summer 1996 and summer 1998 and was discovered in April 2001 using interferometric synthetic aperture radar (InSAR). Swelling is centered about 5 km west of the summit of South Sister, a composite basaltic-andesite to rhyolite volcano that last erupted between 2200 and 2000 yr ago, and it affects an area &sim;20 km in diameter within the Three Sisters Wilderness. Yearly InSAR observations show that the average maximum displacement rate was 3–5 cm/yr through summer 2001, and the velocity of a continuous GPS station within the deforming area was essentially constant from June 2001 to June 2004. The background level of seismic activity has been low, suggesting that temperatures in the source region are high enough or the strain rate has been low enough to favor plastic deformation over brittle failure. A swarm of about 300 small earthquakes ( M max = 1.9) in the northeast quadrant of the deforming area on March 23–26, 2004, was the first notable seismicity in the area for at least two decades. The U.S. Geological Survey (USGS) established tilt-leveling and EDM networks at South Sister in 1985–1986, resurveyed them in 2001, the latter with GPS, and extended them to cover more of the deforming area. The 2001 tilt-leveling results are consistent with the inference drawn from InSAR that the current deformation episode did not start before 1996, i.e., the amount of deformation during 1995–2001 from InSAR fully accounts for the net tilt at South Sister during 1985–2001 from tilt-leveling. Subsequent InSAR, GPS, and leveling observations constrain the source location, geometry, and inflation rate as a function of time. A best-fit source model derived from simultaneous inversion of all three datasets is a dipping sill located 6.5 ± 2.5 km below the surface with a volume increase of 5.0 × 10 6 ± 1.5 × 10 6 m 3 /yr (95% confidence limits). The most likely cause of tumescence is a pulse of basaltic magma intruding the upper crust along the brittle–ductile interface — a process that must occur episodically beneath the Cascade Range but in the past would have escaped detection in the absence of unusual seismicity. We speculate that such intrusive episodes last from days to years and are separated by quiescent periods of decades to centuries. The likelihood that the current episode at Three Sisters will culminate in an eruption is judged to be low, but the impact of an eruption could be great. The USGS has updated its volcano hazards assessment for the Three Sisters region, notified appropriate agencies and the public, and is helping to prepare an emergency coordination and communication plan.

Central Oregon Cascade Range, Three Sisters volcan↗