key: cord-0786790-a7t7t83b
authors: Coronado, Yaxk’in; Navarro, Roberto; Mosqueda, Carlos; Valenzuela, Valeria; Pérez, Juan Pablo; González-Mendoza, Víctor; de la Torre, Mayra; Rocha, Jorge
title: SARS-CoV-2 in wastewater from Mexico City used for irrigation in the Mezquital Valley: quantification and modeling of geographic dispersion
date: 2021-08-09
journal: Environ Manage
DOI: 10.1007/s00267-021-01516-4
sha: 1ad43b511ef2d45fc92ee8fb5e7ab89e3255657f
doc_id: 786790
cord_uid: a7t7t83b

Quantification of SARS-CoV-2 in urban wastewaters has emerged as a cheap, efficient strategy to follow trends of active COVID-19 cases in populations. Moreover, mathematical models have been developed that allow the prediction of active cases following the temporal patterns of viral loads in wastewaters. In Mexico, no systematic efforts have been reported in the use of these strategies. In this work, we quantified SARS-CoV-2 in rivers and irrigation canals in the Mezquital Valley, Hidalgo, an agricultural region where wastewater from Mexico City is distributed and used for irrigation. Using quantitative RT-PCR, we detected the virus in six out of eight water samples from rivers and five out of eight water samples from irrigation canals. Notably, samples showed a general consistent trend of having the highest viral loads in the sites closer to Mexico City, indicating that this is the main source that contributes to detection. Using the data for SARS-CoV-2 concentration in the river samples, we generated a simplified transport model that describes the spatial patterns of dispersion of virus in the river. We suggest that this model can be extrapolated to other wastewater systems where knowledge of spatial patterns of viral dispersion, at a geographic scale, is required. Our work highlights the need for improved practices and policies related to the use of wastewater for irrigation in Mexico and other countries.

The ongoing global pandemic of COVID-19 disease, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is a public health emergency of international concern (Organization and Fund (UNICEF) 2020a, 2020b). SARS-CoV-2 ribonucleic acid (RNA) has been detected in feces from both symptomatic and asymptomatic patients (Chen et al. 2020; Holshue et al. 2020; Jiehao et al. 2020; Tang et al. 2020; W. Wang et al. 2020; Zhang et al. 2020 ) and in wastewater (Ahmed et al. 2020; Lodder and Husman 2020; Medema et al. 2020) . For this reason, quantification of SARS-CoV-2 RNA in wastewater has emerged as a cheap, efficient method for monitoring active cases in large populations (Ahmed et al. 2020; S. Wang et al. 2020) , small towns (Kitajima et al. 2020; Randazzo et al. 2020) , or campuses (Harris-Lovett et al. 2021) . Notably, since both symptomatic and asymptomatic individuals contribute to the viral load in wastewaters (Vallejo et al. 2020; Wu et al. 2020) , this strategy allows for * Jorge Rocha jorge.rocha@ciad.mx 1 Conacyt-Unidad Regional Hidalgo, Centro de Investigación en Alimentación y Desarrollo. Ciudad del Conocimiento y la Cultura de Hidalgo, Blvd. Santa Catarina S/N, San Agustín Tlaxiaca, Hidalgo, México42163 a one-week anticipation in the active cases, compared to health systems registries; i.e. an increase in recently infected individuals who remain asymptomatic may not be registered by health systems, but can be detected in wastewater viral loads (Lesimple et al. 2020; Medema et al. 2020; Harris-Lovett et al. 2021) . Mexico City has more than 661,446 confirmed cases of COVID-19 (Mexico City Government 2021; consulted on June 1st, 2021); this is the city with the highest number of cases in the country. In Hidalgo, a state north of Mexican Valley Metropolitan Area (MVMA, Fig. 1a ), more than 39,012 positive cases have accumulated (State of Hidalgo Government 2021), of which more than 11,000 cases (28%) correspond to the Mezquital Valley, a highly productive agricultural region.

To date, no systematic efforts have been reported about detection of SARS-CoV-2 in wastewater in Mexico.

