Between February 21 and March 6, 2009 slope sediment samples were collected during a cruise on board of R/V Atlantis (Leg AT 15-44) off Costa Rica’s Pacific margin. This was complementary to a larger study focused on macrofauna inhabiting methane seep carbonates. The area is characterized by an intense, well-developed, mid-water OMZ which impinges on the margin at depths between 300 and 700 m (Levin et al., 2015). Along the Costa Rica margin there is a high degree of local and regional heterogeneity caused by carbonates, submarine fluid venting sites, methane seeps, hydrothermal seeps, mud mounds, seamount subduction scars, and slope failures as well-strong vertical hydrographic gradients (Han et al., 2004; Mau et al., 2006; Sahling et al., 2008; Füri et al., 2010; Levin et al., 2012, 2015). Surface net primary production (NPP) is relatively high, as a result of macronutrient supply associated with thermocline shoaling, mixing and wind-driven upwelling, and exhibits intra-annual variability with highest primary production in September–October and lowest in January–February (Nielsen-Muñoz and Quesada-Alpízar, 2006; Pennington et al., 2006). Thus, the OMZ along with both topography and hydrography result in strong gradients in bottom-water oxygen concentration, temperature, and pH over a relatively short depth range (Levin et al., 2015).

Sediment samples were collected with a multicorer (9.6 cm core liners i.d.) at six open slope sites (400, 620, 800, 1000, 1200, and 1800 m water depth; Figure 1). At each depth, three independent multicorer deployments were conducted (except the 1200 m site n = 2). Only undisturbed cores with clear overlying water were used. Subsamples were taken simultaneously from the same corer, for meiofauna and sediment properties, using Perspex tubes (3.6 cm inner diameter, 10 cm⁻²). All sediment cores were sectioned vertically on board ship at 0–1, 1–2, 2–3, 3–5, and 5–10 cm. Sections were preserved in 4% buffered formalin (Pfannkuche and Thiel, 1988) for meiofauna analysis, while fractions for sediment properties were frozen until analysis. Water content (WC) was determined by weight loss after freeze-drying sediment aliquots. Sediment total organic carbon (TOC) and total nitrogen (TN) were determined from freeze-dried, homogenized sediment using a Costech CHN analyzer (Valencia, CA, United States) after removal of inorganic carbonates. Chl a and phaeopigments were determined spectrophotometrically from freeze-dried sediment (Hagerthey et al., 2006) after extraction with 90% acetone (Danovaro, 2010). Chloroplastic pigment equivalent (CPE), i.e., the sum of chl a and phaeopigments was used as a measure of OM of phytodetrital origin. The ratio chl a to phaeopigments was used as a proxy of ‘freshness’ of the OM derived from primary production (García and Thomsen, 2008). The chl a to TOC ratio was used to indicate bioavailability of the bulk organic matter (OM; i.e., reflecting contribution of organic carbon of phytoplanktonic origin, Ramalho et al., 2014). Net primary production values for study sites were provided by M. Kahru (SIO) and were estimated from satellite data of chl-a concentration, photosynthetically active radiation and sea-surface temperature using an empirical algorithm (Kahru et al., 2012). Hydrographic data including temperature, salinity, dissolved oxygen, and pH were obtained from CTD casts.

Metazoan meiofauna were extracted from sediment by suspension-decantation (Pfannkuche and Thiel, 1988) using a 40 μm-sieve as the lower limit. The procedure was repeated at least five times for the 1-cm sections, 8–10 times for the larger fractions. Residual sediment was examined for extraction efficiency (>98%, Neira et al., 2013). All meiofaunal organisms were counted and classified to higher taxon level under a stereoscopic microscope after staining with Rose Bengal (0.5 g L⁻¹) (Pfannkuche and Thiel, 1988; Giere, 2009). Nauplii were counted separately given their ecological differences in size and diet compared to adults and copepodites (Decho and Fleeger, 1988).

Meiofaunal density, taxon richness (number of higher taxa) and diversity indices [Pielou’s evenness (J′), Shannon–Wiener (H′ log10) and Rank 1 dominance (R1D, the proportion of the most abundant taxon, which is an inverse indicator of assemblage evenness)] were used to describe meiofauna assemblage structure. Taxon richness was also examined through rarefaction curves (number of taxa per number of individuals sampled – Hurlbert, 1971).

