Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Pure & Applied Geophysics”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5Linked to original sources
USGS2006

Puhimau thermal area: a window into the upper east rift zone of Kilauea Volcano, Hawaii?

We report the results of two soil CO 2 efflux surveys by the closed chamber circulation method at the Puhimau thermal area in the upper East Rift Zone (ERZ) of Kilauea volcano, Hawaii. The surveys were undertaken in 1996 and 1998 to constrain how much CO 2 might be reaching the ERZ after degassing beneath the summit caldera and whether the Puhimau thermal area might be a significant contributor to the overall CO 2 budget of Kilauea. The area was revisited in 2001 to determine the effects of surface disturbance on efflux values by the collar emplacement technique utilized in the earlier surveys. Utilizing a cutoff value of 50 g m −2 d −1 for the surrounding forest background efflux, the CO 2 emission rates for the anomaly at Puhimau thermal area were 27 t d −1 in 1996 and 17 t d −1 in 1998. Water vapor was removed before analysis in all cases in order to obtain CO 2 values on a dry air basis and mitigate the effect of water vapor dilution on the measurements. It is clear that Puhimau thermal area is not a significant contributor to Kilauea's CO 2 output and that most of Kilauea's CO 2 (8500 t d −1 ) is degassed at the summit, leaving only magma with its remaining stored volatiles, such as SO 2 , for injection down the ERZ. Because of the low CO 2 emission rate and the presence of a shallow water table in the upper ERZ that effectively scrubs SO 2 and other acid gases, Puhimau thermal area currently does not appear to be generally well suited for observing temporal changes in degassing at Kilauea.

Hawaii↗
USGS2007

Generating an image of dispersive energy by frequency decomposition and slant stacking

We present a new algorithm for calculating an image of dispersive energy in the frequency-velocity (f-v) domain. The frequency decomposition is first applied to a shot gather in the offset-time domain to stretch impulsive data into pseudo-vibroseis data or frequency-swept data. Because there is a deterministic relationship between frequency and time in a sweep used in the frequency decomposition, the first step theoretically completes the transform from time to frequency. The slant stacking is then performed on the frequency-swept data to complete the transform from offset to velocity. This simple two-step algorithm generates an image of dispersive energy in the f-v domain. The straightforward transform only uses offset information of data so that this algorithm can be applied to data acquired with arbitrary geophone-acquisition geometry. Examples of synthetic and real-world data demonstrate that this algorithm generates accurate images of dispersive energy of the fundamental as well as higher modes. ?? Birkha??user Verlag, Basel, 2007.

Pure and Applied Geophysics↗
USGS2008

Rayleigh-wave dispersive energy imaging using a high-resolution linear radon transform

Multichannel Analysis of Surface Waves (MASW) analysis is an efficient tool to obtain the vertical shear-wave profile. One of the key steps in the MASW method is to generate an image of dispersive energy in the frequency-velocity domain, so dispersion curves can be determined by picking peaks of dispersion energy. In this paper, we propose to image Rayleigh-wave dispersive energy by high-resolution linear Radon transform (LRT). The shot gather is first transformed along the time direction to the frequency domain and then the Rayleigh-wave dispersive energy can be imaged by high-resolution LRT using a weighted preconditioned conjugate gradient algorithm. Synthetic data with a set of linear events are presented to show the process of generating dispersive energy. Results of synthetic and real-world examples demonstrate that, compared with the slant stacking algorithm, high-resolution LRT can improve the resolution of images of dispersion energy by more than 50%. ?? Birkhaueser 2008.

Pure and Applied Geophysics↗
USGS2009

The November 15, 2006 Kuril Islands-generated tsunami in Crescent City, California

