Nitrate removal processes in the riverbed during a single-peak flood event

The effect of the interplay between unsteady flow and bedform in a flood event on nitrogen cycling in the hyporheic zone (HZ) remains poorly understood. In this study, a reactive transport groundwater model with different flood hydrographs was proposed to investigate the effect of modified hyporheic flow on nitrate dynamics in the HZ, including nitrate source-sink function, response to the single-peak flood event and removal efficiency. The results demonstrate that there exists an optimal range of river channel gradients that could enhance the biogeochemical reactions (respiration, nitrification and denitrification) in a flood event. The HZ acts as a nitrate sink especially after the flood event, and its source-sink function is independent of the unsteady discharge/stage conditions. The nitrate in the HZ has a hysteretic response to peak stage/discharge, and its removal efficiency is decreased by up to 70% compared to steady flow conditions. These findings not only provide a better understanding of nitrogen dynamics under the effect of unsteady channel flow, but also can be applied for river restoration to efficiently remove nitrate in the HZ by modifying river channel gradients.     

Algorithms for solving the diffusive wave flood routing equation

Generally, the diffusive wave equation, obtained by neglecting the acceleration terms in the Saint-Venant equations, is used in flood routing in rivers. Methods based on the finite-difference discretization techniques are often used to calculate discharges at each time step. A modified form of the diffusive wave equation has been developed and new resolution algorithms proposed which are better adapted to flood routing along a complex river network. The two parameters of the equation, celerity and diffusivity, can then be taken as functions of the discharge. The resolution algorithm allows the use of any distribution of lateral inflow in space and time. The accuracy of the new algorithms were compared with a traditional algorithm by numerical experimentation. Special attention was given to the instability caused by the inflow signal which constitutes the upstream boundary condition. For the fully diffusive wave flood routing problem, all three algorithms tested gave good results. The results also indicate that the efficiency of the new algorithms could be significantly improved if the position of the x-axis is modified by rotation. The new algorithms were applied to flood routing simulation over the Gardon d'Anduze catchment (542 km²) in southern France.     

Hyporheic exchange flows induced by constructed riffles and steps in lowland streams in southern Ontario,Canada

Stream–subsurface water interaction induced by natural riffles and constructed riffles/steps was examined in lowland streams in southern Ontario, Canada. The penetration of stream water into the subsurface was analysed using hydrometric data, and the zone of > 10% stream water was calculated from a chemical mixing equation using tracer injection of bromide and background chloride concentrations. The constructed riffles studied induced more extensive hyporheic exchange than the natural riffles because of their steeper longitudinal hydraulic head gradients and coarser streambed sediments. The depth of > 10% stream water zone in a small and a large constructed riffle extended to > 0·2 m and > 1·4 m depths respectively. Flux and residence time distribution of hyporheic exchange were simulated in constructed riffles using MODFLOW, a finite‐difference groundwater flow model. Hyporheic flux and residence time distribution varied along the riffles, and the exchange occurring upstream from the riffle crest was small in flux and had a long residence time. In contrast, hyporheic exchange occurring downstream from the riffle crest had a relatively short residence time and accounted for 83% and 70% of total hyporheic exchange flow in a small and large riffle respectively. Although stream restoration projects have not considered the hyporheic zone, our data indicate that constructed riffles and steps can promote vertical hydrologic exchange and increase the groundwater–surface water linkage in degraded lowland streams. Copyright © 2006 John Wiley & Sons, Ltd.     

Seasonal variation in cascade‐driven hyporheic exchange,northern Honduras

