A two-dimensional finite difference model is developed for the simulation of saltwater intrusion in wave-induced aquifer system with either a confined or phreatic aquifer. The model considers many important factors, such as the dynamic pressure induced by wave motion, the pressure wave equation, the density-dependent Darcy¡¦s Law, and the salt transport equation. This paper presents numerical study of the effect of wave motion, resulting salinity structure responses and phreatic surface fluctuation on the process of seawater intrusion ¡K etc.

Kyoto University (京都大学) / 0048 / 新制・課程博士 / 博士(工学) / 甲第14599号 / 工博第3067号 / 新制||工||1456(附属図書館) / 26951 / UT51-2009-D311 / 京都大学大学院工学研究科都市環境工学専攻 / (主査)教授 大西 有三, 教授 間瀬 肇, 准教授 西山 哲 / 学位規則第4条第1項該当

Groundwater contamination is a very serious problem as it leads to the depletion of water resources. Seawater intrusion is a special category of groundwater contamination that threatens the health and possibly lives of many people living in coastal areas. The focus of this work is to develop a numerical model to study seawater intrusion and its effects on groundwater quality and develop a control method to effectively control seawater intrusion. Two major approaches are used in this study: the first approach is the development of a finite element model to simulate seawater intrusion; the second is the development of a simulation-optimization model to study the control of seawater intrusion in coastal aquifers using different management scenarios. The simulation-optimization model is based on the integration of a genetic algorithm optimization technique with the transient density-dependent finite element model developed in this research. The finite element model considers the coupled flow of air and water and solute transport in saturated and unsaturated soils. The governing differential equations include two mass balance equations of water and air phases and the energy balance equation for heat transfer, together with a balance equation for miscible solute transport. The nonlinear governing differential equations are solved using the finite element method in the space domain and a finite difference scheme in the time domain. A two dimensional finite element model is developed to solve the governing equations and provide values of solute concentration, pore water pressure, pore air pressure and temperature at different points within the region at different times. The mathematical formulation and numerical implementation of the model are presented. The numerical model is validated by application to standard examples from literature followed by application to a number of case studies involving seawater intrusion problems. The results show good agreement with previous results reported in the literature. The model is then used to predict seawater intrusion for a number of real world case studies. The developed model is capable of predicting, with a good accuracy, the intrusion of seawater in coastal aquifers. In the second approach, a simulation-optimization model is developed to study the control of seawater intrusion using three management scenarios: abstraction of brackish water, recharge of fresh water and combination of abstraction and recharge. The objectives of these management scenarios include minimizing the total costs for construction and operation, minimizing salt concentrations in the aquifer and determining the optimal depths, locations and abstraction/recharge rates for the wells. Also, a new methodology is presented to control seawater intrusion in coastal aquifers. In the proposed methodology ADR (abstraction, desalination and recharge), seawater intrusion is controlled by abstracting brackish water, desalinating it using a small scale reverse osmosis plant and recharging to the aquifer. The simulation-optimization model is applied to a number of case studies. The efficiencies of three different scenarios are examined and compared. Results show that all the three scenarios could be effective in controlling seawater intrusion. However, ADR methodology can result in the lowest cost and salt concentration in aquifers and maximum movement of the transition zone towards the sea. The results also show that for the case studies considered in this work, the amount of abstracted and treated water is about three times the amount required for recharge; therefore the remaining treated water can be used directly for different proposes. The application of ADR methodology is shown to be more efficient and more practical, since it is a cost-effective method to control seawater intrusion in coastal aquifers. This technology can be used for sustainable development of water resources in coastal areas where it provides a new source of treated water. The developed method is regard as an effective tool to control seawater intrusion in coastal aquifers and can be applied in areas where there is a risk of seawater intrusion. Finally, the developed FE model is applied to study the effects of likely climate change and sea level rise on seawater intrusion in coastal aquifers. The results show that the developed model is capable of predicting the movement of the transition zone considering the effects of sea level rise and over-abstraction. The results also indicate that the change of water level in the sea side has a significant effect on the position of the transition zone especially if the effect of sea level rise is combined with the effect of increasing abstraction from the aquifer.

