Groundwater arsenic contamination and its health effects in India

During a 28-year field survey in India (1988–2016), groundwater arsenic contamination and its health effects were registered in the states of West Bengal, Jharkhand, Bihar and Uttar Pradesh in the Ganga River flood plain, and the states of Assam and Manipur in the flood plain of Brahamaputra and Imphal rivers. Groundwater of Rajnandgaon village in Chhattisgarh state, which is not in a flood plain, is also arsenic contaminated. More than 170,000 tubewell water samples from the affected states were analyzed and half of the samples had arsenic >10 μg/L (maximum concentration 3,700 μg/L). Chronic exposure to arsenic through drinking water causes various health problems, like dermal, neurological, reproductive and pregnancy effects, cardiovascular effects, diabetes mellitus, diseases of the respiratory and gastrointestinal systems, and cancers, typically involving the skin, lungs, liver, bladder, etc. About 4.5% of the 8,000 children from arsenic-affected villages of affected states were registered with mild to moderate arsenical skin lesions. In the preliminary survey, more than 10,000 patients were registered with different types of arsenic-related signs and symptoms, out of more than 100,000 people screened from affected states. Elevated levels of arsenic were also found in biological samples (urine, hair, nails) of the people living in affected states. The study reveals that the population who had severe arsenical skin lesions may suffer from multiple Bowens/cancers in the long term. Some unusual symptoms, such as burning sensation, skin itching and watering of eyes in the presence of sun light, were also noticed in arsenicosis patients.     

Bearing capacity of sandy soil treated by Kenaf fibre geotextile

Bio-based materials are widely used recently in order to introduce a more sustainable construction material. Kenaf is a type of bio-based material that can be easily obtained in a tropical country, which could be a potential material to be utilised as a geotextile material because it has good tensile strength. The geotextile could be used to improve the bearing capacity of a loose soil. This paper presents a series of small-scale physical modelling tests to investigate the bearing capacity performance of Kenaf fibre geotextile laid on and inside the sand layer. A rigid footing was used to replicate a strip footing during the loading test, and sand was prepared based on 50% of relative density in a rigid testing chamber for ground model preparation. In order to treat the soil, Kenaf fibre geotextile was laid at four difference locations which are on the soil surface and underneath the ground model surface at 50, 75 and 100 mm deep. It was found that the usage of the Kenaf fibre geotextile has improved the bearing capacity of the sandy soil up to 414.9% as compared to untreated soil. It was also found that the depth of the Kenaf fibre geotextile treated into the soil also affects the soil performance.     

Delineation of groundwater potential zones using remote sensing,GIS, and AHP technique in Tehran–Karaj plain,Iran

Evaluation of groundwater resources in dry areas without enough data is a challenging task in many parts of the world, including Tehran–Karaj plain in Iran, which includes Tehran, the capital city of Iran and Karaj, one of Iran’s biggest cities. Water demand due to increasing agricultural and industrial activities caused many problems in the field of water resources management. In this study, the potential of groundwater resources was evaluated using remote sensing, geographic information system (GIS), and analytic hierarchy process (AHP) for the first time. Digital Elevation Model from Shuttle Radar Topography Mission was used to generate a slope map and drainage density map. Three Landsat-8 satellite images were utilized to provide lineament density and land cover/land use maps. Geological and soil type maps were provided from the Geological Survey and Mineral Explorations of Iran (GSI). Tropical Rainfall Measuring Mission data were used to prepare average annual precipitation map. Discharge values from 102 pumping wells in the time period of 2002–2014 were used to evaluate the results. Seven data layers were prepared, and the geodatabase was made in GIS. The layers and their classes were assigned weights using AHP method. Finally, the layers were overlaid based on their weights, and the potential map of groundwater resources was generated. The area was classified into five zones with very high, high, moderate, low, and very low potentials. The zones covered 5.95, 32.90, 22.70, 10.20, and 28.25% of the study area, respectively. The results showed good agreement with the field data obtained from discharge wells.     

Comparing different statistical models for assessing Fe-contaminated soils based on VNIR/SWIR spectral data

