Frequency‐Dependent Moment Tensors of Induced Microearthquakes

Analysis of 984 induced microearthquakes from The Geysers geothermal reservoir in California reveals that the retrieved moment tensors depend on the frequency band of the inverted waveforms. The observed dependence is more significant for the percentages of the double‐couple, compensated linear vector dipole, and isotropic (ISO) components than for the focal mechanisms. The average root‐mean‐square of the moment tensors obtained in different frequency bands is correlated with spectra of ambient noise. The percentages of double‐couple and ISO components tend to decrease and increase with the upper cutoff frequency (fu), respectively. This suggests that shear rupture radiates energy preferentially in a lower frequency band and tensile rupture in a higher frequency band. Events displaying a strong increase of the ISO with fu are confined within the same depth interval as the injection points. This might be related to the strong thermoelastic effects in the vicinity of injection points that promote opening of small cracks adjacent to the main fractures

Seismotectonic setting at the North Anatolian Fault Zone after the 1999 Mw=7.4 Izmit earthquake based on high-resolution aftershock locations

International audienceThe most recent devastating earthquakes that occurred along the North Anatolian Fault Zone (NAFZ) in northwestern Turkey were the 1999 Izmit (Mw=7.4) and Düzce (Mw=7.1) events. In this study we present a catalog of Izmit aftershock hypocenters that was deduced from a network covering the entire 140 km long rupture of the mainshock. 7348 events with a location accuracy better than 5 km are analysed. Aftershocks were observed along the entire ruptured segment along a 20 km wide band of activity. Events are clustered in distinct regions and dominantly occur at 5 to 15 km depth. The eastern termination of the Izmit rupture is characterized by a sharp and steeply dipping boundary exactly where the Düzce mainshock initiated 87 days after the Izmit event. Relocation of the events using double-difference technology results in 4696 high-resolution hypocenters that allow resolving the internal structure of the seismically active areas with a resolution of 300 m (horizontal) and 400m (vertical). Below the Akyazi Plain, representing a small pull-apart structure at a triple junction of the NAFZ, we identify planes of activity that can be correlated with nodal planes of EW extensional normal faulting aftershocks. Along the easternmost Karadere-Düzce segment we identify the down-dip extension of the Karadere fault that hosted about 1 m of right-lateral coseismic slip. At the easternmost rupture we correlate a cloud-type distribution of seismic activity with the largest aftershocks in this area, a subevent of the Izmit mainshock and the Düzce mainshock that all have an almost identical focal mechanism. This part of the NAFZ is interpreted as a classical example of a seismic barrier along the fault

Gas and seismicity within the Istanbul seismic gap

Understanding micro-seismicity is a critical question for earthquake hazard assessment. Since the devastating earthquakes of Izmit and Duzce in 1999, the seismicity along the submerged section of North Anatolian Fault within the Sea of Marmara (comprising the “Istanbul seismic gap”) has been extensively studied in order to infer its mechanical behaviour (creeping vs locked). So far, the seismicity has been interpreted only in terms of being tectonic- driven, although the Main Marmara Fault (MMF) is known to strike across multiple hydrocarbon gas sources. Here, we show that a large number of the aftershocks that followed the M 5.1 earthquake of July, 25th 2011 in the western Sea of Marmara, occurred within a zone of gas overpressuring in the 1.5–5 km depth range, from where pressurized gas is expected to migrate along the MMF, up to the surface sediment layers. Hence, gas-related processes should also be considered for a complete interpretation of the micro- seismicity (~M < 3) within the Istanbul offshore domain

Gas and seismicity within the Istanbul seismic gap

Understanding micro-seismicity is a critical question for earthquake hazard assessment. Since the devastating earthquakes of Izmit and Duzce in 1999, the seismicity along the submerged section of North Anatolian Fault within the Sea of Marmara (comprising the “Istanbul seismic gap”) has been extensively studied in order to infer its mechanical behaviour (creeping vs locked). So far, the seismicity has been interpreted only in terms of being tectonic-driven, although the Main Marmara Fault (MMF) is known to strike across multiple hydrocarbon gas sources. Here, we show that a large number of the aftershocks that followed the M 5.1 earthquake of July, 25th 2011 in the western Sea of Marmara, occurred within a zone of gas overpressuring in the 1.5–5 km depth range, from where pressurized gas is expected to migrate along the MMF, up to the surface sediment layers. Hence, gas-related processes should also be considered for a complete interpretation of the micro-seismicity (~M < 3) within the Istanbul offshore domain

Metre-scale damage zone characterization using S-coda waves from active ultrasonic transmission measurements in the STIMTEC project, URL Reiche Zeche, Germany

