High-resolution Digital Terrain Models Research Articles

Groundwater flooding hazard assessment in a semi-urban aquifer through probability modelling of surrogate data

Many areas of the world are experiencing groundwater flooding (GF) events with serious impacts on the environment and human health. In Europe, GF has received specific attention with the Directive 2007/60/EC, which requires Member States to map GF hazard and risk and propose mitigation measures. In Italy, despite the commitment assigned by the Legislative Decree 49/2010 to River Basin District Authorities, no models have been developed for groundwater flooding hazard (GFH) evaluation also for the lack of standard methods based on commonly available data. In this paper, a new method has been developed through a specific procedure that accounts for i) probability of GF events estimation based on available time series of groundwater levels and ii) probability mapping through a newly developed MATLAB™ code, high-resolution topographic data and groundwater table model. The methodology has been tested and validated in the Metropolitan City of Naples (Italy), affected by GF since 2010, for which i) sixty-two-year long time series of groundwater levels, ii) high-resolution Digital Elevation and Digital Terrain Models and iii) basement/building geometry, are available. The results show that about 59% of the buildings present in the study area have potentially floodable underground structures. In the 2010–2021 period, only 2% of these buildings were flooded and most of them were built in the period 1980–2000, when the piezometric levels were deeper due to over-pumping of the groundwater. The obtained GFH map shows that only 2% of the built area is involved by 1- to 2- year events, while 6% is affected by 1- to 10- year events. A good spatial correspondence between the return period and locations of GF events occurred during 2010–2021 period was observed, thus indicating the significant performance of the model. The obtained GFH map output can support decisions of authorities and policy-makers to reduce the GF risk.

Reconstructing Digital Terrain Models from ArcticDEM and WorldView-2 Imagery in Livengood, Alaska

ArcticDEM provides the public with an unprecedented opportunity to access very high-spatial resolution digital elevation models (DEMs) covering the pan-Arctic surfaces. As it is generated from stereo-pairs of optical satellite imagery, ArcticDEM represents a mixture of a digital surface model (DSM) over a non-ground areas and digital terrain model (DTM) at bare grounds. Reconstructing DTM from ArcticDEM is thus needed in studies requiring bare ground elevation, such as modeling hydrological processes, tracking surface change dynamics, and estimating vegetation canopy height and associated forest attributes. Here we proposed an automated approach for estimating DTM from ArcticDEM in two steps: (1) identifying ground pixels from WorldView-2 imagery using a Gaussian mixture model (GMM) with local refinement by morphological operation, and (2) generating a continuous DTM surface using ArcticDEMs at ground locations and spatial interpolation methods (ordinary kriging (OK) and natural neighbor (NN)). We evaluated our method at three forested study sites characterized by different canopy cover and topographic conditions in Livengood, Alaska, where airborne lidar data is available for validation. Our results demonstrate that (1) the proposed ground identification method can effectively identify ground pixels with much lower root mean square errors (RMSEs) (<0.35 m) to the reference data than the comparative state-of-the-art approaches; (2) NN performs more robustly in DTM interpolation than OK; (3) the DTMs generated from NN interpolation with GMM-based ground masks decrease the RMSEs of ArcticDEM to 0.648 m, 1.677 m, and 0.521 m for Site-1, Site-2, and Site-3, respectively. This study provides a viable means of deriving high-resolution DTM from ArcticDEM that will be of great value to studies focusing on the Arctic ecosystems, forest change dynamics, and earth surface processes.

LiDAR-Derived Relief Typology of Loess Patches (East Poland)

