Bi-normalized response spectra and seismic intensity in Bucharest for 1986 and 1990 Vrancea seismic events

The Vrancea subcrustal earthquakes of August 30,1986 and May 30,1990 are the two most recent seismic events that have occurred in Romania with moment magnitudes M W ≥ 6.9.The spectral analysis of the strong ground motions recorded in Bucharest reveals that despite small differences in magnitude between the 1986 and 1990 earthquakes,their frequency contents are very different,sometimes even opposing.The main focus of this study is to conduct a comparative analysis of the response spectra in terms of the bi-normalized response spectra(BNRS) proposed by Xu and Xie(2004 and 2007) for strong ground motions recorded in Bucharest during these two seismic events.The mean absolute acceleration and relative velocity response spectra for the two earthquakes are discussed and compared.Furthermore,the mean bi-normalized absolute acceleration and normalized relative velocity response spectra with respect to the control period T C are computed for the ground motions recorded in Bucharest in 1986 and 1990.The predominant period T P is also used in this study for the normalization of the spectral period axis.Subsequently,the methodology proposed by Martinez-Perreira and Bommer(1998) is applied in order to estimate the seismic intensity of the two events.The results are discussed and several conclusions regarding the possibility of using the bi-normalized response spectra(BNRS) are given.     

高土石坝振动台模型试验理论、技术与应用

对土石坝振动台模型试验理论和技术进行系统阐述,提出基于原型和模型坝料静、动力特性试验的模型相似设计方法和不同强度地震动递进输入(白噪声微振-设计地震-校核地震-破坏试验)的振动试验方法。基于1g大型振动台和ng超重力离心机振动台设备性能现状,结合高土石坝的结构特点和动力试验相似模拟要求,对土石坝振动台模型试验的优势及局限进行深入讨论。结合已有的工程实践,对土石坝振动台模型试验在工程中的应用进行总结,并以某实际高面板堆石坝为例研究面板坝生命周期内经历多次地震情况下结构动力特性的演化规律。     

A time probabilistic approach to seismic landslide hazard estimates in Iran

Understanding where seismically induced landslides are most likely to occur is crucial in land use planning and civil protection actions aimed at reducing property damage and loss of life in future earthquakes. For this purpose an approach proposed by Del Gaudio et al. [1] has been applied to the whole Iranian territory to provide the basis to assess location and temporal recurrence of conditions of seismic activation of slope failures, according to the Newmark's model [2]. Following this approach, occurrence probabilities for different levels of seismic shaking in a time interval of interest (50 years) were first obtained through a standard hazard estimate procedure. Then, empirical formulae in the form proposed by Jibson et al. [3] and calibrated for the main seismogenic Iranian regions were used to evaluate the slope critical acceleration (Ac)_x for which a prefixed probability exists that, under seismic shakings, Newmark's displacement D_N exceeds a threshold×corresponding to landslide triggering conditions. The obtained (Ac)_x values represent the minimum slope resistance required to limit the probability of landslide seismic triggering within the prefixed value. A map reporting the spatial distribution of these values gives comparative indications on regional different exposure of slopes to shaking capable of inducing failures and provides a reference for hazard estimate at local scale. The obtained results show that the exposure to landslide seismic induction is maximum in the Alborz Mountains region, where critical accelerations up to ∼0.1 g are required to limit the probability of seismic triggering of coherent type landslides within 10% in 50 years.     

A substructure shaking table test for reproduction of earthquake responses of high‐rise buildings

When subjected to long‐period ground motions, high‐rise buildings' upper floors undergo large responses. Furniture and nonstructural components are susceptible to significant damage in such events. This paper proposes a full‐scale substructure shaking table test to reproduce large floor responses of high‐rise buildings. The response at the top floor of a virtual 30‐story building model subjected to a synthesized long‐period ground motion is taken as a target wave for reproduction. Since a shaking table has difficulties in directly reproducing such large responses due to various capacity limitations, a rubber‐and‐mass system is proposed to amplify the table motion. To achieve an accurate reproduction of the floor responses, a control algorithm called the open‐loop inverse dynamics compensation via simulation (IDCS) algorithm is used to generate a special input wave for the shaking table. To implement the IDCS algorithm, the model matching method and the H_∞ method are adopted to construct the controller. A numerical example is presented to illustrate the open‐loop IDCS algorithm and compare the performance of different methods of controller design. A series of full‐scale substructure shaking table tests are conducted in E‐Defense to verify the effectiveness of the proposed method and examine the seismic behavior of furniture. The test results demonstrate that the rubber‐and‐mass system is capable of amplifying the table motion by a factor of about 3.5 for the maximum velocity and displacement, and the substructure shaking table test can reproduce the large floor responses for a few minutes. Copyright © 2009 John Wiley & Sons, Ltd.     

