T288 Challenges in Tectonics: Fault Zone Behavior through Time, from Earth’s Surface to the Upper Mantle (Posters)

Tuesday, 24 October 2017
9:00 AM - 6:30 PM
Halls 4EF

USING MAXIMUM-LIKELIHOOD ESTIMATES AND BAYESIAN STATISTICS TO STUDY SHEAR ZONES AND RIDGE-TRANSFORM SYSTEMS

Paper number: T288-9
Program: Technical Sessions
Time: 9:00 AM – 6:30 PM

Authors

Abstract

There are two common approaches for quantitative modeling of shear zones and other fault-related systems. One approach is to construct dynamic models that are grounded in continuum mechanics. The modeler makes assumptions about geometry, rheology, etc., and then uses numerical methods to approximate the velocity and stress fields of the system. The other approach begins with the collection of field data, such as fabric orientations, deformed markers, and clast populations. The geologist builds a kinematic model and uses the field data to constrain its kinematic vorticity, amount of shortening, and other geometric parameters. Both methods have advantages and disadvantages. Dynamic models provide a wealth of physical information, but the results are often only loosely compared to field studies. Kinematic models often make use of large field data sets but provide only broad constraints on deformation.

We describe a statistically rigorous method that integrates numerical models and field data. The method relies on our ability to write down mathematical expressions for many different types of data. Following the classic dynamic modeling method, we choose a set of model assumptions and parameters, which gives us an approximation of the velocity field. By numerically integrating this velocity field, we are able to make predictions of structural data, which can then be compared to real field observations to build a likelihood function. Then, through maximum likelihood estimation and Bayesian Markov chain Monte Carlo simulation, we can deduce the values of the deformation parameters from the data.

In a series of numerical experiments, we have tested the method on an idealized transpressional shear zone. In some experiments, the velocity field is known analytically; in others, it is known only at discrete points in the domain, simulating output from a finite element method. Then, using a dynamic model of a ridge-transform intersection, we have applied the method to field data from a fossil ridge-transform system in the Troodos ophiolite in Cyprus.

Geological Society of America Abstracts with Programs. Vol. 49, No. 6

doi: 10.1130/abs/2017AM-299714

© Copyright 2017 The Geological Society of America (GSA), all rights reserved. Permission is hereby granted to the author(s) of this abstract to reproduce and distribute it freely, for noncommercial purposes. Permission is hereby granted to any individual scientist to download a single copy of this electronic file and reproduce up to 20 paper copies for noncommercial purposes advancing science and education, including classroom use, providing all reproductions include the complete content shown here, including the author information. All other forms of reproduction and/or transmittal are prohibited without written permission from GSA Copyright Permissions.