Stormflow concentration–discharge dynamics of suspended sediment and dissolved phosphorus in an agricultural watershed

Concentration–discharge (C-Q) relationships are an effective tool for identifying watershed biogeochemical source and transport dynamics over short and long timescales. We examined stormflow C-Q, hysteresis, and flushing patterns of total suspended sediment (TSS) and soluble reactive phosphorus (SRP) in two stream reaches of a severely impaired agricultural watershed in northeastern Wisconsin, USA. The upper watershed reach—draining a relatively flat, row crop-dominated contributing area—showed predominantly anti-clockwise TSS hysteresis during storms, suggesting that particulate materials were mobilized more from distal upland sources than near- and in-channel areas. In contrast, the incised lower watershed reach produced strong TSS flushing responses on the rising limb of storm hydrographs and clockwise hysteresis, signalling rapid mobilization of near- and in-channel materials with increasing event flows. C-Q relationships for SRP showed complex patterns in both the upper and lower reaches, demonstrating largely non-linear chemodynamic C-Q behaviour during events. As with TSS, anti-clockwise SRP hysteresis in the upper reach suggested a delay in the hydrologic connectivity between SRP sources and the stream, with highly variable SRP concentrations during some events. A broad range of clockwise, anti-clockwise, and complex SRP hysteresis patterns occurred in the lower watershed, possibly influenced by in-channel legacy P stores and connection to tile drainage networks in the lower watershed area. Total suspended sediment and SRP responses were also strongly related to precipitation event characteristics including antecedent precipitation, recovery period, and precipitation intensity, highlighting the complexity of stormflow sediment and phosphorus responses in this severely impaired agricultural stream.     

Modelling catchment inflows into Lake Victoria: uncertainties in rainfall–runoff modelling for the Nzoia River

Abstract

Climate and soil characteristics vary considerably around the Lake Victoria basin resulting in high spatial and temporal variability in catchment inflows. However, data for estimating the inflows are usually sparsely distributed and error-prone. Therefore, modelled estimates of the flows are highly uncertain, which could explain early difficulties in reproducing the lake water balance. The aim of this study was to improve the estimates of catchment flow to Lake Victoria. The WASMOD model was applied to the Nzoia River, one of the major tributaries to Lake Victoria. Uncertainty was assessed within the GLUE framework. During calibration, log-transformation was performed on both simulated and observed flows. The results showed that, despite its simple structure, WASMOD produces acceptable results for the basin. For a Nash-Sutcliffe efficiency (NS) threshold of 0.6, the percentage of observations bracketed by simulations (POBS) was 74%, the average relative interval length (ARIL) was 0.93, and the maximum NS value was 0.865. The residuals were shown to be homoscedastic, normally distributed and nearly independent. When the NS threshold was increased to 0.8, POBS decreased to 54% with an improvement of ARIL to 0.49, highlighting the effect of the subjective choice of likelihood threshold.

Citation Kizza, M., Rodhe, A., Xu, C.-Y. & Ntale, H. K. (2011) Modelling catchment inflows into Lake Victoria: uncertainties in rainfall–runoff modelling for the Nzoia River. Hydrol. Sci. J. 56(7), 1210–1226.     

Rainfall–interception–evaporation–runoff relationships in a semi-arid catchment,northern Limpopo basin,Zimbabwe

Abstract

Two types of monthly water balance models at basin scale are used: PE models use precipitation and potential evapotranspiration (PET) as their observed input data, whereas P models need only precipitation. Calibration proceeds by comparing model runoff and observed runoff. Calibration is entirely automatic with the exclusion of subjective elements. All models differ only by their actual evapotranspiration equations. PE models from previous papers are generalized essentially by replacing the constant evapotranspiration parameter by a periodic one, thus increasing the number of parameters by two (a “parameter” is an unknown constant to be estimated, and which is a characteristic of the river basin to be described). P models use a periodic “driving force”, which is intended to represent periodicity of hydrological phenomena, normally originating in the (unavailable) PET time series. These eight PE models and three P models are then applied to 55 river basins in 10 countries with widely diverging climates and soil conditions. A marked improvement of model performance in about one third of the basins is due to the introduction of the above mentioned periodic functions. Even when PET data are available it is sometimes useful to consider P models. P models scarcely perform less well than PE models. An engineer, wanting to try out as few models as possible on a given river basin, can restrict his attention to the optimization of two or three models. The paper is an extension of a long effort towards monthly water balance models, and is believed to give a solution in most circumstances.     

