Agua Amarga coastal aquifer in southern Spain has been the subject of chemical and physical measurements since May 2008 in order to monitor the potential effects of water withdrawal for the Alicante desalination plants on the salt marsh linked to the aquifer. Electrical conductivity contour maps and depth profiles, piezometric-head contour maps, hydrochemical analyses, isotopic characterizations and temperature depth profiles show not only the saltwater intrusion caused by water abstraction, but also the presence of a pronounced convective density-driven flow below the salt marsh; this flow was a consequence of saltwork activity in the early 1900s which generated saline groundwater contamination. The influence of a seawater recharge programme, carried out over the salt marsh in 2009–2010, on the diminishing groundwater salinity and the recovery of groundwater levels is also studied. Based on collected field data, the project provides a deeper understanding of how these successive anthropogenic interventions have modified flow and mixing processes in Agua Amarga aquifer.

Analysis of geothermal and hydrogeologic characteristics of Paleozoic interior basins of Brazil has identified an association between the geographic distribution of thermal springs and areas of occurrences of groundwater flow. Specifically, thermal springs are found to be absent in regions inferred to have lateral flows of groundwater. This trend is evident in the basins of the Amazon region, in the central parts of the Parnaíba basin and in the west-central parts of the Paraná basin. Model studies help to elucidate mutual exclusion of regions of thermal springs and sub-horizontal flows of groundwater. Numerical simulations indicate that groundwater flows with velocities ≥ 1 cm/year are capable of masking the occurrence of thermal anomalies. Also, down flow through distributed recharge zones can lead to development of large zones of relatively low temperature. The observational data sets of temperature gradients and Peclet numbers have been employed outlining advection-convection domains of subsurface strata in the sedimentary basins of the Amazon region, Parnaíba and Paraná. Results obtained indicate that thermal buoyancy forces are incapable of overcoming advective flows in basins of the Amazon region. Similar conditions are also found to prevail in the central parts of the Parnaíba and Paraná basins.

Among the advances made in analytical and numerical analysis methods to quantify groundwater/surface-water interaction, one methodology that stands out is the use of heat as an environmental tracer. A large data set of river and riverbed temperature profiles from the Aa River in Belgium has been used to examine the spatial-temporal variations of groundwater/surface-water interaction. Exchange fluxes were calculated with the numerical heat-transport code STRIVE. The code was applied in transient mode to overcome previous limitations of steady-state analysis, and allowed for the calculation of model quality. In autumn and winter the mean exchange fluxes reached −90 mm d⁻¹, while in spring and early summer fluxes were −42 mm d⁻¹. Predominantly gaining conditions occurred along the river reach; however, in a few areas the direction of flow changed in time. The river banks showed elevated fluxes up to a factor of 3 compared to the center of the river. Higher fluxes were detected in the upstream section of the reach. Due to the influence of exchange fluxes along the river banks, larger temporal variations were found in the downstream section. The exchange fluxes at the river banks seemed more driven by variable local exchange flows, while the center of the river was dominated by deep and steady regional groundwater flows. These spatial and temporal differences in groundwater/surface-water exchange show the importance of long-term investigations on the driving forces of hyporheic processes across different scales.

Investigating the interaction of groundwater and surface water is key to understanding the hyporheic processes. The vertical water fluxes through a streambed were determined using Darcian flux calculations and vertical sediment temperature profiles to assess the pattern and magnitude of groundwater/surface-water interaction in Beiluo River, China. Field measurements were taken in January 2015 at three different stream morphologies including a meander bend, an anabranching channel and a straight stream channel. Despite the differences of flux direction and magnitude, flux directions based on vertical temperature profiles are in good agreement with results from Darcian flux calculations at the anabranching channel, and the Kruskal-Wallis tests show no significant differences between the estimated upward fluxes based on the two methods at each site. Also, the upward fluxes based on the two methods show similar spatial distributions on the streambed, indicating (1) that higher water fluxes at the meander bend occur from the center of the channel towards the erosional bank, (2) that water fluxes at the anabranching channel are higher near the erosional bank and in the center of the channel, and (3) that in the straight channel, higher water fluxes appear from the center of the channel towards the depositional bank. It is noted that higher fluxes generally occur at certain locations with higher streambed vertical hydraulic conductivity (K ᵥ) or where a higher vertical hydraulic gradient is observed. Moreover, differences of grain size, induced by stream morphology and contrasting erosional and depositional conditions, have significant effects on streambed K ᵥ and water fluxes.

