Regional Stress Research Articles

Paleostress reconstruction of the Chaling-Yongxing basin in the central South China Block: Bulk (areal or multisite) vs. site stress tensors

A critical problem in fault-slip data stress analysis is separating the heterogeneous fault-slip data into homogeneous subsets. Recent attempts to solve this problem use tectonostratigraphic criteria based on the age of the sites where site stress tensors are calculated with a best-fit stress inversion method. The present study calculates bulk (areal or multisite) stress tensors from fault-slip data published in the Chaling-Yongxing Basin (CYB) in central South China with the Tensor Ratio Method (TRM). These tensors are compared with published site stress tensors to examine whether bulk or site stress tensors define regional or local stress regimes and similar or not paleostress history. Two regional stress regimes for the Late Cretaceous-Paleogene period, describing an older Late Cretaceous stage with transpression-strike-slip (TRP-SS) tectonics (NNE-SSW contraction and WNW-ESE extension) and a younger Paleogene stage with radial-pure (RE-PE) extensional tectonics (WNW-ESE extension), are defined. These stages are interpreted to be transitive through a transtension (TRN) tectonics of similar WNW-ESE extension and responsible for the formation and evolution of the CYB under tectonic processes related to the India-Asia collision and the retreat of the (Paleo)Pacific plate. The present analysis shows that tectonostratigraphic criteria based on the age of the exposed rocks might be insufficient for the coherent correlation of stress tensors resolved in sites of limited size and, thus, for unraveling the paleostress history of a region.

Slip partitioning and seismogenic stress buildup adjacent to the junction of the Lenglongling fault and the Tuolaishan fault: Insights from kinematically constrained numerical simulations

Understanding fault kinematics plays a crucial role in estimating seismic hazards and linking the relationship between seismogenic stress buildup and tectonic deformation. The northwestern segment of the Qilian-Haiyuan fault zone (QHFZ) serves as a natural laboratory for investigating how the oblique block motion decomposes into a complex fault system. In this study, we constructed a 3D FEM model incorporating the active faults as contact discontinuities to obtain a continuous numerical image of fault slip rates. Our findings revealed that deformation partitioning within the context of oblique block convergence generates a spatially heterogeneous distribution of fault slip rates. The inferred strike-slip rates of the Tuolaishan fault and the Lenglongling fault (LLLF) range from approximately 2 to 4 mm/yr and 4 to 7 mm/yr, respectively. The Lenglongling fault's slip behavior exhibits a spatial transition, with a dominant strike-slip component dominating in the western part and gradually changing to a combination of reverse faulting and strike-slip in the eastern segment. The Northern Lenglongling fault (NLLLF), characterized by a listric shape, merges with the LLLF in the deep crust, resulting in the formation of a positive flower structure that enhances the accumulation of compressional seismogenic stress. The region with high shear strain rates within this area has experienced significant earthquakes, including the 2022 Menyuan MS 6.9 earthquake, the 2016 Menyuan MS 6.4 earthquake, and the 1986 Menyuan MS 6.4 earthquake. Our numerical study gained valuable insights into the relationship between slip partitioning, fault structures, and regional stress.

A Weighted k-Nearest-Neighbors-Based Spatial Framework of Flood Inundation Risk for Coastal Tourism—A Case Study in Zhejiang, China