Meiofauna community composition (multivariate, square root transformed, Bray–Curtis Similarity) and total meiofauna abundance (univariate, square root transformed, Euclidean Distance) were analyzed with PERMANOVA tests (three-factor design: water depth, sediment depth, core). The random factor ‘core ‘was inserted as nested in water depth to account for the dependency of the sediment layers within each core; this approach has been used in several studies as a good solution to the sediment layer dependence issue (e.g., Ingels and Vanreusel, 2013; Rosli et al., 2016). This design allows us to investigate differences between water depths and sediment layers simultaneously. Additional PERMANOVA tests were performed on the meiofauna community matrix and the individual diversity indices (based on pooled 0–10 cm sediment data). Pairwise tests were performed for main factors and interactions when significant results were obtained. PERMDISPs were performed when PERMANOVA returned significant results, to investigate heterogeneity in dispersions. Monte Carlo tests were applied when the number of available permutations was <100. Non-parametric multidimensional scaling ordination (nMDS, square root transformation, Bray–Curtis Similarity) was used to visualize the dissimilarity in meiofaunal community structure among sites. Similarity percentage (SIMPER) analysis was used to determine which taxa were responsible for observed differences (Clarke, 1993).

Environmental data were subjected to univariate PERMANOVA analysis, testing differences between water depths for each variable separately (no transformation, Euclidean distance). Spearman-rank correlations were performed to examine collinearity between environmental variables, and possible relationships between environmental variables and meiofaunal taxa abundance. Redundancy analysis (RDA) was used to examine correlations between meiofaunal community composition and environmental variables (Ter Braak and Šmilauer, 2012) in order to assess the extent of the contribution of each variable to variance in community composition. Prior to RDA, a forward stepwise selection procedure was performed to minimize redundancy of variables (between-variable correlations) and select a set of explanatory environmental variables that account for the composition variation (Lepš and Šmilauer, 2003). This was followed by a Monte Carlo permutation test to determine significant relationships between the taxa and environment.

All above analyses were performed with PRIMER v6 (Clarke and Gorley, 2006) and the PERMANOVA+ add-on (Anderson et al., 2008), except for the Spearman Rank correlations (JMP Pro 12.0.1., SAS Institute, Inc.) and the RDA ordination and gradient analysis [CANOCO for Windows 5 (Ter Braak and Šmilauer, 2012)].

Environmental and sediment properties are summarized in Table 1. Bottom-water temperature decreased progressively with depth along the entire transect, while bottom-water dissolved oxygen (DO) values progressively increased with depth, with a minimum of 0.05 ml L⁻¹ at the 400 m site, and maximum of 1.61 ml L⁻¹ at the 1800 m site. Bottom-water pH fluctuated between 7.68 and 7.77 along the transect. The OMZ 400 (core) and 620 m sites were enriched in TOC (∼2.5%DW) and CPE (147 μg g⁻¹) compared to the deepest oxygenated site (∼1.1%DW, 8.2 μg g⁻¹, respectively). The chl a value at the OMZ core site (400 m) was very high, as was the phaeopigment concentration, leading to a lower ‘freshness’ (chl a/phaeo) value than was observed at 800 m and deeper.

Quality of available OM was least refractory at the deepest site (1800 m), with lowest C/N (7.7) and highest chl a/phaeo (1.54). The vertical profile of chl a/phaeo values along the transect is shown in Figure 2. Sediment environmental variables were significantly different between study sites (PERMANOVA; Supplementary Table S1). Surface water chl a has been considered an indicator of the flux of phytodetrital OM to the seabed and thereby food availability to benthic animals (Rex and Etter, 2010). The relatively high and homogeneous NPP in surface waters for the study area (271–332 mg C m⁻² d⁻¹) calculated for the sampling period (Table 1) was reflected in the partially degraded phytopigments found in sediments. Indeed, sediment C/N ratio compared with δ¹³C signatures indicate that deposited sedimentary OM is predominantly of marine origin (Figure 3).

Several correlations between environmental variables were significant (Supplementary Table S2), but most prominent were the significant relations between DO, bottom-water temperature, and chl a. Bottom temperature was inversely related to DO (r = −1.000), and positively to chl a (r = 1.000). In general, variables related to food quantity, freshness, and availability correlated strongly with each other. As a result, DO correlated significantly, negatively with TOC, TN, TOC/TN, Phaeo, CPE, chl a/TOC and positively with chl a/Phaeo.

