Local Site Effects Research Articles

Variation in the contribution of macroinvertebrates to wood decomposition as it progresses

AbstractAlthough necromass decay rates are limited by the slowest portions to decompose, most decomposition studies examine only the earliest stage of decay. As such, these studies run the risk of yielding misleading results regarding the relative contributions of different decomposers. For example, the contributions of macroinvertebrates to wood decomposition remain mostly unknown beyond the first 50% of mass lost, despite drastic changes in substrate conditions over time. We sought to clarify how the macroinvertebrate contribution to decay changes over the course of wood decomposition in the Southeastern United States—a region with a long history of wood decomposition research. To this end, we (1) compiled data from published studies comparing wood decay with and without macroinvertebrates; and (2) conducted a field study assessing wood mass loss, with and without macroinvertebrate access, at three sites across the region over four years. With these combined data, we analyzed macroinvertebrate contribution as decay progressed, revealing a quadratic relationship, wherein macroinvertebrate contribution increased early in decomposition and then began to decline as decay progressed. Strong local site effects, particularly the abundance and activity of termites, determine the time required for wood to reach this point of mass loss.

Local site effects and seismic microzonation around Suban Area, Curup Rejang Lebong, Bengkulu deduced by ambient noise measurements

BackgroundThe Suban area of Curup Rejang Lebong is a tourist region in Bengkulu Province, Indonesia, close to the active Ketaun and Musi faults, which are segments of the Sumatra Fault System (SFS). However, no studies have been conducted in this area to assess how geological structures affect seismic ground motions and contribute to seismic hazard and risk assessment.MethodsThe first study of seismic microzonation in the Suban area of Curup City by ambient noise measurements was conducted at 100 sites, spaced ~ 1 km apart, with 60 min of data acquisition for each site. All microseismic data were processed using the Horizontal to Vertical Spectral Ratios (HVSR) method.ResultsThe HVSR method revealed the amplification factors (A0) ranging from 1.23 to 8.26 times, corresponding to natural frequency (f0) variations between 1.24 and 9.67 Hz. About 13% and 55% of the sites show high (6 ≤ A0 ≤ 9) and medium (3 ≤ A0 ≤ 6) amplifications, respectively, predominantly in the western parts of the study area, consistent with a high seismic vulnerability index (Kg). Furthermore, we also estimated the ground shear strain (GSS) of the region using the Kanai method with two large historical earthquakes at the Ketahun segment in 1943 (Mw 7.4) and the Musi segment in 1979 (Mw 6.0). The Kg value is consistent with the GSS values and indicates areas of severe damage during the historic earthquakes.ConclusionsThus, the western parts of the Suban region are vulnerable to severe damage from an earthquake. These findings could provide valuable insights for future planning and risk management efforts aimed at minimizing the impact of earthquakes in the Suban region.

Unraveling the Seismic Source in Archaeoseismology: A Combined Approach on Local Site Effects and Geochemical Data Integration

Archaeoseismological research often deals with two unresolved questions: the magnitude and level of damage caused by past earthquakes, and the precise location of the seismic source. We propose a comprehensive review of an integrated approach that combines site effects with the analysis of geochemical data in the field of archaeoseismology. This approach aims to identify active buried faults potentially related to the causative seismic source and provide insights into earthquake parameters. For each integrated method, we report the foundational principles, delineation of theoretical field procedures, and exemplification through two case studies. Site effects analysis in archaeoseismology assumes a pivotal role in unraveling historical seismic occurrences. It enables estimating the earthquake magnitude, assessing the seismotectonic patterns, and determining the resulting damage level. Valuable data related to earthquake parameters can be extracted by analyzing vibration frequencies and acceleration measurements from structures within archaeological sites. This information is instrumental in characterizing seismic events, evaluating their impact on ancient structures, and enhancing our understanding of earthquake hazards within the archaeological context. Geochemical investigations supply indispensable tools for identifying buried active faults. The analysis of fluids and gases vented in proximity to faults yields valuable insights into their nature, activity, and underlying mechanisms. Faults often manifest distinctive geochemical imprints, enabling the differentiation between tectonically active and volcanically related fault systems. The presence of specific gases can further serve as indicators of the environmental conditions surrounding these fault networks. Integrating site effects analysis and geochemical investigations within archaeoseismological research is crucial to improving our understanding of unknown past earthquakes. Moreover, it enhances the seismic hazard assessment of the region under study.

