ALS-Based Detection of Past Human Activities in the Białowie ˙za Forest—New Evidence of Unknown Remains of Past Agricultural Systems

: The Białowie˙za Forest (BF), a unique ecosystem of historical signiﬁcance in central Europe, has a long history of assumed human settlement, with at least 200 known archaeological sites (until 2016). This study uncovers new evidence of the cultural heritage of this unique forest area using Airborne Laser Scanning (ALS) technology combined with traditional archaeological ﬁeld assessment methods to verify the ALS data interpretations and to provide additional evidence about the function and origin of the newly detected archaeological sites. The results of this study include (1) a scientiﬁc approach for an improved identiﬁcation of archaeological resources in forest areas; (2) new evidence about the history of the human use of the BF based on ALS data, covering the entire Polish part of the BF; and (3) an improved remote sensing infrastructure, supporting existing GIS (Geographic Information System) systems for the BF, a famous UNESCO Heritage site. Our study identiﬁed numerous locations with evidence of past human agricultural activities known in the literature as “ﬁeld systems”, “lynchets” and “Celtic ﬁelds”. The initial identiﬁcation included more than 300 km of possible ﬁeld boundaries and plough headlands, many of which we have veriﬁed on the ground. those 13th The results of this study demonstrate that past activities in Polish part of the had been than As a practical of the a geodatabase was created; this has practical applications for the system of monument protection in Poland, as well as for local communities and the BF’s management and conservation. The more widely achieved results are in line with the implementation of the concept of a cultural heritage inventory in forested and protected areas—the actions taken specify (built globally) the forms of protection and management of cultural and environmental goods.

Remote Sens. 2020, 12, 2657 3 of 23 from the 5th century BC to the beginning of the 1st century AD) and Roman influence (1st century AD to approximately the 4th century AD)) in this area represent the foundations of the culture of barium ceramics, evidenced by either the Przeworsk culture or the Wielbark culture. In ancient times, the BF was used in a variety of ways: for hay-mowing, bar-planting, foraging, fishing, firing tar and potash production, as well as for grazing of animals, hunting, gathering, mining, honey production, mushroom collecting and finally forest management of various kinds, including logging [20][21][22]. Most of the remnants of human activities are cemeteries or single graves [23], including cremation graves from the 8th to the 9th century and inhumation graves from the 11th to the 12th century. From the 14th to the end of the 18th century, these areas were the royal hunting ground of the Jagiellonian dynasty [24]. As a result, even if the exploitation of BF has varied during the past few centuries, it is still one of the few lowland forests in Europe that has never been completely cleared of forest cover [25].

Study Area
The Białowieża Forest, a remnant of the Central European lowland forests that once occupied much of the central European lowland plains [47], covers approximately 150,000 ha with 42% of its area located in Poland (Figure 1). The Polish part of the BF (approximately 63,000 ha) comprises four distinct spatial entities: the Białowieża National Park (BNP), 10,517 ha; the Białowieża Forest District, 12,593 ha; the Browsk Forest District, 20,417 ha; and the Hajnówka Forest District, 19,654 ha. In 1979, the BNP area was included in the United Nations Organization for Education, Science and Culture (UNESCO) World Heritage List. The area of the UNESCO site was extended in 2014, and now covers an area of 141,885 ha, with a buffer zone of 166,708 ha [48], including the Polish and Belarusian parts of the BF.

