Impact of river damming on sediment texture and trace metals distribution along the watershed and the coastal zone of Nestos River (NE Greece)

This study examines the sediment particle size distribution and the trace metal concentrations from a dammed-river watershed (Nestos River) to its deltaic zone in NE Greece. The study area is relatively unpolluted. The distribution of trace metals (Cu, Cr, Cd, Ni, Pb, Hg) in sediments throughout the catchment area showed selective “trapping” of certain elements behind the two artificial dams (Thissavros and Platanovrisi dams) in the watershed and a sudden reduction downstream (83% for Cd, 81% for Cr, 94% for Cu, 90% for Ni, 86% for Hg and 33% for Pb). Marked sediment particle separation is observed at the upstream dam (Thissavros), where coarse material including sand is trapped (coarse fraction 12.9–49.3%). Fine-grained material (<63 μm) is trapped behind the Platanovrisi dam (68.1%), and the reservoir showed elevated metal concentrations, especially for Cu and Cd (16.3 and 0.5 μg/g, respectively). Lead exhibited a homogenous distribution throughout the watershed (20.1–32.3 μg/g). All other trace metals (Cu, Cr, Cd, Ni and Hg) decline sharply downstream of the dam complex. In the delta system, nearshore sediments consist of shallow deposits in the vicinity of river mouth and are enriched in Cr (4.4–53.0 μg/g) and Ni (2.6–44.3 μg/g), while the further offshore and slightly deeper (20–40 m) sediments illustrate elevated Hg (0–0.07 μg/g), Cd (0.09–0.18 μg/g), Cu (11.5–18.3 μg/g) and Ni (38–54.5 μg/g).     

The displacement coefficient method in near‐source conditions

The use of nonlinear static procedures for performance‐based seismic design (PBSD) and assessment is a well‐established practice, which has found its way into modern codes for quite some time. On the other hand, near‐source (NS) ground motions are receiving increasing attention, because they can carry seismic demand systematically different and larger than that of the so‐called ordinary records. This is due to phenomena such as rupture forward directivity (FD), which can lead to distinct pulses appearing in the velocity time‐history of the ground motion. The framework necessary for taking FD into account in probabilistic seismic hazard analysis (PSHA) has recently been established. The objective of the present study is to discuss the extension of nonlinear static procedures, specifically the displacement coefficient method (DCM), with respect to the inelastic demand associated with FD. In this context, a methodology is presented for the implementation of the DCM toward estimating NS seismic demand, by making use of the results of NS‐PSHA and a semi‐empirical equation for NS‐FD inelastic displacement ratio. An illustrative application of the DCM, with explicit inclusion of NS‐pulse‐like effects, is given for a set of typical plane R/C frames designed under Eurocode provisions. Different scenarios are considered in the application and nonlinear dynamic analysis results are obtained and discussed with respect to the static procedure estimates. Conclusions drawn from the results may help to assess the importance of incorporating NS effects in PBSD. Copyright © 2014 John Wiley & Sons, Ltd.     

Seismicity Variations in the Southern Aegean,Greece, Before and After the Large (M7.7) 1956 Amorgos Earthquake Due to Evolving Stress

The largest earthquake (M₀=4.9·10²⁷ dyn·cm) of the 20^th century in the territory of Greece occurred south of Amorgos Island, causing extensive destruction in the southern Aegean area. It occurred on an ENE–trending normal fault that is seated parallel to the Islands southern coastline. Changes in the rates of moderate–size earthquakes (M 5.0) that occurred before and after the Amorgos earthquake, within circular regions centered on its epicenter with radii of 100, 150 and 200 km, are investigated. The rate for moderate–size events just before the main shock appears to be considerably increased when compared to those of either preceding or subsequent periods. Further inspection reveals that more evident seismicity fluctuations are attributed to distances exceeding 100 km. These changes may be indicative of a broad region that is approaching a high stress state prior to an eventual large earthquake. Close to the main event, that is, within the 100–km radius, a remarkable quiescence period lasting about two decades before its occurrence was observed. Changes in seismicity are discussed in combination with static stress changes calculated by the application of the stress evolutionary model that takes into account the coseismic slip associated with the larger events (M 6.5) since the beginning of the 20^th century and the tectonic loading on the major faults in the study area. These larger events, as with the intermediate magnitude seismicity taking place at distances exceeding 100 km and which encircled the quiescent area observed during the last 22 years before the Amorgos earthquake, are well correlated with stress-enhanced areas in each stage of the evolutionary model.     

