Next Article in Journal
Comparative Evaluation of Crash Hotspot Identification Methods: Empirical Bayes vs. Potential for Safety Improvement Using Variants of Negative Binomial Models
Previous Article in Journal
How Do Existing Organizational Theories Help in Understanding the Responses of Food Companies for Reducing Food Waste?
Previous Article in Special Issue
Analysis of the Water Resource Carrying Capacity in Guyuan City
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact of the Citrus Industry on the Water Quality of the Filobobos River in Veracruz, Mexico

by
Luis Carlos Sandoval Herazo
1,
Humberto Raymundo González-Moreno
2,3,
Luis Guillermo Carreto-Hernandez
4,
Florentina Zurita
5,
Graciela Nani
1,
Sergio Zamora
6,
Mayerlin Sandoval-Herazo
1,* and
Georgina Martínez-Reséndiz
1,7,*
1
Wetlands and Environmental Sustainability Laboratory, Division of Graduate Studies and Research, Tecnológico Nacional de México/ITS de Misantla, Km. 1.8 Carretera a Loma del Cojolite, Misantla 93821, Veracruz, Mexico
2
Division of Graduate Studies and Research, Tecnológico Nacional de Mexico/ITS de Misantla, Km. 1.8 Carretera a la Loma del Cojolite, Misantla 93821, Veracruz, Mexico
3
Departamento de Ingeniería Civil, Tecnológico Nacional de Mexico/ITS de Misantla, Km. 1.8 Carretera a la Loma del Cojolite, Misantla 93821, Veracruz, Mexico
4
Tecnológico Nacional de México, Centro Nacional de Investigación y Desarrollo Tecnológico, Interior Internado Palmira s/n, Colonia Palmira, Cuernavaca 62490, Morelos, Mexico
5
Environmental Quality Research Center, Centro Universitario de la Cienega, University of Guadalajara, Ocotlan 47820, Jalisco, Mexico
6
Facultad de Ingeniería, Construcción y Hábitat, Universidad Veracruzana, Bv. Adolfo Ruíz Cortines 455, Costa Verde, Boca del Rio 94294, Veracruz, Mexico
7
Postdoctoral Position by CONAHCYT (Consejo Nacional de Humanidades Ciencias y Tecnologias) in Tecnológico Nacional de México Campus Misantla, Misantla 93821, Veracruz, Mexico
*
Authors to whom correspondence should be addressed.
Sustainability 2024, 16(4), 1536; https://doi.org/10.3390/su16041536
Submission received: 28 December 2023 / Revised: 30 January 2024 / Accepted: 8 February 2024 / Published: 11 February 2024
(This article belongs to the Special Issue Advances and Challenges in the Sustainable Water Management)

Abstract

:
Veracruz, Mexico, boasts abundant water resources. However, the region is facing challenges stemming from population growth and industrial expansion, leading to a notable increase in wastewater discharge into surface water bodies. This discharge has resulted in significant alterations to water quality. In this study, we analyzed the water quality in the Bobos River, Veracruz, with a focus on sustainability. A total of 12 samples were collected from various points along the reservoir. Parameters including temperature, pH, dissolved oxygen (DO), conductivity, total nitrogen (TN), nitrate, chemical oxygen demand (COD), total phosphorus (TP), phosphate, dissolved solids (DSs), and fecal coliforms (FCs) were measured and analyzed. The pH levels ranged from slightly acidic (6.40 ± 0.71) to slightly alkaline (8.65 ± 0.07), with the lowest dissolved oxygen concentration recorded at 0.15 ± 0.07 mg L−1. Conductivity varied between 0.26 and 3.81 mS. Total nitrogen concentrations ranged from 0.21 ± 3.8 to 0.491 ± 1.3 mg L−1. Dissolved solid concentrations were measured at 0.39 ± 0.05 mg L−1, while the final sampling point exhibited elevated values for temperature, nitrate, chemical oxygen demand, total phosphorus, and coliforms (30.35 ± 0.21 °C, 0.35, 389.00 ± 11.31, 13.20 ± 0.85 mg L−1, and 3.2 × 105 ± 3.707107 × 104, respectively). To address these concerning trends and safeguard public health and environmental integrity, the implementation of continuous monitoring and stringent control measures is imperative.

1. Introduction

Urbanization, population growth, and industrial and agricultural development can significantly impact river quality [1]. Research indicates a strong correlation between urbanization [2] and water quality deterioration. As cities expand, water pollutants, including chemicals, heavy metals, and nutrients [3], increase. In many developing countries, insufficient public policies fail to adequately protect natural resources like water [4,5]. This is particularly evident in rivers receiving untreated wastewater, which diminishes their dilution and self-purification capacities. The economic repercussions of river pollution affect vulnerable sectors, including fishing [6], livestock [7], and agriculture [8], and pose extinction risks to ecosystem organisms [9].
In Mexico, major population centers are often located along riverbanks, leading to over 70% of water bodies exhibiting some contamination and severe issues regarding availability and access to this resource [10]. Approximately 40% of generated wastewater, both industrial and domestic, is discharged into surface water bodies without prior treatment [11]. Veracruz State, Mexico, boasts significant water resources, with 17 hydrographic basins and 24 rivers. However, these rivers are impacted by wastewater from population centers and industries, making it unclear if they have any treatment processes. Mexican standards NOM-001-SEMARNAT-1996, NOM-002-SEMARNAT-1996, and NOM-003-SEMARNAT-1997 set limits for wastewater discharges.
All water, for public or industrial use, must undergo treatment before reintroduction into the cycle. Unfortunately, Veracruz has only a few registered treatment plants, with 74 urban wastewater treatment plants identified in 2002, treating only 30% of wastewater. The number of treatment plants in the 212 municipalities of Veracruz is insufficient, leading to inadequate treatment of urban wastewater that only focuses on certain parameters such as solids, carbon, and coliforms and can be considered as non-functional.
Access to clean drinking water is indispensable for fostering human health and promoting sustainable development. Regrettably, human activities pose a threat to water quality across diverse regions globally, and the Bobos River basin in southeastern Mexico stands as no exception to this concerning trend. Specifically, the discharge of waste originating from the citrus industry significantly compromises the river’s water quality. Hence, this study undertakes a comprehensive physicochemical and microbiological assessment of water quality within the Bobos River, with a primary focus on analyzing the detrimental impact caused by the citrus industry in the area. Furthermore, it endeavors to address environmental sustainability concerns surrounding industrial operations by advocating for the implementation of sustainable measures aimed at mitigating adverse effects and fostering the preservation and conservation of the Bobos River and other water bodies alike.

2. Materials and Methods

2.1. Description of the Bobos River Basin

The study area encompasses the Nautla River basin, traversing the municipalities of Plan de Arroyos, Tlapacoyan, Martínez de la Torre, San Rafael, and Nautla, located in Veracruz, Mexico, until its convergence with the Gulf of Mexico. The basin exhibits a maximum elevation of 3299 m and a minimum of 6 m, spanning approximately 121 km in its main channel. Its average slope measures 3.25%, with a retention time of 520 min. The region sustains a population of 366,500 inhabitants, residing in approximately 97,500 occupied dwellings.
It is noteworthy that roughly 92,500 homes are equipped with piped water infrastructure, while 3400 lack such access. Additionally, 90,500 households benefit from sewerage systems, while 6500 do not. Notably, there are settlements close to the Bobos River (Figure 1).

2.2. Sampling

2.2.1. Area of Study

The study area encompasses the Nautla River basin, situated in the states of Puebla and Veracruz, Mexico, spanning a total area of 3544 km2. This basin extends across several municipalities, including Puebla, Teziutlán, Martínez de la Torre, Misantla, and Nautla. The Nautla River, with a length of 121 km, eventually merges into the Gulf of Mexico. The relief of the basin is predominantly mountainous, featuring altitudes ranging from 50 to 3299 m above sea level. The region exhibits a warm and humid climate, characterized by annual precipitation levels varying from 1000 to 3500 mm [12].
On the other hand, the Bobos River basin, situated in Veracruz, Mexico, encompasses an area of 1417 km2 and extends over a length of 77 km. The Bobos River itself measures 58 km in length, with an average slope of 1.5%. The region experiences an average annual precipitation of 2229 mm, with 75% of rainfall concentrated during the June to October season. However, intensified agricultural and livestock activities, coupled with deforestation and urbanization, have heightened the basin’s susceptibility to soil erosion and sedimentation [13].

2.2.2. Water Samples

Six water samples were systematically collected at twelve distinct points along the Bobos River (Figure 1) in the summer season (June) of 2022. The sampling points were sequentially numbered in ascending order as the river flowed downstream, ranging from SW1 to SW12. High-density polyethylene containers, each with a capacity of 2 L, were employed to gather the water samples. Concentrated nitric acid was added to the containers to attain a pH level below 2, ensuring preservation during transportation and subsequent analysis.
Following collection, each container was meticulously labeled with the corresponding date, name, and location of the sampling site. Subsequently, the samples were securely stored in coolers, adhering to established transportation protocols and guidelines [14,15]. Sampling was conducted between 6:30 am and 9:00 am.
pH was measured at each site immediately after sampling using a portable pH meter (HANNA Instruments, Smithfield, RI, USA, HI98130).
Electrical conductivity was measured at each site immediately after sampling using a portable meter (HANNA Instruments, Smithfield, RI, USA, HI98130).
Temperature was measured at each site immediately after sampling using a portable meter (HANNA Instruments, Smithfield, RI, USA, HI98130).
Dissolved oxygen was measured at each site immediately after sampling using a portable meter (Milwaukee, Mulund West, Bombay, Maharashtra, India, model 600).
Parameters such as chemical oxygen demand (COD) [16], total nitrogen, NO3, total phosphorus (TP), PO4, and fecal coliforms were determined at the Wetlands and Environmental Sustainability Laboratory located at the Instituto Tecnológico Nacional de México Campus Misantla, Veracruz. Samples were preserved at 4 °C until processing; the techniques used were taken from the Mexican Standards which are based on the standard methods for water and wastewater analysis [15].
The research point was determined using a “sample survey method”, which is a sampling method carried out by dividing the research area into points that are expected to represent the study population.
The determination of sampling points was based on consideration of ease of access, costs, and time so that the points were considered to represent the water of the Bobos River.
The initial step in determining the location of river water sampling is to know the geography of the river and the activities around the watershed [17,18,19].

