Role of water flow in modeling methane emissions from flooded paddy soils

Methane (CH₄) is a potent greenhouse gas that is emitted from paddy fields, and the large CH₄ fluxes represent a worldwide issue for the rice production eco-compatibility. In this work a model is proposed to investigate the role of water flows on CH₄ emissions from flooded paddy soils. The model is based on a system of partial differential mass balance equations of the chemical species affecting CH₄ fate, and water flows are modeled by the Darcy equation. Moreover, in order to properly model the dynamics of CH₄, a number of physico-chemical processes and features not included in currently available CH₄ emission models are considered: paddy soil stratigraphy; nutrient adsorption and root water uptake; gas transport and respiration within root aerenchyma compartment. The proposed model allows to simulate the spatio-temporal dynamics of chemical compounds within paddy soil as well as to quantify the influence of different processes on nutrient input/output budgets. Simulations without water flow have shown a considerable overestimation of CH₄ emissions due to a different spatio-temporal dynamics of dissolved organic matter (DOC – source of energy for CH₄ production). In particular, when water fluxes have not been modeled the overestimation can reach 54%, 41% and 67% of daily minimum, daily maximum, and total over the whole growing season CH₄ emission, respectively. Moreover, the model results suggest that roots influence CH₄ dynamics principally due to their nutrient uptake, while root effect on advective flow plays a minor role. Finally, the analysis of CH₄ transport fluxes has shown the limiting effect of upward dispersive transport fluxes on the downward CH₄ percolation.     

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.     

Changes in groundwater bacterial community during cyclic groundwater-table variations

Column experiments containing an aquifer sand were subjected to static and oscillating water tables to investigate the impact of natural fluctuations and rainfall infiltration on the groundwater bacterial community just below the phreatic surface, and its association with the geochemistry. Once the columns were established, the continuously saturated zone was anoxic in all three columns. The rate of soil organic matter (SOM) mineralization was higher when the water table varied cyclically than when it was static due to the greater availability of NO₃⁻ and SO₄²⁻. Natural fluctuations in the water table resulted in a similar NO₃⁻ concentration to that observed with a static water table but the cyclic wetting of the intermittently saturated zone resulted in a higher SO₄²⁻ concentration. Rainfall infiltration induced cyclic water-table variations resulted in a higher NO₃⁻ concentration than those in the other two columns, and a SO₄²⁻ concentration intermediate between those columns. As rainwater infiltration resulted in slow downward displacement of the groundwater, it is inferred that NO₃⁻ and SO₄²⁻ were being mobilized from the vadose zone. NO₃⁻ was mainly released by SOM mineralization (which was enhanced by the infiltration of oxygenated rainwater), but the larger amount of SO₄²⁻ release required a second mechanism (possibly desorption). Different groundwater bacterial communities evolved from initially similar populations due to the different groundwater histories.     

Rapid nutrient load reduction during infiltration of managed aquifer recharge in an agricultural groundwater basin: Pajaro Valley,California

Artificial recharge of groundwater is an increasingly important method for augmenting groundwater supply and can have a positive or negative influence on the quality of water resources. We instrumented a managed aquifer recharge (MAR) pond in central coastal California to assess how patterns of infiltration and recharge affect the load of nitrate delivered to the underlying aquifer. The concentration of nitrate in infiltrating water consistently decreased during passage through the first metre of subsurface soils. Enrichment of ¹⁸O and ¹⁵ N in the residual nitrate in infiltrating water proceeded in a ratio of 1:2, indicating that denitrification plays a significant role in the quantitative reduction of nutrients exported during infiltration through shallow soils. The extent and rate of nitrate removal was spatially and temporally variable across the bottom of the recharge pond, with 30% to 60% of the nitrate load being removed over the first 6 weeks of managed aquifer recharge operation. During the period of highest N loading to the system, when the average infiltration rate was > 1 m/day, the recharge pond achieved a load reduction efficiency of 7 kg NO₃^?‐N/day/ha, which compares favourably to nitrate load reductions achieved by treatment wetlands. Groundwater mounding and water composition below the recharge pond suggest that recharge and subsequent lateral transport occur heterogeneously in the underlying aquifer. Nitrate concentrations in the aquifer following infiltration were lowered primarily by dilution, with little evidence for additional denitrification occurring in the aquifer in comparison to high rates documented during shallow infiltration. Copyright © 2011 John Wiley & Sons, Ltd.     

