Discussion of “Interfacing the geographic information system,remote sensing,and the soil conservation service-curve number method to estimate curve number and runoff volume in the ASIR region of Saudi Arabia” by Fawzi S. Mohammad,Jan Adamowski

The soil conservation service curve number (SCS-CN) method is one of the most commonly used methods to compute the direct runoff from a rainfall event. Since the method was established, numerous researches were undertaken to improve the method through accurate estimation of its parameter and especially the curve number (CN). However, the essence of the SCS method, as an event-based Hortonian mechanism method, remained unchanged. The main assumption of the method related to the rainfall input is that the rainfall is continuous in time and uniform over the watershed. Mohammad and Adamowski (2015) paper apparently used the SCS method to estimate the annual runoff using the annual rainfall as one cumulative rainfall input value, which is a violation of the event-based principle of the method and of the assumption of the continuity of the rainfall event.To re-estimate the average annual runoff more realistically for the Asir region, Saudi Arabia, daily rainfall data from 14 rainfall stations are used for calculating the resulting runoff depths, on a daily event-by-event rainfall basis, throughout the whole simulation period. The resulting runoff depths are added for each year, and the total cumulative annual runoff values for each year are averaged to get the average annual runoff. The runoff values based on the previously mentioned procedure are an upper limit of the actual average annual runoff as the underlying SCS equations discard evaporation and similar long-term losses. Nevertheless, the average runoff values obtained in the discussion paper are an order of magnitude (at least five to tenfold) lower than the ones of the original paper. An equation is proposed to obtain a more realistic estimate of the average annual runoff, to be used with the average annual rainfall as an input, if the annual value is the only available rainfall information.     

Dead Sea Basin Imaged by Ambient Seismic Noise Tomography

In the framework of the Dead Sea Integrated Research project (DESIRE), 59 seismological stations were deployed in the region of the Dead Sea Basin. Twenty of these stations recorded data of sufficiently high quality between May and September 2007 to be used for ambient seismic noise analysis. Empirical Green’s functions are extracted from cross-correlations of long term recordings. These functions are dominated by Rayleigh waves, whose group velocities can be measured in the frequency range from 0.1 to 0.5 Hz. Analysis of positive and negative correlation lags of the Green’s functions makes it possible to identify the direction of the source of the incoming energy. Signals with frequencies higher than 0.2 Hz originate from the Mediterranean Sea, while low frequencies arrive from the direction of the Red Sea. Travel times of the extracted Rayleigh waves were measured between station pairs for different frequencies, and tomographically inverted to provide independent velocity models. Four such 2D models were computed for a set of frequencies, all corresponding to different sampling depths, and thus together giving an indication of the velocity variations in 3D extending to a depth of 10 km. The results show low velocities in the Dead Sea Basin, consistent with previous studies suggesting up to 8 km of recent sedimentary infill in the Basin. The complex structure of the western margin of the Basin is also observed, with sedimentary infill present to depths not exceeding 5 km west of the southern part of the Dead Sea. The high velocities associated with the Lisan salt diapir are also observed down to a depth of ~5 km. The reliability of the results is confirmed by checkerboard recovery tests.     

利用瞬变电磁和多频电磁方法监测地下水位及其对Hawara遗址的影响

利用瞬变电磁（TEM）和地面电磁（GEM）方法,本文研究了埃及Hawara遗址区的地下水对遗址的影响.这个遗址包括Hawara金字塔、北区墓地、以及被称作“迷宫”的南区墓葬场所.现今Hawara遗址已经完全荒废,被耕作区和Bahr Wahba运河包围.地表水和地下水会破坏Hawara金字塔和“迷宫”遗址的地基,现在金字塔入口已被淹没至地面下约6 m深处.本文在耕作区、金字塔以及“迷宫”等三个地区进行了TEM勘察;在“迷宫”地区进行了GEM勘察.综合分析以上两种勘察数据,我们发现,农业灌溉是浸入到地下的水的主要来源;本地区地下水位随着观测点及海拔的不同,在地面以下2~7 m之间变化.     

Neoproterozoic nascent island arc volcanism from the Nubian Shield of Egypt: Magma genesis and generation of continental crust in intra-oceanic arcs

