Influence of vertical ground motion on seismic-induced displacements of gravity structures

A method enabling assessment of seismic-induced movements of gravity block is proposed. This method is based on dynamic analysis of the block subjected to both horizontal and vertical ground accelerations, and resultant forces acting on the block. There are two kinds of the motion. The first one is trivial, as it is a common motion of the block and subsoil. The second is a relative motion of the block with respect to subsoil, depending on combination of ground accelerations and the resultant of horizontal forces. Attention is focused on the influence of vertical ground motion on seismic-induced displacements of the structure. Experimental investigations performed in order to validate a numerical algorithm are described, and then experimental data are compared with numerical predictions. A good agreement is obtained. Then, the dynamics analysis of the block is presented for harmonic and pseudo-stochastic ground motions in order to illustrate the algorithm proposed. The results of computations show that, for some typical data and ground accelerations corresponding to strong earthquakes, large permanent displacements of the block take place, as observed in real conditions.     

A deterministic seismic hazard analysis for shallow earthquakes in Greece

The maximum expected ground motion in Greece is estimated for shallow earthquakes using a deterministic seismic hazard analysis (DSHA). In order to accomplish this analysis the input data include an homogeneous catalogue of earthquakes for the period 426 BC–2003, a seismogenic source model with representative focal mechanisms and a set of velocity models. Because of the discrete character of the earthquake catalogue and of errors in location of single seismic events, a smoothing algorithm is applied to the catalogue of the main shocks to get a spatially smoothed distribution of magnitude. Based on the selected input parameters synthetic seismograms for an upper frequency content of 1 Hz are computed on a grid of 0.2° × 0.2°. The resultant horizontal components for displacement, velocity, acceleration and DGA (Design Ground Acceleration) are mapped. The maps which depict these results cannot be compared with previously published maps based on probabilistic methodologies as the latter were compiled for a mean return period of 476 years. Therefore, in order to validate our deterministic analysis, the final results are compared with PGA estimated from the maximum observed macroseismic intensity in Greece during the period 426 BC–2003.Since the results are obtained for point sources, with the frequency content scaled with moment magnitude, some sensitivity tests are performed to assess the influence of the finite extent of fault related to large events. Sensitivity tests are also performed to investigate the changes in the peak ground motion quantities when varying the crustal velocity models in some seismogenic areas. The ratios and the relative differences between the results obtained using different models are mapped and their mean value computed. The results highlight the importance in the deterministic approach of using good and reliable velocity models.     

An intelligent approach to prediction and control ground vibration in mines

Drilling and Blasting are still considered to be the most economical method for rock excavation either on surface or underground. The explosive energy, which breaks the rockmass, is not fully utilized for this purpose. Only 20–30% of explosive energy is utilized for fragmenting the rockmass and the rest wasted away in the form of ground vibration, air blast, noise, fly rock, back breaks, etc. Among them, ground vibration is considered to have the most damaging effect. A number of predictor equations have been proposed by various researchers to predict ground vibration prior to blasting. Still, it is difficult to recommend any one predictor for a particular ground condition because ground vibration is influenced by a number of parameters. These parameters are either controllable or non-controllable like blast geometry, explosive types, rock strength properties, joints patterns, etc. In the present paper, an attempt has been made to predict the ground vibration using an Artificial Neural Network incorporating large number of parameters, which affect the ground vibration. Results are also compared with the values obtained from regression analysis and observed field data sets. Finally, it is found that the neural network approach is more accurate and able to predict the value of blast vibration without increasing error with increasing number of inputs and non-linearity among these.     

Method of anthropoecological stabilization of landscape and its interpretation

The study deals with the theoretical approach to anthropoecological stabilization of landscape (ASL). It describes a method consisting of four stages including detailed differentiation from the viewpoint of aim, objectives of investigation, and means of evaluation and processing. Application of the method requires to provide: the preparatory, analytical ground — I. stage (analytical part), the structural — functional evaluation of the territory under solution — II. stage (analytical-synthetical part), the solution of ecooptimization relations and bonds by way of material and energy flows — III. stage (regulative part), the monitoring and signalling activities, determination of all-system biogeoindications — IV. stage (ensuring part).The results are exemplified by application to model territories in South Bohemia region.     

