Numerous studies have illustrated that mineral transformations have the capability to induce faulting at elevated pressure and temperature (PT), circumstances in which ductile flow would typically dominate. This mechanism, commonly known as transformational faulting, emerges as a plausible explanation for the puzzling phenomenon of deep-focus earthquakes occurring at depths up to 700 km. Currently, the debate partly revolves around determining why certain phase transformations lead to faulting while others do not. To better understand this phenomenon, we can compare different transformations taking place in similar experimental conditions and see how they do or do not cause strain localization and faulting. In this regard, we conducted a series of five deformation experiments in the large volume press at the PB61 beamline at DESY synchrotron. Two of these experiments involved deforming germanium-olivine samples as they transformed into ringwoodite (the high-pressure phase). The other three experiments were carried out on quartzite (novaculite) samples while they were transforming to coesite. Throughout the experiments, we collected X-ray diffraction patterns and images concurrently with the collection of Acoustic Emissions (AEs).

The results indicate, in both quartz and olivine experiments, the growth of the high-pressure phase at various rates depending on PT conditions and equilibrium overstep. Specifically, we observed rapid olivine-ringwoodite kinetics at elevated PT, far from equilibrium, while slower kinetics were noted for the quartz-coesite transformation. Thousands of AEs were collected in each experiment, and their locations reconstructed using arrival times on the six transducers used. Interestingly, the spatial distribution of these AEs revealed that for some quartz-coesite experiments, AEs originated from fault planes that formed within the initially intact rock cores. Furthermore, an analysis of the AE catalogues, focusing on the magnitude-frequency distribution, revealed a wide range of b-values influenced by varying PT conditions and transformation kinetics. This observation underscores the different underlying mechanisms since the obtained b-values are high when transformation and strain are distributed and lower when strain is localized (i.e., when a fault plane develops).

Our study supports the major role of mineral transformations in inducing faulting under high PT. These findings will help better quantify the intricate relationships between mineral transformations and faulting and in turn contribute to a better understanding of the fundamental geological processes behind deep and intermediate earthquakes.

How to cite: Mingardi, G., Gasc, J., Farla, R., and Schubnel, A.: Combining synchrotron and acoustic emission techniques to reveal the secrets of high PT faulting, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-21427, https://doi.org/10.5194/egusphere-egu24-21427, 2024.

Seismicity studies on both fast- and slow-spreading ridge systems have found along-strike variations in mantle mechanical behavior on oceanic transform faults (OTFs), at pressure and temperature conditions above the long-term brittle-ductile transition of peridotites at ~700°C. Plate motion on some sections of the fault is accommodated by aseismic slip only (ductile deformation), whereas motion on other sections is by slip and deep swarms of microearthquakes (semi-brittle deformation) of mantle rocks (e.g., McGuire et al., 2012; Yu et al., 2021). To explore the mechanisms responsible for lateral variations in mantle mechanical behavior and the occurrence of this deep mantle microseismicity, we carried out an integrated study on peridotite mylonites dredged from two OTFs on the Southwest Indian Ridge that record deformation at 700-1000°C.

The samples show variable degrees of deformation, ranging from proto- to ultra-mylonitic textures. The most deformed zones of the mylonites are characterized by an increase in the proportion of fine grained (<10 micron) mylonitic shear bands compared to coarse grained (millimeter) porphyroclasts inherited from the protolith. These shear bands contain syn-deformation chlorine-rich amphibole indicating seawater-peridotite interaction during shear band formation.

Olivine and pyroxene porphyroclasts in protomylonites contain evidence for intense brittle deformation. The presence of subgrain walls, high aspect ratios, and internal misorientations crosscut by fractures imply that they deformed by low-temperature plasticity before brittle deformation. Fractures are sealed by the fine-grained shear bands present in the samples. In (ultra)mylonites, porphyroclasts also show evidence of fracturing after flowing through low-T plasticity. Fracturing was coeval with viscous flow of surrounding weak and hydrated mylonitic shear bands and triggered by hardening of larger grains due to dislocation accumulation. From existing flow laws, such brittle deformation of peridotite minerals necessitates high strain rate deformation, from 10^-9 to 10^-5s^-1, similar to strain rates associated with slow slip events.

