Seismicity and seismic gap in the Lesser Antilles arc and earthquake hazard in Guadeloupe

Summary. A seismic study of the Lesser Antilles arc has been carried out, first for the period 1950–1978, for which we can use local seismic networks to draw maps of instrumental seismicity, then for the period 1530–1950, for which we have catalogues of felt earthquakes. The striking feature of the spatial distribution of foci is the cluster of epicentres in the northern half of the arc; all large earthquakes ( M > 7.5) are located north of 14° latitude. Seismicity cross-sections through the arc show a variable dipping subduction zone along the arc; the deep seismic zone is steeper in the centre of the arc than on the extremity.
The time-space diagram for historical seismicity, and the evidence of a seismic gap at the east of Guadeloupe lead us to consider the northern half arc as a likely site for a large earthquake in the near future.
The seismic slip rate calculated from all major earthquakes since 1530 is of much greater value than that obtained from recent plate tectonic models, suggesting that the recurrence rate of earthquakes is more than many hundreds of years with a possible aseismic creep.     

Constraints on the frequency–magnitude relation and maximum magnitudes in the UK from observed seismicity and glacio-isostatic recovery rates

Earthquake populations have recently been shown to have many similarities with critical-point phenomena, with fractal scaling of source sizes (energy or seismic moment) corresponding to the observed Gutenberg–Richter (G–R) frequency–magnitude law holding at low magnitudes. At high magnitudes, the form of the distribution depends on the seismic moment release rate M˙ and the maximum magnitude m _max . The G–R law requires a sharp truncation at an absolute maximum magnitude for finite M˙ . In contrast, the gamma distribution has an exponential tail which allows a soft or 'credible' maximum to be determined by negligible contribution to the total seismic moment release. Here we apply both distributions to seismic hazard in the mainland UK and its immediate continental shelf, constrained by a mixture of instrumental, historical and neotectonic data. Tectonic moment release rates for the seismogenic part of the lithosphere are calculated from a flexural-plate model for glacio-isostatic recovery, constrained by vertical deformation rates from tide-gauge and geomorphological data. Earthquake focal mechanisms in the UK show near-vertical strike-slip faulting, with implied directions of maximum compressive stress approximately in the NNW–SSE direction, consistent with the tectonic model. Maximum magnitudes are found to be in the range 6.3–7.5 for the G–R law, or 7.0–8.2 m _L for the gamma distribution, which compare with a maximum observed in the time period of interest of 6.1 m _L . The upper bounds are conservative estimates, based on 100 per cent seismic release of the observed vertical neotectonic deformation. Glacio-isostatic recovery is predominantly an elastic rather than a seismic process, so the true value of m _max is likely to be nearer the lower end of the quoted range.     

Recurrence-time and frequency-slip statistics of slip events on the creeping section of the San Andreas fault in central California

An important question in seismic hazard assessments is the frequency-size and recurrence interval statistics at a point on a fault. Does a point on a fault obey the same statistics as earthquakes in a region do? This is a difficult question to answer because the number of repetitive earthquakes on a particular fault that have been observed is small. In order to overcome this difficulty we consider slip events on the creeping section of the San Andreas fault in central California. Sequences of up to 100 events are obtained from creepmeter records. We compare the statistical distribution of recurrence times with the Brownian passage-time, lognormal, and Weibull distributions and using goodness-of-fit tests find that the Weibull is the preferred distribution. We also consider the frequency–amplitude distribution of slip events. We find that the data clearly do not obey a Gutenberg–Richter distribution. Instead there is a uniform distribution of event sizes for a large fraction of events.     

Earthquake deformations in the Lisan deposits and seismotectonic implications

