Groundwater-Surface Water Interaction in the Nera River Basin (Central Italy): New Insights after the 2016 Seismic Sequence

The highest part of the Nera River basin (Central Italy) hosts significant water resources for drinking, hydroelectric, and aquaculture purposes. The river is fed by fractured large carbonate aquifers interconnected by Jurassic and Quaternary normal faults in an area characterized by high seismicity. The October 30, 2016, seismic sequence in Central Italy produced an abrupt increase in river discharge, which lasted for several months. The analysis of the recession curves well documented the processes occurring within the basal aquifer feeding the Nera River. In detail, a straight line has described the river discharge during the two years after the 2016 seismic sequence, indicating that a turbulent flow characterized the emptying process of the hydrogeological system. A permeability enhancement of the aquifer feeding the Nera River—due to cleaning of fractures and the co-seismic fracturing in the recharge area—coupled with an increase in groundwater flow velocity can explain this process. The most recent recession curves (2019 and 2020 periods) fit very well with the pre-seismic ones, indicating that after two years from the mainshock, the recession process recovered to the same pre-earthquake conditions (laminar flow). This behavior makes the hydrogeological system less vulnerable to prolonged droughts, the frequency and length of which are increasingly affecting the Apennine area of Central Italy.

The use of Radon (Rn222) isotopes to detect groundwater discharge in streams draining Table Mountain Group (TMG) aquifers

Environmental isotopes have been used for decades as natural tracers in studies aimed at understanding complex hydrogeological processes such as groundwater and surface water interactions. Radon (Rn222) is a naturally occurring, radioactive isotope which is produced from radium (Ra226) during the radioactive decay series of uranium (U238). Since U238 is present in most geological substrates, Rn222 is produced in various lithological structures and subsequently transported with groundwater through fractures and pore spaces in an aquifer towards surface water discharge points in rivers and springs. This study aimed to determine (i) the concentration of Rn222 within both surface water and groundwater in Table Mountain Group (TMG) aquifer systems, and (ii) the feasibility of using Rn222 isotopes as a natural tracer in groundwater-surface water interaction studies. This study was conducted in a highly fractured TMG aquifer system near Rawsonville, South Africa. Surface water from two perennial rivers (i.e. Gevonden and Molenaars), together with groundwater from a nearby borehole, were sampled and their corresponding Rn222 concentrations measured. Our study found median Rn222 concentrations in the Gevonden River of 76.4 Bq∙L-1 and 67.2 Bq∙L-1 in the dry and wet seasons, respectively. Nearly 12% of surface water samples exceeded 100 Bq∙L-1. These abnormally high Rn222 concentrations can only be attributed to the influx of groundwater with extremely high Rn222 concentrations. Under ambient (no pumping) conditions, Rn222 concentrations in groundwater range between 130–270 Bq∙L-1. However, when the borehole was pumped, and inflowing water from the surrounding aquifer was sampled, even higher Rn222 concentrations (391–593 Bq∙L-1) were measured. These extremely high Rn222 concentrations in groundwater are believed to be attributed to the underlying granitic geology and the prevalence of faults. The use of Rn222 isotopes as an environmental tracer in groundwater–surface water interaction studies is therefore regarded as a feasible option in similar highly fractured aquifer systems.

Modelling of Groundwater–Surface Water Interaction Applying the Hyporheic Flux Model

The objective of the present paper is to develop a methodology that could allow the representation of the analytical hyporheic flux equation model (AHF) in a numerical model done in MODFLOW. Therefore, the scope of the research is to show the viability of the methodology suggested in a real case (Biebrza river, Poland, Europe). Considering that the model requires extensive manipulation in the creation of the packages, a test phase through the seepage package of MODFLOW is carried out with the aim of representing the river package of MODFLOW. FloPy is the tool chosen to develop this implementation due to the versatility of manipulating the packages available in MODFLOW through coding. The obtained results showed a correct implementation of the AHF model using the example of the Biebrza River. The results obtained will enable a better understanding regarding the modelling of the interaction between the river and the aquifer, considering streams with specific geometries where the depth is dimensionally higher than the width.

