Search USGS⌕ Search

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 109 records · Page 6Linked to original sources
USGS2025

New paleomagnetic constraints on the eruption timing, stratigraphy, and post-emplacement deformation of the Picture Gorge Basalt within the Columbia River Basalt Group

The Picture Gorge Basalt (PGB) is part of the Miocene Columbia River Basalt Group (CRBG). The PGB, which outcrops in eastern Oregon, is considered coincident in time with the voluminous Grande Ronde Basalt. New radiometric ages have expanded the age‐range of the PGB, suggesting it began erupting prior to the Steens Basalt to its south and continued until after cessation of the Grande Ronde Basalt eruptions, an interval of 1.5 Ma. However, the existing paleomagnetism of the PGB implies this eruption timeline is an overestimate. To reconcile the radiometric and paleomagnetic timescales for the PGB, we conducted a paleomagnetic study on sections of the PGB to construct a detailed, high‐quality magnetostratigraphy. Our data indicate the stratigraphically lowest lava flows in the PGB are of reversed polarity, revealing a new paleomagnetic transition with the PGB and a reversed (R)–normal (N)–reversed (R) sequence. This suggests one of two timeline possibilities for PGB volcanism: (a) eruptions began and during through CRBG polarity chrons R0–N0–R1, penecontemporaneous with Steens Basalt, or (b) eruptions began and persisted during CRBG polarity chrons R1–N1–R2. Our work supports a longer interval of PGB volcanism than was suggested by previous paleomagnetic data but is at odds with the suggestion that PGB eruptions lasted through the entire main CRBG. We favor a scenario wherein PGB eruptions begin with R0 and continue into the R1 paleomagnetic interval. The paleomagnetic results also record a ∼18° vertical‐axis rotation of east‐central Oregon after ∼16 Ma with respect to stable North America.

Oregon↗
USGS2025

Rupture continuity through intermittent pauses in Cascadia slow slip events

Cascadia slow slip events (SSEs) are often envisioned as smooth, continuous ruptures, progressively activating tremor asperities as they propagate. Macroscopically, geodetic inversions and spatiotemporal maps of tremor epicenters show steady, uniform migration. In detail tremor is more chaotic and discontinuous. Larger long-term SSEs observed in daily geodetic solutions are inferred to exhibit intermittent pauses that reflect temporary re-locking of the fault, but this temporal resolution limits tests for similar re-locking on shorter timescales. We use temporal measurements of the areal growth and radiated energy of tremor clusters to investigate SSE intermittence. We find that ruptures mirror tremor pauses. Areal growth rate, however, does not reset, and removing the pauses results in smoother and more similar growth measurements among all SSEs. The rupture similarity occurs regardless of size or location and hints at an underlying uniformity and lack of predeterminism in eventual SSE size. Epicentral uncertainty precludes quantifying early rupture stages, but for larger events areal growth follows a power-law and slows with increasing size. Temporal correlations in tremor energy with inferred SSE propagation velocities and tremor rates suggest its use as a proxy for slip velocity. We find that tremor energy is tidally modulated at daily and sub-daily frequencies, and this modulation is continuous through pauses, suggesting a memory of slip state is sustained through them. We argue these pauses reflect unsteady propagation of the slip front, marked by rapid re- and un-locking, and excluding them removes rupture complexity to reveal a diffusive-like slip process and underlying universality in growth.

California, Oregon, Washington↗
USGS2026

Role of volcano-tectonic interactions during early-phase magma-assisted continental rifting

We study volcanic-tectonic interactions during the early phase of rifting in the Natron Rift, a magmatic segment of the East African Rift System. We use the open-source crustal deformation modeling code, PyLith, to model deflation during the 2007–2008 explosive eruptions and inflation during 2016–2023. We compare predictions of surface displacements resulting from magma reservoir inflation from 3D models with and without topography and homogeneous and heterogeneous material properties against Global Navigation Satellite System displacements from 2016 to 2023. We prefer the model with heterogeneous material properties and topography for fitting the observed displacements. Including topography improves the fit to GNSS data and modulates stress changes from shallow magma sources. Our numerical modeling of Coulomb stress changes suggests that deflation of a magma reservoir under the active Ol Doinyo Lengai volcano could trigger oblique normal slip on the Natron Fault with Coulomb stress increases exceeding ∼0.1 MPa on the volcano-adjacent fault segment. Fault slip on the Natron Fault is less likely during inflation of the magma reservoir because normal compressive stress increases on the fault. We also demonstrate that Coulomb stress changes are sensitive to the inclusion of topography in our model, particularly because the depth of the magma reservoir is similar to the amplitude of the topographic relief.

