Jingjing Huang, Hang Lyu, Chao Wang, Xiaosi Su, Weihong Dong, Yuyu Wan, Teijun Song, Xiaofang Shen
In this study, we developed a 2.5-D microfluidic experimental platform that enables nondestructive visualization of pore-throat microstructures and the migration processes of residual non-aqueous phase liquids (NAPL) within a three-phase system (water–ice–NAPL) under freeze–thaw conditions. Simulated experiments were then used to investigate the freeze–thaw migration of typical NAPL components. We analyzed the forces within the system during freezing and proposed an expulsion factor (<span data-altimg="/cms/asset/9472be0b-6031-4ca2-99a3-4f3073e667af/wrcr70777-math-0001.png"></span><mjx-container ctxtmenu_counter="145" ctxtmenu_oldtabindex="1" jax="CHTML" role="application" sre-explorer- style="font-size: 103%; position: relative;" tabindex="0"><mjx-math aria-hidden="true" location="graphic/wrcr70777-math-0001.png"><mjx-semantics><mjx-mrow><mjx-msub data-semantic-children="0,1" data-semantic- data-semantic-role="latinletter" data-semantic-speech="upper N Subscript f" data-semantic-type="subscript"><mjx-mi data-semantic-annotation="clearspeak:simple" data-semantic-font="italic" data-semantic- data-semantic-parent="2" data-semantic-role="latinletter" data-semantic-type="identifier"><mjx-c></mjx-c></mjx-mi><mjx-script style="vertical-align: -0.15em; margin-left: -0.085em;"><mjx-mi data-semantic-annotation="clearspeak:simple" data-semantic-font="italic" data-semantic- data-semantic-parent="2" data-semantic-role="latinletter" data-semantic-type="identifier" size="s"><mjx-c></mjx-c></mjx-mi></mjx-script></mjx-msub></mjx-mrow></mjx-semantics></mjx-math><mjx-assistive-mml display="inline" unselectable="on"><math altimg="urn:x-wiley:00431397:media:wrcr70777:wrcr70777-math-0001" display="inline" location="graphic/wrcr70777-math-0001.png" xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><msub data-semantic-="" data-semantic-children="0,1" data-semantic-role="latinletter" data-semantic-speech="upper N Subscript f" data-semantic-type="subscript"><mi data-semantic-="" data-semantic-annotation="clearspeak:simple" data-semantic-font="italic" data-semantic-parent="2" data-semantic-role="latinletter" data-semantic-type="identifier">N</mi><mi data-semantic-="" data-semantic-annotation="clearspeak:simple" data-semantic-font="italic" data-semantic-parent="2" data-semantic-role="latinletter" data-semantic-type="identifier">f</mi></msub></mrow>${N}_{f}$</annotation></semantics></math></mjx-assistive-mml></mjx-container>) to evaluate the migration potential of residual NAPL. The results indicated that the residual NAPL exhibited three primary mechanisms during freezing: expulsion, snap-off, and compression. Expulsion significantly drove the residual NAPL movement during the early stages with fewer freeze–thaw cycles. As the number of cycles increased, snap-off and compression became dominant. This substantially restricted the expulsion of the residual NAPL. These findings have significant theoretical and practical value to understand and predict contaminant mo
{"title":"Residual NAPL Intermittent Expulsion With Increasing Freeze–Thaw Cycles in Saturated Porous Media: Findings From a 2.5-D Microfluidic Platform","authors":"Jingjing Huang, Hang Lyu, Chao Wang, Xiaosi Su, Weihong Dong, Yuyu Wan, Teijun Song, Xiaofang Shen","doi":"10.1029/2025wr042230","DOIUrl":"https://doi.org/10.1029/2025wr042230","url":null,"abstract":"In this study, we developed a 2.5-D microfluidic experimental platform that enables nondestructive visualization of pore-throat microstructures and the migration processes of residual non-aqueous phase liquids (NAPL) within a three-phase system (water–ice–NAPL) under freeze–thaw conditions. Simulated experiments were then used to investigate the freeze–thaw migration of typical NAPL components. We analyzed the forces within the system during freezing and proposed an expulsion factor (<span data-altimg=\"/cms/asset/9472be0b-6031-4ca2-99a3-4f3073e667af/wrcr70777-math-0001.png\"></span><mjx-container ctxtmenu_counter=\"145\" ctxtmenu_oldtabindex=\"1\" jax=\"CHTML\" role=\"application\" sre-explorer- style=\"font-size: 103%; position: relative;\" tabindex=\"0\"><mjx-math aria-hidden=\"true\" location=\"graphic/wrcr70777-math-0001.png\"><mjx-semantics><mjx-mrow><mjx-msub data-semantic-children=\"0,1\" data-semantic- data-semantic-role=\"latinletter\" data-semantic-speech=\"upper N Subscript f\" data-semantic-type=\"subscript\"><mjx-mi data-semantic-annotation=\"clearspeak:simple\" data-semantic-font=\"italic\" data-semantic- data-semantic-parent=\"2\" data-semantic-role=\"latinletter\" data-semantic-type=\"identifier\"><mjx-c></mjx-c></mjx-mi><mjx-script style=\"vertical-align: -0.15em; margin-left: -0.085em;\"><mjx-mi data-semantic-annotation=\"clearspeak:simple\" data-semantic-font=\"italic\" data-semantic- data-semantic-parent=\"2\" data-semantic-role=\"latinletter\" data-semantic-type=\"identifier\" size=\"s\"><mjx-c></mjx-c></mjx-mi></mjx-script></mjx-msub></mjx-mrow></mjx-semantics></mjx-math><mjx-assistive-mml display=\"inline\" unselectable=\"on\"><math altimg=\"urn:x-wiley:00431397:media:wrcr70777:wrcr70777-math-0001\" display=\"inline\" location=\"graphic/wrcr70777-math-0001.png\" xmlns=\"http://www.w3.org/1998/Math/MathML\"><semantics><mrow><msub data-semantic-=\"\" data-semantic-children=\"0,1\" data-semantic-role=\"latinletter\" data-semantic-speech=\"upper N Subscript f\" data-semantic-type=\"subscript\"><mi data-semantic-=\"\" data-semantic-annotation=\"clearspeak:simple\" data-semantic-font=\"italic\" data-semantic-parent=\"2\" data-semantic-role=\"latinletter\" data-semantic-type=\"identifier\">N</mi><mi data-semantic-=\"\" data-semantic-annotation=\"clearspeak:simple\" data-semantic-font=\"italic\" data-semantic-parent=\"2\" data-semantic-role=\"latinletter\" data-semantic-type=\"identifier\">f</mi></msub></mrow>${N}_{f}$</annotation></semantics></math></mjx-assistive-mml></mjx-container>) to evaluate the migration potential of residual NAPL. The results indicated that the residual NAPL exhibited three primary mechanisms during freezing: expulsion, snap-off, and compression. Expulsion significantly drove the residual NAPL movement during the early stages with fewer freeze–thaw cycles. As the number of cycles increased, snap-off and compression became dominant. This substantially restricted the expulsion of the residual NAPL. These findings have significant theoretical and practical value to understand and predict contaminant mo","PeriodicalId":23799,"journal":{"name":"Water Resources Research","volume":"53 1","pages":""},"PeriodicalIF":5.4,"publicationDate":"2026-03-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147361026","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Döndü Sarışen, Fayyaz Ali Memon, Hossein Mohammadi, Raziyeh Farmani
Many researchers have used deterministic models to address issues in intermittent water supply (IWS) systems, such as inequitable distribution, but these models often overlook key uncertainties. IWS systems vary widely in practices and operations, introducing uncertainties in areas like water consumption pattern, supply characteristics, and household tank sizes. To address these challenges, this study proposes a novel framework for assessing uncertainty in model input parameters and applies it to an IWS network using an EPA-SWMM-based hydraulic simulation. In the first phase, uncertainty analysis (UA) uses probabilistic methods and Monte Carlo simulations to evaluate the impact of uncertainties on performance indicators. The second phase applies global sensitivity analysis (SA) (Sobol's method) to identify the most influential parameters. The findings reveal that system performance is primarily governed by supply characteristics, while household tank size exerts a secondary but nonlinear influence on both pressure and consumption based indicators. Excessive household tank size is shown to reduce the pressure and supply equity, whereas moderate tank size improves the fairness of the water distribution. This work provides a foundation for more accurate and reliable IWS modeling approaches.
{"title":"Analyzing and Quantifying Key Sources of Uncertainty in Intermittent Water Supply Simulation: Supply Characteristics, Household Tank Size and Time Series of Water Consumption From the Household Tanks","authors":"Döndü Sarışen, Fayyaz Ali Memon, Hossein Mohammadi, Raziyeh Farmani","doi":"10.1029/2024wr039282","DOIUrl":"https://doi.org/10.1029/2024wr039282","url":null,"abstract":"Many researchers have used deterministic models to address issues in intermittent water supply (IWS) systems, such as inequitable distribution, but these models often overlook key uncertainties. IWS systems vary widely in practices and operations, introducing uncertainties in areas like water consumption pattern, supply characteristics, and household tank sizes. To address these challenges, this study proposes a novel framework for assessing uncertainty in model input parameters and applies it to an IWS network using an EPA-SWMM-based hydraulic simulation. In the first phase, uncertainty analysis (UA) uses probabilistic methods and Monte Carlo simulations to evaluate the impact of uncertainties on performance indicators. The second phase applies global sensitivity analysis (SA) (Sobol's method) to identify the most influential parameters. The findings reveal that system performance is primarily governed by supply characteristics, while household tank size exerts a secondary but nonlinear influence on both pressure and consumption based indicators. Excessive household tank size is shown to reduce the pressure and supply equity, whereas moderate tank size improves the fairness of the water distribution. This work provides a foundation for more accurate and reliable IWS modeling approaches.","PeriodicalId":23799,"journal":{"name":"Water Resources Research","volume":"14 1","pages":""},"PeriodicalIF":5.4,"publicationDate":"2026-03-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147360150","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Veerle C. Bril, Jens de Bruijn, Hans de Moel, Tarun Sadana, Tim Busker, W. J. Wouter Botzen, Jeroen C. J. H. Aerts
Floods are expected to increase in frequency and severity due to climate change. Recent floods have shown that many catchments worldwide are vulnerable to floods, highlighting the need for additional adaptation measures. This study extends the Geographical, Environmental, and Behavioral (GEB) model by coupling it to a hydrodynamic and a flood risk model to assess the effects of dry-proofing, wet-proofing, retention ponds, reforestation, and the creation of natural grassland. A key innovation is the integration of all local-scale models, thereby allowing for a catchment-wide assessment of the impacts of various measures on interlinked hydrological conditions, flood extents and depths, damages, and risk. We apply our method to the Geul catchment (shared between the Netherlands, Belgium and Germany), which was heavily flooded in July 2021. Our results show that reforestation and creation of natural grassland (both 10 km2) reduce flood extent by 12% and average water depth by 10%. Damage is decreased up to 38%. Larger retention ponds (1 m deeper) have a much smaller reduction in flood extent (3%), depth (0.5%) and damage (1.6%), due to limited storage capacity compared to excess rainfall. The building-level adaptation scenarios outperform all nature-based solutions, with dry-proofing reducing more damage (up to 95%) than wet-proofing (around 55%). A cost-benefit analysis shows that several adaptation measures are economically attractive. Overall, our findings show a coupled model is essential for comparing the relative effectiveness of different flood adaptation measures and supporting informed risk management decisions. The open-source model is transferable to other catchments worldwide to guide decision-making.
{"title":"Assessing the Effectiveness of Nature-Based Solutions and Building-Level Flood Risk Reduction Measures: An Open-Source Coupled Model","authors":"Veerle C. Bril, Jens de Bruijn, Hans de Moel, Tarun Sadana, Tim Busker, W. J. Wouter Botzen, Jeroen C. J. H. Aerts","doi":"10.1029/2025wr041436","DOIUrl":"https://doi.org/10.1029/2025wr041436","url":null,"abstract":"Floods are expected to increase in frequency and severity due to climate change. Recent floods have shown that many catchments worldwide are vulnerable to floods, highlighting the need for additional adaptation measures. This study extends the Geographical, Environmental, and Behavioral (GEB) model by coupling it to a hydrodynamic and a flood risk model to assess the effects of dry-proofing, wet-proofing, retention ponds, reforestation, and the creation of natural grassland. A key innovation is the integration of all local-scale models, thereby allowing for a catchment-wide assessment of the impacts of various measures on interlinked hydrological conditions, flood extents and depths, damages, and risk. We apply our method to the Geul catchment (shared between the Netherlands, Belgium and Germany), which was heavily flooded in July 2021. Our results show that reforestation and creation of natural grassland (both 10 km<sup>2</sup>) reduce flood extent by 12% and average water depth by 10%. Damage is decreased up to 38%. Larger retention ponds (1 m deeper) have a much smaller reduction in flood extent (3%), depth (0.5%) and damage (1.6%), due to limited storage capacity compared to excess rainfall. The building-level adaptation scenarios outperform all nature-based solutions, with dry-proofing reducing more damage (up to 95%) than wet-proofing (around 55%). A cost-benefit analysis shows that several adaptation measures are economically attractive. Overall, our findings show a coupled model is essential for comparing the relative effectiveness of different flood adaptation measures and supporting informed risk management decisions. The open-source model is transferable to other catchments worldwide to guide decision-making.","PeriodicalId":23799,"journal":{"name":"Water Resources Research","volume":"5 1","pages":""},"PeriodicalIF":5.4,"publicationDate":"2026-03-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147361027","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Daniel J. Short Gianotti, Meriah J. Gannon, Dara Entekhabi
Meteorological droughts (persistent precipitation deficits) often, but not always, transition into agricultural droughts (persistent soil moisture deficits). The intensity of agricultural drought, however, can vary for a given precipitation deficit due to a number of catalyzing co-factors beyond precipitation such as atmospheric evaporative demand and temperature. In this study we use Earth System Model data to quantify (a) how warm temperature anomalies affect this evolution from meteorological-to-agricultural drought and (b) how the evolution of droughts from historical and future climate scenarios differ. We benchmark these results against observational data and use a multi-model ensemble to quantify agreement on future drought propagation. Broadly speaking, drought temperatures in the upper third of local distributions correspond with shifts on the order of 5 percentile of the soil moisture distribution. We would expect today's meteorological droughts to propagate into agricultural droughts roughly one drought classification more severe in the SSP3-7.0 scenario in most regions. Even regions with increases in precipitation are likely to see more intense meteorological-to-agricultural drought propagation by the end of the 21st century. Models disagree on drought propagation changes in Africa for the same precipitation deficit, but suggest that all historical droughts would have had worse agricultural droughts in Europe and Eastern North America if they happened under SSP3-7.0. When accounting for precipitation changes—which tend toward more frequent accumulated precipitation deficits—the increased severity of meteorological-to-agricultural drought evolution leads to predictions of major increases in moderate to extreme (D1–D3) drought events in all regions globally by the end of the century.
{"title":"Meteorological to Agricultural Drought Transitions Compounded by Heat Waves in Historical and Future Climates","authors":"Daniel J. Short Gianotti, Meriah J. Gannon, Dara Entekhabi","doi":"10.1029/2024wr039170","DOIUrl":"https://doi.org/10.1029/2024wr039170","url":null,"abstract":"Meteorological droughts (persistent precipitation deficits) often, but not always, transition into agricultural droughts (persistent soil moisture deficits). The intensity of agricultural drought, however, can vary for a given precipitation deficit due to a number of catalyzing co-factors beyond precipitation such as atmospheric evaporative demand and temperature. In this study we use Earth System Model data to quantify (a) how warm temperature anomalies affect this evolution from meteorological-to-agricultural drought and (b) how the evolution of droughts from historical and future climate scenarios differ. We benchmark these results against observational data and use a multi-model ensemble to quantify agreement on future drought propagation. Broadly speaking, drought temperatures in the upper third of local distributions correspond with shifts on the order of 5 percentile of the soil moisture distribution. We would expect today's meteorological droughts to propagate into agricultural droughts roughly one drought classification more severe in the SSP3-7.0 scenario in most regions. Even regions with increases in precipitation are likely to see more intense meteorological-to-agricultural drought propagation by the end of the 21st century. Models disagree on drought propagation changes in Africa for the same precipitation deficit, but suggest that all historical droughts would have had worse agricultural droughts in Europe and Eastern North America if they happened under SSP3-7.0. When accounting for precipitation changes—which tend toward more frequent accumulated precipitation deficits—the increased severity of meteorological-to-agricultural drought evolution leads to predictions of major increases in moderate to extreme (D1–D3) drought events in all regions globally by the end of the century.","PeriodicalId":23799,"journal":{"name":"Water Resources Research","volume":"48 1","pages":""},"PeriodicalIF":5.4,"publicationDate":"2026-03-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147360151","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The filling process of large hydropower reservoirs is an often-overlooked process that can attract most of the attention and criticisms of a dam's lifetime. Understanding what drove past filling efforts and exploring trade-off solutions is therefore critical for regions where hydropower development is flourishing. Yet, there is a lack of methodological frameworks applicable across different contexts (e.g., existing and planned systems) that can characterize how trade-offs evolve under varying hydro-climatological conditions. Here, we focus on the filling strategies of the Upper Mekong hydropower system, one of the largest mega-dam cascades (capacity