High diatom species turnover in a Baltic Sea rock pool metacommunity

Different metacommunity perspectives have been developed to describe the relationship between environmental and spatial factors and their relative roles for local communities. However, only little is known about temporal variation in metacommunities and their underlying drivers. We examined temporal variation in the relative roles of environmental and spatial factors for diatom community composition among brackish-watered rock pools on the Baltic Sea coast over a 3-month period. We used a combination of direct ordination, variation partition, and Mantel tests to investigate the metacommunity patterns. The studied communities housed a mixture of freshwater, brackish, and marine species, with a decreasing share of salinity tolerant species along both temporal and spatial gradients. The community composition was explained by both environmental and spatial variables (especially conductivity and distance from the sea) in each month; the joint effect of these factors was consistently larger than the pure effects of either variable group. Community similarity was related to both environmental and spatial distance between the pools even when the other variable group was controlled for. The relative influence of environmental factors increased with time, accounting for the largest share of the variation in species composition and distance decay of similarity in July. Metacommunity organization in the studied rock pools was probably largely explained by a combination of species sorting and mass effect given the small spatial study scale. The found strong distance decay of community similarity indicates spatially highly heterogeneous diatom communities mainly driven by temporally varying conductivity gradient at the marine-freshwater transition zone.


Introduction
Microbial communities contribute to a variety of biogeochemical processes vital for the functioning of ecosystems. Variation among these communities thus potentially affects the whole ecosystem, and better knowledge of the key factors determining microbial communities is urgently needed (Van der Gucht et al. 2007). Due to their high diversity, relatively well-defined habitat preferences, short generation times, and generally-yet not invariably-geographically wide distribution, microbial algae such as diatoms are useful model taxa in testing different ecological theories over short timescales (Smucker and Vis 2011). So far, studies of community dynamics and especially spatiotemporal variation in species composition for microbial eukaryotic organisms are still rare compared with those of macroorganisms (Martiny et al. 2006;Meier and Soininen 2014). A metacommunity consists of local communities occupying discrete habitat patches connected by potentially interacting species through dispersal (Hanski and Gilpin 1991;Leibold et al. 2004). Through the application of metacommunity theory, community-level patterns in species abundance, distribution and biotic interactions are tightly coupled and explained against the theoretical background of spatial and environmental processes on a scale of metacommunities (Leibold et al. 2004;Logue and Lindström 2008).
The current metacommunity theory comprises four conceptual archetypes, each reflecting variable roles of environmental and spatial factors for metacommunity dynamics. We note here that these archetypes form recognizable dynamics in Communicated by B. Beszteri Electronic supplementary material The online version of this article (https://doi.org/10.1007/s12526-019-01016-z) contains supplementary material, which is available to authorized users. a continuum of metacommunity parameter space rather than discrete paradigms (Logue et al. 2011). Neutral theory (Hubbell 2001) stresses the importance of spatial processes in regulating metacommunities (Chave 2004;Cottenie 2005). Species are thought to be equal in their environmental demands, and random stochastic events related to dispersal or births and deaths are the only determinants of species distribution. According to species sorting, each species has specific tolerance toward environmental conditions expressed by ecological niches, separating species to environmentally varying habitat patches (Leibold et al. 2004). Species dispersal rate is low enough not to hold the dominating role in community regulation (Logue et al. 2011) and acts in the shadow of more powerful environmental processes, enabling species to reach suitable sites while not resulting in independent spatial patterns in a metacommunity (Cottenie and De Meester 2004;Cottenie 2005). The mass effect perspective (Shmida and Wilson 1985) involves both species dispersal and environmental factors as being responsible for a community organization (Leibold et al. 2004). As environmental preferences divide species to their local habitats, dispersal between patches maintains the viability of communities by enabling temporary species occurrence outside their preferred distribution ranges (Pulliam 2000;Logue et al. 2011). The fourth archetype, that of patch dynamics (Leibold et al. 2004), emphasizes dispersal limitation as in control of local species diversity among identical patches, either occupied or unoccupied according to local extinction-colonization dynamics.
Although many studies have now addressed the metacommunity structure in different types of ecosystems, temporal variation in metacommunity structure and the underlying factors are still poorly studied (but see, e.g., Azeria and Kolasa 2008;Vanschoenwinkel et al. 2010;Langenheder et al. 2012) and occasionally weakly connected to variation in species composition (Padial et al. 2014). Specifically, the dependence of spatial variables such as pool connectivity on temporal variability in local precipitation patterns and consequent changes in dispersal rates between sites may alter species composition and turnover rate especially in coastal communities (Langenheder et al. 2012). Ecosystems are expected to face rapid changes in the current era of global warming, highlighting the need for metacommunity-level information about temporal and spatial variation in and their effects on species composition and turnover rate not limited by snapshot sampling. Two metrics-halving distance and initial similarity-developed to describe patterns of distance decay between communities give cues about species turnover more effectively than traditional correlation coefficients. Limited dispersal or high environmental heterogeneity between sites results in faster decrease in community similarity along environmental and spatial distance and higher turnover rate, characterized by low initial similarity and short halving distance Wetzel et al. 2012).
For investigating metacommunity dynamics, granitic rock pools provide an ideal model system. These are geologically old formations occurring as ubiquitously scattered depressions on rocky outcrops typically surrounded by sea (Jocque et al. 2006). Despite their small size, these physically and biologically diverse waterbodies form a disproportionately large share of regional-scale aquatic biodiversity (Blaustein and Schwartz 2001;Brendonck et al. 2016). As hierarchically nested and closely located depressions, each rock pool serves as a small-scale community connected by species dispersal within a metacommunity (Blaustein and Schwartz 2001). During the last decade, the potential of rock pools as ecological model systems has been shown (Vanschoenwinkel et al. 2007). Representative sampling and a quantification of community dynamics and dispersal processes is more effortless and more easily replicated in time than in larger and more complex aquatic systems and study results are also usually comparable and easily interpolated to different spatial scales (Srivastava et al. 2004;De Meester et al. 2005).
This study aimed to examine temporal variation in the relative roles of environmental and spatial factors for a metacommunity composition. Specifically, this study aimed to answer to the following research questions: First, which environmental and spatial variables affect the community composition in rock pools over time? Second, how do rock pool diatom community similarities vary along the environmental and spatial gradients over time? Third, which metacommunity archetype or a combination of archetypes best explain the observed variation among the rock pool diatom communities?

