An anomalous high-velocity layer at shallow crustal depths in the Narmada zone, India

The Narmada zone in central India is a zone of weakness that separates the region of Vindhyan (Meso-Neoproterozoic) deposition to the north from Gondwana (Permo-Carboniferous–lower Cretaceous) deposits to the south. The reinterpretation of analogue seismic refraction data, acquired during the early 1980s, using 2-D ray-tracing techniques reveals a basement (velocity 5.8–6.0 km s⁻¹ ) topography suggesting that the Narmada zone, bounded by the Narmada North and Narmada South faults is a region of basement uplift. A layer of anomalously high velocity (6.5–6.7 km s⁻¹ ) at depths between 1.5 and 9.0 km appears to be present in the entire region. Within the Narmada zone this layer occurs at shallower depths than outside the Narmada zone. At two places within the Narmada zone this layer is at a depth of about 1.5 km. This layer cannot be considered as the top of the lower crust because in this case it should have produced large positive gravity anomalies at the shallowest parts. Instead, these parts correspond to Bouguer gravity lows. Furthermore, lower crust at such shallow depths has not been reported from any other part of the Indian shield. Therefore, this layer is likely to represent the top of a high-velocity mafic body that has different thicknesses in different places.     

Magnetotelluric interpretation of crustal and mantle structure in the Grenville Province

Summary Accurate determinations of depths and conductivities of electrical structures in shield regions are often difficult because of the inhomogeneity of the uppermost crust. A magnetotelluric (MT) station (BAT) in the Grenville Province of the Precambrian Shield in eastern Canada has been in operation since 1975 for time-dependency studies of electrical resistivity changes related to earthquakes. The MT response of the station displays low skew with small to moderate anisotropy. One-dimensional inversion of the apparent resistivity and phase reveals two well-defined conductors in the crust, one at 10 km and the second at the base of the crust. The latter has a resistivity less than 50 Ω m. These results are substantiated by three additional MT stations located up to 40 km distant.
Data from other new MT stations and from stations previously published in the literature are compared with two-dimensional computer model results and with the three-dimensional analogue scale model results of Dosso et al. While additional data for periods less than 100 s would be desirable the results from a number of the MT stations are not inconsistent with a widespread occurrence of a conducting zone at the base of the crust in the Grenville. The inversion analysis also indicates the existence of a conductor at some depth greater than 100 km with a resistivity less than 30 Ω m. This may coincide with a seismic low-velocity zone observed in the mantle under the Canadian Shield.     

Crustal thickness variations in the margins of the Gulf of California from receiver functions

Receiver functions (RFs) from teleseismic events recorded by the NARS-Baja array were used to map crustal thickness in the continental margins of the Gulf of California, a newly forming ocean basin. Although the upper crust is known to have split apart simultaneously along the entire length of the Gulf, little is known about the behaviour of the lower crust in this region. The RFs show clear P -to- S wave conversions from the Moho beneath the stations. The delay times between the direct P and P -to- S waves indicate thinner crust closer to the Gulf along the entire Baja California peninsula. The thinner crust is associated with the eastern Peninsular Ranges batholith (PRB). Crustal thickness is uncorrelated with topography in the PRB and the Moho is not flat, suggesting mantle compensation by a weaker than normal mantle based on seismological evidence. The approximately W–E shallowing in Moho depths is significant with extremes in crustal thickness of ∼21 and 37 km. Similar results have been obtained at the northern end of the Gulf by Lewis et al., who proposed a mechanism of lower crustal flow associated with rifting in the Gulf Extensional Province for thinning of the crust. Based on the amount of pre-Pliocene extension possible in the continental margins, if the lower crust did thin in concert with the upper crust, it is possible that the crust was thinned during the early stages of rifting before the opening of the ocean basin. In this case, we suggest that when breakup occurred, the lower crust in the margins of the Gulf was still behaving ductilely. Alternatively, the lower crust may have thinned after the Gulf opened. The implications of these mechanisms are discussed.     

The deep seismicity of the Tyrrhenian Sea

The study reappraises the deep seismicity of the Tyrrhenian Sea. Careful examination of the quality of reported hypocentres shows that the earthquakes define a zone dipping NW, about 200 km along strike, 50 km thick, and reaching a depth of about 500 km. The zone is slightly concave to the NW at a depth of 300 km, but, contrary to many previous reports, is not tightly concave, nor are there significant spatial gaps in the seismicity, which is effectively continuous with depth. Seismicity is, however, concentrated in the depth interval 250–300 km, where the dip of the seismic zone changes from 70° (above 250 km) to a more gentle dip of 45° at greater depths. Seven fault-plane solutions are available for the largest earthquakes in this depth interval, all of them consistent with a P -axis down the dip of the seismic zone, and all of them requiring movement on faults out of the plane of the subducting slab.
Two deep earthquakes near Naples lie well outside the main zone of activity; for one of which a fault-plane solution is available that has a P -axis not aligned with the dip of the seismic zone. The tightly concave slab-geometry favoured by other reports is supported mainly by the location of these events near Naples, which we think may represent deformation in a separate, probably shallower dipping, piece of subducted lithosphere.
The lack of shallow seismicity, and particularly of thrust faulting earthquakes, at the surface projection of the Benioff zone suggests that active subduction has ceased. Estimates of the convergence rate responsible for subduction in the last 10 Myr far exceed the present convergence rate of Africa and Eurasia, suggesting that the subduction was related instead to the stretching and thinning of the crust in the Tyrrhenian Sea.     

The crustal structure beneath the northwest fjords, Iceland, from receiver functions and surface waves

Five broad-band seismic stations were operated in the northwest fjords area of Iceland from 1996 to 1998 as part of the Iceland Hotspot project. The structures of the upper 35  km or so beneath these stations were determined by the modelling and joint inversion of receiver functions and regional surface wave phase velocities. More than 40 teleseismic events and a few regional events containing high-quality surface wave trains were used. Although the middle period passband of the seismograms is corrupted by oceanic microseismic noise, which hinders the interpretation of structural details, the inversions reveal the overall features. Many profiles obtained exhibit large velocity gradients in the upper 5  km or so, smaller zero gradients below this, and, at ~23  km depth, a zone 2–4  km thick with higher velocity gradients. The two shallower intervals are fairly consistent with the 'upper' and 'lower' crust, defined by Flovenz (1980 ). The deep zone of enhanced velocity gradient seems to correspond to the sharp reflector first reported by Bjarnason et al . (1993 ) and identified by them as the 'Moho'. However, this type of structure is not ubiquitous beneath the northwest fjords area. The distinctiveness of the three intervals is variable, and in some cases a structure with velocity gradient increasing smoothly with depth is observed. We term these two end-members structures of the first and second types respectively. Structures of the second type correlate with older areas. Substantial variation in fundamental structure is to be expected in Iceland because of the great geological heterogeneity there.     

A seismic study along the East Greenland margin from 72°N to 77°N

Around 4370 km of new seismic reflection data, collected along the East Greenland margin between 71°30'N and 77°N in 2003, provide a first detailed view of the sediment distribution and tectonic features along the East Greenland margin. After processing and converting the data to depth, we correlated ODP-Site 913 stratigraphy into the new seismic network. Unit GB-2 shows the greatest glacial sediment deposits beneath the East Greenland continental shelf. This unit is characterized by the beginning of prograding sequences and has, according to our stratigraphic correlation, a Middle Miocene age. It might have been caused by rapid changes in sea level and/or glacial erosion by an early ice sheet or glaciers along the coast. A basement high, presumably a 360 km long basement structure at 77°N–74°54'N, prevents continuous sediment transport from the shelf into the deep sea area in times before 15 Myr. The origin of this prominent structure remains speculative since no rock sample from this structure is available. Seaward dipping reflectors at the eastern flank of this structure strongly support that it is a volcanic construction and is most likely emplaced on continental or transitional crust. The compilation of sediment thickness provide an insight into the regional sediment distribution in the Greenland Basin. An average sediment thickness of 1 km is observed. The north bordering Boreas Basin has a sediment thickness of 1.8 km close to the Greenland fracture zone (GFZ).     

Anisotropy of magnetic susceptibility investigations of the St Malo dyke swarm (Brittany, France): emplacement mechanism of doleritic intrusions

