First seismic record for intra-arc strike-slip tectonics along the Liquiñe-Ofqui fault zone at the obliquely convergent plate margin of the southern Andes

A temporal seismic network recorded local seismicity along a 130 km long segment of the transpressional dextral strike-slip Liquiñe-Ofqui fault zone (LOFZ) in southern Chile. Seventy five shallow crustal events with magnitudes up to M_w 3.8 and depths shallower than 25 km were observed in an 11-month period mainly occurring in different clusters. Those clusters are spatially related to the LOFZ, to the volcanoes Chaitén, Michinmahuida and Corcovado, and to active faulting on secondary faults. Further activity along the LOFZ is indicated by individual events located in direct vicinity of the surface expression of the LOFZ. Focal mechanisms were calculated using deviatoric moment tensor inversion of body wave amplitude spectra which mostly yield strike-slip mechanisms indicating a NE–SW direction of the P-axis for the LOFZ at this latitude. The seismic activity reveals the present-day activity of the fault zone. The recent M_w 6.2 event near Puerto Aysén, Southern Chile at 45.4°S on April 21, 2007 shows that the LOFZ is also capable of producing large magnitude earthquakes and therefore imposing significant seismic hazard to this region.     

Focal mechanism inversion in the Giudicarie–Lessini seismotectonic region (Southern Alps, Italy): Insights on tectonic stress and strain

A set of 41 focal mechanisms (1989–2006) from P-wave first polarities is computed from relocated seismic events in the Giudicarie–Lessini region (Southern Alps). Estimated hypocentral depths vary from 3.1 to 20.8 km, for duration magnitudes (M_D) in the range 2.7–5.1. Stress and strain inversions are performed for two seismotectonic zones, namely G (Giudicarie) and L (Lessini). This subdivision is supported by geological evidence, seismicity distribution, and focal mechanism types. The available number of data (16 in G, 22 in L) does not make possible any further subdivisions. Seismotectonic zones G and L are undergoing different kinematic regimes: thrust with strike-slip component in G, and strike-slip in L. Principal stress and strain axes in each sub-region show similar orientations. The direction of maximum horizontal compressive stress is roughly perpendicular to the thrust fronts along the Giudicarie Belt in zone G, and compatible with right-lateral strike-slip reactivation of the faults belonging to the Schio-Vicenza system in zone L. On the whole, kinematic regimes and horizontal stress orientations show a good fit with other stress data from focal mechanisms and breakouts and with geodetic strain rate axes.     

Implications of north pacific plate tectonics in central Alaska: Focal mechanisms of earthquakes

The focal mechanisms for 86 selected earthquakes (3.0 m_b 5.5) located in central Alaska have been investigated from P-wave first motions; the data were gathered by local seismic networks. The results show a depth-dependent characteristic to the fault-plane solutions. For earthquakes having focal depths shallower than 60–70 km, the focal mechanisms indicate either strike-slip or normal faults, while for earthquakes with foci at intermediate depths the focal mechanisms correspond to thrust faults. The nature of the seismicity indicates the hinge line of the Pacific lithospheric plate under the study area to be striking N17°E from Cook Inlet towards interior Alaska. The comparison of the focal mechanisms with the seismicity shows that the strike-slip and normal faults are the predominant processes of stress release along the shallow section of the plate. The earthquakes with intermediate foci systematically occur along the inclined section of the plate. If the gently dipping nodal planes for these earthquakes are chosen as the fault planes, the focal mechanisms correspond to underthrust motions at the foci. In these, the slip vectors are oriented either to the west or north with the resultant being in the N30°W direction. The tension axes for the underthrust solutions are also found to be parallel to the local dip of the plate, indicating that the subducted plate in interior Alaska is undergoing gravitational sinking.     

Spatial variations in the similarity of earthquake populations: The case of the 1997 Colfiorito–Sellano (northern Apennines, Italy) seismic sequence

