Ground water recharge and flow characterization using multiple isotopes

Stable isotopes of delta(18)O, delta(2)H, and (13)C, radiogenic isotopes of (14)C and (3)H, and ground water chemical compositions were used to distinguish ground water, recharge areas, and possible recharge processes in an arid zone, fault-bounded alluvial aquifer. Recharge mainly occurs through exposed stream channel beds as opposed to subsurface inflow along mountain fronts. This recharge distribution pattern may also occur in other fault-bounded aquifers, with important implications for conceptualization of ground water flow systems, development of ground water models, and ground water resource management. Ground water along the mountain front near the basin margins contains low delta(18)O, (14)C (percent modern carbon [pmC]), and (3)H (tritium units [TU]), suggesting older recharge. In addition, water levels lie at greater depths, and basin-bounding faults that locally act as a flow barrier may further reduce subsurface inflow into the aquifer along the mountain front. Chemical differences in ground water composition, attributed to varying aquifer mineralogy and recharge processes, further discriminate the basin-margin and the basin-center water. Direct recharge through the indurated sandstones and mudstones in the basin center is minimal. Modern recharge in the aquifer is mainly through the broad, exposed stream channel beds containing coarse sand and gravel where ground water contains higher delta(18)O, (14)C (pmC), and (3)H (TU). Spatial differences in delta(18)O, (14)C (pmC), and (3)H (TU) and occurrences of extensive mudstones in the basin center suggest sluggish ground water movement, including local compartmentalization of the flow system.     

Sr–Nd isotopes and geochemistry of the infrastructural rocks in the Meatiq and Hafafit core complexes, Eastern Desert, Egypt: Evidence for involvement of pre-Neoproterozoic crust in the growth of Arabian–Nubian Shield

Abstract Meatiq and Hafafit core complexes are large swells in the Eastern Desert of Egypt, comprising two major tectono‐stratigraphic units or tiers. The lower (infrastructure) unit is composed of variably cataclased gneissose granites and high‐grade gneisses and schists. It is structurally overlain by Pan–African ophiolitic mélange nappes (the higher unit). The two units are separated by a low‐angle sole thrust, along which mylonites are developed. Major and trace element data indicate formation of the gneissose granites in both volcanic arc and within‐plate settings. Nevertheless, all analyzed gneissose granites and other infrastructural rocks, exhibit low initial ratios (Sr_i) (<0.7027), positive ε_Nd(t) (+4.9 to +10.3) and Neoproterozoic Nd model age (T_DM) (592–831 Ma for the gneissose granite samples). Although these values are compatible with other parts of the Arabian– Nubian Shield considered to be juvenile, the ε_Nd(t) values and several incompatible element ratios of the gneissose granites are too low to be derived from a mantle source without contribution from an older continental crust. Our geological, Sr–Nd isotopic and chemical data combined with the published zircon ages indicate the existence of a pre‐Neoproterozoic continent in the Eastern Desert that started to break up at ca 800 Ma. Rifting and subsequent events caused the formation of oceanic crust and emplacement within‐plate alkali basalts in the hinterland domains of the old continent. The emplacement of basaltic magma might have triggered melting of lower crust in the old continent and resulted in emplacement of the within‐plate granite masses between 700 Ma and 626 Ma. The granite masses and other rocks in the old continent have been subjected to deformation during the over‐thrusting of Pan–African nappes, probably because of the oblique convergence between East and West Gondwanaland. Rb–Sr isotopes of the gneissose granites in both Meatiq and Hafafit core complexes defines an isochron age of 619 ± 25 Ma with Sr_i of 0.7009 ± 0.0017 and mean squares of weighted deviates = 2.0. We interpret this age as the date of thrusting of the Pan–African nappes in the Eastern Desert. Continued oblique convergence between East and West Gondwanaland could have resulted in the formation northwest–southeast‐trending Meatiq and Hafafit anticlinoriums.     

