Low B-values Research Articles

Morphometry and active tectonics of the Konkan coast, western India

The Indian western coast is petroliferous drawing attention from academicians and industrial experts. In this work, geomorphic indices have been calculated to decode neotectonics of the Konkan onshore region on India’s west coast. Earthquake and Bouguer anomaly data have been used along with the present-day stress data from the World Stress Map project. Gravity modeling was performed in order to gain seismotectonic insights. The b-values were determined using Z-MAP 7.1 (2021). Morphometric analysis at both linear and spatial scales were performed using the digital elevation model from the data derived from the Shuttle Radar Topography Mission, analysed with ArcGIS 10.3 (2014) software. Maps depicting slopes, aspects, and reliefs were created. NW-SE trending lineaments in the Konkan plain guided the major stream courses. Two of the five watersheds, watersheds 4 and 5, reveal high tectonic activity, are landslide-prone and host hot springs. Interestingly, watersheds 4 and 5 show high b-values (except near the rivers’ sources), low Bouguer anomalies, and higher Hypsometric integral values (0.18523 and 0.16698) than the other watersheds. A low b-value in watershed 3 indicates stress accumulation. Over a larger area, the gravity trend varies from ∼ -80 to 30 mGal. The lineaments diagram deduced from the first vertical derivative technique shows that the structural fabrics mostly trend ∼ NW-SE at the west of the Western Ghat Escarpment (WGE) while it is NE-SW at the east. The tilt derivative ratio technique reveals a major NE-SW trend to the west of the WGE and an E-W trend to the east. Structural interpretations based on drill-cores around Koyna combined with geophysical studies for deep crust will be required are required for a better understanding of the blind (active) structures in the region.

Microseismicity and fault structure in the Daliangshan block within the Southeastern Tibetan Plateau

The Daliangshan block is located between the Tibetan Plateau and the South China block and has accommodated several M > 6.5 damaging earthquakes in the past ∼600 years, as well as intense tectonic deformation and complex fault structures. In this study, we analyze more than two years of continuous seismic data recorded by a recently deployed dense seismic array. We used a recently developed machine learning-based earthquake location workflow (ESPRH) to construct a high-precision earthquake catalog for the region and obtained 3539 earthquakes, which is approximately three times as many as the National Earthquake Data Center (NEDC) catalog contains. The seismicity distribution not only confirms the nature of the faults marked on the map but also delineates the detailed geometry of the unmapped faults, including the en échelon faults at the northern end of the Zemuhe Fault and the “V”-shaped conjugate fault within the Mabian Fault Zone. The Zemuhe Fault and Lianfeng Fault are prone to hosting large earthquakes according to the derived low b-value. The western side of the Daliangshan block is dominated by strike-slip faults. Combining the fault geometry presented in this paper, we observed that the fault properties on the eastern side are complex. This tectonic phenomenon is attributed to the fact that during the lateral extrusion of the southeastern edge of the Chuandian fragments, the northeastern part of the Daliangshan block was squeezed by the South China block more strongly than its southwestern part. We provide the first precise earthquake catalog for Daliangshan block, which can be used as important seismological data for regional hazard assessment and research on the southeastern (SE) margin of the Tibetan Plateau.

Spatial–temporal variations of b-values prior to medium-to-large earthquakes in Taiwan and the feasibility of real-time precursor monitoring

This study explores the b-value variations prior to M ≥ 6.0 earthquakes in Taiwan, examining their potential as earthquake precursors. Focusing on the 2018 Hualien earthquake and others between 1999 and 2021, we found that many large earthquakes occurred in areas with low b-values a year prior, although there were no significant temporal changes near the epicenters. However, for more accurate earthquake precursors, incorporating additional factors is recommended to minimize uncertainty.Graphical

Source parameters and aftershock pattern of the October 7, 2021, M5.9 Harnai earthquake, Pakistan

