Application of Airborne LiDAR Data to the Archaeology of Agrarian Land Use: The Case Study of the Early Medieval Microregion of Bled (Slovenia)

The use of topographic airborne LiDAR data has become an essential part of archaeological prospection, particularly as a tool for detecting archaeological features in the landscape. However, its use for landscape reconstruction and understanding archaeological sites in their environmental context is still underutilised. To this end, we took an innovative approach to using LiDAR data as a means of discovering, documenting, and interpreting agricultural land use systems by looking for significant environmental variation within a microregion. We combined information from LiDAR-derived DEM derivatives with archaeological, geological, and soil data. We introduced two methodological innovations. The first is the modified wetness index, which combines the LiDAR-derived precision with the accuracy of the effective field capacity of the soil to obtain a very realistic predictor of soil quality. The second is the modified landform classification, a combination of topographic position index and visual geomorphological analysis, which amalgamates two of the most important predictive variables for the distribution of plant species. Our approach is demonstrated by a case study focusing on early medieval settlements in the context of agricultural land use in the subalpine microregion of Bled (Slovenia). It revealed that early medieval settlers were drawn to light soils with high water retention capacity. Such soils were particularly suitable for the cultivation of barley, which is known to have been one of the most important staple crops of the period, especially in colder climate such as subalpine. Soils with lower water retention capacity were not colonized until the eleventh century, which may signify the transition at that time to a higher level of agricultural organisation and wheat as a staple cereal food.


