Parameter uncertainty of a dynamic multispecies size spectrum models

Dynamic size spectrum models have been recognized as an effective way of describing how size-based interactions can give rise to the size structure of aquatic communities. They are intermediate-complexity ecological models that are solutions to partial differential equations driven by the size-dependent processes of predation, growth, mortality, and reproduction in a community of interacting species and sizes. To be useful for quantitative fisheries management these models need to be developed further in a formal statistical framework. Previous work has used time-averaged data to calibrate the model using optimization methods with the disadvantage of losing detailed time-series information. Using a published multispecies size spectrum model parameterized for the North Sea comprising 12 interacting fish species and a background resource, we fit the model to time-series data using a Bayesian framework for the first time. We capture the 1967-2010 period using annual estimates of fishing mortality rates as input to the model and time series of fisheries landings data to fit the model to output. We estimate 38 key parameters representing the carrying capacity of each species and background resource, as well as initial inputs of the dynamical system and errors on the model output. We then forecast the model forward to evaluate how uncertainty propagates through to population-and community-level indicators under alternative management strategies.