Oceanographic Structure and Light Levels Drive Patterns of Sound Scattering Layers in a Low-Latitude Oceanic System

Boswell, Kevin M.; D’Elia, Marta; Johnston, Matthew W.; Mohan, John A.; Warren, Joseph D.; Wells, R. J. David; Sutton, Tracey T.

doi:10.3389/fmars.2020.00051

ORIGINAL RESEARCH article

Front. Mar. Sci., 19 February 2020

Sec. Deep-Sea Environments and Ecology

Volume 7 - 2020 | https://doi.org/10.3389/fmars.2020.00051

This article is part of the Research Topic Deep Pelagic Ecosystem Dynamics in a Highly Impacted Water Column: The Gulf of Mexico After Deepwater Horizon View all 14 articles

Oceanographic Structure and Light Levels Drive Patterns of Sound Scattering Layers in a Low-Latitude Oceanic System

$\r\nKevin M. Boswell*$ Kevin M. Boswell^1*

R. J. David Wells^3,5

Tracey T. Sutton²

¹Department of Biological Sciences, Marine Science Program, Florida International University, North Miami, FL, United States
²Halmos College of Natural Sciences and Oceanography, Nova Southeastern University, Dania Beach, FL, United States
³Department of Marine Biology, Texas A&M University at Galveston, Galveston, TX, United States
⁴School of Marine and Atmospheric Sciences, Stony Brook University, Southampton, NY, United States
⁵Department of Wildlife and Fisheries Sciences, Texas A&M University, College Station, TX, United States

Several factors have been reported to structure the spatial and temporal patterns of sound scattering layers, including temperature, oxygen, salinity, light, and physical oceanographic conditions. In this study, we examined the spatiotemporal variability of acoustically detected sound scattering layers in the northern Gulf of Mexico to investigate the drivers of this variability, including mesoscale oceanographic features [e.g., Loop Current-origin water (LCOW), frontal boundaries, and Gulf Common Water]. Results indicate correlations in the vertical position and acoustic backscatter intensity of sound scattering layers with oceanographic conditions and light intensity. LCOW regions displayed consistent decreases, by a factor of two and four, in acoustic backscatter intensity in the upper 200 m relative to frontal boundaries and Gulf Common Water, respectively. Sound scattering layers had greater backscatter intensity at night in comparison to daytime (25x for frontal boundaries, 17x for LCOW, and 12x for Gulf Common Water). The importance of biotic (primary productivity) and abiotic (sea surface temperature, salinity) factors varied across oceanographic conditions and depth intervals, suggesting that the patterns in distribution and behavior of mesopelagic assemblages in low-latitude, oligotrophic ecosystems can be highly dynamic.

Introduction

The oceanic biome is approximately 71% of the planet’s area and much more of the planet’s living space by volume, yet it remains vastly understudied (Childress, 1983; Webb et al., 2010). Perhaps the most conspicuous features of this biome are the persistent and ubiquitous sound scattering layers (Marshall, 1954; Barham, 1966; Gjosaeter and Kawaguchi, 1980; Irigoien et al., 2014; Cade and Benoit-Bird, 2015; Davison et al., 2015) formed by zooplankton and micronekton (Kloser et al., 2002; Irigoien et al., 2014; Béhagle et al., 2017). These organisms are responsible for the Earth’s largest animal migration, a process known as diel vertical migration (DVM) (Marshall, 1954; Pearre, 2003; Brierley, 2014; Aksnes et al., 2017; Behrenfeld et al., 2019). Recently, the fish component of the global mesopelagic micronekton community (crustaceans, cephalopods, and fishes, ∼2–10 cm in length) was estimated to exceed 5 billion tons (Irigoien et al., 2014; Klevjer et al., 2016; Aksnes et al., 2017).

The migrating layers serve as important trophic pathways linking meso- and bathypelagic habitats with the epipelagic through active vertical movement of animals. In general, micronekton actively swim toward the surface at dusk, seeking foraging opportunities (Merrett and Roe, 1974; Brodeur et al., 2005; Bianchi et al., 2013; Sutton, 2013; Sutton et al., 2020), and descend at dawn into the deep ocean. An important consequence of DVM is that it facilitates trophic interactions and biogeochemical exchange, vertically integrating the world’s oceans (Sutton and Hopkins, 1996a; Hidaka et al., 2001; Davison et al., 2013, Davison et al., 2015; Schukat et al., 2013; Hudson et al., 2014; Trueman et al., 2014; Ariza et al., 2016; Sutton et al., 2020) through extensive, coordinated animal movement. In addition to significant contributions to the biological pump, mesopelagic communities also serve an important role in oceanic food webs by facilitating linkages among secondary producers (zooplankton) and higher-level consumers, including oceanic apex predators (Robertson and Chivers, 1997; Potier et al., 2007; Spear et al., 2007; Benoit-Bird et al., 2017). In spite of the fundamental ecological importance for open-ocean functioning, and increasing interest in commercial exploitation, the mesopelagic community remains one of the least-studied components of oceanic systems (Handegard et al., 2013; Irigoien et al., 2014; Davison et al., 2015).

The spatial and temporal variability observed in the sound scattering layers are known to fluctuate horizontally and vertically across abiotic and biotic gradients, primarily temperature (Kumar et al., 2005; Brierley, 2014; Béhagle et al., 2017; Proud et al., 2017), oxygen content (Devol, 1981; Bertrand et al., 2010; Brierley, 2014; Béhagle et al., 2017), salinity (Forward, 1976; Wang et al., 2014), and light intensity (Frank and Widder, 1997; Aksnes et al., 2009; Lebourges-Dhaussy et al., 2014; Last et al., 2016; Aksnes et al., 2017; Kaartvedt et al., 2017). The importance of these factors in structuring sound scattering layers can vary and is dependent on the location, community dynamics, and physical setup of the oceanic system. For example, the processes that structure high-latitude systems may vary in scale relative to mid-latitude or tropical systems (Godø et al., 2012; Peña et al., 2014; Røstad et al., 2016; Aksnes et al., 2017).