Collapse assessment of moment frame buildings,considering vertical ground shaking

While many cases of structural damage in past earthquakes have been attributed to strong vertical ground shaking, our understanding of vertical seismic load effects and their influence on collapse mechanisms of buildings is limited. This study quantifies ground motion parameters that are capable of predicting trends in building collapse because of vertical shaking, identifies the types of buildings that are most likely affected by strong vertical ground motions, and investigates the relationship between element level responses and structural collapse under multi‐directional shaking. To do so, two sets of incremental dynamic analyses (IDA) are run on five nonlinear building models of varying height, geometry, and design era. The first IDA is run using the horizontal component alone; the second IDA applies the vertical and horizontal motions simultaneously. When ground motion parameters are considered independently, acceleration‐based measures of the vertical shaking best predict trends in building collapse associated with vertical shaking. When multiple parameters are considered, Housner intensity (SI), computed as a ratio between vertical and horizontal components of a record (SI_V/SI_H), predicts the significance of vertical shaking for collapse. The building with extensive structural cantilevered members is the most influenced by vertical ground shaking, but all frame structures (with either flexural and shear critical columns) are impacted. In addition, the load effect from vertical ground motions is found to be significantly larger than the nominal value used in US building design. Copyright © 2016 John Wiley & Sons, Ltd.     

Soil–pile-bridge structure interaction in liquefying ground using shake table testing

The evaluation of seismic pile response is particularly useful for geotechnical engineers involved in the design of foundations in liquefying site. Shake table testing was performed to study the dynamic interactive behavior of soil–pile foundations in liquefying ground under different shaking frequency and amplitude. The soil profile consisted of a clayey layer over liquefiable sand over clay. The model was tested with a series of El Centro earthquake motions with peak accelerations ranging from 0.15g to 0.50g, and time step from 0.006 to 0.02 s. Representative data, including time histories of accelerations and excess pore pressure ratios that characterize the important aspects of soil–pile interaction in liquefying ground are presented. The shaking frequency has no significant effect on the magnitudes of excess pore pressure ratio, ground and pile accelerations and pile bending moments. Excess pore pressure ratio, ground acceleration and pile acceleration, and pile bending moment largely depend on the shaking amplitude.     

Shaking table testing of multidrum columns and portals

A common type of ancient monuments around the Mediterranean is the ancient Greek temple. Unfortunately, very few remain intact; most of them surviving in the form of free‐standing multidrum columns. Composed of stones resting on top of each other without any connection, such columns are considered vulnerable to earthquakes. The paper presents an experimental study of such structures, aiming to explore their seismic vulnerability and derive insights on the key factors affecting their response. Reduced scale models of a single multidrum column and of a portal were tested at the shaking table of the National Technical University of Athens Laboratory of Soil Mechanics. The models, constructed of marble just as the originals, were excited by idealized Ricker pulses and real seismic records. Single columns exhibit a remarkable earthquake resistance. Subjected to the strongest motions ever recorded in Greece, where many such monuments are situated, the columns hardly suffered any permanent deformation. Collapse is probable only for extremely harsh directivity‐affected seismic motions. Portals proved even more robust, surviving extreme seismic excitations. Their superior performance is related to the beneficial role of the epistyle, which adds energy dissipation and restoring force to the system. Their performance is very sensitive to minor changes in geometry or input motion. The complexity increases exponentially with the number of drums, being directly associated with the number of drum‐to‐drum interfaces and the increased probability of interface imperfections. In contrast to PGA, the maximum spectral displacement SD_max and the length scale L_p have turned out to be effective intensity measures. Copyright © 2014 John Wiley & Sons, Ltd.     

Phasing issues in the seismic response of yielding,gravity-type earth retaining walls – Overview and results from a FEM study