From these results we propose that swarms of microseismicity on OTFs are triggered by deformation of a heterogeneous mantle. Seismic rupture occurs in lenses of coarse-grained peridotites, possibly driven by aseismic creep of surrounding hydrated mylonitic shear zones. Importantly, observations also suggest plastic flow of brittle (seismic) patches before rupture.

How to cite: Cécile, P. and Jessica, W.: Origin of brittle deformation and microseismicity in the ‘ductile’ mantle on oceanic transform faults, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-16752, https://doi.org/10.5194/egusphere-egu24-16752, 2024.

Dehydration reactions can influence the occurrence of earthquakes at a range of depths, highlighting the importance of understanding these reactions in the study of seismic activity. We have performed several pilot experiments that included the construction of a controllable fast pressure drop unit attached to the piston-cylinder apparatus. This setup makes it possible to simulate conditions that represent a fast pressure drop during an earthquake event. We focused on serpentinite dehydration because 1/ it plays an important link between the deep geodynamic processes occurring in subduction zones and the seismic and volcanic activity and 2/ the interplay between serpentinite dehydration and deformation during the earthquake cycle is not yet fully understood. To test the experimental setting, we first performed a series of static experiments under the conditions that are already in the olivine stability field. After the static experiments at high pressure (1.1 GPa), we performed the controlled fast pressure drop experiments to 0.3 GPa as well as the ramping experiments, in which a series of pressure build-ups and drops were performed maintaining the high temperatures (570 to 640 °C) to simulate the earthquake cycle. In these experiments, olivine was an order of magnitude more abundant than in the one-hour low-pressure static experiment. The pressure drop occurs in seconds. The ramping experiment lasted only 10 mins before cooling down. The results may challenge conventional wisdom about the timescales of mineral reactions under extreme conditions, such as during earthquakes.

How to cite: Tajčmanová, L., Cionoiu, S., and Schuppenhauer, D.: Could earthquakes cause rapid dehydration of serpentinite?, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-19186, https://doi.org/10.5194/egusphere-egu24-19186, 2024.

Understanding the mechanisms controlling brittle rock failure at the grain to sub-grain scale is a fundamental challenge in geosciences. Recent advances in triaxial compression and dynamic shock experiments combined with dynamic X-ray microtomography provide unparalleled insights into the 3D strain field evolution within deforming rocks. However, these methods do not accurately predict the heterogeneous internal stress field prior to failure, which is crucial for predicting microfracture initiation and propagation, leading to macroscopic failure. In the past decade, efforts have focused on developing synchrotron X-ray diffraction techniques leveraging the high penetrative capacity of hard X-rays from the last generations of synchrotron light sources. These techniques offer spatially resolved information on crystal phase orientation and elastic strain within a 3D volume. The local orientation and elastic strain tensor is reconstructed grain-by-grain, with precision down to approximately 10^-3 radian for orientation and 10^-4 for strain. Stress is then calculated using Hooke's law for anisotropic materials and the elastic constants of the crystal phases. We employed 3D X-ray diffraction to investigate the internal stress field evolution in a rock core sample deformed under triaxial compression in the Hades apparatus. A 5mm-diameter core of Berea sandstone was subjected to axial step loading under constant radial stress of 10 MPa, reaching brittle failure at around 90 MPa differential stress. Elastic strain of individual quartz grains were measured at different load steps, and elastic stresses were calculated, providing maps of the internal strain and stress field in the sample. Results reveal progressive elastic shortening of quartz grains parallel to the compression axis and elongation in orthogonal directions due to the Poisson’s effect. Reorientation of principal stress components is also observed with increasing axial stress, which tend to align with the macroscopic stress field. Internal stresses distribution varies within a range of ca. 300 MPa, suggesting local stress amplifications occurred interpreted as force chains, potentially favoring crack nucleation. This experiment is among the first ones to characterize in-situ the stress distribution in a natural rock under compressive loading, and demonstrates the potential of synchrotron diffraction techniques for investigating strain and stress in geological materials.

