Numerical modelling of stream–aquifer interaction: Quantifying the impact of transient streambed permeability and aquifer heterogeneity

Stream–aquifer interaction plays a vital role in the water cycle, and a proper study of this interaction is needed for understanding groundwater recharge, contaminants migration, and for managing surface water and groundwater resources. A model‐based investigation of a field experiment in a riparian zone of the Schwarzbach river, a tributary of the Rhine River in Germany, was conducted to understand stream–aquifer interaction under alternative gaining and losing streamflow conditions. An equivalent streambed permeability, estimated by inverting aquifer responses to flood waves, shows that streambed permeability increased during infiltration of stream water to aquifer and decreased during exfiltration. Aquifer permeability realizations generated by multiple‐point geostatistics exhibit a high degree of heterogeneity and anisotropy. A coupled surface water groundwater flow model was developed incorporating the time‐varying streambed permeability and heterogeneous aquifer permeability realizations. The model was able to reproduce varying pressure heads at two observation wells near the stream over a period of 55 days. A Monte Carlo analysis was also carried out to simulate groundwater flow, its age distribution, and the release of a hypothetical wastewater plume into the aquifer from the stream. Results of this uncertainty analysis suggest (a) stream–aquifer exchange flux during the infiltration periods was constrained by aquifer permeability; (b) during exfiltration, this flux was constrained by the reduced streambed permeability; (c) the effect of temporally variable streambed permeability and aquifer heterogeneity were found important to improve the accurate capture of the uncertainty; and (d) probabilistic infiltration paths in the aquifer reveal that such pathways and the associated prediction of the extent of the contaminant plume are highly dependent on aquifer heterogeneity.     

Scales and magnitude of hyporheic,river–aquifer and bank storage exchange fluxes

Many studies have investigated the exchange processes that occur between rivers and groundwater systems and have successfully quantified the water fluxes involved. Specifically, these exchange processes include hyporheic exchange, river–aquifer exchange (groundwater discharge and river loss) and bank storage exchange. Remarkably, there are relatively few examples of field studies where more than one exchange process is quantified, and as a consequence, the relationships between them are not well understood. To compare the relative magnitudes of these common exchange processes, we have collected data from 54 studies that have quantified one or more of these exchange flux types. Each flux value is plotted against river discharge at the time of measurement to allow the different exchange flux types to be compared. We show that there are positive relationships between the magnitude of each exchange flux type and increasing river discharge across the different studies. For every one order of magnitude increase in river discharge, the hyporheic, river–aquifer and bank storage exchange fluxes increase by factors of 2.7, 2.9 and 2.5, respectively. On average, hyporheic exchange fluxes are almost an order of magnitude greater than river–aquifer exchange fluxes, which are, in turn, approximately four times greater than bank storage exchange fluxes for the same river discharge. Unless measurement approaches that can distinguish between different types of exchange flux are used, there is potential for hyporheic exchange fluxes to be misinterpreted as river–aquifer exchange fluxes, with possible implications for water resource management decisions. Copyright © 2015 John Wiley & Sons, Ltd.     

Variations in dissolved CO2 and CH4 in a first‐order stream and catchment: an investigation of soil–stream linkages

Spatial and seasonal variations in CO₂ and CH₄ concentrations in streamwater and adjacent soils were studied at three sites on Brocky Burn, a headwater stream draining a peatland catchment in upland Britain. Concentrations of both gases in the soil atmosphere were significantly higher in peat and riparian soils than in mineral soils. Peat and riparian soil CO₂ concentrations varied seasonally, showing a positive correlation with air and soil temperature. Streamwater CO₂ concentrations at the upper sampling site, which mostly drained deep peats, varied from 2·8 to 9·8 mg l^?1 (2·5 to 11·9 times atmospheric saturation) and decreased markedly downstream. Temperature‐related seasonal variations in peat and riparian soil CO₂ were reflected in the stream at the upper site, where 77% of biweekly variation was explained by an autoregressive model based on: (i) a negative log‐linear relationship with stream flow; (ii) a positive linear relationship with soil CO₂ concentrations in the shallow riparian wells; and (iii) a negative linear relationship with soil CO₂ concentrations in the shallow peat wells, with a significant 2‐week lag term. These relationships changed markedly downstream, with an apparent decrease in the soil–stream linkage and a switch to a positive relationship between stream flow and stream CO₂. Streamwater CH₄ concentrations also declined sharply downstream, but were much lower (<0·01 to 0·12 mg l^?1) than those of CO₂ and showed no seasonal variation, nor any relationship with soil atmospheric CH₄ concentrations. However, stream CH₄ was significantly correlated with stream flow at the upper site, which explained 57% of biweekly variations in dissolved concentrations. We conclude that stream CO₂ can be a useful integrative measure of whole catchment respiration, but only at sites where the soil–stream linkage is strong. Copyright © 2004 John Wiley & Sons, Ltd.     

Sedimentary and geochemical characterization of Middle–Late Miocene formations in the Neogene Tsugaru Basin,Japan by means of DTH27‐1 well sediment analysis

Middle–Late Miocene age siliceous formations outcropping along the northwestern side of Honshu Island are considered prospective source rocks for hydrocarbons. An analysis of geophysical, sedimentological, and geochemical properties is essential to evaluate the formations' source potential, and to understand the factors that determined the accumulation and preservation of organic matter. This study investigates the Middle–Late Miocene geological record of the Tsugaru back‐arc basin, located in the western part of Aomori prefecture, through an analysis of a 200 m long portion of a core from the DTH27‐1 well; this core is composed of the diatomaceous siltstones of the Akaishi Formation and the siliceous mudstones of the Odoji Formation. Sedimentological and geophysical characterization showed that the Akaishi Formation's diatomaceous siltstones are mostly massive and bioturbated, have low magnetic susceptibility, and demonstrate moderate natural radioactivity. Although the Odoji Formation's siliceous mudstones are massive, they have exceedingly low magnetic susceptibility and high natural radioactivity. Geochemical data from a Rock‐Eval Pyrolysis such as total organic carbon and generative potential (S1 + S2) revealed that, in the Tsugaru area, only the Odoji Formation is a likely prospective source rock for hydrocarbons. On the other hand, T_max values indicate that both the formations are thermally immature for generating hydrocarbons. The difference between the Akaishi and Odoji Formation in the sedimentological facies, in terms of the degree of bioturbation and the organic carbon content, indicates variations in lithological properties, such as porosity and grain size; moreover, this difference indicates a variation in the paleo‐oxygenation of bottom waters, with the transition from oxygen‐deficient conditions in the Middle Miocene to the more oxygenated conditions in the Late Miocene. Both the lithological and paleo‐environmental factors possibly influenced the organic richness in the two formations.     

Timing and development of sedimentation of the Cleaverville Formation and a post‐accretion pull‐apart system in the Cleaverville area,coastal Pilbara Terrane,Pilbara, Western Australia

In the Cleaverville area of Western Australia, the Regal, Dixon Island, and Cleaverville Formations preserve a Mesoarchean lower‐greenschist‐facies volcano‐sedimentary succession in the coastal Pilbara Terrane. These formations are distributed in a rhomboidal‐shaped area and are unconformably overlain by two narrowly distributed shallow‐marine sedimentary sequences: the Sixty‐Six Hill and Forty‐Four Hill Members of the Lizard Hills Formation. The former member is preserved within the core of the Cleaverville Syncline and the latter formed along the northeast‐trending Eighty‐Seven Fault. Based on the metamorphic grade and structures, two deformation events are recognized: D₁ resulted in folding caused by a collisional event, and D₂ resulted in regional sinistral strike‐slip deformation. A previous study reported that the Cleaverville Formation was deposited at 3020 Ma, after the Prinsep Orogeny (3070–3050 Ma). Our SHRIMP U–Pb zircon ages show that: (i) graded volcaniclastic–felsic tuff within the black shale sequence below the banded iron formation in the Cleaverville Formation yields an age of (3 114 ±14) Ma; (ii) the youngest zircons in sandstones of the Sixty‐Six Hill Member, which unconformably overlies pillow basalt of the Regal Formation, yield ages of 3090–3060 Ma; and (iii) zircons in sandstones of the Forty‐Four Hill Member show two age peaks at 3270 Ma and 3020 Ma. In this way, the Cleaverville Formation was deposited at 3114–3060 Ma and was deformed at 3070–3050 Ma (D₁). Depositional age of the Cleaverville Formation is at least 40–90 Myr older than that proposed in previous studies and pre‐dates the Prinsep Orogeny (3070–3050 Ma). After 3020 Ma, D₂ resulted in the formation of a regional strike‐slip pull‐apart basin in the Cleaverville area. The lower‐greenschist‐facies volcano‐sedimentary rocks are distributed only within this basin structure. This strike‐slip deformation was synchronous with crustal‐scale sinistral shear deformation (3000–2930 Ma) in the Pilbara region.     

