Microstructural Analysis for Dynamic Pulverization and Asymmetric Damage at the Base of Seismogenic Strike-slip Faults

Although the mechanics of continental, seismogenic strike-slip faults have been primarily studied around active faults near Earth’s surface, large earthquakes on these faults commonly extend to depths between 10 and 20 km. At the base of seismogenic strike-slip faults, interaction and feedback between coseismic brittle fracturing and post- and interseismic viscous flow affect transient and long-term changes in stress cycling, fluid and heat transport, fault strength, and associated strain localization and deformation mechanisms. A primary goal of my dissertation is to explore the deeper structures of damage zones near the base of the seismogenic zone and to better understand the influence of the damaged rocks on rupture dynamics, by examining microstructures of exhumed fault rocks. My study area, the Sandhill Corner shear zone that is the longest strand of the Paleozoic Norumbega fault system in Maine, USA, represents large-displacement, seismogenic strike-slip faults at frictional-to-viscous transition depths (corresponding to temperatures of ~400–500 °C). The shear zone contains mutually overprinting pseudotachylyte and mylonite, and juxtaposes quartzofeldspathic mylonites and mica-rich schists. I analyzed fractured and fragmented garnet grains using particle size distributions, microfracture patterns, and electron backscatter diffraction fabrics. Microstructural studies of fragmented garnets reveal asymmetric distribution of dynamic pulverization with a width of ~70 m in the Sandhill Corner shear zone, and these results imply that the same damage processes observed around active seismogenic strike-slip faults operate at the base of the seismogenic zone. Garnet microstructures formed during earthquake cycles at the frictional-viscous transition can also provide evidence for dynamic pulverization even though the particle size distribution is modified by quasi-static fragmentation during post- and interseismic shearing. Elastic and seismic properties of the quartzofeldspathic rock and the mica-rich schist are quantified using the Thermo-Elastic and Seismic Analysis (TESA) numerical toolbox. The results illustrate how elastic contrast across bimaterial faults separating two different anisotropic materials affects preferred rupture propagation and asymmetric damage distribution. Strong anisotropy occurs in fault zones where preferentially aligned phyllosilicate minerals are a major component of the modal mineralogy. My findings suggest that the orientation and proportion of preferentially aligned phyllosilicates, or other highly anisotropic minerals, should be considered when investigating fault ruptures in anisotropic rocks

Ground movement associated with microtunneling

Microtunneling is a trenchless technology for construction of pipelines. Its process is a cyclic pipe jacking operation. Microtunneling has been typically used for gravity sewer systems in urban areas. Despite its good success record overall, several large ground settlement cases caused by microtunneling have been reported. Also, in contrast with large diameter urban tunneling, there are few research projects about the ground settlement caused by microtunneling. In this dissertation, the ground settlement caused by microtunneling is studied using a theoretical approach, empirical approach, numerical simulation approach, and artificial intelligence approach. In the theoretical approach, the equivalent ground loss and settlement caused by concentrated ground loss have been used to drive the ground settlement profile. In the empirical approach, the ground settlement caused by large diameter tunneling case histories is used. In the numerical approach, FLAC 3D software, a commercially available finite difference code, is used to simulate the ground settlement caused by microtunneling. In the artificial intelligence approach, a three-layer back propagation neural network is developed to predict the ground settlement caused by microtunneling using the numerical simulation results. It is found that the neural network developed as part of this thesis work provides a means of rapid prediction of the surface ground settlement curve based on the soil parameters, project geometry and estimated ground loss. This prediction matches FLAC3D results very well over the full range of parameters studied and has a reasonable correspondence to the field results with which it was compared

A Second Replicated Quantitative Analysis of Fault Distributions in Complex Software Systems

Background. Software engineering is in search for general principles that apply across contexts, for example to help guide software quality assurance. Fenton and Ohlsson presented such observations on fault distributions, which have been replicated once. Objectives.We intend to replicate their study a second time in a new environment. Method.We conducted a close replication, collecting defect data from five consecutive releases of a large software system in the telecommunications domain, and conducted the same analysis as in the original study. Results. The replication confirms results on un-evenly distributed faults over modules, and that fault proneness distribution persist over test phases. Size measures are not useful as predictors of fault proneness, while fault densities are of the same order of magnitude across releases and contexts. Conclusions. This replication confirms that the un-even distribution of defects motivates un-even distribution of quality assurance efforts, although predictors for such distribution of efforts are not sufficiently precise

Strain Partitioning and Frictional Behavior of Opalinus Clay During Fault Reactivation

