Next Article in Journal
Hypochlorite Generation from a Water Softener Spent Brine
Next Article in Special Issue
Halloysite Nanotubes as Adsorptive Material for Phosphate Removal from Aqueous Solution
Previous Article in Journal
Evaluating Glacier Volume Changes since the Little Ice Age Maximum and Consequences for Stream Flow by Integrating Models of Glacier Flow and Hydrology in the Cordillera Blanca, Peruvian Andes
Previous Article in Special Issue
Review of Heavy Metal Adsorption Processes by Several Organic Matters from Wastewaters
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Indication of Groundwater Contamination Using Acesulfame and Other Pollutants in a Rural Area of Korea

1
Department of Geological Sciences, Pusan National University, Busan 46241, Korea
2
Department of Civil and Environmental Engineering, Pusan National University, Busan 46241, Korea
3
Soil and Groundwater Division, National Institute of Environmental Research, Incheon 22689, Korea
*
Author to whom correspondence should be addressed.
Water 2018, 10(12), 1731; https://doi.org/10.3390/w10121731
Submission received: 1 October 2018 / Revised: 3 November 2018 / Accepted: 21 November 2018 / Published: 26 November 2018
(This article belongs to the Special Issue Water Pollution and Treatment: Challenges and Opportunities)

Abstract

:
Approximately 40,000 chemical products are currently used in Korea; these products can contaminate the groundwater/soil, the surrounding environment, and organisms for extended periods of time. In this study, a hydrological field survey, a water quality analysis, and groundwater modeling were performed to identify the source and transport path of pollution that was caused by inorganic matter and artificial sweeteners, especially acesulfame, in the groundwater of an agricultural area in Chungnam Province, Korea. In the study area, a higher concentration of acesulfame displayed a spatial distribution similar to nitrate-nitrogen concentration. The characteristics of the groundwater flow and the distribution of the acesulfame were simulated using the Visual MODFLOW Classic Interface ver. 2014.1 and the MT3DMS module, respectively. The modeled area was divided into hilly (southern), residential (northwest), and agricultural (northeast) zones. The stream’s boundary was set to be the drainage channel in the southern hilly zone. From the simulation, we found that acesulfame spread actively from the source for 1–3 years before it reached equilibrium in the northern part of the model domain (the area downstream of the stream’s boundary). The concentration of acesulfame in the agricultural zone of the model domain decreased after five years, and it reached the steam boundary and residential zone within 10 years. After 10 years, most of the acesulfame was discharged from the agricultural zone and the hilly zone, while the concentration in the residential zone was approximately the same. Acesulfame is considered to be a potential indicator of man-made contamination for use in the management of groundwater quality.

1. Introduction

With the development of industries, the use of a diverse range of chemicals has increased and presents an environmental contamination risk. Since 2009, the residues of artificial sweeteners in the environment have been actively investigated using various water samples, including sewage water, wastewater, river water, and groundwater [1,2]. Research on artificial sweeteners as indicators of wastewater has continuously been carried out [3]. Table 1 shows the concentrations of artificial sweeteners, such as acesulfame, cyclamate, saccharin, and sucralose, in groundwater in some countries [4,5,6,7,8]. Acesulfame, which has low adsorption in soil due to its high solubility and low adsorption [2] and is highly persistent due to the low biodegradation during sewage treatment [5,9,10,11], is known to be a pollutant in groundwater and surface water due to sewage inflow [12,13]. The Swiss Federal Research Center for Agroecology and Agriculture in Zurich has reported that acesulfame is a highly persistent material that is detected in high concentrations in the environment and can be used to track wastewater. In Switzerland, the acesulfame that had not been eliminated from the wastewater treatment plants, and was quite persistent in surface waters, showed an increased concentration with a high density of population in the catchment area and a decreased concentration in the water throughflow area, and it was consistently detected in untreated and treated wastewater samples (12–46 μg/L), in most surface waters samples, in 65% of the investigated groundwater samples, and even in several tap water samples (up to 2.6 μg/L).
Medicinal substances in groundwater have been studied in European countries, the United States, China, etc. [14,15,16,17,18]. Some medicinal materials have also been studied in Korea [19,20]. The use of veterinary medicines, which leave residues in the environment, is increasing [21]. Numerous studies [22,23,24] have focused on antibiotics, since antibiotic-resistant bacteria have continuously been detected in the environment. In 2011, the World Health Organization (WHO) [25] highlighted the need for studies monitoring the spread of antibiotics and their propagation pathways. Antibiotics in animals account for approximately 50% of the total antibiotic use worldwide [26]. In the U.S., between 6350 and 13,380 tons of antibiotics are used on farm animals each year [27]. Bartelt-Hunt et al. [28], Batt et al. [29], Hu et al. [22], and Watanabe et al. [30] studied the animal medicinal substances found in groundwater that was collected from animal husbandry farms and detected antibiotics, such as sulfonamides, microlides, and tetracyclines. In particular, sulfamethazine, a sulfonamide antibiotic, was detected at a relatively high concentration along with monensin, a feed additive, and hormonal drugs. Barnes et al. [31] found various medicinal substances, including antibiotics and hormones, in 80% of river water samples and 95% of groundwater samples. In particular, in 2002, a high concentration of sulfamethazine was found in the groundwater in agricultural and livestock areas.
Nitrate, which is one of the most common inorganic pollutants in groundwater, creates a significant water-quality problem in diverse land-use areas, especially in agricultural regions [32]. Nitrate in groundwater originates from non-point sources like inorganic fertilizers and point sources like sewage from septic tanks and animal/human excrement [33,34,35]. Kim et al. [36] reported the significant potential of nitrate pollution in groundwater by the application of nitrogen fertilizer (>250 N kg/ha) in agricultural areas of Korea. Cheong et al. [37] characterized nitrate contamination in groundwater in the Gimpo agricultural area in Korea, by using a risk assessment method.
In Korea, the residual status, the behavioral characteristics, and the risk assessment of medicinal substances have been studied by investigating livestock and marine products, as well as the environment [38,39]. In particular, the National Institute of Environmental Research (NIER) has studied the residual state and behavior of medicinal substances in the environment several times. For instance, between 2008 and 2011, they monitored medicinal substances in water systems and surveyed the production/use of medicines in Korea through an ‘Investigation on the residual status and the chemical analysis methodology of residual medicinal substances’. In the context of this investigation, 40 medicinal chemicals, including analgesics, antipyretics, anti-inflammatories, disinfectants, anti-depressants, and beta-blockers for human use, antibiotics, antimicrobials, analgesics, antipyretics, and anti-inflammatories for human/animal use, and antibiotics and antimicrobials for animal use, were surveyed in groundwater [40,41].
Korea’s Groundwater Act prescribes 19 water quality standards for 19 chemicals including nitrate nitrogen (nitrate-N) [42]; however, it needs to include more water quality standards for medicinal substances, artificial sweeteners, pesticides, etc. The NIER [43] found two representative indicators among the drinking water quality standards for agricultural areas in Korea: nitrate-N and the total coliform group. Korea’s Animal and Plant Quarantine Agency conducted a survey on the behavior of drugs for animals in the soil and groundwater around livestock farms for two years (2009–2010) and, in the groundwater around pig and poultry farms, found a high concentration of medicinal substances for livestock, such as sulfamethazine and sulfonamide-based antibiotics [20].
This study examined the characteristics and distribution of acesulfame and other pollutants in the groundwater in a rural area in Chungnam Province, Korea. For this study, groundwater samples from eight wells were collected three times in 2017. Water quality analyses (cations, ions, and artificial sweeteners), a hydrogeological survey (a slug test and a tracer test), a grain size analysis, and groundwater modeling were carried out. Temperature, pH, and electrical conductivity (EC) were measured in the field. K+, Na+, Ca2+, Mg2+, HCO3, CO32−, SO42−, Cl, acesulfame, cyclamate, and saccharin were analyzed in the laboratory. For this study, nitrate-N concentrations by the previous study [43] were also compared with the acesulfame concentrations.

