Locating a strand of the Seattle Fault Zone in southeast Seattle, Washington through topographic analysis, borehole analysis and ground penetrating radar

Contributing to the evaluation of seismic hazards, a previously unmapped strand of the Seattle Fault Zone (SFZ), cutting across the southwest side of Lake Washington and southeast Seattle, is located and characterized on the basis of bathymetry, borehole logs, and ground penetrating radar (GPR). Previous geologic mapping and geophysical analysis of the Seattle area have generally mapped the locations of some strands of the SFZ, though a complete and accurate understanding of locations of all individual strands of the fault system is still incomplete. A bathymetric scarp-like feature and co-linear aeromagnetic anomaly lineament defined the extent of the study area. A 2-dimensional lithology cross-section was constructed using six boreholes, chosen from suitable boreholes in the study area. In addition, two GPR transects, oblique to the proposed fault trend, served to identify physical differences in subsurface materials. The proposed fault trace follows the previously mapped contact between the Oligocene Blakeley Formation and Quaternary deposits, and topographic changes in slope. GPR profiles in Seward Park and across the proposed fault location show the contact between the Blakeley Formation and unconsolidated glacial deposits, but it does not constrain an offset. However, north-dipping beds in the Blakely Formation are consistent with previous interpretations of P-wave seismic profiles on Mercer Island and Bellevue, Washington. The profiles show the mapped location of the aeromagnetic lineament in Lake Washington and the inferred location of the steeply-dipping, high-amplitude bedrock reflector, representing a fault strand. This north-dipping reflector is likely the same feature identified in my analysis. I characterize the strand as a splay fault, antithetic to the frontal fault of the SFZ. This new fault may pose a geologic hazard to the region.

River and Hillslope Response to Localized Uplift Along a Left Bend of the San Andreas Fault, Santa Cruz Mountains, CA

Studying landscape evolution of the Earthís surface is difficult because both tectonic forces and surface processes control its response to perturbation, and ultimately, its shape and form. Researchers often use numerical models to study erosional response to deformation because there are rarely natural settings in which we can evaluate both tectonic activity and topographic response over appropriate time scales (103-105 years). In certain locations, however, geologic conditions afford the unique opportunity to study the relationship between tectonics and topography. One such location is along the Dragonís Back Pressure Ridge in California, where the landscape moves over a structural discontinuity along the San Andreas Fault and landscape response to both the initiation and cessation of uplift can be observed. In their landmark study, Hilley and Arrowsmith (2008) found that geomorphic metrics such as channel steepness tracked uplift and that hillslope response lagged behind that of rivers. Ideal conditions such as uniform vegetation density and similar lithology allowed them to view each basin as a developmental stage of response to uplift only. Although this work represents a significant step forward in understanding landscape response to deformation, it remains unclear how these results translate to more geologically complex settings. In this study, I apply similar methodology to a left bend along the San Andreas Fault in the Santa Cruz Mountains, California. At this location, the landscape is translated through a zone of localized uplift caused by the bend, but vegetation, lithology, and structure vary. I examine the geomorphic response to uplift along the San Andreas Fault bend in order to determine whether predicted landscape patterns can be observed in a larger, more geologically complex setting than the Dragonís Back Pressure Ridge. I find that even with a larger-scale and a more complex setting, geomorphic metrics such as channel steepness index remain useful tools for evaluating landscape evolution through time. Steepness indices in selected streams of study record localized uplift caused by the restraining bend, while hillslope adjustment in the form of landsliding occurs over longer time scales. This project illustrates that it is possible to apply concepts of landscape evolution models to complex settings and is an important contribution to the body of geomorphological study.

Automatic message handling for a national counselling service

An automatic email handling system (AutoRouter) was introduced at a national counselling service in Australia. In 2003, counsellors responded to a total of 7421 email messages. Over nine days in early May 2004 the administrator responsible for the management of the manual email counselling service recorded the time spent on managing email messages. The AutoRouter was then introduced. Since the implementation of the AutoRouter the administrator's management role has become redundant, an average of 12 h 5 min per week of staff time has been saved. There have been further savings in supervisor time. Counsellors were taking an average of 6.2 days to respond to email messages (n=4307), with an average delay of 1.2 days from the time counsellors wrote the email to when the email was sent. Thus the response was sent on average 7.4 days after receipt of the original client email message. A significant decrease in response time has been noted since implementation of the AutoRouter, with client responses now taking an average of 5.4 days, a decrease of 2.0 days. Automatic message handling appears to be a promising method of managing the administration of a steadily increasing email counselling service.

