Nitrate Contamination in Brazilian Urban Aquifers: A Tenacious Problem

: This study follows the geochemistry of nitrogen in a Cretaceous and unconﬁned sedimentary aquifer in the city of Ur â nia (Brazil) over 20 years. Although the sewer network was built in the 1970s, the nitrate contamination problem ( > 45 mg / L-NO 3 − ) persists to this day. The oldest urbanization areas located in the north of the city initially used cesspits for wastewater and currently present the highest nitrate concentrations ( > 120 mg / L-NO 3 − ), with the plume reaching the deeper aquifer portions (up to 100 m). The contamination is not as dramatic in the south part of the city, where urbanization including installation of the sewage network with PVC pipes that are more resistant to leak than the old ceramic networks occurred after 1985. Based on the constructive well proﬁles, three hydrogeochemical zones were established: shallow ( < 20 m deep), with average nitrate of 63 mg / L-NO 3 − ; intermediate (20–60 m), with 30 mg / L-NO 3 − ; and deep ( > 60 m), with 17 mg / L-NO 3 − . The current total nitrate mass in the aquifer exceeds 731 kg-NO 3 − . Numerical ﬂow (Modﬂow) and transport (MT3D) model scenarios support the hypothesis that the nitrate contamination is caused by substantial inﬁltration of nitrogen through the cesspits until the 1970s, but after the 1990s, leaks from the sewer network should be responsible for the maintenance of the recently observed high concentrations of nitrate.


Introduction
Nitrate is one of the most ubiquitous pollutants in groundwater in the world. At concentrations higher than 45 mg/L NO 3 − (10 mg/L N-NO 3 − ), it causes methemoglobinemia and cancer [1,2].
The origins of this contamination have been widely studied and are associated with the use of inorganic fertilizers in agriculture [3][4][5] and in intensive animal husbandry [6], in addition to inefficient sanitation systems (septic tanks or cesspits) [7,8] and leakage from public sewer mains [9][10][11].
The presence of nitrate in groundwater in cities is related to the management of its effluents; however, its distribution in the aquifer is also associated with urban infrastructure, which changes recharge patterns. Cities alter the circulation of water completely; thus, it is essential to relate sewage leaks to the urban occupation and the various sources of water that enter the aquifer. In Liverpool, England, losses in the water supply network resulted in leaks of 36.5 Mm 3 /y, with a refill equivalent to 180 mm/y [12]. In Hat Yai, Thailand, there was an increase of 97 mm/y for an urban refill of 240 mm/y (40.4%) [13]. In Toluca, Mexico, 27% of public water is lost, which is equivalent to 10.6 Mm 3 /y [14]. Regarding leaks from sewage networks, there are far fewer studies. In Nottingham and Doncaster, in the UK, the values were 0.33 Mm 3 /y and 0.13 Mm 3 /y, respectively [15,16]). In central London, the UK, it is estimated that there are losses of 0.9 Mm 3 /y.
In developing countries, where there is less maintenance of sewage networks, leakage is expected to be higher, reaching 5-15% of the total effluent circulating in the mains [17,18]. The contamination observed in Brazilian cities such as Natal (Rio Grande do Norte), Belém (Pará) and Sete Lagoas (Minas Gerais) [17][18][19][20][21] show that these values are realistic. The old sewage network, lacking maintenance and built with ceramic pipes, was quite common until the end of the 1990s; it is one of the causes of contamination of cities in the State of São Paulo [18], as is the case in the urban areas of São José do Rio Preto [22], Bauru [23,24], Marília [24], Presidente Prudente [25] and Urânia [11].
Nitrate contamination in an unconfined sedimentary aquifer has been studied in Urânia (Figure 1) since the 1990s. The city has a population of around 9000 people, a typical residential occupation, a common pattern to small urban centers in the interior of São Paulo. Dug wells and shallow tubular wells (up to 60 m) have shown high nitrate concentrations in plumes that extend over almost all of its urban area, although the sewage network has been installed since the 1970s and today reaches more than 98% of the population.
This article aims to establish the extent and evolution of the nitrate problem that has affected the Adamantina Aquifer under the city of Urânia since the founding in 1950. The work also quantifies the nitrate mass in the aquifer, using temporal and spatial analysis of land use and occupation. To accomplish this, potential sources of contamination were mapped during visits to almost all residences in the city. This study's originality is in monitoring the nitrate contamination evolution and the geochemical and isotopic groundwater for more than 20 years, which allowed us to evaluate the contaminating plume's future development using a numerical model of flow and transport.

