Modification of the liquefaction potential index and liquefaction susceptibility mapping for a liquefaction-prone area (Inegol,Turkey)

Although some liquefaction assessment methods were proposed to evaluate the liquefaction potential of sandy soils, the conventional method based on the standard penetration test (SPT) has been commonly used in most countries and in Turkey. However, it alone is not a sufficient tool for the evaluation of liquefaction potential. The liquefaction potential index was proposed to quantify the severity of liquefaction. Nevertheless, the liquefaction potential index and the severity categories do not answer the question: "Which areas will not liquify?" Besides, the categories do not include a "moderate" category; on the other hand, the "high" and "low" categories are included. This situation is also contrary to the nature of classification schemes. In this study, the liquefaction potential index and the liquefaction potential categories were modified by considering the existing form of the categories based on the liquefaction potential index. While the category of low was omitted, the categories of moderate and "non-liquefied" were adopted. A factor of safety of 1.2 was assumed as the lowest value for the liquefaction potential category of non-liquefied. In addition, the town of Inegol in the Marmara region became the case study for checking the performance of the liquefaction potential categories suggested in this study.     

A numerical model for the transient analysis of offshore foundations under cyclic loading

A comprehensive numerical model for the analysis of offshore foundations under a general transient loading is presented here. The theoretical basis of the model lies on the Swansea formulation of Biot’s equations of dynamic poroelasticity combined with a constitutive model that reproduces key aspects of cyclic soil behaviour in the frame of the theory of generalised plasticity. On the practical side, the adoption of appropriate finite element formulations may prevent the appearance of spurious numerical instabilities of the pore pressure field. In this respect, the use of a coupled enhanced-strain element is here proposed. On the other hand, the practicality of the presented model depends ultimately on its computational efficiency. Some practical recommendations concerning the solution strategies, the matrix storage/handling procedures and the parallel multi-processor computation are here provided. Finally, the performance of the model with a benchmark study case and its practical application to analyse the soil–structure interaction of an offshore monopile under a realistic transient storm loading are discussed.     

Engineering geologic and geotechnical analysis of paleoseismic shaking using liquefaction effects: field examples

The greatest impediments to the widespread acceptance of back-calculated ground motion characteristics from paleoliquefaction studies typically stem from three uncertainties: (1) the significance of changes in the geotechnical properties of post-liquefied sediments (e.g., “aging” and density changes), (2) the selection of appropriate geotechnical soil indices from individual paleoliquefaction sites, and (3) the methodology for integration of back-calculated results of strength of shaking from individual paleoliquefaction sites into a regional assessment of paleoseismic strength of shaking. Presented herein are two case studies that illustrate the methods outlined by Olson et al. [Engineering Geology, this issue] for addressing these uncertainties.

The first case study is for a site near Memphis, Tennessee, wherein cone penetration test data from side-by-side locations, one of liquefaction and the other of no liquefaction, are used to readily discern that the influence of post-liquefaction “aging” and density changes on the measured in situ soil indices is minimal. In the second case study, 12 sites that are at scattered locations in the Wabash Valley and that exhibit paleoliquefaction features are analyzed. The features are first provisionally attributed to the Vincennes Earthquake, which occurred around 6100 years BP, and are used to illustrate our proposed approach for selecting representative soil indices of the liquefied sediments. These indices are used in back-calculating the strength of shaking at the individual sites, the results from which are then incorporated into a regional assessment of the moment magnitude, M, of the Vincennes Earthquake. The regional assessment validated the provisional assumption that the paleoliquefaction features at the scattered sites were induced by the Vincennes Earthquake, in the main, which was determined to have M7.5. The uncertainties and assumptions used in the assessment are discussed in detail.     

A semi‐empirical method for the liquefaction analysis of offshore foundations

