Next Article in Journal
Ecological Quality Status Evaluation of Port Sea Areas Based on EW-GRA-TOPSIS Model
Next Article in Special Issue
Pollen and Flora as Bioindicators in Assessing the Status of Polluted Sites: The Case Study of the Mantua Lakes (SIN “Laghi di Mantova e Polo Chimico”; N Italy)
Previous Article in Journal
Revolutionizing IC Genset Operations with IIoT and AI: A Study on Fuel Savings and Predictive Maintenance
Previous Article in Special Issue
Status and Trend of the Main Allergenic Pollen Grains and Alternaria Spores in the City of Rome (2003–2019)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Sustainability
Volume 15
Issue 11
10.3390/su15118807
sustainability-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Human-Impact Gradients through Anthropogenic Pollen Indicators in a Mediterranean Mosaic Landscape (Balearic Islands)

by
Gabriel Servera-Vives
1,2,*,
Maurici Mus Amezquita
3,
Grant Snitker
4,*,
Assunta Florenzano
2,5,
Paola Torri
2,
Maurici Ruiz
6,7 and
Anna Maria Mercuri
2,5
1
ArqueoUIB, Department of Historical Sciences and Theory of Art, University of the Balearic Islands, Carretera de Valldemossa Km 7.5, 07122 Palma, Mallorca, Spain
2
Laboratory of Palynology and Paleobotany, Department of Life Sciences, Università degli Studi di Modena e Reggio Emilia, Via Campi 287, 41125 Modena, Italy
3
Department of Biology, University of the Balearic Islands, Carretera de Valldemossa Km 7.5, 07122 Palma, Mallorca, Spain
4
Center for Applied Fire and Ecosystem Science, New Mexico Consortium 4200 W, Jemez Rd., Suite 301, Los Alamos, NM 87544, USA
5
NBFC, National Biodiversity Future Center, 90133 Palermo, Italy
6
Mediterranean Hydrological and Ecogeomorphological Connectivity Research Team, Department of Geography, University of the Balearic Islands, Carretera de Valldemossa Km 7.5, 07122 Palma, Mallorca, Spain
7
Institute of Agro-Environmental and Water Economy Research—INAGEA, University of the Balearic Islands, Carretera de Valldemossa Km 7.5, 07122 Palma, Mallorca, Spain
*
Authors to whom correspondence should be addressed.
Sustainability 2023, 15(11), 8807; https://doi.org/10.3390/su15118807
Submission received: 19 April 2023 / Revised: 18 May 2023 / Accepted: 26 May 2023 / Published: 30 May 2023
(This article belongs to the Special Issue Palynology for Sustainability: A Classical and Versatile Tool for New Challenges)
Browse Figures
Versions Notes

Abstract

:
This paper proposes new anthropogenic pollen indicators for the Balearic Islands and attempts to assess gradients of human impact on vegetation in Mediterranean islands. A combination of modern pollen analogue studies, complemented by phytosociological descriptions and ordination techniques using quantitative and presence/absence data was used. Redundancy analysis allowed us to evaluate the relationships between pollen types and significant environmental variables and propose regional (e.g., Centaurea, Rubus, Plantago lanceolata-t) and local/microregional anthropogenic pollen indicators (e.g., Cerealia, Poygonum aviculare, Matricaria-t). Additionally, an anthropogenic index score (AIS) for each sampled location was calculated to correlate each pollen type to a specific degree of human impact: (a) low (e.g., Cerastium-t, Erica arborea-t, Cistus albidus), (b) moderate (e.g, Sinapis-t, Sanguisorba minor-t, Plantago bellardii-t), (c) high (e.g., Papaveraceae undiff., Dipsacaceae, Secale-t). This paper contributes to a further understanding of land-use dynamics and to defining the degree of impact, which is especially necessary to assess colonization and anthropization rhythms in Mediterranean island environments.

Graphical Abstract

1. Introduction

The long-term history of human impacts on vegetation and the development of cultural landscapes have been a major research topic in Europe since the second half of the 20th century [1]. In this context, anthropogenic pollen indicators (API) and the indicator-species approach have been valuable tools for understanding human-induced vegetation changes in the past. Karl-Ernst Behre’s [2] seminal approach proposed a list of taxa indicators of anthropogenic vegetation types and land use, providing a way to identify human impact through time via pollen assemblages. While regional and local-specific anthropogenic pollen indicators have proven to be essential in evaluating human impacts, classical indicators, such as those proposed by Behre, are still of great value and common in the palynological praxis for characterizing main trends in anthropogenic dynamics recorded on pollen diagrams, especially for detecting ancient human activities occurring in a forest-dominated environment [3]. Nevertheless, several authors suggest the necessity of establishing region-specific anthropogenic pollen indicators to further decipher past human practices through pollen analysis [4]. In this sense, significant research on modern pollen analogues has been developed in many areas of the world to further understand source areas (regional vs. local indicators) and human activities, namely pollen types that are indicators of specific practices such as agriculture or grazing [4,5,6]. Instead, Mercuri et al. [7] observed which synanthropic plants are always or recurrently represented in the historical and prehistoric contexts they studied, being evidence of human presence and activities in a certain area. While this bioindicator-based approach is common in many Central and Northern European areas, it has scarcely been applied in the Mediterranean region. Recent research also proposes new indices such as LUP (agricultural land use probability) as a general index to assess and quantify human impact intensity on European ecosystems [8]. Detailed reviews on different human impact indices and ratios have been published in Deza-Araujo et al. [8,9].
Overall, Mediterranean landscapes are characterized by a mosaic-like landscape structure, characterized by alternating forests, shrubby vegetation, and open-lands where a wide range of human activities have traditionally been developed (e.g., forestry, agriculture, grazing, mining, etc.). This type of landscape physiognomy in the Mediterranean region results in highly valuable environments due to its great biodiversity [10,11,12]. The origin of mosaic landscapes is linked to the onset of sylvoagropastoral activities during the Neolithic period, and especially since the mid-Holocene [13,14]. These legacy landscapes result from long-term sustainable practices through time and require rural landholders to preserve them [15].
Mediterranean island environments are particularly interesting for the study of Holocene plant–cultural interactions due to their particular biogeographical and sociocultural histories and because islands are more vulnerable to human and climate-induced changes [16]. Islands may be considered fragile environments where human impact and other environmental triggers may cause a tipping point, breaking their resilience status due to specific evolutionary adaptations that make to island species particularly vulnerable to human pressure [17]. For instance, the Balearic Archipelago seems to have been settled by humans relatively late when compared to other Mediterranean islands, with the scholarly accepted date for human arrival at ca. 2500 cal BC. Moreover, early agropastoral activities occurring in a context of forested environments and small-scale shifting human activities are difficult to detect from solely a palynological point of view, and therefore require multiproxy and interdisciplinary research to evaluate robust evidence of land use. In this sense, it has been suggested that the transition “from influence to impact”, during the Late Holocene and the development of Chalcolithic and Bronze Age cultures (III-II millennium cal BC for the Balearics), has been detected through pollen analysis in different Mediterranean contexts [14]. In general, to assess human impact on past environments, especially to further assess cultural change and island colonization processes in Mediterranean island environments, it is necessary to obtain reliable bioindicators.
In this paper, we present research on modern pollen analogues conducted in the Balearic Islands to determine pollen types related to local/microregional and regional anthropogenic activities, which in turn provide useful indicators of human land use for Mediterranean and circum-Mediterranean sites. In previous research, we analyzed pollen–vegetation relationships for specific vegetation types from the Balearic mosaic landscape, proposing their indicator value and degree of representativity [18]. However, the identification of anthropogenic pollen indicators still needs an exhaustive evaluation to assess their role and significance as apophytes in Mediterranean environments. We investigate the potential indicator value of pollen types to interpret pollen-derived human impact gradients in Mediterranean island environments.

2. Study Area

