Effects of environmental filtering on taxonomic and functional diversity patterns : When spiders and plants provide complementary information to water level management in the Seine estuary

doi:10.21203/rs.3.rs-4598956/v1

Download PDF

Research Article

Effects of environmental filtering on taxonomic and functional diversity patterns : When spiders and plants provide complementary information to water level management in the Seine estuary

https://doi.org/10.21203/rs.3.rs-4598956/v1

This work is licensed under a CC BY 4.0 License

Version 1

posted

You are reading this latest preprint version

In the Seine estuary in northern France, many artificial structures limit the effect of the tide on associated alluvial zones. Consequently, this affects natural environmental filtering mechanisms linked to tidal regimes and water salinity, which directly influences the structure of organism assemblies in adjacent ecosystems. Here, we propose to study the influence of these filters' modifications on spiders and plants, two compartments recognized as complementary in terms of bioindication. However, this complementarity has only been studied to a limited extent and rarely in estuarine environments. To highlight this complementarity, we studied the taxonomic and functional patterns obtained across different topographical levels of two sites with contrasting water level managements. Moreover, particular attention was paid to the influence of the weight attributed to rare species (Q order) to shed light on processes affecting species dominance dynamics between taxa. Overall, spider communities appeared less influenced by environmental filtering than plants on both study sites, with taxonomic and functional diversity both demonstrating a low sensitivity to salinity. Spider community assemblies also demonstrated compositional shifts across study sites, mainly driven by changes in abundance and dominance. In contrast, plant communities appeared more sensitive to environmental constraints and water level management, with responses in terms of composition and species turnover rather than unbalanced abundance, suggesting responses at distinct spatial scales between plant and spider communities.

estuarine ecosystem

spiders

plants

taxonomic diversity

functional diversity

rare species weight

Estuaries are transitional ecosystems between rivers and adjacent seas with high functional and conservation values worldwide (Pétillon et al., 2023). As a result of their geographical positioning, they have unique properties including a longitudinal salinity gradient due to the mix of fresh- and saltwater, as well as a lateral salinity and flooding gradient due to tidal and seasonal cycles (Desender & Maelfait, 1999). Because of the abiotic conditions’ high daily and seasonal variability, estuaries cannot be defined as stable ecosystems but rather as an interconnected mosaic of changing habitats (Davidson, 1991). Despite this diversity of habitats, estuaries are characterized by low species richness due to the limited number of species adapted to flooding and salinity, but with locally abundant populations of specialist species (McLusky & Elliott, 2004; Meire et al., 2005). Therefore, estuaries act as a reserve habitat for unique fauna (Little, 2000). In addition, strong fluxes in biogeochemical cycles support numerous functions such as high primary production and carbon fixation (Costanza, Kemp, & Boynton, 1993). Due to this high functionality, estuaries provide various ecosystem services from a local scale (e.g., food production, tourism) to a more global and drastically essential level (e.g., nutrient cycling, climate regulation; Barbier et al., 2011; Heckbert et al., 2011; Thrush et al., 2014; Hambäck et al., 2022). For centuries, estuaries have been managed to optimize commercial shipping, fishing, and economic activity, resulting in high anthropogenic levels and degradation of riverbanks (Freeman et al., 2019). In addition, land reclamation for agricultural development reduces natural habitat areas and influences natural connectivity between habitats due to water level management (McLusky & Elliott, 2004). Furthermore, modifications of the natural flow of river systems primarily affect estuarine ecosystems by altering natural freshwater inputs (Gillanders & Kingsford, 2002). This, in turn, can result in modifications of salinity and temperature regimes, affecting the abundances, distribution, and composition of resident organisms (Kennish, 2002), and could lead to functional homogenization (Cavalcante et al., 2023). All these factors mainly affect the integrity of estuaries and the ecosystem services they provide. To reconcile economic development with the ecological integrity of estuaries, specific management strategies can be employed (Boerema & Meire, 2017). In order to develop management solutions, a thorough understanding of the target ecosystem is necessary, especially because land use issues are often site dependent.

