Spatial and Temporal Changes in Groundwater Salinity in South Florida

A combination of statistical and interpolation methods and Geographic Information System (GIS) spatial analysis was used to evaluate the spatial and temporal changes in groundwater Cl− concentrations in Collier and Lee Counties (southwestern Florida), and Miami-Dade and Broward Counties (southeastern Florida), since 1985. In southwestern Florida, the average Cl− concentrations in the shallow wells (0–43 m) in Collier and Lee Counties increased from 132 mg L−1 in 1985 to 230 mg L−1 in 2000. The average Cl− concentrations in the deep wells (>43 m) of southwestern Florida increased from 392 mg L−1 in 1985 to 447 mg L−1 in 2000. Results also indicated a positive correlation between the mean sea level and Cl− concentrations and between the mean sea level and groundwater levels for the shallow wells. Concentrations in the Biscayne Aquifer (southeastern Florida) were significantly higher than those of southwestern Florida. The average Cl− concentrations increased from 159 mg L−1 in 1985 to 470 mg L−1 in 2010 for the shallow wells (<33 m) and from 1360 mg L−1 in 1985 to 2050 mg L−1 in 2010 for the deep wells (>33 m). In the Biscayne Aquifer, wells showed a positive or negative correlation between mean sea level and Cl− concentrations according to their location with respect to the saltwater intrusion line. Wells located inland behind canal control structures and west of the saltwater intrusion line showed negative correlation values, whereas wells located east of the saltwater intrusion line showed positive values. Overall, the results indicated that since 1985, there was a potential decline in the available freshwater resources estimated at about 12–17% of the available drinking-quality groundwater of the southeastern study area located in the Biscayne Aquifer.

The Influence of Hydrologic Restoration on Groundwater-Surface Water Interactions in a Karst Wetland, the Everglades (FL, USA)

Efforts to rehydrate and restore surface water flow in karst wetlands can have unintended consequences, as these highly conductive and heterogeneous aquifers create a close connection between groundwater and surface water. Recently, hydrologic restoration efforts in the karstic Taylor Slough portion of the Everglades has changed from point source delivery of canal water (direct restoration), to the use of a series of surface water recharge retention basins (diffuse restoration). To determine the influence of restoration on groundwater-surface water interactions in the Taylor Slough headwaters, a water budget was constructed for 1997–2011 using 70 hydro-meteorological stations. With diffuse restoration, groundwater seepage from the Everglades toward the urban boundary increased, while the downstream delivery of surface water to the main portion of the slough declined. The combined influence of diffuse restoration and climate led to increased intra-annual variability in the volume of groundwater and surface water in storage but supported a more seasonally hydrated wetland compared to the earlier direct tactics. The data further indicated that hydrologic engineering in karst wetland landscapes enhances groundwater-surface water interactions, even those designed for restoration purposes.

Essays on intergenerational and regional aspects of water management

This dissertation consists of three essays on different aspects of water management. The first essay focuses on the sustainability of freshwater use by introducing the notion that altruistic parents do bequeath economic assets for their offspring. Constructing a two-period, over-lapping generational model, an optimal ratio of consumption and pollution for old and young generations in each period is determined. Optimal levels of water consumption and pollution change according to different parameters, such as, altruistic degree, natural recharge rate, and population growth. The second essay concerns water sharing between countries in the case of trans-boundary river basins. The paper recognizes that side payments fail to forge water-sharing agreement among the international community and that downstream countries have weak bargaining power. An interconnected game approach is developed by linking the water allocation issue with other non-water issues such as trade or border security problems, creating symmetry between countries in bargaining power. An interconnected game forces two countries to at least partially cooperate under some circumstances. The third essay introduces the concept of virtual water (VW) into a traditional international trade model in order to estimate water savings for a water scarce country. A two country, two products and two factors trade model is developed, which includes not only consumers and producer's surplus, but also environmental externality of water use. The model shows that VW trade saves water and increases global and local welfare. This study should help policy makers to design appropriate subsidy or tax policy to promote water savings especially in water scarce countries.^

