PAPATHANASSIOU_corretto_Layout 6 735 Liquefaction phenomena along the paleo-Reno River caused by the May 20, 2012, Emilia (northern Italy) earthquake Giorgos Papathanassiou1, Riccardo Caputo2,*, Dimitra Rapti-Caputo2 1 Aristotle University of Thessaloniki, Department of Geology, Thessaloniki, Greece 2 Università di Ferrara, Dipartimento di Fisica e Scienze della Terra, Ferrara, Italy ANNALS OF GEOPHYSICS, 55, 4, 2012; doi: 10.4401/ag-6147 1. Introduction On the May 20, 2012 (04:03:52 local time; 02:03:52 UTC), a moderate earthquake (ML 5.9) [Scognamiglio et al. 2012, this volume] with a focal mechanism showing E-W-trending, S-dipping, reverse-faulting occurred in the eastern sector of the alluvial plain of the Po River, close to the border between the Regions of Emilia-Romagna and Lombardia (northern Italy). The tectonic structure is completely blind, but it was well known from a dense grid of seismic profiles for hydro- carbon explorations [e.g., Pieri and Groppi 1981, Toscani et al. 2009]. The earthquake triggered extensive liquefaction-in- duced ground effects at the surface, and caused severe struc- tural damage to nonreinforced masonry and precast industrial buildings within the broader epicentral area. The hypocenter was at 44.89 ˚N, 11.23 ˚E, at a depth of 6.3 km [Scognamiglio et al. 2012], while the maximum acceleration was recorded in Mirandola, with peak ground acceleration 310 cm/s2 and 264 cm/s2 along the vertical and horizontal components, re- spectively [Bozzoni et al. 2012, this volume]. In this report, we focus on a zone including the Sant'A- gostino, San Carlo and Mirabello villages (west Ferrara Province), which were built along an abandoned reach of the Reno River and where liquefaction phenomena were partic- ularly diffuse, with very intense local effects (Figure 1). The map in Figure 1 is based on original data collected in the days following the mainshock (and before the event of May 29, 2012) and locally integrated with data published by ISPRA [2012]. Similar maps have also been presented by other geo- logical teams [e.g., Decanini et al. 2012, EMERGEO Working Group 2012], including of the western sector of the epicen- tral area and the manifestations that occurred during the sec- ond major earthquake. However, for the aims of this report that focuses on the Sant'Agostino, San Carlo and Mirabello case study, they will be not discussed. Based on stratigraphic reconstructions from boreholes and penetrometric tests, and on new direct measurements within numerous water wells, a preliminary conceptual hy- drogeological model for the first ca. 30 m is also proposed for the investigated area. This hydrogeological model pro- vides essential information for the liquefaction susceptibility discussed in this report. Grain-size analyses of ejected mate- rial collected a few days after the event in several locations were also carried out, and are here discussed in relation to the classical susceptibility criteria proposed by Tsuchida [1971]. In addition, a preliminary evaluation of the macro- seismic intensity is performed based on the Environmental Seismic Intensity (ESI) scale [Michetti et al. 2007]. 2. Earthquake-induced ground effects A few hours after the event, a post-earthquake field sur- vey was organized to report the coseismic effects and to com- pile a map of their spatial distribution (Figure 1). As this survey began a few hours after the mainshock, this provided the op- portunity to observe, measure, photographically record, and map the sites where the sandy material was ejected, before the residents and local authorities started removing the liquefied sandy materials from the courtyards and roads. During this post-event investigation, more than 150 lo- cations were surveyed, with the focus on a zone 8 km long and 0.5 km wide, with the parameters that characterize the surface liquefaction manifestations were measured and recorded carefully. In particular, the length and width of vent fractures, and the diameter of the sand craters were meas- ured, and samples were collected from ejected materials at representative sites. The length of the liquefaction-induced surface ruptures varied, as expected, from a few centimetres to more than 10 m, and locally they were reported as up to 50 m [Lai et al. 2012]. The width variation of the cracks ranged from a few centimetres to 30 cm. The ejection of sandy material occurred through conduits and fractures in farmed fields, although the most spectacular effects affected pavements and large-diameter water wells within the urban- Article history Received July 23, 2012; accepted September 18, 2012. Subject classification: Seismic hazard assessment, Po Plain, Morphogenic earthquakes. 2012 EMILIA EARTHQUAKES ized areas (Figure 2). As a consequence, the number of sand craters and volcanoes were significantly lower than the vent fractures, but above all, their location and geometry was strongly controlled by the artificial discontinuities within the paving and other infrastructures. The diameter in these cases ranged from 4 cm up to 40 cm. Lateral spreading due to liquefaction and the conse- quent upward injection of sandy material is also well docu- mented by electrical resistivity tomographies [Abu-Zeid et al. 2012, this volume], and this phenomenon induced struc- tural damage within the urbanized areas, including broken pipelines and differential settling, causing minor damage to buildings locally. In addition, short interviews with the villagers were car- ried out relating to the ejection of fluids and sandy material. The questions posed focused on the timing of the ejection phenomena, relative to the occurrence of the mainshock, the duration of the outflow, and the height of the ejection with respect to the ground surface. Important outcomes of these interviews were: (i) most of the sandy material within the three urban areas that was focused on (Mirabello, San Carlo, Sant'Agostino) was ejected through large diameter water wells; (ii) the ejection of sandy material started almost at the same time as the mainshock; and (iii) in the first seconds/min- utes of the ejection process, clear water came to the surface first, only later mixed with sand. The duration of the phe- nomenon commonly lasted up to a few minutes, although at one site (44˚ 48' 03" N, 11˚ 24' 34" E) an eyewitness stated that the outpouring of water after the major sand ejection phe- nomenon lasted for several hours. Differences in the amounts of material and timing are likely to be due to local geological, hydrogeological and anthropogenic conditions. For the height of the ejection, the most characteristic re- port is that from a well at the center of San Carlo (44 ̊48' 15.39"N, 11˚ 24' 35.14"E), where the liquefied material was ejected up to 2 m above the ground surface, and with such intensity as to damage the well. Taking into account that the depth to the water at the site was ca. 5 m (see next section), it can be as- sumed that the total height of the ascent was more than 7 m. Considering the average density of the outflowing material (a mixture of water and sand), this implies a fluid pressure of probably more than one bar (0.1 MPa). Similar ascent and ejection of sand through wells was triggered by the Friuli 1976, northeastern Italy, earthquake [Sirovich 1996]. To assess the macroseismic intensity in the area studied, we applied the ESI scale [Michetti et al. 2007], as the gener- PAPATHANASSIOU ET AL. 736 Figure 1. Map of the spatial distribution of the liquefaction-induced surface evidence within the west Ferrara Province from original data (yellow circles) and from the ISPRA [2012] report (white circles). Red star, May 20, 2012, earthquake epicenter. Base map from Google Inc. (2012). 737 ated geological effects were widespread and extensive. Thus, taking into account the length and width of the surface rup- tures, and the dimensions of vent fractures and sand boils of liquefaction phenomena, we conclude that the macroseismic intensity should be evaluated as VIII in the area of San Carlo. In particular, according to the ESI definition of intensity de- grees [Michetti et al. 2007], a value of VII-VIII is appropriate in areas where fractures up to 50 cm wide and up to hundreds of meters long are commonly observed in loose alluvial de- posits. Decametric cracks are common in paved roads, and liquefaction can be frequent in the epicentral area. Localized lateral spreading and settling (up to ca. 30 cm) can result in the most typical effects of sand boiling up to ca. 1 m in diameter. 3. Hydrogeological model The conceptual shallow of the hydrogeological model of the broader Ferrara Province is characterized by the pres- ence of two aquifer systems. The first is superficial and semi- confined or locally unconfined (A0 according to Molinari et al. [2007]), and the second is deeper and always confined (A1; Regione Emilia-Romagna and ENI-AGIP [1998]). Within the investigated area, and based on published [Emilia-Romagna Region Open file] and unpublished strati- graphic data, the A0 aquifer system consists of silty-sandy and locally peat and sandy lenses, with variable grain-sized distributions containing a semi-confined and locally uncon- fined aquifer. According to interviews with the villagers, the piezometric level in many wells showed limited variations even after a long exploitation period, thus indicating that the aquifer bodies belonging to A0 are commonly hydraulically inter-connected. Only minor lenses were probably isolated, as suggested by their complete draining and a longer recharge time following a few hours of pumping. The exploitation of this aquifer system is performed using boreholes of about 1 m in diameter, which were often excavated many years ago, commonly to a depth of 7 m to 13 m. This traditional water-well excavation and their architec- ture was maintained by local people up to the last decades. Before the aqueduct network reached these areas, these large diameter wells were the only means to provide the water for every domestic use. The spatial distribution and the geometric characteris- tics of these highly heterogeneous aquifer units were largely influenced by the latest Quaternary geomorphological and sedimentary evolution of the territory, and especially by the presence and path of active and abandoned channels and out- flow cones [e.g. Rapti-Caputo 2000]. Within the investigated area, the thickness of the aquifer varies from between a few meters to 6 m to 8 m, in agreement with the paleo-reach of the Reno River (Figure 3). The hydraulic conductivity of the A0 aquifer system was preliminarily estimated based on the available stratigraphic data and the grain-size analyses we performed from collected samples of ejected material (see next section). According to Hazen [1911], this parameter varies from between 10–5 and 10–4 m/s, thus indicating highly permeable layers. The prin- 2012 EMILIA EVENT LIQUEFACTION PHENOMENA Figure 2. Example of sand ejection from a large-diameter water well at the base of the levee associated with an abandoned reach of the Reno River (San Carlo-Sant'Agostino road). Figure 3. Topographic profile (black line) drawn from a lidar-based DEM, and piezometric level of the unconfined aquifer system (A0), as reconstructed from water-well measurements a few weeks after the seismic event. cipal recharge of the aquifer was primarily due to vertical in- filtration from the precipitation, and secondarily to recharges via lateral infiltration from the active hydrographic network, as documented in nearby areas with a comparable hydroge- ological–geomorphological setting [e.g. Rapti-Caputo 2000]. Indeed, the local alimentation is shown by high values of tri- tium (3-6 TU) that show that the mean residence time of the water in the underground reservoir is only 2 years to 5 years [Rapti-Caputo 2008, Rapti-Caputo and Martinelli 2009]. The piezometric level varies by about 1.5 m to 2.5 m during the hydrological year, due to the variability of both the seasonal precipitation and the hydrometric level within the hydro- graphic network [Rapti-Caputo 2008, ARPA GdL 2010]. Based on a campaign of measurements specifically car- ried out for 45 wells during the weeks following the earth- quake, the water depths were about 1.4 m to 3.0 m in the interfluvial plain, showing deeper values (4 m to 5 m from the surface) in the boreholes located along the levees of the paleo-Reno River. This hydrogeological setting is shown in Figure 3, where the influence of the paleo-channel on the di- rection of the groundwater flow is evident for generating a marked groundwater divide. A low permeability layer that consists prevalently of clay and silt locally alternating with peat lenses separates the un- confined (A0) from the underlying confined aquifer (A1) lo- cated about 16 m to 18 m below the surface [Rapti-Caputo 2008, Regione Emilia-Romagna GdL-liquefazione 2012]. The latter aquifer body consists of high permeability sediments (mainly sand), and has a thickness of between 4 m and 9 m. 4. Liquefaction-related effects To evaluate the susceptibility to liquefaction of the soil deposits within the investigated area, we followed the pro- cedure recommended by the Manual for Zonation on Seis- mic Geotechnical Hazards, prepared by The Technical Committee for Earthquake Geotechnical Engineering, TC4, of the International Society for Soil Mechanics and Geot- echnical Engineering, and published by The Japanese Geot- echnical Society [TC4, 1999]. Initially, the relationship between the epicentral distance of liquefied sites and earth- quake magnitude was examined. The epicentral distances of the villages of Sant'Agostino, San Carlo and Mirabello are 16 km, 17 km and 19 km, respectively (Figure 4). In addition, the epicentral distances of the liquefied area of Scortichino- Burana, documented by ISPRA [2012], varied from 3 km to 11 km, while the zone situated to the north of Mirabello is about 17 km from the epicenter. According to several em- pirical relationships between magnitude and epicentral dis- tance (Figure 4), these areas should be clearly considered as susceptible to liquefaction for an earthquake of ML 5.9, as they lie above most of the proposed curves, with the only possible exceptions of some sites in the Mirabello area. Differences between the proposed empirical relation- ships (and hence the curves represented in Figure 4) are due to the different datasets used in these studies. For example, Ambraseys [1988] used a worldwide set of liquefaction data, while Papadopoulos and Lefkopoulos [1993] updated the dataset collected by Ambraseys with 30 new cases from Greek earthquakes and with three other ones. Wakamatsu PAPATHANASSIOU ET AL. 738 Figure 4. Distribution of the liquefied sites based on their epicentral distances and comparison with the triggering curves proposed by several studies for the assessment of liquefaction susceptibility. 739 [1993] extended the work of Kuribayashi and Tatsuoka [1975] by adding 46 new data points from Japanese earthquakes, while Galli [2000] and Aydan et al. [2000] re-evaluated the seismic parameters of Italian and Turkish earthquakes, re- spectively, using data from after 1900. The curve suggested by Papathanassiou et al. [2005] was constructed based on lique- faction phenomena that were reported after 1900 in Greece. The sites where liquefaction-induced ground effects are documented generally consist of Holocene fluvial sediments and locally man-made deposits. The villages of Sant'A- gostino, San Carlo and Mirabello were founded upon the lev- ees of a former reach of the Reno River [Graziani et al. 1992, Castiglioni et al. 1999, Caputo et al. 2012 (this volume)]. This was for the obvious reason that this site was (and still is) 3 m to5 m higher with respect to the surrounding alluvial plain, which at that time (XV-XVIII century) was repeatedly flooded. Moreover, the presence of a locally sandy substra- tum associated with the crevasse deposits ascertained better stability for house foundations, in contrast to other sectors of the plain where the muddy material was not suitable for con- struction. Notwithstanding these historical and social rea- sons, these morphological, stratigraphic and hydrogeological conditions were characterized by: (i) topographic gradients; (ii) the presence of cohesionless deposits; and (iii) a depth of the water table of less than 1 m locally. This represented the perfect preparatory factors for liquefaction, just waiting for sufficient shaking to trigger the physical process. Eventually, on May 20, 2012, sufficient shaking occurred, and this caused locally important and widespread liquefaction phenomena. At the surface, this produced small-scale structural damage, mainly due to lateral spreading (44˚ 48' 22" N, 11˚ 24' 46" E) and differential settling (44˚ 49' 29"N - 11˚ 27' 26"E). Similar coseismic effects within old and/or former river channel zones are well documented for at least three recent earthquakes: Luzon Philippines in 1990, Niigaten Chuetsu- oki, Japan, in 2007 and Darfield, New Zealand, in 2010 [Wotherspoon et al. 2012]. This confirms that this type of morphological setting is favorable for earthquake-induced surface ruptures and liquefaction phenomena. The above conclusion is in agreement with the characterization of this type of deposit as highly susceptible to liquefaction during strong seismic shaking, based on the geological criteria that were suggested by Youd and Perkins [1978]. The samples collected during the field survey were an- alyzed in the Laboratory of Engineering Geology and Hy- drogeology of the Department of Geology at the Aristotle University of Thessaloniki (Greece), to obtain the grain-size distribution, the liquid limit, and the plasticity index [e.g., Kramer 1996] (Figure 5). Based on the Unified Soil Classifi- cation System (USCS), samples are classified (see Table 1) as silty sand (SM) with non-plastic fine particles, and fine sand (SP-SM). The fines content of the samples analyzed varied from 6% to 12% for the SP-SM, and from 20% to 34% for the SM. The results suggest that the majority of the ejected sam- ples had probably come from a unique sedimentary layer. The probable source level was identified using cone pene- tration tests and cores, which was carried out within the in- vestigated area in a similar morphological, hydrogeological 2012 EMILIA EVENT LIQUEFACTION PHENOMENA Figure 5. Comparison between the grain-size distribution curves of the collected samples and the boundary curves for liquefiable (dashed lines) and probable liquefiable (dotted lines) soils, as suggested by Tsuchida [1971]. Similar critical curves were proposed by EN 1998-1 [2004] and NTC [2008]. and stratigraphic setting [Regione Emilia-Romagna GdL-liq- uefazione 2012]. This source unit consists of a medium-to- coarse sand layer at a depth of 7 m to 8 m from the surface, which varies in thickness between 1.5 m and 4 m, and which probably recorded large crevasse splays and overbank flood episodes. Due to the stratigraphic setting and the hydrogeo- logical framework, this layer represents a semi-confined aquifer that is characterized by a high degree of liquefaction susceptibility. The presence of a non-liquefied silty-clayey cap layer and the occurrence of some erosion during the upward flow along the dikes and conduits might have caused an over- estimation of the fines content of the ejected material. Table 1 lists the results of the collected samples that were provided by the laboratory testing, including the fines content (%), the classification (USCS), and their geographic location. The results of the grain-size analyses of the samples are also compared with the curves proposed by Tsuchida [1971] for soils at uniform grading, as shown in Figure 5. The grain-size distribution curves obtained fall within the suggested range for liquefaction conditions, therefore char- acterizing these deposits as susceptible to liquefaction in terms of their grain size. 5. Concluding remarks The May 20, 2012, Emilia earthquake triggered exten- sive liquefaction-induced ground effects at the surface and caused structural damage to unreinforced masonry and pre- cast industrial buildings within the broader epicentral area. The post-earthquake field survey organized a few hours after the event focused on a zone 8 km long and 0.5 km wide, and it was devoted to the reporting of the co-seismic effects and to compile a map of their spatial distribution. The lengths of the liquefaction-induced surface ruptures varied from a few centimeters to more than 10 m, while the width variation of the cracks was not so great, ranging from a few centimeters to 30 cm. The ejection of sandy material occurred through conduits and fractures in farmed fields, although the most spectacular effects affected roads, pavements and water wells within the urbanized areas. In addition, short interviews with the villagers were car- ried out concerning the ejection of fluids and sandy mate- rial. The outcomes of these interviews are that most of the sandy material within the urban areas of Mirabello, San Carlo and Sant'Agostino was ejected through water wells of large diameter, which are typical of the region, and exploited the first unconfined aquifer. The phenomenon started almost at the same time as the mainshock, and in the first sec- onds/minutes of the ejection process clear water came to the surface first, only afterwards mixed with sand. Furthermore, samples of the ejected material were col- lected, to obtain grain-size distribution, liquid limits and plas- ticity indices. Based on the USCS, the materials were classified as SM with non-plastic fine particles, and SP-SM. The fines content of these deposits varied from 6% to 12% for the SP-SM and from 20% to 34% for the SM. The results suggest that the majority of the ejected samples had proba- bly come from the same sedimentary layer. A preliminary outcome of these field surveys is that the studied area has suffered large-scale liquefaction-induced manifestations that should be considered as secondary co- seismic effects (i.e., not directly associated with the upwards propagation of the fault rupture). Although liquefaction phe- PAPATHANASSIOU ET AL. 740 a/a Longitude Latitude USCS classification Fines content (%) S1 44˚ 48' 16.85" 11˚ 25' 08.30" SM 31 S2 44˚ 48' 17.19" 11˚ 25' 07.42" SM 29 S3 44˚ 48' 17.73" 11˚ 24' 37.49" SP-SM 7 S4 44˚ 48' 12.68" 11˚ 24' 25.74" SP-SM 6 S5 44˚ 48' 10.61" 11˚ 24' 30.99" SM 18 S6 44˚ 48' 08.99" 11˚ 24' 33.14" SM 17 S7 44˚ 48' 25.16" 11˚ 24' 45.67" SM 32 S8 44˚ 48' 20.28" 11˚ 24' 57.20" SM 34 S9 44˚ 47' 58.54" 11˚ 24' 09.06" SM 23 S10 44˚ 47' 55.64" 11˚ 24' 04.41" SP-SM 8 S11 44˚ 47' 50.16" 11˚ 24' 03.68" SM 18 S12 44˚ 48' 35.81" 11˚ 25' 47.42" SM 30 S13 44˚ 49' 37.12" 11˚ 27' 41.68" SM 20 S14 44˚ 49' 29.57" 11˚ 27' 27.13" SM 30 Table 1. Grain size characteristics of the collected material. 741 nomena were diffuse in the epicentral area and in general the liquefaction potential might be characterized as 'high' for re- cent alluvial plains, only small-scale liquefaction-induced structural damage was documented [Galli et al. 2012, Ter- tulliani et al. 2012, this volume]. Taking into account that the geological (stratigraphic, morphological, hydrogeological seismotectonic) setting of the investigated area is comparable with those of many other sectors of the broader Po Plain, as well as worldwide, and due to the abundance of liquefaction-related features pro- duced here, the investigated case study represents a valuable example for future investigations devoted to better evalua- tion of liquefaction hazard. Acknowledgements. The authors thank the Servizio Geologico, Sis- mico e dei Suoli, Regione Emilia-Romagna for providing the lidar-based DEMs and preliminary access to the sedimentary cores drilled in the in- vestigated area. References Abu Zeid, N., S. Bignardi, R. Caputo, G. Santarato and M. Stefani (2012). Electrical resistivity tomography investi- gation of coseismic liquefaction and fracturing at San Carlo, Ferrara Province, Italy, Annals of Geophysics, 55 (4); doi:10.4401/ag-6149. Ambraseys, N. (1988). Engineering seismology, J. Earthq. Eng. Struct. Dyn., 17, 1-105. ARPA GdL (2010). Supporto tecnico alla Regione Emilia- Romagna per la ricostruzione del quadro conoscitivo relativo ai principali acquiferi freatici di pianura; http:// www.arpa.emr.it/cms3/documenti/_cerca_doc/acqua/ piacenza/sotterranee/freatico_pc2010_allegati.pdf Aydan, O., R. Ulusay, H. Kumsar and E. Tuncay (2000). Site investigation and engineering evaluation of the Duzce- Bolu earthquake of November 12, 1999, Turkish Earth- quake Foundation, Instanbul, Report No. TDV/DR 09-51, 307 pp. Bozzoni, F., C.G. Lai and L. Scandella (2012): Preliminary re- sults of ground-motion characteristics, Annals of Geo- physics, 55 (4); doi:10.4401/ag-6121. Caputo, R., K. Iordanidou, L. Minarelli, G. Papathanassiou, M.E. Poli, D. Rapti-Caputo, S. Sboras, M. Stefani and A. Zanferrari (2012). Geological evidence of pre-2012 seis- mic events, Emilia-Romagna, Italy, Annals of Geophysics, 55 (4); doi:10.4401/ag-6148. Castiglioni, G., A. Biancotti, M. Bondesan, G.C. Cortemiglia, C. Elmi, V. Favero, G. Gasperi, G. Marchetti, G. Orombelli, G.B. Pellegrini and C. Tellini (1999). Geomorphological map of the Po Plain, Italy, at a scale of 1:250.000, Earth Surf. Proc. Land., 24, 1115-1120. Decanini, L.D., D. Liberatore, L. Liberatore and L. Sor- rentino (2012). Preliminary Report on the 2012, May 20, Emilia Earthquake, vers. 1, http://www.eqclearing- house.org/2012-05-20-italy-it/ EMERGEO working group (2012). Rilievi geologici di ter- reno effettuati nell'area epicentrale della sequenza sismica della Pianura Padana Emiliana del maggio-giugno 2012, Rapporto 04.06.2012; http://emergeo.ingv.it/images/ emergeo/Emergeo04062012.pdf EN 1998-1 (2004). Eurocode 8 – Design of Structures for earthquake resistance – Part 1: General rules, seismic ac- tions and rules for buildings, EN 1998-1, CEN, Brussels. Galli, P. (2000). New empirical relationships between magni- tude and distance for liquefaction, Tectonophysics, 324, 169-187. Galli, P., S. Castenetto and E. Peronace (2012). The MCS macroseismic survey of the Emilia 2012 earthquakes, An- nals of Geophysics, 55 (4); doi:10.4401/ag-6163. Graziani, S., R. Ferri and M. Bondesan (1992). Aspetti geo- morfologici e problemi paleogeografici della zona fra Bondeno, Finale Emilia e Mirabello nel quadro degli an- tichi domini idrografici del Secchia, del Panaro e del Reno, In: S. Gelichi (ed.), Un mito e un territorio: Ansalaregina e l'Alto Ferrarese nel Medioevo, All'Insegna del Giglio, Firenze, 13-44. Hazen, A. (1911). Discussion: dams on sand foundations, T. Am. Soc. Civ. Eng., 73, 199 pp. ISPRA (2012). Geological effects induced by the seismic se- quence started on May 20, 2012, in Emilia (Mw=5.9) – Preliminary report, http://www.isprambiente.gov.it/files/ progetti/inqua/inqua-primo-rapporto-emilia-2012.pdf Kramer, S.L. (1996). Geotechnical Earthquake Engineering, Prentice Hall, New Jersey, 348-355. Kuribayashi, E., and F. Tatsuoka (1975). Brief review of liq- uefaction during earthquakes in Japan, Soils Found., 15, 81-92. Lai, C.G., F. Bozzoni, M.-D. Mangriotis, M. Martinelli, L. Scandella and F. Dacarro (2012). Geotechnical aspects of the May 20, 2012 M5.9 Emilia earthquake, Italy, version 1, July 13, 2012; http://www.eqclearinghouse.org/2012-05- 20-italy/files/2012/07/Report_GEO_Emilia_FINAL.pdf Michetti, A.M., E. Esposito, L. Guerrieri, S. Porfido, L. Serva, R. Tatevossian, E. Vittori, F. Audemard, T. Azuma, J. Clague, V. Comerci, A. Gurpinar, J. McCalpin, B. Mo- hammadioun, N.A. Mörner, Y. Ota and E. Roghozin (2007). Intensity Scale ESI 2007, In: L. Guerrieri and E. Vittori (eds.), Memorie Descrittive della Carta Geologica d'Italia, Roma, 74, 53 pp. Molinari, F.C, G. Boldrini, P. Severi, G. Duroni, D. Rapti-Ca- puto and G. Martinelli (2007). Risorse idriche sotterranee della Provincia di Ferrara. Regione Emilia Romagna, DB MAP Ed., Firenze, 61 pp. NTC (2008). Norme tecniche per le Costruzioni, D.M. 14.1.2008 (Italian Building Code, in Italian). Papadopoulos, G.A., and G. Lefkopoulos (1993). Magnitude- distance relations for liquefactionin soil from earth- quakes, B. Seismol. Soc. Am., 83 (3), 925-938. 2012 EMILIA EVENT LIQUEFACTION PHENOMENA Papathanassiou, G., S. Pavlides, B. Christaras and K. Pitilakis (2005). Liquefaction case histories and empirical relations of earthquake magnitude versus distance from the broader Aegean Region, J Geodyn., 40, 257-278. Pieri, M., and G. Groppi (1981). Subsurface geological struc- ture of the Po Plain, Italy, Consiglio Nazionale delle Ricerche, Progetto Finalizzato Geodinamica, Sottoprog- etto Modello Strutturale, 414, Roma, 13 pp. Rapti-Caputo, D. (2000). Risorse idriche sotterranee a est di Ferrara: indagini sul comportamento idrogeologico e idrochimico. Proposte per una gestione ottimale, Ph.D. thesis, University of Ferrara, 215 pp. Rapti-Caputo, D. (2008). Hydrodynamic and hydrochemical monitoring of the Sant'Agostino landfill, University of Ferrara, Internal report, 180 pp. (unpublished). Rapti-Caputo, D., and G. Martinelli (2009). Geochemical and isotopic composition of the groundwater resources of the Po River delta (Northern Italy): implications for the environmental impact, J. Hydrogeology, 17, 467-480. Regione Emilia-Romagna (Open file). La banca dati geog- nostica; http://ambiente.reg ione.emilia-romagna.it/ geolog ia/cartograf ia/webg is-banchedati/banca-dati- geognostica Regione Emilia-Romagna and ENI-AGIP (1998). Riserve idriche sotterranee della Regione Emilia-Romagna, G. Di Dio (ed.), S.EL.CA., Firenze, 120 pp. Regione Emilia-Romagna GdL-liquefazione (2012). Primo rapporto sugli effetti della liquefazione osservati a S. Carlo, frazione di S. Agostino (Provincia di Ferrara); http://ambiente.regione.emilia-romagna.it/geologia/ temi/sismica/liquefazione-gruppo-di-lavoro/rapporto_ sancarlo.pdf Scognamiglio, L., L. Margheriti, F.M. Mele, E. Tinti, A. Bono, P. De Gori, V. Lauciani, F.P. Lucente, A.G. Mandiello, C. Marcocci, S. Mazza, S. Pintore and M. Quintiliani (2012). The 2012 Pianura Padana Emiliana seimic sequence: lo- cations, moment tensors and magnitudes, Annals of Geo- physics, 55 (4); doi:10.4401/ag-6159. Sirovich, L. (1996). Repetitive liquefaction at a gravelly site and liquefaction in verconsolidated sands, Soils Found., 36 (4), 23-34. TC4 (1999). Manual for Zonation on Seismic Geotechnical Hazards (Revised Version), The Japanese Geotechnical Society. Tertulliani, A., L. Arcoraci, M. Berardi, F. Bernardini, B. Brizuela, C. Castellano, S. Del Mese, E. Ercolani, L. Graziani, A. Maramai, A. Rossi, M. Sbarra and M. Vecchi (2012). The Emilia 2012 sequence: a macroseismic survey, Annals of Geophysics, 55 (4); doi:10.4401/ag-6140. Toscani, G., P. Burrato, D. Di Bucci, S. Seno and G. Valensise (2009). Plio-Quaternary tectonic evolution of the North- ern Apennines thrust fronts (Bologna-Ferrara section, Italy): seismotectonic implications, B. Soc. Geol. Ital., 128 (2), 1-10. Tsuchida, H. (1971). Estimation of liquefaction potential of sandy soils, In: Proc. 3rd Joint Meeting, US-Japan, UNJR. Youd, T.L., and D.M. Perkins (1978). Mapping of Liquefac- tion induced Ground Failure Potential, J. Geotech. Eng. Div., 104, 433-446. Wakamatsu, K. (1993). History of Soil liquefaction in Japan and Assessment of Liquefaction Potential based on Geo- morphology, Ph.D. thesis, Waseda University, Tokyo, Japan, 245 pp. Wotherspoon, L., M. Pender and R. Orense (2012). Rela- tionship between observed liquefaction at Kaiapoi fol- lowing the 2010 Darfield earthquake and former channels of the Waimakariri River, Eng. Geol., 125, 45-55. *Corresponding author: Riccardo Caputo, Università di Ferrara, Dipartimento di Fisica e Scienze della Terra, Ferrara, Italy; email: rcaputo@unife.it. © 2012 by the Istituto Nazionale di Geofisica e Vulcanologia. All rights reserved. PAPATHANASSIOU ET AL. 742 << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles false /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Warning /CompatibilityLevel 1.3 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.1000 /ColorConversionStrategy /LeaveColorUnchanged /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams true /MaxSubsetPct 100 /Optimize false /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments false /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile (None) /AlwaysEmbed [ true /AndaleMono /Apple-Chancery /Arial-Black /Arial-BoldItalicMT /Arial-BoldMT /Arial-ItalicMT /ArialMT /CapitalsRegular /Charcoal /Chicago /ComicSansMS /ComicSansMS-Bold /Courier /Courier-Bold /CourierNewPS-BoldItalicMT /CourierNewPS-BoldMT /CourierNewPS-ItalicMT /CourierNewPSMT /GadgetRegular /Geneva /Georgia /Georgia-Bold /Georgia-BoldItalic /Georgia-Italic /Helvetica /Helvetica-Bold /HelveticaInserat-Roman /HoeflerText-Black /HoeflerText-BlackItalic /HoeflerText-Italic /HoeflerText-Ornaments /HoeflerText-Regular /Impact /Monaco /NewYork /Palatino-Bold /Palatino-BoldItalic /Palatino-Italic /Palatino-Roman /SandRegular /Skia-Regular /Symbol /TechnoRegular /TextileRegular /Times-Bold /Times-BoldItalic /Times-Italic /Times-Roman /TimesNewRomanPS-BoldItalicMT /TimesNewRomanPS-BoldMT /TimesNewRomanPS-ItalicMT /TimesNewRomanPSMT /Trebuchet-BoldItalic /TrebuchetMS /TrebuchetMS-Bold /TrebuchetMS-Italic /Verdana /Verdana-Bold /Verdana-BoldItalic /Verdana-Italic /Webdings ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 150 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.10000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 150 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.10000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.08250 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile (None) /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName (http://www.color.org) /PDFXTrapped /Unknown /CreateJDFFile false /SyntheticBoldness 1.000000 /Description << /ENU (Use these settings to create PDF documents with higher image resolution for high quality pre-press printing. The PDF documents can be opened with Acrobat and Reader 5.0 and later. These settings require font embedding.) /JPN /FRA /DEU /PTB /DAN /NLD /ESP /SUO /NOR /SVE /KOR /CHS /CHT /ITA >> >> setdistillerparams << /HWResolution [2400 2400] /PageSize [595.000 842.000] >> setpagedevice