Fluid flow and heat transfer in porous media have been extensively studied due to their importance in numerous industrial and environmental applications. This review provides a comprehensive analysis of numerical and experimental approaches, presenting a multiscale perspective that bridges molecular, pore, and macroscopic levels. This study emphasizes the importance of understanding the underlying principles governing these processes, as this knowledge is essential for optimizing and innovating applications ranging from energy systems to environmental engineering. The review synthesizes key theoretical frameworks, including Darcy’s law, the Brinkman equation, and volume-averaging methods, offering a robust foundation for interpreting complex interactions in porous media. A novel aspect of this work is the integration of experimental and numerical insights to address challenges such as heterogeneity, anisotropy, and scale effects, demonstrating their complementary roles in advancing this field. Additionally, the review highlights emerging methodologies, including advanced pore-scale modeling, the lattice Boltzmann method, and machine learning, as transformative tools for overcoming existing limitations and exploring future directions. By identifying critical knowledge gaps and proposing innovative solutions, this article serves as a vital resource for researchers and practitioners, fostering interdisciplinary approaches and paving the way for cutting-edge advancements in the study of fluid flow and heat transfer in porous media.

Contamination of groundwater by naturally occurring arsenic has recently become a disturbing environmental problem in Viterbo area, Central Italy. Arsenic concentrations in most of the public supply networks exceed the maximum allowable limit of 10 μg/l (WHO) for drinking water. The primary purpose of this paper is to obtain a better understanding of the factors contributing to the high levels of As in water supply networks. This study focuses on (a) the determination of basic hydrochemical characteristics of groundwater, (b) the identification of the major sources and processes controlling the As contamination in public supply networks, (c) to find out possible relationships among the As and other trace elements through principal component analysis (PCA). Groundwater samples from public water supply wells and springs were collected and analysed for physico-chemical parameters and trace elements. Springs and well water samples are predominantly of the Na–HCO3, Na –Ca–HCO3 and Ca–HCO3 types and the highest arsenic concentrations were observed in Na–HCO3 type water. Eh-pH diagrams reveal that H2AsO4− and HAsO42−, As(V) arsenate, are the dominating As species highlighting slightly to moderately oxidizing conditions. Geochemical modeling indicates that arsenic-bearing phases were undersaturated in the groundwater, however most of the samples were saturated with respect to Fe (i.e. magnetite, hematite and goethite) and Al (diaspore and boehmite) oxide and hydroxide minerals. Concentrations of As, Li, B, Co, Sr, Mo, U and Se are highly correlated (r > 0.7) with each other, however in some groundwater samples As show also good correlations (r > 0.5) with Fe and Mn elements reflecting the relationships among the trace elements result from different geochemical processes. Evaluation of the principal component (PCA) analysis and geochemical modeling suggest that the occurrence of As and other trace element concentrations in groundwater are probably derived from (i) weathering and/or dissolution of volcanic source aquifer materials and (ii) adsorption/desorption processes on the Fe and Al oxide and hydroxide minerals.

In southern Latium region, Central Italy, groundwater and spring water resources in the carbonate aquifers are the major contributors of drinking and irrigation water supply. The aim of this study was to review hydrochemical processes that control the groundwater chemistry and to determine the suitability of springs and groundwater for irrigation and drinking purposes on the basis of the water quality indices. Physical (pH, electrical conductivity, total dissolved solids) and hydrochemical characteristics (Na⁺, K⁺, Ca²⁺, Mg²⁺, HCO₃⁻, Cl⁻, and SO₄⁻) of springs and groundwater were determined. To assess the water quality, chemical parameters like sodium adsorption ratio (SAR), total hardness, Mg-hazard (MH), sodium percentage (Na %), salinity hazard, permeability index, and Kelly’s ratio were calculated based on the analytical results. A Durov diagram plot revealed that the groundwater has been evolved from Ca to HCO₃ recharge water, followed by mixing and reverse ion exchange processes, due to the respective dominance of Na–Cl and Ca–Cl water types. According to Gibbs’s diagram plots, chemical weathering of rock forming minerals is the major driving force controlling water chemistry in this area. Groundwater and spring samples were grouped into six categories according to irrigation water quality assessment diagram of US Salinity Laboratory classification and most of the water samples distributed in category C2–S1 and C3–S1 highlighting medium to high salinity hazard and low sodium content class. The results of hydrochemical analyses and the calculated water quality parameters suggest that most of the water samples are suitable for irrigation and drinking purposes, except for the samples influenced by seawater and enhanced water–rock interaction. High values of salinity, Na %, SAR, and MH at certain sites, restrict the suitability for agricultural uses.

