Petrological features of boseti mountains,a complex volcanic system in the axial portion of the main Ethiopian rift

A petrographical and chemical research carried out mainly for the rocks belonging to the 2nd cycle of volcanic activity of Boseti Mts., shows the predominance of the salic over the mafic and the intermediate products. The magmatic evolution, which chiefly took place in a thermodynamic environment with initial low \(P_{H_2 O} \) and \(f_{O_2 } \) , has brought about subaluminiferous mugearitic types, also with trachytic affinity. Among peralkali-rhyolites, an anorthoclase fractionation can explain both their chemical variation and the production of pantelleritic liquids.     

Neogene calc-alkaline extrusive and intrusive rocks of the karalar-yesiller area

The Karalar-Ye?iller area lies on the estern flank of the Kazda \(\bar g\) i massif in northwestern Anatolia (Asiatic Turkey), and includes about 600 square kilometers, of which approximately 80% is covered by rhyodacite-quartz latite extrusive rocks, and a comagmatic granodiorite-quartz monzonite batholith, and stocks, all of middle Miocene age. Extrusive rocks consist chiefly of rhyodacite lava flows of the Hallaçlar Formation and quartz latite-rhyodacite ash-flow tuff, lava, and mudflow deposits of the overlying Dede Tepe Formation. These volcanic rocks lie on a basement composed, in asceding stratigraphic order, of: 1) pre-Permian Kalabak sequence, 2) Upper Triassic Halilar Formation, 3) post-Upper Triassic Ba \(\bar g\) burun Formation, and 4) allochthons of middle Permian: and Upper Jurassic limestone. Intrusion of a granodiorite-quartz monzonite batholith during the middle Miocene was accompanied or shortly followed by extensive alteration and base metal mineralization of the Hallaçlar and older rocks. Intrusion of six rhyodacitequartz latite stocks followed the main phase of hydrothermal alteration. At least one of these stocks may have contributed material to the Dede Tepe Formation of middle Miocene age. Field relations, petrologic, and geochemical data as well as radiometric age dates suggest that the intrusive and extrusive rocks are comagmatic. Parent magma is probably derived from partial melting of subducted oceanic crust and accompanying oceanic sediments. Geologic relations locally indicate assimilation of sialic crust by the grandiorite-quartz monzonite batholith.     

Quasi-Static Strain and Stress Fields due to a Moment Tensor in Elastic–Viscoelastic Layered Half-Space

We derived explicit expressions in the time domain for 3-D quasi-static strain and stress fields, due to a point moment tensor source in an elastic surface layer overlying viscoelastic half-space under gravity. The expressions of strain in the elastic surface layer were directly obtained from the expressions of displacement in our previous paper. The conversion of strain into stress is easy, because the stress–strain relation of elastic material is linear. In the viscoelastic substratum, the expressions of strain were obtained by applying the correspondence principle of linear viscoelasticity to the associated elastic solution. The strain–stress conversion is not straightforward, as the stress–strain relation of viscoelastic material is usually given in a differential form. To convert strain into stress, we used an integral form of the stress–strain relation instead of the usual differential form. The expressions give the responses of elastic half-space at \( t = 0 \) , and the responses of an elastic plate floating on non-viscous liquid at \( t = \infty \) . The moment tensor is rationally decomposed into the three independent force systems, corresponding to isotropic expansion, shear faulting and crack opening, and so the expressions include the strain and stress fields for these force systems as special cases. As the first numerical example, we computed the temporal changes in strain and stress fields after the sudden opening of an infinitely long vertical crack cutting the elastic surface layer. Here, we observe that the stress changes caused by the sudden crack opening gradually decay with time and vanish at \( t = \infty \) everywhere. After the completion of stress relaxation, a characteristic pattern of shear strain remains in the viscoelastic substratum. Since the strain and stress fields at \( t = \infty \) can be read as the strain- and stress-rate fields caused by steady crack opening, respectively, this numerical example demonstrates the realization of a steady stress state supported by steady viscous flow in the asthenosphere, associated with steady seafloor spreading at mid-ocean ridges. For the second numerical example, we computed the temporal changes in strain and stress fields after the 2011 Tohoku-oki mega-thrust earthquake, which occurred at the North American-Pacific plate interface. In this numerical example, the stress changes caused by coseismic fault slip vanish at \( t = \infty \) in the viscoelastic substratum, but remain in the elastic surface layer. The coseismic stress changes (and also strain changes) in the elastic surface layer diffuse away from the source region with time, due to gradual stress relaxation in the viscoelastic substratum.     