In the OMZ core (400 m), most meiofauna was concentrated in the upper 2 cm, almost divided equally between the 0–1 and 1–2-cm horizons, while densities at the adjacent OMZ site (620 m) peaked in the 1–2 cm layer, but with relatively high numbers at the surface (0–1 cm) (Figure 8). At the deeper, more oxygenated sites, however, meiofauna are mostly concentrated in the top 1 cm, declining consistently with depth in the sediment (Figure 8). In strongly oxygen-depleted sites (∼0.05–0.16 ml L⁻¹) about 20–40% of meiofauna abundance is concentrated in the uppermost 1 cm, while at deeper slope sites (>800–1800 m water depth), 55–70% of meiofauna were found in the top 1 cm sediment layer (Figure 9A). The relative importance of copepods and nauplii increased (and nematodes decreased) toward the deeper (620 to 1800 m), more oxygenated sites (Figures 4, 9B). The PERMANOVA results on density and community composition reported above support these observations statistically, with sediment layer differences being more evident at the OMZ core and adjacent sites (620 and 800 m) where abundances do not significantly change in the upper fraction (∼0–3 cm). In contrast, the more oxygenated sites (1000 and 1800 m) show significant differences between the 0–1 cm fraction and deeper layers, reflecting the higher abundances concentrated at the sediment surface (Supplementary Tables S3, S4).

Spearman’s correlation coefficients are shown in Supplementary Table S6. Copepod and nauplii densities were positively correlated with DO (p < 0.05) and freshness of the OM (chl a/phaeo) (p < 0.05) (Figure 9B). At the same time, their densities were negatively correlated with chl a and phaeo (p < 0.05) and bioavailability of the OM (chl a/TOC) (p < 0.001). Polychaetes and other taxa showed no significant relationships with the measured environmental variables. Nematodes on the other hand correlated significantly and positively with sediment TOC (p < 0.05).

The forward stepwise selection procedure identified a suite of environmental variables that in combination offer the best explanatory value for meiofauna community structure (cf. RDA ordination model visualized for meiofauna community in Figure 10). Those environmental variables that were selected are shown in Supplementary Table S7 and include chl a/TOC, phaeopigments, chl a, water depth, and DO. The RDA analysis yielded four axes that explained 98.02% of the variance between meiofaunal community structure and environmental variables, which were significant (p = 0.002) for all canonical axes. The direction of vectors indicates that OM related descriptors (chl a, phaeo, chl a/TOC) increased along the first axis, and are associated with presence of nematodes, which are markedly dominant at the OMZ core site. Copepods, nauplii as well as other crustaceans such as amphipods, ostracods and cumaceans, mostly occur in the uppermost sediment layer, where they are likely associated with increased DO availability. Other taxa such as kinorhynchs, tanaids, and oligochaetes do not appear to be associated with these environmental variables over the range tested (Figure 10).

Sediment TOC levels were not significantly different within the different OMZ sites (400-m core, 620-m adjacent site, and 800-m lower boundary site) but they were higher in the OMZ compared to those at the deeper, more oxygenated sites (Supplementary Table S1). The strong inverse correlation between TOC (and other food-related variables) and DO follows a similar trend to those observed previously in other OMZs, such as off Peru (Neira et al., 2001b), off central Chile both shelf (Sellanes and Neira, 2006) and slope (Veit-Köhler et al., 2009; Neira et al., 2013), and in the Arabian Sea (Gooday et al., 2000; Levin et al., 2000). This strong correlation suggests that DO may play a primary role in controlling TOC distribution on the Costa Rica margin, potentially through regulation of biological activity resulting in variable consumption and degradation of TOC. Also, very high levels of sinking OM may lead to rapid consumption and aerobic degradation (Stukel et al., 2016), quickly consuming the available oxygen in low-ventilation waters and leaving relatively large amounts of OM in the sediment under hypoxic conditions where anaerobic degradation by microbial organisms takes over (McKew et al., 2013). In such a scenario other factors such as bioturbation, burial, sediment sorption, as well as sediment transport, sedimentation rates and local topography may additionally influence margin TOC distributions (Cowie and Levin, 2009; Cowie et al., 2009).