Introduction
Airborne Laser Scanning (ALS) or airborne Light Detection and Ranging data (hereafter LiDAR) are used in archaeology for visualization and detailed morphological analysis of the archaeological landscape. First and foremost, LiDAR has become an essential part of archaeological prospection as a tool for detecting archaeological features [1][2][3][4][5][6][7][8][9][10]. The free availability of LiDAR data in Slovenia since 2015 [11], for example, has led to the discovery of numerous archaeological sites and features-such as prehistoric settlements, prehistoric and Roman field systems, Roman military camps, and Late Antique settlements [12][13][14][15][16][17][18]particularly in densely forested areas. Moreover, LiDAR data allows each site or feature to be observed at different scales [19][20][21]. From the large "human" scale that provides overwhelming layout detail at the intra-site level to the small landscape scale where patterns of site distribution can be readily observed, it has broadened our understanding of archaeological and historical landscapes. LiDAR data, however, are only suitable for the detection of those archaeological features that are recognizable in the terrain morphology (either embedded, partially embedded, or standing features or standing objects [22]). Thus, the impact of LiDAR data on archaeology as a discipline was uneven. One area of low impact was the detection of early medieval settlements in Eastern Alpine region (hereafter

Archaeological Context
To introduce the case study, a brief overview of the relevant archaeological context is required. The settlement of the early medieval Slavs is a controversial topic in European archaeology [48][49][50][51], for there are no first-hand, written sources before the end of the ninth century and archaeological evidence is sparse compared to many other early medieval peoples. As a result, there is a propensity for sweeping explanations on, for example, how could large parts of eastern Europe have been Slavicized in a relatively short period of time between the sixth and seventh centuries AD [52]. A small but by no means insignificant piece of this puzzle is the settlement of the Slavs in the Eastern Alpine region. Given the technical nature of this paper, we use the term "Slavs" hereafter as a terminus technicus describing a particular archaeological record and material culture typical of the Eastern Alpine region from the 6th century AD to the eleventh century AD, for example, typical pottery, pit huts, burial rites; we understand these as regional features and consider ethnic identities and habitus as complex research questions to be addressed in specialized studies [53,54].
Previous attempts to understand the landscape context of EMS in the Eastern Alpine region often reduced observations to height above sea level and soil type. One early analysis found that Slavs in Slovenia settled mainly in upland areas with dry soils and tended to avoid plains, narrow valleys, and wet soils [36]. In a preceding analysis of the Bled microregion, the reconstruction of the field system located the most suitable areas for early medieval agriculture and concluded that local topography had a direct influence on the EMS location choice model [37,55,56]. A similar attempt to define the landscape type and soil type in which EMS occurs was made in Lower Austria. Under the term mesoregion, 36 EMS were analysed within their respective 5 km radii. Soil type and geomorphological

Archaeological Context
To introduce the case study, a brief overview of the relevant archaeological context is required. The settlement of the early medieval Slavs is a controversial topic in European archaeology [48][49][50][51], for there are no first-hand, written sources before the end of the ninth century and archaeological evidence is sparse compared to many other early medieval peoples. As a result, there is a propensity for sweeping explanations on, for example, how could large parts of eastern Europe have been Slavicized in a relatively short period of time between the sixth and seventh centuries AD [52]. A small but by no means insignificant piece of this puzzle is the settlement of the Slavs in the Eastern Alpine region. Given the technical nature of this paper, we use the term "Slavs" hereafter as a terminus technicus describing a particular archaeological record and material culture typical of the Eastern Alpine region from the 6th century AD to the eleventh century AD, for example, typical pottery, pit huts, burial rites; we understand these as regional features and consider ethnic identities and habitus as complex research questions to be addressed in specialized studies [53,54].
Previous attempts to understand the landscape context of EMS in the Eastern Alpine region often reduced observations to height above sea level and soil type. One early analysis found that Slavs in Slovenia settled mainly in upland areas with dry soils and tended to avoid plains, narrow valleys, and wet soils [36]. In a preceding analysis of the Bled microregion, the reconstruction of the field system located the most suitable areas for early medieval agriculture and concluded that local topography had a direct influence on the EMS location choice model [37,55,56]. A similar attempt to define the landscape type and soil type in which EMS occurs was made in Lower Austria. Under the term mesoregion, 36 EMS were analysed within their respective 5 km radii. Soil type and geomorphological context, which provided a description of the predominant landform types, were considered. The results showed that the EMS occur in two landscape types: (flood) plains and mountainous regions. Approximately half of EMS were located on alluvial river terraces, at least some of them within coeval floodplains on naturally elevated land. The other half of EMS was located in upland and hilly areas above 300 m a.s.l. In these areas, loess and brown earth soils were clearly preferred [38].
In the archaeologically relevant neighbourhood, river terraces and hills were also recognized as the predominant locations for EMS in Bohemia [39]. Similar conclusions regarding landscape preference, habitat description, and soil conditions were also drawn for Great Moravia in Czech Republic [57] and Slovakia [58] and for several microregions in Slovenia (Krško polje [59], Prekmurje and Podravje [60,61], and Bled [62,63]). A somewhat different situation was detected for the sixth-century Slavs in the Northern Danube region (present-day Slovakia, Moravia, Czech Republic, and Upper Austria), who settled the lowlands in strategic locations along roads and at river fords, while mountainous terrain was avoided [64] (Figure 2). context, which provided a description of the predominant landform types, were considered. The results showed that the EMS occur in two landscape types: (flood) plains and mountainous regions. Approximately half of EMS were located on alluvial river terraces, at least some of them within coeval floodplains on naturally elevated land. The other half of EMS was located in upland and hilly areas above 300 m a.s.l. In these areas, loess and brown earth soils were clearly preferred [38].
In the archaeologically relevant neighbourhood, river terraces and hills were also recognized as the predominant locations for EMS in Bohemia [39]. Similar conclusions regarding landscape preference, habitat description, and soil conditions were also drawn for Great Moravia in Czech Republic [57] and Slovakia [58] and for several microregions in Slovenia (Krško polje [59], Prekmurje and Podravje [60,61], and Bled [62,63]). A somewhat different situation was detected for the sixth-century Slavs in the Northern Danube region (present-day Slovakia, Moravia, Czech Republic, and Upper Austria), who settled the lowlands in strategic locations along roads and at river fords, while mountainous terrain was avoided [64] (Figure 2). As a note, it should be mentioned that EMS within floodplains would have severely restricted access to agricultural land. This suggests that the exploitation of riparian vegetation and other resources must have played an important and hitherto neglected role in early medieval economic life. The riparian zone was able to provide for fish, freshwater crabs, various edible plants; wild vines and similar could be gathered without having to invest in cultivation. Reeds for covering houses, but possibly also for making vessels, and willow twigs for building wattle walls in house construction could be gathered in the floodplain forests, as well as wood for timber construction [38,59].
Our particular interest in this paper is EMS in the Bled microregion in the context of agricultural land use. The Bled microregion (80 km 2 ) is located in the northwest of Slovenia, in the subalpine area of Julian Alps. The microregion is bounded by the confluence of As a note, it should be mentioned that EMS within floodplains would have severely restricted access to agricultural land. This suggests that the exploitation of riparian vegetation and other resources must have played an important and hitherto neglected role in early medieval economic life. The riparian zone was able to provide for fish, freshwater crabs, various edible plants; wild vines and similar could be gathered without having to invest in cultivation. Reeds for covering houses, but possibly also for making vessels, and willow twigs for building wattle walls in house construction could be gathered in the floodplain forests, as well as wood for timber construction [38,59].
Our particular interest in this paper is EMS in the Bled microregion in the context of agricultural land use. The Bled microregion (80 km 2 ) is located in the northwest of Slovenia, in the subalpine area of Julian Alps. The microregion is bounded by the confluence of the rivers Sava Bohinjka and Sava Dolinka in the east, and by the high mountain plateaus of Pokljuka and Mežakla in the west and north ( Figure 3). The area is notable for its intensive fluvio-glacial geomorphology. The archaeological significance of this microregion lies in the fact that it encompasses the entire territory of župa, which was the smallest administrative entity of the early medieval Slavs [37,65]. Bled has long been the focus of both archaeological and historical research, and from the point of view of early medieval archaeology, it is the best researched microregion in Slovenia. Since the 1880s, and most intensively in the 1970s and 1980s, 17 noteworthy early medieval archaeological sites have been documented by archaeological excavations [27,37,63,[66][67][68][69][70][71][72] (Table 1; Figure 1). the rivers Sava Bohinjka and Sava Dolinka in the east, and by the high mountain plateaus of Pokljuka and Mežakla in the west and north ( Figure 3). The area is notable for its intensive fluvio-glacial geomorphology. The archaeological significance of this microregion lies in the fact that it encompasses the entire territory of župa, which was the smallest administrative entity of the early medieval Slavs [37,65]. Bled has long been the focus of both archaeological and historical research, and from the point of view of early medieval archaeology, it is the best researched microregion in Slovenia. Since the 1880s, and most intensively in the 1970s and 1980s, 17 noteworthy early medieval archaeological sites have been documented by archaeological excavations [27,37,63,[66][67][68][69][70][71][72] (Table 1; Figure 1).  Table 1) and settlements (letters refer to Table 2) in the Bled microregion. The colours refer to the century of foundation (labelled as year AD in the legend).  Table 1) and settlements (letters refer to Table 2) in the Bled microregion. The colours refer to the century of foundation (labelled as year AD in the legend).   [37,73]); indirect-inferred indirectly, based on the landscape analysis (source [37]); written sources-terminus ante quem from written sources (after [37]).
Only one settlement in Bled area has been fully excavated (Pristava at Bled) and further two (Grad, Zasip) have been confirmed by excavations, but the chronology of several others could be inferred from their respective cemeteries. Remaining settlements were dated by terminus ante quem, provided in written sources or inferred indirectly from the landscape analysis and retrograde analysis of the historical cadastre (Table 2; Figure 3). However, no detailed and systematic archaeobotanical research has been carried out in the Bled microregion to date, and there are no published palynological results dealing with the early medieval vegetation in this area yet. Similarly, extensive underwater archaeological investigations of the Lake Bled yielded minor early medieval finds [74] but as yet no significant findings of relevance to this study.
Three decades have passed since the last comprehensive analysis of the Bled microregion, in which A. Pleterski combined archaeology, written sources, and retrograde analysis Remote Sens. 2021, 13, 3228 7 of 30 of historical cadastres [37,55]. He reconstructed the arable areas, which occurred in small patches scattered in the valley plains (Appendix A: Figure A1). His key conclusions were that most settlements were continuously inhabited from the early medieval period to the present time; the economic model was dominated by agriculture, with little developed crafts [72]. Therefore, each settlement was located adjacent to soils suitable for agriculture. Moreover, most, if not each, settlements had a cemetery nearby. The validity of the original study was subsequently confirmed with archaeological excavations on three separate locations in Žale near Zasip [75], Zasip, and Došca [76]. He was therefore able to infer where and when the settlement took place with great level of confidence, but not why and how.
These conclusions confirmed the theory of central land cores put forward for the medieval settlement of present-day Slovenia by Ilešič in the 1950s [77]. He noted that each medieval settlement initially had relatively little cultivated land on particularly favourable soils in the immediate vicinity of the settlement. As the settlement grew, the existing fields were divided up and new ones further from the village were asserted. Thus, the central land core became increasingly fragmented and the total area of cultivated land increased.
The theory of central land cores has good parallels with the site-catchment analysis proposed in the 1970s [78]. The similarities are not coincidental, as both are based on mid 20th century human geography. The site catchment is defined as an area within which the exploitation of natural resources is economically justified. The area was defined as 5 km or an hour's walk for sedentary farming communities and the share of arable land was estimated to be between 5% and 10%. Flannery [79], Rossman [80], and Zarky [81] empirically tested the model on Mesoamerican villages and found that the site catchment area was at least half and the share of arable land up to ten times smaller than in the original theoretical estimates. They concluded that the distance between villages was determined by social rather than ecological factors. Similar conclusion was reached for the early medieval Bled microregion, where the site catchment for the field was estimated to be 7 min walking distance [82]. Modern studies of the site catchment reinforce the distinction between the exploitation area and its social status, i.e., direct exploitation is not the same as the area that is claimed to define the political status of a settlement [83].
The key advantage of modern studies is that the catchment area is no longer forcefully simplified into circles but is much more realistically estimated in terms of time of walking or energy expended. This is achieved in GIS by computing the time distance based on DEM and realistic formulas obtained through experiments [84][85][86][87][88].

