SEARCH · USGS Science

Results for “Journal of Seismology”

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 73 records · Page 4Linked to original sources

USGS2005

Earthquake fracture energy inferred from kinematic rupture models on extended faults

We estimate fracture energy on extended faults for several recent earthquakes by retrieving dynamic traction evolution at each point on the fault plane from slip history imaged by inverting ground motion waveforms. We define the breakdown work (Wb) as the excess of work over some minimum traction level achieved during slip. Wb is equivalent to "seismological" fracture energy (G) in previous investigations. Our numerical approach uses slip velocity as a boundary condition on the fault. We employ a three-dimensional finite difference algorithm to compute the dynamic traction evolution in the time domain during the earthquake rupture. We estimate Wb by calculating the scalar product between dynamic traction and slip velocity vectors. This approach does not require specifying a constitutive law and assuming dynamic traction to be collinear with slip velocity. If these vectors are not collinear, the inferred breakdown work depends on the initial traction level. We show that breakdown work depends on the square of slip. The spatial distribution of breakdown work in a single earthquake is strongly correlated with the slip distribution. Breakdown work density and its integral over the fault, breakdown energy, scale with seismic moment according to a power law (with exponent 0.59 and 1.18, respectively). Our estimates of breakdown work range between 4 ?? 105 and 2 ?? 107 J/m2 for earthquakes having moment magnitudes between 5.6 and 7.2. We also compare our inferred values with geologic surface energies. This comparison might suggest that breakdown work for large earthquakes goes primarily into heat production. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2012

Dynamic rupture modeling with laboratory-derived constitutive relations

A laboratory-derived state variable friction constitutive relation is used in the numerical simulation of the dynamic growth of an in-plane or mode II shear crack. According to this formulation, originally presented by J. H. Dieterich, frictional resistance varies with the logarithm of the slip rate and with the logarithm of the frictional state variable as identified by A. L. Ruina. Under conditions of steady sliding, the state variable is proportional to (slip rate) −1 . Following suddenly introduced increases in slip rate, the rate and state dependencies combine to produce behavior which resembles slip weakening. When rupture nucleation is artificially forced at fixed rupture velocity, rupture models calculated with the state variable friction in a uniformly distributed initial stress field closely resemble earlier rupture models calculated with a slip weakening fault constitutive relation. Additional rupture models are calculated in which rupture nucleation is achieved naturally, with numerical simulations of the quasi-static response of the fault leading to the onset of unstable, dynamic rupture. When rupture nucleation with the state variable friction law takes place naturally, a large fraction of the fault accelerates before accelerating slip is concentrated in what ultimately becomes the rupture nucleation patch. The state evolution accompanying this accelerating slip leads to higher average rupture speeds or a more rapid rupture acceleration to near P wave rupture speeds. Rupture models are also calculated for the seismological asperity problem, that is, the failure of a highly stressed fault patch surrounded by a region of zero stress drop. Dynamic overshoot of slip into the region of zero stress drop roughly agrees with a simple energy balance analysis; the final size of the rupture is proportional to the square of the size of the high stress patch. Earlier frictional stability analyses have led to the definition of a critical fault patch size for rupture nucleation. This critical patch size is generally different from critical crack lengths determined from crack tip energy balance considerations applied to a simpler slip weakening law. In the model calculations, dynamic rupture does not nucleate if the starting patch size is less than the critical patch size. This is consistent with the frictional stability analyses. Thus these model calculations suggest that dynamic rupture following a state variable friction relation is similar to that following a simpler fault slip weakening law. However, when modeling the full cycle of fault motions, rate-dependent frictional responses included in the state variable formulation are important at low slip rates associated with rupture nucleation. The critical rupture nucleation dimension appropriate for a slip weakening fault does not predict the critical nucleation dimension for a state variable fault.

Journal of Geophysical Research Solid Earth↗

USGS2009

Validation of the rupture properties of the 2001 Kunlun, China (Ms = 8.1), earthquake from seismological and geological observations

