The reduction of plankton biomass induced by mesoscale stirring: A modeling study in the Benguela upwelling
Introduction
Marine ecosystems of the eastern boundary upwelling zones are well-known for their major contribution to the world ocean productivity. They are characterized by wind-driven upwelling of cold nutrient-rich waters along the coast that supports elevated plankton and pelagic fish production (Mackas et al., 2006). Variability is introduced by strong advection along the shore, physical forcings by local and large scales winds, and high submeso- and mesoscale activities over the continental shelf and beyond, linking the coastal domain with the open ocean.
The Benguela Upwelling System (BUS) is one of the four major Eastern Boundary Upwelling Systems (EBUS) of the world. The coastal area of the Benguela ecosystem extends from southern Angola (around ) along the west coast of Namibia and South Africa (). It is surrounded by two boundary currents, the warm Angola current in the north and the temperate Agulhas current in the south. The BUS can itself be subdivided into two subdomains by the powerful Luderitz upwelling cell (Hutchings et al., 2009). Most of the biogeochemical activity occur within the upwelling front and the coast, although it can be extended further offshore toward the open ocean by the numerous filamental structures developed in offshore (Monteiro, 2009). In the BUS, as in the other major upwelling areas, high mesoscale activity due to eddies and filaments impacts strongly a marine planktonic ecosystem over the continental shelf and beyond (Brink and Cowles, 1991, Martin, 2003, Sandulescu et al., 2008, Rossi et al., 2009).
The purpose of this study is to analyze the impact of horizontal stirring on phytoplankton dynamics in the BUS within an idealized two dimensional modeling framework. Based on satellite data of the ocean surface, Rossi et al., 2008, Rossi et al., 2009 recently suggested that mesoscale activity has a negative effect on chlorophyll standing stocks in the four EBUS. This was obtained by correlating remote sensed chlorophyll data with a Lagrangian measurement of lateral stirring in the surface ocean (see Methods section). This result was unexpected since mesoscale physical structures, particularly mesoscale eddies, have been related to higher planktonic production and stocks in the open ocean (McGillicuddy et al., 2007) as well as off a major EBUS (Correa-Ramirez et al., 2007). A more recent and thorough study performed by Gruber et al. (2011) in the California and the Canary current systems extended the initial results from Rossi et al., 2008, Rossi et al., 2009. Based on satellite derived estimates of net primary production, of upwelling strength and of Eddy Kinetic Energy (EKE) as a measure of the intensity of mesoscale activity, they confirmed the suppressive effect of mesoscale structures on biological production in upwelling areas. Investigating the mechanism behind this observation by means of on 3D eddy-resolving coupled models, Gruber et al. (2011) showed that mesoscale eddies tend to export offshore and downward a certain pool of nutrients not being effectively used by the biology in the coastal areas. The process they called “nutrients leakage” is also having a negative feedback by diminishing the pool of deep nutrients available in the surface waters being re-upwelled continuously.
In our work, we focused on the Benguela area, being the most contrasting area of all EBUS in terms of stirring intensity (Rossi et al., 2009). Although the mechanisms studied by Gruber et al. (2011) seem to involve 3D dynamics, the initial observation of this suppressive effect was essentially based on two-dimensional (2D) datasets (Rossi et al., 2008). In this work we use 2D numerical analysis in a semi-realistic framework to better understand the effects of a 2D turbulent flow on biological dynamics, apart from the complex 3D bio-physical processes. The choice of this simple horizontal numerical approach is indeed supported by other theoretical 2D studies that also displayed a negative correlation between stirring and biomass (Tél et al., 2005, MacKiver and Neufeld, 2009, Neufeld and Hernández-García, 2009). Meanwhile, since biological productivity in upwelling areas relies on the (wind-driven) vertical uplift of nutrients, we introduced in our model a nutrient source term with an intensity and spatial distribution corresponding to the upwelling characteristics. Instead of the commonly used EKE, which is an Eulerian diagnostic tool, we used here a Lagrangian measurement of mesoscale stirring that has been demonstrated as a powerful tool to study patchy chlorophyll distributions influenced by dynamical structures at mesoscale, such as upwelling filaments (Calil and Richards, 2010). The Lagrangian perspective provides a complementary insight to transport phenomena in the ocean with respect to the Eulerian one. In particular, the concept of Lagrangian Coherent Structure may give a global idea of transport in a given area, separating regions with different dynamical behavior, and signaling avenues and barriers to transport, which are of great relevance for the marine biological dynamics. While the Eulerian approach describes the characteristics of the velocity field, the Lagrangian one addresses the effects of this field on transported substances, which is clearly more directly related to the biological dynamics. For example the work by Hernández-Carrasco et al. (2012) describes currents in the world ocean having the same level of Eddy Kinetic Energy but having two different stirring characteristics, as quantified by Lagrangian tools. Further discussions comparing Lagrangian and Eulerian diagnostics can be found, for example, in d'Ovidio et al. (2009) and the above cited Hernández-Carrasco et al. (2012). To consider velocity fields with different characteristics and to test the effect of the spatial resolution, different flow fields are used, one derived from satellite and two produced by numerical simulations at two different spatial resolutions. Our modeled chlorophyll a concentrations are compared with observed distributions of chlorophyll a (a metric for phytoplankton) obtained from the SeaWiFS satellite sensor.
This paper is organized as follows. Section 2 is a brief description of the different datasets used in this study. Section 3 depicts the methodology, which includes the computation of the finite-size Lyapunov exponents, and the numerical plankton-flow 2D coupled model. Then, our results are analyzed and discussed in Section 4. Finally in Section 5, we summed-up our main findings.
Section snippets
Satellite and simulated data
We used three different 2D surface velocity fields of the Benguela area. Two are obtained from the numerical model Regional Ocean Model System (ROMS), and the other one from a combined satellite product.
Finite-Size Lyapunov Exponents (FSLEs)
FSLEs (Artale et al., 1997, Aurell et al., 1997, Boffetta et al., 2001) provide a measure of dispersion, and thus of stirring and mixing, as a function of the spatial resolution. This Lagrangian tool allows isolating the different regimes corresponding to different length scales of the oceanic flows, as well as identifying Lagrangian Coherent Structures (LCSs) present in the data (Tew Kai et al., 2009). FSLEs are computed from , the time required for two particles of fluid (one of them placed
Horizontal stirring
We compute the FSLE with an initial separation of particles equal to the spatial resolution of each velocity field ( for Satellite1/4 and ROMS1/4, and for ROMS1/12), and a final distance of to focus on transport processes by mesoscale structures at mid latitudes. The areas of more intense horizontal stirring due to mesoscale activity can be identified by large values of temporal averages of backward FSLEs (see Fig. 2). While there are visible differences between the
Conclusions
We have studied the biological dynamics in the Benguela area by considering a simple biological NPZ model coupled with different velocity fields (satellite and model). Although in a simple framework, a reduction of phytoplankton concentrations in the coastal upwelling for increasing mesoscale activity has been successfully simulated. Horizontal stirring was estimated by computing the FSLEs and was correlated negatively with chlorophyll stocks. Similar correlations are found, though not
Acknowledgments
I.H.-C. was supported by a FPI grant from MINECO to visit LEGOS. We acknowledge support from MINECO and FEDER through projects FISICOS (FIS2007-60327) and ESCOLA (CTM2012-39025-C02-01). V.G. thanks CNES funding through Hiresubcolor project. We are also grateful to J. Sudre for providing us velocity datasets both from ROMS and from the combined satellite product. Ocean color data were produced by the SeaWiFS project at GES and were obtained from DAAC.
References (50)
- et al.
Detecting barriers to transporta review of different techniques
Physica D
(2001) - et al.
A new coupled, one-dimensional biological physical model for the upper ocean application to the JGOFS Bermuda Atlantic Time-series Study (BATS) site
Deep Sea Research II
(1996) - et al.
Comparison between Eulerian diagnostics and finite-size Lyapunov exponents computed from altimetry in the Algerian basin
Deep Sea Research I
(2009) - et al.
Lagrangian coherent structures and mixing in two-dimensional turbulence
Physica D
(2000) - et al.
How reliable are finite-size Lyapunov exponents for the assessment of ocean dynamics?
Ocean Modelling
(2011) - et al.
A pelagic ecosystem model calibrated with BATS data
Deep Sea Research II
(1996) - et al.
The Benguela currentan ecosystem of four components
Progress in Oceanography
(2009) Phytoplankton patchinessthe role of lateral stirring and mixing
Progress in Oceanography
(2003)- et al.
Assessment of satellite ocean color products at a coastal site
Remote Sensing Environment
(2007) Oceanic diffusion diagrams
Deep Sea Research
(1971)
Biological activity in the wake of an island close to a coastal upwelling
Ecological Complexity
Chemical and biological activity in open flowsa dynamical system approach
Physics Reports
The Benguelaa laboratory for a comparative modeling studies
Progress in Oceanography
Dispersion of passive tracers in closed basinsbeyond the diffusion coefficient
Physics of Fluids
Predictability in the largean extension of the Lyapunov exponent
Journal of Physics A
The coastal transition zone program
Journal of Geophysical Research
Transient upwelling hot spots in the oligotrophic North Pacific
Journal of Geophysical Research
Estimation of potential productivity in Eastern Boundary Currents using remote sensing
Deep Sea Research II
Seasonal characteristics of the surface mixed layer in the Australasian regionimplications for primary production regimes and biogeography
Marine and Freshwater Research
Mesoscale eddies and high chorophyll concentrations off central Chile ()
Geophysical Research Letters
Climatology and variability of sea surface temperature and surface chlorophyll in the Benguela and Agulhas ecosystems as observed by satellite
African Journal of Marine Science
Mixing structures in the Mediterranean sea from finite-size Lyapunov exponents
Geophysical Research Letters
Eddy-induced reduction of biological production in eastern boundary upwelling systems
Nature Geoscience
Coupled physical/biogeochemical modeling including 02-dependent processes in the eastern boundary upwelling systemsapplication in the Benguela
Biogeosciences
Cited by (37)
Lagrangian Coherent Structures in the Mediterranean Sea: Seasonality and basin regimes
2023, Progress in OceanographyModelling coastal upwelling off NW Iberian Peninsula: New insights on the fate of phytoplankton blooms
2023, Science of the Total EnvironmentEffects of the seasonality of mesoscale eddies on the planktonic dynamics off eastern Brazil
2022, Dynamics of Atmospheres and OceansCitation Excerpt :Indeed, offshore advection of shelf material caused by eddies is commonly observed in the world ocean. Other studies reported analogous processes: off southeastern Brazil (Yoshinaga et al., 2010); in the Mozambique Channel (Tew-Kai and Marsac, 2009); in the eastern Gulf of Alaska (Crawford et al., 2005); off Chile (Correa-Ramirez et al., 2007); and in the Benguela (Hernández-Carrasco et al., 2014) and the Canary upwelling systems (Lovecchio et al., 2017; Bonino et al., 2020). However, we should mention that not all the Chl-a variability is controlled by this process.
Preferential selection of zooplankton and emergence of spatiotemporal patterns in plankton population
2021, Chaos, Solitons and FractalsCitation Excerpt :Different biophysical properties of the aquatic body are considered to be the regulating factors of such spatial patterns. For example, non-diffusive advection [20], upwelling [21], turbulent stirring [22–24] and predation [25] have been shown as the causal factors of spatial heterogeneity of marine plankton. It is however unknown whether a preferential selection of zooplankton can create patchiness by the spectra of species concentrations.
Stability analysis of the plankton community with advection
2021, Chaos, Solitons and FractalsCitation Excerpt :In this article, we studied the influence of advection on the self-organization of plankton. Plankton drift with ocean currents, and the horizontal advection velocity of each plankton may be different due to their size, which makes a difference in the stability of the plankton community [37–39]. Through many aspects of comparison, we choose the Arditi-Ginzburg type-II functional response as the predation function [21] and consider allelopathic effect, the self-defense of phytoplankton.