Evolution of stress deficit and changing rates of seismicity in cellular automaton models of earthquake faults

We investigate the internal dynamics of two cellular automaton models with heterogeneous strength fields and differing nearest neighbour laws. One model is a crack-like automaton, transferring ail stress from a rupture zone to the surroundings. The other automaton is a partial stress drop automaton, transferring only a fraction of the stress within a rupture zone to the surroundings. To study evolution of stress, the mean spectral density. f(k(r)) of a stress deficit held is: examined prior to, and immediately following ruptures in both models. Both models display a power-law relationship between f(k(r)) and spatial wavenumber (k(r)) of the form f(k(r)) similar tok(r)(-beta). In the crack model, the evolution of stress deficit is consistent with cyclic approach to, and retreat from a critical state in which large events occur. The approach to criticality is driven by tectonic loading. Short-range stress transfer in the model does not affect the approach to criticality of broad regions in the model. The evolution of stress deficit in the partial stress drop model is consistent with small fluctuations about a mean state of high stress, behaviour indicative of a self-organised critical system. Despite statistics similar to natural earthquakes these simplified models lack a physical basis. physically motivated models of earthquakes also display dynamical complexity similar to that of a critical point system. Studies of dynamical complexity in physical models of earthquakes may lead to advancement towards a physical theory for earthquakes.

Psychosocial sequelae of the 1989 Newcastle earthquake .2. Exposure and morbidity profiles during the first 2 years post-disaster

Background. A sample of 1089 Australian adults was selected for the longitudinal component of the Quake Impact Study, a 2-year, four-phase investigation of the psychosocial effects of the 1989 Newcastle earthquake. Of these, 845 (78%) completed a survey 6 months post-disaster as well as one or more of the three follow-up surveys. Methods. The phase I survey was used to construct dimensional indices of self-reported exposure to threat the disruption and also to classify subjects by their membership of five 'at risk' groups (the injured; the displaced; owners of damaged small businesses; helpers in threat and non-threat situations). Psychological morbidity was assessed at each phase using the 12-item General Health Questionnaire (GHQ-12) and the Impact of Event Scale (IES). Results. Psychological morbidity declined over time but tended to stabilize at about 12 months post-disaster for general morbidity (GHQ-12) and at about 18 months for trauma-related (IES) morbidity. Initial exposure to threat and/or disruption were significant predictors of psychological morbidity throughout the study and had superior predictive power to membership of the targeted 'at risk' groups. The degree of ongoing disruption and other life events since the earthquake were also significant predictors of morbidity. The injured reported the highest levels of distress, but there was a relative absence of morbidity among the helpers. Conclusions. Future disaster research should carefully assess the threat and disruption experiences of the survivors at the time of the event and monitor ongoing disruptions in the aftermath in order to target interventions more effectively.

Psychosocial sequelae of the 1989 Newcastle earthquake .3. Role of vulnerability factors in postdisaster morbidity

Background. This paper examines the contributions of dispositional and non-dispositional factors to post-disaster psychological morbidity. Data reported are from the 845 participants in the longitudinal component of the Quake Impact Study. Methods. The phase 1 survey was used to construct dimensional indices of threat and disruption exposure. Subsequently, a range of dispositional characteristics were measured, including neuroticism, personal hopefulness and defence style. The main morbidity measures were the General Health Questionnaire (GHQ-12) and Impact of Event Scale (IES). Results. Dispositional characteristics were the best predictors of psychological morbidity throughout the 2 years post-disaster, contributing substantially more to the variance in morbidity (12-39%) than did initial exposure (5-12%), but the extent of their contribution was greater for general (GHQ-12) than for post-traumatic (IES) morbidity. Among the non-dispositional factors, avoidance coping contributed equally to general and post-traumatic morbidity (pr = 0.24). Life events since the earthquake (pr = 0.18), poor social relationships (pr = -0.25) and ongoing earthquake-related disruptions (pr = 0.22) also contributed to general morbidity, while only the latter contributed significantly to post-traumatic morbidity (pr = 0.15). Conclusions. Medium-term post-earthquake morbidity appears to be a function of multiple factors whose contributions vary depending on the type of morbidity experienced and include trait vulnerability, the nature and degree of initial exposure, avoidance coping and the nature and severity of subsequent events.