Flood inundation causes socioeconomic losses for coastal tourism under climate extremes, progressively attracting global attention. Predicting, mapping, and evaluating the flood inundation risk (FIR) is important for coastal tourism. This study developed a spatial tourism-aimed framework by integrating a Weighted k Nearest Neighbors (WkNN) algorithm, geographic information systems, and environmental indexes, such as precipitation and soil. These model inputs were standardized and weighted using inverse distance calculation and integrated into WkNN to infer the regional probability and distribution of the FIR. Zhejiang province, China, was selected as a case study. The evaluation results were mapped to denote the likelihood of an FIR, which was then validated by the historical Maximum Inundation Extent (MIE) extracted from the World Environment Situation Room. The results indicated that 80.59% of the WkNN results reasonably confirmed the MIE. Among the matched areas, 80.14%, 90.13%, 65.50%, and 84.14% of the predicted categories using WkNN perfectly coincided with MIE at high, medium, low, and very low risks, respectively. For the entire study area, approximately 2.85%, 64.83%, 10.8%, and 21.51% are covered by a high, medium, low, and very low risk of flood inundation. Precipitation and elevation negatively contribute to a high-medium risk. Drainage systems positively alleviate the regional stress of the FIR. The results of the evaluation illustrate that in most inland areas, some tourism facilities are located in high-medium areas of the FIR. However, most tourism facilities in coastal cities are at low or very low risk, especially from Hangzhou-centered northern coastal areas to southern Wenzhou areas. The results can help policymakers make appropriate strategies to protect coastal tourism from flood inundation. Moreover, the evaluation accuracy of WkNN is higher than that of kNN in FIR. The WkNN-based framework provides a reasonable method to yield reliable results for assessing FIR. The framework can also be extended to other risk-related research under climate change.

3-D crustal shear wave velocity and azimuthally anisotropic structure of the Pearl River Delta onshore-offshore area and its tectonic implications

The Pearl River Delta (PRD) onshore-offshore area lies in the transition zone between South China and northern continental margin of the South China Sea. To better understand the tectonic and deformation patterns at depth, we developed a 3-D high-resolution crustal shear wave velocity (Vs) and azimuthal anisotropy model by ambient noise tomography in this region. The data sets comprise ∼30 days of seismic ambient noise recorded by 124 land seismic stations and 28 Ocean Bottom Seismometers (OBS). Our results reveal a significant lateral contrast in azimuthal anisotropy in the crust. The fast axes of the azimuthal anisotropy in land show an approximately NW-SE direction, which may be mainly controlled by the regional stress field. In the offshore PRD, however, the fast axes of anisotropy strike nearly in the NE-SW direction, likely dominated by the shearing of the littoral fault zone. Combined with other previous results, the anisotropic model reveals three major crustal deformation processes in this region since the Yanshanian. Furthermore, an upper-middle crustal frozen mafic magma chamber with high Vs (∼3.85–4.0 km/s) is imaged beneath Shenzhen, Zhuhai, and Hong Kong. It is interpreted as a source of intense magmatic activities in and around Hong Kong during the late Mesozoic. Our new model provides important constraints for understanding the crustal architecture and the multistage deformation processes in and around the PRD onshore-offshore area.

Deep Seismic Reflection Imaging of Mesozoic Kachchh Rift, NW India: Implications for Evolution

AbstractDeep crustal seismic reflection profiling is carried out for first time across various tectonic domains of seismically active Mesozoic Kachchh rift basin, formed during the breakup of Gondwanaland. The seismic data, processed using the common reflection surface stack approach, provided maiden images of the shallow and deep sub‐surface structures in the region. These images reveal a 15 km thick subhorizontal lower crustal reflection fabric and crustal‐scale domal‐type structure extending from the surface to the Moho. We interpret the earlier structure represent magmatic underplating and the latter as the Kachchh Mainland uplift. We find large variations in the thickness of sediments from 150 m to 6.5 km and crustal thickness from 45 to 35 km from north to south, with a Moho up warp of 4 km beneath the Kachchh Mainland fault. The Kachchh rift basin exhibits an unusually thick crust of 45 km, contrary to many rift basins. We interpret the syn‐rifting and the Reunion mantle plume activity, manifested as Deccan volcanics, are responsible for magmatic underplating and crustal thickening. Uplift in the region is multi‐genetic in origin. Present study illuminated new faults and nature of various other faults. Moderate to large earthquakes in the region are attributed to the regional and local stresses resulting from the plate boundary and heterogeneous crustal structure. Based on the mafic lower crust and distribution of aftershocks through the entire crust, we interpret the lower crust is brittle, contrary to most models of continental rheology.

Tectonic stress regime and faulting style analysis due to sequence of earthquakes in Western Nepal

