Biostimulation and Bioaugmentation to Enhance Reductive Dechlorination of TCE in a Long‐Term Flow Through Column Study

Large laboratory columns (15.2 cm diameter, 183 cm long) were fed with groundwater containing trichloroethylene (TCE), were biostimulated and bioaugmented, and were monitored for over 7.5 years. The objective of the study was to observe how the selection of the carbon and energy source, i.e., whey, Newman Zone^® standard surfactant emulsified oil and Newman Zone nonionic surfactant emulsified oil, affected the rate and extent of dechlorination. Column effluent was monitored for TCE and its degradation products, redox indicators (nitrate‐N, Fe(II), sulfate), and changes in iron mineralogy. Total bacteria and Dehalococcoides mccartyi strains were quantified using q‐PCR. Complete dechlorination was only observed in the whey treated columns, occurring 1 year after bioaugmentation with addition of a culture known to dechlorinate TCE to ethene, and 3 years later in the non‐bioaugmented column. The addition of the emulsified oils with or without bioaugmentation resulted in dechlorination only through cis‐DCE and vinyl chloride. While Dehalococcoides mccartyi strains are the only known bacteria that can fully dechlorinate TCE, their presence, either natural or augmented, was not the sole determiner of complete dechlorination. The establishment of a supporting microbial community and biogeochemistry that developed with continuous feeding of whey, in addition to the presence of D. mccartyi, were necessary to support complete reductive dechlorination. Results confirm that careful selection of a biostimulant is critical to the success of TCE dechlorination in complex soil environments.     

Simultaneous Organic Matter and Nitrogen Removal from Piggery Slaughterhouse Wastewater Using Bioaugmented Gaslift Membrane Bioreactor

A single bioaugmentation reactor and a side-stream gaslift membrane bioreactor combined with bioaugmentation are conducted to treat real wastewater from a centralized piggery slaughterhouse in Vietnam. The bioaugmented reactor is inoculated with heterotrophic microorganisms (Bacillus sp.) isolated from piggery slaughterhouse wastewater. The results of a single bioaugmentation experiment show high removal efficiency of chemical oxygen demand (COD) (84.8%–97.5%) and total nitrogen (TN) (69.9%–87.2%) at loading rates of 1.28–3.89 and 0.14–0.37 kg m⁻³ d⁻¹, respectively. The combined system demonstrates a significantly higher TN removal efficiency (89.0%–96.1%) (p < 0.001), more stable flux (36.0–38.4 L m⁻² h⁻¹), and transmembrane pressure (0.95–1.05 bar), and better capacity of separation of solid–liquid phases compared to the single bioaugmentation. High COD and TN removal efficiency is possibly due to assimilation and simultaneous nitrification and denitrification processes. The results of this study also indicate the feasibility and propitious efficiency of the bioaugmented gaslift membrane bioreactor for piggery slaughterhouse wastewater treatment.     

Bioaugmentation and Propane Biosparging for In Situ Biodegradation of 1,4‐Dioxane

Propane biosparging and bioaugmentation were applied to promote in situ biodegradation of 1,4‐dioxane at Site 24, Vandenberg Air Force Base (VAFB), CA. Laboratory microcosm and enrichment culture testing demonstrated that although native propanotrophs appeared abundant in the shallow water‐bearing unit of the aquifer (8 to 23 ft below ground surface [bgs]), they were difficult to be enriched from a deeper water‐bearing unit (82 to 90 feet bgs). Bioaugmentation with the propanotroph Rhodococcus ruber ENV425, however, supported 1,4‐dioxane biodegradation in microcosms constructed with samples from the deep aquifer. For field testing, a propane‐biosparging system consisting of a single sparging well and four performance monitoring wells was constructed in the deep aquifer. 1,4‐dioxane biodegradation began immediately after bioaugmentation with R. ruber ENV425 (36 L; 4 × 10⁹ cells/mL), and apparent first‐order decay rates for 1,4‐dioxane ranged from 0.021 day^?1 to 0.036 day^?1. First‐order propane consumption rates increased from 0.01 to 0.05 min^?1 during treatment. 1,4‐dioxane concentrations in the sparging well and two of the performance monitoring wells were reduced from as high as 1090 µg/L to <2 µg/L, while 1,4‐dioxane concentration was reduced from 135 µg/L to 7.3 µg/L in a more distal third monitoring well. No 1,4‐dioxane degradation was observed in the intermediate aquifer control well even though propane and oxygen were present. The demonstration showed that propane biosparging and bioaugmentation can be used for in situ treatment of 1,4‐dioxane to regulatory levels.     