A synthesis of the findings from the Quake Impact Study: A two year investigation of the psychosocial sequelae of the 1989 Newcastle earthquake

This paper summarises the major findings from the Quake Impact Study (QIS), a four-phase longitudinal project that was conducted in the aftermath of the 1989 Newcastle (Australia) earthquake. A total of 3,484 subjects participated in at least one component of the QIS, comprising a stratified sample of 3,007 drawn from community electoral rolls and 477 from specially targeted supplementary samples (the injured, the displaced, the owners of damaged businesses, and the helpers). Subjects' initial earthquake experiences were rated in terms of weighted indices of exposure to threat and disruption. Psychological morbidity was measured at each phase using the General Health Questionnaire (GHQ-12) and the Impact of Event Scale (IES). Selected findings and key conclusions are presented for each of six areas of investigation: service utilisation during the first 6 months post-disaster; patterns of earthquake experience and short-term (6-month) psychosocial outcome; earthquake exposure and medium term (2-year) psychosocial outcome; vulnerability factors and medium-term psychosocial outcome: specific community groups at increased risk (e.g., the elderly and immigrants from non-English-speaking backgrounds); the effects of stress debriefing for helpers. Threshold morbidity (i.e., likely caseness) rates are also presented for a broad range of subgroups. In addition to presenting an overview of the QIS, this paper synthesises the major findings and discusses their implications for future disaster management and research from a mental health perspective.

Stress correlation function evolution in lattice solid elasto-dynamic models of shear and fracture zones and earthquake prediction

It has been argued that power-law time-to-failure fits for cumulative Benioff strain and an evolution in size-frequency statistics in the lead-up to large earthquakes are evidence that the crust behaves as a Critical Point (CP) system. If so, intermediate-term earthquake prediction is possible. However, this hypothesis has not been proven. If the crust does behave as a CP system, stress correlation lengths should grow in the lead-up to large events through the action of small to moderate ruptures and drop sharply once a large event occurs. However this evolution in stress correlation lengths cannot be observed directly. Here we show, using the lattice solid model to describe discontinuous elasto-dynamic systems subjected to shear and compression, that it is for possible correlation lengths to exhibit CP-type evolution. In the case of a granular system subjected to shear, this evolution occurs in the lead-up to the largest event and is accompanied by an increasing rate of moderate-sized events and power-law acceleration of Benioff strain release. In the case of an intact sample system subjected to compression, the evolution occurs only after a mature fracture system has developed. The results support the existence of a physical mechanism for intermediate-term earthquake forecasting and suggest this mechanism is fault-system dependent. This offers an explanation of why accelerating Benioff strain release is not observed prior to all large earthquakes. The results prove the existence of an underlying evolution in discontinuous elasto-dynamic, systems which is capable of providing a basis for forecasting catastrophic failure and earthquakes.

Load-unload response ratio and accelerating moment/energy release critical region scaling and earthquake prediction

The main idea of the Load-Unload Response Ratio (LURR) is that when a system is stable, its response to loading corresponds to its response to unloading, whereas when the system is approaching an unstable state, the response to loading and unloading becomes quite different. High LURR values and observations of Accelerating Moment/Energy Release (AMR/AER) prior to large earthquakes have led different research groups to suggest intermediate-term earthquake prediction is possible and imply that the LURR and AMR/AER observations may have a similar physical origin. To study this possibility, we conducted a retrospective examination of several Australian and Chinese earthquakes with magnitudes ranging from 5.0 to 7.9, including Australia's deadly Newcastle earthquake and the devastating Tangshan earthquake. Both LURR values and best-fit power-law time-to-failure functions were computed using data within a range of distances from the epicenter. Like the best-fit power-law fits in AMR/AER, the LURR value was optimal using data within a certain epicentral distance implying a critical region for LURR. Furthermore, LURR critical region size scales with mainshock magnitude and is similar to the AMR/AER critical region size. These results suggest a common physical origin for both the AMR/AER and LURR observations. Further research may provide clues that yield an understanding of this mechanism and help lead to a solid foundation for intermediate-term earthquake prediction.

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.