2. Materials and Methods

The study area, S village in O Myeon, Buyeo-gun, Chungcheongnam-do province, is composed of a hilly zone (southern part of study area) at an altitude of ~230 m above the mean seal level (a m.s.l.), a residential zone that is 20–40 m a m.s.l., and an agricultural zone that is 0–20 m a m.s.l. Streams flow from the south to the northeast. The main pollution sources are the farmland, surrounding sloped forest orchards that mostly produce chestnuts and one livestock farm that uses composted manure. The total watershed has an area of 458,567 m2. Farmland covers 226,466 m2, which is 49.6% of the total area; this farmland is 70% paddies, 18% dry agricultural land, and 12% orchard [43]. Rice, garlic, sesame, onions, potatoes, sweet potatoes, etc., are cultivated in the agricultural area. Eight wells (BS-3, BS-5, BS-6, BS-8, BS-9, BS-M1, BS-M2, and BS-M3) were sampled three times in May, August, and October in 2017 (Figure 1).

2.1. Water Quality Analysis

Groundwater samples from eight wells were collected for a laboratory analysis with the physicochemical components (pH, EC, and water temperature) in a stabilized state and at the same locations for an in situ water quality analysis (Figure 1). The samples were filtered using 45-μm filter paper. In addition, the samples that were used for a cation analysis were treated with a 0.05 N nitric acid solution in order to prevent the adsorption of the cations on the sample bottle. The filtered samples were stored in 200-mL sterilized bottles and refrigerated at or below 4 °C as they were transported to the laboratory for analysis. Four cations (K+, Na+, Ca2+, and Mg2+) and four anions (HCO3, CO32−, SO42−, and Cl) in the groundwater samples were analyzed by using an inductively-coupled plasma atomic emission spectrometer (Model ICP-IRIS, Thermo Jarrell Ash Co., Cambridge, MA, USA) and ion chromatography (Model DX-500, Dionex Co., Sunnyvale, CA, USA), respectively, in the Core Research Facilities at Pusan National University (PNU). We also analyzed the artificial sweeteners acesulfame, cyclamate, and saccharin. The artificial sweeteners and the triclosan were analyzed by using an Agilent 1200 HPLC (Agilent Technologues, Santa Clara, CA, USA) for liquid chromatography, and an Agilent 6460 triple-quadrupole mass spectrometer (Agilent Technologues, Santa Clara, CA, USA) for detection, in the laboratory of the Department of Environmental Engineering at PNU (Table 2).

2.2. Hydrogeological Survey

Hydraulic conductivity can be determined by using hydraulic field tests for bedrock aquifers, such as a pumping test, a slug test, and a water pressure test, and hydraulic field tests for unconsolidated formations, such as a con penetration test [44,45,46], a pumping test [47,48], a slug test, and a grain size analysis [49]. Specifically, the hydraulic conductivity of contaminated soils can be evaluated by [50,51]. In this study, the hydraulic conductivity values were determined by a using grain size analysis and a slug test. Besides this, a tracer test was carried out to estimate the longitudinal dispersivity and the velocity of solutes.

2.2.1. The Grain Size Analysis and the Slug Test

A sieve analysis of dry grains was performed using sieve numbers of up to 200 (0.074 mm), and the size of grains that were smaller than 0.074 mm was determined using a laser grain size analysis. First, the samples were sifted in order of grain size, and then, the material that passed through the sieves was weighed. Table 3 shows the sieve numbers and the mesh sizes that were used in this study.
The hydraulic conductivity, K (LT−1), can be estimated using four empirical formulas (Table 4) that follow from the grain size analysis based on Equation (1) [52]:
  K = ( g μ ) × β   × v ( n ) × d e 2  
where g is the gravitational acceleration (LT−2); μ is the kinematic viscosity (L2T−1); β is a dimensionless coefficient that depends on the texture, particle shape, rock composition, anisotropy, etc.; v(n) is the porosity function; n is the porosity (dimensionless); and de is the effective diameter of the porous medium (L). In this study, the μ value of groundwater was applied based on the average temperature of groundwater (15 °C).
A slug test is a cost-effective hydraulic field test that is used to obtain the hydraulic parameters. During the test, a dummy or water is instantaneously inserted into the water of the borehole, and the change in the water level is measured. From the moment when the dummy rises above the water, the water level recovery is measured. Hydraulic parameters are commonly estimated using the Hvorslev [57], Bower and Rice [58], or Cooper–Bredehoeft–Papadopulos [59] method.

2.2.2. Tracer Test

Common tracers in groundwater surveys include chloride (Cl), bromide (Br), lithium (Li+), ammonium (NH4+), magnesium (Mg2+), potassium (K+), iodide (I), and sulfate (SO42−). The dye tracers fluorescein, pyranine, lissamine FF, rhodamine B, rhodamine WT, and sulfo-rhodamine B have high detectability, enable rapid field analysis, are low cost, and have low toxicity [60].
In this study, a natural gradient tracer test was used, which is one of the multi-well tracer tests. The test needs one injection well and at least one observation well. For the test, a small amount of tracer is inserted into the aquifer in order to monitor its transport through the natural groundwater flow. In principle, the observation and injection wells should be arranged perpendicular to the direction of flow. However, this method is not suitable if there is a long distance between the injection and observation wells. The equation used for the tracer test analysis using TRAC ver. 1.7 [61] is:
  C ( x ,   t ) =   M 0 A ω   [ 2 4 π D L T exp ( ( u t ) 2 4 D L t )     μ 2 D L exp ( u x D L ) e r f c ( ( x   +   u t ) 4 D L T ) ] exp ( λ t )  
where M is the injection tracer’s mass (kg), u is the pore water velocity (m/s), A is the groundwater flow section (m2), ω is the effective porosity (0 < ω < 1), DL is the longitudinal dispersion coefficient (m2/s), erfc is the error function, t is time, r is the radial distance (m), x is the distance between the injection and monitoring wells along the x-axis (m), and λ is the decay constant (1/s).

2.3. Groundwater Modeling

Groundwater modeling simulates the groundwater flow in porous or fractured saturated media in order to estimate or predict the groundwater quantity and/or quality in aquifers in steady and transient states. The world’s most popular software for groundwater modeling is Visual MODFLOW. It uses the finite difference method that is equipped by MODFLOW [62] and simulates the solute transport of non-reactive and reactive chemical constituents in 3D using the MT3DMS [63] and RT3D [64] packages. The MT3DMS package can precisely describe the transport of miscible contaminants in groundwater and takes into account advection, dispersion, and diffusion, as well as some basic chemical reactions. Various boundary conditions and external sources or sinks can be applied, and it is able to interpret the spatial concentrations of bio-chemical reactive multi-species pollutants such as BTEX (Benzene, Toluene, Ethylbenzene, Xylene) and chlorinated organic materials. The basic equation for non-reactive and reactive solute transport is:
  C k t =   x i ( D i , j C k x i )     x i ( v i C k )   +   q s Φ C s k + k = 1 N r c , k = 1 ,   2 ,   3    
