Groundwater Flow Research Articles

ОЦЕНКА ВЗАИМОСВЯЗИ ПОВЕРХНОСТНЫХ И ПОДЗЕМНЫХ ВОД РЕКИ СЫРДАРИЯ (В ПРЕДЕЛАХ КАЗАХСТАНА)

This study analyzes and evaluates the relationship of surface and groundwater of the Syrdarya River within Kazakhstan. The share of groundwater recharge of the river was calculated by two methods: by genetically partitioning the total runoff hydrograph and by calculating low f low. The difference of groundwater flow into the river between the selected hydrostations is also calculated. A comparison of the obtained results for the quasi-natural and disturbed periods is made. In natural conditions, the infiltration of river runoff into underground horizons prevailed along the length of the Syrdarya River, however, under the influence of anthropogenic activity, the picture changed to the opposite: the discharge of groundwater prevails along the length of the Syrdarya River, including the emergence of filtration waters below the dam of the Shardara HPP, collector-drainage and discharge waters from irrigation f ields of the adjacent territory. Also, apparently, the truncation of spring flood peaks leads to an increase in the groundwater level and an increase in the share of underground river supply. Recommendations for the application of the methods used in the work are given.

Adaptive management of borehole heat exchanger fields under transient groundwater flow conditions

Uncontrolled heat extraction by multiple interacting borehole heat exchangers (BHEs) in high-density energy-use districts can lead to undesirable thermal conditions in the subsurface which can affect both system performance and regulatory compliance. The difficulty in controlling heat extraction arises in particular from predictive uncertainties, such as when forecasting trends in energy demand or groundwater flow. In this study, a combined simulation-calibration-optimization framework is introduced to consider BHE fields with the presence of a transient groundwater flow regime. In the first part, a semi-analytical modeling technique is proposed based on temporal superpositioning of variable flow conditions. Two synthetic case studies verify its accuracy under different groundwater fluctuation patterns. The mean absolute error of the proposed model in comparison to numerical calculation does not exceed 0.18 K over ten years of operation. In the second part, the model is augmented by a parameter estimation algorithm that is employed for continuous model updating. The benefit of resolving transient flow conditions is demonstrated by using this approach for monthly optimization of individual BHE heat extraction. The result of dynamic optimization compared to a synthetic case without calibration shows a 10 % lower imposed temperature change in the subsurface.

Seismically-induced groundwater flow into confining clays: An experimental and numerical study

This study investigates seismically-induced hydraulic responses in an aquifer system focusing on the hydraulic interaction between a model aquifer and confining clays following shaking as experienced during seismic events. These responses are being investigated by both experimental and numerical methods. In contrast to the conventional assumption assuming a clay layer as a hydraulic barrier to seismically-induced hydrological responses in an aquifer, this study highlights that the excess pore pressure in the aquifer generated by shaking may facilitate water flow crossing the overlying clay layer. This hypothesis is supported by experimental results that show measurable surface settling, excess pore pressure and step-like wellbore water level increase following shaking-induced sand consolidation, and notably, upward flow of water from an aquifer into the clay layer indicated by sustained water level drops in the well after shaking experiments. Numerical simulations further reveal that perturbations influence water flow from the aquifer to the clay layer, modifying pressure distribution within both layers. Key factors such as aquifer flow regime, hydraulic conductivity of the aquifer and clay layer, and their effects on flow dynamics are addressed by numerical modeling. This study provides insights into hydraulic responses and contributes to risk assessments in the context of seismic perturbations.

Ditch cleaning in boreal catchments: Impacts on water chemistry and dissolved greenhouse gases in runoff