General Concept
The concept of the inventory of archaeological heritage in the BF area is closely related to a broader undertaking-with the natural and cultural inventories of the State Forests National Forest Holding in Poland-which is an initiative to create digital resources like spatial data for other forest areas in Poland. Since 2016, these works have also been carried out in south-eastern Poland-mainly In the present study, the first of its kind on such a large scale in the BF, we collected data for the entire area by utilizing Light Detection and Ranging (LiDAR) technology mounted on an aircraft, also known as an Airborne Laser Scanning (ALS) system. In many recent archaeological projects, ALS was one of the leading methods [26][27][28][29][30][31][32][33][34][35]. In Europe, numerous archaeological projects use public ALS data [13], but their utility for archaeological purposes was often low (approximately 1 pts/m 2 (points per square metre) or less), which means that other solutions, such as those dedicated to Remote Sens. 2020, 12, 2657 4 of 23 archaeology but also for natural research data acquisition, are necessary. It was no different in the Polish part of the Białowieża Forest, where two datasets were available. There are ALS data from the ISOK (Informatyczny System Osłony Kraju przed nadzwyczjnymi zagrożeniami-the computer system covers the country against extraordinary threats) project, covering almost the total area of Poland [36,37]. Only a part of the BF was available. However, in 2015, additional ALS data with higher density were obtained from the LIFE+ ForBioSensing project (www.forbiosensing.pl). The ForBioSensing ALS dataset was of much better quality and covers the whole BF, which is why it was used in the presented research.
ALS is an active remote sensing technology that gives users the ability of surveying large areas with high precision. ALS is capable of obtaining ground echoes even under forest canopies [38]. It was introduced to the archaeological community at the beginning of the 21st century and became an essential tool in the recognition of past human activity remains. There are many successful examples of using ALS data for inventories of archaeological objects under forest cover [7,39,40]. ALS data are exceptionally suitable for producing a Digital Terrain Model (DTM) of the ground under forest cover [41], and because of this functionality, ALS technology complements archaeological research in forested areas [42,43]. The DTM interpolated from ALS is extremely useful for identifying small-scale topographic anomalies [44]. Based on local changes, areas with probable anthropogenic origin are determined [45]. Visual DTM interpretation requires knowledge-based pattern recognition skills and the ability to identify and classify complex landforms based on experience and previous archaeological knowledge [46]. After detecting possible human-made features based on DTM interpretation, the field part of the research is carried out. In the field, previously marked areas are verified to confirm their actual origin and condition. Despite this complexity, it is possible to pinpoint potential human activity residues quite precisely, even in forest areas.
Based on the experience of other scientists worldwide, as well as the availability of ALS data, it was decided to verify whether such good results can also be achieved for the Polish part of the BF. The aim of this study was to detect the remains of past human activities and verify them in the field in the whole Polish part of the BF. The main output was recognized as a base for future deep studies of all human activity remains in order to verify the scale and degree of human influence on the BF. In this study, we are focusing on the identification of linear objects, as well as their detection and verification.

Study Area
The Białowieża Forest, a remnant of the Central European lowland forests that once occupied much of the central European lowland plains [47], covers approximately 150,000 ha with 42% of its area located in Poland ( Figure 1). The Polish part of the BF (approximately 63,000 ha) comprises four distinct spatial entities: the Białowieża National Park (BNP), 10,517 ha; the Białowieża Forest District, 12,593 ha; the Browsk Forest District, 20,417 ha; and the Hajnówka Forest District, 19,654 ha. In 1979, the BNP area was included in the United Nations Organization for Education, Science and Culture (UNESCO) World Heritage List. The area of the UNESCO site was extended in 2014, and now covers an area of 141,885 ha, with a buffer zone of 166,708 ha [48], including the Polish and Belarusian parts of the BF.

General Concept
The concept of the inventory of archaeological heritage in the BF area is closely related to a broader undertaking-with the natural and cultural inventories of the State Forests National Forest Holding in Poland-which is an initiative to create digital resources like spatial data for other forest areas in Poland. Since 2016, these works have also been carried out in south-eastern Poland-mainly the Remote Sens. 2020, 12, 2657 5 of 23 area of the Regional Directorate of State Forests in Krosno that largely coincides with the area of the Podkarpackie Voivodeship.
Starting in 2016, the initiative to conduct extensive research in the study area towards the identification of archaeological (cultural) resources has been a large-scale interdisciplinary effort aimed at investigating various forms of evidence for past human activities; it includes the recognition, registration and protection of remains, as well as the public dissemination of information about these findings. The ongoing identification of archaeological heritage resources, including previously unknown objects and relics of field systems, is based on multi-element and multi-stage research procedures applied to the polish part of the Białowieża Forest. The procedures consist of (1a) acquisition and processing of archival and remote sensing data-mainly ALS-data and aerial images; (1b) analysis and interpretation of the data's correlations with existing research results from 2016 and with the Białowieża Forest's current parameters, such as contemporary topography and activities; (1c) analysis, interpretation and verification of the data about the identified objects during and after the field inspection of those objects; (2) multi-site field verifications of point probes and wells, and field surveys and excavations; and (3) specialized analyses, such as dating and analysis of the chemical soil composition, ecofacts and artefacts.

