DOI: 10.20937/ATM.53272

Received: March 12, 2023; Accepted: November 13, 2023

Trend analysis and forecast of annual precipitation and temperature series in the Eastern Mediterranean region

RESUMEN

El objetivo del presente estudio es examinar a nivel nacional la tendencia actual de las series anuales de precipitación y temperatura referidas a la cuenca del Mediterráneo Oriental, incluyendo la variabilidad anual de las anomalías de la temperatura superficial del mar y el efecto de conexión con los índices de la Oscilación del Atlántico Norte y la Oscilación del Mediterráneo (NAOI y MOI por sus siglas en inglés, respectivamente). El periodo considerado abarca principalmente los últimos 32 años, desde 1990 hasta 2021. Adicionalmente, la predicción de las series de temperatura y precipitación mensuales y anuales, basadas en el Modelo Autorregresivo Integrado de Media Móvil sobre una base nacional, para los próximos cuatro u ocho años, es otra característica interesante del presente estudio. Los resultados indican tendencias crecientes en la precipitación anual de 1990 a 2021, en su mayoría tendencias de calentamiento significativas y no significativas en la temperatura anual, incluida la temperatura superficial del mar. El índice NAOI está altamente correlacionado principalmente con la temperatura anual, mientras que el índice MOI no afecta la variación de la precipitación y temperatura anuales. Los resultados del presente estudio coincidieron en general con aquellos de los estudios disponibles en la literatura, y podrían suponer un alto interés para autoridades nacionales y gremios u organizaciones ambientalistas, como instrumento para ayudar a los encargados de la toma de decisiones a enfrentar el cambio climático.

ABSTRACT

The present study aims to examine the current trend of annual precipitation and temperature series referred to the Eastern Mediterranean basin on a national basis, including the annual variability of sea surface temperature anomalies and the connection effect with the North Atlantic Oscillation and the Mediterranean Oscillation indices (NAOI and MOI, respectively). The period under consideration is mainly the last 32 years, from 1990 to 2021. Additionally, the prediction of monthly and yearly temperature and precipitation series, based on the autoregressive integrated moving average model on a national base, for the next four or eight years, is another interesting feature of the present study. Results indicate rising trends in the annual precipitation during 1990-2021, mostly non-significant and significant warming trends in the annual temperature, including sea surface temperature. The NAOI is highly correlated mostly with annual temperature, whereas the MOI does not affect the variation of annual precipitation and temperature. The results of the present study are in general agreement with the results of available studies in the literature, and they could be of high interest to national authorities and environmental unions/organizations, to help decision-makers face climate change.

Keywords: time series, ARIMA, CRU TS 4.05, climate change impact, climate variability, Greece, Turkey, Cyprus, Israel, Lebanon.

1. Introduction

The eastern part of the Mediterranean Sea and its continental areas extend from 18.80-43.5º N and 29.9-41.7º E, including Greece, Turkey, Cyprus, Israel, Lebanon, Egypt, Syria and Libya. The study area (Fig. 1) is economically important, being exploited for fisheries, oil and gas extraction, and offshore renewable energy; it is also used for shipping, which links its many harbors. In the coastal zone, tourism, road transportation, recreation, agricultural and livestock production are major uses. The Mediterranean basin is characterized by strong spatiotemporal contrasts in the precipitation pattern, with high variability at multiple timescales often enhanced by orography, associated with large inter-annual fluctuations and a high heterogeneity of precipitation regimes, terrain, and land surface properties (de Luis et al., 2010; Piras et al., 2016; Cid et al., 2017). The eastern part of Mediterranean area is considered one of the most vulnerable regions to climate change (Giorgi, 2006; IPCC, 2007; Lu et al., 2007; Zittis et al., 2022). The Mediterranean Sea was possibly the first of the world’s seas where changes in temperature were attributed to global warming (Bethoux et al., 1990).

Fig. 1. The Eastern Mediterranean Region.

Precipitation and temperature are the most important climatic variables in the hydrology/climate/atmosphere fields and in water resources due to their critical effect on the temporal and spatial patterns of water availability, transpiration, and evaporation. Studies to detect trends in temperature and precipitation series are essential to analyze the impact of climate change in such areas. Several studies have analyzed temperature, precipitation, or extreme temperature-precipitation trends in the Eastern Mediterranean area (Sahsamanoglou and Makrogiannis, 1992; Xoplaki et al., 2004; Kostopoulou and Jones, 2005; Zhang et al., 2005; Kuglitsch et al., 2010; Philandras et al., 2011; Lelieveld et al., 2012; Caloiero et al., 2018; Hafez, 2018; Mathbout et al., 2017), and also sea surface temperature (SST) trends and variability in the Mediterranean basin/sub-basins (Nykjaer, 2009; Shaltout and Omstedt, 2014; Androulidakis and Krestenitis, 2022).

On a national basis, these trends have been addressed by the studies of Pnevmatikos and Katsoulis (2006), Feidas et al. (2007), Philandras et al. (2008), Kambezidis et al. (2010), Nastos et al. (2011), Mamara et al. (2015), Feidas (2017), Tzanis et al. (2019), Kastridis et al. (2022), Sakalis (2022), Varlas et al. (2022), Stathi et al. (2023), referred to Greece; Kadioğlu (1997), Türkeş (1999), Türkeş et al. (2002), Türkeş and Erlat (2003), Partal and Kahya (2006), Hadi and Tombul (2018), Topuz et al. (2020), Mersin et al. (2022), Todaro et al. (2022), Aalijahan et al. (2023), referred to Turkey; Ben-Gai et al. (1999), al-Rimmawi et al. (2010), Givati and Rosenfeld (2013), Shadeed (2013), Yosef et al. (2019), Ajjur and al-Ghamdi (2021), Salameh et al. (2022), referred to Israel; Price et al. (1999), Giannakopoulos et al. (2010), Zaifoğlu et al. (2017), Zittis et al. (2021), referred to Cyprus; and Shaban (2011) and Ramadan et al. (2013), referred to Lebanon.

In some of the above-mentioned studies, yearly teleconnection indices (the North Atlantic Oscillation Index, NAOI, and the Mediterranean Oscillation index, MOI) were both analyzed in accordance with the temperature and precipitation variability on a national basis or within the Mediterranean region/sub-regions.

Although there are numerous studies on a national or regional basis referred to trend detection of the main weather variables, straight comparisons between these are not always representative. One of the main purposes is the examined periods, which differed between each other from some years to more than a decade. Statistical methods for trend detection, such as the non-parametric Mann-Kendal (Mann, 1945) and the Sens slope estimator (Sen, 1968) are affected by the length of the existing data series (Burn and Elnur, 2002). Otherwise, collecting data from weather stations is often a difficult and time-consuming procedure, especially when these stations belong to different countries, and different authorities are responsible for their operation. Additionally, the quality of ground-based weather station series sometimes is low, often with the absence of the relevant metadata and possible data gaps (missing values), relocations, calibration issues, etc., and sometimes homogenization (Gofa et al., 2019) and quality control is needed.

An issue that may affect the accuracy of results in the Eastern Mediterranean basin is the low spatial density of the existing ground-based weather stations that were used in several studies, especially in the mountainous areas of Greece and Turkey, a possibility which led to a range of unknown scale of microclimates. A way to overcome these possible issues is the use of results provided by several existing and reliable weather variables databases, such as the Climate Research Unit Time series (CRU TS) database from the University of East Anglia in the United Kingdom (the current version is 4.05), which covers the global region for the period 1901 to the present (Harris et al., 2020). Relevant data are based on high-quality weather stations, interpolated in 0.5º latitude-longitude grid cells. By using gridded data, climate models can be evaluated more directly since they represent spatial area averages rather than point data; moreover, averaging over regions is straightforward and, finally, estimated data become available for non-observed locations. Several studies have used the CRU-TS datasets for trend detection of climatic variables in the Eastern Mediterranean basin on a national or regional basis, such as those from Philandras et al. (2011), Ramadan et al. (2013), Hadi and Tombul (2017, 2018), and Sakalis (2022).

Accurately predicting trends of precipitation and temperature can play an important role in irrigation scheduling, use of water resources, social planning, human health, production programming, environmental and civil protection, etc. Forecasting methods such as the autoregressive integrated moving average (ARIMA; Box and Jenkins, 1976) are often used to forecast several variables in meteorology, atmospheric physics, hydrology, and also in marketing, economical sciences, trading, human health, epidemiology, medicine, etc. Non-seasonal ARIMA and/or seasonal ARIMA (SARIMA) methods were recently successfully implemented to forecast future precipitation, temperature, or evapotranspiration variables by Soltani et al. (2006), Wang et al. (2014, 2021), Feng et al. (2016), Mahmud et al (2017), Rahman et al. (2017), al-Balasmeh et al. (2019), Coban et al. (2021), and Tadesse and Dinka, (2022).

The aim of the present study is to examine the current trend of annual precipitation and temperature series in the Eastern Mediterranean basin on a national basis, including the annual variability of SST anomalies and the connection effect with the NAOI and MOI. The period under consideration is mainly the last 32 years, from 1990 to 2021. Another interesting feature of the present study is the prediction of monthly and yearly temperature and precipitation series on a national basis for the next four or eight years. In the available literature there is not a recent and complete study obtained on a national basis, that either includes trend detection in common periods, predicted annual and monthly results, and the examination of the correlation with the teleconnection patterns for the Eastern Mediterranean basin. The results of the present study could be of high interest to national authorities and environmental unions/organizations aiming to understand the current and future influence of climate change phenomena, by aiding them to organize and improve the existing corrective actions for reducing its negative effects on the society and environment.