We determine the finite-fault slip distribution of the 2001 Kunlun earthquake ( M s = 8.1) by inverting teleseismic waveforms, as constrained by geological and remote sensing field observations. The spatial slip distribution along the 400-km-long fault was divided into five segments in accordance with geological observations. Forward modelling of regional surface waves was performed to estimate the variation of the speed of rupture propagation during faulting. For our modelling, the regional 1-D velocity structure was carefully constructed for each of six regional seismic stations using three events with magnitudes of 5.1–5.4 distributed along the ruptured portion of the Kunlun fault. Our result shows that the average rupture velocity is about 3.6 km s −1 , consistent with teleseismic long period wave modelling. The initial rupture was almost purely strike-slip with a rupture velocity of 1.9 km s −1 , increasing to 3.5 km s −1 in the second fault segment, and reaching a rupture velocity of about 6 km s −1 in the third segment and the fourth segment, where the maximum surface offset, with a broad fault zone, was observed. The rupture velocity decelerated to a value of 3.3 km s −1 in the fifth and final segment. Coseismic slip on the fault was concentrated between the surface and a depth of about 10 km. We infer that significant variations in rupture velocity and the observed fault segmentation are indicative of variations in strength along the interface of the Kunlun fault, as well as variations in fault geometry.

Qinghai Province↗

USGS2020

Towards understanding relationships between atmospheric pressure variations and long-period horizontal seismic data: A case study

Variations in atmospheric pressure have long been known to introduce noise in long-period (>10 s) seismic records. This noise can overwhelm signals of interest such as normal modes and surface waves. Generally, this noise is most pronounced on the horizontal components where it arises due to tilting of the seismometer in response to changes in atmospheric pressure. Several studies have suggested methodologies for correcting unwanted pressure-induced noise using collocated microbarograph records. However, how applicable these corrections are to varying geologic settings and installation types (e.g. vault versus post-hole) is unclear. Using coefficients obtained by solving for the residuals of these corrections, we can empirically determine the sensitivity of instruments in a specific location to the influences of pressure. To better understand how long-period, pressure-induced noise changes with time and emplacement, we examine horizontal seismic records along with barometric pressure at five different Global Seismographic Network stations, all with multiple broadband seismometers. We also analyse three Streckeisen STS-2 broadband seismometers, which are collocated with a microbarograph, at the Albuquerque Seismological Laboratory. We observe periods of high magnitude-squared-coherence (γ 2 -coherence; γ 2 > 0.8) between the seismic and pressure signals which fluctuate through time, frequency, and even between seismic instruments in the same vault. These observations suggest that these tilt-generated signals are highly sensitive to very local (<10 m) site effects. However, we find that in cases where instruments are not located at a large depth (<100 m), the pressure-induced noise is polarized in a nearly constant direction that is consistent with local topographic features or the geometry of the vault. We also find that borehole instruments at a large depth (>100 m) appear to be unaffected by pressure-loading mechanisms outlined by Sorrells (1971) but possibly by Newtonian attraction. Correlating the induced-noise polarization direction with times of high coherence, we work to identify sensors that are ultimately limited by pressure-induced horizontal noise as well as period bands that can benefit from pressure corrections. We find that while the situation is complex, each sensor appears to have its own unique response to pressure. Our findings suggest that we can determine empirical relationships between pressure and induced tilt on a case by case basis.

Geophysical Journal International↗

USGS2012

W phase source inversion for moderate to large earthquakes (1990-2010)

