SEARCH · USGS Science

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 289 records · Page 16Linked to original sources

USGS2002

Crustal structure of central Lake Baikal: Insights into intracontinental rifting

The Cenozoic rift system of Baikal, located in the interior of the largest continental mass on Earth, is thought to represent a potential analog of the early stage of breakup of supercontinents. We present a detailed P wave velocity structure of the crust and sediments beneath the Central Basin, the deepest basin in the Baikal rift system. The structure is characterized by a Moho depth of 39–42.5 km; an 8-km-thick, laterally continuous high-velocity (7.05–7.4 km/s) lower crust, normal upper mantle velocity (8 km/s), a sedimentary section reaching maximum depths of 9 km, and a gradual increase of sediment velocity with depth. We interpret the high-velocity lower crust to be part of the Siberian Platform that was not thinned or altered significantly during rifting. In comparison to published results from the Siberian Platform, Moho under the basin is elevated by <3 km. On the basis of these results we propose that the basin was formed by upper crustal extension, possibly reactivating structures in an ancient fold-and-thrust belt. The extent and location of upper mantle extension are not revealed by our data, and it may be offset from the rift. We believe that the Baikal rift structure is similar in many respects to the Mesozoic Atlantic rift system, the precursor to the formation of the North Atlantic Ocean. We also propose that the Central Baikal rift evolved by episodic fault propagation and basin enlargement, rather than by two-stage rift evolution as is commonly assumed.

Lake Baikal↗

USGS2016

Volcano deformation source parameters estimated from InSAR: Sensitivities to uncertainties in seismic tomography

The eruption cycle of a volcano is controlled in part by the upward migration of magma. The characteristics of the magma flux produce a deformation signature at the Earth's surface. Inverse analyses use geodetic data to estimate strategic controlling parameters that describe the position and pressurization of a magma chamber at depth. The specific distribution of material properties controls how observed surface deformation translates to source parameter estimates. Seismic tomography models describe the spatial distributions of material properties that are necessary for accurate models of volcano deformation. This study investigates how uncertainties in seismic tomography models propagate into variations in the estimates of volcano deformation source parameters inverted from geodetic data. We conduct finite element model-based nonlinear inverse analyses of interferometric synthetic aperture radar (InSAR) data for Okmok volcano, Alaska, as an example. We then analyze the estimated parameters and their uncertainties to characterize the magma chamber. Analyses are performed separately for models simulating a pressurized chamber embedded in a homogeneous domain as well as for a domain having a heterogeneous distribution of material properties according to seismic tomography. The estimated depth of the source is sensitive to the distribution of material properties. The estimated depths for the homogeneous and heterogeneous domains are 2666 ± 42 and 3527 ± 56 m below mean sea level, respectively (99% confidence). A Monte Carlo analysis indicates that uncertainties of the seismic tomography cannot account for this discrepancy at the 99% confidence level. Accounting for the spatial distribution of elastic properties according to seismic tomography significantly improves the fit of the deformation model predictions and significantly influences estimates for parameters that describe the location of a pressurized magma chamber.

Journal of Geophysical Research B: Solid Earth↗

USGS2021

Exploring GPS observations of postseismic deformation following the 2012 MW7.8 Haida Gwaii and 2013 MW7.5 Craig, Alaska Earthquakes: Implications for viscoelastic Earth structure

The Queen Charlotte-Fairweather Fault (QC-FF) system off the coast of British Columbia and southeast Alaska is a highly active dextral strike-slip plate boundary that accommodates ∼50 mm/yr of relative motion between the Pacific and North America plates. Nine M W ≥ 6.7 earthquakes have occurred along the QC-FF system since 1910, including a M S (G-R) 8.1 event in 1949. Two recent earthquakes, the October 28, 2012 Haida Gwaii ( M W 7.8) and January 5, 2013 Craig, Alaska ( M W 7.5) events, produced postseismic transient deformation that was recorded in the motions of 25 nearby continuous Global Positioning System (cGPS) stations. Here, we use 5+ yr of cGPS measurements to characterize the underlying mechanisms of postseismic deformation and to constrain the viscosity structure of the upper mantle surrounding the QC-FF. We construct forward models of viscoelastic deformation driven by coseismic stress changes from these two earthquakes and explore a large set of laterally heterogeneous viscosity structures that incorporate a relatively weak back-arc domain; we then evaluate each model based on its fit to the postseismic signals in our cGPS data. In determining best-fit model structures, we additionally incorporate the effects of afterslip following the 2012 event. Our results indicate the occurrence of a combination of temporally decaying afterslip and vigorous viscoelastic relaxation of the mantle asthenosphere. In addition, our best-fit viscosity structure (transient viscosity of 1.4–2.0 × 10 18 Pa s; steady-state viscosity of 10 19 Pa s) is consistent with the range of upper mantle viscosities determined in previous studies of glacial isostatic rebound and postseismic deformation.

