Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Geophysical Research B: Solid Earth”

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

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

At least 217 records · Page 12Linked to original sources
USGS2003

The crustal structure from the Altai Mountains to the Altyn Tagh fault, northwest China

[1] We present a new crustal section across northwest China based on a seismic refraction profile and geologic mapping. The 1100‐km‐long section crosses the southern margin of the Chinese Altai Mountains, Junggar Accretional Belt and eastern Junggar basin, easternmost Tianshan Mountains, and easternmost Tarim basin. The crustal velocity structure and Poisson's ratio (σ), which provide a constraint on crustal composition, were determined from P and S wave data. Despite the complex geology, the crustal thickness along the entire profile is nearly uniform at 50 km. The thickest crust (56 km) occurs at the northern end of the profile beneath the Altai Mountains and the thinnest (46 km) crust is beneath the Junggar basin. Beneath surficial sediments, the crust is found to have three layers with P wave velocities ( V p ) of 6.0–6.3, 6.3–6.6, and 6.9–7.0 km/s, respectively. The southern half of the profile, including the eastern Tianshan Mountains and eastern margin of the Tarim basin, shows low P wave velocities and σ = 0.25 to a depth of 30 km, which suggests a quartz‐rich, granitic upper crustal composition. The northern half of the profile below the Altai Mountains and Junggar Accretional Belt has a higher Poisson's ratio of σ = 0.26–0.27 to a depth of 30 km, indicative of an intermediate crustal composition. The entire 1100‐km‐long profile is underlain by a 15–30 km thick high velocity (6.9–7.0 km/s; σ = 0.26–0.28) lower‐crustal layer that we interpret to have a bulk composition of mafic granulite. At the southern end of the profile, a 5‐km‐thick midcrustal low‐velocity layer ( V p = 5.9 km/s, σ = 0.25) underlies the Tianshan and the region to the south, and may be indicative of a near‐horizontal detachment interface. P n velocities are ∼7.7–7.8 km/s between the Tianshan and the Junggar basin, and ∼7.9–8.0 km/s below the Altai Mountains and eastern margin of the Tarim basin. We interpret the consistent three‐layer stratification of the crust to indicate that the crust has undergone partial melting and differentiation after Paleozoic terrane accretion. The thickness (50 km) of the crust appears to be related to compression resulting from the Indo‐Asian collision.

Altai Mountains↗
USGS2002

The 12 September 1999 Upper East Rift Zone dike intrusion at Kilauea Volcano, Hawaii

Deformation associated with an earthquake swarm on 12 September 1999 in the Upper East Rift Zone of Kilauea Volcano was recorded by continuous GPS receivers and by borehole tiltmeters. Analyses of campaign GPS, leveling data, and interferometric synthetic aperture radar (InSAR) data from the ERS-2 satellite also reveal significant deformation from the swarm. We interpret the swarm as resulting from a dike intrusion and model the deformation field using a constant pressure dike source. Nonlinear inversion was used to find the model that best fits the data. The optimal dike is located beneath and slightly to the west of Mauna Ulu, dips steeply toward the south, and strikes nearly east-west. It is approximately 3 by 2 km across and was driven by a pressure of ∼15 MPa. The total volume of the dike was 3.3 × 10 6 m 3 . Tilt data indicate a west to east propagation direction. Lack of premonitory inflation of Kilauea's summit suggests a passive intrusion; that is, the immediate cause of the intrusion was probably tensile failure in the shallow crust of the Upper East Rift Zone brought about by persistent deep rifting and by continued seaward sliding of Kilauea's south flank.

Hawaii↗
USGS2002

Strain accumulation and rotation in western Oregon and southwestern Washington

Velocities of 75 geodetic monuments in western Oregon and southwestern Washington extending from the coast to more than 300 km inland have been determined from GPS surveys over the interval 1992–2000. The average standard deviation in each of the horizontal velocity components is ∼1 mm yr −1 . The observed velocity field is approximated by a combination of rigid rotation (Euler vector relative to interior North America: 43.40°N ± 0.14°, 119.33°W ± 0.28°, and 0.822 ± 0.057° Myr −1 clockwise; quoted uncertainties are standard deviations), uniform regional strain rate (ε EE = −7.4 ± 1.8, ε EN = −3.4 ± 1.0, and ε NN = −5.0 ± 0.8 nstrain yr −1 , extension reckoned positive), and a dislocation model representing subduction of the Juan de Fuca plate beneath North America. Subduction south of 44.5°N was represented by a 40‐km‐wide locked thrust and subduction north of 44.5°N by a 75‐km‐wide locked thrust.

