Refine
Year of publication
Document Type
- Article (61)
- Other (9)
- Postprint (2)
- Monograph/Edited Volume (1)
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)
Records from ocean bottom seismometers (OBSs) are highly contaminated by noise, which is much stronger
compared to data from most land stations, especially on the horizontal components. As a consequence, the high energy of the oceanic noise at frequencies below 1 Hz considerably complicates the analysis of the teleseismic earthquake signals recorded by OBSs.
Previous studies suggested different approaches to remove low-frequency noises from OBS recordings but mainly focused on the vertical component. The records of horizontal components, which are crucial for the application of many methods in passive seismological analysis of body and surface waves, could not be much improved in the teleseismic frequency band. Here we introduce a noise reduction method, which is derived from the harmonic–percussive separation algorithms used in Zali et al. (2021), in order to separate long-lasting narrowband signals from broadband transients in the OBS signal. This leads to significant noise reduction of OBS records on both the vertical and horizontal components and increases the earthquake signal-to-noise ratio (SNR) without distortion of the broadband earthquake waveforms. This is demonstrated through tests with synthetic data. Both SNR and cross-correlation coefficients showed significant improvements for different realistic noise realizations. The application of denoised signals in surface wave analysis and receiver functions is discussed through tests with synthetic and real data.
Records from ocean bottom seismometers (OBSs) are highly contaminated by noise, which is much stronger compared to data from most land stations, especially on the horizontal components. As a consequence, the high energy of the oceanic noise at frequencies below 1 Hz considerably complicates the analysis of the teleseismic earthquake signals recorded by OBSs.
Previous studies suggested different approaches to remove low-frequency noises from OBS recordings but mainly focused on the vertical component. The records of horizontal components, which are crucial for the application of many methods in passive seismological analysis of body and surface waves, could not be much improved in the teleseismic frequency band. Here we introduce a noise reduction method, which is derived from the harmonic–percussive separation algorithms used in Zali et al. (2021), in order to separate long-lasting narrowband signals from broadband transients in the OBS signal. This leads to significant noise reduction of OBS records on both the vertical and horizontal components and increases the earthquake signal-to-noise ratio (SNR) without distortion of the broadband earthquake waveforms. This is demonstrated through tests with synthetic data. Both SNR and cross-correlation coefficients showed significant improvements for different realistic noise realizations. The application of denoised signals in surface wave analysis and receiver functions is discussed through tests with synthetic and real data.
A 3-D crustal shear wave velocity model and Moho map below the Semail Ophiolite, eastern Arabia
(2022)
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.
Studies of seismic tomography have been highly successful at imaging the deep structure of subduction zones. In a study complementary to these tomographic studies, we use array seismology and reflected waves to image a stagnant slab in the mantle transition zone. Using P and S (SH) waves we find a steeply dipping reflector centred at ca. 400 km depth and ca. 550 km west of the present Mariana subduction zone (at 20N, 140E). The discovery of this anomaly in tomography and independently in array seismology (this paper) helps in understanding the evolution of the Mariana margin. The reflector/stagnant slab may be the remains of the hypothetical North New Guinea Plate, which was theorized to have subducted ca. 50 Ma.
Earthquakes often rupture across more than one fault segment. If such rupture segmentation occurs on a significant scale, a simple point-source or one-fault model may not represent the rupture process well. As a consequence earthquake characteristics inferred, based on one-source assumptions, may become systematically wrong. This might have effects on follow-up analyses, for example regional stress field inversions and seismic hazard assessments. While rupture segmentation is evident for most M-w > 7 earthquakes, also smaller ones with 5.5 < M-w < 7 can be segmented. We investigate the sensitivity of globally available data sets to rupture segmentation and their resolution to reliably estimate the mechanisms in presence of segmentation. We focus on the sensitivity of InSAR (Interferometric Synthetic Aperture Radar) data in the static near-field and seismic waveforms in the far-field of the rupture and carry out non-linear and Bayesian optimizations of single-source and two-sources kinematic models (double-couple point sources and finite, rectangular sources) using InSAR and teleseismic waveforms separately. Our case studies comprises of four M-w 6-7 earthquakes: the 2009 L'Aquila and 2016 Amatrice (Italy) and the 2005 and 2008 Zhongba (Tibet) earthquakes. We contrast the data misfits of different source complexity by using the Akaike informational criterion (AIC). We find that the AIC method is well suited for data-driven inferences on significant rupture segmentation for the given data sets. This is based on our observation that an AIC-stated significant improvement of data fit for two-segment models over one-segment models correlates with significantly different mechanisms of the two source segments and their average compared to the single-segment mechanism. We attribute these modelled differences to a sufficient sensitivity of the data to resolve rupture segmentation. Our results show that near-field data are generally more sensitive to rupture segmentation of shallow earthquakes than far-field data but that also teleseismic data can resolve rupture segmentation in the studied magnitude range. We further conclude that a significant difference in the modelled source mechanisms for different segmentations shows that an appropriate choice of model segmentation matters for a robust estimation of source mechanisms. It reduces systematic biases and trade-off and thereby improves the knowledge on the rupture. Our study presents a strategy and method to detect significant rupture segmentation such that an appropriate model complexity can be used in the source mechanism inference. A similar, systematic investigation of earthquakes in the range of M-w 5.5-7 could provide important hazard-relevant statistics on rupture segmentation. In these cases single-source models introduce a systematic bias. Consideration of rupture segmentation therefore matters for a robust estimation of source mechanisms of the studied earthquakes.
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
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.
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
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.
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.
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.
Fluid intelligence is the ability to think flexibly and to understand abstract relations. People with high fluid intelligence (hi-fluIQ) perform better in analogical reasoning tasks than people with average fluid intelligence (ave-fluIQ). Although previous neuroimaging studies reported involvement of parietal and frontal brain regions in geometric analogical reasoning (which is a prototypical task for fluid intelligence), however, neuroimaging findings on geometric analogical reasoning in hi-fluIQ are sparse. Furthermore, evidence on the relation between brain activation and intelligence while solving cognitive tasks is contradictory. The present study was designed to elucidate the cerebral correlates of geometric analogical reasoning in a sample of hi-fluIQ and ave-fluIQ high school students. We employed a geometric analogical reasoning task with graded levels of task difficulty and confirmed the involvement of the parieto-frontal network in solving this task. In addition to characterizing the brain regions involved in geometric analogical reasoning in hi-fluIQ and ave-fluIQ, we found that blood oxygenation level dependency (BOLD) signal changes were greater for hi-fluIQ than for ave-fluIQ in parietal brain regions. However, ave-fluIQ showed greater BOLD signal changes in the anterior cingulate cortex and medial frontal gyrus than hi-fluIQ. Thus, we showed that a similar network of brain regions is involved in geometric analogical reasoning in both groups. Interestingly, the relation between brain activation and intelligence is not mono-directional, but rather, it is specific for each brain region. The negative brain activation-intelligence relationship in frontal brain regions in hi-fluIQ goes along with a better behavioral performance and reflects a lower demand for executive monitoring compared to ave-fluIQ individuals. In conclusion, our data indicate that flexibly modulating the extent of regional cerebral activity is characteristic for fluid intelligence.
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.
We present the first image of the Madeira upper crustal structure, using ambient seismic noise tomography. 16 months of ambient noise, recorded in a dense network of 26 seismometers deployed across Madeira, allowed reconstructing Rayleigh wave Green's functions between receivers. Dispersion analysis was performed in the short period band from 1.0 to 4.0 s. Group velocity measurements were regionalized to obtain 20 tomographic images, with a lateral resolution of 2.0 km in central Madeira. Afterwards, the dispersion curves, extracted from each cell of the 2D group velocity maps, were inverted as a function of depth to obtain a 3D shear wave velocity model of the upper crust, from the surface to a depth of 2.0 km. The obtained 3D velocity model reveals features throughout the island that correlates well with surface geology and island evolution. (C) 2015 Elsevier B.V. All rights reserved.
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.