An Integrated Approach for the Assessment of Environmental Sustainability in Agro-Industrial Waste Management Practices: The Case of the Tequila Industry

Sustainable waste management is a key component needed to achieve multiple of the Sustainable Development Goals established by the United Nations for the year 2030. Such waste management can reduce the degradation of superficial water sources and can contribute to the sustainable and efficient use of resources. Due to their usually high nutrient and organic loads, agro-industrial wastes pose a threat to soil and water resources and are a major contributor to greenhouse gas generation. In this study, a novel approach was proposed to assess the environmental sustainability of vinasse management practices, integrating an extensive physicochemical characterization of tequila vinasses, a GIS-based weighted overlay analysis and a scenario analysis. Mathematical models were applied for the determination of discharged pollutants (nitrogen, phosphorus, BOD5), as well as greenhouse gas emissions. This methodology was applied to an environmental sustainability assessment of the management practices of tequila vinasses, which are one of the main waste streams generated by the agro-industry in the Mexican state of Jalisco and are one of the main contributors to organic and nutrient loads affecting the water quality of Jalisco’s superficial water sources. Through this integrated approach, critical regions for the management of tequila vinasses were determined and an extensive physicochemical characterization of tequila vinasses was performed and applied to assess the environmental sustainability of four management scenarios for tequila vinasses. These results can be used by decision-makers for the implementation of public policy and infrastructure for the improvement of local and regional waste management systems. Additionally, these data may be used to increase the environmental sustainability of the tequila industry. The proposed methodology has the potential to be applied to different waste-intensive regions as well as different organic waste streams for the assessment of the environmental sustainability of specific management practices considering the local and regional context.


INTRODUCTION
The management of agro-industrial waste poses a great challenge for production-intensive (and waste-intensive) regions, as the inappropriate disposal of large volumes of solid organic waste and wastewater can cause many environmental and health problems (Madurwar et al., 2013;Bedoić et al., 2019). Consequently, waste management is a key component to achieving multiple of the United Nations (UN) Sustainable Development Goals (SDG) established for 2030. The sustainable management of waste contributes to integrated water resource management for the protection and restoration of water-related ecosystems and to achieving the sustainable and efficient use of natural resources (UN, 2017).
Tequila is a fermented alcoholic beverage exclusively made from Agave tequilana Weber var. azul (DOF, 2012b). Tequila has been given a protected geographical status that restrains where the agave plants, and tequila itself, can be produced. As a result, tequila and agave plants can only be produced in certain municipalities in the Mexican states of Jalisco, Michoacán, Tamaulipas, Nayarit, and Guanajuato, as shown in Figure 1 (DOF, 2012a; CRT, 2019). Tequila is one of the most commonly consumed alcoholic beverages both in Mexico and globally (Bowen and Zapata, 2009;CRT, 2021). Mexico produces a yearly estimate of 1,501,081 tons of A. tequilana Weber var. azul over an approximate surface area of 87,000 ha. The state of Jalisco alone contributes 75% of the national agave production (SIAP, 2020). The total tequila production in 2019 reached an all-time high with 374 million liters. Seventy percent of the produced tequila is destined for exportation and the remainder is distributed nationally (CRT, 2021).
Tequila production can be divided into four stages: 1) agave processing and sugar extraction, 2) fermentation, 3) distillation, and 4) maturation and bottling (Cedeño, 1995). Tequila vinasses are the liquid sludge generated during tequila production and are usually made at a 7-15:1 ratio of vinasse to tequila volume FIGURE 1 | Distribution of tequila-producing states in Mexico.
Frontiers in Environmental Science | www.frontiersin.org June 2021 | Volume 9 | Article 682093 (España-Gamboa et al., 2011). While all four production stages have a direct effect on the quality, character and sensory characteristics of the produced tequila, only the first three can affect the physicochemical characteristics of the vinasses generated (Cedeño, 1995;López-López et al., 2010). Although vinasses are mainly generated at the distillation stage, effluents from other production stages may also be included (Cedeño, 1995). Tequila vinasses present high concentrations of total solids (TS) as a result of the sugar extraction processes, like mill extraction and diffusion extraction, prior to the fermentation stage (Íñiguez-Covarrubias and Peraza-Luna, 2007), as well as a result of extremely high chemical oxygen demand (COD) concentrations and high acidity (España-Gamboa et al., 2011). Tequila vinasses also contain high concentrations of phenolic compounds that give them a characteristic dark brown color (López-López et al., 2010). These characteristics make vinasses a great environmental concern, as their incorrect disposal can cause significant environmental impacts (López-López et al., 2010). This is especially a concern in regions with high rates of Tequila production, such as in the state of Jalisco (Macías, 1997). The high agricultural and agro-industrial activities of the state of Jalisco have led to increased exploitation and pollution of its water resources (Montoya et al., 1997;Abrol and Sharma, 2012). Eutrophic conditions caused by high nutrient and organic matter loads have been reported for multiple waterbodies across the state. This is mainly due to excessive uncontrolled solid waste and wastewater discharges by agro-industrial activities, leading to anoxic conditions and toxic cyanobacterial blooms, which can damage wildlife and even pose a threat to nearby settlements (de Anda et al., 2019;Díaz-Torres et al., 2020;Gradilla-Hernández et al., 2020). In order to reduce this eutrophication, the management practices of both solid waste and wastewater must be improved to decrease the amounts of organic matter and nutrient discharge going into Jalisco's hydrographic systems (Montoya et al., 1997;Abrol and Sharma, 2012;Jayme-Torres and Hansen, 2018). Additionally, Jalisco has officially declared an aim for the year 2030 to reduce its current GHG (greenhouse gas) emissions by 22%, of which 14% have been attributed to waste and wastewater management (INECC, 2018;SEMADET, 2018).
The evaluation of environmental sustainability of waste management practices has been accomplished based on the following inputs: public and stakeholder opinion (Sahar and Ahmad, 2020), overall waste production estimation and management costs analysis (Daian and Ozarska, 2009), combined Life Cycle Analysis (LCA) and Analytical Hierarchy Process (AHP) (Lijó et al., 2018), combined LCA and Resource Conservation Efficiency (RCE) determination (Kaufman et al., 2010), waste flow mapping (Kurdve et al., 2015), and carbon footprint assessment (Ibanescu et al., 2018), among others. However, most waste management environmental sustainability studies usually lack practical application, as they do not provide an integral approach based on the local geographic, environmental and socioeconomical conditions, while maintaining replicability for other regions or waste streams (Bani et al., 2009). New integrated methodologies are required to generate adaptative environmental sustainability assessments that can be easily applied by decision-makers to establish public policy and infrastructure to improve local and regional waste management systems, especially for wasteintensive regions.
As a contribution to the literature, this study presents a novel approach to spatially characterize the environmental impacts (nutrient and organic matter discharge and GHG production) of the tequila vinasses generated in the state of Jalisco. This work herein also quantitively assesses different management scenarios, combining the physicochemical characterization of vinasses collected in field, GIS-based multicriteria analysis and scenario analysis. Due to the lack of academic literature regarding the spatial variations of the physicochemical characteristics of tequila vinasses, an extensive physicochemical characterization of Tequila vinasses was carried out including samples from Tequila Production Units (TPU) across the state of Jalisco. In order to identify the most critical regions to establish vinasse management strategies, a weighted overlay analysis (WOA) was implemented to evaluate the potential pollution risk for the entire territory of the state, integrating five productive and environmental parameters: 1) vinasse production at the nearest tequila production unit (TPU), 2) total vinasse production for the municipality, 3) distance to the nearest TPU, 4) distance to the nearest water body and 5) distance to the nearest urban community. Furthermore, the environmental impacts generated as a result of the current management practices of tequila vinasses were evaluated in addition to the environmental benefits that could be gained if other management practices were implemented. This was accomplished using information generated through the physicochemical characterization and focusing on the critical regions previously determined through the WOA. Nutrient (nitrogen and phosphorus) and organic matter discharge estimations, as well as GHG production estimations were performed for the critical regions as well as for the entire state, based on the existing models reported for the management conditions of each scenario.

