The large scale extraction of geothermal energy can reduce COsub.2 emissions. For hot dry rocks, the key to successful utilization depends on the efficiency of reservoir reconstruction. The chemical ...and thermal stimulation methods are always used in geothermal reservoir reconstruction except in hydraulic fracturing with high fluid injection pressure, which is believed to reduce the seismic hazard by applying before the high-pressure hydraulic fracturing stimulation. However, at the laboratory scale, there are still very limited experimental studies illustrating the combined effects of chemical and thermal stimulation on the permeability and mechanical properties of granite, which is regarded as the main type of hot dry rock. In this paper, comparative stimulation experiments were carried out, including thermal/cold stimulation, COsub.2 bearing solution hydro-chemical stimulation, combined thermal and COsub.2 bearing fluid stimulation. By means of nuclear magnetic resonance analysis, permeability test and triaxial compression test, the changes of the micro-structure, permeability and mechanical properties of granite under various stimulation conditions were analyzed. The experimental results show that, compared with the single thermal stimulation and COsub.2 bearing fluid hydro-chemical stimulation, the superposition effect of thermal and COsub.2 bearing fluid hydro-chemical stimulation can increase the number of micro-fractures in granite more effectively, thus increasing the permeability, while the elastic modulus and compressive strength decrease. Moreover, the cooling mode on the granite also has a certain influence on the stimulation effect. After water-cooling on the heated granite (300 °C), combined with the COsub.2 bearing fluid stimulation (240 °C, 20 MPa), the permeability of granite is the highest, increasing by 17 times that of the initial state, and the porosity also increases by 144.4%, while the elastic modulus and compressive strength decrease by 14.3% and 18.4%, respectively. This implies that the deterioration of mechanical properties due to the micro-fractures increased by the thermal and chemical stimulation can enhance the fluid conductivity and heat extraction of granite. The methods in this paper can provide a reference for the combined application of thermal and chemical stimulation technology in artificial reservoir reconstruction of hot dry rocks.
SUMMARY
The International Seismological Centre (ISC) is a non‐governmental, non‐profit organization with the primary mission of producing the definitive account of the Earth's seismicity. The ISC ...Bulletin covers some 50 yr (1960–2011) of seismicity. The recent years have seen a dramatic increase both in the number of reported events and especially in the number of reported phases, owing to the ever‐increasing number of stations worldwide. Similar ray paths will produce correlated traveltime prediction errors due to unmodelled heterogeneities in the Earth, resulting in underestimated location uncertainties, and for unfavourable network geometries, location bias. Hence, the denser and more unbalanced the global seismic station coverage becomes, the less defensible is the assumption (that is the observations are independent), which is made by most location algorithms.
To address this challenge we have developed a new location algorithm for the ISC that accounts for correlated error structure, and uses all IASPEI standard phases with a valid ak135 traveltime prediction to obtain more accurate event locations. In this paper we describe the new ISC locator, and present validation tests by relocating the ground truth events in the IASPEI Reference Event List, as well as by relocating the entire ISC Bulletin.
We show that the new ISC location algorithm provides small, but consistent location improvements, considerable improvements in depth determination and significantly more accurate formal uncertainty estimates. We demonstrate that the new algorithm, through the use of later phases and testing for depth resolution, considerably clusters event locations more tightly, thus providing an improved view of the seismicity of the Earth.
SUMMARY
The uneven distribution of earthquakes and stations in seismic tomography leads to slower convergence of nonlinear inversions and spatial bias in inversion results. Including dense regional ...arrays, such as USArray or Hi-Net, in global tomography causes severe convergence and spatial bias problems, against which conventional pre-conditioning schemes are ineffective. To save computational cost and reduce model bias, we propose a new strategy based on a geographical weighting of sources and receivers. Unlike approaches based on ray density or the Voronoi tessellation, this method scales to large full-waveform inversion problems and avoids instabilities at the edges of dense receiver or source clusters. We validate our strategy using a 2-D global waveform inversion test and show that the new weighting scheme leads to a nearly twofold reduction in model error and much faster convergence relative to a conventionally pre-conditioned inversion. We implement this geographical weighting strategy for global adjoint tomography.
SUMMARY
Earthquakes come in clusters formed of mostly aftershock sequences, swarms and occasional foreshock sequences. This clustering is thought to result either from stress transfer among faults, a ...process referred to as cascading, or from transient loading by aseismic slip (pre-slip, afterslip or slow slip events). The ETAS statistical model is often used to quantify the fraction of clustering due to stress transfer and to assess the eventual need for aseismic slip to explain foreshocks or swarms. Another popular model of clustering relies on the earthquake nucleation model derived from experimental rate-and-state friction. According to this model, earthquakes cluster because they are time-advanced by the stress change imparted by the mainshock. This model ignores stress interactions among aftershocks and cannot explain foreshocks or swarms in the absence of transient loading. Here, we analyse foreshock, swarm and aftershock sequences resulting from cascades in a Discrete Fault Network model governed by rate-and-state friction. We show that the model produces realistic swarms, foreshocks and aftershocks. The Omori law, characterizing the temporal decay of aftershocks, emerges in all simulations independently of the assumed initial condition. In our simulations, the Omori law results from the earthquake nucleation process due to rate and state friction and from the heterogeneous stress changes due to the coseismic stress transfers. By contrast, the inverse Omori law, which characterizes the accelerating rate of foreshocks, emerges only in the simulations with a dense enough fault system. A high-density complex fault zone favours fault interactions and the emergence of an accelerating sequence of foreshocks. Seismicity catalogues generated with our discrete fault network model can generally be fitted with the ETAS model but with some material differences. In the discrete fault network simulations, fault interactions are weaker in aftershock sequences because they occur in a broader zone of lower fault density and because of the depletion of critically stressed faults. The productivity of the cascading process is, therefore, significantly higher in foreshocks than in aftershocks if fault zone complexity is high. This effect is not captured by the ETAS model of fault interactions. It follows that a foreshock acceleration stronger than expected from ETAS statistics does not necessarily require aseismic slip preceding the mainshock (pre-slip). It can be a manifestation of a cascading process enhanced by the topological properties of the fault network. Similarly, earthquake swarms might not always imply transient loading by aseismic slip, as they can emerge from stress interactions.
SUMMARY
The Gutenberg–Richter (GR) b-value describes the relative proportion of small to large earthquakes in a scale-free population and is a critical parameter for probabilistic estimation of ...seismic hazard. At low magnitudes, the scale-free behaviour breaks down below the magnitude of completeness mc due to censoring of the data, when the instrumentation used to construct the catalogue is incapable of completely recording all earthquakes in the study region above the background noise. At high magnitudes, it must also break down because natural tectonic and volcanic processes are incapable of an infinite release of energy. This breakdown at large magnitudes is commonly modelled as an exponential roll-off to either the incremental or cumulative GR distribution. This introduces an extra parameter and hence requires relatively more data to justify the additional model complexity. For tectonic seismicity, the estimated b-value is commonly close to unity. In contrast, studies of volcanic and induced seismicity often report significantly higher estimates of the b-value, albeit using relatively small data sets—both in sample size and dynamic (magnitude) range for data above mc. Here, using synthetic data, we show that when we have low dynamic range, it is statistically challenging to test whether the sample is representative of the scale-free GR behaviour or whether it is controlled primarily by the finite size roll-off. We then explore the potential biases that arise when the data quality does not allow this distinction to be made and what the implications are for interpreting studies that have high estimated b-values. We find that systematically higher b-values than those used to generate the synthetic data are regularly obtained when assuming the wrong model and when having a too high mc, resulting in too small catalogues. This is important because it changes our understanding of the accuracy of elevated or variable b-values in catalogues of different dynamic ranges, and quantifies the likely bias in the inferred b-value compared to the underlying true distribution and its associated uncertainty. Finally, we recommend steps to minimize this bias.
SUMMARY
The Gulf of Guinea exhibits a continuous emission of narrow-band and long-period signals (16, 26 and 27 s) on teleseismic records, yet the underlying excitation mechanism remains unclear. ...This study establishes a connection between these tremors and the vibration of thin, decoupled crustal plates at unexplored volcanoes in the gulf. We first formulate the damped plate oscillation equation, by incorporating the vibration of the thin surface crustal plate and magma flow in the subsurface sill. The findings reveal that a fundamental-mode vibration with a period of several dozen seconds can be induced by a crustal plate that is less than 1.0 km thick but extends over tens of kilometres in both length and width, given a subsurface sill depth exceeding 10.0 cm. The thin plate hypothesis also allows for excitation of a few overtone modes, but such waves in higher frequencies diminish over long distances, leaving only the monotonous fundamental-mode vibration at teleseismic stations. The long duration of Guinea tremors at each recurrence is attributed to the presence of low viscosity basaltic magma, which influences the damping factor. Direct wave loads at the shallow gulf serve as the primary vibration source, accounting for seasonal variations and recurring patterns. Sporadic energy bursts may also occur due to large storms. Radiation patterns of Guinea tremors are linked to the geometric structure of the thin plate. Our theoretical estimates of tremor spectra closely align with observed data, confirming the model’s accuracy in capturing reported Guinea tremor characteristics. This study provides valuable insights into the origins of very long-period tremors at continental volcanoes.
SUMMARY
This study examines the feature space of seismic waveforms often used in machine learning applications for seismic event detection and classification problems. Our investigation centres on ...the southern Alaska region, where the seismic record captures diverse seismic activity, notably from the calving of marine-terminating glaciers and tectonic earthquakes along active plate boundaries. While the automated discrimination of earthquakes and glacier quakes is our nominal goal, this data set provides an outstanding opportunity to explore the general feature space of regional seismic phases. That objective has applicability beyond ice quakes and our geographic region of study. We make a noteworthy discovery that features rooted in the spectral content of seismic waveforms consistently outperform statistical and temporal features. Spectral features demonstrate robust performance, exhibiting resilience to class imbalance while being minimally impacted by factors such as epicentral distance and signal-to-noise ratio. We also conduct experiments on the transferability of the model and find that transferability primarily depends on the appearance of the waveforms. Finally, we analyse misclassified events and find examples that are identified incorrectly in the original regional catalogue.