Application and evaluation of kriging and cokriging methods on groundwater depth mapping

Groundwater and water resources management play a key role in conserving the sustainable conditions in arid and semi-arid regions. Applying some techniques that can reveal the critical and hot conditions of water resources seem necessary. In this study, kriging and cokriging methods were evaluated for mapping the groundwater depth across a plain in which has experienced different climatic conditions (dry, wet, and normal) and consequently high variations in groundwater depth in a 12 year led in maximum, minimum, and mean depths. During this period groundwater depth has considerable fluctuations. Results obtained from geostatistical analysis showed that groundwater depth varies spatially in different climatic conditions. Furthermore, the calculated RMSE showed that cokriging approach was more accurate than kriging in mapping the groundwater depth though there was not a distinct difference. As a whole, kriging underestimated the real groundwater depth for dry, wet, and normal conditions by 5.5, 2.2, and 5.3%, while cokriging underestimations were 3.3, 2, and 2.2%, respectively; which showed the unbiasedness in estimations. Results implied that in the study area farming and cultivation in dry conditions needs more attention due to higher variability in groundwater depth in short distances compared to the other climate conditions. It is believed that geostatistical approaches are reliable tools for water resources managers and water authorities to allocate groundwater resources in different environmental conditions.     

SPATIAL ANALYSIS AND EUTROPHICATION ASSESSMENT FOR CHLOROPHYLL a IN TAIHU LAKE

Spatial structure analysis and kriging analysis have been identified to be useful tools in illustrating the spatial patterns of variables. Taihu Lake is one of the largest fresh water lakes in China, and has suffered serious eutrophication in recent years due to the rapid economic development and growing environmental pollution in the Taihu Catchment. In this paper, spatial structural analysis, kriging interpolation and eutrophication assessment were carried out for chlorophyll a in the lake. Studies show that spherical model could be applied to fit all experimental variograms. Positive nuggets were observed for three directions except NE–SW direction. The variograms show some anisotropy with anisotropic ratio falling within 1.76. The spatial structural patterns of chlorophyll a in the lake were affected by factors such as distribution of pollution sources, water flow and wind. Two-dimensional ordinary block kriging was applied for interpolation process. An eutrophication assessment map was also made based on a water-quality evaluation standard. Results show that the content of chlorophyll a in Taihu Lake was quite high. The whole lake has suffered serious eutrophication. However, the eutrophic situation varied in space. Higher contents of chlorophyll a appeared mainly in the northern part of the lake.     

Groundwater levels time series sensitivity to pluviometry and air temperature: a geostatistical approach to Sfax region,Tunisia

In this paper, the pattern of groundwater level fluctuations is investigated by statistical techniques for 24 monitoring wells located in an unconfined coastal aquifer in Sfax (Tunisia) for a time period from 1997 to 2006. Firstly, a geostatistical study is performed to characterize the temporal behaviors of data sets in terms of variograms and to make predictions about the value of the groundwater level at unsampled times. Secondly, multivariate statistical methods, i.e., principal component analysis (PCA) and cluster analysis (CA) of time series of groundwater levels are used to classify groundwater hydrographs regard to identical fluctuation pattern. Three groundwater groups (A, B, and C) were identified. In group “A,” water level decreases continuously throughout the study periods with rapid annual cyclic variation, whereas in group “B,” the water level contains much less high-frequency variation. The wells of group “C” represents a steady and gradual increase of groundwater levels caused by the aquifer artificial recharge. Furthermore, a cross-correlation analysis is used to investigate the aquifer response to local rainfall and temperature records. The result revealed that the temperature is more affecting the variation of the groundwater level of group A wells than the rainfall. However, the second and the third groups are less affected by rainfall or temperature.     

Indicator and probability kriging methods for delineating Cu, Fe, and Mn contamination in groundwater of Najafgarh Block, Delhi, India

