Geology Reports⌕ Search

SEARCH · Geology Reports

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

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

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

At least 577 records · Page 32Linked to original sources

Chemistry of hydrothermal solutions from the southern Juan de Fuca Ridge

Fluids from three vent fields on the southern Juan de Fuca Ridge were sampled in September 1984 using the DSRV Alvin. The fluids are uniquely high in both chloride, which ranges up to twice the seawater value, and in metal content. Simple vapor-liquid phase separation could not have produced both the high chlorinity and gas concentrations observed in these fluids. The cause of the elevated chlorinity can not be uniquely identified but may be the result of either or a combination of two processes: (1) subsurface formation of a degassed brine and subsequent mixing of a small amount of this brine with a hydrothermal seawater which has not undergone a phase separation and (2) dissolution of a chloride-rich phase combined with a possible small loss of gas during sampling procedures. Although measured temperatures were all less than 300°C, quartz geothermometry suggests that the fluids have equilibrated at greater than 340°C. Quartz geobarometry is also in agreement with geophysical estimates of depth to the local magma chamber.

Journal of Geophysical Research B: Solid Earth↗

Friction of ice

The frictional strength of ice seems to be well below that for all other rocks. Triaxial testing of cylinders of pure water ice containing a 45° inclined sawcut, at temperatures of 77 ≤ T ≤ 115 K and confining pressures 0.1 ≤ P ≤ 250 MPa, reveals the frictional laws τ = 0.20 σ n + 8.3 MPa for P ≥ 10 MPa and τ = 0.55 σ n + 1.0 MPa for P ≤ 5 MPa, where τ and σ n are the shear and normal stresses on the sawcut surface. Friction is independent of T and, over the one order of magnitude tested, of average sliding velocity. The sliding behavior is invariably stick slip. Examination of sliding surfaces shows only minor gouge development. A curious arrangement of densely packed short vertical fractures occurs primarily in samples with anomalously low strength. These results apply to models of near‐surface tectonic activity on Jupiter's three large icy moons, Europa, Callisto, and Ganymede. They indicate, for example, that a global expansion on Ganymede of 3 linear percent will cause extensional movement on preexisting faults at depths to 7±3 km.

Journal of Geophysical Research B: Solid Earth↗

Structure of the lower crust beneath the Carolina Trough, U.S. Atlantic continental margin

Data from three large-offset seismic profiles provide information on the crustal structure beneath the Carolina trough. The profiles, obtained by the U.S. Geological Survey, the Naval Oceanographic Research Development Agency, and the Scripps Institution of Oceanography in 1985, were oriented parallel to the trough and were located (1) seaward of the East Coast Magnetic Anomaly (ECMA), which is generally thought to represent the boundary between oceanic and continental crust; (2) along the axis of the trough between the ECMA and the hinge zone, which is thought to reflect the landward limit of highly stretched and altered transitional crust; and (3) along the Carolina platform landward of the basement hinge zone on crust thought to have been thinned only slightly during rifting. These data constrain the velocity structure of the lower crust and provide evidence for a thick lens of high-velocity (>7.1 km/s) lower crustal material that extends beneath the Carolina trough and the adjacent ocean basin. This lens reaches a maximum thickness of about 13 km beneath the deepest part of the trough, thins to about 5 km seaward of the ECMA, and is either very thin or absent landward of the hinge zone. It is interpreted to represent material that was underplated beneath and/or intruded into the crust during the late stage of continental rifting and that led to an anomalously thick plutonic layer during the early seafloor spreading phase. These data thus support the recent conclusions of White et al. (1987b) and Mutter et al. (1988) that the initiation of seafloor spreading is attended in many, if not most, cases by the generation of an anomalously large volume of melt.

Carolina Trough↗

The competition between thermal contraction and differentiation in the stress history of the Moon

