Modelling the effects of climatic and hydrological regime changes on the sediment dynamics of the Fraser River Basin,British Columbia,Canada

Future sediment dynamics may be affected by changing climates or hydrological regimes because of the close link between hydrology and sediment erosion, deposition, and transport. Previously, investigations of these potential changes have been constrained by a combination of limited observational data, hydrological drivers, and appropriate mechanistic models. Additionally, there is often ambiguity regarding how to disentangle the impacts of climate and hydrology from direct human factors such as reservoirs and land‐use change, which often exert more control over sediment dynamics. In this study, we utilize a recently developed, large‐scale, distributed, mechanistic sediment transport model to project future sediment erosion, deposition, and transportation within the Fraser River Basin in British Columbia, Canada—a basin with historical water flux and sediment load observations and limited anthropogenic influences upstream of its delta. The sediment model is driven by synthetic land‐surface hydrology derived from Scenarios A1B, A2, and B1 of the Special Report on Emissions Scenarios, which were provided by the Pacific Climate Impacts Consortium. Resulting simulations of water flux and sediment load from 1965 to 1994 are first validated against observational data then compared with future projections. Future projections show an overall increase in annual hillslope erosion and in‐channel transportation, a shift towards earlier spring peak erosion and transportation, and longer persistence of the sediment signal through the year. These shifts in timing and annual yield may have deleterious effects on spawning sockeye salmon and are insufficient to counteract future coastal retreat caused by sea‐level rise.     

Modelling of sediment dynamics in a laboratory‐scale experimental catchment

The variability of hillslope form and function is examined experimentally using a simple model catchment in which most landscape development parameters are either known or controlled. It is demonstrated that there is considerable variability in sediment output from similar catchments, subjected to the same hydrological processes, and for which the initial hillslope profiles are the same. The results demonstrate that, in the case of catchments with a linear initial hillslope profile, the sediment output is initially high but reduces through time, whereas for a concave initial profile the sediment output was smaller and relatively constant. Concave hillslope profiles also displayed reduced sediment output when compared with linear slopes with the same overall slope. Using this experimental model catchment data, the SIBERIA landscape evolution model was tested for its ability to predict temporal sediment transport. When calibrated for the rainfall and erodible material, SIBERIA is able to simulate mean temporal sediment output for the experimental catchment over a range of hillslope profiles and rainfall intensities. SIBERIA is also able to match the hillslope profile of the experimental catchments. The results of the study provide confidence in the ability of SIBERIA to predict temporal sediment output. The experimental and modelling data also demonstrate that, even with all geomorphic and hydrological variables being known and/or controlled, there is still a need for long‐term stream gauging to obtain reliable assessments of field catchment hydrology and sediment transport. Copyright © 2005 John Wiley & Sons, Ltd.     

A synthesis of morphological evolutions and Holocene stratigraphy of a wave-dominated estuary: The Arcachon lagoon,SW France

