Qanat is not a hazard

Comments are presented on the article by Abbasnejad et al. (Environ Earth Sci 75:1306, 2016) dealing with qanat and hazards in Iran. My rebuttal starts from the direct attribution of the described hazards to qanat and addresses the importance in the correct use of terminology for geological hazards. All of the problems pointed out in Iran (subsidence, sinkholes, pollution) have, to me, an anthropogenic origin and cannot be directly ascribed to qanat. Eventually, I present some additional elements to highlight the remarkable importance of qanat systems and their influence on the development of similar underground structures in many countries of the Mediterranean Basin. This remarks the relevance of qanat as cultural heritage sites and the need for their preservation and valorization.     

2D probabilistic prediction of sparsely measured earth properties constrained by geophysical imaging fully accounting for tomographic reconstruction ambiguity

Many hydrological, environmental, or engineering exploration tasks require predicting spatially continuous scenarios of sparsely measured borehole logging data. We present a methodology to probabilistically predict such scenarios constrained by ill-posed geophysical tomography. Our approach allows for transducing tomographic reconstruction ambiguity into the probabilistic prediction of spatially continuous target parameter scenarios. It is even applicable to data sets where petrophysical relations in the survey area are non-unique, i.e., different facies related petrophysical relations may be present. We employ static two-layer artificial neural networks (ANNs) for prediction and additionally evaluate, whether the training performance of the ANNs can be used to rank geophysical tomograms, which are mathematically equal reconstructions of physical parameter distributions in the ground. We illustrate our methodology using a realistic synthetic database for maximal control about the prediction performance and ranking potential of the approach. For doing so, we try to link geophysical radar and seismic tomography as input parameters to porosity of the ground as target parameter of ANN. However, the approach is flexible and can cope with any combination of geophysical tomograms and hydrologic, environmental or engineering target parameters. Ranking of equivalent geophysical tomograms based on additional borehole logging data is found to be generally possible, but risks remain that the ranking based on the ANN training performance does not fully coincide with the closeness of geophysical tomograms to ground truth. Since geophysical field data sets do usually not offer control options similar to those used in our synthetic database, we do not recommend the utilization of recurrent ANNs to learn weights for the individual geophysical tomograms used in the prediction procedure.     

An evaluation of the relationship between physical/mechanical properties and mineralogy of landscape rocks as determined by hyperspectral reflectance

We investigated the relationships between mineral content and the physical and mechanical properties of landscape rock using a non-destructive remote sensing method applied in the laboratory. Using this technique, the spectral properties of the landscape rock could be collected at different wavelengths without harming the samples. Differences in spectral reflectance were compared with the physical and mechanical properties of the stone. Significant correlations were observed between reflectance values and the rocks’ mechanical and physical properties, with correlation coefficients of 95 to 99 %. However, establishing a correlation between two variables is not a sufficient condition to establish a causal relationship. Mineral densities and mineral content are characteristics used for the classification of landscape rock. We have concluded that although spectral signatures from landscape rock can be used for predicting which stones might have similar features when comparing two batches of stone, the high correlations we discovered cannot confirm a cause and effect relationship that would allow for the prediction of a rock’s physical and mechanical properties. Although this conclusion is disappointing, the mineral content and the significant correlations discovered by hyperspectral reflectance scanning can be used as supplementary information when comparing two samples of landscape rock.     

Artificial neural network inversion of magnetic anomalies caused by 2D fault structures

A new approach is proposed to interpret magnetic anomalies caused by 2D fault structures. This approach is based on the artificial neural network inversion, utilizing particularly modular neural network algorithm. The inversion process is implemented to estimate the parameters of 2D fault structures where it has been verified first on synthetic models. The results of the inversion show that the parameters derived from the inversion agree well with the true ones. The analysis of noise has been studied in order to investigate the stability of the approach where it has been tested for contaminated anomalies with 5 and 10 % of white Gaussian noise. The results of the inversion provide satisfactory results even with contaminated signals.The validity of the approach has been demonstrated through real data taken from New South Wales, Australia. A comparable and satisfactory agreement is shown between the inversion results of the neural network and those from techniques published in literatures.     

Using inverse modeling and hierarchical cluster analysis for hydrochemical characterization of springs and Talkhab River in Tang-Bijar oilfield,Iran

Hierarchical cluster analysis (HCA) and inverse modeling (PH REdox EQuilibrium (in C language) (PHREEQC)) were simultaneously useful approaches in interpreting surface water hydrochemistry within Talkhab River in the Tang-Bijar oilfield, Iran, where large uncertainties exist in the understanding of the water quality system. Q-mode HCA applied to the data revealed three major surface water associations distinguished on the basis of the major causes of variation in the hydrochemistry. The three water groups were classified as upstream waters (group 1: Ca–SO₄), intermediate waters (group 2: Ca–SO₄–Cl), and downstream waters (group 3: Na–Cl). Geochemical reaction models were constructed using PHREEQC to establish the reactions associated with the different mineral phases through inverse modeling. The hydrochemical compositions of the water groups and the mass balance calculations indicate that the dominant processes and reactions responsible for the hydrochemical evolution in the system are (1) dissolution of evaporites, (2) precipitation of carbonate minerals, (3) silicate weathering reactions, (4) limited mixing with saline water, and (5) ion exchange.     

