Shear modulus reduction and damping ratios curves joined with engineering geological units in Italy

Numerical simulations of seismic site response require the characterization of the nonlinear behaviour of shallow subsoil. When extensive evaluations are of concern, as in the case of seismic microzonation studies, funding problems prevent from the systematic use of laboratory tests to provide detailed evaluations. For this purpose, 485 shear modulus reduction, G\G0(γ) and damping ratio, D(γ) curves were collected from multiple literature sources available in Italy. Each curve was associated with the related engineering geological units considered in seismic microzonation studies. A statistical analysis of the data was carried out with the aim of shedding light on the significant difference between the laboratory classification of samples and the macroscopic/engineering geological one, provided during seismic microzonation studies. Since the engineering geological classification plays a prominent role in extensive site response evaluations, the outcomes of the present work may be of help at least when preliminary seismic response estimates are of concern. The dataset provides reference information that can serve as key data for large-scale hazard assessments worldwide.


Background & Summary
Simulations of waves propagation are a recurrent practice for the quantification of the ground motion expected at a site.Since they are performed considering the uncertainties of one or more parameters that play a role in seismic hazard, they may be included among the probabilistic seismic hazard applications (PSHA).As reported in the existing literature, PSHA contemplate different levels of increasing complexity 1 .In level 0, a fully probabilistic seismic hazard estimate requires the convolution of the hazard defined on bedrock virtually ignoring site-specific information [2][3][4] , and of amplification factors derived from a Ground Motion Prediction Equations logic tree 5 .In the subsequent levels, fully probabilistic hazard estimates are retrieved adopting a site-specific approach, introducing single-station standard deviation values and validations at seismological networks 6 , or the convolution of the hazard curve on rock with the probability distribution of the amplification functions obtained from analytical soil response analyses.In this latter case, estimates on the variability of ground motion, S-waves velocities Vs and nonlinear properties associated with the local seismo-stratigraphical configuration are necessary 7 .
It is now worth noting that the geographical scale of the PSHA fits well with the scale of the seismic microzonation (SM) studies.Therefore, information obtained at large (local) or even regional scale (<1:25,000) through SM studies may be especially useful for PSHA.SM is a practice that is able to account for site effects at different levels of detail (known as 1st, 2nd or 3rd level), highlighting on 1:10,000 scale or more detailed specific maps the areas most prone to seismic hazard 8 .It is commonly accepted that SM is a tool of fundamental importance for land use, planning and to maintain engineering infrastructures 9,10 , although SM should be considered only the first step towards a comprehensive seismic hazard assessment of the total site-specific hazard.Starting from the results of extensive numerical modelling performed over 138 municipalities in central Italy 11 , the attempt to represent absolute estimates of the seismic hazard has already been followed by Mori et al. 12 .Also, Barani et al. 13 provided another example in Italy of incorporating the results of 2nd level SM studies into probabilistic seismic hazard analysis.
Still further analyses should be performed, however, and advanced methods developed to extend the results over wider areas.This can be achieved using already available data such as the significant amount of morphological and/or geological-geotechnical data, which together may provide the possibility of defining appropriate proxies (Vs30) that are suitable to catch the site amplification, with a partially probabilistic/hybrid approach 14 .It should be noted that, in common practice, due to the fact that increasing the SM level is a function of the available economic resources, SM studies are often only qualitative (i.e. if the SM is conducted at a 1st level) or, if they contain quantitative estimates (i.e. in terms of amplification factors, if the SM is conducted at a 2nd or 3rd level), they do not provide synthetic parameters that express the local and reference hazard together computed probabilistically.Only at the building/structure scale, several analyses have been carried out by one-dimensional Monte Carlo ground response analyses using a fully probabilistic approach [15][16][17] .
Thus, the rationale behind this study rests on the proposition of a robust harmonization of the data available from SM studies to support the need for fully probabilistic seismic site-specific hazard studies, of which still today no example exists.Otherwise, it should be considered that in Italy up to 2,000 1st level SM studies are available (https://www.webms.it/servizi/stats.php) and the integration of the subsoil engineering geological model 18 with with the association with Vs has already been attempted by Romagnoli et al. 19 .
One of the basic products of a SM study is an engineering geological map and the conceptual interpretation of the subsoil under investigation in terms of engineering geological units.This kind of information differs from that contained in a basic geological map because the latter does not represent the dynamic nature of the subsoil 20 .
Starting from the engineering geological setting, it is possible to aggregate samples used for laboratory tests to the same engineering geological unit: a dataset associated with this kind of engineering geological unit classification constitutes an untried method in existing research.
The paper is organized in three main parts.After a general description of the data and methods, statistics are shown and, finally, comparisons with literature data are discussed in detail.

Methods
As the first step, shear modulus reduction, G\G 0 (γ) and damping ratio, D(γ) curves were collected and associated with the related engineering geological units considered in SM studies (eg-units hereafter).The goal of the engineering-geological classification proposed for Italian SM studies is to group together soils and rocks in two main categories, the "Cover terrains" and "Geological bedrock" units respectively, considering their geological and geotechnical properties or attributes in order to analyse the seismic local effects at urban municipality scale 21 .The Cover terrains units collectively represent all kinds of loose, incoherent and unconsolidated superficial deposits such as, gravel, sand, clay, and organic material that originated generally in the Quaternary era.These include slope deposits, terraced and recent alluvial deposits, terraced marine deposits, polygenic detritic coverings, ancient glacial deposits, lacustrine sediments, eluvial-colluvial and landslide deposits.The cover units are classified, according to the Unified Soil Classification System 22,23 (Table 1), into coarse grain and fine-grained soils.The coarse grain soils consist of gravels (G) and sands (S).Each class is further subdivided into four units depending upon the grading and inclusion of other grain sized materials, combining the "G" or "S" acronyms with "W" for well graded, "P" for poorly graded, "M" for containing fine materials and "C" for clay binder.The coarse grain soils also contain anthropic deposits (RI).The fine-grained soils include silts and clays and are divided into three classes named with the acronyms "M" for inorganic silts and very fine sands, "C" for inorganic clays and "O" for organic silts and clays.These are combined, on the basis of their liquid limit and plasticity index, with the acronyms "L" for low plasticity and "H" for middle and high plasticity.The fine-grained soils also contain peat and other highly organic soils (PT).The designation of unconsolidated units refers to the dominant grain size of clastic material mixtures of different sizes.The Cover terrains are thus classified in 16 eg-units (Table 1).The Geological bedrock units consist of lithoid and consolidated deposits of geological formations, comprising weathered and fractured portions, classified following lithostratigraphic criteria, structural features and facies 19,24,25 .Examples of Geological bedrock units are limestones, sandstones, siltstone, dolomites, chert, marly calcareous and marly bedrock, pelitic and arenaceous bedrock, brecciated and conglomeratic bedrock.This category also includes 16 e-g units starting from 4 main types of rocks: lapideous rocks "LP" (e.g.limestone, dolomites), grainy cemented rocks "GR" (e.g.sandstones, conglomerates), cohesive over-consolidated rocks "CO" (e.g.over-consolidated clays) and deposits characterized by alternations of the contrasting lithotypes "AL" (e.g.flysch deposits).All the other units derive from these four main units.If they are stratified, the acronym "S" is added to form another 4 units ("LPS", "GRS", "COS" and "ALS").If the previous 8 eg-units are fractured and/ or weathered, the prefix "SF" is added to the beginning of the acronym (e.g."SFLP", "SFGRS", "SFALS", "SFCO"; Table 1).The Italian SM classification of the geological bedrock units considers the intact, stratified and weathered and/or fractured rock properties 21 .
All the G\G 0 (γ) and D(γ) curves were singularly regularized according to the procedure proposed by Yokota et al. 26 .This latter allows the relationships between G\G 0 and the strain γ to be found as the simplified formula in Eq. ( 1), and D as the simplified formula in Eq. ( 2).Substantially, G\G 0 and D values are defined by means of the three constants: namely, α and β for G\G 0 and λ for D: The constants α, β and λ are obtained through a double-step procedure of adaptation of the experimental data to the analytical linearized expressions of Eqs.(1, 2): firstly the parameters of Eq. ( 1) are obtained and then they are used to calibrate λ.
Curves were regularized up to a value of the γ level of 0.0001%.The raw data was archived As-Is in the original pdf file, while curves were each individually regularized using a codified procedure.Before defining parameters of the Yokota modeling adaptation, points recognized outliers were manually deleted.A further and parallel regularization was also performed considering all the curves for each eg-unit with the aim of representing the behaviour of soils in macroscopic terms.

Data Records
The collection, available at the link: https://doi.org/10.5281/zenodo.8134979 27, has been carried out nationwide considering the available data from SM studies, public databases and published works, according to European Commission principles 28 .
Primarily, in order to have a Findable, Accessible, Interoperable, Re-usable (so called FAIR) dataset (developed according to the European Open Science Cloud -EOSC policies, https://ec.europa.eu/research/openscience/index.cfm?pg=open-science-cloud), each set of laboratory test results has been saved in a standard file.The dataset will be useful for accomplishing the following purposes: to access and interoperate research data throughout web-accessible services (for instance, by means of the webpage: https://www.webms.it/servizi/catalog.php) and guarantee public access to subsoil information in the perspective of data integration with already existing web-based databases (i.e., European Geotechnical Database service -EGD, http://egd-epos.civil.auth.gr/; the New Zealand Geotechnical Database -NZGD, https://www.nzgd.org.nz/).This approach will ensure the continuous maintenance of the dataset, which can be updated every time new information is available.
Each file in the dataset is described in its accompanying metadata file, which can be seen as a complementary footnotes sheet.The metadata file contains the following information: • Rootfilename: basename of the raw archived file; • Macroarea: name of the macroarea (region or area) of the SM study where the sample was collected; • Municipality: the municipality where the sample was taken; • Type of laboratory test; • Depth top and bot (m): depth of the top and bottom sampling computed from the surface level; • γ (kN\m 3 ): unit weight; • WL (%): water content at the liquid limit; • PI (%): Plasticity Index; • USCS code: code according to the USCS classification; • eg-unit SM: code of the eg-unit retrieved from SM study; • X and Y coordinates from a WGS84/UTM-33N datum; • Ref: link or references to the data source.
• Namely, the dataset consists of 485 G\G 0 (γ) and D(γ) curves obtained from: • the third level of the SM studies carried out following the 2016-2017 Central Italy seismic sequence (https://sis-ma2016data.it/), in which the dynamic behaviour of silty and clayey soils was first studied by Ciancimino et al. 29 ; • SM studies carried out in the Emilia Romagna region 30 ; • SM studies of Roma Palatino 16,31,32   • SM studies carried out on the eastern flank of the Mount Etna volcano following the 2002 Santa Venerina earthquake (Protezione Civile Catania Working Group 35    The geographic location of the sites where laboratory tests were collected is shown in Fig. 1. We analysed the experimental G\G 0 (γ) and D(γ) curves obtained from different types of geotechnical laboratory tests: Double Specimen Direct Simple Shear, DSDSS; Resonant Column, RC; Cyclic Triaxial, TXC; Cyclic Torsional Test, CT; Cyclic Torsional Shearing, CTS; Resonant Column and Cyclic Torsional Test, RCT.In several sites, for each sampled layer, different laboratory tests were performed to enlarge the range of deformations analysed.In these cases, the results are reported in the dataset in separate rows.Figure 2 graphically visualises the similarities among available samples: samples were taken mostly in unconsolidated clastic deposits of cover terrain units, although about 10% of the tests were carried out for geological bedrock units.Fig. 3 Violin plots for the depth of sampling for each eg-unit 53,54 .Samples were taken mostly at depths ranging from 1 m to 70 m below ground level (Fig. 3).The initial (small strain) values of the damping ratio, D 0 , range between 0.2% and 63%.The smallest and the largest shear strain values obtained from laboratory tests are 1.0 × 10 −5 % and 5.2 × 10 −1 %, respectively.The Fig. 7 G\G 0 (γ) curves (a) and D(γ) curves (b), adapted from Yokota et al. 26 , for α CLSM , α MLSM and α SMSM egunits.Each unit is represented by the mean and by the Darendeli confidence levels (±95%) 37 .dataset also includes some samples of a few organic clays with very high water content and low unit weight.The structure of the data array is depicted in Fig. 4.Moreover, each case history includes identifying information (e.g.ID, geographic coordinates) in a metadata file: the compiled post-processed data is presented in a single file suitably archived.

technical Validation
All the data was associated with an engineering geological unit and, if possible, with the USCS classification.At this point, the discussion deserves a focus on the representativeness of the samples, since differences exist between the USCS units obtained from the laboratory certificate, and the eg-units in the dataset.This constitutes a crucial aspect intrinsic to the process of extending results that are available for a few centimetres (i.e. the dimensions of the samples) to meters (i.e. the layer thickness) and from a few verticals (i.e.boreholes) to larger areas (i.e.cross-sections and seismically homogeneous microzones, SHM 8 ).Statistical analysis was performed to investigate the correspondence between the USCS and eg-unit in all those cases where the two classifications are available for at least 7 samples.Figure 5 shows the distribution of the USCS codes among each seismic microzonation code as vectors going from the origin of the plot to the percentage of data availability (each angle direction is represented by a different USCS code).
No code attributed in the seismic microzonation coincides with the USCS code for more than 50% of the population.From Fig. 5, it can be seen that for CL SM only 46% of the samples (44 of a population of 95 samples) are described by the same acronym.Meanwhile, for ML SM only 12% of the samples (14 out of 115 samples) are described by the same acronym, while for SM SM only 15.4% of the samples (8 of 52 samples).The heterogeneity and anisotropy in the materials and geological formations seem more marked in the cases of SC SM and GM SM than in the cases of CL SM , ML SM and SM SM .This behaviour may be due to the nature of the materials that constitute the specimens and to difficulties in the sampling operations: in all these cases, the specimens contain some finer levels of the main coarse deposits.In general, three concomitant aspects should be considered as sources of bias: 1) heterogeneity and anisotropy in the materials and geological formations; 2) unavailability of samples according to regular meshes of investigation, due to the cost of a theoretical massive-invasive exploitation; 3) subjectivity in the visual inspections of the sample.This latter may induce different attributions of the code attributed in SM.In the case of the finer materials, the differences between the two classifications were extrapolated also on a Casagrande chart for the most populated eg-unit classes (Fig. 6).Contextually, the correlation coefficients for CL SM and ML SM were also computed.The values (equal to 0.86 and 0.85, respectively) denote a higher variability in the case of ML SM than CL SM for the two variables: water content, WL and Plasticity Index, IP.
As a further step, a representation of these in terms of median formulations of the G\G 0 (γ) and D(γ) curves were obtained (Fig. 7).The laws of variation of G\G 0 (γ) and D(γ) curves for each eg-unit were determined through the formulation of Darendeli 37 , which describes the standard deviation for the normalized modulus reduction and the damping curves in the form of equations based on statistically retrieved parameters.Mean and mean ± standard deviation curves for each unit were also made available in the archive in the folder "average".The standard deviations have the form indicated in Eqs. 3, 4, respectively, for G\G 0 (γ) and D(γ): ( ) In Table 2 the maximum and minimum values of the σ G\G0 and σ D with the strain are reported for each eg-group.The values quantify the apparent aleatory randomness of the G\G 0 and D values with a confidence level of 95%.
Comparing the curves in terms of mean values, only significant differences for the SM SM curves may be distinguished with respect to the ML SM and CL SM curves at low γ values (lower than 0.03%): the CL SM and ML SM eg-units follow approximately the same behaviour both for G\G 0 and for D.
We tested the null hypothesis of the pairwise difference between data vectors of the α and β parameters used to smooth the curves according to Yokota et al. 26 .The results are synthetized in Table 3.At the 5% significance level, the returned value of h = 1 for α CLSM vs α SM SM indicates that the t-Student test rejects the null hypothesis, and thus suggests the presence of significant differences between the two populations.Symbol p in Table 2 is the probability of observing a test statistic to be as extreme as, or more extreme than, the observed value under the null hypothesis.For α MLSM and α SM SM the returned value of p is equal to 0.9; otherwise, for α CLSM and α MLSM , p = 0.13.No significant differences are identified among the three populations of β used for D regularization.
As stated in Wasserstein et al. 38 , conclusions should not be based solely on whether an association was found to be statistically significant.According to this consideration, the most commonly used in numerical modelling curves were investigated.The variation with the Vucetic and Dobry 39 and Darendeli and Stokoe 40 models were simulated for a Plasticity Index ranging from 15 to 50% (Fig. 8).The highest differences of the means curves for CL SM from those in the literature may be recorded at very high strain levels (0.3 and 0.8-0.9%respectively for G\G 0 and D).Generally, the behaviour of the seismic microzonation curves ± the standard deviations is able to include the predicted variability based on the Plasticity Index recorded in the dataset.Despite this evidence, the SM SM curves show the highest standard deviations compared to the literature data.As a remark, the rationale of this study is that by extending seismic microzonation data, it is possible to account for uncertainty in a coherent framework, where subsurface geometries and buried morphologies also have a similar amount of uncertainty.As a result, the predictions made by this study are larger than those in the literature are.

Usage Notes
The curves shown in this work and identified using the laboratory data of the seismic microzonation studies can be adopted as input in 1D calculation codes to carry out local seismic response studies, as shown in Fig. 9.
The results shown before suggest that a further merge of the eg-units is possible.This was previously confirmed also in terms of the S-waves velocity, Vs by Romagnoli et al. 19 .In practice, from the point of view of the non-linear behaviour of soils, a macro-group of eg-units may be constructed including all the eg-units relating to clays and inorganic silts (ML, CL, MH and CH) in one single macro-group, while two other macro-groups may be defined for: 2. OH + OL and 3. SM + SC + SP + SW.It is worth noticing that all the curves defined for each macro-group may be adopted only to reproduce the response of soils located whitin the first 15 m.The laws of variations of G\G 0 and D have the forms indicated by Eqs. 5, 6, respectively, and the parameters of Table 4.
The G\G 0 and D curves may be described using the aggregated variation laws defined ad hoc for seismic microzonation through the parameters reported in Table 4.For completeness, Table 5 reports also the maximum and minimum values of the σ for G\G 0 and D with the strain for the three previously introduced macro-groups.
Thus, the parameters of the hyperbolic model for eg-unit groups and macro groups, respectively, are shown in Tables 2, 5. Neverthless, the outcomes of this study can be used in any code that simulates 1D propagating waves by using the parameters provided in Table 4 and the formulation in Eqs.(5, 6), when Darendeli's model is not implemented.
The present work fits in the field of fully probabilistic seismic hazard assessment.The level at which these results feature in the entire process is indicated in Fig. 9 with a red asterisk.
It is outside the scope of the work to suggest a correlation model that examines the connection between the variation in G\G 0 reduction and the variation in D increase 41,42 , but it may be a topic for future research.
The dataset may be used to adapt models from the laboratory to the regional/local scale, similarly to what happens for analogous models in the laboratory to real-scale models 43 .In other words, the eg-unit definition allows the modelling of the dynamical properties of a geological body when changes of scale are applied.This scaling operation is then even more important considering that at least four sources of uncertainties may affect the numerical modelling results when using laboratory tests data as inputs: 1) loading directionality; 2) simplified schemes of application of the cyclic loading; 3) drainage conditions and 4) representativeness of samples.The present lack of knowledge of engineering geology at a regional scale has until now limited the interpretation of the available data.Thus, these results aim to provide new insights about this topic, consequentially looking at seismic prevention at a regional scale, rather than at single municipality scale.This scale is even more important since agglomerates of adjacent hamlets strictly interact with each other.This approach is one of the principles at the base of the new Italian code of Civil Protection 44 .
Moreover, this study illustrates relevant information in the perspective of performing 3D numerical modelling at a local/regional scale, which was described as one of the grand challenges by Forsyth et al. 45 .3D numerical modelling is recently being even more widely diffused and adopted because of its ability to explain the complex pattern of strong ground motions after or before an earthquake event, but nowadays only linear simulations are performed due to the computational cost and lack of data.Therefore, starting from the average values of the defined curves of this study, future simulations could be run where eg-unit models are available 46,47 .
The cascade effect resulting from this analysis can also provide new data suitable for achieving a detailed physical understanding of the nonlinear processes of waves propagation after events causing damage.It is generally assumed as a rule of thumb that the damping ratio D may be related to the a-dimensional Q factor using the expression: D = 1\ (2Q) 48 .On this subject, Dimitriu et al. 49 and Lacave-Lachet et al. 50also showed that an important contribution to κ (the seismological measure of wave attenuation) is the inelastic attenuation (D) in the site's subsurface geology.The data described in this study can provide further information suitable for comparisons with seismological data 51 and it consequentially has the potential to be used for several purposes (i.e., stochastic ground-motion prediction calculations; nonlinearity and attenuation of seismic waves relationships).
Moreover, the present dataset may allow the relaxing of the ergodicity hypothesis on the nonlinearity among all the parameters, which regulates the seismic response which also includes stratigraphy, shear wave velocities and Vs.

Fig. 1
Fig. 1 Location of the investigated sites.

Fig. 2
Fig. 2 Dendrogram of the number of available samples.Each node represents a single eg-unit.The size of the nodes is proportional to the number of available samples.

Fig. 5
Fig. 5 Vector distribution among USCS classes for the eg-units CL SM (a), ML SM (b), GM SM (c), SM SM (d), SC SM (e), and SW SM (f).Units on the vector axes are synched and expressed in %.

Fig. 6
Fig.6 Scatter plot of the classification of the finer soils according to the Casagrande chart.SM subscript in legend stands for "seismic microzonation perspective".

Fig. 8
Fig. 8 Comparisons with existing literature curves.G\G 0 (γ) and D(γ) curves for CL SM eg-unit (a and b, c and d) ± Darendeli conference levels (adapted from Yokota et al. 26 ), for a mean confining effective pressure σ' of about 200 kPa, compared respectively with Vucetic and Dobry 39 and Darendeli and Stokoe 40 , confining effective pressure σ' = 200 kPa; G\G 0 (γ) and D(γ) curves for SM SM eg-unit (e and f) ± Darendeli conference levels (adapted from Yokota et al. 26 ), for a mean confining effective pressure σ' of 180 kPa, compared with Seed and Idriss 56 curves -mean, upper and lower bound.

Fig. 9
Fig. 9 Flowchart of the overall process devoted to the probabilistic hazard assessment.Red asterisk indicates the present research positioning.

Table 1 .
21gineering-geological classification adopted in SM studies by following Italian standards21 ;• SM studies of Nocera Umbra 33 ;• SM studies carried out following the 2009 L' Aquila earthquake (MS AQ Working Group 34 );

Table 2 .
Minimum and maximum σ values for each eg-group.