Modelling of deformation around magmatic intrusions with application to gold-related structures in the Yilgarn Craton, Western Australia

Abstract

This study simulates rock deformation around high temperature granite intrusions and explores how gold bearing shear zones near intrusions were developed in the Yilgarn, using a new continuum damage mechanics algorithm that considers the temperature and time dependent elastic-visco-plastic constitutive behaviour of crustal materials. The results demonstrate that strain rates have the most significant effects on structural patterns for both extensional and compressional cases. Smaller strain rates promote the formation of narrow high-strain shear zones and strong strain localisation along the flank or shoulder areas of the intrusion and cold granite dome. Wider diffuse shear zones are developed under higher strain rates due to strain hardening. The cooling of the intrusion to background temperatures occurred over a much shorter time interval when compared to the duration of deformation and shear zones development. Strong strain localisation near the intrusion and shear zone development in the crust occurred under both extensional and compressional conditions. There is always clear strain localisation around the shoulders of the intrusion and the flanks of the “cold” granitic dome in early deformation stages. In the models containing a pre-existing fault, strain localisation near the intrusion became asymmetric with much stronger localisation and the development of a damage zone at the shoulder adjacent to the reactivated fault. At higher deformation stages, the models produced a range of structural patterns including graben and half graben basin (extension), “pop-up” wedge structures (compression), tilted fault blocks and switch of shear movement from reverse to normal on shear zones. The model explains in part why a number of gold deposits (e.g. Wallaby and Paddington deposits) in the Yilgarn were formed near the flank of granite-cored domes and deep “tapping” faults, and shows that the new modelling approach is capable of realistically simulating high strain localisation and shear zone development.