Lattice Boltzmann Modeling of Fluid Flow and Solute Transport in Karst Aquifers

A novel modeling approach is applied to karst hydrology. Long-standing problems in karst hydrology and solute transport are addressed using Lattice Boltzmann methods (LBMs). These methods contrast with other modeling approaches that have been applied to karst hydrology. The motivation of this dissertation is to develop new computational models for solving ground water hydraulics and transport problems in karst aquifers, which are widespread around the globe. This research tests the viability of the LBM as a robust alternative numerical technique for solving large-scale hydrological problems. The LB models applied in this research are briefly reviewed and there is a discussion of implementation issues. The dissertation focuses on testing the LB models. The LBM is tested for two different types of inlet boundary conditions for solute transport in finite and effectively semi-infinite domains. The LBM solutions are verified against analytical solutions. Zero-diffusion transport and Taylor dispersion in slits are also simulated and compared against analytical solutions. These results demonstrate the LBM’s flexibility as a solute transport solver. The LBM is applied to simulate solute transport and fluid flow in porous media traversed by larger conduits. A LBM-based macroscopic flow solver (Darcy’s law-based) is linked with an anisotropic dispersion solver. Spatial breakthrough curves in one and two dimensions are fitted against the available analytical solutions. This provides a steady flow model with capabilities routinely found in ground water flow and transport models (e.g., the combination of MODFLOW and MT3D). However the new LBM-based model retains the ability to solve inertial flows that are characteristic of karst aquifer conduits. Transient flows in a confined aquifer are solved using two different LBM approaches. The analogy between Fick’s second law (diffusion equation) and the transient ground water flow equation is used to solve the transient head distribution. An altered-velocity flow solver with source/sink term is applied to simulate a drawdown curve. Hydraulic parameters like transmissivity and storage coefficient are linked with LB parameters. These capabilities complete the LBM’s effective treatment of the types of processes that are simulated by standard ground water models. The LB model is verified against field data for drawdown in a confined aquifer.

Evolução espaço-temporal da ocupação urbana sobre áreas naturais no município de Natal-RN

In this work the landscape morphodynamics was used to check the strength and importance of the changes carried out by man on the environment over time, in Natal-RN municipality. The occupation of partially preserved natural areas was analyzed, but environmentally fragile, such as riparian forests, vegetation on the banks of waterways, which play regulatory role of the water flow, and the dunes, which guarantee the rapid recharge of aquifers. The impacts of urban sprawl in Natal Southern and West zones Were identified and characterized, through a detailed mapping in the period between 1969 and 2013 the main Permanent Preservation Areas - PPA (banks of rivers and lagoons, and dunes remaining) and their temporal changes. For this were used aerial photographs and satellite imagery, altimetry data, and pre-existing information, which allowed the creation of a spatial database, and evolution of maps of impervious areas, evolution of the use and occupation and Digital Terrain Model (DTM) from contour lines with contour interval of 1 meter. Based on this study presents a diagnosis of the environmental situation and the state of conservation of natural areas, over the last 44 years, compared to human pressures. In general, it was found that the urban settlement has advanced about 60% of studied natural areas. This advance was growing by the year 2006, when there was a slowdown in the process, except for the Environmental Protection Zone (EPZ) 03, where the river Pitimbú and your PPA, which experienced a more significant loss area. The urban occupation affected the natural drainage and contributed to the contamination of groundwater Natal, due to increased sealed area, the release of liquid and solid waste, as well as the removal of riparian vegetation. Changed irreversibly the natural landscape, and reduced the quality and quantity of water resources necessary for the population. Thus, it is necessary to stimulate the adoption of use and protection of PPA planning measures, to the preservation of the San Valley Region inserted into the EPZ 01, and integrate more remaining dunes, in good condition, this EPZ, due to the importance of those remaining on the environment and the maintenance of quality of life. It is suggested, also, protection of catchment areas, such as PPA ponds and Pitimbú River. Finally, it is expected that this study can assist the managers in making decisions in urban and environmental planning of the municipality

