SEARCH · Search USGS

Results for “Journal of Geophysical Research, Solid Earth”

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 595 records · Page 33Linked to original sources

USGS2011

Elevated gas hydrate saturation within silt and silty clay sediments in the Shenhu area, South China Sea

Gas hydrate saturations were estimated using five different methods in silt and silty clay foraminiferous sediments from drill hole SH2 in the South China Sea. Gas hydrate saturations derived from observed pore water chloride values in core samples range from 10 to 45% of the pore space at 190–221 m below seafloor (mbsf). Gas hydrate saturations estimated from resistivity (R t ) using wireline logging results are similar and range from 10 to 40.5% in the pore space. Gas hydrate saturations were also estimated by P wave velocity obtained during wireline logging by using a simplified three-phase equation (STPE) and effective medium theory (EMT) models. Gas hydrate saturations obtained from the STPE velocity model (41.0% maximum) are slightly higher than those calculated with the EMT velocity model (38.5% maximum). Methane analysis from a 69 cm long depressurized core from the hydrate-bearing sediment zone indicates that gas hydrate saturation is about 27.08% of the pore space at 197.5 mbsf. Results from the five methods show similar values and nearly identical trends in gas hydrate saturations above the base of the gas hydrate stability zone at depths of 190 to 221 mbsf. Gas hydrate occurs within units of clayey slit and silt containing abundant calcareous nannofossils and foraminifer, which increase the porosities of the fine-grained sediments and provide space for enhanced gas hydrate formation. In addition, gas chimneys, faults, and fractures identified from three-dimensional (3-D) and high-resolution two-dimensional (2-D) seismic data provide pathways for fluids migrating into the gas hydrate stability zone which transport methane for the formation of gas hydrate. Sedimentation and local canyon migration may contribute to higher gas hydrate saturations near the base of the stability zone.

Journal of Geophysical Research B: Solid Earth↗

USGS2011

A reevaluation of the Pallett Creek earthquake chronology based on new AMS radiocarbon dates, San Andreas fault, California

The Pallett Creek paleoseismic record occupies a keystone position in most attempts to develop rupture histories for the southern San Andreas fault. Previous estimates of earthquake ages at Pallett Creek were determined by decay counting radiocarbon methods. That method requires large samples which can lead to unaccounted sources of uncertainty in radiocarbon ages because of the heterogeneous composition of organic layers. In contrast, accelerator mass spectrometry (AMS) radiocarbon dates may be obtained from small samples that have known carbon sources and also allow for a more complete sampling of the section. We present 65 new AMS radiocarbon dates that span nine ground‐rupturing earthquakes at Pallett Creek. Overall, the AMS dates are similar to and reveal no dramatic bias in the conventional dates. For many layers, however, individual charcoal samples were younger than the conventional dates, leading to earthquake ages that are overall slightly younger than previously reported. New earthquake ages are determined by Bayesian refinement of the layer ages based on stratigraphic ordering and sedimentological constraints. The new chronology is more regular than previously published records in large part due to new samples constraining the age of event R. The closed interval from event C to 1857 has a mean recurrence of 135 years ( σ = 83.2 years) and a quasiperiodic coefficient of variation (COV) of 0.61. We show that the new dates and resultant earthquake chronology have a stronger effect on COV than the specific membership of this long series and dating precision improvements from sedimentation rates.

California↗

USGS2011

The Augustine magmatic system as revealed by seismic tomography and relocated earthquake hypocenters from 1994 through 2009

We incorporate 14 years of earthquake data from the Alaska Volcano Observatory with data from a 1975 controlled‐source seismic experiment to obtain the three‐dimensional P and S wave velocity structure and the first high‐precision earthquake locations at Augustine Volcano to be calculated in a fully three‐dimensional velocity model. Velocity tomography shows two main features beneath Augustine: a narrow, high‐velocity column beneath the summit, extending from ∼2 km depth to the surface, and elevated velocities on the south flank. Our relocation results allow a thorough analysis of the spatio‐temoral patterns of seismicity and the relationship to the magmatic and eruptive activity. Background seismicity is centered beneath the summit at an average depth of 0.6 km above sea level. In the weeks leading to the January 2006 eruption of Augustine, seismicity focused on a NW‐SE line along the trend of an inflating dike. A series of drumbeat earthquakes occurred in the early weeks of the eruption, indicating further magma transport through the same dike system. During the six months following the onset of the eruption, the otherwise quiescent region 1 to 5 km below sea level centered beneath the summit became seismically active with two groups of earthquakes, differentiated by frequency content. The deep longer‐period earthquakes occurred during the eruption and are interpreted as resulting from the movement of magma toward the summit, and the post‐eruptive shorter‐period earthquakes may be due to the relaxation of an emptied magma tube. The seismicity subsequently returned to its normal background rates and patterns.