The Balearic Islands are an archipelago located in the Western Mediterranean Sea formed by a group of 151 islands and islets, with only the four largest being inhabited. The vegetation of the Balearic Islands is characterized by a mosaic-like pattern (Figure 1) which results from environmental and anthropogenic dynamics [19]. The climate of the Balearic Archipelago is typically Mediterranean, with seasonality (dry and hot summers and maximal rainfall during the autumn season). The highest rainfall is recorded in the Serra de Tramuntana mountain range (ca. 1500–1400 mm) in Mallorca, while the lowest values (300–450 mm) are recorded in Formentera [20]. These Mediterranean mosaics have been constructed through natural and social processes since the middle Holocene and are best understood using a longue durée perspective [21,22]. Therefore, the Balearic landscape is primarily composed of heterogeneous forests, shrublands, grasslands, tree orchards, agricultural, and pastoral areas. The prevalent tree vegetation are evergreen forests, woodlands, and sclerophyllous shrublands dominated by holm oaks, pins, olive trees, mastic, and junipers [23]. The primary agropastoral uses are based on dry extensive farming including annual cereal crops, with dispersed orchards and post-cereal-harvest livestock grazing [19]. Arboreal cultivars have traditionally been based on fig, olive, carob, and almond tree agriculture. Since the second half of the 20th century, many coastal areas have been widely urbanized in response to the tourism boom, and former crop fields and pastures are used for non-farming and touristic activities [24].
As a result of this long-term, co-dependent relationship between people and ecosystems, Balearic landscapes and vegetation are well adapted to human-caused disturbances, such as the preference or tolerance of sclerophyllous dry-adapted species for grazing environments [25,26] (Figure 2). The spread of garrigues and maquis during the Late Holocene seems to be correlated to both drier climate conditions and an increase in agropastoral activities [22]. Vegetation community distribution may be explained through heterogeneous environmental conditions and different degrees of impact. Some synanthropic communities include Urtico membranaceae-Smyrnietum olusatri A. & O. Bolòs in O. Bolòs & Molinier 1958, Resedo albae-Chrysanthemetum coronarii O. Bolòs & Molinier 1958 or Reichardio gracilis-Stipetum capensis Costa & Loidi 1992 corr. [23], while key synanthropic taxa from the Balearic landscape include Centaurea cyanus L., Asphodelus fistulosus L., Diplotaxis erucoides (L.) DC, Galactites tormentosa Moench or Plantago lanceolata L. Vegetation community distribution may be explained through heterogeneous environmental conditions and different degrees of impact. Synanthropic vegetation is broadly present in the archipelago, including communities of Artemisietea vulgaris Lohmeyer, Preising & Tüxen ex von Rochow 1951 and Stellarietea medie Tüxen, Lohmeyer & Preising ex von Rochow 1951 classes, xerophytic grasslands of the Lygeo-Stipetea Rivas-Martínez 1978 class or meadow grasslands communities of Molinio-Arrhenatheretea Tüxen 1937 [23]. In addition to some export of synanthropic taxa to other neighboring natural communities with the characteristic mosaic vegetation arrangement of the Balearic Islands, there are species that intrude deeply as companions of plant formations favored by human disturbances; we highlight taxa such Asphodelus aestivus Brot., Arisarum vulgare Targ.-Torz., Carlina corymbosa L., or Plantago lanceolata L. Additionally, many non-native species are listed in Balearic botanical research on natural vegetation communities, resulting from intentional and unintended anthropogenic introduction [26,27,28], which are clearly favored by human activities. Nevertheless, it is still difficult to discern between archaeophytes and neophytes in the current state of botanical-paleoenvironmental research integration.

3. Materials and Methods

Our methods rely on a combination of vegetation surveys, high-morphological resolution pollen identification and numerical analyses to propose anthropogenic indicators to evaluate long-term Holocene human impacts and environmental sustainability in Mediterranean island environments. Further details of fieldwork, laboratory treatment and the selection of environmental variables are reported in Servera-Vives et al. [18].

3.1. Fieldwork

Fifty-six modern pollen samples were collected within the main islands of the Balearic Archipelago: Mallorca, Menorca, Eivissa, Formentera and Cabrera, in 2016 and 2018 [18]. The locations were chosen to document the heterogeneity of major vegetation types of the Balearic mosaic-like landscapes in relation to a wide range of human activities and anthropization intensity. The studied samples consisted of moss cushions, as they may record several years of pollen rain on the sampling location [29,30]. Surface sedimentary samples were collected when local environmental constraints prevented the moss cushions from growing [31]. At least 5 moss cushions or surface sedimentary samples were combined into one sample from each sampling point and phytosociological descriptions were conducted in a 100 m2 square using the Braun-Blanquet [32] system.

3.2. Environmental Variables

Environmental variables were recorded at each sampling site to explore the modern pollen–vegetation–anthropization relationship. Additional variables were also acquired through regional and international geospatial databases of landcover and other landscape metrics. A total of 34 environmental variables were collected for each sampling point, including vegetation type, meteorology, soil, landscape openness, fire activity, topography, landscape metrics and land use information. Detailed information for all environmental variables is described in Table S1 of the Supplementary Material.

3.3. Pollen Analysis

Both moss cushions and surface sedimentary samples were collected to perform pollen analyses on the current pollen rain. Pollen samples were dried overnight. While moss samples were treated following the standard procedures in palynology, including KOH, 200-µm sieving and acetolysis [33,34], surface sedimentary samples were analyzed using a 7-μm nylon handheld sieve and heavy-liquid separation. Pollen counts and identification were performed using a 400×, 600× and 1000× light microscope. Pollen percentages are based on pollen sum, excluding Pinus due to its overrepresentation. The taxonomic identification of pollen grains was based on pollen atlases [35,36], morphological keys [37,38] and the pollen reference collection of the Laboratorio di Palinologia e Paleobotanica (Università degli Studi di Modena e Reggio Emilia). In this work, cereals are recognized as Poaceae pollen grains with diameter larger than 40 µm, a pore annulus ≥ 12 µm [34] and scabrate exine decoration [35,39].

3.4. Statistical Analysis

3.4.1. Data Preparation and Statistical Software

A total of 63 pollen types were selected from all modern pollen sampling locations based on their role as potential anthropogenic indicators, as identified in the specialized literature (see, e.g., [2,4,6,7,8,9,40,41]). A square-root transformation was applied to the selected data to reduce the influence of the most abundant taxa, following Ejarque et al. [6], López-Sáez et al. [42] and Servera-Vives et al. [18]. All data transformations, analyses and visualizations were performed using R, version 4.2.0 [43], the vegan package, version 2.6-4 [44] and the ggplot2 package version 3.4.0 [45].

3.4.2. Redundancy Analysis (RDA) of Pollen Taxa and Environmental Variables

Ordination techniques were used to evaluate the relationship between modern pollen sampling sites, selected pollen taxa and the environmental variables recorded at each sampling site. Detrended correspondence analysis (DCA) was performed to determine the most appropriate ordination techniques for these data, following Smilauer and Lepš [46], Ejarque et al. [6], López-Sáez et al. [47] and others. The resulting gradient length of the selected pollen dataset was 2.58 standard deviation of species turnover units, indicating that a linear ordination method, such as redundancy analysis, is suitable for these data. Additionally, the numerous environmental variables compiled for this analysis (n = 48) increases the likelihood of multicollinearity among each of these predictors, which in turn can obscure relationships between all variables. Variance inflation factors (VIF) were calculated for all environmental variables to assess redundancy in our measures. Environmental variables with VIF values greater than 20.0 were considered highly collinear and removed, leaving eight variables for redundancy analysis. The RDA of pollen, sites, and environmental variables was conducted using a combination forward/backward stepwise model selection, after which ANOVA permutation tests were performed to identify the influential of each predictor on the structure of the pollen dataset. Pearson’s correlation coefficients (r) were calculated to evaluate the relationship between pollen taxa and each environmental variable using the orthogonal linear combinations of the explanatory variables from the RDA.

3.4.3. Constructing the Anthropogenic Intensity Score (AIS) for Sample Sites

We evaluate the degree of anthropization at each sample site by creating an anthropogenic intensity score (AIS), which accounts for the values of the most influential environmental variables at each site, as identified through the RDA. The steps in creating the AIS are outlined below:
(1)
At each sample site, the values of the 8 most influential environmental variables were compiled and transformed to equalize the magnitude of each measure. For example, agropastoral use is a binary variable with values ranging from 0 to 1, while herbivory pressure is an ordinal variable with values ranging from 0 to 2. To give both variables equal magnitude in constructing the AIS, the values of the agropastoral use variable are multiplied by 2 to match the range of values of the herbivory pressure variable. This logic was applied to all variables in our dataset.
(2)
The resulting values of each environmental variable were weighted by the p-values calculated by the ANOVA permutation tests applied to the results of the RDA, following Equation (1). Weighting these variables gave the more influential environmental predictors greater effect in determining the AIS than less influential environmental variables.
W e i g h t e d   V a r i a b l e = T r a n s f o r m e d   V a r i a b l e × 1 − p - value
(3)
Weighted environmental variables were then summed for each sample site to create a raw anthropogenic intensity score value.
(4)
The final AIS was calculated by rescaling the raw anthropogenic intensity values for each site to a 0 to 1 index, with 1 representing the highest anthropogenic intensity score and 0 the lowest anthropogenic intensity score.

3.4.4. Evaluating Anthropogenic Indicator Gradients