Alaska, British Columbia↗

USGS1998

Closing the gap between regional and global travel time tomography

Recent global travel time tomography studies by Zhou [1996] and van der Hilst et al . [1997] have been performed with cell parameterizations of the order of those frequently used in regional tomography studies (i.e., with cell sizes of 1°–2°). These new global models constitute a considerable improvement over previous results that were obtained with rather coarse parameterizations (5° cells). The inferred structures are, however, of larger scale than is usually obtained in regional models, and it is not clear where and if individual cells are actually resolved. This study aims at resolving lateral heterogeneity on scales as small as 0.6° in the upper mantle and 1.2°–3° in the lower mantle. This allows for the adequate mapping of expected small-scale structures induced by, for example, lithosphere subduction, deep mantle upwellings, and mid-ocean ridges. There are three major contributions that allow for this advancement. First, we employ an irregular grid of nonoverlapping cells adapted to the heterogeneous sampling of the Earth's mantle by seismic waves [ Spakman and Bijwaard , 1998]. Second, we exploit the global data set of Engdahl et al . [1998], which is a reprocessed version of the global data set of the International Seismological Centre. Their reprocessing included hypocenter redetermination and phase reidentification. Finally, we combine all data used ( P , pP , and pwP phases) into nearly 5 million ray bundles with a limited spatial extent such that averaging over large mantle volumes is prevented while the signal-to-noise ratio is improved. In the approximate solution of the huge inverse problem we obtain a variance reduction of 57.1%. Synthetic sensitivity tests indicate horizontal resolution on the scale of the smallest cells (0.6° or 1.2°) in the shallow parts of subduction zones decreasing to approximately 2°–3° resolution in well-sampled regions in the lower mantle. Vertical resolution can be worse (up to several hundreds of kilometers) in subduction zones with rays predominantly pointing along dip. Important features of the solution are as follows: 100–200 km thick high-velocity slabs beneath all major subduction zones, sometimes flattening in the transition zone and sometimes directly penetrating into the lower mantle; large high-velocity anomalies in the lower mantle that have been attributed to subduction of the Tethys ocean and the Farallon plate; and low-velocity anomalies continuing across the 660 km discontinuity to hotspots at the surface under Iceland, east Africa, the Canary Islands, Yellowstone, and the Society Islands. Our findings corroborate that the 660 km boundary may resist but not prevent (present day) large-scale mass transfer from upper to lower mantle or vice versa. This observation confirms the results of previous, global mantle studies that employed coarser parameterizations.

Journal of Geophysical Research B: Solid Earth↗

USGS2012

Remarkable isotopic and trace element trends in potassic through sodic Cretaceous plutons of the Yukon-Koyukuk Basin, Alaska, and the nature of the lithosphere beneath the Koyukuk terrane

