Long-range automaton models of earthquakes: Power-law accelerations, correlation evolution, and mode-switching

We introduce a conceptual model for the in-plane physics of an earthquake fault. The model employs cellular automaton techniques to simulate tectonic loading, earthquake rupture, and strain redistribution. The impact of a hypothetical crustal elastodynamic Green's function is approximated by a long-range strain redistribution law with a r(-p) dependance. We investigate the influence of the effective elastodynamic interaction range upon the dynamical behaviour of the model by conducting experiments with different values of the exponent (p). The results indicate that this model has two distinct, stable modes of behaviour. The first mode produces a characteristic earthquake distribution with moderate to large events preceeded by an interval of time in which the rate of energy release accelerates. A correlation function analysis reveals that accelerating sequences are associated with a systematic, global evolution of strain energy correlations within the system. The second stable mode produces Gutenberg-Richter statistics, with near-linear energy release and no significant global correlation evolution. A model with effectively short-range interactions preferentially displays Gutenberg-Richter behaviour. However, models with long-range interactions appear to switch between the characteristic and GR modes. As the range of elastodynamic interactions is increased, characteristic behaviour begins to dominate GR behaviour. These models demonstrate that evolution of strain energy correlations may occur within systems with a fixed elastodynamic interaction range. Supposing that similar mode-switching dynamical behaviour occurs within earthquake faults then intermediate-term forecasting of large earthquakes may be feasible for some earthquakes but not for others, in alignment with certain empirical seismological observations. Further numerical investigation of dynamical models of this type may lead to advances in earthquake forecasting research and theoretical seismology.

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.

Effect of soil composition and dissolved organic matter on pesticide sorption

The effect of the solid and dissolved organic matter fractions, mineral composition and ionic strength of the soil solution on the sorption behaviour of pesticides were studied. A number of soils, chosen so as to have different clay mineral and organic carbon content, were used to study the sorption of the pesticides atrazine (6-chloro-N-2-ethyl-N-4-isopropyl-1,3,5-triazine-2,4-diamine), 2,4-D ((2,4-dichlorophenoxy) acetic acid), isoproturon (3-(4-isopropylphenyl)1,1-dimethylurea) and paraquat (1,1'-dimethyl-4,4'-bipyridinium) in the presence of low and high levels of dissolved organic carbon and different background electrolytes. The sorption behaviour of atrazine, isoproturon and paraquat was dominated by the solid state soil components and the presence of dissolved organic matter had little effect. The sorption of 2,4-D was slightly affected by the soluble organic matter in the soil. However, this effect may be due to competition for adsorption sites between the pesticide and the soluble organic matter rather than due to a positive interaction between the pesticide and the soluble fraction of soil organic matter. It is concluded that the major factor governing the sorption of these pesticides is the solid state organic fraction with the clay mineral content also making a significant contribution. The dissolved organic carbon fraction of the total organic carbon in the soil and the ionic strength of the soil solution appear to have little or no effect on the sorption/transport characteristics of these pesticides over the range of concentrations studied. (C) 2002 Elsevier Science B.V. All rights reserved.