Quantifying the dynamics of microtopography during a snowmelt event

Knowledge of soil microtopography and its changes in space and over time is important to the understanding of how tillage influences infiltration, runoff generation and erosion. In this study, the use of a terrestrial laser scanner (TLS) is assessed for its ability to quantify small changes in the soil surface at high spatial resolutions for a relatively large surface area (100 m²). Changes in soil surface morphology during snow cover and melt are driven by frost heave, slaking, pressure exertion by the snowpack and overland flow (erosion and deposition). An attempt is undertaken to link these processes to observed changes at the soil surface. A new algorithm for soil surface roughness is introduced to make optimal use of the raw point cloud. This algorithm is less scale dependent than several commonly used roughness calculations. The results of this study show that TLSs can be used for multitemporal scanning of large surfaces and that small changes in surface elevation and roughness can be detected. Statistical analysis of the observed changes against terrain indices did not yield significant evidence for process differentiation. © 2019 The Authors. Earth Surface Processes and Landforms Published by John Wiley & Sons Ltd. © 2019 The Authors. Earth Surface Processes and Landforms Published by John Wiley & Sons Ltd.     

A physical model of the high-frequency seismic signal generated by debris flows

We propose a physical model for the high-frequency (>1 Hz) spectral distribution of seismic power generated by debris flows. The modeled debris flow is assumed to have four regions where the impact rate and impulses are controlled by different mechanisms: the flow body, a coarser-grained snout, a snout lip where particles fall from the snout on the bed, and a dilute front composed of saltating particles. We calculate the seismic power produced by this impact model in two end-member scenarios, a thin-flow and thick-flow limit, which assume that the ratio of grain sizes to flow thicknesses are either near unity or much less than unity. The thin-flow limit is more appropriate for boulder-rich flows that are most likely to generate large seismic signals. As a flow passes a seismic station, the rise phase of the seismic amplitude is generated primarily by the snout while the decay phase is generated first by the snout and then the main flow body. The lip and saltating front generate a negligible seismic signal. When ground properties are known, seismic power depends most strongly on both particle diameter and average flow speed cubed, and also depends on length and width of the flow. The effective particle diameter for producing seismic power is substantially higher than the median grain size and close to the 73rd percentile for a realistic grain size distribution. We discuss how the model can be used to estimate effective particle diameter and average flow speed from an integrated measure of seismic power. © 2019 The Authors. Earth Surface Processes and Landforms Published by John Wiley & Sons Ltd. © 2019 The Authors. Earth Surface Processes and Landforms Published by John Wiley & Sons Ltd.     

Fate of pioneering vegetation patches in a dynamic meandering river

Fluvial biomorphodynamics in actively meandering rivers entail interactions between hydromorphodynamics and pioneering tree species that have eco-engineering effects. Here we study spatiotemporal patterns of vegetation patches smaller than 150 m² in a 4 km reach of the river Allier in France in order to unravel causes for tree persistence and mortality and identify spatial trends across the river valley. To this end we analysed aerial photographs by object-based image analysis over a period of 56 years and tracked individual patches through time. Furthermore the cover and surface age of the study reach were classified. The large-scale shifts of channels, bars and vegetation are consistent with the meandering process and chute cutoffs. However, the spatiotemporal patterns of the vegetation patches are surprising in that they are ubiquitous and have ages up to decades on the highly dynamic meander belt, but hardly expand into larger vegetation patches. Patches disappear exponentially as a function of their age, and faster so in the last decades. Causes are amalgamation into the riparian forest flanking the meander belt and mortality likely due to desiccation or erosion. Patches have a higher probability of survival when further away from the active channel and closer to high vegetation patches and valley boundary. The window of opportunity of vegetation settlement widens towards the valley boundaries and in floodplain lows of former channels and chutes. These results imply a gradual cross-valley gradient of riparian vegetation settling, survival and succession. © 2019 The Authors. Earth Surface Processes and Landforms Published by John Wiley & Sons Ltd.     

Complex disturbance-driven reactivation of near-surface sediments in the largest dunefield in North America during the last 200 years

