Vol51,2_3,2008 397 ANNALS OF GEOPHYSICS, VOL. 51, N. 2/3, April/June 2008 Postseismic deformation and body forces shaping the Apennines and adjacent sedimentary basins in Umbria-Marche Riccardo E.M. Riva (1) and Abdelkrim Aoudia (2) (1) Dept. of Earth Observation and Space Systems, Delft University of Technology, Delft, The Netherlands (2) Earth System Physics Section, the Abdus Salam International Centre for Theoretical Physics, Trieste, Italy Abstract The geodynamic complexity of the Apennines and adjacent sedimentary basins in Umbria-Marche (North-Central Italy) makes the dynamics of the present day deformation and its relationships with the seismicity less well under- stood. In this paper, we argue that, further to buoyancy forces, postseismic deformation of earthquakes taking place on the Apennines contributes to the regional deformation. We investigate the interaction between the normal fault- ing system responsible of the 1997 Umbria-Marche earthquake sequence (Colfiorito fault) and the low angle nor- mal faulting system bordering the sedimentary basins, namely the Altotiberina fault. We set-up a 2D finite element model of the lithosphere-asthenosphere accounting for lateral heterogeneities and investigate how this heteroge- neous structure is capable of localizing strain under the Umbria-Marche sedimentary basins, providing a way for the Colfiorito fault to influence the evolution of the Altotiberina fault. We show how the two different length and time scale processes, namely postseismic deformation and buoyancy, are complementary in shaping the Apennines and adjacent sedimentary basins. The high resolution deformation patterns modeled in this study can hardly be re- produced by a model accounting only for external forces such as a rotating or subducting or retreating Adria. Key words Viscoelastic relaxation – Aseismic slip – Body forces – Normal faults 1. Introduction The geodynamic complexity of the Apen- nines and adjacent sedimentary basins in Um- bria-Marche (North-Central Italy) makes the dy- namics of the present day deformation and its re- lationships with the seismicity less well under- stood. The recent 1997-1998 Umbria-Marche earthquake sequence clearly highlighted complex deformation processes both in its coseismic (e.g. Amato et al., 1998; Zollo et al., 1999; Salvi et al., 2000; Chimera et al., 2003) and postseismic (Ri- va et al., 2000; Aoudia et al., 2003; Riva et al., 2007) phases. The Umbria-Marche crustal nor- mal faulting seismic sequence, which consisted of three moderate earthquakes (26/09/97 at 00:33, MW=5.7; 26/09/1997 at 09:40, MW=6.1 and 14/10/1997 at 15:23, MW=5.6), ruptured three distinct fault segments (e.g., Galadini et al., 1999; Meghraoui et al., 1999; Barba and Basili, 2000). In the same region an earthquake occurred in the upper mantle, on March 28, 1998, MW=5.2 (e.g., Olivieri and Ekstrom, 1999; Chimera et al., 20003) highlighting the importance of deep processes in the lithosphere. The lithosphere ve- locity models beneath Central Italy clearly high- lighted the important lateral heterogeneity of the physical properties of the crust-uppermost mantle system along with pronounced variations of the Mailing address: Dr. Riccardo E.M. Riva, Dept. of Earth Observation and Spaces System, Delft University of Technology, Kluyverweg 1, 2629 HS Delft, The Netherlands; e-mail: r.e.m.riva@tudelft.nl Vol51,2_3,2008 4-03-2009 10:28 Pagina 397 398 R.E.M. Riva and A. Aoudia geometry of both the crust-mantle and litho- sphere-asthenosphere boundaries (e.g. Agostinet- ti et al., 2002; Chimera et al., 2003; Mele and Sandvol, 2003). These models exhibit clear evi- dence of lithospheric roots without any continu- ous slab. This would imply that the slab has been eroded and no engine is actually left to drive rel- ative subduction processes. Therefore our fa- vorite model of the dynamics of the ongoing con- tinental deformation is the one reported by Aoudia et al. (2007) that invokes body forces in- stead of external forces to satisfy the recent kine- matics. Aoudia et al. (2007) have shown that the buoyancy forces resulting from the heteroge- neous density distribution in the lithosphere gov- ern the present day deformation within Central Italy and can explain regional coexisting contrac- tion and extension at shallow depth and unusual distribution of sub-crustal earthquakes. Their modeled uppermost mantle flow supports the lithospheric delamination beneath the peninsula and provides a unifying background for petrolog- ic and geochemical studies of recent magmatism and volcanism in Tuscany. Further to the buoyancy forces already report- ed in Aoudia et al. (2007), this paper argues that postseismic deformation of earthquakes taking place on the Apennines contributes to the region- al deformation. We investigate the interaction be- tween the normal faulting system responsible for the 1997 Umbria-Marche earthquake sequence (Colfiorito fault, CF) and the low angle normal faulting system bordering the sedimentary basins, namely the Altotiberina fault (AF) (e.g. Boncio and Lavecchia, 2000; Chiaraluce et al., 2007). By means of a half-space model we investigate the effect on the CF of prolonged slip on the AF. Af- ter that, we set-up a 2D finite element model of the lithosphere-asthenosphere accounting for lat- eral heterogeneities and investigate how this het- erogeneous structure is capable of localizing strain under the Umbria-Marche sedimentary basins, providing a way for the CF to influence the evolution of the AF. 2. Secular postseismic relaxation When approaching the study of postseismic deformation in the area of the 1997 Colfiorito earthquake, different mechanisms have to be taken into account (e.g. viscoelastic relaxation, afterslip, poroelastic relaxation…). Here, we are not interested in the details of present day relaxation (see, e.g., Riva et al., 2007): instead, we want to estimate the cumulative impact of viscoelastic relaxation due to repeated seismic events in the area. According to historical records (Boschi et al., 2000), the largest events that occurred in the direct neighborhoods of the Colfiorito fault (CF) are the 1279 Camerino in the north and the 1328 Norcia in the south, with magnitudes 6.6 and 6.2 respectively. We model the cumulative effect of those two earthquakes by means of a 50-km-long nor- mal fault rupturing at a depth of 8 km and dip- ping 45 degrees to the SW, which approximate- ly corresponds to an extension of the main fault of the 1997 sequence by 15 km to the north and 25 km to the south. The moment magnitude is equivalent to the cumulative effect of a 6.6 and 6.2 event (M0 = 6x1018 Nm). The full relaxation of this fault is modeled by means of a normal modes approach (Sabadini and Vermeersen, 1997), with a stratified earth model comprising a 15-km elastic upper crust, a Moho at 30 km, a 20-km-thick mantle wedge and a homogeneous mantle. Since we discuss the fully relaxed mod- el (i.e. at infinite time), the viscosity profile which controls the temporal variations of the relaxation process is here irrelevant, the key earth parameter being the thickness of the upper elastic layer. Results for vertical relaxation are shown in fig. 1, where the thin solid lines describe the surface projection of the boundaries of the three main patches that represent the Altotiberina fault (AF) (Boncio and Lavecchia, 2000), start- ing at the surface in the SW and dipping to- wards the NE, and smaller open rectangles rep- resent the surface projections of the three main faults of the 1997 sequence (from the north: the 9:40 event on 26 Sept., the 0:33 event on 26 Sept, and the Sellano event on 14 Oct.). We see how the whole area between the CF and the sur- face expression of the AF is dominated by sub- sidence: it is noteworthy; that in the same area we find deep sedimentary basins and most of the spread microseismicity observed in the Um- bria-Marche Apennines (Deschamps et al., Vol51,2_3,2008 4-03-2009 10:28 Pagina 398 399 Postseismic deformation and body forces shaping the Apennines and adjacent sedimentary basins in Umbria-Marche 1989). The maximum deformation amounts to about 22 mm, providing a secular rate of 0.03 mm/a if we assume a recurrence time of 700 a, as suggested by Pantosti et al. (1999) after in- terpreting the 1997 sequence as a reoccurrence of the 1279 earthquake. Since the upper crust has brittle properties, modeled vertical relax- ation at the surface is representative of the dis- placement of the whole seismogenic layer. Therefore, we can project the vertical motion onto the AF to obtain average slip rates between 0.02 mm/a and 0.12 mm/a, which are much smaller than the value of 1 mm/a suggested by Chiaraluce et al. (2007) for the northern seg- ment of the AF. The average slip rates we have computed are meant to represent the long term effect of the CF relaxation on the evolution of the AF. Therefore, we are not interested in a detailed description of the relaxation process, which would have to take into account additional fac- tors such as time-variable relaxation and near- ly-free shear stress conditions on the AF. De- pending on the viscosity profile, deformation rates in the first few years after each earthquake can easily reach values of several millimeters per year, as we discuss in a later section. More- over, allowing the AF to slip under specific conditions would surely have an impact on the stress evolution within the upper crust that we have here modeled as a continuous elastic lay- er. However, the use of a more sophisticated model would not have a strong impact on the purpose of this section, i.e. showing that relax- ation of the CF is capable to load the AF, but it cannot represent the primary loading mecha- nism. 3. Slip on the Altotiberina fault Since it appears that the average rates of postseismic relaxation due to earthquakes in the Colfiorito area account for only a small fraction of the AF slip, here we explore the possibility of an opposite interaction, namely an active role of the AF in the evolution of the Umbria- Marche fault zone. For this purpose, we model co- and post-sesimic deformation due to a 50- km-long NE-dipping normal fault representing the AF. Following Boncio and Lavecchia (2000), we distinguish three fault patches along the dip, with the shallowest (0-4 km depth) and deepest (7-14 km depth) patches dipping 30o and the middle patch dipping 13o. We model the coseismic response by means of a half-space model (Feigl and Dupre, 1998) and viscoelastic relaxation as in the previous section. Figure 2 shows the total (elastic plus full re- laxation) vertical deformation due to 70 cm of slip on the AF, consistent with a rate of 1 mm/a (Chiaraluce et al., 2007) during a 700-year seismic cycle (Pantosti et al., 1999). The whole area under the AF is dominated by subsidence, with a maximum due to the elastic component where the AF intersects the free surface and a second maximum due to viscoelastic relax- ation below the CF. The large subsidence be- low the CF reaches a value of about 35 cm, comparable with the average slip accommodat- ed on the largest faults during the recent 1997 sequence (Zollo et al., 1999), suggesting that the AF is capable of actively loading the CF throughout a seismic cycle. Fig. 1. Vertical postseismic relaxation due to slip on a 50-km-long normal fault, dipping 45 degrees to the SW and representing a prolongation of the largest shock of the 1997 sequence (northernmost open box). The three small open boxes represent the surface pro- jections of largest shocks of the 1997 sequence, while the large boxes with bounded by thin lines represent the three segments of the Altotiberina fault. Vol51,2_3,2008 4-03-2009 10:28 Pagina 399 400 R.E.M. Riva and A. AoudiaR.E.M. Riva and A. Aoudia 4. Postseismic relaxation with a lateral heterogeneous lithosphere The first part of our study focused on the av- erage features of postseismic relaxation after earthquakes on the CF. The purpose of the pres- ent section is to discuss an improved postseismic relaxation model consistent with lateral crustal variations present in the Central Apennines. By means of a 2D finite elements approach (Govers and Wortel, 1995), we model a SW-NE section through the Central Apennines extending for 250 km at each side of the CF, where the do- main size has been chosen to minimize the effect of boundary conditions on the area of interest. Lateral variations are introduced according to the results published by Chimera et al. (2003). A large level of heterogeneity characterizes the to- pography of the Moho, which is located at a depth of 30 km under the CF area and at 38 km under the Adriatic Sea, whereas on the left side it presents a minimum of 28 km west of the AF and finally deepens to 33 km towards the Tuscan domain. According to Chimera et al. (2003), a 20-km-thick layer of hot lithospheric material, called mantle wedge, underlies the Moho be- neath the Tuscan magmatic complex and one large part of the Central Apennines. This layer dies out before reaching the Adriatic coast, where a colder thermal profile is expected is al- so observed by heat flow studies (Della Vedova et al., 2000). In our model domain, the two end members are represented by the Tuscan magmat- ic domain on the left, with a thin lithospheric lid leading to high temperature and low viscosities, and the Adriatic domain at the right side, with a thick lithospheric lid and a consequent cold isotherm and high viscosities. In the middle of the domain, under the CF and AF, the presence of the hot mantle wedge induces low viscosities also in the lower crust, even if the area is located above thick lithospheric roots. We impose a lat- erally varying temperature profile consistent with the lithospheric structure and with the ob- served heat flow, and let the viscosity profile be controlled by empirical flow laws for rocks and by laboratory derived rheological parameters (e.g., Ranalli, 1995) . Since we have a large choice of possible rheologies and stress condi- tions, we have introduced an additional con- straint by looking for those combinations of rhe- ology and stress that provide viscosities any- where not lower than 1017 Pa s at the bottom of the lower crust and 1016 Pa s in the mantle wedge. A viscosity of 1017 Pa s at the bottom of the lower crust is consistent with laboratory rhe- ological parameters for soft lower crustal materi- als (e.g. wet diorite in Hansen and Carter, 1982) under realistic conditions in regions of active tec- tonics, with a surface heat flow of 60-70 mW/m2 and a background stress of 10-100 MPa. In the mantle wedge, the presence of low viscosities is inferred from the very low seismic velocities (VS less than 4.2 km/s), suggesting that the wedge might represent a partially molten mantle (Chimera et al., 2003; Aoudia et al., 2007). With our approach, we aim to represent a plausible viscosity structure, where the choice of different parameters would influence the actual values, but not the relative horizontal viscosity varia- tions, which represent our main point of interest. Figure 3, shows the resulting postseismic strain rates one year after the Colfiorito main shock, where the horizontal axis is centered at Fig. 2. Vertical deformation (elastic plus full relax- ation) due to 70 cm of slip on the Altotiberina fault (AF). The large black boxes represent the three seg- ments of the Altotiberina fault, dipping towards the NE with angles of 30o, 13o and 30o (shallow, middle and deep segment respectively). The three small box- es represent the surface projections of the largest shocks of the 1997 sequence. Vol51,2_3,2008 4-03-2009 10:28 Pagina 400 401 Postseismic deformation and body forces shaping the Apennines and adjacent sedimentary basins in Umbria-Marche the CF and the fault is represented by a small white segment. The two top layers represent the crust (upper and lower), the third layer the man- tle wedge, followed by the rest of the litho- sphere and by the asthenosphere. The Moho is indicated by a black line bounded in white. We can clearly see how the highest strain rates are localized in three distinct regions: in the first few kilometers around the fault, and in the bot- tom part of both the lower crust and the mantle wedge. The high strain rates in the upper crust are mainly due to the very large differential stresses induced by the earthquake, while local- ization of strain rates in the other layers is main- ly controlled by the viscosity profile. The com- bined effect of viscosity jumps at the Moho and at the base of the mantle wedge due to changes in rheology, and lateral viscosity variations driv- en by the thermal structure, leads to the localiza- tion of postseismic relaxation in a confined re- gion below the Umbria-Marche Apennines. Vertical postseismic relaxation is shown in fig. 4 for the first 10 years after the earthquake: the choice of representing this particular time- span is due to the fact that large portions of the lower crust and mantle wedge have viscosities lower than 1019 Pa s, meaning that they almost completely relax in the first few years. The up- per panel shows viscoelastic relaxation for the lateral heterogeneous model, while the bottom panel shows the same result for a lateral homo- geneous model with the same structure as the re- gion below the AF. For the case of the heteroge- neous model, the deformation pattern is domi- nated by a large subsidence in a rather narrow region between 10 km and 75 km west of the CF, bounded by a narrow uplifting region around the CF and by a broader uplift further to the west. The large subsidence occurs where large sedimentary basins such as the Foligno basin are mapped. These sedimentary basins are bordered by low-angle normal faults. The effect of lateral heterogeneities is mainly visible in the eastern side: when approaching the colder Adri- atic domain, we observe a drastic reduction in subsidence with respect to the homogeneous model. A further difference between the two models is visible towards the Tuscan domain, where the combined effect of the thinner crust and very warm geotherm leads to an amplifica- tion of the uplifting bulge. According to those Fig. 3. 2D section through the Central Apennines, showing postseismic strain rates at 1 year after the earthquake. The Colfiorito Fault (CF) is represented by the small white segment in the center of the do- main, and the Moho by the thin black line bounded in white. Fig. 4. Postseismic relaxation for the 2D approxi- mation of the CF for both the lateral heterogeneous (top panel) and the homogeneous (bottom panel) model. The four lines depict vertical deformation at four different times after the earthquake: t = 1 a (sol- id), 2 a (long dashed), 5 a (short dashed) and 10 a (dotted). Vol51,2_3,2008 4-03-2009 10:28 Pagina 401 402 R.E.M. Riva and A. AoudiaR.E.M. Riva and A. Aoudia results, we think that a lateral heterogeneous lithospheric structure has the effect of localizing postseismic relaxation due to earthquakes in the Colfiorito area in a rather narrow region below the AF and the related sedimentary basins. We have verified that our conclusions are not dependent on the adopted 2D approxima- tion, because there is no substantial qualitative difference between relaxation of a 50-km-long fault, representative of multiple ruptures in the Colfiorito area between Camerino and Norcia, and that of an infinitely long fault. The key pa- rameter controlling the wavelength of the re- laxation signal, in fact, is represented by the thickness of the upper elastic layer, which is here about 15 km and therefore much shorter than the extended CF. By comparing relaxation in a vertically layered medium for finite faults of different length, we have further verified that the quantitative effect of the 2D approxi- mation results in the amplification of the mag- nitudes of relaxation by about a factor of two, which is again not relevant for the purpose of this study. 5. Discussion and conclusions Here elaborate further on the possibility of an interaction between the low-angle NE-dip- ping Altotiberina fault (AF) and the SW nor- mal faulting structures responsible for the 1997 Colfiorito sequence. With loose con- straints on the magnitude and location of his- torical earthquakes, as well as on their recur- rence time, we have modeled a possible sce- nario of postseismic relaxation which appears to be consistent that the formation of major sedimentary basins at the western side of the Umbria-Marche faulting zone. The location of the postseismic relaxation signal favours the loading of the AF, even if the slip rates that we have obtained are rather small: however, un- certainties are very large because those values are directly dependent on the average earth- quake magnitude and recurrence time. On the other hand, after modeling the effect of slip on the AF, we have shown that the AF alone could be responsible for providing most of the load necessary to rupture the CF. When we compare average postseismic re- laxation rates due to the CF with slip rates on the AF, the latter appears to provide the domi- nating signal. However, for a realistic viscosity structure where the lower crust and the mantle wedge below the Umbria-Marche Apennines are characterized by low viscosities (in the range 1016-1018 Pa s), about half of the total postseismic relaxation is likely to take place in the first 10 years after the earthquake. This means that, at least locally, the postseismic loading on the AF could easily exceed the aver- age slip rates. In addition, the impact of stress relaxation on the evolution of the AF can be in- creased by the presence of large lateral hetero- geneities in the lithospheric structure, which have the effect of concentrating postseismic re- laxation below the AF. In conclusion, a possible scenario for the in- teraction between the AF and the CF allows aseismic slip on the AF to effectively load the CF, while postseismic relaxation after each rup- ture event on the CF plays the role of a feed- back process, reloading sections of the AF on short time-scales and potentially affecting the evolution of the AF itself. From a regional perspective, the strain pat- tern retrieved from the 2D finite element mod- eling of the postseismic deformation mimics the modeled stress pattern that accounts only for buoyancy forces as reported by Aoudia et al. (2007) along the same section and using the same geometry of the crust upper mantle. In fact the realistic pattern of localized vertical deformation reported in fig. 4 compares fairly well with its equivalent part reported in Aoudia et al. (2007, fig. 3). This is evidence that the two different length and time scale processes, namely, postseismic deformation and buoyancy, are complementary in shaping the Apennines and adjacent sedimentary basins. The high resolution deformation pat- terns modeled in this study as well as the one reported in Aoudia et al. (2003) satisfying the ongoing observed and complex geodynamics across North Central Italy with emphasis to the Umbria-Marche region can hardly be re- produced by a model accounting only for ex- ternal forces such as a rotating or subducting or retreating Adria. Vol51,2_3,2008 4-03-2009 10:28 Pagina 402 403 Postseismic deformation and body forces shaping the Apennines and adjacent sedimentary basins in Umbria-Marche Acknowledgements This work was done in the framework of the SISMA- Italian Space Agency - Project. We thank Rob Govers for the use of the FEM soft- ware GTECTON and two anonymous review- ers for their constructive comments. REFERENCES AGOSTINETTI, N.P., F. P. LUCENTE, G. SELVAGGI and M. DI BONA (2002): Crustal structure and Moho geometry be- neath the Northern Apennines (Italy), Geophys. Res. Lett., 29. AMATO, A., R. AZZARA, C. CHIARABBA, G.B. CIMINI, M. COCCO, M. DI BONA, L. MARGHERITI, S. MAZZA, F. MELE, G. SELVAGGI, A. BASILI, E. BOSCHI, F. COUR- BOULEX, A. DESCHAMPS, S. GAFFET, G. BITTARELLI, L. CHIARALUCE, D. PICCININI and M. RIPEPE (1998): The 1997 Umbria-Marche, Italy, earthquake sequence: a first look at the main shocks and aftershoks, Geophys. Res. Lett., 25 (15), 2861-2864. AOUDIA, A., A.T. ISMAIL-ZADEH and F. ROMANELLI (2007): Buoyancy-Driven Deformation and Contemporary Tectonic Stress in the Lithosphere beneath Central Italy, Terra Nova, 19 (6), 490-495, doi: 10.1111/j.1365-3121.2007.00776.x AOUDIA, A., A. BORGHI, R. RIVA, R. BARZAGHI, B.A.C. AM- BROSIUS, R. SABADINI, L.L.A. VERMEERSEN and G.F. PANZA (2003): Postseismic deformation following the 1997 Umbria-Marche (Italy) moderate normal faulting earthquakes, Geophys. Res. Lett., 70, 1390-1393. BARBA, S., and R. BASILI (2000): Analysis of seismological and geological observations for moderate-size earth- quakes: the Colfiorito Fault System (Central Apen- nines, Italy), Geophys. J. Int., 141, 241-252. BONCIO, P., and G. LAVECCHIA (2000): A geological model for the Colfiorito earthquakes (September-October 1997, central Italy), Journal of Seismology, 4, 345-356. BOSCHI, E., E. GUIDOBONI, G. FERRARI, D. MARIOTTI, G. VALENSISE and P. GASPERINI (eds) (2000): Catalogue of Strong Italian Earthquakes, Ann. Geofis., 43 (4), 268. CHIARALUCE, L., C. CHIARABBA, C. COLLETTINI, D. PICCINI- NI and M. COCCO (2007): Architecture and mechanics of an active low-angle normal fault: Alto Tiberina Fault, northern Apennines, Italy, Journal of Geophysical Re- search, 112, B10310, doi: 10.1029/2007JB005015, CHIMERA, G., A. AOUDIA, A. SARAO and G.F. PANZA (2003): Active tectonics in Central Italy: constraints from surface wave tomography and source moment tensor inversion, Phys. Earth Planet. Inter., 138, 241- 262. DELLA VEDOVA, B., S. BELLANI, G. PELLIS and P. SQUARCI (2000): Deep Temperatures and Surface Heat-Flow Distribution, in Anatomy of an orogen: the Apennines and adjacent Mediterranean basins, edited by G.B. VAI and L.P. MARTINI (Kluwer Academic Publishers, Dor- drecht). DESCHAMPS, A., R. SCARPA and G. SELVAGGI (1989): Anali- si sismologica del settore settentrionale dell’Appenni- no umbro-marchigiano, Atti GNGTS 8th Conference, 1, 9-15. FEIGL, K.L. and E. DUPRE (1998): RNGCHN: a program to calculate displacement components from dislocations in an elastic half-space with applications for modeling geodetic measurements of crustal deformation, revised for Computers and Geosciences. GALADINI, F., P. GALLI, I. LESCHIUTTA, G. MONACHESI and M. STUCCHI (1999): Active tectonics and seismicity in the area of the 1997 earthquake sequence in central Italy: A short review, Journal of Seismology, 3 (2), 167-175. GOVERS, R. and M.J.R. WORTEL (1995): Extension of stable continental lithosphere and the initiation of lithosphere scale faults, Tectonics, 14 (4), 1041-1055. HANSEN, F.D. and N.L. CARTER, (1982): EOS Trans. AGU, 63, 437. MEGHRAOUI, M., V. BOSI and T. CAMELBEECK (1999): Fault fragment control in the 1997 Umbria-Marche, central Italy, earthquake sequence, Geophys. Res. Lett., 26 (8), 1069-1072. MELE, G. and E. SANDVOL (2003): Deep crustal roots be- neath the northern Apennines inferred from teleseismic receiver functions, Earth Planet. Sci. Lett., 211, 69-78. OLIVIERI, M. and G. EKSTROM (1999): Rupture depths and source processes of the 1997-1998 earthquake sequence in Central Italy, Bull. Seism. Soc. Am., 89, 305-310. PANTOSTI, D., P.M. DE MARTINI, P. GALLI, F. GALADINI, P. MESSINA, M. MORO and A. SPOSATO, (1999): Studi pa- leosismologici attraverso la rottura superficiale prodot- ta dal terremoto del 14 ottobre 1997 (Umbria-Marche) Atti 18° Convegno Annuale del GNGTS, sessione 10 “Forti terremoti dell’area mediterranea”. RANALLI, G., (1995): Rheology of the earth (second edition, Kluwer Academic Publishers, Dordrecht), pp. 432. RIVA, R., A. AOUDIA, L.L.A. VERMEERSEN, R. SABADINI and G.F. PANZA (2000): Crustal versus asthenospheric re- laxation and post-seismic deformation for shallow nor- mal faulting earthquakes: the Umbria-Marche (central Italy) case, Geophys. J. Int., 141, F7-F11. RIVA, R.E.M., A. BORGHI, A. AOUDIA, R. BARZAGHI, R. SA- BADINI and G. F. PANZA (2007): Viscoelastic relaxation and long-lasting after-slip following the 1997 Umbria- Marche (Central Italy) earthquakes, Geophys. J. Int., 169, 534-546, doi: 10.1111/j.1365-246X.2007.03315.x SABADINI, R. and L.L.A. VERMEERSEN (1997): Influence of lithospheric adn mantle layering on global post-seismic deformation, Geophys. Res. Lett., 24, 2075-2078. SALVI, S., S. STRAMONDO, M. COCCO, M. TESAURO, I. HUN- STAD, M. ANZIDEI, P. BRIOLE, P. BALDI, E. SANSOSTI, G. FORNARO, R. LANARI, F. DOUMAZ, A. PESCI and A. GALVANI (2000): Modeling coseismic displacements resulting from SAR interferometry and GPS measure- ments during the 1997 Umbria-Marche seismic se- quence, Journal of Seismology, 4, 479-499. ZOLLO, A., MARCUCCI, S., MILANA, G. and P. CAPUANO (1999): The 1997 Umbria-Marche (central Italy) earth- quake sequence: Insights on the mainshock ruptures from near source strong motion records, Geophys. J. Int., 26 (20), 3165-3168. Vol51,2_3,2008 4-03-2009 10:28 Pagina 403