Sediment dispersal from a typical Mediterranean flood: The Têt River,Gulf of Lions

This paper describes an integrated study of a typical Mediterranean flood event in the Gulf of Lions. A flood with a 5-year return interval occurred in the Têt River basin and adjacent inner-shelf in the Gulf of Lions, northwest Mediterranean, during April 2004. Data were collected during this flood as part of event-response investigations of the EU-funded Eurostrataform (European Margin Strata Formation) project. Southeasterly storm winds led to a flood which directly modified the inner-shelf hydrodynamics. Sediment delivery to the coastal zone during this flood represented more than half of the mean annual discharge of the Têt River to the Gulf of Lions. This river transported a large amount of sand in suspension, representing 25% of the total suspended load, and as bedload representing 8% of the total load, during this event. Sand introduced in the nearshore was transported northwards during the peak storm and nourished a small delta. Fine sediments were separated from coarse sediments at the river mouth, and were advected southwards and seawards by the counter-clockwise general circulation. Fine-grained sediments were transported via a hypopycnal plume along the coast towards the southern tip of the Gulf of Lions and the Cap Creus canyon. The along-shore currents, which intensified from north to south of the Gulf of Lions, particularly between the Cap Creus promontory and the Cap Creus canyon, favoured the transfer of fine-grained sediments from the continental shelf of the Gulf of Lions towards the continental slope. Our results show that floods with a few-year return interval in small coastal rivers can play a significant role in the transport of sediments on microtidal continental margins and their export from the shelf through canyons.     

Comparison of horizontal and downward particle fluxes across canyons of the Gulf of Lions (NW Mediterranean): Meteorological and hydrodynamical forcing

In the framework of the European project Eurostrataform, an array of six near-bottom mooring lines was deployed in the heads of the major submarine canyons incising the continental slope of the Gulf of Lions (NW Mediterranean). All moorings were equipped with sediment traps, current meters and turbidity Optical Backscatters Sensors (OBS) situated at few meters above the bottom. Particulate mass fluxes were recorded weekly by the sediment traps between November 2003 and May 2004 and compared with horizontal export fluxes obtained from the turbidity and current meters.     

Calibration of subsurface batch and reactive-transport models involving complex biogeochemical processes

In this study, the calibration of subsurface batch and reactive-transport models involving complex biogeochemical processes was systematically evaluated. Two hypothetical nitrate biodegradation scenarios were developed and simulated in numerical experiments to evaluate the performance of three calibration search procedures: a multi-start non-linear regression algorithm (i.e. multi-start Levenberg–Marquardt), a global search heuristic (i.e. particle swarm optimization), and a hybrid algorithm that combines the particle swarm procedure with a regression-based “polishing” step. Graphical analysis of the selected calibration problems revealed heterogeneous regions of extreme parameter sensitivity and insensitivity along with abundant numbers of local minima. These characteristics hindered the performance of the multi-start non-linear regression technique, which was generally the least effective of the considered algorithms. In most cases, the global search and hybrid methods were capable of producing improved model fits at comparable computational expense. In other cases, the multi-start and hybrid calibration algorithms yielded comparable fitness values but markedly differing parameter estimates and associated uncertainty measures.     

Higher and lowest order mixed finite element approximation of subsurface flow problems with solutions of low regularity

In this work we study mixed finite element approximations of Richards’ equation for simulating variably saturated subsurface flow and simultaneous reactive solute transport. Whereas higher order schemes have proved their ability to approximate reliably reactive solute transport (cf., e.g. [Bause M, Knabner P. Numerical simulation of contaminant biodegradation by higher order methods and adaptive time stepping. Comput Visual Sci 7;2004:61–78]), the Raviart–Thomas mixed finite element method (RT₀) with a first order accurate flux approximation is popular for computing the underlying water flow field (cf. [Bause M, Knabner P. Computation of variably saturated subsurface flow by adaptive mixed hybrid finite element methods. Adv Water Resour 27;2004:565–581, Farthing MW, Kees CE, Miller CT. Mixed finite element methods and higher order temporal approximations for variably saturated groundwater flow. Adv Water Resour 26;2003:373–394, Starke G. Least-squares mixed finite element solution of variably saturated subsurface flow problems. SIAM J Sci Comput 21;2000:1869–1885, Younes A, Mosé R, Ackerer P, Chavent G. A new formulation of the mixed finite element method for solving elliptic and parabolic PDE with triangular elements. J Comp Phys 149;1999:148–167, Woodward CS, Dawson CN. Analysis of expanded mixed finite element methods for a nonlinear parabolic equation modeling flow into variably saturated porous media. SIAM J Numer Anal 37;2000:701–724]). This combination might be non-optimal. Higher order techniques could increase the accuracy of the flow field calculation and thereby improve the prediction of the solute transport. Here, we analyse the application of the Brezzi-Douglas-Marini element (BDM₁) with a second order accurate flux approximation to elliptic, parabolic and degenerate problems whose solutions lack the regularity that is assumed in optimal order error analyses. For the flow field calculation a superiority of the BDM₁ approach to the RT₀ one is observed, which however is less significant for the accompanying solute transport.     