Seismic Hazard Mapping and Microzonation in the Sikkim Himalaya through GIS Integration of Site Effects and Strong Ground Motion Attributes

The seismic ground motion hazard is mapped in the Sikkim Himalaya with local and regional site conditions incorporated through geographic information system. A strong motion network in Sikkim comprising of 9 digital accelerographs recorded more than 100 events during 1998–2002, of which 41 events are selected with signal-to-noise ratio 3 for the estimation of site response (SR), peak ground acceleration (PGA) and predominant frequency (PF) at all stations. With these and inputs from IRS-1C LISS III digital data, topo-sheets, geographical boundary of the State of Sikkim, surface geological maps, soil taxonomy map in 1:50,000 scale and seismic refraction profiles, the seismological and geological thematic maps, namely, SR, PGA, PF, lithology, soil class, %slope, drainage, and landslide layers are generated. The geological themes are united to form the basic site condition coverage of the region. The seismological themes are assigned normalized weights and feature ranks following a pair-wise comparison hierarchical approach and later integrated to evolve the seismic hazard map. When geological and seismological layers are integrated together through GIS, microzonation map is prepared. The overall site response, PGA and predominant frequency show an increasing trend in the NW–SE direction peaking at Singtam in the lesser Himalaya. As Main Boundary Thrust (MBT) is approached, the attribute value increases further. A quasi-probabilistic seismic hazard index has been proposed based on site response, peak ground acceleration and predominant frequency. Six seismic hazard zones are marked with percent probability <22%, 22–37%, 37–52%, 52–67%, 67–82%, >82% at 3 Hz and <20%, 20–34%, 34–48%, 48–61%, 61–75%, >75% at 9 Hz. In the microzonation vector layer of integrated seismological and geological themes also six major zones are mapped, with percent probability <15%, 15–31%, 31–47%, 47–63%, 63–78%, >78% at low frequency end. The maximum risk is attached to the probability greater than 78% in the Singtam and its adjoining area. These maps are generally better spatial representation of seismic hazard including site-specific analysis.     

Site-specific earthquake response study for hazard assessment in Kolkata city,India

The assessment of local site effects on seismic ground motions is of great importance in earthquake engineering practice. Several destructive earthquakes in the past have demonstrated that the amplification of ground motion and associated damage to structures due to local site conditions is a significant consideration in earthquake hazard analysis. A recent paper published in this journal highlights the hazard posed by earthquakes in the megacity of Kolkata in India due to its seismic and geological settings. The seismic hazard assessment study speculates that the deep alluvial deposit in the city may increase the seismic hazard probably due to the amplification of the seismic energies. This paper focuses on the seismic response studies of the various soil strata (i.e. for local subsurface conditions) obtained from various construction sites in the city for predicted earthquake. It is very well recognized that site response studies (a part of seismic microhazard zonation for urban areas) are the first step towards performance-based foundation design or seismic risk analysis and mitigation strategy. One of the problems for carrying out site-specific study in Kolkata is the lack of recorded strong motion data in the city. Hence, this paper outlines a methodology to carry out site-specific study, where no strong motion data or seismic data are available. The methodology uses wavelet-based spectrum compatibility approach to generate synthetic earthquake motions and equivalent linear method for seismic site response analysis. The Mega City of Kolkata has been considered to explain the methodology. Seismic hazard zonation map by the Bureau of Indian Standards classifies the City of Kolkata as moderate seismic zone (Zone III) with a zone factor 0.16. On the other hand, GSHAP(Global Seismic Hazard Assessment Program) map which is based on 10% probability of exceedance in 50 years specifies a maximum peak ground acceleration (PGA) of 1.6 m/s² (0.163 g) for this region. In the present study, the seismic response has been carried out based on GSHAP. The results of the analysis indicate the amplification of ground motion in the range of 4.46–4.82 with the fundamental period ranging from 0.81 to 1.17 s. Furthermore, the maximum spectral accelerations vary in the range of 0.78–0.95 g.     

