River flux uncertainties predicted by hydrological variability and riverine material behaviour

Data on riverine fluxes are essential for calculating element cycles (carbon, nutrients, pollutants) and erosion rates from regional to global scales. At most water‐quality stations throughout the world, riverine fluxes are calculated from continuous flow data (q) and discrete concentration data (C), the latter being the main cause of sometimes large uncertainties. This article offers a comprehensive approach for predicting the magnitude of these uncertainties for water‐quality stations in medium to large basins (drainage basin area > 1000 km²) based on the commonly used discharge‐weighted method. Uncertainty levels – biases and imprecisions – for sampling intervals of 3 to 60 days are correlated first through a nomograph with a flux variability indicator, the quantity of riverine material discharged in 2% of time (M_2%). In turn, M_2% is estimated from the combination of a hydrological reactivity index, W_2% (the cumulative flow volume discharged during the upper 2% of highest daily flow) and the truncated b_50sup exponent, quantifying the concentration versus discharge relationship for the upper half of flow values (C = a q ^b50sup, for q > q₅₀, where q₅₀ is the median flow): M_2% = W_2% + 27.6b_50sup. W_2% can be calculated from continuous flow measurements, and the b_50sup indicator can be calculated from infrequent sampling, which makes it possible to predict a priori the level of uncertainty at any station, for any type of riverine material either concentrated (b_50sup > 0) or diluted (b_50sup > 0) with flow. A large data base of daily surveys, 125 station variables of suspended particulate matter (SPM), total dissolved solids (TDS) and dissolved and particulate nutrients, was used to determine uncertainties from simulated discrete surveys and to establish relationships between indicators. Results show, for example, that for the same relatively reactive basin (W_2% > 25%), calculated fluxes from monthly sampling would yield uncertainties approaching ±100% for SPM (b_50sup > 1.4) fluxes and ±10% for TDS (b_50sup = ?0.2). The application to the nitrate survey of the river Seine shows significant trends for the 1972–2009 records. Copyright © 2012 John Wiley & Sons, Ltd.     

Choosing methods for estimating dissolved and particulate riverine fluxes from monthly sampling

Abstract

In discrete water quality surveys, riverine fluxes are associated with unknown uncertainties (biases and imprecisions). Annual flux errors have been determined from the generation of discrete surveys by Monte Carlo sorting for monthly sampling, from 10 years of daily records (120 records). Eight calculation methods were tested for suspended particulate matter, dissolved solids and dissolved and total nutrients in medium to large basins (10³ to 10⁶ km²) covering a wide range of hydrological conditions and riverine biogeochemistry. The performance of each method was analysed first by type of riverine material, which appeared to be much less pertinent than the flux variability matrix. The latter combines the river flow duration in two percent of time (W_2%) and the truncated exponent (b_50sup) defining the relationship of concentration vs discharge (C–Q) at higher flows (C = aQ^b50sup). As flux variability increases (high W_2% and/or high b_50sup), averaging and rating curve methods become less efficient compared to hydrograph separation methods. Flux biases and imprecisions were plotted in the [W_2%, b_50sup] matrix for discrete monthly surveys.

Editor Z. W. Kundzewicz

Citation Raymond, S., Moatar, F., Meybeck, M., and Bustillo, V., 2013. Choosing methods for estimating dissolved and particulate riverine fluxes from monthly sampling. Hydrological Sciences Journal, 58 (6), 1326–1339.     

Predicting suspended sediment concentrations in the Meuse river using a supply‐based rating curve