This paper presents a relatively simple method for three‐dimensional liquefaction analysis of granular soil under offshore foundations. In this method, the Mohr–Coulomb model, which defines the elasto–plastic stress–strain relationship under monotonic loading, is modified to accommodate the plastic strains generated by cyclic loading. The effects of cyclic loading, evaluated from the results of laboratory tests on saturated samples of soil, are incorporated into the model. The method is implemented in an efficient finite element program for analyses of three‐dimensional consolidating soil. The practicability of the model is demonstrated by analysis of a typical offshore foundation, and the predictions of the numerical analysis are compared with the observed behaviour of the foundation. Copyright © 2000 John Wiley & Sons, Ltd.     

Geotechnical analysis of paleoseismic shaking using liquefaction features: a major updating

A new methodology is proposed for the geotechnical analysis of strength of paleoseismic shaking using liquefaction effects. The proposed method provides recommendations for selection of both individual and regionally located test sites, provides techniques for validation of field data for use in back-analysis, and presents a recently developed energy-based solution to back-calculate paleoearthquake magnitude and strength of shaking. The proposed method allows investigators to qualitatively assess the influence of post-earthquake density change and aging. The proposed method also describes how the back-calculations from individual sites should be integrated into a regional assessment of paleoseismic parameters.     

The 2003 Lefkada earthquake: Field observations and preliminary microzonation map based on liquefaction potential index for the town of Lefkada

A strong earthquake (Mw = 6.2) occurred offshore the island of Lefkada (or Lefkas) on August 14, 2003. The maximum intensity has been evaluated Io = VIII (EMS) at Lefkada municipality, while VI to VII+ intensities were evaluated at many other villages of the island. The offshore NNE-SSW oriented strike-slip right-lateral fault was activated by the main shock. This fault is the northern termination of the Cephalonia transform fault. The most characteristic macroseismic effects were extensive typical ground failures like rock falls, soil liquefactions, subsidence, densification, ground cracks and landslides. These macroseismic effects are remarkably similar to those reported from some historical Lefkada shocks, e.g. 1704, 1914 (Ms = 6.3) and 1948 (Ms = 6.5). Sand boils and ground fissures with ejection of mud were observed at the seaside of the town of Lefkada, and in the villages of Nydri and Vassiliki.In situ soil profiles are obtained based mainly on borings after the earthquake. Boreholes records with SPT values (standard penetration test) are obtained and the “simplified procedure” originally developed by Seed and Idriss was employed to evaluate the liquefaction resistance of soils. The results indicated that the silty sandy layer, which lies beneath the artificial fill in the coastal zone in Lefkada town, had liquefied during the 14 August Earthquake. An attempt was also made to establish a preliminary microzonation map for Lefkada town using the data from Liquefaction Potential Index analyses. Our map was validated by the occurrence of liquefaction phenomena inside the town.     

Artificial neural network and liquefaction susceptibility assessment: a case study using the 2001 Bhuj earthquake data,Gujarat, India

This study pertains to prediction of liquefaction susceptibility of unconsolidated sediments using artificial neural network (ANN) as a prediction model. The backpropagation neural network was trained, tested, and validated with 23 datasets comprising parameters such as cyclic resistance ratio (CRR), cyclic stress ratio (CSR), liquefaction severity index (LSI), and liquefaction sensitivity index (LSeI). The network was also trained to predict the CRR values from LSI, LSeI, and CSR values. The predicted results were comparable with the field data on CRR and liquefaction severity. Thus, this study indicates the potentiality of the ANN technique in mapping the liquefaction susceptibility of the area.     

Comparison of SPT-N-based analysis methods in evaluation of liquefaction potential during the 1999 Chi-chi earthquake in Taiwan

SPT-N-based methods have been adopted for liquefaction assessment of soils during earthquakes for decades. However, there has not been a consistent way of assessing the accuracy and applicability of these methods. The Chi-chi earthquake of 1999, which has been the most serious ground shaking in Taiwan within the century, caused extensive liquefactions in mid-west alluvial deposits of the island. This paper assesses the prediction accuracy of several SPT-N-based methods using liquefaction and non-liquefaction incidents observed during the earthquake. A sensitivity study on commonly adopted parameters shows that the SPT blow count and peak ground acceleration are most sensitive in computing liquefaction potential. By comparing the error in predicting liquefaction and non-liquefaction incidents, this study concludes that Tokimatsu and Yoshimi’s method is more accurate than the other methods. However, the differences between prediction errors of various methods are minimal, indicating all of the methods examined are applicable for the 1999 earthquake in Taiwan.     