Complexity–entropy analysis of daily stream flow time series in the continental United States

Complexity–entropy causality plane (CECP) is a diagnostic diagram plotting normalized Shannon entropy \({\cal H}_S\) versus Jensen–Shannon complexity \({\cal C}_{JS}\) that has been introduced in nonlinear dynamics analysis to classify signals according to their degrees of randomness and complexity. In this study, we explore the applicability of CECP in hydrological studies by analyzing 80 daily stream flow time series recorded in the continental United States during a period of 75 years, surrogate sequences simulated by autoregressive models (with independent or long-range memory innovations), Theiler amplitude adjusted Fourier transform and Theiler phase randomization, and a set of signals drawn from nonlinear dynamic systems. The effect of seasonality, and the relationships between the CECP quantifiers and several physical and statistical properties of the observed time series are also studied. The results point out that: (1) the CECP can discriminate chaotic and stochastic signals in presence of moderate observational noise; (2) the signal classification depends on the sampling frequency and aggregation time scales; (3) both chaotic and stochastic systems can be compatible with the daily stream flow dynamics, when the focus is on the information content, thus setting these results in the context of the debate on observational equivalence; (4) the empirical relationships between \({\mathcal H}_S\) and \({\mathcal C}_{JS}\) and Hurst parameter H, base flow index, basin drainage area and stream flow quantiles highlight that the CECP quantifiers can be considered as proxies of the long-term low-frequency groundwater processes rather than proxies of the short-term high-frequency surface processes; (6) the joint application of linear and nonlinear diagnostics allows for a more comprehensive characterization of the stream flow time series.     

Mantle Transition Zone Structure Beneath Southeastern China and its Implications for Stagnant Slab and Water Transportation in the Mantle

We determined depth variation of the 410- and 660-km discontinuities beneath southeastern China by common-converted-point stacking of \(\rm P\) -wave receiver functions of 121 permanent Chinese seismic stations. We then combined the results with seismic velocity variation to estimate temperature and water content variations in the mantle transition zone of the region. Previous tomographic studies have shown a stagnant slab in the mantle transition zone in eastern Asia that is connected to subduction of the western Pacific. Temperature variations obtained clearly outline the shape of the stagnant slab, with its western edge at 113.5 \(^\circ\) E and the southern edge at 28.5 \(^\circ\) N. The correlation between the location of the stagnant slab and surface tectonics suggests that the Cenozoic extension in eastern China is closely associated with the subduction of the western Pacific and its eastward migration. The water content of the stagnant slab is lower than in surrounding slabs, suggesting that the water has already been released from the subducting slab into the upper mantle.     

Interannual variability of the air–sea CO2 flux in the north Indian Ocean

A simple biogeochemical model coupled to an offline ocean tracer transport model driven by reanalysis ocean data is used to simulate the seasonal and interannual CO $_2$ flux variability in the northern Indian Ocean. The maximum of seasonal and interannual CO $_2$ emission variances in the northern Indian Ocean are located in the coastal Arabian Sea (AS) and Southern Peninsular India (SP) with a basin-wide seasonal amplitude and standard deviation of 0.044 $\pm $ 0.04 Pg C year $^{-1}$ . The area integrated CO $_2$ emissions from these two regions in the model are significantly correlated (above a 95 % level) with the observations of Takahashi et al. (Deep-Sea Res-II, 56:554–577, 2009). The interannual anomalies of CO $_2$ emission from the AS and SP are found as 40 and 30 % of their respective seasonal amplitudes. Both the Arabian Sea (AS) and Southern Peninsular India (SP) interannual CO $_2$ emission anomalies show a 3–4-year variability. The correlations of AS and SP CO $_2$ emission anomalies with the Indian Ocean Dipole/Zonal Mode (IODZM) and Southern Oscillation (SO) indices from 1980 to 1999 are 0.35, 0.21 and 0.32, 0.01 respectively. A 5-year window moving correlation analysis shows that the relationship of AS and SP CO $_2$ emission to the SO and IODZM are complementary to each other. During the years when the correlation of air–sea CO $_2$ emission with the IODZM is stronger, the corresponding correlation with the SO is weaker or opposite. The total change in pCO $_2$ is broken down into changes induced by the individual components such as dissolved inorganic carbon (DIC), sea surface temperature (SST), alkalinity, and salinity and found that (1) the effect of SST in the AS CO $_2$ emission increases (decreases) when the correlation of CO $_2$ emission with the IODZM is positive (negative), and (2) the SP CO $_2$ emission is strongly controlled by the circulation-driven DIC changes; however, this relation is found to be weaker when the SO correlates negatively with the SP CO $_2$ emission.     

