تأثیر جهت دامنه‌های کوهستانی ایران بر تفاوت ارتفاع خط تعادل آب و یخ کواترنر

نوع مقاله : مقاله پژوهشی

نویسندگان

1 استادیار ژئومورفولوژی، دانشگاه اصفهان، ایران

2 استاد ژئومورفولوژی، دانشگاه اصفهان، ایران

3 دانشیار ژئومورفولوژی، دانشگاه تهران، ایران

4 دانشجوی دکتری ژئومورفولوژی، دانشگاه اصفهان، ایران

چکیده

 یکی از عواملی که بر ارتفاع خط تعادل آب و یخ اثر می گذارد، مقدار انرژی دریافتی از خورشید است که به عوامل متعددی، از جمله زاویه تابش، عرض جغرافیایی، مقدار و جهت شیب دامنه ها و ... بستگی دارد. در این جا کوشش شده است با استفاده از زاویه تابش خورشید و روابط فیزیکی، میانگین انرژی دریافتی روزانه در طی سال بر سطوح ارضی، در هر مقدار شیب و در هر جهت جغرافیایی محاسبه شود. در مرحله بعدی، با در نظر گرفتن دمای ثابت صفر درجه سانتی گراد (برای برآورد برف مرز دایمی) و دمای پنج درجه (در برآورد ارتفاع خط تعادل آب و یخ) و مقایسه ( comparative ) مقدار انرژی دریافتی سطوح ارضی که شیب یکسان، ولی جهت جغرافیایی مخالف هم دارند، تفاوت مقدار انرژی دریافتی از خورشید آنها برآورد گردید و سپس مقدار تفاوت انرژی به تفاوت ارتفاعی تبدیل شد؛ و به این ترتیب، تفاوت ارتفاع خط تعادل آب و یخ در سطوح ارضی مختف به دست آمد. هر چند اعداد برآوردی ارتفاع برف مرز دایمی و خط تعادل آب و یخ، رقوم ارتفاعی هستند که با ارتفاع کنونی بسیاری از ناهمواریهای ایران منطبق نیستند، ولی با توجه به پایین بدن ارتفاع برف مرز دایمی و خط تعادل آب و یخ دوره­های یخچالی کواترنری، می توان نتایج حاصله را به گذشته نیز تعمیم داد. از بررسی به عمل آمده، می توان چنین نتیجه گرفت: 1- در صورتی که افت محیطی دما، در هر صد متر افزایش ارتفاع 48 صدم سانتی گراد باشد، با کاهش 5 درجه تا 12 درجه سانتی گراد دمای ایران زمین در دوره کواترنری، ارتفاع برف مرز دایمی 1042 تا 2500 متر پایین تر از ارتفاع برآوردی بوده است؛ 2- بیشترین تأثیر جهت شیب سطوح ارضی بر تفاوت ارتفاعی خط تعادل آب و یخ، در دامنه هایی دیده می شود که امتداد کوهستانی آنها، غربی- شرقی باشد؛ 3- امتداد شمالی- جنوبی ناهمواریها، هیچ اثری بر تفاوت ارتفاع خط تعادل آب و یخ ندارد.   

کلیدواژه‌ها


عنوان مقاله [English]

The impact of mountainous skirts direction of Iran on differences in altitude of water and ice equilibrium line of quaternary

نویسندگان [English]

  • M. Moayeri 1
  • M. H. Ramesht 2
  • A. Saif 1
  • M. Yamani 3
  • Gh. H. Jafari 4
1 Assistant Professor of Geomorphology, University of Isfahan, Iran
2 Professor of Geomorphology, University of Isfahan, Iran
3 Associate Professor of Geomorphology, University of Tehran, Iran
4 Ph. D Student of geomorphology, University of Isfahan, Iran
چکیده [English]