Analysis of ambient noise in Yalova, Turkey: discrimination between artificial and natural excitations

Ambient noise measurements acquired in Yalova, which was highly damaged during the 1999 Izmit earthquake, are analyzed to explore the site characteristics. The region of Yalova is governed by complex geological and geomorphological structures consisting of river beds extending from the mountains to the sea, ridges between them, plains in front of them with different size, and the sea coast. As a result of these shallow geological features, the H/V curves exhibit complex patterns. Clear peaks in the H/V curves, which can be interpreted as reliable site resonance frequency, are observed only at about half of the measurement sites. At the remaining sites industrial peaks, broad peaks, or flat responses dominate the spectral ratio graphs. We observed that man-made noises generated by marble cutting machines in Hersek delta mask the site resonance frequencies or can be misinterpreted as a resonance frequency. In total, we identified three anthropogenic noise sources at fundamental frequencies of 1.3, 1.5, and 1.7 Hz along with their two- and threefold harmonics. The parts of H/V curves showing unusual low scattering can be a clue to identify anthropogenic effects. In the assessment of H/V curves, the site location and the similarity of the near surface geology were taken into account. The Laledere plain with thick and soft sediment sequence surprisingly displays flat responses due to a possible low impedance contrast. The Ciftlikkoy and Hacimehmet plains exhibit clear resonance peaks at nearly 1 Hz possessing the largest amplitudes. These sites experienced the highest damage in Yalova during the Izmit earthquake. In contrast, the Cinarcik region which was also exposed to high damage, do not show any obvious amplifications on the H/V curves. Generally, the H/V curves for valley and ridge sites in Yalova reveals a resonance peak at approximately 1 Hz and almost flat curves, respectively. However, several sites on the ridges and valleys portray different patterns.     

Monthly stream flow forecasting via dynamic spatio-temporal models

In this research, a dynamic linear spatio-temporal model (DLSTM) was developed and evaluated for monthly streamflow forecasting. For parameter estimation, coupled expectation–maximization (EM) algorithm and Kalman filter was adopted. This combination enables the model to estimate the state vector and parameters concurrently. Different forecast scenarios including various combinations of upstream stations were considered for downstream station streamflow forecasting. Several statistical criteria, nonparametric and visual tests were used for model evaluation. Results indicated that the spatio-temporal model performed acceptably in almost all scenarios. The dynamic model was able to capitalize on coupled spatial and temporal information provided that there is spatial connectivity in the studied hydrometric stations network. Moreover, threshold level method was used for model evaluation in drought and wet periods. Results indicated that, in validation phase, the model was able to forecast the drought duration and volume deficit/over threshold, although volume deficit/over threshold could not be accurately simulated.     

Regional gravimetric quasi-geoid model and transformation surface to national height system for Turkey (THG-09)

Turkish regional geoid models have been developed by employing a reference earth gravitational model, surface gravity observations and digital terrain models. The gravimetric geoid models provide a ready transformation from ellipsoidal heights to the orthometric heights through the use of GPS/leveling geoid heights determined through the national geodetic networks. The recent gravimetric models for Turkish territory were computed depending on OSU91 (TG-91) and EGM96 (TG-03) earth gravitational models. The release of the Earth Gravitational Model 2008 (EGM08), the collection of new surface gravity observations, the advanced satellite altimetry-derived gravity over the sea, and the availability of the high resolution digital terrain model have encouraged us to compute a new geoid model for Turkey. We used the Remove-Restore procedure based on EGM08 and applied Residual Terrain Model (RTM) reduction of the surface gravity data. Fast Fourier Transformation (FFT) was then used to obtain the residual quasigeoid from the reduced gravity. We restored the individual contributions of EGM08 and RTM to the whole quasi-geoid height (TQG-09). Since the Helmert orthometric height system is adopted in Turkey, the quasi-geoid model (TQG-09) was then converted to the geoid model (TG-09) by making use of Bouguer gravity anomalies and digital terrain model. After all we combined a gravimetric geoid model with GPS/leveling geoid heights in order to obtain a hybrid geoid model (THG-09) (or a transformation surface) to be used in GPS applications. The RMS of the post-fit residuals after the combination was found to be ± 0.95 cm, which represents the internal precision of the final combination. And finally, we tested the hybrid geoid model with GPS/leveling data, which were not used in the combination, to assess the external accuracy. Results show that the external accuracy of the THG-09 model is ± 8.4 cm, a precision previously not achieved in Turkey until this study.     

