HR: 0800h
AN: NG41B-0516    [Abstracts]
TI: Virtual California: Interplay Between Fault System Complexity, Dynamics, and Numerical Stability
AU: * Rundle, P
EM: p.b.rundle@gmail.com
AF: Stanford University, Department of Geophysics, Stanford, CA 94305, United States
AU: Morein, G
EM: gleb@cse.ucdavis.edu
AF: University of California, Department of Physics One Shields Ave, Davis, CA 95616, United States
AU: Rundle, J
EM: jbrundle@ucdavis.edu
AF: University of California, Department of Physics One Shields Ave, Davis, CA 95616, United States
AU: Tiampo, K
EM: ktiampo@uwo.ca
AF: University of Western Ontario, Department of Earth Sciences, London, ON N6A 5B7, Canada
AU: Donnellan, A
EM: andrea.donnellan@jpl.nasa.gov
AF: Jet Propulsion Laboratory, Division of Earth and Space Science, Pasadena, CA 91109, United States
AU: Turcotte, D
EM: turcotte@geology.ucdavis.edu
AF: University of California, Department of Geology One Shields Ave, Davis, CA 95616, United States
AB: Virtual California is a topologically realistic simulation of the interacting earthquake faults in California. Inputs to the model arise from field data, and typically include realistic fault system topologies, realistic long term slip rates, and realistic frictional parameters. Outputs from the simulations include synthetic earthquake sequences and space-time patterns together with associated surface deformation and strain patterns that are similar to those seen in nature. It is a type of model called a backslip model, so called because the loading in the model arises from negative slip or "backslip" applied to each fault boundary element at its geologically observed long- term rate of offset, V(x). We note that in this context, a fault element is regarded simply as a "degree of freedom", rather than a spatially coherent entity with geological meaning. All of the rectangular fault elements are embedded in an elastic half space, and they interact with each other by means of quasistatic elastic interactions, whose stress Greens functions are computed by means of a boundary element method. Frictional coefficients are assigned to each fault element, along with other frictional parameters, by means of a data assimilation technique. When the model is used to produce a simulation, the result is a history of slip on the fault elements in response to the driving forces. The interactions between the fault elements serve to organize the system so that, instead of a sequence of single elements breaking individually, multiple elements break simultaneously, producing large earthquakes. The dynamical evolution in the model is also produced by means of a stochastic, cellular automaton method, in which a random overshoot or undershoot component is added at the time of sliding of a element. The topology of the model fault system has evolved over time, becoming more complex, and using progressively finer scales of detail. Here we discuss the latest versions of the fault system topology and the associated numerical simulations. We focus in particular on the interplay between fault system complexity, and complexity in the dynamics. We find that the mesh scale of the boundary elements is strongly affected by the scale fault topology, and that the stability properties of the dynamics are likewise strongly affected.
DE: 4430 Complex systems
DE: 4435 Emergent phenomena
DE: 4465 Phase transitions
DE: 5104 Fracture and flow
SC: Nonlinear Geophysics [NG]
MN: 2007 Fall Meeting