Iran is a low-water country located in an arid and semiarid region of the Earth. Therefore, snowfall is considered as the most important factor affecting water supply, especially in the summer, and also snowfall is one of the important factors that controls the hydro-climate of the geographical area. In northwest Iran, precipitation occurs more frequently in the form of snow in the cold seasons, and due to the mountainous characteristics of this area, snow resources persist throughout the year. Consequently, the scientific management of water supplies in the form of snow is necessary to maintain a supply of consumable water.1,2
The development of predictive models of runoff resulting from snowmelt requires the determination of many spatial parameters. Knowledge of the distribution of precipitation, and especially the spatial and temporal distribution of snowfall and soil moisture, is required for the characterization of runoff.3 However, these data are often scattered and of poor quality, but remote sensing has good potential to eliminate these shortcomings.
Snow measurement and assessment in mountainous areas are difficult due to lack of appropriate access. The diverse features of topography and physiography also complicate the estimation of hydrological snow parameters and the inadequate numbers of weather stations in the highlands point to the necessity of using indirect methods applications, such as the analysis of remote sensing data for measuring snow cover area (SCA).4 Satellite images and application of the available algorithms allow extraction of SCA with a reasonable accuracy. The estimation of snow cover surface significantly improves the process of calculation and application of snowmelt models using remote sensing images.56.–7
However, these issues have already been partially resolved by the MODIS instrument onboard the environmental remote sensing satellites of the Earth observing systems (EOSs) Terra (EOS-AM1) and Aqua (EOS-PM1), which were launched successfully in 1999 and in 2002, respectively.89.–10 The MODIS Terra satellite crosses the equator at about 10.30 a.m. (descending orbit) UTC, with a scanning swath of 2330 km (cross-track) by 10 km (along-track at nadir). The MODIS has a total of 36 different wavelength channels suited for a wide range of applications. An algorithm has been developed by Hall et al.8 to map snow cover at 500-m spatial resolution using the MODIS observations. The 500-m spatial resolution results from the use of the MODIS bands with demonstrated capability for detecting snow and separating snow from clouds.11 Assessment of the MODIS snow product accuracy and execution of applied research is now possible,12 and many scientists are now focusing on the study of snow and ice (e.g., Refs. 13188.8.131.52.18.19.–20).
These studies have examined the relative accuracy of snow mapping for various snow cover products or the relationship between snow mapping accuracy and spatial resolution, ruggedness, snow depth, or land cover type. The snowmelt runoff model (SRM) has been successfully tested in numerous mountainous basins around the world (e.g., Refs. 2184.108.40.206.26.27.28.–29). The main advantage of temperature index models (i.e., SRM), over other model types, is that they only require air temperature data to obtain estimates of melt.30 The wide availability of these data (mainly in mountainous areas) has led to the popularity of temperature index models, especially in operational stream flow forecasting.31
The purposes of this research are: (a) to investigate the capability of MODIS-8 daily images at 500-m resolution to provide the map of SCA32 and recognize changes in snow cover using the normalized difference snow index (NDSI); (b) to compare snow cover maps extracted from MODIS with the ground-based snow courses at some selected points in the central Zab basin to yield snow water equivalent (SWE); (c) to perform an SCA validation using the MODIS products during the water year 2010 to 2011 for accumulation and ablation periods; and (d) to use the SRM model to model the snowmelt runoff in the study area.
The study area is located in the southwest mountainous area of West Azerbaijan province, Iran, along the Zab river basin in Sardasht, between the latitudes of 36° 8’ 25”N and 36° 26’ 27”N and the longitudes of 45° 21’ 21”E and 45° 40’ 44”E. The central part of the Zab river basin stretches 30 km from north to south and east to west (Fig. 1).
The whole study area is ; 30.45% is located in the west and 69.55% is located east of the main river.33 Wherever the size of the catchment of waterways is greater, the rate of flow and volume of alluvial cone material will be higher. In western part of central Zab basin, the rate of flow and volume of runoff is low.
The average annual rainfall in the central Zab basin is 865 mm. The highest rainfall in the study area occurs in the months of March and April; about 53% of the annual rainfall occurs during the months of March to May and 9% occurs from June to October. The climate condition of central Zab basin is semiarid with strong vertical gradients in precipitation, so the area is cold. A high percentage of the annual precipitation, especially in the autumn and winter seasons, occurs as snow. Due to the prevailing synoptic weather systems in the region and the local topography, snow is a stable element throughout the year, especially in the highlands area.
Materials and Methods
The application of satellite images to model runoff processes obtained from snowmelt was evaluated using the MODIS satellite images obtained on 8 different days that signal quantization is 12-bit radiometric. These satellite images have been used because of their ability to provide high temporal separation with a one day return period as well as the high spatial resolution similar to other satellite images.
The data used in this study include: Topographic maps and field data recorded at the snow gauge stations, including Sardasht, Brisweh, and Mirabad (see Fig. 2). The SWE was considered for a fixed water year (2010 to 2011), that starts from January until May. Due to closure during the melting season and the limited time for accessing data from the snow gauge stations, the daily snow cover maps received from MODIS were provided by the Iranian Space Agency with a L1B process level, and are used in deriving the snow depletion curve, which is one of the input parameters of the SRM (Table 1).
The dates of images used for the central Zab basin study.
|1||January 14, 2011|
|2||February 6, 2011|
|3||February 14, 2011|
|4||February 26, 2011|
|5||March 8, 2011|
|6||March 21, 2011|
|7||April 4, 2011|
|8||May 3, 2011|
The range of the basin and the snow line were determined using base data and topographic maps with a scale of 1:50,000 in digital format. The DEM data were generated from advanced space borne thermal emission and reflection radiometer (ASTER) images provided by the Iranian Space Agency. Images were acquired near the nadir of the 3N band (0.78 to 0.86 μm), and stereo images were acquired for the same area after examining the 3B band (0.78 to 0.86 μm). The ASTER 3N and 3B images were designed for the DEM creation with a Baseline/Height (B/H) from the nadir and aft-telescope of . One hundred tie points were used to form an epipolar image, which had a . parallax comes from a correlation analysis between the 3N band image of the left-hand side and the 3B image of the right-hand side. Therefore, the final DEM is extracted from visible and near-infrared ( to 0.86 μm) in ASTER with a spatial resolution of 15 m. The root of the mean square error (RMSE) of the DEM was . The DEM was generated using the bilinear interpolation method. Figure 2 shows the central Zab basin DEM extracted from the ASTER image.
The stages of data preparation were performed using Arc GIS 9.3 software as follows: changing the data format, layer integration, applying the geometric correction, and extraction of hydrologic parameters (such as exterminating the basin area range and waterways network). In the processing stage, the MODIS sensor images were extracted during two stages of snow cover: first, in the programming environment of the ERDAS software on the image of the B1 surface, earth operations. The images were referenced with the “Lambert” system and radiometric correction was conducted. The automatic algorithm for snow cover was implemented and the extraction of the SCA was started. In the process of the image processing operation, the following steps were implemented:
In the first stage, reflective criteria were applied if NDSI . In other words, if the reflectance of pixels in band 2 was more than 11% and the reflective value of pixels in band 4 was equal to or more than 10%. Based on above mentioned algorithms, conditional triple tests were conducted to extract the NDSI values.
In the second stage, the majority filter was implemented with the aim of reducing the reflectance role of a single pixel. This filter is a window of and its application classifies the valuable pixels of the specified NDSI as snow cover. Figure 3 shows the snow, snow free, and cloud covers on February 14, 2011, as 185, 238, and , respectively.
Reflective Feature of Snow and the Index of Snow Cover
The reflectance of a snow surface is affected by changes in factors, such as snow grain size, form, water content, surface roughness, snow depth, and impurities as well as the angle of incident solar radiation and reflection. Increasing snow age reduces its reflection in the visible spectrum and near-infrared regions, and the main factor in this reduction is an increase in snow grain size due to its melting and refreezing.34 The long and short waves received at the snow surface can be an important source of the snowmelt process. The most important feature of snow that causes a difference in the spectral reflection of snow is the difference in the VNIR spectral range.35
A study by Klein et al.36 showed that the capability for displaying snow using the MODIS sensor has made significant progress compared to other similar sensors in the unique and high resolution bands. The obtained images have a wide spectral range and the techniques used in the algorithm related to the snow map are based on detection at the local and regional scales.36 Processing is carried out using the images that have characteristics, such as L1B, no clouds, and taken in daylight are selected to extract the snow cover. These pixels must have been recorded during daylight with sufficient light and sky without cloud.37
In this research, the algorithm for the spatial snow cover detection provided by the United States space agency is used so that that the bands 1 to 6 of the MODIS sensor data were studied and analyzed. After applying radiometric and geometric corrections, MODIS/Terra Snow Cover Daily L3 Global 500-m Grid (MOD10A1) is used to snowmelt modeling.
The NDSI is a spectral ratio that uses the spectral difference of the infrared and visible bands in the MODIS sensors to recognize changes in snow cover surface.38 The NDSI index takes advantage of the fact that spectral snow reflection has a high reflection in the visible band and low reflection in the spectral range of infrared as a way to discriminate snow from cloud and snow-free areas.39 This index, like many spectral ratio methods, reduces atmospheric effects.8
The NDSI index can be used to separate snow and ice from each other and also to separate snow from the clouds at a high atmosphere altitude, like cumulonimbus clouds.41,42 This index is insensitive to the range of lighting conditions and can be adjusted for atmospheric effects. In other words, the NDSI index does not just depend on reflection values in a band but it also depends on the digital number of pixel reflections.43 Pellika and Rees44 proved that pure snow has a high NDSI, but when it is mixed with other substances (such as dust and smoke), the percentage of its purity is reduced. The NDSI index is also used as a separator of cloud-snow from clouds that reflect both VNIR light.
Hall et al.37 proved that the algorithm of a snow map is best operated in the condition of full snow cover in areas with thin vegetation, such as grasslands, agricultural fields, and tundra. In these conditions, band 2 of MODIS is mainly processed to detect snow and the NDSI components of the snow map algorithm effectively filter clouds (with the exception of high altitude clouds). These clouds contain ice pieces and may cause a misclassification of snow cover. According to this criterion, the results of the NDSI index can be accepted if the reflectance of band 2 is more than 11%. The second criterion is introduced as dark targets by Klein et al.36. In this mode, 10% reflection in band 4 is set as the low limit of identification and separation of vegetation from snow. For pixels classified as snow, reflection in band 4 must be equal to or more than 10%. Despite the high value of the NDSI index, in some cases, dark targets prevent correct snow classification. Thus, based on these two criteria, the algorithm for snow cover will consider a pixel as snow if the following conditions are achieved.
This algorithm is used along with two other thresholds:
1. 1-band 2 (0.841 to 0.867 μm) has the reflectance of more than 11%.
2. 2-band 4 (0.545 to 0.665 μm) has a reflectance equal to or more than 10%. In total, the NDSI index must be estimated as more than 0.4.
The NDSI index of snow cover is only sensitive to certain types of clouds that have ice particles and some clouds that have reflections similar to that of snow.45 Note that cloud typically has high reﬂectance in the VNIR wavelengths 4 and the value decreases slowly in shortwave infrared wavelengths (SWIR), while snow has high reﬂectance in visible and NIR and the value decreases signiﬁcantly in the SWIR (shown in Fig. 4).3)]: 47 (1) The NDVI values are , (2) the NDSI values are , and (3) the reflectance in MODIS band2 is higher than 0.75. According to these rules, the snow and nonsnow pixels can be extracted.46
Reflectance derived from radiance measured by sensors at the top of the atmosphere (TOA) may be increased or decreased when compared to the surface reflectance as a function of the target’s reflectance, its environment, sensor spectral band, viewing and solar geometry, and atmospheric characteristics.48
The snow cover algorithm developed by Hall et al.37 uses the TOA reflectance in the NDSI and other thresholds. In Hall et al.’s algorithm, the effects of atmosphere have not been corrected for, which may cause some errors for estimating snow cover in mountainous areas.49 Taking into account the computational speed and accurate technique for the atmospheric correction, the surface reflectances are retrieved from the TOA reflectances by using an updated simplified method for the atmospheric correction model in this work.50 If is the spectral surface reflectance of the target, surrounded by a homogeneous environment of spectral reflectance , the TOA spectral reflectance, at the satellite level can be expressed as [Eq. (4)]:4) can be simpliﬁed to [Eq. (5)]:50
Analysis of the Results and Discussion
Geometric Corrections of Satellite Images
Given the 500-m separation capability and considering that each pixel of the image covers a surface of , the geometric correction of the satellite-based MODIS data is very important when estimating the snow cover surface.51 Using appropriate ground control points (GCPs) for both positioning and number,52 it is possible to rectify images in a very simple manner by ensuring that the resultant RMSE is under a certain threshold.53,54
All rectiﬁcation methods produce a residual in the and axes between input and output coordinates; for each GCP, RMSE, e.g., the resultant vector from residuals in and axes is calculated based [Eq. (6)]:55
Total RMSE is then derived as [Eq. (7)]
Nearest neighbor resampling method simply associates the nearest original pixel brightness value to the corrected image pixel. For each area and rectiﬁcation method, total RMSE was calculated, and the dispersion cloud of GCPs residuals was plotted in order to verify the relation between terrain roughness and cloud dimensions.55
The nearest neighbor resampling method with the output pixel dimension of 2 m was performed on the SPOT 5 satellite image taken in 2008. The three-dimensional image of the central Zab Basin, extracted using the SPOT 5 satellite, is shown in Fig. 5.
Figure 6 shows an accuracy test of the process of geometric correction using the layer of basin waterways before (A) and after (B) the process of geometric correction.
The RMSE produced in the earth process of referencing MODIS satellite images is given in Table 2.
The root mean square error (RMSE) produced in the earth process of referencing moderate resolution imaging spectroradiometer (MODIS) satellite images.
|Year||Moon||Day||Number of control points||RMSE (m)|
Reflectance of Elements of Earth and Snow in the MODIS Satellite Images
Snow monitoring due to spatial and temporal snow variables depends on the remelting and reraining conditions. In addition, due to the lack of snow telemetry, the phenomenon mentioned in the result of map preparation, geometric accuracy is not possible when studying the accuracy of digital interpretation.
Providing the Images of Snow-Covered Area
After calculating NDSI and NDVI, the snow cover is determined in pixels according to these two values and other distinct criteria. Note that the reﬂectances used are the atmospheric corrected reﬂectances. To prevent errors in pixels containing very dark targets, such as black spruce forests, the threshold values of the surface reﬂectances in MODIS bands 2 (0.841 to 0.876 μm) and 4 (0.545 to 0.565 μm) are adjusted to . The NDSI was applied as an equation on the image so that the numbers higher than and equal to 0.4 were isolated as a conditional snow pixel and reﬂectance in MODIS bands 2 and 4 are . The currently proposed snow cover algorithm for MODIS data is shown in Table 3.
The proposed snow cover algorithm for MODIS data.
|Type||Proposed snow cover algorithm|
|Input data||MOD02-L1B, MOD10A1, DEM|
|Models||NDSI, NDVI, NDCI|
|Constraints of threshold||For nonforest area: NDCI >0.4 and cld-msk≠0 and NDSI≥0.4 and Pband 2>0.09 and Pband 4>0.09 For forest area: NDCI >0.4 and cld-msk≠0 and NDSI<0.4 and Pband 2>0.09 and Pband 1>0.09 and −0.4NDSI+0.4≤NDVI≤3.5NDSI|
P is surface true reﬂectance.
As depicted in Fig. 7 in red, the basin surface is considered as snow using this index. Therefore, to separate water from snow and also to separate dark and shadowy areas from snow, the reflections of band 4 were used for dark targets and band 2 for water.56
According to operation in the pixel value or by recoding the reflectance of band 2 higher than 11% and the reflectance of band 4 that was equal or , these values were multiplied in the NDSI image to obtain the final image of the surface snow cover (see Fig. 8).
Comparing the real images of MODIS and the image of snow separated by the snow normalized index gave the snow difference. This index separates the snow surface well (Fig. 9). In Iran, mapping of the snow surface is not conducted so that we cannot estimate the accuracy of this index; however, according to a study conducted in the Rio Grande watershed on MODIS snow cover during the years of 2000 to 2001, assessment of the quality of MODIS snow cover maps is accomplished through comparisons of MODIS snow maps with those produced by NOHRSC as well as through comparisons with in situ observations of snow depth calculated from automated measurements of SWE. As the entire snow season is being considered in this study, comparisons of MODIS snow maps were created from higher resolution sensors, such as Landsat ETM+ or ASTER. These comparisons are, however, being conducted as part of ongoing MODIS validation efforts.13
In this study, runoff from snowmelt was modeled by dividing the central Zab basin into 5 elevation zones with around 200-m classes of elevation surfaces in the Arc GIS 9.3, and the snow-covered surface was obtained in each of these separated zones. The amount of SCA in the central Zab basin, using NDSI extracted from the MODIS satellite images, is given in Table 4.
Amount of snow cover area in the central Zab basin using the normalized difference snow index (NDSI) extracted from MODIS satellite images.
|Zone||Area (km2)||Mean elevation (m)||Snow covered area (%)|
|Date of measurement (2011)|
The snow depletion curves belonging to each elevation class obtained from the MODIS snow products for central Zab basin are shown in Fig. 10.
Interpolation of Snow-Covered Surface
The snow-covered surface was then interpolated for days without satellite images. The snow-covered surface in the central Zab basin was calculated through the interpolation of the time series of MODIS images based on the algorithm provided by Malcher and Rott.57 The satellite carrier MODIS sensor takes images from the surface of earth daily but due to the cost of satellite imaging and time-consuming nature of the process and daily extraction of snow-covered surfaces in watersheds, the possibility of the daily use of satellite images was virtually ruled out. The complex nature of the snow and its location and time variation meant that the selection of dates of the satellite images was very important and could affect the accuracy of determining the snow-covered surface. Temperature was an additional factor of vital importance, with regard to the process of snow melting in the melting season, as pointed out earlier. The temperature factor has a greater effect than other phenomena, such as rainfall on the snow ridge or the quantities of radiation, in the process of melting snow on snow ridge. Extrapolation of the snow-covered surface is determined based on the temperature of the snowfall in each of the height zones. According to the Malcher algorithm, the melting factor is a function of critical temperature and a degree–day factor and according to Eq. (8), this can be expressed as58
Finally, taking all these factors into account, the snow-covered surface for three snow gauge stations including Sardasht, Brisweh, and Mirabad are obtained with interpolation of the extracted snow-covered surface derived from the MODIS satellite images (Table 5).
The snow-covered surface for three snow gauge stations obtained from MODIS satellite images.
|Snow gauge station||snow-covered surface (km2)|
Simulation of Snowmelt Runoff
In this research, the snowmelt runoff in the central Zab basin was modeled using the SRM model that provides the snow-covered surface through satellite images and enters it into the model as one of the basic variables. If rain falls on this snow cover, it is assumed that the same amount of water is released from the snowpack so that rain from the entire area is added to snowmelt. The melting effect of rain is neglected because the additional heat supplied by the liquid precipitation is considered to be small. A few models in hydrology also use remote sensing data as input data. The SRM is one of the models that uses remote sensing data directly from hydrological models.60
The model was proposed for the first time by Martinec and Rango60 for the proper management of water resources in the mountainous areas of the Alps, and its Windows version was upgraded in 2003. The main equation of the model, which allows the transformation of water produced from snowmelt and rainfall into the daily discharge, is based on a simple degree-day method.59
In the SRM model, discharge from rainfall and snow melting is calculated on a daily basis. The basic equation of the model is as follows Eq. (10):61
The degree-days are used to determine the melt rate of the snowpack in the area of the basin covered by snow as observed from satellites. The degree-day factor is not a constant, and changes in response to changing snow properties and atmospheric conditions. Martinec59 developed the following relationship for the computation of the degree-day factor (), which takes into account the positive correlation between snow density and the degree–day factor Eq. (11):12):62
Along with the topographic, meteorological, and hydrologic parameters, snow melt simulation requires 11 parameters and variables, as follows: recession coefficient, time lag, critical temperature, runoff coefficient for snow, runoff coefficient for rain, degree-day factor, rainfall contributing area, reference elevation, temperature lapse rate, precipitation lapse rate, and initial discharge; these were all considered as input parameters to the SRM model. The parameter values used for the central Zab basin are given in Table 6.
Parameter values used in the snowmelt runoff model for snow melt simulation in the central Zab basin.
|Recession coefficient (k)||x=1.121, y=0.018|
|Time lag||18 h|
|Critical temperature||0.65°C to 4°C|
|Runoff coefficient for snow||0.5 to 0.83|
|Runoff coefficient for rain||0.5 to 0.96|
|Degree-day factor||0.4 to 0.65|
|Rainfall contributing area||0 (Feb to May)|
|Reference elevation||980 m|
|Temperature lapse rate||0.5 °C/100 m|
|Precipitation lapse rate||6.5 mm|
|Initial discharge||43 m3/s|
Determination of the SRM model parameters for the central Zab basin is related to the basin size to some extent. Usually, it seems more complex in a large basin covering a variety of climatic zones than in a small basin with relatively simple features. Because the quality of the input data is directly related to the model accuracy and efficiency, the calculation and determination of variables and parameters is most important to SRM. Hence, according to previous studies, the SRM research worldwide has been mainly concentrated on the acquisition of variables, optimization of parameters, and proper analyses of the hydrological and physical characteristics of a basin. All the variables and parameters have to be specified for model input on a daily step. The SRM calculates the daily average stream flow on the ()’th day by the addition of snowmelt and precipitation contributing to runoff and discharge on the preceding day.
Validation of the SRM model
Generally, the results of the model are verified by comparison between the computed and measured discharge. The computed and measured discharges were plotted and the graphs matched well with each other. Apart from the visual observation, the model accuracy can also be checked by statistical parameters. The SRM model used two other accuracy criteria based on the statistical coefficients: one is the coefficient of determination () and the second is the volume difference (). The equations for and are as follows Eq. (13):4814):
By calculating and entering these input variables and parameters, a snowmelt runoff hydrological model (SRM) was generated.63 The comparison between the observed discharge and that calculated by the SRM model is shown in Fig. 11.
The statistical criterion for evaluation of the efficiency of fit between the observed and calculated values was applied. The coefficient of determination () for simulation was 0.8953, which shows a good correlation between the observed and calculated values. Another index was the ratio. The value for the simulation run was 0.1498%. The average measured runoff and average computed runoff were obtained as 23.156 and , respectively (see Fig. 11). The error rate for snowmelt runoff modeling using the mean comparison test (t-test) was calculated as 0.384. According to the values obtained from the t-test statistical profile, the measured flow rate showed no significant difference at the 93% confidence interval. These results suggest that the SRM can successfully simulate both diurnal water discharge and SWEs in the mountainous study area.
The use of satellite images to monitor natural phenomena such as snow cover is recommended due to considerable reductions in costs and time. Continuous monitoring of snow requires sensors with a short repeat period and with precise ability for ground separation and good color separation. The results of the present study show that the MODIS-8 daily images are appropriate for extract snow cover and use it for snowmelt estimation. The reason for using 8-day snow cover maps instead of daily ones was the fact that there was no ground truth data that might have increased the accuracy of the daily image. However, 8-day snow cover maps have the possibility that cloud that contaminated grid cells might have been determined as cloud-free and incorrectly classified. The forecasting studies show that the use of the MODIS snow maps appear to be expedient for prediction of the magnitude and timing of snowmelt runoff that contributes large volumes of flows to reservoirs.
The output hydrograph indicates that the major water resources of the central Zab basin are supplied from snowfall. The delayed effect of snow on the basin flow is also seen in the output hydrographs that show the necessity of managing water resources in this basin to reduce the damage caused by snowmelt flooding in early melting season and during spring with rain fall in the warm days of this season where the temperature of raindrops is high. Thus, quick and accurate information about the snow coverage is very important for operational forecasting and simulating water resource applications. For simulation of the runoff resulting from snowmelt in the snow gauge stations that are lacking in information about the snow coverage, the SRM model provides SCAs in the regions by providing these data from satellite images. Considering that snow gauge data collection is not done in most basins in Iran, this method is considered as an excellent alternative. The efficiency of the model’s applications in several mountainous subcatchments of the basin showed that the model performed accurately in snowmelt estimation of most basins.
First and foremost, the first author would like to thank the Universiti Teknologi Malaysia (UTM) and University of Kurdistan, Iran, for providing the facilities for this investigation. The author also would like to thank the two anonymous reviewers for the constructive comments which significantly improved the quality of this paper. The authors are also grateful to the Editor-in-Chief Dr. Wei Gao and Associate Editor Dr. Huadong Guo for their helpful comments and suggestions that enabled us to improve the paper.
Himan Shahabi is currently a PhD student of remote sensing (natural hazards) at the Universiti Teknologi Malaysia (UTM). He received his MSc in geomorphology in 2009 at the University of Tabriz, Iran. He finished his BS in physical geography (geomorphology) in 2007 at the University of Tehran, Iran. His interest research fields are environmental assessment and modeling, landslide susceptibility mapping, natural hazards and disaster management in urban planning.
Saeed Khezri is assistant professor in geomorphology in the Department of Physical Geography at the University of Kurdistan, Iran. He received a PhD degree in geomorphology in 2006 and a master’s degree in hydroclimatology in 1997 at the University of Tabriz.
Baharin Bin Ahmad is associate professor in remote sensing, faculty of Geoinformation and Real Estate in Universiti Teknologi Malaysia (UTM). He finished his PhD degree in geography in New South Wales, Australia. He received a master’s degree in surveying and mapping in Curtin, Australia.