Study area
The study area consisted of 30 granitic rock pools on a rocky outcrop (ca. 1200 m 2 ) on the western island of Pihlajasaari (60°8′ 10″ N, 24°54′ 46″ E), located approximately 2 km south of Helsinki, Finland, and surrounded by the northern part of the Baltic Sea, the Gulf of Finland (salinity ca. 5‰) (Fig. 1). Tidal influence in the area is virtually absent, and changes of sea level are mainly caused by wind and differences in air pressure. The majority of the studied pools were located at or close to the sea surface level (< 5 m a.s.l.). Pools of varying size and location were chosen for the study and named by running numbers from 1 to 35. As pools 31-34 dried up until the first sampling in May and pool 24 was merged with pool 22 by a rainfall-induced overflow from June onwards, the data consist of pools 1-23, 25-30, and 35 (Online Resource 1).
The majority of the pools were rainfall-fed, but the pools closest to the sea were directly influenced by brackish seawater by wind-caused waves and indirectly by salty sprays from the sea (Fig. 2). Despite merging of the pools 22 and 24 after May, the pools were isolated from each other and not permanently interconnected by watercourses. Due to an impermeable granitic pool base, groundwater influence was virtually absent, as was surface outflow from the surroundings in the absence of any nearby freshwater sources. The pools were not shaded by the sparse vegetation, making solar radiation conditions comparable. Nutrient enrichment was likely mostly of a faunal origin, caused by excrements and decaying remains of aquatic insects and their larvae and nutrient inputs from, e.g., tadpoles, bird droppings, and other animal feces (Methratta 2004;Brendonck et al. 2016).

