Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “forward modeling”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS2024

Spectral Correlation in MODIS Water-Leaving Reflectance Retrieval Uncertainty

Spectral remote sensing reflectance, Rrs(λ) (sr−1), is the fundamental quantity used to derive a host of bio-optical and biogeochemical properties of the water column from satellite ocean color measurements. Estimation of uncertainty in those derived geophysical products is therefore dependent on knowledge of the uncertainty in satellite-retrieved R rs . Furthermore, since the associated algorithms require R rs at multiple spectral bands, the spectral (i.e., band-to-band)error covariance in R rs is needed to accurately estimate the uncertainty in those derived properties. This study establishes a derivative-based approach for propagating instrument random noise, instrument systematic uncertainty, and forward model uncertainty into R rs as retrieved using NASA’s multiple-scattering epsilon (MSEPS) atmospheric correction algorithm, to generate pixel-level error covariance in R rs . The approach is applied to measurements from Moderate Resolution Imaging Spectroradiometer (MODIS) on the Aqua satellite and verified using Monte Carlo (MC) analysis. We also make use of this full spectral error covariance in R rs to calculate uncertainty in phytoplankton pigment chlorophyll-a concentration (chl a , mg/m 3 ) and diffuse attenuation coefficient of downwelling irradiance at 490 nm (K d (490), m -1 ). Accounting for the error covariance in R rs generally reduces the estimated relative uncertainty in chl a by ∼1-2% (absolute value) in waters with chl a < 0.25 mg/m 3 where the color index (CI) algorithm is used. The reduction is ∼5-10% in waters with chl a > 0.35 mg/m 3 where the blue-green ratio (OCX) algorithm is used. Such reduction can be higher than 30% in some regions. For K d (490), the reduction by error covariance is generally ∼2%, but can be higher than 20% in some regions. The error covariance in R rs is further verified through forward-calculating chl a from MODIS-retrieved and in situ R rs and comparing estimated uncertainty with observed differences. An 8-day global composite of propagated uncertainty shows that the goal of 35% uncertainty in chl a can be achieved over deep ocean waters (chl a ≤ 0.1 mg/m3). While the derivative-based approach generates reasonable error covariance in R rs some assumptions should be updated as our knowledge improves. These include the inter-band error correlation in top-of-atmosphere reflectance, and uncertainties in the calibration of MODIS 869 nm band, in ancillary data, and in the in situ data used for system vicarious calibration.

Ocean color↗
NASA NTRS2014

The Architecture of the LkCa 15 Transitional Disk Revealed By High-Contrast Imaging

We present four new epochs of Ks-band images of the young pre-transitional disk around LkCa 15 and perform extensive forward modeling to derive the physical parameters of the disk. We find indications of strongly anisotropic scattering (Henyey-Greenstein g = 0.67 (+0.18/−0.11)) and a significantly tapered gap edge ("round wall") but see no evidence that the inner disk, whose existence is predicted by the spectral energy distribution, shadows the outer regions of the disk visible in our images.We marginally confirm the existence of an offset between the disk center and the star along the line of nodes; however, the magnitude of this offset (x = 27 (+19/−20) mas) is notably lower than that found in our earlier H-band images. Intriguingly, we also find an offset of y = 69 (+49/−25) mas perpendicular to the line of nodes at high significance. If confirmed by future observations, this would imply a highly elliptical - or otherwise asymmetric - disk gap with an effective eccentricity of e ≈ 0.3. Such asymmetry would most likely be the result of dynamical sculpting by one or more unseen planets in the system. Finally, we find that the bright arc of scattered light we see in direct imaging observations originates from the near side of the disk and appears brighter than the far side because of strong forward scattering.

dynamical sculpting↗
NASA NTRS2014

The Architecture of the LkCa 15 Transitional Disk Revealed by High-contrast Imaging

