Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geophysical”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Stochastic Inversion of Gaussian Random Media Using Transverse Coherence Functions for Reflected Waves: Theory and Method

The transverse coherence functions (TCFs) of phase and amplitude fluctuations of a seismic wave are powerful to estimate the spatial distribution, length scales, and strength of random heterogeneities. However, TCFs have been formulated for transmitted waves only, not for reflected waves. In this paper, we derive reflection TCFs for Gaussian random media. Furthermore, we propose to invert for Gaussian random media using the reflection TCFs based on the grid search. We validate the new reflection TCF formulas using 2D finite-difference numerical experiments. The numerical example also illustrates the feasibility and efficiency of the inversion. The stochastic inversion using reflected waves can be used in both exploration and global seismology.

58 GEOSCIENCES↗
DOE OSTI2021

The Role of Stratigraphy and Loading History in Generating Complex Compaction Bands in Idealized Field-Scale Settings

The Buckskin Gulch locality in Utah is a landmark example of compaction localization. The outcrop of this locality involves distinct stratigraphic heterogeneity and was exposed to complex loading history. It features multiple sets of deformation bands with different kinematics and orientation. Similar formations were seen in the Valley of Fire, Nevada, and the Orange quarry, France, among other localities. The formation of such complex structures, their propagation mechanisms, and frequency is affected by numerous local and ambient factors whose impacts are not yet fully understood. The simulation of the above-mentioned localities is not feasible because of the limited amount of available information. This work, instead, investigates from a geomechanics standpoint how the interplay among material nonlinearity, outcrop stratigraphy, and loading history interconnects with specific spatiotemporal patterns of compaction band propagation. Our study shows that the system stratigraphy can be responsible for the emergence of coexisting compaction bands with different inclination and kinematics. Specifically, we show that stiffness contrasts induce nonlocal stress changes which may favor the initiation of secondary structures with different compaction localization characteristics. Furthermore, systems of inclined compaction bands induced by burial increase display secondary, noncontemporaneous sets of vertical compaction bands under the effects of postburial shortening. Our results indicate that stages of intermediate burial decrease prior to tectonic shortening can promote the formation of such complex systems. Despite the simplifications involved in our analyses, these findings show how geomechanics computations complement field observations and could provide a mechanics-based validation of site-specific reconstruction hypothesis.

58 GEOSCIENCES↗
DOE OSTI2020

Intrinsic Elastic Anisotropy of Westerly Granite Observed by Ultrasound Measurements, Microstructural Investigations, and Neutron Diffraction

Westerly granite (WG) has been generally accepted as an isotropic homogeneous rock. Here, we return to WG and observe significant elastic anisotropy using multidirectional ultrasonic sounding on spherical samples at pressures up to 400 MPa. Thermal treatment of WG leads to formation of microcracks that reduce elastic wave velocities and increase its elastic anisotropy. The 3D distribution of P-wave velocities at low pressure is close to orthorhombic symmetry. Application of hydrostatic pressure closes most of thermally induced microcracks and decreases elastic anisotropy of WG, but at high pressure the anisotropy is practically reversed compared to low pressure: maximum P-wave velocity direction at low pressures is near minimum velocity direction at high pressure and vice versa. To understand this effect, microstructures of the rock were investigated by optical and scanning electron microscopy. Preferred orientations of four major rock-forming minerals—quartz, orthoclase, plagioclase, and biotite—were measured by time-of-flight neutron diffraction, which confirms significant crystal alignment. All these data were used to numerically model anisotropic elastic properties of WG. It is shown that WG possesses weak intrinsic elastic anisotropy related mainly to the preferred orientation of feldspars formed during igneous crystallization. Observed microcracks are mostly related to the cleavage planes of feldspars and biotite, and thus also demonstrate preferred orientation. Higher preheating temperatures produce larger quantity of longer microcracks. A numerical model shows that these microcracks act against the weak intrinsic elastic anisotropy of WG, and define the elastic anisotropy at low pressures.

58 GEOSCIENCES↗
DOE OSTI2021

Methane Hydrate Formation and Evolution During Sedimentation