The Opalinus Clay (OPA) formation is considered a suitable host rock candidate for nuclear waste storage. However, the sealing integrity and long-term safety of OPA are potentially compromised by pre-existing natural or artificially induced faults. Therefore, characterizing the mechanical behavior and microscale deformation mechanisms of faults and the surrounding rock is relevant for predicting repository damage evolution. In this study, we performed triaxial tests using saw-cut samples of the shaly and sandy facies of OPA to investigate the influence of pressure and mineral composition on the deformation behavior during fault reactivation. Dried samples were hydrostatically pre-compacted at 50 MPa and then deformed at constant strain rate, drained conditions and confining pressures (pc) of 5–35 MPa. Mechanical data from triaxial tests was complemented by local strain measurements to determine the relative contribution of bulk deformation and fault slip, as well as by acoustic emission (AE) monitoring, and elastic P-wave velocity measurements using ultrasonic transmissions. With increasing pc, we observe a transition from brittle deformation behavior with highly localized fault slip to semi-brittle behavior characterized by non-linear strain hardening with increasing delocalization of deformation. We find that brittle localization behavior is limited by pc at which fault strength exceeds matrix yield strength. AEs were only detected in tests performed on sandy facies samples, and activity decreased with increasing pc. Microstructural analysis of deformed samples revealed a positive correlation between increasing pc and gouge layer thickness. This goes along with a change from brittle fragmentation and frictional sliding to the development of shear zones with a higher contribution of cataclastic and granular flow. Friction coefficient at fault reactivation is only slightly higher for the sandy (µ ~ 0.48) compared to the shaly facies (µ ~ 0.4). Slide-hold-slide tests performed after ~ 6 mm axial shortening suggest stable creeping and long-term weakness of faults at the applied conditions. Our results demonstrate that the mode of fault reactivation highly depends on the present stress field and burial history

Geomechanical modelling of sinkhole development using distinct elements: model verification for a single void space and application to the Dead Sea area

Mechanical and/or chemical removal of material from the subsurface may generate large subsurface cavities, the destabilisation of which can lead to ground collapse and the formation of sinkholes. Numerical simulation of the interaction of cavity growth, host material deformation and overburden collapse is desirable to better understand the sinkhole hazard but is a challenging task due to the involved high strains and material discontinuities. Here, we present 2-D distinct element method numerical simulations of cavity growth and sinkhole development. Firstly, we simulate cavity formation by quasi-static, stepwise removal of material in a single growing zone of an arbitrary geometry and depth. We benchmark this approach against analytical and boundary element method models of a deep void space in a linear elastic material. Secondly, we explore the effects of properties of different uniform materials on cavity stability and sinkhole development. We perform simulated biaxial tests to calibrate macroscopic geotechnical parameters of three model materials representative of those in which sinkholes develop at the Dead Sea shoreline: mud, alluvium and salt. We show that weak materials do not support large cavities, leading to gradual sagging or suffusion-style subsidence. Strong materials support quasi-stable to stable cavities, the overburdens of which may fail suddenly in a caprock or bedrock collapse style. Thirdly, we examine the consequences of layered arrangements of weak and strong materials. We find that these are more susceptible to sinkhole collapse than uniform materials not only due to a lower integrated strength of the overburden but also due to an inhibition of stabilising stress arching. Finally, we compare our model sinkhole geometries to observations at the Ghor Al-Haditha sinkhole site in Jordan. Sinkhole depth ∕ diameter ratios of 0.15 in mud, 0.37 in alluvium and 0.33 in salt are reproduced successfully in the calibrated model materials. The model results suggest that the observed distribution of sinkhole depth ∕ diameter values in each material type may partly reflect sinkhole growth trends

Thermo-hydro-mechanical simulation of a generic geological disposal facility for radioactive waste

Geological disposal is required for the safe and long-term disposal of legacy radioactive waste. High level waste and spent fuel generate significant heat that will cause thermo-hydro-mechanical coupled processes in the rock mass. The thermal expansion of the fluid will be greater than the grains causing a decrease in mean effective stress with the low permeability restricting Darcy flow and excess pore pressure equilibration. A decrease in mean effective stress can reduce material strength in granular materials, which may be significant near excavations where differential stress is increased. Microseismic monitoring provides cost effective, non-intrusive and three-dimensional data that can be calibrated with the stress and strain behaviour of a rock mass. However, there is no precedent for the microseismic monitoring of heat-producing radioactive waste. Generic concepts, analogue materials and data from in situ experiments are used to demonstrate the potential for the microseismic monitoring of heat-producing radioactive waste in lower strength sedimentary rocks. A mechanism for early post-closure microseismicity is demonstrated, whereby excess pore pressure decreases the mean effective stress towards yielding in shear. The rock and fluid property uncertainties are ranked according to their contribution to the excess pore pressure. Permeability is found to be important as expected, however, Biot's coefficient is demonstrably more important and yet often overlooked. Furthermore, the microseismic event locations, timings and pseudo scalar seismic moments are shown to have statistically significant relationships with the engineered backfill swelling pressure. Therefore, early post-emplacement microseismic monitoring could provide constraints for the engineered backfill swelling pressure and rock property uncertainties whilst the facility is still operational. Insights could prove timely for adapting the engineering designs, if they are not behaving as expected, in further high level waste and spent fuel tunnels

Identification of fault and top seal effectiveness through an integration of hydrodynamic and capillary analysis techniques

