Engineering geological studies and tunnel support design at Sulakyurt dam site, Turkey

This paper presents the results of engineering geological investigations and tunnel support design studies, carried out at the Sulakyurt dam site, northeast of Ankara, Turkey. The Sulakyurt dam will be used for flow control and water storage for irrigation projects. Studies were carried out both in the field and the laboratory. Field studies include engineering geological mapping, intensive discontinuity surveying, core drilling and sampling for laboratory testing. The diversion tunnel will be driven in rock mass, consisting of granite and diorite. Empirical, analytical and numerical methods were combined for safe tunnel design. Rock mass rating (RMR), Rock mass quality (Q) and Geological strength index (GSI) systems were used for empirical rock mass quality determination, site characterization and support design. The convergence–confinement method was used as analytical method and software called Phase², a 2D finite element program, was utilized as numerical method. According to the results acquired from the empirical, analytical and numerical methods, tunnel stability problems were expected in both granite and diorite rock masses. The support system, suggested by empirical methods, was applied and the performance of suggested support system was evaluated by means of numerical modelling. It was concluded that the suggested support systems were adequate, since after applying the suggested support system to granite and diorite, tunnel deformation and the yielded elements around the tunnel decreased significantly. Thus, it is suggested that for more reliable support design empirical, numerical and analytical methods should be combined.     

Testing earthquake prediction methods: «The West Pacific short-term forecast of earthquakes with magnitude MwHRV ≥ 5.8»

Analyzing the tables and probability maps posted by Yan Y. Kagan and David D. Jackson in April 2002–September 2004 at http://scec.ess.ucla.edu/~ykagan/predictions_index.html and the catalog of earthquakes for the same period, the conclusion is drawn that the underlying method could be used for prediction of aftershocks, while it does not outscore random guessing when main shocks are considered.     

A laboratory experiment to monitor the contact state of a fault by transmission waves

We performed a series of laboratory experiments in which elastic waves were transmitted across a simulated fault. Two types of experiments were carried out: (1) Normal Stress Holding Test (NSHT): normal stress was kept constant for about 3 h without shear stress and transmission waves were observed. (2) Shear Stress Increasing Test (SSIT): shear stress was gradually increased until a stick-slip event occurred. Transmission waves were continuously observed throughout the process of stress accumulation. We focused on the change in transmission waves during the application of shear stress and especially during precursory slips.It was found in NSHT that the amplitude of transmission waves linearly increased with the logarithm of stationary contact time. The increase amounted to a few percent after about 3 h. Creep at asperity contacts is responsible for this phenomenon. From a theoretical consideration, it was concluded that the real contact area increased with the logarithm of stationary contact time.We observed in SSIT a significant increase in wave amplitude with shear stress application. This phenomenon cannot be attributed to the time effect observed in NSHT. Instead, it can be explained by the mechanism of “junction growth” proposed by Tabor. Junction growth yields an increase in real contact area. It is required for junction growth to occur that the material in contact is already plastic under a purely normal loading condition. A computer simulation confirmed that this requirement was satisfied in our experiments. We also found that the rate at which the amplitude increased was slightly reduced prior to a stick-slip event. The onset time of the reduction well coincides with the onset of precursory slip. The cause of the reduction is attributed to the reset of stationary contact time due to displacement. This interpretation is supported by the result of NSHT. Taking the time of stationary contact in SSIT into account, we may expect the change in wave amplitude to be, at most, only a few percent. The observed slight reduction in increasing rate is, in this sense, reasonable. The static stiffness of the fault also decreases with precursory slip. It was also found that low frequency waves are a better indicator of precursory slip than high frequency waves. This might suggest that low frequency waves with longer wavelength are a better indicator of average behavior of faults. The problem, however, merits a further investigation. The shifts in phase were also found to be a good indicator of the change in contact state of the fault. The changes in both amplitude and phase of transmission waves are unifyingly understood through the theory of transmission coefficient presented by Pyrak-Nolte et al. Rough surfaces have a tendency to give larger stick-slips than smooth surfaces. The amount of precursory slip is larger for rough surfaces than for smooth surfaces. Although it was confirmed by a computer simulation that rough surfaces have larger contact diameters than smooth surfaces, the rigorous relationship between the surface roughness (contact diameter) and the amount of precursory slips was not established.     

A way to synchronize models with seismic faults for earthquake forecasting: Insights from a simple stochastic model

Numerical models are starting to be used for determining the future behaviour of seismic faults and fault networks. Their final goal would be to forecast future large earthquakes. In order to use them for this task, it is necessary to synchronize each model with the current status of the actual fault or fault network it simulates (just as, for example, meteorologists synchronize their models with the atmosphere by incorporating current atmospheric data in them). However, lithospheric dynamics is largely unobservable: important parameters cannot (or can rarely) be measured in Nature. Earthquakes, though, provide indirect but measurable clues of the stress and strain status in the lithosphere, which should be helpful for the synchronization of the models.The rupture area is one of the measurable parameters of earthquakes. Here we explore how it can be used to at least synchronize fault models between themselves and forecast synthetic earthquakes. Our purpose here is to forecast synthetic earthquakes in a simple but stochastic (random) fault model. By imposing the rupture area of the synthetic earthquakes of this model on other models, the latter become partially synchronized with the first one. We use these partially synchronized models to successfully forecast most of the largest earthquakes generated by the first model. This forecasting strategy outperforms others that only take into account the earthquake series. Our results suggest that probably a good way to synchronize more detailed models with real faults is to force them to reproduce the sequence of previous earthquake ruptures on the faults. This hypothesis could be tested in the future with more detailed models and actual seismic data.     

A multi-method approach to study the stability of natural slopes and landslide susceptibility mapping

