Soil-structure interaction effects on an instrumented building

The paper highlights the use of fem and bi-directional lumped-mass-storey-stiffness numerical models for the study of the soil–structure interaction (ssi) effects on an instrumented building. Data on the structural response have been obtained through the project for seismic instrumentation of a 16-storey r/c cast-in-place dwelling building (Chisinau, Republic of Moldova) during a series of earthquakes (Gutenberg–Richter M _GR = 5.0−6.7). The effect of soil–structure interaction is clearly observed comparing the responses recorded on foundation and free-field. ssi becomes more pronounced for higher level of ground shaking amplifying the natural period of the structure and slightly suppressing high frequences on the foundation in comparison with the free-field motion.     

Seismic response of bridges with massive foundations

The dynamic response of vehicular overpasses with massive foundations built in highly populated earthquake prone regions is studied to assess the massive foundation potential of being a technically sound mean to reduce the structural response during major earthquakes. The study consists on numerical simulations using 3-D finite element models. Two typical supports of a major 23 km long vehicular overpass recently built in the north east part of Mexico City valley were considered in this research. This zone is characterized by the presence of stiff soils comprised by dense and very dense silty sands and sandy silts, randomly intercalated by stiff clays layers of variable thickness. Initially, a conventional raft foundation structurally connected to four precast closed-end concrete piles was considered. Then, the potential beneficial effect of a massive foundation of variable depth was assessed. Sets of 3-D finite element models were developed and the response of the systems was evaluated for a typical seismic scenario such as that prevailing at the zone, assuming a potential 8.1 M_w seismic event, and for a hypothetical 8.7 M_w extreme event. Important attenuations of about 39% to 48% at the upper deck spectral accelerations and of 21% to 30% in maximum lateral foundation displacements were achieved with the massive soil improvement for the cases analyzed. Thus, the massive foundations seem to be a convenient alternative to reduce the overall structural seismic response.     

Dimensional analysis of structures with translating and rocking foundations under near-fault ground motions

This paper evaluates the inertial soil–structure interaction (SSI) effects on linear and bilinear structures supported on foundation that is able to translate and rock when subject to near-fault ground motions. Through rigorous dimensional analysis, the peak structural responses (e.g. structural drift and total acceleration) of the soil–foundation–structure interacting (SFSI) systems are characterized by a set of dimensionless Π-parameters, which can decisively describe the interactive behavior of SFSI systems. By comparing the normalized structural responses of various structure–foundation systems with their fixed-base counterparts, the study reveals that SSI effects highly depend on the structure-to-pulse frequency ratio, Π_ω, the foundation-to-structure stiffness ratio, Π_k, damping coefficient of foundation impedance, Π_c, the foundation rocking, and the development of nonlinearity in structures. For linear structures, the SSI effects are insignificant when the structure-to-pulse frequency ratio (Π_ω) is smaller than 1.5 and can amplify the structural responses when Π_ω is higher than 1.5. Foundation rocking can shift and enlarge the response amplification zone of SSI. For nonlinear structures, SSI tends to reduce the structural responses for Π_ω<3 while can increase the ductility demands for Π_ω≥3. The bilinear structures may experience more significant SSI effects than linear structures in certain frequency ranges. The numerical simulations on ten real building cases exhibiting significant rocking and a detailed case study on a nine-story frame structure demonstrate the applicability of dimensional analysis results to predict the SSI effects on realistic building structures. The study demonstrates that the dimensional analysis provides a concise and systematic way of evaluating SSI effects.     

Activation of seismicity in Central and South Asia after the Makran earthquakes: Possible acceleration of preparation of large seismic events in the Tien Shan region

Two zones of seismicity (ten events with M _w = 7.0–7.7) stretching from Makran and the Eastern Himalaya to the Central and EasternTien Shan, respectively, formed over 11 years after the great Makran earthquake of 1945 (M _w = 8.1). Two large earthquakes (M _w = 7.7) hit theMakran area in 2013. In addition, two zones of seismicity (M ≥ 5.0) occurred 1–2 years after theMakran earthquake in September 24, 2013, stretching in the north-northeastern and north-northwestern directions. Two large Nepal earthquakes struck the southern extremity of the “eastern” zone (April 25, 2015, M _w = 7.8 and May 12, 2015, M _w = 7.3), and the Pamir earthquake (December 7, 2015, M _w = 7.2) occurred near Sarez Lake eastw of the “western” zone. The available data indicate an increase in subhorizontal stresses in the region under study, which should accelerate the possible preparation of a series of large earthquakes, primarily in the area of the Central Tien Shan, between 70° and 79° E, where no large earthquakes (M _w ≥ 7.0) have occurred since 1992.     

