Calibration of local magnitude scale for Hindukush continental subduction zone

Hindukush is an active subduction zone where at least one earthquake occurs on daily basis. For seismic hazard studies, it is important to develop a local magnitude scale using the data of local seismic network. We have computed local magnitude scale for Hindukush earthquakes using data from local network belonging to Center for Earthquake Studies (CES) for a period of three years, i.e. 2015–2017. A total of 26,365 seismic records pertaining to 2,683 earthquakes with magnitude 2.0 and greater, was used with hypocentral distance less than 600 km. Magnitude scale developed by using this data comes to be M_L = logA + 0.929logr + 0.00298r – 1.84. The magnitude determined through formulated relation was compared with that of standard relation for Southern California and relation developed by the same authors for local network for Northern Punjab. It was observed that Hindukush region has high attenuation as compared to that of Southern California and Northern Punjab which implies that Hindukush is tectonically more disturbed as compared to the said regions, hence, seismically more active as well. We have calculated station correction factors for our network. Station correction factors do not show any pattern which probably owes to the geological and tectonic complexity of this structure. Standard deviation and variance of magnitude residuals for CES network determined using Hutton and Boore scale and scale developed in this study were compared, it showed that a variance reduction of 44.1% was achieved. Average of magnitude residuals for different distance ranges was almost zero which showed that our magnitude scale was stable for all distances up to 600 km. Newly developed magnitude scale will help in homogenization of earthquake catalog. It has been observed that b-value of CES catalog decreases when magnitude is calculated by using newly developed magnitude scale.     

Exploration tree approach to estimate historical earthquakes Mw and depth,test cases from the French past seismicity

The estimation of the seismological parameters of historical earthquakes is a key step when performing seismic hazard assessment in moderate seismicity regions as France. We propose an original method to assess magnitude and depth of historical earthquakes using intensity data points. A flowchart based on an exploration tree (ET) approach allows to apply a consistent methodology to all the different configurations of the earthquake macroseismic field and to explore the inherent uncertainties. The method is applied to French test case historical earthquakes, using the SisFrance (BRGM, IRSN, EDF) macroseismic database and the intensity prediction equations (IPEs) calibrated in the companion paper (Baumont et al. Bull Earthq Eng, 2017). A weighted least square scheme allowing for the joint inversion of magnitude and depth is applied to earthquakes that exhibit a decay of intensity with distance. Two cases are distinguished: (1) a “Complete ET” is applied to earthquakes located within the metropolitan territory, while (2) a “Simplified ET” is applied to both, offshore and cross border events, lacking information at short distances but disposing of reliable data at large ones. Finally, a priori-depth-based magnitude computation is applied to ancient or poorly documented events, only described by single/sporadic intensity data or few macroseismic testimonies. Specific processing of “felt” testimonies allows exploiting this complementary information for poorly described earthquakes. Uncertainties associated to magnitude and depth estimates result from both, full propagation of uncertainties related to the original macroseismic information and the epistemic uncertainty related to the IPEs selection procedure.     

Comparison of ground motions from the 2010 Mw 7.2 El Mayor–Cucapah earthquake with the next generation attenuation ground motion prediction equations

