Endah Saptutyningsih1 Faculty of Economics. Universitas Muhammadiyah Yogyakarta e-mail: end_naufal@yahoo. Suryanto Faculty of Economics. Universitas Sebelas Maret e-mail: yanto. rimsy@gmail. Abstract The natural disasters that took place in many parts of Indonesia in the last 10 years have caused an enormous losses, including agricultural sector. The objective of this paper is to estimate the magnitude of the influence of flood disaster on land price changes in Yogyakarta, using a hedonic price The sample are is choosen from the vulnerable flood area mapped by Geographical Information Systems (GIS), from which farmers and land owner are selected as the respondents. The paper finds the evidence of a high level of flood stream coefficient, indicating that the flood significantly reduces the land price. The average of householdAos marginal willingness to pay for the decrease in flood stream level is estimated to be Rp 2. Keywords:A gricultural land, marginal willingness to pay (MWTP), hedonic price JEL classification numbers: Q15. Q19 Abstrak Bencana alam yang terjadi di berbagai daerah Indonesia dalam kurun 10 tahun terakhir telah menyebabkan kerugian yang besar di berbagai sektor ekonomi, termasuk sektor pertanian. Tujuan dari makalah ini adalah untuk memperkirakan besarnya pengaruh bencana banjir pada perubahan harga tanah di Yogyakarta, dengan menggunakan pendekatan harga hedonis. Sampel dipilih dari daerah rawan banjir yang dipetakan oleh Sistem Informasi Geografis. Responden yang dipilih adalah petani dan pemilik tanah. Makalah ini menemukan bukti tingginya koefisien aliran banjir, menunjukkan bahwa banjir secara signifikan mengurangi harga tanah. Rata-rata kesediaan marjinal rumah tangga untuk membayar upaya demi mengurangi tingkat banjir sungai diperkirakan bernilai Rp 2,175. Keywords: Tanah pertanian, keinginan marjinak untuk membayar, harga hedonis JEL classification numbers: Q15. Q19 INTRODUCTION1 Disaster risk is influenced by three factors which are the level of vulnerability, danger, and community capacity. The danger can be minimized if the vulnerability level can be reduced and the community capacity level Corresponding autor can be increased. In flood subscription area, the awareness of public on mitigating the risk might enable to minimize the loss. contrast, the low awareness would cause a very detrimental disaster effect. The natural disasters might affect the sustainability of development processes. The direct impact is the destruction of economic assets including stock of capital, infrastructure, natural resources, and human being. The disasters can also affect government budget. Hofman and Brukoff . states that when the disaster occurred then the government are faced by two urgent financings which are the operating fund and reconstruction fund. In its 1999 report, the United Nations Economic Commission for Latin America and the Caribbean estimated the impact of the direct loss of 25 major disasters occurred during a specific time period. The investigation shows that there have been serious damages caused by the impact of disasters in several states. Several studies have been conducted to assess the impact of natural disaster. The use of a hedonic price method to perform valuation has become a long debate especially since the study conducted by Ridker and Hennings in 1967. The debate is on the identification of functional forms which eventually shifts to a simple functional form The hedonic price method estimates demand and price of goods without traditional economic markets. This method is used to examine its constituent characteristics, and to estimate the value of each characteristic. In the case of land, which is a heterogeneous bundle of goods, researchers have used hedonic price to assess the implicit prices for characteristics associated with various properties, such as structural components, environmental factors, public services, and urban form. Kulon Progo is one of the regencies in Yogyakarta Special Region (Yogyakrta here afte. in which the flood reached 1,944 hectares and the number of inundated houses were 129 units. The data indicates that the substantial losses were mainly in the form of physical damage to homes, lands, farms and forests (Table . Many studies on flood disaster have been done using various methods of analysis such as. Cost and Benefit Analysis, econometrics. Geographic Information System. Data Envelope Analysis, and hedonic Table 2 presents some summary of previous related studies. Dai et al. mention that high rainfall will increase the risk of landslides, slope instability as well as landslides risk of vegetations. It is stressed out that landslides is associated with rainfall levels. Meanwhile Parson et . use GIS to change the mitigation Cowell and Zeng . also use GIS as a modeling-prone area due to weather changes. The GIS method has been widely used to identify areas of potential disasters. Wood and Good . use GIS to identify vulnerabilities at airports and seaports by the earthquake and tsunami. Rashed . measures the environmental context of social insecurity by the earthquake. Zerger . uses GIS to examine the risk disaster model and Cowell and Zeng . integrate the theory of uncertainty by as a modeling-prone area due to weather Table 1: Damage and Casualties of Flood in Yogyakart Special Region Over the Last 10 Years District/city Victims Yogyakarta Bantul Gunung Kidul Kulon Progo Sleman Houses . Damages Land . Source: Data of Information and Business Indonesia. DIY 2009 Garden/forest . Hedonic Price Approach A (Saptutyningsih and Suryant. The purposes of this study are to map flood area especially the agricultural land in Yogyakarta Special Region and analyzes factors that influence agricultural land price in the flood area. Based on the prone area map, the paper will make the economic valuation of loss suffered from flood in the most vulnerable area by using a hedonic price method. This study also aimed to evaluate how far the effects of the flood on land prices are concerned. dict flood inundation caused by river overflow. Water propagation will be modeled by iterative calculations. Iterative calculation using the value of DEM is to identify the nearest neighbor pixels. Nearest neighbor pixels will be recognized as inundation areas based on the scenario. In the scenario of flood inundation as high as 1 meter high, the iteration process will stop at the DEM pixels that is worth 1 meter higher than the initial surface water. METHODS Geographic Information System To identify the flood prone area in Yogyakarta, this paper uses the Geographical Information Systems (GIS). To construct the flood modeling, this paper uses neighborhood model operation with an iterative This method is adopted from Marfai . model in the case of flood modeling study in Semarang. An iteration method is a mathematical model using the calculation result of iteration as input for subsequent calculations. The calculation is done in the form of raster spatial data. Raster data is represented with a square shape with a certain For example, if one zooms in an image on a computer, the image will be difficult to be recognized. The appearance of existing will be only in the form of boxes that have different colors. The color indicates a value. In the field of cartography, raster data is often encountered in remote sensing image. In addition, raster data is also widely used to represent the height of a place which is represented by a DEM (Digital Elevation Mode. In DEM, elevation values will be represented somewhere in each picture box, or often referred to as pixels . icture elemen. The value of each pixel represents the value of altitude on the earth surface. Figure 1 depicts the illustration of neighborhood operations and iteration procedures. In this study. DEM is used to pre- Source: Marfai, 2003 Figure 1: Neighborhood Operations The data calculation in Geographical Information System applications is hdone with the help of ILWIS 3. ArcView 3. and ArcGIS 9. 3 software packages, including a Point Interpolation Altitude. Vector Data to raster conversion, and the Iteration Technique-Neighborhood Analysis Flood Vulnerability Index Scoring Model The emphasis of this second method is to identify flood areas based on several factors that are owned by the region or study The data analysis techniques with this method is similar to the first method. The variables used in this method is not only based on the high scenario puddle and altitude regions, but by using the basic data (Table . , one can map to identify regions most catastrophic floods potential. These basic data which carry out layers of various parameters can indicate the potential for flood. In addition, by scoring and weighting of these parameters, one can get the index of vulnerability to flood. Flood-prone areas can be identified according to the scores and weights that have been made in which Vulnerability = . weighting scor. This formula enables cross-mapping on the basis of the data (Figure . Table 2: Basic Data for Mapping PotentialF lood Areas Type of data Land Use Administrative Boundaries Topographic Landform Infiltration Soil texture Slope Flood events Sources Bakosurtanal Bakosurtanal Bakosurtanal Bapeda DIY Bapeda DIY Bapeda DIY Bapeda DIY BPSDA-PU Scale 1:25. 1:25. 1:25. 1:250. 1:250. 1:250. 1:250. 1:250. Figure 2: Chart of Flood Vulnerability Index Scoring Method Year Hedonic Price Approach A (Saptutyningsih and Suryant. Hedonic Price Model This study uses two hedonic price models. The models use variables related to land The model is LP = 0 1 x1 2x2 3x3 LP is Land price x1 is Vector of property . characteristics x2 is Vector of neighborhood characteristics x3 is Vector of flood risk characteristics x4 is Vector society awareness x5 is Vector socioeconomic characteristics e,u is standard error The data is needed to conduct mapping of the flood prone areas including . Map of Yogyakarta using RBI scale 1:25. Geological Map scale 1:100,000. Data of height point of Yogyakarta Special Region. Information of road. Information of River. Information of land use. Data of flood inundation in Yogyakarta Special Region. Data of River Flow. The other type of data for economic valuation is the primary data resulted from the respondents, which are previously taken by the mean of distributed questionnaire. This paper applies a strategic random sampling method to get the sample from which the primary data are obtained. The respondents consists of agricultural land owners in the research location. The steps for collecting data are focus group discussion, pre-test, and final survey RESULTS DISCUSSION The flood in the province of Yogyakarta was concentrated in Kulonprogo Regency. The areas prone to flood is those located in southern county district of Kulonprogo i. Temon. Wates, and Panjatan. This can be observed from the height of the low area . -10 m above sea leve. Kulonprogo Regency along with Oyo River are also considered as sensitive areas. Although it is not widespreaded, but the area along the Oyo River must be alert to the flood hazard. Results from the tentative mapping of flood inundation by using iteration method is cross checked by the flood inundation data published of Board of Water Resource Management of Opak-ProgoOyo River. The result indicates that there are some areas prone to the flood. However, the are cannot be included as areas prone to flood by the iteration method. the example is the flood caused by the overflow of Opak River soaking about 10 hectares of paddy fields and 12 hectares of The inability of the iteration method to map the flood is caused by the lack of the stretcher on the height data in a high detailed scale. Iteration methods heavily rely on data which is represented by the DEM altitude (Digital Elevation Mode. The more detailed the elevation data is used as input, the more detailed the flood mapping that will be produced. DEM generated in this study cannot represent micro relief thus it cannot represent the localized flood. Most of the residential population occupies with clumped pattern. The characteristics of lowland areas tend to be flat and it facilitate the population to exploit the land for settlement. Based on field observations in Temon. Panjatan, and Wates Districts, the location of densely residential and clustered are in village centers. In Temon District, population are clustered in the Glagah and Kedundang Villages. Population in the Panjatan District are clustered in Tayuban and Wates Villages. Wates District which is the capital of Kulonprogo Regency has the highest concentration of population located in Wates Village and Village Dam. Mapping using Flood Vulnerability Index Scoring Method produces vulnerable area among other areas in south of Kulon- progo Regency area such as the Temon. Panjatan, and Wates Districts. There is also identified areas other than Kulonprogo regency that is areas along the River Oyo of Bantul Regency. Vulnerability to flood risk in each zone are varied and divided into four levels of vulnerability scenarios. The scenarios are water height 25 cm, 50 cm, 100 cm and 150 cm. Based on observations the most vulnerable area are the settlements and rice field area. In the scenario of 25 cm flood inundation. Plumbon and Kalidengen are the area the most prone to flood risk. The characteristics in these two regions are also mostly rice fields and un-irrigated land. the scenario of 50 cm, flood inundation begin to overflow into areas around Kalidengen and Plumbon which are mostly rice field and un-irrigated land. Inundation scenarios of 100 cm and 150 cm will reach Sogan region, eastern Kuwaru and western Palihan. From the analysis results it can be seen that the most vulnerable area in Temon District is the irrigated fields. While that of for Wates and Panjatan District are the moor and garden. From the analysis of population density level and settlement it can be noted that flood in settlements will reach 34 hectares. Of the two methods of analysis, it can be concluded that the highest floodprone areas in the region Yogyakarta Special Region is Kulonprogo Regency especially in the three districts are Temon. Wates, and the Panjatan Districts. So the primary data are collected from these districts. Table 3 provides and overview the descriptive statistic of the variables used in the analysis. Table 3: Description andS ummary statistics of Land Price -Related Variables Variables Land Characteristics Description of Variables Mean Standard PRICE COST YIELD OWNER Price of land . Cost of production . Amount of harvest . Ownership status D_HEALTH D_SCHOOL D_HALTE D_HDRAIN D_DRAIN Distance to health center . Distance to elementary school . Distance to bus stop . Distance to drainage . Distance of drainage . PERCEPT LONG FAILED LOSS HEIGHT Perception of Flood Time period of flood inundation . Failed of harvest Losses caused by flood Height level of flood inundation . AGE SEX EDUCATE CHILD INCOME Age of respondent Sex . ale or femal. Education level Number of children . Income level per month Neighborhood Characteristics Flood Risk Characteristics Socioeconomic Characteristics Source: Data processed Hedonic Price Approach A (Saptutyningsih and Suryant. The price of land in the area of the sample are still a reasonable price. Mostly the land is owned by the farmers. The average yield for one-time harvesting is 2. million rupiah. Characteristics of neighborhood variables are the distance of the area with the health, educational, and bus stop facilities, the distance to the nearest river, and the distance from drainage. On average, the distance of land to health facilities is 731 meters and the distance to education facilities is 1. 38 km. The location of the observed land is not too far away to rivers and drains i. 166 meters and 49 meters. Length of flood inundation is approximately eight days with the 542,549. 00 IDR crop losses. The socioeconomic characteristics of respondents can be reflected by age, sex, educational level, number of children, and income level. In general the education level of respondents is graduates of Junior Secondary School (SMP) and the level of their average income is less than one million rupiah per month. The number of children is also relatively low at one child. Hedonic price method has been used to study natural disasters. The underlying assumption of the revealed preference are based on housing market prices in which the purchase option of a house is based on structural characteristics, as well as characteristics associated with the location of the agricultural land. Each of these characteristics has a value and these values are added into the transaction price that is being observed. Table 4: Hedonic Valuation of Marginal Willingness to Pay for Flood Mitigation Related Variable Constant to Land Value in Kulon Progo Regency Full Model Fit Model 273 . OWNER LAREA 283(-1. 196(-2. LNCOST 089 (. 164 (-2. LNYIELD 014 (-. D_HEALTH 000 (. D_SCHOOL 010 (-. D_BUSSTOP 001 (-. 380 (-2. D_HDRAIN 000 (. D_DRAIN 000 (. PERCEPT 027 (-. LONG 007 (. FAILED 376 (-20. *** 307 (-1. LOSS 000 (. AGE 004 (-. SEX 002 (-. EDUCATE 035 (. CHILD 059 (. INCOME 081 (. HEIGHT 004 (-2. 005 (-2. No. Observation Adjusted R2 Note: ***, **, * indicates significant at 1%, 5%, 10% levels of significance, respectively. Source: Data estimation Column . of Table 4 indicates that the production cost has a negative effect on land prices. The higher the production cost is, the lower the price of land. indicates that the land buyers will consider cost as one of factors to decide the price of Higher production cost will reduce their future profit. Distance to the center of education and bus stops negatively affect the price of Further from the health services and the bus stop will lower the price of land. Meanwhile, crop failures affect the price of This can be ascertained that the success of land to produce crops would affect price of those land. Inundation height also has a negative effect on prices of land. This is because the prospective buyer of land would have to consider the risk of any loss due to flood inundation. Therefore, the higher of flood inundation, the lower the price of those land. The result shows that all coefficients are significantly different from zero. All coefficients have expected sign and their magnitudes are comfortably reasonable. The marginal implicit prices of land characteristics are not constant due to the nonlinearity of the hedonic function. The marginal implicit price of the variable of interest. HEIGHT, is calculated by differentiating the hedonic price function with respect to HEIGHT. Therefore, for a given household, each unit increment in height of flood inundation results an estimated decrease in land unit value as Ae0. 004 times land unit value divided by the associated height level. For example, a household whose land unitAos value 12. 000,00 IDR with associated level of height at 50 centimeters will suffer a marginal damage 000,00 IDR. The same householdAos marginal benefit of a unit decrease in height level of flood inundation is equivalently 1. 000,00 IDR. In other words, this particular household is expected to be willing to pay no more than 1. 000,00 IDR to avoid the damage to their land associ- ated with a one unit increase in height level of flood inundation. Denoting the size of the sample by N, the sampleAos average MWTP is calculated as the average of marginal damages as follows: MWTP = HEIGHT HEIGHT Thus, if the sample is fairly representative of Kulon Progo Regency, an average household should be willing to pay 00 IDR for a unit decrease in height level of flood inundation. While an average household should be willing to pay 175,00 IDR for a unit decrease in height level of flood inundation. Needless to say, the benefit is capitalized in the property, and the associated MWTP is not to be thought of as a payment per some unit of time, but rather as a lump sum payment for moving from the status quo to a lower level of flood inundation. There are several variables that theoretically affect the price of land. However for the following flood-prone areas such variables are not significantly affected. This is possible because people underestimate of the risks caused by flood. Moreover, socioeconomic variables that do not affect the price of property and land might because the perception of the community that think the flood was a common. Therefore, the marginal WTP of their flood mitigation of the impact was quite low. CONCLUSION Using a Geographical Information Systems (GIS) approach, the flood was concentrated in Kulon Progo Regency of the Province of Yogyakarta Special Region. Areas prone to flood located in southern county of Kulon Progo Regency i. Temon. Wates, and Panjatan Districts. These area experiencing severe flood every year need a special treatment such as by optimizing channels Hedonic Price Approach A (Saptutyningsih and Suryant. of existing water. Many channels of water are not optimally functioned. It also required an evaluation about the volume of water that does not fit into drains or causing flood. The serious attention is required to cope the floods in the future. It related to the increasing population while it not accompanied by the improving environmental carrying capacity. The decreasing environmental carrying capacity cause higher vulnerability to flood. The results of this study indicates that the hedonic price method can be applied in flood case of Kulon Progo Regency. It is to obtain the estimates of mar- REFERENCES