Storm-driven sediment transport in Massachusetts Bay

Massachusetts Bay is a semi-enclosed embayment in the western Gulf of Maine about 50 km wide and 100 km long. Bottom sediment resuspension is controlled predominately by storm-induced surface waves and transport by the tidal- and wind-driven circulation. Because the Bay is open to the northeast, winds from the northeast (‘Northeasters’) generate the largest surface waves and are thus the most effective in resuspending sediments. The three-dimensional oceanographic circulation model Regional Ocean Modeling System (ROMS) is used to explore the resuspension, transport, and deposition of sediment caused by Northeasters. The model transports multiple sediment classes and tracks the evolution of a multilevel sediment bed. The surficial sediment characteristics of the bed are coupled to one of several bottom-boundary layer modules that calculate enhanced bottom roughness due to wave–current interaction. The wave field is calculated from the model Simulating WAves Nearshore (SWAN). Two idealized simulations were carried out to explore the effects of Northeasters on the transport and fate of sediments. In one simulation, an initially spatially uniform bed of mixed sediments exposed to a series of Northeasters evolved to a pattern similar to the existing surficial sediment distribution. A second set of simulations explored sediment-transport pathways caused by storms with winds from the northeast quadrant by simulating release of sediment at selected locations. Storms with winds from the north cause transport southward along the western shore of Massachusetts Bay, while storms with winds from the east and southeast drive northerly nearshore flow. The simulations show that Northeasters can effectively transport sediments from Boston Harbor and the area offshore of the harbor to the southeast into Cape Cod Bay and offshore into Stellwagen Basin. This transport pattern is consistent with Boston Harbor as the source of silver found in the surficial sediments of Cape Cod Bay and Stellwagen Basin.     

A transport‐distance approach to scaling erosion rates: 3. Evaluating scaling characteristics of Mahleran

In the two previous papers of this series, we demonstrated how a novel approach to erosion modelling (Mahleran – Model for Assessing Hillslope‐Landscape Erosion, Runoff And Nutrients) provided distinct advantages in terms of process representation and explicit scaling characteristics when compared with existing models. A first evaluation furthermore demonstrated the ability of the model to reproduce spatial and temporal patterns of erosion and their particle‐size characteristics on a large rainfall‐simulation plot. In this paper, we carry out a more detailed evaluation of the model using monitored erosion events on plots of different size. The evaluation uses four plots of 21·01, 115·94, 56·84 and 302·19 m², with lengths of 4·12, 14·48, 18·95 and 27·78 m, respectively, on similar soils to the rainfall‐simulation plot, for which runoff and erosion were monitored under natural rainfall. Although the model produces the correct ranking of the magnitude of erosion events, it performs less well in reproducing the absolute values and particle‐size distributions of the eroded sediment. The implications of these results are evaluated in terms of requirements for process understanding and data for parameterization of improved soil‐erosion models. We suggest that there are major weaknesses in the current understanding and data underpinning existing models. Consequently, a more holistic re‐evaluation is required that produces functional relationships for different processes that are mutually consistent, and that have appropriate parameterization data to support their use in a wide range of environmental conditions. Copyright © 2008 John Wiley & Sons, Ltd.     

One-dimensional numerical simulation of non-uniform sediment transport under unsteady flows

One-dimensional numerical models are popularly used in sediment transport research because they can be easily programmed and cost less time compared with two- and three-dimensional numerical models. In particular, they possess greater capacity to be applied in large river basins with many tributaries. This paper presents a one-dimensional numerical model capable of calculating total-load sediment transport. The cross-section-averaged sediment transport capacity and recovery coefficient are addressed in the suspended load model. This one-dimensional model, therefore, can be applied to fine suspended loads and to hyperconcentrated flows in the Yellow River. Moreover, a new discretization scheme for the equation of unsteady non-uniform suspended sediment transport is proposed. The model is calibrated using data measured from the Yantan Reservoir on the Hongshui River and the Sanmenxia Reservoir on the Yellow River. A comparison of the calculated water level and river bed deformation with field measurements Shows that the improved numerical model is capable of predicting flow, sediment transport, bed changes, and bed-material sorting in various situations, with reasonable accuracy and reliability.     

