A GIS-based approach for gully erosion susceptibility modelling: a test in Sicily, Italy

The aim of this study is to analyze the susceptibility conditions to gully erosion phenomena in the Magazzolo River basin and to test a method that allows for driving the factors selection. The study area is one of the largest (225 km²) watershed of southern Sicily and it is mostly characterized by gentle slopes carved into clayey and evaporitic sediments, except for the northern sector where carbonatic rocks give rise to steep slopes. In order to obtain a quantitative evaluation of gully erosion susceptibility, statistical relationships between the spatial distributions of gullies affecting the area and a set of twelve environmental variables were analyzed. Stereoscopic analysis of aerial photographs dated 2000, and field surveys carried out in 2006, allowed us to map about a thousand landforms produced by linear water erosion processes, classifiable as ephemeral and permanent gullies. The linear density of the gullies, computed on each of the factors classes, was assumed as the function expressing the susceptibility level of the latter. A 40-m digital elevation model (DEM) prepared from 1:10,000-scale topographic maps was used to compute the values of nine topographic attributes (primary: slope, aspect, plan curvature, profile curvature, general curvature, tangential curvature; secondary: stream power index; topographic wetness index; LS-USLE factor); from available thematic maps and field checks three other physical attributes (lithology, soil texture, land use) were derived. For each of these variables, a 40-m grid layer was generated, reclassifying the topographic variables according to their standard deviation values. In order to evaluate the controlling role of the selected predictive variables, one-variable susceptibility models, based on the spatial relationships between each single factor and gullies, were produced and submitted to a validation procedure. The latter was carried out by evaluating the predictive performance of models trained on one half of the landform archive and tested on the other. Large differences of accuracy were verified by computing geometric indexes of the validation curves (prediction and success rate curves; ROC curves) drawn for each one-variable model; in particular, soil texture, general curvature and aspect demonstrated a weak or a null influence on the spatial distribution of gullies within the studied area, while, on the contrary, tangential curvature, stream power index and plan curvature showed high predictive skills. Hence, predictive models were produced on a multi-variable basis, by variously combining the one-variable models. The validation of the multi-variables models, which generally indicated quite satisfactory results, were used as a sensitivity analysis tool to evaluate differences in the prediction results produced by changing the set of combined physical attributes. The sensitivity analysis pointed out that by increasing the number of combined environmental variables, an improvement of the susceptibility assessment is produced; this is true with the exception of adding to the multi-variables models a variable, as slope aspect, not correlated to the target variable. The addition of this attribute produces effects on the validation curves that are not distinguishable from noise and, as a consequence, the slope aspect was excluded from the final multi-variables model used to draw the gully erosion susceptibility map of the Magazzolo River basin. In conclusion, the research showed that the validation of one-variable models can be used as a tool for selecting factors to be combined to prepare the best performing multi-variables gully erosion susceptibility model.     

Comprehensive chemical analyses of natural cordierites: implications for exchange mechanisms

Cordierite samples from pegmatites and metamorphic rocks have been analysed for major [electron microprobe analysis (EMPA)] and trace elements [inductively coupled plasma mass spectrometry (ICP-MS), secondary ion mass spectrometry analyses (SIMS)] as well as for H₂O and CO₂ (coulometric titration), and the results evaluated in conjunction with published data in order to determine which exchange mechanisms are significant. Apart from the homovalent substitutions FeMg₋₁ and MnMg₋₁ on the octahedral site, some minor KNa₋₁ on the Ch0 channel site, and Fe³⁺Al₋₁ on the T₁1 tetrahedral site, the three most important substitution mechanisms are those for the incorporation of Li on the octahedral sites (NaLi□₋₁Mg₋₁), and of Be and other divalent cations on the tetrahedral T₁1 site (NaBe□₋₁Al₋₁ and Na(Mg,Fe²⁺)□₋₁Al₋₁). The dominant role of the last vector is clearly demonstrated. We propose a new generalized formula for cordierite: ^Ch(Na,K)_0–1 ^VI(Mg,Fe²⁺,Mn,Li)₂ ^IVSi₅ ^IVAl₃ ^IV(Al, Be, Mg, Fe²⁺, Fe³⁺)O₁₈ *x^Ch(H₂O, CO₂…). Our results show that the population of (Mg, Fe²⁺) on the T₁1-site is limited to about 0.08 a.p.f.u. Other exchange mechanisms that were encountered in experiments operate only under P–T conditions or in bulk compositions that are rarely realized in nature. Routine analyses by electron microprobe in which Li and Be are not determined can be plotted as (Mg+Fe+Mn) versus (Si+Al) to assess whether significant amounts of Li and Be could be present. These amounts can be calculated as Li (a.p.f.u.)=Al+Na–4 and Be (a.p.f.u.)=10–2Al–M²⁺–Na.     