Modeling Dehalococcoides sp. Augmented Bioremediation in a Single Fracture System

A numerical reactive transport model was developed to simulate the bioremediation processes in a perchloroethene (PCE) contaminated single fracture system augmented with Dehalococcoides sp. (DHC). The model describes multispecies bioreactive transport processes that include bacterial growth and detachment dynamics, biodegradation of chlorinated species, competitive inhibition of various reactive species, and the loss of daughter products because of back‐partitioning effects. Two sets of experimental data, available in the study by Schaefer et al. (2010b) , were used to calibrate and test the model. The model was able to simulate both datasets. The simulation results indicated that the yield coefficient and the DHC maximum utilization rate coefficient were the two important process parameters. A detailed sensitivity study was completed to quantify the sensitivity of the model to variations in these two parameter values. The results show that an increase in yield coefficient increases bacterial growth and thus expedites the dechlorination process, whereas an increase in maximum utilization rate coefficient greatly increased dechlorination rates. The proposed model provides a mathematical framework for simulating remediation systems that employ DHC bioaugmentation for restoring chlorinated‐solvent contaminated groundwater aquifers.     

Determination of aquifer parameters based on measurements of tidal effects on a coastal aquifer near Beihai,China

Measurements of the tide and groundwater levels in coastal zones are of importance in determining the properties of coastal aquifers. The solution to a one‐dimensional unsteady groundwater flow model in a coastal confined aquifer with sinusoidal fluctuation of the tide shows that the tidal efficiency decreases exponentially with distance and the time lag increases linearly with distance from the coast. The aquifer property described by the ratio of storage coefficient to transmissivity is determined if the damping constant of the tidal efficiency or the slope of the time lag with distance are obtained on the basis of tidal measurements. Hourly observations of the tide and groundwater levels at 10 wells on the northern coast near Beihai, China show that with distance from the coast, tidal efficiency decreases roughly exponentially and the time lag increases roughly linearly. The estimated ratio of storage coefficient to transmissivity of the confined aquifer ranges from 1·169 × 10^?6 d m^?2 to 1·83 × 10^?7 d m^?2. For a given transmissivity of 750 m² d^?1, the storage coefficient of the aquifer is 8·7675 × 10^?4 with the tidal efficiency method and 1·3725 × 10^?4 with the time lag method. The damping constant of the tidal efficiency with distance can be defined as the tidal propagation coefficient. The value of the confined aquifer is determined as 0·0018892 m^?1. Copyright © 2007 John Wiley & Sons, Ltd.     

Microcosm Assessment of Polaromonas sp. JS666 as a Bioaugmentation Agent for Degradation of cis-1,2-dichloroethene in Aerobic,Subsurface Environments

Chlorinated ethenes such as tetrachloroethene and trichloroethene have been widely used as dry-cleaning and degreasing solvents. Under anaerobic conditions, microorganisms reduce these parent compounds to less-chlorinated daughter products such as cis-1,2-dichloroethene (cDCE), and often further to ethene. This process can be stalled at cDCE, due to insufficient supply of reductants and/or inadequate microbial-community composition. Recently, a novel bacterium, Polaromonas sp. JS666, was isolated that is able to aerobically oxidize cDCE as sole carbon and energy source. As such, it is a promising candidate for use as a subsurface, bioaugmentation agent at sites where anaerobic bioremediation is inappropriate or has stalled and cDCE has migrated to, and accumulated within, aerobic zones, or where it is practical to impose aerobic conditions. Subsurface sediments or groundwater from six such cDCE-contaminated sites were used to construct microcosms. In every sediment or groundwater inoculated with JS666, the organism was able to degrade cDCE, provided that the pH remained circum-neutral. Even when JS666 was challenged with an alternate carbon source, or in the presence of competitive/predatory microorganisms, there was a measure of success. Collectively, these microcosm studies suggest that JS666 is a viable candidate for the bioaugmentation of aerobic, cDCE-contaminated sites. A minimum inoculation level in excess of 10⁵ cells per mL is recommended for field applications. At this level of inoculation, 100 L of inoculum culture grown to an OD600 of 1.0 should be able to treat a 10-m × 30-m × 80-m (24,000-m³) plot.     

