Application of geostatistics with Indicator Kriging for analyzing spatial variability of groundwater arsenic concentrations in Southwest Bangladesh

This article seeks to explore the spatial variability of groundwater arsenic (As) concentrations in Southwestern Bangladesh. Facts about spatial pattern of As are important to understand the complex processes of As concentrations and its spatial predictions in the unsampled areas of the study site. The relevant As data for this study were collected from Southwest Bangladesh and were analyzed with Flow Injection Hydride Generation Atomic Absorption Spectrometry (FI-HG-AAS). A geostatistical analysis with Indicator Kriging (IK) was employed to investigate the regionalized variation of As concentration. The IK prediction map shows a highly uneven spatial pattern of arsenic concentrations. The safe zones are mainly concentrated in the north, central and south part of the study area in a scattered manner, while the contamination zones are found to be concentrated in the west and northeast parts of the study area. The southwest part of the study area is contaminated with a highly irregular pattern. A Generalized Linear Model (GLM) was also used to investigate the relationship between As concentrations and aquifer depths. A negligible negative correlation between aquifer depth and arsenic concentrations was found in the study area. The fitted value with 95 % confidence interval shows a decreasing tendency of arsenic concentrations with the increase of aquifer depth. The adjusted mean smoothed lowess curve with a bandwidth of 0.8 shows an increasing trend of arsenic concentration up to a depth of 75 m, with some erratic fluctuations and regional variations at the depth between 30 m and 60 m. The borehole lithology was considered to analyze and map the pattern of As variability with aquifer depths. The study has performed an investigation of spatial pattern and variation of As concentrations.