In this paper, a multi-method approach for the assessment of the stability of natural slopes and landslide hazard mapping applied to the Dakar coastal region is presented. This approach is based on the effective combination of geotechnical field and laboratory works, of GIS, and of mechanical (deterministic and numerical) stability analysis. By using this approach, valuable results were gained regarding instability factors, landslide kinematics, simulation of slope failure and coastal erosion. This led to a thorough assessment and strong reduction in the subjectivity of the slope stability and hazard assessment and to the development of an objective landslide danger map of the SW coast of Dakar. Analysis of the results shows that the slides were influenced by the geotechnical properties of the soil, the weathering, the hydrogeological situation, and the erosion by waves. The landslide susceptibility assessment based on this methodological approach has allowed for an appropriate and adequate consideration of the multiple factors affecting the stability and the optimization of planning and investment for land development in the city.     

Towards Formulation of a Space-borne System for Early Warning of Floods: Can Cost-Effectiveness Outweigh Prediction Uncertainty?

The three most important components necessary for functioning of an operational flood warning system are: (1) a rainfall measuring system; (2) a soil moisture updating system; and, (3) a surface discharge measuring system. Although surface based networks for these systems can be largely inadequate in many parts of the world, this inadequacy particularly affects the tropics, which are most vulnerable to flooding hazards. Furthermore, the tropical regions comprise developing countries lacking the financial resources for such surface-based monitoring. The heritage of research conducted on evaluating the potential for measuring discharge from space has now morphed into an agenda for a mission dedicated to space-based surface discharge measurements. This mission juxtaposed with two other upcoming space-based missions: (1) for rainfall measurement (Global Precipitation Measurement, GPM), and (2) soil moisture measurement (Hydrosphere State, HYDROS), bears promise for designing a fully space-borne system for early warning of floods. Such a system, if operational, stands to offer tremendous socio-economic benefit to many flood-prone developing nations of the tropical world. However, there are two competing aspects that need careful assessment to justify the viability of such a system: (1) cost-effectiveness due to surface data scarcity; and (2) flood prediction uncertainty due to uncertainty in the remote sensing measurements. This paper presents the flood hazard mitigation opportunities offered by the assimilation of the three proposed space missions within the context of these two competing aspects. The discussion is cast from the perspective of current understanding of the prediction uncertainties associated with space-based flood prediction. A conceptual framework for a fully space-borne system for early-warning of floods is proposed. The need for retrospective validation of such a system on historical data comprising floods and its associated socio-economic impact is stressed. This proposal for a fully space-borne system, if pursued through wide interdisciplinary effort as recommended herein, promises to enhance the utility of the three space missions more than what their individual agenda can be expected to offer.     

Verification of the motion fields over Southeast Asia by the Florida State University (FSU) superensemble forecast

The Florida State University (FSU) multimodel superensemble forecast is evaluated against several other operational weather models for the Southeast Asia region. The superensemble technique has demonstrated its exceptional skills in forecasting precipitation, motion and mass fields compared to either individual global operational or ensemble mean forecasts. The motion field investigation for the season of 2001 reveals that the superensemble forecasts are closer to the observed data compared to the other global member operational models through its low systematic errors at the 850 hPa level. The FSU multimodel superensemble forecasts exhibit the lowest root mean square errors (RSMEs), the highest correlation against the best observed data and the lowest systematic errors compared to the other operational model members. These forecasts have the potential to provide better daily weather predictions over the Southeast Asia region, particularly during the early northeast monsoon that often causes heavy rainfall in the equatorial part of the Southeast Asia region.     

Natural and man-induced stress evolution of slopes: the Monte Mario hill in Rome

The paper deals with stress-release effects induced by man-made cuts or excavations into natural stiff clay slopes that experienced erosion in response to valley deepening. The study was focused on the Monte Mario hill in Rome (Italy), which formed part of an area of recent urban expansion. The methodology of the study relied on a reference engineering-geology model, which was developed on the basis of site and laboratory data and stress–strain analyses. The latter analyses were carried out with the finite-difference code FLAC 4.0. Numerical modelling was based on a sequential approach, taking into account the main evolutionary stages of the Tiber river valley in Romeȁ9s urban area and then making cuts at the bottom of the slope located south of the Monte Mario Astronomical Observatory. The simulation revealed the stress-release effects that fluvial erosion and excavation fronts have caused on the investigated slopes and their consequent gravitational instabilities. These processes appear with metre-scale displacements, followed by stress-release cracks (actually observed on the slopes under review). In quantifying stress-release deformations, the simulation took into account the possible role of creep in the observed retardation of stress-release effects.Research activities were carried on with the co-operation of Dȁ9Arcangelo A. (Consorzio TREESSE, andrea.darcangelo@libero.it) and Moretti S. (IMG S.r.l. Servizi Tecnici per lȁ9Ingegneria e lȁ9Ambiente, moretti@img-srl.it)     

Occurrence of salt water above fresh water in dynamic equilibrium in a coastal groundwater flow system near De Panne, Belgium

A salt water lens is found above fresh water under the shore between Dunkerque (France) and Nieuwpoort (Belgium). This inverse density distribution is in a dynamic equilibrium. It develops due to the infiltration of salt water on the back shore during high tide. Under this salt water lens, water infiltrated in the adjacent dune area flows towards the sea and discharges at the seabed. This water quality distribution differs from the classic salt water wedge under fresh water described in the literature. Here, the evolution to this water quality distribution is simulated with a density dependent numerical model. A large tidal range, shore morphology and a permeable groundwater reservoir are the main conditions for the observed water quality distribution. By altering these conditions, intermediate water quality distributions between the classic salt water wedge and the one discussed here develop. Based on these simulations, it is expected that similar kinds of inverse density distribution could be present in a number of coastal areas, which have tides, a gently sloping shore and a permeable substratum.