First cosmogenic 10Be age constraint on the timing of Younger Dryas glaciation and ice cap thickness,western Scottish Highlands

We use cosmogenic ¹⁰Be surface exposure age techniques at a locality close to Rannoch Moor, western Scottish Highlands, in order to establish the age and chronology of its most recent glaciation. Glacial erratics and an in situ bedrock quartz vein sampled from this site—the summit of Beinn Inverveigh—have yielded zero‐erosion exposure ages of 12.9 ± 1.5 ka to 11.6 ± 1.0 ka, implying complete ice cover of the mountain during the Younger Dryas, or Loch Lomond Stadial. These results fit closely with published ¹⁴C dates that bracket the maximum (lateral) extent of ice cap outlet glaciers, and are the first internally consistent ages to specifically address this period of glaciation in Scotland. Furthermore, the dates imply that previous palaeoglaciological reconstructions for this area may have underestimated both the thickness of the former ice cap and, by implication, its volume. © British Geological Survey/Natural Environment Research Council copyright 2007. Reproduced with the permission of BGS/NERC. Published by John Wiley & Sons, Ltd.     

Statistical biases and errors inherent in photoclinometric surface slope estimation with natural light

Photoclinometry is the most common method used to obtain high-resolution topographic maps of planetary terrain. We derive the likelihood function of photoclinometric surface slope from (1) the probability distribution of the measured photon count of natural sunlight through a Charge-Coupled Device (CCD) including uncertainty due to camera shot noise, camera read noise, small-scale albedo fluctuation and atmospheric haze, and (2) a photometric model relating photocount to surface orientation. We then use classical estimation theory to determine the theoretically exact biases and errors inherent in photoclinometric surface slope and show when they may be approximated by asymptotic expressions for sufficiently high sample size. We show how small-scale albedo variability often dominates biases and errors, which may become an order of magnitude larger than surface slopes when surface reflectance has a weak dependence on surface tilt. We provide bounds on the minimum possible error of any unbiased photoclinometric surface slope estimate, and compute the sample sizes necessary to constrain errors within desired design thresholds.     

Scales and causes of heterogeneity in bars in a large multi‐channel river: Río Paraná, Argentina

To date, published studies of alluvial bar architecture in large rivers have been restricted mostly to case studies of individual bars and single locations. Relatively little is known about how the depositional processes and sedimentary architecture of kilometre‐scale bars vary within a multi‐kilometre reach or over several hundreds of kilometres downstream. This study presents Ground Penetrating Radar and core data from 11, kilometre‐scale bars from the Río Paraná, Argentina. The investigated bars are located between 30 km upstream and 540 km downstream of the Río Paraná – Río Paraguay confluence, where a significant volume of fine‐grained suspended sediment is introduced into the network. Bar‐scale cross‐stratified sets, with lengths and widths up to 600 m and thicknesses up to 12 m, enable the distinction of large river deposits from stacked deposits of smaller rivers, but are only present in half the surface area of the bars. Up to 90% of bar‐scale sets are found on top of finer‐grained ripple‐laminated bar‐trough deposits. Bar‐scale sets make up as much as 58% of the volume of the deposits in small, incipient mid‐channel bars, but this proportion decreases significantly with increasing age and size of the bars. Contrary to what might be expected, a significant proportion of the sedimentary structures found in the Río Paraná is similar in scale to those found in much smaller rivers. In other words, large river deposits are not always characterized by big structures that allow a simple interpretation of river scale. However, the large scale of the depositional units in big rivers causes small‐scale structures, such as ripple sets, to be grouped into thicker cosets, which indicate river scale even when no obvious large‐scale sets are present. The results also show that the composition of bars differs between the studied reaches upstream and downstream of the confluence with the Río Paraguay. Relative to other controls on downstream fining, the tributary input of fine‐grained suspended material from the Río Paraguay causes a marked change in the composition of the bar deposits. Compared to the upstream reaches, the sedimentary architecture of the downstream reaches in the top ca 5 m of mid‐channel bars shows: (i) an increase in the abundance and thickness (up to metre‐scale) of laterally extensive (hundreds of metres) fine‐grained layers; (ii) an increase in the percentage of deposits comprised of ripple sets (to >40% in the upper bar deposits); and (iii) an increase in bar‐trough deposits and a corresponding decrease in bar‐scale cross‐strata (<10%). The thalweg deposits of the Río Paraná are composed of dune sets, even directly downstream from the Río Paraguay where the upper channel deposits are dominantly fine‐grained. Thus, the change in sedimentary facies due to a tributary point‐source of fine‐grained sediment is primarily expressed in the composition of the upper bar deposits.     

Quantifying flow resistance in mountain streams using computational fluid dynamics modeling over structure-from-motion photogrammetry-derived microtopography

Flow resistance in mountain streams is important for assessing flooding hazard and quantifying sediment transport and bedrock incision in upland landscapes. In such settings, flow resistance is sensitive to grain-scale roughness, which has traditionally been characterized by particle size distributions derived from laborious point counts of streambed sediment. Developing a general framework for rapid quantification of resistance in mountain streams is still a challenge. Here we present a semi-automated workflow that combines millimeter- to centimeter-scale structure-from-motion (SfM) photogrammetry surveys of bed topography and computational fluid dynamics (CFD) simulations to better evaluate surface roughness and rapidly quantify flow resistance in mountain streams. The workflow was applied to three field sites of gravel, cobble, and boulder-bedded channels with a wide range of grain size, sorting, and shape. Large-eddy simulations with body-fitted meshes generated from SfM photogrammetry-derived surfaces were performed to quantify flow resistance. The analysis of bed microtopography using a second-order structure function identified three scaling regimes that corresponded to important roughness length scales and surface complexity contributing to flow resistance. The standard deviation σ_z of detrended streambed elevation normalized by water depth, as a proxy for the vertical roughness length scale, emerges as the primary control on flow resistance and is furthermore tied to the characteristic length scale of rough surface-generated vortices. Horizontal length scales and surface complexity are secondary controls on flow resistance. A new resistance predictor linking water depth and vertical roughness scale, i.e.  H/σ_z, is proposed based on the comparison between σ_z and the characteristic length scale of vortex shedding. In addition, representing streambeds using digital elevation models (DEM) is appropriate for well-sorted streambeds, but not for poorly sorted ones under shallow and medium flow depth conditions due to the missing local overhanging features captured by fully 3D meshes which modulate local pressure gradient and thus bulk flow separation and pressure distribution. An appraisal of the mesh resolution effect on flow resistance shows that the SfM photogrammetry data resolution and the optimal CFD mesh size should be about 1/7 to 1/14 of the standard deviation of bed elevation. © 2019 John Wiley & Sons, Ltd.     

From taxonomies to ontologies: formalizing generalization knowledge for on-demand mapping

Automation of the cartographic design process is central to the delivery of bespoke maps via the web. In this paper, ontological modeling is used to explicitly represent and articulate the knowledge used in this decision-making process. A use case focuses on the visualization of road traffic accident data as a way of illustrating how ontologies provide a framework by which salient and contextual information can be integrated in a meaningful manner. Such systems are in anticipation of web-based services in which the user knows what they need, but do not have the cartographic ability to get what they want.     

What is a flux tube? On the magnetic field topology of buoyant flux structures

We study the topology of field lines threading buoyant magnetic flux structures. The magnetic structures, visually resembling idealized magnetic flux tubes, are generated self-consistently by numerical simulation of the interaction of magnetic buoyancy and a localized velocity shear in a stably stratified atmosphere. Depending on the parameters, the system exhibits varying degrees of symmetry. By integrating along magnetic field lines and constructing return maps, we show that, depending on the type of underlying behaviour, the stages of the evolution, and therefore the degree of symmetry, the resulting magnetic structures can have field lines with one of three distinct topologies. When the x -translational and y -reflectional symmetries remain intact, magnetic field lines lie on surfaces but individual lines do not cover the surface. When the y symmetry is broken, magnetic field lines lie on surfaces and individual lines do cover the surface. When both x and y symmetries are broken, magnetic field lines wander chaotically over a large volume of the magnetically active region. We discuss how these results impact our simple ideas of a magnetic flux tube as an object with an inside and an outside, and introduce the concept of 'leaky' tubes.