The TOC levels (∼2–2.5%DW) of the OMZ sediments were comparable to those reported for the Golfo Dulce (Costa Rica), a nearshore oxygen-depleted basin (∼100–200 m depth) located southeast of the present study region (Nichols-Driscoll, 1976; Thamdrup et al., 1996) and Panama Basin (3900 m) (Aller et al., 1998), but lower than in other eastern Pacific OMZs such as off Chile (∼122–309 m) (Neira et al., 2001c, 2013; Veit-Köhler et al., 2009). However, at the Costa Rica margin, C/N (molar) ratios were much lower (7.7–9.7; Table 1) and similar to more oceanic waters with OM of phytoplanktonic origin, compared to the nearshore Golfo Dulce (6.1–12.9), where terrestrial OM input with enhanced content of N-poor detritus occurs (Thamdrup et al., 1996). Higher C/N ratios can be indicative of preferential use of nitrogen during the remineralization of OM of phytoplanktonic origin (Sánchez and Carriquiry, 2007). The relatively high chl a/phaeo values at deeper sites suggest that material deposited is more fresh relative to that of OMZ sites (400 and 620 m). Furthermore, C/N ratios and carbon isotopic composition shows little spatial variability between sites (Figure 3). This suggests that little degradation of OM or at least only partial degradation occurs in the overlying water column (Cowie et al., 2009), and that the deposited OM is of marine origin, with δ¹³C signatures typically between −18 and −22^o/_oo and C/N ratios of 6–11. Spatial variability of C/N ratios among sites reflects the diagenetic state, reactivity and the quality of the particulate OM as well as bioturbation and bio-irrigation by fauna (Cowie, 2005; Cowie et al., 2009; Lessin et al., 2018).

Metazoan densities along the outer shelf-slope margin generally decrease with water depth and surface primary productivity (Rex et al., 2006). Table 5 summarizes information on mean total meiofauna abundances, number of taxa, and relative abundance of nematodes for the Costa Rica margin and compares it with studies conducted in other regions along the eastern Pacific. The mean densities registered at the OMZ core site (400 m) off Costa Rica (∼3688 ind. 10 cm⁻²) are one the highest compared to other OMZ core sites (independent of water depth) along the eastern Pacific, such as off Concepción (< 2000 ind. 10 cm⁻²) (Neira et al., 2001c, 2013; Sellanes et al., 2003; Sellanes and Neira, 2006), off Callao, Peru (∼1500 ind. 10 cm⁻²), and almost double that reported for an OMZ off Antofagasta (Veit-Köhler et al., 2009). The drastic decline in total abundance of meiofauna from the OMZ core toward the adjacent site (620 m) with increasing oxygen, but still within the OMZ, was a pattern similar to that observed at the Peru margin (Neira et al., 2001b).

Marine benthic diversity usually increases with increasing water depth along the upper continental slope (Menot et al., 2010; Rex and Etter, 2010). In this study, a total of 19 meiofaunal taxa were found, but at the OMZ core site, only two taxa were collected (Nematoda and Polychaeta), with nematodes completely dominating meiofaunal composition (99.9% at the 400 m site, Table 2). This low number was even lower than reported for other Pacific OMZ sites such as off Antofagasta (4 taxa, Veit-Köhler et al., 2009), central Chile (4–8 taxa, Neira et al., 2001c; Sellanes and Neira, 2006) and, off Peru (7 taxa, Neira et al., 2001b) (Table 5). Differing diversity between OMZ vs. non-OMZ sites was also clearly reflected in diversity indices (Table 4 and Supplementary Table S5). Studies off Peru suggest that beyond the OMZ habitat, sediment heterogeneity has more relevance in shaping metazoan meiofaunal diversity compared to within the OMZ core, where oxygen is the main limiting factor (Neira et al., 2001b; Gooday et al., 2010; Sellanes et al., 2010). Oxygen availability in sediments plays a key role in controlling meiofaunal diversity (Gooday et al., 2010; Zeppilli et al., 2015). Reduced conditions may negatively impact meiofaunal diversity (e.g., Neira and Rackemann, 1996), and nematode community structure and trophic diversity (Ingels et al., 2011; Neira et al., 2013). Interestingly, this is illustrated by the changing diversity of Foraminifera from a site off Pakistan, showing consistently higher diversity during the inter-monsoon season (when the study site was outside the OMZ) compared to during the monsoon season (when the site was within the OMZ) (Larkin and Gooday, 2009; Gooday et al., 2010). Similarly, off Pakistan macrofaunal (>300 μm) and megafaunal diversity was strongly reduced where DO was lowest, but diversity recovered rapidly at slightly higher oxygen levels (Gooday et al., 2009; Levin et al., 2009). These observations are likely attributable to an ‘edge effect’ where a release from oxygen stress combined with an abundant food supply occurs (Murty et al., 2009; Gooday et al., 2010). Veit-Köhler et al. (2009) reported 10 taxa at an OMZ site off Chile that had slightly elevated oxygen levels (0.79 ml L⁻¹) compared to 4 taxa from the lower-oxygen sites, suggesting comparable processes may be taking place in that area. For the Costa Rica margin, this is reflected by the increase in the number of meiofauna taxa from 2 (400 m) to 9 (620 m) within a gradient of DO from 0.05 to 0.16 ml L⁻¹ (Table 5), and likely related to the relatively low C/N values and higher food quality (resulting from fresher OM, i.e., higher chl a/phaeo) toward the more oxygenated sites (Figure 2 and Table 1).