Quantification of major earthquake tsunamis of the Japan Sea

The size of major tsunamigenic earthquakes which occurred in the Japan Sea is quantified on the basis of seismic and tsunamigenic source parameters. The tsunami magnitude M_t is determined from the instrumental tsunami-wave amplitudes. The M_t values thus obtained are on average 0.2 units larger than the values of moment magnitude M_w, though the M_t scale has originally been adjusted to agree with M_w. Moreover, the volume of displaced water at the source is on average 2.3 times as large as that for the Pacific events with a comparable M_w. Nevertheless, the observed height of the sea-level disturbance at the source is found consistent with the amount of crustal deformation computed for the seismic fault models. These results indicate that the tsunami source potential itself is large for M_w in comparison with the Pacific events. The large source potential is explained in terms of the effective difference both in the rigidity of the source medium and in the geometry of the fault motion. For the Japan Sea events, the M_t scale still provides the physical measure of the tsunami potential, and M_t minus 0.2 corresponds to M_w. This predicts that the maximum amplitude of tsunami waves from Japan Sea earthquakes is at least two times as large as that from Pacific earthquakes with a comparable M_w.     

Seismic response evaluation of an urban overpass

The seismic response of one section of a 23 km strategic urban overpass to be built in the so‐called transition and hill zones in Mexico City is presented. The subsoil conditions at these zones typically consist on soft to stiff clay and medium to dense sand deposits, randomly interbedded by loose sand lenses, and underlain by rock formations that may outcrop in some areas. Several critical supports of this overpass are going to be instrumented with accelerometers, inclinometers and extensometers, tell tales and end pile cell pressures to assess their seismic performance during future earthquakes and to generate a database to calibrate soil–structure interaction numerical models. This paper presents the seismic performance evaluation of the critical supports located in one section of the overpass. Sets of finite elements models of the soil–foundation–structure systems were developed. Initially, the model was calibrated analyzing the seismic response that an instrumented bridge support exhibited during the June 15th, 1999 Tehuacan (M_w = 7) Earthquake. This bridge is located also within the surroundings of Mexico City, but in the lake zone, where highly compressible clays are found. The computed response was compared with the measured response in the free field, pile‐box foundation and bridge deck. Once the model prediction capabilities were established, the seismic response of the critical supports of the urban overpass was evaluated for the design earthquake in terms of transfer functions and displacement time histories. Copyright © 2010 John Wiley & Sons, Ltd.     

Earthquake ground motion predictive equations for Garhwal Himalaya,India

Predictive equations based on the stochastic approach are developed for earthquake ground motions from Garhwal Himalayan earthquakes of 3.5≤M_w≤6.8 at a distance of 10≤R≤250 km. The predicted ground motion parameters are response spectral values at frequencies from 0.25 to 20 Hz, and peak ground acceleration (PGA). The ground motion prediction equations (GMPEs) are derived from an empirically based stochastic ground motion model. The GMPEs show a fair agreement with the empirically developed ground motion equations from Himalaya as well as the NGA equation. The proposed relations also reasonably predict the observed ground motion of two major Himalayan earthquakes from Garhwal Himalayan region. For high magnitudes, there is insufficient data to satisfactorily judge the relationship; however it reasonably predicts the 1991 Uttarkashi earthquake (M_w=6.8) and 1999 Chamoli earthquake (M_w=6.4) from Garhwal Himalaya region.     

Deep structure and tomographic imaging of strong earthquake source zones

This work generalizes the results of tomographic imaging performed by the authors for epicentral zones. Seismic events in North Africa (the M _w = 5.8 earthquake of 1985 near the town of Constantine), eastern Anatolia (the Erzincan M _w = 6.7 earthquake of 1992), the Lesser and Greater Caucasus (the 1988 Spitak M _w = 6.8 and the 1991 Racha M _w = 7.0 earthquakes), and northern Sakhalin (the 1995 Neftegorsk M _w = 7.1 earthquake) are examined. It is shown how various morphokinematic types of active faults differ in the resulting tomographic images at various depths. A classification of tomographic images of strong earthquake source zones is proposed in accordance with the rank of their generating faults. The sources of the Spitak, Racha, and Erzincan earthquakes are confined to large boundary faults separating tectonic zones. Lower velocity bands are revealed in the tomographic images, and low velocity “pockets” 1–2 km or somewhat more in width penetrating to a depth of up to 15 km are observed near the fault zones. The Constantine and Neftegorsk earthquakes were generated by faults of a lower rank. The source zones of these events are imaged tomographically as narrow gradient zones.     