The current study examines the source mechanism and the regional stress regime of a sequence of moderate-sized earthquakes between October 2022 and October 2023. ISOLA-GUI-a Matlab coded program has been used to estimate Moment Tensor Solution to understand the source mechanisms of the earthquakes in the region. The earthquakes are recorded at the broadband station arrays of the regional seismological network around the Tehri region, installed by the Indian Institute of Technology Roorkee in India. The study focuses on earthquakes with epicentral distances of less than 10°. The investigation aims to comprehend the source process by estimating the focal mechanisms utilizing regional waveform data and then inverting the acquired focal mechanism of these earthquakes to comprehend the active stress regime in the region. The obtained source mechanism reveals a predominance of a double-couple mechanism in faulting and hence, suggests that the majority of the seismic energy is released through shear motion along a defined fault plane. By revealing the complex interaction of geological forces in the Himalayan belt, this investigation will help us to better understand the distribution of stress and deformation patterns beneath the Earth's crust. The main emphasis lies on assessing stress distribution through the utilization of earthquake focal mechanism data and employing the Stress Inversion technique using the WinTensor program. The result of stress inversion suggests that the maximum compressive stress is plunging horizontally, while the minimum compressive stress is oriented vertically. Hence revealing a thrust stress regime in Western Nepal, trending in the NNE- SSW direction corroborating the ongoing tectonic processes associated with the collision of the Indian plate with the Eurasian plate.

Source Parameters of Strong Turkish Earthquakes on February 6, 2023 (<i>M<sub>w</sub></i> = 7.8 and <i>M<sub>w</sub></i> = 7.7) from Surface Wave Data

Abstract—Based on the amplitude spectra of surface waves, the source parameters of the strong Turkish earthquakes of February 6, 2023 (Mw = 7.8 and Mw = 7.7) were calculated in two approximations: an instantaneous point source and an elliptical shear dislocation. As a result, rupture planes were identified, data were obtained on the scalar seismic moment, moment magnitude, focal mechanism, and source depth of the considered seismic events, and the integral parameters characterizing the rupture geometry and its development in time were estimated. It is shown that the sources of the earthquakes under study were formed under the influence of the regional stress field and their focal mechanisms were left lateral faults with a strike direction close to the strike of the East Anatolian fault zone for the first event and close to the strike of the Sürgü-Çardak fault system for the second. For the first earthquake, our estimates of the rupture duration and its length (t = 52.5 s, L = 180 km) probably refer not to the entire rupture, but only to its main phase, confined to the northeastern segments of the East Anatolian fault and characterized by maximum displacements and values of the released seismic moment. The values of t = 30 s and L = 180 km that we obtained for the second earthquake fully characterize the entire rupture.

Bedrock Fractures Control Groundwater‐Driven Mountain Slope Deformations

AbstractSeasonal deformation of mountain rock slopes can be driven by groundwater infiltration and depletion. Such processes could explain our field observation in the Aletsch Valley, Switzerland, where GNSS‐derived 3D annual displacement amplitudes reach 3.4 cm. However, the physical mechanisms behind such groundwater‐driven surface displacements are not well understood. Here, we develop a fully coupled hydromechanical model to simulate the relevant processes in a valley slope embedded with numerous fractures of variable sizes. The magnitude and orientation of transient annual slope surface displacement obtained from our model are in overall agreement with the field observations. The key geological factors controlling the type and magnitude of reversible mountain slope deformations are fracture network geometry, fracture aperture, and regional stress field. We show that the heterogeneity and anisotropy of bedrock hydromechanical responses, originating from depth‐dependent variations of fracture properties, play a critical role in groundwater recharge and valley slope deformation. During recharge events, pore pressure perturbations migrate downward from the groundwater table and toward the receiving stream and the deep subsurface. This process driven by pressure diffusion and poroelastic stressing develops in the subsurface with a great reach of up to a few kilometers, called critical hydromechanical response zone, and controls surface deformation patterns. During groundwater recession, this hydromechanical response zone expands downward and ground surface displacement vectors rotate upwards. Our results suggest that slope surface deformation can inform about subsurface permeability structures and pore pressure fluctuations, which have important implications for understanding groundwater flow in fractured bedrock slopes.

