@article{FreybourgerKruegerAchauer1999, author = {Freybourger, Marion and Kr{\"u}ger, Frank and Achauer, Ulrich}, title = {A 22 degree long seismic profile for the study of the top D"}, year = {1999}, 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{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{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} } @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} } @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{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} } @article{BrechnerKlingeKruegeretal.1998, author = {Brechner, Stefan and Klinge, Klaus and Kr{\"u}ger, Frank and Plenefisch, Thomas}, title = {Backazimuthal variations of splitting parameters of teleseismic SKS phases observed at the broadband stations in Germany}, year = {1998}, language = {en} } @article{DonnerKruegerRoessleretal.2014, author = {Donner, Stefanie and Kr{\"u}ger, Frank and Roessler, Dirk and Ghods, Abdolreza}, title = {Combined Inversion of broadband and short-period waveform data for regional moment tensors: A case study in the Alborz Mountains, Iran}, series = {Bulletin of the Seismological Society of America}, volume = {104}, journal = {Bulletin of the Seismological Society of America}, number = {3}, publisher = {Seismological Society of America}, address = {Albany}, issn = {0037-1106}, doi = {10.1785/0120130229}, pages = {1358 -- 1373}, year = {2014}, abstract = {In this study, we suggest a novel approach for the retrieval of regional moment tensors for earthquakes with small to moderate magnitudes. The first modification is the combined inversion of broadband and short-period waveform data. The broadband waveforms are inverted in a frequency range suitable for surface waves, whereas for the short-period data a frequency range suitable for body waves is applied. The second modification is the use of first-motion body-wave polarities to select the most probable solution out of all solutions from inversion. To combine three different criteria for selecting the most probable solution (i.e., residual from inversion, double-couple content of solution, number of nonmatching first-motion body-wave polarities), the L2 norm is applied to the normalized parameters. We chose five earthquakes within the Alborz mountains, Iran, as a case study (3.1 <= M-w <= 4.1). In this area, several factors exacerbate the difficulty of performing inversion for moment tensors, for example, a heterogeneous station network and large azimuthal gaps. We have demonstrated that our approach supplies reliable moment tensors when inversion from broadband data alone fails. In one case, we successfully retrieved a stable solution from short-period waveform data alone. Thus, our approach enables successful determination of seismic moment tensors wherever a sparse network of broadband stations has thus far prevented it.}, language = {en} } @article{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{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{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{CzubaGradMjeldeetal.2011, author = {Czuba, Wojciech and Grad, Marek and Mjelde, Rolf and Guterch, Aleksander and Libak, Audun and Kr{\"u}ger, Frank and Murai, Yoshio and Schweitzer, Johannes}, title = {Continent-ocean-transition across a trans-tensional margin segment: off Bear Island, Barents Sea}, series = {Geophysical journal international}, volume = {184}, journal = {Geophysical journal international}, number = {2}, publisher = {Oxford Univ. Press}, address = {Oxford}, organization = {IPY Project Grp}, issn = {0956-540X}, doi = {10.1111/j.1365-246X.2010.04873.x}, pages = {541 -- 554}, year = {2011}, abstract = {P>A 410 km long Ocean Bottom Seismometer profile spanning from the Bear Island, Barents Sea to oceanic crust formed along the Mohns Ridge has been modelled by use of ray-tracing with regard to observed P-waves. The northeastern part of the model represents typical continental crust, thinned from ca. 30 km thickness beneath the Bear Island to ca. 13 km within the Continent-Ocean-Transition. Between the Hornsund FZ and the Kn circle divide legga Fault, a 3-4 km thick sedimentary basin, dominantly of Permian/Carboniferous age, is modelled beneath the ca. 1.5 km thick layer of volcanics (Vestbakken Volcanic Province). The P-wave velocity in the 3-4 km thick lowermost continental crust is significantly higher than normal (ca. 7.5 km s-1). We interpret this layer as a mixture of mafic intrusions and continental crystalline blocks, dominantly related to the Paleocene-Early Eocene rifting event. The crystalline portion of the crust within the south-western part of the COT consists of a ca. 30 km wide and ca. 6 km thick high-velocity (7.3 km s-1) body. We interpret the body as a ridge of serpentinized peridotites. The magmatic portion of the ocean crust accreted along the Knipovich Ridge from continental break-up at ca. 35 Ma until ca. 20 Ma is 3-5 km thicker than normal. We interpret the increased magmatism as a passive response to the bending of this southernmost part of the Knipovich Ridge. The thickness of the magmatic portion of the crust formed along the Mohns Ridge at ca. 20 Ma decreases to ca. 3 km, which is normal for ultra slow spreading ridges.}, language = {en} } @article{DahmCescaHainzletal.2015, author = {Dahm, Torsten and Cesca, Simone and Hainzl, Sebastian and Braun, Thomas and Kr{\"u}ger, Frank}, title = {Discrimination between induced, triggered, and natural earthquakes close to hydrocarbon reservoirs: A probabilistic approach based on the modeling of depletion-induced stress changes and seismological source parameters}, series = {Journal of geophysical research : Solid earth}, volume = {120}, journal = {Journal of geophysical research : Solid earth}, number = {4}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9313}, doi = {10.1002/2014JB011778}, pages = {2491 -- 2509}, year = {2015}, abstract = {Earthquakes occurring close to hydrocarbon fields under production are often under critical view of being induced or triggered. However, clear and testable rules to discriminate the different events have rarely been developed and tested. The unresolved scientific problem may lead to lengthy public disputes with unpredictable impact on the local acceptance of the exploitation and field operations. We propose a quantitative approach to discriminate induced, triggered, and natural earthquakes, which is based on testable input parameters. Maxima of occurrence probabilities are compared for the cases under question, and a single probability of being triggered or induced is reported. The uncertainties of earthquake location and other input parameters are considered in terms of the integration over probability density functions. The probability that events have been human triggered/induced is derived from the modeling of Coulomb stress changes and a rate and state-dependent seismicity model. In our case a 3-D boundary element method has been adapted for the nuclei of strain approach to estimate the stress changes outside the reservoir, which are related to pore pressure changes in the field formation. The predicted rate of natural earthquakes is either derived from the background seismicity or, in case of rare events, from an estimate of the tectonic stress rate. Instrumentally derived seismological information on the event location, source mechanism, and the size of the rupture plane is of advantage for the method. If the rupture plane has been estimated, the discrimination between induced or only triggered events is theoretically possible if probability functions are convolved with a rupture fault filter. We apply the approach to three recent main shock events: (1) the M-w 4.3 Ekofisk 2001, North Sea, earthquake close to the Ekofisk oil field; (2) the M-w 4.4 Rotenburg 2004, Northern Germany, earthquake in the vicinity of the Sohlingen gas field; and (3) the M-w 6.1 Emilia 2012, Northern Italy, earthquake in the vicinity of a hydrocarbon reservoir. The three test cases cover the complete range of possible causes: clearly human induced, not even human triggered, and a third case in between both extremes.}, 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{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} } @misc{DonnerStreckerRoessleretal.2009, author = {Donner, Stefanie and Strecker, Manfred and R{\"o}ßler, Dirk and Ghods, Abdolreza and Kr{\"u}ger, Frank and Landgraf, Angela and Ballato, Paolo}, title = {Earthquake source models for earthquakes in Northern Iran}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-32581}, year = {2009}, abstract = {The complex system of strike-slip and thrust faults in the Alborz Mountains, Northern Iran, are not well understood yet. Mainly structural and geomorphic data are available so far. As a more extensive base for seismotectonic studies and seismic hazard analysis we plan to do a comprehensive seismic moment tensor study also from smaller magnitudes (M < 4.5) by developing a new algorithm. Here, we present first preliminary results.}, 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{HaendelOhrnbergerKrueger2016, author = {H{\"a}ndel, Annabel and Ohrnberger, Matthias and Kr{\"u}ger, Frank}, title = {Extracting near-surface Q(L) between 1-4 Hz from higher-order noise correlations in the Euroseistest area, Greece}, series = {Geophysical journal international}, volume = {207}, journal = {Geophysical journal international}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggw295}, pages = {655 -- 666}, year = {2016}, abstract = {Knowledge of the quality factor of near-surface materials is of fundamental interest in various applications. Attenuation can be very strong close to the surface and thus needs to be properly assessed. In recent years, several researchers have studied the retrieval of attenuation coefficients from the cross correlation of ambient seismic noise. Yet, the determination of exact amplitude information from noise-correlation functions is, in contrast to the extraction of traveltimes, not trivial. Most of the studies estimated attenuation coefficients on the regional scale and within the microseism band. In this paper, we investigate the possibility to derive attenuation coefficients from seismic noise at much shallower depths and higher frequencies (> 1 Hz). The Euroseistest area in northern Greece offers ideal conditions to study quality factor retrieval from ambient noise for different rock types. Correlations are computed between the stations of a small scale array experiment (station spacings < 2 km) that was carried out in the Euroseistest area in 2011. We employ the correlation of the coda of the correlation (C-3) method instead of simple cross correlations to mitigate the effect of uneven noise source distributions on the correlation amplitude. Transient removal and temporal flattening are applied instead of 1-bit normalization in order to retain relative amplitudes. The C-3 method leads to improved correlation results (higher signal-to-noise ratio and improved time symmetry) compared to simple cross correlations. The C-3 functions are rotated from the ZNE to the ZRT system and we focus on Love wave arrivals on the transverse component and on Love wave quality factors Q(L). The analysis is performed for selected stations being either situated on soft soil or on weathered rock. Phase slowness is extracted using a slant-stack method. Attenuation parameters are inferred by inspecting the relative amplitude decay of Love waves with increasing interstation distance. We observe that the attenuation coefficient gamma and Q(L) can be reliably extracted for stations situated on soft soil whereas the derivation of attenuation parameters is more problematic for stations that are located on weathered rock. The results are in acceptable conformance with theoretical Love wave attenuation curves that were computed using 1-D shear wave velocity and quality factor profiles from the Euroseistest area.}, language = {en} }