METHODS
The general methodology used in the present study is shown in Figure 2.

Data Collection and Processing
The geographic information and the annual production of tequila production units (TPU), expressed in m 3 year −1 , was obtained from the National Statistic Directory of Economic Units in Mexico (DENUE by its acronym in Spanish) (INEGI, 2021a) and validated with official records provided by the Environment and Land Development Ministry of the state of Jalisco (SEMADET by its acronym in Spanish) (SEMADET, 2019b). When possible, this information was further corroborated directly with the producers and updated. TPUs were then classified into production categories according to the internal classification system of the Tequila Regulatory Council (CRT by its acronym in Spanish). This classification was based on the tequila production volumes presented in Table 1, and production Frontiers in Environmental Science | www.frontiersin.org June 2021 | Volume 9 | Article 682093 records were used to estimate the production of vinasses for each individual TPU. The current management practices of vinasses within the TPUs in the state were also identified from official records (SEMADET, 2019b) in order to ascertain the most common practices by the production category according to the CRT classification system for the TPUs (Table 1). Information related to the hydrographic structure, runoff network, digital elevation model, official municipal limits, and location of urban areas within the study region was obtained from the National Institute of Statistics and Geography (INEGI) (INEGI, 2021b). A geographical information system (GIS) was created on ArcGis 10.5 to integrate the data corresponding to the TPUs and the previously mentioned environmental and geographical parameters. All raster images had a 15 m × 15 m resolution.

Selection of Parameters for the WOA
The critical regions for the management of tequila vinasses were determined based on vinasse production records and environmental parameters. A weighted overlay analysis allows for the integration of multiple classified raster layers into a single weighted resulting raster based on a weighted index matrix to classify the evaluated parameters (Yazdani et al., 2015;Al-Anbari et al., 2018). A weighted index matrix was generated to classify and rank the rasters comprising each individual point within the territory of Jalisco according to their potential risk of environmental degradation as a result of inadequate management of vinasses. The following parameters were transformed into raster format for the weighted overlay analysis for the establishment of the critical regions: 1) vinasse production at the nearest TPU, 2) total vinasse production for the municipality, 3) distance to the nearest TPU, 4) distance to the nearest water body, and 5) distance to the nearest urban community. A similar approach was reported by (Díaz-Vázquez et al., 2020). These parameters were selected as they are commonly used for the evaluation of the pollution risk as well  as the spatial suitability for the implementation of management strategies and infrastructure for organic waste and wastewaters (Comber et al., 2015;Mukherjee et al., 2015;Venier and Yabar, 2017;Díaz-Vázquez et al., 2020). These parameters are further explained in the following paragraphs.
The rate of vinasse production at each individual TPU (V) was determined based on the mean, minimum, and maximum production ratio of vinasses (liter of vinasses per liter of tequila) (R) previously reported in the literature ( Table 2) and the tequila production volume reported for each individual TPU (T), as shown in Eq. 1. The mean production rates of vinasses were used for the WOA based on the mean production ratio of vinasses (R mean ) reported in Table 2.
TPUs can act as both point and non-point sources of pollution to water sources, as the generated vinasses can be discharged directly from the TPUs into water bodies and municipal wastewater systems or they can be transported for field application as fertilizers and then be transported to water bodies by rainfallrunoff processes (Fuess and Garcia, 2014). Therefore, the vinasse production rate at the nearest TPU has a great influence on the pollution risk at every point in the territory.
(2) Total TPUs within the municipality The total count of TPUs for each municipality in Jalisco was determined based on the TPU inventories. The TPU density within each municipality was included in the analysis as it is related to the density of vinasses that require to be disposed of within the municipal area. Additionally, according to current legislation, all waste of special management is the responsibility of the state, including vinasses (DOF, 1988;SADER, 1989). Therefore the state has the authority to set specific guidelines and regulations for the management of special wastes, as well as the right to establish strategies and region-specific programs in order to improve management of high priority waste, like vinasses (DOF, 1988;SADER, 1989). Higher densities of TPU can facilitate the implementation of pollution prevention strategies within specific regions or small clusters of TPUs rather than at a state level. Additionally, high densities of point sources can lead to the accumulation of pollutants in the areas nearby. Therefore, it is important to locate possible high-density regions in order to establish preventive measures and regulations to reduce pollution risk (Nobre et al., 2007).
(3) Distance to the nearest TPU The distance from every point in the territory in Jalisco to the nearest TPU was determined using the Euclidean distance tool in ArcGIS 10.5. High production rates of vinasses, and organic waste in general, imply elevated transport costs due to the high water content of vinasses. As a consequence, vinasses tend to be disposed of and accumulated in areas nearby the TPUs in an attempt to reduce disposal costs (Carral et al., 2011;Shahmoradi and Isalou, 2013;Fuess and Garcia, 2017;Díaz-Vázquez et al., 2020). This leads to soil degradation, waterbody eutrophication and, in extreme cases, health hazards for nearby populations (Fuess and Garcia, 2014;Chowdhary et al., 2018).
(4) Distance to the nearest waterbody The distance from every point in the territory in Jalisco to the nearest water body or water flow were determined using the Euclidean distance tool in ArcGIS 10.5 based on the hydrological data available for the state of Jalisco (INEGI, 2021b). The disposal of vinasses can have a direct effect on the water quality of nearby superficial water bodies (Carral et al., 2011;Shahmoradi and Isalou, 2013). Therefore, it is important to determine how close the evaluated point in the territory is to the nearest superficial waterbody.
(5) Distance to the nearest urban community The distances from every point in the territory of Jalisco to the nearest urban area were determined using the Euclidean distance tool in ArcGIS 10.5 based on the geographical data available for the state of Jalisco (INEGI, 2021b). The disposal of vinasses to local sewage systems can overload the treatment capacity of local wastewater treatment plants due to the high concentration of pollutants in vinasses. While this is a common management strategy used by TPUs located within urban areas (Beltrán et al., 2001), current legislation requires that each TPU is responsible for the management and disposal of their wastewaters and should not transfer this responsibility to municipal wastewater treatment plants (DOF, 1988;SADER, 1989). Additionally, the accumulation of untreated vinasses in soils and water bodies can generate foul odors due to organic matter degradation (Carral et al., 2011;Shahmoradi and Isalou, 2013).