Slope instabilities triggered by the 11th May 2011 Lorca earthquake (Murcia,Spain): comparison to previous hazard assessments and proposition of a new hazard map and probability of failure equation

The Lorca Basin has been the object of recent research aimed at studying the phenomena of earthquake-induced landslides and its assessment in the frame of different seismic scenarios. However, it has not been until the 11th May 2011 Lorca earthquakes when it has been possible to conduct a systematic approach to the problem. In this paper we present an inventory of slope instabilities triggered by the Lorca earthquakes which comprises more than 100 cases, mainly rock and soil falls of small size (1–100  \(\hbox {m}^{3}\) ). The distribution of these instabilities is here compared to two different earthquake-triggered landslide hazard maps: one considering the occurrence of the most probable earthquake for a 475-years return period in the Lorca Basin \((\hbox {M}_{\mathrm{w}}=5.0)\) based on both low- and high-resolution digital elevation model (DEM); and a second one matching the occurrence of the \(\hbox {M}_{\mathrm{w}}=5.2\) 2011 Lorca earthquake, which was performed using the higher resolution DEM. The most frequent Newmark displacements related to the slope failures triggered by the 2011 Lorca earthquakes are lower than 2 cm in both the hazard scenarios considered. Additionally, the predicted Newmark displacements were correlated to the inventory of slope instabilities to develop a probability of failure equation. The fit seems to be very good since most of the mapped slope failures are located on the higher probability areas. The probability of slope failure in the Lorca Basin for a seismic event similar to the \(\hbox {M}_{\mathrm{w}}\) 5.2 2011 Lorca earthquake can be considered as very low (0–4 %).     

Sub-Gaussian model of processes with heavy-tailed distributions applied to air permeabilities of fractured tuff

Earth and environmental variables are commonly taken to have multivariate Gaussian or heavy-tailed distributions in space and/or time. This is based on the observation that univariate frequency distributions of corresponding samples appear to be Gaussian or heavy-tailed. Of particular interest to us is the well-documented but heretofore little noticed and unexplained phenomenon that whereas the frequency distribution of log permeability data often seems to be Gaussian, that of corresponding increments tends to exhibit heavy tails. The tails decay as powers of ? $ \alpha $ where 1 <  $ \alpha $  < 2 is either constant or grows monotonically toward an asymptote with increasing separation distance or lag. We illustrate the latter phenomenon on 1-m scale log air permeabilities from pneumatic tests in 6 vertical and inclined boreholes completed in unsaturated fractured tuff near Superior, Arizona. We then show theoretically and demonstrate numerically, on synthetically generated signals, that whereas the case of constant $ \alpha $ is consistent with a collection of samples from truncated sub-Gaussian fractional Lévy noise, a random field (or process) subordinated to truncated fractional Gaussian noise, the case of variable $ \alpha $ is consistent with a collection of samples from truncated sub-Gaussian fractional Lévy motion (tfLm), a random field subordinated to truncated fractional Brownian motion. Whereas the first type of signal is relatively regular and characterized by Lévy index $ \alpha $ , the second is highly irregular (punctuated by spurious spikes) and characterized by the asymptote of $ \alpha $ values associated with its increments. We describe a procedure to estimate the parameters of univariate distributions characterizing such signals and apply it to our log air permeability data. The latter are found to be consistent with a collection of samples from tfLm with $ \alpha $ slightly smaller than 2, which is easily confused with a Gaussian field (characterized by constant $ \alpha $  = 2). The irregular (spiky) nature of this signal is typical of observed fractured rock properties. We propose that distributions of earth and environmental variable be inferred jointly from measured values and their increments in a way that insures consistency between these two sets of data.     

Empirical ground-motion relations using moderate earthquakes recorded by Medellín–Aburrá Valley (Colombia) strong-motion networks

