SEARCH · USGS Science

Results for “Journal of Geophysical Research- Earth Surface”

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 433 records · Page 24Linked to original sources

USGS2018

Findings and lessons learned from the assessment of the Mexico-United States transboundary San Pedro and Santa Cruz aquifers: The utility of social science in applied hydrologic research

Study Region This study region encompasses the Transboundary San Pedro and Santa Cruz aquifers which are shared between the states of Sonora (Mexico) and Arizona (US). Special regional considerations include a semi-arid climate, basin-fill aquifers with predominantly montane recharge areas, economic drivers in the mining, trade, and military sectors, groundwater-dependent cities with expanding cones of depression, interbasin groundwater transfers, ground- and surface-water contamination, and protected aquatic and riparian habitats that act as significant migration corridors for hundreds of species, including some that are threatened and endangered. Study Focus We focus on lessons learned from the hydrologic assessment of the Transboundary San Pedro and Santa Cruz aquifers. We conducted the work, in two phases: (1) laying the groundwork and (2) implementation. The “laying the groundwork” phase consisted of binational meetings with stakeholders and key actors (agencies and individuals), and the development of an understanding of the physical, institutional, historical, and socio-political context. This led to signing of the binational Transboundary Aquifer Assessment Program (TAAP) agreement in 2009 and detailed the process for cooperation and coordination in the assessment of shared aquifers. The implementation phase began with an agreement to proceed with the study of four “focus” aquifers (Santa Cruz, San Pedro, Mesilla (Conejos-Médanos in Mexico), and Hueco Bolson (Bolsón del Hueco in Mexico)) and development of associated technical teams. Though we do include a brief discussion of the lessons learned from the physical science portion of the study, the results have been described and published elsewhere. The bulk of the paper instead focuses on the findings and lessons learned from the integration of social-science perspectives into a largely physical-science based program, since there is a growing recognition of the need for this type of approach especially in the management and assessment of transboundary aquifers. New Hydrological Insights for the Region The Sonora-Arizona effort succeeded because both countries were adequately represented, and because of flexibility of skills and ability of teams comprising both university and government scientists. Teams included social and earth scientists. Including the social sciences was critical to research design and implementation, and to addressing the cultural, institutional, and socio-political contexts of transboundary aquifer assessment. Significant components of the continuing implementation phase include strategic planning, data compilation and analysis, cross-border integration of datasets, geophysical and geochemical surveys, and internal, peer, and stakeholder engagement.

Arizona, Sonora↗

USGS2022

Quantifying large-scale surface change using SAR amplitude images: Crater morphology changes during the 2019-2020 Shishaldin Volcano eruption

Morphological processes often induce meter-scale elevation changes. When a volcano erupts, tracking such processes provides insights into the style and evolution of eruptive activity and related hazards. Compared to optical remote-sensing products, synthetic aperture radar (SAR) observes surface change during inclement weather and at night. Differential SAR interferometry estimates phase change between SAR acquisitions and is commonly applied to quantify deformation. However, large deformation or other coherence loss can limit its use. We develop a new approach applicable when repeated digital elevation models (DEMs) cannot be otherwise retrieved. Assuming an isotropic radar cross-section, we estimate meter-scale vertical morphological change directly from SAR amplitude images via an optimization method that utilizes a high-quality DEM. We verify our implementation through simulation of a collapse feature that we modulate onto topography. We simulate radar effects and recover the simulated collapse. To validate our method, we estimate elevation changes from TerraSAR-X stripmap images for the 2011–2012 eruption of Mount Cleveland. Our results reproduce those from two previous studies; one that used the same dataset, and another based on thermal satellite data. By applying this method to the 2019–2020 eruption of Shishaldin Volcano, Alaska, we generate elevation change time series from dozens of co-registered TerraSAR-X high-resolution spotlight images. Our results quantify previously unresolved cone growth in November 2019, collapses associated with explosions in December–January, and further changes in crater elevations into spring 2020. This method can be used to track meter-scale morphology changes for ongoing eruptions with low latency as SAR imagery becomes available.

Alaska↗

USGS2012

Fluid expulsion sites on the Cascadia accretionary prism: mapping diagenetic deposits with processed GLORIA imagery