This work presents the first synthesis of secular to millenary morphological evolutions and stratigraphy of a wave-dominated estuary, the Arcachon lagoon, from a combination of unpublished bathymetric maps (1865 and 2001), core results and high-resolution seismic profiles recorded for the first time in this lagoon. The Arcachon lagoon is located on the Atlantic coast of France, facing the wave-dominated shelf of the Bay of Biscay. It is a mesotidal semi-enclosed environment of about 160 km².The sediment budget of the Arcachon lagoon was computed by subtracting the 1865 bathymetric map from that of 2001. The computed volume difference is low (?9.9±35×10⁶ m³ in 136 yrs) and is the result of the balance between erosion and accretion that occurs within tidal channels and tidal flats, respectively. This morphological evolution pattern is explained by low sediment supply and also by the tidal distortion resulting from the morphology of the lagoon. Deep channels connected to the inlet are dominated by ebb currents inducing erosion. Tidal flats and transverse channels display weak or flood-dominated tidal currents leading to the deposition of silts. The areas of tidal flat siltation locally correlate with the presence of oyster farms, suggesting the influence of Man on the lagoon sediment-fill. Transverse channel-infill is related to weak tidal currents resulting from the hydraulically inefficient orientation of these channels which served as an ancient drainage network.Evidence for tidal channel-infill and channel abandonment are also provided by seismic profiling and cores. The upper stratigraphic succession of the lagoon (about 10 m thick) includes four main stratigraphic units dominated by channel-fills. The two lower units (around 7500–2800 yrs BP) display tabular-shape sandy channels interpreted to be records of the open estuarine phase of the Arcachon lagoon. The two upper units (around 2800 yrs BP to present-day) display U-shaped mixed sand-and-mud channel-fills interpreted to be records of the closure of the lagoon. Given that the basal estuarine units are transgressive and the upper lagoonal units are regressive, the main stratigraphic change at around 2800 yrs BP is interpreted as being the maximum flooding surface (MFS). This late MFS is explained by the low sediment supply. It is proposed that the transition from the estuarine to the lagoonal phase is related to the development of the Cap-Ferret spit in response to an increase in the ratio between wave power to tide power. This change in wave-to-tide ratio may be triggered by wave power increase following the Subboreal/Subatlantic climate instability or a decrease in tide power following a decrease in tidal prism related to the lagoon sediment-fill.Thus, the evolution of the Arcachon lagoon over the last millenaries was mainly controlled by its spit development, leading to a wave-dominated estuary in terms of its geomorphology. Once it was partially closed, extensive mud flats developed in the lagoon which became ebb-dominated.     

Modelling sediment dynamics due to hillslope–river interactions: incorporating fluvial behaviour in landscape evolution model LAPSUS

Landscape evolution models (LEMs) simulate the three‐dimensional development of landscapes over time. Different LEMs have different foci, e.g. erosional behaviour, river dynamics, the fluvial domain, hillslopes or a combination. LEM LAPSUS is a relatively simple cellular model operating on timescales of centuries to millennia and using annual timesteps that has had a hillslope focus. Our objective was to incorporate fluvial behaviour in LAPSUS without changing the existing model equations. The model should be able to reproduce alternating aggradation and incision in the floodplains of catchments, depending on simulated conditions. Testing was done using an artificial digital elevation model (DEM) and a demonstration of the ability for fluvial simulation was performed for a real landscape (Torrealvilla catchment, southeast Spain). Model equations to calculate sediment dynamics and water routing were similar for both hillslope and fluvial conditions, but different parameter values were used for these domains, defined based on annual discharge. Parameters changing between the domains are convergence factor p, which is used in the multiple flow algorithm to route water, and discharge and gradient exponents m and n, used in transport capacity calculations. Erodibility and ‘sedimentability’ factors K and P were changed between cold (little vegetation, high erodibility) and warm conditions (more vegetation, lower erodibility). Results show that the adapted parameters reproduced alternating aggradation – due to divergent flow in the floodplain and sediment supply under cold conditions – and incision due to reduced sediment supply and resulting clean water erosion during simulated warm conditions. The simulated results are due to interactions between hillslopes and floodplains, as the former provide the sediments that are deposited in the latter. Similar behaviour was demonstrated when using the real DEM. Sensitivity and resolution analysis showed that the model is sensitive to changes in m, n and p and that model behaviour is influenced by DEM resolution. Copyright © 2012 John Wiley & Sons, Ltd.     

Assessing the potential underestimation of sediment and nutrient loads to the Great Barrier Reef lagoon during floods

