JBCS

INTRODUCTION

Rivers and streams are the main pathways by which terrestrial and atmospheric carbon (particulate carbon-PC), nitrogen (N), phosphorus (P), and other elements are delivered to coastal marine areas.¹ The general properties of river systems tend to reflect a combination of geomorphological attributes modified by a wide range of direct and indirect climatic and human influences.

Human activities have greatly changed the flow of substances into coastal marine areas.² Multiple problems (including, eutrophication, harmful algal blooms, hypoxia, fish kills) are caused by the accumulation of N, P, and other elements in river water in route to the ocean.³ Various activities can potentially generate nutrients along the coast: wastewater production, including the release of untreated sewage, agriculture, livestock, soil erosion, and industrial processes. Quantifying nutrient loads provides important information about biogeochemical processes in rivers, estuaries, and the downstream coastal zone.

Tropical rivers and streams are valuable resources, but pollution from urban and agricultural areas poses a threat to their water quality. There has been little quantification of the nutrient loads from coastal drainage basins in northeastern Brazil.

Quantifying nutrient loads provides important information about biogeochemical processes in estuaries and the coastal zone. The objective of this work was to determine the N and P loads from 12 hydrographic basins (24,403 km²) that flow into the estuaries along the northeastern coast of Brazil, considering both natural and anthropogenic sources from 1990-2016. Additionally, this work's objective was to determine how hydrographic basins have changed with alterations (natural and anthropogenic) over 27 years.

METHODS

Study area

The state of Pernambuco occupies an area of 98,938 km², which represents approximately 6% of the area of northeastern Brazil and 1.2% of the total territory of the country. The state has 187 km of coastline and is home to 12 major coastal rivers that flow toward the Atlantic Ocean.⁴ These rivers, which were the focus of our study, are Goiana-GO, Botafogo-BF, Igarassu-IG, Timbo-TB, Paratibe-PB, Beberibe-BE, Capibaribe-CB, Jaboatão-JB, Pirapama-PP, Ipojuca-IP, Sirinhaem-SI, and Una-UN (Figure 1).

Figure 1. Map showing the locations of the 12 studied river basins in northeastern Brazil. The colored circles indicate monitoring stations in each river. The colors in the sidebar show the geographic location of the sampling stations according to the different rivers. The major cities are highlighted in the left panel

The 12 systems differ by almost two orders of magnitude in area, from the larger Capibaribe (7,757 km²) basin to the smaller Beberibe (80 km²) basin. They are mostly short systems with lower basins located within the "Tabuleiros Costeiros do Nordeste" formation, which is characterized by Tertiary and Quaternary sediments that form coastal plains constituted by sandy soils closer to the coast and yellow-red latosols and podsols inland. Small stretches of alluvial eutrophic soils occur along river valleys.⁵ The longer system (Ipojuca River) crosses semi-arid inner areas where it can be intermittent (depending on the amount of rainfall), but it is perennial along its more coastal reach. Rainfall varies from ~ 2,000 mm yr^-1 on the coast to < 500 mm yr^-1 in the interior valleys (semi-arid).⁶

The natural vegetation in most of the area has been converted into land for non-mechanized subsistence agriculture and pastureland. The major crops in the sandy soils are coconut and banana, whereas the latosols are used mostly for sugar cane and pasture. Population density is generally high (> 1,000 persons per km²) and is concentrated in cities near the sea (see in the left panel of Figure 1; Paulista, Olinda, Recife, Jaboatão). Industrial activities are located mainly in areas near the shore. Most urban waste is not treated, and agriculture and livestock production have recently intensified along the coastal region. Therefore, prior data already suggest that some of these areas are showing signs of incipient eutrophication.⁷ Table 1 shows the environmental characteristics of the 12 rivers studied.

Emission factors methods

The loads associated with the natural edaphic and atmospheric processes were estimated by taking into account the composition and distribution of soils, rates of release and/or retention of the substances of interest in the soil types, and soil loss equations. The main nutrients with high eutrophication potential (N and P) were obtained through coefficients reported in the literature. The contribution of natural N and P emissions through atmospheric deposition was obtained by the average atmospheric deposition,^12,13,14 weighted by the actual annual rainfall on the river basins (mm) and corrected by factors of retention of each element in the soil (N: 0.37; P: 0.30).^15,16