Patterns of significant seismic quiescence on the Mexican Pacific coast

Many authors have proposed that the study of seismicity rates is an appropriate technique for evaluating how close a seismic gap may be to rupture. We designed an algorithm for identification of patterns of significant seismic quiescence by using the definition of seismic quiescence proposed by Schreider (1990). This algorithm shows the area of quiescence where an earthquake of great magnitude may probably occur. We have applied our algorithm to the earthquake catalog on the Mexican Pacific coast located between 14 and 21 degrees of North latitude and 94 and 106 degrees West longitude; with depths less than or equal to 60 km and magnitude greater than or equal to 4.3, which occurred from January, 1965 until December, 2014. We have found significant patterns of seismic quietude before the earthquakes of Oaxaca (November 1978, M_w = 7.8), Petatlán (March 1979, M_w = 7.6), Michoacán (September 1985, M_w = 8.0, and M_w = 7.6) and Colima (October 1995, M_w = 8.0). Fortunately, in this century earthquakes of great magnitude have not occurred in Mexico. However, we have identified well-defined seismic quiescences in the Guerrero seismic-gap, which are apparently correlated with the occurrence of silent earthquakes in 2002, 2006 and 2010 recently discovered by GPS technology.     

Seismic monitoring of source zones of strong earthquakes

A rational choice of the scalar seismic moment and ordering index is proposed that can be advantageously used for the monitoring of source zones of strong earthquakes in order to predict the development of a seismic situation. These parameters are the main characteristics of seismotectonic deformation. The ordering index characterizes a regular change in time of chaotization and ordering phases of the seismic process related to the occurrence of strong aftershocks. Using the December 5, 1997, Kronotskii (M _w = 7.8) and December 26, 2004, Sumatra (M _w = 9.0) earthquakes as an example, temporal variations of the studies parameters in the aftershock zones of these earthquakes are analyzed in detail.     

Ambient vibration tests of a seven-story reinforced concrete building in Van Nuys, California, damaged by the 1994 Northridge earthquake

Results of two detailed ambient vibration surveys of a 7-story reinforced concrete building in Van Nuys, California, are presented. Both surveys were conducted after the building was severely damaged by the 17 January 1994, Northridge earthquake (M_L=5.3, epicenter 1.5 km west from the building site) and its early aftershocks. The first survey was conducted on 4 and 5 February 1994, and the second one on 19 and 20 April 1994, about one month after the 20 March aftershock (M_L=5.3, epicenter 1.2 km north–west from the building site). The apparent frequencies and two- and three-dimensional mode shapes for longitudinal, transverse and vertical vibrations were calculated. The attempts to detect the highly localized damage by simple spectral analyses of the ambient noise data were not successful. It is suggested that very high spatial resolution of recording points is required to identify localized column and beam damage, due to the complex building behavior, with many interacting structural components. The loss of the axial capacity of the damaged columns could be seen in the vertical response of the columns, but similar moderate or weak damage typically would not be noticed in ambient vibration surveys. Previous analysis of the recorded response of this building to 12 earthquakes suggests that, during large response of the foundation and piles, the soil is pushed sideways and gaps form between the foundation and the soil. These gaps appear to be closing during “dynamic compaction” when the building site is shaken by many small aftershocks. The apparent frequencies of the soil–foundation–structure system appear to be influenced significantly by variations in the effective soil–foundation stiffness. These variations can be monitored by a sequence of specialized ambient vibration tests.     

Source characteristics of the deep Philippine earthquake cluster, 23 and 24 July 2010

We determine the rupture velocity, rupture area, stress drop and duration of four strong deep-focus earthquakes in the Philippines by back-projecting the teleseismic P waves. Four deep-focus earthquakes occurred in a totally consumed Molucca microplate; their focal depths were greater than 550 km and their moment magnitudes were between M _w 6.6 and M _w 7.6. By studying this deep-focus cluster, we are able to estimate the rupture velocity, rupture area and stress drop which would assist in constraining the physical mechanism for earthquakes deeper than 500 km. Since the Molucca microplate is totally consumed, little evidence is left on the surface for us to do research. This deep-focus cluster provides us the opportunity to reveal the properties of this totally consumed microplate by using seismic method for the first time. Four earthquakes in this deep-focus cluster all have multiple rupture subevents. The M _w 7.3 event ruptures in two subevents, the M _w 7.6 and M _w 7.4 events both have three subevents. The M _w 6.6 event has single peak on the amplitude as a function of time; however, its energy releases at two spatially separated areas. Our results show that this deep-focus cluster has a slow rupture velocity which is about 0.27 to 0.43 of the shear wave velocity, long-scaled duration, concentrated energy release area, and high stress drop. These source properties are similar to those of other deep earthquakes occurring in warm slabs and indicate that the totally consumed Molucca microplate possibly is a warm plate.     