LiDAR Data
The airborne LiDAR data used in this study was acquired in 2014. These data have a nominal density of 5 points/m 2 and an estimated horizontal and vertical root mean square error of 0.09 m and are distributed via the eVode webservice [11,89] (for correlation between point cloud density and DEM quality see, e.g., [22]). The data were processed using an algorithm developed specifically for archaeology [22]. The relevant metadata and paradata have been presented elsewhere [31]. The main product used in this study is 0.5 m DEM with archaeology-specific off-terrain features included [22].
As already mentioned, in archaeology, processed LiDAR data are mostly used for interpretative mapping of archaeological features, i.e., feature detection. In this case study, however, we have used the data for what is termed integrated multi-scale 'deep' interpretation, which aims to deepen the understanding of archaeological features in their landscape context [31]. In this case, the digital terrain model is treated not just as a set of elevation values, but as an important habitat descriptor. The specific tools to achieve this are described below in more detail.

Geological Data
The Bled area can be divided broadly in four geomorphological areas: (i) the high alpine karst plateaus, (ii) the intramountain area, (iii) the till plain, and (iv) the marshy area. Jelovica (900-1411 m) were formed by glaciers in the Pleistocene (Appendix A: Figure A2). The area is composed of Middle Triassic dolomites and limestones. (ii) Sedimentary deposits on the Quaternary slope cover the intramountain area between Poljšica and Podhom, which slopes gently towards the alpine Radovna River valley and the glacial Lake Bled (lithostratigraphic unit al.-alluvium). (iii) The Bohinj and Radovna glaciers had a particularly strong influence on the geomorphology and postglacial fluvial processes, with strong glacial activity leading to the deposition of a till plain with up to several 10 m of Quaternary sediments and with a small marsh basin in the northeast part of Lake Bled (lithostratigraphic unit, pr.-till; b-marsh deposits). (iv) During glaciation, the marshy area between Lake Bled and River Rečica was formed.
A characteristic feature of the Bled landform is the frontal moraine on the northeast edge of the lake and the dome-shaped monadnocks rising above the general level of glacial deposits [90][91][92]. Bled Castle is located on one such cliff-like dome-shaped monadnocks [93]. The Pleistocene fluvioglacial sediments formed the terraces of Sava Dolinka and Sava Bohinjka Rivers (lithostratigraphic unit al.-alluvium and pr.-till). An important aftereffect of the underlying geological conditions in the study area is the lack of perennial water and permanent water streams (see Figures 1 and 3 for the locations mentioned in the text).
Of particular relevance to our case study is the overall glacial nature of the area, which is clear evidence that the geomorphology has not changed significantly since the Pleistocene, let alone since the beginning of the early medieval period.

