Catchment-scale modelling of flow and nutrient transport in the Chalk unsaturated zone

The unsaturated zone exerts a major control on the delivery of nutrients to Chalk streams, yet flow and transport processes in this complex, dual-porosity medium have remained controversial. A major challenge arises in characterising these processes, both at the detailed mechanistic level and at an appropriate level for inclusion within catchment-scale models for nutrient management. The lowland catchment research (LOCAR) programme in the UK has provided a unique set of comprehensively instrumented groundwater-dominated catchments. Of these, the Pang and Lambourn, tributaries of the Thames near Reading, have been a particular focus for research into subsurface processes and surface water-groundwater interactions. Data from LOCAR and other sources, along with a new dual permeability numerical model of the Chalk, have been used to explore the relative roles of matrix and fracture flow within the unsaturated zone and resolve conflicting hypotheses of response. From the improved understanding gained through these explorations, a parsimonious conceptualisation of the general response of flow and transport within the Chalk unsaturated zone was formulated. This paper summarises the modelling and data findings of these explorations, and describes the integration of the new simplified unsaturated zone representation with a catchment-scale model of nutrients (INCA), resulting in a new model for catchment-scale flow and transport within Chalk systems: INCA-Chalk. This model is applied to the Lambourn, and results, including hindcast and forecast simulations, are presented. These clearly illustrate the decadal time-scales that need to be considered in the context of nutrient management and the EU Water Framework Directive. (C) 2007 Elsevier B.V. All rights reserved.

Fine sediment delivery and transfer in lowland catchments: modelling suspended sediment concentrations in response to hydrological forcing

Fine sediment delivery to and storage in stream channel reaches can disrupt aquatic habitats, impact river hydromorphology, and transfer adsorbed nutrients and pollutants from catchment slopes to the fluvial system. This paper presents a modelling toot for simulating the time-dependent response of the fine sediment system in catchments, using an integrated approach that incorporates both land phase and in-stream processes of sediment generation, storage and transfer. The performance of the model is demonstrated by applying it to simulate in-stream suspended sediment concentrations in two lowland catchments in southern England, the Enborne and the Lambourn, which exhibit contrasting hydrological and sediment responses due to differences in substrate permeability. The sediment model performs well in the Enborne catchment, where direct runoff events are frequent and peak suspended sediment concentrations can exceed 600 mg l(-1). The general trends in the in-stream concentrations in the Lambourn catchment are also reproduced by the model, although the observed concentrations are low (rarely exceeding 50 mg l(-1)) and the background variability in the concentrations is not fully characterized by the model. Direct runoff events are rare in this highly permeable catchment, resulting in a weak coupling between the sediment delivery system and the catchment hydrology. The generic performance of the model is also assessed using a generalized sensitivity analysis based on the parameter bounds identified in the catchment applications. Results indicate that the hydrological parameters contributing to the sediment response include those controlling (1) the partitioning of runoff between surface and soil zone flows and (2) the fractional loss of direct runoff volume prior to channel delivery. The principal sediment processes controlling model behaviour in the simulations are the transport capacity of direct runoff and the in-stream generation, storage and release of the fine sediment fraction. The in-stream processes appear to be important in maintaining the suspended sediment concentrations during low flows in the River Enborne and throughout much of the year in the River Lambourn. Copyright (c) 2007 John Wiley & Sons, Ltd.

Modelling sediment supply and transport in the River Lugg: strategies for controlling sediment loads

The River Lugg has particular problems with high sediment loads that have resulted in detrimental impacts on ecology and fisheries. A new dynamic, process-based model of hydrology and sediments (INCA- SED) has been developed and applied to the River Lugg system using an extensive data set from 1995–2008. The model simulates sediment sources and sinks throughout the catchment and gives a good representation of the sediment response at 22 reaches along the River Lugg. A key question considered in using the model is the management of sediment sources so that concentrations and bed loads can be reduced in the river system. Altogether, five sediment management scenarios were selected for testing on the River Lugg, including land use change, contour tillage, hedging and buffer strips. Running the model with parameters altered to simulate these five scenarios produced some interesting results. All scenarios achieved some reduction in sediment levels, with the 40% land use change achieving the best result with a 19% reduction. The other scenarios also achieved significant reductions of between 7% and 9%. Buffer strips produce the best result at close to 9%. The results suggest that if hedge introduction, contour tillage and buffer strips were all applied, sediment reductions would total 24%, considerably improving the current sediment situation. We present a novel cost-effectiveness analysis of our results where we use percentage of land removed from production as our cost function. Given the minimal loss of land associated with contour tillage, hedges and buffer strips, we suggest that these management practices are the most cost-effective combination to reduce sediment loads.