Ditch cleaning (DC) is a forestry practice commonly conducted in boreal regions that aim to lower groundwater tables (GWT) in waterlogged soils, thereby maintaining or improving forest growth. However, there is limited information on the impact of DC on water chemistry and dissolved greenhouse gases (GHG) in draining ditch networks. Based on a repeated synoptic sampling of a paired catchment design we here evaluated water chemistry and GHG data in ditch waters from 25 cleaned sites (being cleaned 1–4 years prior to sampling) and 25 non-cleaned reference sites (REF). The sampled sites were further selected to test whether there were any differences in the DC effect if the operations were conducted in forested or clear-cut areas. Across all sites, we found that DC sites exhibited higher pH, sulfate and calcium concentrations than REF sites. Also, lower dissolved carbon dioxide and higher nitrous oxide concentrations were observed in DC sites. In forested areas, DC sites exhibited significantly higher calcium and potassium concentrations, along with reduced levels of methylmercury and carbon dioxide. In clear-cuts, sulfate concentrations were significantly elevated in the DC sites. We suggest that the observed differences in water chemistry and GHG´s between DC and REF sites were induced by the DC and largely driven by lower GWT following DC, resulting in deeper groundwater flow paths through more mineral-rich soil layers, and altered redox conditions. Also, the removal of organic rich sediments and vegetation from the ditches themselves may affect water chemistry and GHG´s e.g. by decreasing the formation of methylmercury and carbon dioxide.

Hydrogeochemistry of aquifers in the northern portion of the Lagoa Santa Karst Environmental Protection Area and surroundings, state of Minas Gerais, Brazil

The studied region, located in the central part of Minas Gerais state, Brazil, encompasses the northern portion of the Lagoa Santa Karst Environmental Protection Area (APA) and adjacent areas to the west. It is notable for harboring important karstic aquifers and archaeological and paleontological sites of global significance. The research aims to define the key parameters that characterize the waters of the aquifers, quantifying their most frequent maximum and minimum concentrations. For this purpose, water sampling and hydrochemical analyses were conducted at 51 water points. Hierarchical clustering analysis facilitated the association of samples with aquifer units. Utilizing Piper and Stiff diagrams, boxplot graphs, and proportional symbol maps enabled the classification of water and the determination of its hydrochemical signatures. All sampled waters are bicarbonated, with 87% being calcic. Waters associated with metacalcareous rocks exhibit higher concentrations of calcium (Ca2⁺), bicarbonate (HCO₃⁻), electrical conductivity (EC), total dissolved solids (TDS), and hydrogenionic potential (pH) compared to waters associated with non-carbonate rocks. Waters from the Lagoa Santa Member differ from those of the Pedro Leopoldo Member by higher concentrations of Ca2⁺ and HCO₃⁻ and slightly lower concentrations of sodium (Na⁺), magnesium (Mg2⁺), and silica (SiO₄2⁻). Non-carbonate waters related to granite-gneisses differ from those associated with metapelites by slightly higher concentrations of potassium (K⁺), sodium (Na⁺), silica (SiO₄2⁻), and barium (Ba2⁺). Additionally, hydrochemical variations have been linked to the regional groundwater flow (from west to east). The results obtained have enabled advancements in the hydrogeochemical understanding of the study area and can be valuable for comparative studies in other karst regions, contributing to more effective water resource management aimed at environmental conservation and water sustainability, a growing global concern.

Hydrogeological assessment and aquifer potential of a coastal area, South Nigeria: insights from VES surveys and spatial analysis

ABSTRACT This study provides a comprehensive analysis of aquifer properties using vertical electrical sounding (VES) surveys conducted in Eket, Nigeria, to address the challenges of groundwater management in coastal areas. The VES method was employed for hydrogeo-electrical analysis and to investigate the hydrogeological dynamics of the region. The results reveal varying subsurface layers with distinct resistivity values. Spatial distribution analysis indicates densely settled areas exhibit lower resistivities, while sparsely settled regions show higher resistivities. The spatial location labeled EK2 features the largest aquifer, measuring 246.70 m, predominantly located in the central research area. Aquifer depths range from 87.10 to 250.20 m, with variations in thickness highlighting geographical heterogeneity and its impact on groundwater transport. These spatial differences in groundwater characteristics are crucial for understanding groundwater movement, storage, and extraction potential. The study underscores significant variability in porosity, reflecting differences in the aquifers’ water storage capacity and susceptibility to contamination. These findings have important implications for groundwater flow rates and extraction feasibility. The research provides essential data for hydrogeological investigations and groundwater resource management, emphasizing the intricate relationship among geological formations, aquifer properties, human activities, and the potential risks of contamination in variable porosity zones.

Using Expert Participation to Evaluate the Accuracy of Hand-Drawn Water-Table Maps.

