Refine
Year of publication
Language
- English (19)
Is part of the Bibliography
- yes (19)
Keywords
- Database (2)
- Geomechanical Model (2)
- Long-Term (2)
- Middle Strand (2)
- NW Turkey (2)
- North Anatolian Fault (2)
- Northwestern Anatolia (2)
- Pull-Apart (2)
- Shear Zone (2)
- Slip Distribution (2)
- Stress pattern (2)
- Tectonic Evolution (2)
- Tectonic stress (2)
- in-situ stress (2)
- map project (2)
- 3D geomechanical numerical model (1)
- Alberta (1)
- Appalachian plateau (1)
- Borehole image logs (1)
- Canada (1)
- Chilean subduction zone (1)
- Circular statistics (1)
- Coulomb failure stress (1)
- GNSS (1)
- Geomechanical modelling (1)
- Iceland (1)
- InSAR (1)
- Induced seismicity (1)
- Inverse Theory (1)
- Jura (1)
- Kinematics of crustal (1)
- Lithospheric (1)
- New Zealand (1)
- North Anatolian Fault Zone (1)
- Rheological modeling (1)
- Sea of Marmara (1)
- Seismic gap (1)
- Seismicity (1)
- Statistical Methods (1)
- Stress field (1)
- Stress tensor (1)
- Taranaki Basin (1)
- Thermal modeling (1)
- afterslip inversion (1)
- alpine foreland (1)
- basin (1)
- deep-level mining (1)
- evolution (1)
- fault (1)
- fault reactivation (1)
- geomechanical model (1)
- geothermal-reservoirs (1)
- in situ stress (1)
- induced seismicity (1)
- insitu stress (1)
- km depth (1)
- mantle deformation (1)
- orientation (1)
- perturbations (1)
- plate tectonics (1)
- poroelasticity (1)
- power-law rheology (1)
- static stress change (1)
- strength (1)
- subduction zone (1)
- system (1)
- tabular mining (1)
- tectonics (1)
- topography (1)
Institute
The selection of earthquake focal mechanisms (FMs) for stress tensor inversion (STI) is commonly done on a spatial basis, that is, hypocentres. However, this selection approach may include data that are undesired, for example, by mixing events that are caused by different stress tensors when for the STI a single stress tensor is assumed. Due to the significant increase of FM data in the past decades, objective data-driven data selection is feasible, allowing more refined FM catalogues that avoid these issues and provide data weights for the STI routines. We present the application of angular classification with expectation-maximization (ACE) as a tool for data selection. ACE identifies clusters of FM without a priori information. The identified clusters can be used for the classification of the style-of-faulting and as weights of the FM data. We demonstrate that ACE effectively selects data that can be associated with a single stress tensor. Two application examples are given for weighted STI from South America. We use the resulting clusters and weights as a priori information for an STI for these regions and show that uncertainties of the stress tensor estimates are reduced significantly.
Knowledge of the present-day crustal in-situ stress field is a key for the understanding of geodynamic processes such as global plate tectonics and earthquakes. It is also essential for the management of geo-reservoirs and underground storage sites for energy and waste. Since 1986, the World Stress Map (WSM) project has systematically compiled the orientation of maximum horizontal stress (S-Hmax). For the 30th anniversary of the project, the WSM database has been updated significantly with 42,870 data records which is double the amount of data in comparison to the database release in 2008. The update focuses on areas with previously sparse data coverage to resolve the stress pattern on different spatial scales. In this paper, we present details of the new WSM database release 2016 and an analysis of global and regional stress pattern. With the higher data density, we can now resolve stress pattern heterogeneities from plate-wide to local scales. In particular, we show two examples of 40 degrees-60 degrees S-Hmax rotations within 70 km. These rotations can be used as proxies to better understand the relative importance of plate boundary forces that control the long wave-length pattern in comparison to regional and local controls of the crustal stress state. In the new WSM project phase IV that started in 2017, we will continue to further refine the information on the S-Hmax orientation and the stress regime. However, we will also focus on the compilation of stress magnitude data as this information is essential for the calibration of geomechanical-numerical models. This enables us to derive a 3-D continuous description of the stress tensor from point-wise and incomplete stress tensor information provided with the WSM database. Such forward models are required for safety aspects of anthropogenic activities in the underground and for a better understanding of tectonic processes such as the earthquake cycle.
Iceland is located on the Mid-Atlantic Ridge which is the plate boundary between the Eurasian and the North American plates. It is one of the few places on earth where an active spreading centre is located onshore but the stress pattern has not been extensively investigated so far. In this paper we present a comprehensive compilation of the orientation of maximum horizontal stress (S-Hmax). In particular we interpret borehole breakouts and drilling induced fractures from borehole image logs in 57 geothermal wells onshore Iceland. The borehole results are combined with other stress indicators including earthquake focal mechanism solutions, geological information and overcoring measurements resulting in a dataset with 495 data records for the S-Hmax orientation. The reliability of each indicator is assessed according to the quality criteria of the World Stress Map project The majority of S-Hmax orientation data records in Iceland is derived from earthquake focal mechanism solutions (35%) and geological fault slip inversions (26%). 20% of the data are borehole related stress indicators. In addition minor shares of S-Hmax orientations are compiled, amongst others, from focal mechanism inversions and the alignment of fissure eruptions. The results show that the S-Hmax orientations derived from different depths and stress indicators are consistent with each other.
The resulting pattern of the present-day stress in Iceland has four distinct subsets of S-Hmax orientations. The S-Hmax orientation is parallel to the rift axes in the vicinity of the active spreading regions. It changes from NE-SW in the South to approximately N-S in central Iceland and NNW-SSE in the North. In the Westfjords which is located far away from the ridge the regional S-Hmax rotates and is parallel to the plate motion. (C) 2016 Elsevier B.V. All rights reserved.
The stress field at depth is a relevant parameter for the design of subsurface constructions and reservoir management. Yet the distortion of the regional stress field due to local-scale features such as sedimentary and tectonic structures or topography is often poorly constrained. We conduct a stress sensitivity analysis using 3-D numerical geomechanical modelling with an elasto-plastic material law to explore the impact of such site-specific features on the stress field in a sedimentary sequence of the Swiss Alpine foreland. The model's dimensions are 14 x 14 x 3 km(3) and it contains 10 units with different mechanical properties, intersected by two regional fault zones. An initial stress state is established involving a semi-empirical relationship between the ratio of horizontal to vertical stress and the overconsolidation ratio of argillaceous sediments. The model results indicate that local topography can affect the stress field significantly to depths greater than the relief contrasts at the surface, especially in conjunction with horizontal tectonic loading. The complexity and frictional properties of faults are also relevant. The greatest variability of the stress field arises across the different sedimentary units. Stress magnitudes and stress anisotropy are much larger in stiffer formations such as massive limestones than in softer argillaceous formations. The stiffer formations essentially carry the load of the far-field forces and are therefore more sensitive to changes of the boundary conditions. This general characteristic of stress distribution in the stiff and soft formations is broadly maintained also with progressive loading towards the plastic limit. The stress field in argillaceous sediments within a stack of formations with strongly contrasting mechanical properties like in the Alpine foreland appears to be relatively insensitive to changes in the tectonic boundary conditions and is largely controlled by the maximum stiffness contrast with respect to the load-bearing formations.
Megathrust earthquakes impose changes of differential stress and pore pressure in the lithosphere-asthenosphere system that are transiently relaxed during the postseismic period primarily due to afterslip, viscoelastic and poroelastic processes.
Especially during the early postseismic phase, however, the relative contribution of these processes to the observed surface deformation is unclear.
To investigate this, we use geodetic data collected in the first 48 days following the 2010 Maule earthquake and a poro-viscoelastic forward model combined with an afterslip inversion.
This model approach fits the geodetic data 14% better than a pure elastic model. Particularly near the region of maximum coseismic slip, the predicted surface poroelastic uplift pattern explains well the observations.
If poroelasticity is neglected, the spatial afterslip distribution is locally altered by up to +/- 40%.
Moreover, we find that shallow crustal aftershocks mostly occur in regions of increased postseismic pore-pressure changes, indicating that both processes might be mechanically coupled.
The surface deformation associated with the 2010 M-w 8.8 Maule earthquake in Chile was recorded in great detail before, during and after the event. The high data quality of the continuous GPS (cGPS) observations has facilitated a number of studies that model the postseismic deformation signal with a combination of relocking, afterslip and viscoelastic relaxation using linear rheology for the upper mantle. Here, we investigate the impact of using linear Maxwell or power-law rheology with a 2D geomechanical-numerical model to better understand the relative importance of the different processes that control the postseismic deformation signal. Our model results reveal that, in particular, the modeled cumulative vertical postseismic deformation pattern in the near field (< 300 km from the trench) is very sensitive to the location of maximum afterslip and choice of rheology. In the model with power-law rheology, the afterslip maximum is located at 20-35 km rather than > 50 km depth as suggested in previous studies. The explanation for this difference is that in the model with power-law rheology the relaxation of coseismically imposed differential stresses occurs mainly in the lower crust. However, even though the model with power-law rheology probably has more potential to explain the vertical postseismic signal in the near field, the uncertainty of the applied temperature field is substantial, and this needs further investigations and improvements.
Postseismic surface deformation associated with great subduction earthquakes is controlled by asthenosphere rheology, frictional properties of the fault, and structural complexity. Here by modeling GPS displacements in the 6 years following the 2010 M-w 8.8 Maule earthquake in Chile, we investigate the impact of heterogeneous viscosity distribution in the South American subcontinental asthenosphere on the 3-D postseismic deformation pattern. The observed postseismic deformation is characterized by flexure of the South America plate with peak uplift in the Andean mountain range and subsidence in the hinterland. We find that, at the time scale of observation, over 2 orders of magnitude gradual increase in asthenosphere viscosity from the arc area toward the cratonic hinterland is needed to jointly explain horizontal and vertical displacements. Our findings present an efficient method to estimate spatial variations of viscosity, which clearly improves the fitting to the vertical signal of deformation. Lateral changes in asthenosphere viscosity can be correlated with the thermomechanical transition from weak subvolcanic arc mantle to strong subcratonic mantle, thus suggesting a stationary heterogeneous viscosity structure. However, we cannot rule out a transient viscosity structure (e.g., power law rheology) with the short time span of observation.
On 27 December 2007, a 1.9 seismic event occurred within a dyke in the deep-level Mponeng Gold Mine, South Africa. From the seismological network of the mine and the one from the Japanese-German Underground Acoustic Emission Research in South Africa (JAGUARS) group, the hypocentral depth (3,509 m), focal mechanism and aftershock location were estimated. Since no mining activity took place in the days before the event, dynamic triggering due to blasting can be ruled out as the cause. To investigate the hypothesis that stress transfer, due to excavation of the gold reef, induced the event, we set up a small-scale high-resolution three-dimensional (3D) geomechanical numerical model. The model consisted of the four different rock units present in the mine: quartzite (footwall), hard lava (hanging wall), conglomerate (gold reef) and diorite (dykes). The numerical solution was computed using a finite-element method with a discretised mesh of approximately elements. The initial stress state of the model is in agreement with in situ data from a neighbouring mine, and the step-wise excavation was simulated by mass removal from the gold reef. The resulting 3D stress tensor and its changes due to mining were analysed based on the Coulomb failure stress changes on the fault plane of the event. The results show that the seismic event was induced regardless of how the Coulomb failure stress changes were calculated and of the uncertainties in the fault plane solution. We also used the model to assess the seismic hazard due to the excavation towards the dyke. The resulting curve of stress changes shows a significant increase in the last in front of the dyke, indicating that small changes in the mining progress towards the dyke have a substantial impact on the stress transfer.
The Sea of Marmara is a tectonically active basin that straddles the North Anatolian Fault Zone (NAFZ), a major strike-slip fault that separates the Eurasian and Anatolian tectonic plates. The Main Marmara Fault (MMF), which is part of the NAFZ, contains an approximately 150 km long seismotectonic segment that has not ruptured since 1766. A key question for seismic hazard and risk assessment is whether or not the next rupture along this segment is likely to produce one major earthquake or a series of smaller earthquakes. Geomechanical characteristics such as along-strike variations in rock strength may provide an important control on seismotectonic segmentation. We find that variations in lithospheric strength throughout the Marmara region control the mechanical segmentation of the MMF and help explain its long-term seismotectonic segmentation. In particular, a strong crust that is mechanically coupled to the upper mantle spatially correlates with aseismic patches, where the MMF bends and changes its strike in response to the presence of high-density lower crustal bodies. Between the bends, mechanically weaker crustal domains that are decoupled from the mantle indicate a predominance of creeping. These results are highly relevant for the ongoing debate regarding the characteristics of the Marmara seismic gap, especially in view of the seismic hazard (Mw > 7) in the densely populated Marmara region.
Rotations of the principal stress axes are observed as a result of fluid injection into reservoirs. We use a generic, fully coupled 3-D thermo-hydro-mechanical model to investigate systematically the dependence of this stress rotation on different reservoir properties and injection scenarios. We find that permeability, injection rate, and initial differential stress are the key factors, while other reservoir properties only play a negligible role. In particular, we find that thermal effects do not significantly contribute to stress rotations. For reservoir types with usual differential stress and reservoir treatment the occurrence of significant stress rotations is limited to reservoirs with a permeability of less than approximately 10(-12)m(2). Higher permeability effectively prevents stress rotations to occur. Thus, according to these general findings, the observed principal stress axes rotation can be used as a proxy of the initial differential stress provided that rock permeability and fluid injection rate are known a priori.