USGS Science⌕ Search

SEARCH · USGS Science

Results for “Geothermics”

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

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

At least 1,009 records · Page 56Linked to original sources

Heat flux from a vapor-dominated hydrothermal field beneath Yellowstone Lake

We report results from 149 heat flux measurements made over n ∼2-year interval at sites in and around a vapor-dominated geothermal field located at water depths of ∼100–120 m in Yellowstone Lake, Wyoming. Measurements of both in situ temperature and thermal conductivity as a function of depth were made with a 1 m probe via a remotely operated vehicle, and are combined to compute the vertical conductive heat flux. Inside the ∼55.5 × 10 3 m 2 bathymetric depression demarcating the vapor-dominated field, the median conductive flux is 13 W m −2 , with a conductive output of 0.72 MW. Outside the thermal field, the median conductive flux is 3.5 W m −2 . We observed 49 active vents inside the thermal field, with an estimated mass discharge rate of 56 kg s −1 , a median exit-fluid temperature of 132°C, and a total heat output of 29 MW. We find evidence for relatively weak secondary convection with a total output of 0.09 MW in thermal area lake floor sediments. Our data indicate that vapor beneath the thermal field is trapped by a low-permeability cap at a temperature of ∼189°C and a depth of ∼15 m below the lake floor. The thermal output of the Deep Hole is among the highest of any vapor-dominated field in Yellowstone, due in part to the high boiling temperatures associated with the elevated lake floor pressures.

Wyoming↗

Direct observation of the depth of active groundwater circulation in an alpine watershed

The depth of active groundwater circulation is a fundamental control on stream flows and chemistry in mountain watersheds, yet it remains challenging to characterize and is rarely well constrained. We collected hydraulic conductivity, hydraulic head, temperature, chemical, noble gas, and 3 H/ 3 He groundwater age data from discrete levels in two boreholes 46 and 81 m deep in an alpine watershed, in combination with chemical and age data from shallow groundwater discharge, to discern groundwater flow rates at different depths and directly observe active and inactive groundwater. Vertical head gradients are steep (average of 0.4) and thermal profiles are consistent with typical linear conductive continental geotherms. Groundwater deeper than ∼20 m is distinct from shallow groundwater and creek water in its chemistry, noble gas signature, and age (dominantly >65 years compared to <9 years). Together these results suggest low vertical groundwater flow velocities and a relatively shallow active circulation depth of ∼20 m. This hypothesis is tested with a simple 2‐D numerical fluid flow and heat transport model representing a hillslope transect through the two boreholes. The modeling indicates that the subhorizontally bedded sedimentary rocks underlying the basin are highly anisotropic with low vertical hydraulic conductivity, and at most ∼10% of bedrock recharge (equivalent to <2% of stream baseflow) flows below a depth of 20 m. The study demonstrates the considerable value of discrete‐depth hydrogeologic, chemical, and age data for determining active circulation depth, and illustrates an approach for maximizing the utility of individual boreholes drilled for mountain bedrock aquifer characterization.

Colorado↗

Multidisciplinary constraints on the thermal-chemical boundary between Earth's core and mantle

Heat flux from the core to the mantle provides driving energy for mantle convection thus powering plate tectonics, and contributes a significant fraction of the geothermal heat budget. Indirect estimates of core-mantle boundary heat flow are typically based on petrological evidence of mantle temperature, interpretations of temperatures indicated by seismic travel times, experimental measurements of mineral melting points, physical mantle convection models, or physical core convection models. However, previous estimates have not consistently integrated these lines of evidence. In this work, an interdisciplinary analysis is applied to co-constrain core-mantle boundary heat flow and test the thermal boundary layer (TBL) theory. The concurrence of TBL models, energy balance to support geomagnetism, seismology, and review of petrologic evidence for historic mantle temperatures supports Q CMB ∼15 TW, with all except geomagnetism supporting as high as ∼20 TW. These values provide a tighter constraint on core heat flux relative to previous work. Our work describes the seismic properties consistent with a TBL, and supports a long-lived basal mantle molten layer through much of Earth's history.

Geochemistry, Geophysics, Geosystems↗

Scaling microseismic cloud shape during hydraulic stimulation using in-situ stress and permeability

