Sustenance of phytoplankton in the subpolar North Atlantic during the winter through patchiness

This study investigates the influence of two factors that change the mixed layer depth and can potentially contribute to the phytoplankton sustenance over winter: 1) variability of air-sea fluxes and 2) three-dimensional processes arising from strong fronts. To study the role of these factors, we perform several three-dimensional numerical simulations forced with air-sea fluxes at different temporal averaging frequencies as well as different spatial resolutions. Results show that in the winter, when the average mixed layer is much deeper than the euphotic layer and the days are short, phytoplankton production is relatively insensitive to the high-frequency variability in air-sea fluxes. The duration of upper ocean stratification due to high-frequency variability in air-sea fluxes is short and hence has a small impact on phytoplankton production. On the other hand, slumping of fronts creates patchy, stratified, shallow regions that persist considerably longer than stratification caused by changes in air-sea fluxes. Simulations show that before spring warming, the average MLD with fronts is about 700 m shallower than the average MLD without fronts. Therefore, fronts increase the residence time of phytoplankton in the euphotic layer and contribute to phytoplankton growth. Results show that before the spring warming, the depth-integrated phytoplankton concentration is about twice as large as phytoplankton concentration when there are no fronts. Hence, fronts are important for setting the MLD and sustaining phytoplankton in the winter. Model results also show that higher numerical resolution leads to stronger restratification, shallower mixed layers, greater variability in the MLD and higher production of phytoplankton

Coherent pathways for subduction from the surface mixed layer at ocean fronts

Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 126(5), (2021): e2020JC017042, https://doi.org/10.1029/2020JC017042.In frontal zones, water masses that are tens of kilometers in extent with origins in the mixed layer can be identified in the pycnocline for days to months. Here, we explore the pathways and mechanisms of subduction, the process by which water from the surface mixed layer makes its way into the pycnocline, using a submesoscale-resolving numerical model of a mesoscale front. By identifying Lagrangian trajectories of water parcels that exit the mixed layer, we study the evolution of dynamical properties from a statistical standpoint. Velocity- and buoyancy-gradients increase as water parcels experience both mesoscale (geostrophic) and submesoscale (ageostrophic) frontogenesis and subduct beneath the mixed layer into the stratified pycnocline along isopycnals that outcrop in the mixed layer. Subduction is transient and occurs in coherent regions along the front, the spatial and temporal scales of which influence the scales of the subducted water masses in the pycnocline. An examination of specific subduction events reveals a range of submesoscale features that support subduction. Contrary to the forced submesoscale processes that sequester low potential vorticity (PV) anomalies in the interior, we find that PV can be elevated in subducting water masses. The rate of subduction is of similar magnitude to previous studies (∼100 m/year), but the Lagrangian evolution of properties on water parcels and pathways that are unraveled in this study emphasize the role of submesoscale dynamics coupled with mesoscale frontogenesis.This research was funded by the ONR CALYPSO DRI grant N00014-16-1-3130. MAF was partially funded by a Martin Fellowship from MIT

Decomposition of vertical velocity for nutrient transport in the upper ocean

Author Posting. © American Meteorological Society, 2019. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 49(6), (2019): 1561-1575, doi:10.1175/JPO-D-19-0002.1.Within the pycnocline, where diapycnal mixing is suppressed, both the vertical movement (uplift) of isopycnal surfaces and upward motion along sloping isopycnals supply nutrients to the euphotic layer, but the relative importance of each of these mechanisms is unknown. We present a method for decomposing vertical velocity w into two components in a Lagrangian frame: vertical velocity along sloping isopycnal surfaces and the adiabatic vertical velocity of isopycnal surfaces . We show that , where is the isopycnal slope and is the geometric aspect ratio of the flow, and that accounts for 10%–25% of the total vertical velocity w for isopycnal slopes representative of the midlatitude pycnocline. We perform the decomposition of w in a process study model of a midlatitude eddying flow field generated with a range of isopycnal slopes. A spectral decomposition of the velocity components shows that while is the largest contributor to vertical velocity, is of comparable magnitude at horizontal scales less than about 10 km, that is, at submesoscales. Increasing the horizontal grid resolution of models is known to increase vertical velocity; this increase is disproportionately due to better resolution of , as is shown here by comparing 1- and 4-km resolution model runs. Along-isopycnal vertical transport can be an important contributor to the vertical flux of tracers, including oxygen, nutrients, and chlorophyll, although we find weak covariance between vertical velocity and nutrient anomaly in our model.MAF was supported by a National Defense Science and Engineering Graduate Fellowship and AM by NSF OCE-I434788. The authors thank Glenn Flierl and Ruth Curry for helpful conversations, and three anonymous reviewers for comments that improved the manuscript.2020-06-1

Submesoscale-selective compensation of fronts in a salinity-stratified ocean