A characterization of hyporheic exchange for dry and wet season baseflow, as well as partially dewatered discharge, was done in Prieta Creek, a first‐order cascade in northern Honduras. The cascade had discharges from 1 to 15 l s^?1, had average slopes of 12%, pool spacing of 3 m, and shallow substrate of sand and gravel. Tracer tests were conducted in a 15‐m sub‐reach, a length considered to be adequate for the experiment based on the DaI test, a ratio of exchange and transport processes. In the three tests, between 9 and 18% of tracer was not recovered, possibly due to entrainment in flowpaths passing beneath the downstream monitoring location. Tracer data were analysed by the one‐dimensional transport with inflow and storage (OTIS) transient storage model (TSM) to derive standard exchange parameters, and by the solute transport in rivers (STIR) model to examine hyporheic residence time distributions (RTDs). The best fit of the observed tracer breakthrough curves was obtained by using the STIR model with a combination of two exponential RTDs to represent hyporheic retention. With increasing discharge, the OTIS model predicted increasing storage exchange fluxes and exchange coefficients and decreasing storage zone areas and transient storage times, which are trends supported by riparian and streambed piezometric head data. Riparian water levels rose during the transition from the dry to wet season, which could constrict the hyporheic storage zone. Thirteen of the 19 streambed piezometers recorded seasonal changes in hydraulic gradients and flux direction, with fewer yet stronger upwelling zones during higher discharges. The MODFLOW model missed the observed seasonal changes, possibly due to subtle changes in the seasonal change in water surface profiles. We conclude that partially dewatered dry season exchange, compared to wet season exchange, was initiated and terminated with smaller pressure gradients and, in different streambed locations, was smaller in volume, had longer residence times, and may connect with deeper and longer flow paths. Copyright © 2010 John Wiley & Sons, Ltd.     

In situ quantification of spatial and temporal variability of hyporheic exchange in static and mobile gravel‐bed rivers

Seepage meters modified for use in flowing water were used to directly measure rates of exchange between surface and subsurface water in a gravel‐ and cobble bed river in western Pennsylvania, USA (Allegheny River, Q_mean = 190 m³/s) and a sand‐ and gravel‐bed river in Colorado, USA (South Platte River, Q_mean = 9·7 m³/s). Study reaches at the Allegheny River were located downstream from a dam. The bed was stable with moss, algae, and river grass present in many locations. Median seepage was + 0·28 m/d and seepage was highly variable among measurement locations. Upward and downward seepage greatly exceeded the median seepage rate, ranging from + 2·26 (upward) to ? 3·76 (downward) m/d. At the South Platte River site, substantial local‐scale bed topography as well as mobile bedforms resulted in spatial and temporal variability in seepage greatly in exceedence of the median groundwater discharge rate of 0·24 m/d. Both upward and downward seepage were recorded along every transect across the river with rates ranging from + 2·37 to ? 3·40 m/d. Despite a stable bed, which commonly facilitates clogging by fine‐grained or organic sediments, seepage rates at the Allegheny River were not reduced relative to those at the South Platte River. Seepage rate and direction depended primarily on measurement position relative to local‐ and meso‐scale bed topography at both rivers. Hydraulic gradients were small at nearly all seepage‐measurement locations and commonly were not a good indicator of seepage rate or direction. Therefore, measuring hydraulic gradient and hydraulic conductivity at in‐stream piezometers may be misleading if used to determine seepage flux across the sediment‐water interface. Such a method assumes that flow between the well screen and sediment‐water interface is vertical, which appears to be a poor assumption in coarse‐grained hyporheic settings. Copyright © 2011 John Wiley & Sons, Ltd.     

Impact of a low-permeability lens on dune-induced hyporheic exchange

Hyporheic exchange induced by dunes is a key process controlling water fluxes and biogeochemical processes in river networks. Owing to the limitations of instrumental detection at small spatial scales, previous studies have focused mainly on dune-induced hyporheic exchange in homogeneous systems. A low-permeability lens is a natural, widespread heterogeneity in stream beds, and probably affects the processes of water flow and contaminant transportation significantly. To quantitatively analyse the response mechanism of hyporheic exchange to a low-permeability lens, a two-dimensional dune-generated hyporheic exchange model was developed using the VS2DH model. The results indicate a lens in a stream bed can hinder or enhance hyporheic exchange processes, depending on its relative spatial location to dunes. Both the increase in length and thickness of the lens could strengthen its impacts on hyporheic exchange. Regional groundwater–surface water interactions of higher intensity suppress the flow of hyporheic exchange in a stream bed with a low-permeability lens.     

复杂介质地震波传播的褶积微分算子数值模拟h

将佛尔塞（Forsyte）广义正交多项式微分算子地震波场正演模拟算法运用于复杂非均匀介质模型的波场数值模拟中,并比较了该方法在计算效率和计算精度方面与有限差分方法和伪谱法的差异. 数值结果表明,这种广义正交多项式微分算子法计算速度快、精度高,对计算资源需求低,是一种颇具潜力的数值模拟方法.      