Forecasting microseismic cloud shape as a proxy of stimulated rock volume may improve the design of an energy extraction system. The microseismic cloud created during hydraulic stimulation of geothermal reservoirs is known empirically to extend in the general direction of the maximum principal stress. However, this empirical relationship is often inconsistent with reported results, and the cloud growth process remains poorly understood. This study investigates microseismic cloud growth using data obtained from a hydraulic stimulation project in Basel, Switzerland, and explores its correlation with measured in situ stress. We applied principal component analysis to a time series of microseismicity for macroscopic characterization of microseismic cloud growth in two- and three-dimensional space. The microseismic cloud, in addition to extending in the general direction of maximum principal stress, expanded in the direction of intermediate principal stress. The orientation of the least microseismic cloud growth was stable and almost identical to the minimum principal stress direction. Further, microseismic cloud shape ratios showed good agreement when compared with in situ stress magnitude ratios. The permeability tensor estimated from microseismicity also provided a good correlation in terms of direction and magnitude with the microseismic cloud growth. We show that in situ stress plays a dominant role by controlling the permeability of each existing fracture in the reservoir fracture system. Consequently, microseismic cloud growth can be scaled by in situ stress as a first-order approximation if there is sufficient variation in the orientation of existing faults.

Basel↗

Flow heterogeneity controls dissolution dynamics in topologically complex rocks

Rock dissolution is a common subsurface geochemical reaction affecting pore space properties, crucial for reservoir stimulation, carbon storage, and geothermal energy. Predictive models for dissolution remain limited due to incomplete understanding of the mechanisms involved. We examine the influence of flow, transport, and reaction regimes on mineral dissolution using 29 time-resolved data from 3D rocks. We find that initial pore structure significantly influences the dissolution pattern, with reaction rates up to two orders of magnitude lower than batch conditions, given solute and fluid-solid boundary constraints. Flow unevenness determines the location and rate of dissolution. We propose two models describing expected dissolution patterns and effective reaction rates based on dimensionless metrics for flow, transport, and reaction. Finally, we analyze feedback between evolving flow and pore structure to understand conditions that regulate/reinforce dissolution hotspots. Our findings underscore the major impact of flow arrangement on reaction-front propagation and provide a foundation for controlling dissolution hotspots.

Geophysical Research Letters↗

The United States Magnetotelluric Array and the National Impedance Map

The United States Magnetotelluric Array (USMTArray) data set, collected in the years 2006–2024, consists of more than 1,700 long-period magnetotelluric stations covering the entirety of the contiguous United States on a quasi-regular 70 km grid. Funding across multiple federal agencies was critical to sustaining this effort to its completion. Important components of the project included active guidance and participation from the MT community, the open and timely availability of all data, and the application of consistent instrumentation and robust data processing. Together with parallel advancement in the development of publicly available three-dimensional (3D) inversion codes, the USMTArray has revitalized the US magnetotelluric community and increased the visibility of magnetotellurics within the Earth-science community. Taken as a whole, these data are visualized as the National Impedance Map, which, together with a 3D synthesis conductivity model of the nation, reveals the electrical architecture of the contiguous US. USMTArray data are used by researchers worldwide for fundamental and applied studies, including investigations of continental architecture and evolution, estimation of hazards to critical infrastructure due to geomagnetic storms, and assessment of the nation's undiscovered geothermal and mineral resources. We here review the history and development of the project, discuss the challenges and successes in its execution, present the National Impedance Map and synthesis conductivity model, and highlight the breadth of research stemming from this rich data set.

conterminous United States↗

Compact seismicity bursts have different characteristics from regional seismicity

Earthquakes tend to cluster, developing into sequences driven by stress perturbations and transient fault-zone processes. Depending on the driving process, earthquake sequences show differing behaviors. This variability challenges our ability to observe or distinguish these driving processes in high resolution. Here we systematically identify seismicity bursts throughout southern California using new statistical methods and examine their causes with other independent geophysical observations. Seismicity bursts are defined as statistically significant seismicity rate anomalies. Our approach has the advantage of being data driven, depending on few earthquake occurrence assumptions. We find abundant seismicity bursts across southern California, most frequently occurring along the San Jacinto Fault and in the Salton Sea and Coso geothermal fields. These seismicity bursts are highly compact in space and time, often encompassed by a 5 km radius and have durations less than 10 hr. Many of the seismicity bursts have their largest earthquake near the beginning of the sequence, but the largest earthquake is usually not an obvious mainshock. We utilize a variety of independent geophysical data sets to analyze the characteristics of the seismicity bursts, finding that the seismicity bursts have low b-values, low spectral stress drops, and varied stress ratios compared to regional seismicity. These differences suggest that seismicity bursts are driven by transient processes acting frequently across fault networks.

