Vol51,2_3,2008 491 ANNALS OF GEOPHYSICS, VOL. 51, N. 2/3, April/June 2008 Key words Microzonation, Umbria-Marche earthquake, seismic codes 1. Introduction The working hypothesis behind microzona- tion studies is that we can study a limited num- ber of ideal situations prior to an earthquake oc- currence to forecast the areas where the damage will be increased by amplification of seismic motion or by secondary site effects (liquefac- tion, soil settlement, slope instability, etc.). As often it happens in seismology, also for micro- zonation each earthquake is an occasion to learn new lessons, and the Umbria-Marche seis- mic sequence occurred in September 1997 is no exception. In the 10 years elapsed since that event, several microzonation studies were car- ried out in the area. Immediately after the event, the studies were aimed to the most damaged area and to support emergency intervention and reconstruction plans. In the following years, re- gional and national projects were aimed to transfer the lessons learned to other towns (among the others Spoleto, Cerreto di Spoleto, Città di Castello, Gubbio, Fabriano, Cagli, Treia, Serra de’Conti, Offida, Senigallia). Hav- Codes, models and reality: reductionism vs. holism in a review of microzonation studies in the Umbria-Marche region Marco Mucciarelli Dipartimento di Strutture, Geotecnica e Geologia Applicata, Università della Basilicata, Potenza, Italy Abstract In the 10 years since the Umbria-Marche earthquake, several microzonation studies were carried out in the two regions. In the immediate aftermath of the event, the focus was on the epicentral area and toward emergency in- tervention and reconstruction plans. In the following years, regional and national projects aimed to transfer the lessons learned to other towns. Usually, those two kinds of microzonation studies are referred to as «simplified» and «detailed». The difference is more subtle, and leads to the question of whether a microzonation study can be tackled following a reductionist approach, i.e. leaving different experts taking care of a limited subject (geol- ogy, geophysics, seismology, geotechnics, structural engineering). The impression looking back at 10 years of studies is that a holistic approach would be more appropriate to describe a system (structure-soil-bedrock) that is non-linear, inhomogeneous, and presenting feedback among its components. A second problem that emerges is the link with codes and practitioners. During the past 10 years the seismic code has been changed and a new version is on arrival. The last proposed version of the code is based on a parameter (Vs30) that is discussed in the same country where it was first adopted, and introduces a parameter (acclivity) that appears to be a second- or third order problem with respect to others that are completely disregarded (e.g., 2-d site effects). A possible explanation for this mismatch between codes, models and reality is that our knowledge of distribution and am- plitude of site effects is biased by selective under-sampling. Being driven by damage, and paying less attention to a uniform distribution of studied sites and situations, we act like a drunk man looking of his lost keys under a street lamp, not because he is sure that he lost them there, but because the light is there. Mailing address: Dr. Marco Mucciarelli, Dipartimento di Strutture, Geotecnica e Geologia Applicata, Università della Basilicata, Potenza, Italy; e-mail: marco.mucciarel- li@unibas.it. Vol51,2_3,2008 4-03-2009 10:28 Pagina 491 492 M. Mucciarelli ing more time (and sometimes more resources) available, the second kind of studies were car- ried out performing more experimental work in situ and more detailed modelling, while the first studies used a limited amount of geophysical and geotechnical data and rely mostly upon ge- ological maps and 1-d models. Usually, those two kinds of microzonation studies are referred to as «simplified» and «de- tailed». The difference is more subtle, and leads to the question if a microzonation study can be tackled following a reductionist approach, i.e. leaving different experts taking care of a limit- ed subject (geology, geophysics, seismology, geotechnics, structural engineering) and draw- ing conclusions from that. A further need of microzonation studies is to produce something that can be transferred to codes and practitioners, and here the problem is to (quoting A. Einstein) «make things as simple as possible, but not simpler». This prob- lem was born with modern seismology. The year 2007 marks not only the 10th anniversary of the 1997 Umbria-Marche quake but also the 150th anniversary of the 1857 Val d’Agri quake. After a field trip in the affected area, Robert Mallet wrote its book «First principles of ob- servational seismology». The problem of which geological and morphological situations could influence the observed damage distribu- tion was already tackled by previous studies, but Mallet is the first to understand the prob- lem of biased distribution of experimental evi- dence. He wrote: «Much pain has been devot- ed by earthquake describers (Hamilton, Dolomieu, Spallanzani, and others) to the question whether towns situated upon the plain on loose material, or those on the solid rocks on hill tops, suffered most; and their discus- sions evince much obscurity of thought, having no guiding principle.... Were the whole of the facts ... discussed with regard to this question, it would be found, that more places were de- stroyed upon the rock, than upon loose clay or other such foundations, and more upon hills than upon plains; but this would not conduct us to the whole truth, because in all South Italy, there are a great many more places upon rock and upon hills, than upon clays and low plains». 2. Review of simplified microzonation studies A work that takes into account the problem outlined by Mallet was published by Tertulliani (2000). This paper examines all the localities af- fected by the first two shock of the Umbria- Marche sequence with a damage reaching VI-VII degree in the European Macroseismic Scale. A simplified geo-morphological scheme is then proposed, distinguishing town located either a) on alluvium near the basin edge; (b) on fan de- bris; (c) on debris-alluvium transition; (d–e) on bedrock debris transition; (f) on bedrock-alluvi- um transition; (g) on bedrock; (h) on alluvium far from the basin edge. The conclusion is that «heavy damage occurred in hamlets settled on al- luvial deposits near to the edge of the basin, with a border fault to mark the bedrock limit» while «the severity gradually decreases as we move to villages settled on debris and on a bedrock-debris transition. The lowest intensity is found on bedrock sites and on the inner part of the basin». Stratigraphy appears thus to be a governing factor more important than morphology. The most extensive microzonation study aimed to emergency purposes was performed by Pergalani et al. (1999). In this work «a team of surveyors were trained to collect field informa- tion such as geologic and geomorphologic fea- tures and, where possible, pre-existing geotech- nic or geophysic information. Information was collected and analysed with the aid of dynamic codes to calculate the possible local site effects». Using 60 villages as a training set, the authors produced a generalization of the stratigraphic and morphologic situations, synthesized in a table grouping morphology types (valley, ridges or slopes) and lithologic units and thickness, and finally assigning to each zone a frequency-inde- pendent amplification factor according with the provision of the then enforced seismic code. The conclusion was that «the procedure ... prove to be a suitable tools for site effects assessment in post-event situation when time and budget con- straint are present; in fact about 1000 sites have been surveyed and an estimate of the possible amplification established in less than six months. The average cost was also relatively small: about US$1500 for each site». Some years later, Mar- zorati et al. (2003) performed a work aimed to Vol51,2_3,2008 4-03-2009 10:28 Pagina 492 493 Codes, models and reality: reductionism vs. holism in a review of microzonation studies in the Umbria-Marche region «evaluating the accuracy of that methodology, following a completely different approach. The actual damages, detected by an aerial photo in- terpretation, are compared to the theoretical damages estimated through empirical curves». The conclusions are particularly interesting from the point of view of the knowledge transfer to practitioners: «most of the errors detected in the comparison are... not ascribed to the methodolo- gy, but are simply due to the incorrect applica- tion of the reference table proposed for the zona- tion, or to the inaccurate geologic and geomor- phologic surveys. Those erroneous evaluations could be justified only if considering that the mi- crozoning had to be carried out so rapidly that several surveyors not adequately trained were employed... Other inaccuracies revealed new ge- omorphological and geological settings that have to be included in the reference table, like the sad- dle-like morphology and the surface faulting». Among the single-technique approaches to microzonation, two methodologies were tested on large data sets for the first time in Italy. Goretti and Dolce (2002) used an objective methodology to identify damage enhancements from observed data, called Normalised Damage Index: «the distribution of the seismic intensity that actually affected buildings is obtained by means of well-assessed vulnerability functions, together with a Bayesian approach, and taking into account uncertainties in building type and building damage». Mucciarelli and Monachesi (1998) per- formed a series of measurements using the HVSR technique (Nakamura, 1989) following the suggestion about how to collect and process data provided by Mucciarelli (1998). The meas- urements were taken in 12 sites arranged in cou- ple of places where a strong difference in inten- sity occurred within a short distance. Then a cor- relation was sought between the intensity varia- tion and the amplification ratio in different band of frequency. The result was that the higher cor- relation occurred in the range 1-2 Hz. 3. Review of detailed microzonation studies From 1997 to 2007, several studies for mi- crozonation were performed in towns located in the Umbria and Marche regions. The main studies are listed in the following: Nocera Um- bra (Cattaneo and Marcellini, 2000), Fabriano (Marcellini and Tiberi, 2000), Cagli, Treia, Ser- ra de’Conti, Offida (Mucciarelli and Tiberi, 2004), Città di Castello (Bordoni et al. 2003), Cerreto di Spoleto (Martino et al., 2006) Seni- gallia (Mucciarelli and Tiberi, 2007), Spoleto (Vuan et al., 2007), Gubbio (Fiorini et al., 2007). These works mixed in a different degree geological, geophysical and geotechnical data, but are all united by a common feature: the use of recorded earthquakes. The main results of these works concern three points: 1) the importance of 2-d effects in sedimentary basins; 2) mixed experimental and numerical evidence of strong morphologic ef- fects; 3) the importance of soil-building reso- nance for damage enhancement. 3.1. 2-d effects The strongest amplification recorded for the Umbria-Marche earthquake occurred in the Gubbio basin (Castro et al., 2004). In gen- eral, alluvial plains were affected by high am- plification, as demonstrated also for Città di Castello basin (Bordoni et al., 2003) and Colfiorito basin (Di Giulio et al., 2003). The importance of this kind of effects is well known to seismologist and engineers, but nothing has been done so far for their inclu- sion in the seismic codes, notwithstanding proposals aimed to this goal (Chavez-Garcia and Faccioli, 2000). 3.2. Morphologic effects There has been mixed evidence about the role of morphology in amplifying recorded seismic ground motion in the Umbria-Marche area. Laurenzano et al. (2008) noted no ampli- fication on the Treia ridge but a strong amplifi- cation in the underlying valley. Marsan et al. (2004) analysing the data from two arrays, «ob- served a large amplification in the basin border of Cesi site, in opposite small amplification due to topografic effect in Sellano site». Martino et Vol51,2_3,2008 4-03-2009 10:28 Pagina 493 494 M. Mucciarelli al. (2006) studied the sharp ridge of Cerreto di Spoleto. They deployed a seismological array «in order to highlight possible topographic am- plification effects all along the ridge; the find- ings, however, at least with the low energy lev- el of the recorded ground motion, do not point out any topographic amplification» and after the study of the rock mass fracturation, they concluded that «rock sites can show seismic amplification effects other than the topographic one. In particular, the seismic response in a rock mass ridge can be affected by rock mass joint- ing: adjacent rock masses with significantly dif- ferent jointing and specific geometries seem to favour trapped wave amplification, if persist- ence along depth of these characteristics oc- curs». The amplification due to waves trapped inside a fault zone was already observed and modelled by Rovelli et al. (2002) for Nocera Umbra. The Italian seismic code disregarded the problem of amplification at rock sites, abol- ishing the hard-rock class of the NEHRP classi- fication scheme (Vs30>1500 m/s). On the con- trary, the new provisions introduce topographic amplification for all the slopes with more than 15° inclination. 3.3. Soil-building resonance During the Umbria-Marche sequence, the enhancement of damage due to soil-building resonance was studied for first time in Italy with digital instrumentation. Some of the most peculiar damage patterns observed were ascribed to soil-building reso- nance, as for the case of damage to reinforced concrete buildings in a limited area of Fabriano (Marcellini and Tiberi, 2000; Mucciarelli et al., 2001;) or the damage to retrofitted buildings in the historical centre of Sellano (Mucciarelli et al., 2001; Nunziata and Natale, 2004). The case of these last two works is interesting because they reach similar conclusions starting from in- dependent dataset and models. Mucciarelli et al. (2001) performed empirical measurement of ambient vibration on both soil and buildings, concluding that «the amplifcation is mainly due to a 2-d hill-top effect, with variability due to anthropic fills. It is clear that the damage is greater for buildings whose frequency ap- proaches the site transfer function». Nunziata and Natale (2004) aimed to reproduce accelero- metric recording using numerical modelling and concluded that «taking into account the building typology of 2–3 floors, the high dam- age level might be attributed not only to the de- graded conditions of the brick houses but also to 2D amplification effects, as the frequencies of the maximum peaks of the spectral amplifi- cation correspond to the eigenfrequencies of the buildings». Notwithstanding this evidence, no mention is made in the Italian code about the possibility of taking into account the frequency of soil and building for the existing structures to be retrofitted. 4. Discussion Other events in Italy (like the 2002 Molise earthquake) were tackled using methodologies that were first applied after the Umbria- Marche sequence. The data collected during Umbria-Marche microzonation campaigns were a valuable source of information for fur- ther studies. For example, shear waves velocity estimates contributed with those performed in other re- gions to verify the reliability of Vs30 as a proxy of seismic amplification (Mucciarelli and Gallipoli, 2006; Gallipoli and Mucciarelli, 2009). These studies concluded that «Vs30 ... is not a good proxy of observed amplification ef- fects in about one third of the sites. The reason why in Italy Vs30 does not provide satisfactory estimates is linked to peculiar geological set- tings that are widespread in the country. Vs30 seems to work fine only if a site has a strictly monotonic velocity profile increasing with depth and a strong impedance contrast in the first dozen meters». The data from Umbria-Marche helped also to investigate the effect of selective under-sam- pling mentioned in the introduction (Mucciarel- li and Gallipoli, 2004). From 1997 to 2004 years, our group conducted two kinds of HVSR campaigns: microzonation studies and post earthquake surveys. The aim was different: in the first case, we sampled all the municipalities Vol51,2_3,2008 4-03-2009 10:28 Pagina 494 495 Codes, models and reality: reductionism vs. holism in a review of microzonation studies in the Umbria-Marche region or sites within the study area without a priori selection; in the second case, we performed measurements aimed at investigating possible resonance phenomena between the fundamen- tal frequency of the soil and the one of damaged buildings. Thus in this second case we per- formed a selection driven by the observed dam- age. The existence of hidden selection criteria in a sample may introduce unwanted bias in the outcome. We performed all the measurement using the same instrumentation and processing technique. With a database of 540 HVRS meas- urements, we compared three samples: 407 measurements from microzonation studies in different Italian regions (Umbria, Marche, Basilicata, Molise), 79 from post- earthquake surveys (Umbria-Marche, Slovenia, Calabro- Lucano, Izmit) and 54 from our latest post- earthquake survey (Molise). In this last case, we performed measurements not only at the most damaged sites, but also more uniformly in the investigated area (Gallipoli et al., 2004). The variable we considered for our statistics is the highest HVSR value in the frequency range 0.5-10 Hz. It is clear from fig. 1 that the data from microzonation studies and post-quake sur- vey have similar distributions but very different parametrisation. The main result from fig. 1 is the probability of observing high HVSR values at randomly sampled sites. The data sampled during microzonation studies have a lower probability than the ones resulting from post- quake studies. This means that there are few sites with strong local amplification, and they are those that claim for attention when the dam- age is driving our researches. The distribution observed after the Molise 2002 event shows that a study conducted after an earthquake yield the same distribution of earthquake-independ- ent studies, provided that the sampled location are not pre-selected but the whole area is inves- tigated. Fig. 1. Comparison between the empirical cumulative distribution functions of different sets of HVSRs, as a function of the largest value observed in the frequency range 0.5-10 Hz. See text for details. 100 90 80 70 60 50 40 30 20 10 0 2 4 6 0.5 Hz < Max HVSR <10 Hz Microzonation Post-earthquake Molise 2002 % 80 10 Vol51,2_3,2008 4-03-2009 10:28 Pagina 495 496 M. Mucciarelli 5. Conclusions The Umbria-Marche sequence provided several valuable lessons about microzonation studies. These can be summarised in the follow- ing: 1. no simple method is safe enough if the practitioners are not well-trained to use it; 2. the strongest amplifications recorded for the Umbria-Marche earthquake were due to 2-d basin effects; 3. topographic amplification, if and where present, is a second-order problem with respect to impedance contrast amplification; 4. the most peculiar damage patterns ob- served were due to soil-building resonance; 5. Vs30 is not a good proxy of observed am- plification effects in about one third of the cas- es. Very little of this information was taken in- to account for the formulation of the Italian seismic code. A leitmotiv often repeated is that the codes must be simple or the practitioners could have problems in applying them. One may wonder if any of us would oppose the diffusion of medical ultrasonography or computed tomography scans on the basis that old-fashioned doctors are familiar just with their wooden stetho- scopes. A wrong diagnosis is never good, not only if the illness is worst than assumed, but al- so in case of over-conservatism, causing unnec- essary expense and routing limited resources where they are not needed. Ten years after the Umbria-Marche quake, but also ten years after the introduction of the NEHRP-FEMA provisions, the world outside Italy has treasured the experiences gathered in the meanwhile. – The relationship between acclivity and amplification has been revised by Wald and Allen (2007) who propose using the topograph- ic map as a first order support for determining amplification via Vs30 since «topographic vari- ations should be an indicator of near-surface geomorphology and lithology ... with steep mountains indicating rock, nearly flat basins in- dicating soil, and a transition between the end members on intermediate slopes». – The presence of amplification in sites with 1500