Petrological and geochemical evidence for a hot crystallization path and a recharge filtering bypass at Antimilos, Milos volcanic field, Greece

Antimilos volcano in the South Aegean Volcanic Arc, Greece, comprises an andesite–dacite suite that follows a distinct evolutionary path than the main edifice of the Milos volcanic field, despite their proximity. Petrographic and geochemical analyses reveal that basaltic andesite to low-Si dacite lavas have similar phenocryst assemblages that indicate crystallization from hot, relatively dry magmas in an upper crustal storage region. Rare antecrystic high-Mg# clinopyroxene cores with low Y, low Dy, and high Sr contents record the cryptic involvement of amphibole, a phase nominally absent from the erupted products, in the deeper parts of the plumbing system. Low temperature antecrysts with textures recording various degrees of disequilibrium suggest a protracted history of interaction between the upper crustal reservoir and deeper mafic melts, forming mobile hybrid magmas that consequently erupt as highly mingled, crystal-rich lava domes. Antimilos magmas seem to have escaped recharge filtering in the upper crust and prolonged stalling, which is the process that is probably responsible for the paucity of mafic eruptions in the rest of the Milos volcanic system. Large extensional structures offshore of Antimilos promote rapid ascent of mafic melts, inhibiting prolonged stalling and interaction with the arc crust. This model highlights the dominant role of the regional stress field in generating petrologically distinct suites in the marginal parts of some volcanic fields.

Regional evaluation method of ground stress in shale oil reservoirs-taking the Triassic Yanchang formation in northern Shaanxi area as an example

Estimating in-situ stress in shale oil reservoirs is critical for predicting hydraulic fracture height. However, the varied longitudinal lithology and transverse isotropy challenge the regional evaluation of ground stress. A new method based on the Kaiser effect method and the Sn Model was proposed to evaluate regional ground stress from point to point. Kaiser tests obtained the regional tectonic stress coefficients. Combined with the dynamic and static transformation of elastic parameters, the Sn ground stress calculation model suitable for the Northern Shaanxi Area was constructed. According to the calculation results of ground stress in the Chang 7 section of 600 wells in this region, the plane distribution maps of maximum horizontal ground stress (SH) and the difference between SH and Minimum horizontal stress (Sh) were drawn. It was found that the SH in this area was between 20 and 60 MPa, and the So was between 2.6 and 8.8 MPa. Compared with the 3D finite element method simulation results, the forward modeling results coincided well with the inversion results in this research. The causes of the plane distribution of ground stress in the Chang 7 section were explained by regional tectonic stress and paleo topography. It was found that the movement of plate and shell and deposition were critical factors of regional ground stress distribution. The research results provide a reference for further integration of geologic and engineering in evaluating shale oil reservoirs.

Decoding the state of stress and fluid pathways along the Andean Southern Volcanic Zone

Decoding means decrypting a hidden message. Here, the encrypted messages are the state of stress, fluid pathways, and volcano tectonic processes occurring in volcanoes of the Andean Southern Volcanic Zone (SVZ). To decode these messages, we use earthquake focal mechanisms, fault slip data, and a Monte Carlo simulation that predicts potential pathways for magmatic and hydrothermal fluids. From this analysis, we propose that SVZ volcanoes have three end-member stress patterns: (i) Stress-A, a strike-slip regime coupled with the regional far-field tectonic stress; (ii) Stress-B, an extensional regime that may be promoted by volcanic edifice loading and upward pressure due to magma inflation occurring within the upper brittle-crust; and (iii) Stress-C, a local and transient fluid-driven stress rotated ~90 degrees from Stress-A. Notoriously, Stress-C pattern was observed in most volcanoes with historical eruptions. We propose that volcanoes presenting Stress-B are attractive geothermal targets, while Stress-C could be used as a predicting signal for impending eruptions.

Serial RV wall stress measurements: association with right ventricular function in repaired Tetralogy of Fallot patients