Hydrogeological assessment of the Tabarteh fault zone by physicochemical analysis and multivariate statistical methods

Hydrogeologically, faults may impede, conduit, exert no influence, or may play a combination of these roles on groundwater flow. The object of this paper is to study the hydrogeological role of the Tabarteh fault, which is located on the border of Zagros and Central Iran tectonic zones in an alluvial aquifer. The recorded data of water table levels, chemical parameters, and discharge rate of wells, in addition to geological maps and geophysical results, were collected and evaluated. The outcrop of travertine in limited areas and the emergence of a few small springs within the alluvium show a barrier role of the fault in the groundwater flow. The spatial analysis of chemical components, head time series, and groundwater flow direction assessment demonstrated that the fault acts as both a barrier and a non-barrier in different sections. The multivariate statistical methods of cluster and discriminant analyses also confirm the dual role of the fault.     

Characteristics of groundwater flow field after land creation engineering in the hilly and gully area of the Loess Plateau

In order to study the degree of influence and control mechanism to groundwater flow field caused by land creation engineering in the hilly and gully area of the Loess Plateau, based on the geological and engineering conditions of the first stage project of Yan’an new district in China, numerical simulation of groundwater flow is carried out by the Feflow and GIS technologies. From the simulation, punning measure relatively reduces infiltration recharge and artificial gravel drain increases groundwater seepage. The basic characteristics of groundwater flow field is controlled by the old and new topographies in the whole study area, and artificial gravel drain plays an auxiliary role in accelerating groundwater drawdown upstream and promotes groundwater rise downstream. According to differences of groundwater level and declining percentages of hydraulic gradient in the main and secondary gullies, dewatering of artificial gravel drain in the secondary gully is more effective than that in the main gully, which will yet play an important role in the future. The study results will make contributions to understand groundwater response to land creation engineering and will be beneficial to take necessary measures to prevent collapse of loess and failure of building foundation in the hilly and gully area of the Loess Plateau.     

Impact of spate irrigation of flood waters on agricultural drought and groundwater recharge: case of Sidi Bouzid plain,Central Tunisia

In Sidi Bouzid plain, rainfall alone is insufficient to satisfy crop water requirements. Within this framework, and in order to improve water resources in the region, the Tunisian State adopted non-conventional water mobilization techniques, among which artificial spate irrigation. The objective of the study is to evaluate the impact of spate irrigation of flood water on the mitigation of agricultural drought and the enhancement of groundwater recharge. Annual and monthly rainfall data as well as flood water volumes were monitored. The study focused on the groundwater drawdown monitoring. Results showed a high flood water contribution to crop water requirements that exceeded rainfall. This water prevented drought in the spate perimeters. The groundwater drawdown was found to fluctuate over time, with an average decreasing rate of 0.4 to 0.5 m/year. Groundwater recharge was found to be highly correlated with flood water contribution through spate irrigation (R ² = 84 %). Out of the spate zone, a high decrease in the groundwater level was noted. The lowest rate of 1 m/year was that of the farthest piezometer from the spate perimeters. This is influenced by the excessive pumping out of the spate zone. In 1980, groundwater flew from the west to the east. In 2015, the flow movement from the east to the center of the plain did not change due to the presence of the spate perimeters. Nevertheless, excessive pumping around sabkhas changed the flow directions at the outlet zone. A variation in groundwater salinity was observed in both space and time. In 1975, salinity was very low. The outlet zone was the most affected where the drawdown reached several meters, causing saltwater intrusion from the surrounding sabkhas.     

Three-dimensional supervirtual seismic refraction interferometry

Improved supervirtual interferometry (SVI) method that enhance the signal-to-noise ratio (SNR) of noisy seismic refraction data is presented. This 3D supervirtual algorithm is helpful in removal of side lobes caused by convolution in the conventional 3D SVI method. In this study, we extend the supervirtual seismic refraction interferometry method to 3D geometries commonly used in active seismic exploration. To achieve this objective, synthetic 3D seismic refraction data were created using single patch orthogonal geometry. Simple two-layers velocity model was used. A zero phase wavelet was convolved with the refraction travel time to create the seismic wiggle trace. The refraction arrival is computed using the equation for horizontal refractor case and is not adapted for dipping refractor. Noise of Gaussian distribution with zero mean and 0.25 standard deviation was added to simulate a case of moderate ambient noise. The supervirtual algorithm, consisting of crosscorrelation, alignment, summation, and first arrival calculation is performed. By aligning and summation of all the correlogram, the stationary position of source-receiver pairs as required in the 2D supervirtual method are eliminated in this case. Synthetic data presented in this study shows accurate first arrivals after the application of the 3D SVI and traces with much better SNR than the actual traces.