Understanding the historical activity of desert dune systems is important for identifying both the palaeoenvironmental drivers of change and the likelihood of future reactivation. Dating dune sediments in the Nebraska Sandhills has identified regional-scale dune activity over centennial and millennial timescales during the Holocene, occurring at 9.6–6.5, 3.8 and 2.5 ka, and most recently spanning the Medieval Climatic Anomaly 1050–650 years BP. These periods have been interpreted as palaeoclimatic evidence of intense aridity lasting decadal and centennial timescales. A detailed record of dune activity in the historical period, since EuroAmerican arrival, is lacking however, yet important for interpreting the role of human agency amongst the factors influencing disturbance. Without a high-resolution record of short-term, historical, local sediment mobilization, it is not possible to distinguish the environmental factor(s) responsible for local reactivation. In this paper, the individual drivers of vegetation disturbance are reviewed and presented alongside a luminescence-dated reconstruction of dune sediment deposition ages. This allows an integrated assessment of the relationship between drivers and environmental response over a recorded period. We focused our investigation on the aeolian reactivations of surface dune sediments and blowout features around the Niobrara Valley Preserve in the northern limits of the Nebraska Sandhills. Results show a near-continuous (within uncertainties) timeline of local reactivation across the sites studied, with variation between the individual features indicating that both regional (i.e. climatic) and local (i.e. land use) forcings contribute to surface disturbance. © 2019 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd. © 2019 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Plant root and rhizome strength: Are there differences between and within species and rivers?

The strength and architecture of roots and other below-ground organs of riparian and aquatic plants affect plant resistance to uprooting and contribute to reinforcing river bank, bar and bed materials. Therefore, root properties are an important element in models for estimating river bank stability and such models may focus on the role of plants by using root strength–diameter relationships for the particular plant species that are present. Here we explore the degree to which there appear to be significant differences in strength–diameter relationships between and within species-specific data sets obtained for two riparian tree/shrub (Populus nigra, Salix alba) and two emergent aquatic macrophyte (Sparganium erectum, Phalaris arundinacea) species in different European river environments. While the analysed data sets were not specifically collected to answer these research questions, the results are sufficiently compelling to make the case for the collection of a more comprehensive data set and its rigorous analysis. This would allow recommendations to be made on the degree to which (i) species-specific or more general relationships between root/rhizome strength and diameter are appropriate, (ii) such relationships are applicable within and between rivers in different geographical regions and subject to different local environmental conditions, and (iii) further (minimalist) field observations are needed to calibrate such relationships for investigations of new locales or species. © 2018 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Distinct patterns of bank erosion in a navigable regulated river

Distinct bankline patterns appeared after the removal of protection works along a navigable reach of the Meuse River. A series of oblique embayments now dominate the riverine landscape after ten years of bank erosion, but their location and asymmetry cannot be explained yet. This work analyses and integrates field measurements of flow, ship waves, bank composition, bed topography and historical maps to explain the observed patterns along two reaches of the river. An extraordinary low-water-level event generated by a ship accident provided the unique opportunity to also analyse the subaqueous bank topography. The results indicate that the formation of oblique embayments arises from the combination of floodplain heterogeneity, structured by scroll-bar deposits, and the regulation of water levels, resulting in ship-wave attack at a narrow range of bank elevation for 70% of the time. Substrate erodibility acts on the effectiveness of trees to slow down local bank erosion rates, which is possibly enhanced by a positive feedback between woody roots and cohesive soil. The strong regulation of water levels and the waves generated by the intense ship traffic produce an increasingly long mildly-sloping terrace at the bank toe and progressively dominate the bank erosion process. This study demonstrates the important role of floodplain and scroll bar formation in shaping later bank erosion, which has implications for predictive numerical models, restoration strategies, and understanding the role of vegetation in bank erosion processes. © 2019 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd     

Advances in computational morphodynamics using the International River Interface Cooperative (iRIC) software

Results from computational morphodynamics modeling of coupled flow–bed–sediment systems are described for 10 applications as a review of recent advances in the field. Each of these applications is drawn from solvers included in the public-domain International River Interface Cooperative (iRIC) software package. For mesoscale river features such as bars, predictions of alternate and higher mode river bars are shown for flows with equilibrium sediment supply and for a single case of oversupplied sediment. For microscale bed features such as bedforms, computational results are shown for the development and evolution of two-dimensional bedforms using a simple closure-based two-dimensional model, for two- and three-dimensional ripples and dunes using a three-dimensional large-eddy simulation flow model coupled to a physics-based particle transport model, and for the development of bed streaks using a three-dimensional unsteady Reynolds-averaged Navier–Stokes solver with a simple sediment-transport treatment. Finally, macroscale or channel evolution treatments are used to examine the temporal development of meandering channels, a failure model for cantilevered banks, the effect of bank vegetation on channel width, the development of channel networks in tidal systems, and the evolution of bedrock channels. In all examples, computational morphodynamics results from iRIC solvers compare well to observations of natural bed morphology. For each of the three scales investigated here, brief suggestions for future work and potential research directions are offered. © 2019 The Authors Earth Surface Processes and Landforms Published by John Wiley & Sons Ltd     

Using lake sediment archives to improve understanding of flood magnitude and frequency: Recent extreme flooding in northwest UK

We present the first quantitative reconstruction of palaeofloods using lake sediments for the UK and show that for a large catchment in NW England the cluster of devastating floods from 1990 to present is without precedent in this 558-year palaeo-record. Our approach augments conventional flood magnitude and frequency (FMF) analyses with continuous lake sedimentary data to provide a longer-term perspective on flood magnitude recurrence probabilities. The 2009 flood, the largest in >558 years, had a recurrence interval larger (1:2,200 year) than revealed by conventional flood estimation using shorter duration gauged single station records (1:1,700 year). Flood-rich periods are non-stationary in their correlation with climate indices, but the 1990-2018 cluster is associated with warmer Northern Hemisphere Temperatures and positive Atlantic Multidecadal Oscillation. Monitored records rarely capture the largest floods and our palaeoflood series shows, for this catchment, such omissions undermine evaluations of future risk. Our approach provides an exemplar of how to derive centennial palaeoflood reconstructions from lakes coupled well with their catchments around the world. © 2019 The Authors. Earth Surface Processes and Landforms Published by John Wiley & Sons Ltd.     

Estuarine biofilm patterns: Modern analogues for Precambrian self-organization

This field and laboratory study examines whether regularly patterned biofilms on present-day intertidal flats are equivalent to microbially induced bedforms found in geological records dating back to the onset of life on Earth. Algal mats of filamentous Vaucheria species, functionally similar to microbial biofilms, cover the topographic highs of regularly spaced ridge–runnel bedforms. As regular patterning is typically associated with self-organization processes, indicators of self-organization are tested and found to support this hypothesis. The measurements suggest that biofilm-induced sediment trapping and biostabilization enhance bedform relief, strength and multi-year persistence. This demonstrates the importance of primitive organisms for sedimentary landscape development. Algal-covered ridges consist of wavy-crinkly laminated sedimentary deposits that resemble the layered structure of fossil stromatolites and microbially induced sedimentary structures. In addition to layering, both the morphological pattern and the suggested formation mechanism of the recent bedforms are strikingly similar to microbialite strata found in rock records from the Precambrian onwards. This implies that self-organization was an important morphological process in times when biofilms were the predominant sessile ecosystem. These findings furthermore emphasize that self-organization dynamics, such as critical transitions invoking ecosystem emergence or collapse, might have been captured in fossil microbialites, influencing their laminae. This notion may be important for paleoenvironmental reconstructions based on such strata. © 2019 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd     

Improving bank erosion modelling at catchment scale by incorporating temporal and spatial variability

Bank erosion can contribute a significant portion of the sediment budget within temperate catchments, yet few catchment scale models include an explicit representation of bank erosion processes. Furthermore, representation is often simplistic resulting in an inability to capture realistic spatial and temporal variability in simulated bank erosion. In this study, the sediment component of the catchment scale model SHETRAN is developed to incorporate key factors influencing the spatio‐temporal rate of bank erosion, due to the effects of channel sinuosity and channel bank vegetation. The model is applied to the Eden catchment, north‐west England, and validated using data derived from a GIS methodology. The developed model simulates magnitudes of total catchment annual bank erosion (617–4063 t y^‐1) within the range of observed values (211–4426 t yr^‐1). In addition, the model provides both greater inter‐annual and spatial variability of bank eroded sediment generation when compared with the basic model, and indicates a potential 61% increase of bank eroded sediment as a result of temporal flood clustering. The approach developed within this study can be used within a number of distributed hydrologic models and has general applicability to temperate catchments, yet further development of model representation of bank erosion processes is required. © 2017 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Field techniques for measuring bedrock erosion and denudation

Bedrock erosion rates in natural landscapes are usually slow, of the order of millimeters per year or less, and sophisticated techniques have been developed to measure them. Different techniques have proved to be valuable depending on the spatial and temporal scale on which information is needed, on the environment and on the scientific question that is asked. Here, we give an overview of the various methods that have been developed. We introduce their working principles and outline their advantages and disadvantages. Further, we provide comprehensive references to relevant literature, both on the methods and on scientific examples of their application. © 2016 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Patch‐scale representation of vegetation within hydraulic models

Submerged aquatic vegetation affects flow, sediment and ecological processes within rivers. Quantifying these effects is key to effective river management. Despite a wealth of research into vegetated flows, the detailed flow characteristics around real plants in natural channels are still poorly understood. Here we present a new methodology for representing vegetation patches within computational fluid dynamics (CFD) models of vegetated channels. Vegetation is represented using a Mass Flux Scaling Algorithm (MFSA) and drag term within the Reynolds‐averaged Navier–Stokes Equations, which account for the mass and momentum effects of the vegetation, respectively. The model is applied using three different grid resolutions (0.2, 0.1 and 0.05 m) using time‐averaged solution methods and compared to field data. The results show that the model reproduces the complex spatial flow heterogeneity within the channel and that increasing the resolution leads to enhanced model accuracy. Future applications of the model to the prediction of channel roughness, sedimentation and key eco‐hydraulic variables are presented, likely to be valuable for informing effective river management. © 2016 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Comparing the transitional behaviour of kaolinite and bentonite suspension flows

Past research has demonstrated the dramatic effects that variations in suspended clay can have on the properties of flow by producing a range of transitional flows between turbulent and laminar states, depending on clay concentration and fluid shear. Past studies have been restricted to kaolinite flows, a clay mineral that has relatively weak cohesive properties. This paper extends these studies to suspension flows of bentonite, a clay mineral that attains higher viscosities at far lower volumetric concentrations within a flow. The results show that the types of transitional flow behaviour recognized in past studies can also be found in bentonite suspension flows, but at lower suspended sediment concentrations, thus demonstrating an even more dramatic effect on flow properties, and potentially on sediment transport and resulting bed morphology, than kaolinite flows. The paper proposes new stability diagrams for the phase space of bentonite flows and compares these to past work on kaolinite suspension flows. These new data suggest that the transitional‐flow Reynolds number can be used to delineate the types of transitional flow across different clay types and assess modern and ancient clay‐suspension flows. © 2016 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Controls on the lateral channel-migration rate of braided channel systems in coarse non-cohesive sediment

Lateral movements of alluvial river channels control the extent and reworking rates of alluvial fans, floodplains, deltas, and alluvial sections of bedrock rivers. These lateral movements can occur by gradual channel migration or by sudden changes in channel position (avulsions). Whereas models exist for rates of river avulsion, we lack a detailed understanding of the rates of lateral channel migration on the scale of a channel belt. In a two-step process, we develop here an expression for the lateral migration rate of braided channel systems in coarse, non-cohesive sediment. On the basis of photographic and topographic data from laboratory experiments of braided channels performed under constant external boundary conditions, we first explore the impact of autogenic variations of the channel-system geometry (i.e. channel-bank heights, water depths, channel-system width, and channel slope) on channel-migration rates. In agreement with theoretical expectations, we find that, under such constant boundary conditions, the laterally reworked volume of sediment is constant and lateral channel-migration rates scale inversely with the channel-bank height. Furthermore, when channel-bank heights are accounted for, lateral migration rates are independent of the remaining channel geometry parameters. These constraints allow us, in a second step, to derive two alternative expressions for lateral channel-migration rates under different boundary conditions using dimensional analysis. Fits of a compilation of laboratory experiments to these expressions suggest that, for a given channel bank-height, migration rates are strongly sensitive to water discharges and more weakly sensitive to sediment discharges. In addition, external perturbations, such as changes in sediment and water discharges or base level fall, can indirectly affect lateral channel-migration rates by modulating channel-bank heights. © 2019 The Author. Earth Surface Processes and Landforms published by John Wiley & Sons, Ltd. © 2019 The Author. Earth Surface Processes and Landforms published by John Wiley & Sons, Ltd.     