On November 15, 2006, Crescent City in Del Norte County, California was hit by a tsunami generated by a M w 8.3 earthquake in the central Kuril Islands. Strong currents that persisted over an eight-hour period damaged floating docks and several boats and caused an estimated $9.2 million in losses. Initial tsunami alert bulletins issued by the West Coast Alaska Tsunami Warning Center (WCATWC) in Palmer, Alaska were cancelled about three and a half hours after the earthquake, nearly five hours before the first surges reached Crescent City. The largest amplitude wave, 1.76-meter peak to trough, was the sixth cycle and arrived over two hours after the first wave. Strong currents estimated at over 10 knots, damaged or destroyed three docks and caused cracks in most of the remaining docks. As a result of the November 15 event, WCATWC changed the definition of Advisory from a region-wide alert bulletin meaning that a potential tsunami is 6 hours or further away to a localized alert that tsunami water heights may approach warning- level thresholds in specific, vulnerable locations like Crescent City. On January 13, 2007 a similar Kuril event occurred and hourly conferences between the warning center and regional weather forecasts were held with a considerable improvement in the flow of information to local coastal jurisdictions. The event highlighted the vulnerability of harbors from a relatively modest tsunami and underscored the need to improve public education regarding the duration of the tsunami hazards, improve dialog between tsunami warning centers and local jurisdictions, and better understand the currents produced by tsunamis in harbors.

California↗
USGS2009

Approximation to cutoffs of higher modes of Rayleigh waves for a layered earth model

A cutoff defines the long-period termination of a Rayleigh-wave higher mode and, therefore is a key characteristic of higher mode energy relationship to several material properties of the subsurface. Cutoffs have been used to estimate the shear-wave velocity of an underlying half space of a layered earth model. In this study, we describe a method that replaces the multilayer earth model with a single surface layer overlying the half-space model, accomplished by harmonic averaging of velocities and arithmetic averaging of densities. Using numerical comparisons with theoretical models validates the single-layer approximation. Accuracy of this single-layer approximation is best defined by values of the calculated error in the frequency and phase velocity estimate at a cutoff. Our proposed method is intuitively explained using ray theory. Numerical results indicate that a cutoffs frequency is controlled by the averaged elastic properties within the passing depth of Rayleigh waves and the shear-wave velocity of the underlying half space. ?? Birkh??user Verlag, Basel 2009.

Pure and Applied Geophysics↗
USGS2009

Geometry of the Nojima fault at Nojima-Hirabayashi, Japan - II. Microstructures and their implications for permeability and strength

Samples of damage-zone granodiorite and fault core from two drillholes into the active, strike-slip Nojima fault zone display microstructures and alteration features that explain their measured present-day strengths and permeabilities and provide insight on the evolution of these properties in the fault zone. The least deformed damage-zone rocks contain two sets of nearly perpendicular (60–90° angles), roughly vertical fractures that are concentrated in quartz-rich areas, with one set typically dominating over the other. With increasing intensity of deformation, which corresponds generally to increasing proximity to the core, zones of heavily fragmented rock, termed microbreccia zones, develop between prominent fractures of both sets. Granodiorite adjoining intersecting microbreccia zones in the active fault strands has been repeatedly fractured and locally brecciated, accompanied by the generation of millimeter-scale voids that are partly filled with secondary minerals. Minor shear bands overprint some of the heavily deformed areas, and small-scale shear zones form from the pairing of closely spaced shear bands. Strength and permeability measurements were made on core collected from the fault within a year after a major (Kobe) earthquake. Measured strengths of the samples decrease regularly with increasing fracturing and fragmentation, such that the gouge of the fault core and completely brecciated samples from the damage zone are the weakest. Permeability increases with increasing disruption, generally reaching a peak in heavily fractured but still more or less cohesive rock at the scale of the laboratory samples. Complete loss of cohesion, as in the gouge or the interiors of large microbreccia zones, is accompanied by a reduction of permeability by 1-2 orders of magnitude below the peak values. The core samples show abundant evidence of hydrothermal alteration and mineral precipitation. Permeability is thus expected to decrease and strength to increase somewhat in active fault strands between earthquakes, as mineral deposits progressively seal fractures and fill pore spaces.

Pure and Applied Geophysics↗
USGS2009

Geometry of the Nojima fault at Nojima-Hirabayashi, Japan - I. A simple damage structure inferred from borehole core permeability