California↗

Mid‐Cretaceous extensional tectonics of the Yukon‐Tanana Terrane, Trans‐Alaska Crustal Transect (TACT), east‐central Alaska

Mid‐Cretaceous crustal extension played a fundamental role in the structural evolution of the Yukon‐Tanana terrane (YTT) in the northern Cordilleran interior. In the central portion of the YTT northwest of Delta Junction, Alaska, a mylonitic shear zone juxtaposes greenschist facies rocks in the upper plate against middle to upper amphibolite facies metamorphic rocks in the lower plate, a juxtaposition suggesting elimination of as much as 10 km of crustal section. The mylonites form a partial sheath enveloping a domal footwall structure and kinematic analysis of the mylonite zone yields a uniform transport direction of hanging wall to ESE. These relations suggest analogies to the metamorphic core complexes of the southern Cordillera. However, the YTT structures are entirely ductile, suggesting either a relatively deep erosional level or relatively high geothermal gradients during extension. In the study area remnants of an older preextensional thrusting event are preserved at the highest structural levels at the base of the Seventymile terrane and the leading edge of YTT in the Wickersham terrane. However, most areas display a regional, subhorizontal fabric that is superimposed on older fabrics, and in the study area this latest fabric is subparallel to the mylonitic sheath of the apparent extensional structure. Thus the conventional viewpoint that this latest fabric is related to thrusting needs to be reevaluated and this fabric may be entirely extensional in origin. Further evidence for extension is provided by clear similarities between YTT and characteristic features of other extensional terranes. Thus we suggest that the YTT is a deeply eroded view of highly extended continental crust. The tectonic mechanism for the extensional event and the magnitude of the extension is uncertain because of complications in regional timing relationships and in alternative interpretations of the reconstruction of the crustal section. Three end‐member models based on analogies with Neogene extensional systems are presented as working models to accommodate the alternative interpretations: (1) a Jurassic collision and Cretaceous extension model based on comparisons with the Neogene history of the Mediterranean region; (2) an Early to mid‐Cretaceous syncollisional model analogous to the Carpathian Mountains of eastern Europe; and (3) a syncollisional plateau uplift model with extension driven by gravity spreading.

Alaska↗

Seismic velocity structure and composition of the continental crust: A global view

Seismic techniques provide the highest‐resolution measurements of the structure of the crust and have been conducted on a worldwide basis. We summarize the structure of the continental crust based on the results of seismic refraction profiles and infer crustal composition as a function of depth by comparing these results with high‐pressure laboratory measurements of seismic velocity for a wide range of rocks that are commonly found in the crust. The thickness and velocity structure of the crust are well correlated with tectonic province, with extended crust showing an average thickness of 30.5 km and orogens an average of 46.3 km. Shields and platforms have an average crustal thickness nearly equal to the global average. We have corrected for the nonuniform geographical distribution of seismic refraction profiles by estimating the global area of each major crustal type. The weighted average crustal thickness based on these values is 41.1 km. This value is 10% to 20% greater than previous estimates which underrepresented shields, platforms, and orogens. The average compressional wave velocity of the crust is 6.45 km/s, and the average velocity of the uppermost mantle ( Pn velocity) is 8.09 km/s. We summarize the velocity structure of the crust at 5‐km depth intervals, both in the form of histograms and as an average velocity‐depth curve, and compare these determinations with new measurements of compressional wave velocities and densities of over 3000 igneous and metamorphic rock cores made to confining pressures of 1 GPa. On the basis of petrographic studies and chemical analyses, the rocks have been classified into 29 groups. Average velocities, densities, and standard deviations are presented for each group at 5‐km depth intervals to crustal depths of 50 km along three different geotherms. This allows us to develop a model for the composition of the continental crust. Velocities in the upper continental crust are matched by velocities of a large number of lithologies, including many low‐grade metamorphic rocks and relatively silicic gneisses of amphibolite facies grade. In midcrustal regions, velocity gradients appear to originate from an increase in metamorphic grade, as well as a decrease in silica content. Tonalitic gneiss, granitic gneiss, and amphibolite are abundant midcrustal lithologies. Anisotropy due to preferred mineral orientation is likely to be significant in upper and midcrustal regions. The bulk of the lower continental crust is chemically equivalent to gabbro, with velocities in agreement with laboratory measurements of mafic granulite. Garnet becomes increasingly abundant with depth, and mafic garnet granulite is the dominant rock type immediately above the Mohorovicic discontinuity. Average compressional wave velocities of common crustal rock types show excellent correlations with density. The mean crustal density calculated from our model is 2830 kg/m 3 , and the average SiO 2 content is 61.8%.

