Parallel 3D Simulation of a Fault Gouge using the Lattice Solid Model

Despite the insight gained from 2-D particle models, and given that the dynamics of crustal faults occur in 3-D space, the question remains, how do the 3-D fault gouge dynamics differ from those in 2-D? Traditionally, 2-D modeling has been preferred over 3-D simulations because of the computational cost of solving 3-D problems. However, modern high performance computing architectures, combined with a parallel implementation of the Lattice Solid Model (LSM), provide the opportunity to explore 3-D fault micro-mechanics and to advance understanding of effective constitutive relations of fault gouge layers. In this paper, macroscopic friction values from 2-D and 3-D LSM simulations, performed on an SGI Altix 3700 super-cluster, are compared. Two rectangular elastic blocks of bonded particles, with a rough fault plane and separated by a region of randomly sized non-bonded gouge particles, are sheared in opposite directions by normally-loaded driving plates. The results demonstrate that the gouge particles in the 3-D models undergo significant out-of-plane motion during shear. The 3-D models also exhibit a higher mean macroscopic friction than the 2-D models for varying values of interparticle friction. 2-D LSM gouge models have previously been shown to exhibit accelerating energy release in simulated earthquake cycles, supporting the Critical Point hypothesis. The 3-D models are shown to also display accelerating energy release, and good fits of power law time-to-failure functions to the cumulative energy release are obtained.

On the possibility of elastic strain localisation in a fault

The phenomenon of strain localisation is often observed in shear deformation of particulate materials, e.g., fault gouge. This phenomenon is usually attributed to special types of plastic behaviour of the material (e.g., strain softening or mismatch between dilatancy and pressure sensitivity or both). Observations of strain localisation in situ or in experiments are usually based on displacement measurements and subsequent computation of the displacement gradient. While in conventional continua the symmetric part of the displacement gradient is equal to the strain, it is no longer the case in the more realistic descriptions within the framework of generalised continua. In such models the rotations of the gouge particles are considered as independent degrees of freedom the values of which usually differ from the rotation of an infinitesimal volume element of the continuum, the latter being described for infinitesimal deformations by the non-symmetric part of the displacement gradient. As a model for gouge material we propose a continuum description for an assembly of spherical particles of equal radius in which the particle rotation is treated as an independent degree of freedom. Based on this model we consider simple shear deformations of the fault gouge. We show that there exist values of the model parameters for which the displacement gradient exhibits a pronounced localisation at the mid-layers of the fault, even in the absence of inelasticity. Inelastic effects are neglected in order to highlight the role of the independent rotations and the associated additional parameters. The localisation-like behaviour occurs if (a) the particle rotations on the boundary of the shear layer are constrained (this type of boundary condition does not exist in a standard continuum) and (b) the contact moment-or bending stiffness is much smaller than the product of the effective shear modulus of the granulate and the square of the width of the gouge layer. It should be noted however that the virtual work functional is positive definite over the range of physically meaningful parameters (here: contact stiffnesses, solid volume fraction and coordination number) so that strictly speaking we are not dealing with a material instability.

Finite deformation model of simple shear of fault with microrotations: apparent strain localisation and en-echelon fracture pattern

Strain localisation is a widespread phenomenon often observed in shear and compressive loading of geomaterials, for example, the fault gouge. It is believed that the main mechanisms of strain localisation are strain softening and mismatch between dilatancy and pressure sensitivity. Observations show that gouge deformation is accompanied by considerable rotations of grains. In our previous work as a model for gouge material, we proposed a continuum description for an assembly of particles of equal radius in which the particle rotation is treated as an independent degree of freedom. We showed that there exist critical values of the model parameters for which the displacement gradient exhibits a pronounced localisation at the mid-surface layers of the fault, even in the absence of inelasticity. Here, we generalise the model to the case of finite deformations characteristic for the gouge deformation. We derive objective constitutive relationships relating the Jaumann rates of stress and moment stress to the relative strain and curvature rates, respectively. The model suggests that the pattern of localisation remains the same as in the linear case. However, the presence of the Jaumann terms leads to the emergence of non-zero normal stresses acting along and perpendicular to the shear layer (with zero hydrostatic pressure), and localised along the mid-line of the gouge; these stress components are absent in the linear model of simple shear. These additional normal stresses, albeit small, cause a change in the direction in which the maximal normal stresses act and in which en-echelon fracturing is formed.