The Seine estuary in northern France plays a major ecological role with an area of about 150 km² at high tide and an average discharge of 380 m³ s⁻¹ (Avoine, 1986). It also holds significant economic importance, with a large part of the French population and agriculture concentrated along this stream (i.e., 16 million inhabitants and 40% of the country's industry and agriculture) and supports 50% of the national river traffic (Mouny et al., 1998; Lafite & Romaña, 2001). Consequently, the Seine estuary is a highly anthropized area with numerous artificial structures (dikes) limiting the effect of the tide on the associated alluvial zones. This results in a compartmentalized estuary with direct consequences on adjacent ecosystems, particularly through the modification of the lateral gradient of salinity and flooding. Consequently, natural transitions between subhalophilous meadows to freshwater meadows in alluvial zones have been modified in favor of agricultural freshwater environments (McLusky & Elliott, 2004). Additionally, the Seine estuary presents a national nature reserve including various land uses such as grazed and mowed areas and local hunting activities (Maison de l’Estuaire, 2023; Réserves Naturelles de France, 2023; Maison de l’Estuaire, 2018; ADREE, n.d.). These land uses can locally interact with water level management and impact local organisms either directly (e.g., modification of arthropod communities by grazing and mowing; Pétillon et al., 2007) or indirectly because water level management must comply with the chosen land use types.

Ecological indicators are essential to highlight the impact of water level management, local land uses, and their interaction on the Seine estuary’s biodiversity. Plants appear to be the most used bioindicators due to their high sensitivity to abiotic factors in general, including flooding and salinity stresses (McKee & Mendelssohn, 1989; Reed, 1995; Gough & Grace, 1998; Flindt et al., 1999). In the case of estuaries, plants have been proven to be relevant bioindicators (Weilhoefer, 2011). In order to better understand the link between water level management, land use, and biodiversity of alluvial zones, the use of several bioindicator taxa seems relevant. In many cases, spiders are neglected, even though their bioindication capacity is widely demonstrated (Pearce & Venier, 2006; Borchard et al., 2014). Additionally, they are present in high abundance in estuaries and are also able to indicate changes in flooding and salinity (Pétillon et al., 2003; Pétillon et al., 2014; Fournier et al., 2015; Ridel et al., 2021), but are still rarely used in this type of ecosystem (Desender & Maelfait, 1999; David et al., 2016). The coupled use of plants and spiders has the advantage of employing taxa from different trophic levels, exhibiting complementary responses to environmental changes (Lafage et al., 2015; 2019; Hacala et al., 2020). However, despite this apparent complementarity, few studies deal with these taxa as joint bioindicators (Hacala et al., 2023).

In order to properly assess ecosystem functionality, selecting complementary metrics that reflect variations in ecosystem processes as accurately as possible is an essential step. In fact, studying the local composition of community assemblages may reflect environmental filtering processes in both plants and spiders (Hacala et al., 2023). Moreover, translating this composition into a diversity metric can provide a different way to understand the impact of environmental change on organisms (Santini et al., 2017). For example, species richness combined with Shannon and Simpson indices are traditionally used and considered complementary when assessing taxonomic diversity (Chiarucci et al., 2011). These metrics reflect differences in richness and distribution of species within the assemblages and provide a good understanding of the influence of rare species' weight (Q order) on the patterns obtained in space and time (i.e., species richness: all species are equal weight; Shannon index: proportional abundance weighted; Simpson index: smaller weights attributed to rare species).

On the other hand, partitioning only taxonomic diversity metrics can lead to an incomplete assessment of the functionality of the ecosystem and/or target species (Devictor et al., 2010). To fill this gap, splitting diversity into taxonomic and functional facets can highlight more precise responses to environmental changes because species are considered not only as taxonomic units but also as organisms with their own functional characteristics (Petchey & Gaston, 2002; 2006). In many cases, taxonomic and functional diversity are linked because increasing species richness indirectly captures more functional traits (Pardo et al., 2017; Pavoine et al., 2013). However, the relationship between these metrics can provide new information for understanding biological processes like spatial variation (Devictor et al., 2010; Cadotte & Tucker, 2018). Therefore, combined diversity metrics can provide useful insights into links between biodiversity, ecosystem functions, and habitat characteristics such as environmental filtering (Díaz et al., 2007). To best describe functional diversity, numerous indices have been developed (Mouchet et al., 2010). However, when sampling small and mobile taxa such as spiders, a bias of under-sampling persists (Scharff et al., 2003), and few functional indices take this into account.

In order to properly compare the diversity of taxa without risks of under-sampling, methods have been developed to standardize sampling by ensuring its completeness (Chao et al., 2009; Chao & Jost, 2012). Chao indices have the advantage of partitioning taxonomic and functional diversity in a comparable way while still taking into account Q order (i.e., the weight given to rare species; Chao et al., 2009; 2014; Pavoine et al., 2016). Indeed, rare species can play a significant part in the overall functionality of the ecosystem (Lyons et al., 2005) because of their potential divergence from the rest of the community, supporting vulnerable functions (Mouillot et al., 2013; Jain et al., 2014). As a result, partitioning functional diversity into Q orders can highlight functional divergence between rare and abundant species, particularly in estuarine ecosystems where specialist species are in general highly abundant.

Here, we propose to investigate the influence of water management and local land uses on plants and spider assemblages and, more generally, to test the bioindicator complementarity of these taxa in estuarine environments using a case study in the Seine estuary. To our knowledge, this is the first time the bioindicator role of these taxa is explored simultaneously in a taxonomic and functional way in estuaries.

To achieve this objective, the following hypotheses will be tested:

Different patterns are expected between spider and plant assemblages due to high complementarity in the bioindicator capacity of these taxa (Lafage et al., 2015; Hacala et al., 2020; 2024). On the other hand, similar patterns are expected between these metrics for each taxon because in many cases taxonomic and functional diversity are linked (Pavoine et al., 2013; Hacala et al., 2021; Ridel et al., 2021).

Because salinity and flooding stress generally promote high abundance of specialist species (McLusky & Elliott, 2004; Meire et al., 2005), we expect an effect of Q order on taxonomic and functional diversity patterns on non-managed sites but not on managed sites for both taxa.

Due to water level management and local land uses, significant differences are expected between the composition of spider (Pétillon et al., 2003; Pétillon et al., 2014; Fournier et al., 2015; Ridel et al., 2021) and plant assemblages (McKee & Mendelssohn, 1989; Reed, 1995; Gough & Grace, 1998; Flindt et al., 1999) between sites. Additionally, differences are expected between topographical zones for each site (distance from the Seine on a lateral gradient), more pronounced on non-managed sites due to environmental filtering and indicator species in line with both water level management and local land uses.

Finally, we expect an increasing proportion of halophilic species and individuals in the non-managed sites along the topographic gradient for spiders (Pétillon et al., 2008) and plants (Adam, 1981; Bertness & Ellison, 1987; Brewer & Grace, 1990; Wilson & Stubbs, 2012; Kim & Ohr, 2020), contrasting with a stable proportion in non-managed sites.

Study sites

The study was conducted in the estuary of the river Seine (Le Havre, Normandy, France) (Fig. 1A, B), where three sites were selected for their differences in hydrological management type and lateral salinity gradient. Indeed, two of these sites are located on the north shore of the river and have water levels managed by valves and channels (called Managed 1 and Managed 2, abbreviated as M1 & M2), contrasting with the site located on the river’s south shore which is not managed and therefore directly influenced by the tidal cycle (called Non-Managed and abbreviated as NM) (Fig. 1C). Additionally, the three study sites are positioned on a lateral gradient of proximity to the sea (and therefore theoretically a salinity gradient) with M2, M1, and NM ordered from closest to furthest from the sea. On each site, three zones of 50m² were defined according to their topography: high, middle, and low zones, using a digital elevation model of the Seine estuary (GIP Seine-Aval, 2012) and vegetation community characteristics (e.g., the presence of salt-tolerant species like Bolboschoenus maritimus and Juncus gerardi). High elevation zones are inundated during exceptional events (e.g., spring tides, storm surges), contrasting with the low elevation zones closest to the river, which are regularly exposed to inundations. Additionally, all study sites are influenced by the characteristic winter flooding of the river Seine, creating a seasonal flooding gradient. However, this is mitigated by the hydraulic control of the waters from the river Seine through the valves and channels present on the north shore. In addition, the study sites present different land uses, with the NM site being grazed across all topographical levels. Site M1 is mowed with late season grazing on the high and intermediate elevation zones, while the low zones are grazed all year round. Finally, site M2 is mown across all elevation zones.

Sampling design

Spider sampling

To sample spiders, we used two plots including four sampling traps on each of the three topographical zones defined previously for a total of 24 sample points per site. At each sampling point, 50 aspirations of 2 seconds with a thermic aspirator were performed on the ground. The extremity of the g-vac measures 10 cm in diameter (i.e., 0.39 m² sampled for each point). The content of the sample was immediately stored in 70% ethanol to avoid intraspecific predation during transport. Spiders were then sorted in the lab and stored in 70% ethanol. This sampling design was repeated four times in 2020 (from 14 to 16/06, from 22 to 24/07, from 24 to 26/08, and from 14 to 16/09).

Vegetation surveys

On the same three zones along the topographical gradients of the three study sites previously defined for spider sampling, vegetation surveys were carried out within two square plots per zone, each containing a grid of 16 squares of 1 m² for a total of 96 sampling points per site. In each square, indices of vegetation cover based on Braun-Blanquet have been assigned to each plant species. These surveys were carried out once per site between 27/05 and 02/07/2020, depending on the site. To account for spatial distribution heterogeneity between taxa, the vegetation surveys were pooled (sum of plants covering) to obtain 24 sample points of 4 m² per site, and the Braun-Blanquet coefficient was converted into the Van Der Maarel coefficient, offering an intermediate solution for the weighting of rare species.

Environmental variables

Litter depth and vegetation cover were measured using the same spatial protocol applied to the vegetation surveys, with one value for each of the 16 squares defined in each plot. The thickness of the litter was measured once between 25 and 26/08/2020 with an accuracy of 0.5 cm. The vegetation cover was visually assessed once in 2020 (between 27/05/2020 and 02/07/2020). At the same time, vegetation height was also measured at 10 points regularly distributed on each of the two plots previously defined per zone for a total of 72 sampling points per site. In addition to these structural variables, pH and conductivity were measured on the superficial part of the soil (between 0 and 15 cm deep) with three measurements per zone, carried out one time between 30/06 and 01/07/2020 (using a Mettler Toledo FiveEasy sensor).

Functional traits used

In order to assess the functional diversity of plants and spiders (see below for the method), functional traits were assigned to each species following existing literature. The traits selected here are seen as indicators of essential biological functions and broadly comparable between taxa. Respectively for spiders and plants: global development (maximum size of females and maximum height), annual periodicity (seasonal activity; start of flowering periods and vegetative formation), access to resources (hunting strategies; LDMC and SLA), and dispersal ability (ballooning ability; dissemination). For spiders, traits were extracted from the literature: see (Roberts, 1985; 1987; Uetz, 1999; Harvey, Nellist, & Telfer, 2002; Bell et al., 2005; Cristofoli et al., 2010; Cardoso et al., 2011; Simonneau, Courtial, & Pétillon, 2016). Missing trait values from literature (particularly for ballooning) were completed by linking them to the closest genus available. The patterns obtained with this method remained similar to those obtained with keeping missing values. Average plant leaf traits were extracted from the TRY database (Kattge et al., 2020) and root traits from the Global Root Traits database (GRooT; Guerrero-Ramírez et al., 2021). They were completed with overall life strategies from the baseflor database (Julve, 2018). A summary of the modalities of each trait used here is attached in Appendix 1.

Statistical analyses

Habitat characteristics

In order to highlight whether the environmental variables differed between the three topographical zones of each of the sites, ANOVA tests followed by Tukey post-hoc tests or Kruskal–Wallis followed by Dunn tests were carried out as appropriate. Holm correction was used for post-hoc tests if necessary. All statistical were performed using R STUDIO software (v. 4.3.2). and packages detailed hereafter.

Taxonomic vs functional diversity patterns and Q order influence

Taxonomic diversity was estimated with the INEXT.3D package (v. 1.0.1) (Chao et al., 2021) for each Q diversity order (Q0, Q1, and Q2), corresponding to species richness, the Shannon diversity index, and the Simpson diversity index. These indices were calculated for each topographical zone and for each site for both taxa using the iNEXT3D function with 50 bootstrap replicates. This method considers the sample coverage and therefore avoids the risk of false conclusions due to insufficient sampling. For spiders, the individuals present in the different sampling sessions were pooled. Similarly, functional diversity was estimated using the INEXT.3D package for each Q order previously mentioned, defined as being comparable to the calculation method used for taxonomic diversity (for more details on the methods, see Chao, Chiu, & Jost, 2014; Chao et al., 2021). This method defines the absence of overlap in the confidence interval of produced curves as a significant difference.

Differences in community composition and species richness

To test differences in terms of composition and abundance between sites for each taxon, PERMANOVA was performed with the adonis2 function (vegan package, v.2.6-4) in a Bray-Curtis distance matrix. If tests signaled significant differences, multiple comparison tests were carried out with the pairwise.adonis2 function (pairwiseAdonis package) on the model.

Indicator species

To identify indicator species of each topographic zone per site, the indicator index (IndVal) proposed by Dufrene and Legendre (1997) was calculated using the multipatt function (indicspecies package, v. 1.7.14). According to the described methods, a threshold level of 25% for the index was considered significant, indicating the presence of the species in more than 50% of the traps in targeted topographic zones as well as a relative abundance of more than 50% in this zone.

Proportion of halophilic species

To compare the proportion of halophilic species between topographic zones per site, the halophilic character was attributed to each spider species according to the literature and coded with binary responses (halophilic species vs. non-halophilic species). Vegetation tolerance to salinity was assessed using average Ellenberg indicator values for salt tolerance from the TRY database (Kattge et al., 2020). Afterwards, community-weighted means were calculated with the functcomp function (FD package, v. 1.0-12.3) on a presence/absence matrix (equal weight for each species) and on an abundance matrix for both taxa (abundance proportional weight). Then, significant differences between topographic zones were assessed by ANOVA tests followed by Tukey post-hoc tests or Kruskal–Wallis followed by Dunn tests with Holm correction where appropriate.

Habitat characteristics

As expected, the non-managed site NM featured an increase in salinity from the high to the low zones as well as a slightly lower vegetation cover in the low areas (Table 1). The managed site M1 is characterized by a significant increase in the litter depth from the high to the low areas, high global vegetation cover, and similar conductivity measured between each topographic zone. Salinity standard deviation of the high topographic zone appeared particularly strong on this site. Finally, the managed site M2 is characterized by a significant increase in the litter depth and smaller vegetation height in the middle and low zones. Moreover, salinity appeared similar between each topographic zone and surprisingly, the measured values were low.

Tab.1 Environmental variables (mean ± s.d., n = 8) for each zone and for each site. (Successive letters indicate significant differences. H, high ; M, middle; L, low. NM = Non Managed site, M1 = Managed site one and M2 = Managed site two

		NM				M1				M2
average height of vegetation (cm)	H	38.6	±	14.5	a	97.0	±	14.7	A	101.0	±	22.7	x
	M	34.4	±	8.5	a	90.6	±	10.4	A	61.6	±	22.9	y
	L	41.3	±	12.2	a	60.8	±	6.1	B	56.6	±	16.2	y
total vegetation cover (%)	H	100.0	±	0.0	a	100.0	±	0.0	A	100.0	±	0.0	xy
	M	99.8	±	0.4	a	100.0	±	0.0	A	100.0	±	0.0	x
	L	82.3	±	10.5	b	100.0	±	0.0	A	99.2	±	2.5	y
litter depth (cm)	H	0.0	±	0.0	a	0.7	±	0.5	A	0.2	±	0.3	x
	M	0.0	±	0.0	a	0.3	±	0.4	B	1.5	±	0.8	y
	L	0.0	±	0.0	a	1.1	±	0.7	C	2.0	±	0.7	z
soil conductivity (µS/cm)	H	524.8	±	101.2	a	1361.7	±	1527.3	A	427.9	±	164.8	x
	M	1010.9	±	212.6	b	1186.3	±	305.3	A	225.6	±	40.0	x
	L	1917.3	±	148.5	c	2531.3	±	618.9	A	744.8	±	365.0	x
pH	H	7.8	±	0.1	a	7.9	±	0.3	A	7.9	±	0.1	x
	M	8.1	±	0.0	b	8.0	±	0.1	A	8.0	±	0.0	x
	L	8.2	±	0.1	b	7.7	±	0.2	A	7.9	±	0.1	x

Taxonomic vs functional diversities patterns and Q order influence

On the non-managed site NM, species richness of spiders were similar for each topographic zone due to the high variability of estimated values (especially for high and low zones) (Fig. 2A). The Shannon and Simpson diversity indices were found to be highest in the low zones and lowest in the middle and low zones, without a statistically significant difference between the two indices. For plants, a similar pattern was obtained for each Q order, with diversity significantly decreasing from the upper zone to the lower (Fig. 2B). On this site, the functional richness of spiders showed a similar pattern to species richness with no significant difference between each topographic zone (Fig. 3A). Similarly, functional diversity for the other Q orders showed maximal values in the low zones and minimal values for middle and low zones without significant difference. For plants, functional richness appeared significantly higher in the upper zone but similar in middle and low zones. This pattern remained true for each Q order (Fig. 3B).

In managed site M1, spider species richness (q = 0) appeared higher in the low zone and similar between the middle and high zones (Fig. 4A). This pattern was also found for the other Q orders. In contrast, the species richness of plants appeared maximal for the high zone, intermediate for the low zone, and minimal for the middle zone (Fig. 4B). Moreover, for this taxon, specific diversity remained at a maximum in the upper zone according to all the Q orders, but the middle zone was respectively lower than the low zone for species richness, equal for Shannon diversity, and higher for Simpson diversity.

Spider functional richness was recorded as similar between each topographic zone, contrasting with species richness (Fig. 5A). For other Q orders, higher values were obtained for the lower zone and similar values for high and middle zones, similar to those obtained for taxonomic diversity. For plants, all functional diversity values between Q orders appeared highest for the high elevation zone, intermediate for the low elevation zone, and lowest for the middle elevation zone (Fig. 5B). This pattern is congruent with those obtained for species richness but different from those obtained for Shannon and Simpson diversity.

On managed site M2, spider species richness (q = 0) in the high elevation zone appeared higher than in the low elevation zone. The middle zone displayed intermediate and similar values (Fig. 6A). In contrast, for Shannon diversity (q = 1), maximal diversity was still observed in the upper elevation zone while the diversity of the middle and lower elevation zones was lower and similar to one another. Regarding Simpson diversity (q = 2), the high and middle elevation zones showed different values but no significant difference was noted between low vs high or middle elevation zones. For plants, a similar pattern was obtained for each diversity Q order with significant differences between all topographic zones. The highest value was obtained for the middle zone, intermediate value for the high elevation zone, and minimal value for the low elevation zone (Fig. 6B).

On this site, spider functional richness (q = 0) appeared similar for each topographic zone (Fig. 7A). For the other Q orders, the diversity values of the middle elevation zones stood out, being significantly different from the high and low elevation zones. For plants, the global pattern does not change between diversity Q orders, with the highest value obtained for the middle zone, intermediate value for the high zone, and minimal value for the low zone, consistent with taxonomic diversity (Fig. 7B).

Composition

The composition and abundance of both spiders and plants showed significant differences between sites (p-value systematically lower than 0.001; see appendix 2 for more details on test values). More precisely, site by site, all topographic zones appeared significantly different from each other for spider composition (p-value systematically lower than 0.001 except for the managed site M1 high vs middle zones p = 0.003). Similarly, all topographic zones appeared significantly different from each other for plant composition (p-value systematically lower than 0.001 except for managed site M2 high vs middle zones p = 0.002) (see appendices 3, 4, and 5 for more details on test values).

Indicator species

In the non-managed site NM, two spider species appeared as significant indicator species for the high elevation zone (Agyneta mollis and Pachygnatha clercki), none for the middle zone, and just one (Pardosa purbeckensis) in the low elevation zone (Table 2). For plants, six species (Hordeum secalinum; Phleum pratense; Trifolium repens; Potentilla anserina; Cynosurus cristatus; Holcus lanatus) were identified as significantly indicative for the high zone and none for the middle and low zones (Table 3).

Tab.2 Significant indicator spider’s species foreach topographic zones (H= Hight, M = Middle and L = low) per study site, and associated p-value. Halophilic species are in bold. NM = Non Managed site, M1 = Managed site one and M2 = Managed site two. Signif. codes: *** p<0.001 ** p<0.01 * p<0.05

	NM		M1		M2
	species	p	species	p	species	p
H	Agyneta mollis	0.030 *	Ozyptila sanctuaria	0.015 *	Agyneta mollis	0.001 ***
	Pachygnatha clercki	0.032 *			Tenuiphantes tenuis	0.001 ***
					Pardosa proxima	0.001 ***
					Agyneta rurestris	0.001 ***
					Oedothorax retusus	0.013 *
					Pardosa palustris	0.024 *
					Enoplognatha mordax	0.009 **
M	/		/		Gnathonarium dentatum	0.012 *
L	Pardosa purbeckensis	0.001 ***	Pardosa proxima	0.001 ***	Pirata piraticus	0.003 **
			Pardosa prativaga	0.001 ***
			Piratula latitans	0.003 **
			Arctosa leopardus	0.018 *

Tab.3 Significant indicator plants species for each topographic zone (H= Hight, M = Middle and L = low) per study site, and associated p-value. Ellenberg salinity value are reported for each species. NM = Non Managed site, M1 = Managed site one and M2 = Managed site two. Signif. codes: *** p<0.001 ** p<0.01 * p<0.05

	NM			M1			M2
	species	p	Ellenberg salinity	species	p	Ellenberg salinity	species	p	Ellenberg salinity
H	Hordeum secalinum	0.001***	2.00	Arrhenatherum elatius	0.001***	0.50	Bromus racemosus	0.001***	0.50
	Phleum pratense	0.001***	0.33	Holcus lanatus	0.001***	0.67	Festuca arundinacea	0.001***	1.83
	Trifolium repens	0.003**	1.25	Plantago lanceolata	0.001***	0.50	Holcus lanatus	0.001***	0.67
	Potentilla anserina	0.001***	2.25	Potentilla anserina	0.001***	2.25	Lolium perenne	0.001***	0.33
	Cynosurus cristatus	0.006**	0.00	Carex distans	0.001***	4.00	Cynosurus cristatus	0.001***	0.00
	Holcus lanatus	0.019*	0.67	Trifolium pratense	0.001***	0.60
				Plantago major	0.002**	0.60
				Juncus gerardii	0.007**	6.00
				Crepis biennis	0.008**	0.33
				Cynosurus cristatus	0.014*	0.00
				Lotus corniculatus	0.017*	0.75
M	/			Phleum pratense	0.001***	0.33	Juncus articulatus	0.001***	1.00
							Carex otrubae	0.001***	2.00
							Juncus gerardii	0.001***	6.00
							Myosotis laxa	0.008**	0.00
							Alopecurus geniculatus	0.004**	1.25
L	/			Ranunculus sardous	0.001***	2.20	/
				Atriplex prostrata	0.001***	1.80
				Oenanthe fistulosa	0.002**	0.33

In managed site M1, only one spider species (Ozyptila sanctuaria) was noted as significantly indicative for the high elevation zone, none for the middle zone, and four for the low elevation zone (Pardosa proxima; Pardosa prativaga; Piratula latitans; Arctosa leopardus). Interestingly, all indicator spider species of the low elevation zone belong to the same family, the Lycosidae. Concerning plants, indicator species showed a different pattern to that obtained for spiders on this site, with a large number of indicator species in the upper elevation zone (11 species, see Table 3 for details), then a single species (Phleum pratense) for the middle zone and three species for the lower zone (Ranunculus sardous; Atriplex prostrata; Oenanthe fistulosa).

Concerning managed site M2, a large number of spiders species were recorded as indicators of the upper elevation zone with seven species distributed across three families (see Table 2 for details). For the middle and lower elevation zones of this site, only one spider species, was recorded as an indicator species (respectively: Gnathonarium dentatum and Pirata piraticus). For plants, an equal number of five indicator species were found for the high and middle zones (respectively: Bromus racemosus; Festuca arundinacea; Holcus lanatus; Lolium perenne; Cynosurus cristatus and Juncus articulates; Carex otrubae; Juncus gerardii; Myosotis laxa; Alopecurus geniculatus), while no species were identified as indicative of the lower zone.

Salinity affinity

On the non-managed site NM, the proportion of halophilic spider species showed increasing values from the high elevation to the low elevation zone, with significant differences observed among all values obtained for the topographical zones (Fig. 8.A1). Interestingly, when weighted by abundance, only the lower elevation zone appeared to have significantly fewer species compared to the middle and high elevation zones, with no significant difference between the latter (Fig. 8.A2). Similarly, for plants, the mean Ellenberg salinity value followed the same pattern as observed for spiders, with significantly increasing values from the high elevation to the low elevation zones (Fig. 8.B1). The results obtained with abundance-weighted values appeared consistent with this pattern (Fig. 8.B2).

On managed site M1, there were no significant differences observed in the proportion of halophilic spider species present in the assemblages among the various topographical zones (Fig. 8.A1). This pattern was also evident for the proportion of halophilic individuals when weighted by species abundance (Fig. 8.A2). However, for plants, the mean Ellenberg salinity value of the assemblage appeared maximal for the high and low elevation zones, with no significant difference between them, while significantly fewer halophytes were found in the middle zone (Fig. 8.B1). A similar pattern was observed for the abundance-weighted metric (CWM) (Fig. 8.B2).

On managed site M2, the proportion of halophilic spider species appeared to be higher in the high elevation zone, lower in the low elevation zones, and intermediate but not significantly different from each other in the middle zone (Fig. 8.A1). This pattern remained consistent when abundance-weighted metrics were considered (Fig. 8.A2). For plants, the mean Ellenberg salinity value showed significantly higher values for the middle zone compared to the others, with similar values observed between the high and low elevation zones (Fig. 8.B1). However, when abundance-weighted values were considered, the middle zone still appeared significantly higher, with the high elevation zone showing an intermediate value and the low elevation zone showing a significantly lower value compared to the upper zone (Fig. 8.B2).

Habitat characteristics

For the non-managed site NM, the observed salinity gradient aligns with the absence of water level management practices. The lower vegetation cover in the low topographical zones may to be a consequence of grazing on wetlands, as livestock has been shown to increase the frequency and duration of waterlogging through trampling and soil compaction (Dausse et al. 2012). In managed site M1, the significant increase in litter depth could be attributed to the water management practices. During the winter period, water valves are closed (Reserve Naturelle Estuaire de Seine 2023), leading to water retention and consequently accumulation of litter that cannot be discharged into the Seine River. The similarity of salinity levels across each topographic zone in this site corresponds to the water management practices. Moreover, the variability in salinity measurements in the high elevation zones could be attributed to dredging activity in this area, resulting in soil with heterogeneous properties, as supported by the high variability in granulometry measured in another research work (Neupert et al., in press).

Furthermore, the higher litter depth observed in the lower topographic zone of managed site M2 aligns with the previous observation of water retention during winter, and the uniformity of salinity throughout the site is consistent with water level management. However, the low soil salinity values measured on this site were unexpected given its proximity to the river mouth, and could be indicative of an hydrological anomaly, possibly linked to the upwelling of freshwater from the karst aquifers of the surrounding chalk cliffs, as suggested by Neupert et al. (in press) and Soueid-Ahmed et al. (2017).

Taxonomic vs functional diversity patterns

For the non-managed site NM, a contrasting pattern was observed between the taxonomic diversity of spiders and plants, supporting our initial hypothesis of high complementarity between indicator taxa (Lafage et al., 2015; Hacala et al., 2020; Hacala et al., 2024). Specifically, while a similar pattern was observed for spider species richness across topographical levels, decreasing richness values were obtained for plants with decreasing elevation, typically indicative of a stress gradient (i.e., reduction in species richness due to increasing abiotic filter strengths). This result could reflect a turnover process in spider assemblages without influencing taxonomic richness, possibly due to insufficient salinity to filter only halophilic species, as observed in harsher environments like salt marshes (Pétillon et al., 2003). These results suggest a weaker specific selection on spider species compared to plant species on this site, possibly due to their higher mobility compared to sessile organisms (Lafage et al., 2015). The pattern of spider functional richness, similar to species richness, is consistent with the high correlation between these two metrics (Pavoine et al., 2013), supporting species selection based on functional characteristics, as previously observed for spiders in harsher environments (Ridel et al., 2021). Conversely, the similar functional richness of plants in the middle and low elevation zones, contrasting with previous results, suggests functional trait convergence among present species (Meinzer, 2003).

In managed site M1, the complementarity of spider versus plant taxonomic diversity was also evident. Specifically, the higher spider species richness observed in the lower zone could be linked to the increased litter depth, resulting from limited discharge of organic matter into the Seine River, which may enhance spider species diversity by reducing interspecific competition (Döbel et al., 1990). Conversely, the highest species richness was recorded in the high elevation zones, possibly due to the high salinity variability, suggesting co-occurrence of species with different ecological affinities by reducing interspecific competition (Crain et al., 2004). This heterogeneity could be linked to dredging activity, as previously mentioned. Contrasting patterns were also noted in functional diversity between plants and spiders, highlighting their complementary functional roles (Hacala et al., 2024). Interestingly, spider functional richness showed no difference between topographic zones, contrasting with taxonomic richness results, suggesting a diversification of species that are functionally similar, possibly linked to similar environmental structures driving spider functional assemblages (Leroy et al., 2014). Conversely, for plants, the pattern obtained for functional richness aligned with previous results, supporting a global redundancy between these metrics for this taxon (Pavoine et al., 2013).

Managed site M2 remained consistent with other study sites regarding the observed complementarity between taxa. Specifically, the reduction in spider species richness from high to low topographic zones is consistent with environmental filtering processes (Pétillon et al., 2008). This result contrasts with the water level management of this site by valves but could be linked to freshwater resurgence (i.e., upwelling of freshwater from karst aquifers), as spider assemblages are sensitive to flooding and/or hydric soil conditions (Fournier et al., 2015). Conversely, plant species richness appeared higher in the middle zone, possibly linked to freshwater upwelling coupled with restrictions on tidal influence, generating intermediate environmental conditions allowing the co-occurrence of different species types (i.e., halophilic and hygrophilous species). In terms of functionality, the spider pattern remained generally coherent with the taxonomic one, with minor differences observed between pairs of topographical zones but no shifting patterns (Pavoine et al., 2013; Ridel et al., 2021). In contrast to this result, plant diversity in the lower zone appeared to shift from lower taxonomic to intermediate functional richness values. This result could be linked to the partial selection of flooding specialist species in the lower zone, reducing taxonomic but increasing functional diversity due to high divergence between functional hygrophilous and salt-tolerant species.

Q order influence

As expected, on the non-managed site NM, the patterns of taxonomic and functional diversity of spiders changed between Q orders. Surprisingly, for abundance-weighted metrics, diversity in the low elevation zone appeared higher. This suggests balanced proportions of halophilic and non-halophilic species near the Seine River, possibly due to an increasing proportion of specialist species in this zone. This supports lightly filtered spider assemblages (i.e. no exclusion of generalist species) compared to harsher environments (Pétillon et al., 2008; Ridel et al., 2021), which is consistent with the site's position furthest from the sea on the lateral gradient. For plants, taxonomic and functional diversity patterns did not significantly change with Q order on this site, suggesting a high turnover process (Yuan et al., 2012). This interesting contrast between spiders and plants highlights different spatial responses to environmental filtering.

On managed site M1, the patterns observed for spider diversity remain consistent for each Q order, in line with the initial hypothesis (Hacala et al., 2023). In contrast, functional patterns for spiders diverge from this previous stability, supporting the hypothesis of functional convergence mentioned above. However, for plants, the shifting of taxonomic diversity values between middle and low elevation zones as Q order increases supports unbalanced abundance in the lower zones, possibly due to winter flooding in that area (Fournier et al., 2015). On the other hand, functional plant diversity remains stable between each Q order, suggesting functional redundancy between some rare and abundant species (Meinzer, 2003).

On managed site M2, the shifting position of spider taxonomic and functional diversity in the middle zone indicates unbalanced abundance. In this site, the upwelling of freshwater from the water table creates longer periods of waterlogging and reduced soil salinity, which could favor the presence of hygrophilous species. Conversely, taxonomic and functional patterns of plants remain stable between Q orders, suggesting a higher influence of environmental filtering on composition than on abundance for this taxon at the spatial scale of the defined topographical gradient. Interestingly, a general taxonomic and functional convergence was observed for spiders on this site, coupled with an effect of Q order on functional metrics, contrasting with the taxonomic and functional divergence observed for plants, with an effect of Q order on taxonomic metrics. This underscores the strong complementarity between the bioindicator taxa used (Lafage et al., 2015; Hacala et al., 2020; 2024).

Composition of assemblages

According to our initial hypothesis, there are differences in the specific composition and abundance between sites, indicating the sensitivity of spiders and plants to water management types (non-managed vs managed), consistent with their sensitivity to flooding and salinity (respectively Pétillon et al., 2003; Pétillon et al., 2014; Fournier et al., 2015; Ridel et al., 2021, and McKee & Mendelssohn, 1989; Reed, 1995; Gough & Grace, 1998; Flindt et al., 1999). However, differences in composition between managed sites highlight the sensitivity of these taxa to local parameters and exploitation types (mowing and grazing). Site by site, the differences in composition partially contrast with the results on diversity (e.g., different composition but similar diversity obtained in spider taxonomic diversity between high and middle elevation zones on managed site M1). This opposition can indicate shifting assemblages, sometimes based on composition rather than on species proportion, highlighting high turnover, as observed in other coastal systems (see for spiders: Pétillon et al., 2008; Coccia & Fariña, 2019 and for plants: Janousek & Folger, 2014; Lawrence et al., 2022).

Indicator Species

On the non-managed site NM, the two spiders indicator species of the high elevation zone (Agyneta mollis and Pachygnatha clercki) were common and generalist species (Hänggi et al., 1995), indicating low-filtered zones. For the lower elevation zones, the presence of Pardosa purbeckensis, a species occurring in shores (Puzin et al., 2014), is consistent with the management type and environmental salinity gradient measured here. For plants, the presence of sub-hygrophilic species indicators of high elevation zones (e.g., Hordeum secalinum, Potentilla anserina) with the absence of strict halophilic species aligns with the results obtained for spiders.

On managed site M1, the only spider indicator species found for high elevation zone assemblages, Ozyptila sanctuaria, is typically found in open habitats (Dawson et al., in prep), consistent with the management type (i.e., mowing with late-season grazing). In the lower elevation zones, the presence of four Lycosidae species supports the diversification of functionally close species linked to an increase in litter depth by adding new prey guilds, favoring ground-hunting species (Döbel et al., 1990, Uetz, 1991). Finally, the presence of Pirata piraticus indicates wet conditions in this zone (Harvey et al., 2002). For plants, assemblages show numerous indicator species of high zones, with sub-halophilic species (e.g., Juncus gerardi) mixed with non-halophilic species (e.g., Trifolium pratense) (Julve, 1998). This co-occurrence of species with different life strategies is congruent with the high taxonomic and functional diversity previously obtained, supporting the strong heterogeneity of this zone. In the lower elevation zone, the three species found (i.e., Ranunculus sardous, Atriplex prostata, and Oenanthe fistulosa) are hygrophilic species typical of flooded meadows (Julve, 1998), consistent with the wet conditions mentioned above.

Finally, for managed site M2, spiders indicator species in higher zones are principally generalist species, except Enoplognatha mordax, a coastal shore species. However, caution is necessary regarding the presence of this halophilic species alone, as it is sometimes found in very different habitats such as agroecosystems (Djoudi et al., 2018), and it is possible that E. mordax forms a specific complex (Bosmans & Van Keer, 1999) encompassing strictly halophilic and non-halophilic forms. For middle and low zones, the presence of Gnathonarium dentatum and Pirata piraticus, two species occurring in flooded habitats (Harvey et al., 2002), supports a high hydric stress gradient, as previously discussed (i.e., freshwater resurgence from the water table). For plants, only mesophilic meadow-characteristic plants were found in the high elevation zones (e.g., Holcus lanatus, Lolium perenne) (Julve, 1998), illustrating their low environmental constraint. However, in the middle elevation zone, the mixture of hygrophilous plant species (e.g., Carex otrubae) (Julve, 1998) and halophilic plants like Juncus gerardi contrasts with the absence of halophilic spider species, highlighting the higher sensitivity of this taxa, as previously observed. Moreover, this co-occurrence of different species types (i.e., halophilic and hygrophilous) is consistent with the higher taxonomic and functional diversity previously obtained.

Overall, the absence of indicator species from some topographical zones for each taxa is due to the absence of some exclusive species, indicating transitional assemblages, as suggested by the presence of indicative species of pooled zones (Appendix 6 and 7).

Proportion of halophilic species

For spiders, the results obtained are generally consistent with expected patterns, with constant proportions of halophilic species on each managed site and an increasing proportion in the non-managed site NM. An exception should be noted for managed site M2, with higher rates of halophilic species and individuals in the high topographic zones. This result is exclusively linked to the presence of E. mordax, listed as a halophilic species. As suggested above, this could be a non-halophilic form of this species here. Further work is needed on this genus given the contrasting ecologies and distributions described for this species (Bosmans & Van Keer, 1999). As expected, when metrics are abundance-weighted, the pattern changed only on the non-managed site NM. However, the results obtained reflect a specific composition with more specialist species in lower elevation zones (Pétillon et al., 2008), but with an equal proportion of specialist individuals, corresponding to lightly filtered spider assemblages without the exclusion of generalist species.

Similarly to spiders, the salinity tolerance index for plants (Ellenberg indicator values) increases near the river Seine on the non-managed site NM, but significant variations were also found on the two managed sites. This more precise response of plants in terms of halophilic species proportion confirms the idea of a generally higher sensitivity to salinity and flooding than spiders. However, this result could be due to the still limited knowledge of spider traits (Pekár et al., 2021), allowing for a qualitative implementation of halophilic affinity for this taxon. Additionally, and contrasting with spiders, the plant salinity tolerance index is broadly consistent with taxonomic and functional diversity values, supporting high heterogeneity and co-occurrence of different species types in some places (e.g., high elevation zone of managed site M1 and middle elevation zone of managed site M2).

In conclusion, on the unmanaged site NM, spider assemblages appeared lightly influenced by salinity, exhibiting a turnover of species along the stress gradient and balanced abundance between halophilic and non-halophilic species. Conversely, plant assemblages appeared more strongly influenced by filtering processes, with species selection and functional homogenization from the middle elevation zone towards the river Seine. On managed site M1, spiders seemed to be influenced by local factors, leading to species diversification but functional convergence in the lower elevation zone. The flooding constraints present on this site appeared insufficient to strongly influence these taxa, as evidenced by the unbalanced taxonomic abundance. In contrast, plants appeared more sensitive, exhibiting the co-occurrence of different ecotypes possibly linked with winter flooding, along with functional homogenization in the lower elevation zone indicating stronger environmental filtering. Finally, on managed site M2, resurgence from the water table impacted spider assemblages, showing a response to flooding. Conversely, plants exhibited mixed halophilic and hygrophilous assemblages on this site, illustrating the vegetaiton’s higher sensitivity to low salinity levels.

Author Contribution

Conception and design: J.P. Acquisition of data: A.R. and M.NAnalysis and interpretation: A.R. and J.P Drafting: A.R. Revising the article: J.P., M.N and E.L.

Acknowledgement

We would like to thank the staff of the Seine estuary nature reserve and Olivier Jambon help during fieldwork. This project was funded by the GIP 'Seine-AVal' (Projet FEREE « Comparaison du Fonctionnement Écologique de secteurs intertidaux contrastés pour la compréhension de leurs connectivités et la Restauration des fonctions Écologiques Estuariennes »).

Adam P (1981) The vegetation of British saltmarshes. New Phytol 88(1):143–196
ADREE s. d. Cartographie des sols de la Réserve Naturelle Nationale de l’estuaire de la Seine
Avoine J (1986) Evaluation des apports fluviatiles dans l’estuaire de la Seine. Act Colloq IFREMER 4:117–124
Barbier EB, Hacker SD, Kennedy C, Koch EW, Stier AC, Silliman BR (2011) The value of estuarine and coastal ecosystem services. Ecol Monogr 81:169–193. https://doi.org/10.1890/10-1510.1
Bell JR, Bohan DA, Shaw EM, et, Weyman GS (2005) Ballooning Dispersal Using Silk: World Fauna, Phylogenies, Genetics and Models. Bull Entomol Res 95(2):69–114. https://doi.org/10.1079/BER2004350
Bertness MD, Aaron ME (1987) Determinants of Pattern in a New England Salt Marsh Plant Community. Ecol Monogr 57(2):129–147. https://doi.org/10.2307/1942621
Boerema A, Meire P (2017) Management for Estuarine Ecosystem Services: A Review. Ecol Eng 98:172–182. https://doi.org/10.1016/j.ecoleng.2016.10.051
Borchard F, Buchholz S, Helbing F, Fartmann T (2014) Carabid Beetles and Spiders as Bioindicators for the Evaluation of Montane Heathland Restoration on Former Spruce Forests. Biol Conserv 178:185–192. https://doi.org/10.1016/j.biocon.2014.08.006
Bosman R, Van Keer J, J (1999) The Genus Enoplognatha Pavesi 1880 in the Mediterranean Region Araneae: Theridiidae. Bull Br Arachnol Soc 11(6):209–241
Brewer JS, Grace JB (1990) Plant community structure in an oligohaline tidal marsh. Vegetation 90(2):93–107. https://doi.org/10.1007/BF00033019
Cadotte MW, Tucker MT (2018) Difficult Decisions: Strategies for Conservation Prioritization When Taxonomic, Phylogenetic and Functional Diversity Are Not Spatially Congruent. Biol Conserv 225:128–133. https://doi.org/10.1016/j.biocon.2018.06.014
Cardoso P, Pekár S, Jocqué R, Jonathan A, Coddington (2011) Global Patterns of Guild Composition and Functional Diversity of Spiders. PLoS ONE 6(6):e21710. https://doi.org/10.1371/journal.pone.0021710
Cavalcante LL, Salete Daga V, Rennó Braga R, Andrian Padial A (2023) Functional Homogenization in Aquatic Ecosystems: A Review and Framework Proposal. Hydrobiologia 850(6):1283–1302. https://doi.org/10.1007/s10750-022-04919-4
Chao A, Chiu CH, Jost L (2014) Unifying Species Diversity, Phylogenetic Diversity, Functional Diversity, and Related Similarity and Differentiation Measures Through Hill Numbers. Annu Rev Ecol Evol Syst 45(1):297–324. https://doi.org/10.1146/annurev-ecolsys-120213-091540
Chao A, Colwell RK, Lin CH, Gotelli J, N (2009) Sufficient Sampling for Asymptotic Minimum Species Richness Estimators. Ecology 90(4):1125–1133. https://doi.org/10.1890/07-2147.1
Chao A, Gotelli NJ, Hsieh TC, Sander EL, Ma KH, Colwell RK, Ellison AM (2014) Rarefaction and Extrapolation with Hill Numbers: A Framework for Sampling and Estimation in Species Diversity Studies. Ecol Monogr 84(1):45–67. https://doi.org/10.1890/13-0133.1
Chao A, Henderson PA, Chiu C, Moyes F, Hu K, Dornelas M, Magurran AE (2021) Measuring temporal change in alpha diversity: A framework integrating taxonomic, phylogenetic and functional diversity and the iNEXT.3D standardization. Methods Ecol Evol 12:1926–1940. https://doi.org/10.1111/2041-210X.13682
Chao A, Jost L (2012) Coverage-based rarefaction and extrapolation: standardizing samples by completeness rather than size. Ecology 93:2533–2547. https://doi.org/10.1890/11-1952.1
Chiarucci A, Bacaro G, Scheiner SM (2011) Old and new challenges in using species diversity for assessing biodiversity. Phil Trans R Soc B 366:2426–2437. https://doi.org/10.1098/rstb.2011.0065
Coccia C, Fariña JM (2019) Partitioning the effects of regional, spatial, and local variables on beta diversity of salt marsh arthropods in Chile. Ecol Evol 9:2575–2587. https://doi.org/10.1002/ece3.4922
Costanza R, Kemp WM, Boynton WR (1993) Predictabiliw, Scale, and Biodiversity in Coastal and Estuarine Ecosystems: Implications for Management. 88–96
Crain CM, Silliman BR, Bertness SL, Bertness MD (2004) PHYSICAL AND BIOTIC DRIVERS OF PLANT DISTRIBUTION ACROSS ESTUARINE SALINITY GRADIENTS. Ecology 85:2539–2549. https://doi.org/10.1890/03-0745
Cristofoli S, Mahy G, Kekenbosch R, Lambeets K (2010) Spider communities as evaluation tools for wet heathland restoration. Ecol Ind 10:773–780. https://doi.org/10.1016/j.ecolind.2009.11.013
Dausse A, Garbutt A, Norman L, Papadimitriou S, Jones LM, Robins PE et al (2012) Biogeochemical functioning of grazed estuarine tidal marshes along a salinity gradient. Estuarine, Coastal and Shelf Science, Recent advances in biogeochemistry of coastal seas and continental shelves. 100:83–92
David AT, Goertler PAL, Munsch SH, Jones BR, Simenstad CA, Toft JD, Cordell JR, owe ER, Gray A, Hannam MP, Matsubu W, Morgan EE (2016) Influences of Natural and Anthropogenic Factors and Tidal Restoration on Terrestrial Arthropod Assemblages in West Coast North American Estuarine Wetlands. Estuaries Coasts 39:1491–1504. https://doi.org/10.1007/s12237-016-0091-3
Davidson NC (1991) Nature Conservation and Estuaries in Great Britain. Note. The Full Book Is in 4 Pdf Parts, Downloadable from: Http://Jncc.Defra.Gov.Uk/Page-2563. Nature Conservancy Council/JNCC. https://doi.org/10.13140/2.1.3522.9448
Desender K, Maelfait J-P (1999) Diversity and conservation of terrestrial arthropods in tidal marshes along the River Schelde: a gradient analysis. Biol Conserv 87:221–229. https://doi.org/10.1016/S0006-3207(98)00058-5
Devictor V, Mouillot D, Meynard C, Jiguet F, Thuiller W, Mouquet N (2010) Spatial mismatch and congruence between taxonomic, phylogenetic and functional diversity: the need for integrative conservation strategies in a changing world: Spatial mismatch between diversity facets. Ecol Lett 13:1030–1040. https://doi.org/10.1111/j.1461-0248.2010.01493.x
Díaz S, Lavorel S, Chapin FS, Tecco PA, Gurvich DE, Grigulis K (2007) Functional Diversity at the Crossroads between Ecosystem Functioning and Environmental Filters. In: Canadell JG, Pataki DE, Pitelka LF (eds) Terrestrial Ecosystems in a Changing World. Global Change — The IGBP Series. Springer, Berlin Heidelberg, Berlin, Heidelberg, pp 81–91. https://doi.org/10.1007/978-3-540-32730-1_7
Djoudi EA, Marie A, Mangenot A, Puech C, Aviron S, Plantegenest M, Pétillon J (2018) Farming system and landscape characteristics differentially affect two dominant taxa of predatory arthropods. Agric Ecosyst Environ 259:98–110. https://doi.org/10.1016/j.agee.2018.02.031
DÖbel HG, Denno RF, Coddington JA (1990) Spider (Araneae) Community Structure in an Intertidal Salt Marsh: Effects of Vegetation Structure and Tidal Flooding. Environ Entomol 19:1356–1370. https://doi.org/10.1093/ee/19.5.1356
Dufrene M, Legendre P (1997) Species assemblages and indicator species. Ecol Monogr 67:345–366
Flindt MR, Pardal MÂ, Lillebø AI, Martins I, Marques JC (1999) Nutrient cycling and plant dynamics in estuaries: A brief review. Acta Oecol 20:237–248. https://doi.org/10.1016/S1146-609X(99)00142-3
Fournier B, Gillet F, Le Bayon R, Mitchell EAD, Moretti M (2015) Functional responses of multitaxa communities to disturbance and stress gradients in a restored floodplain. J Appl Ecol 52:1364–1373. https://doi.org/10.1111/1365-2664.12493
Freeman LA, Corbett DR, Fitzgerald AM, Lemley DA, Quigg A, Steppe CN (2019) Impacts of Urbanization and Development on Estuarine Ecosystems and Water Quality. Estuaries Coasts 42:1821–1838. https://doi.org/10.1007/s12237-019-00597-z
Gillanders B, Kingsford M (2002) Impact of Changes in Flow of Freshwater on Estuarine and Open Coastal Habitats and the Associated Organisms. In: Atkinson R, Barnes M (eds) Oceanography and Marine Biology, An Annual Review. Oceanography and Marine Biology - An Annual Review, vol 40. CRC, pp 233–309. https://doi.org/10.1201/9780203180594.ch5
GIP Seine-Aval (2012) Levé topographique haute résolution de l’estuaire de la seine. Acquisition par lidar aéroporté et exploitation des données. [dataset]
Gough L, Grace JB (1998) Effects of fooding, salinity and herbivory on coastal plant communities, Louisiana, United States
Guerrero-Ramírez NR, Mommer L, Freschet GT, Iversen CM, McCormack ML, Kattge J et al (2021) Global root traits (GRooT) database. Glob Ecol Biogeogr 30:25–37. https://doi.org/10.1111/geb.13179
Hacala A, Lafage D, Prinzing A, Sawtschuk J, Pétillon J (2021) Drivers of taxonomic, functional and phylogenetic diversities in dominant ground-dwelling arthropods of coastal heathlands. Community Ecol 197:511–522
Hacala A, Le Roy M, Sawtschuk J, Pétillon J (2020) Comparative responses of spiders and plants to maritime heathland restoration. Biodivers Conserv 29:229–249. https://doi.org/10.1007/s10531-019-01880-y
Hacala A, Salgueiro-Simon M, Bain A, Marguerie D, Pétillon J (2023) Spider and vascular plant assemblages in subarctic peat bogs are complementary ecological indicators of variation in local and landscape factors. Ecol Ind 158:111389. https://doi.org/10.1016/j.ecolind.2023.111389
Hambäck P, Dawson L, Geranmayeh (Kynkäänniemi) P, Jarsjö J, Kačergytė I, Peacock M, Collentine D, Destouni G, Futter M, Hugelius G, Hedman S, Jonsson S, Klatt BK, Lindström A, Nilsson J, Pärt T, Schneider L, Strand J, Cordero P, Blicharska M (2022) Tradeoffs and synergies in wetland multifunctionality: A scaling issue. Sci Total Environ 862:160746. https://doi.org/10.1016/j.scitotenv.2022.160746
Hänggi A, Stockli E, Nentwig W (1995) Habitats of central European spiders. Miscellanea Faunistica Helvetiae 4. Neuchatel
Harvey PR, Nellist DR, Telfer MG (2002) Provisional atlas of British spiders (Arachnida, Araneae), Volumes 1 & 2. Biological Records Centre, Huntingdon
Heckbert S, Costanza R, Poloczanska ES, et, Richardson AJ (2011) Climate Regulation as a Service from Estuarine and Coastal Ecosystems. In Treatise on Estuarine and Coastal Science, 199–216. Elsevier. https://doi.org/10.1016/B978-0-12-374711-2.01211-0
Jain M, Flynn DFB, Prager CM, Hart GM, DeVan CM, Ahrestani FS, Palmer MI, Bunker DE, Knops JMH, Jouseau CF, Naeem S (2014) The importance of rare species: a trait-based assessment of rare species contributions to functional diversity and possible ecosystem function in tall‐grass prairies. Ecol Evol 4:104–112. https://doi.org/10.1002/ece3.915
Janousek CN, Folger CL (2014) Variation in Tidal Wetland Plant Diversity and Composition within and among Coastal Estuaries: Assessing the Relative Importance of Environmental Gradients. J Veg Sci 25(2):534–545. https://doi.org/10.1111/jvs.12107
Julve P (1998) Baseflor. Index botanique, écologique et chorologique de la flore de France. Version 31 décembre 2002. http://perso.wanadoo.fr/philippe.julve/catminat.html
Julve P (2018) Listes départementales des plantes de France. Version 2018.04 du 24 avril 2018. Programme chorologie départementale de Tela Botanica. https://www.tela-botanica.org/bdtfx-nn-55763-repartition
Jens K, Bönisch G, Díaz S, Lavorel S, Prentice IC, Leadley P, Tautenhahn S et al (2020) TRY Plant Trait Database – Enhanced Coverage and Open Access. Glob Change Biol 26(1):119–188. https://doi.org/10.1111/gcb.14904
Kennish MJ (2002) Environmental Threats and Environmental Future of Estuaries. Environ Conserv 29(1):78–107. https://doi.org/10.1017/S0376892902000061
Kim D, Ohr S (2020) Coexistence of plant species under harsh environmental conditions: An evaluation of niche differentiation and stochasticity along salt marsh creeks. J Ecol Environ 44(1):19. https://doi.org/10.1186/s41610-020-00161-y
Lafage D, Djoudi EA, Perrin G, Gallet S, Pétillon J (2019) Responses of ground-dwelling spider assemblages to changes in vegetation from wet oligotrophic habitats of Western France. Arthropod-Plant Interact 13:653–662. https://doi.org/10.1007/s11829-019-09685-0
Lafage D, Maugenest S, Bouzillé J-B, Pétillon J (2015) Disentangling the influence of local and landscape factors on alpha and beta diversities: opposite response of plants and ground-dwelling arthropods in wet meadows. Ecol Res 30:1025–1035. https://doi.org/10.1007/s11284-015-1304-0
Lafite R, Romaña L-A (2001) A man-altered macrotidal estuary: The Seine estuary (France): Introduction to the special issue. Estuaries 24:939–939. https://doi.org/10.1007/BF02691262
Lawrence PJ, Sullivan MJP, Mossman HL (2022) Restored Saltmarshes Have Low Beta Diversity Due to Limited Topographic Variation, but This Can Be Countered by Management. J Appl Ecol 59(7):1709–1720. https://doi.org/10.1111/1365-2664.14179
Leroy B, Le Viol I, Pétillon J (2014) Complementarity of rarity, specialisation and functional diversity metrics to assess community responses to environmental changes, using an example of spider communities in salt marshes. Ecol Ind 46:351–357. https://doi.org/10.1016/j.ecolind.2014.06.037
Little C (2000) The Biology of Soft Shores and Estuaries. OUP Oxford
Lyons KG, Brigham CA, Traut BH, et, Schwartz MW (2005) Rare Species and Ecosystem Functioning. Conserv Biol 19(4):1019–1024. https://doi.org/10.1111/j.1523-1739.2005.00106.x
Maison de l’Estuaire (2018) 4ème Plan de gestion de la réserve naturelle nationale de l’estuaire de la Seine - Tome I: Diagnostic. http://www.normandie.developpement-durable.gouv.fr/plan-de-gestion-de-la-rnnes-a77.html#sommaire_1
Maison de l’Estuaire (2023) La réserve naturelle nationale de l’estuaire de la Seine. 2023. https://maisondelestuaire.org/reserve.html
McKee KL, Mendelssohn IA (1989) Response of a Freshwater Marsh Plant Community to Increased Salinity and Increased Water Level. Aquat Bot 34(4):301–316. https://doi.org/10.1016/0304-3770(89)90074-0
McLusky DS, Elliott M (2004) The estuarine ecosystem: ecology, threats, and management, 3rd edn. ed. Oxford University Press, Oxford; New York
Meinzer F (2003) Functional Convergence in Plant Responses to the Environment. Oecologia 134(1):1–11. https://doi.org/10.1007/s00442-002-1088-0
Meire P, Ysebaert T, Damme SV, Bergh EVD, Maris T, Struyf E (2005) The Scheldt estuary: a description of a changing ecosystem. Hydrobiologia 540:1–11. https://doi.org/10.1007/s10750-005-0896-8
Mouchet MA, Villéger S, Mason NWH, Mouillot D (2010) Functional diversity measures: an overview of their redundancy and their ability to discriminate community assembly rules: Functional diversity measures. Funct Ecol 24:867–876. https://doi.org/10.1111/j.1365-2435.2010.01695.x
Mouillot D, Loiseau N, Grenié M, Algar AC, Allegra M, Cadotte MW, Casajus N, Denelle P, Guéguen M, Maire A, Maitner B, McGill BJ, McLean M, Mouquet N, Munoz F, Thuiller W, Villéger S, Violle C, Auber A (2021) The dimensionality and structure of species trait spaces. Ecol Lett 24:1988–2009. https://doi.org/10.1111/ele.13778
Mouny P, Dauvin JC, Bessineton C, Elkaim B, Simon S (1998) Biological components from the Seine estuary: first results. In: Amiard J-C, Le Rouzic B, Berthet B, Bertru G (eds) Oceans, Rivers and Lakes: Energy and Substance Transfers at Interfaces. Springer Netherlands, Dordrecht, pp 333–347. https://doi.org/10.1007/978-94-011-5266-2_26
Neupert M, Aubert M, Langlois E (2024) Tidal restriction does not always supplant local conditions shaping estuarine plant communities. Estuarine, Coastal and Shelf Science
Pardo I, Roquet C, Lavergne S, Olesen JM, Gómez D, García MB (2017) Spatial congruence between taxonomic, phylogenetic and functional hotspots: true pattern or methodological artefact? Divers Distrib 23:209–220. https://doi.org/10.1111/ddi.12511
Pavoine S, Gasc A, Bonsall MB, Mason NWH (2013) Correlations between Phylogenetic and Functional Diversity: Mathematical Artefacts or True Ecological and Evolutionary Processes? Édité par Andreas Prinzing. J Veg Sci 24(5):781–793. https://doi.org/10.1111/jvs.12051
Pavoine S, Marcon E, Ricotta C (2016) Equivalent numbers’ for species, phylogenetic or functional diversity in a nested hierarchy of multiple scales. Methods Ecol Evol 7:1152–1163. https://doi.org/10.1111/2041-210X.12591
Pearce JL, Venier LA (2006) The use of ground beetles (Coleoptera: Carabidae) and spiders (Araneae) as bioindicators of sustainable forest management: A review. Ecol Ind 6:780–793. https://doi.org/10.1016/j.ecolind.2005.03.005
Pekár S, Wolff JO, Černecká Ľ, Birkhofer K, Mammola S, Lowe EC, Fukushima CS, Herberstein ME, Kučera A, Buzatto BA, Djoudi EA, Domenech M, Enciso AV, Piñanez Espejo YMG, Febles S, García LF, Gonçalves-Souza T, Isaia M, Lafage D, Líznarová E, Macías-Hernández N, Magalhães I, Malumbres-Olarte J, Michálek O, Michalik P, Michalko R, Milano F, Munévar A, Nentwig W, Nicolosi G, Painting CJ, Pétillon J, Piano E, Privet K, Ramírez MJ, Ramos C, Řezáč M, Ridel A, Růžička V, Santos I, Sentenská L, Walker L, Wierucka K, Zurita GA, Cardoso P (2021) The World Spider Trait database: a centralized global open repository for curated data on spider traits. Database 2021: baab064. https://doi.org/10.1093/database/baab064
Petchey OL, Gaston KJ (2002) Functional diversity (FD), species richness and community composition. Ecol Lett 5:402–411. https://doi.org/10.1046/j.1461-0248.2002.00339.x
Petchey OL, Gaston KJ (2006) Functional diversity: back to basics and looking forward. Ecol Lett 9:741–758. https://doi.org/10.1111/j.1461-0248.2006.00924.x
Pétillon J, Georges A, Canard A, Lefeuvre J-C, Bakker JP, Ysnel F (2008) Influence of abiotic factors on spider and ground beetle communities in different salt-marsh systems. Basic Appl Ecol 9:743–751. https://doi.org/10.1016/j.baae.2007.08.007
Pétillon J, Georges A, Canard A, Ysnel F (2007) Impact of cutting and sheep grazing on ground–active spiders and carabids in intertidal salt marshes (Western France). Anim Biodivers Conserv 30(2):201–209
Pétillon J, McKinley E, Alexander M, Adams JB, Angelini C, Balke T, Griffin JN, Bouma T, Hacker S, He Q, Hensel MJS, Ibáñez C, Macreadie PI, Martino S, Sharps E, Ballinger R, De Battisti D, Beaumont N, Burdon D, Daleo P, D’Alpaos A, Duggan-Edwards M, Garbutt A, Jenkins S, Ladd CJT, Lewis H, Mariotti G, McDermott O, Mills R, Möller I, Nolte S, Pagès JF, Silliman B, Zhang L, Skov MW (2023) Top ten priorities for global saltmarsh restoration, conservation and ecosystem service research. Sci Total Environ 898:165544. https://doi.org/10.1016/j.scitotenv.2023.165544
Pétillon J, Potier S, Carpentier A, Garbutt A (2014) Evaluating the success of managed realignment for the restoration of salt marshes: Lessons from invertebrate communities. Ecol Eng 69:70–75. https://doi.org/10.1016/j.ecoleng.2014.03.085
Pétillon J, Ysnel F, Gleut SL (2003) Responses of spider communities to salinity and flooding in a tidal salt marsh (Mont St.-Michel Bay, France). (Proceedings of the 21st European Colloquium of Arachnology, St.-Petersburg, August 2003). 236–247
Puzin C, Leroy B, Pétillon J (2014) Intra- and inter-specific variation in size and habitus of two sibling spider species (Araneae: Lycosidae): taxonomic and biogeographic insights from sampling across Europe: Morphological variations in sibling wolf spiders. Biol J Linn Soc Lond 113:85–96. https://doi.org/10.1111/bij.12303
Reed DJ (1995) The Response of Coastal Marshes to Sea-Level Rise: Survival or Submergence? Earth Surf Proc Land 20(1):39–48. https://doi.org/10.1002/esp.3290200105
Reserve Naturelle Estuaire de Seine (2023) 4ème Plan de gestion de la réserve naturelle nationale de l'estuaire de la Seine TOME III. OPERATIONS ET ANNEXES Version révisée
Réserves Naturelles de France (2023) Estuaire de la Seine. https://www.reserves-naturelles.org/estuaire-de-la-seine
Ridel A, Lafage D, Devogel P, Lacoue-Labarthe T, Pétillon J (2021) Habitat filtering differentially modulates phylogenetic and functional diversity relationships between predatory arthropods. R Soc open sci 8:202093. https://doi.org/10.1098/rsos.202093
Roberts MJ (1985) The Spiders of Great Britain and Ireland, Volume 1: Atypidae to Theridiosomatidae (1985) Harley Books, Colchester, UK, ; ISBN 978-0-946589-18-0
Roberts MJ (1987) The Spiders of Great Britain and Ireland, Volume 2: Linyphiidae and Check List (1985) Harley Books, Colchester, UK, ; ISBN 978-90-04-61178-8
Santini L, Belmaker J, Costello MJ, Pereira HM, Rossberg AG, Schipper AM, Ceaușu S, Dornelas M, Hilbers JP, Hortal J, Huijbregts MAJ, Navarro LM, Schiffers KH, Visconti P, Rondinini C (2017) Assessing the suitability of diversity metrics to detect biodiversity change. Biol Conserv 213:341–350. https://doi.org/10.1016/j.biocon.2016.08.024
Scharff N, Coddington JA, Griswold CE, Hormiga G, Bjørn PDP, IN A NORTHERN EUROPEAN DECIDUOUS FOREST (2003) WHEN TO QUIT? ESTIMATING SPIDER SPECIES RICHNESS. J Arachnology 31:246–273. https://doi.org/10.1636/0161-8202(2003)031[0246:WTQESS]2.0.CO;2
Simonneau M, Courtial C, Pétillon J (2016) Phenological and meteorological determinants of spider ballooning in an agricultural landscape. CR Biol 339:408–416. https://doi.org/10.1016/j.crvi.2016.06.007
Soueid-Ahmed A, Jardani A, Dupont JP (2017) Rapport Scientifique de la thèse HYDROMAR. AESN, GPMH, Maison de l’Estuaire
Thrush SF, Townsend M, Hewitt JE, Davies K, Lohrer AM, Lundquist C, Cartner K (2014) THE MANY USES AND VALUES OF ESTUARINE ECOSYSTEMS. In book: Ecosystem Services in New Zealand – Condition and TrendsPublisher. Manaaki Whenua PressEditors: J. Dymond
Uetz GW (1991) Habitat structure and spider foraging. In: Bell SS, McCoy ED, Mushinsky HR (eds) Habitat Structure. Population and Community Biology Series, vol 8. Springer, Dordrecht. https://doi.org/10.1007/978-94-011-3076-9_16
Uetz GW (1999) Guild Structure of Spiders in Major Crops. J Arachnology 27:270–280
Weilhoefer CL (2011) A review of indicators of estuarine tidal wetland condition. Ecol Ind 11:514–525. https://doi.org/10.1016/j.ecolind.2010.07.007
Wilson JB, Stubbs WJ (2012) Evidence for assembly rules: Limiting similarity within a saltmarsh. J Ecol 100(1):210–221. https://doi.org/10.1111/j.1365-2745.2011.01891.x
Yuan X, Ma K, Wang D (2012) Partitioning the effects of environmental and spatial heterogeneity on distribution of plant diversity in the Yellow. River Estuary Sci China Life Sci 55:542–550. https://doi.org/10.1007/s11427-012-4338-3

No competing interests reported.

appendixrideletal.docx

Download PDF

Version 1

posted

You are reading this latest preprint version

Effects of environmental filtering on taxonomic and functional diversity patterns : When spiders and plants provide complementary information to water level management in the Seine estuary

Status:

Version 1

Abstract

Figures

Introduction

Material and methods

Study sites

Sampling design

Vegetation surveys

Environmental variables

Functional traits used

Statistical analyses

Habitat characteristics

Taxonomic vs functional diversity patterns and Q order influence

Differences in community composition and species richness

Indicator species

Proportion of halophilic species

Results

Habitat characteristics

Taxonomic vs functional diversities patterns and Q order influence

Composition

Indicator species

Discussion

Q order influence

Composition of assemblages

Indicator Species

Proportion of halophilic species

Conclusion

Declarations

Author Contribution

Acknowledgement

References

Additional Declarations

Supplementary Files

Status:

Version 1