Oregon, Washington↗
USGS2002

High-resolution image of Calaveras Fault seismicity

By measuring relative earthquake arrival times using waveform cross correlation and locating earthquakes using the double difference technique, we are able to reduce hypocentral errors by 1 to 2 orders of magnitude over routine locations for nearly 8000 events along a 35-km section of the Calaveras Fault. This represents ∼92% of all seismicity since 1984 and includes the rupture zone of the M 6.2 1984 Morgan Hill, California, earthquake. The relocated seismicity forms highly organized structures that were previously obscured by location errors. There are abundant repeating earthquake sequences as well as linear clusters of earthquakes. Large voids in seismicity appear with dimensions of kilometers that have been aseismic over the 30-year time interval, suggesting that these portions of the fault are either locked or creeping. The area of greatest slip in the Morgan Hill main shock coincides with the most prominent of these voids, suggesting that this part of the fault may be locked between large earthquakes. We find that the Calaveras Fault at depth is extremely thin, with an average upper bound on fault zone width of 75 m. Given the location error, however, this width is not resolvably different from zero. The relocations reveal active secondary faults, which we use to solve for the stress field in the immediate vicinity of the Calaveras Fault. We find that the maximum compressive stress is at a high angle, only 13° from the fault normal, supporting previous interpretations that this fault is weak.

California↗
USGS2002

Crustal structure and relocated earthquakes in the Puget Lowland, Washington, from high-resolution seismic tomography

The availability of regional earthquake data from the Pacific Northwest Seismograph Network (PNSN), together with active source data from the Seismic Hazards Investigation in Puget Sound (SHIPS) seismic experiments, has allowed us to construct a new high-resolution 3-D, P wave velocity model of the crust to a depth of about 30 km in the central Puget Lowland. In our method, earthquake hypocenters and velocity model are jointly coupled in a fully nonlinear tomographic inversion. Active source data constrain the upper 10–15 km of the model, and earthquakes constrain the deepest portion of the model. A number of sedimentary basins are imaged, including the previously unrecognized Muckleshoot basin, and the previously incompletely defined Possession and Sequim basins. Various features of the shallow crust are imaged in detail and their structural transitions to the mid and lower crust are revealed. These include the Tacoma basin and fault zone, the Seattle basin and fault zone, the Seattle and Port Ludlow velocity highs, the Port Townsend basin, the Kingston Arch, and the Crescent basement, which is arched beneath the Lowland from its surface exposure in the eastern Olympics. Strong lateral velocity gradients, consistent with the existence of previously inferred faults, are observed, bounding the southern Port Townsend basin, the western edge of the Seattle basin beneath Dabob Bay, and portions of the Port Ludlow velocity high and the Tacoma basin. Significant velocity gradients are not observed across the southern Whidbey Island fault, the Lofall fault, or along most of the inferred location of the Hood Canal fault. Using improved earthquake locations resulting from our inversion, we determined focal mechanisms for a number of the best recorded earthquakes in the data set, revealing a complex pattern of deformation dominated by general arc-parallel regional tectonic compression. Most earthquakes occur in the basement rocks inferred to be the lower Tertiary Crescent formation. The sedimentary basins and the eastern part of the Olympic subduction complex are largely devoid of earthquakes. Clear association of hypocenters and focal mechanisms with previously mapped or proposed faults is difficult; however, seismicity, structure, and focal mechanisms associated with the Seattle fault zone suggest a possible high-angle mode of deformation with the north side up. We suggest that this deformation may be driven by isostatic readjustment of the Seattle basin.

Washington↗
USGS2002

Prehistoric ruptures of the Gurvan Bulag fault, Gobi Altay, Mongolia

The 1957 Gobi Altay M8.3 earthquake in southern Mongolia was associated with the simultaneous rupture of several faults, including the Gurvan Bulag reverse fault, which is located about 25 km south of the main strike-slip Bogd fault. Our study of paleoseismic excavations across the Gurvan Bulag fault suggests that the penultimate surface rupture occurred after 6.0 ka, most likely between 2.6 and 4.4 ka, and a possible earlier rupture occurred after 7.3 ka. Our interpretation of the stratigraphic relations in one of the exposures suggests that at least five earthquakes have generated surface rupture of the Gurvan Bulag fault since the abandonment of an ancient alluvial fan surface. Luminescence dating of sediment associated with this surface indicates that it formed either 26.6 ± 2.1 ka or 16.1 ± 2.0 ka. These data imply that the recurrence intervals for surface faulting on the Gurvan Bulag and Bogd faults are similar, on the order of several thousands of years, but that the penultimate surface ruptures of the two faults did not occur during the same earthquake.