The 2023 Turkey–Syria Earthquake Sequence: Ground-Motion and Local Site-Effect Analyses for Kahramanmaras City

The 2023 Turkey–Syria Earthquake Sequence: Ground-Motion and Local Site-Effect Analyses for Kahramanmaras City

Quaternary Deposit Response to Earthquakes in Pemalang City Based on Peak Ground Acceleration, Earthquake Intensity, and Microtremor Method

The northern part of Pemalang City consists of Quaternary deposits, having the potential for earthquake amplification effect. This amplification effect amplifies the ground shaking because of an earthquake (the local site effect) that has the potential to cause damage. This study investigated the amplification factor from the HVSR curve of microtremor measurements due to soil response based on ground shear strain, the risk level of the earthquake based on peak ground acceleration (PGA), and earthquake intensity. The microtremor data from five locations in Pemalang were used to calculate the amplification factor and predominant frequency. The damaging earthquake parameters around Java during 2010-2020 were used to calculate the PGA. The microtremor data were processed using the HVSR method, and PGA was calculated using the Kanai equation. The HVSR result shows that Pemalang has an amplification factor ranging from 6.23 to 19.59 and ground shear strain varying between 0.86 x 10-4 and 6.67 x 10-4, which shows that Pemalang only experiences the vibration when an earthquake occurs. The PGA results using the Kanai equation (19.71-54.56 gal) were included in the low vulnerability category, and MMI earthquake intensity (3.08-4.70) was included in the felt earthquake category (II SIG BMKG scale). Therefore, the amplification factor from the HVSR curve of microtremor measurement, ranging from 6.23 to 19.59, showed low soil response and low-risk vulnerability based on the damaging earthquake parameter around Java during 2010-2020. Keywords: peak ground acceleration, amplification, microtremor, PGA

Stochastic simulation of earthquake ground motions based on improved finite-fault model

In the seismic reliability analysis of engineering structures, it is often necessary to adopt stochastic ground motions as inputs. However, the current methods used to generate stochastic samples, such as the spectral representation method, the spectral transfer method, the dynamic source method, the physical source method, and the comprehensive method, etc., are not very convincing since they do not contain the complete information of seismic occurrence and propagation. In this study, an improved finite-fault model for stochastic simulation of ground motions is proposed. Compared to existing stochastic finite-fault model, the improved model considers the propagation law of phase spectrum, local site effect, and the impact of wave velocity on the arrival time of seismic waves. Besides, identification and statistics of model parameters are carried out based on ground motion records. Reasonable samples of stochastic ground motions are then obtained. In addition, parameter sensitivity analysis is conducted to determine critical random variables, which reduces the complexity of stochastic simulation. Simulation results are compared with recorded accelerograms from the 1999 Chi-Chi earthquake. It is shown that the improved stochastic finite-fault model is reliable, and the acceleration response spectra of simulated samples are consistent with those of station records collected at various site classes and the design spectra.

Full-waveform simulation of DAS records, response and cable-ground coupling

SUMMARY Over the past several years, the use of optical fibre cables as ground motion sensors has become a central topic for seismologists, with successful applications of distributed acoustic sensing (DAS) in various key fields such as seismic monitoring, structural imaging and source characterization. DAS response is a combination of both instrument response and cable-ground coupling, with the latter having a strong, spatially variable, but yet largely unquantified effect. This limits the application of a large number of staple seismological techniques (e.g. earthquake magnitude estimation, waveform tomography) that can require accurate knowledge of a signal's amplitude and frequency content. Here we present a method for accurately simulating a DAS cable and its ground coupling. The scheme is based on molecular dynamic-like particle-based numerical modelling, allowing the investigation of the effect of varying DAS-ground coupling scenarios. We start by computing the full strain field directly, for each pair of neighbouring particles in the model. We then define a virtual DAS cable, embedded within the model and formed by a single string of interconnected particles. This allows us to control all aspects of the cable-ground coupling and their properties at an effective granular level through changing the bond stiffness and bond types (e.g. non-linearity) for both the cable and the surrounding medium. Arbitrary cable geometries and heterogeneous materials can be accommodated at the desired scale of investigation. We observe that at the meter scale, the cable-ground coupling and local site effects can substantially alter the recorded signal. We find that the stiffness of the thin layer of material to which the cable is coupled has the strongest effects, selectively amplifying portions of the wave train and contributing to substantial phase delays. These differences show that cable coupling and local site effects should be considered both when designing a DAS deployment and analysing its data when either true or along-cable relative amplitudes and/or frequencies are considered. The codes developed herein for calculating full waveform DAS responses and coupling are made publicly available.