The spatial properties of events in the 1997 Colfiorito–Sellano seismic sequence (Northern Apennines, Italy) were investigated using coherence, a parameter derived from seismic moment tensors that quantifies the kinematic similarity between focal mechanisms. The 1997 Colfiorito–Sellano seismic sequence predominantly consists of normal faulting earthquakes, with a few strike-slip and reverse faulting episodes. This kinematic heterogeneity is possibly related to the contemporaneous activity of two different sets of faults: NW–SE normal faults and NNE–SSW sub-vertical faults, the latter inherited from the previous Miocene compressional phase. The study used two independently-derived data sets of the same seismic sequence characterized by a different number of events and by different precision of spatial localisation. Their statistical significances, assessed through a reshuffling procedure, reveal that data sets with at least some hundreds of events and good positional precision are required to obtain significant results through coherence analysis. Results from the better quality data set indicate that this seismic sequence is characterized by a rapid decrease in the kinematic similarity between earthquake pairs within 2 km of separation, particularly along directions sub-perpendicular to the normal fault strike. The decrease rate seems to be controlled by the geometric characteristics of the normal faults, given that the mean along-dip distance between fault segments is 2 km. In proximity to pre-existing tectonic lineaments the relative abundance of strike-slip and reverse faults tends to decrease the kinematic similarity between events but does not influence the coherence decrease rate. The presence of mixed focal mechanisms (normal, reverse and strike-slip) in a single seismic phase implies that mixed fault types are not restricted to polyphase tectonic histories: such heterogeneous kinematics during a single phase may be induced by the presence of inherited discontinuities.     

Hierarchy of earthquakes in the Baikal rift system: implications for lithospheric stress

We investigated space-time lithospheric stress patterns of the Baikal rift system according to the hierarchy of earthquakes using mechanisms of 265 K_P ≥ 10 events recorded from 1950 to 1998 and seismic moments of 802 K_P ≥ 11 events from 1968 to 1994. The lithosphere of the region was confirmed to undergo rifting with mostly normal-slip events, while local areas of frequent strike-slip and reverse motions may record stress heterogeneity. The dominance of rifting, although being evident in the stress dynamics, is unstable, which is indicated by increase in strike-slip and reverse motions to as many as normal slip events in the latest 1980s–earliest 1990s. The lithospheric stress patterns inferred from seismic-moment data are generally consistent with those derived from the classical focal mechanism method. The suggested approach of seismic zoning according to earthquake slip geometry may provide a more reliable background for successful mitigation of seismic hazard in the region.     

A study of microearthquake seismicity and focal mechanisms within the Sea of Marmara (NW Turkey) using ocean bottom seismometers (OBSs)

We have carried out seismological observations within the Sea of Marmara (NW Turkey) in order to investigate the seismicity induced after Gölcük–İzmit (Kocaeli) earthquake (M_w 7.4) of August 17, 1999, using ocean bottom seismometers (OBSs). High-resolution hypocenters and focal mechanisms of microearthquakes have been investigated during this Marmara Sea OBS project involving deployment of 10 OBSs within the Çınarcık (eastern Marmara Sea) and Central-Tekirdağ (western Marmara Sea) basins during April–July 2000. Little was known about microearthquake activity and their source mechanisms in the Marmara Sea. We have detected numerous microearthquakes within the main basins of the Sea of Marmara along the imaged strands of the North Anatolian Fault (NAF). We obtained more than 350 well-constrained hypocenters and nine composite focal mechanisms during 70 days of observation. Microseismicity mainly occurred along the Main Marmara Fault (MMF) in the Marmara Sea. There are a few events along the Southern Shelf. Seismic activity along the Main Marmara Fault is quite high, and focal depth distribution was shallower than 20 km along the western part of this fault, and shallower than 15 km along its eastern part. From high-resolution relative relocation studies of some of the microearthquake clusters, we suggest that the western Main Marmara Fault is subvertical and the eastern Main Marmara Fault dips to south at 45°. Composite focal mechanisms show a strike-slip regime on the western Main Marmara Fault and complex faulting (strike-slip and normal faulting) on the eastern Main Marmara Fault.     

Seismological study of the central Ecuadorian margin: Evidence of upper plate deformation

A seismic study of a segment of the convergent margin of Ecuador is presented. During the SISTEUR campaign a network of 24 Ocean Bottom Seismometers (OBS) was deployed on the Carnegie Ridge, one line along the main axes of the ridge and two lines across the strike of the edge of the ridge, during one month. This marine network was complemented with a land network of 20 stations distributed in two lines: one parallel to the margin and the other perpendicular to it.The seismic event recorded by these networks, were located using different crustal models defined from the wide-angle seismic data modeling. Relative location techniques were used to improve earthquake locations. Seismogram waveform modeling allowed us to constrain hypocentral location for events farther than ～50 km from the network. This modeling also provided additional information to constrain the focal mechanisms of these events. The upper limit of the Interplate Seismogenic Zone (ISZ) is estimated to be at a 10 km depth in the region. The background seismic activity of the upper plate provided new insights:1) A seismic cluster that reaches the base of the overriding plate is linked to the Jipijapa-Portoviejo fault. The reactivation of this Quaternary fault is confirmed by focal mechanisms that provide rupture planes parallel to its superficial projection (N10°–N25°).2) The focal mechanisms presented in this study are compatible with a homogeneous regional stress field corresponding to an E–W to ESE–WNW compression and an NNE–SSW extension. The presence of strike-slip deformation, with a reverse component, corresponds to the NNE escape of the North Andean Block. Normal faulting accommodating this movement suggests that this part of the North Andean Block cannot be considered as a rigid block.     

