Holopelagic Sargassum aggregations provide warmer microhabitats for associated fauna

Drifting aggregations of Sargassum algae provide critical habitat for endemic, endangered, and commercially important species. They may also provide favorable microclimates for associated fauna. To quantify thermal characteristics of holopelagic Sargassum aggregations, we evaluated thermal profiles of 50 aggregations in situ in the Sargasso Sea. Sea surface temperature (SST) in the center of aggregations was significantly higher than in nearby open water, and SST differential was independent of aggregation volume, area, and thickness. SST differential between aggregation edge and open water was smaller than those between aggregation center and aggregation edge and between aggregation center and open water. Water temperature was significantly higher inside and below aggregations compared to open water but did not vary inside aggregations with depth. Holopelagic Sargassum aggregations provide warmer microhabitats for associated fauna, which may benefit marine ectotherms, though temperature differentials were narrow (up to 0.7 °C) over the range of aggregation sizes we encountered (area 0.01–15 m2). We propose a hypothetical curve describing variation in SST differential with Sargassum aggregation size as a prediction for future studies to evaluate across temporal and geographic ranges. Our study provides a foundation for investigating the importance of thermal microhabitats in holopelagic Sargassum ecosystems.


Methods
Study area and sampling design.We evaluated thermal profiles of 50 Sargassum aggregations by assessing temperature differentials between aggregations and nearby open water during a research expedition to the Sargasso Sea from 30 July to 12 August 2019 aboard the Greenpeace vessel MV Esperanza.After departing from Bermuda, vessel headings were set using ocean surface temperature and ocean color imagery obtained from Roffer's Ocean Fishing Forecasting Service (https:// www.roffs.com) to target convergence zones and increase the probability of encountering Sargassum aggregations.All sampling occurred southeast of Bermuda (30.3-31.5°N,61.7-63.9°W)within the boundaries of the Sargasso Sea (22-38°N, 43-76°W) defined by the clockwise flow of major ocean currents 5 .Daily, from 0800 to 1800 h, observers with binoculars on the outer bridge wings of the vessel (circa 11.3 m above sea level) monitored the port and starboard sides of the vessel for Sargassum.When Sargassum was sighted, a rigid inflatable boat (RIB) was deployed from which to conduct sampling (Fig. 1).www.nature.com/scientificreports/Sampling was conducted during daylight hours when the solar disk was not obstructed by clouds and in sea states of Beaufort Force ≤ 3. Sampling was restricted to these conditions to minimize variation in incoming solar radiation and water movement and control for temporal variation in environmental conditions that affect seawater temperature 32 and Sargassum aggregation state 33 .Drifting Sargassum occurs as dispersed fragments and individual thalli or "clumps" 17,21,33,34 typically 0.1-0.5 m across 3 ; lines or "windrows" of Sargassum aggregated parallel with the wind direction and stabilized by Langmuir circulation 2,3 varying in width from < 0.5 m to several m 22,33,34 ; and "mats" or "rafts" of aggregated Sargassum several m in extent up to tens or hundreds of m across 3,17,22,33,34 , which may be associated with windrows 33,34 or occur as discrete aggregations more likely to form in calm conditions 2,3 .Sargassum distribution in the Sargasso Sea is irregular and patchy 3 , large aggregations are uncommon, and clumps and windrows are most frequently observed in groups of 1-5 33 .Because patterns of Sargassum aggregation and distribution observed during the study period were consistent with previous reports, we sampled Sargassum aggregations as they were encountered unless aggregations of multiple size strata were visible from the RIB.If so, to ensure the range of aggregation sizes sampled represented the range of aggregation sizes present, we sampled aggregations ≥ 1 m wide (n = 28) as they were encountered and sampled every fourth aggregation < 1 m wide (n = 22) encountered.Sargassum aggregations were considered independent units if they were at least 5 m apart.
Aggregation length and maximum width (perpendicular to the length axis) were visually estimated to the nearest 0.25 m from an eye height above sea level of 2.0 m.Using a YSI Pro20i temperature probe rigged to an extendable pole (Fig. 1), we measured SST at 0.5-m intervals along a transect from aggregation edge at the aggregation's widest point toward aggregation center.At each horizontal position, we also inserted the probe into the aggregation to measure temperature at 0.2-m increments to a depth of 1 m.After recording temperature measurements, we snorkeled into the aggregation to measure aggregation thickness at each horizontal position using a PVC pole marked at 0.2-m increments.For aggregations < 0.5 m wide, only SST at aggregation center was measured, and thickness was visually estimated to the nearest 0.1 m.For aggregations ≥ 0.5 and < 1 m wide, only temperatures at edge and center were measured, and length and width were recorded to the nearest 0.05 m.For each aggregation, we measured temperatures at 0.2-m depth increments, as described above, in nearby open water at a single point circa 5 m away.For aggregations in which we only measured SST, only SST was measured at the paired open water point.Using an extendable pole (up to 4.9 m) allowed for data collection across the entire range of aggregation sizes encountered while minimizing water disturbance by the RIB.See Supplementary Information for a video of our methods.
Because the greater proportion of aggregations were elongate (length > width), we calculated area for each aggregation as maximum width × length rather than as circular area.We calculated aggregation volume as area × thickness at aggregation center.We classified temperature measurements as inside aggregation or below aggregation if depth of the measurement was ≤ or > aggregation thickness, respectively.

