Lessons from April 6, 2009 L’Aquila earthquake to enhance microzoning studies in near-field urban areas

This study focuses on two weak points of the present procedure to carry out microzoning study in near-field areas: (1) the Ground Motion Prediction Equations (GMPEs), commonly used in the reference seismic hazard (RSH) assessment; (2) the ambient noise measurements to define the natural frequency of the near surface soils and the bedrock depth. The limitations of these approaches will be discussed throughout the paper based on the worldwide and Italian experiences performed after the 2009 L’Aquila earthquake and then confirmed by the most recent 2012 Emilia Romagna earthquake and the 2016–17 Central Italy seismic sequence. The critical issues faced are (A) the high variability of peak ground acceleration (PGA) values within the first 20–30 km far from the source which are not robustly interpolated by the GMPEs, (B) at the level 1 microzoning activity, the soil seismic response under strong motion shaking is characterized by microtremors’ horizontal to vertical spectral ratios (HVSR) according to Nakamura’s method. This latter technique is commonly applied not being fully compliant with the rules fixed by European scientists in 2004, after a 3-year project named Site EffectS assessment using AMbient Excitations (SESAME). Hereinafter, some “best practices” from recent Italian and International experiences of seismic hazard estimation and microzonation studies are reported in order to put forward two proposals: (a) to formulate site-specific GMPEs in near-field areas in terms of PGA and (b) to record microtremor measurements following accurately the SESAME advice in order to get robust and repeatable HVSR values and to limit their use to those geological contests that are actually horizontally layered.