However, the Mezquital Valley has several characteristics that make the study of SARS-CoV-2 in wastewater relevant. This region is located 80 km north of Mexico City (Fig. 1) . It is inhabited mainly by the indigenous Otomí group that speaks a native language called hñähñu. The average annual precipitation is 409 mm (Moreno Alcántara et al. 2006 ) and the main economic activity is agriculture, but the arid climate contributes to the extreme poverty, historically associated to the region (SAGARPA 2003; Instituto Nacional de Lenguas Indígenas (Mexico) 2009; Rossette 2017) . For this reason, wastewater is widely used for irrigation; it is received from Mexico City through the West Interceptor Tunnel, Central Interceptor Tunnel, and Grand Drainage Canal (Fig. 1b) . Wastewater is partially processed in the wastewater treatment plant (WWTP) in Atotonilco de Tula, Hidalgo before feeding the Tula River and a complex system of irrigation canals (Fig. 1c) (Contreras et al. 2017 This represents a health and environmental risk (Lesser et al. 2018) , as farmers, inhabitants and consumers in the Mezquital Valley are in contact with water, soil or agricultural products. However, the use of wastewater for irrigation is one of the main causes that allowed that this Valley is no longer in extreme poverty (García-Salazar and García-Salazar 2019).

Using data for SARS-CoV-2 concentration in wastewater, several models have been proposed for to find a correlation between the temporal patterns of viral concentration and the number of active COVID-19 cases (Hart and Halden 2020a) . Likewise, a previous work presented the first model of spatial and temporal patterns of viral loads in seaway systems, showing the dispersion of the virus in an urban region (Hart and Halden 2020b) . However, models for spatial patterns of viral dispersion are needed in open waterbodies that contain wastewater, in order to understand the virus transport and to identify the zones with higher risk of infection.

In this study we sampled water from the Tula River, Salado River and irrigation canals in the Mezquital Valley ( Fig. 1) , which receive wastewater from Mexico City, to assess the presence of SARS-CoV-2 and generate a mathematical model that describes its spatial dispersion. We propose that our results are highly relevant not only to follow the epidemic in Mexico City and municipalities in Hidalgo, but also to evaluate a possible risk of transmission through environmental matrices and measure the stability of SARS-CoV-2 with geographic resolution.

Between September 22nd and October 6th 2020, we collected water samples in eight locations of the Tula River, representing the main flow of wastewater from Mexico City in the Mezquital Valley, this also included two locations in the Salado River, which receives wastewater from the Grand Drainage Canal (sewage stream indicated with a yellow arrow in Fig. 1b) , and Tepeji river, which contains wastewaters from local municipalities (Fig. 1c, Table 1 ). These samples were collected between the WWTP and Ixmiquilpan ( Fig. 1) . Water samples were also collected from irrigation canals adjacent to agricultural fields as well as soil samples from these agricultural fields (Fig. 1c , Table 1 ). Produce samples were also collected at three of the agricultural fields: a sample of coriander at site 3, and samples of lettuce at sites 6 and 7. In total, our sampling covered ≈80 km in the Mezquital Valley, and included the municipalities of Tepeji del Río, Atotonilco de Tula, Tula de Allende, Tezontepec de Aldama, Mixquiahuala de Juárez, Tlaxcoapan, Tetepango and Ixmiquilpan (Table 1) .

Appropriate safety equipment was always used, consisting of a cotton lab coat, latex gloves, KN95 respiratormask, rubber boots, disposable cap, and safety glasses. A simple sampling technique was used (NOM, 1980; EPA, 2017) , locating sites where the wastewater is well mixed near to the center of the flow channel, approximately between 40 and 60 percent of water depth, where turbulence is maximum, and there is minimum sedimentation of solids. For each location, three water samples were collected as follows: two samples of 400 ml in glass bottles, for SARS-CoV-2 detection and for microbiological analyses; and one 4 L sample in a plastic bottle, for physicochemical analysis. In the case of samples for SARS-CoV-2 detection, a single sample was collected in each location of the Tula River, while duplicate samples were collected in each location for irrigation canals. For each water sample collected at an irrigation canal, we also sampled 500 g of soil from an adjacent agricultural field. All samples were immediately placed in ice until arrival at the laboratory. Water samples for SARS-CoV-2 detection were inactivated by incubation at 60˚C for 1 h upon arrival at the laboratory.

Samples for physicochemical and microbiological analyses were stored in ice from sampling until delivery to external laboratories where determinations were carried out. Due to the remoteness of sampling sites, analyses were performed between 18 to 40 h after sampling. Water samples were analyzed for physicochemical parameters at Laboratorio de Análisis Químicos Analíticos (LAQA), at Facultad de Química, Universidad Autónoma de Querétaro (Querétaro, México). Physicochemical parameters measured were: Biochemical oxygen demand (BOD), chemical oxygen demand (COD); pH, total dissolved solids (TDS), suspended solids (SS), total suspended solids (TSS) and total solids (TS). Microbiological analysis of total coliforms (TC), fecal coliforms (FC), Escherichia coli and Salmonella (PCR detection) were analyzed (according to norm NOM-210-SSA-2014) at Laboratorio para la Evaluación y Control de Riesgos Microbianos en Alimentos (LECRIMA) at UAQ (Querétaro, México). The 'most probable number' method was used for quantification of TC, FC and E. coli.

Duplicate RNA extractions were performed for each sample collected for SARS-CoV-2 detection. For soil samples, viral particles were first eluted from 15 g of sample in 150 ml of 10% beef extract buffer (pH 7) by stirring for 30 min (Horm et al. 2012 ). Next, soil suspensions and inactivated wastewater samples (150 ml) were processed following previously reported protocols (Kocamemi et al. 2020; Kumar et al. 2020; Wu et al. 2020 ) with slight adaptations. First, samples were filtered through a No. 5 Whatman filter (Millipore), and then through 0.22 µm polyethersulfone membrane (Corning).

Next, viral particles were precipitated by adding NaCl (100 g/L) and Polyethilenglicol 6000 (22 g/L) and stirring for 20 min. Then, 105 ml samples were centrifuged at 12,000 × g for 2 h, and the resulting pellet was suspended in 400 µl RNAse free water. Viral nucleic acids were extracted from these samples using the Purelink Viral DNA/RNA kit (Invitrogen), with an elution volume of 50 µl. Purified nucleic acid samples were quantified by spectrophotometry in a Nanodrop One (Thermo Scientific). RNAs from soil samples were further processed with the OneStep PCR Inhibitor Removal Kit (Zymo Research).

For detection of SARS-CoV-2, the kit Decov2 Triplex (Gen-es2Life, Irapuato, Mexico) was used, which includes reverse transcription and PCR detection of three targets of N ORF (N1-FAM, N2-HEX and N3-TexasRed) in a single multiplex reaction. For each reaction, 5 µl of viral RNA were added as template in a total volume of 25 µl. The RT-qPCR reaction was performed in a QuantStudio 5 (Thermo Scientific) instrument, and settings for passive reference was changed to 

The transport model for dispersion of virus in unsaturated media (Torkzaban et al. 2006 ) was adapted. The original model describes the viral concentration with respect to the distance from the virus input site and the decrease over time, and considers terms of adsorption and detachment of virus with respect to the soil-water interface and air-water interface. We considered only dispersion in the water bulk; we also considered that dispersion of virus is independent on time, and it depends on the distance. Hence, the attachment and detachment rate coefficients and the inactivation rate coefficient in soil and air media were equal to zero (κ a att = 0, κ a det = 0, κ s att = 0, κ s det = 0, μ s = 0 and μ a = 0). We also assume that Darcy flux (q) and inactivation rate coefficient (μ w ) are constant along the river. The variables amount of water (θ) and coefficient of dispersion (D) were fitted to the experimental data. Finally, we excluded the contributions of the adsorption from the general model and considered two entrances corresponding to the Tepeji River at 5.88 km and the Salado River at 57.88 km (with respect to RW2; Figs. 1c and 2a) , both of which affect the viral concentration (Supplementary File 1) . Detailed derivation of the equations is shown in Supplementary File 1. The resulting Eq. (1) is shown below:

where Cw, is the viral concentration and the other variables are previously decribed.

The equation was solved numerically with initial and boundary conditions of initial distance of zero at the point RW2 corresponding to the Tepeji River (as this is the initial entrance of water to the Tula River) and an initial normalized viral concentration at the point RW2, according to the maximum value obtained by the experimental measurements. We used the software package Scipy Odeint in Python. The concentration values for the virus in the Tula River were obtained experimentally as described above. 

We detected RNA from SARS-CoV-2 in water samples collected in the Tula River, Tepeji River and Salado River (Fig. 2) . Detection was achieved from the three targets included in the kit (N1-FAM, N2-HEX and N3-TexasRed) in most cases; however, quantification using N2-HEX Ct values and its corresponding standard curve yielded more consistent results, and absolute values were always intermediate to those obtained with N1-FAM and N3-TexasRed (not shown). From these criteria, in this study we used N2-HEX data; for this probe and following the protocols detailed in methods, we calculated a limit of detection of <1 copy per ml of water sample. SARS-CoV-2 RNA was detected in six out of eight water samples from river ( Fig. 2) . As expected, the highest concentration was found in sample RW2, at the entrance flow of the wastewater treatment plant (WWTP, Figs. 1c and 2), with a concentration of 79 RNA copies per ml; from the sampled locations in the Tula River, this is the closest to the Mexican Valley Metropolitan Area (MVMA; Figs. 1 and 2) . This value decreased through the flow of the river (which consists mostly of untreated wastewater) in samples RW3 and RW4, and was undetectable in RW5, downstream of the Endhó Dam. SARS-CoV-2 was also detected at the entrance fluxes of Tepeji River (RW1, 18 copies per ml) and Salado River (RW6, 24 copies per ml), both of which contribute to viral load in the Tula River. Finally, SARS-CoV-2 was not detected in sample RW8, which was collected at Ixmiquilpan, ≈40 km downstream of sample RW7 (Fig. 2) .

RNA from SARS-CoV-2 was also detected in water samples collected in irrigation canals throughout the Mezquital Valley. Irrigation canals are mostly fed from treated water, however, untreated wastewaters from municipalities in Hidalgo are also connected to these canals. Since the irrigation canal network in the Mezquital Valley is highly complex and a detailed map is unavailable, we only show the general known source of water in each location (Fig. 3a) . RNA From SARS-CoV-2 was detected in 5 out of 8 water samples collected from irrigation canals (Fig. 3b) . The highest concentration of 112 copies per ml was found in sample CW1 (Fig. 3b) . This location is the closest from MVMA, however, it is reportedly fed solely from wastewater treated in the WWTP (Fig. 3a, Table 1 ). SARS-CoV-2 RNA was also detected in samples CW2, CW4, CW4-II and CW5, all of which are fed from the WWTP, but also from wastewater from municipalities in Hidalgo, and from the Endhó Dam (Fig. 3) . SARS-CoV-2 was not detected in samples CW3, CW6 and CW7; notably, irrigation canals at Ixmiquilpan (CW6 and CW7) are fed from freshwater bodies (not wastewater) such as groundwater wells or local spas. In contrast, CW3 is reportedly fed from the WWTP and local municipalities (Fig. 3) and lack of detection in this specific location may be related to degradation or dilution along the flow through the complex canal system, where distance from urban sources is mostly unknown.

Since SARS-CoV-2 RNA was detected in water samples from irrigation canals, we hypothesized that it may be detectable in agricultural soil from adjacent fields (Fig. 1c , Table 1 ). However, SARS-CoV-2 was not detected in these samples (Supplementary Table S1 and S2). Similarly, Fig. 3 Quantification of SARS-CoV-2 ARN in water from irrigation canals in the Mezquital Valley. a Schematic representation showing the reported source of water in irrigation canals sampled. b RT-qPCR results of SARS-CoV-2 quantification in water samples from irrigation canals SARS-CoV-2 RNA was not detected in any of the produce samples, which were collected in the field or at the Ixmiquilpan Market (Table 1, Supplementary Table S2 ). Since the internal control PMMV was detected in all of these samples (Supplementary Table S2 ), we discard the presence of PCR inhibitors as the reason for lack of detection, but further studies are needed to confirm the absence of SARS-CoV-2 since viral particles can adsorb to soil particles and other solids.

Results on physicochemical analyses are shown in Supplementary Table S3 and results on microbiological analyses are shown on Supplementary Table S4 . Based on standards established by Mexican laws, all water samples had values of BOD, TSS and SS below the permitted limits for water used for irrigation (NOM-001-SEMARNAT-1996) . Regarding limits for microbiological parameters (NOM-001-ECOL-1996) , five out of eight water samples from the Tula River, and 5 out of 8 water samples from irrigation canals, presented total coliform values above the allowed limits for wastewater used for irrigation (Supplementary Table S4 ).

We assessed if the physicochemical characteristics of wastewater, organic matter and microbial concentration correlated with the concentration of the virus. In samples from the river water, correlation analyses showed that several variables were associated with SARS-CoV-2 RNA concentration (Table 2 , left). In contrast, analyses using parameters obtained from irrigation canals, showed no significant correlations between SARS-CoV-2 RNA concentration and physicochemical or microbiological variables (Table 2, right). Significant p values were obtained for correlation analyses with physicochemical parameters BOD, COD, pH, and SS (Table 2 , left), and also for all three microbiological parameters in river samples (Table 2 , Supplementary Fig. S1 ). Further inspection of correlation trends revealed correlations with pH and SS could be spurious, since they were highly grouped in discrete values ( Supplementary Fig. S1 ), and a higher number of samples should be analyzed to confirm these observations (Table 1 , Suplementary Fig. S1 ).

COD and BOD are measurements of oxygen required for chemical oxidation of organic material in water. From these values, the Biodegradability index can be calculated (COD/ BOD ratio) as an indicator of the non-biodegradable organic-chemicals. We assessed if the Biodegradability index correlated with the virus concentration of the water system (Abdalla and Hammam 2014); interestingly, we found a clear correlation between Biodegradability index and SARS-CoV-2 concentration in samples from the Tula River ( Supplementary Fig. S2 ).

Using data for SARS-CoV-2 concentration in the Tula River, we sought to adapt a model that describes viral dispersion in the wastewater bulk. Our model only considers the dispersion of the virus in water and the modification of the flux in the mainstream. The dispersion coefficient was calculated through least square optimization to the experimental data (Suplementary File 1) (Molugaram et al. 2017) . The concentration in the initial sampling site RW2 (Fig. 1c ) was normalized to one (Ratio of Concentration C/initial Concentration Ci); since it was the maximum value, the experimental data were normalized to this value. Additionally, we modeled two entrances at the points RW2 and RW6, as a result the dispersion coefficient increases at these points. The model fitted to the experimental data along the Tula River (Fig. 4a) . Therefore, with only three variables, dispersion coefficient (D), amount of water (Θ) and a distance function of D [D(x)], we could define the dispersion of the virus in the Tula River. Finally, we used QGIS3 to obtain a geographic representation of our model and the experimental data (Fig. 4b) , which show the general trends of SARS-CoV-2 dispersion in the Tula River.

In this work, we quantified SARS-CoV-2 RNA in wastewaters from Mexico City, which are distributed throughout the Tula River, Salado River and irrigation canals and are used for irrigation of crops in the Mezquital Valley. Notably, we detected the virus in geographically representative water samples from both rivers and the irrigation canals. With our data, we generated a viral dispersion model which can be used to assess concentration of SARS-CoV-2 to other river systems if the dispersion coefficient (D), amount of water (Θ) and a distance function of D [D(x)] of the river are fitted. Samples from both the Tula River and irrigation canals followed the same general trend of presenting the highest concentration in the sites located near the WWTP. Viral concentration decreased in the remaining sites and was undetectable in the northernmost sites, located near the town of Ixmiquilpan. This result indicates that the most important source of SARS-CoV-2 RNA in these water systems is the wastewater from Mexican City Valley, and that the viral particles tend to degrade or dilute along the flow of water through the Mezquital Valley. However, we do not discard that municipalities in the Mezquital Valley contribute to viral load and fecal contamination of Tula River, Salado River and the irrigation canals.

Our results showed that the SARS-CoV-2 virus was not detected on soil samples. While this could be related to PCR inhibitors in the samples, we found that the PMMV virus was detected and hence this hypothesis was discarded. Virus associate with soil and other particulate matter and both adsorption and detachment of virus depends on temperature, moisture, pH, and the physicochemical characteristics of the virus capsid surface and the particles. Similarly, the microbes and organic matter may have a protective effect on virus (AzadpourKeeley et al. 2003) . Further studies are needed to asses if lack of detection of SARS-CoV-2 is related to degradation or dispersion of the virus, or alternatively, if viral particles adsorb in the soil matrix, which would make our RNA extraction method not appropriate.

We found that physicochemical variables COD and BOD and microbiological variables showed significant correlations with SARS-CoV-2 quantification, but only in samples from the Tula River. However, these observations should be confirmed with measurements in a higher number of samples. This result suggests that in the river, anthropological sources from a single source of contamination predominate, i.e. the MVMA. However, irrigation canalswhere no correlations were foundmay receive other sources of contamination (from agriculture, livestock, or industry) (Guédron et al. 2014; Lesser et al. 2018) . Furthermore, our correlation analyses suggest that BOD, COD, and the Biodegradability Index (BOD/COD), as well as fecal coliforms may be good indicators of SARS-CoV-2 RNA concentration in municipal wastewaters. Fig. 4 The transport model generated for SARS-CoV-2 dispersion in the Tula River, fits the data obtained experimentally. a Graphic representation of viral concentration in the Tula River, comparing the experimental data and the transport model obtained in this study. b Map representation of the geographic dispersion of SARS-COV-2 presence in the Tula River, comparing the transport model (left) and the experimental data (right)

The model proposed to explain the dispersion of SARS-CoV-2 virus along the Tula River, as well as the experimental data, indicate that the viral concentration decreases from south to north, with an initial dispersion point at the entrance of the Tula River at the WWTP and a minimal contribution from the communities along the Mezquital Valley (Fig. 4a, b) . Our model predicts viral dispersion in the Mezquital Valley, and may be useful for the selection of sampling locations prior to temporal monitoring to quantify the evolution of SARS-CoV-2 (or other virus) epidemics in populations. Previously available models for viral dispersion in environmental matrices (Bivins et al. 2020; Farkas et al. 2020; Kitajima et al. 2020) are useful for detailed dissection of variables governing dispersion patterns of viral particles in at a smaller scale in controlled environments. Our model is, to our knowledge, the first to allow prediction of viral dispersion in linear water bodies at a geographic scale. This model could be extrapolated to other wastewater streams, rivers, canals, etc.

One important purpose of this work was to assess the possible risk of environmental transmission of SARS-CoV-2 through environmental matrices. Indeed, some recent reports suggest that the virus may be infective (Kitajima et al. 2020) , and hence, direct contact with water bodies in the Mezquital Valley may represent a risk of transmission, especially for farmers during flooding irrigation (Lüneberg et al. 2018) . While our results suggest that improved wastewater treatment processes are needed as the WWTP normally works at 30% of its capacity (Chamizo-Checa et al. 2020) , it is important to note that the producers refuse to use the treated water in their irrigation systems. Indeed, wastewaters contribute with large amounts of organic matter decreasing the costs of fertilization. In turn, increased wastewater treatment could directly elevate production costs, price of water, and needs for fertilizers, which could negatively impact land rentability (Pérez Camarillo 2002; Pérez et al. 2006; SIAP 2018) .

This work does not seek to eliminate wastewater usage in the Mezquital Valley, as this practice is one of the key factors that allowed the economic development of the region. In fact, the Mezquital Valley is considered a study model for the agricultural development of arid and semiarid rural regions (Anderson 2020; Bonvehi Rosich and Seth Denizen 2021; Seth Denizen 2021) . We suggest that our work should be consideredalong with social, economic and agronomic factors--for the creation and modification of practices and policies related to water treatment (Belhadi et al. 2020) , pathogen monitoring, agricultural practices, and safety measures for farmers in the Mezquital Valley. Finally, our results support the widespread notion that wastewater surveillance is an efficient method to follow transmission of SARS-CoV-2 in populations, especially in low-income regions where massive testing is not economically viable. Furthermore, the information generated through these strategies could also feed the decision-making processes during vaccination programs, by focusing efforts in populations with high-incidence.

Correlation between biochemical oxygen demand and chemical oxygen demand for various wastewater treatment plants in Egypt to obtain the biodegradability indices

First confirmed detection of SARS-CoV-2 in untreated wastewater in Australia: a proof of concept for the wastewater surveillance of COVID-19 in the community

The Right to Sewage: Agriculture, climate change, and the growing need for cities to embrace wastewater reuse, Harvard Graduate School of Design

Infectious waste management strategy during COVID-19 Pandemic in Africa: an integrated decision-making framework for selecting sustainable technologie

Movement and Longevity of Viruses in the Subsurface. Environmental Protection Agency, EPA

Wastewater-based epidemiology: global collaborative to maximize contributions in the fight against COVID-19'

New Visions for Wastewater Equity in The Mezquital Valley

Megacity wastewater poured into a nearby basin: looking for sustainable scenarios in a case study

Environmental Management

The presence of SARS-CoV-2 RNA in the feces of COVID-19 patients

Health risks from exposure to untreated wastewater used for irrigation in the Mezquital Valley, Mexico: A 25-year update

Wastewater and public health: the potential of wastewater surveillance for monitoring COVID-19

El agua residual como generadora del espacio de la actividad agrícola en el Valle del Mezquital

Methyl)mercury, arsenic, and lead contamination of the world's largest wastewater irrigation system: the Mezquital Valley (Hidalgo State-Mexico)