Archival Data
The analysis of archival resources included literature and other archival materials (e.g., aerial photos from World War II, archival maps) and documentation, including resources from the Provincial Office for the Protection of Monuments in Białystok, the Institute of Archaeology and Ethnology of the Polish Academy of Sciences and the Podlasie Museum in Białystok. In addition, archival resources acquired as part of the activities of the Forest Culture Centre in Gołuchów were analysed. This work involved the correlation of features of known archaeological objects (shape, size and location) with those of the newly recognized anthropogenic objects. Data related to land use, along with data on the occurrence and condition of potential archaeological resources, were also collected. The result of these activities is compiled into a database of 200 archaeological sites ( Figure 1). Correlation with historical (archival) data indicates the probability of the destruction of some sites or single archaeological features. This can happen where, for forest renewal purposes, the top layer of soil is subjected to mechanical preparation for planting, or due to other activities or environmental processes.

ALS Data
Primary recognition of the potential archaeological objects was based on three critical elements associated with the ALS data: (i) the ALS data acquisition; (ii) data processing; and (iii) data interpretation. Raster data for the visual interpretation of potential archaeological objects in the Białowieża Forest were generated based on ALS data obtained in the leaf-off season between 25 November and 7 December 2015. The ALS data were collected with a full-waveform sensor Riegl LMS-Q680i scanner, with an average point density above 11 pts/m 2 at a flying height of 500 m above ground. To cover the study area, 135 individual flight paths were flown with a 40% overlap. The horizontal accuracy was 0.20 m and the vertical accuracy was 0.15 m. The point cloud was automatically classified into five classes using TerraSolid software based on the ASPRS 1.2 LAS format specification [49]: Positioning System-Real Time Kinematic (GPS-RTK) method. The DTM constituted the base for generating the derived products used in the inventory.

DTM Processing
Recognition of potential archaeological objects was based on different visualizations generated from the DTM [44,50]. These products were produced using the Relief Visualization Toolbox [42]. The following results of the ALS-based DTM processing were used in the interpretation (main visualizations): We used the open-source software QGIS (version 2.14) to standardize the inventory results of all the interpreters' work for viewing data and searching objects of interest. ESRI Shapefiles were created to store the final version of the interpretation. Uniform schematics made it easier to combine data from different interpreters, e.g., for the entire forest district. For each catalogued point of interest verified by an archaeologist, a part of a point cloud, raster image, terrain profile and 3D view were created to show the best representation of the terrain. We used the Profile Tool in the QGIS Toolbox to generate the elevation profile for the selected section. The basic inventory unit was the forest compartment. We used a convex polygon tool to calculate the area of the groups of linear objects, which were clearly detected and formed closed structures. Convex polygons create the smallest possible border around the marked object/clusters in our case. Afterwards, the area calculation was based on the borders around separate groups of detected linear objects.
In total, 2108 forest compartments were analysed in the whole research area, including 444 compartments in the Białowieża Forest District, 396 in Białowieża National Park, 598 in the Browsk Forest District and 670 in the Hajnówka Forest District. The average size of the forest compartment in the BF is approximately 3 ha. Military areas were excluded from the inventory.

Primary Analysis and Interpretation
The main focus of the analysis and interpretation of the ALS data was on the visual recognition of topographic differences, their locations, and their correspondence with other elements of the environment and with contemporary objects, as well as the correlation of their shapes, sizes and complexities with those of the objects they are found with in conjunction. The interpreters used many supplementary data to support their work, including The next stage of the inventory processing was to control the data accuracy and the integrity of the prepared shapefiles and processed attributes. The merged and validated file for the entire Forest District was passed to the archaeologist for verification. Figure 2 shows the workflow diagram for this process. Two approaches to analyse and interpret the ALS data were used in the research to identify cultural resources: (1) office works and (2) office-terrain works. The first one was based solely on Remote Sens. 2020, 12, 2657 7 of 23 the analysis of data in the form of "from the desk", and the second was based on the analysis and interpretation of data both as part of office and field work (beginning with terrain verification).