Non-double-couple mechanisms of microearthquakes induced during the 2000 injection experiment at the KTB site, Germany: A result of tensile faulting or anisotropy of a rock?

Moment tensors of microearthquakes induced during the 2000 injection experiment at the KTB deep drilling borehole at a depth level of 5.4 km are studied. A family of 37 most reliable moment tensors contains significant non-double-couple (non-DC) components. The DC is on average 60% and the non-DC is 40%. Fault plane solutions computed from the DC part show preferred strike-slip mechanisms with small normal or reverse components. A predominant azimuth of P and T axes is in the range of N320°–340°E and of N230°–250°E, respectively. The non-DC components contain both the isotropic (ISO) and compensated linear vector dipole (CLVD) components. The mean value of ISO is 1.5%, the mean value of CLVD is − 5.7%. The predominantly negative CLVD components are inconsistent with the concept of the non-DC mechanisms as a result of tensile faulting due to fluid injection into the rock. The main origin of the non-DC components is probably anisotropy in the focal area. The other origins are errors produced by mismodelling of the medium when calculating the Green functions, and numerical errors produced by noise and limitations of input data. Adopting four alternative models of anisotropy obtained by other seismic measurements at the KTB, we have employed the non-DC components for estimating an optimum orientation of anisotropy in the focal area. The optimum orientation of the symmetry plane of anisotropy is nearly vertical with a strike of N335°–340°E. This strike coincides well with the strike of 330° typical for many major lithological units and faults and with the orientation of the transversely isotropic model inferred by other authors. After removing the anisotropy effects from the moment tensors by calculating the source tensors, the distribution of ISO is significantly narrowed. This indicates predominantly shear, but not tensile faulting.     

The Woods Reef (New South Wales) earthquake of 14 November 1990: Focal mechanism derived from amplitude ratios and synthetic seismograms

On the morning of 15 November 1990 local time, Armidale and the area to the west of Armidale was shaken by a magnitude 3.2 earthquake. The epicentre was located at 30.39° S, 150.88° E and the depth of focus at 12 ± 7 km. As the epicentre was close to the Peel Fault an attempt was made to constrain the focal mechanism of this earthquake. The conventional method, which is based on the analysis of P wave polarities, was not applicable because the event was not strong enough. In an alternative method, the amplitudes of various seismic phases recorded at a number of stations well distributed in azimuth were compared with theoretical amplitudes calculated with the reflectivity method for a point shear dislocation in a layered medium. The differences between observed and calculated amplitudes were minimized as a function of fault strike, fault dip and direction of the slip vector. The analysis indicates that none of the possible fault planes had the strike of the Peel Fault. The solution suggests predominantly strike slip motion along two possible, steeply dipping fault planes. The inferred direction of the maximum compressional stress. is east‐west which is in good agreement with other estimates of the stress field for eastern Australia.     

Fault plane solutions of the january 26th, 2001 bhuj earthquake sequence

A 12-station temporary microearthquake network was established by the Geological Survey of India for aftershock monitoring of the January 26th, 2001 Bhuj earthquake (M _w 7.6) in the Kutch district of Gujarat state, western India. The epicentres of the aftershocks show two major trends: one in the NE direction and the other in the NW direction. Fault-plane solutions of the best-located and selected cluster of events that occurred along the NE trend, at a depth of 15–38 km, show reverse faulting with a large left-lateral strike-slip motion, which are comparable with the main-shock solution. The NW trending upper crustal aftershocks at depth <10 km, on the other hand, show reverse faulting with right-lateral strike-slip motion, and the mid crustal and lower crustal aftershocks, at a depth of 15–38 km, show pure reverse faulting as well as reverse faulting with right-lateral and left-lateral strike-slip motions; these solutions are not comparable with the main-shock solution. It is inferred that the intersection of two faults has been the source area for stress concentration to generate the main shock and the aftershocks.     

