SEARCH · USGS Science

Results for “Journal of Geophysical Research - Atmospheres”

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 577 records · Page 32Linked to original sources

USGS2007

Sulfur species behavior in soil organic matter during decomposition

Soil organic matter (SOM) is a primary re??servoir of terrestrial sulfur (S), but its role in the global S cycle remains poorly understood. We examine S speciation by X-ray absorption near-edge structure (XANES) spectroscopy to describe S species behavior during SOM decomposition. Sulfur species in SOM were best represented by organic sulfide, sulfoxide, sulfonate, and sulfate. The highest fraction of S in litter was organic sulfide, but as decomposition progressed, relative fractions of sulfonate and sulfate generally increased. Over 6-month laboratory incubations, organic sulfide was most reactive, suggesting that a fraction of this species was associated with a highly labile pool of SOM. During humification, relative concentrations of sulfoxide consistently decreased, demonstrating the importance of sulfoxide as a reactive S phase in soil. Sulfonate fractional abundance increased during humification irrespective of litter type, illustrating its relative stability in soils. The proportion of S species did not differ systematically by litter type, but organic sulfide became less abundant in conifer SOM during decomposition, while sulfate fractional abundance increased. Conversely, deciduous SOM exhibited lesser or nonexistent shifts in organic sulfide and sulfate fractions during decomposition, possibly suggesting that S reactivity in deciduous litter is coupled to rapid C mineralization and independent of S speciation. All trends were consistent in soils across study sites. We conclude that S reactivity is related to spqciation in SOM, particularly in conifer forests, and S species fractions in SOM change, during decomposition. Our data highlight the importance of intermediate valence species (sulfoxide and sulfonate) in the pedochemical cycling of organic bound S. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research G: Biogeosciences↗

USGS2009

Three‐dimensional model of Hellenic Arc deformation and origin of the Cretan uplift

[1] The Hellenic Arc of Greece is the most seismically active part of Europe, but little is know about its mechanics. We modeled deformation along the arc using a finite element model. The model was intended to capture large‐scale 3‐D structure of Nubian plate subduction beneath the Aegean block and its deformational consequences. The shape of the interface was developed using mapped traces at the surface and earthquake hypocenters at depth. Model block motions were constrained by recent compilations of GPS velocity vectors. We simulated a 10 ka period of convergence between Nubia and the Aegean and calculated the strain field in the overriding plate as well as the spatial distribution and orientation of differential stress (∣ σ 1 − σ 3 ∣). From these calculations we derived testable quantities such as the expected seismic moment rate on the interplate contact, uplift pattern, and distribution of strain modes. Our relatively simple model broadly reproduced observed uplift patterns, earthquake activity, and loci of extension and contraction. The model showed a localization of uplift near the island of Crete, where the fastest Aegean uplift rates are well documented. Comparison of calculated expected seismic moment and observed earthquake catalogs implies a nearly fully coupled interplate contact. On the basis of our modeling results, we suggest that south Aegean deformation is driven primarily by the fast moving (∼33 mm a −1 ) Aegean upper plate overriding a nearly stalled (∼5 mm a −1 ) Nubian lower plate. This tectonic setting thus more closely resembles a continental thrust than it does a typical oceanic subduction zone.

Aegean Sea, Hellenic Arc↗

USGS2009

Premonitory acoustic emissions and stick-slip in natural and smooth-faulted Westerly granite

A stick-slip event was induced in a cylindrical sample of Westerly granite containing a preexisting natural fault by loading at constant confining pressure of 150 MPa. Continuously recorded acoustic emission (AE) data and computer tomography (CT)-generated images of the fault plane were combined to provide a detailed examination of microscale processes operating on the fault. The dynamic stick-slip event, considered to be a laboratory analog of an earthquake, generated an ultrasonic signal that was recorded as a large-amplitude AE event. First arrivals of this event were inverted to determine the nucleation site of slip, which is associated with a geometric asperity on the fault surface. CT images and AE locations suggest that a variety of asperities existed in the sample because of the intersection of branch or splay faults with the main fault. This experiment is compared with a stick-slip experiment on a sample prepared with a smooth, artificial saw-cut fault surface. Nearly a thousand times more AE were observed for the natural fault, which has a higher friction coefficient (0.78 compared to 0.53) and larger shear stress drop (140 compared to 68 MPa). However at the measured resolution, the ultrasonic signal emitted during slip initiation does not vary significantly between the two experiments, suggesting a similar dynamic rupture process. We propose that the natural faulted sample under triaxial compression provides a good laboratory analogue for a field-scale fault system in terms of the presence of asperities, fault surface heterogeneity, and interaction of branching faults. ?? 2009.

