Molluscs Community as a Keystone Group for Assessing the Impact of Urban Sprawl at Intertidal Ecosystems

Mollusc communities are getting endangered as a result of urban sprawl because articial structures do not surrogate natural substrates. In this study, we compared diversity and community and trophic structures of molluscs among different models of articial substrate and its adjacent natural rock, to detect relationships between some abiotic variables and the molluscs community. Complexity, chemical composition and age were tested as potential drivers of the community. Diversity, community and trophic structure differed between natural and articial substrates. Complexity at the scale of cm was detected as the most important factor driving community structure. Besides, a chemical composition based on silica and/or scarce calcium carbonates seems to be relevant for molluscs, as well as for the secondary substrate where they inhabit. However, age did not seem to be a driving factor. Among the different articial structures, macroscale complexity was detected as the main factor diverging a drastically poor community at seawalls from other articial structures.


Introduction
Biodiversity on natural coastal habitats is under threat by many causes, mainly: coastal arti cialization, exploitation of renewable ( sheries) and non-renewable (mineral and energy extraction) resources, pollutant discharge and marine debris (Dulvy et al. 2003;Jacob et al. 2018).
In the Bay of Sydney (Australia), 50% of the natural coastline is replaced by arti cial substrates (Champan, 2006;Dafforn, 2015) and around 22.000 Km² of European coasts are covered with concrete or asphalt (EEA 2006;Airoldi and Beck, 2007). This coastal transformation, the so-called 'urban sprawl' (Firth et al. 2016), is being boosted by shore erosion due to more frequent stormy events and the sea-level rise (Bouma et al. 2014;Bulleri and Chapman, 2010), altogether, threatening intertidal ecosystems.
Intertidal communities are diverse and complex due to the broad range of biotic and abiotic interactions that occur on intertidal natural substrates (Chapman, 2013).
For example, wave and tide action (Southward and Orton, 1954), desiccation or top-down processes (predation, competition, grazing, etc.) modulate both the sessile and vagile biota, promoting the development of rich and ecologically important communities.
Molluscs are one of the most abundant taxa in the intertidal zone, providing important ecosystem services (see Table 2 in Firth et al. 2016). They are considered early colonizers of substrates (Underwood and Chapman, 2013) and play important roles in C and Si cycles (Meysman and Montserrat, 2017).
Sessile lter molluscs can contribute to clean water and improve nutrient uptake for algae (Eriksson et al. 2017) and together with other sessile organisms, they serve as engineers (Simboura et al. 2017) setting up a secondary substrate for many different species. Also, mobile grazers can feed on macrophytes, cleaning areas for subsequent colonization of many species (Firth et al. 2016). Previous studies have reported a negative impact of arti cial substrates on intertidal molluscs. For example, Moreira et al. (2006) suggested that seawall do not sustain viable populations of limpets. Furthermore, complexity/heterogeneity (e.g. micro-roughness) can affect the abundance of chitons (Moreira et al. 2007) or limpets (Rivera-Ingraham et al. 2011) on arti cial substrates. In fact, substrate complexity is one of the biggest drivers of intertidal biodiversity. Concrete-made arti cial substrates usually lack microhabitats (crevices, rock pools, etc.), preventing refuge from stressful conditions such as desiccation or predation, and are largely responsible for the biodiversity de cit of arti cial substrates compared with the natural rocky shore (Firth et al. 2016 and references therein). Therefore, eco-engineering actions that added habitat complexity at different scales have been related with a higher number of taxa (Ke et al. 2015;Strain et al. 2018) and enhanced recruitment and survival of sessile and mobile macrofauna (Atilla and Finelli, 2005).
Substrate composition, like minerals and elements, is also known to be an important factor affecting communities developing on arti cial structures (Coombes et al. 2015;Sempere-Valverde et al. 2018). For example, acidic siliceous quartz from sandstone may cause oxidative stress and hold less diverse and mature community when compared to limestone (Bavestrello et al. 2000;Cattaneo-vietti et al. 2005). The mineralogical composition usually varies from arti cial substrates (normally made from concrete) to natural substrates (Ido and Shimrit, 2015;Ponti et al. 2015). Concrete may liberate toxic metals and carbonates that enhance alkalinity (pH~13) producing stress on individuals (Ido and Shimrit, 2015). In the case of molluscs, higher saturation of aragonite can facilitate a higher occurrence of burrowing bivalves (Mos et al. 2019) and alkaline concrete surfaces may increase oysters' recruitment (Anderson, 1996).
Although ecological succession may not occur in a parallel manner on arti cial and natural substrates , the age of substrates has been considered as an important factor explaining the differences between arti cial and natural substrates (Glasby and Connel, 1999). Some authors have estimated that it takes from 5 to 20 years for arti cial structures to reach climax communities (Coombes, 2011;Hawkins et al. 1983;Pinn et al. 2005), while others suggest that communities on low crested structures never reach climax (Gacia et al. 2007) or take more than 100 years (Perkol-Finkel et al. 2005).
Consequently, we decided to consider the date of substrates deployment in the present study, together with substrate composition and complexity, to study mollusc diversity associated with arti cial substrates.
Furthermore, changes in the community structure of epifaunal organisms associated with arti cial substrates can cause trophic shifts (Sedano et al. 2020a). Arti cial substrates are known to affect prey resources (Munsch et al. 2015), limiting the diet of some mollusc species (Burgos-Rubio et al. 2015) and ultimately restricting the diversity of trophic strategies. For example, the reduced primary productivity on seawalls has been related to the scarcity of herbivore grazers (Lai et al. 2018). These effects, among others, call for an ecological evaluation of coastal arti cial substrates to prevent the decline of intertidal habitats (Dafforn et al. 2015;Firth et al. 2016) and promote other ecological services (García-Gómez et al. 2014;Dearborn and Kark, 2010). Taking into account that molluscs are diverse and contribute highly to this habitat (Ricciardi et al. 1997), we decided to study the community of molluscs as a model to detect relationships between the abiotic features of the man-made intertidal substrate and the associated fauna. We focused on habitat complexity, substrates composition and age intending to identify which factors are driving the differences in molluscs taxonomic and trophic structure between arti cial substrates and natural substrates. In this regard, we hypothesized that: 1. Substrate complexity and composition would be the main drivers differentiating arti cial from natural substrates, given the differences in complexity and composition between arti cial and natural substrates in our study area. Additionally, we hypothesized that the mollusc community at rip-raps (an arti cial substrate made from natural rock) would be the most similar to natural substrates.
2. Trophic community structure would vary among different arti cial substrates and between arti cial and natural substrates.
3. Age will be a driver structuring intertidal molluscs' community on arti cial substrates.

