Fluid Geochemistry as Indicator of Tectonically-Related,Deep Water Circulations in the Sardinian Rift-Campidano Graben (Italy): New Insights from Environmental Isotopes

The EC funded Geochemical Seismic Zonation program (EEC GSZ Project 1996–1998) chose Sardinia as a low-seismicity site, in which the relationships between fluid geochemistry and seismo-tectonics had to be investigated and results compared with outcomes from other selected high-seismicity sites. A first article, examining the role of fault segmentation and seismic quiescence on the geochemical composition of groundwaters and gases, has already been presented (Angelone et al. 2005). This article deals with environmental isotopes which, together with selected hydrochemical data, give hints on tectonically-related fluid circulations. Four water-dominated hydrothermal systems were considered, all located along regional fault systems and discharging groundwaters belonging to the Na–HCO₃ and Na–Cl facies. In the considered systems, groundwater circulation takes place, principally, in the Palaeozoic Crystalline Basement (PCB), with the exception of the Logudoro system, where hydrological circuits develop in the Mesozoic Carbonate Platform (MCP). The high CO₂ contents, the non-attainment of fluid-rock equilibrium and the large lithological variability prevent the construction of a unique hydrogeological–geochemical conceptual model. In this case, stable isotopes provide a useful tool to describe the origin of fluids and their subterranean movements. Stable isotopes of water, integrated with hydrochemical data, indicate that fluids are derived from three main end members. The dominant component is a relatively recent local meteoric water; the second one is marine water; and the third one is a fossil freshwater, depleted in heavy isotopes with respect to modern rains. The latter end member entered the aquifer system in the past, when climatic conditions were greatly different from today. At least two circulation systems can be recognised, namely a shallow cold system and a deep hydrothermal system, as well as two distinct hydrological processes: (1) gravity-controlled descent of cold water towards greater depths and (2) convection linked to a thermal gradient, causing deep fluids to rise up from the hydrothermal reservoir towards the surface. The highly variable δ¹³C_TDIC values suggest the presence of two distinct CO₂ sources, namely, a biogenic one and a thermogenic one. The relation between the isotopic compositions of CO₂ and He indicates an increased mantle signature in uprising CO₂-rich fluids.     

Evaluation of bank filtration as a pretreatment method for the provision of hygienically safe drinking water in Norway: results from monitoring at two full-scale sites

Two case studies were carried out in central Norway in order to assess the performance of bank filtration systems in cold-climate fluvial aquifers relying on recharge from humic-rich surface waters with moderate microbial contamination. Three municipal wells and two surface-water sources at operative bank filtration systems were monitored for naturally occurring bacteriophages, fecal indicators, natural organic matter (NOM) and physico-chemical water quality parameters during a 4-month period. Aquifer passage effectively reduced the microorganism and NOM concentrations at both study sites. Bacteriophages were detected in 13 of 16 (81%) surface-water samples and in 4 of 24 (17%) well-water samples, and underwent 3 ± 0.3 log₁₀ reduction after 50–80-m filtration and 20–30 days of subsurface passage. NOM reductions (color: 74–97%; dissolved organic carbon: 54–80%; very hydrophobic acids: 70%) were similar to those achieved by conventional water-treatment processes and no further treatment was needed. Both groundwater dilution and sediment filtration contributed to the hygienic water quality improvements, but sediment filtration appeared to be the most important process with regard to microbial and NOM reductions. A strengths-weaknesses-opportunities-threats analysis showed that bank filtration technology has a high potential as a pretreatment method for the provision of hygienically safe drinking water in Norway.     

Investigating the impact of the properties of pilot points on calibration of groundwater models: case study of a karst catchment in Rote Island,Indonesia

A robust configuration of pilot points in the parameterisation step of a model is crucial to accurately obtain a satisfactory model performance. However, the recommendations provided by the majority of recent researchers on pilot-point use are considered somewhat impractical. In this study, a practical approach is proposed for using pilot-point properties (i.e. number, distance and distribution method) in the calibration step of a groundwater model. For the first time, the relative distance–area ratio (d/A) and head-zonation-based (HZB) method are introduced, to assign pilot points into the model domain by incorporating a user-friendly zone ratio. This study provides some insights into the trade-off between maximising and restricting the number of pilot points, and offers a relative basis for selecting the pilot-point properties and distribution method in the development of a physically based groundwater model. The grid-based (GB) method is found to perform comparably better than the HZB method in terms of model performance and computational time. When using the GB method, this study recommends a distance–area ratio of 0.05, a distance–x-grid length ratio (d/X _grid) of 0.10, and a distance–y-grid length ratio (d/Y _grid) of 0.20.     