A rating curve provides a reasonable estimate of the suspended sediment concentration at a given discharge. However, analysis of a detailed 9‐year time‐series of suspended sediment concentration (SSC) and discharge Q of the Meuse River in The Netherlands indicates that SSC is (besides discharge) controlled by exhaustion and replenishment of different sediment sources. Clockwise hysteresis and other effects of sediment exhaustion can be observed during and after flood events, and the effects of stockpiling of sediment in the river bed during low‐discharge periods are obvious in the SSC of the next flood. In a single regression equation we have implemented a parameter that represents the presence or absence of stock for sediment uptake. In comparison with a rating curve of SSC and Q, adding this parameter is shown to be a more reliable and comprehensive method to predict SSCs at all discharge regimes with all preceding discharge conditions, for single‐peaked and multi‐peaked runoff events as well as for low flow conditions. The method is probably applicable to other small‐ to medium‐scaled river basins. Copyright © 2007 John Wiley & Sons, Ltd.     

Prediction of depth‐integrated fluxes of suspended sediment in the Amazon River: particle aggregation as a complicating factor

Large rivers have been previously shown to be vertically heterogeneous in terms of suspended particulate matter (SPM) concentration, as a result of sorting of suspended solids. Therefore, the spatial distribution of suspended sediments within the river section has to be known to assess the riverine sedimentary flux. Numerous studies have focused on the vertical distribution of SPM in a river channel from a theoretical or experimental perspective, but only a few were conducted so far on very large rivers. Moreover, a technique for the prediction of depth‐integrated suspended sediment fluxes in very large rivers based on sediment transport dynamics has not yet been proposed. We sampled river water along depth following several vertical profiles, at four locations on the Amazon River and its main tributaries and at two distinct water stages. Depending on the vertical profile, a one‐ to fivefold increase in SPM concentration is observed from river channel surface to bottom, which has a significant impact on the ‘depth‐averaged’ SPM concentration. For each cross section, a so‐called Rouse profile quantitatively accounts for the trend of SPM concentration increase with depth, and a representative Rouse number can be measured for each cross section. However, the prediction of this Rouse number would require the knowledge of the settling velocity of particles, which is dependent on the state of aggregation affecting particles within the river. We demonstrate that in the Amazon River, particle aggregation significantly influences the Rouse number and renders its determination impossible from grain‐size distribution data obtained in the lab. However, in each cross section, the Rouse profile obtained from the fit of the data can serve as a basis to model, at first order, the SPM concentration at any position in the river cross section. This approach, combined with acoustic Doppler current profiler (ADCP) water velocity transects, allows us to accurately estimate the depth‐integrated instantaneous sediment flux. Copyright © 2010 John Wiley & Sons, Ltd.     

Analysis of flow-sediment rating curve hysteresis based on flow and sediment travel time estimations

A flow-sediment rating curve is used to describe the relation between flow discharge and suspended-sediment concentration for a specific location. Five types of flow-sediment rating curves - single-valued line, clockwise loop, counterclockwise loop, single-valued line plus loop, and figure eight - were found to rely on the flow and available sediment arriving at the measuring site. In this study, equations for flow and sediment travel time were derived according to soil, rainfall, and watershed geomorphologic characteristics. The hysteresis of the rating curve was related to the travel times by a series of numerical tests. Field data collected from the Goodwin Creek Experimental Watershed, Mississippi, United States were used to verify the proposed rating curve hysteresis analysis. The results indicate that when the flow travel time is more extended than the sediment travel time, the rating curve shows a clockwise loop. A counterclockwise loop in the rating curve shows that the flow travel time is less extended than the sediment travel time. If the flow travel time exceeds the sediment travel time in specific runoff states and is less than the sediment travel time in other runoff states, then a single line plus a loop rating curve or a figure-eight rating curve is observed. The criterion for the model parameters to obtain equalization of the flow and sediment travel times was derived, which can identify the type of flow-sediment rating curve in a specific watershed.     

The effects of sample scheduling and sample numbers on estimates of the annual fluxes of suspended sediment in fluvial systems