We explored methane hydrate formation with sedimentation with a newly developed one-dimensional, multiphase flow, multicomponent transport numerical model. Our model couples methane hydrate formation from in situ microbial methane generation within the hydrate stability zone (HSZ), methane recycling, and microbial methane generation below the base of the hydrate stability zone (BHSZ). Both recycled methane and deeply generated methane are transported into the HSZ by buoyancy-driven free gas flow. Free gas flows through the HSZ by both the processes of capillary-dependent pore fillings and by salt exclusion during hydrate formation, with the former being the dominant mechanism. We quantitively illustrated the formation of enriched hydrate in muddy sediments above, and interconnected free gas below, the BHSZ, which are common features along the world's continental margin. In addition, we showed two ways to form concentrated methane hydrate above the BHSZ. The first mechanism is local free gas flow during methane recycling. This happens at sites with sufficient methane generation above the BHSZ. The second mechanism is deep microbial methane generation which is transported into the HSZ by free gas flow. This mechanism plays a more important role at sites with high sedimentation rates. This study provides new insights into methane hydrate formation and distribution below the seafloor. It is important for understanding the carbon cycle and carbon storage below the seafloor and for resource evaluation and exploitation.

58 GEOSCIENCES↗
DOE OSTI2021

Regional Subsidence Analysis Through a Multi-Scale Modeling Framework Based on Breakage Mechanics

Although poroelastic models are often used to explain the delay between subsurface fluid depletion and ground subsidence, inelastic compaction involving permanent changes of rock microstructure may exacerbate hydro-mechanical coupling, thus influencing the interpretation of measurements and long-term forecasts. Here, a multi-scale modeling approach is discussed, which accounts for the inherent connection between rock microstructure, hydraulic conductivity, and pore compaction. A constitutive model built within the framework of breakage mechanics is proposed to link the hydraulic conductivity of granular rocks with inelastic deformations and changes in grading caused by injection-depletion cycles at stress levels far from yielding. The proposed model has been incorporated into large-scale simulation frameworks, thus enabling the spatiotemporal mapping of regional subsidence through a hybrid, semi-analytical approach. Numerical results based on this strategy show that the model allows isolating near-field and far-field effects into the computation of land subsidence and can generate forecasts for different modeling scenarios (e.g., elastic and inelastic compaction, constant permeability, and concurrent change of compressibility and permeability). In particular, examples of simulations for the case of the Groningen gas field are discussed, showing the model capabilities to use both field measurements and laboratory tests for the generation of reasonable subsidence maps, without expensive computational costs. Results indicate that ignoring coupled inelastic effects has major consequences on the predicted timescale of subsidence. Specifically, while all the model scenarios produced similar long-term ground settlements, those ignoring breakage-dependent permeability changes result in a variation of the temporal window of residual subsidence of the order of several decades.

58 GEOSCIENCES↗
DOE OSTI2021

Precambrian Crystalline Basement Properties From Pressure History Matching and Implications for Induced Seismicity in the US Midcontinent

Injection-induced seismicity across the US midcontinent has almost exclusively occurred in the crystalline basement that underlies the Arbuckle Group aquifer and its equivalents, the primary wastewater disposal zone in this region. However, the properties of the basement are not well known. Newly compiled data, from Class I wells in Kansas, provide a unique record of pressures in the Arbuckle and an opportunity to constrain the reservoir-scale properties of the basement such as permeability, diffusivity, and specific storage. Constraints on these parameters are critical for modeling fluid flow and pressures across the entire Arbuckle-basement system, and are necessary for accurate evaluation and prediction of injection-induced earthquakes. Here, we present a detailed, three-dimensional geological and pressure history-matched numerical model for the Arbuckle and basement, based on data from >400 wells covering a large region in south-central Kansas, where injection-induced seismicity has been concentrated since 2014. Simulations of dynamic data from 319 wells indicate that Arbuckle pressures have increased by 1.1 MPa in high injection rate areas and an overpressure of <0.1 MPa may be the cause of seismicity in the basement. Pressure-history matching also yields the likely range in porosity (0.3%–7%), permeability (0.1–0.7 mD), and diffusivity (0.004–0.07 m2 /s) for the basement. The resulting estimates suggest reservoir-scale properties of the basement are enhanced by faults and fractures. Importantly, the diffusivities determined in this study are lower than estimates derived from Kansas earthquake triggering fronts, and suggest that such seismicity-based techniques may have limitations, particularly where spacetime patterns between injection and seismicity are complex.

