HYSOGs250m, global gridded hydrologic soil groups for curve-number-based runoff modeling

Soils have a fundamental role in the global hydrologic cycle by governing rainfall infiltration and groundwater recharge, which ultimately affects the lateral transport of water and subsequent runoff potential. Knowledge of soil hydraulic properties is therefore of interest to ecologists, hydrologists, and soil scientists, and is critical for parameterization of a variety of empirical and physically-based hydrologic models, dynamic-vegetation models, and land-surface models^1–3.

The U.S. Department of Agriculture (USDA) curve-number (CN) method provides a simplified approach to the estimation of key hydrologic processes while being grounded in a physical understanding of saturated flow and runoff processes^4–6. The CN method avoids the problems inherent to parameterizing and running more complex models due to its simplicity and relatively low data input requirements, and has been implemented in a variety of hydrologic, erosion, and water-quality models^7–9. CN selection is derived from the hydrologic response of various combinations of soil types and land cover classes^2,10. Particularly relevant to the subject of this analysis, and the data product we make available, is the classification and development of soil parameters for CN-based runoff modeling. The lack of globally consistent data derived from contemporary soil information served as the overarching motivation for this analysis.

CN-based runoff estimates require information regarding the minimum infiltration rate of rainfall into the soil and the transmission rate of groundwater through the soil profile after prolonged wetting. Runoff occurs when the rainfall rate exceeds the infiltration capacity of soils. The rate at which these processes occur is primarily affected by the physical nature of soils (e.g., texture, compaction), in addition to land cover, antecedent moisture, and rainfall intensity. For example, coarse-textured sandy soils have larger pore spacing, allowing water to infiltrate quickly relative to fine-textured clay soils.

Soils are thus classified into four hydrologic soil groups (HSGs) to infer runoff potential (Table 1)¹¹. HSG-A has the lowest runoff potential (typically contains more than 90% sand and less than 10% clay), HSG-B has moderately low runoff potential (typically contains between 10 to 20% clay and 50 to 90% sand), HSG-C has moderately high runoff potential (typically contains between 20 to 40% clay and less than 50% sand), and HSG-D has high runoff potential (typically contains more than 40% clay and less than 50% sand). Classification is determined by the least transmissive soil layer—often measured as saturated hydraulic conductivity (K_s)—depth to water table or depth to an impermeable layer (e.g., duripan, bedrock). If K_s is unknown or not available, infiltration and transmission rates can be inferred from soil texture, with the underlying assumption that soils with similar content of sand, silt, and clay have analogous hydraulic properties^12–14. Wet soils have high runoff potential (regardless of texture) due to the presence of a groundwater table within 60 cm of the surface. These soils are assigned dual HSGs, as a less restrictive group can be assigned (according to texture or K_S) if they can be adequately drained.

We derived HSGs from texture classes in accordance with USDA¹¹ specifications (Table 1). The resulting data product—HYSOGs250m—represents typical soil runoff potential suitable for regional, continental, and global scale analyses and is available in a gridded format at a spatial resolution of 250 m (Fig. 1).

Hydrologic soil groups A, B, C, and D correspond to low, moderately low, moderately high, and high runoff potential, respectively. Wet soils are assigned a dual HSG (e.g., HSG A/D) and have high runoff potential due to the presence of a water table within 60 cm of the surface. A less restrictive group can be assigned if these soils are drained (e.g., HSG-A).

Full size image

Our analysis indicates that soils with moderately high runoff potential dominate the global distribution (57.4%), followed by soils with moderately low (HSG-B 12.2%), high (HSG-D 10.1%), and low runoff potential (HSG-A 3.0%) (Table 2). Dual HSGs A/D, B/D, C/D, and D/D accounted for 0, 1.4, 13.5, and 2.4% of the global distribution, respectively. Some global trends were observed for soils with high and low runoff potential. Low runoff potential soils are found predominantly in parts of the Sahara and Arabian Deserts, which are characterized by very deep and well-drained sandy soils. High runoff potential soils occur predominantly within tropical and sub-tropical zones (with notable additions occurring in the Alaska-Yukon Arctic and Canadian Taiga and Boreal Shield) and are characterized by soils with high clay content or shallow soils (<50 cm to bedrock). No clear pattern could be discerned for soils with moderately low runoff potential at the global scale, as these HSGs occur in arid and humid environments and at both high and low elevations.

