T142 Structural Geology in the 21st Century: Rheology of Naturally Deformed Rocks

Monday, 23 September 2019
1:30 PM - 5:30 PM
Room 229A, North Building

CALCULATION OF THE BEST-FIT HOMOGENEOUS STRAIN FROM A DISPLACEMENT FIELD USING DISPLACEMENT STRETCH PLOT AND STRAIN PROBE THEORY

Paper number: T142-8
Program: Technical Sessions
Time: 3:15 PM – 3:30 PM

Authors

Abstract

Studies of glacial flow, continental deformation, slope creep, and analog experiments use the positions of survey markers at successive intervals in time to determine displacement fields. Relationships between displacements and strain are given by continuum mechanics (Mase 1970; Ramsay 1967; Elliott 1970; Means 1976; Robin 2019).

A "displacement stretch plot", implemented in FORTRAN by Vollmer (in Means 1983) for analog experiments, is quantified here. The plot is created by calculating stretches and final orientations of the line segments defined by combinations of n markers. For i in [1..n], initial positions are [xi, yi], and final positions are [x'i, y'i], giving c = n!/(2(n-2)!) line segments. For k in [1..c] the initial lengths are given by dk² = dxk² + dyk², final lengths by d'k² = dx'k² + dy'k², stretches Sk = d'k/dk, and final angles α'k = atan(dy'k/dx'k). A symmetric polar plot of the stretches defines the estimated strain ellipse. The best-fit ellipse is calculated by minimizing the sum of the squares of the radial residuals (Fitzgibbon et al. 1996; Halır and Flusser 1998), and solved for strain axes and dilation, Δ.

Robin (2019) proposed a "strain probe" method to estimate the deformation tensor U, and displacement gradient tensor, E, from the initial and final positions. For finite strain, U can be decomposed into a spin tensor R and a symmetric deformation tensor. Left decomposition gives the ellipse in a Lagrangian reference frame, the spin tensor R, rotation angle ω, and symmetric tensor UL. Eigenvectors of UL give the principal strains. A coefficient of determination, CD on [0..1], is calculated to evaluate the solution.

Data from an analogue experiment by Schweinberger and Fueten is used to compare results. Displacement stretch plot: A = 1.0914, B = 0.9749, R = 1.1194, ϕ = 28.155°, Δ = 0.0640, CD = 0.9996. Strain probe: A = 1.0943, B = 0.9680, R = 1.1304, ϕ = 28.623°, Δ = 0.0592, CD = 0.9909. The methods are implemented in EllipseFit 3.7.0.

Geological Society of America Abstracts with Programs. Vol. 51, No. 5

doi: 10.1130/abs/2019AM-340051

© Copyright 2019 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.