Evidence of episodic brittle faulting in the cratonic part of the Peninsular India and its implications for seismic hazard in slow deforming regions

The region around Waclakkancheri, in the province of Kerala, India, which lies in the vicinity of Palghat-Cauvery ;hear zone (within the Precambrian crystalline terrain), has been a site of microseismic activity since 1989. Earlier studies had identified a prominent WNW-ESE structure overprinting on the E-W trending lineaments associated with Palghat-Cauvery shear zone. We have mapped this structure, located in a chamockite quarry near Desamangalam, Waclakkancheri, which we identify as a ca. 30 km-long south dipping reverse fault. This article presents the characteristics of this fault zone exposed on the exhumed crystalline basement and discusses its significance in understanding the earthquake potential of the region. This brittle deformation zone consists of fracture sets with small-scale displacement and slip planes with embedded fault gouges. The macroscopic as well as the microscopic studies of this fault zone indicate that it evolved through different episodes of faulting in the presence of fluids. The distinct zones within consolidated gouge and the cross cutting relationship of fractures indicate episodic fault activity. At least four faulting episodes can be recognized based on the sequential development of different structural elements in the fault rocks. The repeated ruptures are evident along this shear zone and the cyclic behavior of this fault consists of co-seismic ruptures alternating with inter-seismic periods, which is characterized by the sealed fractures and consolidated gouge. The fault zone shows a minimum accumulated dip/oblique slip of 2.1 m in the reverse direction with a possible characteristic slip of 52 cm (for each event). The ESR dating of fault gouge indicates that the deformation zone records a major event in the Middle Quaternary. The empirical relationships between fault length and slip show that this fault may generate events M >= 6. The above factors suggest that this fault may be characterized as potentially active. Our study offers some new pointers that can be used in other slow deforming cratonic hinterlands in exploring the discrete active faults.

Structural and clumped-isotope constraints on the mechanisms of displacement along low-angle detachments

Despite years of research on low-angle detachments, much about them remains enigmatic. This thesis addresses some of the uncertainty regarding two particular detachments, the Mormon Peak detachment in Nevada and the Heart Mountain detachment in Wyoming and Montana.

Constraints on the geometry and kinematics of emplacement of the Mormon Peak detachment are provided by detailed geologic mapping of the Meadow Valley Mountains, along with an analysis of structural data within the allochthon in the Mormon Mountains. Identifiable structures well suited to constrain the kinematics of the detachment include a newly mapped, Sevier-age monoclinal flexure in the hanging wall of the detachment. This flexure, including the syncline at its base and the anticline at its top, can be readily matched to the base and top of the frontal Sevier thrust ramp, which is exposed in the footwall of the detachment to the east in the Mormon Mountains and Tule Springs Hills. The ~12 km of offset of these structural markers precludes the radial sliding hypothesis for emplacement of the allochthon.

The role of fluids in the slip along faults is a widely investigated topic, but the use of carbonate clumped-isotope thermometry to investigate these fluids is new. Faults rocks from within ~1 m of the Mormon Peak detachment, including veins, breccias, gouges, and host rocks, were analyzed for carbon, oxygen, and clumped-isotope measurements. The data indicate that much of the carbonate breccia and gouge material along the detachment is comminuted host rock, as expected. Measurements in vein material indicate that the fluid system is dominated by meteoric water, whose temperature indicates circulation to substantial depths (c. 4 km) in the upper crust near the fault zone.

Slip along the subhorizontal Heart Mountain detachment is particularly enigmatic, and many different mechanisms for failure have been proposed, predominantly involving catastrophic failure. Textural evidence of multiple slip events is abundant, and include multiple brecciation events and cross-cutting clastic dikes. Footwall deformation is observed in numerous exposures of the detachment. Stylolitic surfaces and alteration textures within and around “banded grains” previously interpreted to be an indicator of high-temperature fluidization along the fault suggest their formation instead via low-temperature dissolution and alteration processes. There is abundant textural evidence of the significant role of fluids along the detachment via pressure solution. The process of pressure solution creep may be responsible for enabling multiple slip events on the low-angle detachment, via a local rotation of the stress field.

Clumped-isotope thermometry of fault rocks associated with the Heart Mountain detachment indicates that despite its location on the flanks of a volcano that was active during slip, the majority of carbonate along the Heart Mountain detachment does not record significant heating above ambient temperatures (c. 40-70°C). Instead, cold meteoric fluids infiltrated the detachment breccia, and carbonate precipitated under ambient temperatures controlled by structural depth. Locally, fault gouge does preserve hot temperatures (>200°C), as is observed in both the Mormon Peak detachment and Heart Mountain detachment areas. Samples with very hot temperatures attributable to frictional shear heating are present but rare. They appear to be best preserved in hanging wall structures related to the detachment, rather than along the main detachment.

Evidence is presented for the prevalence of relatively cold, meteoric fluids along both shallow crustal detachments studied, and for protracted histories of slip along both detachments. Frictional heating is evident from both areas, but is a minor component of the preserved fault rock record. Pressure solution is evident, and might play a role in initiating slip on the Heart Mountain fault, and possibly other low-angle detachments.

Correlação entre contexto morfoestrutural e sismicidade nas regiões de João Câmara e São Rafael (RN)

This MSc thesis describes brittle deformation in two seismic zones located in north-eastern Brazil: João Câmara and São Rafael, Rio Grande do Norte State. Both areas show seismogenic faults, Samambaia and São Rafael, indicated by narrow zones of epicentres with a strike of 040o, a lenght of 30 km and 4 km, and a depth of 1-12 and 0,5-4 km, respectively. The first seismological and geological studies suggested blind faults or faults that were still in the beginning of the nucleation process. The region is under E-W-oriented compression and is underlain by Precambrian terrains, deformed by one or more orogenic cycles, which generated shear zones generally marked by strong pervasive foliation and sigmoidal shapes. The crystalline basement is capped by the Cretaceous Potiguar basin, which is also locally capped by Pliocene continental siliciclastic deposits (Barreiras Formation), and Quaternary alluvium. The main aim of this study was to map epicentral areas and find whether there are any surface geological or morphotectonic expression related to the seismogenic faults. A detailed geological map was carried out in both seismic areas in order to identify brittle structures and fault-related drainage/topographic features. Geological and morphotectonic evidence indicate that both seismogenic faults take place along dormant structures. They either cut Cenozoic rocks or show topographic expression, i.e., are related to topographic heights or depressions and straight river channels. Faults rocks in the Samambaia and São Rafael faults are cataclasite, fault breccia, fault gouge, pseudotachylyte, and quartz veins, which point to reactivation processes in different crustal levels. The age of the first Samambaia and the São Rafael faulting movement possibly ranges from late Precambrian to late Cretaceous. Both fault cut across Precambrian fabric. They also show evidence of brittle processes which took place between 4 and 12 km deep, which probably have not occurred in Cenozoic times. The findings are of great importance for regional seismic hazard. They indicate that fault zones are longer than previously suggested by seismogenic studies. According to the results, the methodology used during this thesis may also be useful in other neotectonic investigation in intraplate areas

Solute transport in crystalline rocks at Äspö — I: Geological basis and model calibration

Water-conducting faults and fractures were studied in the granite-hosted A¨ spo¨ Hard Rock Laboratory (SE Sweden). On a scale of decametres and larger, steeply dipping faults dominate and contain a variety of different fault rocks (mylonites, cataclasites, fault gouges). On a smaller scale, somewhat less regular fracture patterns were found. Conceptual models of the fault and fracture geometries and of the properties of rock types adjacent to fractures were derived and used as input for the modelling of in situ dipole tracer tests that were conducted in the framework of the Tracer Retention Understanding Experiment (TRUE-1) on a scale of metres. After the identification of all relevant transport and retardation processes, blind predictions of the breakthroughs of conservative to moderately sorbing tracers were calculated and then compared with the experimental data. This paper provides the geological basis and model calibration, while the predictive and inverse modelling work is the topic of the companion paper [J. Contam. Hydrol. 61 (2003) 175]. The TRUE-1 experimental volume is highly fractured and contains the same types of fault rocks and alterations as on the decametric scale. The experimental flow field was modelled on the basis of a 2D-streamtube formalism with an underlying homogeneous and isotropic transmissivity field. Tracer transport was modelled using the dual porosity medium approach, which is linked to the flow model by the flow porosity. Given the substantial pumping rates in the extraction borehole, the transport domain has a maximum width of a few centimetres only. It is concluded that both the uncertainty with regard to the length of individual fractures and the detailed geometry of the network along the flowpath between injection and extraction boreholes are not critical because flow is largely one-dimensional, whether through a single fracture or a network. Process identification and model calibration were based on a single uranine breakthrough (test PDT3), which clearly showed that matrix diffusion had to be included in the model even over the short experimental time scales, evidenced by a characteristic shape of the trailing edge of the breakthrough curve. Using the geological information and therefore considering limited matrix diffusion into a thin fault gouge horizon resulted in a good fit to the experiment. On the other hand, fresh granite was found not to interact noticeably with the tracers over the time scales of the experiments. While fracture-filling gouge materials are very efficient in retarding tracers over short periods of time (hours–days), their volume is very small and, with time progressing, retardation will be dominated by altered wall rock and, finally, by fresh granite. In such rocks, both porosity (and therefore the effective diffusion coefficient) and sorption Kds are more than one order of magnitude smaller compared to fault gouge, thus indicating that long-term retardation is expected to occur but to be less pronounced.

Solute transport in crystalline rocks at Äspö — II: Blind predictions, inverse modelling and lessons learnt from test STT1

Based on the results from detailed structural and petrological characterisation and on up-scaled laboratory values for sorption and diffusion, blind predictions were made for the STT1 dipole tracer test performed in the Swedish A¨ spo¨ Hard Rock Laboratory. The tracers used were nonsorbing, such as uranine and tritiated water, weakly sorbing 22Na+, 85Sr2 +, 47Ca2 +and more strongly sorbing 86Rb+, 133Ba2 +, 137Cs+. Our model consists of two parts: (1) a flow part based on a 2D-streamtube formalism accounting for the natural background flow field and with an underlying homogeneous and isotropic transmissivity field and (2) a transport part in terms of the dual porosity medium approach which is linked to the flow part by the flow porosity. The calibration of the model was done using the data from one single uranine breakthrough (PDT3). The study clearly showed that matrix diffusion into a highly porous material, fault gouge, had to be included in our model evidenced by the characteristic shape of the breakthrough curve and in line with geological observations. After the disclosure of the measurements, it turned out that, in spite of the simplicity of our model, the prediction for the nonsorbing and weakly sorbing tracers was fairly good. The blind prediction for the more strongly sorbing tracers was in general less accurate. The reason for the good predictions is deemed to be the result of the choice of a model structure strongly based on geological observation. The breakthrough curves were inversely modelled to determine in situ values for the transport parameters and to draw consequences on the model structure applied. For good fits, only one additional fracture family in contact with cataclasite had to be taken into account, but no new transport mechanisms had to be invoked. The in situ values for the effective diffusion coefficient for fault gouge are a factor of 2–15 larger than the laboratory data. For cataclasite, both data sets have values comparable to laboratory data. The extracted Kd values for the weakly sorbing tracers are larger than Swedish laboratory data by a factor of 25–60, but agree within a factor of 3–5 for the more strongly sorbing nuclides. The reason for the inconsistency concerning Kds is the use of fresh granite in the laboratory studies, whereas tracers in the field experiments interact only with fracture fault gouge and to a lesser extent with cataclasite both being mineralogically very different (e.g. clay-bearing) from the intact wall rock.

(Table 2) Trace element composition of seleceted clast and matrix materials from ODP Site 125-786