We tested attenuation relations obtained for different regions of the world to verify their suitability to predict strong-motion data recorded by Medellín and Aburrá Valley Accelerographic Networks. We used as comparison criteria, the average of the difference between the observed and the predicted data as a function of epicenter distance and its standard deviation. We also used the approach developed by Sherbaum et al. (Bull Seism Soc Am 94:2164–2185, 2004) that provides a method to evaluate the overall goodness-of-fit of ground-motion prediction equations. The predictive models selected use a generic focal depth. We found that this parameter has an important influence in the ground-motion predictions and must be taken into account as an independent variable. We also found important to characterize the local soil amplification to improve the attenuation relations. We found empirical relations for peak horizontal acceleration PGA and velocity PGV based on the Kamiyama and Yanagisawa (Soils Found 26:16–32, 1986) approach. $$\begin{aligned} \log _{10} (PGA)=0.5886M_L -1.0902\log _{10}(R)-0.0035H+C_{st}\pm 0.\text{29} \end{aligned}$$ $$\begin{aligned} \log _{10} (PGV)=0.7255M_L -1.8812\log _{10}(R)-0.0016H+C_{st}\pm 0.36 \end{aligned}$$ where PGA is measured in cm/s $^{2}$ and PGV in cm/s, $M_{L}$ is local magnitude in the range 2.8–6.5, $R$ is epicentral distance up to 290 km, $H$ is focal depth in km and $C_{st}$ is a coefficient that accounts for the site response due to soil conditions of each recording station. The introduction of focal depth and local site conditions as independent variables, minimize the residuals and the dispersion of the predicted data. We conclude that $H$ and $C_{st}$ are sensitive parameters, having a strong influence on the strong-motion predictions. Using the same functional form, we also propose an empirical relation for the root mean square acceleration a $_\mathrm{rms}$ : $$\begin{aligned} \log _{10} \left( {a_{rms} } \right)=0.4797M_L -1.1665\log _{10} (R)-0.00201H+C_{st}\pm 0.40 \end{aligned}$$ where a $_\mathrm{rms}$ is measured in cm/s $^{2}$ , from the S-wave arrival and using a window length equal to the rupture duration. The other variables are the same as those for PGA and PGV. The site correction coefficients $C_{st}$ found for PGA, PGV and a $_\mathrm{rms}$ show a similar trend indicating a good correlation with the soil conditions of the recording sites.     

On the ability of Moho reflections to affect the ground motion in northeastern Italy: a case study of the 2012 Emilia seismic sequence

It has been observed that post-critically reflected S-waves and multiples from the Moho discontinuity could play a relevant role on the ground motion due to medium to strong size earthquakes away from the source. Although some studies investigated the correlation between the Moho reflections amplitudes and the damage in the far field, little attention was given to the frequency content of these specific phases and their scaling with magnitude. The 2012 Emilia seismic sequence in northern Italy, recorded by velocimetric and accelerometric networks, is here exploited to investigate Moho reflections and multiples (SmSM). A single station method for group velocity-period estimation, based on the multiple filter technique, is applied to strong motion data to detect SmSM. Amplitude and frequency scaling with magnitude is defined for earthquakes from \(\hbox {Mw}=3.9\) to \(\hbox {Mw}=5.9\) . Finally, the ability of SmSM to affect the ground motion for a maximum credible earthquake within the Po plain is investigated by extrapolating observed engineering parameters. Data analysis shows that high amplitude SmSM can be recognized within the Po plain, and at the boundaries between the Po plain and the Alpine chain, at epicentral distances larger than 80 km, in the period range from 0.25 to 3 s and in the group velocity window from about 2.6 to 3.2 km/s. 5 % damped pseudo-spectral accelerations at different periods (0.3, 1.0 and 2.0 s), and Housner intensities, are obtained from data characterized by large amplitude SmSM. A scaling relationship for both pseudo-spectral accelerations and Housner intensities is found for the earthquakes of the 2012 Emilia seismic sequence. \(\hbox {I}_{\mathrm{MCS}}\) from VII to VIII is estimated, as a result of SmSM amplitude enhancement, at about 100 km for a maximum credible earthquake ( \(\hbox {Mw}=6.7\) ) in the Po plain, showing that moderate to high damage cloud be caused by these specific phases.     

The role of site effects in the comparison between code provisions and the near field strong motion of the Emilia 2012 earthquakes