The 1995 Kobe (Hyogo-ken Nanbu) earthquake, M = 7.2, ruptured the Nojima fault in southwest Japan. We have studied core samples taken from two scientific drillholes that crossed the fault zone SW of the epicentral region on Awaji Island. The shallower hole, drilled by the Geological Survey of Japan (GSJ), was started 75 m to the SE of the surface trace of the Nojima fault and crossed the fault at a depth of 624 m. A deeper hole, drilled by the National Research Institute for Earth Science and Disaster Prevention (NIED) was started 302 m to the SE of the fault and crossed fault strands below a depth of 1140 m. We have measured strength and matrix permeability of core samples taken from these two drillholes. We find a strong correlation between permeability and proximity to the fault zone shear axes. The half-width of the high permeability zone (approximately 15 to 25 m) is in good agreement with the fault zone width inferred from trapped seismic wave analysis and other evidence. The fault zone core or shear axis contains clays with permeabilities of approximately 0.1 to 1 microdarcy at 50 MPa effective confining pressure (10 to 30 microdarcy at in situ pressures). Within a few meters of the fault zone core, the rock is highly fractured but has sustained little net shear. Matrix permeability of this zone is approximately 30 to 60 microdarcy at 50 MPa effective confining pressure (300 to 1000 microdarcy at in situ pressures). Outside this damage zone, matrix permeability drops below 0.01 microdarcy. The clay-rich core material has the lowest strength with a coefficient of friction of approximately 0.55. Shear strength increases with distance from the shear axis. These permeability and strength observations reveal a simple fault zone structure with a relatively weak fine-grained core surrounded by a damage zone of fractured rock. In this case, the damage zone will act as a high-permeability conduit for vertical and horizontal flow in the plane of the fault. The fine-grained core region, however, will impede fluid flow across the fault.

Pure and Applied Geophysics↗
USGS2009

Constructing constitutive relationships for seismic and aseismic fault slip

For the purpose of modeling natural fault slip, a useful result from an experimental fault mechanics study would be a physically-based constitutive relation that well characterizes all the relevant observations. This report describes an approach for constructing such equations. Where possible the construction intends to identify or, at least, attribute physical processes and contact scale physics to the observations such that the resulting relations can be extrapolated in conditions and scale between the laboratory and the Earth. The approach is developed as an alternative but is based on Ruina (1983) and is illustrated initially by constructing a couple of relations from that study. In addition, two example constitutive relationships are constructed; these describe laboratory observations not well-modeled by Ruina's equations: the unexpected shear-induced weakening of silica-rich rocks at high slip speed (Goldsby and Tullis, 2002) and fault strength in the brittle ductile transition zone (Shimamoto, 1986). The examples, provided as illustration, may also be useful for quantitative modeling.

Pure and Applied Geophysics↗
USGS2010

Finite-difference modeling and dispersion analysis of high-frequency love waves for near-surface applications

Love-wave propagation has been a topic of interest to crustal, earthquake, and engineering seismologists for many years because it is independent of Poisson's ratio and more sensitive to shear (S)-wave velocity changes and layer thickness changes than are Rayleigh waves. It is well known that Love-wave generation requires the existence of a low S-wave velocity layer in a multilayered earth model. In order to study numerically the propagation of Love waves in a layered earth model and dispersion characteristics for near-surface applications, we simulate high-frequency (>5 Hz) Love waves by the staggered-grid finite-difference (FD) method. The air-earth boundary (the shear stress above the free surface) is treated using the stress-imaging technique. We use a two-layer model to demonstrate the accuracy of the staggered-grid modeling scheme. We also simulate four-layer models including a low-velocity layer (LVL) or a high-velocity layer (HVL) to analyze dispersive energy characteristics for near-surface applications. Results demonstrate that: (1) the staggered-grid FD code and stress-imaging technique are suitable for treating the free-surface boundary conditions for Love-wave modeling, (2) Love-wave inversion should be treated with extra care when a LVL exists because of a lack of LVL information in dispersions aggravating uncertainties in the inversion procedure, and (3) energy of high modes in a low-frequency range is very weak, so that it is difficult to estimate the cutoff frequency accurately, and "mode-crossing" occurs between the second higher and third higher modes when a HVL exists. ?? 2010 Birkh??user / Springer Basel AG.

Pure and Applied Geophysics↗
USGS2010

Seismic hazard assessment: Issues and alternatives