Journal of Geophysical Research B: Solid Earth↗

USGS2009

Inducing in situ, nonlinear soil response applying an active source

[1] It is well known that soil sites have a profound effect on ground motion during large earthquakes. The complex structure of soil deposits and the highly nonlinear constitutive behavior of soils largely control nonlinear site response at soil sites. Measurements of nonlinear soil response under natural conditions are critical to advancing our understanding of soil behavior during earthquakes. Many factors limit the use of earthquake observations to estimate nonlinear site response such that quantitative characterization of nonlinear behavior relies almost exclusively on laboratory experiments and modeling of wave propagation. Here we introduce a new method for in situ characterization of the nonlinear behavior of a natural soil formation using measurements obtained immediately adjacent to a large vibrator source. To our knowledge, we are the first group to propose and test such an approach. Employing a large, surface vibrator as a source, we measure the nonlinear behavior of the soil by incrementally increasing the source amplitude over a range of frequencies and monitoring changes in the output spectra. We apply a homodyne algorithm for measuring spectral amplitudes, which provides robust signal-to-noise ratios at the frequencies of interest. Spectral ratios are computed between the receivers and the source as well as receiver pairs located in an array adjacent to the source, providing the means to separate source and near-source nonlinearity from pervasive nonlinearity in the soil column. We find clear evidence of nonlinearity in significant decreases in the frequency of peak spectral ratios, corresponding to material softening with amplitude, observed across the array as the source amplitude is increased. The observed peak shifts are consistent with laboratory measurements of soil nonlinearity. Our results provide constraints for future numerical modeling studies of strong ground motion during earthquakes.

Journal of Geophysical Research B: Solid Earth↗

USGS2009

Tracking magma volume recovery at Okmok Volcano using GPS and an unscented kalman filter

Changes beneath a volcano can be observed through position changes in a GPS network, but distinguishing the source of site motion is not always straightforward. The records of continuous GPS sites provide a favorable data set for tracking magma migration. Dense campaign observations usually provide a better spatial picture of the overall deformation field, at the expense of an episodic temporal record. Combining these observations provides the best of both worlds. A Kalman filter provides a means for integrating discrete and continuous measurements and for interpreting subtle signals. The unscented Kalman filter (UKF) is a nonlinear method for time‐dependent observations. We demonstrate the application of this technique to deformation data by applying it to GPS data collected at Okmok volcano. Seven years of GPS observations at Okmok are analyzed using a Mogi source model and the UKF. The deformation source at Okmok is relatively stable at 2.5 km depth below sea level, located beneath the center of the caldera, which means the surface deformation is caused by changes in the strength of the source. During the 7 years of GPS observations more than 0.5 m of uplift has occurred, a majority of that during the time period January 2003 to July 2004. The total volume recovery at Okmok since the last eruption in 1997 is ∼60–80%. The UKF allows us to solve simultaneously for the time‐dependence of the source strength and for the location without a priori information about the source. 2009 by the American Geophysical Union.

Alaska↗

USGS2009

Seismic wave triggering of nonvolcanic tremor, episodic tremor and slip, and earthquakes on Vancouver Island

