Gullies of two Hungarian regions – a case study 325 Hungarian Geographical Bulletin 60 (4) (2011) 325–342. Gullies of two Hungarian regions – a case study Gergely JAKAB1, Balázs MADARÁSZ1, Anna ŐRSI1, Zoltán SZALAI1 and Ádám KERTÉSZ1 Abstract Gully erosion plays a decisive role in removing the fertile layer of the soil and it has an important long term eff ect in relief formation. The objective of this paper is to compare gully development and distribution in two pilot areas in Hungary and to reveal the factors controlling gully formation. The pilot areas are natural micro-regions, i.e. the Börzsöny Mountains and the Zselic hills. The results presented in this paper are based on the data of a nation-wide gully cadastre recently under compilation by the authors. The cadastre contains the gullies shown on digital maps (1: 10,000), gully lengths, mean gradients, land use and main properties of parent materials and soils. The results point to the develop- ment and formation of signifi cantly diff erent gullies in the pilot areas as a consequence of diff erent environmental conditions. The most important result is the introduction of the concept of equivalent gully length characterizing the gullies of the given category. Topography is the main driving force in gully formation followed by land use type, par- ent rock and soil properties. Keywords: gully erosion, gully cadastre, Zselic, Börzsöny Introduction Soil erosion is one of the most important agents in contemporary landscape formation. This statement applies also for the subhumid regions of Central Europe. Soil erosion att acks the uppermost, fertile soil layer and the eroded soil contains valuable nutrients (Farsang, A. et al. 2011; Borcsik, Z. et al. 2011). If the eroded soil will be transported into lakes eutrophication will be accel- erated (Csathó, P. et al. 2007). Sheet erosion processes aff ect extended areas, the result is, however, a relatively slow change in topography. Gully erosion appears on relatively limited portions of the surface but it leads to the removal of a huge amount of soil and it makes rapid and remarkable change in surface topography (Pécsi, M. 1955; Kertész, Á. 2009). Both sheet and gully erosion contribute to relief formation (Jakab, G. et al. 2009). Gully erosion is a threshold 1 Geographical Research Institute Hungarian Academy of Sciences, H-1112 Budaörsi út 45. Budapest, Hungary, e-mail: jakabg@mtafk i.hu 326 phenomenon (Poesen, J. et al. 2003) but the identifi cation of the threshold value is a very complex procedure (Poesen, J. et al. 2003; Kirkby, M. and Bull, L. 2000). The occurrence of gully, or sheet erosion may change also periodically. The glacial periods are presumed to have been ruled by gully erosion while the interglacials characterized by sheet erosion (Pécsi, M. 1997). Climate is not the only environmental condition that regulates soil erosion (Kertész, Á. 2006; Smolska, E. 2007) the role of soil properties, parent material (Buzek, L. 2007), topography and land use is also important (Kertész, Á. 2008). Many authors refer to the signifi cance of land use changes (Grace, J.M. 2004; Gábris, Gy. et al. 2003; Galang, M.A. et al. 2007; Centeri, Cs. et al. 2009). The development of an already incised gully can only be stopped by radical changes and it is extremely expensive to remove it from the landscape (Kirkby, M. and Bracken, L.J. 2009) Gullies develop mainly on loose, unconsolidated sediments such as loess and loess-like deposits (Poesen, J. et al. 2003; Zglobicki, W. and Baran- Zglobicka, B. 2011) as well as on marine sediments in the Mediterranean (Poesen, J. et al. 2006), and on sandstones (Hegedűs, K. et al. 2008) etc. Gullies developed on unconsolidated volcanic rocks represent a special type (Pintér, Z. et al. 2009). An example is the badland developed on rhyolite tuff in Kazár, Hungary (Horváth, G. et al. 2010). Gully erosion has received less att ention due to the complexity and dif- fi culty of its investigation (Valentine, C. et al. 2005), nevertheless adequate gul- ly susceptibility prediction would be necessary to control soil loss (Conforti, M. et al. 2010). There is no soil type, rock type, land use or topography which could alone launch gully initiation. The interaction of the factors controlling gully erosion is needed for gully formation (Muńoz-Robles, C. et al. 2010). The aim of this study is to identify the spatial properties of gully ero- sion in two pilot areas in Hungary with the analysis of the recently compiled gully cadastre. The hypothesis to be tested is whether diff erent environmental conditions (topography, land use, parent material and soil) generate distinct spatial patt erns and distributions of gully systems. An additional goal is to identify the relationship between gully formation and land use change. The study sites In this study two pilot areas are compared, one from a mountain range and another from a hilly country, namely from the Börzsöny Mountains and from the Zselic hills, respectively. The Börzsöny Mountains are of volcanic origin, forested, with steep slopes and with high relative relief values. This area is se- lected because results of previous research are available from here (Madarász, B. 2009; Madarász, B. and Jakab, G. 2009). 327 Most of the Zselic hilly country pilot area is covered by loess and cultivated. Detailed investigations have been carried out in the Somogy hilly country for almost two decades (Jakab, G. 2008; Jakab, G. et al. 2005, 2009, 2010a) so the reason for choosing the nearby Zselic was to have a new sample area similar to the Somogy hills. The pilot areas were identifi ed according to the Inventory of Natural Micro-regions of Hungary (Dövényi, Z. ed. 2010). The Börzsöny pilot area cov- ers 447 km2 and is situated in the northern part of Hungary (Figure 1). It is a paleo-volcano formed in the Miocene between 16.5 and 13.5 Ma B.P. (Pécskay, Z. et al. 1995; Karátson, D. 2007). The older volcanic rocks of are mostly covered by younger volcanic deposits, however, at the margins of the mountain at some spots older deposits outcrop. During the Pliocene and Quaternary times ero- sion and tectonic movements were the main landscape forming agents in the Börzsöny (Láng, S. 1955). The recent surface and valley system was formed by the erosion of the nearly 800 m thick volcanic strata (Karátson, D. 2007). The el- evation of the study area varies between 120 m and 939 m a.s.l. The characteristic surface forms are debris fl ows as well as elongated and steep hillsides dissected by relatively young V-shaped valleys. Relative relief values gradually decrease towards the periphery of the mountains, from 350–370 m km-2 to 100–150 m km-2 (Madarász, B. 2009). Most of the slopes are exposed to west–north-west and to east–south-east. The study area belongs to the cool and wet climate do- main. Mean annual temperature does not exceed 8–8.5 °C with only 7–8 °C in the highest region. Annual precipitation varies between 600 and 800 mm with higher values in the peak region (Dövényi, Z. ed. 2010). The Zselic hills (1170 km2) are situated in the southern part of the Transdanubian Hills macro-region (Figure 1). It is a relatively plane area with gentle slopes similar to a pediment (Sebe, K. et al. 2008). Almost the whole area is covered by loess, mainly of Wurm origin, and can be classifi ed mostly as slope loess (Kapronczay, J. 1965). At some deeper lying spots sandy Pannonian sediments occur on the surface (Lóczy, D. and Gyenizse, P. 2003). Forest man- agement was the traditional main occupation and income source in the area (Gyenizse, P. et al. 2008) supported by a rather high amount of annual precipi- tation (700–750 mm, Kapronczay, J. 1966). With the expansion of arable fi elds soil erosion became an increasingly serious problem aft er Second World War (Lóczy, D. and Gyenizse, P. 2003). The area can be characterized as fragmented from geological and social aspects. Methods The gullies of the study areas were digitized using the 1:10,000 scale maps of the Unifi ed National Map System (EOTR). A GIS was organized which con- 328 Fi g. 1 . L oc at io n of th e st ud y si te s w ith th e su rv ey ed g ul lie s 329 sisted of the digital map of the gullies, the map of soil properties, the land use map and the topographic map of the pilot areas. The map of soil properties is a map series derived from the AGROTOPO map at a scale of 1:100,000 (RIS- SAC 1991). This database contains information on soil type, parent material, texture, clay mineral composition, hydrology, pH, soil organic matt er (SOM) and soil depth. The land use map was compiled on the basis of the CORINE Land Cover database (1:50,000) (Büttner, G. et al. 2002). Gully sections with intermitt ent or permanent water fl ows were excluded. The original CORINE database has more than 60 land use types which had to be generalized for the purposes of the pilot study. As a result six land use types were created focusing on soil conservation, ranked according to their endangerment by soil erosion, i.e. cultivated area including arable land, vineyard and orchard followed by permanent vegetation, forest and wetland/ scrubland. Elevation data are taken from SRTM (Shutt le Radar Topography Mission, Rabus, B. et al. 2003). The resolution is 90 m. The soil gradient map with fi ve category classes (see e.g. Pécsi, M. 1991) was derived from the eleva- tion model. Soil, land use and topographic data are given for each gully. If a gully extends over two or more pixels then the gully is identifi ed with the value/ property occupying the largest area in the gully. Concerning numerical data the gully is identifi ed with the arithmetical mean of the values of the pixels in question. Results and discussion The distribution of soil parent material and soil type is shown in Table 1. The Börzsöny is dominated by volcanic rocks mainly covered by Luvisol. On the remaining area Cambisol and Phaeozem are to be found. Phaeozem stands for a special type of soil developed on volcanic rocks named “erubáz soil” in Table 1. Soil type and parent material distribution Area Soil type (%) Parent material (%) Ph ae oz em Lu vi so l C am bi so l C he rn oz em G le ys ol G la ci al a nd al lu vi al d ep os its Lo es s Te rt ia ry a nd ol de r d ep os its A nd es ite , b as al t, rh yo lit e Börzsöny Zselic 4 – 81 63 15 19 0 5 0 13 0 13 1 87 33 1 66 – 330 Fi g. 2 . M ai n so il ty pe s of th e Zs el ic (l eft ) an d Bö rz sö ny (r ig ht ) s tu dy a re as 331 the Hungarian nomenclature, see Madarász, B. 2009. In the valley bott oms and near the Danube Fluvisol occurs (Figure 2). In the Zselic area loess is the dominant parent material. Texture and organic matt er distribution data are shown in Table 2, soil depth in Table 3. The Phaeozems of the Börzsöny site are shallow (soil depth < 40 cm), other soil types are thicker, up to 70 cm (Table 3). Presumably, deep gullies cannot be formed on volcanic rocks covered by shallow soils. There are gul- lies developed on this soil types, too. During the fi eld campaigns there were detected deep gullies with steep slopes cut into the hard volcanic rock. These gullies can be older, developed presumably during the Pleistocene. There are also shallow gullies running parallel with each other and brought about by man induced activities (e.g. roads used for wood transport or for military training). The central part of the Börzsöny was used by the army as a training area in the second half of the last century. These shallow gullies are not shown in the digital maps. The Zselic pilot area has deep soils. Their high proportion (95%) can be explained by map generalization, the data are taken from the AGROTOPO map series (1:100,000). Small spots of Leptosol, however, could be observed all over in the area. The deep soils are mainly colluvia accumulated by erosion. Due to the porous parent rock it is hard to determine the boundary between the rock and the soil, especially in case of arable fi elds. The two study sites have distinct land use structure (Table 4 and Figure 3). The Börzsöny is almost completely covered by forest. The Zselic has large arable fi elds and the forest spots are fragmented. The average land area of the individual land use categories show similarity in the two pilot Table 2. Texture and organic matt er (OM) distribution Area Texture (%) OM content distribution (%) Loam Clay loam 50–100 t ha-1 100–200 t ha-1 200–300 t ha-1 300–400 t ha-1 400 < t ha-1 Börzsöny 23 77 0 96 4 – – Zselic 100 0 63 24 1 7 5 Table 3. Soil depth distribution Area Soil depth (%, cm) 20–40 40–70 70–100 100 < Börzsöny 4 62 – 34 Zselic 0 – 5 95 332 Fi g. 3 . S im pl ifi ed la nd u se m ap o f t he Z se lic (l eft ) an d Bö rz sö ny (r ig ht ) s tu dy a re as 333 Table 4. Land use distribution of the study sites derived form the CLC50 (2000) database Land use Zselic Börzsöny Total area (ha) Rate (%) Average area (ha) Total area (ha) Rate (%) Average area (ha) Artifi cial surface Cultivated area Forest Pasture Swamp Scrubland 3,585 46,524 43,688 19,076 340 3,304 3 40 37 16 0 3 50 358 299 64 43 49 985 4,235 33,928 4,134 1 1,447 2 9 76 9 0 3 47 132 1,696 63 0 39 areas except for cultivated area and forest. Artifi cial areas include built-up areas, open-cast mines etc. The distinct land use structure of the sites can be explained by diff erent parent materials and relief conditions. In Börzsöny the largest homogeneous areas are covered by forest and they fall into the highest slope gradient class. In Zselic the steepest slopes have a very limited spatial extension, they are covered by forest, but the woodland is fragmented (fi gures 4 and 5). Gully properties are summarized in Table 5. The gully dissection index value is almost the same in the two pilot areas. The average values are three times larger than the lower limit value (0.5) of the highest category of the Hungarian classifi cation underlying the necessity of introducing an additional gully dissection category as already suggested (Jakab, G. et al. 2010b). Average gully length is slightly higher in the Börzsöny. The diff erence between the average and median values point to the anomalous distribution of the data similar to the case of the Tetves catchment (Jakab, G. et al. 2005). The diff erence between the median values is also very small. The minimum gully length values (2 m, see Table 5) do not refer to real gully lengths as they represent only parts of gullies which are longer but they are cut by the border of the natural micro-regions. Comparing maximum gully lengths of the pilot areas the value of Zselic is twice as large (16 km) as that of the Börzsöny. The explanation is the porous parent material and smaller gradient values in the Zselic area. Comparing tables 6 and 7 a conspicuous ob- servation can be made, i.e. looking at the properties (e.g. soil type, land use etc.) there are only small diff er- ences between the percent- age values calculated from the number of gullies, com- Table 5. Main gully properties Indicator Börzsöny Zselic Number of gullies Dissection index (km km-2) Total length (m) Average length (m) Minimum length (m) Maximum length (m) Median (m) 2,260 1.43 638,309 282 2 7,308 126 6,579 1.45 1,693,374 258 2 16,220 133 334 Fi g. 4 . S lo pe g ra di en t d is tr ib ut io ns o f t he s tu dy a re as d er iv ed fr om th e SR TM d at ab as e (R ab us , B . e t a l. 20 03 ) 335 Fig. 5. Slope gradient classes of the Zselic and the Börzsöny study sites Table 6. Average properties of the surveyed gullies in the Börzsöny Property Rate according to the number of gullies (%) Rate according to gully length (%) Average gully length (m) Soil type Phaeozem Luvisol Cambisol Chernozem 1.1 86.9 11.5 0.4 0 88 11 1 118 283 279 808 Parent rock Loess Tertiary and older deposits Volcanic rocks 1.9 36.4 61.7 3 31.0 66.0 441 241 302 Soil texture Sandy loam Loam Clay loam 0 30 70 1 26 73 110 242 297 OM content t ha-1 100–200 200–300 98.9 1.1 100 0 284 118 Soil depth (cm) 20–40 40–70 > 100 1 61 38 0 66 34 118 305 251 336 pared to those from gully length, e.g. the percentage of gullies on Luvisol is 86.9% versus 88.0%. Analyzing the distribution of gullies according to soil types of the Börzsöny the fi rst interesting phenomenon is that there are relatively few gul- lies on Phaeozem. The reason for this is the shallow soil on solid parent rock as mentioned above. Another reason may be the geomorphological position of these soils, i.e. they develop on the highest parts of the mountain, at the caldera fringe which is not a favourable location for gully formation because of the lack of suffi cient catchment area. It is striking that the percentage of Luvisol is slightly higher than the percentage of its area. Compared with Cambisol the share of the latt er is less than its territorial percentage. Luvisol is more resist- ant to gully erosion than Cambisol (Nachtergaele, J. and Poesen, J. 2002). The diff erence can be explained by the morphological position of the gullies. Analyzing the role of soil texture the problem has to be dealt with that the AGROTOPO texture data give only one single value for the total soil profi le and textural diff erences between the soil horizons are not taken into account. Loam is overrepresented, with a higher ratio of gullies on it than the territorial distribution of loam. In Zselic the percentage of gullies on Luvisol is also high (Table 8) compared with the areal representation of this soil type. The reason for this Table 7. Average properties of the surveyed gullies in the Zselic Property Rate according to the number of gullies (%) Rate according to gully length (%) Average gully length (m) Soil type Luvisol Cambisol Gleysol 79 18 3 81 17 1 265 236 195 Parent rock Glacial and alluvial deposits Loess Tercier and older deposits 2 96 2 2 97 1 224 260 144 Soil texture Loam Clay loam 98.3 1.0 99 1 259 144 OM content t ha-1 50–100 100–200 > 400 79 20 1 81 18 0 264 232 191 Soil depth (cm) 70–100 > 100 1 99 0 99 191 258 337 is the same as in the Börzsöny, i.e. the morphological position of the gullies. The share of gullies developed on Cambisol is proportional with the territo- rial extension. On Gleysol and Chernozem there are hardly any gullies. This is normal concerning Gleysol but there should have been more gullies on Chernozem. This soil type is used as arable land where ephemeral gullies de- velop being not always shown on topographic maps (Jakab, G. et al. 2010a). Gully distribution data on various parent rocks refl ect the well known fact that gully development favours loess environment (Poesen, J. et al. 2003, 2005; Valentin, C. et al. 2005). Comparing the areal percentage of parent rocks with the percentage occurrence of gullies on them we can see that loess is over- represented while glacial and alluvial deposits are underrepresented. Analyzing the role of organic matt er the well known positive eff ect of organic matt er in preventing gully erosion can be recognized, i.e. with low OM content more gullies develop. There is a close relationship between gully erosion and land use (Table 8). In both pilot areas most of the surveyed gullies were in the forest when the map was prepared. If gullies are deeply incised arable cultivation must be stopped and these gullies will not be classifi ed into arable land any more (Jakab, G. 2006). These gullies will soon be covered by forest. The overwhelming area of Börzsöny is covered by forest, but the pro- portion of gullies outside the forest area is higher than in the forest. In the Zselic most of the gullies are located in the forest (90%) in spite of the fact that the percentage of forests is less. Gullies on cultivated land generally are longer as shown by Vanwalleghem, T. et al. (2003). Diff erences in average gully length among the land use categories in the Börzsöny are negligible, in the Zselic area they are bigger. Land use types of the gullies in 2000 are presented in Table 9. Since 1985 the percentage of the number of gullies on arable land had increased in both areas, in Börzsöny from 2 to 6%, in Zselic from 3 to 8% and average gully length from 2 to 13% and 3 to 22%, respectively. Table 8. Land use type distribution of the gully areas in 1985 Indicator Börzsöny Zselic Arable land Pasture Forest Arable land Pasture Forest Total length (m) Total length (%) Number of gullies Ratio of gullies (%) Average length (m) Median length (m) Minimum length (m) Maximum length (m) 11,343 2 37 2 307 161 14 2,558 76,619 12 266 12 290 131 22 6,990 550,347 86 1,959 86 281 124 2 7,308 45,032 3 201 3 226 184 2 1,235 135,941 8 867 13 157 110 2 2,430 1,512,401 89 5,512 84 274 135 2 22,124 338 Intensive land use increases the rate of both sheet and gully erosion (Centeri, Cs. 2002; Gábris, Gy. et al. 2000). A similar trend is connected with chang- ing land use, i.e. with the transformation of the former forest and pasture into ar- able land. The dissection index values vary between 0.8 and 2.3. Arable land and pas- ture have similar values in both areas. Much higher values should have been on arable land. The reason why this is not the case is that ephemeral gullies are not surveyed as mentioned above (Jakab, G. et al. 2010b) The forest can also be dis- sected due to the dirt roads running in them and because of the eff ect of rill and gully erosion taking place in the arable fi eld upward the slope (Jakab, G. et al. 2010a). Analyzing the trends of changes the following statements can be made. In Zselic 66% of the gullies which were classifi ed as forest in 1985 had become arable land by 2000, in the Börzsöny this value is 50%. It is diffi cult to evaluate the ef- fect of gradient on gully formation as an elongated form of sometimes more than 10 km length is characterized by only one gradient value. Average gradient in the Börzsöny Mountains is 2.98, in the Zselic hills it is 2.69. Conclusions Two pilot areas were analyzed and com- pared in detail. As it was expected the gullies and gully systems of the pilot ar- eas were diff erent. It is diffi cult to assess the eff ects of the environmental factors Ta bl e 9 . L an d us e t yp e d ist rib ut io n of th e g ul ly a re as d er iv ed fr om th e 2 00 0 CL C5 0 da ta ba se In di ca to r Bö rz sö ny Zs el ic C ul tiv at ed ar ea Fo re st Pa st ur e Sc ru bl an d C ul tiv at ed ar ea Fo re st Pa st ur e Sc ru bl an d To ta l l en gt h (m ) To ta l l en gt h (% ) N um be r o f g ul lie s Ra tio o f g ul lie s (% ) A ve ra ge le ng th (m ) M ed ia n le ng th (m ) M in im um le ng th (m ) M ax im um le ng th (m ) D is se ct io n in de x (k m k m -2 ) 82 ,9 68 13 12 6 6 65 8 26 5 22 6, 99 0 2. 0 45 5, 93 3 73 1, 81 2 82 25 2 11 9 2 5, 56 1 1. 3 81 ,7 18 13 24 2 11 33 9 15 0 17 7, 30 8 2. 0 6, 24 7 1 33 1 18 9 11 5 20 1, 07 8 – 37 0, 86 2 22 1, 12 8 18 32 9 18 0 2 16 ,2 20 0. 8 1, 08 8, 53 1 64 4, 39 8 67 24 8 10 8 2 8, 23 7 2. 5 19 7, 87 2 12 84 7 13 23 4 14 3 6 3, 08 1 1. 0 30 ,4 91 2 16 5 3 18 5 10 3 18 2, 53 1 – 339 because they are interrelated, not independent from each other, i.e. they consti- tute a complex system and this system as a whole controls gully development. Being aware of this statement it is assumed that soil properties exert the small- est impact on gully development. This conclusion is confi rmed by the high share of gullies formed on Luvisol and it is known that Luvisol is resistant to gully erosion. The role of parent material in the process is very important. High gully erosion rate can be found on loess as the parent material. Relief and land use play the most important role in gully development. The eff ects of these two environmental factors cannot be treated separately, except in very small areas. In accordance with the results of Menéndez-Duarte, R., et al. (2007) the eff ect of topography is more decisive of the two because it controls also land use. The methodology of gully identifi cation is not perfect as such formation is classifi ed as a gully in the forest even if a considerable part of its catchment is on arable land. The method applied in the paper for the identifi cation of gully gradient is suitable in small scale only. In the Börzsöny Mountains very deep gullies can be found that have developed into valleys in some cases. They must be of Pleistocene, early Holocene origin as they are deeply cut into hard rock. The analysis of the two pilot area revealed important characteristics of the gullies and pointed to some features of gully development. Future research will be devoted to the classifi ed survey of gullies in the country (i.e. a detailed country-wide survey of ephemeral gullies). Acknowledgement: Research activities reported in this paper were funded by the Hungarian Scientifi c Research Fund (OTKA K 76434) and the support is gratefully acknowledged here. The data input was carried out by Varga, E. and this activity is also gratefully ac- knowledged by the authors. REFERENCES Borcsik, Z., Farsang, A., Barta, K. and Kitka, G. 2011. Humuszanyagok mennyiségi és minőségi eróziójának mérése a Tolna megyei Szálka település mellett i vízgyűjtőn. (Measurement of qualitative and quantitative caracteristics of soil erosion in a catch- ment, near Szálka, Tolna County.) Talajvédelem (Supplement) 127–137. Buzek, L. 2007. Water erosion in the watershed of the Upper Ostravice river from 1976 to 2000. Moravian Geographical Reports 15. (4): 2–14. Büttner, G., Feranec, J. and Jaffrain, G. 2002. Corine land cover update 2000. European Environment Agency, Copenhagen, Denmark Centeri, Cs. 2002. The role of vegetation cover in the control of soil erosion on the Tihany peninsula Acta Botanica Hungarica 44. (3–4): 285–295. Centeri, Cs., Herczeg, E., Vona, M., Balázs, K. and Penksza, K. 2009. The eff ects of land- use change on plant-soil-erosion relations, Nyereg Hill, Hungary. Journal of Plant Nutrition and Soil Science 172. (4): 586–592. 340 Conforti, M., Aucelli, P., Robustelli G. and Scarciglia, F. 2011. Geomorphology and GIS analysis for mapping gully erosion susceptibility in the Turbolo stream catchment (Northern Calabria, Italy) Natural Hazards 56. 881–898. Csathó, P., Sisák, I., Radimszky, L., Lushaj, S., Spiegel, H., Nikolova, M. T., Nikolov, N., Čermák, P., Klir, J., Astover, A., Karklins, A., Lazauskas, S., Kopiński, J., Hera, C., Dumitru, E., Manojlovic, M., Bogdanović, D., Torma, S., Leskošek, M. and Khristenko, A. 2007. Agriculture as a source of phosphorus causing eutrophication in Central and Eastern Europe. Soil Use and Management 23. 36–56. Dövényi, Z. (ed.) 2010. Magyarország kistájainak katasztere. (Inventory of Natural Micro-re- gions of Hungary.) Budapest, Geographical Research Institute HAS, 876 p. Farsang, A., Kitka, G. and Barta, K. 2011. Mezőgazdaságilag hasznosított kisvízgyűjtők talajerózióhoz kötődő elemdinamikája. (Element dinamics of cultivated watersheds related to soil erosion) Talajvédelem (Supplement) 339–349. Gábris, Gy., Kertész, Á. and Zámbó, L. 2003. Land use change and gully formation over the last 200 years in a hilly catchment. Catena 50. 151–164. Galang, M.A., Markewitz, D., Morris, L.A. and Bussell, P. 2007. Land use change and gully erosion in the Piedmont region of South Carolina. Journal of Soil and Water Conservation 62. 122–129. Grace, J.M. 2004. Soil erosion following forest operations in the Southern Piedmont of central Alabama. Journal of Soil and Water Conservation 59. 180–185. Gyenizse, P., Elekes, T., Nagyváradi, L. and Pirkhoffer, E. 2008. A tájhasználat-változás értékelése a Zselicben, az Árpád-kortól a 20. század végéig. (Assessment of ladscape change in the Zselic micro-region from the 12th to 20th century.) Földrajzi Közlemények 132. (1): 83–88. Hegedűs, K., Horváth, G., Karancsi, Z. and Prakfalvi, P. 2008. Eróziós vizsgálatok a Medves-vidék egy homokkőszurdokában. (Soil erosion investigations in a sandstone gully of the Medves Area.) Földrajzi Közlemények 132. (2): 157–173. Horváth, G., Pintér, Z. and Prakfalvi, P. 2010. Erosional micro-, meso- and macroforms of the Medves Area with special respect to the development of the sandstone gorges. Budapest, Report for Hungarian Scientifi c Research Fund (OTKA). htt p://real.mtak.hu/2254/ Jakab, G., Madarász, B. and Szalai, Z. 2009. Gully- or sheet erosion? A case study at catch- ment scale. Hungarian Geographical Bulletin 58 (3): 151–161. Jakab, G. 2006. A vonalas erózió megjelenésének formái és mérésének lehetőségei. (Gully types and posssibilities of their investigation.) Hungarian Journal of Landscape Ecology 4. (1): 17–33. Jakab, G. 2008. Természeti tényezők hatása a talajpusztulás vonalas formáinak kialakulására (The in- fl uence of environmental conditions on gully development). Budapest, Geographical Research Institute HAS. 109 p. (PhD dissertation) Jakab, G., Kertész, Á. and Papp, S. 2005. Az árkos erózió vizsgálata a Tetves-patak vízgyűjtő- jén. (Gully erosion in the Tetves catchment.) Földrajzi Értesítő/Hungarian Geographical Bulletin 54. (1–2): 149–165. Jakab, G., Kertész, Á. and Szalai, Z. 2010b. Scale dependence of gully investigations. Hungarian Geographical Bulletin 59 (3): 319–330. Jakab, G., Kertész, Á., Madarász, B., Ronczyk, L. and Szalai, Z. 2010a. Az erózió és a dom- borzat kapcsolata szántóföldön, a tolerálható talajveszteség tükrében. (The role of relief in soil erosion with special emphasis on tolerable soil loss.) Hungarian Journal of Landscape Ecology 8. (1): 35–45. Kapronczay, J 1966. Adatok a Zselic éghajlatához. (Contributions to the climate of Zselic.) Földrajzi Értesítő/ Hungarian Geographical Bulletin 15. (2): 179–197. 341 Kapronczay, J. 1965. Adatok a Zselic geomorfológiájához. (Data to the geomorphology of Zselic.) Földrajzi Értesítő/Hungarian Geographical Bulletin 14. (1): 29–45. Karátson, D. 2007. A Börzsönytől a Hargitáig (From Börzsöny to Harghita). Budapest, Typotex, 463 p. Kertész, Á. 2006. Az éghajlati tényezők szerepe a lineáris erózióban. (The role of climatic conditions in rill and gully erosion.) Földrajzi Közlemények 130. (3–4): 115–122. Kertész, Á. 2008. Water and wind erosion in Hungary. In Dimensions and trends in Hungarian geography. Eds. Kertész, Á. and Kovács, Z. Studies in Hungarian Geography 33. Budapest, Geographical Research Institute HAS. 47–54. Kertész, Á. 2009. Environmental conditions of gully erosion in Hungary. Hungarian Geographical Bulletin 58. (2): 79–89. Kirkby, M.J. and Bracken, L.J. 2009. Gully processes and gully dynamics. Earth Surface Processes and Landforms 34. 1841–1851. Kirkby, M.J. and Bull, L.J. 2000. Some factors controlling gully growth in fi ne-grained sedi- ments: a model applied to southeast Spain. Catena 40. 127–146. Láng, S. 1955. A Mátra és a Börzsöny természeti földrajza (Physical geography of the Mátra and Börzsöny mountains). Budapest, Akadémiai Kiadó, 512 p. Lóczy, D. and Gyenizse, P. 2003. A dunántúli dombságok változó hasznosítása és értéke a Zselic példáján. (Varying land uses and changing values of Transdanubian hilly countries on the example of the Zselic micro region.) In A Dunántúl és a Kisalföld történeti földrajza. Eds. Frisnyák, S. and Tóth, J. Nyíregyháza–Pécs, Nyíregyházi Főiskola Földrajz Tanszéke – PTE Földrajzi Intézete, 165–180. Madarász, B. 2009. A magyarországi erubáz talajok komplex talajtani vizsgálata, különös tekintet- tel agyagásvány-összetételükre. (Complex pedologic analysis of the Hungarian erubáz soils with special reference to their clay-mineral composition.) Budapest, Eötvös Loránd University, 134 p. (PhD dissertation). Madarász, B. and Jakab, G. 2009. Volcanic soils of the High Börzsöny and their relationship with geomorphological conditions. Hungarian Geographical Bulletin 58. (4): 227–242. Menéndez-Duarte, R., Marquínez, J., Fernández-Menéndez, S. and Santos, R. 2007. Incised channels and gully erosion in Northern Iberian Peninsula: Controls and geomorphic sett ing. Catena 71. 267–278. Muńoz-Robles, C., Reid, N., Frazier P., Tighe, M., Briggs, S. and Wilson, B. 2010. Factors related to gully erosion in woody encroachment in south-eastern Australia. Catena 83. 148–157. Nachtergaele, J. and Poesen, J. 2002. Spatial and temporal variations in resistance of loess-derived soils to ephemeral gully erosion. European Journal of Soil Science 53. (3): 449– 463. Pécsi, M. 1991. Geomorfológia és domborzatminősítés. (Geomorphology and relief classifi ca- tion.) Budapest, Geographical Research Institute HAS. 296 p. Pécsi, M. 1997. Szerkezeti és váztalaj képződés Magyarországon. (Eff ect of the Quaternary peri- glacial processes and structural soil formation in Hungary) Budapest, Geographical Research Institute HAS. 99 p. Pécskay, Z., Lexa, J., Szakács, A., Balogh, K., Seghedi, I., Konecny, V., Kovacs, M., Márton, E., Kaliciak, M., Széky-Fux, V., Póka, T., Gyarmati, P., Edelstein, O., Rosu, É. and Zec, B. 1995. Space and time distribution of Neogene-Quaternary volcanism in the Carpatho–Pannonian Region. Acta Vulcanologica 7. (2): 15–28. Pintér, Z., Prakfalvi, P., Karancsi, Z. and Horváth, G. 2009. Medves-vidéki riolitt ufák eróziós formakincse. (Erosional forms on rhyolite tuff s in Medves Region.) Földrajzi Közlemények 133 (3): 219–238. 342 Poesen, J., Nachtergaele, J., Verstraeten, G. and Valentin, C. 2003. Gully erosion and environmental change: importance and research needs. Catena. 50. 91–133. Poesen, J., Vanwalleghem, T., de Vente, J., Knapen, A., Verstraeten, G. and Martínez- Casanovas, J. 2006. Gully erosion in Europe. In Soil erosion in Europe. Eds. Bordman, J. and Poesen, J. New Jersey, John Wiley and Sons, 515–536. Rabus, B., Eineder, M., Roth, A. and Bamler, R. 2003. The shutt le radar topography mission- a new class of digital elevation models acquired by spaceborne radar, Photogrammetric Engineering and Remote Sensing 57. 241–262. RISSAC 1991. htt p://www.mta-taki.hu/hu/osztalyok/gis-labor/agrotopo Sebe, K., Csillag, G. and Konrád, Gy. 2008. The role of neotectonics in fl uvial landscape development in the Western Mecsek Mountains and related foreland basins (SE Transdanubia, Hungary) Geomorphology 102 (1): 55–67. Smolska, E. 2007. Extreme Rainfalls and their Impact on Slopes — Evaluation Based on Soil Erosion Measurements (As Exemplifi ed by the Suwałki Lakeland, Poland) Geographia Polonica 80. (2): 151–164. Valentin, C., Poesen, J. and Li, Y. 2005. Gully erosion: impacts, factors and control. Catena 63. 132–153. Vanwalleghem, T., Eeckhaut, Van Den M., Poesen, J., Deckers, J., Nachtergaele, J., Oost, Van K. and Slenters, C. 2003. Characteristics and controlling factors of old gullies under forest in a temperate humid climate: a case study from the Meerdaal Forest (Central Belgium). Geomorphology 56. (1–2): 15–29. Zglobicki, W. and Baran–Zglobicka, B. 2011. Gullies as an indicator of human impact on loess landscape (Case study: North Western Part of Lublin Upland, Poland) Zeitschrift für Geomorphologie 55. 119–137. << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles true /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Error /CompatibilityLevel 1.3 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.0000 /ColorConversionStrategy /LeaveColorUnchanged /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams false /MaxSubsetPct 100 /Optimize false /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments true /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile () /AlwaysEmbed [ true ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 300 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.50000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 300 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.50000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.50000 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile (None) /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName () /PDFXTrapped /False /CreateJDFFile false /Description << /ARA /BGR /CHS /CHT /CZE /DAN /DEU /ESP /ETI /FRA /GRE /HEB /HRV (Za stvaranje Adobe PDF dokumenata najpogodnijih za visokokvalitetni ispis prije tiskanja koristite ove postavke. Stvoreni PDF dokumenti mogu se otvoriti Acrobat i Adobe Reader 5.0 i kasnijim verzijama.) /ITA /JPN /KOR /LTH /LVI /NLD (Gebruik deze instellingen om Adobe PDF-documenten te maken die zijn geoptimaliseerd voor prepress-afdrukken van hoge kwaliteit. De gemaakte PDF-documenten kunnen worden geopend met Acrobat en Adobe Reader 5.0 en hoger.) /NOR /POL /PTB /RUM /RUS /SKY /SLV /SUO /SVE /TUR /UKR /ENU (Use these settings to create Adobe PDF documents best suited for high-quality prepress printing. Created PDF documents can be opened with Acrobat and Adobe Reader 5.0 and later.) /HUN >> /Namespace [ (Adobe) (Common) (1.0) ] /OtherNamespaces [ << /AsReaderSpreads false /CropImagesToFrames true /ErrorControl /WarnAndContinue /FlattenerIgnoreSpreadOverrides false /IncludeGuidesGrids false /IncludeNonPrinting false /IncludeSlug false /Namespace [ (Adobe) (InDesign) (4.0) ] /OmitPlacedBitmaps false /OmitPlacedEPS false /OmitPlacedPDF false /SimulateOverprint /Legacy >> << /AddBleedMarks false /AddColorBars false /AddCropMarks false /AddPageInfo false /AddRegMarks false /ConvertColors /ConvertToCMYK /DestinationProfileName () /DestinationProfileSelector /DocumentCMYK /Downsample16BitImages true /FlattenerPreset << /PresetSelector /MediumResolution >> /FormElements false /GenerateStructure false /IncludeBookmarks false /IncludeHyperlinks false /IncludeInteractive false /IncludeLayers false /IncludeProfiles false /MultimediaHandling /UseObjectSettings /Namespace [ (Adobe) (CreativeSuite) (2.0) ] /PDFXOutputIntentProfileSelector /DocumentCMYK /PreserveEditing true /UntaggedCMYKHandling /LeaveUntagged /UntaggedRGBHandling /UseDocumentProfile /UseDocumentBleed false >> ] >> setdistillerparams << /HWResolution [2400 2400] /PageSize [612.000 792.000] >> setpagedevice