Predicting soil bulk density for incomplete databases

Authors

Cleiton H. Sequeira, University of Nebraska-LincolnFollow
Skye A. Wills, USDA-NRCSFollow
Cathy A. Seybold, USDA-NRCSFollow
Larry T. West, USDA-NRCSFollow

Date of this Version

2014

Citation

Published in Geoderma 213 (2014) 64–73. DOI: 10.1016/j.geoderma.2013.07.013

Comments

This article is a U.S. government work and is not subject to copyright in the United States.

Abstract

Soil bulk density (ρ_b) is important because of its direct effect on soil properties (e.g., porosity, soilmoisture availability) and crop yield. Additionally, ρ_b measurements are needed to express soil organic carbon (SOC) and other nutrient stocks on an area basis (kg ha−1). However, ρbmeasurements are commonlymissing fromdatabases for reasons that include omission due to sampling constraints and laboratory mishandling. The objective of this study was to investigate the performance of novel pedotransfer functions (PTFs) in predicting ρ_b as a function of textural class and basic pedon description information extracted from the horizon of interest (the horizon for which ρb is being predicted), and ρ_b, textural class, and basic pedon description information extracted from horizons above or below and directly adjacent or not adjacent to the horizon of interest. A total of 2,680 pedons (20,045 horizons) were gathered from the USDA-NRCS National Soil Survey Center characterization database. Twelve ρ_b PTFs were developed by combining PTF types, database configurations, and horizon limiting depths. Different PTF types were created considering the direction of prediction in the soil profile: upward and downward prediction models. Multiple database configurations were used to mimic different scenarios of horizons missing ρ_b values: random missing (e.g., ρ_b sample lost in transit) and patterned or systematic missing (e.g., no ρ_b samples collected for horizons N 30 cm depth). For each database configuration scenario, upward and downward models were developed separately. Three limiting depths (20, 30, and 50 cm) were tested to identify any threshold depth between upward and downward models. For both PTF types, validation results indicated thatmodels derived from the database configuration mimicking randomhorizonsmissing ρ_b performed better than those derived from the configuration mimicking clear patterns of missing ρ_b measurements. All 12 PTFs performed well (RMSPE: 0.10–0.15 g cm−3). The threshold depth of 50 cm most successfully split the database between upward and downward models. For all PTFs, the ρ_b of other horizons in the soil profile was the most important variable in predicting ρ_b. The proposed PTFs provide reasonably accurate ρb predictions, and have the potential to help researchers and other users to fill gaps in their database without complicated data acquisition.