Field sampling
The field data were collected once a month on 17 May, 22 June, and 22 July 2016. The environmental data consist of water physicochemical characteristics (pH, conductivity, temperature, total P, and total N) measured in the field with YSI field meter and pool size measurements (max length, width, and depth) measured with a meter stick to the nearest centimeter. Pool area was calculated by multiplying pool length by pool width. For two very irregularly shaped pools (pools 10 and 12), means of multiple measurements of both pool width and length were used to better reflect the real pool surface area. A 0.5-L water sample was collected from each pool in a plastic container and preserved at 4°C until the determination of total P following standard SFS- EN ISO 6878 (2004) and total N following standard SFS-EN ISO 11905-1 (1998) in a laboratory. The determination of total N was occasionally disturbed by high organic matter content, which, like nitrogen, also absorbs radiation on the UV wavelength of 220 nm measured by the UV/VIS spectrometer. Due to such a high disturbance in 21 samples (15 in June and 6 in July) and two results exceeding the determination limit in May, all total N values were omitted from the data.
The spatial data consisted of pool X and Y coordinates and mean isolation. Due to the very small spatial scale of the study and consequent difficulties in determining the exact pool coordinates by GPS or from aerial photographs, the coordinates were estimated from a drawn grid map of the study area showing the relative location of the pools to each other and to the seashore. X coordinates were based on perpendicular pool distance from the shore, while Y coordinates described horizontal pool distance from the map origin (the south end) parallel to the shoreline. The relative isolation was calculated for each pool as a mean Euclidean distance, i.e., the sum of distances to five closest pools divided by five (Vanschoenwinkel et al. 2007).
The benthic diatoms were sampled following EN 13946 standard (2003). Ten epilithic subsamples of c. 25-cm 2 area were collected by a toothbrush from each pool and combined as a single composite sample in a plastic test tube in the field. The subsamples were collected from different sides of each pool from the pool bottom. From the two steepest-walled pools (pools 10 and 14) the pool bottom was unreachable and only the walls were sampled deep enough to ensure the frustules were permanently submerged. Between each sampling, the toothbrush was rinsed in pool water to remove any  (2018) attached cells as this effectively reduces the possibility of contamination and taxonomical bias between the samples (Kelly et al. 1998). The samples were stored at 4°C for 24 h until microscopical analyses. The samples were treated with 30% H 2 0 2 to remove organic material and mounted on slides with Naphrax. A total of 500 frustules per slide were counted and identified to the lowest taxonomic level possible (mostly species level) with a light microscope with 1000-fold magnification following Krammer andLange-Bertalot (1986-1991). From the four most sparsely slides celled, less than 500 frustules could be counted (pool 21 (n = 476) in June and pools 9 (n = 414), 11 (n = 385), and 20 (n = 156) in July). The total, absolute cell abundances were converted to relative percentages to make the cell numbers comparable with each other. The identified species were classified according to their preferences for water salinity into marine, brackish, and freshwater species, as well as species favoring both marine and freshwater environments (excluding 8 specimens identified to genus level and 1 Nitzschia specimen with an unknown salinity preference) mainly after Cantonati et al. (2017).