How to cite: Jacob, J., Cordonnier, B., Wright, J., and Renard, F.: Mapping internal stress field in deforming rocks using synchrotron high-energy X-ray diffraction, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-3863, https://doi.org/10.5194/egusphere-egu24-3863, 2024.

Strength-depth profiles for ductile portions of continental crust are derived from either extrapolation of flow laws from deformation experiments or paleopiezometric estimates in deformed and nominally hydrated plate margins. Lower continental crust in intracontinental settings, in contrast, is relatively dry and should be considerably stronger than the lower crust of hydrated plate margins. The relative strengths of dry quartz and feldspar are poorly constrained by experiments and paleopiezometric estimates from such rocks are sparse. As such, the strength of intracratonic lower crust is difficult to ascertain. Here, we use a recently calibrated subgrain-size piezometer to estimate paleostresses from feldspar and quartz deformed in relatively dry (<20 ppm H2O) lower continental crust of the Musgrave Ranges in central Australia. Neocrysts of plagioclase, K-feldspar, and quartz mantle partially recrystallized porphyroclasts, which is indicative of bulging and subgrain-rotation recrystallization. Using crystallographic preferred orientations and plotting misorientation axes of subgrain boundaries of each phase, we infer that dislocation creep involved the slip systems (010)[100] and (010)[001] for plagioclase, (010)[101] for K-feldspar, and (0001)<11-20> and {01-10}<0001> for quartz. Titanium in quartz and gradients in concentration of Ca and K in feldspars within neocrysts and along subgrain boundaries verify that subgrains in all three phases were formed at a temperature of ~650°C under dry, eclogite-facies conditions. Subgrain sizes of 10.6–18.1 µm in quartz, 11.5–16.9 µm in plagioclase, and 12.0–17.5 µm in K-feldspar correspond to differential paleostresses between 22–36 MPa and are consistent with a single mean paleostress of 28 MPa. Our results demonstrate that there is minimal stress partitioning between dry quartz, plagioclase and K-feldspar under typical crustal thermal gradients. Moreover, the differential stress accommodated by felsic rocks in the Davenport shear zone is lower than predicted by previous strength-depth profiles of lower cratonic crust.

How to cite: Osinchuk, A., Dyck, B., Wallis, D., and Camacho, A.: Subgrain-size piezometry of feldspar and quartz records a single paleostress from dry lower continental crust, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-16849, https://doi.org/10.5194/egusphere-egu24-16849, 2024.

Field studies established that seismicity in the lower crust is linked to brittle failure of dry, strong rocks. Failure of these strong rocks implies build-up of differential stresses to gigapascal (GPa) levels, but this requirement contrasts with current models of continental lithosphere deformation, which favour distributed flow of weak viscous lower crust. Although several mechanisms have been proposed to generate transiently high stresses, direct measurements are lacking. Recent advancements in microanalytical techniques (i.e., high-angular resolution electron backscatter diffraction, HR-EBSD) have proven successful at measuring the residual stress resulting from elastic strain retained in mineral grains.

We investigated with HR-EBSD the residual stresses in garnet and diopside from exhumed faults containing pseudotachylytes (quenched frictional melts produced during seismic slip). The samples come from Lofoten (Norway), Holsnøy (Bergen Arcs, Norway), and Musgrave Ranges (Central Australia). Pseudotachylytes from all three localities represent single earthquake events and formed at lower-crustal conditions (T = 500–720 °C, P = 0.5–1.0 GPa). Pseudotachylytes from Holsnøy show an asymmetric damage distribution, where host-rock garnet is pulverized nearby the fault on the side subjected to predominantly tensional stresses during rupture propagation, while garnet is intact on the other side. This asymmetry provides an opportunity to compare the residual stresses on both sides of the fault.