Mesoscale oceanographic features (e.g., eddies, frontal boundaries) have also been identified as important in mediating the dynamics of sound scattering layers (Owen, 1981; Sabarros et al., 2009; Godø et al., 2012; Scales et al., 2014; Ternon et al., 2014; Gaube et al., 2018). These features operate across multiple spatial scales (10–100s km) and produce areas of physical and biological heterogeneity and are thought to play an important role in mediating the transport and accumulation of biological material (i.e., larvae, eggs, plankton) as well as nutrients and heat (Sabarros et al., 2009; Chelton et al., 2011). Given that these features exist across a wide range of oceanic geographies, the ubiquitous vertically migrating sound scattering layers are also likely to be influenced by these mesoscale features (Fennell and Rose, 2015).

Within the Gulf of Mexico (GoM), the dominant mesoscale features are eddies and frontal boundaries associated with the Loop Current, interspersed among larger regions of Gulf Common Water (Johnston et al., 2019). The Loop Current is formed by warm, highly saline Caribbean water entering the GoM through the Yucatan Channel. The Loop Current’s position within the GoM varies and is dependent upon its retracted or extended state. When retracted, the Loop Current flows directly east from the Yucatan Channel, bypassing the GoM proper, flanking the Florida Keys, eventually forming the Gulf Stream in the North Atlantic. When extended, the Loop Current protrudes into the far north and eastern GoM as far as 28° north latitude. The Loop Current drives environmental heterogeneity in the upper 1000 m in the pelagic GoM (Cardona and Bracco, 2016) through the shedding of energetic eddies, both cyclonic and anticyclonic. Loop Current eddies are large (10–100 kms in diameter) and persistent (average lifespan of 8–9 months; Hall and Leben, 2016) anticyclonic (downwelling) features. They are characterized by elevated mean sea surface height anomalies, clockwise rotation, elevated temperatures extending to c. 1000 m water depth (Biggs, 1992; Vukovich, 2007; Herring, 2010), and low surface chlorophyll a concentrations. Cyclonic eddies can be formed as well, although usually on the periphery of the large anticyclones, especially when large eddies are first sloughed off the Loop Current. Cyclones are typically much more ephemeral than anticyclones. Loop Current eddies typically form in the eastern GoM and ebb westward, eventually mixing with resident GoM water to form Gulf common water. Loop Current eddies and anticyclonic regions in the GoM can be distinguished from Gulf Common Water by the presence of the Subtropical Underwater water mass, which originates in the Caribbean (Rivas et al., 2005). The boundaries between these two water types are gradients, herein termed frontal boundaries (‘mixed water’ of Johnston et al., 2019), exhibit intermediate characteristics. These boundaries are known as important regions that concentrate or attract prey for pelagic organisms and may affect faunal distributions from surface waters to the benthos (Richards et al., 1993).

Previous examination of micronekton through acoustic-based surveys has indicated that mesoscale features may serve to structure mesopelagic organism distribution in oceanic systems (Drazen et al., 2011; Godø et al., 2012). The intent of this study was to examine the variability in acoustic backscatter associated with mesopelagic sound scattering layers among major oceanographic features in the GoM, a dynamic, oligotrophic, oceanic bioregion. The GoM represents an excellent model system for a study such as this, as high-resolution, taxon-specific vertical distribution data exist for the numerically dominant fishes (Hopkins and Lancraft, 1984; Gartner et al., 1987; Sutton and Hopkins, 1996a; Hopkins et al., 1996; Sutton et al., 2017; Milligan and Sutton, 2020), macrocrustaceans (Heffernan and Hopkins, 1981; Flock and Hopkins, 1992; Kinsey and Hopkins, 1994; Burdett et al., 2017; Frank et al., 2020), and cephalopods (Passarella and Hopkins, 1991; Judkins et al., 2016; Judkins and Vecchione, 2020). These data, developed during multi-decadal research programs (see references in Hopkins et al., 1996; Sutton et al., 2020), make the GoM one of the best-known deep-pelagic ecosystems in the World Ocean with respect to micronekton/nekton faunal composition and vertical distribution.

In this study, we sought to directly integrate existing biotic data, vessel-based acoustic surveys, remotely sensed oceanographic data, and predictive hydrographic ocean modeling to characterize the dominant mesoscale patterns of sound scattering layer distribution and intensity. Specifically, we examined how mesopelagic sound scattering layers respond to gradients in oceanographic conditions, light intensity, primary productivity, and temperature and salinity in order to better understand the drivers of pelagic ecosystem in a highly speciose, low-latitude pelagic ecosystem.

Materials and Methods

Study Area

Four acoustic surveys were conducted in the northern GoM aboard the R/V Point Sur (Figure 1) during the boreal late spring (DP01: 1–7 May 2015, Boswell, 2016; DP03: 1–15 May 2016, Boswell, 2017b) and summer (DP02: 9–23 August 2015, Boswell, 2017a; DP04: 6–21 August 2016, Boswell, 2017c). Sampling sites were an offshore extension of the standard Southeast Area Monitoring and Assessment Program (SEAMAP) plankton-sampling grid, which extends from the Texas shelf to the West Florida Shelf. Grid cells that comprise the survey design are 55.6 × 55.6 km, with sampling stations located at the mid-point of each grid cell. Cruise tracks were designed during the cruises to sample multiple oceanographic features. Paired multi-frequency acoustic and biological catch (10 m² MOCNESS net) data were collected at each site. The acoustic methodologies are described below. Net-sampling methodologies and subsequent data are described in detail in Kupchik et al. (2018), Milligan et al. (2019), Sutton et al. (2020), and Milligan and Sutton (2020).