Two statistical models including partial least squares regression (PLSR) and principal component regression were comparatively utilized to determine the predictive accuracy of visible–near-infrared and short-wave infrared reflectance spectroscopy in quantifying the Fe concentration in contaminated soils. Two scenarios were applied to select the best model: Scenario I included all wavelengths (400–2450 nm) and Scenario II encompassed characteristic bands of Fe. Pre-processing techniques used to select the best model included: first and second derivatives (FD and SD), multiplicative scatter correction (MSC) and standard normal variate. The abilities of the predictive models were evaluated by splitting soil samples into two random groups (80 and 20%). The first group (80%) was used to evaluate calibration and validation sets by employing the cross‐validation method, and the second group (20%) was applied to test the models. The coefficient of determination (R ²), root mean square error and residual prediction deviation were calculated to evaluate the models. Applying Scenario I indicated that the PLSR model with SD pre-processing was a more accurate technique for predicting the Fe concentration, whereas in the Scenario II, the PLSR model with MSC pre-processing had a better performance. Comparing Scenarios I and II indicated that the more reliable models for predicting the soil Fe content could be constructed by the PLSR model with the SD pre-processing techniques and all wavelengths. The modeling results produced by the PLSR model with the SD pre-processing could be used to detect, map and monitor Fe-contaminated soils by proximal and remote sensing in the mining areas.     

Modeling and mapping of solar radiation using geostatistical analysis methods in Iran

Because of economic and technical limitations, measuring solar energy received at ground level (R _s ) isn’t possible in all parts of the country, and in only 12% of synoptic stations is this parameter measured and recorded. Thus, it should be estimated and modeled spatially based on other climatic variables using mathematical methods. In this research, many attempts have been made to introduce an air temperature-based model for Rs estimation, and then, based on the output of the mentioned models, several geostatistical methods have been tested, and finally an elegant spatial model is proposed for (Rs) zoning in Iran. In this regard, the relationships between the measured amounts of monthly solar radiation and other climatic parameters, such as a monthly average, maximum and minimum temperature, precipitation, relative humidity, and the number of sunny hours during the period 1970–2010, are examined and modeled. It was revealed that based on the linear relationship between the monthly average air temperatures and solar radiation values recorded in each of the stations, that the best-fit linear model, with R ²  = 0.822, MAE = 1.81, RMSE = 2.51%, and MAPE = 10.08, can be introduced for Rs estimation. Then, using the outputs of the proposed model, the amounts of (R _s ) are estimated in another 171 meteorological stations (a total of 192 stations), and eight geostatistical methods (IDW, GPI, RBF, LPI, OK, SK, UK, and EBK) were investigated for zoning. Comparing the resulting variograms showed that in addition to proof of spatial correlation between solar radiation data, they can be applied for modeling changes in various directions. Analyzing the ratio of the nugget effect on the roof of the variograms showed that the Gaussian model with the lowest ratio (Co/Co + C = 0.883) and (R ²  = 0.972), could model the highest correlation between the data and, therefore, it was used for data interpolation. To select the best geostatistical model, R2, MAE, and RMSE were used. On this basis, it was found that the RBF method with R ²  = 0.904, MAE = 3.02, RMSE = 0.39% is the most effective. Also, the IDW method with R ²  = 0.90, MAE = 3.08, RMSE = 0.391%, compared to other methods is the most effective. In addition, for data validation, correlations between observed and estimated values of solar radiation were studied and found R ²  = 0.86.     

Sorptive Potential of Beach Sand to Remove Ni(II) Ions: An Equilibrium Isotherm Study

The potential to remove Ni(II) ions from aqueous solutions using sea beach sand, a carbonate‐quartz mineral, was thoroughly investigated. The effects of relevant parameters such as solution pH, adsorbent dose, metal ions concentration, and temperature on Ni(II) sorption onto beach sand were examined. The sorption data followed the Langmuir, Freundlich and Dubinin‐Radushkevich (D‐R) isotherms. The adsorption was endothermic in nature at ambient temperature and the computation of the parameters, ΔH, ΔS and ΔG indicated the interactions between sorbate and sorbent to be thermodynamically favorable. Equilibrium was achieved very quickly within 30 min of shaking. A pseudo‐first order Lagergren equation was used to test the adsorption kinetics. Other kinetic models, e. g., the Morris‐Weber and Reichenberg equations, were used to calculate the rate constant of intraparticle diffusion and the fate of the diffusion process, respectively. The influence of some of the common cations and anions were also a subject of this study.     

GPS observations at quasi-conjugate points under disturbed conditions

The conjugacy effects of the GPS scintillation activities during the geomagnetic storms of October 2003, November 2003 and July 2004 have been investigated at the approximately geomagnetically conjugate stations: Scott Base, Antarctica (SBA) and Resolute Cornwallis Island (RESO) in the high Arctic region. The measurements aim at investigation of the similarities and differences of the scintillation activities occurring at the conjugate points in the polar regions under storm conditions and examine the relationship between the Storm Enhanced Density (SED) and scintillation activity. The measurements of the scintillation activities obtained from total scintillation index during these storm periods at both hemispheres showed asymmetry in the ionospheric scintillation occurrence at the conjugate points. Pronounced scintillation activity was observed at the nightside hemisphere with the total scintillation index higher than at the dayside hemisphere. The results also show that the durations of severe scintillation activity were longer at the nightside hemisphere. The measurements showed that the intense scintillation periods were corresponding to the presence of the SED events where more pronounced SED events were observed over the nightside hemisphere.     