58 GEOSCIENCES↗
DOE OSTI2022

Multidisciplinary Constraints on the Thermal‐Chemical Boundary Between Earth's Core and Mantle

Abstract 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.

58 GEOSCIENCES↗
DOE OSTI2021

Isolating Detrital and Diagenetic Signals in Magnetic Susceptibility Records From Methane‐Bearing Marine Sediments

Abstract Volume‐dependent magnetic susceptibility (κ) is commonly used for paleoenvironmental reconstructions in both terrestrial and marine sedimentary environments where it reflects a mixed signal between primary deposition and secondary diagenesis. In the marine environment, κ is strongly influenced by the abundance of ferrimagnetic minerals regulated by sediment transport processes. Post‐depositional alteration by H 2 S, however, can dissolve titanomagnetite, releasing reactive Fe that promotes pyritization and subsequently decreases κ. Here, we provide a new approach for isolating the detrital signal in κ and identifying intervals of diagenetic alteration of κ driven by organoclastic sulfate reduction (OSR) and the anaerobic oxidation of methane (AOM) in methane‐bearing marine sediments offshore India. Using the correlation of a heavy mineral proxy from X‐ray fluorescence data (Zr/Rb) and κ in unaltered sediments, we predict the primary detrital κ signal and identify intervals of decreased κ, which correspond to increased total sulfur content. Our approach is a rapid, high‐resolution method that can identify overprinted κ resulting from pyritization of titanomagnetite due to H 2 S production in marine sediments. In addition, total organic carbon, total sulfur, and authigenic carbonate δ 13 C measurements indicate that both OSR and AOM can drive the observed κ loss, but AOM drives the greatest decreases in κ. Overall, our approach can enhance paleoenvironmental reconstructions and provide insight into paleo‐positions of the sulfate‐methane transition zone, past enhancements of OSR or paleo‐methane seepage, and the role of detrital iron oxide minerals on the marine sediment sulfur sink, with consequences influencing the development of chemosynthetic biological communities at methane seeps.

58 GEOSCIENCES↗
DOE OSTI2021

Episodic Earthquake Swarms in the Mineral Mountains, Utah Driven by the Roosevelt Hydrothermal System

Over 1,000 earthquakes (-2.0 < M < 2.0), identified using a matched-filter method, occurred in the Mineral Mountains, Utah between 2016 and 2019. The enhanced catalog is complete down to M -0.9 and contains roughly 15 times more events than originally cataloged. Earthquake relocation of ~800 earthquakes shows that activity is concentrated in a <2 km long E-W striking narrow zone, ~4 km east of the Roosevelt hydrothermal system. Two fault orientations, both N-S and E-W parallel to the Opal Mound and Mag Lee faults, respectively, are observed after computing composite focal mechanisms of highly similar earthquakes. Looking solely at the temporal distribution of the seismicity, we identify 15 periods of swarm-like activity, with two major clusters occurring in December 2016, recorded by three stations, and in October 2019 recorded by eight stations. The October 2019 swarm, the best recorded sequence in the area, provides evidence for the underlying triggering mechanism. We show that a complex mechanism of fluid diffusion and aseismic slip is responsible for the swarm evolution with migration velocities reaching 10 km/day. We hypothesize that these episodic swarms in the Mineral Mountains are primarily driven by migrating fluids that originate within the Roosevelt hydrothermal system.

58 GEOSCIENCES↗
DOE OSTI2021

Comparison of Biotite Elastic Properties Recovered by Spherical Nanoindentations and Atomistic Simulations — Influence of Nano-Scale Defects in Phyllosilicates

Phyllosilicate minerals, due to their sheets structure and morphology, are known to cause anisotropy in bulk rock properties and make the bulk rock more compliant. Accurately characterizing the micromechanical behavior of phyllosilicate minerals from laboratory observations, which eventually translates to the bulk rock behavior, is still challenging due to their fine-grained nature. Recent advances in atomistic simulations open the possibility of theoretically investigating such mineral mechanical behavior. We compare the elastic properties of biotites recovered by spherical nanoindentation with those predicted from density functional theory (DFT) simulations to investigate to what extent theoretical predictions reproduce actual phyllosilicate properties. Spherical nanoindentation was conducted using schist rocks from Poorman Formation, South Dakota, USA, to recover continuous indentation stress-strain curves. Loading in the layer-normal orientation shows an average indentation modulus (M) of about 35 GPa, while loading in the layer-parallel orientation gives a higher average of about 95 GPa. To facilitate comparison, the elastic stiffness constants (cij) determined from DFT were converted to indentation modulus (M) using solutions proposed in this study. The majority of the nanoindentation modulus results are below the values inferred from the simulation results representing ideal defect-free minerals. We suggest that crystal defects present at the nano-scale, potentially ripplocations, are the dominant cause of the lower indentation modulus recovered from nanoindentation compared to those inferred from DFT simulations. Results highlight the importance of acknowledging the defects that exist down to the nano-scale as it modifies the mechanical properties of phyllosilicates compared to its pure defect-free form.

58 GEOSCIENCES↗
DOE OSTI2022

Monitoring Fracture Saturation With Internal Seismic Sources and Twin Neural Networks

Seismic coda-wave analysis is a well-developed method for detecting subtle physical changes in complex media by measuring arrival times in the late-arriving energy from multiply scattered or reflected waves. However, a challenge arises when multiply scattered waves are not sufficiently separated in time from the direct arrivals to provide a clear coda wave train. Additional complications for monitoring changes in fracture systems arise when the signals originate from unsynchronized internal sources, such as natural or induced seismicity, from acoustic emission, or from transportable intra-fracture sources (chattering dust), that generate uncontrolled signals that vary in arrival time, amplitude and frequency content. Here, we use a twin neural network (TNN also known as a Siamese neural network) for dimensionality reduction to analyze signals from chattering dust to classify the fluid saturation state of a synthetic fracture system. The TNN with shared weights generates a low-dimensional representation of the data input by minimizing contrastive loss, serving as the input to a multiclass classifier that accurately classifies whether multiple fractures in a fracture system are fully saturated or partially saturated, or whether a change in saturation has occurred in different fractures in the system. Furthermore, these results show that information buried in unresolved codas from uncontrolled sources can be extracted using machine learning to monitor the evolution of fracture systems caused by physical and chemical processes even when the scattered and direct wave fields overlap.

58 GEOSCIENCES↗
DOE OSTI2022

An End-To-End Earthquake Detection Method for Joint Phase Picking and Association Using Deep Learning

Earthquake monitoring by seismic networks typically involves a workflow consisting of phase detection/picking, association, and location tasks. In recent years, the accuracy of these individual stages has been improved through the use of machine learning techniques. Here, in this study, we introduce a new end-to-end approach that improves overall earthquake detection accuracy by jointly optimizing each stage of the detection pipeline. We propose a neural network architecture for the task of multi-station processing of seismic waveforms recorded over a seismic network. This end-to-end architecture consists of three sub-networks: a backbone network that extracts features from raw waveforms, a phase picking sub-network that picks P- and S-wave arrivals based on these features, and an event detection sub-network that aggregates the features from multiple stations to associate and detect earthquakes across a seismic network. We use these sub-networks together with a shift-and-stack module based on back-projection that introduces kinematic constraints on arrival times, allowing the neural network model to generalize to different velocity models and to variable station geometry in seismic networks. We evaluate our proposed method on the STanford EArthquake Dataset (STEAD) and on the 2019 Ridgecrest, CA earthquake sequence. The results demonstrate that our end-to-end approach can effectively pick P- and S-wave arrivals and achieve earthquake detection accuracy rivaling that of other state-of-the-art approaches. Because our approach preserves information across tasks in the detection pipeline, it has the potential to outperform approaches that do not.

58 GEOSCIENCES↗
DOE OSTI2022