[1] We explore the physical conditions that enable triggering of nonvolcanic tremor and earthquakes by considering local seismic activity on Vancouver Island, British Columbia during and immediately after the arrival of large-amplitude seismic waves from 30 teleseismic and 17 regional or local earthquakes. We identify tremor triggered by four of the teleseismic earthquakes. The close temporal and spatial proximity of triggered tremor to ambient tremor and aseismic slip indicates that when a fault is close to or undergoing failure, it is particularly susceptible to triggering of further events. The amplitude of the triggering waves also influences the likelihood of triggering both tremor and earthquakes such that large amplitude waves triggered tremor in the absence of detectable aseismic slip or ambient tremor. Tremor and energy radiated from regional/local earthquakes share the same frequency passband so that tremor cannot be identified during these smaller, more frequent events. We confidently identify triggered local earthquakes following only one teleseism, that with the largest amplitude, and four regional or local events that generated vigorous aftershock sequences in their immediate vicinity. Earthquakes tend to be triggered in regions different from tremor and with high ambient seismicity rates. We also note an interesting possible correlation between large teleseismic events and episodic tremor and slip (ETS) episodes, whereby ETS events that are “late” and have built up more stress than normal are susceptible to triggering by the slight nudge of the shaking from a large, distant event, while ETS events that are “early” or “on time” are not.

Vancouver Island↗

USGS2009

Constraints on the stress state of the San Andreas fault with analysis based on core and cuttings from SAFOD drilling phases 1 and 2

Analysis of field data has led different investigators to conclude that the San Andreas Fault (SAF) has either anomalously low frictional sliding strength ( μ < 0.2) or strength consistent with standard laboratory tests ( μ > 0.6). Arguments for the apparent weakness of the SAF generally hinge on conceptual models involving intrinsically weak gouge or elevated pore pressure within the fault zone. Some models assert that weak gouge and/or high pore pressure exist under static conditions while others consider strength loss or fluid pressure increase due to rapid coseismic fault slip. The present paper is composed of three parts. First, we develop generalized equations, based on and consistent with the Rice (1992) fault zone model to relate stress orientation and magnitude to depth‐dependent coefficient of friction and pore pressure. Second, we present temperature‐ and pressure‐dependent friction measurements from wet illite‐rich fault gouge extracted from San Andreas Fault Observatory at Depth (SAFOD) phase 1 core samples and from weak minerals associated with the San Andreas Fault. Third, we reevaluate the state of stress on the San Andreas Fault in light of new constraints imposed by SAFOD borehole data. Pure talc ( μ ≈0.1) had the lowest strength considered and was sufficiently weak to satisfy weak fault heat flow and stress orientation constraints with hydrostatic pore pressure. Other fault gouges showed a systematic increase in strength with increasing temperature and pressure. In this case, heat flow and stress orientation constraints would require elevated pore pressure and, in some cases, fault zone pore pressure in excess of vertical stress.

California↗

USGS2009

Rupture parameters of the 2003 Zemmouri (M w 6.8), Algeria, earthquake from joint inversion of interferometric synthetic aperture radar, coastal uplift, and GPS

We study the surface deformation associated with the 21 May 2003 (M w = 6.8) Zemmouri (Algeria) earthquake, the strongest seismic event felt in the Algiers region since 1716. The thrust earthquake mechanism and related surface deformation revealed an average 0.50 m coastal uplift along ??55-km-long coastline. We obtain coseismic interferograms using Envisat advanced synthetic aperture radar (ASAR) (IS2) and RADARSAT standard beam (ST4) data from both the ascending and descending orbits of Envisat satellite, whereas the RADARSAT data proved useful only in the descending mode. While the two RADARSAT interferograms cover the earthquake area, Envisat data cover only the western half of the rupture zone. Although the interferometric synthetic aperture radar (InSAR) coherence in the epicenter area is poor, deformation fringes are observed along the coast in different patches. In the Boumerdes area, the maximum coseismic deformation is indicated by the high gradient of fringes visible in all interferograms in agreement with field measurements (tape, differential GPS, leveling, and GPS). To constrain the earthquake rupture parameters, we model the interferograms and uplift measurements using elastic dislocations on triangular fault patches in an elastic and homogeneous half-space. We invert the coseismic slip using first, a planar surface and second, a curved fault, both constructed from triangular elements using Poly3Dinv program that uses a damped least square minimization. The best fit of InSAR, coastal uplift, and GPS data corresponds to a 65-km-long fault rupture dipping 40?? to 50?? SE, located at 8 to 13 km offshore with a change in strike west of Boumerdes from N60??-65?? to N95??-105??. The inferred rupture geometry at depth correlates well with the seismological results and may have critical implications for the seismic hazard assessment of the Algiers region. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗

USGS2009

Volcano-tectonic implications of 3-D velocity structures derived from joint active and passive source tomography of the island of Hawaii

We present a velocity model of the onshore and offshore regions around the southern part of the island of Hawaii, including southern Mauna Kea, southeastern Hualalai, and the active volcanoes of Mauna Loa, and Kilauea, and Loihi seamount. The velocity model was inverted from about 200,000 first-arrival traveltime picks of earthquakes and air gun shots recorded at the Hawaiian Volcano Observatory (HVO). Reconstructed volcanic structures of the island provide us with an improved understanding of the volcano-tectonic evolution of Hawaiian volcanoes and their interactions. The summits and upper rift zones of the active volcanoes are characterized by high-velocity materials, correlated with intrusive magma cumulates. These high-velocity materials often do not extend the full lengths of the rift zones, suggesting that rift zone intrusions may be spatially limited. Seismicity tends to be localized seaward of the most active intrusive bodies. Low-velocity materials beneath parts of the active rift zones of Kilauea and Mauna Loa suggest discontinuous rift zone intrusives, possibly due to the presence of a preexisting volcanic edifice, e.g., along Mauna Loa beneath Kilauea's southwest rift zone, or alternatively, removal of high-velocity materials by large-scale landsliding, e.g., along Mauna Loa's western flank. Both locations also show increased seismicity that may result from edifice interactions or reactivation of buried faults. New high-velocity regions are recognized and suggest the presence of buried, and in some cases, previously unknown rift zones, within the northwest flank of Mauna Loa, and the south flanks of Mauna Loa, Hualalai, and Mauna Kea. Copyright 2009 by the American Geophysical Union.

Hawaii↗

USGS2008

Seismicity and deformation induced by magma accumulation at three basaltic volcanoes

We analyzed the evolution of volcano‐tectonic (VT) seismicity and deformation at three basaltic volcanoes (Kilauea, Mauna Loa, Piton de la Fournaise) during phases of magma accumulation. We observed that the VT earthquake activity displays an accelerating evolution at the three studied volcanoes during the time of magma accumulation. At the same times, deformation rates recorded at the summit of Kilauea and Mauna Loa volcanoes were not accelerating but rather tend to decay. To interpret these observations, we propose a physical model describing the evolution of pressure produced by the accumulation of magma into a reservoir. This variation of pressure is then used to force a simple model of damage, where damage episodes are equivalent to earthquakes. This model leads to an exponential increase of the VT activity and to an exponential decay of the deformation rate during accumulation phases. Seismicity and deformation data are well fitted by such an exponential model. The time constant, deduced from the exponential increase of the seismicity, is in agreement with the time constant predicted by the model of magma accumulation. This VT activity can thus be a direct indication of the accumulation of magma at depth, and therefore can be seen as a long‐term precursory phenomenon, at least for the three studied basaltic volcanoes. Unfortunately, it does not allow the prediction of the onset of future eruptions, as no diverging point (i.e., critical time) is present in the model.

Journal of Geophysical Research B: Solid Earth↗

USGS2009

A search in strainmeter data for slow slip associated with triggered and ambient tremor near Parkfield, California