FIGURE 1

Figure 1. Survey locations (filled symbols) overlaid on the feature classification map illustrating the gradients among the three primary features examined (Loop Current-origin water = red; frontal boundary = green; Gulf Common Water = blue). Lower panels represent survey tracks during each cruise (DP01–DP04). Contour lines represent 500 m depth intervals.

Acoustic Data Collection and Processing

Acoustic data were collected during the day and nighttime periods (defined by local sunrise and sunset) when the transducer was deployed, allowing for continuous surveys (∼8 h) during each transect at each station. A multiple-frequency echosounder system (Simrad EK60/Simrad EK80) was used and operated transducers at 18, 38, 70, 120 kHz. The transducers were mounted in a faring and suspended 2.5 m below the water surface. Given the limitations of using a pole-mounted system, transects were conducted at an approximate vessel speed of 2 knots. Transducers were calibrated according to the standard sphere method (Demer et al., 2015). For this paper, we examined the acoustic backscatter from the sound scattering layers using only the 38 kHz frequency due to: the widespread use of this frequency to study pelagic biomass (Davison et al., 2013; D’Elia et al., 2016; Aksnes et al., 2017; Kaartvedt et al., 2017), the complicating factor of resonance effects from gas-bearing organisms at 18 kHz, and low signal-to-noise ratios in the higher frequencies (70, 120 kHz) (Godø et al., 2009; Fennell and Rose, 2015; Davison et al., 2015). The pulse duration for 38 kHz echosounder was 4 ms with a power setting of 2000 W, and ping repetition rate of 0.2 pings s^–1. Sound speed profiles and absorption coefficient were computed from bin-averaged CTD data using the Ocean Toolbox (McDougall and Barker, 2011) in Matlab.

Raw acoustic backscatter data were imported and manually scrutinized in Echoview (v8, Myriax). Data from the transducer face to 15 m depth were excluded from the analysis to account for beam formation and to eliminate surface-associated interference (e.g., bubble sweep down). Data beyond 1000 m were not included in the analysis due to range dependent losses in attenuation and signal strength. Compromised data due to interference from other shipboard sonar systems (intermittent or spike noise), false bottom, and background noise were excluded from the analysis. False bottoms were manually excluded. To remove occurrences of spike noise, each sample was compared to the preceding and successive sample. If the single ping-to-ping difference was greater than 10 dB the sample was considered a spike candidate and replaced with the mean S_V of four neighboring samples (D’Elia et al., 2016). Background noise was identified and removed following a modified process described by De Robertis and Higginbottom (2007). A minimum signal-to-noise ratio of 15 dB was applied to data collected at 38 kHz. Samples that did not satisfy this threshold were considered indistinguishable from the background noise and flagged as ’no data.’ The measurements of Nautical Area Scattering Coefficient (NASC; m² nmi^–2) were derived from the echo integral in 500-m along-track x 5-m vertical bins with a −80 dB re 1 m² integration threshold (MacLennan et al., 2002). NASC is considered to be proportional to the abundance of biological scatterers and serves as a comparable index of organism biomass (Hazen et al., 2009; Zwolinski et al., 2010; Fennell and Rose, 2015). Integrated backscatter was further binned into three depth intervals: 15–200 m (epipelagic), 200–600 m (upper mesopelagic) and 600–1000 m (lower mesopelagic). The center of mass (m) was derived for each of the three depth intervals using the approach of Urmy et al. (2012) to describe depth of the statistical center of the backscatter within each depth interval.

Oceanographic Feature Identification Methods

Oceanographic feature classes were identified following Johnston et al. (2019); these include Loop Current-origin water (LCOW), Gulf Common Water, and frontal boundaries. These feature classes were derived from the GoM HYbrid Coordinate Ocean Model (HYCOM + NCODA Gulf of Mexico 1/25° Analysis, GoM l0.04/expt_32.5) (Chassignet et al., 2007), which is a three-dimensional, eddy-resolving circulation model that assimilates satellite- and in situ-derived measures to depict ocean conditions (e.g., sea surface height, zonal velocity, meridional velocity, temperature, and salinity) in near real time, from surface waters to the benthos. In the GoM, HYCOM data are available at 1/25° (c. 4 km²) horizontal resolution, in hourly intervals from 1993 to the present day¹ (Johnston et al., 2018, 2019). Velocity fronts were calculated as the difference between the minimum and maximum water speed within a 0.10 arc degree radius (∼11 km) of each location, derived from the HYCOM and measured in m s^–1.

LCOW is generally characterized by increased SSHA, increased water temperatures (extending down to ca. 1000 m), and a reduction in surface water chlorophyll concentrations. Based on the values used in Johnston et al. (2019), we derived an index that normalized the response of the LCOW and represented a derived quantity that utilized location-specific HYCOM output, given the equation:

LCOW index = SSHA {}_{i}- (SSHA {}_{GOM}+ 0.067) + T {}_{i}- 15.922 (1)

where SSHA_i represents the location-specific SSHA, SSHA_GOM is the daily mean SSHA in the GoM and T_i represents the location-specific temperature at 300 m. The output from Eq. 1 was scaled between 1.0 and 2.0, where 1.0 is the weakest condition and 2.0 represents the most intense condition measured for the entire GoM. The LCOW index. As such, the lower SSHA threshold for LCOW is SSHA_GOM + 0.067 m and the lower temperature threshold for the is 15.92°C, following Johnston et al. (2019).

Gulf Common Water is generally characterized by decreased SSHA and water column temperature and increased surface water chlorophyll concentrations when compared to LCOW. An index was computed for Gulf Common Water as the difference between the Johnston et al. (2019) temperature threshold at 300 m depth (13.46°C) and the location specific temperature at 300 m – i.e., the colder the water at 300 m, the greater the Gulf Common Water index value. The common water index was normalized to range from 0.0 to −1.0.

