Dynamic rupture in a 3-D particle-based simulation of a rough planar fault

An appreciation of the physical mechanisms which cause observed seismicity complexity is fundamental to the understanding of the temporal behaviour of faults and single slip events. Numerical simulation of fault slip can provide insights into fault processes by allowing exploration of parameter spaces which influence microscopic and macroscopic physics of processes which may lead towards an answer to those questions. Particle-based models such as the Lattice Solid Model have been used previously for the simulation of stick-slip dynamics of faults, although mainly in two dimensions. Recent increases in the power of computers and the ability to use the power of parallel computer systems have made it possible to extend particle-based fault simulations to three dimensions. In this paper a particle-based numerical model of a rough planar fault embedded between two elastic blocks in three dimensions is presented. A very simple friction law without any rate dependency and no spatial heterogeneity in the intrinsic coefficient of friction is used in the model. To simulate earthquake dynamics the model is sheared in a direction parallel to the fault plane with a constant velocity at the driving edges. Spontaneous slip occurs on the fault when the shear stress is large enough to overcome the frictional forces on the fault. Slip events with a wide range of event sizes are observed. Investigation of the temporal evolution and spatial distribution of slip during each event shows a high degree of variability between the events. In some of the larger events highly complex slip patterns are observed.

Efficient accommodation of local minima in watershed model calibration

The Gauss-Marquardt-Levenberg (GML) method of computer-based parameter estimation, in common with other gradient-based approaches, suffers from the drawback that it may become trapped in local objective function minima, and thus report optimized parameter values that are not, in fact, optimized at all. This can seriously degrade its utility in the calibration of watershed models where local optima abound. Nevertheless, the method also has advantages, chief among these being its model-run efficiency, and its ability to report useful information on parameter sensitivities and covariances as a by-product of its use. It is also easily adapted to maintain this efficiency in the face of potential numerical problems (that adversely affect all parameter estimation methodologies) caused by parameter insensitivity and/or parameter correlation. The present paper presents two algorithmic enhancements to the GML method that retain its strengths, but which overcome its weaknesses in the face of local optima. Using the first of these methods an intelligent search for better parameter sets is conducted in parameter subspaces of decreasing dimensionality when progress of the parameter estimation process is slowed either by numerical instability incurred through problem ill-posedness, or when a local objective function minimum is encountered. The second methodology minimizes the chance of successive GML parameter estimation runs finding the same objective function minimum by starting successive runs at points that are maximally removed from previous parameter trajectories. As well as enhancing the ability of a GML-based method to find the global objective function minimum, the latter technique can also be used to find the locations of many non-global optima (should they exist) in parameter space. This can provide a useful means of inquiring into the well-posedness of a parameter estimation problem, and for detecting the presence of bimodal parameter and predictive probability distributions. The new methodologies are demonstrated by calibrating a Hydrological Simulation Program-FORTRAN (HSPF) model against a time series of daily flows. Comparison with the SCE-UA method in this calibration context demonstrates a high level of comparative model run efficiency for the new method. (c) 2006 Elsevier B.V. All rights reserved.

Fault evaluation for security-critical communications devices

Communications devices for government or military applications must keep data secure, even when their electronic components fail. Combining information flow and risk analyses could make fault-mode evaluations for such devices more efficient and cost-effective.

Holocene depositional history of the Burdekin River delta of northeastern Australia: A model for a low-accommodation, highstand delta

The Burdekin River of northeastern Australia has constructed a substantial delta during the Holocene (delta plain area 1260 km2). The vertical succession through this delta comprises (1) a basal, coarse-grained transgressive lag overlying a continental omission surface, overlain by (2) a mud interval deposited as the coastal region was inundated by the postglacially rising sea, in turn overlain by (3) a generally sharp-based sand unit deposited principally in channel and mouth-bar environments with lesser volumes of floodplain and coastal facies. The Holocene Burdekin Delta was constructed as a series of at least thirteen discrete delta lobes, formed as the river avulsed. Each lobe consists of a composite sand body typically 5-8 m thick. The oldest lobes, formed during the latter stages of the postglacial sea-level rise (10-5.5 kyr BP), are larger than those formed during the highstand (5.5-3 kyr BP), which are in turn larger than those formed during the most recent slight sea-level lowering and stillstand (3-0 kyr BP). Radiocarbon ages and other stratigraphic data indicate that inter-avulsion period has decreased through time coincident with the decrease in delta lobe area. The primary control on Holocene delta architecture appears to have been a change from a pluvial climate known to characterize the region 12-4 kyr BP to the present drier, ENSO-dominated climate. In addition to decreasing the sediment supply via lower rates of chemical weathering, this change may have contributed to the shorter avulsion period by facilitating extreme variability of discharge. More frequent avulsion may also have been facilitated by the lengthening of the delta-plain channels as the system prograded seaward. Copyright © 2006, SEPM (Society for Sedimentary Geology).