Modelling of groundwater–vegetation interactions in a tidal marsh

Wetting and drying due to tidal fluctuations affect soil conditions and hence plant growth in tidal marshes. Here, a coupled one-dimensional model was developed to simulate interacting groundwater flow and plant growth in these wetlands. The simulation results revealed three characteristic zones of soil conditions for plant growth along a cross-creek section subjected to the combined influences of spring-neap tides and evapotranspiration: (1) a near-creek zone affected by semi-diurnal tides over the whole spring-neap cycle, where the soil is well aerated although the plant growth could be slightly limited by the local water content dropping periodically below the wilting point on the ebb tide; (2) a less well-drained zone where drainage occurs only during neap tides (for which the daily inundation is absent) and plant growth is aeration-limited; and (3) an interior zone where evapotranspiration determines the soil–water saturation. Plant growth dynamics, which depend on these soil conditions, lead to spatial biomass distributions that are consistent with the characteristic zonation. The simulations shed light on the feedback mechanism for groundwater–vegetation interactions in the marsh system. It was demonstrated that the growth of pioneer plants can improve the soil aeration condition as a result of transpiration. The strength of this feedback varies spatially in accordance with the three characteristic zones of soil–water saturation. However, the development of another species in the marsh system is likely to be more complicated than suggested by the “positive feedback” mechanism proposed previously, due to the influence of inter-species competition. The feedback effects are generally more complex, involving both plant growth enhancement and inhibition depending on the combined influence of the intra- and inter-species competition, the ecosystem’s carrying capacity and plant transpiration. These findings demonstrate the interplay of ecological and hydrological processes in tidal marshes, and provide guidance for future research, including field investigations that aim to establish the principle relationship between marsh morphology and plant zonation.     

Evaluation of storage–discharge relationships and recession analysis-based distributed hourly runoff simulation in large-scale,mountainous and snow-influenced catchment

ABSTRACT

Evaluation of a recession-based “top-down” model for distributed hourly runoff simulation in macroscale mountainous catchments is rare in the literature. We evaluated such a model for a 3090 km² boreal catchment and its internal sub-catchments. The main research question is how the model performs when parameters are either estimated from streamflow recession or obtained by calibration. The model reproduced observed streamflow hydrographs (Nash-Sutcliffe efficiency up to 0.83) and flow duration curves. Transferability of parameters to the sub-catchments validates the performance of the model, and indicates an opportunity for prediction in ungauged sites. However, the cases of parameter estimation and calibration excluding the effects of runoff routing underestimate peak flows. The lower end of the recession and the minimum length of recession segments included are the main sources of uncertainty for parameter estimation. Despite the small number of calibrated parameters, the model is susceptible to parameter uncertainty and identifiability problems.

EDITOR D. Koutsoyiannis; ASSOCIATE EDITOR A. Carsteanu     

Load or concentration,logged or unlogged? Addressing ten years of uncertainty in neural network suspended sediment prediction

This paper addresses current inconsistencies in methodological approaches for neural network modelling of suspended sediment. An expansion in the number of case studies being published over the last decade has yet to result in agreed guidelines on whether suspended sediment load or concentration should be modelled, and whether log‐transformation of data is either necessary or potentially beneficial. This contrasts with the well‐recognized guidelines that direct traditional sediment rating curve studies. The paper reports a comprehensive set of single‐input single‐output neural network suspended sediment modelling experiments performed on two catchments in Puerto Rico. It examines the impact of internal complexity, input variable choice and data transformation on the form, consistency and physical rationality of model outputs, the existence of localized overfitting and the usefulness of global performance metrics. Sound guidance on whether to model sediment load or concentration, and whether to model log‐transformed data is provided. Copyright © 2011 John Wiley & Sons, Ltd.     

Anhydrite in the 1989–1990 lavas and xenoliths from Redoubt Volcano,Alaska