Water-table maps are fundamental to hydrogeological studies and a manual, hand-drawn method is still commonly used to produce them. Despite this, the accuracy and variability of such maps have received little attention in international literature. In a unique experiment, 63 groundwater professionals drew water-table equipotential contours based on the same dataset of point measurements and were asked to infer flow directions and predict groundwater elevations at predefined locations. The root mean squared error (RMSE) for the average map calibration data was 10.5 m, which is accuracy comparable to numerical groundwater models. This study confirmed that to produce hand-drawn water-table maps, practitioners seek to not only fit the spatial data, but also to conform to their own cognitive model of hydrogeological concepts and processes. The calibration accuracy increased with experience; from a RMSE of 13.3 m for practitioners with 0-3 years of experience to a RMSE of 9.2 m for those with four or more years. Despite considerable variability in the style of the hand-drawn water-table maps, the maps were consistent in their representation of the dominant regional groundwater flow directions. There was less consensus, however, in predicting the direction of surface water-groundwater interaction for a stream reach. Hand-drawn water-table mapping remains useful and valid, especially as a starting point for hydrogeological conceptualization, yet further work is required to resolve issues around transparency, repeatability, and reproducibility.

Hydrogeophysical Evaluation of the Karst Aquifer near the Western Edge of the Ring of Cenotes, Yucatán Peninsula

In this work, electrical resistivity tomography was carried out together with physical hydrogeology techniques to evaluate the karst aquifer in the northwest region of the Yucatán Peninsula in a study area near the western edge of the Ring of Cenotes of the Chicxulub Crater. In addition, based on a systematic compilation of open-access data of water levels reported for the peninsular aquifer, maps of groundwater isolines and groundwater flows were generated using IDW interpolation, Empirical Bayesian Kriging, and the Flow Net method. From these results, a shallow aquifer is observed, with the presence of heterogeneities such as possible dissolution conduits and/or flooded caverns, approximately 20 m below ground level, formed by the dissolution processes of limestone rocks. On a regional scale, the geomorphological influence of the Ring of Cenotes on groundwater flows was observed. In general, the flow directions observed from these maps coincide with those conceptualized for this region of the peninsular aquifer. Nevertheless, some differences were observed depending on the interpolation method used. Our results contribute to hydrogeological studies carried out in the periphery of this ring, where the vulnerability of the aquifer to anthropogenic contamination has been highlighted due to the intrinsic features of the karst environment.

Assessing the Impact of Groundwater Recharge on Underground Reservoir Replenishment in Eastern Uttar Pradesh, India

Groundwater is considered a fresh resource of water and its uses have tremendously increased in the recent past due to an increase in population, rapid urbanization, and industrialization. In India, the groundwater level is declining in some parts of the country due to over-exploitation, low or negligible recharge of aquifer systems, and unsustainable development of groundwater resources. The groundwater modeling is an important tool for studying the past and present groundwater behavior and in the development of future strategies for sustainable groundwater management plans. To study the Impact of groundwater recharge on the replenishment of underground reservoir, Ballia district of Uttar Pradesh has been selected which is one of the districts of the most populous state of India, Uttar Pradesh. An attempt has been made to develop a groundwater model using Modflow software to simulate the groundwater trends and predict future groundwater heads. The calibration and validation of the model were done for 5 years and 3 years respectively. The correlation coefficient for calibration and validation was found 0.85 and root mean square errors vary from 2.89 to 3.2m variation in future trends of groundwater heads. The results of the study show that the developed model can be effectively used to predict the future groundwater heads. The groundwater flow was observed from the northwest to southeast direction. It was predicted from the study that groundwater draft will increase by 10% with a decrease in groundwater level by approximately 0.24 m in the north-west direction by the year 2025. However, no impact was observed in the south side of the district and it was predicted that the groundwater level would remain the same in this zone during the next 3 years.

Uncertainty modeling and propagation for groundwater flow: a comparative study of surrogates

We compare sparse grid stochastic collocation and Gaussian process emulation as surrogates for the parameter-to-observation map of a groundwater flow problem related to the Waste Isolation Pilot Plant in Carlsbad, NM. The goal is the computation of the probability distribution of a contaminant particle travel time resulting from uncertain knowledge about the transmissivity field. The latter is modelled as a lognormal random field which is fitted by restricted maximum likelihood estimation and universal kriging to observational data as well as geological information including site-specific trend regression functions obtained from technical documentation. The resulting random transmissivity field leads to a random groundwater flow and particle transport problem which is solved realization-wise using a mixed finite element discretization. Computational surrogates, once constructed, allow sampling the quantities of interest in the uncertainty analysis at substantially reduced computational cost. Special emphasis is placed on explaining the differences between the two surrogates in terms of computational realization and interpretation of the results. Numerical experiments are given for illustration.

