T111 From Outcrops to the Base of the Crust: Integrated Geophysical and Geological Illumination of Earth’s Tectonic Fabrics to Refine Understanding of Crustal Evolution and Natural Resources

Monday, 10 October 2022
1:30 PM - 5:30 PM
504

AN EIGENVALUE BASED K-MEANS APPROACH TO GEOLOGIC CLUSTER ANALYSIS

Paper number: T111-15
Program: Technical Sessions
Time: 4:45 PM – 5:00 PM

Authors

Abstract

Geologic directional data is often multimodal, and it is desirable to objectively partition it into clusters. These include fractures and faults, whose partitioning may assist in identifying overprinted deformation events, or coeval slip partitioning along a major fault system. Cluster analysis is a form of unsupervised learning, which includes parametric mixture and k-means models. In a k-means model each data point is assigned to one of k distributions based on a measure of distance. This nonparametric model is widely applied to geologic data.

The present approach maximizes the sum of an eigenvalue, or resultant vector, index over the clusters. This allows weighting, and is implemented for axis, vector, and girdle distributions. Cluster centers are seeded using random rotation matrices, from which iterative solutions converge by minimizing a cost function. This requires selection of the number of clusters, k, which can be subjective, so methods are given for the evaluation of k. A scalar distance-based minimization is also implemented, where cluster centers are found by minimizing the sum of spherical distances from the data. This is less flexible, but is more stable for large k and small n. Finally, a density-based cluster analysis is implemented.

Three data sets illustrate these methods. The first is fracture data from the San Manual copper mine, Arizona. Contouring the data on a Schmidt plot suggests k=3 clusters. The second are magnetic remanence vectors, which on theoretical grounds should have k=4 clusters. A third example is fault slip partitioning. Slip partitioning has been suggested to occur on transpressional fault systems, such as the Denali fault, and is a well established tenet of plate tectonic theory along extensional plate boundaries, such as the Mid-Atlantic Ridge (MAR). This example uses k=3 to distinguish three fault slip modes along the southern MAR using cluster analysis.

These implementations are available in the free Orient software.

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

doi: 10.1130/abs/2022AM-383667

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