Two non-parametric kriging methods such as indicator kriging and probability kriging were compared and used to estimate the probability of concentrations of Cu, Fe, and Mn higher than a threshold value in groundwater. In indicator kriging, experimental semivariogram values were fitted well in spherical model for Fe and Mn. Exponential model was found to be best for all the metals in probability kriging and for Cu in indicator kriging. The probability maps of all the metals exhibited an increasing risk of pollution over the entire study area. Probability kriging estimator incorporates the information about order relations which the indicator kriging does not, has improved the accuracy of estimating the probability of metal concentrations in groundwater being higher than a threshold value. Evaluation of these two spatial interpolation methods through mean error (ME), mean square error (MSE), kriged reduced mean error (KRME), and kriged reduced mean square error (KRMSE) showed 3.52% better performance of probability kriging over indicator kriging. The combined result of these two kriging method indicated that on an average 26.34%, 65.36%, and 99.55% area for Cu, Fe, and Mn, respectively, are coming under the risk zone with probability of exceedance from a cutoff value is 0.6 or more. The groundwater quality map pictorially represents groundwater zones as ??desirable?? or ??undesirable?? for drinking. Thus the geostatistical approach is very much helpful for the planners and decision makers to devise policy guidelines for efficient management of the groundwater resources so as to enhance groundwater recharge and minimize the pollution level.     

Spatial assessment of water quality parameters in Jhelum city (Pakistan)

In this study, we assess the drinking water quality of Jhelum city. Two hundred and ninety-two drinking water samples were randomly collected in the study area. These samples were chemically analyzed for three key toxic (in excess) elements such as pH, total dissolved solids (TDS), and calcium. Geostatistical techniques such as variogram and kriging were used to investigate the spatial variations of these minerals across the city. The spatial structure for each element was found to be anisotropic, and thus, anisotropic variograms were used. The kriging predictions revealed significant concentrations of the above-stated elements at some locations in the study area. While comparing with the World Health Organization, United States Environmental Protection Agency, and Pakistan Environmental Protection Agency standards, the water samples were found to be unsatisfactory for drinking. We conclude that the drinking water in this region is of poor quality and needs proper treatment to make it palatable.     

Spatial evaluation of the risk of groundwater quality degradation. A comparison between disjunctive kriging and geostatistical simulation

In some previous papers a probabilistic methodology was introduced to estimate a spatial index of risk of groundwater quality degradation, defined as the conditional probability of exceeding assigned thresholds of concentration of a generic chemical sampled in the studied water system. A crucial stage of this methodology was the use of geostatistical techniques to provide an estimation of the above-mentioned probability in a number of selected points by crossing spatial and temporal information. In this work, spatial risk values were obtained using alternatively stochastic conditional simulation and disjunctive kriging. A comparison between the resulting two sets of spatial risks, based on global and local statistical tests, showed that they do not come from the same statistical population and, consequently, they cannot be viewed as equivalent in a statistical sense. At a first glance, geostatistical conditional simulation may appear to represent the spatial variability of the phenomenon more effectively, as the latter tends to be smoothed by DK. However, a close examination of real case study results suggests that disjunctive kriging is more effective than simulation in estimating the spatial risk of groundwater quality degradation. In the study case, the potentially ‘harmful event’ considered, threatening a natural ‘vulnerable groundwater system,’ is fertilizer and manure application.     

Comparison of Groundwater Level Estimation Using Neuro-fuzzy and Ordinary Kriging

Water level in aquifer plays the main role in groundwater modeling as one of the input data. In practice, due to aspects of time and cost, data monitoring of water levels is conducted at a limited number of sites, and interpolation technique such as kriging is widely used for estimation of this variable in unsampled sites. In this study, the efficiency of the ordinary kriging (OK) and adaptive network-based fuzzy inference system (ANFIS) was investigated in interpolation of groundwater level in an unconfined aquifer in the north of Iran. The results showed that ANFIS model is more efficient in estimating the groundwater level than OK.     

Estimation of spatial distrubition of groundwater level and risky areas of seawater intrusion on the coastal region in Çarşamba Plain,Turkey, using different interpolation methods

