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 199 records · Page 11Linked to original sources

Tectonic setting of the Yukon-Koyukuk basin and its borderlands, western Alaska

The Yukon-Koyukuk basin of western Alaska is composed of an arcuate belt of Jurassic and Early Cretaceous subduction-related volcanic and plutonic rocks (Koyukuk terrane) flanked by deep subbasins filled with mid-Cretaceous terrigenous sedimentary rocks. The basin is bordered on three sides by metamorphosed Proterozoic and Paleozoic continental rocks (Seward, Arctic Alaska, and Ruby terranes) and is separated from the metamorphic borderlands by a narrow, highly tectonized belt of oceanic crust and mantle rocks (composite Angayucham-Tozitna terrane). The oceanic and mantle rocks, which dip inward beneath the basin and are thrust outward onto the borderlands, are divided into three separate thrust panels: (1) a structurally lowest panel (Slate Creek) composed of phyllite and metagraywacke of probable Devonian age, (2) a middle panel (Narvak) composed of imbricated basalt, chert, and gabbro of Devonian to Early Jurassic age, and (3) a structurally highest panel (Kanuti) composed of gabbro and peridotite of probable Middle and Late Jurassic age. The three thrust panels appear to represent a reversely stacked sequence that progresses from continental slope deposits in the lower panel to cumulus and mantle peridotites in the upper. Metamorphic mineral K-Ar ages from garnet amphibolite on the sole of the upper panel suggest that the upper panel was emplaced on the middle panel in the Middle to Late Jurassic. Subsequent accretion of all three panels to the continental rocks of the borderlands occurred in the latest Jurassic and Early Cretaceous, synchronous with arc volcanism within the basin. Arc volcanism waned and the accretionary phase ended in the middle of Early Cretaceous time. Uplift and erosion of the metamorphic borderlands and the obducted oceanic rocks began in late Early Cretaceous and was accompanied by the rapid filling of two flanking subbasins with turbiditic sediments. In the latest Early Cretaceous and early Late Cretaceous, shallow marine and nonmarine conglomerates were deposited around the margins of the basin, and a prograding delta was built out from the southeast margin of the basin across the turbiditic subbasins and the remnant volcanic arc. In the Late Cretaceous, western Alaska was subjected to strong east-west compression which severely deformed both the Yukon-Koyukuk basin and the borderlands. Several widespread magmatic episodes in the mid- and Late Cretaceous and in early Tertiary transgress the basin boundaries and stitch together the accreted arc and oceanic terranes and the continental borderlands.

Journal of Geophysical Research Solid Earth

The Kanuti ophiolite, Alaska

The Kanuti ophiolite is a mafic-ultramafic thrust sheet of probable Jurassic age, formerly considered to be the upper part of the Yukon-Koyukuk ophiolite belt (Angayucham terrane). It is here called the Kanuti ophiolite after the Kanuti River region on the southeastern flank of the Yukon-Koyukuk Basin. The thrust sheet crops out discontinuously for a distance of more than 900 km along the northern and southeastern margins of the basin. It is probably correlative with similar ophiolite thrust sheets to the north in the western Brooks Range and to the south in the Ruby geanticline. Technically, the ophiolite is considered to be the Kanuti thrust panel of the Angayucham-Totzitna terrane. The Kanuti consistently overlies another extensive thrust sheet, consisting mostly of pillow basalt and radiolarian chert of Devonian to Jurassic age (Narvak thrust panel). This sheet is thrust over a third sheet consisting of probable Devonian phyllite and metagraywacke, which is in turn thrust over older metamorphic rocks (Slate Creek thrust panel). The Kanuti ophiolite is a partial ophiolite that consists of a lower residual mantle suite and an upper magmatic suite, but dikes, extrusives, and sediments are absent. The residual mantle suite is composed of harzburgite and dunite with refractory mineral compositions. The harzburgite is attributed to partial melting and extraction of basaltic magma; residual dunite is attributed to partial melting or to reaction of orthopyroxene out of harzburgite in contact with ascending melt diapirs. The magmatic suite consists of layered ultramafic and gabbroic rocks, containing minerals having limited iron enrichment. The absence of large volumes of magmatic rocks intermediate in composition between cumulus ultramafics and evolved gabbros favors periodic introduction of magma, rather than closed system fractional crystallization. The ultramafic rocks of both the residual mantle and magmatic suites are tectonites, which have undergone high-temperature deformation involving isoclinal folding on all scales and related syntectonic recrystallization. The olivine fabric is consistent with the glide system {0kl} [100], which has been produced experimentally at 800°–1190°C at 20 kbar. Olivine Z axes and subparallel isoclinal fold axes have consistent, northeast trends throughout the Kanuti region (>100 km NE-SW) and may be close to the original upper mantle flow direction, despite later low-angle thrust faulting. The order of crystallization in the cumulus ultramafic rocks of the magmatic suite is olivine, clinopyroxene, plagioclase, and orthopyroxene. The high Mg numbers of clinopyroxene (0.85–0.93) coexisting with olivine suggest that the cumulus ultramafic rocks crystallized at relatively high pressures (>10 kbar). The effects of parental magma composition cannot be evaluated, but the small difference in Mg numbers of coexisting olivine and clinopyroxene in the cumulus ultramafic rocks and in residual harzburgite suggests that regardless of absolute pressures, the pressure difference between the melting that produced the basalt magma and the initial fractional crystallization of the magma is small. Because of the limited range in rock types in the ophiolite, the tectonic environment cannot be interpreted unambiguously. However, the structural and petrological data are best reconciled with an origin in a volcanic arc tectonic setting.

