New variogram modeling method using MGGP and SVR

A critical step for kriging in geostatistics is estimation of the variogram. Traditional variogram modeling comprise of the experimental variogram calculation, appropriate variogram model selection and model parameter determination. Selecting of the variogram model and fitting of model parameters is the most controversial aspect of geostatistics. Shapes of valid variogram models are finite, and sometimes, the optimal shape of the model can not be fitted, leading to reduced estimation accuracy. In this paper, a new method is presented to automatically construct a model shape and fit model parameters to experimental variograms using Support Vector Regression (SVR) and Multi-Gene Genetic Programming (MGGP). The proposed method does not require the selection of a variogram model and can directly provide the model shape and parameters of the optimal variogram. The validity of the proposed method is demonstrated in a number of cases.     

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.     

Characteristics of deuterium excess parameters for geothermal water in Beijing

This study investigates the characteristics of geothermal water in 10 geothermal fields in Beijing. The relationships between the deuterium excess parameter (d) and temperature, depth, age of geothermal groundwater, groundwater flow field, and Eh were investigated using geothermal groundwater samples. Results showed that (1) the average d value of geothermal water is 5.4, whereas that of the groundwater in normal temperature is 6.04. The differences are induced by the oxygen isotope exchange during the water–rock interaction, which may be more easily completed in geothermal water than in cold groundwater. (2) The d value increases remarkably with the age of the geothermal groundwater. The d value increases from 11.2 to 14.6 when the age of the geothermal water is 12,760 ± 130 a and 38,960 ± 630 a, respectively. Moreover, the isotope heat exchange for composition of the hydrogen and oxygen isotopes in the geothermal groundwater proceeds sufficiently with time. (3) The d value decreases from 5.72 to 3.03 when the depth increases from 125.13 to 3221 m. Generally, in the same area, the d value decreases with depth because the temperature is increasing. (4) The d value of the groundwater gradually reduces from the northern recharge area to the southern discharge area. The average d value is 7.31 in the northern recharge area and 5.68 in the middle Beijing Depression, whereas the d value in the southern area of Fengheying is ?9.20. The larger difference in d values between the recharge and discharge areas is due to the slower velocity of underwater flow, which induces longer time for oxygen exchange. (5) The relationship between the d and Eh is complex. When Eh is <200 mV, the d value of the geothermal water decreases with the decrease in Eh. When Eh is higher than 200 mV, the d value increases slightly with the decrease in Eh. The study of the characteristics of deuterium excess parameters for geothermal water could provide a scientific isotopic evidence for assessment and exploitation measures in geothermal groundwater systems.