Anthropogenic contamination of groundwater with nitrate in arid region: case study of southern Hodna (Algeria)

In southern Hodna, an arid region of Algeria, the extension of irrigated agriculture and overexploitation of groundwaters have led to a rupture with traditional land use by pastoralism. Due to the arid conditions, groundwaters are the only water resources for man and irrigation. This study focuses on nitrate contamination of groundwaters and its increase since the last 40 years, on the basis of two archives and recent analyses. The chemical facies of waters is sulphate-chloride, the risk of salinity is moderate to large. The sodicity risk, as estimated by the SAR, is low due to the large salinity and the sandy soil texture. Nitrate contamination of deep groundwater is a clear evidence when recent data are compared to the historical data. From 1996 to 2008, the proportion of samples with NO ₃ ^? concentrations larger than 50 mg/L increases from 24 to 61 %, the proportion of samples with NO ₃ ^? concentrations larger than 100 mg/L increases from 12 to 27 %. The study points at a general contamination of aquifers from anthropic origin: phreatic water was already contaminated when pastoralism was dominant (1967); as population increased and irrigation agriculture developed with large use of N-fertilizers and organic amendments, contamination extended spatially and vertically to deep aquifer. To remediate this contamination, fertilizers should be fractionated and limited to plant requirements on the basis of soil and plant analyses. Excess of irrigation should be avoided. Soil surface should be continuously covered to prevent nitrate leaching. Secondly, more efficient irrigation techniques, such as localized irrigation, should be implemented, with adaptation of crop systems and fertilizer inputs to arid conditions. Farmers should then be trained; simultaneously, local agronomic references that are presently lacking should be acquired.     

Estimation of anisotropic hydraulic conductivity using geophysical data in a coastal aquifer of Karnataka,India

Estimation of hydraulic parameters is essential to understand the interaction between groundwater flow and seawater intrusion. Though several studies have addressed hydraulic parameter estimation, based on pumping tests as well as geophysical methods, not many studies have addressed the problem with clayey formations being present. In this study, a methodology is proposed to estimate anisotropic hydraulic conductivity and porosity values for the coastal aquifer with unconsolidated formations. For this purpose, the one-dimensional resistivity of the aquifer and the groundwater conductivity data are used to estimate porosity at discrete points. The hydraulic conductivity values are estimated by its mutual dependence with porosity and petrophysical parameters. From these estimated values, the bilinear relationship between hydraulic conductivity and aquifer resistivity is established based on the clay content of the sampled formation. The methodology is applied on a coastal aquifer along with the coastal Karnataka, India, which has significant clayey formations embedded in unconsolidated rock. The estimation of hydraulic conductivity values from the established correlations has a correlation coefficient of 0.83 with pumping test data, indicating good reliability of the methodology. The established correlations also enable the estimation of horizontal hydraulic conductivity on two-dimensional resistivity sections, which was not addressed by earlier studies. The inventive approach of using the established bilinear correlations at one-dimensional to two-dimensional resistivity sections is verified by the comparison method. The horizontal hydraulic conductivity agrees with previous findings from inverse modelling. Additionally, this study provides critical insights into the estimation of vertical hydraulic conductivity and an equation is formulated which relates vertical hydraulic conductivity with horizontal. Based on the approach presented, the anisotropic hydraulic conductivity of any type aquifer with embedded clayey formations can be estimated. The anisotropic hydraulic conductivity has the potential to be used as an important input to the groundwater models.     

3-D seismic velocities calculated from lattice-preferred orientation and reflectivity of a lower crustal section: examples of the Val Sesia section (Ivrea zone, northern Italy)

