SEARCH · Geology Reports

Results for “Geophysical Journal”

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

USGS2014

Seismicity, faulting, and structure of the Koyna-Warna seismic region, Western India from local earthquake tomography and hypocenter locations

Although seismicity near Koyna Reservoir (India) has persisted for ~50 years and includes the largest induced earthquake ( M 6.3) reported worldwide, the seismotectonic framework of the area is not well understood. We recorded ~1800 earthquakes from 6 January 2010 to 28 May 2010 and located a subset of 343 of the highest-quality earthquakes using the tomoDD code of Zhang and Thurber (2003) to better understand the framework. We also inverted first arrivals for 3-D Vp , Vs , and Vp / Vs and Poisson's ratio tomography models of the upper 12 km of the crust. Epicenters for the recorded earthquakes are located south of the Koyna River, including a high-density cluster that coincides with a shallow depth (<1.5 km) zone of relatively high Vp and low Vs (also high Vp / Vs and Poisson's ratios) near Warna Reservoir. This anomalous zone, which extends near vertically to at least 8 km depth and laterally northward at least 15 km, is likely a water-saturated zone of faults under high pore pressures. Because many of the earthquakes occur on the periphery of the fault zone, rather than near its center, the observed seismicity-velocity correlations are consistent with the concept that many of the earthquakes nucleate in fractures adjacent to the main fault zone due to high pore pressure. We interpret our velocity images as showing a series of northwest trending faults locally near the central part of Warna Reservoir and a major northward trending fault zone north of Warna Reservoir.

Maharashtra↗

USGS2014

Slip rates and spatially variable creep on faults of the northern San Andreas system inferred through Bayesian inversion of Global Positioning System data

Fault creep, depending on its rate and spatial extent, is thought to reduce earthquake hazard by releasing tectonic strain aseismically. We use Bayesian inversion and a newly expanded GPS data set to infer the deep slip rates below assigned locking depths on the San Andreas, Maacama, and Bartlett Springs Faults of Northern California and, for the latter two, the spatially variable interseismic creep rate above the locking depth. We estimate deep slip rates of 21.5 ± 0.5, 13.1 ± 0.8, and 7.5 ± 0.7 mm/yr below 16 km, 9 km, and 13 km on the San Andreas, Maacama, and Bartlett Springs Faults, respectively. We infer that on average the Bartlett Springs fault creeps from the Earth's surface to 13 km depth, and below 5 km the creep rate approaches the deep slip rate. This implies that microseismicity may extend below the locking depth; however, we cannot rule out the presence of locked patches in the seismogenic zone that could generate moderate earthquakes. Our estimated Maacama creep rate, while comparable to the inferred deep slip rate at the Earth's surface, decreases with depth, implying a slip deficit exists. The Maacama deep slip rate estimate, 13.1 mm/yr, exceeds long-term geologic slip rate estimates, perhaps due to distributed off-fault strain or the presence of multiple active fault strands. While our creep rate estimates are relatively insensitive to choice of model locking depth, insufficient independent information regarding locking depths is a source of epistemic uncertainty that impacts deep slip rate estimates.

California↗

USGS2014

Latest Quaternary paleoseismology and evidence of distributed dextral shear along the Mohawk Valley fault zone, northern Walker Lane, California

The dextral-slip Mohawk Valley fault zone (MVFZ) strikes northwestward along the eastern margin of the Sierra Nevada in the northern Walker Lane. Geodetic block modeling indicates that the MVFZ may accommodate ~3 mm/yr of regional dextral strain, implying that it is the highest slip-rate strike-slip fault in the region; however, only limited geologic data are available to constrain the system’s slip rate and earthquake history. We mapped the MVFZ using airborne lidar data and field observations and identified a site near Sulphur Creek for paleoseismic investigation. At this site, oblique dextral-normal faulting on the steep valley margin has created a closed depression that floods annually during spring snowmelt to form an ephemeral pond. We excavated three fault-perpendicular trenches at the site and exposed pond sediment that interfingers with multiple colluvial packages eroded from the scarp that bounds the eastern side of the pond. We documented evidence for four surface-rupturing earthquakes on this strand of the MVFZ. OxCal modeling of radiocarbon and luminescence ages indicates that these earthquakes occurred at 14.0 ka, 12.8 ka, 5.7 ka, and 1.9 ka. The mean ~4 kyr recurrence interval is inconsistent with slip rates of ~3 mm/yr; these rates imply surface ruptures of more than 10 m per event, which is geologically implausible for the subdued geomorphic expression and 60 km length of the MVFZ. We propose that unidentified structures not yet incorporated into geodetic models may accommodate significant dextral shear across the northern Walker Lane, highlighting the role of distributed deformation in this region.

