Beschreibung: A growing population on a planet with limited resources demands finding new sources of protein. Hence, fisheries are turning their perspectives towards mesopelagic fish, which have, so far, remained relatively unexploited and poorly studied. Large uncertainties are associated with regards to their biomass, turn-over rates, susceptibility to environmental forcing and ecological and biogeochemical role. Models are useful to disentangle sources of uncertainties and to understand the impact of different processes on the biomass. In this study, we employed two food-web models – OSMOSE and the model by Anderson et al. (2019, or A2019) – coupled to a regional physical–biogeochemical model to simulate mesopelagic fish in the Eastern Tropical South Pacific ocean. The model by A2019 produced the largest biomass estimate, 26 to 130% higher than OSMOSE depending on the mortality parameters used. However, OSMOSE was calibrated to match observations in the coastal region off Peru and its temporal variability is affected by an explicit life cycle and food web. In contrast, the model by A2019 is more convenient to perform uncertainty analysis and it can be easily coupled to a biogeochemical model to estimate mesopelagic fish biomass. However, it is based on a flow analysis that had been previously applied to estimate global biomass of mesopelagic fish but has never been calibrated for the Eastern Tropical South Pacific. Furthermore, it assumes a steady-state in the energy transfer between primary production and mesopelagic fish, which may be an oversimplification for this highly dynamic system. OSMOSE is convenient to understand the interactions of the ecosystem and how including different life stages affects the model response. The combined strengths of both models allow us to study mesopelagic fish from a holistic perspective, taking into account energy fluxes and biomass uncertainties based on primary production, as well as complex ecological interactions.

Beschreibung: In geoscience and other fields, researchers use models as a simplified representation of reality. The models include processes that often rely on uncertain parameters that reduce model performance in reflecting real-world processes. The problem is commonly addressed by adapting parameter values to reach a good match between model simulations and corresponding observations. Different optimization tools have been successfully applied to address this task of model calibration. However, seeking one best value for every single model parameter might not always be optimal. For example, if model equations integrate over multiple real-world processes which cannot be fully resolved, it might be preferable to consider associated model parameters as random parameters. In this paper, a random parameter is drawn from a wide probability distribution for every singe model simulation. We developed an optimization approach that allows us to declare certain parameters random while optimizing those that are assumed to take fixed values. We designed a corresponding variant of the well known Covariance Matrix Adaption Evolution Strategy (CMA-ES). The new algorithm was applied to a global biogeochemical circulation model to quantify the impact of zooplankton mortality on the underlying biogeochemistry. Compared to the deterministic CMA-ES, our new method converges to a solution that better suits the credible range of the corresponding random parameter with less computational effort.