A Recent Study of the Relationship between the Precipitation Rate over Saudi Arabia in the Fall and Climatic Indices

,is paper presents a recent study of the relationship between precipitation rate (PR) over Saudi Arabia (SA) within the months of the fall season and climatic indices. ,e fall monthly PR data spanning the study period between 1948 and 2018 is considered. In addition, the monthly climatic index records (arctic oscillation (AO), global surface air temperature (GSAT), multivariate ENSO index (MEI), North Atlantic Oscillation (NAO) index, Nino 3.4 index, and Southern Oscillation Index (SOI)) for the fall months were also considered. ,e statistical trend, anomaly, and correlation analyses are applied in this study. ,e results reveal that the sweeping changes in PR show generally positive trends throughout the fall seasons of the past decades. Moreover, the climatic indices have an effect on the PR over SA within the fall months and season. During the study period, the most substantial relationship recorded, with an inverse correlation of −0.7, is between the PR over SA and the climatic index of GSATfor September and October. Moreover, there is a clear correlation of +0.5 between the PR over SA and the ENSO and Nino 3.4 index for October and November.


Introduction
Recently, Saudi Arabia (SA) has suffered from an extreme precipitation variation, mainly during the fall season. However, the SA climatic zone is that of a dry climate. e distribution, quantity, and duration of precipitation in arid and semiarid zones were considered to be highly variable [1][2][3][4][5][6]. e northern regions of SA, the region around the Central North, and the southwest areas of SA recorded the highest level of precipitation (171.5 mm) during the year of 1982.
Meanwhile, the southeast areas of SA are arid. e southwest region of SA experiences precipitation throughout the twelve months of each year [7]. In fact, the SA climate varies significantly, according to the geography and climate yeartide of the annual cycle, and is influenced by circulations in the atmosphere [4,6,[8][9][10]. Actually, the regional circulation pattern of precipitation for the duration of the fall over SA was precisely related to the same regional circulation and distribution during the spring season [11].
Several global scientific works presented in the literature consider the teleconnections between precipitation and climatic indices. In particular, the authors in [12] used the Southern Oscillation Index (SOI) to enhance precipitation forecast and manage agricultural water in the Republic of Ghana. e authors of [13] studied how to predict rainfall over Australia by the use of the SOI phase. Recently, there are several scientific literatures challenge the teleconnections between the precipitations and large-scale phenomena like as atmospheric blocking systems [14][15][16].
From previous research works in [17][18][19][20], it was found that there is robust proof for the essence of the link between climatic conditions and the variability of the Arctic oscillation (AO). It appears that there is a connection between AO and the climatic conditions for many areas in the Northern Hemisphere [21][22][23][24][25][26][27][28]. Besides, weather variations and climatic conditions and North Atlantic Oscillation (NAO) variations are significantly and positively correlated [29][30][31]. Both AO and NAO variations play significant roles in understanding the physical and dynamical mechanisms of the variations of the Northern Hemisphere [32]. e authors of [33] realized that the seasonal and interannual changes in the average air temperature of SA show a warming period that began in the last years of the twentieth century. is temperature variability follows global warming signs.
Furthermore, it appears that the sturdy variability in temperature over SA was mostly related to the North Atlantic Oscillation indices, majorly throughout the winter season. Mostly, El Nino Southern Oscillation (ENSO) index, AO, and NAO indices take significant roles in the continuous temperature rise over SA [4,6,7,34]. Additionally, global mean surface temperature (GSAT) change was analyzed by many organizations, including the NASA Goddard Institute for Space Studies. ese studies found that global warming is enduring persistently. e authors of [31] noted that there is a significant, substantial connection between the global mean surface air temperature and precipitation over Europe. In the same manner, the authors of [35] found that the ENSO phenomenon generates a significant proportion of short-term climate variability globally. e relationship between ENSO and precipitation in South East Queensland, Australia, has been studied. e authors of [36] found that the Southern Oscillation Index (SOI) and rainfall amounts could be predicted in the short term using the analysis of time series methods. e authors of [6] found that SA climate parameters, temperature, and PR were correlated with and controlled by the Oceanic Nino Index (ONI). e present work aims to uncover the relationship between the PR variability over SA in the fall months and season and the climatic indices throughout the study period (1948-2018).

Materials and Methods
2.1. Data. National Center for Atmospheric Research (NCAR) and National Centers for Environmental Prediction (NCEP) reanalysis program makes use of a perfect analysis system to assimilate records with a resolution n of 2.5 × 2.5 degree lat/long grid, by making use of previous data between 1948 and 2018. e NCEP/NCAR reanalysis 1 is the source of Saudi Arabian PR used in the present work. However, the NCEP/NCAR reanalysis 1 project is using a state-of-the-art analysis/forecast system to perform data assimilation using past data from 1948 to the present time. A large subset of this data is available from the Physical Sciences Division (PSD) in its original four times daily format and as daily averages.
is data was obtained eight times in a day as per model, since the inputs available in that era were available at 3Z, 9Z, 15Z, and 21Z, whereas the four times daily data has been available as 0Z, 6Z, 12Z, and 18Z. ese latter times were forecasted and the combined result of this early era is eight times daily. e NCEP/NCAR reanalysis 1 data was obtained from https://www.esrl.noaa.gov/psd/data/gridded/data. ncep.reanalysis.html. e monthly statistics set of the PR over the region of SA for the study period of 1948 to 2018 for the fall months (September-October-November) was considered. Moreover, the monthly datasets for the climatic indices (AO, GSAT, MEI, NAO, Nino 3.4, and SOI) for the fall months of the study period (1948-2018) are used. is data was provided with the support of the NOAA, ESRL PSD, OAR, Boulder, Colorado, USA, and [37]. In the present work, the domain of the NCEP data for the PR over the Arabian Peninsula, including SA, was considered to extend over the area of 12°N-32°N, 33°E-60°E. Figure 1 shows a map, including the Arabian Peninsula and Saudi Arabia.

Methodology.
Time series and anomaly analysis methods were used to analyze the monthly and seasonal PR over SA throughout the study period . e climatic average values of the PR were considered throughout the time between 1981 and 2010. However, this data used for the time between 1981 and 2010 is the recent climatic mean taken from climate-related research works. e anomaly of precipitation (A ' ) was calculated as the difference between the precipitation amount (A) and its climatic mean (Â) for each grid point in the domain of the area of study during that time . Additionally, linear correlation methods were employed to study the teleconnections between the PR over SA in the fall season and the climatic indices through the period of 1948-2018 used for the current evaluations. e methodology of Monte Carlo has been used [38]. Statistically, based on this methodology, a correlation of 0.2 would be significant for 66 degrees of freedom. In addition to that, the linear correlation method according to [38,39] is used. e statistical significance was determined using the Kendall-tau test. According to that test, a trend was deemed to be "statistically significant" if it has at least 90% significance.

Computation of the Fall Seasonal Means.
A fall seasonal mean of PR for every grid point of the domain of SA (9 × 12-degree lat/long grids) was calculated, whereas 9 grid points from 12°north to 32°north for latitudes and 12 grid points from 33°east to 60°east for longitudes have been used. Fall PR seasonal means over SA were calculated. Both were computed by the use of NCEP/NCAR reanalysis monthly figures for the fall season ((September + October + November)/ three). e climatic PR means for the period of study (1981-2010) were also calculated. e monthly datasets of the climatic indices (AO, GSAT, MEI, NAO, Nino 3.4, and SOI) over the fall months of the study period  were analyzed in the similar method. e interactive plotting and analysis NCEP and NCAR software packages were applied in these computations.

Local Significance and Computation of Correlation
Patterns.
e set data for monthly and seasonal PR averages at every grid point correlated with the earlier climatic indices. t-tests were performed on the correlation of every grid point for local significance using [40], allowing sequential autocorrelation as per [38] technique. To evaluate the correlation, we apply the Monte Carlo methodology. e field significance statistic is deemed to be area-weighted average of the absolute correlation of a given correlation map. e field-significance threshold equals the 95th percentile of a 1000-member Monte Carlo population.

Variability of PR over SA in the Fall Season throughout the Period of 1948-2018.
e NCEP/NCAR reanalysis shows that all datasets of the fall PR over SA for between the years 1948 and 2018 were scrutinized using the anomaly technique. e results showed the following: (1) Time series analysis of the September monthly PR over SA shows that the PR is above its normal values from 1950 to 1966. e maximum positive anomaly is 1.1 mm/day in 1960. After that period, the amounts of PR oscillate around its normal value until 2008. For the period of 2009-2018, the PR remains above its normal values. e polynomial trend analysis of the September PR shows that it increased over the past 10 years with a positive trend, as illustrated in Figure 2.   Advances in Meteorology the north-central regions of SA during the study period, as shown in Figure 6(e) and Table 1.

Discussion.
e present work provides a teleconnection study between PR over SA in the fall season and the climatic indices. According to [8], there are extreme changes in temperature and precipitation in the Arab region. Some significant relationships were found between climate extremes with ENSO over the eastern part of the Arab region. e authors of [4] studied the rainfall features and variations over SA. It was found that the climate of SA was instructed by ENSO, NAO, and SOI. e authors of [6] studied the relationship between the Oceanic Nino Index (ONI) and temperature and precipitation rates over SA. It was found that the SA climatic parameters, temperature, and precipitation rates were controlled by ONI. From the present analysis of PR variability of the fall through the period of study , it is clear that the PR over SA varies dramatically from month to month throughout the fall season during the study period. It is obvious that each month (September, October, and November) has unique PR variations. Generally speaking, the trend of PR variations is a positive trend through the last decade. rough the present work, the teleconnection between the PR over SA and the climatic indices (AO, GSAT, MEI, NAO, Nino 3.4, and SOI) has been studied. It becomes clear that the PR over SA in fall season was affected significantly by the global warming and

Conclusions
e teleconnection between the precipitation and climatic indices has attracted the attention of researchers of climate and climatic changes for decades. is paper investigated a teleconnection study between PR over SA in the fall season and the climatic indices (AO, GSAT, MEI, NAO, Nino 3.4, and SOI) throughout the period of study . e results revealed that PR over SA has an outstanding significant correlation with the climatic indices. Finally, one can conclude that the climatic indices affect and control the PR over several regions of SA in fall season. However, the results obtained through the present work aid the future research regarding rainfall forecast in SA by using the assessed indices for prediction models.

Data Availability
NCEP/NCAR reanalysis datasets of monthly precipitation rate are available at https://psl.noaa.gov/data/composites/ datasets.html. A monthly/seasonal meantime series of the precipitation rate from the NCEP reanalysis dataset is available at https://psl.noaa.gov/cgi-bin/data/timeseries/ timeseries1.pl. e datasets of monthly climate indices (AO, GSAT, MEI, NAO, Nino 3.4, and SOI) are available at https://psl.noaa.gov/data/climateindices/. e datasets of linear correlations in Atmospheric Seasonal/Monthly Averages of climatic indices with weather elements are available at NOAA physical science laboratory.

Conflicts of Interest
e authors declare that they have no conflicts of interest.