SEARCH · Search USGS

Results for “Journal of the Mechanics and Physics of Solids”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

25 records · Page 2Linked to original sources

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↗

USGS1996

Quantitative investigations of the Missouri gravity low: A possible expression of a large, Late Precambrian batholith intersecting the New Madrid seismic zone

Analysis of gravity and magnetic anomaly data helps characterize the geometry and physical properties of the source of the Missouri gravity low, an important cratonic feature of substantial width (about 125 km) and length (> 600 km). Filtered anomaly maps show that this prominent feature extends NW from the Reelfoot rift to the Midcontinent Rift System. Geologic reasoning and the simultaneous inversion of the gravity and magnetic data lead to an interpretation that the gravity anomaly reflects an upper crustal, 11-km-thick batholith with either near vertical or outward dipping boundaries. Considering the modeled characteristics of the batholith, structural fabric of Missouri, and relations of the batholith with plutons and regions of alteration, a tectonic model for the formation of the batholith is proposed. The model includes a mantle plume that heated the crust during Late Precambrian and melted portions of lower and middle crust, from which the low-density granitic rocks forming the batholith were partly derived. The batholith, called the Missouri batholith, may be currently related to the release of seismic energy in the New Madrid seismic zone (earthquake concentrations occur at the intersection of the Missouri batholith and the New Madrid seismic zone). Three qualitative mechanical models are suggested to explain this relationship with seismicity. Copyright 1996 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2007

Geology and geochronology of the Spirit Mountain batholith, southern Nevada: Implications for timescales and physical processes of batholith construction

The Spirit Mountain batholith (SMB) is a ∼ 250 km2 composite silicic intrusion located within the Colorado River Extensional Corridor in southernmost Nevada. Westward tilting of 40–50° has exposed a cross-section from the roof through deep levels of the batholith. Piecemeal construction is indicated by zircon geochronology, field relations, and elemental geochemistry. Zircon U/Pb data (SHRIMP) demonstrates a ∼ 2 million year (17.4–15.3 Ma) history for the SMB. Individual samples contain zircons with ages that span the lifetime of the batholith, suggesting recycling of extant zircon into new magma pulses. Field relations reveal several distinct intrusive episodes and suggest a common injection geometry of stacked horizontal sheets. The largest unit of the SMB is a gradational section (from roof downward) of high-silica leucogranite through coarse granite into foliated quartz monzonite. Solidification of this unit spans most of the history of the batholith. The 25 km × 2 km leucogranite was emplaced incrementally as subhorizontal sheets over most or all of the history of this section, suggesting repeated fractional crystallization and melt segregation events. The quartz monzonite and coarse granite are interpreted to be cumulate residuum of this fractionation. Age data from throughout this gradational unit show multiple zircon populations within individual samples. Subsequent distinct intrusions that cut this large unit, which include minor populations of zircons that record the ages of earlier events in construction of the batholith, preserve a sheeted, sill-on-sill geometry. We envision the SMB to have been a patchwork of melt-rich, melt-poor, and entirely solid zones throughout its active life. Preservation of intrusion geometries and contacts depended on the consistency of the host rock. Zircons recycled into new pulses of magma document remobilization of previously emplaced crystal mush, suggesting the mechanisms by which evidence for initial construction of the batholith became blurred.

Nevada↗

USGS2012

Laboratory triggering of stick-slip events by oscillatory loading in the presence of pore fluid with implications for physics of tectonic tremor

The physical mechanism by which the low-frequency earthquakes (LFEs) that make up portions of tectonic (also called non-volcanic) tremor are created is poorly understood. In many areas of the world, tectonic tremor and LFEs appear to be strongly tidally modulated, whereas ordinary earthquakes are not. Anomalous seismic wave speeds, interpreted as high pore fluid pressure, have been observed in regions that generate tremor. Here we build upon previous laboratory studies that investigated the response of stick-slip on artificial faults to oscillatory, tide-like loading. These previous experiments were carried out using room-dry samples of Westerly granite, at one effective stress. Here we augment these results with new experiments on Westerly granite, with the addition of varying effective stress using pore fluid at two pressures. We find that raising pore pressure, thereby lowering effective stress can significantly increase the degree of correlation of stick-slip to oscillatory loading. We also find other pore fluid effects that become important at higher frequencies, when the period of oscillation is comparable to the diffusion time of pore fluid into the fault. These results help constrain the conditions at depth that give rise to tidally modulated LFEs, providing confirmation of the effective pressure law for triggering and insights into why tremor is tidally modulated while earthquakes are at best only weakly modulated.

Journal of Geophysical Research B: Solid Earth↗

USGS2005

Applying petrophysical models to radar travel time and electrical resistivity tomograms: Resolution-dependent limitations