Transport in the low-latitude tropopause zone diagnosed using particle trajectories

Recent literature has described a “transition zone” between the average top of deep convection in the Tropics and the stratosphere. Here transport across this zone is investigated using an offline trajectory model. Particles were advected by the resolved winds from the European Centre for Medium-Range Weather Forecasts reanalyses. For each boreal winter clusters of particles were released in the upper troposphere over the four main regions of tropical deep convection (Indonesia, central Pacific, South America, and Africa). Most particles remain in the troposphere, descending on average for every cluster. The horizontal components of 5-day trajectories are strongly influenced by the El Niño–Southern Oscillation (ENSO), but the Lagrangian average descent does not have a clear ENSO signature. Tropopause crossing locations are first identified by recording events when trajectories from the same release regions cross the World Meteorological Organization lapse rate tropopause. Most crossing events occur 5–15 days after release, and 30-day trajectories are sufficiently long to estimate crossing number densities. In a further two experiments slight excursions across the lapse rate tropopause are differentiated from the drift deeper into the stratosphere by defining the “tropopause zone” as a layer bounded by the average potential temperature of the lapse rate tropopause and the profile temperature minimum. Transport upward across this zone is studied using forward trajectories released from the lower bound and back trajectories arriving at the upper bound. Histograms of particle potential temperature (θ) show marked differences between the transition zone, where there is a slow spread in θ values about a peak that shifts slowly upward, and the troposphere below 350 K. There forward trajectories experience slow radiative cooling interspersed with bursts of convective heating resulting in a well-mixed distribution. In contrast θ histograms for back trajectories arriving in the stratosphere have two distinct peaks just above 300 and 350 K, indicating the sharp change from rapid convective heating in the well-mixed troposphere to slow ascent in the transition zone. Although trajectories slowly cross the tropopause zone throughout the Tropics, all three experiments show that most trajectories reaching the stratosphere from the lower troposphere within 30 days do so over the west Pacific warm pool. This preferred location moves about 30°–50° farther east in an El Niño year (1982/83) and about 30° farther west in a La Niña year (1988/89). These results could have important implications for upper-troposphere–lower-stratosphere pollution and chemistry studies.

Wave refraction and longshore transport patterns along the southern Santa Catarina coast

Com base no clima de ondas para o sul da costa brasileira, a refração de ondas foi modelada para a obtenção de alturas e direção das ondas ao longo da costa sul de Santa Catarina. Essas informações permitem a estimativa da deriva litorânea potencial para a região. De acordo com a orientação da linha de costa, diferentes setores do litoral apresentam variados padrões de deriva. As estimativas foram realizadas com base nos dados médios anuais e para cada estação do ano, demonstrando assim a variabilidade dos padrões ao longo do ano. A configuração da linha de costa e os resultados das estimativas de intensidade e direção da deriva indicam o predomínio da deriva litorânea na porção sul da área e um equilíbrio entre processos de espraiamento e de deriva na porção norte. Padrões contrastantes de deriva entre a porção sul e norte da linha de costa indicam um balanço positivo de sedimentos na porção central, disponibilizando sedimentos para os processos de transporte através da costa. Considerando escalas temporais maiores, os padrões de deriva concordam com o processo de preenchimento costeiro que é controlado principalmente pelo aporte de sedimentos pela deriva litorânea.

Hydrogeology, Chemical Characteristics, and Transport Processes in the Zone of Contribution of a Public-Supply Well in York, Nebraska

In 2001, the U.S. Geological Survey, as part of the National Water Quality Assessment (NAWQA) Program, initiated a topical study of Transport of Anthropogenic and Natural Contaminants (TANC) to PSW (public-supply wells). Local-scale and regional-scale TANC study areas were delineated within selected NAWQA study units for intensive study of processes effecting transport of contaminants to PSWs. This report describes results from a local-scale TANC study area at York, Nebraska, within the High Plains aquifer, including the hydrogeology and geochemistry of a 108-square-kilometer study area that contains the zone of contribution to a PSW selected for study (study PSW), and describes factors controlling the transport of selected anthropogenic and natural contaminants to PSWs.

Zum Transport von partikelgebundenen Schadstoffen in der ungesättigten Zone