BackgroundOptimal timing of pulmonary valve replacement (PVR) in Tetralogy of Fallot (TOF) patients remains challenging. Ventricular wall stress is considered to be an early marker of right ventricular (RV) dysfunction.ObjectivesTo investigate the association of RV wall stresses and their change over time with functional parameters in TOF patients.MethodsTen TOF patients after surgical repair with moderate/severe pulmonary regurgitation were included. At two timepoints (median follow-up time 7.2 years), patient-specific computational biventricular models for wall stress assessment were created using CMR short-axis cine images and echocardiography-based RV pressures. RV ejection fraction (RVEF), NT-proBNP and cardiopulmonary exercise tests were used as outcome measures reflecting RV function. Associations between regional RV diastolic wall stress and RV function were investigated using linear mixed models.ResultsIncreased wall stress correlated with lower RV mass (rrm = −0.70, p = 0.017) and lower RV mass-to-volume (rrm = −0.80, p = 0.003) using repeated measures. Wall stress decreased significantly over time, especially in patients with a stable RVEF (p < 0.001). Higher wall stress was independently associated with lower RVEF, adjusted for left ventricular ejection fraction, RV end-diastolic volume and time since initial surgery (decrease of 1.27% RVEF per kPa increase in wall stress, p = 0.029) using repeated measurements. No association was found between wall stress, NT-proBNP, and exercise capacity.ConclusionsUsing a computational method to calculate wall stress locally in geometrically complex ventricles, we demonstrated that lower wall stress might be important to maintain ventricular function. RV wall stress assessment can be used in serial follow-up, and is potentially an early marker of impending RV dysfunction.

The controls of strike-slip faults on fracture systems: Insights from 3D seismic data in the central Tarim Basin, northwest China

The central Tarim Basin has gained wide attention for its petroleum reserves, especially the recent commercial discovery of carbonate reservoirs of the Ordovician age in the Shunbei oil and gas field. A systematic analysis is conducted based on 3D seismic interpretation in the central Tarim Basin. The results indicate the presence of several major strike-slip faults and associated fracture systems. The characteristics of major strike-slip faults indicate a lower positive or negative flower structure in the Lower-Middle Ordovician and faults with a normal sense of movement in the Upper Ordovician due to the regional stress field. The major strike-slip faults commonly cause fractures in multiple scales. The fracture systems in different segments of major strike-slip faults contain significant differences in characteristics. In addition, two models are established with regard to strike-slip fault-associated fracture systems in a compressional setting and an extensional setting. The development of fracture systems varies to a large extent in relation to the scales of a strike-slip fractured zone. The fracture systems in a small-scale shear zone or strike-slip fractured zone usually develop along the fracture plane. In contrast, the fracture systems in a large-scale strike-slip fractured zone commonly develop along the fault zone on both sides.

