Validation and Application of the Second Generation Three Dimensional Hydrodynamic Model of Chesapeake Bay

The validation and subsequent application of the current three-dimensional numerical hydrodynamic model of Chesapeake Bay is presented. The numerical model solves conservation equations for water mass, momentum, salinity, and heat on a boundary-fitted grid in the horizontal plane and a Cartesian z-grid in the vertical. A generalized ADI finite difference scheme is employed in conjunction with mode splitting technique, solving external and the internal modes. The 10-year boundary conditions including tide, slinity, temperature, wind, heat exchange coefficient, river and non-point source flows were constructed. Model validation was accomplished by demonstrating the model's ability to reproduce observed data over time scales ranging from tidal to seasonal periods. The major parameters compared include tidal elevation, intra-tidal and residual velocities, salinity, temperature, stratification, and flux calculated through the Bay mouth.After validation, the model was applied to simulate bay hydrodynamics for the 10 years of 1985–94. These results were used to drive the three-dimensional water quality model of Chesapeake Bay, which is discussed in a companion paper.     

The impact of pro‐glacial lakes on downstream sediment delivery in Norway

Glacier recessions caused by climate change may uncover pro‐glacial lakes that form important sedimentation basins regulating the downstream sediment delivery. The impact of modern pro‐glacial lakes on fluvial sediment transport from three different Norwegian glaciers: Nigardsbreen, Engabreen and Tunsbergdalsbreen, and their long‐term development has been studied. All of these lakes developed in modern times in overdeepened bedrock basins. The recession of Nigardsbreen uncovered a 1.8 km long and on average 15 m deep pro‐glacial lake basin during 1937 to 1968. Since then the glacier front has been situated entirely on land, and the sediment input and output of the lake has been measured. The suspended sediment transport into and out of the lake averaged 11 730 t yr^?1 and 2340 t yr^?1 respectively. Thus, 20% remained in suspension at the outlet. The measured mean annual bedload supplied to the lake was 11 800 t yr^?1, giving a total transport of 23 530 t yr^?1 which corresponds to a specific sediment yield of 561 t km^?2 yr^?1. A 1.9 km long and up to 90 m deep pro‐glacial lake basin downstream from Engabreen glacier was uncovered during 1890 to 1944. The average suspended sediment load delivered from the glacier during the years 1970–1981 amounted to 12 375 t yr^?1and the transport out of the lake was 2021 t yr^?1, giving an average of 16% remaining in suspension. The mean annual bedload was 8000 t yr^?1, thus the total transport was 20 375 t yr^?1, giving a specific sediment yield of 566 t km^?2 yr^?1. For Tunsbergdalsbreen glacier, measurements in the early 1970s indicated that the suspended sediment transport was on average 44 000 t yr^?1. From 1987 to 1993 the recession of the glacier uncovered a small pro‐glacial lake, 0.3 km long and around 9 m deep. Downstream from this, the suspended sediment load measured in 2009 was 28 000 t yr^?1, indicating that as much as 64% remained in suspension. Flow velocity, grain size of sediment, and morphology of the lake are important factors controlling the sedimentation rate in the pro‐glacial lakes. A survey of the sub‐glacial morphology of Tunsbergdalsbreen revealed that there are several overdeepened basins beneath the glacier. The largest is 4 km long and 100 m deep. When the glacier melts back they will become lakes and act as sedimentation basins. Despite an expected increase in sediment yield from the glacier, little sediment will pass these lakes and downstream sediment delivery will be reduced markedly. Beneath Nigardsbreen there was only a small depression that may form a lake and the sediment delivery will not be significantly affected. © 2014 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

The responses of hydrological processes and sediment yield to land‐use and climate change in the Be River Catchment,Vietnam