The 2012–2013 earthquake swarm in the eastern Guadalquivir basin (South Spain): A case of heterogeneous faulting due to oroclinal bending

From October 2012 to October 2013, a seismic swarm released more than 7000 microearthquakes beneath the eastern Guadalquivir foreland basin. From double-difference relocations of 501 events (m_d > 1.5), we can image the active structures associated with this swarm. Most of the events occurred along two ~ N–S trending lineaments separated ~ 1 km. Relocation places most events at 4–6.5 km depth in the Iberian-massif basement below the basin. Moment tensor inversion yields strike-slip mechanisms consistent with the hypocenter alignments, attributing left-lateral motion to the N–S structures and right-lateral motion to the ESE–WNW ones, in compliance with the ~ NNW direction of the main compressive stress field in the central Betics. These structures respond to a vertical-axis bend in the mountain front associated with the protrusion of Sierra Cazorla east of the epicentral area. This bend is mimicked by concordant, gentle bends in the foreland units, which are evident from the surface geology as well as through structural elements like strike-slip faults, crisscrossing joints. In this context, the right-lateral shear zone responsible for the Torreperogil sequence is taking up deformation in the western limb of the foreland bend.     

The 10 June 2012 Fethiye Mw 6.0 aftershock sequence and its relation to the 24–25 April 1957 Ms 6.9–7.1 earthquakes in SW Anatolia,Turkey

The 10 June 2012 M_w 6.0 aftershock sequence in southwestern Anatolia is examined. Centroid moment tensors for 23 earthquakes with moment magnitudes (M_w) between 3.7 and 6.0 are determined by applying a waveform inversion method. The mainshock is a shallow focus strike-slip with reverse component event at a depth of 30 km. The seismic moment (M_o) of the mainshock is estimated as 1.28 × 10¹⁸ Nm and rupture duration of the Fethiye mainshock is 38 s. The focal mechanisms of the aftershocks are mainly strike-slip faulting with a reverse component. The geometry of the focal mechanisms reveals a strike-slip faulting regime with NE–SW trending direction of T-axis in the entire activated region. A stress tensor inversion of focal mechanism data is performed to obtain a more accurate picture of the Fethiye earthquake stress field. The stress tensor inversion results indicate a predominant strike-slip stress regime with a NW–SE oriented maximum horizontal compressive stress (S_H). According to variance of the stress tensor inversion, to first order, the Fethiye earthquake area is characterized by a homogeneous interplate stress field. The Coulomb stress change associated with the mainshock and the largest aftershock are also investigated to evaluate any significant enhancement of stresses along the Gulf of Fethiye and surrounding region. Positive lobes with stress more than 0.4 bars are obtained, indicating that these values are large enough to increase the Coulomb stress failure towards NNW–SSE and E–W directions.     

Upper mantle anisotropy inferred from shear wave splitting beneath the Eastern Indian Shield region

We estimate the shear wave splitting parameters vis-à-vis the thicknesses of the continental lithosphere beneath the two permanent seismic broadband stations located at Dhanbad (DHN) and Bokaro (BOKR) in the Eastern Indian Shield region. Broadband seismic data of 146 and 131 teleseismic earthquake events recorded at DHN and BOKR stations during 2007–2014 were analyzed for the present measurements. The study is carried out using rotation-correlation and transverse component minimization methods. We retain our “Good”, “Fair” and “Null” measurements, and estimate the splitting parameters using 13 “Good” results for DHN and 10 “Good” results for BOKR stations. The average splitting parameters (ϕ, δt) for DHN and BOKR stations are found to be 50.76°±5.46° and 0.82 ± 0.2 s and 56.30°±5.07° and 0.95 ± 0.17 s, and the estimated average thicknesses of the anisotropic layers beneath these two stations are ∼ 94 and ∼109 km, respectively. The measured deviation of azimuth of the fast axis direction (ϕ) from the absolute motion of the Indian plate ranges from ∼8° to 14°. The measured deviation of azimuth of the fast axis direction (ϕ) from the absolute motion of the Indian plate ranges from ∼8° to 14°. The eastward deviation of the fast axis azimuths from absolute plate motion direction is interpreted to be caused by induced outflow from the asthenosphere. Further, the delay time found in the present analysis is close to the global average for continental shield areas, and also coherent with other studies for Indian shield regions. The five “Null” results and the lower delay time of ∼0.5–0.6 s might be indicating multilayer anisotropy existing in the mantle lithosphere beneath the study area.     

