Site effect assessment using KiK-net data: Part 1. A simple correction procedure for surface/downhole spectral ratios

When data is available, the estimation of site effects is usually performed using the “standard spectral ratio” (SSR) technique with respect to an outcropping, reference rock site. This study uses the Japanese KiK-net network, which has more than 600 pairs of surface-downhole stations allowing the computation of empirical borehole transfer functions, consisting of mean spectral ratios of surface over downhole recordings. The borehole transfer function deviates from the SSR in two respects: the reference is located at depth, and the downhole velocity varies from one site to another. These differences bias the estimation of the transfer function with reference to a standard outcrop rock site. The goal of this paper is to develop a simple and robust methodology to correct for such bias. The proposed correction procedure consists of two steps: a depth correction designed to account, in a simplified and physically acceptable way, for the existence at depth of destructive interferences and the absence of free-surface effects in the high-frequency range; and an impedance correction designed to normalize the shear wave velocity at depth. The depth correction involves a simple, frequency-dependent curve to be adapted for each site as a function of the first destructive interference frequency at depth. The impedance normalization combines the use of “generic” rock velocity profiles and a quarter-wavelength approach, resulting in a smooth frequency-dependent amplitude correction. The proposed methodology is applied on a large subset of KiK-net data in view of analysing the correlation between site amplification factors and site parameters in a companion paper.     

Ground motion prediction in Beirut: a multi-step procedure coupling empirical Green’s functions,ground motion prediction equations and instrumental transfer functions

The strike slip Yammouneh fault is the longest fault in Lebanon, crossing the territory from South to North. It was responsible for major historical earthquakes like the 1202 A.D. earthquake, estimated to M_s7.6. This paper presents a site-specific estimation of the ground motion caused by a potential M_w7.5 earthquake on the Yammouneh fault, similar to the 1202 event, for various sites within the Beirut area. The empirical Green’s function technique EGF is used to estimate the median and the standard deviations of the seismic ground motion at the reference station BHL, taking into account epistemic and aleatory uncertainties related to source parameters. These uncertainties were quantified through a sensitivity analysis of the position of the rupture nucleation X_nuc, the slip roughness parameter K, the corner frequency f_c and the magnitude M_c of the EGF. The rock ground motion is then transferred to various other sites within the Beirut area, using instrumental Fourier transfer functions. Site amplification factors are next deduced by computing the ratio between response spectra at sediment sites and at a reference rock station. Considering the limits of the EGF method in the near field of extended sources, the EGF approach is considered only up to a magnitude M_w of 6.5. Selected Ground Motion Predictive Equations are then used to simulate a M_w7.5 event at a reference station. By applying the amplification factors, the response spectra at the different sites of Beirut are also calculated and compared with the actual response spectra used in the Lebanese regulations.     

Modeling coupled erosion and filtration of fine particles in granular media

One of the major causes of instability in geotechnical structures such as dikes or earth dams is the phenomenon of suffusion including detachment, transport and filtration of fine particles by water flow. Current methods fail to capture all these aspects. This paper suggests a new modeling approach under the framework of the porous continuous medium theory. The detachment and transport of the fine particles are described by a mass exchange model between the solid and the fluid phases. The filtration is incorporated to simulate the filling of the inter-grain voids created by the migration of the fluidized fine particles with the seepage flow, and thus, the self-filtration is coupled with the erosion process. The model is solved numerically using a finite difference method restricted to one-dimensional (1-D) flows normal to the free surface. The applicability of the model to capture the main features of both erosion and filtration during the suffusion process has been validated by simulating 1-D internal erosion tests and by comparing the numerical with the experimental results. Furthermore, the influence of the coupling between erosion and filtration has been highlighted, including the development of material heterogeneity induced by the combination of erosion and filtration.

     

Structural mapping in basin-and-range-like geology by electromagnetic methods: a powerful aid to seismic

A case history is presented where electromagnetic (EM) methods were applied as a complement to seismic, for structural mapping in basin-and-range-like geology: 366 five-component magnetotelluric (MT) soundings were carried out together with 331 transient soundings (TDEM) along seismic lines. Due to high structural complexity, seismic shows a number of limitations. For the same reasons, MT is highly perturbed and three specific interpretation techniques were comprehensively applied: 1. a classical correction of static effect using TDEM sounding, to determine the high-frequency nondistorted apparent resistivities and thus the corrected tensor; 2. a so-called regional correction based upon the same concept as the static effect, to transform distorted resistivity curves due to the horst/graben situation into plausible 1D curves, through the use of nomograms built for 2D H-polarization situations; 3. a stripping technique which made it possible to map areas where a deep conductive Mesozoic shale was present below carbonates, at a depth of 3 km. After the best MT interpretation was obtained along each line, it was integrated with seismic and with the results from two boreholes. A crude empirical law relating resistivity and acoustic velocity was established and the MT horizons were plotted on the two-way traveltime seismic sections. The final integrated cross-sections obtained are undoubtedly of greater use to the explorationist than the initial seismic sections alone and two wells were accurately predicted.     

Temperature distribution in karst systems: the role of air and water fluxes

A better understanding of heat fluxes and temperature distribution in continental rocks is of great importance for many engineering aspects (tunnelling, mining, geothermal research, etc.). This paper aims at providing a conceptual model of temperature distribution in karst environments which display thermal ‘anomalies’ as compared with other rocks. In temperate regions, water circulation is usually high enough to ‘drain‐out’ completely the geothermal heat flux at the bottom of karst systems (phreatic zone). A theoretical approach based on temperature measurements carried out in deep caves and boreholes demonstrates, however, that air circulation can largely dominate water infiltration in the karst vadose zone, which can be as thick as 2000 m. Consequently, temperature gradients within this zone are similar to the lapse rate of humid air (～0.5 °C 100 m^?1). Yet, this value depends on the regional climatic context and might present some significant variations.     

Seismic Microzonation: A Comparison between Geotechnical and Seismological Approaches in Pointe-à-Pitre (French West Indies)

Pointe-à-Pitre, the main city of Guadeloupe in the French West Indies, has on several occasions been partially destroyed by major historical earthquakes. Moreover, a post-seismic assessment of the damage from the 1985 Montserrat earthquake indicates that the town is prone to site effects. Consequently, from 1996 to 1998, BRGM conducted a seismic microzonation study based on geotechnical and geological data. At the same time, three seismological studies were being conducted – two based on earthquake recordings using a time-series analysis and the classical spectral ratio (CSR) method (CETE/LCPC and BRGM), and the third based on noise measurement at 400 points using the horizontal-to-vertical noise ratio (HVNR) method (CETE/LCPC). The objective of this paper is not to carry out a new microzonation study by taking into account all the results, but rather to show in what respects the results of these different methods are in agreement or not. A comparison of the results of the seismological studies with the geotechnical microzonation shows that they are in fairly good agreement, albeit with some discrepancies. The results indicate that the seismological methods and the geotechnical data are highly complementary and should be used together in compiling seismic transfer-function microzonation maps. This revised version was published online in July 2006 with corrections to the Cover Date.     

Internal erosion in dike-on-foundation modeled by a coupled hydromechanical approach

One of the major causes of instability in geotechnical structures such as dikes or earth dams is internal erosion, an insidious process that occurs over a long period of time. Research on this topic is still fairly new and much more needs to be understood in order to solve the problems posed by this phenomenon. This paper proposes a hydromechanical model based on porous continuous medium theory to assess how internal erosion impacts the safety of earthen structures. The saturated soil is considered as a mixture of four interacting constituents: soil skeleton, erodible fines, fluidized fine particles, and fluid. The detachment and transport of the fine particles are described by a mass exchange model between the solid and the fluid phases. An elastoplastic constitutive model for sand-silt mixtures has been developed to monitor the effect of the evolution of both porosity and fines content induced by internal erosion upon the behavior of the soil skeleton. The model has been numerically solved with the finite element method. It has then been applied to the specific case study of a dike foundation subjected to internal erosion induced by the presence of a karstic cavity beneath the alluvium layer. The numerical results show the onset of erosion, the time-space evolution of the eroded zone, and the hydromechanical response of the soil constituting the dike, all of which highlights the effects of the cavity location, the erosion rate, and the fines content.     

Geophysical characterisation of karstic networks – Application to the Ouysse system (Poumeyssen,France)

In the framework of the management of karstic aquifers, geophysical reconnaissance can be used to locate conduits and caves, and to characterise the surrounding limestone matrix. Suitable characterisation of heterogeneities in the karstic environment is, however, challenging for ground-based geophysical methods. The present article describes the results, and evaluates the response and accuracy of combined geophysical measurements carried out at the Poumeyssen test site in France, involving electrical resistivity imaging (ERI), magnetic resonance sounding (MRS), “mise-à-la-masse” electrical mapping, and seismic tomography. This site provides the opportunity to study a relatively wide, shallow, water-filled conduit whose location and shape are known from topographic work carried out by cave divers. Seismic and MRS provided the exact location and width of the conduit, to within a few meters. The seismic and electrical data suggest that the limestone medium surrounding the conduit is not homogeneous.