Alaska↗

USGS2011

A model for seasonal changes in GPS positions and seismic wave speeds due to thermoelastic and hydrologic variations

It is known that GPS time series contain a seasonal variation that is not due to tectonic motions, and it has recently been shown that crustal seismic velocities may also vary seasonally. In order to explain these changes, a number of hypotheses have been given, among which thermoelastic and hydrology-induced stresses and strains are leading candidates. Unfortunately, though, since a general framework does not exist for understanding such seasonal variations, it is currently not possible to quickly evaluate the plausibility of these hypotheses. To fill this gap in the literature, I generalize a two-dimensional thermoelastic strain model to provide an analytic solution for the displacements and wave speed changes due to either thermoelastic stresses or hydrologic loading, which consists of poroelastic stresses and purely elastic stresses. The thermoelastic model assumes a periodic surface temperature, and the hydrologic models similarly assume a periodic near-surface water load. Since all three models are two-dimensional and periodic, they are expected to only approximate any realistic scenario; but the models nonetheless provide a quantitative framework for estimating the effects of thermoelastic and hydrologic variations. Quantitative comparison between the models and observations is further complicated by the large uncertainty in some of the relevant parameters. Despite this uncertainty, though, I find that maximum realistic thermoelastic effects are unlikely to explain a large fraction of the observed annual variation in a typical GPS displacement time series or of the observed annual variations in seismic wave speeds in southern California. Hydrologic loading, on the other hand, may be able to explain a larger fraction of both the annual variations in displacements and seismic wave speeds. Neither model is likely to explain all of the seismic wave speed variations inferred from observations. However, more definitive conclusions cannot be made until the model parameters are better constrained. Copyright ?? 2011 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2011

Seismic and geodetic signatures of fault slip at the Slumgullion Landslide Natural Laboratory

We tested the hypothesis that the Slumgullion landslide is a useful natural laboratory for observing fault slip, specifically that slip along its basal surface and side-bounding strike-slip faults occurs with comparable richness of aseismic and seismic modes as along crustal- and plate-scale boundaries. Our study provides new constraints on models governing landslide motion. We monitored landslide deformation with temporary deployments of a 29-element prism array surveyed by a robotic theodolite and an 88-station seismic network that complemented permanent extensometers and environmental instrumentation. Aseismic deformation observations show that large blocks of the landslide move steadily at approximately centimeters per day, possibly punctuated by variations of a few millimeters, while localized transient slip episodes of blocks less than a few tens of meters across occur frequently. We recorded a rich variety of seismic signals, nearly all of which originated outside the monitoring network boundaries or from the side-bounding strike-slip faults. The landslide basal surface beneath our seismic network likely slipped almost completely aseismically. Our results provide independent corroboration of previous inferences that dilatant strengthening along sections of the side-bounding strike-slip faults controls the overall landslide motion, acting as seismically radiating brakes that limit acceleration of the aseismically slipping basal surface. Dilatant strengthening has also been invoked in recent models of transient slip and tremor sources along crustal- and plate-scale faults suggesting that the landslide may indeed be a useful natural laboratory for testing predictions of specific mechanisms that control fault slip at all scales.

Colorado↗

USGS2011

InSAR observations of aseismic slip associated with an earthquake swarm in the Columbia River flood basalts

In 2009 a swarm of small shallow earthquakes occurred within the basalt flows of the Columbia River Basalt Group (CRBG). The swarm occurred within a dense seismic network in the U.S. Department of Energys Hanford Site. Data from the seismic network along with interferometric synthetic aperture radar (InSAR) data from the European Space Agencys (ESA) ENVISAT satellite provide insight into the nature of the swarm. By modeling the InSAR deformation data we constructed a model that consists of a shallow thrust fault and a near horizontal fault. We suggest that the near horizontal lying fault is a bedding-plane fault located between basalt flows. The geodetic moment of the modeled fault system is about eight times the cumulative seismic moment of the swarm. Precise location estimates of the swarm earthquakes indicate that the area of highest slip on the thrust fault, ???70mm of slip less than ???0.5km depth, was not located within the swarm cluster. Most of the slip on the faults appears to have progressed aseismically and we suggest that interbed sediments play a central role in the slip process. Copyright 2011 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2011