城市景观湖泊不同时间尺度CH4通量特征及影响因素分析——以南京市莫愁湖为例

城市景观湖泊对温室气体的收支发挥着重要作用。本文以南京市莫愁湖为研究对象,采用静态箱—温室气体分析仪法实时监测湖泊水—气界面CH₄通量,分析湖泊主要温室气体CH₄在日尺度和季节尺度上因冒泡和扩散排放方式不同对其通量的影响,探究影响湖泊CH₄通量的因素。结果表明:(1)在日尺度上,四季24 h内CH₄均呈排放状态,受白天冒泡影响,四季CH₄总通量均存在白天高于夜间的日变化特征。(2)在季节尺度上,莫愁湖CH₄排放通量呈现显著的时空异质性,受冒泡通量的影响夏季CH₄通量明显高于春、秋、冬三季;B区的CH₄总通量(6.04 nmol/(m²·s))显著高于A区(3.82 nmol/(m²·s)),水体的营养化程度和离岸距离是空间变化的主要影响因素。A、B两区CH₄排放夏季以冒泡排放为主,春、秋、冬以扩散排放为主。(3)在日尺度上,CH₄     

Identification of groundwater recharge sources and processes in a heterogeneous alluvial aquifer: results from multi‐level monitoring of hydrochemistry and environmental isotopes in a riverside agricultural area in Korea

We evaluated sources and pathways of groundwater recharge for a heterogeneous alluvial aquifer beneath an agricultural field, based on multi‐level monitoring of hydrochemistry and environmental isotopes of a riverside groundwater system at Buyeo, Korea. Two distinct groundwater zones were identified with depth: (1) a shallow oxic groundwater zone, characterized by elevated concentrations of NO₃^? and (2) a deeper (>10–14 m from the ground surface) sub‐oxic groundwater zone with high concentrations of dissolved Fe, silica, and HCO₃^?, but little nitrate. The change of redox zones occurred at a depth where the aquifer sediments change from an upper sandy stratum to a silty stratum with mud caps. The δ¹⁸O and δ²H values of groundwater were also different between the two zones. Hydrochemical and δ¹⁸O? δ²H data of oxic groundwater are similar to those of soil water. This illustrates that recharge of oxic groundwater mainly occurs through direct infiltration of rain and irrigation water in the sandy soil area where vegetable cropping with abundant fertilizer use is predominant. Oxic groundwater is therefore severely contaminated by agrochemical pollutants such as nitrate. In contrast, deeper sub‐oxic groundwater contains only small amounts of dissolved oxygen (DO) and NO₃^?. The ³H contents and elevated silica concentrations in sub‐oxic groundwater indicate a somewhat longer mean residence time of groundwater within this part of the aquifer. Sub‐oxic groundwater was also characterized by higher δ¹⁸O and δ²H values and lower d‐excess values, indicating significant evaporation during recharge. We suggest that recharge of sub‐oxic groundwater occurs in the areas of paddy rice fields where standing irrigation and rain water are affected by strong evaporation, and that reducing conditions develop during subsequent sub‐surface infiltration. This study illustrates the existence of two groundwater bodies with different recharge processes within an alluvial aquifer. Copyright © 2009 John Wiley & Sons, Ltd.     

Dissolved Organic Carbon in Groundwater Overlain by Irrigated Sugarcane