The eruption of Redoubt Volcano in Alaska produced a moderate sulfur emission (estimated at 1 × 10 tons SO₂), but relatively small volume of lava (0.11 km ) with pre-eruption estimates of 840–950 °C and fO₂1.5 to 2.0 log units above NNO (Swanson, S.E., Nye, C.J., Miller, T.P., Avery, V.F., 1994. Magma mixing in the 1989–1990 eruption of Redoubt Volcano: Part II. Evidence from mineral and glass chemistry. Journal of Volcanology and Geothermal Research 62, 453–468). Petrologic estimates of sulfur production (Sigurdsson, H., Devine, J.D.,Davis, A.N., 1985. The petrologic estimation of volcanic degassing. Jokull 35, 1–8) from this eruption (Gerlach, T., Westrich, H.R., Casadevall, T.J., Finnegan, D.L., 1994. Vapor Saturation and accumulation in magmas of the 1989–1990 eruption of Redoubt Volcano, Alaska. Journal of Volcanology and Geothermal Research 62, 317–337) are considerably less than the measured sulfur emission, leading workers to propose the existence of a pre-eruption vapor phase to explain the “excess” sulfur.     

The relative impact of future storminess versus offshore dredging on suspended sediment concentration in a shallow coastal embayment: Rødsand lagoon,western Baltic Sea

Natural and human-induced disturbances challenge the resilience of coastal ecosystems. The relative impact of a projected increase in storm intensity versus a temporary sediment spill from the construction of the planned Fehmarnbelt Fixed Link on the suspended sediment concentrations was estimated in the shallow coastal embayment of Rødsand lagoon. The investigations were performed through a combination of field observations and numerical modeling. The numerical model was calibrated and validated with the field observations. The lagoonal response to a regional increase in wind speed showed temporarily amplified suspended particulate matter concentration (SPMC) levels (26–37%) and amplified exceedance duration above wash load (30–32%) in the entire lagoon, but especially in the shallow regions. The lagoonal response to a sediment spill outside Rødsand lagoon showed temporarily amplified SPMC levels (2–4%) and amplified exceedance duration above wash load (1–2%) in the entire lagoon, but especially in the regions in close proximity to the spill source. The regional response to increasing wind speed was thus 10-fold larger than the sediment spill where only advective sediment transport took place. This highlights the importance of the internal lagoonal processes as opposed to exterior disturbances. The study contributed to an improved understanding of the potential natural and human-induced stressors on the sediment dynamics in lagoonal ecosystems.

     

Magma–ice–sediment interactions and the origin of lava/hyaloclastite sequences in the Síða formation, South Iceland

Products of subglacial volcanism can illuminate reconstructions of paleo-environmental conditions on both local and regional scales. Competing interpretations of Pleistocene conditions in south Iceland have been proposed based on an extensive sequence of repeating lava-and-hyaloclastite deposits in the Síða district. We propose here a new eruptive model and refine the glacial environment during eruption based on field research and analytical data for the Síða district lava/hyaloclastite units. Field observations from this and previous studies reveal a repeating sequence of cogenetic lava and hyaloclastite deposits extending many kilometers from their presumed eruptive source. Glasses from lava selvages and unaltered hyaloclastites have very low H₂O, S, and CO₂ concentrations, indicating significant degassing at or close to atmospheric pressure prior to quenching. We also present a scenario that demonstrates virtual co-emplacement of the two eruptive products. Our data and model results suggest repeated eruptions under thin ice or partially subaerial conditions, rather than eruption under a thick ice sheet or subglacial conditions as previously proposed.     

Application of ~(137)Cs fingerprinting technique to interpreting sediment production records from reservoir deposits in a small catchment of the Hilly Loess Plateau, China

According to variations of 137Cs and clay contents, 44 flood couplets were identified in a profile of res- ervoir deposit with a vertical length of 28.12 m in the Yuntaishan Gully. Couplet 27 at the middle of the profile had the highest average 137Cs content of 12.65 Bq·kg?1, which indicated the 1963s' deposits, then 137Cs content decreased both downward and upward in the profile. The second top and bottom couplets had average 137Cs contents of 2.15 Bq·kg?1 and 0.92 Bq·kg?1, respectively. By integrated analysis of reservoir construction and management history, variations of 137Cs contents over the profile, sediment yields of flood couplets and rainfall data during the period of 1958-1970, individual storms related to the flood couplets were identified. 44 floods with a total sediment yield of 2.36×104 m3 occurred and flood events in a year varied between 1 and 10 times during the period of 1960-1970. 7-10 flood events occurred during the wet period of 1961-1964 with very wet autumn, while only 1-2 events during the dry period of 1965-1969. Average annual specific sediment yield was 1.29×104 t·km?2·a?1 for the Yuntaishan Gully during the period of 1960-1970, which was slightly higher than 1.11 ×104 t·km?2·a?1 for the Upper Yanhe River Basin above the Ganguyi Hydrological Station and slightly lower than 1.40 ×104 t·km?2·a?1 for the nearby Zhifang Gully during the same period. Annual specific sediment yields for the Yuntaishan Gully were correlated to the wet season's rainfalls well.     

