Abstract
The first large-scale network of 79 tree-ring chronologies in the Eastern Mediterranean and Near East (EMNE; 33°N–42°N, 21°E–43°E) is described and analyzed to identify the seasonal climatic signal in indices of annual ring width. Correlation analysis and cluster analysis are applied to tree-ring data and gridded climate data to assess the climate signal embedded in the network in preparation for climate field reconstructions and formal proxy/model intercomparison experiments. The lengths of the 79 combined chronologies range from 89 to 990 years. The monthly correlations and partial correlations reveal a pervasive positive association with May, June, and sometimes July precipitation, positive correlations with winter and spring (December through April) temperatures, and negative relationships with May through July temperature, although as expected, there are site-to-site exceptions to these general patterns. Cluster analysis suggests three groups of sites based on their association with climate. The chronologies for the EMNE have coherent seasonal precipitation and temperature signals across a fairly broad geographical domain. The predominant signal is a positive growth response to May–June precipitation. Collectively, the findings suggest that the network can be exploited to develop season-specific field reconstructions of precipitation and drought history in the EMNE.
Keywords
Introduction
The Eastern Mediterranean and Near East (EMNE) region is influenced by multiple local and remote climate processes. The region is exposed to the distal influences of the Asian Monsoon (Raicich et al., 2003; Rodwell and Hoskins, 1996) and the proximal effects of the East Atlantic Jet (Duenkeloh and Jacobeit, 2003; Touchan et al., 2005b; Xoplaki et al., 2004) in summer. The Siberian High Pressure System (e.g. Xoplaki et al., 2001) and North Atlantic Oscillation (Corte-Real et al., 1995; Duenkeloh and Jacobeit, 2003; Xoplaki et al., 2004) also influence regional climate. A semi-objective classification of daily synoptic maps identified no fewer than six large synoptic groups important to climatic variation in the eastern Mediterranean: Cyprus lows, the Persian trough, the Red Sea trough, Sharav lows, and Siberian/subtropical highs (Alpert et al., 2004). Precipitation trends in the last half of the 20th century appear to be related to changes in the frequency or annual number of days of these synoptic types, notably the Cyprus lows and Red Sea trough (Alpert et al., 2004). Orography and land-sea interactions (including distance from the sea) and smaller scale processes (Lolis et al., 1999; Xoplaki et al., 2001, 2003a, 2003b, 2004) are also important and can control local-scale pattern climate heterogeneity. Mediterranean climate is further influenced by the Mediterranean Sea itself (e.g. Mariotti et al., 2002; Trigo et al., 1999), which represents an important source of energy and moisture for cyclone development. Furthermore, the complex land topography around the Mediterranean plays a crucial role in steering airflow (e.g. Bartzokas et al., 1994; Trigo et al., 1999). Perhaps, as a consequence of these multiple and superimposed influences, hydrologic variability spans a broad range of timescales and is unlikely to be fully described by the modern instrumental record. Meanwhile, the population of the eastern Mediterranean grows by 3.5% annually, while irrigation practices consume at least 80% of the available water supply. As a result, precipitation is a key variable affecting public health and political stability (De’Donato and Michelozzi, 2014; Rubio, 2009). One example of this is that drought conditions have been associated with devastating fire seasons across the region, resulting in the destruction of hundreds of hectares of forests and crops even in the relatively more humid and colder regions of the Mediterranean Basin (e.g. Dimitrakopoulos et al., 2011; Pausas et al., 2008). Drought also has profound implications for regional food security, societal upheaval, and economic stability (Kaniewski et al., 2012; Sowers et al., 2011). Quantifying and understanding climatic changes at these regional scales are among the most important and uncertain issues in the study of global change. For example, extreme regional-scale droughts exhibit much larger amplitudes than global averages, and affect regional societies, economies, water supplies, and agricultural ecosystems (e.g. Luterbacher et al., 2004; Xoplaki et al., 2005). Consequently, understanding the spatiotemporal details of drought is of critical importance in this region, where even currently, the consequences can be severe, and an increased occurrence of such events is projected.
Refined knowledge and improved understanding of the full range of past hydroclimatic variability in the EMNE is critical for identifying possible causative factors, evaluating the potential for spatially extensive and temporally persistent anomalies, and for assessing the ability of general circulation models (GCMs) to reproduce variability at timescales from decades to centuries. Meteorological data are sparse in the EMNE and are typically not long enough to effectively capture the potential range of multidecadal to century-scale variability and the spatiotemporal response to radiative forcing. In view of existing uncertainties, longer records of natural hydroclimatology are necessary for assessing the causes of variability and trends in the instrumental record and evaluating the accuracy of the forced response in forecast GCMs. Touchan et al. (1999, 2003, 2005a, 2005b, 2007, 2008a, 2008b, 2011, 2012), Akkemik and Aras (2005), Akkemik et al. (2008), Köse et al. (2011), Griggs et al. (2007), Hughes et al. (2001), and D’Arrigo and Cullen (2001) have demonstrated that paleoclimate proxy records developed from trees in this region offer a longer term perspective on episodic drought in the EMNE and have identified the steps needed to improve the development and applicability of such records.
In this paper, a network of tree-ring chronologies in the EMNE (33°N–42°N, 21°E–43°E) is described and is analyzed to identify their seasonal climatic signal. This analysis is a preliminary but necessary step toward application of this network to the study of long-term climate variability, the association between regional climate and atmosphere–ocean circulations, and the ability of GCMs to reproduce important drought-related features of regional climate. We apply correlation and cluster analyses to tree-ring and gridded climate data in order to assess the climate signal embedded in the network as preparation for climate field reconstructions and formal proxy/model intercomparison experiments. The objective of the analysis is the improvement and refinement of reconstruction protocols, including the selection of appropriate target reconstruction and climate model fields, all with the larger goal of understanding Mediterranean climate variability at interannual to centennial scales, providing out-of-sample assessment of GCMs, and evaluating the myriad interacting influences on EMNE drought. These analyses also provide complementary data about the variability of climate controls on tree growth across the region that may be useful for predicting the future spatiotemporal response of key forest species to climate change.
Materials and methods
Tree-ring data and chronology development
This study represents the first large-scale systematic dendroclimatic sampling campaign focused on developing a network of drought-sensitive chronologies from Turkey, Syria, Lebanon, Cyprus, and Greece. Fieldwork conducted over the period 2000–2011 has resulted in development of 79 chronologies from 82 tree-ring sites (Table 1 and Figure 1). Samples were collected from species known to form annual rings and those that demonstrated a high degree of common variation driven by strong climate. Increment cores were taken at all sites, and full cross sections were taken from stumps of cedar and juniper. Samples were surfaced and cross-dated following standard dendrochronological techniques (Stokes and Smiley, 1996). The width of each annual ring on the cores and cross sections was measured to the nearest 0.01 mm.
Site information for the eastern Mediterranean region.
Species key: 1. Pinus nigra, 2. Pinus sylvestris, 3. Juniperus excelsa, 4. Pinus brutia, 5. Cedrus libani, 6. Pinus heldreichii, 7. Abies cephalonica, 8. Cedrus brevifolia, 9. Abies cilicica and 10. Pinus pinea.