Study area
Our study area was located in the Algeciras Bay (Cadiz, Spain), which achieves 400 m in depth and occupies 73 Km² of area. This deep bay is found next to one of the most relevant marine regions in the world, the Strait of Gibraltar. It is a marine area with high biodiversity due to its location and structure, which is placed between Africa and Europe and between two water bodies, the Atlantic Ocean and the Mediterranean Sea (Usero et al. 2016). Algeciras Bay contains ve different substrates (four arti cial and their nearest natural rocky shore), very close to each other and under very similar environmental conditions. We selected four nearby arti cial substrates (acropods, cubes, rip-raps and seawall) and compared the molluscan assemblages and trophic structure among them and with the nearest natural substrate. Given the di culty to nd different arti cial substrates next to each other, we limited our study area to this single Bay (Fig. 1).

Abiotic analysis
To identify possible drivers of the differences between substrates, we measured the physicochemical features of each substrate. The variables included macro and microscale complexity, elemental composition, minerals, crystallinity, calcination percentage (C.P) and age. Complexity measures were divided into macroscale complexity (m) and microscale complexity (cm). In both cases (macro and microscale), substrate roughness was calculated as in Rivera-Ingraham et al. (2011) using the equation by Blanchard and Bourget (1999): Roughness or topographical heterogeneity index (THI) = Tr / Ts, where Tr is the "effective" distance between two points "A-B" (measuring the contour between A-B) and Ts is the linear distance between A-B. Macroscale roughness was calculated over 15 m length transects. Three transects were selected at each substrate and a exible meter was laid directly over it, trying to conform as closely as possible to all contours of the bare substrate. Regarding microscale roughness, three 15 cm pro le gauges with 0.5 mm pins were pushed onto the bare rock to record the surface of each substrate (Frost et al. 2005). The resulting pro les were photographed, and the images were digitally processed with Adobe Photoshop to obtain two coloured images. The length of the contour of the pro le was obtained with ImageJ software. The elemental composition and calcination percentage, mineralogical absorption spectra, crystallinity and lithology composition of each sample were obtained from Sedano et al. 2019. All chemical composition was characterized using three powdered fragments of each substrate.
Age of the substrate was based on the date of construction and resulting from the di culty of dating age of the natural substrate, the oldest possible date in the same order of magnitude compared to the oldest arti cial substrate was used instead. Also, wave exposure was quanti ed at each substrate using a combination of the maximum fetch and the modi ed effective fetch (Fe) index developed by (Howes et al. 1994): Fe = [∑ (cos i ) X F i /∑cos i ], where i is the angle between the shore-normal, and the directions 0º, 45º left and 45º right and F i is the fetch distance in Km along the relevant vector. To determine if substrates differed physico-chemically and to detect the most relevant abiotic components that separate the substrates, we performed a Principal Component Analyses (PCA) using macro, microscale complexity, elemental composition (calcium, silicon and magnesium), crystallinity and age. Data were normalized before analyses.