To quantify the seismic properties of lower crustal rocks and to better constrain the origin of the lower crustal seismic reflectivity, we determined the complete 3-D seismic properties of a lower crustal section. Eight representative samples of the main lithologic and structural units outcropping in the Val Sesia (Ivrea zone) were studied in detail. The seismic velocities were calculated using the single crystal stiffness coefficients and the lattice preferred orientation (LPO) of each mineral in all samples. The 21 stiffness coefficients characterizing the elastic behaviour of each rock are determined. Mafic and ultramafic rocks such as pyroxenite and pyroxene-bearing gabbros display complex shear wave properties. These rocks are weakly birefringent (maximum 0.1 kms⁻¹) and it is difficult to find consistent relationships between the seismic properties and the rock structure. On the other hand, seismic properties of deformed felsic rocks are essentially controlled by mica. They display strong S -wave birefringence (0.3 km s⁻¹) and relatively high V _p anisotropy (7.6 per cent). Amphibole also strongly influences the rock birefringence patterns. For both kind of rocks, the foliation is highly birefringent and the fast polarized shear wave is systematically oriented parallel to the foliation. We show that the number of mineral phases in the rock strongly controls the anisotropy. The seismic anisotropy has a complex role in the P -wave reflectivity. Compared to the isotropic case, anisotropy enhances the reflection coefficient for about 60 per cent of the possible lithological interfaces. For 40 per cent of the interfaces, the reflection coefficient is much lower when one considers the medium as anisotropic.     

Rotational failure mechanisms for the face stability analysis of tunnels driven by a pressurized shield

The aim of this paper is to determine the collapse and blow‐out face pressures of a circular tunnel driven by a pressurized shield. The analysis is performed in the framework of the kinematical approach of the limit analysis theory. Two rotational failure mechanisms are proposed for the active and passive cases. These mechanisms have two significant advantages with respect to the available ones: (i) they take into account the entire circular tunnel face instead of an inscribed ellipse to this circular area, and (ii) they are more consistent with the rotational rigid‐block movement observed in the experimental tests. For both the active and passive cases, the three‐dimensional failure surface was generated ‘point by point’ instead of simple use of the existing standard geometric shapes such as cones or cylinders. This was achieved by employing a spatial discretization technique. The numerical results have shown that the present rotational mechanisms provide, in the case of frictional soils (with or without cohesion), a significant improvement with respect to the translational mechanisms. Finally, an extension of the proposed collapse mechanism to include a tension cut‐off in the classical Mohr–Coulomb failure criterion is presented and discussed. Copyright © 2010 John Wiley & Sons, Ltd.     

The stability of gypsum in marine sediments using the entire ODP/IODP porewater composition database

The stability of gypsum in marine sediments has been investigated through the calculation of its saturation index at the sediment in situ temperature and pressure, using the entire ODP/IODP porewater composition database (14416 samples recovered from sediments collected during 95 ODP and IODP Legs). Saturation is reached in sediment porewaters of 26 boreholes drilled at 23 different sites, during 12 ODP/IODP Legs. As ocean bottom seawater is largely undersaturated with respect to gypsum, the porewater Ca content or its SO₄ concentration, or both, must increase in order to reach equilibrium. At several sites equilibrium is reached either through the presence of evaporitic gypsum layers found in the sedimentary sequence, and/or through a salinity increase due to the presence of evaporitic brines with high concentrations of Ca and SO₄. Saturation can also be reached in porewaters of seawater-like salinity (≈ 35‰), provided sulfate reduction is limited. In this case, saturation is due to the alteration of volcanogenic material which releases large amounts of Ca to the porewaters, where the Ca concentration can reach 55 times its seawater value as for example at ODP Leg 134 site 833. At a few sites, saturation is reached in hydrothermal environments, or as a consequence of the alteration of the basaltic basement. In addition to the well known influence of brines on the formation of gypsum, these results indicate that the alteration of sediments rich in volcanogenic material is a major process leading to gypsum saturation in marine sediment porewaters. Therefore, the presence of gypsum in ancient and recent marine sediments should not be systematically interpreted as due to hypersaline waters, especially if volcanogenic material is present.     

Dolomitic concretions in the Eocene Sobrarbe delta (Spanish Pyrenees): Fluid circulation above a submarine slide scar infilling