California↗

USGS2014

Frequency-dependent effects of rupture for the 2004 Parkfield mainshock, results from UPSAR

The frequency-dependent effects of rupture propagation of the Parkfield, California earthquake (Sept. 28, 2004, M6) to the northwest along the San Andreas fault can be seen in acceleration records at UPSAR (USGS Parkfield Seismic Array) in at least two ways. First, we can see the effects of directivity in the acceleration traces at UPSAR, which is about 11.5 km from the epicenter. Directivity or the seismic equivalent of a Doppler shift has been documented in many cases by comparing short duration, high-amplitude pulses (P or S) in the forward direction with longer duration body waves in the backward direction. In this case we detect a change from a relatively large amplitude, coherent, high-frequency signal at the start of rupture to a low-amplitude, low-coherent, low-frequency signal at about the time the rupture front transfers from the forward azimuth to the back azimuth at about 34-36 s (time is UTC and are the seconds after day 272 and 17 hours and 15 minutes. S arrival is just after 30s) for rays leaving the fault and propagating to UPSAR. The frequency change is obvious in the band about 5 to 30 Hz, which is significantly above the corner frequency of the earthquake (about 0.11Hz). From kinematic source models, the duration of faulting is about 9.2 s and the change in frequency is during faulting as the rupture extends to the northwest. Understanding the systematic change in frequency and amplitude of seismic waves in relation to the propagation of the rupture front is important for predicting strong ground motion. Second, we can filter the acceleration records from the array to determine if the low frequency energy emerges from the same part of the fault as the high frequency signal (e.g. has the same back azimuth and apparent velocity at UPSAR) an important clue to the dynamics of rupture. Analysis of sources of strong motion (characterized by relatively high frequencies) compared to kinematic slip models (relatively low frequency) for the March 11, 2011 Tohoku earthquake as well as Maule (Feb. 27, 2010) and Chi-Chi (Sept. 20, 1999) earthquakes show that high- and low-frequency sources do not have the same locations on the fault. In this paper we filter the accelerograms from UPSAR for the 2004 mainshock in various passbands and then re-compute the cross correlations to determine the vector slowness of the incoming waves. At Parkfield, it appears that for seismic waves with frequencies above 1 Hz there is no discernible frequency-dependent difference in source position (up to 8 Hz) based on estimates of back azimuth and apparent velocity. However at lower frequencies, sources appear to be from shallower depths and trail the high frequencies as the rupture proceeds down the fault. This result is greater than one standard deviation of an estimate of error, based on a new method of estimating error that is a measure of how broad the peak in correlation is and an estimate of the variance of the correlation values. These observations can be understood in terms of a rupture front that is more energetic and coherent near the front of rupture (radiating higher frequencies) and less coherent and less energetic (radiating in a lower frequency band) behind the initial rupture front. This result is a qualitative assessment of changes in azimuth and apparent velocity with frequency and time and does not include corrections to find the source location on the fault.

California↗

USGS2014

High-frequency imaging of elastic contrast and contact area with implications for naturally observed changes in fault properties