On October 7, 2021, a magnitude 5.9 earthquake struck the Harnai (Baluchistan) region of Pakistan, causing several fatalities and injuries within the epicentral area. First-order tectonic deformation in this region is caused by the convergence of the Indian Plate with respect to the Eurasian Plate. The Katwaz Block hinders the motion of the Indian Plate, resulting in the formation of strike-slip faults. In this study, the P-wave first-motion polarity technique was used to determine the mainshock faulting style. Cyclic scanning of the polarity solutions was applied to determine the most suitable focal mechanism solution among the available solutions generated by the FOCMEC (focal mechanism) software. The nodal planes correspond to different faulting styles (i.e., thrust and strike-slip faulting). A nodal plane oriented in the NW-SE direction corresponded to a strike-slip mechanism, which was considered to be the fault plane. Tectonically, this earthquake was associated with the Harnai-Karahi strike-slip fault zone owing to the fault strike and direction of slip. The apparent stress drop, fault length, and moment magnitude of the Harnai earthquake were 35.4 bar, 6.1 km, and 5.9, respectively. A lower b-value for the Gutenberg-Richter law was observed prior to the earthquake. Higher α- than b-values (α >b) indicate that this earthquake was governed by large events as opposed to small-magnitude events. The Harnai sequence had a decay exponent close to unity, lasted for 145 days, and produced few aftershocks. The study will help the future hazard mitigation in the region.

Stable estimation of the Gutenberg–Richter b-values by the b-positive method: a case study of aftershock zones for magnitude-7 class earthquakes

The traditional approach for estimating the b-value of the Gutenberg–Richter law, which is posited to inversely correlate with differential stress, has historically relied on the maximum likelihood technique, utilizing data from earthquakes exceeding a magnitude cutoff, Mc. This traditional approach is significantly influenced by the value of Mc, leading to extensive research focused on methods for determining Mc with greater accuracy. However, a recent study introduced a novel method based on the frequency distribution of magnitude difference, termed the b-positive method. This innovative method could enable more robust b-value estimations, even in scenarios where Mc may vary spatially and temporally. Our study concentrated on analyzing aftershocks, related to 25 magnitude-7 class earthquakes surrounding the Japanese archipelago. We estimated the b-values using both the goodness-of-fit test for Mc, a traditional approach, and the b-positive method. The aftershock data were examined over two distinct time frames: the initial 10 days following each mainshock and an extended period of 1000 days. Our findings indicated that the estimates produced by the b-positive method showed negligible variation between the 10-day and 1000-day aftershock periods (correlation coefficient of 0.95), whereas the traditional approach tended to yield lower b-values for the 10-day aftershocks compared to those from the 1000-day period. Variations in b-values, when analyzed using the traditional approach, could be inaccurately ascribed to temporal fluctuations in differential stress that may not actually be present. The b-positive method offers a vital solution to prevent these erroneous interpretations, serving as an essential alternative.Graphical abstract

Precise relative magnitude measurement improves fracture characterization during hydraulic fracturing

SUMMARY Microseismic monitoring is an important technique to obtain detailed knowledge of in-situ fracture size and orientation during stimulation to maximize fluid flow throughout the rock volume and optimize production. Furthermore, considering that the frequency of earthquake magnitudes empirically follows a power law (i.e. Gutenberg–Richter), the accuracy of microseismic event magnitude distributions is potentially crucial for seismic risk management. In this study, we analyse microseismicity observed during four hydraulic fracture treatments of the legacy Cotton Valley experiment in 1997 at the Carthage gas field of East Texas, where fractures were activated at the base of the sand-shale Upper Cotton Valley formation. We perform waveform cross-correlation to detect similar event clusters, measure relative amplitude from aligned waveform pairs with a principal component analysis, then measure precise relative magnitudes. The new magnitudes significantly reduce the deviations between magnitude differences and relative amplitudes of event pairs. This subsequently reduces the magnitude differences between clusters located at different depths. Reduction in magnitude differences between clusters suggests that some attenuation-related biases could be effectively mitigated with relative magnitude measurements. The maximum likelihood method is applied to understand the magnitude frequency distributions and quantify the seismogenic index of the clusters. Statistical analyses with new magnitudes suggest that fractures that are more favourably oriented for shear failure have lower b-value and higher seismogenic index, suggesting higher potential for relatively larger earthquakes, rather than fractures subparallel to maximum horizontal principal stress orientation.

Arterial pulsation dependence of perivascular cerebrospinal fluid flow measured by dynamic diffusion tensor imaging in the human brain