Meiofaunal community distribution is influenced by numerous factors such as physical sediment properties (Giere, 2009), interactions with macrofauna (Pfannkuche and Thiel, 1987; Lambshead et al., 1994; Van Colen et al., 2009) or microbiota (Danovaro et al., 1995, 2000), habitat heterogeneity (Van Gaever et al., 2009), and oxygen availability (Levin et al., 2002; Sellanes and Neira, 2006). The potential role of oxygen is reflected in the RDA plot (Figure 10) where oxygen is positively correlated with a number of taxa. In fact, there are distinctly divergent meiofaunal responses across OMZ-influenced settings compared to deeper, more oxygenated ones (Figures 6, 7). The most oxygen-depleted sites were characterized by lower meiofaunal taxon richness, lower diversity (H′) and lower taxonomic evenness (J′), but higher rank 1 dominance (R₁D) (Table 4). The more similar evenness (J′) and diversity (H′) exhibited at deeper, better-oxygenated sites suggests a more equitable distribution of meiofaunal assemblages relative to the OMZ core site (400 m) (Figure 6 and Table 4). Nematode densities at our Costa Rica sites were enhanced at low-oxygen sites, consistent with other OMZs studies (e.g., Levin et al., 1991; Cook et al., 2000; Neira et al., 2001b,c, 2013; Sellanes and Neira, 2006). The very high densities and total numerical dominance of nematodes at the site with lowest oxygen concentration, and the proportional decline of nematode abundance with increasing oxygen (Figure 9B) supports previous observations of a potential positive effect of low-oxygen levels on nematode abundance (e.g., Neira et al., 2001b), including the indirect effect of decreasing hypoxia-intolerant predators and competitors (Neira et al., 2001b; Levin et al., 2002; Levin, 2003). Time series data at a shelf site (98 m) off Callao, Peru has shown fluctuations in nematode density and a community 100% dominated by nematodes. Such extreme dominance, termed a “nematode state” was suggested to take place during extreme La Niña conditions, when bottom-water oxygen concentrations are at their lowest and macrofaunal densities are strongly reduced (Gutiérrez et al., 2008). This is likely reflected by the consistent decline in abundance of meiofaunal-size (and larger) polychaetes toward the most oxygen-depleted sites (Figure 4). Also crustaceans, such as copepods and nauplii, isopods, amphipods and ostracods are sensitive to low-oxygen conditions (Josefson and Widbom, 1988; Wetzel et al., 2001; De Troch et al., 2013) and compete with nematodes for space and resources, so their absolute and relative decline with decreasing oxygen levels likely allows nematodes to dominate the sediment system. Copepods are active swimmers and their metabolism requires more oxygen than other meiofaunal groups (Moodley et al., 1997; Modig and Olafsson, 1998; Giere, 2009). Summer hypoxic conditions (1985/86) on the Louisiana shelf resulted in the loss of copepods for a month, with moderate loss of nematode and kinorhynch densities. After the habitat returned to normal oxygen conditions, the copepods recovered slowly (Murrell and Fleeger, 1989). Oxygen depletion may lead to a feeding stop and unsuccessful reproduction in adult copepods (Veit-Köhler et al., 2009, 2018; De Troch et al., 2013) and hence absence of nauplii and adult copepods at the OMZ core. Recovery, speed and extent of colonization of meiofauna after hypoxic events are variable. For instance, the high mobility of copepods and ostracods promotes a rapid recolonization of azoic sediments, compared to nematodes which require more time to cover similar distances (Chandler and Fleeger, 1983). Experiments in the Wadden Sea reduced sediments (“black spots”) showed that meiofauna recolonized sediments only after redevelopment of an oxidized sediment layer (Neira and Rackemann, 1996). Van Colen et al. (2009) found strong changes in nematode community composition after exposing them to experimentally-induced hypoxia, but this did not result in complete nematode mortality. For instance, under El Niño conditions, which cause enhanced bottom-water oxygen concentrations, there was an increased representation of copepods and nauplii by one order of magnitude at sites influenced by the OMZ (Neira et al., 2001c). We suspect that the most important driver of meiofaunal communities is bottom-water oxygen, but that the co-occurrence of low oxygen and chl a/phaeo, and high chl a, phaeo, and chl a/TOC (and hence food quality and availability) lie behind some of the causality of the correlations.

