Linear interpolation of irregularly spaced LIDAR elevation data sets is needed to develop realistic spatial models. We evaluated inverse distance weighting (IDW) and ordinary kriging (OK) interpolation techniques and the effects of LIDAR data density on the statistical validity of the linear interpolators. A series of 10 forested 1000‐ha LIDAR tiles on the Lower Coastal Plain of eastern North Carolina was used. An exploratory analysis of the spatial correlation structure of the LIDAR data set was performed. Weighted non‐linear least squares (WNLS) analysis was used to parameterize best‐fit theoretical semivariograms on the empirical data. Tile data were sequentially reduced through random selection of a predetermined percentage of the original LIDAR data set, resulting in data sets with 50%, 25%, 10%, 5% and 1% of their original densities. Cross‐validation and independent validation procedures were used to evaluate root mean square error (RMSE) and kriging standard error (SE) differences between interpolators and across density sequences. Review of errors indicated that LIDAR data sets could withstand substantial data reductions yet maintain adequate accuracy (30 cm RMSE; 50 cm SE) for elevation predictions. The results also indicated that simple interpolation approaches such as IDW could be sufficient for interpolating irregularly spaced LIDAR data sets.

Detailed soil property maps are increasingly important for land management decisions and environmental modeling. The US Soil Survey is investing in production of the Soil Landscapes of the United States (SOLUS), a new set of national predictive soil property maps. This paper documents initial 100‐m resolution maps of 20 soil properties that include various textural fractions, physical parameters, chemical parameters, carbon, and depth to restrictions. Many of these properties have not been previously mapped at this resolution. A hybrid training strategy helped increase training data by roughly 10‐fold over previous similar studies by combining commonly used laboratory data with underutilized field descriptions tied to soil survey map unit component property estimates (to help represent within polygon variability) as well as randomly selected soil survey map unit weighted average property estimates. Relative prediction intervals were used to help select which training data sources improved model performance. Conventional and spatial cross‐validation strategies yielded generally strong coefficients of determination between 0.5 and 0.7, but with substantial variability and outliers among the various properties, types of training data, and depths. Internal review of the maps highlighted both strengths and weaknesses of the maps, but most of the critical comments were in areas with high model uncertainty that can be used to guide future improvements. Generally, previously glaciated areas and complex large alluvial basins were harder to model. The new SOLUS 100‐m maps will be updated in the future to address identified issues and feedback as users interact with the data.

The probabilistic mapping of landslide occurrence at a high spatial resolution and over a large geographic extent is explored using random forests (RF) machine learning; light detection and ranging (LiDAR)-derived terrain variables; additional variables relating to lithology, soils, distance to roads and streams and cost distance to roads and streams; and training data interpreted from high spatial resolution LiDAR-derivatives. Using a large training set and all predictor variables, an area under the receiver operating characteristic (ROC) curve (AUC) of 0.946 is obtained. Our findings highlight the value of a large training dataset, the incorporation of a variety of terrain variables and the use of variable window sizes to characterize the landscape at different spatial scales. We also document important variables for mapping slope failures. Our results suggest that feature selection is not required to improve the RF modeling results and that incorporating multiple models using different pseudo absence samples is not necessary. From our findings and based on a review of prior studies, we make recommendations for high spatial resolution, large-area slope failure probabilistic mapping.