6 Sensitivity Index (SI)

Calculating a sensitivity index allows the impact of changing an independent variable (e.g. a model parameter) to be assessed on a dependent variable. Another way to think of this is how much does an output change if I change an input? Given all other variables are held constant, the SI is computed by dividing the percentage change in the simulated result (i.e. the output) by the percentage change of the input parameter:

\[\begin{equation} SI = \frac{(Output_{new}-Output_{original})/Output_{original}}{(Parameter_{new}-Parameter_{original})/Parameter_{original}} \tag{6.1} \end{equation}\]

If there is a large change in the output variable, we say that the variable is sensitive to the parameter. If there is a small change in the output variable, we say the variable is not sensitive to the parameter.

In this exercise we will compare the sensitivity of the water temperature and algae concentration to the mean wind speed, and the water clarity.

Calculate the sensitivity of the modelled temperature to changes in water clarity (the light extinction coefficient, \(K_w\)) and wind speed (wind_factor). These can be found in glm3.nml: &light and &meteorology, respectively.

Try increasing and decreasing the default parameter value by 0.2 (Kw) and 0.3 (wind) and see how much the output changes.

Using the SI calculation on the simulation output, assess how sensitive the water temperature and phytoplankton biomass is to water clarity (Kw):

GLM results Water clarity Decrease 0.37 Water clarity Original 0.57 Water clarity Increase 0.77
Average WQ_35 temperature
Average WQ_35 phytoplankton (green, crypto, diatom) biomass


Assess how sensitive the temperature and phytoplankton biomass is to wind speed (wind_factor):


GLM results Wind Speed Decrease 0.6 Wind Speed Original 0.9 Wind Speed Increase 1.2
Average WQ_35 temperature
Average WQ_35 phytoplankton (green, crypto, diatom) biomass

Are the water variables more sensitive to the change in wind, or the change water clarity?