Estuarine ecosystems are under enormous stress due to rapid coastal developments and climate change. Proper management of these important ecosystems requires a good understanding of their key processes. In this thesis, riparian groundwater-surface water interaction is explored for an aquifer-estuary system primarily by a series of numerical experiments. The work focuses on riparian-scale groundwater flow and salinization. The overall aim of the study was to extend our understanding of aquifer-estuary exchange, which is currently centered on the lower marine estuarine reach, to middle estuaries (i.e., the estuary reach that has variable salinity). The numerical experiments were guided by previous studies and observations made from an exploratory field investigation conducted in and next to Sandy Creek, a macro-tidal estuary incised in the alluvial aquifer of the Pioneer Valley, North-eastern Australia (Longitude 49.11°, Latitude -21.27°). The following observations were made from the field investigation: Sandy Creek estuary experiences a variable salinity regime in its mid reaches that consists of periods of 1) freshwater flushing due to up catchment-derived flooding, 2) persistent freshwater conditions for at least 2 months following the flooding, 3) tidal salinity fluctuations and 4) constant near-seawater salinity; laterally extensive and disconnected aquitards were found to occur at the field site; Sandy Creek had an essentially ‘vertical’ bank slope. Numerical simulations were conducted using the finite element modeling code FEFLOW for saturated unsaturated, variable-density groundwater flow and solute transport, to examine the influence of the following factors on aquifer-estuary exchange: a tidally varying estuarine salinity and hydraulic head, a seasonal freshwater flush (i.e., estuary with freshwater and an elevated stage due to an up catchment sourced flood), near estuary aquitard layers, lateral asymmetry (about the estuary centerline) in hydraulic conductivity and regional hydraulic gradients. The simulations neglected seepage face development after numerical experiments showed that for a vertical bank estuary interacting with a sandy loam aquifer, seepage face effects on groundwater flow and associated salinity distribution were minimal. The following observations were drawn from the range of numerical experiments considered. Tidal salinity fluctuations in the estuary (varying between 0 and 1 - i.e., using a relative salinity scale where a salinity of 1 is seawater) produced flow paths and residence times that were distinctly different to the constant seawater salinity case. While the constant average 0.5 salinity case and the corresponding tidally-varying salinity case (i.e., salinity varying between 0 and 1) produced somewhat comparable results in terms of RUC and RLC (RUC represents groundwater discharge to the estuary that originated from recharge to the estuary bank and RLC groundwater discharge to the estuary that originated from recharge through the estuary bed), whereas flow paths and the total salt mass in the aquifer differed. Freshwater flushing simulations indicated that the near-estuary aquifer responds rapidly to a 2-day ‘wet season’ flushing event with a short-lived freshwater lens created through freshening of the hyporheic zone. Annual cycling of the seasonal flushing led to significant disruption of the estuary water circulation in the aquifer thereby impacting on residence times, transport pathways, and RUC and RLC, and acting to potentially remobilize groundwater and contaminants previously trapped in continuous and semi-continuous re-circulation cells. Although groundwater flow paths determined using tide-averaged velocity vectors were representative of flow paths from transient tidally driven flow vector field, residence times calculated from the two flow fields were markedly different. The influence of riparian scale aquitards and lateral asymmetry (about the estuary centreline) in hydraulic gradients and hydraulic conductivity on groundwater flow and associated salinity distribution was also found to be sensitive to estuarine salinity conditions. The results indicate that observations made about aquifer-estuary interaction in the lower estuary may not be directly applicable to the middle estuary. According to the simulations, tidal salinity variations in the estuary are important factors that affect hyporheic-riparian salt transport processes and that the use of a time averaged estuarine salinity as an approximation to variable salinity conditions is unsuitable for the accurate prediction of the near-estuary dynamics in middle estuaries. This study was based on a two dimensional representation of the riparian scale interaction and it is clear that future research needs to focus on the three-dimensionality of the aquifer-estuary system, incorporating spatially and temporally varying flow and transport characteristics. That is, many estuaries are tortuous and the aquifer geology spatially complex such that assumptions required for the two-dimensional section will most likely restrict application to the field. The tidal dynamics in the middle estuary is also expected to generate three dimensional aspects to the aquifer-estuary interaction. Thus further investigation that explicitly models the hydrodynamics and salt transport in the estuary and estuarine morphology is required to refine the insight provided by the simple conceptual model adopted in this study.