Bei einer Risikoabschätzung bezüglich einer Gefährdung des Schutzgutes Grundwasser müssen alle relevanten Transportpfade, auf denen Schadstoffe durch die Bodenzone bis ins Grundwasser verlagert werden, identifiziert und quantifiziert werden. Die Verlagerung von Schadstoffen gebunden an mobile Partikel im Sickerwasser wird dabei oft vernachlässigt. In dieser Arbeit wurden sowohl experimentelle Untersuchungen zum Partikeltransport in der Bodenzone als auch Szenarienmodellierungen hinsichtlich der Wechselwirkung Partikel/Schadstoff durchgeführt. Die experimentellen ungesättigten Säulenversuche wurden unter naturnahen stationären und instationären hydraulischen und hydrochemischen Bedingungen durchgeführt. Dabei wurde der Einfluss der Parameter Durchmesser Bodenmatrix, Partikelgröße, Beregnungsintensität, Oberflächenspannung und Hydrochemie auf den Transport von natürlichen und synthetischen Partikeln untersucht. Des Weiteren wurden Untersuchungen zur partikelgebundenen Verlagerung von Phenanthren durchgeführt. In einer numerischen Szenarienmodellierung mit dem Modell SMART wurde untersucht, unter welchen Randbedingungen der Transport von Partikeln gleichzeitig zu signifikanten partikelgebundenen Schadstoffkonzentrationen im Grundwasser führt. Dabei wurden die Parameter Lithologie Partikel/Boden, Hydrophobizität Schadstoff, Partikelkonzentration, Partikeldurchmesser sowie Körnung Bodenmatrix variiert. Die Ergebnisse dieser Arbeit zeigen, dass der partikelgebundene Schadstofftransportpfad in der ungesättigten Bodenzone in verschiedenen Szenarien den Anteil mobiler Schadstoffe, die mit dem Sickerwasser ins Grundwasser gelangen, signifikant erhöht. Auf Basis der experimentellen und theoretischen Untersuchungen wurde ein zweistufiges Bewertungsschema entwickelt, das bereits im Vorfeld einer Risikoabschätzung als Entscheidungshilfe hinsichtlich der Relevanz einer Mobilisierung, eines Transports und des Rückhalts von partikelgebundenen Schadstoffen in der ungesättigten Zone dient.

Reactive mass transport modelling of a three-dimensional vertical fault zone with a finger-like convective flow regime

For a three-dimensional vertically-oriented fault zone, we consider the coupled effects of fluid flow, heat transfer and reactive mass transport, to investigate the patterns of fluid flow, temperature distribution, mineral alteration and chemically induced porosity changes. We show, analytically and numerically, that finger-like convection patterns can arise in a vertically-oriented fault zone. The onset and patterns of convective fluid flow are controlled by the Rayleigh number which is a function of the thermal properties of the fluid and the rock, the vertical temperature gradient, and the height and the permeability of the fault zone. Vigorous fluid flow causes low temperature gradients over a large region of the fault zone. In such a case, flow across lithological interfaces becomes the most important mechanism for the formation of sharp chemical reaction fronts. The degree of rock buffering, the extent and intensity of alteration, the alteration mineralogy and in some cases the formation of ore deposits are controlled by the magnitude of the flow velocity across these compositional interfaces in the rock. This indicates that alteration patterns along compositional boundaries in the rock may provide some insights into the convection pattern. The advective mass and heat exchanges between the fault zone and the wallrock depend on the permeability contrast between the fault zone and the wallrock. A high permeability contrast promotes focussed convective flow within the fault zone and diffusive exchange of heat and chemical reactants between the fault zone and the wallrock. However, a more gradual permeability change may lead to a regional-scale convective flow system where the flow pattern in the fault affects large-scale fluid flow, mass transport and chemical alteration in the wallrocks

Hazard index maps for woody material recruitment and transport in alpine catchments

A robust and reliable risk assessment procedure for hydrologic hazards deserves particular attention to the role of transported woody material during flash floods or debris flows. At present, woody material transport phenomena are not systematically considered within the procedures for the elaboration of hazard maps. The consequence is a risk of losing prediction accuracy and of underestimating hazard impacts. Transported woody material frequently interferes with the sediment regulation capacity of open check dams and moreover, when obstruction phenomena at critical crosssections of the stream occur, inundations can be triggered. The paper presents a procedure for the determination of the relative propensity of mountain streams to the entrainment and delivery of recruited woody material on the basis of empirical indicators. The procedure provided the basis for the elaboration of a hazard index map for all torrent catchments of the Autonomous Province of Bolzano/Bozen. The plausibility of the results has been thoroughly checked by a backward oriented analysis on natural hazard events, documented since 1998 at the Department of Hydraulic Engineering of the aforementioned Alpine Province. The procedure provides hints for the consideration of the effects, induced by woody material transport, during the elaboration of hazard zone maps.