The thermal profile of a streambed is affected by a number of factors including: temperatures of stream water and groundwater, hydraulic conductivity, thermal conductivity, heat capacity of the streambed, and the geometry of hyporheic flow paths. Changes in these parameters over time cause changes in thermal profiles. In this study, temperature data were collected at depths of 30, 60, 90 and 150 cm at six streambed wells 5 m apart along the thalweg of Little Kickapoo Creek, in rural central Illinois, USA. This is a third-order low-gradient baseflow-fed stream. A positive temperature gradient with inflection at 90-cm depth was observed during the summer period. A negative temperature gradient with inflection at 30 cm was observed during the winter period, which suggests greater influence of stream-water temperatures in the substrate during the summer. Thermal models of the streambed were built using VS2DHI to simulate the thermal profiles observed in the field. Comparison of the parameters along with analysis of temperature envelopes and Peclet numbers suggested greater upwelling and stability in temperatures during the winter than during the summer. Upwelling was more pronounced in the downstream reach of the pool in the riffle and pool sequence.

A hydrogeological conceptual model of the source, circulation pathways and temporal variation of a low-enthalpy thermal spring in a fractured limestone setting is derived from a multidisciplinary approach. St. Gorman’s Well is a thermal spring in east-central Ireland with a complex and variable temperature profile (maximum of 21.8 °C). Geophysical data from a three-dimensional(3D)audio-magnetotelluric(AMT) survey are combined with time-lapse hydrogeological data and information from a previously published hydrochemical analysis to investigate the operation of this intriguing hydrothermal system. Hydrochemical analysis and time-lapse measurements suggest that the thermal waters flow within the fractured limestones of the Carboniferous Dublin Basin at all times but display variability in discharge and temperature. The 3D electrical resistivity model of the subsurface revealed two prominent structures: (1) a NW-aligned faulted contact between two limestone lithologies; and (2) a dissolutionally enhanced, N-aligned, fault of probable Cenozoic age. The intersection of these two structures, which has allowed for karstification of the limestone bedrock, has created conduits facilitating the operation of relatively deep hydrothermal circulation (likely estimated depths between 240 and 1,000 m) within the limestone succession of the Dublin Basin. The results of this study support a hypothesis that the maximum temperature and simultaneous increased discharge observed at St. Gorman’s Well each winter is the result of rapid infiltration, heating and recirculation of meteoric waters within a structurally controlled hydrothermal circulation system.

Temperature distribution and heat transport are studied in a coastal aquifer at De Panne in the western Belgian coastal plain. Field observations include temperature profiles of groundwater in the dunes and temperature measurements at the water table in a profile on the shore. Freshwater–saltwater distribution is known from previous studies. These are used to constrain a density-dependent model simulating the freshwater–saltwater distribution and heat transport using the SEAWAT code. The yearly fluctuation of the groundwater temperature in the phreatic aquifer under the dunes, shore and sea, and the influence of a tidal inlet in the dunes are simulated. The observations show that seawater temperature variations determine the temperature variations on the shore whereas atmospheric temperature changes determine this in the dunes. Yearly temperature fluctuations imposed at the water table propagate mainly vertically in the aquifer with only limited lateral influence. Heat transport is mainly convection dominated. Thickness of the surficial zone is determined by the amplitude of the groundwater temperature at the water table and the groundwater flow. Establishment of a tidal inlet in the dunes results in asymmetric temperature profiles under and in the vicinity of it.