Numerical simulation of seismic wave propagation in complex media by convolutional differentiator

We apply the forward modeling algorithm constituted by the convolutional Forsyte polynomial differentiator pro-posed by former worker to seismic wave simulation of complex heterogeneous media and compare the efficiency and accuracy between this method and other seismic simulation methods such as finite difference and pseudospec-tral method. Numerical experiments demonstrate that the algorithm constituted by convolutional Forsyte polyno-mial differentiator has high efficiency and accuracy and needs less computational resources, so it is a numerical modeling method with much potential.     

地震体波斜入射情形下台阶地形引起的波型转换

本文首次对地震体波斜入射情形下竖直、倾斜台阶地形引起的波型转换进行了数值模拟，结果表明：地震波的斜入射会使台阶上角点引起更强的转换面波；转换Rayleigh面波最大振幅可达弹性半空间表面自由场位移的1.1倍左右。     

Identification of groundwater exfiltration,interflow discharge,and hyporheic exchange flows by fibre optic distributed temperature sensing supported by electromagnetic induction geophysics

Water exchange across the sediment–water interface of streams impresses a characteristic thermal pattern at the interface. The use of fibre optic distributed temperature sensing at the sediment–water interface in a small sand‐bed stream identifies such temperature patterns. Groundwater and interflow can be differentiated based on the temporal evolution of temperature patterns. Additionally, sudden temperature changes at the sediment–water interface observed during the transit of floods enable spatial identification of local up and downwelling. Electromagnetic induction geophysics can detect subsurface texture structures that support groundwater–surface water exchange. Our results show that areas of permanent temperature anomalies observed with fibre optic distributed temperature sensing match areas of comparatively homogeneous electrical conductivity. This indicates groundwater discharge and enables differentiating groundwater discharge from interflow and local downwelling.     

Effects of river discharge on hyporheic exchange flows in salmon spawning areas of a large gravel‐bed river

The flow magnitude and timing from hydroelectric dams in the Snake River Basin of the Pacific north‐western US is managed in part for the benefit of salmon. The objective of this research was to evaluate the effects of Hells Canyon Dam discharge operations on hydrologic exchange flows between the river and riverbed in Snake River fall Chinook salmon spawning areas. Interactions between river water and pore water within the upper 1 m of the riverbed were quantified through the use of self‐contained temperature and water level data loggers suspended inside of piezometers. The data were recorded at 20 min intervals over a period of 200 days when the mean daily discharge was 218–605 m³ s^?1, with hourly stage changes as large as 1·9 m. Differences in head pressure between the river and riverbed were small, often within ± 2 cm. Measured temperature gradients in the riverbed indicated significant interactions between the surface and subsurface water. At the majority of sites, neither hydraulic nor temperature gradients were significantly affected by either short‐ or long‐term changes in discharge operations from Hells Canyon Dam. Only 2 of 14 study sites exhibited acute flux reversals between the river and riverbed resulting from short‐term, large magnitude changes in discharge. The findings suggest that local scale measurements may not be wholly explanatory of the hydrological exchange between the river and riverbed. The processes controlling surface water exchange at the study sites are likely to be bedform‐induced advective pumping, turbulence at the riverbed surface, and large‐scale hydraulic gradients along the longitudinal profile of the riverbed. By incorporating the knowledge of hydrological exchange processes into water management planning, regional agencies will be better prepared to manage the limited water resources among competing priorities that include salmon recovery, flood control, irrigation supply, hydropower production, and recreation. Copyright © 2007 John Wiley & Sons, Ltd.     

Reducing monitoring gaps at the aquifer–river interface by modelling groundwater–surface water exchange flow patterns