Investigating the impact of irrigation practices on hydrologic fluxes in a highly managed river basin

Irrigation practices and sources can have significant impacts on water resources and the hydrologic fluxes that control these resources. To better manage water resources and future water supply, the influence of irrigation practices and management on these hydrologic fluxes should be quantified in time and space at varying scales, under potential irrigation management practices. To fulfill this objective, a surface-subsurface modeling approach was applied to simulate watershed-scale hydrologic processes in the Cache la Poudre River Basin, Colorado, USA (4824 km2), in which both surface water irrigation and groundwater irrigation are prevalent. The model chosen for this study is the watershed model SWAT+, using the spatially distributed, physically based groundwater module gwflow, in which unconfined groundwater storage, flows, and interaction with land surface features are simulated using a collection of grid cells that represent control volumes of the aquifer. Major groundwater inflows and outflows include pumping, recharge, groundwater-channel exchange, groundwater-lake exchange, and tile drainage outflow. To investigate the impact of irrigation practices, detailed surface and groundwater irrigation routines and canal-aquifer interactions were added to the SWAT+ source code, requiring information of irrigation sources and irrigation canal locations throughout the river basin. Model calibration and testing was performed using monthly stream discharge and groundwater head. The calibrated model is used to quantify the impact of surface water and groundwater irrigation scenarios on water availability and hydrologic fluxes within the river basin. A total of 22 scenarios were conducted and grouped into five main groups: irrigation source, irrigation amount, irrigation type, canal bed thickness, and partial or full sealing of earthen irrigation canals. Using groundwater as the only irrigation source decreases groundwater discharge to streams (by 14 %) due to lowering groundwater levels; converting flood irrigation to sprinkler irrigation throughout the basin decreases surface runoff by 22 %; and sealing earthen canals leads to a lowering of groundwater levels, which decreases groundwater discharge to streams by 9 %, leading to an overall decrease in streamflow in the Cache la Poudre River and changes to temporal patterns in streamflow. Overall, irrigation amount and type and canal sealing have a small impact on total groundwater storage, compared to changes in the percent of fields irrigated by groundwater pumping. Results are helpful for informed decision-making in agriculture water management and can lead to sustainable, efficient, and equitable use of water resources, helping to address the challenges posed by water scarcity and environmental conservation.

Drivers of Daily Water Level Fluctuation of Shallow Groundwater in the Inner Delta of the River Danube

Groundwater flow systems are influenced by the changes in surface waters as well as climatic factors. These teleconnections significantly increase in cases of extreme weather conditions. To prepare and mitigate the effect of such phenomena, the background factors that create and influence natural processes must be recognized. In the present study, 94 shallow groundwater (SGW) wells’ water level time series were analyzed in the inner delta of the River Danube (Europe) the Szigetköz region to explore which factors contribute to the development of diurnal periodicity of SGW and what its drivers are. The relationship between surface meteorological processes and SGW dynamics in the Szigetköz region was investigated using hourly data from monitoring wells. Hourly water temperature data exhibited weak correlations with meteorological parameters. However, daily averaged data revealed stronger correlations, particularly between SGW levels and air temperature and potential evapotranspiration. Diurnal periodicity in SGW fluctuations correlated strongly with potential evapotranspiration. The study also demonstrated the role of capillary fringe dynamics in linking surface evapotranspiration with SGW fluctuations. Changes in groundwater levels, even small, can significantly affect soil moisture, vegetation, and ecosystem functioning, highlighting the sensitivity of the unsaturated zone to SGW fluctuations driven by surface processes.

Predicting heavy metal transport in groundwater around Lemna dumpsite: implications for residence utilizing borehole water in Cross River State, Nigeria

