Integrated geophysical methods for the characterisation of an archaeological site (Massenzio Basilica — Roman forum, Rome, Italy)

A geophysical study that involved different techniques was carried out with the aim to improve the knowledge of the archaeological site where the Basilica of Maxentius was founded and to discern individual covered structures (foundations).Vertical Electrical Sounding (VES), seismic refraction and Electrical Resistivity Tomography (ERT) studies were performed in the archaeological site. VES and seismic refraction allowed to characterise the main geological formations of the hill where the Basilica was built and to distinguish the concrete floor and backfilling. Electrical data were processed using different algorithms; their results were compared to appraise the inverted models' robustness.ERT inversion algorithms were used to delineate shape and size of a much more complex structure, that were originally expected from archaeological excavation plan. The results of the commercial program were used as a posteriori information to include them in the algorithm proposed by the authors; the sequential use of the programs defined a processing procedure.The integrated use of different geophysical techniques reduced a great deal the intrinsic ambiguities of each method. Direct explorations (boreholes and archaeological excavations) confirmed the geophysical results.     

Mapping soil effects in Barcelona, Spain, through an integrated GIS environment

Maps of soil response for the city of Barcelona have been for the first time obtained through a GIS environment that integrates the different procedures for soil response estimation within a single tool. These maps constitute part of the results of the local scale application of the computer prototype for seismic risk assessment that was developed within the European project SERGISAI. The approach involves: collection of available data relevant to local geology, implementation of geotechnical models, estimation of the reference seismic action, generation of synthetic strong-motion time histories, and soil response calculation through 1D analytical method. The resulting predictive hazard maps of predominant period and amplification ratio delineate potential variations on ground shaking and constitute a first approximation towards an integrated approach to Barcelona urban area microzonation. Analysis of the observed differences, when comparing the analytical results in this study with previous empirical studies, provides a useful feedback to establish site dependence suitability and reliability of methods, to extract information on at-present inaccessible parameters needed for the characterisation of physical properties of soil, and also to delimit those areas where further in-depth survey research is needed for a proper seismic hazard assessment.     

王家沟断层组与九家湾断层组转换区的地球物理研究

针对王家沟断层组与九家湾断层组转换区基岩走势及目标地层的变化情况,确定适合转换区物探方法探测的观测系统参数,并辅以现场踏勘,确定合理的探测测线,利用wseis15和RES2DINV软件对在探测测线上采集到的原始数据进行处理,得到时间剖面和地层视电阻率剖面,比对两条断层组已有物探资料中土层变化的特征,总结归纳出在转换区内土层变化情况。     

A stochastic based approach for a new site classification method: application to the Algerian seismic code

Building codes have widely considered the shear wave velocity to make a reliable subsoil seismic classification, based on the knowledge of the mechanical properties of material deposits down to bedrock. This approach has limitations because geophysical data are often very expensive to obtain. Recently, other alternatives have been proposed based on measurements of background noise and estimation of the H/V amplification curve. However, the use of this technique needs a regulatory framework before it can become a realistic site classification procedure. This paper proposes a new formulation for characterizing design sites in accordance with the Algerian seismic building code(RPA99/ver.2003), through transfer functions, by following a stochastic approach combined to a statistical study. For each soil type, the deterministic calculation of the average transfer function is performed over a wide sample of 1-D soil profiles, where the average shear wave(S-W) velocity, V s, in soil layers is simulated using random field theory. Average transfer functions are also used to calculate average site factors and normalized acceleration response spectra to highlight the amplification potential of each site type, since frequency content of the transfer function is significantly similar to that of the H/V amplification curve. Comparison is done with the RPA99/ver.2003 and Eurocode8(EC8) design response spectra, respectively. In the absence of geophysical data, the proposed classification approach together with micro-tremor measures can be used toward a better soil classification.     

Joint Audio-Magnetotelluric and Passive Seismic Imaging of the Cerdanya Basin