Effect of variable‐density groundwater flow on nitrate flux to coastal waters

A two‐dimensional variable‐density groundwater flow and transport model was developed to provide a conceptual understanding of past and future conditions of nitrate (NO₃) transport and estimate groundwater nitrate flux to the Gulf of Mexico. Simulation results show that contaminant discharge to the coast decreases as the extent of saltwater intrusion increases. Other natural and/or artificial surface waters such as navigation channels may serve as major sinks for contaminant loading and act to alter expected transport pathways discharging contaminants to other areas. Concentrations of NO₃ in the saturated zone were estimated to range between 30 and 160 mg?L^?1 as NO₃. Relatively high hydraulic vertical gradients and mixing likely play a significant role in the transport processes, enhancing dilution and contaminant migration to depth. Residence times of NO₃ in the deeper aquifers vary from 100 (locally) to about 300 years through the investigated aquifer system. NO₃ mass fluxes from the shallow aquifers (0 to 5.7 × 10⁴ mg?m^?2?day^?1) were primarily directed towards the navigation channel, which intersects and captures a portion of the shallow groundwater flow/discharge. Direct NO₃ discharge to the sea (i.e. Gulf of Mexico) from the shallow aquifer was very low (0 to 9.0 × 10¹ mg · m^?2?day^?1) compared with discharge from the deeper aquifer system (0 to 8.2 × 10³ mg?m^?2?day^?1). Both model‐calibrated and radiocarbon tracer‐determined contaminant flux estimates reveal similar discharge trends, validating the use of the model for density‐dependent flow conditions. The modelling approach shows promise to evaluate contaminant and nutrient loading for similar coastal regions worldwide. Copyright © 2015 John Wiley & Sons, Ltd.     

Bacterial and Virus‐Like Particle Abundances in Purged and Unpurged Groundwater Depth Profiles

Bacteria and viruses are ubiquitous in subterranean aquatic habitats. Bacterial abundance is known to vary with depth in aquifers; however, whether viral abundance varies with depth is less well known. Here we use flow cytometry (FCM) to enumerate bacteria and virus‐like particles (VLP) from groundwater depth profiles. Groundwater samples were obtained from a set of nested piezometers from depths of 15, 30, 45, 60, 80, and 90 m and bacteria and VLP abundances were determined in purged aquifer water and unpurged water at each slot depth. Mean bacterial abundance (cells / mL) was not significantly different in unpurged water (3.2 × 10⁵) compared to purged water (1.4 × 10⁵); however, mean VLP abundance (particles / mL) was significantly greater in unpurged water (4.4 × 10⁵) compared to purged water (2.3 × 10⁵). Purged water was used to investigate the aquifer depth profile and bacterial and VLP abundances were observed to vary significantly between depths. The virus‐bacteria ratio was determined and was observed to steadily increase with depth. Overall, our data indicate the dynamic nature of bacterial and viral abundances in subsurface environments which should be considered when designing groundwater microbial sampling methodologies.     

Hydrogeology of Wadi Al-Yammaniyah, Saudi Arabia

Quaternary alluvium, ranging in thickness from a few to 100 meters underlain by Precambrian rocks of metamor-phic and igneous origin, constitutes an important source of ground water in Wadi Al-Yammaniyah, Saudi Arabia. The purpose of this report is to assess the hydraulic properties, quality of water, and estimated change in storage in waterbearing rocks in the area. The results of eight pumping tests carried out in hand-dug, large-diameter wells, indicate that the hydraulic conductivity of the alluvial aquifer ranges from 5.6 × 10⁻⁵ to 1.85 × 10⁻³ cm/second (3.36 × 10⁻⁵ to 1.11 × 10⁻³ m/minute) and that its storativity varies from 8.23 × 10⁻² to 1.17 × 10⁻¹. The aquifer is replenished by sporadic but intensive rainfall of short duration. The present withdrawal is only about 10 percent of the annual recharge which is estimated at 52 × 10⁶ m³. It is shown that there is a substantial potential for the future development of potable ground water which would be required for the development of the area.     

The Impact of the Degree of Aquifer Confinement and Anisotropy on Tidal Pulse Propagation

