USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research B: Solid Earth”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8Linked to original sources
USGS2018

Using stereo satellite imagery to account for ablation, entrainment, and compaction in volume calculations for rock avalanches on Glaciers: Application to the 2016 Lamplugh Rock Avalanche in Glacier Bay National Park, Alaska

The use of preevent and postevent digital elevation models (DEMs) to estimate the volume of rock avalanches on glaciers is complicated by ablation of ice before and after the rock avalanche, scour of material during rock avalanche emplacement, and postevent ablation and compaction of the rock avalanche deposit. We present a model to account for these processes in volume estimates of rock avalanches on glaciers. We applied our model by calculating the volume of the 28 June 2016 Lamplugh rock avalanche in Glacier Bay National Park, Alaska. We derived preevent and postevent 2‐m resolution DEMs from WorldView satellite stereo imagery. Using data from DEM differencing, we reconstructed the rock avalanche and adjacent surfaces at the time of occurrence by accounting for elevation changes due to ablation and scour of the ice surface, and postevent deposit changes. We accounted for uncertainties in our DEMs through precise coregistration and an assessment of relative elevation accuracy in bedrock control areas. The rock avalanche initially displaced 51.7 ± 1.5 Mm3 of intact rock and then scoured and entrained 13.2 ± 2.2 Mm3 of snow and ice during emplacement. We calculated the total deposit volume to be 69.9 ± 7.9 Mm3. Volume estimates that did not account for topographic changes due to ablation, scour, and compaction underestimated the deposit volume by 31.0–46.8 Mm3. Our model provides an improved framework for estimating uncertainties affecting rock avalanche volume measurements in glacial environments. These improvements can contribute to advances in the understanding of rock avalanche hazards and dynamics.

Alaska↗
USGS2013

Missing great earthquakes

The occurrence of three earthquakes with moment magnitude (Mw) greater than 8.8 and six earthquakes larger than Mw 8.5, since 2004, has raised interest in the long-term global rate of great earthquakes. Past studies have focused on the analysis of earthquakes since 1900, which roughly marks the start of the instrumental era in seismology. Before this time, the catalog is less complete and magnitude estimates are more uncertain. Yet substantial information is available for earthquakes before 1900, and the catalog of historical events is being used increasingly to improve hazard assessment. Here I consider the catalog of historical earthquakes and show that approximately half of all Mw ≥ 8.5 earthquakes are likely missing or underestimated in the 19th century. I further present a reconsideration of the felt effects of the 8 February 1843, Lesser Antilles earthquake, including a first thorough assessment of felt reports from the United States, and show it is an example of a known historical earthquake that was significantly larger than initially estimated. The results suggest that incorporation of best available catalogs of historical earthquakes will likely lead to a significant underestimation of seismic hazard and/or the maximum possible magnitude in many regions, including parts of the Caribbean.

Journal of Geophysical Research B: Solid Earth↗
USGS2013

The utility of atmospheric analyses for the mitigation of artifacts in InSAR

The numerical weather models (NWMs) developed by the meteorological community are able to provide accurate analyses of the current state of the atmosphere in addition to the predictions of the future state. To date, most attempts to apply the NWMs to estimate the refractivity of the atmosphere at the time of satellite synthetic aperture radar (SAR) data acquisitions have relied on predictive models. We test the hypothesis that performing a final assimilative routine, ingesting all available meteorological observations for the times of SAR acquisitions, and generating customized analyses of the atmosphere at those times will better mitigate atmospheric artifacts in differential interferograms. We find that, for our study area around Mount St. Helens (Amboy, Washington, USA), this approach is unable to model the refractive changes and provides no mean benefit for interferogram analysis. The performance is improved slightly by ingesting atmospheric delay estimates derived from the limited local GPS network; however, the addition of water vapor products from the GOES satellites reduces the quality of the corrections. We interpret our results to indicate that, even with this advanced approach, NWMs are not a reliable mitigation technique for regions such as Mount St. Helens with highly variable moisture fields and complex topography and atmospheric dynamics. It is possible, however, that the addition of more spatially dense meteorological data to constrain the analyses might significantly improve the performance of weather modeling of atmospheric artifacts in satellite radar interferograms.

