Numerical modeling of groundwater flow based on explicit and fully implicit schemes of finite volume method

This paper documents a novel numerical model for calculating the behavior of unsteady, one-dimensional groundwater flow by using the finite volume method. The developed model determined water table fluctuations for different scenarios as follows: Drainage and recession from an unconfined aquifer, and water table fluctuations above an inclined leaky layer due to ditch recharge with a constant and variable upper boundary condition. The Boussinesq equation, which is the governing equation in this domain, is linearized and solved numerically in both of the explicit and fully implicit conditions. Meanwhile, the upwind scheme is employed to discretize the governing equation. The computed outcomes of both the explicit and implicit approaches agreed well with the results of analytical solution and laboratory experiments. The main reason is that in the first half of simulation process explicit scheme obtains slightly better results and in the second half of the simulation process fully implicit scheme predicts more reliable outcomes that are hidden in the neighbor node points. As a final point, the numerical outcomes confirm that the developed model is capable of calculating satisfactory outcomes in engineering and science applications.

Prediction of Slow-Moving Landslide Mobility Due to Rainfall Using a Two-Wedges Model

In the present study, the landslides cyclically reactivated by water-table oscillations due to rainfall are dealt with. The principal kind of motion that usually characterizes such landslides is a slide with rather small velocity. As another feature, soil deformations are substantially accumulated inside a narrow shear zone situated below the landslide body so that the latter approximately slides rigidly. Within this framework, a new approach is developed in this paper to predict the mobility of this type of landslides due to rainfall. To this end, a two-wedges model is used to schematize the moving soil mass. Some analytical solutions are derived to link rain recordings with water-table fluctuations and in turn to landslide displacements. A well-documented landslide frequently activated by rainfall is studied to check the forecasting capacity of the proposed method.

Analyzing Trends in Water Table Elevations at the Marcell Experimental Forest, Minnesota, U.S.A.

Water table fluctuations in peatlands are closely coupled with the local climate setting and drive critical ecosystem processes such as nutrient cycling. In Minnesota, USA, peatlands cover ten percent of the surface area, approximately 2.5 million hectares, some of which are actively managed for forest products. To explore the relationship between peatland water tables and precipitation, long-term data (1961 to 2019) were used from the Marcell Experimental Forest in northern Minnesota. Starting in 1961, water table data from seven peatlands, including two types of peatlands (bogs and fens), were measured. We used the Theil-Sen estimator to test for monotonic trends in mean monthly water table elevations for individual peatlands and monthly precipitation. Water levels in bogs were both more variable and had mean water table elevations that were closer to the surface. Individual trends of water table elevations differed among peatlands. Water table elevations increased over time in three of the bogs studied and decreased over time in two of the bogs studied. Trends within fens were notably nonlinear across time. No significant linear trend was found for mean monthly precipitation between 1961 and 2019. These results highlight differences in peatlands types, local physiography, and the importance of understanding how changes in long-term dynamics coupled with changing current conditions will influence the effects of water table fluctuations on ecosystem services. The variability of water table elevations in bogs poses potential difficulties in modeling these ecosystems or creating adaptive management plans. KEYWORDS: Peatlands; Hydrology; Water tables; Bogs; Fens; Monitoring; Minnesota; Climate Change

Predicting preferential flow and water table fluctuations in karst systems using film-flow theory and source-responsive models

<p>The locally focused dissolution of the rock material (e.g., below dolines and dry valleys) in karst systems and in general percolating clusters of fractures in consolidated aquifer systems trigger the development of preferential flow paths in the vadose zone. Rainfall events may initiate rapid mass fluxes via macropores and fractures (e.g., as gravitationally-driven films) that lead to source-responsive water table fluctuations and comparably short residence times within the vadose zone. The degree of partitioning into a slow diffuse infiltration component and a rapid localized part depends, amongst others, on the hydraulic interaction of porous matrix and fracture domain as well as the geometrical characteristics of the fracture systems (e.g., persistence, connectivity) that are often difficult to obtain or unknown under most field conditions. Given their importance in water-resource management, specifically in arid and semi-arid regions (e.g., Mediterranean), it is desirable to recover such infiltration dynamics in porous-fractured systems with physically-based yet not overparameterized models. Here, we simulate water table fluctuations in a karst catchment in southwest Germany (Gallusquelle) using a source-responsive film flow model based on borehole and precipitation data. The model takes into account interfacial connectivity between slow and fast domain as well as phreatic zone discharge via classical recession analysis. This case study shows the potential importance of preferential flows while modeling water table responses in karst systems and recognizes the need for formulations other than those applied for a diffuse bulk fractured domain where infiltration patterns are assumed to be homogeneous without formation of infiltration instabilities along preferential pathways.</p>

General considerations for modeling water table dynamics in peatlands

<p>A basic and universal characteristic of peatlands is that the water table frequently rises near the surface of the soil profile. Surface peat is naturally loose and open-structured, and often has microtopographic features; the water table frequently rises above the level of local depressions. Therefore, water table fluctuations in peatlands cause rapid changes in the permeability and effective porosity of the medium through which flow occurs. We use a simple model based on Boussinesq's equation to explore the challenges that arise from these basic and universal physical aspects of peatland hydrology. We show that simulation of water table fluctuations in peatlands requires precipitation data with a high temporal resolution, and careful attention to the time derivative for accuracy of the mean water tables and correct water balance for two reasons. First, large vertical gradients in specific yield can result in large mass balance errors analogous to errors from naive discretization of the Richards equation; a change of variables from water table elevation to water storage can eliminate these errors and also speed up calculations by allowing larger time steps. Second, large vertical gradients in permeability near the peat surface cause a strongly nonlinear response to precipitation, so that time-averaged precipitation data or neglect of diurnal cycles of evapotranspiration results in erroneously high water levels, and careful time stepping is required around rain storms.  Consideration of these features of peatland hydrology results in efficient hydrologic models that can be used to predict spatial and temporal patterns in greenhouse gas uptake and emissions in peatlands.</p>