where, [N; P] atm.dep is the total load of N or P via atmospheric deposition in t yr^-1; E is the average deposition of [N; P] in mg m² yr^-1, for rainfall of 1,000 mm; F is the retention factor [0.37; 0.30]; R is the rainfall over the period for the hydrographic basin in mm yr^-1, and A is the area of the river basin in km². The contribution of natural N and P emissions from physical and chemical denudation of soils was obtained by the factors indicated at for tropical soils.^13,14

where, [N; P] s.r. is the total load of N or P via soil runoff in t yr^-1; A is the area of the hydrographic basin in km²; M is the soil loss (12 kg yr^-1), and E is the loss factor of the element [N; P] in t yr^-1 with the loss factors for each type of soil.^13,14

To estimate loads from anthropogenic sources, the relationships between productive activities and their respective load of elements released into the environment are sought, such relationships are obtained using data on the production and use of agricultural inputs, volumes, and concentrations available for effluents, and retention rates in soils in the basin.¹³ Nutrient loads to industrialized coastal watersheds are mostly from diffuse sources and therefore are difficult to directly measure. For this study, we instead used indirect approaches based on emission factors as previously described.^7,14 We used emission factors available in the literature for each separate activity or process. However, all emission factors were adapted to local conditions as needed. For example, correction factors were used to adjust for semi-arid conditions to better estimate actual water consumption rates by the population in each basin.¹⁴

Emission factors for anthropogenic activities are based on the latest population data available (1990 census; 2016 population estimates) and factors indicated by Bidone and Lacerda, and Bidone.^17,18 Urban runoff depends on the area of impermeable surfaces, the number and size of buildings, annual rainfall and the urban area.^6,19

Nutrients emitted from livestock are released into soils as animal excrement, and their nutrient contributions depend on the types of animals raised in a region and the concentration of nutrients in their feed. The average values and ranges of excretion for horses, cattle, pigs, poultry, and sheep, as well as the rates of retention of nitrogen and phosphorus in soil and bodies of water were estimated through coefficients.¹⁴

The industrial loads of N and P were calculated using the Load Organic Remaining method (t BDO d^-1)⁶ and the coefficients were found by SanDiego-McGlone et al.²⁰ N and P contributions from sugarcane cultivation were estimated using the area of cultivation, fertilization practices, and loss of N and P. Data on the area of cultivation,⁹ fertilization information,¹⁶ and loss of nutrients¹¹ were derived from previous studies and coefficients for tropical areas. The potential for environmental degradation-PED is defined as the deterioration in environmental quality from ambient concentrations of pollutants and other activities and processes such as improper land use and natural disasters.⁷ The PED in the region was determined by calculating the sum of loads from atmospheric deposition and soil erosion and leaching based on soil physics and chemistry (Level-Based Natural) (NBN). Three levels were proposed: high (> 20 times NBN), medium (2-20 times NBN), and low (< 2 times NBN).^7,14

Statistical analysis

The statistical analysis used in this work were as follows:

