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.     

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.     

Cr-spinel/olivine and Cr-spinel/liquid nickel partition coefficients from natural samples

The Nernst partition coefficient of nickel (D^Ni) between Cr-spinel and silicate melt in natural systems has been investigated using mid-ocean ridge basalts (MORB) and other volcanic rocks. The Cr-spinel/olivine D^Ni values in volcanic rocks are between 1.2 and 0.3, indicating that the Cr-spinel/liquid D^Ni values vary from slightly higher to significantly lower than the olivine/liquid D^Ni values in natural systems. The Cr-spinel/liquid D^Ni values from the MORB samples vary between 6 and 11, slightly higher than those from the S-bearing experiments of Satari et al. [Satari P., Brenan J. M., Horn I. and McDonough W. F. (2002) Experimental constraints on the sulfide- and chromite-silicate melt partitioning behavior of rhenium and platinum-group elements. Economic Geology97, 385-398]. The results of the MORB samples and the experiments of Satari et al. (2002) indicate a negative correlation between the Cr-spinel/liquid D^Ni and the X_Cr values in Cr-spinels (Cr cation number on the basis of 3 total cations in the spinel structure). Variations of Cr-spinel/liquid D^Ni values with Cr-spinel compositions can be estimated from an empirical equation based on the results of the MORB samples and the experiments by Satari et al. (2002). The choice of Cr-spinel/liquid D^Ni = 10 for numerical modeling by Righter et al. [Righter K., Leeman W. P. and Hervig R. L. (2006) Partitioning of Ni, Co, and V between spinel-structured oxides and silicate melts: importance of spinel composition. Chemical Geology227, 1-25] is reasonable for basaltic systems. For picritic and komatiitic systems a lower value of ∼5 is more appropriate.     

Reaction mechanism of uranyl in the presence of zero-valent iron nanoparticles

The current study provides an investigation of abiotic reduction of an oversaturated uranyl solution driven by iron nanoparticle oxidation. The reactivity of nano-scale zero-valent iron (ZVI) under mildly oxic conditions (1.2% O₂ and 0.0017% CO₂) was studied in 1000 ppm uranyl solution in the pH range 3-7, at reaction times from 10 min to 4 h. Reductive precipitation of UO₂ was observed as the main process responsible for the removal of uranium from solution with the kinetics of reaction becoming increasingly favourable at higher pH. Despite working with an oversaturated uranium solution, the precipitation of UO₂ occurred in preference to precipitation of UO₃·2H₂O (metaschoepite) at reaction times between 1 and 4 h and for uranyl solutions initially set up at pH ?5. Characterisation of both solid and solution phases was performed using X-ray photoelectron spectroscopy (XPS), focused ion beam (FIB) imaging, X-ray diffraction (XRD) and inductively coupled plasma atomic emission spectroscopy (ICP-AES).     

Biotite weathering and nutrient uptake by ectomycorrhizal fungus, Suillus tomentosus, in liquid-culture experiments

Ectomycorrhiza-forming fungi (EMF) alter the nutrient-acquisition capabilities of vascular plants, and may play an important role in mineral weathering and the partitioning of products of weathering in soils under nutrient-limited conditions. In this study, we isolated the weathering function of Suillus tomentosus in liquid-cultures with biotite micas incubated at room temperature. We hypothesized that the fungus would accelerate weathering by hyphal attachment to biotite surfaces and transmission of nutrient cations via direct exchange into the fungal biomass. We combined a mass-balance approach with scanning electron microscopy (SEM) and atomic force microscopy (AFM) to estimate weathering rates and study dissolution features on biotite surfaces. Weathering of biotite flakes was about 2-3 orders of magnitude faster in shaken liquid-cultures with fungus compared to shaken controls without fungus, but with added inorganic acids. Adding fungus in nonshaken cultures caused a higher dissolution rate than in inorganic pH controls without fungus, but it was not significantly faster than organic pH controls without fungus. The K⁺, Mg²⁺ and Fe²⁺ from biotite were preferentially partitioned into fungal biomass in the shaken cultures, while in the nonshaken cultures, K⁺ and Mg²⁺ was lost from biomass and Fe²⁺ bioaccumulated much less. Fungal hyphae attached to biotite surfaces, but no significant surface changes were detected by SEM. When cultures were shaken, the AFM images of basal planes appeared to be rougher and had abundant dissolution channels, but such channel development was minor in nonshaken conditions. Even under shaken conditions the channels only accounted for only 1/100 of the total dissolution rate of 2.7 × 10⁻¹⁰ mol of biotite m⁻² s⁻¹. The results suggest that fungal weathering predominantly occurred not by attachment and direct transfer of nutrients via hyphae, but because of the acidification of the bulk liquid by organic acids, fungal respiration (CO₂), and complexation of cations which accelerated dissolution of biotite. Results further suggest that both carbohydrate source (abundant here) and a host with which nutrients are exchanged (missing here) may be required for EMF to exert an important weathering effect in soils. Unsaturated conditions and physical dispersal of nutrient-rich minerals in soils may also confer a benefit for hyphal growth and attachment, and promote the attachment-mediated weathering which has been observed elsewhere on soil mineral surfaces.     

Isotopic studies of leaf water. Part 1: A physically based two-dimensional model for pine needles

Understanding of isotopic variations in leaf water is important for reconstruction of paleoclimate and assessment of global biochemical processes. We report here a study of isotopic distributions within a single needle of two pine species, Pinus resinosa Ait and Pinus strobes L., with the objective of understanding how isotopic compositions of leaf water are controlled by environmental and physiological variables. A 2D model was developed to simulate along-leaf isotopic variations and bulk leaf water isotopic compositions. In addition to variables common to all leaf water isotopic models, this 2D model also takes into account the specific geometry and dimensions of pine needles and the isotopic transport in xylem and mesophyll. The model can successfully simulate oxygen isotopic variations along a single needle and averaged over a leaf (bulk leaf water). The simulations suggest that isotopic composition of the bulk leaf water does not always depend only upon the average transpiration rate, which in turn raises questions about using leaf water isotopic values to estimate transpiration rates. An unsuccessful attempt to simulate along-needle hydrogen isotopic variations suggests that certain unknown biological process(es) may not have been incorporated into our 2D model, and if so, it calls for a reevaluation of all other models for hydrogen isotopic simulations of leaf water since they too lack these processes.Existing leaf water isotopic models are reviewed in this work. In particular, we evaluate the most frequently used model, the stomatal boundary layer model (also referred to as the Craig-Gordon model). We point out that discrepancy between the boundary layer model and the measured bulk leaf water seems to depend upon relative humidity. Using our 2D model, we show that this humidity dependency is a result of an interplay between environmental and physiological conditions: if the transpiration rate of plant leaves decreases with increasing relative humidity, our 2D model can reproduce the pattern of isotopic discrepancy between boundary layer model predictions and observations, enabling us to understand better the reason behind this discrepancy.