# USLE output generated from InVEST SDR

**URL:** <https://community.naturalcapitalalliance.org/t/usle-output-generated-from-invest-sdr/1820>\
**Category:** Software Support\
**Tags:** sdr\
**Created:** [April 5, 2021, 8:20pm UTC](https://community.naturalcapitalalliance.org/t/usle-output-generated-from-invest-sdr/1820 "2021-04-05T20:20:58Z")\
**Posts on this page:** 8\
**Page:** 1

<div class="post-metadata">

**Author:** ![Habtamu\_Assaye](https://avatars.discourse-cdn.com/v4/letter/h/b2d939/32.png) [@Habtamu\_Assaye](https://community.naturalcapitalalliance.org/u/Habtamu_Assaye)\
**Post date:** [April 5, 2021, 8:20pm UTC](https://community.naturalcapitalalliance.org/t/usle-output-generated-from-invest-sdr/1820/1 "2021-04-05T20:20:58Z")

</div>

I have found the usle output generated from InVEST SDR model to be smaller than the one that I run directly in ArcGIS. Would there be any possible reason?

---

<div class="post-metadata">

**Author:** ![esoth](https://avatars.discourse-cdn.com/v4/letter/e/9d8465/32.png) [@esoth](https://community.naturalcapitalalliance.org/u/esoth)\
**Post date:** [April 5, 2021, 9:15pm UTC](https://community.naturalcapitalalliance.org/t/usle-output-generated-from-invest-sdr/1820/2 "2021-04-05T21:15:41Z")

</div>

Hi @Habtamu_Assaye,  
Do you mean the raster extent is smaller, or the values are smaller? What is your ArcGIS method for calculating USLE? Could you please upload the two versions? Thanks.

---

<div class="post-metadata">

**Author:** ![Habtamu\_Assaye](https://avatars.discourse-cdn.com/v4/letter/h/b2d939/32.png) [@Habtamu\_Assaye](https://community.naturalcapitalalliance.org/u/Habtamu_Assaye)\
**Post date:** [April 5, 2021, 11:51pm UTC](https://community.naturalcapitalalliance.org/t/usle-output-generated-from-invest-sdr/1820/3 "2021-04-05T23:51:01Z")

</div>

Dear esoth,  
I have run the RUSLE using the raster calculator in ArcGIS using same input rasters of R, K, LS, C and P that I used for InVEST. However the usle output from InVEST was quite different from the one I found from calculated using the raster calculator. I have attache the two rasters. usle\_SDR3 is the InVEST SDR output and A3 is using rastrer calculator. I was expecting similar result but there is big difference. I can send you all other files if needed.  
Thank you for your support.  
[usle\_SDR3.tif](https://community.naturalcapitalalliance.org/uploads/short-url/ty4ymxs1WS2SBj4LScVorGYKw3o.tif) (146.7 KB)

---

<div class="post-metadata">

**Author:** ![swolny](https://yyz1.discourse-cdn.com/flex031/user_avatar/community.naturalcapitalalliance.org/swolny/32/181_2.png) [@swolny](https://community.naturalcapitalalliance.org/u/swolny)\
**Post date:** [April 6, 2021, 2:51am UTC](https://community.naturalcapitalalliance.org/t/usle-output-generated-from-invest-sdr/1820/4 "2021-04-06T02:51:59Z")

</div>

Hi @Habtamu_Assaye -

It looks like only one of the rasters was attached (usle\_SDR3.tif), would you try to upload the other one again? Thanks!

When you’re doing the RUSLE calculation using raster calculator, are you using the [equations given in the InVEST User Guide](https://storage.googleapis.com/releases.naturalcapitalproject.org/invest-userguide/latest/sdr.html#annual-soil-loss), or different equations? In particular, we do the LS calculation a little differently. I don’t know if this explains your results, but might.

~ Stacie

---

<div class="post-metadata">

**Author:** ![Habtamu\_Assaye](https://avatars.discourse-cdn.com/v4/letter/h/b2d939/32.png) [@Habtamu\_Assaye](https://community.naturalcapitalalliance.org/u/Habtamu_Assaye)\
**Post date:** [April 6, 2021, 7:53am UTC](https://community.naturalcapitalalliance.org/t/usle-output-generated-from-invest-sdr/1820/5 "2021-04-06T07:53:13Z")

</div>

Hi@swolny,  
I am sorry it is not uploading the other raster, a bit strange. I used the Dezmet and Govers 1996 method of LS factor determination which is integrated in the SAGA gis within the qgis environment which is similar to what InVEST SDR uses. My e-mail is [habtamuasaye@gmail.com](mailto:habtamuasaye@gmail.com) and I can send you all the factor files if you can help me figure it out. Thank so much indeed!

---

<div class="post-metadata">

**Author:** ![esoth](https://avatars.discourse-cdn.com/v4/letter/e/9d8465/32.png) [@esoth](https://community.naturalcapitalalliance.org/u/esoth)\
**Post date:** [April 6, 2021, 4:08pm UTC](https://community.naturalcapitalalliance.org/t/usle-output-generated-from-invest-sdr/1820/6 "2021-04-06T16:08:54Z")

</div>

Hi @Habtamu_Assaye,  
It would be helpful if you can send the R, K, LS, C, P files and the ArcGIS output as well so we can check if it’s working as expected. You can email them to me: esoth @ stanford.edu Thanks!

---

<div class="post-metadata">

**Author:** ![esoth](https://avatars.discourse-cdn.com/v4/letter/e/9d8465/32.png) [@esoth](https://community.naturalcapitalalliance.org/u/esoth)\
**Post date:** [April 12, 2021, 10:38pm UTC](https://community.naturalcapitalalliance.org/t/usle-output-generated-from-invest-sdr/1820/7 "2021-04-12T22:38:08Z")

</div>

Hi @Habtamu_Assaye,  
Thank you for sharing your data with me! I’m going to reply here in the forum in case this is useful to others. The output from InVEST looks correct to me.  
There is a small mistake in the user’s guide that might have led you astray:

 ![image](https://canada1.discourse-cdn.com/flex031/uploads/naturalcapitalproject/original/2X/5/599aad4bfa00977258968e4ec315c190e95df2e2.png)  
It incorrectly compares θ to percent values, which really should be radians: arctan(9/100) ≈ 0.08975, so this should say θ \< 0.08975 and θ \>= 0.08975. I will update it today.

However I doubt that would cause the large difference in USLE you are seeing. I assume you are using the same USLE equation in ArcGIS: `USLE = R x K x LS x C x P`.

- For the `LS`, `C`, and `P` inputs in ArcGIS, are you using those intermediate outputs from InVEST (`intermediate/ls.tif`, `intermediate/cp.tif`), or are you calculating them independently?
- If you are calculating `LS` independently, are you sure you’re using the same method (described [here](https://storage.googleapis.com/releases.naturalcapitalproject.org/invest-userguide/latest/sdr.html#annual-soil-loss))? There are some different formulas for `LS` that could account for the difference.

Thanks for your patience working this out!

---

<div class="post-metadata">

**Author:** ![Habtamu\_Assaye](https://avatars.discourse-cdn.com/v4/letter/h/b2d939/32.png) [@Habtamu\_Assaye](https://community.naturalcapitalalliance.org/u/Habtamu_Assaye)\
**Post date:** [April 13, 2021, 8:43am UTC](https://community.naturalcapitalalliance.org/t/usle-output-generated-from-invest-sdr/1820/8 "2021-04-13T08:43:23Z")

</div>

Dear Emily,  
Thank you for your detailed response. I calculated LS, C and P independently. For the LS, I used Desmet and Govers 1996, directly from SAGA gis integrated with the Qgis environment which uses similar equations described in the InVEST guide. I will check with the updated version of the program and come back if I will still have further questions. Thank you so much indeed for your time, willingness and support!

With best regards,
