An analysis of river bank slope and unsaturated flow effects on bank storage

Recognizing the underlying mechanisms of bank storage and return flow is important for understanding streamflow hydrographs. Analytical models have been widely used to estimate the impacts of bank storage, but are often based on assumptions of conditions that are rarely found in the field, such as vertical river banks and saturated flow. Numerical simulations of bank storage and return flow in river-aquifer cross sections with vertical and sloping banks were undertaken using a fully-coupled, surface-subsurface flow model. Sloping river banks were found to increase the bank infiltration rates by 98% and storage volume by 40% for a bank slope of 3.4° from horizontal, and for a slope of 8.5°, delay bank return flow by more than four times compared with vertical river banks and saturated flow. The results suggested that conventional analytical approximations cannot adequately be used to quantify bank storage when bank slope is less than 60° from horizontal. Additionally, in the unconfined aquifers modeled, the analytical solutions did not accurately model bank storage and return flow even in rivers with vertical banks due to a violation of the dupuit assumption. Bank storage and return flow were also modeled for more realistic cross sections and river hydrograph from the Fitzroy River, Western Australia, to indicate the importance of accurately modeling sloping river banks at a field scale. Following a single wet season flood event of 12 m, results showed that it may take over 3.5 years for 50% of the bank storage volume to return to the river.     

A review of sediment and nutrient concentration data from Australia for use in catchment water quality models

Land use (and land management) change is seen as the primary factor responsible for changes in sediment and nutrient delivery to water bodies. Understanding how sediment and nutrient (or constituent) concentrations vary with land use is critical to understanding the current and future impact of land use change on aquatic ecosystems. Access to appropriate land-use based water quality data is also important for calculating reliable load estimates using water quality models. This study collated published and unpublished runoff, constituent concentration and load data for Australian catchments. Water quality data for total suspended sediments (TSS), total nitrogen (TN) and total phosphorus (TP) were collated from runoff events with a focus on catchment areas that have a single or majority of the contributing area under one land use. Where possible, information on the dissolved forms of nutrients were also collated. For each data point, information was included on the site location, land use type and condition, contributing catchment area, runoff, laboratory analyses, the number of samples collected over the hydrograph and the mean constituent concentration calculation method. A total of ～750 entries were recorded from 514 different geographical sites covering 13 different land uses. We found that the nutrient concentrations collected using "grab" sampling (without a well defined hydrograph) were lower than for sites with gauged auto-samplers although this data set was small and no statistical analysis could be undertaken. There was no statistically significant difference (p<0.05) between data collected at plot and catchment scales for the same land use. This is most likely due to differences in land condition over-shadowing the effects of spatial scale. There was, however, a significant difference in the concentration value for constituent samples collected from sites where >90% of the catchment was represented by a single land use, compared to sites with <90% of the upstream area represented by a single land use. This highlights the need for more single land use water quality data, preferably over a range of spatial scales. Overall, the land uses with the highest median TSS concentrations were mining (～50,000mg/l), horticulture (～3000mg/l), dryland cropping (～2000mg/l), cotton (～600mg/l) and grazing on native pastures (～300mg/l). The highest median TN concentrations are from horticulture (～32,000μg/l), cotton (～6500μg/l), bananas (～2700μg/l), grazing on modified pastures (～2200μg/l) and sugar (～1700μg/l). For TP it is forestry (～5800μg/l), horticulture (～1500μg/l), bananas (～1400μg/l), dryland cropping (～900mg/l) and grazing on modified pastures (～400μg/l). For the dissolved nutrient fractions, the sugarcane land use had the highest concentrations of dissolved inorganic nitrogen (DIN), dissolved organic nitrogen (DON) and dissolved organic phosphorus (DOP). Urban land use had the highest concentrations of dissolved inorganic phosphorus (DIP). This study provides modellers and catchment managers with an increased understanding of the processes involved in estimating constituent concentrations, the data available for use in modelling projects, and the conditions under which they should be applied. Areas requiring more data are also discussed.     

Alignment time-scale of the microquasar GRO J1655−40

The microquasar GRO J1655−40 has a black hole with spin angular momentum apparently misaligned to the orbital plane of its companion star. We analytically model the system with a steady-state disc warped by Lense–Thirring precession and find the time-scale for the alignment of the black hole with the binary orbit. We make detailed stellar evolution models so as to estimate the accretion rate and the lifetime of the system in this state. The secondary can be evolving at the end of the main sequence or across the Hertzsprung gap. The mass-transfer rate is typically 50 times higher in the latter case but we find that, in both the cases, the lifetime of the mass-transfer state is at most a few times the alignment time-scale. The fact that the black hole has not yet aligned with the orbital plane is therefore consistent with either model. We conclude that the system may or may not have been counter aligned after its supernova kick but that it is most likely to be close to alignment rather than counter alignment now.     

In situ characterization of grain‐scale fluvial morphology using Terrestrial Laser Scanning

The grain‐scale morphology of fluvial sediments is an important control on the character and dynamics of river systems; however current understanding of its role is limited by the difficulties of robustly quantifying field surface morphology. Terrestrial Laser Scanning (TLS) offers a new methodology for the rapid acquisition of high‐resolution and high‐precision surface elevation data from in situ sediments. To date, most environmental and fluvial applications of TLS have focused on large‐scale systems, capturing macroscale morphologies. Application of this new technology at scales necessary to characterize the complexity of grain‐scale fluvial sediments therefore requires a robust assessment of the quality and sources of errors in close‐range TLS data. This paper describes both laboratory and field experiments designed to evaluate close‐range TLS for sedimentological applications and to develop protocols for data acquisition. In the former, controlled experiments comprising high‐resolution scans of white, grey and black planes and a sphere were used to quantify the magnitude and source of three‐dimensional (3D) point errors resulting from a combination of surface geometry, reflectivity effects and inherent instrument precision. Subsequently, a methodology for the collection and processing of grain‐scale TLS data is described through an application to a coarse grained gravel system, the River Feshie (D₅₀ 32 to 63 mm). This stepwise strategy incorporates averaging repeat scans and filtering scan artefact and non‐surface points using local 3D search algorithms. The sensitivity of the results to the filter parameter values are assessed by careful internal validation of Digital Terrain Models (DTMs) created from the resulting point cloud data. The transferability of this methodology is assessed through application to a second river, Bury Green Brook, dominated by finer gravel (D₅₀ 18 to 33 mm). The factor limiting the resolution of DTMs created from this second dataset was found to be the relative sizes of the laser footprint and smallest grains. Copyright © 2009 John Wiley & Sons, Ltd.     

Changes in sediment and organic carbon accumulation in a highly-disturbed ecosystem: The Sacramento-San Joaquin River Delta (California, USA)

We used the Sacramento-San Joaquin River Delta CA (Delta, hereafter) as a model system for understanding how human activities influence the delivery of sediment and total organic carbon (TOC) over the past 50–60 years. Sediment cores were collected from sites within the Delta representing the Sacramento River (SAC), the San Joaquin River (SJR), and Franks Tract (FT), a flooded agricultural tract. A variety of anthropogenic tracers including ¹³⁷Cs, total DDE (∑DDE) and brominated diphenyl ether (BDE) congeners were used to quantify sediment accumulation rates. This information was combined with total organic carbon (TOC) profiles to quantify rates of TOC accumulation. Across the three sites, sediment and TOC accumulation rates were four to eight-fold higher prior to 1972. Changes in sediment and TOC accumulation were coincident with completion of several large reservoirs and increased agriculture and urbanization in the Delta watershed. Radiocarbon content of TOC indicated that much of the carbon delivered to the Delta is “pre-aged” reflecting processing in the Delta watershed or during transport to the sites rather than an input of predominantly contemporary carbon (e.g., 900–1400 years BP in surface sediments and 2200 yrs BP and 3610 yrs BP at the base of the SJR and FT cores, respectively). Together, these data suggest that human activities have altered the amount and age of TOC accumulating in the Delta since the 1940s.     

Using MODFLOW 2000 to model ET and recharge for shallow ground water problems

In environments with shallow ground water elevation, small changes in the water table can cause significant variations in recharge and evapotranspiration fluxes. Particularly, where ground water is close to the soil surface, both recharge and evapotranspiration are regulated by a thin unsaturated zone and, for accuracy, must be represented using nonconstant and often nonlinear relationships. The most commonly used ground water flow model today, MODFLOW, was originally designed with a modular structure with independent packages representing recharge and evaporation processes. Systems with shallow ground water, however, may be better represented using either a recharge function that varies with ground water depth or a continuous recharge and evapotranspiration function that is dependent on depth to water table. In situations where the boundaries between recharging and nonrecharging cells change with time, such as near a seepage zone, a continuous ground water flux relationship allows recharge rates to change with depth rather than having to calculate them at each stress period. This research article describes the modification of the MODFLOW 2000 recharge and segmented evapotranspiration packages into a continuous recharge-discharge function that allows ground water flux to be represented as a continuous process, dependent on head. The modifications were then used to model long-term recharge and evapotranspiration processes on a saline, semiarid floodplain in order to understand spatial patterns of salinization, and an overview of this process is given.     

Formation and erosion of sediment cover in an experimental bedrock‐alluvial channel

Sediment grains in a bedrock‐alluvial river will be deposited within or adjacent to a sediment patch, or as isolated grains on the bedrock surface. Previous analysis of grain geometry has demonstrated that these arrangements produce significant differences in grain entrainment shear stress. However, this analysis neglected potential interactions between the sediment patches, local hydraulics and grain entrainment. We present a series of flume experiments that measure the influence of sediment patches on grain entrainment. The flume had a planar bed with roughness that was much smaller than the diameters of the mobile grains. In each experiment sediment was added either as individual grains or as a single sediment pulse. Flow was then increased until the sediment was entrained. Analysis of the experiments demonstrates that: (1) for individual grains, coarse grains are entrained at a higher discharge than fine grains; (2) once sediment patches are present, the different in entrainment discharge between coarse and fine grains is greatly reduced; (3) the sheltering effect of patches also increases the entrainment discharge of isolated grains; (4) entire sediment patches break‐up and are eroded quickly, rather than through progressive grain‐by‐grain erosion; (5) as discharge increases there is some tendency for patches to become more elongate and flow‐aligned, and more randomly distributed across the bed. One implication of this research is that the critical shear stress in bedrock‐alluvial channels will be a function of the extent of the sediment cover. Another is that the influence of sediment patches equalizes critical shear stresses between different grain sizes and grain locations, meaning that these factors may not need to be accounted for. Further research is needed to quantify interactions between sediment patches, grain entrainment and local hydraulics on rougher bedrock surfaces, and under different types of sediment supply. Copyright © 2016 John Wiley & Sons, Ltd.     

Evapotranspiration and canopy characteristics of two lodgepole pine stands following mountain pine beetle attack

Over the past decade, British Columbia (BC), has experienced the largest mountain pine beetle (MPB) outbreak on record. This study used the eddy‐covariance (EC) technique to examine the impact of the MPB attack on evapotranspiration (E) and associated canopy characteristics of two lodgepole pine stands with secondary structure (trees, saplings and seedlings surviving the attack) located in central BC. MPB‐06, an 85‐year‐old almost pure stand of pine trees, was first attacked in 2006, and by 2010, ~80% of the trees had been killed. MPB‐03, a 110‐year‐old stand with an overstory consisting of over 90% pine and a developed sub‐canopy, was first attacked in 2003 and by 2007 had > 95% pine canopy mortality. EC measurements began in August 2006 at MPB‐06 and in March 2007 at MPB‐03, and continued for four years. Annual total E ranged from 226 mm to 237 mm at MPB‐06, and from 280 to 297 mm at MPB‐03, showing relatively little year‐to‐year change at both sites over the four years. Increased E from the accelerated growth of the surviving vegetation (secondary structure, shrubs and herbs) compensated for reduction in E due to the death of the overstory. Monthly average daytime canopy conductance, the Priestley–Taylor (α), and the canopy–atmosphere decoupling coefficient (Ω) steadily increased during the growing season reaching approximate maximum values of 5 mm s^?1, 0.75 and 0.12, respectively. Potential evapotranspiration was approximated using a vapour pressure deficit‐dependent α obtained at high soil water content. Calculated water deficits indicated some water‐supply limitation to the surviving trees and understory at both sites. Rates of root zone drainage during the growing season were low relative to precipitation. Copyright © 2013 John Wiley & Sons, Ltd.     

Carbon losses in soils previously exposed to elevated atmospheric CO2 in a chaparral ecosystem: Potential implications for a sustained biospheric C sink

Between 1996 and 2001 an experimental set up in a chaparral community near San Diego, CA, examined various plant and ecosystem responses to CO₂ concentrations ranging from 250 to 750 μl l^− 1. These experiments indicated a significant increase in soil C sequestration as CO₂ rose above the ambient levels. In 2003, two years after the cessation of the CO₂ treatments, we returned to this site to examine soil C dynamics with a particular emphasis on stability of specific pools of C. We found that in as little as two years, C content in the surface soils (0–15 cm) of previously CO₂ enriched plots had dropped to levels below those of the ambient and pretreatment soils. In contrast, C retained in response to CO₂ enrichment was more durable in the deeper soil layers (> 25 cm deep) where both organic and inorganic C were on average 26% and 55% greater, respectively, than C content of ambient plots. Using stable isotope tracers, we found that treatment C represented 25% of total soil C and contributed to 55% of soil CO₂ efflux, suggesting that most of treatment C is readily accessible to decomposers. We also found that, C present before CO₂ fumigation was decomposed at a faster rate in the plots that were exposed to elevated CO₂ than in those exposed to ambient CO₂ levels. To our knowledge, this is the first report that allows for a detail accounting of soil C after ceasing CO₂ treatments. Our study provides a unique insight to how stable the accrued soil C is as CO₂ increases in the atmosphere.