Journal of Geophysical Research Solid Earth

Historic surface slip along the San Andreas Fault near Parkfield, California

The Parkfield Earthquake Prediction Experiment is focusing close attention on the 44-km-long section of the San Andreas fault that last ruptured seismically in 1966 ( M s 6.0). The 20-km-long central segment of the 1966 Parkfield rupture, extending from the mainshock epicenter at Middle Mountain southeastward to Gold Hill, forms a 1- to 2-km salient northeastward away from the dominant N40°W strike. Following the 1966 earthquake afterslip, aseismic slip has been nearly constant. Moderate Parkfield earthquakes have recurred on average every 21 years since 1857, when a great earthquake ( M ≈8) ruptured at least as far north as the southern Parkfield segment. Many measurements of slip have been made near Parkfield since 1966. Nevertheless, much of the history of surface slip remained uncertain, especially the total amount associated with the 1966 event. In 1985 we measured accumulated slip on the four oldest cultural features offset by the fault along the 1966 Parkfield rupture segment. We interpret net slip on each feature as a sum of event slip (sum of coseismic and rapid preseismic and postseismic slip) from Parkfield earthquakes and steady interseismic slip as measured over the last 20 years on nearby alinement arrays, creep meters, and trilateration lines. We assumed for each site that event slip was identical for the 1922, 1934, and 1966 Parkfield events and that long-term average rate of interseismic slip was constant between all events. Two fences on the southern segment, southeast of Gold Hill, indicate event slip of 13 and 15 cm and interseismic slip rate of 0.36 and 0.30 cm/yr since 1959 and 1908, respectively. At these sites, redundant independent data support our assumption that both event and interseismic slip occur uniformly. On the central segment, near Parkfield, both the 1934 and the 1966 ruptures offset a bridge built in 1932. Interseismic slip rate near the bridge has been about 1.1 cm/yr since 1966; thus we deduce an average event slip of 31 cm for the 1934 and 1966 earthquakes. On a parallel fault trace, 1 km to the southwest, slip was about 8 cm in 1966; thus total event slip summed across the entire fault zone near Parkfield was nearly 40 cm. On Middle Mountain, 4 km north of the 1966 mainshock epicenter, an offset fence indicates 17 cm of slip in 1966 and a 2.26-cm/yr interseismic slip rate since circa 1946. Thus the central segment of the 1966 rupture is characterized by much larger event slip (∼40 cm) than both distal segments (∼15 cm). This amount of surface slip per event is about twice what had been previously assumed. Larger 1966 surface slip in the central part of the rupture is geodetically compatible with a coseismic slip of 65±10 cm slip on a narrow, buried asperity between Middle Mountain and Gold Hill that has been inferred from the depth distribution of early aftershocks. Assuming our characteristic surface slip model, one can further deduce a deficit in slip since the great 1857 earthquake. Taking the long-term slip rate as 3.3 cm/yr, the surface slip deficit is 3±0.2 m south of Gold Hill but only 0.3±0.3 m northward from Parkfield.

