Journal of Hydrology最新文献

Water distribution systems (WDS) and power grids (PG) are critical infrastructure systems that are vital to all human activity. As such, their quality of service is of great importance for economic, environmental, and human welfare reasons. Although traditionally being analyzed separately, the two systems are interconnected and can mutually affect one another. In order to utilize the potential benefits that the two systems can produce for each other, their design and operation should be analyzed conjunctively. In this paper, a conjunctive optimal design approach for water and power networks is presented, with the objective of finding the dimensions of the systems’ facilities that will result in minimal overall costs, for both design and operation. The model is formulated and implemented on two example applications using an off-the-shelf nonlinear solver by MATLAB and compared to the optimal design of the independent WDS. A sensitivity analysis is performed to provide validity to the obtained results. The conjunctive design is compared to the design of an independent WDS to emphasize the effect of including the PG in the optimization problem. Results show a clear link between the availability of renewable energy and sizing of WDS components. The design of the independent WDS leads to the violation of PG constraints, which are satisfied when including both systems under a single optimization model, demonstrating the importance of a holistic design approach.

To address the challenge of escalating urban flood risk and the deficiency in effective flood emergency management, this study introduces a novel Coupled Human and Natural Systems (CHANS) modelling framework that employs hierarchical reinforcement learning to optimise mobile pump scheduling and placement for urban flood risk mitigation. The CHANS framework integrates hydrodynamic and agent-based models within a multi-GPU computing environment for high-resolution, real-time flood inundation modelling and risk assessment to enrich Reinforcement Learning (RL) training. In the application to Ninh Kieu District in Can Tho City, Vietnam, the new RL-enabled modelling framework is used to evaluate optimal mobile pumping strategies for concurrent pluvial flooding and post-flooding events against the traditional deployment approaches. Results demonstrate that RL-based strategies can significantly enhance flood risk reduction, outperforming traditional methods by achieving 2× and 4× improvements in the concurrent and post-flooding periods/events, respectively. Incorporating human factors and adapting to local conditions, the RL agent provides valuable insights into mobile pump scheduling and deployment strategies. Sensitivity analysis confirms the robustness of the CHANS modelling framework and underscores the role of RL in optimising mobile pump scheduling and placement where traditional rule-based strategies are challenging.

Land use change can significantly affect soil hydrology in arid and semi-arid regions, making it crucial to understand the relationship between vegetation roots and soil moisture. Current models often fail to predict root growth and its impacts on water dynamics accurately. Our work presents a novel model that seamlessly integrates the Community Land Model (CLM) with the Soil & Water Assessment Tool (SWAT). Furthermore, it enhances the root module within the CLM, enabling more accurate simulations of dynamic root depth and distribution across varying tree ages. This improvement particularly considers the crucial processes of dormancy and plant maturity. Soil moisture and root patterns under apple trees of varying ages and in wheat fields on the Loess Plateau was analyzed. Our findings indicate that our dynamic root depth model outperforms traditional static models, and can accurately reflect soil moisture levels with high precision (R² = 0.80–0.81; Nash-Sutcliffe efficiency (NSE) = 0.65–0.75). In contrast to methods that utilize fixed root depths, dynamic root simulation can provide new insights. As apple orchards mature, the roots of 22-year-old apple trees have been found to reach a depth of 21 m in the soil. Conversely, the maximum root depth of wheat is limited to 1.9 m. This latter finding aligns more closely with the measured root depths, highlighting the accuracy of dynamic simulations. This model reveals that older apple orchards show decreased soil moisture at greater depths (>20 m), contrasting with wheat fields that affect moisture mostly within the top 2 m. Our results underscore the crucial role of dynamic modeling in comprehending root-soil water interactions. Furthermore, they imply that extended orchard cultivation practices can lead to a substantial depletion of deep soil moisture. Specifically, over a period of 1 to 22 years, a water deficit of up to 85 mm yr⁻¹ has been observed. For a 22-year-old forest, the D-D (dynamic distributions of coarse and fine roots) method calculates a significant cumulative deep Soil Water Storage loss. Over the course of 22 years, this loss amounts to 1664 mm, which is almost three times compare to the annual rainfall recorded. Such a large loss has the potential to significant impact on groundwater recharge. This highlights the need for careful consideration in future afforestation efforts to prevent increased soil aridity.