An overview of past and recent developments on the subject of seismic earth pressures on yielding, gravity-type walls, retaining cohesionless backfill, is first presented, focusing on available data on the issue of phase difference that develops between the peak values of wall inertia and seismic earth thrust increment. The results of a FEM parametric study are next presented regarding the dependence on the resulting dynamic earth thrust reduction – acting on the time of peak wall inertia – on backfill rigidity, wall height, and shaking characteristics. The reliability of the numerical analyses was verified by modeling centrifuge tests reported by Nakamura [24] and successfully comparing measured vs. computed behavior. The results of the parametric analyses indicate that the seismic active earth thrust, acting on the wall at the time of maximum wall inertia, is significantly reduced (compared to its peak value) with increasing shaking intensity of backfill, increasing wall displacements, increasing wall height, and decreasing backfill rigidity. No systematic dependence on the ratio of input motion frequency to the natural frequency of the backfill (f/f₁) was observed. The above findings: (1) verify earlier experimental and numerical results, (2) explain the reported lack of damage to retaining walls under strong ground shaking, and (3) indicate the need for revising the pertinent provisions of current seismic codes. Graphs summarizing the results of the numerical analyses are presented which may be used as a guide for selecting the magnitude of seismic active earth thrust that needs to be taken into account in the design of the examined type of earth retaining walls.     

Shaking table test and numerical simulation of an RC frame‐core tube structure for earthquake‐induced collapse

The reinforced concrete frame‐core tube structure is a common form of high‐rise building; however, certain vertical components of these structures are prone to be damaged by earthquakes, debris flow, or other accidents, leaving no time for repair or retrofit. This study is motivated by a practical problem—that is, the seismic vulnerability and collapse resistant capability under future earthquakes when a vertical member has failed. A reduced scale model (1:15 scale) of a typical reinforced concrete frame‐core tube with a corner column removed from the first floor is designed, fabricated, and tested. The corner column is replaced by a jack, and the failure behavior is simulated by manually unloading the jack. The model is then excited by a variety of seismic ground motions on the shaking table. Experimental results concerning the seismic responses and actual process of collapse are presented herein. Finally, the earthquake‐induced collapse process is simulated numerically using the software program ANSYS/LS‐DYNA. Validation and calibration of the model are carried out by comparison with the experimental results. Furthermore, based on both experimental investigations and numerical simulations, the collapse mechanism is discussed, and some suggestions on collapse design are put forward. Copyright © 2016 John Wiley & Sons, Ltd.     

Shaking table comparative test and associated study of a stepped wall-frame structure

A new structural system called a stepped wall-frame structure is proposed in this study to solve the bottom yielding problem of RC frames, which widely occurred during previous earthquakes such as the Wenchuan and Yushu earthquakes in China. A 1/5 scale ordinary RC frame model and a stepped wall-frame model were subjected to shake table motions together to study the seismic behavior of the new structural system. This paper presents the dynamic characteristics, the seismic responses and the failure and collapse mechanism of the two models under low, moderate and high intensity shaking. The test results and further analysis demonstrate that the seismic performance of stepped wall-frame structures is superior to ordinary RC frames in terms of the well-controlled deformation pattern and more uniformly distributed damage. The stepped wall can effectively suppress the bottom yielding mechanism, and is simple, economical and practical for engineering practice.     

Performance of base‐isolated reinforced concrete buildings under bidirectional seismic excitation considering pounding with retaining walls including friction effects

Seismic pounding of base‐isolated buildings has been mostly studied in the past assuming unidirectional excitation. Therefore, in this study, the effects of seismic pounding on the response of base‐isolated reinforced concrete buildings under bidirectional excitation are investigated. For this purpose, a three‐dimensional finite element model of a code‐compliant four‐story building is considered, where a newly developed contact element that accounts for friction and is capable of simulating pounding with retaining walls at the base, is used. Nonlinear behavior of the superstructure as well as the isolation system is considered. The performance of the building is evaluated separately for far‐fault non‐pulse‐like ground motions and near‐fault pulse‐like ground motions, which are weighted scaled to represent two levels of shaking viz. the design earthquake (DE) level and the risk‐targeted maximum considered earthquake (MCE_R) level. Nonlinear time‐history analyses are carried out considering lower bound as well as upper bound properties of isolators. The influence of separation distance between the building and the retaining walls at the base is also investigated. It is found that if pounding is avoided, the performance of the building is satisfactory in terms of limiting structural and nonstructural damage, under DE‐level motions and MCE_R‐level far‐fault motions, whereas unacceptably large demands are imposed by MCE_R‐level near‐fault motions. In the case of seismic pounding, MCE_R‐level near‐fault motions are found to be detrimental, where the effect of pounding is mostly concentrated at the first story. In addition, it is determined that considering unidirectional excitation instead of bidirectional excitation for MCE_R‐level near‐fault motions provides highly unconservative estimates of superstructure demands. Copyright © 2014 John Wiley & Sons, Ltd.     