All samples preserve intragrain residual stress heterogeneities reaching 100s of MPa to GPa levels due to local high density of unrecovered lattice defects (dislocations). However, the timing of formation of lattice defects with respect to the seismic event differs. In samples from the Musgrave Ranges, the absence of any later deformation along with the sluggish mobility of dislocations in garnet at the ambient deformation conditions (500 °C, 0.5 GPa) allowed preservation of the high dislocation density produced during the earthquake rupture propagation, recording stress heterogeneities of as much as 6 GPa. In Holsnøy, residual stress heterogeneities of up to 1 GPa are only measured in pulverized grains and are also associated with unrelaxed dislocations generated during the earthquake rupture propagation. Intact garnet grains from the less damaged side of the fault show a limited range of intragrain stress heterogeneities, generally within 100 MPa, and a low density of dislocations. Residual stresses in diopside from Lofoten are only elevated (600 MPa) within 200 µm of the pseudotachylyte. Diopside recorded the progressive build-up of stress during interseismic loading, as suggested by the presence of coseismic cracks crosscutting lattice undulations that preserve the greatest stress heterogeneities. However, the ability of diopside to build up stress is limited, as stress is efficiently dissipated by the development of deformation twins.

In conclusion, great stress heterogeneities can be preserved in mineral grains that experienced the earthquake cycle in the lower crust. Different mineral phases can preserve stress heterogeneities to different extents, depending on the mobility of dislocations after their formation and on other relaxation mechanisms (e.g., twinning). Information on residual stress have important implications for the energy budget of an earthquake, the earthquake cycle deformation, and crustal rheology.

How to cite: Menegon, L., Toffol, G., van Schrojenstein Lantman, H. W., Wallis, D., Pennacchioni, G., and Jamtveit, B.: Earthquake induced residual stresses preserved in fault rocks exhumed from the lower crust, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-9228, https://doi.org/10.5194/egusphere-egu24-9228, 2024.

Understanding the stress conditions of active subduction zones has been a longstanding hurdle with critical implications for natural disasters considering stress/strain orientations and magnitudes can control shallow earthquakes and tsunamigenesis. The Sestola-Vidiciatico Unit (SVU) in the Northern Apennines is an exhumed subduction channel with exposures of up to 9 km paleodepth, having reached up to 200°C. This unit experienced a relatively limited deformation history and serves as a rare analog to the shallowest portions of active subduction megathrusts. We use calcite twin data from shear veins along mineralized faults surrounding the exhumed subduction interface to reconstruct paleostress orientations through calcite twin stress inversion. Combining orientation data with calcite twin paleopiezometry and geothermometry, we are able to reconstruct the stress state of the SVU during peak subduction and subsequent exhumation.

During subduction, the maximum principal stress axis was oriented at a low angle to the subduction interface and the minimum principal stress axis oriented at a high angle, indicating N/NE directed compression. As subduction ceased and exhumation initiated, stress orientations inverted with the maximum principal stress axis becoming oriented at a high angle to the subduction interface and the minimum principal stress axis oriented at a low angle, indicating N/NE directed extension driven by primarily the weight of overburden material. These findings are consistent with theoretical orientations for both of these tectonic regimes and agree with previous studies interpreting subduction zone stress orientations. Calcite twin paleopiezometry and geothermometry suggests the rotation of principal stresses coincides with higher differential stresses during early exhumation. Based on the interpreted differential stresses and the reconstructed paleostress orientations, we model different possible explanations including contrasting mechanical strength between the contractional and extensional faults or changes in pore fluid pressure conditions between the two different tectonic regimes.

How to cite: Williams, S., French, M., and Rubin, C.: Determining the deformation temperatures and paleostress conditions of the Sestola-Vidiciatico Unit in the Northern Apennines, an exhumed shallow subduction zone, using calcite deformation twins, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-14306, https://doi.org/10.5194/egusphere-egu24-14306, 2024.

Dykes are essentially magma filled fractures within the earth’s crust often formed by the pressure imparted by the intruding magma. Magnitude of the magma overpressure has been traditionally determined utilizing elastic properties of the host rock and the complete dimension i.e., full length and maximum width of the fractures. Full exposures of dykes (from tip to tip) are rare, however, as most of the dyke bodies encountered in the field are subject to erosion or disruption along its length as a result of geological time, making estimation of aspect ratios challenging.

