NDR model - Nutrient runoff proxy

Hello everyone,

I’m using the NDR model to assess nutrient export in relation to changes in precipitation.

I’m using a raster showing the average annual precipitation over my study area. This is my Nutrient runoff proxy. In QGIS, I created a new raster file that shows an increase of rainfall of 25%.

In other words, I ran the NDR model twice: using the “average annual precipitation” as proxy and the “average annual precipitation +25%” as another proxy.

The problem is that the results of nutrient exports are almost identical between the two runs. It seems to me that an increase in precipitation does not change nutrient export, probably because this proxy represents “the spatial variability in runoff potential”, but I wonder if I’m correct.

I utilized different DEMs with different resolutions, the coordinate systems are correct, and I also changed the Threshold flow accumulation values.

Thank you for your help.


Hi Marco -

First, just wondering if only the patterns look the same, or if the underlying values are also the same.

Did you increase all pixels in the precipitation raster by 25%? If so, I think you’re correct to note that the runoff proxy represents spatial variability. If you don’t have spatial variability in the change in precipitation, then the model results won’t change. This is related to what the User Guide says about the calculation of RPI:

RPIi is the runoff potential index on pixel i. It is defined as:
RPIi = RPi / RPav
where RPIi is the nutrient runoff proxy for runoff on pixel i, and RPav is the average RP over the raster.

So since it’s creating an index in this way, if everything changes by 25%, it seems like you’d still end up with the same index. You can check this by looking at the model output intermediate_outputs/runoff_proxy_index.tif.

~ Stacie

1 Like