Coastal aquifers, the transition zone between freshwater and saltwater, show large salinity contrasts in the subsurface. Salinity is a key parameter to understand coastal groundwater flow dynamics and consequently also geochemical and microbial processes. For mapping porewater salinity, a variety of methods exists, mainly using electrical conductivity as a proxy. We investigate methods including hydrological/geochemical (well sampling, fluid logger) as well as geophysical method (direct push, geoelectrics) utilizing measurements near the high-water line of a high-energy beach at the North Sea island of Spiekeroog. We compare the methods, discuss their benefits and limitations and assess their spatial and temporal resolution. One key to enable a comparison is the estimation of formation factors transforming bulk conductivity measured by geophysical tools in to fluid conductivities obtained from direct measurements. We derive depth-dependent formation factors derived from time-series measurements of fluid loggers and a vertical electrode installation. Using these formation factors, the vertical electrode chain proves to provide reliable salinities at high spatial and temporal dimension. Direct-push profiling data provide the highest vertical resolution. However, a careful calibration is needed to allow for salinity quantification. On the other hand, electrical resistivity tomography (ERT) exhibits the lowest spatial resolution, but can image two-dimensional salinity distributions. We found ERT to fit very well to all other methods, but the data analysis should be aimed at salinities instead of bulk conductivities, i.e. including formation factors and temperature models into the inversion process.

Precipitation is an important piece of information needed in many areas such as transportation, military and agriculture, and microwave links have proven to be an effective means of acquiring it and can be used as a complementary means for professional precipitation measurement instruments such as rain gauges, weather radar, rain measuring satellites, etc. In this paper, two microwave links (at 26 GHz, both vertically polarised) located across a river in eastern China and a nearby OTT PARSIVEL disdrometer were used to carry out a precipitation observation experiment from October 2020 to March 2022 (March to November 2021 for liquid precipitation, others for possible non-liquid precipitation). The applicability of the dry and rainy period identification methods (standard deviation method and correlation coefficient method) based on this link data is analyzed first. In addition, a wet antenna attenuation correction model based on a Long Short-Term Memory (LSTM) neural network was applied to non-winter precipitation inversion. The results show that the rain rate and cumulative rainfall calculated from the two links are in good agreement with that calculated by the disdrometer data. Furthermore, quantitative inversion of winter precipitation using two microwave links is performed, and the variation characteristics of the link levels during non-liquid precipitation periods are also analyzed.

To minimize the loss of life caused by earthquakes, it is crucial to have early warning tools that provide sufficient warning time. Although various methods are available, their accuracy is uncertain. This review explores alternative indicators related to hydrogeochemical anomalies that appear before earthquakes, which could potentially offer earlier warnings. It presents a theoretical basis for ionic and molecular changes initiated by the generation of hydroxyl radicals (•OH) and hydrogen peroxide (H₂O₂) in groundwater before an earthquake. The review analyses 32 earthquake events by measuring ionic anomalies through standard reduction potential. The study also identifies significant pre-seismic erosional activity in minerals exposed to hydrogen peroxide. Additionally, the oxidation of organic compounds is shown to predict notable geochemical changes. The changes in the groundwater chemistry can also be validated by examining simultaneous microbial/ zooplankton responses. The findings suggest that the ionic, mineral, organic, and microbiological changes triggered by hydroxyl radicals and hydrogen peroxide could serve as a multi-parametric earthquake early-warning system. Implementation of such a system could improve the accuracy and timeliness of earthquake warning systems potentially reducing the associated risks and damages.

