In the marine ecological system, the prime role of water management and durability of an ecosystem is being played by the vegetation patches. The vegetation patches in open channels can significantly affect the flow velocity, discharge capacity and hinder energy fluxes, which ultimately helps in controlling catastrophic floods. In this study, the numerical simulation for turbulent flow properties, i.e. velocity distribution, Reynolds stresses and Turbulent Intensities (TI) near the circular vegetation patches with progressively increasing density, were performed using the computational fluid dynamics (CFD) code ANSYS FLUENT. For examination of the turbulent flow features in the presence of circular patches with variable densities, Reynolds averaged Navier-Stokes equations, and Reynolds stress model (RSM) were employed. The numerical investigation was performed in the presence of in-line emergent and submerged patches having variable vegetation density in the downstream direction. Two of the cases were investigated with three circular patches having a clear gap to patch diameter ratio of La/D = 1 (where La is the clear spacing between the vegetation patches and D is the diameter of the circular patch), and the other two cases were analyzed with two patches having a clear gap ratio of La/D = 3. The case with a clear gap ratio (La/D = 3) showed 10.6% and 153% inflation in the magnitude of longitudinal velocity at the downstream of the sparse patch (aD = 0.8) and upstream of the dense patch (aD = 3.54), respectively (where aD is the flow blockage, in which “a” represents the patch frontal area and “D” represents the patch diameter). The velocity was reduced to 94% for emergent and 99% for submerged vegetation due to successive increase in vegetation density made by introducing a middle patch which reduced the clear gap ratio (La/D = 1). For La/D = 1, the longitudinal velocities at depth z = 15cm were increased by 319% than at depth z = 6cm at the downstream of the dense patch (aD = 3.54). Whereas it was observed to 365% higher in the case of La/D = 3. The magnitude of turbulent characteristics was observed 36% higher for submerged vegetation cases having a clear gap ratio of La/D = 1. The successive increase in the patch density reduced the Reynolds stresses, turbulent kinetic energy and turbulent intensities significantly within the gap region. The major reduction in the flow velocities and turbulent properties in the gaps provides a stable environment for aquatic ecosystems nourishment and fosters sediment deposition, and supports further vegetation growth.
Variable rate technology (VRT) in nutrient management has been developed in order to apply crop inputs according to the required amount of fertilizers. Meteorological conditions rarely differ within one field; however, differences in soil conditions responding to precipitation or evaporation results within field variations. These variations in soil properties such as moisture content, evapotranspiration ability, etc. requires site-specific treatments for the produced crops. There is an ongoing debate among experts on how to define management zones as well as how to define the required amount of fertilizers for phosphorus and nitrogen replenishment for winter wheat (Triticum aestivum L.) production. For management zone delineation, vegetation based or soil based data collection is applied, where various sensor technology or remote sensing is in help for the farmers. and The objective of the study reported in this paper was to investigate the effect of soil moisture data derived from Sentinel-2 satellite images moisture index and variable rate phosphorus and nitrogen fertilizer by means of variable rate application (VRA) in winter wheat in Mezőföld, Hungary. Satellite based moisture index variance at the time of sowing has been derived, calculated and later used for data comparison. Data for selected points showed strong correlation (R2 = 0.8056; n = 6) between moisture index and yield, however generally for the whole field correlation does not appear. Vegetation monitoring has been carried out by means of NDVI data calculation. On the field level, as indicated earlier neither moisture index values at sowing nor vegetation index data was sufficient to determine yield. Winter wheat production based on VRA treatment resulted significant increase in harvested crop: 5.07 t/h in 2013 compared to 8.9 t/ha in 2018. Uniformly managed (control) areas provided similar yield as VRA treated areas (8.82 and 8.9 t/ha, respectively); however, the input fertilizer was reduced by 108 kg/ha N and increased by 37 kg/ha P.
The paper describes a mathematical and physical modelling of flow of complex slurries in pipelines, i.e. a flow of slurries composed of solids covering a very broad range of particle sizes that overlaps more than one flow patterns – non-Newtonian, pseudohomogeneous, heterogeneous and fully-stratified. A typical examples are residual products (“tailings”) from mining industry with normal average particle sizes of 20 to 100 μm or more. Experimental results of flows of complex slurries composing of non-Newtonian carrier fluid and three fractions of glass particles in 50 mm pipe are presented. Depending on the particle size, particles show different flow patterns and therefore considerable differences in pressure drops. Fine particles tend behave as a homogeneous matter, while coarser particles exhibit heterogeneous behaviour and even coarser particles form a sliding bed. A mathematical 3-component predictive model for turbulent flow of complex slurries is presented based on well-established semi-empirical formulae developed originally for flows with Newtonian carrier. The predicted values of pressure drops show very reasonable agreement with experimental results and indicate suitability of the model for engineering practice.
The two-dimensional particle image velocimetry (PIV) data are inevitably contaminated by noise due to various imperfections in instrumentation or algorithm, based on which the well-established vortex identification methods often yield noise or incomplete vortex structure with a jagged boundary. To make up this deficiency, a novel method was proposed in this paper and the efficiency of the new method was demonstrated by its applications in extracting the twodimensional spanwise vortex structures from 2D PIV data in open-channel flows. The new method takes up a single vortex structure by combining model matching and vorticity filtering, and successfully locates the vortex core and draws a streamlined vortex boundary. The new method shows promise as being more effective than commonly used schemes in open-channel flow applications.
Backward erosion piping is driven by seepage forces acting on the soil grains at the downstream end of the seepage path. A new device for the laboratory testing of backward erosion progression was developed and tested. The device consists of a plexiglass prism at which the seepage path has been predefined. The prism was equipped with an inflow consisting of gravel separated from tested sand by a strainer. The hydraulic gradient along the seepage pipe was observed by a set of piezometers and pressure cells, and the seepage discharge was measured volumetrically. The transported sediment was trapped in a vertical cone located downstream from the device. The progression of the seepage path, the piezometric heads and the trapped material was observed by two synchronous cameras. 15 trial tests have been carried out to date, and from these, the interim results are presented.
Fishway design not only takes into account the swimming abilities of target fishes, but also considers the hydrodynamic characteristics within the fishway. In this study, the flow fields of one vertical-slot fishway (i.e. VSF), five T-shape fishways (i.e. TSF-1~TSF-5) and two H-shape fishways (i.e. HSF-1 and HSF-2) are numerically simulated by solving the three-dimensional Reynolds-averaged Navier-Stokes equations and the K-Omega-SST turbulence model. The numerical results clearly indicate that the hydrodynamic properties of HSF-2 are overall superior to the remaining seven cases, in terms of the time-averaged flow pattern, the time-averaged velocity magnitude, the depth-mean timeaveraged velocity magnitude along the vertical-slot section, the volume percentages of the time-averaged velocity magnitude less than some critical values, and the distribution of the time-averaged turbulent kinetic energy. Therefore, HSF-2 is more friendly for fishes with relatively smaller sizes and weaker swimming capacities to transfer upstream. The novel HSF-2 is firstly proposed in this paper, and it is naturally designed during the process of improving the flow regime. Furthermore, the generalizability of the superiority of HSF-2 over VSF and the original T-shape fishway (i.e. TSF-1) has been exhibited with the aid of the numerical results of four operating conditions (i.e. Q = 400 L/s, 600 L/s, 800 L/s and 1000 L/s).
Accurate estimation of precipitation in mountain catchments is challenging due to its high spatial variability and lack of measured ground data. Weather radar can help to provide precipitation estimates in such conditions. This study investigates the differences between measured and radar-estimated daily precipitation in the mountain catchment of the Jalovecký Creek (area 22 km2, 6 rain gauges at altitudes 815–1900 m a.s.l.) in years 2017–2020. Despite good correlations between measured and radar-based precipitation at individual sites (correlation coefficients 0.68–0.90), the radar-estimated precipitation was mostly substantially smaller than measured precipitation. The underestimation was smaller at lower altitude (on average by –4% to –17% at 815 m a.s.l.) than at higher altitudes (–35% to –59% at 1400–1900 m a.s.l.). Unlike measured data, the radar-estimated precipitation did not show the differences in precipitation amounts at lower and higher altitudes (altitudinal differences). The differences between the measured and radar-estimated precipitation were not related to synoptic weather situations. The obtained results can be useful in preparation of more accurate precipitation estimates for the small mountain catchments.
The design and evaluation of algorithms for adaptive stochastic control of the reservoir function of a water reservoir using an artificial intelligence method (learned fuzzy model) are described in this article. This procedure was tested on the Vranov reservoir (Czech Republic). Stochastic model results were compared with the results of deterministic management obtained using the method of classical optimisation (differential evolution). The models used for controlling of reservoir outflow used single quantile from flow duration curve values or combinations of quantile values from flow duration curve for determination of controlled outflow. Both methods were also tested on forecast data from real series (100% forecast). Finally, the results of the dispatcher graph, adaptive deterministic control and adaptive stochastic control were compared. Achieved results of adaptive stochastic management were better than results provided by dispatcher graph and provide inspiration for continuing research in the field.
We analyse water balance, hydrological response, runoff and snow cover characteristics in the Jalovecký Creek catchment (area 22 km2, mean elevation 1500 m a.s.l.), Slovakia, in hydrological years 1989–2018 to search for changes in hydrological cycle of a mountain catchment representing hydrology of the highest part of the Western Carpathians. Daily air temperature data from two meteorological stations located in the studied mountain range (the Tatra Mountains) at higher elevations show that the study period is 0.1°C to 2.4°C warmer than the climatic standard period 1951–1980. Precipitation and snow depth data from the same stations do not allow to conclude if the study period is wetter/drier or has a decreasing snow cover. Clear trends or abrupt changes in the analysed multivariate hydrometric data time series are not obvious and the oscillations found in catchment runoff are not coherent to those found in catchment precipitation and air temperature. Several time series (flashiness index, number of flow reversals, annual and seasonal discharge maxima, runoff coefficients) indicate that hydrological cycle is more dynamic in the last years of the study period and more precipitation runs off since 2014. The snow cover characteristics and climatic conditions during the snow accumulation and melting period do not indicate pronounced changes (except the number of days with snowfall at the Kasprowy Wierch station since 2011). However, some data series (e.g. flow characteristics in March and June, annual versus summer runoff coefficients since 2014) suggest the changes in the cold period of the year.
δ18O in precipitation at station Liptovský Mikuláš (about 8.5 km south from the outlet of the Jalovecký Creek catchment) remains constantly higher since 2014 that might be related to greater evaporation in the region of origin of the air masses bringing precipitation to the studied part of central Europe. Increased δ18O values are reflected also in the Jalovecký Creek catchment runoff. Seasonality of δ18O in the Jalovecký Creek became less pronounced since 2014. The most significant trends found in annual hydrological data series from the catchment in the study period 1989–2018 have the correlation coefficients 0.4 to 0.7. These trends are found in the number of flow reversals (change from increasing to decreasing discharge and vice versa), June low flow, number of simple runoff events in summer months (June to September) and the flashiness index. The attribution analysis suggests that drivers responsible for the changes in these data series include the number of periods with precipitation six and more days long, total precipitation amount in February to June, number of days with precipitation in June to September and total precipitation in May on days with daily totals 10 mm and more, respectively. The coefficients of determination show that linear regressions between the drivers and supposedly changed data series explain only about 31% to 36% of the variability. Most of the change points detected in the time series by the Wild Binary Segmentation method occur in the second and third decades of the study period. Both hydrometric and isotopic data indicate that hydrological cycle in the catchment after 2014 became different than before.
Study about the mechanical energy balance and the energy loss of 3-D turbulent flows in open-channels has its own complexities. The governing equation of the mechanical energy in turbulent flows has been previously known and includes turbulence parameters that their calculations or measurements are not easy. In this study, a form of the total mechanical energy equation that leads to a number of significant physical insights is analytically investigated, from which analytical relationships for the energy loss estimation in 3-D turbulent flows are defined. The effect of different turbulence parameters is reflected on the new relationships and analyzed by equalizations replacing unknown correlations with closure approximations using the numerical turbulence simulation. In order to investigate the application of the analytical relationships, numerical simulations are performed by using OpenFOAM software to solve the Navier-Stokes equations with the RSM turbulence model in open-channels with different geometries. Then, the contribution of the turbulence parameters to the total mechanical energy balance is evaluated in uniform and nonuniform turbulent flows and their difference is analyzed, that leads to identify the parameters affecting the friction and local losses. The results demonstrate that the magnitudes of the turbulent diffusion, the work done by the viscous stresses pertaining to the mean motion and the viscous diffusion of the turbulence energy are substantially smaller than the other terms of the total energy equation for turbulent flows in open-channels with different geometries, while the effect of the variations of the turbulence kinetic energy and the work done by the turbulence stresses, that has not been considered in the previous mechanical energy equations, is more important in complex flows. From a practical viewpoint, in order to study the details of the total mechanical energy balance and the energy loss in 3-D turbulent flows with the presence of the secondary currents, the proposed method can be useful.
Authors propose a beneficial methodology for hydrological planning in their study. Prospective evaluations of the basins' net capacity can be done using the technique presented. The HEC-HMS (Hydrologic Modelling System) software can be used to estimate in a basin, the sediment emitted. For a certain precipitation, this methodology allows estimating, within a certain range, the gradual blockage of a reservoir, and even a projected date for total blockage. This has some applications to adopt corrective measures that prevent or delay the planned blockage deadlines. The model is of the semi-distributed type, estimating the generation and emission of sediments by sub-basins. The integration of different return periods in HEC-HMS with a semi-distributed model by sub-basins and the application of a mathematical model are the differentiating element of this research. The novelty of this work is to allow prognosing the reservoir sedimentation rate of basins in a local and regional scale with a medium and large temporary framework. The developed methodology allows public institutions to take decisions concerning hydrological planning. It has been applied to the case of "Charco Redondo" reservoir, in Cádiz, Andalusia, in southern Spain. Applying the methodology to this case, an average soil degradation of the reservoir basin has been estimated. Therefore, it is verified that in 50 years the reservoir is expected to lose 8.4% of its capacity.
Vapor condensation, whether due to dew or fog, may add a stable and important source of water to deserts. This was also extensively assessed in the Negev, regarded as a dew desert. Dew deserts necessitate a large reservoir of vapor, and are therefore confined to near oceans or seas. Yet, examples of such deserts are scarce. Here we try to assess whether the Tabernas Desert in SE Spain can be regarded as a dew desert, and may therefore facilitate the growth of certain organisms that otherwise would not survive the dry season. We analyze some of the abiotic conditions of four relatively dry months (June, July, August, September) in the Tabernas and Negev deserts (with the Negev taken as an example of a dew desert) during 2003–2012. The analysis showed substantially lower values of relative humidity (by 10–13%) in the Tabernas in comparison to the Negev, with RH ≥95% being on average only 0.9–1.1 days a month in the Tabernas in comparison to 9.7–13.9 days in the Negev. Our findings imply that the Tabernas Desert cannot be regarded as a dew desert, suggesting that rain will be the main factor responsible for the food web chain in the Tabernas.
The presence of biocrusts changes water infiltration in the Mu Us Desert. Knowledge of the hydraulic properties of biocrusts and parameterization of soil hydraulic properties are important to improve simulation of infiltration and soil water dynamics in vegetation-soil-water models. In this study, four treatments, including bare land with sporadic cyanobacterial biocrusts (BL), lichen-dominated biocrusts (LB), early-successional moss biocrusts (EMB), and latesuccessional moss biocrusts (LMB), were established to evaluate the effects of biocrust development on soil water infiltration in the Mu Us Desert, northwest of China. Moreover, a combined Wooding inverse approach was used for the estimation of soil hydraulic parameters. The results showed that infiltration rate followed the pattern BL > LB > EMB > LMB. Moreover, the LB, EMB, and LMB treatments had significantly lower infiltration rates than the BL treatment. The saturated soil moisture (θs) and shape parameter (αVG) for the EMB and LMB treatments were higher than that for the BL and LB treatments, although the difference among four treatments was insignificant. Water retention increased with biocrust development at high-pressure heads, whereas the opposite was observed at low-pressure heads. The development of biocrusts influences van Genuchten parameters, subsequently affects the water retention curve, and thereby alters available water in the biocrust layer. The findings regarding the parameterization of soil hydraulic properties have important implications for the simulation of eco-hydrological processes in dryland ecosystems.
Agriculture faces several challenges to use the available resources in a more environmentally sustainable manner. One of the most significant is to develop sustainable water management. The modern Internet of Things (IoT) techniques with real-time data collection and visualisation can play an important role in monitoring the readily available moisture in the soil. An automated Arduino-based low-cost capacitive soil moisture sensor has been calibrated and developed for data acquisition. A sensor- and soil-specific calibration was performed for the soil moisture sensors (SKU:SEN0193 - DFROBOT, Shanghai, China). A Repeatability and Reproducibility study was conducted by range of mean methods on clay loam, sandy loam and silt loam soil textures. The calibration process was based on the data provided by the capacitive sensors and the continuously and parallelly measured soil moisture content by the thermo-gravimetric method. It can be stated that the response of the sensors to changes in soil moisture differs from each other, which was also greatly influenced by different soil textures. Therefore, the calibration according to soil texture was required to ensure adequate measurement accuracy. After the calibration, it was found that a polynomial calibration function (R2 ≥ 0.89) was the most appropriate way for modelling the behaviour of the sensors at different soil textures.
The objectives of the study were to: (1) assess the strength of associations of direct CO2 and N2O emissions with the seasonal variations in the relevant soil properties under both tillage systems; 2) evaluate how CT and RT affect magnitudes of seasonal CO2 and N2O fluxes from soil. Field studies were carried out on plots for conventional tillage (up to 0.22–0.25 m) and reduced tillage (up to 0.10–0.12 m) during the growing season and post-harvest period of red clover. The results showed that daily CO2 emissions significantly correlated only with soil temperature during the growing season under conventional and reduced tillage. Soil temperature demonstrated its highest influence on daily N2O emissions only at the beginning of the growing season in both tillage systems. There were no significant inter-system differences in daily CO2 and N2O emissions from soil during the entire period of observations. Over the duration of post-harvest period, water-filled pore space was a better predictor of daily CO2 emissions from soils under CT and RT. The conventional and reduced tillage did not cause significant differences in cumulative N2O and CO2 fluxes from soil.
In an open channel with a mobile bed, intense transport of bed load is associated with high-concentrated sediment-laden flow over a plane surface of the eroded bed due to high bed shear. Typically, the flow exhibits a layered internal structure in which virtually all sediment grains are transported through a collisional layer above the bed. Our investigation focuses on steady uniform turbulent open-channel flow with a developed collisional transport layer and combines modelling and experiment to relate integral quantities, as the discharge of solids, discharge of mixture, and flow depth with the longitudinal slope of the bed and the internal structure of the flow above the bed. A transport model is presented which considers flow with the internal structure described by linear vertical distributions of granular velocity and concentration across the collisional layer. The model employs constitutive relations based on the classical kinetic theory of granular flows selected by our previous experimental testing as appropriate for the flow and transport conditions under consideration. For given slope and depth of the flow, the model predicts the total discharge and the discharge of sediment. The model also predicts the layered structure of the flow, giving the thickness of the dense layer, collisional layer, and water layer. Model predictions are compared with results of intense bed-load experiment carried out for lightweight sediment in our laboratory tilting flume.
Sand-water slurry was investigated on an experimental pipe loop of inner diameter D = 100 mm with the horizontal, inclined, and vertical smooth pipe sections. A narrow particle size distribution silica sand of mean diameter 0.87 mm was used. The experimental investigation focused on the effects of pipe inclination, overall slurry concentration, and mean velocity on concentration distribution and deposition limit velocity. The measured concentration profiles showed different degrees of stratification for the positive and negative pipe inclinations. The degree of stratification depended on the pipe inclination and on overall slurry concentration and velocity. The ascending flow was less stratified than the corresponding descending flow, the difference increasing from horizontal flow up to an inclination angle of about +30°. The deposition limit velocity was sensitive to the pipe inclination, reaching higher values in the ascending than in the horizontal pipe. The maximum deposition limit value was reached for an inclination angle of about +25°, and the limit remained practically constant in value, about 1.25 times higher than that in the horizontal pipe. Conversely, in the descending pipe, the deposition limit decreased significantly with the negative slopes and tended to be zero for an inclination angle of about −30°, where no stationary bed was observed.
The event runoff coefficient (Rc) and the recession coefficient (tc) are of theoretical importance for understanding catchment response and of practical importance in hydrological design. We analyse 57 event periods in the period 2013 to 2015 in the 66 ha Austrian Hydrological Open Air Laboratory (HOAL), where the seven subcatchments are stratified by runoff generation types into wetlands, tile drainage and natural drainage. Three machine learning algorithms (Random forest (RF), Gradient Boost Decision Tree (GBDT) and Support vector machine (SVM)) are used to estimate Rc and tc from 22 event based explanatory variables representing precipitation, soil moisture, groundwater level and season. The model performance of the SVM algorithm in estimating Rc and tc is generally higher than that of the other two methods, measured by the coefficient of determination R2, and the performance for Rc is higher than that for tc. The relative importance of the explanatory variables for the predictions, assessed by a heatmap, suggests that Rc of the tile drainage systems is more strongly controlled by the weather conditions than by the catchment state, while the opposite is true for natural drainage systems. Overall, model performance strongly depends on the runoff generation type.