Studies of controlled hydraulic stimulation experiments with active and passive seismic monitoring conducted in Underground Research Laboratories (URLs) benefit from specific knowledge of hydraulic parameters, close by microseismic monitoring revealing structural details of the rock mass, and detailed evolution of seismicity in response to injection operations. Microseismic monitoring is commonly used to characterize a stimulated reservoir volume, for example, in terms of damage evolution of the rock mass. Since seismic attenuation is affected by damage of the rock volume, active seismic sources covering sizes from the centimetre to decimetre scale may help us to investigate space–time varying attenuation properties in a reservoir. This may allow us to monitor damage evolution of the stimulated rock volume in more detail, also since active seismic sources produce stronger signals leading to a broader frequency range that can be analysed compared to passive seismic signals. Within the STIMTEC project in the URL Reiche Zeche (URL-RZ) in Freiberg (Germany), more than 300 active Ultrasonic Transmission (UT) measurements were performed before and after hydraulic stimulations in two boreholes in the targeted rock volume, an anisotropic metamorphic gneiss. The signal-frequency content ranges between 1 and 60 kHz. Assuming scattering attenuation to dominate over intrinsic attenuation, we here apply the single isotropic scattering model. S-coda waves of 88 spatially representative UT measurements are used to estimate the coda quality factor (QC). We obtain stable QC estimates for centre frequencies of octave-width frequency bands between 3 and 21 kHz. We group neighbouring UT measurements to stabilize the observations and form eight UT groups in total, covering different depth intervals in three boreholes and four different time periods to investigate scattering attenuation changes in a spatiotemporal manner. Our final mean QC (⁠QC¯¯¯¯¯¯¯⁠) estimates show characteristic frequency-dependence as observed at the field scale in geological reservoirs. We find temporal variations of QC are strongly connected to hydraulic stimulation, and these variations are more significant than those resolved from velocity changes. QC¯¯¯¯¯¯¯ estimates at frequencies above 15 kHz indicate healing of injection-induced small-scale fractures during a two-months post-stimulation phase. Larger fractures, mostly sampled by lower frequencies (<15 kHz), seem to be more persistent with time (over 15 months). We observe spatial differences of QC¯¯¯¯¯¯¯ values near the mine galleries (driftway and vein drift) and relate these observations to different extents and characteristics of the galleries’ excavation damage zones. Our results further support previous assumptions based on borehole televiewer logs and mapped structures of an existing fault with larger damage zone that crosses the stimulated rock volume NW-SE between the galleries. We conclude that the coda analysis of active UT measurements complements established imaging methods used during experiments in URLs. In particular, coda analysis is a powerful tool for the detection of damage zones and for monitoring local fracture networks with immediate application for imaging georeservoirs considered for exploitation or underground storage of gases and liquids

Sea Level Changes Affect Seismicity Rates in a Hydrothermal System Near Istanbul

Small stress changes such as those from sea level fluctuations can be large enough to trigger earthquakes. If small and large earthquakes initiate similarly, high-resolution catalogs with low detection thresholds are best suited to illuminate such processes. Below the Sea of Marmara section of the North Anatolian Fault, a segment of urn:x-wiley:00948276:media:grl65397:grl65397-math-0001150 km is late in its seismic cycle. We generated high-resolution seismicity catalogs for a hydrothermal region in the eastern Sea of Marmara employing AI-based and template matching techniques to investigate the link between sea level fluctuations and seismicity over 6 months. All high resolution catalogs show that local seismicity rates are larger during time periods shortly after local minima of sea level, when it is already rising. Local strainmeters indicate that seismicity is promoted when the ratio of differential to areal strain is the largest. The strain changes from sea level variations, on the order of 30–300 nstrain, are sufficient to promote seismicity

Crustal Thickness Variation Across the Sea of Marmara Region, NW Turkey: A Reflection of Modern and Ancient Tectonic Processes

The Marmara region in Turkey is an important geological setting, both from a tectonic and a seismic hazard/risk perspective. We present a new map of crustal thickness variation across this complex region to better understand the interplay of past and present tectonic processes that have formed present‐day structure. Maps of crustal thickness are created using Ps converted phases and receiver function (RF) analysis of earthquakes recorded at all publicly available seismic stations and stations in the national monitoring network (run by AFAD Disaster and Emergency Management Authority Turkey). RFs are converted from time to depth using a local 3‐D full‐waveform tomographic model and are combined in multiphase common conversion point stacks. Direct P to S converted arrivals and associated multiples are mapped to produce continuous maps of the Moho discontinuity. Results show Moho depths ranging from 26–41 km with a regional trend of westward thinning reflecting the effects of the extensional regime in western Anatolia and the neighboring Aegean Sea. The thinnest crust is observed beneath the western end of the Sea of Marmara, attributed to transtensional basin opening. A distinct region of increased crustal thickness bounded by the West Black Sea Fault in the west, and the northern strand of the North Anatolian Fault in the south, defines the ancient crustal terrane of the Istanbul Zone. Isostatic arguments indicate that the thickened crust and lower elevation in the Istanbul Zone require it to be underlain by thicker lithosphere, a conclusion that is consistent with its hypothesized origin near the Odessa shelf

Anatomy of the high-frequency ambient seismic wave field at the TCDP borehole.

International audienceThe Taiwan Chelungpu-fault Drilling Project (TCDP) installed a vertical seismic array between 950 and 1270 m depth in an active thrust fault environment. In this paper we analyze continuous noise records of the TCDP array between 1 and 16 Hz. We apply multiple array processing and noise correlation techniques to study the noise source process, properties of the propagation medium, and the ambient seismic wave field. Diurnal amplitude and slowness patterns suggest that noise is generated by cultural activity. The vicinity of the recording site to the excitation region, indicated by a narrow azimuthal distribution of propagation directions, leads to a predominant ballistic propagation regime. This is evident from the compatibility of the data with an incident plane wave model, polarized direct arrivals of noise correlation functions, and the asymmetric arrival shape. Evidence for contributions from scattering comes from equilibrated earthquake coda energy ratios, the frequency dependent randomization of propagation directions, and the existence of correlation coda waves. We conclude that the ballistic and scattered propagation regime coexist, where the first regime dominates the records, but the second is weaker yet not negligible. Consequently, the wave field is not equipartitioned. Correlation signal-to-noise ratios indicate a frequency dependent noise intensity. Iterations of the correlation procedure enhance the signature of the scattered regime. Discrepancies between phase velocities estimated from correlation functions and in-situ measurements are associated with the array geometry and its relative orientation to the predominant energy flux. The stability of correlation functions suggests their applicability in future monitoring efforts