We propose a new method of estimating total length and maximum width of dykes from their partial outcrops featuring at least one exposed tip. Taking into account the fact that majority of dykes form as dominantly opening mode fractures with an elliptical shape of opening, the method involves solving the equation of this ellipse using every conceivable combination of a pair of ground points recorded on the dyke margin considering the visible tip as the origin. Validity of the method has been checked using published data obtained from incomplete dyke outcrops exposed in the caldera walls of Miyake-jima volcano in Japan. The calculated estimates are in line with the results acquired through a previous published method. The present method has been effectively utilized to calculate the aspect ratios of partially exposed mafic dykes emplaced within the younger granite of the Chitradurga Schist Belt in the western Dharwar craton of peninsular India. We discuss the ranges of their magma overpressure and depths of origin as well as the stress intensity factors associated with the host granite.

How to cite: Mondal, T. K. and Biswas, S. K.: Dykes and their magma overpressure, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-392, https://doi.org/10.5194/egusphere-egu24-392, 2024.

Faults are normally thought to present shear fractures that develop at an angle to the main principal stresses, so that they have shear stresses active parallel to the fault plane and thus move. Here we present two “fault” features that deviate from this principle, they develop not due to stress but during kinematic movement, they are both oriented parallel to two of the main principle stresses and as such have no shear stresses in their planes. On the large tectonic plate scale one of these features are the well known transform faults between mid ocean ridges. The ridges themselves are extensional features with the lowest principle stress perpendicular to the ridge. Transform faults are oriented perpendicular to the ridges and show movement only or mainly between the ridges where the plates move in opposite directions. These are faults that do not develop due to shear stress, they develop only because of differential movement and are therefore only or mainly kinematic. On the small scale a very similar feature is the side of a stylolite tooth. Stylolites are dissolution features, they are thus in a way the opposite to mid ocean ridges and have the largest principal stress oriented perpendicular to the stylolite plane. Due to differential movement and growth of the stylolite roughness they develop steep teeth where the sides of teeth become oriented perpendicular to the stylolite plane. These are also movement surfaces or “faults” that have no shear stress.

What does this mean for stress inversion analysis? At least stylolite teeth show a quite pronounced set of striations on their sides and slikolites are also often developed on fault planes, at least in limestone. How do we separate a purely kinematic from a stress-related fault? This discussion and the potential consequences for stress inversion studies is not new, but remains to be very important and should be debated.

How to cite: Koehn, D., Hafermaas, D., and Koehler, S.: The development of kinematic shear-stress free faults , EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-10900, https://doi.org/10.5194/egusphere-egu24-10900, 2024.

Knowledge of the present‐day stress field in the Earth's crust is key for understanding the mechanical behaviour of rocks and structures under tectonic forces. The study of stress fields therefore remains a pivotal area for understanding the mechanical behaviour of rocks, fluid flow at depth and in revealing mechanisms that cause tectonic plates to creep, fail, or rupture. stress patterns in the Earth's crust appear on different scales: first order (plate scale), second order (regional scale), and third order (local scale). The latter is mainly controlled by basin geometry, topography, local inclusions, density contrasts, and active faults and can mask regional and plate stress patterns.

In this contribution, a couple of examples of stress states at the local, regional and large scale are presented using borehole breakouts as main stress indicator for the current stress field orientation.

In order to understand the influence of stress field evolution at local scale, a case study in Hawai´i concerns the effects of the two large overlapping shield volcanoes on the stress field at depth. The analysis reveals that the two-competing gravitational loads primary control the orientation of the present-day stress field, which deviates significantly from the plate and regional tectonic stress field. Therefore, knowledge of local and shallow stress fields can have a significant impact on future borehole planning. From a more regional point of view, an example of current stress orientation in Sweden is presented. The main objectives are to constrain the orientation of horizontal stresses using borehole data, and to discuss implications for geothermal exploration. Thus, obtaining detailed and accurate data on the stress state is of paramount significance to optimise the design of underground installations in order to maximise fluid flow and minimise the risks of wellbore instability. Finally, a large-scale study in Italy investigates the stress field at plate scale to reveal whether the orientation of horizontal stresses may change with depth or laterally indicating stress perturbations and heterogeneities related to areas with complex geo-tectonic setting.

In conclusion, this contribution aims to illustrate and emphasise the relevance of determining the horizontal stress orientation at depth in order to improve the understanding of subsurface stress fields and their applications in different fields of geosciences and different geological settings.