We present four new epochs of Ks-band images of the young pre-transitional disk around LkCa 15, and perform extensive forward modeling to derive the physical parameters of the disk. We find indications of strongly anisotropic scattering (Henyey-Greenstein g = 0.67+0.18 −0.11) and a significantly tapered gap edge ('round wall'), but see no evidence that the inner disk, whose existence is predicted by the spectral energy distribution, shadows the outer regions of the disk visible in our images. We marginally confirm the existence of an offset between the disk center and the star along the line of nodes; however, the magnitude of this offset (x = 27+19 −20 mas) is notably lower than that found in our earlier H-band images (Thalmann et al. 2010). Intriguingly, we also find, at high significance, an offset of y = 69+49 −25 mas perpendicular to the line of nodes. If confirmed by future observations, this would imply a highly elliptical- or otherwise asymmetric-disk gap with an effective eccentricity of e ≈ 0.3. Such asymmetry would most likely be the result of dynamical sculpting by one or more unseen planets in the system. Finally, we find that the bright arc of scattered light we see in direct imaging observations originates from the near side of the disk, and appears brighter than the far side because of strong forward scattering.

anisotropic scattering↗
DOE OSTI2021

Predicting images for the dynamics of stellar clusters ( π-DOC ): a deep learning framework to predict mass, distance, and age of globular clusters

ABSTRACT Dynamical mass estimates of simple systems such as globular clusters (GCs) still suffer from up to a factor of 2 uncertainty. This is primarily due to the oversimplifications of standard dynamical models that often neglect the effects of the long-term evolution of GCs. Here, we introduce a new approach to measure the dynamical properties of GCs, based on the combination of a deep-learning framework and the large amount of data from direct N-body simulations. Our algorithm, π-DOC (Predicting Images for the Dynamics Of stellar Clusters) is composed of two convolutional networks, trained to learn the non-trivial transformation between an observed GC luminosity map and its associated mass distribution, age, and distance. The training set is made of V-band luminosity and mass maps constructed as mock observations from N-body simulations. The tests on π-DOC demonstrate that we can predict the mass distribution with a mean error per pixel of 27 per cent, and the age and distance with an accuracy of 1.5 Gyr and 6 kpc, respectively. In turn, we recover the shape of the mass-to-light profile and its global value with a mean error of 12 per cent, which implies that we efficiently trace mass segregation. A preliminary comparison with observations indicates that our algorithm is able to predict the dynamical properties of GCs within the limits of the training set. These encouraging results demonstrate that our deep-learning framework and its forward modelling approach can offer a rapid and adaptable tool competitive with standard dynamical models.

Chardin, Jonathan↗
DOE OSTI2022

Evaluating seasonal and regional distribution of snowfall in regional climate model simulations in the Arctic

In this study, we investigate how the regional climate model HIRHAM5 reproduces the spatial and temporal distribution of Arctic snowfall when compared to CloudSat satellite observations during the examined period of 2007–2010. For this purpose, both approaches, i.e., the assessments of the surface snowfall rate (observation-to-model) and the radar reflectivity factor profiles (model-to-observation), are carried out considering spatial and temporal sampling differences. The HIRHAM5 model, which is constrained in its synoptic representation by nudging to ERA-Interim, represents the snowfall in the Arctic region well in comparison to CloudSat products. The spatial distribution of the snowfall patterns is similar in both identifying the southeastern coast of Greenland and the North Atlantic corridor as regions gaining more than twice as much snowfall as the Arctic average, defined here for latitudes between 66 and 81°N. Excellent agreement (difference less than 1%) in the Arctic-averaged annual snowfall rate between HIRHAM5 and CloudSat is found, whereas ERA-Interim reanalysis shows an underestimation of 45% and significant deficits in the representation of the snowfall rate distribution. From the spatial analysis, it can be seen that the largest differences in the mean annual snowfall rates are an overestimation near the coastlines of Greenland and other regions with large orographic variations as well as an underestimation in the northern North Atlantic Ocean. To a large extent, the differences can be explained by clutter contamination, blind zone or higher resolution of CloudSat measurements, but clearly HIRHAM5 overestimates the orographic-driven precipitation. The underestimation of HIRHAM5 within the North Atlantic corridor south of Svalbard is likely connected to a poor description of the marine cold air outbreaks which could be identified by separating snowfall into different circulation weather type regimes. By simulating the radar reflectivity factor profiles from HIRHAM5 utilizing the Passive and Active Microwave TRAnsfer (PAMTRA) forward-modeling operator, the contribution of individual hydrometeor types can be assessed. Looking at a latitude band at 72–73°N, snow can be identified as the hydrometeor type dominating radar reflectivity factor values across all seasons. The largest differences between the observed and simulated reflectivity factor values are related to the contribution of cloud ice particles, which is underestimated in the model, most likely due to the small sizes of the particles. The model-to-observation approach offers a promising diagnostic when improving cloud schemes, as illustrated by comparison of different schemes available for HIRHAM5.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Ingrained: An Automated Framework for Fusing Atomic-Scale Image Simulations into Experiments

