TY - JOUR A1 - Dahm, Torsten A1 - Fischer, Tomas T1 - Velocity ratio variations in the source region of earthquake swarms in NW Bohemia obtained from arrival time double-differences JF - Geophysical journal international N2 - Crustal earthquake swarms are an expression of intensive cracking and rock damaging over periods of days, weeks or month in a small source region in the crust. They are caused by longer lasting stress changes in the source region. Often, the localized stressing of the crust is associated with fluid or gas migration, possibly in combination with pre-existing zones of weaknesses. However, verifying and quantifying localized fluid movement at depth remains difficult since the area affected is small and geophysical prospecting methods often cannot reach the required resolution. We apply a simple and robust method to estimate the velocity ratio between compressional (P) and shear (S) waves (upsilon(P)/upsilon(S)-ratio) in the source region of an earthquake swarm. The upsilon(P)/upsilon(S)-ratio may be unusual small if the swarm is related to gas in a porous or fractured rock. The method uses arrival time difference between P and S waves observed at surface seismic stations, and the associated double differences between pairs of earthquakes. An advantage is that earthquake locations are not required and the method seems lesser dependent on unknown velocity variations in the crust outside the source region. It is, thus, suited for monitoring purposes. Applications comprise three natural, mid-crustal (8-10 km) earthquake swarms between 1997 and 2008 from the NW-Bohemia swarm region. We resolve a strong temporal decrease of upsilon(P)/upsilon(S) before and during the main activity of the swarm, and a recovery of upsilon(P)/upsilon(S) to background levels at the end of the swarms. The anomalies are interpreted in terms of the Biot-Gassman equations, assuming the presence of oversaturated fluids degassing during the beginning phase of the swarm activity. KW - Tomography KW - Earthquake source observations KW - Volcano seismology Y1 - 2014 U6 - https://doi.org/10.1093/gji/ggt410 SN - 0956-540X SN - 1365-246X VL - 196 IS - 2 SP - 957 EP - 970 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Maghsoudi, Samira A1 - Hainzl, Sebastian A1 - Cesca, Simone A1 - Dahm, Torsten A1 - Kaiser, Diethelm T1 - Identification and characterization of growing large-scale en-echelon fractures in a salt mine JF - Geophysical journal international N2 - The spatiotemporal seismicity of acoustic emission (AE) events recorded in the Morsleben salt mine is investigated. Almost a year after backfilling of the cavities from 2003, microevents are distributed with distinctive stripe shapes above cavities at different depth levels. The physical forces driving the creation of these stripes are still unknown. This study aims to find the active stripes and track fracture developments over time by combining two different temporal and spatial clustering techniques into a single methodological approach. Anomalous seismicity parameters values like sharp b-value changes for two active stripes are good indicators to explain possible stress accumulation at the stripe tips. We identify the formation of two new seismicity stripes and show that the AE activities in active clusters are migrated mostly unidirectional to eastward and upward. This indicates that the growth of underlying macrofractures is controlled by the gradient of extensional stress. Studying size distribution characteristic in terms of frequency-magnitude distribution and b-value in active phase and phase with constant seismicity rate show that deviations from the Gutenberg-Richter power law can be explained by the inclusion of different activity phases: (1) the inactive period before the formation of macrofractures, which is characterized by a deficit of larger events (higher b-values) and (2) the period of fracture growth characterized by the occurrence of larger events (smaller b-values). KW - Earthquake source observations KW - Statistical seismology Y1 - 2014 U6 - https://doi.org/10.1093/gji/ggt443 SN - 0956-540X SN - 1365-246X VL - 196 IS - 2 SP - 1092 EP - 1105 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Cesca, Simone A1 - Sen, Ali Tolga A1 - Dahm, Torsten T1 - Seismicity monitoring by cluster analysis of moment tensors JF - Geophysical journal international N2 - We suggest a new clustering approach to classify focal mechanisms from large moment tensor catalogues, with the purpose of automatically identify families of earthquakes with similar source geometry, recognize the orientation of most active faults, and detect temporal variations of the rupture processes. The approach differs in comparison to waveform similarity methods since clusters are detected even if they occur in large spatial distances. This approach is particularly helpful to analyse large moment tensor catalogues, as in microseismicity applications, where a manual analysis and classification is not feasible. A flexible algorithm is here proposed: it can handle different metrics, norms, and focal mechanism representations. In particular, the method can handle full moment tensor or constrained source model catalogues, for which different metrics are suggested. The method can account for variable uncertainties of different moment tensor components. We verify the method with synthetic catalogues. An application to real data from mining induced seismicity illustrates possible applications of the method and demonstrate the cluster detection and event classification performance with different moment tensor catalogues. Results proof that main earthquake source types occur on spatially separated faults, and that temporal changes in the number and characterization of focal mechanism clusters are detected. We suggest that moment tensor clustering can help assessing time dependent hazard in mines. KW - Persistence KW - memory KW - correlations KW - clustering KW - Earthquake source observations Y1 - 2014 U6 - https://doi.org/10.1093/gji/ggt492 SN - 0956-540X SN - 1365-246X VL - 196 IS - 3 SP - 1813 EP - 1826 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Cesca, Simone A1 - Grigoli, Francesco A1 - Heimann, Sebastian A1 - Gonzalez, Alvaro A1 - Buforn, Elisa A1 - Maghsoudi, Samira A1 - Blanch, Estefania A1 - Dahm, Torsten T1 - The 2013 September-October seismic sequence offshore Spain: a case of seismicity triggered by gas injection? JF - Geophysical journal international N2 - A spatially localized seismic sequence originated few tens of kilometres offshore the Mediterranean coast of Spain, close to the Ebro river delta, starting on 2013 September 5, and lasting at least until 2013 October. The sequence culminated in a maximal moment magnitude M-w 4.3 earthquake, on 2013 October 1. The most relevant seismogenic feature in the area is the Fosa de Amposta fault system, which includes different strands mapped at different distances to the coast, with a general NE-SW orientation, roughly parallel to the coastline. However, no significant known historical seismicity has involved this fault system in the past. The epicentral region is also located near the offshore platform of the Castor project, where gas is conducted through a pipeline from mainland and where it was recently injected in a depleted oil reservoir, at about 2 km depth. We analyse the temporal evolution of the seismic sequence and use full waveform techniques to derive absolute and relative locations, estimate depths and focal mechanisms for the largest events in the sequence (with magnitude mbLg larger than 3), and compare them to a previous event (2012 April 8, mbLg 3.3) taking place in the same region prior to the gas injection. Moment tensor inversion results show that the overall seismicity in this sequence is characterized by oblique mechanisms with a normal fault component, with a 30A degrees low-dip angle plane oriented NNE-SSW and a subvertical plane oriented NW-SE. The combined analysis of hypocentral location and focal mechanisms could indicate that the seismic sequence corresponds to rupture processes along shallow low-dip surfaces, which could have been triggered by the gas injection in the reservoir, and excludes the activation of the Amposta fault, as its known orientation is inconsistent with focal mechanism results. An alternative scenario includes the iterated triggering of a system of steep faults oriented NW-SE, which were identified by prior marine seismics investigations. KW - Earthquake dynamics KW - Earthquake source observations Y1 - 2014 U6 - https://doi.org/10.1093/gji/ggu172 SN - 0956-540X SN - 1365-246X VL - 198 IS - 2 SP - 941 EP - 953 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Grigoli, Francesco A1 - Cesca, Simone A1 - Amoroso, Ortensia A1 - Emolo, Antonio A1 - Zollo, Aldo A1 - Dahm, Torsten T1 - Automated seismic event location by waveform coherence analysis JF - Geophysical journal international N2 - Automated location of seismic events is a very important task in microseismic monitoring operations as well for local and regional seismic monitoring. Since microseismic records are generally characterized by low signal-to-noise ratio, automated location methods are requested to be noise robust and sufficiently accurate. Most of the standard automated location routines are based on the automated picking, identification and association of the first arrivals of P and S waves and on the minimization of the residuals between theoretical and observed arrival times of the considered seismic phases. Although current methods can accurately pick P onsets, the automatic picking of the S onset is still problematic, especially when the P coda overlaps the S wave onset. In this paper, we propose a picking free earthquake location method based on the use of the short-term-average/long-term-average (STA/LTA) traces at different stations as observed data. For the P phases, we use the STA/LTA traces of the vertical energy function, whereas for the S phases, we use the STA/LTA traces of a second characteristic function, which is obtained using the principal component analysis technique. In order to locate the seismic event, we scan the space of possible hypocentral locations and origin times, and stack the STA/LTA traces along the theoretical arrival time surface for both P and S phases. Iterating this procedure on a 3-D grid, we retrieve a multidimensional matrix whose absolute maximum corresponds to the spatial coordinates of the seismic event. A pilot application was performed in the Campania-Lucania region (southern Italy) using a seismic network (Irpinia Seismic Network) with an aperture of about 150 km. We located 196 crustal earthquakes (depth < 20 km) with magnitude range 1.1 < M-L < 2.7. A subset of these locations were compared with accurate manual locations refined by using a double-difference technique. Our results indicate a good agreement with manual locations. Moreover, our method is noise robust and performs better than classical location methods based on the automatic picking of the P and S waves first arrivals. KW - Time-series analysis KW - Inverse theory KW - Earthquake source observations KW - Seismicity and tectonics KW - Early warning Y1 - 2014 U6 - https://doi.org/10.1093/gji/ggt477 SN - 0956-540X SN - 1365-246X VL - 196 IS - 3 SP - 1742 EP - 1753 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Hannemann, Katrin A1 - Krüger, Frank A1 - Dahm, Torsten T1 - Measuring of clock drift rates and static time offsets of ocean bottom stations by means of ambient noise JF - Geophysical journal international N2 - Marine seismology usually relies on temporary deployments of stand alone seismic ocean bottom stations (OBS), which are initialized and synchronized on ship before deployment and re-synchronized and stopped on ship after recovery several months later. In between, the recorder clocks may drift and float at unknown rates. If the clock drifts are large or not linear and cannot be corrected for, seismological applications will be limited to methods not requiring precise common timing. Therefore, for example, array seismological methods, which need very accurate timing between individual stations, would not be applicable for such deployments. We use an OBS test-array of 12 stations and 75 km aperture, deployed for 10 months in the deep sea (4.5-5.5 km) of the mid-eastern Atlantic. The experiment was designed to analyse the potential of broad-band array seismology at the seafloor. After recovery, we identified some stations which either show unusual large clock drifts and/or static time offsets by having a large difference between the internal clock and the GPS-signal (skew). We test the approach of ambient noise cross-correlation to synchronize clocks of a deep water OBS array with km-scale interstation distances. We show that small drift rates and static time offsets can be resolved on vertical components with a standard technique. Larger clock drifts (several seconds per day) can only be accurately recovered if time windows of one input trace are shifted according to the expected drift between a station pair before the cross-correlation. We validate that the drifts extracted from the seismometer data are linear to first order. The same is valid for most of the hydrophones. Moreover, we were able to determine the clock drift at a station where no skew could be measured. Furthermore, we find that instable apparent drift rates at some hydrophones, which are uncorrelated to the seismometer drift recorded at the same digitizer, indicate a malfunction of the hydrophone. KW - Time-series analysis KW - Interferometry KW - Broad-band seismometers Y1 - 2014 U6 - https://doi.org/10.1093/gji/ggt434 SN - 0956-540X SN - 1365-246X VL - 196 IS - 2 SP - 1034 EP - 1042 PB - Oxford Univ. Press CY - Oxford ER -