The structure of Cerdanya Basin (north-east of Iberian Peninsula) is partly known from geological cross sections, geological maps and vintage geophysical data. However, these data do not have the necessary resolution to characterize some parts of Cerdanya Basin such as the thickness of soft soil, geometry of bedrock or geometry of geological units and associated faults. For all these reasons, the main objective of this work is to improve this deficiency carrying out a detailed study in this Neogene basin applying jointly the combination of passive seismic methods (H/V spectral ratio and seismic array) and electromagnetic methods (audio-magnetotelluric and magnetotelluric method). The passive seismic techniques provide valuable information of geometry of basement along the profile. The maximum depth is located near Alp village with a bedrock depth of 500 m. The bedrock is located in surface at both sites of profile. The Neogene sediments present a shear-wave velocity between 400 and 1000 m/s, and the bedrock basement presents a shear-wave velocity values between 1700 and 2200 m/s. These results are used as a priori information to create a 2D resistivity initial model which constraints the inversion process of electromagnetic data. We have obtained a 2D resistivity model which is characterized by (1) a heterogeneous conductivity zone (<40 Ohm m) that corresponds to shallow part of the model up to 500 m depth in the centre of the profile. These values have been associated with Quaternary and Neogene sediments formed by silts, clays, conglomerates, sandstones and gravels, and (2) a deeper resistive zone (1000–3000 Ohm m) interpreted as Palaeozoic basement (sandstones, limestones and slates at NW and conglomerates and microconglomerates at SE). The resistive zone is truncated by a discontinuity at the south-east of the profile which is interpreted as the Alp-La Tet Fault. This discontinuity is represented by a more conductive zone (600 Ohm m approx.) and is explained as a combination of fractured rock and a fluid network. The result highlights that the support between different geophysical methods is essential in producing geophysical meaningful models.     

Soil characterization for seismic damage scenarios purposes: application to Angra do Heroísmo (Azores)

Angra do Heroismo, the main town of Terceira Island in the Azores Archipelago, was hit in 1980 by a 7.2 magnitude earthquake that caused great destruction in the central part of the town. Taking into consideration the high seismic hazard of the region and the cultural and social importance of Angra do Heroísmo, the elaboration of damage scenarios is of particular importance to implement measures for preserve and protect the town against future earthquakes. The first step is to perform microzonation studies in order to characterize the soil seismic behaviour. Taking into consideration the available geologic, geotechnical and geophysical information, a detailed soil characterization was performed based on the results from numerical modelling and the analysis of microtremor experimental measurements. Nine different soil profiles were identified, characterized and classified. Discussion on the detailed soil classification and the Eurocode 8 soil classification is presented. This study shows that even with an available code, microzonation studies must be developed in order to identify differences on soil behaviour inside the interested area. It shows also that the use of experimental measurements presents a great help on soil characterization. The obtained detailed classification will be used on the estimation of damage scenarios for Angra do Heroísmo.     

Structural study of the Ballandaz landslide (French Alps) using geophysical imagery

This study of the Ballandaz landslide (Savoie, French Alps) was carried out as part of the ANR/PGCU-SIGMA research project. Its aim was to characterize the slope by different methods of geophysical imagery, and then use this characterization to provide a combined interpretation of the morpho-structures in order to simulate the geomechanical behavior of the sliding mass. First, electromagnetic mapping was used to identify the variations of the surface lithology and moisture from which one can more precisely locate the active zone of the slope. Then, electrical and seismic 2D imagery methods were used along several transverse and longitudinal profiles in order to produce electrical resistivity and seismic-velocity depth sections showing up the slope's structures. The H/V method was also tested locally to determine the depth to basement so as to complement the profiles: it reveals potential seismic site effects from the deeper structures. Helped by drilling and inclinometer surveys, the geophysical interpretations have revealed the various units structuring the landslide: (i) the active, very heterogeneous shallow level overlying (ii) a more rigid, less porous, and probably stable shallow bedrock, (iii) boulders, and (iv) the sound geological basement of quartzite and gypsiferous facies. The study has shown the usefulness of combining such different sounding techniques for studying complex environments like landslides; it has also revealed the limitations of each method when used for studying very heterogeneous environments.     

A seismic multi-approach method for characterizing contaminated sites

A geophysical campaign to characterize the subsurface of a contaminated site down to a depth of several tens of meters was carried out under the HYGEIA-CEE project. On this site, seismic techniques were combined to image the geological structures; i.e. seismic reflection, P-wave tomography and spectral analysis of surface waves. Because these techniques consider different wave components in the processing, they can be expected to provide complementary information concerning the site lithology. The special feature of this experiment is the fact that the same seismic acquisition device, consisting of a mobile central unit, a drop-weight seismic source, and a sensor line of gimbal mounted geophones, was used for each of the techniques. Two perpendicular seismic lines were set up in the field for testing two geophone spacings. Three processing procedures, one each for the seismic reflection, P-wave tomography and spectral analysis of surface waves, were developed for producing seismic images from the P-wave reflectivity, the first P-wave arrivals and the dispersion of Rayleigh waves, respectively. The images show good complementarity in terms of investigation depth. The results are also in good agreement with available borehole data: the sandy layers seem to be related to low velocities, since the high velocities are better explained by the presence of clayey and gravelly intervals. The contribution and the limits of this seismic multi-approach method is discussed.     

