What causes short-term oscillations in ecosystem models of the ocean mixed layer?

Abstract An ecosystem model of the oceanic upper mixed layer (UML) that simulates the differences in the annual production cycle observed between the North Atlantic and North Pacific is proposed. It is based on the model of Fasham (1993) modified as follows: (1) the diffusive mixing between the UML and the seasonal pycnocline that is due to internal wave breaking is not included; (2) the nitrate concentration in the seasonal pycnocline is considered to vary with depth rather than be constant; (3) an alternative parameterisation of Zooplankton grazing is used and (4) a number of model parameters are assigned other values. The modelled and observed ecosystem characteristics at Ocean Weather Stations (OWS) I and P are in satisfactory agreement with each other. The model solution obtained shows different regimes of seasonal variability: a periodic regime having a simple annual oscillation and a so-called biharmonic regime that also has oscillations with a period of order one month. The latter oscillations are of interest in that there are observational data showing a succession of phytoplankton blooms at OWS I, although the data are not of sufficient quality to show that the mechanisms causing the oscillations in the model are the same as those producing the oscillations in the data. A qualitative mathematical analysis of the model equations with stationary forcing shows that in the case of entrainment the sole equilibrium solution with non-zero concentrations of all of seven ecosystem components exists, while in the case of detrainment no such solution exists. Analysis of the system of the equations linearised in the vicinity of the former equilibrium state shows it to be stable or unstable and usually corresponds to monotonic changes or damped oscillations at periodic forcing. The model short-term oscillations are connected with feedbacks in the ecosystem initiated by either abrupt changes of the UML thickness in spring or strong deepening of the UML in autumn. A necessary condition of their existence is a combination of intense photosynthetically active radiation, small thickness of the UML, significant mean annual entrainment velocity and a high nitrogen concentration in the seasonal pycnocline. Such combinations are uncommon in the world ocean and this is the reason why such oscillations of the UML ecosystem are only to be expected in some regions.

[1]  H. Ducklow,et al.  Plankton succession and carbon cycling at 47°N 20°W during the JGOFS North Atlantic Bloom Experiment , 1993 .

[2]  P. Hartman.  Ordinary Differential Equations , 1965 .

[3]  R. Goericke,et al.  Primary production in the subarctic Pacific Ocean: Project SUPER , 1993 .

[4]  C. L. Leonard,et al.  Coastal zone color scanner pigment concentrations in the Southern Ocean and relationships to geophysical surface features , 1993 .

[5]  P. Burkill,et al.  Microzooplankton and their herbivorous activity in the northeastern Atlantic Ocean , 1993 .

[6]  C. S. Wong.  Analysis of trends in primary productivity and chlorophyll-a over two decades at Ocean Station P (50°N, 145°W) in the subarctic northeast Pacific Ocean , 1995 .

[7]  Thomas M. Powell,et al.  Ecological dynamics in the subarctic Pacific, a possibly iron-limited ecosystem , 1991 .

[8]  Donald L. DeAngelis,et al.  The global carbon cycle. , 1990 .

[9]  B. Frost.  The role of grazing in nutrient-rich areas of the open sea , 1991 .

[10]  W. Richard,et al.  TEMPERATURE AND PHYTOPLANKTON GROWTH IN THE SEA , 1972 .

[11]  S. Levitus.  Climatological Atlas of the World Ocean , 1982 .

[12]  K. Banse.  Grazing, Temporal Changes of Phytoplankton Concentrations, and the Microbial Loop in the Open Sea , 1992 .

[13]  H. Ducklow,et al.  A nitrogen-based model of plankton dynamics in the oceanic mixed layer , 1990 .

[14]  G. Evans,et al.  Towards a model of ocean biogeochemical processes , 1993 .

[15]  M. Fasham,et al.  Modelling the Marine Biota , 1993 .

[16]  John H. Steele,et al.  The Significance of Interannual Variability , 1993 .

[17]  I. Totterdell.  An Annotated Bibliography of Marine Biological Models , 1993 .

[18]  J. Steele,et al.  The role of predation in plankton models , 1992 .

[19]  Paul G. Falkowski,et al.  Primary Productivity and Biogeochemical Cycles in the Sea , 1992 .

[20]  R. Weller,et al.  Bio‐optical variability associated with phytoplankton dynamics in the North Atlantic Ocean during spring and summer of 1991 , 1995 .

[21]  William H. Press,et al.  Numerical Recipes: FORTRAN , 1988 .

[22]  J. Sarmiento,et al.  Ecosystem behavior at Bermuda Station “S” and ocean weather station “India”: A general circulation model and observational analysis , 1993 .

[23]  G. Evans,et al.  The Use of Optimization Techniques to Model Marine Ecosystem Dynamics at the JGOFS Station at 47 degrees N 20 degrees W [and Discussion] , 1995 .

[24]  B. Frost.  A modelling study of processes regulating plankton standing stock and production in the open Subarctic Pacific Ocean , 1993 .

[25]  C. Carlson,et al.  Nitrate utilization in surface waters of the Iceland Basin during spring and summer of 1989 , 1993 .