Soil Conditions
The underlying lithology (bedrock) described above is one of the diagnostic criteria for the variety of soil types in the study area, which are briefly summarized here. Rendzinas formed on limestone, dolomite, moraines, and talus deposits. Dystric brown soils formed on carbonate and siliciclastic rocks. All are mostly suitable for forest and alpine pastures. Eutric brown soils formed on moraine and talus deposits and on fluvioglacial sandy gravel sediments. Small patches of rendzinas formed on limestone mostly support forests and meadows. Brown soils on fluvioglacial sandy gravel sediments are among the most fertile soils in subalpine areas. They occur in the plains, are well drained, sufficiently deep, and have favourable physical and chemical properties for intensive cropland. However, the brown soils formed on moraine and talus deposits are of limited use as arable land for modern agriculture, as the soil skeleton consists of moraine loam and stones. A notable depression with hydromorphic soils (hypogley) formed on a Pleistocene clay and loam northeast of Lake Bled; it is mostly suitable for grassland. The areas adjacent to the riverbeds of Sava Bohinjka and Sava Dolinka are dominated by undeveloped soils on alluvial river deposits that have been frequently flooded in the past. Suitable land uses here are riparian forests and grassland [94] (Appendix A: Figure A3).

Effective Field Capacity of Soil
For agricultural use, arguably the most important soil property is its ability to retain water. This quality is defined as the soil's effective field capacity (hereafter FC). FC depends on soil texture, depth, and organic matter content and is measured as the water content of a soil after gravity has drained as much water from the soil as possible [95]. The higher the FC value of a soil, the more water it is able to retain and the less susceptible it is to drought.
For mapping purposes, soil types are defined as discrete pedocartographic units, and FC is one of the criteria used. In the Soil map of Slovenia [94], which holds the best available data for the Bled microregion, FC is part of the description of pedocartographic units and is presented in 5 classes (Table 3; Appendix A: Figure A4). From the perspective of archaeology, the problem with soil maps is that they are produced on a small or medium scale. This is also the case with the soil map of Slovenia, which is designed for use at 1:25,000 scale, which is somewhat coarse for our purposes. To improve this, further analyses can be undertaken, such as the wetness index described below.

Modified Landform Classification Method
The landform or morphological classification of DEM, also termed geomorphology or morphometry, provides an objective and quantitative description of landform shapes, defined as specific geomorphic features, for example, plains, mountain ranges, hills, and valleys. The available methods have mostly been developed for geomorphological analysis of the terrain and are based on advanced spatial statistics [96][97][98]. We applied an automated landform classification method, topographic position index-based landform classification (hereafter TPI), implemented as a module in SAGA GIS (Hamburg, Germany) [99][100][101]. TPI provides a simple and powerful means of classifying the landscape into morphological classes. It is calculated as the difference between the elevation of a cell and the average elevation in large-and small-scale neighbourhoods. Positive values indicate that the cell is higher than its neighbours, while negative values indicate that the cell is lower (Figure 4) [97,101,102].
able data for the Bled microregion, FC is part of the description of pedocartographic units and is presented in 5 classes (Table 3; Appendix A: Figure A4). From the perspective of archaeology, the problem with soil maps is that they are produced on a small or medium scale. This is also the case with the soil map of Slovenia, which is designed for use at 1:25,000 scale, which is somewhat coarse for our purposes. To improve this, further analyses can be undertaken, such as the wetness index described below.

Modified Landform Classification Method
The landform or morphological classification of DEM, also termed geomorphology or morphometry, provides an objective and quantitative description of landform shapes, defined as specific geomorphic features, for example, plains, mountain ranges, hills, and valleys. The available methods have mostly been developed for geomorphological analysis of the terrain and are based on advanced spatial statistics [96][97][98]. We applied an automated landform classification method, topographic position index-based landform classification (hereafter TPI), implemented as a module in SAGA GIS (Hamburg, Germany) [99][100][101]. TPI provides a simple and powerful means of classifying the landscape into morphological classes. It is calculated as the difference between the elevation of a cell and the average elevation in large-and small-scale neighbourhoods. Positive values indicate that the cell is higher than its neighbours, while negative values indicate that the cell is lower (Figure 4)   TPI has proven to be one of the most important predictive variables for vegetation species distribution. For example, in a study of plant distribution in the Spring Mountains of Nevada (USA), TPI was second only to elevation as the most important predictive variable [102]. In other words, in a typical landscape, TPI classes are informative not only of TPI has proven to be one of the most important predictive variables for vegetation species distribution. For example, in a study of plant distribution in the Spring Mountains of Nevada (USA), TPI was second only to elevation as the most important predictive variable [102]. In other words, in a typical landscape, TPI classes are informative not only of landform classes but indirectly also of plant communities. This demonstrates the importance of TPI for all landscape-aware human decisions, including the choice of EMS location.
In our application to archaeology, the results of TPI have presented significant challenges to analysis ( Figure 5). The areas of moderately steep slopes and till plain were clearly defined, but the mountainous plateau and river terraces were not. Therefore, an additional visual geomorphological analysis was carried out. For this purpose, hypsometric tinting of DEM, transparently (60%) superimposed on a hillshade visualisation of the same DEM, was used to improve terrain classification and visualise relief differences more clearly. The most important criterion was the height above sea level. Applying this additional analytical step, we were able to precisely describe the mountainous plateau and the Holocene river terraces (Appendix A: Figure A5). clearly defined, but the mountainous plateau and river terraces were not. Therefore, an additional visual geomorphological analysis was carried out. For this purpose, hypsometric tinting of DEM, transparently (60%) superimposed on a hillshade visualisation of the same DEM, was used to improve terrain classification and visualise relief differences more clearly. The most important criterion was the height above sea level. Applying this additional analytical step, we were able to precisely describe the mountainous plateau and the Holocene river terraces (Appendix A: Figure A5). Our modified landform classification is thus a combination of TPI and visual geomorphological analysis that incorporates height above sea level. It allowed us to define Our modified landform classification is thus a combination of TPI and visual geomorphological analysis that incorporates height above sea level. It allowed us to define quantified catchment descriptors of landscape morphology, which we termed Zones. Defined in this way, Zones represent two of the most important predictive variables of plant species distribution: TPI and height above sea level [102].