Development of Classified Raster Images
Raster images were constructed for the five selected parameters. The raster pixels within each raster image of each selected parameter were classified into non-dimensional values ranging from 1 to 5 according to the ranges presented in Table 3. These values were assigned based on how each of the selected parameters increases the risk of environmental pollution from the inadequate management of vinasses at every point in the territory of Jalisco. 1 was assigned to the areas with the lowest pollution risk, and 5 represents the highest pollution risk (Al-Anbari et al., 2018;Díaz-Vázquez et al., 2020).
To evaluate the overall risk of pollution as a result of the inadequate management of vinasses, at every point in the territory, weights were assigned to each parameter based on the literature and studies that evaluate similar parameters (Yager, 1988;Al-Hanbali et al., 2011;Shahmoradi and Isalou, 2013;Vaissi and Sharifi, 2019;Díaz-Vázquez et al., 2020). A WOA was performed using the five generated classified raster images according to Eq. 2, where RF represents the risk factor for every individually analyzed pixel (1-5 scale), W i represents the assigned weight for the selected parameter and S i represents the reclassified raster value of the selected parameter (1-5) for the analyzed pixel. Weights were distributed evenly between all parameters based on previous studies (Yager, 1988;Al-Anbari et al., 2018;Vaissi and Sharifi, 2019;Díaz-Vázquez et al., 2020). The WOA was performed on ArcGIS 10.5 using the Weighted Overlay Analysis tool.

RF
The resulting classified raster, with a 15 m × 15 m resolution, represents the overall pollution risk in each point of the territory (each pixel on the map) considering all five selected parameters expressed on a 1-5 scale. Contiguous municipalities showing zones with a pollution risk of 5 were classified as critical regions.

Collection of Vinasses and Physicochemical Characterization
The physicochemical characterization of vinasses was required due to the lack of information regarding the spatial variability of tequila vinasse physicochemical characteristics in the existing literature. Vinasses were collected from 20 TPUs distributed across the previously identified critical regions. A stratified sampling was performed to include all TPU size categories within the study. 5 L grab samples were taken the same day that these were produced, transported in refrigeration (4°C) for less than 5 h and kept at 4°C until further processing within the next 24 h, as established by the applicable regulation (DOF, 1997). No composite samples were prepared as all samples were analyzed individually in duplicate. Hydrogen potential (pH) was measured at the collection site by the potentiometric method while total suspended solids (TSS), settleable solids (SS), biological oxygen demand (BOD5), total Kjeldahl nitrogen (TKN), total phosphorus (TP), and fat oil and grease (FOG) were determined upon arrival by the Analytical and Metrologic Services Unit (USAM by its by its acronym in Spanish) of the Centro de Investigación y Asistencia Tecnológica del Estado de Jalisco (CIATEJ). TSS and SS were determined by gravimetric methods, BOD5 was determined by the respirometric method, TKN was determined using the Kjeldahl method, TP was determined by the persulfate digestion-molybdovanadate method, and FOG was determined by the Soxhlet extraction method. All determinations were made following the Federation, W. E., and American Public Health Association methods (Eaton et al., 2005).
A physicochemical characterization was required for the evaluation of the management scenarios. One-way ANOVAs and Tukey tests were performed on Minitab 19.2020.1 to determine significant spatial differences between the sampled vinasses collected from different regions of the state of Jalisco considering each of the water quality parameters (pH, BOD 5 , FOG, TKN, TP, and TSS). The same analyses were performed to determine possible differences in the vinasses produced by different TPU size categories ( Table 1).  Table 4.

Management Scenarios Analysis
The potential nutrient and organic matter release into the environment, as well as the potential GHG emissions (CO 2 eq) Vinasses production at the nearest TPU to the site m 3 year −1 0-1,000 1,000-3,000 3,000-10,000 10,000-30,000 >30,000 20 Total TPUs in the municipality Units 0-1 1 -5 5 -10 10-15 >15 20 Distance from the point in the territory to the nearest TPU M 0-2,500 2,500-5,000 5,000-7,500 7,500-10,000 >10,000 20 Distance from the point in the territory to the nearest superficial waterbody or water flow M 0 -250 250-500 500-750 750-1,000 >1,000 20 Distance from the point in the territory to the nearest urban area M 0-500 500-1,000 1,00-1,500 1,500-2,000 >2,000 20 Frontiers in Environmental Science | www.frontiersin.org June 2021 | Volume 9 | Article 682093 from the disposal of the vinasses within the critical regions, were determined based on the characterization of vinasses performed in the present work for the four proposed scenarios (A, B, C, and D).

Nutrient Release to the Environment
The total vinasse production, TN production, and TP production were estimated based on the total production volumes reported for each individual TPU and the physicochemical characterization of the sampled vinasses using the minimum, mean and maximum values (for both vinasse production rates and physicochemical parameters) as shown in Table 1 and Table 6, for each of the critical regions. These estimations are total nutrient productions without the differentiation of the bioavailability. The total nitrogen (N T ) (Mg year −1 ) potential release was determined according to Eq. 3 for each scenario. This was based on the total vinasses generated (V) (m 3 year −1 ) in each critical region and the mean, minimum, and maximum nitrogen concentration in the sampled vinasses (N) (Mg m −3 ) displayed in Table 6. The nitrogen load (N D ) (Mg year −1 ) in the discharged effluent was determined according to Eq. 4, where R N (nondimensional) represents the nitrogen removal rate under each scenario.
The total phosphorus (P T ) (Mg year −1 ) potential release was determined following the same method used for the calculation of the total nitrogen potential release, as shown in Eq. 5, where V is the total vinasse generation (m 3 year −1 ) in each critical region and P is the phosphorus (mean, minimum, and maximum) composition of the sampled vinasses (Mg m −3 ) as displayed in Table 6. The phosphorus load in the discharged effluent (P D ) (Mg year −1 ) was determined according to Eq. 6, where R P (Nondimensional) represents the phosphorus removal rate under each scenario.