Summary. Earthquake deformations and induced sedimentary structures preserved in Quaternary sediments include faults, folds, fissures, slumps, sand boils and other effects of liquefaction. Such deformations and structures are well preserved in the Lisan deposits of the Dead Sea. Of most importance are the fold-type deformations known as décollement structures which are present all along the eastern side of the Lisan and seem to decrease gradually westwards to disappear in the middle of the Lisan. These may indicate that palaeoearthquakes originating along the Araba fault have triggered such structures due to shaking of elastoplastic unconsolidated sediments over gentle slopes dipping to the west.
Preliminary results from studies on décollement structures preserved in a section representing some 1733 years of continuous deposition in the uppermost? Pleistocene, in the vicinity of Wadi Araba, indicate that: (1) seismic activity has fluctuated with time. Average recurrence period is about 340 ± 20yr for earthquakes with magnitudes greater than or equal to 6.5, Earthquakes with magnitude greater than 7 seem to have occurred along the Araba fault. (2) Deduced earthquake magnitudes conform to the frequency–magnitude relationship: log N = 5.24–0.68 M . (3) The deduced seismic slip rate along the Araba fault seems to be not less than 0.64 ± 0.04 cm yr⁻¹.     

Coseismic slip model of the 2007 August Pisco earthquake (Peru) as constrained by Wide Swath radar observations

The Pisco earthquake ( M _w 8.0; 2007 August 15) occurred offshore of Peru's southern coast at the subduction interface between the Nazca and South American plates. It ruptured a previously identified seismic gap along the Peruvian margin. We use Wide Swath InSAR observations acquired by the Envisat satellite in descending and ascending orbits to constrain coseismic slip distribution of this subduction earthquake. The data show movement of the coastal regions by as much as 85 cm in the line-of-sight of the satellite. Distributed-slip model indicates that the coseismic slip reaches values of about 5.5 m at a depth of ∼18–20 km. The slip is confined to less than 40 km depth, with most of the moment release located on the shallow parts of the interface above 30 km depth. The region with maximum coseismic slip in the InSAR model is located offshore, close to the seismic moment centroid location. The geodetic estimate of seismic moment is 1.23 × 10²¹ Nm ( M _w 8.06), consistent with seismic estimates. The slip model inferred from the InSAR observations suggests that the Pisco earthquake ruptured only a portion of the seismic gap zone in Peru between 13.5° S and 14.5° S, hence there is still a significant seismic gap to the south of the 2007 event that has not experienced a large earthquake since at least 1687.     

Discrimination Analysis of Earthquakes and Man-Made Events Using ARMA Coefficients Determination by Artificial Neural Networks

A Quadratic Neural Networks (QNNs) model has been developed for identifying seismic source classification problem at regional distances using ARMA coefficients determination by Artificial Neural Networks (ANNs). We have devised a supervised neural system to discriminate between earthquakes and chemical explosions with filter coefficients obtained by windowed P-wave phase spectra (15 s). First, we preprocess the recording's signals to cancel out instrumental and attenuation site effects and obtain a compact representation of seismic records. Second, we use a QNNs system to obtain ARMA coefficients for feature extraction in the discrimination problem. The derived coefficients are then applied to the neural system to train and classification. In this study, we explore the possibility of using single station three-component (3C) covariance matrix traces from a priori-known explosion sites (learning) for automatically recognizing subsequent explosions from the same site. The results have shown that this feature extraction gives the best classifier for seismic signals and performs significantly better than other classification methods. The events have been tested, which include 36 chemical explosions at the Semipalatinsk test site in Kazakhstan and 61 earthquakes (mb = 5.0–6.5) recorded by the Iranian National Seismic Network (INSN). The 100% correct decisions were obtained between site explosions and some of non-site events. The above approach to event discrimination is very flexible as we can combine several 3C stations.     

Traveltime measurements from noise correlation: stability and detection of instrumental time-shifts