The zoogeomorphology of case-building caddisfly: Quantifying sediment use

Caddisfly (Trichoptera) larvae are an abundant and widespread aquatic insect group characterized by the construction of silk structures, including nets and cases. Case-building caddisfly have the potential to modify the sorting and mobility of sand and fine gravel via: (1) case construction, resulting in altered sediment properties; (2) transporting sediment incorporated into cases over the river bed; and (3) changing the structure of river beds via burrowing activity. To investigate these mechanisms, it is necessary to understand the mass, size distribution and spatial variability of sediment use by case-building caddisfly larvae. We quantified the mineral sediment used by individuals and communities of case-building caddisfly in 27 samples, from three sites on a gravel-bed stream. The mass and size distribution of sediment in individual cases varied between taxa (mass = 0.001–0.83 g, D₅₀ = 0.17–4 mm). The mean mass of sediment used by the caddisfly community was 38 g m⁻² and varied locally. Sediment use was predominantly coarse sand (D₅₀ = 1 mm). 64% of sediment use was attributable to Agapetus fuscipes (Glossosomatidae). Due to within-species variability in case mass, the abundance of most taxa, including A. fuscipes, was only weakly associated with the mass of sediment used by this species, at the river scale. Whilst the caddisfly community used a small percentage of the total sediment available (average 2.99% of the 1–1.4 mm size fraction), A. fuscipes used more fine sediment in their cases at sites where it was more available. Despite variability in local habitat, all sites supported diverse case-building caddisfly communities utilizing mineral sediment. Consequently, geomorphological effects of case-building caddisfly are potentially widespread. The results provide novel insights into the specific grain sizes and quantities of fine sediment used by caddisfly larvae, which represents an important step towards understanding their zoogeomorphic activities. © 2019 The Authors. Earth Surface Processes and Landforms Published by John Wiley & Sons Ltd.     

Velocity profiles and basal stresses in natural debris flows

The internal deformation within debris flows holds essential information on dynamics and flow resistance of such mass-wasting processes. Systematic measurements of velocity profiles in real-scale debris flows are not yet available. Additionally, data on basal stresses of the solid and the fluid phase are rare. Here, we present and analyse measurements of vertical velocity profiles in two debris flows naturally occurring in the Gadria Creek, Italy. The method is based on cross-correlation of paired conductivity signals from an array of sensors installed on a fin-shaped wall located in the middle of the channel. Additionally, we measure normal stress and pore fluid pressure by two force plates with integrated pressure transducers. We find internal deformation throughout the flows. Only at the very front was some en-bloc movement observed. Velocity profiles varied from front to tail and between flows. For one debris flow, pore fluid pressure close to normal stress was measured, whereas the other flow was less liquefied. The median shear rates were mostly less than 5 s⁻¹ and Savage numbers at the basal layer ranged from 0.01 to 1. Our results highlight the variable nature of debris flows and provide quantitative data on shear rate and basal stress distribution to help guide model development for hazard assessment and landscape evolution. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Deep-seated gravitational slope deformation (DSGSD) and slow-moving landslides in the southern Tien Shan Mountains: new insights from InSAR,tectonic and geomorphic analysis