Point-discharge fluid expulsion on accretionary prisms is commonly indicated by diagenetic deposition of calcium carbonate cements and gas hydrates in near-surface (<10 m below seafloor; mbsf) hemipelagic sediment. The contrasting clastic and diagenetic lithologies should be apparent in side scan images. However, sonar also responds to variations in bottom slope, so unprocessed images mix topographic and lithologic information. We have processed GLORIA imagery from the Oregon continental margin to remove topographic effects. A synthetic side scan image was created initially from Sea Beam bathymetric data and then was subtracted iteratively from the original GLORIA data until topographic features disappeared. The residual image contains high-amplitude backscattering that we attribute to diagenetic deposits associated with fluid discharge, based on submersible mapping, Ocean Drilling Program drilling, and collected samples. Diagenetic deposits are concentrated (1) near an out-of-sequence thrust fault on the second ridge landward of the base of the continental slope, (2) along zones characterized by deep-seated strikeslip faults that cut transversely across the margin, and (3) in undeformed Cascadia Basin deposits which overlie incipient thrust faults seaward of the toe of the prism. There is no evidence of diagenetic deposition associated with the frontal thrust that rises from the dècollement. If the dècollement is an important aquifer, apparently the fluids are passed either to the strike-slip faults which intersect the dècollement or to the incipient faults in Cascadia Basin for expulsion. Diagenetic deposits seaward of the prism toe probably consist dominantly of gas hydrates.

Cascadia Basin↗

USGS2014

Ice sheet load cycling and fluid underpressures in the Eastern Michigan Basin, Ontario, Canada

Strong fluid underpressures have been detected in Paleozoic strata in the eastern Michigan Basin, with hydraulic heads reaching ~400 m below land surface (~4 MPa underpressure) and ~200 m below sea level in strata where unusually low permeabilities (~10 −20 –10 −23  m 2 ) were measured in situ. Multiple glaciations, including three with as much as 3 km of ice cover at the site in the last 120 ka, suggest a causal link with the underpressures. We examined this possibility using a one-dimensional groundwater flow model incorporating mechanical loading from both ice weight and lithospheric flexure. Because hydrologic and mechanical changes during glaciation are not well characterized and subsurface properties are imperfectly known, the model was used inversely to estimate flexural loads and loosely constrained permeabilities by matching observed pressures. Acceptable matches were obtained for a surprisingly wide range of scenarios with permeabilities close to measured values and plausible flexural loads. Matches were not obtained when too many parameters were preselected, or when permeabilities were constrained to be significantly larger than measured values. In successful model runs groundwater expulsion under glacial-mechanical loads caused the underpressuring, and flexural loads were important if aquifer and sub-glacial pressures were significantly elevated during glaciation. Simulated fluid pressures in the low-permeability strata fluctuated by 30–40 MPa during glacial cycles but resulted in advective transport of only tens of meters or less. Although other mechanisms cannot be ruled out, we conclude that glacial-mechanical forcing of a water-saturated system can explain the observed underpressures.

Ontario↗

USGS2012

Heat flow in vapor dominated areas of the Yellowstone Plateau volcanic field: implications for the thermal budget of the Yellowstone Caldera

Characterizing the vigor of magmatic activity in Yellowstone requires knowledge of the mechanisms and rates of heat transport between magma and the ground surface. We present results from a heat flow study in two vapor dominated, acid-sulfate thermal areas in the Yellowstone Caldera, the 0.11 km 2 Obsidian Pool Thermal Area (OPTA) and the 0.25 km 2 Solfatara Plateau Thermal Area (SPTA). Conductive heat flux through a low permeability layer capping large vapor reservoirs is calculated from soil temperature measurements at >600 locations and from laboratory measurements of soil properties. The conductive heat output is 3.6 ± 0.4 MW and 7.5 ± 0.4 MW from the OPTA and the SPTA, respectively. The advective heat output from soils is 1.3 ± 0.3 MW and 1.2 ± 0.3 MW from the OPTA and the SPTA, respectively and the heat output from thermal pools in the OPTA is 6.8 ± 1.4 MW. These estimates result in a total heat output of 11.8 ± 1.4 MW and 8.8 ± 0.4 MW from OPTA and SPTA, respectively. Focused zones of high heat flux in both thermal areas are roughly aligned with regional faults suggesting that faults in both areas serve as conduits for the rising acid vapor. Extrapolation of the average heat flux from the OPTA (103 ± 2 W·m −2 ) and SPTA (35 ± 3 W·m −2 ) to the ~35 km 2 of vapor dominated areas in Yellowstone yields 3.6 and 1.2 GW, respectively, which is less than the total heat output transported by steam from the Yellowstone Caldera as estimated by the chloride inventory method (4.0 to 8.0 GW).

