Indonesian Journal of Forestry Research Vol. 8 No. April 2021, 23-35 ISSN: 2355-7079/E-ISSN: 2406-8195 RAINFALL VARIABILITY IN GUNUNG SEWU KARST AREA, JAVA ISLAND. INDONESIA Ahmad Cahyadi1*. Eko Haryono2. Tjahyo Nugroho Adji2. Margaretha Widyastuti2. Indra Agus Riyanto3. Dzakwan Taufiq Nur Muhammad1, and Naufal Fattah Tastian1 Department of Environmental Geography. Universitas Gadjah Mada. Sekip Utara Jalan Kaliurang. Bulaksumur. Yogyakarta, 55281 Indonesia Karst Research Group. Universitas Gadjah Mada. Sekip Utara. Jalan Kaliurang Bulaksumur. Yogyakarta. 55281 Indonesia Lentera Geosains Yogyakarta. Jl. Tongkol V. Minomartani. Ngaglik. Sleman. Yogyakarta Indonesia Received: 21 October 2020. Revised: 5 April 2021. Accepted: 6 April 2021 RAINFALL VARIABILITY IN GUNUNG SEWU KARST AREA. JAVA ISLAND. INDONESIA. Karst area is highly susceptible to changes to climate parameters. One of the parameters is rainfall variability. addition to shaping the condition of water resources, rainfall in the Gunung Sewu karst area determines the nature of crop and livestock of the agriculture sectorsAithe local population's main economic activities, warranting the significance of the rainfall variability studies. Rainfall variability in karst areas also affects disaster conditions such as drought and floods. However, due to insufficient meteorological data in quality and quantity, there has been no rainfall variability studies conducted in this locality. The research intended to analyze rainfall variability in the Gunung Sewu karst area in 1979Ae2013 by utilizing rainfall predictions from satellite images that many scholars had tested in different locations and recognized as having good quality. In the analysis, mean monthly rainfall was calculated, and the trends of annual rainfall and average rainfall intensity, dry and rainy seasons, the number of rainy days, and the effect of ENSO (El Niyo Southern Oscillatio. on rainfall were analyzed. The research data were 35 years of daily rainfall records derived from the National Centers for Environmental Prediction (NCEP) Climate Forecast System Reanalysis (CFSR). The analysis results showed that the mean rainfall, number of rainy days, and rainfall intensity had an increasing trend. Also. El Niyo quantitatively influenced the rain in the Gunung Kidul karst area. Keywords: Rainfall variability, karst. ENSO. Gunung Sewu. NCEP. CFSR VARIABILITAS CURAH HUJAN DI KAWASAN KARST GUNUNG SEWU. PULAU JAWA, INDONESIA. Kawasan karst merupakan kawasan yang sangat sensitif terhadap perubahan parameter iklim, salah satunya adalah pada variabilitas curah hujan. Variabilitas curah hujan menjadi sangat penting dikaji di kawasan karst Gunung Sewu karena selain sangat mempengaruhi kondisi sumber daya air, curah hujan akan berpengaruh terhadap sektor pertanian dan peternakan yang menjadi sektor ekonomi utama bagi masyarakat di kawasan karst Gunung Sewu. Variabilitas curah hujan di kawasan karst juga mempengaruhi kondisi bencana seperti kekeringan dan banjir. Namun demikian, analisis variabilitas di kawasan karst Gunung Sewu belum dilakukan karena ketersediaan data meteorologis yang kurang baik dari sisi kualitas ataupun kuantitas. Penelitian ini bertujuan untuk menganalisis variabilitas curah hujan di kawasan karst Gunung Sewu pada tahun 1979Ae2013 dengan memanfaatkan data prediksi dari citra satelit yang telah diuji oleh banyak peneliti di lokasi lain dan memiliki kualitas yang baik. Analisis variabilitas hujan yang dilakukan dalam penelitian ini meliputi perhitungan rerata hujan bulanan, analisis trend hujan tahunan, analisis hujan musim kemarau dan penghujan, analisis jumlah hari hujan, analisis trend rata-rata intensitas hujan dan analisis pengaruh ENSO (El Niyo Southern Oscillatio. terhadap curah hujan di kawasan karst Gunung Sewu. Data yang digunakan dalam penelitian ini berasal dari The National Centers for Environmental Prediction (NCEP) Climate Forecast System Reanalysis (CFSR). Data yang digunakan berupa data harian selama 35 tahun. Hasil analisis menunjukkan bahwa rerata curah hujan, jumlah hari hujan dan intensitas hujan memiliki trend kenaikan. Selain itu diketahui bawah hujan di kawasan karst Gunung Kidul secara kualitatif nampak sangat dipengaruhi oleh peristiwa EL Niyo. Kata kunci: Variabilitas hujan, karst. ENSO. Gunung Sewu. NCEP. CFSR Corresponding author: ahmad. cahyadi@geo. A2021 IJFR All rights reserved. Open access under CC BY-NC-SA license. doi:10. 20886/ijfr. Indonesian Journal of Forestry Research Vol. No. April 2021, 23-35 INTRODUCTION Karst landscapes are fragile and highly sensitive to environmental changes in which ecological restoration can be challenging to achieve quickly (Quine et al. , 2017. Kang et al. , 2. The changes in question are those of hydrometeorological or climatic variables, which currently are prevalent issues. Dissolution as the dominant process in karst areas has led to a hydrologically dry surface because all water landing on this landscape is drained immediately into underground rivers (Cahyadi, 2014. Here, water storage highly depends on the epikarst layer, i. , the top layer of a karst landscape that consists of soil and the widening zone of fissures and conduits in carbonate rocks. In other words, soil thickness and the presence of loose rock grains strongly determine its capacity. However, resistance to rock disintegration has limited the formation of soils in high volume in a relatively short time of this landscape (Green et al. , 2. The epikarst layer requires a very long time to form. Gunung Sewu is one of Indonesia's karst areas that reportedly experiences recurrent droughts (Cahyadi. Marfai. Rahmadana, & Nucifera, 2012. Fatchurohman & Cahyadi. Nevertheless, in the event of high rainfall, floods can still occur and affect the Gunung Sewu is highly susceptible to the variability of climatic parameters, especially Apart from the strong influence on water resources, rainfall is also a crucial factor in crop and livestock agricultural practices, which are the main economic sectors in Gunung Sewu. Many of its inhabitants still rely on rain as a primary source of clean water throughout the rainy season, along with lakes and epikarst Water availability is much dependent on rainfall. Even for merely a few months, the absence of rains can lower the water levels and dry out most of these lakes and epikarst springs. For these reasons, a study of rainfall variability is relatively important. Gunung Sewu karst areais part of the Indonesian Maritime Continent that regulates many determinants of its rain variability, for ISSN: 2355-7079/E-ISSN: 2406-8195 example, the Asia-Australia monsoon. El Niyo-Southern Oscillation (ENSO), eastwest circulations (Zonal/Walker circulatio. , north-south circulations (Meridional/Hadley Circulatio. and some circulations that are heavily influenced by local factors (Hermawan. Aldrian. Gates, & Widodo, 2. Indonesia, the complex climate conditions are attributed to many factors, to which every part of the country shows varying responses (Gutman. Csiszar, & Romanov, 2000. Haylock & McBride, 2001. Boer, 2003. Athoillah. Sibarani, & Doloksaribu, 2. Other studies have confirmed that southern Indonesia, including the Gunung Sewu Karst Area, is inclined to have higher rainfall variability than the northern side (Juaeni, 2. The climate of Java is significantly controlled by the Asian-Australian monsoon (Ramage. Aldrian, 2. History has documented that ENSO contributes to shaping the rainfall on this island (Boer. Faqih, & Ariani 2. Further research stated that hills, mountains, and volcanoes in the southern and central parts of the island cause the southern side to be wetter than its northern counterpart (Qian. Robertson, & Moron, 2. In tropical studies, rainfall variability often takes precedence over other climatic parameters, given that it is the most influencing factor of plant productivity and water resources. The tropical condition may be linked to parameters other than rainfall that rarely experience similar large fluctuations (Naylor. Battisti. Vimont. Falcon, & Burke. Furthermore, comprehensive knowledge of rainfall variability helps manage water resources, understand and reduce disaster risks, assist in infrastructure work schedule, and plan the agriculture, fishery, livestock, and tourism sectors (Stefanidis & Stathis, 2. Rainfall variability analysis in the Gunung Sewu karst areais vital, but long-term daily rainfall data required for the study remain insufficiently In total, the Gunung Kidul Regency has 13 rain gauge stations, three of which are in Gunung Sewu and since 1979 have long timeseries of records with low quality as too many Rainfall Variability in Gunung Sewu Karst Area. Java Island. Indonesia. (Ahmad Cahyadi et al. data are missing: 42% of data is missing at Panggang, 46% at Tepus, and 38% at Rongkop. However, their daily rainfall availability . atio of filled to missin. shows that the data quality is At this state, rainfall analysis instead should use satellite data with relatively complete long-term records. Besides, many studies have relied on these global data, which also have easy accessibility and high accuracy (Christanto. Setiawan. Nurkholis. Sartohadi, & Hadi, 2. This research was intended to analyze the rainfall variability in the Gunung Sewu karst area from 1979 until 2013. A local-scale analysis is vital in that rainfall also varies depending on local factors (Hermawan, 2010. Hamada. Yamanaka. Matsumoto. Fukao. Winarso, & Sribimawati, 2. This research is expected to contribute to climate variability at the study site to formulate the most appropriate plans to meet water needs in karst areas in the future. Furthermore, it can be used to understand rainfall characteristics, especially in temporal studies, as reference material for mitigating meteorological disasters. II. MATERIAL AND METHOD This report's climatological data are rainfall data obtained from the Global Weather Data for the SWAT (Soil and Water Assessment Too. The data are accessible at https:// Analyses were then conducted to produce the desired data. Daily rainfall data were generated from the National Centers for Environmental Prediction (NCEP) Climate Forecast System Reanalysis (CFSR), global coverage and adequately high resolution. These data cover the years from 1979 to mid2013 . ata length= 35 year. The spatial resolution of the CFSR data is 38 km or 0. 3125 degrees. Therefore. Fuka. Walter. Macalister. Degaetano. Steenhuis, and Easton, 2. stated that this data is very suitable for studies on an area of about 40 square km2. This study only represents the Gunung Sewu karst area in the Gunung Kidul Regency. Yogyakarta. Dile and Srinivasan . stated that the CFSR data's accuracy compared to the rain recording data on the ground station is high. The difference between the CFSR data and the ground station is between 5. 8% and 57%, with an average of 16. However, the studies of data utilization from Fuka. Macalister. Degaetano. Steenhuis, & Easton . in the USA and Ethiopia. Dile and Srinivasan . in the Blue Nile River Basin. Ethiopia and Christanto. Setiawan. Setiawan. Sartohadi, & Hadi . in Serayu Watershed. Indonesia, showed an excellent performance even as an input to further hydrological models such as runoff modelling or water balance in a Here, the rainfall variability was determined by calculating mean monthly rainfall and analyzing trends of annual rainfall and mean rainfall intensity, dry and rainy seasons, the number of rainy days, and the effect of ENSO on rains in the Gunung Sewu karst Rainfall trends, intensity, and the number of rainy days were analyzed using linear In contrast, the effect of ENSO on rainfall variability was analyzed descriptively from the rain characteristics during the El Niyo and La Niya months or years. Data on these El Niyo and La Niya events were obtained from the National Oceanic and Atmospheric Administration website. US Department of Commerce, i. , https://noaa. In this Table 1. El Niyo and La Niya classification based on ONI values ONI Value Range 5 - 1. >1. 0 - 1. >1. 5 - 2. >2. Classification Week El Niyo Moderate El Niyo Strong El Niyo Very Strong El Niyo ONI Value Range 0 - <-0. 5 - <-1. 0 - <-1. < -2. Classification Week La Niya Moderate La Niya Strong La Niya Very Strong La Niya Indonesian Journal of Forestry Research Vol. No. April 2021, 23-35 research. El Niyo and La Niya occurrences were classified according to the Oceanic Niyo Index (ONI) values, i. NOAA's main parameter used to determine each. ONI values were measured from sea surface temperature in the Niyo 3. 4 region, where ONI > 0. 5 marks an El Niyo occurrence, whereas ONI = -0. 5 indicates La Niya. Table 1 shows the El Niyo-La Niya classification details used in this research. This research's data coverage represents the western part of the Gunung Sewu karst area, which includes Gunung Kidul Regency. Nevertheless, the study is expected to provide various related parties with a representative overview of the area as an input to their decision-making processes. This data coverage is considered adequate, given that groundwater resource and disaster management in Indonesia are currently part of provincial and local Discussion on water resources in the western part of the karst areais believed to give more focused details and thoughts as a contribution to Gunung Kidul Regency's governments and the Yogyakarta Province. RESULT AND DISCUSSION The analysis of rainfall data from 1979 to 2013 revealed that the Gunung Sewu karst ISSN: 2355-7079/E-ISSN: 2406-8195 area's monthly rainfall showed a monsoonal It shows the horse saddle-like shape in Figure 1. High rain lasted from January to March and November to December, whereas low rainfall occurred from MayAeSeptember. April and October were transitional months between these two seasons. The highest rainfall, 258 mm/month, was found in February, and the lowest, 15 mm/month, occurred in August. These characteristics support the argument that the monsoon winds and the Intertropical Convergence Zone (ITCZ) control the wet period in NovemberAeMarch (Asnani, 1. It is consistent with Aldrian . , who categorized the study area into the climate region AAialso known as the Australian Monsoon RegionAiwhere monsoon winds strongly affect the rainfall pattern. Another corroborating opinion is that Java Island, where the Gunung Sewu karst area is located, is the center of the Australian-Asian Monsoon Region (Ramage. Because of this, the monsoon wind tends to exert a substantial effect on the climate of the study area. The annual rainfall in the Gunung Sewu karst area was, on average, 1,831 mm/year. The highest rainfall, 3,885 mm/year, occurred in 2010, while the lowest was 909 mm/year, which appeared in 1982. Figure 2 shows an Figure 1. Mean monthly rainfall in the Gunung Sewu karst area Rainfall Variability in Gunung Sewu Karst Area. Java Island. Indonesia. (Ahmad Cahyadi et al. Figure 2. The trend of the annual rainfall in the Gunung Sewu karst area Figure 3. The trend of rainfall in rainy and dry seasons in the Gunung Sewu karst area increasing trend of annual rainfall by A 650 mm in 35 years. Prior scholars have also reported similar trends in neighbouring regions. Sleman Regency. Yogyakarta City, and Bantul Regency (Dipayana. Cahyadi. Mutaqin, & Nurjani. Dipayana. Nurjani, & Adji, 2012. , and Magelang Regency that is located relatively close to the study site (Suprayogi et al. , 2. However. Figure 3 shows that rainfall in the rainy and dry seasons increased at different magnitudes . he determination of the dry and rainy seasons follows (Brunsch. Stoffel. Ikhwan. Oberle, & Nestmann, 2. , in which the rainy season starts in November and ends in April, while the dry season lasts from May until Octobe. Rainfall increased more significantly in the rainy season than in the dry season, as indicated by the gradient of the higher in the former than the latter. The monsoonal pattern is observable from the amount of rainfall and the number of rainy days per month (Figure . From November through March, the study area had many rainy days, but this number was much fewer in Indonesian Journal of Forestry Research Vol. No. April 2021, 23-35 ISSN: 2355-7079/E-ISSN: 2406-8195 Figure 4. The number of rainy days in the Gunung Sewu karst area Figure 5. The trend of the number of rainy days in the Gunung Sewu karst area MayAeSeptember. Although the highest average rainfall was in February, the highest number of rainy days was in March, i. , 24 days. August was the peak of the dry season and had the lowest number of rainy days. It also corresponds to the mean monthly rainfall, which was the least in this month. The number of rainy days is crucial in the research, especially concerning agricultural activities and flood events. Rain lasting for many days is most likely to sustain plants' lives that depend on soil moisture, such as food crops and horticulture. However, too many consecutive days of rain can oversaturate the soil and lead to landslides and floods. Gunung Sewu karst area had 193 rainy days per year. The most extended wet condition occurred in 2019, with 313 rainy days/year, while the shortest one was in 1987 and 1984, with 140 rainy days/year. In detail, the rainy season had, on average, 132 rainy days, while throughout the dry season, rains still occurred during 61 days. Figure 5 shows an increasing trend in the number of rainy days in the Gunung Sewu karst The gently sloping trendline indicates no significant change for 35 years . n increase by Rainfall Variability in Gunung Sewu Karst Area. Java Island. Indonesia. (Ahmad Cahyadi et al. Figure 6. The trend of the number of rainy days during the rainy and dry seasons in the Gunung Sewu karst area merely A10 rainy days in this perio. However, after a further study of every season, it turned out that the number of rainy days increased only during the rainy season but tended to decrease in the dry season (Figure . The rate of increase was A 20 rainy days per 35 years in the former, while in the latter, it was A ten rainy days per 35 years. The Gunung Sewu karst area showed a steep increasing trend of rainfall but a less significant elevation in the number of rainy days . s evident from the gently sloping trendlin. This condition led to an increase in the daily rainfall intensity (Figure . High rain in the area has been reported to cause floods, such as the 2017 flood that coincided with cyclone cempaka (Cahyadi, & Mardiatno, 2019. Haryono et , 2020. Samodra et al. , 2. and the 2019 flood due to cyclone savannah-induced rainfall (Riyanto et al. , 2. In Gunung Sewu, some landforms are prone to this disaster, namely areas surrounding the ponor where the allogenic river enters the underground river system, springs or resurgences that are controlled by large conduit flows, and karst windows and closed basins . that have poor drainage because they are covered by sediments or have small-size sinkholes (Cahyadi & Mardiatno. Rainfall variability in the Gunung Sewu karst areais inseparable from the effect of ENSO. Figure 8 depicts the annual rainfall in Indonesia . lue lin. and its relation to El Niyo. From 1979 to 2013, there were eleven El Niyo years, including two very strong events in 1982/83 and 1997/98. In the other nine El Niyo events, rains were below the mean annual rainfall in the Gunung Sewu karst area . ed lin. other words. El Niyo strongly influenced the rainfall in the Gunung Sewu karst area, which is consistent with Boer et al. , who has found a close link between 41 droughts in southern Indonesia (Java. Bali. Nusa Tenggara. South Sulawesi island. from 1844 until 2009 due to El Niyo. El Niyo also affects Indonesia's agricultural sector. for example, in 1997/98. El Niyo caused damages amounting to USD 2. billion and a total economic loss of USD 9 billion nationwide (Bappenas, 1999. Kirono & Tapper, 1. Low rainfall as an impact of El Niyo will have a devastating effect on regions where most of the populations work in the crop and agricultural livestock sectors (Ayanlade. Radeny. Morton, & Muchaba 2. like the Gunung Sewu karst area (Cahyadi et al. , 2012. Cahyadi, 2014. This area is mainly used for dry crop cultivations and rainfed rice farming Indonesian Journal of Forestry Research Vol. No. April 2021, 23-35 ISSN: 2355-7079/E-ISSN: 2406-8195 Figure 7. The Trend of daily rainfall intensity in Gunung Sewu karst area Figure 8. Rainfall data of the Gunung Sewu karst area from 1979 until 2013 with the weak, moderate, strong, and very strong El Niyo years as the reference dependent on rainwater (Naylor et al. , 2007. Lestariningsih. Cahyadi. Rahmat, & Zein. Also, very low rainfall is most likely to induce secondary disasters, namely droughts, decreased agricultural production, and food shortages (Cahyadi et al. , 2013. Abaje. Sawa, & Ati, 2. La Niya also shapes the rainfall variability in the Gunung Sewu karst area. Figure 9 shows the relationship between annual rainfall . ed lin. and La Niya events in Indonesia. From 1979 through 2013. La Niya happened seven times or fewer than El Niyo. Figure 9 shows that during strong La Niya, the rainfall increased considerably compared to the average rainfall . ellow lin. However, unlike El Niyo, which has a strong influence, the weak and moderate La Niya years do not appear to be Rainfall Variability in Gunung Sewu Karst Area. Java Island. Indonesia. (Ahmad Cahyadi et al. Figure 9. Rainfall data of the Gunung Sewu from 1979 until 2013 with the weak, moderate, and strong La Niya years as a reference Table 2. Rainfall with different return periods in the Gunung Sewu karst area Probabilities (%) Return Periods (Yea. strong enough to affect the study area's annual As evidence, the annual rainfall during these events was still below the mean yearly rainfall line. Hamada et al. and Hidayat and Ando . also confirmed that El Niyo is more likely to affect Indonesia's rainfall As presented in detail in Table 2, the highest rainfall in 35 years of recording in the Gunung Sewu karst area was 161 mm/day, with a return period of 200 years. The return period for 50 mm/day rainfall was five years, meaning that extreme rainfall events of this magnitude can occur every five years. Consequently, extreme events pose high hydrometeorological hazards in the area. Rainfall . m/da. IV. CONCLUSION There is generally an increase in the mean annual rainfall trend, the number of rainy days, and daily rainfall intensity based on the analysis of rainfall data in the Gunung Sewu karst areafrom 1979 until 2013. In more detail, the number of rainy days in the rainy season tends to increase, while that in the dry season is decreasing. It creates a challenge to the area because the possibility of drought in the dry season will be higher. The yearly rainfall is always below the annual average during the El Niyo events, indicating a strong influence of El Niyo on the rain. On the contrary. La Niya has less effect, which is apparent from the absence Indonesian Journal of Forestry Research Vol. No. April 2021, 23-35 of high rainfall following this climatological Based on the analysis results, there are less rainfall and fewer rainy days in the dry season, but each rainy event occurs in higher intensity. this state, drought and flood prevention should be the primary concern of future improvements of relevant disaster management infrastructure. Considering that rainy conditions, particularly in dry seasons, will be increasingly unreliable in the future, providing infrastructure for clean water extraction from underground rivers that caters to all residents of the Gunung Sewu karst areais, in this case, is a suggested option. The analysis results show that for 35 years . 9Oe2. , there has been an increase in the mean annual rainfall of 650 mm. If this amount is properly distributed temporally, it will cause more water resources reserves in the karst area. However, this study's results also indicate that there is a trend of increasing daily rainfall intensity. This condition also means that the possibility of extreme rainfall is also getting higher. This means that the increase in the average annual rainfall that occurs does mean more water resources and the possibility of extreme rainfalls and floods that can occur more frequently. In the future, to complement this research, it is necessary to conduct further studies related to drought susceptibilities, climate change predictions, and analyses of extreme rain Studies surrounding their applications, , flood and drought impact analysis and food security, are also essential. They can at least contribute to long-term environmental management and disaster management. Although satellite data meet the requirements of rainfall variability analysis, field-recorded details remain necessary. Weather satellite data are not without limitations in that information with small to medium resolution is not sufficient for large-scale analysis. Field data available in long time series will facilitate and create a better Besides, the vast area of Gunung Sewu suggests that there need to be more rain gauge For an even spatial distribution, at least ISSN: 2355-7079/E-ISSN: 2406-8195 one station should be added to each district, particularly at a location selected based on the sea's distance to provide distinct variations between rainfall conditions along the coast and far inland. ACKNOWLEDGEMENT This research is a part of the first Author's dissertation in the Geography Study Program. Faculty of Geography. Universitas Gadjah Mada. It received funding assistance from the Faculty of Geography. Universitas Gadjah Mada through the 2020 Independent Research Grant Scheme for the Faculty Lecturers (Hibah Penelitian Mandiri Dosen Fakultas Geografi Universitas Gadjah Mad. under the Decree of the Faculty Dean No. 1291/UN1/FGE/KPT/ SETD/2020. REFERENCES