Hydrologic Characterization Research Articles

Modeling and Optimization of Detention Basins For Cost Minimization and Flood Mitigation in Urban Areas

Objective: The objective of this study was to develop and apply an optimization model for the sizing of detention reservoirs in an urban watershed, aiming to minimize implementation costs and mitigate flood risks in critical areas. Theoretical Framework: The study is based on flood control models and urban drainage infrastructure, focusing on optimization techniques applied to water resource management. References include the SCS Curve Number method for surface runoff, the PCSWMM hydrological model for precipitation simulations, and the HEC-RAS hydraulic model for flood area delineation. Method: The research was divided into three main stages: (i) physical and hydrological characterization of the watershed using PCSWMM; (ii) delineation of flood-prone areas with HEC-RAS 2D; and (iii) development of a nonlinear optimization model in GAMS for the sizing of detention reservoirs. The decision variables included the height, surface area, and outlet diameter of the reservoirs, with implementation cost as the main objective function. Results and Discussion: The optimization modeling determined the ideal characteristics of six reservoirs distributed throughout the watershed to maximize retention and flow control. The total flooded area was estimated at 270,989 m², with maximum depths reaching 1.914 meters. The total cost of the reservoir system was calculated at approximately R$ 1,164,845.91, providing an optimized solution for flood control in the watershed’s vulnerable areas. Research Implications: The results suggest a practical and feasible solution for flood control in urban watersheds, offering a replicable model that can aid in planning drainage infrastructure in other regions facing similar challenges. Originality/Value: This study is original in its integrated application of PCSWMM and HEC-RAS models with a nonlinear optimization approach in GAMS, resulting in an efficient and economical solution for flood mitigation in urbanized areas.

Automatic topology and capacity generation framework for urban drainage systems with deep learning-based land use segmentation and hydrological characterization

This paper proposed a comprehensive and automated framework that integrates the entire urban drainage design process, from refined land use segmentation, to catchment subdivision and characterization, to network topology generation, and finally to hydraulic estimation of pipeline capacity, based on multiple open-source datasets. The method incorporated deep learning-based semantic segmentation, GIS-based hydrological characterization, topology generation and hydraulic calculation algorithms. The performance and applicability of the method were validated through case studies of two cities. The results showed that the deep learning model with prior knowledge could accurately segment urban land use to obtain the detailed shape and boundary description, and achieved higher detection accuracy for water bodies, buildings, and roads. The stormwater drainage topology and conveyance capacity were automatically generated based on subcatchment division and hydrological characterization under a set of topological parameter conditions. In the optimal Critical Success Index (CSI) scenario, only 8% and 21% of the generated pipelines were inconsistent with the expert-designed networks. Compared to the traditional manual design process, the automated approach can save significant time and resources and reduce the reliance on field surveys and empirical/expert work. At the same time, the tool can be employed for system layout and capacity optimization, providing important support for the automation, design, and rehabilitation of stormwater drainage systems.

Investigating a surface water quality monitoring approach for QwaQwa, South Africa, by combining biological in vitro tests and chemical analyses

IntroductionWater quality assessment is a fundamental requirement for comprehensive water management. Access to water is not exclusively a quantitative issue, as human activities often lead to negative impacts on water quality due to pollution. This results in a need for the systematic monitoring of water bodies to surveil ongoing pollution dynamics and, if needed, facilitate the implementation of suitable solutions for water quality management.MethodsTo that effect, in March 2022 and July 2023, we conducted two sampling campaigns in four headwater streams in QwaQwa, South Africa, to assess the water quality under high (summer) and low (winter) flow conditions. The overarching aim was to support local decision-makers with relevant information on water quality in a semirural and less studied area. Beside the hydrological characterization of the catchment, which drains the mountainous areas around Phuthaditjhaba (the main town of the region), physicochemical and biological monitoring were conducted, including in vitro tests detecting endocrine disruption, dioxin-like impacts, and genotoxic potentials in the water.ResultsThe elemental analysis showed that at several sampling sites, Al, Cu, Mn, Fe, and Zn were above the Target-Water-Quality-Range (TWQR) for aquatic ecosystems and Fe, Al, and Mn were additionally above the TWQR for domestic use. Interestingly, the organic micropollutants (OMPs) revealed a region with a high incidence of HIV infections, given the occurrence of the drugs efavirenz, its transformation product 8-OH-efavirenz and nevirapine, among other micropollutants. Possibly, the signals from the in vitro tests for endocrine disruption are also induced by the detected micropollutants. Our results ultimately show that even though anthropogenic impacts increase from upstream to downstream, none of the sampling sites are without concern. This indicates the urgent need for measures to increase the water quality of these headwater streams in this region.

