Refine
Year of publication
Language
- English (72) (remove)
Keywords
- Body waves (7)
- Seismologie (6)
- Seismology (6)
- Earthquake source observations (5)
- Erdbeben (5)
- Seismicity and tectonics (5)
- Theoretical seismology (4)
- Array Seismology (3)
- Arrayseismologie (3)
- Crustal structure (3)
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.
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.
The correct estimation of site-specific attenuation is crucial for the assessment of seismic hazard. Downhole instruments provide in this context valuable information to constrain attenuation directly from data. In this study, we apply an interferometric approach to this problem by deconvolving seismic motions recorded at depth with those recorded at the surface. In doing so, incident and surface-reflected waves can be separated. We apply this technique not only to earthquake data but also to recordings of ambient vibrations. We compute the transfer function between incident and surface-reflected waves in order to infer frequency-dependent quality factors for S waves. The method is applied to a 87m deep borehole sensor and a colocated surface instrument situated at a hard-rock site in West Bohemia/Vogtland, Germany. We show that the described method provides comparable attenuation estimates using either earthquake data or ambient noise for frequencies between 5 and 15 Hz. Moreover, a single hour of noise recordings seems to be sufficient to yield stable deconvolution traces and quality factors, thus, offering a fast and easy way to derive attenuation estimates from borehole recordings even in low- to mid-seismicity regions.
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.
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.
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).
Tensile source components of swarm events in West Bohemia in 2000 by considering seismic anisotropy
(2006)
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.
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.
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.
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.
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.
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.
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.
An der Universitä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 Öffentlichkeit via Internet zur Verfügung zu stellen. Es ist vorgesehen, das System in das Deutsch-Indonesische Tsunamifrühwarnsystem (GITEWS) zu integrieren, fü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ömmliche Arrayverfahren in teilweise modifizierter Form. Die Semblance wir als Ähnlichkeitsmaß verwendet, um Seismogramme eines Stationsnetzes zu vergleichen. Im Falle eines Erdbebens ist die Semblance unter Berücksichtigung des Hypozentrums zur Herdzeit und während des Bruchvorgangs deutlich zeitlich und räumlich erhöht und konzentriert. Indem wir die Ergebnisse verschiedener Stationsnetzwerke kombinieren, erreichen wir Unabhängigkeit von der Herdcharakteristik und eine raum-zeitliche Auflö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ösungsmöglichkeiten und vergleichen die Ergebnisse der automatisierten Echtzeitanwendung mit nachträglichen Berechnungen. Weiterhin stellen wir eine Internetseite zur Verfügung, die die Ergebnisse prä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
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.
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.
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.
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.
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
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.