Damage data analysis and vulnerability estimation following the August 14, 2003 Lefkada Island, Greece, Earthquake

Earthquake damage data is of vital importance both as an indicator of the likely performance of buildings in future earthquakes, as well as for the calibration of existing theoretical-analytical models regarding building vulnerability. The analysis of damage data collected shortly after the August 14, 2003 Lefkada Island, Greece earthquake revealed the higher vulnerability of masonry buildings vis-à-vis all other building typologies on the island. However, promising means in strengthening the existing masonry stock emerge, when considering the improved performance of buildings of a dual-system of stone-masonry and timber frame – a construction practice uniquely adopted in the island. Based on the parameterless seismic intensity scale (psi) proposed by Spence et al. (1991) a set of preliminary vulnerability curves for the typologies of the island’s buildings are proposed.     

Implementation of ADI-FDTD subgrids in ground penetrating radar FDTD models

Realistic numerical modeling of ground penetrating radar (GPR) using the finite-difference time-domain (FDTD) method could greatly benefit from the implementation of subgrids – supporting finer spatial resolution – into the conventional FDTD mesh. This is particularly important, when parts of the computational domain need to be modeled in detail or when there are features or regions in the overall computational mesh with values of high dielectric constant supporting propagation of waves at very short wavelengths. A scheme that simplifies the process of implementing these subgrids into the traditional FDTD method is presented. This scheme is based on the combination of the standard FDTD method and the unconditionally stable alternating-direction implicit (ADI) FDTD technique. Because ADI-FDTD is unconditionally stable its time-step can be set to any value that facilitates the accurate calculation of the fields. By doing so, the two grids can efficiently communicate information across their boundary without requiring to use a costly time-interpolation scheme. This paper discusses the performance of ADI-FDTD subgrids when implemented into the traditional FDTD method, using different communication schemes for the information exchange at the boundary of the two grids. The developed algorithm, can handle cases where the subgrid crosses dielectrically inhomogeneous media. In addition, results from the comparison between the proposed scheme and a commonly employed purely FDTD subgridding technique are presented.     

Using detrital garnets to determine provenance: a case study from the Vertiskos Unit (Serbomacedonian Massif,N. Greece)

Garnet single crystals of several millimeters in diameter were collected from the uppermost horizon of a soil profile developing immediately on the gneissic rocks of the Vertiskos Unit of the Serbomacedonian Massif in northern Greece. The garnets were analyzed for major elements by EDS analyzer mounted on a scanning electron microscope, and the obtained data were utilized to determine their source rocks. Bivariate diagrams, spider diagrams as well as statistical analysis were used in order to correlate and compare the garnet composition of the basement rocks of the Vertiskos Unit with the existing reference data. This case study demonstrates the difficulty in assigning a source rock to sediment, using only the chemical compositional of detrital garnet. Direct linking of the detrital garnets and the outcropping rocks is not always possible despite well documented outcrop lithologies. This is largely due to a complex metamorphic evolution that leads to overlapping compositions between garnets originating from different lithologies that have undergone similar metamorphic processes and alteration effects.     

Natural radioactivity distribution and gamma radiation exposure of beach sands from Sithonia Peninsula

This study aims to evaluate the activity concentrations of ²³⁸U, ²²⁶Ra, ²³²Th, ²²⁸Th and ⁴⁰K along the beaches of Sithonia Peninsula which are adjacent to the rock-types of the Sithonia Plutonic Complex. These range from 6–673, 5–767, 5–1750, 6–1760 and 185–875 Bq/kg respectively. The (% wt.) heavy magnetic (HM) fraction (epidote, allanite, hornblende, biotite and garnet), the heavy non-magnetic (HNM) fraction (monazite, zircon, titanite and apatite) and the total heavy (TH) fraction, were correlated with the concentrations of the measured radionuclides in the bulk samples. The HNM fraction seems to control the activity concentrations of ²³⁸U in all samples, while the HM fraction, at least for the heavy mineral rich samples bearing high amounts of epidote crystals with allanite cores, controls their ²³²Th content. The measured radionuclides in beach sands were normalized to the respective values measured in the granitic rocks, which are their most probable parent rocks, in order to provide data on their enrichment or depletion. The annual effective dose varies from 0.013 to 0.688 mSv y^?1 for local people working on the beach, while for tourists the annual external effective dose ranges between 0.003 and 0.165 mSv y^?1.     

Influence of deterioration modelling on the seismic response of steel moment frames designed to Eurocode 8

This paper assesses the influence of cyclic and in‐cycle degradation on seismic drift demands in moment‐resisting steel frames (MRF) designed to Eurocode 8. The structural characteristics, ground motion frequency content, and level of inelasticity are the primary parameters considered. A set of single‐degree‐of‐freedom (SDOF) systems, subjected to varying levels of inelastic demands, is initially investigated followed by an extensive study on multi‐storey frames. The latter comprises a large number of incremental dynamic analyses (IDA) on 12 frames modelled with or without consideration of degradation effects. A suite of 56 far‐field ground motion records, appropriately scaled to simulate 4 levels of inelastic demand, is employed for the IDA. Characteristic results from a detailed parametric investigation show that maximum response in terms of global and inter‐storey drifts is notably affected by degradation phenomena, in addition to the earthquake frequency content and the scaled inelastic demands. Consistently, both SDOF and frame systems with fundamental periods shorter than the mean period of ground motion can experience higher lateral strength demands and seismic drifts than those of non‐degrading counterparts in the same period range. Also, degrading multi‐storey frames can exhibit distinctly different plastic mechanisms with concentration of drifts at lower levels. Importantly, degrading systems might reach a “near‐collapse” limit state at ductility demand levels comparable to or lower than the assumed design behaviour factor, a result with direct consequences on optimised design situations where over‐strength would be minimal. Finally, the implications of the findings with respect to design‐level limit states are discussed.     