Journal of Geophysical Research Solid Earth

Pb, Sr, and Nd isotopes in seamount basalts from the Juan de Fuca Ridge and Kodiak-Bowie seamount chain, northeast Pacific

Pb, Sr, and Nd isotopic ratios and their parent/daughter element concentrations for 28 basalts from 10 hotspot and nonhotspot seamounts are reported. Nd and Sr isotopic compositions ( 143 Nd/ 144 Nd = 0.51325–0.51304; 87 Sr/ 86 Sr = 0.70237–0.70275) plot in the envelope for Juan de Fuca-Gorda ridge basalts with tholeiitic basalts showing more depleted sources and a better negative correlation than transitional to alkalic basalts. Pb isotopic ratios in tholeiitic and alkalic basalts overlap ( 206 Pb/ 204 Pb = 18.29–19.44) and display a trend toward more radiogenic Pb in alkalic basalts. The isotopic data for hotspot and nonhotspot basalts are indistinguishable and correlate broadly with rock composition, implying that they are controlled by partial melting. The isotopic variation in the seamount basalts is about 60% (Nd-Sr) to 100% (Pb) of that in East Pacific Rise basalts and is interpreted as a lower limit for the magnitude of mantle heterogeneity in the northeast Pacific. The data indicate absence of a chemically distinct plume component in the linear seamount chains and strongly suggest an origin from mid-ocean ridge basalt-like east Pacific mantle.

Journal of Geophysical Research Solid Earth

Late Cretaceous paleomagnetism and clockwise rotation of the Silver Bell Mountains, south central Arizona

Late Cretaceous ash flow volcanism in the Silver Bell Mountains of southern Arizona (32.3°N, 248.5°E) was associated with caldera formation and porphyry copper mineralization. Oriented samples from 34 sites in volcanic, volcaniclastic, and intrusive units related to this episode of igneous activity (73–69 Ma) yield a mean paleomagnetic direction of I = 61.2°, D = 24.0°, α 95 = 7.6°. Primary remanent magnetizations are indicated by the presence of both normal and reversed polarities and by the significantly improved grouping (95% confidence level) of site-mean directions for a widespread tuff unit after structural correction. A comparison of this direction with a direction at Silver Bell calculated from contemporaneous units for stable North America ( I = 58.8°, D = 342.3°, α 95 = 7.7°) indicates a Late Cretaceous paleolatitude anomaly of −2.4°±7.6° (not significant at the 95% confidence level) and a declination anomaly of 41.7°±14.3° for the Silver Bell Mountains. Previously determined paleomagnetic data for southeastern Arizona suggest that this apparent clockwise rotation in the Silver Bell Mountains is a local phenomenon. Although preliminary, the average paleomagnetic direction for Oligocene and lower Miocene rocks in the Silver Bell area ( I = 43.8°, D = 357.3°, α 95 = 16.5°) is similar to that calculated for stable North America ( I = 50.2°, D = 352.2°, α 95 = 3.9°), implying that the observed rotation in the Silver Bell Mountains occurred before 26 Ma and was most likely associated with dextral strike-slip movement along the Ragged Top and related WNW trending faults bounding the Silver Bell Mountain block. These data, in conjunction with plate reconstructions and other paleomagnetic data from southwestern North America, imply that WNW trending strike-slip faults may have played an important role during Late Cretaceous to early Tertiary deformation in southern Arizona.

Journal of Geophysical Research Solid Earth

Crustal strain near the Big Bend of the San Andreas Fault: Analysis of the Los Padres-Tehachapi Trilateration Networks, California