Groundwater level plays a significant role in coastal plains. Heavy pumping and excessive use of near-coast groundwater can increase the intrusion of seawater into the aquifers. In the present study, groundwater levels were measured at 59 groundwater wells at different times during pre- and post-irrigation seasons (April and September of the year 2012) in Çar?amba Plain, Turkey. To select the best method, two deterministic interpolation methods (inverse distance weighing (IDW) with the weights of 1, 2, and 3 and radial basis function (RBF) with spline with tension (SPT) and completely regularized spline (CRS)) and two stochastic methods (ordinary kriging (OK) with spherical, exponential, and Gaussian variograms) and cokriging (COK)) were compared and then the best interpolation method was used to evaluate the spatial distribution of groundwater levels in different seasons and seasonal changes. A total of nine different techniques were tested. Also, risky areas of seawater intrusion in coastal area were determined using the best methods for two periods. The performance of these interpolation methods is evaluated by using a validation test method. Statistical indices of correlation (R ²), mean absolute error (MAE), and root-mean-square error (RMSE) were used to select and validate the best methods. Comparisons between predicted and observed values indicated RBF as the optimal method for groundwater level estimation in April and September. When the best method RBF and the worst method IDW were compared, significant differences were observed in the spatial distribution of groundwater. Results of the study also revealed that excessive groundwater withdrawals during the post-irrigation season dropped the groundwater levels up to 2.0 m in some sections. With regard to seawater intrusion, 9,103 ha of land area was determined to be highly risky and risky.     

GIS-based assessment and characterization of groundwater quality in a hard-rock hilly terrain of Western India

The growing population, pollution, and misuse of freshwater worldwide necessitate developing innovative methods and efficient strategies to protect vital groundwater resources. This need becomes more critical for arid/semi-arid regions of the world. The present study focuses on a GIS-based assessment and characterization of groundwater quality in a semi-arid hard-rock terrain of Rajasthan, western India using long-term and multi-site post-monsoon groundwater quality data. Spatio-temporal variations of water quality parameters in the study area were analyzed by GIS techniques. Groundwater quality was evaluated based on a GIS-based Groundwater Quality Index (GWQI). A Potential GWQI map was also generated for the study area following the Optimum Index Factor concept. The most-influential water quality parameters were identified by performing a map removal sensitivity analysis among the groundwater quality parameters. Mean annual concentration maps revealed that hardness is the only parameter that exceeds its maximum permissible limit for drinking water. GIS analysis revealed that sulfate and nitrate ions exhibit the highest (CV?>?30%) temporal variation, but groundwater pH is stable. Hardness, EC, TDS, and magnesium govern the spatial pattern of the GWQI map. The groundwater quality of the study area is generally suitable for drinking and irrigation (median GWQI?>?74). The GWQI map indicated that relatively high-quality groundwater exists in northwest and southeast portions of the study area. The groundwater quality parameter group of Ca, Cl, and pH were found to have the maximum value (6.44) of Optimum Index factor. It is concluded that Ca, Cl, and pH are three prominent parameters for cost-effective and long-term water quality monitoring in the study area. Hardness, Na, and SO₄, being the most-sensitive water quality parameters, need to be monitored regularly and more precisely.     

Spatiotemporal variation of groundwater quality using integrated multivariate statistical and geostatistical approaches in Amol–Babol Plain,Iran

In recent years, groundwater quality has become a global concern due to its effect on human life and natural ecosystems. To assess the groundwater quality in the Amol–Babol Plain, a total of 308 water samples were collected during wet and dry seasons in 2009. The samples were analysed for their physico-chemical and biological constituents. Multivariate statistical analysis and geostatistical techniques were applied to assess the spatial and temporal variabilities of groundwater quality and to identify the main factors and sources of contamination. Principal component analysis (PCA) revealed that seven factors explained around 75 % of the total variance, which highlighted salinity, hardness and biological pollution as the dominant factors affecting the groundwater quality in the Plain. Two-way analysis of variance (ANOVA) was conducted on the dataset to evaluate the spatio-temporal variation. The results showed that there were no significant temporal variations between the two seasons, which explained the similarity between six component factors in dry and wet seasons based on the PCA results. There are also significant spatial differences (p?>?0.05) of the parameters under study, including salinity, potassium, sulphate and dissolved oxygen in the plain. The least significant difference (LSD) test revealed that groundwater salinity in the eastern region is significantly different to the central and western side of the study area. Finally, multivariate analysis and geostatistical techniques were combined as an effective method for demonstrating the spatial structure of multivariate spatial data. It was concluded that multiple natural processes and anthropogenic activities were the main sources of groundwater salinization, hardness and microbiological contamination of the study area.     

Evaluating Kriging as a Tool to Improve Moderate Resolution Maps of Forest Biomass