A geohazard reconnaissance study based on geoscientific information for development needs of the western region of Istanbul (Turkey)

The role of geoenvironmental information is becoming increasingly important as legislative changes have forced developers and planning authorities to consider more implications and impact on the environment of large-scale development initiatives. Therefore, integration of surface and subsurface geoscientific information for development needs has prime importance and provides a means of identifying potential problems and opportunities at an early stage in any planned development. However, from the experience of recent natural disasters, it is evident that this was not case the taken into consideration in many countries. In addition to thousands of casualties, many urbanized areas, industrial districts and large-scale engineering structures suffered severe damages from the natural hazards due to many reasons including the lack of preliminary engineering geological maps and zoning maps of the settlement areas. Turkey is one of the countries which is exposed to natural hazards such as earthquakes, landslides and floods. In particular, the devastating 1999 Kocaeli earthquake, which affected the Marmara Region of Turkey, focused the attention on densely urbanized and industrialized metropolitan areas such as Istanbul. The rapid growth of Istanbul, particularly towards west with minimal geoscientific information resulted in an overwhelming pressure on the natural environment. In addition, a large earthquake, which is expected to occur in the Marmara Sea within the next 30 years, also pose a threat to the city and its surroundings. In this study, on the basis of the geological, geomorphological and geophysical reconnaissance study, an integrated geoscientific data were collected from the western region of Istanbul and evaluated for geohazards. The paper focuses on the geological and geomorphological aspects that control the occurrence of some geohazards such as earthquake-induced liquefaction, landslides and flooding. In this context, the geological map of the region was revised and Quaternary deposits were classified into 11 units, in detail. Liquefaction-prone areas were evaluated by using geomorphological criteria based on field investigation, by the examination of the available records from 88 boreholes drilled on recent deposits and by the data from resistivity profiles. The landslides within the region were classified according to their type, relative depth and activity. In addition, fluvial and marine flood-prone areas were also delimited within the region. Finally, a series of maps such as landslide inventory maps, and maps showing liquefaction- and flood-prone areas were produced with the aid of Geographic Information Systems (GIS) to assist in designing further detailed site investigations and to reduce costs by ensuring a more focused approach to strategic planning and site selection.     

The great 1950 Assam Earthquake revisited: Field evidences of liquefaction and search for paleoseismic events

Extensive field investigations were carried out for the first time in the meizoseismal area of the great 1950 Assam Earthquake aimed at exploring the paleoseismic history of the NE Indian region through documentation of liquefaction features and radiocarbon (¹⁴C) dating. Trenching at more than a dozen locations along the Burhi Dihing River valley and within the alluvial fans adjoining the Brahmaputra and Dibang Rivers resulted in the identification of more than a dozen very prominent liquefaction features (sand dykes, sills, sand blows etc.) as evidences of large to great earthquakes. ¹⁴C dating of the organic material associated with some of the features indicates a paleoseismic record of about 500 yrs archived by the sediments in this region. Compelling geological evidence(s) of the great 1950 earthquake are well constrained by ¹⁴C dating. Out of the two historically reported seismic events (1548 AD and 1697 AD) from this region, ¹⁴C dating could constrain the 1548 AD event though not distinctly. Further studies using combined ¹⁴C and OSL dating may better constrain the seismo-chronology of the study region.     

An experimental,analytical and numerical study on the thermal efficiency of energy piles in unsaturated soils

Energy piles are bi-functional foundation elements used as structural support as well as ground heat exchangers for shallow geothermal energy systems. Because they are relatively short, energy piles may be partially embedded in unsaturated soils. Saturation conditions influence the thermal properties of the ground and therefore the heat exchange rate, which in turn affects the efficiency of energy piles. This paper combines analytical, experimental and numerical investigations to evaluate the heat exchange rate of energy piles partially or fully embedded in unsaturated soils. The proposed analytical solution is based on the cylindrical heat source theory that treats the soil as a semi-infinite, homogeneous, and isotropic medium. The solution from this theory is multiplied by a function developed analytically in this paper and the outcome is the heat exchange rate for energy piles in unsaturated soils. The proposed function depends on soil saturation, soil and pile thermal properties, and pile geometry. The analytical solution was compared against a finite element solution; which was in turn validated against results from laboratory experiments. Very satisfactory agreements between the analytical, numerical and experimental outputs were observed. The proposed method can be used for a quick and simple evaluation of the efficiency of energy piles in unsaturated conditions. The proposed analytical solution can also be a useful tool for the verification of numerical codes developed for the design of energy piles in unsaturated soils.     

Use of backpropagation neural network for landslide monitoring: a case study in the higher Himalaya

Static and dynamic rock slope stability analyses were performed using a numerical discontinuum modelling technique for a 700-m high rock slope in western Norway. The rock slope has been investigated by the Geological Survey of Norway (NGU), which has been carrying out rock slide studies for the county Møre and Romsdal in western Norway. The purpose of numerical modelling was to estimate the volume of the rock mass that could potentially slide under static and dynamic forces. This estimation was required to assess the run-up heights (tsunami) in a fjord that could potentially be caused by the rockslide. Three cases have been simulated for predicting the behaviour of the rock slope. First, an initial static loading is applied in the numerical model to simulate the prevailing rock mass conditions at the site. Second, saturated and weathered joint conditions are modelled by reducing the residual friction angle along the discontinuities of the rock mass. In doing so, the model simulates the effect of degradation of discontinuities in the rock slope. Third, a dynamic loading, based on peak ground accelerations expected in the area, is applied to simulate dynamic earthquake conditions.

These numerical studies have provided some useful insights into the deformation mechanisms in the rock slope. Both sliding and rotation of blocks start to occur once the residual friction angle along the discontinuities is reduced and when the region is shaken by a strong earthquake. The results indicate that, due to variations in the inclination of discontinuities, the entire slope does not become unstable and that down-slope sliding and rotation of blocks occur mainly on the top layers of the slope. Within the range of parameter values considered for this study, it is unlikely that the whole rock slope can be destabilised. The study provides an illustration of how the geo-mechanical properties of a rock mass can be integrated in a discontinuum rock slope model, which is used for predicting the behaviour of the slope under existing environmental and earthquake conditions. This model has helped not only to better understand the dynamics of the rockslide but also to estimate the potential rock volume that can become unstable when subjected to static and dynamic loads.     

关于喷反装置效能评判的分析与讨论

在开展水平孔某型喷反装置实验研究及测试中发现,吸水效率η这单一因素不能全面、正确地评判或比较其喷反装置的性能。针对喷反效能的评判问题开展了分析与讨论,并应用功率转换因素k值进行喷反效能评判实例分析,提出了用系统压力、吸水效率及功率转换系数综合衡量喷反装置效能的评估观点。     

A statistical procedure for the analysis of seismotectonically induced hydrochemical signals: A case study from the Eastern Carpathians, Romania

Changes in the stress field of an aquifer system induced by seismotectonic activity may change the mixing ratio of groundwaters with different compositions in a well, leading to hydrochemical signals which in principle could be related to discrete earthquake events. Due to the complexity of the interactions and the multitude of involved factors the identification of such relationships is a difficult task. In this study we present an empiric statistical approach suitable to analyse if there is an interdependency between changes in the chemical composition of monitoring wells and the regional seismotectonic activity of a considered area. To allow a rigorous comparison with hydrochemistry the regional earthquake time series was aggregated into an univariate time series. This was realized by expressing each earthquake in form of a parameter “e”, taking into consideration both energetic (magnitude of a seismic event) and spatial parameters (position of epi/hypocentrum relative to the monitoring site). The earthquake and the hydrochemical time-series were synchronised aggregating the e-parameters into “earthquake activity” functions E, which takes into account the time of sampling relative to the earthquakes which occurred in the considered area. For the definition of the aggregation functions a variety of different “e” parameters were considered. The set of earthquake functions E was grouped by means of factor analysis to select a limited number of significant and representative earthquake functions E to be used further on in the relation analysis with the multivariate hydrochemical data set. From the hydrochemical data a restricted number of hydrochemical factors were extracted. Factor scores allow to represent and analyse the variation of the hydrochemical factors as a function of time. Finally, regression analysis was used to detect those hydrochemical factors which significantly correlate with the aggregated earthquake functions.This methodological approach was tested with a hydrochemical data set collected from a deep well monitored for two years in the seismically active Vrancea region, Romania. Three of the hydrochemical factors were found to correlate significantly with the considered earthquake activities. A screening with different time combinations revealed that correlations are strongest when the cumulative seismicity over several weeks was considered. The case study also showed that the character of the interdependency depends sometimes on the geometrical distribution of the earthquake foci. By using aggregated earthquake information it was possible to detect interrelationships which couldn't have been identified by analysing only relations between single geochemical signals and single earthquake events. Further on, the approach allows to determine the influence of different seismotectonic patterns on the hydrochemical composition of the sampled well. The method is suitable to be used as a decision instrument in assessing if a monitoring site is suitable or not to be included in a monitoring net within a complex earthquake prediction strategy.     

A comparative error analysis of manual versus automated methods of data acquisition for algebraic strain estimation

This paper presents a statistical analysis of the algebraic strain estimation algorithm of Shimamoto and Ikeda [Shimamoto, T., Ikeda, Y., 1976. A simple algebraic method for strain estimation from deformed eillipsoidal objects: 1. Basic theory. Tectonophysics 36, 315–337]. It is argued that the error in their strain estimation procedure can be quantified using an expected discrepancy measure. An analysis of this measure demonstrates that the error is inversely proportional to the number of clasts used. The paper also examines the role of measurement error, in particular that incurred under (i) a moment based and (ii) manual data acquisition methods. Detailed analysis of these two acquisition methods shows that in both cases, the effect of measurement error on the expected discrepancy is small relative to the effect of the sample size (number of objects). Given their relative speed advantage, this result favours the use of automated measurement methods even if they incur more measurement error on individual objects. Validation of these results is carried out by means of a simulation study, as well as by reference to studies appearing in previous literature. The results are also applied to obtain an upper bound on the error of strain estimation for various studies published in the literature on strain analysis.     

Dynamic causes for the opening of the Baikal Rift Zone: a numerical modelling approach

Previous dynamic models of the Baikal Rift Zone (BRZ) are mostly two-dimensional on vertical plane. In this study, a numerical model of neotectonics in the region on map view was constructed using the adapted PLATES program. The present work is an attempt to test different mechanisms for opening Baikal Rift by comparing the modelled and observed stress and strain rate fields. The following rifting scenarios were tested: (1) pure northwest–southeast extension, (2) pure northeast–southwest compression, (3) oblique rift opening and (4) combined northwest–southeast extension and northeast–southwest compression. The models are calibrated using geologically and GPS-derived strain rates and stress-tensor determinations from fault-slip data and earthquake focal mechanisms. The most successful model requires a combination of NE–SW compression and orthogonal extension. The model results indicate that the present extensional regime in BRZ can be explained by combining the India plate indentation northward into Eurasia, east–west convergence between the North America and Eurasia plates and southeastward extrusion of the Amur plate in northeastern Asia. Predicted fault-slip rates for the best-fit model are consistent with the observed Holocene fault-slip rates in the Lake Baikal region. The generally accepted rotation of the Amur and Mongolia microplates are used as independent constraints for the choice of the best-fit model. These data correlate well with the predicted direction of rotation in our best model.     

Evaluation of earthquake potential along the Northern Anatolian Fault Zone in the Marmara Sea using comparisons of GPS strain and seismotectonic parameters