Statistical analyses
Prior to statistical analyses, normality of the explanatory variables was tested with the Shapiro-Wilk's normality test (Shapiro and Wilk 1965) and all the non-normally distributed explanatory variables (except pH) were log 10 -transformed to reduce their skewness. Spearman's rank correlation coefficients were calculated to detect statistical dependence between the variables. In June, conductivity and Y coordinates were strongly correlated (r s = − 0.79) and hence Y coordinates were excluded from the June's statistical models (but included in the combined RDA model as the correlation between these variables was < |0.7| in the combined dataset). Temporal fluctuations in pool physicochemical characteristics were compared with a short-term climate data from Kumpula meteorological station, Helsinki, including daily precipitation (mm), mean air temperature (°C), and mean wind speed (m s −1 ) covering a 7-day period before each sampling (FMI 2017; Online Resource 2).
Spatial autocorrelation was calculated for the environmental variables by Moran's I (Moran 1950). Water salinity is an important factor regulating pool biota especially in coastal areas under significant seawater influence, and even minor changes along a pool conductivity gradient may have a major impact on communities (Vanormelingen et al. 2008). Thus, the distance to the sea (X) was considered more important gradient than the horizontal distance along the shoreline (Y). Since Moran's I automatically calculates coefficients for the longest geographical gradient (i.e., Y coordinates), the value of each Y coordinate was set to zero while keeping X coordinates unmodified. The data were divided into five distance classes at 2.5-m intervals. The p values of the correlation coefficients were Bonferroni-corrected with a formula α′ = α/k, where k is the number of distance classes and α the significance level of 0.05 (Legendre and Legendre 2012). Every correlation coefficient below 0.05/5 ≈ 0.01 was considered statistically significant.
The most influential explanatory variables for diatom community composition were examined with redundancy analyses (RDA; ter Braak 1994) with a Hellinger-transformed species data. An RDA is a direct gradient analysis in which species observations are directly related to a set of explanatory variables. An RDA was conducted separately with each monthly dataset; in addition, a fourth RDA was run for the combined dataset from May to July. The statistical Variation partitioning (Borcard et al. 1992) was conducted for each monthly dataset and for the combined data to share the explained variation in community composition to fractions explained by pure spatial, pure environmental, and pure temporal factors (only for the combined data), as well as by their joint effect. The statistical significance of the fractions was assessed by RDAs run separately for each variable set using the F test with 999 permutations (Legendre 1993).
Mantel tests (Mantel 1967) were performed to reveal the roles of spatial and environmental distance on community similarity (expressed using the Bray-Curtis dissimilarity index; Bray and Curtis 1957) in each month. Three distance matrices were created: a biotic dissimilarity matrix using species relative abundances, an environmental Euclidean distance matrix using water physicochemical characteristics and pool morphometrics, and a spatial Euclidean distance matrix using pool X and Y coordinates. Halving distances of community similarities were calculated with a formula (β − α)/2β, where β is the regression coefficient and α the intercept, and initial similarities were determined at a 1-m distance ) which was considered sufficiently long distance given the small study scale and short inter-pool distances. A partial Mantel test (Legendre and Legendre 2012) was conducted for each month with the same explanatory variables as in the Mantel tests. The partial Mantel test compares all three distance matrices simultaneously and divides the variation in distances among sites explained by pure environmental and pure spatial variables, while controlling for the effects of the other variable set. The significance of correlations for the Mantel and partial Mantel tests was assessed with 9999 permutations.