A high-quality aeromagnetic survey of northern Brittany (line spacing 250  m; flight elevation 150  m) has been used to delineate the Lower Carboniferous St Malo dyke swarm in detail. The dyke swarm fans at its northern and southern ends, and is affected by N60° sinistral transcurrent faults. After restoration of these offsets, the full structure trends in a N–S direction. Small dykes are not imaged, and only one-third of the swarm is evidenced by magnetism. Gravity and magnetic modelling shows that the swarm overlies a single N–S elongated magma chamber.
  The distribution of K _max inclinations of anisotropy of magnetic susceptibility (AMS) suggests that the dykes display a fanning magma flow in section. Computed K _max inclinations are usually shallower than the measured geological dips, probably because the flow becomes more disturbed as the dyke becomes shallower. We observe that the mean susceptibility values increase when the magma flow is steeper than about 70°.
  A detailed cross-section of St Briac dyke, which is part of the St Malo dyke swarm, reveals that the main carrier of magnetization is magnetite in the centre of the dyke and magnetite + maghaemite on its rims. The emplacement of the St Malo dyke swarm could have been contemporaneous with the N60° shearing which displaced the dyke swarm by about 20  km. The dyke swarm is cut at its southern end by an E–W-orientated fault which probably acted during Late Carboniferous times.     

Seismotectonics and analysis of earthquakes from the Hindukush region

Summary. Over 80 earthquakes, exclusively from the Hindukush focal region, which were recorded at the Gauribidanur seismic array (GBA) have been used in this study. These events have similar epicentral distances and a narrow azimuthal range from GBA but varying focal depths from 10 to 240 km. A fault plane dipping steeply (75°) in the north-west direction and striking N 66° E has been investigated on the basis of the spatial distribution of earthquakes in two vertical planes through 68° E and 32° N. Short period P -wave recordings up to 30 s were processed using the adaptive cross-correlation filtering technique. Slowness and azimuthal anomalies were obtained for first arrivals. These anomalies show positive as well as negative bias and are attributed to a steep velocity gradient in the upper mantle between the 400–700 km depth range where the seismic rays have their maximum penetration. Relative time residuals between the stations of GBA owe their origin very near to the surface beneath the array. A search of the signals across the array revealed that most of the events occurring at shallower depths had complex signatures as compared to the deeper events. The structure near the source region, complicated source functions and the scattering confined to the crust—upper mantle near source are mainly responsible for the complexity of the Hindukush earthquakes as the transmission zone of the ray tubes from turning point to the recording station is practically the same.     

Shear wave anisotropy in the crust around the San Andreas fault near Parkfield: spatial and temporal analysis

We systematically analysed shear wave splitting (SWS) for seismic data observed at a temporary array and two permanent networks around the San Andreas Fault (SAF) Observatory at Depth. The purpose was to investigate the spatial distribution of crustal shear wave anisotropy around the SAF in this segment and its temporal behaviour in relation to the occurrence of the 2004 Parkfield M 6.0 earthquake. The dense coverage of the networks, the accurate locations of earthquakes and the high-resolution velocity model provide a unique opportunity to investigate anisotropy in detail around the SAF zone. The results show that the primary fast polarization directions (PDs) in the region including the SAF zone and the northeast side of the fault are NW–SE, nearly parallel or subparallel to the SAF strike. Some measurements on the southwest side of the fault are oriented to the NNE–SSW direction, approximately parallel to the direction of local maximum horizontal compressive stress. There are also a few areas in which the observed fast PDs do not fit into this general pattern. The strong spatial variations in both the measured fast PDs and time delays reveal the extreme complexity of shear wave anisotropy in the area. The top 2–3 km of the crust appears to contribute the most to the observed time delays; however substantial anisotropy could extend to as deep as 7–8 km in the region. The average time delay in the region is about 0.06 s. We also analysed temporal patterns of SWS parameters in a nearly 4-yr period around the 2004 Parkfield main shock based on similar events. The results show that there are no appreciable precursory, coseismic, or post-seismic temporal changes of SWS in a region near the rupture of an M 6.0 earthquake, about 15 km away from its epicentre.     

Crustal structure of Atlantic fracture zones–II. The Vema fracture zone and transverse ridge