Numerous dolomite concretions have been discovered in marls of the Eocene Sobrarbe deltaic complex as part of the Ainsa Basin (Spain). This paper presents the first analyses of the shapes, the spatial relationships, the mineralogical, chemical and isotopic compositions of these concretions.The concretions are located above a major fossil submarine slide scar. They are mainly perpendicular to the sedimentary layers. Four distinct shapes of concretions have been distinguished: horizontal flat, sub-vertical cylindrical or cylindrical-complex and stocky. Three main mineral phases are associated with most of the concretions: calcite, celestite and barite. Concretion shapes and mineral occurrences are organized vertically in the marls from bottom to top: (i) at the bottom, flat shapes with septarian cracks filled by calcite and celestite, (ii) in the middle and at the top, cylindrical and cylindrical-complex concretions associated with prismatic barite, calcite and celestite filling conduits related to bioturbations, and (iii) at the top, cylindrical and cylindrical-complex concretions associated with calcite and celestite filling conduits related to bioturbations, and stocky shape concretions.We postulate that concretions have formed by dolomite cementation of the surrounding marls during early diagenesis in the zone of methanogenesis. The high sedimentation rate of the infilling seems to be a factor controlling the mineralogical composition of the concretions. Brown calcite precipitated in voids and fractures of the concretions. Celestite precipitated during burial, completing the filling of voids and fractures. Barite precipitated before celestite, but its time of precipitation relative to brown calcite remains unknown.     

Upper-mantle flow beneath French Polynesia from shear wave splitting

Upper-mantle flow beneath the South Pacific is investigated by analysing shear wave splitting parameters at eight permanent long-period and broad-band seismic stations and 10 broad-band stations deployed in French Polynesia from 2001 to 2005 in the framework of the Polynesian Lithosphere and Upper Mantle Experiment (PLUME). Despite the small number of events and the rather poor backazimuthal coverage due to the geographical distribution of the natural seismicity, upper-mantle seismic anisotropy has been detected at all stations except at Tahiti where two permanent stations with 15 yr of data show an apparent isotropy. The median value of fast polarization azimuths (N67.5°W) is parallel to the present Pacific absolute plate motion direction in French Polynesia (APM: N67°W). This suggests that the observed SKS fast polarization directions result mainly from olivine crystal preferred orientations produced by deformation in the sublithospheric mantle due to viscous entrainment by the moving Pacific Plate and preserved in the lithosphere as the plate cools. However, analysis of individual measurements highlights variations of splitting parameters with event backazimuth that imply an actual upper-mantle structure more complex than a single anisotropic layer with horizontal fast axis. A forward approach shows that a two-layer structure of anisotropy beneath French Polynesia better explains the splitting observations than a single anisotropic layer. Second-order variations in the measurements may also indicate the presence of small-scale lateral heterogeneities. The influence of plumes or fracture zones within the studied area does not appear to dominate the large-scale anisotropy pattern but may explain these second-order splitting variations across the network.     

A Tertiary asthenospheric flow beneath the southern French Massif Central indicated by upper mantle seismic anisotropy and related to the west Mediterranean extension

Upper mantle flow beneath the French Massif Central is investigated using teleseismic shear wave splitting induced by seismic anisotropy. About 25 three-component stations (short period, intermediate and broadband) were installed during the period 1998-1999 in the southern Massif Central, from the Clermont Ferrand volcanic area to the Mediterranean Sea. Teleseismic shear waves (SKS, SKKS and PKS) were used to determine the splitting parameters: the fast polarization direction and the delay time. Delay times ranging between 0.7 and 1.5 s have been observed at most of the sites. The azimuths of the fast split shear waves trend homogeneously NW-SE in the southern Massif Central suggesting a homogeneous mantle flow beneath this area. The observed NW-SE direction differs from the N100°E Pyrenean anisotropy further south. It does not appear to be correlated to Hercynian structures nor to the present-day motion of the plate but is well correlated to the Tertiary extension direction. We propose that the opening of the western Mediterranean induced by the rotation of the Corsica-Sardinia lithospheric block and the roll-back to the SE of the Tethys slab may have generated a large asthenospheric mantle flow beneath the southern Massif Central and a deflection of the up going plume centered beneath the northern Massif Central toward the SE.