Locations of tree-ring chronology sites, with symbol color and size indicating genus and length, respectively.
For the purposes of this study, each series of tree-ring width measurements was fit with a 67% cubic smoothing spline (frequency response 0.50 at a wavelength 2/3 the sample length), to remove nonclimatic trends due to age, size, and the effects of stand dynamics (Cook and Kairiukstis, 1990). The detrended series were then prewhitened with low-order autoregressive models to remove persistence not related to climatic variations. The individual indices were combined into single averaged chronologies for each combination of site and species using a bi-weight robust estimate of the mean (Cook, 1985). Geographic location, confirmed by visual cross-dating and computer-based quality control suggested a strong similarity among several of the resulting 82 residual chronologies. Accordingly, prior to correlation and correlation analysis, the following sites were combined to form single chronologies: (1) PIBR Forest-SYR combines Atera (ATEP), Bait Hamik (BAIP), and Mafrak Bait Ablak (MBAP) and (2) ABCE Forest-SYR combines Bedayat Al Khandak (BKTA) and Rawisat Almedeki (RWIA).
We used expressed population signal (EPS), calculated from data on sample size and between-tree correlation, to assess the adequacy of replication of the chronology (Wigley et al., 1984). The EPS estimates the ability of a finite sample to approximate the true, unknown, population tree-ring signal at the site. Following recommended guidelines (Cook and Kairiukstis, 1990; Wigley et al., 1984), we judged chronologies as adequately sampled for use in climatic reconstruction (e.g. precipitation and drought index) only if the chronology EPS reaches 0.85 before the start of the instrumental climate record.
Climate data and cluster analysis
This analysis was done using 77 of the 79 residual chronologies. We excluded from the analysis two chronologies from Lebanon (HAD, Cedrus libani, and KFR, Pinus pinea; see Table 2) as they did not cover the full span of the instrumental data (beginning only in 1919 and 1923, respectively). The remaining 77 sites were compared against local monthly gridded climate data. We used 1.0° gridded monthly precipitation data from the Global Precipitation Climatology Centre (GPCC; version 6) dataset covering 1901 through 2010 (Becker et al., 2013; Schneider et al., 2011, 2013) and 0.5° gridded temperature data from the Climatic Research Unit (CRU) TS3.1 dataset (Mitchell and Jones, 2005) covering the period from 1901 to 2009. We also examined the correlation between sites in our network and the Palmer Drought Severity Index (PDSI; Palmer, 1965) developed by Van der Schrier et al. (2013), which covers the period from 1901 to 2009.
Summary statistics for the 79 chronologies for the ARSTAN program (Cook and Holmes, 1999; Cook and Krusic, 2005).
MSSL: mean sample segment length; Std: standard deviation; SK: skewness; KU: Kurtosis; EPS: Expressed Population Statistic (Wigley et al., 1984); MCAR: mean correlation among radii; EV: explained variance.
Brutia pine forest combines ATEP, BAIP, and MBAP.
Abies cilicica forest combines BKTA and RWIA.
More samples are needed.
We conducted a site-by-site correlation analysis of each residual ring-width chronology against the local (nearest) gridded climate data using the seasonal correlation (Seascorr) procedure developed by Meko et al. (2011) with exact simulation (Percival and Constantine, 2006) for significance testing. We used individual monthly as well as seasonal values integrating 2, 3, or 4 months. We considered a 14-month window starting in the August prior to the growth year and ending in the following September. We then performed a k-means cluster analysis on the resulting 112 monthly and seasonal correlation and partial coefficients for precipitation and temperature for each site. This allows us to analyze groups of chronologies based on their association with their local climate. We estimated the optimal number of clusters using silhouette plots and the gap statistic (Kaufman and Rousseeuw, 1990; Tibshirani et al., 2001).
Results
Tree-ring chronologies
Data for individual tree-ring chronologies are summarized in Tables 1 and 2. The lengths of the 79 chronologies range from 89 (Hadath, Lebanon) to 990 years (Elmali, Turkey; Table 1). Statistical analyses of each chronology are summarized in Table 2. The mean correlation among individual radii at each site (the interseries correlation) represents the strength of their common signal and ranges from 0.17 to 0.59. The highest interseries correlation is for the NESJ site in Turkey, and the lowest is for the chronology developed from the JRA site in Lebanon. The EPS of three chronologies (including the two already eliminated from consideration due to their limited time span) fails to reach 0.85 by the start of the climatic record, leaving a total of 76 chronologies potentially useful for climatic reconstruction without additional sampling. The mean sample segment length (MSSL) of all 79 chronologies ranges from 68 to 408 years. Half of these chronologies have a MSSL greater than 200 years in length, and several have MSSL exceeding 400 years.
Climate and cluster analysis
Analysis of monthly correlations and partial correlations between gridded climate data and our tree-ring width chronologies reveals a pervasive positive association with May, June, and sometimes July precipitation; positive correlations with winter and spring (December through April) temperatures; and negative relationships with May through July temperature, although as expected, there are site-to-site exceptions to these general patterns (Figure 2).