Rapid characterization of the earthquake source and of its effects is a growing field of interest. Until recently, it still took several hours to determine the first-order attributes of a great earthquake (e.g. M w ≥ 7.5), even in a well-instrumented region. The main limiting factors were data saturation, the interference of different phases and the time duration and spatial extent of the source rupture. To accelerate centroid moment tensor (CMT) determinations, we have developed a source inversion algorithm based on modelling of the W phase, a very long period phase (100–1000 s) arriving at the same time as the P wave. The purpose of this work is to finely tune and validate the algorithm for large-to-moderate-sized earthquakes using three components of W phase ground motion at teleseismic distances. To that end, the point source parameters of all M w ≥ 6.5 earthquakes that occurred between 1990 and 2010 (815 events) are determined using Federation of Digital Seismograph Networks, Global Seismographic Network broad-band stations and STS1 global virtual networks of the Incorporated Research Institutions for Seismology Data Management Center. For each event, a preliminary magnitude obtained from W phase amplitudes is used to estimate the initial moment rate function half duration and to define the corner frequencies of the passband filter that will be applied to the waveforms. Starting from these initial parameters, the seismic moment tensor is calculated using a preliminary location as a first approximation of the centroid. A full CMT inversion is then conducted for centroid timing and location determination. Comparisons with Harvard and Global CMT solutions highlight the robustness of W phase CMT solutions at teleseismic distances. The differences in M w rarely exceed 0.2 and the source mechanisms are very similar to one another. Difficulties arise when a target earthquake is shortly (e.g. within 10 hr) preceded by another large earthquake, which disturbs the waveforms of the target event. To deal with such difficult situations, we remove the perturbation caused by earlier disturbing events by subtracting the corresponding synthetics from the data. The CMT parameters for the disturbed event can then be retrieved using the residual seismograms. We also explore the feasibility of obtaining source parameters of smaller earthquakes in the range 6.0 ≤M w < 6.5. Results suggest that the W phase inversion can be implemented reliably for the majority of earthquakes of M w = 6 or larger.

Geophysical Journal International↗

USGS2014

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

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

Journal of Geophysical Research B: Solid Earth↗

USGS2010

Historical seismograms for unravelling a mysterious earthquake: The 1907 Sumatra Earthquake

History of instrumental seismology is short. Seismograms are available only for a little more than 100 years; high-quality seismograms are available only for the last 50 years and the seismological database is very limited in time. To extend the database, seismograms of old events are of vital importance. Many unusual earthquakes are known to have occurred, but their seismological characteristics are poorly known. The 1907 Sumatra earthquake is one of them (1907 January 4, M = 7.6). Gutenberg and Richter located this event in the outer-rise area of the Sunda arc. This earthquake is known to be anomalous because of its extensive tsunami, which is disproportionate of its magnitude. The tsunami affected the coastal areas over 950 km along the Sumatran coast. We investigated this earthquake using the historical seismograms we could collect from several seismological observatories. We examined the P -wave arrival times listed in the Strassburg Bulletin (1912) and other station bulletins. The scatter of the Observed−Computed traveltime residuals ranges from –30 to 30 s, too large to locate the event accurately. The uncertainty of the epicentre estimated from an S-P grid-search relocation study is at least 1° (∼110 km). We interpreted the Omori seismograms from Osaka, Mizusawa and Tokyo, and the Wiechert seismograms from Göttingen and Uppsala by comparing them with the seismograms simulated from modern broad-band seismograms of the 2002, 2008 and two 2010 Sumatra earthquakes which occurred near the 1907 earthquake. From the amplitude of Rayleigh waves recorded on the Omori seismograms we conclude that the magnitude of the 1907 earthquake at about 30 to 40 s is about 7.8 (i.e. 7.5 to 8.0). The SH waveforms recorded on the Göttingen and Uppsala seismograms suggest that the 1907 earthquake is a thrust earthquake at a shallow depth around 30 km. The most likely scenario is that the 1907 earthquake initiated on the subduction interface, and slowly ruptured up-dip into the shallow sediments and caused the extensive tsunami. Although their quantity and quality are limited, historical seismograms provide key quantitative information about old events that cannot be obtained otherwise. This underscores the importance of preserving historical seismograms.

Geophysical Journal International↗

USGS2022

Three Mw ≥ 4.7 earthquakes within the Changning (China) shale gas field ruptured shallow faults intersecting with hydraulic fracturing wells