During the period from 110 to 80 m.y. ago, a 450-km-long magmatic belt was active along the northern margin of Yukon-Koyukuk basin and on eastern Seward Peninsula. The plutons intruded Upper Jurassic(?) and Lower Cretaceous volcanic arc rocks and Cretaceous sedimentary rocks in Yukon-Koyukuk basin and Proterozoic and lower Paleozoic continental rocks in Seward Peninsula. Within Yukon-Koyukuk basin, the plutons vary in composition from calc-alkalic plutons on the east to potassic and ultrapotassic alkalic plutons on the west. Plutons within Yukon-Koyukuk basin were analyzed for trace element and isotopic compositions in order to discern their origin and the nature of the underling lithosphere. Farthest to the east, the calc-alkalic rocks of Indian Mountain pluton are largely tonalite and sodic granodiorite, and have low Rb (average 82 ppm), high Sr (>600 ppm), high chondrite-normalized (cn) Ce/Yb (16–37), low δ 18 O (+6.5 to +7.1), low initial 87 Sr/ 86 Sr (SIR) (0.704), and high initial 143 Nd/ 144 Nd (NIR) (0.5126). These rocks resemble those modelled elsewhere as partial melts and subsequent fractionates of basaltic or gabbroic metaigneous rocks, and may be products of melting in the deeper parts of the Late Jurassic(?) and Early Cretaceous volcanic arc. Farthest to the west, the two ultrapotassic bodies of Selawik and Inland Lake are high in Cs (up to 93 ppm), Rb (up to 997 ppm), Sr, Ba, Th, and light rare earth elements, have high (Ce/Yb)cn (30, 27), moderate to low δ 18 O (+8.4, +6.9), high SIR (0.712, 0.710), and moderate NIR (0.5121–0.5122). These rocks resemble rocks of Australia and elsewhere that were modelled as melts of continental mantle that had been previously enriched in large cations. This mantle may be Paleozoic or older. The farthest west alkalic pluton of Selawik Hills is largely monzonite, quartz monzonite, and granite; has moderate Rb (average 284 ppm), high Sr (>600 ppm), high (Ce/Yb)cn (15–25), moderate δ 18 O (+8.3 to +8.6), high SIR (0.708–0.712), and moderate NIR (0.5121–0.5122). These rocks may be the product of interaction of magma derived from old continental mantle and magma derived from old continental crust. Plutons between eastern and western extremes show completely gradational variations in the concentration of K and Rb and in the isotopic compositions of Sr, Nd, and O. These plutons probably originated either by melting in a mixed source composed of a Paleozoic or older continental section (mantle + crust) overlain by Mesozoic mafic arc rocks, or by mixing of ultrapotassic to potassic magmas from continental sources (mantle + crust), and tonalitic magmas from arc sources. We infer from these results that the northwest portion of Yukon-Koyukuk basin is underlain by a substantial continental basement of Paleozoic or greater age. This basement probably thins out to the east. There is no geochemical evidence for continental basement east of about longitude 157°, or along a belt of at least 50 km width flanking Ruby Geanticline as far to the southwest as about longitude 161°. These areas are probably underlain by oceanic and Mesozoic arc rocks.

Alaska↗

USGS1996

Failure in laboratory fault models in triaxial tests

A model of a fault in the Earth is a sand-filled saw cut in a granite cylinder subjected to a triaxial test. The saw cut is inclined at an angle α to the cylinder axis, and the sand filling is intended to represent gouge. The triaxial test subjects the granite cylinder to a constant confining pressure and increasing axial stress to maintain a constant rate of shortening of the cylinder. The required axial stress increases at a decreasing rate to a maximum, beyond which a roughly constant axial stress is sufficient to maintain the constant rate of shortening. Such triaxial tests were run for saw cuts inclined at angles α of 20°, 25°, 30°, 35°, 40°, 45°, and 50° to the cylinder axis, and the apparent coefficient of friction μ a (ratio of the shear stress to the normal stress, both stresses resolved onto the saw cut) at failure was determined. Subject to the assumption that the observed failure involves slip on Coulomb shears (orientation unspecified), the orientation of the principal compression axis within the gouge can be calculated as a function of μ a for a given value of the coefficient of internal friction μ i . The rotation of the principal stress axes within the gouge in a triaxial test can then be followed as the shear strain across the gouge layer increases. For μ i ∼ 0.8, an appropriate value for highly sheared sand, the observed values μ a imply that the principal axis of compression within the gouge rotates so as to approach being parallel to the cylinder axis for all saw cut angles (20° < α < 50°). In the limiting state (principal compression axis parallel to cylinder axis) the stress state in the gouge layer would be the same as that in the granite cylinder, and the failure criterion would be independent of the saw cut angle.

Journal of Geophysical Research B: Solid Earth↗

USGS2013

Triggered tremor sweet spots in Alaska