Fecal pellets flux at DYFAMED sediment traps

Because zooplankton feces represent a potentially important transport pathway of surface-derived organic carbon in the ocean, we must understand the patterns of fecal pellet abundance and carbon mobilization over a variety of spatial and temporal scales. To assess depth-specific water column variations of fecal pellets on a seasonal scale, vertical fluxes of zooplankton fecal pellets were quantified and their contribution to mass and particulate carbon were computed during 1990 at 200, 500, 1000, and 2000 m depths in the open northwestern Mediterranean Sea as part of the French-JGOFS DYFAMED Program. Depth-averaged daily fecal pellet flux was temporally variable, ranging from 3.04 * 10**4 pellets m**2/d in May to a low of 6.98 * 10**2 pellets m**2/d in September. The peak flux accounted for 50% of the integrated annual flux of fecal pellets and 62% of pellet carbon during only two months in mid-spring (April and May). Highest numerical fluxes were encountered at 1000 m, suggesting fecal pellet generation well below the euphotic zone. However, there was a trend toward lower pellet carbon with increasing depth, suggesting bacterial degradation or in situ repackaging as pellets sink through the water column. At 500 m, both the lowest pellet numerical abundance and carbon flux were evident during the spring peak. Combined with data indicating that numerical and carbon fluxes are dominated at 500 m by a distinct type of pellet found uniquely at this depth, these trends suggest the presence of an undescribed mid-water macro-zooplankton or micro-nekton community. Fecal pellet carbon flux was highest at 200 m and varied with depth independently of overall particulate carbon, which was greatest at 500 m. Morphologically distinct types of pellets dominated the numerical and carbon fluxes. Small elliptical and spherical pellets accounted for 88% of the numerical flux, while larger cylindrical pellets, although relatively rare (<10%), accounted for almost 40% of the overall pellet carbon flux. Cylindrical pellets dominated the pellet carbon flux at all depths except 500 m, where a large subtype of elliptical pellet, found only at that depth, was responsible for the majority of pellet carbon flux. Overall during 1990, fecal pellets were responsible for a depth-integrated annual average flux of 1.03 mgC/m**2/d, representing 18% of the total carbon flux. The proportion of vertical carbon flux attributed to fecal pellets varied from 3 to 35%, with higher values occurring during periods when the water column was vertically mixed. Especially during these times, fecal pellets are a critical conveyor of carbon to the deep sea in this region.

Sediment metabolism in samples off Northwest Africa

Sediments off northwest Africa were assayed for activities of the respiratory electron transport system (ETS) and for primary amino nitrogen. ETS activities were used to compute respiratory oxygen consuption, carbon oxidation, and nitrate reduction rates. Activities were correlated with depth of the water column, and their longshore distribution resembled that of euphotic zone phytoplankton productivity. Protein concentrations were closely correlated with ETS activities. Carbon biomass was calculated from protein and compared with other computed values. The carbon oxidation rate accounted for 13 % of the region's primary production.

Multi-proxy analyses of sediment cores from the Pakistan continental margin

Late Holocene laminated sediments from a core transect centred in the oxygen minimum zone (OMZ) impinging at the continental slope off Pakistan indicate stable oxygen minimum conditions for the past 7000 calendar years. High SW-monsoon-controlled biological productivity and enhanced organic matter preservation during this period is reflected in high contents of total organic carbon (TOC) and redox-sensitive elements (Ni, V), as well as by a low-diversity, high-abundance benthic foraminiferal Buliminacea association and high abundance of the planktonic species Globigerina bulloides indicative of upwelling conditions. Surface-water productivity was strongest during SW monsoon maxima. Stable OMZ conditions (reflected by laminated sediments) were found also during warm interstadial events (Preboreal, Bølling-Allerød, and Dansgaard-Oeschger events), as well as during peak glacial times (17-22.5 ka, all ages in calendar years). Sediment mass accumulation rates were at a maximum during the Preboreal and Younger Dryas periods due to strong riverine input and mobilisation of fine-grained sediment coinciding with rapid deglacial sea-level rise, whereas eolian input generally decreased from glacial to interglacial times. In contrast, the occurrence of bioturbated intervals from 7 to 10.5 ka (early Holocene), in the Younger Dryas (11.7-13 ka), from 15 to 17 ka (Heinrich event 1) and from 22.5 to 25 ka (Heinrich event 2) suggests completely different conditions of oxygen-rich bottom waters, extremely low mass and organic carbon accumulation rates, a high-diversity benthic fauna, all indicating lowered surface-water productivity. During these intervals the OMZ was very poorly developed or absent and a sharp fall of the aragonite compensation depth favoured the preservation of pteropods. The abundance of lithogenic proxies suggests aridity and wind transport by northwesterly or northeasterly winds during these periods coinciding with the North Atlantic Heinrich events and dust peaks in the Tibetan Loess records. The correlation of the monsoon-driven OMZ variability in the Arabian Sea with the rapid climatic fluctuations in the high northern latitudes suggests a close coupling between the climates of the high and low latitudes at a global scale.