Organic Matter Release and Greenhouse Gas Emissions
The mean, minimum and maximum organic matter content present in the untreated vinasses and the discharged effluent for each respective scenario was calculated according to Eqs. 7, 8, respectively. BOD T (Mg year −1 ) represents the total BOD 5 content in the untreated vinasses according to BOD 5 (Mg m −3 ) values shown in Table 4. These were obtained as a result of the physicochemical characterization of the sampled vinasses. BOD D (Mg year −1 ) represents the remaining BOD 5 present in the discharged effluent after the proposed treatment, according to the specific organic matter removal factors of each process as reported in the literature ( Table 4).
The mean, minimum and maximum GHG emissions (GHG i ) (Mg CO 2 eq year −1 ) generated though the proposed scenarios were calculated for each scenario (A, B and C, and D, respectively) according to Eqs. 9-11, using the mean, minimum and maximum values from the vinasse production rates and the physicochemical characterization. These equations were adapted from the IPCC Guidelines for National Greenhouse gas inventories based on process specific factors, summarized in Table 4 (Doorn et al., 2006). The factors proposed by the IPCC are region specific in most cases and, thus, the Latin-America values were used when available. When not available a mean parameter was used, as suggested by the IPCC (De Klein, 2006;Doorn et al., 2006).
B o (nondimensional) is the most commonly used methane production capacity in the calculation of GHG inventories, which indicates the fraction of organic matter that can be transformed to CH 4 and tends to be independent of the degradation process and its conditions, assumed as 0.6 (Doorn et al., 2006). MCF (nondimensional) represents the methane conversion factor used to estimate the conversion of BOD 5 into methane under specific degradation conditions (Doorn et al., 2006). V (m 3 year −1 ) represents the generation of vinasses in the individual TPU. N T (Mg year −1 ) represents the nitrogen concentration in the untreated effluent, and N D (Mg year −1 ) represents the nitrogen present in the discharged effluent after the proposed treatment. The fraction of nitrogen transformed and released into the atmosphere as N 2 O is expressed as EF (nondimensional). Both CF CH4 (nondimensional) and CF N2O (nondimensional) represent the conversion factors used to transform the CH 4 and N 2 O into CO 2eq , respectively. Eqs. 9-11 were adapted from the equations proposed by De Klein, (2006) and Doorn et al. (2006), as they are the most common methods for estimating GHG resulting from waste management and agricultural practices (De Klein, 2006;Doorn et al., 2006). The GHG estimation for scenarios C and D includes only the direct GHG generated by organic matter and nutrient degradation; it does not include possible GHG reductions due to energetic need reduction or even energetic surplus due to the energetic revalorization of CH 4 .

Characterization of the Tequila Production Sector and of Current Producer Management Practices
The total production of tequila in Jalisco is accomplished by 158 TPUs, from which 109, 17, 14, and 18 units are categorized as micro, small, medium, and large, respectively. This is according to the internal classification system used by the CRT, which is based on the tequila production volume shown in Table 1. The mean estimated vinasse production rate in Jalisco was found to be 2'645,655.69 m 3 year −1 . The TPUs classified as large contribute 84.3% of the state vinasse production, while the remaining 15.7% is generated by micro, small and medium TPUs contributing with 1.7, 4.1, and 9.8% respectively. Almost 80% of the state's vinasse production is accomplished by only seven municipalities: Tequila, Zapotlanejo, Atotonilco el Alto, Amatitan, Arandas, Tlaquepaque, and Tototlan. The complete vinasse production and size distribution for the municipalities responsible for 80% of the vinasses within the state is displayed in Table 5. Tequila and Amatitan are located within the middle-western region of the state, Atotonilco, Arandas and Tototlan belong to the eastern region, and Zapotlanejo and Tlaquepaque are located in the central region of the state. Most small, micro, and medium TPUs used pH neutralization and direct field application as their management practices for vinasses, while large TPUs used a variety of treatment alternatives, mainly categorized as biological treatment (aerobic or anaerobic) and, with a lower frequency, sequential anaerobic-aerobic treatment.

Critical Region Determination
Five classification maps were generated, each corresponding to one of the selected parameters ( Figure 3) used to evaluate the risk of environmental degradation due to inadequate management of tequila vinasses. High values (5), which indicated a higher risk of pollution, were located within the middle-western and eastern regions of the state of Jalisco. This can be attributed to the higher density of TPUs, and their high tequila production rates, within these regions as well as the abundance of urban settlements across these areas.
As a result of the WOA, a classification map was obtained to identify the critical regions ( Figure 4) weighting all five selected parameters. Two main critical regions were found to be located in the western region of the state (Valles region) and in the eastern region of the state (Altos region). These critical regions, Valles and Altos, have the highest possible risk of environmental pollution as a result of nutrient and organic matter release and the highest potential to produce GHG emissions resulting from poor management of the tequila vinasses. The contiguous municipalities displaying critical zones (Figure 4) were selected and clustered, as shown in Figure 5. The Valles critical region comprises the municipalities of Tequila, Amatitan, and El Arenal, while the Altos critical region consists of the municipalities of Tepatitlan de Morelos, Atotonilco el Alto, Arandas, and Jesus Maria, as displayed in Figure 5. The municipalities of Tequila, Amatitan, Atotonilco el Alto, and Arandas belong to the municipalities responsible for 80% of vinasse production in the state, as shown in Table 5.

Characterization of Vinasses
The results of the physicochemical characterization, summarized by descriptive statistics, of vinasses collected at 20 TPUs distributed across the critical regions are displayed in Table 6. Additionally, based on the current Mexican legislation (DOF, 1997), the newly proposed Mexican legislation (DOF, 2018), the European Union water directive (Council directive, 2014) and the EPA guidelines for centralized waste treatment point discharges (EPA, 2021), the permissible limits are shown to compare the results to the applicable federal discharge guidelines as well as to international standards. While the physicochemical characteristics of vinasses from other alcoholic production processes (molasses, sugarcane, sugar beet, etc.) are well documented (España-Gamboa et al., 2011;Christofoletti et al., 2013;Chowdhary et al., 2018;Mikucka and Zielińska, 2020), there is a lack of information regarding tequila vinasses because most of the existing reports on the physicochemical characteristics have been obtained from single production units and thus fail to represent the variability in vinasse characteristics resulting from a wide variety of productive processes used across the state (Íñiguez-Covarrubias and Peraza-Luna, 2007;López-López et al., 2010).
The pH values of the sampled vinasses exhibited the lowest coefficient of variation (CV 5%) while FOG displayed the highest coefficient of variation of all analyzed parameters (CV 158%). This indicates that variations within the tequila production process have a great influence on the FOG content of the discharged vinasses. TP and SS also presented elevated coefficients of variation, at 156 and 133%, respectively.
According to Table 6, vinasses surpass the permissible discharge limits established in the current Mexican legislation for all of the regulated parameters. Furthermore, neither of the vinasse samples were within acceptable limits for any of the seven analyzed parameters according to current legislation. Therefore, a treatment process would be required in order to fulfill the current local legislation standards (DOF, 1997). The local permissible limits for wastewater discharge are currently under a revision process since a new legislation was recently proposed (DOF, 2018). New lower permissible limits may be established within the short term due to the new legislation proposal, which would increase the need for effective wastewater management strategies in order to fulfill the new standards that not only reduce current permissible limits but add new parameters to be regulated, as is the case of color and toxicity (DOF, 2018). Both international discharge guidelines presented in this study ( Table 6) exhibit lower limits than the current local legislation, making it unacceptable to discharge untreated vinasses under those guidelines (Council directive, 2014; EPA, 2021).
The mean concentrations of phosphorus and nitrogen (42.01 mg L −1 and 131.63 mg L −1 ) were closer to the current permissible limits, 40 and 110% higher than the current legislation, in comparison to the mean concentrations for other parameters. However, when compared to international standards, tequila vinasses display concentrations of phosphorus and nitrogen that were 2,100 and 400% higher than the permissible limits set for the European Community (Council directive, 2014). BOD 5 is the water quality parameter that most highly exceeded the current local and international guidelines, reflecting the extreme organic load in the vinasses and emphasizing the need for effective vinasse management practices in order to prevent environmental degradation.
From the seven analyzed parameters, significant mean differences between the vinasses produced by TPUs of different size categories were only found when considering SS (p 0.011). No significant differences were found between the vinasses produced in the different regions (Altos and Valles) of the state of Jalisco, indicating there is no difference between the vinasses generated in either of the tequila producing regions.