2.2.3. Quantification of NaOH (Sodium Hydroxide) (%)

The quantification of NaOH percentage was performed using the titration method, where drops of phenolphthalein were added, and titration was carried out with sulfuric acid (0.1 N). The NaOH percentage was determined using the following formula:
%NaOH = Volume of acid used × Acid normality×meq of NaOH (0.04) × 100/Weight of the analyzed sodium sample.

2.3. Statistical Analysis

To investigate the relationship between water quality factors, correlation analysis and principal component analysis (PCA) were performed using Statistical Package for the Social Sciences software (SPSS, ver. 25.0; SPSS Inc., Chicago, IL, USA).

3. Results and Discussion

3.1. Temperature and Dissolved Oxygen

Temperature and dissolved oxygen (DO) stand as pivotal parameters in the assessment of water quality within aquatic systems. In this study, we observed significant variability in water temperature across the Bobos River basin, influenced by both geographic location and human activity (Table 1 and Figure 2).
Dissolved oxygen (DO) emerges as a critical parameter defining water quality and the vitality of aquatic ecosystems [20]. Aquatic flora and fauna, including fish, rely on oxygen for survival. For instance, fish cannot endure prolonged periods in water with dissolved oxygen levels below 5 mg L−1 [21]. Reduced levels of dissolved oxygen in water signal contamination and serve as crucial determinants in water quality assessment, pollution management, and treatment procedures [22].
The increment in water temperature along the river can be attributed to the convergence of tributaries with elevated temperatures and the reduction in water velocity as the river broadens. In certain instances, notable increases in water temperature stem from the discharge of industrial wastewater [22]. Regular monitoring of temperature levels and other water quality variables is imperative to ensure the sustainability of the river ecosystem in the Bobos River basin, even though water temperature values fall within normal ranges.
Our study unveiled that the highest dissolved oxygen (DO) value in the river was documented at SW4, accompanied by an average temperature of 28.68 ± 0.32 °C. Notably, there exists a temperature differential of 10 °C between the lowest and highest points. Rajesh et al. [23] note that for every 1 °C elevation in river water temperature, there is an approximate 2.3% decline in DO saturation level concentrations. However, our study did not discern this relationship, potentially influenced by climate-change-induced warming of river water [24] and supplementary anthropogenic heat emissions [25]. In a study conducted by Sousissi et al. [26], analyzing 25 water bodies, the average annual temperature was determined to be 28 °C.
The optimal DO concentration in a healthy water body ranges from 8 to 12 mg L−1, with concentrations below 8 mg L−1 posing adverse effects on the survival of aquatic species [27]. As depicted in Table 1, values observed at the assessed points are <6 mg L−1 (excluding SW4), limiting the presence of aquatic life in general. Another factor influencing DO levels is the turbulence of the current due to air, which aids in its dissolution in water. A water body registering a dissolved oxygen level ≤ 2 mg L−1 or ≤30% saturation is categorized as hypoxic [28]. Oxygen depletion or low oxygen content typically arises from a combination of heightened biological productivity and diminished water exchange [29].
The influence of water temperature on dissolved oxygen (DO) levels, particularly in polluted conditions, has been extensively documented [30]. Elevated temperatures can exacerbate low DO levels by accelerating the decomposition of organic matter, thereby inducing alterations in ecological and environmental conditions, such as heightened eutrophication through increased nutrient discharge into sediments [31]. Furthermore, when dissolved oxygen levels plummet, fish and other aquatic organisms struggle to survive. Cold water has a higher capacity to retain oxygen, whereas warm water tends to dissolve less oxygen. In this study, no significant correlation was observed between these variables. Samples SW1 and SW7 exhibited identical DO values despite temperature variations.
It has been established that water temperature can significantly influence the amount of dissolved oxygen in water, particularly under polluted conditions. Elevated temperatures can exacerbate low dissolved oxygen (DO) levels by intensifying the decomposition of organic matter, thereby inducing changes in ecological and environmental conditions, such as an increase in eutrophication through enhanced nutrient release into sediments. Additionally, when dissolved oxygen levels plummet, fish and other aquatic organisms face survival challenges. Colder water possesses a higher capacity to retain oxygen, whereas warmer water tends to dissolve less oxygen. However, in this study, no correlation was identified between these parameters.
The samples SW1 and SW7 exhibited identical dissolved oxygen (DO) values despite temperature variations. The results obtained in this study ranged between 0.15 ± 0.07 and 6.45 ± 1.06. In a comparison between these findings and those from a study on the Bagmati River, where the dissolved oxygen fluctuated between 2.1 mg L−1 and 8.4 mg L−1, such a wide variation can be attributed to a higher organic load originating from urban areas. There exists an approximate difference of 6 units between the data obtained in this study and the findings of Giri et al. [32], where the main source of contamination is industrial.

3.2. pH and Conductivity

pH stands as a critical determinant influencing water quality and the viability of aquatic life. When pH levels deviate from the optimal range, adverse impacts on aquatic flora and fauna can occur [33], potentially indicating the presence of contaminants. The pH of most natural waters typically falls within the range of 6.5 to 8.0 [29], primarily regulated by the carbonate system. However, the ideal pH for river water should hover around 7.4 [34].
In this study, pH levels generally remained within the anticipated range for freshwater bodies, albeit exhibiting slight alkalinity or acidity at certain points (Table 2). Notably, Point SW4, situated near population (Figure 2), registered the highest pH levels, hinting at potential human-induced alterations in water quality at this locale, with variations observed from slightly acidic to slightly alkaline. It is noteworthy that certain aquatic organisms struggle to thrive in acidic environments. Except for SW6, pH values trended slightly above optimal levels, tending towards alkalinity; this point is in close relation to the companies located in the area. The presence of human settlements around SW3 may contribute to this pH elevation owing to wastewater discharge with higher pH levels than natural waters.
Therefore, sustained monitoring of pH levels within the watershed remains imperative to avert adverse impacts on aquatic life and the broader aquatic ecosystems (Figure 2). Extreme pH conditions can render a water body inhospitable to life, with low pH proving particularly detrimental to fish and immature insects. Additionally, acidic water accelerates the leaching of heavy metals, posing toxicity risks to fish.
Pure water exhibits low electrical conductivity due to the absence of ions in the solution. Conversely, river water harbors dissolved ions, rendering its conductivity proportional to the concentration of dissolved mineral salts and temperature [35]. Table 2 delineates the conductivity values documented in this study, revealing notable variations among the different sampling points. Notably, SW12 exhibited the highest and significantly disparate conductivity value compared to other points (p < 0.05). Elevated electrical conductivity in natural water bodies is primarily attributed to household and commercial waste, salts dissolved via natural runoff, regional geology, and human activities. A substantial alteration in conductivity signals potential discharge or pollution sources in a water body [36,37].
River conductivity typically fluctuates between 50,000 and 1,500,000 mS cm−1 [38]. Conversely, Rios-Villamizar [39] assessed the water quality of various Amazonian rivers (Tapajós, Incoloros streams, Cuparí, Jutaí, Tefé, and Juruá), noting diverse electrical conductivity levels, with values of 0.014331, 0.01512, 0.0599, 0.00871, 0.007367, and 0.19114 (mS), respectively. Compared to the findings of this study, these values appear diminutive.
Low conductivity (ranging from 0 to 0.2 mS cm−1) serves as an indicator of pristine or baseline conditions. Mid-range conductivity (ranging from 0.2 to 1 mS cm−1) represents the typical background for most rivers and lakes, as evidenced by the observations at SW1 to SW9. Conductivity levels beyond this spectrum could suggest unsuitability for certain animal species. Elevated conductivity (ranging from 1 to 10 mS cm−1) signifies saline conditions [40], as observed at SW12.

3.3. Total Nitrogen, NO3, and Total Phosphorus, PO4