Perivascular cerebrospinal fluid (pCSF) flow is a key component of the glymphatic system. Arterial pulsation has been proposed as the main driving force of pCSF influx along the superficial and penetrating arteries; however, evidence of this mechanism in humans is limited. We proposed an experimental framework of dynamic diffusion tensor imaging with low b-values and ultra-long echo time (dynDTIlow-b) to capture pCSF flow properties during the cardiac cycle in human brains. Healthy adult volunteers (aged 17–28 years; seven men, one woman) underwent dynDTIlow-b using a 3T scanner (MAGNETOM Prisma, Siemens Healthcare, Erlangen, Germany) with simultaneously recorded cardiac output. The results showed that diffusion tensors reconstructed from pCSF were mainly oriented in the direction of the neighboring arterial flow. When switching from vasoconstriction to vasodilation, the axial and radial diffusivities of the pCSF increased by 5.7 % and 4.94 %, respectively, suggesting that arterial pulsation alters the pCSF flow both parallel and perpendicular to the arterial wall. DynDTIlow-b signal intensity at b=0 s/mm2 (i.e., T2-weighted, [S(b=0 s/mm2)]) decreased in systole, but this change was ∼7.5 % of a cardiac cycle slower than the changes in apparent diffusivity, suggesting that changes in S(b=0 s/mm2) and apparent diffusivity arise from distinct physiological processes and potential biomarkers associated with perivascular space volume and pCSF flow, respectively. Additionally, the mean diffusivities of white matter showed cardiac-cycle dependencies similar to pCSF, although a delay relative to the peak time of apparent diffusivity in pCSF was present, suggesting that dynDTIlow-b could potentially reveal the dynamics of magnetic resonance imaging-invisible pCSF surrounding small arteries and arterioles in white matter; this delay may result from pulse wave propagation along penetrating arteries. In conclusion, the vasodilation-induced increases in axial and radial diffusivities of pCSF and mean diffusivities of white matter are consistent with the notion that arterial pulsation can accelerate pCSF flow in human brain. Furthermore, the proposed dynDTIlow-b technique can capture various pCSF dynamics in artery pulsation.

The Stress State before the MS 6.8 Luding Earthquake on 5 September 2022 in Sichuan, China: A Retrospective View Based on the b-Value

On 5 September 2022 (BJT), Luding, located in southwestern Sichuan Province, China, experienced an MS 6.8 earthquake. This earthquake occurred within the historical rupture zone of the 1786 MS 7.75 event, part of the southern section of the Xianshui He Fault belt. Given the average 155-year recurrence interval for strong earthquakes in this area, the 236 years since the last event made this earthquake somewhat expected. However, prior to this event, we did not detect any anomalies indicating low surface b-values, which are often indicative of a high-stress state in the source area before strong earthquakes, as highlighted by numerous studies. Our research focused on the northern section of the eastern boundary of the Sichuan–Yunnan sub-block, encompassing the Xianshui He, Anning He, Zemu He, and Daliang Shan fault belts. We meticulously located earthquakes of ML ≥ 1.5 from 2009 to May 2022. The catalog was divided into two periods: 2009–2014 and 2015–May 2022. Using an AIC-constraint method, we analyzed the changes in b-values (Δb) in the latter period compared to the former. Our findings revealed a significant abnormal Δb zone (Δb < −0.3), with a radius of approximately 50 km, when ΔAIC ≥ 2 was selected. Intriguingly, the epicenter of the recent Luding MS 6.8 earthquake fell within this abnormal zone. Furthermore, we calculated the b-value cross-section for the southern section of the Xianshui He fault belt using a directory of precisely located small earthquakes. This revealed that the location, scale, and shape of the abnormally low-b-value area corresponded with the large displacement co-seismic area of the main earthquake, affirming the b-value’s effectiveness in identifying asperities. The b-value’s temporal evolution prior to the mainshock exhibited a nearly decade-long continuous decrease, signifying a long-term stress-loading process akin to that observed before many strong earthquakes. The b-value anomalies observed from different profiles before the Luding earthquake underline the necessity of a comprehensive, multi-dimensional analysis of such anomalies. Finally, our analysis indicates that nine earthquakes with MS ≥ 6.5, including the Luding MS 6.8 event, have contributed to increased Coulomb Failure Stress change (ΔCFS) in the Daofu (DF)–Kangding (KD) section of the Xianshui He fault belt and the northern section of the Anning He fault belt south of Shimian (SM), with amplitudes surpassing the 0.01 MPa threshold. This suggests the potential for strong earthquakes in these zones.