2. Data and methods

2.1 Study area

The study area (Fig. 1) comprises the eastern part of the Mediterranean Sea. Its continental regions extend from 18.80 to 43.5º N and from 29.9 to 41.7º E. The area includes Greece, Turkey, Cyprus, Israel, and Lebanon. The selection criterion for countries that belong to the area under consideration is the percentage (> 30%) of their perimeter washed by the Mediterranean Sea. Based on this criterion, Egypt, Syria, and Libya were excluded from the study. Another reason for this exclusion was that the main territories of these particular states are extended far away from the Eastern Mediterranean (EM) basin and they are mostly affected by different climatic phenomena (Northern Africa and Western Asia areas have a mostly dry climate).

According to the Köppen-Geiger’s climate classification system, the climate of the EM basin (Fig. 1) is characterized mainly as hot-summer Mediterranean climate (temperate Csa, Cfa, Cfb), with the north-eastern parts of Turkey as continental (Dsa, Dfa) and the southern part of Israel as dry (Bwh).

2.2 Historical weather data

2.2.1 Precipitation and Temperature

Country-based long-term monthly precipitation and temperature data records were obtained from the gridded monthly data of the CRU TS v. 4.05 (Harris et al., 2020), downloaded from the Climate Change Knowledge Portal (CCKP, 2021). for the period 1990-2021 (the last 32 years). Considering that data records form CRU TS v. 4. are available from 1901 to the present, and for comparison with the available results of other studies, the period under consideration spans appropriately towards the past.

2.2.2 North Atlantic Oscillation Index (NAOI)

According to Feidas et al. (2007), the North Atlantic Oscillation (NAO) is one of the large-scale modes of climate variability in the Northern Hemisphere. It defines the large-scale meridional oscillation of atmospheric mass between the center of the subtropical high surface pressure located near the Azores and the sub-polar low surface pressure near Iceland.

Strong positive phases of the NAO tend to be associated with above-normal temperatures in the eastern United States and across northern Europe and below-normal temperatures in Greenland and oftentimes across southern Europe and the Middle East. They are also associated with above-normal precipitation over northern Europe and Scandinavia and below-normal precipitation over southern and central Europe. Opposite patterns of temperature and precipitation anomalies are typically observed during strong negative phases of the NAO. During particularly prolonged periods dominated by one particular phase of the NAO, abnormal height and temperature patterns are also often seen extending well into central Russia and north-central Siberia (Walker, 1924; Walker and Bliss, 1932; van Loon and Rogers, 1978; Rogers and van Loon, 1979).