This study investigates spatial patterns and temporal dynamics of aquifer–river exchange flow at a reach of the River Leith, UK. Observations of sub‐channel vertical hydraulic gradients at the field site indicate the dominance of groundwater up‐welling into the river and the absence of groundwater recharge from surface water. However, observed hydraulic heads do not provide information on potential surface water infiltration into the top 0–15 cm of the streambed as these depths are not covered by the existing experimental infrastructure. In order to evaluate whether surface water infiltration is likely to occur outside the ‘window of detection’, i.e. the shallow streambed, a numerical groundwater model is used to simulate hydrological exchanges between the aquifer and the river. Transient simulations of the successfully validated model (Nash and Sutcliff efficiency of 0·91) suggest that surface water infiltration is marginal and that the possibility of significant volumes of surface water infiltrating into non‐monitored shallow streambed sediments can be excluded for the simulation period. Furthermore, the simulation results show that with increasing head differences between river and aquifer towards the end of the simulation period, the impact of streambed topography and hydraulic conductivity on spatial patterns of exchange flow rates decreases. A set of peak flow scenarios with altered groundwater‐surface water head gradients is simulated in order to quantify the potential for surface water infiltration during characteristic winter flow conditions following the observation period. The results indicate that, particularly at the beginning of peak flow conditions, head gradients are likely to cause substantial increase in surface water infiltration into the streambed. The study highlights the potential for the improvement of process understanding of hyporheic exchange flow patterns at the stream reach scale by simulating aquifer‐river exchange fluxes with a standard numerical groundwater model and a simple but robust model structure and parameterization. Copyright © 2011 John Wiley & Sons, Ltd.     

人工边界高频振荡失稳机理的一点注记

进一步完善了人工边界高频振荡失稳机理的解释. 解析证明了对于一维波动模型的离散模型, 当人工波速大于1.5倍的离散空间步距与时间步距的比值时, 则在某一高频段内, 其稳态波动解在人工边界上反射系数的模大于1的命题. 在此基础上对失稳的高频段作了进一步的讨论, 揭示了高频振荡失稳发生于对波动数值模拟无意义的高频段.      

Ocean wave imaging mechanism by imaging radar

Analytical representations of the high frequency spectra of ocean wave and its variation due to the variation of ocean surface current are derived from the wave-number spectrum balance equation. The ocean surface imaging formulation of real aperture radar (RAR) is given using electromagnetic wave backscattering theory of ocean surface and the modulations of ocean surface winds, currents and their variations to RAR are described. A general representation of the phase modulation induced by the ocean surface motion is derived according to standard synthetic aperture radar (SAR) imaging theory. The detectability of ocean current and sea bottom topography by imaging radar is discussed. The results constitute the theoretical basis for detecting ocean wave fields, ocean surface winds, ocean surface current fields, sea bottom topography, internal wave and so on.     

Implication of the flow resistance formulae on the prediction of flood wave propagation

ABSTRACT

This study examines the difference in the predictions of flood wave propagation in open channels depending on the flow resistance formulae, such as the Chézy and Manning’s equation. The celerity and diffusion coefficient are functions of the channel geometry, slope, roughness as well as the resistance formulae. The results suggest that substituting the Chézy equation with Manning’s equation results in different characteristics of flood propagation, which are consistent regardless of the cross-sectional geometry except for a circular cross-section: increasing celerity and decreasing diffusion coefficient. The celerity is more sensitive to the selection of resistance formulae than the diffusion coefficient. Geometry has a greater effect on the celerity and diffusion coefficient, and consequently on the resulting hydrographs. Manning’s equation results in a larger difference in celerity and diffusion coefficient compared to Chézy equation regardless of the water depth. Overall, this study shows that the selection of resistance formulae is important in terms of the resulting hydrographs and peak flow.

EDITOR Z.W. Kundzewicz ASSOCIATE EDITOR not assigned     

Vertical surface water–groundwater exchange processes within a headwater floodplain induced by experimental floods

Restoring hydrologic connectivity between channels and floodplains is common practice in stream and river restoration. Floodplain hydrology and hydrogeology impact stream hydraulics, ecology, biogeochemical processing, and pollutant removal, yet rigorous field evaluations of surface water–groundwater exchange within floodplains during overbank floods are rare. We conducted five sets of experimental floods to mimic floodplain reconnection by pumping stream water onto an existing floodplain swale. Floods were conducted throughout the year to capture seasonal variation and each involved two replicate floods on successive days to test the effect of varying antecedent moisture. Water levels and specific conductance were measured in surface water, soil, and groundwater within the floodplain, along with surface flow into and out of the floodplain. Vegetation density varied seasonally and controlled the volume of surface water storage on the floodplain. By contrast, antecedent moisture conditions controlled storage of water in floodplain soils, with drier antecedent moisture conditions leading to increased subsurface storage and slower flood wave propagation across the floodplain surface. The site experienced spatial heterogeneity in vertical connectivity between surface water and groundwater across the floodplain surface, where propagation of hydrostatic pressure, preferential flow, and bulk Darcy flow were all mechanisms that may have occurred during the five floods. Vertical connectivity also increased with time, suggesting higher frequency of floodplain inundation may increase surface water–groundwater exchange across the floodplain surface. Understanding the variability of floodplain impacts on water quality noted in the literature likely requires better accounting for seasonal variations in floodplain vegetation and antecedent moisture as well as heterogeneous exchange flow mechanisms. Copyright © 2016 John Wiley & Sons, Ltd.     