Journal of Geophysical Research B: Solid Earth↗

Thermal and unroofing history of a thick, tilted Basin-and-Range crustal section in the Tortilla Mountains, Arizona

We estimate here a geothermal gradient of only 17 ± 5°C km −1 for the tilted Grayback fault block in southeastern Arizona when extension began ∼25 Ma. This gradient is lower than preextension gradients estimated elsewhere in the Basin and Range, is only about 50% of typical gradients in the Basin and Range today, and needs to be accounted for in models of continental extension. The Grayback block exposes a 12-km-thick crustal section of Proterozoic and Cretaceous granitoids, which was tilted 90° during extension between 25 and 15 Ma. Zircon fission-track ages decrease structurally downward (westward) across the block and were all within a zone of partial track annealing prior to tilting and quenching. The zircon age gradient suggests that the 220°–240°C isotherm migrated downward 5–6 km during Paleogene erosion and regional cooling. Apatite fission-track ages decrease westward from ∼83 Ma in the structurally highest crystalline rocks to ∼24 Ma at ∼6-km paleodepth and then to ∼15 Ma another 6 km farther west. Track-length analysis confirms that apatites above the break in slope in age at ∼5.7-km paleodepth resided in a zone of partial annealing prior to tilting, and deeper apatites record rapid cooling upon tilting and unroofing beginning ∼25 Ma. At that time the 110 ± 10°C isotherm determined by the depth at which tracks in apatite were fully erased was at a basement paleodepth of ∼5.7 km, and the 220 ± 30°C isotherm as estimated from zircon data resided at a pretilting basement depth of ∼12.15 km. From consistent values of paleogeothermal gradient for two depth intervals we estimate the pretilt gradient was 17 ± 5°C km −1 . From 25 to 15 Ma the rotating Grayback block cooled rapidly as higher, westward moving blocks unroofed it tectonically at a rate of ∼1 km m.y. −1 .

Journal of Geophysical Research B: Solid Earth↗

Stress/strain changes and triggered seismicity at The Geysers, California

The principal results of this study of remotely triggered seismicity in The Geysers geothermal field are the demonstration that triggering (initiation of earthquake failure) depends on a critical strain threshold and that the threshold level increases with decreasing frequency, or, equivalently, depends on strain rate. This threshold function derives from (1) analyses of dynamic strains associated with surface waves of the triggering earthquakes, (2) statistically measured aftershock zone dimensions, and (3) analytic functional representations of strains associated with power production and tides. The threshold is also consistent with triggering by static strain changes and implies that both static and dynamic strains may cause aftershocks. The observation that triggered seismicity probably occurs in addition to background activity also provides an important constraint on the triggering process. Assuming the physical processes underlying earthquake nucleation to be the same, Gomberg [this issue] discusses seismicity triggered by the M W 7.3 Landers earthquake, its constraints on the variability of triggering thresholds with site, and the implications of time delays between triggering and triggered earthquakes. Our results enable us to reject the hypothesis that dynamic strains simply nudge prestressed faults over a Coulomb failure threshold sooner than they would have otherwise. We interpret the rate-dependent triggering threshold as evidence of several competing processes with different time constants, the faster one(s) facilitating failure and the other(s) inhibiting it. Such competition is a common feature of theories of slip instability. All these results, not surprisingly, imply that to understand earthquake triggering one must consider not only simple failure criteria requiring exceedence of some constant threshold but also the requirements for generating instabilities.

Journal of Geophysical Research B: Solid Earth↗

Crustal shear-wave splitting from local earthquakes in the Hengill triple junction, southwest Iceland