The NAOI data used for the present study were based on Hurrell et. al. (2003), and provided by the Climate Analysis Section of the National Center of Atmospheric Research (https://www.cgd.ucar.edu/sections/cas). The data analyzed span mainly from 1990 to 2021. NAOI data records are available from 1899 to present.

2.2.3 Mediterranean Oscillation Index (MOI)

The MOI has been defined by Conte et al. (1989) and Palutikof et al. (1996) as the normalized pressure difference between Algiers (36.4º N, 3.1º E) and Cairo (30.1º N, 31.4º E). It is described as a teleconnection pattern with opposite pressure and rainfall anomalies between the Central-Western and Eastern Mediterranean areas. It explains the annual precipitation variability across the Central-Western Mediterranean basin better than the NAOI. The MOI is constructed from the 500 hPa surface height monitored by soundings at stations in Algiers and Cairo is often assumed to be representative of the Western and Eastern basin, respectively. When the pressure increases in the Western basin a decrease is found in the East, and vice versa.

The MOI data used for the present study were provided by the CRU, and analyzed mainly for the period from 1990 to 2021. The MOI data records are available from 1948 to the present.

2.2.4 Sea Surface Temperature (SST)

SST data from the Mediterranean Basin based on satellite measurements for the period 1990-2014 were provided by the European Environmental Agency (EEA, 2012). SST data records referred to the Mediterranean basin are available from 1870 to 2014.

2.3 Methodology

2.3.1 Homogeneity test of the data

Meteorological observations are often subject to influences caused by instrument changes or other reasons (i.e., relocation of the observatory, change in the observation times, missing data, methods used to calculate the daily means, and increased urbanized and or industrialized areas), which can introduce artificial trends in time series and render their data unrepresentative of the actual climate variation. The four well-known nonparametric methods selected to test the departure from homogeneity in time series are the standard normal homogeneity test (SNHT) for a single break (Alexanderson, 1986), The Buishand range test (Buishand, 1982), the Pettitt test (Pettitt, 1979) and the von Neumann ratio test (von Neumann, 1941).

The CRU TS v. 4.05 data analyzed in the present study (see section 2.2.1) are based on results from selected weather stations with high quality standards, so the possible detected inhomogeneities could probably occur due to real climatic conditions rather than artificial influences.

2.3.2 Trend detection

The auto-correlated Mann-Kendall test (Mann, 1945; Kendall, 1955; Hamed and Ramachandra, 1998) with Sen’s slope estimator method (Sen, 1968) is used to eliminate the autocorrelation and to verify the presence of a trend in precipitation and temperature time series data. The method was used for annual precipitation and temperature data as well as for SST and NAOI/MOI data.

2.3.2.1 Auto-correlated Mann-Kendall test

The Mann-Kendall test is a nonparametric test for monotonic trend detection. It does not assume the data to be normally distributed and is flexible to outliers in the data. The test assumes a null hypothesis (H₀) of no trend and an alternate hypothesis (H_a) of increasing or decreasing monotonic trend. For a time series X_i = x₁, x₂,..., x_n, the Mann-Kendall test statistic S is calculated as

(1)

where n is the number of data points, x_i and x_j are the data values in time series i and j (j > i), respectively, and sign (x_j – x_i) is the sign function calculated as

(2)

Statistic S is normally distributed with parameters E(S) and variance V(S) as given below:

(3)

where n is the number of data points, m is the number of tied groups, and t_k denotes the number of ties of extent k. The standardized test statistic Z is calculated using Eq. (4):

(4)

To test for a monotonic trend at an α significance level, the alternate hypothesis of trend is accepted if the absolute value of standardized test statistic Z is greater than the Z_{1 – α/2} value obtained from the standard normal cumulative distribution tables. A positive sign of the test indicates an increasing trend and a negative sign indicates a decreasing trend.

According to the variance correction approach proposed by Hamed and Ramachandra (1998), recently used by Patakamuri et al. (2020), when the time series exhibit a significant serial correlation, the following modified variance V(S)^* of the Mann-Kendall test is used:

(5)

The proposed correction factor CF is

(6)

where r_k and r_k^R are the lag − k serial correlation coefficients of data and ranks of data, respectively and n is the total length of the series. In the case of the correction factor CF, only significant correlation coefficients are used.

2.3.2.2 Sen’s slope estimator

The Sen’s slope method (Sen, 1968) is a robust nonparametric method of estimating the magnitude of a trend slope. For a given time series X_i = x₁, x₂,..., x_n, with N pairs of data, the slope is calculated as:

(7)

Median of N values of β_i gives the Sen’s estimator of slope β

(8)

2.3.2.3 The (S)ARIMA forecast models

The ARIMA model, one of the Box et al. (1976) methods, is an important statistical approach in forecasting and analyzing time series data related to precipitation and temperature (Soltani et al., 2006; Wang et al., 2014, 2021; Feng et al., 2016; Rahman et al., 2017; Mahmud et al., 2017; al-Balasmeh et al., 2019; Coban et al., 2021; Tadesse and Dinka, 2022). It is more common to label it as ARIMA(p, d, q), where p is the order of the autoregressive process, q is the order of the moving average process, and d is the differentiation process. In general, the ARIMA(p, d, q) model for time series data can be expressed as:

(9)

where φ and θ are the autoregressive and moving average parameters, respectively, B is the backward shift operator, Z_t the observation at time t, θ the moving average parameter, and ε_t the residual of the model at time t.

The Seasonal ARIMA model (SARIMA), proposed by Box and Jenkins (1976), was used for model building and forecasting monthly rainfall and temperature. It can be labeled as ARIMA(p, d, q) (P, D, Q), where (p, d, q) is the nonseasonal part and (P, D, Q) is the seasonal part of the model, which could be written as:

(10)

where Φ and Θ are the seasonal autoregressive and moving average parameters, respectively, S is the seasonal difference, and the differencing operator.

In the present study, the appropriate combination of parameters above was performed in order to choose the best model for all cases, as described in the next section.

The evaluation of methods was performed using the XLSTAT package. Diagrams were performed with the Origin Pro V.8 package.

2.3.2.4 Model performance criteria

The performance efficiency of the predicted models can be evaluated using the statistical indices mean absolute error (MAE), root mean square error (RMSE), and the Nash-Sutcliffe efficiency coefficient (NSE), defined as

(11)

(12)

(13)

where N is the number of the observations, and y_i, and are the observed, the predicted and the mean value, respectively.

3. Results and discussion

3.1 Homogeneity test of the data

In the present study, homogeneity tests were performed on both precipitation and temperature data series. The time period in which these tests were implemented is the complete available set of yearly precipitation and temperature of the CRU TS 4.05 data series (from 1901 to 2021). Whenever at least two of the homogeneity tests analyzed in section 2.3.1 detect an inhomogeneity (95% confidence level [c.l.]) in the same year, this year is evaluated as a possible change point of the series (Schonwieseand Rapp, 1997; Wijngaard et al., 2003).

In the case of the CRU TS 4.05 precipitation series, no possible change point was detected, whereas, for the CRU TS 4.05 temperature series, the possible change points (1997 for Greece and 1993 for all the other countries) are presented in Table I. It must be noted that the SST series were also subjected to the homogeneity test, without any possible change point being detected.

Table I. Trend detection results for precipitation, temperature series and MOI and NAOI.

Trend detection results
		Period	Year of change point	Mann-Kendal coefficient Z	Significance (p)	Sen’s slope (mm yr-1)
Precipitation	GREECE	1990-2021		2.29	(*)	4.356
	CYPRUS	1990-2021	-	0.76	0.223	1.523
	ISRAEL	1990-2021	-	0.18	0.403	0.310
	TURKEY	1990-2021	-	0.73	0.203	1.121
	LEBANON	1990-2021	-	0.34	0.307	0.738
Indices	NAOI	1990-2021 1955-2021 1955-2013	-	–0.99 2.57 1.66	0.161 (*) (+)	–0.022 0.031 0.014
	MOI	1990-2021	-	-0.70	-	–0.001
Temperature	GREECE	1990-2021	1997	4.69	(***)	0.047
	CYPRUS	1990-2021	1993	4.75	(***)	0.059
	ISRAEL	1990-2021	1993	4.69	(***)	0.047
	TURKIE	1990-2021	1993	4.65	(***)	0.047
	LEBANON	1990-2021	1993	4.49	(***)	0.047
	SST	1990-2014	-	5.00	(***)	0.042

MOI: Mediterranean Oscillation Index; NAOI: North Atlantic Oscillation Index.
Bold characters denote significant trends at a 90% confidence level (+), 95% (*), 99% (**), and 99.9% (***).

According to Hadi and Tombul (2018) 1997 was also detected in their study as a change point in the Turkey temperature series. In the case of Greece, 1997 and 1998 were frequently denoted as years of starting trends by Mamara et al. (2015) and Feidas (2017).

As stated in section 2.3.1, the detected inhomogeneities of the CRU data series could occur due to real climatic conditions rather than artificial influences.

3.2 Precipitation trends

Table I presents the series of annual trend precipitation for the EM region countries, together with the NAOI and MOI. Results are referred mainly to the period 1990-2021. For comparison with the available results of other studies, the period under consideration spans appropriately towards the past for Greece, Israel, Turkey, Cyprus (Table II), and the NAOI (Table I).

Table II. Results of annual precipitation and temperature trends compared to other studies available in the literature.

Present study Literature
		Period	Mann-Kendal coefficient Z	Significance	Trend (Sen’s slope) (mm y–1 or ºC y–1)	Available studies	Trend (mm y–1 or ºC y–1)	Area-averaged data based on
Precipitation	GREECE	1955-2001 1951-2009 1950-2020	–3.39 –2.97 –1.91	(***) (***) (*)	–3.512 (5.002 % decade–1) –2.23 –1.086	Feidas (2007) Filandras (2011) Varlas et al. (2022)	(5±2% decade–1) (*) – –1.02 (+)	WS*(22) CRU 3.01 ERA 5
	CYPRUS	1985-2015	0.37	–	0.916	Zaifoğlu et al. (2017)	–	WS (37)
	ISRAEL	1950-2009	–1.18	(p=0.129)	–0.595	Givati (2013)	–0.701 (*)	NCEP/NCAR Reanalysis
	TURKEY	1955-2013	0.51	–	+0.354	Topuz (2020)	+0.9	W.S (29)
	LEBANON	–	–	–	–	–	–	–
Temperature	GREECE	1979-2017		(***)	0.038	Tzanis et al (2019)	0.034 (**)	ERA INTERIM Reanalysis
	CYPRUS	1901-2000	2.37	(*)	0.005 (0.5 ºC century–1)	Price et al. (1999)	(≈1 ºC century–1)	WS (2)
	ISRAEL	1988-2017	4.42	(***)	0.048 (0.48 ºC decade–1)	Yosef et al. (2019)	(0.53 decade–1) (*)	WS (34)
	TURKIE	1955–1989	-1.16	(**)	-0.018 (-0.63 ºC)	Cadioglu (1997)	(-0.81 ºC)	WS (17)
	LEBANON	1963–2009	3.83	(***)	0.030 (1.38ºC)	Shaban (2011)	(1.8 ºC)	WS

W.S* referred to Weather Station.
Bold characters denote significant trends at a 90% confidence level (+), 95% (*), 99% (**), and 99.9% (***).

The annual variation of precipitation series, together with the relative linear trend (Sen’s slope), including the NAOI and MOI, is visualized in Figure 2.

Fig. 2. Precipitation series (a-e) with relative trends (Sen’s slope), (f) NAOI, and (g) MOI during the period 1990-2020. (NAOI: North Atlantic Oscillation Index; MOI: Mediterranean Oscillation Index).

During the period 1990-2021 (Table I), all the observed annual precipitation trends in the EM region are positive, whereas only in Greece the trend is significant at the 95% c.l. These positive annual precipitation trends span form 0.310 mm yr^–1 (Israel) to 4.36 mm yr^–1 (Greece).

The results of the present study in various periods (Table II) are in close agreement with the results of the available literature referred to gridded monthly data (Givati and Rosenfeld, 2013; Varlas, 2022), and in a general qualitative agreement with the results of Philandras et al. (2011), who based their study on CRU 3.01 gridded data, but did not provide an area-averaged value of the relative trend.

The results of the present study are in general agreement with the result of the available literature referred to area-averaged weather station data (Feidas et al., 2007) and in general qualitative agreement with the results of Zaifoğlu et al. (2017), which are based on 37 weather stations of northern Cyprus (these authors do not provide the area-averaged data). However, the results of the present study (referred to Turkey) are in slight agreement with those of Topuz et al. (2020).

According to the present study, during the period 1990-2021 the observed annual trend of the NAOI (Table I, Fig. 2) is negative but not significant (p = 0.161). This result indicates that during the period 1990-2021, the NAOI trend is weakened, resulting possibly in warmer and wetter conditions within the influenced areas. The increased positive precipitation trend in Greece (2.29 mm yr^–1, p = 0.161) and Turkey (0.73 mm yr^–1, p = 0.203) during this period could partially be explained by the relative negative trend of the NAOI. The reversing sign of the NAOI trend for the last 30 years is a remarkable observation, in spite of the strong positive trends of the periods 1955-2021 and 1955-2013.

Table III presents the correlation coefficients between the precipitation series of the EM with the NAOI and MOI. The series are referred mainly to the period 1990-2021. Additionally, sub-periods between 1951 and 2020 are also examined for comparison with other studies.

Table III. Correlation between the series of temperature and precipitation with NAOI and MOI.

Correlation between series
	Country	Period	ΝΑΟΙ		MOI		SST*
			Correlation coefficient (R)	Probability (p)	Correlation coefficient (R)	Probability (p)	Correlation coefficient (R)	Probability (p)
Precipitation	Greece	1990-2021 1955-2001 1951-2010 1955-1995	–0.295 –0.44 –0.373	0.101 0.002 0.0033	–0.268 –0.19	0.138 0.230	0.50	0.011
	Cyprus	1990-2021 1955-2001	–0.017 –0.193	0.929 0.19	0.183	0.316	0.17	0.42
	Israel	1990-2021 1955-2001 1970-2020	0.331 0.266 0.214	0.065 0.0703 0.132	0.237 +0.27	0.192 0.059	0.24	0.26
	Turkey	1990-2021 1955-2013	–0.278 –0.40	0.124 0.00158	–0.159 –0.06	0.383 0.65	0.24	0.24
	Lebanon	1990-2021 1955-2001	0.209 –0.045	0.251 0.765	0.327	0.068	0.05	0.81
Indices	NAOI	1990-2014					–0.29	0.17
	MOI	1990-2014					–0.08	0.7
Temperature	Greece	1990-2021 1951-2007	–0.19 0.27	0.31 0.043	–0.07 0.063	0.71 0.64	0.80	2*10–6
	Cyprus	1990-2021	–0.44	0.013	0.23	0.21	0.79	2.4*10–6
	Israel	1990-2021	–0.50	0.004	0.32	0.08	–0.72	5.2*10–5
	Turkey	1990-2021	–0.41	0.02	–0.26	0.15	0.75	1.6*10–5
	Lebanon	1990-2021	–0.49	0.005	0.25	0.17	–0.76	1.16*10–5

NAOI: North Atlantic Oscillation Index; MOI: Mediterranean Oscillation Index.
Bold characters denote significant trends at a 95% confidence level.

Selected correlation diagrams between NAOI and MOI and the precipitation series are presented in Figure 3.

Fig. 3. (a-e) Correlation diagrams between NAOI and MOI, and the precipitation series. (NAOI: North Atlantic Oscillation Index; MOI: Mediterranean Oscillation Index).

Based on the results of Table III and Figure 3, it is observed that during the period from 1990 to 2021, precipitation series form Greece, Cyprus, and Turkey have a negative correlation with the NAOI, whereas Israel and Lebanon have a positive correlation. None of the correlations are significant at the 95% c.l. The MOI has an increased positive correlation with Lebanon’s precipitation series (p = 0.068, non-significant at the 95% c.l.), a slightly increased negative correlation with Greece’s precipitation series (p = 0.138, non-significant at the 95% c.l.), and a slightly increased positive correlation with Israel’s precipitation series (p = 0.192, non-significant at 95% c. l.). As stated by Xoplaki et al (2004) and confirmed by the present study, the NAOI and MOI are not completely independent. For the examined period (1990-2021), the correlation coefficient between them is R = 0.635 (p = 0.00009, highly significant). The correlation between NAOI and Greece’s precipitation (R = –0.44, period 1955-2001) is in close agreement with the result (R = –0.41+/–1) of Feidas et al. (2007) for the same period, and in general qualitative agreement with the result of Philandras et al (2011) for the period between 1951 and 2010. The correlation between the NAOI and Turkey’s precipitation (R = –0.40, p = 0.00158, significant at the 95% c.l., period 1955-2013) differs from the result (R = –0.22, non-significant at the 95% c. l.) of Topuz et al. (2020) for the same period. The correlation referred to the MOI for the same period is also not confirmed by the results of Topuz et al. (2020).

According to Salameh et al. (2022), of the two indices MOI is correlated significantly (R = 0.27-0.38) between 1970 and 2020 with most of the precipitation series based on weather stations inf Israel, which is generally confirmed with the results of the present study (Table II, R = 0.27, p = 0.059) referred to the same period. However, the resulting correlation between MOI and Greece (R = –0.19, p = 0.23 for the period 1955-1995, non-significant at the 95% c.l.) differs gradually from the result (–0.4 ± 0.1, significant at the 95% c.l.) calculated by Feidas et al (2007) for the same period. It should be noted that both the results of Topuz et al. (2020) and Feidas et al. (2007) were referred to weather stations areal-averaged rainfall data.

Another interesting observation based on the last 30-yr period is the reduction of the influence (lower correlation coefficient) of the annual NAOI index (which changes from R = –0.40 to –0.278 from 1955-2013 to 1990- 2021, respectively, and from R = –0.44 to –0.295 from 1955-2001 to 1990-2021, respectively) in the Turkish and Greek territories. Taking the precipitation series of Greece as an example, this means that during the period from 1955 to 2001, 44% of the total rainfall could be explained by the relative NAOI, whereas this amount is reduced to 29.5% from 1990 to 2021.

The mean annual SST anomalies in the Mediterranean region have been studied to assess their simultaneous relationship with the EM ground-based annual precipitation. The correlation between mean annual SST anomalies and precipitation series (Table II) is generally positive but significant (at a 95% c.l.) only regarding the precipitation series of Greece (R = 0.50, p = 0.011). However, this single significant correlation seems to be temporary and is valid only for this examined period; however, no significant correlation is valid during the sub-periods 1960-2014, 1970-2014, and 1980-2014 (results are not presented in Table II). As a result, the Mediterranean annual SST anomalies did not govern the precipitation anomalies in the neighboring continental regions, except for the Greece territory, between 1990 and 2014.

3.3 Temperature trends

Table I also presents the trend of annual temperature series in the EM region, together with mean annual SST anomalies in the Mediterranean region. Results in Table II refer mainly to the period from 1990 to 2021 and other common periods for comparison. The annual variation of the series together with the relative linear trend (Sen’s slope) are visualized in Figure 4.

Fig. 4. (a-e) Variations of the mean annual temperature series, (f) mean SST series, together with the relative trend line (Sen’s slope). (SST: sea surface temperature).

Highly significant (at a 99% c.l.) annual temperature trends were detected for all the mean temperature series, including the SST series. These warming temperature trends span between 0.047 ºC yr^–1 (Greece, Israel, Lebanon) and 0.059 ºC yr^–1 (Cyprus, Turkey).

The results of the present study for annual temperature in various periods (Table II) are in close agreement with results found in the available literature (Tzanis et al., 2019; Yosef et al., 2019) and in general agreement with the results of Kadioğlu (1997), Price et al. (1999), and Shaban (2011). However, the area-averaged value of the trend calculated by Kadioğlu (1997) and Price et al. (1999) may not be representative of the Turkish (based on only 17 weather stations) and Cypriote (based on only two weather stations) territories.

Warming annual trends of the same order of magnitude as those found in the present study were also observed for the SST since 1985 by several studies (Nykjaer 2009; Shaltout and Omstedt, 2014; Androulidakis and Krestenitis, 2022). Selected correlation diagrams between the NAOI and the temperature series (including the SST) are presented in Figure 5.

Fig. 5. (a-e) Correlation diagrams between the NAOI and the temperature series, including the (f) SST. (NAOI: North Atlantic Oscillation Index; SST: sea surface temperature).

Significant correlation coefficients (greater than 95% c.l.), during 1990-2021 between the temperature series and the NAOI (except for the series of Greece) were also detected (Table III), in contrast with the relevant results of annual precipitation series. The observed significant correlation coefficients span between 0.41 (Turkey) and 0.50 (Israel). Results in Table III also indicate that the MOI did not play an important role in the annual temperature variation of the EM territory during the period from 1990 to 2021. Regarding the annual SST anomalies, highly significant correlations (greater than 95% c.l.) were observed for all the annual temperature series, ranging between 0.72 (Israel) and 0.80 (Greece). These findings show that between 1990 and 2014, both seawater and land responded synchronously to atmospheric forcing. However, an increasing negative correlation (non-significant at a 95% c.l.) is observed between SST anomalies and the NAOI, whereas almost no correlation is observed between SST and the MOI.

The correlation results of the present study (Table III, Fig. 5) between the temperature series of Greece and the MOI are in agreement with the findings of Nastos et al. (2011), which are mostly positive with non-statistical significance at a 95% c.l. However, in the case of the NAOI, the results of the present study do not confirm the non-significant correlation results of Nastos et al (2011), which significantly are referred to area-averaged temperature based on 29 weather stations within the Greek territory.

Finally, it is worth mentioning that an overall significant or non-significant yearly trend on a national basis does not mean that the regional trend remains the same within the national area. Precipitation trends are relative to several factors, such as the altitude of the region, orientation, distance from the sea, etc. Taking as an example the case of Greece, regional precipitation trends in recent years are mostly stationary, with a few of them positive/negative and significant (Margaritidis, 2021; Sakalis, 2022; Stathi et al., 2023). On the other hand, in the case of temperature, a statistically significant increase in air temperature in recent decades is more pronounced (Proutsos et al., 2020 Kastridis et al., 2022; Stathi et al., 2023). However, the investigation of regional trends is beyond the scope of the present study.

3.4 Forecast results

Several ARIMA and SARIMA models were tested in order to forecast the yearly and monthly precipitation and temperature series for the next four and eight years, respectively. Before the models’ identification, the stationary test of the time series was evaluated by the Kwiatolsky-Philips-Scmidt-Shin (KPSS) (Kwiatkowsky et al., 1992), the Dickey and Fuller (DF) (Dickey and Fuller, 1979), and the Phillips and Peron (PP) (Phillips and Peron, 1988) tests at a 95% c.l. When at least one of these three tests indicated the possible existence of a unit root in the time series, then a model with q = 1 or even 2 was entered in the model test list (differentiation of the time series). The initial models were identified by the inspection of the autocorrelation function (ACF) and partial autocorrelation function (PACF) diagrams. The parameter D of the SARIMA model, which is used for the monthly (seasonal) forecast, was set to D = 12 due to the obvious 12-month periodicity. All the possible combinations using several values of p, q and P, Q were examined in order to determine the best model. The model that gave the best combination of minimum RMSE (Eq. [12]) and MAE (Eq. [11]), the minimum value of the AIC criterion (Akaike, 1994), and the highest Nash-Sutcliffe efficiency coefficient (NSE; Eq. [13]) was selected as the best-fit model (Table III). At the end of the procedure, the identification of the shape of the ACF and PACF functions of the residuals was essential (Fig. 6). When the majority of the ACF and PACF functions lie within the confidence limits, this indicates that there is no significant correlation among them. Regarding yearly forecasts, the period 1990-2020 is used as training period and 2021-2028 as the forecast period. Considering that the yearly values of 2021 are already known, this year is used as validity check of the models. In the case of monthly (seasonal) forecasts (SARIMA models), the years 2021-2024 were used as forecast period and the monthly values of 2021 were used for validation of the seasonal models.

Fig. 6. Selected autocorrelograms and partial autocorrelograms of the models’ residual: (a) ARIMA (15,0,2) for Israel precipitation series, (b) SARIMA (2,0,7) (5,1,2) for Turkey precipitation series, (c) ARIMA (24,2,5) for Lebanon temperature series, and (d) SARIMA (1,0,0) (3,1,1) for Greece temperature series. (ARIMA: autoregressive integrated moving average; SARIMA: seasonal ARIMA).

3.4.1 Precipitation

According to Table IV, of the five selected ARIMA models (one for each country) referred to the yearly forecast of the training period, the MAE ranges from 12.93 to 65.13 mm yr^–1, with a mean value of 30.67 mm yr^–1, while the RMSE ranges from 17.76 to 76.52 mm yr^–1 (mean: 37.62 mm yr^–1) and the NSE from 0.76 to 0.919 mm yr^–1 (mean: 0.87 mm yr^–1).

Table IV. Efficiency of the ARIMA and SARIMA models.

	Country	(S)ARIMA	MAE 1980-2020	RMSE 1980-2020	NSE 1980-2020	SD* (2021)	MAE* (2021)	MAE*/SD*	SD (2021)	RMSE (2021)	RMSE/ SD	MAE (2021)	MAE/ SD	R (2021)
Precipitation	Greece	(2,0,6) (5,1,2)	23.57	30.6	0.386	104.2	40.7	0.39	49.0	29.66	0.50	24.51	0.61	0.84
	Cyprus	(2,0,6) (5,1,2)	16.36	26.51	0.648	110.5	129	1.17	44.4	20.39	0.32	14.20	0.41	0.916
	Israel	(2,0,1) (3,1,1)	10.336	17.68	0.66	76.8	8.8	0.11	30.8	13.22	0.25	7.66	0.43	0.895
	Turkey	(2,0,7) (5,1,2)	14.67	19.43	0.53	63.5	60.5	0.95	28.6	17.08	0.47	13.45	0.60	0.80
	Lebanon	(2,0,5) (7,1,2)	22.53	36.54	0.712	159.7	89.94	0.56	68.7	22.26	0.24	16.34	0.32	0.95
	Greece	(24,1,5)	22.59	29.10	0.919	104.2	192.6	1.85
	Cyprus	(15,2,2)	26.76	34.20	0.901	110.5	142.58	1.29
	Israel	(15,0,2)	25.93	30.48	0.834	76.8	143.66	1.87
	Turkey	(24,2,6)	12.93	17.76	0.919	63.5	2.611	0.04
	Lebanon	(16,0,2)	65.13	76.52	0.76	159.7	11.84	0.07
Temperature	Greece	(1,0,0) (3,1,1)	0.8212	1.070	0.976	0.59	0.1012	0.17	6.89	0.877	0.12	0.79	0.13	0.9916
	Cyprus	(6,0,3) (6,1,5)	0.826	1.075	0.972	0.72	0.0384	0.05	6.39	1.007	0.15	0.931	0.16	0.987
	Israel	(1,0,0) (3,1,2)	0.780	0.995	0.962	0.62	0.080	0.13	5.1	0.710	0.12	0.622	0.14	0.99
	Turkey	(1,0,0) (0,1,2)	0.969	1.263	0.976	0.76	0.054	0.07	8.11	0.962	0.10	0.837	0.12	0.993
	Lebanon	(1,0,0) (0,1,5)	0.849	1.092	0.971	0.68	0.0036	0.01	6.42	0.849	0.11	0.696	0.13	0.991
	Greece	(23,2,5)	0.166	0.225	0.85	0.59	0.327
	Cyprus	(23,2,4)	0.253	0.329	0.782	0.72	0.301
	Israel	(23,1,12)	0.127	0.165	0.928	0.62	0.04
	Turkey	(23,2,5)	0.284	0.37	0.76	0.76	0.452
	Lebanon	(24,2,5)	0.239	0.32	0.77	0.68	0.258

SARIMA: seasonal autoregressive integrated moving average; MAE: mean absolute error; RMSE: root mean square error; NSE: Nash-Sutcliffe efficiency coefficient; SD: standard deviation; R: Correlation coefficient.
MAE*, SD*, and MAE, SD refer to yearly and monthly values, respectively.
Forecasted results in bold type may have issues of acceptance.

The MAE between forecasted and measured values in 2021 ranges from 2.611 (Turkey) to 192.6 mm yr^–1 (Greece), with an average of 98.7 mm yr^–1. Despite the sufficient accuracy of the precipitation forecast referred to the training period, the validation of the models in 2021 indicates medium accuracy, owing to the high variance of the precipitation series. According to the NSE criterion (Eq. [13]), values between 0.0 and 1.0 are generally viewed as acceptable levels of performance, with a value of 1 being the optimal. Providing the division of MAE and RMSE with the standard deviation of the observed data (Table IV), stricter criteria may occur. According to this, values of MAE/SD and RMSE/SD less than 0.65 may be considered low and acceptable (Singh et al. 2005; Moriasi et al., 2007; Margiorou et al., 2022). Based on the criteria analyzed before, ARIMA forecasted results for Greece, Cyprus, and Israel may be recognized as non-acceptable (underlined values in Table IV). The same order in accuracy results (referred to the R and the RMSE parameters) is observed in other studies (Mahmud et al., 2017; Wang et al., 2021). The appropriateness of the yearly ARIMA models for forecasting precipitation is partly confirmed.

Regarding the SARIMA models, the MAE of the monthly forecast of the training period indicates a mean value of 22.01 mm yr^–1, the RMSE a value of 31.97 mm yr^–1 and the NSE is 0.59. The mean values of the MAE, RMSE and correlation coefficient (R) between the forecasted monthly values and the true existing value of 2021 are 15.23, 20.52, and 0.88 mm yr^–1, respectively. According to the criteria of MAE and RMSE divided by de SD of the observed data, forecasted results based on SARIMA for Turkey and Cyprus may not be acceptable (higher than 0.65) for yearly precipitation (bold values of Table IV). However, all the monthly forecasted results of the SARIMA model were sufficiently higher than the acceptance values. These results confirm the general appropriateness of the seasonal (SARIMA) models in monthly forecasting. It should be noted also that the accuracy of the forecasted results of the validation period (12 months of the year 2021) are slightly higher than in the training period. Another interesting result is the ability of the SARIMA model to forecast the yearly precipitation of 2021 (average MAE value: 65.79 mm yr^–1) better than the relevant ARIMA model (98.7 mm yr^–1). Finally, the overall performance of the seasonal SARIMA models seems to be higher than the non-seasonal ARIMA models.

Table V presents the values of the average yearly rainfall during the period from 2022 to 2024, based on the SARIMA model. It is observed that between the period 2022-2024 and the reference period 1990-2021, the average forecasted precipitation is increased by 76.53, 24.72, and 10.26 mm yr^–1 for Greece, Turkey, and Cyprus, respectively, and reduced by 22.16 and 2.54 mm yr^–1 for Israel and Lebanon, respectively.

Table V. Averaged yearly precipitation and temperature forecast values between 2022 and 2024 based on SARIMA models.

	Country	Average 1990-2021	Predicted average 2022-2024	Predicted average 2022	Predicted average 2023	Predicted average 2024
Precipitation (mm yr–1)	Greece	683.81	760.64	762.574	747.242	772.086
	Cyprus	456.39	466.65	494.32	435.152	470.47
	Israel	265.52	243.36	192.061	261.95	256.825
	Turkey	598.78	623.50	624.108	614.163	632.22
	Lebanon	652.16	649.32	704.092	610.53	633.336
Temperature (ºC)	Greece	14.04	14.73	14.70825	14.73217	14.73375
	Cyprus	18.86	19.68	19.66	19.71	19.67
	Israel	19.97	20.67	20.65	20.68	20.69
	Turkey	11.62	12.49	12.48	12.49	12.49
	Lebanon	15.37	16.11	16.08	16.13	16.11

SARIMA: seasonal autoregressive integrated moving average.

The observed and predicted rainfall data, based on ARIMA and SARIMA, are visualized in Figures 7 and 8, respectively.

Fig. 7. (a-e) Observed precipitation data, and ARIMA predicted and forecasted data for the period 1990-2028. (ARIMA: autoregressive integrated moving average).

Fig. 8. (a-e) Precipitation data: observed, predicted with SARIMA, forecasted, and real (2021) for the period 1990-2024. (SARIMA: seasonal autoregressive integrated moving average).

3.4.2 Temperature

The MAE, RMSE, and NSE of the five selected ARIMA models referred to the yearly forecast of the training period indicate a mean value of 0.214, 0.28, and 0.82 ºC, respectively. Correspondingly, the MAE between the forecasted and measured values of the year 2021 ranges from 0.04 and 0.452 ºC for Israel and Turkey, respectively, with an average of only 0.276 ºC (Table III).

Regarding the SARIMA models, the MAE and the RMSE of the monthly forecast of the training period indicate mean values of 0.849 and 1.099 ºC, respectively, and NSE of 0.972. The mean values of the MAE, RMSE and correlation coefficient (R) between the forecasted and measured monthly values of the year 2021 are 0.775 ºC, 0.881 ºC, and 0.991, respectively. As it is indicated in the case of the precipitation forecast, the accuracy of the forecasted mean temperature of 2021 by the seasonal SARIMA model (average MAE: 0.056 ºC) is better than the relevant ARIMA model (0.19 ºC), and the accuracy of the forecasted results of the validation period (12 months during 2021) is slightly higher than in the training period. Based on the NSE criteria, both the ARIMA and SARIMA models used for temperature forecasting are acceptable based on the criteria of MAE/SD and RMSE/SD (Nash and Sutcliffe., 1970; Singh et al., 2005; Moriasi et al., 2007; Margiorou et al., 2022). The overall performance of the seasonal SARIMA models seems to be higher than the non-seasonal ARIMA models.

It is obvious that the ARIMA and SARIMA models used for the prediction of the average temperature are more efficient and reliable than the relevant models used for precipitation, mostly due to the lower variance of temperature series as compared to precipitation series.

The plot of the observed and predicted annual temperature verified that ARIMA models (Fig. 9) and SARIMA models (Fig. 10) are generally appropriate for forecasting.

Fig. 9. (a-e) Temperature data: observed, predicted with ARIMA, and forecasted for the period 1990-2028. (ARIMA: autoregressive integrated moving average).

Fig. 10. (a-e) Temperature data: observed, predicted with SARIMA, forecasted, and real (2021) for the period 1990-2024. (SARIMA: seasonal autoregressive integrated moving average.)

Finally, based on Table V the average forecasted temperature between the periods 1990-2021 and 2022-2024 is increased by 0.69, 0.87, 0.7, 0.74, and 0.82 ºC for Greece, Turkey, Israel, Lebanon, and Cyprus, respectively.

4. Conclusions

Temperature and precipitation are of the most interesting weather variables affected by climate change. In the EM region under consideration, only positive trends were observed for precipitation series on a national basis during the period from 1990 to 2021. Between them, the positive trend of Greece is the only significant one at a 95% c.l. These increasing positive trends were observed mostly during the last two decades. During the period from 1990 to 2021, the annual NAOI trend weakened, resulting possibly in warmer and wetter conditions in the influenced areas. The increased positive precipitation trend in Greece and Turkey during this period could be partially explained by the relative negative trend of the NAOI, which led to warmer and wetter conditions. The reversed negative sign of the NAOI trend during the last 30 years is a remarkable observation, in spite of the strong positive trends of the periods from 1955 to 2021 and from 1955 to 2013. None of the examined countries has a significant correlation with the NAOI. The MOI has a strong positive and significant correlation (95% c.l.) only with the Lebanon precipitation series and a slight correlation with the Greece and Israel series. Although Greece and Turkey are sufficiently large countries in order to study their influences based on large-scale circulation patterns, the limited areas referred to Lebanon, Cyprus, and Israel may indicate a possible loss of consistency in response to circulation patterns.

The Mediterranean annual SST anomalies do not govern the precipitation anomalies in the neighboring continental regions, except for the Greek territory during the period from 1990 to 2014.

Highly significant (at a 99% c.l.) annual temperature trends were detected for all the available mean temperature series, including the SST series, during the period from 1990 to 2021. Significant correlation coefficients (greater than 95% c.l.) during 1990-2021 between the temperature series and the NAOI (except for the series of Greece) were also detected, in contrast with the relevant results of annual precipitation series. According to the observed significant correlation between temperature series and the SST, it is observed that the annual mean SST of the Mediterranean region is driven by the ground-based mean annual temperature of the EM and vice-versa. The MOI did not play an important role in the temperature variation of the EM territory during the period from 1990 to 2021 (correlation is not observed).

Despite the different types of weather data used in this and other relevant studies (gridded CRU, ERA5, ERA interim, NCEP/NCAR reanalysis, or ground-based weather station data), the results of the present study are mostly in agreement with the available results found in the literature.

The influence of climate change in the EM seems to be higher regarding mean annual temperature (including the SST) as compared to mean annual precipitation. However, several relatively extreme indices need to be considered in order to document a complete study.

Based on this study and compared to the results found in the available literature, the (S)ARIMA models are appropriate in general for forecasting future temperature and precipitation. Especially, forecasted results based on temperature series indicate higher efficiency (NSE, MAE/SD, RMSE/SD, and coefficient R) than the results of precipitation series. For all the series under consideration, the forecast results of the SARIMA model, referred both to precipitation and temperature in the validating year (2021), are more accurate than the relevant non-seasonal ARIMA models.

Despite the satisfactory accuracy of the predicted precipitation results of the non-seasonal ARIMA models referred to the training period, the models’ validation in the validating year (2021) is sometimes unacceptable, mostly due to the high variance of the precipitation series compared to the temperature series, and also due to the nature of the problems. The single year (2021) used for the validation of the non-seasonal ARIMA models could give an appropriate order of magnitude of the model’s accuracy. However, in the present study, it is imperative to wait for a longer validation period in order to reach more definite and safe conclusions. On the other hand, the validation period of the relevant SARIMA models (12 months) is supposed to be sufficient to extract safer conclusions referred to the departure from the real values.

Finally, it should be mentioned that cases in which an accurately validated ARIMA or SARIMA model does not lead to equivalent accurate relevant forecasted results, are often observed in the available literature. This means that such cases are representative of the limit of the ability of the ARIMA/SARIMA models to produce accurate and reliable forecasts. As a result of these possible issues, several advanced ARIMA models or even different methods (ANN, decomposition, spectral methods, etc.) are already being used by the scientific community to enhance the accuracy of the interesting but difficult and complicated problem of forecasting.

A possible improvement of the present study, which is the analysis of trend detection of temperature and precipitation on a seasonal basis, could be an interesting subject, especially because seasonal large-scale circulation patterns (winter and summer NAOI and MOI) have generally more pronounced climate forcing than annual large-scale circulation patterns. The patterns in winter and summer may depend on different factors, and if their annual average deviates from a long-term average it may not be related to any significant climate forcing and could be the result of pure coincidence.

Acknowledgments

The author wishes to thank the Northern Kynouria Municipal Water supply and Sewerage Company for the opportunity to conduct the relevant research study.

Special thanks to Climate Change Knowledge Portal CCKP, World Bank Group, for providing the ability to download the CRU TS 4.05 data, to the European Environmental Agency for the SST data, to the Climate Analysis Section of the National Center of Atmospheric Research (NCAR Climate Data Guide) for the NAOI data, and finally to the Climate Research Unit CRU, for the MOI data.

References

Aalijahan M, Karataş A, Lupo AR, Efe B, Khosravichenar A. 2023. Analyzing and modeling the spatial-temporal changes and the impact of GLOTI index on precipitation in the Marmara region of Türkiye. Atmosphere 14: 489. https://doi.org/10.3390/atmos14030489

Ajjur SB, al-Ghamdi SG. 2021. Variation in seasonal precipitation over Gaza (Palestine) and its sensitivity to teleconnection patterns. Water 13: 667. https://doi.org/10.3390/w13050667

Akaike H. 1974. A new look at the statistical model identification. IEEE Transactions on Automatic Control 19: 716-723. https://doi.org/10.1109/TAC.1974.1100705

Al-Balasmeh O, Babbar R, Karmaker T. 2019. Trend analysis and ARIMA modeling for forecasting precipitation pattern in Wadi Shueib catchment area in Jordan. Arabian Journal of Geosciences 12: 27. https://doi.org/10.1007/s12517-018-4205-z

Alexandersson H. 1986. A homogeneity test applied to precipitation data. Journal of Climatology 6: 661-675. https://doi.org/10.1002/joc.3370060607

Al-Rimmawi Η, Ghanem Μ, Shalash Ι. 2010. Rainfall trends in the district of Ramallah and al-Bireh, Palestine. Journal of Water Resource and Protection 2: 345-352. https://doi.org/10.4236/jwarp.2010.24040

Androulidakis, YS, Krestenitis, YN. 2022. Sea surface temperature variability and marine heat waves over the Aegean, Ionian, and Cretan seas from 2008-2021. Journal of Marine Science and Engineering 10: 42. https://doi.org/10.3390/jmse10010042

Ben-Gai T, Bitan A, Manes A, Alpert P, Rubin S. 1999. Temporal and spatial trends of temperature patterns in Israel. Theoretical and Applied Climatology 64: 163-177. https://doi.org/10.1007/s007040050120

Bethoux JP, Gentili B, Raunet J, Tailliez D. 1990. Warming trend in the western Mediterranean deep water. Nature 347: 660-662. https://doi.org/10.1038/347660a0

Box GEP, Jenkins GM. 1976. Time series analysis: Forecasting and control. Holden-Day, San Francisco, USA, 575 pp.

Buishand TA. 1982. Some methods for testing the homogeneity of rainfall records. Journal of Hydrology 58: 11-27. https://doi.org/10.1016/0022-1694(82)90066-X

Burn DH, Elnur MAH. 2002. Detection of hydrological trends and variability. Journal of Hydrology; 255: 107-122. https://doi.org/10.1016/S0022-1694(01)00514-5

Caloiero Τ, Caloiero P, Frustaci F. 2018. Long-term precipitation trend analysis in Europe and in the Mediterranean basin. Water and Environment Journal 32: 1-13. https://doi.org/10.1111/wej.12346

Cid N, Bonada N, Carlson SM, Grantham TE, Gasith A, Resh VH. 2017. Variability is a defining component of Mediterranean-climate rivers and their biota. Water 9: 52. https://doi.org/10.3390/w9010052

CCKP. 2021. Climate Change Knowledge Portal. World Bank Group. Available at https://climateknowledgeportal.worldbank.org/ (accessed 13 Nov 2022)

Coban V, Guler E, Kilic T, Kandemir SY. 2021. Precipitation forecasting in Marmara region of Turkey. Arab Journal of Geosciences 14, 86. https://doi.org/10.1007/s12517-020-06363-x

Conte M, Giuffrida S, Tedesco S. 1989. The Mediterranean oscillation. Impact on precipitation and hydrology in Italy. Proceedings of the Conference on Climate and Water, vol. 1. Academy of Finland 9: 121-137.

De Luis M, Brunetti M, González-Hidalgo JC, Longares LA, Martín-Vide J. 2010. Changes in seasonal precipitation in the Iberian Peninsula during 1946-2005. Global and Planetary Change 74: 27-33. https://doi.org/10.1016/j.gloplacha.2010.06.006

Dickey DA, Fuller WA. 1979. Distribution of the estimators for autoregressive time series with a unit root. Journal of the American Statistical Association 74: 427-431. https://doi.org/10.2307/2286348

EEA. 2012. European Environmental Agency. Maps and graphs. Available at https://www.eea.europa.eu/data-and-maps/figures/annual-average-sea-surface-temperature (accessed 17 Nov 2022).

Feidas H, Noulopoulou N, Makrogiannis T, Bora-Senta E. 2007. Trend analysis of precipitation time series in Greece and their relationship with circulation using surface and satellite data: 1955-2001. Theoretical and Applied Climatology 87: 155-177. https://doi.org/10.1007/s00704-006-0200-5

Feidas H. 2017. Trend analysis of air temperature time series in Greece and their relationship with circulation using surface and satellite data: recent trends and an update to 2013. Theoretical and Applied Climatology 129: 1383-1406. https://doi.org/10.1007/s00704-016-1854-2

Feng G, Cobb S, Abdo Z, Fisher DK, Ouyang Y, Adeli A, Jenkins JN. 2016. Trend analysis and forecast of precipitation, reference evapotranspiration, and rainfall deficit in the Blackland Prairie of eastern Mississippi. Journal of Applied Meteorology and Climatology 55: 1425-1439. https://doi.org/10.1175/JAMC-D-15-0265.1

Giannakopoulos C, Hadjinicolaou P, Kostopoulou E, Varotsos KV, Zerefos C. 2010. Precipitation and temperature regime over Cyprus as a result of global climate change. Advances in Geosciences 23: 17-24. https://doi.org/10.5194/adgeo-23-17-2010

Giorgi F. 2006. Climate change hot-spots. Geophysical Research Letters 33, 101029. https://doi.org/10.1029/2006GL025734

Givati A, Rosenfeld D. 2013. The Arctic Oscillation, climate change and the effects on precipitation in Israel. Atmospheric Research 132-133: 114-124. https://doi.org/10.1016/j.atmosres.2013.05.001

Gofa F, Mamara A, Anadranistakis M, Flocas H. 2019. Developing gridded climate data sets of precipitation for Greece based on homogenized time series. Climate 7, 68. https://doi.org/10.3390/cli7050068

Hadi SJ, Tombul M. 2017. Conversion of CruTS 3.23 data and evaluation of precipitation and temperature variables in a local scale. In: Proceedings of the International Conference on Advances in Sustainable Construction Materials and Civil Engineering Systems (ASCMCES-17), April 18-20, Sharjah, UAE.

Hadi SJ, Tombul M. 2018. Long-term spatiotemporal trend analysis of precipitation and temperature over Turkey. Meteorological Applications 25: 445-455. https://doi.org/10.1002/met.1712

Hafez Y. 2018. A recent study of seasonal and interannual climate variability over the Eastern Mediterranean region. Journal of Geoscience and Environment Protection 6: 132-151. https://doi.org/10.4236/gep.2018.61009

Hamed KH, Ramachandra Rao A. 1998. A modified Mann-Kendall trend test for autocorrelated data. Journal of Hydrology 204: 182-196. https://doi.org/10.1016/S0022-1694(97)00125-X

Harris I, Osborn TJ, Jones P, Lister D. 2020. Version 4 of the CRU TS monthly high-resolution gridded multivariate climate dataset. Scientific Data 7: 109. https://doi.org/10.1038/s41597-020-0453-3

Hurrell JW, Kushnir Y, Visbeck M, Ottersen G. 2003. An overview of the North Atlantic Oscillation. In: The North Atlantic Oscillation: Climate significance and environmental impact (Hurrell JW, Kushnir Y, Ottersen G, Visbeck M, Eds.). Geophysical Monograph Series 134: 1-35. https://doi.org/10.1029/134GM01

IPCC. 2007. Climate change 2007: The physical science basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change (Solomon S, Qin D, Manning M, Chen Z, Marquis M, Averyt KB, Tignor M, Miller HL, Eds.). Cambridge University Press, Cambridge, USA, 996 pp.

Kadioğlu M. 1997. Trends in surface air temperature data over Turkey. International Journal of Climatology 17: 511-520. https://doi.org/10.1002/(SICI)1097-0088(199704)17:5<511::AID-JOC130>3.0.CO;2-0

Kambezidis HD, Larissi IK, Nastos PT, Paliatsos AG. 2010. Spatial variability and trends of the rain intensity over Greece. Advances in Geosciences 26: 65-69. https://doi.org/10.5194/adgeo-26-65-2010

Kastridis A, Kamperidou V, Stathis D. 2022. Dendroclimatological analysis of fir (A. borisii-regis) in Greece in the frame of climate change investigation. Forests 13, 879. https://doi.org/10.3390/f13060879

Kendall MG. 1955. Rank correlation methods. Charles Griffin, London, UK.

Kostopoulou E, Jones P. 2005. Assessment of climate extremes in the Eastern Mediterranean. Meteorology and Atmospheric Physics 89: 69-85. https://doi.org/10.1007/s00703-005-0122-2

Kuglitsch FG, Toreti A, Xoplaki E, Della-Marta EPM, Zerefos CS, Türkeş M, Luterbacher MJ. 2010. Heat wave changes in the eastern Mediterranean since 1960. Geophysical Research Letters 37): 1-5. https://doi.org/10.1029/2009GL041841