Elevated dissolved organic carbon (DOC) has been detected in groundwater beneath irrigated sugarcane on the Burdekin coastal plain of tropical northeast Australia. The maximum value of 82 mg/L is to our knowledge the highest DOC reported for groundwater beneath irrigated cropping systems. More than half of the groundwater sampled in January 2004 (n = 46) exhibited DOC concentrations greater than 30 mg/L. DOC was progressively lower in October 2004 and January 2005, with a total decrease greater than 90% indicating varying load(s) to the aquifer. It was hypothesized that the elevated DOC found in this groundwater system is sourced at or near the soil surface and supplied to the aquifer via vertical recharge following above average rainfall. Possible sources of DOC include organic‐rich sugar mill by‐products applied as fertilizer and/or sugarcane sap released during harvest. CFC‐12 vertical flow rates supported the hypothesis that elevated DOC (>40 mg/L) in the groundwater results from recharge events in which annual precipitation exceeds 1500 mm/year (average = 960 mm/year). Occurrence of elevated DOC concentrations, absence of electron acceptors (O₂ and NO₃^–) and both Fe²⁺ and Mn²⁺ greater than 1 mg/L in shallow groundwater suggest that the DOC compounds are chemically labile. The consequence of high concentrations of labile DOC may be positive (e.g., denitrification) or negative (e.g., enhanced metal mobility and biofouling), and highlights the need to account for a wider range of water quality parameters when considering the impacts of land use on the ecology of receiving waters and/or suitability of groundwater for irrigated agriculture.     

Hydrological control of dissolved organic carbon dynamics in a forested headwater catchment,Kiryu Experimental Watershed,Japan

To evaluate the influence of hydrological processes on dissolved organic carbon (DOC) dynamics in a forested headwater catchment, DOC concentration was observed along the flow path from rainfall to stream water via throughfall, soil water, groundwater, and spring water for 4 years, and DOC flux through the catchment was calculated. The spatial and temporal variations in DOC concentration and flux were compared with physical hydrological observations and the mean residence time of water. In the upslope soil layer, DOC concentrations were not significantly correlated with water fluxes, suggesting that DOC concentrations were not strictly controlled by water fluxes. In the upslope perennial groundwater, DOC concentration was affected by the change in the amount of microbial degradation of DOC produced by changes in the mean residence time of water. In stream water, the temporal variation in DOC concentration was usually affected by changes in DOC concentration of the inflow component via vertical infiltration from above the perennial groundwater. During dry periods, however, the component from inflow via vertical infiltration was negligible and DOC in the upslope perennial groundwater became the major component of stream water DOC. The temporal variation in stream water DOC concentration during baseflow was affected by rainfall patterns over several preceding months. Therefore, records of rainfall over several preceding months are one of the most important factors for predicting changes in DOC concentration on a catchment scale. Copyright © 2007 John Wiley & Sons, Ltd.     

Modelling water infiltration in cracked paddy field soil

Periodic paddy field flooding is a major source of groundwater recharge. Many paddy fields thus are used as groundwater recharge ponds after harvesting the first crop of the summer. Following rice harvesting, paddy field surfaces may crack into fissures as a result of drainage and exposure to sunlight. Field observation indicates that applying precipitation to the paddy field can increase the rate of infiltration. To quantitatively evaluate the amount of infiltration in a cracked paddy field, this study sets up a simple soil crack model to simulate the field infiltration process. A three‐dimensional groundwater model FEMWATER is adopted to simulate water movement in the paddy field subjected to various crack conditions. Using the field and laboratory data of irrigation water requirements, soil physical properties, hydraulic conductivities and soil profiles obtained from Ten‐Chung, FEMWATER simulates the water movement in the dry cracked paddy. Simulation results show that if the cracks develop extensively and penetrate the ploughed soil, the infiltration rate may increase significantly. The infiltration fluxes of crack with depths of 80, 60 and 27·5 cm are 18·77, 14·50 and 8·06 times higher than that of 20 cm, respectively. The simulation results of cracks with 80 cm depth correlated closely with field observations. The results of the study elucidate the processes of unsaturated water movement in a dry cracked paddy field. Copyright © 2004 John Wiley & Sons, Ltd.     

Estimating groundwater recharge in lowland watersheds