Journal of Geophysical Research B: Solid Earth↗

USGS2001

Resolution analysis of finite fault source inversion using one- and three-dimensional Green's functions 1. Strong motions

We develop a methodology to perform finite fault source inversions from strong motion data using Green's functions (GFs) calculated for a three-dimensional (3-D) velocity structure. The 3-D GFs are calculated numerically by inserting body forces at each of the strong motion sites and then recording the resulting strains along the target fault surface. Using reciprocity, these GFs can be recombined to represent the ground motion at each site for any (heterogeneous) slip distribution on the fault. The reciprocal formulation significantly reduces the required number of 3-D finite difference computations to at most 3NS, where NS is the number of strong motion sites used in the inversion. Using controlled numerical resolution tests, we have examined the relative importance of accurate GFs for finite fault source inversions which rely on near-source ground motions. These experiments use both 1-D and 3-D GFs in inversions for hypothetical rupture models in order (1) to analyze the ability of the 3-D methodology to resolve trade-offs between complex source phenomena and 3-D path effects, (2) to address the sensitivity of the inversion results to uncertainties in the 3-D velocity structure, and (3) to test the adequacy of the 1-D GF method when propagation effects are known to be three-dimensional. We find that given "data" from a prescribed 3-D Earth structure, the use of well-calibrated 3-D GFs in the inversion provides very good resolution of the assumed slip distribution, thus adequately separating source and 3-D propagation effects. In contrast, using a set of inexact 3-D GFs or a set of hybrid 1-D GFs allows only partial recovery of the slip distribution. These findings suggest that in regions of complex geology the use of well-calibrated 3-D GFs has the potential for increased resolution of the rupture process relative to 1-D GFs. However, realizing this full potential requires that the 3-D velocity model and associated GFs should be carefully validated against the true 3-D Earth structure before performing the inverse problem with actual data.

Journal of Geophysical Research B: Solid Earth↗

USGS1998

Paleoseismic investigations in the Santa Cruz mountains, California: Implications for recurrence of large-magnitude earthquakes on the San Andreas fault

Trenching, microgeomorphic mapping, and tree ring analysis provide information on timing of paleoearthquakes and behavior of the San Andreas fault in the Santa Cruz mountains. At the Grizzly Flat site alluvial units dated at 1640–1659 A.D., 1679–1894 A.D., 1668–1893 A.D., and the present ground surface are displaced by a single event. This was the 1906 surface rupture. Combined trench dates and tree ring analysis suggest that the penultimate event occurred in the mid-1600 s, possibly in an interval as narrow as 1632–1659 A.D. There is no direct evidence in the trenches for the 1838 or 1865 earthquakes, which have been proposed as occurring on this part of the fault zone. In a minimum time of about 340 years only one large surface faulting event (1906) occurred at Grizzly Flat, in contrast to previous recurrence estimates of 95–110 years for the Santa Cruz mountains segment. Comparison with dates of the penultimate San Andreas earthquake at sites north of San Francisco suggests that the San Andreas fault between Point Arena and the Santa Cruz mountains may have failed either as a sequence of closely timed earthquakes on adjacent segments or as a single long rupture similar in length to the 1906 rupture around the mid-1600 s. The 1906 coseismic geodetic slip and the late Holocene geologic slip rate on the San Francisco peninsula and southward are about 50–70% and 70% of their values north of San Francisco, respectively. The slip gradient along the 1906 rupture section of the San Andreas reflects partitioning of plate boundary slip onto the San Gregorio, Sargent, and other faults south of the Golden Gate. If a mid-1600 s event ruptured the same section of the fault that failed in 1906, it supports the concept that long strike-slip faults can contain master rupture segments that repeat in both length and slip distribution. Recognition of a persistent slip rate gradient along the northern San Andreas fault and the concept of a master segment remove the requirement that lower slip sections of large events such as 1906 must fill in on a periodic basis with smaller and more frequent earthquakes.