To fully leverage the power of image simulation to corroborate and explain patterns and structures in atomic resolution microscopy, an initial correspondence between the simulation and experimental image must be established at the outset of further high accuracy simulations or calculations. Furthermore, if simulation is to be used in context of highly automated processes or high-throughput optimization, the process of finding this correspondence itself must be automated. In this work, "ingrained," an open-source automation framework which solves for this correspondence and fuses atomic resolution image simulations into the experimental images to which they correspond, is introduced. Here, the overall "ingrained" workflow, focusing on its application to interface structure approximations, and the development of an experimentally rationalized forward model for scanning tunneling microscopy simulation are described.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Using ultrasonic attenuation in cortical bone to infer distributions on pore size

Here, in this work we infer the underlying distribution on pore radius in human cortical bone samples using ultrasonic attenuation data. We first discuss how to formulate polydisperse attenuation models using a probabilistic approach and the Waterman Truell model for scattering attenuation. We then compare the Independent Scattering Approximation and the higher-order Waterman Truell models’ forward predictions for total attenuation in polydisperse samples. Following this, we formulate an inverse problem under the Prohorov Metric Framework coupled with variational regularization to stabilize this inverse problem. We then use experimental attenuation data taken from human cadaver samples and solve inverse problems resulting in nonparametric estimates of the probability density function on pore radius. We compare these estimates to the “true” microstructure of the bone samples determined via microCT imaging. We find that our methodology allows us to reliably estimate the underlying microstructure of the bone from attenuation data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

MADLens, a python package for fast and differentiable non-Gaussian lensing simulations

Here, we present MADLens a python package for producing non-Gaussian lensing convergence maps at arbitrary source redshifts with unprecedented precision. MADLens is designed to achieve high accuracy while keeping computational costs as low as possible. A MADLens simulation with only particles produces convergence maps whose power agrees with theoretical lensing power spectra up to within the accuracy limits of HaloFit. This is made possible by a combination of a highly parallelizable particle-mesh algorithm, a sub-evolution scheme in the lensing projection, and a machine-learning inspired sharpening step. Further, MADLens is fully differentiable with respect to the initial conditions of the underlying particle-mesh simulations and a number of cosmological parameters. These properties allow MADLens to be used as a forward model in Bayesian inference algorithms that require optimization or derivative-aided sampling. Another use case for MADLens is the production of large, high resolution simulation sets as they are required for training novel deep-learning-based lensing analysis tools. We make the MADLens package publicly available under a Creative Commons License

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

High-dimensional isotomics, part 2: Observations of over 100 constraints on methionine's isotome

The abundances of different isotopic forms of a compound, or isotopologues, will vary based on its physical and chemical history. The number of isotopologues increases combinatorically with the size of a molecule, and even small molecules such as amino acids have thousands of potentially observable isotopic variants. However, due to the analytical challenges of separating and observing isotopologues, only a few dimensions of isotopic diversity are routinely measured. Overcoming these challenges requires both an experimental method to observe many isotopic properties and a theoretical framework for interpreting these experiments. In Part 1, we presented such a theoretical framework; here, we demonstrate an experimental method, which we apply to methionine. Our approach uses a Q Exactive HF Orbitrap to perform several “M + N experiments”, where a sample is ionized, a subset of its isotopologues with cardinal mass N daltons greater than the unsubstituted isotopologue is selected and fragmented, and the proportions of all detectable isotopic forms of those fragment ions are quantified. We perform M + 1, M + 2, M + 3, and M + 4 experiments of a methionine sample and standard where the sample has a 100 ‰ enrichment of 13C at the methyl carbon relative to the natural 13C abundance at that position in the standard, and is otherwise identical to the standard. We observe isotopic forms of 8 fragment ion species for each version of the M + N experiment. With the assistance of a forward model of expected mass spectra, we identify isotopic peaks for each fragment ion based on observed mass and abundance, screen these for data quality, and quantify abundances for 146 unique isotopic peaks at precisions of ≈ 0.3–3 ‰. We present our direct observations and use them to reconstruct the concentrations of 19 individual singly, doubly, and triply-substituted isotopologues; doing so gives fewer constraints and broader error bars than working with the direct observations, but may be more interpretable for some applications. We also examine possibilities for measuring additional peaks, which are primarily limited by the detection limit of the Orbitrap-IRMS method. We then suggest some possible uses of our direct measurements for chemical forensics and hypothesis testing. Furthermore, our results demonstrate the diversity of isotopic constraints currently observable and interpretable for organic molecules.