One of the factors influencing the altitude ofwater and ice equilibrium line is the amount of the energy received from sun which depends on the radiation angles of sun and the slope. Thus, as the surfaces slopes vary, the absorption of energy will also change. Initially, this article makes an attempt to measure the radiation angle and the energy absorbed daily by the surface of land with respect to different slopes and directions. Then, the differences between the amounts of energy absorbed by two different surfaces varying in terms of slope and direction will be estimated, considering a constant temperature. Consequently, the energy differences will be converted to altitude differences, using the statistical relations. This paves the way for the estimation of the differences in altitude of water and ice equilibrium line. Although this article has used a hypothetical altitude, the resulting finding can also be used in the context of the past ages because the altitude of water and ice equilibrium line was low at the ice age of quaternary. This study had the following findings; 1- Averaged reduction of temperature of environment is .48 per 100 meters increasing in height with reduction of temperature in Iran about 5-12 in quaternary period height of snow is downer about 1042-2500 then current circumstance. 2- The highest impact of the land surfaces directions on the altitude of ice and water equilibrium line can be seen in the skirts whose mountainous ridges expansion is west and east. 3- The north-south expansion of mountainous ridges doesn’t have any effect on the difference in altitude of ice and water equilibrium line.

کلیدواژه‌ها [English]

  • altitude
  • water and ice equilibrium
  • quaternary
  • energy
  • sun radiation
  • slope

What is problem?

Today, Iran (and other countries) wants to expand his economic by tourism. Rural tourism is very important for economic development because it have important role on rural development and therefore planners and experts focus to it. Iran has many places with beautiful landscape (natural and building) and different culture that are attractive for tourists. These help to economic development especially for villages and decrease problems that are in villages (example out-immigration, poverty…). This paper study tourism effects on local people.

The rural tourism has economic, social and environmental dimension and tourism experts focus on it because the tourism has many profit for villages. But there are a questions; what groups use from tourism profits? Do local people use from tourism profit or none local people use from tourism profit? This is the important question that this paper investigates about it.

The rural tourism is a kind of tourism and has many profits (for example decreasing poverty or making job) for remote villages, especially villages that have wealthy culture (WTO: 2005).

 

Introduction

In this article has been attempted by help of SPSS and surfer soft ware's and implement of densitizing and kriging to present a realistic analysis on structure of differences between the altitude of water and ice equilibrium line of quaternary in Iran. The level of energy role, temperature and geographical latitude can be specified on the regional structure of Iran latitude of water and ice equilibrium line. Such analysis are very useful, for example it shows that by increasing 1 degree in environment temperature, the permanent height of snow rise up to 208 m (in case of decrease in environment temperature to 4.8 leads to increase height up to 1000 meters) It makes possible to estimate isotherm line in every place of country at any slope and direction. The Iran altitude of water and ice equilibrium increase from was to east due to influence of thermal, cooling and humid cells position and from north to south because of the geographical latitude effects. A one degree increase in geographical latitude in slope of 13º will decline the height of permanent snow from 140m (Northern skirts) to 204 meters (eastern north skirts).

 

The historical Back ground of the topic