Since the 1970s, there has been both continuing and growing interest in developing accurate estimates of the annual fluvial transport (fluxes and loads) of suspended sediment and sediment‐associated chemical constituents. This study provides an evaluation of the effects of manual sample numbers (from 4 to 12 year^?1) and sample scheduling (random‐based, calendar‐based and hydrology‐based) on the precision, bias and accuracy of annual suspended sediment flux estimates. The evaluation is based on data from selected US Geological Survey daily suspended sediment stations in the USA and covers basins ranging in area from just over 900 km² to nearly 2 million km² and annual suspended sediment fluxes ranging from about 4 Kt year^?1 to about 200 Mt year^?1. The results appear to indicate that there is a scale effect for random‐based and calendar‐based sampling schemes, with larger sample numbers required as basin size decreases. All the sampling schemes evaluated display some level of positive (overestimates) or negative (underestimates) bias. The study further indicates that hydrology‐based sampling schemes are likely to generate the most accurate annual suspended sediment flux estimates with the fewest number of samples, regardless of basin size. This type of scheme seems most appropriate when the determination of suspended sediment concentrations, sediment‐associated chemical concentrations, annual suspended sediment and annual suspended sediment‐associated chemical fluxes only represent a few of the parameters of interest in multidisciplinary, multiparameter monitoring programmes. The results are just as applicable to the calibration of autosamplers/suspended sediment surrogates currently used to measure/estimate suspended sediment concentrations and ultimately, annual suspended sediment fluxes, because manual samples are required to adjust the sample data/measurements generated by these techniques so that they provide depth‐integrated and cross‐sectionally representative data. Published 2014. This article is a U.S. Government work and is in the public domain in the USA.     

Monitoring urban impacts on suspended sediment,trace element,and nutrient fluxes within the City of Atlanta,Georgia, USA: program design,methodological considerations,and initial results

Atlanta, Georgia (City of Atlanta, COA), is one of the most rapidly growing urban areas in the US. Beginning in 2003, the US Geological Survey established a long‐term water‐quantity/quality monitoring network for the COA. The results obtained during the first 2 years have provided insights into the requirements needed to determine the extent of urban impacts on water quality, especially in terms of estimating the annual fluxes of suspended sediment, trace/major elements, and nutrients. During 2004/2005, suspended sediment fluxes from the City of Atlanta (COA) amounted to about 150 000 t year^?1; ≥ 94% of the transport occurred in conjunction with storm‐flow, which also accounted for ≥ 65% of the annual discharge. Typically, storm‐flow averaged ≤20% of the year. Normally, annual suspended sediment fluxes are determined by summing daily loads based on a single calculation step using mean‐daily discharge and a single rating curve‐derived suspended sediment concentration. Due to the small and ‘flashy’ nature of the COAs streams, this approach could produce underestimates ranging from 25% to 64%. Accurate estimates ( ± 15%) require calculation time‐steps as short as every 2–3 h. Based on annual median base‐flow/storm‐flow chemical concentrations, the annual fluxes of ≥ 75% of trace elements (e.g. Cu, Pb, Zn), major elements (e.g. Fe, Al), and total P occur in association with suspended sediment; in turn, ≥ 90% of the transport of these constituents occur in conjunction with storm‐flow. As such, base‐flow sediment‐associated and dissolved contributions represent relatively insignificant portions of the total annual load. An exception is total N, whose sediment‐associated fluxes range from 50% to 60%; even so, storm‐related transport typically exceeds 80%. Hence, in urban environments, non‐point‐sources appear to be the dominant contributors to the fluxes of these constituents. Published in 2007by John Wiley & Sons, Ltd.     

Effect of sediment resuspension on underwater light field in shallow lakes in the middle and lower reaches of the Yangtze River: A case study in Longgan Lake and Taihu Lake

