# Vulnerability of Mediterranean Ecosystems to Long-Term Changes along the Coast of Israel

David Kaniewski, Elise van Campo, Christophe Morhange, Joel Guiot, Dov

Zviely, Sabine Le Burel, Thierry Otto, Michal Artzy

► To cite this version:

David Kaniewski, Elise van Campo, Christophe Morhange, Joel Guiot, Dov Zviely, et al.. Vulnerability of Mediterranean Ecosystems to Long-Term Changes along the Coast of Israel. PLoS ONE, 2014, 9 (7), 10.1371/journal.pone.0102090. [hal-01457741](https://hal.science/hal-01457741v1)

HAL Id: hal-01457741

https://hal.science/hal-01457741v1

Submitted on 5 Oct 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
# Vulnerability of Mediterranean Ecosystems to Long-Term Changes along the Coast of Israel

David Kaniewski¹,²,³*, Elise Van Campo¹,², Christophe Morhange³,⁴, Joël Guiot⁴, Dov Zviely⁵, Sabine Le Burel¹,², Thierry Otto¹,², Michal Artzy⁵,⁶

¹Laboratoire d'Ecologie Fonctionnelle et Environnement (EcoLab), Université Paul Sabatier-Toulouse 3, Toulouse, France, ²Laboratoire d'Ecologie Fonctionnelle et Environnement (EcoLab), Centre National de la Recherche Scientifique, Toulouse, France, ³Institut Universitaire de France, Paris, France, ⁴Centre de Recherche et d'Enseignement de Géosciences de l'Environnement (CEREGE), Unité Mixte 34, Aix-Marseille Université, Centre National de la Recherche Scientifique, Aix-en-Provence, France, ⁵Leon Recanati Institute for Maritime Studies, University of Haifa, Haifa, Israel, ⁶Hatter Laboratory, University of Haifa, Haifa, Israel

<figure><img src="image_2.png" /></figure>

## Abstract

Although human activity is considered to be a major driving force affecting the distribution and dynamics of Mediterranean ecosystems, the full consequences of projected climate variability and relative sea-level changes on fragile coastal ecosystems for the next century are still unknown. It is unclear how these waterfront ecosystems can be sustained, as well as the services they provide, when relative sea-level rise and global warming are expected to exert even greater pressures in the near future (drought, habitat degradation and accelerated shoreline retreat). Haifa Bay, northern Israel, has recorded a landward sea invasion, with a maximum sea penetration 4,000 years ago, during an important period of urban development and climate instability. Here, we examine the cumulative pressure of climate shifts and relative sea-level changes in order to investigate the patterns and mechanisms behind forest replacement by an open-steppe. We provide a first comprehensive and integrative study for the southern Levant that shows that (i) human impact, through urbanization, has been the main driver behind ecological erosion in the past 4,000 years; (ii) climate pressures have reinforced this impact; and (iii) local coastal changes have played a decisive role in eroding ecosystem resilience. These three parameters, which have closely interacted during the last 4,000 years in Haifa Bay, clearly indicate that for an efficient management of the coastal habitats, anthropogenic pressures linked to urban development must be reduced in order to mitigate the predicted effects of Global Change.

**Citation:** Kaniewski D, Van Campo E, Morhange C, Guiot J, Zviely D, et al. (2014) Vulnerability of Mediterranean Ecosystems to Long-Term Changes along the Coast of Israel. PLoS ONE 9(7): e102090. doi:10.1371/journal.pone.0102090

**Editor:** Grant Ballard, Point Blue Conservation Science, United States of America

**Received:** March 18, 2014; **Accepted:** June 13, 2014; **Published:** July 8, 2014

**Copyright:** © 2014 Kaniewski et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability:** The authors confirm that all data underlying the findings are fully available without restriction. All data are included within the Supporting Information files.

**Funding:** Support was provided by the Institut Universitaire de France, CLIMSORIENT and Geoarchaeology of ancient harbours projects, the Mistrals-ENVIMED-GEOISRAEL program, and by the Hatter Laboratory, Recanati Institute for Maritime Studies, University of Haifa. This work is also a contribution to the Labex OT-Med (ANR-11-LABEX-0061) funded by the « Investissements d'Avenir », French Government project of the French National Research Agency (ANR) through the A\*Midex project (ANR-11-IDEX-0001-02). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

* Email: david.kaniewski@univ-tlse3.fr

## Introduction

Mitigating the impact of ongoing climate change and increasing human activities on ecosystems is one of the biggest challenges of the coming decades, at all geographical scales and across all economic spheres [1–5]. Coastal ecosystems are among the most threatened because relative sea-level rise associated with global warming renders coastal habitats particularly vulnerable [6–7]. High-resolution climate simulations conducted worldwide predict that temperature increases will severely affect the Eastern Mediterranean (EM) over the next century [8–10]. An intensification of extreme weather events is also predicted to interact with patterns of species distribution even if non-climatic influences may dominate local short-term biological shifts [11]. Moreover, regional patterns of warming-induced changes in hydrology seem to be more complex and less certain than those in temperature, chiefly in the EM's seabord areas [8,12–13]. Shoreline changes, which may result from a range of interacting physical and

chemical processes, will also have important consequences for biological features of the EM coastlines [14–15].

Despite the importance of Mediterranean coastal vegetation in earth-atmosphere energy budgets [16] and ecosystem functioning [17], there is still a paucity of empirical data concerning long-term ecosystem distributions in the face of human impacts, climate pressures and coastline changes. This gap is most dramatic for the southeastern Mediterranean coasts, an outstanding center of biodiversity but also one of the most threatened, by economic activity, coastal modifications and drought [6]. Here we investigate 6000 years of coastal ecosystem dynamics in Haifa Bay (Fig. 1), the northernmost depositional sediment sink of the Nile littoral cell [18]. Pollen-derived information on vegetation changes from Akko (Acre), Akko Plain, was used to reconstruct biological trends and assess the role of the local landward sea invasion and subsequent coastal progradation in ecological erosion at the millennial/centennial scale. Previous research conducted at Akko has highlighted human impacts since the early phases of urban
development [19]. Here, the data are used to decipher the relative pressures resulting from the sea ingression and climate shifts. The long-term changes along the coasts of Israel [7,20] that have caused strong biological stresses and major ecological alterations during the past 6000 years, can be seen as a model to understand the consequences of Mediterranean relative sea-level rise for the next century (up to 61 cm) [14].

## Haifa Bay

The coastal part of Israel is subject to the basic pattern of Mediterranean climate with mild rainy winter and hot summers, with decreasing amount of rainfall and increasing temperature southward. At Akko, annual precipitation and mean temperature are 499 mm and 20°C, respectively. The Mediterranean forests of Israel historically comprised evergreen, sclerophyllous and pine tree species, which were destroyed over several millennia by ancient civilizations for agriculture, shipbuilding and fuel. Today the predominant species in the low altitude coastal woodlands are the carob (*Ceratonia siliqua*), Pistachio (*Pistacia lentiscus*), oak (*Quercus calliprinos*) and Aleppo pine (*Pinus halepensis*) [21]. In the Zevulun Plain, a rift valley between Haifa and Akko, underground flow of seawater inland and its mixing with freshwater enables the development of salt marshes populated by *Arthrocnemum macrostachyum*, *Sarcocornia fruticosa*, *Limonium narbonense*, *Atriplex portulacoides* and *Tamarix tetragyna* [21].

## Methods

### Core and modern samples

Biological indicators, extracted from a 215-cm continuous core (32°54′N, 35°05′E; +3 m above sea-level, about 600 m from the sea) drilled in the northern part of Haifa Bay, were used to reconstruct the long-term ecosystem dynamics and plant distribution in the coastal area. Work permits were obtained through to the University of Haifa (the Leon Recanati Institute for Maritime Studies and the Hatter Laboratory). The field studies did not involve endangered or protected species. The chronology is based on accelerator mass spectrometry (AMS) radiocarbon (¹⁴C) ages of short-lived terrestrial samples (seeds; Fig. 2). Modern plant species distributions (Fig. 3) were assessed using pollen deposition. Sample locales (Table S1) were chosen to cover a variety of vegetation types from the sandy shore (32°58′N, 35°04′E; +1 m above sea-level, about 10 m from the sea) to the coastal hinterland (32°54′N, 35°06′E; +8 m above sea-level, about 2.4 km from the sea; Table 1). They were also selected within the broad city limits, in order to ensure that vegetation communities were subjected to similar anthropogenic pressure and that the “distance to the sea” was their main source of ecological stress (inversely related to the distance). The distances to the sea for the modern samples are indicated in Figure 3 and in Table 1.

### Biological data

A total of 86 samples were prepared for pollen analysis using the standard palynological procedure for clay and surface samples

<figure><img src="image_3.png" /></figure>

**Figure 1. Map of Haifa Bay in northern Israel.** The study site is denoted by a black star. The maximum marine ingression during the Bronze Age is shown by a dotted line, with the present coastline as reference.
doi:10.1371/journal.pone.0102090.g001
<table><caption>**Table 1.** Geographical details of the modern samples.</caption><thead><tr><th rowspan="2">Modern samples</th><th colspan="2">Coordinates</th><th>Altitude</th><th>Distance to the sea</th></tr><tr><th>North</th><th>East</th><th>(m a.s.l.)</th><th>(m)</th></tr></thead><tbody><tr><td>1</td><td>32°55′50.51″</td><td>35°04′12.99″</td><td>5</td><td>13</td></tr><tr><td>2</td><td>32°55′18.15″</td><td>35°04′20.88″</td><td>7</td><td>17</td></tr><tr><td>3</td><td>32°55′15.92″</td><td>35°04′31.05″</td><td>5</td><td>31</td></tr><tr><td>4</td><td>32°56′11.52″</td><td>35°04′22.43″</td><td>4</td><td>40</td></tr><tr><td>5</td><td>32°54′39.70″</td><td>35°04′54.80″</td><td>2</td><td>59</td></tr><tr><td>6</td><td>32°54′43.08″</td><td>35°04′59.63″</td><td>3</td><td>210</td></tr><tr><td>7</td><td>32°54′30.28″</td><td>35°05′14.73″</td><td>3</td><td>519</td></tr><tr><td>8</td><td>32°54′41.44″</td><td>35°05′23.18″</td><td>2</td><td>775</td></tr><tr><td>9</td><td>32°53′33.73″</td><td>35°05′14.02″</td><td>6</td><td>813</td></tr><tr><td>10</td><td>32°52′45.17″</td><td>35°05′12.96″</td><td>8</td><td>1138</td></tr><tr><td>11</td><td>32°51′50.38″</td><td>35°04′59.01″</td><td>10</td><td>1294</td></tr><tr><td>12</td><td>32°52′43.44″</td><td>35°05′26.49″</td><td>7</td><td>1493</td></tr><tr><td>13</td><td>32°57′20.58″</td><td>35°05′30.38″</td><td>12</td><td>1591</td></tr><tr><td>14</td><td>33°03′52.51″</td><td>35°07′25.82″</td><td>16</td><td>1828</td></tr><tr><td>15</td><td>33°00′20.61″</td><td>35°06′45.12″</td><td>20</td><td>2187</td></tr><tr><td>16</td><td>32°59′38.26″</td><td>35°06′42.90″</td><td>27</td><td>2512</td></tr></tbody></table>

The coordinates (°N, °E), the altitudes and the distances to the sea are given.
doi:10.1371/journal.pone.0102090.t001

[22]. Pollen grains and dinoflagellate cysts were counted under x400 and x1000 magnification. Pollen frequencies (%) are based on the terrestrial pollen sum excluding local hygrophytes and the spores of non-vascular cryptogams. The arboreal pollen was summed and is depicted as a single curve (using the “subtract mean” function) to provide visualisation of the deviation from the

average value (Fig. 2). Dinoflagellate cysts were counted on pollen-slides and are displayed as concentrations (cysts per cm⁻³). The dinoflagellate cysts (resting stage of marine plankton), were used as an independent proxy for the sea incursion [23]. The biological data were resampled (AnalySeries 2.0; sampling grid: one sample

<figure><img src="image_4.png" /><figcaption>**Figure 2. AMS <sup>14</sup>C calibrated ages (1σ-68%) and suggested age-depth curve.** The lithology of the core sampled at Akko is detailed and plotted on a linear depth-scale. The arboreal pollen curve from the core is shown as “subtract mean” function and plotted on a linear depth-scale.
doi:10.1371/journal.pone.0102090.g002</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Modern pollen deposition from the northern Haifa Bay expressed as pollen-derived clusters. The polynomial curves indicate the general trend. A variety of vegetation types from the sandy shore (open steppe) to the coastal hinterland (woodland) were selected. doi:10.1371/journal.pone.0102090.g003</figcaption></figure>

every 50 years) in order to obtain an identical chronological scale for all the proxies.

## Multivariate analyses

Our aim is to integrate each sample from the core into a modern vegetation transect, based on its "distance to the sea". Our goal is to study the vegetation dynamics during the last 6000 years, using present distributions as analogues. All numerical analyses were performed with the software PAST.

All pollen data (fossil and modern) were transformed into pollen-derived clusters ($^{Pd}$Clusters) using Neighbour Joining (NJ) analysis (Fig. 4). The NJ method is an alternative process for hierarchical cluster analysis, finding hierarchical groupings in

multivariate data sets. Here, it is based on pollen-type time-series (presence/absence and abundance). NJ analysis was used to compute the lengths of branches of a tree, using branches as ecological distances between groups of taxa (descending type). NJ was computed using correlation as similarity measure and final branch as root. The pollen-types from each cluster were summed (Tables S1 and S2) to create local pollen-derived vegetation patterns (Figs 3 and 4).

Each sample from the core, expressed as PdCluster-scores, was integrated in the modern pollen deposition matrix (also expressed as PdCluster-scores), and then assigned to a modern distribution using Principal Components Analysis (PCA). The method is based on the PCA axes 1 and 2 (that load the maximum amount of total variance) and only selects the closest score between one modern sample and the fossil sample. The fossil sample was, according to the outcome of the “closest score” method, affiliated to one modern sample. Each affiliation was controlled by cluster analysis (with paired group as algorithm and correlation as the similarity measure) and by k-means clustering (number of clusters: 16), using the same pollen deposition matrix as the PCA and testing the fossil samples one by one. The palaeo-distributions were then recomputed in one unique curve in order to create a proxy for vegetation dynamics ($^{P}$VD) along the Mediterranean coast of Israel for the last 6000 years (Fig. 5).

The $^{Pd}$Clusters were also analysed using CABFAC analysis (Fig. 6) to test the impact of environmental factors on the $^{P}$VD. The CABFAC factor analysis (Q-mode factor analysis) implements the classical method of factor analysis and environmental regression (CABFAC and REGRESS). Selected environmental data are regressed on the CABFAC factors using the second-order (parabolic) method, with cross terms. The environmental regression model (RMA) reports the observed values against the values reconstructed from the factors. The selected environmental data for this study [anthropogenic pressures; EM precipitation; EM temperature; sea incursion; ratio of actual and potential evaporation; transpiration (E/PE) for the EM; difference between annual precipitation and actual evapotranspiration (P-E) for the EM] cover a wide range of pressures. The consistency of each reconstructed model is indicated by the $R^2$. The $R^2$ values were ordered from the highest to the lowest, and used as a hierarchy of pressures on ecosystems.

The link between these independent proxies and the $^{Pd}$Clusters was established using NJ analysis (Fig. 7). This clustering analysis shows the highest affinities between environmental pressures and the $^{Pd}$Clusters, suggesting a close relationship during the last 6000 years.

## Smoothing

In figures 3 and 5, polynomial smoothing curves were selected to synthesise the main changes observed at the site and in the EM. The best smoothing was chosen according to the highest $R^2$.

## Vegetation model

The vegetation model BIOME4 [24] is used in an inverse mode to reconstruct climate along the pollen diagram taking into account variations in the atmospheric concentration of CO$_2$ [25]. The reconstruction is based on a Bayesian approach where a priori information on the climate at a given site and time period is refined given the corresponding pollen assemblages and the BIOME4 model, assuming that the differences between model simulations and data follow a Gaussian law. The results are validated using regional modern pollen assemblages. The reconstructions for the Holocene are expressed as anomalies from the estimated modern climate (top core). In this paper, we focus on the
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Neighbour Joining analysis of the main taxa (modern and fossil) computed with *correlation* as similarity measure and *final branch* as root. The pollen-types from each cluster were summed to create pollen-derived vegetation patterns.
doi:10.1371/journal.pone.0102090.g004</figcaption></figure>

ratio of actual and potential evapotranspiration (E/PE), which is an index of drought felt by vegetation, and on the difference between annual precipitation and actual evapotranspiration (P-E), which is an estimate of the site's contribution to run-off in the catchment basin (Fig. 5 and Table S3). The two other proxies (temperature and precipitations) for the Eastern Mediterranean are derived from Kaniewski *et al.* [19] (Table S4). All reconstructions concern the Eastern Mediterranean (data from Turkey, Cyprus, Syria, and Israel) and are considered as independent proxies.

## Results

The PVD shows that, during the sea incursion deduced from the presence of Dinoflagellate cysts [4000-2900 calibrated year Before the Present (cal yr BP)], the coastal vegetation evolved from a wet

woodland through a shrubland to an open-steppe (Fig. 5) similar to those that can be found nowadays along the coast of Akko, western Galilee (Fig. 3). Before the sea incursion, from c. 6000 to 4000 cal yr BP, the area was dominated by wet woodland (Fig. 4), with a limited influence of the salt spray brought inland by the northwesterly sea breeze. The relative absence of dinoflagellate cysts suggests no sea incursion and no direct pressures on ecosystems by saline waters. The P-E polynomial curve shows that run-offs were important, yielding higher local freshwater budget/outflows in the EM. Sea-level rise, that started around 8000 cal yr BP in the Haifa Bay and reached the coring site at 4000 cal yr BP, acted like an erosive agent that cumulated with human pressures and climate instability (Fig. 5), leading to an alteration of ecosystems (Fig. 2) and to a shift from a wet woodland to a thorny shrub-steppe (Fig. 5). The sea incursion is directly
<figure><img src="image_7.png" /><figcaption>**Figure 5.** A 6000-year environmental reconstruction of Haifa Bay. P-E, E/PE, precipitation, temperatures and anthropogenic activities are displayed as scores, mm. yr⁻¹, °C and %. A polynomial curve was established for each component and plotted on a linear timescale. The coastal ecosystem dynamics and the dinoflagellate cysts (marine phytoplankton) were also plotted on a linear timescale. The numbers associated with the coastal ecosystem dynamics are related to the modern samples (Fig. 3). The marine ingression is underlined in blue and is bracketed between horizontal dotted lines.</figcaption></figure>

doi:10.1371/journal.pone.0102090.g005

associated with the development of the PdClusters dry steppe and phrygana in the coastal area (Fig. 7). The P-E values indicate a decrease of freshwater budget/outflows, mainly during three periods: 4200-3800, 3500 and 3200-2900 cal yr BP. The E/PE drops after 4000 cal yr BP, reaching its lowest value at 2950 cal yr BP (0.43). Coastal progradation started around 3400-3300 cal yr BP and ended at 2900 cal yr BP at the coring site (disappearance of dinoflagellate cysts) and is concomitant with one of the driest/warmest phases from 3200 to 2900 cal yr BP. Owing to the cumulative pressures, the shrub-steppe turned into an open-steppe and human activity decreased (Fig. 5). The dry wooded components then recovered and dominated until 2100 cal yr BP, mainly led by anthropogenic pressures (Fig. 7). Both steppe and shrubs have remained dominant since 2100 cal yr BP (Fig. 2), and have increased in sub-recent times (Fig. 5), despite more important run-off and a higher local freshwater budget.

## Discussion

A full understanding of the ecological consequences of future relative sea-level variations not only needs to integrate numerous amplifiers (e.g. human impact, urban development and increasing drought), but also to consider the complex mechanisms at work at various timescales [26]. At Akko, the RMAs show that the anthropogenic pressures, with an R² value of 0.81 (Fig. 7), best explain the vegetation dynamics. The high values associated with climate pressures [precipitation (R² = 0.75) and temperature (R² = 0.71)] and sea level changes (R² = 0.70), suggest that these three parameters have interacted over the last 6000 years.

## Sea ingression

The altimetry data of the Topex/Poseidon satellite have shown that sea level rose continuously at a rate of up to 20 mm per yr⁻¹ over the period 1993–1999 in the Levantine basin [27]. The predicted maximum value for the year 2100 in the EM is 60 cm higher than the 1990 value, providing an analogue with Haifa Bay where the sea transgressed eastwards (>2.5 km) and flooded the

Zevulum Plain (Fig. 1) [28–29]. Regardless of the physical mechanisms responsible for the sea ingression on the Zevulum Plain and the local influence of the hydro-sedimentary budget from the Na’aman River [28], the added pressure of the sea invasion between 4000 and 2900 cal yr BP on ecosystem dynamics already affected by human activity and climate stress, led to a deep ecological erosion (Fig. 5).

The main consequences of the sea invasion were a retreat of the coastal forest (drop in forested/wooded PdClusters) (Fig. 5), a loss of resilience and a disappearance of the initial local biogeographic zonation. The forest replacement by a thorny shrub-steppe and then by an open-steppe appeared to follow, rather than cause, failure of tree regeneration. To sum up, several causes cumulatively mediated forest dynamics between 4000 and 2900 cal yr BP: (i) the lack of regenerating stands of different tree species responsible for the retreat of the coastal forest [26,30–31] which probably resulted from the increases in saltwater intrusions, and coastal erosion that washed away entire tree stands and provoked seedling burial by shifting sands along the sandy coast of Haifa Bay [18]; (ii) changes of the physical characteristics of the pedologic substrate following coastal progradation, with increasing salt concentration and sand deposition in the Bay; and (iii) increased human pressures and climate forcing [19].

## Exposure to salt

Due to the intrusion of the saline water table in fresh water streams, the freshwater wetland associated with the Na'aman River was deeply impacted after 4000 cal yr BP, with a fall in hygrophilous-hydrophilous herbs (Fig. 8). The subsequent coastal progradation that started at 3400-3300 cal yr BP (Fig. 5) left an eroded sandy-salty area, colonized by a steppic vegetation that became dominant until the end of the shoreline retreat (2900 cal yr BP). A similar process was observed in the Salt Lake of Larnaca, Cyprus, where a shift from sheltered marine to lagoonal environments between 3400 and 3300 cal yr BP produced an ecological change with a strong increase in xerophytic vegetation-
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Environmental regression models from Haifa Bay. The selected environmental data are indicated below each graph. The consistency of each reconstructed model is indicated by the R². The lower graph corresponds to a ranking of the environmental data according to the R² values.
doi:10.1371/journal.pone.0102090.g006</figcaption></figure>

types colonizing the shores that were no longer washed by seawater [32].

## The end of the sea incursion

During the marine stage at Akko, the main impact on ecosystems was identified at 3200 cal yr BP, when the physical processes linked with the regression of the sea were added to higher climate pressures, leading to decreased human activities (Fig. 5). The 3.2 kyr BP drought event is now well identified in different locations from Syria [33–35], Cyprus [32], Israel [36–41], and Egypt [42]. The signature of this climate event is particularly strong at Akko because the effects of drought were reinforced by the physical processes linked to the coastal changes,

<figure><img src="image_9.png" /><figcaption>**Figure 7.** Neighbour Joining analysis of the selected pressures and pollen-derived vegetation patterns. The dendrogram shows the statistical link between one pressure and pollen-derived vegetation patterns, suggesting an environmental relationship between the pressure and the resulting vegetation. Each numerical link between a pressure and a pollen-derived vegetation pattern is denoted by a different colour.
doi:10.1371/journal.pone.0102090.g007</figcaption></figure>

causing a dramatic demise in wooded ecosystems (Fig. 5). After 4000 cal yr BP, the wet woodland, the open-forest, and the fen trees never regained their former importance in the coastal area. The wooded components recovered after 2850 cal yr BP, but only dry shrubland (Fig. 4) was able to develop, until 2100 cal yr BP (Fig. 5).

## Coastal retreat

Shoreline migration analyses of the coast of Akko, northern Haifa Bay, during the last 200 years, show that the coast retreated by 15 meters from 1926 to 2006, primarily due to human-induced modification of the Na'aman River's hydro-sedimentary budget, changes in sediment transport, and coastal structures built during the last 50 years [43]. This new phase of sea-level invasion is well attested by high values of steppic vegetation (Fig. 5) and a renewed drop in hygrophilous-hydrophilous herbs (Fig. 8) at the very end of the sequence. This suggests that similar stresses generate analogous biological processes, whatever the period under consideration.

## Clues for the future

The long-term environmental evolution of the northern part of Haifa Bay suggests that the future of coastal ecosystems and wetlands in the southern Levant will mainly be mediated by
<figure><img src="image_10.png" /><figcaption>**Figure 8.** Hygrophilous and hydrophilous herbs from Haifa Bay. The hygrophilous and hydrophilous herbs, expressed as percentages, are displayed as a single curve plotted on a linear timescale. The average value is indicated by a vertical dotted line. The maximum sea ingression is bracketed between the two lines.
doi:10.1371/journal.pone.0102090.g008</figcaption></figure>

complex interactions between human impacts, climate change and
changing coastal geomorphology. This questions human behavior
in these fragile coastal systems: how will humans respond to rising
relative sea level and increased drought, and how will they
produce a sustainable management of coastal ecosystems and
wetlands while their encroachment already affects the landscape,
through an overexploitation of resources and by an overuse of
freshwater. In Haifa Bay, human activities have played an
important role during the last 4000 years, by reinforcing and
gradually supplanting the natural stressors. Efficient management
of the coastal habitats of southern Levant would imply a strong
reduction in anthropogenic pressures linked to urban develop-
ment, in order to mitigate the predicted effects of Global Change.

References

1. Sala OE, Chapin III FS, Armesto JJ, Berlow E, Bloomfield J, et al. (2000) Global biodiversity scenarios for the year 2100. Science 287: 1770-1774.

2. Rahmstorf S (2007) A semi-empirical approach to projecting future sea-level rise. Science 315: 368-370.

3. Chenoweth J, Hadjinicolaou P, Bruggeman A, Lelieveld J, Levin Z, et al. (2011) Impact of climate change on the water resources of the eastern Mediterranean and Middle East region: Modeled 21st century changes and implications. Water Resources Research 47: W06506; DOI:10.1029/wr010269.

4. Soltani S, Saboohi R, Yaghmaei L (2012) Rainfall and rainy days trend in Iran. Climatic Change 110: 187-213.

5. Milano M, Ruelland D, Dezetter A (2013) Anticipating water stress in the Mediterranean: a decision support approach for better water resources management. Water Resources and natural environment Plan Bleu Notes 26: 1-4.

6. Nicholls RJ, Wong PP, Burkett VR, Codignotto JO, Hay JE, et al. (2007) Coastal systems and low-lying areas. Climate Change 2007: Impacts, Adaptation and Vulnerability. In: Parry ML, Canziani OF, Palutikof JP, van der Linden PJ, Hanson CE, editors. Contribution of working group II to the fourth assessment

Modern cities in the Levant are undergoing strong demographic
growth, but their ongoing expansion needs to carefully assess the
environmental limits of local natural resources and ecosystems'
capacity to respond.

Conclusions

Haifa Bay's biogeographical evolution at millennial- to centen-
nial time-scales suggests that projected sea-level rise scenarios for
the next century will deeply impact the coastal ecosystems,
especially within the context of more intense drought stress,
increased human activities and diminished sediment supply by
local rivers. The major physical processes that have shaped the
ecosystems since 6000 cal yr BP can be attributed to climate
fluctuations (Fig. 6) and to the sea-level transgression 4000 cal yr
BP and subsequent coastal retreat since 3400-3300 cal yr BP.
Stand-level regeneration failures of the forest-woodland vegeta-
tion, as well as a high rate of freshwater wetland loss, might be
expected in the coastal southern Levant in the near future.

**Supporting Information**

<caption>**Table S1 Data - Modern pollen samples.**
(XLSX)</caption>

<caption>**Table S2 Data - Samples from the core.**
(XLSX)</caption>

<caption>**Table S3 Data - Reconstructed climate values for the Eastern Mediterranean, E/PE and P-E.**
(XLSX)</caption>

<caption>**Table S4 Data - Reconstructed climate values for the Eastern Mediterranean, annual precipitation, mean temperature of the coldest month, mean annual temperature, mean temperature of the warmest month.**
(XLSX)</caption>

Acknowledgments

We wish to thank Patricia L. Fall, José S. Carrión and an anonymous
reviewer for their critical remarks and useful recommendations.

Author Contributions

Conceived and designed the experiments: DK EVC. Performed the experiments: DK EVC CM JG SLB TO. Analyzed the data: DK EVC CM JG DZ SLB TO MA. Contributed reagents/materials/analysis tools: DK EVC CM JG DZ SLB TO MA. Contributed to the writing of the manuscript: DK EVC JG. Obtained permission for all the drillings: MA.

report of the intergovernmental panel on climate change. Cambridge:
Cambridge University Press. 315-356.

7. Anthony EJ (2009) Shore processes and their palaeoenvironmental applications. Developments in marine geology vol. 4. Amsterdam: Elsevier Science, 519 p.

8. Milly PCD, Dunne KA, Vecchia AV (2005) Global pattern of trends in streamflow and water availability in a changing climate. Nature 438: 347-350.

9. Suppan P, Kunstmann H, Heckel A, Rimmer A (2008) Impact of climate change on water availability in the Near East. In: Zereini F, Hotzl H, editors. Climate changes and water resources in the Middle East and North Africa. Berlin: Springer. 47-58.

10. Evans JP (2009) 21st century climate change in the Middle East. Climatic Change 92: 417-432.

11. Parmesan C, Yohe G (2003) A globally coherent fingerprint of climate change impacts across natural systems. Nature 424: 37-42.

12. Kafle HK, Bruins HJ (2009) Climatic trends in Israel 1970-2002: warmer and increasing aridity inland. Climatic Change 96: 63-77.

13. Sowers J, Vengosh A, Weinthal E (2011) Climate change, water resources, and the politics of adaptation in the Middle East and North Africa. Climatic Change 104: 599-627.
14. Marcos M, Tsimplis MN (2008) Comparison of results of AOGCMs in the Mediterranean Sea during the 21st century. Journal of Geophysical Research 113: C12028; DOI:10.1029/jc004820.

15. Lichter M, Zviely D, Klein M, Sivan D (2010) Sea level changes in the Mediterranean: past, present and future - a review. In: Israel A, Einav R, Seckbach J, editors. Seaweeds and their role in globally changing environments. Berlin: Springer. 3-17.

16. Sellers PJ, Tucker CJ, Collatz GJ, Los SO, Justice CO, et al. (1996) A revised land surface parameterization (SiB2) for atmospheric GCMs - part II: the generation of global fields of terrestrial biophysical parameters from satellite data. Journal of Climate 9: 706-737.

17. Sternberg M, Shoshany M (2001) Influence of slope aspect on Mediterranean woody formations: Comparison of a semiarid and an arid site in Israel. Ecological Research 16: 335-345.

18. Zviely D, Kit E, Klein M (2007) Longshore sand transport estimates along the Mediterranean coast of Israel in the Holocene. Marine Geology 237: 61-73.

19. Kaniewski D, Van Campo E, Morhange C, Guiot J, Zviely D, et al. (2013a) Early urban impact on Mediterranean coastal environments. Scientific Reports 3: 3540; DOI:10.1038/srep03540.

20. Stewart IS, Morhange C (2009) Coastal geomorphology and sea-level change. In: Woodward J, editor. The Physical Geography of the Mediterranean. Oxford: Oxford University Press. 385-414.

21. Danin A (1981) The impact of geomorphologic and climatic conditions on the vegetation of salt marshes along the Mediterranean coast of Israel and Sinai. Actas III Congr. Optima. Anales Jard. Bot. Madrid 37: 269-275.

22. Faegri K, Iversen I (1989) Textbook of Pollen Analysis, fourth edition. New York: J. Wiley and Sons, 340 p.

23. Gräfe KU, Wendler JE (2003) Foraminifers and calcareous dinoflagellate cysts as proxies for deciphering sequence stratigraphy, sea-level change, and paleocea- nography of Cenomanian-Turonian hemipelagic sediments in Western Europe. Society for Sedimentary Geology 75: 229-262.

24. Kaplan JO, Bigelow NH, Prentice IC, Harrison SP, Bartlein PJ, Christensen TR, Cramer W, Matveyeva NV, McGuire AD, Murray DF (2003) Climate change and arctic ecosystems: 2. Modeling, paleodata-model comparisons, and future projections. Journal of Geophysical Research 108: 8171.

25. Guiot J, Wu H, Garreta V, Hatté C, Magny M (2009) A few prospective ideas on climate reconstruction: from a statistical single proxy approach towards a multi- proxy and dynamical approach. Climate of the Past 5: 571-583.

26. Williams K, Ewel KC, Stumpf RP, Putz FE, Workman TW (1999) Sea level rise and coastal forest retreat on the west coast of Florida, USA. Ecology 80: 2045- 2063.

27. Cazenave A, Cabanes C, Dominh K, Mangiarotti S (2001) Recent sea level change in the Mediterranean Sea revealed by Topex-Poseidon satellite altimetry. Geophysical Research Letters 28: 1607-1610.

28. Zviely D, Sivan D, Ecker A, Bakler N, Rohlrich V, et al. (2006) Holocene evolution of the Haifa Bay area, Israel, and its influence on ancient tell settlements. The Holocene 16: 849-861.

29. Porat N, Sivan D, Zviely D (2008) Late Holocene embayment infill and shoreline migration, Haifa Bay, Eastern Mediterranean. Israel Journal of Earth Sciences 57: 21-31.

30. Williams K, MacDonald M, da Silveira Lobo Sternberg L (2003) Interactions of storm, drought, and sea-level rise on coastal forest: a case study. Journal of Coastal Research 19: 1116-1121.

31. Desantis LRG, Bhotika S, Williams K, Putz FE (2007) Sea-level rise and drought interactions accelerate forest decline on the Gulf Coast of Florida, USA. Global Change Biology 13: 2349-2360.

32. Kaniewski D, Van Campo E, Guiot J, LeBurel S, Otto T, et al. (2013b) Environmental roots of the Late Bronze Age crisis. Plos One 8: e71004; DOI:10.1371/journal.pone.0071004.

33. Fiorentino G, Caracuta V, Calcagnile L, D'Elia M, Matthiae P, et al. (2008) Third millennium B.C. climate change in Syria highlighted by Carbon stable isotope analysis of $^{14}$C-AMS dated plant remains from Ebla. Palaeogeography, Palaeoclimatology, Palaeoecology 266: 51-58.

34. Kaniewski D, Paulissen E, Van Campo E, Al-Maqdissi M, Bretschneider J, et al. (2008) Middle East coastal ecosystem response to middle-to-late Holocene abrupt climate changes. Proceedings of the National Academy of Sciences of the USA 105: 13941-13946.

35. Kaniewski D, Paulissen E, Van Campo E, Weiss H, Otto T, et al. (2010) Late Second-Early First Millennium BC abrupt climate changes in coastal Syria and their possible significance for the history of the Eastern Mediterranean. Quaternary Research 74: 207-215.

36. Schilman B, Bar-Matthews M, Almogi-Labin A, Luz B (2001) Global climate instability reflected by Eastern Mediterranean marine records during the Late Holocene. Palaeogeography, Palaeoclimatology, Palaeoecology 176: 157-176.

37. Schilman B, Ayalon A, Bar-Matthews M, Kagan EJ, Almogi-Labin A (2002) Sea-Land paleoclimate correlation in the Eastern Mediterranean region during the Late Holocene. Israel Journal of Earth Sciences 51: 181-190.

38. Bar-Matthews M, Ayalon A, Gilmour M, Matthews A, Hawkesworth CJ (2003) Sea-land oxygen isotopic relationships from planktonic foraminifera and speleothems in the Eastern Mediterranean region and their implication for paleorainfall during interglacial intervals. Geochimica et Cosmochimica Acta 67: 3181-3199.

39. Litt T, Ohlwein C, Neumann FH, Hense A, Stein M (2012) Holocene climate variability in the Levant from the Dead Sea pollen record. Quaternary Science Reviews 49: 95-105.

40. Langgut D, Finkelstein I, Litt T (2013) Climate and the Late Bronze Collapse: new evidence from the Southern Levant. Tel Aviv 40: 149-175.

41. Langgut D, Neumann FH, Stein M, Wagner A., Kagan EJ, et al. (2014) Dead Sea pollen record and history of human activity in the Judean Highlands (Israel) from the Intermediate Bronze into the Iron Ages (~2500-500 BCE). Palynology, DOI:10.1080/01916122.2014.906001.

42. Bernhardt CE, Horton BP, Stanley JD (2012) Nile Delta vegetation response to Holocene climate variability. Geology 40: 615-618.

43. Zviely D, Kit E, Rosen B, Galili E, Klein M (2009) Shoreline migration and beach-nearshore sand balance over the last 200 years in Haifa Bay (SE Mediterranean). Geo-Marine Letters 29: 93-110.