A Sponge Village Flood Response Method Based on GIS and RS Analysis Formation—A Case Study of Jiangou Village

This study was conducted in response to the Beijing–Tianjin–Hebei mega heavy rainfall event at the end of July 2023, and the severely affected and representative Jiangou village in Beijing was selected as the study area. A variety of methods were used to synthesize and analyze the situation and propose an adaptive response to heavy rainfall and flooding in the village. Based on multi-source remote sensing (RS) data, a comprehensive topographic and hydrological characterization was carried out, and the precipitation before and after the disaster was analyzed; the flood inundation area was extracted using the improved normalized water body index (MNDWI) and OTSU thresholding methods, and the changes of water bodies during the flooding period were quantitatively analyzed; and an improved convolutional-neural-network-based building identification and extraction model was constructed to extract the research distribution of buildings in the area. The sponge city construction (SPCC) method was improved to obtain a method that can mitigate flood risk and adapt to villages by constructing small artificial lakes and local topographic buffers to improve the water storage and drainage capacity of villages. The study shows that these methods are innovative in flood hazard analysis and mitigation but still need further improvement in data accuracy, simulation depth, and system evaluation.

Hydraulic analysis for assessing environmental flow selection and ecological model formulation

AbstractEcosystem management depends on transforming qualitative observations (e.g., slow‐moving shallow conditions provide nursery refugia for silvery minnow larvae) into management actions to increase habitat quantity or improve habitat quality. To be effective, decision metrics that are developed for management objectives should be validated with field observations. Model assumptions, precision and parameter importance can be refined by comparing the fidelity of selected parameters computed as habitat quality metrics and the correlation of these metrics to real‐world observations. Validated environmental metrics are more credible for management and can be compatible with ecosystem monitoring and project design processes. In this study, streamflow monitoring data and hydraulic modelling are used to quantify fish habitat extent for 15 years of spring runoff. The spring runoff event coincides with larval maturation to a free‐swimming juvenile phase for the silvery minnow, a critical period in Rio Grande habitat management. Different methods to estimate habitat availability (i.e., hydrology statistics, inundation extents based on hydraulic modelling and areal habitat availability based on different formulations of a habitat suitability index curve) were used to test the efficacy of different metrics relative to species population monitoring. This analysis finds that flow–ecology relationships based on hydraulic modelling or hydrology statistics are both effective and highly correlated to larval production. The investigation shows how seasonal hydrologic characterization and hydraulic discretization have varying levels of correlation with seasonal fish production. This study demonstrates how hydraulic modelling data and hydrologic characterization of riverine environments can be used to validate or develop conceptual ecological models.

Geotechnical Characterisation of Flysch-Derived Colluvial Soils from a Pre-Alpine Slope Affected by Recurrent Landslides

Heterogeneous rock masses that include rhythmic alternations of marl, shale, marly limestone, sandstone, siltstone, and argillite, such as Flysch, are particularly prone to generating colluvial deposits on gentle slopes, which are often subject to failures triggered by heavy rainfall. Flysch-derived colluvial soils are made up of highly heterogeneous sediments ranging from clayey loam to rock fragments, and they have been studied more rarely than homogeneous soils. In this work, we present a geotechnical and hydraulic characterisation performed both in situ and in the laboratory on flysch-derived colluvial soils that were involved in a channelised landslide in the pre-alpine area of the Friuli Venezia Giulia region (NE Italy). The investigated soils were characterised by the average values of the grain size composition of about 25% gravel, 20% sand, 30% silt, and 25% clay. The loamy matrix presented low-to-medium values of the liquid and plastic limits, as well as of the plasticity index (LL = 40%, PL = 23%, and PI = 17%, respectively). The values of the peak friction angle for natural intact samples were 33° < ϕ’p < 38°, whereas the residual friction angle fell to 23–24° at great depths and high vertical stresses, for a prevailing silty–clayey matrix. Variable head permeability tests were performed both in situ and in the laboratory, showing that the values of the vertical and horizontal permeability were very close and in the range 1 × 10−4–1 × 10−6 m/s. The soil permeability measured in the field was generally higher than the hydraulic conductivity calculated on laboratory samples. The proposed geotechnical and hydrological characterisation of flysch-derived colluvial soils can be of fundamental importance before the use of more thorough analyses/models aimed at forecasting the possible occurrence of slope failures and evaluating the related landslide hazard. The reported geotechnical and hydraulic parameters of flysch-derived colluvial materials can represent a useful reference for rainfall infiltration modelling and slope stability analyses of colluvial covers that are subject to intense and/or prolonged precipitation. However, when facing engineering problems involving colluvial soils, particularly those coming from flysch rock masses, the intrinsic variability in their grain size composition, consistency, and plasticity characteristics is a key feature and attention should be paid to the proper assumption of the corresponding geotechnical and hydraulic parameters.