Kwiatkowsky D, Phillips PC, Scmidt P, Shin Y. 1992. Testing the null hypothesis of stationarity against the alternative of a unit root: How sure are we that economic time series have a unit root? Journal of Econometrics 54: 159-178. https://doi.org/10.1016/0304-4076(92)90104-Y

Lelieveld J, Hadjinicolaou P, Kostopoulou E, Chenoweth J, El Maayar M, Giannakopoulos C, Hannides C, Lange MA, Tanarhte M, Tyrlis E, Xoplaki E. 2012. Climate change and impacts in the Eastern Mediterranean and the Middle East. Climatic Change 114: 667-687. https://doi.org/10.1007/s10584-012-0418-4

Lu J, Vecchi GA, Reichler T. 2007. Expansion of the Hadley cell under global warming. Geophysical Research Letters 34: L06805. https://doi.org/10.4491/eer.2016.075

Mahmud I, Bari SH, Tauhid U, Rahman Μ. 2017. Monthly rainfall forecast of Bangladesh using autoregressive integrated moving average method. Environmental Engineering Research 22: 162-168. https://doi.org/10.4491/eer.2016.075

Mamara A, Argiriou AA, Anadranistakis M. 2015. Recent trend analysis of mean air temperature in Greece based on homogenized data. Theoretical and Applied Climatology 126: 543-573. https://doi.org/10.1007/s00704-015-1592-x

