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 433 records · Page 24Linked to original sources
USGS2012

Rapid extension recorded by cooling‐age patterns and brittle deformation, Naxos, Greece

The metamorphic core complex exposed as the island of Naxos in the Aegean Sea records an unusually complete sequence of structures developed as a result of continental extension. The structures formed during Miocene rise and cooling from ductile, upper amphibolite facies and anatectic conditions to brittle near‐surface conditions beneath the Naxos detachment fault. Top‐to‐the‐north ductile fabrics in the footwall, which initially developed during amphibolite facies prograde metamorphism, were overprinted by a succession of north directed, normal sense lower‐temperature brittle structures as the footwall was tectonically unloaded and unroofed. Pseudotachylite and cataclasite formation, brittle faulting, alteration, and erosion of the footwall occurred during continued slip and tectonic denudation. Neogene conglomerate and megabreccia, in part derived from exhumation of this footwall, lie structurally above the peripheral Naxos fault. Published K‐Ar and 40 Ar/ 39 Ar ages for hornblende, white mica, and biotite in the footwall decrease northwestward; apparent ages 17–50 Ma in the southeast correspond to areas of low metamorphic grade where preextension argon was partially retained. Published ages 16–10 Ma in higher‐grade rocks of the domal core in the north are cooling ages that for each of the three minerals show a component of younging in the NNE direction of extension. Assuming this is the direction of unroofing, we interpret the rate of this younging as the fault slip rate as the footwall rocks moved >20 km SSW relative to their hanging wall along the base of the Naxos detachment fault. The calculated rates of slip average 5–8 mm/yr, comparable to maximum rates reported in the Basin and Range province.

Naxos↗
USGS2012

Effects of geothermal development on deformation in the Long Valley Caldera, eastern California, 1985-1994

Long Valley caldera in east central California has been the site of crustal unrest in the form of seismicity and ground deformation in response to magmatic inflation since 1980. Uplift of the resurgent dome has totaled ∼0.6 m over the 1975–1992 period. Within this region of uplift, and near the southwestern edge of the resurgent dome, is the Casa Diablo area, which experienced relative subsidence of 0.17 m between 1985 and 1992. Geothermal fluid production began at Casa Diablo in 1985 to supply a binary electric power plant; currently, three such plants generate a total of about 40 MW of electricity. The plants are supplied by 170°C water pumped from depths near 150m; all the produced water is reinjected at cooler temperatures at depths near 600 m. Analyses of data from (1) regional leveling lines and a network of bench marks in the geothermal well field area, (2) tilt observations from an L‐shaped array of bench marks, and (3) calculations relating pressure and temperature changes to subsidence indicate that relative subsidence at Casa Diablo results mainly from reductions in pressure in the shallow production reservoir and overlying formations and reductions in temperature in the underlying injection reservoir. Secondary effects on deformation are caused by offsets along the Taylor‐Bryant fault bounding the west side of the well field.

California↗
USGS1996

Identifying coseismic subsidence in tidal-wetland stratigraphic sequences at the Cascadia subduction zone of western North America

Tidal-wetland stratigraphy reveals that great plate boundary earthquakes have caused hundreds of kilometers of coast to subside at the Cascadia subduction zone. However, determining earthquake recurrence intervals and mapping the coastal extent of past great earthquake ruptures in this region are complicated by the effects of many sedimentologic, hydrographic, and oceanographic processes that occur on the coasts of tectonically passive as well as active continental margins. Tidal-wetland stratigraphy at many Cascadia estuaries differs little from that at similar sites on passive-margin coasts where stratigraphic sequences form through nonseismic processes unrelated to coseismic land level changes. Methods developed through study of similar stratigraphic sequences in Europe provide a framework for investigating the Cascadia estuarine record. Five kinds of criteria must be evaluated when inferring regional coastal subsidence due to great plate boundary earthquakes: the suddenness and amount of submergence, the lateral extent of submerged tidal-wetland soils, the coincidence of submergence with tsunami deposits, and the degree of synchroneity of submergence events at widely spaced sites. Evaluation of such criteria at the Cascadia subduction zone indicates regional coastal subsidence during at least two great earthquakes. Evidence for a coseismic origin remains equivocal, however, for the many peat-mud contacts in Cascadia stratigraphic sequences that lack (1) contrasts in lithology or fossils indicative of more than half a meter of submergence, (2) well-studied tsunami deposits, or (3) precise ages needed for regional correlation. Paleoecologic studies of fossil assemblages are particularly important in estimating the size of sudden sea level changes recorded by abrupt peat-mud contacts and in helping to distinguish erosional and gradually formed contacts from coseismic contacts. Reconstruction of a history of great earthquakes for the Cascadia subduction zone will require rigorous application of the above criteria and many detailed investigations.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Strain accumulation across the central Nevada seismic zone, 1973–1994