Deformation patterns in the southwestern part of the Mediterranean Ridge (South Matapan Trench,Western Greece)

Seismic reflection data and bathymetry analyses, together with geological information, are combined in the present work to identify seabed structural deformation and crustal structure in the Western Mediterranean Ridge (the backstop and the South Matapan Trench). As a first step, we apply bathymetric data and state of art methods of pattern recognition to automatically detect seabed lineaments, which are possibly related to the presence of tectonic structures (faults). The resulting pattern is tied to seismic reflection data, further assisting in the construction of a stratigraphic and structural model for this part of the Mediterranean Ridge. Structural elements and stratigraphic units in the final model are estimated based on: (a) the detected lineaments on the seabed, (b) the distribution of the interval velocities and the presence of velocity inversions, (c) the continuity and the amplitudes of the seismic reflections, the seismic structure of the units and (d) well and stratigraphic data as well as the main tectonic structures from the nearest onshore areas. Seabed morphology in the study area is probably related with the past and recent tectonics movements that result from African and European plates’ convergence. Backthrusts and reverse faults, flower structures and deep normal faults are among the most important extensional/compressional structures interpreted in the study area.     

Sensitivity of tidal characteristics in double inlet systems to momentum dissipation on tidal flats: a perturbation analysis

In a tidal channel with adjacent tidal flats, along–channel momentum is dissipated on the flats during rising tides. This leads to a sink of along–channel momentum. Using a perturbative method, it is shown that the momentum sink slightly reduces the M₂ amplitude of both the sea surface elevation and current velocity and favours flood dominant tides. These changes in tidal characteristics (phase and amplitude of sea surface elevations and currents) are noticeable if widths of tidal flats are at least of the same order as the channel width, and amplitudes and gradients of along–channel velocity are large. The M₂ amplitudes are reduced because stagnant water flows from the flats into the channel, thereby slowing down the current. The M₄ amplitudes and phases change because the momentum sink acts as an advective term during the fall of the tide, such a term generates flood dominant currents. For a prototype embayment that resembles the Marsdiep–Vlie double–inlet system of the Western Wadden Sea, it is found that for both the sea surface elevation and current velocity, including the momentum sink, lead to a decrease of approximately 2% in M₂ amplitudes and an increase of approximately 25% in M₄ amplitudes. As a result, the net import of coarse sediment is increased by approximately 35%, while the transport of fine sediment is hardly influenced by the momentum sink. For the Marsdiep–Vlie system, the M₂ sea surface amplitude obtained from the idealised model is similar to that computed with a realistic three–dimensional numerical model whilst the comparison with regard to M₄ improves if momentum sink is accounted for.     

Estimation of seawater retreat timescales in homogeneous and confined coastal aquifers based on dimensional analysis

Reliable estimation of the timescales of seawater retreat (SWR) in coastal aquifers has significant implications for sustainable groundwater management in coastal areas. Due to the complexity of coastal dynamics, analytical estimates of SWR timescales are limited in number, and numerical models are mainly utilized. Although numerical models allow for accurate analysis of the coastal systems, they concern particular aquifer cases, thus their results are not generally applicable. Here, we perform a dimensional analysis for the depth-integrated freshwater/seawater continuity equations and identify the dimensionless parameters that affect SWR in confined, homogeneous and flux-controlled aquifers. Based on numerical simulations, we gain insight into how these parameters affect SWR, and produce type curves which can be used as a simple tool for estimating SWR timescales in idealized coastal systems, independent of the hydraulic and geometric characteristics of the aquifer. The reliability of our estimates in real-world applications is also discussed.     

A hybrid optimization scheme for self-potential measurements due to multiple sheet-like bodies in arbitrary 2D resistivity distributions

Self-potential is a passive geophysical method that can be applied in a straightforward manner with minimum requirements in the field. Nonetheless, interpretation of self-potential data is particularly challenging due to the inherited non-uniqueness present in all potential methods. Incorporating information regarding the target of interest can facilitate interpretation and increase the reliability of the final output. In the current paper, a novel method for detecting multiple sheet-like targets is presented. A numerical framework is initially described that simulates sheet-like bodies in an arbitrary 2D resistivity distribution. A scattered field formulation based on finite differences is employed that allows the edges of the sheet to be independent of the grid geometry. A novel analytical solution for two-layered models is derived and subsequently used to validate the accuracy of the proposed numerical scheme. Lastly, a hybrid optimization is proposed that couples linear least-squares with particle-swarm optimization in order to effectively locate the edges of multiple sheet-like bodies. Through numerical and real data, it is proven that the hybrid optimization overcomes local minimal that occurs in complex resistivity distributions and converges substantially faster compared to traditional particle-swarm optimization.