1. Introduction
Climate is the main driver of natural vegetation and agricultural landscape evolution [
1,
2]. Understanding the effects of climate, both currently and in the near and far future, is one of the most critical aspects of effective planning for landscape conservation, management and agricultural sustainability. The warm-temperate regions of the world, distinctive biomes characterized by dry summers and wet winters (Mediterranean climate), are especially vulnerable to the effects of climatic variability and the change of climate [
3]. This is critically important in the Mediterranean, where a long-term co-evolution of natural and cultural landscape is documented [
4]. Mediterranean Europe provides examples of some of the most complex and heterogenous landscapes and local climatic conditions in the Continent, rendering intriguing problems in climate-driven landscape change prediction [
5] and climate-change adaptation planning. In fact, in recent times, effort has been invested into re-evaluating the historical impact on traditional landscape development; it has long been said that the euro-Mediterranean is too complex to allow for simple generalizations [
6].
The anthropogenic induced climate change (CC) ranks high in the hierarchy of scientific interests owing to its expeditious spatiotemporal development and menacing nature [
7]. The associated modified bioclimate poses considerable threats to the survival of living organisms, the sustainability of natural landscapes, the viability of agricultural activity, the conservation of natural resources, the maintenance, improvement, deterioration or loss of ecosystem services, and the efficiency of management policies [
8,
9,
10,
11,
12,
13,
14].
In the coming decades, Greece will appear as climatically pressurized owing to the marked transformation of the bioclimate that is already arising over the entire territory [
15,
16,
17]. CC projections demonstrate notable warming under the impact of the extreme RCP8.5 climatic scenario [
18,
19]. By the year 2100, the near-surface temperature is expected to rise by an average value of 4.3 °C, resulting, thus, in a significant increase of the night frosts, continual dry spell days, the annual number of hot days and tropical nights, the growing season length and a reduction of frost days. The foreseen increase in the annual number of consecutive dry days by 30% (15.4 days) and the concomitant decline by 16% to 40% of the seasonal precipitation, underline a crucially drier bioclimatic footprint in the years to come [
20,
21,
22,
23]. Additionally, in the overall concept of both the projected and observed CC, the present and expected extreme weather events appear as more frequent and intensified (e.g., more severe and prolonged drought events; earlier start and later ending of heatwaves period; increasing trends of heavy precipitation), highlighting the climate crisis impacts on human security [
24,
25,
26,
27].
Climate projections over Greece in conjunction with its present rather dry-thermal climatic regime give grounds for the country’s escalating susceptibility to the phenomenon of CC which may induce far-reaching impacts on its extensive natural ecosystems and highly heterogeneous agricultural sector [
28,
29].
Several present and foreseen effects of CC on the natural landscapes of Greece include significant habitat alterations and elevated extinction risk for autochthonous species and vegetation communities (see e.g.,
Juniperus drupacea forests of Greece) [
30], increase in drought-related tree dieback [
31], declining tree productivity and development [
32], reduction in habitat-suitable areas [
33], projected altitudinal and elevational shifts and extinctions of endemic species [
34,
35], spatial redistribution of species, invasions and biodiversity loss [
36], increased likelihood of severe wildfires and floods [
37,
38] and projected altered fire behavior [
39].
Some of the CC impacts related to rural and agricultural areas may involve the northward shift of the agro-climatic zones [
40], crop production quality decline and quantity limitation [
17,
41], effects on the capacity of crops to adapt [
42], increased susceptibility of crops to extreme weather [
43], differentiations in the geographical expansion of cultivations [
44], alterations in the cultivation area suitability [
45], soil erosion and risk for land desertification [
46], impacts on groundwater quantity and quality, limitation of surface water availability [
47,
48], changes in crop phenology [
49], increased growing season length [
50] and negative social and economic development trends in rural areas [
51].
Temperature and precipitation are fundamental climatic parameters and are decisive inputs for the investigation of CC. However, the overall apprehension and magnitude of CC and the comprehension of its impacts on vegetation and vegetation types’ distribution is far more expressed and achieved by applying bioclimatic indices [
52,
53,
54]. The bioclimatic indices’ major significance as tools for categorizing bioclimates is indisputably exhibited by their broad exploitation in climatological, bioclimatological, and agricultural research studies [
55,
56,
57].
Approximately one hundred years ago, the de Martonne bioclimatic index (DMI) was proposed for the evaluation of a specific environment’s dryness degree by categorizing it into seven (7) classes, from “dry” to “extremely humid” based on the fundamental climatic parameters of air temperature and precipitation [
58,
59]. The DMI has been extensively implemented for the classification of the bioclimate due to its reliability, effectiveness, and validity [
60,
61]. Scientific investigations, that frequently employ the DMI, fall within the fields of climatology/bioclimatology and agricultural and land or water resources management, while the index also appears as a serviceable tool for achieving environmental assessment reports [
60,
62,
63,
64,
65].
The Emberger index (EI), commonly described as the “pluviothermic quotient Q”, categorizes the Mediterranean area’s bioclimate zones corresponding to a scheme extending from the “Per-Humid” to the “Per-Arid” type of bioclimate or bioclimatic category established on the parameters of temperature, precipitation, and evaporation. Estimations of the EI include the representation of the annual temperature by the average maximum temperature value of the hottest month (M) and the average minimum value of the coldest month (m), considering that the growth of vegetation appears dependable on these thermal limits. Precipitation (P) is expressed on an annual basis. At the same time, evaporation is indirectly represented by the difference between the two temperature values (M—m), owing to the EI’s increase with the latter parameter [
66,
67,
68]. Also, a simplified algorithm based on the minimum winter temperature (m) falling within the range of the “Very Hot” to the “Very Cold” temperature characterizations is applied, serving the purpose of the phytoclimate classification in bioclimatic subtypes, often termed as the ‘Q2’ [
68,
69]. Thus, an area’s phytoclimatic footprint is characterized through the combination of the bioclimatic types of characterizations derived from the estimates of the Q values along with the temperature conditions based on the estimates of the m values, which results in the Emberger’s Q2 bioclimatic subtypes; for example, a Q2 subtype may be described as the “Semi-arid, Mild winter” subtype.
Within the changing climate’s research framework, studies involving the implementation of the DMI and EI indexes are rather limited in Greece. For example, Baltas [
70] has denoted the variability of aridity in Northern Greece from 1965 to 1995 as derived from the DMI characterization range of the Semidry to the Very humid bioclimate. The same author estimated the bioclimatic footprint of the entire country which was described by the Dry to the Very humid classes. Another study in Northern Greece conducted by Mattas et al. [
71] on the upper agricultural area of the Gallikos river basin has revealed bioclimatic variations from the Semiarid to the Humid conditions based on climatic data of a 27-year period (1980–2006). Also, mountainous regions have been classified by Sidiropoulou et al. [
72], who have characterized respectively Northern Greece’s Mt Vermio and Southern Greece’s Mt Zireia as moderately humid and humid, respectively on the basis of climatic data during the years 1990–2019. Lappas et al. [
73] have revealed the bioclimatic conditions of a water district in eastern–central Greece illustrated by the Arid to Very humid classes over the 1980 to 2001 timeframe. A moderate trend, between the years 1958 and 2011, towards a more dry-thermal regime over the traditional agroforest system of Thriasio Plain (Northwest Athens) has been exhibited by Mavrakis et al. [
74]. Further south, Beloiu et al. [
75] pinpoints the transition of Crete Island’s bioclimate from Very humid to Humid conditions over the 1979–2013 period.
According to the recent findings of Charalampopoulos et al., 2023 [
19] for the Greek territory, by 2100 and under the RCP8.5, the semidry bioclimatic conditions are expected to prevail over a significant percentage natural landscape area of approximately 34%, followed by the Mediterranean regime (percentage area of nearly 19%) and the most xerothermic dry bioclimate (very limited area of 0.1%). Also, profound trends towards the more dry-thermal regime are exhibited for the agricultural areas where the semidry conditions are projected to influence substantial parts under both the RCP7 (area of 53.9%) and RCP8.5 (area of 60%), by the end of the century. The authors point out that the sustainability of the ‘Semidry’ agricultural areas will depend entirely on water supply through appropriate irrigation schemes, while the areas with Mediterranean bioclimatic characteristics (nearly 20% for both scenarios) are expected to satisfy their water needs through supplementary irrigation. These evolutionary situations raise uncertainty for the future of the agricultural sector by considering severe impacts on Greece’s production capacity.
Evidently, most utilizations of the DMI concern bioclimatic studies based on past timeframes, while there has been no application of the EI until recently and over the entire country involving past, present, but also future time-periods [
18,
19]. This highlights the limited scientific knowledge on the future bioclimatic development on a local scale over the natural and agricultural areas, with no applications of the DMI and EI as tools for bioclimatic categorization and for projections on conditions foreseen on a local scale in Greece, owing to CC.
Motivation, therefore, for the conduction of the present investigation lies in the relatively limited scientific outcomes related to the DMI and EI applications, along with the critical need for the protection and preservation of the natural and agricultural landscapes already threatened by the changing bioclimatic conditions.
More specifically, in this study, both the bioclimatic DMI and phytoclimatic EI indices are computed. Their categories are spatiotemporally illustrated, for the first time, for the present (1970–2000) and future (2030–2060; 2070–2100) climate conditions at very high resolution (~500m) and under two specific greenhouse gas emission scenarios, namely the RCP4.5 and RCP8.5.
Overall, and under the framework of the LIFE-IP AdaptInGR National Project, the main objective of the present study is to outline the major short-term and long-term bioclimatic alteration trends that are foreseen to occur and shall undoubtedly impact the already pressurized natural and anthropogenic ecosystems in Northern Greece. The resulting mapping of the spatial distributions of the indices’ classes, coupled with the spatial statistics per class, period, and emission scenario, may serve as useful tools for conducting further environmental research in the broader framework on the conservation of landscapes threatened by CC. The final outcomes aim to support decision and policy making for regional natural capital, ecosystem services, and agricultural development management to adapt in, and mitigate future predicted conditions.
4. Conclusions
The study outcomes outline the bioclimatic change over the natural and agricultural landscapes of the Region of Eastern Macedonia and Thrace, Northeastern Greece. Thus, the present study’s innovative findings resulting from the application of both bioclimatic indices (de Martonne index; Emberger index) may be summarized by the following major conclusions:
(a) Profound altered bioclimatic regimes predict trends towards more xerothermic environments under the two examined emission scenarios RCP4.5 and RCP8.5. Comparisons between the three study time frames (1971–2000; 2031- 2060; 2071–2100) highlight the temporal development in the direction of the more dry-thermal conditions, as predicted under both emission scenarios.
(b) Significant changes following the implementation of the de Martonne index concern the long-term timeframe of the extreme RCP8.5 scenario. Owing to high emission regimes and associated global warming by the end of the century, the more xerothermal Mediterranean and Semi-humid conditions will account for approximately 40% of the investigated area. The parts to be particularly influenced consist of the eastern, southern and western agricultural areas which are projected to require increased supplementary irrigation. Natural areas may also be increasingly stressed due to reduced humidity and precipitation regimes.
(c) As documented by the Emberger index, the spatiotemporal xerothermic evolution over the investigated area follows a North-South direction. An intensified xerothermic trend under the extreme RCP8.5 is also projected by the Emberger index, in the long run. By the year 2100, a substantial part of the natural and agricultural areas amounting up to nearly 42% is projected to be impacted by the Sub-humid, Mild winter conditions and up to nearly 33% by the Sub-humid, Cool winter conditions. The possible existence of the Warm winter combinations expected to impact on both agricultural and natural landscapes (distribution of nearly 5%) underlines the peaking of the bioclimatic shift towards warmer and drier conditions.
(d) Finally, a multi-scaled approach at the regional level should be helpful in planning for natural capital conservation, restoration and sustainability initiatives at the landscape level. This particular studied area is so diverse and complex and has several high-profile and well-studied protected natural areas or other designated areas of outstanding interest, both in its natural/semi-natural and cultural landscapes. However, new pressures from climate-associated stress may create severe problems associated primarily with water-related ecosystem services, i.e., water stress in agricultural and aquatic, wetland and riparian areas accompanied by wildfire frequence increase. Climate adaptation analysis would require a scientifically-informed and guided landscape approach and should not be relegated only to specific landscapes or high-profile protected areas as has already been done in some recent initiatives. Also, it is important to carefully assess where adaptation and mitigation measures may conflict with nature conservation initiatives. It is recommended that landscape-scale approaches should be taken in promoting climate adaptation and climate-smart planning initiatives throughout the Region of Eastern Macedonia and Thrace.