Collapse behavior of a brick masonry house using a shaking table and numerical simulation based on the extended distinct element method

Damage caused by devastating earthquakes has occurred in many developing countries. In order to mitigate such damage by promoting the study of adequate seismic design strategies, the authors conducted a dynamic collapse test on 3 m × 3 m × 3 m brick masonry house constructed with Pakistani bricks, using a one-direction horizontal large-scale shaking table. In order to analyze and simulate seismic performance of the masonry structures, the authors applied a new numerical simulating method based on the Extended Distinct Element Method (EDEM) and conducted collapse simulations of the brick masonry house behavior during the shaking table tests. In the numerical simulation model, bricks were assumed to be rigid bodies, and mortar was modeled using a mortar spring that consists of a normal spring and a shear spring. The parameters of each mortar spring were defined based on the results of material tests. Simulated results showed various collapsing processes, and the simulated aspects were found to be similar to the results of the shaking table tests.     

Effects of near-fault ground motions on the nonlinear behaviour of reinforced concrete framed buildings

The design provisions of current seismic codes are generally not very accurate for assessing effects of near-fault ground motions on reinforced concrete(r.c.)spatial frames,because only far-fault ground motions are considered in the seismic codes.Strong near-fault earthquakes are characterized by long-duration(horizontal)pulses and high values of the ratio α_(PGA)of the peak value of the vertical acceleration,PGA_V,to the analogous value of the horizontal acceleration,PGA_H,which can become critical for girders and columns.In this work,six- and twelve-storey r.c.spatial frames are designed according to the provisions of the Italian seismic code,considering the horizontal seismic loads acting(besides the gravity loads)alone or in combination with the vertical ones.The nonlinear seismic analysis of the test structures is performed using a step-by-step procedure based on a two-parameter implicit integration scheme and an initial stress-like iterative procedure.A lumped plasticity model based on the Haar-Karman principle is adopted to model the inelastic behaviour of the frame members.For the numerical investigation,five near-fault ground motions with high values of the acceleration ratio α_(PGA) are considered.Moreover,following recent seismological studies,which allow the extraction of the largest(horizontal) pulse from a near-fault ground motion,five pulse-type(horizontal)ground motions are selected by comparing the original ground motion with the residual motion after the pulse has been extracted.The results of the nonlinear dynamic analysis carried out on the test structures highlighted that horizontal and vertical components of near-fault ground motions may require additional consideration in the seismic codes.     

Population dynamics of earthquakes and mathematical modeling

We present a mathematical model that describes temporal variations of earthquakes. This model is represented as $$dn(t)/dt = n(t)\left[ {\alpha - \beta n(t) - \int_{ - \infty }^t {n(s)h(t - s)ds} } \right].$$ Heren(t) shows the numberof earthquakes per unit time in a certain region. α and β are constants. The functionh(t) denotes the hysteresis effect of the earthquake occurrences and can take the following forms depending on the physical conditions of the crusts; (A)h(t)=0: the equation represents a logistic type increase or decrease and approaches a stationary state asymptotically. This describes aftershock series of large earthquakes and earthquake swarms of large scale such as the Wakayama and Matsushiro swarms in Japan; (B)h(t)=constant (β=0): frequencyn(t) increases initially and then decreases gradually and shows some kind of volcanic swarms; (C)h(t) = κ · {exp(?γ₁ t) ? exp(γ₂ t)}, (γ₂ > γ₁): this denotes time delay effects and the model shows periodic patterns of bursts or “rhythms” of earthquakes, which are observed in earthquake swarms. When external effects are taken into consideration, the model is further generalized and can describe various seismic patterns. These effects represent various influences of the circumstances like the earth tide and fluctuations of plate motions, etc. Whenh(t) takes type (A) and the external effect is random, the equation displays repetitive random patterns with bursts. Particularly interesting cases may be those whenh(t) is type (C) and the external force is periodic like the earth tide. Various nonperiodic as well as periodic patterns of earthquakes appear. These are the phenomena of “chaos” and “entrainment”, etc. and can be commonly observed. Varieties of actual earthquake patterns seem to be, at least partly, explained by the nonlinear coupling between the tidal forces and autonomous rhythms of earthquakes.     