Journal of Geophysical Research B: Solid Earth↗
USGS2013

Heterogeneous rupture in the great Cascadia earthquake of 1700 inferred from coastal subsidence estimates

Past earthquake rupture models used to explain paleoseismic estimates of coastal subsidence during the great A.D. 1700 Cascadia earthquake have assumed a uniform slip distribution along the megathrust. Here we infer heterogeneous slip for the Cascadia margin in A.D. 1700 that is analogous to slip distributions during instrumentally recorded great subduction earthquakes worldwide. The assumption of uniform distribution in previous rupture models was due partly to the large uncertainties of then available paleoseismic data used to constrain the models. In this work, we use more precise estimates of subsidence in 1700 from detailed tidal microfossil studies. We develop a 3-D elastic dislocation model that allows the slip to vary both along strike and in the dip direction. Despite uncertainties in the updip and downdip slip extensions, the more precise subsidence estimates are best explained by a model with along-strike slip heterogeneity, with multiple patches of high-moment release separated by areas of low-moment release. For example, in A.D. 1700, there was very little slip near Alsea Bay, Oregon (~44.4°N), an area that coincides with a segment boundary previously suggested on the basis of gravity anomalies. A probable subducting seamount in this area may be responsible for impeding rupture during great earthquakes. Our results highlight the need for more precise, high-quality estimates of subsidence or uplift during prehistoric earthquakes from the coasts of southern British Columbia, northern Washington (north of 47°N), southernmost Oregon, and northern California (south of 43°N), where slip distributions of prehistoric earthquakes are poorly constrained.

British Columbia;Washington;Oregon;California↗
USGS2013

Evolution of dike opening during the March 2011 Kamoamoa fissure eruption, Kīlauea Volcano, Hawai`i

The 5–9 March 2011 Kamoamoa fissure eruption along the east rift zone of Kīlauea Volcano, Hawai`i, followed months of pronounced inflation at Kīlauea summit. We examine dike opening during and after the eruption using a comprehensive interferometric synthetic aperture radar (InSAR) data set in combination with continuous GPS data. We solve for distributed dike displacements using a whole Kīlauea model with dilating rift zones and possibly a deep décollement. Modeled surface dike opening increased from nearly 1.5 m to over 2.8 m from the first day to the end of the eruption, in agreement with field observations of surface fracturing. Surface dike opening ceased following the eruption, but subsurface opening in the dike continued into May 2011. Dike volumes increased from 15, to 16, to 21 million cubic meters (MCM) after the first day, eruption end, and 2 months following, respectively. Dike shape is distinctive, with a main limb plunging from the surface to 2–3 km depth in the up‐rift direction toward Kīlauea's summit, and a lesser projection extending in the down‐rift direction toward Pu`u `Ō`ō at 2 km depth. Volume losses beneath Kīlauea summit (1.7 MCM) and Pu`u `Ō`ō (5.6 MCM) crater, relative to dike plus erupted volume (18.3 MCM), yield a dike to source volume ratio of 2.5 that is in the range expected for compressible magma without requiring additional sources. Inflation of Kīlauea's summit in the months before the March 2011 eruption suggests that the Kamoamoa eruption resulted from overpressure of the volcano's magmatic system.

Hawai'i↗
USGS2013

InSAR Evidence for an active shallow thrust fault beneath the city of Spokane Washington, USA