Nitrate, a chemical compound resulting from the combination of nitrogen and oxygen, occurs naturally in trace amounts in drinking water. However, elevated nitrate levels detected in surface waters may signify contamination stemming from the presence of industrial and municipal fertilizers [41].
In this study, total nitrogen concentrations exhibited an ascending trend from SW1 to SW12, with a notable disparity (p < 0.05) among sampling points. Furthermore, it was noted that both total nitrogen and nitrate (N-NO3) concentrations, along with the NT/N-NO3 ratio, were notably higher at points SW6 to SW12. This suggests either unfavorable conditions for the oxidation of nitrogen to nitrate in the riverbed or the presence of an additional nitrogen source, such as nearby human settlements or industrial discharges proximate to the sampling sites.
Elevated nitrogen concentrations in surface waters can adversely impact the aquatic ecosystem, leading to diminished dissolved oxygen levels and toxicity to aquatic species. Moreover, intensified agricultural, livestock, and industrial activities, coupled with population growth and industrial centers within a watershed, can substantially alter nitrogen and phosphorus concentrations in water [42] (Figure 2). It is worth noting that nitrate also poses a health risk, with the permissible limit in water intended for human consumption set at 10 mg L−1 [43]; this limit was not exceeded in our study.
Phosphorus originates from various sources, both natural and anthropogenic, including soil, rocks, wastewater treatment plants, runoff from fertilized areas, and industrial activities [44]. Surface water sources of phosphorus encompass agricultural and urban runoff, soil erosion, wastewater treatment plant discharges, atmospheric deposition, and industrial operations, according to a 2021 report by the Minnesota Natural Resources Council. The report highlights that phosphorus can also leach from sediments accumulated at the water body’s bottom, subsequently fueling algal blooms and compromising water quality [45].
Orthophosphate, the oxidized mineral form of phosphorus, mirrors total nitrogen’s pattern, indicating suboptimal conditions for the oxidation of phosphorus to orthophosphate from SW1 to SW5 and from SW6 to P12 (Figure 3b).
The concentration of phosphorus in domestic wastewater typically ranges from approximately 4 to 50 mg L−1 [46]. Within this study, values falling within this range were observed, specifically at points SW6 to SW12. These findings suggest potential wastewater discharges into the river, likely attributed to the proximity of human settlements along the water body’s banks, as depicted in Figure 3b. The highest phosphate concentration was detected at SW12, surpassing values reported in a study by Fadiran et al. [47], ranging between 0.01 and 0.09 mg L−1. Upon comparison with studies conducted in rivers within the Manzini and Lubombo regions of Swaziland, it becomes apparent that the average phosphate levels exceed the recommended threshold of 2.0 mg L−1 set by the SWSC (Swaziland Water Service Corporation). Specifically, the values obtained in this study are up to 23 times higher than those reported by Fadiran et al. [47]. This discrepancy may arise from material originating from higher elevations reaching the site and being retained due to flatter relief.
Furthermore, total nitrogen levels in water tend to increase downstream, attributed to the cumulative total nitrogen transported by water from its source to its terminus. The EPA (Environmental Protection Agency) has established recommended limits of 0.05 mg L−1 of total phosphate in lake inlets and 0.1 mg L−1 of total phosphorus in runoff. It is essential to acknowledge that these limits are established to uphold water quality standards and safeguard the environment.

3.4. Total Nitrogen/Total Phosphorus Ratio

The nitrogen/phosphorus (N/P) ratio in river water serves as a pivotal indicator of aquatic ecosystem health and profoundly influences the growth of aquatic organisms such as algae, plants, and phytoplankton. In this study, the dynamics of these two parameters exhibit similar trends (Figure 4a). An imbalanced N/P ratio can precipitate excessive growth of algae and aquatic plants, culminating in eutrophication and adverse impacts on water quality and dependent organisms. Typically, an N/P ratio of 16:1 is deemed optimal for the growth of aquatic organisms in freshwater. However, N/P ratios may vary contingent on geographical location and specific environmental conditions [48], as evidenced in Figure 4b of this study.
The influx of phosphorus (P) and nitrogen (N) into freshwater systems can trigger the overgrowth of aquatic plants, oxygen depletion, and deleterious alterations in aquatic fauna diversity [49]. Total nitrogen (TN) and total phosphorus (TP) stand as primary nutrients instigating water body eutrophication [50]. Phosphate and nitrate emerge as critical limiting nutrients, fostering eutrophication in freshwater systems [50].

3.5. Chemical Oxygen Demand (COD) and Total Dissolved Solids (TDSs)

Sources of COD include leaves and woody debris, dead plants and animals, animal manure, effluent from pulp and paper mills, wastewater treatment plants, feedlots, feed processing plants, failing septic systems, and urban stormwater runoff. Discharges of wastes with high COD levels can cause water quality problems such as severe dissolved oxygen depletion and fish kills in receiving waters [51].
High COD levels not only pose water quality concerns but also raise public health apprehensions. The decomposition of organic matter in water can engender the formation of toxic organic compounds, such as trihalomethanes, which have been associated with an escalated risk of cancer and other health maladies. Thus, COD emerges as a pivotal parameter for gauging the presence of organic compounds in contaminated wastewater [52].
A study by Lee et al. [53] assessed COD levels in four rivers in South Korea, revealing values ranging from 0.8 mg L−1 to 11.6 mg L−1. In contrast, in the present study, values ranged from 0 mg L−1 to 381 mg L−1 (Figure 5a), denoting higher contamination levels in the examined samples. Another study by Perez-Castresana et al. [54] scrutinized the contamination of the Atoyac River in the metropolitan area of Puebla, Mexico, unveiling average COD values of 130.91 ± 39.52 mg L−1 in the wet season and 260.60 ± 75.32 mg L−1 in the dry season. Conversely, the average value obtained in this study was 185.06 ± 169.23 mg L−1 (Figure 5a), indicating substantial variability among sampling points and an intermediate concentration compared to the Atoyac River study. The COD concentration significantly escalated downstream (SW7-SW12). In another study that evaluated only six points, COD concentrations were documented at 118.8 mg L−1 at the first point and 138.6 mg L−1, implying mild contamination [17]. In this research, values ranged from 0 to 389 ± 11.31 mg L−1, considered severe.
The variation in dissolved solids was noteworthy (p < 0.05) for SW6, 7, 12 (Figure 5b). Dissolved solids’ composition may encompass salts, metals, metalloids, and dissolved organic matter [55]. Processes contributing to organic dissolved solids in streams involve the release of organic molecules during biological growth (e.g., by plant roots and microbes) and the decomposition of organic matter within the stream or on its banks. As per the NOM-001-SEMARNAT-1996 standard, the TSS content should be less than 150 mg L−1 for agricultural use, less than 75 mg L−1 for urban public use, and less than 40 mg L−1 to preserve aquatic life [56].
The most substantial concentration with a significant difference (p < 0.05) was noted at SW12 (refer to Figure 5b). Dissolved solids in water are a focal point in environmental research due to their potential impact on water quality and human health.
According to a study by Akpan et al. [57], average values of 156.1 mg L−1 were reported, approximately half of the average obtained in this study (350 ± 20 mg L−1). Dissolved solids may encompass salts, metals, and dissolved organic matter, stemming from biological processes such as organic matter decomposition and plant and microbe growth. In another study in the Usumacinta River, an average of 265.7 mg L−1 was documented, a value close to that found in this study.
Elevated concentrations of total dissolved solids can diminish water clarity, contributing to reduced photosynthesis and elevated water temperature. Excessively high or low dissolved solid concentrations may impede growth and lead to the demise of numerous aquatic organisms. Dissolved solids induce toxicity through salinity increments, alterations in water ionic composition, and toxicity of individual ions. Salinity rises have been shown to prompt shifts in biotic communities, curtail biodiversity, exclude less-tolerant species, and induce acute or chronic effects at specific life stages [58].

3.6. Fecal Coliforms (FCs)