A Geophysical Study of the Carcavai Fault Zone, Portugal

The Algarve province is located a few hundred kilometres north of the crossing of the E–W Eurasia–Africa plate boundary in an area of diffuse seismicity and broad deformation. It is characterised by a moderate seismicity, with some important historical and instrumental earthquakes causing loss of lives and significant material damages. The area is affected not only by plate boundary earthquakes but also by moderate to large events generated by local sources. The assessment of onshore local sources is, therefore, of vital importance for an evaluation of the regional seismic hazard. This paper discusses the application of geophysical data to the study of the Carcavai fault zone, an outcropping structure more than 20?km long which is seen to deform sediments of the Plio-Quaternary age. The location of some sectors of the fault zone, as well as the vertical offsets of the structure, are still to be confirmed. In order to estimate these and to study the geometry of the fault zone at depth, geophysical data were acquired together with new geological data. Where the location of the fault was less certain, EM and seismic reflection profiles with coarse spatial sampling were carried out. After the detailed location of the fault zone, seismic reflection profiles with a more dense spatial resolution were acquired. The integrated interpretation of the geological and geophysical data confirmed the presence of a large fault zone. The total fault length is still unknown as its extension offshore is still being studied. Together with estimated values of the throw obtained, this data set has improved understanding the seismic hazard in the area by providing more refined estimates of co-seismic rupture, maximum expected earthquake and return periods.     

Characterization of seismic hazard and structural response by energy flux

Seismic safety of structures depends on the structure's ability to absorb the seismic energy that is transmitted from ground to structure. One parameter that can be used to characterize seismic energy is the energy flux. Energy flux is defined as the amount of energy transmitted per unit time through a cross-section of a medium, and is equal to kinetic energy multiplied by the propagation velocity of seismic waves. The peak or the integral of energy flux can be used to characterize ground motions. By definition, energy flux automatically accounts for site amplification. Energy flux in a structure can be studied by formulating the problem as a wave propagation problem. For buildings founded on layered soil media and subjected to vertically incident plane shear waves, energy flux equations are derived by modeling the building as an extension of the layered soil medium, and considering each story as another layer. The propagation of energy flux in the layers is described in terms of the upgoing and downgoing energy flux in each layer, and the energy reflection and transmission coefficients at each interface. The formulation results in a pair of simple finite-difference equations for each layer, which can be solved recursively starting from the bedrock. The upgoing and downgoing energy flux in the layers allows calculation of the energy demand and energy dissipation in each layer. The methodology is applicable to linear, as well as nonlinear structures.     

Seismic technique to determine the allowable bearing pressure for shallow foundations in soils and rocks

Based on a variety of case histories of site investigations, including extensive bore-hole data, laboratory testing and geophysical prospecting, an empirical formulation is proposed for the rapid determination of allowable bearing capacity of shallow foundations. The proposed expression consistently corroborates the results of the classical theory and is proven to be rapid and reliable. It consists of only two soil parameters, namely, the in situ measured shear wave velocity, and the unit weight. The unit weight may be also determined with sufficient accuracy by means of another empirical expression using the P-wave velocity. It is indicated that once the shear and P-wave velocities are measured in situ by an appropriate geophysical survey, the allowable bearing capacity as well as the coefficient of subgrade reaction and many other elasticity parameters may be determined rapidly and reliably through a single step operation, not only for soils, but also for rock formations. Such an innovative approach, using the seismic wave velocities only, is considerably cost- and time-saving in practice.     

Relationships between water flow rate and geophysical measurements in an alluvial aquifer

Our study focuses on the potential usefulness of surface geophysical data to constrain the water content within an alluvial aquifer. On a study area where two wells have been drilled, we have performed several geophysical measurements, including ground penetrating radar, DC resistivity prospecting, seismic refraction survey and magnetic resonance soundings. From these data, we estimated several parameters, namely, the water height in the deposits, the effective porosity, the water content, the permeability, and the transmissivity of alluvial deposits. These physical parameters allow us to characterize the alluvial deposits in order to constrain the estimation of the potential water flow. The lithology and water flow rate known from the wells enabled us to compare geophysical results obtained in a high water flow rate zone to those in a low water flow rate zone. Correlation has been found between the water flow rate observed in both wells and the geophysical data obtained in the vicinity of these wells.     

