Hi @michelle.vanhilten -
I looked up the layer you referenced from ArcGIS online. It sounds like it was derived from the same underlying data that I’ve used in the US (gNATSTO and gSSURGO). But, it also says this about units: “The values are in units of Micrometers per second (μm/s).” This is odd, since the description explicitly says that this layer is used in the USLE, but the USLE requires different units. And I’m really not sure how to convert micrometers per second to the complicated annual units (ton⋅ha⋅hr⋅(ha⋅MJ⋅mm)−1) that are required.
That said, I’m looking back at my notes about how I derived this layer recently in the US, and was reminded about the eternal problem of K factor units. Sometimes people use imperial units in place of metric units when they’re using the USLE and don’t seem to care. Sometimes no units are listed at all. Case in point, I have not yet found any documentation that explicitly says what the units of K factor are, as derived from either the old tool I used to use (Soil Data Viewer) or the new tool (Soil Data Development Toolbox), or in the underlying metadata for gSSURGO. And people use them all the time in USLE-based analyses. Based on the values the tools return, I have been assuming that the units are imperial, and convert them to metric accordingly. But this may be a naïve assumption!
If you are using the results in a relative way (identifying high/low values, not needing absolute values), then the units might not matter as much, since they still give an indication of whether the soil is more or less erodible. I’m sorry this all just muddies the waters, but it is a long-standing issue, and I don’t work in the US very often.
If you’re using ArcGIS, here are quick instructions I wrote up recently for how to use the Soil Data Development Toolbox, which is what I’d recommend if you want to derive this layer yourself.
gNATSGO from USDA:
Download gNATSGO data itself for whatever state(s) are needed (or CONUS, if working with a larger area, but it’s coarser resolution), from buttons on this page.
Also linked on this page, download “ArcTools for Working with gNATSGO”. Unzip the download, open an ArcMap session and add the Soil Data Development Toolbox to ArcToolbox.
Also linked on this page, look at the document “How to Create Soil Property or Interpretation Grid”. What follows is how I followed this doc to create the Kfactor.
From the gNATSGO geodatabase, add MapunitRaster_10m to the Arc session.
- Note that this is a huge file (32GB), but if you clip it before creating the soil map, you lose options in the tool interface, so it looks like you need to clip the result instead.
In the ArcGIS toolbox, open gSSURGO Mapping Toolset > Create Soil Map.
- Map Unit Layer = MapunitRaster_10m
- SDV Folder = Soil Erosion Factors
- SCV Attribute = K Factor, Whole Soil
- Aggregation Method = Dominant Condition (you can play with this)
- Top Depth = 0cm
- Bottom Depth = 5cm (play with this too, really we just need the surface, but it’s unclear how far down is recommended to consider)
- Leave the rest of the settings unchanged (many you won’t be able to change anyway)
- Run the tool
Export the resulting layer to a TIFF. Clip the result as needed.
The result will have an attribute table, including a Kfact column. Note that (I think - even the SSURGO table metadata doesn’t list units!) these values are in US units, they must be converted to metric for use in SDR. Thus sayeth the User Guide: “multiplication by 0.1317 is needed to convert from US customary units to ton⋅ha⋅hr⋅(ha⋅MJ⋅mm)−1, as detailed in Appendix A of the USDA RUSLE handbook (Renard et al., 1997).”