Oceanic tidal fluctuations which propagate long distances up coastal rivers can be exploited to constrain hydraulic properties of riverbank aquifers. These estimates, however, may be sensitive to degree of aquifer confinement and aquifer anisotropy. We analyzed the hydraulic properties of a tidally influenced aquifer along the Meghna River in Bangladesh using: (1) slug tests combined with drilling logs and surface resistivity to estimate Transmissivity (T); (2) a pumping test to estimate T and Storativity (S) and thus Aquifer Diffusivity (D_PT); and (3) the observed reduction in the amplitude and velocity of a tidal pulse to calculate D using the Jacob‐Ferris analytical solution. Average Hydraulic Conductivity (K) and T estimated with slug tests and borehole lithology were 27.3 m/d and 564 m²/d, respectively. Values of T and S determined from the pumping test ranged from 400 to 500 m²/d and 1 to 5 × 10^?4, respectively with D_PT ranging from 9 to 40 × 10⁵ m²/d. In contrast, D estimated from the Jacob‐Ferris model ranged from 0.5 to 9 × 10⁴ m²/d. We hypothesized this error resulted from deviations of the real aquifer conditions from those assumed by the Jacob‐Ferris model. Using a 2D numerical model tidal pulses were simulated across a range of conditions and D was calculated with the Jacob‐Ferris model. Moderately confined (K_top/K_aquifer < 0.01) or anisotropic aquifers (K_x/K_z > 10) yield D within a factor of 2 of the actual value. The order of magnitude difference in D between pumping test and Jacob‐Ferris model at our site argues for little confinement or anisotropy.     

Measuring a Low Horizontal Hydraulic Gradient in a High Transmissivity Aquifer

High transmissivity aquifers typically have low hydraulic gradients (i.e., a flat water table). Measuring low gradients using water levels can be problematic because measurement error may be greater than the true difference in water levels (i.e., a low signal-to-noise ratio). In this study, the feasibility of measuring a hydraulic gradient in the range of 10⁻⁶ to 10⁻⁵ m/m was demonstrated. The study was performed at a site where the depth to water from land surface ranged from 40.1 to 94.2 m and the aquifer transmissivity was estimated at 41,300 m²/d (hydraulic conductivity of 18,800 m/d). The goals of the study were to reduce measurement error as much as practicable and assess the importance of factors affecting water level measurement accuracy. Well verticality was the largest source of error (0.000 to 0.168 m; median of 0.014 m), and geodetic survey of casing elevations was the next most important source of error (0.002 to 0.013 m; median of 0.005 m). Variability due to barometric pressure fluctuations was not an important factor at the site. Hydraulic heads were measured to an accuracy of ±0.0065 m, and the average hydraulic gradient was estimated to be 8.0 × 10⁻⁶ (±0.9 × 10⁻⁶) m/m. The improvement in accuracy allowed for two reversals in the groundwater flow direction to be identified, after which the gradient averaged 2.5 × 10⁻⁵ (±0.4 × 10⁻⁵) m/m. This study showed it is possible to sufficiently control sources of error to measure hydraulic gradients in the 10⁻⁶ to 10⁻⁵ m/m range.     

Leakage properties of the multi-aquifer system underlying the Greater Onitsha water scheme borehole field,Anambra State,Nigeria

Abstract

Leakage properties and the potential for land subsidence due to groundwater withdrawal from a multi-aquifer water supply system were investigated by applying leaky type curve and one dimensional consolidation models to drawdown data that were obtained during a pumping test experiment in an aquifer-aquitard system. The producing aquifer has transmissivity and storativity values of 5.3 × 10^?3 m² s^?1 and 9.54 × 10^?4 respectively. It is recharged through leakage at a rate of 5.67 × 10^?8 m s^?1, giving a leakage amount of more than 0.007 m³ s^?1. Drainage of the aquifer-aquitard system could result in aquitard compaction of between 50 and 180 mm year^?1 for pumping periods of 6 and 22 h day^?1, respectively. The observed leakage has important implications for land subsidence problems and waste disposal practices in the area.     

Aquifer Studies Using Flow Simulations

Flow simulations supplemented resource evaluation studies of the Lower Greensand aquifer. The annual recharge has hitherto been considered to be the resource available for development although a large confined storage exists. Pumping is principally in urban areas but recent monitoring suggests that the aquifer may have suffered excessive drawdowns. Therefore, simulation studies were carried out to refine the annual recharge estimate, to define the past effects of pumping and to obtain an indication of future trends. Catchment water balance indicates average annual recharge of 28 × 10⁶ m³. One-dimensional (ID) ground-water flow simulation suggests that regional abstractions affect outcrop-water levels. The study showed that the aquifer reached a new equilibrium with the current pumping regime. An improved insight to the aquifer has permitted the formulation of strategies for future resource development.     