In this study, we investigated the responses of hydrology and sediment yield with impacts of land‐use and climate change scenarios in the Be River Catchment, using the Soil and Water Assessment Tool (SWAT) hydrological model. The calibration and validation results indicated that the SWAT model is a powerful tool for simulating the impact of environmental change on hydrology and sediment yield in this catchment. The hydrologic and sediment yield responses to land‐use and climate changes were simulated based on the calibrated model. The results indicated that a 16.3% decrease in forest land is likely to increase streamflow (0.2 to 0.4%), sediment load (1.8 to 3.0%), and surface runoff (SURQ) (4.8 to 10.7%) and to decrease groundwater discharge (GW_Q) (3.5 to 7.9%). Climate change in the catchment leads to decreases in streamflow (0.7 to 6.9%) and GW_Q (3.0 to 8.4%), increase in evapotranspiration (0.5 to 2.9%), and changes in SURQ (?5.3 to 2.3%) and sediment load (?5.3 to 4.4%). The combined impacts of land‐use and climate changes decrease streamflow (2.0 to 3.9%) and GW_Q (12.3 to 14.0%), increase evapotranspiration (0.7 to 2.8%), SURQ (8.2 to 12.4%), and sediment load (2.0 to 7.9%). In general, the separate impacts of climate and land‐use changes on streamflow, sediment load, and water balance components are offset each other. However, SURQ and some component of subsurface flow are more sensitive to land‐use change than to climate change. Furthermore, the results emphasized water scarcity during the dry season and increased soil erosion during the wet season. Copyright © 2012 John Wiley & Sons, Ltd.     

Wind effects on the lateral structure of density-driven circulation in Chesapeake Bay

The response of the density-driven circulation in the Chesapeake Bay to wind forcing was studied with numerical experiments. A model of the bay with realistic bathymetry was first applied to produce the density-driven flow under average river discharge and tidal forcing. Subsequently, four spatially uniform wind fields (northeasterly, northwesterly, southwesterly, and southeasterly) were imposed to examine the resulting cross-estuary structure of salinity and flow fields. In general, northeasterly and northwesterly winds intensified the density-driven circulation in the upper and middle reaches of the bay, whereas southeasterly and southwesterly winds weakened it. The response was different in the lower bay, where downwind flow from the upper and middle reaches of the bay competed with onshore/offshore coastal flows. Wind remote effects were dominant, over local effects, on volume transports through the bay entrance. However, local effects were more influential in establishing the sea-level slopes that drove subtidal flows and salinity fields in most of the bay. The effect of vertical stratification on wind-induced flows was also investigated by switching it off. The absence of stratification allowed development of Ekman layers that reached depths of the same order as the water depth. Consequently, bathymetric effects became influential on the homogeneous flow structure causing the wind-induced flow inside the bay to show a marked transverse structure: downwind over the shallow areas and upwind in the channels. In the presence of stratification, Ekman layers became shallower and the wind-induced currents showed weaker transverse structure than those that developed in the absence of stratification. In essence, the wind-driven flows were horizontally sheared under weak stratification and vertically sheared under stratified conditions.     

Lateral dynamics and associated transport of sediment in the upper reaches of a partially mixed estuary,Chesapeake Bay,USA

Data from time series of transects made over a tidal period across a section of the upper Chesapeake Bay, USA, reveal the influence of lateral dynamics on sediment transport in an area with a deep channel and broad extents of shallower flanks. Contributions to lateral momentum by rotation (Coriolis plus channel curvature), cross channel density gradients and cross channel surface slope were estimated, and the friction and acceleration terms needed to complete the balance were compared to patterns of observed lateral circulation. During ebb, net rotation effects were larger because of river velocity and reinforcement of Coriolis by curvature. During flood, stratification was greater because of landward advection of strong vertical density gradients. Together, the ebb intensified lateral circulation and flood intensified stratification focused sediment and sediment transport along the left side of the estuary (looking seaward). The tendency for greater stratification on flood and net sediment flux toward the left-hand shoal are contrary to more common models which, in the northern hemisphere, predict greater resuspension on flood and move sediment toward the right-hand shoal. These tidal asymmetries interact with the lateral circulation to focus net sediment flux on the left side of the estuary, and to produce net ebb directed sediment transport at the surface of the same order of magnitude as net flood directed sediment transport at the bottom.     

Computing spatially distributed sediment delivery ratios: inferring functional sediment connectivity from repeat high‐resolution digital elevation models