Análise computacional de imagens de lâminas delgadas na predição dos parâmetros hidráulicos porosidade e condutividade hidráulica do aquífero Barreiras-RN

The understanding of the occurrence and flow of groundwater in the subsurface is of fundamental importance in the exploitation of water, just like knowledge of all associated hydrogeological context. These factors are primarily controlled by geometry of a certain pore system, given the nature of sedimentary aquifers. Thus, the microstructural characterization, as the interconnectivity of the system, it is essential to know the macro properties porosity and permeability of reservoir rock, in which can be done on a statistical characterization by twodimensional analysis. The latter is being held on a computing platform, using image thin sections of reservoir rock, allowing the prediction of the properties effective porosity and hydraulic conductivity. For Barreiras Aquifer to obtain such parameters derived primarily from the interpretation of tests of aquifers, a practice that usually involves a fairly complex logistics in terms of equipment and personnel required in addition to high cost of operation. Thus, the analysis and digital image processing is presented as an alternative tool for the characterization of hydraulic parameters, showing up as a practical and inexpensive method. This methodology is based on a flowchart work involving sampling, preparation of thin sections and their respective images, segmentation and geometric characterization, three-dimensional reconstruction and flow simulation. In this research, computational image analysis of thin sections of rocks has shown that aquifer storage coefficients ranging from 0,035 to 0,12 with an average of 0,076, while its hydrogeological substrate (associated with the top of the carbonate sequence outcropping not region) presents effective porosities of the order of 2%. For the transport regime, it is evidenced that the methodology presents results below of those found in the bibliographic data relating to hydraulic conductivity, mean values of 1,04 x10-6 m/s, with fluctuations between 2,94 x10-6 m/s and 3,61x10-8 m/s, probably due to the larger scale study and the heterogeneity of the medium studied.

Ocorrência de compostos de interesse emergente no aquífero Dunas-Barreiras e nos esgotos de Natal/RN

The detection of emerging interest microcontaminants in environmental samples of surface water, groundwater, drinking water, wastewater and effluents from water and sewage treatment plants (WTP and STP), in many countries, suggests these pollutants are widespread in the environment, mainly in urban areas. This is a reason for great concern, since many of these compounds are potentially harmful for humans other living beings, and they are not efficiently removed in the majority of WTP and STP, which is exacerbated by precariousness of water supply and sanitation services. In Natal, like other Brazilian cities, the sewage system serves only part of the urban area (about 30%), so that the rest of the wastewater is infiltrated in the sandy soil of the region in cesspool-dry well systems. This has resulted in contamination of groundwater in the area (sand-dune barrier aquifer, which supplies more than 50% of the city population), which has been observed by the increase in nitrate concentration in supply wells. The vulnerability of the sanddune barrier aquifer, combined with reports of the presence of emerging interest microcontaminants in Brazil and worldwide, led to this research, which investigated the occurrence of fifteen microcontaminants in Natal groundwater and sewage. Samples were collected at five wells used for water supply, the raw sewage and the effluents from biological reactors from STP (UASB and activated sludge reactors). Two samples of each sample were taken, with one week apart between the samples. To determine the contaminants, extraction of aquifer water, and raw and treated sewage samples were performed, through the technique of using SPE Strata X cartridge (Phenomenex®) to the aquifer water, and Strata SAX and Strata X (Phenomenex® ) for samples of raw and treated sewage. Subsequently the extracts were analyzed using GC-MS technique. Much of the analyzed microcontaminants were detected in groundwater and sewage. The concentrations in groundwater are generally lower than those found in the sewers. Some of the compounds (estrone, estradiol, bisphenol A, caffeine, diclofenac, naproxen, paracetamol and ibuprofen) are partially removed at STP.

Análise do grau de proteção do Aquífero Barreiras em perímetros sob fertilização química a partir de dados hidrogeofísicos - área da bacia do Rio Catu-RN

