Flood is one of the most catastrophic natural disasters in the United States, particularly in the Southeast states where hurricanes and tropical storms are most prevalent, causing billions of dollars in damage annually and significant losses of life and property. The Weather Research and Forecasting Hydrological model (WRF-Hydro) is a community-based hydrologic model designed to improve the skill of

Safe and cost-effective design of infrastructures, such as dams, bridges, highways, often requires knowing the magnitude and frequency of peak floods. The Generalized Extreme Value distribution (GEV) prevailed in flood frequency analysis along with distributions comprising location, scale, and shape parameters. Here we explore alternative models and propose power-type models, having one scale and two

Urban pluvial flooding is a threatening natural hazard in urban areas all over the world, especially in recent years given its increasing frequency of occurrence. In order to prevent flood occurrence and mitigate the subsequent aftermath, urban water managers aim to predict precipitation characteristics, including peak intensity, arrival time and duration, so that they can further warn inhabitants

A procedure combining experiments and deep learning is demonstrated to acquire pore-scale images of oil- and water-wet surfaces over a large field of view in microfluidic devices and to classify wettability based upon these pore scale images. Deep learning supplants the manual, time-consuming, error-prone investigation and categorization of such images. Image datasets were obtained by visualizing the

The morphodynamics of tidal marshes is investigated using a novel model (MarshMorpho2D) that couples 2D tidal dynamics, wave dynamics, mass-conserving sediment transport, edge erosion, pond dynamics, and marsh upland migration. By only changing tidal range, sediment supply, and RSLR rate, the model reproduces a great variety of spatial patterns in marsh loss. For all tidal ranges, an increase in RSLR

Estimation of the permeability tensor is crucial for many natural and industrial porous media applications. Homogenization provides a rigorous framework to calculate effective parameters from pore-scale images of a representative unit cell of the porous medium by solving a boundary value problem, also known as closure problem, subject to global constraints. However, the latter are hard to satisfy for

In this article we consider a two-phase flow model in a highly heterogeneous porous column. The porous column consists of homogeneous blocks, where the porosity and permeability vary from one block to the other. The flow direction is perpendicular to the layering of the porous column, and hence can be approximated by one-dimensional model equations. The periodic change in porosity and absolute permeability

A series of micro-scale core-flooding experiments were conducted in a miniature water-wet Berea sandstone sample to examine how changes in the desaturation (water-displacing-oil) protocol impact the non-wetting phase recovery and to provide, to the best of our knowledge, the first pore-scale evidences of the associated displacement subtleties. We studied (i) the continuous displacement mode (CDM),

In this work we develop a multiscale model for flow and transport problem in CO2-enhanced coalbed methane recovery. The coalbed methane reservoir is characterized by two levels of porosity associated with nanopores in the matrix and cleat network. Mass conservation equations for fluid mixture (CH4 and CO2) in the matrix at the microscale are rigorously derived by using the formal homogenization technique

The effect of emergent stems on the transport of downstream floating particles (e.g., buoyant seeds) is explored theoretically and experimentally at moderate to high Reynolds number (Rd=2Ubas/υ>300) in an open channel, where Ub is the bulk velocity, as is stem radius, and ν is the kinematic viscosity. Longitudinal dispersion (Dl) of such seeds is shown to be given by Dl=Up3τ02η(1−η)/2/(S1+ητ0Up) thereby

Hydraulic barriers using either two or three injection wells are studied analytically to evaluate their performance on improving groundwater extraction in a coastal aquifer under the threat of seawater intrusion. The multiple well system is operated by extracting fresh groundwater through a landward well and partially reinjecting into the aquifer through two or three wells located parallel to the coastline

We explore the coupling of surface and subsurface flows on fully unstructured meshes that conform to complex soil structures. To accommodate the distorted meshes that inevitably result from explicit representation of complex soil structures, we leverage the structure of the Mimetic Finite Difference (MFD) spatial discretization scheme to couple surface and subsurface flows. The MFD method achieves