Journal of Geophysical Research B: Solid Earth↗

USGS1999

Formation of natural gas hydrates in marine sediments 1. Conceptual model of gas hydrate growth conditioned by host sediment properties

The stability of submarine gas hydrates is largely dictated by pressure and temperature, gas composition, and pore water salinity. However, the physical properties and surface chemistry of deep marine sediments may also affect the thermodynamic state, growth kinetics, spatial distributions, and growth forms of clathrates. Our conceptual model presumes that gas hydrate behaves in a way analogous to ice in a freezing soil. Hydrate growth is inhibited within fine-grained sediments by a combination of reduced pore water activity in the vicinity of hydrophilic mineral surfaces, and the excess internal energy of small crystals confined in pores. The excess energy can be thought of as a "capillary pressure" in the hydrate crystal, related to the pore size distribution and the state of stress in the sediment framework. The base of gas hydrate stability in a sequence of fine sediments is predicted by our model to occur at a lower temperature (nearer to the seabed) than would be calculated from bulk thermodynamic equilibrium. Capillary effects or a build up of salt in the system can expand the phase boundary between hydrate and free gas into a divariant field extending over a finite depth range dictated by total methane content and pore-size distribution. Hysteresis between the temperatures of crystallization and dissociation of the clathrate is also predicted. Growth forms commonly observed in hydrate samples recovered from marine sediments (nodules, and lenses in muds; cements in sands) can largely be explained by capillary effects, but kinetics of nucleation and growth are also important. The formation of concentrated gas hydrates in a partially closed system with respect to material transport, or where gas can flush through the system, may lead to water depletion in the host sediment. This "freeze-drying" may be detectable through physical changes to the sediment (low water content and overconsolidation) and/or chemical anomalies in the pore waters and metastable presence of free gas within the normal zone of hydrate stability.

Journal of Geophysical Research B: Solid Earth↗

USGS1996

Type of faulting and orientation of stress and strain as a function of space and time in Kilauea's south flank, Hawaii

Earthquake focal mechanisms of events occurring between 1972 and 1992 in the south flank of Kilauea volcano, Hawaii, are used to infer the state of stress and strain as a function of time and space. We have determined 870 fault plane solutions from P wave first motion polarities for events with magnitudes M L ≥ 2.5 and depth ranging between 6 and 12 km. Faulting is characterized by a mixture of decollement, reverse, and normal faults. Most large earthquakes with magnitude M < 7 slip on reverse faults striking NE at 40° and dipping SE between 60° and 70°. In Hawaii, the earthquakes with M > 7 rupture the decollement plane, since it is the only surface large enough to generate magnitude 7 or larger earthquakes. The percentage of reverse faulting events is high compared to the decollement and normal faulting mechanisms for the period 1972–1983. The percentage of decollement type focal mechanisms becomes dominant after 1983. This pattern of faulting activity suggests that pressure was building up within Kilauea's rift zone prior to the 1983 Puu'Oo eruption. Overall, a single stress orientation with the maximum compressive stress oriented SE perpendicular to the rift and dipping at 45° is compatible with the coeval existence of decollement, reverse, and normal faults. However, in a crustal volume east of longitude 155°10′W, we find a change of the orientation of σ 1 from nearly horizontal to plunging 45° SE occurring in 1979. This stress rotation suggests magma movements within the aseismic part of Kilauea's east rift zone. The strain and stress orientations are coaxial in the south flank except within the volume where the stress rotation is observed. We observe a change in the relationship between stress and strain directions caused either by the shifting of seismic activity from reverse faults to decollements, while stress stays constant, or by a rotation of stress, while strain remains constant. Assuming that the model of a noncohesive Coulomb wedge is appropriate for Kilauea's south flank, we find that high pore pressures are prevalent along the decollement and within the wedge for a coefficient of friction equal to 0.85.

Hawaii↗

USGS2012

Effects of compression direction on the plasticity and rheology of hydrolytically weakened synthetic quartz crystals at atmospheric pressure

A hydrothermally grown synthetic quartz crystal with 370±60 ppm hydroxyl impurity was cut into right rectangular prisms in eight crystallographic orientations. We compressed the prisms under constant axial force corresponding to a uniaxial stress of 140.0±0.5 MPa, and temperatures of 510° and 750°C. All but one of the samples sustained permanent axial strains of 2–3%. We established the operating slip systems from specimen shape change, slip bands and dislocation etch pits on polished surfaces, crystallographic orientation changes, stress optical features in thin sections, and transmission electron microscopy. The observed creep behavior and plasticity divided the samples into three groups: (1) Crystals compressed at 45° to [0001] and [ ] and those compressed ⊥( ) and ⊥( ) deformed principally by slip parallel to [0001]. Creep rates were relatively high and were not strongly sensitive to test temperature. Dislocation arrays approximately parallel to ( ) are common. Dislocation loops are elongate parallel to [0001], indicating that the edge segments were more mobile than the screw segments. (2) The second group of samples were loaded normal to [0001] in three orientations: ⊥( ), ⊥( ), and at 45° to ( ). These samples deformed primarily by { } 〈a〉 slip with some evidence for secondary slip on the other systems. They were more creep resistant than the first group and displayed a much higher sensitivity of creep rate to test temperature. Dislocation loops are very elongate parallel to [0001], indicating that the screw dislocation segments were much more mobile than the edge segments. (3) A sample compressed parallel to [0001] at 750°C crept at a barely detectable rate (∼10 s) and no optical scale slip features were observed. These results confirm our earlier work on one orientation each from groups 1 and 2, which indicated a strong creep anisotropy for this same crystal. This creep anisotropy parallels a remarkably similar anisotropy in the diffusivity of impurities in quartz, suggesting a causal relationship between impurity diffusion and creep associated with hydrolytic weakening. Appendix Tables A1‐A5 are available with entire article on microfiche. Order from American Geophysical Union, 2000 Florida Avenue, N.W., Washington, DC 20009. Document B84‐004; $2.50. Payment must accompany order.

Journal of Geophysical Research B: Solid Earth↗

USGS1997

Continuous monitoring of surface deformation at Long Valley Caldera, California, with GPS

Continuous Global Positioning System (GPS) measurements at Long Valley Caldera, an active volcanic region in east central California, have been made on the south side of the resurgent dome since early 1993. A site on the north side of the dome was added in late 1994. Special adaptations for autonomous operation in remote regions and enhanced vertical precision were made. The data record ongoing volcanic deformation consistent with uplift and expansion of the surface above a shallow magma chamber. Measurement precisions (1 standard error) for “absolute” position coordinates, i.e., relative to a global reference frame, are 3–4 mm (north), 5–6 mm (east), and 10–12 mm (vertical) using 24 hour solutions. Corresponding velocity uncertainties for a 12 month period are about 2 mm/yr in the horizontal components and 3–4 mm/yr in the vertical component. High precision can also be achieved for relative position coordinates on short (less than 10 km) baselines using broadcast ephemerides and observing times as short as 3 hours, even when data are processed rapidly on site. Comparison of baseline length changes across the resurgent dome between the two GPS sites and corresponding two-color electronic distance measurements indicates similar extension rates within error (∼2 mm/yr) once we account for a random walk noise component in both systems that may reflect spurious monument motion. Both data sets suggest a pause in deformation for a 3.5 month period in mid-1995, when the extension rate across the dome decreased essentially to zero. Three dimensional positioning data from the two GPS stations suggest a depth (5.8±1.6 km) and location (west side of the resurgent dome) of a major inflation center, in agreement with other geodetic techniques, near the top of a magma chamber inferred from seismic data. GPS systems similar to those installed at Long Valley can provide a practical method for near real-time monitoring and hazard assessment on many active volcanoes.

Journal of Geophysical Research B: Solid Earth↗

USGS2011

Spatiotemporal evolution of dike opening and décollement slip at Kīlauea Volcano, Hawai'i

