SANICLAY: simple anisotropic clay plasticity model

SANICLAY is a new simple anisotropic clay plasticity model that builds on a modification of an earlier model with an associated flow rule, in order to include simulations of softening response under undrained compression following K_o consolidation. Non‐associativity is introduced by adopting a yield surface different than the plastic potential surface. Besides, the isotropic hardening of the yield surface both surfaces evolve according to a combined distortional and rotational hardening rule, simulating the evolving anisotropy. Although built on the general premises of critical state soil mechanics, the model induces a critical state line in the void ratio–mean effective stress space, which is a function of anisotropy. To ease interpretation, the model formulation is presented firstly in the triaxial stress space and subsequently, its multiaxial generalization is developed systematically, in a form appropriate for implementation in numerical codes. The SANICLAY is shown to provide successful simulation of both undrained and drained rate‐independent behaviour of normally consolidated sensitive clays, and to a satisfactory degree of accuracy of overconsolidated clays. The new model requires merely three constants more than those of the modified Cam clay model, all of which are easily calibrated from well‐established laboratory tests following a meticulously presented procedure. Copyright © 2006 John Wiley & Sons, Ltd.     

Multicausal environmental severity: A flexible framework for microtidal sandy beaches and the role of polychaetes as an indicator taxon

The aim of this paper is to study the macrofaunal community dynamics and the biological–environmental interactions in the mid- and sublittoral ecosystems of the microtidal Mediterranean sandy shores. Four sandy beaches, three on the island of Crete and one on the northwest coast of Italy were selected to investigate the spatial and temporal changes in the community structure and the associated environmental variables. The littoral zone, which has not been adequately studied in the Eastern Mediterranean, presents special interest not only from the scientific point of view but also for practical reasons of ecological management. The multivariate techniques revealed that the community pattern of the sandy beach macrofauna is mainly spatial rather than temporal. There are pronounced differences in species composition and abundance of the macrofaunal assemblages of the mid- and sublittoral zone. The multicausal environmental severity hypothesis appears to be valid for the sandy beach macrofaunal communities of the Mediterranean. The abundance and composition of the macrofaunal assemblages are highly variable and are affected by the synergistic effects of many environmental variables. The polychaete taxonomic assemblage structure closely follows the macrofaunal community pattern. Differences between the two patterns may arise from the different responses that polychaetes may show to the environmental stress.     

The Ms = 6.2, June 15, 1995 Aigion earthquake (Greece): evidence for low angle normal faulting in the Corinth rift

We present the results of a multidisciplinary study of the M_s = 6.2, 1995, June 15, Aigion earthquake (Gulf of Corinth, Greece). In order to constrain the rupture geometry, we used all available data from seismology (local, regional and teleseismic records of the mainshock and of aftershocks), geodesy (GPS and SAR interferometry), and tectonics. Part of these data were obtained during a postseismic field study consisting of the surveying of 24 GPS points, the temporary installation of 20 digital seismometers, and a detailed field investigation for surface fault break. The Aigion fault was the only fault onland which showed detectable breaks (< 4 cm). We relocated the mainshock hypocenter at 10 km in depth, 38 ° 21.7 N, 22 ° 12.0 E, about 15 km NNE to the damaged city of Aigion. The modeling of teleseismic P and SH waves provides a seismic moment M_o = 3.4 10¹⁸ N.m, a well constrained focal mechanism (strike 277 °, dip 33 °, rake – 77°), at a centroidal depth of 7.2 km, consistent with the NEIC and the revised Harvard determinations. It thus involved almost pure normal faulting in agreement with the tectonics of the Gulf. The horizontal GPS displacements corrected for the opening of the gulf (1.5 cm/year) show a well-resolved 7 cm northward motion above the hypocenter, which eliminates the possibility of a steep, south-dipping fault plane. Fitting the S-wave polarization at SERG, 10 km from the epicenter, with a 33° northward dipping plane implies a hypocentral depth greater than 10 km. The north dipping fault plane provides a poor fit to the GPS data at the southern points when a homogeneous elastic half-space is considered: the best fit geodetic model is obtained for a fault shallower by 2 km, assuming the same dip. We show with a two-dimensional model that this depth difference is probably due to the distorting effect of the shallow, low-rigidity sediments of the gulf and of its edges. The best-fit fault model, with dimensions 9 km E–W and 15 km along dip, and a 0.87 m uniform slip, fits InSAR data covering the time of the earthquake. The fault is located about 10 km east-northeast to the Aigion fault, whose surface breaks thus appears as secondary features. The rupture lasted 4 to 5 s, propagating southward and upward on a fault probably outcropping offshore, near the southern edge of the gulf. In the shallowest 4 km, the slip – if any – has not exceeded about 30 cm. This geometry implies a large directivity effect in Aigion, in agreement with the accelerogram aig which shows a short duration (2 s) and a large amplitude (0.5 g) of the direct S acceleration. This unusual low-angle normal faulting may have been favoured by a low-friction, high pore pressure fault zone, or by a rotation of the stress directions due to the possible dip towards the south of the brittle-ductile transition zone. This fault cannot be responsible for the long term topography of the rift, which is controlled by larger normal faults with larger dip angles, implying either a seldom, or a more recently started activity of such low angle faults in the central part of the rift.     

VHF double-peaked spectra at negative Doppler shifts in the morning sector of radio aurora

This paper presents a new Doppler spectral type of VHP (42 MHz) radio auroral backscatter. This spectrum, which has a double-peaked structure, was observed repeatedly during the morning sector of an exceptionally strong event (A_p = 48) and is due to irregularities moving northwards with quite different velocities. The stronger spectral component, which has a smaller Doppler shift, is centred in frequency at ～?130 Hz, corresponding to the ion-acoustic velocity range in the medium; the weak component, which has a greater frequency shift, usually is centred at about ?300 Hz (～ 1050 m s^?1). Evidence based on spectral analysis of sequential short time sequences shows that the spectral power alternates in time between the two distinct frequency bands where the peaks are located, suggesting that the double-peaked spectrum may result from two competing processes which cannot operate simultaneously. The possibility exists that the theoretical model proposed by Sato (1977), which predicts two different quasi-linear stabilization mechanisms for the two-stream instability, could explain the observed double-peaked spectral type.     

Earthquake generation in Cyprus revealed by the evolving stress field

Strong earthquake occurrence (M ≥ 6.0) onshore and offshore the Cyprus Island constitutes significant seismic hazard because they occur close to populated areas. Seismicity is weak south of the Island along the Cyprean Arc and strong events are aligned along the Paphos transform fault and Larnaka thrust fault zone that were already known and the Lemessos thrust fault zone that defined in the present study. By combining the past history of strong (M ≥ 6.0) events and the long-term tectonic loading on these major fault zones, the evolution of the stress field from 1896 until the present is derived. Although uncertainties exist in the location, magnitude and fault geometries of the early earthquakes included in our stress evolutionary model, the resulting stress field provides an explanation of later earthquake triggering. It was evidenced that the locations of all the strong events were preceded by a static stress change that encouraged failure. The current state of the evolved stress field may provide evidence for the future seismic hazard. Areas of positive static stress changes were identified in the southwestern offshore area that can be considered as possible sites of future seismic activity.     

The Mw = 6.0, 7 September 1999 Athens Earthquake

On 7 September 1999 at 11:56 GMT a destructive earthquake (M_w = 6.0) occurred close to Athens (Greece). The rupture process is examined using data from the Cornet local permanent network, as well as teleseismic recordings. Data recorded by a temporary seismological network were analyzed to study the aftershock sequence. The mainshock was relocated at 38.105°N, 23.565°E, about 20 km northwest of Athens. Four foreshocks were also relocated close to the mainshock. The modeling of teleseismic P and SH waves provides a well-constrained focal mechanism of the mainshock (strike = 105°, dip = 55° and rake = -80°) at a depth of 8 km and a seismic moment M₀ = 1.010²⁵ dyn·cm. The obtained fault plane solution represents normal faulting indicating an almost north-south extension. More than 3500 aftershocks were located, 1813 of which present RMS < 0.1 s and ERH, ERZ < 1.0 km. Two main clusters were distinguished, while the depth distribution is concentrated between 2 and 11 km. Over 1000 fault plane solutions of aftershocks were constrained, the majority of which also correspond to N–S extension. No surface breaks were observed but the fault plane solution of the mainshock is in agreement with the tectonics of the area and with the focal mechanisms obtained by aftershocks. The hypocenter of the mainshock is located on the deep western edge of the fault plane. The relocated epicenter coincides with the fringe that represents the highest deformation observed on the differential interferometric image. The calculated source duration is 5 sec, while the estimated dimensions of the fault are 15 km length and 10 km width. The source process is characterized by unilateral eastward rupture propagation, towards the city of Athens. An evident stop phase observed in the recordings of the Cornet local stations is interpreted as a barrier caused by the Aegaleo Mountain.     

A critical state interpretation for the cyclic liquefaction resistance of silty sands

Contrary to many laboratory investigations, common empirical correlations from in situ tests consider that the increase in the percentage of fines leads to an increase of the cyclic liquefaction resistance of sands. This paper draws upon the integrated Critical State Soil Mechanics framework in order to study this seemingly not univocal effect. Firstly the effect of fines on the Critical State Line (CSL) is studied through a statistical analysis of a large data set of published monotonic triaxial tests. The results show that increasing the content of non-plastic fines practically leads to a clockwise rotation of the CSL in (e–ln p) space. The implication of this effect on cyclic liquefaction resistance is subsequently evaluated with the aid of a properly calibrated critical state elasto-plastic constitutive model, as well as a large number of published experimental results and in situ empirical correlations. Both sets of data show clearly that a fines content, less than about 30% by weight, may prove beneficial at relatively small effective stresses (p₀<50–70 kPa), such as the in situ stresses prevailing in most liquefaction case studies, and detrimental at larger confining stresses, i.e. the stresses usually considered in laboratory tests. To the extent of these findings, a correction factor is proposed for the practical evaluation of liquefaction resistance in terms of the fines content and the mean effective confining stress.     

Intrinsic groundwater vulnerability determination at the aquifer scale: a methodology coupling travel time estimation and rating methods

Groundwater vulnerability has been subject of much research due to the valuable information it provides concerning groundwater protection and exploitation potential. Up to now, most groundwater vulnerability studies adopt subjective systems of rating the various factors and subsequently, their results are often ambiguous and contradicting. Within the present study a methodology for the estimation of intrinsic groundwater vulnerability at the aquifer scale is presented. The methodology is based on travel time estimation from specified sources of pollution to the aquifer. Besides the deterministic calculation of travel times, the methodology provides a rating system for each pollution source, based on its relative severity and the estimated threat that it poses to the aquifer. Therefore, it can be regarded as a hybrid method that couples the advantages provided by the physically based methods with those of the subjective rating systems. The methodology is applied to the Neon Sidirochorion aquifer, Northeastern Greece, an overexploited aquifer where river waters, sea waters and lake waters interact, causing groundwater quality deterioration to the aquifer. The results indicated that the proposed groundwater vulnerability assessment methodology is well capturing pollution related to saltwater intrusion and agricultural activities, while it is concluded that the conceptual model is significantly affecting the vulnerability assessment results and therefore has to be previously developed.     

Testing optically stimulated luminescence dating on sand-sized quartz of deltaic deposits from the Sperchios delta plain,central Greece

This study reports on the first investigation into the potential of luminescence dating to establish a chronological framework for the depositional sequences of the Sperchios delta plain, central Greece. A series of three borehole cores(20 m deep) and two shallow cores(4 m deep), from across the delta plain, were extracted, and samples were collected for luminescence dating. The luminescence ages of sand-sized quartz grains were obtained from small aliquots of quartz, using the Single-Aliquot Regenerative-dose(SAR) protocol.The equivalent dose determination included a series of tests and the selection of the Minimum Age Model(MAM) as the most appropriate statistical model. This made it possible to confirm the applicability of quartz Optically Stimulated Luminescence(OSL) dating to establish absolute chronology for deltaic sediments from the Sperchios delta plain.Testing age results of the five cores showed that the deltaic sediments were deposited during the Holocene.A relatively rapid deposition is implied for the top ~14 m possibly as a result of the deceleration in the rate of the sea-level rise and the transition to terrestrial conditions, while on the deeper parts, the reduced sedimentation rate may indicate a lagoonal or coastal environment.