58 GEOSCIENCES↗
DOE OSTI2024

Predicting non-linear stress–strain response of mesostructured cellular materials using supervised autoencoder

Recent breakthroughs in advanced manufacturing capabilities have made it possible to design and print sophisticated topologies of cellular structures using diverse engineering materials such as metals, polymers, and ceramics. In these architectured materials, it is often desirable to tailor the mechanical properties by altering the unit cell topology. This necessitates an in-depth understanding of how the topology of the unit cell structure affects the macroscopic behavior of the material in both the linear and the non-linear regimes encountered under large compression. Here, we have developed a machine learning (ML) approach capable of accelerating the prediction of the stress–strain response of a polymer-based cellular structure under uniaxial confined compression. As part of generating the training data for ML, 60,000 mesostructures were generated using a relatively novel approach based on cellular automata, and their corresponding stress–strain responses were obtained from the finite element simulations. Principal component analysis (PCA) was used to reduce the dimensionality of the stress–strain curves. With only 20 principal components, PCA captured 99.89% of the variance in the stress–strain curves while reducing the dimensionality by 5X. ML using supervised autoencoder was able to successfully speed up the prediction of the non-linear stress–strain response of a unit cell by up to 4600X. The proposed method can serve as an efficient data generation tool and a rapid means for predicting the structure–property relationship through accelerated forward modeling of cellular materials under compaction, in cases where the macroscopic stress–strain response is governed by the unit-cell topology.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Noble gas insights into early impact delivery and volcanic outgassing to Earth's atmosphere: A limited role for the continental crust

Earth's atmosphere, crust and mantle have evolved together through continuous geochemical exchange throughout Earth's history. Constraints on the transport of volatile elements and compounds between these reservoirs are crucial for understanding how Earth could have stayed habitable for extended periods of time. In this work, we present a new forward model of He, Ne and Ar in the mantle, crust and atmosphere. We explore concentrations of noble gases at the end of accretion, bulk silicate Earth K/U ratios, crustal growth scenarios, and upper mantle processing rates throughout Earth's history. We search for parameter combinations that simultaneously satisfy observational constraints on present-day mantle 4 He/ 3 He and 40 Ar/ 36 Ar (sensitive to mantle outgassing and continental crust growth, which depletes the mantle of U, Th and K), atmospheric 20 Ne/ 22 Ne (which tracks the mix of outgassed vs. delivered Ne in the atmosphere), and atmospheric 40 Ar/ 36 Ar in the past and today (sensitive to volatile delivery, mantle outgassing and crustal growth). Leveraging this intertwined set of noble gas abundances and isotopic compositions yields new constraints on initial noble gas abundances in the mantle and on the proportions of atmospheric volatiles originating from delivery by impact degassing, mantle outgassing, and degassing of the continental crust. We find that atmospheric Ar isotopic evolution is primarily sensitive to the mantle processing rate history; the atmospheric Ar isotopic record should therefore not be used to reconstruct continental crust growth, but instead provides valuable insights into mantle processing rates. Our model predicts a measurably low 20 Ne/ 22 Ne ratio of ~9.7 in Archean atmospheric samples. Most of atmospheric primordial 36 Ar was directly delivered by chondritic bodies and not transferred to the atmosphere during an intense early episode of mantle outgassing. Nitrogen delivered by impact degassing could account for the present-day atmospheric nitrogen inventory.

58 GEOSCIENCES↗
DOE OSTI2020