Gobi Altay↗
USGS2002

Preeruptive inflation and surface interferometric coherence characteristics revealed by satellite radar interferometry at Makushin Volcano, Alaska: 1993-2000

Pilot reports in January 1995 and geologic field observations from the summer of 1996 indicate that a relatively small explosive eruption of Makushin, one of the more frequently active volcanoes in the Aleutian arc of Alaska, occurred on 30 January 1995. Several independent radar interferograms that each span the time period from October 1993 to September 1995 show evidence of ???7 cm of uplift centered on the volcano's east flank, which we interpret as preeruptive inflation of a ???7-km-deep magma source (??V = 0.022 km3). Subsequent interferograms for 1995-2000, a period that included no reported eruptive activity, show no evidence of additional ground deformation. Interferometric coherence at C band is found to persist for 3 years or more on lava flow and other rocky surfaces covered with short grass and sparsely distributed tall grass and for at least 1 year on most pyroclastic deposits. On lava flow and rocky surfaces with dense tall grass and on alluvium, coherence lasts for a few months. Snow and ice surfaces lose coherence within a few days. This extended timeframe of coherence over a variety of surface materials makes C band radar interferometry an effective tool for studying volcano deformation in Alaska and other similar high-latitude regions.

Journal of Geophysical Research B: Solid Earth↗
USGS2002

Undrained poroelastic response of sandstones to deviatoric stress change

Deformation of porous crustal rock through diagenesis, tectonic loading, or other processes can change pore volume and affect fluid pressure. The largest stress-induced pore pressure changes occur when fluid is trapped in pores in an “undrained” condition. We have measured the undrained poroelastic response of two sandstones to changes in mean and deviatoric stress. Pore pressure was found to respond to mean stress σ m in the usual manner: Δ p = B Δσ m ( B ranging from 0.4 to 0.7), nearly independent of the ambient deviatoric stress state. However, variations in deviatoric stress (σ d = (σ 1 − σ 3 )/2) at constant mean stress were also found to induce a reversible (elastic) pore pressure response to stress levels up to and exceeding 80% failure stress (i.e., Δ p = = const.). The coefficient η became more negative with increasing deviatoric stress level in sandstone and Ottawa sand samples. That is, η represents a dilatant response where increased deviatoric stress causes a decrease in pore pressure. The poroelastic response to deviatoric stress is explained in terms of anisotropic matrix stiffening due to closure of crack-like pore space or flattening of grain contacts at high ambient stress levels and can be important in calculations of earthquake stress transfer.

Journal of Geophysical Research B: Solid Earth↗
USGS2002

Concordant paleolatitudes for Neoproterozoic ophiolitic rocks of the Trinity Complex, Klamath Mountains, California

New paleomagnetic results from the eastern Klamath Mountains of northern California show that Neoproterozoic rocks of the Trinity ophiolitic complex and overlying Middle Devonian volcanic rocks are latitudinally concordant with cratonal North America. Combining paleomagnetic data with regional geologic and faunal evidence suggests that the Trinity Complex and related terranes of the eastern Klamath plate were linked in some fashion to the North American craton throughout that time, but that distance between them may have varied considerably. A possible model that is consistent with our paleomagnetic results and the geologic evidence is that the Trinity Complex formed and migrated parallel to paleolatitude in the basin between Laurasia and Australia–East Antarctica as the Rodinian supercontinent began to break up. It then continued to move parallel to paleolatitude at least through Middle Devonian time. Although the eastern Klamath plate served as a nucleus against which more western components of the Klamath Mountains province amalgamated, the Klamath superterrane was not accreted to North America until Early Cretaceous time.

California↗
USGS2003

Importance of a stochastic distribution of floods and erosion thresholds in the bedrock river incision problem