where Ck is the concentration of pollutants in the groundwater; r is the species reaction term; t is time (T); x i is the coordinate distance (L); Di is the hydraulic dispersion coefficient (L2T−1); v i is the void flow rate (LT−1); q s is the volumetric flux, with a positive value for sources and a negative value for sinks; C s k is the concentration of sources and sinks; and r c is the chemical reaction term (adsorption, biodegradation, etc.).
Particle tracking is a simple method that is commonly used to estimate the groundwater flow and contaminant pathways and to approximate the transport of contaminants by advection without using the complex MT3DMS and RT3D packages. Traditionally, particle tracking by the first-order Euler algorithm is given by:
  x n + 1 = x n   +   Δ t R v x ( x n , y n ,   z n ) ,   y n + 1 = y n   +   Δ t R v y ( x n ,   y n ,   z n ) ,   z n + 1 = z n   +   Δ t R v z ( x n ,   y n ,   z n )  
where x n + 1 , y n + 1 , and z n + 1 , as well as x n , y n , and z n are the coordinates x, y, and z of the particles at the (n + 1)th step and the nth step, respectively; v x , v y , and v z are the linear velocities of the particles at the nth step; R is the retardation factor of the adsorption reaction in the solute transport equation; and Δt is the time step range of the solute behavior, which is usually determined by Equation (5), the Courant condition:
  | Δ t |     γ c R × M I N ( Δ t v x ,   Δ y v y ,   Δ z v z )  
where γ c is the Courant number, which is the number of cells and represents the transport distance of one particle in any direction during one time step of forward particle tracking for positive Δt values and backward particle tracking for negative Δt values. Forward particle tracking chiefly predicts the range across which the contaminant will spread from a known source, while backward particle tracking is used to identify the source and transport pathway of a contaminant when the source is unidentified.
In this study, Visual MODFLOW software was used to simulate the groundwater flow and the acesulfame transport in the study area.

3. Results

3.1. Water Quality

The physicochemical component pH displayed a mean values of 6.01 with minimum (5.43) at BS-M2 and the maximum (6.86) at BS-M3. The physicochemical components EC displayed a mean value of 468 μS/cm with the minimum (172 μS/cm) at BS-3 and maximum (919 μS/cm) at BS-9. The groundwater types were classified into the [Ca-HCO3] type, which was dominant in the residential zone, and the [Na-HCO3] type, which was dominant in the agricultural zone (Table 5).
Groundwater samples were collected from the eight wells BS-3, BS-5, BS-6, BS-8, BS-9, BS-M1, BS-M2, and BS-M3 for analyses in May, August, and October 2017. Acesulfame was detected in all eight samples from the first, second, and third surveys, with mean values of 102 ng/L (8.1–345 ng/L), 577 ng/L (17.3–2695 ng/L), and 214 ng/L (10.3–764 ng/L), respectively (Table 6). Cyclamate was detected in all eight samples from the first survey with a mean value of 249 ng/L (13.5–860 ng/L); however, it was only detected in the BS-M1 sample from the second survey and in the BS-9 and BS-M1 samples from the third survey. Saccharin was detected in some of the samples from the first survey, while it was detected in all eight samples from the second and third surveys with mean values of 257 ng/L (68.1–1156 ng/L) and 134 ng/L (12.7–236 ng/L), respectively.
For comparing with the acesulfame concentrations (8.1–2690 ng/L), the nitrate-N concentrations (1.8–29.2 mg/L) were examined, which were detected in 29 wells during October 2014 [43]. Both acesulfame and nitrate-N showed similarly high concentrations in the residential zone, moving towards the stream in the northern region (Figure 2). Spatially, a higher concentration of acesulfame was similar to the nitrate-nitrogen concentration [65]; this showed levels in excess of Korea’s drinking water criterion for nitrate-N concentration (10 mg/L) at 17 of the 29 wells, with a mean value of 12.72 mg/L and a median value of 12.6 mg/L (Table 7).

3.2. Hydrogeological Results

For the grain size analysis, soil samples were collected from within a 1-m radius of eight wells. After removing more than 30 cm of the topsoil, the soil, except for the organic matter, was collected. Using the four empirical equations given in Table 5 and based on a de between 0.047 mm (BS-8) and 0.661 mm (BS-6) with a mean of 0.194 mm, we obtained mean K values that ranged from 2.678 × 10−5 m/s (BS-8) to 1.124 × 10−2 m/s (BS-6), with a mean K value of 1.65 × 10−3 m/s (Figure 3, Table 8).
According to the slug test analysis at the five wells BS-3, BS-5, BS-6, BS-M2, and BS-M3, the maximum K values were obtained at the BS-5 well: 4.353 × 10−3 m/s according to the Hvorslev method and 5.31 × 10−3 m/s according to the Bouwer and Rice method, with an average of 4.92 × 10−3 m/s. On the other hand, at the BS-M3 well, the minimum K values were estimated to be 1.991 × 10−6 m/s using the Hvorslev method and 2.067 × 10−6 m/s using the Bouwer and Rice method, with an average of 2.029 × 10−6 m/s (Table 9).
A multi-well tracer test was performed in order to evaluate the connectivity between the two test wells and the diffusion coefficient in the geological formations. A NaCl injection solution was instantaneously injected into the injection well (BS-6), and changes in the EC and salinity were monitored at the observation well (BS-8) using a CTD (Conductivity, Temperature, Depth)-diver; the wells were ~110 m apart. Changes in the EC and salinity were detected 72 min after the injection, and they persisted for 1032 min. The maximum EC was 0.004 mS/cm, and the maximum salinity was 0.003 PSU (Practical Salinity Unit). After 1104 min, the EC and salinity rapidly decreased and returned to the background concentration. As a result of the tracer test analysis, the longitudinal dispersivity and solute velocity were estimated to be 12.3 m and 2.76 × 10−3 m/s, respectively.

3.3. Groundwater Modeling

In order to simulate the behavior of groundwater pollutants, a numerical model was generated according to the following steps: data acquisition and hydrogeological characterization, conceptualization, calibration, and prediction. First, the steady-state groundwater flow was modeled, then the transport of acesulfame in the groundwater was included using the MT3DMS module in Visual MODFLOW ver. 2014.1. This revealed the groundwater pollution in the recharge, middle, and discharge areas.

3.3.1. The Conceptual Model and Input Parameters

The conceptual model was built with 43,776 (= 152 × 72 × 4) cells. The geological profile was based on data from right wells (Table 10): The study area consists of a landfill layer, agricultural soil, alluvial deposit, weathered soil, weathered rock, and soft rock. Based on the geological profile, the model’s layers were defined as Layer 1 (a landfill layer, agricultural soil, an alluvial deposit, and weathered soil), Layer 2 (weathered rock), Layer 3 (soft rock), and Layer 4 (fresh rock). Irrigation drainage occurs through a drainage channel from the southern hilly zone to the northeastern region.
The external boundary was set to be the Neman no-flow boundary condition. For the internal boundary conditions, the stream boundary was set to be the drainage channel in the southern hilly zone (that runs along the village and discharge area (Figure 4a)). Each stream boundary’s thickness, width, and stream conductance was 0.2 m, 0.5 m, and 1 × 10−5 m/s, respectively. The observed groundwater levels at the eight observation wells (that were converted to an average mean sea level) were used as the inputs for the model.
Based on the results of the slug test, the K, storativity, effective porosity, and recharge rate values were input differently for the hilly zone, the residential zone, and the agricultural zone (Figure 4b, Table 11). The longitudinal dispersivity (that was obtained by the tracer test analysis) of 12.7 m, as well as the concentration of acesulfame on 22 May 2017, was also entered into the model.