Excessive pumping in coastal aquifers results in seawater intrusion where optimal and efficient planning is essential. In this study, numerical solution of single potential solution by Strack is combined with genetic algorithm (GA) to find the maximum extraction amount in a coastal aquifer. Seawater intrusion is tracked with the potential value at the extraction well locations. A code is developed by combining GA and a subroutine repeatedly calling MODFLOW as a numerical solver to calculate the potential distribution for different configurations of solution (trial solutions). Potential distributions are used to evaluate the fitness values for GA. The developed model is applied to a previous work by Mantoglou. Another heuristic method, simulated annealing (SA) is utilized to compare the results of GA. Different seawater prevention methods (i.e. injection wells, canals) and decision variables related to those methods (i.e. location of the injection wells or canals) are added to model to further prevent the seawater intrusion and improve the coastal aquifer benefit. A method called &ldquo / Alternating Constraints Method&rdquo / is introduced to improve the solution for the cases with variable location. The results show that both proposed method and the regular solution with GA or SA prove to be successful methods for the optimal management of coastal aquifers.

In the era of intense pressures on water resources, the loss of groundwater by increased seawater intrusion (SWI), driven by climate, sea level and landscape changes, may be critical for many people living in commonly populous coastal regions. Analytical solutions have been derived here for interface flow in coastal aquifers, which allow for simple quantification of SWI under extended conditions from previously available such solutions and are suitable for first-order regional vulnerability assessment and mapping of the implications of climate- and landscape-driven change scenarios and related comparisons across various coastal world regions. Specifically, the derived solutions can account for the hydraulically significant aquifer bed slope in quantifying the toe location of a fresh-seawater sharp interface in the present assessments of vulnerability and safe exploitation of regional coastal groundwater.  Results show high nonlinearity of SWI responses to hydro-climatic and groundwater pumping changes on the landside and sea level rise on the marine side, implying thresholds, or tipping points, which, if crossed, may lead abruptly to major SWI of the aquifer. Critical limits of coastal groundwater change and exploitation have been identified and quantified in direct relation to prevailing local-regional conditions and stresses, defining a safe operating space for the human use of coastal groundwater. Generally, to control SWI, coastal aquifer management should focus on adequate fresh groundwater discharge to the sea, rather than on maintaining a certain hydraulic head at some aquifer location. First-order vulnerability assessments for regional Mediterranean aquifers of the Nile Delta Aquifer, the Israel Coastal Aquifer  and the Cyprus Akrotiri Aquifer show that in particular the first is seriously threatened by advancing seawater. Safe operating spaces determined for the latter two show that the current pumping schemes are not sustainable under declining recharge. / <p>The thesis was founded by two research programmes: NEO private-academic sector partnership and Ekoklim, a strategic governmental funding through Stockholm University</p><p>At the time of the doctoral defense, the following paper was unpublished and had a status as follows: Paper 3: Manuscript.</p><p> </p>