We test the hypothesis that, as in subduction zones, slow slip facilitates triggered and ambient tremor in the transform boundary setting of California. Our study builds on the study of Peng et al. (2009) of triggered and ambient tremor near Parkfield, California during time intervals surrounding 31, potentially triggering, M ≥ 7.5 teleseismic earthquakes; waves from 10 of these triggered tremor and 29 occurred in periods of ambient tremor activity. We look for transient slow slip during 3-month windows that include 11 of these triggering and nontriggering teleseisms, using continuous strain data recorded on two borehole Gladwin tensor strainmeters (GTSM) located within the distribution of tremor epicenters. We model the GTSM data assuming only tidal and “drift” signals are present and find no detectable slow slip, either ongoing when the teleseismic waves passed or triggered by them. We infer a conservative detection threshold of about 5 nanostrain for abrupt changes and about twice this for slowly evolving signals. This could be lowered slightly by adding analyses of other data types, modeled slow slip signals, and GTSM data calibration. Detection of slow slip also depends on the slipping fault's location and size, which we describe in terms of equivalent earthquake moment magnitude, M. In the best case of the GTSM above a very shallow slipping fault, detectable slip events must exceed M~2, and if the slow slip is beneath the seismogenic zone (below ~15 km depth), even M~5 events are likely to remain hidden.

Journal of Geophysical Research B: Solid Earth↗

USGS2009

Kilauea slow slip events: Identification, source inversions, and relation to seismicity

Several slow slip events beneath the south flank of Kilauea Volcano, Hawaii, have been inferred from transient displacements in daily GPS positions. To search for smaller events that may be close to the noise level in the GPS time series, we compare displacement fields on Kilauea's south flank with displacement patterns in previously identified slow slip events. Matching displacement patterns are found for several new candidate events, although displacements are much smaller than previously identified events. One of the candidates, 29 May 2000, is coincident with a microearthquake swarm, as are all of the previously identified slow slip events. The microearthquakes follow the onset of slow slip, implying that they are triggered by stress changes during slip. The new slow slip event brings the total number of events on Kilauea, between 1997 and 2007, to eight, the smallest having M W = 5.3, and the largest having M W = 6.0. While the recurrence time between the four largest events is 2.11 ± 0.01 years, the repeat time for all eight events is 0.9 ± 0.6 years. We invert for the fault geometry and distribution of slip during the slow slip events. The optimal source depths of 5 km, assuming uniform slip dislocations in an elastic half‐space, are considerably shallower than the accompanying swarm earthquakes (6.5–8.5 km), which would place the earthquakes in a zone of decreased Coulomb stress. Inversions including the effects of topography and layered elastic structure in the forward models favor depths comparable to microearthquake depths, such that the earthquakes are located in a region of increased Coulomb stress. We also invert for time‐dependent fault slip directly from the 30 s GPS phase observations, constraining the source to the optimal uniform slip geometry. On the basis of these inversions, the larger events last between 1.5–2.2 days. The data are unable to resolve migration of slip along the fault. The temporal pattern of accompanying microearthquakes is consistent with the fault slip history assuming a seismicity rate theory based on rate and state‐friction, making the swarm earthquakes coshocks and aftershocks of the slow slip events.

Hawaii↗

USGS2009

Elastic wave speeds and moduli in polycrystalline ice Ih, si methane hydrate, and sll methane-ethane hydrate

We used ultrasonic pulse transmission to measure compressional, P, and shear, S, wave speeds in laboratory-formed polycrystalline ice Ih, si methane hydrate, and sll methane-ethane hydrate. From the wave speed's linear dependence on temperature and pressure and from the sample's calculated density, we derived expressions for bulk, shear, and compressional wave moduli and Poisson's ratio from -20 to 15??C and 22.4 to 32.8 MPa for ice Ih, -20 to 15??C and 30.5 to 97.7 MPa for si methane hydrate, and -20 to 10??C and 30.5 to 91.6 MPa for sll methane-ethane hydrate. All three materials had comparable P and S wave speeds and decreasing shear wave speeds with increasing applied pressure. Each material also showed evidence of rapid intergranular bonding, with a corresponding increase in wave speed, in response to pauses in sample deformation. There were also key differences. Resistance to uniaxial compaction, indicated by the pressure required to compact initially porous samples, was significantly lower for ice Ih than for either hydrate. The ice Ih shear modulus decreased with increasing pressure, in contrast to the increase measured in both hydrates ?? 2009.