Ol Doinyo Lenga↗
USGS2012

Application of continuum models to deformation of the Aleutian Island Arc

Continuum models were constructed to describe large‐scale deformation of the Aleutian Island Arc over the past 5 m.y. These models consider the island arc as a continuum in the horizontal plane with the velocity boundary condition at the Pacific edge stated as a fraction of Pacific plate convergence transferred to the arc. First, a simple model of uniformly distributed strain is formulated to illustrate the mechanics of continuous deformation. Lineaments along the arc massif rotated about a vertical axis are matched by small‐element rotation calculated from the model. However, this model does not predict across‐arc variations in deformation and produces an unrealistic amount of crustal thickening after 5 m.y. A physically more meaningful model of deformation is the thin viscous sheet model based on averages of stress and rheology throughout the lithosphere. The amount of motion transferred from the Pacific plate to the arc is constrained by the rotated lineaments, while the effective stress‐strain exponent ( n ) and the ability the lithosphere has to sustain crustal thickness contrasts (the Argand number) are independent variables. Primarily, bathymetry, earthquake focal mechanisms, and styles of faulting are used to evaluate the models. The preferred model is one where the amount of motion transferred from the Pacific plate is greater in an arc‐parallel direction than in an arc‐normal direction, producing stresses consistent with strike‐slip faulting at the far western end of the arc and tensional stresses consistent with transverse normal faulting elsewhere in the arc massif. This model agrees with observations of slip vectors by Ekström and Engdahl (1989), who conclude that a portion of the arc‐parallel component of relative plate motion is taken up in the overriding plate. This model implies that compressive stress transferred to the arc is small in comparison to along‐arc shear stress and that stresses conducive to strikeslip faulting are prevalent throughout the arc.

Alaska↗
USGS2012

Compositions, growth mechanisms, and temporal relations of hydrothermal sulfide‐sulfate‐silica chimneys at the northern Cleft segment, Juan de Fuca Ridge

Three active hydrothermal vents forming sulfide mounds and chimneys (Monolith, Fountain, and Pipe Organ) and more widely distributed inactive chimneys are spatially related to a system of discontinuous fissures and young sheet flow lavas at the northern Cleft segment, Juan de Fuca Ridge. The formation of zoned tubular Curich chimneys (type I) on the Monolith sulfide mound is related to focused flow of high‐temperature (to 328°C) fluid. Bulbous chimneys (type II or “beehives”) at the Monolith and Fountain vents are products of diffuse high‐temperature (to 315°C) discharge. A broader zone of vigorous mixing between the hydrothermal fluid and seawater results in quench crystallization of anhydrite‐rich shells. Columnar Zn‐sulfide‐rich chimneys with narrow channelways (type III) are constructed where focused and relatively low‐temperature (261°C) fluid vents directly from the basalt substrate. The bulk chemistry (low Cu; high Pb, Ag, and SiO 2 contents), mineralogy (pyrite‐marcasite‐wurtzite‐amorphous silica‐anglesite), colloform and filamentous textures, and oxygen isotope characteristics of inactive (type IV) chimneys indicate a low‐temperature (<250°C) origin involving diffuse and sluggish flow patterns and conductive cooling. Seafloor observations and 210 Pb data indicate that (1) type IV chimneys are products of an earlier period of hydrothermal activity that ended no more than 60 years ago but prior to the sheet flow eruption, (2) the high‐temperature Monolith and Fountain vents are manifestations of the same heating event (shallow emplacement of magma) that led to the sheet flow eruption and recent megaplumes, and (3) the Pipe Organ Vent is in a very youthful stage of development and chimney deposition postdates the sheet flow eruption.