Little attention has been given to the role of groundwater in the hydrological cycle of lowland watersheds. Our objective in this study was to estimate total recharge to groundwater by analysing water table response to storm events and the rate at which water was transferred into the shallow aquifer. This was conducted at three sites in a rural watershed in the lower Atlantic coastal plain near Charleston, South Carolina, USA. A novel version of the water table fluctuation method was used to estimate total recharge to the shallow aquifer by comparing hourly data of water table position following storm events and measuring water table recession behavior, rather than subjective graphical analysis methods. Also, shallow aquifer recharge rates (vertical fluxes) were estimated using Darcy's Law by comparing static water levels in a water table well and in a shallow piezometer during dry periods. The total annual recharge estimated ranged from 107 ± 39 mm·yr^–1 (5–10% of annual precipitation) at a poorly drained topographic low area to 1140 ± 230 mm·yr^–1 (62–94% of annual precipitation) for a moderately well‐drained upland site. The average aquifer recharge rate was 114 ± 60 mm·yr^–1, which is similar to previous estimations of base flow for the ephemeral third‐order streams in this watershed. The difference in the two methods may have been caused by processes not accounted for in the Darcy flux method, soil moisture deficits, and average evapotranspiration demand, which is about 1100 mm·yr^–1 for this region. Although other factors also can affect partitioning of recharge, an integrated approach to inspecting easily gathered groundwater data can provide information on an often neglected aspect of water budget estimation. We also discuss the effects of land use change on recharge reduction, given a typical development scenario for the region. Copyright © 2011 John Wiley & Sons, Ltd.     

Disconnected Surface Water and Groundwater: From Theory to Practice

When describing the hydraulic relationship between rivers and aquifers, the term disconnected is frequently misunderstood or used in an incorrect way. The problem is compounded by the fact that there is no definitive literature on the topic of disconnected surface water and groundwater. We aim at closing this gap and begin the discussion with a short introduction to the historical background of the terminology. Even though a conceptual illustration of a disconnected system was published by Meinzer (1923) , it is only within the last few years that the underlying physics of the disconnection process has been described. The importance of disconnected systems, however, is not widely appreciated. Although rarely explicitly stated, many approaches for predicting the impacts of groundwater development on surface water resources assume full connection. Furthermore, management policies often suggest that surface water and groundwater should only be managed jointly if they are connected. However, although lowering the water table beneath a disconnected section of a river will not change the infiltration rate at that point, it can increase the length of stream that is disconnected. Because knowing the state of connection is of fundamental importance for sustainable water management, robust field methods that allow the identification of the state of connection are required. Currently, disconnection is identified by showing that the infiltration rate from a stream to an underlying aquifer is independent of the water table position or by identifying an unsaturated zone under the stream. More field studies are required to develop better methods for the identification of disconnection and to quantify the implications of heterogeneity and clogging processes in the streambed on disconnection.     

Impact of Hurricane Irene and Tropical Storm Lee on riparian zone hydrology and biogeochemistry

Although riparian zones are well known to reduce nitrogen (N) and phosphorus (P) runoff to streams, they also have the potential to affect greenhouse gas (CO₂, N₂O, and CH₄) fluxes to the atmosphere. Following large storms, soil biogeochemical conditions often become more reduced, especially in oxbow depressions and side channels, which can lead to hot moments of greenhouse gas production. Here, we investigate the impact of the remnants of Hurricane Irene and Tropical Storm Lee on riparian zone hydrology (water table: WT), and biogeochemistry (oxidation‐reduction potential [ORP], dissolved oxygen [DO], NO₃^?, PO₄^3?, CO₂, N₂O, CH₄). Results indicate that large storms have the potential to reset WT levels for weeks to months. Overbank flooding at our site following Irene and Lee led to the infiltration of well‐oxygenated water at depth (higher DO and ORP) while promoting the development of anoxic conditions within soil aggregates near the soil surface (increased N₂O and CH₄ fluxes). A short‐term increase in CO₂ emission was observed following Irene at our site where aerobic respiration was water‐limited. Over a 2‐year period, an oxbow depression exhibited higher WT, higher N₂O and CH₄ fluxes (hot moment), higher CO₂ fluxes (seasonal), and lower NO₃^? concentrations (seasonal) than the rest of the riparian zone. However, neither Irene, nor Lee, nor the oxbow depression significantly impacted PO₄^3?. Dissolved organic carbon, ORP, and DO data illustrate the time‐lag (>20 years) between the creation of an oxbow depression and the development of reducing conditions despite clear differences in riparian zone and oxbow WT dynamics.     

