nuts: Compute the Nutrient Balance and Ecosystem Services of...
In BenjaminDelory/heathland: Analysis of Heathland Management Trade-offs

Description Usage Arguments Details Value Author(s) See Also Examples

This function computes the cumulative balance of one essential nutrient (N, P, K, Mg, or Ca) or one ecosystem service (groundwater recharge, groundwater quality, carbon removal, cost, or appreciation) for many possible management scenarios.

1
2
3

nuts(scenarios, data=import_nutrients(), 
     sqa=list(grazing=1, mowing=5, burning=5, choppering=7),
     nutrient="N", service="waterrecharge", init=0)

`scenarios`	A `big.matrix` object created by `scenarios` or `import_scenarios`.
`data`	A `nutrients` object created by `import_nutrients` or a `services` object created by `import_services`.
`sqa`	A list giving the time period until leaching rates reach status quo ante (years) for each management.
`nutrient`	The nutrient to focus on. Only used if data is a nutrients object. Could be either `“N”`, `“P”`, `“K”`, `“Mg”`, or `“Ca”`.
`service`	The ecosystem service to focus on. Only used if data is a services object. Could be either `“waterrecharge”`, `“waterquality”`, `“carbon”`, `“cost”`, or `“appreciation”`.
`init`	The groundwater recharge/groundwater quality balance at the start of the simulation. Default to zero.

The following convention is used: 0 for no management, 1 for grazing, 2 for mowing, 3 for burning, and 4 for choppering.

Returns a big.matrix object containing the cumulative nutrient or ecosystem service balance for each scenario (row) each year (column). Values are given in kg/ha for the nutrient balance, mm/year for groundwater recharge, kg N/ha for groundwater quality, kg/ha for carbon removal, and euros/ha for cost. Appreciation is adimensional.

Benjamin M. Delory

import_nutrients, import_services, scenarios, import_scenarios.

#Create management scenarios
sc<-scenarios(management=c("none", "grazing", "mowing", "burning"), n=15)

#Compute the N balance for each scenario
balance<-nuts(sc)

#Plot the results
boxplot(balance[,ncol(balance)], ylab="Cumulative N balance (kg/ha)", las=1)
abline(h=0, lty=2)