Thermal models of arcs and ridges — A “Source-span-sink” model

A theory of two-dimensional geothermic problems is elaborated by the active temperature function at the vertical contact of two horizontally layered media. The approach offered before for oceanic ridges is extended to the case of continental margins and the upper part of a descending slab, i.e. “sink”, in island-arc areas. It is assumed that the plate motion in the oceanic area exists; in a descending area it is directed downward but remains zero on a continental side. Mathematically it symbolizes a “source—span—sink” thermal model. Numerical parameters are given for a theoretical thermal model of the heat-flow profile across the Kuril island arc, from the trench through Iturup Island, Sakhalin Island and the Tatarian Trough.     

微差爆破模拟天然地震的数值分析与效果评价

微差爆破可以创造出与天然地震相类似的振动环境。为了更好地对爆破模拟天然地震进行研究和评价,在LS-DYNA软件框架内,建立了适合于爆炸荷载作用下土体动力分析的实用模型,对原位土体场地爆破地震动试验进行了全尺寸数值模拟,并将模拟结果与现场实测数据进行比较分析,验证数值计算的准确性。计算结果表明,模拟得到的场地质点加速度-时程与实测数据比较吻合,因此,可以预先利用数值手段为试验方案设计提供指导;同时,数值模拟还总结了现场试验中未得到的场地压力变化规律和振动规律,分析了爆破地震动与天然地震的异同,丰富了爆破地震动的研究,表明利用爆破手段能够近似地模拟天然地震,为微差爆破模拟天然地震的研究工作提供了经济、简便而实用的方法。     

近断层速度脉冲地震动引起的边坡永久位移预测模型

基于小波分析方法,从NGA数据库的3 551条地震记录中选取189条速度脉冲地震动,地震动均转换成发生最强脉冲的方向。基于Newmark方法,分析了近断层速度脉冲地震动作用引起的边坡永久位移值。结果表明：近断层速度脉冲地震动对边坡产生特殊的破坏作用,表现在滑动位移值大、滑动体破坏力强等方面;边坡永久位移值与速度脉冲地震动的峰值速度具有高度相关性,位移值较大时尤为明显。建立了基于单变量形式的峰值速度及双变量形式的峰值速度、峰值加速度两种边坡永久位移预测模型,模型简单实用,与回归数据具有很好的相关性,前者更适用于预测对实际工程影响较大的永久位移值,且离散性较小。提出的预测模型为考虑近断层地震动速度脉冲特性影响的边坡永久位移值的概率地震灾害分析提供了基础。     

The externality fields of football grounds: a case study of the Dell, Southampton

Football grounds generate a negative externality — or nuisance— effect which, according to previous work, displays a distance decay pattern. However, on the basis of an exploratory investigation of the negative externality field produced by Southampton Football Club's ground it was found that the various elements which together comprise the general nuisance field varied both in intensity and in spatial extent. The noise nuisance, not surprisingly, is the most localized whilst the parked cars and traffic nuisances are the most extensive. In addition, whereas the specific nuisances of noise, pedestrians and ‘hooliganism’ each exhibit a distance decay pattern the incidence of the parked cars nuisance and, to a lesser extent, the traffic nuisance are less in the immediate vicinity of the ground than further away. Finally, despite the impression conveyed by the media, ‘hooliganism’ was perceived by local residents as much less of a football-generated nuisance than traffic congestion and parked cars.     

Perceptible earthquakes in the broad Aegean area

A probabilistic estimate of seismic hazard can be obtained from the spatial distribution, of earthquake sources, their frequency–magnitude distribution and the rate of attenuation of strong ground motion with distance. We calculate the earthquake perceptibility, i.e. the annual probability that a particular level of ground shaking will be generated by earthquakes of particular magnitude, by weighting frequency–magnitude data with the predicted felt area for a given level of ground shaking at a particular magnitude. This provides an earthquake selection criterion that can be used in the anti-seismic design of non-critical structures. We calculate the perceptibility, at a particular value of isoseismal intensity, peak ground acceleration and velocity, as a function of source magnitude and frequency for the broad Aegean area using local attenuation laws. We use frequency–magnitude distributions that were previously obtained by combining short-term catalogue data with tectonic moment rate data for 14 tectonic zones in Greece with sufficient earthquake data, and where contemporary strain rates are available from satellite data. Many of the zones show a ‘characteristic earthquake’ distribution with the most perceptible earthquake equal to the maximum magnitude earthquake, but a relatively flat perceptibility between magnitudes 6 and 7. The maximum perceptible magnitude is in the fastest-deforming region in the middle of the Aegean sea, and tends to be systematically low on the west in comparison to the east of the Aegean sea. The tectonic data strongly constrain the long-term recurrence rates and lead to low error estimates (±0.2) in the most perceptible magnitudes.     

The role of the Kazakhstan orocline in the late Paleozoic amalgamation of Eurasia

After the 2005 Kashmir earthquake, we mapped surface ground fractures in Tangdhar, Uri, Rajouri and Punch sectors and liquefaction features in Jammu area lying close to the eastern side of the Line of Control (LOC) in Kashmir, India. The NW trending ground fractures occurred largely in the hanging wall zone of the southeastern extension of the causative fault in Tangdhar and Uri sectors. The principal compressive stress deduced from the earthquake induced ground fractures is oriented at N10°, whereas the causative Balakot–Bagh fault strikes 330°. The fault-plane solution indicates primarily SW thrusting of the causative fault with a component of strike–slip motion. The ground fractures reflect pronounced strike–slip together with some tensile component. The Tangdhar area showing left-lateral strike–slip motion lies on the hanging wall, and the Uri region showing right-lateral strike–slip movement is located towards the southeastern extension of the causative fault zone. The shear fractures are related to static stress that was responsible for the failure of causative fault. The tensile fractures with offsets are attributed to combination of both static and dynamic stresses, and the fractures and openings without offsets owe their origin due to dynamic stress. In Punch–Rajouri and Jammu area, which lies on the footwall, the fractures and liquefactions were generated by dynamic stress. The occurrence of liquefaction features in the out board part of the Himalayan range front near Jammu is suggestive of stress transfer  230 km southeast of the epicenter. The Balakot–Bagh Fault (BBF), the Muzaffarabad anticline, the rupture zone of causative fault and the zone of aftershocks — all are aligned in a  25 km wide belt along the NW–SE trending regional Himalayan strike of Kashmir region and lying between the MBT and the Riasi Thrust (Murree Thrust), suggesting a seismogenic zone that may propagate towards the southeast to trigger an earthquake in the eastern part of the Kashmir region.     

Estimation of regionalized compositions: A comparison of three methods

A regionalized composition is a random vector function whose components are positive and sum to a constant at every point of the sampling region. Consequently, the components of a regionalized composition are necessarily spatially correlated. This spatial dependence—induced by the constant sum constraint—is a spurious spatial correlation and may lead to misinterpretations of statistical analyses. Furthermore, the cross-covariance matrices of the regionalized composition are singular, as is the coefficient matrix of the cokriging system of equations. Three methods of performing estimation or prediction of a regionalized composition at unsampled points are discussed: (1) the direct approach of estimating each variable separately; (2) the basis method, which is applicable only when a random function is available that can he regarded as the size of the regionalized composition under study; (3) the logratio approach, using the additive-log-ratio transformation proposed by J. Aitchison, which allows statistical analysis of compositional data. We present a brief theoretical review of these three methods and compare them using compositional data from the Lyons West Oil Field in Kansas (USA). It is shown that, although there are no important numerical differences, the direct approach leads to invalid results, whereas the basis method and the additive-log-ratio approach are comparable.     

Effects of open gaps on particle velocity measurements