The USDA Forest Service, Forest Inventory and Analysis program (FIA) recently produced a nationwide map of forest biomass by modeling biomass collected on forest inventory plots as nonparametric functions of moderate resolution satellite data and other environmental variables using Cubist software. Efforts are underway to develop methods to enhance this initial map. We explored the possibility of modeling spatial structure to make such improvements. Spatial structure in the field biomass data as well as in residuals from the map was investigated across 18 ecological zones in the Interior Western U.S. Exploratory tools included directional graphs of summary statistics, three dimensional maps, Moran’s I correlograms, and variograms. Where spatial pattern was present, field and residual biomass were kriged, and predictions made for an independent test set were evaluated for improvement over predictions in the initial biomass map. While kriging has some potential benefit when analyzing the field data and exploring spatial structure, kriging residuals resulted in little or no improvement in the initial biomass map developed using Cubist software. Stationarity assumptions, variogram behavior, and appropriate model fitting strategies are discussed.     

Spatio-temporal variations in water quality of Nullah Aik-tributary of the river Chenab, Pakistan

This study reports the spatio-temporal changes in water quality of Nullah Aik, tributary of the Chenab River, Pakistan. Stream water samples were collected at seven sampling sites on seasonal basis from September 2004 to April 2006 and were analyzed for 24 water quality parameters. Most significant parameters which contributed in spatio-temporal variations were assessed by statistical techniques such as Hierarchical Agglomerative Cluster Analysis (HACA), Factor Analysis/Principal Components Analysis (FA/PCA), and Discriminant Function Analysis (DFA). HACA identified three different classes of sites: Relatively Unimpaired, Impaired and Less Impaired Regions on the basis of similarity among different physicochemical characteristics and pollutant level between the sampling sites. DFA produced the best results for identification of main variables for temporal and spatial analysis and separated eight parameters (DO, hardness, sulphides, K, Fe, Pb, Cr and Zn) that accounted 89.7% of total variations of spatial analysis. Temporal analysis using DFA separated six parameters (E.C., TDS, salinity, hardness, chlorides and Pb) that showed more than 84.6% of total temporal variation. FA/PCA identified six significant factors (sources) which were responsible for major variations in water quality dataset of Nullah Aik. The results signify that parameters identified by statistical analyses were responsible for water quality change and suggest the possibility of industrial, municipal and agricultural runoff, parent rock material contamination. The results suggest dire need for proper management measures to restore the water quality of this tributary for a healthy and promising aquatic ecosystem and also highlights its importance for objective ecological policy and decision making process.     

Groundwater quality assessment and its correlation with gastroenteritis using GIS: a case study of Rawal Town,Rawalpindi, Pakistan