Groundwater is considered the most important natural resource to mankind. Groundwater constitutes an important part of the hydrological cycle and is more prone to pollution. Dumpsite located in close proximity to groundwater resources is highly susceptible to leachates pollution. Predicting the susceptibility of groundwater pollution is crucial to address industry-standard codes for groundwater flow, contaminant transport, local to regional-scale water quality, and source water protection issues. Therefore, predicting heavy metal transport in groundwater around Lemna dumpsite in Cross River State, Nigeria, was examined. Soil samples were purposively collected with a soil Auger, along a straight line at (5 m, 25 m and 50 m) in the dumpsite. Water samples were purposively collected from five (5) boreholes close to Lemna dumpsite. The study utilized pumping test method to obtain data for the analysis of heavy metal transport in groundwater. Data analysis of the laboratory results of soil and borehole water quality focuses on arsenic, lead, cadmium, chromium, nickel, and mercury. Paired sample t test was used to analyse the soil and borehole water quality. Visual Modflow was also used to analyse the solute transport of heavy metals in groundwater around Lemna dumpsite. The paired sample t test of the analysis of heavy metals in soil exhibited a significant difference (p < 0.05) compared to National Environmental Standard Regulation and Enforcement Agency limits. The paired sample t test of the analysis of heavy metals in borehole water exhibited a significant difference (p < 0.05) compared to World Health Organization limits. The significant level indicates contamination of the soil and borehole water. The findings revealed a spatial spread of 259.2000 m2/day, with the contaminant travelling up to 94,608 m2/year. The extent of heavy metals concentration exhibited a maximum of 0.991 mg/l to a minimum of (− 6.72 × 10–18 mg/l), with concentrations decreasing as the plume extend. The study recommends the need for remediation and stringent monitoring to mitigate heavy metal contamination of boreholes near Lemna dumpsite.

Modeling Contaminant Transport from a Red Mud Pond – A Case Study from South India

The present study focused on impact assessment of redmud pond on groundwater using numerical groundwater flow and contaminant transport modeling using pertinent data including water quality and 2D electrical resistivity imaging data. The observed groundwater depths range from 1.2 to 25.8 meters below ground level with deeper levels observed in upstream and higher elevated regions in the study area. Elevated pH and electrical conductivity (&gt;3000 μS/cm) characterize the Redmud ponds, while other groundwater samples meet BIS drinking water standards. Major ion concentrations contouring indicates high total dissolved solids (TDS) around Redmud ponds due to leakage and in the downstream due to domestic sewage pollution. Most locations adhere to BIS drinking water limits except for one downstream site. The resistivity data indicated hard formations (&gt;18 m depth) with resistivity &gt;500 Ohm.m, weathered basalt (12-18 m, 60-150 Ohm.m), and saturated water (12 m, 1-5 Ohm.m). Contaminated aquifers (&lt;1 Ohm.m) are detected up to 3 m depth, and noticed accumulating contaminated water (27 m depth) in the NW corner of the Redmud pond from seepage. Predominant groundwater flow from the Redmud pond towards public water bodies and streams is highlighted by simulated contaminant transport model. Collaboration and comprehensive management are recommended to protect the watershed’s environmental integrity and public health.

Groundwater inputs could be a significant but often overlooked source of phosphorus in lake ecosystems

Freshwater lakes are severely threatened, due largely to excess inputs of nutrients and other contaminants. Phosphorus (P) is receiving renewed attention due to recent increases in toxic cyanobacteria blooms in lakes worldwide. We investigated groundwater seepage for its role in P loading dynamics at Oneida Lake, New York, USA—one of the most well-studied lakes globally. P loading was measured at representative sites along the 88 km shoreline over three summers by directly measuring groundwater flow using seepage meters and porewater samplers. Groundwater seepage was a continuous and significant source of dissolved P over the summer months, comparable to tributary sources to the lake during that time. This constant input has enriched the concentrations of P in the nearshore surface waters, significantly above levels in the pelagic zone. Pore Total Phosphorus (TP) concentrations and loads reached extremely high values (up to 100 mg/L), with inorganic P representing only ~ 10% of TP per site. Groundwater seepage flows and P loadings were highly variable across space and time, partially explained by adjacent land uses and precipitation. Our research concludes that groundwater seepage is a significant, but overlooked, source of dissolved P and a crucial factor driving summer primary production at Oneida Lake, and likely other temperate lakes.

Optimum identification of aquifer parameter zone structures using the SHuffled Ant Lion Optimization approach considering general form of the Voronoi tessellation

A new simulation–optimization approach is proposed for solving the inverse parameter structure identification problems in groundwater systems. In the simulation part of the proposed approach, the numerical groundwater flow process is simulated by using MODFLOW. Parameter structures of the aquifer system are identified by partitioning the flow domain into a finite number of zones employing the general form of the Voronoi tessellation (VT). This MODFLOW and VT-based simulation model is then integrated into an optimization model where the SHuffled Ant Lion Optimization approach (SHALO) is used. This work is the first application of SHALO by considering the general form of the VT zonation approach for solving the inverse groundwater parameter structure identification problems. The applicability of the proposed simulation–optimization approach is evaluated on a hypothetical aquifer model in the literature. This evaluation is conducted by solving the same problem using Ant Lion Optimization (ALO) and its improved version, the self-adaptive ALO (saALO) for the same conditions. Furthermore, the identified results are compared with those obtained using harmony search (HS) and hybrid genetic algorithm (hybrid GA) based solution approaches in the literature. The obtained results indicated that the proposed approach provided better identification results than ALO, saALO, HS, and hybrid GA in terms of different evaluation metrics.

Influence of precipitation infiltration recharge on hydrological processes of the karst aquifer system and adjacent river

Understanding the hydrological processes of the karst aquifer system (KAS) and adjacent river can further accurately evaluate karst water storage. However, the mechanism of the KAS and adjacent river hydrological processes under precipitation infiltration recharge conditions, the micro-flow behavior of groundwater in the KAS, and changes in media saturation remain unclear. This study coupled the KAS and adjacent river through the Volume of Fluid (VOF), simulated the flow behavior of groundwater in different medium of the KAS and coupled the free flow process of the adjacent river by using the multiphase Darcy-Brinkman-Stokes (DBS) equation, applied the van Genuchten-Mualem and van Genuchten-modified Mualem equations to characterize the seepage processes in soil. The results indicate that the saturation of different porous media controlled by recharge intensity directly influences the flow behavior and flow paths of the porous media. The main groundwater flow path under low precipitation intensity (b=3) is primarily the lower porous medium, while under moderate precipitation intensity (3<b≤8), the main flow path of groundwater is the karst conduit. During heavy precipitation (8<b≤12), due to the limited capacity of karst conduit, they cannot bear the large amount of groundwater supplied from upper porous medium. The peak conduit saturation in laminar flow is 0.3 lower than that in turbulent flows. Under laminar flow conditions, the conduit drainage leads to a rapid decrease in saturation, making it easier for porous media to recharge to the conduit and discharge towards the adjacent river. The turbulent state conduit maintains a high saturation recession, but the discharge is consistently lower than laminar flow. The variation process line of average flow velocity and interfacial flow velocity are obtained using the DBS method and N-S simulation method respectively. There is a certain discrepancy between the simulation results of the two methods since the N-S simulation scheme dose not consider the water retention characteristics of the porous medium.

Groundwater modelling in SWAT+: Considerations for a realistic baseflow simulation

Hydrological models are valuable tools that support the evaluation and management of water resources. Among catchment scale models, SWAT has been widely used, and a revision of this model, SWAT+, was recently released. Groundwater simulation has been pointed out as one of the main weaknesses of these models, and despite the introduction of some improvements in the new version, an accurate simulation of groundwater flow is still a pending task. One of the most common issues is the inability of the model to maintain baseflow during long dry periods. Baseflow is common in aquifer areas and has great relevance for water quality and ecosystems. Although some authors proposed different solutions to its inaccurate simulation (e.g., editing aquifers configuration, coupling with other models, etc.) this work aimed to determine if a realistic groundwater simulation is possible using SWAT+ and determine the reasons why this is not being achieved. For this purpose, a groundwater dominated catchment was modelled, and different calibration approaches were carried out and compared, including strategies such as hard calibration, soft calibration, or sensitivity analysis. A comprehensive analysis of model outputs was crucial to achieve a model with a realistic simulation of groundwater contribution to the streamflow, both in amount and timing. The properly parameterization of some parameters (perco, latq_co, lat_ttime, and especially alpha) was key. Limitations of hard calibration were evidenced, such as the necessity of complementing this strategy with an evaluation of the simulation of hydrological processes. For those parameters governing streamflow components simulation, a detailed description of its functioning is included in this manuscript, together with some guidelines to achieve a statistically accurate model with a realistic groundwater flow simulation. Follow these recommendations will be useful for SWAT users, and will serve to create more reliable models.

Elucidating biogeochemical characterization of nitrogen in the vadose zone integrating geochemistry, microorganism, and numerical simulation

A thorough comprehension of nitrogen biogeochemical processes in the vadose zone is crucial for the effective prevention and remediation of soil-groundwater system contamination. Despite the growing research on this subject, the full scope of nitrogen biogeochemical characterization in different geological environments remains poorly understood. This study addresses this knowledge gap by integrating geochemical, microbiological and numerical simulation approaches to gain a deeper insight into nitrogen biogeochemistry in agriculture. Our findings indicate the biogeochemical behavior of nitrogen in the vadose zone is mediated by microorganisms, driven by hydraulics, influenced by geological conditions and environmental factors. Along the groundwater flow, NH4+-N was found to be heavily accumulated in the topsoil of 0–40 cm, while NO3−-N was transported and driven by hydrodynamics from both vertical and horizontal directions. Microbial diversity, species composition and functional microorganisms were significantly influenced by soil depth, rather than geomorphological types. Oxidation-reduction potential (ORP), total organic carbon (TOC), soil moisture (MOI), bicarbonate (HCO3−), and ferrous (Fe2+) were identified as the principal environmental factors that regulate nitrogen metabolism and the dominant biochemical processes, encompassing nitrogen fixation, nitrification, and denitrification. Driven by hydrodynamics, NH4+-N, NO2−-N and NO3−-N tend to form distinct biochemical reaction zones in the vertical vadose zone. These areas are dynamic and subject to geomorphologies. It should be noted that NO3−-N can migrate towards groundwater from the clayey sand in the Alluvial Plain, which presents a potential risk of groundwater contamination. The fissure structure of loess may serve as the major transport pathway for groundwater nitrogen contamination in the Loess Tableland. This finding highlights the importance of integrating microbiology, geochemistry and hydraulics to elucidate the biogeochemical processes of nitrogen in the vadose zone with a dynamic mindset.

Theoretical analysis of using multiple borehole heat exchangers for production of heating and cooling energy in shallow geothermal reservoirs with underground water flow

In this work the feasibility of using the double-U borehole heat exchangers (BHE) in multiple arrangements for production of heating and cooling energy has been investigated.Theoretical case studies are given for a few hypothetical lithological cases corresponding to the lithology that can be found in City of Zagreb. The results show interference between the borehole heat exchangers and limits of use of shallow reservoir with respect to slow convection and loss of accumulated energy in the reservoir. Special emphasis is placed on the description of the groundwater flow determined by the parameters of the shallow geothermal reservoir, such as permeability, water saturation, thermal conductivity, specific heat capacity and soil density, topography, precipitation, etc. Simulations have been performed in FeFlow coupled with Python user-scripts for individual control of the BHE’s, and the surface equipment is modelled using RES2GEO. To assess the accuracy of numerical predictions in simulating thermal responses, FeFlow software was utilized to conduct Thermal Response Tests (TRT) for individual Borehole Heat Exchangers (BHEs). These numerical results were then juxtaposed with empirical TRT data collected from a site near Velika Gorica, a town within the City of Zagreb. For each case simulation is done for one year with one hour resolution. The temperature distribution analysis highlights notable differences between scenarios, where the temperature field with lower permeability experiences significant variations, particularly near the BHE. Here, temperatures can soar to 20 °C during cooling periods and plummet to 4 °C in heating seasons. These fluctuations cause overheating and subcooling of the reservoir, impairing the heat pump’s efficiency and necessitating increased electricity consumption. By using Multivariate linear regression and Multi-layer perceptron regressor method, the main parameter for determination of COP of ground source heat pump. The analysis reveals that the majority of errors fall within the range of ± 5 %, although the MLA prediction method exhibits slightly more frequent errors. Additionally, the R2 scores demonstrate strong performance: the MLA procedure achieves 0.95 for heating and 0.99 for cooling, while the MLP procedure attains 0.99 for both heating and cooling phases.