Shear-wave velocity determination by combining data from passive and active source field investigations in Kumamoto city, Japan

We present the 1D subsoil structure and local site effects at KUMA strong ground motion station in Kumamoto City, Japan. We analyze data from a field campaign conducted in the framework of the Blind Prediction BP1 test of the 6th IASPEI/IAEE International Symposium: Effects of Surface Geology on Seismic Motion. In parallel with other participants of the BP1 test, we process data from passive and active source measurements aiming to determine the shear-wave velocity, Vs, structure and the site response at KUMA station. Passive measurements are associated to five microtremor arrays. In each array, seven seismometers have been deployed in a common-center triangle shape, recording microtremors simultaneously for 1 to 2 h. The vertical component of microtremors was analyzed using the spatial autocorrelation (SPAC) method. Cross-correlation coefficients were computed for all station pairs available for each array. By fitting the average SPAC’s coefficients to the first-kind zero-order Bessel function, J0, and assuming that microtremors primarily comprise fundamental mode Rayleigh waves, phase velocity dispersion curves were determined. Phase velocity values for frequencies > 15 Hz were obtained from data of a close-by active source geophone profile. We integrated the data with those of the passive measurements and obtained an experimental phase velocity dispersion curve. The resulting curve shows low velocity variation, from 150 to 200 m/s, in the surface layers, whereas significant dispersion appears in frequencies below 2.5 Hz. By inverting this curve, we achieved to determine the 1D shear-wave velocity structure at KUMA station. Site response characteristics were determined by applying the Horizontal-to-Vertical-Spectral-Ratios method. Significantly amplified peaks in the frequency range between 0.3 to 1.5 Hz dominate HVSR spectral ratios. These peaks correspond to resonant frequencies of soils and originate from different impedance contrasts within the substratum of the site.Graphical

A Simplified Approach to Estimate Seismic Vulnerability and Damage Scenarios Including Site Effects. Application to the Historical Centre of Horta, Azores, Portugal

ABSTRACT The research attempts to provide a method for evaluating the susceptibility of heritage buildings and potential damage considering local site effects, using as a case study an urban sector situated in Horta, Faial, Azores, Portugal. The physical vulnerability of the investigated structures was appraised using an index-based method devised for structural masonry aggregates to emphasize their likelihood of seismic damage. As a result, a damage scenario based on real accelerograms collected after the 1998 Azores seismic event was developed by using a reliable seismic intensity prediction equation. In addition, the soil amplification, which increases the extent of seismic damage to the examined structures, is assessed by varying the soil conditions using EN 1998–1 Code. Thus, this research emphasises the ground impact in forecasting the expected damage scenarios, enabling optimised risk mitigation strategies in urban settings.

Vs30 and depth to bedrock estimates from integrating HVSR measurements and geology-slope approach in the Oslo area, Norway

In order to estimate well-constrained seismic hazard and risk on local scales, the knowledge of site amplification factors is one of several important requirements. Seismic hazard studies on national or regional scales generally provide the level of earthquake shaking only at bedrock conditions, thereby avoiding the difficulties that are caused through local site effects. Oftentimes, local site conditions are not well understood or even non-existent. In this study we investigate an efficient and non-invasive methodology to derive the local average shear wave velocity in the uppermost 30 m of the ground (Vs30). The Vs30 value is a useful parameter to define soil classes and soil amplification used in seismic hazard assessment and to extend the knowledge of the site to include the depth to basement rock. At the level of the municipality of Oslo, there is currently no map available that describes the Vs30, and as such any seismic risk study is lacking potentially critical information on local site amplification. The new proposed methodology includes the use of existing well databases (with knowledge on minimum basement depth), topographic slope derived from Digital Elevation Models (as a proxy for both depth to basement and Vs30, integrated with geological maps) and near-surface Quaternary geological maps. The Horizontal to Vertical Spectral Ratio (HVSR) method and a statistics-based geological mapping tool (COHIBA) are used to integrate the various sources of data estimates. Finally, we demonstrate our new methodology and workflow with data from three different regions within the Oslo municipality and propose an approach to conduct cost-efficient mapping for seismic site amplification on a general municipality scale.

Seismic microzonation mapping of Greater Vancouver based on various site classification metrics

The goal of the multi-year seismic microzonation mapping project for Greater Vancouver, British Columbia, Canada, is to produce seismic hazard maps inclusive of local site effects, in particular seismic hazard specific to one-dimensional site response and three-dimensional Georgia sedimentary basin amplification, as well as liquefaction and landslide hazard potential. We explore the variability in key seismic site characterization measures most often used for seismic microzonation mapping to evaluate the impact on mapping and communication of seismic microzonation of Greater Vancouver. This study focuses on the comparison of seismic microzonation maps of Greater Vancouver based on up to three seismic site term parameters and their associated classification schemes: 1) the time-averaged shear-wave velocity (Vs) of the upper 30 m (Vs30) and associated Canadian National Building Code (NBC) site class; 2) Vs30-based site classification proposed for the updated Eurocode 8; 3) site period (T0) determined from microtremor site amplification spectra; and 4) a hybrid site classification based on T0 and the average Vs and thickness of soil. 810 Vs30 and 2,200 T0 values are determined to evaluate sub-regional differences in these important seismic site parameters in Greater Vancouver. We find that the seismic microzonation of Greater Vancouver depends on the chosen seismic site parameter (Vs30, T0, or a combination of parameters) and that classification schemes with greater class divisions are beneficial to communicating the great variability in seismic site conditions in Greater Vancouver. We recommend that either one hybrid classification map or two classification maps of Vs30 and T0 together are required for effective communication of the seismic microzonation of Greater Vancouver.

Insights on the role of local site effects on damage distribution in the Izmir metropolitan area induced by the October 30, 2020 Samos earthquake

On the 30th of October 2020, a 6.6 magnitude earthquake occurred 14 km north of Samos Island, causing 119 casualties (117 in Izmir, Türkiye, and 2 in Samos, Greece) and significant damage in the 3rd biggest city of Türkiye, Izmir. Although the city is roughly 70 km far away from the epicenter, the damage was significant and concentrated in the city center settled on alluviums. This paper aims to analyze the distribution of damage in Izmir province, by crosschecking the recorded motions, the subsoil conditions and the evidence of damage as collected by an ad-hoc on-site reconnaissance. The intrinsic behavior of the Samos earthquake was investigated by employing three different ground-motion prediction equations. The results of the analyses revealed that site effects play a significant role in the amplification of ground motions, and valley effects are responsible for the concentration of damage. The damage in buildings was classified in terms of the intensity and structural typologies for the 30 districts of Izmir metropolitan area. In-depth analysis of the distribution of damages revealed that the earthquake caused damage all over the boundaries of Izmir province, and the concentration of damage in Bornova and Karşıyaka districts has a clear correlation with double resonance effects.

Evaluation of the Local Site Effects and Their Implication to the Seismic Risk of the UNESCO World Heritage Site Old City of Dubrovnik (Croatia)

ABSTRACT UNESCO World Heritage Site Old City of Dubrovnik (Croatia) has been devastated by historical earthquake in 1667 of M ~ 7 as well 1979 M6.8 Montenegro earthquakes causing strong damage effects. The aim of this research is the evaluation of the geological structure and local site effects under the Old City, using non-invasive seismic survey methods to define the areas that can be subjected to stronger damage in the case of an earthquake. The results present basis for seismic hazard and seismic risk assessment of cultural heritage – which are necessary activities for their preservation as an important witness of the past.

Study on multidimensional and multipoint seafloor spatial ground motion simulation under an ice-water layer and the dynamic response of a cross-sea bridge

A combined ice-water reflection coefficient is developed to measure the effect of saltwater and floating ice on seabed ground motions to explore the effect of floating ice on seafloor ground motions. The combined reflection coefficient is made up of three parts: the reflection coefficient of the seawater layer, the transmission coefficient of the water layer, and the reflection coefficient of the ice layer. The combined reflection coefficient analyzes the reflection, transmission, and refraction pathways of seismic waves from the soil layer to the water layer, then to the ice layer and air, and then reflects the attenuation principle of seismic waves in an ice-water layer based on the mechanical parameters of the soil layer, water layer, ice layer, and air. The spatial coherence effect, traveling wave effect, and local site effect are then addressed using the trigonometric series method, and the ice-water layer combined reflection coefficient is utilized to analyze the influence of offshore complex sites on seafloor spatial ground motion. Under the floating ice and seawater layer, an artificial synthetic approach of multidimensional multipoint spatial ground motion is described. The random ice load time history was artificially produced using the ice load spectrum theory and the recorded sea ice parameters. The random response of a long-span cable-stayed bridge is investigated under the combined action of seafloor spatial ground motion and random ice load, and the structural responses of the bridge and tower vary substantially under the combined action of ground motion and ice load.

Stochastic Ground-Motion Simulation of the 2021 Mw 5.9 Woods Point Earthquake: Facilitating Local Probabilistic Seismic Hazard Analysis in Australia

ABSTRACT The 2021 Mw 5.9 Woods Point event is the largest onshore earthquake that has occurred in the recorded history of southeastern Australia since European settlement. To study its source and ground-motion characteristics and to extract information for local seismic hazard analysis, we employ a stochastic finite-fault simulation approach to simulate ground motions for this event based on the observations collected from 36 onshore stations. We determine the regional distance-dependent attenuation parameters using the horizontal Fourier acceleration amplitude spectrum in the frequency range of 0.1–20 Hz. We parameterize path parameters using different models to consider uncertainties and sensitivities. To investigate local site effects, we construct a VS30-based site amplification model. Source parameters are then determined by fitting the theoretical Brune’s ω2 model with a reference Fourier source spectrum at 1.0 km. The κ0 value for the reference rock site is estimated as κ0=0.01 s, and dynamic stress drop is found to be 41.0 MPa by minimizing the overall absolute residual of 5% damped pseudospectral acceleration. We validate the simulations by comparing simulated and observed ground motions in terms of various intensity measurements; analyses of residuals show that the simulations are in good agreement with observations (average residual is close to 0). To facilitate future probabilistic seismic hazard analysis, six selected ground-motion models are ranked using the deviance information criteria based on an independent data set consisting of field observations and simulated ground motions.

Shallow seismic structure around the Vrancea Seismic Zone from joint inversion of ambient noise H/V ratios and surface wave dispersion

The South-East Carpathians in Romania host the Vrancea Seismic Zone (VSZ), an intraplate seismic nest of frequent high magnitude earthquakes that cause unusually high damage over extended areas. Estimating the amplification of seismic signals strongly depends on the depth to the geophysical bedrock and the complex near surface geological structures. To improve seismic hazard estimates and ultimately mitigate seismic risk, we determine the frequency of resonance at 81 broadband seismic stations in southeast Romania using the horizontal-to-vertical spectral ratio (HVSR) of ambient seismic noise. We jointly invert HVSR functions with Rayleigh wave phase velocity dispersion curves using a Markov-chain Monte Carlo algorithm under the assumption of a diffuse field. Our results show excellent correlation between the fundamental frequency of resonance and local surface geology: stations located on foreland sedimentary basins show HVSR peaks in the range 0–1 Hz, while those located in orogenic nappes reach values as high as 10 Hz. Seismic bedrock topography mirrors the crystalline bedrock depth inferred from gravity modeling, but is consistently ∼2000 m shallower, corresponding to the interface between Paleozoic-Triassic to Jurassic-Cretaceous sedimentary levels on the Moesian Platform constrained by borehole data. The near surface seismic model shows anomalously low velocities in the Focsani forearc basin, adjacent to the VSZ, consistent with extremely low resonance frequencies and deep sedimentary layer successions (>5 km). In the back-arc of the Carpathian Orogen, at the contact with the Transylvanian Basin, average resonance frequencies are <2 Hz and inferred seismic bedrock depths vary laterally, likely affected by buried volcanic structures and diapiric salt domes. Our study presents the largest endeavor to study anomalous intensity patterns and local site effects in and around the VSZ and has fundamental implications for seismic risk evaluation and future civil engineering design planning.

Seismic bearing capacity of shallow foundations subjected to inclined and eccentric loading using modified pseudo-dynamic method

In this paper, a comprehensive parametric study is conducted to evaluate the seismic bearing capacity of shallow strip foundations overlying dry and cohesionless granular soil under the actions of inclined and eccentric loadings. For this purpose, a systematic combination of the lower-bound theorems of the limit analysis, the finite element method, and the nonlinear programming is implemented. The seismic loading condition is simulated by the well-established modified pseudo-dynamic (MPD) approach by accounting for the significant influence of phase difference as well as the primary and shear waves propagation through applying non-uniform inertia forces along the vertical and horizontal directions, respectively. Unlike the pseudo-static approach, in the so-called modified pseudo-dynamic loading, a more realistic dynamic nature of the seismic excitation is taken into account using various soil dynamic properties and spectral parameters of the problem under study. In this way, a more viable response of the shallow foundation subjected to seismic loading could be captured by simulating a visco-elastic material through the well-defined Kelvin-Voigt model. The modified pseudo-dynamic (MPD) method considers the influences of shear and primary wave velocities, amplification of earthquake loading, and period of lateral shaking. The current MPD approach is capable of considering the local site effects and damping characteristics of the soil deposit underneath shallow footing by adopting a non-uniform and nonlinear acceleration field inside the soil medium. The results of the seismic bearing capacity are presented in the forms of spectral responses and failure envelopes for the eccentric and inclined loadings. The results show that the failure envelopes of the shallow foundation subjected to either inclined or eccentric loading significantly shrink with the increase in the earthquake accelerations and decrease in the material damping of the underlying soil mass. The amount of changes in the size of failure envelopes in the normalized V-H and V-M spaces due to the variation of seismic intensity and material damping depends directly on the non-dimensional frequency of the earthquake excitation.

Calibration of the local magnitude scale (Ml) for Central Southern Africa

The Central Southern Africa is an area of moderate seismic activity generally caused by the presence of the East Africa Rift System. To improve the quantification of seismicity in this region, we propose a local magnitude scale (Ml), based on the original Richter definition to be used by the Zimbabwe and Botswana national networks. The magnitude scale is developed using 621 seismic events that occurred between 1997–2000 and 2013–2018. These events were recorded by 61 broadband seismic stations located in Botswana, Zimbabwe, and South Africa. We evaluated 5306 traces of zero to peak maximum amplitude, recorded on the horizontal channel from simulated Wood-Anderson seismograms. All 5306 maximum amplitude measurements were inverted simultaneously to determine the attenuation constants. The resultant Ml is defined as $ M_{l} = \log _{10}A+0.80 \times \log _{10}(R)+0.00086 \times R-1.37 \pm S $ in which A is ground displacement amplitude determined from instrument corrected synthetic Wood-Anderson seismograms in nanometres, R is the epicentral distance (km), and S is the station correction. Station corrections were determined for all stations during analysis resulting in values ranging between − 0.44 and + 0.31. The range in station corrections suggests a strong influence of local site effects on the amplitude of the seismic signal. Without station correction, the overall standard deviation of the magnitude residuals using our magnitude relation is 0.32, while using station corrections, the standard deviation 0.17. A comparison of our local magnitude relation with published mblg relationships for Southern Africa shows that these two magnitude scales are almost equivalent. The Ml equation for Central Southern Africa derived in this study has good correlation with mb values reported by the ISC and NEIC and attains systematically lower magnitude values than the South African relationship.