Figure 2.
General workflow with ALS data acquisition, processing and interpretation, as well as the subsequent field-based verification.

Preparations for Field Verification
Figures 3 and 4 show the locations of potential archaeological objects (Section 3.1) indicated by the office team of interpreters in the Polish part of the Białowieża Forest. Based on layers with potential archaeological objects corrected by an archaeologist, cartographic materials and illustrations were prepared for the operational teams to find ALS-based identifications in the field. The review materials were created from the hillshade, local relief model, and other raster transformations. For proper interpretation of the data in the field, it was assumed that the scale of the visual materials should be 1:1000 for point objects, and 1:2500 for linear objects. Materials were saved in georeferenced TIFF format in two versions, with vector objects marked and unmarked. Such data can be displayed on the backdrop of other cartographic and vector materials and is also easier to use on tablets with built-in GPS when locating objects on the ground. For archaeologists working in the field, georeferenced maps in TIFF or JPG format were made that contained detailed locations of the forest compartments on a contour map. For preparing field-based materials, additional data were used, such as

Preparations for Field Verification
Figures 3 and 4 show the locations of potential archaeological objects (Section 3.1) indicated by the office team of interpreters in the Polish part of the Białowieża Forest. Based on layers with potential archaeological objects corrected by an archaeologist, cartographic materials and illustrations were prepared for the operational teams to find ALS-based identifications in the field. The review materials were created from the hillshade, local relief model, and other raster transformations. For proper interpretation of the data in the field, it was assumed that the scale of the visual materials should be 1:1000 for point objects, and 1:2500 for linear objects. Materials were saved in georeferenced TIFF format in two versions, with vector objects marked and unmarked. Such data can be displayed on the backdrop of other cartographic and vector materials and is also easier to use on tablets with built-in GPS when locating objects on the ground. For archaeologists working in the field, georeferenced maps in TIFF or JPG format were made that contained detailed locations of the forest compartments on a contour map. For preparing field-based materials, additional data were used, such as

Systematic and Multi-Site Field Verification and Specialist Analyses
The systematic and multi-site field verification and specialist analyses of the anthropogenic remains of agrarian systems consisted of the following steps: (1) field surveys; (2) geological drilling and/or soil coring; (3) excavations surveys; (4) soil pits; and (5) specialist analyses/laboratory analysis. As part of the fieldwork, a comparison of LiDAR data with the situation in the field was conducted.

Primary Recognition
As a result of the primary recognition process (Section 2.6) conducted over the Polish part of the Białowieża Forest, more than 20,000 potential archaeological features (single archaeological objects-features such as barrows, etc.) or sites (for example several features in a group) were identified ( Figure 3). In addition, approximately 4000 so-called "potential linear objects" were found ( Figure 4)-and were mostly field system remnants. Table 1 lists the potential archaeological objects identified in each administrative unit by the team of interpreters. We cannot treat these first

Systematic and Multi-Site Field Verification and Specialist Analyses
The systematic and multi-site field verification and specialist analyses of the anthropogenic remains of agrarian systems consisted of the following steps: (1) field surveys; (2) geological drilling and/or soil coring; (3) excavations surveys; (4) soil pits; and (5) specialist analyses/laboratory analysis. As part of the fieldwork, a comparison of LiDAR data with the situation in the field was conducted.