Fluvial erosion of bedrock occurs during occasional flood events when boundary shear stress exceeds a critical threshold to initiate incision. Therefore efforts to model the evolution of topography over long timescales should include an erosion threshold and should be driven by a stochastic distribution of erosive events. However, most bedrock incision models ignore the threshold as a second‐order detail. In addition, climate is poorly represented in most landscape evolution models, so the quantitative relationship between erosion rate and measurable climatic variables has been elusive. Here we show that the presence of an erosion threshold, when combined with a well‐constrained, probabilistic model of storm and flood occurrence, has first‐order implications for the dynamics of river incision in tectonically active areas. First, we make a direct calculation of the critical shear stress required to pluck bedrock blocks for a field site in New York. Second, we apply a recently proposed stochastic, threshold, bedrock incision model to a series of streams in California, with known tectonic and climatic forcing. Previous work in the area has identified a weak relationship between channel gradient or relief and rock uplift rate that is not easily explained by simpler detachment‐limited models. The results with the stochastic threshold model show that even low erosion thresholds, which are exceeded in steep channels during high‐frequency flood events, fundamentally affect the predicted relationship between gradient and uplift rate in steady state rivers, in a manner consistent with the observed topography. This correspondence between theory and data is, however, nonunique; models in which a thin alluvial cover may act to inhibit channel incision in the low uplift rate zone also provide plausible explanations for the observed topography. Third, we explore the broader implications of the stochastic threshold model to the development of fluvial topography in active tectonic settings. We suggest that continued field applications of geomorphic models, including physically meaningful thresholds and stochastic climate distributions, are required to advance our knowledge of interactions among surficial, climatic, and crustal processes.

California↗
USGS2002

Waveform inversion of oscillatory signatures in long-period events beneath volcanoes

The source mechanism of long-period (LP) events is examined using synthetic waveforms generated by the acoustic resonance of a fluid-filled crack. We perform a series of numerical tests in which the oscillatory signatures of synthetic LP waveforms are used to determine the source time functions of the six moment tensor components from waveform inversions assuming a point source. The results indicate that the moment tensor representation is valid for the odd modes of crack resonance with wavelengths 2 L / n , 2 W / n , n = 3, 5, 7, …, where L and W are the crack length and width, respectively. For the even modes with wavelengths 2 L / n , 2 W / n , n = 2, 4, 6, …, a generalized source representation using higher-order tensors is required, although the efficiency of seismic waves radiated by the even modes is expected to be small. We apply the moment tensor inversion to the oscillatory signatures of an LP event observed at Kusatsu-Shirane Volcano, central Japan. Our results point to the resonance of a subhorizontal crack located a few hundred meters beneath the summit crater lakes. The present approach may be useful to quantify the source location, geometry, and force system of LP events, and opens the way for moment tensor inversions of tremor.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Slip history and dynamic implications of the 1999 Chi-Chi, Taiwan, earthquake

[1] We investigate the rupture process of the 1999 Chi‐Chi, Taiwan, earthquake using extensive near‐source observations, including three‐component velocity waveforms at 36 strong motion stations and 119 GPS measurements. A three‐plane fault geometry derived from our previous inversion using only static data [ Ji et al. , 2001 ] is applied. The slip amplitude, rake angle, rupture initiation time, and risetime function are inverted simultaneously with a recently developed finite fault inverse method that combines a wavelet transform approach with a simulated annealing algorithm [ Ji et al. , 2002b ]. The inversion results are validated by the forward prediction of an independent data set, the teleseismic P and SH ground velocities, with notable agreement. The results show that the total seismic moment release of this earthquake is 2.7 × 10 20 N m and that most of the slip occurred in a triangular‐shaped asperity involving two fault segments, which is consistent with our previous static inversion. The rupture front propagates with an average rupture velocity of ∼2.0 km s −1 , and the average slip duration (risetime) is 7.2 s. Several interesting observations related to the temporal evolution of the Chi‐Chi earthquake are also investigated, including (1) the strong effect of the sinuous fault plane of the Chelungpu fault on spatial and temporal variations in slip history, (2) the intersection of fault 1 and fault 2 not being a strong impediment to the rupture propagation, and (3) the observation that the peak slip velocity near the surface is, in general, higher than on the deeper portion of the fault plane, as predicted by dynamic modeling.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Three-dimensional velocity structure of crust and upper mantle in southwestern China and its tectonic implications