The Neoproterozoic Wadi Ranga metavolcanic rocks, South Eastern Desert of Egypt, constitute a slightly metamorphosed bimodal sequence of low-K submarine tholeiitic mafic and felsic volcanic rocks. The mafic volcanic rocks are represented by massive and pillow flows and agglomerates, composed of porphyritic and aphyric basalts and basaltic andesites that are mostly amygdaloidal. The felsic volcanic rocks embrace porphyritic dacites and rhyolites and tuffs, which overlie the mafic volcanic rocks. The geochemical characteristics of Wadi Ranga volcanic rocks, especially a strong Nb depletion, indicate that they were formed from subduction-related melts. The clinopyroxene phenocrysts of basalts are more akin to those crystallizing from island-arc tholeiitic magmas. The tholeiitic nature of the Wadi Ranga volcanics as well as their LREE-depleted or nearly flat REE patterns and their low K₂O contents suggest that they were developed in an immature island arc setting. The subchondritic Nb/Ta ratios (with the lowest ratio reported for any arc rocks) and low Nb/Yb ratios indicate that the mantle source of the Wadi Ranga mafic volcanic rocks was more depleted than N-MORB-source mantle. Subduction signature was dominated by aqueous fluids derived from slab dehydration, whereas the role of subducted sediments in mantle-wedge metasomatization was subordinate, implying that the subduction system was sediment-starved and far from continental clastic input. The amount of slab-derived fluids was enough to produce hydrous magmas that follow the tholeiitic but not the calc-alkaline differentiation trend. With Mg# > 64, few samples of Wadi Ranga mafic volcanic rocks are similar to primitive arc magmas, whereas the other samples have clearly experienced considerable fractional crystallization.The low abundances of trace elements, together with low K₂O contents of the felsic metavolcanic rocks indicate that they were erupted in a primitive island arc setting. The felsic volcanic rocks are characterized by lower K/Rb ratios compared to the mafic volcanic rocks, higher trace element abundances (~ 2 to ~ 9 times basalt) on primitive arc basalt-normalized pattern and nearly flat chondrite-normalized REE patterns, which display a negative Eu anomaly. These features are largely consistent with fractional crystallization model for the origin of the felsic volcanic rocks. Moreover, SiO₂-REE variations for the Wadi Ranga volcanic rocks display steadily increasing LREE over the entire mafic to felsic range and enriched La abundances in the felsic lavas relative to the most mafic lavas, features which are consistent with production of the felsic volcanic rocks through fractional crystallization of basaltic melts. The relatively large volume of Wadi Ranga silicic volcanic rocks implies that significant volume of silicic magmas can be generated in immature island arcs by fractional crystallization and indicates the significant role of intra-oceanic arcs in the production of Neoproterozoic continental crust. We emphasize that the geochemical characteristics of these rocks such as their low LILE and nearly flat REE patterns can successfully discriminate them from other Egyptian Neoproterozoic felsic volcanic rocks, which have higher LILE, Zr and Nb and fractionated REE patterns.     

Use of remote sensing, GIS and groundwater monitoring to estimate artificial groundwater recharge in Riyadh, Saudi Arabia

This study aims to estimate artificial recharge of groundwater by using remote sensing technology, geographical information systems, and groundwater surveys. This study is part of the King Fahd project for rainfall and runoff water harvesting, within the premises of Alilb Dam in Diriyah to the west of Riyadh. Digital elevation models were obtained with the help of aerial photography from the year 2007. These models were used to delineate watershed. Average rainfall was calculated using isoheytal method, and the area of each of the storage lakes was measured using SPOT 5 satellite images from 2007. Fluctuations in groundwater levels, evaporation, and infiltration rates were used to determine the water balance for the purpose of estimating of artificial recharge. Artificial recharge rates were found to surpass natural recharge from rainfall. Recharge wells caused a reduction in the effect of evaporation on storage lakes and helped in supplying water to the groundwater reservoir. Moreover, 80% and 86% of the rainwater was found to be available for artificial recharge in Alilb at 2005 and 2007, respectively. The study recommends the establishment of strategic projects for water storage using artificial recharge wells, an increase in the number of monitoring wells around the dams, and the monitoring of hydrochemical changes in groundwater both before and after the artificial recharge. It also recommends the erection of a weather station in the northwest of Wadi Hanifa.     

Neoproterozoic contaminated MORB of Wadi Ghadir ophiolite,NE Africa: Geochemical and Nd and Sr isotopic constraints

The ophiolitic metabasalts (pillowed and sheeted dikes) of Wadi Ghadir area, Eastern Desert, Egypt, were analyzed for their major, trace and rare earth elements, Nd and Sr isotopes and the chemistry of their plagioclase, amphibole and chlorite was also reported. Geochemically these rocks range from tholeiitic basalt to basaltic andesite. The generally low MgO, Cr and Ni and high Zr contents are consistent with derivation of these rocks from an evolved magma. The high TiO₂ contents (mostly between 1.76% and 2.23%) classify Wadi Ghadir ophiolitic metabasalts as MORB ophiolite. The chondrite-normalized REE patterns of most samples display small LREE-enrichment with (La/Yb)_n ranging from 1.44 to 2.56. The MORB-normalized spider diagram shows variable LILE abundances, which are either similar to or enriched relative to MORB, and most samples display small Nb depletion. The abundances of some LILE (Ba, Rb and K) as well as Na and Si were modified by post-magmatic seafloor hydrothermal alteration. Enrichment of the least mobile LILE (Th & U) indicates that Wadi Ghadir ophiolitic metabasalts are akin to C (contaminated)-MORB. These geochemical characteristics are similar to BABB modified by contamination.Wadi Ghadir metabasalts have low initial Sr ratios (0.7010–0.7034) which are similar to those of MORB, while their ε_Nd(t) values (+7.7–+4.5) are either more or less positive than the value of depleted mantle (DM). The more positive ε_Nd(t) values indicate DM source for these basalts, while the less positive ε_Nd(t) values reflect the involvement of slightly older component in Wadi Ghadir ophiolite. We suggest that the parent magma of metabasalts was contaminated by slightly older material, most probably oceanic-arc crustal rocks, which caused enrichment in LREE, and by analogy LILE, but did not significantly affect Nd isotopic systematics or modify Sr isotopes.Such contaminated MORB character also revealed by other ophiolitic metavolcanics in the Central Eastern Desert, contrasting the N-MORB character of the Gerf ophiolite in the South Eastern Desert. Moreover, the present work suggests the increase of the degree of contamination of the ophiolitic metabasalts from south to north in the Central Eastern Desert.     

