Circular Image

M. Hrachowitz

info

Please Note

98 records found

Journal article (2026) - Kailong Li, Saman Razavi, Holger R. Maier, Markus Hrachowitz, Ehsan Nabavi, Natasha Harvey, Khaled Akhtar, Fisaha Unduche
The development of AI models is increasing at a rapid rate. However, when are they ready to be deployed in real-world operational settings? In this paper, we introduce a framework to support such assessments and apply it to Google’s recently released AI-based flood prediction system, which is claimed to achieve “reliability in predicting extreme riverine events” and provide “accurate and timely warnings” that are available “earlier and over larger and more impactful events in ungauged basins”. The system has been integrated into an operational early-warning platform producing open, real-time forecasts in more than 80 countries. While this development promises to usher in a new and exciting age in global flood forecasting, the supporting evidence relies heavily on several subjective choices, the implications of which have not been acknowledged or assessed. Here, we evaluate the consequences of these choices on claims of operational deployment readiness across four dimensions: predictive accuracy, forecast timeliness, the characterization of extreme events, and benchmarking against state-of-the-art models. Our assessment reveals that the system’s actual predictive accuracy is likely to be substantially lower than reported—particularly for extreme events—raising concerns about responsible practices across modelling and publicity in high-stakes applications. The deployment of the Google AI model therefore risks misinforming those who depend on its outputs for evacuation and preparedness decisions, particularly in less-developed countries such as those targeted by the enterprise, given its alarmingly high (>90%) rates of false positives and false negatives. Beyond the immediate operational consequences, if left unaddressed, these outcomes may erode public trust in AI within hydrological sciences. We conclude by calling for greater transparency, accountability, and methodological rigor in the integration of AI into flood forecasting. ...
Journal article (2026) - Hatice Turk, Christine Stumpp, Markus Hrachowitz, Peter Strauss, Günter Blöschl, Michael Stockinger
Preferential flow paths (e.g., macropores or subsurface pipe networks) in hydrological systems facilitate the rapid transmission of precipitation and solutes to streams, resulting in streamflow responses characterized by the release of younger water (i.e., recent precipitation) from the catchment and correspondingly short transit times (on the order of days). While preferential flow paths are documented in both the unsaturated zone and groundwater aquifers, it remains uncertain whether catchment-scale isotope-based transport models can adequately represent preferential flow using tracer measurements in streamflow. In this study, we hypothesize that the preferential release of young water from both the unsaturated zone and groundwater aquifers can be isolated from the streamflow tracer signal. This can be studied with StorAge Selection (SAS) functions, which describe how young or old water leaves a storage. We systematically compared multiple parameterizations of SAS functions describing how water of different ages is released from the unsaturated zone and groundwater aquifer within a single catchment-scale transport model using long-term measurements of hydrogen isotopes in water ( δ2H) from two headwater catchments (the Hydrological Open Air Laboratory (HOAL) in Austria and the Wüstebach catchment in Germany). The results show that δ2H measurements in streamflow exhibited sufficient variability to isolate the preferential release of younger water through preferential flow paths in the unsaturated zone. In contrast, the variability of δ2H in streamflow was insufficient to isolate the preferential release of younger water from the groundwater aquifer, as any seasonal variations in pore water δ2H were largely damped by substantial passive groundwater storage (water that mixes with the tracer signal of the active groundwater volume). Consistent with this interpretation, the degree of attenuation in the simulated streamflow isotope signal increased with increasing passive groundwater storage volumes and became pronounced when passive storage was orders of magnitude larger than active groundwater storage. The size of passive groundwater storage, in combination with groundwater SAS function parametrizations, regulated the long tails (100<T<1000 d) of transit time distributions, resulting in considerable uncertainty (± 20 % for HOAL and ± 23 % for Wüstebach) in the fraction of streamflow older than 100 d. The findings demonstrate that stable water isotope measurements from streamflow outlets is insufficient to constrain preferential groundwater flow in the two study catchments and plausibly in similar catchments characterized by large passive groundwater storage. The variability in streamflow TTD estimates arising from different groundwater storage SAS function parametrizations is considerable. Reducing uncertainty in groundwater transit time estimates and preferential flow contributions to streamflow requires complementary data sources, including multiple tracers, high-frequency tracer analysis, and groundwater-level observations, to improve catchment-scale transit time modelling. ...
Journal article (2026) - Thiago V.M. do Nascimento, Julia Rudlang, Sebastian Gnann, Jan Seibert, Markus Hrachowitz, Fabrizio Fenicia
Although large-sample hydrology data sets are increasingly used to advance predictions in ungauged basins, the influence of landscape data quality on model regionalization remains insufficiently explored. This study investigates whether geological catchment attributes derived from maps of increasing detail—global, continental, and regional—improve parameter transfer and model regionalization. To ensure robustness across model approaches, we applied both a semi-distributed process-based hydrological model using hydrological response units (HRUs) and a data-driven Long Short-Term Memory (LSTM) model. The analysis covered a total of 130 catchments in the Moselle (Central Europe) and Garonne (southwestern France) basins. We conducted five model experiments differing only in the representation of geological information: a benchmark without geology, a benchmark with random geology classes, and configurations based on the global-, continental-, and regional-scale geological maps. Model performance was evaluated using a modified Nash-Sutcliffe (NSE) metric for daily streamflow, as well as Pearson correlation and relative bias for three streamflow signatures: baseflow index, slope of the flow duration curve, and half-flow date. Across both basins and modeling frameworks, increasing geological detail consistently improved predictive performance under space–time evaluation. While differences in NSE were modest, improvements were pronounced for streamflow signatures: only models using the more detailed geological information, especially the regional map, consistently reproduced spatial variability in baseflow and flow regime characteristics. These findings highlight the importance of integrating high-quality geological data into hydrological modeling, particularly for improving predictions in ungauged basins through more reliable parameter transfer and regionalization. ...
Journal article (2026) - M. Ibrahim, Fransje van Oorschot, R.J. van der Ent, M. Hrachowitz, Miriam Coenders-Gerrits
Quantification of long-term partitioning of precipitation into evaporation and runoff is a fundamental pursuit in catchment hydrology. The Budyko framework provides a theoretical basis for this and estimates the evaporative fraction based on the aridity index. However, deviations from the global-average Budyko curve point to additional controls on precipitation partitioning beyond the aridity index. We hypothesized that root zone storage capacity (Sr,max), defined as maximum subsurface water volume accessible to vegetation roots, is a key driver of these deviations. The relationship between Sr,max and precipitation partitioning in the Budyko space was investigated globally across >5000 catchments. Sr,max was calculated using the memory method based on runoff observations and the water balance. The ω-parameter from Fu’s equation, which was used here to construct parametric Budyko curves, reflects deviations from the global-average Budyko curve and hence precipitation partitioning. Results revealed a globally stronger correlation (Spearman’s ρ= 0.68) of ω with Sr,max, than with other potential controls, indicating Sr,max as a dominant driver of precipitation partitioning. Further analysis based on Köppen–Geiger climatic zone classification revealed variations in the Sr,max–ω relationship, with the strongest correlations observed in cold (ρ= 0.87) and Mediterranean (ρ = 0.83) climates, followed by temperate (ρ = 0.76), tropical (ρ = 0.64) and arid climates (ρ = 0.61). Regional differences in Sr,max indicate that, at a given aridity, EA/P largely reflects vegetation adaptation to the seasonal interplay between water supply and atmospheric water demand. This study provides strong empirical evidence on a global scale for Sr,max as a governing factor in modulating catchment precipitation partitioning, as evident in the Budyko space. As a major implication our results provide a theoretical basis for the maximum values of Sr,max found in nature, as constrained by the water and energy limits of the Budyko framework. ...
Journal article (2025) - Qiaojuan Xi, Hongkai Gao, Lan Wang-Erlandsson, Jianzhi Dong, Fabrizio Fenicia, Hubert H.G. Savenije, Markus Hrachowitz
Adaptation of ecosystems’ root zones to climate change critically affects drought resilience and vegetation productivity. However, a global quantitative assessment of this mechanism is missing. In this study, we analyzed high-quality observation-based data to find that the global average root zone water storage capacity (SR) increased by 11%, from 182 to 202 mm in 1982–2020. The total increase of SR equals to 1652 billion m3 over the past four decades. SR increased in 9 out of 12 land cover types, while three relatively dry types experienced decreasing trends, potentially suggesting the crossing of ecosystems’ tipping points. Our results underscore the importance of accounting for root zone dynamics under climate change to assess drought impacts. ...
Journal article (2025) - Hatice Türk, Christine Stumpp, Markus Hrachowitz, Karsten Schulz, Peter Strauss, Günter Blöschl, Michael Stockinger
The rainfall-runoff transformation in catchments usually follows a variety of slower and faster flow paths, leading to a mixture of "younger"and "older"water in streamflow. Previous studies have investigated the time-variable distribution of water ages in streamflow (transit time distribution, TTD) using stable isotopes of water (δ 18O, δ 2H) together with transport models based on Storage Selection (SAS) functions. These functions are traditionally formulated based on soil moisture to mimic the preferential release of younger water as the system becomes wetter. In this study, we hypothesized that, in a heterogeneous catchment with a significant fast-runoff response component, precipitation intensity, in addition to soil moisture, plays a critical role in the preferential release of younger water. To test this hypothesis, we used high-resolution δ 18O data (weekly and event-based streamflow δ18O samples) in a 66 ha agricultural catchment. We tested two scenarios of the SAS function parameterization for the preferential-flow age selection: one as a function of soil moisture only and one as a function of both soil moisture and precipitation intensity. The results showed that accounting for both soil moisture and precipitation intensity to define the shape of SAS functions for preferential flow improved the tracer simulation in streamflow (increasing the Nash-Sutcliffe efficiency from 0.31 to 0.51). This also led to a higher percentage of streamflow (an increase from 2.87 % to 4.38 %) with shorter transit times (TTs younger than 7 d), with the largest differences occurring during the summer and autumn months. This was due to the fact that incorporating both soil wetness and precipitation intensity in the SAS formulation accounts for rapid flow pathways such as infiltration excess overland flow, preferential flow through macropores, and tile drain flow - allowing precipitation water to bypass much of the soil matrix and to reach the stream with minimal storage or mixing, even under dry soil conditions. We showed for the agricultural study catchment that a significant portion of event water bypasses the soil matrix through fast-flow paths, resulting in younger water reaching the stream for both low- and high-intensity precipitation. Thus, in catchments where preferential flows and overland flow are the dominant flow processes, soil-wetness-dependent and precipitation-intensity-conditional SAS functions may be required to better describe the timescale of solute transport in modelling, which has implications for stream water quality and agricultural management practices such as the timing of fertilizer application. ...
Accurate rainfall observations with high spatial and temporal resolutions are key for hydrological applications, in particular for reliable flood forecasts. However, rain gauge networks operated by regional or national environmental agencies are often sparse, and weather radars tend to underestimate rainfall. As a complementary source of information, rain gauges from personal weather stations (PWSs), which have a network density 100 times higher than dedicated rain gauge networks in the Netherlands, can be used. However, PWSs are prone to additional sources of error compared to dedicated gauges, because they are generally not installed and maintained according to international guidelines. A systematic long-term analysis involving PWS rainfall observations across different seasons, accumulation intervals, and rainfall intensity classes has been missing so far. Here, we quantitatively compare rainfall estimates obtained from PWSs with rainfall recorded by automatic weather stations (AWSs) from the Royal Netherlands Meteorological Institute (KNMI) over the 2018–2023 period, including a sample of 1760 individual rainfall events in the Netherlands. This sample consists of the 10 highest rainfall accumulations per season and accumulation intervals (1, 3, 6, and 24 h) over a 6-year period. It was found that the average of a cluster of PWSs severely underestimates rainfall (around 36 % and 19 % for 1 h and 24 h intervals, respectively). By adjusting the data with areal reduction factors to account for the spatial variability of rainfall extremes and applying a bias correction factor of 1.22 to compensate for instrumental bias, the average relative bias decreases to −5 % for 1 h intervals or almost zero for intervals of 3 h and longer. The highest correlations (0.85 and 0.86) and lowest coefficients of variation (0.14 and 0.18) were found for 24 h intervals during winter and autumn, respectively. We show that most PWSs are able to capture high rainfall intensities up to around 30 mm h−1, indicating that these can be utilized for applications that require rainfall data with a spatial resolution of the order of kilometres, such as for flood forecasting in small, fast-responding catchments. PWSs did not observe the most intense rainfall events, which were associated with return periods exceeding 10 or 50 years (above approximately 30 mm h−1) and occurred in spring and summer. However, the spatial distribution of rainfall likely played a large role in the observed differences rather than instrumental limitations. This emphasizes the importance of having a dense rain gauge network. In addition, the variation in undercatch is likely partly due to the disproportional underestimation of tipping bucket rain gauges with increasing intensities. Outliers during winter were likely caused by solid precipitation and can potentially be removed using a temperature module from the PWS. We recommend additional research on dynamic calibration of the tipping volumes to improve this further. ...
On 29 October 2024, torrential rainfall locally exceeding 300 mm within less than 24 h, caused devastating floods in the province of Valencia in Spain. In this study we quantify and describe the spatial and temporal structure of the rainfall event on this day using rainfall observations from approximately 225 personal weather stations (PWSs), low-cost commercial devices primarily operated by citizens. The network density of PWSs is ∼7 times higher compared to the dedicated rain gauge network operated by the Spanish Meteorological Agency (AEMET) in the province of Valencia, allowing a more detailed analysis of the spatial and temporal rainfall dynamics. In addition, PWS observations are available in near real-time to the public with a temporal resolution of 5 min, whereas the data from AEMET are not available in real time for the public and at a lower publicly available temporal resolution (1 h). Daily rainfall sums recorded by the PWSs showed a high correlation (r=0.94) and low bias (underestimation of 4 %) compared to rainfall reported by AEMET. In the upstream parts of the Magro catchment (1661 km2), a first burst of extreme rainfall, reaching up to 180 mm of rainfall in a few hours, started in the morning, leading to the generation of a first flood wave in the upstream parts of the catchment. While the resulting flood wave was propagating downstream through the channel network, a second rainfall peak occurred, which moved downstream along with the flood wave. This spatial and temporal coincidence has likely exacerbated the devastating power of this event. Based on the PWS data, it could have been anticipated that the extreme rainfall already occurring early in the morning would likely result in flooding in the Magro catchment. Areal rainfall maps based on interpolating PWS data indicated catchment average rainfall exceeding 150 mm d−1 across an area of more than 2500 km2. However, the total accumulated rainfall remains uncertain due to interrupted measurements likely caused by power outage and inherent uncertainty associated with interpolating point measurements. For the Rambla de Poyo catchment, the resulting average discharge was around 900 m3 s−1. The estimated return period of the catchment-average rainfall and resulting discharge from this event exhibits large uncertainties, with on average exceeding 10 000 and 900 years, respectively. This study shows the potential of PWSs for real-time rainfall monitoring and potentially flood early warning systems, by complementing dedicated rain gauge networks in order to reduce the uncertainty from areal rainfall estimates and to localize potential flooding more accurately. ...
Journal article (2025) - Thiago V.M. do Nascimento, Julia Rudlang, Sebastian Gnann, Jan Seibert, Markus Hrachowitz, Fabrizio Fenicia
Large-sample hydrology datasets have advanced hydrological research, yet the impact of landscape map details on identifying dominant streamflow generation processes remains underexplored. This study investigates the role of geology using maps of increasing detail – global, continental, and regional – each reclassified into four permeability classes. These geological attributes were used along with topography, soil, land use, and climate attributes to identify dominant controls on streamflow signatures across 4469 European catchments. To distinguish landscape influences from the otherwise dominant influence of climate, we conducted separate analyses on nested basins. Three scales were considered to assess scale-dependent patterns: large (63 nested basins), intermediate (the Moselle nested basin), and small (five nested catchments within the Moselle). The large-scale study used geology information from global and continental maps, while the others also incorporated regional maps. At the large scale, dominant controls varied widely between nested basins, but landscape generally outweighed climate, highlighting the value of our nested basin design. At this scale, continental and global geology maps produced different correlation patterns, with neither consistently superior. At the intermediate scale, increased geological detail led geology to shift from the least to the most correlated variable for certain streamflow signatures. The small-scale experiment reinforced these findings, as the regional map highlighted controls more consistent with process understanding. This study underscores the benefit of integrating detailed, region-specific geological data into large sample hydrology studies, and demonstrates the utility of a nested basins design. These findings have important implications for hydrological regionalization and streamflow prediction in ungauged basins. ...
Review (2025) - Paul D. Wagner, Doris Duethmann, Maria Staudinger, Larisa Tarasova, Stephan Thober, Nicola Fohrer, Doerthe Tetzlaff, Thorsten Wagener, Björn Guse, Jens Kiesel, Sandra Pool, Markus Hrachowitz, Serena Ceola, Anna Herzog, Tobias Houska, Ralf Loritz, Diana Spieler
While measured streamflow is commonly used for hydrological model evaluation and calibration, an increasing amount of data on additional hydrological variables is available. These data have the potential to improve process consistency in hydrological modeling and consequently for predictions under change, as well as in data-scarce or ungauged regions. Here, we show how these hydrological data beyond streamflow are currently used for model evaluation and calibration. We consider storage and flux variables, namely snow, soil moisture, groundwater level, terrestrial water storage, evapotranspiration, and altimetric water level. We aim at summarizing the state-of-the-art and providing guidance for the use of additional hydrological variables for model evaluation and calibration. Based on a review of the current literature, we summarize observation methods and uncertainties of currently available data sets, challenges regarding their implementation, and benefits for model consistency. The focus is on catchment modeling studies with study areas ranging from a few km2 to ~500,000 km2. We discuss challenges for implementing alternative variables that are related to differences in the spatio-temporal resolution of observations and models, as well as to variable-specific features, for example, discrepancy between observed and simulated variables. We further discuss advancements required to deal with uncertainties of the hydrological data and to integrate multiple, potentially inconsistent datasets. The increased model consistency and improvement shown by most reviewed studies regarding the additional variables often come at the cost of a slight decrease in streamflow model performance. ...
Journal article (2025) - Magali Ponds, Sarah Hanus, Harry Zekollari, Marie-claire ten Veldhuis, Gerrit Schoups, Roland Kaitna, Markus Hrachowitz
Hydrological models play a vital role in projecting future changes in streamflow. Despite the strong awareness of non-stationarity in hydrological system characteristics, model parameters are typically assumed to be stationary and derived through calibration on past conditions. Integrating the dynamics of system change in hydrological models remains challenging due to uncertainties related to future changes in climate and ecosystems.