How to cite: Pierdominici, S.: Reconstruction of the state of stress in the upper crust by borehole breakouts stress indicators, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-5744, https://doi.org/10.5194/egusphere-egu24-5744, 2024.

Transient stress perturbations caused by passing waves of distant earthquakes have been observed to activate fault slip. Observations of remotely triggered earthquakes at distances greater than ~2-3 mainshock fault lengths suggest that certain conditions promote fault activation, including large-amplitude shaking at periods below ~ 10 seconds within a geothermal setting and extensional and/or transtensional tectonics. Yet, it is still unclear if remote dynamic triggering is ubiquitous. An additional complication to determining the prevalence of triggering is that there are likely many small-magnitude earthquakes within sparsely instrumented regions that are uncataloged. Additionally, large mainshock signals can mask smaller local events that also are missing from the local catalogs. As a result, possible triggering mechanism(s) remain enigmatic. Bounding the necessary physical conditions for remote triggering, such as determining the upper or lower bounds of stress or strain amplitudes, the orientation of the seismic wave’s traversal with respect to the local stress field or fault geometry, or the geologic properties conducive to triggering can help provide clues about the physics of remote triggering. 

            The northern Chilean subduction margin provides an ideal setting to study remote dynamic triggering. Its dense instrumentation provides a long history of both seismic and aseismic deformation in both the subduction system and forearc faults, including the Atacama fault system. Our investigation combines a new, detailed regional earthquake catalog (2007-2021) from Sippl et al., (2023) and documented cases of triggered aseismic slip in the Atacama fault system (Victor et al., 2018). We use a twofold approach to determine the prevalence of earthquake triggering by candidate mainshocks that produce strains at our target location ranging from 1 to ~140 microstrain. The approach uses 1) a difference-of-means test of cataloged seismicity outside of the mainshock cluster (including foreshocks and aftershocks), and 2) a waveform-based approach to look for earthquake triggering at seismic stations located close to creepmeters that recorded triggered aseismic slip events.  We find a lack of evidence of persistent, statistically significant seismicity increases outside of the mainshock cluster associated with any of the candidate mainshocks.  Notably, there is an absence of significant seismicity changes outside of clustered foreshock or aftershock seismicity associated with the series of 11 M6.2-8.2 earthquakes that produced high-strain-rate events during the 2014 Iquique sequence. Seismic recordings of the 2011 M9.1 Japan earthquake at stations CX.PB01-CX.PB14 located near creep meter stations CAR3 and CH01 on the Mejillones peninsula near the Chomache fault reveal evidence of remote triggering. We observe local, uncataloged earthquakes that are only visible after applying a high pass filter that removes the mainshock signal that otherwise overprinted and swamped the local signals.  The uniformity of particle motions (circular or oblong) generated by local earthquakes on multiple stations (N=9) is lacking in most non-triggering mainshocks. This uniformity suggests that the orientation of transient stress perturbations imparted by the mainshock waveforms, in relation to the local fault orientations, may play a role in the triggering process. 

How to cite: Harrington, R. M., Kilb, D., Verdecchia, A., and Victor, P.: Putting faults into motion by remote dynamic triggering: local ground-motion orientations seem to eclipse strain in the northern Chilean subduction forearc crust, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-6650, https://doi.org/10.5194/egusphere-egu24-6650, 2024.