Descriptive statistics: The mean, standard deviation, minimum, maximum, and percentages were calculated. Nonparametric statistics: The Mann-Whitney test was used to verify significant differences between 2 sample sets. The statistics for the time series were calculated using the Mann-Kendall trend test. In the Mann-Kendall test, a trend is considered negative or positive, indicating a decrease or increase in the attributes of the historical series analyzed, if the Mann-Kendall (Kendall's tau) score is negative or positive, respectively. In addition, the trend indicated by this methodology is considered significant when the p-value has a value less than α = 0.05. All statistical analysis were performed using mStatGraph v1.0 software.²¹

RESULTS

Rainfall and population density

As indicated above, rainfall varies from ~ 2,000 mm yr^-1 on the coast to < 500 mm yr^-1 in the interior valleys (semi-arid),⁶ while, the population density in the river basins varies between 100-8000 inhabitants km^-2 (Table 1). Figure 2 indicates the time series increment rate (1990-2016) for all basins in this study, for rainfall, and population density. The results indicated an annual average increase of +1.1% in population density, while precipitation showed no significant trend (Kendall test; p > 0.05). The time series indicated a significant trend for population density (Mann-Kendall test; p = 0.000006; α = 0.05) and not significant for rainfall (Mann-Kendall test; p = 0.2; α = 0.05). Differences between historical rainfall (1976-2016) and the study period (1990-2016) were not observed (Mann-Whitney test; p > 0.05; α: 0.05).

Figure 2. Time series increment rate (1990-2016) for all basins in this study, for rainfall, and population density. The p symbol indicates the result of the Kendall trend test. The values used correspond to the average annual values of rainfall and population density for all the river basins in this study

Natural emissions

Soil erosional and leaching loss dramatically increases with agriculture and depends on the soil type and climate, while, atmospheric deposition is a function of the basin area, annual rainfall, and the concentration of N and P in bulk deposition (dry and wet).⁷ Losses of N and P (atmospheric deposition + soil runoff) in basins ranged from 8.8 to 726.0 t yr^-1 and 1.2 to 28.0 t yr^-1, respectively (Table 2). Across all basins, 72% of N was due to soil and 28% was due to atmospheric deposition. For P, the distribution was 52% and 48%, respectively. Inputs of N and P were higher in the CB (Capibaribe River) and UN (Una River) basins, which have the largest areas (Tables 1 and 2).

Anthropogenic loads

The total anthropogenic loads of N and P for each basin studied are presented in Table 2. These estimates show N inputs varying from 797 to 43,650 t yr^-1 and P inputs varying from 215 to 1,685 t yr^-1, with maximum values at the CB (Capibaribe) and GO (Goiana) basins. Wastewater was a major source of nutrients to coastal areas, accounting for 50% of N and 44% of P among all natural and anthropogenic sources combined. The second-highest source of contributions was agriculture (sugar cane cultivation; 15,352 and 4,315 t yr^-1 of N and P, respectively), accounting for ~ 34% and 41% of all natural and anthropogenic contributions in the region, respectively (Table 2).

Total natural loads were 2,069 and 93 t yr^-1 for N and P respectively, while total anthropogenic contributions were 45,719 and 10,568 t yr^-1 for N and P, respectively. These values showed that anthropogenic inputs were 21 and 113 times higher than natural inputs for N and P, respectively. Anthropogenic sources accounted for ~ 96% of all loads.

The main anthropogenic nutrient source was urban use (wastewater, urban runoff) with 0.9 and 0.2 t km^-2 yr^-1 of N and P, respectively; followed by agriculture (sugar cane) with 0.6 and 0.2 t km^-2 yr^-1 of N and P, respectively.

Long-term load trends - N

The time series of natural and anthropogenic N loads are shown in Figures 345. The results of the time series of atmospheric deposition + soil runoff (Figure 3) did not show a statistically significant trend for the period 1990-2016 in the watersheds of this study (Mann-Kendall test; p = 0.08; α = 0.05). The anthropogenic N contributions that showed significant positive trends according to the Mann-Kendall test (p < 0.05; τ = + or −; α = 0.05) were: wastewater (p = 0.0001; τ = +0.98), industry (p = 0.0001; τ = +0.94) and urban runoff (p = 0.001; τ = +0.43); while negative trends were observed in livestock (p = 0.004; τ = −0.39) and agriculture (p = 0.0001; τ = −0.68).

Figure 3. Temporal series (1990-2016) of atmospheric N deposition + soil runoff (black line) and rainfall (red line) in the hydrographic basins (a)-(l). Arrows indicate the watersheds with the highest and lowest contributions

Figure 4. Temporal series (1990-2016) of N wastewater (black line), industry (red line) and urban runoff (blue line) in the hydrographic basins (a)-(l). Arrows indicate the watersheds with the highest and lowest contributions

Figure 5. Temporal series (1990-2016) of N livestock (black line) and agriculture (red line) in the hydrographic basins (a)-(l). Arrows indicate the watersheds with the highest and lowest contributions

In Figure 3c and Figure 3l, the largest N contributions by atmospheric deposition and soil runoff were in the Una River (UN) and Capibaribe River (CB) watersheds, while the largest N loads of anthropic origin were: wastewater and agriculture. The high loads via wastewater originated mainly through the Capibaribe River (CB), Pirapama River (PP) and Beberibe River (BE) basins. The N load through wastewater in the Beberibe River basin (BE) reached 92% of the total loads of natural and anthropic origin (Figure 4 and Table 2).

N contributions from agricultural sources (Figure 5) showed high loads in the watersheds of the Goiana River (GO) and Sirinhaem River (SI), while the smallest contribution was from the Paratibe River (PA), as indicated in Figure 5d, Figure 5j, and Figure 5h, respectively. N emissions through livestock showed higher inputs in the Jaboatão (Figure 5g) and Pirapama (Figure 5i) rivers. According to the statistical analysis (indicated above), agriculture and livestock showed negative trends in the N time series (Figure 5).

Long-term emissions trends - P

The time series of natural and anthropogenic P loads are shown in Figures 678.

Figure 6. Temporal series (1990-2016) of atmospheric P deposition + soil runoff (black line) and rainfall (red line) in the hydrographic basins (a)-(l). Arrows indicate the watersheds with the highest and lowest contributions

Figure 7. Temporal series (1990-2016) of P wastewater (black line), industry (red line) and urban runoff (blue line) in the hydrographic basins (a)-(l). Arrows indicate the watersheds with the highest and lowest contributions

Figure 8. Temporal series (1990-2016) of P livestock (black line) and agriculture (red line) in the hydrographic basins (a)-(l). Arrows indicate the watersheds with the highest and lowest contributions

The results of the time series of atmospheric deposition + soil runoff (Figure 6) did not show a statistically significant trend for the period 1990-2016 in the watersheds of this study (Mann-Kendall test; p = 0.07; α = 0.05). The anthropogenic P contributions that showed significant positive trends according to the Mann-Kendall test (p < 0.05; τ = + or −; α = 0.05) were: wastewater (p = 0.0001; τ = +0.99), industry (p = 0.0001; τ = +0.95) and urban runoff (p = 0.0001; τ = +0.42); while negative trends were observed in livestock (p = 0.002; τ = −0.42) and agriculture (p = 0.0001; τ = −0.70).

Wastewater had the highest P input (44%), where the Capibaribe and Pirapama river basins (Figure 7c and Figure 7i) showed the highest inputs, while the lowest input was observed in the Igarassu River (Figure 7e). P input via urban runoff, in a smaller proportion than wastewater, showed higher contributions from the Capibaribe and Beberibe river basins (Figure 7c and Figure 7a), while P input from industries also represented a low percentage of the total (8.5%; Table 2), concentrating important contributions in the basins of the Sirinhaem and Goiana rivers (Figure 7j and Figure 7d).

Similar to P inputs via industries, livestock loads represented < 10% of the total (Table 2). The watersheds of the Jaboatão and Igarassu rivers showed the highest P inputs to the adjacent estuarine systems (Figure 8g and Figure 8e). P inputs via agriculture represented the second-largest input of anthropogenic contributions (41%; Table 2). The watersheds of the Goiana and Sirinhaem rivers recorded the largest contributions (Figure 8d and Figure 8j). The smallest contributions corresponded to the Paratibe River (Figure 8h) which has a lower ground cover for agricultural activity.

DISCUSSION

Natural emissions

Soil loss and atmospheric deposition are the two major natural processes contributing to N and P loads in the studied basins. Losses from agricultural land in tropical regions average 130 t km^-2 yr^-1 for flat land without mechanized agriculture,²² which characterizes the coastal plains in northeastern Brazil.^7,14 N contributions due to soil loss in temperate climates without mechanized agriculture range from 75 to 230 kg km^-2 yr^-1, with an average of 133 kg km^-2 yr^-1, whereas P loads range from 5 to 50 kg km^-2 yr^-1.³

However, recent experiments carried out in southern Bahia-Brazil measured annual rates ranging from 48 to 123 t km^-2. The soils in this region (latosols, podzsols) are similar to those observed in our study region (Table 1).

Losses of N and P in the studied basins were within the lower ranges reported for areas with mechanized agriculture (average of 95 and 9 kg km^-2 yr^-1 for N and P, respectively). These values are similar to those reported for the Contas and Gongogi river basins in the state of Bahia, Brazil.¹³ The loads vary depending on basin size, as was reported for northeastern Brazil by Lacerda.¹⁴ Capibaribe and Una basins had the greatest losses of N and P, respectively. These two basins are the largest in Pernambuco state and together represent 58% of the total area studied.

Along the Brazilian coast, total N and P atmospheric deposition ranges from 80 to 300 mg N m^-2 yr^-1 and 4 to 10 mg P m^-2 yr^-1 in pristine vs. heavy industrial areas, respectively, with an annual rainfall of about 1,000 mm.²² Deposition also depends on the degree of coastal urbanization and industrialization. Inputs to estuaries from the atmosphere estimated using these parameters were low (average of 30 and 2 mg m^-2 yr^-1 for N and P, respectively), and similar values were obtained in the basins of northeastern Brazil (35 and 5.6 mg m^-2 yr^-1 for N and P, respectively).¹⁴ Natural source loads studied here are small compared to other urban areas but are similar to values reported in other areas of northeast Brazil. Estimated loads from natural sources represented ~ 4% of total emissions. Additionally, statistical analysis of the time series of N and P did not show significant trends for the watersheds of this included in this study.

Anthropogenic emissions

The most important anthropogenic sources of N and P were studied here, and their relative contributions varied depending on the degree of urbanization, population, and extent of agricultural lands.

Anthropogenic sources showed similar temporal trends for N and P in the 1990-2016 series. Within the anthropic sources; wastewater and industries recorded the highest contributions of N and P to estuaries in the adjacent coastal region. Accordingly, we can stipulate that the main factor responsible for these positive trends is population growth. The increase in N and P loads were a direct consequence of the population annual growth rate in the study region (1.1%). Additionally, N and P loads via natural sources accounted for < 4% of total contributions. Thus, the anthropogenic loads were 21 and 112 times higher than natural inputs for N and P, respectively. Anthropogenic sources accounted for ~ 96% of all contributions.

The intense population growth, mainly in hydrographic basins with high urban coverage and a low rate of sewage treatment, has increased loads of N and P in the region. The positive trends of these loads (mainly wastewater) represented an annual average increase of 1.1% for N and 1.2% for P, for all river basins included in this study. These percentages represent average annual loads of 251 t N yr^-1 and 54.8 t P yr^-1, respectively. Agriculture is a major source of excess N and P for the coastal region in non-urbanized areas. Generally, clay soils lose 10 to 40% of nutrients applied as fertilizers, while in sandy soils the loss may reach up to 25 to 80%.²³

The need for N fertilization is particularly high for sugar cane (90 to 275 kg ha^-1).¹⁴ Nutrient losses from this crop are 26-32% for N and 6-20% for P.¹⁶ Small basins may receive more N from agriculture. For example, the relatively small basins (Botafogo, Igarassu, Paratibe, Beberibe and Timbo) with 751 km², receive more N from agriculture (1,686 t yr^-1) than the larger Capibaribe basin (7,557 km²), which received only 1,241 t yr^-1 of N, due to the dominance of sugar cane.

The estimated loads, normalized by the total area of the hydrographic basins, showed that N loads reached 1,600 kg N km^-2 yr^-1 and 417 kg P km^-2 yr^-1, respectively. These values represent almost double the N reported at the Rio Grande do Norte basins for emissions from agricultural activity;¹⁴ and 7 times that registered in the Ceará region (Figure 9).²⁴ Additionally, P loads via agriculture also showed high rates when compared to other hydrographic basins in the tropical region. According to our estimates, P loads represented 2.4 times (172 kg P km^-2 yr^-1) the estimates in the Contas River (Brazil),¹³ and 1.5 times (266 kg P km^-2 yr^-1) in the Rio Grande do Norte basins.¹⁴ Thus, the contributions of N and P normalized by the area of agricultural activity were higher than other basins in northeastern Brazil. On the other hand, livestock showed similar values when compared to other estimates from tropical regions. The values estimated in this study (368 kg N km^-2 yr^-1 and 207 kg P km^-2 yr^-1) are similar to those reported by other authors.^14,24,25

Figure 9. Comparison of N and P loads (kg km-2 yr-1) (logarithmic scale) from basins in 4 sites in northeastern Brazil; Ceará State,14 Rio Grande do Norte State,25 Contas River,13 and Pernambuco State

Nutrients loaded from livestock activity are released to the soil as animal excrement, and their concentration depends on nutrient concentrations in animal feed and on the types of animals raised in a region.⁷ Contributions of N and P due to livestock were independent of the basin area. Jaboatão and Igarassu basins are 17 and 52 times smaller than the Capibaribe basin respectively but received more N and P from livestock due to the predominance of poultry.

Inputs due to industrial activity represent 10% of total contributions, and basins with higher industrialization had the highest inputs of N and P, independent of the size of the basin. Significant positive trends of N and P were observed in the study region, including the 12 hydrographic basins. These basins will show increases of ~ 1% (at t yr^-1) for N and P, respectively.

Major parameters determining nutrient loss by runoff from urban areas are the area of impermeable surfaces, the number and size of buildings, and annual rainfall.

The results show a direct relationship between N and P loads and the extent of the urban area. The highest estimated N and P loads corresponded to basins with a high urbanization rate (Capibaribe, Beberibe, Paratibe rivers). Rainfall did not contribute much to variation between basins due to the little difference in rainfall totals.

The total urban runoff contribution corresponded to < 1%, similar to the values estimated in Contas River.¹³

Wastewater is one of the major sources of nutrients to coastal areas, particularly in urbanized estuaries.^7,14,24 When no treatment plants exist, nutrient loads from this source are directly proportional to the human population and the amount of water used per inhabitant.^7,25

The estimates show that N and P inputs due to wastewater are independent of basin size and proportional to the population. The inputs of N and P from wastewater in basins in the state of Pernambuco were higher than other basins in northeastern Brazil (Figure 9).

Basins with a high population, such as Capibaribe, Beberibe and Pirapama, had more emissions than basins with smaller populations. The population density was also highly correlated with N and P loads due to wastewater (r = 0.5; r = 0.7 and r = 0.4, respectively).

The estimated trends showed a positive annual increment of 1.2% for N and P, respectively.

There is a ~ 1% annual population growth rate for this region, which results in an increase in the volume of wastewater every year and an increased nutrient input.

The potential for environmental degradation-PED

The PED (estimated through the levels of natural base and the anthropogenic loads of wastewater and urban runoff) showed high values for N and P in the basins of the study region. Estimated PED levels for N showed 50% of catchments with high PED (n = 6) and 50% of catchments with a medium PED level (n = 6). Low levels of PED (< 2 NBN) were not recorded. Additionally, the P showed 92% of basins (n = 11) with a high level of PED and 8% (n = 1) with medium levels.

The results showed that in basins with high population density, PED levels were higher (Paratibe, Jaboatão, Pirapama, Beberibe) for N and P, respectively. When compared to other regions of northeastern Brazil, the watersheds in this study showed higher levels of PED.

Watersheds in the state of Ceará showed low levels of PED for N and P, respectively (N_PED = 1; P_PED = 0.5), while in basins of Rio Grande do Norte high and medium PED levels were estimated (N_PED = 2.2; P_PED = 2). In the Contas River basin, PED levels were low and high for N and P, respectively (N_PED = 0.3; P_PED = 10).

As indicated above, the estimated PED levels in Pernambuco were high (N_PED = 20; P_PED = 77), indicating a strong anthropogenic pressure in this region, mainly in river basins with high population density and low levels of sewage treatment.

The timeline of N and P loads showed that loads have steadily increased since 1990, mainly wastewater, as seen in Figure 10. N and P contributions through wastewater have increased by 37% in 27 years; while natural loads increased by 4% and 6%, respectively (Figure 10a and Figure 10b).

Figure 10. Sectoral comparison of N(a) and P(b) loads in the study region. (AD + SR): atmospheric deposition + soil runoff; (WW): wastewater; (IN): industry; (UR): urban runoff; (LI): livestock; (AG): agriculture. Values in t yr-1 (1990 and 2016) and % (1990-2016)

The reduction of N and P loads from agricultural activity (~ 50%) is indicative of changes in new agricultural technologies that reduce N and P contributions; but also, in the reduction of the land occupation area for this activity. However, N and P loads via agricultural activities still represent high values when compared to other river basins in the northeast region of Brazil. Additionally, population growth has generated greater urban runoff due to the growth of cities inserted in the watersheds of this study.

CONCLUSION

These results demonstrate that anthropogenic sources (wastewater, industry, livestock, urban runoff, and agriculture) produce sufficient nutrient levels to cause eutrophication and environmental degradation. The results indicated an average annual increase of +1.1% in population density, while precipitation showed no significant trend.

Anthropogenic sources accounted for ~ 96% of all contributions. The time series of N and P in the study region showed trends that point to a constant increment, implying a greater contribution of these loads to the adjacent coastal region year after year. N and P emissions from wastewater have increased by 37% in 27 years, while natural loads increased by 4% and 6%, respectively. High population density basins such as Beberibe, Jaboatão, Paratibe, for example, have 2 and 10 times the upper acceptable limit for urban N and P emissions, respectively. Estimated PED levels showed high levels of N in 50% of the basins studied, while P showed high levels of PED in 92% of the basins in this study. The comparative analysis of the study region with other coastal regions of northeastern Brazil showed higher loads of N and P from wastewater to the Atlantic Ocean in the period 1990- 2016. Nutrient levels in these basins are sufficient to generate serious environmental problems in the rivers and adjacent estuarine region.

ACKNOWLEDGMENTS

INCT AmbTropic, Brazilian National Institute of Science and Technology for Tropical Marine Environments, and CNPq/FAPESB provided financial support (Grants 565054/2010-4 and 8936/2011). M. A. and C. N. received support from the Brazilian Research Network on Global Climate Change (Rede CLIMA) (FINEP Grant 01.13.0353-00). C. N. was supported by the National Oil Agency- ANP and the PRH-ANP 38.1 program with grant (process: 044819). L. V. acknowledges the TRIATLAS project, which has received funding from the European Union's Horizon 2020 Research and Innovation Program under grant agreement nº 817578 and the International Joint Laboratory TAPIOCA (IRD-UFPE-UFRPE).

REFERENCES

1. Huang, T. H.; Chen, C. T. A.; Tseng, H. C.; Lou, J. Y.; Wang, S. L.; Yang, L.; Biogeosciences 2017, 122, 1239. [Crossref]

2. da Cunha, L. C.; Buitenhuis, E. T.; Le Quéré, C.; Giraud, X.; Ludwig, W.; Global Biogeochem. Cycles 2007, 21, 1. [Crossref]

3. Howarth, R.; Anderson, D.; Cloern, J.; Elfring, C.; Hopkinson, C.; Lapointe, B.; Issues in Ecology 2000, 7, 1. [Link] accessed in April 2023

4. Agência Estadual de Meio Ambiente - CPRH; Monitoramento das Bacias Hidrográficas de Pernambuco 2001-2005. [Link] accessed in May 2023

5. Silva, L. F.; Solos Tropicais: Aspectos Pedológicos, Ecológicos e de Manejo, 1ª ed.; Terra brasilis: São Paulo, 1996.

6. Agência Estadual de Meio Ambiente - CPRH; Report on the Monitoring of Water Quality in the Watersheds of the State of Pernambuco in 2016. [Link] accessed in April 2023

7. Noriega, C.; Araujo, M.; J. Coastal Res. 2009, 56, 871. [Link] accessed in April 2023

8. National Meteorological Institute (INMET); Meteorological Database for Teaching and Research - BDMEP. [Link] accessed in April 2023

9. Brazilian Institute of Geography and Statistics (IBGE), available at http://www.ibge.gov.br/cidades, accessed in April 2023.

10. Noriega, C.; Araujo, M.; Sci. Rep. 2014, 4, 1. [Crossref]

11. Santos, H. G.; Jacomine, P. K.; Anjos, L. H.; Oliveira, V. A.; Lumbreras, J. F.; Coelho, M. R.; Sistema Brasileiro de Classificação de Solos, 5ª ed.; Embrapa Solos: Brasília, 2018.

12. Silva-Filho, E. V.; Sella, S. M.; Spinola, E. C.; Santos, I. R.; Machado, W.; Lacerda, L. D.; Microchem. J. 2006, 82, 196. [Crossref]

13. de Paula, F. C. F.; Lacerda, L. D.; Marins, R.; Aguiar, J.; Ovalle, A.; Falcão Filho, C.; Quim. Nova 2010, 33, 70. [Link] accessed in April 2023

14. Lacerda, L. D.; Braz. J. Aquat. Sci. Technol. 2006, 10, 13. [Link] accessed in April 2023

15. Bouwman, A. F.; Lee, D. S.; Asman, W. A. H.; Dentener, F. J.; Hoek, K.W.; Van Der Olivier, J. G. J.; Global Biogeochem. Cycles 1997, 11, 561. [Crossref]

16. Malavoltas, E.; Dantas, J. In Melhoramento e Produção do Milho no Brasil; Paterniari, E., ed.; Fundação Cargill: São Paulo, 1980, p. 429.

17. Bidone, E. D.; Lacerda, L. D. In A Preliminary Approach of the Link Between Socio-Economic and Natural Indicators into a Driver- Pressure-Impact-Response Framework Case Study: Guanabara Bay Basin, Rio de Janeiro, Brazil, LOICZ Reports & Studies 21; Lacerda, L. D.; Kremer, H. H.; Kjerfve, B.; Solomons, W.; Marshall-Crossland, J. I.; Crossland, J. C., eds.; LOICZ: Textel, 2002, p. 142.

18. Bidone, E. D. In Análise Econômica-Ambiental Aplicada à Contaminação de Águas Fluviais de Pequenas Bacias Costeiras do Estado do Rio de Janeiro; Esteves, F. A.; Lacerda, L. D., eds.; Ed. UFRJ: Rio de Janeiro, 2000, p. 371-394.

19. Gianessi, L. P.; Peskin, H. M.; An Overview of the RFF Environmental Data Inventory: Methods, Sources and Preliminary Results, 1st ed.; Resources for the Future: Washington, 1984.

20. SanDiego-McGlone, M.; Smith, S. V.; Nicolas, V.; Mar. Pollut. Bull. 2000, 40, 325. [Crossref]

21. Varona, H. L.; Noriega, C.; Araujo, M.; mStatGraph version 1.0, Zenodo, Switzerland, 2021.

22. Greeland, D. J.; Lal, R.; Soil Conservation and Management in the Humid Tropics, 1st ed.; Wiley: Chichester, 1977.

23. Howarth, R. W.; Billen, D.; Swaney, A.; Townsed, N.; Janarski, K.; Lajtha, K.; Downing, J. A.; Elmgren, R.; Caraco, N.; Jordan, T.; Berendse, F.; Freney, J.; Kudeyorov, V.; Murdcoh, P.; Liang, Z. Z.; Biogeochemistry 1996, 35, 75. [Crossref]

24. Lacerda, L. D.; Molisani, M. M.; Sena, D. L.; Maia, L. P.; Environ. Monit. Assess. 2008, 141, 149. [Crossref]

25. Vasconcelos, V. H. F.: Emissões Naturais e Antrópicas de Nitrogênio e Fósforo para os Principais Açudes da Bacia Hidrográfica do Rio Seridó, RN; Dissertação de Mestrado, Universidade Federal do Rio Grande do Norte, Natal, Brasil, 2011. [Link] accessed in May 2023