Hydrological Characterization of Mediterranean Catchments

Hydrological Characterization of Mediterranean Catchments

ERT-based experimental integrated approach for soil hydrological characterization in rainfall-induced shallow landslides prone areas

Shallow landslides triggered by heavy rainfalls are slope instabilities, developed in the most superficial eluvial layers, involving the first 2 m from the ground level. A crucial predisposing factor in shallow landslides occurrence is the soil water content, generally measured trough sensors installed in the first soil layers. However, despite being a very precise approach, this monitoring technique provides for a site-specific dataset. An integrated method to extend the hydrological characterization from site-specific to slope scale is presented, combining geotechnical analyses, field data monitoring, and geophysical investigations, in two experimental test sites located on Italian Apennines. Ten Electrical Resistivity Tomographies (ERT) of the first soil horizons were performed through different array geometries (2D-3D-Time-Lapse), calibrated and interpreted basing on stratigraphic logs, trenches, and monitored soil water content field data. The test sites colluvial covers composition was analyzed and compared to resistivity values to build conceptual hydrogeological models of the deep-water circulation. In addition, two time-lapse (4D) ERT surveys were performed in both test sites simulating very intense precipitations, to determine the resistivity variations at different soil drainage conditions, thus estimating the average bulk permeability. Bulk permeability can be also a useful input parameter for slope stability models, widely employed in engineering practices. This integrated method proved to be very useful for the hydrogeological characterization of the subsoil at slope scale, where it is susceptible to slope instability, improving the knowledge of water circulation, as well as the bulk permeability heterogeneities, which are shallow landslides triggering parameters.

EDITH: An R‐package to spatially project eDNA‐based biodiversity across river networks with minimal prior information

Abstract Ecological and ecosystem monitoring is rapidly shifting towards using environmental DNA (eDNA) data, particularly in aquatic systems. This approach enables a combined coverage of biodiversity across all major organismal groups and the assessment of ecological indices. Yet, most current approaches are not exploiting the full potential of eDNA data, largely interpreting results in a localized perspective. In riverine networks, by explicitly modelling hydrological transport and associated DNA decay, hydrology‐based models enable upscaling eDNA‐based diversity information, providing spatially integrated inference. To capitalize on these unprecedented biodiversity data and translate it into space‐filling biodiversity projections, a streamlined implementation is needed. Here, we introduce the eDITH R‐package, implementing the eDITH model to project biodiversity across riverine networks with minimal prior information. eDITH couples a species distribution model relating a local taxon's eDNA shedding rate in streamwater to environmental covariates, a mass balance expressing the eDNA concentration at a river's cross‐section as a weighted sum of upstream contributions, and an observational model accounting for uncertainties in eDNA measurements. By leveraging on spatially replicated eDNA measurements and minimal hydro‐morphological data, eDITH enables disentangling the various upstream eDNA sources, and produces space‐filling maps of a taxon's spatial distribution at any chosen resolution. eDITH is applicable to both eDNA concentration and metabarcoding data, and to any taxon whose DNA can be retrieved in streamwater. The eDITH package provides user‐friendly functions for single‐run execution and fitting of eDITH to eDNA data with both Bayesian methods (via the BayesianTools package) and non‐linear optimization. An interface to the DHARMa package allows model validation via posterior predictive checks. Necessary preliminary steps such as watershed delineation and hydrological characterization are implemented via the rivnet package. We illustrate eDITH's workflow and functionalities with two case studies from published fish eDNA data. The eDITH package provides a user‐friendly implementation of eDITH, specifically intended for ecologists and conservation biologists. It can be used without previous modelling knowledge but also allows customization for experienced users. Ultimately, eDITH allows upscaling eDNA biodiversity data for any river globally, transforming how state and change in biodiversity in riverine systems can be tracked at high resolution in a highly versatile manner.