Journal of Geophysical Research B: Solid Earth↗

USGS2009

Moderate rates of late Quaternary slip along the northwestern margin of the Basin and Range Province, Surprise Valley fault, northeastern California

The 86-km-long Surprise Valley normal fault forms part of the active northwestern margin of the Basin and Range province in northeastern California. We use trench mapping and radiocarbon, luminescence, and tephra dating to estimate displacements and timing of the past five surface-rupturing earthquakes on the central part of the fault near Cedarville. A Bayesian OxCal analysis of timing constraints indicates earthquake times of 18.2 ± 2.6, 10.9 ± 3.2, 8.5 ± 0.5, 5.8 ± 1.5, and 1.2 ± 0.1 ka. These data yield recurrence intervals of 7.3 ± 4.1, 2.5 ± 3.2, 2.7 ± 1.6, and 4.5 ± 1.5 ka and an elapsed time of 1.2 ± 0.1 ka since the latest surface-rupturing earthquake. Our best estimate of latest Quaternary vertical slip rate is 0.6 ?? 0.1 mm/a. This late Quaternary rate is remarkably similar to long-term (8-14 Ma) minimum vertical slip rates (>0.4-0.5 ± 0.3 mm/a) calculated from recently acquired seismic reflection and chronologic and structural data in Surprise Valley and the adjacent Warner Mountains. However, our slip rate yields estimates of extension that are lower than recent campaign GPS determinations by factors of 1.5-4 unless the fault has an unusually shallow (30°-35°) dip as suggested by recently acquired seismic reflection data. Coseismic displacements of 2-4.5 ± 1 m documented in the trench and probable rupture lengths of 53-65 km indicate a history of latest Quaternary earthquakes of M 6.8-7.3 on the central part of the. Surprise Valley fault.

California↗

USGS2009

Thermal conductivity of hydrate-bearing sediments

A thorough understanding of the thermal conductivity of hydrate-bearing sediments is necessary for evaluating phase transformation processes that would accompany energy production from gas hydrate deposits and for estimating regional heat flow based on the observed depth to the base of the gas hydrate stability zone. The coexistence of multiple phases (gas hydrate, liquid and gas pore fill, and solid sediment grains) and their complex spatial arrangement hinder the a priori prediction of the thermal conductivity of hydrate-bearing sediments. Previous studies have been unable to capture the full parameter space covered by variations in grain size, specific surface, degree of saturation, nature of pore filling material, and effective stress for hydrate-bearing samples. Here we report on systematic measurements of the thermal conductivity of air dry, water- and tetrohydrofuran (THF)-saturated, and THF hydrate–saturated sand and clay samples at vertical effective stress of 0.05 to 1 MPa (corresponding to depths as great as 100 m below seafloor). Results reveal that the bulk thermal conductivity of the samples in every case reflects a complex interplay among particle size, effective stress, porosity, and fluid-versus-hydrate filled pore spaces. The thermal conductivity of THF hydrate–bearing soils increases upon hydrate formation although the thermal conductivities of THF solution and THF hydrate are almost the same. Several mechanisms can contribute to this effect including cryogenic suction during hydrate crystal growth and the ensuing porosity reduction in the surrounding sediment, increased mean effective stress due to hydrate formation under zero lateral strain conditions, and decreased interface thermal impedance as grain-liquid interfaces are transformed into grain-hydrate interfaces.

Journal of Geophysical Research B: Solid Earth↗

USGS2009

Gas hydrate saturations estimated from fractured reservoir at Site NGHP-01-10, Krishna-Godavari Basin, India

