ZoopModelSizeSpectra/ZooMSS_User_Manual.Rmd at master · MathMarEcol/ZoopModelSizeSpectra · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
---
title: "ZooMSS User Manual"
author: "The UQ ZooMSS Team"
date: "18/03/2020"
output: pdf_document
---

```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE)
```

Note: This is still very much a work in progress.

## Overview of ZooMSS

The Zooplankton Model of Size Spectra (ZooMSS) is a functional size-spectrum model of the marine ecosystem (following Heneghan et al. 2016, ZooMSSv1) to resolve phytoplankton, nine zooplankton functional groups (heterotrophic flagellates and ciliates, omnivorous and carnivorous copepods, larvaceans, euphausiids, salps, chaetognaths and jellyfish) and three size-based fish groups. Zooplankton functional groups are resolved using their body-size ranges, size- based feeding characteristics and carbon content, and the zooplankton community emerges from the model across global environmental gradients, depending on the functional traits of the different groups.

We developed the Zooplankton Model of Size Spectra version 2 (ZooMSSv2) based on the prototype of Heneghan et al. (2016). ZooMSSv2 uses the functional size-spectrum framework (Blanchard et al., 2017) to resolve the body size ranges, size-based feeding characteristics and carbon content of nine zooplankton groups and three fish groups. The model is run on 5° grid squares across the global ocean. For each region, the model is forced with the long- term mean satellite sea surface temperature and chlorophyll a from MODIS-Aqua.

ZooMSSv2 represents the marine ecosystem as three communities: phytoplankton, zooplankton and fish. The zooplankton community consists of nine of the most abundant zooplankton groups, and the fish community was made up of a small, medium and large group. Dynamics of the phytoplankton are not explicitly resolved in the model, rather the mean size structure of the phytoplankton community is estimated directly from satellite chlorophyll a observations (Brewin et al., 2010; Barnes et al., 2011; Hirata et al., 2011). Abundances of the zooplankton and fish communities are driven by size-dependent processes of growth and mortality, with the temporal dynamics of each functional group governed by separate second- order McKendrick-von Foerster equations.


## How to run ZooMSS using `R`

### Model Structure

The model is split into 5 main files.

1. Setup_RUN_NAME.R - Set up your own specific simulation. This is the only file you usually need to edit.

2. fZooMSS_Model.R - The top-level model function which seuqentially runs the following 3 functions, and saves the output.

    + fZooMSS_Params.R - Set all the parameter values for the model.

    + fZooMSS_Setup.R - Do all calculations which can be done in advance of the time loop.

    + fZooMSS_Run.R - Contains the time-dependent loop.

### Input files

ZooMSS requires two input files:

1. TestGroups.csv which contains all the taxa-specific parameter values, including size ranges.

2. And RDS file which contains all environmental data for each grid cell. This includes chlorophyll a (chlo), sea surface temperature (sst), bathymetry (bathy), phytoplankton slope and intercept, and latitude/longitude which is useful for plotting.


### Parameter Values

Insert table from Heneghan et al (in review)

### Output

TBC

### Running
To run the model, you simply need to provide the relevent information in Setup_RUN_NAME.R and run the script.

### Plotting

TBC

## Useful reading

Datta Paper
The other thing to look at would be some of the numerical methods we use
  - Simpson’s method for the growth, mortality and diffusion integrals and
  - a semi-implicit upward finite difference scheme to solve the PDE at each time step.
  For the PDE numerical method, check out Appendix 1 of my thesis (attached)
  Info on the difference scheme in Press’ 2007 book: Numerical Recipes, the Art of Scientific Computing (Available at UQ Library)


On 1 Apr 2020, at 06:50, Richardson, Anthony (O&A, St. Lucia) <Anthony.Richardson@csiro.au> wrote:

<Datta et al. - 2010 - A Jump-Growth M odel for Predator–Prey Dynamics De.pdf>

## Publications

1. Heneghan, R.F., Everett, J.D., Blanchard, J.L., Richardson, A.J., 2016. Zooplankton Are Not Fish: Improving Zooplankton Realism in Size-Spectrum Models Mediates Energy Transfer in Food Webs. Front. Mar. Sci. 3, 1–15. https://doi.org/10.3389/fmars.2016.00201

2. Heneghan, R.F., Everett, J.D., Sykes, P., Batten, S.D., Edwards, M., Takahashi, K., Suthers, I.M.,  Blanchard, J.L., Richardson, A.J., in review, A global size-spectrum model of the marine ecosystem that resolves zooplankton composition. Ecological Modelling

## References

TBC

## Getting Help

If you are running the model come across a problem, raise an issue on GitHub: https://github.com/MathMarEcol/ZooMSS/issues

If you find some errors, or just want to improve the model, we'd love you to go ahead and make the changes and then submit a pull request for us to check and approve.

Same for this manual. If there is information lacking, errors etc, please fix it rather than just reporting it. This should be a team effort.