Hydrological modeling for the discharge production in an experimental area of the Guayas river basin, Ecuador

Hydrological simulation is a vital issue within the comprehensive management of hydrographic basins, with estimates being an input for decision-making by territorial managers. The objective of the current study is to conduct the hydrological modeling of an experimental basin of the Guayas River, where the discharge of the Baba River will be simulated using the HEC-HMS hydrological modeler. Hereby, the morphometric and hydrological modeling of the basin was performed from a global DEM (ASTER GDEM) in an environment of Geographic Information Systems (GIS) complemented with the analysis of thematic coverage of land types and uses as well as of precipitation data for the sector. The simulation of the discharge of the experimental basin determined that twenty percent of the precipitation that fell in the time window of analysis generated direct runoff and the residual eighty percent is attributed to abstractions. This is a condition that is given by the presence of agricultural crops corresponding to seventy percent forest masses, native forest and shrubby vegetation, as well as by the interception processes as well as water catchment systems in the urban area. The adequate results obtained in the application of the HEC - HMS model are based on the previous work of physical and hydrological characterization of the basin expressed through the physiographic structure of the basin and the valorization of the parameters that characterize the hydrological processes. This research tends to be useful for small hydrographic basins in Ecuador that play an important role in the ecosystem and lack of hydrometeorological information.

Automated Surface Runoff Estimation with the Spectral Unmixing of Remotely Sensed Multispectral Imagery

This work presents a methodology for the hydrological characterization of natural and urban landscapes, focusing on accurate estimations of infiltration capacity and runoff characteristics. By combining existing methods from the literature, we created a systemic process that integrates satellite-based vegetation maps, topography, and soil permeability data. This process generates a detailed vegetation classification and slope-corrected composite curve number (CNcα) map using information at the subpixel level, which is crucial for estimating excess runoff during intense precipitation events. The algorithm designed with this methodology is automated and utilizes freely accessible multispectral imagery. Leveraging the vegetation–impervious–soil (V-I-S) model, it is assumed that land cover comprises V-I-S components at each pixel. Automated Music and spectral Separability-based Endmember Selection is employed on a generic spectral library to obtain the most relevant V-I-S endmember spectra for a particular image, which is then employed in multiple endmember spectral mixture analysis to obtain V-I-S fraction maps. The derived fractions are utilized in combination with the Normalized Difference Vegetation Index and the Modified Normalized Difference Water Index to adapt the CNcα map to different seasons and climatic conditions. The methodology was applied to Esch-sur-Alzette, Luxembourg, over a four-year period to validate the methodology and quantify the increase in the impervious surface area in the commune and the relationship with the runoff dynamics. This approach provides valuable insights into infiltration and runoff dynamics across diverse temporal and geographic ranges.

Land Use and Land Cover Trends and Their Impact on Streamflow and Sediment Yield in a Humid Basin of Brazil’s Atlantic Forest Biome

Understanding the trends in land use and land cover (LULC) is crucial for modeling streamflow and sediment yield, particularly in hydrological basins. This study examined the impact of LULC on the dynamics of streamflow and sediment yield within a humid tropical basin of the Atlantic Forest biome in Brazil, focusing on the period from 2000 to 2016. Changes in LULC were analyzed using annual MapBiomas data products for the same period. The Soil and Water Assessment Tool (SWAT) model was deployed to simulate streamflow and sediment yield based on LULC changes. To investigate temporal trends in LULC, a suite of non-parametric statistical tests, including the Mann–Kendall, Pettitt, and Sen’s slope estimator tests, was employed. Ecological diversity indices such as Shannon–Weaver, Simpson, and Pielou were applied to assess forest fragmentation, along with the Forest Fragmentation Index. The results revealed a growing trend in urban and sugarcane areas, coupled with a decline in dense vegetation, mangroves, and other forms of dense vegetation. With regard to the correlation between land uses and hydrological variables, the findings indicate minor variations in hydrological balance, attributable to the not-so-significant changes among the studied land-use scenarios, except for sediment yield estimates, which showed more considerable alterations. Notably, the estimates for 2000 and 2013–2016 were the most divergent. In a broader scientific context, this research conclusively establishes that the incorporation of dynamic LULC data into the SWAT model augments the precision and robustness of simulations pertaining to agricultural watersheds, thereby enabling a more comprehensive hydrological characterization of the study area.