Improvement in the Determination of the Marine Geoid by Estimating the Bathymetry from Altimetry and Depth Soundings

The contribution of bathymetry to the estimation of gravity field related quantities is investigated in an extended test area in the Mediterranean Sea. The region is located southwest of the island of Crete, Greece, bounded between 33^ˆ ≤ ϕ ≤ 35^ˆ and 15^ˆ ≤ λ ≤ 25^ˆ. Gravity anomalies from the KMS99 gravity field and shipborne depth soundings are used with a priori statistical characteristics of depths in a least-squares collocation procedure to estimate a new bathymetry model. Two different global bathymetry models, namely JGP95E and Sandwell and Smith V8, are used to derive the depth a priori statistical information, while the estimated model is compared against both the global ones and the shipborne depth soundings to assess whether there is an improvement. Various marine geoid models are estimated using ERS1 and GEOSAT Geodetic Mission altimetry and shipborne gravity data. In that process, the effect of the bathymetry is computed using both the estimated and the original depths through a residual terrain modeling reduction. The TOPEX/Poseidon Sea Surface Heights, known for their high accuracy and precision, and the GEOMED solution for the geoid in the Mediterranean are used as control for the validation of the new geoid models and to assess the improvement that the estimated depths offer to geoid modeling. The results show that the newly estimated bathymetry agrees better (by about 30 to 300 m) with the shipborne depth soundings and provides smoother residual geoid heights and gravity anomalies (by about 8-20%) than those from global models. Finally, the achieved accuracy in geoid modeling ranges between 6 and 10 cm (1σ).     

Assessing water erosion in Mediterranean tree crops using GIS techniques and field measurements: the effect of climate change

In this work, a dynamic GIS modeling approach is presented that incorporates: a) geoinformatic techniques, b) 55-year historical meteorological data, and c) field measurements, in order to estimate soil erosion risk in intensively cultivated regions. The proposed GIS-based modeling approach includes the estimation of soil erosion rates due to surface water flow under current and future climate change scenarios A2 and B1 for the years 2030 and 2050. The soil erosion was estimated using the Universal Soil Loss Equation (USLE). The proposed soil erosion model was validated using field measurements at different sites of the study area. The results show that an extended part of the study area is under intense erosion with the mean annual loss to be 4.85 t/ha year⁻¹. Moreover, an increase in rainfall intensity, especially for scenario B1, can generate a significant increase (32.44 %) in soil loss for the year 2030 and a much more (50.77 %) for the year 2050 in comparison with the current conditions. Regarding the scenario A2, a slight decrease (1.85 %) in soil loss was observed for the year 2030, while for 2050 the results show an adequate increase (7.31 %) in comparison with the present. All these approaches were implemented at one of the most productive agricultural areas of Crete in Greece dominated by olive and citrus crops.

     

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.     

An analysis of the linear fixed altimetry-gravimetry boundary-value problem

In this paper the linear fixed altimetry-gravimetry boundary-value problem is analyzed with respect to the existence and uniqueness of the solution. Nowadays, it is possible to determine very precisely points on the physical surface of the Earth by 3D satellite positioning and the problem is to determine the disturbing potential in an unbounded domain representing the exterior of the Earth. In order to establish realistic boundary conditions, a Dirichlet condition is imposed at seas and an oblique derivative condition on land. Then, mathematical methods are used within the frame of functional analysis for attacking the problem under consideration. Specifically, the Stampacchia theorem is used to decide upon the existence and uniqueness of the weak solution of the problem in a weighted Sobolev space. Finally, we confirm that the condition of validity for such a theorem has a geometrical interpretation.     

Earthquake magnitude or seismic moment in seismic hazard evaluation?

Seismic hazard analysis requires the estimation of the probabilities that earthquakes will take place within a region of interest, and the expected level of ground motion which will be received at a site during the nextt years. The earthquake magnitude has been used as a basic parameter, because it is available, under the assumption that the earthquake occurrence is a compound Poisson process with exponential or multinomial distribution of magnitude.For improving the hazard prediction, we used the seismic moment as a basic parameter to estimate the mean rate, , of occurrence of earthquakes in a function of seismic moment rate and slip rate released in a seismogenic region.As an illustration of the model, the seismic hazard analysis at different sites in and around the Gulf of Corinth, central Greece, is presented on the basis of the earthquake magnitude and the seismic moment. Comparison of the results shows that determination of the mean rate of earthquake occurrence, using the conventional Gutenberg-Richter recurrence model, underestimates the seismic hazard at a site.