Delineation of groundwater protection zones by the backward particle tracking method: theoretical background and GIS-based stochastic analysis

The backward particle tracking method, an effective and powerful tool that can be used to delineate groundwater protection zones, is presented. The theoretical background and insights on the applicability of this method are provided. Moreover, the present work enriches the backward particle tracking method with an uncertainty analysis concerning the porosity values, applying a Monte Carlo (MC) approach, coupled with the use of geographical information systems (GIS). As an application example, a wellfield in the Komotini area, Greece, is investigated. The present study may serve as a potential guideline for wellfield delineation, particularly in areas like Greece where lack of data related to the hydrogeological system is often a problem.     

The double solid reactant method: II. An application to the shallow groundwaters of the Porto Plain,Vulcano Island (Italy)

This paper presents an example of application of the double solid reactant method (DSRM) of Accornero and Marini (Environmental Geology, 2007a), an effective way for modeling the fate of several dissolved trace elements during water–rock interaction. The EQ3/6 software package was used for simulating the irreversible water–rock mass transfer accompanying the generation of the groundwaters of the Porto Plain shallow aquifer, starting from a degassed diluted crateric steam condensate. Reaction path modeling was performed in reaction progress mode and under closed-system conditions. The simulations assumed: (1) bulk dissolution (i.e., without any constraint on the kinetics of dissolution/precipitation reactions) of a single solid phase, a leucite-latitic glass, and (2) precipitation of amorphous silica, barite, alunite, jarosite, anhydrite, kaolinite, a solid mixture of smectites, fluorite, a solid mixture of hydroxides, illite-K, a solid mixture of saponites, a solid mixture of trigonal carbonates and a solid mixture of orthorhombic carbonates. Analytical concentrations of major chemical elements and several trace elements (Cr, Mn, Fe, Ni, Cu, Zn, As, Sr and Ba) in groundwaters were satisfactorily reproduced. In addition to these simulations, similar runs for a rhyolite, a latite and a trachyte permitted to calculate major oxide contents for the authigenic paragenesis which are comparable, to a first approximation, with the corresponding data measured for local altered rocks belonging to the silicic, advanced argillic and intermediate argillic alteration facies. The important role played by both the solid mixture of trigonal carbonates as sequestrator of Mn, Zn, Cu and Ni and the solid mixture of orthorhombic carbonates as scavenger of Sr and Ba is emphasized.

Timing of UHP metamorphism in the Hong’an area,western Dabie Mountains,China: evidence from zircon U–Pb age,trace element and Hf isotope composition

The Hong’an area (western Dabie Mountains) is the westernmost terrane in the Qinling-Dabie-Sulu orogen that preserves UHP eclogites. The ages of the UHP metamorphism have not been well constrained, and thus hinder our understanding of the tectonic evolution of this area. LA-ICPMS U–Pb age, trace element and Hf isotope compositions of zircons of a granitic gneiss and an eclogite from the Xinxian UHP unit in the Hong’an area were analyzed to constrain the age of the UHP metamorphism. Most zircons are unzoned or show sector zoning. They have low trace element concentrations, without significant negative Eu anomalies. These metamorphic zircons can be further subdivided into two groups according to their U–Pb ages, and trace element and Lu–Hf isotope compositions. One group with an average age of 239 ± 2 Ma show relatively high and variable HREE contents (527 ≥ Lu_N ≥ 14) and ¹⁷⁶Lu/¹⁷⁷Hf ratios (0.00008–0.000931), indicating their growth prior to a great deal of garnet growth in the late stage of continental subduction. The other group yields an average age of 227 ± 2 Ma, and shows consistent low HREE contents and ¹⁷⁶Lu/¹⁷⁷Hf ratios, suggesting their growth with concurrent garnet crystallization and/or recrystallization. These two groups of age are taken as recording the time of prograde HP to UHP and retrograde UHP–HP stages, respectively. A few cores have high Th/U ratios, high trace element contents, and a clear negative Eu anomaly. These features support a magmatic origin of these zircon cores. The upper intercept ages of 771 ± 86 and 752 ± 70 Ma for the granitic gneiss and eclogite, respectively, indicate that their protoliths probably formed as a bimodal suite in rifting zones in the northern margin of the Yangtze Block. Young Hf model ages (T _DM1) of magmatic cores indicate juvenile (mantle-derived) materials were involved in their protolith formation. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Mapping ambivalence: Exploring the geographies of community change and rails-to-trails development using photo-based Q method and PPGIS