Mann HB. 1945. Nonparametric tests against trend. Econometrica 13: 245-259. http://dx.doi.org/10.2307/1907187

Margaritidis AK. 2021. Site and regional trend analysis of precipitation in central Macedonia, Greece. Computational Water, Energy, and Environmental Engineering 10: 49-70. https://doi.org/10.4236/cweee.2021.102004

Margiorou S, Kastridis A, Sapountzis M. 2022. Pre/post-fire soil erosion and evaluation of check-dams effectiveness in Mediterranean suburban catchments based on field measurements and modeling. Land 11: 1705. https://doi.org/10.3390/land11101705

Mersin D, Tayfur G, Vaheddoost B, Safari MJS. 2022. Historical trends associated with annual temperature and precipitation in Aegean Turkey. Where are we heading? Sustainability 14: 13380. https://doi.org/10.3390/su142013380

Mathbout S, López-Bustins JA, Roye D, Martín-Vide J, Bech J, Rodrigo FS. 2017. Observed changes in daily precipitation extremes at annual timescale over the Eastern Mediterranean during 1961-2012. Pure and Applied Geophysics 175: 3875-3890. https://doi.org/10.1007/s00024-017-1695-7

Moriasi DN, Arnold JG, van Liew MW, Bingner RL, Harmel RD, Veith TL. 2007. Model evaluation guidelines for systematic quantification of accuracy in watershed simulations. Transactions of the ASABE 50: 885-900. https://doi.org/10.13031/2013.23153