From 2017 to 2019, three destructive earthquakes (27 January 2017 Mw 4.7, 16 December 2018 Mw 5.2, and 3 January 2019 Mw 4.8) occurred in the Changning shale gas field in the southwest Sichuan Basin, China. Previous seismological studies attributed these events to hydraulic fracturing (HF), but were unable to identify the causative seismogenic faults and their slip behaviors. Here, we use Sentinel-1 synthetic aperture radar data to measure surface deformation triggered by the three events and conduct geodetic inversions to characterize their rupture models. The resulting coseismic interferograms show prominent surface deformation with the maximum line-of-sight displacements of up to 4 cm. The inversion results show that all three earthquakes mainly ruptured sedimentary formations above the shale gas bed, in the upper 3 km of the crust, with slip magnitudes ranging from 8.5 to 15 cm, and stress drops ranging from ∼1.8 to ∼3.3 MPa. Their source faults intersect with horizontal HF wells, but do not root in the crystalline basement. Combined with the reported difficulty of increasing HF operation pressures prior to the three events, we argue that they were most likely induced by direct injection of pressurized fluids into the fault zones. Crustal deformation patterns inferred from regional topography and GPS velocities highlight that the Changning field is located within a triple junction region near the southeastern margin of the Tibetan Plateau with large deformation gradients; such conditions are not only favorable to the development of critically stressed faults, but also facilitate the occurrence of at least moderate magnitude earthquakes.

Changning shale gas field, Sichuan Basin↗

USGS2020

Relocated aftershocks and background seismicity in eastern Indonesia shed light on the 2018 Lombok and Palu earthquake sequences

High seismicity rates in eastern Indonesia occur due to the complex interaction of several tectonic plates which resulted in two deadly, destructive earthquake sequences that occurred in Lombok Island and the city of Palu, Sulawesi in 2018. The first sequence began in July with an M w 6.4 event near Lombok, culminating in an M w 7.0 event 8 d later. This was then followed by a nearby M w 6.9 event 12 d later. Approximately 1000 km to the northeast, a separate sequence began several weeks later near Palu where an M w 7.5 event occurred that triggered a tsunami. In this study, we present hypocentre relocations for both earthquake sequences as well as all other regional earthquakes in eastern Indonesia. The relocations were performed using a teleseismic double-difference relocation method and arrival times for P and S waves from stations at local, regional, and teleseismic distances. The catalogue and phase data were taken from the Agency for Meteorology, Climatology and Geophysics (BMKG) of Indonesia and the International Seismological Centre (ISC) for the period of April 2009 through November 2018. The relocated catalogue provides an improved view of seismicity in eastern Indonesia over the study period, sharpening locations and interpretations of seismogenic features throughout the region. In the Lombok area, the relocated earthquakes clearly show a backarc thrust to the north of the Sunda-Banda Arc transition zone. The relocated aftershocks show that the destructive M w 7.0 and M w 6.9 earthquakes of the Lombok sequence ruptured two different regions: The M w 7.0 earthquake propagated westward, whereas the M w 6.9 earthquake propagated eastward. The entire sequence of Lombok earthquakes was most likely started by the M w 6.4 event as the initial event or foreshock, which then triggered backarc thrusts on both sides. Several weeks later and far to the northeast, the M w 7.5 Palu earthquake occurred along the Palu-Koro Fault, filling a seismic gap that had not ruptured in an M w 6.0 event or larger since at least 1900. The distribution of aftershocks indicates that the northern part of the Palu-Koro Fault has lower relative seismicity rates than the southern part at shallow depths, and that off fault aftershocks are mostly located to the east of the Palu-Koro Fault.

Lombok, Palu↗

USGS2012

Evolution of the Precambrian lithosphere: Seismological and geochemical constraints

