Numerical modeling assessment of three-gate structures for capturing contaminated groundwater

This modeling study evaluated the capability of alternative funnel-and-gate structures with three gates for capturing contaminated groundwater in a hypothetical unconfined aquifer. Simulated interceptor structures were linear and 45 m wide, consisting of three gates and two funnels (walls). One gate occupied the center and two gates occupied the ends of the interceptor structures. The structures, positioned perpendicular to regional groundwater flow, traversed the entire thickness of the aquifer. A total of four structures were evaluated (numbers designate widths of end, center, and end gates, respectively, in meters): 3-3-3, 2-5-2, 1-7-1, and 4-1-4. Particle tracking and zonal water budgets identified shapes of capture zones and discharge patterns for each interceptor structure. A mass transport model, accounting for advection and hydrodynamic dispersion, tested the capability of each structure for capturing a contaminant plume. Results suggest that: time-dependent capture zones underestimate the amount of time to capture a contaminant plume, wide center gates facilitate plume capture, and wide end gates facilitate lateral containment of contaminants. Of the structures simulated, the 2-5-2 configuration was relatively efficient at processing and containing the simulated contaminant plume.     

Groundwater recharge of carbonate aquifers of the Silesian-Cracow Triassic (southern Poland) under human impact

A Triassic carbonate unit has been intensively drained by zinc and lead ore mines and numerous borehole fields since the nineteenth century. Its groundwater recharge has increased due to: pumping of water from boreholes, mining activity, and urbanization. An approach to determine the amounts of the recharge at a variety of spatial scales is presented in the paper. Different methods were used to identify and quantify recharge components on a regional and local scale: mathematical modelling was performed for four aquifers included in an aquifer system, an analytical estimation based on the assumption that an average recharge is equal to the average discharge of the hydrogeological system—for six man-made drainage centres, and the method of water level fluctuation (WLF) was applied in one observation borehole. Results of modelling have been supplemented by observation of environmental tracers (δ¹⁸O, δ²H, ³H), noble gases temperatures, and ⁴He_exc in groundwater. The regional aquifer’s current recharge according to estimations performed by means of modelling varies from 39 to 101 mm/year on average. Depending on the aquifer site the average precipitation ranges from 779 to 864 mm/year. In the confined part of the aquifer average recharge ranges from 26 to 61 mm/year. Within outcrops average recharge varies from 96 to 370 mm/year. Current recharge estimated by the analytical method for man-made drainage centres varies from 158 up to 440 mm/year. High values are caused by different recharge sources like precipitation, induced leakage from shallow aquifers, and water losses from streams, water mains and sewer systems. Pumping of water, mining and municipal activities constitute additional factors accounting for the intensified recharge.     

Mapping susceptibility of rainfall-triggered shallow landslides using a probabilistic approach

To prepare a landslide susceptibility map is essential to identify hazardous regions, construct appropriate mitigation facilities, and plan emergency measures for a region prone to landslides triggered by rainfall. The conventional mapping methods require much information about past landslides records and contributing terrace and rainfall. They also rely heavily on the quantity and quality of accessible information and subjectively of the map builder. This paper contributes to a systematic and quantitative assessment of mapping landslide hazards over a region. Geographical Information System is implemented to retrieve relevant parameters from data layers, including the spatial distribution of transient fluid pressures, which is estimated using the TRIGRS program. The factor of safety of each pixel in the study region is calculated analytically. Monte Carlo simulation of random variables is conducted to process the estimation of fluid pressure and factor of safety for multiple times. The failure probability of each pixel is thus estimated. These procedures of mapping landslide potential are demonstrated in a case history. The analysis results reveal a positive correlation between landslide probability and accumulated rainfall. This approach gives simulation results compared to field records. The location and size of actual landslide are well predicted. An explanation for some of the inconsistencies is also provided to emphasize the importance of site information on the accuracy of mapping results.     

Experimental study on the deterioration and natural remediation of the Ariake Sea tidal mud caused by the sea laver treatment acid practice and the upward seepage of pore water liquid

This study presents a laboratory study of the following two aspects: (1) the influence of sea laver treatment acid on the geoenvironmental properties of Ariake Sea tidal mud, and (2) the natural remediation effect on the sea laver treatment acid contaminated Ariake Sea tidal mud caused by the upward seepage of pore water liquid in the mud. Firstly, the mechanisms of the transport of sea laver treatment acid in the Ariake Sea tidal mud and the generation mechanisms of the upward seepage flow in the Ariake Sea tidal mud are discussed. Secondly, a series of one-dimensional laboratory infiltration tests were carried out to investigate the deterioration of the Ariake Sea tidal mud caused by the sea laver acid treatment practice. Test results reveal that the acid treatment practice caused considerable change in the geochemical properties of the mud in terms of increase in sulfide content and decrease in pH value. After the treatment by the sea laver treatment acid, the sulfide content of the mud even exceeded the safe limit value for the benthos, which represents undesirable living condition for benthos. Thirdly, series of laboratory fresh seawater infiltration tests for the deteriorated Iida site mud were conducted to illustrate this natural remediation efficiency. It is found that with the infiltration of the fresh seawater, the sulfide content of the Iida site mud was considerably reduced and pH value increased to an acceptable range for benthos living in the tidal flat mud. With the increase in the infiltration time and the hydraulic gradient, the remediation efficiency could be increased.     

Inclusions of silicate and sulfate melts in chrome diposide from the Inagli deposit,Yakutia, Russia

Melt inclusions were studied in chrome diopside from the Inagli deposit of gemstones in the Inagli massif of alkaline ultrabasic rocks of potassic affinity in the northwestern Aldan shield, Yakutia, Russia. The chrome diopside is highly transparent and has an intense green color. Its Cr₂O₃ content varies from 0.13 to 0.75 wt %. Primary and primary-secondary polyphase inclusions in chrome diopside are dominated by crystal phases (80–90 vol %) and contain aqueous solution and a gas phase. Using electron microprobe analysis and Raman spectroscopy, the following crystalline phases were identified. Silicate minerals are represented by potassium feldspar, pectolite [NaCa₂Si₃O₈(OH)], and phlogopite. The most abundant minerals in the majority of inclusions are sulfates: glaserite (aphthitalite) [K₃Na(SO₄)₂], glauberite [Na₂Ca(SO₄)₂], aluminum sulfate, anhydrite (CaSO₄), gypsum (CaSO₄ × 2H₂O), barite (BaSO₄), bloedite [Na₂Mg(SO₄)₂ × 4H₂O], thenardite (NaSO₄), polyhalite [K₂Ca₂Mg(SO₄)₄ × 2H₂O], arcanite (K₂SO₄), and celestite (SrSO₄). In addition, apatite was detected in some inclusions. Chlorides are probably present among small crystalline phases, because some analyses of aggregates of silicate and sulfate minerals showed up to 0.19–10.3 wt % Cl. Hydrogen was identified in the gas phase of polyphase inclusions by Raman spectroscopy. The composition of melt from which the chrome diopside crystallized was calculated on the basis of the investigation of silicate melt inclusions. This melt contains 53.5 wt % SiO₂, considerable amounts of CaO (16.3 wt %), K₂O (7.9 wt %), Na₂O (3.5 wt %), and SO₃ (1.4 wt %) and moderate amounts of Al₂O₃ (7.5 wt %), MgO (5.8 wt %), FeO (1.1 wt %), and H₂O (0.75 wt %). The content of Cr₂O₃ in the melt was 0.13 wt %. Many inclusions were homogenized at 770–850°C, when all of the crystals and the gas phase were dissolved. The material of inclusions heated up to the homogenization temperature became heterogeneous even during very fast quenching (two seconds) producing numerous small crystals. This fact implies that most of the inclusions contained a salt (rather than silicate) melt of sulfate-dominated composition. Such inclusions were formed from salt globules (with a density of about 2.5 g/cm³) occurring as an emulsion in the denser (2.6 g/cm³) silicate melt from which the chrome diopside crystallized.     

Metamorphic rocks of the Murmansk domain (Kola Peninsula) as compared with the oldest rock associations of the northeastern Baltic shield,Canada, and Greenland

Rock associations characterized by heterogeneous sets of petrogeochemical parameters were compared by quantifying the degree of their similarity-dissimilarity and searching for discrimination trends between them. Using procedures specially developed for this purpose, it was demonstrated for the first time that the lithotectonic complexes of the Murmansk domain are fundamentally different from those of typical granulite-gneiss terrains and resemble the granite-greenstone terrains of the Baltic shield, Greenland, and Canada. Based on the whole data set, the Murmansk domain can be considered as a deeply eroded Archean granite-greenstone terrain retaining only the tonalite-trondhjemite-gneiss basement with abundant supracrustal enclaves. A trend of the compositional difference between the older and younger rock associations is similar to that between the tholeiitic and boninitic volcanic series. It was suggested that the petrogeochemical “age” trend reflects the initial stage of the compositional evolution of the metamagmatic rocks of the region from metamorphic rocks similar to tholeiites at the early stages to the Paleoproterozoic boninite-like rocks, which are believed to be linked to the unique PGE-bearing province of the northeastern Baltic shield. This implies that the specific metallogenic features of the region emerged already in the Archean, which supported the suggestion on the long duration of geological processes in the Early Precambrian.     

Conditions of Quaternary magmatism at Spitsbergen Island

Petrological and geochemical data obtained on the Quaternary lavas of volcanoes at Spitsbergen Island indicate that the rocks were produced via the deep-seated crystallization of parental alkaline magmas at 8–10 kbar. The character of clinopyroxene enrichment in incompatible elements indicates that the mineral crystallized from more enriched melts than those inferred from the composition of the host lavas. These melts were close to the parental melts previously found as veinlets in mantle hyperbasite xenoliths in the lavas. According to the character of their enrichment in Pb and Sr radiogenic isotopes and depletion in Nd, the basalts from Spitsbergen Island define a single trend with the weakly enriched tholeiites of the Knipovich Ridge, a fact suggesting the closeness of the enriched sources beneath the continental margin of Spitsbergen and beneath the spreading zone. Magmatic activity at Spitsbergen was related to the evolution of the Norwegian-Greenland basin, which evolved in pulses according to the shift of the spreading axes. The most significant of the latter events took place in the Neogene, when the Knipovich Ridge obtained its modern position near the western boundary of Spitsbergen. Early in the course of the evolution, the emplacement of alkaline melts generated at Spitsbergen into the oceanic mantle could form the enriched mantle, which was later involved in the melting process beneath the spreading zone.     

Stable isotope study of the Archaean rocks of the Vredefort impact structure,central Kaapvaal Craton,South Africa

The Vredefort dome in the Kaapvaal Craton was formed as a result of the impact of a large meteorite at 2.02 Ga. The central core of Archaean granitic basement rocks is surrounded by a collar of uplifted and overturned strata of the Witwatersrand Supergroup, exposing a substantial depth section of the Archaean crust. Orthogneisses of the core show little variation in whole-rock δ ¹⁸O value, with the majority being between 8 and 10‰, with a mean of 9.2‰ (n = 35). Quartz and feldspar have per mil differences that are consistent with O-isotope equilibrium at high temperatures, suggesting minimal interaction with fluids during subsequent cooling. These data refute previous suggestions that the Outer Granite Gneiss (OGG) and Inlandsee Leucogranofels (ILG) of the core represent middle and lower crust, respectively. Granulite-facies greenstone remnants from the ILG have δ ¹⁸O values that are on average 1.5‰ higher than the ILG host rocks and are unlikely, therefore, to represent the residuum from the partial melting event that formed the host rock. Witwatersrand Supergroup sedimentary rocks of the collar, which were metamorphosed at greenschist-to amphibolite-facies conditions, generally have lower δ ¹⁸O values than the core rocks with a mean value for metapelites of 7.7‰ (n = 45). Overall, through an ∼20 km thick section of crust, there is a general increase in whole-rock δ ¹⁸O value with increasing depth. This is the reverse of what is normal in the crust, largely because the collar rocks have δ ¹⁸O values that are unusually low in comparison with metamorphosed sedimentary rocks worldwide. The collar rocks have δD values ranging from −35 to −115‰ (average −62‰, n = 29), which are consistent with interaction with water of meteoric origin, having a δD of about −25 to −45‰. We suggest that fluid movement through the collar rocks was enhanced by impact-induced secondary permeability in the dome structure. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.