Refine
Has Fulltext
- no (51)
Year of publication
Document Type
- Article (51) (remove)
Language
- English (51)
Is part of the Bibliography
- yes (51)
Keywords
Shallowly situated evaporites in built-up areas are of relevance for urban and cultural development and hydrological regulation. The hazard of sinkholes, subrosion depressions and gypsum karst is often difficult to evaluate and may quickly change with anthropogenic influence. The geophysical exploration of evaporites in metropolitan areas is often not feasible with active industrial techniques. We collect and combine different passive geophysical data as microgravity, ambient vibrations, deformation and hydrological information to study the roof morphology of shallow evaporites beneath Hamburg, Northern Germany. The application of a novel gravity inversion technique leads to a 3-D depth model of the salt diapir under study. We compare the gravity-based depth model to pseudo-depths from H/V measurements and depth estimates from small-scale seismological array data. While the general range and trend of the diapir roof is consistent, a few anomalous regions are identified where H/V pseudo-depths indicate shallower structures not observed in gravity or array data. These are interpreted by shallow residual caprock floaters and zones of increased porosity. The shallow salt structure clearly correlates with a relative subsidence in the order of 2 mm yr(-1). The combined interpretation of roof morphology, yearly subsidence rates, chemical analyses of groundwater and of hydraulic head in aquifers indicates that the salt diapir beneath Hamburg is subject to significant ongoing dissolution that may possibly affect subrosion depressions, sinkhole distribution and land usage. The combined analysis of passive geophysical data may be exemplary for the study of shallow evaporites beneath other urban areas.
Reliable estimations of magnitude of completeness (M-c) are essential for a correct interpretation of seismic catalogues. The spatial distribution of M-c may be strongly variable and difficult to assess in mining environments, owing to the presence of galleries, cavities, fractured regions, porous media and different mineralogical bodies, as well as in consequence of inhomogeneous spatial distribution of the seismicity. We apply a 3-D modification of the probabilistic magnitude of completeness (PMC) method, which relies on the analysis of network detection capabilities. In our approach, the probability to detect an event depends on its magnitude, source receiver Euclidian distance and source receiver direction. The suggested method is proposed for study of the spatial distribution of the magnitude of completeness in a mining environment and here is applied to a 2-months acoustic emission (AE) data set recorded at the Morsleben salt mine, Germany. The dense seismic network and the large data set, which includes more than one million events, enable a detailed testing of the method. This method is proposed specifically for strongly heterogeneous media. Besides, it can also be used for specific network installations, with sensors with a sensitivity, dependent on the direction of the incoming wave (e.g. some piezoelectric sensors). In absence of strong heterogeneities, the standards PMC approach should be used. We show that the PMC estimations in mines strongly depend on the source receiver direction, and cannot be correctly accounted using a standard PMC approach. However, results can be improved, when adopting the proposed 3-D modification of the PMC method. Our analysis of one central horizontal and vertical section yields a magnitude of completeness of about M-c approximate to 1 (AE magnitude) at the centre of the network, which increases up to M-c approximate to 4 at further distances outside the network; the best detection performance is estimated for a NNE-SSE elongated region, which corresponds to the strike direction of the low-attenuating salt body. Our approach provides us with small-scale details about the capability of sensors to detect an earthquake, which can be linked to the presence of heterogeneities in specific directions. Reduced detection performance in presence of strong structural heterogeneities (cavities) is confirmed by synthetic waveform modelling in heterogeneous media.
On 2012 May 20 and 29, two damaging earthquakes with magnitudes M-w 6.1 and 5.9, respectively, struck the Emilia-Romagna region in the sedimentary Po Plain, Northern Italy, causing 26 fatalities, significant damage to historical buildings and substantial impact to the economy of the region. The earthquake sequence included four more aftershocks with M-w, >= 5.0, all at shallow depths (about 7-9 km), with similar WNW-ESE striking reverse mechanism. The timeline of the sequence suggests significant static stress interaction between the largest events. We perform here a detailed source inversion, first adopting a point source approximation and considering pure double couple and full moment tensor source models. We compare different extended source inversion approaches for the two largest events, and find that the rupture occurred in both cases along a subhorizontal plane, dipping towards SSW Directivity is well detected for the May 20 main shock, indicating that the rupture propagated unilaterally towards SE. Based on the focal mechanism solution, we further estimate the co-seismic static stress change induced by the May 20 event. By using the rate-and-state model and a Poissonian earthquake occurrence, we infer that the second largest event of May 29 was induced with a probability in the range 0.2-0.4. This suggests that the segment of fault was already prone to rupture. Finally, we estimate peak ground accelerations for the two main events as occurred separately or simultaneously. For the scenario involving hypothetical rupture areas of both main events, we estimate M-w = 6.3 and an increase of ground acceleration by 50 per cent. The approach we propose may help to quantify rapidly which regions are invested by a significant increase of the hazard, bearing the potential for large aftershocks or even a second main shock.
Seismicity induced by coal mining in the Ruhr region, Germany, has been monitored continuously over the last 25 yr. In 2006, a dense temporary network (HAMNET) was deployed to locally monitor seismicity induced by longwall mining close to the town of Hamm. Between 2006 July and 2007 July, more than 7000 events with magnitudes M-L from -1.7 to 2.0 were detected. The spatiotemporal distribution of seismicity shows high correlation with the mining activity. In order to monitor rupture processes, we set up an automated source inversion routine and successfully perform double couple and full moment tensor (MT) inversions for more than 1000 events with magnitudes above M-L -0.5. The source inversion is based on a full waveform approach, both in the frequency and in the time domain, providing information about the centroid location, focal mechanism, scalar moment and full MT. Inversion results indicate a strong dominance of normal faulting focal mechanisms, with a steeper plane and a subhorizontal one. Fault planes are oriented parallel to the mining stopes. We classify the focal mechanisms based on their orientation and observe different frequency-magnitude distributions for families of events with different focal mechanisms; the overall frequency-magnitude distribution is not fitting the Gutenberg-Richter relation. Full MTs indicate that non-negligible opening tensile components accompanied normal faulting source mechanisms. Finally, extended source models are investigated for largest events. Results suggest that the rupture processes mostly occurred along the subvertical planes.
Salt diapirs are common features of sedimentary basins. If close to the surface, they can bear a significant hazard due to possible dissolution sinkholes, karst formation and collapse dolines or their influence on ground water chemistry. We investigate the potential of ambient vibration techniques to map the 3-D roof morphology of shallow salt diapirs. Horizontal-to-vertical (H/V) spectral peaks are derived at more than 900 positions above a shallow diapir beneath the city area of Hamburg, Germany, and are used to infer the depth of the first strong impedance contrast. In addition, 15 small-scale array measurements are conducted at different positions in order to compute frequency-dependent phase velocities of Rayleigh waves between 0.5 and 25 Hz. The dispersion curves are inverted together with the H/V peak frequency to obtain shear-wave velocity profiles. Additionally, we compare the morphology derived from H/V and array measurements to borehole lithology and a gravity-based 3-D model of the salt diapir. Both methods give consistent results in agreement with major features indicated by the independent data. An important result is that H/V and array measurements are better suited to identify weathered gypsum caprocks or gypsum floaters, while gravity-derived models better sample the interface between sediments and homogeneous salt. We further investigate qualitatively the influence of the 3-D subsurface topography of the salt diapir on the validity of local 1-D inversion results from ambient vibration dispersion curve inversion.
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.
Identification and characterization of growing large-scale en-echelon fractures in a salt mine
(2014)
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).
Aseismic transient driving the swarm-like seismic sequence in the Pollino range, Southern Italy
(2015)
Tectonic earthquake swarms challenge our understanding of earthquake processes since it is difficult to link observations to the underlying physical mechanisms and to assess the hazard they pose. Transient forcing is thought to initiate and drive the spatio-temporal release of energy during swarms. The nature of the transient forcing may vary across sequences and range from aseismic creeping or transient slip to diffusion of pore pressure pulses to fluid redistribution and migration within the seismogenic crust. Distinguishing between such forcing mechanisms may be critical to reduce epistemic uncertainties in the assessment of hazard due to seismic swarms, because it can provide information on the frequency-magnitude distribution of the earthquakes (often deviating from the assumed Gutenberg-Richter relation) and on the expected source parameters influencing the ground motion (for example the stress drop). Here we study the ongoing Pollino range (Southern Italy) seismic swarm, a long-lasting seismic sequence with more than five thousand events recorded and located since October 2010. The two largest shocks (magnitude M-w = 4.2 and M-w = 5.1) are among the largest earthquakes ever recorded in an area which represents a seismic gap in the Italian historical earthquake catalogue. We investigate the geometrical, mechanical and statistical characteristics of the largest earthquakes and of the entire swarm. We calculate the focal mechanisms of the M-l > 3 events in the sequence and the transfer of Coulomb stress on nearby known faults and analyse the statistics of the earthquake catalogue. We find that only 25 per cent of the earthquakes in the sequence can be explained as aftershocks, and the remaining 75 per cent may be attributed to a transient forcing. The b-values change in time throughout the sequence, with low b-values correlated with the period of highest rate of activity and with the occurrence of the largest shock. In the light of recent studies on the palaeoseismic and historical activity in the Pollino area, we identify two scenarios consistent with the observations and our analysis: This and past seismic swarms may have been 'passive' features, with small fault patches failing on largely locked faults, or may have been accompanied by an 'active', largely aseismic, release of a large portion of the accumulated tectonic strain. Those scenarios have very different implications for the seismic hazard of the area.
In this study, we analyze acoustic emission (AE) data recorded at the Morsleben salt mine, Germany, to assess the catalog completeness, which plays an important role in any seismicity analysis. We introduce the new concept of a magnitude completeness interval consisting of a maximum magnitude of completeness (M-c(max)) in addition to the well-known minimum magnitude of completeness. This is required to describe the completeness of the catalog, both for the smallest events (for which the detection performance may be low) and for the largest ones (which may be missed because of sensors saturation). We suggest a method to compute the maximum magnitude of completeness and calculate it for a spatial grid based on (1) the prior estimation of saturation magnitude at each sensor, (2) the correction of the detection probability function at each sensor, including a drop in the detection performance when it saturates, and (3) the combination of detection probabilities of all sensors to obtain the network detection performance. The method is tested using about 130,000 AE events recorded in a period of five weeks, with sources confined within a small depth interval, and an example of the spatial distribution of M-c(max) is derived. The comparison between the spatial distribution of M-c(max) and of the maximum possible magnitude (M-max), which is here derived using a recently introduced Bayesian approach, indicates that M-max exceeds M-c(max) in some parts of the mine. This suggests that some large and important events may be missed in the catalog, which could lead to a bias in the hazard evaluation.