Primary Recognition
As a result of the primary recognition process (Section 2.6) conducted over the Polish part of the Białowieża Forest, more than 20,000 potential archaeological features (single archaeological objects-features such as barrows, etc.) or sites (for example several features in a group) were identified ( Figure 3). In addition, approximately 4000 so-called "potential linear objects" were found ( Figure 4)-and were mostly field system remnants. Table 1 lists the potential archaeological objects identified in each administrative unit by the team of interpreters. We cannot treat these first recognition results as archaeological remains. There are potential archaeological sites that were subject to additional field verification. *-Potential archaeological features (such as barrows, charcoal kilns, etc.) or contemporary anthropogenic objects or unspecified objects). **-Potential archaeological features or contemporary anthropogenic objects or unspecified objects with their own landscape sculpture (with different widths of elevation, ridges or faults-especially potential remains of field boundaries (among other banks and lynchets demarcating the fields)).
Among these, clusters-groups of linear ridges-were specified, defining the polygonal (usually irregular quadrilateral) surfaces (Figure 4). Similar arrangements are known from the literature and are identified with ancient agricultural systems-"field systems", "lynchets" or "Celtic fields" [53][54][55][56]-whose existence has also been confirmed by other research teams for the Białowieża Forest [57]. A similar diagnosis of old agricultural systems (and similar elements of this type of construction, e.g., copper) is shown in numerous studies, including those that use ALS data [58][59][60][61][62][63]. There were at least two types of systems observed with the linear elevations: (1) with an irregular concentric system and (2) with more regularly arranged forms. Numerous objects occurred separately, or apart from clusters, which in many situations could be associated with the destruction of more substantial deposits. Some of these objects are probably the remains of otherwise undetected human activities in the distant past.