Reduced cerebrospinal fluid motion in patients with Parkinson’s disease revealed by magnetic resonance imaging with low b-value diffusion weighted imaging

BackgroundParkinson’s disease is characterized by dopamine-responsive symptoms as well as aggregation of α-synuclein protofibrils. New diagnostic methods assess α-synuclein aggregation characteristics from cerebrospinal fluid (CSF) and recent pathophysiologic mechanisms suggest that CSF circulation disruptions may precipitate α-synuclein retention. Here, diffusion-weighted MRI with low-to-intermediate diffusion-weightings was applied to test the hypothesis that CSF motion is reduced in Parkinson’s disease relative to healthy participants.MethodsMulti-shell diffusion weighted MRI (spatial resolution = 1.8 × 1.8 × 4.0 mm) with low-to-intermediate diffusion weightings (b-values = 0, 50, 100, 200, 300, 700, and 1000 s/mm2) was applied over the approximate kinetic range of suprasellar cistern fluid motion at 3 Tesla in Parkinson’s disease (n = 27; age = 66 ± 6.7 years) and non-Parkinson’s control (n = 32; age = 68 ± 8.9 years) participants. Wilcoxon rank-sum tests were applied to test the primary hypothesis that the noise floor-corrected decay rate of CSF signal as a function of b-value, which reflects increasing fluid motion, is reduced within the suprasellar cistern of persons with versus without Parkinson’s disease and inversely relates to choroid plexus activity assessed from perfusion-weighted MRI (significance-criteria: p < 0.05).ResultsConsistent with the primary hypothesis, CSF decay rates were higher in healthy (D = 0.00673 ± 0.00213 mm2/s) relative to Parkinson’s disease (D = 0.00517 ± 0.00110 mm2/s) participants. This finding was preserved after controlling for age and sex and was observed in the posterior region of the suprasellar cistern (p < 0.001). An inverse correlation between choroid plexus perfusion and decay rate in the voxels within the suprasellar cistern (Spearman’s-r=-0.312; p = 0.019) was observed.ConclusionsMulti-shell diffusion MRI was applied to identify reduced CSF motion at the level of the suprasellar cistern in adults with versus without Parkinson’s disease; the strengths and limitations of this methodology are discussed in the context of the growing literature on CSF flow.

Radon concentration in seawater as a geochemical indicator of submarine fault activity in the Yatsushiro Sea, Japan

This study examined the relationship between radon (222Rn) concentrations in seawater and crustal activity in the Yatsushiro Sea by investigating the submarine fault zone situated at the southern end of the Futagawa–Hinagu fault zone, activated by the 2016 Kumamoto earthquake (M7.3). We conducted an analysis of 222Rn concentration in samples of bottom water just above the seafloor and pore water in sediments, utilizing multiple and piston cores from the Hakuho Maru Expedition KH18-3. The findings revealed significantly elevated 222Rn concentrations in the central sites of the Yatsushiro Sea, coinciding with a high-stress field exhibiting dense active faults. Seismicity analysis revealed heightened moment release and a low b-value post the 2016 Kumamoto earthquake, indicative of increased seismic activity and the potential for substantial earthquakes in the Yatsushiro Sea vicinity. Our results indicate that heightened concentrations of 222Rn in seawater can serve as an effective tracer for identifying and estimating submarine fault activities. Moreover, our research highlights the utility of 222Rn concentrations in detecting active submarine faults and assessing their activity. It contributes to a comprehensive understanding of the potential for significant earthquakes in the Yatsushiro Sea in the future.

PadGAN: An End-to-End dMRI Data Augmentation Method for Macaque Brain