Summary Joints and bedding planes are present in nearly all geologic media. These discontinuities vary in the extent to which they are open and the material with which they are filled. The presence of such discontinuities affect the ground motions that result from the detonation of an explosive source. The effects of these discontinuities are especially important if it is desired to predict the size of the explosive source from measurements of ground motion — such as particle velocity or stress. In this paper we present results from a series of model tests conducted to investigate the effect of open gaps on the particle velocity that can be transmitted across them. Electromagnetic velocity gages were used to record particle velocities at several locations on either side of an open gap of known width. The study was restricted to the transmission of particle velocities normal to the gap and to the case where the gap was located in a region for which the stress was high but not greatly above the elastic limit of the material. The experimental results were compared qualitatively to results obtained from WONDY — a one dimensional Lagrangian finite difference code. Although the tests were conducted in laboratory models, the results can be scaled to compare to the full scale situation.     

Probabilistic Analysis of Ground Surface Vibration Due to Train Movement, a Case Study on Tehran Metro Line 4

Ground surface vibration produced by moving train is one of the most important aspects in urban areas. The purpose of this study is the probabilistic analysis of ground surface vibration, which is generated by metro transportations. For this reason, Tehran metro line 4 is considered as a case study. In this paper, at first, a new procedure is used to simulate train dynamic load. In the second step, based on the variation of geomechanical properties and train characteristics in Tehran metro line 4, more than 60 numerical models are simulated. The results of numerical simulations are analyzed by multivariate statistical technique and an equation for prediction of peak particle velocity (PPV) in the ground surface is presented. In the next step, probabilistic analysis is done using Monte Carlo Simulation (MCS). Finally, sensitivity of input data on ground surface vibration is discussed and the impact of geomechanical properties and train characteristics on the surface vibration is considered. Based on the probabilistic analysis, PPV in the surface region of Tehran metro line 4 is <2.76?mm/s with 95?% probability.     

Seismic hazard and site-specific ground motion for typical ports of Gujarat

Economic importance of major ports is well known, and if ports are located in seismically active regions, then site-specific seismic hazard studies are essential to mitigate the seismic risk of the ports. Seismic design of port sites and related structures can be accomplished in three steps that include assessment of regional seismicity, geotechnical hazards, and soil structure interaction analysis. In the present study, site-specific probabilistic seismic hazard analysis is performed to identify the seismic hazard associated with four typical port sites of Gujarat state (bounded by 20°–25.5°N and 68°–75°E) of India viz. Kandla, Mundra, Hazira, and Dahej ports. The primary aim of the study is to develop consistent seismic ground motion for the structures within the four port sites for different three levels of ground shaking, i.e., operating level earthquake (72 years return period), contingency level earthquake (CLE) (475 year return period), and maximum considered earthquake (2,475 year return period). The geotechnical characterization for each port site is carried out using available geotechnical data. Shear wave velocities of the soil profile are estimated from SPT blow counts using various empirical formulae. Seismicity of the Gujarat region is modeled through delineating the 40 fault sources based on the seismotectonic setting. The Gujarat state is divided into three regions, i.e., Kachchh, Saurashtra, and Mainland Gujarat, and regional recurrence relations are assigned in the form of Gutenberg-Richter parameters in order to calculate seismic hazard associated with each port site. The horizontal component of ground acceleration for three levels of ground shaking is estimated by using different ground motion attenuation relations (GMAR) including one country-specific GMAR for Peninsular India. Uncertainty in seismic hazard computations is handled by using logic tree approach to develop uniform hazard spectra for 5% damping which are consistent with the specified three levels of ground shaking. Using recorded acceleration time history of Bhuj 2001 earthquake as the input time motion, synthetic time histories are generated to match the developed designed response spectra to study site-specific responses of port sites during different levels of ground shaking. It is observed that the Mundra and Kandla port sites are most vulnerable sites for seismic hazard as estimated CLE ground motion is in order of 0.79 and 0.48 g for Mundra and Kandla port sites, respectively. Hazira and Dahej port sites have comparatively less hazard with estimated CLE ground motion of 0.17 and 0.11 g, respectively. The ground amplification factor is observed at all sites which ranges from 1.3 to 2.0 for the frequency range of 1.0–2.7 Hz. The obtained spectral accelerations for the three levels of ground motions and obtained transfer functions for each port sites are compared with provisions made in Indian seismic code IS:1893-Part 1 (2002). The outcome of present study is recommended for further performance-based design to evaluate the seismic response of the port structures with respect to various performance levels.     