Five trilateration networks extending for 280 km along the central Nevada seismic zone (1915 Pleasant Valley, M = 7.3; 1954 Dixie Valley, M = 6.8; 1954 Stillwater, M = 6.8; 1954 Rainbow Mountain, M = 6.6; 1954 Fairview Peak, M = 7.1; and 1932 Cedar Mountain, M = 7.2) have been surveyed 6 times since 1973 to determine deformation along the zone. Within the precision of measurement the deformation appears uniform along the zone and is described by the principal strain rates 0.036±0.008 μstrain/yr N60°W±3° and −0.031±0.008 μstrain/yr N30°E±3°, extension reckoned positive. The observed strain rates are consistent with simple, right‐lateral, tensor shear at the rate of 0.033 μstrain/yr across a shear zone striking N15°W. This central Nevada shear zone appears to be the northward continuation of the eastern California shear zone. The orientation of the strike‐slip and normal‐slip ruptures within the central Nevada seismic zone are consistent with principal stress axes parallel to the measured principal strain rate axes. Space‐based geodetic measurements (very long baseline interferometry) indicate that the relative motion accommodated across the Basin and Range province west of Ely, Nevada, is about 9.1±1.5 mm/yr N16°W±8° (Dixon et al., 1995.) Notice that the right‐lateral shear zone postulated to explain deformation in the central Nevada seismic zone is properly oriented to accommodate that relative motion. However, a 135‐km effective width of the shear zone would be required to accommodate all of the 9.1 mm/yr relative motion at the strain rates observed in the Nevada seismic zone; only about 3 mm/yr of that relative motion is accommodated within the span of the trilateration networks.

California, Nevada↗
USGS2012

Tholeiitic‐alkalic transition at subglacial volcanoes, Tuya region, British Columbia, Canada

Ash Mountain, South Tuya, and Tuya Butte are three small basaltic volcanoes in the Stikine volcanic belt of northern British Columbia. The volcanoes rise 700, 500, and 400 m above their bases and are about 3.2, 1.6, and 2.6 km 3 in volume, respectively. They began eruptive activity under several hundred meters of overlying glacial ice, or water in an ice‐impounded lake, and undegassed pillow lava was erupted and forms the bases of all three. Later, as the vents grew into shallow water, explosive phreatomagmatic activity erupted partly degassed glassy tuffs. Finally, when the volcano emerged through the surface of the ice or water (or the water was drained), degassed subaerial lava flows were erupted and were converted to assemblages of foreset‐bedded pillow breccia and pillow lava when subaerial flows crossed a shoreline and flowed into meltwater lakes. The undegassed subglacial pillow base of Ash Mountain is overlain by partly degassed pillows and hyaloclastite tuff cut by dikes; at South Tuya the pillow base is overlain by hyaloclastite tuffs and lenses of pillow lava; at Tuya Butte the pillow base is overlain by foreset‐bedded pillow lava, pillow breccias, and hyaloclastite tuffs, which in turn are overlain by subaerial lava flows composing a small shield volcano. The undegassed basal subglacial pillow lava of the three volcanoes contain 0.10 ± 0.01 wt % sulfur and ∼0.5 wt % H 2 O. The overlying partly degassed assemblages contain 0.06 ± 0.02% sulfur and ∼0.2% H 2 O at Ash Mountain, 0.07±0.01% sulfur at South Tuya, and 0.03±0.01% sulfur at Tuya Butte. The differences in the degree of degassing can be related to the nature of eruption and quenching and the distance of flow of the subaerial lava. When the volcanoes switched from subglacial to shallow water or subaerial eruptions, as shown by change to more explosive activity and then to subaerial lava flows (and by a marked reduction of sulfur in volcanic glass), the magma shifted from tholeiitic to alkalic composition. This transition occurs at each of the three volcanoes. The tholeiitic and alkalic magmas cannot be related by shallow crystal fractionation and apparently originated by differing degrees of deep melting at a mantle source. Prior to eruption the tholeiitic melts overlay alkalic melts in shallow chambers underlying each of the volcanoes because of their lower density and were, therefore, the first to erupt under subglacial conditions. As the volcano grew through the ice (or ice‐impounded water), the volcanic conduit vented to the atmosphere, producing a partial depressurization of the conduit and the subsurface chamber. This sudden reduction in confining pressure caused enhanced vesiculation of volatile saturated melts, particularly of the more volatile‐rich alkalic melts, causing them to rise to the top of the chamber and erupt.

British Columbia↗
USGS1996

The 1992 M=7 Cape Mendocino, California, earthquake: Coseismic deformation at the south end of the Cascadia megathrust

We invert geodetic measurements of coseismic surface displacements to determine a dislocation model for the April 25, 1992, M = 7 Cape Mendocino, California, earthquake. The orientation of the model slip vector, which nearly parallels North America-Juan de Fuca relative plate convergence, and the location and orientation of the model fault relative to the offshore Cascadia megathrust, suggest that the 1992 Cape Mendocino earthquake is the first well-recorded event to relieve strain associated with the Cascadia subduction zone. We use data from three geodetic techniques: (1) the horizontal and vertical displacements of 13 monuments surveyed with the Global Positioning System, corrected for observed horizontal interseismic strain accumulation, (2) 88 section-elevation differences between leveling monuments, and (3) the uplift of 12 coastal sites observed from the die-off of intertidal marine organisms. Maximum observed displacements are 0.4 m of horizontal movement and 1.5 m of uplift along the coast. We use Monte Carlo techniques to estimate an optimal uniform slip rectangular fault geometry and its uncertainties. The optimal model using all the data resolves 4.9 m of slip on a 14 by 15 km fault that dips 28° SE. The fault extends from 1.5 to 8.7 km in depth and the main-shock hypocenter is close to the downdip projection of the fault. The shallowly dipping fault plane is consistent with the observed aftershock locations, and the estimated geodetic moment is 3.1 × 10 19 N m, 70% of the seismic moment. Other models that exclude leveling data collected in 1935 and 1942 are more consistent with seismological estimates of the fault geometry. If the earthquake is characteristic for this segment, the estimated horizontal slip vector compared with plate convergence rates suggests a recurrence interval of 140 years, with a 95% confidence range of 100–670 years. The coseismic uplift occurred in a region that also has high Quaternary uplift rates determined from marine terrace studies. If repeated ruptures of this southernmost segment of the Cascadia megathrust are responsible for the Quaternary uplift, a comparison of the coseismic uplift with coastal uplift rates suggests a recurrence interval of 200–400 years. Thus comparing horizontal and vertical coseismic to long-term deformation suggests a recurrence interval of about 100–300 years for M = 7 events at the south end of the Cascadia megathrust.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Equatorial origin for Lower Jurassic radiolarian chert in the Franciscan Complex, San Rafael Mountains, southern California

Lower Jurassic radiolarian chert sampled at two localities in the San Rafael Mountains of southern California (∼20 km north of Santa Barbara) contains four components of remanent magnetization. Components A , B ′, and B are inferred to represent uplift, Miocene volcanism, and subduction/accretion overprint magnetizations, respectively. The fourth component ( C ), isolated between 580° and 680°C, shows a magnetic polarity stratigraphy and is interpreted as a primary magnetization acquired by the chert during, or soon after, deposition. Both sequences are late Pliensbachian to middle Toarcian in age, and an average paleolatitude calculated from all tilt-corrected C components is 1° ± 3° north or south. This result is consistent with deposition of the cherts beneath the equatorial zone of high biologic productivity and is similar to initial paleolatitudes determined for chert blocks in northern California and Mexico. This result supports our model, in which deep-water Franciscan-type cherts were deposited on the Farallon plate as it moved eastward beneath the equatorial productivity high, were accreted to the continental margin at low paleolatitudes, and were subsequently distributed northward by strike-slip faulting associated with movements of the Kula, Farallon, and Pacific plates. Upper Cretaceous turbidites of the Cachuma Formation were sampled at Agua Caliente Canyon to determine a constraining paleolatitude for accretion of the Jurassic chert sequences. These apparently unaltered rocks, however, were found to be completely overprinted by the A component of magnetization. Similar in situ directions and demagnetization behaviors observed in samples of other Upper Cretaceous turbidite sequences in southern and Baja California imply that these rocks might also give unreliable results.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Recognizing and dating prehistoric liquefaction features: Lessons learned in the New Madrid seismic zone, central United States

The New Madrid seismic zone (NMSZ), which experienced severe liquefaction during the great New Madrid, Missouri, earthquakes of 1811 and 1812 as well as during several prehistoric earthquakes, is a superb laboratory for the study of world-class, earthquake-induced liquefaction features and their use in paleoseismology. In seismically active regions like the NMSZ, frequent large earthquakes can produce a complex record of liquefaction events that is difficult to interpret. Lessons learned studying liquefaction features in the NMSZ may help to unravel the paleoseismic record in other seismically active regions. Soil characteristics of liquefaction features, as well as their structural and stratigraphic relations to Native American occupation horizons and other cultural features, can help to distinguish prehistoric liquefaction features from historic features. In addition, analyses of artifact assemblages and botanical content of cultural horizons can help to narrow the age ranges of liquefaction features. Future research should focus on methods for defining source areas and estimating magnitudes of prehistoric earthquakes from liquefaction features. Also, new methods for dating liquefaction features are needed.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Thermal and unroofing history of a thick, tilted Basin-and-Range crustal section in the Tortilla Mountains, Arizona

We estimate here a geothermal gradient of only 17 ± 5°C km −1 for the tilted Grayback fault block in southeastern Arizona when extension began ∼25 Ma. This gradient is lower than preextension gradients estimated elsewhere in the Basin and Range, is only about 50% of typical gradients in the Basin and Range today, and needs to be accounted for in models of continental extension. The Grayback block exposes a 12-km-thick crustal section of Proterozoic and Cretaceous granitoids, which was tilted 90° during extension between 25 and 15 Ma. Zircon fission-track ages decrease structurally downward (westward) across the block and were all within a zone of partial track annealing prior to tilting and quenching. The zircon age gradient suggests that the 220°–240°C isotherm migrated downward 5–6 km during Paleogene erosion and regional cooling. Apatite fission-track ages decrease westward from ∼83 Ma in the structurally highest crystalline rocks to ∼24 Ma at ∼6-km paleodepth and then to ∼15 Ma another 6 km farther west. Track-length analysis confirms that apatites above the break in slope in age at ∼5.7-km paleodepth resided in a zone of partial annealing prior to tilting, and deeper apatites record rapid cooling upon tilting and unroofing beginning ∼25 Ma. At that time the 110 ± 10°C isotherm determined by the depth at which tracks in apatite were fully erased was at a basement paleodepth of ∼5.7 km, and the 220 ± 30°C isotherm as estimated from zircon data resided at a pretilting basement depth of ∼12.15 km. From consistent values of paleogeothermal gradient for two depth intervals we estimate the pretilt gradient was 17 ± 5°C km −1 . From 25 to 15 Ma the rotating Grayback block cooled rapidly as higher, westward moving blocks unroofed it tectonically at a rate of ∼1 km m.y. −1 .

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Gravity field over the Sea of Galilee: Evidence for a composite basin along a transform fault