Historical perspective on seismic hazard to Hispaniola and the northeast Caribbean region

We evaluate the long-term seismic activity of the North-American/Caribbean plate boundary from 500 years of historical earthquake damage reports. The 2010 Haiti earthquakes and other earthquakes were used to derive regional attenuation relationships between earthquake intensity, magnitude, and distance from the reported damage to the epicenter, for Hispaniola and for Puerto Rico and the Virgin Islands. The attenuation relationship for Hispaniola earthquakes and northern Lesser Antilles earthquakes is similar to that for California earthquakes, indicating a relatively rapid attenuation of damage intensity with distance. Intensities in Puerto Rico and the Virgin Islands decrease less rapidly with distance. We use the intensity-magnitude relationships to systematically search for the location and intensity magnitude MI which best fit all the reported damage for historical earthquakes. Many events occurred in the 20th-century along the plate-boundary segment from central Hispaniola to the NW tip of Puerto Rico, but earlier events from this segment were not identified. The remaining plate boundary to the east to Guadeloupe is probably not associated with M > 8 historical subduction-zone earthquakes. The May 2, 1787 earthquake, previously assigned an M 8–8.25, is probably only M I 6.9 and could be located north, west or SW of Puerto Rico. An M I 6.9 earthquake on July 11, 1785 was probably located north or east of the Virgin Islands. We located M I < 8 historical earthquakes on April 5, 1690, February 8, 1843, and October 8, 1974 in the northern Lesser Antilles within the arc. We speculate that the December 2, 1562 (M I 7.7) and May 7, 1842 (M I 7.6) earthquakes ruptured the Septentrional Fault in northern Hispaniola. If so, the recurrence interval on the central Septentrional Fault is ~300 years, and only 170 years has elapsed since the last event. The recurrence interval of large earthquakes along the Hispaniola subduction segment is likely longer than the historical record. Intra-arc M ≥ 7.0 earthquakes may occur every 75–100 years in the 410-km-long segment between the Virgin Islands and Guadeloupe.

Caribbean Sea↗

USGS2011

Estimates of stress drop and crustal tectonic stress from the 27 February 2010 Maule, Chile, earthquake: Implications for fault strength

The great 27 February 2010 M w 8.8 earthquake off the coast of southern Chile ruptured a ∼600 km length of subduction zone. In this paper, we make two independent estimates of shear stress in the crust in the region of the Chile earthquake. First, we use a coseismic slip model constrained by geodetic observations from interferometric synthetic aperture radar (InSAR) and GPS to derive a spatially variable estimate of the change in static shear stress along the ruptured fault. Second, we use a static force balance model to constrain the crustal shear stress required to simultaneously support observed fore‐arc topography and the stress orientation indicated by the earthquake focal mechanism. This includes the derivation of a semianalytic solution for the stress field exerted by surface and Moho topography loading the crust. We find that the deviatoric stress exerted by topography is minimized in the limit when the crust is considered an incompressible elastic solid, with a Poisson ratio of 0.5, and is independent of Young's modulus. This places a strict lower bound on the critical stress state maintained by the crust supporting plastically deformed accretionary wedge topography. We estimate the coseismic shear stress change from the Maule event ranged from −6 MPa (stress increase) to 17 MPa (stress drop), with a maximum depth‐averaged crustal shear‐stress drop of 4 MPa. We separately estimate that the plate‐driving forces acting in the region, regardless of their exact mechanism, must contribute at least 27 MPa trench‐perpendicular compression and 15 MPa trench‐parallel compression. This corresponds to a depth‐averaged shear stress of at least 7 MPa. The comparable magnitude of these two independent shear stress estimates is consistent with the interpretation that the section of the megathrust fault ruptured in the Maule earthquake is weak, with the seismic cycle relieving much of the total sustained shear stress in the crust.

Maule↗

USGS2012

Slab1.0: A three-dimensional model of global subduction zone geometries

