[1] Oceanic plankton plays an important role in the marine food chain and through its significant contribution to the global carbon cycle can also influence the climate. Plankton bloom is a sudden rapid increase of the population. It occurs naturally in the North Atlantic as a result of seasonal changes. Ocean fertilization experiments have shown that supply of iron, an important trace element, can trigger a phytoplankton bloom in oceanic regions with low natural phytoplankton density. Here we use a simple mathematical model of the combined effects of stirring by ocean eddies and plankton evolution to consider the impact of a transient local perturbation, e.g., in the form of iron enrichment as in recent 'ocean fertilization' experiments. The model not only explains aspects of the bloom observed in such experiments but predicts the unexpected outcome of a large scale bloom that in its extent could be comparable to the spring bloom in the North Atlantic. [Truscott and Brindley, 1994a] suggested that plankton blooms might be explained by the excitability of the dynamical system describing biological interactions in the plankton ecosystem. The characteristic feature of excitable systems [Meron, 1992] is that perturbations exceeding a certain threshold can induce a temporary large deviation from the equilibrium state. The plankton bloom produced artificially by iron fertilization experiments has been shown to be consistent with excitable system behaviour [Pitchford and Brindley, 1999]. However, previous work has considered only the spatially homogeneous situation, and need not be relevant to the iron fertilization experiments, in which localised perturbations produce a localised bloom that is subsequently strongly affected by stirring by mesoscale ocean eddies.[3] The role of stirring has been recognised and investigated in the context of the SOIREE experiment [Boyd et al., 2000;Abraham et al., 2000] and the spatial characteristics of the bloom filament was found to be similar to that seen in tracer release experiments [Sundermeyer and Price, 1998], in having a roughly constant width and exponentially increasing length. However, the biological evolution in this analysis was modelled by a simple linear growth rate, which cannot be adequate for predicting the long term behavior of the bloom when nonlinearities, e.g. interactions with other species, become important.[4] Here we consider a, perhaps minimal, model that unifies the two approaches described above by taking into account spatial distribution, stirring and non-linear interactions between biological components. The model equations arewhere P, Z represent the phytoplankton and zooplankton concentration (in nitrogen currency) and F is the deviation of the iron concentration from the original homogeneous background. All three components are advected by a two-dimensional flow, v(r, t), obtained from a so-called seeded eddy model [Abraham, 1998], and diffusive terms represent the effect of small-scale mixing (k = 4 m 2 s À1 ). The biological terms on the right-han...