Estuarine ecosystems are under enormous stress due to rapid coastal developments and climate change. Proper management of these important ecosystems requires a good understanding of their key processes. In this thesis, riparian groundwater-surface water interaction is explored for an aquifer-estuary system primarily by a series of numerical experiments. The work focuses on riparian-scale groundwater flow and salinization. The overall aim of the study was to extend our understanding of aquifer-estuary exchange, which is currently centered on the lower marine estuarine reach, to middle estuaries (i.e., the estuary reach that has variable salinity). The numerical experiments were guided by previous studies and observations made from an exploratory field investigation conducted in and next to Sandy Creek, a macro-tidal estuary incised in the alluvial aquifer of the Pioneer Valley, North-eastern Australia (Longitude 49.11°, Latitude -21.27°). The following observations were made from the field investigation: Sandy Creek estuary experiences a variable salinity regime in its mid reaches that consists of periods of 1) freshwater flushing due to up catchment-derived flooding, 2) persistent freshwater conditions for at least 2 months following the flooding, 3) tidal salinity fluctuations and 4) constant near-seawater salinity; laterally extensive and disconnected aquitards were found to occur at the field site; Sandy Creek had an essentially ‘vertical’ bank slope. Numerical simulations were conducted using the finite element modeling code FEFLOW for saturated unsaturated, variable-density groundwater flow and solute transport, to examine the influence of the following factors on aquifer-estuary exchange: a tidally varying estuarine salinity and hydraulic head, a seasonal freshwater flush (i.e., estuary with freshwater and an elevated stage due to an up catchment sourced flood), near estuary aquitard layers, lateral asymmetry (about the estuary centerline) in hydraulic conductivity and regional hydraulic gradients. The simulations neglected seepage face development after numerical experiments showed that for a vertical bank estuary interacting with a sandy loam aquifer, seepage face effects on groundwater flow and associated salinity distribution were minimal. The following observations were drawn from the range of numerical experiments considered. Tidal salinity fluctuations in the estuary (varying between 0 and 1 - i.e., using a relative salinity scale where a salinity of 1 is seawater) produced flow paths and residence times that were distinctly different to the constant seawater salinity case. While the constant average 0.5 salinity case and the corresponding tidally-varying salinity case (i.e., salinity varying between 0 and 1) produced somewhat comparable results in terms of RUC and RLC (RUC represents groundwater discharge to the estuary that originated from recharge to the estuary bank and RLC groundwater discharge to the estuary that originated from recharge through the estuary bed), whereas flow paths and the total salt mass in the aquifer differed. Freshwater flushing simulations indicated that the near-estuary aquifer responds rapidly to a 2-day ‘wet season’ flushing event with a short-lived freshwater lens created through freshening of the hyporheic zone. Annual cycling of the seasonal flushing led to significant disruption of the estuary water circulation in the aquifer thereby impacting on residence times, transport pathways, and RUC and RLC, and acting to potentially remobilize groundwater and contaminants previously trapped in continuous and semi-continuous re-circulation cells. Although groundwater flow paths determined using tide-averaged velocity vectors were representative of flow paths from transient tidally driven flow vector field, residence times calculated from the two flow fields were markedly different. The influence of riparian scale aquitards and lateral asymmetry (about the estuary centreline) in hydraulic gradients and hydraulic conductivity on groundwater flow and associated salinity distribution was also found to be sensitive to estuarine salinity conditions. The results indicate that observations made about aquifer-estuary interaction in the lower estuary may not be directly applicable to the middle estuary. According to the simulations, tidal salinity variations in the estuary are important factors that affect hyporheic-riparian salt transport processes and that the use of a time averaged estuarine salinity as an approximation to variable salinity conditions is unsuitable for the accurate prediction of the near-estuary dynamics in middle estuaries. This study was based on a two dimensional representation of the riparian scale interaction and it is clear that future research needs to focus on the three-dimensionality of the aquifer-estuary system, incorporating spatially and temporally varying flow and transport characteristics. That is, many estuaries are tortuous and the aquifer geology spatially complex such that assumptions required for the two-dimensional section will most likely restrict application to the field. The tidal dynamics in the middle estuary is also expected to generate three dimensional aspects to the aquifer-estuary interaction. Thus further investigation that explicitly models the hydrodynamics and salt transport in the estuary and estuarine morphology is required to refine the insight provided by the simple conceptual model adopted in this study.