Numerical studies of vertical Cl−, δ2H and δ18O profiles in the aquifer–aquitard system in the Pearl River Delta,China

Sedimentation may have a significant effect on the transport of solutes and environmental isotopes in sediment. The depth profiles of the Cl^?, δ²H and δ¹⁸O in a borehole in the aquifer–aquitard system in the Pearl River Delta (PRD), China, were obtained by centrifuging the core sediment samples. A one‐dimensional model based on the sedimentation and sea level changes of the PRD during the Holocene was built to investigate numerically the transport mechanisms of Cl^?, δ²H and δ¹⁸O. The sedimentation process was modelled as a moving boundary problem with the moving rate equal to the sedimentation rate. The model was calibrated and the parameters were obtained by comparing simulated and measured data. Very good agreement between all the three observed profiles and the simulated ones demonstrates the reliability of the model and the parameters. Simulation results show that the shapes of the curves are controlled by the combination of sedimentation and upper boundary conditions. Diffusion solely is adequate to reconstruct the observed profiles, which indicates that diffusion is the dominant vertical transport mechanism. The effective diffusion coefficients of the aquitard and the aquifer equal to 5.0 × 10^?11 and 2.0 × 10^?10 m²/s, respectively. The results of this study will help in understanding the transport mechanisms of solutes and environmental tracers in deltas with geology and hydrogeology similar to the PRD. Copyright © 2015 John Wiley & Sons, Ltd.     

Investigations of local Bouguer gravity anomalies of the Apennines and Taurus Mountains of the Moon

Local line-of-sight (LOS) Bouguer gravity anomalies of the Apennines and Taurus Mountains of the Moon have been calculated form low-altitude LOS free-air Doppler gravity profiles. The topography of the mountain areas is reflected by free-air gravity highs indicating no complete isostatic compensation. The resultant Bouguer gravity shows no anomalies for the Apennines, indicating lack of isostatic compensation. For the older Taurus Mountains significant local Bouguer minima of about ? 15 mgal indicate at least partial compensation.If a viscoelastic compensation mechanism (bending of a viscoelastic plate overlying a fluid half-space) is assumed, models for the crustal viscosity as a function of time give limits of the range of possible models from 10²⁴ to 5 × 10²⁵ P at 4.4 × 10⁹ y BP, 10²⁶ – 10²⁷ P at 3.9 × 10⁹ y, and 5 × 10²⁶ – 10²⁸ P at 3.0 × 10⁹ y. For earlier times only a lower bound of 10²⁷ P can be given.Two profiles of the Taurus area have been investigated; they show no significant Bouguer anomalies across the mare basalt patches of Lacus Bonitatis and Sinus Amoris and thus can be used to estimate an upper limit for the basalt thicknesses. For Lacus Bonitatis this limit is 1.3 km; the limit is reached for Sinus Amoris at an average thickness of 0.3 km, with 1.5 km in the centre. Earlier results from DeHon and Waskom are consistent with the gravity data.     

Groundwater potential of the Egbe-Mopa basement area,central Nigeria

Abstract

An investigation on the groundwater potentials of the Egbe-Mopa area in central Nigeria, underlain by the Basement Complex, is presented. The investigation involved mapping of the subsurface by use of vertical electrical soundings; measurement of depth to groundwater; and evaluation of hydraulic conductivity, transmissivity and yield by means of pumping test interpretation. The results indicate subsurface units that range from three to five resistivity layers; depth to groundwater of 0–10 m; overburden thickness of 3–16 m; hydraulic conductivity of 6.2?×?10^?6 to 3.4?×?10^?4 m/s; transmissivity of 4.3?×?10^?7 to 2?×?10^?3 m²/s; and groundwater yield of 0.2–2.5 L/s. The hydraulic head assessments revealed a general northward groundwater flow direction. The study identified three aquifer potential types, of high, medium and low productivity, respectively. Based on the longitudinal conductance of the overburden units, four distinct Aquifer Protective Capacity zones were delineated, namely, poor, weak, moderate and good.

Citation Okogbue, C.O. and Omonona, O.V., 2013. Groundwater potential of Egbe-Mopa basement area, central Nigeria. Hydrological Sciences Journal, 58 (4), 826–840.     

