Refine
Year of publication
Document Type
- Article (61)
- Other (9)
- Postprint (2)
- Monograph/Edited Volume (1)
Keywords
- Body waves (7)
- Seismologie (6)
- Seismology (6)
- Earthquake source observations (5)
- Erdbeben (5)
- Seismicity and tectonics (5)
- Theoretical seismology (4)
- Array Seismology (3)
- Arrayseismologie (3)
- Crustal structure (3)
We compute seismic velocity profiles by a combined inversion of surface-wave phase-velocity dispersion curves together with the full spectrum of the microtremor horizontal-to-vertical (H/V) spectral ratio at two sediment-covered sites in Germany. The sediment deposits are approximately 100 m thick at the first test site and approximately 400 m thick at the second test site. We have used an extended physical model based on the diffuse wavefield assumption for the interpretation of the observed microtremor H/V spectral ratio. The extension includes the interpretation of the microtremor H/V spectral ratio observed at depth (in boreholes). This full-wavefield approach accounts for the energy contribution from the body and surface waves, and thus it allows for inverting the properties of the shallow subsurface. We have obtained the multimode phase velocity dispersion curves from an independent study, and a description of the extracted branches and their interpretation was developed. The inversion results indicate that the combined approach using seismic ambient noise and actively generated surface-wave data will improve the accuracy of the reconstructed near-surface velocity model, a key step in microzonation, geotechnical engineering, seismic statics corrections, and reservoir imaging.
We focus on the relation between seismic and total postseismic afterslip following the Maule M-w 8.8 earthquake on 2010 February 27 in central Chile. First, we calculate the cumulative slip released by aftershock seismicity. We do this by summing up the aftershock regions and slip estimated from scaling relations. Comparing the cumulative seismic slip with afterslip modelswe showthat seismic slip of individual aftershocks exceeds locally the inverted afterslip model from geodetic constraints. As the afterslip model implicitly contains the displacements from the aftershocks, this reflects the tendency of afterslip models to smear out the actual slip pattern. However, it also suggests that locally slip for a number of the larger aftershocks exceeds the aseismic slip in spite of the fact that the total equivalent moment of the afterslip exceeds the cumulative moment of aftershocks by a large factor. This effect, seen weakly for the Maule 2010 and also for the Tohoku 2011 earthquake, can be explained by taking into account the uncertainties of the seismicity and afterslip models. In spite of uncertainties, the hypocentral region of the Nias 2005 earthquake is suggested to release a large fraction of moment almost purely seismically. Therefore, these aftershocks are not driven solely by the afterslip but instead their slip areas have probably been stressed by interseismic loading and the mainshock rupture. In a second step, we divide the megathrust of the Maule 2010 rupture into discrete cells and count the number of aftershocks that occur within 50 km of the centre of each cell as a function of time. We then compare this number to a time-dependent afterslip model by defining the 'afterslip to aftershock ratio' (ASAR) for each cell as the slope of the best fitting line when the afterslip at time t is plotted against aftershock count. Although we find a linear relation between afterslip and aftershocks for most cells, there is significant variability in ASAR in both the downdip and along-strike directions of the megathrust. We compare the spatial distribution of ASAR with the spatial distribution of seismic coupling, coseismic slip and Bouguer gravity anomaly, and in each case we find no significant correlation.
We use recently deployed seismological arrays in Africa to sample a 2D cross section through the mantle down to the core-mantle boundary (CMB). By making use of travel-time residuals of S, ScS, and SKS phases, a new shear-velocity model of the African low-velocity zone (ALVZ) is derived. Our model suggests between 1.2% shear-velocity reduction at the top and 5% at the bottom with respect to 1D reference models. The average reduction over the whole low-velocity zone (LVZ) amounts to 2% in the presented model and is therefore about twice as strong as values found in global tomographic models. The top of the LVZ reaches up to 1200-km depth, and its lateral extent at the CMB is about 35 degrees. We propose the existence of a gap of 300 km, splitting the structure into two blocks. Our results are based on remarkable differences in SK(K) S travel-time residuals over a few degrees distance. The complexity of the structure could provide a key to an improved understanding of the deep-mantle LVZ dynamics and composition by comparison to geodynamic models. The gap in the model might suggest that the 2D cross section is cutting through a 3D indentation in the boundary of the ALVZ but may also be interpreted as a sign of two individual plumes, rather than one large homogeneous upwelling.
Continent-ocean-transition across a trans-tensional margin segment: off Bear Island, Barents Sea
(2011)
P>A 410 km long Ocean Bottom Seismometer profile spanning from the Bear Island, Barents Sea to oceanic crust formed along the Mohns Ridge has been modelled by use of ray-tracing with regard to observed P-waves. The northeastern part of the model represents typical continental crust, thinned from ca. 30 km thickness beneath the Bear Island to ca. 13 km within the Continent-Ocean-Transition. Between the Hornsund FZ and the Kn circle divide legga Fault, a 3-4 km thick sedimentary basin, dominantly of Permian/Carboniferous age, is modelled beneath the ca. 1.5 km thick layer of volcanics (Vestbakken Volcanic Province). The P-wave velocity in the 3-4 km thick lowermost continental crust is significantly higher than normal (ca. 7.5 km s-1). We interpret this layer as a mixture of mafic intrusions and continental crystalline blocks, dominantly related to the Paleocene-Early Eocene rifting event. The crystalline portion of the crust within the south-western part of the COT consists of a ca. 30 km wide and ca. 6 km thick high-velocity (7.3 km s-1) body. We interpret the body as a ridge of serpentinized peridotites. The magmatic portion of the ocean crust accreted along the Knipovich Ridge from continental break-up at ca. 35 Ma until ca. 20 Ma is 3-5 km thicker than normal. We interpret the increased magmatism as a passive response to the bending of this southernmost part of the Knipovich Ridge. The thickness of the magmatic portion of the crust formed along the Mohns Ridge at ca. 20 Ma decreases to ca. 3 km, which is normal for ultra slow spreading ridges.
Earthquakes occurring close to hydrocarbon fields under production are often under critical view of being induced or triggered. However, clear and testable rules to discriminate the different events have rarely been developed and tested. The unresolved scientific problem may lead to lengthy public disputes with unpredictable impact on the local acceptance of the exploitation and field operations. We propose a quantitative approach to discriminate induced, triggered, and natural earthquakes, which is based on testable input parameters. Maxima of occurrence probabilities are compared for the cases under question, and a single probability of being triggered or induced is reported. The uncertainties of earthquake location and other input parameters are considered in terms of the integration over probability density functions. The probability that events have been human triggered/induced is derived from the modeling of Coulomb stress changes and a rate and state-dependent seismicity model. In our case a 3-D boundary element method has been adapted for the nuclei of strain approach to estimate the stress changes outside the reservoir, which are related to pore pressure changes in the field formation. The predicted rate of natural earthquakes is either derived from the background seismicity or, in case of rare events, from an estimate of the tectonic stress rate. Instrumentally derived seismological information on the event location, source mechanism, and the size of the rupture plane is of advantage for the method. If the rupture plane has been estimated, the discrimination between induced or only triggered events is theoretically possible if probability functions are convolved with a rupture fault filter. We apply the approach to three recent main shock events: (1) the M-w 4.3 Ekofisk 2001, North Sea, earthquake close to the Ekofisk oil field; (2) the M-w 4.4 Rotenburg 2004, Northern Germany, earthquake in the vicinity of the Sohlingen gas field; and (3) the M-w 6.1 Emilia 2012, Northern Italy, earthquake in the vicinity of a hydrocarbon reservoir. The three test cases cover the complete range of possible causes: clearly human induced, not even human triggered, and a third case in between both extremes.
Earthquake source arrays
(2020)
A collection of earthquake sources recorded at a single station, under specific conditions, are considered as a source array (SA), that is interpreted as if earthquake sources originate at the station location and are recorded at the source location. Then, array processing methods, that is array beamforming, are applicable to analyse the recorded signals. A possible application is to use source array multiple event techniques to locate and characterize near-source scatterers and structural interfaces. In this work the aim is to facilitate the use of earthquake source arrays by presenting an automatic search algorithm to configure the source array elements. We developed a procedure to search for an optimal source array element distribution given an earthquake catalogue including accurate origin time and hypocentre locations. The objective function of the optimization process can be flexibly defined for each application to ensure the prerequisites (criteria) of making a source array. We formulated four quantitative criteria as subfunctions and used the weighted sum technique to combine them in one single scalar function. The criteria are: (1) to control the accuracy of the slowness vector estimation using the time domain beamforming method, (2) to measure the waveform coherency of the array elements, (3) to select events with lower location error and (4) to select traces with high energy of specific phases, that is, sp- or ps-phases. The proposed procedure is verified using synthetic data as well as real examples for the Vogtland region in Northwest Bohemia. We discussed the possible application of the optimized source arrays to identify the location of scatterers in the velocity model by presenting a synthetic test and an example using real waveforms.
The complex system of strike-slip and thrust faults in the Alborz Mountains, Northern Iran, are not well understood yet. Mainly structural and geomorphic data are available so far. As a more extensive base for seismotectonic studies and seismic hazard analysis we plan to do a comprehensive seismic moment tensor study also from smaller magnitudes (M < 4.5) by developing a new algorithm. Here, we present first preliminary results.
We develop an amplitude spectral ratio method for event couples from clustered earthquakes to estimate seismic wave attenuation (Q-1) in the source volume. The method allows to study attenuation within the source region of earthquake swarms or aftershocks at depth, independent of wave path and attenuation between source region and surface station. We exploit the high-frequency slope of phase spectra using multitaper spectral estimates. The method is tested using simulated full wave-field seismograms affected by recorded noise and finite source rupture. The synthetic tests verify the approach and show that solutions are independent of focal mechanisms but also show that seismic noise may broaden the scatter of results. We apply the event couple spectral ratio method to northwest Bohemia, Czech Republic, a region characterized by the persistent occurrence of earthquake swarms in a confined source region at mid-crustal depth. Our method indicates a strong anomaly of high attenuation in the source region of the swarm with an averaged attenuation factor of Qp < 100. The application to S phases fails due to scattered P-phase energy interfering with S phases. The Qp anomaly supports the common hypothesis of highly fractured and fluid saturated rocks in the source region of the swarms in northwest Bohemia. However, high temperatures in a small volume around the swarms cannot be excluded to explain our observations.
We develop an amplitude spectral ratio method for event couples from clustered earthquakes to estimate seismic wave attenuation (Q-1) in the source volume. The method allows to study attenuation within the source region of earthquake swarms or aftershocks at depth, independent of wave path and attenuation between source region and surface station. We exploit the high-frequency slope of phase spectra using multitaper spectral estimates. The method is tested using simulated full wave-field seismograms affected by recorded noise and finite source rupture. The synthetic tests verify the approach and show that solutions are independent of focal mechanisms but also show that seismic noise may broaden the scatter of results. We apply the event couple spectral ratio method to northwest Bohemia, Czech Republic, a region characterized by the persistent occurrence of earthquake swarms in a confined source region at mid-crustal depth. Our method indicates a strong anomaly of high attenuation in the source region of the swarm with an averaged attenuation factor of Qp < 100. The application to S phases fails due to scattered P-phase energy interfering with S phases. The Qp anomaly supports the common hypothesis of highly fractured and fluid saturated rocks in the source region of the swarms in northwest Bohemia. However, high temperatures in a small volume around the swarms cannot be excluded to explain our observations.
Knowledge of the quality factor of near-surface materials is of fundamental interest in various applications. Attenuation can be very strong close to the surface and thus needs to be properly assessed. In recent years, several researchers have studied the retrieval of attenuation coefficients from the cross correlation of ambient seismic noise. Yet, the determination of exact amplitude information from noise-correlation functions is, in contrast to the extraction of traveltimes, not trivial. Most of the studies estimated attenuation coefficients on the regional scale and within the microseism band. In this paper, we investigate the possibility to derive attenuation coefficients from seismic noise at much shallower depths and higher frequencies (> 1 Hz). The Euroseistest area in northern Greece offers ideal conditions to study quality factor retrieval from ambient noise for different rock types. Correlations are computed between the stations of a small scale array experiment (station spacings < 2 km) that was carried out in the Euroseistest area in 2011. We employ the correlation of the coda of the correlation (C-3) method instead of simple cross correlations to mitigate the effect of uneven noise source distributions on the correlation amplitude. Transient removal and temporal flattening are applied instead of 1-bit normalization in order to retain relative amplitudes. The C-3 method leads to improved correlation results (higher signal-to-noise ratio and improved time symmetry) compared to simple cross correlations. The C-3 functions are rotated from the ZNE to the ZRT system and we focus on Love wave arrivals on the transverse component and on Love wave quality factors Q(L). The analysis is performed for selected stations being either situated on soft soil or on weathered rock. Phase slowness is extracted using a slant-stack method. Attenuation parameters are inferred by inspecting the relative amplitude decay of Love waves with increasing interstation distance. We observe that the attenuation coefficient gamma and Q(L) can be reliably extracted for stations situated on soft soil whereas the derivation of attenuation parameters is more problematic for stations that are located on weathered rock. The results are in acceptable conformance with theoretical Love wave attenuation curves that were computed using 1-D shear wave velocity and quality factor profiles from the Euroseistest area.
Fluid intelligence is the ability to think flexibly and to understand abstract relations. People with high fluid intelligence (hi-fluIQ) perform better in analogical reasoning tasks than people with average fluid intelligence (ave-fluIQ). Although previous neuroimaging studies reported involvement of parietal and frontal brain regions in geometric analogical reasoning (which is a prototypical task for fluid intelligence), however, neuroimaging findings on geometric analogical reasoning in hi-fluIQ are sparse. Furthermore, evidence on the relation between brain activation and intelligence while solving cognitive tasks is contradictory. The present study was designed to elucidate the cerebral correlates of geometric analogical reasoning in a sample of hi-fluIQ and ave-fluIQ high school students. We employed a geometric analogical reasoning task with graded levels of task difficulty and confirmed the involvement of the parieto-frontal network in solving this task. In addition to characterizing the brain regions involved in geometric analogical reasoning in hi-fluIQ and ave-fluIQ, we found that blood oxygenation level dependency (BOLD) signal changes were greater for hi-fluIQ than for ave-fluIQ in parietal brain regions. However, ave-fluIQ showed greater BOLD signal changes in the anterior cingulate cortex and medial frontal gyrus than hi-fluIQ. Thus, we showed that a similar network of brain regions is involved in geometric analogical reasoning in both groups. Interestingly, the relation between brain activation and intelligence is not mono-directional, but rather, it is specific for each brain region. The negative brain activation-intelligence relationship in frontal brain regions in hi-fluIQ goes along with a better behavioral performance and reflects a lower demand for executive monitoring compared to ave-fluIQ individuals. In conclusion, our data indicate that flexibly modulating the extent of regional cerebral activity is characteristic for fluid intelligence.
The correct estimation of site-specific attenuation is crucial for the assessment of seismic hazard. Downhole instruments provide in this context valuable information to constrain attenuation directly from data. In this study, we apply an interferometric approach to this problem by deconvolving seismic motions recorded at depth with those recorded at the surface. In doing so, incident and surface-reflected waves can be separated. We apply this technique not only to earthquake data but also to recordings of ambient vibrations. We compute the transfer function between incident and surface-reflected waves in order to infer frequency-dependent quality factors for S waves. The method is applied to a 87m deep borehole sensor and a colocated surface instrument situated at a hard-rock site in West Bohemia/Vogtland, Germany. We show that the described method provides comparable attenuation estimates using either earthquake data or ambient noise for frequencies between 5 and 15 Hz. Moreover, a single hour of noise recordings seems to be sufficient to yield stable deconvolution traces and quality factors, thus, offering a fast and easy way to derive attenuation estimates from borehole recordings even in low- to mid-seismicity regions.
The H/V spectral ratio has emerged as a single station method within the seismic ambient noise analysis field by its capability to quickly estimate the frequency of resonance at a site and through inversion the average profile information. Although it is easy to compute from experimental data, its counter theoretical part is not obvious when building a forward model which can help in reconstructing the derived H/V spectrum. This has led to the simplified assumption that the noise wavefield is mainly composed of Rayleigh waves and the derived H/V often used without further correction. Furthermore, only the right (and left) flank around the H/V peak frequency is considered in the inversion for the subsurface 1-D shear wave velocity profile. A new theoretical approach for the interpretation of the H/V spectral ratio has been presented by Sanchez-Sesmaet al. In this paper, the fundamental idea behind their theory is presented as it applies to receivers at depth. A smooth H/V(z, f) spectral curve on a broad frequency range is obtained by considering a fine integration step which is in turn time consuming. We show that for practical purposes and in the context of inversion, this can be considerably optimized by using a coarse integration step combined with the smoothing of the corresponding directional energy density (DED) spectrum. Further analysis shows that the obtained H/V(z, f) spectrum computed by the mean of the imaginary part of Green's function method could also be recovered using the reflectivity method for a medium well illuminated by seismic sources. Inversion of synthetic H/V(z, f) spectral curve is performed for a single layer over a half space. The striking results allow to potentially use the new theory as a forward computation of the H/V(z, f) to fully invert the experimental H/V spectral ratio at the corresponding depth for the shear velocity profile (Vs) and additionally the compressional velocity profile (Vp) using receivers both at the surface and in depth. We use seismic ambient noise data in the frequency range of 0.2-50 Hz recorded at two selected sites in Germany where borehole information is also available. The obtained 1-D Vs and Vp profiles are correlated with geological log information. Results from shallow geophysical experiment are also used for comparison.
We track a bilateral rupture propagation lasting similar to 160 s, with its dominant branch rupturing northeastwards at about 3 kms(-1). The area of maximum energy emission is offset from the maximum coseismic slip but matches the zone where most plate interface aftershocks occur. Along dip, energy is preferentially released from two disconnected interface belts, and a distinct jump from the shallower belt to the deeper one is visible after about 20 s from the onset. However, both belts keep on being active until the end of the rupture. These belts approximately match the position of the interface aftershocks, which are split into two clusters of events at different depths, thus suggesting the existence of a repeated transition from stick-slip to creeping frictional regime.
Microseismic data from observatories in Europe, which have been continuously recorded since about 100 years, contain information on the wave-climate in the North Atlantic. They can potentially be used as additional constraints in high-resolution temporal and spatial reconstructions of the storminess and oceanic waveheights in the past. To resolve spatial patterns data from observatories in different regions are needed. While previous recent studies analyzed only few observatory archives and relatively short time ranges, this is a first attempt to process the whole available data archive from different observatories. We correct and compare smoothed microseismic data from different stations and discuss their correlation and possible use for studies of storminess variability. Microseismic amplitudes at four seismic stations in northern Europe show amplitude peaks in 1920 and 1925, a slow decline in amplitudes till the middle of the 1930's followed by a steady increase of amplitudes till about 1990. From 1990 on microseismic amplitudes decrease. We find a good correlation between the average surface wind velocity in the North Atlantic and microseismic amplitudes at inland stations far away from the coast. Coastal stations are more influenced by local swell and are thus potentially useful to recover regional changes in wind and ocean wavefields with time. The study demonstrates that the analysis of microseismic has the potential to assess climate changes during the last 100 years
The majority of original seismograms recorded at the very beginning of instrumental seismology (the early 1900s) did not survive till present. However, a number of books, bulletins, and catalogs were published including the seismogram reproductions of some, particularly interesting earthquakes. In case these reproductions contain the time and amplitude scales, they can be successfully analyzed the same way as the original records. Information about the Atushi (Kashgar) earthquake, which occurred on August 22, 1902, is very limited. We could not find any original seismograms for this earthquake, but 12 seismograms from 6 seismic stations were printed as example records in different books. These data in combination with macroseismic observations and different bulletins information published for this earthquake were used to determine the source parameters of the earthquake. The earthquake epicenter was relocated at 39.87A degrees N and 76.42A degrees E with the hypocenter depth of about 18 km. We could further determine magnitudes m (B) = 7.7 +/- 0.3, M (S) = 7.8 +/- 0.4, M (W) = 7.7 +/- 0.3 and the focal mechanism of the earthquake with strike/dip/rake - 260A degrees +/- 20/30A degrees +/- 10/90A degrees +/- 10. This study confirms that the earthquake likely had a smaller magnitude than previously reported (M8.3). The focal mechanism indicates dominant thrust faulting, which is in a good agreement with presumably responsible Tuotegongbaizi-Aerpaleike northward dipping thrust fault kinematic, described in previous studies.
The new in situ geodynamic laboratory established in the framework of the ICDP Eger project aims to develop the most modern, comprehensive, multiparameter laboratory at depth for studying earthquake swarms, crustal fluid flow, mantle-derived CO2 and helium degassing, and processes of the deep biosphere. In order to reach a new level of high-frequency, near-source and multiparameter observation of earthquake swarms and related phenomena, such a laboratory comprises a set of shallow boreholes with high-frequency 3-D seismic arrays as well as modern continuous real-time fluid monitoring at depth and the study of the deep biosphere.
This laboratory is located in the western part of the Eger Rift at the border of the Czech Republic and Germany (in the West Bohemia–Vogtland geodynamic region) and comprises a set of five boreholes around the seismoactive zone. To date, all monitoring boreholes have been drilled. This includes the seismic monitoring boreholes S1, S2 and S3 in the crystalline units north and east of the major Nový Kostel seismogenic zone, borehole F3 in the Hartoušov mofette field and borehole S4 in the newly discovered Bažina maar near Libá. Supplementary borehole P1 is being prepared in the Neualbenreuth maar for paleoclimate and biological research. At each of these sites, a borehole broadband seismometer will be installed, and sites S1, S2 and S3 will also host a 3-D seismic array composed of a vertical geophone chain and surface seismic array. Seismic instrumenting has been completed in the S1 borehole and is in preparation in the remaining four monitoring boreholes. The continuous fluid monitoring site of Hartoušov includes three boreholes, F1, F2 and F3, and a pilot monitoring phase is underway. The laboratory also enables one to analyze microbial activity at CO2 mofettes and maar structures in the context of changes in habitats. The drillings into the maar volcanoes contribute to a better understanding of the Quaternary paleoclimate and volcanic activity.
The new in situ geodynamic laboratory established in the framework of the ICDP Eger project aims to develop the most modern, comprehensive, multiparameter laboratory at depth for studying earthquake swarms, crustal fluid flow, mantle-derived CO2 and helium degassing, and processes of the deep biosphere. In order to reach a new level of high-frequency, near-source and multiparameter observation of earthquake swarms and related phenomena, such a laboratory comprises a set of shallow boreholes with high-frequency 3-D seismic arrays as well as modern continuous real-time fluid monitoring at depth and the study of the deep biosphere.
This laboratory is located in the western part of the Eger Rift at the border of the Czech Republic and Germany (in the West Bohemia-Vogtland geodynamic region) and comprises a set of five boreholes around the seismoactive zone. To date, all monitoring boreholes have been drilled. This includes the seismic monitoring boreholes S1, S2 and S3 in the crystalline units north and east of the major Novy Kostel seismogenic zone, borehole F3 in the Hartousov mofette field and borehole S4 in the newly discovered Bazina maar near Liba. Supplementary borehole P1 is being prepared in the Neualbenreuth maar for paleoclimate and biological research. At each of these sites, a borehole broadband seismometer will be installed, and sites S1, S2 and S3 will also host a 3-D seismic array composed of a vertical geophone chain and surface seismic array. Seismic instrumenting has been completed in the S1 borehole and is in preparation in the remaining four monitoring boreholes. The continuous fluid monitoring site of Hartousov includes three boreholes, F1, F2 and F3, and a pilot monitoring phase is underway. The laboratory also enables one to analyze microbial activity at CO2 mofettes and maar structures in the context of changes in habitats. The drillings into the maar volcanoes contribute to a better understanding of the Quaternary paleoclimate and volcanic activity.