The application of the automated analysis of remote sensing data processed into high-resolution digital terrain models (DTMs) using geographic information systems (GIS) tools provides a geomorphometric characterization of the diversity of the relief of loess patches over large areas. Herein, a quantitative classification of 79 loess patches with a total area of 3361 km2, distributed within the eastern part of the Polish Uplands belt, is carried out. A high-resolution 1 × 1 m DTM was generated from airborne laser scanning (ALS) data with densities ranging from 4 pts/m2 to 12 pts/m2, which was resampled to a resolution of 5 × 5 m for the study. This model was used to classify landform surfaces using the r.geomorphon (geomorphon algorithm) function in GRASS GIS software. By comparing the values in the neighborhood of each cell, a map of geomorphometric features (geomorphon) was obtained. The classification and typology of the relief of the studied loess patches was performed using GeoPAT2 (Geospatial Pattern Analysis Toolbox) software. Pattern signatures with a resolution of 100 × 100 m were extracted from the source data grid, and the similarity of geomorphological maps within the signatures was calculated and saved as a signature file and segment map using the spatial coincidence method. The distance matrix between each pair of segments was calculated, and the heterogeneity and isolation of the maps were generated. R system was used to classify the segments, which generated a dendrogram and a heat map based on the distance matrix. This made it possible to distinguish three main types and eight subtypes of relief. The morphometric approach used will contribute to a better understanding of the spatial variation in the relief of loess patches.

Assessing the Extent of Flood-Prone Areas in a South-American Megacity Using Different High Resolution DTMs

Current forecasts estimate that almost 68% of the global population will be living in urban centers by 2050. As a result, the increase in impermeable surface area can result in severe hydrological impacts, such as the increase in surface runoff and the frequency of floods and their magnitude. Thus, this work analyzes the performance of the hydrodynamic model HEC-RAS for assessing the extent of flood-prone areas, using two digital terrain models (DTM) with different spatial resolutions (5 and 0.50 m). Four different computing intervals (1, 15, 30, and 60 s) were adopted aiming to evaluate the simulations outputs performance. Additionally, reported data by the civil defense are used for calibration and validation. In general, the model showed to be a powerful tool in the identification of susceptible areas to urban flooding. The simulated results in this work provide crucial geographic information when identifying spots with the highest risk of flooding, which should receive priority attention during such events. The simulations with a spatial resolution of 5 m showed the flood maps with the largest coverage of the flooded points (278 points out of 286—97.20%), within the shortest computation times. We highlight that the more refined DTM derived from spatial images did not produce the best flood simulation compared to the DTM with a spatial resolution of 5 m derived from orthoimages.

Rediscovering the Lost Roman Landscape in the Southern Trieste Karst (North-Eastern Italy): Road Network, Land Divisions, Rural Buildings and New Hints on the Avesica Road Station

An interdisciplinary study of the ancient landscape of the Trieste Karst (north-eastern Italy) is presented in this paper. Airborne Laser Scanning (ALS) has been applied to obtain high-resolution topography of the 25 km2 investigated area in order to identify potential archaeological anomalies. The ALS-derived high-resolution Digital Terrain Models have been visualized and managed using QGIS and Relief Visualization Toolbox. Possible archaeological anomalies have been verified through field surveys and interpreted using a multidisciplinary approach mainly based on the collection of associated archaeological materials and geomorphological and stratigraphic evidence. From a methodological perspective, the elaboration and study of ALS-derived images, and in particular the local relief model visualization, combined with the collection of Roman shoe hobnails, have proven to be effective approaches for the certain identification and dating of Roman roads in karst environments. The obtained results have revealed an almost completely unknown Roman landscape: the investigated area was crossed by important public roads, whose layout has been accurately reconstructed for a total length of over 10 km, and occupied by large country estates, sometimes enclosed within boundary walls perfectly fitting the Roman land division grid. One of the identified buildings could correspond to a road station, perhaps the Avesica known from ancient itinerary documents—i.e., the itinerarium Antonini Augusti—due to its position and proximity to a major road junction.

Deciphering complexity of block-faulted topography – Insights from morphometric variables, drainage patterns and river profiles (Orlickie-Bystrzyckie Mountains Block, Central Europe)

The Orlickie-Bystrzyckie Mountains Block (OBMB) is a medium-size, transboundary morphostructure within the Sudetes range in Central Europe, subject to non-uniform uplift in the late Cenozoic. Geologically, it consists of Variscan basement and Permian to Cretaceous sedimentary cover. This paper attempts to unravel its morphotectonic history and decipher the spatial pattern of uplift and erosional response, using several complementary, quantitative approaches to characterize topography. The focus is on: (1) general relief properties captured by selected topographic variables; (2) fault-generated escarpments, which define the eastern and northern limits of the study area; (3) drainage patterns and its peculiarities; (4) diversity of longitudinal river profiles and valley morphology. A LiDAR-based high-resolution digital terrain model was the primary source of elevation data. An important component of the study was examination of the role of lithological diversity on relief properties. A statistical approach generally revealed a secondary role of the rock factor and provided a foundation to interpret the macro-topography of the OBMB in terms of tectonic controls. Based on geomorphometric characteristics and information about the distribution of sedimentary rocks, which can be used as marker horizons to identify vertical displacements, four major and one secondary centres of surface uplift have been recognized. They are all asymmetric, being bounded by distinct fault-generated escarpments on one side and sloping gently in the opposite direction. Erosional response to uplift has been variable and includes, apart from simple incision into rising crustal blocks, the origin of epigenetic-antecedent water gaps, drainage diversions, headward erosion, and captures. This study also reveals that morphometric indices routinely used to extract the tectonic signal from landform patterns do not necessarily provide consistent information.

Resolving Ambiguities in SHARAD Data Analysis Using High-Resolution Digital Terrain Models

The SHAllow RADar (SHARAD) onboard Mars Reconnaissance Orbiter (MRO) is a 20 MHz Synthetic Aperture Radar (SAR) that probes the first hundreds of meters of the Martian subsurface. In order to interpret the detection of subsurface interfaces with ground penetrating radars, simulations using Digital Terrain Models (DTM) are necessary. This methodology paper focuses on the analysis of the first tens of meters of the Martian subsurface with SHARAD, comparing the use of different high-resolution DTMs for radar simulation, namely, from the High-Resolution Stereo Camera (HRSC) onboard the Mars Express and from the Context Camera (CTX) onboard MRO. The region of Terra Cimmeria was chosen as a demonstration area. It is a highly cratered southern midlatitude region, where, as will be discussed, the higher resolution of the aforementioned terrain models is mandatory to describe the surface at an acceptable level of detail for shallow subsurface radar interpretation. With a DTM corrected by photoclinometry using CTX imagery, we show that a reflector that was visible on SHARAD data but not on the simulation made with an HRSC DTM is, in fact, a surface echo that was not reproduced by the HRSC surface model. We also show that, unlike laser altimetry DTMs, optical DTMs are prone to artifacts that can make radar analysis more complicated for some scenarios. Reciprocally, we show that the comparison between radar and its corresponding simulated data is a way of assessing a DTM’s quality, which is especially useful in missions where ground control points are lacking, unlike Martian observations.

Building a High-resolution Digital Terrain Model of Bennu from Laser Altimetry Data

The Origins, Spectral Interpretation, Resource Identification, and Security-Regolith Explorer (OSIRIS-REx) spacecraft orbited the near-Earth asteroid (101955) Bennu to characterize the asteroid prior to sampling. One important aspect of this characterization was the creation of a high-resolution (5–7 cm) global shape model using the OSIRIS-REx Laser Altimeter (OLA). We describe the data collected by OLA, along with the approach used to register overlapping topography using keypoints and keypoint descriptors in order to produce a globally self-consistent set of data. These globally registered sets of topographic scans were used to generate digital terrain models at both global and regional scales. We also describe efforts to correct for a change in behavior of the scanning mirror after the launch and highlight the improvements to the data after implementing an updated calibration of the mirror. The resulting model represents the highest-fidelity global OLA data set.

Power laws for accurate determination of landslide volume based on high-resolution LiDAR data