Boundary crossing and non‐linear theory in earth‐system sciences – a proof of concept based on tsunami and post‐eruption scenarios on Java Island,Indonesia

Although volcanic eruptions are well‐known to be the trigger of some weather and climatic changes, land‐cover changes by pyroclastic‐flows and lahars do not get this recognition, neither do major hazards such as tsunami. These two earth processes are even lesser considered as being able to modify other earth processes they are not directly connected to, such as landslides or river discharge in non‐connected basins more than a hundred kilometres away. In this contribution the authors argue that these ideas are mainly driven by the process of being ‘educated’ in a single academic discipline and once put to the test interactions and retroactions between earth processes and atmospheric processes are far more reaching than commonly thought. For this study, the site of Java Island (Indonesia) was chosen to conduct (1) an analysis of a major tsunami impact – in the same area as the 2006 Java tsunami and (2) an analysis of the post‐eruption impacts of Merapi Volcano after a major eruption – excluding any ejecta in the atmosphere for the sake of the demonstration. The atmospheric feedback simulations were conducted using the regional climate model (RegCM‐4) with calibration from weather stations in Java Island. As a result, both simulations have proven that large scale deposits of pyroclasts (not introducing the ejectas sent in the atmosphere) and tsunamis can have outstanding impacts on the atmospheric situation and the bio‐geomorphologic evolution of the landscape in the following weeks to months. Interestingly enough these impacts are not limited to the area impacted by the earth process and the effect are not linear in time as they work following thresholds. These rainfalls ‘tele‐impacts’ are important enough to, in turn, modify earth‐surface processes in areas remote from the original phenomenon. This system acts in the same manner as a famous butterfly in Africa that could trigger a hurricane on the other side of the Atlantic Ocean. Copyright © 2012 John Wiley & Sons, Ltd.     

Water and air dynamics within a deep vadose zone of a karst massif: Observations from the Lukina jama–Trojama cave system (−1,431 m) in Dinaric karst (Croatia)

Extreme heterogeneity of karst systems makes them very challenging to study. Various processes within the system affect its global response, usually measured at karst springs. Research conducted in caves provides a unique opportunity for in situ analysis of separate processes in karst underground. The aim of the present study was to research the water and air dynamics within a deep karst system. Air and water basic physical parameters across the Lukina jama–Trojama cave system (?1,431 m) were continuously monitored during a 1‐year period. Recorded hydrograph of the siphon lake at the bottom of the cave was used to interpret the characteristics of an unexplored phreatic/epiphreatic conduit network. Water origin in the siphon was determined based on temperature and electrical conductivity. Air temperature and humidity monitoring revealed a strong inflow of air of sub‐zero temperature into the upper portion of the cave during winter. Cave passage morphology was interpreted as the main determinant of air dynamics, which caused ice to accumulate extensively in the upper portions of the cave and caused the temperature on the top of the homothermic zone to be significantly below the mean outside temperature. Air dynamics also lowered the temperature of water flowing through the cave vadose zone and feeding the phreatic zone of the massif. The pronounced temperature difference between the phreatic zone and the top of the homothermic zone probably contributed to the thermal gradient observed in the cave, which is steeper than in ice‐free caves in the area. Our results enabled the development of a conceptual model that describes coupling between air and water dynamics in the cave system and its surroundings.