We test the feasibility of using Green's functions extracted from records of ambient seismic noise to monitor temporal changes in the Earth crust properties by repeated measurements at regional distances. We use about 11 yr of continuous recordings to extract surface waves between three pairs of stations in California. The correlations are computed in a moving 1-month window and we analyse the temporal evolution of measured interstation traveltimes. The comparison of the arrival times in the positive and negative correlation time of Rayleigh and Love waves allows us to separate time-shifts associated with any form of physical change in the medium, those resulting from clock drift or other instrumental errors, and those due to change in the localization of the noise sources. This separation is based on the principle of time symmetry. When possible, we perform our analysis in two different period bands: 5–10 and 10–20 s. The results indicate that significant instrumental time errors (0.5 s) are present in the data. These time-shifts can be measured and tested by closure relation and finally corrected independently of any velocity model. The traveltime series show a periodic oscillation that we interpret as the signature of the seasonal variation of the region of origin of the seismic noise. Between 1999 and 2005, the final arrival time fluctuations have a variance of the order of 0.01 s. This allows us to measure interstation traveltimes with errors smaller than 0.3 per cent of the interstation traveltime and smaller than 1 per cent of the used wave period. This level of accuracy was not sufficient to detect clear physical variation of crustal velocity during the considered 11 yr between the three stations in California. Such changes may be more easily detectable when considering pairs of stations more closely located to each other and in the vicinity of tectonically active faults or volcanoes.     

Quaternary normal faulting in southeastern Sicily (Italy):a seismic source for the 1693 large earthquake

We present geological and morphological data, combined with an analysis of seismic reflection lines across the Ionian offshore zone and information on historical earthquakes, in order to yield new constraints on active faulting in southeastern Sicily. This region, one of the most seismically active of the Mediterranean, is affected by WNW–ESE regional extension producing normal faulting of the southern edge of the Siculo–Calabrian rift zone. Our data describe two systems of Quaternary normal faults, characterized by different ages and related to distinct tectonic processes. The older NW–SE-trending normal fault segments developed up to ≈400  kyr ago and, striking perpendicular to the main front of the Maghrebian thrust belt, bound the small basins occurring along the eastern coast of the Hyblean Plateau. The younger fault system is represented by prominent NNW–SSE-trending normal fault segments and extends along the Ionian offshore zone following the NE–SW-trending Avola and Rosolini–Ispica normal faults. These faults are characterized by vertical slip rates of 0.7–3.3  mm  yr ⁻¹ and might be associated with the large seismic events of January 1693. We suggest that the main shock of the January 1693 earthquakes ( M ~ 7) could be related to a 45  km long normal fault with a right-lateral component of motion. A long-term net slip rate of about 3.7  mm  yr ⁻¹ is calculated, and a recurrence interval of about 550 ± 50  yr is proposed for large events similar to that of January 1693.     

A multidisciplinary study of a slow-slipping fault for seismic hazard assessment: the example of the Middle Durance Fault (SE France)

Assessing seismic hazard in continental interiors is difficult because these regions are characterized by low strain rates and may be struck by infrequent destructive earthquakes. In this paper, we provide an example showing that interpretations of seismic cross sections combined with other kinds of studies such as analysis of microseismicity allow the whole seismogenic source area to be imaged in this type of region. The Middle Durance Fault (MDF) is an 80-km-long fault system located southeastern France that has a moderate but regular seismicity and some palaeoseismic evidence for larger events. It behaves as an oblique ramp with a left-lateral-reverse fault slip and has a low strain rate. MDF is one of the rare slow active fault system monitored by a dedicated dense velocimetric short period network. This study showed a fault system segmented in map and cross section views which consists of staircase basement faults topped by listric faults ramping off Triassic evaporitic beds. Seismic sections allowed the construction of a 3-D structural model used for accurate location of microseismicity. Southern part of MDF is mainly active in the sedimentary cover. In its northern part and in Alpine foreland, seismicity deeper than 8 km was also recorded meaning active faults within the crust cannot be excluded. Seismogenic potential of MDF was roughly assessed. Resulting source sizes and estimated slip rates imply that the magnitude upper limit ranges from 6.0 to 6.5 with a return period of a few thousand years. The present study shows that the coupling between 3-D fault geometry imaging and accurate location of microseismicity provides a robust approach to analyse active fault sources and consequently a more refined seismic hazard assessment.     

Source investigation of a small event using empirical Green's functions and simulated annealing