Landslide volume is critical in landslide risk management due to its close association with the landslide mobility and damage to elements at risk. Limited by accessibility, site conditions and availability of equipment, directly measuring the landslide volume on site is challenging. Instead, many empirical area-volume power-law models have been developed to estimate the landslide volume. Yet the validity of these power-law models for small-scale shallow landslides is in question considering the substantial uncertainty of field-estimated landslide volume data for developing these empirical power laws. In this study, two high-resolution LiDAR-derived digital terrain models taken in 2010 and 2020 are leveraged to interpret the volumes of 1326 shallow landslides that occurred in Hong Kong during the same period. The high precision of the digital terrain models ensures the accurate interpretation of landslides and the identification of landslide source and transport areas. On the basis of the interpreted landslide characteristics, new power-law models and a multivariate distribution model are proposed to accurately estimate the volumes of small-scale shallow landslides based on easy-to-obtain landslide source features such as source area, length, width or maximum depth. The multivariate distribution model rigorously expresses the conditional probability distribution of the landslide volume and quantifies the uncertainty of volume estimation. The models therefore enhance the capability of landslide risk assessment and management.

Modeling the thermal environment of Mercury’s north pole using MLA. Implications for locations of water ice

Radar observations made by the Arecibo Observatory have revealed radar-bright deposits in the north polar region of Mercury which are consistent with water ice. These deposits are generally associated with floors of deep craters in permanent shadow. In this paper, we revisit topographic data, illumination, and associated thermal conditions of the north polar region of Mercury to investigate these radar-bright deposits. We used Mercury Laser Altimeter (MLA) data from the MESSENGER mission to derive a dedicated, adjusted high-resolution Digital Terrain Model (DTM) of the north polar region at 250 m/pixel. Based on this DTM, illumination conditions and temperature were derived for the first time on the fully available MLA dataset. Generally, craters equatorward of 86 ∘N are found to be too hot to contain surface water ice due to effective self-heating. For those warm craters that do contain radar-bright deposits, temperatures within the upper 10–20 cm of regolith are cold enough for water ice to be stable. Using our improved DTM, we confirm previously MLA-modeled illumination conditions in Prokofiev crater and argue that radar-bright patches outside the PSR in the area may be explained by icy deposits in cold and shallow (within 5 cm) sub-surface layers. In the radar observations we find a strong, unphysical correlation between radar-bright pixels and the number of radar observations per pixel. Consequently, several PSRs which were classified as free of radar-bright material (mainly at longitudes 120–340 ∘E) are most likely subject to a bias in the averaging process in deriving the composite radar image and might contain water ice after all.

Enhancing the Identification and Mapping of Fluvial Terraces Combining Geomorphological Field Survey with Land-Surface Quantitative Analysis

A methodological approach to refining the identification and mapping of fluvial terraces has been applied, combining geomorphological field surveys with the computation and assessment of different morphometric parameters (local, statistical, and object-oriented), derived from a high-resolution digital terrain model (DTM) obtained from a LiDAR survey. The mid-sector floodplain of the Misa River basins was taken as a valid example of the main river valleys draining the northern Marche Apennines (Italy) and was considered an ideal site to test a combination of different geomorphological techniques for enhancing fluvial terraces’ detection and mapping. In this area, late Pleistocene–Holocene fluvial terraces are well exposed, and their geomorphological and geochronological characteristics have largely already been studied. However, a reliable distinction of the different Holocene terrace levels, including a detailed geomorphological mapping of different terrace features, is still lacking due to the very complex terrace geometry and the lack of good-quality deposit outcrops. Land-surface quantitative (LSQ) analysis has been coupled with the available outcomes of previous studies and ad-hoc geomorphological field surveys to enhance the identification and mapping of fluvial terraces. The results of this work provided information for the discernment of terrace remnants belonging to the full-glacial fill terrace generation (late Pleistocene) as well as reconstruction of the terrace top–surface, and can be used to distinguish the inner terrace limits coinciding with the margin of the floodplain. It has also been possible to identify and delimit the late Pleistocene terrace from a staircase of three younger strath terraces formed during the Holocene. The results of this study demonstrated that the investigation of fluvial landforms, at different scales, can strongly benefit from the integration of field surveys and quantitative geomorphic analysis based on high-resolution digital topographic datasets. In particular, the integration of LSQ analysis with ground-truth geomorphological data can be dramatically helpful for the identification and mapping of fluvial terraces.

A semi-automated, GIS-based framework for the mapping of supraglacial hydrology

