TY - JOUR A1 - Blanke, Aglaja A1 - Kwiatek, Grzegorz A1 - Goebel, Thomas H. W. A1 - Bohnhoff, Marco A1 - Dresen, Georg T1 - Stress drop-magnitude dependence of acoustic emissions during laboratory stick-slip JF - Geophysical journal international / the Royal Astronomical Society, the Deutsche Geophysikalische Gesellschaft and the European Geophysical Society N2 - Earthquake source parameters such as seismic stress drop and corner frequency are observed to vary widely, leading to persistent discussion on potential scaling of stress drop and event size. Physical mechanisms that govern stress drop variations arc difficult to evaluate in nature and are more readily studied in controlled laboratory experiments. We perform two stick-slip experiments on fractured (rough) and cut (smooth) Westerly granite samples to explore fault roughness effects on acoustic emission (AE) source parameters. We separate large stick-slip events that generally saturate the seismic recording system from populations of smaller AE events which are sensitive to fault stresses prior to slip. AE event populations show many similarities to natural seismicity and may be interpreted as laboratory equivalent of natural microseismic events. We then compare the temporal evolution of mechanical data such as measured stress release during slip to temporal changes in stress drops derived from Alis using the spectral ratio technique. We report on two primary observations: (1) In contrast to most case studies for natural earthquakes, we observe a strong increase in seismic stress drop with AE size. (2) The scaling of stress drop with magnitude is governed by fault roughness, whereby the rough fault shows a more rapid increase of the stress drop magnitude relation with progressing large stick-slip events than the smooth fault. The overall range of AE sizes on the rough surface is influenced by both the average grain size and the width of the fault core. The magnitudes of the smallest AE events on smooth faults may also be governed by grain size. However, AEs significantly grow beyond peak roughness and the width of the fault core. Our laboratory tests highlight that source parameters vary substantially in the presence of fault zone heterogeneity (i.e. roughness and narrow grain size distribution), which may affect seismic energy partitioning and static stress drops of small and large AE events. KW - Acoustic properties KW - Body waves KW - Earthquake dynamics KW - Earthquake source observations KW - Dynamics and mechanics of faulting Y1 - 2021 U6 - https://doi.org/10.1093/gji/ggaa524 SN - 0956-540X SN - 1365-246X VL - 224 IS - 2 SP - 1372 EP - 1381 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Weidle, Christian A1 - Wiesenberg, Lars A1 - El-Sharkawy, Amr A1 - Krüger, Frank A1 - Scharf, Andreas A1 - Agard, Philippe A1 - Meier, Thomas T1 - A 3-D crustal shear wave velocity model and Moho map below the Semail Ophiolite, eastern Arabia JF - Geophysical journal international N2 - 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. KW - Composition and structure of the continental crust KW - Asia KW - Body waves KW - Seismic anisotropy KW - Seismic tomography KW - Surface waves and free oscillations Y1 - 2022 U6 - https://doi.org/10.1093/gji/ggac223 SN - 0956-540X SN - 1365-246X VL - 231 IS - 2 SP - 817 EP - 834 PB - Oxford University Press CY - Oxford ER - TY - JOUR A1 - Hannemann, Katrin A1 - Eulenfeld, Tom A1 - Krüger, Frank A1 - Dahm, Torsten T1 - Seismic scattering and absorption of oceanic lithospheric S waves in the Eastern North Atlantic JF - Geophysical journal international N2 - The scattering and absorption of high-frequency seismic waves in the oceanic lithosphere is to date only poorly constrained by observations. Such estimates would not only improve our understanding of the propagation of seismic waves, but also unravel the small-scale nature of the lithosphere and its variability. Our study benefits from two exceptional situations: (1) we deployed over 10 months a mid-aperture seismological array in the central part of the Eastern North Atlantic in 5 km water depth and (2) we could observe in total 340 high-frequency (up to 30 Hz) Po and So arrivals with tens to hundreds of seconds long seismic coda from local and regional earthquakes in a wide range of backazimuths and epicentral distances up to 850 km with a travel path in the oceanic lithosphere. Moreover, the array was located about 100 km north of the Gloria fault, defining the plate boundary between the Eurasian and African plates at this location which also allows an investigation of the influence of an abrupt change in lithospheric age (20 Ma in this case) on seismic waves. The waves travel with velocities indicating upper-mantle material. We use So waves and their coda of pre-selected earthquakes to estimate frequency-dependent seismic scattering and intrinsic attenuation parameters. The estimated scattering attenuation coefficients are between 10(-4) and 4 x 10(-5) m(-1) and are typical for the lithosphere or the upper mantle. Furthermore, the total quality factors for So waves below 5 Hz are between 20 and 500 and are well below estimates from previous modelling for observations in the Pacific Ocean. This implies that the Atlantic Ocean is more attenuative for So waves compared to the Pacific Ocean, which is inline with the expected behaviour for the lithospheric structures resulting from the slower spreading rates in the Atlantic Ocean. The results for the analysed events indicate that for frequencies above 3 Hz, intrinsic attenuation is equal to or slightly stronger than scattering attenuation and that the So-wave coda is weakly influenced by the oceanic crust. Both observations are in agreement with the proposed propagation mechanism of scattering in the oceanic mantle lithosphere. Furthermore, we observe an age dependence which shows that an increase in lithospheric age is associated with a decrease in attenuation. However, we also observe a trade-off of this age-dependent effect with either a change in lithospheric thickness or thermal variations, for example due to small-scale upwellings in the upper mantle in the southeast close to Madeira and the Canaries. Moreover, the influence of the nearby Gloria fault is visible in a reduction of the intrinsic attenuation below 3 Hz for estimates across the fault. This is the first study to estimate seismic scattering and absorption parameters of So waves for an area with several hundreds of kilometres radius centred in the Eastern North Atlantic and using them to characterize the nature of the oceanic lithosphere. KW - Body waves KW - Seismic attenuation KW - Wave scattering and diffraction Y1 - 2021 U6 - https://doi.org/10.1093/gji/ggab493 SN - 0956-540X SN - 1365-246X VL - 229 IS - 2 SP - 948 EP - 961 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Letort, Jean A1 - Retailleau, Lise A1 - Boue, Pierre A1 - Radiguet, Mathilde A1 - Gardonio, Blandine A1 - Cotton, Fabrice A1 - Campillo, Michel T1 - Lateral variations of the Guerrero-Oaxaca subduction zone (Mexico) derived from weak seismicity (M(b)3.5+) detected on a single array at teleseismic distance JF - Geophysical journal international N2 - Detections of pP and sP phase arrivals (the so-called depth phases) at teleseismic distance provide one of the best ways to estimate earthquake focal depth, as the P-pP and the P-sP delays are strongly dependent on the depth. Based on a new processing workflow and using a single seismic array at teleseismic distance, we can estimate the depth of clusters of small events down to magnitude M-b 3.5. Our method provides a direct view of the relative variations of the seismicity depth from an active area. This study focuses on the application of this new methodology to study the lateral variations of the Guerrero subduction zone (Mexico) using the Eielson seismic array in Alaska (USA). After denoising the signals, 1232 M-b 3.5 + events were detected, with clear P, pP, sP and PcP arrivals. A high-resolution view of the lateral variations of the depth of the seismicity of the Guerero-Oaxaca area is thus obtained. The seismicity is shown to be mainly clustered along the interface, coherently following the geometry of the plate as constrained by the receiver-function analysis along the Meso America Subduction Experiment profile. From this study, the hypothesis of tears on the western part of Guerrero and the eastern part of Oaxaca are strongly confirmed by dramatic lateral changes in the depth of the earthquake clusters. The presence of these two tears might explain the observed lateral variations in seismicity, which is correlated with the boundaries of the slow slip events. KW - North America KW - Time-series analysis KW - Body waves KW - Earthquake source observations KW - Seismicity and tectonics Y1 - 2018 U6 - https://doi.org/10.1093/gji/ggy035 SN - 0956-540X SN - 1365-246X VL - 213 IS - 2 SP - 1002 EP - 1012 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Hendriyana, Andri A1 - Bauer, Klaus A1 - Muksin, Umar A1 - Weber, Michael T1 - AIC-based diffraction stacking for local earthquake locations at the Sumatran Fault (Indonesia) JF - Geophysical journal international N2 - We present a new workflow for the localization of seismic events which is based on a diffraction stacking approach. In order to address the effects from complex source radiation patterns, we suggest to compute diffraction stacking from a characteristic function (CF) instead of stacking the original waveform data. A new CF, which is called in the following mAIC (modified from Akaike Information Criterion) is proposed. We demonstrate that both P- and S-wave onsets can be detected accurately. To avoid cross-talk between P and S waves due to inaccurate velocity models, we separate the P and S waves from the mAIC function by making use of polarization attributes. Then, the final image function is represented by the largest eigenvalue as a result of the covariance analysis between P-and S-image functions. Results from synthetic experiments show that the proposed diffraction stacking provides reliable results. The workflow of the diffraction stacking method was finally applied to local earthquake data from Sumatra, Indonesia. Recordings from a temporary network of 42 stations deployed for nine months around the Tarutung pull-apart basin were analysed. The seismic event locations resulting from the diffraction stacking method align along a segment of the Sumatran Fault. A more complex distribution of seismicity is imaged within and around the Tarutung basin. Two lineaments striking N-S were found in the centre of the Tarutung basin which support independent results from structural geology. KW - Time-series analysis KW - Body waves KW - Computational seismology KW - Earthquake source observations KW - Seismicity and tectonics Y1 - 2018 U6 - https://doi.org/10.1093/gji/ggy045 SN - 0956-540X SN - 1365-246X VL - 213 IS - 2 SP - 952 EP - 962 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Kulikova, Galina A1 - Schurr, Bernd A1 - Krüger, Frank A1 - Brzoska, Elisabeth A1 - Heimann, Sebastian T1 - Source parameters of the Sarez-Pamir earthquake of 1911 February 18 JF - Geophysical journal international N2 - The Ms ∼ 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 (∼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. KW - Earthquake source observations KW - Seismicity and tectonics KW - Body waves KW - Theoretical seismology Y1 - 2016 U6 - https://doi.org/10.1093/gji/ggw069 SN - 0956-540X SN - 1365-246X VL - 205 SP - 1086 EP - 1098 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Hannemann, Katrin A1 - Krüger, Frank A1 - Dahm, Torsten A1 - Lange, Dietrich T1 - Oceanic lithospheric S-wave velocities from the analysis of P-wave polarization at the ocean floor JF - Geophysical journal international N2 - Our knowledge of the absolute S-wave velocities of the oceanic lithosphere is mainly based on global surface wave tomography, local active seismic or compliance measurements using oceanic infragravity waves. The results of tomography give a rather smooth picture of the actual S-wave velocity structure and local measurements have limitations regarding the range of elastic parameters or the geometry of the measurement. Here, we use the P-wave polarization (apparent P-wave incidence angle) of teleseismic events to investigate the S-wave velocity structure of the oceanic crust and the upper tens of kilometres of the mantle beneath single stations. In this study, we present an up to our knowledge new relation of the apparent P-wave incidence angle at the ocean bottom dependent on the half-space S-wave velocity. We analyse the angle in different period ranges at ocean bottom stations (OBSs) to derive apparent S-wave velocity profiles. These profiles are dependent on the S-wave velocity as well as on the thickness of the layers in the subsurface. Consequently, their interpretation results in a set of equally valid models. We analyse the apparent P-wave incidence angles of an OBS data set which was collected in the Eastern Mid Atlantic. We are able to determine reasonable S-wave-velocity-depth models by a three-step quantitative modelling after a manual data quality control, although layer resonance sometimes influences the estimated apparent S-wave velocities. The apparent S-wave velocity profiles are well explained by an oceanic PREM model in which the upper part is replaced by four layers consisting of a water column, a sediment, a crust and a layer representing the uppermost mantle. The obtained sediment has a thickness between 0.3 and 0.9 km with S-wave velocities between 0.7 and 1.4 km s(-1). The estimated total crustal thickness varies between 4 and 10 km with S-wave velocities between 3.5 and 4.3 km s(-1). We find a slight increase of the total crustal thickness from similar to 5 to similar to 8 km towards the South in the direction of a major plate boundary, the Gloria Fault. The observed crustal thickening can be related with the known dominant compression in the vicinity of the fault. Furthermore, the resulting mantle S-wave velocities decrease from values around 5.5 to 4.5 km s(-1) towards the fault. This decrease is probably caused by serpentinization and indicates that the oceanic transform fault affects a broad region in the uppermost mantle. Conclusively, the presented method is useful for the estimation of the local S-wave velocity structure beneath ocean bottom seismic stations. It is easy to implement and consists of two main steps: (1) measurement of apparent P-wave incidence angles in different period ranges for real and synthetic data, and (2) comparison of the determined apparent S-wave velocities for real and synthetic data to estimate S-wave velocity-depth models. KW - Time-series analysis KW - Body waves KW - Theoretical seismology KW - Oceanic transform and fracture zone processes Y1 - 2016 U6 - https://doi.org/10.1093/gji/ggw342 SN - 0956-540X SN - 1365-246X VL - 207 SP - 1796 EP - 1817 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Cristiano, L. A1 - Meier, T. A1 - Krüger, F. A1 - Keers, H. A1 - Weidle, C. T1 - Teleseismic P-wave polarization analysis at the Grafenberg array JF - Geophysical journal international N2 - P-wave polarization at the Grafenberg array (GRF) in southern Germany is analysed in terms of azimuthal deviations and deviations in the vertical polarization using 20 yr of broad-band recordings. An automated procedure for estimating P-wave polarization parameters is suggested, based on the definition of a characteristic function, which evaluates the polarization angles and their time variability as well as the amplitude, linearity and the signal-to-noise ratio of the P wave. P-wave polarization at the GRF array is shown to depend mainly on frequency and backazimuth and only slightly on epicentral distance indicating depth-dependent local anisotropy and lateral heterogeneity. A harmonic analysis is applied to the azimuthal anomalies to analyse their periodicity as a function of backazimuth. The dominant periods are 180A degrees A and 360A degrees. At low frequencies, between 0.03 and 0.1 Hz, the observed fast directions of azimuthal anisotropy inferred from the 180A degrees A periodicity are similar across the array. The average fast direction of azimuthal anisotropy at these frequencies is N20A degrees E with an uncertainty of about 8A degrees A and is consistent with fast directions of Pn-wave propagation. Lateral velocity gradients determined for the low-frequency band are compatible with the Moho topography of the area. A more complex pattern in the horizontal fast axis orientation beneath the GRF array is observed in the high-frequency band between 0.1 and 0.5 Hz, and is attributed to anisotropy in the upper crust. A remarkable rotation of the horizontal fast axis orientation across the suture between the geological units Moldanubicum and Saxothuringicum is observed. In contrast, the 360A degrees A periodicity at high frequencies is rather consistent across the array and may either point to lower velocities in the upper crust towards the Bohemian Massif and/or to anisotropy dipping predominantly in the NE-SW direction. Altogether, P-wave polarization analysis indicates the presence of layered lithospheric anisotropy in the area of the GRF array. Seismic anisotropy is more variable in the brittle upper crust compared to lower crustal and subcrustal depths. KW - Body waves KW - Seismic anisotropy KW - Wave propagation Y1 - 2016 U6 - https://doi.org/10.1093/gji/ggw339 SN - 0956-540X SN - 1365-246X VL - 207 SP - 1456 EP - 1471 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Kulikova, Galina A1 - Krüger, Frank T1 - Source process of the 1911 M8.0 Chon-Kemin earthquake: investigation results by analogue seismic records JF - Geophysical journal international N2 - 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. KW - Earthquake source observations KW - Seismicity and tectonics KW - Body waves KW - Theoretical seismology Y1 - 2015 U6 - https://doi.org/10.1093/gji/ggv091 SN - 0956-540X SN - 1365-246X VL - 201 IS - 3 SP - 1891 EP - 1911 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Knapmeyer-Endrun, Brigitte A1 - Krüger, Frank T1 - Moho depth across the Trans-European Suture Zone from P- and S-receiver functions JF - Geophysical journal international N2 - 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. KW - Body waves KW - Cratons KW - Crustal structure KW - Europe Y1 - 2014 U6 - https://doi.org/10.1093/gji/ggu035 SN - 0956-540X SN - 1365-246X VL - 197 IS - 2 SP - 1048 EP - 1075 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Alinaghi, Alireza A1 - Kruger, Frank T1 - Seismic array analysis and redetermination of depths of earthquakes in Tien-Shan: implications for strength of the crust and lithosphere JF - Geophysical journal international N2 - We have redetermined focal depths of moderate and major earthquakes with reported lower-crust and upper-mantle depths that have occurred in Tien-Shan, since the availability of broad-band array data. Records of earthquakes at global arrays have been used for identification and modelling of depth phases in order to make accurate estimation of focal depths. Our results show that half of the purportedly deep earthquakes are indeed originating from depths attributable to middle-crust and lower-crust regions. Also one exceptional event in the northern foreland of Tien-Shan in Junggar Basin is located in the upper mantle at the depth of 64 km. Such unusually deep earthquakes for intraplate continental tectonic domain are all located at the margin of Tien-Shan with its adjacent stable blocks and at least some of them have occurred where the brittle behaviour of continental rocks is not highly expected. The reverse mechanisms of all these earthquakes and their proximity to formerly subducting and later colliding and underplating stable blocks and their interactions with overlying Tien-Shan are clues to explain this extremity. KW - Earthquake source observations KW - Seismicity and tectonics KW - Body waves Y1 - 2014 U6 - https://doi.org/10.1093/gji/ggu141 SN - 0956-540X SN - 1365-246X VL - 198 IS - 2 SP - 1111 EP - 1129 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Di Giacomo, Domenico A1 - Bindi, Dino A1 - Parolai, Stefano A1 - Oth, Adrien T1 - Residual analysis of teleseismic P-wave energy magnitude estimates: inter- and intrastation variability JF - Geophysical journal international N2 - P>Computing the magnitude of an earthquake requires correcting for the propagation effects from the source to the receivers. This is often accomplished by performing numerical simulations using a suitable Earth model. In this work, the energy magnitude M(e) is considered and its determination is performed using theoretical spectral amplitude decay functions over teleseismic distances based on the global Earth model AK135Q. Since the high frequency part (above the corner frequency) of the source spectrum has to be considered in computing M(e), the influence of propagation and site effects may not be negligible and they could bias the single station M(e) estimations. Therefore, in this study we assess the inter- and intrastation distributions of errors by considering the M(e) residuals computed for a large data set of earthquakes recorded at teleseismic distances by seismic stations deployed worldwide. To separate the inter- and intrastation contribution of errors, we apply a maximum likelihood approach to the M(e) residuals. We show that the interstation errors (describing a sort of site effect for a station) are within +/- 0.2 magnitude units for most stations and their spatial distribution reflects the expected lateral variation affecting the velocity and attenuation of the Earth's structure in the uppermost layers, not accounted for by the 1-D AK135Q model. The variance of the intrastation error distribution (describing the record-to-record component of variability) is larger than the interstation one (0.240 against 0.159), and the spatial distribution of the errors is not random but shows specific patterns depending on the source-to-station paths. The set of coefficients empirically determined may be used in the future to account for the heterogeneities of the real Earth not considered in the theoretical calculations of the spectral amplitude decay functions used to correct the recorded data for propagation effects. KW - Time series analysis KW - Earthquake source observations KW - Body waves KW - Site effects KW - Wave propagation Y1 - 2011 U6 - https://doi.org/10.1111/j.1365-246X.2011.05019.x SN - 0956-540X VL - 185 IS - 3 SP - 1444 EP - 1454 PB - Wiley-Blackwell CY - Malden ER - TY - JOUR A1 - Muksin, Umar A1 - Haberland, Christian A1 - Bauer, Klaus A1 - Weber, Michael H. T1 - Three-dimensional upper crustal structure of the geothermal system in Tarutung (North Sumatra, Indonesia) revealed by seismic attenuation tomography JF - Geophysical journal international N2 - The geothermal potential in Tarutung is controlled by both the Sumatra Fault system and young arc volcanism. In this study we use the spatial distribution of seismic attenuation, calculated from local earthquake recordings, to image the 3-D seismic attenuation of the area and relate it with the temperature anomalies and the fluid distribution of the subsurface. A temporary seismic network of 42 stations was deployed around Tarutung and Sarulla (south of Tarutung) for a period of 10 months starting in 2011 May. Within this period, the network recorded 2586 local events. A high-quality subset of 229 events recorded by at least 10 stations was used for the attenuation inversion (tomography). Path-average attenuation (t(p)*) was calculated by using a spectral inversion method. The spread function, the contour lines of the model resolution matrix and the recovery test results show that our 3-D attenuation model (Q(p)) has good resolution around the Tarutung Basin and along the Sarulla graben. High attenuation (low Q(p)) related to the geothermal system is found in the northeast of the Tarutung Basin suggesting fluid pathways from below the Sumatra Fault. The upper part of the studied geothermal system in the Tarutung district seems to be mainly controlled by the fault structure rather than by magmatic activities. In the southwest of the Tarutung Basin, the high attenuation zone is associated with the Martimbang volcano. In the Sarulla region, a low-Q(p) anomaly is found along the graben within the vicinity of the Hopong caldera. KW - Seismicity and tectonics KW - Body waves KW - Seismic attenuation KW - Seismic tomography Y1 - 2013 U6 - https://doi.org/10.1093/gji/ggt383 SN - 0956-540X SN - 1365-246X VL - 195 IS - 3 SP - 2037 EP - 2049 PB - Oxford Univ. Press CY - Oxford ER -