The groundwater samples were clustered into four groups. The clustering of the samples led to the finding that streamflows play a significant role in the hydrological balance as a source of local recharge to the aquifer.

In this research, vegetation indices (VIs) were analyzed as indicators of the spatio-temporal variation of vegetation in a semi-arid region. For a better understanding of this dynamic, interactions between vegetation and climate should be studied more widely. To this end, the following methodology was proposed: (1) acquire the NDVI, EVI, SAVI, MSAVI, and NDMI by classification of vegetation and land cover categories in a monthly period from 2014 to 2018; (2) perform a geostatistical analysis of rainfall and temperature; and (3) assess the application of ordinary and uncertainty least squares linear regression models to experimental data from the response of vegetation indices to climatic variables through the BiDASys (bivariate data analysis system) program. The proposed methodology was tested in a semi-arid region of Zacatecas, Mexico. It was found that besides the high values in the indices that indicate good health, the climatic variables that have an impact on the study area should be considered given the close relationship with the vegetation. A better correlation of the NDMI and EVI with rainfall and temperature was found, and similarly, the relationship between VIs and climatic factors showed a general time lag effect. This methodology can be considered in management and conservation plans of natural ecosystems, in the context of climate change and sustainable development policies.

The presence of arsenic in groundwater constitutes a hazard for the environment and human health, and the determination of its source has become a global challenge, which can be approached by defining the natural background levels (NBL) in conjunction with the indicator kriging method, with the aim of delineating anthropogenically contaminated areas. However, having a unique value of NBL for large areas can generate interpretation errors. This research integrates the determination of the flow systems present in the Calera Aquifer, and the definition of the natural background levels in each flow system by making estimation maps in ArcGIS using two databases, 10 years apart, to evaluate the spatio-temporal variation of arsenic in groundwater. The results indicate a notable increase in the probability of exceeding the arsenic NBL, mainly in the intermediate flow, which may be due to movement resulting from mining activities as well as a mixture of regional and intermediate flows caused by the extraction of water for agriculture and drinking water supplies. The presented values exceed the maximum limits allowed for human consumption, as stated by the World Health Organization.

The estimation of the hydraulic parameters of an aquifer such as the hydraulic conductivity is somehow complicated due to its heterogeneity, on the other hand field and laboratory tests are both time consuming and costly. The use of geostatistical-based techniques for data assimilation could represent an alternative tool that allows the use of space-time aquifer behaviour to characterize hydraulic conductivity heterogeneity. In this paper, a spatiotemporal bivariate methodology was implemented combining historical hydraulic head data with hydraulic conductivity sparse data in order to obtain an estimate of the spatial distribution of the latter variable. This approach takes advantage of the correlation between the hydraulic conductivity (K) and the hydraulic head (H) behaviour through time. In order to evaluate this approach, a synthetic experiment was constructed through a transitory numerical flow-model that simulates hydraulic head values in a horizontally-heterogeneous aquifer. Geostatistical tools were used to describe the correlation between simulated spatiotemporal data of hydraulic head and the spatial distribution of the hydraulic conductivity in a group of model nodes. Subsequently, the Kalman filter was used to estimate the hydraulic conductivity values at nonsampled sites. The results showed acceptable differences between estimated and synthetic hydraulic conductivity data, with low estimate error variances (predominating the 1 m2/day2 value for K for all the cases, however, the smallest number of cells with values above 2 m2/day2 correspond to the bivariate spatiotemporal case) and the best agreement between the estimated errors and the selected model variance (SMSE values of 0.574 and 0.469) were found for the bivariate cases, which suggests that the implemented methodology could be used for reducing calibration efforts, particularly when the hydraulic parameters data are scarce.