The process for producing HYSOGs250m consisted of five primary steps (Fig. 2). We classified HSGs from USDA-based soil texture classes (Fig. 3), depth to bedrock (Fig. 4), and groundwater table depth (Fig. 5) as specified by the USDA-Natural Resources Conservation Service (USDA-NRCS) National Engineering Handbook (NEH)¹¹. Soil texture classes and depth to bedrock were obtained from the SoilGrids predictions (soilgrids.org) Food and Agriculture Organization (FAO) soilGrids250m system¹⁵. These data and associated meta-data are available for download as GeoTiffs at ftp://ftp.soilgrids.org/data/recent. Groundwater table depth¹⁶ and associated meta-data are available for download as NetCDF at https://glowasis.deltares.nl/thredds/catalog/opendap/opendap/Equilibrium_Water_Table/catalog.html. All computations were performed within the R open source environment for statistical computing¹⁷ and functions from the raster package¹⁸.

textStack represents USDA-based soilGrids250m texture classes¹⁵ for six depth intervals (0, 5, 15, 30, 60, and 100 cm); hsgStack represents hydrologic soil group (HSGs) classified from each texture class, maxHSG represents HSGs defined by the most restrictive layer (0 to 1 m), R represents bedrock depth¹⁵, maxHSG_R represents HSGs re-classified to the bedrock depth criteria, and waterTable represents the HSGs reclassified to account for both the depth to bedrock and the water table criteria.

Full size image

Adapted from SoilGrids predictions (soilgrids.org)¹⁵. Cl is clay, SiCl is silty clay, SaCl is sandy clay, ClLo is clay loam, SiClLo is silty clay loam, SaClLo is sandy clay loam, Lo is loam, SiLo is silty loam, SaLo is sandy loam, Si is silt, LoSa is loamy sand, Sa is sand. Note that mapped texture classes represent the soil surface (0 cm).

Full size image

Adapted from SoilGrids predictions (soilgrids.org)²⁰. Note that individual grid cells (bedrock occurrence) may not be visible at the global scale.

Full size image

Adapted from Fan et al.¹⁶ Dual hydrologic soil groups were assigned to grid cells (pedons) based upon the presence of a water table (<60 cm of the surface).

Full size image

Soil texture to 1 m depth was represented with SoilGrids predictions (soilgrids.org) soilGrids250m texture classes at six depths: 0, 5, 15, 30, 60, and 100 cm. The soilGrids were stacked into a multi-band raster (textStack) using the raster::stack function (Fig. 2a). For the purpose of this analysis, we refer to individual grid cells (~250 m×250 m) in the raster stack (1 m depth) as soil pedons. Each grid cell in the raster stack (or pedon) was re-classified into one of four HSGs (hsgStack) using the classification scheme reported in Table 1 (Fig. 2b). This allowed us to infer the water transmissivity of each layer in the profile from the stacked texture classes. Note that integers 1, 2, 3, and 4 were used to represent HSGs A, B, C, and D, respectively. The raster::max function (Fig. 2c) was then used to determine the largest value of each grid cell in the raster stack, allowing us to infer the most restrictive layer in the pedon. This value (maxHSG) was used to assign HSGs for each pixel in the stack, thus representing soil runoff potential for each pedon. Shallow soils (bedrock within 50 cm of the surface, Fig. 4) were re-classified to HSG-D (maxHSG_R, Fig. 2d). Dual HSGs were assigned to pedons with shallow water tables (<60 cm from the surface) using the depth to groundwater table dataset¹⁶ (Fig. 2e). Integers 11, 12, 13, and 14 were used to denote dual HSGs A/D, B/D, C/D, and D/D, respectively.

Code availability

The R code used to develop HYSOGs250m, described in Fig. 2, is available for download from the Oak Ridge National Laboratory (ORNL) Distributed Active Archive Center (DAAC) (Data Citation 1).

HYSOGs250m (Data Citation 1) is available for download as an un-projected GeoTiff at 7.5 arc-second (approximately 250 m resolution). The value column variables 1, 2, 3, 4, 11, 12, 13, and 14 correspond to HSG A, B, C, D, A/D B/D, C/D, and D/D, respectively.