Probabilistic analysis of masked loads with aggregated photovoltaic production

In this paper we present a probabilistic analysis framework to estimate behind-the-meter photovoltaic generation in real time. We develop a forward model consisting of a spatiotemporal stochastic process that represents the photovoltaic generation and a stochastic differential equation with jumps that represents the demand. Here, we employ this model to disaggregate the behind-the-meter photovoltaic generation using net load and irradiance measurements.

14 SOLAR ENERGY↗
DOE OSTI2024

Development of a scintillator based fast-ion loss detector for the Wendelstein 7-X stellarator

A new scintillator based fast-ion loss detector (FILD) system has been designed for the Wendelstein 7-X (W7-X) stellarator. The mechanical design of the system is presented here along with engineering analyses of the system. This includes an assessment of the structural loads which also considers electromagnetic forces calculated for possible disruption scenarios as well as a mode frequency analysis of the system. Furthermore, an analysis of the thermal characteristics of the actively cooled probe head under prescribed steady state and transient heat loads is presented. Finally, the optical relay system is described and its properties used to forward model synthetic camera signals using the FILDSIM code and based on theoretical fast-ion losses incident on the probe head calculated using the ASCOT5 code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The solubility of thorium in carbonate-bearing solutions at hydrothermal conditions

Thorium mineralization is frequently hosted in carbonate-bearing rocks, and thorium commonly substitutes into the structures of carbonate-bearing minerals that have precipitated from or been modified by hydrothermal fluids. Given this common association, it is reasonable to consider the hypothesis that the presence of carbonate ligands in hydrothermal solutions promotes the transport of Th through the formation of stable aqueous complexes. Our ability to evaluate this hypothesis, however, is hindered by the lack of experimental data for Th-carbonate species at conditions beyond ambient. The low-temperature data indicate that carbonate is a strong complexing agent for Th. In this contribution, we investigate the solubility of Th in carbonate-bearing fluids relevant to natural systems (0.05–0.5 m NaHCO 3 /Na 2 CO 3 ; pH T ~ 7.8–9.8) at elevated temperature (175–250 °C). We demonstrate that, in contrast to the behavior of Th at low temperature, the stability of Th-carbonate complexes is not sufficient for them to predominate at these conditions. Instead, the solubility of Th is governed by hydrolysis reactions. Under the experimental conditions investigated, the predominant hydroxyl complexes are Th(OH) 4 0 and Th(OH) 5 -. Furthermore, thermodynamic formation constants were derived for these species at the temperatures considered in our experiments (log β 4 = 43.34 and 44.31 at 175 and 200 °C, respectively, and log β 5 = 46.15 and 47.9 at 225 and 250 °C, respectively) to permit forward modeling of Th mobility in natural systems. Our study indicates that carbonate ions are unlikely to play a role in transporting Th in hydrothermal fluids. Summarizing the results of this study and our previous studies of the solubility of Th in hydrothermal fluids, we conclude that SO 4 2- is the primary ligand responsible for the hydrothermal transport of Th.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Proxy-based Bayesian inversion of strain tensor data measured during well tests

Recent instrument developments have made it possible to measure the strain tensor caused by injecting or pumping fluid from aquifers or reservoirs, but the full value of these data is limited because the long runtimes of poroelastic forward models makes it impractical to use many inversion schemes. This limits the interpretation of strain data for managing the recovery of resources or storage of wastes in the subsurface. This paper describes a method of inverting deformation data using a poroelastic numerical simulator so the results can be used to manage reservoirs or aquifers. We developed a workflow designed to reduce the number of simulations sufficiently to make it feasible to use DREAMzs, an advanced Bayesian inversion method that translates the uncertainties from different sources into unbiased posterior parameter distributions and uncertainty envelopes around the field data. Using a KNN proxy model for the poroelastic simulator is key to reducing the overall computations, and the workflow includes a strategy for ensuring the proxy model results converge on the results from the simulator. The workflow is tested using an idealized example that verifies the ability to correctly identify parameters and characterize noise used to perturb the data. Field data from an injection test at an oil reservoir near Tulsa, Oklahoma, are also used to evaluate the efficacy of the workflow with a real dataset. The workflow identified 265 history matching solutions out of 1240 total simulation runs (21% acceptance ratio), where the results were used to characterize posterior parameter distribution and evaluate the prediction uncertainty. Furthermore, this workflow is significant because it enables strain tensor, or other geomechanical measurements to be interpreted to guide decision-making during energy and environmental processes in the subsurface.