The use of chemical fertilization in arable perimeters provides increased productivity, though it can eventually lead to a qualitative depreciation of groundwater sources, especially if such sources are unconfined in nature. In this context, this thesis presents results from an analysis of the level of natural protection of the Barreiras Aquifer in an area located on the eastern coast of the Rio Grande do Norte State - Brazil. Such an aquifer is clastic in nature and has an unconfined hydraulic character, which clearly makes it susceptible to contamination from surface ground loads with contaminants associated with the leaching of excess fertilizers not absorbed by ground vegetation. The methodology used was based on the use of hydro-geophysical data, particularly inverse models of vertical electrical soundings (VES) and information from well profiles, allowing the acquisition of longitudinal conductance cartographies (S), data in mili-Siemens (mS), and the vulnerability of the aquifer. Such maps were prepared with emphasis to the unsaturated overlying zone, highlighting in particular its thickness and occurrence of clay lithologies. Thus, the longitudinal conductance and aquifer vulnerability reveal areas more susceptible to contamination in the northeast and east-central sections of the study area, with values equal to or less than 10mS and greater than or equal to 0,50, respectively. On the other hand, the southwestern section proved to be less susceptible to contamination, whose longitudinal conductance and vulnerability indices are greater than or equal to 30mS and less than or equal to 0,40, respectively.

Análise do grau de proteção do Aquífero Barreiras em perímetros sob fertilização química a partir de dados hidrogeofísicos - área da bacia do Rio Catu-RN

The use of chemical fertilization in arable perimeters provides increased productivity, though it can eventually lead to a qualitative depreciation of groundwater sources, especially if such sources are unconfined in nature. In this context, this thesis presents results from an analysis of the level of natural protection of the Barreiras Aquifer in an area located on the eastern coast of the Rio Grande do Norte State - Brazil. Such an aquifer is clastic in nature and has an unconfined hydraulic character, which clearly makes it susceptible to contamination from surface ground loads with contaminants associated with the leaching of excess fertilizers not absorbed by ground vegetation. The methodology used was based on the use of hydro-geophysical data, particularly inverse models of vertical electrical soundings (VES) and information from well profiles, allowing the acquisition of longitudinal conductance cartographies (S), data in mili-Siemens (mS), and the vulnerability of the aquifer. Such maps were prepared with emphasis to the unsaturated overlying zone, highlighting in particular its thickness and occurrence of clay lithologies. Thus, the longitudinal conductance and aquifer vulnerability reveal areas more susceptible to contamination in the northeast and east-central sections of the study area, with values equal to or less than 10mS and greater than or equal to 0,50, respectively. On the other hand, the southwestern section proved to be less susceptible to contamination, whose longitudinal conductance and vulnerability indices are greater than or equal to 30mS and less than or equal to 0,40, respectively.

Essays on Intergenerational and Regional Aspects of Water Management

This dissertation consists of three essays on different aspects of water management. The first essay focuses on the sustainability of freshwater use by introducing the notion that altruistic parents do bequeath economic assets for their offspring. Constructing a two-period, over-lapping generational model, an optimal ratio of consumption and pollution for old and young generations in each period is determined. Optimal levels of water consumption and pollution change according to different parameters, such as, altruistic degree, natural recharge rate, and population growth. The second essay concerns water sharing between countries in the case of trans-boundary river basins. The paper recognizes that side payments fail to forge water-sharing agreement among the international community and that downstream countries have weak bargaining power. An interconnected game approach is developed by linking the water allocation issue with other non-water issues such as trade or border security problems, creating symmetry between countries in bargaining power. An interconnected game forces two countries to at least partially cooperate under some circumstances. The third essay introduces the concept of virtual water (VW) into a traditional international trade model in order to estimate water savings for a water scarce country. A two country, two products and two factors trade model is developed, which includes not only consumers and producer’s surplus, but also environmental externality of water use. The model shows that VW trade saves water and increases global and local welfare. This study should help policy makers to design appropriate subsidy or tax policy to promote water savings especially in water scarce countries.

