UDK 903:911.52:903.4(497.12)"634/£36" Documenta Praehistorica XXVI Landscape dynamics on the Ljubljana Moor Dimitrij Mlekuz Department of Archeology. University of Ljubljana, Slovenia dimitrij.mlekuz® kiss.uni-lj.si ABSTRACT - A geographical information System (GIS) based dynamical model of streams migration iras constructed in order to gain insight into the dynamics of the palaeo landscape Model implications related to the settlement patterns and perception of the landscape by the mid-Holocene men and women dwelling in il are discussed IZMEČEK - Članek predstavlja uporabo geografskega informacijskega sistema (GIS) pri rekonstrukciji dinamike poplavne ravnice na Ljubljanskem barju Predstavljen je model dinamike krajine Ljubljanskega barja v srednjem holocenu in model nastanka poselitvenega vzorca. KEY ViORDS - Ljubljansko barje; MesoHthic: SeoHthic, I: neolithic, Geographical Information System, allutiul geoarchaeology; landscape dynamics The purpose of modeling is insight, not numbers (Hamming 1%2) INTRODUCTION The aim of the present paper is to access the palaeo-dynamics of the Ljubljana Moor landscape in the period between 6500 BP and -«(KM) BP. My intention is to take one step beyond the traditional debate concerned mainly with the description of the static site environs and the detection of the human impact on the landscape and which is - at least in my opinion -unable to access long-term (»evolutionary dy namics of human groups and landscape on the Ljubljana Moor. The following text aims to suggest that Geographical Information Systems (GIS) can provide not only a sophisticated cartographic tool, but a flexible environment within which the dynamics of the past landscapes can be modelled and explored. Using an on-going case-study focused upon the landscape dynamics of the Ljubljana Moor, a dynamic GIS model was constructed in order to explore the dy namics of the floodplain and its implications in the context of the Mesolithic/Neolithic/Eneolithic occupation of the Ljubljana Moor. BACKGROUND TO THE STUDY AREA Ljubljana Moor or Ljubljansko barje (Fig. 1) Ls a large wetland in the extreme south of the Ljubljana basin in the central part of Slovenia. It is a tectonic depression. with an extensive alluvial floor and flat surface, fringed by three major topographical regions, the Alpine foothills, carstic Dolenjska region and northern tip of carstified Dinaric mountains. The bottom can be generally divided into three basic topographic units, the most extensive being a marshy flood-plain, then alluvial fans and. finally, isolated hills. The whole area Ls characterised by dynamic tectonic activity taking forms from long-term subsidence and uplift to catastrophic earthquakes, the last one in 1895. The main water-courses are Ljubljanica, Izica and I.v ka rivers, with extensive carstic watersheds in the Dinaric Mountains, and the River Grada&ica from the Alpine foothills. Prior to the major regulation works and commercial peat extraction in the 19th century , parts of the area were covered by up to 6 m of peat and characterised by predictable cyclical annual floods nested within an unpredictable cy cle 185 Dimitrij Mtekui Fig. 1. False colour sattelite image of Ljubljana Moor with the places mentioned in the text. Landsat TM. channels 2.5.7. Source: ElK IMAGE. Frascatl. Italy. Owner of the data: Joint Research Centre. Ispra. Italy. Owner of the geo-referenced mosaic: Zarod Republike Slovenije za Statistiko. Slovenia. of major flood events. Today, as a result of flood control works, the landscape is effectively controllable and stable. RECENT ARCHAEOLOGICAL RESEARCH Archaeological interest in the region started in the second half of the 19'1' century when the first pile dw ellings were discovered. This discovery was a part of the wider cultural historical context in the second half of the 19'h century , when a chain reaction of pile-dwelling discoveries was triggered by the initial discovery of the pile dwelling on the Lake Zurich (Greif 1997.72). The peaceful landscape of a pile-dwelling situated at the edge of a lake surrounded by mountains served as a model for interpreting Neolithic settlement patterns and palaeo environment [Keller 1854). In the case of the Ljubljana Moor this romantic picture was supported by the generalised stratigraphic sequence of chalky lacustrine clays and marls followed by organogenic muds and finally peat interpreted as a result of linear hydroseral succession from lake to swamp and raised bog (Melik 1946: Sercelj 1966.443) A vicious circle w as formed when archaeological data was used to support environmental data and vice versa. The model created in the mid-19th century was applied to the archaeological and palaeo environmental data at the Ljubljana Moor until the 1990s. This was a period when the search for pile-dwellings was fashionable and intensive - although poorly documented - excavations were important for the insti- tutional promotion of Slovenian prehistoric archaeology (Budja 1995.176). The floor of Ljubljana Moor was priv ileged and other parts of the Ljubljana basin and background hills were neglected until the 1980s, when the discovery of Mesolithic lithic scatters (JosipoviC ¡985) focused research interest on -and only on - the fringes of the Moor and at edges of isolated hills. But new data was still collected un-systematically and researchers ignored the Moor's hinterland. The Mesolithic settlement pattern w as interpreted is an indicator of lakeside activities and w as used to support the model of gradual and linear succession from lake to swamp and raised bog. In 1995 Budja proposed an alternative scenario for the landscape dynamics and settlement patterns based on Sifrer {1984) flood-plain hypothesis. Sifter interprets chalky silt and clay deposits as die sign of erosion of v arious soils in the Dinaric Karst region at the beginning of the Holocene. These sediments are interpreted as indicators of climate changes and past flooding of the Ljubljana Moor, rather than as proofs of a lake existence. Budja (/995) explained Neolithic and Eneolithic sites as riparian settlements and their pattern as a human response to flood-plain evolution. ACCESSING PALAEO DYNAMICS Detailed study area The eastern part of the Ljubljana Moor was chosen for intensive research into floodplain dynamics. A 186 Landscape dynamics on the L)ut>l|ana Moor detailed digital elevation model (DEM) was generated and floodplain features such as lakes, oxbow lakes and palaeo channels were identified from aero-photos and DEM. The pre-regulation hydrological network was reconstructed from historical maps (Fig. 2). Settlement patterns The remnant settlement pattern shows a clear preference for the floodplain by the Neolithic, Eneolithic "pile-dwellings" and a position just next to the palaeo channels. The distribution of the Mesolithic/ Early Neolithic lithic scatters is limited to the edges of the floodplain. Spatial statistics tests were performed in order to observe the relationship between sites and environmental variables, in this case, distance to the nearest river and palaeo channel. Firstly, the relationship between sites and the pre-regulation hydrological network was tested using Kol-mogorov -Smirnov two sample test {Hodder and Orion 1976.226-229). Samples of the 32 pile dwellings and 19 lithic scatters was used in the analysis. The difference between distributions of pile-dwellings and streams proved to be significant at well below the 1% level. Pile-dwellings show a marked preference for proximity to water, with about 90% of the sites located less than 500 m from the nearest stream. Lithics scatters display a similar, but less marked preference, and the difference between distributions was proven not to be significant at the 5% level (Fig. 3). Another test was performed to observe the relationship between the location of pile-dwellings and palaeo channels in the detailed study area (Fig. 4). 13 pile-dwellings w ere used in analysis and the difference between two distributions was proven to be significant well below the 1% level. Possible interpretations of the observ ed pattern are: • There were no large scale river migrations since the Neolithic/Eneolithic. 0 500 1000 1500 2000 2500 300C Distance horn nearest river [m] Fig. 3- Relationship betu een number of sites and background (expected number of sites) at distances from nearest rirer (obtained from historical map) on a regional (Ljubljana Moor) scale. 187 Dim--lri| MlekuZ Fig. 4. Relationship between number of pi led uel-lings and background (expected number of sites) at distances from nearest palaeo channel on a cat-chement (Detailed study area) scale. • Pile-dwellings show a preference for riparian environments. • Proximity to running water was not the key factor for choosing the location for activities, observed today as lithic scatters or • drastic landscape restructuring (drying up of the lake and formation of the floodplain) took place in the Late Mesolithic. How ever, this is a modern, quantifiable and first of all static landscape, and the next step is to access the dynamics which shaped it in the past. Processes that shaped the floodplain Looking at the long-term structuring of the Ljubljana Moor landscape it is obvious that floods are the main factor which contributed to the shaping of the flood-plain. However, it was high magnitude, unpredictable floods that changed the floodplain dramatically. There are four interlinked processes w hich are the consequences of flood events of different magnitudes and which can be read from the modern landscape. Levee formation Levees (Fig. 5) are linear features located just next to a channel. They are formed by the deposition of sediment when flood-water exceeds bank height and its velocity drops abruptly (van Atidel and Runnels /995.485, Brown /997, Table I.I; Knighton 1998.145). Meander migration In meandering systems (as most of Ljubljana Moor rivers are) the most common type of channel change is meander migration (Fig. 6). This can be classified into a number of different forms, such as rotation, translation, extension, enlargement and combinations of these (Brown 1997.26-27; Knigltlon 1998. 226-227). Meander neck cutoff High energy flood water can find a shorter route and cut-offs meander into neck (Fig. 7). Result is a distinctive oxbow lake or swamp (Brown 1997.28). Meander neck cut-off can also be caused in a low magnitude, cyclical flood event. Avulsion The most dramatic process is called avulsion (Fig. 8). The gradual rise of the river bed and adjacent levees causes the river to seek a new course lower down in the flood-plain. Avulsion can be triggered by a high magnitude flood. Ribbon lakes or palaeochannels are typical result of this process (Brown 1997.28; van Andet and Runnels 1995. Fig. 2. 484). The model A simple GIS based dynamic model of stream behaviour was constructed in order to access the long-term structuring of the floodplain. It simulates the Levee Backsuamp Flood level ^ \ Channel \ Fig. 5. Levee formation. 188 Landscape dynamo on the Ljubljana Moot Oxbow lake or swamp Fig. 7. Meander neck cutoff. effect of various permanent alterations inflicted upon a landscape by the action of floods events of different magnitudes. G1S was then used to apply a group of hydrological modelling algorithms to the post-flood topographical data in order to simulate the movement and behaviour of water flows across the post-flood landscape. I"he modelled landscape has it own history as topographic changes of previous flood events are inherited and contribute to the shaping of subsequent events. Stream course data obtained after each run of the model were overlaid to observe the long-term effects of floodplain dynamics. The resultant map (Fig. 9) can be seen as a fuzzy set, with different levels of probability (uncertainty) for fluvial activity at each grain of the floodplain. A by-product of the model is an animation of stream migration which offers a truly dynamic view of the floodplain and replaces the static abstraction of the map. Fig. 9. Model results after 200 runs. Dark shaded areas indicate areas of high probability of fluvial activity. »««ncn ♦ Pitr-doHliniP Mew channel Ribbon . * take *«•• (Old channel) Backswamp Old channel \ New channel Fig. H. Avulsion. RESULTS What strikes me most how dynamici the simulated floodplain is. Simulated streams change their course after each flood event. Small-scale stream migrations and large scale migrations aLso took place, but not 189 Dimitni Mlekuz Fig. 10. Very detailed digital elevation model (DEM) of a small area of the floodplain. Several palaeochan-nel.% and adjoining teeves can he observed. Note the position of the pile-duellings next to the pataeochan nets. randomly. Several possible corridors for each stream occur and flow oscillates between them. If we translate this into the palaeo landscape, it means that the major consequences of the floods were not simply effects of the inundation, of high w ater, but the migration of w ater streams w hich reshaped the flood-plain and changed Its structure. The floodplain was probably a dynamic landscape even on a human time scale. Another feature is die enormous sensitivity of the simulated floodplain to disturbance. Simulating even the smallest changes in die topography caused major changes in the shape and structure of the floodplain. A corollary of this observation could he that human intervention in the drainage catchement of the flood plain or in the floodplain itself - forest burning, deforestation. grazing, cultivation - could have had a profound effect on the landscape and the future actions of the people inhabiting it. An important feature of the simulated floodplain is the occurrence of isolated patches almost completely unaffected by fluvial activity. These are areas slightly elevated from the floodplain. Lithic scatters are located on some of these slightly elev ated areas. Absence of ev idence of lithic scatters on the floodplain is, in my opinion, not evidence of 190 Landscape dynamics on ihe Ljubljana Moor Backxwanip lev« IIg. II. Model of site relocation. Site located on a levee next to Ihe active channel before flood. After high magnitude flood event, uhich restructures sites environ and cuts off site from the running water, site is relocated to Ihe active channel. absence. It is a result of a selective preservation of the sites with permanent structures (pile dwellings), whereas fragile Mesollthic/Earfy Neolithic sites of activities were washed away or rendered invisible by fluvial processes. Neolithic and Eneolithic "pile-dwellings" display a different pattern. They are located just next to or in the stream corridors of die simulated flood-plain and just next to the identified palaeo channels in the real floodplain. As seen above, this pattern can be observed all over Ljubljana Moor. .As suggested also by the stratigraphy of the excavated pileHlwellings (Bregard 1975: Budja 1995). sites were obviously located next to active channels on natural levees (Fig. 10), which were dry most of the time and never deeply submerged by flood-waters. Flood can have several effects on the site: • Inundation, which could have caused structural damage to the site. A human answer to that would be rebuilding of the site, probably at the same location. • Bank erosion, which could have promoted undercutting of the levee and eventual bank collapse and damage to the site. A structure interpreted as a wave-break which w as discovered at the Mahar-ski prekop site could have served to protect the site from bank erosion (Bregant 197. Fig. I). • Stream migration and av ulsion, as a result of catastrophic flooding which would have effectively changed the shape of the floodplain, restructured site environs and cut the site off from the active channel and running water. Human action would be the relocation of the site to the next active channel (Fig. 11). I believe that proximity to running w ater was a key factor in choosing a place to live. • High density and number of sites can be seen in this perspective as a result of short duration of occupation and frequent relocation of the sites due to the effects of the catastrophic floods. Short duration of occupation was proved also by den-drochronological research on two sites to be about 50-80 years per site (Cufar et al. 1998). On the other side, long spans of radiocarbon data from sites Maharski prekop (Budja 1995174) and Parti (Harej 1978.74; Harej 1982.46) suggest either long duration of occupation or more possibly, several re-occupation of the sites. DISCUSSION What can this exercise tell us about the perception of the flood-plain landscape by the people dw elling in it? I believe that cyclical, annual, predictable floods formed an integral part of the everyday lives of the people inhabiting flood-plain, just like the presence 191 Dimitrij Mleku* of enriched soils, reed beds, wild-fowl, fish, water chestnut and other malign and benign things that die flood-plain offers. The floods would have served to coastruct and maintain temporalities (GiUings 1998). However, it was rare and unpredictable, catastrophic floods which contributed to the long term structuring of the landscape. Their not being embedded in people's taxonomy of natural events and their being thought of as being external to the natural order of events, would have caused a major change in the way of life of mid-Holocene women and men. ACKNOWLEDGEMENTS This paper is a part of a diploma dissertation. Landscape dynamics and the first cultural landscapes at the Ljubljana Moor', at the Department of Archaeology, University of Ljubljana, under supervision by Professor Mihael Budja and is currently a work in progress. I would like to thank Professor Mihael Budja for his support during preparation of the dissertation and this paper. All mistakes, however, remain my own. REFERENCES van ANDEL T. H. and RUNNELS C. N. 1995. The earliest farmers in Europe. Antiquity 69: 481-500. BREGANT T. 1975. Kolišče ob Maharskem prekopu pri Igu - raziskovanja 1973- in 1974. leta (Der Pfah bau am Maharski-Kanal - Forschungen aus dem Jahren 1973 und 1974). Poročilo o raziskovanju neo-iita in eneolita v Sloveniji 4: 7-107. BROWN A. G. 1997. Alluvia! geoarchaeology. Cambridge Manuals in Archaeology. Cambridge University Press. BUDJA M. 1995. Spreminjanje naravne in kulturne krajine v neolitiku in eneolitiku na Ljubljanskem barju I (Landscape changes in the Neolithic and Eneolithic in Slovenia. Case study: Ljubljansko barje I). Poročilo o raziskovanju paleolitika, neolitika in eneolitika v Sloveniji 23: 163-181. ClJFAR K.. LEVANTlC T. and VELUŠČEK A. 1998. Dendrokronološke raziskave na koliščih Spodnje mostišče 1 in 2 ter Hočevarica (Dendrochronological investigations in the pile dwellings Spodnje Mostišče 1 and 2 and Hočevarica from Ljubljana moor). Arheološki vestnik 49. 75-92. G1LLINGS M. 1998. Embracing uncertainty and challenging dualism in the GIS-based study of a palaeo-flood plain. European Journal of Archaeology 1: 117-144. GREIF T. 1997. Prazgodovinska kolišča na Ljubljanskem barju (The Prehistoric pile dwellings in the Ljubljansko barje). Arheo 18. HAREJ Z. 1978. Kolišče v Partih pri Igu na Ljubljanskem barju (Der Pfahlbau in Parte bei Ig auf dem Moor von Ljubljana). Poročilo o raziskovanju pa-leolita. neolita in eneolita v Sloveniji 6. 61-75 1982. Kolišče v Partih pri Igu na Ljubljanskem barju - Raziskovanja 1978. in 1979. leta (Der Pfahlbau in Parte bei Ig auf dem Moor von Ljubljana -Forschungen in den Jahren 1978 und 1979) Poročilo o raziskovanju /mleoliUi. neolita in eneolita v Sloveniji 9-10:31-101. HODDERI. and ORTON C. 1976. Spatial analyisis in Archaeology. Cambridge University Press. JOSIPOVIČ D. 1985. Ob odkritju najdišč sileksov na Ljubljanskem barju. Kronika 33• /-5. KELLER F. 1854. Die keltischen Pfahlbauten in den Schweizerseen. Mittelungen der Antiquarischen Gesellschaß in Zürich (Bericht 1)1. 65-100. KNIGHTON D. 1998. Fluvial Forms and Processes. Arnold. MELIK A. 1946. Ljubljansko koliščarsko jezero in dediščina po njem. Dela SA7AJ. SF.RCELJ A. 1966. Pelodne analize pleistocenskih in holocenskih sedimentov Ljubljanskega barja. Razprave S AZU IX/9: 431-427. ŠIFRER M. 1984. Nova dognanja o geomorfološkem razvoju Ljubljanskega barja (New findings on the Geomorphological Development on the Ljubljansko barje). Acta Geographica 23/2. 9-45. 192