Sharif Digital Repository

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... 

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... 

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... 

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... 

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... 

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... 

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... 

The fresh groundwater lenses (FGLs) of small islands can be highly vulnerable to climate change impacts, including sea-level rise (SLR). Many real cases of atoll or sandy islands involve two-layer hydrogeological conceptualizations. In this paper, the influential factors that affect FGLs in two-layer small islands subject to SLR are investigated. An analytical solution describing FGLs in circular islands, composed of two geological layers, is developed for the simplified case of steady-state and sharp-interface conditions. An application of the developed model is demonstrated to estimate the FGL thickness of some real-world islands by comparison with existing FGL thickness data. Furthermore,... 

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... 

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...