In 2001, a nearly five month long sequence of shallow, mostly small magnitude earthquakes occurred beneath the city of Spokane, a city with a population of about 200,000, in the state of Washington. During most of the sequence, the earthquakes were not well located because seismic instrumentation was sparse. Despite poor-quality locations, the earthquake hypocenters were likely very shallow, because residents near the city center both heard and felt many of the earthquakes. The combination of poor earthquake locations and a lack of known surface faults with recent movement make assessing the seismic hazards related to the earthquake swarm difficult. However, the potential for destruction from a shallow moderate-sized earthquake is high, for example Christchurch New Zealand in 2011, so assessing the hazard potential of a seismic structure involved in the Spokane earthquake sequence is important. Using interferometric synthetic aperture radar (InSAR) data from the European Space Agency ERS2 and ENVISAT satellites and the Canadian Space Agency RADARSAT-1, satellite we are able to show that slip on a shallow previously unknown thrust fault, which we name the Spokane Fault, is the source of the earthquake sequence. The part of the Spokane Fault that slipped during the 2001 earthquake sequence underlies the north part of the city, and slip on the fault was concentrated between ~0.3 and 2 km depth. Projecting the buried fault plane to the surface gives a possible surface trace for the Spokane Fault that strikes northeast from the city center into north Spokane.

Washington↗
USGS2013

Chemical controls on fault behavior: weakening of serpentinite sheared against quartz-bearing rocks and its significance for fault creep in the San Andreas system

The serpentinized ultramafic rocks found in many plate-tectonic settings commonly are juxtaposed against crustal rocks along faults, and the chemical contrast between the rock types potentially could influence the mechanical behavior of such faults. To investigate this possibility, we conducted triaxial experiments under hydrothermal conditions (200-350°C), shearing serpentinite gouge between forcing blocks of granite or quartzite. In an ultramafic chemical environment, the coefficient of friction, µ, of lizardite and antigorite serpentinite is 0.5-0.6, and µ increases with increasing temperature over the tested range. However, when either lizardite or antigorite serpentinite is sheared against granite or quartzite, strength is reduced to µ ~ 0.3, with the greatest strength reductions at the highest temperatures (temperature weakening) and slowest shearing rates (velocity strengthening). The weakening is attributed to a solution-transfer process that is promoted by the enhanced solubility of serpentine in pore fluids whose chemistry has been modified by interaction with the quartzose wall rocks. The operation of this process will promote aseismic slip (creep) along serpentinite-bearing crustal faults at otherwise seismogenic depths. During short-term experiments serpentine minerals reprecipitate in low-stress areas, whereas in longer experiments new Mg-rich phyllosilicates crystallize in response to metasomatic exchanges across the serpentinite-crustal rock contact. Long-term shear of serpentinite against crustal rocks will cause the metasomatic mineral assemblages, which may include extremely weak minerals such as saponite or talc, to play an increasingly important role in the mechanical behavior of the fault. Our results may explain the distribution of creep on faults in the San Andreas system.

California↗
USGS2013

Bayesian inversion of data from effusive volcanic eruptions using physics-based models: Application to Mount St. Helens 2004--2008

Physics-based models of volcanic eruptions can directly link magmatic processes with diverse, time-varying geophysical observations, and when used in an inverse procedure make it possible to bring all available information to bear on estimating properties of the volcanic system. We develop a technique for inverting geodetic, extrusive flux, and other types of data using a physics-based model of an effusive silicic volcanic eruption to estimate the geometry, pressure, depth, and volatile content of a magma chamber, and properties of the conduit linking the chamber to the surface. A Bayesian inverse formulation makes it possible to easily incorporate independent information into the inversion, such as petrologic estimates of melt water content, and yields probabilistic estimates for model parameters and other properties of the volcano. Probability distributions are sampled using a Markov-Chain Monte Carlo algorithm. We apply the technique using GPS and extrusion data from the 2004–2008 eruption of Mount St. Helens. In contrast to more traditional inversions such as those involving geodetic data alone in combination with kinematic forward models, this technique is able to provide constraint on properties of the magma, including its volatile content, and on the absolute volume and pressure of the magma chamber. Results suggest a large chamber of >40 km 3 with a centroid depth of 11–18 km and a dissolved water content at the top of the chamber of 2.6–4.9 wt%.

Washington↗
USGS2013

Comment on “Historical perspective on seismic hazard to Hispaniola and the northeast Caribbean region” by U. ten Brink et al.