Nash JE, Sutcliffe JV. 1970. River flow forecasting through conceptual models. Part I-A. Discussion of principles. Journal of Hydrology 10: 282-290. https://doi.org/10.1016/0022-1694(70)90255-6

Nastos PT, Philandras CM, Founda D, Zerefos CS. 2011. Air temperature trends related to changes in atmospheric circulation in the wider area of Greece. International Journal of Remote Sensing 32: 737-750. https://doi.org/10.1080/01431161.2010.517796

Nykjaer L. 2009. Mediterranean Sea surface warming 1985-2006. Climate Research 39: 11-17. https://doi.org/10.3354/cr00794

Palutikof JP, Trigo RM, Adcock ST. 1996. Scenarios of future rainfall over the Mediterranean: Is the region drying? Proceedings of the conference on Mediterranean Desertification, Research Results and Policy Implications. Crete, Greece. 20 October-1 November.

Partal T, Kahya E. 2006. Trend analysis in Turkish precipitation data. Hydrological Processes 20: 2011-2006. https://doi.org/10.1002/hyp.5993

Patakamuri SK, Muthiah K, Sridhar V. 2020. Long-term homogeneity, trend, and change-point analysis of rainfall in the arid district of Ananthapuramu, Andhra Pradesh state, India. Water 12: 211. https://doi.org/10.3390/w12010211