Field assessment of surface water–groundwater connectivity in a semi‐arid river basin (Murray–Darling,Australia)

In semi‐arid and arid river basins, understanding the connectivity between rivers and alluvial aquifers is one of the key challenges for the management of groundwater resources. The type of connection present (gaining, losing‐connected, transitional and losing‐disconnected) was assessed at 12 sites along six Murray–Darling Basin river reaches. The assessments were made by measuring the hydraulic head in the riparian zone near the rivers to evaluate if the water tables intersected the riverbeds and by measuring fluid pressure (ψ) in the riverbeds. The rationale for the latter was that ψ will always be greater than or equal to zero under connected conditions (either losing or gaining) and always lesser than or equal to zero under losing‐disconnected conditions. A mixture of losing‐disconnected, losing‐connected and gaining conditions was found among the 12 sites. The losing‐disconnected sites all had a riverbed with a lower hydraulic conductivity than the underlying aquifer, usually in the form of a silty clay or clay unit 0.5–2 m in thickness. The riparian water tables were 6 to 25 m below riverbed level at the losing‐disconnected sites but never lower than 1 m below riverbed level at the losing‐connected ones. The contrast in water table depth between connected and disconnected sites was attributed to the conditions at the time of the study, when a severe regional drought had generated a widespread decline in regional water tables. This decline was apparently compensated near losing‐connected rivers by increased infiltration rates, while the decline could not be compensated at the losing‐disconnected rivers because the infiltration rates were already maximal there. Copyright © 2012 John Wiley & Sons, Ltd.     

Benzene Dynamics and Biodegradation in Alluvial Aquifers Affected by River Fluctuations

The spatial distribution and temporal dynamics of a benzene plume in an alluvial aquifer strongly affected by river fluctuations was studied. Benzene concentrations, aquifer geochemistry datasets, past river morphology, and benzene degradation rates estimated in situ using stable carbon isotope enrichment were analyzed in concert with aquifer heterogeneity and river fluctuations. Geochemistry data demonstrated that benzene biodegradation was on‐going under sulfate reducing conditions. Long‐term monitoring of hydraulic heads and characterization of the alluvial aquifer formed the basis of a detailed modeled image of aquifer heterogeneity. Hydraulic conductivity was found to strongly correlate with benzene degradation, indicating that low hydraulic conductivity areas are capable of sustaining benzene anaerobic biodegradation provided the electron acceptor (SO₄^2–) does not become rate limiting. Modeling results demonstrated that the groundwater flux direction is reversed on annual basis when the river level rises up to 2 m, thereby forcing the infiltration of oxygenated surface water into the aquifer. The mobilization state of metal trace elements such as Zn, Cd, and As present in the aquifer predominantly depended on the strong potential gradient within the plume. However, infiltration of oxygenated water was found to trigger a change from strongly reducing to oxic conditions near the river, causing mobilization of previously immobile metal species and vice versa. MNA appears to be an appropriate remediation strategy in this type of dynamic environment provided that aquifer characterization and targeted monitoring of redox conditions are adequate and electron acceptors remain available until concentrations of toxic compounds reduce to acceptable levels.     

Towards understanding organic matter fluxes and reactivity in surface waters: Filtering impact on DOC and POC degradation

