Abstract The rheology of the upper mantle impacts a variety of geodynamic processes, including postseismic deformation following great earthquakes and post-glacial rebound. The deformation of upper mantle rocks is controlled by the rheology of olivine, the most abundant upper mantle mineral. The mechanical properties of olivine at steady state are well constrained. However, the physical mechanism underlying transient creep, an evolutionary, hardening phase converging to steady state asymptotically, is still poorly understood. Here, we constrain a constitutive framework that captures transient creep and steady state creep consistently using the mechanical data from laboratory experiments on natural dunites containing at least 94% olivine under both hydrous and anhydrous conditions. The constitutive framework represents a Burgers assembly with a thermally activated nonlinear stress-versus-strain-rate relationship for the dashpots. Work hardening is obtained by the evolution of a state variable that represents internal stress. We determine the flow law parameters for dunites using a Markov chain Monte Carlo method. We find the activation energy $$430\pm 20$$ 430 ± 20 and $$250\pm 10$$ 250 ± 10 kJ/mol for dry and wet conditions, respectively, and the stress exponent $$2.0\pm 0.1$$ 2.0 ± 0.1 for both the dry and wet cases for transient creep, consistently lower than those of steady-state creep, suggesting a separate physical mechanism. For wet dunites in the grain-boundary sliding regime, the grain-size dependence is similar for transient creep and steady-state creep. The lower activation energy of transient creep could be due to a higher jog density of the corresponding soft-slip system. More experimental data are required to estimate the activation volume and water content exponent of transient creep. The constitutive relation used and its associated flow law parameters provide useful constraints for geodynamics applications. Graphical Abstract
more »
« less
Modeling and Quantifying Parameter Uncertainty of Co‐Seismic Non‐Classical Nonlinearity in Rocks
Dynamic perturbations reveal unconventional nonlinear behavior in rocks, as evidenced by field and laboratory studies. During the passage of seismic waves, rocks exhibit a decrease in elastic moduli, slowly recovering after. Yet, comprehensive physical models describing these moduli alterations remain sparse and insufficiently validated against observations. Here, we demonstrate the applicability of two physical damage models—the internal variable model (IVM) and the continuum damage model (CDM)—to provide quantitative descriptions of nonlinear co‐seismic elastic wave propagation observations. While the IVM uses one internal variable to describe the evolution of elastic material moduli, the CDM damage variable is a mathematical representation of microscopic defects. We recast the IVM and CDM models as nonlinear hyperbolic partial differential equations and implement 1D and 2D numerical simulations using an arbitrary high‐order discontinuous Galerkin method. We verify the modeling results with co‐propagating acousto‐elastic experimental measurements. Subsequently, we infer the parameters for these nonlinear models from laboratory experiments using probabilistic Bayesian inversion and 2D simulations. By adopting the Adaptive Metropolis Markov chain Monte Carlo method, we quantify the uncertainties of inferred parameters for both physical models, investigating their interplay in 70,000 simulations. We find that the damage variables can trade off with the stress‐strain nonlinearity in discernible ways. We discuss physical interpretations of both damage models and that our CDM quantitatively captures an observed damage increase with perturbation frequency. Our results contribute to a more holistic understanding of co‐seismic damage and post‐seismic recovery after earthquakes bridging the worlds of theoretical analysis and laboratory findings.
more » « less- NSF-PAR ID:
- 10484651
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 129
- Issue:
- 1
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)Loss of operation or devastating damage to buildings and industrial structures, as well as equipment housed in them, has been observed due to earthquake-induced vibrations. A common source of operational downtime is due to the performance reduction of vital equipment, which are sensitive to the total transmitted acceleration. A well-known method of protecting such equipment is seismic isolation of the equipment itself (or a group of equipment), as opposed to the entire structure due to the lower cost of implementation. The first objective of this dissertation is assessing a rolling isolation system (RIS) based on existing design guidelines for telecommunications equipment. A discrepancy is observed between the required response spectrum (RRS) and the one and only accelerogram recommended in the guideline. Several filters are developed to generate synthetic accelerograms that are compatible with the RRS. The generated accelerograms are used for probabilistic assessment of a RIS that is acceptable per the guideline. This assessment reveals large failure probability due to displacement demands in excess of the displacement capacity of the RIS. When the displacement demands on an isolation system are in excess of its capacity, impacts result in spikes in transmitted acceleration. Therefore, the second objective of this dissertation is to design impact prevention/mitigation mechanisms. A dual-mode system is proposed where the behavior changes when the displacement exceeds a predefined threshold. A new piecewise optimal control approach is developed and applied to find the best possible mechanism for the region beyond the threshold. By utilizing the designed curves obtained from the proposed optimal control procedure, a Kelvin-Voigt device is tuned for illustrative purposes. On the other hand, the preference for protecting equipment decreases as the earthquake intensity increases. In extreme seismic loading, the response mitigation of the primary structure (i.e., life safety and collapse prevention) is of greater concern than protecting isolated equipment. Therefore, the third objective of this dissertation is to develop an innovative dual-mode system that can behave as equipment isolation under low to moderate seismic loading and passively transition to behave as a vibration absorber for the primary structure under extreme seismic loading. To reduce the computational cost of simulating a large linear elastic structure with nonlinear attachments (i.e., equipment isolation with cubic hardening nonlinearity), a reduced order modeling method is introduced that can capture the behavior of such nonlinear coupled systems. The method is applied to study the feasibility of dual-mode vibration isolation/absorber. To this end, nonlinear transmissibility curves for the roof displacement and isolated mass total acceleration are developed from the steady-state responses of dual-mode systems using the harmonic balanced method. The final objective of this dissertation is to extend the reduced order modeling method developed for linear elastic structure with nonlinear attachment to inelastic structures (without attachments). The new inelastic model condensation (IMC) method uses the modal properties of the full structural model (in the elastic range) to construct a linear reduced order model in conjunction with a hysteresis model to capture the hysteretic inter-story restoring forces. The parameters of these hysteretic forces are easily tuned, in order to fit the inelastic behavior of the condensed structure to that of the full model under a variety of simple loading scenarios. The fidelity of structural models condensed in this way is demonstrated via simulation for different ground motion intensities on three different building structures with various heights. The simplicity, accuracy, and efficiency of this approach could significantly alleviate the computational burden of performance-based earthquake engineering.more » « less
-
null (Ed.)Earthquake prediction is the holy grail of seismology. Many previous studies have searched for robust precursory signals to inform us of imminent earthquakes, the most significant of which are seen in laboratory experiments as temporal changes in pressure and shear wave velocities during the seismic cycle. Similar changes are seen in natural faults and the surrounding structurally complex network of fractures with nested hierarchy of localized deformation, referred to as fault damage zone. However, little is known whether such temporal changes in material properties contains any precursory signals for imminent earthquakes.Conversely, the effect of precursory velocity changes on the seismic cycle is not well understood. By imposing shear wave velocity changes in fault damage zones, we investigate the effects of these precursors on multiple stages of the seismic cycle, including nucleation, coseismic, postseismic, and interseismic stages. We perform 2D fully dynamic earthquake cycle simulations with a fault-parallel damage zone for strike-slip fault systems with antiplane geometry. The fault is governed by rate-state-dependent friction laws, and the fault damage zone material is considered elastic. Our preliminary results show that the temporal onset of shear wave velocity drop causes a reduction in earthquake recurrence intervals over the seismic cycle. Furthermore, a dynamic earthquake rupture within the seismic cycle terminates much faster and abruptly in models with precursory velocity changes. We will also discuss how the precursory velocity changes affect the fault-slip behavior, including fast-slip, slow-slip, and aseismic creep, for different amplitudes of shear wave velocity changes at different compliance contrast of the fault damage zones. Our results highlight the importance of short and long-term monitoring of fault zone structures for better assessment of regional seismic hazard.more » « less
-
more » « less
ABSTRACT: Due to rock mass being commonly subjected to compressive or shear loading, the mode II fracture toughness is an important material parameter for rocks. Fracturing in rocks is governed by the behavior of a nonlinear region surrounding the crack tip called the fracture process zone (FPZ). However, the characteristics of mode II fracture are still determined based on the linear elastic fracture mechanics (LEFM), which assumes that a pure mode II loading results in a pure mode II fracture. In this study, the FPZ development in Barre granite specimens under mode II loading was investigated using the short beam compression (SBC) test. Additionally, the influence of lateral confinement on various characteristics of mode II fracture was studied. The experimental setup included the simultaneous monitoring of surface deformation using the two-dimensional digital image correlation technique (2D-DIC) to identify fracture mode and characterize the FPZ evolution in Barre granite specimens. The 2D-DIC analysis showed a dominant mixed-mode I/II fracture in the ligament between two notches, irrespective of confinement level on the SBC specimens. The influence of confinement on the SBC specimens was assessed by analyzing the evolution of crack displacement and changes in value of mode II fracture toughness. Larger levels of damage in confined specimens were observed prior to the failure than the unconfined specimens, indicating an increase in the fracture resistance and therefore mode II fracture toughness with the confining stress.
1. INTRODUCTION The fracturing in laboratory-scale rock specimens is often characterized by the deformation of the inelastic region surrounding the crack tips, also known as the fracture process zone (FPZ) (Backers et al., 2005; Ghamgosar and Erarslan, 2016). While the influence of the FPZ on mode I fracture in rocks has been extensively investigated, there are limited studies on FPZ development in rocks under pure mode II loading (Ji et al., 2016; Lin et al., 2020; Garg et al., 2021; Li et al., 2021).
-
Toward a Self‐Consistent Characterization of Lithospheric Plates Using Full‐Spectrum Viscoelasticitymore » « less
Determining the thickness of the lithosphere in any given setting combines uncertainty in both the observational method and laboratory‐derived understanding of mantle rheology. The many observational and modeling criteria across geophysical subfields for plate thickness lead to significant differences in plate thickness estimates depending on the process of interest, be it seismic wave propagation or relaxation in response to changes in loads—from earthquakes, ice sheets to volcanoes—or convection. This paper proposes a framework in which to model and interpret upper mantle mechanical structure smoothly across the full spectrum of geophysical timescales. We integrate viscous, elastic, and linear anelastic constitutive models and calculate the mechanical response from convective to seismic wave timescales (i.e., 0 to infinite frequency or, in practice,
10−15 to 1 Hz). We apply these calculations to 1‐D thermal structures and determine the normalized complex viscosity, a quantity that shows clearly the role of transient creep in weakening rock relative to the associated Maxwell rheology. Using various criteria for the lithosphere‐asthenosphere boundary, we show that the apparent plate thickness will be thicker at higher frequencies than at lower frequencies. Additional calculations for nonlinear Maxwell behavior (dislocation mechanisms) demonstrate significant changes in the apparent plate structure, decreasing the long‐term plate thickness, consistent with observations. Other effects such as dislocation damping (associated with a steady‐state dislocation structure), melt, water, major element composition, and grain size are not included here but, when incorporated into this framework, will significantly change the full‐spectrum plate thickness predictions.