In the region of the Los Padres-Tehachapi geodetic network, the San Andreas fault (SAF) changes its orientation by over 30° from N40°W, close to that predicted by plate motion for a transform boundary, to N73°W. The strain orientation near the SAF is consistent with right-lateral shear along the fault, with maximum shear rate of 0.38±0.01 μrad/yr at N63°W. In contrast, away from the SAF the strain orientations on both sides of the fault are consistent with the plate motion direction, with maximum shear rate of 0.19±0.01 μrad/yr at N44°W. The strain rate does not drop off rapidly away from the fault, and thus the area is fit by either a broad shear zone below the SAF or a single fault with a relatively deep locking depth. The fit to the line length data is poor for locking depth d less than 25 km. For d of 25 km a buried slip rate of 30 ± 6 mm/yr is estimated. We also estimated buried slip for models that included the Garlock and Big Pine faults, in addition to the SAF. Slip rates on other faults are poorly constrained by the Los Padres-Tehachapi network. The best fitting Garlock fault model had computed left-lateral slip rate of 11±2 mm/yr below 10 km. Buried left-lateral slip of 15±6 mm/yr on the Big Pine fault, within the Western Transverse Ranges, provides significant reduction in line length residuals; however, deformation there may be more complicated than a single vertical fault. A subhorizontal detachment on the southern side of the SAF cannot be well constrained by these data. We investigated the location of the SAF and found that a vertical fault below the surface trace fits the data much better than either a dipping fault or a fault zone located south of the surface trace.

Journal of Geophysical Research Solid Earth

Order and diversity in the modes of circum-Pacific earthquake recurrence

Recurrence characteristics of great circum-Pacific earthquakes and determinations of spatial distribution of seismic moment release are surveyed in order to delineate their general features and relate earthquake slip distribution to models of recurrent rupture. As noted by others, the pattern of moment release is typically very irregular, with strong concentrations in a few isolated regions of a much larger aftershock zone. Despite this complexity, rupture nucleation is notably systematic, with mainshock epicenter showing a strong tendency to locate in or immediately adjacent to identified regions of high moment release. This generalization suggests that earthquake recurrence is more likely to be controlled by maximum rather than average fault slip. Well-documented case histories from twelve plate boundary regions demonstrate that seismic strain release tends to be temporally well-ordered, while source dimensions, slip, and cumulative moment release vary considerably from cycle to cycle. On none of the segments studied do earthquakes consistently recur in nearly identical events. Instead, individual great earthquakes differ significantly from cycle to cycle or rerupture takes place in a sequence of two or more smaller events. Despite these differences the duration of the seismic cycle is approximately uniform. Furthermore, when strain release occurs in a sequence of large earthquakes, these events take place towards the end of the cycle and occur with increasing rupture length and magnitude. These ordered and irregular features of earthquake recurrence argue for the existence of corresponding elements on plate boundary faults. The ordered characteristics are identified with zones of concentrated moment release that slip comparable amounts in each cycle and have high shear strength. The irregular features are associated with intervening weaker regions that move in response to the stress concentration of dynamic rupture and slip by differing amounts in each event. Although these mechanistic associations are indirect and tentative, the features of recurrent behavior documented here have implications for long-term earthquake hazard assessment that are not dependent on the models proposed to explain them. Chief among these are the cycle-to-cycle differences among gap-filling events and the absence of shocks that fill major slip-deficient regions early in the seismic cycle.

Journal of Geophysical Research Solid Earth

Cooling rate and thermal structure determined from progressive magnetization of the dacite dome at Mount St. Helens, Washington

Our study of a magnetic anomaly associated with the recently active dacite dome at Mount St. Helens suggests that the dome consists of a hot, nonmagnetized core surrounded by a cool, magnetized carapace and flanking talus. The talus does not contribute to the anomaly because its constituent blocks are randomly oriented. Temporal changes in the magnetic anomaly indicate that the magnetized carapace thickened at an average rate of 0.03±0.01 m/d from 1984 to 1986. Petrographic and rock magnetic properties of dome samples indicate that the dominant process responsible for these changes is magnetization of extensively oxidized rock at progressively deeper levels within the dome as the rock cools through its blocking temperature, rather than subsequent changes in magnetization caused by further oxidation. Newly extruded material cools rapidly for a short period as heat is conducted outward in response to convective heat loss from its surface. The cooling rate gradually declines for several weeks, and thereafter the material cools at a relatively constant rate by convective heat loss from its interior along fractures that propagate inward. The rate of internal convective heat loss through fractures varies with rainfall, snowmelt, and large-scale fracturing during subsequent eruptive episodes. In accordance with a model for solidification of the 1959 lava lake at Kilauea Iki, Hawaii, we picture the dome's magnetized carapace as being a two-phase, porous, convective zone separated from the nonmagnetized core of the dome by a thin, single-phase conductive zone. As a consequence of the heat balance between the conductive and convective zones, the blocking-temperature isotherm migrates inward at a relatively constant rate. If the dome remains inactive, the time scale for its complete magnetization is estimated to be 18–36 years, a forecast which can be refined by shallow drilling into the dome and by continuing studies of its growing magnetic anomaly.