The supply of drinking water to the population is an important challenge facing humanity, since both surface and underground sources present a great variability of water storage with respect to space and time. This problem is further aggravated in arid and semi-arid areas where rainfall is low and torrential, which makes groundwater the main source of supply; therefore, it is necessary to carry out studies that allow evaluating the evolution of the quantity and quality of water. This study addresses the behavior of groundwater in a semi-arid region, considering the theory of flow systems to identify movement as well as water quality, es determined by a water quality index (WQI), calculated considering arsenic and fluorine. In addition, a quality irrigation classification is used, employing the norms in accordance with international standards and the Mexican Norm, which allows for a comparison. Local, regional, intermediate and mixed flow systems are identified, and the evolution of cations and anions in addition to temperature is examined. It is observed that the drinking water quality index classifies them as excellent in most of the monitored wells (<50), but with a negative evolution. Regarding irrigation, most of the water samples are classified without restriction for the establishment of any type II crop (C2S1) and with restrictions for horticultural crops. It is observed that arsenic had values between 0.49 and 61.40 (µg/L) in 2005, while in 2015 they were between 0.10 and 241.30 (µg/L). In addition, fluoride presented values between 0.00 and 2.6 (mg/L) in 2005, while in 2015 they were between 0.28 and 5.40 (mg/L). The correlations between arsenic and fluorine are noted as well as WQI and SAR. A finding in this research was to include arsenic and fluorine in the calculation of the WQI allowing a better interpretation of the quality of water for both human consumption and for agricultural use to based on this make the best decision to control any harmful effects for the population, in addition to identifying the appropriate purification treatment required to control pollutants. It is concluded that arsenic is an element of utmost importance when considering water quality, so it is necessary to examine its evolution and continue to monitor its levels constantly.

One-third of the global population depends on groundwater for drinking, which is an even larger proportion for arid regions. The integration of isotopic and geophysical applications has been very useful in understanding the process of groundwater recharge. The aim of this study is to define a conceptual model that describes groundwater functions within an aquifer located in a semi-arid region by identifying recharge patterns based on the isotopic characteristics of: Rainfall, surface water, shallow and deep groundwater, and incorporating regional geophysical data. We demonstrated that rainfall was affected by sub-cloud evaporation and altitude. Shallow and deep modern groundwater samples were clustered and exhibited similar evolution from rainfall. However, different groups recharged from different precipitation sources compared to the local one. In the current study, we analyzed the isotopic evolution of deep groundwater over a 10-year period, which was mainly affected by the incorporation of different flows with different isotopic signatures and the hydrodynamics of the area. We performed two geoelectrical sections in the study area to improve the understanding of the hydrogeological setting and water movement patterns. The new conceptual model should help stakeholders in the context of water management policies for the study area.

The increasing population demands a greater quantity of food. In order to satisfy the world’s demand, one of the main challenges worldwide consists of modernizing the current irrigation systems. This research shows an experience carried out in a modernized irrigation module in Central-North Mexico following these objectives: to evaluate the impact of the modernization of the irrigation module, to analyze the agricultural productivity, and to assess a group of parameters related with the agricultural production (system conduction and distribution efficiencies, water productivity, among others) and the water volumes after and before the modernization. After a drought period, a methodology was performed in commercial parcels in 2013 to increase the yield of different crops. Some of the activities were: soil leveling, estimation of the soil properties (field capacity, wilting point, bulk density, pH, and organic matter), optimum fertilization applications, use of a model to scheduling irrigation, measure volumes extracted at the parcel level. With the modernization and the method used, around 1800 and 2000 m³ ha–1 were saved with respect to the initial granted volume by the Comisión Nacional del Agua and increase in the global efficiency was also achieved (from 55% to 85%). All crop yields increased, i.e., for corn from 2.5 kg/m³ to 3.8 kg/m³. The impact of modernization accompanied with an effective operation allowed a significant increase of the crop yield and water productivity. Despite a controlled distribution of water being carried out, future research should contemplate free water demand scenarios and automation irrigation for improving the module operation.

Infiltration estimation is made by tests such as concentric cylinders, which are prone to errors, such as the lateral movement under the ring. Several possibilities have been developed over the last decades to compensate these errors, which are based on physical, electronic, and mathematical principles. In this research, two approaches are proposed to measure the water infiltration rate in a silty loam soil by means of the mass values of a lysimeter weighing under rainfall conditions and different moisture contents. Based on the fact that with the lysimeter it is possible to determine acting soil flows very precisely, then with the help of mass conservation and assuming a downward vertical movement, 12 rain events were analyzed. In addition, it was possible to monitor the behavior of soil moisture and to establish the content at field capacity from the values of the weighing lysimeter, from which both approach are based. The infiltration rate of these events showed a variable rate at the beginning of the rainfall until reaching a maximum, to descend to a stable or basic rate. This basic infiltration rate was 1.49 ± 0.36 mm/h, and this is because soils with fine textures have reported low infiltration capacity. Four empirical or semi-empirical models of infiltration were calibrated with the values obtained with our approaches, showing a better fit with the Horton’s model.