A New Method for C k -Surface Approximation from a Set of Curves, with Application to Ship Track Data in the Marianas Trench

We introduce a surface approximation technique to address the problem of fitting a surface to a given set of curves. The originality of the method lies in its ability to take into account the continuous aspect of the data, and also in the possibility to arbitrarily select the regularity (C⁰, C¹, or higher) of the approximant obtained. We demonstrate the efficiency of the approach by constructing a bathymetry map of the Marianas trench based upon a set of SONAR (SOnic Navigation And Ranging) bathymetry ship track data.     

Global change revealed by palaeolimnological records from remote lakes: a review

Over recent decades, palaeolimnological records from remote sites have provided convincing evidence for the onset and development of several facets of global environmental change. Remote lakes, defined here as those occurring in high latitude or high altitude regions, have the advantage of not being overprinted by local anthropogenic processes. As such, many of these sites record broad-scale environmental changes, frequently driven by regime shifts in the Earth system. Here, we review a selection of studies from North America and Europe and discuss their broader implications. The history of investigation has evolved synchronously with the scope and awareness of environmental problems. An initial focus on acid deposition switched to metal and other types of pollutants, then climate change and eventually to atmospheric deposition-fertilising effects. However, none of these topics is independent of the other, and all of them affect ecosystem function and biodiversity in profound ways. Currently, remote lake palaeolimnology is developing unique datasets for each region investigated that benchmark current trends with respect to past, purely natural variability in lake systems. Fostering conceptual and methodological bridges with other environmental disciplines will upturn contribution of remote lake palaeolimnology in solving existing and emerging questions in global change science and planetary stewardship.     

Modeling surface water–groundwater interaction in New Zealand: Model development and application

Most rivers worldwide have a strong interaction with groundwater when they leave the mountains and flow over alluvial plains before flowing into the seas or disappearing in the deserts, and in New Zealand, typically, rivers lose water to the groundwater in the upper plains and generally gain water from the groundwater in the lower plains. Aiming at simulating surface water–groundwater interaction nationally in New Zealand, we developed a conceptual groundwater module for the national hydrologic model TopNet to simulate surface water–groundwater interaction, groundwater flow, and intercatchment groundwater flow. The developed model was applied to the Pareora catchment in South Island of New Zealand, where there are concurrent spot gauged flows. Results show that the model simulations not only fit quite well to flow measurement but also to concurrent spot gauged flows, and compared to the original TopNet, it has a significant improvement in the low flows. Sensitivity analysis shows river flow is sensitive to the river losing/gaining rate instead of groundwater characteristic, while groundwater storage is sensitive to both river losing/gaining rate and groundwater characteristic. This indicates our conceptual approach is promising for nationwide modeling without the large amount of geology and aquifer data typically required by physically‐based modeling approaches.     

Recovering physical parameters from galaxy spectra using MOPED

We derive physical parameters of galaxies from their observed spectra using MOPED, the optimized data compression algorithm of Heavens, Jimenez & Lahav. Here we concentrate on parametrizing galaxy properties, and apply the method to the NGC galaxies in Kennicutt's spectral atlas. We focus on deriving the star formation history, metallicity and dust content of galaxies. The method is very fast, taking a few seconds of CPU time to estimate ∼17 parameters, and is therefore specially suited to studying large data sets, such as the Anglo-Australian two-degree-field (2dF) galaxy survey and the Sloan Digital Sky Survey (SDSS). Without the power of MOPED, the recovery of star formation histories in these surveys would be impractical. In Kennicutt's atlas, we find that for the spheroidals a small recent burst of star formation is required to provide the best fit to the spectrum. There is clearly a need for theoretical stellar atmospheric models with spectral resolution better than 1 Å if we are to extract all the rich information that large redshift surveys contain in their galaxy spectra.     

Linked fluid and tectonic evolution in the High Himalaya mountains (Nepal)

Fluid inclusions were studied in a quartz lens from the structurally highest unit of the Himalaya mountains in Nepal from a textural, geometrical, chemical and isotopic point of view. Six types of fluid inclusions were distinguished. One of these types consists of annular inclusions; this shape is attributed to a confining pressure increase in a non-isotropic stress field. Two successive stress fields were deduced from the orientation of the inclusion planes relative to the schistosity. The bulk composition of the fluid was dominated by CO₂ (>84 mol%) and H₂O. The composition remained constant during the whole history of the sample indicating that it was buffered by the carbonaceous host rock and/or that one single fluid was reworked in situ by decrepitation. Stable isotope of fluids and minerals indicate (1) that fluids were buffered by surrounding rocks for O and C and (2) that at least two types of water (metamorphic and meteoric) were involved. Finally, a P-T-t-- path is proposed for the sample, taking into account the southward thrusting along the Main Central Thrust, the northward tectonic denudation of the Himalaya mountains inducing tectonic burying below the Annapurna Range, and lastly, rapid uplift.