AbstractSupraglacial drainage networks play an integral role in both glacier dynamics and run-off timing, and mapping them provides insight into their role in glacial systems. Here we present a reproducible approach for semi-automated mapping of supraglacial hydrologic features, which complements existing work in automated and manual mapping by providing clear definitions for identification of features. This framework uses a digital terrain model (DTM) to identify potential flow routes on the glacier surface, which are then classified using a set of standardized rules based on the DTM and an orthomosaic. We found that the normalized difference water index calculated from digital imagery was influenced by image brightness and introduce a new approach using average RGB values to correct for this. Using this framework we mapped supraglacial drainage networks at Nàłùdäy and Thores Glacier, Canada. The framework was easier to implement with high-resolution (0.5 m) imagery and DTMs, compared to data with lower resolution (10 m), due to the increased detail in topography and feature boundaries at high-resolution. Lower-resolution data captured larger streams (>2 pixels wide), however, indicating that the framework can still be used at this resolution. Mapping supraglacial hydrology using standardized methods opens possibilities for investigating many questions relating to changes in supraglacial hydrology over time.

The methodology of determining optimum access routes to remote areas for the purposes of crisis management

ABSTRACT The article presents a methodology to determine access routes to places with difficult access, situated outside the regular transport network. To determine the optimum route, the existing road network was used, along with terrain passability maps, as well as high-resolution digital terrain models together with the traction parameters of the vehicle to enable conducting a detailed analysis of micro-relief and excluding inaccessible areas. The research was extended to include a comprehensive analysis of terrain accessibility for areas located in the vicinity of the capital cities of three European countries on various levels of economic development: Germany, Poland, and Ukraine. The degree of development of the analysed area was determined with the use of the gross domestic product (GDP), human development index (HDI), and the road quality coefficient. The case study presented in the paper illustrates the practical application of the developed methodology that was implemented to a real crisis situation. Access routes were determined for a MAN TGM 15.290 vehicle that went to take part in emergency action in buildings situated approximately 500 m away from paved roads. The developed methodology may be widely used in crisis management, where the emergency actions often require rescue units to access remote places.

Subalpine woody vegetation in the Eastern Carpathians after release from agropastoral pressure

The subalpine vegetation in the Eastern Carpathians has been under agropastoral influence as a high-mountain open pasture for about five centuries. Today, the subalpine zone released by human intervention is growing as thickets. In this study, we use a numerical model of tree crowns (CHM, Canopy Height Model) based on laser scanning (LiDAR) and a high-resolution digital terrain model (DTM) to delineate the subalpine thicket distribution. Anselin ‘Local Moran's I’ statistic was used to find hot and cold spots in vegetation cover. We used a logistic generalized linear model (GLM) and Principal Component Analysis (PCA) to set for the historical, climatic and terrain conditions candidates as the predictors of the present-day distribution of vegetation hot spots. We use variance partitioning to assess the interaction of climate and terrain variables. The resulting model suggests key environmental controls that underlie the vegetation pattern. Namely, snow in terrain depressions protects woody vegetation against abrasion and winter drought and increased insolation reduces the site humidity in the summer on S-E exposure hampering re-vegetation. In addition, the increasing distance from the treeline declines the rate of secondary succession. In all, the spatial model predicts the 35% coverage by thickets as a theoretical maximum of available climatic-terrain niches. The results suggest that the growth of the subalpine thicket, in the face of growing global temperature, may be restricted due to the limited number of niches available.

Remote-sensing detection of points for the heavy vehicles crossing in the large lowland riverbeds

Abstract The erosion-resistant culminations of the alluvial subsoil occurring in the riverbeds of some rivers may play the role of “geological bridges” under favorable hydrological conditions, constituting places convenient for forcing by heavy vehicles. The article discusses the possibility of remote-sensing identification of such places based on a common analysis of aerial images, high-resolution satellite images and Digital Terrain Models based on Airborne Laser Scanning. The result of the work was the identification of relief forms in the valley bottoms of large lowland indicating the available locations for heavy vehicles forcing the bed. This allowed for the formulation of general assumptions of the SYSTEM OF DESIGNING ZONES CONVENIENT FOR CROSSROADS USING AIR PLATFORMS.