Based on three continuous in situ underwater light field measurement under different wind waves conditions in Longgan Lake, Meiliang Bay of Taihu Lake in July 2003 and littoral zone near TLLER in July 2004, respectively, the effects of sediment resuspension caused by wind waves on PAR diffuse attenuation, absorption coefficients and euphotic depths are analyzed. In Longgan Lake, PAR diffuse attenuation coefficients during small, middle and large wind waves were 1.74, 2.02 and 2.45 m-1, respectively, and the corresponding PAR spectral diffuse attenuations ranged from 0.98 to 2.97, 1.34 to 3.95 and 1.80 to 5.40 m-1, respectively. In Meiliang Bay, PAR diffuse attenuation coefficients were 2.63, 3.72, 4.37 m-1 during small, middle and large wind waves. PAR diffuse attenuation coefficients increased by 41% and 66% from small to middle, large wind waves, respectively. Absorption coefficients integrated over the range of PAR of CDOM, phytoplankton were 0.26, 0.28 m-1; 0.76, 0.49 m-1, respectively during middle and large wind waves. Absorption coefficients integrated over the range of PAR of non-algal particulate matter and total suspended particulate matter increased from 0.94 to 1.73 m-1, and from 1.70 to 2.22 m-1, respectively during middle and large wind waves. Relative contributions of absorption coefficients of non-algal particulate matter to total absorption coefficient integrated over the range of PAR were 44.14%, 65.05%, respectively, during middle and large wind waves. PAR euphotic depths decreased by 0.40, 0.19, 0.20 m from middle to large wind waves in Longganhu Lake, Meliang Bay and littoral zone near TLLER. Significant correlations were found between transparency, PAR diffuse attenuation coefficients, euphotic depths and total suspended paniculate matter, wind velocity, wave height. Most significant correlations were found between transparency, PAR diffuse attenuation coefficients, euphotic depths and inorganic suspended paniculate matter but low correlations for chlorophyll a, dissolved organic carbon. Increase of total suspended paniculate matter, especially inorganic suspended paniculate matter caused by wind waves was the dominant factor affecting underwater light field in shallow lakes in the middle and lower reaches of the Yangtze River based on observations at three stations.     

Modeling the impact of wind and waves on suspended particulate matter fluxes in the East Frisian Wadden Sea (southern North Sea)

Suspended particulate matter (SPM) fluxes and dynamics are investigated in the East Frisian Wadden Sea using a coupled modeling system based on a hydrodynamical model [the General Estuarine Transport Model (GETM)], a third-generation wave model [Simulating Waves Nearshore (SWAN)], and a SPM module attached to GETM. Sedimentological observations document that, over longer time periods, finer sediment fractions disappear from the Wadden Sea Region. In order to understand this phenomenon, a series of numerical scenarios were formulated to discriminate possible influences such as tidal currents, wind-enhanced currents, and wind-generated surface waves. Starting with a simple tidal forcing, the considered scenarios are designed to increase the realism step by step to include moderate and strong winds and waves and, finally, to encompass the full effects of one of the strongest storm surges affecting the region in the last hundred years (Storm Britta in November 2006). The results presented here indicate that moderate weather conditions with wind speeds up to 7.5 m/s and small waves lead to a net import of SPM into the East Frisian Wadden Sea. Waves play only a negligible role during these conditions. However, for stronger wind conditions with speeds above 13 m/s, wind-generated surface waves have a significant impact on SPM dynamics. Under storm conditions, the numerical results demonstrate that sediments are eroded in front of the barrier islands by enhanced wave action and are transported into the back-barrier basins by the currents. Furthermore, sediment erosion due to waves is significantly enhanced on the tidal flats. Finally, fine sediments are flushed out of the tidal basins due to the combined effect of strong erosion by wind-generated waves and a longer residence time in the water column because of their smaller settling velocities compared to coarser sediments.

山东省东平湖菹草(Potamogeton crispus)腐烂过程中水体温室气体溶存浓度及扩散通量变化

