1. Introduction
Posidonia oceanica (L) Delile is a Mediterranean endemic seagrass and, in optimal conditions, is considered the main habitat-forming species of one of the climax communities of some Mediterranean infralittoral areas due to the high productivity, complexity and stability of this vascular plant [1,2]. Meadows of this seagrass have an estimated coverage of 25–50% of infralittoral bottoms of the Mediterranean Sea (down to 50 m in certain areas), corresponding to the most extensive seagrass meadows for this basin [2,3]. These meadows support a highly diverse associated biotic community [1], providing an elevated number of ecosystem services such as the maintenance of marine biodiversity (including species of commercial interest), buffering of strong water movement leading to protection of the coastline against erosion, nutrient cycling and carbon sequestration, among others [2,4,5,6,7]. Despite all these services, P. oceanica has suffered a strong decline over the last centuries [7,8,9], with an estimated regression of meadows at 34% in the last 50 years as a consequence of the cumulative effects of multiple local stressors [10].
The decline and negative impacts suffered by this ecologically important species have led the EC to include P. oceanica meadows as priority habitat Posidonia beds (Posidonion oceanicae) (Habitat 1120) for conservation purposes within the EU Habitats Directive (92/43/EC). In addition, it is an indicator species and habitat for evaluating the “Good Environmental Status” of European waters in line with the Marine Strategy Framework Directive (MSFD, 2008/56/EC). Therefore, it is necessary to study the phenological parameters and the spatio-temporal variability of P. oceanica meadows across its biogeographical distribution, including the boundaries of its geographical distribution, as a preliminary step for further monitoring programs in order to improve the knowledge and management of this important marine ecosystem element.
Posidonia oceanica has its westernmost distribution limit close to the Strait of Gibraltar (SoG) in the northern Alboran Sea, being generally absent along the Moroccan coastline except for the Chafarinas Islands (35°05′ N, 02°25′ E) [10,11]. At its easternmost limit, P. oceanica occurs in a few localities along the southern Levantine part of Turkey close to the Kizilliman marine protected area and Alexandria (Egypt) on the southern coast [10,12,13]. In the Alboran Sea, P. oceanica forms extensive beds in its easternmost sector, close to the Almeria-Oran Front. In the central and western sectors, these beds are reduced from usually fragmented meadows to small patches [5,11]. Despite the fragmentation experienced by P. oceanica meadows close to the SoG, several authors have documented highly diverse associated invertebrate assemblages (e.g., molluscs [14] and crustacean decapods [15]). This fosters the idea that P. oceanica patches have an important ecological and conservation value whatever their size, as even small patches comprising a low number of shoots still support a higher abundance and diversity of invertebrates than adjacent sedimentary bottoms, as observed in other seagrass species [14,15].
The Alboran Sea displays transitional characteristics between the Atlantic Ocean and the Mediterranean Sea, with a predominance of colder and less salty water masses than other Mediterranean areas [16]. Moreover, the entrance of surface Atlantic waters through the SoG generates a series of anticyclonic gyres due to several physical factors (northeast orientation of the SoG, the Coriolis effect, the topography of the Alboran basin, the local climatology) that produce upwellings of cold and nutrient enriched Mediterranean deep waters in the northwestern Alboran coasts [17,18] (Figure 1). This sector of the Mediterranean Sea supports one of the highest productive areas within the region [19,20], which, together with the complex seafloor heterogeneity and the geographical location of the Alboran basin, promotes a high diversity of habitats and species [21,22], including ones of commercial importance. The latter has allowed the development of an intensive fishing activity targeting demersal, small pelagic and benthic species [23,24,25,26], with potentially adverse effects on vegetated benthic assemblages [27,28].
Meadows of P. oceanica have been widely studied in the Mediterranean Sea, with most studies being focused on seasonal and monthly variations [28,29,30,31,32,33], whereas few of them have considered interannual differences [34,35,36]. Despite the ecological importance of P. oceanica, studies on phenological aspects are scarce at the limits of the seagrass distribution range, including the Alboran Sea, where hydrological characteristics of the basin are very different from those in other parts of the Mediterranean Sea due to the Atlantic influence. Within this context, peripheral or edge populations are important for a species’ long-term survival and evolution [37,38], as these populations usually present a high genetic differentiation compared to central populations due to reduced population size and genetic drift in fragmented habitats, which results in a lower genetic diversity [39]. Ecotypes derived from this genetic differentiation respond differently to environmental changes, being able to show even better tolerance to thermal shock than central populations [40]; however, their response depends on local population traits, as has been observed in some marine macrophytes (e.g., Laminaria digitata, Fucus vesiculosus, Zostera marina, P. oceanica) [40,41,42]. This highlights the importance of studying and conserving edge populations of marine macrophytes.
The main objectives of this study are (1) to characterize P. oceanica meadows at its westernmost gradient of distribution in the Mediterranean Sea; (2) to analyze the seasonal and interannual variation of different phenological parameters in the Mediterranean–Atlantic transition region in comparison to other Mediterranean areas; and (3) to evaluate flowering events in different P. oceanica meadows of the northwestern Alboran coasts.
2. Materials and Methods
2.1. Study Area
Characterization of phenological parameters of P. oceanica meadows, together with the assessment of flowering events, was carried out at four sites along the northwestern Alboran margin (Figure 1). These were, from east to west: the marine protected area ‘‘Paraje Natural Acantilados de Maro-Cerro Gordo” (hereafter Maro MPA) (Natura 2000 site code: ES6170002), the Special Area of Conservation “Calahonda” (hereafter Calahonda SAC) (Natura 2000 site code: ES6170030), the Special Area of Conservation “Fondos Marinos de la Bahía de Estepona” (hereafter Estepona SAC) (Natura 2000 site code: ES6170036) and Punta Chullera (hereafter Chullera) located just 40 km from the SoG.
Maro MPA is located between Nerja (Málaga) and Almuñécar (Granada) (Figure 1) and includes a narrow coastal strip (12 km) with a series of small beaches, soft and rocky bottoms, coves and cliffs, underwater seawalls and submerged caves [43]. This site covers 15.29 km2 of the MPA with a depth range from the shore down to ca. 70 m. The studied P. oceanica meadow is located in a small cove with mixed pebble, rocky and sandy bottoms called “Molino de Papel” at a depth between 2 m (on rocky bottoms) and 12 m (on sandy bottoms). Calahonda SAC is located ca. 10 km further west between Fuengirola and Marbella, covering an area of 14 km2 from the shore down to ca. 30 m depth (Figure 1). This site is characterized by the presence of different soft bottoms and rocky outcrops in the infralittoral zone, as well as intertidal sandy and rocky shores [43]. The studied meadows correspond to patches of variable sizes located on shallow rocky bottoms at 2–4 m depth. Further west, Estepona SAC is located between the harbour of Estepona and Punta de la Plata (Figure 1). This 6 km2 MPA is separated by a distance of ca. 150 m from the shore and covers a depth range from 1 to ca. 50 m [43]. The studied meadow occurs on shallow rocky bottoms at depths not exceeding 4 m. Finally, the westernmost site was Chullera, a rocky cape with several coves and low cliffs located along the coast of Manilva and very close to the SoG and Atlantic Ocean (Figure 1). Soft bottoms predominate at Chullera, with rocky outcrops appearing in the shallowest area [43]. Here, P. oceanica occurs as very small patches on rocky bottoms at a depth of around 1–2 m. More information about these MPAs and on the general distribution maps of P. oceanica meadows can be found in Mateo-Ramírez et al. [43] and Ruiz et al. [44], respectively. Sea-surface temperature and salinity, although with temporal variations, generally increase from Chullera to Maro MPA due to the greater influence of Atlantic waters (less salty and warm) on the western sites, with values ranging from 16.9 °C to 18.3 °C and from 36.7 to 37.2 psu, respectively [45,46,47]. Wave action and tidal variation in the northwestern Alboran Sea are somehow influenced by its proximity to the Atlantic Ocean. For example, in Calahonda SAC, maximum wave height values are reported in autumn–winter (1–1.5 m) and, due to its proximity to the SoG, the frequent westerly and easterly winds (velocities ranging between 3.6–5.7 m·s−1) favour constant wave action during the whole year (data taken from
2.2. Data Collection
The general characterization of the selected P. oceanica meadows was done using SCUBA diving during different years from 2006 to 2010 and 2015 when information about bathymetric range, substrate types where P. oceanica occurs, patch size, coverage, shoot density and other phenological parameters (e.g., leaf height, number of leaves per shoot, number of inflorescences per square meter), as well as environmental variables was collected.
The line intercept method was used for cover estimates and included the use of eight 25 m transects in each sampling station. Cover estimates corresponded to the percentage of the substrate covered by live P. oceanica in relation to the whole surface area indicated as “% coverage”. On the other hand, a non-destructive technique was used for collecting shoot density and phenological data among the studied meadows within the 2–7 m depth range (depth of maximum coverage in most of the sites) from mid-November to mid-December of 2015. At each station, patches of different sizes were sampled, including three small patches (patch size <1 m2), three medium patches (1–5 m2) and three large patches (>5 m2). Three quadrats of 50 × 50 cm were deployed for shoot density measurements (number of shoots per m2) in each of these patches (total = 27 quadrats per station). In the case of small patches, quadrats were deployed in directly adjacent patches as the patch size was too small for three replicates. The sampling area of the quadrats is similar to that used by other authors in previous studies on shallow Mediterranean seagrass beds [31,33,49,50] and is also used for seagrass monitoring [51]. In addition, 10 haphazardly selected shoots were analyzed in situ in each quadrat for leaf parameters (total = 270 shoots per station) following the methodology suggested by Gobert [35], such as the number of leaves per shoot (considering adult, intermediate and juvenile leaves) and the height (based on leaf length) of the highest (longest) leaf, which was measured to the nearest mm from the basal part of the sheath to the blade tip using a meter rule. Finally, comparative flowering data (number of flowers in the quadrats, i.e., flowering density) of P. oceanica was recorded across the northwestern Alboran Sea (Maro MPA, Calahonda SAC, Estepona SAC and Chullera) in November–December 2015. Flowering frequency was calculated as the percentage of flowering shoots in relation to the total number of shoots within the quadrats, according to Pergent et al. [30].
Seasonal and interannual variation of phenological parameters of P. oceanica was studied in greater detail at Calahonda SAC as it represents a site with a good representation of meadows located in an intermediate area within the seagrass westernmost geographical boundary. Moreover, this site is still highly influenced by surface Atlantic waters, which promotes the occurrence of a large number of algae species with Atlantic affinity [52,53] (Figure 1). Phenological and environmental data were recorded at two sampling stations within this site, “Punta de Calaburras” (hereafter Calaburras; 36°30.4′ N–04°38.3′ W) and “Calahonda” (36°29.2′ N–04°42.04′ W), the latter located 7 km to the west and ca. 100 km to the SoG. Shoot density and leaf parameters (the same mentioned above) were recorded seasonally from July 2007 to March 2010 in patches with sizes ranging between 1–130 m2 at a depth of 2–3 m. Shoot density was estimated from 5 quadrats (50 × 50 cm) at each sampling station.
Assessment of leaf parameters at Calahonda SAC included quantification of the number of leaves per shoot counted in 10 shoots randomly selected in situ within each quadrat, together with measurements of the highest leaf height in 50 shoots per season and sampling station. When inflorescences were observed, flowering data (flowering density and frequency) were collected for studying the seasonal and interannual variation of flowering events. Additionally, the area of those patches containing inflorescences was measured.
In order to explore the relationships of seasonal and interannual changes of seagrass phenological parameters with those of environmental parameters, surface seawater temperature (SST), solar irradiance and concentration of chlorophyll a (Chl a) were measured 1–2 weeks before, during and after plant phenology data collection at Calaburras and Calahonda. SST was measured at midday (12:00) with an alcohol thermometer. Two replicates of 1 litre of surface seawater were collected per day at the site and transported in darkness at low temperatures to the laboratory for Chl a determination. Pigment analyses were carried out by filtering through Whatman GF/C glass filters, and extraction made using 100% acetone for 12h in cool and dark conditions. Measurements were made using a spectrophotometer at wavelengths of 630, 647, 664 and 750 nm. Concentrations of Chl a were obtained using the equation proposed by Jeffrey and Humphrey [54]. Solar irradiance measurements were obtained using HOBO Data Loggers (model UA-002-08, Onset) that were deployed seasonally from winter 2009 to winter 2010 (an annual cycle). Each HOBO was placed daily close to the studied seagrass meadows at 4–5 m depth and set to record data with 10 min intervals between 10:00 and 14:00. Additionally, the organic matter content (% OM) in sediment was also included in the analyses and measured from 5 replicates of sediment samples collected within the studied P. oceanica meadows or adjacent to patches in each season. The % OM in samples of dry sediment (80 g) was obtained from the weight loss after ignition at 560 °C for 1 h.
Finally, the biotic index based on P. oceanica (BiPo) [55] was calculated in order to evaluate the conservation status of the different P. oceanica meadows studied and documented here. The BiPo index was developed to evaluate the ecological status of this seagrass, and it is in accordance with the EU Water Framework Directive requirements. The quantification of the ecological status is made according to the degree of deviation from reference conditions and expressed as a numerical value ranging from 1 (best conditions) to 0 (worst conditions), known as the ecological quality ratio (EQR) (see Lopez y Royo et al. [55] for more details).
2.3. Statistical Analyses
Differences in shoot density and phenological parameters (number of leaves, leaf height) were analyzed using a two-way PERMANOVA design 1 with two fixed factors: site (4 levels: Maro MPA, Calahonda SAC, Estepona SAC, Chullera); and patch size (3 levels: small, medium, large). Seasonal and interannual differences in P. oceanica dynamics were analyzed using a two-way PERMANOVA design 2 with two fixed factors: seasons (4 levels: Summer [July–September], Autumn [November–December], Winter [January–March], Spring [May–June]); and years (3 levels: A1 [July 2007–May 2008], A2 [September 2008–June 2009], A3 [September 2009–March 2010]). The two-way PERMANOVA design 2 was also used to test for statistical differences in environmental variables between seasons and years. Differences in inflorescence density could only be studied in Calahonda SAC meadows using a one-way PERMANOVA design 3 with one fixed factor: patch size (3 levels: small, medium, large). Analyses were based on Euclidean distances [56], and the significance of p values was determined through 9999 permutations of raw data and residuals under a reduced model for one-way and two-way analysis, respectively. When the numbers of permutations were low, the Monte Carlo method was used to contrast the significance of p values. Two-way PERMANOVA pair-wise post hoc tests (p < 0.05) were used for posterior multiple comparisons across months and years in design 2. In order to know the magnitude of temporal fluctuations in shoot density and phenological parameters, the coefficient of variation (CV, %) was calculated with all collected data.
When seagrass parameters, sediment and environmental variables displayed a similar pattern, forward stepwise multiple linear regression analyses were used to investigate the corresponding relationships. Solar irradiance was analyzed independently because data was only available from March 2009 to March 2010. Prior to carrying out the regression analyses, potential significant relationships between the independent variables (seagrass parameters, sediment and water column variables) were investigated using Pearson’s correlation coefficient, and those that were significantly and highly correlated were eliminated from the analyses. Data used in the regression analyses were checked for normality using Kolmogorov and Smirnov test and transformed to natural logarithms when necessary. These statistical procedures were performed using the software SPSS 20 and PRIMER 6.0, and PERMANOVA statistical package.
3. Results
3.1. Spatial Variability along the Northwestern Alboran Sea
The bathymetric distribution of P. oceanica meadows within each site displayed a decreasing trend towards the Strait of Gibraltar (SoG), with those of the Maro MPA (easternmost station) displaying a wide depth range (2–12 m depth); those of Calahonda SAC displaying an intermediate depth range (1–5 m depth); and those of Estepona SAC and Chullera (westernmost stations) displaying a very narrow depth range (1–2.5 m and 1–1.5 m depth, respectively).
The studied P. oceanica meadows were generally fringing or growing on rocky outcrops, displaying a gradient in the meadow structure from continuous meadows to small patches towards the SoG. The meadows at Maro MPA presented small patches (0.25–0.5 m2) at shallow depths (2–3 m depth) and medium to large patches (up to 8 m2) at a depth of 6 m depth, where it generally changed to a more or less continuous bed (Figure 1). Meadows of P. oceanica at Calahonda SAC and Estepona SAC were structured as medium to large patches, sometimes becoming fragmented meadows (Figure 1). Finally, the meadows at Chullera were reduced to a few small to medium patches (0.5–1.5 m2) located on very shallow rocky outcrops (Figure 1).
Coverage of P. oceanica within the 1 to 6 m depth range displayed a decreasing trend towards the SoG, ranging from ca. 50% in Maro MPA to ca. 45% in Calahonda SAC and Estepona SAC and less than 5% in Chullera with a low number of patches of small size. In Maro MPA, P. oceanica coverage also decreased with depth, showing values lower than 10% at its lower limit (12 m depth). Dead matte heavily covered by macroalgae (mainly Halopteris scoparia) was detected between patches of living matte in most study sites (Figure 2), especially at Calahonda SAC, resulting in isolated and/or fragmented seagrass populations.
An increasing trend of shoot density was detected towards the SoG, with the largest shoot density mean values detected at Chullera (753 ± 63 shoots m–2; mean ± SE) and Estepona SAC (662 ± 27 shoots m–2), whereas similar values were reported at Maro MPA (550 ± 44 shoots m−2) and Calahonda SAC (514 ± 42 shoots m–2) (Figure 3A). Mean maximum leaf height was highest at Chullera (25.8 ± 3.8 cm) and Estepona SCI (25.3 ± 1.2 cm), followed by Calahonda SAC (22.1 ± 1.4 cm) (Figure 3B; Table 1) and Maro MPA (17.6 ± 0.9 cm) (Figure 3B). The mean number of leaves per shoot showed a significant decreasing trend from Maro MPA (6.6 ± 0.3 leaves shoot−1) to Chullera (5.9 ± 0.4 leaves shoot−1), with intermediate values at Calahonda SAC and Estepona SAC (6.1 ± 0.2 and 6.3 ± 0.3 leaves shoot−1), respectively (Figure 3C; Table 1).
Regarding patch size, and despite non-significant differences detected in shoot density with patch size (Table 1), the pairwise analysis showed statistical differences for the interaction term of site and patch size (Figure 3A). The highest shoot density values were recorded in small-size patches of Chullera (823 ± 27 shoots m−2), whereas the lowest values were recorded in small-size patches of Maro MPA (406 ± 20 shoots m−2) (Figure 3A). On the other hand, the highest leaf height values observed at the time of the study (November–December) were generally detected in medium-size patches at all sites, except those of Maro MPA, with the highest mean values observed at Chullera (31.2 ± 1.9 cm). Meadows at Maro MPA showed similar leaf height independent of patch size and had the lowest value from all the studied sites (17.7 ± 1.9 cm) (Figure 3B; Table 1). Finally, the number of leaves per shoot was significantly higher in small-size and medium-size patches, ranging between 6.7 ± 0.7 and 6.1 ± 1 leaves shoot−1 (Figure 3C; Table 1).
3.2. Flowering
Three flowering events were detected in the P. oceanica meadows of Calahonda SAC during the six years of monitoring (2006–2010 and 2015). Mature inflorescences were observed at Calaburras in March 2009 on orthotropic shoots of several patches with areas ranging between 4 and 10 m2 at a depth of ca. 2 m, with a mean value of 12.7 ± 4.1 inflorescences per m2 (infl. m−2) that represented a 1.9% flowering frequency (flow. freq.), with a maximum density of 28 infl. m−2. A second and major flowering event was observed in November 2009 in almost all patches studied in Calahonda SAC, with mean values of 34.0 ± 7.4 infl. m−2 at Calaburras (7.0 ± 3.0% flow. freq.) and of 69.0 ± 15.3 infl. m−2 at Calahonda (6.7 ± 3.4% flow. freq.), and with maximum densities of 144 and 184 infl. m−2, respectively. Finally, a massive flowering event was detected at Calahonda in November 2015, with a mean value of 115.2 ± 98.2 infl. m−2 (22.3 ± 6.6% flow. freq.) and a maximum density of 330 infl. m−2 (Figure 2A). Here, PERMANOVA analyses did not reveal significant differences in inflorescence density among patch sizes (Table 1). On this last occasion, no inflorescences were detected at Maro MPA, Estepona SAC and Chullera.
3.3. Seasonal and Interannual Variability
Mean surface seawater temperatures (SST) showed a clear temporal trend at both Calahonda and Calaburras sampling stations, with values ranging from 22 °C in warm months (with exceptional values of 27 °C in September 2009) to 15 °C in cold months (Figure 4A). Seasonal differences in SST were significant at both sampling stations, whereas annual differences were only significant at Calaburras (Table 2). The Chl a concentration displayed a similar temporal trend, with the highest values in warm months and the lowest ones in cold months, with Calahonda generally displaying higher values (Figure 4B). Regarding Chl a, seasonal differences were significant at both sampling stations, with significant annual differences only at Calahonda (Table 2). Solar irradiance displayed a temporal trend with high values in warm months and low ones in cold months (Figure 4C), but no significant seasonal differences were detected (Table 2). The percentage of organic matter in sediments collected within the seagrass patches showed a similar seasonal trend at both sampling stations, with high values in warm months (ca. 2.4%) and low ones in cold months (ca. 1.6%) (Table 2).
The mean shoot density of P. oceanica was not significantly different amongst years (annual cycles) and sampling stations of Calahonda SAC (one-way PERMANOVA; F = 0.031, p = 0.862), with very similar values of 850 ± 15 shoots m−2 at Calaburras and 845 ± 20 shoots m−2 at Calahonda (Table 3). Seasonal (monthly) variability was only significant at Calahonda, where shoot density displayed a temporal trend with higher values in cold months (maximum: 1055 shoots m−2, November 2007) during the first two years (Figure 5A; Table 3), whereas at Calaburras shoot density was very similar throughout the year (Figure 5A; Table 3, Supplementary Table S1). Finally, shoot density displayed negative correlations with the number of leaves (RPearson = −0.387, p < 0.01), as well as with maximum leaf height (RPearson = −0.268, p < 0.01).
Maximum leaf height displayed similar values in both sampling stations, with mean values ranging between 26.6 ± 1.4 cm at Calahonda and 24.9 ± 1.2 cm at Calaburras (one-way PERMANOVA: F = 0.863, p = 0.357). A temporal trend with significantly higher leaf height in cold months (maximum 39 cm in both sites, May 2008) was observed in both sampling stations (Figure 5B; Table 3; Supplementary Table S1), as well as significant interannual differences with mean annual values ranging between 22.5 ± 1.7 and 27.6 ± 2.5 cm (Figure 5B; Table 3; Supplementary Table S1).
The mean number of leaves per shoot was similar at both sampling stations (Calaburras: 5.7 ± 0.1 leaves shoot−1; Calahonda: 5.6 ± 0.1 leaves shoot−1) (one-way PERMANOVA; F = 0.422, p = 0.519), with no significant differences amongst annual cycles. On the other hand, the number of leaves displayed a temporal trend with significant seasonal (monthly) differences and maximum values in cold months (maxima: 6.3 ± 0.1 leaves shoot−1 at Calahonda, March 2009; 6.1 ± 0.1 leaves shoot−1 at Calaburras, January 2008) (Figure 5C; Table 3; Supplementary Table S1).
3.4. Relationships between Phenological, Sediment and Environmental Variables
This analysis could be performed in the seasonal and intra-annual study carried out at Calahonda SAC. Here, results for the irradiance showed the highest values of R squared and explained 47.9% and 59.4% of the variability in relation to the number of leaves and shoot density, respectively (Table 4).
3.5. Status
Most of the studied meadows presented a good conservation status during the sampling period (2007–2015), as reflected in the EQR value of the BiPo index ranging between 0.57 and 0.64 (Table 5). The exception was the meadow studied at Maro MPA, where a lower EQR value indicates a moderate conservation status (Table 5). The majority of meadows presented a sharp lower limit, which was located on shallow bottoms (5 m depth or less), and elevated values of shoot density and leaf height (Table 5). In the case of Maro MPA, meadows showed a sparse and deeper lower limit with lower values of shoot density and leaf height (Table 5).
4. Discussion
The distribution of Posidonia oceanica along the northern Alboran Sea, where its western limit of distribution is located, is not uniform. Here, the most extensive and continuous beds occur off El Ejido (Almería) at the Site of Community Importance, “Fondos Marinos de Punta Entinas-Sabinar” (code ES6110009) [11]. From here, towards the Strait of Gibraltar (SoG) and the Atlantic Ocean, P. oceanica meadows are becoming increasingly fragmented, and their extent is considerably reduced. The P. oceanica meadows studied along the northwestern Alboran Sea display a decreasing trend of patch size, bathymetric range and the number of leaves per shoot on moving toward the SoG, as well as an increasing trend of shoot density and leaf height. Similarly, the maximum depth reached by P. oceanica meadows decreases westwards from Punta Entinas-Sabinar (ca. 15 m depth) [57] to Chullera (1–1.5 m depth). The influence of surface Atlantic water entering eastwards through the SoG, together with the characteristic oceanographic and hydrological conditions of the northern Alboran Sea (e.g., higher turbidity, stronger wave action) [5,44,58], might represent a developmental limiting factor for P. oceanica in this part of the Mediterranean. In this respect, other studies have documented the presence of other habitat types at unusually shallow bottoms in the northwestern Alboran Sea, such as coralligenous assemblages found at a depth of 15 m, which is usually characteristic of the circalittoral biocenosis [59,60]. This indicates lower water transparency in the northwestern Alboran Sea compared to other Mediterranean areas, as has also been detected in the adjacent Gulf of Cádiz. Moreover, the analyzed P. oceanica meadows are located in one of the highest biological productivity areas within the Mediterranean [61,62], which may also induce low water transparency. This was reflected in the low irradiance values obtained at a depth of 4–5 m at the Calahonda SAC, which are similar to those found at a depth of 10 m in the eastern Adriatic Sea [50].
Solar irradiance and seawater temperature are key factors modulating the abundance and productivity of seagrasses [8,12,63]. Indeed, light and surface seawater temperature explained a high percentage of variability of leaf number and shoot density of P. oceanica at the Calahonda SAC throughout the year. On the other hand, a narrower bathymetric distribution (ranging between 1–5 m depth) exposes the meadows to continuous wave action and turbulence, especially in the areas close to the SoG, where patches were very close to the surf zone (1–1.5 m depth). Regarding this aspect, shallow patchy meadows of different seagrasses have been reported to be affected by such extreme conditions (e.g., low water transparency and strong wave action), including P. oceanica meadows at its northern distributional boundary in the Adriatic Sea [64,65,66,67]. Within this context, studies carried out on P. oceanica, and other seagrasses show a negative relationship between seagrass biomass, leaf height, epiphyte biomass, leaf width and shoot area with wave energy [68,69]. Moreover, high-energy environments support P. oceanica meadows landscapes with higher patch numbers, lower coverage, more complex patch shapes and lower intra-patch architecture than sites exposed to lower wave action [69].
Shoot densities of the studied P. oceanica meadows were generally higher than those located in other Mediterranean areas at similar depths, including the north Adriatic Sea [66]) and the northeastern Levantine Sea (eastern distributional boundary [12]), but similar to those reported in meadows in Tunisia growing on a rocky substratum, as the ones studied here [70] (see Figure 6). Regarding leaf height, values recorded in the present study were lower than those documented in other Mediterranean areas (see Figure 6). Regarding this, Marbà et al. [8] reported P. oceanica meadows with high shoot densities and low leaf height values in other northern Alboran Sea locations [8]. Several physical factors acting together seem to highly influence the westernmost P. oceanica meadows, such as the shallow location and the rocky nature of the substratum, together with higher hydrodynamic exposure due to wave action. Furthermore, they may represent key drivers structuring ecological features and phenological parameters of these meadows along the northern Alboran Sea [11,68,69], being particularly evident at the westernmost meadows of Estepona SAC and especially of Chullera. This could represent a stress reduction strategy that involves increasing rhizome growth for anchorage as well as self-shading effects [35,71,72], as suggested by the negative correlations observed at Calahonda SAC between shoot densities, leaf height and the number of leaves per shoot.
The temporal trend observed for leaf height and number of leaves per shoot is similar to that reported for other P. oceanica meadows of both European and African coasts [30,31,32,36,73,74] Regarding shoot density and despite being a more stable descriptor, higher values were observed in cold months in Calahonda SAC. Moreover, an unusual seasonal trend was observed during the first two years, similar to that previously documented in other meadows, some of which were located close to the northern and eastern distribution edges of P. oceanica in Tunisia and Turkey [36,71,75]. Nevertheless, shoot density and number of leaves per shoot displayed lower intra-annual fluctuations when compared to other P. oceanica meadows (see Table 6), probably linked to the particular environmental characteristics of the northwestern Alboran Sea, including lower intra-annual temperature changes (15.3–22.5 °C) in comparison to other Mediterranean areas (e.g., Adriatic Sea: 9–24 °C; Tunisian coasts: 14–31 °C; southern Italian coasts: 11.8–24.8 °C [33,73,76]), which may highly influence these meadows due to their shallow location. On the other hand, maximum leaf height displayed high differences amongst years, as observed previously [8,35,36]. Interannual studies of phenological parameters are scarce, but all of them have indicated some interannual variability, highlighting the need for performing long-term studies of meadow population characters and phenological parameters [8,35,36].
The inflorescences followed by fruit development observed at Calahonda SAC would indicate that sexual reproduction of P. oceanica may potentially occur in meadows at their western limit of geographical distribution [80]. Flowering of P. oceanica is highly variable within and between meadows, with differences even at small spatial scales of a few meters [81]. Despite this, the flowering and recruitment of this Mediterranean seagrass seem to be more frequent than previously expected [82,83]. The available information seems to indicate that flowering represents an irregular event mainly related to high seawater temperature [84,85]. In fact, three flowering events were registered in the meadows at Calahonda SAC between 2006 and 2015, with the highest flowering episode coinciding with the warmest seawater temperature values ever recorded for the area at that time. Regarding this, it is remarkable that the high flowering frequency observed in 2015, with similar or even higher values than those reported in meadows of the central and eastern Mediterranean basin [75,86,87]. Nevertheless, no flowering events were detected in Maro MPA, Estepona SAC, or Chullera during that year. This is of interest as Calahonda SAC meadows are located close to the westernmost distributional limit of P. oceanica, which would indicate (1) that these meadows present an adequate health status for sexual reproduction; or (2) a survival mechanism in response to natural or anthropogenic disturbances [85,88]. Furthermore, the intensity of the flowering and the absence of significant differences amongst patches of different sizes in relation to inflorescence density may be related to the latter notion. Fruit production and seed germination still need to be monitored in these P. oceanica meadows in long-term studies in order to evaluate the extent of sexual reproduction in the formation of new patches. This is of importance to evaluate the survival capacity of P. oceanica, as edge populations could play an important role in the long-term survival of this seagrass; in fact, a recent study shows that thermal edge populations of P. oceanica have a greater resilience against thermal stress than central populations [42].
Disentangling natural, anthropogenic impacts on different organisms and habitats is still very difficult to achieve. In spite of being exposed to high anthropogenic pressure (intense coastal development and high human population), most of the studied P. oceanica meadows were in a good conservation status as reflected in the BiPo index, with EQR values similar to those of other meadows located in areas that are more pristine and where the human population density is low [55]. This could be related to the peculiar characteristics of the P. oceanica meadows at our study sites compared to other sites in the Mediterranean Sea (see Figure 6). Nevertheless, the Maro MPA meadow showed evidence of decline, which could be partly related to large flooding linked to a torrential rainfall (200 mm per m2) recorded in the area in September 2007. This had serious consequences for the P. oceanica meadow located here as several streams flow into the small cove of “Molino de Papel”. Flooding of these watercourses led to the transport of sand and mud, plant debris, and waste of different kinds into the sea, which were deposited in a large part of the meadow. This caused the burial and subsequent death of approximately half the existing meadow [89] (see Figure 2D). The Calahonda SAC meadows also show a high presence of old and dead rhizomes, with no leaves and the matte covered by macroalgae, which may indicate that P. oceanica has suffered a regression and it was more abundant in the past (see Figure 2C).
5. Conclusions
The P. oceanica meadows of the northwestern Alboran Sea represent the westernmost located ones ever studied. These meadows have an important ecological role in the area, supporting high-diversity communities [4,5,14,15,43,44,59] and showing some unique characteristics (e.g., higher shoot density and lower leaf height) that differentiate them from other Mediterranean meadows present at similar depths. Nevertheless, the seasonal dynamics of phenological parameters were similar to those recorded in other areas that fall within the biogeographical distribution of P. oceanica. Although the conservation status of the studied meadows is good and flowering events have been reported, the present study highlights recent impacts on some of them linked to different anthropogenic pressures occurring at the studied coastal area of southern Spain (e.g., coastal infrastructures, elevated urbanism, insufficient wastewater treatment plants and illegal fishing). For these reasons, it is recommended to carry out detailed long-term monitoring of these meadows in order to improve the management of the westernmost populations of this EU-priority habitat. Further monitoring of these seagrass meadows as part of the Marine Strategy Framework Directive may provide crucial information on the improvement or decline in the state of these singular meadows located in a very different environmental situation compared to other Mediterranean ones. Furthermore, studies on the potential restoration of P. oceanica, and other seagrasses, will contribute to balancing the alarming and continuous regression of meadows detected in the northern Alboran Sea during the last decades [90].
Á.M.-R., P.M., J.L.R. and J.U. contributed to the conceptualization and design of the study. Á.M.-R., P.M., J.L.R., E.B.-E. and J.U. contributed to data acquisition. Á.M.-R., P.M., A.M.-A., E.B.-E., J.E.G.R., J.L.R. and J.U. contributed to writing, editing, and reviewing the manuscript. J.L.R. and J.U. supervised and coordinated the study. All authors have read and agreed to the published version of the manuscript.
Not applicable.
We would like to express our sincere gratitude to the “Consejería de Medio Ambiente de la Junta de Andalucía” and RNM-0141 Research Group from the Universidad de Málaga for their institutional support, the staff of the “Centro de Buceo Benalmádena”, for their help and continuous interest, and to Terence W. Edwards from Sunny View School for the English revision of earlier drafts of this manuscript.
The authors declare no conflict of interest.
Footnotes
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Enlarge this image.Figure 1. Map of the study area in the northwestern Alboran Sea and images of the general configuration of Posidonia oceanica patches in the studied meadows. (A) General map. Orange areas indicate the location of the different sites with P. oceanica meadows sampled: (B) Maro MPA, (C) Calahonda SAC, (D) Estepona SAC and (E) Chullera. Green lines indicate areas influenced by nutrient-rich upwellings; yellow arrows indicate the direction of the Surface Atlantic Waters and pink arrows indicate the direction of the Mediterranean Intermediate and Deep. Redrawn from http://www.juntadeandalucia.es/medioambiente/rediam (accessed on: 14 April 2022).
Enlarge this image.Figure 2. (A) Posidonia oceanica patch with several inflorescences (arrows) and (B) detail of one inflorescence from Calahonda SAC meadows, November 2015. (C) Dead matte covered by macroalgae (mainly Halopteris scoparia) from Calahonda SAC meadows. (D) Dead matte from Maro MPA meadows after burial due to sediment transport linked to a torrential rainfall in September 2007.
Enlarge this image.Figure 3. Shoot density and phenological parameters of Posidonia oceanica as a function of patch sizes and sites. (A) Shoot density (shoots m–2), (B) maximum leaf height (cm) and (C) number of leaves per shoot (leaves shoot–1). Different letters indicate significant differences within each patch size category (p < 0.05). Mean + standard error.
Enlarge this image.Figure 4. Temporal changes of environmental variables in Posidonia oceanica meadows in Calaburras (empty bars) and Calahonda (solid bars) sampling stations within the Calahonda SAC in different months/seasons throughout the first (A1), second (A2) and third (A3) annual cycles: (A) Surface seawater temperature, (B) concentration of Chlorophyll a and (C) solar irradiance in different seasons. July, Jl; November, Nv; January, Ja; May, My; September, Sp; December, Dc; March, Mr; June, Ju. Different letters indicate significant differences between months, with uppercase letters for Calaburras and lowercase letters for Calahonda (p < 0.05). Mean + standard error.
Enlarge this image.Figure 5. Temporal changes in shoot density, leaf height and number of leaves of Posidonia oceanica meadows at Calaburras empty rhombus and Calahonda solid square within the Calahonda SAC, in different months/seasons throughout the first (A1), second (A2) and third (A3) annual cycles: (A) Shoot density (shoots m−2), (B) maximum leaf height (cm) and (C) number of leaves shoot−1). July, Jl; November, Nv; January, Ja; May, My; September, Sp; December, Dc; March, Mr; June, Ju. Mean ± standard error.
![View Image - Figure 6. Shoot density (A) and leaf height (B) of Posidonia oceanica at the studied sites of the northwestern Alboran Sea (in the case of Calahonda SAC, data were clustered from different seasons and years considering Calaburras and Calahonda), as well as at other locations of the Mediterranean Sea. The bold black line in the box represents the median value. The upper and lower ends of the box represent the third (Q3) and first (Q1) quartiles, respectively. The ends of the lines extending vertically from the bars (whiskers) are the minimum and maximum values. Outliers are plotted as individual points. CHU, Punta Chullera; ESTE, Estepona SAC; CALA, Calahonda SAC; MARO, Maro MPA; TAB, Tabarca Island (E Spain) [32]; AL, Algeria [77]; FR, France [30,34]; TUS and TUD, Tunisian shallow (<5 m depth) and deep (>5 m) meadows, respectively [73]; IT, Italy [31,33,36]; LI, Libya [78]; MA, Malta [79]; CR, Croatia [50]; GRS and GRD, Greek shallow (<5 m) and deep (>5 m) meadows, respectively [49]; TKS and TKD, Turkish shallow (<5 m) and deep (>5 m) meadows, respectively [12].](https://media.proquest.com/media/hms/THUM/12/l1CzQ?_a=ChgyMDI1MDcyNTEwMzQyNzk4MzoxNDYyNjUSBzEzODMzMjMaCk9ORV9TRUFSQ0giDTIxNi43My4yMTYuNzYqBzUwNDY4NzQyCjI3OTE2NzM5NzE6C0lubGluZUltYWdlQgEwUgZPbmxpbmVaA0lMSWIEVEhVTWoKMjAyMy8wMS8wMXIKMjAyMy8xMi8zMXoAggEgUC0xMDAwMDg0LTIwMTM5NS1TVEFGRi1udWxsLW51bGySAQZPbmxpbmXKAWdNb3ppbGxhLzUuMCBBcHBsZVdlYktpdC81MzcuMzYgKEtIVE1MLCBsaWtlIEdlY2tvOyBjb21wYXRpYmxlOyBDbGF1ZGVCb3QvMS4wOyArY2xhdWRlYm90QGFudGhyb3BpYy5jb20p0gESU2Nob2xhcmx5IEpvdXJuYWxzmgIHUHJlUGFpZKoCG09TOkVNUy1GdWxsVGV4dC1nZXRGdWxsVGV4dMoCD0FydGljbGV8RmVhdHVyZeICAPICAPoCAVmCAwNXZWKKAxxDSUQ6MjAyNTA3MjUxMDM0Mjc5ODM6NDc5MjY0&_s=tJ9Ko8hkajUKdZCc4R2NIxcXL9g%3D "View Image - Figure 6. Shoot density (A) and leaf height (B) of Posidonia oceanica at the studied sites of the northwestern Alboran Sea (in the case of Calahonda SAC, data were clustered from different seasons and years considering Calaburras and Calahonda), as well as at other locations of the Mediterranean Sea. The bold black line in the box represents the median value. The upper and lower ends of the box represent the third (Q3) and first (Q1) quartiles, respectively. The ends of the lines extending vertically from the bars (whiskers) are the minimum and maximum values. Outliers are plotted as individual points. CHU, Punta Chullera; ESTE, Estepona SAC; CALA, Calahonda SAC; MARO, Maro MPA; TAB, Tabarca Island (E Spain) [32]; AL, Algeria [77]; FR, France [30,34]; TUS and TUD, Tunisian shallow (<5 m depth) and deep (>5 m) meadows, respectively [73]; IT, Italy [31,33,36]; LI, Libya [78]; MA, Malta [79]; CR, Croatia [50]; GRS and GRD, Greek shallow (<5 m) and deep (>5 m) meadows, respectively [49]; TKS and TKD, Turkish shallow (<5 m) and deep (>5 m) meadows, respectively [12].")
Enlarge this image.Figure 6. Shoot density (A) and leaf height (B) of Posidonia oceanica at the studied sites of the northwestern Alboran Sea (in the case of Calahonda SAC, data were clustered from different seasons and years considering Calaburras and Calahonda), as well as at other locations of the Mediterranean Sea. The bold black line in the box represents the median value. The upper and lower ends of the box represent the third (Q3) and first (Q1) quartiles, respectively. The ends of the lines extending vertically from the bars (whiskers) are the minimum and maximum values. Outliers are plotted as individual points. CHU, Punta Chullera; ESTE, Estepona SAC; CALA, Calahonda SAC; MARO, Maro MPA; TAB, Tabarca Island (E Spain) [32]; AL, Algeria [77]; FR, France [30,34]; TUS and TUD, Tunisian shallow (<5 m depth) and deep (>5 m) meadows, respectively [73]; IT, Italy [31,33,36]; LI, Libya [78]; MA, Malta [79]; CR, Croatia [50]; GRS and GRD, Greek shallow (<5 m) and deep (>5 m) meadows, respectively [49]; TKS and TKD, Turkish shallow (<5 m) and deep (>5 m) meadows, respectively [12].
Results of two-way PERMANOVA test for differences in shoot density, maximum leaf height, and number of leaves of Posidonia oceanica amongst sites and patch size. Results of one-way PERMANOVA test for differences in inflorescence density of P. oceanica meadows of Calahonda SAC amongst patch sizes. Significant differences are indicated by: ** p < 0.01; *** p < 0.001.
| Source of Variation | df | MS | F | p |
|---|---|---|---|---|
| Shoot density | ||||
| Site | 3 | 77,026 | 5.7217 | ** |
| Patch size | 2 | 9290.4 | 0.69011 | 0.5228 |
| Site × Patch size | 5 | 30,653 | 2.2769 | 0.0767 |
| Residual | 21 | 13,462 | ||
| Total | 31 | |||
| Leaf height | ||||
| Site | 3 | 1215.5 | 66.207 | *** |
| Patch size | 2 | 397.53 | 21.653 | *** |
| Site × Patch size | 5 | 143.97 | 7.8415 | *** |
| Residual | 309 | 18.359 | ||
| Total | 319 | |||
| Number of leaves | ||||
| Site | 3 | 8.7974 | 13.739 | *** |
| Patch size | 2 | 3.1323 | 4.8916 | ** |
| Site × Patch size | 5 | 2.4972 | 38,998 | ** |
| Residual | 309 | 0.64035 | ||
| Total | 319 | |||
| Inflorescences density | ||||
| Patch size | 2 | 15,900 | 2.1032 | 0.2107 |
| Residual | 6 | 7559.8 | ||
| Total | 8 |
Results of one and two-way PERMANOVA test for spatial and temporal changes of environmental variables in Posidonia oceanica meadows at Calaburras and Calahonda (Calahonda SAC, northwestern Alboran Sea) in relation to sampling stations, months from different seasons and years (annual cycles). Significant differences are indicated by: * p < 0.05; ** p < 0.01; *** p < 0.001.
| Source of Variation | Surface Seawater Temperature | Solar Irradiance | Chlorophyll a | Percentage of Organic Matter | |
|---|---|---|---|---|---|
| Calaburras | Year | F = 4.033 * | -- | F = 2.896 | F = 0.898 |
| p = 0.068 | p = 0.414 | ||||
| Season | F = 75.677 *** | F = 2.075 | F = 5.324 ** | F = 6.474 *** | |
| p = 0.164 | |||||
| Year × Season | F = 19.292 *** | -- | F = 4.971 ** | F = 1.315 | |
| p = 0.263 | |||||
| Calahonda | Year | F = 1.172 | -- | F = 6.575 ** | F = 1.822 |
| p = 0.344 | p = 0.171 | ||||
| Season | F = 8.872 ** | F = 1.410 | F = 4.410 * | F = 10.323 *** | |
| p = 0.307 | |||||
| Year × Season | F = 0.646 | -- | F = 3.947 * | F = 2.324 * | |
| p = 0.639 | |||||
| Among sampling stations | F = 0.195 | F = 3.410 | F = 7.136 ** | F = 4.498 * | |
| p = 0.657 | p = 0.078 | ||||
Results of two-way PERMANOVA tests for differences in the shoot density, maximum leaf height and number of leaves per shoot of Posidonia oceanica between seasons and years. Significant differences * p < 0.05; *** p < 0.001.
| Source of Variation | Calaburras | Calahonda | ||||||
|---|---|---|---|---|---|---|---|---|
| df | MS | F | p | df | MS | F | p | |
| Shoot density | ||||||||
| Year | 2 | 12,801 | 1.178 | 0.324 | 2 | 14,726 | 1.616 | 0.206 |
| Season | 3 | 11,404 | 1.05 | 0.382 | 3 | 162,810 | 17.867 | *** |
| Year × Season | 5 | 23,106 | 2.127 | 0.080 | 5 | 50,290 | 5.519 | *** |
| Residual | 43 | 10,865 | 43 | 9112.4 | ||||
| Total | 53 | 53 | ||||||
| Leaf height | ||||||||
| Year | 2 | 136.59 | 8.005 | * | 2 | 214.8 | 20.952 | *** |
| Season | 3 | 1449 | 84.921 | *** | 3 | 1052.8 | 102.7 | *** |
| Year × Season | 5 | 97.508 | 5.715 | *** | 5 | 32.055 | 3.1267 | * |
| Residual | 43 | 17.062 | 43 | 10.252 | ||||
| Total | 53 | 53 | ||||||
| Number of leaves | ||||||||
| Year | 2 | 0.361 | 2.583 | 0.084 | 2 | 0.098 | 0.751 | 0.473 |
| Season | 3 | 1.39 | 9.943 | *** | 3 | 1.713 | 13.176 | *** |
| Year × Season | 5 | 0.232 | 1.663 | 0.167 | 5 | 0.175 | 1.349 | 0.267 |
| Residual | 43 | 0.14 | 43 | 0.13 | ||||
| Total | 53 | 53 | ||||||
Results of regression analyses of Posidonia oceanica for shoot density and number of leaves with irradiance (for only one annual cycle) and surface seawater temperature. Only environmental variables accepted in each regression model (p < 0.05) are listed. * p < 0.05.
| Coefficient | SE | F Ratio | R2 | p | |
|---|---|---|---|---|---|
| Number of leaves | |||||
| Constant | 6.308 | 0.311 | |||
| Temperature | −0.037 | 0.17 | 4.573 | ||
| Number of leaves | 0.479 | * | |||
| Constant | 10.635 | 1.929 | |||
| Irradiance | −1.319 | 0.52 | 6.424 | ||
| Shoot density | 0.594 | * | |||
| Constant | −571.8 | 431.4 | |||
| Irradiance | 372.4 | 116.3 | 10.255 |
Evaluation of the different Posidonia oceanica meadows studied in the northwestern Alboran Sea according to the BiPo index (Biotic Index Posidonia oceanica) [
| Site-Sampling Station (Year) | Lower Limit Depth (m) | Lower Limit Type | Shoot Density (Shoots m−2) | Leaf Height (mm) | EQR | Class | |
|---|---|---|---|---|---|---|---|
| Chullera | (2015) | 2 | Sharp | 753 | 258 | 0.63 | Good |
| Estepona SAC | (2015) | 4 | Sharp | 662 | 253 | 0.61 | Good |
| Calahonda SAC | Calaburrras (2007–2010) | 5 | Sharp | 849 | 250 | 0.61 | Good |
| Calahonda (2007–2010) | 5 | Sharp | 945 | 237 | 0.64 | Good | |
| Calahonda (2015) | 5 | Sharp | 514 | 221 | 0.57 | Good | |
| Maro MPA | (2015) | 12 | Sparse | 550 | 176 | 0.52 | Moderate |
Magnitude of intra-annual fluctuations (coefficient of variation, %) in shoot density and phenological parameters of Posidonia oceanica meadows in Calaburras and Calahonda (NW Alboran Sea), as well as in other locations of the Mediterranean Sea.
| Studied Area | Country | Reference | Shoot Density | Leaf Height | Number of Leaves |
|---|---|---|---|---|---|
| Calaburras | Spain | Present study | 13.1 | 38.7 | 8.5 |
| Calahonda | Spain | Present study | 17.4 | 34.7 | 8.3 |
| Tabarca | Spain | [ |
58.4 | 10.7 | |
| Banyuls-sur-Mer | France | [ |
23.5 | 6.7 | |
| Port-Cros | France | [ |
39.6 | 19.4 | |
| La Revellata Bay (Corsica) | France | [ |
20.6 | 30 | 18.4 |
| Monterosso al Mare | Italy | [ |
42.4 | 43.7 | 20.6 |
| Lacco Ameno (Ischia Island) | Italy | [ |
32.3 | 10.4 | |
| Otranto (Apulia) | Italy | [ |
41.2 | 34.5 | 3.6 |
Supplementary Materials
The following are available online at
References
1. Mazzella, L.; Buia, M.C.; Gambi, M.C.; Lorenti, M.; Russo, G.F.; Scipione, M.B.; Zupo, V. Plant-animal trophic relationships in the Posidonia oceanica ecosystem of the Mediterranean Sea: A review. Plant-Animal Interactions in the Marine Benthos; John, D.M.; Howkins, S.J.; Price, J.H. Clarendon Press: Oxford, UK, 1992; pp. 165-187.
2. Boudouresque, C.F.; Bernard, G.; Bonhomme, P.; Charbonnel, E.; Diviacco, G.; Meinesz, A.; Pergent, G.; Pergent-Martini, C.; Ruitton, S.; Tunesi, L. Protection and Conservation of Posidonia oceanica Meadows; RAMOGE and RAC/SPA: Tunis, Tunisia, 2012; 202. ISBN 2905540311
3. Pasqualini, V.; Pergent-Martini, C.; Clabaut, P.; Pergent, G. Mapping of Posidonia oceanica using Aerial Photographs and Side Scan Sonar: Application off the Island of Corsica (France). Estuar. Coast. Shelf Sci.; 1998; 47, pp. 359-367. [DOI: https://dx.doi.org/10.1006/ecss.1998.0361]
4. Hemminga, M.A.; Duarte, C.M. Seagrass Ecology; Cambridge University Press: Cambridge, UK, 2000; [DOI: https://dx.doi.org/10.1017/CBO9780511525551]
5. Luque, Á.A.; Templado, J. Praderas y Bosques Marinos de Andalucía; Consejería de Medio Ambiente Junta de Andalucía: Sevilla, Spain, 2004; 334. ISBN 9788496329225
6. Pergent-Martini, C.; Pergent, G.; Monnier, B.; Boudouresque, C.-F.; Mori, C.; Valette-Sansevin, A. Contribution of Posidonia oceanica meadows in the context of climate change mitigation in the Mediterranean Sea. Mar. Environ. Res.; 2020; 165, 105236. [DOI: https://dx.doi.org/10.1016/j.marenvres.2020.105236] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33360312]
7. Duarte, C.M. The future of seagrass meadows. Environ. Conserv.; 2002; 29, pp. 192-206. [DOI: https://dx.doi.org/10.1017/S0376892902000127]
8. Marbà, N.; Duarte, C.; Cebrián, J.; Gallegos, M.; Olesen, B.; Sand-Jensen, K. Growth and population dynamics of Posidonia oceanica on the Spanish Mediterranean coast: Elucidating seagrass decline. Mar. Ecol. Prog. Ser.; 1996; 137, pp. 203-213. [DOI: https://dx.doi.org/10.3354/meps137203]
9. Marbà, N.; Duarte, C.M. Mediterranean warming triggers seagrass (Posidonia oceanica) shoot mortality. Glob. Chang. Biol.; 2009; 16, pp. 2366-2375. [DOI: https://dx.doi.org/10.1111/j.1365-2486.2009.02130.x]
10. Telesca, L.; Belluscio, A.; Criscoli, A.; Ardizzone, G.; Apostolaki, E.T.; Fraschetti, S.; Gristina, M.; Knittweis, L.; Martin, C.S.; Pergent, G. Seagrass meadows (Posidonia oceanica) distribution and trajectories of change. Sci. Rep.; 2015; 5, 12505. [DOI: https://dx.doi.org/10.1038/srep12505]
11. Arroyo, M.C.; Barrajón, A.; Brun, F.; del Castillo, F.; de la Rosa, J.; Díaz-Almela, E.; Fernández-Casado, M.; Hernández, I.; Moreno, D.; Pérez-Lloréns, J.L. et al. Praderas de angiospermas marinas de Andalucía. Atlas de las Praderas Marinas de España; Ruiz, J.M.; Guillén, J.E.; Ramos-Segura, A.; Otero, M.M. IEO/IEL/UICN: Murcia-Alicante-Málaga, Spain, 2015; pp. 313-398.
12. Celebi, B. A study on Posidonia oceanica (L.) Delile, 1813, seagrass meadows in the Levantine Sea. Master’s Thesis; METU Institute of Marine Science: Mersin, Turkey, 2007; 124.
13. Pérez-Lloréns, J.L.; Vergara, J.J.; Olivé, I.; Mercado, J.M.; Conde-Álvarez, R.; Pérez-Ruzafa, Á.; Figueroa, F.L. Autochthonous seagrasses. The Mediterranean Sea: Its history and Present Challenges; Goffredo, S.; Dubinsky, Z. Springer: Dordrecht, The Netherlands, 2014; pp. 137-158. [DOI: https://dx.doi.org/10.1007/978-94-007-6704-1_9]
14. Urra, J.; Gofas, S.; Rueda, J.L.; Marina, P. Molluscan assemblages in littoral soft bottoms of the Alboran Sea (Western Mediterranean Sea). Mar. Biol. Res.; 2010; 7, pp. 27-42. [DOI: https://dx.doi.org/10.1080/17451001003660301]
15. Mateo-Ramírez, Á.; Urra, J.; Marina, P.; Rueda, J.L.; García Raso, J.E. Crustacean decapod assemblages associated with fragmented Posidonia oceanica meadows in the Alboran Sea (Western Mediterranean Sea): Composition, temporal dynamics and influence of meadow structure. Mar. Ecol.; 2016; 37, pp. 344-358. [DOI: https://dx.doi.org/10.1111/maec.12284]
16. Vargas-Yáñez, M.; García-Martínez, M.; Moya, F.; Balbín, R.; López-Jurado, J.; Serra, M.; Zunino, P.; Pascual, J.; Salat, J. Updating temperature and salinity mean values and trends in the Western Mediterranean: The RADMED project. Prog. Oceanogr.; 2017; 157, pp. 27-46. [DOI: https://dx.doi.org/10.1016/j.pocean.2017.09.004]
17. Kinder, T.H.; Parrilla, G. Yes, some of the Mediterranean outflow does come from great depth. J. Geophys. Res. Earth Surf.; 1987; 92, 2901. [DOI: https://dx.doi.org/10.1029/JC092iC03p02901]
18. Lafuente, J.G.; Vargas, J.M.; Sarhan, T.; Candela, J.; Bascheck, B.; Plaza, F. Tide at the eastern section of the Strait of Gibraltar. J. Geophys. Res. Earth Surf.; 2000; 105, pp. 14197-14213. [DOI: https://dx.doi.org/10.1029/2000JC900007]
19. Rodríguez, J. Las reservas marinas en el marco ecológico y oceanográfico del Mediterráneo Occidental. La gestión de los Espacios Marinos en el Mediterráneo Occidental; Guirado, J. Instituto de Estudios Almerienses: Almería, Spain, 1995; pp. 13-28. Available online: https://isni.org/isni/000000005921782X (accessed on 27 December 2022).
20. Sarhan, T. Upwelling mechanisms in the northwestern Alboran Sea. J. Mar. Syst.; 2000; 23, pp. 317-331. [DOI: https://dx.doi.org/10.1016/S0924-7963(99)00068-8]
21. Templado, J. La diversidad marina en España. Biodiversidad: Aproximación a la Diversidad Botánica y Zoológica en España, Segunda época, Tomo IX; Viejo, J.L. Memorias de la Real Sociedad Española de Historia Natural: Madrid, Spain, 2011; pp. 343-362.
22. Rueda, J.L.; Gofas, S.; Aguilar, R.; de la Torriente, A.; Raso, J.E.G.; Iacono, C.L.; Luque, A.; Marina, P.; Mateo-Ramírez,; Moya-Urbano, E. et al. Benthic Fauna of Littoral and Deep-Sea Habitats of the Alboran Sea: A Hotspot of Biodiversity. Alboran Sea—Ecosystems and Marine Resources; Báez, J.C.; Vázquez, J.T.; Camiñas, J.A.; Malouli-Idrissi, M. Springer Nature: Geneva, Switzerland, 2021; pp. 285-358. [DOI: https://dx.doi.org/10.1007/978-3-030-65516-7_9] ISBN 9783030655167
23. Agostini, V.N.; Bakun, A. Ocean triads’ in the Mediterranean Sea: Physical mechanisms potentially structuring reproductive habitat suitability (with example application to European anchovy, Engraulis encrasicolus). Fish. Oceanogr.; 2002; 11, pp. 129-142. [DOI: https://dx.doi.org/10.1046/j.1365-2419.2002.00201.x]
24. García, A. Post-larval research of small pelagic species in the Alboran Sea. A call for safeguarding fry concentration sites. Proceedings of the CGPM—SAC Subcommittee of Marine Environment and Ecosystems (GFCM:SAC13/2011/Inf.5.); Malta, 26 November 29–December 2010.
25. Quintanilla, J.M.; Laiz-Carrión, R.; García, A.; Quintanilla, L.F.; Cortés, D.; Gómez-Jakobsen, F.; Yebra, L.; Salles, S.; Putzeys, S.; León, P. et al. Early life trophodynamic influence on daily growth patterns of the Alboran Sea sardine (Sardina pilchardus) from two distinct nursery habitats (bays of málaga and almería) in the western Mediterranean Sea. Mar. Environ. Res.; 2020; 162, 105195. [DOI: https://dx.doi.org/10.1016/j.marenvres.2020.105195]
26. Baro, J.; García, T.; Serna-Quintero, J.M. Description of artisanal fisheries in northern Alboran Sea. Alboran Sea—Ecosystems and Marine Resources; Báez, J.C.; Vázquez, J.T.; Camiñas, J.A.; Malouli-Idrissi, M. Springer Nature: Geneva, Switzerland, 2021; pp. 521-542. [DOI: https://dx.doi.org/10.1007/978-3-030-65516-7_14] ISBN 9783030655167
27. Rueda, J.L.; Marina, P.; Urra, J.; Salas, C. Changes in the composition and structure of a molluscan assemblage due to eelgrass loss in southern Spain (Alboran Sea). J. Mar. Biol. Assoc. United Kingd.; 2009; 89, pp. 1319-1330. [DOI: https://dx.doi.org/10.1017/S0025315409000289]
28. Rueda, J.L.; Marina, P.; Urra, J.; Salas, C. Effects of eelgrass loss on coastal biodiversity. JMBA Glob. Mar. Environ.; 2009; 10, pp. 12-13.
29. Bay, D. A field study of the growth dynamics and productivity of Posidonia oceanica (L.) delile in Calvi Bay, Corsica. Aquat. Bot.; 1984; 20, pp. 43-64. [DOI: https://dx.doi.org/10.1016/0304-3770(84)90026-3]
30. Pergent, G.; Pergent Martini, C. Phénologie de Posidonia oceanica (L.) Delile dans le bassin Méditerranéen. Ann. De L’institut Océanograph.; 1988; 64, pp. 79-100.
31. Buia, M.C.; Zupo, V.; Mazzella, L. Primary Production and Growth Dynamics in Posidonia oceanica. Mar. Ecol.; 1992; 13, pp. 2-16. [DOI: https://dx.doi.org/10.1111/j.1439-0485.1992.tb00336.x]
32. Sánchez Lizaso, J.L. Estudio de la Pradera de Posidonia oceanica (L.) Delile de la Reserva Marina de Tabarca (Alicante): Fenología y Producción Primaria. Ph.D. Thesis; University of Alicante: Alicante, Spain, 1993; 121.
33. Guidetti, P.; Lorenti, M.; Buia, M.C.; Mazzella, L. Temporal Dynamics and Biomass Partitioning in Three Adriatic Seagrass Species: Posidonia oceanica, Cymodocea nodosa, Zostera marina. Mar. Ecol.; 2002; 23, pp. 51-67. [DOI: https://dx.doi.org/10.1046/j.1439-0485.2002.02722.x]
34. Marbà, N.; Duarte, C.M. Interannual changes in seagrass (Posidonia oceanica) growth and environmental change in the Spanish Mediterranean littoral zone. Limnol. Oceanogr.; 1997; 42, pp. 800-810. [DOI: https://dx.doi.org/10.4319/lo.1997.42.5.0800]
35. Gobert, S. Variations Spatiale et Temporelle de L’herbier à Posidonia oceanica (L.) Delile. (Baie de La Revellata, Calvi, Corse). Ph.D. Thesis; University of Liège: Liege, Belgium, 2002; 207.
36. Peirano, A.; Cocito, S.; Banfi, V.; Cupido, R.; Damasso, V.; Farina, G.; Lombardi, C.; Mauro, R.; Morri, C.; Roncarolo, I. et al. Phenology of the Mediterranean seagrass Posidonia oceanica (L.) Delile: Medium and long-term cycles and climate inferences. Aquat. Bot.; 2011; 94, pp. 77-92. [DOI: https://dx.doi.org/10.1016/j.aquabot.2010.11.007]
37. Lesica, P.; Allendorf, F.W. When Are Peripheral Populations Valuable for Conservation?. Conserv. Biol.; 1995; 9, pp. 753-760. [DOI: https://dx.doi.org/10.1046/j.1523-1739.1995.09040753.x]
38. Gibson, S.Y.; Van Der Marel, R.C.; Starzomski, B.M. Climate Change and Conservation of Leading-Edge Peripheral Populations. Conserv. Biol.; 2009; 23, pp. 1369-1373. [DOI: https://dx.doi.org/10.1111/j.1523-1739.2009.01375.x]
39. Eckert, C.G.; Samis, K.E.; Lougheed, S.C. Genetic variation across species’ geographical ranges: The central–marginal hypothesis and beyond. Mol. Ecol.; 2008; 17, pp. 1170-1188. [DOI: https://dx.doi.org/10.1111/j.1365-294X.2007.03659.x]
40. Bennett, S.; Alcoverro, T.; Kletou, D.; Antoniou, C.; Boada, J.; Buñuel, X.; Cucala, L.; Jorda, G.; Kleitou, P.; Roca, G. et al. Resilience of seagrass populations to thermal stress does not reflect regional differences in ocean climate. New Phytol.; 2021; 233, pp. 1657-1666. [DOI: https://dx.doi.org/10.1111/nph.17885] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34843111]
41. Mota, C.F.; Engelen, A.H.; Serrao, E.A.; Coelho, M.A.G.; Marbà, N.; Krause-Jensen, D.; Gareth, A.P. Differentiation in fitness-related traits in response to elevated temperatures between leading and trailing edge populations of marine macrophytes. PLoS ONE; 2018; 13, e0203666. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18302683][DOI: https://dx.doi.org/10.1371/journal.pone.0203666]
42. King, N.G.; McKeown, N.J.; Smale, D.A.; Wilcockson, D.C.; Hoelters, L.; Groves, E.A.; Stamp, T.; Moore, P.J. Evidence for different thermal ecotypes in range centre and trailing edge kelp populations. J. Exp. Mar. Biol. Ecol.; 2019; 514–515, pp. 10-17. [DOI: https://dx.doi.org/10.1016/j.jembe.2019.03.004]
43. Mateo-Ramírez,; Marina, P.; Moreno, D.; Valero, A.F.A.; Aguilar, R.; Báez, J.C.; Bárcenas, P.; Baro, J.; Caballero-Herrera, J.A.; Camiñas, J.A. et al. Marine Protected Areas and Key Biodiversity Areas of the Alboran Sea and Adjacent Areas. Alboran Sea, Ecosystems and Marine Resources; Báez, J.C.; Vázquez, J.T.; Camiñas, J.A.; Malouli-Idrissi, M. Springer Nature: Geneva, Switzerland, 2021; pp. 819-923. [DOI: https://dx.doi.org/10.1007/978-3-030-65516-7_25] ISBN 9783030655167
44. Ruiz, J.M.; Guillén, E.; Ramos-Segura, A.; Otero, M. Atlas de las Praderas Marinas de España; IEO/IEL/UICN: Murcia-Alicante-Málaga, Spain, 2015; 681.
45. Ramírez, T.; Cortés, D.; Mercado, J.; Vargas-Yañez, M.; Sebastián, M.; Liger, E. Seasonal dynamics of inorganic nutrients and phytoplankton biomass in the NW Alboran Sea. Estuarine, Coast. Shelf Sci.; 2005; 65, pp. 654-670. [DOI: https://dx.doi.org/10.1016/j.ecss.2005.07.012]
46. García-Gómez, C.; Yebra, L.; Cortés, D.; Sánchez, A.; Alonso, A.; Valcárcel-Pérez, N.; Gómez-Jakobsen, F.; Herrera, I.; Johnstone, C.; Mercado, J.M. Shifts in the protist community associated with an anticyclonic gyre in the Alboran Sea (Mediterranean Sea). FEMS Microbiol. Ecol.; 2020; 96, fiaa197. [DOI: https://dx.doi.org/10.1093/femsec/fiaa197]
47. Junta de Andalucía. Red de Información Ambiental de Andalicía. REDIAM, pagina de Descargas. Available online: https://portalrediam.cica.es/descargas (accessed on 4 June 2022).
48. Ministerio de Transporte, Movilidad y Agencia Urbana, Gobierno de España. Datos Oceanográficos. Available online: https://www.puertos.es/es-es/oceanografia/Paginas/portus.aspx (accessed on 15 February 2022).
49. Amoutzpoulou-Schina, H.; Haritonidis, S. Distribution and phenology of the marine phanerogam Posidonia oceanica in the Pagassitikos Gulf, Greece. J. Biol. Res.; 2005; 4, pp. 203-211.
50. Kruzic, P. Variations in Posidonia oceanica meadow structure along the coast of the Dugi Otok Island (eastern Adriatic Sea). J. Mar. Biol. Assoc.; 2008; 88, pp. 883-892. [DOI: https://dx.doi.org/10.1017/S0025315408001410]
51. Pergent-Martini, C.; Leoni, V.; Pasqualini, V.; Ardizzone, G.D.; Balestri, E.; Bedini, R.; Belluscio, A.; Belsher, T.; Borg, J.; Boudouresque, C.F. et al. Descriptors of Posidonia oceanica meadows: Use and application. Ecol. Indic.; 2005; 5, pp. 213-230. [DOI: https://dx.doi.org/10.1016/j.ecolind.2005.02.004]
52. Conde, F.; Seoane, J.A. Corología de las especie de algas en relación a ciertos factores ecológicos en el litoral malagueño. Collect. Bot.; 1982; 13, pp. 783-802.
53. Flores-Moya, A.; Vera-Gonzalez, J.J.; Conde, F. Contribución a la corología de las macroalgas marinas bentónicas del litoral malagueño. Acta Bot. Malac.; 1989; 14, pp. 199-201. [DOI: https://dx.doi.org/10.24310/abm.v14i.9375]
54. Jeffrey, S.; Humphrey, G.T. New spectrophotometric equations for determining chlorophylls a, b, c1 and c2 in higher plants, algae and phytoplankton. Biochem. Physiol. Pflanzen; 1975; 167, pp. 191-194. [DOI: https://dx.doi.org/10.1016/S0015-3796(17)30778-3]
55. Lopez y Royo, C.; Casazza, G.; Pergent-Martini, C.; Pergent, G. A biotic index using the seagrass Posidonia oceanica (BiPo), to evaluate ecological status of coastal waters. Ecol. Indic.; 2010; 1, pp. 380-389. [DOI: https://dx.doi.org/10.1016/j.ecolind.2009.07.005]
56. Anderson, M.J. A new method for non-parametric multivariate analysis of variance. Austral Ecol.; 2001; 26, pp. 32-46. [DOI: https://dx.doi.org/10.1111/j.1442-9993.2001.01070.pp.x]
57. Moreno, D.; Guirado, J. Nuevos datos sobre la distribución de las fanerógamas marinas en las provincias de Almería y Granada (SE España). Acta Bot. Malacit.; 2003; 28, pp. 105-120. [DOI: https://dx.doi.org/10.24310/abm.v28i0.7270]
58. Báez, J.C.; Vázquez, J.T.; Camiñas, J.A.; Malouli, M. Alboran Sea—Ecosystems and Marine Resources; Springer Nature: Geneva, Switzerland, 2021; 939. [DOI: https://dx.doi.org/10.1007/978-3-030-65516-7] ISBN 9783030655167
59. García Raso, J.E.; Gofas, S.; Salas Casanova, C.; Manjón-Cabeza, E.; Urra, J.; García Muñóz, J.E. El Mar Más rico de Europa: Biodiversidad del Litoral Occidental de Málaga Entre Calaburras y Calahonda; Consejería de Medio Ambiente, Junta de Andalucía: Sevilla, Spain, 2010; 138. ISBN 9788492807499
60. Urra, J.; Rueda, J.L.; Gofas, S.; Marina, P.; Salas, C. A species-rich molluscan assemblage in a coralligenous bottom of the Alboran Sea (south-western Mediterranean): Intra-annual changes and ecological considerations. J. Mar. Biol. Assoc. United Kingd.; 2011; 92, pp. 665-677. [DOI: https://dx.doi.org/10.1017/S0025315411001391]
61. Lazzari, P.; Solidoro, C.; Ibello, V.; Salon, S.; Teruzzi, A.; Béranger, K.; Colella, S.; Crise, A. Seasonal and inter-annual variability of plankton chlorophyll and primary production in the Mediterranean Sea: A modelling approach. Biogeosciences; 2012; 9, pp. 217-233. [DOI: https://dx.doi.org/10.5194/bg-9-217-2012]
62. Ramírez, T.; Muñoz, M.; Reul, A.; García-Martínez, M.C.; Moya, F.; Vargas-Yáñez, M.; Bautista, B. The Biogeochemical Context of Marine Planktonic Ecosystems. Alboran Sea, Ecosystems and Marine Resources; Báez, J.C.; Vázquez, J.T.; Camiñas, J.A.; Malouli-Idrissi, M. Springer Nature: Geneva, Switzerland, 2021; pp. 207-246. [DOI: https://dx.doi.org/10.1007/978-3-030-65516-7_7] ISBN 9783030655167
63. Sand-Jensen, K. Biomass, net production and growth dynamics in an eelgrass (Zostera marina L.) population in Vellerup Vig, Denmark. Ophelia; 1975; 14, pp. 185-201. [DOI: https://dx.doi.org/10.1080/00785236.1975.10422501]
64. Fonseca, M.S.; Zieman, J.C.; Thayer, G.W.; Fisher, J.S. The role of current velocity in structuring eelgrass (Zostera marina) meadows. Estuar. Coast. Shelf Sci.; 1983; 17, pp. 367-380. [DOI: https://dx.doi.org/10.1016/0272-7714(83)90123-3]
65. Reyes, J.; Sansón, M.; Carrillo, A. Leaf phenology, growth and production of the seagrass Cymodocea nodosa at El Médano (South of Tenerife, Canary Islands). Bot. Mar.; 1995; 38, pp. 457-465. [DOI: https://dx.doi.org/10.1515/botm.1995.38.1-6.457]
66. Turk, R.; Vukovic, A. Phenology of Posidonia oceanica (L.) Delile in the Gulf of Koper (Gulf of Trieste), North Adriatic. Rapp. Comm. Int. Pour L’exploration Sci. De La Mer Méditerranée; 1998; 35, pp. 592-593.
67. Ballesta, L.; Pergent, G.; Pasqualini, V.; Pergent-Martini, C. Distribution and Dynamics of Posidonia oceanica Beds along the Albères Coastline; Elsevier: Amsterdam, The Netherlands, 2000; Volume 323, pp. 407-414. [DOI: https://dx.doi.org/10.1016/s0764-4469(00)00135-9]
68. Rubegni, F.; Franchi, E.; Lenzi, M. Relationship between wind and seagrass meadows in a non-tidal eutrophic lagoon studied by a Wave Exposure Model (WEMo). Mar. Pollut. Bull.; 2013; 70, pp. 54-63. [DOI: https://dx.doi.org/10.1016/j.marpolbul.2013.02.012] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23522697]
69. Pace, M.; Borg, J.A.; Galdies, C.; Malhotra, A. Influence of wave climate on architecture and landscape characteristics of Posidonia oceanica meadows. Mar. Ecol.; 2016; 38, e12387. [DOI: https://dx.doi.org/10.1111/maec.12387]
70. Mabrouk, L.; Hamza, A.; Sahraoui, H.; Bradai, M.N. Donnees sur les caracteristiques et la phenology de l’herbier de Posidonia oceanica (L.) Delile sur les côtes de mahdia (región est de la Tunisie). Bull. Inst. Natl. Sci. Tech. Oceanogr. Peche Salammbo; 2009; 36, pp. 139-147.
71. Abbate, M.; Peirano, A.; Ugolini, U. Structural changes in Posidonia oceanica leaves along the coast of Liguria (Italy): Response to environmental stress?. Biol. Mar. Medit.; 2000; 7, pp. 320-323.
72. Giovannetti, E.; Lasagna, R.; Montefalcone, M.; Bianchi, C.; Albertelli, G.; Morri, C. Inconsistent responses to substratum nature in Posidonia oceanica meadows: An integration through complexity levels?. Chem. Ecol.; 2008; 24, pp. 83-91. [DOI: https://dx.doi.org/10.1080/02757540801966439]
73. Sghaier, Y.R.; Zakhama, R.S.; Charfi, F.C. Status of Posidonia oceanica meadows along the Eastern coast of Tunisia. Biol. Mar. Medit.; 2006; 13, pp. 85-91.
74. Pergent, G.; Gucu, A.C.; Sakinan, S.; Pergent-Martini, C. Cartographie et surveillance de l´herbier de Mersin (Turquie). Programme MedPosidonia /SPA/RAC, Fondation d´entrepise. TOTAL pour La Biodiversité Et La Mer, Memorandun dáccord 01/SPA/CAR, MedPosidonia. 2008; pp. 1-16.
75. Dural, B. Phenological observations on Posidonia oceanica (L.) Delile medows along the coast of Akkum (Sıgacık Bay, Aegean Sea, Turkey). J. Black Sea/Medit. Environ.; 2010; 16, pp. 133-144.
76. Peduzzi, P.; Vukovic, A. Primary production of Cymodocea nodosa in the Gulf of Trieste (Northern Adriatic Sea): A comparison of methods. Mar. Ecol. Prog. Ser.; 1990; 64, pp. 197-207. [DOI: https://dx.doi.org/10.3354/meps064197]
77. Semroud, R.; Benkorteby, N.; Tamouza, F.Z. Phenologie de Posidonia oceanica dans la region d´Alger (Algerie): Donnes preliminaires sur la biometriefoliare. Rapp. Comm. Int. Pour L’exploration Sci. De La Méditerranée; 1990; 32, B-I16: 10-11.
78. Pergent, G.; Djellouli, A.; Hamza, A.A.; Ettayeb, K.S.; El Mansouri, A.A.; Talha, F.M.; Hamza, M.A.A.; Pergent-Martini, C.; Platini, F. Characterization of the benthic vegetation in the FarwàLaggon (Libya). J. Coast. Conser.; 2002; 8, pp. 119-126. [DOI: https://dx.doi.org/10.1652/1400-0350(2002)008[0119:COTBVI]2.0.CO;2]
79. Borg, J.A.; Attrill, M.J.; Rowden, A.A.; Schembri, P.J.; Jones, M.B. Architectural characteristics of two bed types of the seagrass Posidonia oceanica over different spatial scales. Estuar. Coast. Shelf Sci.; 2005; 62, pp. 667-678. [DOI: https://dx.doi.org/10.1016/j.ecss.2004.10.003]
80. Urra, J.; Mateo,; Marina, P.; Rueda, J.L.; Raso, J.E.G. First records of Posidonia oceanica flowering at its westernmost distributional limit (Málaga, Alboran Sea). Bot. Mar.; 2011; 54, pp. 101-104. [DOI: https://dx.doi.org/10.1515/bot.2011.010]
81. Balestri, E. Flowering of the seagrass Posidonia oceanica in a north-western Mediterranean coastal area: Temporal and spatial variations. Mar. Biol.; 2004; 145, pp. 61-68. [DOI: https://dx.doi.org/10.1007/s00227-004-1301-2]
82. Montefalcone, M.; Giovannetti, E.; Morri, C.; Peirano, A.; Bianchi, C.N. Flowering of the seagrass Posidonia oceanica in NW Mediterranean: Is there a link with solar activity?. Mediterr. Mar. Sci.; 2013; 14, 416. [DOI: https://dx.doi.org/10.12681/mms.529]
83. Balestri, E.; Vallerini, F.; Lardicci, C. Recruitment and patch establishment by seed in the seagrass Posidonia oceanica: Importance and conservation implications. Front Plant Sci.; 2017; 8, pp. 1-12. [DOI: https://dx.doi.org/10.3389/fpls.2017.01067] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28670323]
84. Díaz-Almela, E.; Marbà, N.; Duarte, C.M. Consequences of Mediterranean warming events in seagrass (Posidonia oceanica) flowering records. Glob. Chang. Biol.; 2006; 13, pp. 224-235. [DOI: https://dx.doi.org/10.1111/j.1365-2486.2006.01260.x]
85. Ruiz, J.; Marín-Guirao, L.; García-Muñoz, R.; Ramos-Segura, A.; Bernardeau-Esteller, J.; Pérez, M.; Sanmartí, N.; Ontoria, Y.; Romero, J.; Arthur, R. et al. Experimental evidence of warming-induced flowering in the Mediterranean seagrass Posidonia oceanica. Mar. Pollut. Bull.; 2017; 134, pp. 49-54. [DOI: https://dx.doi.org/10.1016/j.marpolbul.2017.10.037] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29102072]
86. Balestri, E.; Cinelli, F. Sexual reproductive success in Posidonia oceanica. Aquat. Bot.; 2003; 75, pp. 21-32. [DOI: https://dx.doi.org/10.1016/S0304-3770(02)00151-1]
87. Balestri, E.; Vallerini, F.; Lardicci, C. On the unusual flowering of plagiotropic shoots in the seagrass Posidonia oceanica. Aquat. Bot.; 2005; 82, pp. 82-88. [DOI: https://dx.doi.org/10.1016/j.aquabot.2005.03.001]
88. Cabaço, S.; Santos, R. Seagrass reproductive effort as an ecological indicator of disturbance. Ecol. Indic.; 2012; 23, pp. 116-122. [DOI: https://dx.doi.org/10.1016/j.ecolind.2012.03.022]
89. Consejería de Agricultura, Ganadería, Pesca y Desarrollo Sostenible, Junta de Andalucía. Informe Final de Resultados del Programa de Gestión Sostenible del Medio Marino Andaluz, 2020. 168.Available online: https://www.juntadeandalucia.es/medioambiente/portal/landing-page-documento/-/asset_publisher/jXKpcWryrKar/content/informes-regionales-sobre-gesti-c3-b3n-sostenible-del-medio-marino-andaluz-2008-2018-/20151 (accessed on 14 April 2022).
90. Urra, J.; Marina, P.; Rueda, J.L.S.O.S. de las praderas de fanerógamas marinas en la Costa del Sol. Quercus; 2008; 270, pp. 28-37.
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Details
; Pablo, Marina 1 ; Martín-Arjona, Alejandro 1 ; Bañares-España, Elena 2 ; García Raso, José E 3
; Rueda, José L 1 ; Urra, Javier 1 1 Centro Oceanográfico de Málaga, Instituto Español de Oceanografía-Consejo Superior de Investigaciones Científicas, Puerto Pesquero s/n, 29640 Málaga, Spain
2 Departamento de Botánica y Fisiología Vegetal, Campus de Teatinos s/n, Universidad de Málaga, 29071 Málaga, Spain
3 Departamento de Biología Animal, Campus de Teatinos s/n, Universidad de Málaga, 29071 Málaga, Spain