@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{WeidleWiesenbergElSharkawyetal.2022, author = {Weidle, Christian and Wiesenberg, Lars and El-Sharkawy, Amr and Kr{\"u}ger, Frank and Scharf, Andreas and Agard, Philippe and Meier, Thomas}, title = {A 3-D crustal shear wave velocity model and Moho map below the Semail Ophiolite, eastern Arabia}, series = {Geophysical journal international}, volume = {231}, journal = {Geophysical journal international}, number = {2}, publisher = {Oxford University Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggac223}, pages = {817 -- 834}, year = {2022}, abstract = {The Semail Ophiolite in eastern Arabia is the largest and best-exposed slice of oceanic lithosphere on land. Detailed knowledge of the tectonic evolution of the shallow crust, in particular during and after ophiolite obduction in Late Cretaceous times is contrasted by few constraints on physical and compositional properties of the middle and lower continental crust below the obducted units. The role of inherited, pre-obduction crustal architecture remains therefore unaccounted for in our understanding of crustal evolution and the present-day geology. Based on seismological data acquired during a 27-month campaign in northern Oman, Ambient Seismic Noise Tomography and Receiver Function analysis provide for the first time a 3-D radially anisotropic shear wave velocity (V-S) model and a consistent Moho map below the iconic Semail Ophiolite. The model highlights deep crustal boundaries that segment the eastern Arabian basement in two distinct units. The previously undescribed Western Jabal Akhdar Zone separates Arabian crust with typical continental properties and a thickness of similar to 40-45 km in the northwest from a compositionally different terrane in the southeast that is interpreted as a terrane accreted during the Pan-African orogeny in Neoproterozoic times. East of the Ibra Zone, another deep crustal boundary, crustal thickness decreases to 30-35 km and very high lower crustal V-S suggest large-scale mafic intrusions into, and possible underplating of the Arabian continental crust that occurred most likely during Permian breakup of Pangea. Mafic reworking is sharply bounded by the (upper crustal) Semail Gap Fault Zone, northwest of which no such high velocities are found in the crust. Topography of the Oman Mountains is supported by a mild crustal root and Moho depth below the highest topography, the Jabal Akhdar Dome, is similar to 42 km. Radial anisotropy is robustly resolved in the upper crust and aids in discriminating dipping allochthonous units from autochthonous sedimentary rocks that are indistinguishable by isotropic V-S alone. Lateral thickness variations of the ophiolite highlight the Haylayn Ophiolite Massif on the northern flank of Jabal Akhdar Dome and the Hawasina Window as the deepest reaching unit. Ophiolite thickness is similar to 10 km in the southern and northern massifs, and <= 5 km elsewhere.}, language = {en} } @article{WeberWicksKruegeretal.1998, author = {Weber, Michael H. and Wicks, Chuck and Kr{\"u}ger, Frank and Jahnke, Gunnar and Schlittenhardt, J{\"o}rg}, title = {Asymmetric radiation of seismic waves from an atoll : nuclear tests in French Polynesia}, year = {1998}, language = {en} } @article{WeberWicksKruegeretal.2000, author = {Weber, Michael H. and Wicks, Chuck and Kr{\"u}ger, Frank and Jahnke, Gunnar and Baumann, M.}, title = {Reply to comment on "Asymmetric radiation of seismic waves from an atoll : Nuclear tests in French Polynesia" by A. Douglas}, year = {2000}, 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{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{ScherbaumKruegerWeber1997, author = {Scherbaum, Frank and Kr{\"u}ger, Frank and Weber, Michael H.}, title = {Double beam imaging : mapping lower mantle heterogeneities using combinations of source and receiver arrays}, year = {1997}, language = {en} } @article{RoesslerRumpkerKrueger2004, author = {R{\"o}ßler, Dirk and Rumpker, Georg and Kr{\"u}ger, Frank}, title = {Ambiguous moment tensors and radiation patterns in anisotropic media with applications to the modeling of earthquake mechanisms in W-Bohemia}, issn = {0039-3169}, year = {2004}, abstract = {Anisotropic material properties are usually neglected during inversions for source parameters of earthquakes. In general anisotropic media, however, moment tensors for pure-shear sources can exhibit significant non-double-couple components. Such effects may be erroneously interpreted as an indication for volumetric changes at the source. Here we investigate effects of anisotropy on seismic moment tensors and radiation patterns for pure-shear and tensile-type sources. Anisotropy can significantly influence the interpretation of the source mechanisms. For example, the orientation of the slip within the fault plane may affect the total seismic moment. Also, moment tensors due to pure- shear and tensile faulting can have similar characteristics depending on the orientation of the elastic tensor. Furthermore, the tensile nature of an earthquake can be obscured by near-source anisotropic properties. As an application, we consider effects of inhomogeneous anisotropic properties on the seismic moment tensor and the radiation patterns of a selected type of micro-earthquakes observed in W-Bohemia. The combined effects of near-source and along- path anisotropy cause characteristic amplitude distortions of the P, S1 and S2 waves. However, the modeling suggests that neither homogeneous nor inhomogeneous anisotropic properties alone can explain the observed large non-double-couple components. The results also indicate that a correct analysis of the source mechanism, in principle, is achievable by application of anisotropic moment tensor inversion}, 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} } @misc{RoesslerKruegerOhrnbergeretal.2008, author = {R{\"o}ßler, Dirk and Kr{\"u}ger, Frank and Ohrnberger, Matthias and Ehlert, Lutz}, title = {Automatic near real-time characterisation of large earthquakes}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-20191}, year = {2008}, abstract = {An der Universit{\"a}t Potsdam wird seit 2008 ein automatisiertes Verfahren angewandt, um Bruchparamter großer Erdbeben in quasi-Echtzeit, d.h. wenige Minuten nachdem sich das Beben ereignet hat, zu bestimmen und der {\"O}ffentlichkeit via Internet zur Verf{\"u}gung zu stellen. Es ist vorgesehen, das System in das Deutsch-Indonesische Tsunamifr{\"u}hwarnsystem (GITEWS) zu integrieren, f{\"u}r das es speziell konfiguriert ist. Wir bestimmen insbesondere die Dauer und die Ausdehnung des Erdbebens, sowie dessen Bruchgeschwindigkeit und -richtung. Dabei benutzen wir die Seismogramme der zuerst eintreffenden P Wellen vom Breitbandstationen in teleseimischer Entfernung vom Beben sowie herk{\"o}mmliche Arrayverfahren in teilweise modifizierter Form. Die Semblance wir als {\"A}hnlichkeitsmaß verwendet, um Seismogramme eines Stationsnetzes zu vergleichen. Im Falle eines Erdbebens ist die Semblance unter Ber{\"u}cksichtigung des Hypozentrums zur Herdzeit und w{\"a}hrend des Bruchvorgangs deutlich zeitlich und r{\"a}umlich erh{\"o}ht und konzentriert. Indem wir die Ergebnisse verschiedener Stationsnetzwerke kombinieren, erreichen wir Unabh{\"a}ngigkeit von der Herdcharakteristik und eine raum-zeitliche Aufl{\"o}sung, die es erlaubt die o.g. Parameter abzuleiten. In unserem Beitrag skizzieren wir die Methode. Anhand der beiden M8.0 Benkulu Erdbeben (Sumatra, Indonesien) vom 12.09.2007 und dem M8.0 Sichuan Ereignis (China) vom 12.05.2008 demonstrieren wir Aufl{\"o}sungsm{\"o}glichkeiten und vergleichen die Ergebnisse der automatisierten Echtzeitanwendung mit nachtr{\"a}glichen Berechnungen. Weiterhin stellen wir eine Internetseite zur Verf{\"u}gung, die die Ergebnisse pr{\"a}sentiert und animiert. Diese kann z.B. in geowissenschaftlichen Einrichtungen an Computerterminals gezeigt werden. Die Internetauftritte haben die folgenden Adressen: http://www.geo.uni-potsdam.de/arbeitsgruppen/Geophysik_Seismologie/forschung/ruptrack/openday http://www.geo.uni-potsdam.de/arbeitsgruppen/Geophysik_Seismologie/forschung/ruptrack}, language = {en} } @misc{RoesslerKruegerOhrnberger2007, author = {R{\"o}ßler, Dirk and Kr{\"u}ger, Frank and Ohrnberger, Matthias}, title = {Rupture propagation of recent large TsE off-coast Sumatra and Java}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-13039}, year = {2007}, abstract = {The spatio-temporal evolution of the three recent tsunamogenic earthquakes (TsE) off-coast N-Sumatra (Mw9.3), 28/03/2005 (Mw8.5) off-coast Nias, on 17/07/2006 (Mw7.7) off-coast Java. Start time, duration, and propagation of the rupture are retrieved. All parameters can be obtained rapidly after recording of the first-arrival phases in near-real time processing. We exploit semblance analysis, backpropagation and broad-band seismograms within 30°-95° distance. Image enhancement is reached by stacking the semblance of arrays within different directions. For the three events, the rupture extends over about 1150, 150, and 200km, respectively. The events in 2004, 2005, and 2006 had source durations of at least 480s, 120s, and 180s, respectively. We observe unilateral rupture propagation for all events except for the rupture onset and the Nias event, where there is evidence for a bilateral start of the rupture. Whereas average rupture speed of the events in 2004 and 2005 is in the order of the S-wave speed (≈2.5-3km/s), unusually slow rupturing (≈1.5 km/s) is indicated for the July 2006 event. For the July 2006 event we find rupturing of a 200 x 100 km wide area in at least 2 phases with propagation from NW to SE. The event has some characteristics of a circular rupture followed by unilateral faulting with change in slip rate. 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 generally increase with the size of the earthquake making possible an automatic detection and classification of large and small earthquakes.}, 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} } @misc{RoesslerKruegerOhrnberger2008, author = {R{\"o}ßler, Dirk and Kr{\"u}ger, Frank and Ohrnberger, Matthias}, title = {Automatic near real-time characterisation of large earthquakes}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-18382}, year = {2008}, abstract = {We use seismic array methods (semblance analysis) to image areas of seismic energy release in the Sunda Arc region and world-wide. Broadband seismograms at teleseismic distances (30° ≤ Δ ≤ 100°) are compared at several subarrays. Semblance maps of different subarrays are multiplied. High semblance tracked over long time (10s of second to minutes) and long distances indicate locations of earthquakes. The method allows resolution of rupture characteristics important for tsunami early warning: start and duration, velocity and direction, length and area. The method has been successfully applied to recent and historic events (M>6.5) and is now operational in real time. Results are obtained shortly after source time, see http://www.geo.uni-potsdam.de/Forschung/Geophysik/GITEWS/tsunami.htm). Comparison of manual and automatic processing are in good agreement. Computational effort is small. Automatic results may be obtained within 15 - 20 minutes after event occurrence.}, language = {en} } @misc{RoesslerKruegerOhrnberger2008, author = {R{\"o}ßler, Dirk and Kr{\"u}ger, Frank and Ohrnberger, Matthias}, title = {Rupture Propagation of the 2008/05/12 Ms8.0 Wenchuan Earthquake}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-29195}, year = {2008}, abstract = {We study the rupture propagation of the 2008/05/12 Ms8.0 Wenchuan Earthquake. We apply array techniques such as semblance vespagram analysis to P waves recorded at seismic broadband station within 30-100° epicentral distance. By combination of multiple large aperture station groups spatial and temporal resolution is enhanced and problems due source directivity and source mechanism are avoided. We find that seismic energy was released for at least 110 s. Propagating unilaterally at sub-shear rupture velocity of about 2.5 km/s in NE direction, the earthquake reaches a lateral extent of more than 300 km. Whereas high semblance during within 70 s from rupture start indicates simple propagation more complex source processes are indicated thereafter by decreases coherency in seismograms. At this stage of the event coherency is low but significantly above noise level. We emphasize that first result of our computations where obtain within 30 minutes after source time by using an atomized algorithm. This procedure has been routinely and globally applied to major earthquakes. Results are made public through internet.}, language = {en} } @misc{RoesslerHiemerBachetal.2009, author = {R{\"o}ßler, Dirk and Hiemer, Stephan and Bach, Christoph and Delavaud, Elise and Kr{\"u}ger, Frank and Ohrnberger, Matthias and Sauer, David and Scherbaum, Frank and Vollmer, Daniel}, title = {Small-aperture seismic array monitors Vogtland earthquake swarm in 2008/09}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-29185}, year = {2009}, abstract = {The most recent intense earthquake swarm in the Vogtland lasted from 6 October 2008 until January 2009. Greatest magnitudes exceeded M3.5 several times in October making it the greatest swarm since 1985/86. In contrast to the swarms in 1985 and 2000, seismic moment release was concentrated near swarm onset. Focal area and temporal evolution are similar to the swarm in 2000. Work hypothysis: uprising upper-mantle fluids trigger swarm earthquakes at low stress level. To monitor the seismicity, the University of Potsdam operated a small aperture seismic array at 10 km epicentral distance between 18 October 2008 and 18 March 2009. Consisting of 12 seismic stations and 3 additional microphones, the array is capable of detecting earthquakes from larger to very low magnitudes (M<-1) as well as associated air waves. We use array techniques to determine properties of the incoming wavefield: noise, direct P and S waves, and converted phases.}, language = {en} } @article{PreussevanderMeerDeshpandeetal.2011, author = {Preusse, Franziska and van der Meer, Elke and Deshpande, Gopikrishna and Kr{\"u}ger, Frank and Wartenburger, Isabell}, title = {Fluid intelligence allows flexible recruitment of the parieto-frontal network in analogical reasoning}, series = {Frontiers in human neuroscienc}, volume = {5}, journal = {Frontiers in human neuroscienc}, number = {3}, publisher = {Frontiers Research Foundation}, address = {Lausanne}, issn = {1662-5161}, doi = {10.3389/fnhum.2011.00022}, pages = {14}, year = {2011}, abstract = {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.}, language = {en} } @article{PaloTilmannKruegeretal.2014, author = {Palo, Mauro and Tilmann, Frederik and Kr{\"u}ger, Frank and Ehlert, Lutz and Lange, Dietrich}, title = {High-frequency seismic radiation from Maule earthquake (M-w 8.8, 2010 February 27) inferred from high-resolution backprojection analysis}, 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/ggu311}, pages = {1058 -- 1077}, year = {2014}, abstract = {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.}, 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{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{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{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{LontsiGarciaJerezCamiloMolinaVillegasetal.2019, author = {Lontsi, Agostiny Marrios and Garcia-Jerez, Antonio and Camilo Molina-Villegas, Juan and Jose Sanchez-Sesma, Francisco and Molkenthin, Christian and Ohrnberger, Matthias and Kr{\"u}ger, Frank and Wang, Rongjiang and Fah, Donat}, title = {A generalized theory for full microtremor horizontal-to-vertical [H/V(z,f)] spectral ratio interpretation in offshore and onshore environments}, series = {Geophysical journal international}, volume = {218}, journal = {Geophysical journal international}, number = {2}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggz223}, pages = {1276 -- 1297}, year = {2019}, abstract = {Advances in the field of seismic interferometry have provided a basic theoretical interpretation to the full spectrum of the microtremor horizontal-to-vertical spectral ratio [H/V(f)]. The interpretation has been applied to ambient seismic noise data recorded both at the surface and at depth. The new algorithm, based on the diffuse wavefield assumption, has been used in inversion schemes to estimate seismic wave velocity profiles that are useful input information for engineering and exploration seismology both for earthquake hazard estimation and to characterize surficial sediments. However, until now, the developed algorithms are only suitable for on land environments with no offshore consideration. Here, the microtremor H/V(z, f) modelling is extended for applications to marine sedimentary environments for a 1-D layered medium. The layer propagator matrix formulation is used for the computation of the required Green's functions. Therefore, in the presence of a water layer on top, the propagator matrix for the uppermost layer is defined to account for the properties of the water column. As an application example we analyse eight simple canonical layered earth models. Frequencies ranging from 0.2 to 50 Hz are considered as they cover a broad wavelength interval and aid in practice to investigate subsurface structures in the depth range from a few meters to a few hundreds of meters. Results show a marginal variation of 8 per cent at most for the fundamental frequency when a water layer is present. The water layer leads to variations in H/V peak amplitude of up to 50 per cent atop the solid layers.}, 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} } @misc{LipkeKruegerRoessler2008, author = {Lipke, Katrin and Kr{\"u}ger, Frank and R{\"o}ßler, Dirk}, title = {Subduction zone structure along Sumatra from receiver functions}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-18260}, year = {2008}, abstract = {Receiver functions are a good tool to investigate the seismotectonic structure beneath the a seismic station. In this study we apply the method to stations situated on or near Sumatra to find constraints on a more detailed velocity model which should improve earthquake localisation. We estimate shallow Moho-depths (~ 21 km) close to the trench and depths of ~30 km at greater distances. First evidences for the dip direction of the slab of ~60° are provided. Receiver functions were calculated for 20 stations for altogether 110 earthquakes in the distance range between 30° and 95° from the receiver. However the number of receiver functions per station is strongly variable as it depends on the installation date, the signal-to-noise-ratio of the station and the reliability of the acquisition.}, 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{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{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{KulikovaKrueger2017, author = {Kulikova, Galina and Kr{\"u}ger, Frank}, title = {Historical seismogram reproductions for the source parameters determination of the 1902, Atushi (Kashgar) earthquake}, series = {Journal of seismology}, volume = {21}, journal = {Journal of seismology}, publisher = {Springer}, address = {Dordrecht}, issn = {1383-4649}, doi = {10.1007/s10950-017-9683-z}, pages = {1577 -- 1597}, year = {2017}, abstract = {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.}, 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} } @misc{KruegerOhrnbergerRoessler2008, author = {Kr{\"u}ger, Frank and Ohrnberger, Matthias and R{\"o}ßler, Dirk}, title = {Rupture imaging of large earthquakes with a poststack isochrone migration method}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-18395}, year = {2008}, abstract = {Rapid and robust characterization of large earthquakes in terms of their spatial extent and temporal duration is of high importance for disaster mitigation and early warning applications. Backtracking of seismic P-waves was successfully used by several authors to image the rupture process of the great Sumatra earthquake (26.12.2004) using short period and broadband arrays. We follow here an approach of Walker et al. to backtrack and stack broadband waveforms from global network stations using traveltimes for a global Earth model to obtain the overall spatio-temporal development of the energy radiation of large earthquakes in a quick and robust way. We present results for selected events with well studied source processes (Kokoxili 14.11.2001, Tokachi-Oki 25.09.2003, Nias 28.03.2005). Further, we apply the technique in a semi-real time fashion to broadband data of earthquakes with a broadband magnitude >= 7 (roughly corresponding to Mw 6.5). Processing is based on first automatic detection messages from the GEOFON extended virtual network (GEVN).}, 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{KruegerKulikovaLandgraf2018, author = {Kr{\"u}ger, Frank and Kulikova, Galina and Landgraf, Angela}, title = {Magnitudes for the historical 1885 (Belovodskoe), the 1887 (Verny) and the 1889 (Chilik) earthquakes in Central Asia determined from magnetogram recordings}, series = {Geophysical journal international}, volume = {215}, journal = {Geophysical journal international}, number = {3}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggy377}, pages = {1824 -- 1840}, year = {2018}, abstract = {Six large magnitude earthquakes in Central Asia which occurred at the end of the 19th century were recorded on early magnetographs in Great Britain. Scalar seismic moment estimates of the 1911 Chon-Kemin, the 1902 Atushi and the 1907 Karatag earthquakes in Central Asia were recently determined by historical seismogram modelling. For those events, we find agreement between moment magnitudes estimated from seismograms and from magnetograms. This supports the assumption of linear scaling of magnetogram amplitudes as function of M-0, which we then use to estimate the moment magnitudes for earlier large-magnitude events, that is, the 1885 Belovodskoe, 1887 Verny and 1889 Chilik earthquakes. The magnetometer data imply that the Chilik earthquake had M(W)7.9, slightly smaller than the Chon-Kemin event with M(W)8.0. The Verny earthquake, however, for which we estimate M(W)7.7, is likely larger than listed in catalogues (M7.3). Similarly, we find a larger magnitude M(W)7.6 (instead of the previous M6.9) for the Belovodskoe earthquake, but this remains uncertain due to measurement imprecision.}, language = {en} } @article{KruegerKulikovaLandgraf2017, author = {Kr{\"u}ger, Frank and Kulikova, Galina and Landgraf, Angela}, title = {Instrumental magnitude constraints for the 11 July 1889, Chilik earthquake}, series = {Seismicity, fault rupture and earthquake hazards in slowly deforming regions}, volume = {432}, journal = {Seismicity, fault rupture and earthquake hazards in slowly deforming regions}, publisher = {The Geological Society}, address = {London}, isbn = {978-1-86239-745-3}, issn = {0305-8719}, doi = {10.1144/SP432.8}, pages = {41 -- 72}, year = {2017}, abstract = {A series of large-magnitude earthquakes above 6.9 occurred in the northern Tien-Shan between 1885 and 1911. The Chilik earthquake of 11 July 1889, has been listed with a magnitude of 8.3, based on sparse macroseismic intensities, constrained by reported damage. Despite the existence of several juvenile fault scarps in the epicentral region, that are possibly associated with the 1889 earthquake, no through-going surface rupture having the dimensions expected for a magnitude 8.3 earthquake has been located - a puzzling dilemma. Could the magnitude have been overestimated? This would have major implications not only for the understanding of the earthquake series, but also for regional hazard estimates. Fortunately, a fragmentary record from an early Rebeur-Paschwitz seismometer exists for the Chilik event, recorded in Wilhelmshaven (Germany). To constrain the magnitude, we compare the late coda waves of this record with those of recent events from Central Asia, recorded on modern instruments in Germany and filtered with Rebeur-Paschwitz instrument characteristics. Additional constraints come from disturbances of historic magnetograms that exist from the Chilik and the 1911 Chon-Kemin earthquakes. Scaling of these historic records confirm a magnitude of about 8 for the 1889 Chilik earthquake, pointing towards a lower crustal contribution to the fault area.}, 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{KruegerDahmHannemann2020, author = {Kr{\"u}ger, Frank and Dahm, Torsten and Hannemann, Katrin}, title = {Mapping of Eastern North Atlantic Ocean seismicity from Po/So observations at a mid-aperture seismological broad-band deep sea array}, series = {Geophysical journal international}, volume = {221}, journal = {Geophysical journal international}, number = {2}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggaa054}, pages = {1055 -- 1080}, year = {2020}, abstract = {A mid-aperture broad-band test array (OBS array DOCTAR) was deployed from June 2011 to April 2012 about 100 km north of the Gloria fault in the Eastern North Atlantic in about 5000 m water depth. In addition arrays were installed on Madeira Island and in western Portugal mainland. For the first time in the Eastern North Atlantic, we recorded a large number of high frequency Po and So waves from local and regional small and moderate earthquakes (M-L < 4). An incoherent beamforming method was adapted to scan continuous data for such Po and So arrivals applying a sliding window waveform migration and frequency-wavenumber technique. We identify about 320 Po and 1550 So arrivals and compare the phase onsets with the ISC catalogue (ISC 2015) for the same time span. Up to a distance of 6 degrees to the DOCTAR stations all events listed in the ISC catalogue could be associated to Po and So phases. Arrivals from events in more than 10 degrees distance could be identified only in some cases. Only few Po and/or So arrivals were detected for earthquakes from the European and African continental area, the continental shelf regions and for earthquakes within or northwest of the Azores plateau. Unexpectedly, earthquake clusters are detected within the oceanic plates north and south of the Gloria fault and far from plate boundaries, indicating active intraplate structures. We also observe and locate numerous small magnitude earthquakes on the segment of the Gloria fault directly south of DOCTAR, which likely coincides with the rupture of the 25 November 1941 event. Local small magnitude earthquakes located beneath DOCTAR show hypocentres up to 30 km depth and strike-slip focal mechanisms. A comparison with detections at temporary mid-aperture arrays on Madeira and in western Portugal shows that the deep ocean array performs much better than the island and the continental array regarding the detection threshold for events in the oceanic plates. We conclude that sparsely distributed mid-aperture seismic arrays in the deep ocean could decrease the detection and location threshold for seismicity with M-L < 4 in the oceanic plate and might constitute a valuable tool to monitor oceanic plate seismicity.}, language = {en} } @article{KruegerDahm1999, author = {Kr{\"u}ger, Frank and Dahm, Torsten}, title = {Higher degree moment inversin using far-field broad-band recordings : theory and evaluation of the method with application to the 1994 Bolivia deep earthauke}, year = {1999}, language = {en} } @book{Krueger2009, author = {Kr{\"u}ger, Frank}, title = {Seismologische Arrays in der teleseismischen Struktur und Herdprozessabbildung : Antrittsvorlesung 2009-06-18}, publisher = {Univ.-Bibl.}, address = {Potsdam}, year = {2009}, abstract = {Frank Kr{\"u}ger, der seit 1997 am Institut f{\"u}r Geowissenschaften in der Seismologie t{\"a}tig ist, wurde eine außerplanm{\"a}ßige Professur f{\"u}r Geophysik {\"u}bertragen. Seine Vorlesung zum Thema "Seismologische Arrays in der teleseismischen Struktur- und Herdprozessabbildung" gibt einen {\"U}berblick {\"u}ber den Einsatz arrayseismologischer Verfahren in der Erdstrukturerkundung und der detaillierten Abbildung von Erdbebenherden. Dichte Netzwerke von seismologischen Messstationen erm{\"o}glichen die Anwendung spezieller hochaufl{\"o}sender Auswertungsverfahren. Diese wurden zun{\"a}chst im Kontext der Spionage entwickelt, finden heutzutage aber breite Anwendung in vielen seismologischen Forschungsbereichen, von der Erkundung der Struktur an der Grenze zum Erdkern, {\"u}ber Zusammenh{\"a}nge von seismischem Wellenfeld und Wetterph{\"a}nomenen bis hin zum Einsatz bei Tsunamifr{\"u}hwarnsystemen zur schnellen Analyse sehr großer Erdbebenherde.}, language = {de} } @article{KrugerOhrnberger2005, author = {Kruger, Frank and Ohrnberger, Matthias}, title = {Tracking the rupture of the M-w=9.3 Sumatra earthquake over 1,150 km at teleseismic distance}, issn = {0028-0836}, year = {2005}, abstract = {On 26 December 2004, a moment magnitude M-w = 9.3 earthquake occurred along Northern Sumatra, the Nicobar and Andaman islands, resulting in a devastating tsunami in the Indian Ocean region(1). The rapid and accurate estimation of the rupture length and direction of such tsunami-generating earthquakes is crucial for constraining both tsunami wave- height models as well as the seismic moment of the events. Compressional seismic waves generated at the hypocentre of the Sumatra earthquake arrived after about 12 min at the broadband seismic stations of the German Regional Seismic Network (GRSN)(2,3), located approximately 9,000 km from the event. Here we present a modification of a standard array- seismological approach and show that it is possible to track the propagating rupture front of the Sumatra earthquake over a total rupture length of 1,150 km. We estimate the average rupture speed to be 2.3-2.7 km s(-1) and the total duration of rupture to be at least 430 s, and probably between 480 and 500 s.}, language = {en} } @misc{KriegerowskiCescaOhrnbergeretal.2018, author = {Kriegerowski, Marius and Cesca, Simone and Ohrnberger, Matthias and Dahm, Torsten and Kr{\"u}ger, Frank}, title = {Event couple spectral ratio Q method for earthquake clusters}, series = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, number = {683}, doi = {10.25932/publishup-42602}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-426029}, pages = {12}, year = {2018}, abstract = {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.}, language = {en} } @article{KriegerowskiCescaOhrnbergeretal.2019, author = {Kriegerowski, Marius and Cesca, Simone and Ohrnberger, Matthias and Dahm, Torsten and Kr{\"u}ger, Frank}, title = {Event couple spectral ratio Q method for earthquake clusters}, series = {Solid Earth}, journal = {Solid Earth}, number = {10}, publisher = {Copernicus Publications}, address = {G{\"o}ttingen}, issn = {1869-9529}, doi = {10.5194/se-10-317-2019}, pages = {317 -- 328}, year = {2019}, abstract = {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.}, 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{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{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{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{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{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{KastleWeberKrueger2017, author = {Kastle, Emanuel D. and Weber, Michael and Kr{\"u}ger, Frank}, title = {Complex Deep Structure of the African Low-Velocity Zone}, series = {Bulletin of the Seismological Society of America}, volume = {107}, journal = {Bulletin of the Seismological Society of America}, publisher = {Seismological Society of America}, address = {Albany}, issn = {0037-1106}, doi = {10.1785/0120160215}, pages = {1688 -- 1703}, year = {2017}, abstract = {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.}, language = {en} } @article{KaramzadehToularoudHeimannDahmetal.2020, author = {Karamzadeh Toularoud, Nasim and Heimann, Sebastian and Dahm, Torsten and Kr{\"u}ger, Frank}, title = {Earthquake source arrays}, series = {Geophysical journal international}, volume = {221}, journal = {Geophysical journal international}, number = {1}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggaa002}, pages = {352 -- 370}, year = {2020}, abstract = {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.}, language = {en} } @article{KaramzadehToularoudHeimannDahmetal.2018, author = {Karamzadeh Toularoud, Nasim and Heimann, Sebastian and Dahm, Torsten and Kr{\"u}ger, Frank}, title = {Application based seismological array design by seismicity scenario modelling}, series = {Geophysical journal international}, volume = {216}, journal = {Geophysical journal international}, number = {3}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, pages = {1711 -- 1727}, year = {2018}, abstract = {The design of an array configuration is an important task in array seismology during experiment planning. Often the array response function (ARF), which depends on the relative position of array stations and frequency content of the incoming signals, is used as the array design criterion. In practice, additional constraints and parameters have to be taken into account, for example, land ownership, site-specific noise levels or characteristics of the seismic sources under investigation. In this study, a flexible array design framework is introduced that implements a customizable scenario modelling and optimization scheme by making use of synthetic seismograms. Using synthetic seismograms to evaluate array performance makes it possible to consider additional constraints. We suggest to use synthetic array beamforming as an array design criterion instead of the ARF. The objective function of the optimization scheme is defined according to the monitoring goals, and may consist of a number of subfunctions. The array design framework is exemplified by designing a seven-station small-scale array to monitor earthquake swarm activity in Northwest Bohemia/Vogtland in central Europe. Two subfunctions are introduced to verify the accuracy of horizontal slowness estimation; one to suppress aliasing effects due to possible secondary lobes of synthetic array beamforming calculated in horizontal slowness space and the other to reduce the event's mislocation caused by miscalculation of the horizontal slowness vector. Subsequently, a weighting technique is applied to combine the subfunctions into one single scalar objective function to use in the optimization process.}, language = {en} }