Frontal boundaries in the GoM are typically areas where significant mixing between LCOW and Common Water occur. To grade the strength of these boundaries, a frontal boundary index was calculated based on the difference between the SSHA_i and the SSHA_GOM and scaled to range from 0.0 to 1.0, with a value of 1.0 representing conditions nearest LCOW, and 0.0 being closest to Common Water (Johnston et al., 2019).

These oceanographic intensity indices were generated to span as a continuum following the classifications of Johnston et al. (2019) and then standardized and scaled to their respective ranges (i.e., LCOW: 1 to 2; frontal boundary: 0 to 1; Gulf Common Water: 0 to −1) based on the strength of the feature for the duration of the sampled period. Along-track positions for each acoustic survey were then used to extract and quantify the oceanographic conditions for each echo integration cell.

Hydrographic Properties of the Water Column

A calibrated SeaBird CTD (SBE 911+; SeaBird Electronics, Inc.) was used to characterize the water column properties. Data were collected during each night and day period as conditions allowed to characterize the diel structure of the water column. The raw instrument data were processed in the SeaBird processing software (v. 7.23), to compute 1-m bin-averaged estimates of salinity (PSU), temperature (°C), dissolved oxygen concentration (mg L^–1) and chl a (mg L^–1).

Approximating Surface Light Intensity and Primary Production

We examined the effect of relative light availability at 5 m depth by computing the instantaneous photosynthetically active radiation (IPAR; W m^–2) along each transect:

I P A R {}_{5 m}= I P A R {}_{s u r f}* 2.72^{(- K d (490) * z)} (2)

where the surface light intensity IPAR_surf was approximated by NOAA’s Geostationary Satellite Server (GOES); Kd(490) represents NASA’s Moderate Resolution Imaging Spectroradiometer (MODIS) derived diffuse attenuation coefficient at 490 nm, and z represents water depth (5 m). Observations at night have IPAR_surf values of 0. Hourly estimates of solar elevation were derived from NOAA’s Earth System Research Lab². The maximum estimates of solar elevation were 71.4° at local-noon (DP02) and minimum was −65.5° (DP04) at local-midnight (Figure 2). Sixty-day net primary production was compiled from the MODIS observations and estimated from the Vertically Generalized Production Model (Behrenfeld and Falkowski, 1997) made available from the Oregon State Ocean Productivity standard products³. Integrated net primary production estimates were extracted for each cruise (Supplementary Figure S1).

FIGURE 2

Figure 2. Estimated mean (black line) and 95% confidence intervals (gray lines) of light intensity at 5 m depth (IPAR₅_m, left panel) and solar elevation (right panel) across the four surveys conducted in the northern GoM during 2015 and 2016.

Net Collection

Micronekton were sampled with a 10 m² Multiple Opening and Closing Net Sampling System (MOCNESS) conducted synchronously with acoustic data. Briefly, the MOCNESS was used to sample discrete depth intervals from 0 to 1500 m water depth at each station. The MOCNESS was configured with 9 identical nets with 333 μm mesh (see Wiebe et al., 1985 for full system description). Samples were sorted, identified to lowest taxonomic level possible, enumerated, and weighed (either individually or in groups depending on size) onboard the vessel. Organisms were preserved in formalin for long-term storage and later analyses.

Data Analysis

Patterns in acoustic backscatter at 38 kHz (NASC, m² nmi^–2) and the center of mass (m) of sound scattering layers were examined by time of day (day and night), and across the depth intervals (D’Elia et al., 2016). We investigated these patterns relative to the three oceanographic feature classes (LCOW, frontal boundary, and Common Water), using a linear mixed effects model, implemented in R (R Core Team, 2013) with the library “nlme.” Since the variation in the residuals differed by day and night and across the three intervals of depth for both NASC and center of mass, a weighting option was added to the model using the varComb and varIdent structure to allow for different variances by time of day and depth domain. NASC values were log₁₀(x) transformed prior to analysis to meet the assumptions of normality. The responses in log-NASC and center of mass were examined relative to the interactions of time of day, depth interval and feature class. The cruise number was included as random effect to allow the magnitude of NASC and center of mass to vary by cruise. Tukey’s post hoc comparisons were used to identify significant differences with respect to log mean NASC and mean center of mass.

Generalized additive mixed models (GAMMs) were used to analyze the relationships of NASC and center of mass with the environmental and oceanographic drivers for each of the three depth intervals. HYCOM-derived sea surface temperature (°C), HYCOM-derived surface salinity (PSU), CTD-derived maximum chlorophyll concentration, and an index representing the gradient of the oceanographic feature classes were included as main effects. We also examined the interactions of surface (5 m) light intensity and the feature class indices as differences in light intensity at depth would be mediated by factors at the surface and likely reflect oceanographic differences.

Generalized additive mixed models were implemented with the “mgcv” library in R using a Gaussian distribution with an identity link function. Model selection was conducted using a null space penalization. To avoid overfitting, the spline fitting process of the main effect was restricted to 5 knots within the GAMM. We included the survey in the model as a random factor. In addition, the varIdent function was inserted as a weighting factor to allow for variance in fit between day and night. A Spearman’s correlation matrix was calculated to determine collinearity among the environmental variables. Variables used in the model were selected by using a cut off value of 0.80. The autocorrelation of residuals was modeled using a first-order autoregressive error structure nested within each deployment.

Results

Properties of Sound Scattering Layers