Seismic spectrograms of an anelastic layer with different source-receivers configurations

Thin and relatively thin anelastic layers (compared to the signal wavelength) generally represent hydrocarbon reservoirs, where the rock is a sandstone or a source rock saturated with brine, oil and gas. The study of the seismic response of these layers is important to detect the hydrocarbons on the basis of the reflection and transmission coefficients and the wave velocity and attenuation properties. Different seismic experiments (source-receiver configurations) can provide useful information to characterise its properties. In this work, we consider varying thicknesses and Q values of the layer and analyse the reflection and transmission coefficients. Moreover, we obtain spectrograms of surface seismic profiles and vertical and horizontal well profiles (VSP and HSP, respectively) to analyse their frequency content with offset due to variations of the attenuation properties of the layer. In addition, we compare the effects due to NMO stretching and intrinsic attenuation related to the low-frequency shadows (LFS) observed in real data after stacking, since LFS can have several causes. Ambiguity is present in this case, indicating that non-stretch NMO is required, otherwise an offset mute of the data may remove useful information regarding the intrinsic (physical) loss.     

Ensemble smoothing of land subsidence measurements for reservoir geomechanical characterization

Geomechanical models are often used to predict the impact on land surface of fluid withdrawal from deep reservoirs, as well as investigating measures for mitigation. The ability to accurately simulate surface displacements, however, is often impaired by limited information on the geomechanical parameters characterizing the geological formations of interest. In this study, we employ an ensemble smoother, a data assimilation algorithm, to provide improved estimates of reservoir parameters through assimilation of measurements of both horizontal and vertical surface displacement into geomechanical model results. The method leverages the demonstrated potential of remote sensing techniques developed in the last decade to provide accurate displacement data for large areas of the land surface. For evaluation purposes, the methodology is applied to the case of a disk‐shaped reservoir embedded in a homogeneous, isotropic, and linearly elastic half space, subject to a uniform change in fluid pressure. Multiple sources of uncertainty are investigated, including the radius, R, the thickness, h, and the depth, c, of the reservoir; the pore pressure change, Δp; porous medium's vertical uniaxial compressibility, c_M, and Poisson's ratio, ν, and the ratio, s, between the compressibilities of the medium during loading and unloading cycles. Results from all simulations show that the ensemble smoother has the capability to effectively reduce the uncertainty associated with those parameters to which the variability and the spatial distribution of land surface displacements are most sensitive, namely, R, c, c_M, and s. These analyses demonstrate that the estimation of these parameters values depends on the number of measurements assimilated and the error assigned to the measurement values. Copyright © 2014 John Wiley & Sons, Ltd.     

Simplified seismic fatigue evaluation for rigid steel connections

A simplified fatigue-life model is proposed for assessing the seismic inelastic rotational capacity of steel connections. First relations are developed for rigid steel connections under lateral loading. Next this is extended to account for the effects of the welded steel moment frame (WSMF) connections of the so-called pre-Northridge type. The seismic fatigue theory is validated against experimental results. The experiments were conducted under increasing ductility amplitudcs until the onset of fracture. Miner‘ rule was used to convert the test results to given an equivalent constant amplitude cyclic fatigue life. Satisfactory agreement is obtained when comparing the experimental observations with the theoretical predictions.     

Line-based modified iterated Hough transform for automatic registration of multi-source imagery

Image registration aims at combining imagery from multiple sensors to achieve higher accuracy and derive more information than that obtained from a single sensor. The enormous increase in the volume of remotely sensed data that is being acquired by an ever-growing number of earth observation satellites mandates the development of accurate, robust, and automated registration procedures. An effective automatic image registration has to deal with four issues: registration primitives, transformation function, similarity measure, and matching strategy. This paper introduces a new approach for automatic image registration using linear features as the registration primitives. Linear features have been chosen because they can be reliably extracted from imagery with significantly different geometric and radiometric properties. The modified iterated Hough transform (MIHT), which manipulates the registration primitives and similarity measure, is used as the matching strategy for automatically deriving an estimate of the parameters involved in the transformation function as well as the correspondence between conjugate primitives. The MIHT procedure follows an optimal sequence for parameter estimation that takes into account the contribution of linear features with different orientations at various locations within the imagery towards the estimation of the transformation parameters in question. Experimental results using real data proved the feasibility and robustness of the suggested approach.