Several recent models of crustal evolution are based on the belief that the thickness of the continental crust is proportional to its age, with ancient crust being the thickest. A worldwide review of seismic structure contradicts this belief and falsifies these models, at least for the Archean. Proterozoic crust has a thickness of 40–55 km and a substantial high‐velocity (>7 km/s) layer at its base, while Archean crust is only 27–40 km thick (except at the site of younger rifts and collisional boundaries) and lacks the basal high‐velocity layer. Seismology also provides evidence that the lithosphere is thickest beneath Archean cratons, while diamond ages show that this lithospheric keel must have already existed in the Archean. Geochemical data also indicate significant differences between Archean and Proterozoic lithosphere. Major and trace element studies of sediments show a change in upper crustal composition between the Archean and Proterozoic. Archean rocks are depleted in Si and K and enriched in Na, Ca, and Mg. There is also a marked change in the Eu/Eu * ratio. Mantle xenoliths and continental flood basalts show that the mantle lithosphere beneath Archean crust is ultradepleted in FeO compared to that beneath post‐Archean crust. The secular change in the crust‐forming process is attributed to a decline in mantle temperature, leading to a change in the composition of the lithospheric mantle. The higher temperature of the Archean mantle led to the eruption of komatiitic lavas, producing a refractory lithospheric mantle which is ultradepleted in FeO and volatiles. The resultant lithospheric keel is intrinsically less dense than the surrounding mantle and thus not susceptible to delamination. It was sufficiently thick and cool for diamonds to form during the Archean. In contrast, Proterozoic crust developed above fertile mantle. The eruption of continental flood basalts and underplating of basaltic sills is attributed to subsequent heating and partial melting of the lithospheric mantle. Consequently, Proterozoic crust is thickened and has a high‐velocity basal layer.

Journal of Geophysical Research B: Solid Earth↗

USGS2017

Reexamination of the magnitudes for the 1906 and 1922 Chilean earthquakes using Japanese tsunami amplitudes: Implications for source depth constraints

Far-field tsunami records from the Japanese tide gauge network allow the reexamination of the moment magnitudes ( M w ) for the 1906 and 1922 Chilean earthquakes, which to date rely on limited information mainly from seismological observations alone. Tide gauges along the Japanese coast provide extensive records of tsunamis triggered by six great ( M w >8) Chilean earthquakes with instrumentally determined moment magnitudes. These tsunami records are used to explore the dependence of tsunami amplitudes in Japan on the parent earthquake magnitude of Chilean origin. Using the resulting regression parameters together with tide gauge amplitudes measured in Japan we estimate apparent moment magnitudes of M w 8.0–8.2 and M w 8.5–8.6 for the 1906 central and 1922 north-central Chile earthquakes. The large discrepancy of the 1906 magnitude estimated from the tsunami observed in Japan as compared with those previously determined from seismic waves ( M s 8.4) suggests a deeper than average source with reduced tsunami excitation. A deep dislocation along the Chilean megathrust would favor uplift of the coast rather than beneath the sea, giving rise to a smaller tsunami and producing effects consistent with those observed in 1906. The 1922 magnitude inferred from far-field tsunami amplitudes appear to better explain the large extent of damage and the destructive tsunami that were locally observed following the earthquake than the lower seismic magnitudes ( M s 8.3) that were likely affected by the well-known saturation effects. Thus, a repeat of the large 1922 earthquake poses seismic and tsunami hazards in a region identified as a mature seismic gap.

Journal of Geophysical Research B: Solid Earth↗

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↗

USGS2019

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

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

Nevada↗

USGS1997

Long-period seismicity at Redoubt Volcano, Alaska, 1989-1990 related to magma degassing

The mass of exsolved magmatic H2O is estimated and compared to the mass of superheated steam (25-50 Mtons) released through the resonating crack producing the December 13-14, 1989 swarm of long-period seismic events at Redoubt Volcano. Results indicate degassing of a H2O-CO2-SO2-saturated magma upon ascending from at least 12 km to 3-4 km beneath the crater as the source of the superheated steam. The mass of exsolved H2O (3.2-250 Mtons) is estimated from solubility diagrams of H2O-CO2-saturated silicate melts for the ascent history of the Redoubt magmas. Crystal size distribution, seismological, petrological, and geochemical data are used to constrain the ascent history of the two andesitic magmas prior to the eruption. Two stages of crystallization are inferred from crystal size distributions of plagioclase crystals in andesites erupted in December 1989. The first stage occurred 30-150 years before the eruption in both magmas and the second stage occurred at least 8 years and 15 years before the eruption in the dacitic andesite and rhyolitic andesite, respectively. The depths of crystallization are constrained from the spatial and temporal variations of volcano-tectonic earthquakes locations (Lahr et al., 1994) and from the P-wave and S-wave velocity structures (Benz et al., 1996). These data suggest that the rhyolitic andesite magma ascended to a depth of 7-8 km within at least 15 years of the eruption. Within at least 8 years of the eruption, the dacitic andesite magma migrated to a depth just below the other magma body where it resided until hours to days of the eruption. At this time, the dacitic andesite magma mixed with the rhyolitic andesite magma and established the reservoir for the eruption. Near the top of the reservoir, some of the mixed magma was displaced into fractures which extended 4-5 km toward the surface. This displaced magma created the eruption conduit and released the fluids related to the resonating crack. This scenario is consistent with the trends in major-and trace-element chemistry, and the stability of hornblende in the pre-eruption Redoubt magmas. It also provides a source for the SO2 and CO2 emissions measured during the eruption.

