Integrated water resource management (IWRM) is facing great challenges due to growing uncertainties caused by climate change (CC), rapid socio-economic and technological changes, and population growth. In the present study, we have developed different indices to assess the availability of water using an IWRM approach. These indices evaluate supply to demands, surface availability, groundwater availability, reservoirs, and environmental flow. Moreover, reliability, resilience, and vulnerability were determined. Sustainability index (SI) and sustainability index by groups (SG) were determined based on the five indices (all indices vary from 0 to 1). The impacts of climate change affect surface and groundwater availability, as do the agricultural, urban, and industrial requirements on the different supplies. We used the generalized AQUATOOL Decision Support System Shell (DSSS) to evaluate the IWRM in the Rio Grande Basin (Morelia, México). Various emission scenarios from representative concentration pathways (RCPs) were applied to the basin for the years 2015–2039 and 2075–2099. The results indicate increases in agricultural and urban demand, and decreases in surface runoff, as well as groundwater recharge. The proposed indices are useful for different approaches (decision-makers, water policy, and drought risks, among others). CC significantly affects the different proposed indices and indicates a decrease of the SI, SG1, and SG2 (i.e., less availability). For example, we found that SG2 decreased from 0.812 to 0.195 under the RCP 8.5 2075–2099 scenario, and SG2 equal to 0.252 and 0.326 for the RCP 6.0 2075–2099 and RCP 4.5 2070–2099 scenarios, respectively (values close to 0 indicate worst drought conditions).

In recent years, due to various anthropogenic activities, such as agriculture and livestock, the presence of nitrogen-associated contaminants has been increasing in surface- and groundwater resources. Among these, the main compounds present in groundwater are ammonia, nitrite, and nitrate. However, it is sometimes difficult to assess such effects given the scarcity or lack of information and the complexity of the system. In the current study, a methodology is proposed to assess nitrate in groundwater from diffuse sources considering spatiotemporal patterns of hydrological systems using a coupled SWAT/MODFLOW/MT3DMS model. The application of the model is carried out using a simplified simulation scheme of hydrological and agricultural systems because of the limited spatial and temporal data. The study area includes the Cuitzeo Lake basin in superficial flow form and the Morelia–Querendaro aquifer in groundwater flow form. The results within the methodology are surface runoff, groundwater levels, and nitrate concentrations present in surface- and groundwater systems. The results indicate that the historical and simulated nitrate concentrations were obtained within acceptable values of the statistical parameters and, therefore, are considered adequate.

Calibration of conceptual rainfall–runoff models (CRRM) for water-resource assessment (WRA) is a complicated task that contributes to the reliability of results obtained from catchments. In recent decades, the application of automatic calibration techniques has been frequently used because of the increasing complexity of models and the considerable time savings gained at this phase. In this work, the traditional Rosenbrock (RNB) algorithm is combined with a random sampling method and the Latin hypercube (LH) to optimize a multi-start strategy and test the efficiency in the calibration of CRRMs. Three models (the French rural-engineering-with-four-daily-parameters (GR4J) model, the Swedish Hydrological Office Water-balance Department (HBV) model and the Sacramento Soil Moisture Accounting (SAC-SMA) model) are selected for WRA at nine headwaters in Spain in zones prone to long and severe droughts. To assess the results, the University of Arizona’s shuffled complex evolution (SCE-UA) algorithm was selected as a benchmark, because, until now, it has been one of the most robust techniques used to solve calibration problems with rainfall–runoff models. This comparison shows that the traditional algorithm can find optimal solutions at least as good as the SCE-UA algorithm. In fact, with the calibration of the SAC-SMA model, the results are significantly different: The RNB algorithm found better solutions than the SCE-UA for all basins. Finally, the combination created between the LH and RNB methods is detailed thoroughly, and a sensitivity analysis of its parameters is used to define the set of optimal values for its efficient performance.

Groundwater models serve the function of predicting and analyzing aquifer behavior. They require input information, such as hydrogeological parameters like hydraulic conductivity and storage coefficient, which are used to calibrate the model, and elementary actions that include recharge and extracted volumes. There are cases in which it is insufficient to know the homogeneous recharge entering through the surface basin, referred to as traditional recharge, since, in many instances, the distribution is altered by changes in land use. For this reason, based on the geomorphological characteristics of the basin, weighting is proposed for sites with greater recharge capacity. The present work shows a solution to the recharge distribution using the potential groundwater recharge (PGR) map, which is formed by weighting spatially distributed information: (i) drainage, (ii) precipitation, (iii) land use, (iv) geological faults, (v) soil type, (vi) slope, and (vii) hydrogeology. A comparison is made between groundwater modeling using traditional recharge and PGR recharge. It is noted that the modeling perform similarly for both recharges, and the errors do not exceed 5% absolute error, which validates the model’s reliability. This manuscript demonstrates how to model and calibrate groundwater in aquifers with scarce information and variable recharge, making it reproducible.

Temperature is one of the most influential weather variables necessary for numerous studies, such as climate change, integrated water resources management, and water scarcity, among others. The temperature and precipitation are relevant in river basins because they may be particularly affected by modifications in the variability, for example, due to climate change. We developed a stochastic model for daily precipitation occurrences and their influence on maximum and minimum temperatures with a straightforward approach. The Markov model has been used to determine everyday occurrences of rainfall. Moreover, we developed a multisite multivariate autoregressive model to represent the short-term memory of daily temperature, called MASCV. The reduction of parameters is an essential factor addressed in this approach. For this reason, the normalization of the temperatures was performed through different nonparametric transformations. The case study is the Jucar River Basin in Spain. The multisite multivariate stochastic model of two states and a lag-one accurately represents both occurrences as well as maximum and minimum temperature. The simulation and generation of occurrences and temperature is considered a continuous multivariate stochastic process. Additionally, time series of multiple correlated climate variables are completed. Therefore, we simplify the complexity and reduce the computational time for the simulation.