We describe and present a new model of global subduction zone geometries, called Slab1.0. An extension of previous efforts to constrain the two-dimensional non-planar geometry of subduction zones around the focus of large earthquakes, Slab1.0 describes the detailed, non-planar, three-dimensional geometry of approximately 85% of subduction zones worldwide. While the model focuses on the detailed form of each slab from their trenches through the seismogenic zone, where it combines data sets from active source and passive seismology, it also continues to the limits of their seismic extent in the upper-mid mantle, providing a uniform approach to the definition of the entire seismically active slab geometry. Examples are shown for two well-constrained global locations; models for many other regions are available and can be freely downloaded in several formats from our new Slab1.0 website, http://on.doi.gov/d9ARbS . We describe improvements in our two-dimensional geometry constraint inversion, including the use of ‘average’ active source seismic data profiles in the shallow trench regions where data are otherwise lacking, derived from the interpolation between other active source seismic data along-strike in the same subduction zone. We include several analyses of the uncertainty and robustness of our three-dimensional interpolation methods. In addition, we use the filtered, subduction-related earthquake data sets compiled to build Slab1.0 in a reassessment of previous analyses of the deep limit of the thrust interface seismogenic zone for all subduction zones included in our global model thus far, concluding that the width of these seismogenic zones is on average 30% larger than previous studies have suggested.

Journal of Geophysical Research B: Solid Earth↗

USGS2011

Shallow conduit system at Kilauea Volcano, Hawaii, revealed by seismic signals associated with degassing bursts

Eruptive activity at the summit of Kilauea Volcano, Hawaii, beginning in March, 2008 and continuing to the present time is characterized by episodic explosive bursts of gas and ash from a vent within Halemaumau Pit Crater. These bursts are accompanied by seismic signals that are well recorded by a broadband network deployed in the summit caldera. We investigate in detail the dimensions and oscillation modes of the source of a representative burst in the 1−10 s band. An extended source is realized by a set of point sources distributed on a grid surrounding the source centroid, where the centroid position and source geometry are fixed from previous modeling of very-long-period (VLP) data in the 10–50 s band. The source time histories of all point sources are obtained simultaneously through waveform inversion carried out in the frequency domain. Short-scale noisy fluctuations of the source time histories between adjacent sources are suppressed with a smoothing constraint, whose strength is determined through a minimization of the Akaike Bayesian Information Criterion (ABIC). Waveform inversions carried out for homogeneous and heterogeneous velocity structures both image a dominant source component in the form of an east trending dike with dimensions of 2.9 × 2.9 km. The dike extends ∼2 km west and ∼0.9 km east of the VLP centroid and spans the depth range 0.2–3.1 km. The source model for a homogeneous velocity structure suggests the dike is hinged at the source centroid where it bends from a strike E 27°N with northern dip of 85° west of the centroid, to a strike E 7°N with northern dip of 80° east of the centroid. The oscillating behavior of the dike is dominated by simple harmonic modes with frequencies ∼0.2 Hz and ∼0.5 Hz, representing the fundamental mode ν 11 and first degenerate mode ν 12 = ν 21 of the dike. Although not strongly supported by data in the 1–10 s band, a north striking dike segment is required for enhanced compatibility with the model elaborated in the 10–50 s band. This dike provides connectivity between the east trending dike and the new vent within Halemaumau Pit Crater. Waveform inversions with a dual-dike model suggest dimensions of 0.7 × 0.7 km to 2.6 × 2.6 km for this segment. Further elaboration of the complex dike system under Halemaumau does not appear to be feasible with presently available data.

Hawai'i↗

USGS2012

Geometry and subsidence history of the Dead Sea basin: A case for fluid-induced mid-crustal shear zone?

Pull‐apart basins are narrow zones of crustal extension bounded by strike‐slip faults that can serve as analogs to the early stages of crustal rifting. We use seismic tomography, 2‐D ray tracing, gravity modeling, and subsidence analysis to study crustal extension of the Dead Sea basin (DSB), a large and long‐lived pull‐apart basin along the Dead Sea transform (DST). The basin gradually shallows southward for 50 km from the only significant transverse normal fault. Stratigraphic relationships there indicate basin elongation with time. The basin is deepest (8–8.5 km) and widest (∼15 km) under the Lisan about 40 km north of the transverse fault. Farther north, basin depth is ambiguous, but is 3 km deep immediately north of the lake. The underlying pre‐basin sedimentary layer thickens gradually from 2 to 3 km under the southern edge of the DSB to 3–4 km under the northern end of the lake and 5–6 km farther north. Crystalline basement is ∼11 km deep under the deepest part of the basin. The upper crust under the basin has lower P wave velocity than in the surrounding regions, which is interpreted to reflect elevated pore fluids there. Within data resolution, the lower crust below ∼18 km and the Moho are not affected by basin development. The subsidence rate was several hundreds of m/m.y. since the development of the DST ∼17 Ma, similar to other basins along the DST, but subsidence rate has accelerated by an order of magnitude during the Pleistocene, which allowed the accumulation of 4 km of sediment. We propose that the rapid subsidence and perhaps elongation of the DSB are due to the development of inter‐connected mid‐crustal ductile shear zones caused by alteration of feldspar to muscovite in the presence of pore fluids. This alteration resulted in a significant strength decrease and viscous creep. We propose a similar cause to the enigmatic rapid subsidence of the North Sea at the onset the North Atlantic mantle plume. Thus, we propose that aqueous fluid flux into a slowly extending continental crust can cause rapid basin subsidence that may be erroneously interpreted as an increased rate of tectonic activity.