Site Description
The municipality of Urânia is in the northwest region of the State of São Paulo, with an area of 210 km 2 and a population of 9135 inhabitants [26]. The urban center occupies 2.14 km 2 , where 80% of the population lives ( Figure 1). Since its founding in 1953 until the present, its main economic activity is agriculture. Between the 1950s and 1990s, cotton grew, followed by coffee and, to a lesser extent, corn, rice, and peanuts. The cultivation of oranges and grapes began in the 1980s and 1990s, respectively, when the use of inorganic nitrogen fertilizers was introduced systematically.
Urânia has a humid tropical climate [27], with two well-defined seasons: hot summers, with monthly averages above 25 • C (September-April); and cold winters, with monthly averages below 22 • C (May-August). Between 1988 and 2018, the Urânia rain station (B7-006) recorded average annual rainfall of 1317 mm. The months December to March are predominantly humid, with an average rainfall of 206 mm/month, while the months June to August are relatively drier, with an average rainfall of 21 mm/month. Geomorphologically, the region is in the Planalto Centro-Ocidental Paulista, characterized by slightly rolling, wide and low hills, with convex and flattened tops and average topographic levels of 458 masl. The drainage is not very dense and has, in general, a parallel pattern, mostly organized by consequent rivers. The main streams are the Comprido and the Matadouro (Figure 1), both belonging to the Turvo river basin. Geomorphologically, the region is in the Planalto Centro-Ocidental Paulista, characterized by slightly rolling, wide and low hills, with convex and flattened tops and average topographic levels of 458 masl. The drainage is not very dense and has, in general, a parallel pattern, mostly organized by consequent rivers. The main streams are the Comprido and the Matadouro (Figure 1), both belonging to the Turvo river basin.
Geologically, the region is inserted in the volcanic-sedimentary Basin of Paraná, filled with continental and marine sediments and basaltic lava from the Serra Geral Formation. In Urânia, Cretaceous sedimentary rocks of the Adamantina Formation (Bauru Group) outcrop and are Geologically, the region is inserted in the volcanic-sedimentary Basin of Paraná, filled with continental and marine sediments and basaltic lava from the Serra Geral Formation. In Urânia, Cretaceous sedimentary rocks of the Adamantina Formation (Bauru Group) outcrop and are deposited on the Serra Geral Formation. The Adamantina Formation, the unit with the most extensive surface distribution in the State of São Paulo, is composed of fine silty to very clayey sandstones, with a reddish-brown color. The grains are sub-rounded and well-selected, distributed in layers with a thicknesses of from 66 to 165 m. The Serra Geral Formation is a set of spills of acidic and basic basaltic lavas characterized by aphanitic dark gray to black rocks, and in Urânia, they are melanocratic and thin phaneritic, and have a thicknesses greater than 1000 m [28,29].
Two hydro-stratigraphic units occur in Urânia, the Adamantina Aquifer, which is the focus of the study. It is of the unconfined to semi-confined type, is slightly heterogeneous and isotropic, and is limited in its lower portion by the basalts. The Serra Geral Aquifer when outcropped is fractured, unconfined, anisotropic and heterogeneous. The average hydrodynamic parameters for the Adamantina Aquifer in regional values for hydraulic conductivity (K) are between 10 −6 and 10 −5 m/s. Its recharge occurs naturally through the rain infiltration [30]. The average annual rainfall varied from 1000 to 1500 mm, with potential monthly evapotranspiration between 50

Materials and Methods
The same group of researchers has studied Urânia since the 1990s, allowing monitoring of the territorial evolution and its urban infrastructure, as well as the geochemical changes of groundwater. Many of the benefits from this long period, including evaluations of hydrogeochemical data, have not been published.