We briefly describe uncertainty assessments of the SoilGrids predictions (soilgrids.org)¹⁵ and groundwater table depth¹⁶ data that were used as input for our analysis; however, readers are referred to the corresponding publications for a detailed description of the methods and uncertainty analysis.

SoilGrids

Soil profile data was compiled by the FAO from approximately 150,000 unique sites covering every continent; however, the tropics, semi-arid to hyper-arid regions, and mountain regions were underrepresented¹⁵. Furthermore, soils with high runoff potential are likely under-estimated due to the uncertainty associated with depth to bedrock¹⁵. However, their depth to bedrock models performed reasonably well, and explained more than 50% of the global variation (R²=0.54).

Accuracy assessment was performed with 10-fold repeated cross-validation using soil profile data from ca. 150 000 globally distributed sites used to develop soilGrids250m¹⁵. In all instances, the amount of variation explained by the soil texture models was higher than 72.6%; root mean square error (RMSE) was lowest for clay (9.5%), followed by silt (9.8%), and sand (13.1%)¹⁵.

Groundwater table depth

A total of 1,603,781 well sites were compiled from government archives and published literature to generate predictions of global groundwater table depth¹⁶. On average, the modeled groundwater table was 1.62 m (±17.91 m) lower than observations at the global scale. Note that local, perched aquifers were not modeled¹⁶. Groundwater pumping, drainage, and irrigation were not represented, thus neglecting the local complexity of human influence and only capturing the broad-scale patterns of groundwater¹⁶.

Comparison with other datasets

Hong and Adler¹⁹ reported that the global distribution of soils was dominated by moderately low runoff potential (36.8%), followed by high (25.3%), low (20.5%), and moderately high (17.4%) runoff potential. Although this is in stark contrast with what we report, these discrepancies are largely attributed to different classification schemes (Table 1), and to a lesser extent, different methodologies.

For comparative purposes only, we used the same classification scheme reported by Hong and Adler^12,19. This comparison revealed that the distribution of the two datasets were in closer agreement, and that soils are dominated by moderately low runoff potential (37%), followed by high (32%), low (17%), and moderately high (15%) runoff potential. However, it is important to note that the classification scheme reported by Hong and Adler was based on earlier work by Musgrave¹³ using rainfall, runoff, and infiltrometer measurements¹³, a practice that has since been abandoned by the USDA¹¹. Furthermore, the deprecated classification scheme does not account for the presence of impermeable layers (e.g., bedrock) or depth to groundwater table.

Other considerations

Note that substantial variation can exist within and between soil texture classes and their respective hydraulic properties (Fig. 6). According to the revised NEH¹¹, HSG-A typically consists soils classified as sand (e.g, more than 90% sand and less than 10% clay content), but can include loamy sand, sandy loam, loam, or silt loam. Likewise, HSG-B typically consists of loamy sand and sandy loam, but can contain loam, silt loam, silt, or sandy clay loam, while HSG-C typically consists of loam, silt loam, sandy clay loam, clay loam, and silty clay loam, but can include clay, silty clay, and sandy clay textures¹¹.

(a) Southeast US, (b) Southeast Asia, (c) northern South America, (d) Australia.

Full size image

Users of this dataset should be aware that HYSOGs250m represents general patterns of soil runoff potential appropriate for regional- to global-scale analyses and may not capture the local variance suitable for fine-scale applications. Although originally developed to support CN-based computations of rainfall runoff, HYSOGs250m can be used as a covariate for empirical analyses investigating various soil-environmental relationships. For example, plant and/or animal species distributions are often related to soil texture, plant available water, and groundwater. HYSOGs250m may be a useful covariate to further explain such relationships, as these data were produced by incorporating depth to bedrock, depth to groundwater table, and soil texture classes. These data can also be used for flood risk assessment and suitability mapping. End-users who are not interested in dual HSGs may simply re-classify HSGs A/D, B/D, C/D, and D/D to HSG-D.

How to cite this article: Ross, C. W. et al. HYSOGs250m, global gridded hydrologic soil groups for curve-number-based runoff modeling. Sci. Data 5:180091 doi: 10.1038/sdata.2018.91 (2018).

Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.