The Importance of Assessing the Geological Site Effects of Ancient Earthquakes from the Archaeoseismological Point of View

Earthquakes have and continue to, occur worldwide, though some places are affected more than others by earthquake-induced ground shaking and the same earthquake can cause more damage in one area than in nearby locations due to site-specific geological site conditions, also known as local site effects. Depending on the chronology of the earthquakes, various disciplines of seismology include instrumental and historical seismology, archaeoseismology, palaeoseismology and neotectonics, each focusing on using specific sources of information to evaluate recent or ancient earthquakes. Past earthquakes are investigated to expand the pre-instrumental and instrumental earthquake catalog and better evaluate a region’s seismic hazard. Archaeoseismology offers a way to achieve these goals because it links how ancient civilizations and their environment might have interacted and responded to past earthquake-induced ground motion and soil amplification. Hence, archaeoseismology explores pre-instrumental (past) earthquakes that might have affected sites of human occupation and their nearby settings, which have left their co-seismic marks in ancient manufactured constructions exhumed by archaeological excavations. However, archaeoseismological observations are often made on a limited epicentral area, poorly constrained dated earthquakes and occasionally on unclear evidence of earthquake damage. Archaeological excavations or field investigations often underestimate the critical role that an archaeological site’s ancient geological site conditions might have played in causing co-seismic structural damage to ancient anthropogenic structures. Nevertheless, the archaeological community might document and inaccurately diagnose structural damage by ancient earthquake shaking to structures and even estimate the size of past earthquakes giving little or no consideration to the role of geological site effects in addressing the causative earthquake. This mixture of factors frequently leads to imprecise estimates of the size of ancient earthquakes and unlikely earthquake environmental impacts, leaving unexplained the location and the moment magnitude of the causative earthquake. Hence, it is essential not to rely solely on earthquake intensities based on archaeologically documented co-seismic damage without assessing the nature of the observed structural damage and the contribution of the geological site effects. This paper explains the geological site effects concept to archaeologists unfamiliar with the notion. It clarifies its role in assessing ground shaking, soil amplification and earthquake intensity by past earthquakes and how and why the geological site effects can be estimated when a site is thought to have been struck by an earthquake. Hence, the geological site effects must be considered when archaeological excavations describe and interpret destruction layers. Conversely, engineers and seismologists dealing with seismic hazard risk assessment must pay close attention to archaeological investigations assessing earthquake intensities and locations based on field evidence of damage to structures attributed to past earthquakes, because the geological site effects might have been factored in inaccurately or not at all.

Fragility curves of masonry buildings in aggregate accounting for local mechanisms and site effects

The seismic evaluation of masonry buildings in aggregate, largely diffused within the existing Italian and European building stock, represents a difficult and open task that has not been exhaustively investigated so far. The study proposes a procedure aimed at evaluating the potential impact of the combination of local mechanisms and site-amplification in terms of fragility curves on an existing unreinforced masonry (URM) aggregate which is made of five adjacent structural units mutually interacting with each other during seismic sequences. The case study is inspired by built heritage of the historic centre of Visso struck by the Central Italy 2016/2017 earthquakes. The in-plane (IP) response of URM buildings was simulated through nonlinear dynamic analyses performed on a 3D equivalent frame model of the structure, whereas out-of-plane (OOP) mechanisms were analysed by adopting the rigid-block assumption but assuming, as seismic input, the floor accelerations derived from the post-processing of data derived from the global 3D model. An innovative procedure considering the pounding effect to the global response of the building is also presented. Two soil conditions were assumed with (freefield) and without (bedrock) site amplification. The results showed that site effects strongly affected the seismic vulnerability of the aggregate, also altering the combination between IP and OOP mechanisms. In fact, for bedrock condition, especially for medium–high damage levels, local mechanisms were prevailing with respect to the IP response. Conversely, for freefield condition, IP mainly governed the overall behaviour for all the damage levels, consistently with the field evidence.