California, Oregon, Washington↗
USGS2012

Surface seismic and electrical methods to detect fluids related to faulting

In the absence of drilling, surface-based geophysical methods are necessary to observe fault zones and fault zone physical properties at seismogenic depths. These in situ physical properties can then be used to infer the presence and distribution of fluids along faults, although such observations are by nature indirect and become less exact with greater depth. Multiple observations of a range of such geophysical properties as compressional and shear seismic velocity ( V p and V s ), V p /V5 ratio (related to Poisson's ratio), resistivity and attenuation in and adjacent to fault zones offer the greatest hope of making inferences of the fault zone geometry, fluids in the fault zone, and fluid reservoirs in the surrounding crust. For simple geometries, fault zone guided waves can provide information on fault zone width and velocities for faults of the order of 200 m wide. To address the question of whether a narrow fault zone can be imaged well enough at depths of seismic rupture to infer the presence of anomalously high fluid/rock ratios, we present synthetic seismic tomography and magnetotelluric examples for an ideal case of a narrow fault zone with a simple geometry, large changes in material properties, and numerous earthquakes within the fault zone. A synthetic 0.5-km wide fault zone with 20% velocity reduction is well imaged using local earthquake tomography. When sequential velocity inversions are done, the true fault width is found, even to 9 km depth, although the calculated amplitude of the velocity reduction is lower than the actual amplitude. V p / V s is as well determined as V p . Magnetotelluric imaging of a synthetic fault zone shows that a conductive fault zone can be well imaged within the upper 10 km. Further, a narrow (1 km) very low resistivity (3 ohm m) fault core can be imaged within a broad (5 km) low resistivity (10 ohm m) fault zone, illustrating that regions of a fault containing large quantities of interconnected fluids within a broader, conductive fault zone should be detectable. Thus variations in fluid content and fluid pressure can be inferred from electrical and seismic methods but there will always be uncertainty in these inferences due to the trade-off with other factors, such as intrinsic variations in porosity, mineralogy, and pore geometry. The best approach is combined modeling of varied seismic and electrical data.

Journal of Geophysical Research Solid Earth↗
USGS1998

The coseismic slip distributions of the 1940 and 1979 Imperial Valley, California, earthquakes and their implications

Geodetic arrays observed by the U.S. Coast and Geodetic Survey span the Imperial fault in southern California. For the 1940 M 7.1 Imperial Valley earthquake, a 1934–1941 triangulation network has sufficient resolution to allow inversion for the coseismic slip distribution on fault segments 5 to 25 km long extending from the surface to a depth of 9 km. The estimated right-lateral slip is 0.8 to 1.7 m on the northern 30 km of the main trace of the Imperial fault, 4.8±0.2 m on a 10-km-long segment straddling the United States - Mexico border, and 1.3±0.4 m on a southern 25-km-long segment in Mexico. Fixing this strike-slip model and inverting 1940 leveling data only for dip slip yields 0.1 m of east-side-down dip slip. The seismic moment for this model is M 0 = (3.2±0.3) ×10 19 N m. The 1979 geodetic data set, mostly elevation changes from leveling routes, has insufficient resolution for inversion. However, it is possible to use this geodetic data set and results published by others to infer that the 1940 and 1979 earthquakes may be similar on the rupture zone common to both events. Our preferred 1940 model is similar to the 1979 geodetic results of Crook [1984] on the segments where both networks have good resolution. Elevation changes from 1940 and 1979 leveling data are very similar. Thus the geodetic data corroborate the surface slip evidence of Sharp [1982b] that the 1940 and 1979 slip distributions are examples of “characteristic slip” on the northern Imperial fault.

California↗
USGS1998

Suppression of large earthquakes by stress shadows: A comparison of Coulomb and rate-and-state failure

Stress shadows generated by California's two most recent great earthquakes (1857 Fort Tejon and 1906 San Francisco) substantially modified 19th and 20th century earthquake history in the Los Angeles basin and in the San Francisco Bay area. Simple Coulomb failure calculations, which assume that earthquakes can be modeled as static dislocations in an elastic half-space, have done quite well at approximating how long the stress shadows, or relaxing effects, should last and at predicting where subsequent large earthquakes will not occur. There has, however, been at least one apparent exception to the predictions of such simple models. The 1911 M >6.0 earthquake near Morgan Hill, California, occurred at a relaxed site on the Calaveras fault. We examine how the more complex rate-and-state friction formalism based on laboratory experiments might have allowed the 1911 earthquake. Rate-and-state time-to-failure calculations are consistent with the occurrence of the 1911 event just 5 years after 1906 if the Calaveras fault was already close to failure before the effects of 1906. We also examine the likelihood that the entire 78 years of relative quiet (only four M ≥6 earthquakes) in the bay area after 1906 is consistent with rate-and-state assumptions, given that the previous 7 decades produced 18 M ≥6 earthquakes. Combinations of rate-and-state variables can be found that are consistent with this pattern of large bay area earthquakes, assuming that the rate of earthquakes in the 7 decades before 1906 would have continued had 1906 not occurred. These results demonstrate that rate-and-state offers a consistent explanation for the 78-year quiescence and the 1911 anomaly, although they do not rule out several alternate explanations.

California↗
USGS1998

Strain accumulation rates in the San Francisco Bay area, 1972–1989

Maps of the strain accumulation rate in the San Francisco Bay area have been constructed from trilateration observations extending from about 1972 until the Loma Prieta earthquake in late 1989. The observations were corrected to remove offsets imposed by shallow fault creep and by four M ∼6 earthquakes that occurred in the Bay area during that time interval. The Bay area was divided into 32 contiguous polygons, and the uniform (in both space and time) strain rates that best explain the changes in the corrected (earthquake and shallow fault creep offsets removed) distances within each polygon were calculated. In a coordinate system with the 1 axis directed N58°E and the 2 axis N32°W (perpendicular and parallel to the local tangent to the small circle drawn about the Pacific-Sierra Nevada pole of rotation) the averages of these 32 strain rates (each weighted by the area of the polygon) are ε11=9.2±7.4, ε12=−160.7±4.6, and ε22=8.2±6.2 nanostrain/yr, where extension is reckoned positive and quoted uncertainties are standard deviations. As expected from the Pacific-Sierra Nevada relative plate motion, the overall strain rate is predominantly right-lateral shear across a vertical plane striking N32°W. The net increase in the 12,225 km2 area of the trilateration network is only 212±110 m2/yr, which arises from almost equal extensions in the N32°W and N58°E directions. Within the network the strain rates vary from polygon to polygon. Significant areal dilatation rates are observed in almost ⅓ of the individual polygons and the N32°W extension rates tend to be negative to the west of the Hayward-Rodgers Creek fault trend and positive east of it. The N58E extension rate is positive in 22 out of 32 polygons, a proportion that is significantly larger than would be expected by chance if the N58°E extension rate were zero or negative. The pre-1989 strain accumulation across the eventual site of the Loma Prieta rupture involves fault normal contraction as well as right-lateral shear, consistent with the rupture mechanism.

California↗
USGS2012

Paleozoic paleomagnetism and northward drift of the Alexander terrane, southeastern Alaska

Paleozoic limestone, graywacke, sandstone, milestone, red beds and volcanic rocks of the Alexander terrane, southeastern Alaska, have yielded six paleomagnetic pole positions after thermal and alternating-field demagnetization. These poles are from sample groups of late Middle Ordovician, Late Ordovician, Devonian, Late Devonian, and early and late Carboniferous age. To test various tectonic models for the structural development of this part of western North America, the paleomagnetic results are compared to those for the North American craton. It is found that the observed inclination and declination values deviate significantly from the values predicted for the present-day position of the Alexander terrane (55.5N, 133.5W). Better matching can be obtained for a paleoposition of the terrane at about 40N, 120W, in the present position of western Nevada and northeastern California. In addition, an in situ 25° clockwise rotation of the terrane is required to restore it to its original position.

Journal of Geophysical Research Solid Earth↗
USGS2012

Stress measurements at depth in the vicinity of the San Andreas fault. Implications for the magnitude of shear stress at depth

Using the hydraulic fracturing technique, we have made a systematic series of in situ stress measurements in wells drilled near the San Andreas fault. In an attempt to provide constraints for the magnitude of shear stress on the San Andreas fault at depth we have measured both the variation of stress with distance from the fault in relatively shallow (∼230 m) wells and the variation of stress with depth in a ∼1-km-deep well located 4 km from the fault. The shallow wells are located along profiles roughly perpendicular to the fault in the western Mojave desert near Palmdale and in central California where the fault is creeping. In both areas the direction of maximum compression was found to be approximately 45° from the local trend of the San Andreas. The two stress profiles show very similar results: (1) shear stress (on planes parallel to the San Andreas) increases with distance from the fault, more markedly in the western Mojave, (2) the far-field shear stress at ∼200 m depth is ∼50 bars, and (3) the horizontal principal stresses as well as shear stress increase with depth more rapidly in the wells farthest from the fault. The ∼1-km-deep well, also located in the western Mojave desert, shows increases of both horizontal principal stresses and shear stress with depth. Shear stress increases from about 25 bars at 150–300 m to about 80 bars at 750–850 m. Although this rapid increase of shear stress with depth suggests that the mean shear stress on the fault at seismogenic depths exceeds several hundred bars, the principal stresses increase with depth in a steplike manner. As this may be a near-surface effect, extrapolation of the measurements to much greater depths may not be warranted.

Journal of Geophysical Research Solid Earth↗
USGS2012

Heat flow and energetics of the San Andreas fault zone

Approximately 100 heat flow measurements in the San Andreas fault zone indicate (1) there is no evidence for local factional heating of the main fault trace at any latitude over a 1000-km length from Cape Mendocino to San Bernardino, (2) average heat flow is high (∼2 HFU, ∼80 mW m −2 ) throughout the 550-km segment of the Coast Ranges that encloses the San Andreas fault zone in central California; this broad anomaly falls off rapidly toward the Great Valley to the east, and over a 200-km distance toward the Mendocino Triple Junction to the northwest. As others have pointed out, a local conductive heat flow anomaly would be detectable unless the frictional resistance allocated to heat production on the main trace were ≲100 bars. Frictional work allocated to surface energy of new fractures is probably unimportant, and hydrologic convection is not likely to invalidate the conduction assumption, since the heat discharge by thermal springs near the fault is negligible. Explanations for the low dynamic friction fall into two intergradational classes: those in which the fault is weak all of the time and those in which it is weak only during earthquakes (possibly just large ones). The first class includes faults containing anomalously weak gouge materials and faults containing materials with normal frictional properties under near-lithostatic steady state fluid pressures. In the second class, weakening is caused by the event (for example, a thermally induced increase in fluid pressure, dehydration of clay minerals, or acoustic fluidization). In this class, unlike the first, the average strength and ambient tectonic shear stress may be large, ∼1 kbar, but the stress allocated to elastic radiation (the apparent stress) must be of similar magnitude, an apparent contradiction with seismic estimates. Unless seismic radiation is underestimated for large earthquakes, it is difficult to justify average tectonic stresses on the main trace of the San Andreas fault in excess of ∼200 bars. The development of the broad Coast Range heat flow anomaly southward from Cape Mendocino suggests that heat flow increases by a factor of 2 within 4 m.y. after the passage of the Mendocino Triple Junction. This passage leaves the San Andreas transform fault zone in its wake; the depth of the anomalous sources cannot be much greater than the depth of the seismogenic layer. Some of the anomalous heat may be supplied by conduction from the warmer mantle that must occur south of the Mendocino transform (where there is no subducting slab), and some might be supplied by shear heating in the fault zone. With no contribution from shear heating, extreme mantle upwelling would be required, and asthenosphere conditions should exist today at depths of only ∼20 km in the northernmost Coast Ranges. If there is an appreciable contribution from shear heating, the heat flow constraint implies that the seismogenic layer is partially decoupled at its base and that the basal traction is in the sense that resists right lateral motion on the fault(s). As a result of these basal tractions, the average shearing stress in the seismogenic layer would increase with distance from the main fault, and the seismogenic layer would offer substantial resistance to plate motion even though resistance on the main fault might be negligible. These speculative models have testable consequences.

Journal of Geophysical Research Solid Earth↗
USGS2012

Tectonic stresses in the lithosphere: constraints provided by the experimental deformation of rocks.

The strengths of rocks clearly place an upper limit on the stress that can be sustained by the upper half of the lithosphere. Laboratory data on rock rheology are generally lacking at intermediate temperatures and pressures on the important rock types expected in the lithosphere, so a definitive accounting of the strength distribution with depth in the upper lithosphere is presently unattainable. Analogies are drawn between the fragmentary strength data on silicates at intermediate temperature and the more extensive experimental data on marble and limestone, and several tentative conclusions are drawn: First, brittle processes, such as faulting and cataclasis, are expected to control rock strength at low pressures and temperatures. The strengths associated with these brittle mechanisms increase rapidly with increasing effective pressure and are relatively insensitive to temperature and strain rate. Second, the transitions between brittle and ductile processes occur at critical values of the least principal stress σ 3 . I suggest that the concept of the deformation mechanism map of Ashby (1972) be extended to brittle-ductile transitions by normalizing the applied differential stress σ by σ 3 , i.e., the transitions occur at critical values of σ/σ 3 . Third, the high temperature flow law of olivine, widely applied to problems involving flow in the asthenosphere, is inappropriate to the conditions of temperature found in the lithosphere, because different dislocation flow mechanisms dominate at low to intermediate temperatures. The fragmentary rheological data suggest the following rheological structure of the lithosphere where it is inelastically deforming: A rapid pressure-driven increase in rock strength with depth culminates with a shear strength maximum of up to 8 kbar at depths that depend on the state of stress and on the temperature distribution. A review of the mechanisms of weakening associated with water suggests that water weakening effects are probably not important in the oceanic lithosphere but are likely to be controlling in the continental crust.

Journal of Geophysical Research Solid Earth↗
USGS2012

Lithospheric loading by the 1896 Riku-u earthquake, northern Japan: Implications for plate flexure and asthenospheric rheology

Under favorable circumstances the time-dependent aseismic deformation resulting from the loading of the lithosphere by the stress drop of large dip slip earthquakes can be used to determine both the effective elastic plate thickness and the asthenospheric viscosity. The deformation has several similarities with the deflection of the lithosphere by surface loads and with movements due to postglacial rebound. Level changes obtained in the 80 years since the M = 7.5, 1896 Riku-u earthquake, an intraplate thrust event in northern Honshu, provide convincing evidence that asthenospheric readjustments are responsible for the observed movements. Leveling surveys crossing the zone of surface faulting have been repeated five times since 1900 and delineate a localized depression that has subsided at a continually decreasing rate. The depression is centered close to the 1896 faulting, and its shape and width, about 75 km, are matched by our model using a plate thickness of 30 km. The decaying subsidence rate constrains the viscosity of the uppermost asthenosphere to be 1×10 20 P. A linear viscous rheology matches the observed decay quite well, although measurements are sparse during the several decades following the earthquake.

Journal of Geophysical Research Solid Earth↗
USGS2012

Elastic moduli of rock glasses under pressure to 8 kilobars and geophysical implications

Shear and longitudinal velocities were measured by the ultrasonic phase comparison method as a function of pressure to 8 kbar on synthetic glasses of basalt, andesite, rhyolite, and quartz composition and on natural obsidian. Velocities of most of the glasses decrease anomalously with pressure, but increasingly more-normal behavior occurs with decrease in SiO 2 content. The pressure derivatives of rigidity and bulk modulus increase linearly, from −3.39 to −0.26 and from −5.91 to +2.09, respectively, with decrease in SiO 2 content from 100 to 49%. The change from negative to positive in the pressure derivatives of both moduli and observed at Poisson's ratio of about 0.25 is consistent with the Smyth model for the anomalous elastic behavior of glass. If the temperature in the upper mantle is about 1500°C, tholeiitic basalt would be molten in accordance with the partial melt explanation for the low-velocity zone; at 1300°C and below, basalt would be in the glassy state, especially if more felsic than tholeiite. At a temperature of 1370°C and at 30 kbar, reasonable values for the upper mantle at 100 km depth, the basalt glass of this study would have a viscosity of about 10 13 P. On the basis of the theory of viscoelasticity the glass would support shear wave propagation at frequencies above 0.1 Hz. Under such conditions of PT, 10 to 30% basalt glass in a matrix either of eclogite or olivine would reduce the seismic velocities by 3 to 9% and could also account for the values observed in the low-velocity zone.

Journal of Geophysical Research Solid Earth↗
USGS2012

Electrical properties of granite with implications for the lower crust

The electrical properties of granite appear to be dominantly controlled by the amount of free water in the granite and by temperature. Minor contributions to the electrical properties are provided by hydrostatic and lithostatic pressure, structurally bound water, oxygen fugacity, and other parameters. The effect of sulfur fugacity may be important but is experimentally unconfirmed. In addition to changing the magnitude of electrical properties, the amount and chemistry of water in granite significantly changes the temperature dependence of the electrical properties. With increasing temperature, changes in water content retain large, but lessened, effects on electrical properties. Near room temperature, a monolayer of water will decrease the electrical resistivity by an order of magnitude. Several weight-percent water may decrease the electrical resistivity by as much as 9 orders of magnitude and decrease the thermal activation energy by a factor of 5. At elevated temperatures just below granitic melting, a few weight-percent water may still decrease the resistivity by as much as 3 orders of magnitude and the activation energy by a factor of 2. Above the melting temperature (650° to 1100°C depending upon water pressure), a few weight-percent water will decrease the resistivity by less than an order of magnitude and will barely change the activation energy. Remarkably, the few weight-percent water must be present as free water. Experiments with hydrated hornblende schist (with structural water) indicate an electrical resistivity very similar to that for dry granite. The implications of these results, together with the findings of deep magnetic sounding and magnetotelluric surveys, suggest much more free water than is commonly associated with the lower crust and possibly into the upper mantle.

Journal of Geophysical Research Solid Earth↗
USGS2012

The Galapagos Spreading Centre at 86° W: A detailed geothermal field study

We report here measurements of the heat flow field of the Galapagos Spreading Center on crust of age less than 1.0 m.y. The 443 measurements in an area of about 570 km 2 reveal the general planform of the geothermal flux and permit the first truly areal estimate of the near-axis conductive heat flux. The intrusion process and associated hydrothermal circulation dominate the surface heat flow pattern, with circulation apparently continuing beyond the limits of our survey. The areal average of the conductive heat flux is 7.1 ± 0.8 HFU (295 ± 33 mW/m 2 ), about one-third the heat flux predicted by plate models. The remaining heat is apparently removed by venting of hydrothermal waters at the spreading axis and through basalt outcrops and hydrothermal mounds off axis. The pattern of surface heat flux is lineated parallel to the axis and the strongly lineated topography. Sharp lateral gradients in the heat flow, greater than 10 HFU/km near escarpments and commonly expressed as high heat flow at the tops of the scarps and lower heat flow in the valleys, may indicate a local concentration of the circulation by surface fault systems and/or variable sediment thickness.

Journal of Geophysical Research Solid Earth↗
USGS2012

The accommodation of relative motion at depth on the San Andreas fault system in California

Plate motion below the seismogenic layer along the San Andreas fault system in California is generally assumed to occur by aseismic slip along a deeper extension of the fault. It is also possible that below the seismogenic layer, deformation is distributed laterally over a zone. Several observed features of the San Andreas fault in California have implications about the mode of accommodation of relative motion along the plate boundary beneath the seismogenic zone: the shallow depth of all earthquakes in California, the depth to which coseismic slip occurred during the 1906 San Francisco earthquake, the broad zone of strain accumulation, the broad heat flow anomaly, and the existence of widely separated parallel faults. The observations strongly imply that below the seismogenic zone, relative motion is distributed over a zone and occurs by inelastic flow rather than by aseismic slip on discrete fault planes. The existence of multiple faults further suggests that tractions at the base of the brittle layer are significant over time periods of years to hundreds of years.

Journal of Geophysical Research Solid Earth↗