USGS Science⌕ Search

SEARCH · USGS Science

Results for “Geophysical Journal”

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 1,675 records · Page 93Linked to original sources
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

The July 12, 1993, Hokkaido-Nansei-Oki, Japan, earthquake: Coseismic slip pattern from strong-motion and teleseismic recordings

We employ a finite fault inversion scheme to infer the distribution of coseismic slip for the July 12, 1993 Hokkaido-Nansei-Oki earthquake using strong ground motions recorded by the Japan Meteorological Agency within 400 km of the epicenter and vertical P waveforms recorded by the Global Digital Seismograph Network at teleseismic distances. The assumed fault geometry is based on the location of the aftershock zone and comprises two fault segments with different orientations: a northern segment striking at N20°E with a 30° dip to the west and a southern segment with an N20°W strike. For the southern segment we use both westerly and easterly dip directions to test thrust orientations previously proposed for this portion of the fault. The variance reduction is greater using a shallow west dipping segment, suggesting that the direction of dip did not change as the rupture propagated south from the hypocenter. This indicates that the earthquake resulted from the shallow underthrusting of Hokkaido beneath the Sea of Japan. Static vertical movements predicted by the corresponding distribution of fault slip are consistent with the general pattern of surface deformation observed following the earthquake. Fault rupture in the northern segment accounts for about 60% of the total P wave seismic moment of 3.4 × 10 20 N m and includes a large circular slip zone (4-m peak) near the earthquake hypocenter at depths between 10 and 25 km. Slip in the southern segment is also predominantly shallower than 25 km, but the maximum coseismic displacements (2.0–2.5 m) are observed at a depth of about 5 km. This significant shallow slip in the southern portion of the rupture zone may have been responsible for the large tsunami that devastated the small offshore island of Okushiri. Localized shallow faulting near the island, however, may require a steep westerly dip to reconcile the measured values of ground subsidence.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Mapping playa evaporite minerals and associated sediments in Death Valley, California, with multispectral thermal infrared images

Efflorescent salt crusts and associated sediments in Death Valley, California, were studied with remote-sensing data acquired by the NASA thermal infrared multispectral scanner (TIMS). Nine spectral classes that represent a variety of surface materials were distinguished, including several classes that reflect important aspects of the playa groundwater chemistry and hydrology. Evaporite crusts containing abundant thenardite (sodium sulfate) were mapped along the northern and eastern margins of the Cottonball Basin, areas where the inflow waters are rich in sodium. Gypsum (calcium sulfate) crusts were more common in the Badwater Basin, particularly near springs associated with calcic groundwaters along the western basin margin. Evaporite-rich crusts generally marked areas where groundwater is periodically near the surface and thus able to replenish the crusts though capillary evaporation. Detrital silicate minerals were prevalent in other parts of the salt pan where shallow groundwater does not affect the surface composition. The surface features in Death Valley change in response to climatic variations on several different timescales. For example, salt crusts on low-lying mudflats form and redissolve during seasonal-to-interannual cycles of wetting and desiccation. In contrast, recent flooding and erosion of rough-salt surfaces in Death Valley probably reflect increased regional precipitation spanning several decades. Remote-sensing observations of playas can provide a means for monitoring changes in evaporite facies and for better understanding the associated climatic processes. At present, such studies are limited by the availability of suitable airborne scanner data. However, with the launch of the Earth Observing System (EOS) AM-1 Platform in 1998, multispectral visible/near-infrared and thermal infrared remote-sensing data will become globally available.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Holocene paleoseismicity, temporal clustering, and probabilities of future large (M > 7) earthquakes on the Wasatch fault zone, Utah

The chronology of M > 7 paleoearthquakes on the central five segments of the Wasatch fault zone (WFZ) is one of the best dated in the world and contains 16 earthquakes in the past 5600 years with an average repeat time of 350 years. Repeat times for individual segments vary by a factor of 2, and range from about 1200 to 2600 years. Four of the central five segments ruptured between ∼620 ± 30 and 1230 ± 60 calendar years B.P. The remaining segment (Brigham City segment) has not ruptured in the past 2120 ± 100 years. Comparison of the WFZ space-time diagram of paleoearthquakes with synthetic paleoseismic histories indicates that the observed temporal clusters and gaps have about an equal probability (depending on model assumptions) of reflecting random coincidence as opposed to intersegment contagion. Regional seismicity suggests that for exposure times of 50 and 100 years, the probability for an earthquake of M > 7 anywhere within the Wasatch Front region, based on a Poisson model, is 0.16 and 0.30, respectively. A fault-specific WFZ model predicts 50 and 100 year probabilities for an M > 7 earthquake on the WFZ itself, based on a Poisson model, as 0.13 and 0.25, respectively. In contrast, segment-specific earthquake probabilities that assume quasi-periodic recurrence behavior on the Weber, Provo, and Nephi segments are less (0.01–0.07 in 100 years) than the regional or fault-specific estimates (0.25–0.30 in 100 years), due to the short elapsed times compared to average recurrence intervals on those segments. The Brigham City and Salt Lake City segments, however, have time-dependent probabilities that approach or exceed the regional and fault specific probabilities. For the Salt Lake City segment, these elevated probabilities are due to the elapsed time being approximately equal to the average late Holocene recurrence time. For the Brigham City segment, the elapsed time is significantly longer than the segment-specific late Holocene recurrence time.

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

