@article{KnapmeyerEndrunKruegerGeissler2017, author = {Knapmeyer-Endrun, Brigitte and Kr{\"u}ger, Frank and Geissler, Wolfram H.}, title = {Upper mantle structure across the Trans-European Suture Zone imaged by S-receiver functions}, series = {Earth \& planetary science letters}, volume = {458}, journal = {Earth \& planetary science letters}, publisher = {Elsevier}, address = {Amsterdam}, organization = {PASSEQ Working Grp}, issn = {0012-821X}, doi = {10.1016/j.epsl.2016.11.011}, pages = {429 -- 441}, year = {2017}, abstract = {We present a high-resolution study of the upper mantle structure of Central Europe, including the western part of the East European Platform, based on S-receiver functions of 345 stations. A distinct contrast is found between Phanerozoic Europe and the East European Craton across the Trans-European Suture Zone. To the west, a pronounced velocity reduction with depth interpreted as lithosphere-asthenosphere boundary (LAB) is found at an average depth of 90 km. Beneath the craton, no strong and continuous LAB conversion is observed. Instead we find a distinct velocity reduction within the lithosphere, at 80-120 km depth. This mid-lithospheric discontinuity (MLD) is attributed to a compositional boundary between depleted and more fertile lithosphere created by late Proterozoic metasomatism. A potential LAB phase beneath the craton is very weak and varies in depth between 180 and 250 km, consistent with a reduced velocity contrast between the lower lithosphere and the asthenosphere. Within the Trans-European Suture Zone, lithospheric structure is characterized by strong heterogeneity. A dipping or step-wise increase to LAB depth of 150 km is imaged from Phanerozoic Europe to 20-22 degrees E, whereas no direct connection to the cratonic LAB or MLD to the east is apparent. At larger depths, a positive conversion associated with the lower boundary of the asthenosphere is imaged at 210-250 km depth beneath Phanerozoic Europe, continuing down to 300 km depth beneath the craton. Conversions from both 410 km and 660 km discontinuities are found at their nominal depth beneath Phanerozoic Europe, and the discontinuity at 410 km depth can also be traced into the craton. A potential negative conversion on top of the 410 km discontinuity found in migrated images is analyzed by modeling and attributed to interference with other converted phases.}, 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{LipkeZitzmannAmbergeretal.2007, author = {Lipke, Katrin and Zitzmann, Max and Amberger, Manuel and Ehlert, Carsten and R{\"o}ßler, Dirk and Kr{\"u}ger, Frank and Ohrnberger, Matthias}, title = {Traveltime residuals at regional and teleseismic distances for SE-Asia}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-14117}, year = {2007}, abstract = {Traveltime residuals for worldwide seismic stations are calculated. We use P and S waves from earthquakes in SE-Asia at teleseismic and regional distances. The obtained station residuals help to enhance earthquake localisation. Furthermore we calculated regional source dependent station residuals. They show a systematic dependence of the locality of the source. These source dependent residuals reflect heterogenities along the path and can be used for a refinement of earthquake localisation.}, language = {en} } @article{KnapmeyerEndrunKruegerLegendreetal.2013, author = {Knapmeyer-Endrun, Brigitte and Kr{\"u}ger, Frank and Legendre, C. P. and Geissler, Wolfram H.}, title = {Tracing the influence of the trans-european suture zone into the mantle transition zone}, series = {Earth \& planetary science letters}, volume = {363}, journal = {Earth \& planetary science letters}, publisher = {Elsevier}, address = {Amsterdam}, organization = {PASSEQ Working Grp}, issn = {0012-821X}, doi = {10.1016/j.epsl.2012.12.028}, pages = {73 -- 87}, year = {2013}, abstract = {Cratons with their thick lithospheric roots can influence the thermal structure, and thus the convective flow, in the surrounding mantle. As mantle temperatures are hard to measure directly, depth variations in the mantle transition zone (MTZ) discontinuities are often employed as a proxy. Here, we use a large new data set of P-receiver functions to map the 410 km and 660 km discontinuities beneath the western edge of the East European Craton and adjacent Phanerozoic Europe across the most fundamental lithospheric boundary in Europe, the Trans-European Suture Zone (TESZ). We observe significantly shorter travel times for conversions from both MTZ discontinuities within the craton, caused by the high velocities of the cratonic root. By contrast, the differential travel time across the MTZ is normal to only slightly raised. This implies that any insulating effect of the cratonic keel does not reach the MTZ. In contrast to earlier observations in Siberia, we do not find any trace of a discontinuity at 520 km depth, which indicates a rather dry MTZ beneath the western edge of the craton. Within most of covered Phanerozoic Europe, the MTZ differential travel time is remarkably uniform and in agreement with standard Earth models. No widespread thermal effects of the various episodes of Caledonian and Variscan subduction that took place during the amalgamation of the continent remain. Only more recent tectonic events, related to Alpine subduction and Quarternary volcanism in the Eifel area, can be traced. While the East European craton shows no distinct imprint into the MTZ, we discover the signature of the TESZ in the MTZ in the form of a linear region of about 350 km width with a 1.5 s increase in differential travel time, which could either be caused by high water content or decreased temperature. Taking into account results of recent S-wave tomographies, raised water content in the MTZ cannot be the main cause for this observation. Accordingly, we explain the increase, equivalent to a 15 km thicker MTZ, by a temperature decrease of about 80 K. We discuss two alternative models for this temperature reduction, either a remnant of subduction or an indication of downwelling due to small-scale, edge-driven convection caused by the contrast in lithospheric thickness across the TESZ. Any subducted lithosphere found in the MTZ at this location is unlikely to be related to Variscan subduction along the TESZ, though, as Eurasia has moved significantly northward since the Variscan orogeny.}, language = {en} } @article{KruegerGrosserBaumbachetal.1998, author = {Kr{\"u}ger, Frank and Grosser, H. and Baumbach, M. and Berckhemer, Hans}, title = {The Erzincan (Turkey) earthquake (Ms 6.8) of March 13, 1992 and its aftershock sequence}, year = {1998}, 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{RoesslerKruegerRuempkeretal.2006, author = {R{\"o}ßler, Dirk and Kr{\"u}ger, Frank and R{\"u}mpker, Georg and Psencik, Ivan}, title = {Tensile source components of swarm events in West Bohemia in 2000 by considering seismic anisotropy}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-12975}, year = {2006}, abstract = {Earthquake swarms occur frequently in West Bohemia, Central Europe. Their occurrence is correlated with and propably triggered by fluids that escape on the earth's surface near the epicentres. These fluids raise up periodically from a seemingbly deep-seated source in the upper mantle. Moment tensors for swarm events in 1997 indicate tensile faulting. However, they were determined under assumption of seismic isotropy although anisotropy can be observed. Anisotropy may obscure moment tensors and their interpretation. In 2000, more than 10,000 swarm earthquakes occurred near Novy Kostel, West Bohemia. Event triggering by fluid injection is likely. Activity lasted from 28/08 until 31/12/00 (9 phases) with maximum ML=3.2. High quality P-wave seismograms were used to retrieve the source mechanisms for 112 events between 28/08/00 and 30/10/00 using > 20 stations. We determine the source geometry using a new algorithm and different velocity models including anisotropy. From inversions of P waves we observe ML<3.2, strike-slip events on steep N-S oriented faults with additional normal or reverse components. Tensile components seem to be evident for more than 60\% of the processed swarm events in West Bohemia during the phases 1-7. Being most significant at great depths and at phases 1-4 during the swarm they are time and location dependent. Although tensile components are reduced when anisotropy is assumed they persist and seem to be important. They can be explained by pore-pressure changes due to the injection of fluids that raise up. Our findings agree with other observations e.g. correlation of fluid transport and seismicity, variations in b-value, forcing rate, and in pore pressure diffusion. Tests of our results show their significance.}, language = {en} } @article{HannemannKruegerDahmetal.2017, author = {Hannemann, Katrin and Kr{\"u}ger, Frank and Dahm, Torsten and Lange, Dietrich}, title = {Structure of the oceanic lithosphere and upper mantle north of the Gloria Fault in the eastern mid-Atlantic by receiver function analysis}, series = {Journal of geophysical research : Solid earth}, volume = {122}, journal = {Journal of geophysical research : Solid earth}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9313}, doi = {10.1002/2016JB013582}, pages = {7927 -- 7950}, year = {2017}, abstract = {Receiver functions (RF) have been used for several decades to study structures beneath seismic stations. Although most available stations are deployed on shore, the number of ocean bottom station (OBS) experiments has increased in recent years. Almost all OBSs have to deal with higher noise levels and a limited deployment time (approximate to 1year), resulting in a small number of usable records of teleseismic earthquakes. Here we use OBSs deployed as midaperture array in the deep ocean (4.5-5.5km water depth) of the eastern mid-Atlantic. We use evaluation criteria for OBS data and beamforming to enhance the quality of the RFs. Although some stations show reverberations caused by sedimentary cover, we are able to identify the Moho signal, indicating a normal thickness (5-8km) of oceanic crust. Observations at single stations with thin sediments (300-400m) indicate that a probable sharp lithosphere-asthenosphere boundary (LAB) might exist at a depth of approximate to 70-80km which is in line with LAB depth estimates for similar lithospheric ages in the Pacific. The mantle discontinuities at approximate to 410km and approximate to 660km are clearly identifiable. Their delay times are in agreement with PREM. Overall the usage of beam-formed earthquake recordings for OBS RF analysis is an excellent way to increase the signal quality and the number of usable events.}, language = {en} } @article{KruegerOhrnberger2005, author = {Kr{\"u}ger, Frank and Ohrnberger, Matthias}, title = {Spatio-temporal source characteristics of the 26 December 2004 Sumatra earthquake as imaged by teleseismic broadband arrays}, year = {2005}, abstract = {We test the capability of broadband arrays at teleseismic distances to image the spatio-temporal characteristics of the seismic energy release during the Dec 26, 2004 Sumatra earthquake at early observation times. Using a non-plane-wave array location technique previously reported values for rupture length (about 1150 km), duration (about 480 s), and average rupture velocity (2.4-2.7 km/s) are confirmed. Three dominant energy releases are identified: one near the hypocenter, a second at 6 degrees N94 degrees E about 130 s later and a third one after 300 s at 9 degrees N92-93 degrees E. The spatio-temporal distribution of the radiated seismic energy in the source region is calculated from the stacked broadband recordings of two arrays in Germany and Japan and results in rough estimates of the total seismic energy of 0.55.10(18) Nm (GRSN) and 1.53.10(18) Nm (FNET) respectively. Changes in the relative ratio of energy as function of spatio-temporal location indicate a rotation of the focal mechanism during the rupture process}, 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{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{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{SteinbergSudhausHeimannetal.2020, author = {Steinberg, Andreas and Sudhaus, Henriette and Heimann, Sebastian and Kr{\"u}ger, Frank}, title = {Sensitivity of InSAR and teleseismic observations to earthquake rupture segmentation}, series = {Geophysical journal international}, volume = {223}, journal = {Geophysical journal international}, number = {2}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggaa351}, pages = {875 -- 907}, year = {2020}, abstract = {Earthquakes often rupture across more than one fault segment. If such rupture segmentation occurs on a significant scale, a simple point-source or one-fault model may not represent the rupture process well. As a consequence earthquake characteristics inferred, based on one-source assumptions, may become systematically wrong. This might have effects on follow-up analyses, for example regional stress field inversions and seismic hazard assessments. While rupture segmentation is evident for most M-w > 7 earthquakes, also smaller ones with 5.5 < M-w < 7 can be segmented. We investigate the sensitivity of globally available data sets to rupture segmentation and their resolution to reliably estimate the mechanisms in presence of segmentation. We focus on the sensitivity of InSAR (Interferometric Synthetic Aperture Radar) data in the static near-field and seismic waveforms in the far-field of the rupture and carry out non-linear and Bayesian optimizations of single-source and two-sources kinematic models (double-couple point sources and finite, rectangular sources) using InSAR and teleseismic waveforms separately. Our case studies comprises of four M-w 6-7 earthquakes: the 2009 L'Aquila and 2016 Amatrice (Italy) and the 2005 and 2008 Zhongba (Tibet) earthquakes. We contrast the data misfits of different source complexity by using the Akaike informational criterion (AIC). We find that the AIC method is well suited for data-driven inferences on significant rupture segmentation for the given data sets. This is based on our observation that an AIC-stated significant improvement of data fit for two-segment models over one-segment models correlates with significantly different mechanisms of the two source segments and their average compared to the single-segment mechanism. We attribute these modelled differences to a sufficient sensitivity of the data to resolve rupture segmentation. Our results show that near-field data are generally more sensitive to rupture segmentation of shallow earthquakes than far-field data but that also teleseismic data can resolve rupture segmentation in the studied magnitude range. We further conclude that a significant difference in the modelled source mechanisms for different segmentations shows that an appropriate choice of model segmentation matters for a robust estimation of source mechanisms. It reduces systematic biases and trade-off and thereby improves the knowledge on the rupture. Our study presents a strategy and method to detect significant rupture segmentation such that an appropriate model complexity can be used in the source mechanism inference. A similar, systematic investigation of earthquakes in the range of M-w 5.5-7 could provide important hazard-relevant statistics on rupture segmentation. In these cases single-source models introduce a systematic bias. Consideration of rupture segmentation therefore matters for a robust estimation of source mechanisms of the studied earthquakes.}, language = {en} } @article{KinscherKruegerWoithetal.2013, author = {Kinscher, Jannes and Kr{\"u}ger, Frank and Woith, H. and L{\"u}hr, B. G. and Hintersberger, E. and Irmak, T. Serkan and Baris, S.}, title = {Seismotectonics of the Armutlu peninsula (Marmara Sea, NW Turkey) from geological field observation and regional moment tensor inversion}, series = {Tectonophysics : international journal of geotectonics and the geology and physics of the interior of the earth}, volume = {608}, journal = {Tectonophysics : international journal of geotectonics and the geology and physics of the interior of the earth}, number = {46}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0040-1951}, doi = {10.1016/j.tecto.2013.07.016}, pages = {980 -- 995}, year = {2013}, abstract = {The Armutlu peninsula, located in the eastern Marmara Sea, coincides with the western end of the rupture of the 17 August 1999, Izmit M-W 7.6 earthquake which is the penultimate event of an apparently westward migrating series of strong and disastrous earthquakes along the NAFZ during the past century. We present new seismotectonic data of this key region in order to evaluate previous seismotectonic models and their implications for seismic hazard assessment in the eastern Marmara Sea. Long term kinematics were investigated by performing paleo strain reconstruction from geological field investigations by morphotectonic and kinematic analysis of exposed brittle faults. Short term kinematics were investigated by inverting for the moment tensor of 13 small to moderate recent earthquakes using surface wave amplitude spectra. Our results confirm previous models interpreting the eastern Marmara Sea Region as an active transtensional pull-apart environment associated with significant NNE-SSW extension and vertical displacement. At the northern peninsula, long term deformation pattern did not change significantly since Pliocene times contradicting regional tectonic models which postulate a newly formed single dextral strike slip fault in the Marmara Sea Region. This area is interpreted as a horsetail splay fault structure associated with a major normal fault segment that we call the Waterfall Fault. Apart from the Waterfall Fault, the stress strain relation appears complex associated with a complicated internal fault geometry, strain partitioning, and reactivation of pre-existing plane structures. At the southern peninsula, recent deformation indicates active pull-apart tectonics constituted by NE-SW trending dextral strike slip faults. Earthquakes generated by stress release along large rupture zones seem to be less probable at the northern, but more probable at the southern peninsula. Additionally, regional seismicity appears predominantly driven by plate boundary stresses as transtensional faulting is consistent with the southwest directed far field deformation of the Anatolian plate. (C) 2013 Elsevier B.V. All rights reserved.}, language = {en} } @article{HannemannEulenfeldKruegeretal.2021, author = {Hannemann, Katrin and Eulenfeld, Tom and Kr{\"u}ger, Frank and Dahm, Torsten}, title = {Seismic scattering and absorption of oceanic lithospheric S waves in the Eastern North Atlantic}, series = {Geophysical journal international}, volume = {229}, journal = {Geophysical journal international}, number = {2}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggab493}, pages = {948 -- 961}, year = {2021}, abstract = {The scattering and absorption of high-frequency seismic waves in the oceanic lithosphere is to date only poorly constrained by observations. Such estimates would not only improve our understanding of the propagation of seismic waves, but also unravel the small-scale nature of the lithosphere and its variability. Our study benefits from two exceptional situations: (1) we deployed over 10 months a mid-aperture seismological array in the central part of the Eastern North Atlantic in 5 km water depth and (2) we could observe in total 340 high-frequency (up to 30 Hz) Po and So arrivals with tens to hundreds of seconds long seismic coda from local and regional earthquakes in a wide range of backazimuths and epicentral distances up to 850 km with a travel path in the oceanic lithosphere. Moreover, the array was located about 100 km north of the Gloria fault, defining the plate boundary between the Eurasian and African plates at this location which also allows an investigation of the influence of an abrupt change in lithospheric age (20 Ma in this case) on seismic waves. The waves travel with velocities indicating upper-mantle material. We use So waves and their coda of pre-selected earthquakes to estimate frequency-dependent seismic scattering and intrinsic attenuation parameters. The estimated scattering attenuation coefficients are between 10(-4) and 4 x 10(-5) m(-1) and are typical for the lithosphere or the upper mantle. Furthermore, the total quality factors for So waves below 5 Hz are between 20 and 500 and are well below estimates from previous modelling for observations in the Pacific Ocean. This implies that the Atlantic Ocean is more attenuative for So waves compared to the Pacific Ocean, which is inline with the expected behaviour for the lithospheric structures resulting from the slower spreading rates in the Atlantic Ocean. The results for the analysed events indicate that for frequencies above 3 Hz, intrinsic attenuation is equal to or slightly stronger than scattering attenuation and that the So-wave coda is weakly influenced by the oceanic crust. Both observations are in agreement with the proposed propagation mechanism of scattering in the oceanic mantle lithosphere. Furthermore, we observe an age dependence which shows that an increase in lithospheric age is associated with a decrease in attenuation. However, we also observe a trade-off of this age-dependent effect with either a change in lithospheric thickness or thermal variations, for example due to small-scale upwellings in the upper mantle in the southeast close to Madeira and the Canaries. Moreover, the influence of the nearby Gloria fault is visible in a reduction of the intrinsic attenuation below 3 Hz for estimates across the fault. This is the first study to estimate seismic scattering and absorption parameters of So waves for an area with several hundreds of kilometres radius centred in the Eastern North Atlantic and using them to characterize the nature of the oceanic lithosphere.}, language = {en} } @article{KitoKruegerNegishi2004, author = {Kito, Tadashi and Kr{\"u}ger, Frank and Negishi, H.}, title = {Seismic heterogeneous structure in the lowermost mantle beneath the southwestern Pacific}, year = {2004}, abstract = {The P and S wave velocity structure of the D" layer beneath the southwestern Pacific was investigated by using short-period data from 12 deep events in the Tonga-Fiji region recorded by the J-Array and the Hi-net (two large- aperture seismic arrays) in Japan. Reflected wave beam forming (RWB) and a migration method were used to extract weak signals originating from heterogeneities in the lowermost mantle. In order to acquire high resolution a double-array method was applied to the data. The results of the RWB method indicate that seismic energy is reflected at discontinuities near the depths of 2520 and 2650 km, which have a negative P wave velocity contrast of 1\% at the most. In addition, there is a positive seismic discontinuity at a depth of 2800 km. In the case of the S wave, reflected energy is produced almost at the same depth (2550 km depth). An apparent depth shift (50 km) of the discontinuity at the depth of 2850 km may indicate that the S wave velocity reduction in the lowermost mantle is similar to2-3 times stronger than that of P. A two-dimensional cross section, constructed with the RWB method, suggests that the observed discontinuities can be characterized as intermittent lateral heterogeneities whose lateral extent is a few hundred kilometers. The migration shows weak evidence of scattering objects which belong to the seismic discontinuities detected by the RWB method. These anomalous structures may represent a part of hot plume generated beneath the southwestern Pacific in the lowermost mantle}, 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} } @article{DonnerRoesslerKruegeretal.2013, author = {Donner, Stefanie and R{\"o}ßler, Dirk and Kr{\"u}ger, Frank and Ghods, Abdolreza and Strecker, Manfred}, title = {Segmented seismicity of the M (w) 6.2 Baladeh earthquake sequence (Alborz Mountains, Iran) revealed from regional moment tensors}, series = {Journal of seismology}, volume = {17}, journal = {Journal of seismology}, number = {3}, publisher = {Springer}, address = {Dordrecht}, issn = {1383-4649}, doi = {10.1007/s10950-013-9362-7}, pages = {925 -- 959}, year = {2013}, abstract = {The M (w) 6.2 Baladeh earthquake occurred on 28 May 2004 in the Alborz Mountains, northern Iran. This earthquake was the first strong shock in this intracontinental orogen for which digital regional broadband data are available. The Baladeh event provides a rare opportunity to study fault geometry and ongoing deformation processes using modern seismological methods. A joint inversion for hypocentres and a velocity model plus a surface-wave group dispersion curve analysis were used to obtain an adapted velocity model, customised for mid- and long-period waveform modelling. Based on the new velocity model, regional waveform data of the mainshock and larger aftershocks (M (w) a parts per thousand yen3.3) were inverted for moment tensors. For the Baladeh mainshock, this included inversion for kinematic parameters. All analysed earthquakes show dominant thrust mechanisms at depths between 14 and 26 km, with NW-SE striking fault planes. The mainshock ruptured a 28A degrees south-dipping area of 24 x 21 km along a north-easterly direction. The rupture plane of the mainshock does not coincide with the aftershock distribution, neither in map view nor with respect to depth. The considered aftershocks form two main clusters. The eastern cluster is associated with the mainshock. The western cluster does not appear to be connected with the rupture plane of the mainshock but, instead, indicates a second activated fault plane dipping at 85A degrees towards the north.}, 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{RoesslerKruegerOhrnberger2007, author = {R{\"o}ßler, Dirk and Kr{\"u}ger, Frank and Ohrnberger, Matthias}, title = {Rupture propagation of the TsE (Mw7.7) on 17 July 2006 off-coast Java}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-12964}, year = {2007}, abstract = {The Mw=7.7 tsunamogenic earthquake (TsE) on 17 July 2006, 08:19:28 shock the Indian Ocean at about 15 km depth off-coast Java, Indonesia. It caused a local tsunami with wave heights exceeding 2 m. The death toll reached several hundred. Thousands of people were displaced. By means of standard array methods, we have investigated the propagation and the extent of the rupture front of the causative earthquake. Waveform similarity is expressed by means of the semblance. We back-propagate the semblance for first-arrival phases recorded at broad-band stations within teleseismic distances (30°-95°). Image enhancement is realised by stacking the semblance of 8 arrays within different epicentral and azimuthal directions. From teleseismic observations we find rupturing of a 200 x 100 km wide area in at least 2 phases with propagation from NW to SE and source duration >125 s. The event has some characteristics of a circular rupture followed by unilateral faulting with change in slip rate. Unusually slow rupturing (≈1.5 km/s) is indicated. Fault area and aftershock distribution coincide. Spatial and temporal resolution are frequency dependent. Studies of a Mw6.0 earthquake on 2006/09/21 and one synthetic source show a ≈1° limit in resolution. Retrieved source area, source duration as well as peak values for semblance and beam power increase with the size of the earthquake making possible an automatic detection and classification of large and small earthquakes.}, language = {en} }