Much of the sediment and nutrient load to the Great Barrier Reef (GBR) lagoon happens during over bank floods, when discharge can be significantly underestimated by standard river gauges. This paper assesses the potential need for a flood load correction for 28 coastal rivers that discharge into the GBR lagoon. For each river, daily discharge was divided into flows above and below a 'flood' threshold to calculate the mean annual percentage flow above this threshold. Most GBR rivers potentially need a flood load correction as over 15% of their mean annual flow occurs above the minor flood level; only seven rivers need little/no correction as their flood flows were less than 5% of the mean annual flow. Improved assessment of the true load of materials to the GBR lagoon would be an important contribution to the monitoring and reporting of progress towards Reef Plan and associated marine load targets.     

Modelling the transport of sediment discharged by Colombian rivers to the southern Caribbean Sea

Ocean Dynamics - The three-dimensional transport of sediments released by the main rivers in the Colombian Caribbean basin is investigated using numerical model simulations. Different types of...     

Historical development of the Venice lagoon contamination as recorded in radiodated sediment cores

Three cores (length 60 cm, i.d. 10 cm) collected in the most contaminated area of the lagoon of Venice (Italy) were sectioned (14, 12 and 15 sections each respectively) and analysed. Concentrations of heavy metals (Hg, Pb, Cd, Co, Cu, Zn, Cr, Ni, Fe), organic carbon, total nitrogen and total phosphorus were determined. Eh and pH profiles are discussed. Background (pre-industrial) levels for metals were estimated from the deepest sections. The principal component analysis performed on the metal concentrations, normalized on the backgrounds, revealed that two components can account for 99% of the variability. The first clearly represents the mercury pollution, the second weights the contribution of zinc and, to a lesser degree, of other metals (Cd, Cu, Co, Pb) to a contamination mainly caused by discharges of industrial tailings. Radiodating of two cores permitted us to reconstruct the historical evolution of the sediment pollution as related to the development of specific activities in the industrial zone of Porto Marghera. Total phosphorus and organic carbon accumulation was significant in the upper 10–20 cm of the cores.     

Introduction to special issue on connectivity in water and sediment dynamics

Connectivity has emerged in recent years as a significant conceptual framework within which to address the spatial and temporal variability in runoff and sediment transport. This special issue draws together several of the papers that were presented in the session “Connectivity in water and sediment dynamics: how do we move forwards?” at the 2012 General Assembly of the European Geosciences Union in Vienna, Austria. The papers submitted for this special issue fall into three groups: empirical studies, modelling studies, and conceptual development. Although it may be evident that the concept of connectivity helps us to express the complexity (in terms of water and sediment yields) of landscape responses to rainfall inputs, does it improve our ability to understand or predict those responses? There would still seem to be some way to go in connectivity research before this nagging concern can be assuaged. That it can be will undoubtedly be an important task for a number of ongoing research initiatives. Copyright © 2015 John Wiley & Sons, Ltd.     

Anthropogenic effects on sediment quality offshore southwestern Taiwan: Assessing the sediment core geochemical record

Slag material was dumped in two sites off southwestern Taiwan by the China Steel Corporation during 1984–1995. By geochemically analyzing four sediment cores, we investigated the impact of slag on the sediment chemistry. Elemental profiles from the dumping sites show a strong depletion in detrital elements and enrichment in alkaline and redox-sensitive elements, especially for the top ～20 cm of sediments. The relative enrichment factor (EF_rel) exhibits wider range for Ca (0.73–13.33), Fe (0.88–3.03), Mg (0.94–9.58) and Mn (1.22–33.30) due to contamination of sediments by slag. Sediment weathering indicators also show a distinct change with lower values in sediments influenced by dumping due to Ca and Sr addition. Higher EF_rel for As, Cd and Pb in the top sections of the cores indicate an additional modern/industrial input of these elements since ～1950. Our study suggests that slag dumping may have a localized effect on biogeochemical processes by inducing the diagenetic remobilization of certain elements.     

The large penumbra: Long-distance effects of artificial beach nourishment on Posidonia oceanica meadows