High‐resolution digital elevation models (DEMs) from repeat LiDAR (light detection and ranging) or SfM (structure from motion) surveys have become an important tool in process geomorphology. The spatial pattern of negative and positive changes of surface elevation on raster DEMs of difference (DoD) can be interpreted in terms of geomorphic processes, and has been used for morphological budgeting. We show how the application of flow routing algorithms and flow accumulation opens new opportunities for the analysis of DoD. By accumulating the values of the DoD along downslope flowpaths delineated on a DEM, these algorithms lend themselves to computing the net balance, i.e. sediment yield (SY), for the contributing area of each cell. Doing the same for the negative subset of the DoD yields a minimum estimate of erosion (E) within the contributing area. The division of SY by E yields (a maximum estimate of) the sediment delivery ratio (SDR), that is the proportion of material eroded within the contributing area of each cell that has been exported from that area. The resulting SDR is a spatially distributed measure of functional sediment connectivity. In this letter, we develop the computationally simple approach by means of an example DoD from a lateral moraine section in the Upper Kaunertal Valley, Austrian Central Alps. We also discuss advantages, assumptions and limitations, and outline potential applications to connectivity research using field‐, laboratory‐, and model‐based DoD. Copyright © 2018 John Wiley & Sons, Ltd.     

Evaluation of a dynamic multi‐class sediment transport model in a catchment under soil‐conservation agriculture

Soil erosion models are essential tools for the successful implementation of effective and adapted soil conservation measures on agricultural land. Therefore, models are needed that predict sediment delivery and quality, give a good spatial representation of erosion and deposition and allow us to account for various soil conservation measures. Here, we evaluate how well a modified version of the spatially distributed multi‐class sediment transport model (MCST) simulates the effectiveness of control measures for different event sizes. We use 8 year runoff and sediment delivery data from two small agricultural watersheds (0·7 and 3·7 ha) under optimized soil conservation. The modified MCST model successfully simulates surface runoff and sediment delivery from both watersheds; one of which was dominated by sheet and the other was partly affected by rill erosion. Moreover, first results of modelling enrichment of clay in sediment delivery are promising, showing the potential of MCST to model sediment enrichment and nutrient transport. In general, our results and those of an earlier modelling exercise in the Belgian Loess Belt indicate the potential of the MCST model to evaluate soil erosion and deposition under different agricultural land uses. As the model explicitly takes into account the dominant effects of soil‐conservation agriculture, it should be successfully applicable for soil‐conservation planning/evaluation in other environments. Copyright © 2008 John Wiley & Sons, Ltd.     

Modelling the hydrologic effects of dynamic land‐use change using a distributed hydrologic model and a spatial land‐use allocation model

This study develops a novel approach for modelling and examining the impacts of time–space land‐use changes on hydrological components. The approach uses an empirical land‐use change allocation model (CLUE‐s) and a distributed hydrological model (DHSVM) to examine various land‐use change scenarios in the Wu‐Tu watershed in northern Taiwan. The study also uses a generalized likelihood uncertainty estimation approach to quantify the parameter uncertainty of the distributed hydrological model. The results indicate that various land‐use policies—such as no change, dynamic change and simultaneous change—have different levels of impact on simulating the spatial distributions of hydrological components in the watershed study. Peak flow rates under simultaneous and dynamic land‐use changes are 5·71% and 2·77%, respectively, greater than the rate under the no land‐use change scenario. Using dynamic land‐use changes to assess the effect of land‐use changes on hydrological components is more practical and feasible than using simultaneous land‐use change and no land‐use change scenarios. Furthermore, land‐use change is a spatial dynamic process that can lead to significant changes in the distributions of ground water and soil moisture. The spatial distributions of land‐use changes influence hydrological processes, such as the ground water level of whole areas, particularly in the downstream watershed. Copyright © 2010 John Wiley & Sons, Ltd.     

Developing a new stream metric for comparing stream function using a bank–floodplain sediment budget: a case study of three Piedmont streams