The Sea of Galilee (Lake Kinneret) is located at the northern portion of the Kinneret-Bet Shean basin, in the northern Dead Sea transform. Three hundred kilometers of continuous marine gravity data were collected in the lake and integrated with land gravity data to a distance of more than 20 km around the lake. Analyses of the gravity data resulted in a free-air anomaly map, a variable density Bouguer anomaly map, and a horizontal first derivative map of the Bouguer anomaly. These maps, together with gravity models of profiles across the lake and the area south of it, were used to infer the geometry of the basins in this region and the main faults of the transform system. The Sea of Galilee can be divided into two units. The southern half is a pull-apart that extends to the Kinarot Valley, south of the lake, whereas the northern half was formed by rotational opening and transverse normal faults. The deepest part of the basinal area is located well south of the deepest bathymetric depression. This implies that the northeastern part of the lake, where the bathymetry is the deepest, is a young feature that is actively subsiding now. The pull-apart basin is almost symmetrical in the southern part of the lake and in the Kinarot Valley south of the lake. This suggests that the basin here is bounded by strike-slip faults on both sides. The eastern boundary fault extends to the northern part of the lake, while the western fault does not cross the northern part. The main factor controlling the structural complexity of this area is the interaction of the Dead Sea transform with a subperpendicular fault system and rotated blocks.

Sea of Galilee↗
USGS1996

Three-dimensional P and S wave velocity structure of Redoubt Volcano, Alaska

The three‐dimensional P and S wave structure of Redoubt Volcano, Alaska, and the underlying crust to depths of 7–8 km is determined from 6219 P wave and 4008 S wave first‐arrival times recorded by a 30‐station seismograph network deployed on and around the volcano. First‐arrival times are calculated using a finite‐difference technique, which allows for flexible parameterization of the slowness model and easy inclusion of topography and source‐receiver geometry. The three‐dimensional P wave velocity structure and hypocenters are determined simultaneously, while the three‐dimensional S wave velocity model is determined using the relocated seismicity and an initial S wave velocity model derived from the P wave velocity model assuming an average Vp / Vs ratio of 1.78. Convergence is steady with approximately 73% and 52% reduction in P and S wave arrival time RMS, respectively, after 10 iterations. The most prominent feature observed in the three‐dimensional velocity models derived for both P and S waves is a relative low‐velocity, near‐vertical, pipelike structure approximately 1 km in diameter that extends from 1 to 6 km beneath sea level. This feature aligns axially with the bulk of seismicity and is interpreted as a highly fractured and altered zone encompassing a magma conduit. The velocity structure beneath the north flank of the volcano between depths of 1 and 6 km is characterized by large lateral velocity variations. High velocities within this region are interpreted as remnant dikes and sills and low velocities as regions along which magma migrates. No large low‐velocity body suggestive of a magma chamber is resolved in the upper 7–8 km of the crust.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Coral ages and island subsidence, Hilo drill hole

A 25.8-m-thick sedimentary section containing coral fragments occurs directly below a surface lava flow (the ∼1340 year old Panaewa lava flow) at the Hilo drill hole. Ten coral samples from this section dated by accelerator mass spectrometry (AMS) radiocarbon and five by thermal infrared multispectral scanner (TIMS) 230 Th/U methods show good agreement. The calcareous unit is 9790 years old at the bottom and 1690 years old at the top and was deposited in a shallow lagoon behind an actively growing reef. This sedimentary unit is underlain by a 34-m-thick lava flow which in turn overlies a thin volcaniclastic silt with coral fragments that yield a single 14 C date of 10,340 years. The age-depth relations of the dated samples can be compared with proposed eustatic sea level curves after allowance for island subsidence is taken. Island subsidence averages 2.2 mm/yr for the last 47 years based on measurements from a tide gage near the drill hole or 2.5–2.6 mm/yr for the last 500,000 years based on the ages and depths of a series of drowned coral reefs offshore from west Hawaii. The age-depth measurements of coral fragments are more consistent with eustatic sea levels as determined by coral dating at Barbados and Albrolhos Islands than those based on oxygen isotopic data from deep sea cores. The Panaewa lava flow entered a lagoon underlain by coral debris and covered the drill site with 30.9 m of lava of which 11 m was above sea level. This surface has now subsided to 4.2 m above sea level, but it demonstrates how a modern lava flow entering Hilo Bay would not only change the coastline but could extensively modify the offshore shelf.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Stress/strain changes and triggered seismicity at The Geysers, California