A damaging seismic sequence hit a wide area mainly located in the Emilia-Romagna region (Northern Italy) during 2012 with several events of local magnitude \(\hbox {M}_\mathrm{l} \ge 5\) , among which the \(\hbox {M}_\mathrm{l}\) 5.9 May 20 and the \(\hbox {M}_\mathrm{l}\) 5.8 May 29 were the main events. Thanks to the presence of a permanent accelerometric station very close to the epicentre and to the temporary installations performed in the aftermath of the first shock, a large number of strong motion recordings are available, on the basis of which, we compared the recorded signals with the values provided by the current Italian seismic regulations, and we observed several differences with respect to horizontal components when the simplified approach for site conditions (based on Vs30 classes) is used. On the contrary, when using the more accurate approach based on the local seismic response, we generally obtain a much better agreement, at least in the frequency range corresponding to a quarter wavelength comparable with the depth of the available subsoil data. Some unresolved questions still remain, such as the low frequency behaviour ( \(<\) 1 Hz) that could be due either to complex propagation at depth larger than the one presently investigated or to near source effects, and the behaviour of vertical spectra whose recorded/code difference is too large to be explained with the information currently available.     

Scaling Transition in Earthquake Sources: A Possible Link Between Seismic and Laboratory Measurements

We estimate the corner frequencies of 20 crustal seismic events from mainshock–aftershock sequences in different tectonic environments (mainshocks 5.7 < M _W < 7.6) using the well-established seismic coda ratio technique (Mayeda et al. in Geophys Res Lett 34:L11303, 2007; Mayeda and Malagnini in Geophys Res Lett, 2010), which provides optimal stability and does not require path or site corrections. For each sequence, we assumed the Brune source model and estimated all the events’ corner frequencies and associated apparent stresses following the MDAC spectral formulation of Walter and Taylor (A revised magnitude and distance amplitude correction (MDAC2) procedure for regional seismic discriminants, 2001), which allows for the possibility of non-self-similar source scaling. Within each sequence, we observe a systematic deviation from the self-similar \( M_{0} \propto \mathop f\nolimits_{\text{c}}^{ - 3} \) line, all data being rather compatible with \( M_{0} \propto \mathop f\nolimits_{\text{c}}^{ - (3 + \varepsilon )} \) , where ε > 0 (Kanamori and Rivera in Bull Seismol Soc Am 94:314–319, 2004). The deviation from a strict self-similar behavior within each earthquake sequence of our collection is indicated by a systematic increase in the estimated average static stress drop and apparent stress with increasing seismic moment (moment magnitude). Our favored physical interpretation for the increased apparent stress with earthquake size is a progressive frictional weakening for increasing seismic slip, in agreement with recent results obtained in laboratory experiments performed on state-of-the-art apparatuses at slip rates of the order of 1 m/s or larger. At smaller magnitudes (M _W < 5.5), the overall data set is characterized by a variability in apparent stress of almost three orders of magnitude, mostly from the scatter observed in strike-slip sequences. Larger events (M _W > 5.5) show much less variability: about one order of magnitude. It appears that the apparent stress (and static stress drop) does not grow indefinitely at larger magnitudes: for example, in the case of the Chi–Chi sequence (the best sampled sequence between M _W 5 and 6.5), some roughly constant stress parameters characterize earthquakes larger than M _W ~ 5.5. A representative fault slip for M _W 5.5 is a few tens of centimeters (e.g., Ide and Takeo in J Geophys Res 102:27379–27391, 1997), which corresponds to the slip amount at which effective lubrication is observed, according to recent laboratory friction experiments performed at seismic slip velocities (V ~ 1 m/s) and normal stresses representative of crustal depths (Di Toro et al. in Nature in press, 2011, and references therein). If the observed deviation from self-similar scaling is explained in terms of an asymptotic increase in apparent stress (Malagnini et al. in Pure Appl Geophys, 2014, this volume), which is directly related to dynamic stress drop on the fault, one interpretation is that for a seismic slip of a few tens of centimeters (M _W ~ 5.5) or larger, a fully lubricated frictional state may be asymptotically approached.     

Discretization Errors in the Hybrid Finite Element Particle-in-cell Method