Water Quality Assessment and Resource Potentials: the Case of Aba–Urban and its Environs,Niger Delta Basin

The quest for improved water supply to cater for the ever increasing population has given rise to the assessment of water quality and resource potentials in Aba-Urban and its environments. The area, which lies within the Niger Delta Basin is underlain by the Benin Formation that is highly aquiferous. Samples of sands and water were subjected to sedimentological and hydrochemical analysis, respectively. Result reveals that the aquifer is thick (over 100 m) and unconfined. The computed aquifer parameters indicate high yielding clean-sands with hydraulic conductivity values ranging from 1.13 × 10^–4 to 5.70 × 10^–3 m/s. The specific discharge is about 14.2 m/year while the average linear groundwater velocity is calculated to be 53.46 m/year. Hydrochemical investigations carried out on water samples from Aba River and the groundwater system revealed low dissolved geochemical constituents. Although, there is a slight increase from the north to the southern part along the flow path. In general, the groundwater is relatively enriched in Ca, Mg?HCO^-₃ions and is predominantly of bicarbonate constituents. However, in some places Cl^? ions dominate over HCO^-₃ions. This is common where the water has come in contact with domestic sewage. This kind of contamination is also accompanied by elevated concentration of NO^-₃ions. Generally, the groundwater in most cases meets the standard for human consumption and is a better alternative to surface water. Consequently an increase exploitation of the prolific Benin aquifer through more hygienic and safe methods will surely be the best way of improving the domestic water supply situation in Aba City.     

On the use of the groundwater replenishment periods approach to determining aquifer parameters

Chalk crops out from a wide belt around the Paris basin, France, covering an area of about 70 000 km². In this region, the chalk presents the most important unconfined aquifer because of its extent and the size of its resources ((11–12) × 10⁹ m³ year^?1). The assessment of underground outflow depends on the vertical feeding, the infiltration and the hydrometry. This paper analyses the regional structural map, interprets the groundwater reaction under rainfall, explains the water circulation in such media where the reservoir geometry plays an essential role on aquifer response, and determines different aquifer physical parameters. Attempts are made to identify the rapid transfer of groundwater at the level of faults and their important fissures. Based on a quantitative study of the seasonal and interannual piezometry fluctuations, it is noted that the mode and the piezometric chronology events are controlled by major parameters of geological and hydrogeologic contexts, aquifer hydraulic characteristics, the position of upstream and downstream basin limits, groundwater depth and replenishment time. This paper ends with determination of groundwater physical parameters of the diffusivity (T/S, T and S) values by examination of the groundwater replenishments periods. Copyright © 2005 John Wiley & Sons, Ltd.     

Coupled groundwater flow and heat transport simulation for estimating transient aquifer–stream exchange at the lowland River Spree (Germany)

Subsurface flow and heat transport near Freienbrink, NE Germany, was simulated in order to study groundwater–surface water exchange between a floodplains aquifer and a section of the lowland River Spree and an adjacent oxbow. Groundwater exfiltration was the dominant process, and only fast surface water level rises resulted in temporary infiltration into the aquifer. The main groundwater flow paths are identified based on a 3D groundwater flow model. To estimate mass fluxes across the aquifer–surface water interfaces, a 2D flow and heat transport modelling approach along a transect of 12 piezometers was performed. Results of steady‐state and transient water level simulations show an overall high accuracy with a Spearman coefficient ρ = 0.9996 and root mean square error (RMSE) = 0.008 m. Based on small groundwater flow velocities of about 10^?7 to 10^?6 ms^?1, mean groundwater exfiltration rates of 233 l m^?2 d^?1 are calculated. Short periods of surface water infiltration into the aquifer do not exceed 10 days, and the infiltration rates are in the same range. The heat transport was modelled with slightly less accuracy (ρ = 0.8359 and RMSE = 0.34 °C). In contrast to the predominant groundwater exfiltration, surface water temperatures determine the calculated temperatures in the upper aquifer below both surface water bodies down to 10 m during the whole simulation period. These findings emphasize prevailing of heat conduction over advection in the upper aquifer zones, which seems to be typical for lowland streams with sandy aquifer materials and low hydraulic gradients. Moreover, this study shows the potential of coupled numerical flow and heat transport modelling to understand groundwater–surface water exchange processes in detail. Copyright © 2013 John Wiley & Sons, Ltd.