An integrated reservoir operation framework for enhanced water resources planning

Climate change induced spatiotemporal variation in global water availability modifies the proposed design criteria of water infrastructure structures like dams and reservoirs. Although reservoir operation is treated as a potential adaptation option, obsolescence of existing operation rules in the climate change scenarios could cause devastating situation through faulty water management practices. Presently onboard simulation–optimization based reservoir operation schemes fail to capture the uncertainty involved in the climate change scenario. Hence, there is a need to identify the limiting application scenario of the existing reservoir operation rule, and subsequently, revise the operation framework to address the future supply–demand uncertainty adequately. This research develops an integrated Soil and Water Assessment Tool (SWAT) (hydrologic), HEC-ResSim (hydraulic), and genetic algorithm (GA) (optimization) based adaptive reservoir operation framework, which is competent enough in accounting the future supply–demand uncertainty. Incorporation of the newly proposed environmental flow assessment approach in the reservoir operation would assist the decision makers in guiding the reservoir release for maintaining the water quality and sustenance of the downstream aquatic species. Certainly, corresponding to the existing operation rules under both the baseline and future climate change scenarios of RCP 4.5 and 8.5, the developed SWAT-HEC-ResSim-GA based reservoir operation scheme could improve the performance of the Kangsabati reservoir with the time and volume reliability estimates of 0.631 and 0.736, respectively. Conclusively, the developed approach in this study could be the best feasible alternative for hydrologic characterization in complex reservoir catchment-command regions with the option for enhanced reservoir planning in global catchment-command regions.

Calibrating Complexity: A Comprehensive Approach to Developing Stage–Storage–Discharge Relationships for Geographically Isolated Wetlands (GIWs) in W-C Florida

Geographically isolated wetlands (GIWs) play a critical role in regional hydrology, streamflow, groundwater recharge, evapotranspiration, and water quality. Despite their importance, research on quantifying discharges from these wetlands remains scarce. This study focuses on the hydrological characterization of GIWs in W-C Florida, a region where they cover a significant proportion of the land surface. The paper introduces a new procedure for developing Stage–Storage–Discharge (rating) characteristics essential for deterministic hydrological modeling in larger geographic areas. The approach synthesizes and extends previous research methodologies and offers simplified key modeling coefficients (m and n), which act as intelligent calibration knobs. These coefficients, when coupled with easily derivable physical attributes such as areas, slopes, and elevations, allow for the accurate simulation of downstream discharge timings and magnitudes, including flood behavior. The proposed method was tested using observational data from well-calibrated models. The results indicated a relative error of −7.2% for stage–storage and 15% for stage–discharge. GIS-based techniques were evaluated against surveyed data, and the results showed an accuracy of 0.16 m (0.52 ft) in estimating both the invert elevation and the maximum depth of GIWs. This effort will ultimately contribute to a better understanding and management of these critical ecosystems.

Seasonal hydrology and gas transport in a composite cover on sulfide tailings

This study presents the field performance of a five-layer composite cover to mitigate acid mine drainage in legacy sulfide tailings in northern Ontario, Canada. Installed in 2008, this cover comprised sand, clay, geosynthetic clay liner, sand, and waste rock layers. To evaluate the effectiveness of the cover in reducing water and oxygen ingress, groundwater and vadose zone hydrological characterization, stable water isotope analysis, pore-gas measurements, oxygen flux calculations, and variably saturated flow modeling were conducted. The results indicate that the clay layer stayed nearly saturated in the spring, fall, and winter, but temporary desiccation occurred during the summer. Compared to uncovered tailings, the cover significantly lowered diffusive oxygen flux. In the summer, fall, and winter, the capillary barrier effect of the cover functioned effectively and inhibited percolation. Atmospheric pore-gas oxygen concentrations at one out of the three monitoring locations indicate potential cover imperfections that enabled oxygen transport into the tailings. In the spring and early summer, snowmelt infiltration resulted in percolation that compromised the capillary barrier effect, as well as lateral drainage. The resulting increase in water saturation in the cover limited oxygen transport. Despite potential cover imperfections, this composite cover reduced oxygen and water ingress a decade after installation.

Rainfall intensity model with spatialization of intensity-duration-frequency curve parameters - A case study for the state of Maranhão, Brazil

The characterization of intense rainfall in engineering projects is fundamental, mainly regarding the estimate of design flows necessary for designing hydraulic works. Intense rainfall events are commonly measured by Equations and curves that relate their intensity, duration, and frequency. Such relations, known as IDF, enable the hydrological characterization of a given region. This article presents a methodological design and results from both determination and spatialization of IDF curve parameters for the state of Maranhão. Historical series of maximum daily rainfalls obtained from National Water and Sanitation Agency (ANA) were used in 126 rainfall gauge stations and the Gumbel probability distribution estimated the maximum daily rainfall for 5, 10, 15, 25, 50, and 100 return periods. The Isozonal Method obtained the IDF correlations of intense rainfall events for 0,1. 1, and 24 h durations, and their performance were conducted by Nash-Sutcliffe R2 coefficient and Root Mean Square Relative Error (RMSE). “K, a, b, and c” parameters of intense rainfall equations were determined by optimization and convergence processes and their spatialization was carried out by interpolation by Inverse Distance Weighted (IDW), which enabled to determine the values of each parameter in regions without physical measurements of rainfall. Similarly, rainfall intensity was spatialized for the entire state. According to the results, the rainfall distribution in the state of Maranhão shows a variation in the indexes of precipitation, with the highest values found in areas located in central-southern, southwestern, and southeastern regions.

Hydrological Characterization and Statistical Analysis of Groundwater Quality for the Ogallala Aquifer in the Permian Basin, Texas, USA

Hydrological Characterization and Statistical Analysis of Groundwater Quality for the Ogallala Aquifer in the Permian Basin, Texas, USA

How does flow connection path and vertical spatial layout of LIDs affect urban runoff? A new LID construction method based on refined landuse and hydrologic characterization

The techniques and modeling methods of Low Impact Development (LID, a more decentralized approach that uses or mimics natural processes to mange stormwater), has been extensively studied due to their positive impacts on urban flooding and water environment. Traditional LID modeling has been typically allocated on a plane layout based on rough land use classification and hydrological characterization. As the LIDs cannot be combined in the same space due to modeling restrictions, the LID effects on runoff control are often underestimated. This paper proposed a vertical LID model construction method based on a refined land use classification and improved hydrological characterization. By establishing a flow connection flux in the LID controller calculations, the overflows of different LIDs can be interconnected and transferred for further infiltration, storage and conveyance. Results showed that the refined approach provided more accurate spatial allocation guidance for LIDs, and the obtained catchment outflow was reduced by 25%-34% in comparison to conventional land use classification in a case study in northern China. The vertical construction method allowed a multi-dimensional LID allocation in space and function and achieved a higher control effect on total outflow volume and peak value. The combined LIDs could cope with the rainfall (i.e., no outflow from the LID plots) within the 50-year event, and the total outflow volume was 13–26% lower than the one from the plane layout method.

Biochar Physical and Hydrological Characterization to Improve Soil Attributes for Plant Production

Biochar Physical and Hydrological Characterization to Improve Soil Attributes for Plant Production

Evaluation of hydrological responses to climate change for a data-scarce mountainous watershed in Taiwan

Abstract Despite the mountainous watersheds being important for the ecosystem, water resources, and hydropower, little hydrological data has been collected. This data scarcity makes it difficult to evaluate their hydrologic response to climate change. This study integrated short-term hydrological data with physics-based meteorological and hydrological models to measure the impact of climate change on future water scarcity in the Wuling Mountain Watershed, Taiwan. Twenty-eight-month hydrological datasets from 2013 to 2015 were used for the hydrological characterization. Scenarios from CMIP5 were selected for the climate projection for the period 2021–2040 based on the inferred 1986–2005 baseline data. The results showed that precipitation, water percolation, and streamflow will decrease by about 10% and increase by about 20–25% in the dry and wet seasons, respectively. The evapotranspiration is lower than that of the baseline in January, March, and December, whereas it can be as high as 4% during the other months. The increase in the annual amount and change in the water distribution for all studied water components indicate the possible acceleration of the water cycle. Spectrum analysis showed that surface water tends toward becoming more irregular. Groundwater remains mildly persistent and thus may serve as a buffer for the impact of climate change on water resources.