From Lake to River: Documenting an Environmental Transition Across the Jura/Knockfarril Hill Members Boundary in the Glen Torridon Region of Gale Crater (Mars)

Between January 2019 and January 2021, the Mars Science Laboratory team explored the Glen Torridon (GT) region in Gale crater (Mars), known for its orbital detection of clay minerals. Mastcam, Mars Hand Lens Imager, and ChemCam data are used in an integrated sedimentological and geochemical study to characterize the Jura member of the upper Murray formation and the Knockfarril Hill member of the overlying Carolyn Shoemaker formation in northern GT. The studied strata show a progressive transition represented by interfingering beds of fine–grained, recessive mudstones of the Jura member and coarser–grained, cross–stratified sandstones attributed to the Knockfarril Hill member. Whereas the former are interpreted as lacustrine deposits, the latter are interpreted as predominantly fluvial deposits. The geochemical composition seen by the ChemCam instrument show K 2 O–rich mudstones (~1–2 wt.%) versus MgO–rich sandstones (>6 wt.%), relative to the average composition of the underlying Murray formation. We document consistent sedimentary and geochemical data sets showing that low–energy mudstones of the Jura member are associated with the K–rich endmember, and that high–energy cross–stratified sandstones of the Knockfarril Hill member are associated with the Mg–rich endmember, regardless of stratigraphic position. The Jura to Knockfarril Hill transition therefore marks a significant paleoenvironmental change, where a long–lived and comparatively quiescent lacustrine setting progressively changes into a more energetic fluvial setting, as a consequence of shoreline regression due to either increased sediment supply or lake–level drop.

58 GEOSCIENCES↗
DOE OSTI2021

Concerning the Aims and Scope for JAMES

In this editorial, we present general principles as well as specific notions that guide the strategy of JAMES' editors in realizing the journal's mission.

Journal, aims and scope, American Geophysics Union↗
DOE OSTI2023

Equilibrated Gas and Carbonate Standard‐Derived Dual (Δ 47 and Δ 48 ) Clumped Isotope Values

Abstract Carbonate clumped isotope geochemistry has primarily focused on mass spectrometric determination of m/z 47 CO 2 for geothermometry, but theoretical calculations and recent experiments indicate paired analysis of the m/z 47 ( 13 C 18 O 16 O) and m/z 48 ( 12 C 18 O 18 O) isotopologues (referred to as Δ 47 and Δ 48 ) can be used to study non‐equilibrium isotope fractionations and refine temperature estimates. We utilize 5,448 Δ 47 and 3,400 Δ 48 replicate measurements of carbonate samples and standards, and 183 Δ 47 and 195 Δ 48 replicate measurements of gas standards from 2015 to 2021 from a multi‐year and multi‐instrument data set to constrain Δ 47 and Δ 48 values for 27 samples and standards, including Devils Hole cave calcite, and study equilibrium Δ 47 ‐Δ 48 , Δ 47 ‐temperature, and Δ 48 ‐temperature relationships. We compare results to previously published findings and calculate equilibrium regressions based on data from multiple laboratories. We report acid digestion fractionation factors, Δ* 63‐47 and Δ* 64‐48 , and account for their dependence on the initial clumped isotope values of the mineral.

58 GEOSCIENCES↗
DOE OSTI2022

Size Ranges of Magnetic Domain States in Tetrataenite

Paleomagnetic studies of meteorites provide unique constraints on the evolution of magnetic fields in the early solar system. These studies rely on the identification of magnetic minerals that can retain stable magnetizations over ≳4.5 billion years (Ga). The ferromagnetic mineral tetrataenite (γ''-Fe 0.5 Ni 0.5 ) is found in iron, stony-iron and chondrite meteorite groups. Nanoscale intergrowths of tetrataenite have been shown to carry records of paleomagnetic fields, although the effect of magnetostatic interactions on their magnetic remanence acquisition remains to be fully understood. Tetrataenite can also occur as isolated, non-interacting, nanoscale grains in many meteorite groups, although the paleomagnetic potential of these grains is particularly poorly understood. Here, we aim to improve our understanding of tetrataenite magnetization to refine our knowledge of existing paleomagnetic analyses and broaden the spectrum of meteorite groups that can be used for future paleomagnetic studies. We present the results of analytical calculations and micromagnetic modeling of isolated tetrataenite grains with various geometries. We find that tetrataenite forms a stable single domain state at grain lengths between 6 and ~160 nm dependent on its elongation. It also possesses a magnetization resistant to viscous remagnetization over the lifetime of the solar system at 293 K. At larger grain sizes, tetrataenite's lowest energy state is a lamellar two-domain state, stable at Ga-scale timescales. Unlike many other magnetic minerals, tetrataenite does not form a single-vortex domain state due to its large uniaxial anisotropy. Our results show that single domain and two-domain tetrataenite grains carry an extremely stable magnetization and therefore are promising for paleomagnetic studies.