© The Author(s), 2018. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Science Advances 4 (2018): e1701504, doi:10.1126/sciadv.1701504.Salinity, rather than temperature, is the leading influence on density in some regions of the world’s upper oceans. In the Bay of Bengal, heavy monsoonal rains and runoff generate strong salinity gradients that define density fronts and stratification in the upper ~50 m. Ship-based observations made in winter reveal that fronts exist over a wide range of length scales, but at O(1)-km scales, horizontal salinity gradients are compensated by temperature to alleviate about half the cross-front density gradient. Using a process study ocean model, we show that scale-selective compensation occurs because of surface cooling. Submesoscale instabilities cause density fronts to slump, enhancing stratification along-front. Specifically for salinity fronts, the surface mixed layer (SML) shoals on the less saline side, correlating sea surface salinity (SSS) with SML depth at O(1)-km scales. When losing heat to the atmosphere, the shallower and less saline SML experiences a larger drop in temperature compared to the adjacent deeper SML on the salty side of the front, thus correlating sea surface temperature (SST) with SSS at the submesoscale. This compensation of submesoscale fronts can diminish their strength and thwart the forward cascade of energy to smaller scales. During winter, salinity fronts that are dynamically submesoscale experience larger temperature drops, appearing in satellite-derived SST as cold filaments. In freshwater-influenced regions, cold filaments can mark surface-trapped layers insulated from deeper nutrient-rich waters, unlike in other regions, where they indicate upwelling of nutrient-rich water and enhanced surface biological productivity.This work was carried out under the Office of Naval Research’s ASIRI (grants N000141612470 and N000141310451) in collaboration with the Indian Ministry of Earth Science’s OMM initiative supported by the Monsoon Missio

Sinking flux of particulate organic matter in the oceans: Sensitivity to particle characteristics

© The Author(s), 2020. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Omand, M. M., Govindarajan, R., He, J., & Mahadevan, A. Sinking flux of particulate organic matter in the oceans: Sensitivity to particle characteristics. Scientific Reports, 10(1), (2020): 5582, doi:10.1038/s41598-020-60424-5.The sinking of organic particles produced in the upper sunlit layers of the ocean forms an important limb of the oceanic biological pump, which impacts the sequestration of carbon and resupply of nutrients in the mesopelagic ocean. Particles raining out from the upper ocean undergo remineralization by bacteria colonized on their surface and interior, leading to an attenuation in the sinking flux of organic matter with depth. Here, we formulate a mechanistic model for the depth-dependent, sinking, particulate mass flux constituted by a range of sinking, remineralizing particles. Like previous studies, we find that the model does not achieve the characteristic ‘Martin curve’ flux profile with a single type of particle, but instead requires a distribution of particle sizes and/or properties. We consider various functional forms of remineralization appropriate for solid/compact particles, and aggregates with an anoxic or oxic interior. We explore the sensitivity of the shape of the flux vs. depth profile to the choice of remineralization function, relative particle density, particle size distribution, and water column density stratification, and find that neither a power-law nor exponential function provides a definitively superior fit to the modeled profiles. The profiles are also sensitive to the time history of the particle source. Varying surface particle size distribution (via the slope of the particle number spectrum) over 3 days to represent a transient phytoplankton bloom results in transient subsurface maxima or pulses in the sinking mass flux. This work contributes to a growing body of mechanistic export flux models that offer scope to incorporate underlying dynamical and biological processes into global carbon cycle models.We thank NSF (OCE 1260080), NASA (NNX16AR48G), and the Ministry of Earth Sciences, Government of India (Monsoon Mission Project on the Bay of Bengal) for support. This work was largely done in 2012 while MMO was a postdoctoral associate at WHOI, during a visit by RG supported by The Mary Sears visiting scholar program to the Woods Hole Oceanographic Institution. Thanks also to Benjamin Hodges for many thoughtful contributions

Monami as an oscillatory hydrodynamic instability in a submerged sea grass bed

The onset of monami ~-- the synchronous waving of sea grass beds driven by a steady flow -- is modeled as a linear instability of the flow. Unlike previous works, our model considers the drag exerted by the grass in establishing the steady flow profile, and in damping out perturbations to it. We find two distinct modes of instability, which we label Mode 1 and Mode 2. Mode 1 is closely related to Kelvin-Helmholtz instability modified by vegetation drag, whereas Mode 2 is unrelated to Kelvin-Helmholtz and arises from an interaction between the flow in the vegetated and unvegetated layers. The vegetation damping, according to our model, leads to a finite threshold flow for both these modes. Experimental observations for the onset and frequency of waving compare well with model predictions for the instability onset criteria and the imaginary part of the complex growth rate respectively, but experiments lie in a parameter regime where the two modes can not be distinguished. % The inclusion of vegetation drag differentiates our mechanism from the previous linear stability analyses of monami.Comment: 4 figures, 13 page

Large-scale alignment of oceanic nitrate and density