The analysis of historical earthquakes in the northeastern Caribbean by ten Brink et al. [2011, hereafter TB11] addresses the occurrence of large and destructive historical earthquakes associated with the North American-Caribbean plate boundary. One conclusion presented in TB11 is that the recurrence interval for large earthquakes on the left-lateral, strike-slip Septentrional Fault (SF) (Figure 1a) is approximately 300 years. Their Figure 7 shows rupture of the SF across the entire island of Hispaniola in CE 1200, 1542, and 1842. Our comment challenges this model for SF earthquake recurrence because it is inconsistent with our published paleoseismic data that show no large historical earthquake is associated with surface rupture along the SF east of Santiago (Figure 1a)[Prentice et al., 1993; Mann et al., 1998; Prentice et al., 2003].

Hispaniola↗
USGS2013

The role of farfield tectonic stress in oceanic intraplate deformation, Gulf of Alaska

An integration of geophysical data from the Pacific Plate reveals plate bending anomalies, massive intraplate shearing and deformation, and a lack of oceanic crust magnetic lineaments in different regions across the Gulf of Alaska. We argue that farfield stress from the Yakutat Terrane collision with North America is the major driver for these unusual features. Similar plate motion vectors indicate that the Pacific plate and Yakutat Terrane are largely coupled along their boundary, the Transition Fault, with minimal translation. Our study shows that the Pacific Plate subduction angle shallows toward the Yakutat Terrane and supports the theory that the Pacific Plate and Yakutat Terranemaintain coupling along the subducted region of the Transition Fault. We argue that the outboard transfer of collisional stress to the Pacific Plate could have resulted in significant strain in the NE corner of the Pacific Plate, which created pathways for igneous sill formation just above the Pacific Plate crust in the Surveyor Fan. A shift in Pacific Plate motion during the late Miocene altered the Yakutat collision with North America, changing the stress transfer regime and potentially terminating associated strain in the NE corner of the Pacific Plate. The collision further intensified as the thickest portion of the Yakutat Terrane began to subduct during the Pleistocene, possibly providing the impetus for the creation of the Gulf of Alaska Shear Zone, a>200 km zone of intraplate strike-slip faults that extend from the Transition Fault out into the Pacific Plate. This study highlights the importance of farfield stress from complex tectonic regimes in consideration of large-scale oceanic intraplate deformation.

Gulf of Alaska↗
USGS2013

Annual modulation of non-volcanic tremor in northern Cascadia

Two catalogs of episodic tremor events in northern Cascadia, one from 2006 to 2012 and the other from 1997 to 2011, reveal two systematic patterns of tremor occurrence in southern Vancouver Island: (1) most individual events tend to occur in the third quarter of the year; (2) the number of events in prolonged episodes (i.e., episodic tremor and slip events), which generally propagate to Vancouver Island from elsewhere along the Cascadia subduction zone, is inversely correlated with the amount of precipitation that occurred in the preceding 2 months. We rationalize these patterns as the product of hydrologic loading of the crust of southern Vancouver Island and the surrounding continental region, superimposed with annual variations from oceanic tidal loading. Loading of the Vancouver Island crust in the winter (when the land surface receives ample precipitation) and unloading in the summer tends to inhibit and enhance downdip shear stress, respectively. Quantitatively, for an annually variable surface load, the predicted stress perturbation depends on mantle viscoelastic rheology. A mechanical model of downdip shear stress on the transition zone beneath Vancouver Island—driven predominantly by the annual hydrologic cycle—is consistent with the 1997–2012 tremor observations, with peak-to-peak downdip shear stress of about 0.4 kPa. This seasonal dependence of tremor occurrence appears to be restricted to southern Vancouver Island because of its unique situation as an elongated narrow-width land mass surrounded by ocean, which permits seasonal perturbations in shear stress at depth.

Journal of Geophysical Research B: Solid Earth↗
USGS2013

The role of viscous magma mush spreading in volcanic flank motion at Kīlauea Volcano, Hawai‘i