Short-term earthquake prediction remains one of the primary goals of seismotectonics. Here detailed observations of unusual fault kinematic behaviour and near-surface crustal stress variations are presented from before, during, and shortly after an earthquake series which culminated with two Mw 4.6 and 4.4 events near Breitenau, Vienna Basin, Austria, on 30 March 2021 and 19 April 2021, respectively. The oblique normal NNE-SSW trending Pitten Fault is exposed in Altaquelle Cave close to the southern margin of the Vienna Basin in the eastern Alps, which is known to have hosted several historical earthquakes of Mw = > 5. This cave has developed in Triassic marbles of the Central Alpine Permomesozoic. The observed branch of this active steeply dipping fault is associated with the seismogenic sinistral Vienna Basin Fault and the NE-SW trending Mur-Mürz Fault. To investigate the fault activity, TM71 moiré extensometers have been used to obtain precise three-dimensional records of fault kinematic behaviour at the micron scale while the recently developed SMB2018 protocol has been used to define the stress state associated with each fault reactivation event. The observations were then compared to the Copernicus European Ground Motion Service InSAR time series derived from Sentinel-1 data. From late 2018 to early 2021, the three-dimensional kinematic behaviour of the fault comprised a variety of different on-plane as well as out-of-plane hanging block displacements ranging in magnitude from 3 to 19 μm. Then, around the time of the earthquake series in 2021, four significant displacement events were recorded: (i) 0.186 mm along a vector of 186/-12° (i.e. upward) on 16 March; (ii) 0.615 mm along a vector of 177/-88° (upward) on 26 March; (iii) 0.066 mm along a vector of 013/26° (downward) on 30 March; and (iv) 0.022 mm along a vector of 308/54° (downward) on 11 May. The third of these events occurred on the same day as the largest earthquake. These events are all much larger than any other record of fault displacement recorded in the Eastern Alps since 2013. This contribution details this unusual fault displacement behaviour and compares the calculated stress states with both the focal solutions for each earthquake and InSAR maps of E-W and vertical ground motion. A comprehensive understanding of this important seismotectonic event helps to shed further light on potential earthquake precursory phenomena.

How to cite: Melichar, R., Baroň, I., Rowberry, M., Jelének, J., Sokol, Ľ., del Puy Papí Isaba, M., Freudenthaler, C., Hausmann, H., Plan, L., Grasemann, B., Stemberk, J., Schultz, R. A., and Bürgmann, R.: Unusual fault kinematic behaviour and near-surface crustal stress variations before and during an earthquake series in the Vienna Basin (Austria) in spring 2021, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-7664, https://doi.org/10.5194/egusphere-egu24-7664, 2024.

The 1D Geomechanical Models (1DGM) were done for four vertical boreholes in the Early Paleozoic shale sequences of the Baltic Basin in Poland. The models assumed an elastic rock behaviour with anisotropy in VTI symmetry. The far-field horizontal stresses were calculated as the sum of two components: the vertical stress derivative acting on the horizontally constrained rock column and the effect of elastic tectonic strains. Local stresses in the borehole wall were assumed to induce breakouts (BBs) and drilling-induced tensile fractures (DITFs). Hydraulic fracturing tests additionally validated the stress modelling results.

Micro-resistivity images from XRMI logging revealed irregular BBs, usually confined to individual layers, often non-symmetrical, and with a tendency to encompass the entire borehole wall. Despite their irregularity, statistical analysis of their orientation provides a good quality and stable stress orientation.

The initial modelling results, balancing the cumulative length of the modelled (BB_M) and observed (BB_O), revealed a systematic misfit between BBm and BBo locations. A detailed comparison between the BBm and BBo intervals concluded that the artificial degradation of Young's modulus and Poisson's ratio is caused by the perturbation of the velocity of the acoustic wave from the dipole acoustic tool passing through the intervals with irregular BBs. This makes it impossible to model the BBs in the places where they are present. To deal with this, the initial stress models were recalculated to the final models in which the BBm were avoided in intervals where there were no BB_O.

The initial and final stress models differ significantly in terms of the tectonic strain values and the combined length of the BBm. Still, their stress profiles are similar due to the small contribution of tectonic strain to the far-field stresses. We concluded that irregular BBs developed due to small differential horizontal stresses, causing abrupt BB failure with rapidly growing angular width. The stress layering between lithostratigraphic units was obtained with a dominance of the normal faulting stress regime in the lower borehole sections and the reverse faulting present in the upper sections. The minor regional elastic tectonic strain value for the shale sequence was determined to be an order of magnitude lower than the strain in the crystalline basement, as determined from the satellite geodetic strain rate. We expect that this discrepancy could be explained by a higher rate of viscous relaxation in the shale sequence with > 60% of the clay mineral content. This suggests the need to implement the viscous relaxation into the 1DGM of sedimentary sequences.