TECTONIC STRESSES IN THE CHERSKY FAULT ZONE (BAIKAL RIFT SYSTEM

Based on the kinematic and structural-paragenetic approaches, there were performed the reconstructions of tectonic stresses in the Chersky fault zone. Rank analysis as applied to specialized mapping of fault zones was made for heterochronous rock complexes: Paleozoic, Upper Oligocene – Lower Pliocene and Upper Pliocene – Quaternary. There are distinguished three regional stress fields: Pre-Cenozoic compressive stress field with the NW-SE striking principal axis and two Cenozoic tensile stress fields with different axial orientations. The Cenozoic tensile stress field with a sub-meridionally striking minimum compression axis may be an earlier occurrence because it was found to exist in the Paleozoic and Oligocene – Early Pliocene rocks, and the tensile stress field, with the NW-SE axial orientation, may be later occurrence as reconstructed for all rock complexes including Quaternary pebbles. No evidence was found of tectonic shear stresses within the studied segment of the Chersky fault zone. Consideration is being given to probable reasons for the absence of residual shear deformations in the region and to the problems in the study of hierarchy of the upper crustal tectonic stresses.

The Impact of Pre-Existing Faults on Fault Geometry during Multiphase Rifts: The Jiyang Depression, Eastern China

The combination of multi-phase extension and pre-existing fault reactivation results in a complex fault pattern within hydrocarbon-bearing basins, affecting hydrocarbon exploration at different stages. We used high-resolution 3D seismic data and well data to reveal the impact of multi-phase extension and pre-existing fault reactivation on Cenozoic fault pattern changes over time in the Jiyang Depression of eastern China. The results show that during the Paleocene, a portion of NW-striking pre-existing faults reactivated under NS extension and controlled the basin structure (type 1). Other parts of the NW-striking pre-existing faults stopped activity and served as weak surfaces, and a series of NNE-striking faults were distributed in an en-echelon pattern along the NW direction at shallow depths (type 2). In areas unaffected by pre-existing faults, NE-striking faults formed perpendicular to regional stresses. During the Eocene, the regional stresses shifted clockwise to near-NS extension, and many EW-striking faults developed within the basin. The NE-striking faults and the EW-striking faults were hard-linked, forming the ENE-striking curved faults that controlled the structure in the basin (type 3). The NNE-striking faults were distinctly strike-slip at this time, with the ENE-striking faults forming a horsetail pattern at their tails. Many ENE-striking faults perpendicular to the extension direction were formed in areas where the basement was more stable and pre-existing faults were not developed (type 4). There were also developing NS-striking faults that were small in scale and appeared in positions overlapping different main faults (type 5). Additionally, different fault patterns can guide different phases of hydrocarbon exploration. Type 1, type 2, and type 3 faults are particularly suitable for early-stage exploration. In contrast, type 4 and type 5 faults are more appropriate for mature exploration areas, where they may reveal smaller hydrocarbon reservoirs.

Fracture characterisation using 3-D seismic reflection data for advanced deep geothermal exploration in the NE German Basin

At the geothermal research platform Gross Schönebeck (NE German Basin), we analysed 3-D seismic reflection data to determine the degree and direction of azimuthal velocity anisotropy which is interpreted as the effect of sub-vertical fracturing. Above the Zechstein salt, the observed anisotropy roughly correlates to fault structures formed by an upwelling salt pillow. Below the salt, faults are not obvious and the direction of less pronounced anisotropy and interpreted fracturing follows the trend of the regional stress field. The fracturing in an extensional setting above salt pillows may cause higher permeability and better conditions for geothermal exploitation.

Threshold effects of regional fiscal stress index on employment

Economic shocks and structural budget imbalances, when combined with ongoing uncertainty, can lead to fiscal stress in governments. This fiscal stress, along with the resulting volatility in the financing of local governments, can worsen their ability to meet short-term and long-term financial commitments and increase their reliance on the central government. Consequently, the effects of this stress, whether positive or negative, are closely tied to the actions and responses of both central and local governments. This highlights the crucial need for policymakers in central and local governments to respond accurately and promptly, through constant monitoring and assessment of fiscal stress indices. This study aims to illustrate the fiscal situation in the 31 provinces of Iran by calculating the local fiscal stress index based on fiscal and budgetary variables specific to each province. Furthermore, it seeks to estimate the threshold and spatial effects of this index on employment during the period of 2005–2017 using the panel smooth transition regression method. The findings reveal that initially, financial stress has an immediate and positive impact on employment. However, once the threshold of financial stress is crossed, and the subsequent pressures accumulate, the ability to control this imbalance diminishes, resulting in a decline in employment. Additionally, the ability or inability of local governments to manage income and expenses not only affects the economic indicators of the region but also spills over to neighbouring regions, leading to capital outflow and workforce migration, which are two major contributing factors to economic growth.

Earthquakes Induced by Rapid Loading of Faults During Pulichintala Reservoir Impoundment in the Stable Continental Region of India

AbstractThe first stage impoundment of the Pulichintala reservoir in the Proterozoic Cuddapah basin in the southeastern Indian shield began in August 2014. The second stage of impoundment started in 2019 when the height of the dam was increased and the filling of the reservoir was rapid. Soon after earthquakes started occurring which continued for a few months with a maximum magnitude of ML 4.6 on 25 January 2020. Two distinct clusters within the shallow crust were identified, which follow the faults/shear zones in the region. The estimated focal mechanisms with predominant strike slip motion on steep planes and the derived stress state are consistent with the regional stress and plate motion of the Indian plate. We simulated the influence of reservoir impoundment on the inferred seismogenic faults and found that the stress and pore pressure due to reservoir load indeed promoted failure, thus inducing earthquakes. We ascribe the induced earthquakes to rapid filling, longer and higher retention of high‐water levels during 2019–2020, which caused a stress change of 4–8 kPa. This stress change was at least two orders of magnitude higher than the typical tectonic loading rate in the stable continental region of the Indian Shield.

A study of earthquake potential analysis in central java and Yogyakarta provinces, Indonesia based on the calculation of a and B-Values with the maximum likelihood method

There are several active faults in Central Java and Yogyakarta (Indonesia), such as the Opak Fault, Baribis- Kendeng Fault, Merapi-Merbabu Fault, and others. In addition, Central Java and Yogyakarta are located near the subduction zone in the south of Java Island due to the collision of the Eurasian Plate and the Indo-Australian Plate. This has contributed to the fact that earthquakes can occur in Central Java and Yogyakarta. The purpose of this research is analyze the earthquake potential in Central Java, which includes an analysis of the a-value (seismicity level) and b-value (regional stress level) as an earthquake hazard mitigation effort. The Gutenberg-Richter equation expresses the relationship between earthquake frequency and magnitude as a straight-line equation containing the a-value and b-value parameters. The calculation of the b-value of the Gutenberg-Richter equation can be conducted with a modification of the maximum likelihood method. From 122 main earthquake events in the IRIS catalog, the values for depth (0-33) km are b = 1.01 ± 0.01 and a = 6.31, for depth (0-10) km, b = 1.44 ± 0.2 and a = 7.58; while for depth (10-33) km, b = 0.91 ± 0.1 and a = 5.73. The depth of (0-33) km and the depth of (0-10) km shows high a-value and b-value, while the depth of (10-33) km shows low a-value and b-value.

Variations and causes of in-situ stress orientations in the Dibei-Tuziluoke Gas Field in the Kuqa Foreland Basin, western China

The Kuqa Foreland Basin, adjacent to the southern Tianshan Mountains, is an important gas-producing basin in northwestern China. The development of the basin is dominated by a regional stress field with orientations of the maximum horizontal stress (SHmax) in the north-south direction, as a result of the collision between the Indian and Eurasian Plates. Previous studies have noted variations in the stress state in the Kuqa Foreland Basin; however, little attention has been paid to understanding the mechanisms causing such variations. In this study, a systematic analysis of in-situ SHmax orientations and fractures has been carried out using image logs from eight wells and XMAC/DSI logs from five wells. These wells are located in an East-West (E-W) trending fold-thrust in the foreland basin. The SHmax orientations (having quality ranked between A and D) obtained from borehole breakouts and drilling-induced fractures are highly variable and occur as two evidently different stress patterns on a scale of less than 10 km. The first stress pattern with SHmax orientations between NNW and NNE (165.08 and 188.59 °N) is broadly comparable to the regional stress field inferred from earthquake focal mechanism solutions in northwestern China from the World Stress Map (WSM). The stress pattern is interpreted to be controlled by plate boundary forces related to the collision between the Indian and Eurasian Plates. The SHmax orientations for the second stress pattern, however, ranging from NNE to ENE (26.65–57.77 °N), are clearly inconsistent with the regional stress filed and have been influenced by local geological factors. Faults and hinge-parallel fractures (trending in an approximately W-E direction) formed during the syn-folding stage probably resulted in heterogeneity of elastic properties and thus deflected SHmax orientations, both laterally and vertically, to generate the second stress pattern. The NE-striking hinge-oblique fractures are interpreted to have formed in response to the stress reorientation. Observations in this study suggest that there may exist a close relationship between fracturing (faulting) and stress variations, and pre-existing faults and fractures can influence subsequent fracturing by regulating local stress fields on a scale of several kilometers. The coupling between fracturing or faulting and stress produced fracture swarms, which can significantly enhance fluid flow and hence petroleum production.