Multiple mechanisms have been suggested to explain seaward motion of the south flank of Kīlauea Volcano, Hawai‘i. The consistency of flank motion during both waxing and waning magmatic activity at Kīlauea suggests that a continuously acting force, like gravity body force, plays a substantial role. Using finite element models, we test whether gravity is the principal driver of long-term motion of Kīlauea's flank. We compare our model results to geodetic data from Global Positioning System and interferometric synthetic aperture radar during a time period with few magmatic and tectonic events (2000-2003), when deformation of Kīlauea was dominated by summit subsidence and seaward motion of the south flank. We find that gravity-only models can reproduce the horizontal surface velocities if we incorporate a regional décollement fault and a deep, low-viscosity magma mush zone. To obtain quasi steady state horizontal surface velocities that explain the long-term seaward motion of the flank, we find that an additional weak zone is needed, which is an extensional rift zone above the magma mush. The spreading rate in our model is mainly controlled by the magma mush viscosity, while its density plays a less significant role. We find that a viscosity of 2.5 × 1017–2.5 × 1019 Pa s for the magma mush provides an acceptable fit to the observed horizontal surface deformation. Using high magma mush viscosities, such as 2.5 × 1019 Pa s, the deformation rates remain more steady state over longer time scales. These models explain a significant amount of the observed subsidence at Kīlauea's summit. Some of the remaining subsidence is probably a result of magma withdrawal from subsurface reservoirs.

Hawai'i↗
USGS2013

A domain decomposition approach to implementing fault slip in finite-element models of quasi-static and dynamic crustal deformation

We employ a domain decomposition approach with Lagrange multipliers to implement fault slip in a finite-element code, PyLith, for use in both quasi-static and dynamic crustal deformation applications. This integrated approach to solving both quasi-static and dynamic simulations leverages common finite-element data structures and implementations of various boundary conditions, discretization schemes, and bulk and fault rheologies. We have developed a custom preconditioner for the Lagrange multiplier portion of the system of equations that provides excellent scalability with problem size compared to conventional additive Schwarz methods. We demonstrate application of this approach using benchmarks for both quasi-static viscoelastic deformation and dynamic spontaneous rupture propagation that verify the numerical implementation in PyLith.

Journal of Geophysical Research B: Solid Earth↗
USGS2013

Seismic evidence for a slab tear at the Puerto Rico Trench

The fore-arc region of the northeast Caribbean plate north of Puerto Rico and the Virgin Islands has been the site of numerous seismic swarms since at least 1976. A 6 month deployment of five ocean bottom seismographs recorded two such tightly clustered swarms, along with additional events. Joint analyses of the ocean bottom seismographs and land-based seismic data reveal that the swarms are located at depths of 50–150 km. Focal mechanism solutions, found by jointly fitting P wave first-motion polarities and S/P amplitude ratios, indicate that the broadly distributed events outside the swarm generally have strike- and dip-slip mechanisms at depths of 50–100 km, while events at depths of 100–150 km have oblique mechanisms. A stress inversion reveals two distinct stress regimes: The slab segment east of 65°W longitude is dominated by trench-normal tensile stresses at shallower depths (50–100 km) and by trench-parallel tensile stresses at deeper depths (100–150 km), whereas the slab segment west of 65°W longitude has tensile stresses that are consistently trench normal throughout the depth range at which events were observed (50–100 km). The simple stress pattern in the western segment implies relatively straightforward subduction of an unimpeded slab, while the stress pattern observed in the eastern segment, shallow trench-normal tension and deeper trench-normal compression, is consistent with flexure of the slab due to rollback. These results support the hypothesis that the subducting North American plate is tearing at or near these swarms. The 35 year record of seismic swarms at this location and the recent increase in seismicity suggest that the tear is still propagating.

Puerto Rico↗
USGS2013

Foreshocks during the nucleation of stick-slip instability

