Residual relief separation: digital elevation model enhancement for geomorphological mapping

Geomorphologically mapped data form a primary set of observations that can be used to infer former environmental conditions. Thus, objective and consistent mapping of landforms from remotely sensed data (e.g. satellite imagery, digital elevation models (DEMs)) is paramount for reconstructing palaeoenvironments. This paper proposes a technique, ‘residual relief separation’, to enhance landforms in DEMs prior to visualization and digital mapping. This is applied to a ～600 km² region surrounding Lough Gara, Ireland, where drumlins (～200 m wide) overlie a regional relief of hills (～10 km wide). Here, residual relief separation uses this difference in width‐scale. Regional relief is approximated by a 1 km wide median filter, then subtracted to leave the drumlins in a ‘residual’ topography. In a second step, the residual relief is normalized to allow for amplitude variations in the drumlins across the area (～5–40 m high). Finally, visualization uses a simple black‐to‐white colour scale for height. Although not numerically outperforming other visualization techniques, this method performs equally well, and as the data are not ‘illuminated’ there is no azimuthal bias. Additional benefits include the relatively simple calculation, intuitive visual comprehension, no emphasis of noise, and the possibility of using any desired visualization technique after the landscape has been topographically manipulated. Copyright © 2008 John Wiley & Sons, Ltd.     

Hyperscale terrain modelling of braided rivers: fusing mobile terrestrial laser scanning and optical bathymetric mapping

Quantifying the morphology of braided rivers is a key task for understanding braided river behaviour. In the last decade, developments in geomatics technologies and associated data processing methods have transformed the production of precise, reach‐scale topographic datasets. Nevertheless, generating accurate Digital Elevation Models (DEMs) remains a demanding task, particularly in fluvial systems. This paper identifies a threefold set of challenges associated with surveying these dynamic landforms: complex relief, inundated shallow channels and high rates of sediment transport, and terms these challenges the ‘morphological’, ‘wetted channel’ and ‘mobility’ problems, respectively. In an attempt to confront these issues directly, this paper presents a novel survey methodology that combines mobile terrestrial laser scanning and non‐metric aerial photography with data reduction and surface modelling techniques to render DEMs from the resulting very high resolution datasets. The approach is used to generate and model a precise, dense topographic dataset for a 2.5 km reach of the braided Rees River, New Zealand. Data were acquired rapidly between high flow events and incorporate over 5 x 10⁹ raw survey observations with point densities of 1600 pts m^‐2 on exposed bar and channel surfaces. A detailed error analysis of the resulting sub‐metre resolution is described to quantify DEM quality across the entire surface model. This reveals unparalleled low vertical errors for such a large and complex surface model; between 0.03 and 0.12 m in exposed and inundated areas of the model, respectively. Copyright © 2013 John Wiley & Sons, Ltd.     

Application of globally available,coarse-resolution digital elevation models for delineating valley bottom segments of varying length across a catchment

The worldwide availability of digital elevation models (DEMs) has enabled rapid (semi-)automated mapping of earth surface landforms. In this paper, we first present an approach for delineating valley bottom extent across a large catchment using only publicly available, coarse-resolution DEM input. We assess the sensitivity of our results to variable DEM resolution and find that coarse-resolution datasets (90 m resolution) provide superior results. We also find that LiDAR-derived DEMs produce more realistic results than satellite-derived DEMs across the full range of topographic settings tested. Satellite-derived DEMs perform more effectively in moderate topographic settings, but fail to capture the subtleties of valley bottom extent in mild gradient, low-lying topography and in narrow headwater reaches. Second, we present a semi-automated technique within ArcGIS for delineating valley bottom segments using DEM-derived network scale metrics of valley bottom width and slope. We use an unsupervised machine-learning technique based on the k-means clustering algorithm to solve a conundrum in GIS-based geomorphic analysis of rivers: the delineation of valley bottom segments of variable length. The delineation of valley bottom segments provides a coarse-scale entry point into automated geomorphic analysis and characterization of river systems. © 2020 John Wiley & Sons, Ltd.     

Geomorphic response to an active transpressive regime: a case study along the Chaman strike‐slip fault,western Pakistan