Pettitt AN. 1979. A non-parametric approach to the change-point detection. Applied Statistics 28: 126-135. https://doi.org/10.1016/0022-1694(82)90066-X

Philandras CM, Nastos PT, Repapis CC. 2008. Air temperature variability and trends over Greece. Global NEST Journal 10: 273-285. https://doi.org/10.30955/gnj.000470

Philandras CM, Nastos PT, Kapsomenakis J, Douvis KC, Tselioudis G, Zerefos CS. 2011. Long term precipitation trends and variability within the Mediterranean region. Natural Hazards and Earth System Sciences 11: 3235-3250. https://doi.org/10.5194/nhess-11-3235-2011

Phillips PC, Perron P. 1988. Testing for a unit root in time series regression. Biometrika 75: 335-346. https://doi.org/10.1093/biomet/75.2.335

Proutsos N, Korakaki E, Bourletsikas A, Solomou A, Avramidou EV, Georgiadis C, Kontogianni AB, Tsagari K. 2020. Urban temperature trends in East Mediterranean: The case of Heraklion Crete. European Water 69/70: 3-14.

Piras Μ, Mascaro G, Deidda Ρ, Vivoni ΕR. 2016. Impacts of climate change on precipitation and discharge extremes through the use of statistical downscaling approaches in a Mediterranean basin. Science of The Total Environment 543: 952-964. https://doi.org/10.1016/j.scitotenv.2015.06.088