During localized slip of a laboratory fault we simultaneously measure the contact area and the dynamic fault normal elastic stiffness. One objective is to determine conditions where stiffness may be used to infer changes in area of contact during sliding on nontransparent fault surfaces. Slip speeds between 0.01 and 10 µm/s and normal stresses between 1 and 2.5 MPa were imposed during velocity step, normal stress step, and slide-hold-slide tests. Stiffness and contact area have a linear interdependence during rate stepping tests and during the hold portion of slide-hold-slide tests. So long as linearity holds, measured fault stiffness can be used on nontransparent materials to infer changes in contact area. However, there are conditions where relations between contact area and stiffness are nonlinear and nonunique. A second objective is to make comparisons between the laboratory- and field-measured changes in fault properties. Time-dependent changes in fault zone normal stiffness made in stress relaxation tests imply postseismic wave speed changes on the order of 0.3% to 0.8% per year in the two or more years following an earthquake; these are smaller than postseismic increases seen within natural damage zones. Based on scaling of the experimental observations, natural postseismic fault normal contraction could be accommodated within a few decimeter wide fault core. Changes in the stiffness of laboratory shear zones exceed 10% per decade and might be detectable in the field postseismically.

Journal of Geophysical Research↗

USGS2014

Continuous uplift near the seaward edge of the Prince William Sound megathrust: Middleton Island, Alaska

Middleton Island, located at the seaward edge of the continental shelf 50 km from the base of the inner wall of the Aleutian Trench, affords an opportunity to make land-based measurements of uplift near the toe of the Prince William Sound megathrust, site of the 1964, M  = 9.2, Alaska earthquake. Leveling surveys (1973–1993) on Middleton Island indicate roughly uniform tilting (~1 µrad/a down to the northwest) of the island, and GPS surveys (1993–2012) show an uplift rate of 14 mm/a of the island relative to fixed North America. The data are consistent with a combined (coseismic and postseismic) uplift (in meters) due to the 1964 earthquake as a function of time τ (years after the earthquake) u ( τ ) = (3.5 + 1.21 log 10  [1 + 1.67  τ ]) H ( τ ) where 3.5 is the coseismic uplift and H ( τ ) is 0 for τ  < 0 and 1 otherwise. The current uplift on Middleton Island is attributed to continuous slip on a fault splaying off from the megathrust, and the long-term uplift is the superposition of the effects of past earthquakes, each earthquake being similar to the 1964 event. Then, the predicted uplift at time t due to a sequence of earthquakes at times t i would be . From studies of strandlines associated with the uplifted terraces on Middleton Island, Plafker et al. (1992) estimated the occurrence times of the last six earthquakes and measured the present-day elevations of those strandlines. The predicted uplift is in rough agreement with those measurements. About half of the predicted uplift is due to postseismic relaxation from previous earthquakes.

Alaska↗

USGS2014

Time-averaged discharge rate of subaerial lava at Kīlauea Volcano, Hawai‘i, measured from TanDEM-X interferometry: Implications for magma supply and storage during 2011-2013

Differencing digital elevation models (DEMs) derived from TerraSAR add-on for Digital Elevation Measurements (TanDEM-X) synthetic aperture radar imagery provides a measurement of elevation change over time. On the East Rift Zone (EZR) of Kīlauea Volcano, Hawai‘i, the effusion of lava causes changes in topography. When these elevation changes are summed over the area of an active lava flow, it is possible to quantify the volume of lava emplaced at the surface during the time spanned by the TanDEM-X data—a parameter that can be difficult to measure across the entirety of an ~100 km 2 lava flow field using ground-based techniques or optical remote sensing data. Based on the differences between multiple TanDEM-X-derived DEMs collected days to weeks apart, the mean dense-rock equivalent time-averaged discharge rate of lava at Kīlauea between mid-2011 and mid-2013 was approximately 2 m 3 /s, which is about half the long-term average rate over the course of Kīlauea's 1983–present ERZ eruption. This result implies that there was an increase in the proportion of lava stored versus erupted, a decrease in the rate of magma supply to the volcano, or some combination of both during this time period. In addition to constraining the time-averaged discharge rate of lava and the rates of magma supply and storage, topographic change maps derived from space-based TanDEM-X data provide insights into the four-dimensional evolution of Kīlauea's ERZ lava flow field. TanDEM-X data are a valuable complement to other space-, air-, and ground-based observations of eruptive activity at Kīlauea and offer great promise at locations around the world for aiding with monitoring not just volcanic eruptions but any hazardous activity that results in surface change, including landslides, floods, earthquakes, and other natural and anthropogenic processes.