Acoustic backscatter intensity (log₁₀ NASC) varied among the three feature classes, time of day, and across the three depth intervals. Within the epipelagic (15–200 m), acoustic backscatter intensity was significantly greater during the night among all feature types (p < 0.001), with values greater by a factor of 25-fold in frontal boundaries, 17-fold in LCOW, and 12-fold in Common Water, relative to daytime values (Table 1). During both nighttime and daytime, the lowest backscatter occurred within the LCOW (Figure 3). Daytime backscatter increased significantly (p < 0.001) from the LCOW to frontal boundaries, and from frontal boundaries to Common Water. However, at night backscatter was greatest within the frontal boundaries and slightly less, although not significantly, (p = 0.258), within Common Water. Similarly, within the upper mesopelagic (200–600 m), acoustic backscatter increased significantly (p < 0.001) from LCOW to frontal boundaries, and from frontal boundaries to Common Water during the day; however, LCOW and frontal boundary features displayed significantly greater backscatter (p < 0.001) at night than Common Water (Figure 3). In contrast to the other two depth intervals, lower mesopelagic zone (600–1000 m) LCOW waters had significantly (p < 0.001) greater backscatter than the same depth interval in either frontal boundary or Common Water, with nearly a 9-fold and 19-fold increase relative to the latter two features at night, respectively (Figure 3).

TABLE 1

Table 1. Mean acoustic backscatter by time of day (TOD) and oceanographic feature class (LCOW, Loop Current-origin water; FB, Frontal boundary; CW, Common Water).

FIGURE 3

Figure 3. Least-squared means of acoustic backscatter, log₁₀ NASC (m² nmi^–2), by feature class, time of day (open symbols are day, filled are night) and depth interval. LCOW, Loop Current-origin water; FB, frontal boundary; and CW, Gulf Common Water. Error bars represent standard error from least-squared mean estimates.

In all cases, during the day the center of mass occurred significantly deeper in frontal boundaries than in LCOW, and in Common Water than frontal boundaries (p < 0.001; Figure 4); however, at night the responses were more variable. At night the center of mass within the LCOW was consistently and significantly (p < 0.001) deeper than within frontal boundaries and Common Water. Within the epipelagic zone the centers of mass of frontal boundaries and Common Water were significantly shallower than within LCOW (p < 0.001). The centers of mass within the upper mesopelagic occurred deepened going from LCOW to frontal boundaries to Common Water, respectively (p < 0.001) (Figure 4). The location of the centers of mass within the lower mesopelagic zone was the most variable at night, with frontal boundaries having a significantly shallower center of mass (716.7 m; p < 0.001) than either Common Water (737.7 m) or LCOW (766.2 m) (Figure 4).

FIGURE 4

Figure 4. Least-squared means of center of mass by feature class, time of day (open symbols are day, filled are night) and depth interval. LCOW, Loop Current-origin water; FB, frontal boundary; and CW, Gulf Common Water. Error bars represent standard error from least-squared mean estimates.

Effect of Environmental Drivers on Backscatter

Surface light intensity (IPAR; W m^–2) at 5 m water depth was not significantly different among the three oceanographic feature classes (Kruskal–Wallis One Way ANOVA on Ranks, p = 0.765); although in general LCOW stations displayed greater surface light intensity, followed by frontal boundaries and Common Water (Figure 5). In all GAMMs light intensity and the LCOW index emerged as the most consistently significant variables (Supplementary Table S1) among all three depth intervals. Temperature, chlorophyll and the Common Water index were significant variables with respect to backscatter in certain cases (Supplementary Tables S1, S2). With the exception of only a few stations in cruises DP01 and DP02, transects were seaward of the coastal production plume associated with the Mississippi River, where net primary production exceeded 700 mg C m^–2 d^–1 (Supplementary Figure S1). Below we discuss the significant interactions between the oceanographic index scores and light intensity with backscatter across the three depth intervals.

FIGURE 5

Figure 5. Mean surface light intensity (IPAR) at 5 m water depth for each of the three oceanographic feature classes surveyed during DP01–DP04. LCOW, Loop Current-origin Water; FB, frontal boundary; and CW, Gulf Common Water. Error bars represent standard error of the mean.

Within the epipelagic, the acoustic backscatter was correlated with surface light intensity, sea surface temperature, chlorophyll concentration, and the LCOW index, as well as the interaction between light intensity and the LCOW index (r² = 0.53; Supplementary Table S1). The partial plots indicated an increasing trend in backscatter at light intensity values < 200 W m^–2 (local night). As expected, backscatter decreased in the epipelagic as light levels increased beyond 200 W m^–2, ostensibly a function of DVM. In addition, backscatter decreased as surface temperature increased. Increased surface chlorophyll concentrations were associated with increased backscatter (Supplementary Table S1). As the intensity of the LCOW index increased in the epipelagic, we observed precipitous declines in backscatter, with a nearly 45% decline at night. In the epipelagic, only the LCOW index and light intensity interaction was significant (p < 0.001; Figure 6).

FIGURE 6

Figure 6. Heatmap representing acoustic backscatter (NASC) as a function of the oceanographic index and light intensity (IPAR₅_m). (Left panels) Represent the epipelagic (15–200 m), (center panels) represent upper-mesopelagic (200–600 m), and (right panels) represent lower mesopelagic (600–1000 m).

In the upper mesopelagic, light intensity was the most significant factor (p < 0.001) explaining variance in backscatter intensity, with lesser explanatory power attributed to the Common Water index (p = 0.004), the interaction of the Common Water index with light intensity (Supplementary Table S1). A decrease in backscatter was associated with increasing values of the Common Water index (p = 0.004) and with decreasing light intensity (p < 0.001) (Supplementary Table S1). During daytime and nighttime, a decrease in backscatter was observed with a greater LCOW index, whereas lower values of the Common Water index were associated with increased backscatter during the day (Figure 4). In general, the correlation of backscatter and light intensity within the upper mesopelagic was less variable than either epipelagic or lower mesopelagic depths.