Grain entrainment and transport rates of magnetic minerals of coastal sediment samples from the Bay of Plenty, New Zealand

Heavy (magnetic & non-magnetic) minerals are found concentrated by natural processes in many fluvial, estuarine, coastal and shelf environments with a potential to form economic placer deposits. Understanding the processes of heavy mineral transport and enrichment is prerequisite to interpret sediment magnetic properties in terms of hydro- and sediment dynamics. In this study, we combine rock magnetic and sedimentological laboratory measurements with numerical 3D discrete element models to investigate differential grain entrainment and transport rates of magnetic minerals in a range of coastal environments (riverbed, mouth, estuary, beach and near-shore). We analyzed grain-size distributions of representative bulk samples and their magnetic mineral fractions to relate grain-size modes to respective transport modes (traction, saltation, suspension). Rock magnetic measurements showed that distribution shapes, population sizes and grain-size offsets of bulk and magnetic mineral fractions hold information on the transport conditions and enrichment process in each depositional environment. A downstream decrease in magnetite grain size and an increase in magnetite concentration was observed from riverine source to marine sink environments. Lower flow velocities permit differential settling of light and heavy mineral grains creating heavy mineral enriched zones in estuary settings, while lighter minerals are washed out further into the sea. Numerical model results showed that higher heavy mineral concentrations in the bed increased the erosion rate and enhancing heavy mineral enrichment. In beach environments where sediments contained light and heavy mineral grains of equivalent grain sizes, the bed was found to be more stable with negligible amount of erosion compared to other bed compositions. Heavy mineral transport rates calculated for four different bed compositions showed that increasing heavy mineral content in the bed decreased the transport rate. There is always a lag in transport between light and heavy minerals which increases with higher heavy mineral concentration in all tested bed compositions. The results of laboratory experiments were validated by numerical models and showed good agreement. We demonstrate that the presented approach bears the potential to investigate heavy mineral enrichment processes in a wide range of sedimentary settings.

A time-slice study in Heinrich Stadial 1 and Mid Holocene sediment cores off northern Senegal

Sediment dynamics in limnic, fluvial and marine environments can be assessed by granulometric and rock-magnetic methodologies. While classical grain-size analysis by sieving or settling mainly bears information on composition and transport, the magnetic mineral assemblages reflect to a larger extent the petrology and weathering conditions in the sediment source areas. Here, we combine both methods to investigate Late Quaternary marine sediments from five cores along a transect across the continental slope off Senegal. This region near the modern summer Intertropical Convergence Zone is particularly sensitive to climate change and receives sediments from several aeolian, fluvial and marine sources. From each of the investigated five GeoB sediment cores (494-2956 m water depth) two time slices were processed which represent contrasting climatic conditions: the arid Heinrich Stadial 1 (~ 15 kyr BP) and the humid Mid Holocene (~ 6 kyr BP). Each sediment sample was split into 16 grain-size fractions ranging from 1.6 to 500 µm. Concentration and grain-size indicative magnetic parameters (susceptibility, SIRM, HIRM, ARM and ARM/IRM) were determined at room temperature for each of these fractions. The joint consideration of whole sediment and magnetic mineral grain-size distributions allows to address several important issues: (i) distinction of two aeolian sediment fractions, one carried by the north-easterly trade winds (40-63 µm) and the other by the overlying easterly Harmattan wind (10-20 µm) as well as a fluvial fraction assigned to the Senegal River (< 10 µm); (ii) identification of three terrigenous sediment source areas: southern Sahara and Sahel dust (low fine-grained magnetite amounts and a comparatively high haematite content), dust from Senegalese coastal dunes (intermediate fine-grained magnetite and haematite contents) and soils from the upper reaches of the Senegal River (high fine-grained magnetite content); (iii) detection of partial diagenetic dissolution of fine magnetite particles as a function of organic input and shore distance; (iv) analysis of magnetic properties of marine carbonates dominating the grain-size fractions 63-500 µm.