为了解东平湖菹草(Potamogeton crispus)腐烂分解对水体温室气体溶存浓度和界面扩散通量的影响,于2016年5-7月在东平湖菹草腐烂期采集上覆水和沉积物柱样,测定上覆水和孔隙水中温室气体(N_2O、CH_4和CO_2)的溶存浓度,采用Fick第一定律和双层模型计算沉积物-水-气界面扩散通量,同时分析上覆水和沉积物的理化性质,并采用网袋分解法于现场进行菹草腐烂分解试验,以探究东平湖菹草腐烂过程中温室气体溶存和扩散的主要影响因子及其主要来源.结果表明,菹草腐烂符合二次指数模型,分为快速衰减和慢速分解两个阶段;菹草腐烂过程中上覆水pH和亚硝态氮浓度表现为先降低后升高,而溶解氧、氨氮、硝态氮和可溶性正磷酸盐浓度则为先升高后下降,沉积物中铵态氮含量表现为先升高后降低,硝态氮为先降低后显著升高,有机质和p H呈降低-升高-降低的波动变化;上覆水中各温室气体浓度和水气界面扩散通量均表现为CO_2 CH_4 N_2O,其扩散通量分别为5862.9±5441.4、31.15±41.3和0.15±0.57μmol/(m~2·h),整体表现为大气温室气体的"源",并以碳排放为主;上覆水中N_2O浓度和水-气界面扩散通量均先降低后升高,孔隙水中N_2O浓度在快速和慢速分解阶段分别出现极大值(22.7和55.6 nmol/L),而其沉积物-水界面通量前期持续增加至腐烂结束后迅速降低;上覆水和孔隙水中CH_4浓度及其各界面通量均表现为前期略有降低后持续升高;上覆水中CO_2浓度和水-气界面通量表现为持续升高后降低并趋于稳定,而孔隙水中CO_2呈波动变化,在菹草腐烂初期向孔隙水扩散,后期向上覆水扩散.水温是影响上覆水中温室气体浓度和水-气界面通量的主要因素;沉积物是水体N_2O和CH_4的主要来源,孔隙水中浓度是控制其沉积物-水界面扩散的重要因素;而上覆水中CO_2呈现多源性,但以上覆水中有机物质的矿化为主.     

Effects of climate-related variability in storage on streamwater solute concentrations and fluxes in a small forested watershed in the Southeastern United States

Streamwater quality can be affected by climate-related variability in hydrologic state, which controls flow paths and affects biogeochemical processes. Thirty-one years of input/output solute fluxes at Panola Mountain Research Watershed, a small, forested, seasonally water-limited watershed near Atlanta, Georgia, were used to quantify the effects of climatic-related variability in storage on streamwater solute concentrations and fluxes. Streamwater fluxes were estimated for ten solutes from weekly and event sample concentrations using regression-based methods. The most pertinent storage attribute (current or antecedent watershed, shallow, and deep storage) for each solute was determined by fitting separate concentration relationships. The concentration-discharge relationships varied more for reactive solutes such as potassium, sulfate, and DOC and less for weathering products (base cations and dissolved silica) and conservative chloride. Many solutes exhibited higher concentrations when storage levels were lower or wetting up, which was likely the result of the concentrating effects of evapotranspiration and/or the buildup and flushing of weathering products associated with longer residence times. The impacts of storage modeling on annual fluxes varied by solute, ranging from about 5% (magnesium) to 52% (nitrate) as relative standard deviations, and sufficiently removed climate-related patterns observed in streamwater concentrations. Sulfate was particularly mobilized following growing season droughts but only if deep storage was sufficiently recharged, possibly indicating that sulfides in the deep storage pool were oxidized to sulfate during droughts and mobilized when re-wetted. The lack of streamwater sulfate response to 61% declines in atmospheric deposition indicates the importance of watershed biogeochemical processes on controls of streamwater export of sulfate. The approach of explicitly incorporating storage in the streamwater concentration modeling elucidated the effects of climate on streamwater water-quality and may provide insight into the effects of climatic change on future fluxes.     

Lithological and hydrological controls on water composition: evaporite dissolution and glacial weathering in the south central Andes of Argentina (33°–34°S)