The earthquake swarm of December 2007 in the Mila region of northeastern Algeria

In December 2007, the Mila region of northeastern Algeria experienced thousands of microearthquakes (0.8????Md????3.9) recorded by eight temporary stations, in addition to permanent stations. Most of the events were too small to be located precisely, but a set of 122 precisely located events shows an alignment of epicenters, extending mainly in a horizontal band at about 1?C2?km depth in a NNW?CSSE direction and concentrated in a small area, 3?km southeast of Jebel Akhal, a small rocky hill between the Beni Haroun dam/reservoir and the Oued Athmania reservoir. The reservoirs are connected by pipelines, and a pumping station ensures water can be transferred between them at transient pressures of up to 80?bars. During the pumping in 2007, only 45?% of the transferred water (~600,000?m³ per day) was recovered at the Oued Athmania reservoir, and a large amount of the slightly pressurized water leaked through defective joints in a tunnel that passes through the Jebel. This water penetrated deeply into the soil with the assistance of preexisting fractures, faults, and karsts. Nine days after the first pumping started, a local increase in pore fluid pressures at shallow depths triggered seismicity southeast of Jebel Akhal, where the faults were probably close to failure. The focal mechanisms show a near vertical N?CS strike-slip fault plane under regional NW?CSE tectonic compression. One of the fault plane solutions is consistent with the NNW?CSSE direction along which the seismic events are aligned. Furthermore, a long-term comparison of the seismic activity in the region versus water levels behind the dam and the pumping of water shows that the earthquake swarm was a one-off event related to the pumping operation.     

Waveform inversion and focal mechanisms of two weak earthquakes in Cordillera Principal (Argentina) between 35° and 35.5° S

Only few (six) focal mechanism, in CMT Catalog, have been so far known for intraplate shallow events in the Andean chain close to Chile–Argentina state border at latitudes ∼35° S. We add two more mechanisms, depths and moment magnitudes by carefully analyzing full waveforms of weak events recorded by broad-band stations of the Chile Argentina Geophysical Experiment (southern profile). The moment magnitudes of both events (Mw = 3.6 and 3.7) are lower than the duration magnitudes (Md = 4.0 and 4.29) reported by NEIC. The source depth, constrained by waveforms for one of the studied events (5.5–8.5 km) seems to be considerably shallower than the hypocenter depth located by means of arrival times (∼20 km). The waveform analysis was complemented by first-motion polarities which resulted in an uncertainty assessment of the focal mechanism. Event 1 (2001-11-03) has a strike-slip mechanism with a small normal component and almost vertical nodal planes in the north-south and east-west directions. The north-south nodal plane could be related to the Calabozos faults system. Event 2 (2002-02-16) has a strike-slip mechanism with a small thrust component. The latter event (its subhorizontal nodal plane) could be associated with the El Diablo-El Fierro fault system. Dextral strike-slip solutions are consistent with recent studies in the area.     

Near-surface earthquakes in the Baltic shield, notably Sweden

A recent series of Swedish earthquakes at a focal depth not exceeding 2–3 km, the largest with I_O = V + (MSK scale) andM_L = 3.2 shows that relatively strong seismic activity can occur in the uppermost part of the Baltic Shield. During the last 15 years several near-surface earthquakes have occurred in this region, as indicated by recorded Rg-waves and/or macroseismic data. Many events are located along the coast of central Sweden, suggesting a seismic belt of minor, near-surface activity, which should be considered in connection with the radioactive waste storage in the Swedish bedrock. The appearance of Rg, common in seismic records of explosions and rockbursts, is not a sufficient discriminator between artificial events and earthquakes.     

Detailed fault structure of the Tarutung Pull-Apart Basin in Sumatra,Indonesia, derived from local earthquake data