The principal results of this study of remotely triggered seismicity in The Geysers geothermal field are the demonstration that triggering (initiation of earthquake failure) depends on a critical strain threshold and that the threshold level increases with decreasing frequency, or, equivalently, depends on strain rate. This threshold function derives from (1) analyses of dynamic strains associated with surface waves of the triggering earthquakes, (2) statistically measured aftershock zone dimensions, and (3) analytic functional representations of strains associated with power production and tides. The threshold is also consistent with triggering by static strain changes and implies that both static and dynamic strains may cause aftershocks. The observation that triggered seismicity probably occurs in addition to background activity also provides an important constraint on the triggering process. Assuming the physical processes underlying earthquake nucleation to be the same, Gomberg [this issue] discusses seismicity triggered by the M W 7.3 Landers earthquake, its constraints on the variability of triggering thresholds with site, and the implications of time delays between triggering and triggered earthquakes. Our results enable us to reject the hypothesis that dynamic strains simply nudge prestressed faults over a Coulomb failure threshold sooner than they would have otherwise. We interpret the rate-dependent triggering threshold as evidence of several competing processes with different time constants, the faster one(s) facilitating failure and the other(s) inhibiting it. Such competition is a common feature of theories of slip instability. All these results, not surprisingly, imply that to understand earthquake triggering one must consider not only simple failure criteria requiring exceedence of some constant threshold but also the requirements for generating instabilities.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Hydrogeology of the Hawaii Scientific Drilling Project borehole KP-1 1. Hydraulic conditions adjacent to the well bore

Temperature and formation resistivity logs obtained in borehole KP-1 of the Hawaii Scientific Drilling Project indicate that the adjacent formation is characterized by several zones of distinctly different average temperature and water salinity. A series of hydraulic analyses and water sampling programs were conducted to rule out the possibility of local hydraulic effects associated with the presence of the borehole in the generation of these apparent groundwater zones. Hydraulic tests and sampling with the borehole cased to a depth of 710 m and open below that depth indicate that the deep aquifer contains seawater at a temperature nearly identical to that of the open ocean at the same depth. Various analyses give estimates of aquifer transmissivity of about 10-3 m2/s in the vicinity of the borehole. Isolation of this deeper aquifer from the overlying groundwater zones was investigated by perforating the casing at six locations and then measuring the changes in water level in the borehole, in the salinity of the fluid column, in the temperature profile of the fluid column, and in the rate of flow in the fluid column induced by the perforations. These results positively confirm that the zones of distinctly different formation properties indicated on the temperature and resistivity logs are not caused by flow in or around casing. Flow and fluid column salinity induced by the perforations also confirm significant differences between the hydraulic heads and geochemistry of the different groundwater zones inferred from the well logs.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Thermal and mass implications of magmatic evolution in the Lassen volcanic region, California, and minimum constraints on basalt influx to the lower crust

We have analyzed the heat and mass demands of a petrologic model of basalt-driven magmatic evolution in which variously fractionated mafic magmas mix with silicic partial melts of the lower crust. We have formulated steady state heat budgets for two volcanically distinct areas in the Lassen region: the large, late Quaternary, intermediate to silicic Lassen volcanic center and the nearby, coeval, less evolved Caribou volcanic field. At Caribou volcanic field, heat provided by cooling and fractional crystallization of 52 km 3 of basalt is more than sufficient to produce 10 km of rhyolitic melt by partial melting of lower crust. Net heat added by basalt intrusion at Caribou volcanic field is equivalent to an increase in lower crustal heat flow of ∼7 mW m −2 , indicating that the field is not a major crustal thermal anomaly. Addition of cumulates from fractionation is offset by removal of erupted partial melts. A minimum basalt influx of 0.3 km 3 (km 2 Ma) −1 is needed to supply Caribou volcanic field. Our methodology does not fully account for an influx of basalt that remains in the crust as derivative intrusives. On the basis of comparison to deep heat flow, the input of basalt could be ∼3 to 7 times the amount we calculate. At Lassen volcanic center, at least 203 km 3 of mantle-derived basalt is needed to produce 141 km 3 of partial melt and drive the volcanic system. Partial melting mobilizes lower crustal material, augmenting the magmatic volume available for eruption at Lassen volcanic center; thus the erupted volume of 215 km 3 exceeds the calculated basalt input of 203 km 3 . The minimum basalt input of 1.6 km 3 (km 2 Ma) −1 is >5 times the minimum influx to the Caribou volcanic field. Basalt influx high enough to sustain considerable partial melting, coupled with locally high extension rate, is a crucial factor in development of Lassen volcanic center; in contrast, Caribou volcanic field has failed to develop into a large silicic center primarily because basalt supply there has been insufficient.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Rheology of water ices V and VI