Crenal Habitats: Sources of Water Mite (Acari: Hydrachnidia) Diversity

Many studies emphasized the role that water mites play within the invertebrate communities of spring ecosystems, regarding species diversity and its significance within the crenal food web, as well as the specific preferences water mites exhibit towards spring typology. In pristine natural springs with permanent flow, water mites are nearly always present and usually display high diversity. This study aimed to determine whether significant differences in water mite assemblages between rheocrene (river-forming springs with dominant riffle habitats) and limnocrene (lake-forming springs with dominant pool habitats) karst springs could be detected in terms of species richness, diversity and abundance, but also in different ratios of specific synecological groups: crenobiont (exclusively found in springs), crenophilous (associated with springs) and stygophilous (associated with groundwater) water mite taxa. Our research was carried out on four limnocrenes and four rheocrenes in the Dinaric karst region of Croatia. Seasonal samples (20 sub-samples per sampling) were taken at each spring with a 200-µm net, taking into consideration all microhabitat types with coverage of at least 5%. Water mite abundance was found not to differ between morphological spring types. Significantly higher values of species richness and diversity indices were found in rheocrenes compared to limnocrenes, like those usually reported for this type of springs. However, unlike the studies previously reported, in this case, the higher shares of crenophilous and crenobiont water mite individuals were found in limnocrenes. The differences between stygophilous water mite taxa ratios among spring morphotypes were not statistically significant, indicating that the degree of the groundwater/surface water interaction (and water mite interaction therein) does not seem to be directly influenced by spring morphotype. Within this research, 40% of identified water mite species (eight out of twenty) were recorded for the first time in Croatia, thus highlighting again a huge gap in water mite knowledge of karst springs.

Groundwater–Surface Water Interaction—Analytical Approach

Modelling of water flow in the hyporheic zone and calculations of water exchange between groundwater and surface waters are important issues in modern environmental research. The article presents the Analytical Hyporheic Flux approach (AHF) permitting calculation of the amount of water exchange in the hyporheic zone, including vertical water seepage through the streambed and horizontal seepage through river banks. The outcome of the model, namely water fluxes, is compared with the corresponding results from the numerical model SEEP2D and simple Darcy-type model. The errors of the AHF model, in a range of 11–16%, depend on the aspect ratio of water depth to river width, and the direction of the river–aquifer water exchange, i.e., drainage or infiltration. The AHF model errors are significantly lower compared to the often-used model based on vertical water seepage through the streambed described by Darcy’s law.

Analytical and Numerical Groundwater Flow Solutions for the FEMME-Modeling Environment

Simple analytical and numerical solutions for confined and unconfined groundwater-surface water interaction in one and two dimensions were developed in the STRIVE package (stream river ecosystem) as part of FEMME (flexible environment for mathematically modelling the environment). Analytical and numerical solutions for interaction between one-dimensional confined and unconfined aquifers and rivers were used to study the effects of a 0.5 m sudden rise in the river water level for 24 h. Furthermore, a two-dimensional groundwater model for an unconfined aquifer was developed and coupled with a one-dimensional hydrodynamic model. This model was applied on a 1 km long reach of the Aa River, Belgium. Two different types of river water level conditions were tested. A MODFLOW model was set up for these different types of water level condition in order to compare the results with the models implemented in STRIVE. The results of the analytical solutions for confined and unconfined aquifers were in good agreement with the numerical results. The results of the two-dimensional groundwater model developed in STRIVE also showed that there is a good agreement with the MODFLOW solutions. It is concluded that the facilities of STRIVE can be used to improve the understanding of groundwater-surface water interaction and to couple the groundwater module with other modules developed for STRIVE. With these new models STRIVE proves to be a powerful example as a development and testing environment for integrated water modeling.