Currently, an increasing number of macaque brain MRI datasets are being made publicly accessible. Unlike human, publicly accessible macaque brain datasets suffer from data quality in diffusion magnetic resonance imaging (dMRI) data. Typically, dMRI data require a minimum ratio of 1:10 between low b-value (b &lt; 10) volumes and high b-value (b &gt; 300) volumes. However, the currently accessible macaque datasets do not meet this ratio. Due to site differences in macaque brain images, traditional human brain image-to-image translation models struggle to perform well on macaque brain images. Our work introduces a novel end-to-end primary-auxiliary dual generative adversarial network (PadGAN) for generating low b-value images. The auxiliary generator in the PadGAN is responsible for extracting the latent space features from peak information maps and transmitting them to the primary generator, enabling the primary generator to generate images with rich details. Experimental results demonstrate that PadGAN outperforms existing methods both qualitatively and quantitatively (mean SSIM increased by 0.1139). Diffusion probabilistic tractography using dMRI data augmented by our method yields superior results.

A framework for optimizing the acquisition protocol multishell diffusion-weighted imaging for multimodel assessment.

Complementary aspects of tissue microstructure can be studied with diffusion-weighted imaging (DWI). However, there is no consensus on how to design a diffusion acquisition protocol for multiple models within a clinically feasible time. The purpose of this study is to provide a flexible framework that is able to optimize the shell acquisition protocol given a set of DWI models. Eleven healthy subjects underwent an extensive DWI acquisition protocol, including 15 candidate shells, ranging from 10 to 3500 s/mm2. The proposed framework aims to determine the optimized acquisition scheme (OAS) with a data-driven procedure minimizing the squared error of model-estimated parameters. We tested the proposed method over five heterogeneous DWI models exploiting both low and high b-values (i.e., diffusion tensor imaging [DTI], free water, intra-voxel incoherent motion [IVIM], diffusion kurtosis imaging [DKI], and neurite orientation dispersion and density imaging [NODDI]). A voxel-level and region of interest (ROI)-level analysis was conducted over the white matter and in 48 fiber bundles, respectively. Results showed that acquiring data for the five abovementioned models via OAS requires 14 min, compared with 35 min for the joint recommended acquisition protocol. The parameters derived from the reference acquisition scheme and the OAS are comparable in terms of estimated values, noise, and tissue contrast. Furthermore, the power analysis showed that the OAS retains the potential sensitivity to group-level differences in the parameters of interest, with the exception of the free water model. Overall, there is a linear correspondence (R2= 0.91) between OAS and reference-derived parameters. In conclusion, the proposed framework optimizes the shell acquisition scheme for a given set of DWI models (i.e., DTI, free water, IVIM, DKI, and NODDI), combining low and high b-values while saving acquisition time.

The relationship between GRACE gravity and the seismic b-value: a case study of the Northern Chile Triple Junction (25° S–40° S)

SUMMARY The northern Chile Triple Junction (CTJ) is characterized by the ongoing subduction of the Nazca plate beneath the South American plate. The geological structures within the subduction zone undergo complex changes, resulting in significant tectonic activities and intense seismicity along the western margin of South America. Based on the Gravity Recovery and Climate Experiment (GRACE) data and earthquake catalogues, this study selects the northern CTJ area (25° S–40° S, 75° W–65° W) as the research object, adopts the mathematical methods of independent component analysis (ICA) and principal component analysis (PCA) to separate the earthquake-related signals within the GRACE data, and fits the changes of seismic b-values through the frequency–magnitude relationship. The characteristics of gravity changes before and after seismic events, the seismic activity parameter b-values, and the relationship between the gravity signals and b-values are discussed. The results show that mathematical methods can effectively extract seismic-related gravity components from the GRACE data. ICA, compared to PCA, provides better results in capturing the temporal variations associated with b-value time-series, which exhibit good consistency in long-term trend changes. The average change of b-values in the study area is 0.66 ± 0.003, fluctuating over time. Generally, prior to larger seismic events, b-values tend to decrease. Along the western margin of South America, b-values are low; this aligns with the active tectonic activities between subducting plates. Additionally, a certain correlation between b-values and gravity changes is observed, but due to the influence of tectonic activities, the correspondence between b-values and gravity anomalies may not be consistent across different areas. The b-value is highly consistent with the strain rate model. Low b-values correspond to high strain rates along the western edge of South America, which is in line with the tectonic characteristics of frequent seismic activity in this area. A gradual concentration of gravity anomalies before major earthquakes is observed, accompanied by the gradual accumulation of smaller seismic events. Meanwhile, several months before the two major earthquakes, the spatial distribution of gravity appears to be similar to the coseismic signals, but the nature of its generation remains to be explored. These methods and results not only add to the applications of GRACE in seismic studies but also raise questions for further exploration.

