A transdisciplinary analysis of water problems in the mountainous karst areas of Morocco

Water management and engineering in the karstic High Atlas of Morocco are difficult tasks under the prevailing geological, hydrogeological, geomorphological, vegetational and climatic conditions. It is important to be able to understand and predict the characteristics and availability of water for future water planning in the region under changing climatic and agricultural conditions. An interdisciplinary analysis of problems and adequate hydrological modelling tools developed by geologists, hydrologists and biologists are necessary. The karst areas of the High Atlas Mountains are characterised by impermeable triassic basalt underlying substantial subsurface reservoirs with high potential discharge rates. The karst groundwater aquifers are extensive but largely unknown in dimension, probably with a hierarchical network of groundwater flow paths. It is estimated that approximately 70% of the surface water is directly lost to groundwater. Steep landslide- and debris flow prone slopes exist next to coarse-grained, highly porous river beds. Infrequent, high intensity rainfall or snowmelt causes a particularly high flood risk to these karst areas. In addition, agriculture and land use changes have degraded the karst areas. The most important driving forces for degradation include permanent overgrazing even during droughts and the use of firewood by a continually growing population. Large scale degradation of vegetation has occurred in the oro-mediterranean (mountainous Mediterranean) zone, between 2600 and 3400 m which coincides with the most important zone for karstic groundwater creation. The combination of high amounts of groundwater flow and rapid surface flow due to sparse vegetation has increased the problems of flood flow.     

Shear strength characterization of municipal solid waste at the Suzhou landfill, China

The current practice of slope stability analysis for a municipal solid waste (MSW) landfill usually overlooks the dependence of waste properties on the fill age or embedment depth. Changes in shear strength of MSW as a function of fill age were investigated by performing field and laboratory studies on the Suzhou landfill in China. The field study included sampling from five boreholes advanced to the bottom of the landfill, cone penetration tests and monitoring of pore fluid pressures. Twenty-six borehole samples representative of different fill ages (0 to 13 years) were used to perform drained triaxial compression tests. The field and laboratory study showed that the waste body in the landfill can be sub-divided into several strata corresponding to different ranges of fill age. Each of the waste strata has individual composition and shear strength characteristics. The triaxial test results showed that the MSW samples exhibited a strain-hardening and contractive behavior. As the fill age of the waste increased from 1.7 years to 11 years, the cohesion mobilized at a strain level of 10% was found to decrease from 23.3 kPa to 0 kPa, and the mobilized friction angle at the same strain level increasing from 9.9° to 26°. For a confinement stress level greater than 50 kPa, the shear strength of the recently-placed MSW seemed to be lower than that of the older MSW. This behavior was consistent with the cone penetration test results. The field measurement of pore pressures revealed a perched leachate mound above an intermediate cover of soils and a substantial leachate mound near the bottom of the landfill. The measurements of shear strength properties and pore pressures were utilized to assess the slope stability of the Suzhou landfill.     

Assessing soil erosion in Mediterranean karst landscapes of Lebanon using remote sensing and GIS

Maps showing the potential for soil erosion at 1:100,000 scale are produced in a study area within Lebanon that can be used for evaluating erosion of Mediterranean karstic terrain with two different sets of impact factors built into an erosion model. The first set of factors is: soil erodibility, morphology, land cover/use and rainfall erosivity. The second is obtained by the first adding a fifth factor, rock infiltration. High infiltration can reflect high recharge, therefore decreasing the potential of surface runoff and hence the quantity of transported materials. Infiltration is derived as a function of lithology, lineament density, karstification and drainage density, all of which can be easily extracted from satellite imagery. The influence of these factors is assessed by a weight/rate approach sharing similarities between quantitative and qualitative methods and depending on pair-wise comparison matrix.The main outcome was the production of factorial maps and erosion susceptibility maps (scale 1:100,000). Spatial and attribute comparison of erosion maps indicates that the model that includes a measure of rock infiltration better represents erosion potential. Field investigation of rills and gullies shows 87.5% precision of the model with rock infiltration. This is 17.5% greater than the precision of the model without rock infiltration. These results indicate the necessity and importance of integrating information on infiltration of rock outcrops to assess soil erosion in Mediterranean karst landscapes.     

Oxygen isotopic compositions of Allende Type C CAIs: Evidence for isotopic exchange during nebular melting and asteroidal metamorphism