The Use of High-Resolution Remote Sensing Data in Preparation of Input Data for Large-Scale Landslide Hazard Assessments

The objective of the study is to show that landslide conditioning factors derived from different source data give significantly different relative influences on the weight factors derived with statistical models for landslide susceptibility modelling and risk analysis. The analysis of the input data for large-scale landslide hazard assessment was performed on a study area (20.2 km2) in Hrvatsko Zagorje (Croatia, Europe), an area highly susceptible to sliding with limited geoinformation data, including landslide data. The main advantage of remote sensing technique (i.e., LiDAR, Light Detection and Ranging) data and orthophoto images is that they enable 3D surface models with high precision and spatial resolution that can be used for deriving all input data needed for landslide hazard assessment. The visual interpretation of LiDAR DTM (Digital Terrain Model) morphometric derivatives resulted in a detailed and complete landslide inventory map, which consists of 912 identified and mapped landslides, ranging in size from 3.3 to 13,779 m2. This inventory was used for quantitative analysis of 16 input data layers from 11 different sources to analyse landslide presence in factor classes and thus comparing landslide conditioning factors from available small-scale data with high-resolution LiDAR data and orthophoto images, pointing out the negative influence of small-scale source data. Therefore, it can be concluded that small-scale landslide factor maps derived from publicly available sources should not be used for large-scale analyses because they will result in incorrect assumptions about conditioning factors compared with LiDAR DTM derivative factor maps. Furthermore, high-resolution LiDAR DTM and orthophoto images are optimal input data because they enable derivation of the most commonly used landslide conditioning factors for susceptibility modelling and detailed datasets about elements at risk (i.e., buildings and traffic infrastructure data layers).

Transverse aeolian ridges in the landing area of the Tianwen-1 Zhurong rover on Utopia Planitia, Mars

The enigmatic transverse aeolian ridges (TARs), with distinct morphology and albedo, are among the key geological features investigated by China's Tianwen-1 Zhurong rover on southern Utopia Planitia, Mars. Their morphologies and morphometrics are investigated through high-resolution imaging science experiment (HiRISE) orthoimage and Digital Terrain Model (DTM) products. A total of 5089 TARs are identified, with barchan TARs being predominant (97.6%). Morphometric analysis shows these TARs to be small and symmetrical aeolian landforms, with an average crest–ridge lengths of 33.9 ± 20.5 m, profile widths of 9.4 ± 3.8 m, profile heights of 0.4 ± 0.4 m, profile-height-width ratios of 0.04 ± 0.02, and profile symmetry ratios of −0.01 ± 0.13. In-situ observations from the Navigation and Terrain Camera (NaTeCam) show the crests of the TARs to be dark and sharp, while the flanks are interlaced by dark and bright materials. Close-up Multispectral Camera (MSCam) images reveal the TARs to be coated by granules of ∼1.5 mm in diameter. Given the morphometric characteristics and the presence of coating granules, the TARs in the landing area could be categorized as megaripples. Buffered crater counting (BCC) technique-derived absolute model age (AMA) reveals the formation time, or the last active period of the TARs, could be as recent as 1 Ma in the Late Amazonian. The morphometrics and direction of the horns of the barchan TARs suggest the winds for the formation of TARs blew mostly from the north. During the spring-summer transition period (Ls: 50°–93°), the Mars Climate Station (MCS) had recorded local bimodal winds in the landing area, with the speed of the northerly wind in the afternoon being a little stronger than the speed of the southerly wind in the morning. These observations are consistent with the wind fields described in the Mars Climate Database (MCD), which imply the northerly winds during the northern winter season to be responsible for the net sediment transport to the south. Two TARs observed in-situ with secondary NW-SE trending crest-ridges indicate that forked TARs might form given sufficient time (i.e., in the order of millions of years) under modern wind conditions, i.e., the TARs may be currently reworked, if only extremely weakly.

The role of microtopography and resident species in post-disturbance recovery of arid habitats in Hawai'i.