Hawai'i↗

USGS2014

The 2010 slow slip event and secular motion at Kilauea, Hawai`i inferred from TerraSAR-X InSAR data

We present here an Small BAseline Subset (SBAS) algorithm to extract both transient and secular ground deformations on the order of millimeters in the presence of tropospheric noise on the order of centimeters, when the transient is of short duration and known time, and the background deformation is smooth in time. We applied this algorithm to study the 2010 slow slip event as well as the secular motion of Kīlauea's south flank using 49 TerraSAR-X images. We also estimate the tropospheric delay variation relative to a given reference pixel using an InSAR SBAS approach. We compare the InSAR SBAS solution for both ground deformation and tropospheric delays with existing GPS measurements and confirm that the ground deformation signal andtropospheric noise in InSAR data are successfully separated. We observe that the coastal region on the south side of the Hilina Pali moves at a higher background rate than the region north side of the Pali. We also conclude that the 2010 SSE displacement is mainly horizontal and the maximum magnitude of the 2010 SSE vertical component is less than 5 mm.

Hawaii↗

USGS2014

Laboratory constraints on models of earthquake recurrence

In this study, rock friction ‘stick-slip’ experiments are used to develop constraints on models of earthquake recurrence. Constant-rate loading of bare rock surfaces in high quality experiments produces stick-slip recurrence that is periodic at least to second order. When the loading rate is varied, recurrence is approximately inversely proportional to loading rate. These laboratory events initiate due to a slip rate-dependent process that also determines the size of the stress drop [Dieterich, 1979; Ruina, 1983] and as a consequence, stress drop varies weakly but systematically with loading rate [e.g., Gu and Wong, 1991; Karner and Marone, 2000; McLaskey et al., 2012]. This is especially evident in experiments where the loading rate is changed by orders of magnitude, as is thought to be the loading condition of naturally occurring, small repeating earthquakes driven by afterslip, or low-frequency earthquakes loaded by episodic slip. As follows from the previous studies referred to above, experimentally observed stress drops are well described by a logarithmic dependence on recurrence interval that can be cast as a non-linear slip-predictable model. The fault’s rate dependence of strength is the key physical parameter. Additionally, even at constant loading rate the most reproducible laboratory recurrence is not exactly periodic, unlike existing friction recurrence models. We present example laboratory catalogs that document the variance and show that in large catalogs, even at constant loading rate, stress drop and recurrence co-vary systematically. The origin of this covariance is largely consistent with variability of the dependence of fault strength on slip rate. Laboratory catalogs show aspects of both slip and time predictability and successive stress drops are strongly correlated indicating a ‘memory’ of prior slip history that extends over at least one recurrence cycle.

Journal of Geophysical Research↗

USGS2014

Preslip and cascade processes initiating laboratory stick slip

Recent modeling studies have explored whether earthquakes begin with a large aseismic nucleation process or initiate dynamically from the rapid growth of a smaller instability in a “cascade-up” process. To explore such a case in the laboratory, we study the initiation of dynamic rupture (stick slip) of a smooth saw-cut fault in a 76mm diameter cylindrical granite laboratory sample at 40–120MPa confining pressure. We use a high dynamic range recording system to directly compare the seismic waves radiated during the stick-slip event to those radiated from tiny (M _6) discrete seismic events, commonly known as acoustic emissions (AEs), that occur in the seconds prior to each large stick slip. The seismic moments, focal mechanisms, locations, and timing of the AEs all contribute to our understanding of their mechanics and provide us with information about the stick-slip nucleation process. In a sequence of 10 stick slips, the first few microseconds of the signals recorded from stick-slip instabilities are nearly indistinguishable from those of premonitory AEs. In this sense, it appears that each stick slip begins as an AE event that rapidly (~20 μs) grows about 2 orders of magnitude in linear dimension and ruptures the entire 150mm length of the simulated fault. We also measure accelerating fault slip in the final seconds before stick slip. We estimate that this slip is at least 98% aseismic and that it both weakens the fault and produces AEs that will eventually cascade-up to initiate the larger dynamic rupture.

Journal of Geophysical Research B: Solid Earth↗

USGS2014

Seismological and geodetic constraints on the 2011 M w 5.3 Trinidad, Colorado earthquake and induced deformation in the Raton Basin

The Raton Basin of southern Colorado and northern New Mexico is an actively produced hydrocarbon basin that has experienced increased seismicity since 2001, including the August 2011 M w 5.3 Trinidad normal faulting event. Following the 2011 earthquake, regional seismic observations were used to relocate 21 events, including the 2011 main shock, two foreshocks, and 13 aftershocks. Additionally, interferometric synthetic aperture radar (InSAR) observations of both the 2011 event and preevent basin deformation place constraint on the spatial kinematics of the 2011 event and localized basin subsidence due to ground water or gas withdrawal. We find that the 2011 earthquake ruptured an 8–10 km long segment of a normal fault at depths of 1.5–6.0 km within the crystalline Precambrian basement underlying the Raton Basin sedimentary rocks. The earthquake also nucleated within the crystalline basement in the vicinity of an active wastewater disposal site. The ensuing aftershock sequence demonstrated statistical properties expected for intraplate earthquakes, though the length of the 2011 earthquake is unexpectedly long for an M w 5.3 event, suggesting that wastewater disposal may have triggered a low stress drop, otherwise natural earthquake. Additionally, preevent and postevent seismicity in the Raton Basin spatially correlates to regions of subsidence observed in InSAR time series analysis. While these observations cannot discern a causal link between hydrocarbon production and seismicity, they constrain spatial relationships between active basin deformation and geological and anthropogenic features. Furthermore, the InSAR observations highlight the utility of space-based geodetic observations for monitoring and assessing anthropogenically induced and triggered deformation.

Journal of Geophysical Research B: Solid Earth↗

USGS2014

Widespread gas hydrate instability on the upper U.S. Beaufort margin

The most climate-sensitive methane hydrate deposits occur on upper continental slopes at depths close to the minimum pressure and maximum temperature for gas hydrate stability. At these water depths, small perturbations in intermediate ocean water temperatures can lead to gas hydrate dissociation. The Arctic Ocean has experienced more dramatic warming than lower latitudes, but observational data have not been used to study the interplay between upper slope gas hydrates and warming ocean waters. Here we use (a) legacy seismic data that constrain upper slope gas hydrate distributions on the U.S. Beaufort Sea margin, (b) Alaskan North Slope borehole data and offshore thermal gradients determined from gas hydrate stability zone thickness to infer regional heat flow, and (c) 1088 direct measurements to characterize multidecadal intermediate ocean warming in the U.S. Beaufort Sea. Combining these data with a three-dimensional thermal model shows that the observed gas hydrate stability zone is too deep by 100 to 250 m. The disparity can be partially attributed to several processes, but the most important is the reequilibration (thinning) of gas hydrates in response to significant (~0.5°C at 2σ certainty) warming of intermediate ocean temperatures over 39 years in a depth range that brackets the upper slope extent of the gas hydrate stability zone. Even in the absence of additional ocean warming, 0.44 to 2.2 Gt of methane could be released from reequilibrating gas hydrates into the sediments underlying an area of ~5–7.5 × 10 3 km 2 on the U.S. Beaufort Sea upper slope during the next century.

Alaska↗

USGS2014

Attenuation and scattering tomography of the deep plumbing system of Mount St. Helens

We present a combined 3-D P wave attenuation, 2-D S coda attenuation, and 3-D S coda scattering tomography model of fluid pathways, feeding systems, and sediments below Mount St. Helens (MSH) volcano between depths of 0 and 18 km. High-scattering and high-attenuation shallow anomalies are indicative of magma and fluid-rich zones within and below the volcanic edifice down to 6 km depth, where a high-scattering body outlines the top of deeper aseismic velocity anomalies. Both the volcanic edifice and these structures induce a combination of strong scattering and attenuation on any seismic wavefield, particularly those recorded on the northern and eastern flanks of the volcanic cone. North of the cone between depths of 0 and 10 km, a low-velocity, high-scattering, and high-attenuation north-south trending trough is attributed to thick piles of Tertiary marine sediments within the St. Helens Seismic Zone. A laterally extended 3-D scattering contrast at depths of 10 to 14 km is related to the boundary between upper and lower crust and caused in our interpretation by the large-scale interaction of the Siletz terrane with the Cascade arc crust. This contrast presents a low-scattering, 4–6 km 2 “hole” under the northeastern flank of the volcano. We infer that this section represents the main path of magma ascent from depths greater than 6 km at MSH, with a small north-east shift in the lower plumbing system of the volcano. We conclude that combinations of different nonstandard tomographic methods, leading toward full-waveform tomography, represent the future of seismic volcano imaging.

Washington↗

USGS2014

Gravity changes and deformation at Kīlauea Volcano, Hawaii, associated with summit eruptive activity, 2009-2012

Analysis of microgravity and surface displacement data collected at the summit of Kīlauea Volcano, Hawaii (USA), between December 2009 and November 2012 suggests a net mass accumulation at ~1.5 km depth beneath the northeast margin of Halema‘uma‘u Crater, within Kīlauea Caldera. Although residual gravity increases and decreases are accompanied by periods of uplift and subsidence of the surface, respectively, the volume change inferred from the modeling of interferometric synthetic aperture radar deformation data can account for only a small portion (as low as 8%) of the mass addition responsible for the gravity increase. We propose that since the opening of a new eruptive vent at the summit of Kīlauea in 2008, magma rising to the surface of the lava lake outgasses, becomes denser, and sinks to deeper levels, replacing less dense gas-rich magma stored in the Halema‘uma‘u magma reservoir. In fact, a relatively small density increase (<200 kg m −3 ) of a portion of the reservoir can produce the positive residual gravity change measured during the period with the largest mass increase, between March 2011 and November 2012. Other mechanisms may also play a role in the gravity increase without producing significant uplift of the surface, including compressibility of magma, formation of olivine cumulates, and filling of void space by magma. The rate of gravity increase, higher than during previous decades, varies through time and seems to be directly correlated with the volcanic activity occurring at both the summit and the east rift zone of the volcano.

Hawaii↗

USGS2014

Eruptions at Lone Star geyser, Yellowstone National Park, USA: 2. Constraints on subsurface dynamics

We use seismic, tilt, lidar, thermal, and gravity data from 32 consecutive eruption cycles of Lone Star geyser in Yellowstone National Park to identify key subsurface processes throughout the geyser's eruption cycle. Previously, we described measurements and analyses associated with the geyser's erupting jet dynamics. Here we show that seismicity is dominated by hydrothermal tremor (~5–40 Hz) attributed to the nucleation and/or collapse of vapor bubbles. Water discharge during eruption preplay triggers high-amplitude tremor pulses from a back azimuth aligned with the geyser cone, but during the rest of the eruption cycle it is shifted to the east-northeast. Moreover, ~4 min period ground surface displacements recur every 26 ± 8 min and are uncorrelated with the eruption cycle. Based on these observations, we conclude that (1) the dynamical behavior of the geyser is controlled by the thermo-mechanical coupling between the geyser conduit and a laterally offset reservoir periodically filled with a highly compressible two-phase mixture, (2) liquid and steam slugs periodically ascend into the shallow crust near the geyser system inducing detectable deformation, (3) eruptions occur when the pressure decrease associated with overflow from geyser conduit during preplay triggers an unstable feedback between vapor generation (cavitation) and mass discharge, and (4) flow choking at a constriction in the conduit arrests the runaway process and increases the saturated vapor pressure in the reservoir by a factor of ~10 during eruptions.

Wyoming↗

USGS2014

The impact of static stress change, dynamic stress change, and the background stress on aftershock focal mechanisms

The focal mechanisms of earthquakes in Southern California before and after four M  ≥ 6.7 main shocks provide insight into how fault systems respond to stress and changes in stress. The main shock static stress changes have two observed impacts on the seismicity: changing the focal mechanisms in a given location to favor those aligned with the static stress change and changing the spatial distribution of seismicity to favor locations where the static stress change aligns with the background stress. The aftershock focal mechanisms are significantly aligned with the static stress changes for absolute stress changes of ≥ 0.02 MPa, for up to ~20 years following the main shock. The dynamic stress changes have similar, although smaller, effects on the local focal mechanisms and the spatial seismicity distribution. Dynamic stress effects are best observed at long periods (30–60 s) and for metrics based on repeated stress cycling in the same direction. This implies that dynamic triggering operates, at least in part, through cyclic shear stress loading in the direction of fault slip. The background stress also strongly controls both the preshock and aftershock mechanisms. While most aftershock mechanisms are well oriented in the background stress field, 10% of aftershocks are identified as poorly oriented outliers, which may indicate limited heterogeneity in the postmain shock stress field. The fault plane orientations of the outliers are well oriented in the background stress, while their slip directions are not, implying that the background stress restricts the distribution of available fault planes.

California↗

USGS2015

Relationships among seismic velocity, metamorphism, and seismic and aseismic fault slip in the Salton Sea Geothermal Field region

The Salton Sea Geothermal Field is one of the most geothermally and seismically active areas in California and presents an opportunity to study the effect of high-temperature metamorphism on the properties of seismogenic faults. The area includes numerous active tectonic faults that have recently been imaged with active source seismic reflection and refraction. We utilize the active source surveys, along with the abundant microseismicity data from a dense borehole seismic network, to image the 3-D variations in seismic velocity in the upper 5 km of the crust. There are strong velocity variations, up to ~30%, that correlate spatially with the distribution of shallow heat flow patterns. The combination of hydrothermal circulation and high-temperature contact metamorphism has significantly altered the shallow sandstone sedimentary layers within the geothermal field to denser, more feldspathic, rock with higher P wave velocity, as is seen in the numerous exploration wells within the field. This alteration appears to have a first-order effect on the frictional stability of shallow faults. In 2005, a large earthquake swarm and deformation event occurred. Analysis of interferometric synthetic aperture radar data and earthquake relocations indicates that the shallow aseismic fault creep that occurred in 2005 was localized on the Kalin fault system that lies just outside the region of high-temperature metamorphism. In contrast, the earthquake swarm, which includes all of the M > 4 earthquakes to have occurred within the Salton Sea Geothermal Field in the last 15 years, ruptured the Main Central Fault (MCF) system that is localized in the heart of the geothermal anomaly. The background microseismicity induced by the geothermal operations is also concentrated in the high-temperature regions in the vicinity of operational wells. However, while this microseismicity occurs over a few kilometer scale region, much of it is clustered in earthquake swarms that last from hours to a few days and are localized near the MCF system.

Salton Sea↗

USGS2015

Stress and mass changes at a “wet” volcano: Example during the 2011–2012 volcanic unrest at Kawah Ijen volcano (Indonesia)

Since 2010, Kawah Ijen volcano has been equipped with seismometers, and its extremely acid volcanic lake has been monitored using temperature and leveling sensors, providing unprecedented time resolution of multiparametric data for an acidic volcanic lake. The nature of stress and mass changes of the volcano is studied by combining seismic analyses and volcanic lake measurements that were made during the strongest unrest ever recorded by the seismic network at Kawah Ijen. The distal VT earthquake swarm that occurred in May 2011 was the precursor of volcanic unrest in October 2011 that caused an increase in shallow earthquakes. The proximal VT earthquakes opened pathways for ﬂuids to ascend by increasing the permeability of the rock matrix. The following months were characterized by two periods of strong heat and mass discharge into the lake and by the initiation of monochromatic tremor (MT) activity when steam/gases interacted with shallow portions of the aquifer. Signiﬁcant seismic velocity variations, concurrent with water level rises in which water contained a large amount of steam/gas, were associated with the crises, that caused an although the unrest did not aﬀect the shallow hydrothermal system at a large scale. Whereas shallow VT earthquakes likely reﬂect a magmatic intrusion, MT and relative seismic velocity changes are clearly associated with shallow hydrothermal processes. These results will facilitate the forecast of future crises.

Kawah Ijen↗