The Chaman left‐lateral strike‐slip fault bounds the rigid Indian plate boundary at the western end of the Himalayan‐Tibetan orogen and is marked by contrasting topographic relief. Deformed landforms along the fault provide an excellent record for understanding this actively evolving intra‐continental strike‐slip fault. The geomorphic response of an active transpessional stretch of the Chaman fault was studied using digital elevation model (DEM) data integrated with Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) Visible and Near Infrared/Short Wave Infrared (VNIR/SWIR) and images from GeoEye‐1. Geologic and geomorphic mapping helped in reconstructing the Late Quaternary landscape history of this transpessional strand of the Chaman strike‐slip fault and the associated Spinatizha thrust fault in western Pakistan. Topographic analysis of a part of the transpression (the thrust bounded Roghani ridge) revealed northward growth of the Spinatizha fault with the presence of three water gaps and two corresponding wind gaps. Geomorphic indices including stream length‐gradient index, mountain front sinuosity, valley floor width to valley height ratios, and entrenchment of recent alluvial fan deposits were used to define the lateral growth and direction of propagation of the Spinatizha fault. Left‐lateral displacement along Chaman fault and uplift along the Spinatizha fault was defined using topographic analysis of the Roghani ridge and geomorphic mapping of an impressive alluvial fan, the Bostankaul fan. The landforms and structures record slip partitioning along the Indian plate boundary, and account for the convergence resulting from the difference in the Chaman fault azimuth and orientation of the velocity vector of the Indian plate. Copyright © 2012 John Wiley & Sons, Ltd.     

Flow network derivation from a high resolution DEM in a low relief,agrarian landscape

Digital flow networks derived from digital elevation models (DEMs) sensitively react to errors due to measurement, data processing and data representation. Since high‐resolution DEMs are increasingly used in geomorphological and hydrological research, automated and semi‐automated procedures to reduce the impact of such errors on flow networks are required. One such technique is stream‐carving, a hydrological conditioning technique to ensure drainage connectivity in DEMs towards the DEM edges. Here we test and modify a state‐of‐the‐art carving algorithm for flow network derivation in a low‐relief, agricultural landscape characterized by a large number of spurious, topographic depressions. Our results show that the investigated algorithm reconstructs a benchmark network insufficiently in terms of carving energy, distance and a topological network measure. The modification to the algorithm that performed best, combines the least‐cost auxiliary topography (LCAT) carving with a constrained breaching algorithm that explicitly takes automatically identified channel locations into account. We applied our methods to a low relief landscape, but the results can be transferred to flow network derivation of DEMs in moderate to mountainous relief in situations where the valley bottom is broad and flat and precise derivations of the flow networks are needed. Copyright © 2013 John Wiley & Sons, Ltd.     

Comparison of geomorphological field mapping and 2D‐InSAR mapping of periglacial landscape activity at Nordnesfjellet,northern Norway

The ability to continuously monitor the dynamic response of periglacial landforms in a climate change context is of increasing scientific interest. Satellite radar interferometry provides information on surface displacement that can be related to periglacial processes. Here we present a comparison of two‐dimensional (2D) surface displacement rates and geomorphological mapping at periglacial landform and sediment scale from the mountain Nordnesfjellet in northern Norway. Hence, 2D Interferometric Synthetic Aperture Radar (InSAR) results stem from a 2009–2014 TerraSAR‐X dataset from ascending and descending orbits, decomposed into horizontal displacement vectors along an east–west plane, vertical displacement vectors and combined displacement velocity. Geomorphological mapping was carried out on aerial imagery and validated in the field. This detailed landform and sediment type mapping revealed an altitudinal distribution dominated by, weathered bedrock blockfields, surrounded primarily by slightly, to non‐vegetated solifluction landforms at the mountain tops. Below, an active rockslide and associated rockfall deposits are located on the steep east‐facing side of the study area, whereas glacial sediments dominate on the gentler western side. We show that 2D InSAR correctly depicts displacement rates that can be associated with typical deformation patterns for flat‐lying or inclined landforms, within and below the regional permafrost limit, for both wet and dry areas. A net lowering of the entire landscape caused by general denudation of the periglacial landforms and sediments is here quantified for the first time using radar remote sensing. Copyright © 2018 John Wiley & Sons, Ltd.     

The normalized topographic method: an automated procedure for gully mapping using GIS

Gully delineation is a critical aspect of accurately determining soil losses but associated methodologies are rarely detailed. Here, we describe a new gully mapping method, the normalized topographic method (NorToM), based on processing digital elevation model (DEM) data, and we assess associated errors when it is applied over a range of geomorphological scales. The NorToM is underpinned by two gully detection variables (normalized slope and elevation) calculated over local windows of prescribed size, and a group of filtering variables. For four study sites, DEMs of gullies were obtained using field and airborne photo‐reconstruction and evaluated using total station and differential global positioning system (dGPS) survey. NorToM provided accurate areal and volume estimates at the individual gully scale but differences increased at the larger gully system and gully network scales. We were able to identify optimal parameters for using the NorToM approach and so confirm that is represents a useful scale‐independent means of gully mapping that is likely to be valid in other environments. Its main limitations are that the normalization process might be time‐consuming at regional scales and the need for a fixed window size when applied to landforms with extreme variations in dimensions. Copyright © 2014 John Wiley & Sons, Ltd.     

Linking landscape morphological complexity and sediment connectivity

Connectivity relates to the coupling of landforms (e.g. hillslopes and channels) and the transfer of water and sediment between them. The degree to which parts of a catchment are connected depends largely on the morphological complexity of the catchment's landscape. Landscapes can have very different and distinct morphologies, such as terraces, V‐shaped valleys or broad floodplains. The objective of this study is to better understand and quantify the relation between landscape complexity and catchment connectivity. We hypothesize that connectivity decreases with increasing landscape morphological complexity. To quantify the connectivity–complexity relationship virtual digital elevation models (DEMs) with distinct morphologies were used as inputs into the landscape evolution model LAPSUS to simulate the sediment connectivity of each landscape. Additionally, the hypothesis was tested on six common real DEMs with widely different morphologies. Finally, the effects of different rainfall time series on catchment response were explored. Simulation results confirm the hypothesis and quantify the non‐linear relation. Results from the exploration of sediment connectivity in response to sequences of rainfall events indicate that feedback between erosion and deposition are more important for certain landscape morphologies than for others: for a given rainfall input, a more effective sediment connectivity and erosion response may be expected from rolling or V‐shaped catchments than from dissected or stepped landscapes. Awareness of the differences in the behaviour and response of different morphologies to catchment processes provides valuable information for the effective management of landscapes and ecosystems through efficiently designed soil and water conservation measures. Copyright © 2013 John Wiley & Sons, Ltd.     

Efficient hybrid breaching‐filling sink removal methods for flow path enforcement in digital elevation models

Digital elevation models (DEMs) that are used in hydrological applications must be processed to remove sinks, mainly topographic depressions. Flow enforcement techniques include filling methods, which raise elevations within depressions, breaching, which carves channels through blockages, and hybrid methods. Despite previous research demonstrating the large impact to DEMs and subsequent analyses of depression filling, it is common practice apply this technique to flow enforcement. This is partly because of the greater efficiency of depression filling tools compared to breaching counterparts, which often limits breaching to applications of small‐ to moderate‐sized DEMs. A new hybrid flow enforcement algorithm is presented in this study. The method can be run in complete breaching, selective breaching (either breached or filled), or constrained breaching (partial breaching) modes, allowing for greater flexibility in how practitioners enforce continuous flow paths. Algorithm performance was tested with DEMs of varying topography, spatial extents, and resolution. The sites included three moderate sized DEMs (52 000 000 to 190 000 000 cells) and three massive DEMs of the Iberian Peninsula, and the Amazon and Nile River basins, the largest containing nearly one billion cells. In complete breaching mode, the new algorithm required 87% of the time needed by a filling method to process the test DEMs, while the selective breaching and constrained breaching modes, operating with maximum breach depth constraints, increased run times by 8% and 27% respectively. Therefore, the new algorithm offers comparable performance to filling and the ability to process massive topographic data sets, while giving practitioners greater flexibility and lowering DEM impact. Copyright © 2015 John Wiley & Sons, Ltd.     

TerEx Toolbox for semi‐automated selection of fluvial terrace and floodplain features from lidar