Statistical analyses.
Analyses were performed in R v. 4.0.0 35 using the dplyr package 36 .SST in the center of Sargassum aggregations was compared to nearby open water using a non-parametric paired-samples Wilcoxon test.Linear regression was used to evaluate relationships between aggregation dimensions (volume, area, thickness) and SST differential (aggregation center-open water), with log-transformed (dimension) and (differential + 1) to meet parametric assumptions.Non-parametric Kruskal-Wallis and post-hoc Dunn's tests with Bonferroni correction were used to compare SST differentials among open water, aggregation edge, and aggregation center; temperatures inside aggregations, below aggregations, and open water; and temperatures among 0.2-m depth intervals.Differences in sample size among groups are due to taking fewer temperature measurements for smaller aggregations, as described above.The assumption of normality was checked for all comparisons using Shapiro-Wilk tests.Homogeneity of variance was checked using Levene's tests for comparisons with a categorical predictor variable and by visual inspection of residual plots for comparisons with a continuous predictor variable.When parametric assumptions were not met, a non-parametric alternative was used, as specified above.All statistical tests were two-tailed, and significance was evaluated at an alpha level of 0.05.All means are reported ± standard error (s.e.m.).
Vol:.( 1234567890 22 and the potential benefit of this thermal niche for associated organisms 23 .Our study is the first to evaluate thermal profiles of this habitat in situ and at a scale relevant to conditions experienced by Sargassumassociated fauna.www.nature.com/scientificreports/Light attenuation due to high macrophyte [37][38][39] or phytoplankton biomass 40 results in surface water warming.Sargassum aggregations likely affect water temperature through the same mechanism 22,23 , altering the absorption profile of solar radiation 41,42 .Temperature differentials are likely driven by increased heating of water inside Sargassum aggregations mediated by environmental factors (incoming solar radiation, wind speed, current velocity) and aggregation characteristics like size and density 37,38,40 .Sargassum aggregations range in size from clumps < 0.5 m across to dense mats hundreds of m across 2,3,33,34 .support the prediction that temperature differential will vary with aggregation size and suggest there may be a size threshold above which temperature differential increases.Marmorino et al. 22 reported that Sargassum aggregations > 10 m across (area > 80 m 2 ) appeared 0.1-0.5 °C warmer than adjacent seawater in airborne imagery, and a differential was not detected for smaller aggregations.Our study provides direct field observations supporting solar warming of smaller Sargassum aggregations (area ≤ 15 m 2 ).The range in aggregation size sampled in the current study is consistent with the aggregation characteristics most frequently encountered in the Sargasso Sea 3,33 and in the tropical Atlantic and Caribbean Sea 33,34 .Future studies should quantify thermal characteristics of larger Sargassum aggregations in situ, given continued expansion of Sargassum biomass 10 and the recurrence of Sargassum aggregations > 50 m across in the Caribbean and off West Africa 33,34 .Though the geographic scope of sampling in the current study was limited to an area southeast of Bermuda in the Sargasso Sea, solar warming of the ocean surface occurs in all ocean basins 32 , thus we expect the mechanisms driving warming of Sargassum aggregations are consistent within and among regions.We predict SST differential increases with aggregation size above a threshold sufficient to reduce water movement inside the aggregation (Fig. 6(i)), with the greatest increase in temperature relative to open water at the center of aggregations large enough to minimize water movement and slow heat dispersal (Fig. 6(ii)).
Water movement driven by currents and wind could substantially affect temperature differentials in two ways.First, water movement through aggregations may transport warmed water out and cooler water in.Natural water movement, apparent in our methods video (see Supplementary Information), likely increased water transport and reduced temperature differentials in the current study (Fig. 6(iii-iv)).Second, water movement and wind acting directly on emergent Sargassum alter shapes and densities of Sargassum aggregations 3,22 , thereby affecting thermal stability, which we predict would be greater in large, dense aggregations (Fig. 6(v)) 38,41,43 .Sargassum aggregations have been described as "wave-subduers" 44 .The stabilizing effect of Sargassum may be responsible for reduced water temperature variance at the surface of aggregations (Fig. 5b) compared to SST in open water (Fig. 5a) potentially via reduced air-sea heat flux 32 .
Diel variation in SST 32 may also affect temperature differentials.In the Sargasso Sea, diurnal SST ranges 0.1-2.6 °C during the diurnal thermal cycle, with minimum and maximum temperatures occurring at 0400-0500 and 1500-1600 UT, respectively 45 .Diurnal temperature fluctuation in dense benthic Sargassum canopies lags hours behind water outside the canopy, with daytime temperatures maintained inside canopies at night 46,47 .This was recorded in a Sargassum forest 50 m wide with a canopy height of 3 m 47 , which is comparable to the dimensions of the largest holopelagic Sargassum aggregations reported in the eastern Caribbean (50-100s of m wide and up to 7 m thick) 34 .Thus, it is unlikely that temperatures in typical Sargassum aggregations (up to several m wide 33,34 and ≤ 0.5 m thick 34 ), including those sampled in the current study, lag behind open water.Diurnal variation in temperature differentials should be evaluated across depth in larger Sargassum aggregations.In dense macrophyte canopies, daytime temperatures are cooler and more stable below the surface, while at night, surface water cooling combined with reduced water movement in dense structure can cause temperature to be higher with depth 43 .Greater thermal stability in large Sargassum aggregations should result in slower warming after sunrise than in open water and daytime temperature differentials that decrease with depth, as well as longer retention of daytime heat after sunset and nighttime temperature differentials that increase with depth.
Biological significance of Sargassum microclimates.Temperature affects the rate at which biochemical reactions occur, and metabolic rates scale with temperature 30,48,49 .Across diverse taxa and habitats, most rates of biological activity increase exponentially with temperature over the temperature range in which an organism is normally active 30 .Small marine ectotherms, which dominate the holopelagic Sargassum community 4,20 , have body temperatures that track ambient water temperature 29 .Thermal microclimates directly impact ectotherm performance and influence fitness 29,30 .Fish and other mobile ectotherms are able to regulate body temperature by moving among thermal microhabitats 28,31 .Behavioral thermoregulation can reduce metabolic costs, improve performance, and increase fitness 28,50 .Similar to other aquatic macroalgae [24][25][26] , holopelagic Sargassum should provide favorable thermal microclimates for associated fauna 23 .Due to a lack of knowledge related to species thermal tolerance ranges for those associated with holopelagic Sargassum, it is difficult to assess how species' thermal tolerance coincides with the temperature differentials we documented in Sargassum aggregations, or how climate change may affect the viability of thermal microclimates.Rising atmospheric CO 2 concentrations are resulting in ocean warming and increased acidity [51][52][53] .The complexity of factors mediating thermal thresholds and fitness outcomes [54][55][56] , as well as limited knowledge of species' thermal tolerance ranges, genetic adaptation potentials, and acclimation abilities 57,58 , make predicting impacts of warming on individual species difficult.Under future climate conditions, temperatures in holopelagic Sargassum aggregations would exceed upper thermal limits of associated species before temperatures in open water.If selecting Sargassum habitat becomes maladaptive despite other fitness advantages 58 , Sargassum aggregations may become ecological traps.Tropical species are probably more vulnerable to warming than temperate species because most have narrower thermal tolerance breadths and limited acclimation ability 56,[58][59][60] .Thus, risk must be evaluated across latitudes in which holopelagic Sargassum is found.Predicting the effects of climate change on thermal microhabitats of holopelagic Sargassum aggregations is outside of the scope of our study.However, our theoretical model of the relationship between temperature differential and aggregation size (Fig. 6) provides a conceptual framework for exploring potential impacts of ocean warming on thermal microclimates, in addition to evaluating thermal profiles across broader geographic and temporal continua.
Because temperature differentials have important biological implications for Sargassum ecosystems, future studies should evaluate temporal and geographic variation in thermal microhabitats in situ, particularly of larger aggregations.Our study provides a foundation for investigating the importance of thermal microclimates for holopelagic Sargassum-associated fauna, which has important implications for anticipating impacts of climate change on the suitability of this thermal habitat for endemic, endangered, and commercially important species that depend on Sargassum ecosystems.

Figure 1 .
Figure 1.Researchers measure water temperatures in a Sargassum aggregation using a YSI Pro20i temperature probe rigged to an extendable pole (left), using black bands on the PVC pipe as a depth reference (right).Photo credits: Tavish Campbell/Greenpeace (left); N.C.(right).

Figure 6 .
Figure 6.Predicted effect of Sargassum aggregation size on the surface temperature differential between Sargassum aggregation center and open water.

Table 1 .
Post-hoc pairwise comparisons of SST differentials (°C) among open water, Sargassum aggregation edge, and aggregation center using Dunn's test with Bonferroni correction.

Table 2 .
Post-hoc pairwise comparisons of temperature (°C) in open water, inside Sargassum aggregations, and below aggregations using Dunn's test with Bonferroni correction.