On the possibility of elastic strain localisation in a fault

The phenomenon of strain localisation is often observed in shear deformation of particulate materials, e.g., fault gouge. This phenomenon is usually attributed to special types of plastic behaviour of the material (e.g., strain softening or mismatch between dilatancy and pressure sensitivity or both). Observations of strain localisation in situ or in experiments are usually based on displacement measurements and subsequent computation of the displacement gradient. While in conventional continua the symmetric part of the displacement gradient is equal to the strain, it is no longer the case in the more realistic descriptions within the framework of generalised continua. In such models the rotations of the gouge particles are considered as independent degrees of freedom the values of which usually differ from the rotation of an infinitesimal volume element of the continuum, the latter being described for infinitesimal deformations by the non-symmetric part of the displacement gradient. As a model for gouge material we propose a continuum description for an assembly of spherical particles of equal radius in which the particle rotation is treated as an independent degree of freedom. Based on this model we consider simple shear deformations of the fault gouge. We show that there exist values of the model parameters for which the displacement gradient exhibits a pronounced localisation at the mid-layers of the fault, even in the absence of inelasticity. Inelastic effects are neglected in order to highlight the role of the independent rotations and the associated additional parameters. The localisation-like behaviour occurs if (a) the particle rotations on the boundary of the shear layer are constrained (this type of boundary condition does not exist in a standard continuum) and (b) the contact moment-or bending stiffness is much smaller than the product of the effective shear modulus of the granulate and the square of the width of the gouge layer. It should be noted however that the virtual work functional is positive definite over the range of physically meaningful parameters (here: contact stiffnesses, solid volume fraction and coordination number) so that strictly speaking we are not dealing with a material instability.

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.

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.

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.

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.

Postharvest characteristics and handling of litchi fruit - an overview

Litchi ( Litchi chinensis Sonn.) is a tropical to subtropical crop that originated in South-East Asia. Litchi fruit are prized on the world market for their flavour, semi-translucent white aril and attractive red skin. Litchi is now grown commercially in many countries and production in Australia, China, Israel, South Africa and Thailand has expanded markedly in recent years. Increased production has made significant contributions to economic development in these countries, especially those in South-East Asia. Non-climacteric litchi fruit are harvested at their visual and organoleptic optimum. They are highly perishable and, consequently, have a short life that limits marketability and potential expansion of demand. Pericarp browning and pathological decay are common and important defects of harvested litchi fruit. Postharvest technologies have been developed to reduce these defects. These technologies involve cooling and heating the fruit, use of various packages and packaging materials and the application of fungicides and other chemicals. Through the use of fungicides and refrigeration, litchi fruit have a storage life of about 30 days. However, when they are removed from storage, their shelf life at ambient temperature is very short due to pericarp browning and fruit rotting. Low temperature acclimation or use of chitsoan as a coating can extend the shelf life. Sulfur dioxide fumigation effectively reduces pericarp browning, but approval from Europe, Australia and Japan for this chemical is likely to be withdrawn due to concerns over sulfur residues in fumigated fruit. Thus, sulfur-free postharvest treatments that maintain fruit skin colour are increasingly important. Alternatives to SO2 fumigation for control of pericarp browning and fruit rotting are pre-storage pathogen management, anoxia treatment, and dipping in 2% hydrogen chloride solution for 6-8 min following storage at 0 degrees C. Insect disinfestation has become increasingly important for the expansion of export markets because of quarantine issues associated with some fruit fly species. Thus, effective disinfestation protocols need to be developed. Heat treatment has shown promise as a quarantine technology, but it injures pericarp tissue and results in skin browning. However, heat treatment can be combined with an acid dip treatment that inhibits browning. Therefore, the primary aim of postharvest litchi research remains the achievement of highly coloured fruit which is free of pests and disease. Future research should focus on disease control before harvest, combined acid and heat treatments after harvest and careful temperature management during storage and transport.

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.