@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{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}
}