Person: Rice, James
Loading...
Email Address
AA Acceptance Date
Birth Date
Research Projects
Organizational Units
Job Title
Last Name
Rice
First Name
James
Name
Rice, James
38 results
Search Results
Now showing 1 - 10 of 38
Publication Heating, weakening and shear localization in earthquake rupture(The Royal Society, 2017-08-21) Rice, JamesField and borehole observations of active earthquake fault zones show that shear is often localized to principal deforming zones of order 0.1-10 mm width. This paper addresses how frictional heating in rapid slip weakens faults dramatically, relative to their static frictional strength, and promotes such intense localization. Pronounced weakening occurs even on dry rock-on-rock surfaces, due to flash heating effects, at slip rates above approximately 0.1 m s-1 (earthquake slip rates are typically of the order of 1 m s-1). But weakening in rapid shear is also predicted theoretically in thick fault gouge in the presence of fluids (whether native ground fluids or volatiles such as H2O or CO2 released by thermal decomposition reactions), and the predicted localizations are compatible with such narrow shear zones as have been observed. The underlying concepts show how fault zone materials with high static friction coefficients, approximately 0.6-0.8, can undergo strongly localized shear at effective dynamic friction coefficients of the order of 0.1, thus fitting observational constraints, e.g. of earthquakes producing negligible surface heat outflow and, for shallow events, only rarely creating extensive melt. The results to be summarized include those of collaborative research published with Nicolas Brantut (University College London), Eric Dunham (Stanford University), Nadia Lapusta (Caltech), Hiroyuki Noda (JAMSTEC, Japan), John D. Platt (Carnegie Institution for Science, now at *gramLabs), Alan Rempel (Oregon State University) and John W. Rudnicki (Northwestern University). This article is part of the themed issue 'Faulting, friction and weakening: from slow to fast motion.Publication Finite Element Modeling of Branched Ruptures Including Off-Fault Plasticity(Seismological Society of America (SSA), 2012) DeDontney, Nora Lee; Rice, James; Dmowska, RenataFault intersections are a geometric complexity that frequently occurs in nature. Here we focus on earthquake rupture behavior when a continuous planar main fault has a second fault branching off of it. We use the finite element (FE) method to examine which faults are activated and how the surrounding material responds for both elastic and elastic–plastic off-fault descriptions. Compared to an elastic model, a noncohesive elastic–plastic material, intended to account for zones of damaged rock bordering maturely slipped faults, will inhibit rupture on compressional side branches and promote rupture of extensional side branches. Activation of extensional side branches can be delayed and is triggered by continued rupture propagation on the main fault. We examine the deformation near the branching junction and find that fault opening is common for elastic materials, especially for compressional side branches. An elastic–plastic material is more realistic because elevated stresses around the propagating rupture tip and at the branching junction should bring the surrounding material to failure. With an elastic–plastic material model, fault opening is inhibited for a range of realistic material parameters. For large cohesive strengths, opening can occur, but with material softening, a real feature of plastically deforming rocks, open- ing can be prevented. We also discuss algorithmic artifacts that may arise due to the presence of such a triple junction. When opening does not occur, the behavior at the triple junction is simplified and standard contact routines in FE programs are able to properly represent the physical situation.Publication Deformation-Induced Melting in the Margins of the West Antarctic Ice Streams(Wiley-Blackwell, 2014) Suckale, Jenny; Platt, John; Perol, Thibaut; Rice, JamesFlow of glacial ice in the West Antarctic Ice Sheet localizes in narrow bands of fast-flowing ice streams bordered by ridges of nearly stagnant ice, but our understanding of the physical processes that generate this morphology is incomplete. Here we study the thermal and mechanical properties of ice-stream margins, where flow transitions from rapid to stagnant over a few kilometers. Our goal is to explore under which conditions the intense shear deformation in the margin may lead to deformation-induced melting. We propose a 2-D model that represents a cross section through the ice stream margin perpendicular to the downstream flow direction. We limit temperature to the melting point to estimate melt rates based on latent heat. Using rheology parameters as constrained by laboratory data and observations, we conclude that a zone of temperate ice is likely to form in active shear margins.Publication Stability and Localization of Rapid Shear in Fluid-Saturated Fault Gouge: 2. Localized Zone Width and Strength Evolution(Wiley-Blackwell, 2014) Platt, John; Rudnicki, John W.; Rice, JamesField and laboratory observations indicate that at seismic slip rates most shearing is confined to a very narrow zone, just a few tens to hundreds of microns wide, and sometimes as small as a few microns. Rice et al. (2014) analyzed the stability of uniform shear in a fluid-saturated gouge material. They considered two distinct mechanisms to limit localization to a finite thickness zone, rate-strengthening friction, and dilatancy. In this paper we use numerical simulations to extend beyond the linearized perturbation context in Rice et al. (2014), and study the behavior after the loss of stability. Neglecting dilatancy we find that straining localizes to a width that is almost independent of the gouge layer width, suggesting that the localized zone width is set by the physical properties of the gouge material. Choosing parameters thought to be representative of a crustal depth of 7 km, this predicts that deformation should be confined to a zone between 4 and 44 μm wide. Next, considering dilatancy alone we again find a localized zone thickness that is independent of gouge layer thickness. For dilatancy alone we predict localized zone thicknesses between 1 and 2 μm wide for a depth of 7 km. Finally, we study the impact of localization on the shear strength and temperature evolution of the gouge material. Strain rate localization focuses frictional heating into a narrower zone, leading to a much faster temperature rise than that predicted when localization is not accounted for. Since the dynamic weakening mechanism considered here is thermally driven, this leads to accelerated dynamic weakening.Publication Stability and Localization of Rapid Shear in Fluid-Saturated Fault Gouge: 1. Linearized Stability Analysis(Wiley-Blackwell, 2014) Rice, James; Rudnicki, John W.; Platt, JohnField observations of major earthquake fault zones show that shear deformation is often confined to principal slipping zones that may be of order 1–100 μm wide, located within a broader gouge layer of order 10–100 mm wide. This paper examines the possibility that the extreme strain localization observed may be due to the coupling of shear heating, thermal pressurization, and diffusion. In the absence of a stabilizing mechanism shear deformation in a continuum analysis will collapse to an infinitesimally thin zone. Two possible stabilizing mechanisms, studied in this paper, are rate-strengthening friction and dilatancy. For rate-strengthening friction alone, a linear stability analysis shows that uniform shear of a gouge layer is unstable for perturbations exceeding a critical wavelength. Using this critical wavelength we predict a width for the localized zone as a function of the gouge properties. Taking representative parameters for fault gouge at typical centroidal depths of crustal seismogenic zones, we predict localized zones of order 5–40 μm wide, roughly consistent with field and experimental observations. For dilatancy alone, linearized strain rate perturbations with a sufficiently large wavelength will undergo transient exponential growth before decaying back to uniform shear. The total perturbation strain accumulated during this transient strain rate localization is shown to be largely controlled by a single dimensionless parameter E, which is a measure of the dilatancy of the gouge material due to an increase in strain rate.Publication Tsunami Wave Analysis and Possibility of Splay Fault Rupture During the 2004 Indian Ocean Earthquake(Springer, 2012) DeDontney, Nora; Rice, JamesThe 2004 Indian Ocean tsunami was observed by two satellites, close in space and time, that traversed the Indian ocean two hours after the Sumatra-Andaman earthquake, but which observed different tsunami lead wave morphologies. The earlier satellite, Jason-1, recorded a lead wave with two peaks of similar amplitude and wavelength, while the later satellite, TOPEX/Poseidon, recorded a lead wave with only one longer wavelength uplift. To resolve this disparity, we examine the travel paths of long wavelength waves over the seafloor bathymetry. Waves traveling from the margin will traverse significantly different paths to arrive at the two satellite transects. The result is that the satellites are sensitive to different parts of the margin; Jason-1 is highly sensitive to the margin in the area of the epicenter, while TOPEX is sensitive to a more northerly section. By developing solutions of the ocean gravity wave equations, accounting for dispersion, we show that the double peak of the Jason-1 satellite observations are consistent with coseismic rupture of a splay fault of limited along-strike extent, located north of Simeulue Island. The doubly peaked morphology can be reproduced with co-activation of the subduction zone interface and the splay fault, which creates a seafloor uplift pattern with two distinct areas of uplift. The Jason-1 satellite is sensitive to a splay fault in this portion of the margin, whereas the TOPEX satellite would not be significantly affected by this uplift pattern. By back-projecting satellite observation points to the margin, we constrain the location of the proposed splay fault and find that it correlates with a bathymetric high. The aftershock locations, uplift of corals on Simeulue Island and a fault scarp on Pulau Salaut Besar are also consistent with the activation of a splay fault in the area delimited by the back-projection. Our work also shows that it is critical to fully capture gravity wave dispersion in order to represent features of the lead wave profile that may not be as well characterized by the shallow water (long-wavelength) model. It is also necessary to account for dispersion so as to precisely assess wavefront travel times; this leads us to conclude that the rupture must have reached very near to the trench and propagated with an updip rupture velocity of order 2.0 km/s or more.Publication Finite Element Model of Branched Ruptures Including Off-Fault Plasticity(Seismological Society of America, 2012) DeDontney, Nora; Rice, James; Dmowska, RenataFault intersections are a geometric complexity that frequently occurs in nature. Here we focus on earthquake rupture behavior when a continuous, planar main fault has a second fault branching off of it. We use the finite element method to examine which faults are activated and how the surrounding material responds for both elastic and elastic-plastic off-fault descriptions. Compared to an elastic model, a non-cohesive, elastic-plastic material, intended to account for zones of damaged rock bordering maturely slipped faults, will inhibit rupture on compressional side branches and promote rupture of extensional side branches. Activation of extensional side branches can be delayed and is triggered by continued rupture propagation on the main fault. We examine the deformation near the branching junction and find that fault opening is common for elastic materials, especially for compressional side branches. An elastic-plastic material is more realistic since elevated stresses around the propagating rupture tip and at the branching junction should bring the surrounding material to failure. With an elastic-plastic material model, fault opening is inhibited for a range of realistic material parameters. For large cohesive strengths opening can occur, but with material softening, a real feature of plastically deforming rocks, opening can be prevented. We also discuss algorithmic artifacts that may arise due to the presence of such a triple junction. When opening does not occur, the behavior at the triple junction is simplified and standard contact routines in finite element programs are able to properly represent the physical situation.Publication Nucleation of Slip-Weakening Rupture Instability in Landslides by Localized Increase of Pore Pressure(American Geophysical Union, 2012) Viesca, Robert C.; Rice, JamesWe model landslide initiation as slip surface growth driven by local elevated pore pressure, with particular reference to submarine slides. Assuming an elastic medium and friction that weakens with slip, solutions exist in which the slip surface may dynamically grow, without further pore pressure increases, at a rate of the order of the sediment shear wave speed, a situation comparable to earthquake nucleation. The size of the rupture at this transition point depends weakly on the imposed pore pressure pro file; however, the amount of slip at the transition depends strongly on whether the pore pressure was broadly or sharply elevated. Sharper profiles may result in pore pressures reaching the total slope-normal stress before dynamic rupture is nucleated. While we do not account for modes of failure other than pure slip on a failure surface, this may be an indication that additional modes involving liquefaction or hydraulic cracking may be factors in the initiation of shallow slope failure. We identify two lengthscales, one geometrical (h, depth below the free surface) and one material (l, determined by the frictional weakening rate) and a transition in nucleation behavior between effectively "deep" and "shallow" limits dependent on their ratio. Whether dynamic propagation of failure is indefinite or arresting depends largely on whether the background shear stress is closer to nominal peak or residual frictional strength. This is determined in part by background pore pressures, and to consider the submarine case we simplify a common sedimentation/consolidation approach to reflect interest in near-seafloor conditions.Publication Modeling Turbulent Hydraulic Fracture Near a Free Surface(American Society of Mechanical Engineers, 2012) Tsai, Victor C.; Rice, JamesMotivated by observations of subglacial drainage of water, we consider a hydraulic fracture problem in which the crack fluid flow. Using a hybrid Chebyshev/series-minimization numerical approach, we solve for the pressure profile, crack opening displacement, and crack growth rate for a crack that begins relatively short but eventually becomes long compared with the distance to the free surface. We plot non-dimensionalized results for a variety of different times, corresponding with different fracture lengths, and find substantial differences when free-surface effects are important.Publication Earthquake Ruptures with Thermal Weakening and the Operation of Major Faults at Low Overall Stress Levels(American Geophysical Union, 2009) Noda, Hiroyuki; Dunham, Eric M.; Rice, JamesWe model ruptures on faults that weaken in response to flash heating of microscopic asperity contacts (within a rate-and-state framework) and thermal pressurization of pore fluid. These are arguably the primary weakening mechanisms on mature faults at coseismic slip rates, at least prior to large slip accumulation. Ruptures on strongly rate-weakening faults take the form of slip pulses or cracks, depending on the background stress. Self-sustaining slip pulses exist within a narrow range of stresses: below this range, artificially nucleated ruptures arrest; above this range, ruptures are crack-like. Natural quakes will occur as slip pulses if faults operate at the minimum stress required for propagation. Using laboratory-based flash heating parameters, propagation is permitted when the ratio of shear to effective normal stress on the fault is 0.2–0.3; this is mildly influenced by reasonable choices of hydrothermal properties. The San Andreas and other major faults are thought to operate at such stress levels. While the overall stress level is quite small, the peak stress at the rupture front is consistent with static friction coefficients of 0.6–0.9. Growing slip pulses have stress drops of ~3 MPa; slip and the length of the slip pulse increase linearly with propagation distance at ~0.14 and ~30 m/km, respectively. These values are consistent with seismic and geologic observations. In contrast, cracks on faults of the same rheology have stress drops exceeding 20 MPa, and slip at the hypocenter increases with distance at ~1 m/km.