@article{KielingRoesslerKrueger2011, author = {Kieling, Katrin and R{\"o}ßler, Dirk and Kr{\"u}ger, Frank}, title = {Receiver function study in northern Sumatra and the Malaysian peninsula}, series = {Journal of seismology}, volume = {15}, journal = {Journal of seismology}, number = {2}, publisher = {Springer}, address = {Dordrecht}, issn = {1383-4649}, doi = {10.1007/s10950-010-9222-7}, pages = {235 -- 259}, year = {2011}, abstract = {In this receiver function study, we investigate the structure of the crust beneath six seismic broadband stations close to the Sunda Arc formed by subduction of the Indo-Australian under the Sunda plate. We apply three different methods to analyse receiver functions at single stations. A recently developed algorithm determines absolute shear-wave velocities from observed frequency-dependent apparent incidence angles of P waves. Using waveform inversion of receiver functions and a modified Zhu and Kanamori algorithm, properties of discontinuities such as depth, velocity contrast, and sharpness are determined. The combination of the methods leads to robust results. The approach is validated by synthetic tests. Stations located on Malaysia show high-shear-wave velocities (V (S)) near the surface in the range of 3.4-3.6 km s (-aEuro parts per thousand 1) attributed to crystalline rocks and 3.6-4.0 km s (-aEuro parts per thousand 1) in the lower crust. Upper and lower crust are clearly separated, the Moho is found at normal depths of 30-34 km where it forms a sharp discontinuity at station KUM or a gradient at stations IPM and KOM. For stations close to the subduction zone (BSI, GSI and PSI) complexity within the crust is high. Near the surface low V (S) of 2.6-2.9 km s (-aEuro parts per thousand 1) indicate sediment layers. High V (S) of 4.2 km s (-aEuro parts per thousand 1) are found at depth greater than 6 and 2 km at BSI and PSI, respectively. There, the Moho is located at 37 and 40 km depth. At station GSI, situated closest to the trench, the subducting slab is imaged as a north-east dipping structure separated from the sediment layer by a 10 km wide gradient in V (S) between 10 and 20 km depth. Within the subducting slab V (S) a parts per thousand aEuro parts per thousand 4.7 km s (-aEuro parts per thousand 1). At station BSI, the subducting slab is found at depth between 90 and 110 km dipping 20A degrees +/- 8A degrees in approximately N 60A degrees E. A velocity increase in similar depth is indicated at station PSI, however no evidence for a dipping layer is found.}, language = {en} } @article{BlaserKruegerOhrnbergeretal.2010, author = {Blaser, Lilian and Kr{\"u}ger, Frank and Ohrnberger, Matthias and Scherbaum, Frank}, title = {Scaling relations of earthquake source parameter estimates with special focus on subduction environment}, issn = {0037-1106}, doi = {10.1785/0120100111}, year = {2010}, abstract = {Earthquake rupture length and width estimates are in demand in many seismological applications. Earthquake magnitude estimates are often available, whereas the geometrical extensions of the rupture fault mostly are lacking. Therefore, scaling relations are needed to derive length and width from magnitude. Most frequently used are the relationships of Wells and Coppersmith (1994) derived on the basis of a large dataset including all slip types with the exception of thrust faulting events in subduction environments. However, there are many applications dealing with earthquakes in subduction zones because of their high seismic and tsunamigenic potential. There are no well-established scaling relations for moment magnitude and length/width for subduction events. Within this study, we compiled a large database of source parameter estimates of 283 earthquakes. All focal mechanisms are represented, but special focus is set on (large) subduction zone events, in particular. Scaling relations were fitted with linear least-square as well as orthogonal regression and analyzed regarding the difference between continental and subduction zone/oceanic relationships. Additionally, the effect of technical progress in earthquake parameter estimation on scaling relations was tested as well as the influence of different fault mechanisms. For a given moment magnitude we found shorter but wider rupture areas of thrust events compared to Wells and Coppersmith (1994). The thrust event relationships for pure continental and pure subduction zone rupture areas were found to be almost identical. The scaling relations differ significantly for slip types. The exclusion of events prior to 1964 when the worldwide standard seismic network was established resulted in a remarkable effect on strike-slip scaling relations: the data do not show any saturation of rupture width of strike- slip earthquakes. Generally, rupture area seems to scale with mean slip independent of magnitude. The aspect ratio L/W, however, depends on moment and differs for each slip type.}, language = {en} } @article{BlaserOhrnbergerKruegeretal.2012, author = {Blaser, Lilian and Ohrnberger, Matthias and Kr{\"u}ger, Frank and Scherbaum, Frank}, title = {Probabilistic tsunami threat assessment of 10 recent earthquakes offshore Sumatra}, series = {Geophysical journal international}, volume = {188}, journal = {Geophysical journal international}, number = {3}, publisher = {Wiley-Blackwell}, address = {Malden}, issn = {0956-540X}, doi = {10.1111/j.1365-246X.2011.05324.x}, pages = {1273 -- 1284}, year = {2012}, abstract = {Tsunami early warning (TEW) is a challenging task as a decision has to be made within few minutes on the basis of incomplete and error-prone data. Deterministic warning systems have difficulties in integrating and quantifying the intrinsic uncertainties. In contrast, probabilistic approaches provide a framework that handles uncertainties in a natural way. Recently, we have proposed a method using Bayesian networks (BNs) that takes into account the uncertainties of seismic source parameter estimates in TEW. In this follow-up study, the method is applied to 10 recent large earthquakes offshore Sumatra and tested for its performance. We have evaluated both the general model performance given the best knowledge we have today about the source parameters of the 10 events and the corresponding response on seismic source information evaluated in real-time. We find that the resulting site-specific warning level probabilities represent well the available tsunami wave measurements and observations. Difficulties occur in the real-time tsunami assessment if the moment magnitude estimate is severely over- or underestimated. In general, the probabilistic analysis reveals a considerably large range of uncertainties in the near-field TEW. By quantifying the uncertainties the BN analysis provides important additional information to a decision maker in a warning centre to deal with the complexity in TEW and to reason under uncertainty.}, language = {en} } @article{GassnerThomasKruegeretal.2015, author = {Gassner, Alexandra and Thomas, Christine and Kr{\"u}ger, Frank and Weber, Michael H.}, title = {Probing the core-mantle boundary beneath Europe and Western Eurasia: A detailed study using PcP}, series = {Physics of the earth and planetary interiors}, volume = {246}, journal = {Physics of the earth and planetary interiors}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0031-9201}, doi = {10.1016/j.pepi.2015.06.007}, pages = {9 -- 24}, year = {2015}, abstract = {We use PcP (the core reflected P phase) recordings of deep earthquakes and nuclear explosions from the Grafenberg (Germany) and NORSAR (Norway) arrays to investigate the core-mantle boundary region beneath Europe and western Eurasia. We find evidence for a previously unknown ultra-low velocity zone 600 km south-east of Moscow, located at the edge of a middle-size low shear- velocity region imaged in seismic tomography that is located beneath the Volga river region. The observed amplitude variations of PcP can be modelled by velocity reductions of P and S-waves of -5\% and -15\%, respectively, with a density increase of +15\%. Travel time delays of pre-and postcursors are indicating a thickness of about 13 km for this ultra-low velocity region (ULVZ). However, our modelling also reveals highly ambiguous amplitude variations of PcP and a reflection off the top of the anomaly for various ULVZs and topography models. Accordingly, large velocity contrasts of up to -10\% in V-P and -20\% in Vs cannot be excluded. In general, the whole Volga river region shows a complex pattern of PcP amplitudes caused most likely by CMB undulations. Further PcP probes beneath Paris, Kiev and northern Italy indicate likely normal CMB conditions, whereas the samples below Finland and the Hungary-Slovakia border yield strongly amplified PcP signals suggesting strong CMB topography effects. We evaluate the amplitude behaviour of PcP as a function of distance and several ULVZ models using the 1D reflectivity and the 2D Gauss beam method. The influence of the velocity and density perturbations is analysed as well as the anomaly thickness, the dominant period of the source wavelet and interface topographies. Strong variation of the PcP amplitude are obtained as a function of distance and of the impedance contrast. We also consider two types of topographies: undulations atop the CMB in the presence of flat ULVZs and vice versa. Where a broad range of CMB topography dimensions lead to large PcP amplitude variations, only large ULVZ undulations generate significant amplitude scattering. Consequently, this indicates that topography effects of anomalies may mask the true medium parameters as well as the ULVZ thickness. Moreover, there might be a possibility of misinterpreting the precursor as PcP, in particular for thin ULVZs. (C) 2015 Elsevier B.V. All rights reserved.}, language = {en} } @article{LontsiJoseSanchezSesmaCamilloMolinaVillegasetal.2015, author = {Lontsi, Agostiny Marrios and Jose Sanchez-Sesma, Francisco and Camillo Molina-Villegas, Juan and Ohrnberger, Matthias and Kr{\"u}ger, Frank}, title = {Full microtremor H/V(z,f) inversion for shallow subsurface characterization}, series = {Geophysical journal international}, volume = {202}, journal = {Geophysical journal international}, number = {1}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggv132}, pages = {298 -- 312}, year = {2015}, abstract = {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.}, language = {en} } @article{KulikovaKrueger2015, author = {Kulikova, Galina and Kr{\"u}ger, Frank}, title = {Source process of the 1911 M8.0 Chon-Kemin earthquake: investigation results by analogue seismic records}, series = {Geophysical journal international}, volume = {201}, journal = {Geophysical journal international}, number = {3}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggv091}, pages = {1891 -- 1911}, year = {2015}, abstract = {Several destructive earthquakes have occurred in Tien-Shan region at the beginning of 20th century. However, the detailed seismological characteristics, especially source parameters of those earthquakes are still poorly investigated. The Chon-Kemin earthquake is the strongest instrumentally recorded earthquake in the Tien-Shan region. This earthquake has produced an approximately 200 km long system of surface ruptures along Kemin-Chilik fault zone and killed about similar to 400 people. Several studies presented the different information on the earthquake epicentre location and magnitude, and two different focal mechanisms were also published. The reason for the limited knowledge of the source parameters for the Chon-Kemin earthquake is the complexity of old analogue records processing, digitization and analysis. In this study the data from 23 seismic stations worldwide were collected and digitized. The earthquake epicentre was relocated to 42.996NA degrees and 77.367EA degrees, the hypocentre depth is estimated between 10 and 20 km. The magnitude was recalculated to m(B) 8.05, M-s 7.94 and M-w 8.02. The focal mechanism, determined from amplitude ratios comparison of the observed and synthetic seismograms, was: str = 264A degrees, dip = 52A degrees, rake = 98A degrees. The apparent source time duration was between similar to 45 and similar to 70 s, the maximum slip occurred 25 s after the beginning of the rupture. Two subevents were clearly detected from the waveforms with the scalar moment ratio between them of about 1/3, the third subevent was also detected with less certainty. Taking into account surface rupture information, the fault geometry model with three patches was proposed. Based on scaling relations we conclude that the total rupture length was between similar to 260 and 300 km and a maximum rupture width could reach similar to 70 km.}, language = {en} } @article{MatosSilveiraMatiasetal.2015, author = {Matos, Catarina and Silveira, Graca and Matias, Luis and Caldeira, Rita and Ribeiro, M. Luisa and Dias, Nuno A. and Kr{\"u}ger, Frank and Bento dos Santos, Telmo}, title = {Upper crustal structure of Madeira Island revealed from ambient noise tomography}, series = {Journal of volcanology and geothermal research}, volume = {298}, journal = {Journal of volcanology and geothermal research}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0377-0273}, doi = {10.1016/j.jvolgeores.2015.03.017}, pages = {136 -- 145}, year = {2015}, abstract = {We present the first image of the Madeira upper crustal structure, using ambient seismic noise tomography. 16 months of ambient noise, recorded in a dense network of 26 seismometers deployed across Madeira, allowed reconstructing Rayleigh wave Green's functions between receivers. Dispersion analysis was performed in the short period band from 1.0 to 4.0 s. Group velocity measurements were regionalized to obtain 20 tomographic images, with a lateral resolution of 2.0 km in central Madeira. Afterwards, the dispersion curves, extracted from each cell of the 2D group velocity maps, were inverted as a function of depth to obtain a 3D shear wave velocity model of the upper crust, from the surface to a depth of 2.0 km. The obtained 3D velocity model reveals features throughout the island that correlates well with surface geology and island evolution. (C) 2015 Elsevier B.V. All rights reserved.}, language = {en} } @article{WeberWicksLeStunffetal.2015, author = {Weber, Michael H. and Wicks, Charles and Le Stunff, Yves and Romanowicz, Barbara and Kr{\"u}ger, Frank}, title = {Seismic evidence for a steeply dipping reflector-stagnant slab in the mantle transition zone}, series = {Geophysical journal international}, volume = {200}, journal = {Geophysical journal international}, number = {2}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggu438}, pages = {1235 -- 1251}, year = {2015}, abstract = {Studies of seismic tomography have been highly successful at imaging the deep structure of subduction zones. In a study complementary to these tomographic studies, we use array seismology and reflected waves to image a stagnant slab in the mantle transition zone. Using P and S (SH) waves we find a steeply dipping reflector centred at ca. 400 km depth and ca. 550 km west of the present Mariana subduction zone (at 20N, 140E). The discovery of this anomaly in tomography and independently in array seismology (this paper) helps in understanding the evolution of the Mariana margin. The reflector/stagnant slab may be the remains of the hypothetical North New Guinea Plate, which was theorized to have subducted ca. 50 Ma.}, language = {en} } @misc{ZaliReinKruegeretal.2023, author = {Zali, Zahra and Rein, Teresa and Kr{\"u}ger, Frank and Ohrnberger, Matthias and Scherbaum, Frank}, title = {Ocean bottom seismometer (OBS) noise reduction from horizontal and vertical components using harmonic-percussive separation algorithms}, series = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {1320}, issn = {1866-8372}, doi = {10.25932/publishup-58882}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-588828}, pages = {15}, year = {2023}, abstract = {Records from ocean bottom seismometers (OBSs) are highly contaminated by noise, which is much stronger compared to data from most land stations, especially on the horizontal components. As a consequence, the high energy of the oceanic noise at frequencies below 1 Hz considerably complicates the analysis of the teleseismic earthquake signals recorded by OBSs. Previous studies suggested different approaches to remove low-frequency noises from OBS recordings but mainly focused on the vertical component. The records of horizontal components, which are crucial for the application of many methods in passive seismological analysis of body and surface waves, could not be much improved in the teleseismic frequency band. Here we introduce a noise reduction method, which is derived from the harmonic-percussive separation algorithms used in Zali et al. (2021), in order to separate long-lasting narrowband signals from broadband transients in the OBS signal. This leads to significant noise reduction of OBS records on both the vertical and horizontal components and increases the earthquake signal-to-noise ratio (SNR) without distortion of the broadband earthquake waveforms. This is demonstrated through tests with synthetic data. Both SNR and cross-correlation coefficients showed significant improvements for different realistic noise realizations. The application of denoised signals in surface wave analysis and receiver functions is discussed through tests with synthetic and real data.}, language = {en} } @article{ZaliReinKruegeretal.2023, author = {Zali, Zahra and Rein, Teresa and Kr{\"u}ger, Frank and Ohrnberger, Matthias and Scherbaum, Frank}, title = {Ocean bottom seismometer (OBS) noise reduction from horizontal and vertical components using harmonic-percussive separation algorithms}, series = {Solid earth}, volume = {14}, journal = {Solid earth}, number = {2}, publisher = {Coepernicus Publ.}, address = {G{\"o}ttingen}, issn = {1869-9529}, doi = {10.5194/se-14-181-2023}, pages = {181 -- 195}, year = {2023}, abstract = {Records from ocean bottom seismometers (OBSs) are highly contaminated by noise, which is much stronger compared to data from most land stations, especially on the horizontal components. As a consequence, the high energy of the oceanic noise at frequencies below 1 Hz considerably complicates the analysis of the teleseismic earthquake signals recorded by OBSs. Previous studies suggested different approaches to remove low-frequency noises from OBS recordings but mainly focused on the vertical component. The records of horizontal components, which are crucial for the application of many methods in passive seismological analysis of body and surface waves, could not be much improved in the teleseismic frequency band. Here we introduce a noise reduction method, which is derived from the harmonic-percussive separation algorithms used in Zali et al. (2021), in order to separate long-lasting narrowband signals from broadband transients in the OBS signal. This leads to significant noise reduction of OBS records on both the vertical and horizontal components and increases the earthquake signal-to-noise ratio (SNR) without distortion of the broadband earthquake waveforms. This is demonstrated through tests with synthetic data. Both SNR and cross-correlation coefficients showed significant improvements for different realistic noise realizations. The application of denoised signals in surface wave analysis and receiver functions is discussed through tests with synthetic and real data.}, language = {en} } @article{LontsiOhrnbergerKrueger2016, author = {Lontsi, Agostiny Marrios and Ohrnberger, Matthias and Kr{\"u}ger, Frank}, title = {Shear wave velocity profile estimation by integrated analysis of active and passive seismic data from small aperture arrays}, series = {Journal of applied geophysics}, volume = {130}, journal = {Journal of applied geophysics}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0926-9851}, doi = {10.1016/j.jappgeo.2016.03.034}, pages = {37 -- 52}, year = {2016}, abstract = {We present an integrated approach for deriving the 1D shear wave velocity (Vs) information at few tens to hundreds of meters down to the first strong impedance contrast in typical sedimentary environments. We use multiple small aperture seismic arrays in 1D and 2D configuration to record active and passive seismic surface wave data at two selected geotechnical sites in Germany (Horstwalde \& Lobnitz). Standard methods for data processing include the Multichannel Analysis of Surface Waves (MASW) method that exploits the high frequency content in the active data and the sliding window frequency-wavenumber (f-k) as well as the spatial autocorrelation (SPAC) methods that exploit the low frequency content in passive seismic data. Applied individually, each of the passive methods might be influenced by any source directivity in the noise wavefield. The advantages of active shot data (known source location) and passive microtremor (low frequency content) recording may be combined using a correlation based approach applied to the passive data in the so called Interferometric Multichannel Analysis of Surface Waves (IMASW). In this study, we apply those methods to jointly determine and interpret the dispersion characteristics of surface waves recorded at Horstwalde and Lobnitz. The reliability of the dispersion curves is controlled by applying strict limits on the interpretable range of wavelengths in the analysis and further avoiding potentially biased phase velocity estimates from the passive f-k method by comparing to those derived from the SPatial AutoCorrelation method (SPAC). From our investigation at these two sites, the joint analysis as proposed allows mode extraction in a wide frequency range (similar to 0.6-35 Hz at Horstwalde and similar to 1.5-25 Hz at Lobnitz) and consequently improves the Vs profile inversion. To obtain the shear wave velocity profiles, we make use of a global inversion approach based on the neighborhood algorithm to invert the interpreted branches of the dispersion curves. Within the uncertainty given by the apparent spread of forward models we find that besides a well defined sediment velocity range also a reasonable minimum estimate of bedrock depth and bedrock velocity can be achieved. The Vs estimate for the best model in Horstwalde ranges from similar to 190 m/s at the surface up to similar to 390 m/s in the bottom of the soft sediment column. The bedrock starts earliest around 200 m depth and bedrock velocities are higher than 1000 m/s. In Lobnitz, we observe slightly lower velocities for the sediments (similar to 165-375 m/s for the best model) and a minimum thickness of 75 m. (C) 2016 Elsevier B.V. All rights reserved.}, language = {en} } @article{HaendelOhrnbergerKrueger2016, author = {H{\"a}ndel, Annabel and Ohrnberger, Matthias and Kr{\"u}ger, Frank}, title = {Extracting near-surface Q(L) between 1-4 Hz from higher-order noise correlations in the Euroseistest area, Greece}, series = {Geophysical journal international}, volume = {207}, journal = {Geophysical journal international}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggw295}, pages = {655 -- 666}, year = {2016}, abstract = {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.}, language = {en} } @article{LontsiOhrnbergerKruegeretal.2016, author = {Lontsi, Agostiny Marrios and Ohrnberger, Matthias and Kr{\"u}ger, Frank and S{\´a}nchez-Sesma, Francisco Jos{\´e}}, title = {Combining surface-wave phase-velocity dispersion curves and full microtremor horizontal-to-vertical spectral ratio for subsurface sedimentary site characterization}, series = {Interpretation : a journal of subsurface characterization}, volume = {4}, journal = {Interpretation : a journal of subsurface characterization}, publisher = {Society of Exploration Geophysicists}, address = {Tulsa}, issn = {2324-8858}, doi = {10.1190/INT-2016-0021.1}, pages = {SQ41 -- SQ49}, year = {2016}, abstract = {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.}, language = {en} } @article{KulikovaSchurrKruegeretal.2016, author = {Kulikova, Galina and Schurr, Bernd and Kr{\"u}ger, Frank and Brzoska, Elisabeth and Heimann, Sebastian}, title = {Source parameters of the Sarez-Pamir earthquake of 1911 February 18}, series = {Geophysical journal international}, volume = {205}, journal = {Geophysical journal international}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggw069}, pages = {1086 -- 1098}, year = {2016}, abstract = {The Ms \&\#8764; 7.7 Sarez-Pamir earthquake of 1911 February 18 is the largest instrumentally recorded earthquake in the Pamir region. It triggered one of the largest landslides of the past century, building a giant natural dam and forming Lake Sarez. As for many strong earthquakes from that time, information about source parameters of the Sarez-Pamir earthquake is limited due to the sparse observations. Here, we present the analysis of analogue seismic records of the Sarez-Pamir earthquake. We have collected, scanned and digitized 26 seismic records from 13 stations worldwide to relocate the epicentre and determine the event's depth (\&\#8764;26 km) and magnitude (mB7.3 and Ms7.7). The unusually good quality of the digitized waveforms allowed their modelling, revealing an NE-striking sinistral strike-slip focal mechanism in accordance with regional tectonics. The shallow depth and magnitude (Mw7.3) of the earthquake were confirmed. Additionally, we investigated the possible contribution of the landslide to the waveforms and present an alternative source model assuming the landslide and earthquake occurred in close sequence.}, language = {en} } @article{LangeBedfordMorenoetal.2014, author = {Lange, Dietrich and Bedford, J. R. and Moreno, M. and Tilmann, F. and B{\´a}ez, Juan Carlos and Bevis, M. and Kr{\"u}ger, Frank}, title = {Comparison of postseismic afterslip models with aftershock seismicity for three subduction-zone earthquakes: Nias 2005, Maule 2010 and Tohoku 2011}, series = {Geophysical journal international}, volume = {199}, journal = {Geophysical journal international}, number = {2}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggu292}, pages = {784 -- 799}, year = {2014}, abstract = {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.}, language = {en} } @article{KnapmeyerEndrunKrueger2014, author = {Knapmeyer-Endrun, Brigitte and Kr{\"u}ger, Frank}, title = {Moho depth across the Trans-European Suture Zone from P- and S-receiver functions}, series = {Geophysical journal international}, volume = {197}, journal = {Geophysical journal international}, number = {2}, publisher = {Oxford Univ. Press}, address = {Oxford}, organization = {PASSEQ Working Grp}, issn = {0956-540X}, doi = {10.1093/gji/ggu035}, pages = {1048 -- 1075}, year = {2014}, abstract = {The Mohorovicic discontinuity, Moho for short, which marks the boundary between crust and mantle, is the main first-order structure within the lithosphere. Geodynamics and tectonic evolution determine its depth level and properties. Here, we present a map of the Moho in central Europe across the Teisseyre-Tornquist Zone, a region for which a number of previous studies are available. Our results are based on homogeneous and consistent processing of P-and S-receiver functions for the largest passive seismological data set in this region yet, consisting of more than 40 000 receiver functions from almost 500 station. Besides, we also provide new results for the crustal vP/vS ratio for the whole area. Our results are in good agreement with previous, more localized receiver function studies, as well as with the interpretation of seismic profiles, while at the same time resolving a higher level of detail than previous maps covering the area, for example regarding the Eifel Plume region, Rhine Graben and northern Alps. The close correspondence with the seismic data regarding crustal structure also increases confidence in use of the data in crustal corrections and the imaging of deeper structure, for which no independent seismic information is available. In addition to the pronounced, stepwise transition from crustal thicknesses of 30 km in Phanerozoic Europe to more than 45 beneath the East European Craton, we can distinguish other terrane boundaries based on Moho depth as well as average crustal v(P)/v(S) ratio and Moho phase amplitudes. The terranes with distinct crustal properties span a wide range of ages, from Palaeoproterozoic in Lithuania to Cenozoic in the Alps, reflecting the complex tectonic history of Europe. Crustal thickness and properties in the study area are also markedly influenced by tectonic overprinting, for example the formation of the Central European Basin System, and the European Cenozoic Rift System. In the areas affected by Cenozoic rifting and volcanism, thinning of the crust corresponds to lithospheric updoming reported in recent surface wave and S-receiver function studies, as expected for thermally induced deformation. The same correlation applies for crustal thickening, not only across the Trans-European Suture Zone, but also within the southern part of the Bohemian Massif.}, language = {en} } @article{CzubaGradMjeldeetal.2011, author = {Czuba, Wojciech and Grad, Marek and Mjelde, Rolf and Guterch, Aleksander and Libak, Audun and Kr{\"u}ger, Frank and Murai, Yoshio and Schweitzer, Johannes}, title = {Continent-ocean-transition across a trans-tensional margin segment: off Bear Island, Barents Sea}, series = {Geophysical journal international}, volume = {184}, journal = {Geophysical journal international}, number = {2}, publisher = {Oxford Univ. Press}, address = {Oxford}, organization = {IPY Project Grp}, issn = {0956-540X}, doi = {10.1111/j.1365-246X.2010.04873.x}, pages = {541 -- 554}, year = {2011}, abstract = {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.}, language = {en} } @article{DonnerKruegerRoessleretal.2014, author = {Donner, Stefanie and Kr{\"u}ger, Frank and Roessler, Dirk and Ghods, Abdolreza}, title = {Combined Inversion of broadband and short-period waveform data for regional moment tensors: A case study in the Alborz Mountains, Iran}, series = {Bulletin of the Seismological Society of America}, volume = {104}, journal = {Bulletin of the Seismological Society of America}, number = {3}, publisher = {Seismological Society of America}, address = {Albany}, issn = {0037-1106}, doi = {10.1785/0120130229}, pages = {1358 -- 1373}, year = {2014}, abstract = {In this study, we suggest a novel approach for the retrieval of regional moment tensors for earthquakes with small to moderate magnitudes. The first modification is the combined inversion of broadband and short-period waveform data. The broadband waveforms are inverted in a frequency range suitable for surface waves, whereas for the short-period data a frequency range suitable for body waves is applied. The second modification is the use of first-motion body-wave polarities to select the most probable solution out of all solutions from inversion. To combine three different criteria for selecting the most probable solution (i.e., residual from inversion, double-couple content of solution, number of nonmatching first-motion body-wave polarities), the L2 norm is applied to the normalized parameters. We chose five earthquakes within the Alborz mountains, Iran, as a case study (3.1 <= M-w <= 4.1). In this area, several factors exacerbate the difficulty of performing inversion for moment tensors, for example, a heterogeneous station network and large azimuthal gaps. We have demonstrated that our approach supplies reliable moment tensors when inversion from broadband data alone fails. In one case, we successfully retrieved a stable solution from short-period waveform data alone. Thus, our approach enables successful determination of seismic moment tensors wherever a sparse network of broadband stations has thus far prevented it.}, language = {en} } @article{KruegerScherbaum2014, author = {Kr{\"u}ger, Frank and Scherbaum, Frank}, title = {The 29 September 1969, Ceres, South Africa, Earthquake: full waveform moment tensor inversion for point source and kinematic source parameters}, series = {Bulletin of the Seismological Society of America}, volume = {104}, journal = {Bulletin of the Seismological Society of America}, number = {1}, publisher = {Seismological Society of America}, address = {Albany}, issn = {0037-1106}, doi = {10.1785/0120130209}, pages = {576 -- 581}, year = {2014}, abstract = {The Ceres earthquake of 29 September 1969 is the largest known earthquake in southern Africa. Digitized analog recordings from Worldwide Standardized Seismographic Network stations (Powell and Fries, 1964) are used to retrieve the point source moment tensor and the most likely centroid depth of the event using full waveform modeling. A scalar seismic moment of 2.2-2.4 x 10(18) N center dot m corresponding to a moment magnitude of 6.2-6.3 is found. The analysis confirms the pure strike-slip mechanism previously determined from onset polarities by Green and Bloch (1971). Overall good agreement with the fault orientation previously estimated from local aftershock recordings is found. The centroid depth can be constrained to be less than 15 km. In a second analysis step, we use a higher order moment tensor based inversion scheme for simple extended rupture models to constrain the lateral fault dimensions. We find rupture propagated unilaterally for 4.7 s from east-southwest to west-northwest for about 17 km ( average rupture velocity of about 3: 1 km/s).}, language = {en} } @article{AlinaghiKruger2014, author = {Alinaghi, Alireza and Kruger, Frank}, title = {Seismic array analysis and redetermination of depths of earthquakes in Tien-Shan: implications for strength of the crust and lithosphere}, series = {Geophysical journal international}, volume = {198}, journal = {Geophysical journal international}, number = {2}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggu141}, pages = {1111 -- 1129}, year = {2014}, abstract = {We have redetermined focal depths of moderate and major earthquakes with reported lower-crust and upper-mantle depths that have occurred in Tien-Shan, since the availability of broad-band array data. Records of earthquakes at global arrays have been used for identification and modelling of depth phases in order to make accurate estimation of focal depths. Our results show that half of the purportedly deep earthquakes are indeed originating from depths attributable to middle-crust and lower-crust regions. Also one exceptional event in the northern foreland of Tien-Shan in Junggar Basin is located in the upper mantle at the depth of 64 km. Such unusually deep earthquakes for intraplate continental tectonic domain are all located at the margin of Tien-Shan with its adjacent stable blocks and at least some of them have occurred where the brittle behaviour of continental rocks is not highly expected. The reverse mechanisms of all these earthquakes and their proximity to formerly subducting and later colliding and underplating stable blocks and their interactions with overlying Tien-Shan are clues to explain this extremity.}, language = {en} }