Seismic hazard and risk are two very important concepts in engineering design and other policy considerations. Although seismic hazard and risk have often been used inter-changeably, they are fundamentally different. Furthermore, seismic risk is more important in engineering design and other policy considerations. Seismic hazard assessment is an effort by earth scientists to quantify seismic hazard and its associated uncertainty in time and space and to provide seismic hazard estimates for seismic risk assessment and other applications. Although seismic hazard assessment is more a scientific issue, it deserves special attention because of its significant implication to society. Two approaches, probabilistic seismic hazard analysis (PSHA) and deterministic seismic hazard analysis (DSHA), are commonly used for seismic hazard assessment. Although PSHA has been pro-claimed as the best approach for seismic hazard assessment, it is scientifically flawed (i.e., the physics and mathematics that PSHA is based on are not valid). Use of PSHA could lead to either unsafe or overly conservative engineering design or public policy, each of which has dire consequences to society. On the other hand, DSHA is a viable approach for seismic hazard assessment even though it has been labeled as unreliable. The biggest drawback of DSHA is that the temporal characteristics (i.e., earthquake frequency of occurrence and the associated uncertainty) are often neglected. An alternative, seismic hazard analysis (SHA), utilizes earthquake science and statistics directly and provides a seismic hazard estimate that can be readily used for seismic risk assessment and other applications. ?? 2010 Springer Basel AG.

Pure and Applied Geophysics↗
USGS2010

CyberShake: A Physics-Based Seismic Hazard Model for Southern California

CyberShake, as part of the Southern California Earthquake Center’s (SCEC) Community Modeling Environment, is developing a methodology that explicitly incorporates deterministic source and wave propagation effects within seismic hazard calculations through the use of physics-based 3D ground motion simulations. To calculate a waveform-based seismic hazard estimate for a site of interest, we begin with Uniform California Earthquake Rupture Forecast, Version 2.0 (UCERF2.0) and identify all ruptures within 200 km of the site of interest. We convert the UCERF2.0 rupture definition into multiple rupture variations with differing hypocenter locations and slip distributions, resulting in about 415,000 rupture variations per site. Strain Green Tensors are calculated for the site of interest using the SCEC Community Velocity Model, Version 4 (CVM4), and then, using reciprocity, we calculate synthetic seismograms for each rupture variation. Peak intensity measures are then extracted from these synthetics and combined with the original rupture probabilities to produce probabilistic seismic hazard curves for the site. Being explicitly site-based, CyberShake directly samples the ground motion variability at that site over many earthquake cycles (i.e., rupture scenarios) and alleviates the need for the ergodic assumption that is implicitly included in traditional empirically based calculations. Thus far, we have simulated ruptures at over 200 sites in the Los Angeles region for ground shaking periods of 2 s and longer, providing the basis for the first generation CyberShake hazard maps. Our results indicate that the combination of rupture directivity and basin response effects can lead to an increase in the hazard level for some sites, relative to that given by a conventional Ground Motion Prediction Equation (GMPE). Additionally, and perhaps more importantly, we find that the physics-based hazard results are much more sensitive to the assumed magnitude-area relations and magnitude uncertainty estimates used in the definition of the ruptures than is found in the traditional GMPE approach. This reinforces the need for continued development of a better understanding of earthquake source characterization and the constitutive relations that govern the earthquake rupture process.

California↗
USGS2010

Aftershock distribution as a constraint on the geodetic model of coseismic slip for the 2004 Parkfield earthquake

Several studies of the 2004 Parkfield earthquake have linked the spatial distribution of the event’s aftershocks to the mainshock slip distribution on the fault. Using geodetic data, we find a model of coseismic slip for the 2004 Parkfield earthquake with the constraint that the edges of coseismic slip patches align with aftershocks. The constraint is applied by encouraging the curvature of coseismic slip in each model cell to be equal to the negative of the curvature of seismicity density. The large patch of peak slip about 15 km northwest of the 2004 hypocenter found in the curvature-constrained model is in good agreement in location and amplitude with previous geodetic studies and the majority of strong motion studies. The curvature-constrained solution shows slip primarily between aftershock “streaks” with the continuation of moderate levels of slip to the southeast. These observations are in good agreement with strong motion studies, but inconsistent with the majority of published geodetic slip models. Southeast of the 2004 hypocenter, a patch of peak slip observed in strong motion studies is absent from our curvature-constrained model, but the available GPS data do not resolve slip in this region. We conclude that the geodetic slip model constrained by the aftershock distribution fits the geodetic data quite well and that inconsistencies between models derived from seismic and geodetic data can be attributed largely to resolution issues.

California↗
USGS2010

Field survey of the March 28, 2005 Nias-Simeulue earthquake and Tsunami