Introduction
On April 6, 2009, at 1:32 a.m. (local time) an M w 6.3 earthquake with shallow hypocentral depth (8.3 km) hit the city of L'Aquila and several municipalities within the Aterno Valley. This earthquake can be considered one of the most mournful seismic event in Italy since 1980 although its magnitude was moderately-high: 308 fatalities and 60.000 people displaced (data source http://www. protezionecivile.it) and estimated damages for 1894 M€ (data source http://www.ngdc.noaa.gov). These numbers showed how dangerous can be an unexpected seismic event in urbanized territories where no preventive actions have been addressed to reduce seismic risk. Hence, meanwhile, some actions were implementing in the post-earthquake time such as updating the reference Italian hazard map (by the Decree OPCM n. 3519 on 28 April 2006) and drawing microzoning maps to be used in the reconstruction stage, two major earthquakes struck the Emilia Romagna Region (in the Northern part of Italy), causing 27 deaths and widespread damage. The first, with M w 6.1, occurred on 20 May at 04:03 local time (02:03 UTC) and was located at about 36 km north of the city of Bologna. Then, a second major earthquake (M w 5.9) occurred on 29 May 2012, in the same area, causing widespread damage, particularly to buildings already weakened by the 20 May earthquake. Later on, the 2016-17 Central Italy Earthquake Sequence occurred, consisting of several moderately-high magnitude earthquakes between M w 5.5 and M w 6.5, from Aug 24, 2016, to Jan 18, 2017, each centered in a different but close location and with its own sequences of aftershocks, spanning several months. The seismic sequence killed about 300 people and injured the other 396. Worldwide, several other strong earthquakes (eg. 1998Northridge earthquake, 2004Parkfield earthquake, 2010Canterbury and 2011and 2017Christchurch earthquakes, 2018 Sulawesi earthquake) produced devastating effects in the same time span. All these events show the need to carry out efficient microzoning studies to plan vulnerability reductions of urban structures and promoting the resilience of the human communities in seismic territories. After the 2009 L'Aquila earthquake, the microzoning studies have been introduced in Italy by law and the Guidelines for Seismic Microzonation (ICMS 2008) have been issued to accomplish these studies according to the most updated international scientific findings. These guidelines provide the local administrators with an efficient tool for seismic microzoning study to predicting the subsoil behavior under seismic shaking. Unfortunately, the ICMS does not give special recommendations for urbanized near field areas (NFAs). The microzoning activity concerning urbanized territories as suggest by ICMS (2008) is made up of four steps: 1-Estimating the reference seismic hazard to provide the input peak horizontal ground acceleration (PGA) at each point on the national territory and the normalized response spectrum at each site. 2-Dynamic characterization of soil deposits overlaying the seismic bedrock at each urban center in order to draw the microzoning maps (MM) at three main knowledge levels. 3-The Level 1 MM consists of geo-lithological maps of the surficial deposits that show typical successions and the amplified frequency map drawn through the measurements of microtremors elaborated by horizontal to vertical spectral ratio HVSR Nakamura's technique. Nakamura's method (1989), the horizontal to vertical noise components are calculated to derive the natural frequency of surficial soft deposits and their thickness. 4-The Level 2/3 MM consists of drawing maps after performing the numerical analyses of (a) seismic local amplification factors in terms of acceleration FA and velocity FV; (b) liquefaction potential LP and (c) permanent displacements due to seismically induced slope instability.
After 10 years of training the ICMS and the related methods, it is now the time to start analyzing some arisen weak points. Starting from the large data acquired worldwide on recent strong motion earthquakes, the experiences developed in seismic hazard assessment, and the sitespecific seismic response characterization carried out by the writing authors after the 2009 L'Aquila earthquake, the aforementioned weak points (related to some aspects of the steps 1 and 3) are hereinafter discussed and some proposals are made to improve the efficiency of the microzoning studies especially in NFAs.
In this paper, after a brief background section on the procedures to accomplish the reference seismic hazard assessment (background section), the methods to calculate the Ground Motion Prediction Equations (GMPEs) and the HVSR (Nakamura's method) are briefly recalled in section 2. Then in section 3, the results from observations of recorded peak ground acceleration (PGA) and pseudospectral acceleration (PSA) values within the NFAs from the L'Aquila earthquake and other worldwide strong earthquakes have been discussed. In addition, some applications of Nakamura's procedure to characterize the natural frequency of the sites throughout the Aterno Valley have been discussed. Finally, in the conclusion section, some relevant points drawn from the discussed microzoning experiences have been highlighted to improve the efficiency of the microzonation studies in urban centers especially located in NFAs.

Background on seismic hazard assessment
Several theoretical and experimental studies performed worldwide in the last 50 years (see Kramer 1996 and the reference herein), highlighted that seismic shaking intensity is due to the magnitude of the earthquake generated at the source, to the travel paths of the seismic waves from the source to the buried or outcropping bedrock (that is called reference seismic hazard RSH) and the additional phenomena of local amplification or de-amplification take place where soil deposits overlay the rocky bedrock, named local seismic response LSR (Paolucci 2002;Vessia and Venisti 2011;Vessia and Russo 2013;Vessia et al. , 2017Boncio et al. 2018, among others). The RSH maps drawn worldwide on national territories do not take into account the results of LSR studies.
The pioneering work by Signanini et al. (1983) after the 1979 Friuli earthquake confirmed the observations on the ground: local seismic effects could enlarge the referenced hazard at a site by 2-3 times in terms of MCS scale Intensity but also in PGA values owing to the local morphological and stratigraphic settings. Such RSL is particularly evident in near field areas, from then on named NFAs. The NFAs have been defined among others by Boore (2014a) as the Fault Damage Zones. These areas cannot be uniquely identified depending on the source rupture mechanisms, the surficial soil deposits and the multiple calculation methods used for measuring the distance between the seismic stations and the source. Especially in these areas, about the first 30 km aside the source, spot-like amplifications are the common amplification pattern captured through the Maximum Intensity Felt maps. These maps estimate the differentiated damages suffered by buildings and urban structures by means of the Macroseismic Intensity scale (e.g. Mercalli-Cancani-Sieberg MCS scale, European Macroseismic scale EMS, Modified Mercalli Intensity MMI scale). One of the Maximum Intensity Felt maps on the Italian territory was drawn by Boschi et al. (1995). They took into account the seismic events that occurred from 1 to 1992 AD with a minimum intensity felt of VI MCS. This latter value is the one commonly used to highlight those areas where seismic events caused relevant damages to dwellings and infrastructures, ranging from severe damages to collapse. Boschi et al. (1995) map is reported in Fig. 1: it showed IX-X MCS at L'Aquila district based on historical earthquakes that are in agreement with the seismic intensity map drawn by Galli and Camassi (2009) after the mainshock of 2009 L'Aquila earthquake. This map is also in very good agreement with other recent earthquakes such as the 2012 Emilia Romagna and 2016-  Boschi et al. 1995, modified) with the areas of two recent Italian earthquake sequences considered in the present work 17 Central Italy earthquakes (Fig. 1). Moreover, Boschi et al. (1995), Midorikawa (2002 and more recently, Paolini et al. (2012) proposed a direct use of the Maximum Felt Intensity maps to highlight those areas where the reference seismic hazard is largely increased by the local amplificated responses of soil deposits, that is the NFAs.
The most used method to perform the reference seismic hazard assessment has been conceived in the late '60s. It is Cornell's method (1968) that was implemented into a numerical code by Mc Guire (1978). Cornell (1968) introduced the Probabilistic Seismic Hazard Assessment (PSHA) method to carry out the reference seismic hazard at a site considering the contribution of the seismic source and the travel path of the seismic waves by considering the uncertainties related to these estimations. This method consists of four steps (Kramer 1996), as illustrated in Fig. 2: STEP 1. To identify the seismogenic sources as single faults and faulting regions in terms of magnitude amplitude generated at different time spans. The probabilistic approach to such a characterization needs to know the rate of the earthquake at different magnitudes at the site and the spatial distribution of the fault segment or the source volume that can be activated. STEP 2_1. To calculate the seismic rate in a region the Gutenberg-Richter law is used, where a and b coefficients are drawn by interpolating numerous data from a database of seismic events (instrumental and noninstrumental) available for a limited number of source areas and affected by the lack of completeness distortions. The earthquake occurrence probability is estimated by means of a Poisson distribution over time that is independent of the time span of the last strong seismic event. STEP 2_2. To calculate the spatial distribution of the seismic events alongside a fault zone is a character difficult to get known and the spatial distribution of earthquake sources within a seismogenic area is commonly assumed uniformly distributed. STEP 3. To define the ground motion prediction equations GMPEs that enable to predict, at different magnitude ranges, the decrease with the distance from the seismic source of the strong motion parameter assumed to be representative of the earthquake at a site. STEP 4. To calculate the probability of exceedance of a target shaking value of the considered ground motion parameter, i.e. PGA, in a time span at a chosen site, due to the contribution of different seismogenic sources.
The previous 4 steps attempt to take into account several sources of uncertainties, such as the limited knowledge about the fault activity, the qualitative and documental estimations of the past earthquake effects at the sites, the lack of completeness of the seismic catalogs (meaning that the database of the seismic events is populated by several data related to both low and high magnitude ones) and the dependency among the recorded strong seismic events. In addition, the distortions in Gutenberg-Richter law, defined for different regions worldwide and the uncertainties related . Logic trees are commonly used to take into account different formulations of GMPEs and several Gutenberg-Richter rates of magnitude occurrence (Kramer 1996). Molina et al. (2001) pointed out that the PSHA has its strength in the systematic parameterization of seismicity and the way in which also epistemic uncertainties are carried out through the computations into the final results.
Recently, alternative approaches to PSHA calculation have been suggested, such as through extensions of the zonation method (Frankel 1995;Frankel et al. 1996Frankel et al. , 2000Perkins 2000) where multiple source zones, parameter smoothing and quantification of geology and active faults have been successfully applied. The Frankel et al. (1996) method applied a Gaussian function to smooth a-values (within Gutenberg-Richter law) from each zone, thereby being a forerunner for the later zonation-free approaches of Woo (1996). This latter approach tries to amalgamate statistical consistency with the empirical knowledge of the earthquake catalogue (with its fractal character) into the computation of seismic hazard. Furthermore, Jackson and Kagan (1999) developed a non-parametric method with a continuous rate-density function (computed from earthquake catalogues) used in earthquake forecasting. Nonetheless, all these methods need a function to propagate the strong motion parameter values from the source to the site under study. To this end, the GMPEs are built by interpolating large databases of seismic records (related to specific geographical and tectonic environments worldwide), taking into account the contributions of the earthquake magnitude M and the distance to the seismic source R, according to the following form (Kramer 1996): where Y is the ground motion parameter, commonly the peak ground horizontal acceleration PGA or the spectral acceleration SA at fixed period; S i is related to the source and site: they are the refinement terms due to the enlargement of the seismic databases and the possibility of drawing specific regional GMPEs.

The uncertainty of GMPEs in near field areas
Several examples of GMPEs are provided in literature (Kramer 1996 among others) while a recent throughout review of several possible formulations of GMPEs used in the USA can be found at the Pacific Earthquake Engineering Research center PEER website http://peer.berkeley. edu/publications/peer_reports_complete.html. In Italy, the GMPEs are built based on the PGAs drawn from the Italian shape wave database of strong motion events (Faccioli 2012;Bindi et al. 2011Bindi et al. , 2014Cauzzi et al. 2014). Commonly, the PGA values represent the strong motion parameter used in microzoning studies but the related GMPEs are highly uncertain especially in the first tens of kilometers as shown in Fig. 3a (Faccioli 2012) and Fig. 4a (Boore 2013). According to Boore (2013Boore ( , 2014a, fault zone records show significant variability in amplitude and polarization of PGA, SA especially at low periods (as shown in Fig. 4a) and magnitude saturation beyond M w 6, although the causes of this variability are not easy to be unraveled. The main drawback of the GMPEs is the weakness of their predictivity at a short distance from the seismic source due to two main issues affecting the NFAs worldwide: 1) a few seismic stations installed; 2) highly scattered measures of strong motion parameters, especially in terms of accelerations (i.e. PGA, SA, etc) ( Fig. 2, step 3), that do not show any decreasing trend with distance.
The PGA spatial uncertainties have been observed after several recent strong earthquakes, such as 2009 L'Aquila earthquake (Lanzo et al. 2010;Bergamaschi et al. 2011;Di Giulio et al. 2011Christchurch and 2010 Darfield earthquakes in New Zealand (Bradley and Cubrinovski 2011) (Fig. 5), 2012 Emilia Romagna earthquake in Italy, 1994 Northridge earthquake in USA (Boore 2004) and 2013 Fivizzano earthquake (Fig. 5). In the case of the 2009 L'Aquila earthquake (Fig. 4b), the areal distribution of PGAs around the source seems to be highly random although they show that the most dramatic increase occurs where thick soft sediments are met over rigid bedrocks or where bedrock basin shapes can be recognized. This latter traps the seismic waves and caused longer duration accelerograms with increased amplitudes at short and moderate periods (lower than 2 s) (Rainone et al. 2013).
The GMPEs based on PGAs tend to saturate for large earthquakes as the distance from the fault rupture to the observation point decreases. Boore (2014b) showed that the PGA parameter is a poor measure of the groundmotion intensity due to its non-unique correspondence to the frequency and acceleration content of the shaking waves (Fig. 4), especially at high frequencies. Furthermore, Bradley and Cubrinovski (2011) and Boore (2004) stated that the influence on the amplitude and shape response by local surface geology and geometrical conditions is noted to be much more relevant than the forward directivity and the source-site path on spectral accelerations in near field areas and for periods shorter than 3 s. Thus, the GMPEs of PGAs within the NFAs are highly uncertain and cumbersome to be predicted even when fitted on single seismic events as shown in  (Fig. 3a). This GMPE was built based on the ITACA 2010 database (Luzi et al. 2008, http://itaca.mi.ingv.it/ItacaNet_30/#/home) that collects Italian strong motion shape waves. It is worth noticing  that within the first tens of kilometers from the source, these data indicate that the GMPEs are not accurate in predicting the PGA values. To avoid the pitfalls in GMPEs based on peak parameters, integral ground motion parameters have been proposed in the literature (Kempton and Stewart 2006;Abrahamson and Silva 2008;Campbell and Bozorgnia 2012), such as Arias intensity (AI) and cumulative absolute velocity (CAV). In addition, Hollenback et al. (2015) and Stewart et al. (2015) formulated new generation GMPEs based on median ground-motion models as part of the Next Generation Attenuation for Central and Eastern North America project. They provided a set of adjustments to median GMPEs that are necessary to incorporate the source depth effects and the rupture distances in the range from 0 to 1500 Km. Moreover, the preceding authors suggest a distinct expression for the GMPE at short-distance to the source (by 10 km), that is: where R RUP is the rupture distance, that is the closest distance to the earthquake rupture plane (km); c 1 and c 2 are the regression coefficients and h is a "fictitious depth" used for ground-motion saturation at close distances.

Ambient noise measures elaborated by means of the Nakamura horizontal to vertical ratio HVSR
In 1989 Nakamura proposed to use the ambient noise measurements to derive a seismic property of a site, that is the frequency range of amplification, through the spectral ratio of horizontal H and vertical V ambient vibration (microtremors) components of the recorded signals. If the site does not amplify, the ratio H/V is equal to 1. The Nakamura method shows the advantage to solve the troublesome issue to find out a reference site. In fact, it considers the vertical component as the one that is not modified by the site where horizontal subsoil layers are set and SH seismic waves represent the ambient noise signal content in a quite site (far from urban or industrial areas). This latter is the reference signal whereas the horizontal component is the only one that can be affected by the amplifying properties of the soils. As a matter of fact, Nakamura assumed that: locally random distributed sources of microtremors generate not directional signals almost made up of shear horizontal or Rayleigh waves; the microtremors are confined in the surficial layers because the subsoil is made up of soft layered sediments overlaying a rigid seismic bedrock.
A relevant implication of the Nakamura method is that the peaks of the ratio H/V are related to the presence of high acoustic impedance contrast at the depth h that can be derived by the following expression: where V S is the mean value of the measured shear wave velocity profile and f o is the amplified frequency measured by means of the noise measurement. The fundamental rules to perform a correct ambient noise recording was provided by the European research project named SESAME (Bard and the WG 2004) that analyzed the possible drawbacks of the simple model introduced by the Nakamura method and issued guidelines that offer important recommendations regarding the places where the method can be successfully used in urbanized areas. The given recommendations are based on a rather strict set of criteria, that are essentially composed of (1) experimental conditions and (2) criteria for gaining reliable results (Table 1).
As can be seen from Table 1, the recommendations are focused on the weather conditions that influence the quality of the noise measurements and they highlight the need to record at distance from structures, trees, slopes because all these items affect the records. Unfortunately, it is not possible to quantify the minimum distance from the structure where the influence is negligible, as this distance depends on too many external factors (structure type, wind strength, soil type, etc.). Furthermore, related to the measurement spacing, SESAME guidelines suggest to never use a single measurement point to derive f 0 value, make at least three measurement points. This latter advice is often disregarded. An interesting alternative way to apply the HVSR method is to investigate the "heavy tails" of its statistical distribution that is like that of a critical system (Signanini and De Santis 2012). This is likely indicative of the strong non-linear properties of rocks forming the uppermost crust resulting in a power-law trend.
However, the Nakamura technique has been introduced in microzoning studies at Level 1 in Italy (ICMS 2008) to draw the natural frequency map of urban sites but limitations to suitable sites have not been prescribed. Although the Nakamura method seems to be simple, lost cost and short time consuming, the suitable sites where it can be applied are few, especially in urban centers. This type of indirect investigation method is not applicable in complex geological contexts (e.g. buried inclined fold settings) and is not easily handled for the difficulties in reproducing the same measurements under variable site conditions and noise sources and acquisitions performed by different operators even at the same site. Rainone et al. (2018) undertook a thorough study on the effectiveness of HVSR in predicting amplification frequencies at two Italian urban areas characterized by different subsoil setting and noise distribution. Results from this study show that HVSR works well only where horizontally layered sediments overlay a rigid bedrock: these conditions are the most relevant and the most influential on the predictivity of the actual f o measured values.

A new proposal for GMPEs in near field areas
Recently, a study to formulate ad hoc GMPEs for PGAs within NFAs of the Central Italy Apennine sector has been performed. Only horizontal PGA values measured from seismic stations set on A and A* soil category (V s ≥ 800 m/s), generated by seismic events with M w ranging between 5.0 and 6.5 and normal fault mechanism, have been extracted from ITACA shape wave database (Luzi et al. 2008). The selected events cover the period from 1997 to 2017 and consider 25 seismic events from three strong seismic sequences generated by normal faults: the 1997 Umbria-Marche, 2009 L'Aquila, 2016-2017 Central Recording parameters Level the sensor as recommended by the manufacturer. Fix the gain level at the maximum possible without signal saturation.
In situ soil-sensor coupling Set the sensor down directly on the ground, whenever possible. Avoid setting the sensor on "soft grounds" (mud, plowed soil, tall grass, etc.) or soil saturated after rain.
Artificial soil-sensor coupling Avoid plates from "soft" materials such as foam rubber, cardboard, etc. On steep slopes that do not allow correct sensor leveling, install the sensor in a sand pile or in a container filled with sand. On snow or ice, install a metallic or wooden plate or a container filled with sand to avoid sensor tilting due to local melting.

Nearby structures
Avoid recording near structures such as buildings, trees, etc. in case of wind blowing (faster than about 5 m/s). It may strongly influence H/V results by introducing some low frequencies in the curves. Avoid measuring above underground structures such as car parks, pipes, sewer lids, etc.

Weather conditions
Wind: protect the sensor from the wind (faster than approx. 5 m/s). This only helps if there are no nearby structures. Rain: avoid measurements under heavy rain. Slight rain has no noticeable influence. Temperature: check sensor and recorder manufacturer's instructions. Meteorological perturbations: indicate on the field sheet whether the measurements are performed during a low-pressure meteorological event.

Disturbances
Monochromatic sources: avoid measurements near construction machines, industrial machines, pumps, generators, etc. Transients: in case of transients (steps, cars, etc) increase the recording duration to allow for enough windows for the analysis, after transient removal.
Italy. The studied source area is a quadrant whose edges' coordinates are (43.5°, 12.3°) and (42.2°, 13.6°) in decimal degrees. The PGA measures within the first 35 km from the seismic source have been taken into account. The hypocentral distance has been used to define the source to site distance. Two GMPEs within the first 35 km have been drawn for two ranges of moment magnitude: 5 ≤ M w1 < 5.5 and 5.5 ≤ M w2 ≤ 6.5. These ranges represent the injurious magnitudes of the Italian moderately-high magnitude earthquakes (Fig. 6). As can be noted from Fig. 6 the PGA values seem not to be highly different in the two magnitude ranges and they do not show a clear trend with the hypocentral distance. Thus, these two datasets have been kept distinct and a box and whisker plot has been used to calculate their medians, quartiles, and interquartile distance. Figure 7a, b show the two datasets with a different number of bins of hypocentral distance: it is due to the circumstance that for higher magnitudes (Fig. 7b) the seismic stations within the first 10 km are that few that cannot be considered a distinct bin. Thus, through Fig.  7a, b the outliers are evidenced and eliminated. Then the 95th percentiles of the PGAs within each bin of the two datasets have been calculated. The mean value of the preceding percentiles has been considered as the representative constant value of the first 30 km of the hypocentral distance: 0.27 g for 5 ≤ M w1 < 5.5 and 0.37 g for 5.5 ≤ M w2 ≤ 6.5. It is worthy to be noted that this proposal is related to the PGA values at the rigid ground to be used in microzonation studies at the sites located in the Central Italy Apennine sector within the first 30 km hypocentral distance and in the two ranges of magnitudes of moderately high earthquakes. The disaggregation pairs at each site within NFAs can be determined according to the Ingv study issued at the website: esse1. mi.ingv.it. Then, to select the reference PGA can be used the abovementioned method and the two values found by this study. Further studies must be accomplished to Fig. 6 Datasets of PGA values recorded at NFAs in the Central Italy Apennine Sector from 1997 to 2017 divided into two ranges of M w : a 5 ≤ M w1 < 5.5; b 5.5 ≤ M w2 ≤ 6.5 characterize the PGA of different seismic regions within Italian territory. The same proposed approach or several other proposals can be conceived and applied worldwide within the NFAs taking into account that the surficial soil response there is not dependent on the distance from the source but it is much more dependent on the non-linearity of the soil response combined with the complex geological conditions that cannot be easily modelled.

HVSR measurements addressed in the Aterno Valley
The seismic characterization of surface geology by means of microtremors was introduced by ICMS (2008) and it was then applied in the aftermath of the 2009 L'Aquila earthquake. Many research groups started to record ambient vibrations and process them through Nakamura's method ignoring, in details, the surface geology of each testing point. At Villa Sant'Angelo and Tussillo sites (falling into the Macroarea 6 of the Aterno Valley named L'Aquila crater), we performed several microtremors acquisitions at one station through two devices: Tromino and DAQLink III. The following acquisition parameters have been used: (1) time windows longer than 30′; (2) the sampling frequency higher than 125 Hz; (3) the sampling time lower than 8 ms. Furthermore, Fast Fourier Transform FFT has been used to calculate the ratio H/V. Finally, the spectral smoothing has been performed by means of the Konno-Ohmachi smoothing window. The HVSR values have been calculated for each sub-windows of 20s, then the mean and the standard deviation of all ratios have been calculated and plotted. Further details on the technical aspects of the acquisitions by both devices can be found in Vessia et al. (2016). Figure 8a shows the HVSR measurements acquired at two neighboring points in Tussillo center, where geological characters were similar, by two research groups: T1 (the writing authors) and M5 (the Italian Department of Civil Protection DPC). These acquisitions have been done by the Tromino equipment. As can be noted, the two plots are different: T1 evidences peaks at 2.5 Hz and 8 Hz; on the contrary, the M5 shows the main peak at 2 Hz and minor peaks at 10-20 Hz, 40 Hz and 55 Hz. In the presence of these peaks, the operator would select the most representative one: of course, this selection is highly subjective although the SESAME rules suggest to take into account the highest peaks, such as 2 Hz in both cases (T1 and M5) and disregard the peaks higher than 20 Hz.
On the contrary, Fig. 8b shows the HVSRs measured at two nearby points on a different type of ground type compared with the previous points: T5 (the writing authors) and S4 (DPC group).
In this latter case, the two plots show an evident peak at 2 Hz although the peak amplitude is double at T5 with respect to S4. This difference could be due to the presence of disregarded Love waves that do not have vertical components contributing to the amplification of the horizontal components. Figure 9a, b compare HVSR measured at NE of T5, at Villa Sant'Angelo historical center. In this case, the signals are recorded by two devices used by us: the Tromino and the DAQLink devices. As can be seen, they show similar peaks although no unique peak values can be drawn from each HVSR. In this case, the operator choices can affect the results in terms of the natural frequency of the site. However, the SESAME rule of three acquisitions at 3 different points to assess HVSR could be useful to get to a robust assessment of f o .
Another weak point in the calculation of the amplified frequency f 0 is the systematic differences in calculated amplified frequencies coming from the noise measurements, that induce very small deformation in soil deposits and f 0 drawn from the weak motion tails of the strong Fig. 7 Box and whistler plots of the two datasets of PGA values recorded at NFAs in the Central Italy Apennine Sector from 1997 to 2017 divided into two ranges of M w : a 5 ≤ M w1 < 5.5; b 5.5 ≤ M w1 ≤ 6.5. The bins of hypocentral distance are: a) 1(0-9.95 km), 2(10-19.95 km), 3(20-29.95 km); b) 1(10-19.95 km), 2(20-29.95 km). The void circles are the outliers identified by the box and whistler method motion signals generated by the strong motion events at the site that cause medium to large deformation levels in the ground. From our experience, the f 0 drawn from noise measurements are rarely confirmed by amplified frequencies from actual records. From the field experience, the amplified frequencies f 0 have been measured through the HVSR function from the noise tracks acquired at the Tussillo site, at the point T1 and the weak motion tails of a seismic event recorded on July 7, 2009, at 10:15 local time at the same site. Figure 10 shows the HVSR functions. It is easy to notice that the calculated f 0 related to the noise (Fig. 10a) Fig. 8 a H/V measurements at two neighbor points at Tussillo center: T1 (this study) and M5 (DPC). b H/V measurements performed by the Tromino at T5 (this study) and at S4 (DPC) under similar ground conditions is shown at 8 Hz whereas and the one related to the weak motion (Fig. 10b) is calculated at 2 Hz: the peak frequencies do not match and the weak tail, after the strong motion excitation shows a lower amplification frequency due to the non-linear response of the soil compared to the peak related to the noise measures. These results have been confirmed by other comparisons accomplished in several other places within the Aterno Valley .
Finally, from the abovementioned experiences, three issues can be pointed out: (1) Nakamura's method often provides more than one peak corresponding to different natural frequencies; (2) the peaks are heavily affected by many external factors, especially in urban areas, that are not easy to be disregarded by filtering the measurements; (3) the peaks in HVSR functions are not commonly related to both weak and strong motion amplified frequencies.
Thus, the use of the noise measurements in microzoning activities to derive the bedrock depth should be discouraged especially when the geological conditions of the site are not known, such as the shear wave velocity profile of the soil deposits up to the bedrock depth. In addition, the amplified frequency of the site should be determined through more than one measurement, according to the SESAME rules, in order to check the possible differences induced by the different time of the day and weather conditions at the site. However, the amplified frequency measured at a very low deformation level is modified at medium and large deformations induced during the strong and even weak motion seismic events. Thus the calculation of the f 0 from noise measures can only be used to determine the bedrock depth through Eq. (3) but the shear wave velocity profile is needed as well as the buried geological conditions to guarantee the applicability of the Nakamura method.

Conclusions
In this paper two weak points in microzoning studies have been discussed starting from the authors' experience in microzonation in Italy, that is: (a) the lack of predictivity of GMPEs of PGA measurements in NFAs and (b) the ability of noise measurements to capture the amplified frequency at site even in a complex geological conditions. Throughout the paper, some past experiences of microzoning activity by the present authors are discussed and two proposals have neem put forward. On one hand, concerning the GMPEs of PGAs according to the reference seismic hazard assessment performed in Italy, the need for specific GMPE values in NFAs have been highlighted by several scientists. Here, the proposal of using the 95 percentile of the scattered values recorded within the first 30 km from the hypocentral distance has been provided for the Central Italy Apennine Sector. These values have been drawn from the ITACA database limited to seismic events ranging from 5 to 6.5 M w occurred from 1997 to 2017. On the other hand, after a large experience gained in the noise measurements recorded in the Aterno Valey after the 2009 L'Aquila earthquake, it can be concluding that the noise measurements are not inherently repeatable, thus at least two or three measurements, according to the SESAME guidelines (this represents the European standard to perform ambient vibration measurements) must be requested to calculate the f 0 by means of Nakamura's method. Although noise measurements can provide relevant differences in amplified frequencies according to the operator or the weather conditions, following the SESAME rules guarantee both technical standards to an unregulated geophysical technique that relies on a simple buried geological model of horizontally layered soil deposits. This model, when not applicable, can make the HVSR function from noise measurements totally misleading. Thus, even at level 1 microzoning studies, the use of direct and indirect measures is needed in order to confirm the layered planar setting of the subsurface geolithological model and measuring shear wave velocity profiles to enable a robust prediction of the bedrock depth by means of the amplified frequency f 0 estimation. Finally, when the noise measurements are compared with the weak motion tails of actual seismic events, they show different amplified ranges of frequencies. Accordingly, it must be kept in mind that soil behavior is strain-dependent: this means that their natural frequencies at small strain levels (microtremors) differ from the ones at medium strain level (weak motions) and at high strain level (strong motions). Then, depending on the purpose of the natural frequency measurement, different strain levels will be investigated to do an adequate characterization of the site response under seismic excitation.