Fault and top seal effectiveness has proved to be a significant risk in exploration success, and creates a large uncertainty in predicting reservoir performance. This is particularly true in the Australian context, but equally applies to exploration provinces worldwide. Seals can be broadly classified into fault, intraformational, and top seal. For geological time-scale processes, intraformational and top seals are typically characterised by their membrane seal capacity and fracture threshold pressure. Fault seals are typically characterised by fault geometry, juxtaposition, membrane seal capacity, and reactivation potential. At the production time scale, subtle variations in the permeability distribution within a reservoir can lead to compartmentalization. These are typically characterised by dynamic reservoir models which assume hydrostatic conditions prior to commencement of production. There are few references in the seals literature concerning the integration of hydrodynamic techniques with the various aspects of seal evaluation. The research for this PhD thesis by published papers includes: Methodology for characterising formation water flow systems in faulted strata at exploration and production time scales; a new theory of hydrodynamics and membrane (capillary) seal capacity; and case study evaluations demonstrating integrated multidisciplinary techniques for the evaluation of seal capacity (fault, intraformational and top seal) that demonstrate the new theory in practice. By incorporating hydrodynamic processes in the evaluation of total seal capacity, the evidence shows that existing shale gouge ratio – across fault pressure difference (SGR-AFPD) calibration plots need adjustment resulting in the calibration envelopes shifting to the centre of the plot.This adjustment sharpens the predictive capacity for membrane seal analysis in the pre-drill scenario. This PhD thesis presents the background and rationale for the thesis topic, presents each published paper to be included as part of the thesis and its contribution to the body of work addressing the thesis topic, and presents related published papers that are not included in the thesis but which support the body of published work on the thesis topic. The result of the thesis is a new theory and approach to characterising membrane seal capacity for the total seal thickness, and has implications for an adjusted SGR-AFPD calibration to be applied in pre-drill evaluations of seal capacity. A large portion of the resources and data required to conduct the research were made available by CSIRO and its associated project sponsors including the CO2CRC

Large Impact Craters and Basins: Mechanics of Syngenetic and Postgenetic Modification

The impact crater is the ubiquitous landform of the solar system. Theoretical mechanical analyses are applied to the modification stage of crater formation, both syngenetic (immediate or short term) and postgenetic (long term). The mechanical stability of an impact crater is analyzed via a quasi-static, axisymmetric slip line theory of plasticity. The yield model incorporated is Mohr-Coulomb and a simplified rectangular profile is used for the transient cavity. The degree of stability (or instability) is described as a function of internal friction angle, depth/diameter ratio, and a dimensionless parameter ρgH/c (ρ = density, g = acceleration of gravity, H = depth, and c = cohesion strength). To match the observed slumping of large lunar craters the cohesion strength of the lunar surface material must be low (&lt;20 bars) and the angle of internal friction must be less than 2°. It is not implausible that these failure strength characteristics are realized by freshly shocked rock. A theoretical description of impact crater collapse is evolved which accounts for the development of wall scallops, slump terraces, and flat floors. A preliminary set of scale model experiments performed in a centrifuge corroborate the theory. The strength of terrestrial planet surfaces under impact is seen to vary by as much as a factor of two. Shortly after the excavation of a large impact crater the transient cavity collapses, driven by gravity. It is shown that at least one concentric fault scarp forms about the crater, if the strength of the target material decreases sufficiently rapidly with increasing depth. This is demonstrated by two classes of models: extrusion flow models which assume a weak layer underlying a strong layer, and plastic flow models which assume a continuous decrease of cohesion strength with depth. Both classes predict that the ratio of the radius of the scarp to the transient crater radius is between 1.2 and 2 for large craters. Large impact basins on Ganymede and Callisto are characterized by one or more concentric rings or scarps. The number, spacing, and morphology of the rings is a function of the thickness and strength of the lithosphere, and crater diameter. When the lithosphere is thin and weak, the collapse is regulated by flow induced in the asthenosphere. The lithosphere fragments in a multiply concentric pattern (e.g., Valhalla, Asgard, Galilee Regio, and a newly discovered ring system on Callisto). The thickness and viscosity of a planetary lithosphere increases with time as the mantle cools. A thicker lithosphere leads to the formation of one (or very few) irregular normal faults concentric to the crater (e.g., Gilgamesh). A gravity wave or tsunami induced by impact into a liquid mantle would result in both concentric and radial extension features. Since these are not observed , this process cannot be responsible for the generation of the rings around the basins on Ganymede and Callisto. The appearance of Galilee Regio and portions of Valhalla is best explained by ring graben, and though the Valhalla system is older, the lithosphere was 1.5-2.0 times as thick at the time of formation. The present lithosphere thickness is too great to permit development of any rings. It has been proposed that a mascon may be in the form of an annulus surrounding the Caloris basin on Mercury, associated with the smooth plains. The effects (stresses, deformation, surface tectonic style, gravity anomalies, etc.) of such a ring load on a floating elastic lithosphere of variable thickness are investigated. The main characteristics of the surface tectonic pattern are normal faulting within the basin and thrust faulting beneath the ring load both in agreement with observation Moreover, the dominant concentric trend of the basin normal faults is consistent with the ring load hypothesis provided the mercurian lithosphere was ≤125 km thick at the time of faulting. Simple updoming within the basin would produce normal faults of predominantly radial orientation.</p