Mauna Loa lava accumulation rates at the Hilo drill site: Formation of lava deltas during a period of declining overall volcanic growth

Accumulation rates for lava flows erupted from Mauna Loa, as sampled in the uppermost 280 m of the Hilo drill hole, vary widely for short time intervals (several thousand years), but overall are broadly similar to those documented elsewhere on this volcano since 100 ka. Thickness variations and accumulation rates for Mauna Loa lavas at the Hilo drill site have been strongly affected by local paleotopography, including funneling and ponding between Mauna Kea and Kilauea. In addition, gentle submerged slopes of Mauna Kea in Hilo Bay have permitted large shoreline displacements by Mauna Loa flows. Ages of eruptive intervals have been determined from published isotopic data and from eustatic sea level curves modified to include the isostatic subsidence of the island of Hawaii at 2.2-2.6 mm/yr. Prior to 10 ka, rates of Mauna Loa lava accumulation at the drill site varied from 0.6 to 4.3 mm/yr for dateable intervals, with an overall rate of 1.8 mm/yr. Major eruptive pulses at about 1.3 and 10 ka, each probably representing a single long-lived eruption based on lack of weathering between flow units, increase the overall accumulation rate to 2.4 mm/yr. The higher rate since 10 ka reflects construction of thick near-shoreline lava deltas as postglacial sea levels rose rapidly. Large lava deltas form only along coastal segments where initially subaerial slopes have been submerged by the combined effects of eustatic sea level rise, isostatic subsidence, or spreading of volcano flanks. Overall accumulation of 239 m of lava at the drill site since 100-120 ka closely balances submergence of the Hilo area, suggesting that processes of coastal lava deposition have been modulated by rise in sea level. The Hilo accumulation rate is slightly higher than average rates of 1-2 mm/yr determined elsewhere along the Mauna Loa coast, based on rates of shoreline coverage and dated sea cliff and fault scarp exposures. Low rates of coastal lava accumulation since 100 ka, near or below the rate of island-wide isostatic subsidence, indicate that Mauna Loa is no longer growing vigorously or even maintaining its size above sea level.

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

Stress perturbation associated with the Amazonas and other ancient continental rifts

The state of stress in the vicinity of old continental rifts is examined to investigate the possibility that crustal structure associated with ancient rifts (specifically a dense rift pillow in the lower crust) may modify substantially the regional stress field. Both shallow (2.0–2.6 km depth) breakout data and deep (20–45 km depth) crustal earthquake focal mechanisms indicate a N to NNE maximum horizontal compression in the vicinity of the Paleozoic Amazonas rift in central Brazil. This compressive stress direction is nearly perpendicular to the rift structure and represents a ∼75° rotation relative to a regional E‐W compressive stress direction in the South American plate. Elastic two‐dimensional finite element models of the density structure associated with the Amazonas rift (as inferred from independent gravity modeling) indicate that elastic support of this dense feature would generate horizontal rift‐normal compressional stresses between 60 and 120 MPa, with values of 80–100 MPa probably most representative of the overall structure. The observed ∼75° stress rotation constrains the ratio of the regional horizontal stress difference to the rift‐normal compressive stress to be between 0.25 and 1.0, suggesting that this rift‐normal stress may be from 1 to 4 times larger than the regional horizontal stress difference. A general expression for the modification of the normalized local horizontal shear stress (relative to the regional horizontal shear stress) shows that the same ratio of the rift‐normal compression relative to the regional horizontal stress difference, which controls the amount of stress rotation, also determines whether the superposed stress increases or decreases the local maximum horizontal shear stress. The potential for fault reactivation of ancient continental rifts in general is analyzed considering both the local stress rotation and modification of horizontal shear stress for both thrust and strike‐slip stress regimes. In the Amazonas rift case, because the observed stress rotation only weakly constrains the ratio of the regional horizontal stress difference to the rift‐normal compression to be between 0.25 and 1.0, our analysis is inconclusive because the resultant normalized horizontal shear stress may be reduced (for ratios >0.5) or enhanced (for ratios <0.5). Additional information is needed on all three stress magnitudes to predict how a change in horizontal shear stress directly influences the likelihood of faulting in the thrust‐faulting stress regime in the vicinity of the Amazonas rift. A rift‐normal stress associated with the seismically active New Madrid ancient rift may be sufficient to rotate the horizontal stress field consistent with strike‐slip faults parallel to the axis of the rift, although this results in a 20–40% reduction in the local horizontal shear stress within the seismic zone. Sparse stress data in the vicinity of the seismically quiescent Midcontinent rift of the central United States suggest a stress state similar to that of New Madrid, with the local horizontal shear stress potentially reduced by as much as 60%. Thus the markedly different levels of seismic activity associated with these two subparallel ancient rifts is probably due to other factors than stress perturbations due to dense rift pillows. The modeling and analysis here demonstrate that rift‐normal compressive stresses are a significant source of stress acting on the lithosphere and that in some cases may be a contributing factor to the association of intraplate seismicity with old zones of continental extension.

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↗