Final anthropogenic intensity scores were ranked by dividing the values into three equal quantiles to represent low (0.00–>0.33), moderate (≤0.33–>0.66) and high (≤0.66–1.00) degrees of anthropogenic influence on vegetation communities at each sample site (Figure 1). Pollen taxa were then summed by their abundance (%) in sample sites belonging to each AIS quantile rank category. Principal components analysis (PCA) and a biplot were used to visually illustrate species abundance by AIS rank.
The relationship between each pollen taxon and AIS rank category provides information on its strength as a potential indicator of diverse anthropogenic activities in the fossil pollen record. We constructed an initial categorization of the taxa that are the most correlated to high, moderate and low AIS using K-means clustering (K = 4) of the summed pollen taxa by gradient category. Four clusters were used to group these data to allow for taxa that were weakly correlated with high, moderate and low AIS to be captured in an intermediate category. While K-means provides a general indicator of a pollen taxon’s relationship to anthropogenic activities, it does not specifically evaluate a taxon’s indicator strength. Subsequently, each pollen taxon’s indicator strength was evaluated by mapping its distance in the PCA factorial plane to an “idealized” taxon that is perfectly correlated with either high, moderate or low AIS values. The distances between taxa and the idealized taxon were calculated and scaled between 0 and 1, which illustrates the indicator strength of each taxon for its AIS group. These values were then displayed as a color gradient in a PCA biplot and in a series of silhouette plots for each group. Finally, pollen sums of key ruderal and nitrophilous taxa were plotted in the PCA biplot as supplementary passive variables (e.g., Plantago S U M ) to assess the potential role of some pollen aggregations in lower taxonomic resolution than pollen types as potential API. This is especially important to be able to evaluate human impact through pollen analysis in on-site sedimentary records, where palynomorph’s preservation generally do not allow high taxonomic resolution [7,48].
Additionally, the local/microregional vs. regional indicator value of each pollen type was evaluated by considering the Davis [49] and fidelity/dispersibility [50,51] indices obtained in Servera-Vives et al. 2022. More specifically, in this paper, local vs. regional indicators are proposed either considering pollen types classed as strongly associated taxa (SAT) or associated taxa (AT) in the Davis indices calculation or high fidelity (>50) + low dispersibility (<50) values. Davis indices and fidelity/dispersibility indices were calculated for all taxa based on their presence or absence in both the pollen and observed vegetation datasets [18].

4. Results

4.1. Pollen Type-Environmental Variables Correlation

Pearson’s correlation coefficients (r) and redundancy analysis values indicating moderate to high (>0.5) correlation between pollen types and environmental variables are shown in Table 1. Pasture area (sqm) and landscape openness show similar correlations with PC1 and PC2 and includes comparable pollen types such as Papaveraceae undiff., Dipsacaceae, Vitis and Trifolium-t. Distance to forest shows the highest correlation with Plantago cf. crassifolia, Verbascum-t, Cistus salviifolius and Hornungia-t, while agropastoral use correlates with Avena-Triticum-group, Plantago lanceolata-t, Rubus and Galium-t. Here, we also propose potential local/microregional and regional indicators based on fidelity/dispersibility, and Davis indices obtained in Servera-Vives et al. [18] and discussed in Section 5.1. Potential local indicators of open and anthropized habitats include Poaceae undiff., Polygonum aviculare-t, Avena/Triticum-group, Matricaria-t, and Urtica dioica-t.

4.2. Redundancy Analysis (RDA)

After the calculation of VIF to test for collinearity, a total of eight environmental variables were retained for redundancy analysis. The two first RDA axes explain a total of 18.63% of the variance within the dataset. Forward–backward selection results in three statistically significant environmental variables (p-value ≤ 0.05), namely distance to forest, pasture size (square meters) and landscape openness. While less significant, agropastoral use shows noticeable importance in structuring our dataset, while trampling, fire, distance to urban and herbivory pressure demonstrate low significance (Table 2).
The first RDA axis shows a gradient from closed to open environments (Figure 3). It contrasts most of the sites, from maquis, garrigues, oak forest and Buxus formations demonstrating positive values to coastal saltmarshes, some anthropized habitats, dune habitats and wet pastures in tall humid grassland and freshwater riverine scrubs showing negative values on the axis. The second RDA axis is most related to anthropogenic environments, with garrigues and maquis demonstrating positive values and oak forest, box formations and coastal salt environments with negative values. Pollen types in the negative RDA1 and positive RDA2 values quadrant (upper left) are related to pastoral activity (pasture square meters and agropastoral use) and landscape openness. In this sense, Plantago lanceolata-t, Cichorieae, Avena/Triticum-group, Plantago albicans, Rumex, Plantago sp. and, to a lesser extent, Hordeum-t, Cistus monspeliensis, Asphodelus, Erica arborea-t, Matricaria-t and Sinapis-t are correlated to agropastoral use. Other taxa placed in the negative RDA1 and RDA2 axis values (lower left quadrant), such as Verbascum-t, Plantago cf. crassifolia and Plantago coronopus-t, are related to open areas (distance to forest) and, less significatively, to trampling. Ranunculaceae undiff. Artemisia, Linaria-t and Erica arborea seem to be correlated with fire occurrence (lower right quadrant), and Plantago SUM and Plantago afra to distance to urban areas and herbivory pressure (upper right quadrant).

4.3. K-means Clustering and Unconstrained PCA with Anthropogenic Impact Gradient Categories

K-means clustering was used to create four groups to describe the types and intensity of anthropogenic activities characterizing the sample sites, which we have termed (i) low, (ii) moderate, (iii) high and (iv) intermediate anthropogenic impacts. (i) The first group is comprised of indicators of low-levels of anthropogenic impacts on vegetation and includes a total of 10 pollen types, such as Erica, Cerastium-t and Medicago. (ii) Moderate impact indicators are represented by a group of 15 pollen types, including Sinapis-t, Sanguisorba minor-t and Plantago bellardii-t. (iii) High anthropogenic impact indicators include 12 pollen types, such as Papaveraceae, Dipsacaceae and Secale-t. (iv) The final group of taxa are not clearly associated with a specific degree of human impact, and include a total of 24 pollen types (e.g., Galium-t, Urtica dioica-t., Cistus, etc.). Pollen sums of key botanic genera and families, also included in this analysis, are for instance sums (SUM) of Brassicaceae, Cerealia-t and Plantago, which may be considered good indicators of high human pressure, while Erica SUM is well represented in low impact habitats. The results obtained from K-means clustering are plotted in an unconstrained PCA, where pollen types are drawn as observations and the degrees of impact obtained by the AIS calculation are plotted as variables (Figure 4A). To give more details about the performance of each pollen type within each group, a complementary PCA showing the indicator strength is also shown (Figure 4B). An extensive list of K-means clustering for each impact degree category is presented in Figure 5.

5. Discussion

5.1. Identifying Local/Microregional and Regional Anthropogenic Pollen Indicators

The anthropogenic indicators for the Balearic Islands can be identified based on results obtained from redundancy analysis and Pearson correlation results (Figure 3 and Table 1). The correlation between pollen types and three variables (landscape openness, pasture size and agropastoral use) from our study is especially useful for proposing anthropogenic pollen indicators. Interestingly, Plantago lanceolata-t, Avena/Triticum-group, Rubus, Galium-t, Caryophyllaceae undiff., Centaurea, Echium and Plantago major/media are the pollen types with the best performance as agropastoral indicators, as shown by the moderate/high correlation with “agropastoral use”. Other taxa, such as Poaceae undiff., Polygonum aviculare-t, Matricaria-t, Papaveraceae undiff, Dipsacaceae, Vitis, Trifolium-t, Plantago sp., Cirsium-t, Urtica dioica-t, Senecio-t and Cichorieae, are more correlated with large pastures and the degree of landscape openness.
Other heliophilous taxa such as Verbascum-t, Plantago cf. crassifolia, Plantago coronopus-t and Chenopodiaceae are correlated to high values of the distance to forest variable, suggesting their preference for open environments. These taxa are especially abundant in samples from coastal, salt-tolerant and dune vegetation in our study, where salt and soil compaction are frequent. In this sense, while not statistically significant, trampling appears correlated to these taxa. Plantago crassifolia pollen is included within the P. coronopus-t [52]. As the P. crassifolia plant mostly grows as a halophyte on shore marshlands [53,54], Plantago coronopus-t must be interpreted carefully as an API in coastal paleoenvironmental records [18,22]. This is the same for Chenopodiaceae. Species from this family are salt-tolerant but they are also frequent in ruderal and nitrophilous environments. Nevertheless, its overrepresentation in the modern pollen rain from the Balearic Islands hinders its consideration as an API, when not accompanied by other APIs in pollen assemblages.
A detailed list of potential regional and local indicators is shown in Table 1 and Figure 4. This list is based on the fidelity/dispersibility and Davis indices obtained in Servera-Vives et al. [18] and the statistical work presented in this work. Potential local indicators of open and anthropized habitats include Poaceae undiff., Polygonum aviculare-t, Avena/Triticum-group, Matricaria-t and Urtica dioica-t. Otherwise, pollen types identified as potential regional indicators include Plantago lanceolata-t, Rubus, Galium-t, Caryophyllaceae undiff., Centaurea, Echium and Plantago major/media. Modern analogue research carried out in mountain areas such as the Pyrenees also highlighted that Plantago lanceolata-t and Plantago major/media indicate regional human activity [4]. Even though Poaceae are anemophilous taxa, the high association values in the Davis indices suggest a good plant–pollen association, as has been recorded in similar research [4,6,55,56].

5.2. Pollen Types Related to Different Degrees of Human Impact Intensity

Despite the many studies on modern pollen analogues that have been carried out to understand how pollen types respond to specific human practices, attempts to specify pollen’s role in differentiating the degree of human impacts has not yet been explored systematically. Consequently, we explored the potential of pollen indicators to evaluate the degree of human impacts based on an ad hoc AIS index and statistical analysis on modern pollen data from the Balearic Islands (Figure 4 and Figure 6). This work provides insights into comprehensive pollen indicators that may also be applied to other Mediterranean island environments.
The group of indicators most related to environments with low anthropogenic impacts aggregates several Ericaceae pollen types (Erica sp., Erica arborea-t, Erica cf. multiflora), as well as Cerastium-t, Medicago and Rubia. Ericaceae pollen curves have regularly been related to early successional vegetation and high fire frequency in pollen analysis from the Mediterranean region [22,57,58]. This may be due to the adaptation of many species of this family to regular fire disturbance through developing lignotubers that allow resprouting after burning [59,60]. In the Balearic Islands, all three Erica species are currently present in the landscape, including E. multiflora L., E. arborea L. and E. scoparia L., present lignotubers and are common in sclerophyllous maquis and garrigues formations [59,60,61]. Erica pollen may be considered a local pollen indicator (associated taxa) with low dispersibility and strong fidelity in the Balearic Islands [18]. Moderate to high values of Erica pollen types can be associated with repeated burning for pastoral activities [58].
Moderate human impact indicators include Cistus pollen types (C. salvifolius and C. monspeliensis), several Plantago types (P. cf. crassifolia, P. bellardii-t and P. albicans-t), several Asteroideae types (Centaurea, C. jacea, Carduus-t, Cirsium-t and Matricaria-t), Verbascum-t, Phlomis and Sanguisorba minor. Centaurea and Plantago types are among the most frequent taxa recorded in Italian archaeopalynological studies, proving their value as cultural pollen indicators [7], and they are both interpreted as pastoral indicator with increased prevalence in grazed areas [5,62]. This study suggests that Sinapis-t and Matricaria-t are also likely important agropastoral indicators that should be considered in addition to Centaurea and Plantago types. Balearic Verbascum species are recurrent in croplands and roadsides, and therefore Verbascum pollen type may be also considered an API based on our research. Matricaria-t includes five plant species recorded during the Braun-Blanquet descriptions, among which three develop in ruderal and agricultural environs, namely Anthemis arvensis L., Chrysanthemum segetum L. and Glebionis coronaria (L.) Cass ex Sprach. In previous modern analogue research from the Balearic Islands [18], Matricaria-t pollen values higher than 1% were only found in crop fields and pastures. Similarly, Sinapis-t has also previously been considered an indicator of weed flora in other modern analogue studies [5], and its role as a cultural indicator is confirmed by the results of both RDA and PCA analyses presented here.
Indicators of high human impact include 12 pollen types: Papaveraceae undiff., Dipsacaceae, Secale-t, Polygonum aviculare-t, Trifolium-t, Centaurea nigra-t, Cerealia SUM, Rumex, Plantago coronopus-t, Plantago SUM, Lotus-t and Brassicaceae SUM. All these pollen types are frequently referred to as cultural pollen indicators in modern pollen analogue research (e.g, [5,6,7]. The family of Dipsacaceae includes several synanthropic taxa in the Balearic Islands such as Centranthus calcitrapae (L.) Dufresne and Scabiosa atropurpurea L., both of which were reported in the Braun-Blanquet descriptions collected during the fieldwork. Centranthus and Dipsacaceae pollen were reported with greater values in sites categorized as high impact in the AIS classification. Although Dipsacaceae is not usually considered an anthropogenic pollen indicator, RDA and PCA results suggest it is most related to sites with high AIS values and that it should be considered an indicator of high levels of anthropization. Moreover, Dipsacaceae pollen has also been related to human practices in other palynological research [62,63,64].
Brassicaceae is also an indicator of high human impact in this study and are commonly related to human-induced environments in pollen-analytic studies [65], but the highest values are related to cultivated fields in modern analogue research [62]. Interestingly, Brassicaceae pollen types are found in high abundance in prehistoric settlements and funerary archaeological sites from the Balearics [21,66], suggesting that humans interacted with Brassicaceae plants through both propagation and cultural uses due to their economic, dietary and symbolic value. In the historical period, Brassicaceae have been found in stratigraphies from the ancient imperial harbors of Rome [67,68] and Naples [69], as well as in rural contexts in Sicily [70]. These taxa have been interpreted as evidence for crop cultivation of cabbage during the Roman period, as attested to by classical authors (e.g., Pliny the Elder) or for oil extraction, as is documented by recovered macroremains [71].
Cereal pollen types are confirmed to be good indicators of both human impact and agropastoral uses by ordination techniques. Secale-t and Cerealia SUM (Avena/Triticum-group + Secale-t + Cerealia-t) are related to high impact, but Avena/Triticum-group is included in the intermediate indicator category. Nevertheless, its scores in the two first axes of the constrained PCA were between high and moderate impact, suggesting that this pollen type may potentially be related to both moderate and high impact sites. In this sense, the RDA logically confirms Avena/Triticum-group as a good indicator of agropastoral activities. Both Plantago SUM and P. coronopus pollen types are included in the high impact category. As explained previously, Plantago cf. crassifolia pollen is included within the P. coronopus type, being the first classed as a moderate impact indicator, therefore pointing out that these pollen taxonomical classifications may be associated with cultural activities, especially when they are associated with other anthropogenic pollen indicators.
Finally, the intermediate impact category includes pollen types that are difficult to assign to a specific degree of AIS. This is likely caused by the ubiquity of such taxa in all the studied sites. This should not be understood as their lack of value as anthropogenic pollen indicators, but rather that they should be considered general indicators of a variety of anthropogenic activities. In this sense, some pollen types clearly associated with agropastoral use in the RDA, such as Urtica dioica-t or Plantago lanceolata-t, do not perform well as reliable indicators for detecting human impact gradients. The intermediate group includes well known API based both on palynological and ecological knowledge. For instance, Plantago lanceolata-t, Avena/Triticum-t and Cichorieae are included in this general category, but much pollen-analysis research has proved their strong value in detecting human impact [5,8]. Plantago lanceolata-t has commonly been used as an anthropogenic pollen indicator in several palynological analyses, especially in North European areas [2] while the interpretations of low values in Mediterranean and circum-Mediterranean areas should be considered carefully [4,6]. Plantago species are found in ruderal habitats subjected to soil compaction, therefore being related to human and herd presence [7].

5.3. Reflection on Pollen Morphology Resolution and Productive Cultural Practices

In the palynological interpretations, the co-occurrence of several anthropogenic indicators is needed to detect human impacts on vegetation. Indeed, our interpretations should rely on the accurate integration of off-site bio-indicators (e.g., pollen, non-pollen palynomorphs, charcoals, diatoms, etc.) with archaeobotanical, historical and archaeological data. Disregarding such important sources of information would diminish our ability to understand the inter-relationship between landscape processes and human behavior through time. Nevertheless, pollen morphology is a complex task and occasionally it is not possible to achieve the high-resolution level needed to link pollen types to specific environments or human practices [5,9,62]. This is particularly true in archaeopalynological research, where taphonomy and post-depositional processes amplify the increased deterioration of some palynomorphs, resulting in the frequent morphological identification limited to genus or family-level morphological identifications.
Nevertheless, despite the fact that the species level would always be suitable for palaeoecological reconstruction, this does not mean that lower taxonomical levels (family, subfamily, genus, etc.) cannot be valuable anthropogenic pollen indicators. In this study, we include some morphological aggregations such as Brassicaceae or Plantago SUM, which perform even better than the pollen type taxonomic level in the ordination analysis. Another taxon performing well as API is Cichorieae, a tribe of the sub-family Cichorioideae included in the Asteraceae family. Cichorieae has proven to be an indicator of secondary pastures and some types of primary open habitats in the Mediterranean area [41]. In the constrained PCA and k-means clustering, it is included in the intermediate category, but based on the scores on the first PCA axis it seems to be related both to medium and high human impact. Furthermore, Cichorieae in the RDA is clearly linked to larger pastures and agropastoral use. Caryophyllaceae undiff. is also linked to agropastoral use, but not manifestly linked to a determinate degree of landscape impact. Otherwise, other family-level morphological classifications such as Ranunculaceae undiff. do not seem to be clearly associated with either human impact gradients or specific environmental variables in the ordination analysis.

6. Conclusions

This paper constitutes the first proposal of anthropogenic pollen indicators for the Balearic Islands, and a new attempt to assess the human-impact gradient in Mediterranean islands. To do so, we used a combination of modern pollen analogues, phytosociological descriptions and ordination techniques employing quantitative and presence/absence data. Redundancy analysis allowed us to evaluate the relation of pollen types to significant environmental variables, from which we devised an anthropogenic index score (AIS) for each site based on the degree to which it exhibited the most impactful anthropogenic environmental variables.
Comparing the redundancy analysis data and the pollen representation indices (Davis’ and fidelity dispersibility indices) used in previous research [18], we proposed local/microregional and regional anthropogenic indicators that will allow us to further spatialize human activities during the Holocene. Local/microregional indicators of human activities are Poaceae undiff., Polygonum aviculare-t, Avena/Triticum-group, Cerealia SUM, Matricaria-t and Urtica dioica-t, while regional indicators are Plantago lanceolata-t, Rubus, Galium-t, Caryophyllaceae undiff., Centaurea, Echium, Plantago major/media among others. Additionally, we identified indicators of low anthropogenic impact (i.e., Erica, Cerastium-t and Medicago); moderate anthropogenic impact, (i.e., Sinapis-t, Sanguisorba minor-t Plantago bellardii-t); and high anthropogenic impact (i.e., Papaveraceae, Dipsacaceae and Secale-t). Some taxa are not clearly related to specific impact gradients, such as Galium-t, Urtica dioica-t. and Cistus, and have been classified as intermediate anthropogenic indicators through this study. Our research also confirmed that lower taxonomical levels in the pollen analysis (genera or family) may also be good indicators of human activities, as for instance with Brassicaceae, Plantago SUM or Cichorieae.
Some potential limitations of our research are that current human activities are not easily and directly comparable to those that took place in the past. Therefore, when constructing general indices, possible biases must be considered. In this sense, quality and resolution of environmental variable datasets can also make the interpretation of pollen type-environmental variable correlations difficult. To tackle this issue, it is mandatory to generate narratives on past landscape dynamics, socioenvironmental interactions and anthropization intensity by putting anthropogenic indicators in a broader context including other palaeoecological, historical and archaeological information.
This research may have implications for the environmental, archaeological and palaeoecological fields. As insights into modern pollen analogues provide new perspectives on pollen–vegetation relationships, this paper extends this understanding by leveraging land-use dynamics as a means to define the degree of anthropogenic impact in these landscapes. This process is especially necessary to assess colonization and anthropization rhythms in Mediterranean island environments. The novel methods applied through this study can also be applied to further interdisciplinary research pertaining to environmental restoration and management that are interested in assessing anthropogenic drivers of change in areas with differing ecological and cultural histories. Forthcoming research will further apply the results presented here to more fully understand Holocene anthropogenic dynamics in the Balearic Islands.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/su15118807/s1, Table S1: Table with the description and source of environmental variables selected for this research; S2-Code: Rmarkdown with the analytical code used in this paper.

Author Contributions

Conceptualization, G.S.-V. and M.M.A.; methodology, G.S.-V., M.M.A. and G.S.; software, G.S. and M.R.; validation, P.T.; formal analysis, G.S.-V., G.S. and M.R.; investigation, G.S.-V., M.M.A. and P.T.; resources, P.T. and A.M.M.; data curation, G.S.-V., A.F. and A.M.M.; writing—original draft preparation, G.S.-V.; writing—review and editing, M.M.A., G.S., A.F., P.T., M.R. and A.M.M.; visualization, G.S.-V.; supervision, G.S.-V., M.M.A., A.F. and A.M.M.; project administration, A.M.M.; funding acquisition, G.S.-V. All authors have read and agreed to the published version of the manuscript.

Funding

This project received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie [grant agreement 895735; OLEAproject); Movilidad y conectividad de las comunidades prehistóricas en el Mediterraneo occidental durante la prehistoria reciente: El caso de las Islas Baleares [PID2019-108692GBI00], which is financed by the Spanish government’s Ministry of Science, Innovation and Universities MCIN/AEI/10.13039/501100011033.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Specific data and code used for this paper can be found as RMarkdown in the Supplementary Material (S2). The whole pollen counts dataset can be found in the NEOTOMA specialized repository.

Acknowledgments

The authors thank Maria Bofill Martínez and the Cabrera National Park for their support during the fieldwork, as well as Alejandra Galmés and Joan Estrany for their technical support. Additionally, special thanks to Marta Mazzanti for her help with pollen morphology doubts. The authors also thank Max Fancourt for his contribution with drone photography, and Sath (Sandro Thomàs Soler) for his collaboration with the artistic figures.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gaillard, M.J. Pollen Methods and Studies: Archaeological Applications. In Encyclopedia of Quaternary Sciences; Elias, S.A., Mock, C.J., Eds.; Elsevier: Amsterdam, The Netherlands, 2013; Volume 3, pp. 880–9903. [Google Scholar]
  2. Behre, K.-E. The Interpretation of Anthropogenic Indicators in Pollen Diagrams. Pollen Et Spores 1981, 23, 225–245. [Google Scholar]
  3. Gauthier, E. Pollen et Microfossiles: Des Crêts Du Jura Aux Fjords Du Groenland; Habilitation à Diriger des Recherches, Université de Franche-Comté: Besançon, France, 2012. [Google Scholar]
  4. Mazier, F.; Galop, D.; Brun, C.; Buttler, A. Modern Pollen Assemblages from Grazed Vegetation in the Western Pyrenees, France: A Numerical Tool for More Precise Reconstruction of Past Cultural Landscapes. Holocene 2006, 16, 91–103. [Google Scholar] [CrossRef]
  5. Brun, C.; Dessaint, F.; Richard, H.; Bretagnolle, F. Arable-Weed Flora and Its Pollen Representation: A Case Study from the Eastern Part of France. Rev. Palaeobot. Palynol. 2007, 146, 29–50. [Google Scholar] [CrossRef]
  6. Ejarque, A.; Miras, Y.; Riera, S. Pollen and Non-Pollen Palynomorph Indicators of Vegetation and Highland Grazing Activities Obtained from Modern Surface and Dung Datasets in the Eastern Pyrenees. Rev. Palaeobot. Palynol. 2011, 167, 123–139. [Google Scholar] [CrossRef]
  7. Mercuri, A.M.; Bandini, M.M.; Florenzano, A.; Montecchi, M.C.; Rattighieri, E.; Torri, P. Anthropogenic Pollen Indicators (API) from Archaeological Sites as Local Evidence of Human-Induced Environments in the Italian Peninsula. Ann. Di Bot. 2013, 3, 143–153. [Google Scholar] [CrossRef]
  8. Deza-Araujo, M.; Morales-Molino, C.; Conedera, M.; Henne, P.D.; Krebs, P.; Hinz, M.; Heitz, C.; Hafner, A.; Tinner, W. A New Indicator Approach to Reconstruct Agricultural Land Use in Europe from Sedimentary Pollen Assemblages. Palaeogeogr. Palaeoclim. Palaeoecol. 2022, 599, 111051. [Google Scholar] [CrossRef]
  9. Deza-Araujo, M.; Morales-Molino, C.; Tinner, W.; Henne, P.D.; Heitz, C.; Pezzatti, G.B.; Hafner, A.; Conedera, M. A Critical Assessment of Human-Impact Indices Based on Anthropogenic Pollen Indicators. Quat. Sci. Rev. 2020, 236, 106291. [Google Scholar] [CrossRef]
  10. Medail, F.; Quezel, P. Biodiversity Hotspots in the Mediterranean Basin: Setting Global Conservation Priorities. Conserv. Biol. 1999, 13, 1510–1513. [Google Scholar] [CrossRef]
  11. Blondel, J. On Humans and Wildlife in Mediterranean Islands. J. Biogeogr. 2008, 35, 509–518. [Google Scholar] [CrossRef]
  12. Grove, A.T.; Rackham, O. The Nature of Mediterranean Europe. An Ecological History; Yale University Press: New Haven, CT, USA; London, UK, 2001. [Google Scholar]
  13. Woodbridge, J.; Roberts, N.; Fyfe, R. Pan-Mediterranean Holocene Vegetation and Land-Cover Dynamics from Synthesized Pollen Data. J. Biogeogr. 2018, 45, 2159–2174. [Google Scholar] [CrossRef]
  14. Mercuri, A.M.; Florenzano, A.; Burjachs, F.; Giardini, M.; Kouli, K.; Masi, A.; Picornell-Gelabert, L.; Revelles, J.; Sadori, L.; Servera-Vives, G.; et al. From Influence to Impact: The Multifunctional Land Use in Mediterranean Prehistory Emerging from Palynology of Archaeological Sites (8.0-2.8 Ka BP). Holocene 2019, 29, 830–846. [Google Scholar] [CrossRef]
  15. Murray, I.; Jover-Avellà, G.; Fullana, O.; Tello, E. Biocultural Heritages in Mallorca: Explaining the Resilience of Peasant Landscapes within a Mediterranean Tourist Hotspot, 1870–2016. Sustainability 2019, 11, 1926. [Google Scholar] [CrossRef]
  16. Dawson, H. Island Archaeology. In Encyclopedia of Global Archaeology; Springer: Cham, Switzerland, 2019; pp. 1–8. [Google Scholar] [CrossRef]
  17. Nogué, S.; Santos, A.M.C.; John, H.; Björck, S.; Castilla-Beltrán, A.; Connor, S.; de Boer, E.J.; de Nascimento, L.; Felde, V.A.; Fernández-Palacios, J.M.; et al. The Human Dimension of Biodiversity Changes on Islands. Science 2021, 372, 488–491. [Google Scholar] [CrossRef]
  18. Servera-Vives, G.; Mus Amezquita, M.; Snitker, G.; Florenzano, A.; Torri, P.; Estrany Bertos, J.; Mercuri, A.M. Modern Analogs for Understanding Pollen-Vegetation Dynamics in a Mediterranean Mosaic Landscape (Balearic Islands, Western Mediterranean). Holocene 2022, 32, 716–734. [Google Scholar] [CrossRef]
  19. Morey, M.; Ruiz-Pérez, M. The Balearic Islands. In Mediterranean Island Landscapes: Natural and Cultural Approaches; Vogiatzakis, I., Pungetti, G., Mannion, A.M., Eds.; Springer: Dordrecht, The Netherlands, 2008; pp. 271–296. ISBN 978-1-4020-5064-0. [Google Scholar]
  20. Guijarro, J.A. Contribución a La Bioclimatología de Baleares; Universitat de les Illes Balears: Palma, Spain, 1986. [Google Scholar]
  21. Picornell-Gelabert, L.; Servera-Vives, G. Landscape Practices and Everyday Life in Domestic Spaces in Bronze Age Mallorca (Balearic Islands): Perspectives for and Archaeology of Fuel and Firewood. Quat. Int. 2017, 431, 73–89. [Google Scholar] [CrossRef]
  22. Servera-Vives, G.; Riera, S.; Picornell-Gelabert, L.; Moffa-Sánchez, P.; Llergo, Y.; Garcia, A.; Mus-Amezquita, M.; García Álvarez, S.; Calvo Trías, M. The Onset of Islandscapes in the Balearic Islands: A Study-Case of Addaia (Northern Minorca, Spain). Palaeogeogr. Palaeoclim. Palaeoecol. 2018, 498, 9–23. [Google Scholar] [CrossRef]
  23. Llorens, L.; Gil, L. The Balearic Islands. In The Vegetation of the Iberian Peninsula; Loidi, J., Ed.; Series “Plant and Vegetation”; Springer: Cham, Switzerland, 2017; Volume 2, pp. 3–33. ISBN 978-3-319-54866-1. [Google Scholar]
  24. Pons Esteva, A.; Rullán Salamanca, O. Artificialization and Islandness on the Spanish Tourist Coast. Misc. Geogr. 2014, 18, 5–16. [Google Scholar] [CrossRef]
  25. Capó, M.; Roig-Oliver, M.; Cardona, C.; Cursach, J.; Bartolomé, J.; Rita, J.; Baraza, E. Historic Exposure to Herbivores, Not Constitutive Traits, Explains Plant Tolerance to Herbivory in the Case of Two Medicago Species (Fabaceae). Plant Sci. 2021, 307, 110890. [Google Scholar] [CrossRef]
  26. Dante Cerrato, M.; Ribas-Serra, A.; Miquel Mir-Rosselló, P.; Cardona Ametller, C.; Gil-Vives, L. Time Pattern Variation of Alien Plant Introductions in an Insular Biodiversity Hotspot: The Balearic Islands as a Case Study for the Mediterranean Region. Biodivers. Conserv. 2023. [Google Scholar] [CrossRef]
  27. Podda, L.; Fraga, I.; Arguimbau, P.; Mascia, F.; Mayoral García-Berlanga, O.; Bacchetta, G. Comparison of the Invasive Alien Flora in Continental Islands: Sardinia (Italy) and Balearic Islands (Spain). Rend. Lincei 2011, 22, 31–45. [Google Scholar] [CrossRef]
  28. Moragues, E.; Rita, J. Espècies Introduïdes Balears. In Conselleria de Medi Ambient; Govern de les Illes Balears: Palma, Spain, 2005; Volume 11. [Google Scholar]
  29. Räsänen, S. Tracing and Interpreting Fine-Scale Human Impact in Northern Fennoscandia with the Aid of Modern Pollen Analogues. Veg. Hist. Archaeobotany 2001, 10, 211–218. [Google Scholar] [CrossRef]
  30. Räsänen, S.; Suutari, H.; Nielsen, A.B. A Step Further towards Quantitative Reconstruction of Past Vegetation in Fennoscandian Boreal Forests: Pollen Productivity Estimates for Six Dominant Taxa. Rev. Palaeobot. Palynol. 2007, 146, 208–220. [Google Scholar] [CrossRef]
  31. Florenzano, A.; Mercuri, A.M.; Rinaldi, R.; Rattighieri, E.; Fornaciari, R.; Messora, R.; Arru, L. The Representativeness of Olea Pollen Fromolive Groves and the Late Holocene Landscape Reconstruction in Central Mediterranean. Front. Earth Sci. 2017, 5, 85. [Google Scholar] [CrossRef]
  32. Braun-Blanquet, J. Fitosociologia Bases Para El Estudio de Las Comunidades Vegetales; Blume: Madrid, Spain, 1979; ISBN 8472141748/9788472141742. [Google Scholar]
  33. Erdtman, G. The Acetolysis Method—A Revised Description. Vensk. Bot. Tidskr. 1960, 54, 561–564. [Google Scholar]
  34. Faegri, K.; Kaland, P.E.; Krzywinski, K. Textbook of Pollen Analysis, 4th ed.; The Blackburn Press: Caldwell, NJ, USA, 1989; ISBN 978-1930665019. [Google Scholar]
  35. Beug, H.J. Leitfaden der Pollenbestimmung Für Mitteleuropa und Angrenzende Gebiete; Verlag Dr. Friedrich Pfeil: München, Germany, 2004; ISBN 3-89937-043-0. [Google Scholar]
  36. Reille, M. Pollen et Spores d’Europe et d’Afrique Du Nord. In Atlas Photographique; Laboratoire de Botanique Historique et Palynologie, URA 1152/CNRS: Marseille, France, 1992. [Google Scholar]
  37. Chester, P.I.P.I.; Raine, J.I.; Raine, I. Pollen and Spore Keys for Quaternary Deposits in the Northern Pindos Mountains, Greece. Grana 2001, 40, 299–387. [Google Scholar] [CrossRef]
  38. Punt, W.; Hoen, P.P. The Northwest European Pollen Flora, 70. Asteraceae—Asteroideae. Rev. Palaeobot. Palynol. 2009, 157, 22–183. [Google Scholar] [CrossRef]
  39. Tweddle, J.C.; Edwards, K.J.; Fieller, N.R.J.J. Multivariate Statistical and Other Approaches for the Separation of Cereal from Wild Poaceae Pollen Using a Large Holocene Dataset. Veg. Hist. Archaeobot. 2005, 14, 15–30. [Google Scholar] [CrossRef]
  40. Kouli, K. Plant Landscape and Land Use at the Neolithic Lake Settlement of Dispilió (Macedonia, Northern Greece). Plant Biosyst. Int. J. Deal. All Asp. Plant Biol. 2015, 149, 195–204. [Google Scholar] [CrossRef]
  41. Florenzano, A.; Marignani, M.; Rosati, L.; Fascetti, S.; Mercuri, A.M. Are Cichorieae an Indicator of Open Habitats and Pastoralism in Current and Past Vegetation Studies? Plant Biosyst. Int. J. Deal. All Asp. Plant Biol. 2015, 149, 154–165. [Google Scholar] [CrossRef]
  42. López-Sáez, J.A.; Glais, A.; Tsiftsis, S.; Lezpez, L. Modern Pollen–Vegetation Relationships along an Altitudinal Transect in the Lefka Ori Massif (Western Crete, Greece). Rev. Palaeobot. Palynol. 2018, 259, 159–170. [Google Scholar] [CrossRef]
  43. R Core Team. R: A Language and Environment for Statistical Computing, version 4.2.0; R Foundation for Statistical Computing; R Core Team: Viena, Austria, 2022.
  44. Oksanen, J.; Legendre, P.; O’Hara, B.; Stevens, M.H.H.; Oksanen, M.J.; Suggests, M. Vegan: Community Ecology Package. R Package, version 2.6-4; Community Ecology Package; DataCamp: New York, NY, USA, 2020.
  45. Wickham, H. Ggplot2: Elegant Graphics for Data Analysis; Springer International Publishing: Berlin/Heidelberg, Germany, 2016; ISBN 9783319242750. [Google Scholar]
  46. Smilauer, P.; Lepš, J. Multivariate Analysis of Ecological Data Using CANOCO 5; Cambridge University Press: Cambridge, UK, 2014; ISBN 9781139627061. [Google Scholar]
  47. López-Sáez, J.A.; Camarero, J.J.; Abel-Schaad, D.; Luelmo-Lautenschlaeger, R.; Pérez-Díaz, S.; Alba-Sánchez, F.; Carrión, J.S. Don’t Lose Sight of the Forest for the Trees! Discerning Iberian Pine Communities by Means of Pollen-Vegetation Relationships. Rev. Palaeobot. Palynol. 2020, 281, 104285. [Google Scholar] [CrossRef]
  48. López Sáez, J.; Burjachs i Casas, F.; López García, P.; Sáez, L.; García, L. Arqueopalinología: Síntesis Crítica. Polen 2003, 35, 5–35. [Google Scholar]
  49. Davis, O.K. Pollen Frequencies Reflect Vegetation Patterns in a Great Basin (USA) Mountain Range. Rev. Palaeobot. Palynol. 1984, 40, 295–315. [Google Scholar] [CrossRef]
  50. Fall, P.L. Modern Vegetation, Pollen and Climate Relationships on the Mediterranean Island of Cyprus. Rev. Palaeobot. Palynol. 2012, 185, 79–92. [Google Scholar] [CrossRef]
  51. McGlone, M.S.; Meurk, C.D. Modern Pollen Rain, Subantarctic Campbell Island, New Zealand. N. Z. J. Ecol. 2000, 24, 181–194. [Google Scholar]
  52. Ubera, J.L.; Galán, C.; Guerrero, F.H. Palynological Study of the Genus Plantago in the Iberian Peninsula. Grana 1988, 27, 1–15. [Google Scholar] [CrossRef]
  53. Llorens, L.; Gil, L.; Tébar, F.J. La Vegetació de l’illa de Mallorca. Bases per a La Definició i Gestió Dels Hàbitats. In Conselleria de Medi Ambient; Govern de les Illes Balears: Palma, Spain, 2007; ISBN 978-84-612-0488-5. [Google Scholar]
  54. Vicente, O.; Boscaiu, M.; Naranjo, M.A.Á.; Estrelles, E.; Bellés, J.M.M.; Soriano, P. Responses to Salt Stress in the Halophyte “Plantago crassifolia” (Plantaginaceae). J. Arid. Environ. 2004, 58, 463–481. [Google Scholar] [CrossRef]
  55. Cañellas-Boltà, N.; Rull, V.; Vigo, J.; Mercadé, A. Modern Pollen-Vegetation Relationships along an Altitudinal Transect in the Central Pyrenees (Southwestern Europe). Holocene 2009, 19, 1185–1200. [Google Scholar] [CrossRef]
  56. Court-Picon, M.; Buttler, A.; De Beaulieu, J.L. Modern Pollen-Vegetation Relationships in the Champsaur Valley (French Alps) and Their Potential in the Interpretation of Fossil Pollen Records of Past Cultural Landscapes. Rev. Palaeobot. Palynol. 2005, 135, 13–39. [Google Scholar] [CrossRef]
  57. Beffa, G.; Pedrotta, T.; Colombaroli, D.; Henne, P.D.; van Leeuwen, J.F.N.; Süsstrunk, P.; Kaltenrieder, P.; Adolf, C.; Vogel, H.; Pasta, S.; et al. Vegetation and Fire History of Coastal North-Eastern Sardinia (Italy) under Changing Holocene Climates and Land Use. Veg. Hist. Archaeobot. 2016, 25, 271–289. [Google Scholar] [CrossRef]
  58. Pedrotta, T.; Gobet, E.; Schwörer, C.; Beffa, G.; Butz, C.; Henne, P.D.; Morales-Molino, C.; Pasta, S.; van Leeuwen, J.F.N.; Vogel, H.; et al. 8000 Years of Climate, Vegetation, Fire and Land-Use Dynamics in the Thermo-Mediterranean Vegetation Belt of Northern Sardinia (Italy). Veg. Hist. Archaeobot. 2021, 30, 789–813. [Google Scholar] [CrossRef] [PubMed]
  59. Paula, S.; Naulin, P.I.; Arce, C.; Galaz, C.; Pausas, J.G. Lignotubers in Mediterranean Basin Plants. Plant Ecol. 2016, 217, 661–676. [Google Scholar] [CrossRef]
  60. Pausas, J.G.; Lamont, B.B.; Paula, S.; Appezzato-da-Glória, B.; Fidelis, A. Unearthing Belowground Bud Banks in Fire-Prone Ecosystems. New Phytol. 2018, 217, 1435–1448. [Google Scholar] [CrossRef]
  61. Paula, S.; Ojeda, F. Resistance of Three Co-Occurring Resprouter Erica Species to Highly Frequent Disturbance. Plant Ecol. 2006, 183, 329–336. [Google Scholar] [CrossRef]
  62. Court-Picon, M.; Buttler, A.; De Beaulieu, J.L. Modern Pollen/Vegetation/Land-Use Relationships in Mountain Environments: An Example from the Champsaur Valley (French Alps). Veg. Hist. Archaeobot. 2006, 15, 151–168. [Google Scholar] [CrossRef]
  63. Hjelle, K.L. Modem Pollen Assemblages from Mown and Grazed Vegetation Types in Western Norway. Rev. Palaeobot. Palynol. 1999, 107, 55–81. [Google Scholar] [CrossRef]
  64. Marinova, E.; Atanassova, J. Anthropogenic Impact on Vegetation and Environment during the Bronze Age in the Area of Lake Durankulak, NE Bulgaria: Pollen, Microscopic Charcoal, Non-Pollen Palynomorphs and Plant Macrofossils. Rev. Palaeobot. Palynol. 2006, 141, 165–178. [Google Scholar] [CrossRef]
  65. Behre, K.-E. The rôle of man in European vegetation history. In Vegetation History. Handbook of Vegetation Science; Huntley, B., Webb, T., Eds.; Springer: Dordrecht, The Netherlands, 1988; Volume 7. [Google Scholar] [CrossRef]
  66. Riera Mora, S.; Servera-Vives, G.; Picornell-Gelabert, L.; Cabanis, M.; Boi, M.; Miras, Y. Pollen Signatures of a Ritual Process in the Collective Burial Cave of Cova des Pas (Late Bronze Age, Minorca, Balearic Islands, Spain). In The Bioarchaeology of Ritual and Religion; Livarda, A., Madgwick, R., Santiago, R., Eds.; Oxbow: Oxford, UK, 2018; pp. 28–43. ISBN 978-1-78570-828-2. [Google Scholar]
  67. Sadori, L.; Allevato, E.; Bellini, C.; Bertacchi, A.; Boetto, G.; Di Pasquale, G.; Giachi, G.; Giardini, M.; Masi, A.; Pepe, C.; et al. Archaeobotany in Italian Ancient Roman Harbours. Rev. Palaeobot. Palynol. 2015, 218, 217–230. [Google Scholar] [CrossRef]
  68. Sadori, L.; Giardini, M.; Giraudi, C.; Mazzini, I. The Plant Landscape of the Imperial Harbour of Rome. J. Archaeol. Sci. 2010, 37, 3294–3305. [Google Scholar] [CrossRef]
  69. Russo Ermolli, E.; Romano, P.; Ruello, M.R.; Barone Lumaga, M.R. The Natural and Cultural Landscape of Naples (Southern Italy) during the Graeco-Roman and Late Antique Periods. J. Archaeol. Sci. 2014, 42, 399–411. [Google Scholar] [CrossRef]
  70. Montecchi, M.C.; Mercuri, A.M. When Palynology Meets Classical Archaeology: The Roman and Medieval Landscapes at the Villa Del Casale Di Piazza Armerina, UNESCO Site in Sicily. Archaeol. Anthropol. Sci. 2016, 10, 743–757. [Google Scholar] [CrossRef]
  71. Bandini Mazzanti, M.; Bosi, G.; Mercuri, A.M.; Accorsi, C.A.; Guarnieri, C. Plant Use in a City in Northern Italy during the Late Medieval and Renaissance Periods: Results of the Archaeobotanical Investigation of “The Mirror Pit” (14th–15th Century AD) in Ferrara. Veg. Hist. Archaeobotany 2005, 14, 442–452. [Google Scholar] [CrossRef]
Sustainability 15 08807 g001 550
Figure 1. Map showing sampling locations of modern pollen samples with the estimated anthropogenic rank score (discussed in Section 3). Vegetation categories are based on Corine Land Cover 2018 classes.
Figure 1. Map showing sampling locations of modern pollen samples with the estimated anthropogenic rank score (discussed in Section 3). Vegetation categories are based on Corine Land Cover 2018 classes.
Sustainability 15 08807 g001
Sustainability 15 08807 g002 550
Figure 2. Example of several landscape types from the Balearic Islands: (A) xeric thyme-garrigue coastal vegetation (frigolar) from the Pytiusic Islands; (B) holm oak formations from the Serra de Tramuntana (Mallorca); (C) Rumex pasture in a ravine wetland from Barranc de Trebalúger (Menorca); (D) mountain olive orchards from the Serra de Tramuntana (Mallorca); (E) aerial view of a mosaic landscape with macchia and crop fields (north-eastern Mallorca); (F) aerial view of extensive agricultural fields from the interior of Mallorca.
Figure 2. Example of several landscape types from the Balearic Islands: (A) xeric thyme-garrigue coastal vegetation (frigolar) from the Pytiusic Islands; (B) holm oak formations from the Serra de Tramuntana (Mallorca); (C) Rumex pasture in a ravine wetland from Barranc de Trebalúger (Menorca); (D) mountain olive orchards from the Serra de Tramuntana (Mallorca); (E) aerial view of a mosaic landscape with macchia and crop fields (north-eastern Mallorca); (F) aerial view of extensive agricultural fields from the interior of Mallorca.
Sustainability 15 08807 g002
Sustainability 15 08807 g003 550
Figure 3. Redundancy analysis results showing the correlation between environmental variables with sites (left) and between environmental variables with the top 50% pollen types (right). Color-shaded quadrants refers to main interpretative groups.
Figure 3. Redundancy analysis results showing the correlation between environmental variables with sites (left) and between environmental variables with the top 50% pollen types (right). Color-shaded quadrants refers to main interpretative groups.
Sustainability 15 08807 g003
Sustainability 15 08807 g004 550
Figure 4. (A) Constrained PCA showing the relation between the human impact gradient and the assigned pollen types to each category. (B) Constrained PCA based on the indicator strength index for each pollen type.
Figure 4. (A) Constrained PCA showing the relation between the human impact gradient and the assigned pollen types to each category. (B) Constrained PCA based on the indicator strength index for each pollen type.
Sustainability 15 08807 g004
Sustainability 15 08807 g005 550
Figure 5. Indicator strength results for pollen types within the correspondent human-impact gradient category. Orange-colored labels correspond to local/microregional agropastoral indicators in the Pearson correlation (Table 1), while dark cyan represent regional agropastoral indicators.
Figure 5. Indicator strength results for pollen types within the correspondent human-impact gradient category. Orange-colored labels correspond to local/microregional agropastoral indicators in the Pearson correlation (Table 1), while dark cyan represent regional agropastoral indicators.
Sustainability 15 08807 g005
Sustainability 15 08807 g006 550
Figure 6. Artistic illustration showing selected pollen types associated with low, moderate, and high human impact degree. An exhaustive list of pollen types for the anthropogenic impact gradient appears in Figure 5.
Figure 6. Artistic illustration showing selected pollen types associated with low, moderate, and high human impact degree. An exhaustive list of pollen types for the anthropogenic impact gradient appears in Figure 5.
Sustainability 15 08807 g006
Table
Table 1. List of taxa showing positive Pearson correlation r values higher than 0.5 to significant environmental variables. The list of pollen types is listed in decreasing correlation order for each category, while the order of environmental variables follows the RDA p-values. Local indicators vs. regional are proposed either considering pollen types classed as SAT or AT in the Davis indices calculation or high fidelity (>50) + low dispersibility (<50) values resulting from Servera-Vives et al. [18]. Taxa linked to the agropastoral use variable are highlighted in bold in the text.
Table 1. List of taxa showing positive Pearson correlation r values higher than 0.5 to significant environmental variables. The list of pollen types is listed in decreasing correlation order for each category, while the order of environmental variables follows the RDA p-values. Local indicators vs. regional are proposed either considering pollen types classed as SAT or AT in the Davis indices calculation or high fidelity (>50) + low dispersibility (<50) values resulting from Servera-Vives et al. [18]. Taxa linked to the agropastoral use variable are highlighted in bold in the text.
Regional Indicators Local/Microregional Indicators
Distance to forestPlantago cf. crassifolia, Verbascum-t, Hornungia-t, Plantago coronopus-t, Lotus-t, Apiaceae, Plantago bellardii-t, Centaurea jacea-t, Linaria-t, Chenopodiaceae, Brassicaceae SUMCistus salviifolius, Artemisia
Pasture areaPapaveraceae undiff, Dipsacaceae, Vitis, Trifolium-t, Plantago sp., Cirsium-t, Senecio-t, Cichorieae, Brassicaceae SUM, Chenopodiaceae, Rumex, Carduus-t, Apiaceae, Sinapis-t, Lotus-t, Plantago albicans, Plantago coronopus-t, Hornungia-t, Plantago major/media, Echium, CentaureaPoaceae undiff., Polygonum aviculare-t, Cerealia SUM, Matricaria-t, Urtica dioica-t
Landscape opennessTrifolium-t, Vitis, Papaveraceae undiff, Dipsacaceae, Plantago sp., Cirsium-t, Senecio-t, Brassicaceae SUM, Chenopodiaceae, Cichorieae, Apiaceae, Lotus-t, Rumex, Plantago coronopus-t, Carduus-t, Sinapis-t, Hornungia-t, Plantago albicansPoaceae undiff., Polygonum aviculare-t, Cerealia SUM, Matricaria-t
Agropastoral usePlantago lanceolata-t, Rubus, Galium-t, Caryophyllaceae undiff., Centaurea, Echium, Plantago major/mediaAvena/Triticum-group
List of potential local/microregional
indicators of anthropized habitats		Poaceae undiff., Polygonum aviculare-t, Avena/Triticum-group, Cerealia SUM, Matricaria-t, Urtica dioica-t
List of potential regional indicators of anthropized habitatsPlantago lanceolata-t, Rubus, Galium-t, Caryophyllaceae undiff., Centaurea, Echium, Plantago major/media, Papaveraceae undiff, Dipsacaceae, Vitis, Trifolium-t, Plantago sp., Cirsium-t, Senecio-t, Cichorieae, Brassicaceae SUM, Chenopodiaceae, Rumex, Carduus-t, Apiaceae, Sinapis-t, Lotus-t, Plantago albicans, Plantago coronopus-t, Hornungia-t
Table
Table 2. ANOVA permutation test results for the environmental variables selected for the RDA using the forward/backward method for stepwise model selection. Significant p-values (≤0.05) appear in bold. Significance codes: 0.001 = “**”, 0.05 = “*” and 0.05 = “·”
Table 2. ANOVA permutation test results for the environmental variables selected for the RDA using the forward/backward method for stepwise model selection. Significant p-values (≤0.05) appear in bold. Significance codes: 0.001 = “**”, 0.05 = “*” and 0.05 = “·”
Environmental VariableDfVarianceFp-Value
Herbivory pressure10.350.830.603
Agropastoral use10.721.690.057·
Trampling10.501.180.230
Landscape openness10.731.720.036 *
EFFIS Fire Occurrence10.511.210.250
Distance to urban10.280.650.803
Distance to forest11.092.560.009 **
Pasture area10.962.250.021 *
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.

Share and Cite

MDPI and ACS Style

Servera-Vives, G.; Mus Amezquita, M.; Snitker, G.; Florenzano, A.; Torri, P.; Ruiz, M.; Mercuri, A.M. Human-Impact Gradients through Anthropogenic Pollen Indicators in a Mediterranean Mosaic Landscape (Balearic Islands). Sustainability 2023, 15, 8807. https://doi.org/10.3390/su15118807

AMA Style

Servera-Vives G, Mus Amezquita M, Snitker G, Florenzano A, Torri P, Ruiz M, Mercuri AM. Human-Impact Gradients through Anthropogenic Pollen Indicators in a Mediterranean Mosaic Landscape (Balearic Islands). Sustainability. 2023; 15(11):8807. https://doi.org/10.3390/su15118807

Chicago/Turabian Style

Servera-Vives, Gabriel, Maurici Mus Amezquita, Grant Snitker, Assunta Florenzano, Paola Torri, Maurici Ruiz, and Anna Maria Mercuri. 2023. "Human-Impact Gradients through Anthropogenic Pollen Indicators in a Mediterranean Mosaic Landscape (Balearic Islands)" Sustainability 15, no. 11: 8807. https://doi.org/10.3390/su15118807

APA Style

Servera-Vives, G., Mus Amezquita, M., Snitker, G., Florenzano, A., Torri, P., Ruiz, M., & Mercuri, A. M. (2023). Human-Impact Gradients through Anthropogenic Pollen Indicators in a Mediterranean Mosaic Landscape (Balearic Islands). Sustainability, 15(11), 8807. https://doi.org/10.3390/su15118807

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop