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.

A unified friction description and its application to the simulation of frictional instability using the finite element method

This article first summarizes some available experimental results on the frictional behaviour of contact interfaces, and briefly recalls typical frictional experiments and relationships, which are applicable for rock mechanics, and then a unified description is obtained to describe the entire frictional behaviour. It is formulated based on the experimental results and applied with a stick and slip decomposition algorithm to describe the stick-slip instability phenomena, which can describe the effects observed in rock experiments without using the so-called state variable, thus avoiding related numerical difficulties. This has been implemented to our finite element code, which uses the node-to-point contact element strategy proposed by the authors to handle the frictional contact between multiple finite-deformation bodies with stick and finite frictional slip, and applied here to simulate the frictional behaviour of rocks to show its usefulness and efficiency.

ACcESS: Australia's contribution to the iSERVO Institute's development

The new Australian Computational Earth Systems Simulator research facility provides a virtual laboratory for studying the solid earth and its complex system behavior. The facility's capabilities complement those developed by overseas groups, thereby creating the infrastructure for an international computational solid earth research virtual observatory.

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.

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.

Capillary jet instability under the influence of gravity

The focus of this paper is on the effect of gravity stretching on disturbed capillary jet instability. Break-up and droplet formation under low flows are simulated using finite difference solution of a one-dimensional approximation of disturbed capillary jet instability chosen from the work by Eggers and Dupont (J. Fluid Mech. 155 (1994) 289). Experiments were conducted using water and aqueous glycerol solutions to compare with simulations. We use a gravity parameter, G, which quantifies gravity stretching by relating flow velocity, orifice size and acceleration and is the reciprocal of the Fronde number. The optimum disturbance frequency Omega(opt) was found to be inversely proportional to G. However, this relationship appears to be complex for the range of G's investigated. At low G, the relationship between Omega(opt) and G appears to be linear but takes on a weakly decaying like trend as G increases. As flows are lowered, the satellite-free regime decreases, although experimental observation found that merging of main and satellite drops sometimes offset this effect to result in monodispersed droplet trains post break-up. Viscosity did not significantly affect the relationship between the disturbance frequency and G, although satellite drops could be seen more clearly close to the upper limit for instability at high G's. It is possible to define regimes of satellite formation under low flows by considering local wavenumbers at the point of instability. (C) 2004 Elsevier Ltd. All rights reserved.

Complete modulational-instability gain spectrum of nonlinear quasi-phase-matching gratings

We consider plane waves propagating in quadratic nonlinear slab waveguides with nonlinear quasi-phase-matching gratings. We predict analytically and verify numerically the complete gain spectrum for transverse modulational instability, including hitherto undescribed higher-order gain bands. (C) 2004 Optical Society of America.

Theoretical investigation of convective instability in inclined and fluid-saturated three-dimensional fault zones.

The convective instability of pore-fluid flow in inclined and fluid-saturated three-dimensional fault zones has been theoretically investigated in this paper. Due to the consideration of the inclined three-dimensional fault zone with any values of the inclined angle, it is impossible to use the conventional linear stability analysis method for deriving the critical condition (i.e., the critical Rayleigh number) which can be used to investigate the convective instability of the pore-fluid flow in an inclined three-dimensional fault zone system. To overcome this mathematical difficulty, a combination of the variable separation method and the integration elimination method has been used to derive the characteristic equation, which depends on the Rayleigh number and the inclined angle of the inclined three-dimensional fault zone. Using this characteristic equation, the critical Rayleigh number of the system can be numerically found as a function of the inclined angle of the three-dimensional fault zone. For a vertically oriented three-dimensional fault zone system, the critical Rayleigh number of the system can be explicitly derived from the characteristic equation. Comparison of the resulting critical Rayleigh number of the system with that previously derived in a vertically oriented three-dimensional fault zone has demonstrated that the characteristic equation of the Rayleigh number is correct and useful for investigating the convective instability of pore-fluid flow in the inclined three-dimensional fault zone system. The related numerical results from this investigation have indicated that: (1) the convective pore-fluid flow may take place in the inclined three-dimensional fault zone; (2) if the height of the fault zone is used as the characteristic length of the system, a decrease in the inclined angle of the inclined fault zone stabilizes the three-dimensional fundamental convective flow in the inclined three-dimensional fault zone system; (3) if the thickness of the stratum is used as the characteristic length of the system, a decrease in the inclined angle of the inclined fault zone destabilizes the three-dimensional fundamental convective flow in the inclined three-dimensional fault zone system; and that (4) the shape of the inclined three-dimensional fault zone may affect the convective instability of pore-fluid flow in the system. (C) 2004 Published by Elsevier B.V.