Numerical modelling of organic waste dispersion from fjord located fish farms

In this study, a three-dimensional particle tracking model coupled to a terrain following ocean model is used to investigate the dispersion and the deposition of fish farm particulate matter (uneaten food and fish faeces) on the seabed due to tidal currents. The particle tracking model uses the computed local flow field for advection of the particles and random movement to simulate the turbulent diffusion. Each particle is given a settling velocity which may be drawn from a probability distribution according to settling velocity measurements of faecal and feed pellets. The results show that the maximum concentration of organic waste for fast sinking particles is found under the fish cage and continue monotonically decreasing away from the cage area. The maximum can split into two maximum peaks located at both sides of the centre of the fish cage area in the current direction. This process depends on the sinking time (time needed for a particle to settle at the bottom), the tidal velocity and the fish cage size. If the sinking time is close to a multiple of the tidal period, the maximum concentration point will be under the fish cage irrespective of the tide strength. This is due to the nature of the tidal current first propagating the particles away and then bringing them back when the tide reverses. Increasing the cage size increases the likelihood for a maximum waste accumulation beneath the fish farm, and larger farms usually means larger biomasses which can make the local pollution even more severe. The model is validated by using an analytical model which uses an exact harmonic representation of the tidal current, and the results show an excellent agreement. This study shows that the coupled ocean and particle model can be used in more realistic applications to help estimating the local environmental impact due to fish farms.     

Convergence-control approach for rock tunnels reinforced by grouted bolts,using the homogenization concept

The scope of this paper is to introduce a method for the analysis of rock tunnels reinforced by grouted bolts, based on the convergence-control approach. The analytical formulations presented in this paper refers to an elasto-plastic behavior of the rock mass and the latest Hoek and Brown yield criterion (Version 2002). In order to model the reinforced plastic zone, the equivalent material approach was taken into account such that the apparent strength of the rock mass is improved as a consequence of the bolting effect. The general design guides and examples presented are intended to facilitate the comprehension and application of the proposed analytical solution in practice.     

A simple and efficient approach to predict reservoir settling volume: case study of Bin El Ouidane reservoir (Morocco)

The negligence of the cumulative deposit volume over years may lead to uselessness of a pre-established and efficient reservoir management method. In this objective, this paper tries to introduce a simple approach to forecast annually the volume lost from the reservoir capacity. The suggested approach employs the modified universal soil loss equation (MUSLE) and the Dendy and Bolton (J Soil Water Conserv 31:264–266, 1976) methods, in addition to the sediment consolidation process. This process permits to follow the evolution of sediment deposit volume throughout the years. For application, we selected the reservoir of Bin El Ouidane (Morocco), for which we have tried to fit a linear regression between sediment yield and observed settling volumes. The optimization process was assured by the algorithm genetic and the OptQuest methods. The results show a high significant (R²?=?0.9513, RMSE?=?8.5428, and p value?=?0.0009) linear regression between cumulative sediment yield volume and accumulated measured lost volume from the reservoir of Bin El Ouidane. The use of a linear relation is justified by the fact that all nonlinear regressions in long term are composed of multiple linear regressions in short term. Nevertheless, this method can become more efficient, and confirmed for real forecasting applications by performing modifications essentially related to in situ measures of different variables. Finally, and considering its simplicity regarding input data and application, the proposed approach converges to efficient results. Therefore, and after validation in other reservoirs, this method can be used to forecast annual reservoir settling at short term.