Biotic analysis
Community and trophic structure, as well as biodiversity indices (richness, Pielou's Evenness and Shannon's diversity), were compared among arti cial and natural substrates. Three different sites were randomly selected within each of the ve substrates (natural, cubes, acropods, rip-raps and seawall). At each site, three replicate quadrats of 20 x 20 cm were scraped (3 sites x 3 replicates x 5 substrates = 45 samples). The samples were collected during low tide and within the lower intertidal zone (5-30 cm over the lowest tidal level). We scraped the biotic substrate (secondary substrate) and the associated fauna and preserved it in 96% ethanol until laboratory analyses. At the laboratory, associated molluscs were sorted out from the rest of sessile and vagile biota, identi ed down to species level whenever possible and quanti ed in terms of their abundance. Since the secondary substrate (sessile biota developing on the hard primary substrate) can in uence the associated fauna (Chapman et al. 2005), all sessile fauna and ora that conform the secondary substrate were volumetrically quanti ed at each replicated site to control this variable (used as a covariate in the analyses). Besides, percentages of the most abundant species of the secondary substrate were recorded as well, to detect possible differences between substrates.
To identify possible trophic shifts, we grouped the different species into trophic categories and compared the trophic structure among substrates. Species were assigned and grouped according to their trophic strategies (see Table 1 in Donnarumma et al. 2018), with slight modi cations to better represent the feeding strategy of the species in our study. We assigned them into different groups depending on what they feed on, and the way to obtain the food (trophic guilds) (Table 1). Despite being known to feed on larvae of animals and detritus (Burgos-Rubio et al. 2015), limpets were considered herbivores. We also considered omnivores all animals that feed on suspended organic-matter or by ltering particles in the water column.
After species identi cation, three biodiversity indexes were measured (Badalamenti et al. 2002): species richness, Shannon-Wiener diversity and Pielou evenness for each replicated site in all substrates. To test for biodiversity differences among substrates, we performed a nested ANOVA for each biodiversity index using GMAV5 (Underwood et al. 2002). A Student Newman Keuls (SNK) test was also conducted to elucidate differences among pairs of substrates. The statistic design had two factors: "Substrate" ( xed), with ve levels: Natural, Cubes, Acropods, Rip-raps and Seawall, and "Sites" (random), which was nested in the substrate and had three levels (1,2,3). Cochran's test was performed to con rm the homoscedasticity of biotic data.
Metric multidimensional scaling (MDS) based on the similarity of Bray Curtis matrix was performed on the community and trophic structure (see Table 1 with six categories). We made an additional CLUSTER and SIMPROF test to group the substrates depending on their dissimilarity. Furthermore, a permutational multivariate analysis of variance (PERMANOVA) was also carried out to test if the community and trophic structure of molluscs varied signi cantly among substrates and sites. Previously, data were square-root transformed and the analysis was carried out on a Bray Curtis triangular matrix. When a signi cant source of variance was detected, a pair-wise test between pairs of substrates was also computed to obtain the corresponding p-values.