Dead Sea shoreline↗

USGS2012

Fixed recurrence and slip models better predict earthquake behavior than the time- and slip-predictable models 1: repeating earthquakes

The behavior of individual events in repeating earthquake sequences in California, Taiwan and Japan is better predicted by a model with fixed inter-event time or fixed slip than it is by the time- and slip-predictable models for earthquake occurrence. Given that repeating earthquakes are highly regular in both inter-event time and seismic moment, the time- and slip-predictable models seem ideally suited to explain their behavior. Taken together with evidence from the companion manuscript that shows similar results for laboratory experiments we conclude that the short-term predictions of the time- and slip-predictable models should be rejected in favor of earthquake models that assume either fixed slip or fixed recurrence interval. This implies that the elastic rebound model underlying the time- and slip-predictable models offers no additional value in describing earthquake behavior in an event-to-event sense, but its value in a long-term sense cannot be determined. These models likely fail because they rely on assumptions that oversimplify the earthquake cycle. We note that the time and slip of these events is predicted quite well by fixed slip and fixed recurrence models, so in some sense they are time- and slip-predictable. While fixed recurrence and slip models better predict repeating earthquake behavior than the time- and slip-predictable models, we observe a correlation between slip and the preceding recurrence time for many repeating earthquake sequences in Parkfield, California. This correlation is not found in other regions, and the sequences with the correlative slip-predictable behavior are not distinguishable from nearby earthquake sequences that do not exhibit this behavior.

Journal of Geophysical Research B: Solid Earth↗

USGS2012

Holocene faulting in the Bellingham forearc basin: Upper-plate deformation at the northern end of the Cascadia subduction zone

The northern Cascadia forearc takes up most of the strain transmitted northward via the Oregon Coast block from the northward-migrating Sierra Nevada block. The north-south contractional strain in the forearc manifests in upper-plate faults active during the Holocene, the northern-most components of which are faults within the Bellingham Basin. The Bellingham Basin is the northern of four basins of the actively deforming northern Cascadia forearc. A set of Holocene faults, Drayton Harbor, Birch Bay, and Sandy Point faults, occur within the Bellingham Basin and can be traced from onshore to offshore using a combination of aeromagnetic lineaments, paleoseismic investigations and scarps identified using LiDAR imagery. With the recognition of such Holocene faults, the northernmost margin of the actively deforming Cascadia forearc extends 60 km north of the previously recognized limit of Holocene forearc deformation. Although to date no Holocene faults are recognized at the northern boundary of the Bellingham Basin, which is 15 km north of the international border, there is no compelling tectonic reason to expect that Holocene faults are limited to south of the international border.

Washington↗

USGS2012

Why the 2002 Denali fault rupture propagated onto the Totschunda fault: implications for fault branching and seismic hazards

The propagation of the rupture of the M w 7.9 Denali fault earthquake from the central Denali fault onto the Totschunda fault has provided a basis for dynamic models of fault branching in which the angle of the regional or local prestress relative to the orientation of the main fault and branch plays a principal role in determining which fault branch is taken. GeoEarthScope LiDAR and paleoseismic data allow us to map the structure of the Denali-Totschunda fault intersection and evaluate controls of fault branching from a geological perspective. LiDAR data reveal the Denali-Totschunda fault intersection is structurally simple with the two faults directly connected. At the branch point, 227.2 km east of the 2002 epicenter, the 2002 rupture diverges southeast to become the Totschunda fault. We use paleoseismic data to propose that differences in the accumulated strain on each fault segment, which express differences in the elapsed time since the most recent event, was one important control of the branching direction. We suggest that data on event history, slip rate, paleo offsets, fault geometry and structure, and connectivity, especially on high slip rate-short recurrence interval faults, can be used to assess the likelihood of branching and its direction. Analysis of the Denali-Totschunda fault intersection has implications for evaluating the potential for a rupture to propagate across other types of fault intersections and for characterizing sources of future large earthquakes.