Seismic response of rocking isolated railway bridge piers with sacrificial components

In this study, sacrificial components were incorporated into self-centering railway bridge piers to improve the lateral stiffness. The seismic response of this new detail was investigated. First, the method to compute the initial uplift moment of the self-centering pier is given. In addition, shaking table tests were conducted on a free-rocking pier without sacrificial components, which was used to validate a two-spring numerical model. Good agreement was obtained between the numerical results and experimental data. Furthermore, the validated model was employed to investigate the influence of sacrificial components on the seismic response of rocking piers. For this purpose, two models were developed, with and without sacrificial components. Nonlinear response history analysis was then performed on both models under three historical motions. The results showed that compared to the one without sacrificial components, the rocking pier with sacrificial components has comparable displacement at the top of the pier, and maximum uplift moment at high amplitude motion. Therefore, incorporating sacrificial components into the rocking pier can increase the lateral stiffness at service load and low amplitude frequent earthquakes but can produce comparable response at high seismic excitation. These results provide support for performance-based seismic design of self-centering rocking piers.

     

Estimation of strong ground motions from hypothetical earthquakes on the Cascadia subduction zone,Pacific Northwest

Strong ground motions are estimated for the Pacific Northwest assuming that large shallow earthquakes, similar to those experienced in southern Chile, southwestern Japan, and Colombia, may also occur on the Cascadia subduction zone. Fifty-six strong motion recordings for twenty-five subduction earthquakes ofM _s7.0 are used to estimate the response spectra that may result from earthquakesM _w<81/4. Large variations in observed ground motion levels are noted for a given site distance and earthquake magnitude. When compared with motions that have been observed in the western United States, large subduction zone earthquakes produce relatively large ground motions at surprisingly large distances. An earthquake similar to the 22 May 1960 Chilean earthquake (M _w 9.5) is the largest event that is considered to be plausible for the Cascadia subduction zone. This event has a moment which is two orders of magnitude larger than the largest earthquake for which we have strong motion records. The empirical Green's function technique is used to synthesize strong ground motions for such giant earthquakes. Observed teleseismicP-waveforms from giant earthquakes are also modeled using the empirical Green's function technique in order to constrain model parameters. The teleseismic modeling in the period range of 1.0 to 50 sec strongly suggests that fewer Green's functions should be randomly summed than is required to match the long-period moments of giant earthquakes. It appears that a large portion of the moment associated with giant earthquakes occurs at very long periods that are outside the frequency band of interest for strong ground motions. Nevertheless, the occurrence of a giant earthquake in the Pacific Northwest may produce quite strong shaking over a very large region.     

土层平均剪切波速确定及其可靠性验证——基于钻井台阵强震动记录

震害资料显示,场地条件对地震动特性以及工程结构破坏程度影响显著。为减少因场地效应而造成的经济损失和社会影响,在进行场地地震反应分析时,需最大限度地减小因场地土层模型参数的不确定性引起的地震动评估偏差,为工程结构地震反应分析选取并生成适当的地震动输入。随着强震动观测技术的逐渐发展,大量可靠的钻井台阵记录为地震过程中场地观测点的动力反应提供了直接数据。以美国加州地区La Cienega钻井台阵强震动观测数据为基础,利用互相关函数,对不同强度地震作用下场地土层的平均剪切波速进行分析,并在此基础上,以Cyclic 1D为模拟平台,建立一维自由场地地震反应有限元分析模型。分析结果表明：通过钻井台阵地震动观测数据识别,得到场地平均剪切波速,能够反映该场地的动力特性,数值模拟计算结果和台阵地震动记录基本吻合,可为数值模型参数选取提供依据。     

The study of the radiation characteristics and propagation of seismic waves in the North Caucasus by modeling the accelerograms of the recorded earthquakes

