Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Geophysical Research, Earth Surface”

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.

439 records · Page 25Linked to original sources
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↗
USGS2012

Electrical structure of Newberry Volcano, Oregon

From the interpretation of magnetotelluric, transient electromagnetic, and Schlumberger resistivity soundings, the electrical structure of Newberry Volcano in central Oregon is found to consist of four units. From the surface downward, the geoelectrical units are (1) very resistive, young, unaltered volcanic rock, (2) a conductive layer of older volcanic material composed of altered tuffs, (3) a thick resistive layer thought to be in part intrusive rocks, and (4) a lower-crustal conductor. This model is similar to the regional geoelectrical structure found throughout the Cascade Range. Inside the caldera, the conductive second layer corresponds to the steep temperature gradient and alteration minerals observed in the USGS Newberry 2 test hole. Drill hole information on the south and north flanks of the volcano (test holes GEO N-1 and GEO N-3, respectively) indicates that outside the caldera the conductor is due to alteration minerals (primarily smectite) and not high-temperature pore fluids. On the flanks of Newberry the conductor is generally deeper than inside the caldera, and it deepens with distance from the summit. A notable exception to this pattern is seen just west of the caldera rim, where the conductive zone is shallower than at other flank locations. The volcano sits atop a rise in the resistive layer, interpreted to be due to intrusive rocks. The intrusive material has served as a heat source to produce enhanced hydrothermal alteration and, perhaps in the case of the west-flank anomaly, elevated fluid temperatures. While no public drill hole information is available to confirm this hypothesis, the west-flank anomaly appears to be a good geothermal target. In addition to the possibility that a region on the west side of the volcano could be favorable for prospecting, part of the resistive structure under the center of the volcano could be due to a vapor-dominated environment with temperatures above 300°C. In other parts of the Cascades, pervasive alteration has produced mixed layer clays and zeolites, resulting in low-resistivity anomalies. Low resistivities cannot be assumed to indicate high-temperature pore fluids. The use of electrical methods that measure resistivity as a function of excitation frequency, such as spectral induced polarization, may provide a way of obtaining information about the type and extent of alteration.

Journal of Geophysical Research Solid Earth↗
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↗
USGS2010

Coherence of Mach fronts during heterogeneous supershear earthquake rupture propagation: Simulations and comparison with observations

We study how heterogeneous rupture propagation affects the coherence of shear and Rayleigh Mach wavefronts radiated by supershear earthquakes. We address this question using numerical simulations of ruptures on a planar, vertical strike-slip fault embedded in a three-dimensional, homogeneous, linear elastic half-space. Ruptures propagate spontaneously in accordance with a linear slip-weakening friction law through both homogeneous and heterogeneous initial shear stress fields. In the 3-D homogeneous case, rupture fronts are curved owing to interactions with the free surface and the finite fault width; however, this curvature does not greatly diminish the coherence of Mach fronts relative to cases in which the rupture front is constrained to be straight, as studied by Dunham and Bhat (2008a). Introducing heterogeneity in the initial shear stress distribution causes ruptures to propagate at speeds that locally fluctuate above and below the shear wave speed. Calculations of the Fourier amplitude spectra (FAS) of ground velocity time histories corroborate the kinematic results of Bizzarri and Spudich (2008a): (1) The ground motion of a supershear rupture is richer in high frequency with respect to a subshear one. (2) When a Mach pulse is present, its high frequency content overwhelms that arising from stress heterogeneity. Present numerical experiments indicate that a Mach pulse causes approximately an ω −1.7 high frequency falloff in the FAS of ground displacement. Moreover, within the context of the employed representation of heterogeneities and over the range of parameter space that is accessible with current computational resources, our simulations suggest that while heterogeneities reduce peak ground velocity and diminish the coherence of the Mach fronts, ground motion at stations experiencing Mach pulses should be richer in high frequencies compared to stations without Mach pulses. In contrast to the foregoing theoretical results, we find no average elevation of 5%-damped absolute response spectral accelerations (SA) in the period band 0.05–0.4 s observed at stations that presumably experienced Mach pulses during the 1979 Imperial Valley, 1999 Kocaeli, and 2002 Denali Fault earthquakes compared to SA observed at non-Mach pulse stations in the same earthquakes. A 20% amplification of short period SA is seen only at a few of the Imperial Valley stations closest to the fault. This lack of elevated SA suggests that either Mach pulses in real earthquakes are even more incoherent that in our simulations or that Mach pulses are vulnerable to attenuation through nonlinear soil response. In any case, this result might imply that current engineering models of high frequency earthquake ground motions do not need to be modified by more than 20% close to the fault to account for Mach pulses, provided that the existing data are adequately representative of ground motions from supershear earthquakes.

Journal of Geophysical Research B: Solid Earth↗
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↗
USGS2012

Global petrologic variations on the moon: A ternary-diagram approach.

A ternary-diagram approach for determination of global petrologic variations on the lunar surface is presented that incorporates valuable improvements in our previous method of using geochemical variation diagrams. This new approach uses a ternary diagram that is subdivided into equally spaced segments along each of its three sides and has a triangular area in its center. The segments are assigned distinct colors that cover the range of the visible spectrum. The apexes are assigned the three primary colors, and the center triangle, which represents equal proportions of the three apexes, is assigned the color gray. The apexes are assigned the average Fe and (Th/Ti) c (ratio normalized to chondrites) compositions of KREEP/Mg-suite rocks, mare basalts, and ferroan anorthosites. The composition of each picture element (pixel) within these Apollo orbital gamma ray databases is used to determine its position within the ternary diagram. The color corresponding to this ternary position is then placed on a classification map at that pixel's position within the orbital databases. Error databases produced for the Fe and (Th/Ti) c data were used to exclude pixels having high errors. The resultant classification map shows spatial transitions among petrologic units, allows direct determination of the relative proportions of each end-member composition in a pixel, and increases the geologic interpretability of these data over that of the individual elemental databases. The classification units correspond remarkably well to observed geologic units, when we consider the spatial resolution (100 km) of the gamma ray spectrometer. Our results are as follows: (1) The highlands contain large areas of relatively pure ferroan anorthosite. (2) The average composition of the upper lunar crust is represented by an “anorthositic gabbro” composition, which supports the “magma-ocean” hypothesis for lunar crustal genesis. (3) KREEP/Mg-suite rocks are a minor fraction of the upper lunar crust. (4) Within the farside highlands, areas of KREEP/Mg-suite rocks coincide mostly with areas of crustal thinning, which are probably areas of KREEP basalt extrusions or localized excavations of Mg-suite rocks or KREEP-rich rocks. (5) Portions of the east limb and farside highlands have considerable amounts of a mafic, chondritic Th/Ti component (like mare basalt) whose occurrences coincide with mapped concentrations of light plains that display dark-halo craters; the presence of this component supports the hypothesis that mare volcanism had occurred within this highland region before the end of final heavy bombardment.

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↗