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.     

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.     

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.     

Phthalic acid complexation and the dissolution of forsteritic glass studied via in situ FTIR and X-ray scattering

Multiple Internal Reflection Fourier Transform Infra-Red (MIR-FTIR) spectroscopy was developed and used for in situ flow-through experiments designed to study the process of organic acid promoted silicate dissolution. In tandem with the FTIR analysis, ex situ X-ray scattering was used to perform detailed analyses of the changes in the surface structure and chemistry resulting from the dissolution process. Phthalic acid and forsteritic glass that had been Chemically Vapour Deposited (CVD) onto an internal reflection element were used as reactants, and the MIR-FTIR results showed that phthalic acid may promote dissolution by directly binding to exposed Mg metal ion centers on the solid surface. Integrated infrared absorption intensity as a function of time shows that phthalic acid attachment apparently follows a t^1/2 dependence, indicating that attachment is a diffusive process. The diffusion coefficient of phthalic acid was estimated to be approximately 7 × 10⁻⁶ cm² s⁻¹ in the solution near the interface with the glass. Shifts in the infrared absorption structure of the phthalate complexed with the surface compared to the solute species indicate that phthalate forms a seven-membered ring chelate complex. This bidentate complex efficiently depletes Mg from the glass surface, such that after reaction as much as 95% of the Mg may be removed. Surface depletion in Mg causes adsorbate density to fall after an initial attachment stage for the organic ligand. In addition, the infrared analysis shows that silica in the near surface polymerizes after Mg removal, presumably to maintain charge balance. X-ray reflectivity shows that the dissolution rate of forsteritic glass at pH 4 based on Mg removal in such flow-through experiments was equal to 4 × 10⁻¹² mol cm⁻² s⁻¹ (geometric surface area normalized). Reflectivity also shows how the surface mass density decreases during reaction from 2.64 g cm⁻³ to 2.2 g cm⁻³, consistent with preferential loss of Mg from the surface. Auxiliary batch experiments with forsteritic glass films deposited onto soda glass were also completed to add further constraints to the mechanism of reaction. By combining reflectivity with diffuse scatter measurements it is shown that the primary interface changes little in terms of atomic-scale roughness even after removal of several hundred angstroms of material. These measurements unequivocally show how a dicarboxylic acid bonds to and may chelate the dissolution of a magnesium-bearing silicate. At the molecular level the solid surface retreat may best be described by a depinning model where Mg is preferentially removed and residual silica tetrahedra polymerize and act to episodically “pin” the surface.     

Microbial mass-dependent fractionation of chromium isotopes

Mass-dependent fractionation of Cr isotopes occurs during dissimilatory Cr(VI) reduction by Shewanella oneidensis strain MR-1. Cells suspended in a simple buffer solution, with various concentrations of lactate or formate added as electron donor, reduced 5 or 10 μM Cr(VI) to Cr(III) over days to weeks. In all nine batch experiments, ⁵³Cr/⁵²Cr ratios of the unreacted Cr(VI) increased as reduction proceeded. In eight experiments covering a range of added donor concentrations up to 100 μM, isotopic fractionation factors were nearly invariant, ranging from 1.0040 to 1.0045, with a mean value somewhat larger than that previously reported for abiotic Cr(VI) reduction (1.0034). One experiment containing much greater donor concentration (10 mM lactate) reduced Cr(VI) much faster and exhibited a lesser fractionation factor (1.0018). These results indicate that ⁵³Cr/⁵²Cr measurements should be effective as indicators of Cr(VI) reduction, either bacterial or abiotic. However, variability in the fractionation factor is poorly constrained and should be studied for a variety of microbial and abiotic reduction pathways.     

Integrating diverse methods to understand climate-land interactions in East Africa

The questions of how land use change affects climate, and how climate change affects land use, require examination of societal and environmental systems across space at multiple scales, from the global climate to regional vegetative dynamics to local decision making by farmers and herders. It also requires an analysis of causal linkages and feedback loops between systems. These questions and the conceptual approach of the research design of the Climate-Land Interaction Project (CLIP) are rooted in the classical human-environment research tradition in Geography.This paper discusses a methodological framework to quantify the two-way interactions between land use and regional climate systems, using ongoing work by a team of multi-disciplinary scientists examining climate-land dynamics at multiple scales in East Africa. East Africa is a region that is undergoing rapid land use change, where changes in climate would have serious consequences for people’s livelihoods, and requiring new coping and land use strategies. The research involves exploration of linkages between two important foci of global change research, namely, land use/land cover (LULC) and climate change. These linkages are examined through modeling agricultural systems, land use driving forces and patterns, the physical properties of land cover, and the regional climate. Both qualitative and quantitative methods are being used to illustrate a diverse pluralism in scientific discovery.     

Oxygen Isotope Evidence for Chemical Interaction of Ki lauea Historical Magmas with Basement Rocks

Klauea historical summit lavas have a wide range in matrix ¹⁸O_VSMOWvalues (4·9–5·6) with lower values in rockserupted following a major summit collapse or eruptive hiatus.In contrast, ¹⁸O values for olivines in most of these lavasare nearly constant (5·1 ± 0·1). The disequilibriumbetween matrix and olivine ¹⁸O values in many samples indicatesthat the lower matrix values were acquired by the magma afterolivine growth, probably just before or during eruption. BothMauna Loa and Klauea basement rocks are the likely sources ofthe contamination, based on O, Pb and Sr isotope data. However,the extent of crustal contamination of Klauea historical magmasis probably minor (< 12%, depending on the assumed contaminant)and it is superimposed on a longer-term, cyclic geochemicalvariation that reflects source heterogeneity. Klauea's heterogeneoussource, which is well represented by the historical summit lavas,probably has magma ¹⁸O values within the normal mid-ocean ridgebasalt mantle range (5·4–5·8) based on thenew olivine ¹⁸O values. KEY WORDS: Hawaii; Klauea; basalt; oxygen isotopes; crustal contamination