Understanding the spatial distributions of river-tide dynamics in estuaries and their response to intensive human interventions is critical. However, the studies on the characteristics of the spatial distributions of water level and tidal range are insufficient, with inadequate direct established empirical formulas. In this study, we propose a general and analytical water level and tidal range distribution model based on the recently proposed General Unit Hydrograph (GUH) theory and apply it to the Modaomen Estuary in China. Due to the intensive human interventions, the evolution of river-tide dynamics in the Modaomen Estuary was divided into three distinct periods: Pre-human, Transitional and Post-human Periods. The results show that the water level increased at the SZ station, yet decreased at other stations, thereby reducing the maximum water level gradient by approximately 9.08 × 10⁻⁶ on average from the Pre-human to the Post-human Periods. The tidal range generally increased along the channel and the absolute value of the maximum tidal range gradient decreased by approximately 54 % on average. Correspondingly, the positions of both maximum water level gradient and minimum tidal range gradient generally move seaward by approximately 27.25 km and 0.34 km, respectively. Additionally, the spatial distributions of river-tide dynamics tended more linear after human interventions. The satisfactory correspondence of the model outputs with observations at the six gauging stations along the Modaomen Estuary (with maximum RMSE values being 0.05 m for the water level and 0.03 m for the tidal range, respectively) indicate that the newly proposed GUH model could serve as a simple tool to describe the spatial distributions of water level and tidal range in a specific channel together with their response to anthropogenic modifications, which is also particularly useful for the identification of regime shifts in river-tide dynamics and sustainable water resources management in other highly human-modified estuaries worldwide.

Energy balance distributed modelling in High Mountain Asia (HMA) is important to examine glaciological and hydrological processes and assess changes in streamflow in the current and future climate. In this study, the Physically based Distributed Snow Land and Ice Model (PDSLIM) using detailed observed meteorological data at hourly scale is employed to simulate the hydrological response of the Naltar catchment, 242.62 km² in size, (in the Karakoram region in Pakistan to simulate its glaciers’ mass balance as well as daily runoff. The results exhibited overall satisfactory performance in terms of coefficient of determination (R² = 0.96) and Nash-Sutcliffe Efficiency (NSE=0.95) modelled against satellite-based snow cover areas, for internal model verification, in eight years. The results of runoff simulations compared for external model verification, with observed daily discharge resulted in NSE 0.90 and 0.89 for calibration and validation period respectively. Flow composition analysis revealed that the streamflow regime of Naltar catchment is composed to 40 % by glacier runoff, 42 % by sub-surface runoff and 18 % by surface runoff. The eight year mean value of net mass balance exhibited a slightly negative mass balance (−0.810 ± 0.31 m w.e. a^-1) less pronounced than that observed globally in several continental glaciers distinct from the Greenland and Antarctic ice sheets in the current climate. It seems that the ‘Karakoram anomaly’, i.e. the balanced to slightly positive glacier budgets observed in the region in the recent decades, a unique dynamics worldwide, has a moderate impact in the central Karakoram. Overall, the distributed energy-balance model PDSLIM, so far tested in the Alps, results to be a suitable tool to estimate energy and mass balance in the glacierized catchments of Karakoram and Himalaya and to better understand snow and ice melt runoff dynamics and floods in highly complex and glacierized mountain basins in the current and, in our research perspective, in the future climate.

The flow patterns within a longitudinal gap area formed by discontinuous distributions of submerged canopy, as well as the momentum and mass exchange characteristics between the gap area and the overlying free-flow, were studied using high-resolution Large Eddy Simulation (LES). The gap area is located within the fully developed region of submerged canopy flow. The simulations considered four aspect ratios of the gap area, L/h (ranging from 1 to 4), where L and h represent the span of the gap area and the height of the canopy, respectively, and two canopy densities, ϕ = 0.08 and 0.15. Results indicate that recirculation vortices appear only at ϕ = 0.15 and penetrate the downstream canopy patches to varying extents, with the distance of invasion decreasing as L/h increases. Influenced by the recirculation vortices, the bed shear stress in the downstream part of the gap area and the initial section of the downstream canopy patch is significantly increased compared to the fully developed region of the upstream canopy patch. Within the parameter range covered, vertical penetration of the mixing layer into the gap area always occurs, consistently falling into the shear layer growth regime, with significant enhancement of turbulence within the gap area even at L/h = 1. The high momentum entrainment elevates longitudinal velocity and turbulence intensity in the upper corner regions of the downstream canopy patch’s leading edge, combined with localized increases in bed shear stress, potentially destabilizing plants at the forefront of the downstream canopy patch. Although the total momentum exchange across the interface between the gap area and the upper free-flow layer is approximately independent of L/h, larger canopy densities lead to stronger momentum transport, and turbulent transport always dominates. Mass exchange generally increases with L/h, with more efficient vertical mass exchange in ϕ = 0.15 cases at L/h = 3 and 4.