Fieldwork Verification
During 2016 [19,64,65] and especially during 2017-2019 [66][67][68], we carried out an extensive campaign of field work, including archaeological surveys, test-pits and excavations. Detailed results of the abovementioned research will be published in separate publications-in this manuscript we are concentrating only on linear ridges, or remains of field systems. The aim of the field verification work was to confirm the origin and role of the potential archaeological sites detected on the ALS-based DTM.
The most important group of objects was related to clusters of linear elevations (ridges) purely seen in the field ( Figure 5), defining the polygonal (usually irregular quadrilateral) surfaces (Figures 4  and 5). More than 300 km of linear ridges, individual or grouped into polygonal surfaces (or occurring outside of groups), could be clearly identified based on the ALS data ( Figure 5B,D). The objects are mainly present on local hills (raised areas). This can be connected to agricultural land-use patterns since these features are not found within floodplain zones-above the flood areas ( Figure 5A). Linear ridge systems cover a large area of the BF. Based on the ALS data, the total area of clusters is over 2500 ha (Figure 4), which is approximately 4% of the entire Polish part of the Białowieża Forest. The area of object complexes in the Białowieża National Park is approximately 1500 ha, which represents approximately 14.3% of its total area.Numerous objects occurred separately, or apart from clusters, which in many situations could be associated with the destruction of more substantial deposits.
The verification concentrated on selected remains of the field systems based on archaeological surveys ( Table 2). These were identified on the ground primarily based on small linear elevations-ridges from a few to several dozen centimetres in height with widths of up to a few metres and lengths often exceeding 100 m (Figures 5-7). A total of 1232 out of 1328 of the linear objects interpreted as linear ridges were positively verified (potential remains of banks or lynchets demarcating the fields) by recognizing and documenting the course of the objects in the field, also in the area of the BPN [67]. Remote Sens. 2020, 12, x FOR PEER REVIEW 11 of 24 The verification concentrated on selected remains of the field systems based on archaeological surveys ( Table 2). These were identified on the ground primarily based on small linear elevations-ridges from a few to several dozen centimetres in height with widths of up to a few metres and lengths often exceeding 100 m (Figures 5-7). A total of 1232 out of 1328 of the linear objects interpreted as linear ridges were positively verified (potential remains of banks or lynchets demarcating the fields) by recognizing and documenting the course of the objects in the field, also in the area of the BPN [67].    In many places, the readability of the objects was poor, which affected their exact field documentation. This situation is related to the destruction of objects either by human activities or by environmental processes. In many places, clusters of stones were lying on the ground surface, either directly on the line of a ridge or in its close vicinity, or even moved from their original locations ( Figure 8). Field observations of existing trees indicated that these objects had been created before the extension of the forest in those areas. In many places, the readability of the objects was poor, which affected their exact field documentation. This situation is related to the destruction of objects either by human activities or by environmental processes. In many places, clusters of stones were lying on the ground surface, either directly on the line of a ridge or in its close vicinity, or even moved from their original locations ( Figure 8). Field observations of existing trees indicated that these objects had been created before the extension of the forest in those areas.   A detailed presentation of the results of the field is the subject of a separate study in progress. Below we summarize the general results.
Soil pits, which revealed the stratigraphic positions of individual objects, provided materials for special analyses. The first surveys and excavations were carried out in 2016 [19,64,65] Additional activities (mainly from 2017), apart from field verification, related to the recognition of linear objects consisted of A detailed presentation of the results of the field is the subject of a separate study in progress. Below we summarize the general results.
Soil pits, which revealed the stratigraphic positions of individual objects, provided materials for special analyses. The first surveys and excavations were carried out in 2016 [19,64,65]   As part of the invasive intensive fieldwork, 91 test pits were made in transects (3 pits in each transect) in order to carry out physico-chemical soil analyses. Volumetric analysis, aimed at determining differences in the degree of soil compaction, and the measurement of the share and vertical distribution of carbon and nitrogen in soils up to 70 cm deep, were used as indicators of soil transformation under the influence of cultivation or fertilization. Out of the 66 pits from 2017, more than half had confirmed significant changes in the vertical distribution of carbon and nitrogen in relation to soils with confirmed natural regeneration and forest management origins [69].
In contrast with the literature on the subject, remains of former fields can be initially associated mainly with the period of Roman influence. The settlement context, recognizing archaeological sites and other factors, also indicate a chronology within the pre-Roman period from the 5th century BC to about the 1st century AD (e.g., Figure 4D). Recognized field systems are dated mainly from about the 1st/2nd century to the 4th century AD. These field systems can be identified as associated with the activities of the Wielbark culture population, while the objects can be associated with the chronologically younger Slavic communities from the period of the 6th-13th century AD. In relation to some of these linear ridges and polygonal spaces, one should also consider the remains of modern constructions. A detailed presentation of the results of field and laboratory tests is the subject of a separate study (in progress). The results of the work are also the basis for further work on the reconstruction of former agricultural systems in the Polish part of the BF based on the acquired spatial data and results of field verification. In some places, linear ridges lie under mounds dated to the early Middle Ages, which allows determining a relative chronology of the fields-the period before the early Middle Ages (Figure 11). 662 calibrated date AD to 770 calibrated date AD/(68.2% probability) 672 calibrated date AD to 765 calibrated date AD; (2) a second archaeological object (prehistoric hearth located in the linear ridge) (95.4% probability), a 138 calibrated date AD to 339 calibrated date AD/(68.2% probability) 223 calibrated date AD to 329 calibrated date AD [70]; and (3) from a place near the field systems (95.4% probability), a 43 calibrated date BC to 125 calibrated date AD/(64.6% probability) 1 calibrated date BC to 82 calibrated date AD [71]. In contrast with the literature on the subject, remains of former fields can be initially associated mainly with the period of Roman influence. The settlement context, recognizing archaeological sites and other factors, also indicate a chronology within the pre-Roman period from the 5th century BC to about the 1st century AD (e.g., Figure 4D). Recognized field systems are dated mainly from about the 1st/2nd century to the 4th century AD. These field systems can be identified as associated with the activities of the Wielbark culture population, while the objects can be associated with the chronologically younger Slavic communities from the period of the 6th-13th century AD. In relation to some of these linear ridges and polygonal spaces, one should also consider the remains of modern constructions. A detailed presentation of the results of field and laboratory tests is the subject of a separate study (in progress). The results of the work are also the basis for further work on the reconstruction of former agricultural systems in the Polish part of the BF based on the acquired spatial data and results of field verification. In some places, linear ridges lie under mounds dated to the early Middle Ages, which allows determining a relative chronology of the fields-the period before the early Middle Ages (Figure 11). Not all linear elevations in the ALS data have been recognized as potential historical objects if their deterioration were too great to recognize them in the field. Furthermore, it is possible that some land elevations may not be related to the post-agricultural clusters. Up to this point, the study of a larger area than the BF territory has not been possible.

Discussion
The application of remote sensing technology for discovering ancient human settlements has developed into an expanding field of research. Our study presents a comprehensive analysis involving the use of ALS in a UNESCO site. We identified a large body of unknown archaeological objects and discovered numerous traces of past human activities in the Polish part of the BF. Many objects were verified in situ ( Figure 12). Not all linear elevations in the ALS data have been recognized as potential historical objects if their deterioration were too great to recognize them in the field. Furthermore, it is possible that some land elevations may not be related to the post-agricultural clusters. Up to this point, the study of a larger area than the BF territory has not been possible.