Hydrogeological Conceptual Model of the Adamantina Aquifer
The conceptual model involved the survey and interpretation of data available over the years of study in Urânia. The altimetry data were obtained from the topographic sheets of Santa Albertina and Jales, scale 1:50,000 [31,32]. Geological and hydrogeological information was compiled from several authors [7,11,[33][34][35][36][37], as well as lithoconstructive profiles of tubular wells [38,39] and geophysical vertical electrical soundings carried out by Pálmer-Rivera (2002) [37].
The tubular and dug wells in the city were divided into three groups according to the drilling depth and the filter position. Dug wells generally draw shallower water, with surface depths of up to 20 m (shallow zone, SZ). Among the tubular wells, there are shallow wells, with depths up to 60 m (intermediate zone, IZ), and deep wells, with depths from 60 to 120 m (deep zone, DZ). Cagnon (2003) acquired hydraulic parameters of the Adamantina Aquifer by pumping test [7]. Measurements of depth, thickness, and hydraulic level of surface water bodies were obtained by Pálmer-Rivera (2002) [37] and Cagnon (2003) [7]. For the estimation of recharge and evapotranspiration in rural areas, we considered the premises of Bertolo et al. (2006) [34], who monitored the zone using pressure tensiometers, while for the urban area, we used data from Maldaner (2010) [30], who made estimations using water balance and water level fluctuation techniques.

Evolution of Urban Occupation and Assessment of Potential Contamination
First, we performed registration and mapping of dug and tubular wells, as well as septic tanks and cesspits (classified as active or inactive), using visits almost a thousand homes, interviewing owners and consulting the Municipal Secretariat's Saúde de Urânia database and SABESP (The Basic Sanitation Company of the State of São Paulo). Agricultural activities present in urban properties were also identified, including small gardens, using manure and inorganic nitrogen fertilizers. To analyze the history of land use and occupation as well as the evolution of the installation of public water supply and sewage networks, aerial photographs at various dates were used (1962, scale 1:27,000; 1965, 1:60,000; 1971, 1:25,000; 1985, 1:25,000), as well as Landsat-5 satellite images (scene 222/074 of 22 June 1997, with 3, 4, and 5 band-composition), Google Earth (scenes from 2000 to 2019), field verification, and interviews with employees of the city hall and SABESP. The data were organized in a geographic information system (ArcGis 10.2). The coordinate system was Universal Transverse Mercator (UTM) projection, zone 23, South American Datum (SAD) 69, with units in meters.

Sampling, Chemical, Physicochemical, and Stable Isotope Analysis
The water collections were made in three phases: (1) June 1998 to September 2001 (nine campaigns); (2) September 2007 to September 2009 (three campaigns); and (3) July 2018 and September 2019 (two campaigns). In total, 463 samples were collected in 49 dug wells, 40 tubular wells, and 4 multilevel monitoring wells over almost 20 years (Table 1). For a temporal assessment of chemical variations, 43 dug and tubular wells were selected, sampled during four campaigns, distributed over four seasons of a hydrogeological year. In the September 2001 campaign, water samples were collected for stable isotopes analysis (δ 18 O and δ 2 H) [40], dissolved organic carbon (DOC), and nitrate isotope (δ 15 N). Wells were selected based on nitrate concentrations and location concerning the history of use and land occupation. Water sampling was made directly at the wellhead using the well's pumping equipment. Care was taken to pump at least one volume of the well in the tubular wells, and pumping in the dug wells was carried out for 30 min. In the case of multilevel monitoring wells, the sampling method was low-flow using the pneumatic groundwater-sampling pump [41].

Estimating Mass of Nitrate Beneath Urânia
From the evaluation of all nitrate concentration data in the 2001 campaign iso-concentration maps, the relative thicknesses and areas of each nitrate concentration range (mg/L) were determined. The volume of each concentration range was associated with the respective nitrate concentration, resulting in the total mass value, in kg.
To estimate the effective mass of nitrate, the total porosity of 15% was adopted from the results of the pumping tests, interpreted by the method of Neuman (1972) [45]. These calculations were made for the shallow zone (SZ), with an average thickness of 12 m, and for the intermediate zone (IZ), with an average thickness of 48 m. The total mass calculation was made considering: Total mass (kg) = area (m 2 ) × average thickness (m) × total porosity

Numerical Model and Nitrate Transit Time Estimation
The MODFLOW code [46], developed by the US Geological Survey (version 4.3), was used to solve the groundwater flow equation using the finite difference method. For the simulation of nitrate transport in the saturated zone, we used the MT3DMS code (version 5.2) [47] contained in the Visual MODFLOW package [48]. The flow model was calibrated by comparing the calculated hydraulic head data with the actual data measured in the field. Monitoring levels estimated the transient condition in three tubular wells with periodic measurements. There was no calibration for the transport model, and only possible scenarios were simulated based on data from the aquifer nitrate concentration maps. However, simulations were checked with the general behavior of the plume observed in the field. The map described in this nitrate distribution work was used as an initial condition

Hydrogeological Conceptual Model of the Adamantina Aquifer
The Adamantina Aquifer occurs throughout the urban area of Urânia. It is an unconfined and locally semi-confined aquifer. Recharge is the result of rain infiltration, with average values of 260 mm/y [30], occurring from December to March, and due to leaks from the water and sewer mains [7]. The potentiometric surfaces follow the geomorphology of the area, where the groundwater dividers coincide with the surface watershed dividers. The groundwater flow lines converge to the Comprido and Matadouro streams, having an effluent character, with preferential direction from south to northeast. Figure 2 shows the potentiometric map of the Adamantina Aquifer, which was built from the hydraulic head measured in dug wells during the recharge period. Locally, the transmissivity (T) and hydraulic conductivity (K) are between 7 × 10 −3 and 9 × 10 −4 m 2 /s, and 7 × 10 −5 and 8 × 10 −6 m/s, respectively, characterizing permeable and homogeneous sediments, with a flow velocity between 1.2 × 10 −5 and 6.2 × 10 −7 m/s, resulting from hydraulic gradients between 1.0% and 2.5% and specific yield (Sy) of 0.1-0.2.
The aquifer was divided into the shallow (SZ), intermediate (IZ), and deep (DZ) zones according to the differences among the physicochemical characteristics of the waters in relation to the depth. SZ waters are captured by urban wells (up to 20 m deep) and range from chloride-nitrate-calcium to chloride-nitrate-sodium, with electric conductivity (EC) between 100 and 1500 µS/cm, pH between 4.4 and 6.5, and with nitrate reaching 422 mg/L-NO 3 − . The IZ is represented by tubular wells with depths between 20 and 60 m, with magnesian calcium bicarbonate waters, pH close to 7, and EC between 100 to 400 µS/cm, with decreasing concentrations of nitrate in depth. DZ is characterized by wells over 60 m deep, with sodium bicarbonate waters, pH between 8.5 and 11.0, and reduced calcium and magnesium concentrations with increased sodium (Figure 2   The Adamantina Aquifer's potentiometric surface map and conceptual hydrogeology cross-section. There are six groups of groundwater in Urânia (Table 2). Although there is a great dispersion and overlap of the results, it is possible to determine that the waters of the SZ (urban wells) have slightly more negative average signatures during the rainy season (δ 18 O = −7.5% ) than in the dry period (from δ 18 O = 7.1% ), and the waters of the deepest wells, which extract water from IZ and DZ, equal these waters (δ 18 O = 7.1% , from −6.8 to −7.6% ). These results are consistent with the conceptual model, indicating an aquifer with rapid recharges during the summer and with more negative signature waters. That part of this recharge is the result of losses in the water and sewage network, which originates in the mixture of the Guarani Aquifer System (δ 18 O = −8% ), which is supplied from Jales (14 km) and with some waters of the deep zone of the Adamantina Aquifer (δ 18 O = −7.25% ) to a lesser extent. Sewage water has a slightly more positive isotopic signature (δ 18 O = −7.75% ) than public water (δ 18 O = −8.0% ), suggesting that the population still uses water from the Adamantina Aquifer extracted from private tubular and urban wells. Water from Jales (GAS) - Water from Jales (GAS) - Public water service (Urânia)

Evolution of Urban Occupation and Assessment of Potential Contamination
Until 1959, Urânia was a district of the city of Jales (about 10 km to the southeast), during which time it had more than 17,000 inhabitants (including the current cities of Urânia, Santa Salete, and Aspásia). With the separation between 1993 and 1997, the municipality was established with an area of 210 km 2 . The reduction in agricultural importance, especially the cultivation of coffee, caused Urânia to lose its population over the years (Table 3). Urbanization began in the north and expanded to the south of the city. In the 1960s, water was supplied by dug wells, and domestic effluent was discharged into cesspits. In 1977, there was already a public water and sewage network collected from 59% of the houses, but that was concentrated in the north and northwest areas of the city. In the 1980s, when 64% of the domestic effluent was collected, the implementation of the sewer mains began in the south and southwest areas, culminating in 1994 in the collection of 89% of the total effluent. In 1996, the final extension of the mains was built, which would reach 96% of the houses by 2001 and 99% today (Table 3, Figure 4).  Agricultural activity in the urban area was identified and mapped on 21 small properties, all of which used organic fertilizers (manure) and only eight of which also used inorganic fertilizers in 1990s. This profile has not changed much over the years (Figure 4). The cesspits, which are used strictly in the city and extensively in rural areas today, have a diameter of up to 1.5 m, and depths between 5 and 10 m, often reaching the water table [8]. Usually, a cesspit receives all the wastewater from a residence, both from the bathrooms and the kitchen.
At the end of the 1990s, there were 72 cesspits located in the area with a sewage network, 57 of which were active, and 15 had been recently deactivated (<3 years). From interviews with residents at that time, three regions with different densities of cesspit were identified, considering the form of urbanization and hydrogeological characteristic (Figure 4): area A (NW), with 23 cesspit/km 2 ; area B (NE/E), with 48 cesspit/km 2 ; and area C (S), with 17 cesspit/km 2 . Area A was the one that first received the sewer mains in the 1970s. In area B, complete coverage only occurred recently, and area C is that where urbanization was, to a large extent, concomitant with the installation of the sanitation system. The cesspits were gradually deactivated, and there are records of at least 20 cesspit/km 2 in area B (NE/E) today.
Agricultural activity in the urban area was identified and mapped on 21 small properties, all of which used organic fertilizers (manure) and only eight of which also used inorganic fertilizers in 1990s. This profile has not changed much over the years (Figure 4).

Chemical and Physicochemical Analysis
The SZ presents predominantly calcium-chlorinate-nitrate and mixed water. In IZ, the water is calcium bicarbonate and calcium bicarbonate (nitrate), with variable impacts of contamination by domestic wastewater ( Figure 5). In the DZ, the water is calcium and sodium bicarbonate. The concentrations of sodium, chlorides, and nitrate decrease, followed by sulfate and potassium, and those of calcium, bicarbonate, and pH increase from SZ to IZ, showing a clear stratification of these ions. The average concentration of chloride in SZ was 19.68 mg/L and decreased to 14.26 mg/L in IZ, and the nitrate of 140 mg/L-NO 3 − in SZ changed to 81 mg/L-NO 3 − in IZ, which may indicate denitrification. Although dilution by rain plays an important role, as noted by the rapid variation in the signature of 18 O, which becomes more negative after periods of rain (−0.35% ). The concentrations of nitrate and chloride varied in the SZ over the course of a hydrological year due to the recharge of the aquifer that affects SZ more than the other deeper zones ( Figure 5). After the period of rain, the water level increases by 2-3 m.
The intensity of the rains also influences changes in the water composition of the SZ wells. In the summer of 2000-2001, for example, rainfall accumulated in the two months preceding the groundwater sampling (January and February, 521 mm) changing the hydrochemical signature of SZ waters from sodium chloride to sodium bicarbonate, while the 2008-2009 summer rains (February and March, 241 mm) did not alter this signature ( Figure 5).
Calcium and sodium bicarbonate water, with low nitrate (<1 mg/L-NO 3 − ) and chloride (4 mg/L), indicate that there is a small impact of urban contamination on the DZ and that this is the natural geochemical signature. In DZ, there is an increase in pH and sodium, accompanied by a reduction in calcium and magnesium concentrations through cation exchange mechanisms [50]. Another characteristic is that in both SZ and IZ, an oxidizing hydrochemical environment predominates, providing the existence of nitrate abundantly and stably. Waters with nitrite, ammonium (except in three wells), ammonia, and organic nitrogen are virtually non-existent. Dissolved organic carbon (DOC) was detected between 0.89 mg/L and 0.29 mg/L. Eh values range from 290 mV to 721 mV, showing relative homogeneity throughout all sampling campaigns, both between SZ and IZ and between dry and wet seasons. These results corroborate those observed in other regional studies of the Adamantina Aquifer, although it is recognized that there are more specific variations associated with the hydraulic heterogeneity of the aquifer, and where denitrification may occur [51]. Additionally, δ 15 N results did not show variations in their signatures between contaminated waters from SZ (δ 15 N from +8.74 to +11.13% ) and IZ (δ 15 N from +7.75 to +11.93%), not indicating denitrification processes, which would enrich the δ 15 N of the nitrate.
Regarding the territorial distribution, the highest concentrations of nitrate (>80 mg/L), chloride (mean value 20 mg/L) and EC (315 µS/cm) were detected in area B (Figure 6), close to the Comprido stream, where the water table is shallow, and where there is the highest density of cesspit still active or recently abandoned. In this location, IZ has low concentrations of nitrate, confirming the notion that the Comprido stream is the discharge of deeper levels of the aquifer (DZ). In area A (NW), the oldest occupied region (before 1962), nitrate is widely distributed throughout the urban area, especially in IZ. In area C, south of the city, there is less contamination in both SZ and IZ ( Figure 6).

Mass of Nitrate Beneath Urânia
From the concentration distribution maps ( Figure 6) and considering aquifer geological characteristics, the remaining nitrate mass was calculated for the shallow and intermediate zones of the Adamantina Aquifer under the city of Urânia (Table 4). A total of 731 kg-NO3 − is distributed 59%

Mass of Nitrate Beneath Urânia
From the concentration distribution maps ( Figure 6) and considering aquifer geological characteristics, the remaining nitrate mass was calculated for the shallow and intermediate zones of the Adamantina Aquifer under the city of Urânia (Table 4). A total of 731 kg-NO 3 − is distributed 59% in the IZ, and 31% in SZ. The results indicate that, although SZ is more contaminated, there is larger NO 3 − mass content at the IZ due to its greater thickness.

Numerical Model and Nitrate Transit Time Estimation
The simulation of contamination scenarios was based on the distribution of nitrate in the city of Urânia in 2001, using the Modflow code, which discretized the area into 89 rows, 93 columns, and seven layers, totaling 8277 cells, with minimum horizontal sizes of 10 × 10 m. The calibration was performed by trial and error, and the model was calibrated with an RMS error of 8%. Some scenarios included pumping wells, which have been running continuously for up to 100 years. The transport model was MT3D, which considered only longitudinal and transverse hydraulic dispersion, because, to date, there is no evidence of extensive denitrification in the area.
The modeled area was divided into ten zones with different hydraulic conductivities in the range of 6 × 10 −8 m/s and 3 × 10 −5 m/s. Total porosity was 10-15%, and effective porosity was 7-10% (Adamantina Aquifer). It was assumed that the aquifer was horizontally isotropic, and the vertical hydraulic conductivity was ten times lower than the horizontal. The basalt of the Serra Geral Formation was considered a regional aquitard, with a hydraulic conductivity of 6 × 10 −10 m/s. The model had two different recharges: the first associated with the city of 230 mm/y, and the rural one, with 183 mm/y, with no variation in time.
The modeled initial concentration (Co) in all cases was established as Co = 100 mg/L-NO 3 − (simulating 100% of the load), and the model results in a fraction of that value (C/Co). Various scenarios were modeled (Table 5, Figure 7), and in some cases, the sensitivity of the model to the variation in longitudinal dispersivity (αL) and the pumping of existing wells in 2001 was also simulated (the year with the highest number of active wells in the city). Flow and transport modeling start in year zero. As the main source of nitrate was the leakage from the sewage network, it did not occur homogeneously or in a single moment. It must be considered that the entire sewage network in the city was built with ceramic pipes, subject to leaks [18]. After the 1990s, these were replaced by plastic tubes, which have few leaks, as observed in the south of the city (Area C). Thus, the numerical modeling considered here, despite having a relationship with the aquifer's reality, is an exercise in understanding the evolution of the contaminating plume and supporting decision making for territorial occupation and the protection of urban groundwater.
Simulations 1, 2, and 3 ( Figure 7) show a continuous source of contamination throughout the area, representing leakage from the sewer network and the effect of hydraulic dispersion, with estimated values of 107.9 m longitudinal dispersion (αL) and transverse dispersion of 10.79 m. The results in the three situations are regionally very similar. The pumping effect can be seen by comparing simulations 2 and 4, where the contaminant mass was extracted by pumping the wells that also caused a smaller horizontal advance of the contaminating plume, but with greater vertical movements of the nitrate (Table 5 and Figure 7). 100% for the cesspits registered for 100 years [7] No No SZ and IZ: 13% of the plume will reach these areas. DZ: 6% of the plume will reach the zone.
No increase in nitrate concentrations in the aquifer in a fraction greater than 13% of the initial contamination.  These simulations 9 and 10 show that the cause of the contamination observed in Urânia was due to the combination of a large infiltration of nitrogen through the cesspits, from the beginning of the urban occupation until the 1970s, and after the 1990s, by leakage of the sewage network. Simulations 5 and 6 consider that the source of contamination ceases, that is, there is no more leakage from the sewer main and without the presence of cesspits after 10 years of activity. It appears that in 40 years (or 30 after the source ceases), there is still 20% of the DZ concentration and that some contaminating mass would be present after 90 years, however, in low concentration. The pumping of wells is effective a great help in reducing the volume and concentrations of the plume horizontally and vertically.
Simulations 7 and 8 reproduce what it would be like if the active cesspits in 2000 continued to operate, indicating that, although there was some local impact, the plume concentrations would never reach more than 10% of the initial concentration, showing that they do not support a plume of the dimensions of the plumes observed today.
In simulations 9 (without pumping) and 10 (with pumping), the area impacted by the contaminating load was reduced by 50% (one cell with a source and another without), showing that the repair of the sewer network would not have to be complete (100% efficient) to reduce the problems arising from its leakage significantly.
These simulations 9 and 10 show that the cause of the contamination observed in Urânia was due to the combination of a large infiltration of nitrogen through the cesspits, from the beginning of the urban occupation until the 1970s, and after the 1990s, by leakage of the sewage network.

Conclusions
The groundwater of the Adamantina Aquifer in Urânia is contaminated by nitrate in its shallow (up to 20 m deep, SZ) and intermediate (20-60 m deep, IZ) zones, despite the sewer network that currently covers 98% of the city. The contamination mainly occurs in the north-central part of the city (Areas A and B), the region of the first settlements, and where the sewer mains dates from the 1970s. In the southern region of the city (Area C), the sewer network was installed concurrently with the development of the neighborhood after the 1990s, and nitrate concentrations are lower (<3 mg/L) and located in the shallow aquifer. It is estimated that there are more than 731 kg of nitrate (NO 3 − ) in the entire aquifer, with 300 kg in shallow zone and the rest in intermediate zone.
The origin of nitrate contamination in Urânia is associated with residential domestic wastewater, which, prior to the 1970s, was fully infiltrated into cesspits. Even though the sewage network was installed more than 40 years ago, the nitrate contamination plume is still persistent. Numerical flow and transport modeling have shown that the time required for a plume of a persistent contaminant to be completely eliminated from the aquifer is more than 100 years, and that it is possible to achieve water potability for nitrate in 30-40 years (<45 mg/L-NO 3 − ), assuming that the source is completely removed. However, a strong correlation observed between recharge rates and variations in the nitrate concentration in the shallow zone suggests that this is a very dynamic system, and that there are still active sources of contamination. Including the still operational cesspits as nitrate sources in the aquifer's numerical simulations, the results indicate that they are not sufficient to sustain the contamination observed today (2020). Therefore, this study suggests that leaks from the sewage network, resulting from low maintenance, are the current primary sources of nitrogen to the aquifer and can sustain the persistent nitrate concentrations of the Adamantina Aquifer.