We investigated deep-seated gravitational slope deformation (DSGSD) and slow mass movements in the southern Tien Shan Mountains front using synthetic aperture radar (SAR) time-series data obtained by the ALOS/PALSAR satellite. DSGSD evolves with a variety of geomorphological changes (e.g. valley erosion, incision of slope drainage networks) over time that affect earth surfaces and, therefore, often remain unexplored. We analysed 118 interferograms generated from 20 SAR images that covered about 900 km². To understand the spatial pattern of the slope movements and to identify triggering parameters, we correlated surface dynamics with the tectono-geomorphic processes and lithologic conditions of the active front of the Alai Range. We observed spatially continuous, constant hillslope movements with a downslope speed of approximately 71 mm year⁻¹ velocity. Our findings suggest that the lithological and structural framework defined by protracted deformation was the main controlling factor for sustained relief and, consequently, downslope mass movements. The analysed structures revealed integration of a geological/structural setting with the superposition of Cretaceous–Paleogene alternating carbonatic and clastic sedimentary structures as the substratum for younger, less consolidated sediments. This type of structural setting causes the development of large-scale, gravity-driven DSGSD and slow mass movement. Surface deformations with clear scarps and multiple crest lines triggered planes for large-scale deep mass creeps, and these were related directly to active faults and folds in the geologic structures. Our study offers a new combination of InSAR techniques and structural field observations, along with morphometric and seismologic correlations, to identify and quantify slope instability phenomena along a tectonically active mountain front. These results contribute to an improved natural risk assessment in these structures. © 2019 The Authors. Earth Surface Processes and Landforms Published by John Wiley & Sons Ltd     

Tools for gauging the capacity of salmon spawning substrates

We present a set of river management tools based on a recently developed method for estimating the amount of salmon spawning habitat in coarse‐bedded rivers. The method, which was developed from a mechanistic model of redd building by female salmon, combines empirical relationships between fish length, redd area, and the sizes of particles moved by fish during spawning. Model inputs are the grain‐size indices D₅₀ and D₈₄ and an estimate of female fish length, which is used to predict the size of the redd that they will build and the size of the largest particle that they can move on the bed. Outputs include predictions of the fraction of the bed that the fish can use for redd building and the number of redds that they can build within the useable area. We cast the model into easy‐to‐use look‐up tables, charts, an Excel worksheet, a JavaScript web applet, and a MATLAB user interface. We explain how these tools can be used in a new, mechanistic approach to assessing spawning substrates and optimizing gravel augmentation projects in coarse‐bedded rivers. © 2016 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Reconstructing ice-flow fields from streamlined subglacial bedforms: A kriging approach

The orientation of several landforms, e.g. drumlins, flutes, crag-and-tails, and mega-scale glacial lineations, records the direction of the overlying ice flow that created them. Populations of such features are used routinely to infer former ice-flow patterns, which serve as the building blocks of reconstructions of palaeo ice-sheet evolution. Currently, the conceptualisation of flow patterns from these flow-direction records is done manually and qualitatively, so the extractable glaciological information is limited. We describe a kriging method (with Matlab code implementation) that calculates continuous fields of ice-flow direction, convergence, and curvature from the flow-direction records, and which yields quantitative results with uncertainty estimates. We test the method by application to the subglacial bedforms of the Tweed Valley Basin, UK. The results quantify the convergent flow pattern of the Tweed Palaeo-Ice Stream in detail and pinpoint its former lateral shear margins and where ice flowed around basal bumps. Ice-flow parameters retrieved by this method can enrich ice-sheet reconstructions and investigations of subglacial till processes and bedform genesis. © 2018 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Bedforms: views and new perspectives from the third international workshop on Marine and River Dune Dynamics (MARID3)

This Commentary introduces a Special Issue of Earth Surface Processes and Landforms dedicated to the memory of Stephen Coleman (1966–2012) who sadly passed away on July 23, 2012. Stephen was an incredible enthusiast for the study of bedforms and made seminal contributions to our knowledge and understanding of flow, sediment transport and bedform dynamics. Stephen presented a keynote address at MARID 3, which has led to this Special Issue (see Coleman SE, Nikora VI. 2011. Fluvial dunes: initiation, characterization, flow structure. Earth Surface Processes and Landforms 36 : 39–57. DOI. 10.1002/esp.2096) and was always keen to discuss research, from theory to techniques to experiments. He was always generous with his time and touched many with the influence of his thoughts and intellect. Our community has lost an outstanding scientist and Stephen will be greatly missed by us all. Copyright © 2012 John Wiley & Sons, Ltd.