We used modelling and field approaches to assess the influence of beach nourishment on a relatively distant Posidonia oceanica seagrass meadow in the NW Mediterranean. Both sediment transport models and in situ observations showed that, two years after the nourishment and under wave climates prevalent during the study period, sediment movement was restricted to shallow waters (<7 m), above meadow shallow limit. The only significant impact on seagrass meadows during this period was an increase in fine sediments, associated with vertical rhizome growth rates 1.5–1.7 times higher than normal. However, running the model with data of wave climate over several decades indicates that strong storms able to transfer these sediments much deeper, potentially burying meadows, occur with a return period of about 25 years. Taken together, our results suggest that beach nourishment could result in significant sub-lethal and lethal consequences for seagrasses that may go unnoticed with short-term evaluations.     

Modelling to analyse the impacts of animal treading effects on soil infiltration

Animal treading can change soil physical properties, and thus is an important factor in hydrological modelling. We investigated the impacts of animal treading on infiltration by using a series of rainfall simulation experiments at Whatawhata Research Center, Waikato, New Zealand. The study identified significant variables for estimating soil steady‐state infiltration at a micro‐site (0·5 m²) and fitted the Green and Ampt equation by modifying or including variables for soil and water parameters and animal activities on grazing paddocks. A regression function for estimating steady‐state infiltration rate was created for each of four scenarios: between tracks (inter‐track), track, easy slope with ash soil, and easy slope with clay soil. Significant variables included the number of days after treading, antecedent soil moisture, field capacity, percentage of bare ground, bulk density, and the high degree of soil damage (damage not compacted). Regression models explained more than 71% of the variance in steady‐state infiltration for three scenarios, but only 53% for the easy slope with clay soil. The remodified Green and Ampt equation provided satisfactory estimation of infiltration for all scenarios (accuracy > 80%), and thus enables us to use the modified model for Waikato hill country pastures of different topography, soil physical condition, season and grazing management. Copyright © 2006 John Wiley & Sons, Ltd.     

Present-day sediment dynamics of the Sanaga catchment (Cameroon): from the total suspended sediment (TSS) to erosion balance

ABSTRACT

Total suspended sediment (TSS) data for 1960–1970 and from recent investigations (1990–2000) are used to evaluate the variability in sediment yield of the Sanaga catchment (Cameroon) and the equivalent rates of erosion. At the annual and seasonal time scales, total suspended sediment concentrations for the Mbam sub-catchment are three to four times higher than for the Sanaga basin, reflecting the higher sensitivity of the former to erosion. Classical clockwise hysteresis loops are observed in both sub-catchments, despite a significant increase in human-induced catchment changes. At the multi-year time frame (over a 40-year period), it appears that the second half of the 1990s is marked by a downward trend in TSS. This shift is due to the control infrastructures (river impoundments and ponds) installed in certain parts of the whole catchment, combined with a drop in annual rainfall and river discharge.

Editor Z.W. Kundzewicz; Associate editor not assigned     

Monitoring baseline suspended sediment in forested basins: the effects of sampling on suspended sediment rating curves

Abstract

Rating curves are widely used for directly assessing changes in the suspended sediment delivery process and indirectly for estimating total yields. Four sampling methods were simulated over a 31-day record of suspended sediment from the North Fork of the Mad River near Korbel, California. The position and size of the four groups of plotted slope/intercept pairs indicated differences in bias and variance among the methods. Estimates of total yield for the 31-day period and for storms of three sizes were also biased according to sampling method. A standard bias-correcting technique improved yield estimates, but did not remove sampling bias uniformly. Methods of data collection have a large and systematic effect on the estimation of rating-curve parameters and on estimates of suspended sediment yield. Differences attributed to land management may, in fact, result from changes in sampling methods.     

Interactive effects of environmental variability and human impacts on the long-term dynamics of an Amazonian floodplain lake and a South Atlantic coastal lagoon

Human activities are exposing freshwater ecosystems to a wide range of stressors, whose direct and indirect effects can be alleviated or exacerbated through interactive effects with dynamic environmental drivers. This study used long-term data from two Neotropical lacustrine freshwater systems (Batata Lake, an Amazonian floodplain lake and Imboassica lagoon, an Atlantic coastal lagoon) subjected to different kinds of environmental fluctuations (i.e., flood pulse and sandbar opening) and anthropogenic impacts (i.e., siltation and eutrophication). Our objective was to determine whether the effects of human perturbations are contingent on modifications of important biotic and abiotic characteristics through environmental variability. For both ecosystems, environmental variability consistently interacted with anthropogenic perturbations to alter most of the variables analyzed, such as nutrient dynamics, chlorophyll-a concentration, zooplankton and benthic invertebrate species richness, and temporal community stability, which indicates that interactive effects between environmental variability and anthropogenic perturbations may impact a myriad of ecosystem properties. Furthermore, the nature of these interactive effects was highly dependent on the variable considered and on the ecosystem analyzed. For example, at Imboassica lagoon, sandbar openings interacted synergistically with trophic state to increase the phosphorus concentration in the water column. At Batata Lake, flooding generally alleviated the negative effects of siltation on species richness by both diluting inorganic suspended material concentration and by promoting local recruitment from the regional species pool. Such results indicate that our ability to understand and predict the outcome of anthropogenic impacts on inland aquatic systems can be hampered if we consider human stressors as “static” phenomena disconnected from dynamic interactions with major local environmental drivers.     

Fluvial system dynamics derived from distributed sediment budgets: perspectives from an uncertainty‐bounded application

The utility of sediment budget analysis is explored in revealing spatio‐temporal changes in the sediment dynamics and morphological responses of a fluvial system subject to significant human impacts during the recent Anthropocene. Sediment budgets require a data‐intensive approach to represent spatially‐differentiated impacts adequately and are subject to numerous estimation uncertainties. Here, field and topographic surveys, historical data, numerical modelling and a representative‐area extrapolation method are integrated to construct a distributed, process‐based sediment budget that addresses historical legacy factors for the highly regulated Lagunitas Creek (213 km²), California, USA, for the period 1983–2008. Independent corroboration methods and error propagation analysis produce an uncertainty assessment unique to a catchment of this size. Current sediment yields of ~20 000 t a^‐1 ± 6000 t a^‐1 equate to unit rates of ~300 t km^‐2 a^‐1 ± 90 t km^‐2 a^‐1 over the effective sediment contributing area of 64 km². This is comparable with yields associated with early Euro‐American settlement in the catchment, despite loss of sediment supply upstream of the two large dams. It occurs because ~57% of the sediment is now derived from incision‐related channel erosion. Further, the highly efficient routing of channel‐derived sediments in these incised channels suggests an efflux of 84% of contemporary sediment production, contrasting with the efflux of ≈10–30% reported for unregulated agricultural catchments. The results highlight that sediment budgets for regulated rivers must accommodate channel morphological responses to avoid significantly misrepresenting catchment yields, and that volumetric precision in sediment budgets may best be improved by repeat, spatially dense, channel cross‐section surveys. Human activities have impacted every aspect of the sediment dynamics of Lagunitas Creek (production, storage, transfer, rates of movement through storage), confirming that, while distributed sediment budgets are data demanding and subject to numerous error sources, the approach can provide valuable insights into Anthropocene fluvial geomorphology. Copyright © 2017 John Wiley & Sons, Ltd.     

Modelling the dynamics and thermodynamics of volcanic degassing

 The rates of passive degassing from volcanoes are investigated by modelling the convective overturn of dense degassed and less dense gas-rich magmas in a vertical conduit linking a shallow degassing zone with a deep magma chamber. Laboratory experiments are used to constrain our theoretical model of the overturn rate and to elaborate on the model of this process presented by Kazahaya et al. (1994). We also introduce the effects of a CO₂–saturated deep chamber and adiabatic cooling of ascending magma. We find that overturn occurs by concentric flow of the magmas along the conduit, although the details of the flow depend on the magmas' viscosity ratio. Where convective overturn limits the supply of gas-rich magma, then the gas emission rate is proportional to the flow rate of the overturning magmas (proportional to the density difference driving convection, the conduit radius to the fourth power, and inversely proportional to the degassed magma viscosity) and the mass fraction of water that is degassed. Efficient degassing enhances the density difference but increases the magma viscosity, and this dampens convection. Two degassing volcanoes were modelled. At Stromboli, assuming a 2 km deep, 30% crystalline basaltic chamber, containing 0.5 wt.% dissolved water, the ∼700 kg s^–1 magmatic water flux can be modelled with a 4–10 m radius conduit, degassing 20–100% of the available water and all of the 1 to 4 vol.% CO₂ chamber gas. At Mount St. Helens in June 1980, assuming a 7 km deep, 39% crystalline dacitic chamber, containing 4.6 wt.% dissolved water, the ∼500 kg s^–1 magmatic water flux can be modelled with a 22–60 m radius conduit, degassing ∼2–90% of the available water and all of the 0.1 to 3 vol.% CO₂ chamber gas. The range of these results is consistent with previous models and observations. Convection driven by degassing provides a plausible mechanism for transferring volatiles from deep magma chambers to the atmosphere, and it can explain the gas fluxes measured at many persistently active volcanoes. Received: 26 September 1997 / Accepted: 11 July 1998     

Modelling of oxygen and nitrogen cycling as a function of macrophyte community in the Thau lagoon

A three-dimensional model coupling physical and biological processes for the whole Thau lagoon (Mediterranean coast of France) was developed in order to assess the relationships between macrophytes and the oxygen and nitrogen cycles. Ten species have been inserted as forcing variables in the model. Plankton dynamics, shellfish cultivation impact and mineralization of organic matter are also considered, as well as nutrient and oxygen exchanges between the sediment and the water column.Simulations with and without the macrophytes have shown that the system can be characterized as having a highly structured pattern involving lagoon nitrogen and oxygen cycles. This pattern is created by the combined influence of macrophytes, watershed and oyster farming.The model has been also used to assess the total annual macrophyte production at the whole lagoon scale. Comparisons with phytoplankton production and with results from other temperate lagoons have underlined the high productivity of the Thau lagoon supported by active nutrient regeneration.     

Post-glacial sediment dynamics in the Irish Sea and sediment wave morphology: Data–model comparisons

The irregular seafloor of the narrow Irish Sea on the NW European Shelf has been documented over several decades. From recently collected swath bathymetry data, very large trochoidal, nearly symmetrical sediment waves are observed in many parts of the Irish Sea and appear similar to those described from other continental shelf seas in North America that were covered by glacigenic sediments during the Last Glacial Maximum. Swath multibeam and single beam bathymetry data, backscatter intensity, shallow seismic imagery, video footage and sediment cores from the Irish Sea high sediment waves have been integrated to identify their genesis with reference to present and past hydrodynamic variability. From cross-sectional profiles over asymmetrical sediment waves in the Irish Sea the direction of asymmetry is used to map residual bed stress directions and associated bedload transport paths. Irish Sea peak bed stress vectors were generated using a two-dimensional palaeo-tidal model for the NW European shelf seas and compare well with the observations. Tidally induced bed stresses are modelled to have increased between 7–10 ka BP, to be nearly symmetrical in magnitude and to have reversed in dominant direction on a millennial scale. These environmental conditions during the post-glacial marine transgression are suggested here to help comprehend the construction of the very large sediment waves, with local variations due to differences in sediment grain size, sediment supply, water depth and intensified currents due to seafloor slopes. Model parameterisation using an open ocean boundary with time-dependent tidal changes and the implementation of high-resolution bathymetric information will improve future models of small-scale bed shear stress patterns and improve the predictive value of such modelling efforts.