Source partitioning of anthropogenic groundwater nitrogen in a mixed-use landscape,Tutuila, American Samoa

This study presents a modeling framework for quantifying human impacts and for partitioning the sources of contamination related to water quality in the mixed-use landscape of a small tropical volcanic island. On Tutuila, the main island of American Samoa, production wells in the most populated region (the Tafuna-Leone Plain) produce most of the island’s drinking water. However, much of this water has been deemed unsafe to drink since 2009. Tutuila has three predominant anthropogenic non-point-groundwater-pollution sources of concern: on-site disposal systems (OSDS), agricultural chemicals, and pig manure. These sources are broadly distributed throughout the landscape and are located near many drinking-water wells. Water quality analyses show a link between elevated levels of total dissolved groundwater nitrogen (TN) and areas with high non-point-source pollution density, suggesting that TN can be used as a tracer of groundwater contamination from these sources. The modeling framework used in this study integrates land-use information, hydrological data, and water quality analyses with nitrogen loading and transport models. The approach utilizes a numerical groundwater flow model, a nitrogen-loading model, and a multi-species contaminant transport model. Nitrogen from each source is modeled as an independent component in order to trace the impact from individual land-use activities. Model results are calibrated and validated with dissolved groundwater TN concentrations and inorganic δ¹⁵N values, respectively. Results indicate that OSDS contribute significantly more TN to Tutuila’s aquifers than other sources, and thus should be prioritized in future water-quality management efforts.     

Probabilistic domain decomposition for the solution of the two-dimensional magnetotelluric problem

Probabilistic domain decomposition is proposed as a novel method for solving the two-dimensional Maxwell’s equations as used in the magnetotelluric method. The domain is split into non-overlapping sub-domains and the solution on the sub-domain boundaries is obtained by evaluating the stochastic form of the exact solution of Maxwell’s equations by a Monte-Carlo approach. These sub-domains can be naturally chosen by splitting the sub-surface domain into regions of constant (or at least continuous) conductivity. The solution over each sub-domain is obtained by solving Maxwell’s equations in the strong form. The sub-domain solver used for this purpose is a meshless method resting on radial basis function-based finite differences. The method is demonstrated by solving a number of classical magnetotelluric problems, including the quarter-space problem, the block-in-half-space problem and the triangle-in-half-space problem.     

Modeling core-scale permeability anisotropy in highly bioturbated “tight oil” reservoir rocks

A high-resolution simulation model of a heterogeneous low-permeability rock sample is used to investigate the effects of physical and biogenic sedimentary structures on scaling and anisotropy of absolute permeability at the core scale. Several simulation sub-samples with random locations and volumes were also selected for evaluation of the effects of scale and lithological composition on the calculated permeability. Vertical and horizontal permeability values (from whole core simulation) are in good agreement with routine core analysis (RCA) measurements from offsetting cores. Despite relatively good reservoir quality associated with geobodies of biogenic and relic bedding structures, results from the full diameter core simulation demonstrate that their limited volumetric abundance and restricted connectivity prevent these features from controlling fluid flow in these rocks. In fact, permeability seems to be dominated by the tighter encasing matrix, which exhibits average permeability values very close to those reported from RCA. Geometric averaging offers a better representation for the upscaling of horizontal permeability datasets; whereas, both geometric and harmonic averaging work similarly well for the vertical measurements. The methodology used in this work is particularly applicable to the detailed characterization of reservoir rocks with a high degree of heterogeneity caused by biological reworking and diagenesis.     

Effects of Microparameters on Macroparameters of Flat-Jointed Bonded-Particle Materials and Suggestions on Trial-and-Error Method