Estuarine ecosystems are under enormous stress due to rapid coastal developments and climate change. Proper management of these important ecosystems requires a good understanding of their key processes. In this thesis, riparian groundwater-surface water interaction is explored for an aquifer-estuary system primarily by a series of numerical experiments. The work focuses on riparian-scale groundwater flow and salinization. The overall aim of the study was to extend our understanding of aquifer-estuary exchange, which is currently centered on the lower marine estuarine reach, to middle estuaries (i.e., the estuary reach that has variable salinity). The numerical experiments were guided by previous studies and observations made from an exploratory field investigation conducted in and next to Sandy Creek, a macro-tidal estuary incised in the alluvial aquifer of the Pioneer Valley, North-eastern Australia (Longitude 49.11°, Latitude -21.27°). The following observations were made from the field investigation: Sandy Creek estuary experiences a variable salinity regime in its mid reaches that consists of periods of 1) freshwater flushing due to up catchment-derived flooding, 2) persistent freshwater conditions for at least 2 months following the flooding, 3) tidal salinity fluctuations and 4) constant near-seawater salinity; laterally extensive and disconnected aquitards were found to occur at the field site; Sandy Creek had an essentially ‘vertical’ bank slope. Numerical simulations were conducted using the finite element modeling code FEFLOW for saturated unsaturated, variable-density groundwater flow and solute transport, to examine the influence of the following factors on aquifer-estuary exchange: a tidally varying estuarine salinity and hydraulic head, a seasonal freshwater flush (i.e., estuary with freshwater and an elevated stage due to an up catchment sourced flood), near estuary aquitard layers, lateral asymmetry (about the estuary centerline) in hydraulic conductivity and regional hydraulic gradients. The simulations neglected seepage face development after numerical experiments showed that for a vertical bank estuary interacting with a sandy loam aquifer, seepage face effects on groundwater flow and associated salinity distribution were minimal. The following observations were drawn from the range of numerical experiments considered. Tidal salinity fluctuations in the estuary (varying between 0 and 1 - i.e., using a relative salinity scale where a salinity of 1 is seawater) produced flow paths and residence times that were distinctly different to the constant seawater salinity case. While the constant average 0.5 salinity case and the corresponding tidally-varying salinity case (i.e., salinity varying between 0 and 1) produced somewhat comparable results in terms of RUC and RLC (RUC represents groundwater discharge to the estuary that originated from recharge to the estuary bank and RLC groundwater discharge to the estuary that originated from recharge through the estuary bed), whereas flow paths and the total salt mass in the aquifer differed. Freshwater flushing simulations indicated that the near-estuary aquifer responds rapidly to a 2-day ‘wet season’ flushing event with a short-lived freshwater lens created through freshening of the hyporheic zone. Annual cycling of the seasonal flushing led to significant disruption of the estuary water circulation in the aquifer thereby impacting on residence times, transport pathways, and RUC and RLC, and acting to potentially remobilize groundwater and contaminants previously trapped in continuous and semi-continuous re-circulation cells. Although groundwater flow paths determined using tide-averaged velocity vectors were representative of flow paths from transient tidally driven flow vector field, residence times calculated from the two flow fields were markedly different. The influence of riparian scale aquitards and lateral asymmetry (about the estuary centreline) in hydraulic gradients and hydraulic conductivity on groundwater flow and associated salinity distribution was also found to be sensitive to estuarine salinity conditions. The results indicate that observations made about aquifer-estuary interaction in the lower estuary may not be directly applicable to the middle estuary. According to the simulations, tidal salinity variations in the estuary are important factors that affect hyporheic-riparian salt transport processes and that the use of a time averaged estuarine salinity as an approximation to variable salinity conditions is unsuitable for the accurate prediction of the near-estuary dynamics in middle estuaries. This study was based on a two dimensional representation of the riparian scale interaction and it is clear that future research needs to focus on the three-dimensionality of the aquifer-estuary system, incorporating spatially and temporally varying flow and transport characteristics. That is, many estuaries are tortuous and the aquifer geology spatially complex such that assumptions required for the two-dimensional section will most likely restrict application to the field. The tidal dynamics in the middle estuary is also expected to generate three dimensional aspects to the aquifer-estuary interaction. Thus further investigation that explicitly models the hydrodynamics and salt transport in the estuary and estuarine morphology is required to refine the insight provided by the simple conceptual model adopted in this study.

Seawater intrusion (SWI) is a widespread environmental problem, particularly in arid and semi-arid coastal areas. Unplanned prolonged over-pumping of groundwater is the most important factor in SWI that could result in severe deterioration of groundwater quality. Therefore, appropriate management strategies should be implemented in coastal aquifers to control SWI with acceptable limits of economic and environmental costs. This PhD project presents the development and application of a simulation-optimization (S/O) model to assess different management methods of controlling saltwater intrusion while satisfying water demands, and with acceptable limits of economic and environmental costs, in confined and unconfined coastal aquifers. The first S/O model (FE-GA) is developed by direct linking of an FE simulation model with a multi-objective Genetic Algorithm (GA) to optimize the efficiency of a wide range of SWI management scenarios. However, in this S/O framework, several multiple calls of the simulation model by the population-based optimization model, evaluating best individual candidate solutions resulted in a considerable computational burden. To solve this problem the numerical simulation model is replaced by an Evolutionary Polynomial Regression (EPR)-based surrogate model in the next S/O model (EPR-GA). Through these S/O approaches (FE-GA and EPR-GA) the optimal coordinates and rates of the both abstraction and recharge barriers are determined in the studied management scenarios. As a result, a new combined methodology, so far called ADRTWW, is proposed to control SWI. The ADRTWW model consists of deep Abstraction of saline water near the coast followed by Desalination of the abstracted water to a potable level for public uses and simultaneously Recharging the aquifer using a more economic source of water such as treated wastewater (TWW). In accordance to the available recharge options (injection through well or infiltration from surface pond), the general performance of ADRTWW is evaluated in different hydro-geological settings of the aquifers indicating that it offers the least cost and least salinity in comparison with other scenarios. The great capabilities of both developed S/O models in identification of the best management solutions and the optimal coordinates and rates of the abstraction well and recharge well/pond are discussed. Both FE-GA and EPR-GA can be successfully employed by a robust decision support system. In the next phase of the study, the general impacts of sea level rise (SLR), associated with its transgression nature along the coastline surface on the saltwater intrusion mechanism are investigated in different hypothetical and real case studies of coastal aquifer systems. The results show that the rate and the amount of SWI are considerably greater in aquifers with flat shoreline slopes compared with those with steep slopes. The SWI process is followed by a significant depletion in quantity of freshwater resources at the end of the century. The situation is exacerbated with combined action of SLR and groundwater withdrawals. This finding is also confirmed by 3D simulation of SWI in a regional coastal aquifer (Wadi Ham aquifer) in the UAE subjected to the coupled actions of SLR and pumping.

