To save content items to your account,
please confirm that you agree to abide by our usage policies.
If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your account.
Find out more about saving content to .
To save content items to your Kindle, first ensure no-reply@cambridge.org
is added to your Approved Personal Document E-mail List under your Personal Document Settings
on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part
of your Kindle email address below.
Find out more about saving to your Kindle.
Note you can select to save to either the @free.kindle.com or @kindle.com variations.
‘@free.kindle.com’ emails are free but can only be saved to your device when it is connected to wi-fi.
‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.
Among all the numerical methods in seismology, the finite-difference (FD) technique provides the best balance of accuracy and computational efficiency. This book offers a comprehensive introduction to FD and its applications to earthquake motion. Using a systematic tutorial approach, the book requires only undergraduate degree-level mathematics and provides a user-friendly explanation of the relevant theory. It explains FD schemes for solving wave equations and elastodynamic equations of motion in heterogeneous media, and provides an introduction to the rheology of viscoelastic and elastoplastic media. It also presents an advanced FD time-domain method for efficient numerical simulations of earthquake ground motion in realistic complex models of local surface sedimentary structures. Accompanied by a suite of online resources to help put the theory into practice, this is a vital resource for professionals and academic researchers using numerical seismological techniques, and graduate students in earthquake seismology, computational and numerical modelling, and applied mathematics.
Intraplate earthquakes occur away from tectonic plate boundaries: their locations are difficult to predict, risking huge damage and loss of life. The 2001 Bhuj earthquake (featured in this book) was the largest intraplate earthquake for three decades and has provided unique insight into these events. This cutting-edge book brings together research from international leading experts in the field. Each chapter provides a comprehensive review of these earthquakes in a different global location, ranging from Australia, China, India and the Sea of Japan, to Western Europe, Brazil, New Madrid (Central USA), and Eastern Canada. They explore similarities and differences between regional features and the mechanical models required to explain them, as well as assessing geophysical techniques used to investigate them. Providing the first global overview of intraplate earthquakes, this is an essential book for academic researchers and professionals in seismology, tectonics, tectonophysics, geodesy, structural geology, earthquake dynamics, geophysics, and structural engineering.
The Neoproterozoic Hedmark Basin in the Caledonides of South Norway was formed at the western margin of the continent Baltica by rifting 750–600 Ma ago. The margin was destroyed in the Caledonian Orogeny and sedimentary basins translated eastwards. This study uses provenance analysis to map the crustal architecture of the pre-Caledonian SW Baltican margin. Conglomerate clasts and sandstones were sampled from submarine fan, alluvial fan and terrestrial glacigenic sedimentary rocks. Samples were analysed for U–Pb isotopes and clast samples additionally for Lu–Hf isotopes. The clasts are mainly granites c. 960 Ma and 1680 Ma old, coeval with the Sveconorwegian Orogeny and formation of the Palaeoproterozoic Transscandinavian Igneous Belt (TIB). Mesoproterozoic (Sveconorwegian) ages are abundant in the western part of the basin, whereas Palaeoproterozoic ages are common in the east. Lu–Hf isotopes support crustally contaminated source for all clasts linking them to Fennoscandia. Detrital zircon ages of the sandstones can be matched with those from the granitic clasts except for ages within the range 1200–1500 Ma. These ages are typically found in the present-day Telemark, SW Norway. The sandstones and conglomerate clasts in the western part of the Hedmark Basin were sourced from the Sveconorwegian domain in the present SW Norway or its continuation to the present-day NW. The conglomerate clasts in the eastern part of the Hedmark Basin were sourced mainly from the TIB domain or its northwesterly continuation. The Hedmark Basin was initiated within the boundary of two domains in the basement: the TIB and the Sveconorwegian domains.
The role of travertine fissure-ridges in reconstructing tectonics and related earthquakes is a challenging issue of recent debate directed at delineating historical/prehistorical seismic records. Indeed, direct measurements on a travertine fissure-ridge immediately after a seismic event have never been previously performed. We describe the co- and post-seismic effects of a M = 3.6 earthquake on fluid flow and travertine deposition in a geothermal area of Tuscany (Italy). Direct observation allows us to demonstrate that thermal spring (re)activation is directly influenced by transient seismic waves, therefore providing a basis for reconstructing seismic events in the past.
The radiometric method, or radiometrics, measures naturally occurring radioactivity in the form of gamma-rays (γ-rays) (Fig. 4.1). Most of this radiation originates from mineral species containing radioactive isotopes of potassium (K), uranium (U) and thorium (Th). Radiometrics is a passive geophysical method because it measures a natural source of energy. Radiometric surveys for mineral exploration are routinely made from the air, on the ground and within drillholes. Airborne radiometrics is particularly common in mineral exploration where the radiometric data are acquired simultaneously with magnetics during airborne surveying, although it would be unusual for a solely radiometric survey to be conducted. Ground radiometric surveys are usually conducted with hand-held instruments, but γ-ray detectors are sometimes mounted on a moving vehicle for larger-scale surveys. Downhole radiometric measurements are limited to γ-logging; downhole surveying as defined in Section 1.2.1 is not possible in radiometrics. Gamma-logging is a common component of multiparameter logging for rapidly and economically measuring in situ physical properties, correlating stratigraphy between drillholes and assisting in the evaluation of uranium deposits. It is also used extensively by the petroleum industry in well-logging.
Historically, the main use of radiometrics in mineral prospecting was detection of anomalies caused by outcropping, highly radioactive, uranium deposits. With improvements in sensor technology and data-processing algorithms, downhole radiometric techniques were developed for estimating the grade of uranium ores intersected by drillholes. The advent of multichannel detectors for airborne surveying, capable of distinguishing radiation from different radioactive elements, the increased sensitivity and resolution of airborne surveying techniques, and the development of new data reduction algorithms have focused airborne radiometrics more toward geological mapping than purely anomaly detection. Newer applications include detecting and mapping areas of hydrothermal alteration and weakly radioactive mineral deposits, e.g. heavy-mineral sands. Renewed interest in the method has also been driven by the requirement to map the materials making up the near-surface. In this context, the radiometric data are usually interpreted in combination with topographic and satellite-borne remote-sensing data. Some applications include regolith mapping for mineral exploration and the mapping of soil types for environmental studies.
Peter Moczo, Univerzita Komenského v Bratislave, Slovakia,Jozef Kristek, Univerzita Komenského v Bratislave, Slovakia,Martin Gális, King Abdullah University of Science and Technology, Saudi Arabia
Peter Moczo, Univerzita Komenského v Bratislave, Slovakia,Jozef Kristek, Univerzita Komenského v Bratislave, Slovakia,Martin Gális, King Abdullah University of Science and Technology, Saudi Arabia
In the previous chapter we studied the propagation of rupture in a material that we assumed to be otherwise unbroken, continuous and perfectly elastic outside the fracture. A direct consequence of these assumptions is that the stress and slip rate present inverse square root singularities as we approach the tip of the crack. These two singularities result from approximations that need to be corrected, since no material of finite strength can sustain infinite stresses or slip at an infinite rate. Moreover, homogeneous models cannot explain the nucleation and arrest of rupture nor the super-shear rupture speeds that sometime occur in earthquake ruptures. To avoid these problems we need to introduce heterogeneities into the medium. These will explain the beginning of the fracture process and its arrest and will also eliminate the singularities of stress and slip rate at the rupture edge.
The cohesive zone
We start with a study of the singularities at the edge of a propagating fracture. In order that the values of the stress and slip rate near the rupture front be finite it is necessary to introduce a transition zone where the material behaves inelastically. Barenblatt (1959a, b, 1962) proposed that the stresses immediately ahead of the fracture front should be regularized by the action of cohesive forces, which hold together the material while it is in the process of being broken. These forces are distributed in the vicinity of the rupture front and may be considered as a material property. They are analogous to the attractive or cohesive forces exerted on a liquid molecule by neighboring molecules which explains their name. The transition zone where these forces act is called the cohesive zone and its dimension d must be small compared with the size of the fracture (Fig. 11.1). At the front of the cohesive zone the opposite sides at the edge of the crack are smoothly joined and the stresses remain finite. Let us consider a very simple model of a cohesive zone of length d where the cohesive stresses are constant and equal to σc, which is a compressive stress that opposes the opening of the cohesive zone. Generally, the cohesive stress σc is larger than the acting stress σ0 (Fig. 11.1). Behind the cohesive zone, at distance d from the tip, the tensional stresses drop to zero.
Peter Moczo, Univerzita Komenského v Bratislave, Slovakia,Jozef Kristek, Univerzita Komenského v Bratislave, Slovakia,Martin Gális, King Abdullah University of Science and Technology, Saudi Arabia
The seismic method is an active form of geophysical surveying that uses elastic waves to investigate the subsurface. The waves are created by a source and propagate through the subsurface before being recorded by detectors that measure deformation of the ground (Fig. 6.1). The deformation of the ground as a function of time since the waves were created comprises a time series, which is called a seismic trace. The passing of a seismic wave appears as a deflection of the trace, referred to as an arrival or an event. The path of the waves from source to detector is controlled by the elastic properties of the material through which they travel. Discontinuities in the elastic properties deflect and divide the seismic waves so that the detectors record a series of waves that have taken different paths through the subsurface. It is by identifying these different arrivals and analysing their travel times and amplitudes that the nature of the subsurface can be inferred.
Most seismic surveys use sources and detectors located on, or near, the Earth’s surface. Surface seismic surveys are of two types. Surveys that particularly exploit waves reflected at elastic discontinuities are known as seismic reflection surveys, whilst those based on waves deflected so as to travel parallel to the discontinuities are known as seismic refraction surveys. Of all the geophysical methods, seismic reflection surveys provide the most detailed information about the subsurface, albeit at the greatest cost (see Fig. 1.3). Reflection survey data are displayed in a pseudosection/pseudovolume form where the travel times of the seismic waves provide an approximate indication of depth; but accurate time-to-depth conversion is difficult. Refraction surveys provide less detailed information, detecting only major changes in elastic properties, although the positions of these boundaries can usually be accurately determined.
Peter Moczo, Univerzita Komenského v Bratislave, Slovakia,Jozef Kristek, Univerzita Komenského v Bratislave, Slovakia,Martin Gális, King Abdullah University of Science and Technology, Saudi Arabia
from
Part III
-
Finite-element method and hybrid finite-difference–finite-element method
Peter Moczo, Univerzita Komenského v Bratislave, Slovakia,Jozef Kristek, Univerzita Komenského v Bratislave, Slovakia,Martin Gális, King Abdullah University of Science and Technology, Saudi Arabia
from
Part IV
-
Finite-difference modelling of seismic motion at real sites
Peter Moczo, Univerzita Komenského v Bratislave, Slovakia,Jozef Kristek, Univerzita Komenského v Bratislave, Slovakia,Martin Gális, King Abdullah University of Science and Technology, Saudi Arabia
Peter Moczo, Univerzita Komenského v Bratislave, Slovakia,Jozef Kristek, Univerzita Komenského v Bratislave, Slovakia,Martin Gális, King Abdullah University of Science and Technology, Saudi Arabia
Peter Moczo, Univerzita Komenského v Bratislave, Slovakia,Jozef Kristek, Univerzita Komenského v Bratislave, Slovakia,Martin Gális, King Abdullah University of Science and Technology, Saudi Arabia
Peter Moczo, Univerzita Komenského v Bratislave, Slovakia,Jozef Kristek, Univerzita Komenského v Bratislave, Slovakia,Martin Gális, King Abdullah University of Science and Technology, Saudi Arabia