Summary. The crustal structure beneath the Vema fracture zone and its flanking transverse ridge was determined from seismic refraction profiles along the fracture zone valley and across the ridge. Relatively normal oceanic crust, but with an upwarped seismic Moho, was found under the transverse ridge. We suggest that the transverse ridge represents a portion of tectonically uplifted crust without a major root or zone of serpentinite diapirism beneath it. A region of anomalous crust associated with the fracture zone itself extends about 20 km to either side of the central fault, gradually decreasing in thickness as the fracture zone is approached. There is evidence to suggest that the thinnest crust is found beneath the edges of the 20 km wide fracture zone valley. Under the fracture zone valley the crust is generally thinner than normal oceanic crust and is also highly anomalous in its velocity structure. Seismic layer 3 is absent, and the seismic velocities are lower than normal. The absence of layer 3 indicates that normal magmatic accretionary processes are considerably modified in the vicinity of the transform fault. The low velocities are probably caused by the accumulation of rubble and talus and by the extensive faulting and fracturing associated with the transform fault. This same fracturing allows water to penetrate through the crust, and the apparently somewhat thicker crust beneath the central part of the fracture zone valley may be explained by the resultant serpentinization having depressed the seismic Moho below its original depth.     

Crustal and uppermost mantle structure in southern Africa revealed from ambient noise and teleseismic tomography

Rayleigh wave phase velocity maps in southern Africa are obtained at periods from 6 to 40 s using seismic ambient noise tomography applied to data from the Southern Africa Seismic Experiment (SASE) deployed between 1997 and 1999. These phase velocity maps are combined with those from 45 to 143 s period which were determined previously using a two-plane-wave method by Li & Burke. In the period range of overlap (25–40 s), the ambient noise and two-plane-wave methods yield similar phase velocity maps. Dispersion curves from 6 to 143 s period were used to estimate the 3-D shear wave structure of the crust and uppermost mantle on an 1°× 1° grid beneath southern Africa to a depth of about 100 km. Average shear wave velocity in the crust is found to vary from 3.6 km s^–1 at 0–10 km depths to 3.86 km s^–1 from 20 to 40 km, and velocity anomalies in these layers correlate with known tectonic features. Shear wave velocity in the lower crust is on average low in the Kaapvaal and Zimbabwe cratons and higher in the surrounding Proterozoic terranes, such as the Limpopo and the Namaqua-Natal belts, which suggests that the lower crust underlying the Archean cratons is probably less mafic than beneath the Proterozoic terranes. Crustal thickness estimates agree well with a previous receiver function study of Nair et al. . Archean crust is relatively thin and light and underlain by a fast uppermost mantle, whereas the Proterozoic crust is thick and dense with a slower underlying mantle. These observations are consistent with the southern African Archean cratons having been formed by the accretion of island arcs with the convective removal of the dense lower crust, if the foundering process became less vigorous in arc environments during the Proterozoic.     

A wide-angle seismic traverse through the Variscan of southwest Ireland

A wide-angle seismic profile across the western peninsulas of SW Ireland was performed. This region corresponds to the northernmost Variscan thrust and fold deformation. The dense set of 13 shots and 109 stations along the 120  km long profile provides a detailed velocity model of the crust.
  The seismic velocity model, obtained by forward and inverse modelling, defines a five-layer crust. A sedimentary layer, 5–8  km thick, is underlain by an upper-crustal layer of variable thickness, with a base generally at a depth of 10–12  km. Two mid-crustal layers are defined, and a lower-crustal layer below 22  km. The Moho lies at a depth of 30–32  km. A low-velocity zone, which coincides with a well-defined gravity low, is observed in the central part of the region and is modelled as a Caledonian granite which intruded upper-crustal basement. The granite may have acted as a buffer to northward-directed Variscan thrusting. The Dingle–Dungarvan Line (DDL) marks a major change in sedimentary and crustal velocity and structure. It lies immediately to the north of the velocity and gravity low, and shows thickness and velocity differences in many of the underlying crustal layers and even in the Moho. This suggests a deep, pre-Variscan control of the structural development of this area. The model is compatible with thin-skinned tectonics, which terminated at the DDL and which incorporated thrusts involving the sedimentary and upper-crustal layers.     

Axial Intrusion Zone beneath the Median Valley of the Mid-Atlantic Ridge at 37°N Detected by Explosion Seismology

Summary. Four seismic refraction lines, three of which had shots every 250 m, were shot across, along and parallel to the median valley of the Mid-Atlantic Ridge at 37° N. A method has been developed for calculating the effect on the travel times of the rough sea-floor relief beneath the profiles and has been used to correct all the travel times for this effect. Most arrivals were from a main refractor of apparent velocity 5·4 to 6·3 km s⁻¹; only beyond 35 km were faster arrivals observed from an 8·09 ± 36 km s⁻¹ refractor. The main refractor corresponds in depth, at least approximately, to the top of Layer 3 of the ocean basins but its velocity is significantly less than normal for Layer 3, perhaps due to dip. A study of time residuals along two profiles across the median valley indicates the presence of a 2 to 3 km wide low velocity zone (about 3·2 km s⁻¹) beneath the median valley floor. This zone extends over the upper 2·5 km of the crust and is believed to represent a zone of intrusion through which magma passes on its way to the sea floor.     

