Recent progress on combining geomorphological and geochronological data with ice sheet modelling,demonstrated using the last British–Irish Ice Sheet

Palaeo-ice sheets are important analogues for understanding contemporary ice sheets, offering a record of ice sheet behaviour that spans millennia. There are two main approaches to reconstructing palaeo-ice sheets. Empirical reconstructions use the available glacial geological and chronological evidence to estimate ice sheet extent and dynamics but lack direct consideration of ice physics. In contrast, numerically modelled simulations implement ice physics, but often lack direct quantitative comparison with empirical evidence. Despite being long identified as a fruitful scientific endeavour, few ice sheet reconstructions attempt to reconcile the empirical and model-based approaches. To achieve this goal, model-data comparison procedures are required. Here, we compare three numerically modelled simulations of the former British–Irish Ice Sheet with the following lines of evidence: (a) position and shape of former margin positions, recorded by moraines; (b) former ice-flow direction and flow-switching, recorded by flowsets of subglacial bedforms; and (c) the timing of ice-free conditions, recorded by geochronological data. These model–data comparisons provide a useful framework for quantifying the degree of fit between numerical model simulations and empirical constraints. Such tools are vital for reconciling numerical modelling and empirical evidence, the combination of which will lead to more robust palaeo-ice sheet reconstructions with greater explicative and ultimately predictive power.     

Linkages between spatio‐temporal patterns of environmental factors and distribution of plant assemblages across a boreal peatland complex

Here we examine the arrangement of plant species across an oligotrophic bog/poor fen peatland complex in the North American boreal plain and the relationships of these species to their physical and chemical environment. A semi‐uniform spatial sampling approach was utilized to describe the species assemblages, pore‐water chemistry and physical condition of 100 plots throughout a single peatland complex. Regardless of sharing the same ground cover of Sphagnum mosses, the remaining species separated into four distinct assemblages, each with unique indicators. These species groups along with associated chemical and physical factors are organized into four ecosites: bog, margin (edge) and two poor fen ecosites. The plant assemblages of this peatland have a complex relationship with numerous gradients, both physical and chemical, including depth to water table, shade, pH, nutrient and base cation. Rather than being homogenous across the landscape, most environmental variables exhibit distinct spatial patterns and do so in relationship to the plant assemblages, forming spatially distinct ecosites across the complex. Base cation concentrations play a smaller role than previously thought in differentiating these ecosites, and in addition to shade and depth to water table, nitrogen in the form of dissolved organic nitrogen was highly related to the placement of these ecosites. Many significant chemical factors appear related to evaporative water loss within the peatland complex, and these chemical factors are used to differentiate the ecosites. However, the mediation of evaporative water loss is due largely to self‐generated responses of the plant assemblages related to shade through plant morphology and peat acrotelm development related to depth to water table. We conclude that plant species and associated environmental gradients act together to form spatially distinct ecosites. The distribution of these ecosites within this large, environmentally complex peatland is largely controlled by differing self‐generated responses along the hydrotopographical gradient of differential water loss.     

Variability in bedform morphology and kinematics with transport stage

Bedform geometry is widely recognized to be a function of transport stage. Bedform aspect ratio (height/length) increases with transport stage, reaches a maximum, then decreases as bedforms washout to a plane bed. Bedform migration rates are also linked to bedform geometry, in so far as smaller bedforms in coarser sediment tend to migrate faster than larger bedforms in finer sediment. However, how bedform morphology (height, length and shape) and kinematics (translation and deformation) change with transport stage and suspension have not been examined. A series of experiments is presented where initial flow depth and grain size were held constant and the transport stage was varied to produce bedload dominated, mixed‐load dominated and suspended‐load dominated conditions. The results show that the commonly observed pattern in bedform aspect ratio occurs because bedform height increases then decreases with transport stage, against a continuously increasing bedform length. Bedform size variability increased with transport stage, leading to less uniform bedform fields at higher transport stage. Total translation‐related and deformation‐related sediment fluxes all increased with transport stage. However, the relative contribution to the total flux changed. At the bedload dominated stage, translation‐related and deformation‐related flux contributed equally to the total flux. As the transport stage increased, the fraction of the total load contributed by translation increased and the fraction contributed by deformation declined because the bedforms got bigger and moved faster. At the suspended‐load dominated transport stage, the deformation flux increased and the translation flux decreased as a fraction of the total load, approaching one and zero, respectively, as bedforms washed out to a plane bed.     

Response of low‐angle dunes to variable flow

Current understanding of bedform dynamics is largely based on field and laboratory observations of bedforms in steady flow environments. There are relatively few investigations of bedforms in flows dominated by unsteadiness associated with rapidly changing flows or tides. As a consequence, the ability to predict bedform response to variable flow is rudimentary. Using high‐resolution multibeam bathymetric data, this study explores the dynamics of a dune field developed by tidally modulated, fluvially dominated flow in the Fraser River Estuary, British Columbia, Canada. The dunes were dominantly low lee angle features characteristic of large, deep river channels. Data were collected over a field ca 1·0 km long and 0·5 km wide through a complete diurnal tidal cycle during the rising limb of the hydrograph immediately prior to peak freshet, yielding the most comprehensive characterization of low‐angle dunes ever reported. The data show that bedform height and lee angle slope respond to variable flow by declining as the tide ebbs, then increasing as the tide rises and the flow velocities decrease. Bedform lengths do not appear to respond to the changes in velocity caused by the tides. Changes in the bedform height and lee angle have a counterclockwise hysteresis with mean flow velocity, indicating that changes in the bedform geometry lag changes in the flow. The data reveal that lee angle slope responds directly to suspended sediment concentration, supporting previous speculation that low‐angle dune morphology is maintained by erosion of the dune stoss and crest at high flow, and deposition of that material in the dune trough.     

Field investigation of deformation characteristics and stress mobilisation of a soil slope

Stress mobilisation and deformation of a slope are important for engineers to carry out reliable design of retaining systems. However, most case histories reported mainly on the response of pore water pressure (PWP), whereas knowledge about the stress deformation characteristics of slope is limited. In this study, a saprolitic soil slope was instrumented to monitor not only the responses of PWP but also horizontal stress and horizontal displacement. To assist in the interpretation of field data, a series of laboratory tests was conducted to characterise volume change behaviour of the soil taken from the site, under the effects of both net stress and suction. During a rainstorm event when positive PWP built up, a remarkably large displacement of 20 mm was recorded between 5.5- and 6-m depths, and the top 5 m of the slope exhibited translational downslope movement. This caused an increase in Bishop’s effective horizontal stress by 350 %, which reached a peak value close to 40 % of a Bishop’s effective passive stress. During the subsequent dry season when suction was recovered, an upslope rebound of 10 mm was recorded. Comparison of field and laboratory data reveals that the rebound was attributed to suction-induced soil shrinkage. This rebound led to a decrease in the Bishop’s effective horizontal stress previously built up during the storm event.     

Uncertainty analysis of space–time prisms based on the moment-design method

Space–time prism (STP) is an important concept for the modeling of object movements in space and time. An STP can be conceptualized as the result of the potential path of a moving object revolving around in the three-dimensional space. Though the concept has found applications in time geography, research on the analysis and propagation of uncertainty in STPs, particularly under high degree of nonlinearity, is scanty. Based on the efficiency and effectiveness of the moment-design (M-D) method, this paper proposes an approach to deal with nonlinear error propagation problems in the potential path areas (PPAs) of STPs and their intersections. Propagation of errors to the PPA and its boundary, and to the intersection of two PPAs is investigated. Performance of the proposed method is evaluated via a series of experimental studies. In comparison with the Monte Carlo method and the implicit function method, simulation results show the advantages of the M-D method in the analysis of error propagation in STPs.     

Dasymetric Spatiotemporal Interpolation

This research applies the principles of dasymetric mapping to spatiotemporal interpolation by extending the spatial concepts of zone and area to their temporal analogs of interval and duration, respectively. An example application of dasymetric spatiotemporal interpolation using crime event data is presented. Results indicate that dasymetric spatiotemporal interpolation significantly improves the accuracy of estimates over areal or duration weighting. In addition, even when dasymetric interpolation in either the spatial or temporal dimension is relatively weak, combining dasymetric estimation in both space and time dimensions simultaneously has the potential to amplify the accuracy of the overall dasymetric estimation.