Simulation of the micro-physics of rocks using LSMearth

The particle-based Lattice Solid Model (LSM) was developed to provide a basis to study the physics of rocks and the nonlinear dynamics of earthquakes (MORA and PLACE, 1994; PLACE and MORA, 1999). A new modular and flexible LSM approach has been developed that allows different microphysics to be easily included in or removed from the model. The approach provides a virtual laboratory where numerical experiments can easily be set up and all measurable quantities visualised. The proposed approach provides a means to simulate complex phenomena such as fracturing or localisation processes, and enables the effect of different micro-physics on macroscopic behaviour to be studied. The initial 2-D model is extended to allow three-dimensional simulations to be performed and particles of different sizes to be specified. Numerical bi-axial compression experiments under different confining pressure are used to calibrate the model. By tuning the different microscopic parameters (such as coefficient of friction, microscopic strength and distribution of grain sizes), the macroscopic strength of the material and can be adjusted to be in agreement with laboratory experiments, and the orientation of fractures is consistent with the theoretical value predicted based on Mohr-Coulomb diagram. Simulations indicate that 3-D numerical models have different macroscopic properties than in 2-D and, hence, the model must be recalibrated for 3-D simulations. These numerical experiments illustrate that the new approach is capable of simulating typical rock fracture behaviour. The new model provides a basis to investigate nucleation, rupture and slip pulse propagation in complex fault zones without the previous model limitations of a regular low-level surface geometry and being restricted to two-dimensions.

Simulation of the influence of rate- and state-dependent friction on the macroscopic behavior of complex fault zones with the lattice solid model

In order to understand the earthquake nucleation process, we need to understand the effective frictional behavior of faults with complex geometry and fault gouge zones. One important aspect of this is the interaction between the friction law governing the behavior of the fault on the microscopic level and the resulting macroscopic behavior of the fault zone. Numerical simulations offer a possibility to investigate the behavior of faults on many different scales and thus provide a means to gain insight into fault zone dynamics on scales which are not accessible to laboratory experiments. Numerical experiments have been performed to investigate the influence of the geometric configuration of faults with a rate- and state-dependent friction at the particle contacts on the effective frictional behavior of these faults. The numerical experiments are designed to be similar to laboratory experiments by DIETERICH and KILGORE (1994) in which a slide-hold-slide cycle was performed between two blocks of material and the resulting peak friction was plotted vs. holding time. Simulations with a flat fault without a fault gouge have been performed to verify the implementation. These have shown close agreement with comparable laboratory experiments. The simulations performed with a fault containing fault gouge have demonstrated a strong dependence of the critical slip distance D-c on the roughness of the fault surfaces and are in qualitative agreement with laboratory experiments.

Parallel simulation system for earthquake generation: Fault analysis modules and parallel coupling analysis

Solid earth simulations have recently been developed to address issues such as natural disasters, global environmental destruction and the conservation of natural resources. The simulation of solid earth phenomena involves the analysis of complex structures including strata, faults, and heterogeneous material properties. Simulation of the generation and cycle of earthquakes is particularly important, but such simulations require the analysis of complex fault dynamics. GeoFEM is a parallel finite-element analysis system intended for solid earth field phenomena problems. This paper describes recent development in the GeoFEM project for the simulation of earthquake generation and cycles.

Tool support for testing concurrent Java components

Concurrent programs are hard to test due to the inherent nondeterminism. This paper presents a method and tool support for testing concurrent Java components. Too[ support is offered through ConAn (Concurrency Analyser), a too] for generating drivers for unit testing Java classes that are used in a multithreaded context. To obtain adequate controllability over the interactions between Java threads, the generated driver contains threads that are synchronized by a clock. The driver automatically executes the calls in the test sequence in the prescribed order and compares the outputs against the expected outputs specified in the test sequence. The method and tool are illustrated in detail on an asymmetric producer-consumer monitor. Their application to testing over 20 concurrent components, a number of which are sourced from industry and were found to contain faults, is presented and discussed.

Functionally partitioned module-based programmable architecture for wireless base-band processing

A specialised reconfigurable architecture is targeted at wireless base-band processing. It is built to cater for multiple wireless standards. It has lower power consumption than the processor-based solution. It can be scaled to run in parallel for processing multiple channels. Test resources are embedded on the architecture and testing strategies are included. This architecture is functionally partitioned according to the common operations found in wireless standards, such as CRC error correction, convolution and interleaving. These modules are linked via Virtual Wire Hardware modules and route-through switch matrices. Data can be processed in any order through this interconnect structure. Virtual Wire ensures the same flexibility as normal interconnects, but the area occupied and the number of switches needed is reduced. The testing algorithm scans all possible paths within the interconnection network exhaustively and searches for faults in the processing modules. The testing algorithm starts by scanning the externally addressable memory space and testing the master controller. The controller then tests every switch in the route-through switch matrix by making loops from the shared memory to each of the switches. The local switch matrix is also tested in the same way. Next the local memory is scanned. Finally, pre-defined test vectors are loaded into local memory to check the processing modules. This paper compares various base-band processing solutions. It describes the proposed platform and its implementation. It outlines the test resources and algorithm. It concludes with the mapping of Bluetooth and GSM base-band onto the platform.

40Ar/39Ar constraints on a temporal link between gold mineralization, magmatism, and continental margin transtension in the Jiaodong Gold Province, eastern China

The Jiaodong gold province is the largest gold repository in China. Both mineralization and granitoid hosts are spatially related to the crustal-scale Tan-Lu strike-slip fault system, which developed along the Mesozoic continental margin in eastern China. A series of Ar-40/Ar-39 laser incremental heating analyses of hydrothermal sericite/muscovite from three major gold deposits (Jiaojia, Xincheng, and Wangershan) and igneous biotite from the granodiorite hosts were performed to establish a possible temporal link between gold mineralization, magmatism, and movement along the Tan-Lu fault zone. Magmatic biotite crystals yield well-defined and concordant plateau ages between 124.5+/-0.4 Ma and 124.0+/-0.4 Ma (2sigma), whereas sericite and muscovite samples (a total of 30 single separates) give reproducible plateau ages ranging from 121.0+/-0.4 Ma to 119.2+/-0.2 Ma (2sigma). An integration of our Ar-40/Ar-39 results with age data from other major gold deposits in Jiaodong demonstrates that widespread gold mineralization occurred contemporaneously during a 2-3-m.yr. period. Most gold deposits show intimate spatial associations with abundant mafic to intermediate dikes. The mafic dikes have K-Ar ages of 123.5-119.6 Ma, in excellent agreement with those of the gold deposits. These newly obtained Ar-40/Ar-39 ages, in combination with other independent geological and geochronological data on granodioritic intrusions (130-126 Ma), volcanic rocks (1243.6-114.7 Ma), and deformed rocks within strike-slip faults (132-120 Ma) in Jiaodong or adjacent areas, also support the idea that gold mineralization postdated the granodioritic magmatism but was contemporaneous with mafic magmatism and volcanism, all controlled by the transtensional motion along the Tan-Lu fault in the Early Cretaceous.

Early-Middle Jurassic stratigraphy of the Fortrose-Chaslands region, southernmost South Island, New Zealand

The 40 km of coastline from Fortrose to Chaslands Mistake (southeastern South Island, New Zealand) comprises sediments that are part of the Early-Middle Jurassic of the Murihiku Terrane. The sediments are dominantly fluvial with some marine beds and alluvial fan deposition, and display an evolution of fluvial style which progresses from perennial flow to seasonal flow. The McPhee Cove Conglomerate is a prominent unit to the north. It has been used to separate two formations which would otherwise, on inherent lithological grounds, be difficult to distinguish. This paper discusses several similar conglomerates which occur in the south, but which are separated from the type area of the McPhee Conglomerate by major tectonic disruption. Hence, the existing lithostratigraphic nomenclature to the north, including the McPhee Cove Conglomerate, cannot be simply extended southwards. The Fortrose-Chaslands area appears to consist of two tectonic blocks, the Slope Point Block and the Brothers Block, which are separated from each other and from the adjacent Papatowai Block by major strike faults (or fault zones). A change is proposed to the existing stratigraphy which involves recognising all terrestrial sediments as part of the False Island Formation. Four prominent clast-supported conglomerate horizons are named as members of the False Islet Formation: the White Head Conglomerate, Black Bluff Conglomerate. Hoiho Conglomerate, and Slope Point Conglomerate Members. The latter contains five named conglomerate beds.

Rare earth element geochemistry of Late Devonian reefal carbonates, canning basin, Western Australia: Confirmation of a seawater REE proxy in ancient limestones

Rare earth element and yttrium (REE+Y) concentrations were determined in 49 Late Devonian reefal carbonates from the Lennard Shelf, Canning Basin, Western Australia. Shale-normalized (SN) REE+Y patterns of the Late Devonian samples display features consistent with the geochemistry of well-oxygenated, shallow seawater. A variety of different ancient limestone components, including microbialites, some skeletal carbonates (stromatoporoids), and cements, record seawater-like REE+Y signatures. Contamination associated with phosphate, Fe-oxides and shale was tested quantitatively, and can be discounted as the source of the REE+Y patterns. Co-occurring carbonate components that presumably precipitated from the same seawater have different relative REE concentrations, but consistent REE+Y patterns. Clean Devonian early marine cements (n = 3) display REE+Y signatures most like that of modern open ocean seawater and the highest Y/Ho ratios (e.g., 59) and greatest light REE (LREE) depletion (average Nd-SN/Yb-SN = 0.413, SD = 0.076). However, synsedimentary cements have the lowest REE concentrations (e.g., 405 ppb). Non-contaminated Devonian microbialite samples containing a mixture of the calcimicrobe Renalcis and micritic thrombolite aggregates in early marine cement (n = 11) have the highest relative REE concentrations of tested carbonates (average total REE = 11.3 ppm). Stromatoporoid skeletons, unlike modern corals, algae and molluscs, also contain well-developed, seawater-like REE patterns. Samples from an estuarine fringing reef have very different REE+Y patterns with LREE enrichment (Nd-SN/Yb-SN > 1), possibly reflecting inclusion of estuarine colloidal material that contained preferentially scavenged LREE from a nearby riverine input source. Hence, Devonian limestones provide a proxy for marine REE geochemistry and allow the differentiation of co-occurring water masses on the ancient Lennard Shelf. Although appropriate partition coefficients for quantification of Devonian seawater REE concentrations from out data are unknown, hypothetical Devonian Canning Basin seawater REE patterns were obtained with coefficients derived from modern natural proxies and experimental values. Resulting Devonian seawater patterns are slightly enriched in LREE compared to most modem seawaters and suggest higher overall REE concentrations, but are very similar to seawaters from regions with high terrigenous inputs. Our results suggest that most limestones should record important aspects of the REE geochemistry of the waters in which they precipitated, provided they are relatively free of terrigenous contamination and major diagenetic alteration from fluids with high, non-seawater-like REE contents. Hence, we expect that many other ancient limestones will serve as seawater REE proxies, and thereby provide information on paleoceanography, paleogeography and geochemical evolution of the oceans. Copyright (C) 2004 Elsevier Ltd.