We propose a two-step inversion of three-component seismograms that (1) recovers the far-field source time function at each station and (2) estimates the distribution of co-seismic slip on the fault plane for small earthquakes (magnitude 3 to 4). The empirical Green's function (EGF) method consists of finding a small earthquake located near the one we wish to study and then performing a deconvolution to remove the path, site, and instrumental effects from the main-event signal.
The deconvolution between the two earthquakes is an unstable procedure: we have therefore developed a simulated annealing technique to recover a stable and positive source time function (STF) in the time domain at each station with an estimation of uncertainties. Given a good azimuthal coverage, we can obtain information on the directivity effect as well as on the rupture process. We propose an inversion method by simulated annealing using the STF to recover the distribution of slip on the fault plane with a constant rupture-velocity model. This method permits estimation of physical quantities on the fault plane, as well as possible identification of the real fault plane.
We apply this two-step procedure for an event of magnitude 3 recorded in the Gulf of Corinth in August 1991. A nearby event of magnitude 2 provides us with empirical Green's functions for each station. We estimate an active fault area of 0.02 to 0.15 km² and deduce a stress-drop value of 1 to 30 bar and an average slip of 0.1 to 1.6 cm. The selected fault of the main event is in good agreement with the existence of a detachment surface inferred from the tectonics of this half-graben.     

Late Quaternary systematic stream offsets caused by repeated large seismic events along the Kunlun fault, northern Tibet

The Kunlun fault is one of the largest strike-slip faults in northern Tibet, China. In this paper, we focus upon the Kusai Lake–Kunlun Pass segment of the fault to understand the geomorphic development of offset streams caused by repeated large seismic events, based on tectono-geomorphic analysis of high-resolution satellite remote sensing images combined with field studies. The results indicate that systematic left-lateral stream offsets appear at various scales across the fault zone: Lateral offsets of small gullies caused by the 2001 M_w 7.8 Kunlun earthquake vary typically from 3 m to 6 m, meanwhile streams with cumulative offsets of 10 m, 25–30 m, 50–70 m, 250–300 m and 750–1400 m have resulted from repeated large seismic events during the late Quaternary. An average slip rate of 10 ± 1 mm/year has been estimated from the lateral stream offsets and ¹⁴C ages of alluvial fan surfaces incised by the streams. A three-dimensional model showing tectono-geomorphic features along a left-lateral strike-slip fault is also presented. The Kusai Lake–Kunlun Pass segment provides an opportunity to understand the relationship between geomorphic features produced by individual large seismic events and long-term geomorphic development caused by repeated large seismic events along a major strike-slip fault.     

The role of seismicity models in probabilistic seismic hazard estimation: comparison of a zoning and a smoothing approach

Seismic hazard estimations are compared using two approaches based on two different seismicity models: one which models earthquake recurrence by applying the truncated Gutenberg-Richter law and a second one which smoothes the epicentre location of past events according to the fractal distribution of earthquakes in space ( Woo 1996 ). The first method requires the definition of homogeneous source zones and the determination of maximum possible magnitudes whereas the second method requires the definition of a smoothing function. Our results show that the two approaches lead to similar hazard estimates in low seismicity regions. In regions of increased seismic activity, on the other hand, the smoothing approach yields systematically lower estimates than the zoning method. This epicentre-smoothing approach can thus be considered as a lower bound estimator for seismic hazard and can help in decision making in moderate seismicity regions where source zone definition and estimation of maximum possible magnitudes can lead to a wide variety of estimates due to lack of knowledge. The two approaches lead, however, to very different earthquake scenarios. Disaggregation studies at a representative number of sites show that if the distributions of contributions according to source–site distance are comparable between the two approaches, the distributions of contributions according to magnitude differ, reflecting the very different seismicity models used. The epicentre-smoothing method leads to scenarios with predominantly intermediate magnitudes events (5 ≤ M ≤ 5.5) while the zoning method leads to scenarios with magnitudes that increase with the return period from the minimum to the maximum magnitudes considered. These trends demonstrate that the seismicity model used plays a fundamental role in the determination of the controlling scenarios and ways to discriminate between the most appropriate models remains an important issue.     

Comparison of frequency of occurrence of earthquakes with slip rates from long-term seismicity data: the cases of Gulf of Corinth, Sea of Marmara and Dead Sea Fault Zone