In situ oxygen isotopic measurements of primary and secondary minerals in Type C CAIs from the Allende CV3 chondrite reveal that the pattern of relative enrichments and depletions of ¹⁶O in the primary minerals within each individual CAI are similar to the patterns observed in Types A and B CAIs from the same meteorite. Spinel is consistently the most ¹⁶O-rich (Δ¹⁷O = −25‰ to −15‰), followed by Al,Ti-dioside (Δ¹⁷O = −20‰ to −5‰) and anorthite (Δ¹⁷O = −15‰ to 0‰). Melilite is the most ¹⁶O-depleted primary mineral (Δ¹⁷O = −5‰ to −3‰). We conclude that the original melting event that formed Type C CAIs occurred in a ¹⁶O-rich (Δ¹⁷O  −20‰) nebular gas and they subsequently experienced oxygen isotopic exchange in a ¹⁶O-poor reservoir. At least three of these (ABC, TS26F1 and 93) experienced remelting at the time and place where chondrules were forming, trapping and partially assimilating ¹⁶O-poor chondrule fragments. The observation that the pyroxene is ¹⁶O-rich relative to the feldspar, even though the feldspar preceded it in the igneous crystallization sequence, disproves the class of CAI isotopic exchange models in which partial melting of a ¹⁶O-rich solid in a ¹⁶O-poor gas is followed by slow crystallization in that gas. For the typical (not associated with chondrule materials) Type C CAIs as well for as the Types A and B CAIs, the exchange that produced internal isotopic heterogeneity within each CAI must have occurred largely in the solid state. The secondary phases grossular, monticellite and forsterite commonly have similar oxygen isotopic compositions to the melilite and anorthite they replace, but in one case (CAI 160) grossular is ¹⁶O-enriched (Δ¹⁷O = −10‰ to −6‰) relative to melilite (Δ¹⁷O = −5‰ to −3‰), meaning that the melilite and anorthite must have exchanged its oxygen subsequent to secondary alteration. This isotopic exchange in melilite and anorthite likely occurred on the CV parent asteroid, possibly during fluid-assisted thermal metamorphism.     

Iron isotopes constrain biologic and abiologic processes in banded iron formation genesis

The voluminous 2.5 Ga banded iron formations (BIFs) from the Hamersley Basin (Australia) and Transvaal Craton (South Africa) record an extensive period of Fe redox cycling. The major Fe-bearing minerals in the Hamersley-Transvaal BIFs, magnetite and siderite, did not form in Fe isotope equilibrium, but instead reflect distinct formation pathways. The near-zero average δ⁵⁶Fe values for magnetite record a strong inheritance from Fe³⁺ oxide/hydroxide precursors that formed in the upper water column through complete or near-complete oxidation. Transformation of the Fe³⁺ oxide/hydroxide precursors to magnetite occurred through several diagenetic processes that produced a range of δ⁵⁶Fe values: (1) addition of marine hydrothermal , (2) complete reduction by bacterial dissimilatory iron reduction (DIR), and (3) interaction with excess that had low δ⁵⁶Fe values and was produced by DIR. Most siderite has slightly negative δ⁵⁶Fe values of ∼ −0.5‰ that indicate equilibrium with Late Archean seawater, although some very negative δ⁵⁶Fe values may record DIR. Support for an important role of DIR in siderite formation in BIFs comes from previously published C isotope data on siderite, which may be explained as a mixture of C from bacterial and seawater sources.Several factors likely contributed to the important role that DIR played in BIF formation, including high rates of ferric oxide/hydroxide formation in the upper water column, delivery of organic carbon produced by photosynthesis, and low clastic input. We infer that DIR-driven Fe redox cycling was much more important at this time than in modern marine systems. The low pyrite contents of magnetite- and siderite-facies BIFs suggests that bacterial sulfate reduction was minor, at least in the environments of BIF formation, and the absence of sulfide was important in preserving magnetite and siderite in the BIFs, minerals that are poorly preserved in the modern marine record. The paucity of negative δ⁵⁶Fe values in older (Early Archean) and younger (Early Proterozoic) BIFs suggests that the extensive 2.5 Ga Hamersley-Transvaal BIFs may record a period of maximum expansion of DIR in Earth’s history.     

Surface complexation of U(VI) on goethite (α-FeOOH)