Peatlands cover a very small area of the Earth, but store globally significant quantities of carbon and export disproportionate quantities of fluvial organic carbon, especially when the peatlands are degraded or disturbed. Peatland headwater catchments with high concentrations of dissolved and particulate organic carbon (DOC and POC) provide an opportunity to investigate the possibility of competing effects that could lead to enhanced or diminished turnover of DOC in the presence of POC. Both POC and DOC can be degraded by light and microbes, producing smaller molecules and releasing CO₂ and CH₄ to the atmosphere, and POC can inhibit light penetration, stabilize DOC by providing adsorption sites and providing surfaces for microbes to interact with DOC. However, the majority of peatland fluvial carbon studies are conducted using filtered water samples, and measure only the DOC concentration, so the impact of the particulate organic matter (POM) on in-stream processing of organic carbon is relatively unknown. It is therefore possible that studies have underestimated carbon transformations in rivers as they have not considered the interaction of the particulate material on the dissolved concentrations; there could be higher losses than previously estimated, increasing the contribution of peatland headwaters to GHG emissions. In this study, we assessed if the current approach of DOC degradation studies accurately represent the impact of POM on DOC degradation, by quantifying DOC production from POM, and therefore POC, over time in water with manipulated POM concentrations. Both filtered and unfiltered water lost 60% of the DOC over 70 hours, whereas the treatment with additional POM lost only 35%. The results showed that filtering does not significantly impact the DOC degradation rates; however, when the POC concentration was doubled, there was a significant reduction in DOC degradation, suggesting that filtering would still be necessary to get accurate rates of DOC transformations in waters with high POC concentrations.     

Modeling Surface Water-Groundwater Interaction with MODFLOW: Some Considerations

The accuracy with which MODFLOW simulates surface water-groundwater interaction is examined for connected and disconnected losing streams. We compare the effect of different vertical and horizontal discretization within MODFLOW and also compare MODFLOW simulations with those produced by HydroGeoSphere. HydroGeoSphere is able to simulate both saturated and unsaturated flow, as well as surface water, groundwater and the full coupling between them in a physical way, and so is used as a reference code to quantify the influence of some of the simplifying assumptions of MODFLOW. In particular, we show that (1) the inability to simulate negative pressures beneath disconnected streams in MODFLOW results in an underestimation of the infiltration flux; (2) a river in MODFLOW is either fully connected or fully disconnected, while in reality transitional stages between the two flow regimes exist; (3) limitations in the horizontal discretization of the river can cause a mismatch between river width and cell width, resulting in an error in the water table position under the river; and (4) because coarse vertical discretization of the aquifer is often used to avoid the drying out of cells, this may result in an error in simulating the height of the groundwater mound. Conditions under which these errors are significant are investigated.     

Quantifying River‐Groundwater Interactions of New Zealand's Gravel‐Bed Rivers: The Wairau Plain

New Zealand's gravel‐bed rivers have deposited coarse, highly conductive gravel aquifers that are predominantly fed by river water. Managing their groundwater resources is challenging because the recharge mechanisms in these rivers are poorly understood and recharge rates are difficult to predict, particularly under a more variable future climate. To understand the river‐groundwater exchange processes in gravel‐bed rivers, we investigate the Wairau Plain Aquifer using a three‐dimensional groundwater flow model which was calibrated using targeted field observations, “soft” information from experts of the local water authority, parameter regularization techniques, and the model‐independent parameter estimation software PEST. The uncertainty of simulated river‐aquifer exchange flows, groundwater heads, spring flows, and mean transit times were evaluated using Null‐space Monte‐Carlo methods. Our analysis suggests that the river is hydraulically perched (losing) above the regional water table in its upper reaches and is gaining downstream where marine sediments overlay unconfined gravels. River recharge rates are on average 7.3 m³/s, but are highly dynamic in time and variable in space. Although the river discharge regularly hits 1000 m³/s, the net exchange flow rarely exceeds 12 m³/s and seems to be limited by the physical constraints of unit‐gradient flux under disconnected rivers. An important finding for the management of the aquifer is that changes in aquifer storage are mainly affected by the frequency and duration of low‐flow periods in the river. We hypothesize that the new insights into the river‐groundwater exchange mechanisms of the presented case study are transferable to other rivers with similar characteristics.     

Characterising groundwater–surface water interactions in idealised ephemeral stream systems