Seismic behavior of a post-tensioned integral bridge including soil–structure interaction (SSI)

Current practice usually pays little attention to the effect of soil–structure interaction (SSI) on seismic analysis and design of bridges. The objective of this research study is to assess the significance of SSI on the modal with geometric stiffness and seismic response of a bridge with integral abutments that has been constructed using a new bridge system technology. Emphasis is placed on integral abutment behavior, since abutments together with piers are the most critical elements in securing the integrity of bridge superstructures during earthquakes. Comparison is made between analytical results and field measurements in order to establish the accuracy of the superstructure–abutment model. Sensitivity studies are conducted to investigate the effects of foundation stiffness on the overall dynamic and seismic response of the new bridge system.     

Radiated Energy of Great Earthquakes from Teleseismic Empirical Green’s Function Deconvolution

We expand on the empirical Green’s function deconvolution method of Ide et al. (2011) to estimate radiated energy for the six largest earthquakes worldwide over the last 10 years: 2011 M _w 9.0 Tohoku-Oki, 2004 M _w 9.1 Sumatra, 2010 M _w 8.8 Maule, 2005 M _w 8.7 Nias, 2007 M _w 8.5 Bengkulu, and 2012 M _w 8.6 off-Sumatra. Deconvolution of P, SV and SH components gives consistent energy results that are comparable to estimates found independently by other researchers. Apparent stress for the five great thrust earthquakes is between 0.4 and 0.8 MPa, while the 2012 off-Sumatra strike-slip earthquake has a higher apparent stress of 3 MPa, which is consistent with other studies that find a tendency for strike-slip events to be more energetic. Our results are within the spread of apparent stress from the wider global earthquake population over a large magnitude range. The azimuthal distribution of energy in each case shows signs of directivity, and in some cases, shows less energy radiated in the trench-ward direction, which may suggest enhanced tsunami potential. We find that eGfs as small as ~M 6.5 can be used for teleseismic deconvolution, and that an eGf-mainshock magnitude difference of 1.5 units yields stable results. This implies that M 8 is the minimum mainshock size for which teleseismic eGf deconvolution will work well. We propose that a database of eGf events could be used to calculate radiated energy and apparent stress of great, hazardous events in near real time, i.e., promptly enough that it could contribute to rapid response measures.     

The European-Mediterranean Earthquake Catalogue (EMEC) for the last millennium

The catalogue by Grünthal et al. (J Seismol 13:517?C541, 2009a) of earthquakes in central, northern, and north-western Europe with M _w????3.5 (CENEC) has been expanded to cover also southern Europe and the Mediterranean area. It has also been extended in time (1000?C2006). Due to the strongly increased seismicity in the new area, the threshold for events south of the latitude 44°N has here been set at M _w????4.0, keeping the lower threshold in the northern catalogue part. This part has been updated with data from new and revised national and regional catalogues. The new Euro-Mediterranean Earthquake Catalogue (EMEC) is based on data from some 80 domestic catalogues and data files and over 100 special studies. Available original M _w and M ₀ data have been introduced. The analysis largely followed the lines of the Grünthal et al. (J Seismol 13:517?C541, 2009a) study, i.e., fake and duplicate events were identified and removed, polygons were specified within each of which one or more of the catalogues or data files have validity, and existing magnitudes and intensities were converted to M _w. Algorithms to compute M _w are based on relations provided locally, or more commonly on those derived by Grünthal et al. (J Seismol 13:517?C541, 2009a) or in the present study. The homogeneity of EMEC with respect to M _w for the different constituents was investigated and improved where feasible. EMEC contains entries of some 45,000 earthquakes. For each event, the date, time, location (including focal depth if available), intensity I ₀ (if given in the original catalogue), magnitude M _w (with uncertainty when given), and source (catalogue or special study) are presented. Besides the main EMEC catalogue, large events before year 1000 in the SE part of the investigated area and fake events, respectively, are given in separate lists.     

Asperity Distribution of the 1952 Great Kamchatka Earthquake and its Relation to Future Earthquake Potential in Kamchatka

—The 1952 Kamchatka earthquake is among the largest earthquakes of this century, with an estimated magnitude of M _w = 9.0. We inverted tide gauge records from Japan, North America, the Aleutians, and Hawaii for the asperity distribution. The results show two areas of high slip. The average slip is over 3 m, giving a seismic moment estimate of 155×10²⁰Nm, or M _w = 8.8. The 20th century seismicity of the 1952 rupture zone shows a strong correlation to the asperity distribution, which suggests that the large earthquakes (M > 7) are controlled by the locations of the asperities and that future large earthquakes will also recur in the asperity regions.     