On the evening of March 28, 2005 at 11:09 p.m. local time (16:09 UTC), a large earthquake occurred offshore of West Sumatra, Indonesia. With a moment magnitude (Mw) of 8.6, the event caused substantial shaking damage and land level changes between Simeulue Island in the north and the Batu Islands in the south. The earthquake also generated a tsunami, which was observed throughout the source region as well as on distant tide gauges. While the tsunami was not as extreme as the tsunami of December 26th, 2004, it did cause significant flooding and damage at some locations. The spatial and temporal proximity of the two events led to a unique set of observational data from the earthquake and tsunami as well as insights relevant to tsunami hazard planning and education efforts. ?? 2010 Springer Basel AG.

Pure and Applied Geophysics↗
USGS2010

Combined effects of tectonic and landslide-generated Tsunami Runup at Seward, Alaska during the Mw 9.2 1964 earthquake

We apply a recently developed and validated numerical model of tsunami propagation and runup to study the inundation of Resurrection Bay and the town of Seward by the 1964 Alaska tsunami. Seward was hit by both tectonic and landslide-generated tsunami waves during the M W "> M W MW 9.2 1964 megathrust earthquake. The earthquake triggered a series of submarine mass failures around the fjord, which resulted in landsliding of part of the coastline into the water, along with the loss of the port facilities. These submarine mass failures generated local waves in the bay within 5 min of the beginning of strong ground motion. Recent studies estimate the total volume of underwater slide material that moved in Resurrection Bay to be about 211 million m 3 (Haeussler et al. in Submarine mass movements and their consequences, pp 269–278, 2007 ). The first tectonic tsunami wave arrived in Resurrection Bay about 30 min after the main shock and was about the same height as the local landslide-generated waves. Our previous numerical study, which focused only on the local landslide-generated waves in Resurrection Bay, demonstrated that they were produced by a number of different slope failures, and estimated relative contributions of different submarine slide complexes into tsunami amplitudes (Suleimani et al. in Pure Appl Geophys 166:131–152, 2009 ). This work extends the previous study by calculating tsunami inundation in Resurrection Bay caused by the combined impact of landslide-generated waves and the tectonic tsunami, and comparing the composite inundation area with observations. To simulate landslide tsunami runup in Seward, we use a viscous slide model of Jiang and LeBlond (J Phys Oceanogr 24(3):559–572, 1994 ) coupled with nonlinear shallow water equations. The input data set includes a high resolution multibeam bathymetry and LIDAR topography grid of Resurrection Bay, and an initial thickness of slide material based on pre- and post-earthquake bathymetry difference maps. For simulation of tectonic tsunami runup, we derive the 1964 coseismic deformations from detailed slip distribution in the rupture area, and use them as an initial condition for propagation of the tectonic tsunami. The numerical model employs nonlinear shallow water equations formulated for depth-averaged water fluxes, and calculates a temporal position of the shoreline using a free-surface moving boundary algorithm. We find that the calculated tsunami runup in Seward caused first by local submarine landslide-generated waves, and later by a tectonic tsunami, is in good agreement with observations of the inundation zone. The analysis of inundation caused by two different tsunami sources improves our understanding of their relative contributions, and supports tsunami risk mitigation in south-central Alaska. The record of the 1964 earthquake, tsunami, and submarine landslides, combined with the high-resolution topography and bathymetry of Resurrection Bay make it an ideal location for studying tectonic tsunamis in coastal regions susceptible to underwater landslides.

Alaska↗
USGS2011

Book review: Extreme ocean waves

‘‘Extreme Ocean Waves’’ is a collection of ten papers edited by Efim Pelinovsky and Christian Kharif that followed the April 2007 meeting of the General Assembly of the European Geosciences Union. A note on terminology: extreme waves in this volume broadly encompass different types of waves, includ- ing deep-water and shallow-water rogue waves (alternatively termed freak waves), storm surges from cyclones, and internal waves. Other types of waves such as tsunamis or rissaga (meteotsunamis) are not discussed in this volume. It is generally implied that ‘‘extreme’’ has a statistical connotation relative to the average or significant wave height specific to each type of wave. Throughout the book, in fact, the reader will find a combination of theoretical and statistical/ empirical treatment necessary for the complete examination of this subject. In the introduction, the editors underscore the importance of studying extreme waves, documenting several dramatic instances of damaging extreme waves that occurred in 2007. Review info: Extreme Ocean Waves . By E. Pelinovsky and C. Kharif (eds), 2008. ISBN: 978-1402083136, xiii, 196 pp.