Asset Price Instability and Policy Responses: The Legacy of Liberalization

The problem of asset price bubbles, and more generally of instability in the financial system, has been a matter of concern since the 1980s but has only recently moved to the center of the macroeconomic policy debate. The main concern with bubbles arises when they burst, imposing losses on investors holding the bubble assets and potentially on the financial institutions that have extended credit to them. Asset price volatility is an inevitable consequence of financial market liberalization and, in extreme cases, generates asset price bubbles, the bursting of which can impose substantial economic and social costs. Policy responses within the existing liberalized financial system face daunting levels of uncertainty and risk. Given the pattern of increasing asset market volatility over recent decades and the policy issues highlighted in this paper, the future looks uncertain. Another significant cycle of asset price movements, especially in one of the major economies, could see a fundamental revision of thinking about the costs and benefits of liberalized financial systems.

Low level microsatellite instability may be associated with reduced cancer specific survival in sporadic stage C colorectal carcinoma

Background: Colorectal cancers (CRCs) may be categorised according to the degree of microsatellite instability (MSI) exhibited, as MSI-high (MSI-H), MSI-low (MSI-L), or microsatellite stable (MSS). MSI-H status confers a survival advantage to patients with sporadic CRC. Aims: To determine if low levels of MSI are related to the clinicopathological features and prognosis of sporadic stage C CRC. Patients: A total of 255 patients who underwent resection for sporadic stage C CRC were studied. No patient received chemotherapy. Minimum follow up was five years. Methods: DNA extracted from archival malignant and non-malignant tissue was amplified by polymerase chain reaction using a panel of 11 microsatellites. MSI-H was defined as instability at greater than or equal to40% of markers, MSS as no instability, and MSI-L as instability at >0% but,40% of markers. Patients with MSI-H CRC were excluded from analysis as they have previously been shown to have better survival. Results: Thirty three MSI-L and 176 MSS CRCs were identified. There was no difference in biological characteristics or overall survival of MSI-L compared with MSS CRC but MSI-L was associated with poorer cancer specific survival (hazard ratio 2.0 (95% confidence interval 1.1-3.6)). Conclusions: Sporadic MSI-L and MSS CRCs have comparable clinicopathological features. Further studies are required to assess the impact of MSI-L on prognosis.

40Ar/39Ar geochronology of gold mineralization in Brasiliano strike-slip shear zones in the Barborema province, NE Brazil

40Ar/39Ar geochronology of muscovite and biotite grains genetically related to gold and Be–Ta–Li pegmatites from the Seridó Belt (Borborema province, NE Brazil) yield well-defined, reliable plateau ages. This information, combined with data about paragenetic and field relationships, reveals Cambro-Ordovician mineralization ages (520 and 500–506 Ma) for the orogenic gold deposits in the Seridó Belt. Biotite ages of 525±2 Ma, which represent the mean weighted results of the incremental heating analysis of six biotite single crystals, record the time of pegmatite emplacement and reactivation of Brasiliano/Pan-African strike-slip shear zones. These results, along with previous structural evolution studies, suggest that shear zones formed during the Brasiliano/Pan-African event were reactivated in the Upper Cambrian–Lower Ordovician. Mineralization occurs late in the history of the orogen.

A genetically informed study of the processes underlying the association between parental marital instability and offspring adjustment

Parental divorce is associated with problematic offspring adjustment, but the relation may be due to shared genetic or environmental factors. One way to test for these confounds is to study offspring of twins discordant for divorce. The current analyses used this design to separate the mechanisms responsible for the association between parental divorce, experienced either before or after the age of 16, and offspring well-being. The results were consistent with a causal role of divorce in earlier initiation of sexual intercourse and emotional difficulties, in addition to a greater probability of educational problems, depressed mood, and suicidal ideation. In contrast, the increased risk for cohabitation and earlier initiation of drug use was explained by selection factors, including genetic confounds.

On the chaotic instability of a nonsliding liquid-filled top with a small spheroidal base via Melnikov-Holmes-Marsden integrals

Chaotic orientations of a top containing a fluid filled cavity are investigated analytically and numerically under small perturbations. The top spins and rolls in nonsliding contact with a rough horizontal plane and the fluid in the ellipsoidal shaped cavity is considered to be ideal and describable by finite degrees of freedom. A Hamiltonian structure is established to facilitate the application of Melnikov-Holmes-Marsden (MHM) integrals. In particular, chaotic motion of the liquid-filled top is identified to be arisen from the transversal intersections between the stable and unstable manifolds of an approximated, disturbed flow of the liquid-filled top via the MHM integrals. The developed analytical criteria are crosschecked with numerical simulations via the 4th Runge-Kutta algorithms with adaptive time steps.