Journal of Volcanology and Geothermal Research↗

USGS2013

The 2011 Mw 7.1 Van (Eastern Turkey) earthquake

We use interferometric synthetic aperture radar (InSAR), body wave seismology, satellite imagery, and field observations to constrain the fault parameters of the M w 7.1 2011 Van (Eastern Turkey) reverse-slip earthquake, in the Turkish-Iranian plateau. Distributed slip models from elastic dislocation modeling of the InSAR surface displacements from ENVISAT and COSMO-SkyMed interferograms indicate up to 9 m of reverse and oblique slip on a pair of en echelon NW 40 °–54 ° dipping fault planes which have surface extensions projecting to just 10 km north of the city of Van. The slip remained buried and is relatively deep, with a centroid depth of 14 km, and the rupture reaching only within 8–9 km of the surface, consistent with the lack of significant ground rupture. The up-dip extension of this modeled WSW striking fault plane coincides with field observations of weak ground deformation seen on the western of the two fault segments and has a dip consistent with that seen at the surface in fault gouge exposed in Quaternary sediments. No significant coseismic slip is found in the upper 8 km of the crust above the main slip patches, except for a small region on the eastern segment potentially resulting from the M w 5.9 aftershock on the same day. We perform extensive resolution tests on the data to confirm the robustness of the observed slip deficit in the shallow crust. We resolve a steep gradient in displacement at the point where the planes of the two fault segments ends are inferred to abut at depth, possibly exerting some structural control on rupture extent.

Journal of Geophysical Research↗

USGS2014

Investigation of the high-frequency attenuation parameter, κ (kappa), from aftershocks of the 2010 M w 8.8 Maule, Chile earthquake

The Bío Bío region of Chile experienced a vigorous aftershock sequence following the 2010 February 27 M w 8.8 Maule earthquake. The immediate aftershock sequence was captured by two temporary seismic deployments: the Quake Catcher Network Rapid Aftershock Mobilization Program (QCN RAMP) and the Incorporated Research Institutions for Seismology CHile Aftershock Mobilization Program (IRIS CHAMP). Here, we use moderate to large aftershocks ( M L ≥ 4.0) occurring between 2010 March 1 and June 30 recorded by QCN RAMP and IRIS CHAMP stations to determine the spectral decay parameter, kappa ( κ ). First, we compare waveforms and κ estimates from the lower-resolution QCN stations to the IRIS CHAMP stations to ensure the QCN data are of sufficient quality. We find that QCN stations provide reasonable estimates of κ in comparison to traditional seismic sensors and provide valuable additional observations of local ground motion variation. Using data from both deployments, we investigate the variation in κ for the region to determine if κ is influenced primarily by local geological structure, path attenuation, or source properties (e.g. magnitude, mechanism and depth). Estimates of κ for the Bío Bío region range from 0.0022 to 0.0704 s with a mean of 0.0295 s and are in good agreement with κ values previously reported for similar tectonic environments. κ correlates with epicentral distance and, to a lesser degree, with source magnitude. We find little to no correlation between the site kappa, κ 0 , and mapped geology, although we were only able to compare the data to a low-resolution map of surficial geology. These results support an increasing number of studies that suggest κ observations can be attributed to a combination of source, path and site properties; additionally, measured κ are often highly scattered making it difficult to separate the contribution from each of these factors. Thus, our results suggest that contributions from the site, path and source should be carefully considered when interpreting κ values.