First environmental surveillance for the presence of SARS-CoV-2 RNA in wastewater and river water in Japan

Wastewater surveillance for SARS-CoV-2 on college campuses: initial efforts, lessons learned, and research needs

Simulated 2017 nationwide sampling at 13,940 major U.S. sewage treatment plants to assess seasonal population bias in wastewater-based epidemiology

Modeling wastewater temperature and attenuation of sewage-borne biomarkers globally

First case of 2019 novel coronavirus in the United States

Highly pathogenic influenza A(H5N1) virus survival in complex artificial aquatic biotopes

Catálogo de las lenguas indígenas nacionales: variantes lingüísticas de México con sus autodenominaciones y referencias geoestadísticas. 1

SARS-CoV-2 in wastewater: State of the knowledge and research needs

SARS-CoV-2 Detection in Istanbul Wastewater Treatment Plant Sludges, medRxiv, 2020.05.12

First proof of the capability of wastewater surveillance for COVID-19 in India through detection of genetic material of SARS-CoV-2

The role of wastewater treatment plants as tools for SARS-CoV-2 early detection and removal

Survey of 218 organic contaminants in groundwater derived from the world's largest untreated wastewater irrigation system

SARS-CoV-2 in wastewater: potential health risk, but also data source

Drylands soil bacterial community is affected by land use change and different irrigation practices in the Mezquital Valley