Seismotectonic parameters including the Gutenberg-Richter b-value and multifractal dimensions D₂ and D₁₅ of seismicity patterns (both spatial and temporal) were compared to GPS-derived maximum shear and dilatation strains measured in the Marmara Sea region of western Turkey along the Northern Anatolian Fault Zone (NAFZ). Comparisons of seismotectonic parameters and GPS-derived maximum shear and dilatation strain along the NAFZ in the vicinity of the 1999 M7.4 Izmit earthquake reveal a positive correlation (r = 0.5, p = 0.05) between average dilatation and the Gutenberg-Richter b-value. Significant negative correlation (r = − 0.56, p = 0.03 and r = − 0.56, p = 0.02) was also observed between the spatial fractal dimension D₂ and GPS-derived maximum geodetic and shear strain. This relationship suggests that, as maximum geodetic and shear strains increase, seismicity becomes increasingly clustered.Anomalous interrelationships are observed in the Marmara Sea region prior to the Izmit event along a bend in the NAFZ near the eastern end of the Marmara Sea known as the Northern Boundary Fault (NBF). An asperity is located near the northwest end of the NBF. Along the 50-km length of the NBF, GPS strains become slightly compressive. The correlation between b-value and GPS-derived dilatation suggests that regions in compression have increased probability of larger magnitude rupture. The NBF appears to serve as an impediment to the transfer of strain from east to west along the NAFZ. Recurrence times for large earthquakes along the NBF are larger than in surrounding areas. Temporal clustering of seismicity in the vicinity of the NBF may represent foreshocks of an impending rupture.     

On the lag time between internal strain and basement involved thrust induced exhumation: The case of the Colombian Eastern Cordillera

Thrust sheets accumulate internal strain before they start moving along discrete fault planes. However, there are no previous studies evaluating the time difference between initiation of strain and fault displacement. In this paper we use observations from the Eastern Cordillera of Colombia to evaluate this interval. We utilize multiple thermochronometers and paleothermometers to refine the timing of deformation. Based on these new data we build time-temperature path estimates that together with geometric outcrop-based structural analysis and fluid inclusions allow us to assign relative timing to features associated with strain, such as cleavage, veins and certain types of fractures, and compare that with the timing of thrusting. We find that cleavage was only formed close to maximum paleotemperatures, almost coeval with the onset of thrust-induced denudation by the Late Oligocene. The corresponding structural level of fold-related veins suggest that they were formed later but still when the country rocks were at temperatures higher than 160 °C, mostly during the Early Miocene and still coexisted with the latest stages of cleavage formation. Our data show that the main period of strain hardening was short (probably a few million years) and occurred before first-order basement thrusting was dominant, but was associated with second-order folding.     

High-temperature emplacement and liquefaction of shallow-water caldera-forming eruption products: Early Miocene Tateyamazaki Dacite in the Oga Peninsula, NE Japan

The Early Miocene Tateyamazaki Dacite infills a 3.2 km diameter caldera. It comprises poorly sorted, massive, biotite-bearing dacite pumice lapilli tuff, in which huge blocks of densely welded dacite lapilli tuff, basaltic andesite lava, and other lithologies are commonly set. Dense blocks are variably cracked and intruded by the host lapilli tuff. Sparse blocks of bedded lapilli tuff and tuff are variably disaggregated to intermingle with the host rocks or are plastically deformed into irregular shapes. Rootless tuff veins millimeters to 30 cm thick are developed within the host rocks, mainly dipping at 10–30°, and are locally branched and mutually cut to form a network. Where thicker, they are stratified and locally carry accidental fragments. Accidental lapilli up to 2 or 3 cm wide and 30 cm long are locally set in near-vertical and variably sinuous arrays. Although poorly defined they are reminiscent of fluid escape structures. The host pumice lapilli tuff, however, retains in part a thermal remnant magnetization (TRM) vector stable at temperatures above 280 °C. Blocks in the caldera fill also retain TRM but the vectors are rotated significantly from those of the host pumice lapilli tuff and the adjacent volcanic rocks. Tateyamazaki Dacite is thus likely to have been emplaced at high temperatures, and intermingled with shattered basement rocks and ambient water to be partly liquefied within the caldera immediately after or during the caldera-forming eruption.