The Tarutung Basin is located at a right step-over in the northern central segment of the dextral strike-slip Sumatran Fault System (SFS). Details of the fault structure along the Tarutung Basin are derived from the relocations of seismicity as well as from focal mechanism and structural geology. The seismicity distribution derived by a 3D inversion for hypocenter relocation is clustered according to a fault-like seismicity distribution. The seismicity is relocated with a double-difference technique (HYPODD) involving the waveform cross-correlations. We used 46,904 and 3191 arrival differences obtained from catalogue data and cross-correlation analysis, respectively. Focal mechanisms of events were analyzed by applying a grid search method (HASH code). Although there is no significant shift of the hypocenters (10.8 m in average) and centroids (167 m in average), the application of the double difference relocation sharpens the earthquake distribution. The earthquake lineation reflects the fault system, the extensional duplex fault system, and the negative flower structure within the Tarutung Basin. The focal mechanisms of events at the edge of the basin are dominantly of strike-slip type representing the dextral strike-slip Sumatran Fault System. The almost north–south striking normal fault events along extensional zones beneath the basin correlate with the maximum principal stress direction which is the direction of the Indo-Australian plate motion. The extensional zones form an en-echelon pattern indicated by the presence of strike-slip faults striking NE–SW to NW–SE events. The detailed characteristics of the fault system derived from the seismological study are also corroborated by structural geology at the surface.     

Double seismic zone beneath the middle of Hokkaido, Japan, in the southwestern side of the Kurile Arc

The vertical section of microearthquakes, determined accurately by using the Hokkaido University network, shows two dipping zones (the double seismic zone) 25–30 km apart in the depth range of 80–150 km beneath the middle of Hokkaido in the southwestern side of the Kurile arc. Hypocentral distribution of large earthquakes (m_b > 4) based on the ISC (International Seismological Centre) bulletin also shows the double seismic zone beneath the same region. The hypocentral distribution indicates that the frequency of events occurring in the lower zone is four times greater than that in the upper zone. The difference in seismic activity between the two zones beneath Hokkaido is in contrast with the region beneath northeastern Honshu in the northeastern Japan arc.Composite focal mechanisms of microearthquakes and individual mechanisms of large events mainly characterize the down-dip extension for the lower zone as is observed beneath northeastern Honshu. For the upper zone, however, the stress field is rather complex and not necessarily similar to that beneath northeastern Honshu. This may be considered to indicate the influence of slab contortion or transformation in the Hokkaido corner between the Kurile and the northeastern Japan arcs.     

The 2002 Molise earthquake sequence: What can we learn about the tectonics of southern Italy?

On October 31, 2002 a M_L = 5.4 earthquake occurred in southern Italy, at the margin between the Apenninic thrust belt (to the west) and the Adriatic plate (to the east). In this area, neither historical event nor seismogenic fault is reported in the literature. In spite of its moderate magnitude, the earthquake caused severe damage in cities close to the epicenter and 27 people, out of a total of 29 casualties, were killed by the collapse of a primary school in S. Giuliano di Puglia. By inverting broadband regional waveforms, we computed moment tensor solutions for 15 events, as small as M_L = 3.5 (M_w = 3.7). The obtained focal mechanisms show pure strike-slip geometry, mainly with focal planes oriented to NS (sinistral) and EW (dextral). In several solutions focal planes are rotated counterclockwise, in particular for later events, occurring west of the mainshock. From the relocated aftershock distribution, we found that the mainshock ruptured along an EW plane, and the fault mechanisms of some aftershocks were not consistent with the mainshock fault plane. The observed stress field, resulting from the stress tensor inversion, shows a maximum principal stress axis with an east–west trend (N83°W), whereas the minimum stress direction is almost N–S. Considering both the aftershock distribution and moment tensor solutions, it appears that several pre-existing faults were activated rather than a single planar fault associated with the mainshock. The finite fault analysis shows a very simple slip distribution with a slow rupture velocity of 1.1 km/s, that could explain the occurrence of a second mainshock about 30 h after. Finally, we attempt to interpret how the Molise sequence is related to the normal faulting system to the west (along the Apennines) and the dextral strike-slip Mattinata fault to the east.     

Atlas of macroseismic maps for French earthquakes with their principal characteristics

The SIRENE macroseismic database has been utilized to draw isoseismal maps for the 140 best-documented French earthquakes, characterized by epicentral intensities of at least V (MSK) and located in all parts of the country. A study of focal depths derived from available local intensity data using an intensity versus distance decay law (Sponheuer) shows that the focal depths of most of the events considered do not exceed about 10 km. Their distribution correlates fairly well with regional dynamic geology features. A relationship is then computed between magnitude, intensity and focal distance, based on 73 instrumenta]ly recorded earthquakes (M _L between 3.3 and 6.3) and on 217 mean radius values (from 2 to 380 km) for isoseismals of intensity VIII to III (MSK). This relationship is applied to historical earthquakes contained in the database SIRENE which are characterised by their intensity only. These results are used in the evaluation as well deterministic as probabilistic of the seismic hazard on the national territory.