Habitat-suitability indices (HSI) have been employed in restoration to identify optimal sites for planting native species. Often, HSI are based on abiotic variables and do not include biotic interactions, even though similar abiotic conditions can favor both native and nonnative species. Biotic interactions such as competition may be especially important in invader-dominated habitats because invasive species often have fast growth rates and can exploit resources quickly. In this study, we test the utility of an HSI of microtopography derived from airborne LiDAR to predict post-disturbance recovery and native planting success in native shrub-dominated and nonnative, invasive grass-dominated dryland habitats in Hawai'i. The HSI uses high-resolution digital terrain models to classify sites' microtopography as high, medium, or low suitability, based on wind exposure and topographic position. We used a split-plot before-after-control-impact design to implement a disturbance experiment within native shrub (Dodonaea viscosa) and nonnative, invasive grass (Cenchrus clandestinus)-dominated ecosystems across three microtopography categories. In contrast to previous studies using the same HSI, we found that microtopography was a poor predictor of pre-disturbance conditions for soil nutrients, organic matter content, or foliar C:N, within both Dodonaea and Cenchrus vegetation types. In invader-dominated Cenchrus plots, microtopography helped predict cover, but not as expected (i.e., highest cover would be in high-suitability plots): D. viscosa had the greatest cover in low-suitability and C. clandestinus had the greatest cover in medium-suitability plots. Similarly, in native-dominated Dodonaea plots, microtopography was a poor predictor of D. viscosa, C. clandestinus, and total plant cover. Although we found some evidence that microtopography helped inform post-disturbance plant recovery of D. viscosa and total plant cover, vegetation type was a more important predictor. Important for considering the success of plantings, percent cover of D. viscosa decreased while percent cover of C. clandestinus increased within both vegetation types 20 months after disturbance. Our results are evidence that HSIs based on topographic features may prove most useful for choosing planting sites in harsh habitats or those already dominated by native species. In more productive habitats, competition from resident species may offset any benefits gained from "better" suitability sites.

Spatial Modelling of Vineyard Erosion in the Neszmély Wine Region, Hungary Using Proximal Sensing

Intense soil erosion in the northern part of the Gerecse Hills, Hungary, is causing significant damage to vineyards in the area. Three vineyards in the Neszmély Wine Region were investigated to quantify the amount of eroded soils. The method was based on monitoring vineyards for one-year between July 2019 and June 2020. Every season, a set of photographs of the vineyards were taken from an unmanned aerial vehicle. The images were processed in a photogrammetric workflow to produce high-resolution digital terrain models (DTMs) and orthophotos, which were used to estimate the soil loss using the Universal Soil Loss Equation (USLE) model. Particular attention was paid to the effect of seasonal variation in vegetation cover and rainfall, and the erosion control effect of the inter-row grassing already applied in the vineyards was also modelled. The results confirm and quantify the extent to which intense summer rainfall has a more significant effect on erosion compared to autumn or winter rainfall.

Environmental factors controlling the Last Glacial multi-phase development of the Moravian Sahara dune field, Lower Moravian Basin, Central Europe

Aeolian sediments and landforms represent a unique source of proxy data allowing the detailed research of palaeoenvironmental changes. For example, dune morphology and dune field patterns could reflect an influence of a wide range of conditions from palaeowind directions to topography. The Moravian Sahara dune field located in southern Moravia represents a unique archive reflecting the effects of these environmental variables, lying far away from the European sand belt which is the main area of interest for aeolian research in Europe. Here, we present the complex geospatial analysis of the dunes and their topography performed on Digital Terrain Models (DTM) from which their basic morphometrics (crest orientation, spacings between the dunes, and others) were obtained. Two groups of differently oriented transverse dunes have been distinguished, suggesting a significant twist in atmospheric circulation connected with the retreat of ice sheets following the Last Glacial Maximum. In addition, other local conditions such as the topography and sediment supply were found to exert a strong influence on the formation of dune field patterns. Thus, the results of the high-resolution DTMs presented here show their applicability to the highly detailed study of the influence of environmental conditions on the development of the aeolian dune fields.