An analysis of soil crust strength in relation to potential abrasion by saltating particles

Wind erosion depends on the ease with which particles can be detached from the soil surface, but suitable tests to characterize this property are not available. Two possible methods to determine surface soil strength in the field were therefore compared on a range of artificially ‘crusted’ surfaces. These were made by spraying or tension wetting aggregates (10–2, 2–0.5 and <0.5 mm) from a structurally unstable sandy loam, followed by drying. Each test involved measuring the force exerted on a probe driven at a steady rate into the surface, using either a flat-tipped 0.6 mm diameter penetrometer or a flat-ended cylindrical punch with inner and outer diameters of 5 and 6 mm, respectively. Both probes showed that crusts could be produced reproducibly. Depending on the probe and aggregate size, penetration mainly occurred either as a result of aggregates being deflected out of the pathway of the probe or by genuine rupture of aggregates or of the crusted surface. The penetrometer, because it was comparable to the size of sand grains, gave results that can be used to characterize surface erodibility to saltating particles. The punch gave results that would be unsuitable for this purpose, as would other strength tests that are on too large a scale. Penetrometer results were analyzed to calculate the energy required for penetration. It was thus possible to demonstrate that only the spray-wetted fine aggregates had a surface that could undergo large-scale rupture by saltating sand grains. For all other surfaces, saltating particles would be unable to supply sufficient energy to rupture aggregates or the crusted surface. Erosion could only occur by a slower process of abrasion in which smaller particles or aggregates are chipped away from the surface. However, it is shown that saltating particles could rupture the interaggregate bonding in the 2–0.5 mm aggregate surfaces, thus permitting creep. An alternative and potentially simpler way of characterizing surface erodibility by using a surface modulus of elasticity is also discussed. Our results demonstrate that the small diameter penetrometer is a promising technique for characterizing erodibility of aggregated and crusted surfaces. © 1997 by John Wiley & Sons, Ltd.     

Signal crayfish burrowing,bank retreat and sediment supply to rivers: A biophysical sediment budget

Burrowing into riverbanks by animals transfers sediment directly into river channels and has been hypothesised to accelerate bank erosion and promote mass failure. A field monitoring study on two UK rivers invaded by signal crayfish (Pacifastacus leniusculus) assessed the impact of burrowing on bank erosion processes. Erosion pins were installed in 17 riverbanks across a gradient of crayfish burrow densities and monitored for 22 months. Bank retreat increased significantly with crayfish burrow density. At the bank scale (<6 m river length), high crayfish burrow densities were associated with accelerated bank retreat of up to 253% and more than a doubling of the area of bank collapse compared with banks without burrows. Direct sediment supply by burrowing activity contributed 0.2% and 0.6% of total sediment at the reach (1.1 km) and local bank (<6 m) scales. However, accelerated bank retreat caused by burrows contributed 12.2% and 29.8% of the total sediment supply at the reach and bank scales. Together, burrowing and the associated acceleration of retreat and collapse supplied an additional 25.4 t km⁻¹ a⁻¹ of floodplain sediments at one site, demonstrating the substantial impact that signal crayfish can have on fine sediment supply. For the first time, an empirical relation linking animal burrow characteristics to riverbank retreat is presented. The study adds to a small number of sediment budget studies that compare sediment fluxes driven by biotic and abiotic energy but is unique in isolating and measuring the substantial interactive effect of the acceleration of abiotic bank erosion facilitated by biotic activity. Biotic energy expended through burrowing represents an energy surcharge to the river system that can augment sediment erosion by geophysical mechanisms.     

Biodegradation and environmental behavior of biodiesel mixtures in the sea: An initial study

Biodiesel, a mixture of fatty acid methyl esters (FAMEs) derived from animal fats or vegetable oils, is rapidly moving towards the mainstream as an alternative source of energy. However, the behavior of biodiesel, or blends of biodiesel with fossil diesel, in the marine environment have yet to be fully understood. Hence, we performed a series of initial laboratory experiments and simple calculations to evaluate the microbial and environmental fate of FAMEs. Aerobic seawater microcosms spiked with biodiesel or mixtures of biodiesel and fossil diesel revealed that the FAMEs were degraded at roughly the same rate as n-alkanes, and more rapidly than other hydrocarbon components. The residues extracted from these different microcosms became indistinguishable within weeks. Preliminary results from physical-chemical calculations suggest that FAMEs in biodiesel mixtures will not affect the evaporation rates of spilled petroleum hydrocarbons but may stabilize oil droplets in the water column and thereby facilitate transport.     

Grain size along two gravel-bed rivers: statistical variation,spatial pattern and sedimentary links

A new set of field data facilitates a detailed analysis of variations in bed material grain size within two confluent gravel-bed rivers in northeastern British Columbia, Canada. A preliminary assessment of grain-size variability establishes a basis for examination of the spatial pattern of grain-size change. Standard ANOVA techniques are inappropriate because individual samples have unequal variances and are not normally distributed. Alternative tests for homoscedasticity and comparison of means are therefore utilized. Within-site, between-sample variability is not significant. The grain-size distributions that were obtained at individual sites are therefore representative of the depositional environments that were sampled. In both rivers mean grain size does vary significantly between sites and there is therefore a basis for examining the data for spatial patterns such as downstream fining. Textural variations along the two rivers studied here are complex and show negligible overall fining (in over 100 km). This is the consequence of a large number of tributary inputs and non-alluvial sediment sources which are the legacy of Late Pleistocene glaciation. The identification of lateral sources like these is fundamental for understanding textural changes within rivers. The sedimentary link (a channel reach between significant lateral sediment inputs) provides a means of isolating fluvial maturation processes (abrasion and sorting) from contingent lateral inputs. Strong fining trends are apparent in most links and classification of grain-size measurements according to their location within particular links greatly improves the statistical explanation of textural variation. Identification of sedimentary links provides a means of applying models of fluvial fining processes, so isolation of link networks will aid the development of basin-scale models of textural variation. © 1998 John Wiley & Sons, Ltd.     

An Analysis Platform for Multiscale Hydrogeologic Modeling with Emphasis on Hybrid Multiscale Methods

One of the most significant challenges faced by hydrogeologic modelers is the disparity between the spatial and temporal scales at which fundamental flow, transport, and reaction processes can best be understood and quantified (e.g., microscopic to pore scales and seconds to days) and at which practical model predictions are needed (e.g., plume to aquifer scales and years to centuries). While the multiscale nature of hydrogeologic problems is widely recognized, technological limitations in computation and characterization restrict most practical modeling efforts to fairly coarse representations of heterogeneous properties and processes. For some modern problems, the necessary level of simplification is such that model parameters may lose physical meaning and model predictive ability is questionable for any conditions other than those to which the model was calibrated. Recently, there has been broad interest across a wide range of scientific and engineering disciplines in simulation approaches that more rigorously account for the multiscale nature of systems of interest. In this article, we review a number of such approaches and propose a classification scheme for defining different types of multiscale simulation methods and those classes of problems to which they are most applicable. Our classification scheme is presented in terms of a flowchart (Multiscale Analysis Platform), and defines several different motifs of multiscale simulation. Within each motif, the member methods are reviewed and example applications are discussed. We focus attention on hybrid multiscale methods, in which two or more models with different physics described at fundamentally different scales are directly coupled within a single simulation. Very recently these methods have begun to be applied to groundwater flow and transport simulations, and we discuss these applications in the context of our classification scheme. As computational and characterization capabilities continue to improve, we envision that hybrid multiscale modeling will become more common and also a viable alternative to conventional single‐scale models in the near future.     

Constitutive relations for fault slip and earthquake instabilities

Constitutive relations for fault slip are described and adopted as a basis for analyzing slip motion and its instability in the form of earthquakes on crustal faults. The constitutive relations discussed include simple rate-independent slip-weakening models, in which shear strength degrades with ongoing slip to a residual frictional strength, and also more realistic but as yet less extensively applied slip-rate and surface-state-dependent relations. For the latter the state of the surface is characterized by one or more variables that evolve with ongoing slip, seeking values consistent with the current slip rate. Models of crustal faults range from simple, single-degree-of-freedom spring-slider systems to more complex continuous systems that incorporate nonuniform slip and locked patches on faults of depth-dependent constitutive properties within elastic lithospheric plates that may be coupled to a viscoelastic asthenosphere. Most progress for the rate and state-dependent constitutive relations is at present limited to single-degree-of-freedom systems. Results for stable and unstable slip with the various constitutive models are summarized. Instability conditions are compared for spatially uniform versus nonuniform slip, including the elastic — brittle crack limit of the nonuniform mode. Inferences of constitutive and fracture parameters are discussed, based on earthquake data for large ruptures that begin with slip at depth, concentrating stress on locked regions within a brittle upper crust. Results of nonlinear stability theory, including regimes of complex sustained stress and slip rate oscillations, are outlined for rate and state-dependent constitutive relations, and the manner in which these allow phenomena like time-dependent failure, restrengthening in nearly stationary contact, and weakening in rapidly accelerated slip, is discussed.     

Bed material texture in low order streams on the Queen Charlotte Islands,British Columbia

Low order channels comprise a large proportion of the links of every drainage basin, and are often at the centre of land management concerns. They exhibit hydrological and geomorphological characteristics atypical of higher order links. This paper examines the nature and causes of variations in the bed material texture of two streams on the Queen Charlotte Islands, British Columbia. The extant, functional exponential model is found to be inadequate for explaining observed changes in grain size parameters with distance downstream. Recurrent disruption of sediment transport by large organic debris jams, and the sporadic contamination of the fluvial sediment population by colluvial inputs, preclude the development of longitudinal structure. Rather, grain size varies erratically over short distances. A stochastic model best describes the observed variations, and should be adopted as an alternative to the exponential model in low order links. Characteristic variances are controlled by the degree of hillslope-channel coupling, and the extent and characteristics of non-alluvial storage mechanisms.     