58 GEOSCIENCES↗
DOE OSTI2023

Salt Diapir‐Driven Recycling of Gas Hydrate

Abstract By harnessing both hypothetical, synthetic basin and gas hydrate (GH) system models and real‐world models of well‐studied salt diapir‐associated GH sites at Green Canyon (Gulf of Mexico) and Blake Ridge (U.S. Atlantic coast), we propose and demonstrate salt movement (and in particular, diapirism) to be a new mechanism for the recycling of marine GH. At Green Canyon, for example, we show that by considering this newly proposed diapir‐driven recycling mechanism in conjunction with previously proposed lithological control on sandy‐reservoir‐hosted hydrate at the base of the GH stability zone (BGHSZ; ∼bottom‐simulating reflector, BSR), modeled GH saturations match drilling data. Overall, salt diapir movement‐induced GH recycling provides a temperature‐driven mechanism by which GH saturations at the BGHSZ may reach >90 vol. % and by which GH volumes near and free gas volumes beneath the BGHSZ may be increased significantly through time. Interestingly, comparison of salt diapir‐driven recycling and sediment burial‐driven recycling scenarios suggests notably higher rates of recycling via diapir‐driven versus burial‐driven processes. Our results suggest that GH and associated free gas accumulations above salt diapir crests represent particularly attractive targets for unconventional and conventional hydrocarbon resource exploration and for scientific and academic drilling expeditions aimed at exploiting GH systems. Salt basins containing GH systems—including passive margin basins of the Gulf of Mexico, southeastern Brazil, and southwestern Africa—are therefore compelling localities for studying salt‐driven GH recycling and for salt diapir‐associated natural gas exploration.

58 GEOSCIENCES↗
DOE OSTI2022

Frictional Behavior Downdip Along the Subduction Megathrust: Insights From Laboratory Experiments on Exhumed Samples at In Situ Conditions

Subduction megathrusts exhibit a spectrum of slip modes, including catastrophic earthquakes. Although the mechanical and frictional properties of materials sampled from subduction zones have been studied extensively, few datasets have been collected for compositions and at pressure and temperature conditions representative of those in situ. The Nankai subduction zone in southwest Japan is a well-studied margin, and abundant data provide an opportunity to advance our understanding of fault and earthquake physics. Here, we use samples exhumed in the Shimanto and Sanbagawa Belts on Shikoku Island of southwest Japan that represent analogs for materials along the present-day megathrust at depths of ~5–>25 km, and we shear these at their peak in situ pressure-temperature (P-T) conditions. Effective normal stresses range from 28 to 192 MPa, and temperatures from 105°C to 470°C. We used pore fluid pressures of 45–240 MPa, corresponding to fluid overpressure ratios λ of 0.65 and 0.90. Slip velocities of 0.1–100 μm/s were used, in order to focus on the nucleation of instability and earthquakes. We found predominantly velocity-strengthening (inherently stable) behavior under all conditions for λ = 0.65. For λ = 0.90, velocity-weakening behavior was observed at 350°C, with velocity-strengthening behavior at lower and higher temperatures. The rate/state frictional stability parameter (a-b) increases with slip velocity at temperatures up to ~200°C and remains constant or decreases with slip velocity at higher temperatures. Altogether, our results demonstrate the potentially important roles of both temperature and slip velocity in controlling the distribution of stress and frictional rheology along subduction thrusts.

58 GEOSCIENCES↗