Correlation analyses.
Correlation tests between the community and abiotic matrix were made to explore potential relationships between community structure and the abiotic variables. Multicollinearity among abiotic factors was previously tested with a Draftsman's plot based on Pearson correlations and only one abiotic factor was used when there were high pairs of correlation (Pearson correlation limit was set a 0.80) (See Figure 5 in Sedano et al. 2019). Variance in ance factor (VIF) was also perfomed to avoid multicollinearity. Furthermore, a distance-based redundancy analysis (dbRDA) was computed using a fourth root transformed biotic matrix paired with a normalized abiotic matrix, to give similar weight to variables measured on different units. DbRDA was portrayed into a bidimensional representation. A BIOENV routine (Clarke and Ainsworth, 1993) was done to detect the best set of variables that better suit the response data. This method calculates correlation coe cients between response variables (community matrix) and predictor variables (abiotic matrix) (Balkenhol et al. 2009). RELATE routine (Clarke and Warwick, 2001) was carried out to detect the correlation coe cient between both, community and abiotic dissimilarity matrix. All multivariant and correlation analyses were carried out with PRIMER + PERMANOVA 6 using 9999 permutations (Anderson et al. 2008) Results

Abiotic analyses
The results of the fetch index indicate that all the substrates belong to a semi-exposed wave exposure class, whereas the age of origin was different for each substrate (Table 2). Regarding substrate complexity, microscale complexity was higher at natural substrate and cubes comparing with the rest of substrates, being very low at seawall and acropods. In contrast, macroscale complexity was higher at acropods and cubes than rip-raps and natural substrate (Table 2). From a chemical point of view, elemental composition differed between natural and arti cial substrates and among all substrates. Silica (SiO2) concentration was higher at natural substrate compared to arti cial substrates, which were characterized by a higher concentration of calcium oxide (CaO) at all samples (Table 3).
According to the mineralogical composition, natural substrate was very different from the arti cial ones, and differences were also found within the arti cial substrates. The natural substrate was composed of high percentages of quartz, while cubes and rip-raps were mostly composed by quartz and CaO in a carbonated form, calcite (CaCO3). Acropods presented high levels of magnesium oxide (MgO) and their mineralogical composition was based on dolomite (CaMg(CO 3 ) 2 ) (  Age causes the separation of cubes and natural samples from the rest of substrate samples, correlating with micro-scale complexity.
Regarding trophic groups, natural substrate contained all the groups measured, and scavenger was exclusive for this substrate. This group was formed by a single species, Tritia tingitana, and its abundance was 2. On the other hand, mostly all arti cial substrates lacked the group detritus feeder except for seawall. The differences in the percentages of groups between substrates were also remarkable.
In terms of percentages, more macro and micro grazers appeared at arti cial substrates compared to natural substrate. In contrast, predators and detritus feeders appeared in higher percentage at natural substrate, with the exception of seawall, where the percentage of predators and detritus feeders was higher than that of natural substrate.
There were also differences between arti cial substrates because the percentage of lter feeders was higher at rip-raps and cubes compared to natural substrate and acropods, while more macro and micro grazers appeared at acropods compared to the rest of substrates. Finally, seawall had a very heterogeneous trophic structure among samples ( Figure 5, supplementary material Table 1).
Shannon's diversity and richness varied signi cantly among substrates (p<0,001) (Table 4). According to SNK test, Shannon's diversity was greater on the natural substrate compared to the arti cial substrate. Among arti cial substrates, acropods showed higher Shannon's diversity values than seawall, cubes and rip-raps. In contrast, cubes, rip-raps and seawall did not differ signi cantly on Shannon's diversity.
MDS analysis for the community structure showed three groups: 1 = natural, 2 = seawall and 3 = acropods, rip-raps and cubes. These groups were statistically supported by the SIMPROF test (p<0.05).
The natural group appeared homogeneous and different from the rest of the substrates. Seawall group was heterogeneous, but it also appeared clearly segregated. The third group, formed by acropods, cubes and rip-raps was homogeneous but distinct from natural substrate and seawall groups (Figure 7).
MDS for the trophic structure revealed three groups: 1 = seawall, 2 = seawall and 3= natural, acropods, cubes, rip-raps and seawall. These groups were statistically supported by the SIMPROF test (p<0.05). Seawall was the most heterogeneous substrate diverging into three groups, while acropods, rip-raps, natural substrate and cubes were more similar. However natural substrate was homogeneous and was signi cantly different from the rest at a level of similarity of 80% (Figure 8).
PERMANOVA test indicated signi cant differences in community structure between substrates and among sites. According to the pair-wise tests, the community at natural substrate differed from the rest. When comparing among arti cial substrates, seawall differed from rip-raps and acropods (p<0,01), but not from cubes. Cubes, acropods and rip-raps seemed to have a similar community structure (Table 5).
PERMANOVA test indicated differences in trophic structure among substrates but not among sites. The pair-wise test revealed that trophic structure on the natural substrate was different from the arti cial substrates. Among arti cial substrates, seawall differed from rip-raps (<0,05) and acropods (p<0,01) but not from cubes. Also, acropods, rip-raps and cubes had a similar trophic community ( Table 6). The volume of the secondary substrate using as a covariable was signi cant for taxonomical PERMANOVA but not for the trophic one.

Correlation analyses:
A dbRDA analysis revealed a relationship between the physical-chemical composition and the molluscan community. The associated community of molluscs at natural substrate was highly correlated with a low CaO concentration and carbonated nature, a high microscale complexity and older age. On the other hand, the community at arti cial substrates was correlated with a high CaO concentration and carbonated nature, low microscale complexity and younger age. Because SiO2 appeared negatively collineated with CaO and C.P, these parameters were not included in the analyses. Among arti cial substrates, the seawall community was also correlated with a low macroscale complexity, and therefore they have clustered apart from the more heterogeneous arti cial substrates (acropods, cubes and ripraps) (Figure 9). The BIOENV analyses showed that the most correlated variables (p<0.01) were macro, microscale complexity and crystallinity (Rho=0,669). The RELATE test showed a signi cant correlation of abiotic with taxonomic matrices (Rho = 0,518; p<0,01).

Discussion
The molluscs community structure and diversity seemed to be signi cantly different between arti cial and natural substrates. Among our studied variables, substrate complexity (macro and microscale roughness) and chemical composition appeared to be the main drivers of those differences. Besides, the trophic structure also seemed to be different between arti cial and natural substrates.
By our results, habitat complexity in terms of the relative abundance of microhabitats such as crevices (Evans et al. 2016), rockpools and macrophytes is considered as one of the most in uencing factors on intertidal communities (Warfe et al. 2008). Higher heterogeneity at the scale of centimeters increases recruitment of spores and larvae (Sempere-Valverde et al. 2018) due to a higher number of refugees (Kostylev, 2005;Coombes et al. 2010). This can be particularly important for intertidal molluscs since they can not only nd shelter against environmental stress (Meager, 2011;Harley and Helmuth, 2003;Loke et al. 2015), but also against predation (Warfe et al. 2008) and competition (Huston, 1979) by nding crevices that t their shell size (Loke et al. 2016), determining community structure and diversity.
Moreover, algae turfs that cover the rocky substrate can in uence the abundance and biodiversity of associated fauna, playing biogenic roles, similar to sessile animals such as barnacles or annelids who play the role of "ecological engineers", providing the secondary substrate where many species live (Simboura et al. 1995;Bavestrello et al. 2000). The greater abundance of the calcareous algae Ellisolandia elongata on natural substrate can in uence the associated fauna via increasing both the habitat volume and habitat complexity (Guerra and Garcia et al. 2012;Veiga et al. 2014;Torres et al. 2015). Moreover, it can decrease desiccation and temperature stress by providing shelter for mobile fauna (Singh et al. 2013;Ke et al. 2015).
In our study, 15 species were exclusive of natural substrate. Natural substrate had high microscale complexity, but they were also highly covered by the calcareous algae Ellisolandia elongata, altogether possibly boosting the higher occurrence of more taxa.  (Avila, 2003). Similarly, more abundance of sea snails and bivalves appeared at natural substrate. For example, Cerithiopsis tubercularis is usually restricted to live on algae that are associated with its food (sponges), as it is the case of the branched Ellisolandia spp. and its association with the sponges Halichondria and Hymeniacidon (Fretter and Manly, 1977). Given the close association between algae and molluscs that certain species can present, the absence or very low abundances of these species at arti cial substrates, where the cover of algae was very scarce, highlights the importance of calcareous algae in supporting richer communities of molluscs on arti cial substrates at this area. In contrast, a lower microscale complexity and scarce algae canopy, probably lead to a lower abundance of the less competitive bivalves and sea snails because fewer microhabitats are available Fairweather 1989, Hills et al. 1996;Strain et al.2018) as it happens at seawall. Also, the sandstone porosity of natural substrates increases algae settlement (Green et al. 2012), probably generating positive cascading effects.
However, a higher complexity at the scale of meters increases recruitment of propagules and the dissipation of water energy (Vieira et al. 2020) on cubes, acropods and rip-raps, boosting the abundance, richness and diversity of associated fauna, on these substrates, in comparison with seawall. However, these species were mostly "limpets" as Fissurela nubeluca, Siphonaria pectinata and Patella caerulea and chitons. The increment of this taxa is possibly related to the fact that arti cial substrates are better habitat for sedentary species, such as limpets and chitons, rather than strictly vagile gastropods (Rivera-Ingraham et al. 2011;Cha et al. 2013), probably by suffering lower predation and being more resilient to wave action. In fact, non-native species of Siphonaria and barnacles have been recorded on seawalls at Plymouth and Singapore (Hsiung et al. 2020) Seawalls have a small intertidal area for recruitment but, as it happens in our study, seawalls harbour abundant beds of mussels (Chapman et al. 2005) and barnacles on the secondary substrate, associated with lower biodiversity values in comparison with natural substrates (People, 2006;Sedano et al. 2020b). The community at seawall was very scarce and had the lowest diversity. Chapman, 2006 suggested that seawalls lack microhabitats for many species and limit the life strategies of specialized intertidal fauna, such as limpets and chitons. For example, as it has been recorded in the chiton belonging to Ischiochiton genera that inhabits underneath the boulders as habitat-specialist (Grayson and Chapman, 2004). On the other hand, the pulmonate limpet Siphonaria pectinata was absent at seawall, in accordance with Moreira et al. 2006 who detected a relation among living on seawall and a reduction on the reproductive output of this limpet. However, these results contrast with Hsiung et al. 2020 who recently detected non-native species of Siphonaria guanemensis and barnacles on seawalls at Plymouth and Singapore.
Chemical composition was also identi ed as a possible driver of the community, mainly differentiating communities settled on natural or arti cial substrates, since the natural rock was mainly pure quartz (SiO2), while arti cial substrates had great amount of carbonated minerals, with high levels of the calcite (CaCO3). The effect of quartz on natural substrates and the carbonated mineralogy at arti cial substrates could affect the associated fauna. For example, it has been reported that quarzitic radicals inhibit the settlement of rst recruits of secondary substrates such as the hydroid Eudendrium glomeratum (Bavestrello et al. 2000) or the sponge Clionia sp. (Cerranoet al.2007) while they are neutral for algae settlement. Moreover, in the present study, associated fauna was more abundant and more diverse at natural substrate, where more Ellisolandia elongata appeared, possibly due to a reduction in competition with other sessile biota affected by the toxicity of silicon radicals (Cerrano et al. 1999). In addition, facilitation by calcium hydroxides that are liberated by concrete arti cial substrates to the substrate surface, alkalinizing the pH, also contributes to the settlement of bivalves (Anderson, 1996;Soniat and Burton, 2005;Burt et al. 2009) and barnacles (Guilbeau et al. 2003) on the sessile substrate, as it occurs in the concrete substrates in this study (acropods and seawall).
The concrete substrates are also rich in magnesium oxides and other minerals, which could in uence the presence of exclusive species. For example, aragonite has been related with improving boring bivalves' settlement (Green et al. 2013), being more soluble on water than calcite (Cornelis and Cornelius, 2007), and acropods that are composed by this material showed the presence of the boring species Leisonelus aristatus.
Therefore, a combination of a carbonated nature and a lower microscale complexity at arti cial substrates possibly promotes a different community of molluscs and increase the dominance of the most 'colonizer' species of the secondary substrate (mussels and barnacles) at cubes, acropods and seawall (Miller and Etter, 2008;Underwood and Chapman, 2013), all disturbing the associated fauna of molluscs.
According to the trophic structure, natural substrate seems to be more diverse, mostly because they contain all the trophic groups measured, while acropods, rip-raps and cubes lack of suspension feeder and scavengers. Contrarily, seawall had all the measured groups except scavengers. Another interesting difference was the different percentage of groups among substrates because detritus feeder and predator were more abundant at natural substrate than at acropods, rip-raps and cubes, where more percentage of macro and micro-grazer appeared. Natural substrate hold the majority of detritus feeders, a fact that could be related to higher sediment retention by macrophytes (Melero et al. 2017;Casoli et al. 2019). In contrast, the detritus feeder Barleeia unifasciata appeared at seawall. In fact, littorinid snails have been related to breakwaters, with lower crevice availability (Aguilera et al. 2014). In addition, predators were exclusive from natural substrate, probably because intraguild predation (Janssen et al. 2007) and the number of preys has been reported to be lower at structures with less complexity. In contrast, macro and micro grazers were highly abundant at arti cial substrates, mostly derived from an increment in species of limpets and chitons (see the rst part of Discussion).
Limpets are known to control the volume of macrophytes in concert with sea-urchins by their grazing activity (Piazzi et al. 2016) and, in this area, they have been recorded as omnivorous and very generalist (Burgos-Rubio et al. 2015). This could explain the lower volume of the secondary substrate at cubes, rip-raps and acropods when compared to natural substrate and seawall. In addition, the higher volume of secondary substrate on seawall, a substrate with a low abundance of molluscs, sustains the hypothesis that grazers control these sessile populations at acropods, rip-raps and cubes, in special the Ellisolandia elongata, as has been observed at rip-raps. The idea that these grazers could be controlling the associated fauna at arti cial substrates should be considered.
Several authors had pointed out that biodiversity could be driven by the age of substrate (Perkol-Finkel et al. 2005;Glasby and Connel, 1999), and others have reported that temporal heterogeneity among arti cial and natural substrates is a relevant factor driving communities (Glasby and Connel, 1999). For example, on arti cial substrates, rst recruits as ephemeral algae, sponges and bivalves occur fast in less than a year on but later, as a consequence of a low microscale complexity, dominant species outcompete the rst colonizers . Nevertheless, among the arti cial substrates studied in the present work, age did not appear as a driver of the community, because the community at cubes (80y) and rip-raps (20y) was similar, and community of seawall (20y) and acropods (20y) differed, independently of age.