The Hengill region in SW Iceland is an unstable ridge-ridge-transform triple junction between an active and a waning segment of the mid-Atlantic spreading center and a transform that is transgressing southward. The triple junction contains active and extinct spreading segments and a widespread geothermal area. We evaluated shear-wave birefringence for locally recorded upper-crustal earthquakes using an array of 30 three-component digital seismographs. Fast-polarization directions, ϕ, are mostly NE to NNE, subparallel to the spreading axis and probably caused by fissures and microcracks related to spreading. However, there is significant variability in ϕ throughout the array. The lag from fast to slow S is not proportional to earthquake depth (ray length), being scattered at all depths. The average wave-speed difference between qS1 and qS2 in the upper 2–5 km of the crust is 2–5%. Our results suggest considerable heterogeneity or strong S scattering.

Geophysical Research Letters↗

Ductile shear zones beneath strike-slip faults: Implications for the thermomechanics of the San Andreas fault zone

We have carried out two-dimensional (2-D) numerical experiments on the bulk flow of a layer of fluid that is driven in a strike-slip sense by constant velocities applied at its boundaries. The fluid has the (linearized) conventional rheology assumed to apply to lower crust/upper mantle rocks. The temperature dependence of the effective viscosity of the fluid and the shear heating that accompanies deformation have been incorporated into the calculations, as has thermal conduction in an overlying crustal layer. Two end-member boundary conditions have been considered, corresponding to a strong upper crust driving a weaker ductile substrate and a strong ductile layer driving a passive, weak crust. In many cases of practical interest, shear heating is concentrated close to the axial plane of the shear zone for either boundary condition. For these cases, the resulting steady state temperature field is well approximated by a cylindrical heat source embedded in a conductive half-space at a depth corresponding to the top of the fluid layer. This approximation, along with the application of a theoretical result for one-dimensional shear zones, permits us to obtain simple analytical approximations to the thermal effects of 2-D ductile shear zones for a range of assumed rheologies and crustal geotherms, making complex numerical calculations unnecessary. Results are compared with observable effects on heat flux near the San Andreas fault using constraints on the slip distribution across the entire fault system. Ductile shearing in the lower crust or upper mantle can explain the observed increase in surface heat flux southeast of the Mendocino triple junction and match the amplitude of the regional heat flux anomaly in the California Coast Ranges. Because ductile dissipation depends only weakly on slip rate, faults moving only a few millimeters per year can be important heat sources, and the superposition of effects of localized ductile shearing on both currently active and now inactive strands of the San Andreas system can explain the breadth of the heat flux anomaly across central California.

Journal of Geophysical Research B: Solid Earth↗

Creep of water ices at planetary conditions: A compilation

Many constitutive laws for the flow of ice have been published since the advent of the Voyager explorations of the outer solar system. Conflicting data have occasionally come from different laboratories, and refinement of experimental techniques has led to the publication of laws that supersede earlier ones. In addition, there are unpublished data from ongoing research that also amend the constitutive laws. Here we compile the most current laboratory-derived flow laws for water ice phases I, II, III, V, and VI, and ice I mixtures with hard particulates. The rheology of interest is mainly that of steady state, and the conditions reviewed are the pressures and temperatures applicable to the surfaces and interiors of icy moons of the outer solar system. Advances in grain-size-dependent creep in ices I and II as well as in phase transformations and metastability under differential stress are also included in this compilation. At laboratory strain rates the several ice polymorphs are rheologically distinct in terms of their stress, temperature, and pressure dependencies but, with the exception of ice III, have fairly similar strengths. Hard particulates strengthen ice I significantly only at high particulate volume fractions. Ice III has the potential for significantly affecting mantle dynamics because it is much weaker than the other polymorphs and its region of stability, which may extend metastably well into what is nominally the ice II field, is located near likely geotherms of large icy moons.

Journal of Geophysical Research E: Planets↗

Non-double-couple earthquakes. 1. Theory