Deformations related to a large (M=6.9) earthquake, the magma discharge, and eruptions in the Karymskii Volcanic Center in 1996–2005

A network of interconnected stations was established in the entire area of the Karymskii Volcanic Center and near the active Karymskii Volcano, Kamchatka in 1971–1988 for the purpose of studying ground deformation. Multiple observations by this network yielded quantitative characteristics of the ground deformation related to the following phenomena: the eruption of Karymskii Volcano during the periods 1976–1982 and January 1, 1996, to 2005 (still continuing, written in February 2008); the discharge of basalt on January 2, 1996, in the bottom of Lake Karymskii situated in the caldera of Akademii Nauk Volcano (this volcano had previously been thought to be extinct) and the subsequent phreatomagmatic eruption lasting approximately 24 hours; and the large (M 6.9) earthquake of January 1, 1996, occurring at 21 h 57 min local time in the Karymskii Volcanic Center at a depth of ～10 km. This paper discusses the relationships of ground deformation to volcanic activity and to the abovementioned unique natural occurrences, and their mechanism as deduced from geodetic data.     

Impact of climate and land-use changes on hydrological processes and sediment yield—a case study of the Be River catchment,Vietnam

Abstract

The impact of climate and land-use changes on hydrological processes and sediment yield is investigated in the Be River catchment, Vietnam, using the Soil and Water Assessment Tool (SWAT) hydrological model. The sensitivity analysis, model calibration and validation indicated that the SWAT model could reasonably simulate the hydrology and sediment yield in the catchment. From this, the responses of the hydrology and sediment to climate change and land-use changes were considered. The results indicate that deforestation had increased the annual flow (by 1.2%) and sediment load (by 11.3%), and that climate change had also significantly increased the annual streamflow (by 26.3%) and sediment load (by 31.7%). Under the impact of coupled climate and land-use changes, the annual streamflow and sediment load increased by 28.0% and 46.4%, respectively. In general, during the 1978–2000 period, climate change influenced the hydrological processes in the Be River catchment more strongly than the land-use change.

Editor Z.W. Kundzewicz; Associate editor Q. Zhang

Citation Khoi, D.N. and Suetsugi, T., 2014. Impact of climate and land-use changes on hydrological processes and sediment yield—a case study of the Be River catchment, Vietnam. Hydrological Sciences Journal, 59 (5), 1095–1108.     

Rainfall and streamflow response to El Niño Southern Oscillation: a case study in a semiarid catchment,Australia

Abstract

This paper aims to compare the shift in frequency distribution and skill of seasonal climate forecasting of both streamflow and rainfall in eastern Australia based on the Southern Oscillation Index (SOI) Phase system. Recent advances in seasonal forecasting of climate variables have highlighted opportunities for improving decision making in natural resources management. Forecasting of rainfall probabilities for different regions in Australia is available, but the use of similar forecasts for water resource supply has not been developed. The use of streamflow forecasts may provide better information for decision-making in irrigation supply and flow management for improved ecological outcomes. To examine the relative efficacy of seasonal forecasting of streamflow and rainfall, the shift in probability distributions and the forecast skill were evaluated using the Wilcoxon rank-sum test and the linear error in probability space (LEPS) skill score, respectively, at three river gauging stations in the Border Rivers Catchment of the Murray-Darling Basin in eastern Australia. A comparison of rainfall and streamflow distributions confirms higher statistical significance in the shift of streamflow distribution than that in rainfall distribution. Moreover, streamflow distribution showed greater skill of forecasting with 0–3 month lead time, compared to rainfall distribution.     

From rainfall to spring discharge: Coupling conduit flow,subsurface matrix flow and surface flow in karst systems using a discrete–continuum model

Physics-based distributed models for simulating flow in karst systems are generally based on the discrete–continuum approach in which the flow in the three-dimensional fractured limestone matrix continuum is coupled with the flow in discrete one-dimensional conduits. In this study we present a newly designed discrete–continuum model for simulating flow in karst systems. We use a flexible spatial discretization such that complicated conduit networks can be incorporated. Turbulent conduit flow and turbulent surface flow are described by the diffusion wave equation whereas laminar variably saturated flow in the matrix is described by the Richards equation. Transients between free-surface and pressurized conduit flow are handled by changing the capacity term of the conduit flow equation. This new approach has the advantage that the transients in mixed conduit flow regimes can be handled without the Preissmann slot approach. Conduit–matrix coupling is based on the Peaceman’s well-index such that simulated exchange fluxes across the conduit–matrix interface are less sensitive to the spatial discretization. Coupling with the surface flow domain is based on numerical techniques commonly used in surface–subsurface models and storm water drainage models. Robust algorithms are used to simulate the non-linear flow processes in a coupled fashion. The model is verified and illustrated with simulation examples.     