Scattering and tangential momentum accommodation at a 2D adsorbate-solid interface

We calculate tangential momentum coefficients for the exchange of momentum between molecules in transport and the internal surface of a membrane pore, modelled as a simple atomic structure. We introduce a local specular reflection (LSR) hypothesis, which states that impinging molecules undergo mirror-like reflection in a plane tangent to a surface atom at the point of impact. As a consequence, the components of the velocity, parallel to the direction of flow will (in general) change on impact. The overall effect is a loss of tangential momentum, since more is lost in the upstream direction than is gained in the downstream direction. The loss of tangential momentum is greater when the size ratio of fluid to solid atom is small, allowing more steeply inclined impact planes to become accessible to the fluid phase molecules. (c) 2005 Elsevier B.V. All rights reserved.

Clay mineralogical and isotopic (K-Ar, delta O-18, delta D) constraints on the evolution of the North Anatolian Fault Zone, Turkey

This study presents the first attempt to constrain the evolution of the North Anatolian Fault Zone (NAFZ) by age dating and isotope tracing of clay minerals formed during near-surface faulting. Extensive illitic clay mineralisation occurred along the NAFZ related to hydrothermal alteration of the fault gouges and pseudotachylytes. Samples representing the pre-fault protoliths outside the fault zone do not contain authigenic illitic clay minerals indicating that hydrothermal processes were confined to the areas within the fault zone. K-Ar age data indicate that the hydrothermal system and the associated illite authigenesis initiated at similar to 57 Ma. This process is interpreted to reflect the onset of significant strike-slip or transtensional faulting immediately after the continental collision related to the closure of the Neotethys Ocean. Following the initiation of the fault movements in the latest Paleocene-Early Eocene, displacements along the NAFZ have continued, with probably intensified fault activities at similar to 26 Ma and later than similar to 8 Ma. Oxygen isotope compositions of the illitic clays from different locations along the NAFZ are similar, with narrow ranges in delta O-18 values indicating clay precipitation from fluids with similar oxygen isotope compositions and crystallisation temperatures. The delta O-18 and delta D values of the calculated fluid isotopic composition (delta O-18=5.9 parts per thousand to 11.2 parts per thousand, delta D=-59 parts per thousand to -73 parts per thousand) are consistent with metamorphic and magmatic origin of fluids mobilised during active tectonism. The interpretation of the fluid flow history of the NAFZ is in agreement with that reported previously for some well-known large-scale high-angle fault zones, which similarly developed along collisional-type orogenic belts and are commonly associated with significant mesothermal ore mineralisation. (c) 2005 Elsevier B.V. All rights reserved.

Construction of an Intraplate Fault System Model of South Australia, and Simulation Tool for the iSERVO Institute Seed Project

To foster ongoing international cooperation beyond ACES (APEC Cooperation for Earthquake Simulation) on the simulation of solid earth phenomena, agreement was reached to work towards establishment of a frontier international research institute for simulating the solid earth: iSERVO = International Solid Earth Research Virtual Observatory institute (http://www.iservo.edu.au). This paper outlines a key Australian contribution towards the iSERVO institute seed project, this is the construction of: (1) a typical intraplate fault system model using practical fault system data of South Australia (i.e., SA interacting fault model), which includes data management and editing, geometrical modeling and mesh generation; and (2) a finite-element based software tool, which is built on our long-term and ongoing effort to develop the R-minimum strategy based finite-element computational algorithm and software tool for modelling three-dimensional nonlinear frictional contact behavior between multiple deformable bodies with the arbitrarily-shaped contact element strategy. A numerical simulation of the SA fault system is carried out using this software tool to demonstrate its capability and our efforts towards seeding the iSERVO Institute.

Effect of rolling on dissipation in fault gouges

Sliding and rolling are two outstanding deformation modes in granular media. The first one induces frictional dissipation whereas the latter one involves deformation with negligible resistance. Using numerical simulations on two-dimensional shear cells, we investigate the effect of the grain rotation on the energy dissipation and the strength of granular materials under quasistatic shear deformation. Rolling and sliding are quantified in terms of the so-called Cosserat rotations. The observed spontaneous formation of vorticity cells and clusters of rotating bearings may provide an explanation for the long standing heat flow paradox of earthquake dynamics.

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.

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.