42 ENGINEERING↗
DOE OSTI2026

In situ quantification of fracture slip induced by hydraulic injections in a deep borehole: A comparison of two different borehole techniques

In situ measurements of fracture deformation during fluid injection are rare, yet essential for understanding the mechanical response of fractured rock. In this study, we evaluate the reliability of two methods by comparing their slip vector estimates: high-resolution borehole acoustic televiewer images captured before and after injection tests, and displacement data from a three-component borehole deformation probe recorded during the injections. Acoustic televiewer images capture only final in-plane displacement, whereas three-component borehole deformation measurements provide full 3D, transient fracture movement. Four injection tests in a fractured granitic rock mass along an inclined borehole at the Bedretto Underground Laboratory (Switzerland), beneath more than 1100 m of overburden, were analyzed. The two methods yielded consistent kinematics and comparable slip magnitudes, typically in the range of 0.2 – 0.6 mm. Angular differences between estimated slip directions ranged from 0.2° to 28°. For two test intervals with pre-existing fractures optimally oriented for slip, the measured slip vectors agreed with forward-modeled slip vectors based on the regional in situ stress field, supporting the Wallace–Bott hypothesis. Measurable shear displacement was identified in two test intervals containing fractures oriented nearly perpendicular to the least principal stress, including a newly formed hydraulic fracture. This demonstrates that even fractures unfavorably oriented for slip can undergo shear due to local stress perturbations around an inclined borehole.

Acoustic televiewer↗
DOE OSTI2021

Joint Bayesian Inference for Near-Surface Explosion Yield and Height-of-Burst

Forensic capabilities to understand chemical and nuclear explosions are greatly aided by an accurate estimate of explosive yield with uncertainty. The relationship between explosive size and geophysical observations of seismic, acoustic, and optical waves can be exploited to provide an estimate of yield. Any near-surface yield estimate is complicated by the surface interaction, so an estimate for the explosion height-of-burst is necessarily included in the relationship. Additionally, the relationship dictates a trade-off between estimates of yield and height-of-burst. Fortunately, the surface interaction for each type of observation is different, which breaks the trade-off, and the inclusion of height-of-burst with multiple data types improves yield estimation. We define simple parametric forward models to relate seismoäcoustoöptic observations from a data set of known explosive yields and height-of-bursts. The parameters of the models and a prediction for the yield and height-of-burst of a new event can then be estimated given new observations via Bayesian inference. We report posterior distribution estimates of the parametric models using a Markov chain Monte Carlo sampling technique. These models are then used to predict the yield and height-of-burst of SUGAR, a historical near-surface nuclear explosion, using its reported historical observations. The reported yield of 1.2 ktonne Trinitrotoluene (TNT)-equivalent (Department of Energy, 2015) is within the estimated posterior. Yield uncertainty can be estimated from the spread of the posterior, which is between 0.9 and 2.1 ktonne TNT-equivalent. The posterior for height-of-burst has a wider range between 10 m below and 8 m above ground that includes the true height-of-burst of 1 m.

58 GEOSCIENCES↗
DOE OSTI2023

Time‐Dependent Weakening of Granite at Hydrothermal Conditions

Abstract The evolution of a fault's frictional strength during the interseismic period is a critical component of the earthquake cycle, yet there have been relatively few studies that examine the time‐dependent evolution of strength at conditions representative of seismogenic depths. Using a simulated fault in Westerly granite, we examined how frictional strength evolves under hydrothermal conditions up to 250°C during slide‐hold‐slide experiments. At temperatures ≤100°C, frictional strength generally increases with hold duration but, at 200 and 250°C, an initial increase in strength transitions to rapid time‐dependent weakening for holds longer than 14 hr. Forward modeling of long hold periods at 250°C using the rate and state friction constitutive equations requires a second, strongly negative, state variable with a long evolution distance. This implies that significant hydrothermal alteration is occurring at 250°C, consistent with microstructural observations of dissolution and secondary mineral precipitation.

Geology↗