Majority of the people of Pakistan get drinking water from groundwater source. Nearly 40 % of the total ailments reported in Pakistan are the result of dirty drinking water. Every summer, thousands of patients suffer from acute gastroenteritis in the Rawal Town. Therefore, a study was designed to generate a water quality index map of the Rawal Town and identify the relationship between bacteriological water quality and socio-economic indicators with gastroenteritis in the study area. Water quality and gastroenteritis patient data were collected by surveying the 262 tubewells and the major hospitals in the Rawal Town. The collected spatial data was analyzed by using ArcGIS spatial analyst (Moran’s I spatial autocorrelation) and geostatistical analysis tools (inverse distance weighted, radial basis function, kriging, and cokriging). The water quality index (WQI) for the study area was computed using pH, turbidity, total dissolved solids, calcium, hardness, alkalinity, and chloride values of the 262 tubewells. The results of Moran’s I spatial autocorrelation showed that the groundwater physicochemical parameters were clustered. Among IDW, radial basis function, and kriging and cokriging interpolation techniques, cokriging showed the lowest root mean square error. Cokriging was used to make the spatial distribution maps of water quality parameters. The WQI results showed that more than half of the tubewells in the Rawal Town were providing “poor” to “unfit” drinking water. The Pearson’s coefficient of correlation for gastroenteritis with fecal coliform was found significant (P?P?P?相似文献   

Comparison of stochastic and deterministic methods for mapping groundwater level spatial variability in sparsely monitored basins

In sparsely monitored basins, accurate mapping of the spatial variability of groundwater level requires the interpolation of scattered data. This paper presents a comparison of deterministic interpolation methods, i.e. inverse distance weight (IDW) and minimum curvature (MC), with stochastic methods, i.e. ordinary kriging (OK), universal kriging (UK) and kriging with Delaunay triangulation (DK). The study area is the Mires Basin of Mesara Valley in Crete (Greece). This sparsely sampled basin has limited groundwater resources which are vital for the island’s economy; spatial variations of the groundwater level are important for developing management and monitoring strategies. We evaluate the performance of the interpolation methods with respect to different statistical measures. The Spartan variogram family is applied for the first time to hydrological data and is shown to be optimal with respect to stochastic interpolation of this dataset. The three stochastic methods (OK, DK and UK) perform overall better than the deterministic counterparts (IDW and MC). DK, which is herein for the first time applied to hydrological data, yields the most accurate cross-validation estimate for the lowest value in the dataset. OK and UK lead to smooth isolevel contours, whilst DK and IDW generate more edges. The stochastic methods deliver estimates of prediction uncertainty which becomes highest near the southeastern border of the basin.     

面向水源保护与污染监控的广东地下水环境监测网研究

通过系统整理分析自然资源、水利、生态环境等部门已有地下水监测网,针对广东地下水监测存在的问题和未来地下水环境管理需求,提出较大空间尺度和水质监控的地下水环境监测网建设框架及广东省"三层七级"共237个网点的层级化地下水环境监测网建设方案,并与水文地质单元、浅层地下水功能区划进行了空间叠加分析。结果表明,广东省三大部门监测网在空间上分布总体较为分散,相互补充性好,提出的监测网方案总体上与广东省地下水本底条件、地下水资源开发利用现状、地下水环境功能分区、重点污染源空间分布特征等相匹配,契合地下水饮用水源保护和重点污染源监控两大管理需求,构建的地下水环境监测网框架可为其他地区提供参考,提出的地下水环境监测网建设方案可为广东省地下水环境管理和污染防控提供决策支持。     

Assessment of groundwater pollution in West Delhi, India using geostatistical approach

The exploration, exploitation, and unscientific management of groundwater resources in the National Capital Territory (NCT) of Delhi, India have posed a serious threat of reduction in quantity and deterioration of quality. The objective of the study is to determine the groundwater quality and to assess the risk of groundwater pollution at Najafgarh, NCT of Delhi. The groundwater quality parameters were analyzed from the existing wells of the Najafgarh and the thematic maps were generated using geostatistical concepts. Ordinary kriging and indicator kriging methods were used as geostatistical approach for preparation of thematic maps of the groundwater quality parameters such as bicarbonate, calcium, chloride, electrical conductivity (EC), magnesium, nitrate, sodium, and sulphate with concentrations equal or greater than their respective groundwater pollution cutoff value. Experimental semivariogram values were fitted well in spherical model for the water quality parameters, such as bicarbonate, chloride, EC, magnesium, sodium, and sulphate and in exponential model for calcium and nitrate. The thematic maps of all the groundwater quality parameters exhibited an increasing trend of pollution from the northern and western part of the study area towards the southern and eastern part. The concentration was highest at the southernmost part of the study area but it could not reflect correctly the groundwater pollution status. The indicator kriging method is useful to assess the risk of groundwater pollution by giving the conditional probability of concentrations of different chemical parameters exceeding their cutoff values. Thus, risk assessment of groundwater pollution is useful for proper management of groundwater resources and minimizing the pollution threat.     

Earth-tide-induced fluctuations in the salinity of an inland river,New South Wales,Australia: a short-term study

Wybong Creek discharges salts into the agriculturally and industrially important Hunter River in New South Wales, Australia. Abrupt increases in salinity occur periodically in the mid-Wybong Creek catchment. In order to understand the processes which cause these abrupt increases, changes in surface and groundwater were investigated. It is shown that salinity increases can be attributed to highly discrete groundwater discharge directly into the river from below. Hourly electrical conductivity data measured in the river showed regular, diurnal electrical conductivity fluctuations of up to 350 μS cm^?1. These fluctuations could not be attributed to barometric pressure, temperature, or evapotranspiration. Instead, a similar periodicity in surface water electrical conductivity and groundwater height in nearby groundwater wells was found. Fluctuations were of similar periodicity to the orthotides which cause fluctuations in surface water height and are induced by Earth tides. The geology in the mid-catchment area indicates conditions are optimal for Earth tides to impact groundwater. The reporting of orthotidal changes in water chemistry in this article is believed to be the first of its kind in the scientific literature, with the large fluctuations noted having important implications for water monitoring and management in the catchment. Further research investigating Earth-tide-induced phases of groundwater heights will better constrain the relationships between surface water chemistry and groundwater height.     

Geostatistical Evaluation of Spatial Variation Related to Groundwater Quality Database: Case Study for Arak Plain Aquifer,Iran

Geostatistical methods are one of the advanced techniques to interpolate groundwater quality data. Geostatistical interpolation techniques employ both the mathematical and the statistical properties of the measured points. Compiling the data distribution on spatial and temporal domain is of crucial importance in order to evaluate its quality and safety. The main purpose of this paper is to assess groundwater quality of Arak plain, Iran, by an unbiased interpolated method so called Kriging. Therefore, seven quality variables of Arak plain aquifer including TDS, SAR, EC, Na⁺, TH, Cl^?, and SO₄ ^2? have been analyzed, studied, and interpreted statistically and geostatistically. Utilized data in this study were collected from 97 water well samples in Arak plain, in 2012. After normalizing data, variogram as a geostatistical tool for defining spatial regression was calculated and experimental variograms have been plotted by GS⁺ software, then the best theoretical model was fitted to each variogram based on minimum RSS error. Cross validation was used to determine the accuracy of the estimated data. The uncertainty of the method could be well assessed via this method since the method not only gave the average error (around 0 in this study) but also gave the standard deviation of the estimations. Therefore, more than 3800 points were estimated by ordinary Kriging algorithm in places which have not been sampled. Finally, estimation maps of groundwater quality were prepared and map of estimation variance, EV, has been presented to assess the quality of estimation in each estimated point. Results showed that the Kriging method is more accurate than the traditional interpolation algorithms not honoring the spatial properties of the database.     

Estimation of nested spatial patterns and seasonal variation in the longitudinal distribution of Sicyopterus japonicus in the Datuan Stream, Taiwan by using geostatistical methods

This study attempts to determine the scale-dependent hierarchical spatial variation and longitudinal distributions of Sicyopterus japonicus year round. The distribution of S. japonicus in the Datuan Stream in northern Taiwan was surveyed during the fall and winter 2007, as well as the spring and summer of 2008. The spatial structure of S. japonicus density was modeled using geostatistics. The longitudinal distributions of S. japonicus density were then estimated using kriging and hydrology distance with nested variogram models. Variography results indicate that nested variogram models could reflect the hierarchical structure in the spatial variation of seasonal S. japonicus density, with the small, median, and large ranges representing three nested scales. Models for the four seasons were consistent in that they shared the same shape of variogram models with various ranges and sill values. This model shape consistency implies stationary spatial correlations in the longitudinal fish distribution across the four seasons. The Kriging geostatistical method based on the multiple scales nested variogram models also provided robust estimates of S. japonicus densities at unsampled sections. We conclude that S. japonicus densities exhibit hierarchical patterns and variation in the four seasons along the study stream. Geostatistical methods with a nested variograms and hydrological distance are a highly effective means of delineating the hierarchical structure in longitudinal patterns of S. japonicus density in each season, providing estimates of the S. japonicus density for hierarchically structured spatial distributions and expanding knowledge of S. japonicus beyond the limits imposed by spatial and temporal scales.     

Tidal along-shore groundwater flow in a coastal aquifer

Cross-shore interactions between the ocean and a coastal aquifer have been studied extensively, whereas the corresponding along-shore case has seldom been examined. This paper presents a numerical model that simulates two-dimensional groundwater flow averaged over the thickness of a coastal aquifer. The model is used to examine the essential features of tide-induced, along-shore effects on an aquifer adjacent to a cross-shore river. The results show that the tide, which fluctuates the water level in the river, induces groundwater table fluctuations and oscillating flows in the along-shore direction. This occurs even at locations much further inland than tidal cross-shore fluctuations can propagate. However, the magnitude of along-shore water table fluctuations and flow velocity at a given cross-shore distance decreases with the distance from the river in the same manner as cross-shore tidal fluctuations. The along-shore groundwater flow, together with the cross-shore flow, forms horizontal circulation and increases mixing of solute in the aquifer. Over a tidal period, a large amount of water exchange occurs at the river-aquifer and ocean-aquifer interfaces, leading to increased transfer of chemicals between the three water bodies. These results have implications for the management of waste discharge in estuaries and coastal aquifers.