The increasing occurrence of widespread drought phenomena is a global environmental emergency, especially for the effects of ongoing climate change on groundwater availability. Dry years and extreme temperatures are common drivers of current climate impacts all over the world, including, for example, the freshwater supply for drinking and agriculture purposes, ecosystems, forestry, health, etc. In this frame, to ensure temporal water availability in water-stressed areas, sustainable groundwater management is an increasing challenge. Most groundwater in the South-East Latium Region, Central Italy, as in the whole Apennine Mountains chain, is stored in karst aquifers. In this area important water resources are present, but even here in the last decades they are affected by groundwater depletion as a consequence of occurring drought events, the upward trend in the global average temperature and the increase in groundwater extraction. Due to the lack of flow rate data of springs in many areas of Italy, spring response modeling could be a useful tool for supporting proper water resource management. Several research studies proposed methods based on relationships between spring discharges and rainfall data. The goal of this paper is to propose a simple linear model, based on rainfall-discharge cross correlations, in order to be able to simulate the effects of variation in rainfall on water resource connected Capodacqua di Spigno Spring minimum discharge. The results obtained using the developed model has been compared to an existing method that uses the Standard Precipitation Index (SPI) for the estimation of the minimum annual spring discharge.

Karst aquifer recharge areas are usually difficult to identify because of the complexity of these aquifers’ characteristics. On the other hand, their identification is very important in the aim of protecting the groundwater resources that these aquifers host. Regarding this topic, this paper presents an approach aimed at identifying karst aquifer recharge areas by the application of oxygen-18 and deuterium isotopes composition of groundwater coupled with hydrological features. Oxygen-18 and deuterium isotope composition of Capodacqua di Spigno Spring, in the South of the Latium Region, has been applied with rainfall and discharge values related to the feeding aquifer of this spring. As δ¹⁸O and δ²H values of groundwater samples are natural tracers of the recharge area’s elevation, we propose a model, based on the distribution of the basin surfaces involved as recharge areas, in relation to elevations. The model estimates, for any discharge value, the percentage of the topographic area involved in the aquifer recharge. The setting up of this simulated distribution is supported by a Weibull cumulative probability function. The results show that the measured discharges increase as larger areas with lower elevations are involved in the recharge process.

Water resources management is one of the most important challenges worldwide because water represents a vital resource for sustaining life and the environment. With the aim of sustainable groundwater management, the identification of aquifer recharge areas is a useful tool for water resources protection. In a well-developed karst aquifer, environmental isotopes provide support for identifying aquifer recharge areas, residence time and interconnections between aquifer systems. This study deals with the use of environmental isotopes to identify the main recharge area of a karst aquifer in the Upper Valley of Aniene River (Central Italy). The analysis of ¹⁸O/¹⁶O and ²H/H values and their spatial distribution make it possible to trace back groundwater recharge areas based on average isotope elevations. The Inverse Hydrogeological Balance Method was used to validate spring recharge elevations obtained by the use of stable isotopes. Areas impacted by direct and rapid rainfall recharge into the study area were delineated, showing groundwater flowpaths from the boundaries to the core of the aquifer. The results of this study demonstrate the contribution that spatial and temporal isotope changes can provide to the identification of groundwater flowpaths in a karst basin, taking into account the hydrogeological setting.