Journal of Geophysical Research Solid Earth

Reflected and mode-converted seismic waves within the shallow Aleutian subduction zone, southern Kenai Peninsula, Alaska

Pronounced secondary phases observed in local recordings of quarry shots and earthquakes on the southern Kenai Peninsula, Alaska, are identified as reflected P and S and converted S-to-P phases originating within four depth ranges: in the upper few kilometers of the Cook Inlet Tertiary basin, at midcrustal depths within the overthrust North American plate, at about 35 km depth near the top of the Wadati-Benioff seismic zone in proximity to the inferred interplate megathrust, and at about 5–10 km below the megathrust in the subducted Pacific plate. Beneath the landward margin of the accreted Chugach terrane, the mid-upper plate (MUP) discontinuity dips as steeply as 20°–30° to the west-northwest over the depth range 12–18 km. At shallower depths it flattens and possibly arches with the crest at about 10 km depth. Similar midcrustal reflectors have been imaged about 125 km to the southwest and about 350 km to the northeast along the structural trend of the Chugach terrane. The extensive reflectors may have a common origin, possible caused by the presence of underplated rocks. Relative amplitudes and polarities of the secondary phases originating at MUP depths provide few constraints on the nature of the discontinuity. It appears that the MUP discontinuity is seismically inactive and does not represent a brittle-ductile transition zone within the upper plate. The two converted S -to- P phases generated near the top of the subducted plate could indicate a low-velocity zone associated with subducted oceanic crust.

Journal of Geophysical Research Solid Earth

A tomographic glimpse of the upper mantle source of magmas of the Jemez lineament, New Mexico

The 800-km-long Jemez lineament is the most active volcanic feature in the southwestern United States. It is the southeastern tectonic boundary of the Colorado Plateau and crosses the Rio Grande rift at the Jemez Mountains. The primary volcanism of the lineament is basaltic and has occurred in the last 4.5 m.y. To infer spatial distributions of partial melt in the upper mantle source zones for the Rio Grande rift and the Jemez lineament, we investigated the lateral variations of P wave velocity in the upper mantle beneath these features. We used teleseismic P wave delays recorded at a 22-station network to perform a damped least squares, three-dimensional inversion for these lateral variations. Our technique employed velocity interpolation within a three-dimensional grid of points, rather than using blocks of constant P wave velocity. This method allows highly realistic computation of seismic ray paths as well as accurate computation of the matrix elements in our system of equations. Determinations of resolution of results were done in two independent ways, both of which gave consistent estimates of resolution. In our best resolved volume the inversion showed no significant concentration of relative low velocity for P waves beneath the Rio Grande rift. However, directly beneath the Jemez lineament there is a ∼100-km-wide, 1–2% low-velocity feature in the depth range of 50–160 km. Because of the association of the low P wave velocity with the Jemez volcanic lineament but not with the Rio Grande rift, because lowered P wave velocity can be associated with increased partial melt, and because the volume of recent volcanism at the lineament greatly exceeds that at the rift, we infer that a large magmatic source zone exists beneath the Jemez lineament but not beneath the Rio Grande rift. This implies that the volcanic potential of the Jemez lineaments continues to greatly exceed that of the Rio Grande rift. The mantle source zones for volcanics of the Jemez lineament are not overridden by, but rather track, the motion of the North American plate; this implies that these sources are within the lithospheric plate, as is clarified in the discussion. The magmatic source zones of the Jemez lineament are modeled as due to clockwise rotation of the Colorado Plateau about a pole in northeastern Colorado. This rotation caused extension of the lithosphere beneath the Jemez lineament, permitting concentration there of partially melted rock in the upper mantle.