Alaska↗

USGS2012

Juan de Fuca slab geometry and its relation to Wadati-Benioff zone seismicity

A new model of the subducted Juan de Fuca plate beneath western North America allows first-order correlations between the occurrence of Wadati-Benioff zone earthquakes and slab geometry, temperature, and hydration state. The geo-referenced 3D model, constructed from weighted control points, integrates depth information from earthquake locations and regional seismic velocity studies. We use the model to separate earthquakes that occur in the Cascadia forearc from those that occur within the underlying Juan de Fuca plate and thereby reveal previously obscured details regarding the spatial distribution of earthquakes. Seismicity within the slab is most prevalent where the slab is warped beneath northwestern California and western Washington suggesting that slab flexure, in addition to expected metamorphic dehydration processes, promotes earthquake occurrence within the subducted oceanic plate. Earthquake patterns beneath western Vancouver Island are consistent with slab dehydration processes. Conversely, the lack of slab earthquakes beneath western Oregon is consistent with an anhydrous slab. Double-differenced relocated seismicity resolves a double seismic zone within the slab beneath northwestern California that strongly constrains the location of the plate interface and delineates a cluster of seismicity 10 km above the surface that includes the 1992 M7.1 Mendocino earthquake. We infer that this earthquake ruptured a surface within the Cascadia accretionary margin above the Juan de Fuca plate. We further speculate that this earthquake is associated with a detached fragment of former Farallon plate. Other subsurface tectonic elements within the forearc may have the potential to generate similar damaging earthquakes.

British Columbia;California;Oregon;Washington↗

USGS2012

Laboratory triggering of stick-slip events by oscillatory loading in the presence of pore fluid with implications for physics of tectonic tremor

The physical mechanism by which the low-frequency earthquakes (LFEs) that make up portions of tectonic (also called non-volcanic) tremor are created is poorly understood. In many areas of the world, tectonic tremor and LFEs appear to be strongly tidally modulated, whereas ordinary earthquakes are not. Anomalous seismic wave speeds, interpreted as high pore fluid pressure, have been observed in regions that generate tremor. Here we build upon previous laboratory studies that investigated the response of stick-slip on artificial faults to oscillatory, tide-like loading. These previous experiments were carried out using room-dry samples of Westerly granite, at one effective stress. Here we augment these results with new experiments on Westerly granite, with the addition of varying effective stress using pore fluid at two pressures. We find that raising pore pressure, thereby lowering effective stress can significantly increase the degree of correlation of stick-slip to oscillatory loading. We also find other pore fluid effects that become important at higher frequencies, when the period of oscillation is comparable to the diffusion time of pore fluid into the fault. These results help constrain the conditions at depth that give rise to tidally modulated LFEs, providing confirmation of the effective pressure law for triggering and insights into why tremor is tidally modulated while earthquakes are at best only weakly modulated.

Journal of Geophysical Research B: Solid Earth↗

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↗

USGS2012

Externally triggered renewed bubble nucleation in basaltic magma: the 12 October 2008 eruption at Halema‘uma‘u Overlook vent, Kīlauea, Hawai‘i, USA

From October 2008 until present, dozens of small impulsive explosive eruptions occurred from the Overlook vent on the southeast side of Halema‘uma‘u Crater, at Kīlauea volcano, USA. These eruptions were triggered by rockfalls from the walls of the volcanic vent and conduit onto the top of the lava column. Here we use microtextural observations and data from clasts erupted during the well-characterized 12 October 2008 explosive eruption at Halema‘uma‘u to extend existing models of eruption triggering. We present a potential mechanism for this eruption by combining microtextural observations with existing geophysical and visual data sets. We measure the size and number density of bubbles preserved in juvenile ejecta using 2D images and X-ray microtomography. Our data suggest that accumulations of large bubbles with diameters of >50 μ m to at least millimeters existed at shallow levels within the conduit prior to the 12 October 2008 explosion. Furthermore, a high number density of small bubbles <50 μ m is measured in the clasts, implying very rapid nucleation of bubbles. Visual observations, combined with preexisting geophysical data, suggest that the impact of rockfalls onto the magma free surface induces pressure changes over short timescales that (1) nucleated new additional bubbles in the shallow conduit leading to high number densities of small bubbles and (2) expanded the preexisting bubbles driving upward acceleration. The trigger of eruption and bubble nucleation is thus external to the degassing system.

Hawaii↗