How to cite: Jarosinski, M., Bobek, K., and Pachytel, R.: Geomechanical models of the shale sequence of the Baltic Basin (Poland): possible case of elastic properties degradation and viscous stress relaxation, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-20951, https://doi.org/10.5194/egusphere-egu24-20951, 2024.

Understanding the intricate dynamics of mantle flow and their influence on lithospheric stress patterns is critical for assessing reservoir responses to potential CO2 or nuclear waste storage, as well as for hazard and risk assessment. Stress patterns play a first-order control in the mechanical response within inherited tectonic structures, with varying stress sources governing different spatiotemporal scales. Our understanding of the present-day mantle flow state has much improved over the past decades, reflected in models that are consistent with first-order features. The World Stress Map (WSM) project serves as a primary observational dataset to validate our understanding of Earth's dynamics through a global compilation of crustal stress indicators.

Here we study mantle flow models as a simplified superposition of Couette and Poiseuille flow types, which have been useful in explaining sub-continental scale deformation in the lithosphere. We aim to understand the role of mantle flow as a stress driver by generating stress fields from an analytical representation of upper mantle flow, derived from the superposition of steady-state flow models. Our approach allows us to build first-order expectations and conduct fast hypothesis testing for upper mantle flow states.

How to cite: Hayek Valencia, J. N., Stotz, I. L., Bunge, H.-P., Carena, S., and Ghelichkhan, S.: First-order global stress patterns inferred from hierarchies of upper mantle flow models, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-16807, https://doi.org/10.5194/egusphere-egu24-16807, 2024.

To comprehend the rheology of the Earth's crust and the relevant rock properties, one key approach is to deform rocks and minerals at elevated pressures and temperatures and then extrapolate the measured stress and strain rate values to natural conditions using constitutive equations. Laboratory experiments are mostly conducted on monomineralic rocks, with quartz being considered as the weakest constituent of the middle continental crust. However, field observations suggest that this is an oversimplification, and polymineralic fault rocks may be weaker than monomineralic quartz rocks. This study presents the first experiments on fine-grained, solid, natural rock samples, containing their natural homogeneities and inhomogeneities, demonstrating that granitoid rocks may be weaker than quartz at mid-crustal conditions. It also highlights the importance of pre-existing faults and polymineralic fine-grained zones for strain localisation and proposes values for extrapolation to natural conditions and their use in numerical models of the deformation of the granitoid crust.

Cylindrical granitoid ultramylonite samples, composed of qtz + ab + K-fsp + bt + ep, with grain sizes of 125-15 μm are deformed in a Grigg’s type apparatus at T=650°C, confining P=1.2 GPa, strain rates=10^-3 to 10^-5s^-1, and 0.2 wt% H₂O added. Mechanical data are combined with light microscope, SEM, TEM, and quantitative image analysis to connect microstructures with stress and strain evolution. We show that polymineralic granitoid rocks deform through other mechanisms than monomineralic quartz aggregates at pressure and temperature conditions characteristic for the middle crust: Ultra-fine grain size reduction down to <50nm is developed by nucleation and growth of new grains in a polymineralic mixture. Grain size remains small because of pinning processes. We therefore refer to the deformation mechanism as pinning-controlled dissolution-precipitation creep (P-DPC).

Furthermore, we establish a new constitutive equation for this P-DPC, based on an exponential diffusion creep flow law, to model our experiments and tackle the extrapolation to various natural conditions. This flow law is supported by the microstructural evidence for the deformation mechanisms. Extrapolations show that the shear zones of the granitoid middle crust may be magnitudes weaker than extrapolated so far, and deformation may occur at magnitudes faster rates. The brittle to viscous transition may be shifted to shallower levels. This may have implications for the seismogenic zone and/or stress fields below geothermal reservoirs. Most importantly, we show the necessity to take polymineralic rocks into consideration for various numerical model applications.

How to cite: Nevskaya, N., Berger, A., Stünitz, H., Zhan, W., Plümper, O., Ohl, M., and Herwegh, M.: Implications for the strength of the Earth’s middle crust from novel experiments on natural fine-grained granitoid rocks , EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-17191, https://doi.org/10.5194/egusphere-egu24-17191, 2024.