Using P and S arrival times from 4625 local and regional earthquakes recorded at 174 seismic stations and associated geophysical investigations, this paper presents a three‐dimensional crustal and upper mantle velocity structure of southwestern China (21°–34°N, 97°–105°E). Southwestern China lies in the transition zone between the uplifted Tibetan plateau to the west and the Yangtze continental platform to the east. In the upper crust a positive velocity anomaly exists in the Sichuan Basin, whereas a large‐scale negative velocity anomaly exists in the western Sichuan Plateau, consistent with the upper crustal structure under the southern Tibetan plateau. The boundary between these two anomaly zones is the Longmen Shan Fault. The negative velocity anomalies at 50‐km depth in the Tengchong volcanic area and the Panxi tectonic zone appear to be associated with temperature and composition variations in the upper mantle. The Red River Fault is the boundary between the positive and negative velocity anomalies at 50‐km depth. The overall features of the crustal and the upper mantle structures in southwestern China are a low average velocity, large crustal thickness variations, the existence of a high‐conductivity layer in the crust or/and upper mantle, and a high heat flow value. All these features are closely related to the collision between the Indian and the Asian plates.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Basin-centered asperities in great subduction zone earthquakes: A link between slip, subsidence, and subduction erosion?

Published areas of high coseismic slip, or asperities, for 29 of the largest Circum-Pacific megathrust earthquakes are compared to forearc structure revealed by satellite free-air gravity, bathymetry, and seismic profiling. On average, 71% of an earthquake's seismic moment and 79% of its asperity area occur beneath the prominent gravity low outlining the deep-sea terrace; 57% of an earthquake's asperity area, on average, occurs beneath the forearc basins that lie within the deep-sea terrace. In SW Japan, slip in the 1923, 1944, 1946, and 1968 earthquakes was largely centered beneath five forearc basins whose landward edge overlies the 350??C isotherm on the plate boundary, the inferred downdip limit of the locked zone. Basin-centered coseismic slip also occurred along the Aleutian, Mexico, Peru, and Chile subduction zones but was ambiguous for the great 1964 Alaska earthquake. Beneath intrabasin structural highs, seismic slip tends to be lower, possibly due to higher temperatures and fluid pressures. Kilometers of late Cenozoic subsidence and crustal thinning above some of the source zones are indicated by seismic profiling and drilling and are thought to be caused by basal subduction erosion. The deep-sea terraces and basins may evolve not just by growth of the outer arc high but also by interseismic subsidence not recovered during earthquakes. Basin-centered asperities could indicate a link between subsidence, subduction erosion, and seismogenesis. Whatever the cause, forearc basins may be useful indicators of long-term seismic moment release. The source zone for Cascadia's 1700 A.D. earthquake contains five large, basin-centered gravity lows that may indicate potential asperities at depth. The gravity gradient marking the inferred downdip limit to large coseismic slip lies offshore, except in northwestern Washington, where the low extends landward beneath the coast. Transverse gravity highs between the basins suggest that the margin is seismically segmented and could produce a variety of large earthquakes. Published in 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Finite element model predictions of static deformation from dislocation sources in a subduction zone: Sensitivities to homogeneous, isotropic, Poisson-solid, and half-space assumptions

Dislocation models can simulate static deformation caused by slip along a fault. These models usually take the form of a dislocation embedded in a homogeneous, isotropic, Poisson-solid half-space (HIPSHS). However, the widely accepted HIPSHS assumptions poorly approximate subduction zone systems of converging oceanic and continental crust. This study uses three-dimensional finite element models (FEMs) that allow for any combination (including none) of the HIPSHS assumptions to compute synthetic Green's functions for displacement. Using the 1995 Mw = 8.0 Jalisco-Colima, Mexico, subduction zone earthquake and associated measurements from a nearby GPS array as an example, FEM-generated synthetic Green's functions are combined with standard linear inverse methods to estimate dislocation distributions along the subduction interface. Loading a forward HIPSHS model with dislocation distributions, estimated from FEMs that sequentially relax the HIPSHS assumptions, yields the sensitivity of predicted displacements to each of the HIPSHS assumptions. For the subduction zone models tested and the specific field situation considered, sensitivities to the individual Poisson-solid, isotropy, and homogeneity assumptions can be substantially greater than GPS. measurement uncertainties. Forward modeling quantifies stress coupling between the Mw = 8.0 earthquake and a nearby Mw = 6.3 earthquake that occurred 63 days later. Coulomb stress changes predicted from static HIPSHS models cannot account for the 63-day lag time between events. Alternatively, an FEM that includes a poroelastic oceanic crust, which allows for postseismic pore fluid pressure recovery, can account for the lag time. The pore fluid pressure recovery rate puts an upper limit of 10-17 m2 on the bulk permeability of the oceanic crust. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Magma supply dynamics at Westdahl volcano, Alaska, modeled from satellite radar interferometry