Analysis of IVIM Perfusion Fraction Improves Detection of Pancreatic Ductal Adenocarcinoma.

The purpose of this study was to evaluate whether intravoxel incoherent motion (IVIM) parameters can enhance the diagnostic performance of MRI in differentiating normal pancreatic parenchyma from solid pancreatic adenocarcinomas. This study included 113 participants: 66 patients diagnosed with pancreatic adenocarcinoma and 47 healthy volunteers. An MRI was conducted at 1.5 T MR unit, using nine b-values. Postprocessing involved analyzing both conventional monoexponential apparent diffusion coefficient (ADC) and IVIM parameters (diffusion coefficient D-pure molecular diffusion coefficient, perfusion-dependent diffusion coefficient D*-pseudodiffusion coeffitient, and perfusion fraction coefficient (f)) across four different b-value selections. Significantly higher parameters were found in the control group when using high b-values for the pure diffusion analysis and all b-values for the monoexponential analysis. Conversely, in the study group, the parameters were affected by low b-values. Most parameters could differentiate between normal and cancerous tissue, with D* showing the highest diagnostic performance (AUC 98-100%). A marked decrease in perfusion in the patients with pancreatic cancer, indicated by the significant differences in the D* medians between groups, was found. In conclusion, standard ADC maps alone may not suffice for a definitive pancreatic cancer diagnosis, and incorporating IVIM into MRI protocols is recommended, as the reduced tissue perfusion detected by the IVIM parameters is a promising marker for pancreatic adenocarcinoma.

Spatio-Temporal Variation Seismicity Pattern in East Java Between 2002 and 2022 Based on the b-value and Seismic Quiescence z-value

East Java is one of the seismically active regions in part of the Sunda Arc whose seismicity pattern must be studied. Therefore, a spatio-temporal analysis of seismicity between 2002 - 2022 in East Java has been conducted. This analysis aims to investigate the seismotectonic conditions in the East Java region over the last twenty years, especially the seismicity patterns associated with significant earthquakes in East Java. The analysis was carried out by quantifying the spatial and temporal variations of seismicity parameters using b-value and z-value. The used data used are earthquake data from the USGS earthquake catalog during 2002 - 2022 in the East Java region. The analysis was performed using the ZMAP 6.0 program code using the maximum likelihood method. The results of the temporal analysis of b-values show that almost all significant earthquakes in East Java were characterized by a decrease in b-values of up to 0.2. Spatially, the significant earthquake on 10 April 2021 was associated with a low b-value (0.9). Temporal z-value analysis shows that the significant East Java earthquake was preceded by a quiet period (positive z-value) between 4 - 5 years (2003 - 2008 and 2012 - 2016). Spatial analysis of the z-value in the last 7 shows an increase in the z-value before the significant earthquakes occurred on 16 July 2019, 18 March 2020, and 10 April 2021. Therefore, the results of this study indicate that the b-value and z-value potentially can be used as precursors of the significant East Java earthquake. HIGHLIGHTS This is the first time a correlation of b-value and z-value as seismicity parameters has been carried out in East Java, Indonesia, as an attempt to determine the significant earthquake precursors The significant earthquakes in East Java during 2002 - 2022 were marked by a decrease in b-value of up to 0.2 The quiet period in East Java lasted between 4 - 7 years, i.e. 2002 - 2009 and 2012 - 2016, characterized by a positive z-value peaked until +3.5 The z-value at the location of a significant earthquake epicenters show a positive value 2 - 3 years prior to the main quake GRAPHICAL ABSTRACT

A test on methods for Mc estimation and spatial-temporal distribution of b-value in the eastern Tibetan Plateau

Seismic b-value is one of the most important parameters for seismological research and seismic hazards assessment, while the accuracy of the b-value largely depended on the completeness of seismic catalog. This article compares eight methods for estimating the minimum magnitude of completeness (Mc). The results indicate that the modified maximum curvature method (MMAXC), exhibits greater stability and accuracy, closely approximating the standard Mc obtained from the synthetic seismic catalogs. We then calculate the b-value using the instrumental seismic catalog from 2000–2023 in the eastern Tibetan Plateau. The results indicate that the five major earthquakes occur in regions with lower b-value. In addition, the temporal evolution of b-value before and after major earthquakes exhibits a common trend of decreasing before earthquakes, and increasing after earthquakes, which may reflect the stress accumulation and release during earthquakes. Combining the results of maximum shear strain rate and b-value, we identify five regions characterized by low b-value and high shear strain rate, indicating a higher potential seismic hazard in the future.

Acoustic emission characteristic of sandstone and sandstone like material under multi-path loading.

Using spline interpolation to select proportions of similar materials, a comparative analysis of the fracturing behavior of sandstone specimens and similar material specimens was conducted through Brazilian splitting tests under multi-path loading. The study revealed that during stepwise loading, both sandstone and similar materials exhibited memory effects and plastic deformation. However, under constant velocity loading, the relationship between force and displacement in sandstone showed linearity after compaction. Employing MATLAB optimization algorithms for the inversion of acoustic emission event information, the distribution of fracture points, and the evolution of cracks were analyzed. The findings indicated that under stepwise loading, both sandstone and similar materials exhibited banded distribution of peak frequencies, with sandstone concentrated in the mid-low-frequency range and similar materials leaning towards the low-frequency range. The amplitude-frequency characteristics of acoustic emission signals suggested that initially, sandstone produced low-frequency, low-amplitude signals. As cracks developed, these signals gradually transformed into high-frequency, high-amplitude signals, ultimately leading to macroscopic failure. The ringing counts and b-values of sandstone displayed an approximate "W" shape distribution, with a subsequent decrease in b-values during final failure. In contrast, the acoustic emission counts were inversely related to b-values. Similar materials exhibited slightly more acoustic emission counts than sandstone, with relatively lower b-values. The crack development process of both sandstone and similar materials was confirmed through these observations. From the perspective of section initiation and local damage, sandstone and similar materials exhibited similar failure characteristics. The proportions of quartz sand: cement: water = 9:1:0.9 in similar materials demonstrated the most similar characteristics to sandstone in terms of mechanical loading, acoustic emission features, and failure morphology. This suggests that these similar materials can be used as substitutes for sandstone in analogous simulation experiments. The study provides theoretical support for understanding rock fracture mechanisms, offers guidance for the selection and proportioning of similar materials, and holds significance for predicting and controlling rock fracture behavior in engineering applications.

Months-long seismicity transients preceding the 2023 MW 7.8 Kahramanmaraş earthquake, Türkiye

Short term prediction of earthquake magnitude, time, and location is currently not possible. In some cases, however, documented observations have been retrospectively considered as precursory. Here we present seismicity transients starting approx. 8 months before the 2023 MW 7.8 Kahramanmaraş earthquake on the East Anatolian Fault Zone. Seismicity is composed of isolated spatio-temporal clusters within 65 km of future epicentre, displaying non-Poissonian inter-event time statistics, magnitude correlations and low Gutenberg-Richter b-values. Local comparable seismic transients have not been observed, at least since 2014. Close to epicentre and during the weeks prior to its rupture, only scarce seismic activity was observed. The trends of seismic preparatory attributes for this earthquake follow those previously documented in both laboratory stick-slip tests and numerical models of heterogeneous earthquake rupture affecting multiple fault segments. More comprehensive earthquake monitoring together with long-term seismic records may facilitate recognizing earthquake preparation processes from other regional deformation transients.

Statistical Seismic Analysis by b-Value and Occurrence Time of the Latest Earthquakes in Italy

The study reported in this paper concerns the temporal variation in the b-value of the Gutenberg–Richter frequency–magnitude law, applied to the earthquakes that struck Italy from 2009 to 2016 in the geographical areas of L’Aquila, the Emilia Region, and Amatrice–Norcia. Generally, the b-value varies from one region to another dependent on earthquake incidences. Higher values of this parameter are correlated to the occurrence of low-magnitude events spread over a wide geographical area. Conversely, a lower b-value may lead to the prediction of a major earthquake localized along a fault. In addition, it is observed that each seismic event has a different “occurrence time”, which is a key point in the statistical study of earthquakes. In particular, its results are absolutely different for each specific event, and may vary from years to months or even just a few hours. Hence, both short- and long-term precursor phenomena have to be examined. Accordingly, the b-value analysis has to be performed by choosing the best time windows to study the foreshock and aftershock activities.

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.