3.3.2. The Calibration and Steady-State Simulation

The simulated water levels were calibrated to the observed levels by varying the values of K over the range obtained in the slug test and by the result of the steady-state simulation (Table 12). The groundwater level data from the BS-3, BS-M1, BS-M2, and BS-M3 wells were excluded from the calibration because those wells were located in a residential zone and were being actively pumped. The calibration was completed with a correlation coefficient of 0.903 and a confidence interval within 95%; the residuals from the calibration ranged from −3.4–0.63 m, with an average of −1.28 m (Figure 5a, Table 13).
The steady-state model showed that groundwater principally moves from south to north along the stream boundary and then converges to the northeast. The groundwater levels were 38–28 m a m.s.l. The mean groundwater velocity through the forward particle tracking was 5.2 × 10−6 m/s; the particles were predicted to flow out through the stream boundary to the north and northeast.

3.3.3. Prediction of the Acesulfame’s Spread

According to the particle tracking, the mean velocity of the groundwater over the entire study area is 5.2 × 10−6 m/s. Through forward particle tracking, we see that particles flow through the stream boundary to the north and northeast discharge boundaries (Figure 6). The particles from wells BS-5, BS-8, and BS-3 flow fast (4.267 × 10−8–1.559 × 10−6 m/s), and the particles from the wells BS-M1 and BS-M2 (which are located in the residential zone) have a low velocity (2.158 × 10−9–2.755 × 10−7 m/s). Forward particle tracking was carried out for the central residential zone (where acesulfame is probably derived); hence, the groundwater flow from the residential zone to the stream boundary was estimated. In the eastern part of the residential zone, the particles usually flow into the drain conduit with a relatively high velocity. However, in the western part, where many houses and the town hall are concentrated together in the village, the particles move slowly. Therefore, pollution in the western drain channel that heads east is expected to be more serious than the pollution along the eastern drain channel (Figure 7).
Subsequently, a 10-year transient modeling with the MT3DMS module was performed to estimate the areal distribution of the acesulfame concentration. The concentration of acesulfame at each well was given an initial value of 102.1 ng/L, which was obtained from the water quality analysis of the samples in the first survey in 2017. Using the MT3DMS model, we observed that, in Years 1–3, acesulfame actively spread from the BS-9 and BS-M3 wells and outflowed through the stream boundary (Figure 8). Then, the concentration of acesulfame attained equilibrium downstream of the stream boundary and started to decrease in the agricultural region after approximately five years. After 10 years, the acesulfame in the agricultural zone, and in the hilly zone, was mostly discharged through the stream boundary, while the acesulfame in the central residential zone remained close to the initial state.

4. Discussion

A typical indicator of pollution in groundwater is nitrate-N, which commonly derives from livestock wastewater, manure, domestic wastewater, pesticides, and agricultural fertilizers, mainly in agricultural and livestock areas. The seasonality of nitrogen concentrations in groundwater occurs in agricultural areas in Europe, showing the highest nitrogen concentrations in spring after winter freezing with a small vegetation cover, as well as in the areas of the high level of nitrogen application [32]. According to the NIER’s 2016 survey in Korea [65], sulfamethoxazole (antibiotics), oxfendazole (insect repellent), and carbofuran (pesticides) were only detected in agricultural and livestock areas and may enter the groundwater as a result of agricultural and animal farm activities.
Acesulfame is recognized as a pollutant in groundwater and surface water that originates from sewage inflow [12,13]. The concentration of acesulfame in the groundwater in urban areas was found to be between 3.5- and 67-times higher than that in agricultural and livestock areas [65]. Compared to nitrate-N, acesulfame is highly persistent and has a low removal rate at the treatment site [66]. At present, the toxicity of acesulfame is not clear, but acesulfame may be an appropriate indicator of groundwater contamination originating from man-made sewage.
The comparison of the nitrate-N water quality analysis with the acesulfame water quality analysis, along with the groundwater modeling, verified that acesulfame is a potential indicator for the contamination of groundwater in rural areas in Korea. Hence, acesulfame has been confirmed to be a very useful indicator for the detection of sources of man-made pollution. Acesulfame may be a useful tool to design regulations, management measures, and policies for groundwater in rural areas where people use large amounts of artificial sweeteners.

5. Conclusions

In this study, acesulfame was confirmed to be the most common artificial sweetener that contaminates groundwater. Acesulfame showed high concentrations (8.1–2690 ng/L) in the residential zone, as it is known to be a pollutant in groundwater and surface water that originates from sewage inflow. Besides this, the nitrate-N concentrations (1.8–29.2 mg/L) displayed a similar distribution pattern to the acesulfame concentrations. By the MT3DMS and particle tracking modeling, acesulfame actively spreads from the BS-9 and BS-M3 wells in the residential zone and outflows through the stream boundary, indicating that it discharges in the agricultural and hilly zones.
The acesulfame concentration in the groundwater was evaluated using particle tracking in Visual MODFLOW with the MT3DMS module. Forward particle tracking showed that the groundwater flows relatively rapidly in the agricultural land, flows slowly in the western residential zone, and finally flows through the stream boundary to the northern region. According to the MT3DMS simulation, acesulfame spreads across the model domain over 1–3 years, and the concentration gradually decreases after ~5 years. Over the 10-year simulation period, acesulfame escaped through the stream boundary in the northern part of the study area, while it accumulated in the residential zone and downstream of the stream boundary. Therefore, it is anticipated that the groundwater’ acesulfame contamination level may be higher in the western drainage channel than in the eastern drainage channel.

Author Contributions

Conceptualization, C.-M.L. and S.-Y.H.; methodology, S.-Y.H. and J.-E.O.; software, C.-M.L. and S.-M.Y.; validation M.K. and H.K.K.; formal analysis, C.-M.L., J.-E.O. and S.-M.Y.; investigation, C.-M.L. and S.-M.Y.; resources, M.K.; data curation, H.K.K.; writing, original draft preparation, C.-M.L. and S.-Y.H.; writing, review and editing, S.-Y.H.; visualization, C.-M.L.; supervision, S.-Y.H.; project administration, S.-Y.H.; funding acquisition, S.-Y.H.

Funding