The aims of this paper are to analyze the effects of microparameters on macroparameters of flat-jointed bonded-particle materials of particle flow code in two dimensions (PFC2D) and improve the efficiency of trial-and-error method in calibration of microparameters. The trial-and-error method is always used to calibrate the microparameters for a PFC2D model, but it is time-consuming and empirical. To improve the efficiency of the trial-and-error method, the effects of microparameters on macroparameters need to be analyzed. By using the trial-and-error method, a desirable set of microparameters for a uniaxial compression PFC2D model (with a flat-joint contact model as the contact model) of limestone were calibrated to match the macroscopic properties obtained from laboratory tests. Based on the calibration set of microparameters, effects of microparameters on macroparameters of the PFC2D model were studied by the orthogonal design method. On this basis, a flowchart of calibrating microparameters via trial-and-error method was proposed. This flowchart was used to calibrate the microparameters of granite. The macroparameters and failure characteristics determined from numerical simulation were very similar to those of laboratory tests, confirming the effectiveness of the proposed flowchart.     

Landslide susceptibility mapping of the Sera River Basin using logistic regression model

Of the natural hazards in Turkey, landslides are the second most devastating in terms of socio-economic losses, with the majority of landslides occurring in the Eastern Black Sea Region. The aim of this study is to use a statistical approach to carry out a landslide susceptibility assessment in one area at great risk from landslides: the Sera River Basin located in the Eastern Black Sea Region. This paper applies a multivariate statistical approach in the form of a logistics regression model to explore the probability distribution of future landslides in the region. The model attempts to find the best fitting function to describe the relationship between the dependent variable, here the presence or absence of landslides in a region and a set of independent parameters contributing to the occurrence of landslides. The dependent variable (0 for the absence of landslides and 1 for the presence of landslides) was generated using landslide data retrieved from an existing database and expert opinion. The database has information on a few landslides in the region, but is not extensive or complete, and thus unlike those normally used for research. Slope, angle, relief, the natural drainage network (including distance to rivers and the watershed index) and lithology were used as independent parameters in this study. The effect of each parameter was assessed using the corresponding coefficient in the logistic regression function. The results showed that the natural drainage network plays a significant role in determining landslide occurrence and distribution. Landslide susceptibility was evaluated using a predicted map of probability. Zones with high and medium susceptibility to landslides make up 38.8 % of the study area and are located mostly south of the Sera River Basin and along streams.     

Flash flood vulnerability assessment for small catchments with a material flow approach

Destructions resulted from natural hazards like earthquake, landslide, or flood in the urban roads and lifelines introduce their negative effects including the psychological damage to citizens as well as decreased urban functions that usually last for a long time. Thus, a quick and efficient recovery of infrastructures, lifelines, and service-providing facilities along with reducing reconstruction costs and time are essential. This paper proposed an approach that consists of four models for forming an algorithm in order to quantitating and integrating of the criteria that have decisive influence in the recovery of urban roadways after a natural disaster. Meanwhile, to aggregate and conclude the data that are collected by means of presented functions and formulations, we applied fuzzy VIKOR technique as a compromise ranking method. The model outputs a priority list showing the revival of which urban paths stands in higher priority for recovery operation after a natural disaster. Results show that not only the model is able to precisely quantize the selected criteria and provide an action plan for post-event recovery prioritization, but also it offers an appropriate order of transportation roads priority for recovery operations. Finally, the results from the recovery model application to a roadway system in Tehran area are provided.     

An approach for estimating the largest probable tsunami from far-field subduction zone earthquakes

Following the recent unexpected earthquake events of 2004 and 2011, it can be cautiously extrapolated that all major subduction zones bearing the capacity to produce mega-earthquake events will eventually do so given enough time, irrespective of the lack of such in the relatively short historical record. This notion has led to an effort of assigning maximum earthquake magnitudes to all major subduction zones, either based on geological constraints or based on size–frequency relations, or a combination of both. In this study, we utilize the proposed maximum magnitudes to assess tsunami hazard in Central California in the very long return periods. We also assessed tsunami hazard following an alternative methodology to calculate maximum magnitudes, which uses scaling relations for subduction zone earthquakes and maximum fault rupture scenarios found in literature. A sensitivity analysis is performed for Central California that is applicable to any coastal site in the Pacific Rim and can readily provide a strong indication for which subduction zones beam the most energy toward a study area. The maximum earthquake scenarios are then narrowed down to a few candidates, for which the initial conditions are examined in more detail. The chosen worst-case scenarios for Central California stem from the Alaska–Aleutian subduction zone that beams more energy and generates the biggest amplitude waves toward the study area. The largest tsunami scenario produces maximum free surface elevations of 15 m and run-up heights greater than 20 m.