Nevertheless, there is increasing evidence that vegetation adjusts its root zone storage capacity – considered a critical parameter in hydrological models – to prevailing hydroclimatic conditions. This adaptation of the root zone to moisture deficits can be estimated by the Memory method. When combined with long-term water budget estimates from the Budyko framework, the Memory method offers a promising approach to estimate future climate–vegetation interaction and thus time-variable parameters in process-based hydrological models.

Our study provides an exploratory analysis of non-stationary parameters for root zone storage capacity in hydrological models for projecting streamflow in six catchments in the Austrian Alps, specifically investigating how future changes in root zone storage impact modeled streamflow. Using the Memory method, we derive climate-based parameter estimates of the root zone storage capacity under historical and projected future climate conditions. These climate-based estimates are then implemented in our hydrological model to assess the resultant impact on modeled past and future streamflow.

Our findings indicate that climate-based parameter estimations significantly narrow the parameter ranges linked to root zone storage capacity. This contrasts with the broader ranges obtained solely through calibration. Moreover, using projections from 14 climate models, our findings indicate a substantial increase in the root zone storage capacity parameters across all catchments in the future, ranging from +10 % to +100 %. Despite these alterations, the model performance remains relatively consistent when evaluating past streamflow, independent of using calibrated or climate-based estimations for the root zone storage capacity parameter. Additionally, no significant differences are found when modeling future streamflow when including future climate-induced adaptation of the root zone storage capacity in the hydrological model. Variations in annual mean, maximum and minimum flows remain within a 5 % range, with slight increases found for monthly streamflow and runoff coefficients. Our research shows that although climate-induced changes in root zone storage capacity occur, they do not notably affect future streamflow projections in the Alpine catchments under study. Our findings suggest that incorporating a dynamic representation of the root zone storage capacity parameter may not be crucial for modeling streamflow in humid and energy-limited catchments. However, our observations indicate relatively larger changes in root zone storage capacity within the less humid catchments, corresponding to higher variations in modeled future streamflow. This suggests a potentially higher importance of dynamic representations of root zone characteristics in arid regions and underscores the necessity for further research on non-stationarity in these regions. ...
Quantification of precipitation partitioning into evaporation and runoff is crucial for predicting future water availability. Within the widely used Budyko framework, which relates the long-term aridity index to the long-term evaporative index, curvilinear relationships between these indices (i.e. parametric Budyko curves) allow for the quantification of precipitation partitioning under prevailing climatic conditions. A common assumption is that movement along a specific Budyko curve with changes in the aridity index over time can be used as a predictor for catchment responses to changing climatic conditions. However, various studies have reported deviations around these curves, which raises questions about the usefulness of the method for future predictions. To investigate whether parametric Budyko curves still have predictive power, we quantified the global, regional, and local evolution of deviations of catchments from their parametric Budyko curves over multiple subsequent 20-year periods throughout the last century based on historical long-term water balance data from over 2000 river catchments worldwide. This process resulted in up to four 20-year distributions of annual deviations from the long-term mean parametric curve for each catchment. To use these distributions of deviations to predict future deviations, the temporal stability of these four distributions of deviations was evaluated between subsequent periods of time. On average, it was found that the majority (62 %) of study catchments did not significantly deviate from their expected parametric Budyko curves. Out of the remaining 38 % of catchments that deviated from their expected curves, the long-term magnitude of median deviations remains minor, with 70 % of catchments falling within the range of ±0.025 of the expected evaporative index. When these median deviations were expressed as relative changes in discharge, catchments in arid regions showed higher susceptibility to larger discharge shifts compared to those in humid regions. Furthermore, a significant majority of catchments, constituting around the same percentage, was found to have stable distributions of deviations across multiple time periods, making them well suited to statistically predict future deviations with high predictive power. These findings suggest that while trajectories of change in catchments do not strictly follow the expected long-term mean parametric Budyko curves, the deviations are minor and quantifiable. Consequently, taking into account these deviations, the parametric formulations of the Budyko framework remain a valuable tool for predicting future evaporation and runoff under changing climatic conditions within quantifiable margins of error. ...
Journal article (2025) - Wouter R. Berghuijs, Ross A. Woods, Bailey J. Anderson, Anna Luisa Hemshorn de Sanchez, Markus Hrachowitz
The water balance of catchments will, in many cases, strongly depend on its state in the recent past (e.g. previous days). Processes causing significant hydrological memory may persist at longer timescales (e.g. annual). The presence of such memory could prolong drought and flood risks and affect water resources over long periods, but the global universality, strength, and origin of long memory in the water cycle remain largely unclear. Here, we quantify annual memory in the terrestrial water cycle globally using autocorrelation applied to annual time series of water balance components. These time series of streamflow, global gridded precipitation, and GLEAM potential and actual evaporation, along with a GRACE (Gravity Recovery and Climate Experiment)-informed global terrestrial water storage reconstruction, indicate that, at annual timescales, memory is typically absent in precipitation but strong in terrestrial water stores (root zone moisture and groundwater). Outgoing fluxes (streamflow and evaporation) positively scale with storage, and so they also tend to hold substantial annual memory. As storage mediates flow extremes, such memory often also occurs in annual extreme flows and is especially strong for low flows and in large catchments. Our model experiments show that storage–discharge relationships that are hysteretic and strongly non-linear are consistent with these observed memory behaviours, whereas non-hysteretic and linear drainage fails to reconstruct these signals. Thus, a multi-year slow dance of terrestrial water stores and their outgoing fluxes is common; it is not simply mirroring precipitation memory and appears to be caused by hysteretic storage and drainage mechanisms that are incorporable in hydrological models. ...
Journal article (2025) - Mark C. Drever, Luis M. Bautista-Sopelana, Juan C. Alonso, Juan G. Navedo, Markus Hrachowitz
The migration of birds through a stopover site may be understood as a physical process of hydrological flow through a reservoir whose water levels rise and fall over a migration period. Hydrological flow models show promise as a way of integrating information on storage (daily counts of birds), inflow/outflow (number of birds arriving/departing each day) and transit time (length of stay (LOS)) over a migration period. We used a hydrological flow model to evaluate the relationship between stopover duration and passage population sizes of migrant birds under variable wetland conditions in two case studies. First, we considered the northward migration of Common Cranes Grus grus at Gallocanta Nature Reserve in Spain. We calibrated the model with daily counts recorded in 1984 and 1985, and then used it to predict transit time distributions based on counts of cranes 30 years later (2015–17). The model was calibrated with a mean transit time of 6.5 days observed in 1984/85, consistent with observed values of LOS, and predicted a mean transit time of 5.2 days for the 2015–17 period. The model also predicted an ~6.2× increase of the total migratory passage population of cranes at Gallocanta, which qualitatively agreed with the large increase in the overall population from the 1980s to 2010s. Second, migration dynamics of Eurasian Spoonbills Platalea leucorodia were considered during southward migration at two intertidal coastal wetlands in northern Spain from 2002 to 2005. The model well captured the observed differences in transit time between Urdaibai Biosphere Reserve (median transit time of 1.1 days) and Santoña Marshes Natural Park (median transit time of 2.1 days). Transit times of both species were negatively correlated with estimates of annual population sizes. These results indicate that hydrological flow models can provide insights into the migration ecology of waterbirds (or species where regular counts during migration are available), and that LOS is a dynamic decision that can depend on the population sizes of migratory birds. ...
The temporal dynamics of water ages provide crucial insights into hydrological processes and transport mechanisms, yet there remains a significant gap in quantifying water age variability across different temporal scales. This study utilizes a comprehensive dataset spanning 70 years of hydrological observations and tritium records (1953–2022) with a semi-distributed hydrological model with integrated tracer routing routine based on StorageAge Selection functions SAS, to explore the temporal evolution of water ages in the 4000 km2 Upper Neckar River basin, Germany. Our findings indicate a systematic convergence of the variability of young water fractions and other metrics of water age in riverflow and evaporation towards stable values when averaging over increasing time scales. While at daily scales exhibiting considerable variability with young water fractions in riverflow Fwy,Q ∼ 0.01–0.91 and in evaporation Fwy,E ∼ 0.02–0.75, the variability of Fwy,Q and Fwy,E gradually reduces with increasing averaging time scales and converge to 0.45–0.47 and 0.96–0.97, respectively, between individual decades. Liquid water input (PL), comprising rainfall and snow melt, emerges as the dominant driver of Fwy,Q across all time scales. In contrast, Fwy,E shows varying controls with time scale: soil moisture content governs daily fluctuations, whereas PL dominates at the decadal scale. Overall, water ages demonstrate remarkable stability with only minor deviations in response to climatic variability: a 20% fluctuation in average decadal PL results in only ∼4% variation in Fwy,Q and ∼1% in Fwy,E over the study period. These findings suggest a lack of major long-term dynamics in water ages. Consequently, the results suggest that the physical transport dynamics in the Upper Neckar River basin, and potentially in comparable river basins with similar water age characteristics, can be considered near-stationary over multiple decades. ...
Journal article (2025) - Jordy Salmon-Monviola, Ophélie Fovet, Markus Hrachowitz
The consistency of hydrological models, i.e. their ability to reproduce observed system dynamics, needs to be improved to increase their predictive power. As using streamflow data alone to calibrate models is not sufficient to constrain them and render them consistent, other strategies must be considered, in particular using additional types of data. The aim of this study was to test whether simultaneous calibration of dissolved organic carbon (DOC) and nitrate (NO3-) concentrations along with streamflow improved the hydrological consistency of a parsimonious solute-transport model. A multi-objective approach with four calibration scenarios was used to evaluate the model's predictions for an intensive agricultural headwater catchment. After calibration, the model reasonably simultaneously reproduced the dynamics of discharge and DOC and NO3- concentrations in the stream of the headwater catchment from 2008-2016. Evaluation using independent datasets indicated that the model usually reproduced dynamics of groundwater level and soil moisture in upslope and riparian zones correctly for all calibration scenarios. Using daily stream concentrations of DOC and NO3- along with streamflow to calibrate the model did not improve its ability to predict streamflow for calibration or evaluation periods. The approach significantly improved the representation of groundwater storage and to a lesser extent soil moisture in the upslope zone but not in the riparian zone. Parameter uncertainty decreased when the model was calibrated using solute concentrations, except for parameters related to fast and slow reservoir flow. This study shows the added value of using multiple types of data along with streamflow, in particular DOC and NO3- concentrations, to constrain hydrological models to improve representation of internal hydrological states and flows. With the increasing availability of solute data from catchment monitoring, this approach provides an objective way to improve the consistency of hydrological models that can be used with confidence to evaluate scenarios. ...
Journal article (2024) - Berit Arheimer, Christophe Cudennec, Attilio Castellarin, Soham Adla, Judith Boekee, Thom Bogaard, Markus Hrachowitz, Saket Pande, Hubert Savenije, More Authors...
The new scientific decade (2023-2032) of the International Association of Hydrological Sciences (IAHS) aims at searching for sustainable solutions to undesired water conditions–whether it be too little, too much or too polluted. Many of the current issues originate from global change, while solutions to problems must embrace local understanding and context. The decade will explore the current water crises by searching for actionable knowledge within three themes: global and local interactions, sustainable solutions and innovative cross-cutting methods. We capitalise on previous IAHS Scientific Decades shaping a trilogy; from Hydrological Predictions (PUB) to Change and Interdisciplinarity (Panta Rhei) to Solutions (HELPING). The vision is to solve fundamental water-related environmental and societal problems by engaging with other disciplines and local stakeholders. The decade endorses mutual learning and co-creation to progress towards UN sustainable development goals. Hence, HELPING is a vehicle for putting science in action, driven by scientists working on local hydrology in coordination with local, regional, and global processes. ...
Journal article (2024) - J.M. van den Brink, H. Scharff, B. Steinert, S. Melchior, M. Hrachowitz, J. Heimovaara, J. Gebert
This study quantifies the field hydraulic performance of a dual-functionality landfill cover, combining microbial methane oxidation with water diversion using a capillary barrier. The investigated 500 m2 test field, constructed on a landfill in the Netherlands, consisted of a cover soil optimised for methane oxidation, underlain by a sandy capillary layer and a gravelly capillary block. Outflows from these layers were measured between 2009 and 2023. Average precipitation was 848 mm/a, evapotranspiration, diverted infiltration and breakthrough amounted to 504 (59.4 %), 282 (33.3 %) and 62 (7.3 %) mm/a, respectively. On average, the capillary barrier diverted 82 % of the inflow into the capillary layer. Breakthrough occurred mainly from October to March when evapotranspiration was low and the maximum water storage capacity of the cover soil was reached. During this period, inflow into the capillary barrier exceeded its diversion capacity, caused by the relatively high hydraulic conductivity of the cover soil due to its optimisation for gas transport. The diversion capacity declined drastically in the year after construction and increased again afterwards. This was attributed to suffusion of sand from the capillary layer into the capillary block and subsequent washout to greater depths or the influence of iron precipitates at the bottom of the capillary layer. The effect of a more finely grained methane oxidation layer on the hydraulic and methane oxidation performance should be investigated further. These measures could further improve the combined performance of the dual functionality landfill cover system under the given conditions of a temperate climate. ...

Application to Hoek van Holland (Netherlands)

Coastal flooding events pose a critical risk in delta areas, since they are characterized by population growth and urban expansion. A better understanding of Extreme Water Levels (EWLs), the mechanisms generating them, and their components, i.e., astronomical tide and storm surge is of great importance as they drive the maintenance and design of flood protection systems. Therefore, a statistical investigation of them can provide new insights for more reliable flood risk mitigation infrastructures. In this study, we analyse these components and compare different probabilistic methods i.e., univariate extreme value analysis, copula functions, and Joint Probability Method (JPM) for the better estimation of EWLs. We use Hoek van Holland (NL) as a representative case study, since the dynamic conditions of this deltaic environment with man-made infrastructures render the area of strategic importance. The results indicate that a more accurate estimate of the declustering time between extreme events can be achieved using correlation of high surges and high wind speeds, taking into consideration also the wind direction. In the Southwest Delta this time estimated to be around 4 days. Furthermore, the EWLs components, i.e., surge and astronomical tide, show negative dependence. From the comparison between statistical approaches to model EWLs, results show that EWLs estimated via EVA and JPM do not vary significantly, while copulas’ seems to outperform the other methods. However, the selection of the proper copula to show the dependence is critical. As a conclusion, the analysis of the dependence between tides and storm surges can lead to more robust inferences of EWLs. ...

An integrated dataset and catalogue of streamflow, hydro-climatic and landscape variables for Europe

Journal article (2024) - Thiago V.M. do Nascimento, Julia Rudlang, Marvin Höge, Ruud van der Ent, Máté Chappon, Jan Seibert, Markus Hrachowitz, Fabrizio Fenicia
Large-sample hydrology datasets have become increasingly available, contributing to significant scientific advances. However, in Europe, only a few such datasets have been published, capturing only a fraction of the wealth of information from national data providers in terms of available spatial density and temporal extent. We present “EStreams”, an extensive dataset of hydro-climatic variables and landscape descriptors and a catalogue of openly available stream records for 17,130 European catchments. Spanning up to 120 years, the dataset includes streamflow indices, catchment-aggregated hydro-climatic signatures and landscape attributes (topography, soils, geology, vegetation and landcover). The catalogue provides detailed descriptions that allow users to directly access streamflow data sources, overcoming challenges related to data redistribution policies, language barriers and varied data portal structures. EStreams also provides Python scripts for data retrieval, aggregation and processing, making it dynamic in contrast to static datasets. This approach enables users to update their data as new records become available. Our goal is to extend current large-sample datasets and further integrate hydro-climatic and landscape data across Europe. ...