During the Indian National Gas Hydrate Program Expedition 01 (NGHP-Ol), one of the richest marine gas hydrate accumulations was discovered at Site NGHP-01-10 in the Krishna-Godavari Basin. The occurrence of concentrated gas hydrate at this site is primarily controlled by the presence of fractures. Assuming the resistivity of gas hydratebearing sediments is isotropic, th?? conventional Archie analysis using the logging while drilling resistivity log yields gas hydrate saturations greater than 50% (as high as ???80%) of the pore space for the depth interval between ???25 and ???160 m below seafloor. On the other hand, gas hydrate saturations estimated from pressure cores from nearby wells were less than ???26% of the pore space. Although intrasite variability may contribute to the difference, the primary cause of the saturation difference is attributed to the anisotropic nature of the reservoir due to gas hydrate in high-angle fractures. Archie's law can be used to estimate gas hydrate saturations in anisotropic reservoir, with additional information such as elastic velocities to constrain Archie cementation parameters m and the saturation exponent n. Theory indicates that m and n depend on the direction of the measurement relative to fracture orientation, as well as depending on gas hydrate saturation. By using higher values of m and n in the resistivity analysis for fractured reservoirs, the difference between saturation estimates is significantly reduced, although a sizable difference remains. To better understand the nature of fractured reservoirs, wireline P and S wave velocities were also incorporated into the analysis.

Journal of Geophysical Research B: Solid Earth↗

USGS2009

Microseismicity at the North Anatolian Fault in the Sea of Marmara offshore Istanbul, NW Turkey

The North Anatolian Fault Zone (NAFZ) below the Sea of Marmara forms a “seismic gap” where a major earthquake is expected to occur in the near future. This segment of the fault lies between the 1912 Ganos and 1999 İzmit ruptures and is the only NAFZ segment that has not ruptured since 1766. To monitor the microseismic activity at the main fault branch offshore of Istanbul below the Çınarcık Basin, a permanent seismic array (PIRES) was installed on the two outermost Prince Islands, Yassiada and Sivriada, at a few kilometers distance to the fault. In addition, a temporary network of ocean bottom seismometers was deployed throughout the Çınarcık Basin. Slowness vectors are determined combining waveform cross correlation and P wave polarization. We jointly invert azimuth and traveltime observations for hypocenter determination and apply a bootstrap resampling technique to quantify the location precision. We observe seismicity rates of 20 events per month for M < 2.5 along the basin. The spatial distribution of hypocenters suggests that the two major fault branches bounding the depocenter below the Çınarcık Basin merge to one single master fault below ∼17 km depth. On the basis of a cross-correlation technique we group closely spaced earthquakes and determine composite focal mechanisms implementing recordings of surrounding permanent land stations. Fault plane solutions have a predominant right-lateral strike-slip mechanism, indicating that normal faulting along this part of the NAFZ plays a minor role. Toward the west we observe increasing components of thrust faulting. This supports the model of NW trending, dextral strike-slip motion along the northern and main branch of the NAFZ below the eastern Sea of Marmara.

Istanbul↗

USGS2009

On near-source earthquake triggering

[1] When one earthquake triggers others nearby, what connects them? Two processes are observed: static stress change from fault offset and dynamic stress changes from passing seismic waves. In the near-source region ( r ≤ 50 km for M ∼ 5 sources) both processes may be operating, and since both mechanisms are expected to raise earthquake rates, it is difficult to isolate them. We thus compare explosions with earthquakes because only earthquakes cause significant static stress changes. We find that large explosions at the Nevada Test Site do not trigger earthquakes at rates comparable to similar magnitude earthquakes. Surface waves are associated with regional and long-range dynamic triggering, but we note that surface waves with low enough frequency to penetrate to depths where most aftershocks of the 1992 M = 5.7 Little Skull Mountain main shock occurred (∼12 km) would not have developed significant amplitude within a 50-km radius. We therefore focus on the best candidate phases to cause local dynamic triggering, direct waves that pass through observed near-source aftershock clusters. We examine these phases, which arrived at the nearest (200–270 km) broadband station before the surface wave train and could thus be isolated for study. Direct comparison of spectral amplitudes of presurface wave arrivals shows that M ∼ 5 explosions and earthquakes deliver the same peak dynamic stresses into the near-source crust. We conclude that a static stress change model can readily explain observed aftershock patterns, whereas it is difficult to attribute near-source triggering to a dynamic process because of the dearth of aftershocks near large explosions.

Journal of Geophysical Research B: Solid Earth↗