Rainwater harvesting (RWH) has been recognized as one of the most reliable and efficient methods for water supply, especially in arid and semi-arid regions (ASARs) facing freshwater scarcity. Nevertheless, due to the inherent uncertainty of input data and subjectivity involved in the selection of influential parameters, the identification of RWH potential areas is a challenging procedure. In this study, two approaches for locating potential RWH sites were implemented. In the first approach, a frequently-used method of the multi-criteria decision analysis and geographic information system (MCDA-GIS) was utilized, while, in the second approach, a novel strategy of integrating the soil and water assessment tool (SWAT) model as a hydrology model into an MCDA-GIS method was proposed to evaluate its performance in locating potential RWH sites. The Mashhad Plain Basin (MPB) was selected as a case study area. The developed potential RWH maps of the two approaches indicated similar patterns for potential RWH areas; in addition, the correlation coefficient (CC) between the two obtained maps were relatively high (i.e., CC = 0.914) revealing that integration of SWAT as a comprehensive hydrologic model does not necessarily result in very different outputs from the conventional method of MCDA-GIS for RWH evaluation. The overlap of developed maps of the two approaches indicated that 3394 km² of the study area, mainly located in the northern parts, was identified as high-potential RWH areas. The performed sensitivity analysis indicated that rainfall and slope criteria, with weights of 0.329 and 0.243, respectively, had the greatest sensitivity on the model in the first approach while in the second approach, the criterion of runoff coefficient (with weights of 0.358) had the highest impact. Based on results from the identification of the potential locations for conventional RWH techniques, pond and pan techniques are the most proper options, covering high-potential areas of RWH more effectively than other techniques over MPB.

The intensity–duration–frequency (IDF) curve is a commonly utilized tool for estimating extreme rainfall events that are used for many purposes including flood analysis. Extreme rainfall events are expected to become more intense under the changing climate, and there is a need to account for non-stationarity IDF curves to mitigate an underestimation of the risks associated with extreme rainfall events. Sydney, Australia, has recently started experiencing flooding under climate change and more intense rainfall events. This study evaluated the impact of climate change on altering the precipitation frequency estimates (PFs) used in generating IDF curves at Sydney Airport. Seven general circulation models (GCMs) were obtained, and the best models in terms of providing the extreme series were selected. The ensemble of the best models was used for comparing the projected 24 h PFs in 2031–2060 with historical values provided by Australian Rainfall and Runoff (ARR). The historical PFs consistently underestimate the projected 24 h PFs for all return periods. The projected 24 h 100 yr rainfall events are increased by 9% to 41% for the least and worst-case scenario compared to ARR historical PFs. These findings highlight the need for incorporating the impact of climate change on PFs and IDF curves in Sydney toward building a more prepared and resilient community. The findings of this study can also aid other communities in adapting the same framework for developing more robust and adaptive approaches to reducing extreme rainfall events’ repercussions under changing climates.

Understanding and projecting drought, especially in the face of climate change, is crucial for assessing its impending risks. However, the causes of drought are multifaceted. As the environmental research paradigm pivots towards machine learning (ML) for predictions, our investigation contrasted multiple ML techniques to simulate the Standardized Precipitation Evapotranspiration Index (SPEI) from 2009 to 2022, utilizing various potential evapotranspiration (PET) methods. Our primary focus was Australia, the world’s driest inhabited continent. Given the challenges with ML model interpretation, SHAP (SHapley Additive exPlanations) values were employed to decipher SPEI variations and to gauge the relative importance of precipitation (Prec) and PET in six key Australian cities. Our findings revealed that while different PET methods resulted in distinct mean values, their trends remained consistent. Post the Millennium Drought, Australia experienced several drought events. SPEI discrepancies based on PET methods were minimal in humid regions like Brisbane and Darwin. However, for arid cities, the Priestley–Taylor equation-driven SPEI differed notably from other methods. Ridge regression was the most adept at mirroring SPEI changes among the assessed ML models. Furthermore, the SHAP explainer discerned that PET-related climate variables had a greater impact on SPEI in drier cities, whereas in humid cities, Prec was more influential. Notably, the research emphasised CO₂′s role in influencing drought dynamics in humid cities. These insights are invaluable for enhancing drought mitigation strategies and refining predictive models. Such revelations are crucial for stakeholders aiming to improve drought prediction and management, especially in drought-prone regions like Australia.

Surface Water–Groundwater (SW–GW) interaction is a crucial aspect of the hydrological cycle and requires accurate modeling for reliable predictions. In many basic hydrological models and calculations, it is common to assume that the water density is constant. However, density-dependent flow, which accounts for changes in water density, plays a significant role in various hydrological processes. This study aims to quantify the effects of density-dependent flow on SW–GW interaction and evaluate the sensitivity of dominant hydrological drivers to density-dependent flow. Our simulations using the HydroGeoSphere model revealed that neglecting density-dependent flow in SW–GW interaction can lead to inaccurate estimations of water and solute balances. In particular, including density-dependent flow in the model yielded more realistic salinity distributions under gaining river scenarios and captured the gradual expansion of freshwater lenses under losing river scenarios. The results also indicated that under non-density-dependent flow, more saline groundwater is exposed to evapotranspiration, resulting in higher solute mass storage and a more saline unsaturated zone. Further, surface recharge and pumping rates played crucial roles in salinity distribution. This study highlights the critical importance of incorporating density-dependent flow in simulations, providing valuable insights for improving the accuracy of predictions and effectively managing water and solute balances in floodplain aquifers.