In this article, through the comparison of knowledge relating to historical earthquakes with the understanding of present-day earthquake mechanics and overall GPS slip rates in the eastern Mediterranean region, it has been possible to obtain an idea of how frequently large earthquakes may be expected in some parts of the region. It has also been possible to make an assessment from these early events of slip rates over a long period of time for the Gulf of Corinth in Greece, the Marmara Sea in Turkey and the Dead Sea Fault System, as well as deriving long-term magnitude–frequency relations for these same regions.
It has been demonstrated that slip rates calculated from historical data are in general comparable to those calculated from GPS measurements and field observations, while the size of historical earthquakes and their uncertainty can be quantified. This permits a more reliable estimation of the long-term hazard, the calculation of which is the concern of the engineering seismologist. It has also been shown that in most cases large earthquakes are less frequent when they are estimated from long-term data sets rather than from the instrumental period making the notion of recurrence time and of hazard assessment, questionable.
This study focuses on some of the few areas in the world for which long-term macroseism information exists and which facilitate this kind of analysis.     

Spherical harmonic representation of the magnetic field in the presence of a current density

Summary. We present the results of a systematic study of events with M _s > 6 in northern Chile (20–33°S), for the period between 1963 and 1971. Medium to large earthquakes near the coast of this region are of three types: (1) Interplate events at the interface between the downgoing slab and the overriding South American plate. These events can be very large reaching magnitudes greater than 8. (2) Intra-plate earthquakes 20–30 km inside the downgoing slab. They have fault mechanisms indicating extension along the dip of the slab and may have magnitudes up to 7.5. (3) Less frequent, M _s∼ 6 events that occur near the top of the downgoing slab and have thrust mechanisms with an almost horizontal E-W compressional axis. This type of mechanism is very different from that of the events of type 1 which are due to shallow dipping reverse faulting. There is a rotation of about 30° of the compressional axis in the vertical plane between events of types (1) and (3). Three groups of events near 32.5°, 25.5° and 21°s were studied in detail. Depth and mechanisms were redetermined by P -wave modelling and relative locations were obtained by a master event technique. Near 32.5°S, only events of types 1 and 2 were found in the time period of this study. At the two other sites, the three types of events were identified. This shows clearly that there are compressive stresses at the top of the slab and extension at the centre, a situation which is usually found in the areas where a double Benioff-zone has been identified in the seismicity.     

基于信息扩散理论的中国西部地区地震风险评估

利用全国地震目录数据库,系统梳理建国以来中国西部地区5.0级以上地震信息;利用信息扩散理论,从年最大震级和各级地震灾害年频次两个方面对西部12省、市、自治区的地震灾害进行风险评估。在此基础上,将风险估计离散值进行曲线拟合,系统评价中国西部地震最为多发5省区的地震风险分布趋势及差异,定量揭示中国西部地区地震风险,为西部地区地震防灾减灾提供技术支持和决策建议。     

A catalogue of seismicity in Greece and adjacent areas

Summary. A new earthquake catalogue for Greece has been formed to cover the instrumental period 1901–78, in particular 605 earthquakes for the period 1917–63 inclusive are relocated using first arrival data from the International Seismological Summary. These relocations incorporate macro-seismically and other well-controlled master events into an ensuing joint epicentre determination technique. The largest annual average shift is 165 km for 78 earthquakes during the decade after 1917, decreasing to 17 km for the later years 1957–63. Magnitudes are redetermined mainly using readings from the Swedish network and Uppsala since as early as 1908. Catalogue completeness exists for magnitudes around 5.5 for at least the last 60 years, but magnitudes of about 4.7 are completely reported only during the most recent 15 years.
The new epicentral positions lead to a better delineation of the seismic zones than has previously been achieved. The majority of shallow earthquakes are contained in a belt parallel to the Hellenic Arc which extends north into Albania, in the south-east the epicentral locations give a more diffuse extension of this zone into the west coast of Turkey. Depths of intermediate earthquakes along the Hellenic Arc tend to relocate to shallower depths, in the south-west part of Crete no earthquake focus deeper than 100 km is found, although in the south-eastern section of the arc a tendency to increased depths is observed. A second zone starts at Leukas Island in the west and extends through central Greece to near Volos on the east coast, where it divides into two branches, which are less well defined, but which eventually join the seismicity of western Turkey. A third zone follows the Saronikos and Corinth gulfs.     

