To evaluate the nitrogen pollution load in an aquifer, a water and nitrogen balance analysis was conducted over a thirty-five year period at five yearly intervals. First, we established a two-horizon model comprising a channel/soil horizon, and an aquifer horizon, with exchange of water between the aquifer and river. The nitrogen balance was estimated from the product of nitrogen concentration and water flow obtained from the water balance analysis. The aquifer nitrogen balance results were as follows: 1) In the aquifer horizon, the total nitrogen pollution load potential (NPLP) peaked in the period 1981-1990 at 1800 t·yr-1; following this the NPLP rapidly decreased to about 600 t·yr-1 in the period 2006-2010. The largest NPLP input component of 1000 t·yr-1 in the period 1976-1990 was from farmland. Subsequently, farmland NPLP decreased to only 400 t·yr-1 between 2006 and 2010. The second largest input component, 600 t·yr-1, was effluent from wastewater treatment works (WWTWs) in the period 1986-1990; this also decreased markedly to about 100 t·yr-1 between 2006 and 2010; 2) The difference between input and output in the aquifer horizon, used as an index of groundwater pollution, peaked in the period 1986-1990 at about 1200 t·yr-1. This gradually decreased to about 200 t·yr-1 by 2006-2010. 3) The temporal change in NPLP coincided with the nitrogen concentration of the rivers in the study area. In addition, nitrogen concentrations in two test wells were 1.0 mg·l-1 at a depth of 150 m and only 0.25 mg·l-1 at 50 m, suggesting gradual percolation of the nitrogen polluted water deeper in the aquifer.
In general, people have a great concern about groundwater quality if they use the groundwater for municipal, industrial and agricultural purposes. Thus, the availability of unpolluted water is a very important factor to their daily life.
In the last few decades, the use of chemical fertilizer applied to farmland has markedly increased. Nitrogen from fertilizer not taken-up by crops percolates through the ground eventually reaching, and thus polluting, groundwater. In addition, the use of wastewater treatment works (WWTWs) has also markedly grown in recent decades in Japan, and nitrogen pollution of groundwater has occurred because wastewater treatment was usually limited to secondary stage processes which do not remove nitrogen. Therefore, nitrogen, not removed by WWTWs flows into the environment eventually reaching and polluting groundwater.
Besides nitrogen pollution originating from agricultural land and WWTWs, there are many sources of nitrogen pollutants such as deposition from precipitation, intake water (where water is taken off to feed irrigation canals), and recharge water from rivers. There are sinks for the removal of nitrogen pollutant components, such as surface and subsurface drainage away from the aquifer, denitrification in paddy fields, and effluent transport in the river. Therefore, to discuss nitrogen pollution of groundwater, all components relating to this pollution must be considered from an integrated (comprehensive) perspective.
Retention time of groundwater is usually very long because of the large capacity of the aquifer and the relatively small rate of recharge. Such characteristics demand that groundwater pollution issues be considered over long periods.
Many hundreds of studies of nitrogen contamination of groundwater have been carried out worldwide. We have listed only a few examples (because of the limitation of space): the nitrogen budget for a late successional hillslope tabonuco forest [
In addition, investigations of water and nitrogen cycles have also been conducted by many researchers. For example, Burt et al. [
Groundwater pollution issues have also recently been very intensively studied. Howden [
Many studies on the relationships between groundwater pollution and individual sources of nitrogen pollution are mentioned above. However, integrated approaches including all pollutant components are limited in number. Some long-term risk assessments of groundwater pollution were found [22,23]; these are necessitated by the long residence times associated with water movement in aquifers. Nitrogen balance investigations are also limited in number, despite their importance for the design of effective measures to decrease pollution.
The results of the aforementioned studies could not be applied directly here because water and nitrogen balance issues are strongly dependent upon the characteristics of the site being investigated. The overall goal of this study was to apply a water and nitrogen balance analysis to the Tedori River Alluvial Fan Area, Ishikawa Prefecture, Japan, to clarify both the NPLP over an extended period, and the importance of its source components, to reduce pollutant loads.
The Tedori River alluvial fan area is an important industrial and agricultural region within Ishikawa prefecture that contains many major industries and commercial facilities, as well as an abundance of developed agricultural land, The facilities utilize a lot of groundwater in the area from olden days while the agricultural land, mainly paddy, utilize much surface water from the Tedori River. Therefore, water in the basin is very important for people live in this region.
The study area comprises a semi-enclosed groundwater basin (
area is bound by the Sai River on the north-eastern side, the Kakehashi River on the south-western side, mountainous areas on the south-eastern side and the Japan Sea on the north-western side. Both sides of the Tedori River contain developed farmland consisting mostly of rice paddies and both sides of the Sai River and National Route 8, the main highway in the study region, are de veloped residential areas. The total study area contains 17,682 ha of flat land, of which paddies occupy 42.6%, upland fields 2.8%, residential areas 34.2%, road areas 15.0%, and rivers (including canals) 5.3% [
The elevation of the top of the fan is about 80 m above sea level and the distance from the top to the end of thefan is 11 km. The mean gradient is 1/140, which is relatively steep. The underground portions of the study area show an excellent aquifer horizon consisting of sand and gravel, with partially contained a clay lens. The region has, therefore, experienced extensive groundwater exploitation.
Additionally, to confirm the reliability of the nitrogen balance analysis, the nitrogen content of two small rivers in the study area, the Fushimi and Kurahashi Rivers, was analyzed by collecting water quality data after 1993 from the prefectural office. The nitrogen content of the Fushimi River was investigated every month. The river has a 6200-ha catchment area, a predominantly residential area located in the northern part of the study area. The nitrogen content of the Kurabe River was investigated every two months. The river has a catchment area of 1700 ha and is representative of land use in the central part of the study area.
Prior to nitrogen balance analysis, water balance analysis was conducted with the aim of clarifying the interrelationships among the various hydrologic components and estimating values for non-observed components for nitrogen balance analysis.
A two-horizon water balance model was applied in the study area to simulate the exchange of water between the channel/soil and aquifer horizon fractions and the Tedori River (
In the atmosphere:
In the channel-soil horizon:
In the aquifer horizon:
where Pr is precipitation, Prd is direct runoff, Prg is infiltration, Qin is the intake water to the study area from head works (irrigation canal off-take points), ET is evapotranspiration, Per is the percolation from the soil surface, Qout is the outflow discharge from the drainage canal, Pup is pumped water from the aquifer, Tin is recharge water from the Tedori River, Tout is effluent water into the Tedori River, and Gout is discharge from the aquifer to the lower reach of the drainage canal or the Sea of Japan directly.
The water balance model and its components are reported in detail elsewhere [
2.2.2.1. Precipitation and Evapotranspiration (ET)
Precipitation rain data, measured by tipping type automatic gauge, used in this study covered 35 years (1976- 2010) and were collected at the Kanazawa Local Meteorological Observation Station (hereafter Kanazawa), which is located about 10 km from the center of the study area. The mean annual precipitation for the region over the 35 year study duration was 2425 mm, 877 mm of which fell during the irrigation period (13 April to 10 September) and the remaining 1548 mm fell outside of the irrigation period (non-irrigation period). The precipitation in the non-irrigation period includes snowfall. Daily evapotranspiration (ET) over the same period was estimated by the Penman equation using data for Kanazawa. However, the ET of roads (15.0%) and roofs of houses (about 11.5%) was assumed to be zero; thus, only 73.5% of the total ET calculated by the penman equation was applied to the entire study area.
2.2.2.2. Water Use and Percolation Water use (water requirement), investigated by paddy lot, for rice paddies is a very important factor in the analysis because this differs markedly before and after farm land consolidation. Consolidation resulted in soil compaction caused by heavy machinery, with consequent reduction in soil permeability.
Farm land consolidation took place over thirty years between 1965 and 1995 (
The mean percolation in paddies during the irrigation period was estimated as the water use minus Penman evapotranspiration; the evapotranspiration coefficient was assumed to be 0.85. The percolation rate used for the entire study area (17,682 ha) was based on that for the cropped rice area.
Percolation water from irrigation is assumed to travel directly to the aquifer, because the study area has a deep groundwater level and a gravel and sand horizon beneath the topsoil.
2.2.2.3 Discharge in the Tedori River The amount of inflow to, or discharge from, the Tedori River is a very important water balance component. Therefore, we conducted a survey of discharge in several sections of the Tedori River during the irrigation and non-irrigation periods and calculated the difference in discharge between the two sections. Discharge clearly declined in the river stretch between 2.2 and 16.4 km
from the river mouth, indicating that water from the river was recharging the aquifer underlying the area. The discharge between 1.1 and 2.2 km from the river mouth increased slightly, reflecting the inflow of effluent from the area into the river.
2.2.2.4. Pumping Groundwater from Wells Ishikawa Prefecture has reported on the amount of groundwater pumped, investigated by pump capacity multing operation time, from wells in the area [
2.2.2.5. Intake Water for Paddy Irrigation To assess paddy intake water quantities, we used observed data surveyed at the top of the Shichika irrigation canal from 1976 to 2010. No observations were available for the Miyatake District, thus paddy intake water was estimated based on water right discharges between two the Districts. Five yearly averages observed intake water quantities are shown in
2.2.2.6. Estimation of Direct Runoff and Infiltration from the Ground Surface To estimate percolation and direct runoff, runoff analysis was conducted, using a ratio of precipitation volume to direct runoff for various land uses (
Land use areas were determined from statistical data for Ishikawa Prefecture. Direct runoff from rotational use of upland paddy areas was estimated using of those runoff ratio. Direct runoff from the paddy areas during the irrigation period was estimated using the runoff ratio for the non-irrigation period. The roof area percentage for residential areas (33.7%) was determined using GIS information in which the runoff ratio was assumed to be 100% [
Based on the above procedure, the amount of direct runoff and percolation from the ground surface for various land uses was estimated using daily precipitation and ET over 35 years divided into five-yearly intervals (1976- 2010).
Based on the water balance analysis, the nitrogen pollution load potential (NPLP) was estimated principally from the product of water flow volume and nitrogen concentrations (hereafter total nitrogen, TN). The nitrogen balance model was constructed from the water balance components adding nitrogen loads from farmland, effluxent from WWTWs and denitrification in paddies (
Nitrogen deposited by precipitation was monitored at weekly intervals over an extended period by the Ishikawa Environment and Health Research Institute [
The observed volume of intake water (the water extracted from the river into irrigation canals, investigated H~Q curve) is shown in
The nitrogen concentration of percolation water, equals to surface water, was measured at 53 surface locations four times during the irrigation period and three times during the non-irrigation period [
The nitrogen concentration of recharge water from the river was calculated according to the observed data; 0.36 mg∙l−1 during irrigation and 0.48 mg∙l−1 during the nonirrigation period in the upand mid-stream portions of the river. The effluent (river outflow) water concentration was determined to be 1.18 mg∙l−1 and 1.11 mg∙l−1 during the irrigation and non-irrigation periods, respectively [
The nitrogen concentration of the pumped (aquifer) water was assumed to be the average nitrogen concentration of the groundwater observed from 53 points across in the entire study area [
The NPLP in surface drainage water was based on measured nitrogen concentrations and estimated surface drainage flows. The concentration was measured near the Sea of Japan, four times at 5 - 6 sites during the irrigation period (n = 22, σ = 0.19), and four times at 6 - 8 sites during the non-irrigation period (N = 28, σ = 0.17) [
Similarly, the NPLP in groundwater discharge was estimated from measured nitrogen concentrations and estimated groundwater discharge. The nitrogen concentration was investigated at 63 locations during the irrigation period (σ = 0.60) and at 55 locations (σ = 0.70) during the non-irrigation periods [
In this section, we explain how the water and nitrogen balance analysis of each individual component of the entire system was performed.
Water balance analysis was conducted as outlined above [
In the table, Qout was estimated by the imbalance of channel/soil horizon in Equation (2) and Gout in Equation (3) was also estimated by the imbalance of aquifer horizon.
In the atmosphere zone, Pr, Prd and Prg were relatively constant. In the surface water zone, Qin was relatively constant except in period 5. Per showed a remarkable decrease because of the decrease in water demand caused by land development, in response to which Qout also gradually increased. In the groundwater zone, Tin and Tout were assumed constant throughout study period due to the lack of observational data and Pup was shown to be a relatively constant. Gout showed a marked decrease consistent with the decrease in Per. These results have a significant influence on the following nitrogen balance analysis.
NPLP from farmland was assessed in previous studies [37,38]; NPLP was estimated as the difference in nitrogen between fertilizer application rate and the yields of various crops. The fertilizer input was determined from standard application rates for various crops, as recommended by the local government office in the study area. The cultivated crops were rice, soybean, barley, vegeta-
bles, and orchard fruits. The nitrogen content in the unit yields for crops was determined from standard tables of food composition. The yield of crops was determined from statistical data collected by the prefectural govern ment for crops every five years. NPLP shows remarkable temporal changes (
The peak of NPLP occurred in 1990 at around 1000 t∙yr−1, decreasing gradually thereafter to about 400 t∙yr−1 during the period 2005-2010. The largest NPLP source was rice cultivation, because of the large total cropping area. The second was horticultural crops; vegetables and orchards, due to heavy fertilizer per unit area application rates.
Nitrogen in effluent from WWTWs is a very important factor for NPLP, and existing estimates were available for the study [
NPLP from WWTW was estimated from population equivalent nitrogen loads for various WWTWs; These were 4.8 g∙person−1∙d−1 for activated sludge, 1.6 g∙person−1∙d−1 for the oxidation ditch method (a kind of public sewage treatment), 6.1 g∙person−1∙d−1 for the rural sewage method, 6.5 g∙person−1∙d−1 for combined septic tanks, and 7.6 g∙person−1∙d−1 for night soil septic tanks [
The peak WWTW NPLP was about 800 t∙yr−1 in 1994 (
The nitrogen balance analysis results for the irrigation period and non-irrigation periods are presented in Tables 2 and 3 respectively; the NPLP was obtained from the product of the volume of each component and its nitrogen concentration (acronyms refer to
Tables 4 shows the input and
The difference between the input and output is the pollution amount, i.e. the mass of pollutant stored in each horizon. In addition, denitrification (De-N) from planted paddy was estimated as 17.1 kg∙ha−1 from the test paddy [
The total input of nitrogen peaked 2155 t∙yr−1 at period
(3), after which rapidly decreased 1035 t∙yr−1 at period (7). The dominant input components were Agr-N and SwN mentioned detail in following section. The other components keep relatively constant (
The total output of nitrogen peaked 1141 t∙yr−1 at period (1) after which rapidly decreased. The dominant decreased component is Gout-N, while Qout-N is gradually increased. The De-N gradually decreased because of paddy area decreasing. The other components maintained relatively constant (
The difference of input and output, index of pollution of the study area peaked 1129 t∙yr−1 at period (3), after which markedly decreases 174 t∙yr−1 at period (7). The main reason of the phenomena was remarkable decreasing of Sw-N and Agr-N in the input component.
In this section, we discuss nitrogen balance analysis in the aquifer horizon for all components, including the NPLP of farmland and WWTWs presented in the former section, and whether these indicate pollution of the aquifer.
One objective of our research was to determine the extent of groundwater pollution.
Unit: Amount (mm∙day−1), NPLP (t∙yr−1).
Unit: Amount (mm∙day−1), NPLP (t∙yr−1).
1981-1990 at 1800 t∙yr−1. After that it rapidly decreased to about 600 t∙yr−1 in the period 2006-2010. The largest component was NPLP from farmland at 1000 t∙yr−1 in the period 1976-1990, after which NPLP decreased markedly to about 400 t∙yr−1. The second largest component was effluent from WWTWs (600 t∙yr−1 in the period 1986- 1990), which decreased significantly to about 100 t∙yr−1 in the period 2006-2010. The other items contributed around 100 - 200 t∙yr−1. Therefore, to further reduce Npollution in the study area, NPLP from farmland and WWTW effluents would have the highest priority for reduction, despite their significant decreases in recent times.
The changes in output components are shown in
The difference between input and output, the nitrogen pollution load of the aquifer horizon, is shown in
The remaining issue is the accuracy of the groundwater discharge (Gout), which was estimated during the water balance analysis. Further investigation into the accuracy of Gout is necessary in future research.