Terraces and floodplains are important indicators of near‐channel sediment dynamics, serving as potential sediment sources and sinks. Increasing availability of high resolution topography data over large areas calls for development of semi‐automated techniques for identification and measurement of these features. In this study we introduce a novel tool that accommodates user‐defined parameters including, a local‐relief threshold selected by a variable‐size moving window, minimum area threshold, and maximum distance from the channel to identify and map discrete terrace and floodplain surfaces. Each of the parameters can easily be calibrated for a given watershed or reach. Subsequently, the tool automatically measures planform area, absolute elevation, and height relative to the local river channel for each terrace polygon. We validate the tool in two locations where terrace maps were previously developed via manual digitization from lidar and extensive field mapping campaigns. The tool is also tested on six different types of rivers to provide examples of starting selection parameters, and to test effectiveness of the tool across a wide range of landscapes. Generally, the tool provides a high quality draft map of terrace and floodplain surfaces across the wide range of environmental conditions for which it has been tested. We find that the tool functions best in catchments where the terraces are spatially extensive, with distinct differences between the terrace and floodplain. The most challenging environments for semi‐automated terrace and floodplain mapping include steep catchments with dense riparian vegetation, and very small terraces (~10 m² in areal extent). We then apply the tool to map terraces and floodplains in the Root River watershed, southeastern Minnesota and generate exceedance plots for terrace heights. These plots provide a first pass analysis to indicate the tributaries and reaches of the river where terraces constitute a significant source of sediment. Copyright © 2013 John Wiley & Sons, Ltd.     

Performance evaluation of three DEM‐based fluvial terrace mapping methods

The availability of high‐resolution digital elevation models (DEMs) derived from airborne light detection and ranging (LiDAR) surveys has spurred the development of several methods to identify and map fluvial terraces. The post‐glacial landscape of the Sheepscot River watershed, Maine, where land‐use change has produced fill terraces upstream of historic dam sites, was selected to implement a comparison between terrace mapping methodologies. At four study sites within the watershed, terraces were manually mapped on LiDAR‐DEM‐derived hillshade images to facilitate the comparison among fully and semi‐automated DEM‐based procedures, including: (1) spatial relationships between interpreted terraces and surrounding natural topography, (2) feature classification algorithms, and (3) the TerEx terrace mapping toolbox. Each method was evaluated based on its accuracy and ease of implementation. The four study sites have varying longitudinal slope (0.0008–0.006 m/m), channel width (< 5–30 m), surrounding landscape relief (20–80 m), type and density of surrounding land use, and mapped surficial geologic units. All methods generally overestimate terrace areas (average predicted area 210% of manually defined area) with the most accurate results achieved within confined river valleys surrounded by the steep hillslopes. Accuracy generally decreases for study sites surrounded by low‐relief landscapes (predicted areas ranged 4–953% of manual delineations). We conclude with the advantages and drawbacks of each method tested and make recommendations for the scenarios where the use of each method is most appropriate. Copyright © 2016 John Wiley & Sons, Ltd.     

The influence of elevation error on the morphometrics of channel networks extracted from DEMs and the implications for hydrological modelling

Stream network morphometrics have been used frequently in environmental applications and are embedded in several hydrological models. This is because channel network geometry partly controls the runoff response of a basin. Network indices are often measured from channels that are mapped from digital elevation models (DEMs) using automated procedures. Simulations were used in this paper to study the influence of elevation error on the reliability of estimates of several common morphometrics, including stream order, the bifurcation, length, area and slope ratios, stream magnitude, network diameter, the flood magnitude and timing parameters of the geomorphological instantaneous unit hydrograph (GIUH) and the network width function. DEMs of three UK basins, ranging from high to low relief, were used for the analyses. The findings showed that moderate elevation error (RMSE of 1·8 m) can result in significant uncertainty in DEM‐mapped network morphometrics and that this uncertainty can be expressed in complex ways. For example, estimates of the bifurcation, length and area ratios and the flood magnitude and timing parameters of the GIUH each displayed multimodal frequency distributions, i.e. two or more estimated values were highly likely. Furthermore, these preferential estimates were wide ranging relative to the ranges typically observed for these indices. The wide‐ranging estimates of the two GIUH parameters represented significant uncertainty in the shape of the unit hydrograph. Stream magnitude, network diameter and the network width function were found to be highly sensitive to elevation error because of the difficulty in mapping low‐magnitude links. Uncertainties in the width function were found to increase with distance from outlet, implying that hydrological models that use network width contain greater uncertainty in the shape of the falling limb of the hydrograph. In light of these findings, care should be exercised when interpreting the results of analyses based on DEM‐mapped stream networks. Copyright © 2007 John Wiley & Sons, Ltd.     