A total of 144 free-field ground motions with closest site-to-rupture distances (R_rup) less than 200 km recorded during the 2010 M_w 7.2 El Mayor–Cucapah earthquake are used to investigate predictive capabilities of the next generation attenuation (NGA) ground-motion prediction equations (GMPE). The NGA GMPEs underpredict observed spectral accelerations at sites with shear wave velocity in the upper 30 m of the site (V_s30) between 180 and 366 m/s with R_rup from about 10 to 50 km and overpredict at sites with R_rup from about 50 to 200 km. Intra-event residuals of the NGA GMPEs exhibit a noticeable negative trend for peak ground acceleration and 0.3, 1.0, and 2.0 s periods. Comparison of the inter-event residual between the 2010 M_w 7.2 El Mayor–Cucapah earthquake and the NGA dataset reveals that short-period inter-event residuals from the 2010 M_w 7.2 El Mayor–Cucapah earthquake is within the scatter of inter-event residuals from the NGA dataset but long-period inter-event residuals do not appear within of the scatter of inter-event residuals from the NGA dataset. Spectral accelerations predicted by the NGA GMPEs are generally unbiased against V_s30 and periods of less than 4.0 s. Observed spectral accelerations show a stronger V_s30 dependence for both short and long periods compared with the NGA GMPEs. The Boore and Atkinson (Earthq Spectra 24(1):99–138, 2008) and Chiou and Youngs (Earthq Spectra 24(1):173–215, 2008) GMPEs perform better in predicting observed short-period spectral accelerations at the sites with V_s30 between 180 and 250 m/s than the Abrahamson and Silva (Earthq Spectra 24(1):67–97, 2008) and Campbell and Bozorgnia (Earthq Spectra 24(1):139–171, 2008) GMPEs.     

Estimation of site amplification,attenuation and source spectra of S-waves in the East-Central Iran

Generalized inversion of the S-wave amplitude spectra from the strong-motion network data in the East-Central Iran has been used to estimate simultaneously source parameters, site response and the S-wave attenuation (Q_s). In this regard, 190 three-component records were used corresponded to 40 earthquakes with the magnitudes M3.5–M7.3. These earthquakes were recorded at 42 stations in the hypocentral distance range from 9 to 200 km. The inverse problem was solved in 20 logarithmically equally spaced points in the frequency band from 0.4 to 15 Hz. The frequency-dependent site amplification was imposed, as a constraint, on two reference site responses in order to remove the undetermined degree of freedom in the inversion and obtain a unique inverse solution. Also, a geometrical spreading factor was assumed for removing the trade-off between geometrical spreading and anelastic attenuation. Different source parameters, such as seismic moment (M₀), seismic energy (E_s), corner frequency (f_c) and Brune stress drop (Δσ), were estimated for each event by fitting an ω² model to the spectra obtained from the inversion. The stress drop values of earthquakes, obtained in this research, are in good agreement with those of other studies. Also average site response values were correlated to the average shear wave velocities in the uppermost 30 m, in high and low frequency bands. The peak frequencies of site amplifications, estimated by the generalized inversion method, where in good agreement with those of horizontal to vertical (H/V) spectral ratios for the S-wave portion of records. However, no perfect matching in amplitude was obtained due to the deficiencies of the H/V ratio technique. By supposing a free shape for Q factor, a frequency dependent function was found, the logarithm of which could be approximated by a linear function, Q(f)=151f^0.75. The uncertainties of model parameters have been evaluated by covariance matrix of least-square fit. The residuals were also analyzed in order to assess the validity of the model. The analysis of residuals with respect to magnitude and distance indicates that they are distributed normally with approximately zero mean. The robustness of the results has been studied concerning their sensitivities to the omission of different datasets, selected randomly from original database. The results obtained here can be used in predicting ground-motion parameters applying stochastic methods.     

Ground-motion predictions in Shillong region, northeast India

We deliver ground-motion prediction equations for Shillong region, northeast India, based on a database generated by finite-fault stochastic simulations. An examination of the regional seismic source characteristics is carried out beforehand. Micro/minor earthquakes (M _W?<?5.0) nucleating at hypocentral depth <21?km in the region recorded at broadband seismic stations are observed to have Brune stress-drop ranging between 2.8 and 99.9?bars. Likewise, macroseismic intensity data for the 1897 Shillong Earthquake that nucleated at a hypocentral depth of ~35?km places the associated stress-drop at 100?200?bars. The apparent variation of the stress-drop parameter with depth is considered with two source zones namely lower-crust and upper-crust. Equations for the lower-crust predict higher ground-motion levels and exhibit affinity to those developed for stable continental region of Eastern North America. The ground-motion levels predicted by the equations for the upper-crust are relatively lower but are still higher compared with those predicted for tectonically active regions, viz., the Himalayas and Western North America.     

一种基于有限动态源的烈度估计方法