A bank and floodplain sediment budget was created for three Piedmont streams tributary to the Chesapeake Bay. The watersheds of each stream varied in land use from urban (Difficult Run) to urbanizing (Little Conestoga Creek) to agricultural (Linganore Creek). The purpose of the study was to determine the relation between geomorphic parameters and sediment dynamics and to develop a floodplain trapping metric for comparing streams with variable characteristics. Net site sediment budgets were best explained by gradient at Difficult Run, floodplain width at Little Conestoga Creek, and the relation of channel cross‐sectional area to floodplain width at Linganore Creek. A correlation for all streams indicated that net site sediment budget was best explained by relative floodplain width (ratio of channel width to floodplain width). A new geomorphic metric, the floodplain trapping factor, was used to compare sediment budgets between streams with differing suspended sediment yields. Site sediment budgets were normalized by floodplain area and divided by the stream's sediment yield to provide a unitless measure of floodplain sediment trapping. A floodplain trapping factor represents the amount of upland sediment that a particular floodplain site can trap (e.g. a factor of 5 would indicate that a particular floodplain site traps the equivalent of 5 times that area in upland erosional source area). Using this factor we determined that Linganore Creek had the highest gross and net (floodplain deposition minus bank erosion) floodplain trapping factor (107 and 46, respectively) that Difficult Run the lowest gross floodplain trapping factor (29) and Little Conestoga Creek had the lowest net floodplain trapping factor (–14, indicating that study sites were net contributors to the suspended sediment load). The trapping factor is a robust metric for comparing three streams of varied watershed and geomorphic character, it promises to be a useful tool for future stream assessments. Published 2012. This article is a U.S. Government work and is in the public domain in the USA.     

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.     

Gully catchments as a sediment sink,not just a source: Results from a long‐term (~12 500 year) sediment budget

Sediment delivery from hillslopes to trunk streams represents a significant pathway of mass transfer in the landscape, with a large fraction facilitated by gully systems. The internal gully geomorphic dynamics represent a considerable gap in many landscape and empirical erosion models, therefore a better understanding of these processes over longer timescales (10–10⁴ years) is needed. This study analyses the sediment mass balance and storage dynamics within a headwater gully catchment in central Europe over the last ~12 500 years. Human induced erosion resulted in hillslope erosion rates ~2.3 times higher than under naturally de‐vegetated conditions (during the Younger Dryas), however the total sediment inputs to the gully system (and therefore gully aggradation), were similar. Net gully storage has consistently increased to become the second largest term in the sediment budget after hillslope erosion (storage is ~45% and ~73% of inputs during two separate erosion and aggradation cycles). In terms of the depletion of gully sediment storage, the sediment mass balance shows that export beyond the gully fan was not significant until the last ~500 years, due to reduced gully fan accommodation space. The significance of storage effects on the gully sediment mass balance, particularly the export terms, means that it would be difficult to determine the influences of human impact and/or climatic changes from floodplain or lake sedimentary archives alone and that the sediment budgets of the headwater catchments from which they drain are more likely to provide these mechanistic links. Copyright © 2015 John Wiley & Sons, Ltd.     

Estimation of suspended sediment loads and delivery in an incised upland headwater catchment,south‐eastern Australia

Historically upland headwater catchments in south‐eastern Australia have undergone extensive gully erosion that has removed large amounts of sediment to lowlands. Recent research suggests these upland areas may continue to dominate fine sediment loads in lowland rivers. Improved understanding of sediment transfer through upland headwater catchments may have implications for interpreting downstream sediment supply. In this study a nested catchment design was utilized to examine suspended sediment yields and delivery from a small tributary sub‐catchment (1·64 km²) to the study catchment outlet (53·5 km²). Monitoring of suspended sediment concentration and discharge was undertaken for a period of nearly two years and used to estimate suspended sediment loads. Estimated total suspended sediment exports over the period of monitoring were 24·16 t from the sub‐catchment and 550·3 t from the catchment, which are generally less than previous reported small catchment yields in south‐eastern Australia. The extent of sediment delivery was examined using between‐site ratios of specific sediment yield per unit area and incised channel length. Sediment delivery was high under average rainfall conditions, but seasonally dependent. Both suspended sediment yields and the extent of delivery peaked over spring months, supplemented by remobilization of sediment stored during summer months in the main catchment channel. The findings of this study suggest much of the suspended sediment exported from small incised upland sub‐catchments (1–2 km²) may be delivered to downstream reaches under average rainfall conditions, which, in conjunction with the findings of previous research supports the potential importance of contributions from these areas to suspended sediment loads in lowland rivers during high flow periods. Copyright © 2007 John Wiley & Sons, Ltd.     

Interactions between sediment delivery,channel change,climate change and flood risk in a temperate upland environment

This paper uses numerical simulation of flood inundation based on a coupled one‐dimensional–two‐dimensional treatment to explore the impacts upon flood extent of both long‐term climate changes, predicted to the 2050s and 2080s, and short‐term river channel changes in response to sediment delivery, for a temperate upland gravel‐bed river. Results show that 16 months of measured in‐channel sedimentation in an upland gravel‐bed river cause about half of the increase in inundation extent that was simulated to arise from climate change. Consideration of the joint impacts of climate change and sedimentation emphasized the non‐linear nature of system response, and the possibly severe and synergistic effects that come from combined direct effects of climate change and sediment delivery. Such effects are likely to be exacerbated further as a result of the impacts of climate change upon coarse sediment delivery. In generic terms, these processes are commonly overlooked in flood risk mapping exercises and are likely to be important in any river system where there are high rates of sediment delivery and long‐term transfer of sediment to floodplain storage (i.e. alluviation involving active channel aggradation and migration). Similarly, attempts to reduce channel migration through river bank stabilization are likely to exacerbate this process as without bank erosion, channel capacity cannot be maintained. Finally, many flood risk mapping studies rely upon calibration based upon combining contemporary bed surveys with historical flood outlines, and this will lead to underestimation of the magnitude and frequency of floodplain inundation in an aggrading system for a flood of a given magnitude. Copyright © 2006 John Wiley & Sons, Ltd.     

Modelling artificial channel and land‐use changes and their impact on floods and sediment yield to the Ishikari basin

The potential for flooding and sediment transport is greatly affected by river channel form and changes in land use. Therefore the modelling of channel morphology prior to canalization and of land‐use change is important with respect to the prediction of floods and sediment yield and their consequences. A combination of land‐use transformation maps and soil properties shows certain decision rules for the conversion of forest into arable or vice versa. The model proposed, from this study, was used to simulate possible past and/or future channel and land‐use patterns. Subsequently, the outcome of this simulation was used to assess the risk of flooding, sediment transport and soil‐erosion under different conditions. In this study, channel morphology prior to canalization and land‐use change in the Ishikari basin, Hokkaido, Japan, were analysed by comparing three scenarios using a physical based channel and slope model. The results indicate that pre‐canalization channel morphology has a significant impact on flood peak, but no significant effect on sediment yield. In contrast, land‐use change has a significant effect on soil eroded from hillslopes, but no significant effect on flooding for Ishikari basin. This study also illustrates the challenges that a simple model, such as a physical based channel and slope model, can simulate large‐scale river basin processes using fewer hydrological data resources. Copyright © 2004 John Wiley & Sons, Ltd.     

Spatial patterns of sediment delivery to valley floors: sensitivity to sediment transport capacity and hillslope hydrology relations

There is considerable interest in large‐scale spatial patterns of sediment transport in catchments, and this topic is often approached using terrain‐based modelling. In such models topography influences the discharge of overland flow and its sediment transport capacity. The sediment transport capacity of overland flow is commonly expressed as a power function of slope and discharge (i.e. q_s=k₁q^βS^γ). The relationship between discharge and contributing area can also be expressed as a power function. Several reviews reveal a limited range of values for the two exponents β and γ. In this paper we examine the sensitivity of catchment‐scale patterns of sediment delivery to valley floors to a range of sediment transport capacity and hillslope hydrology parameterizations, using two catchments on the southern tablelands of New South Wales. The results indicate that, over the limited range of β and γ identified within the literature, sediment deliveries to valley floors across the two catchments are similar for all but one of five sediment transport capacity relationships. The patterns are dominated by the trend in slope through each catchment. The sensitivity to hillslope hydrology of predicted sediment delivery patterns is strong in the catchment with systematic variation in unit hillslope area, and weak in the catchment for which there are no systematic trends in unit hillslope area. We believe there is less experimental evidence to restrict choice of hillslope hydrology parameters than there is for sediment transport capacity. Copyright © 2001 John Wiley & Sons, Ltd.     

Runoff and sediment generation on bench‐terraced hillsides: measurements and up‐scaling of a field‐based model

Despite widespread bench‐terracing, stream sediment yields from agricultural hillsides in upland West Java remain high. We studied the causes of this lack of effect by combining measurements at different spatial scales using an erosion process model. Event runoff and sediment yield from two 4‐ha terraced hillside subcatchments were measured and field surveys of land use, bench‐terrace geometry and storage of sediment in the drainage network were conducted for two consecutive years. Runoff was 3·0–3·9% of rainfall and sediment yield was 11–30 t ha⁻¹ yr⁻¹ for different years, subcatchments and calculation techniques. Sediment storage changes in the subcatchment drainage network were less than 2 t ha⁻¹, whereas an additional 0·3–1·5 t ha⁻¹ was stored in the gully between the subcatchment flumes and the main stream. This suggests mean annual sediment delivery ratios of 86–125%, or 80–104% if this additional storage is included. The Terrace Erosion and Sediment Transport (TEST) model developed and validated for the studied environment was parameterized using erosion plot studies, land use surveys and digital terrain analysis to simulate runoff and sediment generation on the terraced hillsides. This resulted in over‐estimates of runoff and under‐estimates of runoff sediment concentration. Relatively poor model performance was attributed to sample bias in the six erosion plots used for model calibration and unaccounted covariance between important terrain attributes such as slope, infiltration capacity, soil conservation works and vegetation cover. Copyright © 2005 John Wiley & Sons, Ltd.     

Influence of anthropogenic land‐use change on hillslope erosion in the Waipaoa River Basin,New Zealand

European settlement of the Poverty Bay Region resulted in deforestation and conversion of > 90% of the landscape to pastureland. The resulting loss of vegetation triggered a rapid increase in hillslope erosion as widespread landslide complexes and gully systems developed on weak lithologic units in the Waipaoa Basin. To quantify the rate and volume of historic hillslope degradation, we used a 1956–2010 sequence of aerial photographs for a ~16 km² catchment to map temporal changes in the spatial extent of active landslides. Then we created a ‘turf index’ based on the extent and style of pastoral ground disruption, which correlates with downslope velocity. Based on the movement of trees and other features, we assigned average velocities to the turf classes as follows: (1) minimal disrupted ground: 0.6 m/yr, (2) a mix of disrupted ground and intact blocks: 3.4 m/yr, and (3) no intact blocks or vegetation: > 6 m/yr. We then calculated the average annual sediment flux using these turf‐derived velocities, the width of the landslide‐channel intersection, and an average toe depth of 4.4 ± 1.3 m (mean ± standard deviation [SD]) from 37 field measurements. The resulting catchment averaged erosion rates are (mean ± SD): 29.9 ± 12.9 mm/yr (1956), 28.8 ± 13.7 mm/yr (1969), 13.4 ± 4.9 mm/yr (1979), 17.0 ± 6.2 mm/yr (1988), and 9.9 ± 3.6 mm/yr (2010). Compared with long‐term (post‐18 ka) erosion rates (1.6 mm/yr) and the long‐term uplift rate (~1 mm/yr) for this site, the 50‐year anthropogenically‐driven rate is an order of magnitude larger (~20 mm/yr). Previously, we measured an increase in erosion over the past 3.4 kyr (2.2 mm/yr), and here, we demonstrate this increase could be primarily due to human land‐use change – showing that a century of rapid erosion superimposed on the background geologic rate can profoundly skew the interpretation of erosion rates. Copyright © 2016 John Wiley & Sons, Ltd.     

Entropy‐based equation to assess hillslope sediment production

A new equation to assess hillslope sediment production, based on physical and probabilistic approaches, is presented. The equation, which allows the computation of the delivery ratio for every event, considers the physical variables of travel distance, stream power, settling velocity and gross erosion. The probability density function that arises from the new formulation is solved using the principle of maximum entropy. Based on data from five watersheds in both tropical and temperate zones, the new delivery parameter K_v is calibrated and associated with vegetation cover and conservation practice. The proposed equation is rationally based in relation to parameter K_v. The entropy‐based equation was applied to assess sediment yield in two other experimental watersheds, showing good predictability for the set (mean absolute error of 20·8%). No systematic error was found in the analysed data. The entropy‐based equation showed good predictability for long‐term sets of data and for high‐erosivity events, but did not perform well for the low‐erosivity ones. Copyright © 2007 John Wiley & Sons, Ltd.     

Effects of subsurface soil characteristics on wetland–groundwater interaction in the coastal plain of the Chesapeake Bay watershed

Ecosystem services provided by depressional wetlands on the coastal plain of the Chesapeake Bay watershed (CBW) have been widely recognized and studied. However, wetland–groundwater interactions remain largely unknown in the CBW. The objective of this study was to examine the vertical interactions of depressional wetlands and groundwater with respect to different subsurface soil characteristics. This study examined two depressional wetlands with a low‐permeability and high‐permeability soil layer on the coastal plain of the CBW. The surface water level (SWL) and groundwater level (GWL) were monitored over 1 year from a well and piezometer at each site, respectively, and those data were used to examine the impacts of subsurface soil characteristics on wetland–groundwater interactions. A large difference between the SWL and GWL was observed at the wetland with a low‐permeability soil layer, although there was strong similarity between the SWL and GWL at the wetland with a high‐permeability soil layer. Our observations also identified a strong vertical hydraulic gradient between the SWL and GWL at the wetland with a high‐permeability soil layer relative to one with a low‐permeability soil layer. The hydroperiod (i.e., the total time of surface water inundation or saturation) of the wetland with a low‐permeability soil layer appeared to rely on groundwater less than the wetland with a high‐permeability soil layer. The findings showed that vertical wetland–groundwater interactions varied with subsurface soil characteristics on the coastal plain of the CBW. Therefore, subsurface soil characteristics should be carefully considered to anticipate the hydrologic behavior of wetlands in this region.     

Influence of storm-related sediment storage on the sediment delivery from tributary catchments in the upper Waipaoa River,New Zealand

Although much is known about overall sediment delivery ratios for catchments as components of sediment production and sediment yield, little is known about the component of temporary sediment storage. Sediment delivery ratios focused on the influence of storm-related sediment storage are measured at Matakonekone and Oil Springs tributaries of the Waipaoa River basin, east coast of New Zealand. The terrace deposits of both tributaries show abundant evidence of storm-related sedimentation, especially sediment delivered from Cyclone Bola, a 50 year return rainfall event which occurred in 1988. The sediment delivery ratio is calculated by dividing the volume of sediment transported from a tributary to the main stream by the volume of sediment generated at erosion sites in the tributary catchment. Because the sediment delivery volume is unknown, it can be calculated as the difference between sediment generation volume and sediment storage volume in the channel reach of the tributary. The volume of sediment generated from erosion sites in each tributary catchment was calculated from measurements made on aerial photographs dating from 1960 (1:44 000) and 1988 (1:27 000). The volume of sediment stored in the tributary can be calculated from measurements of cross-sections located along the tributary channel, which are accompanied by terrace deposits dated by counting annual growth rings of trees on terrace surfaces. Sediment delivery ratios are 0·93 for both Matakonekone catchment and Oil Springs catchment. Results indicate that Oil Springs catchment has contributed more than twice the volume of sediment to the Waipaoa River than the Matakonekone catchment (2·75 × 10⁶ m³ vs 1·22 × 10⁶ m³). Although large volumes of sediment are initially deposited during floods, subsequent smaller flows scour away much of these deposits. The sediment scouring rate from storage is 1·25 × 10⁴ m³ a⁻¹ for Matakonekone stream and 0·83 × 10⁴ m³ a⁻¹ for Oil Springs stream. Matakonekone and Oil Springs channels respond to extreme storms by instantaneously aggrading, then gradually excavating the temporarily stored sediment. Results from Matakonekone and Oil Springs streams suggest a mechanism by which event recurrence interval can strongly influence the magnitude of a geomorphic change. Matakonekone stream with its higher stream power is expected to excavate sediment deposits more rapidly and allow more rapid re-establishment of storage capacity. Copyright © 1999 John Wiley & Sons, Ltd.