ALBERT

All Library Books, journals and Electronic Records Telegrafenberg

Your email was sent successfully. Check your inbox.

An error occurred while sending the email. Please try again.

Proceed reservation?

Export
  • 1
    facet.materialart.
    Unknown
    Oxford University Press
    In:  New York, Oxford University Press, vol. 4, no. 1, pp. 1-40, (ISBN: 1-4020-1348-5 hb, ISBN: 1-4020-1349-3 pb)
    Publication Date: 1997
    Keywords: Textbook of geology ; Seismology ; Tectonics ; Crustal deformation (cf. Earthquake precursor: deformation or strain) ; Earthquake hazard ; Earthquake risk ; Induced seismicity ; Magnitude ; Maximum likelihood
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 2
    Publication Date: 2015-08-08
    Description: The attenuation and velocity dispersion of sonic waves contain valuable information on the mechanical and hydraulic properties of the probed medium. An inherent complication arising in the interpretation of corresponding measurements is, however, that there are multiple physical mechanisms contributing to the energy dissipation and that the relative importance of the various contributions is difficult to unravel. To address this problem for the practically relevant case of terrestrial alluvial sediments, we analyse the attenuation and velocity dispersion characteristics of broad-band multifrequency sonic logs with dominant source frequencies ranging between 1 and 30 kHz. To adequately compensate for the effects of geometrical spreading, which is critical for reliable attenuation estimates, we simulate our experimental setup using a correspondingly targeted numerical solution of the poroelastic equations. After having applied the thus inferred corrections, the broad-band sonic log data set, in conjunction with a comprehensive suite of complementary logging data, allows for assessing the relative importance of a range of pertinent attenuation mechanisms. In doing so, we focus on the effects of wave-induced fluid flow over a wide range of scales. Our results indicate that the levels of attenuation due to the presence of mesoscopic heterogeneities in unconsolidated clastic sediments fully saturated with water are expected to be largely negligible. Conversely, Monte-Carlo-type inversions indicate that Biot's classical model permits to explain most of the considered data. Refinements with regard to the fitting of the observed attenuation and velocity dispersion characteristics are locally provided by accounting for energy dissipation at the microscopic scale, although the nature of the underlying physical mechanism remains speculative.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 3
    Publication Date: 2015-07-30
    Description: We investigate spatiotemporal variations of the crustal stress field orientation along the rupture zones of the 1999 August Izmit M w 7.4 and November Düzce M w 7.1 earthquakes at the North Anatolian Fault zone (NAFZ) in NW Turkey. Our primary focus is to elaborate on the relation between the state of the crustal stress field and distinct seismotectonic features as well as variations of coseismic slip within the seismogenic layer of the crust. To achieve this, we compile an extensive data base of hypocentres and first-motion polarities including a newly derived local hypocentre catalogue extending from 2 yr prior (1997) to 2 yr after (2001) the Izmit and Düzce main shocks. This combined data set allows studying spatial and temporal variations of stress field orientation along distinct fault segments for the pre- and post-seimic phase of the two large earthquakes in detail. Furthermore, the occurrence of two M  〉 7 earthquakes in rapid succession gives the unique opportunity to analyse the 87-d-long ‘inter-seismic phase’ between them. We use the MOTSI (first MOTion polarity Stress Inversion) procedure directly inverting first-motion polarities to study the stress field evolution of nine distinct segments. In particular, this allows to determine the stress tensor also for the pre- and post-seismic phases when no stable single-event focal mechanisms can be determined. We observe significantly different stress field orientations along the combined 200-km-long rupture in accordance with lateral variations of coseismic slip and seismotectonic setting. Distinct vertical linear segments of the NAFZ show either pure-strike slip behaviour or transtensional and normal faulting if located near pull-apart basins. Pull-apart structures such as the Akyazi and Düzce basins show a predominant normal faulting behaviour along the NAFZ and reflect clearly different characteristic from neighbouring strike-slip segments. Substantial lateral stress field heterogeneity following the two main shocks is observed that declines with time towards the post-seismic period that rather reflects the regional right-lateral strike-slip stress field.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 4
    Publication Date: 2015-08-08
    Description: The Piton de la Fournaise basaltic volcano, on La Réunion Island in the western Indian Ocean, is one of the most active volcanoes in the world. This volcano is classically considered as the surface expression of an upwelling mantle plume and its activity is continuously monitored, providing detailed information on its superficial dynamics and on the edifice structure. Deeper crustal and upper mantle structure under La Réunion Island is surprisingly poorly constrained, motivating this study. We used receiver function techniques to determine a shear wave velocity profile through the crust and uppermost mantle beneath La Réunion, but also at other seismic stations located on the hotspot track, to investigate the plume and lithosphere interaction and its evolution through time. Receiver functions (RFs) were computed at permanent broad-band seismic stations from the GEOSCOPE network (on La Réunion and Rodrigues), at IRIS stations MRIV and DGAR installed on Mauritius and Diego Garcia islands, and at the GEOFON stations KAAM and HMDM on the Maldives. We performed non-linear inversions of RFs through modelling of P -to- S conversions at various crustal and upper mantle interfaces. Joint inversion of RF and surface wave dispersion data suggests a much deeper Mohorovičić discontinuity (Moho) beneath Mauritius (~21 km) compared to La Réunion (~12 km). A magmatic underplated body may be present under La Réunion as a thin layer (≤3 km thick), as suggested by a previous seismic refraction study, and as a much thicker layer beneath other stations located on the hotspot track, suggesting that underplating is an important process resulting from the plume–lithosphere interaction. We find evidence for a strikingly low velocity layer starting at about 33 km depth beneath La Réunion that we interpret as a zone of partial melt beneath the active volcano. We finally observe low velocities below 70 km beneath La Réunion and below 50 km beneath Mauritius that could represent the base of the oceanic lithosphere.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 5
    Publication Date: 2015-07-30
    Description: The scattering of seismic waves travelling in the Earth is not only caused by random velocity heterogeneity but also by surface topography. Both factors are known to strongly affect ground-motion complexity even at relatively short distance from the source. In this study, we simulate ground motion with a 3-D finite-difference wave propagation solver in the 0–5 Hz frequency band using three topography models representative of the Swiss alpine region and realistic heterogeneous media characterized by the Von Karman correlation functions. Subsequently, we analyse and quantify the characteristics of the scattered wavefield in the near-source region. Our study shows that both topography and velocity heterogeneity scattering may excite large coda waves of comparable relative amplitude, especially at around 1 Hz, although large variability in space may occur. Using the single scattering model, we estimate average Q C values in the range 20–30 at 1 Hz, 36–54 at 1.5 Hz and 62–109 at 3 Hz for constant background velocity models with no intrinsic attenuation. In principle, envelopes of topography-scattered seismic waves can be qualitatively predicted by theoretical back-scattering models, while forward- or hybrid-scattering models better reproduce the effects of random velocity heterogeneity on the wavefield. This is because continuous multiple scattering caused by small-scale velocity perturbations leads to more gentle coda decay and envelope broadening, while topography abruptly scatters the wavefield once it impinges the free surface. The large impedance contrast also results in more efficient mode mixing. However, the introduction of realistic low-velocity layers near the free surface increases the complexity of ground motion dramatically and indicates that the role of topography in elastic waves scattering can be relevant especially in proximity of the source. Long-period surface waves can form most of the late coda, especially when intrinsic attenuation is taken into account. Our simulations indicate that both topography and velocity heterogeneity scattering may result in large ground-motion variability, characterized by standard deviation values in the range 0.2–0.5 also at short distance from the source. We conclude that both topography and velocity heterogeneity should be considered to correctly assess the ground-motion variability in earthquake scenario studies even at intermediate frequency.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 6
    Publication Date: 2015-07-30
    Description: In this study, we have systematically investigated the influence of the parameters of the slip-weakening law and the size of nucleation asperity on dynamic rupture of a planar fault in full-space and half-space using the boundary integral equation method, in particular, the occurrence conditions for subshear (or sub-Rayleigh for strike-slip rupture) and supershear ruptures. Besides the well-known rupture styles of subshear (or sub-Rayleigh) and supershear, we defined a new kind of rupture style in this study, termed the ‘self-arresting rupture’, for which the rupture process can be autonomously arrested by itself without any outside interference (e.g. a high strength barrier). Based on the vast number of simulations, we obtained rupture phase diagrams for strike-slip and dip-slip ruptures vertically and obliquely embedded in half-space and full-space with different buried depths. The rupture phase diagram clearly illustrates the occurrence conditions of three kinds of rupture styles and the transitions between them. In full-space, the supershear transition is sensitive with the fault width. Owing to the influence of the free surface, the rupture in half-space becomes much more complicated comparing to the one in full-space. For a strike-slip fault with zero buried depth, all ruptures that occur within the parameter range for sub-Rayleigh ruptures in full-space case become supershear ruptures. This means that as long as a rupture is able to grow incessantly, it will always evolve into a supershear rupture. For dip-slip faults, however, ruptures will always propagate with subshear speed, although slip rate could be almost twice that of a strike-slip fault. Although the influence of the free surface is strong, it is limited to very shallow ruptures (i.e. buried depth 〈1 km). The rupture phase diagram discussed in this study could provide a new insight on earthquake rupture mechanics.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 7
    Publication Date: 2015-08-09
    Description: Mainland Portugal, on the southwestern edge of the European continent, is located directly north of the boundary between the Eurasian and Nubian plates. It lies in a region of slow lithospheric deformation (〈5 mm yr –1 ), which has generated some of the largest earthquakes in Europe, both intraplate (mainland) and interplate (offshore). Some offshore earthquakes are nucleated on old and cold lithospheric mantle, at depths down to 60 km. The seismicity of mainland Portugal and its adjacent offshore has been repeatedly classified as diffuse. In this paper, we analyse the instrumental earthquake catalogue for western Iberia, which covers the period between 1961 and 2013. Between 2010 and 2012, the catalogue was enriched with data from dense broad-band deployments. We show that although the plate boundary south of Portugal is diffuse, in that deformation is accommodated along several distributed faults rather than along one long linear plate boundary, the seismicity itself is not diffuse. Rather, when located using high-quality data, earthquakes collapse into well-defined clusters and lineations. We identify and characterize the most outstanding clusters and lineations of epicentres and correlate them with geophysical and tectonic features (historical seismicity, topography, geologically mapped faults, Moho depth, free-air gravity, magnetic anomalies and geotectonic units). Both onshore and offshore, clusters and lineations of earthquakes are aligned preferentially NNE–SSW and WNW–ESE. Cumulative seismic moment and epicentre density decrease from south to north, with increasing distance from the plate boundary. Only few earthquake lineations coincide with geologically mapped faults. Clusters and lineations that do not match geologically mapped faults may correspond to previously unmapped faults (e.g. blind faults), rheological boundaries or distributed fracturing inside blocks that are more brittle and therefore break more easily than neighbour blocks. The seismicity map of western Iberia presented in this article opens important questions concerning the regional seismotectonics. This work shows that the study of low-magnitude earthquakes using dense seismic deployments is a powerful tool to study lithospheric deformation in slowly deforming regions, such as western Iberia, where high-magnitude earthquakes occur with long recurrence intervals.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 8
    Publication Date: 2015-09-11
    Description: The accurate estimation of dispersion curves has been a key issue for ensuring high quality in geophysical surface wave exploration. Many studies have been carried out on the generation of a high-resolution dispersion image from array measurements. In this study, the sparse signal representation and reconstruction techniques are employed to obtain the high resolution Rayleigh-wave dispersion image from seismic wave data. First, a sparse representation of the seismic wave data is introduced, in which the signal is assumed to be sparse in terms of wave speed. Then, the sparse signal is reconstructed by optimization using l 1 -norm regularization, which gives the signal amplitude spectrum as a function of wave speed. A dispersion image in the f – v domain is generated by arranging the sparse spectra for all frequency slices in the frequency range. Finally, to show the efficiency of the proposed approach, the Surfbar-2 field test data, acquired by B. Luke and colleagues at the University of Nevada Las Vegas, are analysed. By comparing the real-field dispersion image with the results from other methods, the high mode-resolving ability of the proposed approach is demonstrated, particularly for a case with strongly coherent modes.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 9
    Publication Date: 2015-09-26
    Description: We analyse daily cross-correlation computed from continuous records by permanent stations operating in vicinity of the Klyuchevskoy group of volcanoes (Kamchatka). Seismic waves generated by volcanic tremors are clearly seen on the cross-correlations between some pairs of stations as strong signals at frequencies between 0.2 and 2 Hz and with traveltimes typically shorter than those corresponding to interstation propagation. First, we develop a 2-D source-scanning algorithm based on summation of the envelops of cross-correlations to detect seismic tremors and to determine locations from which the strong seismic energy is continuously emitted. In an alternative approach, we explore the distinctive character of the cross-correlation waveforms corresponding to tremors emitted by different volcanoes and develop a phase-matching method for detecting volcanic tremors. Application of these methods allows us to detect and to distinguish tremors generated by the Klyuchevskoy and the Tolbachik, volcanoes and to monitor evolution of their intensity in time.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 10
    Publication Date: 2015-09-26
    Description: In geophysical inversion, inferences of Earth's properties from sparse data involve a trade-off between model complexity and the spatial resolving power. A recent Markov chain Monte Carlo (McMC) technique formalized by Green, the so-called trans-dimensional samplers, allows us to sample between these trade-offs and to parsimoniously arbitrate between the varying complexity of candidate models. Here we present a novel framework using trans-dimensional sampling over tree structures. This new class of McMC sampler can be applied to 1-D, 2-D and 3-D Cartesian and spherical geometries. In addition, the basis functions used by the algorithm are flexible and can include more advanced parametrizations such as wavelets, both in Cartesian and Spherical geometries, to permit Bayesian multiscale analysis. This new framework offers greater flexibility, performance and efficiency for geophysical imaging problems than previous sampling algorithms. Thereby increasing the range of applications and in particular allowing extension to trans-dimensional imaging in 3-D. Examples are presented of its application to 2-D seismic and 3-D teleseismic tomography including estimation of uncertainty.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 11
    Publication Date: 2015-09-26
    Description: Rupture properties, such as rupture direction, length, propagation speed and source duration, provide important insights into earthquake mechanisms. One approach to estimate these properties is to investigate the body-wave duration that depends upon the relative location of the station with respect to the rupture direction. Under the assumption that the propagation is unilateral, the duration can be expressed as a function of the dip and azimuth of the rupture. Examination of duration measurements with respect to both the take-off angle and the azimuth is crucial to obtain robust estimates of rupture parameters, especially for nearly vertical rupture propagation. Moreover, limited data coverage, such as using only teleseismic data, can bias the source duration estimate for dipping ruptures, and this bias can map into estimates of other source properties such as rupture extent and rupture speed. Based upon this framework, we introduce an inversion scheme that uses the duration measurements to obtain four parameters: the source duration, a measure of the rupture extent and speed, and dip and azimuth of the rupture propagation. The method is applied to two deep-focus events in the Sea of Okhotsk region, an M w 7.7 event that occurred on 2012 August 14 and an M w 8.3 event from 2013 May 24. The source durations are 26 ± 1 and 37 ± 1 s, and rupture speeds are 49 ± 4 per cent and 26 ± 3 per cent of shear wave speed for the M w 7.7 and 8.3 events, respectively. The azimuths of the two ruptures are parallel to the trench, but are in opposite directions. The dips of the M w 7.7 and 8.3 events are constrained to be 48° ± 8° downdip and 19° ± 8° updip, respectively. The fit to the data is significantly poorer for the M w 8.3 event than the M w 7.7 event, suggesting that the unilateral rupture may not be a good assumption. The analysis is expanded into a multi-episode model, and a secondary episode is determined for the M w 8.3 event in the southeast direction. The two-episode model gives a better fit to the data than the unilateral model and is compatible with the back-projection analysis, demonstrating that the rupture propagation of the M w 8.3 event is complex.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 12
    Publication Date: 2015-08-29
    Description: We take account of the effect of Earth's surface topography in quasi-dynamic earthquake cycle simulations using a boundary integral equation method. While we have so far assumed a homogeneous elastic half-space medium with a flat free surface, Earth's actual surface topography is complicated. Here, we constructed new slip response functions in half-space with an arbitrarily shaped surface topography in which we used slip response functions in full-space by introducing imaginary free surface cells in addition to embedded fault ones. By comparing analytical slip response functions in the case of a flat surface overlying half-space with the new ones, we developed a computationally efficient method for setting the Earth's surface region, which was divided into cells with the appropriate sizes depending on the fault source cell depth to maintain the computational accuracy. With these new slip response functions, we simulated simple interplate earthquake cycles in the region close to the Japan Trench, off Miyagi, Tohoku, in northeast Japan, which has the amplitude of 7 km in depth. Compared with the case where the flat surface level was set at the trench depth, the slip response functions for the case where actual seafloor topography was used had smaller amplitudes. Hence, the actual topography produces smaller recurrence times for earthquake cycles than that for the flat surface case. These effects of the actual Earth's surface topography mainly come from changes in the distance between the surface and the fault compared with the flat surface case. Changes in the slip response function also represent changes in the fault stiffness of the system. Considering the actual topography of the Earth's surface to be convex upwards as opposed to the flat, the fault stiffness becomes larger compared to the case of the flat Earth's surface. This leads to a change in the frictional instability, and sometime leads to the change in the way of rupture.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 13
    Publication Date: 2015-08-30
    Description: The fate of subducted slabs is enigmatic, yet intriguing. We analyse seismic arrivals at ~20–50 s after the direct P wave in an array in northeast China (NECESSArray) recordings of four deep earthquakes occurring beneath the west-central Pacific subduction zones (from the eastern Indonesia to Tonga region). We employ the array analysing techniques of fourth root vespagram and beam-forming analysis to constrain the slowness and backazimuth of later arrivals. Our analyses reveal that these arrivals have a slightly lower slowness value than the direct P wave and the backazimuth deviates slightly from the great circle direction. Along with calculation of 1-D synthetic seismograms, we conclude that the later arrival is corresponding to an energy of S -to- P converted at a scatterer below the sources. Total five scatterers are detected at depths varying from ~700 to 1110 km in the study region. The past subducted oceanic crust most likely accounts for the seismic scatterers trapped in the mid-mantle beneath the west-central subduction zones. Our observation in turn reflects that oceanic crust at least partly separated from subducted oceanic lithosphere and may be trapped substantially in the mid-mantle surrounding subduction zones, in particular in the western Pacific subduction zones.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 14
    Publication Date: 2015-05-31
    Description: We examine spatio-temporal patterns of microseismicity recorded during one month in an underground mine by addressing three key questions: (1) where does the seismicity occur? (2) Why does it occur in these locations? and (3) what triggers it? To obtain accurate locations, we perform a multiplet analysis and use a modified version of the double-difference (DD) relocation method. This approach leads to highly accurate relative event locations and requires groups of multiplets only. Most of the 281 relocated events are close to the main shaft and tunnels; thus we postulate seismicity is facilitated by stresses associated with the potential for subsidence in addition to the hoop stresses acting on the two vertical shafts. Most events occurred during certain hours of the day and there is a 68 per cent correlation with reported rock removal; therefore, it is likely they were triggered by static and dynamic stress perturbations caused by the transportation of debris along tunnels instead of our initial guess that blasting was the principal causative mechanism. Given that seismicity is present around the main shaft but absent close to the second one, we conclude that for seismicity to occur both a favourable stress state and additional external perturbing forces must exist, thus leading to dynamic event triggering in an initially stable stress situation. This analysis provides more insight into anthropogenic processes that might trigger seismicity, thereby facilitating identification of hazardous and potential damage areas in mine settings.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 15
    Publication Date: 2015-05-31
    Description: We present synthetic tests of 2-D adjoint tomography of surface wave traveltimes obtained by the ambient noise cross-correlation analysis across the Czech Republic. The data coverage may be considered perfect for tomography due to the density of the station distribution. Nevertheless, artefacts in the inferred velocity models arising from the data noise may be still observed when weak regularization (Gaussian smoothing of the misfit gradient) or too many iterations are considered. To examine the effect of the regularization and iteration number on the performance of the tomography in more detail we performed extensive synthetic tests. Instead of the typically used (although criticized) checkerboard test, we propose to carry out the tests with two different target models—simple smooth and complex realistic models. The first test reveals the sensitivity of the result on the data noise, while the second helps to analyse the resolving power of the data set. For various noise and Gaussian smoothing levels, we analysed the convergence towards (or divergence from) the target model with increasing number of iterations. Based on the tests we identified the optimal regularization, which we then employed in the inversion of 16 and 20 s Love-wave group traveltimes.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 16
    Publication Date: 2015-05-31
    Description: Least-squares migration (LSM) is a linearized inversion technique for subsurface reflectivity estimation. Compared to conventional migration algorithms, it can improve spatial resolution significantly with a few iterative calculations. There are three key steps in LSM, (1) calculate data residuals between observed data and demigrated data using the inverted reflectivity model; (2) migrate data residuals to form reflectivity gradient and (3) update reflectivity model using optimization methods. In order to obtain an accurate and high-resolution inversion result, the good estimation of inverse Hessian matrix plays a crucial role. However, due to the large size of Hessian matrix, the inverse matrix calculation is always a tough task. The limited-memory BFGS (L-BFGS) method can evaluate the Hessian matrix indirectly using a limited amount of computer memory which only maintains a history of the past m gradients (often m 〈 10). We combine the L-BFGS method with least-squares pre-stack Kirchhoff depth migration. Then, we validate the introduced approach by the 2-D Marmousi synthetic data set and a 2-D marine data set. The results show that the introduced method can effectively obtain reflectivity model and has a faster convergence rate with two comparison gradient methods. It might be significant for general complex subsurface imaging.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 17
    Publication Date: 2015-05-31
    Description: We present a method for high-resolution imaging of lithospheric structures based on full waveform inversion of teleseismic waveforms. We model the propagation of seismic waves using our recently developed direct solution method/spectral-element method hybrid technique, which allows us to simulate the propagation of short-period teleseismic waves through a regional 3-D model. We implement an iterative quasi-Newton method based upon the L-BFGS algorithm, where the gradient of the misfit function is computed using the adjoint-state method. Compared to gradient or conjugate-gradient methods, the L-BFGS algorithm has a much faster convergence rate. We illustrate the potential of this method on a synthetic test case that consists of a crustal model with a crustal discontinuity at 25 km depth and a sharp Moho jump. This model contains short- and long-wavelength heterogeneities along the lateral and vertical directions. The iterative inversion starts from a smooth 1-D model derived from the IASP91 reference Earth model. We invert both radial and vertical component waveforms, starting from long-period signals filtered at 10 s and gradually decreasing the cut-off period down to 1.25 s. This multiscale algorithm quickly converges towards a model that is very close to the true model, in contrast to inversions involving short-period waveforms only, which always get trapped into a local minimum of the cost function.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 18
    Publication Date: 2015-05-26
    Description: We analyse the physics and geometry of trade-offs between Earth structure and noise sources in interstation noise correlations. Our approach is based on the computation of off-diagonal Hessian elements that describe the extent to which variations in noise sources can compensate for variations in Earth structure without changing the misfit beyond the measurement uncertainty. Despite the fact that all ambient noise inverse problems are special in terms of their receiver configuration and data, some general statements concerning source-structure trade-offs can be made: (i) While source-structure trade-offs may be reduced to some extent by clever measurement design, there are inherent trade-offs that can generally not be avoided. These inherent trade-offs may lead to a mispositioning of structural heterogeneities when the noise source distribution is unknown. (ii) When attenuation is weak, source-structure trade-offs in ambient noise correlations are a global phenomenon, meaning that there is no noise source perturbation that does not trade-off with some Earth structure, and vice versa. (iii) The most significant source-structure trade-offs occur within two elliptically shaped regions connecting a potential noise source perturbation to each one of the receivers. (iv) Far from these elliptical regions, only small-scale structure can trade off against changes in the noise source. (v) While source-structure trade-offs mostly decay with increasing attenuation, they are nearly unaffected by attenuation when the noise source perturbation is located near the receiver-receiver line. This work is intended to contribute to the development of joint source-structure inversions of ambient noise correlations, and in particular to an understanding of the extent to which source-structure trade-offs may be reduced. It furthermore establishes the foundation of future resolution analyses that properly quantify trade-offs between noise sources and Earth structure.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 19
    Publication Date: 2015-05-26
    Description: Reactive flow at depth (either related to underground activities, like enhancement of hydrocarbon recovery and CO 2 storage, or to natural flow like in hydrothermal zones) can alter fractures’ topography, which might in turn change their seismic responses. Depending on the flow and reaction rates, instability of the dissolution front can lead to a wormhole-like pronounced erosion pattern. In a fractal structure of rupture process, we question how the perturbation related to well-spaced long channels alters rupture propagation initiated on a weak plane and eventually the statistical feature of rupture appearance in frequency–magnitude distribution (FMD). Contrary to intuition, a spatially uniform dissolution is not the most remarkable case, since it affects all the events proportionally to their sizes leading to a downward translation of FMD: the slope of FMD ( b -value) remains unchanged. The parameter–space study shows that the increase of b -value (of 0.08) is statistically significant for optimum characteristics of the erosion pattern with spacing to length ratio of the order of ~1/40: large-magnitude events are more significantly affected leading to an imbalanced distribution in the magnitude bins of the FMD. The larger the spacing, the lower the channel's influence. Besides, a spatial analysis shows that the local seismicity anomaly concentrates in a limited zone around the channels: this opens perspective for detecting these eroded regions through high-resolution imaging surveys.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 20
    Publication Date: 2016-07-13
    Description: Earth's free oscillations excited by a mega-thrust earthquake were observed by a continent-scale array of groundwater monitoring sites for the first time. After the occurrence of the 2011 Tohoku M w 9.0 earthquake, water level records at 43 out of 216 wells in the China mainland revealed long-period free oscillation signals. In the time domain, these free oscillations exhibit globe circling Rayleigh surface waves. In some single wells, even the globe-circling Rayleigh wave R7 was visible, which travels three times around the Earth after the first arrival and appears about 10 hr after the earthquake occurrence in the present case. The spectral analysis shows that the principal oscillatory fluctuations seen in the water level records correspond to the spheroidal modes 0 S l ( l  = 2–31 for frequencies between 0.3 and 5.0 mHz) of the Earth's free oscillation. Especially at quiet sites, the spheroidal modes at very low frequencies (〈1.5 mHz) can be identified with high signal-to-noise ratios. Using signal enhancement methods (product spectrum over 43 wells), even the gravest modes of these oscillations can be detected. The results suggest that groundwater level arrays can be considered as a low-cost complementary tool to study the Earth's free oscillations excited by great earthquakes. Additionally, the site-specific aquifer response may provide further insight into local hydrogeological conditions.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 21
    Publication Date: 2016-07-13
    Description: We performed numerical simulations of the 2011 deep-seated Akatani landslide in central Japan to understand the dynamic evolution of friction of the landslide. By comparing the forces obtained from numerical simulation to those resolved from seismic waveform inversion, the coefficient of the friction during sliding was investigated in the range of 0.1–0.4. The simulation assuming standard Coulomb friction shows that the forces obtained by the seismic waveform inversion are well explained using a constant friction of μ = 0.3. A small difference between the residuals of Coulomb simulation and a velocity-dependent simulation suggests that the coefficient of friction over the volume is well constrained as 0.3 most of time during sliding. It suggests the sudden loss of shearing resistance at the onset of sliding, that is, sudden drop of the initial coefficient of friction in our model, which accelerates the deep-seated landslide. Our numerical simulation calibrated by seismic data provides the evolution of dynamic friction with a reasonable resolution in time, which is difficult to obtain from a conventional runout simulation, or seismic waveform inversion alone.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 22
    Publication Date: 2016-07-13
    Description: We introduce a technique to compute exact anelastic sensitivity kernels in the time domain using parsimonious disk storage. The method is based on a reordering of the time loop of time-domain forward/adjoint wave propagation solvers combined with the use of a memory buffer. It avoids instabilities that occur when time-reversing dissipative wave propagation simulations. The total number of required time steps is unchanged compared to usual acoustic or elastic approaches. The cost is reduced by a factor of 4/3 compared to the case in which anelasticity is partially accounted for by accommodating the effects of physical dispersion. We validate our technique by performing a test in which we compare the K α sensitivity kernel to the exact kernel obtained by saving the entire forward calculation. This benchmark confirms that our approach is also exact. We illustrate the importance of including full attenuation in the calculation of sensitivity kernels by showing significant differences with physical-dispersion-only kernels.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 23
    Publication Date: 2016-07-13
    Description: Real Earth media are not perfectly elastic. Instead, they attenuate propagating mechanical waves. This anelastic phenomenon in wave propagation can be modeled by a viscoelastic mechanical model consisting of several standard linear solids. Using this viscoelastic model, we approximate a constant Q over a frequency band of interest. We use a four-element viscoelastic model with a trade-off between accuracy and computational costs to incorporate Q into 2-D time-domain first-order velocity–stress wave equations. To improve the computational efficiency, we limit the Q in the model to a list of discrete values between 2 and 1000. The related stress and strain relaxation times that characterize the viscoelastic model are pre-calculated and stored in a database for use by the finite-difference calculation. A viscoelastic finite-difference scheme that is second order in time and fourth order in space is developed based on the MacCormack algorithm. The new method is validated by comparing the numerical result with analytical solutions that are calculated using the generalized reflection/transmission coefficient method. The synthetic seismograms exhibit greater than 95 per cent consistency in a two-layer viscoelastic model. The dispersion generated from the simulation is consistent with the Kolsky–Futterman dispersion relationship.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 24
    Publication Date: 2016-07-13
    Description: After 83 yr, the great normal-faulting earthquake of 1933 March 2, which took place off the Japan Trench and produced a devastating tsunami on the Sanriku coast and damaging waves in Hawaii, remains the largest recorded normal-faulting earthquake. This study uses advanced methods to investigate this event using far-field seismological and tsunami data and complements a sister study by Uchida et al. which used exclusively arrival times at Japanese stations. Our relocation of the main shock (39.22°N, 144.45°E, with a poorly constrained depth of less than 40 km) places it in the outer trench slope, below a seafloor depth of ~6500 m, in a region of horst-and-graben structure, with fault scarps approximately parallel to the axis of the Japan Trench. Relocated aftershocks show a band of genuine shallow aftershocks parallel to the Japan Trench under the outer trench slope and a region of post-mainshock events landward of the trench axis that occur over roughly the same latitude range and are thought to be the result of stress transfer to the interplate thrust boundary following the normal-faulting rupture. Based on a combination of P -wave first motions and inversion of surface wave spectral amplitudes, we propose a normal-faulting focal mechanism ( = 200°, = 61° and = 271°) and a seismic moment M 0 = (7 ± 1) x 10 28 dyn cm ( M w = 8.5). A wide variety of data, including the distribution of isoseismals, the large magnitudes (up to 8.9) proposed by early investigators before the standardization of magnitude scales, estimates of energy-to-moment ratios and the tentative identification of a T wave at Pasadena (and possibly Riverside), clearly indicate that this seismic source was exceptionally rich in high-frequency wave energy, suggesting a large apparent stress and a sharp rise time, and consistent with the behaviour of many smaller shallow normal-faulting earthquakes. Hydrodynamic simulations based on a range of possible sources consistent with the above findings, including a compound rupture on two opposite-facing normal-faulting segments, are in satisfactory agreement with tsunami observations in Hawaii, where run-up reached 3 m, causing significant damage. This study emphasizes the need to include off-trench normal-faulting earthquake sources in global assessments of tsunami hazards emanating from the subduction of old and cold plates, whose total length of trenches exceed 20 000 km, even though only a handful of great such events are known with confidence in the instrumental record.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 25
    Publication Date: 2016-07-15
    Description: To refine the 3-D seismic velocity model in the greater Parkfield, California region, a new data set including regular earthquakes, shots, quarry blasts and low-frequency earthquakes (LFEs) was assembled. Hundreds of traces of each LFE family at two temporary arrays were stacked with time–frequency domain phase weighted stacking method to improve signal-to-noise ratio. We extend our model resolution to lower crustal depth with LFE data. Our result images not only previously identified features but also low velocity zones (LVZs) in the area around the LFEs and the lower crust beneath the southern Rinconada Fault. The former LVZ is consistent with high fluid pressure that can account for several aspects of LFE behaviour. The latter LVZ is consistent with a high conductivity zone in magnetotelluric studies. A new Vs model was developed with S picks that were obtained with a new autopicker. At shallow depth, the low Vs areas underlie the strongest shaking areas in the 2004 Parkfield earthquake. We relocate LFE families and analyse the location uncertainties with the NonLinLoc and tomoDD codes. The two methods yield similar results.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 26
    Publication Date: 2016-07-20
    Description: The Cadzow rank-reduction method can be effectively utilized in simultaneously denoising and reconstructing 5-D seismic data that depend on four spatial dimensions. The classic version of Cadzow rank-reduction method arranges the 4-D spatial data into a level-four block Hankel/Toeplitz matrix and then applies truncated singular value decomposition (TSVD) for rank reduction. When the observed data are extremely noisy, which is often the feature of real seismic data, traditional TSVD cannot be adequate for attenuating the noise and reconstructing the signals. The reconstructed data tend to contain a significant amount of residual noise using the traditional TSVD method, which can be explained by the fact that the reconstructed data space is a mixture of both signal subspace and noise subspace. In order to better decompose the block Hankel matrix into signal and noise components, we introduced a damping operator into the traditional TSVD formula, which we call the damped rank-reduction method. The damped rank-reduction method can obtain a perfect reconstruction performance even when the observed data have extremely low signal-to-noise ratio. The feasibility of the improved 5-D seismic data reconstruction method was validated via both 5-D synthetic and field data examples. We presented comprehensive analysis of the data examples and obtained valuable experience and guidelines in better utilizing the proposed method in practice. Since the proposed method is convenient to implement and can achieve immediate improvement, we suggest its wide application in the industry.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 27
    Publication Date: 2016-07-21
    Description: An exact analytical solution is presented for the effective dynamic transverse shear modulus in a heterogeneous fluid-filled porous solid containing cylindrical inclusions. The complex and frequency-dependent properties of the dynamic shear modulus are caused by the physical mechanism of mesoscopic-scale wave-induced fluid flow whose scale is smaller than wavelength but larger than the size of pores. Our model consists of three phases: a long cylindrical inclusion, a cylindrical shell of poroelastic matrix material with different mechanical and/or hydraulic properties than the inclusion and an outer region of effective homogeneous medium of laterally infinite extent. The behavior of both the inclusion and the matrix is described by Biot's consolidation equations, whereas the surrounding effective medium which is used to describe the effective transverse shear properties of the inner poroelastic composite is assumed to be a viscoelastic solid whose complex transverse shear modulus needs to be determined. The determined effective transverse shear modulus is used to quantify the S -wave attenuation and velocity dispersion in heterogeneous fluid-filled poroelastic rocks. The calculation shows the relaxation frequency and relative position of various fluid saturation dispersion curves predicted by this study exhibit very good agreement with those of a previous 2-D finite-element simulation. For the double-porosity model (inclusions having a different solid frame than the matrix but the same pore fluid as the matrix) the effective shear modulus also exhibits a size-dependent characteristic that the relaxation frequency moves to lower frequencies by two orders of magnitude if the radius of the cylindrical poroelastic composite increases by one order of magnitude. For the patchy-saturation model (inclusions having the same solid frame as the matrix but with a different pore fluid from the matrix), the heterogeneity in pore fluid cannot cause any attenuation in the transverse shear modulus at all. A comparison with the case of spherical inclusions illustrates that the transverse shear modulus for the cylindrical inclusion exhibits more S -wave attenuation than spherical inclusions.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 28
    Publication Date: 2016-07-24
    Description: Secondary microseism sources are pressure fluctuations close to the ocean surface. They generate acoustic P waves that propagate in water down to the ocean bottom where they are partly reflected and partly transmitted into the crust to continue their propagation through the Earth. We present the theory for computing the displacement power spectral density of secondary microseism P waves recorded by receivers in the far field. In the frequency domain, the P -wave displacement can be modeled as the product of (1) the pressure source, (2) the source site effect that accounts for the constructive interference of multiply reflected P waves in the ocean, (3) the propagation from the ocean bottom to the stations and (4) the receiver site effect. Secondary microseism P waves have weak amplitudes, but they can be investigated by beamforming analysis. We validate our approach by analysing the seismic signals generated by typhoon Ioke (2006) and recorded by the Southern California Seismic Network. Backprojecting the beam onto the ocean surface enables to follow the source motion. The observed beam centroid is in the vicinity of the pressure source derived from the ocean wave model WAVEWATCH III R . The pressure source is then used for modeling the beam and a good agreement is obtained between measured and modeled beam amplitude variation over time. This modeling approach can be used to invert P -wave noise data and retrieve the source intensity and lateral extent.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 29
    Publication Date: 2016-07-24
    Description: In the context of the verification of the Comprehensive Nuclear-Test Ban Treaty in the marine environment, we present a new discriminant based on the empirical observation that hydroacoustic phases recorded at T -phase stations from explosive sources in the water column feature a systematic inverse dispersion, with lower frequencies traveling slower, which is absent from signals emanating from earthquake sources. This difference is present even in the case of the so-called ‘hotspot earthquakes’ occurring inside volcanic edifices featuring steep slopes leading to efficient seismic–acoustic conversions, which can lead to misidentification of such events as explosions when using more classical duration-amplitude discriminants. We propose an algorithm for the compensation of the effect of dispersion over the hydroacoustic path based on a correction to the spectral phase of the ground velocity recorded by the T -phase station, computed individually from the dispersion observed on each record. We show that the application of a standard amplitude-duration algorithm to the resulting compensated time-series satisfactorily identifies records from hotspot earthquakes as generated by dislocation sources, and present a full algorithm, lending itself to automation, for the discrimination of explosive and earthquake sources of hydroacoustic signals at T -phase stations. The only sources not readily identifiable consist of a handful of complex explosions which occurred in the 1970s, believed to involve the testing of advanced weaponry, and which should be independently identifiable through routine vetting by analysts. While we presently cannot provide a theoretical justification to the observation that only explosive sources generate dispersed T phases, we hint that this probably reflects a simpler, and more coherent distribution of acoustic energy among the various modes constituting the wave train, than in the case of dislocation sources embedded in the solid Earth.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 30
    Publication Date: 2016-07-25
    Description: This work presents an innovative strategy to enhance the resolution of surface wave tomography obtained from ambient noise cross-correlation ( C 1 ) by bridging asynchronous seismic networks through the correlation of coda of correlations ( C 3 ). Rayleigh wave group dispersion curves show consistent results between synchronous and asynchronous stations. Rayleigh wave group traveltimes are inverted to construct velocity–period maps with unprecedented resolution for a region covering Mexico and the southern United States. The resulting period maps are then used to regionalize dispersion curves in order to obtain local 1-D shear velocity models ( V S ) of the crust and uppermost mantle in every cell of a grid of 0.4°. The 1-D structures are obtained by iteratively adding layers until reaching a given misfit, and a global tomography model is considered as an input for depths below 150 km. Finally, a high-resolution 3-D V S model is obtained from these inversions. The major structures observed in the 3-D model are in agreement with the tectonic-geodynamic features and with previous regional and local studies. It also offers new insights to understand the present and past tectonic evolution of the region.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 31
    Publication Date: 2016-07-29
    Description: Determination of a response of the sea water column to teleseismic plane wave is important to suppress adverse effects of water reverberations in calculating receiver functions (RFs) using ocean-bottom seismometer (OBS) records. We present a novel non-linear waveform analysis method using the simulated annealing algorithm to determine such a water-layer response recorded by an OBS array. We then demonstrate its usefulness for the RF estimation through its application to synthetic and observed data. Synthetic experiments suggest that the water-layer response constrained in this way has a potential to improve RFs of OBS records drastically even in the high-frequency range (to 4 Hz). By applying it to data observed by the OBS array around the Kii Peninsula, southwestern Japan, we identified a low-velocity zone at the top of the subducting Philippine Sea plate. This zone may represent the incoming fluid-rich sediment layer that has been reported by active-source seismic survey.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 32
    Publication Date: 2016-08-05
    Description: We develop an automated strategy for discriminating deep microseismic events from shallow ones on the basis of the waveforms recorded on a limited number of surface receivers. Machine-learning techniques are employed to explore the relationship between event hypocentres and seismic features of the recorded signals in time, frequency and time–frequency domains. We applied the technique to 440 microearthquakes –1.7 〈  M w  〈 1.29, induced by an underground cavern collapse in the Napoleonville Salt Dome in Bayou Corne, Louisiana. Forty different seismic attributes of whole seismograms including degree of polarization and spectral attributes were measured. A selected set of features was then used to train the system to discriminate between deep and shallow events based on the knowledge gained from existing patterns. The cross-validation test showed that events with depth shallower than 250 m can be discriminated from events with hypocentral depth between 1000 and 2000 m with 88 per cent and 90.7 per cent accuracy using logistic regression and artificial neural network models, respectively. Similar results were obtained using single station seismograms. The results show that the spectral features have the highest correlation to source depth. Spectral centroids and 2-D cross-correlations in the time–frequency domain are two new seismic features used in this study that showed to be promising measures for seismic event classification. The used machine-learning techniques have application for efficient automatic classification of low energy signals recorded at one or more seismic stations.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 33
    Publication Date: 2016-08-05
    Description: We conduct a numerical experiment to investigate potential bias in measurements of S -wave splitting (apparent differences between the arrival times of SH and SV phases) for waves propagating close to the core–mantle boundary (CMB) in the D'' layer. The bias is defined as the discrepancy between shear wave splitting measured from finite frequency synthetic seismograms (‘apparent splitting’) and the splitting predicted by ray theory, which is a high-frequency approximation. For simple isotropic models, we find biases which are typically between 0.5 and 4 s, depending on the model, the Q structure and the dominant period of the synthetics. The bias increases for lower frequencies or lower Q values. The epicentral distance at which the bias starts depends on the frequency and the Q structure. We also compute synthetics for models based on mineral physics (using the elastic constants under lower-mantle pressure and temperature conditions, taking into account the phase transition from Mg-perovskite to Mg-post-perovskite) and geodynamics (the thermal boundary layer) and find that the depth of the positive velocity jump associated with the phase transition and the depth range over which the velocity decreases (due to temperature increases) in the thermal boundary layer significantly influence the wavefield in the lowermost mantle. For example, in cold regions beneath subduction zones, wavefields for SH and SV differ greatly due to the steep velocity decrease close to the CMB. For complex models, apparent splitting can also arise from the possibility that low amplitude direct phases might be overlooked, and larger amplitude later phases might instead incorrectly be picked as the direct arrival. Biases of the type investigated in this study combine with other sources of uncertainty for splitting in D'' (e.g. the correction for upper-mantle anisotropy and the difference between SH and SV ray paths) to make a precise evaluation of the anisotropy in D'' difficult.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 34
    Publication Date: 2016-07-16
    Description: A velocity ( Vs ) and structure model is derived for the Los Angeles Basin, California based on ambient-noise surface wave and receiver-function analysis, using data from a low-cost, short-duration, dense broad-band survey (LASSIE) deployed across the basin. The shear wave velocities show lateral variations at the Compton-Los Alamitos and the Whittier Faults. The basement beneath the Puente Hills–San Gabriel Valley shows an unusually high velocity (~4.0 km s –1 ) and indicates the presence of schist. The structure of the model shows that the basin is a maximum of 8 km deep along the profile and that the Moho rises to a depth of 17 km under the basin. The basin has a stretch factor of 2.6 in the centre grading to 1.3 at the edges and is in approximate isostatic equilibrium.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 35
    Publication Date: 2016-08-01
    Description: We present a catalogue of full seismic moment tensors for 63 events from Uturuncu volcano in Bolivia. The events were recorded during 2011–2012 in the PLUTONS seismic array of 24 broad-band stations. Most events had magnitudes between 0.5 and 2.0 and did not generate discernible surface waves; the largest event was M w 2.8. For each event we computed the misfit between observed and synthetic waveforms, and we used first-motion polarity measurements to reduce the number of possible solutions. Each moment tensor solution was obtained using a grid search over the 6-D space of moment tensors. For each event, we show the misfit function in eigenvalue space, represented by a lune. We identify three subsets of the catalogue: (1) six isotropic events, (2) five tensional crack events, and (3) a swarm of 14 events southeast of the volcanic centre that appear to be double couples. The occurrence of positively isotropic events is consistent with other published results from volcanic and geothermal regions. Several of these previous results, as well as our results, cannot be interpreted within the context of either an oblique opening crack or a crack-plus-double-couple model. Proper characterization of uncertainties for full moment tensors is critical for distinguishing among physical models of source processes.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 36
    Publication Date: 2016-08-02
    Description: The 2-D acoustic wave equation is commonly solved numerically by finite-difference (FD) methods in which the accuracy of solution is significantly affected by the FD stencils. The commonly used cross stencil can reach either only second-order accuracy for space domain dispersion-relation-based FD method or (2 M )th-order accuracy along eight specific propagation directions for time–space domain dispersion-relation-based FD method, if the conventional (2 M )th-order spatial FD and second-order temporal FD are used to discretize the equation. One other newly developed rhombus stencil can reach arbitrary even-order accuracy. However, this stencil adds significantly to computational cost when the operator length is large. To achieve a balance between the solution accuracy and efficiency, we develop a new FD stencil to solve the 2-D acoustic wave equation. This stencil is a combination of the cross stencil and rhombus stencil. A cross stencil with an operator length parameter M is used to approximate the spatial partial derivatives while a rhombus stencil with an operator length parameter N together with the conventional second-order temporal FD is employed in approximating the temporal partial derivatives. Using this stencil, a new FD scheme is developed; we demonstrate that this scheme can reach (2 M )th-order accuracy in space and (2 N )th-order accuracy in time when spatial FD coefficients and temporal FD coefficients are derived from respective dispersion relation using Taylor-series expansion (TE) method. To further increase the accuracy, we derive the FD coefficients by employing the time–space domain dispersion relation of this FD scheme using TE. We also use least-squares (LS) optimization method to reduce dispersion at high wavenumbers. Dispersion analysis, stability analysis and modelling examples demonstrate that our new scheme has greater accuracy and better stability than conventional FD schemes, and thus can adopt large time steps. To reduce the extra computational cost resulting from adopting the new stencil, we apply the variable spatial operator length schemes. Adopting our new FD scheme, characterized by new stencil, LS-based optimization, variable operator lengths and larger time step, modelling efficiency is significantly improved.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 37
    Publication Date: 2016-07-08
    Description: Near-surface geophysical imaging is often performed by generating surface waves, and estimating the subsurface properties through inversion, that is, iteratively matching experimentally observed dispersion curves with predicted curves from a layered half-space model of the subsurface. Key to the effectiveness of inversion is the efficiency and accuracy of computing the dispersion curves and their derivatives. This paper presents improved methodologies for both dispersion curve and derivative computation. First, it is shown that the dispersion curves can be computed more efficiently by combining an unconventional complex-length finite element method (CFEM) to model the finite depth layers, with perfectly matched discrete layers (PMDL) to model the unbounded half-space. Second, based on analytical derivatives for theoretical dispersion curves, an approximate derivative is derived for the so-called effective dispersion curve for realistic geophysical surface response data. The new derivative computation has a smoothing effect on the computation of derivatives, in comparison with traditional finite difference (FD) approach, and results in faster convergence. In addition, while the computational cost of FD differentiation is proportional to the number of model parameters, the new differentiation formula has a computational cost that is almost independent of the number of model parameters. At the end, as confirmed by synthetic and real-life imaging examples, the combination of CFEM + PMDL for dispersion calculation and the new differentiation formula results in more accurate estimates of the subsurface characteristics than the traditional methods, at a small fraction of computational effort.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 38
    Publication Date: 2016-07-08
    Description: Powerful subduction zone earthquakes rupture thousands of square kilometres along continental margins but at certain locations earthquake rupture terminates. To date, detailed knowledge of the parameters that govern seismic rupture and aftershocks is still incomplete. On 2015 September 16, the M w 8.3 Illapel earthquake ruptured a 200 km long stretch of the Central Chilean subduction zone, triggering a tsunami and causing significant damage. Here, we analyse the temporal and spatial pattern of the coseismic rupture and aftershocks in relation to the tectonic setting in the earthquake area. Aftershocks cluster around the area of maximum coseismic slip, in particular in lateral and downdip direction. During the first 24 hr after the main shock, aftershocks migrated in both lateral directions with velocities of approximately 2.5 and 5 km hr –1 . At the southern rupture boundary, aftershocks cluster around individual subducted seamounts that are related to the downthrusting Juan Fernández Ridge. In the northern part of the rupture area, aftershocks separate into an upper cluster (above 25 km depth) and a lower cluster (below 35 km depth). This dual seismic–aseismic transition in downdip direction is also observed in the interseismic period suggesting that it may represent a persistent feature for the Central Chilean subduction zone.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 39
    Publication Date: 2016-06-22
    Description: Least-squares inversion of seismic arrivals can provide remarkably detailed models of the Earth's subsurface. However, cycle skipping associated with these oscillatory arrivals is the main cause for local minima in the least-squares objective function. Therefore, it is often difficult for descent methods to converge to the solution without an accurate initial large-scale velocity estimate. The low frequencies in the arrivals, needed to update the large-scale components in the velocity model, are usually unreliable or absent. To overcome this difficulty, we propose a multi-objective inversion scheme that uses the conventional least-squares functional along with an auxiliary data-domain objective. As the auxiliary objective effectively replaces the seismic arrivals by bumps, we call it the bump functional. The bump functional minimization can be made far less sensitive to cycle skipping and can deal with multiple arrivals in the data. However, it can only be used as an auxiliary objective since it usually does not provide a unique model after minimization even when the regularized-least-squares functional has a unique global minimum and hence a unique solution. The role of the bump functional during the multi-objective inversion is to guide the optimization towards the global minimum by pulling the trapped solution out of the local minima associated with the least-squares functional whenever necessary. The computational complexity of the bump functional is equivalent to that of the least-squares functional. In this paper, we describe various characteristics of the bump functional using simple and illustrative numerical examples. We also demonstrate the effectiveness of the proposed multi-objective inversion scheme by considering more realistic examples. These include synthetic and field data from a cross-well experiment, surface-seismic synthetic data with reflections and synthetic data with refracted arrivals at long offsets.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 40
    Publication Date: 2016-05-05
    Description: Owing to the increasing availability of computational resources, in recent years the probabilistic solution of non-linear, geophysical inverse problems by means of sampling methods has become increasingly feasible. Nevertheless, we still face situations in which a Monte Carlo approach is not practical. This is particularly true in cases where the evaluation of the forward problem is computationally intensive or where inversions have to be carried out repeatedly or in a timely manner, as in natural hazards monitoring tasks such as earthquake early warning. Here, we present an alternative to Monte Carlo sampling, in which inferences are entirely based on a set of prior samples—that is, samples that have been obtained independent of a particular observed datum. This has the advantage that the computationally expensive sampling stage becomes separated from the inversion stage, and the set of prior samples—once obtained—can be reused for repeated evaluations of the inverse mapping without additional computational effort. This property is useful if the problem is such that repeated inversions of independent data have to be carried out. We formulate the inverse problem in a Bayesian framework and present a practical way to make posterior inferences based on a set of prior samples. We compare the prior sampling based approach to a Markov Chain Monte Carlo approach that samples from the posterior probability distribution. We show results for both a toy example, and a realistic seismological source parameter estimation problem. We find that the posterior uncertainty estimates obtained based on prior sampling can be considered conservative estimates of the uncertainties obtained by directly sampling from the posterior distribution.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 41
    Publication Date: 2016-05-05
    Description: The Afar Depression and its adjacent areas are underlain by an upper mantle marked by some of the world's largest negative velocity anomalies, which are frequently attributed to the thermal influences of a lower-mantle plume. In spite of numerous studies, however, the existence of a plume beneath the area remains enigmatic, partially due to inadequate quantities of broad-band seismic data and the limited vertical resolution at the mantle transition zone (MTZ) depth of the techniques employed by previous investigations. In this study, we use an unprecedented quantity (over 14 500) of P -to- S receiver functions (RFs) recorded by 139 stations from 12 networks to image the 410 and 660 km discontinuities and map the spatial variation of the thickness of the MTZ. Non-linear stacking of the RFs under a 1-D velocity model shows robust P -to- S conversions from both discontinuities, and their apparent depths indicate the presence of an upper-mantle low-velocity zone beneath the entire study area. The Afar Depression and the northern Main Ethiopian Rift are characterized by an apparent 40–60 km depression of both MTZ discontinuities and a normal MTZ thickness. The simplest and most probable interpretation of these observations is that the apparent depressions are solely caused by velocity perturbations in the upper mantle and not by deeper processes causing temperature or hydration anomalies within the MTZ. Thickening of the MTZ on the order of 15 km beneath the southern Arabian Plate, southern Red Sea and western Gulf of Aden, which comprise the southward extension of the Afro-Arabian Dome, could reflect long-term hydration of the MTZ. A 20 km thinning of the MTZ beneath the western Ethiopian Plateau is observed and interpreted as evidence for a possible mantle plume stem originating from the lower mantle.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 42
    Publication Date: 2016-05-05
    Description: We propose a procedure for uncertainty quantification in Probabilistic Tsunami Hazard Analysis (PTHA), with a special emphasis on the uncertainty related to statistical modelling of the earthquake source in Seismic PTHA (SPTHA), and on the separate treatment of subduction and crustal earthquakes (treated as background seismicity). An event tree approach and ensemble modelling are used in spite of more classical approaches, such as the hazard integral and the logic tree. This procedure consists of four steps: (1) exploration of aleatory uncertainty through an event tree, with alternative implementations for exploring epistemic uncertainty; (2) numerical computation of tsunami generation and propagation up to a given offshore isobath; (3) (optional) site-specific quantification of inundation; (4) simultaneous quantification of aleatory and epistemic uncertainty through ensemble modelling. The proposed procedure is general and independent of the kind of tsunami source considered; however, we implement step 1, the event tree, specifically for SPTHA, focusing on seismic source uncertainty. To exemplify the procedure, we develop a case study considering seismic sources in the Ionian Sea (central-eastern Mediterranean Sea), using the coasts of Southern Italy as a target zone. The results show that an efficient and complete quantification of all the uncertainties is feasible even when treating a large number of potential sources and a large set of alternative model formulations. We also find that (i) treating separately subduction and background (crustal) earthquakes allows for optimal use of available information and for avoiding significant biases; (ii) both subduction interface and crustal faults contribute to the SPTHA, with different proportions that depend on source-target position and tsunami intensity; (iii) the proposed framework allows sensitivity and deaggregation analyses, demonstrating the applicability of the method for operational assessments.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 43
    Publication Date: 2016-05-05
    Description: The expanding fleet of broad-band ocean-bottom seismograph (OBS) stations is facilitating the study of the structure and seismicity of oceanic plates at regional scales. For continental studies, an important tool to characterize continental crust and mantle structure is the analysis of teleseismic P receiver functions. In the oceans, however, receiver functions potentially suffer from several limiting factors that are unique to ocean sites and plate structures. In this study, we model receiver functions for a variety of oceanic lithospheric structures to investigate the possibilities and limitations of receiver functions using OBS data. Several potentially contaminating effects are examined, including pressure reverberations from the water column for various ocean-floor depths and the effects of a layer of low-velocity marine sediments. These modelling results indicate that receiver functions from OBS data are difficult to interpret in the presence of marine sediments, but shallow-water sites in subduction zone forearcs may be suitable for constraining various crustal elements around the locked megathrust fault. We propose using a complementary approach based on transfer function modelling combined with a grid search approach that bypasses receiver functions altogether and estimates model properties directly from minimally processed waveforms. Using real data examples from the Cascadia Initiative, we show how receiver and transfer functions can be used to infer seismic properties of the oceanic plate in both shallow (Cascadia forearc) and deep (Juan de Fuca Ridge) ocean settings.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 44
    Publication Date: 2016-05-07
    Description: A new implementation of indirect boundary element method allows simulating the elastic wave propagation in complex configurations made of embedded regions that are homogeneous with irregular boundaries or flat layers. In an older implementation, each layer of a flat layered region would have been treated as a separated homogeneous region without taking into account the flat boundary information. For both types of regions, the scattered field results from fictitious sources positioned along their boundaries. For the homogeneous regions, the fictitious sources emit as in a full-space and the wave field is given by analytical Green's functions. For flat layered regions, fictitious sources emit as in an unbounded flat layered region and the wave field is given by Green's functions obtained from the discrete wavenumber (DWN) method. The new implementation allows then reducing the length of the discretized boundaries but DWN Green's functions require much more computation time than the full-space Green's functions. Several optimization steps are then implemented and commented. Validations are presented for 2-D and 3-D problems. Higher efficiency is achieved in 3-D.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 45
    Publication Date: 2016-05-11
    Description: Though well known for layer boundaries, the use of amplitude-versus-offset (AVO) variations for non-welded boundaries like fractures is not yet investigated. Depending on the seismic wavelength used, fractures can be regarded as thin, compliant zones in rocks, in different scales. We explore the potential of multiangle AVO inversion of P-P and P-S reflections from a fracture to estimate fracture properties. We conduct laboratory experiments to measure reflection responses of dry and wet fractures. The observed P-P reflections of the wet fracture and the fracture aperture are very well predicted by the non-welded interface model. We invert the angle-dependent P-P reflectivity of the fracture to estimate both normal and tangential fracture compliances. The estimated value of the normal compliance is accurate, and it is also possible to obtain the value of the non-zero tangential compliance. We find that supplementing the information of converted P-S reflections in the AVO inversion greatly improves the estimate of the tangential compliance. The calculated compliance ratio clearly shows the existence of fluid in the fracture. This finding can be crucial for new applications in a wide range of scale—from earthquake seismology, deep and shallow seismic exploration, to non-destructive material testing.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 46
    Publication Date: 2016-05-11
    Description: The present study evaluates the capacity of the Boom Clay as a host rock for disposal purposes, more precisely its seismic characterization, which may assess its long-term performance to store radioactive wastes. Although the formation is relatively uniform and homogeneous, there are embedded thin layers of septaria (carbonates) that may affect the integrity of the Boom Clay. Therefore, it is essential to locate these geobodies. The seismic data to characterize the Boom Clay has been acquired at the Kruibeke test site. The inversion, which allowed us to obtain the anisotropy parameters and seismic velocities of the clay, is complemented with further information such as log and laboratory data. The attenuation properties have been estimated from equivalent formations (having similar composition and seismic velocities). The inversion yields quite consistent results although the symmetry of the medium is unusual but physically possible, since the anisotropy parameter is negative. According to a time-domain calculation of the energy velocity at four frequency bands up to 900 Hz, velocity increases with frequency, a behaviour described by the Zener model. Then, we use this model to describe anisotropy and anelasticity that are implemented into the equation of motion to compute synthetic seismograms in the space–time domain. The technique is based on memory variables and the Fourier pseudospectral method. We have computed reflection coefficients of the septaria thin layer. At normal incidence, the P -wave coefficient vanishes at specific thicknesses of the layer and there is no conversion to the S wave. For example, calculations at 600 Hz show that for thicknesses of 1 m the septarium can be detected more easily since the amplitudes are higher (nearly 0.8). Converted PS waves have a high amplitude at large offsets (between 30° and 80°) and can be useful to identify the target on this basis. Moreover, we have investigated the effect of septaria embedded in the Boom Clay with several simulations, by considering a lateral partial continuity of the calcareous thin inclusions. The simulations with layers of calcareous material show continuity of the reflections even when the percentage of carbonate within the layer is very small (5–15 per cent), while for low content of the calcareous material, isolated septaria boulders generate diffraction events. We have also simulated the stacked seismic section obtained from processing of the field data. The matching between the field and synthetic sections is acceptable.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 47
    Publication Date: 2016-05-27
    Description: A new Matched Filtering Algorithm (MFA) is proposed for detecting and analysing microseismic events recorded by downhole monitoring of hydraulic fracturing. This method requires a set of well-located template (‘parent’) events, which are obtained using conventional microseismic processing and selected on the basis of high signal-to-noise (S/N) ratio and representative spatial distribution of the recorded microseismicity. Detection and extraction of ‘child’ events are based on stacked, multichannel cross-correlation of the continuous waveform data, using the parent events as reference signals. The location of a child event relative to its parent is determined using an automated process, by rotation of the multicomponent waveforms into the ray-centred co-ordinates of the parent and maximizing the energy of the stacked amplitude envelope within a search volume around the parent's hypocentre. After correction for geometrical spreading and attenuation, the relative magnitude of the child event is obtained automatically using the ratio of stacked envelope peak with respect to its parent. Since only a small number of parent events require interactive analysis such as picking P - and S -wave arrivals, the MFA approach offers the potential for significant reduction in effort for downhole microseismic processing. Our algorithm also facilitates the analysis of single-phase child events, that is, microseismic events for which only one of the S - or P -wave arrivals is evident due to unfavourable S/N conditions. A real-data example using microseismic monitoring data from four stages of an open-hole slickwater hydraulic fracture treatment in western Canada demonstrates that a sparse set of parents (in this case, 4.6 per cent of the originally located events) yields a significant (more than fourfold increase) in the number of located events compared with the original catalogue. Moreover, analysis of the new MFA catalogue suggests that this approach leads to more robust interpretation of the induced microseismicity and novel insights into dynamic rupture processes based on the average temporal (foreshock–aftershock) relationship of child events to parents.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 48
    Publication Date: 2016-05-27
    Description: Azimuthal anisotropy is a powerful tool to reveal information about both the present structure and past evolution of the mantle. Anisotropic images of the upper mantle are usually obtained by analysing various types of seismic observables, such as surface wave dispersion curves or waveforms, SKS splitting data, or receiver functions. These different data types sample different volumes of the earth, they are sensitive to different length scales, and hence are associated with different levels of uncertainties. They are traditionally interpreted separately, and often result in incompatible models. We present a Bayesian inversion approach to jointly invert these different data types. Seismograms for SKS and P phases are directly inverted using a cross-convolution approach, thus avoiding intermediate processing steps, such as numerical deconvolution or computation of splitting parameters. Probabilistic 1-D profiles are obtained with a transdimensional Markov chain Monte Carlo scheme, in which the number of layers, as well as the presence or absence of anisotropy in each layer, are treated as unknown parameters. In this way, seismic anisotropy is only introduced if required by the data. The algorithm is used to resolve both isotropic and anisotropic layering down to a depth of 350 km beneath two seismic stations in North America in two different tectonic settings: the stable Canadian shield (station FFC) and the tectonically active southern Basin and Range Province (station TA-214A). In both cases, the lithosphere–asthenosphere boundary is clearly visible, and marked by a change in direction of the fast axis of anisotropy. Our study confirms that azimuthal anisotropy is a powerful tool for detecting layering in the upper mantle.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 49
    Publication Date: 2016-05-27
    Description: We have obtained new results in the statistical analysis of global earthquake catalogues with special attention to the largest earthquakes, and we examined the statistical behaviour of earthquake rate variations. These results can serve as an input for updating our recent earthquake forecast, known as the ‘Global Earthquake Activity Rate 1’ model (GEAR1), which is based on past earthquakes and geodetic strain rates. The GEAR1 forecast is expressed as the rate density of all earthquakes above magnitude 5.8 within 70 km of sea level everywhere on earth at 0.1 x 0.1 degree resolution, and it is currently being tested by the Collaboratory for Study of Earthquake Predictability. The seismic component of the present model is based on a smoothed version of the Global Centroid Moment Tensor (GCMT) catalogue from 1977 through 2013. The tectonic component is based on the Global Strain Rate Map, a ‘General Earthquake Model’ (GEM) product. The forecast was optimized to fit the GCMT data from 2005 through 2012, but it also fit well the earthquake locations from 1918 to 1976 reported in the International Seismological Centre-Global Earthquake Model (ISC-GEM) global catalogue of instrumental and pre-instrumental magnitude determinations. We have improved the recent forecast by optimizing the treatment of larger magnitudes and including a longer duration (1918–2011) ISC-GEM catalogue of large earthquakes to estimate smoothed seismicity. We revised our estimates of upper magnitude limits, described as corner magnitudes, based on the massive earthquakes since 2004 and the seismic moment conservation principle. The new corner magnitude estimates are somewhat larger than but consistent with our previous estimates. For major subduction zones we find the best estimates of corner magnitude to be in the range 8.9 to 9.6 and consistent with a uniform average of 9.35. Statistical estimates tend to grow with time as larger earthquakes occur. However, by using the moment conservation principle that equates the seismic moment rate with the tectonic moment rate inferred from geodesy and geology, we obtain a consistent estimate of the corner moment largely independent of seismic history. These evaluations confirm the above-mentioned corner magnitude value. The new estimates of corner magnitudes are important both for the forecast part based on seismicity as well as the part based on geodetic strain rates. We examine rate variations as expressed by annual earthquake numbers. Earthquakes larger than magnitude 6.5 obey the Poisson distribution. For smaller events the negative-binomial distribution fits much better because it allows for earthquake clustering.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 50
    Publication Date: 2016-05-28
    Description: In order to improve our understanding of hazardous underground cavities, the development and collapse of a ~200 m wide salt solution mining cavity was seismically monitored in the Lorraine basin in northeastern France. The microseismic events show a swarm-like behaviour, with clustering sequences lasting from seconds to days, and distinct spatiotemporal migration. Observed microseismic signals are interpreted as the result of detachment and block breakage processes occurring at the cavity roof. Body wave amplitude patterns indicated the presence of relatively stable source mechanisms, either associated with dip-slip and/or tensile faulting. Signal overlaps during swarm activity due to short interevent times, the high-frequency geophone recordings and the limited network station coverage often limit the application of classical source analysis techniques. To overcome these shortcomings, we investigated the source mechanisms through different procedures including modelling of observed and synthetic waveforms and amplitude spectra of some well-located events, as well as modelling of peak-to-peak amplitude ratios for the majority of the detected events. We extended the latter approach to infer the average source mechanism of many swarming events at once, using multiple events recorded at a single three component station. This methodology is applied here for the first time and represents a useful tool for source studies of seismic swarms and seismicity clusters. The results obtained with different methods are consistent and indicate that the source mechanisms for at least 50 per cent of the microseismic events are remarkably stable, with a predominant thrust faulting regime with faults similarly oriented, striking NW–SE and dipping around 35°–55°. This dominance of consistent source mechanisms might be related to the presence of a preferential direction of pre-existing crack or fault structures. As an interesting byproduct, we demonstrate, for the first time directly on seismic data, that the source radiation pattern significantly controls the detection capability of a seismic station and network.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 51
    Publication Date: 2016-07-16
    Description: We show that higher modes are an important component of high-frequency Rayleigh waves in the cross-correlations over sedimentary basins. The particle motions provide a good test for distinguishing and separating the fundamental from the first higher mode, with the fundamental mode having retrograde and the first higher mode having prograde motion in the 1–10 s period of interest. The basement depth controls the cut-off period of the first higher mode, which coincides with a rapid increase (over period) in the particle-motion ellipticity or H / V ratio of the fundamental mode. The strong higher mode we observed is not only due to the low-velocity sedimentary layer but also due to the noise sources with significant radial component such as the basin edge scattering. It is important to correctly identify the mode order when inverting the dispersion curves because misidentifying the higher mode as fundamental will lead to an anomalous high V SV velocity.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 52
    Publication Date: 2016-07-16
    Description: We introduce a ‘double-difference’ method for the inversion for seismic wave speed structure based on adjoint tomography. Differences between seismic observations and model predictions at individual stations may arise from factors other than structural heterogeneity, such as errors in the assumed source-time function, inaccurate timings and systematic uncertainties. To alleviate the corresponding non-uniqueness in the inverse problem, we construct differential measurements between stations, thereby reducing the influence of the source signature and systematic errors. We minimize the discrepancy between observations and simulations in terms of the differential measurements made on station pairs. We show how to implement the double-difference concept in adjoint tomography, both theoretically and practically. We compare the sensitivities of absolute and differential measurements. The former provide absolute information on structure along the ray paths between stations and sources, whereas the latter explain relative (and thus higher resolution) structural variations in areas close to the stations. Whereas in conventional tomography a measurement made on a single earthquake-station pair provides very limited structural information, in double-difference tomography one earthquake can actually resolve significant details of the structure. The double-difference methodology can be incorporated into the usual adjoint tomography workflow by simply pairing up all conventional measurements; the computational cost of the necessary adjoint simulations is largely unaffected. Rather than adding to the computational burden, the inversion of double-difference measurements merely modifies the construction of the adjoint sources for data assimilation.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 53
    Publication Date: 2016-07-20
    Description: The creation of a magnitude-homogenized catalogue is often one of the most fundamental steps in seismic hazard analysis. The process of homogenizing multiple catalogues of earthquakes into a single unified catalogue typically requires careful appraisal of available bulletins, identification of common events within multiple bulletins and the development and application of empirical models to convert from each catalogue's native scale into the required target. The database of the International Seismological Center (ISC) provides the most exhaustive compilation of records from local bulletins, in addition to its reviewed global bulletin. New open-source tools are developed that can utilize this, or any other compiled database, to explore the relations between earthquake solutions provided by different recording networks, and to build and apply empirical models in order to harmonize magnitude scales for the purpose of creating magnitude-homogeneous earthquake catalogues. These tools are described and their application illustrated in two different contexts. The first is a simple application in the Sub-Saharan Africa region where the spatial coverage and magnitude scales for different local recording networks are compared, and their relation to global magnitude scales explored. In the second application the tools are used on a global scale for the purpose of creating an extended magnitude-homogeneous global earthquake catalogue. Several existing high-quality earthquake databases, such as the ISC-GEM and the ISC Reviewed Bulletins, are harmonized into moment magnitude to form a catalogue of more than 562 840 events. This extended catalogue, while not an appropriate substitute for a locally calibrated analysis, can help in studying global patterns in seismicity and hazard, and is therefore released with the accompanying software.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 54
    Publication Date: 2016-07-20
    Description: This paper introduces a novel approach to constructing an effective pre-conditioner for finite-difference (FD) electromagnetic modelling in geophysical applications. This approach is based on introducing an FD contraction operator, similar to one developed for integral equation formulation of Maxwell's equation. The properties of the FD contraction operator were established using an FD analogue of the energy equality for the anomalous electromagnetic field. A new pre-conditioner uses a discrete Green's function of a 1-D layered background conductivity. We also developed the formulae for an estimation of the condition number of the system of FD equations pre-conditioned with the introduced FD contraction operator. Based on this estimation, we have established that the condition number is bounded by the maximum conductivity contrast between the background conductivity and actual conductivity. When there are both resistive and conductive anomalies relative to the background, the new pre-conditioner is advantageous over using the 1-D discrete Green's function directly. In our numerical experiments with both resistive and conductive anomalies, for a land geoelectrical model with 1:10 contrast, the method accelerates convergence of an iterative method (BiCGStab) by factors of 2–2.5, and in a marine example with 1:50 contrast, by a factor of 4.6, compared to direct use of the discrete 1-D Green's function as a pre-conditioner.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 55
    Publication Date: 2016-07-29
    Description: Reaching the global minimum of a waveform misfit function requires careful choices about the nonlinear optimization, preconditioning and regularization methods underlying an inversion. Because waveform inversion problems are susceptible to erratic convergence associated with strong nonlinearity, one or two test cases are not enough to reliably inform such decisions. We identify best practices, instead, using four seismic near-surface problems, one regional problem and two global problems. To make meaningful quantitative comparisons between methods, we carry out hundreds of inversions, varying one aspect of the implementation at a time. Comparing nonlinear optimization algorithms, we find that limited-memory BFGS provides computational savings over nonlinear conjugate gradient methods in a wide range of test cases. Comparing preconditioners, we show that a new diagonal scaling derived from the adjoint of the forward operator provides better performance than two conventional preconditioning schemes. Comparing regularization strategies, we find that projection, convolution, Tikhonov regularization and total variation regularization are effective in different contexts. Besides questions of one strategy or another, reliability and efficiency in waveform inversion depend on close numerical attention and care. Implementation details involving the line search and restart conditions have a strong effect on computational cost, regardless of the chosen nonlinear optimization algorithm.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 56
    Publication Date: 2016-07-29
    Description: The conventional finite-difference time-domain (FDTD) method for elastic waves suffers from the staircasing error when applied to model a curved free surface because of its structured grid. In this work, an improved, stable and accurate 3-D FDTD method for elastic wave modelling on a curved free surface is developed based on the finite volume method and enlarged cell technique (ECT). To achieve a sufficiently accurate implementation, a finite volume scheme is applied to the curved free surface to remove the staircasing error; in the mean time, to achieve the same stability as the FDTD method without reducing the time step increment, the ECT is introduced to preserve the solution stability by enlarging small irregular cells into adjacent cells under the condition of conservation of force. This method is verified by several 3-D numerical examples. Results show that the method is stable at the Courant stability limit for a regular FDTD grid, and has much higher accuracy than the conventional FDTD method.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 57
    Publication Date: 2015-05-01
    Description: The Anatolian Plate is composed of different lithospheric blocks and ribbon continents amalgamated during the closure of the Paleotethys Ocean and Neotethys Ocean along a subduction margin. Using ambient noise tomography, we investigate the crustal and uppermost mantle shear wave velocity structure of the Anatolian Plate. A total of 215 broad-band seismic stations were used spanning 7 yr of recording to compute 13 778 cross-correlations and obtain Rayleigh wave dispersion measurements for periods between 8 and 40 s. We then perform a shear wave inversion to calculate the seismic velocity structure of the crust and uppermost mantle. Our results show that the overall crustal shear wave velocities of the Anatolian crust are low (~3.4 km s –1 ), indicative of a felsic overall composition. We find that prominent lateral seismic velocity gradients correlate with Tethyan suture zones, supporting the idea that the neotectonic structures of Turkey are exploiting the lithospheric weaknesses associated with the amalgamation of Anatolia. Anomalously slow shear wave velocities (~3.15 km s –1 at 25 km) are located in the western limb of the Isparta Angle in southwestern Turkey. In the upper crust, we find that these low shear wave velocities correlate well with the projected location of a carbonate platform unit (Bey Dağlari) beneath the Lycian Nappe complex. In the lower crust and upper mantle of this region, we propose that the anomalously slow velocities are due to the introduction of aqueous fluids related to the underplating of accretionary material from the underthrusting of a buoyant, attenuated continental fragment similar to the Eratosthenes seamount. We suggest that this fragment controlled the location of the formation of the Subduction-Transform Edge Propagator fault in the eastern Aegean Sea during rapid slab rollback of the Aegean Arc in early Miocene times. Lastly, we observe that the uppermost mantle beneath continental Anatolia is generally slow (~4.2 km s –1 ), indicating higher than usual temperatures consistent with the influx of asthenosphere to shallow depths as a result of the segmentation and break-up of the subducting African lithosphere.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 58
    Publication Date: 2015-05-01
    Description: The H/V spectral ratio has emerged as a single station method within the seismic ambient noise analysis field by its capability to quickly estimate the frequency of resonance at a site and through inversion the average profile information. Although it is easy to compute from experimental data, its counter theoretical part is not obvious when building a forward model which can help in reconstructing the derived H/V spectrum. This has led to the simplified assumption that the noise wavefield is mainly composed of Rayleigh waves and the derived H/V often used without further correction. Furthermore, only the right (and left) flank around the H/V peak frequency is considered in the inversion for the subsurface 1-D shear wave velocity profile. A new theoretical approach for the interpretation of the H/V spectral ratio has been presented by Sánchez-Sesma et al. In this paper, the fundamental idea behind their theory is presented as it applies to receivers at depth. A smooth H/V( z , f ) spectral curve on a broad frequency range is obtained by considering a fine integration step which is in turn time consuming. We show that for practical purposes and in the context of inversion, this can be considerably optimized by using a coarse integration step combined with the smoothing of the corresponding directional energy density (DED) spectrum. Further analysis shows that the obtained H/V( z , f ) spectrum computed by the mean of the imaginary part of Green's function method could also be recovered using the reflectivity method for a medium well illuminated by seismic sources. Inversion of synthetic H/V( z , f ) spectral curve is performed for a single layer over a half space. The striking results allow to potentially use the new theory as a forward computation of the H/V( z , f ) to fully invert the experimental H/V spectral ratio at the corresponding depth for the shear velocity profile ( Vs ) and additionally the compressional velocity profile ( Vp ) using receivers both at the surface and in depth. We use seismic ambient noise data in the frequency range of 0.2–50 Hz recorded at two selected sites in Germany where borehole information is also available. The obtained 1-D Vs and Vp profiles are correlated with geological log information. Results from shallow geophysical experiment are also used for comparison.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 59
    Publication Date: 2015-05-01
    Description: The 1000-km-long left-lateral Dead Sea fault is a major tectonic structure of the oriental Mediterranean basin, bounding the Arabian Plate to the west. The fault is located in a region with an exceptionally long and rich historical record, allowing to document historical seismicity catalogues with unprecedented level of details. However, if the earthquake time series is well documented, location and lateral extent of past earthquakes remain often difficult to establish, if only based on historical testimonies. We excavated a palaeoseismic trench in a site located in a kilometre-size extensional jog, south of the Dead Sea, in the Wadi Araba. Based on the stratigraphy exposed in the trench, we present evidence for nine earthquakes that produced surface ruptures during a time period spanning 5000 yr. Abundance of datable material allows us to tie the five most recent events to historical earthquakes with little ambiguities, and to constrain the possible location of these historical earthquakes. The events identified at our site are the 1458 C.E., 1212 C.E., 1068 C.E., one event during the 8th century crisis, and the 363 C.E. earthquake. Four other events are also identified, which correlation with historical events remains more speculative. The magnitude of earthquakes is difficult to assess based on evidence at one site only. The deformation observed in the excavation, however, allows discriminating between two classes of events that produced vertical deformation with one order of amplitude difference, suggesting that we could distinguish earthquakes that started/stopped at our site from earthquakes that potentially ruptured most of the Wadi Araba fault. The time distribution of earthquakes during the past 5000 yr is uneven. The early period shows little activity with return interval of ~500 yr or longer. It is followed by a ~1500-yr-long period with more frequent events, about every 200 yr. Then, for the past ~550 yr, the fault has switched back to a quieter mode with no significant earthquake along the entire southern part of the Dead Sea fault, between the Dead Sea and the Gulf of Aqaba. We computed the Coefficient of Variation for our site and three other sites along the Dead Sea fault, south of Lebanon, to compare time distribution of earthquakes at different locations along the fault. With one exception at a site located next to Lake Tiberias, the three other sites are consistent to show some temporal clustering at the scale of few thousands years.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 60
    Publication Date: 2015-05-10
    Description: The proposed study deals with synthesis of seismograms by BIEM (boundary integral equation method) taking into account all three base components-seismic source, wave path and local region of interest. Consider a laterally inhomogeneous geological profile situated in a half-plane with non-parallel layers. Seismic load is time-harmonic or transient in time. It is presented by incident SH wave or wave radiating from an embedded line seismic source. Two types of lateral inhomogeneities with arbitrary shape and located in the inhomogeneous half-plane are considered: (i) free-surface relief as a canyon or a hill; (ii) alluvial basin with properties different from those of the layered half-plane. The computational tool is BIEM based on the frequency-dependent elastodynamic fundamental solutions. A relation between displacements and tractions along the free surface and arbitrary interface of the soil stratum is derived, which is applicable for arbitrary geometry of the interfaces between soil layers. Validation and convergence study is presented. All simulations reveal the sensitivity of the synthetic seismic signals on the type and characteristics of the seismic time-harmonic or transient load, on the wave path inhomogeneity and on the specific geotechnical properties of the local geological region.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 61
    Publication Date: 2015-05-01
    Description: Time-shift estimation between arrivals in two seismic traces before and after a velocity perturbation is a crucial step in many seismic methods. The accuracy of the estimated velocity perturbation location and amplitude depend on this time shift. Windowed cross-correlation and trace stretching are two techniques commonly used to estimate local time shifts in seismic signals. In the work presented here we implement Dynamic Time Warping (DTW) to estimate the warping function – a vector of local time shifts that globally minimizes the misfit between two seismic traces. We compare all three methods using acoustic numerical experiments. We show that DTW is comparable to or better than the other two methods when the velocity perturbation is homogeneous and the signal-to-noise ratio is high. When the signal-to-noise ratio is low, we find that DTW and windowed cross-correlation are more accurate than the stretching method. Finally, we show that the DTW algorithm has good time resolution when identifying small differences in the seismic traces for a model with an isolated velocity perturbation. These results impact current methods that utilize not only time shifts between (multiply) scattered waves, but also amplitude and decoherence measurements.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 62
    Publication Date: 2015-05-03
    Description: We discuss several outstanding aspects of seismograms recorded during 〉4 weeks by a spatially dense Nodal array, straddling the damage zone of the San Jacinto fault in southern California, and some example results. The waveforms contain numerous spikes and bursts of high-frequency waves (up to the recorded 200 Hz) produced in part by minute failure events in the shallow crust. The high spatial density of the array facilitates the detection of 120 small local earthquakes in a single day, most of which not detected by the surrounding ANZA and regional southern California networks. Beamforming results identify likely ongoing cultural noise sources dominant in the frequency range 1–10 Hz and likely ongoing earthquake sources dominant in the frequency range 20–40 Hz. Matched-field processing and back-projection of seismograms provide alternate event location. The median noise levels during the experiment at different stations, waves generated by Betsy gunshots, and wavefields from nearby earthquakes point consistently to several structural units across the fault. Seismic trapping structure and local sedimentary basin produce localized motion amplification and stronger attenuation than adjacent regions. Cross correlations of high-frequency noise recorded at closely spaced stations provide a structural image of the subsurface material across the fault zone. The high spatial density and broad frequency range of the data can be used for additional high resolution studies of structure and source properties in the shallow crust.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 63
    Publication Date: 2015-05-16
    Description: Recently A. Gusev suggested and numerically investigated the doubly stochastic earthquake source model. The model is supposed to demonstrate the following features in the far-field body waves: (1) the omega-square high-frequency (HF) behaviour of displacement spectra; (2) lack of the directivity effect in HF radiation; and (3) a stochastic nature of the HF signal component. The model involves two stochastic elements: the local stress drop (SD) on a fault and the rupture time function (RT) with a linear dominant component. The goal of the present study is to investigate the Gusev model theoretically and to find conditions for (1, 2) to be valid and stable relative to receiver site. The models with smooth elements SD, RT are insufficient for these purposes. Therefore, SD and RT are treated as realizations of stochastic fields of the fractal type. The local smoothness of such fields is characterized by the fractional (Hurst) exponent H , 0 〈 H 〈 1. This allows us to consider a wide class of stochastic functions without regard to their global spectral properties. We show that the omega-square behavior of the model is achieved approximately if the rupture time function is almost regular ( H ~1) while the stress drop is rough function of any index H . However, if the rupture front is linear, the local stress drop has to be function of minimal smoothness ( H ~0). The situation with the directivity effect is more complicated: for different RT models with the same fractal index, the effect may or may not occur. The nature of the phenomenon is purely analytical. The main controlling factor for the directivity is the degree of smoothness of the 2-D distributions of RT random function. For this reason the directivity effect is unstable. This means that in practice the opposite conclusions relative to the statistical significance of the directivity effect are possible.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 64
    Publication Date: 2015-05-16
    Description: In general, seismic slip along faults reduces the average shear stress within earthquake source regions, but stress drops of specific earthquakes are observed to vary widely in size. To advance our understanding of variations in stress drop, we analysed source parameters of small-magnitude events in the greater San Gorgonio area, southern California. In San Gorgonio, the regional tectonics are controlled by a restraining bend of the San Andreas fault system, which results in distributed crustal deformation, and heterogeneous slip along numerous strike-slip and thrust faults. Stress drops were estimated by fitting a Brune-type spectral model to source spectra obtained by iteratively stacking the observed amplitude spectra. The estimates have large scatter among individual events but the median of event populations shows systematic, statistically significant variations. We identified several crustal and faulting parameters that may contribute to local variations in stress drop including the style of faulting, changes in average tectonic slip rates, mineralogical composition of the host rocks, as well as the hypocentral depths of seismic events. We observed anomalously high stress drops (〉20 MPa) in a small region between the traces of the San Gorgonio and Mission Creek segments of the San Andreas fault. Furthermore, the estimated stress drops are higher below depths of ~10 km and along the San Gorgonio fault segment, but are lower both to the north and south away from San Gorgonio Pass, showing an approximate negative correlation with geologic slip rates. Documenting controlling parameters of stress-drop heterogeneity is important to advance regional hazard assessment and our understanding of earthquake rupture processes.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 65
    Publication Date: 2014-12-13
    Description: Current earthquake early warning (EEW) systems lack the ability to appropriately handle multiple concurrent earthquakes, which led to many false alarms during the 2011 Tohoku earthquake sequence in Japan. This paper uses a Bayesian probabilistic approach to handle multiple concurrent events for EEW. We implement the theory using a two-step algorithm. First, an efficient approximate Bayesian model class selection scheme is used to estimate the number of concurrent events. Then, the Rao-Blackwellized Importance Sampling method with a sequential proposal probability density function is used to estimate the earthquake parameters, that is hypocentre location, origin time, magnitude and local seismic intensity. A real data example based on 2 months data (2011 March 9–April 30) around the time of the 2011 M 9 Tohoku earthquake is studied to verify the proposed algorithm. Our algorithm results in over 90 per cent reduction in the number of incorrect warnings compared to the existing EEW system operating in Japan.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 66
    Publication Date: 2014-11-07
    Description: Fractured rocks are known to exhibit seismic anisotropy and shear wave splitting (SWS). SWS is commonly used for fractured rock characterization and has been shown to be sensitive to fluid type. The presence of partial liquid/gas saturation is also known to affect the elastic properties of rocks. The combined effect of both fractures and partial liquid/gas saturation is still unknown. Using synthetic, silica-cemented sandstones with aligned penny-shaped voids, we conducted laboratory ultrasonic experiments to investigate the effect fractures aligned at an oblique angle to wave propagation would have on SWS under partial liquid/gas saturation conditions. The result for the fractured rock shows a saturation dependence which can be explained by combining a fractured rock model and a partial saturation model. At high to full water saturation values, SWS decreases as a result of the fluid bulk modulus effect on the quasi-shear wave. This bulk modulus effect is frequency dependent as a result of wave-induced fluid flow mechanisms, which would in turn lead to frequency dependent SWS. This result suggests the possible use of SWS for discriminating between full liquid saturation and partial liquid/gas saturation.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 67
    Publication Date: 2014-11-07
    Description: The effect of network density and geometric distribution on kinematic non-linear source inversion is investigated by inverting synthetic ground motions from a buried strike-slip fault ( M w 6.5), that have been generated by dynamic spontaneous rupture modelling. For the inversion, we use a physics-based regularized Yoffe function as slip velocity function. We test three different cases of station network geometry: (i) single station, varying azimuth and epicentral distance; (ii) multistation circular configurations, that is stations at similar distances from the fault, and regularly spaced around the fault; (iii) irregular multistation configurations using different numbers of stations. Our results show: (1) single station tests suggest that it may be possible to obtain a relatively good source model even using a single station. The best source model using a single station is obtained with stations at which amplitude ratios between three components are not large. We infer that both azimuthal angle and source-to-station distance play an important role in the design of optimal seismic network for source inversion. (2) Multistation tests show that the quality of the inverted source systematically correlates neither with the number of stations, nor with waveform misfit. (3) Waveform misfit has a direct correlation with the number of stations, resulting in overfitting the observed data without any systematic improvement of the source. It suggests that the best source model is not necessarily derived from the model with minimum waveform misfit. (4) A seismic network with a small number of well-spaced stations around the fault may be sufficient to obtain acceptable source inversion.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 68
    Publication Date: 2015-04-04
    Description: We analyze the cross-correlation function (CCF) of coda of earthquakes, which is used to retrieve the Green's function between two stations as well as the CCF of ambient noise. We select 74 Hi-net stations located in eastern Japan and 66 earthquakes to calculate the CCF. For each earthquake, we calculate the CCFs between possible pairs for the frequency bands of 0.1–0.2 Hz, 0.2–0.4 Hz and 0.4–0.8 Hz. Then we stack the CCFs for different earthquakes at each pair to obtain the average CCF. Although the correlation coefficients between the average and each CCFs are lower than 0.5 for most of the earthquakes, we obtain the propagating Rayleigh wave trace from average CCFs. We focus on the ratio of the amplitude in the positive lag time of the CCF to that in the negative lag time. CCFs for different earthquakes show different ratios which depend on the angle between the path of two stations and the epicentre. The amplitude in the lag time corresponding to the signal travelling from the near source station to the far source station is larger than that in the opposite lag time. Therefore the energy flux is not isotropic even in the coda and the energy from the source side is dominant. We average the ratios of pairs whose absolute values of angles are less than 45°. The average ratios are 0.5 at 0.1–0.2 Hz. For higher frequencies, the ratio is not clear because of the bad signal-to-noise ratio. According to the diffusion model, the ratio is predicted as 0.6. Therefore, the coda is represented as the diffusion state in 0.1–0.2 Hz with our observation setting.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 69
    Publication Date: 2015-04-04
    Description: Microseismic monitoring from large arrays using migration-based detection and location techniques is limited by detections of false positive events, which are the interpretation of spurious/noisy signals as real events. Therefore, semblance has been considered to differentiate between false positive and true events. However, semblance by itself is not suitable for variable signals such as those caused by shear source radiation. We present a new methodology for event detection and location using semblance of amplitudes corrected by a source mechanism. Our method is suitable for multichannel processing of microseismic data sets acquired with large arrays. The amplitudes are corrected by the radiation pattern of the inverted source mechanism before the semblance computation. We show that the source mechanism correction is the key factor in maximizing the value of semblance and makes the detection based on semblance superior to simple stacking. We apply this method to a data set recorded by a large surface star-like array on synthetic as well as on field data.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 70
    Publication Date: 2015-03-27
    Description: A new effective model is presented to compute horizontal-to-vertical spectral ratios (HVSR) relative to ambient vibrations, under the assumption that these are originated by a distribution of spatially correlated random surface sources. The major novelty of this model lies in the description of both ground displacement and sources as stochastic fields defined on the Earth's surface, stationary in time and homogeneous in space. In this frame, the power spectral density of the displacement stochastic field can be written as a function of the power spectral density of the force stochastic field and of the subsoil properties, through the relevant Green's function. Spatial correlation between ambient vibration sources is shown to be a necessary condition to warrant convergence of the integrals defining the frequency power spectra of the displacement field that make up the HVSR curve. Furthermore, it is shown that this HVSR curve may be significantly affected by the effective range of the force-field correlation on the Earth's surface. This formalization reduces computational efforts with respect to the previous version of the model based on distributed surface sources and may provide synthetic HVSR-curve patterns that are in line with those given by that computationally more troublesome version, as well as with those deduced under the assumption that the ambient vibrations constitute a diffuse wavefield.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 71
    Publication Date: 2015-02-06
    Description: Since centroid moment tensor (CMT) solution is not enough to describe the mechanism variations of large earthquakes, multiple CMTs were investigated in recent years by inverting seismic waveforms to retrieve the mechanism complexities. However, a problem in these waveform inversions is the trade-off between spatial mechanism variation and temporal rupture process, which may lead to ambiguous results. In this work, we propose a simple but practical approach to invert geodetic coseismic data for spatial mechanism variation. Because there is no need to solve for rupture velocity, geodetic inversion is generally more efficient and robust than waveform inversion. This approach was used to investigate the mechanism variation of the 2008 Wenchuan earthquake. In the case study, a line fault, whose position had been optimized through a 3-D grid search, was used to approximate the rupture fault, aiming to determine the mechanism variations along the strike direction. From the inversion results, there are thrust-slips and strike-slips in the southwest and northeast, respectively. The mechanism change appears at 160–200 km to the northeast of the epicentre, where compensated linear vector dipoles (CLVDs) have a local maximum and also support a significant mechanism change. By testing eight line faults with different positions, it is found that the inversion results do not strongly depend on the fault position, suggesting that the mechanism and its variations are robustly constrained in coseismic data inversion with our approach. Compared to waveform inversions, this work provides a more effective way to explore the mechanism complexities, and lays a foundation for future joint inversion of seismic and geodetic data sets.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 72
    Publication Date: 2015-02-06
    Description: The seismic structure of Earth's inner core is highly complex, displaying strong anisotropy and further regional variations. However, few seismic waves are sensitive to the inner core and fundamental questions regarding the origin of the observed seismic features remain unanswered. Thus, new techniques to observe different types of inner core waves are imperative to improve data coverage. Here, we detail our method for detecting exotic inner core phases such as PKJKP and PKIIKP, using inner core compressional waves as proof of concept. We use phase weighted stacking on long period data from a global distribution of stations, and employ several synthetic methods, including normal mode summation and SPECFEM, to identify and confirm the inner core phases. We present evidence for two observations of exotic inner core compressional waves, and apply the technique to a previously detected inner core shear wave. A possible new inner core shear wave remains unconfirmed. Additionally, we show how our method is important for rejecting potential observations, and distinguishing between waves with similar traveltime and slowness. The method is most successful for detecting exotic inner core compressional waves, and will provide a new approach for studying the compressional wave structures in the upper inner core.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 73
    Publication Date: 2015-03-11
    Description: Temporal changes in seismic anisotropy can be interpreted as variations in the orientation of cracks in seismogenic zones, and thus as variations in the stress field. Such temporal changes have been observed in seismogenic zones before and after earthquakes, although they are still not well understood. In this study, we investigate the azimuthal polarization of surface waves in anisotropic media with respect to the orientation of anisotropy, from a numerical point of view. This technique is based on the observation of the signature of anisotropy on the nine-component cross-correlation tensor (CCT) computed from seismic ambient noise recorded on pairs of three-component sensors. If noise sources are spatially distributed in a homogeneous medium, the CCT allows the reconstruction of the surface wave Green's tensor between the station pairs. In homogeneous, isotropic medium, four off-diagonal terms of the surface wave Green's tensor are null, but not in anisotropic medium. This technique is applied to three-component synthetic seismograms computed in a transversely isotropic medium with a horizontal symmetry axis, using a spectral element code. The CCT is computed between each pair of stations and then rotated, to approximate the surface wave Green's tensor by minimizing the off-diagonal components. This procedure allows the calculation of the azimuthal variation of quasi-Rayleigh and quasi-Love waves. In an anisotropic medium, in some cases, the azimuth of seismic anisotropy can induce a large variation in the horizontal polarization of surface waves. This variation depends on the relative angle between a pair of stations and the direction of anisotropy, the amplitude of the anisotropy, the frequency band of the signal and the depth of the anisotropic layer.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 74
    Publication Date: 2015-03-11
    Description: Conventional traveltime seismic tomography methods with Tikhonov regularization (L 2 norm) typically produce smooth models, but these models may be inappropriate when subsurface structure contains discontinuous features, such as faults or fractures, indicating that tomographic models should contain sharp boundaries. For this reason, we develop a double-difference (DD) traveltime tomography method that uses a modified total-variation regularization scheme incorporated with a priori information on interfaces to preserve sharp property contrasts and obtain accurate inversion results. In order to solve the inversion problem, we employ an alternating minimization method to decouple the original DD tomography problem into two separate subproblems: a conventional DD tomography with Tikhonov regularization and a L 2 total-variation inversion. We use the LSQR linear solver to solve the Tikhonov inversion and the split-Bregman iterative method to solve the total-variation inversion. Through our numerical examples, we show that our new DD tomography method yields more accurate results than the conventional DD tomography method at almost the same computational cost.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 75
    Publication Date: 2015-03-11
    Description: The Rudbar-Tarom earthquake of 1990 June 20 ( M 0  = 1.4 10 27 dyn cm), the largest one in Iran over the past 35 yr, was accompanied by a small tsunami in the Caspian Sea, which produced run-ups of up to 2 m, and inundations of 1 km, reaffirming the existence of tsunami hazard along the Caspian shoreline, as suggested by historical reports. We present the results of a field survey, documenting the concentration of the effects of the tsunami along a section of coastline not exceeding 30 km in length. A hydrodynamic simulation using the earthquake dislocation as the source of the tsunami fails to reproduce both the amplitude of the waves, and especially their concentration between the cities of Kiashahr and Jafrood. Rather, we show that the model of an underwater landslide, presumably triggered by the earthquake, and taking place on the steep slopes of the continental shelf approximately 10 km offshore, can fit the principal inundation characteristics identified during the survey. We suggest that the occurrence of such underwater landslides should become a primary ingredient to the assessment of tsunami risk along the Southern shores of the Caspian Sea.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 76
    Publication Date: 2015-04-09
    Description: In the Umbria–Marche (Central Italy) region an important earthquake sequence occurred in 1997, characterized by nine earthquakes with magnitudes in the range between 5 and 6, that caused important damages and causalities. In the present paper, we separately estimate intrinsic and scattering Q –1 parameters, using the classical multiple lapse time window analysis (MLTWA) approach in the assumption of a half-space model. The results clearly show that the attenuation parameters Q i –1 and Q s –1 are frequency dependent. This estimate is compared with other attenuation studies carried out in the same area, and with all the other MLTWA estimates obtained till now in other tectonic environments in the Earth. The bias introduced by the half-space assumption is investigated through numerical solutions of the energy transport equation in the more realistic assumption of a heterogeneous crust overlying a transparent mantle, with a Moho located at a depth ranging between 35 and 45 km below the surface. The bias introduced by the half-space assumption is significant only at high frequency. We finally show how the attenuation estimates, calculated with different techniques, lead to different peak ground acceleration decay with distance relationships, using the well-known and well proven Boore's method. This last result indicates that care must be used in selecting the correct estimate of the attenuation parameters for seismic risk purposes. We also discuss the reason why MLTWA may be chosen among all the other available techniques, due to its intrinsic stability, to obtain the right attenuation parameters.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 77
    Publication Date: 2015-04-12
    Description: We present a new approach to measuring crustal thickness and bulk properties from teleseismic data. In contrast to the traditional H – k stacking approach, which involves receiver-function deconvolution followed by stacking along expected arrival-time curves, we eliminate the deconvolution step and generate synthetic transfer functions predicting the relationship between the vertical and radial components. Given a catalogue of precalculated transfer functions, we convolve the vertical component with each assumed transfer function and then calculate a misfit between the real and predicted radial component. As a single-layer crust is no longer a necessary assumption, we use the transfer-function approach to extend H – k analysis to models containing sedimentary basins. Applied to a data set in Minnesota, North Dakota and South Dakota, where the Archean Superior Province, Paleoproterozoic Trans-Hudson Orogen and Meso/Neoproterozoic Mid-Continent Rift are largely overlain by sediments of the Williston Basin, we find that the transfer-function approach is able to recover approximate sedimentary thickness, as well as remove the contaminating effects of the sedimentary layers from crustal thickness and basement P/S velocity ratio measurements. We find that the Superior Province has uniformly low P/S ratio, reflecting a highly felsic composition, but has substantial Moho topography, with substantial crustal thinning to the west possibly related to Trans-Hudson accretion. The Trans-Hudson has thick crust with higher (more mafic) P/S ratio; higher ratios in the Superior near the Trans-Hudson contact suggest that the Superior is overthrust on Trans-Hudson basement material. The vicinity of the Mid-Continent Rift shows thinned crust, elevated P/S ratio, and a localized zone of very high P/S ratio associated with a rift offset, possibly related to a former triple junction, that extends northwestward into the Superior following a recently detected linear mantle anomaly.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 78
    Publication Date: 2015-04-12
    Description: Measurements of seismic anisotropy are commonly used to constrain deformation in the upper mantle. Observations of anisotropy at mid-mantle depths are, however, relatively sparse. In this study we probe the anisotropic structure of the mid-mantle (transition zone and uppermost lower mantle) beneath the Japan, Izu-Bonin, and South America subduction systems. We present source-side shear wave splitting measurements for direct teleseismic S phases from earthquakes deeper than 300 km that have been corrected for the effects of upper mantle anisotropy beneath the receiver. In each region, we observe consistent splitting with delay times as large as 1 s, indicating the presence of anisotropy at mid-mantle depths. Clear splitting of phases originating from depths as great as ~600 km argues for a contribution from anisotropy in the uppermost lower mantle as well as the transition zone. Beneath Japan, fast splitting directions are perpendicular or oblique to the slab strike and do not appear to depend on the propagation direction of the waves. Beneath South America and Izu-Bonin, splitting directions vary from trench-parallel to trench-perpendicular and have an azimuthal dependence, indicating lateral heterogeneity. Our results provide evidence for the presence of laterally variable anisotropy and are indicative of variable deformation and dynamics at mid-mantle depths in the vicinity of subducting slabs.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 79
    Publication Date: 2015-04-12
    Description: We estimated spatiotemporal slip distributions from three long-term slow slip events (L-SSEs) that occurred beneath the Bungo Channel at the convergent plate boundary between the subducting oceanic Philippine Sea plate and the continental Amurian plate in southwest Japan between 1997 and 1998, 2002 and 2004 and 2009 and 2011. For this purpose, we employed an inversion method using a Bayesian Information Criterion (ABIC), which included the following three prior constraints: the spatial slip distribution was smooth to some extent, slip directions were mostly oriented in the direction of plate convergence and the temporal change in slip was smooth to some extent. Our results revealed that the three L-SSEs had a common feature: slipped regions expanded southwestward at accelerating slip velocities. We also found that major slipped regions migrated southwestward by approximately 50–100 km yr –1 . In contrast, southwestward and northeastward migration of the slipped regions, whose direction differed from event to event, was also identified before or after the periods when the slip velocities were at their greatest. Comparing the obtained spatiotemporal slip distributions of the three L-SSEs with slip-deficit rate distributions obtained in our previous study, we investigated the accumulation process of the slip deficit caused by slip-deficit rate distributions and the release processes of the slip deficit caused by the obtained spatiotemporal slip distributions of the three L-SSEs. At the western plate interface of the Bungo Channel, as the slip-deficit rate was small and the amounts of slips associated with the three L-SSEs were large, most of the accumulated slip deficit was estimated to have been released. In contrast, at the eastern plate interface, as the slip-deficit rate was large and the amounts of slips associated with the three L-SSEs were small, the slip deficit was estimated to have accumulated effectively. These results suggest that the slipped regions of the three L-SSEs and the strongly coupled region are not spatially complementary; the accumulated slip deficit showed spatial variation even at approximately the same depth range along the arc.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 80
    Publication Date: 2015-04-12
    Description: We present observations of cross coupled spheroidal modes in the Earth's free oscillation spectrum recorded by the vertical component G-ring laser (Geodetical Station Wettzell) of the 2011 M w 9.0 Tohoku-Oki earthquake. In an attempt to determine which are the mechanisms responsible for spheroidal energy in a vertical axes rotational spectra, we first rule out instrumental effects as well as the effect of local heterogeneity. Secondly, we carry out simulations of an ideal rotational sensor taking into account the effects of the Earth's rotation, its hydrostatic ellipticity and structural heterogeneity, which results in a good fit to the data. Simulations considering each effect separately are performed in order to evaluate the sensitivity of rotational motions to various phenomena compared to traditional translation measurements.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 81
    Publication Date: 2015-04-16
    Description: Local earthquake tomography is a non-linear and non-unique inverse problem that uses event arrival times to solve for the spatial distribution of elastic properties. The typical approach is to apply iterative linearization and derive a preferred solution, but such solutions are biased by a number of subjective choices: the starting model that is iteratively adjusted, the degree of regularization used to obtain a smooth solution, and the assumed noise level in the arrival time data. These subjective choices also affect the estimation of the uncertainties in the inverted parameters. The method presented here is developed in a Bayesian framework where a priori information and measurements are combined to define a posterior probability density of the parameters of interest: elastic properties in a subsurface 3-D model, hypocentre coordinates and noise level in the data. We apply a trans-dimensional Markov chain Monte Carlo algorithm that asymptotically samples the posterior distribution of the investigated parameters. This approach allows us to overcome the issues raised above. First, starting a number of sampling chains from random samples of the prior probability distribution lessens the dependence of the solution from the starting point. Secondly, the number of elastic parameters in the 3-D subsurface model is one of the unknowns in the inversion, and the parsimony of Bayesian inference ensures that the degree of detail in the solution is controlled by the information in the data, given realistic assumptions for the error statistics. Finally, the noise level in the data, which controls the uncertainties of the solution, is also one of the inverted parameters, providing a first-order estimate of the data errors. We apply our method to both synthetic and field arrival time data. The synthetic data inversion successfully recovers velocity anomalies, hypocentre coordinates and the level of noise in the data. The Bayesian inversion of field measurements gives results comparable to those obtained independently by linearized inversion, reconstructing the geometry of the main seismic velocity anomalies. The quantification of the posterior uncertainties, a crucial output of Bayesian inversion, allows for visualizing regions where elastic properties are closely constrained by the data and is used here to directly compare our results to the ones obtained with the linearized inversion. In the case we examined the results of two inversion techniques are not significantly different.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 82
    Publication Date: 2015-04-16
    Description: A concept of seismic system (SS), which is responsible for the preparation of an ensemble of strong earthquakes, is considered as an open dissipative system exchanging energy and entropy with the environment. Open dissipative SS allow one to describe the equilibrium and non-equilibrium states of SS, and the lithosphere evolution under different plate tectonic settings on the basis of seismostatistics. Several new seismic parameters (‘seismic temperature’, ‘seismic time’, dissipation function, efficiency, inelastic energy, dynamical probability) are defined and proposed for better understanding and describing the dynamical processes. The Sakhalin SS is considered to illustrate the behaviour of proposed parameters. By analogy to Liouville's equation in thermodynamics, it is shown that there is no criterion of instability in the domain where the Gutenberg-Richter law is true. In the proposed approach, the instability origination and the formation of seismogenic structures in the lithosphere are based on the energy versus information entropy power law; the existence of ‘time arrow’ also proceeds from such a dependence. Application of energy and trajectory diagrams enables to describe the preparation of strong earthquakes within an ensemble in terms of slow and fast timescales. These diagrams help perform the spatiotemporal-energy monitoring of the instability origination in the lithosphere. It is shown that the information entropy parameter can serve as a measure of the unknown external energy flow into the system (this energy is supplied for the elastic radiation energy in the earthquake sources and for inelastic processes in the system volume). The property of the ensemble of strong earthquakes is periodically to restore the SS equilibrium state that enables to describe the SS energy balance. The results offer possibilities to estimate the fraction of inelastic energy released by the SS medium during the preparation and occurrence of seismic catastrophes. The construction of phase diagrams and ‘entropy funnels’ in the virtual space (with time, energy and entropy coordinates) can provide new opportunities for the visualization of undetectable processes leading to disastrous earthquakes.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 83
    Publication Date: 2015-04-16
    Description: The structure of the Earth's inner core is not well known between depths of ~100–200 km beneath the inner core boundary. This is a result of the PKP core phase triplication and the existence of strong precursors to PKP phases, which hinder the measurement of inner core compressional PKIKP waves at epicentral distances between roughly 143 and 148°. Consequently, interpretation of the detailed structure of deeper regions also remains difficult. To overcome these issues we stack seismograms in slowness and time, separating the PKP and PKIKP phases which arrive simultaneously but with different slowness. We apply this method to study the inner core's Western hemisphere beneath South and Central America using paths travelling in the quasi-polar direction between 140 and 150° epicentral distance, which enables us to measure PKiKP–PKIKP differential traveltimes up to greater epicentral distance than has previously been done. The resulting PKiKP–PKIKP differential traveltime residuals increase with epicentral distance, which indicates a marked increase in seismic velocity for polar paths at depths greater than 100 km compared to reference model AK135. Assuming a homogeneous outer core, these findings can be explained by either (i) inner core heterogeneity due to an increase in isotropic velocity or (ii) increase in anisotropy over the studied depth range. Although this study only samples a small region of the inner core and the current data cannot distinguish between the two alternatives, we prefer the latter interpretation in the light of previous work.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 84
    Publication Date: 2015-04-16
    Description: The PP precursors are seismic waves that form from underside reflections of P waves off discontinuities in the upper mantle transition zone (MTZ). These seismic phases are used to map discontinuity topography, sharpness, and impedance contrasts; the resulting structural variations are then often interpreted as evidence for temperature and/or mineralogy variations within the mantle. The PP precursors as well as other seismic phases have been used to establish the global presence of seismic discontinuities at 410 and 660 km depth. Intriguingly, in more than 80 per cent of PP precursor observations the seismic wave amplitudes are significantly weaker than the amplitudes predicted by seismic reference models. Even more perplexing is the observation that 1–5 per cent of all earthquakes (which are 20–25 per cent of earthquakes with clear PP waveforms) do not show any evidence for the PP precursors from the discontinuities even in the presence of well-developed PP waveforms. Non-detections are found in six different data sets consisting of tens to hundreds of events. We use synthetic modelling to examine a suite of factors that could be responsible for the absence of the PP precursors. The take-off angles for PP and the precursors differ by only 1.2–1.5°; thus source-related complexity would affect PP and the precursors. A PP wave attenuated in the upper mantle would increase the relative amplitude of the PP precursors. Attenuation within the transition zone could reduce precursor amplitudes, but this would be a regional phenomenon restricted to particular source receiver geometries. We also find little evidence for deviations from the theoretical travel path of seismic rays expected for scattered arrivals. Factors that have a strong influence include the stacking procedures used in seismic array techniques in the presence of large, interfering phases, the presence of topography on the discontinuities on the order of tens of kilometres, and 3-D lateral heterogeneity in the velocity and density changes with depth across the transition zone. We also compare the observed precursors’ amplitudes with seismic models from calculations of phase equilibria and find that a seismic velocity model derived from a pyrolite composition reproduces the data better than the currently available 1-D earth models. This largely owes to the pyrolite models producing a stronger minimum in the reflection coefficient across the epicentral distances where the reduction in amplitudes of the PP precursors is observed. To suppress the precursors entirely in a small subset of earthquakes, other effects, such as localized discontinuity topography and seismic signal processing effects are required in addition to the changed velocity model.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 85
    Publication Date: 2015-04-16
    Description: A method for combining the asymptotic operator designed by Beylkin (Born migration operator) for the solution of linearized inverse problems with full waveform inversion is presented. This operator is used to modify the standard L 2 norm that measures the distance between synthetic and observed data. The modified misfit function measures the discrepancy of the synthetic and observed data after they have been migrated using the Beylkin operator. The gradient of this new misfit function is equal to the cross-correlation of the single scattering data with migrated/demigrated residuals. The modified misfit function possesses a Hessian operator that tends asymptotically towards the identity operator. The trade-offs between discrete parameters are thus reduced in this inversion scheme. Results on 2-D synthetic case studies demonstrate the fast convergence of this inversion method in a migration regime. From an accurate estimation of the initial velocity, three and five iterations only are required to generate high-resolution P -wave velocity estimation models on the Marmousi 2 and synthetic Valhall case studies.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 86
    Publication Date: 2015-04-16
    Description: At 155 km, the Pärvie fault is the world's longest known endglacial fault (EGF). It is located in northernmost Sweden in a region where several kilometre-scale EGFs have been identified. Based on studies of Quaternary deposits, landslides and liquefaction structures, these faults are inferred to have ruptured as large earthquakes when the latest ice sheet disappeared from the region, some 9500 yr ago. The EGFs still exhibit relatively high seismic activity, and here we present new earthquake data from northern Sweden in general and the Pärvie fault in particular. More than 1450 earthquakes have been recorded in Sweden north of 66° latitude in the years 2000–2013. There is a remarkable correlation between this seismicity and the mapped EGF scarps. We find that 71 per cent of the observed earthquakes north of 66° locate within 30 km to the southeast and 10 km to the northwest of the EGFs, which is consistent with the EGFs’ observed reverse faulting mechanisms, with dips to the southeast. In order to further investigate the seismicity along the Pärvie fault we installed a temporary seismic network in the area between 2007 and 2010. In addition to the routine automatic detection and location algorithm, we devised a waveform cross-correlation technique which resulted in a 50 per cent increase of the catalogue and a total of 1046 events along the Pärvie fault system between 2003 and 2013. The earthquakes were used to establish an improved velocity model for the area, using 3-D local earthquake tomography. The resulting 3-D velocity model shows smooth, minor velocity variations in the area. All events were relocated in this new 3-D model. A tight cluster on the central part of the Pärvie fault, where the rate of seismicity is the highest, could be relocated with high precision relative location. We performed depth phase analysis on 40 of the larger events to further constrain the hypocentral locations. We find that the seismicity on the Pärvie fault correlates very well with the mapped surface trace of the fault. The events do not align along a well-defined fault plane at depth but form a zone of seismicity that dips between 30° and 60° to the southeast of the surface fault trace, with distinct along-strike variations. The seismic zone extends to approximately 35 km depth. Using this geometry and earthquake scaling relations, we estimate that the endglacial Pärvie earthquake had a magnitude of 8.0 ± 0.4.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 87
    Publication Date: 2015-04-17
    Description: The depth of an earthquake is difficult to estimate because of the trade-off between depth and origin time estimations, and because it can be biased by lateral Earth heterogeneities. To face this challenge, we have developed a new, blind and fully automatic teleseismic depth analysis. The results of this new method do not depend on epistemic uncertainties due to depth-phase picking and identification. The method consists of a modification of the cepstral analysis from Letort et al. and Bonner et al. , which aims to detect surface reflected (pP, sP) waves in a signal at teleseismic distances (30°–90°) through the study of the spectral holes in the shape of the signal spectrum. The ability of our automatic method to improve depth estimations is shown by relocation of the recent moderate seismicity of the Guerrero subduction area (Mexico). We have therefore estimated the depth of 152 events using teleseismic data from the IRIS stations and arrays. One advantage of this method is that it can be applied for single stations (from IRIS) as well as for classical arrays. In the Guerrero area, our new cepstral analysis efficiently clusters event locations and provides an improved view of the geometry of the subduction. Moreover, we have also validated our method through relocation of the same events using the new International Seismological Centre (ISC)-locator algorithm, as well as comparing our cepstral depths with the available Harvard–Centroid Moment Tensor (CMT) solutions and the three available ground thrust (GT5) events (where lateral localization is assumed to be well constrained with uncertainty 〈5 km) for this area. These comparisons indicate an overestimation of focal depths in the ISC catalogue for deeper parts of the subduction, and they show a systematic bias between the estimated cepstral depths and the ISC-locator depths. Using information from the CMT catalogue relating to the predominant focal mechanism for this area, this bias can be explained as a misidentification of sP phases by pP phases, which shows the greater interest for the use of this new automatic cepstral analysis, as it is less sensitive to phase identification.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 88
    Publication Date: 2015-04-19
    Description: Regional seismic networks may and in some cases need to be used to monitor teleseismic or network-outside events. For detecting and localizing teleseismic events automatically and reliably in this case, in this paper we present a novel progressive association algorithm for teleseismic signals recorded by a regional seismic network. The algorithm takes triangle station arrays as the starting point to search for P waves of teleseismic events progressively by that, as detections from different stations actually are from the same teleseismic event, their arrival times should be linearly related to the average slowness vector with which the signal propagates across the network, and the slowness of direct teleseismic P wave basically is different from other major seismic phases. We have tested this algorithm using data recorded by Xinjiang Seismic Network of China (XJSN) for 16 d. The results show that the algorithm can effectively and reliably detect and localize earthquakes outside of the network. For the period of the test data, as all mb 4.0+ events with c 〈 30° and all mb 4.5+ events with c 〈 60° referring to the International Data Center-Reviewed Event Bulletin (IDC REB) were detected, where c is the epicentral distance relative to the network's geographical centre, the rate of false events only accounted for 2.4 per cent, suggesting that the new association algorithm has good application prospect for situations when regional seismic networks need to be used to monitor teleseismic events.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 89
    Publication Date: 2016-04-03
    Description: Pressure solution creep (PSC) is an important elementary process in rock friction at high temperatures where solubilities of rock-forming minerals are significantly large. It significantly changes the frictional resistance and enhances time-dependent strengthening. A recent microphysical model for PSC-involved friction of clay–quartz mixtures, which can explain a transition between dilatant and non-dilatant deformation (d-nd transition), was modified here and implemented in dynamic earthquake sequence simulations. The original model resulted in essentially a kind of rate- and state-dependent friction (RSF) law, but assumed a constant friction coefficient for clay resulting in zero instantaneous rate dependency in the dilatant regime. In this study, an instantaneous rate dependency for the clay friction coefficient was introduced, consistent with experiments, resulting in a friction law suitable for earthquake sequence simulations. In addition, a term for time-dependent strengthening due to PSC was added which makes the friction law logarithmically rate-weakening in the dilatant regime. The width of the zone in which clasts overlap or, equivalently, the interface porosity involved in PSC plays a role as the state variable. Such a concrete physical meaning of the state variable is a great advantage in future modelling studies incorporating other physical processes such as hydraulic effects. Earthquake sequence simulations with different pore pressure distributions demonstrated that excess pore pressure at depth causes deeper rupture propagation with smaller slip per event and a shorter recurrence interval. The simulated ruptures were arrested a few kilometres below the point of pre-seismic peak stress at the d-nd transition and did not propagate spontaneously into the region of pre-seismic non-dilatant deformation. PSC weakens the fault against slow deformation and thus such a region cannot produce a dynamic stress drop. Dynamic rupture propagation further down to brittle-plastic transition, evidenced by geological observations, would require even smaller frictional resistance at coseismic slip rate, suggesting the importance of implementation of dynamic weakening activated at coseismic slip rates for more realistic simulation of earthquake sequences. The present models produced much smaller afterslip at deeper parts of arrested ruptures than those with logarithmic RSF laws because of a more significant rate-strengthening effect due to linearly viscous PSC. Detailed investigation of afterslip would give a clue to understand the deformation mechanism which controls shear resistance of the fault in a region of arrest of earthquake ruptures.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 90
    Publication Date: 2016-04-03
    Description: Surface wave tomography routinely uses empirically scaled density model in the inversion of dispersion curves for shear wave speeds of the crust and uppermost mantle. An improperly selected empirical scaling relationship between density and shear wave speed can lead to unrealistic density models beneath certain tectonic formations such as sedimentary basins. Taking the Sichuan basin east to the Tibetan plateau as an example, we investigate the differences between density profiles calculated from four scaling methods and their effects on Rayleigh wave phase velocities. Analytical equations for 1-D layered models and adjoint tomography for 3-D models are used to examine the trade-off between density and S -wave velocity structures at different depth ranges. We demonstrate that shallow density structure can significantly influence phase velocities at short periods, and thereby affect the shear wave speed inversion from phase velocity data. In particular, a deviation of 25 per cent in the initial density model can introduce an error up to 5 per cent in the inverted shear velocity at middle and lower crustal depths. Therefore one must pay enough attention in choosing a proper velocity–density scaling relationship in constructing initial density model in Rayleigh wave inversion for crustal shear velocity structure.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 91
    Publication Date: 2016-04-07
    Description: Application of the ambient noise surface wave tomography method (ANT) for determination of the upper-mantle structure requires data on long-periodic noise ( T 〉 40 s). The ANT technique implies that noise sources are distributed almost uniformly over the surface. This is practically true for short-periodic noise, however, it is not so in the case of long periods. In this paper we show that the main contribution to noise at long periods is caused by signals from earthquakes. In some cases, they may strongly distort noise cross-correlation. This leads to an incorrect determination of surface wave velocity dispersion curves. To minimize such a distortion we propose two means: (1) to use records of noise for the periods when there is no clustering of earthquakes, such as aftershocks of strong events; (2) to stack cross-correlation functions for a period of at least three years in order to achieve sufficient uniformity of earthquake locations. Validity of this approach is demonstrated by ANT results for Europe. Tomographic reconstruction of Rayleigh wave group velocities for 10–100 s measured along interstation paths was carried out in a central part of Western Europe where resolving power of the data was the highest. Locally averaged dispersion curves were inverted to vertical S -wave velocity sections in this area. The results correspond closely to known features of the structure of the region, namely: strong difference of the crust and upper-mantle structure at the opposite sides from the Tornquist–Teisseyre Line down to ~ 250 km, penetration of high-velocity material of East European Platform lithosphere under Carpathians, as well as penetration of low-velocity asthenospheric layer from the Carpathian region towards the northeast.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 92
    Publication Date: 2016-04-07
    Description: Sensitivity analysis with synthetic models is widely used in seismic tomography as a means for assessing the spatial resolution of solutions produced by, in most cases, linear or iterative nonlinear inversion schemes. The most common type of synthetic reconstruction test is the so-called checkerboard resolution test in which the synthetic model comprises an alternating pattern of higher and lower wave speed (or some other seismic property such as attenuation) in 2-D or 3-D. Although originally introduced for application to large inverse problems for which formal resolution and covariance could not be computed, these tests have achieved popularity, even when resolution and covariance can be computed, by virtue of being simple to implement and providing rapid and intuitive insight into the reliability of the recovered model. However, checkerboard tests have a number of potential drawbacks, including (1) only providing indirect evidence of quantitative measures of reliability such as resolution and uncertainty, (2) giving a potentially misleading impression of the range of scale-lengths that can be resolved, and (3) not giving a true picture of the structural distortion or smearing that can be caused by the data coverage. The widespread use of synthetic reconstruction tests in seismic tomography is likely to continue for some time yet, so it is important to implement best practice where possible. The goal of this paper is to develop the underlying theory and carry out a series of numerical experiments in order to establish best practice and identify some common pitfalls. Based on our findings, we recommend (1) the use of a discrete spike test involving a sparse distribution of spikes, rather than the use of the conventional tightly spaced checkerboard; (2) using data coverage (e.g. ray-path geometry) inherited from the model constrained by the observations (i.e. the same forward operator or matrix), rather than the data coverage obtained by solving the forward problem through the synthetic model; (3) carrying out multiple tests using structures of different scale length; (4) taking special care with regard to what can be inferred when using synthetic structures that closely mimic what has been recovered in the observation-based model; (5) investigating the range of structural wavelengths that can be recovered using realistic levels of imposed data noise; and (6) where feasible, assessing the influence of model parametrization error, which arises from making a choice as to how structure is to be represented.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 93
    Publication Date: 2016-04-07
    Description: We develop an approach for simulating acousto-elastic wave phenomena, including scattering from fluid–solid boundaries, where the solid is allowed to be anisotropic, with the discontinuous Galerkin method. We use a coupled first-order elastic strain-velocity, acoustic velocity–pressure formulation, and append penalty terms based on interior boundary continuity conditions to the numerical (central) flux so that the consistency condition holds for the discretized discontinuous Galerkin weak formulation. We incorporate the fluid–solid boundaries through these penalty terms and obtain a stable algorithm. Our approach avoids the diagonalization into polarized wave constituents such as in the approach based on solving elementwise Riemann problems.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 94
    Publication Date: 2016-04-07
    Description: It is well known that large earthquakes generally trigger aftershock sequences. However, the duration of those sequences is unclear due to the gradual power-law decay with time. The triggering time is assumed to be infinite in the epidemic type aftershock sequence (ETAS) model, a widely used statistical model to describe clustering phenomena in observed earthquake catalogues. This assumption leads to the constraint that the power-law exponent p of the Omori-Utsu decay has to be larger than one to avoid supercritical conditions with accelerating seismic activity on long timescales. In contrast, seismicity models based on rate- and state-dependent friction observed in laboratory experiments predict p ≤ 1 and a finite triggering time scaling inversely to the tectonic stressing rate. To investigate this conflict, we analyse an ETAS model with finite triggering times, which allow smaller values of p . We use synthetic earthquake sequences to show that the assumption of infinite triggering times can lead to a significant bias in the maximum likelihood estimates of the ETAS parameters. Furthermore, it is shown that the triggering time can be reasonably estimated using real earthquake catalogue data, although the uncertainties are large. The analysis of real earthquake catalogues indicates mainly finite triggering times in the order of 100 days to 10 years with a weak negative correlation to the background rate, in agreement with expectations of the rate- and state-friction model. The triggering time is not the same as the apparent duration, which is the time period in which aftershocks dominate the seismicity. The apparent duration is shown to be strongly dependent on the mainshock magnitude and the level of background activity. It can be much shorter than the triggering time. Finally, we perform forward simulations to estimate the effective forecasting period, which is the time period following a mainshock, in which ETAS simulations can improve rate estimates after the occurrence of a mainshock. We find that this effective forecasting period is only in the order of 100 days for moderate mainshocks and in the order of a few years for large events, even if the underlying triggering process lasts much longer.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 95
    Publication Date: 2016-01-01
    Description: We present a new upper-mantle tomographic model derived solely from hum seismic data. Phase correlograms between station pairs are computed to extract phase-coherent signals. Correlograms are then stacked using the time–frequency phase-weighted stack method to build-up empirical Green's functions. Group velocities and uncertainties are measured in the wide period band of 30–250 s, following a resampling approach. Less data are required to extract reliable group velocities at short periods than at long periods, and 2 yr of data are necessary to measure reliable group velocities for the entire period band. Group velocities are first regionalized and then inverted versus depth using a simulated annealing method in which the number and shape of splines that describes the S -wave velocity model are variable. The new S -wave velocity tomographic model is well correlated with models derived from earthquakes in most areas, although in India, the Dharwar craton is shallower than in other published models.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 96
    Publication Date: 2016-01-01
    Description: Accurate prediction of ground motion intensity and its variability is an important element in seismic hazard assessment. Simulation-based ground motion prediction has become popular in light of earthquake rupture and wave-propagation modelling methods and the rapid growth in computing power. Earthquake rupture modelling needs to be physics-based and also computationally efficient. It also requires the ability to quantify the variability of finite source models for future scenario events. A generalized pseudo-dynamic source model (e.g. earthquake source statistics model) for M w 6.5–7.0 has been constructed by analysing a number of dynamic rupture models based on 1-point and 2-point statistics of kinematic source parameters. Synthetic broad-band ground motions derived from the pseudo-dynamic source model are validated against empirical ground motion prediction equations (GMPEs). The pseudo-dynamic source model produces ground motion intensities mostly compatible with the empirical GMPEs in the broad frequency range. The perturbation analysis of 1-point and 2-point statistics helps to elucidate the effect of source statistics on ground motions in a systematic way. The constructed pseudo-dynamic source model may be used to generate a feasible range of rupture scenarios for future events and to simulate expected ground motions for seismic hazard assessment.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 97
    Publication Date: 2016-01-01
    Description: We evaluate the bias in parameter estimates of the ETAS model. We show that when a simulated catalogue is magnitude-truncated there is considerable bias, whereas when it is not truncated there is no discernible bias. We also discuss two further implied assumptions in the ETAS and other self-exciting models. First, that the triggering boundary magnitude is equivalent to the catalogue completeness magnitude. Secondly, the assumption in the Gutenberg–Richter relationship that numbers of events increase exponentially as magnitude decreases. These two assumptions are confounded with the magnitude truncation effect. We discuss the effect of these problems on analyses of real earthquake catalogues.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 98
    Publication Date: 2016-01-02
    Description: The frequency dependence of the quality factor Q has long been predicted by mathematical modelling and laboratory measurements; however, in situ evidence from seismic surveys is still lacking. We have conducted the cross-hole seismic surveys to investigate the near-surface seismic attenuation in the Daqing oilfield in northeastern China. The seismic waves were fired in a source hole of 40 m from the bottom to the surface at an interval of 1 m and were recorded in a receiver hole of 40 m by two geophones with one at the surface and the other one at the bottom. The direct waves were extracted to avoid the noise disturbance and the reflection interference, and the attenuations without the effects of the source signature and the receiver coupling were estimated by a method we proposed. The nonlinear attenuations were observed and fitted using the power-law-based Q . The reliability of Q estimate was verified by the high similarity between the real and the simulated attenuations. Therefore, the experiment we have conducted can be treated as a reliable evidence for the frequency dependence of near-surface  Q .
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 99
    Publication Date: 2016-01-08
    Description: The southeastern European cities of Sofia and Thessaloniki are explored as example site-specific scenarios by geographically zoning their individual localized seismic sources based on the highest probabilities of magnitude exceedance. This is with the aim of determining the major components contributing to each city's seismic hazard. Discrete contributions from the selected input earthquake catalogue are investigated to determine those areas that dominate each city's prevailing seismic hazard with respect to magnitude and source-to-site distance. This work is based on an earthquake catalogue developed and described in a previously published paper by the author and components of a magnitude probability density function. Binned magnitude and distance classes are defined using a joint magnitude–distance distribution. The prevailing seismicity to each city—as defined by a child data set extracted from the parent earthquake catalogue for each city considered—is divided into distinct constrained data bins of small discrete magnitude and source-to-site distance intervals. These are then used to describe seismic hazard in terms of uni-variate modal values; that is, M * and D * which are the modal magnitude and modal source-to-site distance in each city's local historical seismicity. This work highlights that Sofia's dominating seismic hazard—that is, the modal magnitudes possessing the highest probabilities of occurrence—is located in zones confined to two regions at 60–80 km and 170–180 km from this city, for magnitude intervals of 5.75–6.00 M w and 6.00–6.25 M w respectively. Similarly, Thessaloniki appears prone to highest levels of hazard over a wider epicentral distance interval, from 80 to 200 km in the moment magnitude range 6.00–6.25 M w .
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 100
    Publication Date: 2016-03-09
    Description: We collect two months of ambient noise data recorded by 35 broad-band seismic stations in a 9 x 11 km area (1–3 km station interval) near Karamay, China, and do cross-correlation of noise data between all station pairs. Array beamforming analysis of the ambient noise data shows that ambient noise sources are unevenly distributed and the most energetic ambient noise mainly comes from azimuths of 40°–70°. As a consequence of the strong directional noise sources, surface wave components of the cross-correlations at 1–5 Hz show clearly azimuthal dependence, and direct dispersion measurements from cross-correlations are strongly biased by the dominant noise energy. This bias renders that the dispersion measurements from cross-correlations do not accurately reflect the interstation velocities of surface waves propagating directly from one station to the other, that is, the cross-correlation functions do not retrieve empirical Green's functions accurately. To correct the bias caused by unevenly distributed noise sources, we adopt an iterative inversion procedure. The iterative inversion procedure, based on plane-wave modeling, includes three steps: (1) surface wave tomography, (2) estimation of ambient noise energy and biases and (3) phase velocities correction. First, we use synthesized data to test the efficiency and stability of the iterative procedure for both homogeneous and heterogeneous media. The testing results show that: (1) the amplitudes of phase velocity bias caused by directional noise sources are significant, reaching ~2 and ~10 per cent for homogeneous and heterogeneous media, respectively; (2) phase velocity bias can be corrected by the iterative inversion procedure and the convergence of inversion depends on the starting phase velocity map and the complexity of the media. By applying the iterative approach to the real data in Karamay, we further show that phase velocity maps converge after 10 iterations and the phase velocity maps obtained using corrected interstation dispersion measurements are more consistent with results from geology surveys than those based on uncorrected data. As ambient noise in high-frequency band (〉1 Hz) is mostly related to human activities or climate events, both of which have strong directivity, the iterative approach demonstrated here helps improve the accuracy and resolution of ANT in imaging shallow earth structures.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
Close ⊗
This website uses cookies and the analysis tool Matomo. More information can be found here...