In computational geodynamics, the Finite Element (FE) method is frequently used. The method is attractive as it easily allows employment of body-fitted deformable meshes and a true free surface boundary condition. However, when a Lagrangian mesh is used, remeshing becomes necessary at large strains to avoid numerical inaccuracies (or even wrong results) due to severely distorted elements. For this reason, the FE method is oftentimes combined with the particle-in-cell (PIC) method, where particles are introduced which track history variables and store constitutive information. This implies that the respective material properties have to be interpolated from the particles to the integration points of the finite elements. In numerical geodynamics, material parameters (in particular the viscosity) usually vary over a large range. This may be due to strongly temperature-dependent rheologies (which result in large but smooth viscosity variations) or material interfaces (which result in viscosity jumps). Here, we analyze the accuracy and convergence properties of velocity and pressure of the hybrid FE-PIC method in the presence of large viscosity variations. Standard interpolation schemes (arithmetic and harmonic) are compared to a more sophisticated interpolation scheme which is based on linear least squares interpolation for two types of elements ( \(Q_1P_0\) and \(Q_2P_{-1}\) ). In the case of a smooth viscosity field, the accuracy and convergence is significantly improved by the new interpolation scheme. In the presence of viscosity jumps, the order of accuracy is strongly decreased.     

Characterization of the Tail of the Distribution of Earthquake Magnitudes by Combining the GEV and GPD Descriptions of Extreme Value Theory

The present work is a continuation and improvement of the method suggested in Pisarenko et al. (Pure Appl Geophys 165:1–42, 2008) for the statistical estimation of the tail of the distribution of earthquake sizes. The chief innovation is to combine the two main limit theorems of Extreme Value Theory (EVT) that allow us to derive the distribution of T-maxima (maximum magnitude occurring in sequential time intervals of duration T) for arbitrary T. This distribution enables one to derive any desired statistical characteristic of the future T-maximum. We propose a method for the estimation of the unknown parameters involved in the two limit theorems corresponding to the Generalized Extreme Value distribution (GEV) and to the Generalized Pareto Distribution (GPD). We establish the direct relations between the parameters of these distributions, which permit to evaluate the distribution of the T-maxima for arbitrary T. The duality between the GEV and GPD provides a new way to check the consistency of the estimation of the tail characteristics of the distribution of earthquake magnitudes for earthquake occurring over an arbitrary time interval. We develop several procedures and check points to decrease the scatter of the estimates and to verify their consistency. We test our full procedure on the global Harvard catalog (1977–2006) and on the Fennoscandia catalog (1900–2005). For the global catalog, we obtain the following estimates: \( \hat{M}_{{\rm max} } \)  = 9.53 ± 0.52 and \( \hat{Q}_{10} (0.97) \)  = 9.21 ± 0.20. For Fennoscandia, we obtain \( \hat{M}_{{\rm max} } \)  = 5.76 ± 0.165 and \( \hat{Q}_{10} (0.97) \)  = 5.44 ± 0.073. The estimates of all related parameters for the GEV and GPD, including the most important form parameter, are also provided. We demonstrate again the absence of robustness of the generally accepted parameter characterizing the tail of the magnitude-frequency law, the maximum possible magnitude M _max, and study the more stable parameter Q _T (q), defined as the q-quantile of the distribution of T-maxima on a future interval of duration T.     

Towards fully data driven ground-motion prediction models for Europe

We have used the Artificial Neural Network method (ANN) for the derivation of physically sound, easy-to-handle, predictive ground-motion models from a subset of the Reference database for Seismic ground-motion prediction in Europe (RESORCE). Only shallow earthquakes (depth smaller than 25 km) and recordings corresponding to stations with measured $V_{s30}$ properties have been selected. Five input parameters were selected: the moment magnitude $M_{W}$ , the Joyner–Boore distance $R_{JB}$ , the focal mechanism, the hypocentral depth, and the site proxy $V_{S30}$ . A feed-forward ANN type is used, with one 5-neuron hidden layer, and an output layer grouping all the considered ground motion parameters, i.e., peak ground acceleration (PGA), peak ground velocity (PGV) and 5 %-damped pseudo-spectral acceleration (PSA) at 62 periods from 0.01 to 4 s. A procedure similar to the random-effects approach was developed to provide between and within event standard deviations. The total standard deviation ( $\sigma $ ) varies between 0.298 and 0.378 (log $_{10}$ unit) depending on the period, with between-event and within-event variabilities in the range 0.149–0.190 and 0.258–0.327, respectively. Those values prove comparable to those of conventional GMPEs. Despite the absence of any a priori assumption on the functional dependence, our results exhibit a number of physically sound features: magnitude scaling of the distance dependency, near-fault saturation distance increasing with magnitude, amplification on soft soils and even indications for nonlinear effects in softer soils.     