Author Posting. © American Geophysical Union, 2013. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 118 (2013): 5322–5332, doi:10.1002/jgrc.20379.By analyzing global data, we find that over large scales, surfaces of constant nitrate are often better aligned with isopycnals than with isobars, particularly below the euphotic zone. This is unexplained by the movement of isopycnal surfaces in response to eddies and internal waves, and is perhaps surprising given that the biological processes that alter nitrate distributions are largely depth dependent. We provide a theoretical framework for understanding the orientation of isonitrate surfaces in relation to isopycnals. In our model, the nitrate distribution results from the balance between depth-dependent biological processes (nitrate uptake and remineralization), and the along-isopycnal homogenization of properties by eddy fluxes (parameterized by eddy diffusivity). Where the along-isopycnal eddy diffusivity is relatively large, nitrate surfaces are better aligned with isopycnals than isobars. We test our theory by estimating the strength of the eddy diffusivity and biological export production from global satellite data sets and comparing their contributions. Indeed, we find that below the euphotic zone, the mean isonitrate surfaces are oriented along isopycnals where the isopycnal eddy diffusivity is large, and deviate where the biological export of organic matter is relatively strong. Comparison of nitrate data from profiling floats in different regions corroborates the hypothesis by showing variations in the nitrate-density relationship from one part of the ocean to another.We acknowledge the support of the National Science Foundation (Grant OCE-0928617) and NASA (Grant NNX- 08AL80G).2014-04-1

Enhancement in vertical fluxes at a front by mesoscale-submesoscale coupling

Author Posting. © American Geophysical Union, 2014. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 119 (2014): 8495–8511, doi:10.1002/2014JC010211.Oceanic frontal instabilities are of importance for the vertical exchange of properties in the ocean. Submesoscale, O(1) Rossby number, dynamics are particularly relevant for inducing the vertical (and lateral) flux of buoyancy and tracers in the mixed layer, but how these couple with the stratified pycnocline is less clear. Observations show surface fronts often persist beneath the mixed layer. Here we use idealized, three-dimensional model simulations to show how surface fronts that extend deeper into the pycnocline invoke enhanced vertical fluxes through the coupling of submesoscale and mesoscale instabilities. We contrast simulations in which the front is restricted to the mixed layer with those in which it extends deeper. For the deeper fronts, we examine the effect of density stratification on the vertical coupling. Our results show deep fronts can dynamically couple the mixed layer and pycnocline on time scales that increase with the peak stratification beneath the mixed layer. Eddies in the interior generate skew fluxes of buoyancy and tracer oriented along isopycnals, thus providing an adiabatic pathway for the interior to interact with the mixed layer at fronts. The vertical enhancement of tracer fluxes through the mesoscale-submesoscale coupling described here is thus relevant to the vertical supply of nutrients for phytoplankton in the ocean. A further implication for wind-forced fronts is that the vertical structure of the stream function characterizing the exchange between the interior and the mixed layer exhibits significant qualitative differences compared to a linear combination of existing parameterizations of submesoscale eddies in the mixed layer and mesoscale eddies in the interior. The discrepancies are most severe within the mixed layer suggesting a potential role for Ekman-layer dynamics absent in existing submesoscale parameterizations.S.R. and A.T. acknowledge financial support from the National Science Foundation (NSF OCE-0928138) and the Office of Naval Research (ONR N00014-09-1-0196, ONR N00014-12-1-0101). A.M. acknowledges funding from the National Science Foundation (NSF OCE-0928617) and the Office of Naval Research (ONR N00014-12-1-0101).2015-06-1

Spontaneous generation of near-inertial waves by the Kuroshio Front

Author Posting. © American Meteorological Society, 2015. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 45 (2015): 2381–2406, doi:10.1175/JPO-D-14-0086.1.While near-inertial waves are known to be generated by atmospheric storms, recent observations in the Kuroshio Front find intense near-inertial internal-wave shear along sloping isopycnals, even during calm weather. Recent literature suggests that spontaneous generation of near-inertial waves by frontal instabilities could represent a major sink for the subinertial quasigeostrophic circulation. An unforced three-dimensional 1-km-resolution model, initialized with the observed cross-Kuroshio structure, is used to explore this mechanism. After several weeks, the model exhibits growth of 10–100-km-scale frontal meanders, accompanied by O(10) mW m−2 spontaneous generation of near-inertial waves associated with readjustment of submesoscale fronts forced out of balance by mesoscale confluent flows. These waves have properties resembling those in the observations. However, they are reabsorbed into the model Kuroshio Front with no more than 15% dissipating or radiating away. Thus, spontaneous generation of near-inertial waves represents a redistribution of quasigeostrophic energy rather than a significant sink.“The Study of Kuroshio Ecosystem Dynamics for Sustainable Fisheries (SKED)” supported by MEXT, MIT-Hayashi Seed Fund, ONR (Awards N000140910196 and N000141210101), NSF (Award OCE 0928617, 0928138) for support.2016-03-0