Mineral variation in anorthositic, troctolitic, and adamellitic rocks of the Barth Island layered structure in the Nain anorthosite complex, Labrador

The Barth Island layered structure is an oval, 6 by 9 km body, consisting of rhythmically layered adamellitic rock in the center which grades outward through jotunite into troctolite. Farther outward the sequence repeats itself in reversed order, strongly reduced in magnitude and finer grained; the adamellitic zone is followed by jotunite which grades into coarse-grained leuconorite and into anorthosite of the Nain complex. The Barth Island structure, having an inverted conical base topped by a hemispherical depression, seems to represent a distorted sequence of rock layers with troctolite at the bottom, grading upward into adamellitic rocks which grade into anorthosite at the top. Trend-surface analysis demonstrates the regional variation of plagioclase and orthopyroxene compositions in the troctolite—adamellite sequence of the central part of the structure. The fits for the second- and third-degree surfaces are good and significant at the 99 percent level. The regression line for compositional variation in coexisting plagioclase and orthopyroxene in all analyzed rocks has a correlation coefficient of r² = 0.78. The difference between the trends in the troctolite—adamellite sequence and the anorthosite—adamellite sequence is insignificant. The regression curve for compositional variation in coexisting orthopyroxene and olivine has a correlation coefficient of r² = 0.98. The curve shows good correlation with the experimentally established partitioning curve of Medaris, which indicates that equilibrium conditions prevailed during formation of the olivine—ortho-pyroxene pairs. The results suggest that the troctolite—adamellite sequence and the anorthosite—adamellite sequence are products of fractional crystallization, possibly from the same parental magma.     

Flood risk assessment in a coastal lagoon under present and future scenarios: Ria de Aveiro case study

Iranian strong motion records as well as detailed conditions of their instrument sites and the characteristics of their causative seismic sources are compiled and processed. The dataset consists of 2286 three-component records from 461 Iranian earthquakes with at least two high-quality records having moment magnitude from 3.9 to 7.3. These records are about 20% of the Iranian database and are suitable for seismic hazard analysis and engineering applications. Perhaps for the first time in the literature, the distance to the surface projection of the fault is reported for a great number of records corresponding to earthquakes with M > 6.0. The raw accelerations are processed using the wavelet de-noising method. Having corrected and filtered these raw data, the pseudospectral accelerations are calculated for each of the three components of time series, separately. In addition to the ground motion parameters, a large and comprehensive list of metadata characterizing the recording conditions of each record is also developed. Moreover, careful revision of the characteristics of the earthquakes such as location, magnitude, style of faulting and fault rupture plane geometry, if available, is carried out using the best available information in a scientifically sound manner. Finally, we also focus on special ground motion records including records with peak ground acceleration (PGA) >300 cm/s² and distances less than 30 km. These are “exceptional” records in the Iranian dataset and include less than 2% of the selected dataset.     

Ground Vibration due to Blasting in Limestone Quarries

A database of ground vibration due to blasting at 27 limestone quarries, located in various parts of India, has been created. The database contains peak particle velocity (PPV), frequency, other vibration related and blast design parameters. Regression analysis of the data is carried out to derive site constants of the USBM predictor equation for individual quarries. It is found that these site constants are correlated with each other. By combining all the data, a generalised predictor equation is developed to assess and control ground vibration. In addition, mean zone of attenuation has been delineated using the predictor equations of the individual quarries. The dominant frequency of ground vibration with respect to distance and the possibility of modifying it by changing delay intervals in production blasts are also examined.