Assessing groundwater storage changes in the Nubian aquifer using GRACE data

Gravity Recovery and Climate Experiment (GRACE) level two (L2) data is used in estimating the groundwater storage changes (GWSC) in the Nubian Sandstone Aquifer System (NSAS). This set of data consists of spherical harmonics coefficients with specific degree and order. The GRACE data is de-correlated using a sixth degree polynomial in order to reduce the effect of the noise error resulting from the correlation between the spherical harmonics coefficients with the same degree parity. The GRACE estimates of GWSC are smoothed using Gaussian filter with half width of 1000 km. This half width is chosen in order to maximize the correlation between the GRACE estimates of GWSC and previous modeling results of the NSAS. The loss in groundwater storage occurring in each of the four countries sharing the NSAS is calculated to assess the sustainability of using the NSAS as a water resource in each country. The overarching finding in this study is that NSAS is losing its groundwater storage at a very high rate. Also, it is found that Egypt is the fastest in losing its groundwater storage from the NSAS. This loss of groundwater storage in Egypt may not necessarily be resulting from in-country extractions because of the trans-boundary nature of this aquifer. The GRACE-based estimates are found to be close to available data and previous modeling results of the NSAS.     

Automatic identification of fake patterns caused by short-width wavelets in seismic data

Despite the popularity of using the Haar wavelet filter in many applications, it sometimes introduces fake patterns into the multi resolution analysis (MRA) of seismic data. In this work, we compared different wavelet filters to demonstrate that these patterns are fake and not part of the original waveforms and to show that they are a result of using the Haar wavelet filter as a short-width wavelet. To achieve this, many seismic waveforms from two different sources: the Egyptian National Seismic Network (ENSN) and the High Sensitivity Seismograph Network Japan (Hi-net) are used with different wavelet filters. We propose an algorithm based on an autoregressive (AR) model to detect these patterns automatically and fully.     

Tectono-sedimentary evolution of the Miocene Oued Dayr Formation (Ghomaride Complex,Internal Domain of the Maghrebian Rif Belt,Morocco)

This work deals with sedimentological, petrographic, and structural analyses of a middle Miocene late-orogenic sedimentary cycle, denoted Oued Dayr Formation, recognized in the Rifian sector of the Maghrebian Chain (Morocco). The analyzed Formation (75 m thick) starts with 15–20 m of light colored polymict conglomerates, with minor sandstone beds, lying on the Paleozoic basement and Mesozoic cover of the Ghomaride Nappe. Facies analysis indicates a fining-upward deposition in a marine environment characterized by increasing deepening, reflecting a subsidence rate that exceeds sedimentary supply. Petrographic analysis points out that sandstones are represented by litharenites originated by erosion of recycled orogen. The conglomerates pebbles and cobbles consist of Alpine low- to high-grade metamorphic rocks as metarenites, phyllites, mylonitic quartzites, micaschists, augen gneisses deriving from the exhumed deep metamorphic basement, the overlying metasedimentary of the Sebtide Nappes and of sedimentary rocks as sandstones, jaspes, limestones, and shales deriving from the Ghomaride Nappes and their sedimentary cover. Data reveal mixed provenance indicating that the Oued Dayr Formation was fed by the Internal Nappes stack of the Maghrebian Chain. Structural analysis shows that the Oued Dayr Formation accumulated in a Thrust-Top basin, during an early extension (D₀ phase), recorded by synsedimentary normal faults within middle Langhian deposits on the rear of the Internal Nappes stack. Subsequent ductile and brittle compressional (D₁, D₂, D₃) and extensional (D₄) deformation phases occurred during and/or after the stacking, exhumation, and early unroofing of Sebtide Complex coeval with the opening of the western Mediterranean back-arc basins since middle Miocene time.     

The influence of physical and biological processes on the ichthyoplankton communities in the Gulf of Sirte (Southern Mediterranean Sea)

This paper reports on the first study on the ichthyoplankton community structure in the Gulf of Sirte and its spatial distribution in relation to environmental conditions and zooplankton abundance. The results make an important contribution to the future management of fisheries in this unexplored, but much exploited, area. Biological samples were collected in July 2008 using a Bongo40 net. In total, 1914 larvae were found and 1652 of these were identified. In particular, bathypelagic taxa were the most abundant, followed by demersal, mesopelagic, pelagic and epipelagic taxa. The ichthyoplankton community had a patchwork distribution influenced by oligotrophic conditions, the bottom depth and oceanographical features. The results suggest that environmental forcings were able to transport the ichthyoplankton to productive areas. Indeed, maximum fish egg densities were found in coastal stations in correspondence with the Atlantic Tunisian Current inflow, whereas larvae were mainly concentrated in the east side of the Gulf, probably as a result of advection by the anticyclonic circulation. Additionally, the distribution patterns of the total larvae density and the different assemblages were well matched with the abundance of the zooplankton, probably determining final larval survival, growth and recruitment.     

Diversity patterns of free‐living marine nematodes in the southwest continental shelf off Bay of Bengal and their link to abiotic variables