To better understand what controls fault slip along plate boundaries, we have exploited the abundance of seismic and geodetic data available from the richly varied tectonic environments composing Alaska. A search for tremor triggered by 11 large earthquakes throughout all of seismically monitored Alaska reveals two tremor “sweet spots”—regions where large-amplitude seismic waves repeatedly triggered tremor between 2006 and 2012. The two sweet spots locate in very different tectonic environments—one just trenchward and between the Aleutian islands of Unalaska and Akutan and the other in central mainland Alaska. The Unalaska/Akutan spot corroborates previous evidence that the region is ripe for tremor, perhaps because it is located where plate-interface frictional properties transition between stick-slip and stably sliding in both the dip direction and laterally. The mainland sweet spot coincides with a region of complex and uncertain plate interactions, and where no slow slip events or major crustal faults have been noted previously. Analyses showed that larger triggering wave amplitudes, and perhaps lower frequencies (<~0.03 Hz), may enhance the probability of triggering tremor. However, neither the maximum amplitude in the time domain or in a particular frequency band, nor the geometric relationship of the wavefield to the tremor source faults alone ensures a high probability of triggering. Triggered tremor at the two sweet spots also does not occur during slow slip events visually detectable in GPS data, although slow slip below the detection threshold may have facilitated tremor triggering.

Alaska↗

USGS2013

Aftershocks of the 2010 M w 7.2 El Mayor-Cucapah earthquake revealcomplex faulting in the Yuha Desert, California

We detect and precisely locate over 9500 aftershocks that occurred in the Yuha Desert region during a 2 month period following the 4 April 2010 M w 7.2 El Mayor-Cucapah (EMC) earthquake. Events are relocated using a series of absolute and relative relocation procedures that include Hypoinverse, Velest, and hypoDD. Location errors are reduced to ~40 m horizontally and ~120 m vertically.Aftershock locations reveal a complex pattern of faulting with en echelon fault segments trending toward the northwest, approximately parallel to the North American-Pacific plate boundary and en echelon, conjugate features trending to the northeast. The relocated seismicity is highly correlated with published surface mapping of faults that experienced triggered surface slip in response to the EMC main shock. Aftershocks occurred between 2 km and 11 km depths, consistent with previous studies of seismogenic thickness in the region. Three-dimensional analysis reveals individual and intersecting fault planes that are limited in their along-strike length. These fault planes remain distinct structures at depth, indicative of conjugate faulting, and do not appear to coalesce onto a throughgoing fault segment. We observe a complex spatiotemporal migration of aftershocks, with seismicity that jumps between individual fault segments that are active for only a few days to weeks. Aftershock rates are roughly consistent with the expected earthquake production rates of Dieterich (1994). The conjugate pattern of faulting and nonuniform aftershock migration patterns suggest that strain in the Yuha Desert is being accommodated in a complex manner.

California↗

USGS2014

Maximum magnitude earthquakes induced by fluid injection

Analysis of numerous case histories of earthquake sequences induced by fluid injection at depth reveals that the maximum magnitude appears to be limited according to the total volume of fluid injected. Similarly, the maximum seismic moment seems to have an upper bound proportional to the total volume of injected fluid. Activities involving fluid injection include (1) hydraulic fracturing of shale formations or coal seams to extract gas and oil, (2) disposal of wastewater from these gas and oil activities by injection into deep aquifers, and (3) the development of enhanced geothermal systems by injecting water into hot, low-permeability rock. Of these three operations, wastewater disposal is observed to be associated with the largest earthquakes, with maximum magnitudes sometimes exceeding 5. To estimate the maximum earthquake that could be induced by a given fluid injection project, the rock mass is assumed to be fully saturated, brittle, to respond to injection with a sequence of earthquakes localized to the region weakened by the pore pressure increase of the injection operation and to have a Gutenberg-Richter magnitude distribution with a b value of 1. If these assumptions correctly describe the circumstances of the largest earthquake, then the maximum seismic moment is limited to the volume of injected liquid times the modulus of rigidity. Observations from the available case histories of earthquakes induced by fluid injection are consistent with this bound on seismic moment. In view of the uncertainties in this analysis, however, this should not be regarded as an absolute physical limit.

Journal of Geophysical Research B: Solid Earth↗

USGS2014

Observations of static Coulomb stress triggering of the November 2011 M5.7 Oklahoma earthquake sequence