Modelling daily suspended sediment of rivers in Turkey using several data-driven techniques / Modélisation de la charge journalière en matières en suspension dans des rivières turques à l'aide de plusieurs techniques empiriques

Abstract

The transport of sediment load in rivers is important with respect to pollution, channel navigability, reservoir filling, longevity of hydroelectric equipment, fish habitat, river aesthetics and scientific interest. However, conventional sediment rating curves cannot estimate sediment load accurately. An adaptive neuro-fuzzy technique is investigated for its ability to improve the accuracy of the streamflow—suspended sediment rating curve for daily suspended sediment estimation. The daily streamflow and suspended sediment data for four stations in the Black Sea region of Turkey are used as case studies. A comparison is made between the estimates provided by the neuro-fuzzy model and those of the following models: radial basis neural network (RBNN), feed-forward neural network (FFNN), generalized regression neural network (GRNN), multi-linear regression (MLR) and sediment rating curve (SRC). Comparison of results reveals that the neuro-fuzzy model, in general, gives better estimates than the other techniques. Among the neural network techniques, the RBNN is found to perform better than the FFNN and GRNN.     

Variations in discharge and temperature of mineral springs at the Františkovy Lázně spa,Czech Republic,during a nearby earthquake swarm in 1985/1986

The western part of the Bohemian Massif (West-Bohemia/Vogtland region at the Czech-German border) is characterized by relatively frequent intraplate earthquake swarms and by other manifestations of present-day geodynamic activity. During the strong earthquake swarm at the turn of the years 1985 and 1986, significant changes in mineral spring parameters were observed at the spa of Františkovy Lázně. In this study, we present all available data on the mineral springs parameters, and we discuss them in terms of relations to seismic activity. Some changes in discharge were very distinct, amounting up to 40%, and had a co-seismic character. The changes in temperature were less noticeable, but preceded the beginning of the swarm by several months. Some hydrological changes persisted for nearly two years after the earthquake swarm. The character of the observed changes seems to support the hypothesis on an injection of mantle fluids, in particular of CO₂, as the main triggering mechanism of the earthquake swarm and the main cause of discharge anomalies.     

Evolution of a Gas–Dust Cloud in the Upper Atmosphere: Analysis of Events October 26, 2017

Geomagnetism and Aeronomy - A sequence of digital images of a gas–dust cloud of combustion products of solid fuel that formed during the separation of rocket stages is analyzed. The images...     

Hertz Range Pulsations during Recovery Phase of the Magnetic Storm on September 7–8, 2017, and Relation between their Dynamics and Changes in the Parameters of the Interplanetary Medium

Geomagnetism and Aeronomy - The shape and dynamics of an unusual disturbance in the Pc1 geomagnetic pulsation range recorded by ground-based induction magnetometers during the morning hours...     

Multiple volcano deformation sources in a post-rifting period: 1989–2005 behaviour of Krafla,Iceland constrained by levelling,tilt and GPS observations

The Krafla rifting episode, which occurred in North Iceland in 1975–1984, was followed by inflation of a shallow magma chamber until 1989. At that time, gradual subsidence began above the magma chamber and has continued to the present at a declining rate. Pressure decrease in a shallow magma chamber is not the only source of deformation at Krafla, as other deformation processes are driven by exploitation of two geothermal fields, together with plate spreading. In addition, deep-seated magma accumulation appears to take place, with its centre ∼ 10 km north of the Krafla caldera. The relative strength of these sources has varied with time. New results from a levelling survey and GPS measurements in 2005 allow an updated view on the deformation field. Deformation rates spanning 2000–2005 are the lowest recorded in the 30-year history of geodetic studies at the volcano. The inferred rate of 2000–2005 subsidence related to processes in the shallow magma chamber is less than 0.3 cm/yr whereas it was ∼ 5 cm/yr in 1989–1992. Currently, the highest rate of subsidence takes place in the Leirbotnar area, within the Krafla caldera, and appears to be a result of geothermal exploitation.