The scarcity of both extension and compression features on the Moon strongly constrains the history of the lunar radius—to variations of less than ±1 km over the past 3.8 Gyr. This limit has traditionally been interpreted as requiring a delicate balance between thermal contraction of the near‐surface and expansion of a substantial cold interior region. Recent theories of lunar origin (e.g., giant impact), in contrast, favor a “hot” initial state. We propose that a reconciliation may be possible by taking account of the volume change Δ V / V | d due to differentiation. We calculate STP densities based on simplified normative mineralogies for a suite of estimates of the bulk lunar composition, of primary lunar basalt, and of the residuum left when the maximum amount of the latter is extracted from the former. Typically Δ V / V | d ≃ 2 to 5%—an expansion equivalent to heating by ∼10 3 K. Provided the timing of differentiation is correct, one might offset the cooling of a magma ocean as much as 630 km deep by differentiation of the remainder of the Moon (which need not start much below the solidus temperature). A large but not impossible amount of gabbroic melt production is implied: ∼100 times the volume of mare basalts known to have been extruded. We do not address the detailed genetic relationship of this melt to the basalts observed on the lunar surface but point out that it need not have reached the surface directly or even have entered the crust in order for the expansion to have occurred. To assess the timing of melt formation, we investigate a simple conductive lunar thermal model which takes account of both Δ V / V | d and thermal contraction. Our initial state is characterized by a central temperature T c and a depth Z 0 above which the material (derived from the magma ocean) is already at the solidus and is not suceptible to volume changes upon further differentiation. We find a range of models satisfying the limits on radius increase and decrease. The hottest has T c = 1210 K, Z 0 = 400 km; without Δ V / V | d , we would need a larger or colder (or both) core, e.g., T c ≲ 700 K for Z 0 = 200–400 km, in agreement with previous investigators. Our modeling thus lends credence to the idea that the Moon could have been initially ≳50% molten (with the remainder relatively close to the solidus) and yet experienced little volume change over the last 3.8 Gyr.

Journal of Geophysical Research B: Solid Earth↗

Crustal structure of the Chugach Mountains, southern Alaska: A study of peg‐leg multiples from a low‐velocity zone

A seismic refraction profile recorded along the geologic strike of the Chugach Mountains in southern Alaska shows three upper crustal high‐velocity layers (6.9, 7.2, and 7.6 km/s) and a unique pattern of strongly focussed echelon arrivals to a distance of 225 km. The group velocity of the ensemble of echelon arrivals is 6.4 km/s. Modeling of this profile with the reflectivity method reveals that the echelon pattern is due to peg‐leg multiples generated from within a low‐velocity zone between the second and third upper crustal high‐velocity layers. The third high‐velocity layer (7.6 km/s) is underlain at 18 km depth by a pronounced low‐velocity zone that produces a seismic shadow zone wherein peg‐leg multiples are seen as echelon arrivals. The interpretation of these echelon arrivals as multiples supersedes an earlier interpretation which attributed them to successive primary reflections arising from alternating high‐ and low‐velocity layers. Synthetic seismogram modeling indicates that a low‐velocity zone with transitional upper and lower boundaries generates peg‐leg multiples as effectively as one with sharp boundaries. No PmP or Pn arrivals from the subducting oceanic Moho at 30 km depth beneath the western part of the line are observed on the long‐offset (90–225 km) data. This may be due to a lower crustal waveguide whose top is the high‐velocity (7.6 km/s) layer and whose base is the Moho. A deep (∼54 km) reflector is not affected by the waveguide and has been identified in the data. Although peg‐leg multiples have been interpreted on some long‐range refraction profiles that sound to upper mantle depths, the Chugach Mountains profile is one of the few crustal refraction profiles where peg‐leg multiples are clearly observed. This study indicates that multiple and converted phases may be more important in seismic refraction/wide‐angle reflection profiles than previously recognized.

Alaska↗

The crustal structure of the Wrangellia Terrane along the East Glenn Highway, eastern‐southern Alaska