The words, glacier equilibrium line, water and ice equilibrium line and permanent height of snow are common in the glacial literature. In glacial equilibrium line, the level accumulation and ice melting are equal meaning that in a glacier two sub system can be explored. One is the ingredient accumulation section including part of glacier that more ingredient accumulated during a year and have less melting thus accumulation exceeds fraction (abrasion) the other is abrasion section which includes inverse effects, meaning that abrasion is much more compare to accumulation level, these two section borders called as glacial equilibrium line (Motamed, 2000). The water and ice equilibrium line is the line between marginal limit of glacier`s latch and level of running water caused by ice melting. Glacier would move due to various factors such as slope gravity and so on. This movement would be possible till glacier melting equal to its movements (Mahmoudi1990).

This border is known as water and ice equilibrium line. For the first the term "time permanent height of snow" had been used by Libutery (Bellreopomlor 1990) and then the term have been defined by (Sehcheglova and klebelsberg) quoted by Schweitzer 1970 as: the limit on horizontal surfaces and upper which for whole year is covered by snow. In Bellreopomlor opinion permanent height of snow was named glacial equilibrium line. Latu (Libutery and Miserly) quoted by (Schweitzer 1970, king 1970, Hawkins 1989, Humlum 1986, Ritter et al. 1993) made difference between permanent height of snow and glacier equilibrium line and have believed that could not be equal because they have various height usually permanent height of snow is higher than glacier equilibrium, in (Hawkins) opinion 1985. There is several meaning for ecological snow line which has general and limited application; this term in spite of its name refers to a level and surface, where equilibrium line is precise term with more application. In glacier equilibrium line separates ice accumulation region from ice melting one (Selby 1986, Humlum 1988). In geomorphologic view permanent height of snow is important in this respect which indicates the lowest level of glacier climate (Humlum 1998). The first condition for glacial formation is presence of permanent height of snow on the earth (0 degree isotherm) the permanent height of snow is a factor that shows relationship between glacier and climate and the dominant glacial climate on the permanent height of snow support the glacier durability, in case height of permanent line of snow exceed some more leads to mountainous glacier regression and destruction of ice cap (calyptras).

Each glacier at us higher glacier equilibrium line, is able to romove gravel from itself, in some cases, the melting snow water from upper section of equilibrium line. In the melting section influencing by freezing weather transformed to complete ice. These conditions providing the glacier aliment from lower parts (Dallaloghli 2003). There glacial latch have been element from high mountainous regions have capability to come down to foot of mountain and impose their effect to whole valley. The extent of progression can be recognized from geomorphic witnesses, almost all the height values given by various researchers (Scholars) in case of moraine observation and glaciers effects in the height less than 3000m indicating the height of ice and water equilibrium line in Iran. In Touchal northern skirts, the moraine came down up to 2200m till Shahrestanak valley bottom (Aivazi 1995) (Hagedorn et al. 1974) in the Shirkooh-Yazd indicating the existence of glaeial moraines at the height 1800 to 2800m that showing 2 glacial phases (Mahmoudi1887). Right has estimated permanent height of snow lowering up to 1800m from sea level in Kurdistan Mountains. He had observed moraines in small valley high 2600m in west south of Aznaa- Oshtorankooh (Aivazi 1995). In the Gorgan Rood valley glacial latch came down to the height of 1200m and its moraines obstructed Gorgan Rood main valley which made a crag of more than 10 meters height (Zomorrodian 2002, Pedrami 1982) refers to a glacial marks in Kashan which had been protected its lateral moraines very well. According to him glacial latchs of Sardab Rood and Chaloos River in old turgidity had been lowered up to those days Caspian water level.

He talks about small cirques with the height of 1700 m in Rood Barak village west side as well. There are about 50 large and small known cirques that most of them placed at 1800m height and more. Regarding one cirque has come down to 1600m height (Mahmoudi 1997). Ramesht (2002) is one of the researchers in glaciology of central Iran. Approved the findings of scientists like Hagedorn and Kuhle. He studied the glaciers of 1600m height and more with relience on geomorphologic observation and abrasion effects. In several point of central Iran such as Zefreh- Shirkooh emphasized on descending of ice letches till plain lands. Obviously the receiving energy from sun play vital role in lowering water and ice equilibrium line, which is affected by direction and slope of land surfaces. In this case few studies are available tough. Ramesht 2002, shoshtari zadeh 2003 have referred to Zefreh and Salafchegan ice and water equilibrium line difference of 200m very shortly.

 

Methodology of Data

To specify permanent height of snow or height of water and ice equilibrium line according to given isotherms by the researchers. 0 degree for permanent height of snow and 5 degree for ice of water equilibrium line are required to specify temperature difference in sloped surfaces. Temperature itself is result of various factors such as receiving energy from sun which is depends on solar radiation angle and its constant. So for this purpose we have to estimate the receiving energy from sun in various slopped surfaces. The degree of energy receiving from the sun is affected by quantity and slope direction, geographical latitude and etc. for estimation of annual radiation angle average in constant slopped surfaces. The equation I has been used (Joe Michal sky 1988) and along with this equation. The following step carried out:

1. Medium month day was used for estimation of monthly radiation angle average.

2. To specify the monthly radiation angle average in medium day of month, every 4 minute. Radiation angle was calculated and by their summation the average has been calculated.

3. Monthly radiation angle average was estimated for slope of 0 up to 90º degree one by one. (Second stop was repeated 90 times)

4. For slope of horizontal surfaces in 360º degree scale may be infinite variable. In this article 16 directions with 12.5 degree distance have been selected as follows:

North, North- northern east, northern east, east- northern east, east, east- Southern east, southern east, south- southern east, south, south- southern west, southern west, west- southern west, west- Northern west, northern west and North- northern west.

The monthly radiation angle average have been calculated for energy direction and slope degree (second step was repeated 16 times)

Equation I

 

 

 

 

Solar radiation angle: the angle between solar direct radiation line and the vertical direction on horizon surface which is sun height angle supplementary.

Slope angle: the angle between concerned plane surface and horizon surface.

: Solar dip angle for any special day of year which approximate value calculated by equation II (Stickler, 1999)

Equation II 

In equation II n represents the number of day in a year.   Calculation for all the days is very tough wok, so as referred, month median day have been used means that the day which solar dip angle in atmosphere outward surface is proximity value of monthly mean (Azad et al 1998).

 Plane direction angle: perpendicularity on plane image deviation with local meridian of longitude. In northern hemispheres where =0 plane direction has been towards south completely and when plane direction is towards eastern south,  has negative (-) value angle plane towards western south  takes positive value. This angle varies up to ±180.

: Solar hour angle, which conventionally is negative value for pre meridian and positive for post meridian and can be achieved by the equation III.

Equation III       =(X-12)*15

 

X considered as hour and indeed constant 12 is the solar meridian. Basis of 15 is the 15º earth relation in an hour, in order to different length of day during the year. For precise calculation X assigned as a variable with maximum range of 5 to 20 for June and July and minimum range of 7 to 17 for December and January. For Cos estimation in equation I digits accordingly between 105- 135 varying by 4 digit interval and substituting  (earth radiation is one degree every 4 minute this variation provide the radiation angle for energy 4 minute). In case region excludes the solar radiation solar height become less than zero (negative) and the estimated quantities do not contain the data and have not any effect on monthly solar radiation angle.  The reason for this calculation is estimation of annual radiation angle in sloped surfaces with different direction. For horizontal plane  value is zero and equation I get simplified to equation IV.

Equation IV       

For estimation of total receiving energy during each unit of time equation V is used (Alizadeh 2006)

Equation V    I=W.  

I: Amount of receiving energy per minute

W: Solar constant at point of space which earth atmosphere is formed (started) its quantity is  per minute.

: Sun height- the angle between sun direct radiation side and horizon direction and  supplementary (radiation angle) so Cos  is equal to Sin  and in equation VI sin can be substituted by Cos.

5. Annual average of radiation in every direction and every slope degree have been specified by repeating the above mentioned steps for 12 months.

6. Annual average of radiation angle for 16 geographical latitude in Iran (25º to 40º) calculated by above steps as well.

7. Value of solar constant on the earth surface is affected by many factors such as geographical latitude, sun radiation duration (time), variegation rate, relative humidity, water vapor pressure, regularity of weather and etc. For solar constant value estimation on earth surface, given values of energy on horizontal balance of radiometer networks were used for higher confidence coefficient.

Data of ecological station equipped with radiometer during index period (1966-91) and including temperate, total radiation on horizon surfaces height and geographical altitude and latitude of the stations. The quantity of energy received by station has been estimated by Khalili (1997) and temperature average and others required information have been collected from metrology organization site. Spatial variation of receiving radiation sun are negligible in compare with some of meteorological factors such as rainfall, except micro ecological conditions which may influence visibility at point because of special cloud coverage situation or related phenomenas. Every station can represent the radiation condition of a region up to 100 km of distance and more (Azad et al. 1972) considering this and 16 radiometry station in country and by illustrating perpendicular line between a bistation networks. Iran was divided into 16 parts (Khalili 1997). Regarding that in the Quaternary period the north region was influenced by Atlantic ocean and Meditrean and Caspian sea (Tahoni 2004) and further more because of 2 different existing climatically region in north Iran, for accurate examination of quantity of received energy in north Alborz skirts, Tehran station information's were used. Although in estimation, Ramsar station data in Caspian littoral regions are more general and its data can be used for receiving energy station in northern skirts (Khalili 1997). Considering lack of high peaks examination of Persian Gulf and Oman sea littoral regions had been ignored (plan I)

8. Then average receiving energy from sun on earth surface horizontal balance during year (per minutes a day) dividing on Cos (annual radiation angle of horizon surface).

9. By dividing outcome value of previous step on total numbers of minutes a day (24 hour- 1440 min) we obtain the value of annual average of solar constant for any station on the earth.

Regarding to these discussions, in 32º geographical latitude (Yazd) the annual average of receiving energy from sun on horizontal balance of earth surface is 542/3  per day (Budget, planning organization) and annual average of sun radiation angle (Cos 0.514708) has been estimated for 59.2º so by the equation V have:

542.3=W*0.514708

Day

Min    

 

 

 

Figure 1: Iran explanting regarding to receiving energy from sun

 

 

Means annual average of solar constant for (altitude 32ºof horizon surface is 0.732  per minute.

10- quantity of receiving energy for every slope and direction of slopped earth surfaces can be arrived from equation 6.

Equation6   

11. for energy transformation into temperature, first between the monthly average of receiving energy (per day) and monthly average of temperature () in every station linear equation established (table 1) then annual average of receiving energy calculated in previous step with help of these equations, had replaced with temperature average.

12. for estimating average of earth surface slope by using the value of slope situation in Iran according to 109 page of Iran slope map (Aivazi 1995) and equation 7, the regions which have possibility of transforming and passing of Quaternary glacier have been carried out.

 

 

Table 1: Temperature and energy equation in selected stations

Correlation co efficient

temperature of energy equation in each station

Station

Correlation co efficient

temperature of energy equation in each station

Station

0.916

0.954

0.930

0.929

0.902

0.975

0.928

T= -9.317+0.0577E

T= -5.796+0.063E

T= -9.315+0.0465E

T=-7.59+0.0613E

T= -11.918+0.058E

T= -7.341+0.0483E

T= -9.309+0.0464E

Shiraz

Tabas

Tabriz

Tehran

Zanjan

Yazd

Oromieh

0.927

0.899

0.955

0.933

0.923

0.926

0.967

0.956

T= -1.933+0.0647E

T= -10.458+0.0623E

T= -10.142+0.058E

T= -8.533+0.0597E

T= -10.525+0.0595E

T=-12.326+0.063E

T= -11.57+0.064E

T=-7.822+0.060E

Ahvaz

Bakhtaran

Birjand

Bojnord

Esfahan

Hamedan

Kerman

Mashhad

 

 

Equation 7           

S: Average of slope for assigned surface.

a: Area percentage of each slope parts to percentage area of total assigned surfaces.

s: slope percentage of each port of surfaces.

A: sum of total percentage of assigned area.

For example average slope of skirts surface or mountainous or rocky surfaces are calculating as follows:

 

And average of interior plateau slope, plateaus and littoral plateaus trace on slope surfaces of diluvia plateaus, glassy, pediment, wild skirts of plateaus, margins and abrasive surface:

 

 

For example at 1230 m height (Yazd) earth surface with slope to eastern north have mean temperature of 16.177 and to southern most mean temperature will be 4. Now regarding to environment temperature tall 4.8 per 1000m in Iran (Masoodian2003) we have:

 

So actual permanent height of snow in northern east skirts is at 4600m height and southern west skirts will be at 5608 m. Estimated permanent height of snow in northern east is very close to the permanent height of snow given by Grouter (1978) for Shirkooh (4600-4700m). So we can generalize directions of northern west- south east with inffluence of earth surface with 13º slope towards northern east (shady skirts) and southern west (sun exposed) for central Iran Mountains. So likewise for the rest of Iran regions after specifying direction and affluent earth surface slope for receiving solar energy, the permanent height of snow and height of water and ice equilibrium line had been calculated, 28 height peak of country were specified by examination of affluent earth surface direction on quantity of receiving energy which affluent direction are as follows: 19 cases the northern east and southern west skits (67.86%) double peak, west- northern west and east- northern east (7.14%) five- peaked, east- northern east and west- southern west (17.86%) and double peak, eastern and western skirts (7.14%). In any part of Iran the affluent surface in permanent height of snow cannot be in south or north direction. By comparing 2 station (Ahvaz of Tabriz) with almost 6º geographical latitude difference and same geographical altitude (48.5º) locations what is resulted is the permanent height of snow with different direction by increasing geographical latitude have intangible variation which can be ignored.

In other words, difference in snow line has least relativity with geographical latitude and at every geographical latitudes this different in most sun exposed or shaded areas is about 1500m. This value is more than 500m for east- southern east and west- northern west skirts and in west or west direction is zero.

Variation average of snow line height between shady and sun exposed skirts is around 1435m, and average of Iran snow line height is 4170m, and its median is 4353m. The table Π shows the water a ice equilibrium line height difference in Yazd.

 

 

Table 2: water a ice equilibrium line height difference in Yazd

EorW

ESEorWNW

WSWorENE

SEorNW

SWorNE

SSEorNNW

SSWorNNE

N  or S

Direction

City

0

540

1008

1429

1608

Yazd

 


Conclusion

The most outstanding direction of earth surface slope of Iran rough skirts are towards northern east and southern west. Such directions are spread from 28º to38º northern latitude 46º to 62º eastern altitude. The permanent height of snow at 13º slopped surface of southern west skirts of Shirkooh- Yazd is 5608m, and for northern east skirts is 4600m and their permanent height of snow difference is more than 1008m, in case of 8.8 environment temperature fall in Pleistocene period the permanent height of snow had been 1834 m less than actual permanent height of snow, means that permanent height of snow for shaded skirts would be 2760 m, its water and ice equilibrium line in case of accordance with Quaternary 5 temperature, had been 1042 m less than permanent height of snow at 1725m height. Some parts of Tallish Mountains (Bagharodagh peak. Iran- Turkey and Iran- Iraq Boundary Mountains (Ghandil peak) have northern- southern expansion. This expansion examination in different geographical latitude representing that in earth surface skirts such expansions towards east and west have no effect on quantity of solar receiving energy and do not cause any different in height of water and ice equilibrium line the slope of earth surface skirt of Peru peak in Kurdistan is towards north- northern east and south- southern west. Such surface can lead to 1436m height difference in water and ice equilibrium line.

The affluent direction a Sabalan peak skirts is west- northern west and east- south east, and for Sahand, Hezar, Taftan, and Shahvar (Kurdistan and Kermanshah border) mountain skirts are toward east- northern east and west- southern west, the earth surface of those two group function convectively regarding to receiving energy and with geographical latitude proportion creating height difference between 523 to 586m in water and ice equilibrium line. In shady skirts any degree increase or decrease in earth surface slope can cause about 40m decrease or increase on height of snow line and for sun exposed skirts any degree decrease or increase in slope leads to 29m decrease or increase in height of snow line.

If assume the water and ice equilibrium line regarding maximum height difference from north to south with 13º  slope and 1608m height, and examine the assumptive height of water and ice equilibrium line of other direction based on this height the results are shown in figure I (geographical latitude 32º north and geographical altitude 54º east). As seen in this figure, height of water and ice equilibrium line of the south direction. Coinciding to 1608m by reducing disposition of earth surface slope direction to south, water and ice equilibrium line height decreased gradually till that become zero for northern slopped earth surfaces meaning that permanent height of snow as well as water and ice equilibrium line reach to the maximum in south and minimum in north.

Figure 2: Assumptive height of water and ice equilibrium line in different direction (Yazd)

Regarding to quantity of solar receiving energy if we divide the direction of earth surface slope of Iran rough skirts into 2 eastern and western halves for both halves by increasing distance from south, water and ice equilibrium line will reduced so because of less received energy water and ice equilibrium line place in lower height and approach us median in west or east direction. Afterward the height of water and ice equilibrium line reduces till that reaches to its minimum in north. Therefore, two eastern and western halves regarding energy receiving function convectively to each other. As whale water and ice equilibrium line in east to south skirt compare to east to north skirts is higher due to more solar receiving energy the same situation can be seen in west to south skirts as well as west to north skirts.

 

  

 

Figure 2: Assumptive height of water and ice equilibrium line in different direction (Yazd)

 

 
-                      Alaee Taleghani, Mahmoud, Geomorphology of Iran, Tehran, Gomes, 2003.
-                      Alizadeh, Amin, Principles of Applied Hydrology, Ferdowsi University Mashhad, 2006.
-                      American Society for Testing Methods (ASTM) Standards (E-891) and (E-892).
-                      Azad,E, et al Solar energy, Tehran university publication, 1990.
-                      Dallaloghli, Ali, Glacial equilibrium line altitude in the Sabalan Mountain
-                      Dallaloghli, Ali, The researchers in morphological of systems effect in the Sabalan mountain north skirt and form sedimentary plan Meshkin shahr, Tabriz University 2002.
-                      Duffy, J. A. and Beckman, W. A. 1991. Solar Engineering of Thermal Processes, 2nd edn. J. Wiley and Sons, New York.
-                      http://solardat.uoregon.edu/SolarRadiationBasics.html
-                      Jedari Avazi, Jamshid, Geomorphology of Iran, P.N university,1995.
-                      Joe Michal sky, J. J. 1988. The astronomical almanac's algorithm for approximate solar position  (1950-2050), Solar Energy, 40, 227-235.
-                      Khalili, ali and Rezai-e sadr, Hassan, Estimation of Global solar Radiations over Iran Based on climatically Data, Geographical Research, Tome 46,PP.15-35,1997.
-                      Mahmoudi, Farajollah, the evolution of reliefs pendent the Quaternary Iran,  the journal geography researchers Tehran university, number 23, Pp: 5-043, 1988.
-                      Mahmoudi, Farajollah, Geomorphology structural et dynamical extern, Tehran university publication, first Edition, 1989.
-                      Masoodian, S.A, Aspatial analysis of monthly temperature in Iran, Journal of research university of Esfahan, Humanities,Vol 15,No1,2,PP87-96,2003.
-                      Pedrami, Manochehr, The Quaternary absolution old, societies collage journal, Tehran University, volume 17, no3& 4, 1988.   
-                      Ramesht, M.H. 2002, glacial evidence of the Zefreh, proposal number 800305, approved / June/ 2002, Esfahan University.
-                      Richard J Chorley. Geomorphology volume 4, Translated by Ahmad motamed, Samt publication, 2000.
-                      Shoshtari zadeh, Najmeh,  Ice Cape Central Iran ( Sarafchegan Region ) “ M.A. “ Thesis on Geomorphology Islamic Azad University, Najaf Abad Branch
-                      Tahooni, Poran, geography researchers, the geomorphologic instances the Pleistocene glacier erosion in the Talesh Mountains, number 47, Pp: 31-55, may 2004.
-                      The planning and budge organization, the standard technical and codify jobs note, climatic standards in structure, energy solar over Iran, the proof central economic-  society and publication, 1997.
-                      Write, H.E, N.Y. 1963 Preliminary Pollen Studies at Lacked Serbia, Zagros Mountains, Southwest Iran, and Science.
-                      Zomorrodian, MJ, Geomorphology of Iran structural processes &endogen dynamics, Ferdowsi University of Mashhsd,2002.