Compositional characteristics and geodynamic significance of late Miocene volcanic rocks associated with the Chah Zard epithermal gold–silver deposit,southwest Yazd,Iran

Whole‐rock geochemical and Sr–Nd isotopic data are presented for late Miocene volcanic rocks associated with the Chah Zard epithermal Au–Ag deposit in the Urumieh‐Dokhtar Magmatic Arc (UDMA), Iran, to investigate the magma source, petrogenesis and the geodynamic evolution of the study area. The Chah Zard andesitic to rhyolitic volcanic rocks are characterized by significant Large Ion Lithophile Element (LILE) and Light Rare Earth Element (LREE) enrichment coupled with High Field Strength Element (HFSE) depletion. Our geochemical data indicate an adakitic‐like signature for the volcanic rocks (e.g. SiO₂ > 62 wt%, Al₂O₃ > 15 wt%, MgO < 1.5 wt%, Sr/Y > 70, La/Yb > 35, Yb < 1 ppm, and Y < 18 ppm, and no significant Eu anomalies), distinguishing them from the other volcanic rocks of the UDMA. The Chah Zard volcanic rocks have similar Sr and Nd isotopic compositions; the ⁸⁷Sr/⁸⁶Sr(i) ratios range from 0.704 902 to 0.705 093 and the ε_Nd(i) values are from +2.33 to +2.70. However, the rhyolite porphyry represents the final stage of magmatism in the area and has a relatively high ⁸⁷Sr/⁸⁶Sr ratio (0.705 811). Our data suggest that the andesitic magmas are from a heterogeneous source and likely to result from partial melting of a metasomatized mantle wedge associated with a mixture of subducted oceanic crust and sediment. These melts subsequently underwent fractional crystallization along with minor amounts of crustal assimilation. Our study is consistent with the model that the volcanic host rocks to epithermal gold mineralization in the UDMA are genetically related to late Miocene Neo‐Tethyan slab break‐off beneath Central Iran.     

Fault slip and identification of the second fault plane in the Varzeghan earthquake doublet

An intraplate earthquake doublet, with 11-min delay between the events, devastated the city of Varzeghan in northwestern Iran on August 11, 2012. The first Mw 6.5 strike-slip earthquake, which occurred after more than 200 years of low seismicity, was followed by an Mw 6.4 oblique thrust event at an epicentral separation of about 6 km. While the first event can be associated with a distinct surface rupture, the absence of a surface fault trace and no clear aftershock signature makes it challenging to identify the fault plane of the second event. We use teleseismic body wave inversion to deduce the slip distribution in the first event. Using both P and SH waves stabilize the inversion and we further constrain the result with the surface rupture extent and the aftershock distribution. The obtained slip pattern shows two distinct slip patches with dissimilar slip directions where aftershocks avoid high-slip areas. Using the estimated slip for the first event, we calculate the induced Coulomb stress change on the nodal planes of the second event and find a preference for higher Coulomb stress on the N-S nodal plane. Assuming a simple slip model for the second event, we estimate the combined Coulomb stress changes from the two events on the focal planes of the largest aftershocks. We find that 90% of the aftershocks show increased Coulomb stress on one of their nodal planes when the N-S plane of the second event is assumed to be the correct fault plane.     

Numerical simulation of sedimentation process in reservoirs and development of a non-coupled algorithm to improve long-term modeling

This paper investigates the validity of a quasi-steady approximation for sediment transport and presents a new algorithm based on this concept. The developed non-coupled algorithm interacts among hydrodynamic, sediment, and morphology modules which are based on depth-averaged Navier-Stokes equations for the flow, the three-dimensional equation of conservation of sediment, and the mass balance between the bed and sediment (Exner equation) to simulate the reservoir sedimentation process. The non-coupled algorithm solves both the short-term scale and the relatively long-term scale problems of reservoir sedimentation. The proposed algorithm is verified using field data and by comparison with other accurate algorithms. Based upon the results of this investigation, the developed algorithm can be used to simulate long-term reservoir sedimentation while considerably decreasing the computational costs and preserving computational accuracy. The computational cost of the non-coupled algorithm is about 97% less than the conventional semi-coupled approach whereas the errors (Root Mean Square Error, Average Relative Error, and Maximum Relative Error of bed level) of the developed algorithm are approximately 15% greater than those for the semi-coupled algorithm for the average value.