Contemporaneous extension and compression in the Northern

Geophys. J . Int. (1997) 129,368-388
Contemporaneous extension and compression in the Northern
Apennines from earthquake fault-plane solutions
A. Frepoli and A. Amato
Istitho Nazionale di Geofisicu, Via di Vigna Murata 605, 00143 Rome, Italy
Accepted 1997 January 8. Received 1997 January 2; in original form 1996 June I
SUMMARY
The Northern Apenninic arc belongs to the deformation zones surrounding the Adriatic
plate, which behaves as a rigid block. It is a NE-verging fold-and-thrust belt, which
developed since the Miocene with the migration of an extensional-compressional pair.
Previous seismological data are roughly in agreement with this deformation picture,
although the outer compressional front was not clearly defined by the few available
solutions. In this work we confirm the existence of two adjacent zones of contemporaneous extension and compression in the Northern Apennines, defining the extent of
these two zones and their interim1 deformation better than was previously done. This
is accomplished through the analysis of 125 new focal mechanisms of earthquakes
(2.6 < Md < 4.8) recorded by the national network of the lstituto Nazionale di Geofisica
in the period 1988-1995. The two deformation zones are clearly separated and lie very
close to each other with a partial overlapping in the Emilian-Tuscan Apennines. Strikeslip events are scattered in most of the study area. Thrust-faulting earthquakes are
located in the internal part of the most recently active thrust front (middle-upper
Pleistocene). There is a general trend of the compressional P-axes of thrust-fault and
strike-slip solutions to be perpendicular to the Apenninic direction in the external zone,
while the tensional axes of the normal-fault and strike-slip solutions in the belt do not
show a uniform orientation. However, clearly evident is the -ENE direction of
extension in the peri-Tyrrhenian region, consistent with the direction of Shmlnin the
Quaternary volcanoes of this region inferred from borehole breakouts and microearthquake fault-plane solutions. We determined the orientation of the stress tensor in
the two main deformation zones by inversion of the 125 fault-plane solutions determined
in this study. We found that the external belt is characterized by NE-SW compression
(olhorizontal, oriented N45"E),while a normal-faulting regime with o3 E-W oriented
is characteristic of the Umbria-Tuscany (back-arc) region.
-
Key words: focal mechanisms, Northern Apennines, seismicity, stress distribution.
INTRODUCTION
The growth of the Apenninic belt falls within the complex
structural environment of the central Mediterranean region,
dominated by the approximately NNW-SSE convergence
between the African and Eurasian plates at a velocity of less
than 1 cm yr-'. The opening of the Tyrrhenian Basin started
with the rifting in an E-W direction of the continental lithosphere of the Corsica-Tuscany block from the Burdigalian
(17 Ma), and then developed towards the SE with the migration
of the Calabrian Arc (Malinverno & Ryan 1986). The lithospheric thinning and the opening of a back-arc basin are associated
with the migration of the foredeep and the subduction of
ancient oceanic lithosphere attached to the continental Adriatic
plate, according to Reutter (1981), Royden, Patacca &
368
Scandone (1987), Amato et ul. (1993) and Royden (1993),
among many others. The slab retreat and the opening of the
Tyrrhenian Basin have driven the eastward migration of the
'extension-compression' system deforming the Meso-Cenozoic
sedimentary unit that forms the Apenninic thrust and fold belt
(Reutter 1981; Malinverno & Ryan 1986; Philip 1987; Royden
et al. 1987; Laubscher 1988).
For the Northern Apennines we can distinguish three main
geodynamic provinces (Bigi et al. 1990) (Fig. I): (1) a backarc basin, which includes the Tyrrhenian Sea and the periTyrrhenian region of the Italian peninsula (dominated by
extension since the upper Miocene); (2) a belt-foredeep region
along the Adriatic margin characterized by crustal shortening
and westward passive subduction of the Adriatic plate (Reutter
1981; Royden et ul. 1987; Selvaggi & Amato 1992; Amato et ul.
0 1997 RAS
Extension and compression in the Northern Apennines
5.---. contour lines of the isobaths
in km of the base of Pliocene
,, - - - _
369
-
Middle Pliocene-Pleistocenethrust front
Early Pliocene thrust front
~
--
vertical, strike-slip and transform
faults (undifferentiatedage)
normal faults (undifferentiated age)
Messinian thrust front
Serravallian to Tortonian (southern Alps) and
Torionian (Apennines) thrust front
Figure 1. Main tectonic features of the Northern Apennines and isobaths of the base of the Pliocene (modified from the 'Structural model of
Italy', Bigi et a / . 1990) and epicentral map of crustal earthquakes (location errors ERH and ERZ less than 8 km) that occurred between 1988 and
1995 (data source: ING bulletin).
1993); (3) a foreland region formed by the great sedimentary
basin of the Padana Plain and the Adriatic Sea, situated
over continental lithosphere. The crustal shortening of the
Apenninic belt occurred in paroxystic phases of 0.1-1 Myr
duration separated by longer periods of breaks (Vai 1987).
This tectonic evolution is also characterized by repeated frontal
retreat. In fact, the most recent (middle-upper Pleistocene)
thrust front is localized in a more internal position, in
agreement with the distribution of the seismic activity, as
discussed later.
A detailed knowledge of the seismicity distribution enables
us to gain a better understanding of the recent tectonic
evolution and the present-day state of stress of the Northern
Apennines. Analysing the seismic activity of the past 10 years
(Fig. l ) , we note that the maximum concentration of events
approximately coincides with the highest elevations. In this
part of the belt the seismogenic brittle crust has a thickness of
about 15-20 km, whereas in the peri-Tyrrhenian region,
characterized by a higher geothermal gradient (up to 60 mW
m-', Mongelli et al. 1989) and positive Bouguer gravity
anomalies, the seismogenic crust is less than 10 km thick.
The existence of subcrustal events as deep as -90 km, which
apparently dip at 40-45" towards the SW, provides evidence
for the location of subducted Adriatic lithosphere beneath the
Northern Apennines (Selvaggi & Amato 1992). In addition,
the tomographic images in the same portion of the Apenninic
region show the existence of a high-velocity body reaching
0 1997 RAS, G J I 129, 368-388
-300 km depth, interpreted as a remnant of former oceanic
or continental lithosphere (Amato et a/. 1993).
A basic tool for the inference of the present state of deformation of a region is the analysis of earthquake fault-plane
solutions. The availability of digital data from the seismic
network of the Istituto Nazionale di Geofisica (ING) since
1988 has allowed us to conduct a systematic investigation of
the seismicity of the Northern Apennines and to gain a more
detailed knowledge of the strain release in this region. In this
paper, we present the fault-plane solutions of 125 earthquakes
that occurred in the Northern Apennines between 1988
and 1995. These original data are compared with previously
published focal mechanisms and with the recent geological
evolution of the thrust belt.
PREVIOUS S T U D I E S
Many focal mechanisms of moderate events with magnitudes
ranging between 4.0 and 5.8 have been determined for the
Northern Apennines. The available fault-plane solutions of the
period 1939-1980, calculated from arrival times and polarities
read from seismic bulletins, are reported by Gasparini,
Iannaccone & Scarpa (1985). Fig. 2 displays the best focal
mechanisms of Gasparini et al. (1985) as selected by Zoback
(1992) for the world stress map (WSM) database (quality B
and C ) , together with three centroid moment tensor (CMT)
solutions of the Parma (1983), Norcia (1979) and Perugia
370
A. Frepoli and A . Amato
Figure 2. Fault-plane solutions from Gasparini et al. (1985) (compressional quadrants in black) selected by Zoback (1992) for the world stress
map (WSM) database (quality B and C). CMT solutions (compressional quadrants in light grey) (from N to S) of the Parma (1983), Perugia
(1984) and Norcia (1979) earthquakes. Fault-plane solutions of the four largest events of the Porto San Giorgio sequence of July 1987 (compressional
quadrants in dark grey) from Riguzzi et al. (1989).
(1984) earthquakes and the focal mechanisms of four events
of the Porto San Giorgio seismic sequence of July 1987
(Riguzzi, Tertulliani & Gasparini 1989). The two latter CMT
solutions show an evident normal-fault regime along the inner
portion of the Apenninic belt, while the 1983 Parma event, as
well as other mechanisms such as the strike-slip solution of
the Ancona (1972) event, suggest the existence of active compression along the outer part of the Apenninic belt (Gasparini
et al. 1985). The Ancona strike-slip solution shows a compressional axis ( P ) with an ENE-WSW orientation, while that
of the Parma event is aligned NNW-SSE. The first-motion
polarity solution for the latter event (Haessler et al. 1988) is
quite different from the CMT solution. The focal mechanism
is still a thrust solution but with an approximately N-S P-axis.
The two more recent seismic sequences of this region, Norcia
and Perugia, have been studied by Deschamps, Iannaccone &
Scarpa (1984) and Haessler et al. (1988), respectively. The
mb = 5.9 Norcia main-shock focal mechanism calculated by
Deschamps et al. (1984) is a normal-fault-plane solution with
a small strike-slip component. The CMT solution for this
event is slightly different, with an almost pure normal-faultplane solution. The tensional axis (T) for both the first-motion
polarities solution and the CMT solution trends ENE-WSW.
The first-motion-polarities focal mechanism of the mb = 5.2
Perugia earthquake shows a normal-fault-plane solution with
a small strike-slip component. The CMT solution in this case
is also slightly different, showing an almost pure normal
solution. The T-axis is oriented in a NNE-SSW direction for
both solutions.
The seismic sequence of Porto San Giorgio of July 1987,
approximately 50 km south of Ancona, was studied by Riguzzi
et al. (1989). This sequence had four events with magnitudes
ranging between 4.0 and 4.5, for which the authors determined
fault-plane solutions, although no indications about the
reliability of these solutions are reported. The four focal
mechanisms are almost pure thrust solutions with the maximum compression axis trending approximately E-W. This
orientation is in agreement with that observed for the Ancona
strike-slip solution.
DATA SELECTION
The 176 earthquakes analysed in this study are localized in
the Northern Apennines from 42" to 45"N and from 9" to
14"E (Fig. 3). Arrival times and polarities were accurately
re-picked from digital seismic waveforms recorded by the
seismic network of the ING (Fig. 4) in the period 1988-1995.
In addition, for events that occurred in the period May-August
1994, we used arrival times and polarities recorded by stations
of the temporary teleseismic transect of the Northern
Apennines (Fig. 4), carried out in the framework of EC project
GeoModAp (contract EV5V-CT94-0464).
The magnitudes Md of these events range between 2.6 and
4.8. The smaller events (Md between 2.6 and 3.0) are mostly
0 1997 RAS, G J I 129, 368-388
Extension and compression in the Northern Apennines
37 1
Figure 3. Epicentral distribution of the 125 events with focal mechanisms selected in this study. The events are grouped in different geographic
sectors as indicated by the dotted lines (see text for details).
7
9
11
13
15
17
19
47
45
43
41
Figure 4. Station distribution map of the Italian seismic network (only central and northern Italy are shown), and of the temporary teleseismic
transect of the Northern Apennines (May-August 1994).
0 1997 RAS, GJI 129, 368-388
372
A . Frepoli and A. Amato
located in the northern part of the studied region, where the
coverage of the seismic network is denser.
The 176 events were located with the computer program
HYPOINVERSE (Klein 1989) using the gradient velocity structure
reported in Table 1. We have obtained good values of rms
residuals for most analysed earthquakes. Fig. 5 shows the
distribution of rms residuals for the 125 selected events (see
below). Most locations have rms residuals below 0.40 s.
The focal solutions were calculated with the program FPFIT
(Reasenberg & Oppenheimer 1985). From the 176 focal solutions, we selected only those with quality Qf and Q e ual to
9
A and B (see Table 2). Qfreflects the solution prediction misfit
to the polarity data Fj (Fj=O.O represents a perfect fit to the
data, while Fj= 1.0 represents a perfect misfit). Q, summarizes
the three parameter uncertainties As, Ad and Ar (ranges of
perturbations to strike, dip and rake, respectively). 5 1 solutions
with one or both quality factors equal to C were rejected. The
remaining 125 fault-plane solutions (Fig. 6) have quality factors
AA, AB, BA and BB. Only 10 solutions were obtained with
less than eight polarities, but we accepted them because they
are well distributed on the focal sphere. The spacing between
Table 1. Gradient velocity model adopted for
locations and focal mechanism computation.
Top of layer (km)
Velocity (km s - l )
5.2
0.0
6.5
15.0
6.7
25.0
30.0
7.6
8.1
31.0
Figure 5. Distribution of rms residuals for the 125 selected
earthquakes.
Table 2. Fault-plane-solutionquality factors Qr
and Q,. Fj is the solution prediction misfit to
the polarity data. As, Ad and Ar are ranges of
perturbations to the strike, dip and rake,
respectively (see text).
A
B
C
Qf
QP
Fj50.025
0.025 < Fj5 0.1
Fj>O.l
As, Ad, Ar 5 20"
20" to 40"
> 40"
the stations of the seismic network, which is on average around
25-30 km, may introduce errors in location, especially for the
hypocentral depth. Different hypocentral locations for an event
could produce a change in the orientation of the two nodal
planes of the focal solution, although we have generally
observed that the orientations of the principal axes (what we
are mainly concerned with in this paper) do not change
dramatically. In fact, discarding the 51 earthquakes with one
or both quality factors equal to C eliminated the most unstable
solutions and thus guaranteed a more robust data set to work
on. Among the 125 selected earthquakes, only 15 per cent have
Q, = B; the remaining 85 per cent have Q,= A (Table 6), which
means the estimated uncertainties on strike, dip and rake are
less than 20".
Nonetheless, in order to test further the 'tectonic' significance
of the solutions, we selected 12 of the 125 events and analysed
them more carefully. Among these 12 selected events, the only
three earthquakes with hypocentral depths fixed at 5 km
(numbers 16, 23 and 83) are included. The focal solutions
obtained for each event fixing the hypocentral depth at the
deepest and at the shallowest values, considering the vertical
standard error ERZ (see Table 5), are shown in Fig. 7a. We
decided to test the focal mechanisms of the three events with
fixed hypocentral depths, shifting the hypocentre to 15 km and
then to 1 km (Fig. 7a). Multiple solutions have been found for
five earthquakes and are shown in Fig. 7a (Reasenberg &
Oppenheimer 1985). In general, we observe that the earthquakes for which we have a good azimuthal coverage and
many polarity data have more stable solutions (e.g. numbers
23,27,67 and 70 in Fig. 7a), whereas those with large azimuthal
gaps show the largest variations (numbers 54, 83 and 108).
Also these three examples, however, show that one of the two
principal axes (T-axes in all cases) does not change substantially
in trend and plunge (less than -20"). In some other cases
(numbers 5, 16, 84, 87 and 115) we observe that the shift
imposed on the hypocentre determination ( Z fERZ) changes
to some extent the polarity distribution on the focal sphere
(changing the take-off angles for some critical distances), and
consequently the orientation of one of the two nodal planes.
Also in these examples, the tectonic interpretation of such
earthquakes would not change much, at least as far as we base
our interpretation on the distributions of P- and T-axes.
The take-off angles from the source to the receivers depend
on the velocity model adopted for locating the earthquakes.
For this reason we also tested the influence of using different
parametrization of the medium on the stability of the faultplane solutions. We located some of the earthquakes with a
layered velocity model (the one routinely used at ING, see
Table 3) and compared the focal mechanisms obtained with
this model with those computed with the gradient velocity
model shown in Table 1. In general, although the polarity
distribution on the focal sphere changes to some degree (see
for instance the southern sector of the sphere for solution
number 71, and solution number 102 in Fig. 7b), the difference
between the two solutions is small. In some other cases, we
observe larger variations in the geometry of the nodal planes,
but a rather stable orientation of the principal axes (for
instance the P-axis of solution number 47).
FAULT-PLANE-SOLUTION CATEGORIES
To describe the distribution of different styles of deformation
within the studied region, we classify our fault-plane solutions
0 1997 RAS, GJI 129, 368-388
Extension and compression in the Northern Apennines
0 0@ @
373
0
1
a
3
2
4
a
Qa
QJ
@
@
@
Q
0
@
7
6 ++
++
11
9
14
16
aQ
@ @
Q Q
@@
25
26
27
@
@@
13
0
19
ttf +
31
12
17
15
+f+
o +
0
20
T++
21
33
32
22
J:
23
24
0
30
28
@ @ @
34
35
36
P O
37
Figure 6 . Schmidt lower-hemisphere projection of the 125 selected fault-plane solutions of this study. Compression and dilatation polarities are
indicated with crosses and circles, respectively.
according to the stress categories of Zoback (1992) for the
WSM (Table 4), which is based on different plunge ranges for
P- and T-axes. This classification is founded on the assumption
that the earthquake focal mechanisms reflect the state of stress
of the region, implying that the compression and tension axes,
to a first approximation, correspond to the principal stress
axes o1 and 03,respectively. It is worth recalling that for pure
normal- or reverse-faulting focal mechanisms the assumption
that the T- and P-axes correspond with o3 and crl is allowed,
particularly if we look at their horizontal distributions.
The 125 solutions are divided as follows: 27 belong to the
pure normal-fault category (NF), three to the normal-fault
with small strike-slip-component category (NS), 40 are
included in the strike-slip category (SS), nine belong to the
thrust-fault with small strike-slip-component category (TS),
and 20 to the pure thrust-fault category (TF). The remaining
26 focal solutions do not fall into any of these categories.
These solutions are included in the unknown-stress-field category (U). For these events the orientations of the P- and Taxes are not described in the next section. Plunges of the axes
0 1997 RAS, G J I 129, 368-388
of four solutions fall outside the ranges defined in Table 4 with
differences of only one or two degrees. We have decided to
assign these solutions to the appropriate category, in this case
two to the NF category and the other two to the TS one. The
categories assigned to these four solutions are given in brackets
in Table 6.
RESULTS
To simplify the characterization of the Northern Apennine
region in terms of areas with different focal-mechanism categories, we have assembled the T-axes of the solutions belonging
to the N F and NS categories in Fig. 8(a), and the P-axes for
the TS and T F categories in Fig. 8(b). T- and P-axes for SS
solutions are shown in Figs 8(c) and (d), respectively.
The most important result achieved with these new data is
the evident separation between an area characterized by normal
(and strike-slip) faulting and one characterized by reverse (and
strike-slip) faulting (see Figs 8a, 8b, 9a, 9b and 9c). Solutions
for N F and NS are distributed over a large portion of the
374
a 0E3 @
0
f
j
J
a
@
A. Frepoli and A . Amato
+
O
+&+
#
't
46
43
*
47
+
45
44
48
Q
0
+%++
52
49
51
53
54
0Q @ @ @ Q
@ 0@
0
@J- Q
Q 00Q
@ Q
@ 0@
56
57
55
o
58
60
59
+
61
63
66
62
00
68
69
67
s
73O
70
PO
71
++
72
78
75
77 O
74
@
83
a
84
Figure 6. (Continued.)
region under study (Fig. 9a), while the reverse solutions are
confined in a narrow belt running from the Emilian region to
the southern Marche (Figs 8b and 9b), which follows the
curvature of the Northern Apenninic arc. The SS solutions are
spread almost everywhere except in central and southern
Tuscany (Fig. 9b). Events with TF, TS and SS solutions,
indicating compression perpendicular to the arc (Figs 9b and
c), are in an inner position with respect to the most recently
active thrust front (middle-upper Pleistocene) (Vai 1987)
(Fig. 1). The solutions belonging to the U category are mostly
distributed in the northern sector (Fig. 9d).
We grouped focal mechanisms in different geographic sectors
(Fig. 3) following both the distribution of event clusters and
the distribution of fault-plane-solution categories. In the following description the solutions belonging to the U categories are
not considered.
Caorso-Cortemaggiore sector
In this sector four solutions were determined: three are T F
and one is SS. The hypocentral depths are in the upper 8 km
and the magnitudes range between 2.9 and 4.0. All the solutions
have P-axes oriented around a NE-SW direction (Figs 8b
and d).
Emilia sector
This sector also is characterized by T F and SS solutions (three
and two, respectively). The magnitudes of these events range
between 3.0 and 3.7 and the depths are between 8 and 17 km.
P-axes for these solutions have an approximately E-W direction (Figs 8b and d). Only one solution has a maximum
compression axis oriented NNE-SSW (Fig. 8b).
Romagna and Marche sector
In this sector, we do not include the numerous Cesena-Forli
focal mechanisms because they are clustered in a well-defined
area of the Romagna region. We describe this cluster of events
in the next section.
In the Romagna-Marche sector there is a predominance of
T F solutions (eight) over SS (three)and N F (two). The presence
0 1997 RAS, GJI 129, 368-388
Extension and compression in the Northern Apennines
B
87
86
@
Q @
91
0 0
92
CF';
T+
85
+ ++
93
375
88
89
90
@
@
@
94
95
96
0Q @ @ @
@ @ 652 @ @ @
@
98
99
+++
97
101
103
104
105
@o
0
100
106
102
107
10
@
@
Q
@
@ @ 69 Q
(Q
Q
110
111
109
115
116
@
121
122
++ T
+ + *
112
113
118
119
aQ
117
123
114
124
Figure 6. (Continued.)
of these two NF solutions and the nearness of the inner
extensional-regime zone allow us to suppose the existence of
a sudden change in the strain regime or a partial superimposition of the two regimes at a different depth. It is difficult to
demonstrate this last hypothesis because the present spacing
between the stations of the national seismic network does not
allow us to constrain the hypocentral depths for crustal events
well. In general, we can say that the hypocentral depths in this
sector are in the upper 20km, with magnitudes ranging
between 2.9 and 3.5. Generally, the P-axes of both the T F and
SS solutions are oriented around a NNE-SSW direction (Figs
8b and d). Only two T F solutions in the western part of this
sector have P-axes with an average direction around NW-SE
(Fig. 8b).
compute 40 fault-plane solutions, 32 per cent of the entire
data set.
Most of these events occurred in the period from October
15 to December 5, 1993. 17 and 15 events have SS and T F
solutions (Figs 10a and 1la), respectively, while only two have
N F solutions (Fig. 9a) and six belong to the unknown stress
category (Fig. 12). The depths of these events are concentrated
in the upper 15 km, with magnitudes between 2.6 and 3.7.
Generally, the T F events have a larger magnitude than the SS
events. The average magnitude for the T F is 3.3, while for the
SS it is 2.9. Taken together SS and T F solutions have P-axes
with an average orientation around the NE-SW direction
(Figs 10b and llb). The two N F solutions show a T-axis with
an ESE-WNW direction (Fig. 8a).
Cesena-Forli sector
North-western Tuscany and Frignano sector
In the analysed period the Cesena-Forli sector was affected by
a higher seismic activity in comparison with the other regions
of the area under study. For this sector it was possible to
This sector is mostly characterized by extension, with five N F
and nine SS solutions, and no T F earthquakes. Considering
both categories (SS and NF), the T-axes are approximately
0 1997 RAS, GJ1 129, 368-388
376
A . Frepoli and A. Amato
@
T +
z = 1.0
Z = 13.8
z = 3.3
54
Z = 15.0
z = 7.4
Z = 17.5
Z = 11.8
Z = 13.3
@Q
Z = 5.0
Z = 15.0
Q
Be
Z = 8.9
87
z = 17.7
(0 Q
Z = 11.6
a4
Caorso-Cortemagg.
z = 1.0
C.Metallifere
z = 1.0
@
67
83
z = 1.0
P
Z = 9.1
70
Southern Umbria
J:
Z = 5.0
27
Southern Umbria
z = 1.0
Z = 5.0
23
Northern Marche
0
@
@
0
@ '@
a
0 pjj
Z = 6.0
16
O
P+
Z = 9.6
108@
Z = 6.6
115J-f
z = 7.9
Z = 17.4
Emilia
0
Q0
z = 5.9
Cesena-Forli'
z = 9.9
Northern Abruzzi
a
z = 1.0
z = 1.0
Q
T+ +
@
Central Tuscany
Northern Umbria
T++
Z = 11.8
z = 7.4
@
Q
@
z = 1.0
Z = 16.2
Z = 13.2
Umbria
Q
+
Northern Tuscany
Z = 2.6
Figure 7. (a) Focal mechanisms tested for stability shifting the hypocentre considering the vertical standard error ERZ (Table 5). The solutions of
Fig. 6 are given in the left column; in the central and right columns the focal mechanisms (some with multiple solutions) calculated with a fixed
hypocentral depth are shown. (b) Focal mechanisms tested using two different velocity models: solutions calculated with the gradient velocity
model (Table 1j on the left; solutions calculated with the layered velocity model IPO (Table 3) on the right. For each event, numerical order, sector
name (see Table 5 ) and focal depth are indicated.
0 1997 RAS, G J I 129, 368-388
0
(b)
+
47
Z = 17.3
71
0
Extension and compression in the Northern Apennines
Colline Metallifere sector
Romagna
2 = 20.5
Cesena-Forli'
PO
377
++
This portion of central western Tuscany is characterized by
five N F solutions. There is no agreement among the orientation
of the T-axes: two are around a ENE-WSW orientation, two
are in a nearly NW-SE direction and one is in a nearly N-S
direction (Fig. 8a). The hypocentral depths are less than 13 km,
while magnitudes are between 2.9 and 3.8.
Z = 14.4
Romagna
102
Z = 12.9
Z = 7.8
125
@
@
Frignano
Z = 13.2
Z = 8.9
Figure 7. (Continued.)
Table3. Layered velocity model IPO, routinely used for bulletin locations at ING.
Velocity (km s-')
5.O
6.0
8.1
Top of layer (km)
0.0
10.0
30.0
Table 4. Fault-plane categories based on different
plunge ranges for P- and T-axes according to Zoback
(1992) for the WSM.
Plunges of axes
P
pl 2 52"
40"5 pl < 52"
pl <40"
PlI20"
PI< 20"
pl I 35"
FPS categories
T
pl535"
PI< 20"
PlI20"
pl <40"
40" I
pl< 52"
p l 2 52"
NF
NS
ss
ss
TS
TF
oriented in an E-W direction in the middle and lower Arno
valley, while in the Apuane Alps and along the TuscanyEmilia border region the orientation is around NE-SW (Figs
8a and c). In this northern part there are also two N F and
two SS solutions with T-axes oriented NW-SE. These two
SS solutions belong to the Eastern Frignano region, lie near
the external compressional Apenninic belt, and have P-axes
oriented consistently with the P-axis orientations found in this
sector (NE-SW). The magnitudes of these events range
between 2.8 and 4.1, while focal depths are in the first 20 km.
The hypocentral depths of events in the SW portion of this
sector are shallower, between 2 and 6 km.
Umbria-Marche Apenninic belt and the northern Abruzzi
sector
This sector also is characterized by extension. There are 11
and seven N F and SS solutions, respectively (no T F solutions).
The T-axes of these solutions do not show a preferential
orientation (Figs 8a and c). Only four N F and three SS
solutions show an orientation approximately perpendicular to
the trend of the Apennines. Magnitudes are between 3.0 and
4.3, while focal depths are concentrated in the first 20 km.
0 1997 RAS, G J I 129, 368-388
Peri-Tyrrhenian sector (Mt. Amiata, Vulsini and Colli
Albani)
In the Vulsini region there are two NF solutions and one SS.
For all these solutions the T-axis is oriented roughly
ENE-WSW (Figs 8a and c). The N F solution located close to
Mt. Amiata volcano has a NE-SW T-axis. In both regions
focal depths are shallow, between 2 and 6 km. Magnitudes are
between 3.0 and 3.7.
In this sector we included an earthquake that occurred in
Rome on June 12, 1995 (Md=3.8; 12 km depth). Its SS focal
mechanism is very similar to the solution calculated with the
S-wave polarization (Basili et al. 1997), showing a NE-trending
T-axis.
The N F solution of the Colli Albani earthquake of April 23,
1989 (Md=3.7; 4.0 km depth) shows a horizontal tension axis
N80"E (Amato et al. 1994). Both these SS and N F solutions
are in agreement with the general NE-SW extension direction deduced in this sector, and with the regional extensional
stress inferred from breakouts and microearthquake focal
mechanisms in the Quaternary volcano belt (Montone et al.
1995a).
STATE OF STRESS F R O M FOCALMECHANISM INVERSION
The orientations of fault planes and the directions of slip on
these planes observed from a population of focal mechanisms
can be used for the determination of the regional stress tensor.
We have applied the inverse technique of Gephart & Forsyth
(1984) for finding the stress field that is consistent with our
data set of focal mechanisms. The inversion is based on the
assumption that the deviatoric stress tensor is uniform over
the region under study. The inversion determines the four
model parameters (the three principal stress directions and the
)
measure of relative stress magnitude R = ~ J - c T ~ / c T ~ - c T ~that
minimize the misfit of a number of observations from the
predictions of the model. The misfit is given by the measure
of the smallest rotation about an axis of any orientation that
brings one of the two nodal planes with its slip direction and
sense of slip into an orientation that is compatible with the
stress model. These measures of rotation between observation
and model reflect either errors in the data or variations in the
stress field within the region under study. The four model
parameters that most closely match all the observations
are found for the two areas characterized by extension and
compression (Fig. 13).
The inversion result for the extension area gives a stresstensor orientation compatible with a normal-stress regime. The
minimum principal stress axis (c3)is horizontal and close to
the E-W direction, while the maximum principal stress axis
(ol)is vertical. In this inversion 54 focal mechanisms were
used. The relatively high average misfit value (13.8") is probably
378
A. Frepoli and A. Amato
Table 5 . Hypocentral parameters, magnitudes, location azimuthal gaps, location errors and geographic sector indications of the 125 events with
selected fault-plane solutions
Orig. time
Lat
Lon
Z
Md
ERH ERZ Sector
Date
No.
Gap rms
Emilia
6.2
174 0.25 2.1
07:35
44 28 09 56 3.1 3.2
880816
1
Colli Albani
1.5
41 46 12 40 4.0 3.7
200 0.09 1.6
23 :32
890423
2
2.0
Cesena-Forli
139 0.23 1.3
44 08 12 06 7.9 3.1
03 :03
891104
3
Cesena-Forli
44 08 12 06 5.9 3.3
139 0.23 1.4 2.0
03:14
891104
4
Northern Marche
7.8
164 0.42 3.8
43 37 12 27 6.0 3.3
11:30
900126
5
Northern Umbria
2.7
0.29 2.1
43 19 12 17 19.4 3.1
123
05:27
900224
6
Emilia
98 0.28 1.0 2.6
44 42 10 02 9.0 3.7
21 : 38
900411
7
Northern Umbria
2.3
0.36 1.2
103
43 38 12 07 1.6 3.7
22:33
8
900508
C. Metallifere
5.6
188 0.32 3.3
43 18 10 40 2.1 3.5
05 : 52
900624
9
Southern Umbria
5.7
0.45 1.8
42 45 12 45 10.0 3.6
119
02:59
900912
10
Emilia
108 0.28 1.4 3.3
19 : 58
44 56 10 54 17.0 3.1
901016
11
Roma gna
1.2
122 0.22 1.8
44 02 11 44 9.4 3.2
18:43
901124
12
Southern Umbria
124 0.25 1.0 2.1
42 36 12 55 17.5 3.0
20:14
910110
13
Romagna
0.9
173 0.01 2.2
17:55
43 50 12 02 9.6 3.0
910114
14
3.3
Umbria
140 0.40 2.9
42 58 12 51 5.9 3.0
23 :46
910114
15
Southern Umbria
136 0.31 1.0 42 41 12 58 5.0 3.1
03 : 29
910331
16
4.6
Northern Tuscany
103
0.37 1.2
02:32
43 41 10 58 2.0 3.3
910518
17
C. Metallifere
150 0.24 1.0 4.8
43 12 11 02 3.3 3.4
00:12
910530
18
Southern Umbria
3.7
170 0.49 2.6
42 54 12 51 1.6 3.6
15:48
910712
19
1.7 Cesena-Forli
127 0.31 0.9
44 07 12 09 11.2 3.2
02:25
911026
20
Cesena-Forli
2.3
146 0.19 0.8
44 07 12 08 7.8 3.0
14:48
911030
21
CaorseCortemag.
3.1
0.30 0.7
44 57 10 01 4.5 4.0
101
09:31
911031
22
Southern Umbria
151
0.33 1.6
21 :45
42 50 13 21 5.0 3.1
911112
23
3.4 Northern Umbria
1.5
151
0.22
4
3
53
17.7
3.0
20:00
12
02
911122
24
Frignano
4.8
44 22 10 31 5.6 3.6
157 0.34 3.9
12:25
911208
25
Northern Marche
1.6
0.34 3.7
225
22:32
43 38 13 04 9.3 3.4
911215
26
CaorseCortemag.
4.1
98 0.31 1.1
10:12
44 57 09 58 3.3 3.4
920101
27
CaorseCortemag.
3.2
12:30
44 58 09 58 1.3 3.2
90 0.27 0.8
920102
28
Vulsini
2.0
42 41 11 58 5.7 3.0
212
0.08 1.1
22:20
920207
29
4.1 Vulsini
0.10
1.5
42 47 12 03 3.4 3.2
203
23 : 37
920207
30
Vulsini
6.0
195 0.33 5.8
42 44 11 55 1.8 3.2
06: 50
920208
31
Romagna
2.3
0.21 1.3
159
44 14 11 56 6.7 3.1
17:10
920220
32
Cesena-Forli
3.3
44 02 12 12 0.7 3.2
204 0.47 2.0
14:56
920224
33
Romagna
1.2
100 0.27 0.9
44 25 12 01 14.1 3.5
21:23
920403
34
Northern Marche
2.4
200 0.30 3.6
43 39 12 40 3.5 3.1
08:19
920412
35
Frignano
1.0
101
0.13 1.1
44 19 10 50 23.1 3.8
11:59
920417
36
Cesena-Forli
1.6
157 0.15 0.6
44 06 12 10 1.0 3.1
19:23
920510
37
Roma gna
2.3
91
0.12 0.9
44 08 11 31 9.1 3.5
18:13
920527
38
Romagna
5.7
93
0.20 0.9
44 09 11 27 1.8 3.4
18:18
920527
39
Mt. Amiata
3.6
112 0.30 1.4
42 58 11 30 6.5 3.7
09 : 28
920612
40
Frignano
2.0
44 21 10 19 2.1 3.6
144 0.22 1.3
05: 38
920614
41
2.0
Frignano
105
0.30 1.4
44 24 10 18 5.0 3.4
06: 05
920614
42
Romagna
4.2
92 0.40 1.3
44 20 11 31 2.6 3.3
19:43
920702
43
Romagna
82 0.33 1.4 2.9
44 15 11 29 18.3 3.4
10:29
920816
44
Northern Abruzzi
1.1
148 0.13 0.9
42 26 13 21 11.0 3.9
02: 25
920825
45
Romagna
1.5
86
0.35 1.2
43 58 11 55 18.6 3.2
22:18
920829
46
Romagna
105 0.29 1.0 0.8
43 58 11 56 17.3 2.9
22:36
921023
47
Northern Tuscany
7.5
43 45 10 52 5.5 2.8
165 0.34 2.5
02:26
921028
48
0.6
Roma gna
102 0.25 0.9
11:19
44 01 11 55 19.8 3.5
921231
49
Umbria
3.9
109 0.34 1.9
10:51
43 34 12 13 12.4 4.0
930117
50
Cesena-Forli
2.3
192 0.23 1.6
44 05 12 05 9.0 2.9
22: 35
93031 1
51
Cesena-Forli
1.5
140 0.26 0.8
44 10 12 06 5.9 3.0
22: 38
930311
52
C. Metallifere
5.8
160 0.32 2.2
20:48
43 21 10 54 11.6 3.2
930320
53
C. Metallifere
8.4
199
0.30 3.4
43 20 10 54 9.1 3.0
04:59
930321
54
CMetallifere
8.2
0.30 3.5
43 20 10 53 12.5 2.9
201
08:13
930321
55
1.4 Northern Marche
129 0.15 1.2
22:55
43 36 12 36 28.0 3.2
930406
56
Caorso-Cortemag.
3.2
79
0.47 1.1
10:46
44 59 09 55 8.1 2.9
930513
57
Umbria
1.4
83
0.30 0.9
21 : 36
43 06 12 42 8.1 3.7
930604
58
Umbria
1.2
85
0.24 1.2
43 07 12 44 7.2 4.3
19:16
930605
59
5.0 C. Metallifere
43 12 10 46 10.0 3.8
184 0.45 4.5
07:51
930806
60
Southern Umbria
7.0
195 0.13 1.1
42 57 13 17 12.9 3.4
09:12
930806
61
Cesena-Forli
1.6
151
0.28 0.9
44 01 12 10 9.1 3.4
02:43
931015
62
Cesena-Forli
44 07 12 10 3.2 3.0
151
0.36 1.0 3.4
03:OO
931015
63
Cesena-Forli
1.3
0.17 1.0
44 03 12 09 9.7 2.8
23 : 52
189
931015
64
Cesena-Forli
1.2
44 12 12 17 10.3 3.2
110 0.33 1.5
01 : 58
931105
65
0 1997 RAS, GJI 129, 368-388
Extension and compression in the Northern Apennines
Table 5 . (Continued.)
No
Date
Orig. time
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
931105
931106
931107
931108
931109
931109
931109
931112
931113
931113
931113
931114
931128
931205
940202
940209
940309
940313
940325
940408
940413
940430
940516
940612
940616
940630
940705
940711
940711
940712
940719
940720
940721
940722
940722
940727
940730
941010
950211
950530
950530
950530
950530
950612
950614
950723
950823
950824
950825
950914
951010
951013
951120
951207
951214
951226
951227
951227
951228
951231
02:Ol
15:12
23:21
03: 59
13:46
13:52
19:08
01:42
09:23
12:49
16:11
08:43
10:41
15:10
10: 12
22:21
18:23
15:14
04:08
22:20
09:38
06:03
02:02
09:07
16:49
04:23
10:20
21 : 31
21 : 57
19 : 26
18:43
17:41
01:ll
04:08
10:25
20:27
03:16
17:59
04:22
03 : 36
05: 39
11:02
22:11
18:13
05:27
22:38
04:31
17:27
07: 09
14:32
06:54
19:03
02:46
20: 20
15 :26
06:15
19:59
23 :44
01 : 24
21:29
Lat
Lon
Z
Md
Gap
44 07
44 38
44 11
44 04
44 04
44 08
44 08
44 08
44 09
44 07
44 08
44 09
43 46
44 02
44 02
42 42
44 39
42 26
43 31
44 04
44 25
43 25
44 19
44 05
44 17
44 48
44 06
44 07
44 08
44 41
44 10
44 07
44 07
44 06
44 09
44 06
44 06
44 42
44 57
43 26
43 24
44 09
43 25
41 48
44 49
42 43
43 58
44 08
43 50
44 02
44 08
43 24
43 43
42 29
43 26
44 04
44 11
44 08
44 07
44 20
12 10
09 58
12 15
12 09
12 10
12 10
12 10
12 09
12 10
12 09
12 10
12 09
12 49
12 08
12 08
12 14
10 24
13 16
11 21
12 09
10 30
12 19
10 56
10 45
10 51
10 45
12 08
12 10
12 09
10 37
12 15
12 10
12 10
12 10
12 13
12 09
11 37
09 47
09 59
12 39
12 42
12 02
12 44
12 31
09 50
14 02
11 02
10 46
10 26
11 03
10 10
12 38
10 40
13 13
10 39
12 03
12 07
12 04
12 04
10 33
8.2
11.8
13.0
8.6
11.6
14.4
9.6
9.2
15.0
4.5
9.6
13.2
26.0
6.3
8.0
2.7
12.3
5.0
8.9
5.7
18.1
9.6
14.2
11.1
12.7
12.7
1.0
2.5
7.6
14.1
8.5
3.2
2.6
2.1
5.7
5.0
7.8
9.8
7.6
2.0
4.0
7.5
6.6
12.6
8.5
19.2
5.8
9.6
4.5
7.9
2.0
24.9
2.1
7.8
3.4
3.0
3.7
3.2
3.6
3.2
3.1
2.8
2.7
2.6
3.1
2.9
3.3
3.2
3.6
3.5
3.4
3.5
3.3
2.6
2.9
3.1
3.2
3.0
2.9
2.7
3.3
2.7
2.8
2.9
3.0
2.8
2.9
3.1
3.2
3.2
3.5
3.1
3.1
3.8
3.8
3.0
3.2
3.8
3.2
3.7
3.7
4.1
3.5
3.1
4.8
3.2
3.2
3.3
3.3
2.9
3.0
3.5
3.1
4.0
128
149
109
193
159
98
117
142
118
150
117
117
205
184
184
209
108
135
150
191
118
144
87
69
69
134
150
151
114
160
105
151
151
151
106
140
60
134
83
113
106
114
133
188
117
190
99
75
182
108
167
107
156
141
250
186
117
107
118
92
1.0
11.3
12.5
7.4
4.5
8.9
due to the existence of local variations in the stress orientations
within the extension area, which is quite large (-200 x 300 km).
The stress regime in the compressional belt is characterized
by a maximum principal stress (al)oriented N45"E, and
0 1997 RAS, G J I 129, 368-388
rms
0.31
0.29
0.34
0.27
0.28
0.22
0.30
0.30
0.19
0.24
0.29
0.23
0.34
0.30
0.32
0.42
0.26
0.38
0.45
0.33
0.28
0.15
0.30
0.32
0.37
0.40
0.36
0.45
0.30
0.31
0.49
0.41
0.39
0.44
0.50
0.37
0.24
0.23
0.48
0.47
0.24
0.27
0.41
0.27
0.26
0.42
0.28
0.23
0.32
0.34
0.28
0.33
0.40
0.42
0.30
0.21
0.28
0.30
0.35
0.40
ERH
0.8
3.6
1.2
1.1
1.2
0.6
0.7
0.8
0.8
0.8
0.7
0.9
4.2
1.2
1.4
4.7
1.0
1.6
2.3
1.7
1.5
1.1
1.0
0.6
1.0
1.9
0.8
1.1
0.7
2.1
2.4
1.3
1.0
1.1
1.3
0.9
0.6
1.0
1.1
2.5
0.8
0.8
4.7
0.9
0.7
3.8
1.0
0.7
1.5
1.2
1.2
1.6
1.9
1.6
6.3
1.5
1.2
0.7
0.8
0.9
ERZ
1.3
5.9
1.2
1.7
1.7
1.2
1.2
2.0
2.6
1.7
1.5
1.9
1.4
2.3
2.2
8.2
2.5
~
8.5
2.7
1.8
2.2
1.1
1.3
1.2
5.1
1.2
2.9
1.5
3.1
2.0
2.7
2.2
2.3
1.9
1.8
1.1
1.8
3.6
8.3
2.6
2.1
9.6
1.1
1.9
3.7
3.9
1.3
2.4
5.3
1.9
2.1
3.6
3.5
5.2
2.1
1.0
1.3
1.5
0.9
319
Sector
Cesena-Forli
Emilia
Cesena-Forli
Cesena-Forli
Cesena-Forli
Cesena-Forli
Cesena-Forli
Cesena-Forli
Cesena-ForIi
Cesena-Forli
Cesena-Forli
Cesena-Forli
Northern Marche
Cesena-Forli
Cesena-Forli
Vulsini
Emilia
Northern Abruzzi
Central Tuscany
Cesena-Forli
Frignano
Northern Umbria
Frignano
Northern Tuscany
Frignano
Emilia
Cesena-Forli
Cesena-Forli
Cesena-Forli
Emilia
Cesena-Forli
Cesena-Forli
Cesena-Forli
Cesena-Forli
Cesena-Forli
Cesena-Forli
Romagna
Emilia
Caorso-Cortemag.
Umbria
Umbria
Cesena-Forli
Umbria
Rome
Emilia
Southern Marche
Northern Tuscany
Northern Tuscany
Northern Tuscany
Northern Tuscany
Northern Tuscany
Umbria
Northern Tuscany
Northern Abruzzi
C. Metallifere
Cesena-Forli
Cesena-Forli
Cesena-Forli
Cesena-Forli
Frignano
oblique a2 and c3 (Fig. 13). The data in this sector are more
consistent with this stress model. In fact, although the number
of focal mechanisms is higher (71 events) than in the extension
area, the average misfit is smaller (8.3"). For both inversions
380
A. Frepoli and A. Amato
Table 6. Strikes, dips and rakes of the first nodal planes dipping to the east going clockwise from the north, azimuths and plunges
axes, fault-plane-solution categories, number of polarities and quality factors for the 125 selected focal mechanisms.
No
Date
Orig. time
Str.
Dip
T-axes
Pol. num.
Rake
P-axes
FPS categ
azim
plun.
azim. plun
1
880816
07: 35
95
85
246
42
120
33
11
- 120
U
2
0
890423
23 : 32
255
45
20
348
82
82
- 100
NF
3
03 :03
891 104
170
90
20
129
9
20
30
- 30
ss
4
891104
23 5
85
9
10
279
12
03: 14
3
- 170
ss
5
900126
20
211
11:30
55
140
11
358
30
40
TF
900224
171
13
73
05 : 27
28
210
80
13
150
ss
,6
7
900411
5
70
300
21 : 38
16
71
35
337
- 150
ss
8
62
247
6
22 : 33
45
95
19
900508
NF
350
- 50
9
54
20
35
180
14
900624
05 : 52
27 1
48
- 160
U
10
55
8
10
45
02:59
155
22
- 140
NF
9009 12
263
11
166
30
55
19 : 58
40
14
901016
17
10
ss
266
12
18:43
901124
68
63
35
13
265
9
- 60
NF
190
13
56
337
26
25
120
20: 14
11
9101 10
NF
- 140
198
14
55
77
215
17:55
9101 14
9
69
T’F
7
100
207
15
25 1
35
120
23 :46
70
9
9101 14
5
30
ss
157
16
75
1
150
03 : 29
11
910331
45
315
NF
52
- 70
17
24
910518
02:32
2
155
11
75
288
196
20
ss
18
28
46
190
12
910530
40
00: 12
80
166
120
U
19
45
910712
44
15:48
68
3
245
9
- 60
NF
144
20
90
210
45
12
911026
180
0
02 :25
90
TF
29
21
21
14:48
165
80
9
911030
299
20
6
ss
206
22
41
270
70
19
911031
136
09:31
10
140
236
TS
23
34
55
10
8
911112
235
21 : 45
- 90
55
55
NF
24
44
65
90
20 :00
8
911122
74
44
235
- 100
U
25
38
150
20
12:25
13
911208
32
44
- 170
U
250
26
66
195
30
22 : 32
9
911215
258
TF
17
36
60
27
51
240
920101
10: 12
35
9
139
150
22
18
u (TS)
28
38
260
40
920102
142
12: 30
10
27
170
28
U
29
920207
22 : 20
13
90
35
248
68
- 60
25
15
NF
30
15
255
30
23 : 37
920207
82
73
6
NF
-100
28 1
31
920208
32
170
90
137
06: 50
32
7
- 50
22
U
32
19
145
75
920220
116
17: 10
4
44
9
- 50
NS
33
20
170
90
14: 56
920224
7
129
20
- 30
30
ss
34
34
345
80
21 : 23
920403
123
18
20
40
19
ss
35
35
100
35
920412
219
08: 19
10
35
340
0
U
36
21
50
920417
11 : 59
20
115
343
- 170
U
33
238
37
46
19:23
165
85
920510
324
36
13
70
182
U
38
14
58
18: 13
185
75
920527
18
100
TF
29
176
39
36
18: 18
210
85
920527
64
6
130
28
178
U
40
6
09 :28
185
65
920612
218
- 140
NS
45
314
17
41
18
105
35
920614
250
05 : 38
14
57
11
- 40
NF
42
20
145
25
920614
06 :05
317
69
124
9
NF
- 80
43
64
175
25
920702
120
19 :43
11
22
332
120
TF
44
64
20
9208 16
10:29
121
14
105
80
25
292
TF
45
2
140
920825
02 :25
292
50
15
60
26
- 50
NF
46
40
35
920829
22: 18
95
341
U
30
22 1
13
170
47
75
921023
90
224
128
9
22:36
31
12
-150
ss
48
80
130
0
921028
02 : 26
14
265
355
-10
ss
8
49
40
11:19
115
921231
346
27
16
232
-170
U
38
50
9301 17
40
165
252
5
8
10: 51
NF
-120
22
69
51
15
65
930311
105
57
22: 35
12
269
60
TF
31
52
30
235
93031 1
132
47
22: 38
5
14
29
160
U
53
75
55
930320
77
24
20:48
22
- 120
20 1
50
u “F)
54
50
120
930321
140
04: 59
232
0
13
67
- 120
NF
55
50
120
930321
134
08: 13
3
12
235
-110
NF
74
56
65
150
1
930406
22: 55
112
21
- 30
11
38
ss
57
70
170
930513
20
10:46
301
34
9
10
ss
7
58
55
120
930604
17
21 : 36
346
246
20
- 170
ss
30
59
930605
65
135
19: 16
356
4
- 160
264
14
31
ss
60
15
165
07: 51
930806
34
36
22
240
- 150
50
u “F)
61
09: 12
930806
70
170
136
10
14
41
36
- 40
NS
62
35
205
02 :43
931015
305
66
13
51
30
22
u tTS)
63
03 :00
931015
85
100
139
- 150
9
16
24
237
ss
64
50
105
23 : 52
931015
8
344
24 1
39
160
15
ss
of P- and TQr
QP
A
A
A
A
A
B
B
A
B
B
A
A
A
A
A
B
A
A
B
B
B
A
A
A
A
A
A
A
A
A
A
A
A
B
B
B
B
A
A
A
A
A
A
B
A
B
A
B
A
B
B
A
B
B
A
A
A
A
B
A
B
B
A
A
A
A
B
A
A
A
A
A
A
A
A
A
B
A
A
B
A
A
A
A
B
A
B
A
A
B
A
B
A
A
B
A
A
A
A
A
A
A
B
B
B
A
A
A
A
A
A
B
A
A
A
A
A
A
A
A
B
A
A
A
A
A
B
A
0 1997 RAS, GJI 129, 368-388
Extension and compression in the Northern Apennines
Table 6. (Continued.)
No.
Date
Orig. time
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
I15
116
117
118
119
120
121
122
123
124
125
931105
931 105
931106
93 1107
931108
931109
931109
931 109
931112
931113
931113
931113
931114
931 128
931205
940202
940209
940309
9403 13
940325
940408
940413
940430
9405 16
940612
940616
940630
940705
94071 1
94071 1
9407 12
9407 19
940720
940721
940722
940722
940727
940730
941010
950211
950530
950530
950530
950530
950612
950614
950723
950823
950824
950825
950914
951010
951013
951120
951207
951214
951226
951227
951227
951228
951231
01 : 58
02 :01
15: 12
23 :21
03 : 59
13:46
13:52
19 :08
01:42
09 : 23
12:49
16:11
08 :43
10:41
15: 10
10: 12
22:21
18:23
15: 14
04 :08
22 :20
09:38
06 :03
02 :02
09 :07
16 : 49
04:23
10:20
21 : 31
21 : 57
19:26
18 :43
17:41
01:11
04 :08
10:25
20 :27
03 : 16
17:59
04: 22
03 : 36
05 : 39
11:02
22: 11
18: 13
05 :27
22:38
04:31
17:27
07 :09
14:32
06:54
19 :03
02 : 46
20: 20
15:26
06: 15
19:59
23 :44
01 : 24
21 : 29
0 1997 RAS, G J l 129, 368-388
Str.
95
205
105
195
225
220
250
70
185
295
185
275
185
110
170
23.5
280
75
135
260
I80
25
150
165
135
250
100
240
10
205
140
205
125
245
180
30
175
225
130
165
185
190
175
205
85
210
190
165
160
235
275
170
170
105
175
170
40
200
105
175
245
Dip
45
50
35
75
55
45
55
35
50
60
85
80
85
35
70
45
85
55
75
25
80
45
85
90
45
65
35
75
85
55
5
50
65
50
55
75
50
70
30
60
60
65
80
40
90
45
40
60
80
65
25
80
80
80
50
55
80
25
70
60
55
Rake
- 180
60
110
50
40
40
80
110
20
-50
30
120
30
0
40
80
30
50
30
- 90
30
- 140
- 120
- 30
- 170
160
170
40
- 10
20
40
- 10
170
80
50
10
- 10
140
- 160
30
-50
- 20
- 40
- 10
10
100
60
- 70
- 150
170
- 100
170
160
- 170
- 40
-110
40
80
-110
40
- 170
P-axes
azim.
220
225
270
223
77
73
257
235
48
168
224
25 1
224
350
203
62
319
102
173
260
218
133
301
25
255
28
226
274
235
57
6
81
259
252
207
254
51
191
232
22
58
60
36
87
129
22
30
25
296
9
295
305
128
239
60
299
74
27
256
207
11
plun.
30
0
11
19
1
10
9
11
15
55
16
28
16
35
10
0
16
2
9
70
13
55
42
20
35
4
30
14
10
11
41
33
10
4
2
3
33
10
47
2
55
31
34
48
7
0
8
68
28
10
69
0
6
14
53
71
18
20
60
2
30
T-axes
azim
329
317
44
335
345
328
35
9
305
267
322
125
322
229
303
328
57
195
269
79
316
238
175
124
4
120
346
16
325
328
194
336
354
12
300
163
306
91
359
290
157
328
140
333
220
116
277
150
198
104
102
35
35
329
321
184
178
220
120
300
111
plun.
30
67
73
44
51
55
77
73
39
6
24
46
24
35
41
82
24
58
31
19
28
10
33
20
24
31
40
38
3
37
48
21
24
81
58
17
21
41
29
41
6
4
18
27
7
82
69
12
13
24
20
14
21
0
6
7
34
69
22
48
17
FPS categ.
Pol. num.
Qf
U
TF
TF
TS
TS
TF
TF
TF
11
16
9
20
18
18
13
12
12
13
9
14
12
22
13
18
16
18
11
14
6
10
15
11
13
9
6
15
12
12
8
10
13
10
11
12
14
18
6
13
12
20
6
8
13
8
19
16
18
11
10
20
13
11
7
9
8
12
18
12
22
A
A
A
A
B
A
A
A
A
A
A
A
A
B
B
ss
NF
ss
U
ss
U
TS
TF
ss
TF
ss
NF
ss
NF
U
ss
U
ss
U
ss
ss
ss
U
U
ss
TF
TF
ss
U
TS
U
TS
NF
ss
ss
U
ss
TF
TF
NF
ss
ss
NF
ss
ss
ss
NF
NF
ss
TF
NF
TS
ss
A
B
A
A
A
A
A
B
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
B
B
A
B
A
A
A
A
A
A
38 1
QP
A
A
A
A
A
A
A
A
A
A
B
A
A
A
A
A
A
A
A
A
B
A
A
A
A
B
A
A
A
A
A
A
A
A
A
B
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
B
A
A
A
A
A
A
382
A . Frepoli and A . Amato
Figure 8. Distributions of P - and T-axes of the 125 selected fault-plane solutions. (a) T-axes of normal-fault solutions (NF-NS categories);
(b) P-axes of thrust-fault solutions (TF-TS); (c) T-axes of strike-slip solutions (SS); (d) P-axes of strike-slip solutions (SS).
the measure of the relative stress magnitude R is close to 0.5,
i.e. the three principal stresses are different from one another.
A more detailed study based on the same stress-inversion
technique applied to smaller subregions for the entire Italian
peninsula is presently in progress.
DISCUSSION
The Northern Apennines belong to the deforming belts that
surround the Adriatic plate. The Adriatic plate is made of
continental lithosphere with thickness increasing from 70 km
to the north of the Gargano Ridge area to more than 110 km
to the south of this area (Calcagnile & Panza 1982). This
plate behaves as a rigid, almost aseismic block, surrounded by
highly seismically active regions (Apennines, southern Alps,
Dinarides) in which seismicity is concentrated in 50,630 km
wide belts.
From slip directions of the strongest earthquakes of a 21 year
period, Anderson & Jackson (1987) proposed a counterclockwise rotation of the Adriatic plate around a pole in northern
Italy. According to Anderson & Jackson, this rotation could
-
explain the present-day deformation around the Adriatic plate,
with N-S compression in the Eastern Alps, NE-SW compression in the southern Dinarides, and parallel extension
along the Apennines. This rotation involves greater deformation in the areas more distant from the rotation pole, i.e.
the southern Dinarides and the southern Apennines, consistent
with the observations from the distribution of large earthquakes. The seismic deformation rates calculated for both the
compression in the Dinarides and the extension in the central
southern Apennines are around 2.5 mm yr-' for a short period
of observation (21 years; Anderson & Jackson 1987) and are
similar to those calculated by Jackson & McKenzie (1988)
for a longer period (70 years) using stronger earthquakes
(2.1mmyr-' in the Apennines and 1.5mmyr-' in the
Dinarides). Pondrelli, Morelli & Boschi (1995) find similar
results using strong earthquakes of the last 16 years with CMT
solutions. All these works are based on moment summation
of the largest earthquakes that occurred in the Apennines in
this century, and consider the Apennines as a continuous belt
characterized only by extension perpendicular to its length.
Jackson & McKenzie (1988) also remark that it is difficult to
0 1997 RAS, GJI 129, 368-388
Extension and compression in the Northern Apennines
383
Figure 9. Focal mechanisms of this study divided into different stress categories according to Zoback (1992): (a) normal-fault solutions (NF-NS
categories); (b) thrust-fault solutions (TF-TS); (c) strike-slip solutions (SS); (d) unknown-category fault-plane solutions (U). Thrust, strike-slip and
unknown category solutions of the Cesena-Forli sector (indicated with a dashed box) are shown in Figs 10, 11 and 12. The numbers in the first
column of Table 5 are indicated next to each focal solution.
estimate the percentage of aseismic deformation for all the
surrounding areas of the Adriatic plate. Such information is
likely to come from geodetic observations, but due to the low
deformation rates of our region it will take several years before
reliable measurements are available from space geodesy.
Contradictory results were found from VLBI and SLR data
(Ward 1994; Noomen et al. 1996). Ongoing work on old
triangulation data will probably constrain the total strain of
the Apennines better (Hunstad & England 1997).
Taking into account both the instrumental and the historical
seismicity from 1650, Westaway (1992) estimated the average
extension and compression rates in the Northern Apennines
as 0.26 mm yr-' for both, while for the southern Apennines,
characterized only by extension, he estimated a much higher
deformation rate with a maximum of 5 mm yr-' near 41"N.
The rotation of the Adriatic plate can give an explanation
for the different rates of deformation in the Northern and
southern Apennines, but does not explain the existence of a
clear separation between the northern and the southern arcs,
-
0 1997 RAS, GJI 129, 368-388
evident from geological data (Patacca & Scandone 1989),
crustal earthquake distribution (Cocco et al. 1993), subcrustal
seismicity (Selvaggi & Amato 1992) and deep structures
(Amato et al. 1993).
It is well known from chronological and analytical studies
of deformation through seismic-reflection profiles and geological data (Castellarin et al. 1986; Vai 1987) that the Northern
Apennines outer margin is characterized by recent compression, while the inner-belt portion is in extension. In
addition to this difference in style of deformation, other
elements discriminate the Northern from the southern
Apennines. Apart from the Calabrian Arc, the southern
Apennines are characterized by the lack of a clear high-velocity
anomaly in the upper -300 km, due to either a slab detachment (Spakman, van der Lee & van der Hilst 1993) or a
thinner, less rigid subducted lithosphere. Only below -300 km
depth is there evidence of a high-velocity body from tomography, while the Northern Apennines are characterized by a
remnant of subducted lithosphere in the upper 300 km accord-
384
A . Frepoli and A . Amato
20
-
20'
15'
15'
10'
10'
5'
5'
441-
44'
55,
12-
5'
10'
IS
20'
k/
I
12'
55'
25
5'
12-
5'
10'
15'
20'
15'
20'
25'
1
Cesena
55'
10'
25'
(c)
Figure 10. (a) Cesena-Forli sector strike-slip solutions; (b) P-axis directions; (c) 7'-axis directions. The numbers in the first column of Table 5 are
indicated next to each focal solution.
ing to Amato et al. (1993), whereas Spakman et al. (1993)
predicted a detachment in the upper 200 km in this region.
Moreover, subcrustal events are located only beneath the
Northern Apennines down to 90 km (Selvaggi & Amato 1992).
These differences between the two Apenninic arcs could be
connected with a lithospheric tear located in the central
Apennines.
The seismological data presented in this paper confirm that
the Northern Apenninic arc is undergoing a markedly different
evolution than the southern Apennines, and provides further
evidence of active compression along its external front. Some
previously published focal mechanisms, like the 1972 Ancona
earthquake strike-slip solution (Gasparini et al. 1985), the
thrust solutions of the 1983 Parma event (Haessler et al. 1988)
and of the 1987 Porto San Giorgio sequence (Riguzzi et al.
1989),show axes of maximum compression ( P )roughly perpendicular to the Apenninic belt. However, the sparsity of these
data did not allow the definition of the areal extent of the
compressional front.
From our results we observe two regions that are clearly
separated. One is characterized by compression along the outer
margin of the belt, and is flanked along the inner part by
another region, which is dominated by extension. There is also
a partial superimposition of the two deformation zones in the
eastern portion of the Emilian-Tuscan Apennines. From the
distribution of our data we note that there is no evidence of
an almost aseismic zone, as proposed by Lavecchia et aE. (1994).
The thrust and strike-slip solutions of the compressional
zone show maximum compression axes ( P ) generally oriented
perpendicular to the trend of the Apennines (Figs 8b, 8d). In
the analysed period we observe a concentration of events with
thrust-fault and strike-slip solutions in the Cesena-Forli sector
(Figs 10, 11 and 12).
The active compression depicted by our data (Figs 8b and d)
in the easternmost sector of the Northern Apenninic belt
apparently affects the entire brittle upper crust. Hypocentral
depths of the earthquakes are generally above 15-20km,
although they are not well constrained by the present national
network, hence no inference can be made on the stress distribution with depth, and on whether the active compression is
confined within the subducting Adriatic plate, or in the
overthrusting Apenninic crustal sheets, or in both. Along the
0 1997 RAS, GJI 129, 368-388
Extension and compression in the Northern Apennines
foredeep of the Northern Apenninic arc the Adriatic plate is
located at a depth of 7-8 km. It seems likely, therefore, that
the seismicity here also affects the underlying Adriatic plate.
Montone et al. (1992) and Montone et al. (1995b) report
results of borehole breakouts along the outer front at a depth
of 1-5 km, showing directions of Shmx perpendicular to the
thrust front, at least in the northernmost sector of the arc,
between latitudes 43"N and 45"N. This is in very good
agreement with the results of this study, and suggests that the
stress distribution is constant with depth. According to Kruse
& Royden (1994) the Adriatic plate was flexing beneath the
Apennines until middle-upper Pliocene times, and was probably uplifting during the Pleistocene due to diminished forces
acting on the subducted slab at depth. According to our results,
the compressional front is still active, implying that the shallow
processes (crustal compression) associated with the deep structures (subducting or sinking slab) are continuing today.
Another indication that confirms these results on the compression in the easternmost part of the Northern Apennines
is given by the study of recent compressional structures in
the Adriatic foreland, observed through seismic-reflection
profiles (Argnani & Frugoni 1997). The authors identified
compressional structures, elongated in an approximately
WNW-ESE direction from the central Adriatic Sea to the
Northern Apennines. From Mt. Conero (Fig. 1) to the west
these structures underlie the Northern Apenninic thrust fronts,
corresponding with a relatively higher seismicity with respect
to the off-shore part of the compressional structures.
In the internal extensional area most of the normal-fault
and strike-slip solutions show minimum compression axes (T)
oriented between E-W and NE-SW (Figs 8a and c). In some
sectors such as the Colline Metallifere and the Umbria-Marche
Apennine, the fault-plane solutions show T-axes with variable
orientations. It is possible for the latter sector that this
inconsistency of axial orientations is connected to a transition
zone between extension to the west and compression to the
east. For the Colline Metallifere, the distribution of 7'-axes
indicates radial extension, possibly due to magmatic processes
active in this region, but the poor network coverage in this
area also suggests caution in interpreting these results.
In the superimposition zone, where we observe both reverse
and normal faulting, we cannot discriminate whether the stress
varies with depth due to the previously described unreliability
of the focal depths typical of locations with regional networks.
Lavecchia et al. (1994) do not take into account this problem
in their seismotectonic zoning. The 'foothills seismic zone' that
they observe between a more internal zone under extension
and the external compressional zone along the Adriatic coast
is characterized, according to these authors, by shallow extension for events between 3 and 7 km depth (no fault-plane
solutions are indicated) and deep compression between 15 and
25 km depth, constrained only by two poorly resolved focal
mechanisms.
From neither the seismicity distribution (Fig. 1) nor the
analysis of our fault-plane solutions do we observe the presence
of transverse elements that disconnect the continuity of the
extensional and compressional zone in the Northern
Apennines, as it is described in the seismotectonic zoning of
Scandone et al. (1990). Some of the clustered distributions of
earthquakes that are visible in the seismicity maps of this
region (Fig. l), with an elongation (NE-SW) perpendicular to
the Apenninic belt, may be an artefact of the network geometry.
We verified that the location errors for events located at the
-
I
55'
10'
5'
12"
-
15'
20
25'
1
I
F55
12'
5'
LO'
is
20
2s
(b>
Figure 11. (a) Cesena-Forli sector thrust solutions; (b) P-axis directions. The numbers in the first column of Table 5 are indicated next
to each focal solution.
15'
I
\
Forli'O
5
@--
76
55'
I
12.
5
10'
15'
20'
2s
Figure 12. Cesena-Forli sector unknown category solutions. The
numbers in the first column of Table5 are indicated next to each
focal solution.
0 1997 RAS, G J I 129,368-388
385
386
A . Frepoli and A . Amato
Gr :82 1356 (pl. I az.)
b r :5 I 4 5 (pl. I az.)
01: 60 I 3 0 6
03:301138
R
:0.4
misfit : 8.3'
No. : 71
R :0.5
misfit : 13.8'
No. : 54
Extension zone:
northern Tuscany
peri-Tyrrhenian zone
Umbria-northern Abruui
(all data)
Compression zone:
Caorso-Emilia
Romagna-Cesena-Marche
(all data)
Figure 13. Stress inversion results for the extension and compression zones. The three principal stress axis directions (plunge and azimuth), the
measure of relative stress magnitudes R, the measure of misfit and the number of fault-plane solutions used for the inversion are indicated.
external boundary of the arc are at a maximum in this direction
(NE-SW) due to a lack of stations in the Adriatic region.
The N F solutions of this study are only concentrated in the
inner belt and along the peri-Tyrrhenian region (Fig. 8a). This
confirms that the Northern Apenninic belt and back-arc region
are undergoing extension, as previously proposed by Gasparini
et al. (1985) based on a few fault-plane solutions of older
events (period 1905-1980), and also suggested by the solutions
of the two recent strongest earthquakes of this area, such as
the 1979 Norcia one (Deschamps et al. 1984) and the 1984
Perugia one (Haessler et al. 1988). Both the first-motion and
the CMT solutions have a tension axis (T) with a NE-SW or
NNE-SSW direction for the Perugia event, while for the
Norcia event the T-axis is ENE-WSW oriented. The rotation
of the extensional axis between these two solutions may reflect
the bending of the tectonic features in this sector of the
Northern Apennines, whose trend changes from NW-SE to
N-S (Fig. 1). This different structural orientation is linked to
the presence of a less pronounced advancement of the outer
thrust front in the southern sector, and to Plio-Pleistocene
clockwise rotation of the foredeep in this southern portion of
the Northern Apennines (Speranza, Sagnotti & Mattei 1997).
In the peri-Tyrrhenian region our results indicate normalfault and strike-slip solutions with minimum compression axes
( T ) in an ENE-WSW direction (Figs 8a and c). The present
extensional regime along the coastal Tyrrhenian region of the
Northern central Apennines is also indicated by borehole breakout analysis and by stress directions inferred from inversion of
microearthquake (A4<4) focal mechanisms (Montone et al.
1995a). Also in this region, the stress directions deduced from
the seismicity are in general in good agreement with those of
the breakout analysis, and indicate a NE to ENE extension
direction. The Tuscania event (M=4.6), which occurred in 1971
and was studied by Gasparini et al. (1985), also shows a normalfault-plane solution with its T-axis oriented NE-SW (Fig. 2 ) .
CONCLUSIONS
The determination of 125 new focal solutions of earthquakes
that occurred in the Northern Apennines from 1988 to 1995
0 1997 RAS, GJI 129, 368-388
Extension and compression in the Northern Apennines
( M < 5) have allowed us to improve our knowledge of the state
of strain within this region. The main result of this study is
the division of the region into two areas with different styles
of deformation: the inner area of the Apenninic belt, which is
characterized by extension, and the outer margin of the belt,
which shows evidence of active compression. In the extensional
area the T-axes of normal and strike-slip solutions do not
have a homogeneous orientation, while in the compressional
area the P-axes of thrust and strike-slip solutions have a
dominant NE-SW direction. The epicentres of the compressional events are in an inner position with respect to the
more recent thrust front of the Northern Apennines. The stresstensor orientations computed for the compressional and extensional zones shows that the maximum principal stress axis in
the external front is horizontal and NE-SW oriented, while in
the internal region the minimum principal stress axis is horizontal and oriented approximately E-W. Considering the
Neogene-Quaternary tectonics of the Northern Apenninic arc
and the persistence of the extensionxompression pair throughout its evolution, our results suggest that the process driving
this mechanism, most probably the slab roll-back, is still
continuing today.
ACKNOWLEDGMENTS
We thank John Gephart who provided the computer code
relating to the stress-inversion problem. We also thank two
anonymous referees for their critical reviews. Helpful comments
were also provided by Andrea Argnani, Salvatore Barba,
Claudio Chiarabba, Massimo Cocco, Massimo Di Bona,
Francesco Frugoni, Paola Montone and Giulio Selvaggi.
This work was partially funded by EC project GeoModAp
(contract EV5V-CT94-0464). We are grateful to Enzo Boschi
for encouraging this study. One of the authors (AF) dedicates
this work to Guglielmo and Lucia.
REFERENCES
Amato. A., Alessandrini, B., Cimini, G.B., Frepoli, A. & Selvaggi, G.,
1993. Active and remnant subducted slabs beneath Italy: evidence
from seismic tomography and seismicity, Ann. G K O ~ S36,
. , 201-214.
Amato, A., Chiarabba, C., Cocco, M., Di Bona, M. & Selvaggi, G.,
1994. The 1989-1990 seismic swarm in the Alban Hills volcanic
area, Central Italy, J. Vok. Geotherm. RKS.,61, 225-237.
Anderson, H. & Jackson, J., 1987. Active tectonics of the Adriatic
region; Geophys. J. R. astr. Soc., 91, 937-983.
Argnani, A. & Frugoni, F., 1997. Foreland deformation in the Central
Adriatic and its bearing on the evolution of the Northern Apennine,
Ann. G K O ~ Ssubmitted.
.,
Basili, A., Cantore, L., Cocco, M., Frepoli, A,, Margheriti, L., Nostro,
C. & Selvaggi, G., 1997. The June 12, 1995, microearthquake
39, 1167-1175.
sequence in the city of Rome, Ann. GKU~S.,
Bigi, G., Cosentino, D., Parotto, M., Sartori, R. & Scandone, P., 1990.
Structural model of Italy, 1:500000, PFG-CNR, Quad. Ric. Sci., 114.
Calcagnile, G. & Panza, G.F., 1982. The main characteristics of the
lithosphere-asthenosphere system in Italy and surrounding regions,
Pageoph, 119, 865-879.
Castellarin, A,, Eva, C., Giglia, G. & Vai, G.B., 1986. Analisi strutturale
del Fronte Appenninico Padano, Giorn. Geol., 47, 47-76.
Cocco, M., Selvaggi, G., Di Bona, M. & Basili, A,, 1993. Recent seismic
activity and earthquake occurrence along the Apennines, in Recent
Evolution and Seismicity of the Mediterranean Region, pp. 295-312,
eds Boschi, E., Mantovani, E. & Morelli, A,, Kluwer Academic, the
Netherlands.
0 1997 RAS, GJI 129, 368-388
387
Deschamps, A., Iannaccone, G. & Scarpa, R., 1984. The Umbrian
earthquake (Italy) of 19 September 1979, Ann. Geophys., 2, 29-36.
Gasparini, C., Iannaccone, G. & Scarpa, R., 1985. Fault-plane solutions
and seismicity of the Italian Peninsula, Tectonophysics, 117, 59-78.
Gephart, J. & Forsyth, W., 1984. An improved method for determining
the regional stress tensor using earthquake focal mechanism data:
Application to the San Fernando earthquake sequence, J. geophys.
Res., 89, 9305-9320.
Haessler, H. et al., 1988. The Perugia (Italy) earthquake of 29 April
1984 a microearthquake survey, Bull. seism. Soc. Am., 78, 1948-1964.
Hunstad, I. & England, P., 1997. Estimates of strain from old Italian
triangulation network, European Geophys. Soc. X X I l General
Assembly, Vienna, 21-25 (abstract).
Jackson, J. & McKenzie, D., 1988. The relationship between plate
motions and seismic moment tensors, and the rates of active
deformation in the Mediterranean and Middle East, Geophys. J. R .
astr. Soc., 93, 45-73.
Klein, F.W., 1989. Hypoinverse, a program for Vax computers to solve
for earthquake location and magnitude; USGS. Open file rept,
89-314.
Kruse, S.E. & Royden, L.H., 1994. Bending and unbending of an
elastic lithosphere: the Cenozoic history of the Apennine and
Dinaride foredeep basins, Tectonics. 13, 278-302.
Laubscher, H.P., 1988. The arcs of the Alps and the Northern
Apennines: an updated view, Tectonophysics, 146, 67-78.
Lavecchia, G., Brozzetti, F., Barchi, M., Menichetti, M. & Keller,
J.V.A., 1994. Seismotectonic zoning in east-central Italy deduced
from an analysis of the Neogene to present deformations and related
stress fields, Geol. SOC.Am. Bull., 106, 1107-1120.
Malinverno, A. & Ryan, W.B.F., 1986. Extension in the Tyrrhenian
Sea and shortening in the Apennines as result of arc migration
driven by sinking of the lithosphere, Tectonics, 5 , 227-245.
Mongelli, F., Zito, G., Ciaranfi, N. & Pieri, P., 1989. Interpretation of
heat flow density of the Apennine chain, Italy, Tectonophysics,
164, 267-280.
Montone, P., Amato, A., Chiulli, R. & Funiciello, R., 1992. Metodologie
per la determinazione del campo di stress attuale da dati di
perforazioni profonde, Proc. 11th Mtng Gruppo Nazionak di
Geofsica della Terra Solida, pp. 337-348, CNR, Rome, Italy.
Montone, P., Amato, A., Chiarabba, C., Buonasorte, G. & Fiordelisi,
A., 1995a. Evidence of active extension in Quaternary volcanoes of
Central Italy from breakout analysis and seismicity, Geophys. RKS.
Lett., 22, 1909-1912.
Montone, P., Amato, A,, Frepoli, A. & Cesaro, M., 1995b. Active stress
and recent geologic evolution of the Italian region, EOS, Trans. Am.
Geophys. Un., 76, F622 (abstract).
Noomen, R., Springer, T.A., Ambrosius, B.A.C., Herzberger, K.,
Kuijper, D.C., Mets, G.J., Overgaauw, B. & Wakker, K.F., 1996.
Crustal deformations in the Mediterranean area computed from
SLR and GPS observations, J. Geodyn., 21, 73-96.
Patacca, E. & Scandone, P., 1989. Post-Tortonian mountain building
in the Apennines: the role of the passive sinking of a relic lithospheric
slab, In The Lithosphere in Italy: Aduances in Earth Science Research,
pp. 157-176, eds Boriani, A,, Bonafede, M., Piccardo, G.B. & Vai
G.B., Accademia Nazionale dei Lincei, Mid-Term Conference, Rome
5-6 May, 1987.
Philip, H., 1987. Plio-Quaternary evolution of the stress field in
Mediterranean zones of subduction and collision, Ann. Geophys.,
5B,301-320.
Pondrelli, S., Morelli, A. & Boschi, E., 1995. Seismic deformation in
the Mediterranean area estimated by moment tensor summation,
Geophys. J. Int., 122, 938-952.
Reasenberg, P. & Oppenheimer, D., 1985. FPFIT, FPPLOT and
FPPAGE: FORTRAN computer programs for calculating and
displaying earthquake fault-plane solutions, USGS Open f i l e Rept,
85-739.
Reutter, K.J., 1981. A trench-forearc model for the Northern
Apennines, in Sedimentary Basins of Mediterranean Margins, ed.
388
A . Frepoli and A. Amato
Wezel, F.C., pp. 433-443, CNR Italian Project of Oceanography,
Tecnoprint, Bologna.
Riguzzi, F., Tertnlliani, A. & Gasparini, C., 1989. Study of the seismic
sequence of Porto San Giorgio: Marche 3 July 1987, I1 Nuouo
Cimento, 12C, 4,453-466.
Royden, L., 1993. The tectonic expression slab pull at continental
convergent boundaries, Tectonics, 12, 303-325.
Royden, L., Patacca, E. €2 Scandone, P., 1987. Segmentation and
configuration of subducted lithosphere in Italy: an important
cbntrol on thrust-belt and foredeep-basin evolution, Geology, 15,
714-717.
Scandone, P., Patacca, E., Meletti, C., Bellatalla, M., Perilli, N. &
Santini, U., 1990. Struttura geologica, evoluzione cinematica e
schema sismotettonico della penisola italiana, Proc. G N D T
Conference, Pisa, 1, 119-135.
Selvaggi, G. & Amato, A., 1992. Subcrustal earthquakes in the
Northern Apennines (Italy): evidence for a still active subduction?,
Geophys. Res. Lett., 19, 2127-2130.
Spakman, W., van der Lee, S. & van der Hilst, R., 1993. Travel-time
tomography of the European-Mediterranean mantle down to
1400 km, Phys. Earth planet Inter., 79, 3-74.
Speranza, F., Sagnotti, L. & Mattei, M., 1997. Tectonics of the UmbriaMarche-Romagna Arc (central northern Apennines, Italy): New
paleomagnetic constraints, J. geophys. Rex, in press.
Vai, G.B., 1987. Migrazione complessa del sistema fronte deformativoavanfossa-cercine periferico: jl caso dell’Appennino settentrionale,
Mem. Sac. Geol. It., 38, 95-105.
Ward, S.N., 1994. Constraints on the seismotectonics of the central
Mediterranean from very long baseline interferometry, Geophys.
J. Int., 117, 441-452.
Westaway, R., 1992. Seismic moment summation for historical earthquakes in Italy: Tectonic implications; J. geophys. Res., 97,
15437-15 464.
Zoback, M.L., 1992. First- and second-order patterns of stress in the
lithosphere: The world stress map project, J. geophys. Rex, 97,
11703-1 1728.
0 1997 RAS, GJI 129, 368-388