Pure and Applied Geophysics↗
USGS2011

Book review: Nonlinear ocean waves and the inverse scattering transform

Nonlinear Ocean Waves and the Inverse Scattering Transform is a comprehensive examination of ocean waves built upon the theory of nonlinear Fourier analysis. The renowned author, Alfred R. Osborne, is perhaps best known for the discovery of internal solitons in the Andaman Sea during the 1970s. In this book, he provides an extensive treatment of nonlinear water waves based on a nonlinear spectral theory known as the inverse scattering transform. The writing is exceptional throughout the book, which is particularly useful in explaining some of the more difficult mathematical concepts. Review info: N onlinear Ocean Waves and the Inverse Scattering Transform. B y Alfred R. Osborne, 2010. ISBN: 978-125286299, 917 pp.

Pure and Applied Geophysics↗
USGS2011

The role of adsorbed water on the friction of a layer of submicron particles

Anomalously low values of friction observed in layers of submicron particles deformed in simple shear at high slip velocities are explained as the consequence of a one nanometer thick layer of water adsorbed on the particles. The observed transition from normal friction with an apparent coefficient near μ = 0.6 at low slip speeds to a coefficient near μ = 0.3 at higher slip speeds is attributed to competition between the time required to extrude the water layer from between neighboring particles in a force chain and the average lifetime of the chain. At low slip speeds the time required for extrusion is less than the average lifetime of a chain so the particles make contact and lock. As slip speed increases, the average lifetime of a chain decreases until it is less than the extrusion time and the particles in a force chain never come into direct contact. If the adsorbed water layer enables the otherwise rough particles to rotate, the coefficient of friction will drop to μ = 0.3, appropriate for rotating spheres. At the highest slip speeds particle temperatures rise above 100°C, the water layer vaporizes, the particles contact and lock, and the coefficient of friction rises to μ = 0.6. The observed onset of weakening at slip speeds near 0.001 m/s is consistent with the measured viscosity of a 1 nm thick layer of adsorbed water, with a minimum particle radius of approximately 20 nm, and with reasonable assumptions about the distribution of force chains guided by experimental observation. The reduction of friction and the range of velocities over which it occurs decrease with increasing normal stress, as predicted by the model. Moreover, the analysis predicts that this high-speed weakening mechanism should operate only for particles with radii smaller than approximately 1 μm. For larger particles the slip speed required for weakening is so large that frictional heating will evaporate the adsorbed water and weakening will not occur.

Pure and Applied Geophysics↗
USGS2011

Earthquake rupture at focal depth, part II: mechanics of the 2004 M2.2 earthquake along the Pretorius Fault, TauTona Mine, South Africa

We analyze here the rupture mechanics of the 2004, M2.2 earthquake based on our observations and measurements at focal depth (Part I). This event ruptured the Archean Pretorius fault that has been inactive for at least 2 Ga, and was reactivated due to mining operations down to a depth of 3.6 km depth. Thus, it was expected that the Pretorius fault zone will fail similarly to an intact rock body independently of its ancient healed structure. Our analysis reveals a few puzzling features of the M2.2 rupture-zone: (1) the earthquake ruptured four, non-parallel, cataclasite bearing segments of the ancient Pretorius fault-zone; (2) slip occurred almost exclusively along the cataclasite-host rock contacts of the slipping segments; (3) the local in-situ stress field is not favorable to slip along any of these four segments; and (4) the Archean cataclasite is pervasively sintered and cemented to become brittle and strong. To resolve these observations, we conducted rock mechanics experiments on the fault-rocks and host-rocks and found a strong mechanical contrast between the quartzitic cataclasite zones, with elastic-brittle rheology, and the host quartzites, with damage, elastic–plastic rheology. The finite-element modeling of a heterogeneous fault-zone with the measured mechanical contrast indicates that the slip is likely to reactivate the ancient cataclasite-bearing segments, as observed, due to the strong mechanical contrast between the cataclasite and the host quartzitic rock.

Tautona Mine↗