We have measured the mechanical strength (σ) of pure water ices V and VI under steady state deformation conditions. Constant displacement rate compressional tests were conducted in a gas apparatus at confining pressures from 400 < P < 800 MPa, temperatures from 209 < T < 270 K, and strain rates 7 × 10 −7 < < 7 × 10 −4 s −1 . Most of the results fit to an empirical flow law of the form = A σ n exp {−( E * + PV *)/ RT }, where the four material constants A , n , E *, and V * are (for in inverse seconds and P and σ in megapascals) 10 23.0 , 6.0, 136 kJ/mol, and 29 cm 3 /mol, respectively, for ice V; and 10 6.7 , 4.5, 66 kJ/mol, and 11 cm 3 /mol, respectively, for ice VI. Ice VI may weaken to a mechanism of higher E * at T > 250 K. Ices V and VI are thus rheologically distinct but by coincidence have approximately the same strength under the conditions chosen for these experiments. To avoid misidentification, these tests are therefore accompanied by careful observations of the occurrences and characteristics of phase changes. One sample each of ice V and VI was quenched at pressure to metastably retain the high-pressure phase and the acquired deformation microstructures; X ray diffraction analysis of these samples confirmed the phase identification. Surface replicas of the deformed and quenched samples suggest that ice V probably deforms largely by dislocation creep, while ice VI deforms by a more complicated process involving substantial grain size reduction through recrystallization.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Cyclic stressing and seismicity at strongly coupled subduction zones

We use the finite element method to analyze stress variations in and near a strongly coupled subduction zone during an earthquake cycle. Deformation is assumed to be uniform along strike (plane strain on a cross section normal to the trench axis), and periodic earthquake slip is imposed consistent with the long-term rate of plate convergence and degree of coupling. Simulations of stress and displacement rate fields represent periodic fluctuations in time superimposed on an average field. The oceanic plate, descending slab, and continental lithosphere are assumed here to respond elastically to these fluctuations, and the remaining mantle under and between plates is assumed to respond as Maxwell viscoelastic. In the first part of the analysis we find that computed stress fluctuations in space and time are generally consistent with observed earthquake mechanism variations with time since a great thrust event. In particular, trench-normal extensional earthquakes tend to occur early in the earthquake cycle toward the outer rise but occur more abundantly late in the cycle in the subducting slab downdip of the main thrust zone. Compressional earthquakes, when they occur at all, have the opposite pattern. Our results suggest also that the actual timing of extensional outer rise events is controlled by the rheology of the shallow aseismic portion of the thrust interface. The second part of the analysis shows the effects of mantle relaxation on the rate of ground surface deformation during the earthquake cycle. Models without relaxation predict a strong overall compressional strain rate in the continental plate above the main thrust zone, with the strain rate constant between mainshocks. However with significant relaxation present, a localized region of unusually low compressional, or even slightly extensional, strain rate develops along the surface of the continental plate above and somewhat inland from the downdip edge of the locked main thrust zone. The low strain rate starts in the middle or late part of the cycle, depending on position. This result suggests that the negligible or small contraction measured on the Shumagin Islands, Alaska, during 1980 to 1991, may not invalidate an interpretation of that region as being a moderately coupled subduction zone. In contrast, mantle relaxation causes only modest temporal nonuniformity of uplift rates in the overriding plate and of extensional stress rates in the subducting plate, even when the Maxwell time is an order of magnitude less than the recurrence interval.

Journal of Geophysical Research B: Solid Earth↗