The use of reference models from a priori data to guide 2D inversion of electrical resistivity tomography data

Clay-rich till plains cover much of the UK. Such sites are attractive locations for landfills, since clay aquitards lower the risk of landfill leachate entering groundwater. However, such tills often contain sand and gravel bodies that can act as leachate flow routes. Such bodies may not be detected by conventional site investigation techniques such as drilling boreholes and trial pitting. A method of guided inversion, where a priori data are used to construct structural reference models for use in inverting electrical resistivity tomography data, was proposed as a tool to improve the detection of sand and gravel bodies within clay-rich till sequences.
Following a successful 2D guided inversion synthetic modelling study, a field study was undertaken. Wenner 2D electrical resistivity tomography lines, resistivity cone penetrometry bores and electromagnetic induction ground resistivity data were collected over a site on the East Yorkshire coast, England, where sand and gravel lenses were known to exist from cliff exposures. A number of equally valid geoelectrical models were constructed using the electromagnetic and resistivity cone data. These were used as structural reference models in the inversion of the resistivity tomography data. Blind inversion using an homogenous reference model was also carried out for comparison.
It was shown for the first time that the best solution model produced by 2D inversion of one data set with a range of structural reference models could be determined by using the l ₂ model misfit between the solution models and associated reference models (reference misfit) as a proxy for the l ₂ misfit between the solution models and the synthetic model or 'best-guess' geoelectrical model (true misfit). The 2D methodology developed here is applicable in clay-rich till plains containing sand and gravel bodies throughout the UK.     

Form and growth of bars in a wandering gravel‐bed river

The changing form of developing alluvial river bars has rarely been studied in the field, especially in the context of the fixed, compound, mainly alternate gravel bars that are the major morphological feature of the wandering style. Century scale patterns of three‐dimensional growth and development, and the consequent scaling relations of such bars, are examined along the gravel‐bed reach of lower Fraser River, British Columbia, Canada. A retrospective view based on maps and aerial photographs obtained through the twentieth century shows that individual bars have a life history of about 100 years, except in certain, protected positions. A newly formed gravel bar quickly assumes its ultimate thickness and relatively quickly approaches its equilibrium length. Growth continues mainly by lateral accretion of unit bars, consistent with the lateral style of instability of the river. Bar growth is therefore allometric. Mature bars approach equilibrium dimensions and volume that scale with the overall size of the channel. Accordingly, the bars conform with several published criteria for the ultimate dimensions of alternate barforms. Sand bars, observed farther downstream, have notably different morphology. Fraser River presents a typical wandering channel planform, exhibiting elements of both meandered and low‐order braided channels. Hydraulic criteria to which the Fraser bars conform illustrate why this planform develops and persists. The modest rate of bed material transfer along the channel – typical of the wandering type – determines a century‐length time scale for bar development. This time scale is consistent with estimates that have been made for change of the macroform elements that determine the overall geometry of alluvial channels. Copyright © 2009 John Wiley & Sons, Ltd.     

A GIS-based multi-proxy analysis of the evolution of subglacial dynamics of the Quebec–Labrador ice dome,northeastern Quebec,Canada

The glacial records of the inner-core regions of the Laurentide Ice Sheet (LIS) document complex yet coherent patterns reflecting ice-sheet change (e.g. ice-divide migration), providing unique insights into past glacial conditions. This study develops a conceptual model of subglacial dynamics evolution within a major ice-dispersal centre of the LIS in northeastern Quebec, Canada using a GIS-based analysis of the surficial geologic record. Multiple proxies of subglacial conditions (subglacial streamlined landforms, lake density and lake area over thin drift/bedrock) were analysed through grid-overlay techniques and then classified based on different proxy variables ranging from highly mobile warm-based to immobile cold-based conditions. An additional proxy (till blanket) was used to identify areas of thick till deposition, but with few proxies (few lake or landform metrics). Based on local ice-flow reconstructions, the most ‘relict’ glacial terrain zone (GTZ1) has warm-based conditions over 66% of its area and is remarkably well preserved, suggesting laterally extensive warm-based conditions during the oldest identified ice-flow phase. This relict glacial terrain is partially overprinted by two subsequent ice-flow phases in spatially restricted zones in the northeast (73% warm-based), east-central (41% warm-based), and northwest (33% warm-based) of the study area. A zone of more sluggish conditions (only 3% warm-based) was identified in the highlands at the centre of the study area, characterized by thin till cover, few landforms, yet with large patches of relatively abundant small lakes, indicative of areal scouring. No clear evidence of sustained cold-based conditions (i.e. high chemical index of alteration values or high ¹⁰Be abundances) was found in the study area. These results suggest that warm-based conditions (active erosion and/or deposition) were uniformly widespread during the earliest ice-flow phase, later becoming more spatially restricted with broader sluggish ice conditions. These spatially restricted regions of warm-based subglacial regimes were likely controlled by surrounding and down-flow ice streaming. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd