Stukel, M. R., & Kelly, T. B. (2019). The carbon: (234) Thorium ratios of sinking particles in the California current ecosystem 2: Examination of a thorium sorption, desorption, and particle transport model. Marine Chemistry , 212 , 1–15.
Abstract: Thorium-234 (234Th) is a powerful tracer of particle dynamics and the biological pump in the surface ocean; however, variability in carbon: thorium ratios of sinking particles adds substantial uncertainty to estimates of organic carbon export. We coupled a mechanistic thorium sorption and desorption model to a one-dimensional particle sinking model that uses realistic particle settling velocity spectra. The model generates estimates of 238U234Th disequilibrium, particulate organic carbon concentration, and the C:234Th ratio of sinking particles, which are then compared to in situ measurements from quasi-Lagrangian studies conducted on six cruises in the California Current Ecosystem. Broad patterns observed in in situ measurements, including decreasing C:234Th ratios with depth and a strong correlation between sinking C:234Th and the ratio of vertically-integrated particulate organic carbon (POC) to vertically-integrated total water column 234Th, were accurately recovered by models assuming either a power law distribution of sinking speeds or a double log normal distribution of sinking speeds. Simulations suggested that the observed decrease in C:234Th with depth may be driven by preferential remineralization of carbon by particle-attached microbes. However, an alternate model structure featuring complete consumption and/or disaggregation of particles by mesozooplankton (e.g. no preferential remineralization of carbon) was also able to simulate decreasing C:234Th with depth (although the decrease was weaker), driven by 234Th adsorption onto slowly sinking particles. Model results also suggest that during bloom decays C:234Th ratios of sinking particles should be higher than expected (based on contemporaneous water column POC), because high settling velocities minimize carbon remineralization during sinking.
Stukel, M. R., Biard, T., Krause, J. W., & Ohman, M. D. (2018). Large Phaeodaria in the twilight zone: Their role in the carbon cycle. Association for the Sciences of Limnology and Oceanography , .
Abstract: Advances in in situ imaging allow enumeration of abundant populations of large Rhizarians that compose a substantial proportion of total mesozooplankton biovolume. Using a quasi-Lagrangian sampling scheme, we quantified the abundance, vertical distributions, and sinking‐related mortality of Aulosphaeridae, an abundant family of Phaeodaria in the California Current Ecosystem. Inter‐cruise variability was high, with average concentrations at the depth of maximum abundance ranging from < 10 to > 300 cells m−3, with seasonal and interannual variability associated with temperature‐preferences and regional shoaling of the 10°C isotherm. Vertical profiles showed that these organisms were consistently most abundant at 100�150 m depth. Average turnover times with respect to sinking were 4.7�10.9 d, equating to minimum in situ population growth rates of ~ 0.1�0.2 d−1. Using simultaneous measurements of sinking organic carbon, we find that these organisms could only meet their carbon demand if their carbon : volume ratio were ~ 1 μg C mm−3. This value is substantially lower than previously used in global estimates of rhizarian biomass, but is reasonable for organisms that use large siliceous tests to inflate their cross‐sectional area without a concomitant increase in biomass. We found that Aulosphaeridae alone can intercept > 20% of sinking particles produced in the euphotic zone before these particles reach a depth of 300 m. Our results suggest that the local (and likely global) carbon biomass of Aulosphaeridae, and probably the large Rhizaria overall, needs to be revised downwards, but that these organisms nevertheless play a major role in carbon flux attenuation in the twilight zone.
Stukel, M. R., Decima, M., & Kelly, T. B. (2018). A new approach for incorporating 15N isotopic data into linear inverse ecosystem models with Markov Chain Monte Carlo sampling. PLoS One , 13 (6), e0199123.
Abstract: Oceanographic field programs often use delta15N biogeochemical measurements and in situ rate measurements to investigate nitrogen cycling and planktonic ecosystem structure. However, integrative modeling approaches capable of synthesizing these distinct measurement types are lacking. We develop a novel approach for incorporating delta15N isotopic data into existing Markov Chain Monte Carlo (MCMC) random walk methods for solving linear inverse ecosystem models. We test the ability of this approach to recover food web indices (nitrate uptake, nitrogen fixation, zooplankton trophic level, and secondary production) derived from forward models simulating the planktonic ecosystems of the California Current and Amazon River Plume. We show that the MCMC with delta15N approach typically does a better job of recovering ecosystem structure than the standard MCMC or L2 minimum norm (L2MN) approaches, and also outperforms an L2MN with delta15N approach. Furthermore, we find that the MCMC with delta15N approach is robust to the removal of input equations and hence is well suited to typical pelagic ecosystem studies for which the system is usually vastly under-constrained. Our approach is easily extendable for use with delta13C isotopic measurements or variable carbon:nitrogen stoichiometry.
Stukel, M. R., Décima, M., Landry, M. R., & Selph, K. E. (2018). Nitrogen and isotope flows through the Costa Rica Dome upwelling ecosystem: The crucial mesozooplankton role in export flux. Global Biogeochemical Cycles , 32 (12), 1815–1832.
Abstract: The Costa Rica Dome (CRD) is an open-ocean upwelling ecosystem, with high biomasses of picophytoplankton (especially Synechococcus), mesozooplankton, and higher trophic levels. To elucidate the food web pathways supporting the trophic structure and carbon export in this unique ecosystem, we used Markov Chain Monte Carlo techniques to assimilate data from four independent realizations of δ15N and planktonic rate measurements from the CRD into steady state, multicompartment ecosystem box models (linear inverse models). Model results present well-constrained snapshots of ecosystem nitrogen and stable isotope fluxes. New production is supported by upwelled nitrate, not nitrogen fixation. Protistivory (rather than herbivory) was the most important feeding mode for mesozooplankton, which rely heavily on microzooplankton prey. Mesozooplankton play a central role in vertical nitrogen export, primarily through active transport of nitrogen consumed in the surface layer and excreted at depth, which comprised an average 36-46% of total export. Detritus or aggregate feeding is also an important mode of resource acquisition by mesozooplankton and regeneration of nutrients within the euphotic zone. As a consequence, the ratio of passively sinking particle export to phytoplankton production is very low in the CRD. Comparisons to similar models constrained with data from the nearby equatorial Pacific demonstrate that the dominant role of vertical migrators to the biological pump is a unique feature of the CRD. However, both regions show efficient nitrogen transfer from mesozooplankton to higher trophic levels (as expected for regions with large fish, cetacean, and seabird populations) despite the dominance of protists as major grazers of phytoplankton.
Stukel, M. R., Ohman, M. D., Kelly, T. B., & Biard, T. (2019). The Roles of Suspension-Feeding and Flux-Feeding Zooplankton as Gatekeepers of Particle Flux Into the Mesopelagic Ocean in the Northeast Pacific. Front. Mar. Sci. , 6 .
Abstract: Zooplankton are important consumers of sinking particles in the ocean's twilight zone. However, the impact of different taxa depends on their feeding mode. In contrast to typical suspension-feeding zooplankton, flux-feeding taxa preferentially consume rapidly sinking particles that would otherwise penetrate into the deep ocean. To quantify the potential impact of two flux-feeding zooplankton taxa [Aulosphaeridae (Rhizaria), and Limacina helicina (euthecosome pteropod)] and the total suspension-feeding zooplankton community, we measured depth-stratified abundances of these organisms during six cruises in the California Current Ecosystem. Using allometric-scaling relationships, we computed the percentage of carbon flux intercepted by flux feeders and suspension feeders. These estimates were compared to direct measurements of carbon flux attenuation (CFA) made using drifting sediment traps and U-238-Th-234 disequilibrium. We found that CFA in the shallow twilight zone typically ranged from 500 to 1000 m mol organic C flux remineralized per 10-m vertical depth bin. This equated to approximately 6-10% of carbon flux remineralized/10 m. The two flux-feeding taxa considered in this study could account for a substantial proportion of this flux near the base of the euphotic zone. The mean flux attenuation attributable to Aulosphaeridae was 0.69%/10 m (median = 0.21%/10 m, interquartile range = 0.04-0.81%) at their depth of maximum abundance (similar to 100 m), which would equate to similar to 10% of total flux attenuation in this depth range. The maximum flux attenuation attributable to Aulosphaeridae reached 4.2%/10 m when these protists were most abundant. L. helicina, meanwhile, could intercept 0.45-1.6% of carbon flux/10 m, which was slightly greater (on average) than the Aulosphaeridae. In contrast, suspension-feeding zooplankton in the mesopelagic (including copepods, euphausiids, appendicularians, and ostracods) had combined clearance rates of 2-81 L m(-3) day(-1) (mean of 19.6 L m(-3) day(-1)). This implies a substantial impact on slowly sinking particles, but a negligible impact on the presumably rapidly sinking fecal pellets that comprised the majority of the material collected in sediment traps. Our results highlight the need for a greater research focus on the many taxa that potentially act as flux feeders in the oceanic twilight zone.
Sun, J., & Wu, Z. (2019). Isolating spatiotemporally local mixed Rossby-gravity waves using multi-dimensional ensemble empirical mode decomposition. Clim Dyn , (3-4), 1383–1405.
Abstract: Tropical waves have relatively large amplitudes in and near convective systems, attenuating as they propagate away from the area where they are generated due to the dissipative nature of the atmosphere. Traditionally, nonlocal analysis methods, such as those based on the Fourier transform, are applied to identify tropical waves. However, these methods have the potential to lead to the misidentification of local wavenumbers and spatial locations of local wave activities. To address this problem, we propose a new method for analyzing tropical waves, with particular focus placed on equatorial mixed Rossby-gravity (MRG) waves. The new tropical wave analysis method is based on the multi-dimensional ensemble empirical mode decomposition and a novel spectral representation based on spatiotemporally local wavenumber, frequency, and amplitude of waves. We first apply this new method to synthetic data to demonstrate the advantages of the method in revealing characteristics of MRG waves. We further apply the method to reanalysis data (1) to identify and isolate the spatiotemporally heterogeneous MRG waves event by event, and (2) to quantify the spatial inhomogeneity of these waves in a wavenumber-frequency-energy diagram. In this way, we reveal the climatology of spatiotemporal inhomogeneity of MRG waves and summarize it in wavenumber-frequency domain: The Indian Ocean is dominated by MRG waves in the period range of 8–12 days; the western Pacific Ocean consists of almost equal energy distribution of MRG waves in the period ranges of 3–6 and 8–12 days, respectively; and the eastern tropical Pacific Ocean and the tropical Atlantic Ocean are dominated by MRG waves in the period range of 3–6 days. The zonal wavenumbers mostly fall within the band of 4–15, with Indian Ocean has larger portion of higher wavenumber (smaller wavelength components) MRG waves.
Timko, P. G., Arbic, B. K., Hyder, P., Richman, J. G., Zamudio, L., O'Dea, E., et al. (2019). Assessment of shelf sea tides and tidal mixing fronts in a global ocean model. Ocean Modelling , 136 , 66–84.
Abstract: Tidal mixing fronts, which represent boundaries between stratified and tidally mixed waters, are locations of enhanced biological activity. They occur in summer shelf seas when, in the presence of strong tidal currents, mixing due to bottom friction balances buoyancy production due to seasonal heat flux. In this paper we examine the occurrence and fidelity of tidal mixing fronts in shelf seas generated within a global 3-dimensional simulation of the HYbrid Coordinate Ocean Model (HYCOM) that is simultaneously forced by atmospheric fields and the astronomical tidal potential. We perform a first order assessment of shelf sea tides in global HYCOM through comparison of sea surface temperature, sea surface tidal elevations, and tidal currents with observations. HYCOM was tuned to minimize errors in M2 sea surface heights in deep water. Over the global coastal and shelf seas (depths <200 m) the area-weighted root mean square error of the M2 sea surface amplitude in HYCOM represents 35% of the 50 cm root mean squared M2 sea surface amplitude when compared to satellite constrained models TPXO8 and FES2014. HYCOM and the altimeter constrained tidal models TPXO8 and FES2014 exhibit similar skill in reproducing barotropic tidal currents estimated from in-situ current meter observations. Through comparison of a global HYCOM simulation with tidal forcing to a global HYCOM simulation with no tides, and also to previous regional studies of tidal mixing fronts in shelf seas, we demonstrate that HYCOM with embedded tides exhibits quite high skill in reproducing known tidal mixing fronts in shelf seas. Our results indicate that the amount of variability in the location of the tidal mixing fronts in HYCOM, estimated using the Simpson-Hunter parameter, is consistent with previous studies when the differences in the net downward heat flux, on a global scale, are taken into account. We also provide evidence of tidal mixing fronts on the North West Australian Shelf for which we have been unable to find references in the existing scientific literature.
Venugopal, T., Ali, M. M., Bourassa, M. A., Zheng, Y., Goni, G. J., Foltz, G. R., et al. (2018). Statistical Evidence for the Role of Southwestern Indian Ocean Heat Content in the Indian Summer Monsoon Rainfall. Sci Rep , 8 (1), 12092.
Abstract: This study examines the benefit of using Ocean Mean Temperature (OMT) to aid in the prediction of the sign of Indian Summer Monsoon Rainfall (ISMR) anomalies. This is a statistical examination, rather than a process study. The thermal energy needed for maintaining and intensifying hurricanes and monsoons comes from the upper ocean, not just from the thin layer represented by sea surface temperature (SST) alone. Here, we show that the southwestern Indian OMT down to the depth of the 26 degrees C isotherm during January-March is a better qualitative predictor of the ISMR than SST. The success rate in predicting above- or below-average ISMR is 80% for OMT compared to 60% for SST. Other January-March mean climate indices (e.g., NINO3.4, Indian Ocean Dipole Mode Index, El Nino Southern Oscillation Modoki Index) have less predictability (52%, 48%, and 56%, respectively) than OMT percentage deviation (PD) (80%). Thus, OMT PD in the southwestern Indian Ocean provides a better qualitative prediction of ISMR by the end of March and indicates whether the ISMR will be above or below the climatological mean value.
Villas Bôas, A. B., Ardhuin, F., Ayet, A., Bourassa, M. A., Brandt, P., Chapron, B., et al. (2019). Integrated Observations of Global Surface Winds, Currents, and Waves: Requirements and Challenges for the Next Decade. Front. Mar. Sci. , 6 .
Abstract: Ocean surface winds, currents, and waves play a crucial role in exchanges of momentum, energy, heat, freshwater, gases, and other tracers between the ocean, atmosphere, and ice. Despite surface waves being strongly coupled to the upper ocean circulation and the overlying atmosphere, efforts to improve ocean, atmospheric, and wave observations and models have evolved somewhat independently. From an observational point of view, community efforts to bridge this gap have led to proposals for satellite Doppler oceanography mission concepts, which could provide unprecedented measurements of absolute surface velocity and directional wave spectrum at global scales. This paper reviews the present state of observations of surface winds, currents, and waves, and it outlines observational gaps that limit our current understanding of coupled processes that happen at the air-sea-ice interface. A significant challenge for the coming decade of wind, current, and wave observations will come in combining and interpreting measurements from (a) wave-buoys and high-frequency radars in coastal regions, (b) surface drifters and wave-enabled drifters in the open-ocean, marginal ice zones, and wave-current interaction �hot-spots,� and (c) simultaneous measurements of absolute surface currents, ocean surface wind vector, and directional wave spectrum from Doppler satellite sensors.
Vinayachandran, P. N., Davidson, F., & Chassignet, E. P. (2020). Towards joint assessments, modern capabilities and new links for ocean prediction systems. Bull. Amer. Meteor. Soc. , 101 (4).
Abstract: Approximately 260 individuals from forecasting centers, research laboratories, academia, and industry representing 40 countries met to discuss recent developments in operational oceanography and brainstorm about the future directions of ocean prediction services.
