skip to main content


Title: Evaluating Landscape Complexity and the Contribution of Non‐Locality to Geomorphometry
Abstract

A long‐standing question in geomorphology concerns the extent that statistical models of terrain elevations have adequate characteristics with respect to the known scaling properties of landscapes. In previous work, it has been challenging to ascribe statistical significance to metrics adopted to measure landscape properties. Here, we use a recently developed surrogate data algorithm to generate synthetic surfaces with identical elevation values to the source data set, while also preserving the value of the Hölder exponent at any point (the underpinning characteristic of a multifractal surface). Our primary source data are from a laboratory experiment on landscape evolution. This allows us to examine how the statistical properties of the surfaces evolve through time and the extent to which they depart from the simple (multi)fractal formalisms. We show that there is a strong departure that is driven by the diffusive processes in operation. The number of sub‐basins of a given channel order (for orders sufficiently small relative to the basin order) exhibits a clear increase in complexity after a steady‐state for sediment flux is established. We also study elevation data from Florida and Washington States, where the relative departure from simple multifractality is even more strongly expressed but is similar for two very different locations. Our results show that at the very least, the minimum complexity for a stochastic model for terrain statistics with appropriate geomorphic scalings needs to incorporate a conditioning between the pointwise Hölder exponents and elevation.

 
more » « less
Award ID(s):
1811909 1839441
NSF-PAR ID:
10387030
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Earth Surface
Volume:
126
Issue:
4
ISSN:
2169-9003
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The effective characterization of topographic surfaces is a central tenet of geomorphology. Differences in land surface properties reveal variations in structural controls and the nature and efficacy of Earth‐shaping processes. In this paper, we employ the Hölder exponents,α, characterizing the local scaling behavior of topography and commonly used in the study of the (multi)fractal properties of landscapes and show that the joint probability distribution of the area of the terrain with a given elevation andαcontains a wealth of information on topographic structure. The conditional distributions of the hypsometric integrals as a function ofα, that is,Ihyp|α, are shown to capture this structure. A multivariate analysis reveals three metrics that summarize these conditional distributions: Strahler's original hypsometric integral, the standard deviation of theIhyp|α, and the nature of any trend of theIhyp|αagainstα. An analysis of five digital elevation models (DEMs) from different regions of the United States shows that only one is truly described by the hypsometric integral (Mettman Ridge from central Oregon). In the other cases, the new metrics clearly discriminate between instances where topographic roughness is more clearly a function of elevation, as captured by the conditional variables. In a final example, we artificially sharpen the ridges and valleys of one DEM to show that while the hypsometric integral and standard deviation ofIhyp|αare invariant to the change, the trend ofIhyp|αagainstαcaptures the changes in topography.

     
    more » « less
  2. Abstract

    Previous studies discovered a spatially heterogeneous expansion of Siberian larch into the tundra of the Polar Urals (Russia). This study reveals that the spatial pattern of encroachment of tree stands is related to environmental factors including topography and snow cover. Structural and allometric characteristics of trees, along with terrain elevation and snow depth were collected along a transect 860 m long and 80 m wide. Terrain curvature indices, as representative properties, were derived across a range of scales in order to characterize microtopography. A density-based clustering method was used here to analyze the spatial and temporal patterns of tree stems distribution. Results of the topographic analysis suggest that trees tend to cluster in areas with convex surfaces. The clustering analysis also indicates that the patterns of tree locations are linked to snow distribution. Records from the earliest campaign in 1960 show that trees lived mainly at the middle and bottom of the transect across the areas of high snow depth. As trees expanded uphill following a warming climate trend in recent decades, the high snow depth areas also shifted upward creating favorable conditions for recent tree growth at locations that were previously covered with heavy snow. The identified landscape signatures of increasing tall vegetation, and the effects of microtopography and snow may facilitate the understanding of treeline dynamics at larger scales.

     
    more » « less
  3. Remote sensing-based Earth Observation plays an important role in assessing environmental changes throughout our planet. As an image-heavy domain, the evaluation of the data strongly focuses on statistical and pixel-based spatial analysis methods. However, considering the complexity of our Earth system, there are some environmental structures and dependencies that are not possible to accurately describe with these traditional image analysis approaches. One example for such a limitation is the representation of (spatial) networks and their characteristics. In this study, we thus propose a computer vision approach that enables the representation of semantic information gained from images as graphs. As an example, we investigate digital terrain models of Arctic permafrost landscapes with its very characteristic polygonal patterned ground. These regular patterns, which are clearly visible in high-resolution image and elevation data, are formed by subsurface ice bodies that are very vulnerable to rising temperatures in a warming Arctic. Observing these networks’ topologies and metrics in space and time with graph analysis thus allows insights into the landscape’s complex geomorphology, hydrology, and ecology and therefore helps to quantify how they interact with climate change. We show that results extracted with this analytical and highly automated approach are in line with those gathered from other manual studies or from manual validation. Thus, with this approach, we introduce a method that, for the first time, enables upscaling of such terrain and network analysis to potentially pan-Arctic scales where collecting in-situ field data is strongly limited. 
    more » « less
  4. Assessment of lakes for their future potential to drain relied on the 2002/03 airborne Interferometric Synthetic Aperture Radar (IFSAR) Digital Surface Model (DSM) data for the western Arctic Coastal Plain in northern Alaska. Lakes were extracted from the IfSAR DSM using a slope derivative and manual correction (Jones et al., 2017). The vertical uncertainty for correctly detecting lake-based drainage gradients with the IfSAR DSM was defined by comparing surface elevation differences of several overlapping DSM tile edges. This comparison showed standard deviations of elevation between overlapping IfSAR tiles ranging from 0.0 to 0.6 meters (m). Thus, we chose a minimum height difference of 0.6 m to represent a detectable elevation gradient adjacent to a lake as being most likely to contribute to a rapid drainage event. This value is also in agreement with field verified estimates of the relative vertical accuracy (~0.5 m) of the DSM dataset around Utqiaġvik (formerly Barrow) (Manley et al., 2005) and the stated vertical RMSE (~1.0 m) of the DSM data (Intermap, 2010). Development of the potential lake drainage dataset involved several processing steps. First, lakes were classified as potential future drainage candidates if the difference between the elevation of the lake surface and the lowest elevation within a 100 m buffer of the lake shoreline exceeded our chosen threshold of 0.6 m. Next, we selected lakes with a minimum size of 10 ha to match the historic lake drainage dataset. We further filtered the dataset by selecting lakes estimated to have low hydrological connectivity based on relations between lake contributing area as determined for specific surficial geology types and presented in Jones et al. (2017). This was added to the future projection workflow to isolate the lake population that likely responds to changes in surface area driven largely by geomorphic change as opposed to differences in surface hydrology. Lakes within a basin with low to no hydrologic connectivity that had an elevation change gradient between the lake surface and surrounding landscape are considered likely locations to assess for future drainage potential. Further, the greater the elevation difference, the greater the drainage potential. This dataset provided a first-order estimate of lakes classified as being prone to future drainage. We further refined our assessment of potential drainage lakes by identifying the location of the point with the lowest elevation within the 100 m buffer of the lake shoreline and manually interpreted lakes to have a high drainage potential based on the location of the likely drainage point to known lake drainage pathways using circa 2002 orthophotography or more recent high resolution satellite imagery available for the Western Coastal Arctic Plain (WACP). Lakes classified as having a high drainage potential typically had the likely drainage location associated with one or more of the following: (1) an adjacent lake, (2) the cutbank of a river, (3) the ocean, (4) were located in an area with dense ice-wedge networks, (5) appeared to coincide with a potentially headward eroding stream, or (6) were associated with thermokarst lake shoreline processes in the moderate to high ground ice content terrain. We also added information on potential lake drainage pathways to the high potential drainage dataset by manually interpreting the landform associated with the likely drainage site to draw comparisons with the historic lake drainage dataset. 
    more » « less
  5. Abstract

    We collected ground‐penetrating radar (GPR) and frequency‐domain electromagnetic induction (FDEM) profiles in 2011 and 2012 to identify the extent of permafrost relative to surface biomass and solar insolation around Twelvemile Lake near Fort Yukon, Alaska. We compared a Landsat‐derived biomass estimate and modeled solar insolation from a digital elevation model to the geophysical measurements. We show correspondence between vegetation type and biomass relative to permafrost extent and seasonal freeze–thaw. Thicker permafrost (≥25 m) was covered by greater biomass, and seasonal thaw depths in these regions were minimal (1 m). Shallow (1–3 m depth) and thin (20–50 cm) newly forming permafrost or frozen layers from the previous winter occurred below northward oriented slopes with thin biomass cover. South‐facing slopes exhibited permafrost when there was enough biomass to shield incoming solar energy. We developed an artificial neural network to predict permafrost extent across the broader region by mapping GPR‐observed instances of permafrost to FDEM, biomass, and terrain observations with 90.2% accuracy. We identified a strong linear correlation (r = −0.77) between permafrost probability and seasonal thaw depth, indicating that our models may also be used to explore thaw patterns and variability in active layer thickness. This study highlights the combined influence of biomass and terrain on the presence of permafrost and the value of evaluating such parameters via remote sensing to predict permafrost spatial or temporal variability. Incorporating diverse geophysical datasets with in‐situ validation into machine learning models demonstrates a useful approach to upscale estimated permafrost extent across large Arctic expanses.

     
    more » « less