Fragility curves for risk-targeted seismic design maps

Seismic design using maps based on “risk-targeting” would lead to an annual probability of attaining or exceeding a certain damage state that is uniform over an entire territory. These maps are based on convolving seismic hazard curves from a standard probabilistic analysis with the derivative of fragility curves expressing the chance for a code-designed structure to attain or exceed a certain damage state given a level of input motion, e.g. peak ground acceleration (PGA). There are few published fragility curves for structures respecting the Eurocodes (ECs, principally EC8 for seismic design) that can be used for the development of risk-targeted design maps for Europe. In this article a set of fragility curves for a regular three-storey reinforced-concrete building designed using EC2 and EC8 for medium ductility and increasing levels of design acceleration \((\hbox {a}_\mathrm{g})\) is developed. These curves show that structures designed using EC8 against PGAs up to about 1 m/s \(^{2}\) have similar fragilities to those that respect only EC2 (although this conclusion may not hold for irregular buildings, other geometries or materials). From these curves, the probability of yielding for a structure subjected to a PGA equal to \(\hbox {a}_\mathrm{g}\) varies between 0.14 ( \(\hbox {a}_\mathrm{g}=0.7\) m/s \(^{2})\) and 0.85 ( \(\hbox {a}_\mathrm{g}=3\) m/s \(^{2})\) whereas the probability of collapse for a structure subjected to a PGA equal to \(\hbox {a}_\mathrm{g}\) varies between 1.7 \(\times 10^{-7}\) ( \(\hbox {a}_\mathrm{g}=0.7\) m/s \(^{2})\) and 1.0 \(\times 10^{-5}\) ( \(\hbox {a}_\mathrm{g}=3\) m/s \(^{2})\) .     

Changes in dynamic characteristics of Lorca RC buildings from pre- and post-earthquake ambient vibration data

This paper estimates fundamental translational period and damping ratio parameters and examines the changes in dynamic characteristics of a set of low-to-medium rise buildings in Lorca town (SE of Spain) affected by the May 11th, 2011 earthquake. These building parameters have been calculated analysing structural dynamic response from ambient vibration measurements recorded at top RC buildings pre- and post earthquake, using the Fast Fourier Transform and the Randomdec technique. The empirical expression relating period \((T)\) and number of floor \((N)\) here obtained analysing ambient noise recorded on 59 healthy buildings before the earthquake is \(T= (0.054\pm 0.002)\, N\) , very similar to others empirical period–height relationships obtained for RC structures in the European built environment but quite different from code provisions. Measurements performed in 34 damaged buildings show a period elongation after the quake according to \(T^{*} =(0.075\pm 0.002)\,N\) expression. Moreover, we found a rise of the fundamental period with the EMS’s grade of damage of buildings. In contrast to natural frequency, damping ratio \((\xi )\) do not shows a significant variation with earthquake damage degree and the product \(\xi \, T\) remains near constant.     

An analytical method for the evaluation of the in-plan irregularity of non-regularly asymmetric buildings

This paper describes a new method for the evaluation of the static eccentricity $e_{s}$ and the ratio $\Omega _{\uptheta } $ of uncoupled torsional to lateral frequencies in real multi-storey buildings. The above-mentioned parameters greatly affect the lateral-to-torsional coupling of the response of asymmetric systems and thus are of paramount importance in the assessment of the in-plan irregularity of buildings. The proposed method, which is a generalization of that suggested by Calderoni et al. (Earthq Spectra 18(2):219–231, 2002), allows the calculation of the static eccentricity $e_{s}$ and the ratio $\Omega _{\uptheta } $ from the structural response to arbitrary distributions of forces and torsional couples. The effectiveness of the method is validated on some regularly and non-regularly asymmetric buildings characterised by different in-plan irregularity. The analyses demonstrate that the results of the method are rigorous in the case of regularly asymmetric systems and only slightly depend upon the heightwise distribution of the forces in the case of non-regularly asymmetric systems. Finally, the values of the static eccentricity $e_{s}$ and the ratio $\Omega _{\uptheta } $ resulting from the proposed method are compared to those obtained by means of the procedure suggested by Makarios and Anastassiadis in (Struct Des Tall Spec Build 7(1):33–55, 1998a; Struct Des Tall Spec Build 7(1):57–71, 1998b) .