Rapid changes in ground tilt and GPS positions on Kīlauea Volcano, Hawai'i, are interpreted as resulting from a shallow, two‐segment dike intrusion into the east rift zone that began at 1217 UTC (0217 HST) on 17 June 2007 and lasted almost 3 days. As a result of the intrusion, a very small volume of basalt (about 1500 m 3 ) erupted on 19 June. Northward tilt at a coastal tiltmeter, subsidence of south flank GPS sites, southeastward displacements at southwestern flank GPS sites, and a swarm of flank earthquakes suggest that a slow slip event occurred on the décollement beneath Kīlauea's south flank concurrent with the rift intrusion. We use 4 min GPS positions that include estimates of time‐dependent tropospheric gradients and ground tilt data to study the spatial and temporal relationships between the two inferred shallow, steeply dipping dike segments extending from the surface to about 2 km depth and décollement slip at 8 km depth. We invert for the temporal evolution of distributed dike opening and décollement slip in independent inversions at each time step using a nonnegative least squares algorithm. On the basis of these inversions, the intrusion occurred in two stages that correspond spatially and temporally with concentrated rift zone seismicity. The dike opening began on the western of the two segments before jumping to the eastern segment, where the majority of opening accumulated. Dike opening preceded the start of décollement slip at an 84% confidence level; the latter is indicated by the onset of northward tilt of a coastal tiltmeter. Displacements at southwest flank GPS sites began about 18 h later and are interpreted as resulting from slow slip on the southwestern flank. Additional constraints on the evolution of the intrusion and décollement slip come from inversion of an Envisat interferogram that spans the intrusion until 0822 UTC on 18 June 2007, combined with GPS and tilt data. This inversion shows that up to 0822 UTC on 18 June, décollement slip is only required in a limited region offshore of Ka'ena Point. A similar inversion of the complete event, which includes GPS and tilt data up to 21 June and a second Envisat interferogram spanning the complete intrusion until 21 June, shows décollement slip spread westward across the south flank. This may suggest westward migration of the décollement slip as the event progressed.

Hawaii↗

USGS2019

Locating shallow seismic sources with waves scattered by surface topography: Validation of the method at the Nevada Test Site

Accurate and robust source location is fundamental in seismology. Previously, we proposed a new full waveform location method using waves scattered by topography near the source, and we demonstrated its potential in obtaining accurate source location using synthetic data. In the work presented here, we validate this method with observed seismic data from the 1993 Non-Proliferation Experiment, a chemical explosion with a precisely known location in a region with moderate topography at the Nevada Test Site. We select the first arrivals (Pn/Pg) and their immediate codas to determine the source location and compare our solution with the known location. We use a collocated-grid finite-difference method to calculate the strain Green’s tensor in a grid-search volume containing the source, and obtain three-component synthetic waveforms at 12 broadband seismic stations at regional distances using source-receiver reciprocity. We assess the solution using a least-squares misfit between the observed and synthetic waveforms. When scattered coda waves are incorporated, the best solution is within a few hundreds of meters of the exact source location, and the estimated uncertainty of the solution is reduced compared to the waveform based solution using only the P waves. The solution is robust to the choice of the frequency content and to the addition of random velocity heterogeneity. We conclude that the full waveform source location method is effective for real seismic data, though more validation tests like this one are needed to further understand its efficacy for source-station geometry, roughness of topography, signal quality and other factors.

Nevada↗

USGS2023

High-resolution InSAR reveals localized pre-eruptive deformation inside the crater of Agung Volcano, Indonesia

During a volcanic crisis, high-rate, localized deformation can indicate magma close to the surface, with important implications for eruption forecasting. However, only a few such examples have been reported, because frequent, dense monitoring is needed. High-resolution Synthetic Aperture Radar (SAR) is capable of achieving <1 m spatial resolution and sub-weekly revisit times, but is under-used. Here we use high-resolution satellite SAR imagery from COSMO-SkyMed, TerraSAR-X, and Sentinel-1 to detect intra-crater uplift preceding the November 2017 onset of eruptive activity at Agung, Indonesia. Processing the SAR imagery with an up-to-date, accurate, high-resolution digital elevation model was crucial for preventing aliasing of the deformation signal and for accurate georeferencing. We show that >15 cm of line-of-sight shortening occurred over a 400-by-400 m area on the crater floor in September-October 2017, accompanying a deep seismic swarm and flank dyke intrusion. We attribute the deformation to the pressurization of a shallow (<200 m deep) hydrothermal system by the injection of magmatic gases and fluids. We also observe a second pulse of intra-crater deformation of 3–5 cm within 4 days to 11 hr prior to the first phreatomagmatic eruption, which is consistent with interaction between the hydrothermal system and the ascending magma. This phreatomagmatic eruption created the central pathway used during the final stages of magma ascent. Our observations have important implications for understanding unrest and eruption forecasting, and demonstrate the potential of monitoring with high-resolution SAR.