Critical evaluation of strong motion in Kocaeli and Düzce (Turkey) earthquakes

Turkey was struck by two major events on August 17th and November 12th, 1999. Named Kocaeli (M_w=7.4) and Düzce (M_w=7.2) earthquakes, respectively, the two earthquakes provided the most extensive strong ground motion data set ever recorded in Turkey. The strong motion stations operated by the General Directorate of Disaster Affairs, the Kandilli Observatory and Earthquake Research Institute of Bogazici University and Istanbul Technical University have produced at least 27 strong motion records for the Kocaeli earthquake within 200 km of the fault. Kocaeli earthquake has generated six motions within 20 km of the fault adding significantly to the near-field database of ground motions for M_w>=7.0 strike–slip earthquakes. The paper discusses available strong motion data, studies their attenuation characteristics, analyses time domain, as well as spectral properties such as spectral accelerations with special emphasis on fault normal and fault parallel components and the elastic attenuation parameter, kappa. A simulation of the Kocaeli earthquake using code FINSIM is also presented.     

Stress Drops of the 1997–1998 Colfiorito,Central Italy Earthquakes: Hints for a Common Behaviour of Normal Faults in the Apennines

Stress drop estimates of moderate-magnitude earthquakes in the Umbria–Marche region, in the northern Apennines, exhibit a large scatter. For the two M _w 5.7 and 6.0 main shocks of 26 September 1997 near Colfiorito, several papers resulted in stress drop estimates of 20 MPa, but values as low as 2–3 MPa were proposed as well. Also for the largest aftershocks (M _w > 4), estimates spread from < 1 MPa up to values ten times larger. We have critically revisited methods and data used in the literature. We have specifically faced the trade-off between source and propagation effects, as we believe that it is responsible for a part of the large scatter. To keep this trade-off under control, we have applied a methodology that combines the best fit of both source spectra after Empirical Green’s Function (EGF) deconvolution and observed ground motion spectra, finding that the results of the two different data sets converge independently at the same solution. We have used ground motions observed in the Colfiorito basin, where an accelerograph and a co-located seismological broad-band station recorded three clusters of earthquakes in a broad magnitude interval (1.7 ≤ M _w ≤ 6.0). We have found that the mainshock–aftershock sequences result in stress drops of 2–5 MPa at M _w ≥ 5.6, with an average tendency to decrease at smaller magnitudes where stress drop variability increases. These findings confirm the source scaling recently assessed through Empirical Green’s Function deconvolution for another well-monitored seismic sequence of normal-faulting earthquakes, which struck the city of L’Aquila in the central Apennines in April 2009. The similar scaling law of the two areas suggests common mechanisms of stress release for the shallow normal faults in the Apennines. The propensity of smaller earthquakes to increase in variability, with a tendency toward smaller stress drops, may reflect an effect of fault strength heterogeneities for smaller size ruptures.     

Physical size of tsunamigenic earthquakes of the northwestern Pacific

The new scale M_t of tsunami magnitude is a reliable measure of the seismic moment of a tsunamigenic earthquake as well as the overall strength of a tsunami source. This M_t scale was originally defined by Abe (1979) in terms of maximum tsunami amplitudes at large distances from the source. A method is developed whereby it is possible to determine M_t at small distances on the basis of the regional tsunami data obtained at 30 tide stations in Japan. The relation between log H, maximum amplitude (m) and log Δ, a distance of not less than 100 km away from the source (km) is found to be linear, with a slope close to 1.0. Using three tsunamigenic earthquakes with known moment magnitudes M_w, for calibration, the relation, M_t = log H + log Δ + D, is obtained, where D is 5.80 for single-amplitude (crest or trough) data and 5.55 for double-amplitude (crest-to-trough) data. Using a number of tsunami amplitude data, M_t is assigned to 80 tsunamigenic earthquakes that occurred in the northwestern Pacific, mostly in Japan, during the period from 1894 to 1981. The M_t values are found to be essentially equivalent to M_w for 25 events with known M_w. The 1952 Kamchatka earthquake has the largest M_t, 9.0. Of all the 80 events listed, at least seven unusual earthquakes which generated disproportionately-large tsunamis for their surface-wave magnitude M_s are identified from the relation. From the viewpoint of tsunami hazard reduction, the present results provide a quantitative basis for predicting maximum tsunami amplitudes at a particular site.