Shipboard examination of volcanic and sedimentary strata at Site 786 suggested that at least four types of breccias are present: flow-top breccias, associated with cooling and breakup on the upper surface of lava flows; autobreccias, formed by in-situ alteration at the base of flows; fault-gouge breccias; and true sedimentary breccias derived from weathering and erosion of underlying flows. It is virtually impossible to assess the origin of breccia matrix by textural and mineralogical analyses alone. However, it is fundamental for our understanding of breccia provenance to determine the source component of the matrix material. Whether the matrix is uniquely clastderived can be determined by geochemical fingerprinting. Trace elements that are immobile during weathering and alteration do not change their relative abundances. A contribution to the matrix from any source with an immobile trace element signature different from that of the clasts would appear as a perturbation of the trace element signature of the matrix. Trace element analysis of bulk samples from clasts and matrix material in individual breccia units was undertaken in a fashion similar to that used by Brimhall and Dietrich (1987, doi:10.1016/0016-7037(87)90070-6) in analyzing soil provenance: (1) to help distinguish between sedimentary and volcanic breccias, (2) to determine the degree of mixing and depth of erosion in sedimentary breccias, and (3) to analyze the local provenance of the individual breccia components (matrix and clasts). The following elements were analyzed by X-ray fluorescence (XRF): Rb, Sr, Ba, U, Zr, Cu, Zn, Ti, Cr, and V. Of these elements, Zr and Ti probably exhibit truly immobile behavior (Humphris and Thompson, 1978, doi:10.1016/0016-7037(78)90222-3 ). The remaining elements are useful as a reference for the extent of compositional change during the formation of matrix material (Brimhall and Dietrich, 1987, doi:10.1016/0016-7037(87)90070-6).

(Table T1) Magnetic anisotropy of samples from ODP Hole 180-1117A

At Site 1117, drilled during Leg 180 of the Ocean Drilling Program in the Woodlark Basin, we cored a fault zone and recovered fault gouge, mylonitized and brecciated gabbros, and undeformed gabbro. We measured the anisotropy of magnetic susceptibility for the rock samples. The susceptibilities of the fault gouge samples were lower than those of the undeformed gabbro, and those of deformed gabbros were lowest. The anisotropy degrees of the fault gouge samples were higher than those of the deformed and undeformed gabbros. Oblate magnetic fabrics were dominant in the samples from the fault zone.

Whole rock and elements concentrations from IODP Hole 343-C0019E

The microstructures, mineralogy and chemistry of four representative samples collected from cores extracted from the Japan Trench during Integrated Ocean Drilling Project Expedition 343, the Japan Trench Fast Drilling Project (JFAST) have been studied using optical microscopy, TEM, SEM, XRF, XRD and microprobe analyses. The samples provide a transect from relatively undeformed marine sediments in the hanging wall, to the undeformed footwall material, crossing the thrust interface between the Pacific and North American plate, where the fault slipped during the March 2011 Tohoku-Oki earthquake. Our preliminary results suggest that the low strength of JFAST fault gouge material is caused by the high amount of clay minerals (~ 60% smectite, ~ 14 illite). Clay minerals in the décollement (gouge) sample are partly replaced by newly formed manganese oxide, which precipitated from hydrothermal fluids. Dauphine twins were found in quartz grains of the décollement sample suggesting local high stress possible during seismic loading. Other microstructures cannot be assigned unambiguously to co-seismic or a-seismic faulting processes. The observed scaly clay fabric is consistent with observations in many other plate-boundary fault zones. Significant grain size reduction was found in the fault (decollement) zone sample. But a change in lithology of the fault material cannot be ruled out. Microstructures typical for a-seismic deformation like dissolution-precipitation features (e.g. dissolved grain boundaries, mineral alteration) occur in all JFAST core samples, but more frequently in the décollement sample.

The lattice solid model to simulate the physics of rocks and earthquakes: Incorporation of friction

The particle-based lattice solid model developed to study the physics of rocks and the nonlinear dynamics of earthquakes is refined by incorporating intrinsic friction between particles. The model provides a means for studying the causes of seismic wave attenuation, as well as frictional heat generation, fault zone evolution, and localisation phenomena. A modified velocity-Verlat scheme that allows friction to be precisely modelled is developed. This is a difficult computational problem given that a discontinuity must be accurately simulated by the numerical approach (i.e., the transition from static to dynamical frictional behaviour). This is achieved using a half time step integration scheme. At each half time step, a nonlinear system is solved to compute the static frictional forces and states of touching particle-pairs. Improved efficiency is achieved by adaptively adjusting the time step increment, depending on the particle velocities in the system. The total energy is calculated and verified to remain constant to a high precision during simulations. Numerical experiments show that the model can be applied to the study of earthquake dynamics, the stick-slip instability, heat generation, and fault zone evolution. Such experiments may lead to a conclusive resolution of the heat flow paradox and improved understanding of earthquake precursory phenomena and dynamics. (C) 1999 Academic Press.