Finite element analysis of fault bend influence on stick-slip instability along an intra-plate fault

Earthquakes have been recognized as resulting from stick-slip frictional instabilities along the faults between deformable rocks. A three-dimensional finite-element code for modeling the nonlinear frictional contact behaviors between deformable bodies with the node-to-point contact element strategy has been developed and applied here to investigate the fault geometry influence on the nucleation and development process of the stick-slip instability along an intra-plate fault through a typical fault bend model, which has a pre-cut fault that is artificially bent by an angle of 5.6degrees at the fault center. The numerical results demonstrate that the geometry of the fault significantly affects nucleation, termination and restart of the stick-slip instability along the intra-plate fault, and all these instability phenomena can be well simulated using the current finite-element algorithm.

Earthquake statistics in a Block Slider Model and a fully dynamic Fault Model

We examine the event statistics obtained from two differing simplified models for earthquake faults. The first model is a reproduction of the Block-Slider model of Carlson et al. (1991), a model often employed in seismicity studies. The second model is an elastodynamic fault model based upon the Lattice Solid Model (LSM) of Mora and Place (1994). We performed simulations in which the fault length was varied in each model and generated synthetic catalogs of event sizes and times. From these catalogs, we constructed interval event size distributions and inter-event time distributions. The larger, localised events in the Block-Slider model displayed the same scaling behaviour as events in the LSM however the distribution of inter-event times was markedly different. The analysis of both event size and inter-event time statistics is an effective method for comparative studies of differing simplified models for earthquake faults.

A parallel implementation of the lattice solid model for the simulation of rock mechanics and earthquake dynamics

The Lattice Solid Model has been used successfully as a virtual laboratory to simulate fracturing of rocks, the dynamics of faults, earthquakes and gouge processes. However, results from those simulations show that in order to make the next step towards more realistic experiments it will be necessary to use models containing a significantly larger number of particles than current models. Thus, those simulations will require a greatly increased amount of computational resources. Whereas the computing power provided by single processors can be expected to increase according to Moore's law, i.e., to double every 18-24 months, parallel computers can provide significantly larger computing power today. In order to make this computing power available for the simulation of the microphysics of earthquakes, a parallel version of the Lattice Solid Model has been implemented. Benchmarks using large models with several millions of particles have shown that the parallel implementation of the Lattice Solid Model can achieve a high parallel-efficiency of about 80% for large numbers of processors on different computer architectures.

Information flow analysis for fail-secure devices

Information security devices must preserve security properties even in the presence of faults. This in turn requires a rigorous evaluation of the system behaviours resulting from component failures, especially how such failures affect information flow. We introduce a compositional method of static analysis for fail-secure behaviour. Our method uses reachability matrices to identify potentially undesirable information flows based on the fault modes of the system's components.

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.

In situ observation of bainite growth during isothermal holding

The mechanism of bainite growth has been investigated using in situ transmission electron microscopy observations. It was found that, in a number of alloys studied, a bainitic embryo is made of basic transformation units. These units are either a group of stacking faults or, in two dimensions, a series of parallelograms of different sizes. Thickening/widening of the bainite embryo takes place through shear along the stacking fault planes or twining planes. The bainite embryo is elongated by the formation of new transformation units at both tips of the bainite plate. The three-dimensional morphology of bainite is a convex tens-like lath. It is believed that the bainite embryo grows by shearing, which is controlled by the diffusion of solute atoms during the transformation. As the growth rate is much lower than that of martensite, it is therefore detectable. (c) 2006 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Multiagent realization of prediction-based diagnosis and loss prevention

A multiagent diagnostic system implemented in a Protege-JADE-JESS environment interfaced with a dynamic simulator and database services is described in this paper. The proposed system architecture enables the use of a combination of diagnostic methods from heterogeneous knowledge sources. The process ontology and the process agents are designed based on the structure of the process system, while the diagnostic agents implement the applied diagnostic methods. A specific completeness coordinator agent is implemented to coordinate the diagnostic agents based on different methods. The system is demonstrated on a case study for diagnosis of faults in a granulation process based on HAZOP and FMEA analysis.