Results
High spatial and temporal environmental heterogeneity was detected among the rock pools throughout the study (Table 1 and Online Resource 3). In general, the studied pools were small (< 4.8 m 2 ), shallow (≤ 0.4 m deep), eutrophic (mean total P > 80 μg L −1 ), and alkaline (mean pH > 7.5). As it should be noted that the calculated pool area does unavoidably deviate from real pool surface area to some extent, the variable well reflected the relative size difference between the pools. Nutrient concentrations, conductivity, temperature, and pH were the highest in May; extremely high maximum values were recorded especially for conductivity (16,741 μS cm −1 ) and total P (1448.3 μg L −1 ). According to a principal component analysis (PCA), the variation in environmental variables was the lowest in July (the results not shown). The majority of the pools were located farther than 5 m from the seashore and not farther than 10 m from the five closest pools (Online Resource 4). Significant spatial autocorrelation (Bonferronicorrected p < 0.01) was detected for water depth in June (Online Resource 5a) and for total P, temperature, and conductivity in July (Online Resource 5b-d).
Monthly fluctuations in pool physicochemistry clearly reflected temporal variations in the local microclimate especially in May and June when compared with the short-term climate data from Kumpula meteorological station, Helsinki (Online Resource 2, 3). Higher air temperatures and low precipitation were followed by shallow and warm water with high conductivity and nutrient concentrations in May, while heavy rainfalls, higher wind speed, and lower air temperature decreased water temperatures and nutrient levels and increased pool depth in June. The studied pools followed a clear gradient in water conductivity along both X (r May = − 0.45, r June = −0.43, r July = − 0.65, p < 0.05) and Y coordinates (r May = − 0.62, r June = − 0.79, r July = − 0.53, p < 0.01) in each month, as pools with the least conductive water were located in the more sheltered northeastern end of the study area (Online Resource 6a-c).
A total of 179 species belonging to 62 genera were detected from May to July. Of these genera, the majority (73%) were represented by a maximum of two species. The most diverse genera were Nitzschia (29 spp.) and Navicula (26 spp.). Onethird (33%) of all detected species were present in a single month only and the majority (62%) of these species occurred as a single cell only (Online Resource 7). Half of the detected species occurred every month. In May and June, a few species occupied every pool, while in July, none of the species occupied all pools (Online Resource 8). The pools housed a mixture of freshwater (57% of all identified species in May-July), brackish (21%) and marine (14%) species, and a few species (8%) typical of both marine and freshwater habitats ( Table 2). The majority (67% in May up to 87% in July) of the pool communities were dominated by freshwater species in each month, whereas the rest of the pools (3 pools in May and 1 pool in June-July) were dominated by brackish species and one pool in May by marine species (Online Resource 9).
The relative proportion of marine, brackish and freshwater species among the pools varied along a spatial gradient (i.e., by either X or Y coordinates) in each month. The proportion of marine species among the pools decreased significantly (r s = − 0.44, p < 0.05) with Y coordinates in June (r s = − 0.3, p < 0.05) and with X coordinates in July (r s = − 0.53, p < 0.01), while that of brackish species decreased significantly along both X (r s = − 0.44, p < 0.05) and Y coordinates in May (r s = − 0.52, p < 0.01) (Online Resource 6a-c). The per-pool proportion of freshwater species increased significantly along both X (r s = 0.49, p < 0.01) and Y coordinates (r s = 0.55, p < 0.01) in May and along X coordinates in July (r s = 0.48, p < 0.01).
The per-pool proportion of brackish species decreased also significantly (r s = − 0.41, p < 0.001) along a temporal gradient while the proportion of species typical of both marine and freshwater habitats increased significantly (r s = − 0.43, p < 0.001) from May to July (Online Resource 6d). The total monthly proportion of marine species increased and that of brackish species decreased from May to July (Table 2). The total proportion of freshwater species was rather comparable between the months, peaking in June. Fig. 3 Results for the monthly redundancy analyses and the one with combined data: a May; b June; c July; d all months. Shown are the location of the studied rock pools (numbered from May to July, symbolled in the combined model) in relation to the ordination axes and explanatory variables Fig. 4 Results of the monthly Mantel tests. Shown is the change in the Bray-Curtis similarity index along environmental (left) and spatial (right) distance, together with a correlation coefficient for the relationship between species and both spatial and environmental variables, respectively. The significance of correlations is based on the p value: *p < 0.05, **p < 0.01, ***p < 0.001 The factors affecting community composition also varied in time, but conductivity was consistently one of the key factors (Table 3). In May, the community composition was significantly explained by conductivity (6.6% of the explained variance, p < 0.001) and X coordinates (3.4%, p < 0.01) in the RDA; the first two ordination axes explained together 28.5% of the variation in community composition (Fig. 3a). In June, conductivity (9.1%, p < 0.001), X coordinates (3.6%, p < 0.05), and pH (3.2%, p < 0.05) were significant, and the first two axes explained 32.5% of the variation (Fig. 3b). In July, conductivity (10.2%, p < 0.001), temperature (3.8%, p < 0.05), total P (3.7%, p < 0.05), and isolation (3.7%, p < 0.05) were significant, while the first and the second ordination axes explained 35.6% of the variation (Fig. 3c). In the combined RDA, several variables were significant; the variation was best explained by conductivity (7.6%, p < 0.001), X coordinates (2.6%, p < 0.001), isolation (2.2%, p < 0.001), and month as a categorical variable (2.2%, p < 0.01) (Fig. 3d). The first two axes explained 25.7% of the variation.
The community similarity decreased significantly with an increase in both spatial and environmental distance in each month even when the other variable group was controlled for. The correlation between similarities and spatial distance was nearly equal (r May = − 0.27, r June = − 0.28, r July = − 0.28) in ). The unexplained proportion of variation in community composition is expressed as a residual variation (R). The significance of pure environmental and pure spatial components is based on the p value of the F test for a redundancy analysis for the two variable groups: *p < 0.05, **p < 0.01, ***p < 0.001 Fig. 6 The results of variation partitioning for the combined monthly data. Shown are the proportion and statistical significance of the explained variation in species community composition by pure environmental (Environment), pure spatial (Space), and pure temporal (Month) variables, together with the joint effects of these three variable sets. The unexplained proportion of the variation is expressed as a residual variation. The significance of pure environmental, pure spatial, and pure temporal components is based on the p value of the F test for a redundancy analysis for each of the three variable groups: *p < 0.05, **p < 0.01, ***p < 0.001 each month, while decrease in similarity along environmental distance increased from May (r = − 0.18) to July (r = − 0.32) ( Fig. 4; Table 4). In May and June, distance decay was stronger for spatial distance than environmental distance, whereas in July, environmental distance was slightly more strongly related with community similarities. The halving distance of community composition was extremely short and decreased from May (108.7 m) to July (58.9 m). The initial similarity was nearly equal in May and June (r = 0.37), decreasing slightly toward July (r = 0.32). The correlation between environmental and spatial distance was rather low in each month (r ≤ 0.43; Table 4). The joint effect of environmental and spatial factors consistently explained the largest share (8.3-13.1%) of the monthly variation in community composition (Fig. 5). The variation explained by pure environmental and pure spatial factors was nearly similar in each month. In May, pure spatial variables accounted for a slightly higher proportion (8.2%) of the explained variation than pure environmental variables (5.5%), while in June and July, pure environmental variables (6.5% and 6.9%) explained more of the variation than pure spatial ones (6.3% and 2.8%). Temporal variation in the community composition was significant, accounting for 3.7% of the variation over the whole 3-month period (Fig. 6). The largest share of the variation was jointly explained by environmental and spatial variables (9.9%), while their separate effects were nearly equally weak (6.8% and 6.5%, respectively).

