Episodes 2019; 42(4): 333-341
Published online December 1, 2019
Copyright © International Union of Geological Sciences.
Haoyu Zhu1,3,4, Dongsheng Liu1,2,3*, Wenqing Shi1,3, and Qiuwen Chen1,3
1State Key Laboratory of Hydrology-Water Resources & Hydraulic Engineering, Nanjing Hydraulic Research Institute, Nanjing 210029, China
2College of Water Conservancy and Hydropower Engineering, Hohai University, Nanjing 210098, China
3Center for Eco-Environmental Research, Nanjing Hydraulic Research Institute, Nanjing 210098, China; *Corresponding author, E-mail: email@example.com
4School of River and Ocean Engineering, Chongqing Jiaotong University, Chongqing 400074, China
Correspondence to:*E-mail: firstname.lastname@example.org
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
By using a 2D-coupled flow and solute transport and reaction model across the stream-to-riparian continuum, this paper systematically studied the nitrogen-cycling processes driven by a flood wave and their spatial-temporal distributions. The influences on hyporheic nitrogen removal of different waves that vary by amplitude (A), duration (T), wave-type parameter (r) and rising duration (tp) were investigated. During the surface-water-infiltration period, the aerobic respiration, nitrification and denitrification in the hyporheic zone were gradually enhanced, and aerobic respiration was dominant. During the groundwater-backflow period, the reactions gradually weakened, and the role of denitrification was improved. The consumption of each solute continuously increased in the whole process, and the groundwater-backflow period corresponded to a larger consumption but at a smaller consumption rate. The reaction rates in space from high to low were as follows: riparian phreatic zone > variable saturated zone > stream bed. Hyporheic nitrogen removal increased with increasing A, T and tp and with decreasing r, which was contrary to the change rule of nitrogen removal efficiency. The hyporheic nitrogen removal exhibited a good linear relationship with the polynomial ‘A*T*tp/r’, which had great predictive significance.
The rapid fluctuations and lateral spreads of flood events greatly change the hydrological regime within the stream-floodplain system (Sawyer et al., 2009; Gerecht et al., 2011; Liu et al., 2018) and affect the hydrochemical indexes in the hyporheic zone (HZ) (Curry et al., 1994; Arntzen et al., 2006; Hanrahan, 2008; Jeon and Lee, 2015). The effect of a flood event on hydrochemical indexes is through the following two ways. (1) The solutes are carried into the HZ during the surface water infiltration, or the HZ behaves as a source of solutes during the groundwater discharging. (2) The solutes in the HZ undergo a series of chemical reactions that change the corresponding hydrochemical indexes.
Current studies on the hydrochemical reactions in the HZ under flood events focus on the nitrogen cycles. In the processes of hyporheic exchange and nitrogen cycles, surface water is usually rich in dissolved oxygen (i.e., O2), dissolved organic carbon (DOC) and nutrients, whereas these contents in groundwater are generally lower (Dahm et al., 1998). Therefore, when surface water level rises and recharges groundwater, it carries much of O2 into the HZ. In aerobic environment, NH4+ is converted into NO3- by the action of nitrobacteria (i.e., nitrification), which contributes to the NO3- source of the aquifer. Simultaneously, O2 and DOC also undergo aerobic respiration, thereby promoting the degradation of organic matter. When surface water level falls, heterotrophic microorganisms in the HZ preferentially use NO3- as an electron acceptor in the absence of oxygen to promote their metabolism and convert NO3- to N2 and N2O (i.e., denitrification). According to Böhlke and Denver (1995) and Marzadri et al. (2011), the threshold concentration of O2 for denitrification in stream sediment is 3 mg/l, assuming NO3- is completely reduced to N2.
Scholars have conducted a lot of researches on the hyporheic nitrogen cycles under flood events through field tracer tests. For example, Stanley and Boulton (1995) studied the nitrogen cycles in streamfloodplain system under alternating wet and dry processes, indicating that the concentrations of nitrates and dissolved oxygen drop sharply during flood attenuation. Moreover, with weakening hyporheic exchange, the HZ reaches an extremely hypoxic state (O2 concentration <0.5 mg/L), and the nitrate nitrogen concentration is also at a lower level (<0.02 mg/L), after which denitrification dominates. Roley et al. (2012) conducted an in situ measurement and analysis of the nitrate concentration and nitrification rate in the HZ under flood events in Shatto Creek, USA. The authors’ results showed that nitrification rate in the floodplain area gradually increases with the expansion of flood-submerged area and has different performances in the process of flood rise-peak-fall. Harvey et al. (2013) investigated denitrification in the HZ by injecting 15NO3- in the upstream of Sugar Creek, USA, and the HZ plays a dominant role in the denitrification process of the river system. Moreover, the denitrification rate is proportional to the magnitude of hyporheic exchange.
At present, researches on the mechanism of hyporheic nitrogen cycles focus on topographic factors and rely on numerical simulation. For example, Bardini et al. (2012) simulated the nitrogen cycles in the HZ driven by streambed dunes under steady flow conditions. The authors found that the availability of DOC determines whether the stream bed is the source or sink of nitrate. Bardini et al. (2013) further clarified that sediment heterogeneity has minimal effect on stream bed nutrient dynamics compared with stream bed morphology. Hu et al. (2014) emphasized that ignoring the micro-morphology in the large stream bed structure will underestimate the discharge of nitrate nitrogen from the HZ to the stream. Gu et al. (2007) studied the effects of flow and biogeochemical conditions on nitrate reduction in stream bed through indoor flume test and numerical simulation. They found that nitrate concentration of groundwater drainage is controlled by groundwater discharge rate and by the rates of aerobic respiration and denitrification.
However, studies on the mechanism of hyporheic nitrogen cycles driven by flood events are lacking. Gu et al. (2012) was the first to simulate the effects of flood waves on hyporheic exchange and nitrogen cycles. The lateral exchange in the riparian zone greatly improves the efficiency of hyporheic nitrogen removal. Furthermore, they revealed positive correlations between nitrogen removal and wave amplitude, duration, bank storage and maximum solute storage. However, their model did not consider nitrification, which might play a lead role in short hyporheic exchange paths and hydrological retention time (Zarnetske et al., 2011). Shuai et al. (2017) built an improved flow-solute transport model and explored the effects of surface-groundwater hydraulic gradient, aquifer hydraulic conductivity and dispersivity on hyporheic nitrogen removal under flood waves. The study’s results indicated that when surface water is in gaining condition, large aquifer hydraulic conductivity and dispersivity lead to great hyporheic nitrogen removal and efficiency. Similarly, Trauth et al. (2018) revealed that aerobic respiration and denitrification efficiency in the HZ increase correspondingly with the increase of amplitude and duration of flood wave through numerical simulations.
In summary, the understanding of mechanism of hyporheic nitrogen cycles driven by flood waves is unclear. Understanding of spatialtemporal distribution of hydrochemical characteristic quantities in a stream bed, riparian phreatic zone and variable saturated zone is lacking. Moreover, although previous studies fully discussed the influence of wave amplitude and duration on nitrogen cycles, the parameters characterizing the roundness and skewness of a wave remain unexamined. Hence, this study aims to solve the following three problems: (1) revealing the dynamics of reaction rates, solute consumption and consumption rates of hyporheic nitrogen cycles in response to a flood wave; (2) illustrating the variations of solute concentrations, solute plumes and reaction rates in stream-to-riparian continuum; and (3) systematically exploring the influence of four typical wave parameters on hyporheic nitrogen removal.
This study adopted the field test and flow model from Liu et al. (2019), based on which the solute transport and reaction model were coupled. The field site (128°12′18″E, 38°13′15″N) is located in the upstream of the Inbuk stream, Korea, where the stream stage amplitude and duration are able to reach 1 m and 5 days respectively due to seasonal flooding. Several piezometers were installed across the stream-to-riparian continuum at intervals of 10m for evaluating the hyporheic exchange between surface water and groundwater. And some chemical indexes in the surface water and groundwater were also measured.
The model coupled Richard’s Equation (Voss, 1984) and convection-dispersion-reaction Equation (1), as follows:
The nitrogen-cycling processes involved in this model included aerobic respiration, nitrification and denitrification. DOC was represented by the chemical formula ‘CH2O’, and all the reaction equations were as follows:
The Multiple-Monod kinetics model was used to characterize the aforementioned nitrogen-cycling processes (Molz et al., 1986).
The concentrations for O2, NO3- and NH4+ on the stream-aquifer boundary were set based on measured values of 5, 5 and 0.05 mg/L, respectively (Fig. 1). It was also set to 5 mg/L for DOC to compare better the degree of reaction of each solute owing to the lack of DOC concentration in surface water. Moreover, a conservative solute boundary condition (Na+) was also set, which was completely consistent with that of NO3-. The initial values of all solutes in the model were set to zero and were uniformly distributed. The longitudinal and transverse dispersivity of the model were set to 0.1 and 0.01m, respectively, according to the model scale and the empirical values. Other biochemical parameters were also determined by basing on the empirical values of previous similar studies (Table 1).
Conceptual flow-transport-reaction coupling model across the stream-to-riparian continuum. The dotted boundary denotes the submerged area. The blue points denote the observation points in stream bed, riparian phreatic zone and riparian variable saturated zone with the same distance from the stream-aquifer interface.
For multi-component reactive solutes (O2, NH4+, NO3- and DOC), the consumption of solute
The consumption rate (
The consumption efficiency (
where l is the length of the stream-aquifer interface.
In addition, the average reaction rate ( ) of solute
As Inbuk Stream is subjected to seasonal flooding and reservoir drainage throughout the year, the surface water level exhibits approximately sinusoidal fluctuation (Liu et al., 2019), as follows:
The range of each parameter was roughly determined based on the measured stream water level data from 2014 to 2016, and a variety of waveforms of different
Schematic diagram of a flood wave.
Fig. 3 showed the dynamics of average reaction rates, solute consumption and consumption rates in response to hyporheic exchange flux. The processes were roughly divided into three periods. (1) Surface-water-infiltration period (0.5 <
Change processes of surface water-groundwater exchange flux (a), average solute reaction rates (b), consumption rates (c) and consumption (d) within the study domain.
Fig. 4 showed the distributions of hydraulic head and solute concentrations (O2, NH4+, NO3-, DOC and Na+) across the stream-to-riparian continuum at typical moments of a flood wave (base case). During the rising limb (0.5 <
Distributions of hydraulic head and solute concentrations across the stream-to-riparian continuum at typical moments of a flood wave.
The shape of the solute plume in the stream-to-riparian continuum was similar to the envelope surface of tracer particles (Chen et al., 2003), which meant that in the vertical direction, the closer to the stream centre, the more shallowly the solute was transported. By comparison, in the lateral direction, the closer to the initial groundwater level, the more widely the solute was transported. To compare the nitrogen-cycling processes for different areas (stream bed, riparian phreatic zone and riparian variable saturated zone), observation points were set in each area. Fig. 5a-d showed that during the surface-water-infiltration period, the concentration for each solute in the riparian phreatic and riparian variable saturated zones almost increased in the same way. The increase in magnitude and speed were larger than that in the streambed. Similar variation trends were found for reaction rate of each chemical process (Fig. 5e-g). During the groundwater-backflow period, the decline rate of concentration in space from high to low was as follows: riparian phreatic zone >variable saturated zone >stream bed. However, the decline of reaction rate in space from high to low was as follows: variable saturated zone >riparian phreatic zone >stream bed.
Variations in solute concentration and reaction rate of each solute at different regional observation point units under a flood wave. SB, RPZ and RVSZ denote stream bed, riparian phreatic zone and riparian variable saturated zone, respectively.
Table 2 showed the variations of
Dimensionless relationship between hyporheic nitrogen removal (Mrem-NO3) and ‘A*T*tp/r’.
During the surface-water-infiltration period, the solute consumption and consumption rate increased with increasing solute reaction rate (Fig. 3b-d), but their trends gradually slowed down in the latter half. This trend arose because the intensity of recharge from surface water to aquifer gradually decreased (Fig. 3a), resulting in a slower rise in solute concentration, which was an important indicator affecting chemical reaction. Thus, with increasing concentration of each solute during the surface-water-infiltration period, the corresponding chemical reactions were strengthened, especially aerobic respiration. Although the reaction rate of NO3- increased during this period, the magnitude was relatively small due to the inhibition of rich O2. During the groundwater-backflow period, the reaction rates for DOC and O2 decreased significantly, whereas the reaction rate of NO3- decreased minimally. These findings indicated that with decreasing O2 concentration, its inhibition decreased, and the role of denitrification improved. From Fig. 3d, the consumption of each solute was constantly increasing during the flood event. Each solute consumed faster when surface water infiltrated into the aquifer, while it consumed more slowly when the groundwater flowed back. However, the consumption of each solute during the surface-water-infiltration period was smaller than that during the groundwater-backflow period. This result was due to the fact that the hydrological retention time during the groundwater-backflow period was much longer than that in the surface-water-infiltration period, which provided sufficient reaction time for each chemical process.
Generally, each solute plume expanded during the rising limb of surface water level and shrank during the falling limb. However, O2, NO3- and DOC plumes showed expansion trends during a preceding period of the falling limb, because the aquifer was still recharged by the surface water although the surface water level declined, and the solute amount being carried into the HZ was much more than the reaction removal. NO3- plume was the largest, because aerobic respiration and nitrification were strong at this period, whereas denitrification was limited by rich O2. Moreover, nitrification produced some extra NO3-. Although both O2 and DOC underwent aerobic respiration, less O2 participated in nitrification, whereas more DOC participated in denitrification, thereby resulting in the plume of O2 being slightly larger than that of DOC.
By comparing the variations in the concentration of O2, NO3- and DOC in three areas, we found that DOC concentration decreased the most, followed by O2 and then NO3-, which indicated that aerobic respiration was dominant in the whole process. This result was supported by comparing Fig. 5(f) and (g).The oxidation rate during a flood wave was about four times than that of denitrification rate. The reason the consumption of DOC was markedly higher than that of O2 was that the denitrification rate in which DOC was involved was significantly higher than the nitrification rate that O2 was involved in, which was supported by comparing Fig. 5(e) and (g). The denitrification rate was about thrice that of the nitrification rate. In addition, the overall reaction rate of each process in stream bed was smaller than that in the riparian zone, which was due to the relatively shallow but wide aquifer, resulting in a smaller amount of solute being carried into the vertical zone than that being carried into the lateral zone. In the riparian zone, the overall reaction rate of each process in the riparian variable saturated zone was smaller than that in the riparian phreatic zone, because a part of the solute in the variable saturated zone was absorbed by soil particles and remained there, thereby decreasing the overall solute concentrations within the region and reducing the reaction rates.
During a flood wave condition, nitrogen removal in the HZ (
This paper systematically studied the hyporheic nitrogen cycles driven by flood waves using a 2D flow-transport-reaction model. The main findings were as follows:
(1) When surface water infiltrated into the aquifer, the solute concentrations in the HZ increased accordingly, thereby resulting in gradual enhancements of aerobic respiration, nitrification and denitrification, of which aerobic respiration was dominant. When groundwater flowed back, all the chemical reaction processes were gradually weakened, of which the role of denitrification was improved.
(2) In the whole flood duration, solutes in the HZ were continuously consumed, with the consumption in the groundwater-backflow period being greater, although the corresponding reaction rates were slower. In space, the reaction rate of each chemical process from high to low was: riparian phreatic zone >variable saturated zone >stream bed.
(3) A good linear relation existed between the hyporheic nitrogen removal and the polynomial of each parameter ‘
This work is supported by the National Research Foundation of Korea (NRF), funded by the Korea government (MSIT) (No. NRF-2015R1A4A1041105), and partly funded by the key project of Nanjing Hydraulic Research Institute (NHRI) (No. Y919005).