Assessing the pollution risk of a groundwater source field at western Laizhou Bay under seawater intrusion.

Coastal areas have great significance for human living, economy and society development in the world. With the rapid increase of pressures from human activities and climate change, the safety of groundwater resource is under the threat of seawater intrusion in coastal areas. The area of Laizhou Bay is one of the most serious seawater intruded areas in China, since seawater intrusion phenomenon was firstly recognized in the middle of 1970s. This study assessed the pollution risk of a groundwater source filed of western Laizhou Bay area by inferring the probability distribution of groundwater Cl(-) concentration. The numerical model of seawater intrusion process is built by using SEAWAT4. The parameter uncertainty of this model is evaluated by Markov Chain Monte Carlo (MCMC) simulation, and DREAM(ZS) is used as sampling algorithm. Then, the predictive distribution of Cl(-) concentration at groundwater source field is inferred by using the samples of model parameters obtained from MCMC. After that, the pollution risk of groundwater source filed is assessed by the predictive quantiles of Cl(-) concentration. The results of model calibration and verification demonstrate that the DREAM(ZS) based MCMC is efficient and reliable to estimate model parameters under current observation. Under the condition of 95% confidence level, the groundwater source point will not be polluted by seawater intrusion in future five years (2015-2019). In addition, the 2.5% and 97.5% predictive quantiles show that the Cl(-) concentration of groundwater source field always vary between 175mg/l and 200mg/l.

[1]  Jasper A. Vrugt,et al.  High‐dimensional posterior exploration of hydrologic models using multiple‐try DREAM(ZS) and high‐performance computing , 2012 .

[2]  Alyssa M. Dausman,et al.  SEAWAT Version 4: A Computer Program for Simulation of Multi-Species Solute and Heat Transport , 2008 .

[3]  D. Rubin,et al.  Inference from Iterative Simulation Using Multiple Sequences , 1992 .

[4]  D. Higdon,et al.  Accelerating Markov Chain Monte Carlo Simulation by Differential Evolution with Self-Adaptive Randomized Subspace Sampling , 2009 .

[5]  Alain Dassargues,et al.  Conceptual model uncertainty in groundwater modeling: Combining generalized likelihood uncertainty estimation and Bayesian model averaging , 2008 .

[6]  H. Haario,et al.  An adaptive Metropolis algorithm , 2001 .

[7]  H. Yazicigil,et al.  Assessing the impacts of climate change on sustainable management of coastal aquifers , 2014, Environmental Earth Sciences.

[8]  D. Reed.  Understanding the Effects of Sea-Level Rise on Coastal Wetlands: The Human Dimension , 2010 .

[9]  C. Simmons,et al.  Sea‐level rise impact on fresh groundwater lenses in two‐layer small islands , 2014 .

[10]  J. Vrugt,et al.  A formal likelihood function for parameter and predictive inference of hydrologic models with correlated, heteroscedastic, and non‐Gaussian errors , 2010 .

[11]  T. Holford,et al.  Estimating daily nitrogen dioxide level: application of a longitudinal model with spatially-correlated random effects , 2015, Environmental and Ecological Statistics.

[12]  S. Qian,et al.  A Hierarchical Model for Estimating Long‐Term Trend of Atrazine Concentration in the Surface Water of the Contiguous U.S. , 2015 .

[13]  Alberto Guadagnini,et al.  Probabilistic assessment of seawater intrusion under multiple sources of uncertainty , 2015 .

[14]  G. C. Tiao,et al.  Bayesian inference in statistical analysis , 1973 .

[15]  Dongmei Han,et al.  Chemical and isotopic constraints on evolution of groundwater salinization in the coastal plain aquifer of Laizhou Bay, China , 2014 .

[16]  J. Jiao,et al.  Seawater intrusion and coastal aquifer management in China: a review , 2014, Environmental Earth Sciences.

[17]  Timothy R. Ginn,et al.  Impact of numerical artifact of the forward model in the inverse solution of density‐dependent flow problem , 2014 .

[18]  Keith Beven,et al.  The future of distributed models: model calibration and uncertainty prediction. , 1992 .

[19]  Ming Ye,et al.  Assessment of parametric uncertainty for groundwater reactive transport modeling , 2014 .

[20]  F. Celico,et al.  Seawater intrusion in the Guanahacabibes Peninsula (Pinar del Rio Province, western Cuba): effects on karst development and water isotope composition , 2015, Environmental Earth Sciences.

[21]  L. Aquilina,et al.  Origins and processes of groundwater salinization in the urban coastal aquifers of Recife (Pernambuco, Brazil): A multi-isotope approach. , 2015, The Science of the total environment.

[22]  Stephen P. Brooks,et al.  Markov chain Monte Carlo method and its application , 1998 .

[23]  Dong Wang,et al.  The development and control of the seawater intrusion in the eastern coastal of Laizhou Bay, China , 2008 .

[24]  Keith Beven,et al.  Equifinality, data assimilation, and uncertainty estimation in mechanistic modelling of complex environmental systems using the GLUE methodology , 2001 .

[25]  D. A. Barry,et al.  Influence of tides and waves on the fate of nutrients in a nearshore aquifer: Numerical simulations , 2014 .

[26]  Jasper A. Vrugt,et al.  UvA-DARE ( Digital Academic Repository ) DREAM ( D ) : An adaptive Markov chain Monte Carlo simulation algorithm to solve discrete , noncontinuous , and combinatorial posterior parameter estimation problems , 2011 .

[27]  B. Anderson,et al.  The rising tide: assessing the risks of climate change and human settlements in low elevation coastal zones , 2007 .

[28]  Hans Kromhout,et al.  Impact of input data uncertainty on environmental exposure assessment models: A case study for electromagnetic field modelling from mobile phone base stations. , 2014, Environmental research.

[29]  Tyler Smith,et al.  Development of a formal likelihood function for improved Bayesian inference of ephemeral catchments , 2010 .

[30]  W. Danskin,et al.  A Geochemical Approach to Determine Sources and Movement of Saline Groundwater in a Coastal Aquifer , 2014, Ground water.

[31]  Heikki Haario,et al.  DRAM: Efficient adaptive MCMC , 2006, Stat. Comput..

[32]  Hailong Li,et al.  Numerical simulations of steady-state salinity distribution and submarine groundwater discharges in homogeneous anisotropic coastal aquifers , 2014 .

[33]  Lachun Wang,et al.  Uncertainty Evaluation of a Groundwater Conceptual Model by Using a Multimodel Averaging Method , 2015 .

[34]  O. Batelaan,et al.  Three-dimensional hydrostratigraphical modelling to support evaluation of recharge and saltwater intrusion in a coastal groundwater system in Vietnam , 2014, Hydrogeology Journal.

[35]  Wenchao Sun,et al.  Estimating daily time series of streamflow using hydrological model calibrated based on satellite observations of river water surface width: Toward real world applications. , 2015, Environmental research.

[36]  C. Doulgeris,et al.  3D Variable Density Flow Simulation to Evaluate Pumping Schemes in Coastal Aquifers , 2014, Water Resources Management.