Critical Regions for the Management of Vinasses
All of the tequila producers registered in the official records in Jalisco were classified in four categories, as shown in Table 1. Similar proportions of all categories (micro, small, medium, and large) were found for both critical regions, Valles and Altos (SEMADET, 2019b) While most of the vinasses generated in Jalisco are produced by large TPUs (83%), a great majority of these producers already report vinasse treatment systems that offer complete or partial treatment. Large producers generally have the resources and technical knowledge to implement effective wastewater management systems. However, medium, small, and micro-producers generally lack these resources and technical knowledge and thus release higher volumes of untreated vinasses discharged to soils and waterbodies across the state (SEMADET, 2019b).
Both the critical regions determined in this study displayed similar overall vinasse production volumes, as shown in Table 7. However, these regions present different spatial configurations. Most of the TPUs within the Valles critical region are in a smaller area (within a 5 km radius, approximately), which favors the implementation of centralized treatment units. This is not only because of the constant supply of vinasses, but due to the relatively short distances between most TPUs. The TPUs within the Altos critical region are more dispersed across the region (within a 30 km radius), which can be a disadvantage for the implementation of large-scale centralized treatment units due to higher transport costs (Flotats et al., 2009;Fuess and Garcia, 2014). Areas with lower densities of wastewater production may require different approaches in order to achieve effective and economical wastewater management like small-scale clustering or other diverse onsite treatment options (Massoud et al., 2009).
Besides being the top Tequila producer, Jalisco is the leading state for the production of crops and livestock in Mexico (SIAP and SAGARPA, 2018). This agricultural productivity increases the overall organic waste and nutrients released to the environment as a result of inadequate management practices (Montoya et al., 1997;Martinez Rivera, 2004;Jayme-Torres and Hansen, 2018;Díaz-Vázquez et al., 2020). Both critical regions identified in this study, Altos and Valles, display high agricultural and industrial activity that influences the soil and water quality of the region. Additionally, the tequila industry in these regions generates larger volumes of agro-industrial waste compared to other industrial categories and therefore has a great influence on the environmental degradation (Montoya et al., 1997;SEMADET, 2017;Jayme-Torres and Hansen, 2018). Furthermore, these regions are located within the Santiago river basin, which is known to be highly polluted with nitrogen, phosphorus, and organic matter (Abrol and Sharma, 2012). The spatial identification of critical regions can be a powerful tool to identify areas that require special policies or management in order to prevent environmental degradation or health hazards to nearby populations (Delivand et al., 2015;Díaz-Vázquez et al., 2020). Although it is out of the scope of the present work, the identification of waste-intensive regions would open up the possibility for the implementation of synergic waste management systems. Different waste streams could then be integrated and revalorized to higher levels than the ones from individual revalorization of waste streams. One example of this is the anaerobic co-digestion of waste, where the biogas yields from mixed waste stream feeds are significatively higher than the ones from individual waste feeds (Muhammad and Chandra, 2021;Rincón-Pérez et al., 2021).

Potential Environmental Effects of Untreated Vinasses
As shown in Table 6, the mean pH of the analyzed samples is within the range previously reported for tequila vinasses (3.4-4.5) (López-López et al., 2010). If vinasses are discharged without previous pH neutralization, the receptive soils and water bodies can be severely affected due to the high acidity of this agroindustrial waste (Fuess and Garcia, 2017). Acidification of soil and water is known to increase nutrient mobility, altering its uptake in crops as well as increasing its transport by runoff or leaching to receptive waterbodies or groundwater, respectively (World Health Organization, 2006). The most common method used by TPUs to regulate this parameter is to utilize neutralizing agents, such as CaO, CaOH 2, and NaOH. However, special care must be taken, as large amounts of neutralizing agents can increase the content of specific ions in the effluent (Fuess and Garcia, 2017;Chowdhary et al., 2018), which over time can cause soil degradation, heavy metal transport, as well as damage to aquatic species due to increased ion concentrations in receptive Frontiers in Environmental Science | www.frontiersin.org June 2021 | Volume 9 | Article 682093 waterbodies (Dickman and Christy, 2002;Fuess and Garcia, 2014;Fuess and Garcia, 2017). The mean BOD 5 value ( Table 6) in the sampled vinasses (119,011.77 ± 104,097.12 mg L −1 ) was greater than previously reported values for tequila vinasses as well as for other vinasses (20,600-60,000 mg L −1 ) (López-López et al., 2010; España-Gamboa et al., 2011). The higher BOD 5 concentrations in the sampled vinasses can be linked to the high TSS content also found in the characterized vinasses (9,875.80 ± 11,958.41). This is due to the fact that TSS in tequila vinasses are mainly composed of organic unfermentable solids like cellulose and pectin, which are present in the agave juice that was not degraded during the alcoholic fermentation stages (Íñiguez-Covarrubias and Peraza-Luna, 2007; España- Gamboa et al., 2011). The mean TSS content in the sampled vinasses (9,875.80 mg L −1 ) is above the higher limit reported by López-López (2000-8,000 mg L −1 ), possibly due to the mechanical extraction processes applied by the specific TPUs. Mechanical sugar extraction processes, such as mill extraction, increase the total unfermentable solid content in agave juice, which is discarded in the vinasse effluent during the distillation stage (Íñiguez-Covarrubias and Peraza-Luna, 2007;Tequila Sauza, 2021). The implementation of other sugar extraction processes, such as diffusion extraction, can reduce the solid load in the produced vinasses as well as increase sugar availability during the fermentation stage (Tequila Sauza, 2021). The TSS content in the vinasses represents an important fraction of a vinasse DBO 5 and can directly have negative effects in both soils and receptive waterbodies if discharged without previous treatment. High BOD 5 values can significantly reduce the dissolved oxygen levels in receptive waterbodies due to microbial proliferation (EPA, 2001). Over time, this increases GHG production due to the production of metabolic byproducts like CH 4 and CO 2 (Doorn et al., 2006). Significant differences were identified in the SS mean values. Micro and small producers generate lower volumes of vinasses compared to medium and large producers and usually store them for longer periods prior to treatment or disposal (SEMADET, 2019b). Therefore, the settleable solid content in the sampled vinasses could have been affected due to natural sedimentation during the time prior the sampling.
The TKN content (Table 6) in the analyzed vinasses (131.63 ± 97.59 mg L −1 ) was approximately three times higher than the previously reported values for tequila vinasses (20-50 mg L −1 ). However, it was still lower than the values reported in the literature for vinasses from mezcal production, which uses different varieties of agave; likewise, the TKN content reported herein was lower than in vinasses obtained from the fermentation of other substrates like corn, beet, and sugarcane (660-5,650 mg L −1 , 2,000-2,680 mg L −1 , 1,800-4,750 mg L −1 , and 1,000-7,000 mg L −1 , respectively) (Robles- González et al., 2012;Fuess and Garcia, 2017). An important fraction of the total nitrogen content in tequila vinasses is in an organic form, which reduces its potential as a crop fertilizer since it requires mineralization by soil microorganisms to be of use for crop fertilization (Fuess and Garcia, 2017). Organic nitrogen is mineralized in soils or water bodies, into ammonium which can be easily used by plants as a nitrogen source. However, high ammonium concentrations may lead to elevated concentrations of nitrites and nitrates in soil and waters, causing eutrophication, or may also cause increased N 2 O emissions through nitrification or denitrification processes (Kirchman, 2018). The mean phosphorus content of the tequila vinasses analyzed in the present study (42.01 mg L −1 ) corresponds to the values reported for tequila vinasses by España-Gamboa (40 mg L −1 ) (España-Gamboa et al., 2011). Tequila vinasses present lower phosphorus contents compared to vinasses obtained from other common substrates (71-1990 mg L −1 ) (sugar cane, sorghum, sugar beet, molasses, etc.), which can be attributed to the natural phosphorus content in these other substrates (España-Gamboa et al., 2011). Because phosphorus exhibits lower volatility compared to nitrogen, there is the risk of accumulation of phosphorus in soil and receptive waterbodies if vinasses are inadequately managed (World Health Organization, 2006) either by the direct discharge of or by the transport of accumulated phosphorus in the soil through superficial runoff and groundwater (Fuess and Garcia, 2017). The disposal of untreated vinasses by field irrigation or direct discharge to water bodies can be a powerful driver for phosphorus accumulation due to its excessive phosphorus content and reduced volatilization rates, leading to eutrophication of water bodies due to surface runoff and groundwater phosphorus transport (Carpenter et al., 1998;Fuess and Garcia, 2017). The FOG mean values ( Table 6) obtained in the present study (301.14 mg L −1 ) were 300% higher than the upper limit of the reported range by López-López for tequila vinasses (10-100 mg L −1 ) (López-López et al., 2010). The high coefficient of variation ( Table 6) (CV 158%) found for FOG in this study is of special interest since the lower measured values of FOG in the analyzed vinasses fall within the acceptable limits for both current and proposed legislation. The cooking process used for the agave heads may be related to the high variation in the FOG values measured in the analyzed vinasses since the agave heads cuticles release waxy compounds under heat treatment. Prior to sugar extraction by milling, agave heads are cooked mostly either by oven cooking or autoclaving (Cedeño, 1995). The use of diffusion extraction, on the other hand, may prevent the release of waxy compounds into the sugar-rich agave juice and can significantly reduce the content of FOG in the vinasses FIGURE 3 | Classification maps generated for the selected parameters.
Frontiers in Environmental Science | www.frontiersin.org June 2021 | Volume 9 | Article 682093 (Cedeño, 1995;Tequila Sauza, 2021). However, only a small fraction of producers have implemented extraction alternatives like diffusion (Cedeño, 1995;Tequila Sauza, 2021). When vinasses are inadequately managed and released to the environment untreated, the FOG present in them not only contributes to the overall BOD 5 but can reduce gaseous transfer rates in both water bodies and soils and can cause clogging in transport and treatment piping (Husain et al., 2014). The differences presented between previous studies and the present work can be attributed to the larger pool of sampled vinasses since the present study characterized 20 different TPUs across the state of Jalisco compared to individual characterizations performed by Lopez-Lopez, Íñiguez-Covarrubias, and España-Gamboa (Íñiguez-Covarrubias and Peraza-Luna, 2007;López-López et al., 2010;España-Gamboa et al., 2011), which may not represent the complete diversity of TPUs and their vinasses. López-López et al. (2010) reported the only physicochemical characterization of Tequila vinasses available in the literature but neither indicated the amount/ location of TPUs sampled nor the descriptive statistics from the results.