Sorption of U(VI) to goethite is a fundamental control on the mobility of uranium in soil and groundwater. Here, we investigated the sorption of U on goethite using EXAFS spectroscopy, batch sorption experiments and DFT calculations of the energetics and structures of possible surface complexes. Based on EXAFS spectra, it has previously been proposed that U(VI), as the uranyl cation , sorbs to Fe oxide hydroxide phases by forming a bidentate edge-sharing (E2) surface complex, >Fe(OH)₂UO₂(H₂O)_n. Here, we argue that this complex alone cannot account for the sorption capacity of goethite (α-FeOOH). Moreover, we show that all of the EXAFS signal attributed to the E2 complex can be accounted for by multiple scattering. We propose that the dominant surface complex in CO₂-free systems is a bidentate corner-sharing (C2) complex, (>FeOH)₂UO₂(H₂O)₃ which can form on the dominant {101} surface. However, in the presence of CO₂, we find an enhancement of UO₂ sorption at low pH and attribute this to a (>FeO)CO₂UO₂ ternary complex. With increasing pH, U(VI) desorbs by the formation of aqueous carbonate and hydroxyl complexes. However, this desorption is preceded by the formation of a second ternary surface complex (>FeOH)₂UO₂CO₃. The three proposed surface complexes, (>FeOH)₂UO₂(H₂O)₃, >FeOCO₂UO₂, and (>FeOH)₂UO₂CO₃ are consistent with EXAFS spectra. Using these complexes, we developed a surface complexation model for U on goethite with a 1-pK model for surface protonation, an extended Stern model for surface electrostatics and inclusion of all known UO₂-OH-CO₃ aqueous complexes in the current thermodynamic database. The model gives an excellent fit to our sorption experiments done in both ambient and reduced CO₂ environments at surface loadings of 0.02-2.0 wt% U.     

Seasonal patterns of carbon chemistry and isotopes in tufa depositing groundwaters of southwestern Japan

Annually laminated carbonates, known as tufas, commonly develop in limestone areas and typically record seasonal patterns of oxygen- and carbon-isotope compositions. δ¹⁸O values are principally controlled by seasonal changes of water temperature, whereas δ¹³C values are the result of complex reactions among the gaseous, liquid, and solid sources of carbon in the system. We examined the processes that cause the seasonal patterns of δ¹³C in groundwater systems at three tufa-depositing sites in southwestern Japan by applying model calculations to geochemical data. Underground inorganic carbon species are exchanged with gaseous CO₂, which is mainly introduced to the underground hydrological system by natural atmospheric ventilation and by diffusion of soil air. These processes control the seasonal pattern of δ¹³C, which is low in summer and high in winter. Among the three sites we investigated, we identified two extreme cases of the degree of carbon exchange between liquid and gaseous phases. For the case with high radiocarbon composition (Δ¹⁴C) and low pCO₂, there was substantial carbon exchange because of a large contribution of atmospheric CO₂ and a small water mass. For the other extreme case, which was characterized by low Δ¹⁴C and high pCO₂, the contribution of atmospheric CO₂ was small and the water mass was relatively large. Our results suggest that at two of the three sites water residence time within the soil profile was longer than 1 year. Our results also suggested a short residence time (less than 1 year) of water in the soil profile at the site with the smallest water mass, which is consistent with large seasonal amplitude of the springwater temperature variations. The Δ¹⁴C value of tufas is closely related to the hydrological conditions in which they are deposited. If the initial Δ¹⁴C value of a tufa-depositing system is stable, ¹⁴C-chronology can be used to date paleo-tufas.     

Temperature effects on non-bridging oxygen and aluminum coordination number in calcium aluminosilicate glasses and melts

Configurational changes with temperature are important for the thermodynamic and transport properties of most aluminosilicate melts, but in general are not well understood. Here, we present high-resolution ²⁷Al and ¹⁷O NMR data on several calcium aluminosilicate glasses prepared with varying quench rates and thus with fictive temperatures that span ranges up to about 200 K. In all compositions the content of five-coordinated aluminum increases with fictive temperature, in agreement with recent high temperature NMR data on melts. In a glass of CaAl₂Si₂O₈ (“anorthite”) composition, the content of non-bridging oxygens also increases with temperature; however this effect was not observed in a sample with a much higher CaO/Al₂O₃ ratio. We present a consistent notation for reactions among structural species in these systems that clarify why in some cases, high-coordinated network cations may appear on the same side of the reaction, while in others they occur on the opposite sides: the key difference is in accounting for all coordination changes for oxygens. Mixing of non-bridging oxygens and of high-coordinated aluminum make significant contributions to the overall configurational entropy and heat capacity of the melts, as does the mixing of various bridging oxygens and of tetrahedral network cations. Other, less well known, types of increase in disorder with temperature may be important as well.     