Table 1. pH, Temperature, conductivity and hydrochemistry of samples

Carbon dioxide deep geological storage, especially in deep saline aquifers, is one of the preferred technological options to mitigate the effects of greenhouse gases emissions. Thus, in the last decade, studies characterising the behaviour of potential CO2 deep geological storage sites along with thorough safety assessments have been considered essential in order to minimise the risks associated with these sites. The study of natural analogues represents the best source of reliable information about the expected hydrogeochemical processes involved in the CO2 storage in such deep saline aquifers. In this work, a comprehensive study of the hydrogeochemical features and processes taking place at the natural analogue of the Alicún de las Torres thermal system (Betic Cordillera) has been conducted. Thus, the main water/CO2/rock interaction processes occurring at the thermal system have been identified, quantified and modelled, and a principle conclusion is that the hydrogeochemical evolution of the thermal system is controlled by a global dedolomitization process triggered by gypsum dissolution. This geochemical process generates a different geochemical environment to that which would result from the exclusive dissolution of carbonates from the deep aquifer, which is generally considered as the direct result of CO2 injection in a deep carbonate aquifer. Therefore, discounting of the dedolomitization process in any CO2 deep geological storage may lead to erroneous conclusions. This process will also influence the porosity evolution of the CO2 storage formation, which is a very relevant parameter when evaluating a reservoir for CO2 storage. The geothermometric calculation performed in this work leads to estimate that the thermal water reservoir is located between 650 and 800 m depth, which is very close to the minimum required to inject CO2 in a deep geological storage. It is clear that the proper characterisation of the features and hydrogeochemical processes taking place at a natural system analogous to a man-made deep geological storage will provide useful conceptual, semi-quantitative and even quantitative information about the processes and consequences that may occur at the artificial storage system.

Table 2: Water stable isotope results and sampling site information of samples taken between 2003 and 2009 on the island of Corsica (western Mediterranean)

The Mediterranean is regarded as a region of intense climate change. To better understand future climate change, this area has been the target of several palaeoclimate studies which also studied stable isotope proxies that are directly linked to the stable isotope composition of water, such as tree rings, tooth enamel or speleothems. For such work, it is also essential to establish an isotope hydrology framework of the region of interest. Surface waters from streams and lakes as well as groundwater from springs on the island of Corsica were sampled between 2003 and 2009 for their oxygen and hydrogen isotope compositions. Isotope values from lake waters were enriched in heavier isotopes and define a local evaporation line (LEL). On the other hand, stream and spring waters reflect the isotope composition of local precipitation in the catchment. The intersection of the LEL and the linear fit of the spring and stream waters reflect the mean isotope composition of the annual precipitation (dP) with values of -8.6(±0.2) per mil for d18O and -58(±2) per mil for d2H. This value is also a good indicator of the average isotope composition of the local groundwater in the island. Surface water samples reflect the altitude isotope effect with a value of -0.17(±0.02) per mil per 100 m elevation for oxygen isotopes. At Vizzavona Pass in central Corsica, water samples from two catchments within a lateral distance of only a few hundred metres showed unexpected but systematic differences in their stable isotope composition. At this specific location, the direction of exposure seems to be an important factor. The differences were likely caused by isotopic enrichment during recharge in warm weather conditions in south-exposed valley flanks compared to the opposite, north-exposed valley flanks.

Hydrochemical properties of groundwater samples in Semarang area, Java Island, Indonesia (1992, 1993, 2003, 2006, and 2007)

