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.
Subrahmanyam, B., Murty, V. S. N., & O'Brien, J. J. (2005). New Sea Surface Salinity Product in the Tropical Indian Ocean . CAS/JSC Working Group on Numerical Experimentation.
Subrahmanyam, B., J.J. O'Brien, V. Manghanani, L. Xie, and J.M. Morrison. (2003). Estimation of Heat Transports in the Indian Ocean uisng Altimetry and MICOMi. In WMO WORLD Climate Research Programme, CAS/JSC Working Group on Numerical Experimentation, Research Activities in Atmospheric and Oceanic Modelling .
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.
Tartaglione, C. A., Smith, S. R., & O'Brien, J. J. (2003). ENSO Impact on Hurricane Landfall Probabilities for the Caribbean. J. Climate , 16 (17), 2925–2931.
Tartaglione, C. S., Hanley, D. E., O'Brien, J. J., & Smith, S. R. (2002). Regional Effects of ENSO on U.S Hurricane Landfalls . COAPS Technical Report 02-5. Tallahassee, FL: Center for Ocean-Atmospheric Prediction Studies, Florida State University.
Taylor, J. P. (2006). Comparison of ECMWF and Quikscat-Derived Surface Pressure Gradients . Master's thesis, Florida State University, Tallahassee, FL.
Abstract: A technique based solely on QuikSCAT data is developed for determining suspect differences between QSCAT and ECMWF pressure gradients. Pressure fields are computed from scatterometer winds using a variational method that applies a gradient wind conversion. Kinematic analysis of the satellite wind field is performed in order to determine which parameters are physically related to the suspect pressure gradients. It is discovered that the likelihood of these suspect occurrences has the greatest dependence on relative vorticity, total deformation, and the curvature Rossby number. A broad range of these values is tested and a single assessment criterion is derived based upon the value of several skill scores. Overall, the assessment criterion is able to correctly identify the majority of suspect pressure gradients; yet considerable over-flagging does occur in many instances. However, the over-flagging is not random: the false alarms are tightly clustered around the suspect areas, resulting in flagged regions that are too large. Identification of the location of suspect areas in pressure products should be useful to forecasters.
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.