Modified Wetness Index Method
Topographic modelling of soil moisture conditions can help alleviate the scale limitations of standard soil maps. Such modelling based on DEM is possible as water tends to flow and accumulate in response to gradients in gravitational potential energy [103]. The algorithms, commonly referred to as topographic wetness index, describe how susceptible specific areas in a study region are to become saturated [103,104]. They calculate for each cell of the grid the relationship between the specific upstream catchment area and the slope [105,106]. The first defines the potential of water intake (rainfall) and the latter the ability to discharge the water downslope (runoff; Formula: TWI = ln [Catchment Area/Slope]). One can think of these as a rainfall-runoff model (Figure 6a-c).
Topographic modelling of soil moisture conditions can help alleviate the scale limitations of standard soil maps. Such modelling based on DEM is possible as water tends to flow and accumulate in response to gradients in gravitational potential energy [103]. The algorithms, commonly referred to as topographic wetness index, describe how susceptible specific areas in a study region are to become saturated [103,104]. They calculate for each cell of the grid the relationship between the specific upstream catchment area and the slope [105,106]. The first defines the potential of water intake (rainfall) and the latter the ability to discharge the water downslope (runoff; Formula: TWI = ln [Catchment Area/Slope]). One can think of these as a rainfall-runoff model (Figure 6a-c). Methods differ primarily in the way the upslope contributing area is calculated [107]. We used the SAGA wetness index (hereafter SWI) because it does not think of the flow as a very thin film and hence it predicts more realistic (higher) potential soil moisture for valley floors [105]. The field tests demonstrated that SWI in combination with LiDAR derived DEM is the best existing predictor of soil wetness [103,108,109].
Another advantage of the SWI is that it can be refined by setting the suction index [110]. Unfortunately, the suction index function is poorly documented in the SAGA GIS software used, and the best available description is in the source code [111]. In addition, the suction cannot be adjusted locally. Therefore, we developed custom modified SWI (hereafter mSWI) by using the FC value extracted from the Soil map of Slovenia [94] as weighting index (Figure 6d). mSWI was calculated with map algebra using SWI and FC classes as an input.
In this way, we obtained mSWI (Appendix A: Figure A6) which combines the accuracy of the FC with the precision of the fine relief resolution of the SWI and is a very realistic predictor of soil quality. Methods differ primarily in the way the upslope contributing area is calculated [107]. We used the SAGA wetness index (hereafter SWI) because it does not think of the flow as a very thin film and hence it predicts more realistic (higher) potential soil moisture for valley floors [105]. The field tests demonstrated that SWI in combination with LiDAR derived DEM is the best existing predictor of soil wetness [103,108,109].
Another advantage of the SWI is that it can be refined by setting the suction index [110]. Unfortunately, the suction index function is poorly documented in the SAGA GIS software used, and the best available description is in the source code [111]. In addition, the suction cannot be adjusted locally. Therefore, we developed custom modified SWI (hereafter mSWI) by using the FC value extracted from the Soil map of Slovenia [94] as weighting index (Figure 6d). mSWI was calculated with map algebra using SWI and FC classes as an input.
In this way, we obtained mSWI (Appendix A: Figure A6) which combines the accuracy of the FC with the precision of the fine relief resolution of the SWI and is a very realistic predictor of soil quality.

General Methodological Remarks
There are two general methodological remarks to be made. First, our method of combining soil data with TPI and mSWI analysis is based on the premise that soil conditions in the Early Medieval period were similar to those of the modern period. This is justified in this particular case study by the fact that hydrological and surface conditions were subject to similar geomorphological processes throughout the Holocene and that the relationship between land surface properties (e.g., soil, vegetation, and lithology) was not very different in the early medieval period. In this particular case study, the stability is the result of the underlying lithology described above. Consequently, this method is only suitable for areas where either soil conditions have not changed significantly between the archaeological period under investigation and the time of soil data collection, or relevant soil data have been obtained through palaeoenvironmental analysis. In our case study, such example are urban areas where soil properties changed significantly (Figure 1: Zagorice, Želeče, Sp. Bohinjska Bela, Pristava), but those are relatively small and did not have significant influence on the results.
Second, the selection of methods used in this case study is indicative but by no means exhaustive. For example, slope and aspect can also be used as predictor variables for plant species distribution. In addition, climate (temperature, precipitation) and human impact are also very important for the distribution of plant species, as are many other factors. Alternative types of similar predictor variables include airborne LiDAR-derived feature detection used to identify landslides [112], spectral parameters of airborne LiDAR data applied for detection of glacial landforms [113], and object-based image analysis applied for volcanic and glacial landforms mapping [114]. TPI and mSWI methods in no way intend to compete with verified and established methods of environmental archaeology, such as archaeopalynology, archaeobotany, or archaeozoology, e.g., [115][116][117][118][119][120]. Rather, the aim is to introduce and test additional methods and, perhaps more importantly, to add LiDAR as a new data source for the archaeological analysis of past human land use. The suggested good practice would be to use TPI and mSWI in combination with other methods. However, in this case study, on the one hand, LiDAR and soil data are the only data currently available to the author, and on the other hand, TPI and mSWI were sufficient to provide new insights into the archaeological landscape in general and EMS in the context of agricultural land use in particular.