Agung Volcano, Bali↗

USGS1999

Abrupt along-strike change in tectonic style: San Andreas fault zone, San Francisco Peninsula

Seismicity and high-resolution aeromagnetic data are used to define an abrupt change from compressional to extensional tectonism within a 10- to 15-km-wide zone along the San Andreas fault on the San Francisco Peninsula and offshore from the Golden Gate. This 100-km-long section of the San Andreas fault includes the hypocenter of the M w = 7.8 1906 San Francisco earthquake as well as the highest level of persistent microseismicity along that ∼470-km-long rupture. We define two distinct zones of deformation along this stretch of the fault using well-constrained relocations of all post-1969 earthquakes based a joint one-dimensional velocity/hypocenter inversion and a redetermination of focal mechanisms. The southern zone is characterized by thrust- and reverse-faulting focal mechanisms with NE trending P axes that indicate “fault-normal” compression in 7- to 10-km-wide zones of deformation on both sides of the San Andreas fault. A 1- to 2-km-wide vertical zone beneath the surface trace of the San Andreas is characterized by its almost complete lack of seismicity. The compressional deformation is consistent with the young, high topography of the Santa Cruz Mountains/Coast Ranges as the San Andreas fault makes a broad restraining left bend (∼10°) through the southernmost peninsula. A zone of seismic quiescence ∼15 km long separates this compressional zone to the south from a zone of combined normal-faulting and strike-slip-faulting focal mechanisms (including a M L = 5.3 earthquake in 1957) on the northernmost peninsula and offshore on the Golden Gate platform. Both linear pseudogravity gradients, calculated from the aeromagnetic data, and seismic reflection data indicate that the San Andreas fault makes an abrupt ∼3-km right step less than 5 km offshore in this northern zone. A similar right-stepping (dilatational) geometry is also observed for the subparallel San Gregorio fault offshore. Persistent seismicity and extensional tectonism occur within the San Andreas right stepover region and at least 15 km along-strike both to the SE and NW. The 1906 San Francisco earthquake may have nucleated within the San Andreas right stepover, which may help explain the bilateral nature of rupture of this event. Our analysis suggests two seismic hazards for the San Francisco Peninsula in addition to the hazard associated with a M = 7 to 8 strike-slip earthquake along the San Andreas fault: the potential for a M ≈ 6 normal-faulting earthquake just 5–8 km west of San Francisco and a M = 6+ thrust faulting event in the southern peninsula.

Journal of Geophysical Research B: Solid Earth↗

USGS2003

Groundwater flow, heat transport, and water table position within volcanic edifices: Implications for volcanic processes in the Cascade Range

The position of the water table within a volcanic edifice has significant implications for volcano hazards, geothermal energy, and epithermal mineralization. We have modified the HYDROTHERM numerical simulator to allow for a free-surface (water table) upper boundary condition and a wide range of recharge rates, heat input rates, and thermodynamic conditions representative of continental volcano-hydrothermal systems. An extensive set of simulations was performed on a hypothetical stratovolcano system with unconfined groundwater flow. Simulation results suggest that the permeability structure of the volcanic edifice and underlying material is the dominant control on water table elevation and the distribution of pressures, temperatures, and fluid phases at depth. When permeabilities are isotropic, water table elevation decreases with increasing heat flux and increases with increasing recharge, but when permeabilities are anisotropic, these effects can be much less pronounced. Several conditions facilitate the ascent of a hydrothermal plume into a volcanic edifice: a sufficient source of heat and magmatic volatiles at depth, strong buoyancy forces, and a relatively weak topography-driven flow system. Further, the plume must be connected to a deep heat source through a pathway with a time-averaged effective permeability ???1 ?? 10-16 m2, which may be maintained by frequent seismicity. Topography-driven flow may be retarded by low permeability in the edifice and/or the lack of precipitation recharge; in the latter case, the water table may be relatively deep. Simulation results were compared with observations from the Quaternary stratovolcanoes along the Cascade Range of the western United States to infer hydrothermal processes within the edifices. Extensive ice caps on many Cascade Range stratovolcanoes may restrict recharge on the summits and uppermost flanks. Both the simulation results and limited observational data allow for the possibility that the water table beneath the stratovolcanoes is relatively deep.

Journal of Geophysical Research B: Solid Earth↗

USGS2001

Rapid fluid disruption: A source for self-potential anomalies on volcanoes

Self-potential (SP) anomalies observed above suspected magma reservoirs, dikes, etc., on various volcanoes (Kilauea, Hawaii; Mount Unzen, Japan; Piton de la Fournaise, Reunion Island, Miyake Jima, Japan) result from transient surface electric fields of tens of millivolts per kilometer and generally have a positive polarity. These SP anomalies are usually attributed to electrokinetic effects where properties controlling this process are poorly constrained. We propose an alternate explanation that contributions to electric fields of correct polarity should be expected from charge generation by fluid vaporization/disruption. As liquids are vaporized or removed as droplets by gas transport away from hot dike intrusions, both charge generation and local increase in electrical resistivity by removal of fluids should occur. We report laboratory observations of electric fields in hot rock samples generated by pulses of fluid (water) through the rock at atmospheric pressure. These indicate the relative amplitudes of rapid fluid disruption (RFD) potentials and electrokinetic potentials to be dramatically different and the signals are opposite in sign. Above vaporization temperatures, RFD effects of positive sign in the direction of gas flow dominate, whereas below these temperatures, effects of negative sign dominate. This suggests that the primary contribution to observed self-potential anomalies arises from gas-related charge transport processes at temperatures high enough to produce vigorous boiling and vapor transport. At lower temperatures, the primary contribution is from electrokinetic effects modulated perhaps by changing electrical resistivity and RFD effects from high-pressure but low-temperature CO2 and SO2 gas flow ripping water molecules from saturated crustal rocks. If charge generation is continuous, as could well occur above a newly emplaced dike, positive static potentials will be set up that could be sustained for many years, and the simplest method for identifying these hot, active regions would be to identify the SP anomalies they generate.

Hawai'i↗

USGS2022

Orbital and in-situ investigation of periodic bedrock ridges in Glen Torridon, Gale Crater, Mars

Wind has been the dominant agent of landscape modification on Mars for the past ~3 billion years. Among the diversity of features formed by aeolian abrasion on the surface of Mars are periodic bedrock ridges (PBRs), a relatively recently recognized class of erosional bedforms on Mars for which Earth analogues are rare. Gale crater, the field site for NASA’s Mars Science Laboratory Curiosity rover since it landed there in 2012, contains a diverse and extensive record of aeolian deposition and erosion. This study focuses on a series of periodic, linear bedrock ridges that occur within the Fe/Mg-smectite clay-bearing Glen Torridon region of Aeolis Mons (informally Mount Sharp). During Curiosity’s exploration of the Glen Torridon region between sols ~2300-3080, the rover drove through this field of ridges, providing the first opportunity for the in situ observation of these enigmatic erosional features. This study characterizes the Glen Torridon ridges using orbiter and rover data to determine their morphology, spatial distribution, compositional and material properties, and association with other aeolian features in the area. Based on these observations, the Glen Torridon ridges are interpreted to be consistent with an origin as wind-eroded periodic bedrock ridges carved during the most recent exhumation of Mount Sharp into the present-day mound. Although there is evidence for multidirectional winds in the Glen Torridon region based on the orientation of modern ripples, megaripples, TARs and other bedrock indicators, the consistent orientation of the Glen Torridon ridges, coupled with morphologic asymmetries within the ridges, support formation and elongation of the Glen Torridon PBRs forms parallel to a net regional northerly wind direction in and around Gale crater.

Journal of Geophysical Research E: Planets↗