Critical Regions and Scenarios Analysis
Scenario A (field application) is based on the most common management in Jalisco, while scenarios B (aerobic treatment), C (anaerobic treatment), and D (sequential anaerobic-aerobic treatment) are not widely implemented in Jalisco but are generally applied for the management of vinasses in other regions (López-López et al., 2010;España-Gamboa et al., 2011).

Scenario A: Field Irrigation
Because scenario A (field application) assumes no treatment previous to the vinasse's disposal, the complete nitrogen, phosphorus, and organic matter loads are released to the soil through irrigation (348.24 Mg year −1 , 111.14 Mg year, −1 , and 314,864.18 Mg year −1 , respectively) ( Table 8). Although the estimated GHG emissions for scenario A (473,823.03 Mg year −1 ) ( Table 8) are the lowest of the proposed scenarios, the GHG production rates in scenario A are highly dependent on the environmental conditions (temperature, oxygen, nutrient concentration, microorganism populations, etc.) on which the organic matter and nutrients present in the soils and waterbodies are degraded. The estimated GHG emissions for scenario A are likely an underestimation since the emission factors used assume complete aerobic degradation. Moreover, due to the elevated organic content, vinasses are prone to generate partially anaerobic conditions upon applications, which would significatively increase overall GHG estimations (De Klein, 2006;Doorn et al., 2006;Fuess and Garcia, 2014). Additionally, GHG estimations from scenarios C and D do not include possible GHG reductions due to external energy demand reduction, which would reduce overall GHG estimations, lowering the GHG estimation difference between scenario A and scenarios C and D. The mean GHG production estimated for scenario A represents 82.2% of the total GHG production estimate through wastewater management reported for the state of Jalisco (SEMADET, 2019a) and 1.5% of the state's total GHG emissions (SEMADET, 2019a). However, due to lack of information, it is not clear if the GHG resulting from wastewater management previously estimated for Jalisco includes only urban wastewaters or if it also includes agro-industrial wastewater, therefore explaining the low reported wastewater management GHG values in the state. Field application, proposed in scenario A, is the most common disposal method of vinasses globally (Moraes et al., 2017). As shown in Table 8, field irrigation of vinasses has several financial advantages, such as low investment and overall disposal costs, in addition to savings in chemical fertilizer or other fertilizer sources   for crops. However, this management practice also comes with strong disadvantages that reduce its environmental sustainability due to possible environmental impacts generated. Plants mainly assimilate nitrogen in the form of nitrate. Therefore, most of the nitrogen (85-90%) present in vinasses cannot be directly assimilated by plants, reducing the fertilization potential of untreated vinasses (Fuess and Garcia, 2017). Increased nitrogen content in untreated vinasses leads to high volatilization rates by denitrification increasing N 2 O generation (World Health Organization, 2006). Phosphorus displays higher stability compared to nitrogen and can be accumulated on the soil´s surface if regularly applied, making it prone to be transported by runoff during rain events (World Health Organization, 2006). Vinasses typically do not contain significant heavy metal concentrations, but the increased soil acidity caused by the regular application of vinasses can increase the mobility of metals present in the soil, causing the accumulation of metals in receptive waterbodies and groundwater (Fuess and Garcia, 2017). To reduce transport costs, a common practice for tequila producers is to dispose of the vinasses in their own crop fields (SEMADET, 2019b). Since agave is one of the most cultivated crops within both of the critical regions, a great majority of the irrigated vinasses are disposed of in agave fields. However, agave plants display higher soil erosion rates compared to other crops due to their rainfall-concentrating leaf structure, which leads to higher nutrient and solid transport in runoff coming from agave fields (Abrol and Sharma, 2012). Therefore, counterproductive effects can be seen due to excessive vinasse application on crops, as proposed in scenario A. This is due to changes in soil structure and salinity, nutrient bioavailability and microbial communities (World Health Organization, 2006;Fuess and Garcia, 2017).
Jalisco has a total cultivated surface of 1,649,784.9 Ha in 2019 (SIAP, 2020), with an estimated fertilizer consumption of 54.79 kg Ha −1 and 25.74 kg Ha −1 of nitrogen and phosphorus, respectively (FAO, 2021). While the fertilizer requirements of Jalisco largely surpass nitrogen and phosphorus production from tequila vinasses, which represent only 0.03 and 2.5% of the total nitrogen and phosphorus demand respectively, the production of vinasses in the Valles and Altos critical regions and the distribution of the generated vinasses for its use in other regions of the state by fertigation may be economically unfavorable due to increased transport costs (Fuess and Garcia, 2014).

Scenario B: Aerobic Treatment
The aerobic treatment of wastewater (scenario B) is one of the most common processes applied globally for the treatment of organically polluted wastewater since it has proven to be an effective way to remove both organic matter and nutrients from organic wastewaters (SEMADET, 2019b). A downside of this process is the intensive energy requirement to achieve the aeration rates required for the system to operate properly (Rantala and Väänänen, 1985;Eckenfelder, 1998). This energetic demand increases overall treatment costs (Rantala and Väänänen, 1985;Cakir and Stenstrom, 2005). From the four proposed scenarios, the aerobic treatment (scenario B) had the highest GHG emission estimates with a mean of 1,451,781.65 Mg year −1 while having the second lowest nitrogen and phosphorus disposed volumes (208.94 and 61.28 Mg, respectively) ( Table 8). GHG production from aerobic systems tends to be proportional to the organic load of the wastewater feed (Cakir and Stenstrom, 2005). Aerobic treatment of highly polluted wastewaters, like vinasses, tends to demand elevated oxygen concentrations in order to allow for the aerobic degradation of organic matter. If operated under non-ideal conditions (low aeration, low mixing, low retention times, etc.), partial anaerobic conditions can be generated within the treatment system, which increases GHG release (Doorn et al., 2006). Since none of the CH 4 generated by the aerobic process is captured and combusted, the complete CH 4 load is released to the atmosphere, increasing overall GHG estimations for the scenario and increasing the GHG production difference against other proposed management scenarios (Díaz-Vázquez et al., 2020).

Scenario C: Anaerobic Treatment
Scenario C had a mean GHG estimate of 511,225.04 Mg year −1 ( Table 9), which represents a 96% reduction compared to the values estimated for scenario B. This was due to CH 4 combustion and transformation into CO 2 rather than direct release as CH 4 . Although scenario C displayed lower GHG and a significant reduction of the discharged nitrogen and phosphorus load (208.94 Mg year −1 and 61.28 Mg year −1 , respectively) ( Table 9), the implementation of individual anaerobic systems requires high investment from producers, limiting its implementation to medium and big producers (Anderson et al., 2013). The implementation of centralized anaerobic digestion systems can offer the possibility to micro and small producers to treat their vinasses through greener alternatives, which bypass the required high initial investments but do not significatively increase treatment costs (Díaz-Vázquez et al., 2020). However, the operation of a centralized anaerobic digestion system can present difficulties due to the heterogeneous nature of the vinasses generated by multiple production units as well as the presence of other pollutants due to non-ideal wastewater management systems within the production units. Alterations in the feed composition can cause changes in the microbial populations within the anaerobic system leading to reduction in treatment efficiency and even complete shutdown of the system (Naik et al., 2014). This heterogenicity in the vinasses can be regulated by strictly monitoring the incoming vinasses in order to prevent undesired pollution and to maintain stable feed conditions to prevent inhibitory or toxic effects on microbial population within the anaerobic digester (Álvarez et al., 2010). The use of effective nutrient (nitrogen and phosphorus) recovery systems, as proposed in scenarios B, C, and D, not only prevents the excessive discharge of nutrients into the soil and water but additionally presents the opportunity to reduce treatment cost by valorizing the recovered nutrients and allowing for a controlled application of nutrients in crop fields. By concentrating the recovered nutrients in both aerobic and anaerobic sludge, the volume of the nutrient mixture is reduced, simultaneously reducing transport and application costs, which in time allows these recovered nutrients to be applied in fields further from the production site, increasing its economical competitivity against other fertilizer sources (Fuess and Garcia, 2014;Díaz-Vázquez et al., 2020).

Scenario D: Sequential Anaerobic-Anaerobic Treatment
Scenario D displayed the maximum nutrient (nitrogen and phosphorus) and BOD 5 reduction and the third lowest GHG

Field irrigation (A)
• Reduced disposal costs Fuess and Garcia, (2014) • The elevated water content of vinasses limits their transportation for them to be used as fertilizer which leads to its accumulation near production sites Fuess and Garcia, (2017) • Does not require a large initial investment in infrastructure or equipment Fuess and Garcia, (2014) • Potential for soil acidification as well as waterbody eutrophication due to nutrient transport in runoff and groundwater Fuess and Garcia, (2014); Fuess and Garcia, (2017) • Low operation complexity Fuess and Garcia, (2014) • Low nitrogen fertilizer potential due to high content of organic nitrogen Mengel, (1996) • If applied according to crop needs, vinasses can be used as a cheap source of nutrients Christofoletti et al. (2013) • Since there is not BOD 5 reduction the complete nutrient and organic load of the vinasses is discharged into soil leading to its uncontrolled degradation and the production of GHG Fuess and Garcia, (2017) • Can help reduce the state dependency on chemical and imported fertilizer Kanagachandran and Jayaratne, (2006) • • High nutrient and organic load reduction Chan et al. (2009) • High volumes of sludge are generated with lower nutrient content compared to anaerobic sludge and require additional processing, if not used as fertilizer requires additional treatment processes in order to prevent further pollution Field et al. (1984); Nguyen and Shima, (2019) • Sludge can be used as fertilizer Nguyen and Shima, (2019) • Aerobic treatment increases overall GHG generation due to higher CH 4 production rates compared to field degradation Cakir and Stenstrom, (2005) • Does not require temperature control Chan et al. (2009) • Its operation increases treatment costs significantly manly due to aeration needs Rosso et al. (2008) • Short start-up periods Chan et al. (2009) • High initial investment required for the construction of the treatment infrastructure • Does not require alkali addition for its proper operation Beltrán et al. (2001) • Typically requires large extensions for the construction of the needed infrastructure Chan et al. (2009) • Low operation complexity -Anaerobic treatment (C) • Negligible N 2 O generation Cakir and Stenstrom, (2005) • Nutrient and organic load removal rates lower than those of aerobic treatment systems Chan et al. (2009) • Higher nutrient content and bioavailability in generated sludge Field et al. (1984) • High initial investment required for the infrastructure required for both anaerobic treatment and CH 4 revalorization Anderson et al. (2013) • CH 4 capitation systems prevent its release into the atmosphere Cakir and Stenstrom, (2005) • Anaerobic systems require long startup periods Chan et al. (2009) • CH 4 can be used as an energy source either for thermal or electric energy conversion Rantala and Väänänen, (1985) ;Eckenfelder, (1998) • Anaerobic microorganisms present high sensitivity to feed changes and alteration in the vinasses feed can alter the removal rates as well as the CH 4 production within the system McLeod et al. (2015) • Low operational costs especially for wastewater with high organic loads Rantala and Väänänen, (1985) ;Eckenfelder, (1998) • Feeds presenting extreme pH values can reduce the productivity of the system and require pH control Chan et al. (2009), which is the case for vinasses • Energetic costs can be greatly reduced and even completely covered by the energetic production using CH 4 Rantala and Väänänen, (1985); De Bere, (2000) • Require temperature control in order for the system to operate • Higher nutrient and orgnic load removal rates than both individual treatment processes Frostell, (1983) • Higher GHG emission than the individual anaerobic treatment due to the addition of the aerobic stage Keller and Hartley, (2003); Cakir and Stenstrom, (2005) • Maximum nutrient recovery and revalorization potential by the production of nutrient-rich sludge Cervantes et al. (2006) • Highest investment required for the implementation of the infrastructure needed for the operation of the anaerobic-aerobic sequential treatment Chan et al. (2009) • Excess aerobic sludge can be reprocessed into the anaerobic system Frostell, (1983) • Anaerobic systems require long startup periods Chan et al. (2009)  Mg, respectively). The higher GHG production through scenario D in comparison to scenario C can be explained by the addition of an aerobic stage, which slightly increases both CH 4 and N 2 O production even when the BOD 5 and nutrient loads entering the aerobic system are reduced by the previous anaerobic stage (Cakir and Stenstrom, 2005). Anaerobic digestion of effluents with high organic loads tends to result in lower overall GHG generated compared to other treatment alternatives, even reaching negative values due to the use of produced biogas for energy generation and the consequent reduction of required energy during the treatment process, contrary to aerobic treatment processes in which wastewaters with a high organic load tend to increase overall GHG generation and energy consumption (Cakir and Stenstrom, 2005). The GHG reduction by CH 4 energetic revalorization is not accounted for in scenarios C and D, meaning both have the potential to produce even lower GHG emission rates. Table 10 displays a comparison of the reported advantages and disadvantages for each of the management strategies proposed in the present work.

Integrated Management Approaches
The implementation of effective treatment systems, either aerobic (scenario B), anaerobic (scenario C), or sequential anaerobic-aerobic (scenario D), usually requires large investments, which limit the producers that can afford such investments and the corresponding operation costs (Díaz-Vázquez et al., 2020). The implementation of centralized treatment units allows micro and small productive units to access effective treatment systems reducing the initial investment needed for its implementation. Additionally, TPUs located within urban areas, which is the case of most productive units located within the municipality of Tequila (Valles critical region), allows access to treatment systems that would require offsite construction due to spatial constraints, reducing the required investment and reducing the strain of local wastewater treatment plants caused by vinasse discharge into municipal wastewater systems (Flotats et al., 2009;Massoud et al., 2009).
The logistics required for the operation of centralized treatment units can significatively increase treatment costs, especially in waste with elevated water content, which is the case of tequila vinasses. This needs to be taken into consideration by decision-makers in order to prevent an unnecessary treatment cost increase due to logistics (Díaz-Vázquez et al., 2020). The establishment of centralized treatment units within critical regions allows the guarantee of a constant supply of wastewater to maintain the economic viability of the treatment systems and protect the natural resources and services of sensitive regions from degradation, which is the case for most basins within Jalisco (de Anda and Shear, 2001;Martinez Rivera, 2004;Abrol and Sharma, 2012;Jayme-Torres and Hansen, 2018). This is not only applicable to tequila vinasses but for most primary production wastes and byproducts generated in the state, as is the case of livestock production waste and dairy processing wastewater (Díaz-Vázquez et al., 2020).
Since vinasses are generated by most alcoholic beverage and bioethanol production processes, regions in Brazil, Japan, the United States as well as many countries within the European Union that are known for intensive alcoholic production, face similar problems to those found in Jalisco (España-Gamboa et al., 2011). Several other practices have been proposed globally for the management of vinasses, in addition to the ones proposed in the previously discussed scenarios. Some examples include use as animal feed, treatment by electrocoagulation, coagulation-flocculation and advanced oxidation (Ozone) as well as biorefinery approaches using either bacteria, fungi, or algae for the simultaneous treatment and recovery of metabolites, biomass or energy (Waliszewski et al., 1997;López-López et al., 2010;España-Gamboa et al., 2011;Christofoletti et al., 2013;Mikucka and Zielińska, 2020). Several novel processes have been developed on a laboratory scale that may offer new alternatives with promising effects on the quality of generated vinasses, as they use hydrolytic enzymes to increase substrate fermentability and decrease the discarded vinasses organic load (Ávila-Fernández et al., 2009;Waleckx et al., 2011). These alternative management practices could be analyzed following the methodology proposed here to assess their environmental sustainability applied to Tequila vinasses. It would be important to pay special attention to the scale and economic viability for the implementation of such practices, as well as to the local and regional context.
The existence of strong producer associations and certification organizations for TPUs at state and national levels can assist in the establishment of centralized treatment units. They can act as intermediaries between policymakers, local authorities, academia, technology suppliers and producers in order to establish the best management practices and coordinate efforts to improve the Tequila production sector and achieve international environmental standards. 10 | (Continued) Comparison between the advantages and disadvantages of the proposed vinasses management scenarios.

Management scenario Advantages Disadvantages
• Energy recovery by CH 4 revalorization by its conversion into thermal or electric energy can help reduce operating costs from both the anaerobic and the aerobic stages of the treatment process Chan et al. (2009) • High operation complexity and higher operation costs

CONCLUSION
The proposed methodology integrates both physicochemical characterization and GIS-based spatial analysis in order to evaluate management scenarios for the determination of possible environmental impacts derived from management practices, either current or proposed. This can be a powerful tool for policymakers and local authorities to evaluate and propose local or regional integral waste management programs, policies and infrastructure. It likewise has the potential to incentivize sustainable practices among productive sectors. This approach can be applied to other regions or other waste-producing sectors in order to determine best management practices in a local and regional context. The proposed methodology could be improved by the generation of site-specific and waste-specific emission, removal, and conversion factors in order to estimate nutrient and organic load discharge as well as greenhouse gas emissions with higher accuracy based on regional and geographical characteristics, either for environmental (soil and water) degradation of waste or for the evaluation of specific treatment alternatives. It would also be useful to include the integration of multiple organic waste streams into the analysis model. This would open up the possibility to evaluate coprocessing alternatives that may prove to have a more positive environmental impact compared to individual management strategies, especially for dense and highly productive regions, such as the Mexican state of Jalisco.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

FUNDING
This research was funded by the Innovation, Science and Technology Ministry of the State of Jalisco (SICYT by its acronym in Spanish) and the Science and Technology State Council (COECyTJAL by its acronym in Spanish) through Jalisco's Scientific Development Found (FODECIJAL by its acronym in Spanish), assigned as the Project 7947: "Vinazas tequileras como sustratos para la producción de aditivos alimenticios de alto valor proteico para animales monogaśtricos."