Mapping landforms of a hilly landscape using machine learning and high-resolution LiDAR topographic data

Landform maps are important tools in assessment of soil- and eco-hydrogeomorphic processes and hazards, hydrological modeling, and natural resources and land management. Traditional techniques of mapping landforms based on field surveys or from aerial photographs can be time and labor intensive, highlighting the importance of remote sensing products based automatic or semi-automatic approaches. In addition, the time-intensive manual labeling can also be subjective rather than an objective identification of the landform. Here we implemented such an objective approach applying a random forest machine learning algorithm to a set of observed landform data and 1m horizontal resolution bare-earth digital elevation model (DEM) developed from airborne light detection and ranging (LiDAR) data to rapidly map various landforms of a hilly landscape. The landform classification includes upland plateaus, ridges, convex slopes, planar slopes, concave slopes, stream channels, and valley bottoms, across a 400 km² hilly landscape of the Ozark Mountains in northeastern Oklahoma. We used 4200 landform observations (600 per landform) and eight topographic indices derived from 2 m, 5 m and 10 m resolution LiDAR DEM in random forest algorithm to develop 2 m, 5 m and 10 m resolution landform models. We test the effectiveness of DEM resolution in mapping landforms via comparison of observed landforms with modeled landforms. Results showed that the approach mapped ∼84% of observed landforms when covariates were at 2 m resolution to ∼89% when they were at 10 m resolution. However, predicted maps showed that the 2 m resolution covariates performed better at capturing accurate landform boundaries and details of small-sized landforms such as stream channels and ridges. The approach presented here significantly reduces the time required for mapping landforms compared to traditional techniques using aerial imagery and field observations and allows incorporation of a wide variety of covariates. The landform map developed using this approach has several potential applications. It could be utilized in hydrological modeling, natural resource management, and characterizing soil-geomorphic processes and hazards in a hilly landscape.