Sharif Digital Repository

Understanding the fresh groundwater lens (FGL) behavior and potential threat of climatic-induced seawater intrusion (SWI) are significant for the future water resources management of many small islands. In this paper, the FGL of Kish Island, an arid-region case in the Persian Gulf, Iran, is modeled using two-dimensional (2D) and three-dimensional (3D) simulations. These simulations are based on the application of SUTRA, a density-dependent groundwater numerical model. Also, the numerical model parameters are calibrated using PEST, an automated parameter estimation code. Firstly a detailed conceptualization of the FGL model is completed to understand the sensitivity of the FGL to some... 

The implementation of Monte Carlo simulations (MCSs) for the propagation of uncertainty in real-world seawater intrusion (SWI) numerical models often becomes computationally prohibitive due to the large number of deterministic solves needed to achieve an acceptable level of accuracy. Previous studies have mostly relied on parallelization and grid computing to decrease the computational time of MCSs. However, another approach which has received less attention in the literature is to decrease the number of deterministic simulations by using more efficient sampling strategies. Sampling efficiency is a measure of the optimality of a sampling strategy. A more efficient sampling strategy requires... 

A number of challenges including instability, nonconvergence, nonuniqueness, nonoptimality, and lack of a general guideline for inverse modelling have limited the application of automatic calibration by generic inversion codes in solving the saltwater intrusion problem in real-world cases. A systematic parameter selection procedure for the selection of a small number of independent parameters is applied to a real case of saltwater intrusion in a small island aquifer system in the semiarid region of the Persian Gulf. The methodology aims at reducing parameter nonuniqueness and uncertainty and the time spent on inverse modelling computations. Subsequent to the automatic calibration of the... 

Real world models of seawater intrusion (SWI) require high computational efforts. This creates computational difficulties for the uncertainty propagation (UP) analysis of these models due the need for repeated numerical simulations in order to adequately capture the underlying statistics that describe the uncertainty in model outputs. Moreover, despite the obvious advantages of moment-independent global sensitivity analysis (SA) methods, these methods have rarely been employed for SWI and other complex groundwater models. The reason is that moment-independent global SA methods involve repeated UP analysis which further becomes computationally demanding. This study proposes the use of... 

Existing closed-form solutions of contaminant transport problems are limited by the mathematically convenient assumption of uniform flow. These solutions cannot be used to investigate contaminant transport in coastal aquifers where seawater intrusion induces a variable velocity field. An adaptation of the Fourier-Galerkin method is introduced to obtain semi-analytical solutions for contaminant transport in a confined coastal aquifer in which the saltwater wedge is in equilibrium with a freshwater discharge flow. Two scenarios dealing with contaminant leakage from the aquifer top surface and contaminant migration from a source at the landward boundary are considered. Robust implementation of... 

The Henry problem (HP) continues to play a useful role in theoretical and practical studies related to seawater intrusion (SWI) into coastal aquifers. The popularity of this problem is attributed to its simplicity and precision to the existence of semi-analytical (SA) solutions. The first SA solution has been developed for a high uniform diffusion coefficient. Several further studies have contributed more realistic solutions with lower diffusion coefficients or velocity-dependent dispersion. All the existing SA solutions are limited to homogenous and isotropic domains. This work attempts to improve the realism of the SA solution of the dispersive HP by extending it to heterogeneous and... 

In this study, a review of the overlay/index methods served for delineation of vulnerable zones in coastal aquifers affected by SWI is provided. Then, a more realistic presentation of the vulnerability mapping of coastal aquifers to SWI through modified GALDIT index method by incorporating the influential factors on SWI is established. The modifications on GALDIT method including incorporating the seaward hydraulic gradient (i) instead of the height of groundwater level above sea level (L) (so-called GAiDIT), and considering hydraulic gradient (i) as an additional parameter to the GALDIT (so-called GALDIT-i). Three GALDIT, GAiDIT, and GALDIT-i methods were evaluated with data from three... 

In this study, the effect of distance-dependent dispersion coefficients on density-driven flow is investigated. The linear asymptotic model, which assumes that dispersivities increase linearly with distance from the source of contamination and reach asymptotic values at a large asymptotic distance, is employed. An in-house numerical model is adapted to handle distance-dependent dispersion. The effect of asymptotic-dispersion on aquifer contamination is analyzed for two tests: (i) a seawater intrusion problem in a coastal aquifer and (ii) a leachate transport problem from a surface deposit site. Global Sensitivity Analysis (GSA) combined with the Polynomial Chaos Expansion (PCE) surrogate... 

Lake Urmia, which was once the second-largest saline lake in the world, has been shrinking dramatically. Moreover, Lake Urmia has become supersaturated with total salinity averaging more than 350 g/l. Salt precipitation and dissolved materials brought by inflowing rivers have formed a layer of sediment with low hydraulic conductivity on the lakebed. Considering the flat bathymetry of Lake Urmia, we conducted a series of numerical simulation scenarios to study the groundwater flow pattern in the vicinity of the hypersaline Lake Urmia using COMSOL Multiphysics®. In the first step, we performed the simulations in steady-state conditions. Secondly, we simulated the lake level fall in 10 years at... 

The interaction between a lake and groundwater has important implications to the quantity and quality of water in both environments. Quantification of lake-groundwater interaction (LGI) has been challenging in regions with limited in-situ data. LGI can be quantified by physically-based models, direct measurement of seepage, measurements of conservative chemical or isotopic tracers, and lake water balance. Despite the accuracy of the methods based on hydrochemical or isotopic measurements and analysis, they require extensive field data that are costly to collect in large lakes. Instead, the data required to quantify LGI by the lake water budget method can be obtained via typical ground... 

This paper presents an evolutionary based approach to achieve optimal management of a coastal aquifer to control saltwater intrusion. An improved Elitist Continuous Ant Colony Optimization (ECACO) algorithm is employed for optimal control variables setting of coastal aquifer management problem. The objectives of the optimal management are; maximizing the total water-pumping rate, while controlling the drawdown limits and protecting the wells from saltwater intrusion. Since present work is one of the first efforts towards the application of an ECACO algorithm, sharp interface solution for steady state problem is first exploited. The performance of the developed optimization model is evaluated... 

Sea-level rise (SLR) influences groundwater hydraulics and in particular seawater intrusion (SWI) in many coastal aquifers. The quantification of the combined and relative impacts of influential factors on SWI has not previously been considered in coastal aquifers. In the present study, a systematic review of the available literature on this topic is first provided. Then, the potential remaining challenges are scrutinized. Open questions on the effects of more realistic complexities such as gradual SLR, parameter uncertainties, and the associated influences in decision-making models are issues requiring further investigation.We assess and quantify the seawater toe location under the impacts... 

In this study we use polynomial chaos expansion (PCE) to perform uncertainty analysis for seawater intrusion (SWI) in fractured coastal aquifers (FCAs) which is simulated using the coupled discrete fracture network (DFN) and variable-density flow (VDF) models. The DFN-VDF model requires detailed discontinuous analysis of the fractures. In real field applications, these characteristics are usually uncertain which may have a major effect on the predictive capability of the model. Thus, we perform global sensitivity analysis (GSA) to provide a preliminary assessment on how these uncertainties can affect the model outputs. As our conceptual model, we consider fractured configurations of the...