Kinetics of Cd sorption, desorption and fixation by calcite: A long-term radiotracer study

Time-dependent sorption and desorption of Cd on calcite was studied over 210 days utilizing ¹⁰⁹Cd as a tracer to distinguish between ‘labile’ and ‘non-labile’ forms of sorbed Cd. Stabilizing the calcite suspensions for 12 months under atmospheric P_CO2 and controlled temperature was necessary to reliably follow Cd dynamics following initial sorption. Results revealed time-dependant Cd sorption and marked desorption hysteresis by calcite under environmentally relevant conditions. Data obtained were fitted to a first-order kinetic model and a concentric shell diffusion model. Both models described the progressive transfer of Cd²⁺ to a less reactive form within calcite and subsequent desorption of Cd subject to different initial contact times. The kinetic model provided a better fit to the combined sorption and desorption data (R² = 0.992). It differentiates between two ‘pools’ of sorbed Cd²⁺ on calcite, ‘labile’ and ‘non-labile’, in which labile sorbed Cd is in immediate equilibrium with the free Cd²⁺ ion activity in solution whereas non-labile Cd is kinetically restricted. For the diffusion model (R² = 0.959), the rate constants describing Cd dynamics in calcite produced a half-life for Cd desorption of ∼175 d, for release to a ‘zero-sink’ solution. Results from this study allow comment on the likely mechanisms occurring at the calcite surface following long-term Cd sorption.     

The age of the martian meteorite Northwest Africa 1195 and the differentiation history of the shergottites

Samarium-neodymium isotopic analyses of unleached and acid-leached mineral fractions from the recently identified olivine-bearing shergottite Northwest Africa 1195 yield a crystallization age of 347 ± 13 Ma and an value of +40.1 ± 0.9. Maskelynite fractions do not lie on the Sm-Nd isochron and appear to contain a martian surface component with low ¹⁴⁷Sm/¹⁴⁴Nd and ¹⁴³Nd/¹⁴⁴Nd ratios that was added during shock. The Rb-Sr system is disturbed and does not yield an isochron. Terrestrial Sr appears to have affected all of the mineral fractions, although a maximum initial ⁸⁷Sr/⁸⁶Sr ratio of 0.7016 is estimated by passing a 347 Ma reference line through the maskelynite fraction that is least affected by contamination. The high initial value and the low initial ⁸⁷Sr/⁸⁶Sr ratio, combined with the geologically young crystallization age, indicate that Northwest Africa 1195 is derived from a source region characterized by a long-term incompatible-element depletion.The age and initial Sr and Nd isotopic compositions of Northwest Africa 1195 are very similar to those of Queen Alexandra Range 94201, indicating these samples were derived from source regions with similar Sr-Nd isotopic systematics. These similarities suggest that these two meteorites share a close petrogenetic relationship and might have been erupted from a common volcano. The meteorites Yamato 980459, Dar al Gani 476, Sayh al Uhaymir 005/008, and Dhofar 019 also have relatively old ages between 474 and 575 Ma and trace element and/or isotopic systematics that are indicative of derivation from incompatible-element-depleted sources. This suggests that the oldest group of meteorites is more closely related to one another than they are to the younger meteorites that are derived from less incompatible-element-depleted sources. Closed-system fractional crystallization of this suite of meteorites is modeled with the MELTS algorithm using the bulk composition of Yamato 980459 as a parent. These models reproduce many of the major element and mineralogical variations observed in the suite. In addition, the rare earth element systematics of these meteorites are reproduced by fractional crystallization using the proportions of phases and extents of crystallization that are calculated by MELTS. Other shergottites that demonstrate enrichments in incompatible-elements and have evolved Sr and Nd isotopic systematics have some geochemical systematics that are similar to those observed in the depleted group. Most notably, although they exhibit a very limited range of incompatible trace element and isotopic compositions, they have highly variable major element compositions. This is also consistent with evolution from a common mantle source region by variable amounts of fractional crystallization. If this scenario is correct, it suggests that the combined effects of source composition and fractional crystallization are likely to account for the major element, trace element, and isotopic diversity of all shergottites.