Discussion
The application of remote sensing technology for discovering ancient human settlements has developed into an expanding field of research. Our study presents a comprehensive analysis involving the use of ALS in a UNESCO site. We identified a large body of unknown archaeological objects and discovered numerous traces of past human activities in the Polish part of the BF. Many objects were verified in situ ( Figure 12).
The application of remote sensing technology for discovering ancient human settlements has developed into an expanding field of research. Our study presents a comprehensive analysis involving the use of ALS in a UNESCO site. We identified a large body of unknown archaeological objects and discovered numerous traces of past human activities in the Polish part of the BF. Many objects were verified in situ ( Figure 12).  Our findings represent new evidence of extended human activity in the history of the BF in addition to that published in 2016 [65]. The most interesting of the identified objects are remains of former settlement and field systems (Figure 11), forming four-sided or polygonal separated areas that are primarily visible in the form of linear deformations or ridges and arranged in compact assemblages. The archaeological objects that could be identified were based on their own landscape sculpture (linear or point) with different widths of elevation, ridges or faults. We assume (based on analogies from the literature [51,55,56,72]) that they represent the remains of ancient field systems (especially as an irregular accreted field or coaxial fields). The linear deformations and ridges are probably the remains of field boundaries, banks and lynchets demarcating the fields or fences, field roads, clearing cairns, ridges or plough headlands. The shape of these structures, known from numerous other studies [52,53,[68][69][70][71][72][73][74], indicates their anthropogenic origin. Objects of this type-remnants of former agricultural systems-are referred to as field systems with distinct linear elevations and specific, closed spaces in the form of multi-dimensional structures varying in size from a few-dozen to over a 100 m long. These objects can also be associated with the results of palynological research, especially in the BF area, which indicates agricultural activity during the period of Roman influence (the 1st century AD to about the 4th century AD) [8].
Our research confirmed in many places the unnatural anthropogenic nature of the stratigraphic units, which is evidenced by, for example, a clear elevation in the central part (visible in the cross-section of objects) and a greater thickness of specific layers (stratigraphic units) in relation to the edge of the elevation. Some objects were also characterized by erratic stone placements found within and in many cases arranged intentionally in the central part of the elevation. Such structure systems show similarity to other objects known from the literature on the subject, which are remnants of former field systems in Europe [75][76][77]. Additionally, artefacts (fragments of ceramic vessels and metal items) were acquired in several excavations, chronologically bound (based on artefact dating) to the Roman period (the 1st century AD to about the 4th century AD) and younger periods. The thickness of the layers occurring directly under the contemporary layer, or even under the bedding, at the elevation of the ridge, is larger than at the edges, and gradually decreases towards the base of the ridge [66].
An important issue regarding this research in forest areas seems to be the quality of the ALS data [30,35], especially the use of publicly available resources, which have become available in many countries, especially in Europe. We need to take into account that some of these data have sparse density, which does not allow for a detailed investigation. Our data were much better than the previously existing resources from national surveys ( Figure 13). The used data were more than two times denser and covered the whole Polish part of the BF area.
artefact dating) to the Roman period (the 1st century AD to about the 4th century AD) and younger periods. The thickness of the layers occurring directly under the contemporary layer, or even under the bedding, at the elevation of the ridge, is larger than at the edges, and gradually decreases towards the base of the ridge [66].
An important issue regarding this research in forest areas seems to be the quality of the ALS data [30,35], especially the use of publicly available resources, which have become available in many countries, especially in Europe. We need to take into account that some of these data have sparse density, which does not allow for a detailed investigation. Our data were much better than the previously existing resources from national surveys ( Figure 13). The used data were more than two times denser and covered the whole Polish part of the BF area.  In a large area analysis, including regions with limited access, such as the Białowieża Forest and other UNESCO sites, it is essential to use non-invasive technologies that allow remote analysis and monitoring. Forest areas require very specific, individual approaches to make the best use of available technology and data, including flights (in areas with low point cloud coverage, dense vegetation or seasonal water bodies) with the use of Unmanned Aerial Vehicle Laser Scanning (UAV-LS) or Mobile Laser Scanning (MLS), in order to fill data gaps or improve data quality [78,79].
We used visual interpretation, which is accurate but time-consuming and requires experienced human resources [45]. Automated methods could provide additional cognitive value. Methods such as machine learning or neural networks are already in use [80][81][82][83][84] and will be explored in the future. Since we have already verified the results in the field, we may have reference materials for automated solutions evaluation. There is still a significant group of potential historic features that have not been either a subject to ALS data interpretation or field verification, nor to a more complete inventory. Many of those do not have strong topographic signatures or are under the surface, which could be detected with ALS, and need more dense point clouds (i.e., UAV-LS) or additional technologies (e.g., geophysics methods) for their detection and recognition.
Another important issue is the primary detection and analysis of resources in sensitive areas where there are real threats to archaeological monuments, whether as a result of natural processes or human activities. Especially in case of large-scale human activities, it is reasonable to initially carry out, for example, ALS/UAV-LS data acquisition and interpretation in such places and to create and share spatial data with possible/verified archaeological heritage detected. Further, priority should be given to areas where natural processes destroying cultural heritage are taking place, such as landslides, windbreaks and emerging excavations; some are found in the Białowieża Forest. ALS has been used recently in many countries to uncover human activity and settlements in the past [85][86][87][88][89]. This is a testimony of often unknown actions in areas that are often considered or were so-called settlements in the past. Our research has provided a multitude of evidence confirming that the Białowieża Forest was modified by human activities during ancient times in ways not known before. Bush et al. [90,91] revealed a 6900-year history of landscape modification by humans in lowland Amazonia, where humans had changed the forest structure and composition as a result of extensive agriculture. Fisher et al. [92] identified ancient settlement patterns using LiDAR in the Mosquitia Region of Honduras and other teams identified ancient settlement patterns in northern Guatemala [93], Central [94] and North America [95].
The evidence gathered in this study contributes significantly to global research of archaeological sites in forested regions. It seems that there is a large gap in our knowledge of the history of many parts of the world that remain under forest cover. The evidence provided in our study contributes to a new understanding of human activity in Białowieża [96] and provides the basis for more detailed analysis in the future.