For evaluating the parameters of the vibrations of the Earth’s surface in the case of strong earthquakes, which are possible in the future, the regular patterns of the emission and propagation of seismic waves in the North Caucasus regions are investigated. The regional parameters of emission and propagation of seismic waves are evaluated by solution of the inverse problems of stochastic modeling of the accelerograms of the earthquakes, recorded by the seismic station in Sochi. The horizontal components of the strongest earthquakes (M _w ～ 3.9?5.6), that occurred in 2002–2006 within a radius of ～300 km from the seismic station, with source depths up to 60 km are modeled. For calculations of accelerograms, estimates of the quality are used, obtained earlier for this region in the form: Q(f) ～ 80 ～ f ^0.9. The parameter settings are carried out, which determine the shapes of the source spectra, the amplification of the seismic waves in the Earth’s crust, the weakening of the waves at high frequencies (κ), the parameters that determine the shape and duration of accelerograms, etc. Sufficiently good agreement of the calculated and recorded accelerograms is obtained, the regional characteristics of emission and propagation of seismic waves, which can be used for prediction of the parameters of strong motions in the North Caucasus, are evaluated; however, in the future these characteristics should be studied in more detail.     

Out-of-plane shaking table tests of full-scale historic adobe corner walls retrofitted with timber elements

Historic adobe structures pose a high seismic risk mainly because of the poor out-of-plane bending response of their walls that may produce fatalities and significant economic, cultural, and heritage losses. In this paper, we propose a retrofitting technique that increases the wall strength for both in-plane and out-of-plane directions. This technique consists of vertical and horizontal timber elements symmetrically installed on each face of the wall to form a confining wood frame, supplemented with vertical tensors that pre-compress the wall. This study evaluates the performance of this retrofitting technique with a two-set experimental program on full-scale historic adobe walls. On the first set, four specimens were subjected to a static overturning test with boundary conditions representing the confinement effect at both ends by orthogonal walls. On the second set, three full-scale specimens, one unretrofitted and two retrofitted, were subjected to four ground motion records on a shaking table to assess the out-of-plane dynamic behavior of typical corner walls. The unretrofitted specimen collapsed during the second motion (peak ground acceleration [PGA] = 0.39 g), while both retrofitted walls survived all four motions (maximum PGA of 0.75 g) proving the high effectiveness of the proposed retrofitting. The addition of base anchors as a variation of the retrofitting technique significantly reduced the rocking effects and the residual drifts of the system, thus improving its overall seismic performance. Further research is needed to develop guidelines for seismic retrofit of heritage buildings including multistory full-scale tests of specimens with various types of openings and retrofitting strategies that minimize their architectural impact.     

Historical intermediate-depth earthquakes in the southern Aegean Sea Benioff zone: modeling their anomalous macroseismic patterns with stochastic ground-motion simulations

We model the macroseismic damage distribution of four important intermediate-depth earthquakes of the southern Aegean Sea subduction zone, namely the destructive 1926 M?=?7.7 Rhodes and 1935 M?=?6.9 Crete earthquakes, the unique 1956 M?=?6.9 Amorgos aftershock (recently proposed to be triggered by a shallow event), and the more recent 2002 M?=?5.9 Milos earthquake, which all exhibit spatially anomalous macroseismic patterns. Macroseismic data for these events are collected from published macroseismic databases and compared with the spatial distribution of seismic motions obtained from stochastic simulation, converted to macroseismic intensity (Modified Mercalli scale, I_MM). For this conversion, we present an updated correlation between macroseismic intensities and peak measures of seismic motions (PGA and PGV) for the intermediate-depth earthquakes of the southern Aegean Sea. Input model parameters for the simulations, such as fault dimensions, stress parameters, and attenuation parameters (e.g. back-arc/along anelastic attenuation) are adopted from previous work performed in the area. Site-effects on the observed seismic motions are approximated using generic transfer functions proposed for the broader Aegean Sea area on the basis of V_S30 values from topographic slope proxies. The results are in very good agreement with the observed anomalous damage patterns, for which the largest intensities are often observed at distances >?100 km from the earthquake epicenters. We also consider two additional “prediction” but realistic intermediate-depth earthquake scenarios, and model their macroseismic distributions, to assess their expected damage impact in the broader southern Aegean area. The results suggest that intermediate-depth events, especially north of central Crete, have a prominent effect on a wide area of the outer Hellenic arc, with a very important impact on modern urban centers along northern Crete coasts (e.g. city of Heraklion), in excellent agreement with the available historical information.