The quality of food resources, in addition to oxygen, can play an important role in structuring faunal communities. For meiofauna, the role of food availability has been highlighted in many studies as one of the most important factors structuring communities and spatial biodiversity patterns (Danovaro et al., 1995; Neira et al., 2001b; Giere, 2009; Pusceddu et al., 2009). The importance of food quantity and availability on nematode abundance and diversity in OMZs has also been suggested by Cook et al. (2000) for the Arabian Sea. Food availability is a regulating factor along the vertical sediment profile (Ingels et al., 2009, 2011; Ingels and Vanreusel, 2013). At the Costa Rica margin, meiofauna densities declined with increasing depth in the sediment, a recurrent gradient observed in marine sediments (Levin et al., 2002; Soltwedel et al., 2005). However, at the lowest oxygen sites of Costa Rica (400 and 620 m), meiofauna, dominated by nematodes, were concentrated in subsurface sediments (∼41–45% in the 1–2 cm layer); whereas toward deeper sites, the surface 1-cm layer harbors most of the meiobenthos (Figure 9A). The vertical distribution of meiofauna was similar to that observed at the OMZ off Callao (transect 305–1210 m, Neira et al., 2001b), where densities declined consistently with depth in the sediment only at the oxygenated sites (Figure 8). The high meiofaunal densities as well as the high dominance of nematodes at the OMZ sites, and their high subsurface densities, are most likely attributable to the extensive availability of food, and relatively fresh phytodetritus (Figure 2) even in the deeper sediment layers (e.g., CPE ranging from 147 μg g⁻¹ at 0–1 cm to 88 μg g⁻¹ at 5–10 cm depth at site 400 m), as well as the relative tolerance of nematodes to low-oxygen levels (Neira et al., 2001a, 2005; Neira and Decraemer, 2009; Sajan et al., 2010), and reduced competition for food by exclusion of the more oxygen-sensitive macro- and meiofauna (Neira et al., 2001b; Levin et al., 2002; Sajan et al., 2010; Jeffreys et al., 2012). Macrofaunal activity may also influence vertical meiofauna patterns in the sediment through bioturbation. For instance, tuna crab (Pleuroncodes sp.), common on the Central American margin (Pineda et al., 2016), can sift out the food material present in the substrate and extract small animals such as nematodes, and other invertebrates (Boyd, 1962, 1967; Lovrich and Thiel, 2011). During Alvin dives in 2017, large masses of galatheid crabs were observed apparently feeding on the sediment bottom of Quepos slide, an active seep area (bacterial mats) influenced by the OMZ and nearby our 400 m site (L. A. Levin personal observation). Galatheid crabs such as Pleuroncodes spp. have often been seen associated with oxygen-depleted waters (Sellanes et al., 2007; Lovrich and Thiel, 2011; Gallardo et al., 2017) and can form dense aggregations swarming the seafloor (Thiel and Lovrich, 2011; Pineda et al., 2016). Bioturbating macrofauna, mainly annelids tolerant to hypoxic conditions found at OMZs, including at very low oxygen sites (e.g., Levin et al., 2000, 2002; Smith et al., 2000; Sellanes et al., 2007), may induce mixing downward, advecting “fresh” and nutritional material deep into the sediment (Levin et al., 1997; Gutiérrez et al., 2000), which, in turn, may support higher nematode survival rates in subsurface and deeper sediments (Braeckman et al., 2010).