Marine sediments in continental shelf ecosystems harbor a rich biodiversity of benthic communities. In this study, the spatial and temporal diversity and community assemblages of free‐living marine nematodes were studied by sampling at six depths and over 3 years from the southwest continental shelf off Bay of Bengal, one of the least explored tropical shelf ecosystems. The dominant marine nematode species were related with abiotic variables as part of this study. The effects of sediment granulometry generally decreased with increasing depth and the highest nematode density and species diversity were recorded on coarse sand (shallower depths). Multivariate analysis of the nematode community data showed that community structure differed significantly among depths as well as among years. Statistical analyses showed significant correlations between the nematode community and abiotic variables. Sediment texture, organic matter, water pressure and depth profile were crucial factors for determining diversity, vertical profile and feeding types of the nematode community. Other environmental factors, including anthropogenic pressure, did not have an effect on nematode diversity except for the presence of some tolerant species (Metachromadora spp., Sabatieria spp. and Siplophorella sp.). This study represents a baseline of knowledge of free‐living marine nematode communities that can be used in the future to compare nematode assemblages from temperate shelf ecosystems.     

Steady critical discharge rates from vertical and horizontal wells in fresh–saline aquifers with sharp interfaces

In this study, the critical (or maximum) discharge rates before saline water enters a well were determined for vertical and horizontal wells in a freshwater aquifer which is separated from a static saline aquifer by a sharp interface. Flow around the well was solved by integration of a point sink solution along the well axis, and both the critical discharge rate and critical interface rise were determined through a comparison of the heads and vertical gradients at the saline–fresh water interface. The rates were determined for vertical and horizontal wells with various lengths and depths for different aquifer salinities. Results were generalized by drawing dimensionless type curves. The results showed that the dimensionless total critical discharge rates are higher for the longer horizontal wells and longer vertical wells with a certain bottom depth, and they almost linearly decrease with well depth at rates of 0.7–0.9. For the dimensionless well length of 0.2, the dimensionless total discharge rate of a horizontal well is about 0.1 more than that of a vertical well with the same length and well-top depth. Also, the critical discharge rates per unit length of well are inversely proportional to well length and remarkably higher for shallower wells. Additionally, the critical pumping rate is proportional to the salinity difference of the aquifers. These results were confirmed by comparison to existing solutions for vertical wells with dimensionless lengths of 0.2, 0.5 and 0.6, and for critical interface rises in the range of 0.75–0.9.     

Improved regional groundwater flow modeling using drainage features: a case study of the central northern karst aquifer system of Puerto Rico (USA)

In northern Puerto Rico (USA), subsurface conduit networks with unknown characteristics, and surface features such as springs, rivers, lagoons and wetlands, drain the coastal karst aquifers. In this study, drain lines connecting sinkholes and springs are used to improve the developed regional model by simulating the drainage effects of conduit networks. Implemented in an equivalent porous media (EPM) approach, the model with drains is able to roughly reproduce the spring discharge hydrographs in response to rainfall. Hydraulic conductivities are found to be scale dependent and significantly increase with higher test radius, indicating scale dependency of the EPM approach. Similar to other karst regions in the world, hydraulic gradients are steeper where the transmissivity is lower approaching the coastline. This study enhances current understanding of the complex flow patterns in karst aquifers and suggests that using a drainage feature improves modeling results where available data on conduit characteristics are minimal.     

Hydrogeological investigation of an oasis-system aquifer in arid southeastern Morocco by development of a groundwater flow model

Groundwater of the Tafilalet oasis system (TOS) is an important water resource in the lower Ziz and Rheris valleys of arid southeastern Morocco. The unconfined aquifer is exploited for domestic consumption and irrigation. A groundwater flow model was developed to assess the impact of climatic variations and development, including the construction of hydraulic structures, on the hydrodynamic behavior of the aquifer. Numerical simulations were performed by implementing a spatial database within a geographic information system and using the Arc Hydro Groundwater tool with the code MODFLOW-2000. The results of steady-state and transient simulations between 1960 and 2011 show that the water table is at equilibrium between recharge, which is mainly by surface-water infiltration, and discharge by evapotranspiration. After the commissioning of the Hassan Addakhil dam in 1971, hydraulic heads became more sensitive to annual variations than to seasonal variations. Heads are also influenced by recurrent droughts and the highest water-level changes are recorded in irrigated areas. The model provides a way of managing groundwater resources in the TOS. It can be used as a tool to predict the impact of different management plans for the protection of groundwater against overexploitation and deterioration of water quality.     

Applying disk-based discontinuous deformation analysis (DDA) to simulate Donghekou landslide triggered by the Wenchuan earthquake

In this paper, for the first time, disk-based discontinuous deformation analysis (DDA) is applied to simulate a real landslide triggered by an earthquake. For this purpose, the kinematic behaviour of the Donghekou landslide triggered by the Wenchuan earthquake is simulated and the results obtained using disk-based DDA are compared with actual data. The comparisons show that there is a good agreement between the results obtained using disk-based DDA and observed data. The simulation results provided an understanding of the failure behaviour and temporal evolution of the landslide. This study shows that disk-based DDA is a practical numerical tool that can be used to simulate the post-failure behaviour of landslides triggered by an earthquake.