The PROWQM physical-biological model with benthic-pelagic coupling applied to the northern North Sea

PROWQM, a 1-D depth resolving model which couples physical and microbiological processes in the water column with sedimentation/resuspension and benthic mineralisation processes, has been used to simulate seasonal changes of chlorophyll, nutrients and oxygen at the PROVESS north site (59degrees20' N 1 degrees00' E) in the North Sea. PROWQM is derived from the 3-D model COHERENS, and improves COHEREN's benthic and pelagic biology. The physical sub-model of PROWQM implicitly solves turbulence closure equations forced by climatological, or realistic high-frequency, meteorological and tidal data. The pelagic biological sub-model 2MPPD includes a 'diatomy' microplankton (mp1) and a 'flagellatey' (or microbial loop) microplankton (mp2), the cycling of silicon and nitrogen, slow-sinking detritus, and fast-sinking phytodetritus. Phytodetritus is formed by shear-driven aggregation of particulate material, using a simple algorithm. for bulk processes that is derived by considering the interactions of single cells. The microplankton compartments include heterotrophic bacteria and protozoa as well as phytoplankton, and most microplankton rates are specified with the aid of a 'heterotroph fraction' parameter, which was 0.125 for mp1 and 0.6 for mp2. The microbiological system is closed by mesozooplankton grazing pressures imposed as time varying series determined from observed zooplankton abundance. The benthic boundary sub-model includes a superficial fluff layer and a nutrient element reservoir in the consolidated sediment. Particulate material in the fluff layer can be resuspended (in response to bed stress by near-bed flows), mineralised or carried by bioturbation into the underlying, consolidated, sediment, where it is mineralised and its nutrients returned to the water-column at rates mainly dependent on (implicit) macrobenthic pumping. Benthic denitrification can occur when mineralisation rates exceed oxygen supply. Verification of the PROWQM numerical implementation used test cases and checks for nutrient element conservation. Simulations with realistic forcing, for a range of parameter values, were compared with historic observations in the NOWESP data set and during FLEX76, and with those made during the PROVESS cruises in autumn 1998. PROWQM provided a good simulation of the seasonal succession from a diatom-dominated spring bloom to summer dominance by small flagellates. The simulations included sedimentation of organic matter from the spring bloom, and qualitatively realistic behaviour of the fluff layer, but decay rates were too slow and there was almost no denitrification. The simulated surface mixed layer was too shallow during,the summer. Simulated annual net microplankton primary production was in between 59 and 91 g C m(-2) y(-1). A large proportion of mineralisation, 28-47% of nitrogen and 40-67% of silicon mineralisation, took place as a result of the decay of sinking and resuspended detritus whilst in the water column. PROWQM is discussed in relation to other models that have been used to simulate this part of the North Sea, in particular the simpler ECOHAM1 and the more complex ERSEM, and in relation to PROWQM's, evolution from COHERENS. (C) 2002 Elsevier Science B.V All rights reserved.