首先,基于地震烈度与震级和震中距或震源距的经验关系,考虑实际地震断层的有限性以及几何学和运动学特性,提出了一种基于有限动态源模型的烈度预估方法.然后,利用数值实验讨论了体现几何学和运动学特性的参数的作用以及断层的倾角和破裂速度的影响;最后,将这种方法应用于近年来发生的两次破坏性地震,讨论了这种方法的实用性与不足之处.在实时获取破裂过程图像成为可能的今天,用这种方法预估烈度对于震后救援十分重要.     

Intermediate-term precursors of large (<Emphasis Type="Italic">M</Emphasis> ≥ 6.6) Kamchatka earthquakes for the period from 1987 to 2004: A retrospective assessment of their information content for prediction

Data from the literature were used to systematize intermediate-term (with advance times of 1 month to ∼2.5 years) precursors to the M ≥ 6.6 Kamchatka earthquakes of 1987–2004. The precursors were observed as changes in seismological, geodetic, geophysical, water-level, and hydrochemical parameters. Retrospective assessment of the information content in these intermediate-term precursors for earthquake prediction is in progress. The focus was on estimating the occurrence times of various precursors as functions of earthquake parameters (magnitude M, hypocentral distance R, and epicenter location). In the conditions of the Kamchatka observing network, precursors can be identified by a combination of methods, mostly before M ∼ 7 earthquakes or greater south of the Kronotskii Peninsula, for which M/logR ≥ 3. It is shown that the relative proportion of earthquakes for which precursors have been identified in the observations considered here is 0.43–0.86.     

Ometepec-Pinotepa Nacional,Mexico Earthquake of 20 March 2012 (Mw 7.5): A preliminary report

An analysis of local and regional data produced by the shallow, thrust Ometepec-Pinotepa Nacional earthquake (M_w 7.5) of 20 March 2012 shows that it nucleated at 16.254°N 98.531°W, about 5 km offshore at a depth of about 20 km. During the first 4 seconds the slip was relatively small. It was followed by rupture of two patches with large slip, one updip of the hypocenter to the SE and the other downdip to the north. Total rupture area, estimated from inversion of near-source strong-motion recordings, is ~25 km × 60 km. The earthquake was followed by an exceptionally large number of aftershocks. The aftershock area overlaps with that of the 1982 doublet (M_w 7.0, 6.9). However, the seismic moment of the 2012 earthquake is ~3 times the sum of the moments of the doublet, indicating that the gross rupture characteristics of the two earthquake episodes differ. The small-slip area near the hypocenter and large-slip areas of the two patches are characterized by relatively small aftershock activity. A striking, intense, linear NE alignment of the aftershocks is clearly seen. The radiated energy to seismic moment ratios, (E_s/M₀), of five earthquakes in the region reveal that they are an order of magnitude smaller for near-trench earthquakes than those that occur further downdip (e.g., 2012 and the 1995 Copala earthquakes). The near-trench earthquakes are known to produce low A_max. The available information suggests that the plate interface in the region can be divided in three domains. (1) From the trench to a distance of about 35 km downdip. In this domain M~6 to 7 earthquakes with low values of (E_s/M₀) occur. These events generate large number of aftershocks. It is not known whether the remaining area on this part of the interface slips aseismically (stable sliding) or is partially locked. (2) From 35 to 100 km from the trench. This domain is seismically coupled where stick-slip sliding occurs, generating large earthquakes. Part of the area is probably conditionally stable. (3) From 100 to 200 km from the trench. In this domain slow slip events (SSE) and nonvolcanic tremors (NVT) have been reported.The earthquake caused severe damage in and near the towns of Ometepec and Pinotepa Nacional. The PGA exceeded 1 g at a soft site in the epicentral region. Observed PGAs on hard sites as a function of distance are in reasonable agreement with the expected ones from ground motion prediction equations derived using data from Mexican interplate earthquakes. The earthquake was strongly felt in Mexico City. PGA at CU, a hard site in the city, was 12 gal. Strong-motion recordings in the city since 1985 demonstrate that PGAs during the 2012 earthquake were not exceptional, and that similar motion occurs about once in three years.     

Earthquake ground motion predictive equations for Garhwal Himalaya,India

Predictive equations based on the stochastic approach are developed for earthquake ground motions from Garhwal Himalayan earthquakes of 3.5≤M_w≤6.8 at a distance of 10≤R≤250 km. The predicted ground motion parameters are response spectral values at frequencies from 0.25 to 20 Hz, and peak ground acceleration (PGA). The ground motion prediction equations (GMPEs) are derived from an empirically based stochastic ground motion model. The GMPEs show a fair agreement with the empirically developed ground motion equations from Himalaya as well as the NGA equation. The proposed relations also reasonably predict the observed ground motion of two major Himalayan earthquakes from Garhwal Himalayan region. For high magnitudes, there is insufficient data to satisfactorily judge the relationship; however it reasonably predicts the 1991 Uttarkashi earthquake (M_w=6.8) and 1999 Chamoli earthquake (M_w=6.4) from Garhwal Himalaya region.     

The First Version of the Amax (MW , R) Relationship for Kamchatka

—To estimate for the first time the typical relation between peak acceleration A _max?, moment magnitude M _W and hypocentral distance R for Kamchatka, 101 analog strong motion records for 1969–1993 were employed as the initial data set. Records of acceleration and velocity meters were obtained at 15 rock to medium-ground Kamchatkan sites from 33 earthquakes with M _W = 4.5–7.8, at R = 30–250?km. A _max values were determined from "true" acceleration time histories calculated by spectral deconvolution of digitized records. The maximum value over the two horizontal components was used as the A _max value in the further analysis. With the scarce data available, there were no chances to determine reliably the whole A _max?(M _W ?,?R) average surface; thus the shape of this trend surface was determined on a theoretical basis and only the level was fitted to the data. The theoretical model employed included: (1) source spectrum: according to the Brune's spectral model; (2) point-source attenuation: as 1/R plus loss specified by Q(f) = 250?f ^0.8; (3) finite-source correction for a disc-shaped incoherent source, its size depending on M _W ?; (4) accelerogram duration: including source-dependent and distance-dependent terms; (5) A _max value: based on random process representation. Distance trends calculated with this model agree with the empirical ones of FUKUSHIMA and TANAKA (1990). To calculate the absolute level for these trends, observed A _max?(M _W ?,?R) values were reduced to M _W = 8, R = 100?km using the theoretical trends as reference. The median of the reduced values, A _max?(8,?100), equal to 188?gal. was taken as the absolute reference level for the relation we sought. Note that in the process of data analysis we were forced to entirely reject relatively abundant data of two particular stations because of their prominent local amplification (×5.5) or deamplification (×0.45).     

Locations and magnitudes of earthquakes in Central Asia from seismic intensity data

We apply the Bakun and Wentworth (Bull Seism Soc Am 87:1502–1521, 1997) method to determine the location and magnitude of earthquakes occurred in Central Asia using MSK-64 intensity assignments. The attenuation model previously derived and validated by Bindi et al. (Geophys J Int, 2013) is used to analyse 21 earthquakes that occurred over the period 1885–1964, and the estimated locations and magnitudes are compared to values available in literature. Bootstrap analyses are performed to estimate the confidence intervals of the intensity magnitudes, as well as to quantify the location uncertainty. The analyses of seven significant earthquakes for the hazard assessment are presented in detail, including three large historical earthquakes that struck the northern Tien-Shan between the end of the nineteenth and the beginning of the twentieth centuries: the 1887, M 7.3 Verny, the 1889, M 8.3 Chilik and the 1911, M 8.2 Kemin earthquakes. Regarding the 1911, Kemin earthquake the magnitude values estimated from intensity data are lower (i.e. MI_LH?=?7.8 and MI_W?=?7.6 considering surface wave and moment magnitude, respectively) than the value M?=?8.2 listed in the considered catalog. These values are more in agreement with the value M _S?=?7.8 revised by Abe and Noguchi (Phys Earth Planet In, 33:1–11, 1983b) for the surface wave magnitude. For the Kemin earthquake, the distribution of the bootstrap solutions for the intensity centre reveal two minima, indicating that the distribution of intensity assignments do not constrain a unique solution. This is in agreement with the complex source rupture history of the Kemin earthquake, which involved several fault segments with different strike orientations, dipping angles and focal mechanisms (e.g. Delvaux et al. in Russ Geol Geophys 42:1167–1177, 2001; Arrowsmith et al. in Eos Trans Am Geophys Union 86(52), 2005). Two possible locations for the intensity centre are obtained. The first is located on the easternmost sub-faults (i.e. the Aksu and Chon-Aksu segments), where most of the seismic moment was released (Arrowsmith et al. in Eos Trans Am Geophys Union 86(52), 2005). The second location is located on the westernmost sub-faults (i.e. the Dzhil'-Aryk segment), close to the intensity centre location obtained for the 1938, M 6.9 Chu-Kemin earthquake (MI_LH?=?6.9 and MI_W?=?6.8).     

Earthquake parameters and spatial distribution of coseismic effects in Southern Siberia and Mongolia

In this work we review earthquakes that happened in Southern Siberia and Mongolia within the coordinates of 42°–62° N and 80°–124° E and first propose relationships between earthquake parameters (a surface-wave earthquake magnitude M _s and an epicentral intensity(I ₀) based on the MSK-64 scale) and maximal distances from an earthquake epicenter (R _{e max}), hypocenter (R _{h max}), and a seismogenic fault (R _{f max}) to the localities of secondary coseismic effects. Special attention was paid to the study of these relationships for the effects of soil liquefaction. Hence, it was shown that secondary deformations from an earthquake were distributed in space away from an earthquake epicenter, than from an associating seismogenic fault. The effects of soil liquefaction are manifested by several times closer to a seismogenic fault, than all other effects, regardless of the type of tectonic movement in a seismic focus. Within the 40 km zone from an earthquake epicenter 44% of the known manifestations of liquefaction process occurred; within the 40 km zone from a seismogenic fault—90%. We propose the next relationship for effects of soil liquefaction: M _s = 0.007 × R _{e max} + 5.168 that increases the limits of the maximum epicentral distance at an earthquake magnitude of 5.2 ≤ M _s ≤ 8.1 as compared to the corresponding relationships for different regions of the world.     

Determination ofQ S,spreading factor and S-wave acceleration by inversion of seismic intensity data

In this paper, on the basis of intensity data from 85 earthquakes occurred in China, the method of resolving overdetermined equations by using the damped least squares method is applied to inversing for the mean value of 1 s. S wave quality factor, the spreading factor and the S-wave acceleration at the foci of earthquakes in 7 regions of China. The relation between S-wave acceleration at the foci of earthquakes and magnitude is discussed. As an example, 2-DQ _S distribution in Taiwan Province is obtained by inversion. It is found that there exists some corresponding relation between this distribution and Bouguer gravitational anomaly and seismicity. Preliminary analysis and discussion on results of the inversion, and estimation of error inQ _S are made.     

Magnitude–intensity and intensity–attenuation relationships for atlas region and Algerian earthquakes

This paper presents the results of an investigation of the magnitude–intensity and intensity–attenuation relationships for earthquakes in the Atlas block and Algeria using macroseismic data. This work is based on a selected sample of isoseismal maps from 32 events which were recently revised. Surface-wave magnitudes, M_s, are recalculated using the Prague formula and range from 4·2 to 7·45. Because the Atlas mountains block is in a collision zone, earthquakes occur in general within a layer 15 km deep. Expressions of general form for the magnitude–intensity and intensity–attenuation correlations are adopted and are, respectively, and where R² = d² + h², d the source distance in km, h the focal depth in km, M_s the revised surface-wave magnitude, M_sc the predicted surface-wave magnitude, I_i the intensity at isoseismal i, I the predicted intensity, σ the standard deviation and P is zero for 50-percentile values and one for 84-percentile, and the coefficients A's and B's are determined by regression analysis. The results of this study show that the intensity–attenuation models are adequate to predict quite well the die-out of intensity with distance in the Atlas zone and coastal Algeria; it is also found that magnitude can be predicted accurately by calibrating isoseismal radii against revised instrumental surface-wave magnitude. Such magnitude–intensity relationships may be used to evaluate the magnitude of historical earthquakes in the region under survey, with no instrumental data, for which isoseismal radii and intensities are available.     

S-wave envelope broadening characteristics of microearthquakes in the Canary Islands

This study analyzes the S-wave envelope broadening characteristics of 290 earthquakes recorded by 14 stations of the Spanish National Seismograph Network in the Canary Islands region. The S-wave peak delay time (t _p ) and envelope duration (t _q ) parameters are evaluated phenomenologically to infer the strength of velocity inhomogeneities of the medium along each seismic ray path. Crustal (0?≤?h?≤?18 km) and upper mantle (18?<?h?≤?80 km) events are analyzed separately. Results in the frequency range 1 to 12 Hz for hypocentral distances from 30 to 600 km show that both t _p and t _q increase according to a power of hypocentral distance and they are independent of frequency. The spatial distribution of the peak delay time reveals weak strength of heterogeneity in most of the region at shallow depths. Relatively strong inhomogeneous zones are generated under the island of Tenerife and Gran Canaria at depths of 11–22 km.     

Moment and duration of shallow earthquakes from Love-wave modelling for regional distances

Love waves from five earthquakes in southwest Germany, Italy, Yugoslavia, Greece and Algeria, as recorded at the broadband Gräfenberg array in eastern Bavaria, are modelled in detail using an extended reflectivity method. The focal mechanisms are taken from the literature or determined from long-period WWNSS data. The theoretical double-couple moment function used for modelling increases smoothly during the risetime T from zero to the seismic moment M₀. T and M₀ are determined by matching the observed and theoretical Love-wave seismograms in the long-period WWNSS frequency band. T is approximately equal to the rupture duration; values from 1.4 to 20 s are found. Seismic moments generally agree within a factor of about two with those found by other authors from larger data sets. The earthquakes in southwest Germany and Greece were simple events and can be modelled very well by one double couple. In the case of the other earthquakes, the theoretical Love-wave seismograms are too short. For these earthquakes there is some independent evidence for a multiple-event nature. This evidence is condensed into source models consisting of two or three double couples representing separate and complete sub-earthquakes, which then explain the observations quite well. In the case of the Algerian earthquake, the long Love-wave duration may also be due to complications in structure along the wavepath across the Mediterranean and the Alps. The present results show that accurate Love-wave modelling for a single station at regional distance from an epicentre is a reliable method for estimation of seismic moment and rupture duration. Use of a digitally recording station with large dynamic range allows events over a large magnitude interval to be studied. The Gräfenberg array is well suited for such an investigation of European and North African earthquakes.     

Ground motion directionality in the 2010–2011 Canterbury earthquakes

This paper examines the observed directionality of ground motions in the Christchurch urban area during the 2010–2011 Canterbury, New Zealand earthquakes. A dataset of ground motions recorded at 20 strong motion stations over 10 different earthquake events is utilized to examine the ratios of various response spectral directionality definitions and the orientation of the maximum direction. Because the majority of previous related studies have utilized overlapping ground motion datasets from the NGA database, the results of this study provide a largely independent assessment of these ground motion aspects. It is found that the directionality ratio between the maximum (100th percentile) and 50th percentile orientation‐independent spectral acceleration is similar to that obtained from recent studies. Ground motions from the 4 September 2010 Darfield earthquake are shown to exhibit strong directionality for source‐to‐site distances up to R_rup = 30 km, notably further than results from a previous study, which suggests that such effects are generally limited to R_rup < 5 km. The adopted dataset also offers the unique potential to consider site‐specific effects on directionality ratios and maximum direction orientations; however, in both cases, site‐specific effects are found not to be significant in the observed empirical results. Copyright © 2014 John Wiley & Sons, Ltd.     

Damage-based inelastic response spectra for seismic design incorporating performance considerations

Modern seismic design allows a structure to develop inelastic response during moderate to severe earthquakes. The emerging performance-based design requires more clearly defined levels of inelastic response, or damage, to be targeted for different earthquake hazard levels. While there are a range of factors that could influence the level of damage and hence the performance, the design strength remains to be a fundamental design parameter that is inherently related to the structural performance. In this paper, the response reduction factor, which is a normalized form of the design strength, is investigated on a direct damage basis. The implications of the damage-based strength reduction factor (SRF), denoted as R_D factor, on multiple performance targets are discussed. A series of R_D spectra are generated from a large set of ground motions in different groupings to examine the effects of local site condition, earthquake magnitude and distance to rupture on the R_D spectra. The overall mean and standard deviation of the R_D spectra for different levels of damage are obtained, and simple empirical formulas are proposed.     

The French seismic CATalogue (FCAT-17)

In regions that undergo low deformation rates, as is the case for metropolitan France (i.e. the part of France in Europe), the use of historical seismicity, in addition to instrumental data, is necessary when dealing with seismic hazard assessment. This paper presents the strategy adopted to develop a parametric earthquake catalogue using moment magnitude M_w, as the reference magnitude scale to cover both instrumental and historical periods for metropolitan France. Work performed within the framework of the SiHex (SIsmicité de l’HEXagone) (Cara et al. Bull Soc Géol Fr 186:3–19, 2015. doi: 10.2113/qssqfbull.186.1.3) and SIGMA (SeIsmic Ground Motion Assessment; EDF-CEA-AREVA-ENEL) projects, respectively on instrumental and historical earthquakes, have been combined to produce the French seismic CATalogue, version 2017 (FCAT-17). The SiHex catalogue is composed of ~40,000 natural earthquakes, for which the hypocentral location and M_w magnitude are given. In the frame of the SIGMA research program, an integrated study has been realized on historical seismicity from intensity prediction equations (IPE) calibration in M_w detailed in Baumont et al. (submitted) companion paper to their application to earthquakes of the SISFRANCE macroseismic database (BRGM, EDF, IRSN), through a dedicated strategy developed by Traversa et al. (Bull Earthq Eng, 2017. doi: 10.1007/s10518-017-0178-7) companion paper, to compute their M_w magnitude and depth. Macroseismic data and epicentral location and intensity used both in IPE calibration and inversion process, are those of SISFRANCE without any revision. The inversion process allows the main macroseismic field specificities reported by SISFRANCE to be taken into account with an exploration tree approach. It also allows capturing the epistemic uncertainties associated with macroseismic data and to IPEs selection. For events that exhibit a poorly constrained macroseismic field (mainly old, cross border or off-shore earthquakes), joint inversion of M_w and depth is not possible, and depth needs to be fixed to calculate M_w. Regional a priori depths have been defined for this purpose based on analysis of earthquakes with a well constrained macroseismic field where joint inversion of M_w and depth is possible. As a result, 27% of SISFRANCE earthquake seismological parameters have been jointly inverted and for the other 73% M_w has been calculated assuming a priori depths. The FCAT-17 catalogue is composed of the SIGMA historical parametric catalogue (magnitude range between 3.5 up to 7.0), covering from AD463 to 1965, and of the SiHex instrumental one, extending from 1965 to 2009. Historical part of the catalogue results from an automatic inversion of SISFRANCE data. A quality index is estimated for each historical earthquake according to the way the events are processed. All magnitudes are given in M_w which makes this catalogue directly usable as an input for probabilistic or deterministic seismic hazard studies. Uncertainties on magnitudes and depths are provided for historical earthquakes following calculation scheme presented in Traversa et al. (2017). Uncertainties on magnitudes for instrumental events are from Cara et al. (J Seismol 21:551–565, 2017. doi: 10.1007/s10950-016-9617-1).