Journal of Geophysical Research Solid Earth

Gabbroic xenoliths from the northern Gorda Ridge: Implications for magma chamber processes under slow spreading centers

Abundant gabbroic xenoliths in porphyritic pillow basalt were dredged from the northern Gorda Ridge. The host lava is a moderately fractionated, normal mid-ocean ridge basalt with a heterogeneous glass rind (Mg numbers 56-60). Other lavas in the vicinity range from near primary (Mg number 69) to fractionated (Mg number 56). On the basis of textures and mineral compositions, the xenoliths are divided into five types. The xenoliths are not cognate to the host lava, but they are genetically related. Chemistry of mineral phases in conjunction with textural features suggests that the xenoliths formed in different parts of a convecting magma chamber that underwent a period of closed system fractionation. The chamber was filled with a large proportion of crystalline mush when new, more primitive, and less dense magma was injected and mixed incompletely with the contents in the chamber, forming the hybrid host lava. -from Authors

Journal of Geophysical Research Solid Earth

Source complexity of the 1987 Whittier Narrows, California, earthquake from the inversion of strong motion records

Strong motion records for the Whittier Narrows earthquake are inverted to obtain the history of slip. Both constant rupture velocity models and variable rupture velocity models are considered. The results show a complex rupture process within a relatively small source volume, with at least four separate concentrations of slip. Two sources are associated with the hypocenter, the larger having a slip of 55–90 cm, depending on the rupture model. These sources have a radius of approximately 2–3 km and are ringed by a region of reduced slip. The aftershocks fall within this low slip annulus. Other sources with slips from 40 to 70 cm each ring the central source region and the aftershock pattern. All the sources are predominantly thrust, although some minor right-lateral strike-slip motion is seen. The overall dimensions of the Whittier earthquake from the strong motion inversions is 10 km long (along the strike) and 6 km wide (down the dip). The preferred dip is 30° and the preferred average rupture velocity is 2.5 km/s. Moment estimates range from 7.4 to 10.0×10 24 dyn cm, depending on the rupture model.

Journal of Geophysical Research Solid Earth

The relation of catastrophic flooding of Mangala Valles, Mars, to faulting of Memnonia Fossae and Tharsis volcanism

Detailed stratigraphic relations indicate two coeval periods of catastrophic flooding and Tharsis centered faulting (producing Memnonia Fossae) in the Mangala Valles region of Mars. Major sequences of lava flows of the Tharsis Montes Formation and local, lobate plains flows were erupted during and between these channeling and faulting episodes. First, Late Hesperian channel development overlapped in time the Tharsis-centered faulting that trends north 75° to 90°E. Next, Late Hesperian/Early Amazonian flooding was coeval with faulting that trends north 55° to 70°E. In some reaches, resistant lava flows filled the early channels, resulting in inverted channel topography after the later flooding swept through. Both floods likely originated from the same graben, which probably was activated during each episode of faulting. Faulting broke through groundwater barriers and tapped confined aquifers in higher regions west and east of the point of discharge. The minimum volume of water required to erode Mangala Valles (about 5×10 12 m 3 ) may have been released through two floods that drained a few percent pore volume from a relatively permeable aquifer. The peak discharges of the floods may have lasted from days to weeks. The perched water discharged from the aquifer may have been produced by hydrothermal groundwater circulation induced by Tharsis magmatism, tectonic uplift centered at Tharsis Montes, and compaction of saturated crater ejecta due to loading by lava flows.

Journal of Geophysical Research Solid Earth

Some additional observations on inclusions in the granitic rocks of the Sierra Nevada