Lithological and hydrological influence on fluvial physical and chemical erosion was studied in a glacierized sedimentary basin with high evaporite presence. Suspended particulate matter (SPM), total dissolved solids (TDS) and major ion concentrations were analysed for 2 years of different hydrologic condition: (i) 2009–2010, Q = 100% average; and (ii) 2010–2011, Q = 60% average. Annual hydrograph was simple regime‐type with one peak in summer related to snow melting. The intra‐annual SPM and TDS variations were directly and inversely associated to Q, respectively. Snow chemistry showed continental influence (Na⁺/Ca²⁺ = 0.17), and atmospheric input of TDS was <1% of the total exported flux. River water was highly concentrated in Ca²⁺ and SO₄²⁻ (~4 mmol l⁻¹) and in Na⁺ and Cl⁻ (~3 mmol l⁻¹). Ca²⁺/SO₄²⁻ and Na⁺/Cl⁻ molar ratios were ~1 and related to Q, directly and inversely, respectively. Major ion relationships suggest that river chemistry is controlled by evaporite (gypsum and halite) dissolution having a summer input from sulfide oxidation and carbonate dissolution, and a winter input from subsurface flow loaded with silicate weathering products. This variation pattern resulted in nearly chemostatic behaviour for Ca⁺, Mg²⁺ and SO₄²⁻, whereas Na⁺, Cl⁻ and SiO₂ concentrations showed to be controlled by dilution/concentration processes. During the 2009–2010 hydrological year, the fluxes of water, SPM and TDS registered in the snow melting–high Q season were, respectively, 71%, 92% and 67% of the annual total, whereas for equal period in 2010–2011, 56% of water, 86% of SPM and 54% of TDS annual fluxes were registered. The SPM fluxes for 2009–2010 and 2010–2011 were 1.19 × 10⁶ and 0.79 × 10⁶ t year⁻¹, whereas TDS fluxes were 0.68 × 10⁶ and 0.55 × 10⁶ t year⁻¹, respectively. Export rates for 2009–2010 were 484 t km² year⁻¹ for SPM and 275 t km² year⁻¹ for TDS. These rates are higher than those observed in glacierized granite basins and in non‐glacierized evaporite basins, suggesting a synergistic effect of lithology and glaciers on physical and chemical erosion. Copyright © 2014 John Wiley & Sons, Ltd.     

Influence of rainfall and soil properties spatial aggregation on extreme flash flood response modelling: An evaluation based on the Sesia river basin,North Western Italy

High resolution radar rainfall fields and a distributed hydrologic model are used to evaluate the sensitivity of flood and flash flood simulations to spatial aggregation of rainfall and soil properties at catchment scales ranging from 75 to 983 km². Hydrologic modeling is based on a Hortonian infiltration model and a network-based representation of hillslope and channel flow. The investigation focuses on three extreme flood and flash flood events occurred on the Sesia river basin, North Western Italy, which are analysed by using four aggregation lengths ranging from 1 to 16 km. The influence of rainfall spatial aggregation is examined by using the flow distance as a spatial coordinate, hence emphasising the role of river network in the averaging of space–time rainfall. The effects of reduced and distorted rainfall spatial variability on peak discharge have been found particularly severe for the flash flood events, with peak errors up to 35% for rainfall aggregation of 16 km and at 983 km² catchment size. Effects are particularly remarkable when significant structured rainfall variability combines with relatively important infiltration volumes due to dry initial conditions, as this emphasises the non-linear character of the rainfall–runoff relationship. In general, these results confirm that the correct estimate of rainfall volume is not enough for the accurate reproduction of flash flood events characterised by large and structured rainfall spatial variability, even at catchment scales around 250 km². However, accurate rainfall volume estimation may suffice for less spatially variable flood events. Increasing the soil properties aggregation length exerts similar effects on peak discharge errors as increasing the rainfall aggregation length, for the cases considered here and after rescaling to preserve the rainfall volume. Moreover, peak discharge errors are roughly proportional to runoff volume errors, which indicates that the shape of the flood wave is influenced in a limited way by modifying the detail of the soil property spatial representation. Conversely, rainfall aggregation may exert a pronounced influence on the discharge peak by reshaping the spatial organisation of the runoff volumes and without a comparable impact on the runoff volumes.