We report on laboratory experiments which investigate interactions between aseismic slip, stress changes, and seismicity on a critically stressed fault during the nucleation of stick-slip instability. We monitor quasi-static and dynamic changes in local shear stress and fault slip with arrays of gages deployed along a simulated strike-slip fault (2 m long and 0.4 m deep) in a saw cut sample of Sierra White granite. With 14 piezoelectric sensors, we simultaneously monitor seismic signals produced during the nucleation phase and subsequent dynamic rupture. We observe localized aseismic fault slip in an approximately meter-sized zone in the center of the fault, while the ends of the fault remain locked. Clusters of high-frequency foreshocks (M w ~ −6.5 to −5.0) can occur in this slowly slipping zone 5–50 ms prior to the initiation of dynamic rupture; their occurrence appears to be dependent on the rate at which local shear stress is applied to the fault. The meter-sized nucleation zone is generally consistent with theoretical estimates, but source radii of the foreshocks (2 to 70 mm) are 1 to 2 orders of magnitude smaller than the theoretical minimum length scale over which earthquake nucleation can occur. We propose that frictional stability and the transition between seismic and aseismic slip are modulated by local stressing rate and that fault sections, which would typically slip aseismically, may radiate seismic waves if they are rapidly stressed. Fault behavior of this type may provide physical insight into the mechanics of foreshocks, tremor, repeating earthquake sequences, and a minimum earthquake source dimension.

Journal of Geophysical Research B: Solid Earth↗
USGS2013

Eruptions at Lone Star Geyser, Yellowstone National Park, USA, part 1: energetics and eruption dynamics

Geysers provide a natural laboratory to study multiphase eruptive processes. We present results from a four&ndash;day experiment at Lone Star Geyser in Yellowstone National Park, USA. We simultaneously measured water discharge, acoustic emissions, infraredintensity, and visible and infrared video to quantify the energetics and dynamics of eruptions, occurring approximately every three hours. We define four phases in the eruption cycle: 1) a 28&thinsp;&plusmn;&thinsp;3 minute phase with liquid and steam fountaining, with maximum jet velocities of 16&ndash;28&thinsp;m s &minus;&thinsp;1 , steam mass fraction of less than &sim;&thinsp;0.01. Intermittently choked flow and flow oscillations with periods increasing from 20 to 40&thinsp;s are coincident with a decrease in jet velocity and an increase of steam fraction; 2) a 26&thinsp;&plusmn;&thinsp;8 minute post&ndash;eruption relaxation phase with no discharge from the vent, infrared (IR) and acoustic power oscillations gliding between 30 and 40&thinsp;s; 3) a 59&thinsp;&plusmn;&thinsp;13 minute recharge period during which the geyser is quiescent and progressively refills, and 4) a 69&thinsp;&plusmn;&thinsp;14 minute pre&ndash;play period characterized by a series of 5&ndash;10&thinsp;minute&ndash;long pulses of steam, small volumes of liquid water discharge and 50&ndash;70&thinsp;s flow oscillations. The erupted waters ascend froma 160&thinsp;&minus;&thinsp;170&deg; C reservoir and the volume discharged during the entire eruptive cycle is 20.8&thinsp;&plusmn;&thinsp;4.1 m 3 . Assuming isentropic expansion, we calculate a heat output from the geyser of 1.4&ndash;1.5&thinsp;MW, which is <&thinsp;0.1% of the total heat output from Yellowstone Caldera.

Wyoming↗
USGS2013

Seismicity around Parkfield correlates with static shear stress changes following the 2003 M w 6.5 San Simeon earthquake

Earthquakes trigger other earthquakes, but the physical mechanism of the triggering is currently debated. Most studies of earthquake triggering rely on earthquakes listed in catalogs, which are known to be incomplete around the origin times of large earthquakes and therefore missing potentially triggered events. Here we apply a waveform matched-filter technique to systematically detect earthquakes along the Parkfield section of the San Andreas Fault from 46 days before to 31 days after the nearby 2003 M w 6.5 San Simeon earthquake. After removing all possible false detections, we identify ~8 times more earthquakes than in the Northern California Seismic Network catalog. The newly identified events along the creeping section of the San Andreas Fault show a statistically significant decrease following the San Simeon main shock, which correlates well with the negative static stress changes (i.e., stress shadow) cast by the main shock. In comparison, the seismicity rate around Parkfield increased moderately where the static stress changes are positive. The seismicity rate changes correlate well with the static shear stress changes induced by the San Simeon main shock, suggesting a low friction in the seismogenic zone along the Parkfield section of the San Andreas Fault.

California↗