[1] Geophysical imaging has traditionally provided qualitative information about geologic structure; however, there is increasing interest in using petrophysical models to convert tomograms to quantitative estimates of hydrogeologic, mechanical, or geochemical parameters of interest (e.g., permeability, porosity, water content, and salinity). Unfortunately, petrophysical estimation based on tomograms is complicated by limited and variable image resolution, which depends on (1) measurement physics (e.g., electrical conduction or electromagnetic wave propagation), (2) parameterization and regularization, (3) measurement error, and (4) spatial variability. We present a framework to predict how core‐scale relations between geophysical properties and hydrologic parameters are altered by the inversion, which produces smoothly varying pixel‐scale estimates. We refer to this loss of information as “correlation loss.” Our approach upscales the core‐scale relation to the pixel scale using the model resolution matrix from the inversion, random field averaging, and spatial statistics of the geophysical property. Synthetic examples evaluate the utility of radar travel time tomography (RTT) and electrical‐resistivity tomography (ERT) for estimating water content. This work provides (1) a framework to assess tomograms for geologic parameter estimation and (2) insights into the different patterns of correlation loss for ERT and RTT. Whereas ERT generally performs better near boreholes, RTT performs better in the interwell region. Application of petrophysical models to the tomograms in our examples would yield misleading estimates of water content. Although the examples presented illustrate the problem of correlation loss in the context of near‐surface geophysical imaging, our results have clear implications for quantitative analysis of tomograms for diverse geoscience applications.

Journal of Geophysical Research B: Solid Earth↗

USGS2012

Tidal triggering of low frequency earthquakes near Parkfield, California: Implications for fault mechanics within the brittle-ductile transition

Studies of nonvolcanic tremor (NVT) have established the significant impact of small stress perturbations on NVT generation. Here we analyze the influence of the solid earth and ocean tides on a catalog of ∼550,000 low frequency earthquakes (LFEs) distributed along a 150 km section of the San Andreas Fault centered at Parkfield. LFE families are identified in the NVT data on the basis of waveform similarity and are thought to represent small, effectively co-located earthquakes occurring on brittle asperities on an otherwise aseismic fault at depths of 16 to 30 km. We calculate the sensitivity of each of these 88 LFE families to the tidally induced right-lateral shear stress (RLSS), fault-normal stress (FNS), and their time derivatives and use the hypocentral locations of each family to map the spatial variability of this sensitivity. LFE occurrence is most strongly modulated by fluctuations in shear stress, with the majority of families demonstrating a correlation with RLSS at the 99% confidence level or above. Producing the observed LFE rate modulation in response to shear stress perturbations requires low effective stress in the LFE source region. There are substantial lateral and vertical variations in tidal shear stress sensitivity, which we interpret to reflect spatial variation in source region properties, such as friction and pore fluid pressure. Additionally, we find that highly episodic, shallow LFE families are generally less correlated with tidal stresses than their deeper, continuously active counterparts. The majority of families have weaker or insignificant correlation with positive (tensile) FNS. Two groups of families demonstrate a stronger correlation with fault-normal tension to the north and with compression to the south of Parkfield. The families that correlate with fault-normal clamping coincide with a releasing right bend in the surface fault trace and the LFE locations, suggesting that the San Andreas remains localized and contiguous down to near the base of the crust. The deep families that have high sensitivity to both shear and tensile normal stress perturbations may be indicative of an increase in effective fault contact area with depth. Synthesizing our observations with those of other LFE-hosting localities will help to develop a comprehensive understanding of transient fault slip below the “seismogenic zone” by providing constraints on parameters in physical models of slow slip and LFEs.

California↗

USGS1998

Effects of slip, slip rate, and shear heating on the friction of granite

The stability of fault slip is sensitive to the way in which frictional strength responds to changes in slip rate and in particular to the effective velocity dependence of steady state friction Δμ ss /Δ ln V . This quantity can vary substantially with displacement, temperature and slip rate. To investigate the physical basis for this behavior and the possible influence of shear heating, we slid initially bare granite surfaces in unconfined rotary shear to displacements of hundreds of millimeters at normal stresses, σ n of 10 and 25 MPa and at room temperature. We imposed step changes in slip rate within the range 10 −2 to 10 3.5 μm/s and also monitored frictional heating with thermistors embedded in the granite. The transient response of μ to slip rate steps was fit to a rate- and state-dependent friction law using two state variables to estimate the values of several parameters in the constitutive law. The first 20 mm of slip shows rising friction and falling Δμ ss /Δ ln V ; further slip shows roughly constant friction, Δμ ss /Δ ln V and parameter values, suggesting that a steady state condition is reached on the fault surface. At V ≤ 10 μm/s, Δμ ss /Δ ln V = −0.004 ± 0.001. At higher rates the response is sensitive to normal stress: At σ n = 25 MPa granite shows a transition to effective velocity strengthening (Δμ ss /Δ ln V = 0.008 ± 0.004) at the highest slip rates tested. At 10 MPa granite shows a less dramatic change to Δμ ss /Δ ln V ≈ 0 at the highest rates. The maximum temperature measured in the granite is ∼60°C at 25 MPa and 10 3.5 μm/s. Temperatures are in general agreement with a numerical model of heat conduction which assumes spatially homogeneous frictional heating over the sliding surface. The simplest interpretation of our measurements of Δμ ss /Δ ln V is that the granite is inherently velocity weakening (∂μ ss /∂ ln V < 0) and temperature strengthening (∂μ ss /∂ T −1 < 0) at all velocities. At high slip rates the response of μ to changes in temperature from shear heating may outweigh the response to changing velocity, such that the net effect Δμ ss /Δ ln V > 0 mimics velocity strengthening. These results have implications for the frictional behavior of faults during earthquakes. High slip rates may cause a switch to effective velocity strengthening which could limit peak coseismic slip rate and stress drop. For fluid-saturated faults, strengthening by this mechanism may be partly or fully offset by weakening due to thermal pressurization of a poorly drained pore fluid.

Journal of Geophysical Research B: Solid Earth↗