Monthly correlations with total GPCC v6 1 precipitation (Becker et al., 2013; Schneider et al., 2013) for each eastern Mediterranean site used in the network analysis and partial correlations with mean CRU TS3.1 0.5 temperature (Mitchell and Jones, 2005) for each eastern Mediterranean site used in the network analysis. Sites are correlated against the closest grid cell in each dataset.
Our cluster analysis suggests three groups of sites based on their association with climate. Iterative testing of the gap statistic (Tibshirani et al., 2001) suggests that the data could reasonably be represented by 2, 3, or 5 clusters. Using only two clusters, however, failed to fully capture differences in the climate response seen in Figure 2, while analyses using five clusters resulted in the apparently spurious development of small clusters consisting of only one or two sites and whose climate response could not be differentiated from other clusters. Therefore, we proceed with our analysis for the results using three clusters.
The average climate response of the three clusters is shown in Figures 3–5. These indicate (1) a cluster characterized by a May–June positive precipitation response and a positive seasonal temperature response throughout most of the late prior and current growing season, (2) a cluster with a positive precipitation response in May–June with a negative temperature response during the summer and a negligible temperature response during the spring and winter, and (3) a cluster with a positive May–June precipitation response, a positive winter–spring temperature response, and a subsequent and abrupt negative temperature response during summer.

Median monthly and seasonal correlation and partial correlation response from SEASCORR of the tree-ring chronologies belonging to Cluster 1. Significance levels shown here are the median of those determined by exact simulation (Meko et al., 2011; Percival and Constantine, 2006) for each individual site.

Median monthly and seasonal correlation and partial correlation response from SEASCORR of the tree-ring chronologies belonging to Cluster 2. Significance levels shown here are the median of those determined by exact simulation (Meko et al., 2011; Percival and Constantine, 2006) for each individual site.

Median monthly and seasonal correlation and partial correlation response from SEASCORR of the tree-ring chronologies belonging to Cluster 3. Significance levels shown here are the median of those determined by exact simulation (Meko et al., 2011; Percival and Constantine, 2006) for each individual site.
A map of the distribution of the site cluster assignments (Figure 6) shows that Clusters 2 and 3 are intermingled across the domain, although Cluster 2 contains fewer sites (24, compared with 38 in Cluster 3), is broadly coastal, and is predominantly located in the Levant and Greece. Cluster 3 is the largest and is composed of the majority of the sites throughout Turkey, Cyprus, and Crete. The majority of the sites within Cluster 1 are found in northeastern Turkey.

Cluster assignment for sites in the eastern Mediterranean tree-ring chronology network.
We can examine the distribution of species across our cluster (Figure 7). Cluster 1 contains all the Pinus sylvestris (PISY) sites, but also one Pinus nigra (PINI) chronology, several Juniperus excelsa (JUEX) sites, and a C. libani (CDLI) chronology. Cluster 1 sites are also encountered at significantly higher elevations than Cluster 2 or 3 (Figure 8). For Clusters 2 and 3, the remaining P. nigra and J. excelsa chronologies are distributed between them. The majority of the C. libani chronologies are in Cluster 2 as are all the Abies sites, while Cluster 3 contains all the Pinus brutia (PIBR) chronologies and the Cedrus brevifolia (CDBR) site.

Cluster assignment for species in the eastern Mediterranean tree-ring chronology network. International Tree-Ring Data Bank species codes are Pinus sylvestris (PISY), Pinus nigra (PINI), Juniperus excelsa (JUEX), Pinus brutia (PIBR), Cedrus libani (CDLI), Abies cilicica (ABCI), Cedrus brevifolia (CDBR), Pinus heldreichii (PIHE), and Abies cephalonica (ABCE). Species collected during sampling and not part of the 77 chronologies analyzed here are not shown.

Boxplots showing elevation for chronologies within each cluster. Horizontal lines through each box indicate the median value. The edges of the box indicate the 25th and 75th percentiles, and the whiskers reflect the full range of the site elevation in each cluster. Triangles indicate the comparison intervals for the median (McGill et al., 1978), indicating Cluster 1 is significantly different from Clusters 2 and 3.
The observed positive correlations with summer precipitation and negative correlations with coeval temperatures suggest that integrated measures of water availability could represent an even stronger control on tree growth at some sites in our network. Figure 9 shows the correlation between the chronologies in our network and the local PDSI from Van der Schrier et al. (2013). PDSI incorporates the influence or precipitation, evaporation, and storage into a metric of water balance. Correlations with spring–summer (May through August) or summer (June through August) PDSI in our network are positive over the region, with the strongest correlations seen in western and northern Turkey.

Correlations between tree-ring chronologies and the corresponding local gridded PDSI from Van der Schrier et al. (2013) over their respective common intervals. Colors indicate Pearson product–moment correlation coefficients and those sites where p < 0.10 have an additional black outline. Inset histograms show the frequency (y-axis) of correlation coefficients (x-axis). Panels correspond to the mean May-through-August PDSI (left), mean May and June PDSI (center), and mean June through August PDSI (right).
Discussion
The chronology statistics and cluster analysis attest to the potential value of this 79-chronology tree-ring network for climate reconstruction. Replication and time span are currently sufficient at 76 sites for the chronologies to characterize the common tree-growth signal at the site. More samples could also render the remaining three sites, all in Lebanon (KAM, JRA, and HAD), useful for climate reconstructions.
One of the major objectives of dendroclimatology is to obtain the longest possible tree-ring records from living and dead trees to investigate climate variability over several centuries or longer. Throughout most of the eastern Mediterranean, we have been able to collect samples that were several hundred years in age. This has been accomplished primarily by applying knowledge originally developed over many decades in the semi-arid North American southwest region in the selection of species (distinct annual rings), sites (a single climatic factor limiting to growth), and individual trees (many rings and large year-to-year variability in ring width) for analysis (e.g. Fritts, 1976). Furthermore, multi-century mean segment lengths ensure that series span sufficient time to be adequate for the investigation of multidecadal and centennial climate variability when combined with conservative de-trending and standardization (Cook et al., 1995).
Some elevational dependence is evident in the clusters. Sites in Cluster 1 are at a significantly higher elevation (Figure 8; cf. McGill et al., 1978), consistent with the location of the P. sylvestris that comprise this cluster, while Clusters 2 and 3 are composed of sites at lower elevations but are not distinguishable from one another. The dependence on elevation may reflect a tendency for the species themselves to be stratified by elevation or could also suggest climate regimes determined by topography or elevation (e.g. high mountain) that can impose similar tree-growth variations in chronologies from widely separate locations (Bhattarai and Vetaas, 2006; Lenoir et al., 2008).
Although individual sites can reflect a diversity of monthly or seasonal climate responses, considered across the entire network of 76 sites, a number of consistent features emerge from our cluster analysis: most of the network has an early summer (May–June) positive response to precipitation. Similar results have been reported by Touchan et al. (2003, 2007), Griggs et al. (2007), Akkemik et al. (2008), and Köse et al. (2011), who likewise investigated the relationship between climate and tree-ring data in the region. Hughes et al. (2001) evaluated specimens from living tree chronologies from southeastern Europe to the eastern Mediterranean. They demonstrated that cross-dating over large distances in Greece and Turkey has a clear climatological basis, with signature years consistently being associated with specific, persistent atmospheric circulation anomalies. They reported that the most consistently significant relationship of tree-ring growth to climate data was the positive growth response to spring and early summer precipitation, particularly April–June in the eastern Mediterranean. Touchan et al. (2003, 2007) found that precipitation in May and/or June is most consistently the controlling factor of tree-ring growth in a regional chronology from southwestern Anatolia in Turkey. Touchan et al. (2005b) conducted the first large-scale systematic dendroclimatic sampling for this region using a diversity set of species. Their response function analysis identified May–August total precipitation as the most appropriate seasonal predictand for reconstruction based on consistent patterns of monthly series with significant response function elements of the same sign. Akkemik and Aras (2005) reported that oak tree-ring growth is influenced by March–June precipitation in western Black Sea region of Turkey. Griggs et al. (2007) developed a May–June precipitation reconstruction for northeastern Greece and northwestern Turkey. They reported that negative mean temperature of May and June is also a growth-limiting factor owing to its effect on the availability of precipitation to the trees but is more difficult to calibrate and reconstruct accurately owing to the trees’ indirect response and the low number of long temperature records available for the interior of northwestern Turkey. Köse et al. (2011) conducted a dendroclimatic study of 17 black pine (P. nigra) chronologies in western Anatolia, Turkey, and reported that the influence of precipitation on tree-ring growth was positive and significant in May in almost all regions. A positive but less significant influence was found for June. Their results revealed that the most important factor affecting tree-ring width of black pine was May–June precipitation in western Anatolia.
Temperature responses across our sites are more diverse than those with precipitation but generally fall into one of the three categories: a positive response to winter–spring temperature, a positive response to response in nearly all seasons, and a negative response to temperature during the summer. The clearest difference among the clusters we have identified here is between Cluster 1, characterized by higher elevation P. sylvestris in northeastern Turkey with a positive temperature response through most of the year, and lower elevation sites throughout Turkey, Greece, and the Levant in Clusters 2 and 3 with May–June precipitation responses dominated by P. nigra. Cluster 2 tends to be found near the Mediterranean and has no positive response to temperature, while Cluster 3 is widespread throughout Anatolia and switches from a positive response to winter–spring temperature to a negative response during the summer and coincident with a positive correlation with precipitation.
The PDSI is a common target for paleoclimate field reconstructions (e.g. Cook et al., 1999, 2004, 2010; Touchan et al., 2011) since it is designed to capture the balance between precipitation and evapotranspiration and may therefore be a better indicator of effective moisture conditions than either variable alone. The combination of positive precipitation correlations and negative summer temperature correlations across a swath of our network suggests that the PDSI is a reasonable variable to reconstruct in this region. Positive correlation with spring–summer (May through August) and summer (June through August) PDSIs is observed across our network (Figure 9). Interestingly, the highest correlations are seen throughout western and northern Turkey, where they largely correspond to the sites in Cluster 3. Sites in Cluster 3 have overall correlations with summer PDSI more than twice those of the other clusters, although even Cluster 1 with its broad positive associations with temperature throughout the year shows positive correlations with PDSI. These observations indicate that, in addition to May–June precipitation, either May–August (MJJA) or June–August (JJA) PDSI is a reasonable target for climate field reconstruction in our region, and that PDSI appears to be strongly controlled by precipitation.
Conclusion
This study represents the first large-scale systematic dendroclimatic sampling focused specifically on developing a network of tree-ring chronologies from multiple species from Turkey, Greece, Cyprus, Syria, and Lebanon. This network of chronologies contains coherent seasonal precipitation and temperature signals across a fairly broad geographical area, in particular a common significant response to May–June precipitation. Collectively, these findings suggest that it may be possible to reconstruct three fields from these data. May–June precipitation is an important control on nearly the entire network, while the negative association with summer temperatures suggests that the PDSI is an additional and reasonable reconstruction target. Finally, the existence of a winter–spring and even in some cases summer temperature signal, particularly in P. sylvestris from northeastern Turkey, suggest the possibility of performing temperature reconstructions over at least part of the eastern Mediterranean domain.
Footnotes
Acknowledgements
We wish to thank the different Forestry Departments in Turkey, Greece, Cyprus, Lebanon, and Syria for their great help and support in making this study possible. We would like to thank Cyprus Meteorological Service for providing us with climate data. We thank Professor Alexandros P Dimitrakopoulos from the Laboratory of Forest Protection, School of Forestry and Natural Environment, Aristotle University of Thessaloniki, Greece, for his help and support. We thank Christopher Baisan, Russell Biggs, and Gurudas C Bock for their valuable assistance in the field. We also thank Russell Biggs, Victoria L Frazier, Alicia Stout, Gurudas C Bock, Jessica L Little, and Anthony P Trujillo for their assistance in sample preparation and measurement. We thank Daniel Griffin for developing
. We wish to thank the anonymous reviewers for their constructive comments and suggestions on the manuscript.
Funding
Funding was provided by the US National Science Foundation under Grant Earth System History (grant no. 0075956) and ATM-GEO/ATM-Paleoclimate Program 0758486.