Presence of SARS-Coronavirus-2 RNA in sewage and correlation with reported COVID-19 prevalence in the early stage of the epidemic in the Netherlands

Statistical Techniques for Transportation Engineering

Comisión Nacional para el Desarrollo de los Pueblos Indígenas (Pueblos indígenas del México contemporáneo)

Statement on the second meeting of the International Health Regulations (2005) Emergency Committee regarding the outbreak of novel coronavirus (2019-nCoV)

Water, sanitation, hygiene, and waste management for the COVID-19 virus: interim guidance

Evolución del rendimiento y de la producción de maíz de riego en el Valle del Mezquital

Carbón de la Espiga de Maíz. Descripción, síntomas y prevención

SARS-CoV-2 RNA in wastewater anticipated COVID-19 occurrence in a low prevalence area

En este pueblo mexicano todo lo que corre, se arrastra o vuela va a la cazuela

Resultados del Estudio del Diagnóstico Sectorial

Thinking Through Soil: Case Study from the Mezquital Valley', Harvard Graduate School of Design

Atlas Agroalimentario

An updated estimation of the risk of transmission of the novel coronavirus (2019-nCov)

Role of air-water interfaces on retention of viruses under unsaturated conditions

Predicting the number of people infected with SARS-COV-2 in a population using statistical models based on wastewater viral load, medRxiv

High-throughput wastewater analysis for substance use assessment in central New York during the COVID-19 pandemic

Detection of SARS-CoV-2 in Different Types of Clinical Specimens

SARS-CoV-2 titers in wastewater foreshadow dynamics and clinical presentation of new COVID-19 cases, medRxiv, 2020.06.15

SARS-CoV-2 titers in wastewater are higher than expected from clinically confirmed cases

WebGestalt: an integrated system for exploring gene sets in various biological contexts

Epidemiological, clinical characteristics of cases of SARS-CoV-2 infection with abnormal imaging findings

Acknowledgements This work was supported by CONACyT, Mexico, grant 312014 for JR. This work was possible thanks to the collaboration of local farmers in the state of Hidalgo, who allowed us to collect samples in their fields. The authors thank Gabriela Rivas (UPFIM), Ana Alcalá (IBT UNAM), Laura Palomares (IBT), Oscar Monroy (UAMI) and Jesús Hernández (CIAD Hermosillo) for advice and technical assistance. Eneas Aguirre guided the sample collections at milpas from San Agustin Tlaxiaca. We thank personnel from governmental institutions in the state of Hidalgo: Jaime Ortega (SEDA-GROH), Sergio Guzman (CONAGUA) and José Alonso Huerta (CITNOVA), for their valuable contribution in the stages of planning, sampling, and discussion of results.

Conflict of Interest The authors declare no competing interests.Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.