Effects of overlying velocity,particle size,and biofilm growth on stream–subsurface exchange of particles

Stream–subsurface exchange strongly influences the transport of contaminants, fine particles, and other ecologically relevant substances in streams. We used a recirculating laboratory flume (220 cm long and 20 cm wide) to study the effects of particle size, overlying velocity, and biofilm formation on stream–subsurface exchange of particles. Sodium chloride was used as a non‐reactive dissolved tracer and 1‐ and 5‐µm fluorescent microspheres were used as particulate tracers. Surface–subsurface exchange was observed with a clean sand bed and a bed colonized by an autotrophic–heterotrophic biofilm under two different overlying velocities, 0·9 and 5 cm s^?1. Hydrodynamic interactions between the overlying flow and sand bed resulted in a reduction of solute and particle concentrations in the water column, and a corresponding accumulation of particles in both the sediments and in the biofilm. Increasing overlying velocity and particle size resulted in faster removal from the overlying water due to enhanced mass transfer to the bed. The presence of the biofilm did not affect solute exchange under any flow condition tested. The presence of the biofilm significantly increased the deposition of particles under an overlying velocity of 5 cm s^?1, and produced a small but statistically insignificant increase at 0·9 cm^?1. The particles preferentially deposited within the biofilm matrix relative to the underlying sand. These results demonstrate that hydrodynamic transport conditions, particle size, and biofilm formation play a key role in the transport of suspended particles, such as inorganic sediments, particulate organic matter, and pathogenic microorganisms in freshwater ecosystems, and should be taken into consideration when predicting the fate and transport of particles and contaminants in the environment. Copyright © 2009 John Wiley & Sons, Ltd.     

Flood routing based on diffusion wave equation using mixing cell method

A one-dimensional non-linear diffusion wave equation is derived from the Saint Venant equations with neglect of the inertia terms. This non-linear equation has no general analytical solution. Numerical schemes are therefore employed to discretize the space and time axes and convert the differential equation to difference form. In this study, the mixing cell method is used to convert the diffusion wave equation to difference form, in which the difference term can be eliminated by selecting an optimal space step size Δx when time step size Δt is given. When the time step size Δt→0, the space step size Δx=Q/(2S₀BC]_k) where Q is discharge, S₀ is bed slope, B is channel width and C_k is kinematic wave celerity, which is the same as the characteristic length proposed by Kalinin and Milyukov. The results of application to two cases show that the mixing cell and linear channel flow routing methods produce hydrographs that are in agreement with the observed flood hydrographs. © 1997 John Wiley & Sons, Ltd.     

地震波散射理论及应用研究综述

回顾了地震波散射的概念、研究历史及分类.对散射理论及应用研究成果进行了系统评述，特别对前人关于散射波数值模拟和应用研究的不足进行了总结：(1)观测系统设计不合理；(2)对散射波的特征研究缺乏系统性；(3)对复杂散射波场的识别方法没有进行系统研究；(4)适合散射成像的处理软件还有待开发.并介绍了作者关于地震散射波场特征数值模拟研究的初步成果.     

爆炸地震波作用下地下结构动力响应数值分析

爆炸地震波荷载类似于天然地震波荷载，但又不完全相同。基于有效应力动力分析法，运用二维显式有限差分程序FLAC对地下结构在竖向和水平爆炸地震波荷载作用下的动力响应进行数值分析。编制了周围土体介质分析模型的程序模块并与FLAC接口。考虑了水平和竖向爆炸地震波荷载对地下结构的耦合效应，得出了一些定性的结论。