Results
Our modified landform classification is, as mentioned, the combination of TPI and visual geomorphological analysis, which resulted in the definition of four Zones. Below, each Zone is described (Figure 7; Table 4).
Zone 1 is defined as a mountainous plateau with steep and very steep slopes (TPI classes: High Ridges, Midslope Ridges, Local Ridges; 931-580 m a.s.l.). Middle Triassic dolomite and limestone bedrock prevail (Table 4: T2/1; T2/2) and two soil types occur. The first are rendzinas and the second dystric brown soils. The latter have a higher FC (FC index 3; mSWI index: 0, −5). Nowadays the area is forested and suitable for alpine pasture. There are no EMS in Zone 1.
Zone 2 consists of gently sloping terrain at the foothills. It occurs mostly in the western part of the study area, on the low hills surrounding the Lake Bled and above the river terraces (TPI classes: Upper Slopes, Open Slope; 580-510 m a.s.l.). The bedrock are mostly Holocene alluvial fan deposits. Prevailing eutric brown soils were formed on talus slopes mixed with moraine material and deposited directly on inactive alluvial fans [121]. These soils have high FC (FC index 3; mSWI index: 0, −5). The area is mostly suitable for arable land and meadows. 16 out of 19 EMS are located within Zone 2.
Zone 3 represents a large till plain formed in postglacial fluvial processes that deposited up to several tens of metres of Quaternary sediments. It is limited by the riverbeds of Sava Dolinka and Sava Bohinjka (TPI class Plains; 510-480 m a.s.l). Over the glaciofluvial sand gravel deposits (Table 4: pr), fertile deposits of brown soils developed. However, due to the high porosity of Holocene sediments, the FC is low (FC index 2; mSWI: 0.7), which means that the entire area is exposed to drought. This is exacerbated by the absence of permanent surface water. Nevertheless, there are small patches of hydromorphic soils (Alluvial soils, Hypogley, Amphigley) with high FC (FC index 4; mSWI: 0, −5). Their formation was possible due the glacial activity and postglacial fluvial processes, which have resulted in deposition of clayed sediments north of the Bled lake [92]. The brown soils in the Zone 3 are the most suitable soils for modern agriculture in the area [94], providing that drought effect can be mitigated (for example, by irrigation or drought-resistant crops). Only three EMS, all established only in the eleventh century, are located in Zone 3.   It can be concluded that the preferred landscape type for EMS was moderately steep slopes and brown soils with high FC, defined here as Zone 2 (Figure 7). This is the case for most EMS in our case study (Figure 1: Višelnica, Zgornje in Spodnje Gorje, Poljšica, Grmišče, Zasip, and Mužje). The location of two other EMS (Figure 1: Zg. in Sp. Bodešče) fits the landform classification criteria but not the soil conditions as depicted on the pedological map. We explain this by the fact that the existing soil map is not detailed enough to show the microlevel differences. Indeed, the area is full of glacial moraines and micro valleys, and under such conditions, water-rich and marshy soils tend to develop. Their presence in this particular area is confirmed by the historical field names ("V blateh", "Curkovca", "Pretaka", "Nad potokam", which means "In the mud", "Stream", Flow", "Above the stream") [37].
The only other landscape context where three EMS exist is large till plain with fertile brown soils with low FC, defined here as Zone 3. However, all three (Figure 1: Zagorice, Grad 3, Koritno) have only been established in the eleventh century.
The above presented focus of EMS on a landscape characterised by moderately steep slopes and brown soils with high FC is consistent with previous research on EMS in similar landscape conditions by Wawruschka. Her mountainous or hilly areas fit well with the description of our Zone 2, although some of the data (e.g., m a.s.l.) cannot be directly compared [38].
The most important result of this analysis is the definition of the ecological niche that was preferred by the EMS and is based on the agricultural land use. The importance of this lies in the scalability, i.e., this result can be directly applied to regional studies of the early medieval settlement in Eastern Alpine region and possibly other regions with subalpine climate.
The results also enable new insights into the early medieval Bled microregion by characterizing the individual EMS. Exclusive preference for Zone 2 prior to the eleventh century strongly suggests two things. First, these are primarily agricultural settlements. There are two exceptions (Figure 1: Grad 2, Mlino) where the landscape morphology does not allow for the presence of significant arable land and non-agricultural function seems probable [37]. Second, the relatively narrow scope of agricultural land use, as can be inferred from the exclusive occupancy of Zone 2, suggests a not overly diversified agricultural land use system, possibly based on a single staple crop.

Discussion
The significance of these results goes beyond the case study area and is important for the EMS of the entire subalpine area settled by the Slavs. To fully appreciate the advantages of the proposed method, however, a brief discussion of the archaeological context must be introduced.
Brown soils with high FC (Zone 2) in the vicinity of which most of the EMS are located are not particularly suitable for agriculture by high-and late medieval or post-medieval standards because they contain abundant moraine loam enriched with stones. This is particularly conspicuous because often more fertile soils which are easier to work but have low FC (Zone 3) are adjacent ( Figure 5). The late expansion into "the area in the plain where the soil was good" (our Zone 3) was already observed by Pleterski; based on historical and archaeological analysis he concluded that the area was covered with forest but could not explain why [37]. Then, why were areas of high FC (Zone 2) the main attractor of EMS? The land use characteristics presented above suggest that the key factor in the EMS location choice model was the proximity of the specific soil type that minimizes the risk of exposing crops to water stress. This conclusion is noteworthy for two reasons. First, it helps to elucidate the agricultural characteristics of EMS. Second, it can be used to predict the location of landscape contexts suitable for EMS throughout the Eastern Alpine region.
While the latter is the subject of ongoing research, the former can be discussed here. A brief overview reveals that the early medieval Slavs had a limited choice of the principal field crops available. Rye (Secale cereale), wheat (Triticum sp.), oats (Avena sp.), barley (Hordeum sp.), and millet (Panicum miliaceum) were the main cereals in the western Slavic settlement area (present-day eastern Germany) [122]. In Southern Russia we know of wheat and barley [123]. The first Slavic settlers in north-western Russia brought with them a great variety of cereals and legumes, but only the crops that guaranteed agricultural success in the colder north were kept in cultivation: barley and rye [124]. Closer to the Bled microregion, wheat, barley, and rye were the most common crops at Roztoky (Czech Republic) [125]. In Thunau (Austria), both in the settlement and in one of the graves, the most common cereal remains were wheat, millet, barley, and rye [126]. In Kleinklein (Austria), the contents of a settlement pit revealed barley, millet, and rye [127]. On Bled Island, the central location of the Bled microregion, charred barley and millet grains were recovered but are dated to a later period in the mid-13th century [128]. Therefore, the early medieval Slavs, like other contemporary Europeans, relied mainly on wheat, barley, rye, and millet as staple cereal foods. Of these, barley has the greatest ecological amplitude and is able to cope with extreme ecological conditions [129]. It was grown, for example, in Highland Britain [130], in Scandinavian north-western Europe [131], and even in the Faroe Islands [132]. Regardless of climate, barley was the dominant crop in western Europe and Britain at the beginning of the early Middle Ages [122], where it was important enough to warrant a special barley tax [133]. Wheat in western and rye in north-western Europe had replaced barley as the dominant cereal by the end of the early Middle Ages [131] ( Figure 8). and is able to cope with extreme ecological conditions [129]. It was grown, for example, in Highland Britain [130], in Scandinavian north-western Europe [131], and even in the Faroe Islands [132]. Regardless of climate, barley was the dominant crop in western Europe and Britain at the beginning of the early Middle Ages [122], where it was important enough to warrant a special barley tax [133]. Wheat in western and rye in north-western Europe had replaced barley as the dominant cereal by the end of the early Middle Ages [131] (Figure 8). There are many differences between wheat, millet, barley, and rye. Most pertinent to this case study is that while millet and rye have exceptional drought tolerance, wheat and barley under rainfed conditions may suffer from drought resulting in significant yield loss [134,135]. Consequently, it is barley and wheat that require soils with high FC. Between There are many differences between wheat, millet, barley, and rye. Most pertinent to this case study is that while millet and rye have exceptional drought tolerance, wheat and barley under rainfed conditions may suffer from drought resulting in significant yield loss [134,135]. Consequently, it is barley and wheat that require soils with high FC. Between the two, barley matures earlier, has greater tillering capacity, and competes better with weeds, but generally yields less [136]. Wheat requires more labour in terms of cultivation and manuring to achieve potentially high yields. An increase in wheat cultivation thus implies a higher degree of agricultural organization [131,137].
Therefore, we infer that barley was most likely the principal field crop in the Bled microregion until the eleventh century. There are three pieces of evidence to support this. First, barley is best suited to a society with a lower level of agricultural organization [131,137], such as the colonist society that was establishing EMS in the Bled microregion between the 7th and 9th centuries (cf., [138]; also see below). Second, barley was the cereal of choice for the Slavs when they settled in what they perceived as colder climates, such as the abovementioned colonization of north-western Russia [123]; due to its height above sea level close to 500 m the subalpine climate of the Bled microregion is colder (under any climatic conditions) than the areas from which the Slavs were arriving [64], for example, from the western edges of the Pannonian plain that they settled already in the 6th century [25]. Third, under rain-fed conditions, barley prefers high FC to all other soil properties [134,135].
To explain the expansion of EMS in the eleventh century into areas of soil with lower FC (Zone 3), another evidence must be introduced. Namely, in 1004 AD, parts of the Bled microregion were bestowed to the bishops of Brixen by a royal charter [139,140]. In addition to the small estate, they received all the lands in the Bled microregion that were not directly worked by existing owners, including the then forested Zone 3 [37]. This deed was more than just an exchange of ownership, as it signified the assertion of direct control over the Bled microregion by the Kingdom of Germany for the first time after more than a century (for historical context see, e.g., [141,142]). The new owners undoubtedly imposed their own agricultural organization on their newly acquired property. It was based on the complex manorial system, which was a subsistence economy oriented towards stability and based on strategies of risk avoidance through diversification of resources and redistribution by means of storage and transport [143]. In the early eleventh century (i.e., at the end of early medieval period), wheat was the main staple crop in this system [131]. Subsequent institutional and social changes in the Bled microregion that reflected a higher degree of agricultural organization have been chronicled as extensive property changes in the following centuries [37,139,140,144]. This was a marked difference from the pre-eleventhcentury system operated by the Slavs, which used similar technology [56], but was far less diversified and almost exclusively embedded in the local community of a single Župa [37,138], which in this case is the Bled microregion.
The political and social changes of the early 11th century thus introduced a higher degree of agricultural organization that is necessary, among other things, for a successful wheat production [131,137]. More importantly, the EMS expansion into areas with low FC soils was possible since the diversified manorial system was designed to routinely cope with local crop failures by transporting supplies from distant estates [143]. In other words, the pre-eleventh-century agricultural system could not cope with occasional crop failures on low FC soils due to dry Summers. The post-eleventh-century manorial system could and was therefore able to take advantage of this otherwise fertile area, most likely to grow wheat. The fact that none of the eleventh century EMS were established directly by the bishops of Brixen [37] does not change the fact that it was the newly introduced manorial agricultural system that enabled the colonisation of this area.
These conclusions are important not only for the understanding of the Bled microregion EMS but also to further our understanding (e.g., [145,146]) of regional agricultural dynamics during the transition from the early to the high medieval period in Eastern Alpine region.

Conclusions
The main objective of the paper was to present and test new methods that offer additional uses of LiDAR for archaeological analysis of past land use. This is important not only in terms of the specific method and research application presented but also because this is an underdeveloped area in archaeology that needs more attention.
This paper utilizes an existing corpus of open access archaeological database Zbiva [73,147], open access remote sensing data and environmental data (geological and soil data), as well as open-source software tools (e.g., QGIS, SAGA) to reassess existing knowledge on the early medieval archaeological landscapes, specifically on agricultural land use. While the importance of free and open-source software in science in general, e.g., [148], and in the field of airborne LiDAR data for archaeology in particular, e.g., [149], is well recognised, we believe that the importance of the increasingly abundant and easily accessible free environmental and archaeological data, e.g., [150], is too often overlooked. Hopefully, this article is a step towards recognizing the importance that this data source can have for archaeology.
A novel objective method was presented and, perhaps more importantly, LiDAR as a new data source for the archaeological analysis of agricultural land use systems. The suggested good practice would be to use this method in combination with existing complementary methods, such as archaeobotanical analyses. However, in this case study, the presented method was sufficient to provide new insights into the archaeological landscape in general and EMS in the context of agricultural land use in particular.
We used the LiDAR data for what is termed integrated multi-scale 'deep' interpretation, which aims to deepen the understanding of archaeological features in their landscape context. It should be reiterated that, in our opinion, such a use of these data in archaeology remains underexploited despite some promising early studies, e.g., [18,21,28]. The Bled case study illustrates such potential contribution of LiDAR data to explore landscape gradients that have influenced human activities. We have clearly demonstrated a preference of early medieval agriculture for terrain on moderately steep slopes with brown soils that have a high capacity to retain water. We have tentatively connected this to a society with a subsistence agricultural organization that is growing barley as a staple cereal food. It is highly likely that cultural and technological changes at the turn of the millennium brought, among other things, a higher level of agricultural organization and wheat as a staple cereal food.
One of the most important contributions is the discussion of scale issues. Since the scale of many soil maps is inadequate for archaeological analysis, a method to overcome this challenge is presented using various indices. The solution presented is scalable to other types of landscape and other archaeological periods, as well as to other types of soil data.
In the wider context of LiDAR methodology in archaeology, we have focused on the potential of LiDAR data to provide a source for very detailed landscape description and observe environmental components using GIS analysis, specifically modified landform classification and mSWI. This approach leads to a more detailed and objective analysis of the environment and spatial context of any observed archaeological phenomena. Given the rise of open access data and open access tools, there is huge potential for this and similar methods in geocomputational archaeology of the near future.
Funding: This research was funded by Austrian Science Fund (FWF) grant number I 3992 (Internationale Projekte, project title "Settlement of the south-eastern Alpine region in the Early Middle Ages").

Appendix A. Additional Maps
Acknowledgments: I thank Mateja Belak for providing the raw data from the Zbiva database and Benjamin Štular and Jernej Rihter for discussion and suggested readings for the Discussion section. Open Access funding by the Austrian Science Fund (FWF).

Conflicts of Interest:
The author declares no conflict of interest. The funder had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results. Figure A1. Bled microregion (decimal longitude and latitude coordinates of the map centre: 14.1139; 46.3752), early medieval settlements (letters refer to Table 2 in the text), and arable land category "good" based on the retrograde analysis of the 19th century Franciscan Cadastre (source data adopted from [37]; North at the top).  Table 2 in the text), and arable land category "good" based on the retrograde analysis of the 19th century Franciscan Cadastre (source data adopted from [37]; North at the top).  (source data adopted from [94]; North at the top). Figure A3. Bled micro-region (decimal longitude and latitude coordinates of the map centre: 14.1139; 46.3752), soil map (source data adopted from [94]; North at the top).  Figure A4. Bled micro-region (decimal longitude and latitude coordinates of the map centre: 14.1139; 46.3752), effective field soil capacity (FC) classes (source data adopted from [94]; North at the top).