In November 2011, a M 5.0 earthquake occurred less than a day before a M 5.7 earthquake near Prague, Oklahoma, which may have promoted failure of the mainshock and thousands of aftershocks along the Wilzetta fault, including a M 5.0 aftershock. The M 5.0 foreshock occurred in close proximity to active fluid injection wells; fluid injection can cause a buildup of pore fluid pressure, decrease the fault strength, and may induce earthquakes. Keranen et al . [ 2013 ] links the M 5.0 foreshock with fluid injection, but the relationship between the foreshock and successive events has not been investigated. Here we examine the role of coseismic Coulomb stress transfer on earthquakes that follow the M 5.0 foreshock, including the M 5.7 mainshock. We resolve the static Coulomb stress change onto the focal mechanism nodal plane that is most consistent with the rupture geometry of the three M ≥ 5.0 earthquakes, as well as specified receiver fault planes that reflect the regional stress orientation. We find that Coulomb stress is increased, e.g., fault failure is promoted, on the nodal planes of ~60% of the events that have focal mechanism solutions, and more specifically, that the M 5.0 foreshock promoted failure on the rupture plane of the M 5.7 mainshock. We test our results over a range of effective coefficient of friction values. Hence, we argue that the M 5.0 foreshock, induced by fluid injection, potentially triggered a cascading failure of earthquakes along the complex Wilzetta fault system.

Oklahoma↗

USGS2014

Crustal earthquake triggering by pre-historic great earthquakes on subduction zone thrusts

Triggering of earthquakes on upper plate faults during and shortly after recent great ( M> 8.0) subduction thrust earthquakes raises concerns about earthquake triggering following Cascadia subduction zone earthquakes. Of particular regard to Cascadia was the previously noted, but only qualitatively identified, clustering of M> ~6.5 crustal earthquakes in the Puget Sound region between about 1200–900 cal yr B.P. and the possibility that this was triggered by a great Cascadia thrust subduction thrust earthquake, and therefore portends future such clusters. We confirm quantitatively the extraordinary nature of the Puget Sound region crustal earthquake clustering between 1200–900 cal yr B.P., at least over the last 16,000. We conclude that this cluster was not triggered by the penultimate, and possibly full-margin, great Cascadia subduction thrust earthquake. However, we also show that the paleoseismic record for Cascadia is consistent with conclusions of our companion study of the global modern record outside Cascadia, that M> 8.6 subduction thrust events have a high probability of triggering at least one or more M> ~6.5 crustal earthquakes.

Cascadia;Puget Sound↗

USGS2014

Volcanic tremor masks its seismogenic source: Results from a study of noneruptive tremor recorded at Mount St. Helens, Washington

On 2 October 2004, a significant noneruptive tremor episode occurred during the buildup to the 2004–2008 eruption of Mount St. Helens (Washington). This episode was remarkable both because no explosion followed, and because seismicity abruptly stopped following the episode. This sequence motivated us to consider a model for volcanic tremor that does not involve energetic gas release from magma but does involve movement of conduit magma through extension on its way toward the surface. We found that the tremor signal was composed entirely of Love and Rayleigh waves and that its spectral bandwidth increased and decreased with signal amplitude, with broader bandwidth signals containing both higher and lower frequencies. Our modeling results demonstrate that the forces giving rise to this tremor were largely normal to conduit walls, generating hybrid head waves along conduit walls that are coupled to internally reflected waves. Together these form a crucial part of conduit resonance, giving tremor wavefields that are largely a function of waveguide geometry and velocity. We find that the mechanism of tremor generation fundamentally masks the nature of the seismogenic source giving rise to resonance. Thus multiple models can be invoked to explain volcanic tremor, requiring that information from other sources (such as visual observations, geodesy, geology, and gas geochemistry) be used to constrain source models. With concurrent GPS and field data supporting rapid rise of magma, we infer that tremor resulted from drag of nearly solid magma along rough conduit walls as magma was forced toward the surface.

Washington↗

USGS2014

A continuous record of intereruption velocity change at Mount St. Helens from coda wave interferometry

In September 2004, Mount St. Helens volcano erupted after nearly 18 years of quiescence. However, it is unclear from the limited geophysical observations when or if the magma chamber replenished following the 1980–1986 eruptions in the years before the 2004–2008 extrusive eruption. We use coda wave interferometry with repeating earthquakes to measure small changes in the velocity structure of Mount St. Helens volcano that might indicate magmatic intrusion. By combining observations of relative velocity changes from many closely located earthquake sources, we solve for a continuous function of velocity changes with time. We find that seasonal effects dominate the relative velocity changes. Seismicity rates and repeating earthquake occurrence also vary seasonally; therefore, velocity changes and seismicity are likely modulated by snow loading, fluid saturation, and/or changes in groundwater level. We estimate hydrologic effects impart stress changes on the order of tens of kilopascals within the upper 4 km, resulting in annual velocity variations of 0.5 to 1%. The largest nonseasonal change is a decrease in velocity at the time of the deep M w  = 6.8 Nisqually earthquake. We find no systematic velocity changes during the most likely times of intrusions, consistent with a lack of observable surface deformation. We conclude that if replenishing intrusions occurred, they did not alter seismic velocities where this technique is sensitive due to either their small size or the finite compressibility of the magma chamber. We interpret the observed velocity changes and shallow seasonal seismicity as a response to small stress changes in a shallow, pressurized system.

Journal of Geophysical Research B: Solid Earth↗

USGS2014

The 17 May 2012 M4.8 earthquake near Timpson, East Texas: An event possibly triggered by fluid injection

This study summarizes our investigation of the 17 May 2012 MW-RMT4.8 earthquake near Timpson, Texas, the largest earthquake recorded historically in eastern Texas. To identify preshocks and aftershocks of the 17 May event we examined the arrivals recorded at Nacogdoches (NATX) 30 km from the 17 May epicenter, at nearby USArray Transportable Array stations, and at eight temporary stations deployed between 26 May 2012 and mid-2013. At NATX we identified seven preshocks, the earliest occurring in April 2008. Reliably located aftershocks recorded by the temporary stations lie along a 6 km long NW-SE linear trend corresponding to a previously mapped basement fault that extends across the highest-intensity (MMI VII) region of the 17 May main shock. Earthquakes in this sequence are relatively shallow—with focal depths ranging from 1.6 to 4.6 km. Evidence supporting these depths include: hypocentral locations of exceptionally well-recorded aftershocks, S-P intervals at the nearest stations, and comparisons of synthetics and observed seismograms. Within 3 km of the linear trend of aftershock activity there are two Class II injection disposal wells injecting at 1.9 km depth beginning in August 2006 and February 2007, with injection rates averaging 42,750 m3/mo and 15,600 m3/mo, respectively. Several observations support the hypothesis that fluid injection triggered the Timpson sequence: well-located epicenters are situated near a mapped basement fault and near high-volume injection wells, focal depths are at or below the depths of injection, and the earliest preshock (April 2008) occurred after the onset of injection in 2006.

Texas↗

USGS2014

Three-dimensional seismic velocity structure of Mauna Loa and Kilauea volcanoes in Hawaii from local seismic tomography

We present a new three-dimensional seismic velocity model of the crustal and upper mantle structure for Mauna Loa and Kilauea volcanoes in Hawaii. Our model is derived from the first-arrival times of the compressional and shear waves from about 53,000 events on and near the Island of Hawaii between 1992 and 2009 recorded by the Hawaiian Volcano Observatory stations. The V p model generally agrees with previous studies, showing high-velocity anomalies near the calderas and rift zones and low-velocity anomalies in the fault systems. The most significant difference from previous models is in V p / V s structure. The high- V p and high- V p / V s anomalies below Mauna Loa caldera are interpreted as mafic magmatic cumulates. The observed low- V p and high- V p / V s bodies in the Kaoiki seismic zone between 5 and 15 km depth are attributed to the underlying volcaniclastic sediments. The high- V p and moderate- to low- V p / V s anomalies beneath Kilauea caldera can be explained by a combination of different mafic compositions, likely to be olivine-rich gabbro and dunite. The systematically low- V p and low- V p / V s bodies in the southeast flank of Kilauea may be caused by the presence of volatiles. Another difference between this study and previous ones is the improved V p model resolution in deeper layers, owing to the inclusion of events with large epicentral distances. The new velocity model is used to relocate the seismicity of Mauna Loa and Kilauea for improved absolute locations and ultimately to develop a high-precision earthquake catalog using waveform cross-correlation data.

Hawaii↗

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

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↗