According to the World Health Organization (WHO), coliform bacteria represent a group of microorganisms found abundantly in human and animal feces. Their presence in water may signify fecal contamination and, consequently, the potential existence of other associated pathogens [59]. Coliform bacteria, originating from untreated sewage, and septic tanks, among other sources, can lead to severe illnesses such as gastroenteritis or diarrhea [60]. Fecal bacteria groups such as Escherichia coli (the most commonly isolated organism in clinical laboratories [61], Klebsiella, and Enterobacter, typically found in the gastrointestinal tract of warm-blooded animals [62], serve as indicators of water quality. Their presence correlates with some of the most challenging-to-detect pathogens, including Salmonella, Shigella, and Vibrio, which can cause gastroenteritis, dysentery, typhoid, and cholera.
Table 3 delineates the distribution of fecal coliforms at the twelve study sites along the Bobos River. A discernible increase in the concentration of these microorganisms from SW1 to SW12 is observed. As depicted in Figure 2, this escalation in coliform density primarily stems from the type of discharge present at each sampling point. The study identifies a significant variation (p < 0.05) in CFU values across different sampling points (refer to Table 3). Vigilance over fecal contamination in water utilized for industrial and recreational purposes remains imperative owing to the economic and public health implications. Total and fecal coliform concentrations exhibit low levels at upstream stations and gradually escalate toward downstream stations, a trend akin to findings in [60].
These investigations underscore the significance of monitoring fecal contamination in surface waters and advocate for the implementation of control measures to curtail disease dissemination, thus safeguarding public health and the environment. Ongoing monitoring endeavors are requisite to evaluate water quality across diverse regions and institute suitable measures aimed at minimizing fecal contamination in water sources.

3.7. Quantification of NaOH (%)

Sodium hydroxide (NaOH) is a caustic substance commonly present in wastewater from drain and oven cleaners, paint and varnish removers, degreasers, dishwasher detergents, and alkaline batteries. It is utilized in the manufacturing processes of aluminum, rayon, and biodiesel and various organic synthesis procedures. Additionally, it is a frequent byproduct in chlorine production [63] and is found in pharmaceuticals such as Clintest® [64]. NaOH can significantly alter water pH [65], and its detection in drinking water may signify contamination or system malfunctions. In wastewater, NaOH detections may indicate industrial discharges or the use of cleaning products.
Proper identification and control of NaOH are paramount for preserving water quality and mitigating any adverse environmental and public health effects. However, NaOH was detected at points SW5, SW6, and SW7 in this study, with concentrations of 6.8%, 5.6%, and 3.9%, respectively. These findings could be attributed to discharges originating from settlements and industries in the study area. While the Mexican Official Standard NOM-001-SEMARNAT-2021 outlines acceptable levels of wastewater contaminants for release into national water bodies, it does not specifically address NaOH.
Concerns regarding the presence of contaminants are particularly pertinent to the local population, which utilizes these bodies of water for recreational and aquacultural purposes. Given the high level of contamination in the water, immersion in it may lead to health issues. Reports from residents suggest that contact with this tributary’s water can result in skin rashes, itching, diarrhea, and respiratory illnesses. To date, there is a lack of research evaluating these conditions, underscoring the urgency of addressing this matter to safeguard the well-being of the community.

3.8. Statistical Analysis

3.8.1. Pearson Correlation

Correlation analysis determines whether a relationship exists between two variables and the degree of association between them. A positive correlation (direct relationship) is indicated when r > 0, with values closer to 1 signifying a stronger correlation. In this study, all correlated variables exhibit positive correlations, with some displaying coefficients greater than 0.90, indicating a robust connection among them. The correlations identified include SS/temp, TN, NO3, TP, PO4, COD, and FCs, as shown in Table 4.
The examination of possible correlations between independent variables confirms a statistically significant correlation between the sampling point and the parameters listed in Table 4. The Pearson correlation coefficient reveals that as the distance from the discharges to the sampling point increases downstream, there is a corresponding rise in parameter concentrations. This suggests that pollutants evaluated are transported by the current and not naturally purified. The correlation between SS and PO4 may be attributed to agricultural areas along the river, where diffuse pollutants, such as fertilizers and livestock feces, introduced by rainfall, contribute to coliform bacteria appearing and increasing in conducive aquatic environments [60].
A correlation between temperature and TN, TP, PO4, COD, and FCs was observed, consistent with the findings of Tee et al. [66], who noted a linear increase in chemical oxygen demand (COD) with rising temperatures. Sources of nitrate and phosphorus include wastewater treatment plants, runoff from fertilized lawns and agricultural lands, faulty septic systems, runoff from animal manure, and industrial waste discharges [67]. Additionally, according to [68], nitrate concentration (NO3-N) increases with stream temperature.
Conductivity (EC) and total dissolved solids (TDSs) are parameters used to characterize water salinity levels. These two parameters are correlated and are typically expressed by the equation TDSs = k EC (at 25 °C). The relationship is influenced not only by salinity content but also by material content [69,70]. Furthermore, a positive correlation was identified between TN/NO3, TP, PO, COD, and FCs in this study, similar to the findings reported by Seo et al. [60].

3.8.2. Principal Component Analysis (PCA)

Principal component analysis (PCA) involves creating linear combinations of variables. In terms of factor loadings, a high loading is determined to be greater than 0.75, while a moderate loading falls within the range of 0.40 to 0.75. Loadings below 0.40 are regarded as insignificant according to the referenced study [71]. These techniques aim to condense the number of variables into a limited set of indices, such as principal components or factors, while endeavoring to maintain the relationships inherent in the original dataset [72,73].
The results of applying PCA for pollutant concentrations in Bobos River are reported in Table 5. According to the results of the initial eigenvalues, the three principal components were considered, which account for over 94% of the total variance. Principal component analysis for different pollutants (Table 5) is presented.
The first component (PC1), with 61.23% of the total variance, shows strong and positive loadings related to TN, NO3, TP, COD, SW, and FCs, and Pearson analysis results indicated a strong correlation between metals (TN, NO3, TP, and COD). Industry is also a significant source of anthropogenic contamination. Anthropogenic sources of COD are significant, arising mainly from industrial activities such as agri-food [74].
In PC1, a close correlation among parameters, including SW (sampling point), suggests that the sampling point exerts influence over the measured parameters. While temperature and DO may not emerge as predominant variables, ref. [75] underscores their significance, as they impact myriad life processes of organisms, along with various abiotic factors within the ecosystem.
The second component (PC2), with 21% of the total variance, is a positive weighting factor of Cond and SD. This correlation was also shown in the Pearson’s correlation.
The third component (PC3) solely comprised DO, with 11% of the total variance. The DO did not correlate with any other parameter, but the presence of oxygen plays a critical role in the transformation of nitrogen and phosphorus compounds in river water [76], underscoring the importance of monitoring oxygen concentration variability in water as a tool for recommending the sustainable use of aquatic environments [52]. While it was anticipated that a correlation or clustering would be discovered between temperature and the presence of FCs in certain components, such a relationship was not observed in this analysis. This implies that FC growth may be associated with nutrients and water pH.
The primary contaminants influencing rivers worldwide include untreated sewage, the proliferation of harmful algal blooms, biodiversity degradation, and oxygen depletion due to high concentrations of chemicals. Industrial wastewater, sewage from metropolitan areas and scientific laboratories, and surface runoff from urban areas constitute the principal sources of urban wastewater [77].

3.9. The Environmental Impact of Martinez de la Torre, Veracruz, on the Bobos River

Demographic growth, urban development, and human activities in surrounding areas have significantly impacted the environment, especially in the upper parts of watersheds, where industry, population, and rainfall directly affect urban areas. The world is urbanizing at a rapid pace, and it is essential to observe the interaction between urban centers, watersheds, rivers, and slope areas to prevent risks and ensure the water supply is suitable for human consumption, which is crucial for achieving sustainable development. An example of the effects of demographic and urban growth can be observed in Martínez de la Torre, Veracruz, where as early as 1970, there was a growing and cumulative housing problem due to population growth [78].
In warm climates such as Veracruz, a population of more than 150,000 inhabitants consumes approximately 350 L of water per day and produces 262.50 L of wastewater per day. The main sources of pollution are the population living along the river (municipal wastewater) and the citrus industry in the municipality, whose wastewater is discharged into the studied water body. In 2020, the population of Martínez de la Torre will reach 108,842 inhabitants, while Tlapacoyan will have 61,337 inhabitants, San Rafael 30,351 inhabitants, and Nautla 10,130 inhabitants, totaling 210,660 inhabitants that will affect the study area [79]. These inhabitants generate wastewater and, in turn, use water from the basin and the studied water body, which can have a significant impact on water quality. Given a daily water consumption of 350 L per person in a population of this size, it is estimated that approximately 38,069,700 L of water is consumed per day, generating approximately 262.50 L of wastewater per person per day in this region. Therefore, it is calculated that the population of Martínez de la Torre generates about 28,542,750 L of wastewater per day.
After China, Mexico is the second country in the world that uses the most wastewater in agricultural activities and is the Latin American country that irrigates the most hectares with untreated wastewater [80]. An example of this is Mexico City, where wastewater is discharged into the Tula River basin, providing water to an area that is scarce due to climatic conditions [81,82]. Although this provides economic benefits to farmers in the region, the same cannot be said for the study area, as wastewater discharged into the Bobos River basin flows into the open sea, and its effects have not yet been evaluated.
The results obtained show that there is an accumulative effect in this area since nature is not able to clean the existing pollutants in time and space. This phenomenon is visible in the points SW6-SW12, where the highest concentration of pollutants is recorded downstream. According to the Water Quality Index (WQI) developed by [83], the presence of COD > 200 mg L−1 in surface waters indicates that they have been heavily affected by both urban and non-urban wastewater discharges. In this study, values up to 389 mg L−1 were detected, indicating a situation of high pollution.
The alteration of river flows and the transport of sediment downstream are responsible for the most significant environmental impacts, endangering aquatic life and riparian ecosystems. This problem poses a threat to the survival of an entire watershed, as evidenced by various studies. Considering that the world’s population is expected to reach 9 billion by 2050, the demand for water, food, and energy is expected to increase significantly, with estimated increases of 55%, 60%, and 80%, respectively. This situation puts tremendous pressure on natural water systems, making it imperative to reduce the use of water resources. One possible solution would be to limit water consumption to 150 L per person per day, which would allow the generation of only 100 L of wastewater per person per day. This approach would help prevent the depletion of water resources and reduce the continuous discharge of large volumes of wastewater into surface water bodies.
Martínez de la Torre is a city in Veracruz, Mexico, known for its citrus production, particularly lemons and oranges [84]. Despite being known for its citrus industry, it has negatively impacted the local environment, particularly water quality and public health. In addition to contamination during the citrus production chain, the industry has been responsible for contamination with organic waste from juice exporting companies such as Citrusper, Citrex, Citrofrut, and Sicar Farms. These companies have been discharging contaminated water into the Bobos River for many years [85] (Figure 6a). However, it was not until 2014 and 2015 that PROFEPA initiated administrative proceedings against Citrofrut, S.A. de C.V., and Cítricos Ex, S.A. de C.V., primarily for irregularities in wastewater discharges and the characterization of the sludge generated in the treatment plant.
So far, there is no additional information on the characterization of the Bobos River. However, it is known that the inhabitants of this community are economically dependent on aquaculture activities. They deplore this serious pollution, which is gradually reducing the fish population. The local population points out that the citrus companies located in this region are the main culprits for the death of crustaceans and fish because they dispose of polluting liquids in the tributary of the river [86]. The production of these companies generates wastes that harm the environment, especially the water body of the Bobos River (Figure 6b), streams, and aquifers near the areas where they are disposed of. In addition, these wastes have a significant negative impact on the local flora and fauna, affecting fishermen and residents of the area. The pollution also has a significant impact on the air and soil quality in this industrial area and is not limited to the city’s waters.
To address these problems, more sustainable policies and practices are needed in citrus production. Some solutions include the implementation of cleaner agricultural practices and the promotion of composting methods for the disposal of organic waste generated by the citrus industry.

3.10. Preventive Measures to Avoid Contamination of the Bobos River and to Maintain Sustainable Development in the Region

The Bobos River in Veracruz is a crucial water source for the region, making it imperative to take preventive measures to avoid its contamination. Some measures that can be implemented include promoting sustainable agricultural practices, as a significant portion of pollution in the Bobos River stems from agriculture. Hence, initiatives should be undertaken to endorse sustainable practices aimed at reducing the use of chemicals and pesticides in crops. Crop rotation is also an option worth considering, as it helps diminish the accumulation of pests and diseases in the soil, even in perennial crop areas.
Large citrus orchards can introduce cover crops during fallow periods to shield the soil from erosion and enhance its health. These crops can serve as natural nutrient reservoirs, thereby reducing the need for chemical fertilizers. Moreover, integrated pest management practices can be adopted, encompassing mechanical, biological, and cultural techniques to control pests, rather than relying solely on pesticides. Organic agriculture, rooted in natural and sustainable practices, can substantially curtail the use of chemicals in farming and the pollution of the Bobos River.
Educating employees on sustainable and responsible production and waste management practices and fostering awareness regarding environmental protection and the significance of the Bobos River are commendable industry practices. Rigorous monitoring and control are indispensable to ensure compliance with environmental standards and forestall river contamination. Environmental monitoring systems and periodic audits can help identify areas for improvement.
Competent authorities must establish regulations and restrictions for industries utilizing the Bobos River for their operations. Quantitative limits should be imposed on the discharge of chemicals into the river, with penalties enforced against violators. Communities neighboring the Bobos River should have sewage treatment systems to prevent contaminated water from entering the river. Encouraging citizen participation through local communities can aid in preventing Bobos River pollution by actively engaging in river protection and care. Public involvement in cleaning and monitoring activities should be encouraged.
Pollution of the Bobos River poses a threat not only to the health of people in the region but also to the sustainability of the local ecosystem. Hence, the community must take preventive measures to avert river pollution and promote sustainable development. Some actions the community can take are as follows:
Refrain from littering the Bobos River and establish appropriate waste disposal sites.
Control the use of pesticides and fertilizers in homes and farmlands, opting for more natural methods. Promote reforestation by planting trees along the riverbanks and nearby areas to maintain water quality and protect the Bobos River ecosystem.
Reforestation is an important measure to maintain water quality and protect the Bobos River ecosystem. The community can promote reforestation by planting trees along the riverbanks and in nearby areas. Trees will help maintain water quality by absorbing nutrients and pollutants and providing habitat for aquatic organisms. Environmental education is another way to raise awareness about the importance of caring for the Bobos River and the environment in general. Environmental education programs should be established in schools and the community to raise awareness about the importance of protecting the river and the local ecosystem. In addition, river cleaning activities should be promoted to encourage community participation in the conservation of the Bobos River.
Establish environmental education programs in schools and the community to raise awareness about the importance of protecting the river and the local ecosystem. Promote river cleaning activities to encourage community participation in Bobos River conservation efforts. Implement wastewater treatment systems to reduce pollutants discharged into the river. Encourage industries to adopt cleaner and more efficient technologies to reduce emissions and pollution. These measures can help preserve the Bobos River and ensure its sustainability for future generations.

4. Conclusions

The water quality of the Bobos River in the state of Veracruz is compromised due to severe contamination, as evidenced by the obtained results and existing regulations. To safeguard the river and mitigate further negative impacts, it is imperative to regulate wastewater discharges effectively. Furthermore, additional studies are warranted to devise appropriate measures, given the limited information available on the water quality of the Bobos River. Ensuring a sustainable supply and equitable distribution of water resources while upholding sustainability mandates stringent regulatory oversight and ongoing environmental monitoring of existing policies and regulations.

Author Contributions

Conceptualization, H.R.G.-M. and L.C.S.H.; methodology, L.C.S.H.; software, S.Z.; validation, L.G.C.-H., F.Z. and G.N.; formal analysis, F.Z.; investigation, L.C.S.H., G.M.-R., M.S.-H.; resources, H.R.G.-M., L.G.C.-H., G.N.; data curation, L.G.C.-H.; writing—original draft preparation, L.C.S.H. and M.S.-H.; writing—review and editing, G.M.-R. and F.Z.; visualization, L.G.C.-H. and G.N.; supervision, M.S.-H. and G.M.-R.; project administration, M.S.-H.; funding acquisition, S.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Strokal, M.; Bai, Z.; Franssen, W.; Hofstra, N.; Koelmans, A.A.; Ludwig, F.; Ma, L.; van Puijenbroek, P.; Spanier, J.E.; Vermeulen, L.C.; et al. Urbanization: An increasing source of multiple pollutants to rivers in the 21st century. NPJ Urban Sustain. 2021, 1, 24. [Google Scholar] [CrossRef]
  2. Jiang, L.; O’neill, B.C. Global urbanization projections for the Shared Socioeconomic Pathways. Glob. Environ. Chang. 2017, 42, 193–199. [Google Scholar] [CrossRef]
  3. Diamond, J.; Altenburger, R.; Coors, A.; Dyer, S.D.; Focazio, M.; Kidd, K.; Koelmans, A.A.; Leung, K.M.Y.; Servos, M.R.; Snape, J.; et al. Use of prospective and retrospective risk assessment methods that simplify chemical mixtures associated with treated domestic wastewater discharges. Environ. Toxicol. Chem. 2018, 37, 690–702. [Google Scholar] [CrossRef]
  4. Singh, J.; Yadav, P.; Pal, A.K.; Mishra, V. Water Pollutants: Origin and Status. In Sensors in Water Pollutants Monitoring: Role of Material; Springer: Singapore, 2020; pp. 5–20. [Google Scholar] [CrossRef]
  5. Coelho, S.T. Existing Barriers for WtE in Developing Countries and Policy Recommendations. In Municipal Solid Waste Energy Conversion in Developing Countries; Elsevier: Amsterdam, The Netherlands, 2020; pp. 219–234. [Google Scholar] [CrossRef]
  6. Calizza, E.; Favero, F.; Rossi, D.; Careddu, G.; Fiorentino, F.; Caputi, S.S.; Rossi, L.; Costantini, M. Isotopic biomonitoring of N pollution in rivers embedded in complex human landscapes. Sci. Total Environ. 2020, 706, 136081. [Google Scholar] [CrossRef]
  7. Ju, H.; Li, S.; Xu, Y.J.; Zhang, G.; Zhang, J. Intensive Livestock Production Causing Antibiotic Pollution in the Yinma River of Northeast China. Water 2019, 11, 2006. [Google Scholar] [CrossRef]
  8. Wu, J.; Lu, J. Landscape patterns regulate non-point source nutrient pollution in an agricultural watershed. Sci. Total Environ. 2019, 669, 377–388. [Google Scholar] [CrossRef]
  9. Castañeda-Chávez, M.d.R.; Lango-Reynoso, F.; Navarrete-Rodríguez, G. Study on Contamination by Heavy Metals in the Cotaxtla-Jamapa Basin with Influence in the Central Zone of the Gulf of Mexico. Water Air Soil Pollut. 2020, 231, 99. [Google Scholar] [CrossRef]
  10. Chalchisa, D.; Megersa, M.; Beyene, A. Assessment of the quality of drinking water in storage tanks and its implication on the safety of urban water supply in developing countries. Environ. Syst. Res. 2018, 6, 12. [Google Scholar] [CrossRef]
  11. Xiao, J.; Wang, L.; Deng, L.; Jin, Z. Characteristics, sources, water quality and health risk assessment of trace elements in river water and well water in the Chinese Loess Plateau. Sci. Total Environ. 2019, 650, 2004–2012. [Google Scholar] [CrossRef]
  12. Gobierno del Estado de Veracruz, Diagnóstico Cuenca del Río Nautla. Available online: http://www.veracruz.gob.mx/proteccioncivil/wp-content/uploads/sites/5/2021/12/Diagn%C3%B3stico-Cuenca-del-rio-Nautla.pdf (accessed on 15 January 2024).
  13. Procuraduría Estatal de Protección al Medio Ambiente de Veracruz, Programa de Ordenamiento Ecológico de la Cuenca del Río Bobos. Available online: http://www.veracruz.gob.mx/medioambiente/poecrb/ (accessed on 30 December 2023).
  14. USEPA. Methods for Chemical Analysis of Water and Wastes; Report No. EPA-600/4-79-020; USEPA: Washington, DC, USA, 1983; p. 544. [Google Scholar]
  15. APHA. Standard Methods for Examination of Water and Wastewater; American Public Health Association, American Water Works Association, Water Pollution Control Federation: Washington, DC, USA, 1995; p. 1035. [Google Scholar]
  16. Norma Mexicana NMX-AA-030/1-SCFI-2012; Análisis de Agua—Medición de la Demanda Química de Oxígeno en Aguas Naturales, Residuales y Residuales Tratadas.—Método de Prueba—Parte 1—Método de Reflujo Abierto—(Cancela a la Nmx-aa-030-scfi-2001). Water Analysis—Determination of the Chemical Oxygen Demand, in Natural Waters, Wastewaters and Treated Wastewaters—Test Method—Part 1—Opened Reflux Method. 2012, México, DF. Available online: https://www.gob.mx/cms/uploads/attachment/file/166774/NMX-AA-030-1-SCFI-2012.pdf (accessed on 25 December 2023).
  17. Prambudy, H.; Supriyatin, T.; Setiawan, F. The testing of Chemical Oxygen Demand (COD) and Biological Oxygen Demand (BOD) of river water in Cipager Cirebon. J. Phys. Conf. Ser. 2019, 1360, 012010. [Google Scholar] [CrossRef]
  18. Dharma, A.D.; Gumiri, S.; Neneng, L.; Ardianoor; Arifin, S. Rapid Assessment of Water Quality on Environmental Health at Riverside Kahayan. Nat. Volatiles Essent. Oils 2021, 8, 21–33. [Google Scholar]
  19. Supardiono, S.; Hadiprayitno, G.; Irawan, J.; Gunawan, L.A. Analysis of River Water Quality Based on Pollution Index Water Quality Status, Lombok District, NTB. J. Penelit. Pendidik. IPA 2023, 9, 1602–1608. [Google Scholar] [CrossRef]
  20. Demetillo, A.T.; Japitana, M.V.; Taboada, E.B. A system for monitoring water quality in a large aquatic area using wireless sensor network technology. Sustain. Environ. Res. 2019, 29, 12. [Google Scholar] [CrossRef]
  21. Bozorg-Haddad, O.; Delpasand, M.; Loáiciga, H.A. Water quality, hygiene, and health. In Economical, Political, and Social Issues in Water Resources; Elsevier: Amsterdam, The Netherlands, 2021; pp. 217–257. [Google Scholar] [CrossRef]
  22. Dębska, K.; Rutkowska, B.; Szulc, W.; Gozdowski, D. Changes in Selected Water Quality Parameters in the Utrata River as a Function of Catchment Area Land Use. Water 2021, 13, 2989. [Google Scholar] [CrossRef]
  23. Rajesh, M.; Rehana, S. Impact of climate change on river water temperature and dissolved oxygen: Indian riverine thermal regimes. Sci. Rep. 2022, 12, 9222. [Google Scholar] [CrossRef]
  24. Liu, S.; Xie, Z.; Liu, B.; Wang, Y.; Gao, J.; Zeng, Y.; Xie, J.; Xie, Z.; Jia, B.; Qin, P.; et al. Global river water warming due to climate change and anthropogenic heat emission. Glob. Planet. Chang. 2020, 193, 103289. [Google Scholar] [CrossRef]
  25. van Vliet, M.T.H.; Ludwig, F.; Zwolsman, J.J.G.; Weedon, G.P.; Kabat, P. Global river temperatures and sensitivity to atmospheric warming and changes in river flow. Water Resour. Res. 2011, 47, 2. [Google Scholar] [CrossRef]
  26. Souaissi, Z.; Ouarda, T.B.; St-Hilaire, A. River water temperature quantiles as thermal stress indicators: Case study in Switzerland. Ecol. Indic. 2021, 131, 108234. [Google Scholar] [CrossRef]
  27. Zhao, N.; Fan, Z.; Zhao, M. A New Approach for Estimating Dissolved Oxygen Based on a High-Accuracy Surface Modeling Method. Sensors 2021, 21, 3954. [Google Scholar] [CrossRef]
  28. Hazen, E.; Craig, J.; Good, C.; Crowder, L. Vertical distribution of fish biomass in hypoxic waters on the Gulf of Mexico shelf. Mar. Ecol. Prog. Ser. 2009, 375, 195–207. [Google Scholar] [CrossRef]
  29. Dai, M.; Guo, X.; Zhai, W.; Yuan, L.; Wang, B.; Wang, L.; Cai, P.; Tang, T.; Cai, W.-J. Oxygen depletion in the upper reach of the Pearl River estuary during a winter drought. Mar. Chem. 2006, 102, 159–169. [Google Scholar] [CrossRef]
  30. Rajendiran, T.; Sabarathinam, C.; Panda, B.; Elumalai, V. Influence of Dissolved Oxygen, Water Level and Temperature on Dissolved Organic Carbon in Coastal Groundwater. Hydrology 2023, 10, 85. [Google Scholar] [CrossRef]
  31. Szewczyk, C.J.; Smith, E.M.; Benitez-Nelson, C.R. Temperature sensitivity of oxygen demand varies as a function of organic matter source. Front. Mar. Sci. 2023, 10, 1133336. [Google Scholar] [CrossRef]
  32. Giri, R.K.C.I.; Khadka, U.R. Water quality status in Bagmati river of Kathmandu valley, Nepal. In Ecological Significance of River Ecosystems; Elsevier: Amsterdam, The Netherlands, 2022; pp. 481–502. [Google Scholar] [CrossRef]
  33. Akhter, F.; Siddiquei, H.R.; Alahi, E.E.; Mukhopadhyay, S.C. Recent Advancement of the Sensors for Monitoring the Water Quality Parameters in Smart Fisheries Farming. Computers 2021, 10, 26. [Google Scholar] [CrossRef]
  34. Mondal, I.; Bandyopadhyay, J.; Paul, A.K. Water quality modeling for seasonal fluctuation of Ichamati river, West Bengal, India. Model. Earth Syst. Environ. 2016, 2, 113. [Google Scholar] [CrossRef]
  35. Mitryasova, O.; Pohrebennyk, V. Hydrochemical Indicators of Water System Analysis as Factors of the Environmental Quality State. In Sustainable Production: Novel Trends in Energy, Environment and Material Systems; Springer International Publishing: Cham, Germany, 2020; pp. 91–104. [Google Scholar] [CrossRef]
  36. Patel, M.P.; Gami, B.; Patel, A.; Patel, P.; Patel, B. Climatic and anthropogenic impact on groundwater quality of agriculture dominated areas of southern and central Gujarat, India. Groundw. Sustain. Dev. 2020, 10, 100306. [Google Scholar] [CrossRef]
  37. Koki, I.B.; Bayero, A.S. Assessment of Water Quality in Rivers and Lakes With Respect to Heavy Metals and General Water Quality Parameters: A Review. Int. J. Sci. Res. 2016, 4, 135–140. [Google Scholar]
  38. Mathur, A. Conductivity: Water Quality Assesment. Int. J. Eng. Res. Technol. 2015, 3, 1–3. [Google Scholar]
  39. Ríos-Villamizar, E.A.; Piedade, M.T.F.; Da Costa, J.G.; Adeney, J.M.; Junk, W.J. Chemistry of different Amazonian water types for river classification: A preliminary review. In Water and Society II; Brebbia, C.A., Ed.; Wessex Institute of Technology: Southampton, UK, 2013; pp. 17–28. [Google Scholar]
  40. Deepa, P.; Raveen, R.; Venkatesan, P.; Arivoli, S.; Samuel, T. Seasonal variations of physicochemical parameters of Korattur lake, Chennai, Tamil Nadu, India. Int. J. Chem. Stud. 2015, 3, 116–123. [Google Scholar]
  41. Carpenter, S.R.; Caraco, N.F.; Correll, D.L.; Howarth, R.W.; Sharpley, A.N.; Smith, V.H. Nonpoint pollution of surface waters with phosphorus and nitrogen. Ecol. Appl. 1998, 8, 559–568. [Google Scholar] [CrossRef]
  42. Rivas, Z.; Sánchez, J.; Troncone, F.; Márquez, R.; Ledo de Medina, H.; Colina, M.; Gutiérrez, E. Nitrógeno y fósforo totales de los ríos tributarios al sistema lago de Maracaibo, Venezuela. Interciencia 2009, 34, 308–314. [Google Scholar]
  43. Secretaría de Salud y Asistencia. Salud Ambiental, Agua Para Uso y Consumo Humano-Lí¬mites Permisibles de Calidad y Tratamientos a Que Debe Someterse el Agua Para su Potabilización. NOM-127-SSA1-1994. Available online: https://dof.gob.mx/nota_detalle.php?codigo=4866379&fecha=18/01/1996#gsc.tab=0 (accessed on 20 January 2024).
  44. Bhateria, R.; Jain, D. Water quality assessment of lake water: A review. Sustain. Water Resour. Manag. 2016, 2, 161–173. [Google Scholar] [CrossRef]
  45. Guidelines for Drinking-Water Quality: Fourth Edition Incorporating the First Addendum; World Health Organization: Geneva, Switzerland, 2017.
  46. Arteaga-Cortez, V.M.; de Postgraduados, M.C.; Quevedo-Nolasco, A.; del Valle-Paniagua, D.H.; Castro-Popoca, M.; Bravo-Vinaja, Á.; Ramírez-Zierold, J.A.; de México, M.U.N.A. Estado del arte: Una revisión actual a los mecanismos que realizan los humedales artificiales para la remoción de nitrógeno y fósforo. Tecnol. Y Cienc. Del Agua 2019, 10, 319–342. [Google Scholar] [CrossRef]
  47. Fadiran, A.; Dlamini, S.; Mavuso, A. A comparative study of the phosphate levels in some surface and ground water bodiesof Swaziland. Bull. Chem. Soc. Ethiop. 2008, 22, 197–206. [Google Scholar] [CrossRef]
  48. Liu, Y.; Zhu, Y.; Qiao, X.; Zheng, B.; Chang, S.; Fu, Q. Investigation of nitrogen and phosphorus contents in water in the tributaries of Danjiangkou Reservoir. R. Soc. Open Sci. 2018, 5, 170624. [Google Scholar] [CrossRef]
  49. Chambers, P.A.; McGoldrick, D.J.; Brua, R.B.; Vis, C.; Culp, J.M.; Benoy, G.A. Development of Environmental Thresholds for Nitrogen and Phosphorus in Streams. J. Environ. Qual. 2012, 41, 7–20. [Google Scholar] [CrossRef]
  50. He, M.; Wu, S.; Huang, B.; Kang, C.; Gui, F. Prediction of Total Nitrogen and Phosphorus in Surface Water by Deep Learning Methods Based on Multi-Scale Feature Extraction. Water 2022, 14, 1643. [Google Scholar] [CrossRef]
  51. Penn, R.M.; Pauer, J.J.; Mihelcic, R.J. Biochemical Oxygen Demand; 2017; Environmental and Ecological Chemistry—Volume 2. Edited by Aleksandar Sabljic. Available online: https://docplayer.net/21041474-Environmental-and-ecological-chemistry-vol-ii-biochemical-oxygen-demand-michael-r-penn-james-j-pauer-james-r-mihelcic.html (accessed on 26 December 2023).
  52. Grzyb, A.; Wolna-Maruwka, A.; Niewiadomska, A. The Significance of Microbial Transformation of Nitrogen Compounds in the Light of Integrated Crop Management. Agronomy 2021, 11, 1415. [Google Scholar] [CrossRef]
  53. Lee, J.; Lee, S.; Yu, S.; Rhew, D. Relationships between water quality parameters in rivers and lakes: BOD5, COD, NBOPs, and TOC. Environ. Monit. Assess. 2016, 188, 252. [Google Scholar] [CrossRef]
  54. Castresana, G.P.; Flores, V.T.; Reyes, L.L.; Aldana, F.H.; Vega, R.C.; Perales, J.L.M.; Suastegui, W.A.G.; Fonseca, A.D.; Silva, A.H. Atoyac River Pollution in the Metropolitan Area of Puebla, México. Water 2018, 10, 267. [Google Scholar] [CrossRef]
  55. Adjovu, G.E.; Stephen, H.; James, D.; Ahmad, S. Measurement of Total Dissolved Solids and Total Suspended Solids in Water Systems: A Review of the Issues, Conventional, and Remote Sensing Techniques. Remote Sens. 2023, 15, 3534. [Google Scholar] [CrossRef]
  56. Salcedo-Garduño, M.G.; Galaviz-Villa, I.; Pérez-Vazquez, A. Determining water quality of the lower basin of the Usumacinta River in Tabasco, Mexico. Agro Prod. 2022, 15, 23–31. [Google Scholar] [CrossRef]
  57. Akpan Emmanuel, F.; Akpan Veronica, M.; Inyang Udeme, U. Geoelectrical investigation of groundwater quality through estimates of total dissolved solids and electrical conductivity in parts of Akwa Ibom state, southern Nigeria. Malays. J. Geosci. 2020, 4, 32–37. [Google Scholar] [CrossRef]
  58. Weber-Scannell, P.K.; Duffy, L.K. Effects of Total Dissolved Solids on Aquatic Organisms: A Review of Literature and Recommendation for Salmonid Species. Am. J. Environ. Sci. 2007, 3, 1–6. [Google Scholar] [CrossRef]
  59. OMS (Organización Mundial de la Salusd). Guías Para la Calidad del Agua Potable; Cuarta edición; OMS: Ginebra, Switzerland, 2011. [Google Scholar]
  60. Seo, M.; Lee, H.; Kim, Y. Relationship between Coliform Bacteria and Water Quality Factors at Weir Stations in the Nakdong River, South Korea. Water 2019, 11, 1171. [Google Scholar] [CrossRef]
  61. Guentzel, M.N. Escherichia, Klebsiella, Enterobacter, Serratia, Citrobacter, and Proteus. In Medical Microbiology, 4th ed.; Baron, S., Ed.; University of Texas Medical Branch at Galveston: Galveston, TX, USA, 1996. [Google Scholar]
  62. Mishra, M.; Arukha, A.P.; Patel, A.K.; Behera, N.; Mohanta, T.K.; Yadav, D. Multi-Drug Resistant Coliform: Water Sanitary Standards and Health Hazards. Front. Pharmacol. 2018, 9, 311. [Google Scholar] [CrossRef] [PubMed]
  63. Dillon, C.P. Process Industries: Corrosion. In Encyclopedia of Materials: Science and Technology; Elsevier: Amsterdam, The Netherlands, 2001; pp. 7856–7860. [Google Scholar] [CrossRef]
  64. Gad, S.E. Lye. In Encyclopedia of Toxicology; Elsevier: Amsterdam, The Netherlands, 2024; pp. 997–1001. [Google Scholar] [CrossRef]
  65. Naylor, M.A.; Kaiser, H.; Jones, C.L.W. The effect of dosing with sodium hydroxide (NaOH−) on water pH and growth of Haliotis midae in an abalone serial-use raceway. Aquac. Int. 2012, 21, 467–479. [Google Scholar] [CrossRef]
  66. Tee, P.-F.; Abdullah, M.O.; Tan, I.A.; Amin, M.A.; Nolasco-Hipolito, C.; Bujang, K. Effects of temperature on wastewater treatment in an affordable microbial fuel cell-adsorption hybrid system. J. Environ. Chem. Eng. 2017, 5, 178–188. [Google Scholar] [CrossRef]
  67. MA. JoVE Science Education Database. Nutrients in Aquatic Ecosystems. Environmental Science. Cambridge. 2023. Available online: https://www.jove.com/es/v/10023/eutrophication-measuring-nitrate-and-phosphate-in-water (accessed on 28 November 2023).
  68. Jacobs, S.R.; Breuer, L.; Butterbach-Bahl, K.; Pelster, D.E.; Rufino, M.C. Land use affects total dissolved nitrogen and nitrate concentrations in tropical montane streams in Kenya. Sci. Total Environ. 2017, 603–604, 519–532. [Google Scholar] [CrossRef] [PubMed]
  69. Taylor, M.; Elliott, H.A.; Navitsky, L.O. Relationship between total dissolved solids and electrical conductivity in Marcellus hydraulic fracturing fluids. Water Sci. Technol. 2018, 77, 1998–2004. [Google Scholar] [CrossRef]
  70. Rusydi, A.F. Correlation between conductivity and total dissolved solid in various type of water: A review. IOP Conf. Ser. Earth Environ. Sci. 2018, 118, 012019. [Google Scholar] [CrossRef]
  71. Evansa, C.; Daviesa, T.; Wigington, P.; Tranter, M.; Kretser, W. Use of factor analysis to investigate processes controlling the chemical composition of four streams in the Adirondack Mountains, New York. J. Hydrol. 1996, 185, 297–316. [Google Scholar] [CrossRef]
  72. Ouyang, Y. Evaluation of river water quality monitoring stations by principal component analysis. Water Res. 2005, 39, 2621–2635. [Google Scholar] [CrossRef]
  73. Janekovi, F.; Novak, T. PCA—A Powerful Method for Analyze Ecological Niches. In Principal Component Analysis—Multidisciplinary Applications; InTech: Rijeka, Croatia, 2012. [Google Scholar] [CrossRef]
  74. He, Y.; Zhang, Y.; Li, T.; Peng, X.; Jia, X. High-concentration COD wastewater treatment with simultaneous removal of nitrogen and phosphorus by a novel Candida tropicalis strain: Removal capability and mechanism. Environ. Res. 2022, 212, 113471. [Google Scholar] [CrossRef]
  75. Riđanović, L.; Riđanović, S.; Jurica, D.; Spasojević, P. Evaluation of Water Temperature and Dissolved Oxygen Regimes in River Neretva; BALWOIS: Ohrid, North Macedonia, 2010. [Google Scholar]
  76. Kim, C.; Nishimura, Y.; Nagata, T. Role of dissolved organic matter in hypolimnetic mineralization of carbon and nitrogen in a large, monomictic lake. Limnol. Oceanogr. 2006, 51, 70–78. [Google Scholar] [CrossRef]
  77. Ibrahim, A.; Ismail, A.; Juahir, H.; Iliyasu, A.B.; Wailare, B.T.; Mukhtar, M.; Aminu, H. Water quality modelling using principal component analysis and artificial neural network. Mar. Pollut. Bull. 2023, 187, 114493. [Google Scholar] [CrossRef] [PubMed]
  78. Silva, G.L.I. The Housing Problem in Veracruz: The Case of Martinez de la Torre; Centro de Estudios Económicos y Sociales del Instituto de Investigaciones y Estudios Superiores Económicos y Sociales de la Universidad Veracruzana: Xalapa, Veracruz, Mexico, 1995; pp. 281–306. [Google Scholar]
  79. INEGI. Martínez de la Torre. Available online: https://www.inegi.org.mx/app/ageeml/#datos_generales/30/30081 (accessed on 20 January 2024).
  80. Garza, A.V. Agricultural reuse of wastewater from Cd. Juarez, (Chih., Mexico). In The Juarez Valley and Its Impact on Public Health; 2020; p. 10. Available online: https://respyn.uanl.mx/index.php/respyn/article/view/49 (accessed on 21 December 2023).
  81. García-Salazar, E.M. Wastewater as a generator of agricultural activity space in the Mezquital Valley, Hidalgo, Mexico. J. Contemp. Food Reg. Dev. 2019, 29, 2–34. [Google Scholar]
  82. García-Salazar, E.M.; Fuentes-Carrasco, M.E. The wastewater dispute in Mexico as a paradoxical ecological-distributive conflicts. Reg. Cohes. 2021, 11, 54–79. [Google Scholar]
  83. SEMARNAT, “D3_R_AGUA05_01. Compendio de Información Ambiental. Available online: https://apps1.semarnat.gob.mx:8443/dgeia/compendio_2020/dgeiawf.semarnat.gob.mx_8080/approot/dgeia_mce/html/RECUADROS_INT_GLOS/D3_AGUA/D3_AGUA04/D3_R_AGUA05_01.htm (accessed on 2 February 2024).
  84. Cauich, I.C.; Fernández, V.G.P.; Ascencio, F.J.; Rodríguez, L.E.S.; García, J.G.R. Análisis de los costos de producción del limón persa en el municipio de Tlapacoyan, Veracruz. Rev. Biológico Agropecu. Tuxpan 2014, 2, 10–18. [Google Scholar] [CrossRef]
  85. Heraldo de México, Citrusper Genera Contaminación. Available online: https://elheraldodemartinez.com.mx/estado/martinez-de-la-torre/63638-citrusper-genera-contaminacion.html (accessed on 20 January 2024).
  86. Sarmiento, J. Mas Noticias. El Rio Bobos Sufre de Contaminación lo Que Genera Una Mortandad de Especies Acuaticas. Available online: https://www.masnoticias.mx/el-rio-bobos-sufre-de-contaminacion-lo-que-genera-mortandad-de-especies-acuaticas/ (accessed on 20 December 2023).
Figure 1. Location of the study area and sampling site in the Bobos River basin.
Figure 1. Location of the study area and sampling site in the Bobos River basin.
Sustainability 16 01536 g001
Figure 2. Localization of the study area and industries.
Figure 2. Localization of the study area and industries.
Sustainability 16 01536 g002
Figure 3. (a) Total nitrogen and nitrate concentrations; (b) total phosphorous in the river channel evaluated.
Figure 3. (a) Total nitrogen and nitrate concentrations; (b) total phosphorous in the river channel evaluated.
Sustainability 16 01536 g003
Figure 4. (a) Ratio between NT and PT; (b) N/P ratio. TN: total nitrogen; TP: total phosphorous.
Figure 4. (a) Ratio between NT and PT; (b) N/P ratio. TN: total nitrogen; TP: total phosphorous.
Sustainability 16 01536 g004
Figure 5. (a) Concentration of COD and (b) concentration of dissolved solids.
Figure 5. (a) Concentration of COD and (b) concentration of dissolved solids.
Sustainability 16 01536 g005
Figure 6. (a) Discharges of wastewater from two major citrus companies (1 and 2); (b) environmental conditions of a transect of the Bobos River (SW5-7).
Figure 6. (a) Discharges of wastewater from two major citrus companies (1 and 2); (b) environmental conditions of a transect of the Bobos River (SW5-7).
Sustainability 16 01536 g006
Table 1. Temperature and DO concentration.
Table 1. Temperature and DO concentration.
ParameterTemperature (°C)DO
(mg L−1)
SW120.21 ± 0.280.15 ± 0.7
SW226.7 ± 0.140.25 ± 0.07
SW327.35 ± 0.073.80 ± 0.42
SW429.60.71 ± 0.126.45 ± 1.06
SW530.25 ± 0.640.30 ± 0.14
SW629.80 ± 0.42.64 ± 0.05
SE729.55 ± 0.070.70 ± 0.14
SW829.40 ± 0.420.25 ± 0.021
SW930.50 ± 0.570.15 ± 0.07
SW1030.20 ± 0.140.25 ± 0.07
SW1130.30 ± 0.140.35 ± 0.21
SW1230.35 ± 0.210.60 ± 0.14
DO: dissolved oxygen; SW: water sample point; ±: standard deviation.
Table 2. pH and conductivity at the sampled points.
Table 2. pH and conductivity at the sampled points.
SitespHConductivity (mS cm−1)
SW17.85 ± 0.020.30 ± 0.14
SW28.60 ± 0.140.27 ± 0.01
SW38.55 ± 0.070.28 ± 0.01
SW48.65 ± 0.070.27 ± 0.01
SW57.50 ± 0.140.26 ± 0.03
SW66.40 ± 0.710.72 ± 0.03
SW77.75 ± 0.080.70 ± 0.01
SW88.25 ± 0.210.34 ± 0.04
SW97.75 ± 0.040.75 ± 0.02
SW107.77 ± 0.050.54 ± 0.04
SW117.78 ± 0.050.37 ± 0.02
SW127.87 ± 0.093.81 ± 0.04
SW: water sample point; ±: standard deviation.
Table 3. Fecal coliform concentration.
Table 3. Fecal coliform concentration.
Water Sample PointFCs (UFC 100 mL−1)
SW10 × 100 ± 0 × 100
SW27.255 × 102 ± 1.0677 × 102
SW31.005 × 103 ± 5.657 × 101
SW41.215 × 103 ± 1.3435 × 102
SW58.16625 × 104 ± 7.425 × 101
SW68 × 104 ± 8.48528 × 103
SW71.175 × 105 ± 3.53553 × 103
SW81.91 × 105 ± 5.65685 × 103
SW92.2 × 105 ± 2.12132 × 104
SW102.322 × 105 ± 4525.48 × 103
SW112.385 × 105 ± 3.53553 × 103
SW123.2 × 105 ± 3.707107 × 104
FCs: fecal coliforms.
Table 4. Correlation matrix of physicochemical and bacteriologic parameters in Bobos River.
Table 4. Correlation matrix of physicochemical and bacteriologic parameters in Bobos River.
ParameterSWTemppHDOCondDSsTNNO3TPPO4CODFCs
SW1
Temp0.713 **1
pH−0.280−0.2061
DO−0.3340.0860.2531
Cond0.5580.250−0.135−0.1511
DSs0.5150.237−0.133−0.1200.991 **1
TN0.909 **0.604 *−0.451−0.3830.4600.4241
NO30.896 **0.546−0.158−0.4520.4150.3840.931 **1
TP0.949 **0.591 *−0.354−0.4080.5310.4860.985 **0.947 **1
PO40.933 **0.589 *−0.353−0.4120.4810.4390.992 **0.960 **0.996 **1
COD0.952 **0.588 *−0.189−0.4580.4700.4200.936 **0.974 **0.970 **0.970 **1
FCs0.956 **0.866 **−0.203−0.4410.617 *0.5610.911 **0.927 **0.967 **0.950 **0.968 **1
** The correlation is significant at the 0.01 level (two-sided). * The correlation is significant at the 0.05 level (bilateral). Cond: conductivity; sampling site: SW; temperature: Temp.
Table 5. Principal component analysis of evaluated variables.
Table 5. Principal component analysis of evaluated variables.
ParameterComponent
PC1PC2PC3
TN0.9310.250−0.108
NO30.9190.185−0.232
TP0.9210.326−0.146
PO40.9360.272−0.152
COD0.9220.257−0.203
SW0.9090.371−0.005
FCs0.9450.1350.216
DO−0.211−0.0240.897
Cond0.2810.957−0.024
SD0.2450.9660.001
Temp0.7400.0430.552
% of Variance61.23021.35011.890
Cumulative Variance (%)61.23082.58094.470
Cond: conductivity; sampling site: SW; temperature: Temp.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sandoval Herazo, L.C.; González-Moreno, H.R.; Carreto-Hernandez, L.G.; Zurita, F.; Nani, G.; Zamora, S.; Sandoval-Herazo, M.; Martínez-Reséndiz, G. Impact of the Citrus Industry on the Water Quality of the Filobobos River in Veracruz, Mexico. Sustainability 2024, 16, 1536. https://doi.org/10.3390/su16041536

AMA Style

Sandoval Herazo LC, González-Moreno HR, Carreto-Hernandez LG, Zurita F, Nani G, Zamora S, Sandoval-Herazo M, Martínez-Reséndiz G. Impact of the Citrus Industry on the Water Quality of the Filobobos River in Veracruz, Mexico. Sustainability. 2024; 16(4):1536. https://doi.org/10.3390/su16041536

Chicago/Turabian Style

Sandoval Herazo, Luis Carlos, Humberto Raymundo González-Moreno, Luis Guillermo Carreto-Hernandez, Florentina Zurita, Graciela Nani, Sergio Zamora, Mayerlin Sandoval-Herazo, and Georgina Martínez-Reséndiz. 2024. "Impact of the Citrus Industry on the Water Quality of the Filobobos River in Veracruz, Mexico" Sustainability 16, no. 4: 1536. https://doi.org/10.3390/su16041536

APA Style

Sandoval Herazo, L. C., González-Moreno, H. R., Carreto-Hernandez, L. G., Zurita, F., Nani, G., Zamora, S., Sandoval-Herazo, M., & Martínez-Reséndiz, G. (2024). Impact of the Citrus Industry on the Water Quality of the Filobobos River in Veracruz, Mexico. Sustainability, 16(4), 1536. https://doi.org/10.3390/su16041536

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