Recently acquired seismic refraction data from eastern‐southern Alaska provide new information on the structure and composition of the Wrangellia and adjacent terranes. The data comprise a 160‐km‐long refraction profile along the East Glenn (Tok‐Cutoff) Highway that was collected as part of the U.S. Geological Survey's multidisciplinary Trans‐Alaska Crustal Transect program. The upper 3 km of the Wrangellia terrane and associated rocks is characterized by low compressional wave velocities (V p = 1.9, 3.3, 4.6, 5.6 km s −1 ) and high‐velocity gradients common to most onshore seismic refraction profiles. There is also clear seismic expression of the West Fork fault system as a steep, down‐to‐the‐southwest fault that separates the Peninsular terrane to the southwest and the metamorphic complex of Gulkana River to the northeast. In contrast, no seismic expression occurs for the Paxson Lake fault, which separates the Wrangellia terrane from the metamorphic complex of Gulkana River. Adjacent to the Denali fault, within the Wrangellia terrane, two high‐velocity bodies (V p = 6.6 km s −1 ) occur in the upper crust. One of these extends to ∼10‐km depth and correlates with a late Paleozoic dioritic complex, suggesting that the Wrangellia terrane is at least 10 km thick in this part of Alaska. From 5 to 23 km depth, the crust appears seismically homogeneous, with velocity increasing from V p = 6.2 to V p = 6.6 km s −1 . Beneath this level, the crust is less well resolved, although evidence exists for a low‐velocity zone between 23 and 26 km and a possible southwest dipping interface at 35 km. No identifiable mantle refraction or reflection is observed, possibly indicating a crust as thick as 55 km. The relatively low seismic velocities in the upper 23 km of the crust compare favorably with laboratory‐measured velocities on pelitic schists and intermediate‐composition plutonic rocks (granites and granodiorites), both of which are recognized in Wrangellia. We interpret the seismic velocities to indicate that silicic‐to‐intermediate‐composition rocks are important constituents of the basement of this part of Wrangellia. Geologic evidence indicates that the Alaskan part of the Wrangellia terrane is a Paleozoic and Mesozoic island arc: our seismic evidence indicates it may have been built mostly on continental crust as opposed to the fragment of Wrangellia from Vancouver Island which was probably built on oceanic crust.

Alaska↗

Tertiary basin development and tectonic implications, Whipple Detachment System, Colorado River Extensional Corridor, California and Arizona

Colorado River extensional corridor. In the Mohave Mountains and Aubrey Hills of Arizona and the eastern Whipple Mountains of California near Parker Dam, these deposits comprise four unconformity-bounded sequences composed of locally derived epiclastic and volcanic rocks and the Peach Springs Tuff. The three older sequences represent syntectonic units that were deposited coeval with detachment faulting, and the fourth is interpreted to be postextensional. The sequences are correlated between four fault-bounded regions, which are the remnants of four different depositional basins. Similar sequences can be correlated over broad areas of the extensional corridor despite the general lack of widespread units. The basins developed in about the same positions, relative to each other and to volcanic sources, as they occupy at present. This is shown by gradational changes of pre-Tertiary rock types between regions, systematic variations in the abundance of magmatic units, and correlative volcanic units that occur in two adjacent regions. The basins formed in the early Miocene from segmentation of the upper crust into blocks bounded by high-angle faults that trended both parallel and perpendicular to the direction of extension and which were terminated at middle crustal depths by a low-angle detachment fault. Extreme rotation of one large crustal block, which constitutes the central Mohave Mountains, is recorded by a major unconformity in the lower Miocene section of one basin. Because coeval sections that formed in adjoining basins do not record this rotation, the underlying crustal blocks must have been separated by transfer faults that allowed them to rotate independently. These proposed transfer faults are represented at present by major faults with trends that parallel the direction of extension on the Whipple detachment system.

Arizona, California↗

Lithology and evolution of the crust-mantle boundary region in the southwestern Basin and Range Province

Mantle and crustal xenoliths from volcanic rocks in the southwestern Basin and Range province and Colorado Plateau Transition Zone reveal histories of episodic magmatism and deformation that have profoundly influenced the crustal structure of this region. Seismic transects in this area show a strongly reflective Moho of generally low relief, which, in the area of modern transects, consists of a thin zone (<2 km thick) of short reflectors. The upper mantle is transparent and has a P n of 7.8–8.0 km/s similar to much of the western United States. A lower crustal zone, 2–13 km thick, has variable internal reflectivity and a relatively low velocity of 6.6–6.8 km/s. Upper mantle peridotite xenoliths show both ductile and brittle deformational features and have structures and compositions affected by magmatic intrusion; intrusions form complex dike systems and extensive zones of grain boundary infiltration in peridotite xenoliths. Whereas melt infiltration preceded and followed ductile deformation, brittle deformation, represented by closely spaced joint systems and faults, followed ductile deformation and is related to the youngest magmatic episodes. These structural characteristics and high uppermost mantle temperature (∼1000°C) may combine to explain the relatively low P n . Alternating layers of ductily deformed and undeformed peridotites, with or without igneous intrusions, may contribute to the reflectivity of the Moho. Lower crustal xenoliths are dominantly igneous-textured pyroxenites and mafic to intermediate gabbros identical to the dikes in peridotite xenoliths. The crustal xenoliths also commonly are jointed, and in addition many show partial melting and have abundant cavities that probably were filled with CO 2 -rich fluids. These rocks are interpreted as products of underplated magmas that were fed through the mantle dike systems and may represent the lowest crustal unit identified in the seismic records. The mafic compositions and high densities of the crustal xenoliths indicate that the low velocity of the lower crust may be caused in part by fracture systems, partial melts, and high temperatures. Garnet granulite xenoliths from a locality with no mantle peridotite xenoliths probably represent crust of the region before late Miocene extension. Felsic granulite xenoliths from two localities have velocities like those of the two lower crustal units identified seismically and could be present in the modern crust as unequilibrated remnants of old crust. The preferred model for the evolution of the lower lithosphere is one in which extension affects the upper mantle as well as the crust and is overlapped in time by multiple magmatic episodes. The earliest magmatic events preceded extension, and later events accompanied and followed extension.

Arizona, California, Nevada↗

Keweenaw hot spot: Geophysical evidence for a 1.1 Ga mantle plume beneath the Midcontinent Rift System

The Proterozoic Midcontinent Rift System of North America is remarkably similar to Phanerozoic rifted continental margins and flood basalt provinces. Like the younger analogues, the volcanism within this older rift can be explained by decompression melting and rapid extrusion of igneous material during lithospheric extension above a broad, asthenospheric, thermal anomaly which we call the Keweenaw hot spot. Great Lakes International Multidisciplinary Program on Crustal Evolution seismic reflection profiles constrain end-member models of melt thickness and stretching factors, which yield an inferred mantle potential temperature of 1500°–1570°C during rifting. Combined gravity modeling and subsidence calculations are consistent with stretching factors that reached 3 or 4 before rifting ceased, and much of the lower crust beneath the rift consists of relatively high density intruded or underplated synrift igneous material. The isotopic signature of Keweenawan volcanic rocks, presented in a companion paper by Nicholson and Shirey (this issue), is consistent with our model of passive rifting above an asthenospheric mantle plume.

Midcontinent Rift System↗

Flank tectonics of Martian volcanoes

On the flanks of Olympus Mons is a series of terraces, concentrically distributed around the caldera. Their morphology and location suggest that they could be thrust faults caused by compressional failure of the cone. In an attempt to understand the mechanism of faulting and the possible influences of the interior structure of Olympus Mons, we have constructed a numerical model for elastic stresses within a Martian volcano. In the absence of internal pressurization, the middle slopes of the cone are subjected to compressional stress, appropriate to the formation of thrust faults. These stresses for Olympus Mons are ∼250 MPa. If a vacant magma chamber is contained within the cone, the region of maximum compressional stress is extended toward the base of the cone. If the magma chamber is pressurized, extensional stresses occur at the summit and on the upper slopes of the cone. For a filled but unpressurized magma chamber, the observed positions of the faults agree well with the calculated region of high compressional stress. Three other volcanoes on Mars, Ascraeus Mons, Arsia Mons, and Pavonis Mons, possess similar terraces. Extending our analysis to other Martian volcanoes, we find that only these three and Olympus Mons have flank stresses that exceed the compressional failure strength of basalt, lending support to the view that the terraces on all four are thrust faults.

Journal of Geophysical Research B: Solid Earth↗

Dark materials in Valles Marineris: Indications of the style of volcanism and magmatism on Mars

Rifting on the equatorial canyon system of Valles Marineris provides a unique view of the interior of the Martian crust to depths reaching 7 km, exposing several in situ bedrock units which testify to past volcanic and magmatic processes on Mars. A thick, regionally extensive deposit observed in Coprates and Juventae chasmata is interpreted on the basis of spectral reflectance, erosional morphology, and tendency for eolian mobilization to be composed of mafic glass, possibly an ancient Martian analogue of the lunar terra mantling deposits. Spectral mapping suggests that the dark floor‐covering materials in the lower canyons are derived from this unit. A series of cliffs in the Ophir Chasma wall rock is interpreted to be exposures of resistant bedrock; the spectral signature of this massive and uniform unit most closely resembles that of terrestrial mafic rocks altered to or coated by crystalline hematite. Application of computer mapping techniques to probable young volcanic materials in the central troughs yields an inferred distribution of volcanic activity consistent with the interpretation of extrusion along faults near the margins of the canyon floors. This result supports the hypothesis that the valles originated through tectonic extension.

Journal of Geophysical Research B: Solid Earth↗

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↗