Pnevmatikos JD, Katsoulis BD. 2006. The changing rainfall regime in Greece and its impact on climatological means. Meteorological Applications 13: 331-345. https://doi.org/10.1017/S1350482706002350

Price C, Michaelides S, Pashiardis S, Alpert P. 1999. Long term changes in diurnal temperature range in Cyprus. Atmospheric Research 51: 85-98. https://doi.org/10.1016/S0169-8095(99)00022-8

Rahman MA, Yunsheng L, Sultana N. 2017. Analysis and prediction of rainfall trends over Bangladesh using Mann-Kendall, Spearman’s rho tests and ARIMA model. Meteorology and Atmospheric Physics 129: 409-424. https://doi.org/10.1007/s00703-016-0479-4

Ramadan HH, Beighley RE, Ramamurthy AS. 2013. Temperature and precipitation trends in Lebanon’s largest river: The Litani basin. Journal of Water Resources Planning and Management 139: 86-95. https://doi.org/10.1061/(ASCE)WR.1943-5452.0000238

Rogers JC, van Loon H. 1979. The seesaw in winter temperatures between Greenland and northern Europe. Part II: Some oceanic and atmospheric effects in middle and high latitudes. Monthly Weather Review 107: 509-519. https://doi.org/10.1175/1520-0493(1979)107<0509:TSIWTB>2.0.CO;2

Sahsamanoglou HS, Makrogiannis TJ. 1992. Temperature trends over the Mediterranean region. 1950-1988. Theoretical and Applied Climatology 45: 183-192. https://doi.org/10.1007/BF00866191

Sakalis, VD. 2022. Trend analysis of annual and seasonal precipitation data in Arcadia region (Greece). Bulletin of Atmospheric Science and Technology 3: 8. https://doi.org/10.1007/s42865-022-00052-1

Salameh AAM, Ojeda MG-V, Esteban-Parra MJ, Castro-Díez Y, Gámiz-Fortis SR. 2022. Extreme rainfall indices in southern Levant and related large-scale atmospheric circulation patterns: A spatial and temporal analysis. Water 14: 3799. https://doi.org/10.3390/w14233799

Schonwiese CD, Rapp J. 1997. Climate trend atlas of Europe-based on observations 1891-1990. Kluwer, Dordrecht. https://doi.org/10.1007/978-94-015-8818-8

Sen PK. 1968. Estimates of the regression coefficient based on Kendall’s tau. Journal of the American Statistical Association 63: 1379-1389. https://doi.org/10.1080/01621459.1968.10480934

Shaban A. 2011. Analyzing climatic and hydrologic trends in Lebanon. Journal of Environmental Science and Engineering 5: 483-492. https://doi.org/10.17265/1934-8932/2011.04.014

Shadeed S. 2013. Climate changes and trends in rainfall and temperature of Nablus meteorological station. Third International Conference on Energy and Environmental Protection in Sustainable Development (ICEEP III), October 9-10. Palestine Polytechnic University, Hebron, West bank, State of Palestine.

Shaltout M, Omstedt A. 2014. Recent sea surface temperature trends and future scenarios for the Mediterranean Sea. Oceanologia 56: 411-443. https://doi.org/10.5697/oc.56-3.411

Singh J, Vernon Knapp H, Arnold JG, Demissie M. 2005. Hydrological modeling of the Iroquois River watershed using HSPF and SWAT. Journal of the American Water Resources Association 41: 343-360. https://doi.org/10.1111/j.1752-1688.2005.tb03740.x

Soltani S, Modarres R, Eslamian S. 2006. The use of time series modeling for the determination of rainfall climates of Iran. International Journal of Climatology 27: 819-829. https://doi.org/10.1002/joc.1427

Stathi E, Kastridis A, Myronidis D. 2023. Analysis of hydrometeorological trends and drought severity in water demanding Mediterranean islands under climate change conditions. Climate 11: 106. https://doi.org/10.3390/cli11050106

Tadesse KB, and Dinka MO. 2022. The SARIMA model-based monthly rainfall forecasting for the Turksvygbult station at the Magoebaskloof dam in South Africa. Journal of Water and Land Development 53: 100-107. https://doi.org/10.24425/jwld.2022.140785

Todaro V, D’Oria M, Secci D, Zanini A, Tanda MG. 2022. Climate change over the Mediterranean region: Local temperature and precipitation variations at five pilot sites. Water 14: 2499. https://doi.org/10.3390/w14162499

Topuz Μ, Feidas Η, Karabulut Μ. 2020. Trend analysis of precipitation data in Turkey and relations to atmospheric circulation (1955-2013). Italian Journal of Agrometeorology (2): 91-107. https://doi.org/10.13128/ijam-887

Türkeş M. 1999. Vulnerability of Turkey to desertification with respect to precipitation and aridity conditions. Turkish Journal of Engineering and Environmental Science 23: 363-380.

Türkeş M, Sümer UM, Demi̇r I. 2002. Re-evaluation of trends and changes in mean, maximum and minimum temperatures of Turkey for the period 1929-1999. International Journal of Climatology 22: 947-977. https://doi.org/10.1002/joc.777

Türkeş M, Erlat E. 2003. Precipitation changes and variability in Turkey linked to the North Atlantic Oscillation during the period 1930-2000. International Journal of Climatology 23: 1771-1796. https://doi.org/10.1002/joc.962

Tzanis CG, Koutsogiannis I, Philippopoulosa K, Deligiorgi D. 2019. Recent climate trends over Greece. Atmospheric Research 230: 104623. https://doi.org/10.1016/j.atmosres.2019.104623

Van Loon H, Rogers JC. 1978. The seesaw in winter temperatures between Greenland and northern Europe. Part I: Winter. Monthly Weather Review 104: 365-380. https://doi.org/10.1175/1520-0493(1978)106<0296:TSIWTB>2.0.CO;2

Varlas G, Stefanidis K, Papaioannou G, Panagopoulos Y, Pytharoulis I, Katsafados P, Papadopoulos A, Dimitriou E. 2022. Unraveling orecipitation trends in Greece since 1950s using ERA5 climate reanalysis data. Climate 10: 12. https://doi.org/10.3390/cli10020012

Von Neumann J. 1941. Distribution of the ratio of the mean square successive difference to the variance. Annals of Mathematical Statistics 13: 367-395. https://doi.org/10.1214/aoms/1177731677

Walker GT. 1924. Correlations in seasonal variations of weather. IX. Further study of world weather. Memoirs of the India Meteorological Department 24: 275-346.

Walker GT, Bliss E. 1932. World weather V. Memoirs of the Royal Meteorological Society 4: 53-84.

Wang HR, Wang C, Lin X, Kang J. 2014. An improved ARIMA model for precipitation simulations. Nonlinear Processes in Geophysics 21: 1159-1168. https://doi.org/10.5194/npg-21-1159-2014

Wang H, Wang W, Du Y, Xu D. 2021. Examining the applicability of wavelet packet decomposition on different forecasting models in annual rainfall prediction. Water 13: 1997. https://doi.org/10.3390/w13151997

Wijngaard JB, Klein Tank AMG, Können GP. 2003. Homogeneity of 20th century European daily temperature and precipitation series. International Journal of Climatology 23: 679-692. https://doi.org/10.1002/joc.906.

Xoplaki E, González-Rouco JF, Luterbacher J, Wanner H. 2004. Wet season Mediterranean precipitation variability: Influence of large-scale dynamics and trends. Climate Dynamics 23: 63-78. https://doi.org/10.1007/s00382-004-0422-0

Yosef Y, Aguilar E, Alpert P. 2019. Changes in extreme temperature and precipitation indices: Using an innovative daily homogenized database in Israel. International Journal of Climatology 39: 5022-5045. https://doi.org/10.1002/joc.6125

Zaifoğlu H, Akıntuğ B, Yamaz AM. 2017. Quality control, homogeneity analysis, and trends of extreme precipitation indices in Northern Cyprus. Journal of Hydrologic Engineering 22: 05017024. https://doi.org/10.1061/(ASCE)HE.1943-5584.0001589

Zhang X, Aguilar E, Sensoy S, Melkonyan H, Tagiyeva U, Ahmed N, Kutaladze N, Rahimzadeh F, Taghipour A, Hantosh TH, Albert P, Semawi M, Karam Ali M, Al-Shabibi MHS, Al-Oulan Z, Zatari T, Khelet IAD, Hamoud S, Sagir R, Demircan M, Eken M, Adiguzel M, Alexander L, Peterson TC, Wallis T. 2005. Trends in Middle East climate extreme indices from 1950 to 2003. Journal of Geophysical Research 110: D22104. https://doi.org/10.1029/2005JD006181

Zittis G, Bruggeman A, Lelieveld J. 2021. Revisiting future extreme precipitation trends in the Mediterranean. Weather and Climate Extremes 34: 100380. https://doi.org/10.1016/j.wace.2021.100380

Zittis G, Almazroui M, Alpert P, Ciais P, Cramer W, Dahdal Y, Fnais M, Francis D, Hadjinicolaou P, Howari F, Jrrar A, Kaskaoutis DG, Kulmala M, Lazoglou G, Mihalopoulos N, Lin X, Rudich Y, Sciare J, Stenchikov G, Xoplaki E, Lelieveld J. 2022. Climate change and weather extremes in the Eastern Mediterranean and Middle East. Reviews of Geophysics 60: e2021RG000762. https://doi.org/10.1029/2021RG000762