Conclusions And Future Approach
As expected in our rst hypothesis, mollusc community and diversity differed between arti cial and natural substrates. Distortion in the bottom-up interactions between, a combination of low microscale complexity and a carbonated nature, rich in calcium, of the arti cial substrates and the mollusc community, seems to impact over many species of molluscs and the common calcareous algae Ellisolandia elongata they inhabit in, comparing with the natural substrate. Moreover, macroscale complexity seems to in uence the community of molluscs, increasing recruitment of species at acropods, rip-raps and cubes in comparison with seawall, but mostly bene ting limpets, chitons and bivalves, and also barnacles on the secondary substrate.
As for our second hypothesis, physico-chemical factors seem to alter the trophic community increasing the percentage of macro and micro grazers and lter feeders on arti cial substrates. In contrast with our third hypothesis, age did not appear as a driver of the mollusc community.
We suggest that, according to previous studies, increasing habitat heterogeneity by means of increasing crevices (Archambault and Bourget 1996) and rock-pools (Ostalé-Valriberas et al. 2018) and microscale complexity is fundamental in the future designs of arti cial substrates. On the other hand, the chemical structure should be included as an important topic of research in new models of arti cial substrates, possibly depending on the geology and chemistry of the surrounding lands.
In the case of Algeciras, rip-raps and cubes were the most similar to natural substrate in relation to abiotic features. At the same time, in relation to community and trophic structure, these substrates seem to be less disturbed. Even though the study was replicated within substrates, the lack of several types of substrates under similar environmental factors on a higher spatial scale represents a limitation. In this sense, further studies re-analyzing already published data and/or meta-analysis at higher spatial scales will provide valuable insights on the role of arti cial substrates in structuring coastal assemblages.
Furthermore, more research should clarify how molluscs or macrophytes recruitment and survivance are in uenced by chemical and other physical issues of the substrate. Finally, biological interactions among secondary substrates and associated fauna should be also explored in future designs.

ACKNOWLEDGEMENTS
We want to thank our colleagues Rafael Espinar for helping in data collection, Carlos Navarro for helping with data analyses, Jose Manuel Guerra for identi cation assistance and Marta Florido for her tips with the writing. We also express our gratitude to the anonymous reviewers that have contributed to improve the early version of the manuscript. FUNDING X-raY analyses were funded by Ministerio de Economía y Competitividad, Spain (Project CGL2017-82739-P co-nanced by Agencia Estatal de Investigacion-AEI-and Fondo Europeo de Desarrollo Regional -FEDER-).

CONFLICTS OF INTERESTS/ETHICS APPROVAL/ CONSENT TO PARTICIPATE
"Not applicable to this study."

CONSENT FOR PUBLICATION
The author gives Urban Ecosystems consent for publication.

Figure 2
Principal Component Analyses (PCA) portraying the differences among substrates according to the abiotic data measured.
Page 27/33 Figure 3 Total volume of secondary substrate (ml) in the Y axis and substrates in the X axis. Percentages of volume of the most abundant species from secondary substrate (species ml/total species ml) at each substrate.

Figure 4
Venn diagram with the percentage of shared species of associated fauna of molluscs among the ve substrates.