In Semarang City, groundwater has been exploited as a natural resource since 1841. The groundwater exploited in deep wells is concentrated in confined aquifers. The previous hydrogeological model was developed in one unit of aquifer and refined then by using several hydrostratigraphical units following a regional hydrogeological map without any further analysis. At present, there is a lack of precise hydrogeological model which integrates geological and hydrogeological data, in particular for multiple aquifers in Semarang. Thus, the aim of this paper is to develop a hydrogeological model for the multiple aquifers in Semarang using an integrated data approach. Groundwater samples in the confined aquifers have been analyzed to define the water type and its lateral distribution. Two hydrogeological cross sections were then created based on several borelog data to define a hydrostratigraphical unit (HSU). The HSU result indicates the hydrogeological model of Semarang consists of two aquifers, three aquitards, and one aquiclude. Aquifer 1 is unconfined, while Aquifer 2 is confined. Aquifer 2 is classified into three groups (2a, 2b, and 2c) based on analyses of major ion content and hydrostratigraphical cross sections.

(Table 1) Sulfur contents and isotope composition in rocks and minerals from DSDP/ODP Hole 504B

DSDP Hole 504B is the only hole in oceanic crust to penetrate through the volcanic section and into hydrothermally altered sheeted dikes. We have carried out petrologic and sulfur isotopic analyses of sulfide and sulfate minerals and whole rocks from the core in order to place constraints on the geochemistry of sulfur during hydrothermal alteration of ocean crust. The nearly 600 m-thick pillow section has lost sulfur to seawater and has net d34S = -1.8 per mil due to degassing of SO2 during crystallization and subsequent low temperature interaction with seawater. Hydrothermally altered rocks in the 200 m-thick transition zone are enriched in S and 34S (4300 ppm and +3.0 +/-1.2 per mil, respectively), whereas the more than 500 m of sheeted dikes contain 720 ppm S with d34S = +0.6 +/-1.4 per mil. These data are consistent with the presence of predominantly basaltic sulfur in hydrothermal fluids deep in the crust: following precipitation of anhydrite during seawater recharge, small amounts of seawater sulfate were reduced at temperatures >250°C through conversion of igneous pyrrhotite to secondary pyrite and minor oxidation of ferrous iron in the crust. The S- and 34S-enrichments of the transition zone are the results of seawater sulfate reduction and sulfide deposition during subsurface mixing between upwelling hot (up to 350°C) hydrothermal fluids and seawater. Seawater sulfate was probably reduced through oxidation of ferrous iron in hydrothermal fluids and in the transition zone rocks. Alteration of the upper crust resulted in loss of basaltic sulfur to seawater, fixation of minor seawater sulfur in the crust and redistribution of magmatic sulfur within the crust. This caused net increases in sulfur content and d34S of the upper 1.8 km of the oceanic crust.

On the feasibility of time-lapse superconducting gravimetry for reservoir monitoring

The feasibility of monitoring fluid flow subsurface processes that result in density changes, using the iGrav superconducting gravimeter, is investigated. Practical targets include steam-assisted gravity drainage (SAGD) bitumen depletion and water pumping from aquifers, for which there is currently a void in low-impact, inexpensive monitoring techniques. This study demonstrates that the iGrav has the potential to be applied to multi-scale and diverse reservoirs. Gravity and gravity gradient signals are forward modeled for a real SAGD reservoir at two time steps, and for surface-fed and groundwater-fed aquifer pumping models, to estimate signal strength and directional dependency of water flow. Time-lapse gravimetry on small-scale reservoirs exhibits two obstacles, namely, a µgal sensitivity requirement and high noise levels in the vicinity of the reservoir. In this study, both limitations are overcome by proposing (i) a portable superconducting gravimeter, and (ii) a pair of instruments under various baseline geometries. This results in improved spatial resolution for locating depletion zones, as well as the cancellation of noise common in both instruments. Results indicate that a pair of iGrav superconducting gravimeters meet the sensitivity requirements and the spatial focusing desired to monitor SAGD bitumen migration at the reservoir scales. For SAGD reservoirs, the well pair separation, reservoir depth, and survey sampling determine the resolvability of individual well pair depletion patterns during the steam chamber rising phase, and general reservoir depletion patterns during the steam chamber spreading phase. Results show that monitoring water table elevation changes due to pumping and tracking whether groundwater or surface water is being extracted are feasible.