Solving image-based pore-scale flow and transport in porous materials is a mainstream of fundamental and industrial research to reveal the pertinent physics in the field of hydrogeology, reservoir engineering, paper and filter engineering. This research discipline requires tremendous integration of multidisciplinary research areas of image processing, computational fluid modelling, and high-performance

Riverine bathymetry is of crucial importance for shipping operations and flood management. However, obtaining direct measurements of depth is not always easy. Conversely, with recent advances in sensor technology, indirect measurements can be obtained and used to estimate high-resolution river bed topography. Physics-based inverse modeling techniques have been used to estimate bathymetry using indirect

We introduce original sensitivity analysis metrics with the aim of assisting diagnosis of the functioning of a given model. We do so by characterizing model-induced dependencies between a target model output and selected model input(s) through the associated bivariate copula density. The latter fully characterizes the dependencies between two random variables at any order (i.e., without being limited

We present empirical evidence of the impact of water markets on production. From 1954 to 2012, we compare the crop composition in counties of the Rio Grande water market with neighboring control counties before and after the establishment of the water market in 1971. What we find is consistent with water markets facilitating a shift from low- to high-value crops, or, in terms of our empirical proxies

In coastal areas, geological heterogeneity influences the geometry and characteristics of the freshwater-saltwater mixing zone. Geophysical electrical and electromagnetic methods are increasingly used to delineate the freshwater-saltwater mixing zone or for groundwater model calibration. However, in practical applications, it is common when assessing pore water salinity from resistivity models to disregard

Groundwater, being a vital component of the natural water resource system, needs continuous monitoring and dynamic management strategies. That said, we require computationally inexpensive groundwater flow models for repetitive solutions with desirable accuracy under budgetary limitation(s). Natural aquifer systems inherit strong heterogeneity at local scales. In this work, we have proposed ordinary

In this study, we have developed a combined lattice Boltzmann-smoothed profile method to explore coupled mechanisms governing transport of colloids and their retention in porous media. We have considered flow in a constricted tube and included hydrodynamic, gravity, buoyancy, van der Waals and electrostatic forces to simulate colloid transport and aggregation. A major advantage of this complete formulation

Classic extreme value theory provides asymptotic distributions of block maxima (BM) Y or peaks over threshold (POT). However, as BM and POT are relatively small subsets of the values assumed by the parent process Z, alternative approaches have been proposed to make inferences on extremes by using intermediate values and non-asymptotic models. In this study, we investigate the finite sample theory of

This article studies the individual contribution of each term in the Navier-Stokes-Darcy (Brinkman) equation to the fluid linear momentum balance and the consequences of opting for the Darcy or Brinkman equations in 3D reactive flow simulations of a vuggy carbonate. A two-scale continuum model was used to simulate the acidification of a carbonate sample using the Darcy and Brinkman equations. Rock

Multiwavelets (MW) enable the compression, analysis and assembly of model data on a multiresolution grid within Godunov-type solvers based on second-order discontinuous Galerkin (DG2) and first-order finite volume (FV1) methods. Multiwavelet adaptivity has been studied extensively with one-dimensional (1D) hydrodynamic models (Kesserwani et al., 2019), revealing that MWDG2 can be 20 times faster than

The impact of phases distribution on mixing and reaction is hardly assessable experimentally. We use a multiple point statistical method, which belongs to the family of machine learning algorithms, to generate simulations of phases distributions from data out of laboratory experiments. The simulations honour the saturation of the laboratory experiments, resemble the statistical distributions of several

We propose a thermodynamic approach to modeling unsaturated flow in porous media, where the liquid saturation is understood as the state variable. The free energy functional is designed as a symmetric expansion of the traditional capillary energy density in Richards equation, therefore removing ambiguities on the interpretation of the higher-order term in the model equation. The proposed definition

Literature shows the application of the entropy concept in open channels flow. No existing research has been addressed to the application of the entropy concept for evaluating the effect of secondary flow which especially develops in curved channels. The purpose of this paper is to explore the possibility to evaluate the effect of the secondary flow on the stream-wise velocity distribution starting

The lattice Boltzmann method (LBM) is a widely-used numerical technique for simulation of single- and two-phase flow in geometries that are obtained using tomographic imaging of natural porous media. Due to ease of implementation and numerical robustness, a vast majority of LBM-based pore-scale simulations employ the so-called bounceback scheme to enforce no-slip velocity boundary condition. Bounceback

Assessments of groundwater and surface water budgets at a large scale, such as the contiguous United States, often separately analyze the complex dynamics linking the surface and subsurface categories of water resources. These dynamics include recharge and groundwater contributions to streamflow. The time-varying simulation of these complex hydrologic dynamics, across large spatial and temporal scales

The Gravity Recovery and Climate Experiment (GRACE) satellites provide unprecedented perspectives to hydrologists and geoscientists for observing and understanding the variation of terrestrial water storage (TWS) at continental to global scales. However, there are few reliable datasets of past TWS variations before GRACE observations were available (pre-2002). To fill this gap, we attempt to develop

Richards equation for variably saturated flow can exhibit stability problems due to its nonlinearity. The challenge is to resolve sharp wetting fronts without introducing spurious oscillations, especially for simulations with very dry initial conditions. Flow instabilities at sharp wetting fronts may arise particularly in anisotropic and heterogeneous media, leading to oscillations or convergence problems

The study of fluid flow in non-Darcy porous medium subjected to Hall current is very important in scientific and engineering applications such as water filters, groundwater discharge in aquifers, petroleum engineering, MHD generators and chemical engineering. Forchheimer model is needed at high flow rate, where, flow will exhibit non-linearity with respect to velocity which make Darcy law inapplicable

Environmental measurements generate great volumes of high-dimensional data (often noisy and with missing values) from which meaningful messages may be extracted through appropriate organisation and summarisation. The self-organizing map (SOM) is an artificial neural network popular for recognizing patterns, relationships and clusters in such data. Through the finetuning of configuration and training

Accurate simulation of multiphase flow in subsurface formations is challenging, as the formations span large length scales (km) with high-resolution heterogeneous properties. To deal with this challenge, different multiscale methods have been developed. Such methods construct coarse-scale systems, based on a given high-resolution fine-scale system. Furthermore, they are amenable to parallel computing

Agriculture is the largest water user and the driving force behind global water stress. While several studies have already assessed global green and blue water use in crop production at high spatial resolution, few studies have explicitly linked water scarcity to water use in the production of individual crops and assessed the sustainability of each crop production. We provide a detailed assessment

Vegetation is a key source of flow resistance in natural channels and floodplains. It is therefore important to accurately model the flow resistance to inform decision makers and managers. However, it is challenging to predict the resistance of real vegetation, because vegetation models are based on relatively small-scale lab experiments with mostly artificial vegetation. Experimental tests of real

Uncertainties related to permeability heterogeneity can be estimated using geostatistical simulation methods. Usually, these methods are applied on regular grids with cells of constant size, whereas unstructured grids are more flexible to honor geological structures and offer local refinements for fluid-flow simulations. However, cells of different sizes require to account for the support dependency

There is increasing attention to the development of a myriad of complex methods for nonstationary frequency analysis (NFA) of floods, droughts and other hydrologic processes. We assume that the need for NFA arises from well understood deterministic mechanisms of change. A common assumption in NFA, questioned here, is that more accurate estimators of hydrologic statistics result when more realistic

Carbonates can be highly heterogeneous formations with large variations in pore size distribution and pore space topology, which results in complex multiphase flow behavior. Here we investigate the spontaneous imbibition behavior of fluids in vuggy carbonates. Glass beads of 1.0 mm diameter, with dissolvable vug placeholders, are sintered to form multiple configurations of heterogeneous vuggy core

Unforeseen clogging and low well productivity during exploitation of oil and gas reservoirs result from a mobilisation of densely deposited particles. The proposed model predicts the mobilisation in unsaturated porous media based on the changes in water level, and hydrophilicity and size of mobilisable particles. Domains of particle attachment, sliding, rolling, and lifting were found from a phase

Turbulence generated by aquatic vegetation in lakes, estuaries, and rivers can significantly alter the flow structure throughout the entire water column, affecting gas transfer mechanisms at the air-water interface, thus modifying indicators of water quality. A series of laboratory experiments with rigid, acrylic cylinder arrays to mimic vegetation was conducted in a recirculating Odell-Kovasznay type

Twelve actual evaporation datasets are evaluated for their ability to improve the performance of the fully distributed mesoscale Hydrologic Model (mHM). The datasets consist of satellite-based diagnostic models (MOD16A2, SSEBop, ALEXI, CMRSET, SEBS), satellite-based prognostic models (GLEAM v3.2a, GLEAM v3.3a, GLEAM v3.2b, GLEAM v3.3b), and reanalysis (ERA5, MERRA-2, JRA-55). Four distinct multivariate

Lateral unsaturated flow was found to contribute little to rainfall-induced shallow landslides, but it can be important when the landslides are triggered under conditions of varying boundary water level. By considering a hillslope with an impermeable bedrock and rising water level at the uphill/downhill boundary, this study explored how lateral unsaturated flow would affect slope stability. Three different

As river flow debouches into the sea, it is affected by the tidal fluctuation at the estuary mouth, which results in a backwater zone, where the rising residual water level (averaged over a lunar day) in the landward direction is changing periodically. This is known as backwater hydrodynamics, especially the variation of the residual water levels that are controlled by the interplay between tide and

The collapse of a submerged granular pile under the action of the gravity force is investigated by means of a depth averaged model (hereinafter DAM). The granular pile is a mixture of solid particles and ambient liquid, totally covered by a layer of ambient liquid. The DAM is obtained by depth-averaging the mass and momentum equations of the mixture and the ambient fluid and is aimed at providing estimates

Pseudo-plastic fluids exhibit a non-linear stress-strain relationship which can provoke large, localized viscosity gradients. For the flow of such fluids in porous media the consequence is a strong variability of the effective permeability with porosity, angle of the macroscopic pressure gradient, and rheological parameters of the fluid. Such a variability is investigated on the basis of adjoint homogenization

Quantifying and mapping regional and global irrigation water consumption have attracted considerable attention from researchers. The existing models cannot accurately estimate the long-term sequence of irrigation water consumption at the regional scale. Despite being among the countries with the largest irrigated areas, China has not been specifically evaluated in terms of high-resolution long-term

Introduction of fine sand into gravel-bed rivers has surged tremendously as a result of human activities. The ecological concerns for river habitats have called for suitable management strategies for cleaning these particles from gravel bed rivers. This paper is concerned with determining movement threshold of fine sand initially entrapped in the spaces between large immobile elements resembling a

Forecasting water level is an extremely important task as it allows to mitigate the effects of floods, reduce and prevent disasters. Physically based models often give good results but they require expensive computational time and various types of hydro-geomorphological data to develop the forecasting system. Alternatively, data driven forecasting models are usually faster and easier to build. During

Quantitatively evaluating the gravity effect on spontaneous imbibition in fractured porous media is challenging due to the lack of a suitable model and the complex fracture distribution. In this paper, new analytical solutions for spontaneous imbibition in a single fracture with and without gravity are first derived. Then, assuming a fractal distribution of fractures and with the tortuous plate fracture

When calving icebergs interact with water, waves of tens of meters in height, so-called iceberg-tsunamis (IBTs), may be generated. Recent examples include an IBT which reached an amplitude of 45 to 50 m in Eqip Sermia, Greenland, in 2014. A novel numerical methodology and unique large-scale laboratory experiments are presented to investigate the generation and propagation of such IBTs. In the laboratory

Regions with vulnerabilities in their food-energy-water nexus may also have economic vulnerabilities. Actions to reduce environmental vulnerabilities may conflict with community concerns about jobs and economic vitality. To account for these potentially competing concerns, this study builds a nexus simulation model that incorporates water, energy, food, and detailed economic data. The model can be

The Generalized Integral Transform Technique (GITT) is employed in combination with a single domain reformulation strategy and a coupled eigenvalue problem expansion base to construct analytical or hybrid numerical-analytical solutions for flow and transport in fractured porous media. The problem formulation and provided formal solutions encompass both continuum (multi-porosity/multi-permeability)