Solute transport by surface runoff from low‐angle slopes: theory and application

The removal of chemicals in solution by overland flow from agricultural land has the potential to be a significant source of chemical loss where chemicals are applied to the soil surface, as in zero tillage and surface‐mulched farming systems. Currently, we lack detailed understanding of the transfer mechanism between the soil solution and overland flow, particularly under field conditions. A model of solute transfer from soil solution to overland flow was developed. The model is based on the hypothesis that a solute is initially distributed uniformly throughout the soil pore space in a thin layer at the soil surface. A fundamental assumption of the model is that at the time runoff commences, any solute at the soil surface that could be transported into the soil with the infiltrating water will already have been convected away from the area of potential exchange. Solute remaining at the soil surface is therefore not subject to further infiltration and may be approximated as a layer of tracer on a plane impermeable surface. The model fitted experimental data very well in all but one trial. The model in its present form focuses on the exchange of solute between the soil solution and surface water after the commencement of runoff. Future model development requires the relationship between the mass transfer parameters of the model and the time to runoff to be defined. This would enable the model to be used for extrapolation beyond the specific experimental results of this study. The close agreement between experimental results and model simulations shows that the simple transfer equation proposed in this study has promise for estimating solute loss to surface runoff. Copyright © 2000 John Wiley & Sons, Ltd.     

Rainfall runoff and erosion in Napa Valley vineyards: effects of slope,cover and surface roughness

The effects of slope, cover and surface roughness on rainfall runoff, infiltration and erosion were determined at two sites on a hillside vineyard in Napa County, California, using a portable rainfall simulator. Rainfall simulation experiments were carried out at two sites, with five replications of three slope treatments (5%, 10% and 15%) in a randomized block design at each site (0%bsol;64 m² plots). Prior to initiation of the rainfall simulations, detailed assessments, not considered in previous vineyard studies, of soil slope, cover and surface roughness were conducted. Significant correlations (at the 95% confidence level) between the physical characteristics of slope, cover and surface roughness, with total infiltration, runoff, sediment discharge and average sediment concentration were obtained. The extent of soil cracking, a physical characteristic not directly measured, also affected analysis of the rainfall–runoff–erosion process. Average cumulative runoff and cumulative sediment discharge from site A was 87% and 242% greater, respectively, than at site B. This difference was linked to the greater cover, extent of soil cracking and bulk density at site B than at site A. The extent of soil cover was the dominant factor limiting soil loss when soil cracking was not present. Field slopes within the range of 4–16%, although a statistically significant factor affecting soil losses, had only a minor impact on the amount of soil loss. The Horton infiltration equation fit field data better than the modified Philip's equation. Owing to the variability in the ‘treatment’ parameters affecting the rainfall–runoff–erosion process, use of ANOVA methods were found to be inappropriate; multiple‐factor regression analysis was more useful for identifying significant parameters. Overall, we obtained similar values for soil erosion parameters as those obtained from vineyard erosion studies in Europe. In addition, it appears that results from the small plot studies may be adequately scaled up one to two orders of magnitude in terms of land areas considered. Copyright © 2000 John Wiley & Sons, Ltd.     

Empirical and conceptual modelling of the suspended sediment dynamics in a large tropical African river: the Upper Niger river basin

A 7-year sediment transport monitoring on the Upper Niger rivers was used to study the relationship between suspended sediment concentration and river discharge. During annual floods, these relationships show positive hysteresis. This paper presents the results of two models that estimate the time evolution of suspended sediment concentration using water discharge data only. The first model is based on a statistical approach using two relationships, one for the rising stage period of the flood and one for the recession period of the annual flood; the second model is a lumped conceptual one; it supposes that the sediment flux observed in the river comes from two different sources of sediment and that these two sources may be regarded as two different reservoirs. The erosion of the first reservoir represents hillslope erosion observed during the runoff season. Sediment supply from this ‘reservoir’ is limited in time because depletion occurs during the runoff season. The second reservoir is unlimited in time and quantity and its erosion represents contributions coming from bank erosion and mobilisation of deposits in the channel network.

Both of the models are compared with a simple rating curve based model. The model results show that the conceptual model has the highest efficiency to reproduce from weekly discharge only the time evolution of weekly suspended sediment concentrations, the time evolution of weekly sediment fluxes, and the global annual sediment yields.