An algorithm for fast 3D inversion of surface electrical resistivity tomography data: application on imaging buried antiquities

In this work a new algorithm for the fast and efficient 3D inversion of conventional 2D surface electrical resistivity tomography lines is presented. The proposed approach lies on the assumption that for every surface measurement there is a large number of 3D parameters with very small absolute Jacobian matrix values, which can be excluded in advance from the Jacobian matrix calculation, as they do not contribute significant information in the inversion procedure. A sensitivity analysis for both homogeneous and inhomogeneous earth models showed that each measurement has a specific region of influence, which can be limited to parameters in a critical rectangular prism volume. Application of the proposed algorithm accelerated almost three times the Jacobian (sensitivity) matrix calculation for the data sets tested in this work. Moreover, application of the least squares regression iterative inversion technique, resulted in a new 3D resistivity inversion algorithm more than 2.7 times faster and with computer memory requirements less than half compared to the original algorithm. The efficiency and accuracy of the algorithm was verified using synthetic models representing typical archaeological structures, as well as field data collected from two archaeological sites in Greece, employing different electrode configurations. The applicability of the presented approach is demonstrated for archaeological investigations and the basic idea of the proposed algorithm can be easily extended for the inversion of other geophysical data.     

Seismic fragility curves for greek bridges: methodology and case studies

This study focusses on the estimation of seismic fragility curves for all common bridge types found in modern greek motorways. At first a classification scheme is developed in order to classify the existing bridges into a sufficient number of classes. A total of 11 representative bridge classes resulted, based on the type of piers, deck, and pier-to-deck connection. Then an analytical methodology for deriving fragility curves is proposed and applied to the representative bridge models. This procedure is based on pushover analysis of the entire bridge and definition of damage states in terms of parameters of the bridge pushover curves. The procedure differentiates the way of defining damage according to the seismic energy dissipation mechanism in each bridge, i.e. bridges with yielding piers of the column type and bridges with bearings (with or without seismic links) and non-yielding piers of the wall type. The activation of the abutment-backfill system due to closure of the gap between the deck and the abutments is also taken into account. The derived fragility curves are subjected to a first calibration against empirical curves based on damage data from the US and Japan.     

The S(c) spectrum machine to visualize the motion in galaxies

The behavior of the orbits in a galaxy model composed of an harmonic core and a strong bar potential is studied. Numerical calculations show that a large number of orbits display chaotic motion. These orbits are low angular momentun orbits. The percentage of chaotic orbits increases as the angular velocity of the system increases or the strength of the harmonic term decreases. A new dynamical parameter, the S(c) spectrum, is introduced and used to detect the island motion and the evolution of the sticky regions. Comparison to previously obtained results reveals the leading role of the new spectrum. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Investigation of waste disposal areas using electrical methods: a case study from Chania,Crete, Greece

A geophysical investigation of a landfill area was conducted as inextricable phase of a preliminary geotechnical, geological and hydrogeological study of the area under investigation. Electrical resistivity tomography (ERT), one of the most promising prospecting techniques mainly concerning its effective contribution to resolve several environmental problems, was applied for the geophysical modeling. ERT is a robust imaging method the theory and implementation of which are well documented in geophysical research literature. The survey site is located in Akrotirion Peninsula, about 7 km east of Chania city in Crete I. The geological setting comprises Trypalion limestones, marly limestones and a clayed weathered layer. Cavities are also identified at various scales which are recent tectono-karstic voids. Due to the above-mentioned geological singularity and in the light of the requirement for an environmentally safe construction of the landfill, an ERT survey was carried out. Specifically, seven geoelectrical tomographies were conducted reaching the prospecting depth of 36 m. The resulted images conduced, to overcome geotechnical problems since the spatial distribution of karstic voids was determined, to plan the future sites for waste disposal as the geological conditions were studied in detail and to reliably estimate the thickness of the already deposited wastes. The resulted images were confirmed using available borehole logs and ambient noise measurements.     

Physical identification of static and axially-symmetric vacuum metric tensors

We present a method of a possible physical identification of the static and axially symmetric Weyl-type vacuum γ and nγ metrics. This method, in which no interior solutions of any kind are involved, is based on the comparison of the far-field forms of the γ and nγ metrics, and of the far-field form of the metric tensor due to a bounded gravitating perfect-fluid source given correctly to post-Newtonian accuracy. The parameters of the vacuum solutions are expressed in terms of physical parameters of a prolate fluid source, namely total mass-energy, semiaxes and eccentricity, defined consistently to post-Newtonian accuracy. The results, based on the otherwise arbitrary fluid source, appear physically general. Possible astrophysical candidates for the far-field γ andnγ metrics are proposed, based on the conditions imposed on them by the identification method. The advantages and deficiencies of the identification method are briefly discussed.     

Mainshock-consistent ground motion record selection for aftershock sequences

In recent years, the additional risk posed to the built environment due to aftershock sequences and triggered events has been brought to attention, and several efforts have been directed towards developing fragility functions for structures in damaged conditions. Despite this rise of interest, a rather fundamental component for such tasks, namely that of aftershock ground motion record selection, has remained under-scrutinized. Herein, we propose a pragmatic procedure that can be applied for the selection of mainshock-aftershock ground motion pairs using consistent causal parameters and accounting for the correlation between their spectral accelerations. In addition, a structural analysis strategy that can be employed for the analytical derivation of damage-dependent fragility functions is outlined and presented through a case study. A more conventional back-to-back IDA analysis is also carried out in order to compare the derived damage-dependent fragility functions with the ones obtained with the proposed procedure. The results indicate that record selection remains a crucial factor even when assessing the structural vulnerability of damaged buildings, and should thus be treated cautiously.     

Restructuring the kiosk trade in St. Petersburg: A new retail trade model for the post-Soviet period

This study explores the commercial kiosk phenomenon as a reflection of the socio-spatial reordering of St. Petersburg following the introduction of market forces. In this article we interpret empirical evidence on the restructuring of commercial opportunity in post-Soviet St. Petersburg as a result of the emergence of the kiosk phenomenon. Kiosk trade and processes that contribute to it have emerged as tools for mass socioeconomic and psychological adaptation, following the dismantling of the communist regime. In our opinion, the emergence and continuing mutation of St. Petersburg kiosks constitutes the backbone of a process of middle-class formation. An extensive field research of 1989, 1994 and 1995 made it possible to study the dynamics of spatial structures of kiosk phenomenon. It is our final conclusion that kiosks will remain an important feature of St. Petersburgs commercial geography in the foreseeable future.     

Evidence for periodic seismicity in the inner Aegean seismic zone

The most complete and reliable data of strong (M _s6.5), shallow (h<70 km) earthquakes which occurred in the inner Aegean seismic zone have been utilized to describe its seismicity time variation during 1800–1986 by two independent statistical models. The first is a sequentially stationary model of seismicity rates which shows that intervals of low seismicity rate, lasting for some 37 years, alternate with high rate intervals of 8–12 years duration. The second model is a statistical model according which seismic energy released within 5-year time windows approximates a harmonic curve within a period of about 50 years. This model is in agreement with the notion that the time series of strong earthquake occurrences in the inner Aegean seismic zone consists of a random (shocks withM _s=6.5–6.8) and a nonrandom component (M _s6.9). Maxima and minima of the harmonic curve coincide with the high and low rate intervals, respectively. A model of regional stationary accumulation of thermal stresses along certain seismic belts and their cyclic relaxation may explain this periodicity.     

Calculating the surface temperature of the solid underlyingsurface by modified “Force-Restore” method

Summary  The calculation of rocky surface temperature encounters the problem of unrealistic results due to its large changes at the interface where the energy balance equation is applied. In order to avoid this, we have modified the “force-restore” method into a self-consistent procedure for simultaneous determination of both surface and the deep ground temperature. The approach is applicable to any celestial body where external radiation can be represented by an arbitrary periodic function. The method is tested with Earth based infrared observation on lunar surface temperature and it showed a high level of accuracy and a rather fast convergence of procedure. Received June 28, 1999/Revised June 19, 2000     

Partial pressure and air–sea CO2 flux in the Aegean Sea during February 2006

Data on the distribution of fCO₂ were obtained during a cruise in the Aegean Sea during February 2006. The fCO₂ of surface water (fCO₂^sw) was lower than the atmospheric fCO₂ (fCO₂^atm) throughout the area surveyed and ΔfCO₂ values varied from ?34 to ?61 μatm. The observed under-saturation suggests that surface waters in the Aegean represent a sink for atmospheric CO₂ during the winter of 2006. Higher fCO₂^sw values were recorded in the ‘less warm’ and ‘less saline’ shallow northernmost part of the Aegean Sea implying that the lower seawater temperature and salinity in this area play a crucial role in the spatial distribution of fCO₂^sw.A first estimate of the magnitude of the air–sea CO₂ exchange and the potential role of the Aegean Sea in the transfer of atmospheric CO₂ was also obtained. The air–sea CO₂ fluxes calculated using different gas transfer formulations showed that during February 2006, the Aegean Sea absorbs atmospheric CO₂ at a rate ranging from ?6.2 to ?11.8 mmol m^?2 d^?1 with the shipboard recorded wind speeds and at almost half rate (?3.5 to ?5.5 mmol m^?2 d^?1) with the monthly mean model-derived wind speed. Compared to recent observations from other temperate continental shelves during winter period, the Aegean Sea acts as a moderate to rather strong sink for atmospheric CO₂.Further investigations, including intensive spatial and temporal high-resolution observations, are necessary to elucidate the role of the Aegean Sea in the process of transfer of atmospheric CO₂ into the deep horizons of the Eastern Mediterranean.