Microgranular quartz diorite and diorite inclusions are widespread in central Sierra Nevada granitoid rocks and are almost exclusively restricted to hornblende-bearing rocks, most commonly felsic tonalites and mafic granodiorites. Strontium 87/strontium 86 values of the inclusions and host granitoids from individual plutons generally plot on single isochrons that agree closely with previously determined ages. Neodymium 143/neodymium 144 values of five inclusions ranged from 0.51229 to 0.51248 and yield ε Nd values that correlate closely with values of their respective hosts. The Nd-Sm and Rb-Sr systematics indicate that most inclusions were in isotopic equilibrium with enclosing materials at the time of formation. Silica contents of inclusions and granitoids are contiguous, but inclusions generally contain less than, and granitoids more than, 60% SiO 2 . Linear continuity of the femic oxide and H 2 O + trends relative to SiO 2 suggests many inclusions formed as concentrations of hydrous mafic minerals. Variation of other major element oxides and trace elements support this inference. The aforementioned data are generally inconsistent with a direct genetic link between the mafic inclusions and the mantle-derived basalt. However, many initial Sierra Nevada magmas were produced by mixing of mafic, mantle-derived melts and felsic crustal materials, resulting in isotopically varied tonalitic and granodioritic magmas. Most inclusions thus represent fragmented crystal accumulations of early-formed, near-liquidus minerals generated from these previously mixed magmas.

Journal of Geophysical Research Solid Earth

Paleomagnetic rotations and the Cenozoic tectonics of the Cascade Arc, Washington, Oregon, and California

Paleomagnetic results from Cenozoic (62–12 Ma) volcanic rocks of the Cascade arc and adjacent areas indicate that moderate to large clockwise rotations are an important component of the tectonic history of the arc. Two mechanisms of rotation are suggested by the regional pattern of paleomagnetic rotations. The progressive increase in rotation toward the coast in arc and forearc rocks results from distributed dextral shear, which is likely driven by oblique subduction of oceanic plates to the west. Simple shear rotation is accommodated in the upper crust by strike-slip faulting. The right-lateral Mount St. Helens seismic zone may be an active manifestation of this process. Dextral shear probably obscures a subequal contribution to arc and forearc rotation that is driven by intraarc or backarc extension. This rotation is suggested by the average southward increase in continental margin rotations into the region outboard of the Basin and Range. The southward increase in rotation parallels a change in the arc tectonic regime from largely compressional in northern Washington to extensional in Oregon. Concomitant with this change is a southward increase in the volume of eruptive rocks and the number of basaltic vents in the arc. A progressive eastward shift of the arc volcanic front with time in the rotated arc terrane is the result of the westward pivoting of the arc block in front of a zone of extension since Eocene time. Westward migration of bimodal Basin and Range volcanism since at least 16 Ma is tracking westward rotation of the frontal arc block and growth of the Basin and Range in its wake.

Journal of Geophysical Research Solid Earth

Geophysical constraints on Washington convergent margin structure

Gravity and magnetic maps of western Washington reveal the lateral structure and fabric of the Washington Coast Range, Puget Basin, and southern Washington Cascade Range. The magnetic and gravity maps show large amplitude positive anomalies associated with the shallow but largely buried section of Washington Coast Range mafic rocks which are separated by negative anomalies over deep sedimentary basins. The positive anomalies indicate that the Coast Range mafic basement extends farther east than previously thought, at least as far east as the longitude of Seattle. Linear and steep gravity and magnetic gradients indicate many unmapped, often buried faults in the Washington Coast Range Province. Magnetic highs are also associated with mapped batholiths in the Cascade arc. Several magnetic highs observed east of the Coast Range rocks and west of these batholiths may be associated with buried Tertiary plutons or ophiolites. Two-dimensional gravity and magnetic modeling constrained with geological and other geophysical data indicate that the Coast Range Province rocks are about 1 km thick at the coast, thickening to as much as 30 km near their postulated eastern edge. A maximum boundary on the average density of the upper 15–20 km of the rocks that compose the Coast Range Province of 2920 kg/m 3 was established by the modeling, suggesting a composition largely of basalt and gabbro with little interbedded sediments. Under these rocks may be mantle or a subduction complex composed of dense mafic, ultramafic, and sedimentary rocks like that proposed to underlie Vancouver Island. Previous gravity models of the Washington margin include lower densities for the proposed subduction complex than for Vancouver Island, suggesting a lower component of mafic and ultramafic rocks than the rocks underlying Vancouver Island. However, my Washington model requires that the proposed subduction complex be more dense than the trench sediments and, therefore, that material denser than sediments be incorporated within it. The absence of continental mantle and the modeled wedge shape of the Coast Range Province upper crust suggest that erosion of the bottom of the overriding plate by subduction processes may have occurred.

Washington

Stratigraphic, lithologic, and major element geochemical constraints on magmatic evolution at Lassen volcanic center, California