A group of satellite radar interferograms that span the time period from 1991 to 2000 shows that Westdahl volcano, Alaska, deflated during its 1991-1992 eruption and is reinflating at a rate that could produce another eruption within the next several years. The rates of inflation and deflation are approximated by exponential decay functions having time constants of about 6 years and a few days, respectively. This behavior is consistent with a deep, constant-pressure magma source connected to a shallow reservoir by a magma-filled conduit. An elastic deformation model indicates that the reservoir is located about 6 km below sea level and beneath Westdahl Peak. We propose that the magma flow rate through the conduit is governed by the pressure gradient between the deep source and the reservoir. The pressure gradient, and hence the flow rate, are greatest immediately after eruptions. Pressurization of the reservoir decreases both the pressure gradient and the flow rate, but eventually the reservoir ruptures and an eruption or intrusion ensues. The eruption rate is controlled partly by the pressure gradient between the reservoir and surface, and therefore it, too, decreases with time. When the supply of eruptible magma is exhausted, the eruption stops, the reservoir begins to repressurize at a high rate, and the cycle repeats. This model might also be appropriate for other frequently active volcanoes with stable magma sources and relatively simple magma storage systems.

Journal of Geophysical Research B: Solid Earth↗
USGS2004

Stress field variations in the Swiss Alps and the northern Alpine foreland derived from inversion of fault plane solutions

This study is devoted to a systematic analysis of the state of stress of the central European Alps and northern Alpine foreland in Switzerland based on focal mechanisms of 138 earthquakes with magnitudes between 1 and 5. The most robust feature of the results is that the azimuth of the minimum compressive stress, S 3 , is generally well constrained for all data subsets and always lies in the NE quadrant. However, within this quadrant, the orientation of S 3 changes systematically both along the structural strike of the Alpine chain and across it. The variation in stress along the mountain belt from NE to SW involves a progressive, counterclockwise rotation of S 3 and is most clear in the foreland, where it amounts to 45°–50°. This pattern of rotation is compatible with the disturbance to the stress field expected from the indentation of the Adriatic Block into the central European Plate, possibly together with buoyancy forces arising from the strongly arcuate structure of the Moho to the immediate west of our study area. Across the Alps, the variation in azimuth of S 3 is defined by a progressive, counterclockwise rotation of about 45° from the foreland in the north across the Helvetic domain to the Penninic nappes in the south and is accompanied by a change from a slight predominance of strike-slip mechanisms in the foreland to a strong predominance of normal faulting in the high parts of the Alps. The observed rotation can be explained by the perturbation of the large-scale regional stress by a local uniaxial deviatoric tension with a magnitude similar to that of the regional differential stress and with an orientation perpendicular to the strike of the Alpine belt. The tensile nature and orientation of this stress is consistent with the “spreading” stress expected from lateral density changes due to a crustal root beneath the Alps.

European Alps↗
USGS2004

Crystallographic controls on the frictional behavior of dry and water-saturated sheet structure minerals

We compare the frictional strengths of 17 sheet structure mineral powders, measured under dry and water-saturated conditions, to identify the factors that cause many of them to be relatively weak. The dry coefficient of friction μ ranges upward from 0.2 for graphite, leveling off at 0.8 for margarite, clintonite, gibbsite, kaolinite, and lizardite. The values of μ (dry) correlate directly with calculated (001) interlayer bond strengths of the minerals. This correlation occurs because shear becomes localized along boundary and Riedel shears and the platy minerals in them rotate into alignment with the shear planes. For those gouges with μ (dry) < 0.8, shear occurs by breaking the interlayer bonds to form new cleavage surfaces. Where μ (dry) = 0.8, consistent with Byerlee's law, the interlayer bonds are sufficiently strong that other frictional processes dominate. The transition in dry friction mechanisms corresponds to calculated surface energies of 2–3 J/m 2 . Adding water causes μ to decrease for every mineral tested except graphite. If the minerals are separated into groups with similar crystal structures, μ (wet) increases with increasing interlayer bond strength within each group. This relationship also holds for the swelling clay montmorillonite, whose water-saturated strength is consistent with the strengths of nonswelling clays of similar crystal structure. Water in the saturated gouges forms thin, structured films between the plate surfaces. The polar water molecules are bonded to the plate surfaces in proportion to the mineral's surface energy, and μ (wet) reflects the stresses required to shear through the water films.

Journal of Geophysical Research B: Solid Earth↗