As the literature on trail development suggests, recreational trail projects can generate conflicts and controversies, particularly when built on abandoned rail corridors through developed areas. These conflicts are often understood as “not in my back yard” (NIMBY) reactions, suggesting a spatial proximity to conflict which increases as one draws closer to the proposed trail. This research seeks to understand local residents’ perceptions and reactions to recreational trail development in the City of Delaware (Ohio, USA). It addresses two spatially infused questions: Does the potential for conflict related to trail development increase as people live closer to a potential trail (the NIMBY factor)? Can important qualitative factors about favorable and unfavorable land uses including potential recreational trail sites be defined using a participatory methodology and then represented in GIS? The study used a mixed-method approach to collect and analyze qualitative data from a group of local residents. Each participant was interviewed and asked to sort 19 pictures related to trail development. After each of the sorts, participants were asked to explain why they ranked the pictures the way they did. Results of the picture sorts were then analyzed using Q method and mapped with GIS. The results show that spatial proximity matters in the context of trail development and potential NIMBY reactions to trails. Significant differences were found in the picture sorts that reveal the importance of proximity and location, although in a manner contrary to the assumptions in the writings on rails-to-trails. Through combining qualitative methods, Q analysis and PPGIS analysis, the research shows that qualitative place-based studies are capable of generating insights about the complexities of situated geographic change such as recreational trail development.     

Geochronology of Pliocene volcanism in the Dzhavakheti Highland (the Lesser Caucasus). Part 1: Western part of the Dzhavakheti Highland

Isotopic-geochronological study of the Pliocene magmatic activity in western part of the Dzhavakheti Highland (northwestern region of the Lesser Caucasus) is carried out. The results obtained imply that the Pliocene magmatic activity lasted in this part of the highland approximately 2 million years from 3.75 to 1.75–1.55 Ma. As is established, the studied volcanic rocks correspond in composition mostly to K-Na subalkaline and more abundant normal basalts. Time constraints of main phases in development of basic volcanism within the study region are figured out. We assume that individual pulses of silicic to moderately silicic volcanism presumably took place in the Dzhavakheti Highland about 3.2 and 2.5 Ma ago.     

Numerical estimation of REV and permeability tensor for fractured rock masses by composite element method

The Monte Carlo method is used to generate parent stochastic discrete fracture network, from which a series of fractured rock samples of different sizes and orientations are extracted. The fracture network combined with a regular grid forms composite element mesh of the fractured rock sample, in which each composite element is composed of sub‐elements incised by fracture segments. The composite element method (CEM) for the seepage is implemented to obtain the nodal hydraulic potential as well as the seepage flow rates through the fractured rock samples. The application of CEM enables a large quantity of stochastic tests for the fractured rock samples because the pre‐process is facilitated greatly. By changing the sizes and orientations of the samples, the analysis of the seepage characteristics is realized to evaluate the variation of the permeability components, the existence of the permeability tensor and the representative element volume. The feasibility and effectiveness are illustrated in a numerical example. Copyright © 2008 John Wiley & Sons, Ltd.     

Finite Element Modeling of Swelling Soil Pressures Using Joint/Interface Elements

There could be some discontinuities in a soil media such as layers, earth structures, cracks, and fissures so that estimation of stresses and deformations in these types of soil masses are somewhat different from continuous masses. The discontinuities in a soil mass could be considered as a special link between two blocks. Transmitted swelling pressures affected by the soil properties of the backfill especially at the interface between the backfill and retaining wall. In order to estimate the transmitted swelling pressure distribution behind a retaining wall, using with interface element, a new finite element model and a code (SWELPRES) have been developed. To evaluate the effect of backfill characteristics and interface element from the transmitted lateral swelling pressures, four cases of thickness of backfills with or without interface elements are investigated.     

江苏潮滩沉积物激光粒度仪与移液管-筛析分析结果的对比

激光粒度仪的广泛应用带来了与历史数据的对比问题,因此需要建立激光粒度仪与早期分析结果之间的关系。根据江苏海岸潮滩沉积物样品的激光粒度仪和移液管-筛析法分析,对粒度参数、粒度组分等进行了对比,对两种方法的差异进行了分析。结果表明,对于江苏潮滩沉积物粒度参数中平均粒径的激光粒度仪与移液管-筛析分析结果之间有良好的线性关系;筛分法测得的粗颗粒物质较激光法偏少,而移液管法法测得的细颗粒物质较激光法偏多;将样品分类之后再进行两种方法所获粒度参数的回归分析,相关性得以提高,说明不同粒度组成的沉积物对分析结果的对比有不同的影响。两种方法之间的换算关系不仅与研究区域有关,而且与沉积物本身的粒度组成有关。