Historically, most quantitative seismological analyses have been based on the assumption that earthquakes are caused by shear faulting, for which the equivalent force system in an isotropic medium is a pair of force couples with no net torque (a “double couple,” or DC). Observations of increasing quality and coverage, however, now resolve departures from the DC model for many earthquakes and find some earthquakes, especially in volcanic and geothermal areas, that have strongly non-DC mechanisms. Understanding non-DC earthquakes is important both for studying the process of faulting in detail and for identifying nonshear-faulting processes that apparently occur in some earthquakes. This paper summarizes the theory of “moment tensor” expansions of equivalent-force systems and analyzes many possible physical non-DC earthquake processes. Contrary to long-standing assumption, sources within the Earth can sometimes have net force and torque components, described by first-rank and asymmetric second-rank moment tensors, which must be included in analyses of landslides and some volcanic phenomena. Non-DC processes that lead to conventional (symmetric second-rank) moment tensors include geometrically complex shear faulting, tensile faulting, shear faulting in an anisotropic medium, shear faulting in a heterogeneous region (e.g., near an interface), and polymorphic phase transformations. Undoubtedly, many non-DC earthquake processes remain to be discovered. Progress will be facilitated by experimental studies that use wave amplitudes, amplitude ratios, and complete waveforms in addition to wave polarities and thus avoid arbitrary assumptions such as the absence of volume changes or the temporal similarity of different moment tensor components.

Reviews of Geophysics↗

Non-double-couple earthquakes 2. Observations

Most studies assume that earthquakes have double-couple (DC) source mechanisms, corresponding to shear motion on planar faults. However, many well-recorded earthquakes have radiation patterns that depart radically from this model, indicating fundamentally different source processes. Seismic waves excited by advective processes, such as landslides and volcanic eruptions, are consistent with net forces rather than DCs. Some volcanic earthquakes also have single-force mechanisms, probably because of advection of magmatic fluids. Other volcanic earthquakes have mechanisms close to compensated linear vector dipoles and may be caused by magmatic intrusions. Shallow earthquakes in volcanic or geothermal areas and mines often have mechanisms with isotropic components, indicating volume changes of either explosive or implosive polarity. Such mechanisms are consistent with failure involving both shear and tensile faulting, which may be facilitated by high-pressure, high-temperature fluids. In mines, tunnels are cavities that may close. Deep-focus earthquakes occur within zones of polymorphic phase transformations in the upper mantle at depths where stick-slip instability cannot occur. Their mechanisms tend to be deviatoric (volume conserving), but non-DC, and their source processes are poorly understood. Automatic global moment tensor services routinely report statistically significant non-DC components for large earthquakes, but detailed reexamination of individual events is required to confirm such results.

Reviews of Geophysics↗

Thermal maturity of tectonostratigraphic terranes within the Franciscan Complex, California

Indicators of organic metamorphism provide valuable tools for analyzing the thermal history of tectonostratigraphic terranes. Paleotemperature estimates derived from vitrinite reflectance, for example, are more precise than values based upon inorganic mineral assemblages in low‐grade rocks. Isothermal geometries must be interpreted within the context of structural and stratigraphic data, but, by doing so, direct comparisons of thermal history can be made between adjacent terranes. Data from the Franciscan Complex of California show that three basic types of thermal anomalies are related to terrane amalgamation: (1) maturity inversion across terrane boundaries, (2) shear heating along terrane boundaries, and (3) resetting of maturity gradients via thrust loading. In addition, thermal overprints result if there is a substantial rise in local or regional geothermal gradient following terrane amalgamation; under these circumstances, homogenization of thermal maturity may be evident across terrane boundaries. If associated magmatic rocks or hydrothermal veins can be dated, an upper limit is established on the suturing event. Thermal data are also useful for resolving structural ambiguities on a local scale; clear distinctions can be made, for example, between faulted contacts, diapirs, and mildly tectonized depositional contacts or unconformities. Citing Literatu

California↗

Variation of depth to the brittle-ductile transition due to cooling of a midcrustal intrusion

The depth to the brittle-ductile transition in the crust is often defined by the intersection of a shear resistance relation in the brittle upper crust that increases linearly with depth and a power law relation for ductile flow in the lower crust that depends strongly on temperature. Transient variation of this depth caused by a magmatic intrusion at a depth near the regional transition can be modeled by a heat conduction model for a rectangular parallelepiped superimposed on a linear geothermal gradient. When parameters appropriate for the southeastern United States are used, a moderate-sized intrusion is found to decrease the transition depth by as much as 7 km; significant variations last about 10 m.y. Since the base of the seismogenic zone is identified with the brittle-ductile transition and since significant depth variations provide loci for strain energy concentration, these results imply that intrusions of late Tertiary age or younger could be important sources of clustered seismicity.

Geophysical Research Letters↗