Transmission losses from the beds of ephemeral streams are thought to be a widespread mechanism of groundwater recharge in arid and semi-arid regions and support a range of dryland hydro-ecology. Dryland areas cover ~40% of the Earth's land surface and groundwater resources are often the main source of freshwater. It is commonly assumed that where an unsaturated zone exists beneath a stream, the interaction between surface water and groundwater is unidirectional and that groundwater does not exert a significant feedback on transmission losses. To test this assumption, we conducted a series of numerical model experiments using idealised two-dimensional channel-transects to assess the sensitivity and degree of interaction between surface and groundwater for typical dryland ephemeral stream geometries, hydraulic properties and flow regimes. We broaden the use of the term ‘stream–aquifer interactions’ to refer not just to fluxes and water exchange but also to include the ways in which the stream and aquifer have a hydraulic effect on one another. Our results indicate that deep water tables, less frequent streamflow events and/or highly permeable sediments tend to result in limited bi-directional hydraulic interaction between the stream and the underlying groundwater which, in turn, results in high amounts of infiltration. With shallower initial depth to the water table, higher streamflow frequency and/or lower bed permeability, greater ‘negative’ hydraulic feedback from the groundwater occurs which in turn results in lower amounts of infiltration. Streambed losses eventually reach a constant rate as initial water table depths increase, but only at depths of 10s of metres in some of the cases studied. Our results highlight that bi-directional stream–aquifer hydraulic interactions in ephemeral streams may be more widespread than is commonly assumed. We conclude that groundwater and surface water should be considered as connected systems for water resource management unless there is clear evidence to the contrary.     

Bio-geological processes of nitrogen transport and transformation in the aeration zone and aquifer

Abstract

Nitrate contamination in groundwater originates mainly from excessive use of fertilizers and uncontrolled discharge to land of incompletely-treated wastewater associated with agricultural activities. A systematic field investigation was carried out in a sub-catchment of Dianchi Lake, Kunming, Yunnan, China, into the hydrological, biological and geological processes of nitrogen transport and transformation in the aeration zone and aquifer system. In situ experiments showed that the quantity of NO₃-N recharged into groundwater was related to fertilization. Nitrification and denitrification behaved quite differently but were affected by moisture content and Eh value. The vertical infiltration rate was controlled by the groundwater table and hydraulic conductivity of the soil. The existence of a zero-flux plane reflected the dynamics of water fluxes in the soil profile and Eh was measured in the aeration zone. In response to these factors, the nitrification rate was greatest in the top soil and reduced with the depth of soil; it was 6.53 mg/(kg·h) in the vegetated plot and 0.2–0.3 mg/(kg·h) in the unvegetated one. The denitrification rate in the unvegetated plot was 6.36 mg/(kg·h), and it was 2.79 mg/(kg·h) in the vegetated one.     

Hydrogeochemical characteristics of surface water and groundwater in the karst basin,southwest China

Groundwater is a very significant water source used for irrigation and drinking purposes in the karst region, and therefore understanding the hydrogeochemistry of karst water is extremely important. Surface water and groundwater were collected, and major chemical compositions and environmental isotopes in the water were measured in order to reveal the geochemical processes affecting water quality in the Gaoping karst basin, southwest China. Dominated by Ca²⁺, Mg²⁺, HCO₃^? and SO₄^2?, the groundwater is typically characterized by Ca? Mg? HCO₃ type in a shallow aquifer, and Ca? Mg? SO₄ type in a deeper aquifer. Dissolution of dolomite aquifer with gypsiferous rocks and dedolomitization in karst aquifers are important processes for chemical compositions of water in the study basin, and produce water with increased Mg²⁺, Ca²⁺ and SO₄^2? concentrations, and also increased TDS in surface water and groundwater. Mg²⁺/Ca²⁺ molar ratios in groundwater decrease slightly due to dedolomitization, while the mixing of discharge of groundwater with high Mg²⁺/Ca²⁺ ratios may be responsible for Mg²⁺/Ca²⁺ ratios obviously increasing in surface water, and Mg²⁺/Ca²⁺ ratios in both surface water and groundwater finally tending to a constant. In combination with environmental isotopic analyses, the major mechanism responsible for the water chemistry and its geochemical evolution in the study basin can be revealed as being mainly from the water–rock interaction in karst aquifers, the agricultural irrigation and its infiltration, the mixing of surface water and groundwater and the water movement along faults and joints in the karst basin. Copyright © 2009 John Wiley & Sons, Ltd.