Assessing the performance of a nucleus estate and smallholder scheme for oil palm production in West Sumatra: a stochastic frontier analysis

In this paper, we study the performance of smallholders in a nucleus estate and smallholder (NES) scheme in oil palm production schemein West Sumatra by measuring their technical efficiency using a stochastic frontier production function. Our results indicate a mean technical efficiency of 66%, which is below what we would have expected given the uniformity of the climate, soils and plantation construction among the sample farmers. The use of progressive farmers as a means of disseminating extension advice does not appear to have been successful, and more rigorous farmer selection procedures need to be put in place for similar schemes and for general agricultural extension in future. No clear relationship was established between technical efficiency and the use of female labour, suggesting there is no need to target extension services specifically at female labourers in the household. Finally, education was found to have an unexpectedly negative impact on technical efficiency, indicating that farmers with primary education may be more important than those with secondary and tertiary education as targets of development schemes and extension programs entailing non-formal education. (C) 2003 Elsevier Ltd. All rights reserved.

Statistical tests of load-unload response ratio signals by lattice solid model: Implication to tidal triggering and earthquake prediction

Statistical tests of Load-Unload Response Ratio (LURR) signals are carried in order to verify statistical robustness of the previous studies using the Lattice Solid Model (MORA et al., 2002b). In each case 24 groups of samples with the same macroscopic parameters (tidal perturbation amplitude A, period T and tectonic loading rate k) but different particle arrangements are employed. Results of uni-axial compression experiments show that before the normalized time of catastrophic failure, the ensemble average LURR value rises significantly, in agreement with the observations of high LURR prior to the large earthquakes. In shearing tests, two parameters are found to control the correlation between earthquake occurrence and tidal stress. One is, A/(kT) controlling the phase shift between the peak seismicity rate and the peak amplitude of the perturbation stress. With an increase of this parameter, the phase shift is found to decrease. Another parameter, AT/k, controls the height of the probability density function (Pdf) of modeled seismicity. As this parameter increases, the Pdf becomes sharper and narrower, indicating a strong triggering. Statistical studies of LURR signals in shearing tests also suggest that except in strong triggering cases, where LURR cannot be calculated due to poor data in unloading cycles, the larger events are more likely to occur in higher LURR periods than the smaller ones, supporting the LURR hypothesis.

Spatio-temporal scanning and statistical test of the accelerating moment release (AMR) model using Australian earthquake data

The Accelerating Moment Release (AMR) preceding earthquakes with magnitude above 5 in Australia that occurred during the last 20 years was analyzed to test the Critical Point Hypothesis. Twelve earthquakes in the catalog were chosen based on a criterion for the number of nearby events. Results show that seven sequences with numerous events recorded leading up to the main earthquake exhibited accelerating moment release. Two occurred near in time and space to other earthquakes preceded by AM R. The remaining three sequences had very few events in the catalog so the lack of AMR detected in the analysis may be related to catalog incompleteness. Spatio-temporal scanning of AMR parameters shows that 80% of the areas in which AMR occurred experienced large events. In areas of similar background seismicity with no large events, 10 out of 12 cases exhibit no AMR, and two others are false alarms where AMR was observed but no large event followed. The relationship between AMR and Load-Unload Response Ratio (LURR) was studied. Both methods predict similar critical region sizes, however, the critical point time using AMR is slightly earlier than the time of the critical point LURR anomaly.