This research was funded by National Research Foundation of Korea (NRF) under the Ministry of Science and ICT (NRF-2017R1A2B2009033), as well as by the National Institute of Environment Research (NIER), funded by the Ministry of Environment (MOE) of Korea (NIER-2017-03-02-019).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Subedi, B.; Kannan, K. Fate of artificial sweeteners in wastewater treatment plants in New York State, USA. Environ. Sci. Technol. 2014, 48, 13668–13674. [Google Scholar] [CrossRef] [PubMed]
  2. Lange, F.T.; Scheurer, M.; Brauch, H.J. Artificial sweeteners—A recently recognized class of emerging environmental contaminants: A review. Anal. Bioanal. Chem. 2012, 403, 2503–2518. [Google Scholar] [CrossRef] [PubMed]
  3. Van Stampvoort, D.R.; Roy, J.W.; Brown, S.J.; Bickerton, G. Artificial sweeteners as potential tracers in groundwater in urban environments. J. Hydrol. 2011, 401, 126–133. [Google Scholar] [CrossRef]
  4. Loos, R.; Gawlik, B.M.; Boettcher, K.; Locoro, R.; Contini, S.; Bidoglio, G. Sucralose screening in European surface waters using a solid-phase extraction-liquid chromatography-triple quadrupole mass spectrometry method. J. Chromatogr. A 2009, 1216, 1126–1131. [Google Scholar] [CrossRef] [PubMed]
  5. Buerge, I.J.; Buser, H.R.; Kahle, M.; Muller, M.D.; Poiger, T. Ubiquitous occurrence of the artificial sweetener acesulfame in the aquatic environment: An ideal chemical marker of domestic wastewater in groundwater. Environ. Sci. Technol. 2009, 43, 4381–4385. [Google Scholar] [CrossRef] [PubMed]
  6. Arbeitsgemeinschaft Wasserwerke Bodensee Rhein (AWBR). 42. Jahresbericht 2010; Geschäftsstelle der Arbeitsgemeinschaft Wasserwerke Bodensee-Rhein (AWBR): Freiburg, Germany, 2010. [Google Scholar]
  7. Mead, R.N.; Morgan, J.B.; Avery, G.B.; Kieber, R.J.; Kirk, A.M.; Skrabal, S.A.; Willey, J.D. Occurrence of the artificial sweetener sucralose in coastal and marine waters of the United States. Mar. Chem. 2009, 116, 13–17. [Google Scholar] [CrossRef]
  8. Storck, F.R.; Skark, C.; Remmier, F.; Brauch, H.J. Environmental fate and behavior of acesulfame in laboratory experiments. Water Sci. Technol. 2016, 74, 2832–2842. [Google Scholar] [CrossRef] [PubMed]
  9. Mawhinney, D.B.; Young, R.B.; Vanderford, B.J.; Borch, T.; Snyder, S.A. Artificial sweetener sucralose in US drinking water systems. Environ. Sci. Technol. 2011, 45, 8716–8722. [Google Scholar] [CrossRef] [PubMed]
  10. Soh, L.; Connors, K.A.; Brooks, B.W.; Zimmerman, J. Fate of sucralose through environmental and water treatment processes and impact on plant indicator species. Environ. Sci. Technol. 2011, 45, 1363–1369. [Google Scholar] [CrossRef] [PubMed]
  11. Dep, P.; Debnath, P.; Denis, A.F.; Lepcha, O.T. Variability of soil physicochemical properties at different agroecological zones of Himalayan region: Sikkim, India. Environ. Dev. Sustain. 2018, 1–19. [Google Scholar] [CrossRef]
  12. Tran, N.H.; Hu, J.; Li, J.; Ong, S.L. Suitability of artificial sweeteners as indicators of raw wastewater contamination in surface water and groundwater. Water Res. 2014, 48, 443–456. [Google Scholar] [CrossRef] [PubMed]
  13. Wolf, L.; Zwiener, C.; Zemann, M. Tracking artificial sweeteners and pharmaceuticals introduced into urban groundwater by leaking sewer networks. Sci. Total Environ. 2012, 430, 8–19. [Google Scholar] [CrossRef] [PubMed]
  14. Morasch, B. Occurrence and dynamics of micropollutants in a karst aquifer. Environ. Pollut. 2013, 173, 133–137. [Google Scholar] [CrossRef] [PubMed]
  15. Cabeza, Y.; Candela, L.; Ronen, D.; Teijon, G. Monitoring the occurrence of emerging contaminants in treated wastewater and groundwater between 2008 and 2010. J. Hazard. Mater. 2012, 239, 32–39. [Google Scholar] [CrossRef] [PubMed]
  16. López-Serna, R.; Jurado, A.; Vázquez-Suñé, E.; Carrera, J.; Petrović, M.; Barceló, D. Occurrence of 95 pharmaceuticals and transformation products in urban groundwaters underlying the metropolis of Barcelona, Spain. Environ. Pollut. 2013, 174, 305–315. [Google Scholar] [CrossRef] [PubMed]
  17. Schaider, L.A.; Rudel, R.A.; Ackerman, J.M.; Dunagan, S.C.; Brody, J.G. Pharmaceuticals, perfluorosurfactants, and other organic wastewater compounds in public drinking water wells in a shallow sand and gravel aquifer. Sci. Total Environ. 2014, 468, 384–393. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Tong, L.; Huang, S.; Wang, Y.; Liu, H.; Li, M. Occurrence of antibiotics in the aquatic environment of Jianghan Plain, central China. Sci. Total Environ. 2014, 497, 180–187. [Google Scholar] [CrossRef] [PubMed]
  19. National Institute of Environmental Research. Investigation and Survey of Residual Drug Substance Analysis Method (IV); National Institute of Environmental Research: Incheon, Korea, 2011; pp. 1–358.
  20. Animal and Plant Quarantine Agency. Investigation of Veterinary Antimicrobial Drug Residues and Their Behavior in Animal Farming Environment; Animal and Plant Quarantine Agency: Gimcheon, Korea, 2010; pp. 683–704.
  21. Future Market Insights (FMI). Animal Healthcare Market: Global Industry Analysis and Opportunity Assessment 2015–2020; Future Market Insights: London, UK, 2015. [Google Scholar]
  22. Hu, X.; Zhou, Q.; Luo, Y. Occurrence and source analysis of typical veterinary antibiotics in manure, soil, vegetables and groundwater from organic vegetable bases, northern China. Environ. Pollut. 2010, 158, 2992–2998. [Google Scholar] [CrossRef] [PubMed]
  23. Yao, L.; Wang, Y.; Tong, L.; Deng, Y.; Li, Y.; Gan, Y.; Guo, W.; Dong, C.; Duan, Y.; Zhao, K. Occurrence and risk assessment of antibiotics in surface water and groundwater from different depths of aquifers: A case study at Jianghan Plain, central China. Ecotoxicol. Environ. Saf. 2017, 135, 236–242. [Google Scholar] [CrossRef] [PubMed]
  24. Balzer, F.; Zühlke, S.; Hannappel, S. Antibiotics in groundwater under locations with high livestock density in Germany. Water Sci. Technol. Water Supply 2016, 16, 1361–1369. [Google Scholar] [CrossRef]
  25. World Health Organization (WHO). Pharmaceuticals in Drinking-Water; World Health Organization: Geneva, Switzerland, 2011; pp. 1–52. [Google Scholar]
  26. Kümmerer, K. Antibiotics in the aquatic environment—A review—Part I. Chemosphere 2009, 75, 417–434. [Google Scholar] [CrossRef] [PubMed]
  27. U.S. Environmental Protection Agency. Contaminants of Emerging Concern including Pharmaceuticals and Personal Care Products; U.S. Environmental Protection Agency, Office of Water: Washington, DC, USA, 2005.
  28. Bartelt-Hunt, S.; Snow, D.D.; Damon-Powell, T.; Miesbach, D. Occurrence of steroid hormones and antibiotics in shallow groundwater impacted by livestock waste control facilities. J. Contam. Hydrol. 2011, 123, 94–103. [Google Scholar] [CrossRef] [PubMed]
  29. Batt, A.L.; Snow, D.D.; Aga, D.S. Occurrence of sulfonamide antimicrobials in private water wells in Washington County, Idaho, USA. Chemosphere 2006, 64, 1963–1971. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Watanabe, N.; Bergamaschi, B.A.; Loftin, K.A.; Meyer, M.T.; Harter, T. Use and environmental occurrence of antibiotics in freestall dairy farms with manured forage fields. Environ. Sci. Technol. 2010, 44, 6591–6600. [Google Scholar] [CrossRef] [PubMed]
  31. Barnes, K.K.; Kolpin, D.W.; Meyer, M.T.; Thurman, E.M.; Furlong, E.T.; Zaugg, S.D.; Barber, L.B. Water-Quality Data for Pharmaceuticals, Hormones, and Other Organic Wastewater Contaminants in U.S. Streams, 1999–2000; U.S. Environmental Protection Agency: Iowa City, IA, USA, 2002.
  32. Lawniczak, A.E.; Zbierska, J.; Nowak, B.; Achtenberg, K.; Grześkowiak, A.; Kanas, K. Impact of agriculture and land use on nitrate contamination in groundwater and running waters in central-west Poland. Environ. Monit. Assess. 2016, 188, 172. [Google Scholar] [CrossRef] [PubMed]
  33. Appelo, C.A.J.; Postma, D. Geochemistry, Groundwater and Pollution, 2nd ed.; A.A. Balkema Publisher: Amsterdam, The Netherlands, 1996; pp. 271–279. [Google Scholar]
  34. Lolan, B.T.; Ruddy, B.C.; Hitt, K.J.; Helsel, D.R. Risk of nitrate in groundwaters of the United States-A national perspective. Environ. Sci. Technol. 1997, 31, 2229–2236. [Google Scholar]
  35. Pastén-Zapara, E.; Ledesma-Ruiz, R.; Harter, T.; Ramírez, A.I.; Maglknecht, J. Assessment of sources and fate of nitrate in shallow groundwater of an agricultural area by using a multi-tracer approach. Sci. Total Environ. 2014, 470–471, 855–864. [Google Scholar] [CrossRef] [PubMed]
  36. Kim, K.H.; Yun, S.T.; Chae, G.T.; Choi, B.Y.; Kim, S.O.; Kim, K.J.; Kim, H.S.; Lee, C.W. Nitrate contamination of alluvial groundwaters in the Keum river watershed area: Source and behaviors of nitrate, and suggestion to secure water supply. J. Eng. Geol. 2002, 12, 471–484. [Google Scholar]
  37. Cheong, J.Y.; Hamm, S.Y.; Lee, J.H.; Lee, K.S.; Woo, N.C. Groundwater nitrate contamination and risk assessment in an agricultural area, South Korea. Environ. Earth Sci. 2012, 66, 1127–1136. [Google Scholar] [CrossRef]
  38. Kim, J.W.; Cho, M.Y.; Jee, B.; Park, M.A.; Kim, N.Y. Administration and use of aquaculture drugs in Korea. J. Fish Pathol. 2014, 27, 67–75. [Google Scholar] [CrossRef] [Green Version]
  39. Yoon, Y.; Ryu, J.; Oh, J.; Choi, B.G.; Snyder, S.A. Occurrence of endocrine disrupting compounds, pharmaceuticals, and personal care products in the Han River (Seoul, South Korea). Sci. Total Environ. 2010, 408, 636–643. [Google Scholar] [CrossRef] [PubMed]
  40. National Institute of Environmental Research. Study on the Emission Source and Behavior of Medicinal Substances in the Environment; National Institute of Environmental Research: Incheon, Korea, 2011; pp. 1–282.
  41. Ji, K.; Kim, S.; Han, S.; Seo, J.; Lee, S.; Park, Y.; Choi, K. Risk assessment of chlortetracycline, oxytetracycline, sulfamethazine, sulfathiazole, and erythromycin in aquatic environment: Are the current environmental concentrations safe? Ecotoxicology 2012, 21, 2031–2050. [Google Scholar] [CrossRef] [PubMed]
  42. Creed, I.F.; Band, L.E.; Foster, N.W.; Morrison, I.K.; Nicolson, J.A.; Semkin, R.S.; Jeffries, D.S. Regulation of nitrate-N release from temperate forests: A test of the N flushing hypothesis. Water Resour. Res. 1996, 32, 3337–3354. [Google Scholar] [CrossRef]
  43. National Institute of Environmental Research. A study of Optimal Management Measures for Securing Safe Groundwater; National Institute of Environmental Research: Incheon, Korea, 2015; pp. 1–302.
  44. Chai, J.C.; Agung, P.M.A.; Hino, T.; Igaya, T.; Cater, J.P. Estimating hydraulic conductivity from piezocone soundings. Géotechnique 2011, 61, 699–708. [Google Scholar] [CrossRef]
  45. Chai, J.C.; Chanmee, N. A modified method for estimating permeability of clayey soils based on piezocone sounding results. Can. Geotech. J. 2018, 55, 1268–1281. [Google Scholar] [CrossRef]
  46. Shen, S.L.; Wang, J.P.; Wu, H.N.; Xu, Y.S.; Ye, G.L.; Yin, Z.Y. Evaluation of hydraulic conductivity for both marine and deltaic deposit based on piezocone test. Ocean Eng. 2015, 110, 174–182. [Google Scholar] [CrossRef]
  47. Wu, Y.X.; Shen, S.L.; Yuan, D.J. Characteristics of dewatering induced drawdown curve under blocking effect of retaining wall in aquifer. J. Hydrol. 2016, 539, 554–566. [Google Scholar] [CrossRef]
  48. Wu, Y.X.; Shen, J.S.; Cheng, W.C.; Hino, T. Semi-analytical solution to pumping test data with barrier, wellbore storage, and partial penetration effects. Eng. Geol. 2017, 226, 44–51. [Google Scholar] [CrossRef]
  49. Cheong, J.Y.; Hamm, S.Y.; Kim, H.S.; Ko, E.J.; Yang, K.; Lee, J.H. Estimating hydraulic conductivity using grain-size analyses, aquifer tests, and numerical modeling in a riverside alluvial system in South Korea. Hydrogeol. J. 2008, 16, 1129–1143. [Google Scholar] [CrossRef]
  50. Du, Y.J.; Jiang, N.J.; Shen, S.L.; Jin, F. Experimental investigation of influence of acid rain on leaching and hydraulic characteristics of cement-based solidified/stabilized lead contaminated clay. J. Hazard. Mater. 2012, 30, 225–226. [Google Scholar] [CrossRef] [PubMed]
  51. Du, Y.J.; Fan, R.D.; Liu, S.Y.; Reddy, K.R.; Jin, F. Workability, compressibility and hydraulic conductivity of zeolite-amended clayey soil/calcium-bentonite backfills for slurry-trench cutoff walls. Eng. Geol. 2015, 195, 258–268. [Google Scholar] [CrossRef]
  52. Kasenow, M. Determination of Hydraulic Conductivity from Grain Size Analysis; Water Resources Publications LLC: Denver, CO, USA, 2002; pp. 47–84. [Google Scholar]
  53. Hazen, A. Some Physical Properties of Sands and Gravels: With Special Reference to Their Use in Filtration; Annual Report; Massachusetts State Board of Health: Boston, MA, USA, 1892; pp. 539–556.
  54. Beyer, W. On the Determination of Hydraulic Conductivity of Gravels and Sands from Grain-Size Distribution; Wasserwirtschaft-Wassertechnik: Berlin, Germany, 1964; Volume 14, pp. 165–169. [Google Scholar]
  55. Kozeny, J. Ueber kapillare leitung des wassers im boden. Sitzungsber Acad. Wiss. Wien 1927, 136, 271–306. [Google Scholar]
  56. Vukovic, M.; Soro, A. Determination of hydraulic conductivity of porous media from grain-size composition; Water Resources Publications: Littleton, CO, USA, 1992; pp. 1–83. [Google Scholar]
  57. Hvorslev, M.J. Time Lag and Soil Permeability in Ground Water Observations; US Army Corps of Engineers, Waterway Experimentation Station: Vicksburg, MI, USA, 1951; pp. 1–50.
  58. Bouwer, H.; Rice, R.C. A slug test for determining hydraulic conductivity of unconfined aquifers with completely or partially penetrating wells. Water Resour. Res. 1976, 12, 423–428. [Google Scholar] [CrossRef]
  59. Cooper, H.H.; Bredehoeft, J.D.; Papadopulos, S.S. Response of a finite-diameter well to an instantaneous charge of water. Water Resour. Res. 1967, 3, 263–269. [Google Scholar] [CrossRef]
  60. Hubbard, E.F.; Kilpatrick, F.A.; Martens, L.A.; Wilson, J.F. Measurement of Time of Travel and Dispersion in Streams by Dye Tracing—Techniques of Water-Resources Investigations of the US Geological Survey; U.S. Geological Survey: Washington, DC, USA, 1982; pp. 1–52.
  61. Gurierrez, A.; Klinka, T.; Thiéry, D.; Buscarlet, E.; Binet, S.; Jozja, N.; Défarge, C.; Leclerc, B.; Fécamp, C.; Ahumada, Y.; Elsass, J. TRAC, a collaborative computer tool for tracer-test interpretation. EPJ Web Conf. 2013, 50, 1–8. [Google Scholar] [CrossRef] [Green Version]
  62. McDonald, M.G.; Harbaugh, A.W. A Modular Three-dimensional Finite-Difference Ground-water Flow Model; Techniques of Water-Resources Investigations of the United States Geological Survey, Book 6, Chapter A1; U.S. Geological Survey: Washington, DC, USA, 1988.
  63. Zheng, P.C.; Wang, P. MT3DMS: A Modular Three-Dimensional Multispecies Transport Model for Simulation of Advection, Dispersion, and Chemical Reactions of Contaminants in Groundwater Systems; Documentation and User’s Guide; US Army Corps of Engineers: Washington, DC, USA, 1999; pp. 1–219.
  64. Clement, T.P. RT3D—A Modular Computer Code for Simulating Reactive Multi-Species Transport in 3-Dimensional Groundwater Aquifer; U.S. Department of Energy: Washington, DC, USA, 1997; pp. 1–59.
  65. National Institute of Environmental Research. A study on the Tracking Technique of Groundwater Contamination in Korea (I); National Institute of Environmental Research: Incheon, Korea, 2016; pp. 1–101.
  66. Buerge, I.J.; Buster, H.R.; Kahle, M.; Müller, M.D.; Poiger, T. The artificial sweetener acesulfame as marker of domestic wastewater in groundwater. In Proceedings of the Groundwater Quality Management in a Rapidly Changing World, 7th International Groundwater Quality Conference, Zurich, Switzerland, 13–18 June 2010; Volume 342, pp. 103–106. [Google Scholar]