Characterization of recharge through complex vadose zone of a granitic aquifer by time-lapse electrical resistivity tomography

The vadose zone is the main region controlling water movement from the land surface to the aquifer and has a very complex structure. The use of non-invasive or minimally invasive geophysical methods especially electrical resistivity imaging is a cost-effective approach adapted for long-term monitoring of the vadose zone. The main aim of this work is to know the fractures in the vadose zone, of granitic terrene, through which the recharge or preferred path recharge to the aquifer takes place and thus to relate moisture and electrical resistivity. Time lapse electrical resistivity tomography (TLERT) experiment was carried out in the vadose zone of granitic terrene at the Indian Geophysical Research Institute, Hyderabad along two profiles to a depth of 18 m and 13 m each. The profiles are 300 m apart. Piezometric, rainfall and soil moisture data were recorded to correlate with changes in the rainfall recharge. These TLERT difference images showed that the conductivity distribution was consistent with the recharge occurring along the minor fractures. We mapped the fractures in hard rock or granites to see the effect of the recharge on resistivity variation and estimation of moisture content. These fractures act as the preferred pathways for the recharge to take place. A good correlation between the soil moisture and resistivity is established in the vadose zone of granitic aquifer. Since the vadose zone exhibits extremely high variability, both in space and time, the surface geophysical investigations such as TLERT have been a simple and useful method to characterize the vadose zone, which would not have been possible with the point measurements alone. The analyses of the pseudosection with time indicate clearly that the assumption of the piston flow of the moisture front is not valid in hard rocks. The outcome of this study may provide some indirect parameters to the well known Richard's equation in studying the unsaturated zone.     

Local site effect microzonation of Lorca town (SE Spain)

Local site effect assessment based on subsurface ground conditions is often the key to evaluate urban seismic hazard. The site effect evaluation in Lorca town (south-eastern Spain) started with a classification of urban geology through the geological mapping at scale 1:10,000 and the use of geotechnical data and geophysical surveys. The 17 geological formations identified were classified into 5 geological/seismic formations according to their seismic amplification capacity obtained from ambient vibration measurements as well as from simultaneous strong motion records. The shear-wave velocity structure of each geological/seismic formation was evaluated by means of inversion of Rayleigh wave dispersion data obtained from vertical-component array records of ambient noise. Nakamura’s method was applied to determine a predominant period distribution map. The spectral amplification factors were fourfold the values recorded in a reference hard-rock site. Finally, the capability of this study for explaining the damage distribution caused by the May 11th, 2011 Lorca destructive earthquake (Mw \(=\) 5.2) was examined. The methods used in this work are of assistance to evaluate ground amplification phenomena in urban areas of complex geology as Lorca town due to future earthquakes with applicability on urban seismic risk management.     

2005年11月26日九江—瑞昌地震场地影响的初步分析

分析了2005年11月26日九江-瑞昌地震震害与当地工程场地条件间的关系,结果表明：当地工程地质条件复杂,水平向变化大,软土发育,存在溶洞,是此次地震灾害严重的一个重要原因。灾后重建应开展地震小区划工作,为抗震设防提供科学依据。     

Delineation of potential seismic sources for seismic zoning of Iran

A total of 235 potential seismic sources in Iran and neighboring regions are delineated based on available geological, geophysical, tectonic and earthquake data for seismic hazard assessment of the country. In practice, two key assumptions are considered; first, the assumption of earthquake repeatedness, implying that major earthquakes occur preferentially near the sites of previous earthquakes; second, the assumption of tectonic analogy, which implies that structures of analogous tectonic setting are capable of generating same size earthquakes. A two-step procedure is applied for delineation of seismic sources: first, demarcation of seismotectonic provinces; second, determination of potential seismic sources. Preferentially, potential seismic sources are modeled as area sources, in which the configuration of each source zone is controlled, mainly, by the extent of active faults, the mechanism of earthquake faultings and the seismogenic part of the crust.     

Coseismic Damage at an Archaeological Site in Sicily,Italy: Evidence of Roman Age Earthquake Surface Faulting

Archaeoseismology can provide a useful chronological tool for constraining earthquakes and documenting significant evidence that would otherwise be lost. In this paper, we report a case of surface faulting on ancient man-made structures belonging to the archaeological site of Santa Venera al Pozzo situated along the eastern flank of Mt. Etna volcano in eastern Sicily (southern Italy), which is affected by well-developed tectonic faults. Geological surveys highlight a set of fractures affecting the archaeological ruins, suggesting the occurrence of a capable fault zone across the area. An integrated geophysical survey was carried out in order to identify the main subsurface tectonic discontinuity ascribable to the fault zone. The information derived from different geophysical techniques, such as electrical resistivity tomography, seismic refraction tomography, ground-penetrating radar, and magnetic surveys allowed us to infer that the fractures observed at the surface could have been produced by coseismic rupture. They are conceivably linked to a strong earthquake that probably occurred in the Roman period, around mid-end of the third-century AD; time constraints are inferred through the dating of buildings of the archaeological site.     

Geophysical investigation of the Sandalp rock avalanche deposits

In the study of rock avalanche phenomena, numerical modelling makes use of back analyses of the rock avalanche propagation for calibration of the modelling assumptions and parameters. The back analyses require knowledge of the run-out area boundaries and the thickness distribution of the deposit. Geophysical methods can be applied to retrieve the thickness distribution, but, due to strong heterogeneities and logistic problems they are seldom applied. The aim of this work is to assess the potential of integrated geophysical methods to recognise and characterise a deposit created by two rock avalanches which occurred in the Sandalp valley (Switzerland) in 1996. The topography of the site before and after the rock avalanche is known and can be used as a benchmark. Resistivity tomography, seismic P-wave tomography, and active and passive surface wave analysis have been applied on several profiles deployed both on the rock avalanche deposit and in the surrounding area. Innovative approaches for surface wave analysis based on laterally constrained inversion and multimodal inversion have been applied to the data. A comparison of the results of the geophysical investigations with the topographic benchmark has shown the capability of the geophysical methods to locate the bottom of the deposit in the areas where the contrast with the host sediments properties is significant. In these areas, the deposit has higher resistivities and lower seismic velocities than the underlying materials. In the areas where the deposit is thicker and richer in fine-grained materials the geophysical parameters are not able to discriminate between the rock avalanche deposit and the underlying sediments. As a secondary task, the geophysical methods also allowed the bedrock pattern to be outlined.     

汶川地震序列早期特征及成因探讨

地震序列的特征和震型判定工作有助于抗震救灾工作的开展,对其发生成因的研究是解决地震预报难题必须面对的科学问题。汶川8．0级地震序列的初步研究表明：①余震沿龙门山断裂带分布于宽100km,长约330km的带状区域内,并侧向于主震震中的北侧;②序列发展初期有2个快速衰减过程;③序列类型为主震一余震型,最大强余震6．4级;④序列的空间演化过程,强余震震源机制结果和地震精确定位结果分析表明,序列具有分段特征;⑤8．0级地震的发震构造是龙门山断裂带,发震构造在剖面上呈现出“犁形”或“铲形”。地球物理勘探和壳、幔结构反演结果表明,自青藏高原穿越龙门山到四川盆地存在地幔阶梯,上地幔阶梯的阻挡作用使得物质东移速率减慢,并蕴积了汶川8．0级地震所需能量。地震的产生正是东西向应力平衡被打破的结果,余震沿龙门山断裂的分布是高原地壳在印度板块的推挤作用下向北北东方向的运动得以继续的表现。整个龙门山断裂带都参与了活动,龙门山断裂带北端作为断裂带的止裂端与南段同期活动。     

Identification of accelerometric stations in ITACA with distinctive features in their seismic response

In this work, we analyze the records of the Italian strong motion database (ITACA, ) with the aim of identifying stations affected by site effects that are not captured by standard seismic classification schemes. In particular, we consider four different site classifications, two of them based on geological/geophysical characteristics and two driven by data. For each classification we develop a ground motion prediction equation using a random effect approach to isolate the between-station and within-station distribution of errors. The site coefficients obtained for the different classes confirm that site amplification effects are significant for both the horizontal and vertical components. The between-station error normalized to the standard deviation of the between-station error distribution is then used to identify stations characterized by large errors, attributable to site effects not accounted for by the classification schemes. The results show that large errors can affect the predictions when the site effects are not uniquely related to the reduction of the seismic impedance in the uppermost layers. For example, amplifications of ground motion over the long period range are observed for stations installed within alluvial closed-shape basins, as consequence of locally generated surface waves. For these stations, classifications based on the horizontal to vertical response spectra ratio are not reliable, since amplifications are also affecting the vertical component. Another interesting feature which emerges from the analysis is the significant de-amplification of short period spectral ordinates that seems to be related to stations typically set in at the foundation level of massive structures. To increase the usefulness of the data set, the most important distinctive features of the strong motion stations are documented in the ITACA database reports containing the instrument information and the available geological-geotechnical data.