A major seismogenic fault in a 'silent area': the Castrovillari fault (southern Apennines, Italy)

Large historical earthquakes in Italy define a prominent gap in the Pollino region of the southern Apennines. Geomorphic and palaeoseismological investigations in this region show that the Castrovillari fault (CF) is a major seismogenic source that could potentially fill the southern part of this gap. The surface expression of the CF is a complex, 10–13 km long set of prominent scarps. Trenches across one scarp indicate that at least four surface-faulting earthquakes have occurred along the CF since Late Pleistocene time, each producing at least 1 m of vertical displacement. The length of the fault and the slip per event suggest M =6.5-7.0 for the palaeoearthquakes. Preliminary radiocarbon dating coupled with historical considerations imply that the most recent of these earthquakes occurred between 380 BC and 1200 AD, and probably soon after 760 AD; no evidence for this event has been found in the historical record. We estimate a minimum recurrence interval of 1170 years and a vertical slip rate of 0.2-0.5 mm yr^-1 for the CF, which indicates that the seismic behaviour of this fault is comparable to other major seismogenic faults of the central-southern Apennines. The lack of mention or the mislocation of the most recent event in the historical seismic memory of the Pollino region clearly shows that even in Italy, which has one of the longest historical records of seismicity, a seismic hazard assessment based solely on the historical record may not be completely reliable, and shows that geological investigations are critical for filling possible information gaps.     

FIR filter effects and nucleation phases

The symmetric impulse response of linear phase Finite Impulse Response (FIR) filters most commonly used in modern seismic recording systems produces precursory signals to impulsive arrivals. These acausal filter-generated artefacts may result in misinterpretations of various onset properties. Prior to any onset interpretation, these effects have to be removed from the seismic record. This can be achieved without loss of bandwidth by post-filtration of the digital seismograms if the filter coefficients and the decimation ratios are known. We have analysed numerous signals from different instruments and sampling rates for precursory phases and found that—in contrast to commonly held beliefs—FIR-filter-related precursory signals are not always easy to recognize visually from their waveform signature. Furthermore, they can exhibit surprisingly similar properties to those reported for nucleation phases, although the majority of nucleation phases reported in the past have been obtained on instruments with a causal response. We demonstrate examples of filter-related precursory signals for events scanning nine orders of moment, from 10¹⁰ N m to 10¹⁹ N m. Surprisingly, the lower bound of the artefact durations as a function of seismic moment scales close to the cube root of the seismic moment. We interpret this as being caused by the fact that above a certain seismic moment, the attenuated source signal acts as a causal lowpass filter of a smaller bandwidth than the FIR filter. Assuming an ω^-2 source model, constant stress drop and an empirical relationship between the maximum artefact duration and the cut-off frequency of the FIR filter, the artefact durations are expected to scale proportional to the 1/2.5 power of the seismic moment, in comparison to 1/3 as proposed for nucleation phases.     

Mitigating Risks From Fracking-Related Earthquakes: Assessing State Regulatory Decisions

Public concern about earthquakes linked to wastewater injection from fracking operations is rising. However, few have examined how “induced seismicity” is acted upon by state officials. For some, an incremental response to smaller quakes can be viewed as an acceptable risk policy orientation because of the sizeable economic benefits that accompany drilling activities while others prefer risk mitigation policies (such as the use of “threshold policies”) as a better way to address quake-related problems. To account for state response to induced seismicity impacts, we examine three factors: the emergence of quakes as focusing events, the economic importance of oil and gas to state jobs and revenue, and selected characteristics of earthquakes as a policy issue, i.e., complexity and categorical precedence. Using information drawn from documentary sources, we consider which factors are most helpful in accounting for agency decisions aimed at reducing seismic risks linked to nearby injection wells.