Backscatter within the lower mesopelagic was significantly related to the velocity of the front, the LCOW index, and the interaction of light intensity in both LCOW and Common Water stations (Supplementary Table S1). The interaction of light intensity and the Common Water and LCOW indices were significant at depths greater than 600 m, with an increase in backscatter during the day at the lowest values of the Common Water index (Figure 6), and greater backscatter during the day and night periods for higher LCOW index (Figure 6).

Vertical Distribution of Backscatter

The GAM model for the center of mass indicated that within the epipelagic zone the vertical distribution of sound scattering layers were significantly related only to the maximum chlorophyll concentration (p < 0.001; r² = 0.16; Supplementary Table S2), suggesting a deepening of the center of mass as chlorophyll concentration increases.

The model for the upper mesopelagic was significant (p < 0.001; r² = 0.33) and selected six terms explaining variability in the center of mass: light intensity, temperature, salinity, maximum chlorophyll concentration, LCOW index and the interaction between the latter and the light intensity (Supplementary Table S2). The effect of light intensity on the center of mass indicated that biomass was deeper in the day and shallower at night, while an increase in sea surface temperature, salinity, and the LCOW index were associated with a deeper center of mass. In contrast, the center of mass was shallower with increases in the maximum chlorophyll concentration (Supplementary Table S2).

Within the lower mesopelagic, the center of mass was significantly related to the LCOW index (p < 0.001) in addition to the surface temperature (p < 0.001) and surface salinity (p = 0.003). A strong relationship was observed with the LCOW index, suggesting that as the LCOW index increases, the center of mass of layers gets deeper. The deepening of the center of mass occurred both during day and nighttime as indicated by the significant interaction between light intensity and LCOW index (p = 0.02) (Supplementary Table S2).

Biological Ground Truthing

Detailed information on the faunal composition, vertical distribution, and standing stocks of the epi- and mesopelagic fauna collected during DEEPEND surveys are reported elsewhere (Judkins et al., 2016; Burdett et al., 2017; Sutton et al., 2017; Frank et al., 2020; Judkins and Vecchione, 2020; Milligan and Sutton, 2020), but will be briefly summarized here. The two dominant taxonomic groups collected with the MOCNESS were fishes and macrocrustaceans (large euphausiids, decapod shrimps, mysids, lophogastrids). Macrocrustaceans were the most abundant group by number, contributing 26.5% of the total abundance of organisms collected, with euphausiids being numerically dominant (e.g., Nematoscelis atlantica, Stylocheiron abbreviatum, and Thysanopoda obtusifrons). The fish assemblage was dominated by the Order Stomiiformes (75% of all fishes collected; Sutton et al., 2020), particularly species of the genus Cyclothone which contributed 18.7% by number of all organisms collected (fishes and invertebrates). The species Cyclothone pallida accounted for over half (56%) of this dominant genus. Myctophid species often dominated the numbers of upper mesopelagic layers and the epipelagic layer at night (Sutton et al., 2017, 2020; Milligan and Sutton, 2020). Other taxa commonly collected in net samples included: gelatinous zooplankton (e.g., siphonophores, medusae, and pyrosomes), shelled pteropods, cephalopods, and a wide variety of “other fishes” (Aulopiformes, Stephanoberyciformes, early life stages of coastal and benthic taxa; Sutton et al., 2017, 2020). In relation to water mass, LCOW and Common Water stations had the greatest number of individuals from net collections, with 40 and 44%, respectively. Each had similar species composition, dominated by the macrocrustaceans and fishes mentioned above.

Discussion

Horizontally continuous sound scattering layers were ubiquitous in the northern GoM, throughout all periods and features surveyed in this study, occupying all three depth intervals examined. As expected, an increase in acoustic backscatter was observed during periods of low surface light intensity (<200 W m^–2) in the epipelagic, followed by a coincident decrease in upper mesopelagic backscatter at night due to the upwardly migrating mesopelagic assemblage. This pattern was observed consistently across all oceanographic feature classes.

Sound Scattering Layer Response to Oceanographic Features

In general, LCOW stations were characterized by the lowest backscatter intensity within the epipelagic, intermediate intensity in the upper mesopelagic, and greatest in the lower mesopelagic zone. An overall reduction in biomass within anticyclones has been reported across many systems, though the manifestations appear to be system-specific (Godø et al., 2012; Béhagle et al., 2014; Fennell and Rose, 2015; Gaube et al., 2018; but see Goldthwait and Steinberg, 2008). The reduction in lower trophic-level (e.g., zooplankton) biomass associated with anticyclonic features, which are similar in structure to LCOW (Johnston et al., 2019), has been observed previously in the northern GoM (Zimmerman and Biggs, 1999; Wormuth et al., 2000; Ressler and Jochens, 2003; Gasca, 2004) as well as in other low-latitude oceanic regions (Huggett, 2014; Lebourges-Dhaussy et al., 2014). In comparison, Godø et al. (2012) observed variation (∼20 dB re 1 m^–1) in backscatter while transiting across oceanographic discontinuities in the Icelandic Basin, including an anticyclonic eddy. However, they noted patchiness in the biomass estimated across those features. In contrast, Fennell and Rose (2015) demonstrated increased backscatter in mesopelagic sound scattering layers in the mid-North Atlantic Ocean associated with mesoscale anticyclonic eddies and attributed the increased backscatter to transport mechanisms associated with eddy fields. While many studies have noted that cold-core cyclonic eddies, characterized by centrally upwelled, nutrient-enriched water, promote increased primary and secondary productivity in the epipelagic (Biggs, 1992; Zimmerman and Biggs, 1999; Seki et al., 2001; Landry et al., 2008), during our survey periods we did not encounter any such biomass enhancement with respect to higher trophic levels due to the ephemeral nature and/or our undersampling of cyclonic features.

Given the pattern of reduced backscatter (and by proxy, micronekton biomass) in Loop Current and anticyclonic features, two ecological explanations can be proposed: (1) these features support less micronekton biomass as a quasi-self-contained habitat unit, or (2) these features influence avoidance behavior of vertical migrators. Of these explanations, we posit that the second is more likely. First, with respect to in situ production, the known generation times of micronekton are longer than the lifespans of shed Loop Current eddies (e.g., Gartner, 1991) and much longer than cyclonic eddies, which are smaller and more ephemeral than Loop Current eddies (Johnston et al., 2019). This differentiates results found for zooplankton and micronekton – the former can be “spun down” within the lifetime of an anticyclonic eddy due to lack of new production as food resources are exhausted. Second, with respect to spatial coherence of micronekton within a mesoscale feature, differential lateral advection from the surface to depth during vertical migration would be a diffusing agent (Milligan and Sutton, 2020). For micronekton to retain spatial coherence with a mesoscale feature, which is itself in motion, micronekton would have to “track” surface features during the daytime while at depth, although the extent of this effect would be dependent on how deep these features propagate at depth. The classical paradigm of daytime behavior of vertically migrating micronekton is that they are quiescent, conserving energy between feeding bouts, not actively tracking features geographically (Sutton, 2013).

Perhaps the reduction of vertical migration into the epipelagic and upper mesopelagic zones at night under LCOW results in a deep accumulation of biomass in the lower mesopelagic under anticyclonic-like conditions relative to the other two oceanographic features. This supposition would support the behavioral argument posited above. Given that the influence of anticyclonic features can be detected well into, and at times below, the 600–1000 m depth interval (Godø et al., 2012; Furey et al., 2018), it is not surprising to see a response in the sound scattering layers. In our case we detected not only an increase in biomass, but also deepening of the layers associated with LCOW.

The effects of frontal boundaries in surface waters can be highly variable with respect to spatial extent, intensity, and persistence (Belkin et al., 2009), and can therefore have variable effects on aggregating biological resources concentrated through entrainment (Owen, 1981) in addition to larger predators that exploit these hotspots (Bakun, 2006; Scales et al., 2014). The greatest acoustic backscatter was observed at night within frontal boundaries. While Loop Current eddies may be associated with reduced productivity, frontal margins of these features are known to be sites of increased faunal abundance/biomass, for example larval and juvenile fishes (Mohan et al., 2017), and may potentially offset the reduced faunal abundance/biomass observed in the adjacent LCOW.

At some sites, backscatter intensity was not conserved between day and night sampling periods. There are multiple explanations for why this may occur, including the advection of organisms into or out of the study region between sampling intervals or organisms’ target strength varying as a function of depth (which is highly probable for resonant scatterers such as swim-bladders in fish). The fact that backscatter varied over diel periods at some sites and not as much at others presents a challenge in terms of interpreting acoustically measured biomass for migrating organisms, unless light cycle (or time of day) is controlled for in the analysis.

Variation in Depth Distribution

An increase in depth of the sound scattering layers as a function of anticyclonic physics is consistent with downwelling processes characteristic of these features (Carton et al., 2010; Chelton et al., 2011). The estimates of the depth of the center of mass suggest that in LCOW, the sound scattering layers are distributed at greater depths at night and that important biophysical interactions may influence the vertical distribution of the layers through either direct action (i.e., downwelling of migrators and/or their planktonic food) or through influences on individual behavior (migration choice) (Pearre, 2003). Moreover, the interaction between the LCOW index and backscatter suggests that an upper threshold in the oceanographic conditions (indicated by a high LCOW index) might mediate how organisms move into the epipelagic at night (as illustrated in Figure 6). This deepening of the sound scattering layers may indicate that organisms inhabiting these features may remain deeper to avoid the dynamic Loop Current waters, possibly because the current’s hydrodynamics add an additional energetic burden that could contribute to a reduction in the vertical movement of the mesopelagic assemblage. Alternatively, the persistent sound scattering layer detected at depth during the day and night may be attributed in part to both asynchronous migration strategies and non-migrators that continuously remain at depth (Sutton and Hopkins, 1996a; Watanabe et al., 1999; Olivar et al., 2012; Sutton, 2013). The phenomenon of asynchronous migration has been observed across many oceanic systems (Clarke, 1974; Badcock and Merrett, 1976; Kenaley, 2008) including the GoM (Sutton and Hopkins, 1996a, b) where dragonfishes (Stomiidae), the dominant mesopelagic predatory fishes, split their time between the epipelagic and mesopelagic depth intervals at night (Sutton and Hopkins, 1996b). Summaries of MOCNESS data (Sutton et al., 2017; Milligan and Sutton, 2020) indicate that these migration strategies are commonplace among the dominant GoM fish species and that is likely to explain at least in part, the persistent sound scattering layers observed at depth.

Influence of Light Regimes

Light level consistently correlated with temporal patterns in the mesopelagic sound scattering layers. We observed predictable patterns in the way the mesopelagic assemblage responded, with consistent increases in backscatter at night in the epipelagic and mesopelagic. Other studies have demonstrated that light intensity is important for controlling the extent of vertical movement and timing (Frank and Widder, 1997; Aksnes et al., 2017; Kaartvedt et al., 2017). While we were unable to empirically measure the light intensity at depth during this study, our results suggest that predicted light intensities from surface-derived estimates can be used as a predictor variable when direct measurements are unavailable. While others have quantitatively examined the animal response to subtle changes in light intensity (Frank and Widder, 1997; Aksnes et al., 2017; Kaartvedt et al., 2017), our estimates are appropriate to examine an integrated timescale that broadly represents the patterns observed with the DVM of the mesopelagic community. Additionally, while we were able to derive satellite-based estimates of light intensity across the three types of water masses examined, it remains unknown how variability in water transparency in LCOW, frontal boundaries and Common Water might differentially mediate light transmission to the mesopelagic region and in particular, whether the estimates we derived could result in significant differences among the three oceanographic feature classes we examined off the continental shelf.

Implications for Trophic Transfer

While we observed differences in backscatter distribution and intensity across the oceanographic conditions studied, the composition of organisms was not substantially different with net catches being dominated by crustacean macrozooplankton, namely euphausiids (e.g., Nematoscelis atlantica, Stylocheiron abbreviatum, and Thysanopoda obtusifrons), as well as fishes dominated by Cyclothone spp., dominant myctophid species (Sutton et al., 2020; Milligan and Sutton, 2020), other gonoistomatids, and hatchetfishes. An exception was noted for species composition of net hauls at frontal boundary stations, where euphausiids and pteropods were more abundant than fishes (Sutton et al., 2017). This suggests that the variation observed in backscatter is likely attributed in large part to the changes in organismal abundance and vertical distribution rather than a significant change in assemblage structure.

The differences in vertically migrating biomass among oceanographic feature types can have important implications for mediating the strength of trophic interactions and ultimately carbon transfer. Hopkins et al. (1996) estimated that 80% of all trophic exchange within the upper 1000 m of the water column in the eastern GoM (within our study area) occurs within the epipelagic zone at night. These authors determined that this consumption was driven by three dominant fish families (Myctophidae, Sternoptychidae and Gonostomatidae). Based on this study, the reduction in backscatter in LCOW suggests that the Loop Current and its associated eddies are likely areas of reduced trophic exchange, which has important implications for spatially explicit models of the GoM, and by proxy, other large marine ecosystems. Given that Loop Current eddies are persistent and dominant features within the GoM, occupying 100’s of square kilometers and with lifespans exceeding a year, the systematic reduction in trophic transfer likely decreases carbon sequestration by the system as a whole (Volk and Hoffert, 1985; Irigoien et al., 2014; Davison et al., 2015). As reported by Volk and Hoffert (1985), nearly 70% of carbon transport in the upper 1000 m is mediated by the biological pump due to the vertically integrated food web, consequently transporting surface production into the deep ocean (Ducklow et al., 2001; Irigoien et al., 2014; Ariza et al., 2015; Davison et al., 2015).

Conclusion

We demonstrate that Loop Current-origin waters in the upper GoM are associated with decreased acoustic backscatter in comparison to the other oceanographic feature classes examined. These patterns were temporally consistent, suggesting that this oceanographic milieu dampens vertical migration by the mesopelagic assemblage. Perhaps equally important, we demonstrate that within the adjacent frontal boundaries along the margins of Loop Current eddies, increased backscatter was measured in both the epipelagic and mesopelagic at night, and we speculate that this enhancement may offset some portion of the reduced standing stocks observed within the nearby LCOW features.

Physical forcing is an important process that operates at various temporal and spatial scales and can act to structure the distribution of organisms and therefore their roles in ecosystems. In this study, we show that in addition to the previously reported relationship in sound scattering layer dynamics relative to light levels (Røstad et al., 2016; Aksnes et al., 2017), we observed a quantifiable correlation with mesoscale oceanographic features, and the nature of this latter correlation is depth-stratum-specific. Given the ubiquitous distribution, immense biomass, and critical role that migrating sound scattering layers play in the global biological pump, understanding how oceanographic processes mediate the distributional patterns of billions of tons of mesopelagic micronekton is necessary to refine global carbon models (sensu Proud et al., 2017). This is especially true of low-latitude, deep-pelagic ecosystems, which are by far the largest component of the World Ocean. With increases in ocean temperature, and associated ecosystem changes (e.g., expanding oxygen minimum zones; Aksnes et al., 2017), particularly in the marginal seas, approaches that leverage the benefits of large-scale observational techniques with fine-scale, process-based methods provide an efficient means to examine the physical dependencies on biological organization in oceanic systems.

Data Availability Statement

Data are publicly available through the Gulf of Mexico Research Initiative Information and Data Cooperative (GRIIDC) at https://data.gulfresearchinitiative.org (HYCOM 10.7266/N7QR4VK0; Classifications from Johnston et al. (2019) 10.7266/N7QR4VK0; DP01 acoustic data: 10.7289/V5WD3XJS; DP02 acoustic data: 10.7289/V5RN35TF; DP03 acoustic data: 10.7289/V5WD3XSX; and DP04 acoustic data: 10.7289/V5RN3620).

Author Contributions

KB and JW collected the data. KB, MD’E, JM, and JW performed the analysis. All authors contributed to data interpretation and manuscript writing.

Funding

This research was made possible by a grant from The Gulf of Mexico Research Initiative.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank the captain and crew of the R/V Point Sur for their incredible support during the DEEPEND research program and at-sea operations. We appreciate the efforts of the National Centers for Environmental Information (www.ncei.noaa.gov) for assistance with archiving the acoustic data from all the DEEPEND cruises. We thank B. Penta and S. deRada at the Naval Research Lab for access to the high-resolution GoM HyCOM model output as well as C. Hu and the Optical Oceanography Laboratory at the University of South Florida for light intensity data included in the analysis. We thank S. Labua and D. Correa for assistance with data organization and processing and appreciate the comments of the two reviewers that improved this manuscript. This is contribution #176 from the Division of Coastlines and Oceans in the Institute of Environment at Florida International University.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2020.00051/full#supplementary-material

Footnotes

^ Publicly available at: http://hycom.org.
^ https://www.esrl.noaa.gov/gmd/grad/solcalc/azel.html
^ http://www.science.oregonstate.edu/ocean.productivity/

References

Aksnes, D., Dupont, N., Staby, A., Fiksen, O., Kaartvedt, S., and Aure, J. (2009). Coastal water darkening and implications for mesopelagic regime shifts in Norwegian fjord. Mar. Ecol Press Ser. 387, 39–49. doi: 10.3354/meps08120