Discussion
Our study revealed high diatom species turnover in a rock pool metacommunity linked to both environmental heterogeneity and pool spatial distance from the sea through variation Table 2 Monthly statistics for the per-pool proportion of diatom species belonging to the four ecological groups according to species-specific salinity tolerances (Taylor et al. 2007;Antón-Garrido et al. 2013;Cantonati et al. 2017 Table 3 Results for the monthly and combined redundancy analyses. Shown is the proportion of explained variance by each statistically significant explanatory variable of the total explained variance in each month, together with their correlation with the species matrix along the first two ordination axes. The significance of the variables is based on the p value of the F test: *p < 0.05, **p < 0.01, ***p < 0.001 in temporal gradient. Given the small spatial scale of the study, the metacommunity was likely mainly structured by mass effect and species sorting, the latter of which was dominated by a conductivity gradient at the marine-freshwater transition zone. The majority of the pool diatoms were freshwater species both temporally and spatially. However, at the ecotone between marine and freshwater habitat, freshwater-tolerant marine species were introduced into the mainly freshwater pools from the nearby sea (Schröder et al. 2015), shifting species composition of the studied communities from freshwater toward more salinity tolerant species along the conductivity gradient.
Although larger-sized species are rarely encountered and hence often underrepresented in microbial samples (Lowe and Pan 1996;Snoeijs et al. 2002) or the presence of some species is sporadic enough to increase the probability for an omission error due to sampling effect (Bottin et al. 2014), the highly uneven local species abundance observed here is a general pattern very typical of diatom communities (Potapova and Charles 2002;Heino et al. 2010) and usually a consequence of the dominance of few species capable of permanently residing in all sites even under frequent environmental disturbance (Dethier 1984). However, due to generally high local species diversity, relative compositional differences between diatom communities are usually captured at a sufficient accuracy with relatively few samples (Eloranta 2004), independent of potential omittance of the most rarely encountered species in a given habitat.

Drivers of diatom community composition
Water conductivity explained the largest share of variance also in species composition in each month. Several studies support the significance of either conductivity or salinity for both freshwater diatom and coastal rock pool communities (Ganning 1971;Metaxas and Scheibling 1993;Virtanen and Soininen 2012). Either the degree of pool isolation or the distance from the sea was significant for the community composition in each month, with importance of the latter decreasing over time. This might be due to the monthly varying weather conditions. Even minor temporal variations in microclimate may have strong effects on local environmental conditions over a given disturbance regime, directly resulting in relative differences in water depth, wave exposure, and other physicochemical conditions among pools with varying hydroperiod length (Ganning 1971;Jocque et al. 2010). Thus, even under seemingly similar environmental conditions, neighboring habitat patches may well differ in their microscale physicochemistry and hence in community composition (Metaxas and Scheibling 1993;Jenkins and Buikema 1998;Brendonck et al. 2016). Seasonal species succession and temporal changes in water physicochemical variables such as total P, temperature, and pH were clearly important for the studied communities, expressed by the significance of the categorical "month" variable in the combined RDA model. Especially in June, stronger winds may have lowered the influence of the transverse gradient reflecting spatial variation in water conductivity, extending the influence of seawater from the more exposed southwest further toward the sheltered northeast and distributing brackish species more evenly between the pools. Contrary to the two other months, the distribution of both marine, brackish, and freshwater species among the pools was unaffected by the pool spatial location in June.
We also note that banks of sedimented propagules may have acted as a resisting force for local extinctions often targeted to small aquatic habitats by maintaining viable local populations, further affecting the species composition in the studied communities . The potential influence of such a seed bank on the evolution of a given community is nevertheless likely better discerned in longerterm studies covering several years and growing seasons (Ellegaard and Ribeiro 2018), whereas our study was restricted to a single growing season at a timescale of only 3 months.

Distance decay of community similarity
Although typically mostly under environmental control, significant distance decay has been documented for diatom assemblages over varying environmental conditions, at multiple temporal and spatial scales, in lentic and lotic waters and along a salinity gradient from freshwater to marine habitats (Soininen 2007;Vanschoenwinkel et al. 2007; Astorga et al. Table 4 Results for the monthly Mantel and partial Mantel tests. Shown are monthly correlation coefficients for pure environmental (Env) and pure spatial variables (Spat), and the corresponding coefficients for these variable groups when the other variable set was controlled for (Env-Spat and Spat-Env, respectively). For a change in dissimilarity index along the spatial distance, the halving distance of community similarity (HD, in meters) is also shown. The initial similarity (IS) refers to the dissimilarity index value at a spatial distance of 1 m. Shown is also the coefficient of correlation between environmental and spatial distance (distance r) in each month. The statistical significance of the variables is based on the p value: *p < 0.05, **p < 0.01, ***p < 0.001 2012; Bottin et al. 2014;Virta and Soininen 2017). Both environmental and spatial distances were related with community similarities in each month also here even when the other variable group was controlled for. In May and June, spatial distance decay was stronger than for environmental distance, while in July, community similarities decayed more strongly with the environmental distance. The increase of distance decay with time was more pronounced along environmental distance than along geographical distance, reflecting growing relative importance of the environment for the communities during the study period. Very short and steadily decreasing halving distances were observed from May to July-such values of 60-110 m are extremely short compared with a global average of over 630 km across all biotic communities ) and much shorter than has been reported earlier for freshwater diatoms (Teittinen et al. 2016). Such high species turnover indicates that the salinity gradient acted as a strong filter on the diatom species at this marinefreshwater transition zone. The proportion of salinity-tolerant species varied the most in May: while the share of marine species exceeded 20% in few pools and the spatial dominance of brackish species was the highest, the overall proportion of marine species was the lowest and that of brackish species the highest of the 3 months. The dominance of spatial factors and consequent higher variation in species ecological composition in May could be attributed to the rather calm and dry weather lowering water-and windaided dispersal rates between the pools. Wind speed and direction affect aerial carrying capacity of passive dispersers especially on coastal areas, while additional rains enhance wateraided dispersal (Smith 1973;Vanschoenwinkel et al. 2008a). Although random and sudden overflows may act as a strong dispersal agent allowing species to disperse among temporarily connected pools (Jocque et al. 2007;Brendonck et al. 2010), very short duration of these connections may limit significant dispersal. It has been suggested that significant changes in the community composition of passive dispersers occur only beyond 60 m from the nearest pool independent of a study scale (Vanschoenwinkel et al. 2007). Here, the shortest observed halving distance was little less than 60 m, indicating that compositional changes do happen even beneath this distance limit.
Given the small study extent, dispersal limitation was probably not a major factor contributing to the species composition in the studied communities. In June and July, higher precipitation and wind speed, perhaps accompanied with maturation of faunal dispersal agents such as aquatic invertebrates observed in the field, have likely promoted passive dispersal (Vanschoenwinkel et al. 2008a;Nabout et al. 2009;Castillo-Escrivà et al. 2017) and suppressed the importance of dispersal barriers, increasing the influence of water physicochemistry on the communities especially in July. Aerial carrying capacity has been documented to decline as far as 10 m from a source community (Vanschoenwinkel et al. 2008a)-a distance which the majority of the studied pools were located within. Contrary to May, brackish species were most evenly distributed among the communities in June and dominated only one pool, likely reflecting elevated dispersal rates. Despite the apparent isolation of the pools, either the among-pool distances may have been short enough or the weather conditions suitable for efficient passive dispersal to take place (Teittinen and Soininen 2015). Late appearance of ecologically more specialized species such as Bacillaria paxillifera, Cocconeis scutellum, Epithemia sorex, and E. turgida (Van Dam et al. 1994;Cantonati et al. 2017) likely further explains the growing importance of environmental factors such as salinity, nutrients, and water pH toward the end of the hydroperiod from June onwards (Potapova and Charles 2002;Vanschoenwinkel et al. 2007;Pandit et al. 2009).

Rock pool metacommunity organization
The significance of both pure environmental and pure spatial components at such a small study scale points toward mass effect and species sorting as the leading metacommunity types for the studied communities (Cottenie 2005;Soininen and Weckström 2009;Bottin et al. 2016). The highly heterogeneous environmental conditions likely triggered by random disturbance events and a variable hydroperiod gradient serve as a ground for a strong abiotic environmental control by setting demands for a specialized biota (Therriault and Kolasa 2001;O'Neill 2016), while the minor influence of pool isolation and consequently effective passive dispersal over the relative short distances maintain viable communities even in the pools with less favorable ecological conditions (Van der Gucht et al. 2007;Vanschoenwinkel et al. 2008b;Langenheder et al. 2012). Here, an increasing disturbance frequency with time possibly restricted competitive exclusions by maintaining several ecological niches and supporting multiple species with various environmental adaptations along the conductivity gradient . Especially the proportion of brackish species with a tolerance toward weakly saline water decreased significantly with time, while the proportion of species favoring both marine and freshwater habitats increased significantly toward July. Likely accompanied with a growing supply of these latter, more generalist species (e.g., Fragilaria capucina) from the surrounding pools, this may have resulted in less similar communities and in higher species turnover (i.e., shorter halving distance) between the pools with time (Soininen and Heino 2007;Vanschoenwinkel et al. 2013). However, given the relatively high proportion of unexplained variance in community composition in each month, we cannot completely rule out possible effects of neutrality and patch dynamics on the pool metacommunity; rather, our results only suggest mass effect and species sorting as being the archetypes admittedly responsible of the species composition, regardless of the role of the other two in the studied communities.
We also note that, although we included the most important environmental variables in the analyses, omission of some unmeasured spatially structured environmental variables may have artificially increased the relative importance of pure spatial variables in explaining the pool communities (Liu et al. 2016;Soininen 2016).

Conclusions
The studied rock pools were environmentally heterogeneous, with high spatial and temporal variation in water physicochemical characteristics. Together with dispersal processes, such high environmental heterogeneity also supported high species turnover in the metacommunity. The results give support for mass effect and species sorting as the possible leading metacommunity types for the studied rock pool communities, along with potential yet likely minor contribution of patch dynamics and neutrality. The highly variable and disturbed physicochemical environment and efficient passive dispersal over the relatively short amongpool distances, assisted with suitable weather conditions, likely serve as an additional species supply from the neighboring pools. These results also highlight the high potential of rock pools as a model system for metacommunity studies along steep environmental gradients over time.