The Lassen volcanic center is the most recent of several long-lived volcanic centers in the southernmost Cascade Range. These centers have erupted products ranging from basaltic andesite to rhyolite and are superimposed on a background of regional basaltic to andesitic volcanism. The evolution of the Lassen volcanic center is described in three stages. Stages I and II comprise the Brokeoff volcano, and 80 km 3 andesitic stratocone, active from 600 to 400 ka. Brokeoff volcano is compositionally equivalent to the regional basaltic andesite to andesite volcanism in the Lassen region and is the result of structurally controlled focusing of the diffuse regional mafic magmatism. Stage III comprises a silicic dome field and adjacent area of hybrid andesites and has a total volume of about 100 km 3 . Volcanism during stage III was episodic and is subdivided into four sequences of lithologically and temporally distinct lavas. Stage III began at 400 ka with a rhyolitic, caldera-forming pyroclastic eruption and chemically related lavas. Additional sequences of dacite erupted between 250–200 ka and 100–0 ka. Hybrid andesites erupted adjacent to the silicic dome field between 300 and 0 ka. Porphyritic andesite and dacite with high Al 2 O 3 , low TiO 2 , medium K 2 O and FeO/MgO ratios of 1.5–2.0 are the most abundant rock types in the Lassen volcanic center. However, the single most voluminous unit is sparsely phyric rhyolite pumice. In general, the lavas of Lassen volcanic center form a single coherent trend on major element variation diagrams and in pseudo-quaternary phase space, consistant with an origin either by fractional crystallization or magma mixing. In detail, however, the lack of systematic temporal change in silica and subtly crossing trends indicate a complex origin. A variety of statistically successful fractional crystallization models can be constructed that derive Brokeoff andesites from regional magmas. An important conclusion of the modeling is that if fractional crystallization is the process responsible for generation of Brokeoff andesite, then the parent magma must be low to medium K in geochemical affinity in order to explain the variation in K 2 O. However, although major element variation can be modeled by fractional crystallization, petrographic and stratigraphic evidence indicates that magma mixing is an important but subtle process in Brokeoff lavas and suggests that lavas evolved in small independent batches. Lavas erupted during stage III, while predominantly silicic, range from 53 to 75% SiO 2 . Disequilibrium mineral assemblages in the stage III lavas indicate that they are not directly derived from Brokeoff andesite by fractional crystallization. Mixing of silicic magma with regional mafic magma and disaggregation of andesitic quenched magmatic inclusions play dominant roles in the compositional diversity of stage III lavas.

Journal of Geophysical Research Solid Earth

Comparison of ground motion from tremors and explosions in deep gold mines

Seismic body waves, from tamped chemical explosions, two with yields of 50 and one of 150 kg, were compared with corresponding data from three mining-induced tremors with a view to testing methods of discriminating between the two types of events. Detonated at depths of about 2 km, all three explosions generated P waves for which the low-frequency spectral asymptotes agree well with corresponding results scaled down from nuclear shots at the Nevada Test Site. For the two smaller explosions, recorded underground at hypocentral distances ranging from 234 to 871 m, the P wave corner frequencies, defined by the intersections of the low- and high-frequency asymptotes, are also in fine agreement with scaled-down results from the Nevada Test Site; for the 150-kg explosion, which was recorded only at the surface, it appears that 2.39 km of upward propagation caused the corner frequency to be at least a factor of 4 lower than anticipated due to attenuation. All three explosions generated S waves that appear to be a consequence of deviatoric stress release in the immediate environs of the explosions. The three tremors analyzed here generated P and S waves whose spectra and source parameters agree well with standard earthquake source models. With regard to discrimination based on P wave spectra, we conclude that for events of fixed low-frequency spectral asymptotes, the explosions typically have higher corner frequencies than tremors or earthquakes, although counterexamples certainly exist. Interestingly, the 150-kg explosion was identified as such on the basis of P and S wave polarities that are incompatible with the normally expected double-couple source model; instead, these initial motions are consistent with an explosion in conjunction with normal faulting. The body wave spectra of this explosion and those of a nearby tremor, however, were indistinguishable.

Journal of Geophysical Research Solid Earth