Conclusions
The extent of the area and the impact of human activities in the Polish part of the BF was substantially greater than previously believed. Substantial evidence was found of human settlement and agriculture during the Roman period (1st century AD to approximately 4th century AD) and/or sometime in the Early Medieval Ages (6th-13th century AD). The remains that have been preserved probably constitute only a part of the multitude of objects that once existed in this area, and potentially also within the Belarusian part of the BF. Some objects have been destroyed, and some were not detected due to the limited quality of the ALS data, which can be related to the specific natural conditions prevailing in some areas of the BF.
An additional product of this study is the established database, the Cultural Heritage-GIS Puszcza Białowieska (CH-GIS Białowieża Forest), as well as a derivative of the abovementioned Web-based geoportal "Natural and cultural portal Puszcza Białowieska"/"Natural and cultural heritage portal of the State Forests" (https://puszcza-bialowieska.lasy.gov.pl/)-allowing for the sharing of resources and promoting the dissemination of research results. This work is part of a broad interdisciplinary effort in Poland, which includes other forests and protected areas, including the Magurski, Roztoczański and Bieszczady National Parks [78][79][80]. The inventory and monitoring of cultural resources in forest areas represents an emerging research area of global dimensions.
We could verify more than 1300 linear objects in the field, thus opening a new chapter in the study of archaeological sites in the Białowieża Forest. The activities confirm the need to undertake further research that takes into account the current state of knowledge about the cultural heritage of the Bialowieza Forest. The goal of subsequent work should be to identify the resources in detail, learn about the history of the area, and effectively manage the UNESCO World Heritage site-for integrated management of the cultural and environmental goods and services.
Author Contributions: K.S. and R.Z. designed the study concept and methodology; K.S. and R.Z. coordinated the research activities leading to this publication; K.S. and R.Z. were responsible for getting the financial support and for the project leading; B.K. and J.W.