Brackish springs are relatively frequent phenomena in coastal carbonate formations and their existence has been extensively reported in Mediterranean coasts. In fact, more than 300 brackish springs have been identified only in the coast of the former Yugoslavia. They essentially consist of inland or submarine karst outlets discharging waters with flow-dependent salinity. The phenomenon is particularly surprising in inland springs, where high flow rates with significant salinities (presumably coming from the sea) may be discharged several meters above sea level. In addition to its scientific interest, brackish springs hold a strategic potential as a source of water in areas with often limited water resources. In order to design their appropriate management a quantitative understanding of their controlling mechanisms both in general and at every particular spring has to be achieved.These mechanisms have been studied for many years but some controversy still remains. It is clear that they are related to deep well developed karst systems. Under these conditions, groundwater flows in a turbulent mode through a network of interconnected conduits immersed in a porous matrix with slow Darcyan flow velocities. Surprisingly, different models to explain the functioning of the system, although based on different conceptual and methodological approaches lead to similar results. This sugests that a global study on the salinitzation mechanisms of brackish springs should be undertaken. Here, we first derive the equations governing turbulent flow for density-dependent fluids and describe different mechanisms of salinization of inland brackish springs, in order to compare with the spring discharge and concentration response for those mechanisms of salinization.The insights gained in this analysis are applied to the study of S'Almadrava spring (Mallorca, Spain). This spring discharges up to 2 m3/s with salinities of 20 mS/cm at an elevation of 8 m.a.s.l. It generally displays an inverse relation between discharge rate and concentration (i.e., discharging higher salinity waters for low flow rates, and vice versa). A hypothetical but geologically feasible dual permeability model is proposed to reproduce observed salinity variations for both the dry and wet seasons but also to explain the secondary salinity peaks observed after every rainfall event. Model results agree with observations, but the lack of geological information at depth impedes model validation. Therefore, a second validation of the conceptual model is undertaken based on high-frequency geochemical observations. Due to the highly dynamic conditions of the system, the geochemical data was analyzed using fully coupled reactive transport modelling. The interpretation of geochemical data not only helps on validating conceptual models but also yields information on the water-rock interaction processes occurring at deep carbonate systems. In fact, one of the processes initially proposed to explain the occurrence of well-developed karst systems at depth, is the enlargement of tectonic fissures by carbonate dissolution due to the mixing of fresh and seawater.The theory of dissolution by mixing waters is based on the fact that when two solutions are mixed, concentrations in the mixture are volume weighted averages of the two end-members, but the thermodynamic activities of the species controlling the water-mineral reactions are non linear functions of the mixing ratio. Therefore, two end-member solutions in equilibrium with a solid phase could lead to an undersaturated mixture depending on several factors, most notably CO2 content and ionic strength. Observation of mixing and carbonate dissolution at depth has not been possible because of technical difficulties. More accessible to observation is the seawater mixing zone in coastal aquifers where calcite undersaturation and/or calcite dissolution have been reported numerous times. Yet, dissolution in coastal environments is not always clear and oversaturation or lack of dissolution in mixing zones have also been described. This apparent inconsistency on field observations around the world prompted the studies of the second part of the thesis. Flow-through laboratory experiments were performed in CO2-controlled atmosphere in order to quantify the dependence of the dissolution of calcite with the mixing ratio, and the role that CO2 variations may have on enhancing the dissolution capacity of the mixture. Results show that, although dissolution occurs, the major carbonate dissolution in aquifers must be considered only in a geological time scale. Sanford and Konikow (1989) predicted the location and magnitude of long term porosity development of coastal aquifers, based on a two step method. We compare their results with a reactive transport model approach in 1D and 2D, showing that reactive transport is required to properly understand the phenomenon because it is found that dissolution is controlled not only by geochemical factors but also by the rate at which fresh and salt water mix (i.e., by dispersion).