New insight into the crust and upper mantle structure under Alaska

To better understand the seismic structure of the subducting Pacific plate under Alaska, we determined the three-dimensional P-wave velocity structure to a depth of approximately 200 km beneath Alaska using 438,146 P-wave arrival times from 10,900 earthquakes. In this study an irregular grid parameterization was adopted to express the velocity structure under Alaska. The number of grid nodes increases from north to south in the study area so that the spacing between grid nodes is approximately the same in the longitude direction. Our results suggest that the subducting Pacific slab under Alaska can be divided into three different parts based on its geometry and velocity structure. The western part has features similar to those in other subduction zones. In the central part a thick low-velocity zone is imaged at the top of the subducting Pacific slab beneath north of the Kenai Peninsula, which is believed to be most likely the oceanic crust plus an overlying serpentinized zone and the coupled Yakutat terrane subducted with the Pacific slab. In the eastern part, significant high-velocity anomalies are visible to 60–90 km depth, suggesting that the Pacific slab has only subducted down to that depth.     

Insights into the lithospheric structure and tectonic setting of the Barents Sea region from isostatic considerations

We study the tectonic setting and lithospheric structure of the greater Barents Sea region by investigating its isostatic state and its gravity field. 3-D forward density modelling utilizing available information from seismic data and boreholes shows an apparent shift between the level of observed and modelled gravity anomalies. This difference cannot be solely explained by changes in crustal density. Furthermore, isostatic calculations show that the present crustal thickness of 35–37 km in the Eastern Barents Sea is greater than required to isostatically balance the deep basins of the area (>19 km). To isostatically compensate the missing masses from the thick crust and deep basins and to adequately explain the gravity field, high-density material (3300–3350 kg m⁻³) in the lithospheric mantle below the Eastern Barents Sea is needed. The distribution of mantle densities shows a regional division between the Western and Eastern Barents and Kara Seas. In addition, a band of high-densities is observed in the lower crust along the transition zone from the Eastern to Western Barents Sea. The distribution of high-density material in the crust and mantle suggests a connection to the Neoproterozoic Timanide orogen and argues against the presence of a Caledonian suture in the Eastern Barents Sea. Furthermore, the results indicate that the basins of the Western Barents Sea are mainly affected by rifting, while the Eastern Barents Sea basins are located on a stable continental platform.     

Mechanism of formation of fold belts and associated seismic anisotropy

Summary. Fold belts form due to shortening of deep basins on oceaic and continental crust. Basins on the oceanic crust should be characterized by a pronounced seismic anisotropy in the mantle lithosphere. Deep basins on the continental crust may develop from the stretching or the destruction of the lower crust under asthenospheric upwelling. These processes can produce seismic anisotropy in both the crust and mantle lithosphere. The character of the anisotropy is different for different basin forming processes. Considerable anisotropy should also arise from compression of the crust and mantle in fold belts. The formation of fold belts produces the original seismic anisotropy in continental lithosphere.     

Three-dimensional image of the Moho undulations beneath the Gulf of Bothnia using wide-angle seismic data

The BABEL marine seismic experiment has been carried out to investigate the lithospheric structure and antecedent tectonic signatures of the Baltic Shield, including the Archaean-Proterozoic collisional structure in the northern part of the Gulf of Bothnia.
Lithospheric seismic-reflection streamer data and simultaneously recorded wide-angle reflection and refraction data collected in the Gulf of Bothnia as part of the BABEL project have been used for 3-D modelling. The distribution of land stations around the Gulf provides a good 3-D ray coverage of the PMP reflection data recorded at the eight stations in the area and allows an estimation of strikes and dips of the Moho boundary in the area. The traveltimes of reflected phases are calculated using a method that utilizes the finite-difference solution of the eikonal equation. The Moho wide-angle-reflection (PMP) traveltimes are modelled using an inversion method. A 2-D model from the Gulf of Bothnia extended into the third dimension is used as an initial model. During the inversion the velocity is kept constant and only the Moho boundary is allowed to vary. To estimate the strike of the Moho boundary and the stability of the inversion, two initial models with different strikes are examined.
The results indicate that the Moho depth in the Gulf of Bothnia undulates and has a maximum depth of 55 km in the south, rising to 42 km in the north. The Moho depth variations seem to be step-like. This change in the Moho depth coincides with the location of the presumed fossil subduction zone in the area. The crustal-thickness variations seem to be well approximated by a nearly 2-D structure striking parallel to a postulated subduction zone immediately to the south of the Skellefte area. The presence of the step at the crust/mantle boundary can be interpreted as a result of a plate-collision event at about 2 Ga.     

