Search USGSSearch

SEARCH · Search USGS

Results for “Journal of Geophysical Research Solid Earth”

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

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

At least 811 records · Page 45Linked to original sources

Basal slip and mechanical anisotropy of biotite

The basal slip systems of biotite and their mechanical expressions have been investigated by shortening single crystals oriented to maximize and minimize shear stresses on (001). Samples loaded at 45° to (001) exhibit gentle external rotations associated with dislocation glide. High‐angle kink bands in these samples, unlike those developed in micas loaded parallel to (001), are limited to sample corners. Samples shortened perpendicular to (001) show no evidence of nonbasal slip and fail by fracture over all conditions tested. The mechanical response of biotite shortened at 45° to (001) is nearly perfectly elastic‐plastic; stress‐strain curves are characterized by a steep elastic slope, a sharply defined yield point, and continued deformation at low (mostly 1%. Stresses measured beyond the yield point are insensitive to confining pressure over the range 200 to 500 MPa and exhibit weak dependencies upon strain rate and temperature. Assuming an exponential relationship between differential stress σ and strain rate of the form , the data collected over strain rates and temperatures of 10−7 to 10−4 s−1 and 20° to 400°C, respectively, are best fit by an exponential constant α of 0.41±0.08 MPa‐1 and an activation energy of 82±13 kJ/mol. A power law fits the data equally well with = 18±4 and = 51±9 kJ/mol. Samples oriented favorably for slip in directions [100] and [110] are measurably weaker than those shortened at 45° to [010] and [310], consistent with the reported Burgers vectors 〈100〉, 1/2 〈110〉, and 1/2 〈110〉. The anisotropy of biotite is further revealed by contrasting these plastic strengths with results of samples deformed parallel and perpendicular to (001). Previous studies have shown that biotite loaded in the (001) plane is strong prior to the nucleation of kink bands. The strength of biotite shortened perpendicular to (001) exceeds that measured parallel to (001) and is pressure dependent. Application of the results to deformation within the continental crust suggests that biotite oriented favorably for slip is much weaker than most other silicates over a wide range of geologic conditions. Its presence within foliated rocks and shear zones may limit locally the stresses that can be supported.

Journal of Geophysical Research B: Solid Earth

The 40Ar/39Ar thermochronology of the eastern Mojave Desert, California, and adjacent western Arizona with implications for the evolution of metamorphic core complexes

Mesozoic thickening and Cenozoic extension resulted in the juxtaposition of upper and middle crustal rocks in the eastern Mojave Desert, southeastern California and western Arizona. The application of 40 Ar/ 39 Ar thermochronology to rocks in this region provides information about the timing and nature of thrusting, plutonism, metamorphism, denudation, and detachment faulting. The 40 Ar/ 39 Ar ages of 175 to 125 Ma from the Clipper, Piute, Turtle, Mohave, Bill Williams, and Hualapai Mountains are interpreted to be the result of a middle Mesozoic thermal event(s) caused by crustal thickening. The 40 Ar/ 39 Ar data from the Clipper and Piute Mountains suggest that this thermal event was followed by a period of cooling at rates of 1°–5°C/m.y. Orogenesis culminated during the Late Cretaceous when rocks exposed in the Old Woman‐Piute, Chemehuevi, and Sacramento Mountains attained temperatures >500°C which reset the K–Ar systems of minerals from Proterozoic rocks. High‐grade metamorphism in the Old Woman Mountains area was caused by the intrusion of the Old Woman‐Piute batholith at 73±1 Ma. Cooling rates following batholith emplacement in the Old Woman Mountains were ∼100°C/m.y. between 73 and 70 Ma and 5°–10°C/m.y. from 70 to ∼30 Ma. Between 65 and 25 Ma the entire eastern Mojave Desert underwent a period of cooling at a rate of 2°–10°C/m.y. By 30 Ma, rocks exposed in the Old Woman‐Piute, Marble, Ship, Clipper, and Turtle Mountains were below ∼100°C. The 40 Ar/ 39 Ar ages from the Sacramento Mountains suggest that mylonitization caused by the onset of regional extension occurred at 23±1 Ma. When extension started in the Chemehuevi Mountains, rocks exposed in the southwestern and northeastern portions of footwall to the Chemehuevi detachment fault were at ∼180°C and ∼350°C, respectively. This suggests that the exposed part of the Chemehuevi detachment fault initiated at a dip of 5°–30° or as a series of higher‐angle faults that cut to a depth of 10–12 km and were later rotated to their present dip. Unroofing of the footwalls to detachment faults in the Sacramento and Chemehuevi Mountains resulted in average cooling rates of 10°–50°C/m.y. between 22 and 15 Ma.

California

Crustal structure of a transform plate boundary: San Francisco Bay and the central California continental margin

Wide-angle seismic data collected during the Bay Area Seismic Imaging Experiment provide new glimpses of the deep structure of the San Francisco Bay Area Block and across the offshore continental margin. San Francisco Bay is underlain by a veneer (<300 m) of sediments, beneath which P wave velocities increase rapidly from 5.2 km/s to 6.0 km/s at 7 km depth, consistent with rocks of the Franciscan subduction assemblage. The base of the Franciscan at-15-18 km depth is marked by a strong wide-angle reflector, beneath which lies an 8- to 10-km-thick lower crust with an average velocity of 6.75??0.15 km/s. The lower crust of the Bay Area Block may be oceanic in origin, but its structure and reflectivity indicate that it has been modified by shearing and/or magmatic intrusion. Wide-angle reflections define two layers within the lower crust, with velocities of 6.4-6.6 km/s and 6.9-7.3 km/s. Prominent subhorizontal reflectivity observed at near-vertical incidence resides principally in the lowermost layer, the top of which corresponds to the "6-s reflector" of Brocher et al. [1994]. Rheological modeling suggests that the lower crust beneath the 6-s reflector is the weakest part of the lithosphere; the horizontal shear zone suggested by Furlong et al. [1989] to link the San Andreas and Hayward/Calaveras fault systems may actually be a broad zone of shear deformation occupying the lowermost crust. A transect across the continental margin from the paleotrench to the Hayward fault shows a deep crustal structure that is more complex than previously realized. Strong lateral variability in seismic velocity and wide-angle reflectivity suggests that crustal composition changes across major transcurrent fault systems. Pacific oceanic crust extends 40-50 km landward of the paleotrench but, contrary to prior models, probably does not continue beneath the Salinian Block, a Cretaceous arc complex that lies west of the San Andreas fault in the Bay Area. The thickness (10 km) and high lower-crustal velocity of Pacific oceanic crust suggest that it was underplated by magmatism associated with the nearby Pioneer seamount. The Salinian Block consists of a 15-km-thick layer of velocity 6.0-6.2 km/s overlying a 5-km-thick, high-velocity (7.0 km/s) lower crust that may be oceanic crust, Cretaceous arc-derived lower crust, or a magmatically underplated layer. The strong structural variability across the margin attests to the activity of strike-slip faulting prior to and during development of the transcurrent Pacific/North American plate boundary around 29 Ma. Copyright 1996 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth

Isotopic and trace element compositions of upper mantle and lower crustal xenoliths, Cima volcanic field, California: Implications for evolution of the subcontinental lithospheric mantle

Ultramafic and mafic xenoliths from the Cima volcanic field, southern California, provide evidence of episodic modification of the upper mantle and underplating of the crust beneath a portion of the southern Basin and Range province. The upper mantle xenoliths include spinel peridotite and anhydrous and hydrous pyroxenite, some cut by igneous-textured pyroxenite-gabbro veins and dikes and some by veins of amphibole ?? plagioclase. Igneous-textured pyroxenites and gabbros like the dike rocks also occur abundantly as isolated xenoliths inferred to represent underplated crust. Mineral and whole rock trace element compositions among and within the different groups of xenoliths are highly variable, reflecting multiple processes that include magma-mantle wall rock reactions, episodic intrusion and it filtration of basaltic melts of varied sources into the mantle wall rock, and fractionation. Nd, Sr, and Pb isotopic compositions mostly of clinopyroxene and plagioclase mineral separates show distinct differences between mantle xenoliths (??Nd = -5.7 to +3.4; 87Sr/86Sr = 0.7051 - 0.7073; 206Pb/204Pb = 19.045 - 19.195) and the igneous-textured xenoliths (??Nd = +7.7 to +11.7; 87Sr/86Sr = 0.7027 - 0.7036 with one carbonate-affected outlier at 0.7054; and 206Pb/204Pb = 18.751 - 19.068), so that they cannot be related. The igneous-textured pyroxenites and gabbros are similar in their isotopic compositions to the host basaltic rocks, which have ??Nd of+5.1 to +9.3; 87Sr/86Sr of 0.7028 - 0.7050, and 206Pb/204Pb of 18.685 - 21.050. The igneous-textured pyroxenites and gabbros are therefore inferred to be related to the host rocks as earlier cogenetic intrusions in the mantle and in the lower crust. Two samples of peridotite, one modally metasomatized by amphibole and the other by plagioclase, have isotopic compositions intermediate between the igneous-textured xenoliths and the mantle rock, suggesting mixing, but also derivation of the metasomatizing magmas from two separate and distinct sources. Sm-Nd two-mineral "isochrons" yield apparent ages for petrographically identical rocks believed to be coeval ranging from -0 to 113 ?? 26 Ma, indicating the unreliability of dating these rocks with this method. Amphibole and plagioclase megacrysts are isotopically like the host basalts and probably originate by mechanical breakup of veins comagmatic with the host basaltic rocks. Unlike other Basin and Range localities, Cima Cr-diopside group isotopic compositions do not overlap with those of the host basalts. Copyright 1997 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth

Fault slip and seismic moment of the 1700 Cascadia earthquake inferred from Japanese tsunami descriptions

The 1700 Cascadia earthquake attained moment magnitude 9 according to new estimates based on effects of its tsunami in Japan, computed coseismic seafloor deformation for hypothetical ruptures in Cascadia, and tsunami modeling in the Pacific Ocean. Reports of damage and flooding show that the 1700 Casscadia tsunami reached 1-5 m heights at seven shoreline sites in Japan. Three sets of estimated heights express uncertainty about location and depth of reported flooding, landward decline in tsunami heights from shorelines, and post-1700 land-level changes. We compare each set with tsunami heights computed from six Cascadia sources. Each source is vertical seafloor displacement calculated with a three-dimensional elastic dislocation model, for three sources the rupture extends the 1100 km length of the subduction zone and differs in width and shallow dip; for the other sources, ruptures of ordinary width extend 360-670 km. To compute tsunami waveforms, we use a linear long-wave approximation with a finite difference method, and we employ modern bathymetry with nearshore grid spacing as small as 0.4 km. The various combinations of Japanese tsunami heights and Cascadia sources give seismic moment of 1-9 ?? 1022 N m, equivalent to moment magnitude 8.7-9.2. This range excludes several unquantified uncertainties. The most likely earthquake, of moment magnitude 9.0, has 19 m of coseismic slip on an offshore, full-slip zone 1100 km long with linearly decreasing slip on a downdip partial-slip zone. The shorter rupture models require up to 40 m offshore slip and predict land-level changes inconsistent with coastal paleoseismological evidence. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth

Toggling of seismicity by the 1997 Kagoshima earthquake couplet: A demonstration of time-dependent stress transfer

Two M ??? 6 well-recorded strike-slip earthquakes struck just 4 km and 48 days apart in Kagoshima prefecture, Japan, in 1997, providing an opportunity to study earthquake interaction. Aftershocks are abundant where the Coulomb stress is calculated to have been increased by the first event, and they abruptly stop where the stress is dropped by the second event. This ability of the main shocks to toggle seismicity on and off argues that static stress changes play a major role in exciting aftershocks, whereas the dynamic Coulomb stresses, which should only promote seismicity, appear to play a secondary role. If true, the net stress changes from a sequence of earthquakes might be expected to govern the subsequent seismicity distribution. However, adding the stress changes from the two Kagoshima events does not fully capture the ensuing seismicity, such as its rate change, temporal decay, or migration away from the ends of the ruptures. We therefore implement a stress transfer model that incorporates rate/state friction, in which seismicity is treated as a sequence of independent nucleation events that are dependent on the fault slip, slip rate, and elapsed time since the last event. The model reproduces the temporal response of seismicity to successive stress changes, including toggling, decay, and aftershock migration. Nevertheless, the match of observed to predicted seismicity is quite imperfect, due perhaps to inadequate knowledge of several model parameters. However, to demonstrate the potential of this approach, we build a probabilistic forecast of larger earthquakes on the expected rate of small aftershocks, taking advantage of the large statistical sample the small shocks afford. Not surprisingly, such probabilities are highly time- and location-dependent: During the first decade after the main shocks, the seismicity rate and the chance of successive large shocks are about an order of magnitude higher than the background rate and are concentrated exclusively in the stress triggering zones. Copyright 2003 by the American Geophysical Uion.

Journal of Geophysical Research B: Solid Earth

Location of long-period events below Kilauea Volcano using seismic amplitudes and accurate relative relocation

We present methods for improving the location of long-period (LP) events, deep and shallow, recorded below Kilauea Volcano by the permanent seismic network. LP events might be of particular interest to understanding eruptive processes as their source mechanism is assumed to directly involve fluid transport. However, it is usually difficult or impossible to locate their source using traditional arrival time methods because of emergent wave arrivals. At Kilauea, similar LP waveform signatures suggest the existence of LP multiplets. The waveform similarity suggests spatially close sources, while catalog solutions using arrival time estimates are widely scattered beneath Kilauea's summit caldera. In order to improve estimates of absolute LP location, we use the distribution of seismic amplitudes corrected for station site effects. The decay of the amplitude as a function of hypocentral distance is used for inferring LP location. In a second stage, we use the similarity of the events to calculate their relative positions. The analysis of the entire LP seismicity recorded between January 1997 and December 1999 suggests that a very large part of the LP event population, both deep and shallow, is generated by a small number of compact sources. Deep events are systematically composed of a weak high-frequency onset followed by a low-frequency wave train. Aligning the low-frequency wave trains does not lead to aligning the onsets indicating the two parts of the signal are dissociated. This observation favors an interpretation in terms of triggering and resonance of a magmatic conduit. Instead of defining fault planes, the precise relocation of similar LP events, based on the alignment of the high-energy low-frequency wave trains, defines limited size volumes. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth

Theoretical analysis of tsunami generation by pyroclastic flows

Pyroclastic flows are a common product of explosive volcanism and have the potential to initiate tsunamis whenever thick, dense flows encounter bodies of water. We evaluate the process of tsunami generation by pyroclastic flow by decomposing the pyroclastic flow into two components, the dense underflow portion, which we term the pyroclastic debris flow, and the plume, which includes the surge and coignimbrite ash cloud parts of the flow. We consider five possible wave generation mechanisms. These mechanisms consist of steam explosion, pyroclastic debris flow, plume pressure, plume shear, and pressure impulse wave generation. Our theoretical analysis of tsunami generation by these mechanisms provides an estimate of tsunami features such as a characteristic wave amplitude and wavelength. We find that in most situations, tsunami generation is dominated by the pyroclastic debris flow component of a pyroclastic flow. This work presents information sufficient to construct tsunami sources for an arbitrary pyroclastic flow interacting with most bodies of water. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth

Groundwater flow, heat transport, and water table position within volcanic edifices: Implications for volcanic processes in the Cascade Range

The position of the water table within a volcanic edifice has significant implications for volcano hazards, geothermal energy, and epithermal mineralization. We have modified the HYDROTHERM numerical simulator to allow for a free-surface (water table) upper boundary condition and a wide range of recharge rates, heat input rates, and thermodynamic conditions representative of continental volcano-hydrothermal systems. An extensive set of simulations was performed on a hypothetical stratovolcano system with unconfined groundwater flow. Simulation results suggest that the permeability structure of the volcanic edifice and underlying material is the dominant control on water table elevation and the distribution of pressures, temperatures, and fluid phases at depth. When permeabilities are isotropic, water table elevation decreases with increasing heat flux and increases with increasing recharge, but when permeabilities are anisotropic, these effects can be much less pronounced. Several conditions facilitate the ascent of a hydrothermal plume into a volcanic edifice: a sufficient source of heat and magmatic volatiles at depth, strong buoyancy forces, and a relatively weak topography-driven flow system. Further, the plume must be connected to a deep heat source through a pathway with a time-averaged effective permeability ???1 ?? 10-16 m2, which may be maintained by frequent seismicity. Topography-driven flow may be retarded by low permeability in the edifice and/or the lack of precipitation recharge; in the latter case, the water table may be relatively deep. Simulation results were compared with observations from the Quaternary stratovolcanoes along the Cascade Range of the western United States to infer hydrothermal processes within the edifices. Extensive ice caps on many Cascade Range stratovolcanoes may restrict recharge on the summits and uppermost flanks. Both the simulation results and limited observational data allow for the possibility that the water table beneath the stratovolcanoes is relatively deep.

Journal of Geophysical Research B: Solid Earth

Location of seismic events and eruptive fissures on the Piton de la Fournaise volcano using seismic amplitudes

We present a method for locating the source of seismic events on Piton de la Fournaise. The method is based on seismic amplitudes corrected for station site effects using coda site amplification factors. Once corrected, the spatial distribution of amplitudes shows smooth and simple contours for many types of events, including rockfalls, long-period events and eruption tremor. On the basis of the simplicity of these distributions we develop inversion methods for locating their origins. To achieve this, the decrease of the amplitude as a function of the distance to the source is approximated by the decay either of surface or body waves in a homogeneous medium. The method is effective for locating rockfalls, long-period events, and eruption tremor sources. The sources of eruption tremor are usually found to be located at shallow depth and close to the eruptive fissures. Because of this, our method is a useful tool for locating fissures at the beginning of eruptions.

Journal of Geophysical Research B: Solid Earth

Near-field postseismic deformation associated with the 1992 Landers and 1999 Hector Mine, California, earthquakes

After the Landers earthquake (Mw = 7.3, 1992.489) a linear array of 10 monuments extending about 30 km N50??E on either side of the earthquake rupture plus a nearby offtrend reference monument were surveyed frequently by GPS until 2003.2. The array also spans the rupture of the subsequent Hector Mine earthquake (Mw = 7.1, 1999.792 . The pre-Landers velocities of monuments in the array relative to interior North America were estimated from earlier trilateration and very long baseline interferometry measurements. Except at the reference monument, the post-Landers velocities of the individual monuments in the array relaxed to their preseismic values within 4 years. Following the Hector Mine earthquake the velocities of the monuments relaxed to steady rates within 1 year. Those steady rates for the east components are about equal to the pre-Landers rates as is the steady rate for the north component of the one monument east of the Hector Mine rupture. However, the steady rates for the north components of the 10 monuments west of the rupture are systematically ???10 mm yr1 larger than the pre-Landers rates. The relaxation to a steady rate is approximately exponential with decay times of 0.50 ?? 0.10 year following the Landers earthquake and 0.32 ?? 0.18 year following the Hector Mine earthquake. The postearthquake motions of the Landers array following the Landers earthquake are not well approximated by the viscoelastic-coupling model of Pollitz et al. [2000]. A similar viscoelastic-coupling model [Pollitz et al., 2001] is more successful in representing the deformation after the Hector Mine earthquake.

Journal of Geophysical Research B: Solid Earth

Imaging the complexity of an active normal fault system: The 1997 Colfiorito (central Italy) case study

Six moderate magnitude earthquakes (5 < Mw < 6) ruptured normal fault segments of the southern sector of the North Apennine belt (central Italy) in the 1997 Colfiorito earthquake sequence. We study the progressive activation of adjacent and nearby parallel faults of this complex normal fault system using ???1650 earthquake locations obtained by applying a double-difference location method, using travel time picks and waveform cross-correlation measurements. The lateral extent of the fault segments range from 5 to 10 km and make up a broad, ???45 km long, NW trending fault system. The geometry of each segment is quite simple and consists of planar faults gently dipping toward SW with an average dip of 40??-45??. The fault planes are not listric but maintain a constant dip through the entire seismogenic volume, down to 8 km depth. We observe the activation of faults on the hanging wall and the absence of seismicity in the footwall of the structure. The observed fault segmentation appears to be due to the lateral heterogeneity of the upper crust: preexisting thrusts inherited from Neogene's compressional tectonic intersect the active normal faults and control their maximum length. The stress tensor obtained by inverting the six main shock focal mechanisms of the sequence is in agreement with the tectonic stress active in the inner chain of the Apennine, revealing a clear NE trending extension direction. Aftershock focal mechanisms show a consistent extensional kinematics, 70% of which are mechanically consistent with the main shock stress field.

Journal of Geophysical Research B: Solid Earth

New insights into Kilauea's volcano dynamics brought by large-scale relative relocation of microearthquakes

We investigated the microseismicity recorded in an active volcano to infer information concerning the volcano structure and long-term dynamics, by using relative relocations and focal mechanisms of microearthquakes. There were 32,000 earthquakes of the Mauna Loa and Kilauea volcanoes recorded by more than eight stations of the Hawaiian Volcano Observatory seismic network between 1988 and 1999. We studied 17,000 of these events and relocated more than 70%, with an accuracy ranging from 10 to 500 m. About 75% of these relocated events are located in the vicinity of subhorizontal decollement planes, at a depth of 8-11 km. However, the striking features revealed by these relocation results are steep southeast dipping fault planes working as reverse faults, clearly located below the decollement plane and which intersect it. If this decollement plane coincides with the pre-Mauna Loa seafloor, as hypothesized by numerous authors, such reverse faults rupture the pre-Mauna Loa oceanic crust. The weight of the volcano and pressure in the magma storage system are possible causes of these ruptures, fully compatible with the local stress tensor computed by Gillard et al. [1996]. Reverse faults are suspected of producing scarps revealed by kilometer-long horizontal slip-perpendicular lineations along the decollement surface and therefore large-scale roughness, asperities, and normal stress variations. These are capable of generating stick-slip, large-magnitude earthquakes, the spatial microseismic pattern observed in the south flank of Kilauea volcano, and Hilina-type instabilities. Rupture intersecting the decollement surface, causing its large-scale roughness, may be an important parameter controlling the growth of Hawaiian volcanoes.

Journal of Geophysical Research B: Solid Earth

When and where the aftershock activity was depressed: Contrasting decay patterns of the proximate large earthquakes in southern California

Seismic quiescence has attracted attention as a possible precursor to a large earthquake. However, sensitive detection of quiescence requires accurate modeling of normal aftershock activity. We apply the epidemic-type aftershock sequence (ETAS) model that is a natural extension of the modified Omori formula for aftershock decay, allowing further clusters (secondary aftershocks) within an aftershock sequence. The Hector Mine aftershock activity has been normal, relative to the decay predicted by the ETAS model during the 14 months of available data. In contrast, although the aftershock sequence of the 1992 Landers earthquake (M = 7.3), including the 1992 Big Bear earthquake (M = 6.4) and its aftershocks, fits very well to the ETAS up until about 6 months after the main shock, the activity showed clear lowering relative to the modeled rate (relative quiescence) and lasted nearly 7 years, leading up to the Hector Mine earthquake (M = 7.1) in 1999. Specifically, the relative quiescence occurred only in the shallow aftershock activity, down to depths of 5-6 km. The sequence of deeper events showed clear, normal aftershock activity well fitted to the ETAS throughout the whole period. We argue several physical explanations for these results. Among them, we strongly suspect aseismic slips within the Hector Mine rupture source that could inhibit the crustal relaxation process within "shadow zones" of the Coulomb's failure stress change. Furthermore, the aftershock activity of the 1992 Joshua Tree earthquake (M = 6.1) sharply lowered in the same day of the main shock, which can be explained by a similar scenario.

Journal of Geophysical Research B: Solid Earth

Irregular recurrence of paleoearthquakes along the central Garlock fault near El Paso Peaks, California

New investigations at the El Paso Peaks paleoseismic site have refined and extended the record of paleoearthquakes for the central Garlock fault. Event evidence, in the form of buried fissures, scarps, folds, and upward terminations, has allowed us to identify six well-resolved earthquakes, designated as events W, U, R, Q, K, and F, within the last ???7000 years. The exposures that we excavated did not reveal additional evidence for three poorly resolved events reported by McGill and Rockwell [1998]. Evidence for event Y, identified in the previous study, was reexamined and found not to be a paleoearthquake. Radiocarbon dates of detrital charcoal, combined with a method of interpolating the event ages with revised sedimentation rates, provide our preferred constraints on the timing of faulting events. The most recent surface-rupturing earthquake, event W, occurred between A.D. 1450 and 1640. The penultimate event U occurred between A.D. 675 and 950. Event R, which was unrecognized in the previous study, occurred between A.D. 250 and 475. Event Q occurred relatively shortly before event R, between A.D. 25 and 275. Event K occurred between 3340 and 2930 B.C. The oldest identified paleoearthquake, event F, occurred between 5300 and 4670 B.C. The event ages indicate that earthquake recurrence is highly irregular at the El Paso Peaks site, with individual preferred intervals ranging from as little as 215 years to as many as 3300 years. The irregular recurrence of surface-rupturing earthquakes along the central Garlock fault may be related to temporal and spatial clustering of earthquakes proposed for faults in the Eastern California Shear Zone.

Journal of Geophysical Research B: Solid Earth

Fault trends on the seaward slope of the Aleutian Trench: Implications for a laterally changing stress field tied to a westward increase in oblique convergence

Normal faults along the seaward trench slope (STS) commonly strike parallel to the trench in response to bending of the oceanic plate into the subduction zone. This is not the circumstance for the Aleutian Trench, where the direction of convergence gradually changes westward, from normal to transform motion. GLORIA side-scan sonar images document that the Aleutian STS is dominated by faults striking oblique to the trench, west of 179??E and east of 172??W. These images also show a pattern of east-west trending seafloor faults that are aligned parallel to the spreading fabric defined by magnetic anomalies. The stress-strain field along the STS is divided into two domains west and east, respectively, of 179??E. Over the western domain, STS faults and nodal planes of earthquakes are oriented oblique (9??-46??) to the trench axis and (69??-90??) to the magnetic fabric. West of 179??E, STS fault strikes change by 36?? from the E-W trend of STS where the trench-parallel slip gets larger than its orthogonal component of convergence. This rotation indicates that horizontal stresses along the western domain of the STS are deflected by the increasing obliquity in convergence. An analytical model supports the idea that strikes of STS faults result from a superposition of stresses associated with the dextral shear couple of the oblique convergence and stresses caused by plate bending. For the eastern domain, most nodal planes of earthquakes strike parallel to the outer rise, indicating bending as the prevailing mechanism causing normal faulting. East of 172??W, STS faults strike parallel to the magnetic fabric but oblique (10??-26??) to the axis of the trench. On the basis of a Coulomb failure criterion the trench-oblique strikes probably result from reactivation of crustal faults generated by spreading. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth

Fault interaction and stress triggering of twentieth century earthquakes in Mongolia

A cluster of exceptionally large earthquakes in the interior of Asia occurred from 1905 to 1967: the 1905 M7.9 Tsetserleg and M8.4 Bolnai earthquakes, the 1931 M8.0 Fu Yun earthquake, the 1957 M8.1 Gobi-Altai earthquake, and the 1967 M7.1 Mogod earthquake (sequence). Each of the larger (M ??? 8) earthquakes involved strike-slip faulting averaging more than 5 m and rupture lengths of several hundred kilometers. Available geologic data indicate that recurrence intervals on the major source faults are several thousands of years and distances of about 400 km separate the respective rupture areas. We propose that the occurrences of these and many smaller earthquakes are related and controlled to a large extent by stress changes generated by the compounded static deformation of the preceding earthquakes and subsequent viscoelastic relaxation of the lower crust and upper mantle beneath Mongolia. We employ a spherically layered viscoelastic model constrained by the 1994-2002 GPS velocity field in western Mongolia [Vergnolle et al., 2003]. Using the succession of twentieth century earthquakes as sources of deformation, we then analyze the time-dependent change in Coulomb failure stress (????f). At remote interaction distances, static ????f values are small. However, modeled postseismic stress changes typically accumulate to several tenths of a bar over time intervals of decades. Almost all significant twentieth century regional earthquakes (M ??? 6) with well-constrained fault geometry lie in positive ????f lobes of magnitude about +0.5 bar. Our results suggest that significant stress transfer is possible among continental faults separated by hundreds of kilometers and on timescales of decades. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth

The 2000 Nemuro-Hanto-Oki earthquake, off eastern Hokkaido, Japan, and the high intraslab seismic activity in the southwestern Kuril Trench

The 2000 Nemuro-Hanto-Oki earthquake (Mw6.8) occurred in the southwestern part of the Kuril Trench. The hypocenter was located close to the aftershock region of the great 1994 Kuril earthquake (Mw8.3), named "the 1994 Hokkaido-Toho-Oki earthquake" by the Japan Meteorological Agency, for which the fault plane is still in debate. Analysis of the 2000 event provides a clue to resolve the fault plane issue for the 1994 event. The hypocenters of the 2000 main shock and aftershocks are determined using arrival times from a combination of nearby inland and submarine seismic networks with an improved azimuthal coverage. They clearly show that the 2000 event was an intraslab event occurring on a shallow-dipping fault plane between 55 and 65 km in depth. The well-focused aftershock distribution of the 2000 event, the relative location of the 1994 event with respect to the 2000 event, and the similarity between their focal mechanisms strongly suggest that the faulting of the great 1994 earthquake also occurred on a shallow-dipping fault plane in the subducting slab. The recent hypocenter distribution around the 1994 aftershock region also supports this result. Large intraslab earthquakes occuring to the southeast of Hokkaido may occur due to a strong coupling on the plate boundary, which generates relatively large stress field within the subducting Pacific plate.

Journal of Geophysical Research B: Solid Earth