Paleotopographic controls on loess deposition in the Loess Plateau of China

The underlying pre‐existing paleotopography directly influences the loess deposition process and shapes the morphology of current loess landforms. An understanding of the controlling effects of the underlying paleotopography on loess deposition is critical to revealing the mechanism of loess‐landform formation. However, these controlling effects exhibit spatial variation as well as uncertainty, depending on a study's data sources, methodologies and particular research scope. In this study, the geological history of a study area in the Loess Plateau of China that is subject to severe soil erosion is investigated using detailed geological information and digital elevation models (DEMs), and an underlying paleotopographic model of the area is constructed. Based on the models of modern terrain and paleotopography, we introduce a watershed hierarchy method to investigate the spatial variation of the loess‐landform inheritance relationship and reveal the loess deposition process over different scales of drainage. The landform inheritance relationships were characterized using a terrain‐relief change index (TRCI) and a bedrock terrain controllability index (BTCI). The results show that the TRCI appears to have an inverse relationship with increasing research scope, indicating that, compared with the paleotopography of the region, modern terrain has lower topographic relief over the entire area, while it has higher topographic relief in the smaller, local areas. The BTCI strengthens with increasing drainage area, which demonstrates a strong controlling effect over the entire study area, but a weak effect in the smaller, local areas because of the effect of paleotopography on modern terrain. The results provide for an understanding of the spatial variation of loess deposition in relation to paleotopography and contribute to the development of a process‐based loess‐landform evolution model. Copyright © 2015 John Wiley & Sons, Ltd.     

Normal moveout velocity for pure‐mode and converted waves in layered orthorhombic medium

We study the azimuthally dependent hyperbolic moveout approximation for small angles (or offsets) for quasi‐compressional, quasi‐shear, and converted waves in one‐dimensional multi‐layer orthorhombic media. The vertical orthorhombic axis is the same for all layers, but the azimuthal orientation of the horizontal orthorhombic axes at each layer may be different. By starting with the known equation for normal moveout velocity with respect to the surface‐offset azimuth and applying our derived relationship between the surface‐offset azimuth and phase‐velocity azimuth, we obtain the normal moveout velocity versus the phase‐velocity azimuth. As the surface offset/azimuth moveout dependence is required for analysing azimuthally dependent moveout parameters directly from time‐domain rich azimuth gathers, our phase angle/azimuth formulas are required for analysing azimuthally dependent residual moveout along the migrated local‐angle‐domain common image gathers. The angle and azimuth parameters of the local‐angle‐domain gathers represent the opening angle between the incidence and reflection slowness vectors and the azimuth of the phase velocity ψ_phs at the image points in the specular direction. Our derivation of the effective velocity parameters for a multi‐layer structure is based on the fact that, for a one‐dimensional model assumption, the horizontal slowness and the azimuth of the phase velocity ψ_phs remain constant along the entire ray (wave) path. We introduce a special set of auxiliary parameters that allow us to establish equivalent effective model parameters in a simple summation manner. We then transform this set of parameters into three widely used effective parameters: fast and slow normal moveout velocities and azimuth of the slow one. For completeness, we show that these three effective normal moveout velocity parameters can be equivalently obtained in both surface‐offset azimuth and phase‐velocity azimuth domains.     

Multi‐scale relief model (MSRM): a new algorithm for the visualization of subtle topographic change of variable size in digital elevation models

Morphological analysis of landforms has traditionally relied on the interpretation of imagery. Although imagery provides a natural view of an area of interest (AOI) images are largely hindered by the environmental conditions at the time of image acquisition, the quality of the image and, mainly, the lack of topographical information, which is an essential factor for a correct understanding of the AOI's geomorphology. More recently digital surface models (DSMs) have been incorporated into the analytical toolbox of geomorphologists. These are usually high‐resolution models derived from digital photogrammetric processes or LiDAR data. However, these are restricted to relatively small areas and are expensive or complex to acquire, which limits widespread implementation. In this paper, we present the multi‐scale relief model (MSRM), which is a new algorithm for the visual interpretation of landforms using DSMs. The significance of this new method lies in its capacity to extract landform morphology from both high‐ and low‐resolution DSMs independently of the shape or scale of the landform under study. This method thus provides important advantages compared to previous approaches as it: (1) allows the use of worldwide medium resolution models, such as SRTM, ASTER GDEM, ALOS, and TanDEM‐X; (2) offers an alternative to traditional photograph interpretation that does not rely on the quality of the imagery employed nor on the environmental conditions and time of its acquisition; and (3) can be easily implemented for large areas using traditional GIS/RS software. The algorithm is tested in the Sutlej‐Yamuna interfluve, which is a very large low‐relief alluvial plain in northwest India where 10 000 km of palaeoriver channels have been mapped using MSRM. The code, written in Google Earth Engine's implementation of JavaScript, is provided as Supporting Information for its use in any other AOI without particular technical knowledge or access to topographical data. © 2017 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Application of Structure‐from‐Motion photogrammetry to river restoration

Structure‐from‐Motion (SfM) photogrammetry is now used widely to study a range of earth surface processes and landforms, and is fast becoming a core tool in fluvial geomorphology. SfM photogrammetry allows extraction of topographic information and orthophotos from aerial imagery. However, one field where it is not yet widely used is that of river restoration. The characterisation of physical habitat conditions pre‐ and post‐restoration is critical for assessing project success, and SfM can be used easily and effectively for this purpose. In this paper we outline a workflow model for the application of SfM photogrammetry to collect topographic data, develop surface models and assess geomorphic change resulting from river restoration actions. We illustrate the application of the model to a river restoration project in the NW of England, to show how SfM techniques have been used to assess whether the project is achieving its geomorphic objectives. We outline the details of each stage of the workflow, which extend from preliminary decision‐making related to the establishment of a ground control network, through fish‐eye lens camera testing and calibration, to final image analysis for the creation of facies maps, the extraction of point clouds, and the development of digital elevation models (DEMs) and channel roughness maps. The workflow enabled us to confidently identify geomorphic changes occurring in the river channel over time, as well as assess spatial variation in erosion and aggradation. Critical to the assessment of change was the high number of ground control points and the application of a minimum level of detection threshold used to assess uncertainties in the topographic models. We suggest that these two things are especially important for river restoration applications. Copyright © 2016 John Wiley & Sons, Ltd.     

面向对象技术在地球物理电磁成像可视化中的应用

面向对象技术作为一门新兴、热门的程序设计技术,当前被广泛应用到大规模数值分析、各种复杂可视化研究等领域.面向对象技术包含的类型继承、函数重载机制、虚函数多态机制成为设计完整、良好的电磁成像系统非常有利的因素.电磁成像系统需要大量的功能单元,数据的可视化需要功能完整的成像功能,这些需要都可以由面向对象技术中的对象得以表示、完成.另一方面,电磁成像系统的可视化优良性决定了地球物理电磁勘探结果的正确性,因此,十分有必要设计良好的电磁可视化系统.基于上述目的,本文基于面向对象技术,采用C++语言讨论了设计地球物理电磁成像系统可视化的方法结构,阐述了成像系统可视化设计中涉及的重要内容,并简述本文设计的可视化系统的应用.     

Digital landscapes of deglaciation: identifying Late Quaternary glacial lake outburst floods using LiDAR

High resolution DEMs obtained from LiDAR topographic data have led to improved landform inventories (e.g. landslides and fault scarps) and understanding of geomorphic event frequency. Here we use airborne LiDAR mapping to investigate meltwater pathways associated with the Tweed Valley palaeo ice‐stream (UK). In particular we focus on a gorge downstream of Palaeolake Milfield, previously mapped as a sub‐glacial meltwater channel, where the identification of abandoned headcut channels, run‐up bars, rock‐cut terrace surfaces and eddy flow features attest to formation by a sub‐aerial glacial lake outburst flood (GLOF) caused by breaching of a sediment dam, likely an esker ridge. Mapping of these landforms combined with analysis of the gorge rim elevations and cross‐section variability revealed a two phase event with another breach site downstream following flow blockage by higher elevation drumlin topography. We estimate the magnitude of peak flow to be 1–3 × 10³ m³/s, duration of the event to range from 16–155 days, and a specific sediment yield of 10⁷–10⁹ m³/km²/yr. We identified other outburst pathways in the lower Tweed basin that help delineate an ice margin position of the retreating Tweed Valley ice stream. The results suggest that low magnitude outburst floods are under‐represented in Quaternary geomorphological maps. We therefore recommend regional LiDAR mapping of meltwater pathways to identify other GLOFs in order to better quantify the pattern of freshwater and sediment fluxes from melting ice sheets to oceans. Despite the relatively low magnitude of the Till outburst event, it had a significant impact on the landscape development of the lower Tweed Valley through the creation of a new tributary pathway and triggering of rapid knickpoint retreat encouraging new regional models of post‐glacial fluvial landscape response. Copyright © 2015 John Wiley & Sons, Ltd.     

Evaluation of remotely‐sensed DEMs and modification based on plausibility rules and initial sediment budgets of an artificially‐created catchment

To quantify landscape change resulting from processes of erosion and deposition and to establish spatially distributed sediment budgets, ‘models of change’ can be established from a time series of digital elevation models (DEMs). However, resolution effects and measurement errors in DEMs may propagate to these models. This study aimed to evaluate and to modify remotely‐sensed DEMs for an improved quantification of initial sediment mass changes in an artificially‐created catchment. DEMs were constructed from photogrammetry‐based, airborne (ALS) and ground‐based laser scanning (TLS) data. Regions of differing morphological characteristics and vegetation cover were delineated. Three‐dimensional (3D) models of volume change were established and mass change was derived from these models. DEMs were modified region‐by‐region for rill, interrill and alluvial areas, based on logical and hydro‐geomorphological principles. Additional DEMs were constructed by combining multi‐source, modified data. Models were evaluated by comparison with d‐GPS reference data and by considering sediment budget plausibility. Comprehensive evaluation showed that DEM usability depends on a relation between the technique used to obtain elevation data, surface morphology and vegetation cover characteristics. Photogrammetry‐based DEMs were suited to quantification of change in interrill areas but strongly underestimated surface lowering in erosion rills. TLS DEMs were best suited to rill areas, while ALS DEMs performed best in vegetation‐covered alluvial areas. Agreement with reference data and budget plausibility were improved by modifications to photogrammetry‐ and TLS‐based DEMs. Results suggest that artefacts in DEMs can be reduced and hydro‐geomorphic surface structures can be better represented by applying region‐specific modifications. Photogrammetry‐based DEMs can be improved by combining higher and lower resolution data in defined structural units and applying modifications based on principles given by characteristic hydro‐geomorphic evolution. Results of the critical comparative evaluation of remotely‐sensed elevation data can help to better interpret DEM‐based quantifications of earth‐surface processes. Copyright © 2012 John Wiley & Sons, Ltd.     

Converted‐wave beam migration with sparse sources or receivers

Gaussian beam depth migration overcomes the single‐wavefront limitation of most implementations of Kirchhoff migration and provides a cost‐effective alternative to full‐wavefield imaging methods such as reverse‐time migration. Common‐offset beam migration was originally derived to exploit symmetries available in marine towed‐streamer acquisition. However, sparse acquisition geometries, such as cross‐spread and ocean bottom, do not easily accommodate requirements for common‐offset, common‐azimuth (or common‐offset‐vector) migration. Seismic data interpolation or regularization can be used to mitigate this problem by forming well‐populated common‐offset‐vector volumes. This procedure is computationally intensive and can, in the case of converted‐wave imaging with sparse receivers, compromise the final image resolution. As an alternative, we introduce a common‐shot (or common‐receiver) beam migration implementation, which allows migration of datasets rich in azimuth, without any regularization pre‐processing required. Using analytic, synthetic, and field data examples, we demonstrate that converted‐wave imaging of ocean‐bottom‐node data benefits from this formulation, particularly in the shallow subsurface where regularization for common‐offset‐vector migration is both necessary and difficult.