Microseismicity and faulting geometry in the Gulf of Corinth (Greece)

During the summer of 1993, a network of seismological stations was installed over a period of 7 weeks around the eastern Gulf of Corinth where a sequence of strong earthquakes occurred during 1981. Seismicity lies between the Alepohori fault dipping north and the Kaparelli fault dipping south and is related to both of these antithetic faults. Focal mechanisms show normal faulting with the active fault plane dipping at about 45° for both faults. The aftershocks of the 1981 earthquake sequence recorded by King et al . (1985 ) were processed again and show similar results. In contrast, the observations collected near the western end of the Gulf of Corinth during an experiment conducted in 1991 ( Rigo et al . 1996 ), and during the aftershock studies of the 1992 Galaxidi and the 1995 Aigion earthquakes ( Hatzfeld et al . 1996 ; Bernard et al . 1997 ) show seismicity dipping at a very low angle (about 15°) northwards and normal faulting mechanisms with the active fault plane dipping northwards at about 30°. We suggest that the 8–12 km deep seismicity in the west is probably related to the seismic–aseismic transition and not to a possible almost horizontal active fault dipping north as previously proposed. The difference in the seismicity and focal mechanisms between east and west of the Gulf could be related to the difference in the recent extension rate between the western Gulf of Corinth and the eastern Gulf of Corinth, which rotated the faults dipping originally at 45° (as in the east of the Gulf) to 30° (as in the west of the Gulf).     

Evidence for underthrusting beneath the Queen Charlotte Margin, British Columbia, from teleseismic receiver function analysis

The Queen Charlotte Fault zone is the transpressive boundary between the North America and Pacific Plates along the northwestern margin of British Columbia. Two models have been suggested for the accommodation of the ∼20 mm yr⁻¹ of convergence along the fault boundary: (1) underthrusting; (2) internal crustal deformation. Strong evidence supporting an underthrusting model is provided by a detailed teleseismic receiver function analysis that defines the underthrusting slab. Forward and inverse modelling techniques were applied to receiver function data calculated at two permanent and four temporary seismic stations within the Queen Charlotte Islands. The modelling reveals a ∼10 km thick low-velocity zone dipping eastward at 28° interpreted to be underthrusting oceanic crust. The oceanic crust is located beneath a thin (28 km) eastward thickening (10°) continental crust.     

Upper-crustal structure of the Benevento area (southern Italy): fault heterogeneities and potential for large earthquakes

The Benevento region is part of the southern Apennines seismogenic belt, which experienced large destructive seismic events both in historical and in recent times. The study area lies at the northern end of the Irpinia fault, which ruptured in 1980 with a M_s = 6.9 normal faulting event, which caused about 3000 casualties. The aims of this paper are to image lateral heterogeneities in the upper crust of the Benevento region, and to try to identify the fault segments that are expected to generate such large earthquakes. This work is motivated by the recognition that lithological heterogeneities along major fault zones, inferred from velocity anomalies, reflect the presence of fault patches that behave differently during large rupture episodes. In this paper, we define the crustal structure of the Benevento region by using the background seismicity recorded during 1991 and 1992 by a local seismic array. These data offer a unique opportunity to investigate the presence of structural discontinuities of a major seismogenic zone before the occurrence of the next large earthquake. The main result that we obtained is the delineation of two NW-trending high-velocity zones (HVZs) in the upper crust beneath the Matese limestone massif. These high velocities are interpreted as high-strength regions that extend for 30-40 km down to at least 12 km depth. The correspondence of these HVZs with the maximum intensity regions of historical earthquakes (1688 AD, 1805 AD) suggests that these anomalies delineate the extent of two fault segments of the southern Apenninic belt capable of generating M = 6.5⁻⁷ earthquakes. The lateral offset observed between the two segments from tomographic results and isoseismal areas is possibly related to transverse right-lateral faults.