Bío Bío Region↗

USGS2002

Imaging the mantle beneath Iceland using integrated seismological techniques

Using a combination of body wave and surface wave data sets to reveal the mantle plume and plume head, this study presents a tomographic image of the mantle structure beneath Iceland to 400 km depth. Data comes primarily from the PASSCAL-HOTSPOT deployment of 30 broadband instruments over a period of 2 years, and is supplemented by data from the SIL and ICEMELT networks. Three sets of relative teleseismic body wave arrival times are generated through cross correlation: S and SKS arrivals at 0.03–0.1 Hz, and P and PKIKP arrivals at 0.03–0.1 and 0.8–2.0 Hz. Prior to inversion the crustal portion of the travel time anomalies is removed using the crustal model ICECRTb. This step has a significant effect on the mantle velocity variations imaged down to a depth of ∼250 km. Inversion of relative arrival times only provides information on lateral velocity variations. Surface waves are therefore used to provide absolute velocity information for the uppermost mantle beneath Iceland. The average wave number for the Love wave fundamental mode at 0.020 and 0.024 Hz is measured and used to invert for the average S velocity. Combination of the body wave and surface wave information reveals a predominantly horizontal low-velocity anomaly extending from the Moho down to ∼250 km depth, interpreted as a plume head. Below the plume head a near-cylindrical low-velocity anomaly with a radius of ∼100 km and peak V P and V S anomalies of −2% and −4%, respectively, extends down to the maximum depth of resolution at 400 km. Within the plume head, in the uppermost mantle above the core of the plume, there is a relatively high velocity with a maximum V P and V S anomaly of +2%. This high-velocity anomaly may be the result of the extreme degree of melt extraction necessary to generate the thick (46 km) crust in central Iceland. Comparison of the plume volumetric flux implied by our images, the crustal generation rate, and the degree of melting suggested by rare earth element inversions, suggests that (1) mantle material must be flowing horizontally away from the plume core faster than the overlying lithosphere and (2) the bulk of the plume material does not participate in melting beneath Iceland.

Journal of Geophysical Research B: Solid Earth↗

USGS1999

The thin hot plume beneath Iceland

We present the results of a seismological investigation of the frequency-dependent amplitude variations across Iceland using data from the HOTSPOT array currently deployed there. The array is composed of 30 broad-band PASSCAL instruments. We use the parameter t(*), defined in the usual manner from spectral ratios (Halderman and Davis 1991), to compare observed S-wave amplitude variations with those predicted due to both anelastic attenuation and diffraction effects. Four teleseismic events at a range of azimuths are used to measure t(*). A 2-D vertical cylindrical plume model with a Gaussian-shaped velocity anomaly is used to model the variations. That part of t(*) caused by attenuation was estimated by tracing a ray through IASP91, then superimposing our plume model velocity anomaly and calculating the path integral of 1/vQ. That part of t(*) caused by diffraction was estimated using a 2-D finite difference code to generate synthetic seismograms. The same spectral ratio technique used for the data was then used to extract a predicted t(*). The t(*) variations caused by anelastic attenuation are unable to account for the variations we observe, but those caused by diffraction do. We calculate the t(*) variations caused by diffraction for different plume models and obtain our best-fit plume, which exhibits good agreement between the observed and measured t(*). The best-fit plume model has a maximum S-velocity anomaly of - 12 per cent and falls to 1/e of its maximum at 100 km from the plume centre. This is narrower than previous estimates from seismic tomography, which are broadened and damped by the methods of tomography. This velocity model would suggest greater ray theoretical traveltime delays than observed. However, we find that for such a plume, wave-front healing effects at frequencies of 0.03-0.175 Hz (the frequency range used to pick S-wave arrivals) causes a 40 per cent reduction in traveltime delay, reducing the ray theoretical delay to that observed.

Geophysical Journal International↗