Figure 1. The location of the study area and the wells that were investigated for acesulfame and nitrate-N.
Figure 1. The location of the study area and the wells that were investigated for acesulfame and nitrate-N.
Water 10 01731 g001
Figure 2. The distribution of the concentration of (a) nitrate-N and (b) acesulfame.
Figure 2. The distribution of the concentration of (a) nitrate-N and (b) acesulfame.
Water 10 01731 g002
Figure 3. The results of the grain size analysis. (a) Five single wells and (b) three multi wells.
Figure 3. The results of the grain size analysis. (a) Five single wells and (b) three multi wells.
Water 10 01731 g003
Figure 4. The (a) boundary conditions and (b) hydraulic conductivity values in the hilly zone (green color), the residential zone (blue color), and the agricultural zone (blue-green color).
Figure 4. The (a) boundary conditions and (b) hydraulic conductivity values in the hilly zone (green color), the residential zone (blue color), and the agricultural zone (blue-green color).
Water 10 01731 g004
Figure 5. The (a) calculated versus observed levels and (b) groundwater level contours in the steady-state.
Figure 5. The (a) calculated versus observed levels and (b) groundwater level contours in the steady-state.
Water 10 01731 g005
Figure 6. Forward particle tracking of the movement of particles from the wells.
Figure 6. Forward particle tracking of the movement of particles from the wells.
Water 10 01731 g006
Figure 7. Particle tracking from the stream boundary.
Figure 7. Particle tracking from the stream boundary.
Water 10 01731 g007
Figure 8. The MT3DMS model of the areal concentration of acesulfame after (a) 1 year, (b) 3 years, (c) 5 years, (d) 7 years, and (e) 10 years.
Figure 8. The MT3DMS model of the areal concentration of acesulfame after (a) 1 year, (b) 3 years, (c) 5 years, (d) 7 years, and (e) 10 years.
Water 10 01731 g008
Table 1. The concentration of artificial sweeteners in groundwater in some countries.
Table 1. The concentration of artificial sweeteners in groundwater in some countries.
E.U. [4]Switzerland [5]Germany [6]U.S. [7]Serbia [8]
Acesulfame-0.2–6.90.4–1.1-0.09–0.31
Cyclamate--0.01–0.13-0.09–0.13
Saccharin--0.01–3.8-0.006–0.08
Sucralose<0.01–0.092<LOD–0.6<LOQ–0.10.0008–2.0<LOQ
LOD: Limit Of Detection, LOQ: Limit Of Quantification.
Table 2. The laboratory analysis of the artificial sweeteners and the triclosan (negative mode).
Table 2. The laboratory analysis of the artificial sweeteners and the triclosan (negative mode).
Conditions
LCColumnZORBAX Eclipse XDB-C18 (4.6 × 150 mm, 3.5 μm)
Column temp.40 °C
Mobile-phaseA: 5 mM ammonium acetate in methanol
B: 5 mM ammonium acetate in water
Flow rate0.4 mL/min
Gradient programMin01416.917
B (%)6060100160
Injection vol.10 μL
MS/MSIonization modeNegative ion electrospray
Gas temp.250 °C
Gas flow rate10 L/min
Sheath gas temp.380 °C
Sheath gas flow rate10 L/min
Capillary voltage4000 V
Nebulizer pressure40 psi
Table 3. The sieve numbers and mesh sizes.
Table 3. The sieve numbers and mesh sizes.
Sieve No.810204060100140200
Mesh size (mm)2.3620.840.4250.250.150.1050.074
Table 4. The empirical formulas for the grain size analysis.
Table 4. The empirical formulas for the grain size analysis.
Authorβυ(n)deGrain Range
Hazen [53]6 × 10−4(1 + 10(n − 0.26))d100.1 mm < de < 3 mm, U < 5
Beyer [54]6 × 10−4log(500/C)1d100.06 mm < de < 0.6 mm, 1 < U < 20
Kozeny [55]8.3 × 10−3n3/(1 − n)2d10Coarse sand
Sauerbrei [56]3.75 × 10−3n3/(1 − n)2d17Sand and sandy clay, de < 0.5 mm
U: d 60 d 10 .
Table 5. The physicochemical properties and the groundwater types in the residential and agricultural zones.
Table 5. The physicochemical properties and the groundwater types in the residential and agricultural zones.
Well No.TpHECWater Type
°C-µS/cm
Residential zoneBS-31st16.05.56172Na-HCO3
2nd19.15.98584Na-Cl
3rd17.26.03231Na-HCO3
BS-M11st18.06.12393Ca-HCO3
2nd19.86.66555Ca-HCO3
3rd16.36.34615Ca-HCO3
BS-M21st17.315.43272Ca-HCO3
2nd18.755.88242Ca-HCO3
3rd16.835.94277Ca-NO3
BS-M31st16.386.09575Ca-HCO3
2nd21.56.86799Ca-HCO3
3rd16.26.02649Ca-HCO3
Agricultural zoneBS-51st17.465.47330Na-HCO3
2nd20.55.87539Na-HCO3
3rd18.85.52373Na-HCO3
BS-61st17.946.31288Na-HCO3
2nd19.085.78412Na-HCO3
3rd17.65.77344Ca-HCO3
BS-81st16.95.99428Na-NO3
2nd19.056.08545Na-NO3
3rd17.45.74390Na-NO3
BS-91st17.056.13597Ca-Cl
2nd19.186.22919Ca-Cl
3rd17.16.51713Ca-HCO3
Min. 16.05.43172
Max.21.56.86919
Mean18.06.01468
Table 6. The concentrations of artificial sweeteners and agrichemicals (unit: ng/L).
Table 6. The concentrations of artificial sweeteners and agrichemicals (unit: ng/L).
ConstituentsResidential ZoneAgricultural Zone
BS-3BS-M1BS-M2BS-M3BS-5BS-M6BS-8BS-9
Acesulfame1st25.840.48.134517.911.738.1330
2nd10889.723.5140812017.31562695
3rd13.976410.347.455.711.388.5720
Cyclamate1st21427.113.577.524.6103860674
2ndN.D.367N.D.N.D.N.D.N.D.N.D.N.D.
3rdN.D.337N.D.N.D.441N.D.N.D.3467
Saccharin1st47.412136.911.9N.D.273N.D.N.D.
2nd96.7115611468.111628687.4135
3rd20497.856.517323619912.788.9
N.D.: Not Detected.
Table 7. The concentrations of nitrate-N (unit: mg/L).
Table 7. The concentrations of nitrate-N (unit: mg/L).
WellS-01S-02S-03S-04S-05S-06S-07S-08S-09S-10
Nitrate-N1.83.84.43.94.76.28.44.911.112.6
WellS-11S-12S-13S-14S-15S-16S-17S-18S-19S-20
Nitrate-N15.914.611.213.214.717.129.27.94.58.2
WellS-22S-26S-27S-28S-29S-30S-31S-32S-33
Nitrate-N14.619.220.917.426.124.023.320.25.0
Table 8. The hydraulic conductivity values based on the empirical grain size analysis (unit: m/s).
Table 8. The hydraulic conductivity values based on the empirical grain size analysis (unit: m/s).
Welld60d50d30d17d10 (de)HazenBeyerSauerbreiKozenyMean
BS-31.8450.7220.3010.1390.0799.038 × 10−54.282 × 10−51.595 × 10−41.128 × 10−57.599 × 10−5
BS-52.0861.1790.4830.2870.1834.842 × 10−42.831 × 10−46.763 × 10−46.046 × 10−53.760 × 10−4
BS-64.7583.8212.4642.0260.6616.348 × 10−34.165 × 10−33.364 × 10−27.926 × 10−41.124 × 10−2
BS-80.6710.4440.2000.0900.0473.040 × 10−51.778 × 10−55.555 × 10−53.396 × 10−62.678 × 10−5
BS-92.1171.1620.4120.2180.1121.834 × 10−49.303 × 10−53.880 × 10−42.290 × 10−51.718 × 10−4
BS-M13.2532.5371.4690.5090.3261.451 × 10−39.333 × 10−41.792 × 10−31.621 × 10−41.085 × 10−3
BS-M22.0501.1350.4670.2490.0871.103 × 10−45.214 × 10−55.078 × 10−41.377 × 10−51.710 × 10−4
BS-M32.2971.8880.4190.1590.0584.659 × 10−51.946 × 10−51.738 × 10−45.204 × 10−66.126 × 10−5
Table 9. The hydraulic conductivity values that were estimated using slug tests (unit: m/s).
Table 9. The hydraulic conductivity values that were estimated using slug tests (unit: m/s).
MethodBS-3BS-5BS-6BS-M2BS-M3
Bouwer and Rice1.475 × 10−55.305 × 10−32.074 × 10−32.022 × 10−52.067 × 10−6
Hvorslev2.261 × 10−64.353 × 10−32.622 × 10−31.963 × 10−51.991 × 10−6
Mean1.868 × 10−54.920 × 10−32.348 × 10−31.992 × 10−52.029 × 10−6
Table 10. The well locations, groundwater levels, and acesulfame concentrations.
Table 10. The well locations, groundwater levels, and acesulfame concentrations.
WellXYObserved Head (a m.s.l.)Acesulfame Concentration (ng/L)
BS-3174,790.728398,990.11739.4425.8
BS-5174,904.162399,102.69736.1917.9
BS-6174,840.516399,078.76336.5211.7
BS-8174,916.284399,193.84935.4738.1
BS-9174,800.598399,266.79928.88330
BS-M1174,720.349399,030.52142.8840.4
BS-M2174,760.159399,100.99836.458.1
BS-M3174,826.636399,180.23834.59345
Table 11. The input values of the hydraulic parameters for the different zones.
Table 11. The input values of the hydraulic parameters for the different zones.
ZoneK (m/s)StorativityRecharge (mm/year)Effective Porosity
Hilly7.926 × 10−4–1.710 × 10−21 × 10−52800.15
Residential1.991 × 10−6–6.126 × 10−55 × 10−51200.15
Agricultural2.678 × 10−5–4.920 × 10−31 × 10−52000.15
Table 12. The hydraulic conductivity values from the steady-state model’s calibration.
Table 12. The hydraulic conductivity values from the steady-state model’s calibration.
ZoneK (m/s)
Hilly9.354 × 10−4 m/s
Residential1.354 × 10−5 m/s
Agricultural3.588 × 10−5 m/s
Table 13. The results of the steady-state model’s calibration (unit: a m.s.l., m).
Table 13. The results of the steady-state model’s calibration (unit: a m.s.l., m).
WellObserved HeadCalculated HeadResidual
BS-536.1939.59−3.40
BS-636.5237.73−1.21
BS-835.4736.60−1.13
BS-928.8828.250.63

Share and Cite

MDPI and ACS Style

Lee, C.-M.; Hamm, S.-Y.; Yun, S.-M.; Oh, J.-E.; Kim, M.; Kim, H.K. Indication of Groundwater Contamination Using Acesulfame and Other Pollutants in a Rural Area of Korea. Water 2018, 10, 1731. https://doi.org/10.3390/w10121731

AMA Style

Lee C-M, Hamm S-Y, Yun S-M, Oh J-E, Kim M, Kim HK. Indication of Groundwater Contamination Using Acesulfame and Other Pollutants in a Rural Area of Korea. Water. 2018; 10(12):1731. https://doi.org/10.3390/w10121731

Chicago/Turabian Style

Lee, Chung-Mo, Se-Yeong Hamm, Sul-Min Yun, Jeong-Eun Oh, MoonSu Kim, and Hyun Koo Kim. 2018. "Indication of Groundwater Contamination Using Acesulfame and Other Pollutants in a Rural Area of Korea" Water 10, no. 12: 1731. https://doi.org/10.3390/w10121731

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop