Analyses of the Pr ˛adnik riverbed Shape Based on Archival and Contemporary Data Sets—Old Maps, LiDAR, DTMs, Orthophotomaps and Cross-Sectional Proﬁle Measurements

: Analyses of riverbed shape evolution are crucial for environmental protection and local water management. For narrow rivers located in forested, mountain areas, it is di ﬃ cult to use remote sensing data used for large river regions. We performed a study of the Pr ˛adnik River, located in the Ojc ó w National Park (ONP), Poland. A multitemporal analysis of various data sets was performed. Light detection and ranging (LiDAR)-based data and orthophotomaps were compared with classical survey methods, and 78 cross-sectional proﬁles were done via GNSS and tachymetry. In order to add an extra time step, the old maps of this region were gathered, and their content was compared with contemporary data. The analysis of remote sensing data suggests that they do not provide su ﬃ cient information on the state and changes of riverbanks, river course or river depth. LiDAR data sets do not show river bottoms, and, due to plant life, do not document riverbanks. The orthophotomaps, due to tree coverage and shades, cannot be used for tracking the whole river course. The quality of old maps allows only for general shape analysis over time. This paper shows that traditional survey methods provide su ﬃ cient accuracy for such analysis, and the resulted cross-sectional proﬁles can and should be used to validate other, remote sensing, data sets. We diagnosed problems with the inventory and monitoring of such objects and proposed methods to reﬁne the data acquisition.


Introduction
River morphology studies, including riverbed analyses, are carried out by researchers in many different fields, as they are crucial for assessing flood risks, maintaining ecological balance and protecting the environment [1,2]. These studies include riverbed shape, which will be described in this paper, velocity of water, transportation and accumulation analysis etc. [3,4]. The riverbed shape may be monitored with the use of classical surveying methods or various remote sensing methods [5]: 3D laser scanning [6] and synthetic-aperture radar (SAR) [7], Light Detection and Ranging (LiDAR) [8]  The Prądnik is a left-bank tributary of the Vistula River, which is one of the major watercourses in Poland (Figure 1e). The Prądnik is 33.4 km long and its meandering course extends over the northsouth distance of 12.2 km in the Ojców National Park. The analysed part of the Prądnik River is approximately 6 km long. The river may be characterised by slight water level and flow variations [39,40]. The Prądnik River was examined in terms of the physicochemical parameters of water springs [41], the contamination of water and sediments with heavy metals [42], the water conditions in the basins [40,43]. The structure of runoff of the Prądnik River in Ojców is dominated by groundwater runoff (90%) [40]. The water levels of the Prądnik River are remarkably stable, with extremely small amplitudes for medium and extreme states. This is a very unique observation for Polish watercourses [44]. The valley was shaped by karst phenomena. Its floor is overlain by young alluvial deposits consisting of sand, gravel, clay, mud, peat and limestone. These sediments have been accumulating in this region for 10 thousand years [36,45]. In the Holocene, sediments of the Prądnik River valley distinguished three types of stratigraphic sequences [46]. In the study area, there are traces of deposits of narrow valleys crossing the elevated part of the upland. Molluscan fauna, in particular types of sediments, reflects the evolution of the environment and the transformation of ecosystems in the Prądnik River valley. The riverbed shape of the Prądnik River was influenced by both natural and anthropogenic processes. Along with the development of industry, mills, fishing sites and roads were created in Ojców and the Prądnik River valley [39]. The Prądnik River valley is affected only sporadically by floods, but they are high in volume and range. During the flood on May 18th 1996, the measured volume of the Prądnik River reached 37 m 3 /s for one km 2 of catchment. This is a record value registered for this part of Europe [47]. During that time, a flood wave of 3-4 m was observed, which left the whole of the valley covered in 0.5-0.8 m of water. Such large floods cause The monitoring of natural resources of the Ojców National Park is a subject of various studies [48][49][50][51]. They concern such environmental features as meteorology, air quality, animal habitats, number of species, plants etc. Such activities aim at providing long-term observations of existing patterns, assessing implemented protection measures and preserving natural character of the ONP [52]. Hence, the research concerning detecting changes in the riverbed of the Prądnik over time meets real needs of the Ojców National Park, and may become a valuable component of the aforementioned environmental studies. The shape of the Prądnik River riverbed and its changes have not been studied so far.

Materials and Methods-Overview of Available Data, Data Preparation and Data Acquisition
The following flowchart ( Figure 2) represents used data sets and methods of processing. In some cases, it shows relations between the data sets and how they were used in the Prądnik riverbed shape analysis.

Materials and Methods-Overview of Available Data, Data Preparation and Data Acquisition
The following flowchart ( Figure 2) represents used data sets and methods of processing. In some cases, it shows relations between the data sets and how they were used in the Prądnik riverbed shape analysis.

LiDAR Based Point Clouds and DTMs of the Prądnik River Valley
During the course of this study, the authors managed to obtain two current airborne LiDAR based point clouds. These products were acquired as .las files that included XYZ coordinates, intensity and the RGB colour of each measured point.
The first point cloud was created during the ISOK program. This program aims to gather the data on the terrain of Poland that could be used for early warning and monitoring of natural disasters, mainly floods, to which the area is susceptible and that in the past has caused a lot of damage. The ISOK data set consists of 3 types of products: a point cloud, a DTM and a digital surface model (DSM) in raster format.
The point cloud was acquired with the use of an airborne platform consisting of an LMS-Q680i RIEGL laser scanner, GPS antenna and a camera installed on a plane flying at about 950 m height. The point clouds were registered with the use of flat surfaces of the rooftops. Four corners measured via the means of tachymetry defined the surfaces. The tachymetric measurements were conducted

LiDAR Based Point Clouds and DTMs of the Prądnik River Valley
During the course of this study, the authors managed to obtain two current airborne LiDAR based point clouds. These products were acquired as .las files that included XYZ coordinates, intensity and the RGB colour of each measured point.
The first point cloud was created during the ISOK program. This program aims to gather the data on the terrain of Poland that could be used for early warning and monitoring of natural disasters, mainly floods, to which the area is susceptible and that in the past has caused a lot of damage. The ISOK data set consists of 3 types of products: a point cloud, a DTM and a digital surface model (DSM) in raster format.
The point cloud was acquired with the use of an airborne platform consisting of an LMS-Q680i RIEGL laser scanner, GPS antenna and a camera installed on a plane flying at about 950 m height. The point clouds were registered with the use of flat surfaces of the rooftops. Four corners measured via the means of tachymetry defined the surfaces. The tachymetric measurements were conducted from two points measured by real-time kinematic GPS (RTK-GPS), one was used as a setup point for tachymetry and the other as a reference. The RGB colours were added to the registered point cloud from an orthomosaic of pictures taken during separate flights. The scanning parameters were as follows: the size of the laser beam on the surface was −0.29 m; the density of the scan was −4.7 pts/m 3 ; the density of points varied, being denser in the direction parallel to the flight route at −0.45 m, and less dense in the perpendicular direction of −0.47 m; the overlap between parallel scanning lines was −35%. The observed density varies within a typical inter-point distance of 0.20 m to 0.45 m. This point cloud Remote Sens. 2020, 12, 2208 6 of 30 was acquired in spring and summer 2013. The accuracy of the registration procedure was reported as 0.25 m.
The second data set was obtained from the Ojców National Park, unfortunately limited metadata were available. The authors of this study managed to find the original order on which a public auction was based. The authorities of the National Park confirm that, even though the acceptance report was lost, the point clouds did fulfil the minimal requirements of the order. The order states that the scanning parameters should be as follows: density of the scan −20 pts/m 3 , overlap between parallel scanning lines-50%. The survey should be carried out before the 15th of June 2012, the registration should be based on 6 flat surfaces and its accuracy should be of at least 0.20 m. RGB colour was added to the registered point cloud based on the orthomosaic acquired during the same flight. This data set was done for the general inventory of the ONP, most likely in late spring, evidenced by the state of vegetation. The observed density varies within a typical inter-point distance of 0.10 m to 0.35 m.
In addition to the point cloud, the ISOK project provided us with DTM and DSM models developed at a resolution of 1 m. The height accuracy of the DTM model is +/−0.15 m in open spaces and +/−0.3 m in forest areas. On the slopes these accuracy values are even smaller [53]. The height accuracy of the DSM model reaches +/−1 m. Due to the nature of the terrain (forested, steep limestone rocks), the DTM is also characterised by low vertical accuracy (explained later by Figure 8(S 52)) [53]. Based on the point cloud from the Ojców National Park, the DTM and the DSM were also produced with a 0.5 m resolution. According to the order specification, the height accuracy of the DTM and the DSM should be +/−0.3 m (no differentiation per type of surface).
LiDAR based data use Polish cartographic projection PL-1992 (EPSG code: 2180). All other data described in this manuscript were transformed to this projection. The PL-1992 is based on a geocentric ellipsoid GRS-80. For the whole territory of Poland, one zone of the Gauss-Kruger projection is used. The central meridian is 19 • with a scale factor of 0.9993 [54].
The characteristics of the data used in the article are described in Table S1 (Table S1-All used data and its characteristic.pdf) in Supplementary Materials.

Survey of the Prądnik River and its Floodplain
The Prądnik can be described as a narrow, alluvial river with a highly visible floodplain that is surrounded by significantly higher limestone outcrops and stacks, characteristic to this area, partly covered with woods. The valley is 40 m wide in its narrowest part and 180 meters in its widest, also partly covered with trees and buildings. The use of GNSS systems is significantly limited, due to these factors [55,56].
The initial plan for measuring cross-sectional profiles of the river and its floodplain involved the usage of real-time network (RTN) mode with ASG-EUPOS (the network of ground reference stations in Poland) corrections or real-time kinematic (RTK) with corrections from a base located in the nearest area with a good sky visibility. This remote sensing method is the most time-effective due to the fact that it does not require measuring a traverse and can be done with a limited number of personnel. However, its implementation was not possible, due to the field conditions-inadequate satellite coverage above the horizon, which was limited by trees and hills, or a lack of the correction transmission from the system or base station [57]. This proves that current remote sensing methods are severely limited in the areas of narrow mountain rivers. Hence, it was decided to place a traverse along the river and conduct a tachymetric survey from its points, integrated with the GNSS points.
It was not possible to find data on enough survey marks from the Polish national survey mark system, thus 2 pairs (Table 1) of traverse points were measured via static GNSS session [58,59]. The measurements were conducted with 2 Leica GS16 GNSS receivers, collecting the observations synchronically. All of the receivers were mounted on geodetic tripods, the heights of which were reduced to the level of measured points. The elevation angle was set to 10 • for all points, each session took two hours with 1 s sampling interval. The points were measured in pairs-a pair at the south end of the traverse was measured first, then a pair in its northern part was surveyed. The GNSS static Remote Sens. 2020, 12, 2208 7 of 30 survey was conducted on the 10th of April 2019. The GNSS measured traverse points were set at the ends of the valley, the only two areas where the satellite coverage was sufficient for a static survey via this remote sensing method. The post-processing of the measurements was carried out using Leica Geo Office 8.4 software. As a reference, observations form 4 ASG-EUPOS stations ( Figure 3) were chosen (KATO-Katowice, KRA1-Kraków, PROS-Proszowice, LELO-Lelów), located 17-57 km from the analysed area.
Remote Sens. 2020, 12, x FOR PEER REVIEW 7 of 30 measurements were conducted with 2 Leica GS16 GNSS receivers, collecting the observations synchronically. All of the receivers were mounted on geodetic tripods, the heights of which were reduced to the level of measured points. The elevation angle was set to 10° for all points, each session took two hours with 1 s sampling interval. The points were measured in pairs-a pair at the south end of the traverse was measured first, then a pair in its northern part was surveyed. The GNSS static survey was conducted on the 10 th of April 2019. The GNSS measured traverse points were set at the ends of the valley, the only two areas where the satellite coverage was sufficient for a static survey via this remote sensing method. The post-processing of the measurements was carried out using Leica Geo Office 8.  The post-processing strategy was based on the precise final orbits from the International GNSS Service (IGS) [60] and GPS+GLONASS observations, Hopfield tropospheric model [61], automatic ionosphere model, stochastic modelling and fix ambiguities solutions. The heights (H) in the normal height system were determined with the use of the current quasi-geoid model for the Polish territory-PL-geoid-2011 [54]. The coordinates of 4 reference points and their accuracy are given in Figure 3.
The cross-sectional profiles were measured during two days of a survey session on the 8 th and 9th of April 2019. The survey included the measurements of distances and bearings in two series and two positions of the telescope (transits of the scope). Topcon OS total station with a paired standard prism was used in this survey. The accuracy of this instrument is 3" (ISO 17123-3:2001) for a single angle measurement, ±(2 mm + 2 ppm × D) for prism measurements and ±(3 mm + 2 ppm × D) for no prism (ISO 14123-4:2001). The range of this tachymeter is 500 m, thus, it was enough for the purpose of this study. There were 31 traverse points, placed in the average distance of 110 m. Two points at The post-processing strategy was based on the precise final orbits from the International GNSS Service (IGS) [60] and GPS+GLONASS observations, Hopfield tropospheric model [61], automatic ionosphere model, stochastic modelling and fix ambiguities solutions. The heights (H) in the normal height system were determined with the use of the current quasi-geoid model for the Polish territory-PL-geoid-2011 [54]. The coordinates of 4 reference points and their accuracy are given in Figure 3.
The cross-sectional profiles were measured during two days of a survey session on the 8th and 9th of April 2019. The survey included the measurements of distances and bearings in two series and two positions of the telescope (transits of the scope). Topcon OS total station with a paired standard prism was used in this survey. The accuracy of this instrument is 3" (ISO 17123-3:2001) for a single angle measurement, ±(2 mm + 2 ppm × D) for prism measurements and ±(3 mm + 2 ppm × D) for no prism (ISO 14123-4:2001). The range of this tachymeter is 500 m, thus, it was enough for the purpose of this study. There were 31 traverse points, placed in the average distance of 110 m. Two points at each end of the traverse were also measured with a static GNSS survey, and were used as a reference to the coordinate system. The accuracy of traverse points can be described by horizontal and vertical root mean square error (RMSE) values, which were 43 mm and 13 mm, respectively. The survey of the cross-section profiles was conducted simultaneously with the traverse measurements. Each measured Remote Sens. 2020, 12, 2208 8 of 30 point was surveyed twice with a standard tachymetric method. The average error of each point in XY was 0.043 m, the largest error did not exceed 0.057 m, and the average error in H was 0.013 m and did not exceed 0.016 m.
As a result, 78 cross-section profiles were done, each of 9 points on average, a total of 718 points were measured. The profiles consisted of 1 point on each side of the river representing the water level, at least 1 point in the deepest part of the river and 2-5 points representing the surrounding terrain on each side of the river.

Orthophotomaps
Orthophotomaps, being another kind of remote sensing data created to provide accurate spatial information on vast areas, were obtained from National Geodetic and Cartographic Resource. A vector layer of the river course was created manually for each orthophotomap, using the ArcGIS ArcMap software. The polygons were constructed in the areas where the river was visible. Some areas were covered by tree crowns or shadows, so it was impossible to identify the water course there. The data set from 2019 provided the least data, because the images were taken during spring, when the vegetation significantly limited visibility. It was also difficult to analyse the data set from 2009, as a photogrammetric flight took place during a sunny day, which caused shadows and inconsistencies in contrast. The most data were derived from the orthophotomap from 2003. The approximate length of the vectorised river is presented in Table 2.
During the research, two additional old maps were considered unsuitable for further analyses of the Prądnik River. These maps are the Map of the Free City of Kraków (pol. Mapa Okręgu Wolnego Miasta Krakowa) from 1833 and the Map of caves in the vicinity of Kraków and Ojców (pol. Mapa jaskiń okolic Krakowa i Ojcowa) from 1929. As a result of the lack of any characteristic objects, it was impossible to georeference them and compare to the current data. However, an approximate location of the riverbed is shown on these maps.
One of the ways to analyse old maps is to compare them by overlaying on the present ones. Problems which can occur while comparing maps from different times are: • Georeferencing-through ages various maps were created using many cartographic projections with different parameters and accuracy [67]. Some maps were presented with established mathematic assumptions, and some were created without any projection. Georeferencing of this type of maps is possible only by fitting based on points, which have not changed since the maps were created. Sometimes, it is difficult to define such points. Furthermore, deviations on non-cartometric maps are not systematic, thus, it is necessary to use many evenly distributed points in order to perform any georeferencing.

•
Interpretation-a map is defined as an ordered and generalised model of reality represented in an arbitrary way using cartographic signs system [68]. Over time, the definitions of symbols have changed. The first complete legend on Polish maps was presented on topographic maps published by the Military Geographic Institute in Warsaw in the 1930s. The first printed maps were created only in one or two colours, so different signs were used to distinguish objects, often very similar to each other and in fact undistinguishable, even within one map sheet. Additionally, a problematic factor is that Polish maps published during the 123 years when Poland was not a sovereign country (from 1795 until the end of the First World War) were created in different countries. Printing techniques, scales, units, sheet formats and used language were different for particular regions occupied by the countries participating in the Partitions of Poland.

•
Scale-different scales result in different importance of errors [67]. For example, georeferencing a map at 1:10,000 scale with an RMS error of 100 m would produce an error of 10 mm on this map and can be significant. However, on a map at 1:100,000 scale, the same terrain error results in 1 mm error. This would probably be neglectable as, most likely, the symbols representing objects are larger than the error value.
The first step was the transformation of map scans to one projection based on characteristic points whose location was assumed unchanged through time. The lack of old buildings in the national park induces the decision of choosing less clearly defined objects, such as crossroads, for this task. The reference data was a layer of roads from the Database of Topographic Objects (BDOT10k). The maps were transformed to this projection using the ArcGIS ArcMap software. A transformation method was set as the first order polynomial (affine), which is based on the least-square fitting algorithm. The minimal number of control points for this transformation is 3. Better results are obtained by using more than 3 points distributed evenly around the area of interest. As a result, a raster dataset is scaled, shifted and rotated and receives a new georeference. The next step was the manual vectorisation of the Prądnik River on each map.

The Josephine map of Galicia
The oldest map, the Josephine map of Galicia ( Figure 4a,b), was created over the 1779-1783 period for military purposes, and includes lands occupied by Austria after the First Partition of Poland. This map is also called the Mieg map of Galicia. It is a good cartographic source for a historical study of urban development [69]. The map of Galicia was created at 1:28,800 scale. The accuracy of these type of maps is often described as 'irregular', as some objects were placed only roughly, and the quality depends on the cartographer's skills [70]. The survey was based neither on precisely defined triangular points, nor on clearly defined projection. Objects were mapped with an unaided eye and the attention was paid to details which were significant from the military planning perspective, such as communication and transport, castles, churches, rivers [71]. The accuracy of georeferencing of this map can be around 100-200 meters, even less if just a small part of the sheet is georeferenced [72]. The reference points were derived from the Detailed Map of Poland at 1:25,000 scale, which was georeferenced earlier. This process simplified the identification of control points. A total of 8 points were chosen, mostly on characteristic crossroads. The total RMSE is 87.77 m. The river course was vectorised manually.

The Map of Ojców surroundings
The map of Ojców surroundings, Edition of the Ojców Bazaar map (Figure 4c,d) was published in 1907 at 1:100,000 scale. The georeferencing of this map was based on 6 control points located around the river. The reference data was the roads layer of the Topographic Objects Database (BDOT10k). The total RMSE is 27.44 m. The river course was vectorised manually.

The Russian and Soviet "two-verst"
The Russian and Soviet map from 1914 is a "two-verst" map (Figure 4e,f) so it was compiled at 1:84,000 scale. It was the main topographic map for civil and military purposes until 1930s. It was created not only for the Russian territory, but also for Poland, Latvia, Estonia, the Crimea and the Western Caucasus [70]. The map was issued in a two-colour version-black for terrain objects and brown for contour lines. The georeferencing of this map was based on 6 control points located around the analysed part of the river. The reference data was a layer of roads from the Topographic Objects Database (BDOT10k). Total RMSE is 11.05 m. The river course was vectorised manually.

WIG25k
The Detailed Map of Poland at 1:25,000 scale was published by the Military Geographical Institute (WIG) in Warsaw (Figure 4g,h). This map series covered around 50% of the current Polish area and was developed in stages, which causes differences in graphical representation or the reference ellipsoid. The sheet which was used for this study was published in 1935 as a "for tourism" version. It was printed in five colours. The source maps were mainly large-scale maps from the period of the Partitions of Poland [73]. The first step of georeferencing of this map was transforming it into its original cartographic projection. Then, the vector layer that represents the Prądnik was created manually. Later, just the vector layer was spatially adjusted to the PL-1992 projection using control points located on bridges. An affine transformation was used, the total RMSE is 19.80 m.

The Godfryd Ossowski map
As part of the research concerning using old maps to study the shape of the Prądnik riverbed, a map at 1:10,000 scale was found [74]. Unfortunately, the map coverage did not include the whole of the study area, only its south part ( Figure 5). This map was created to mark the location of caves in the Prądnik valley. Mills with millraces are also easy to recognise on this map ( Figure 5). The

The Map of Ojców surroundings
The map of Ojców surroundings, Edition of the Ojców Bazaar map (Figure 4c,d) was published in 1907 at 1:100,000 scale. The georeferencing of this map was based on 6 control points located around the river. The reference data was the roads layer of the Topographic Objects Database (BDOT10k). The total RMSE is 27.44 m. The river course was vectorised manually.
The Russian and Soviet "two-verst" The Russian and Soviet map from 1914 is a "two-verst" map (Figure 4e,f) so it was compiled at 1:84,000 scale. It was the main topographic map for civil and military purposes until 1930s. It was created not only for the Russian territory, but also for Poland, Latvia, Estonia, the Crimea and the Western Caucasus [70]. The map was issued in a two-colour version-black for terrain objects and brown for contour lines. The georeferencing of this map was based on 6 control points located around the analysed part of the river. The reference data was a layer of roads from the Topographic Objects Database (BDOT10k). Total RMSE is 11.05 m. The river course was vectorised manually.

WIG25k
The Detailed Map of Poland at 1:25,000 scale was published by the Military Geographical Institute (WIG) in Warsaw (Figure 4g,h). This map series covered around 50% of the current Polish area and was developed in stages, which causes differences in graphical representation or the reference ellipsoid. The sheet which was used for this study was published in 1935 as a "for tourism" version. It was printed in five colours. The source maps were mainly large-scale maps from the period of the Partitions of Poland [73]. The first step of georeferencing of this map was transforming it into its original cartographic projection. Then, the vector layer that represents the Prądnik was created manually. Later, just the vector layer was spatially adjusted to the PL-1992 projection using control points located on bridges. An affine transformation was used, the total RMSE is 19.80 m.

The Godfryd Ossowski map
As part of the research concerning using old maps to study the shape of the Prądnik riverbed, a map at 1:10,000 scale was found [74]. Unfortunately, the map coverage did not include the whole of the study area, only its south part ( Figure 5). This map was created to mark the location of caves in the Prądnik valley. Mills with millraces are also easy to recognise on this map ( Figure 5). The publication describing the caves and the map were made by Godfryd Ossowski, one of the fathers of Polish archaeology and the editor of numerous maps [75].
Remote Sens. 2020, 12, x FOR PEER REVIEW 11 of 30 publication describing the caves and the map were made by Godfryd Ossowski, one of the fathers of Polish archaeology and the editor of numerous maps [75]. The map was georeferenced using the BDOT10k database and a 1:10,000 topographic map. The affine method was used. RMSE is 8.40 m, which corresponds to an error of 0.84 mm on the map at 1:10,000 scale. This result should be considered satisfactory, given the age of this map. The map was probably made using plane table-based survey.

Contemporary Vector Data
Due to the small width of the river in the analysed area, only polyline layers can be found in various spatial databases. The polyline layers used during this study are: • Map of the Hydrographic Division of Poland (MPHP) database. This database was created using the topographic 1:50,000 scale maps and was up to date as of 2010 [76]. Objects are stored in the MPHP database, it is a typical representation of a river network on medium-scale maps [77]. Contrary to other analysed databases, generalised (quantitative and qualitative) objects are stored in the MPHP database.

•
The second vector layer was the ONP watercourses (ONPw). This layer was created by the employees of the ONP using manual vectorisation on a high resolution DTM, based on the point cloud from 2012 from the ONP, described in Section 3.1. We used it as a base layer and other layers were compared to it. It is the most current vector layer, and is derived from the most detailed data set.

•
The third layer was one of the feature classes of the BDOT10k database-SWRS (pol. sieć wodna -rzeki, strumienie, ang. water network -rivers, streams). The SWRS layer stores river and stream axes. The database was based on the 17.11.2011 Regulation on 'the topographic objects database and the database of general geographic objects, standards of cartographic studies' [78]. The BDOT10k database is currently the largest georeferenced Polish database [79]. It was created over many years using various techniques. The analysed fragment of the database is up to date as of 2009 (geometry) and 2013 (attributes). Data in the BDOT10k are in Polish cartographic projection PL-1992 (EPSG code: 2180).

•
The next layer (Topo10k92) is not derived from any database, but from the vectorisation of topographic maps at the 1:10,000 scale. These are digital maps but the source layers that were used to render them are not available. Therefore, it was decided to vectorise GeoTIFF files. The vectorisation was based on two sheets: south M-34-64-B-c-4 "Biały Kościół" and north M-34-64-B-c-2 "Skała". The "Skała" sheet is up to date as of 1996 and "Biały Kościół" sheet as of 2002. It should be noted that these two sheets differ, not only in topicality, but also in technical The map was georeferenced using the BDOT10k database and a 1:10,000 topographic map. The affine method was used. RMSE is 8.40 m, which corresponds to an error of 0.84 mm on the map at 1:10,000 scale. This result should be considered satisfactory, given the age of this map. The map was probably made using plane table-based survey.

Contemporary Vector Data
Due to the small width of the river in the analysed area, only polyline layers can be found in various spatial databases. The polyline layers used during this study are: • Map of the Hydrographic Division of Poland (MPHP) database. This database was created using the topographic 1:50,000 scale maps and was up to date as of 2010 [76]. Objects are stored in the MPHP database, it is a typical representation of a river network on medium-scale maps [77]. Contrary to other analysed databases, generalised (quantitative and qualitative) objects are stored in the MPHP database.

•
The second vector layer was the ONP watercourses (ONPw). This layer was created by the employees of the ONP using manual vectorisation on a high resolution DTM, based on the point cloud from 2012 from the ONP, described in Section 3.1. We used it as a base layer and other layers were compared to it. It is the most current vector layer, and is derived from the most detailed data set.

•
The third layer was one of the feature classes of the BDOT10k database-SWRS (pol. sieć wodna -rzeki, strumienie, ang. water network -rivers, streams). The SWRS layer stores river and stream axes. The database was based on the 17 November 2011 Regulation on 'the topographic objects database and the database of general geographic objects, standards of cartographic studies' [78]. The BDOT10k database is currently the largest georeferenced Polish database [79]. It was created over many years using various techniques. The analysed fragment of the database is up to date as of 2009 (geometry) and 2013 (attributes). Data in the BDOT10k are in Polish cartographic projection PL-1992 (EPSG code: 2180).

•
The next layer (Topo10k92) is not derived from any database, but from the vectorisation of topographic maps at the 1:10,000 scale. These are digital maps but the source layers that were used to render them are not available. Therefore, it was decided to vectorise GeoTIFF files. The vectorisation was based on two sheets: south M-34-64-B-c-4 "Biały Kościół" and north M-34-64-B-c-2 "Skała". The "Skała" sheet is up to date as of 1996 and "Biały Kościół" sheet as of 2002. It should be noted that these two sheets differ, not only in topicality, but also in technical requirements [80]. The sheet covering the southern part of the study area was issued in accordance with the guidelines of 1999, while the northern sheet was issued in accordance with the guidelines of 1989 [81]. The PL-1992 projection was used during creation of these maps.

•
The last layer (Topo10k65), similarly to the previous one, is based on the vectorisation of topographic maps at the 1:10,000 scale, but from an earlier edition. The maps were made in analogue technology. The vectorisation was performed using the WMS layer made available by the Central Office of Geodesy and Cartography in Poland. The study area was covered on two sheets: 163.311 "Wielka Wieś" up-to-date as of 1983-1986 and 163.133 "Skała", up-to-date as of 1978-1979. These were maps made in accordance with the technical instructions guidelines of 1980 [82]. The maps were made in the PUWG 1965 projection (zone 1) (EPSG: 3120).

Point Cloud Comparison
In order to estimate the quality of our two primary remote sensing data sets, point clouds, a cross evaluation was done. Some initial conclusions are listed below: • None of the data sets could be used for measuring the water floor-topographic LiDAR does not allow for accurate measurements through water, and no other application, for example, bathymetric LiDAR, was used to provide this data, the ISOK data set is incomplete regarding flood monitoring and prediction.

•
The ONP point cloud was created relatively fast, within one or two days, judging by the vegetation. It was acquired most likely during late spring, which is unfavourable, due to a large number of deciduous trees that cover the area, thus limiting the density of point cloud representing terrain and creating a large number of blind spots, without any points. However, its short time of creation and relatively high density works in its favour.

•
The ISOK point cloud was created during at least two survey series, one in early spring, evidenced by snow cover still visible in the gullies, and the other in mid to late spring, evidenced by the trees and the colour of the grass. Further analysis showed that this point cloud was 'patched-up' with other point clouds obtained during many separate surveys, this is evidenced by a significant difference of illumination and difference in the density of the point clouds. Some of those things can be attributed to the fact that images for RGB colour were taken separately, but changes in density and in the vegetation cover are too significant to ignore. As a result, this point cloud is less dense and harder to interpret.

•
The ONP point cloud was created with unified density, only changed locally by the amount and type of vegetation. The ISOK point cloud, most likely due to its patched-up origin, consists of square areas of low density, one point every 0.60 m, and areas of high-density, one point every 0.20 m. Additionally, in the ONP, the points are aligned in one direction while in the case of ISOK scan lines change with density. It is worth noting that the metadata do not reflect this information for ISOK. A user might be under the impression that the LiDAR session was done in one set, while following a precisely designed flight plan. No information on changes during scanning or additional data acquisition was given.
Since the point clouds can be viewed as the most recent 3D representation of the valley, we decided to validate them against one another. This was done as follows. The point clouds were loaded into AutoDesk ReCap software. The same, easy to identify objects were localised on both point clouds; these were rooftops, edges of the bridge, corners of the buildings, church tower etc. They were dispersed evenly around the area. Their coordinates were compared, the difference was at 0.34 m at average. Additionally, the −/+ difference was checked, to see if the point clouds were not tilted against each other. This was not observed, the differences were evenly dispersed along the point clouds.
For further processing, the non-terrain elements, trees, buildings, ruins of the castle, and all found artefacts were removed from the point clouds. They were also cut, so that they represent the same parts of the valley and surrounding rock structures. The distances between the point clouds were measured with the use of the CloudCompare software ( Figure 6). The initial analysis shows further proper alignment of the point cloud in the vicinity of the road, the only stable, man-made object not extracted for this analysis. The analysis of the road shows that the difference in distance is about 0.10 m-this is fair less than the accuracy of registration for the ISOK point clouds. A total of 65% of surrounding terrain shows similar differences. A total of 15% of differences larger than 0.30 m can be attributed to plant life, more specifically, all grass or small bushes that were impossible to remove. This difference did not exceed 0.50 m. The remaining 20% of differences cannot be attributed to that. This is interesting, since the largest, up to 0.60 m, changes appeared in the east bank of the river (Figure 8(S 12)), smaller, 0.40-0.50 m differences are visible throughout the valley on the east bank and are usually larger in less urbanised areas and on parts of land between the river and the road (this can be either west or east). If terrains with those 20% of differences attributed to the grass are excluded, the models agree with each other well.
Remote Sens. 2020, 12, x FOR PEER REVIEW 13 of 30 For further processing, the non-terrain elements, trees, buildings, ruins of the castle, and all found artefacts were removed from the point clouds. They were also cut, so that they represent the same parts of the valley and surrounding rock structures. The distances between the point clouds were measured with the use of the CloudCompare software ( Figure 6). The initial analysis shows further proper alignment of the point cloud in the vicinity of the road, the only stable, man-made object not extracted for this analysis. The analysis of the road shows that the difference in distance is about 0.10 m-this is fair less than the accuracy of registration for the ISOK point clouds. A total of 65% of surrounding terrain shows similar differences. A total of 15% of differences larger than 0.30 m can be attributed to plant life, more specifically, all grass or small bushes that were impossible to remove. This difference did not exceed 0.50 m. The remaining 20% of differences cannot be attributed to that. This is interesting, since the largest, up to 0.60 m, changes appeared in the east bank of the river (Figure 8 S12), smaller, 0.40-0.50 m differences are visible throughout the valley on the east bank and are usually larger in less urbanised areas and on parts of land between the river and the road (this can be either west or east). If terrains with those 20% of differences attributed to the grass are excluded, the models agree with each other well. Figure 6. Difference in height between the point clouds, the Prądnik River is not represented by points, the road is constantly represented in blue, difference in height is minimal. A-minimal difference between point clouds, B-largest difference in height between the river and the road, Clargest difference on the east bank in a non-urbanised area (visible difference between urbanised and non-urbanised area).

Comparison between DTMs and Measured Cross-Section Profiles of the Prądnik River.
The authors compared existing spatial data, the ONP and the ISOK based DTMs with measured cross-section profiles, in order to validate their quality and observe if any significant changes occurred between the times when the ISOK data were obtained (2013) and the cross section profiles were measured (2019). A profile analysis was performed as follows. Surveyed points from one profile were connected by a polyline in the ArcGIS ArcMap software (Figure 7a). On this polyline, vertices in 0.5 m sampling were generated. Then, the heights for these vertices were taken from the ONP and the ISOK DTMs (Figure 7b,d). Three height sets-surveyed points and vertices with heights from the DTM ONP and the DTM ISOK models-were compared in Microsoft Excel (Figure 7f). During the

Comparison between DTMs and Measured Cross-Section Profiles of the Prądnik River
The authors compared existing spatial data, the ONP and the ISOK based DTMs with measured cross-section profiles, in order to validate their quality and observe if any significant changes occurred between the times when the ISOK data were obtained (2013) and the cross section profiles were measured (2019). A profile analysis was performed as follows. Surveyed points from one profile were connected by a polyline in the ArcGIS ArcMap software (Figure 7a). On this polyline, vertices in 0.5 m sampling were generated. Then, the heights for these vertices were taken from the ONP and the ISOK DTMs (Figure 7b,d). Three height sets-surveyed points and vertices with heights from the DTM ONP and the DTM ISOK models-were compared in Microsoft Excel (Figure 7f). During the field measurements, we tried to keep straight lines for each profile, but this was not always possible (e.g., due to the visibility of surveyed points from the station). Analyses in MS Excel were performed along the section lines connecting the surveyed points.  Figure 8 shows the comparison between several typical profiles with DTMs. It is important to mention at this point that the point clouds, on which the DTMs were based, did not consist of points from the river bottoms. This is due to the fact that airborne laser scanners usually do not penetrate water or create many artefacts when a laser beam passes through this medium. Everything given as a reflection from the water surface and realised as points within a point cloud can be removed via various algorithms, if the inspection of intensity defines those points as false measurements (Figure 8(S 11)). The ISOK DTM is often higher than the ONP DTM in the area of the river bottom, thus, the ISOK represents the riverbanks ending at the water surface and the ONP DTM ends near the water bottom ( Figure 8(S 1)). This is most likely due to the different laser scanners used and the different algorithms applied for filtering the noise.
This means that, only differences above a threshold of the water level could be analysed. The DTMs seem to be placed higher in all points measured, this is most likely the result of the technology used during scanning, and also the time of the year when the grass was tall enough to be recorded in a point cloud. This last point is the most significant since the steeper the area gets (less vegetation grows there), the smaller the differences are.
The comparison of models carried out on the profiles shows that in places of large slopes (slopes or rocks), models generated from the point clouds are too smooth. This can result from malfunctioning algorithms for the automatic detection of watercourses. This situation can be observed on the right side of profile 52 (DTM ISOK) and 67 (DTM ONP) and on most profiles on the water-land border ( Figure 8). This can be observed more often on the DTM ONP model, despite its greater resolution.  The analysis of profiles allows for concluding that for small rivers, similar to the Prądnik River, waterbed modelling cannot be done accurately by using airborne laser scanning. This applies to both products; point clouds and DTMs. Terrestrial laser scanning could be used for this purpose, but this would require clearing the area, which is impossible due to the National Park regulations of undisturbed nature being left on display. Using UAV with LiDAR or a camera that would allow for structure-from-motion processing could provide sufficient data. However, this would require extensive planning. Airborne LiDAR point clouds, analysed in this study, did not provide enough data, thus, the methodology would have to be customised strictly for this area of study. Another issue is the Polish law, which requires special permissions from the National Park administration and others to perform a UAV flight. This suggests that terrestrial laser scanning would be a reasonable choice for such a survey.
Further profile analysis indicates that, at times, the DTM ONP model shows a wide and flat river bottom (Figure 8(S 11,S 52)). Creating such "lakes" should be described as incorrect. This is one of the problems that may be encountered when using DTMs which were based on airborne laser scanning (ALS) data and created with general purpose algorithms, for hydrological analyses.
In order to give the final check of the quality of the DTMs, a simple statistical analysis was performed. A total of 478 measured points were used (points measured under the water surface were excluded). This was based on calculating the distances between the ISOK and ONP DTMs and measured cross-section profiles. The results are presented in Figure 9 and Table 3. This analysis has shown that the ISOK DTM agrees better with measured data, since the difference is significantly smaller than in the case of the ONP. In addition, the ONP DTM is less consistent, which is evidenced by the outliers. This all proves that the ISOK data set, even with previously described problems with density, is better suited for this kind of analysis.
Remote Sens. 2020, 12, x FOR PEER REVIEW 17 of 30 In order to give the final check of the quality of the DTMs, a simple statistical analysis was performed. A total of 478 measured points were used (points measured under the water surface were excluded). This was based on calculating the distances between the ISOK and ONP DTMs and measured cross-section profiles. The results are presented in Figure 9 and Table 3. This analysis has shown that the ISOK DTM agrees better with measured data, since the difference is significantly smaller than in the case of the ONP. In addition, the ONP DTM is less consistent, which is evidenced by the outliers. This all proves that the ISOK data set, even with previously described problems with density, is better suited for this kind of analysis.  The current shape of the Prądnik River was evaluated based on the measured cross-section profiles and DTMs. The width and depth of the river is presented in Table 4. The maximum and minimum width are 9.39 m and 1.53 m. The average width is 3.75 m. In most areas of the river, the depth is lower than 0.3 m, and only in two areas the depth is greater than 1 m. The difference in height between the highest and the lowest point of the river bottom is 31.09 m. The riverbanks are steeper on the west bank.   The current shape of the Prądnik River was evaluated based on the measured cross-section profiles and DTMs. The width and depth of the river is presented in Table 4. The maximum and minimum width are 9.39 m and 1.53 m. The average width is 3.75 m. In most areas of the river, the depth is lower than 0.3 m, and only in two areas the depth is greater than 1 m. The difference in height between the highest and the lowest point of the river bottom is 31.09 m. The riverbanks are steeper on the west bank.

Orthophotomaps Analysis
We planned to use the orthophotomaps to perform a detailed analysis of changes of the Prądnik River shape. However, it was impossible to vectorise the whole course of the Prądnik River, because of the high and dense vegetation that covered the river, and shadows which occurred in many parts of the analysed area.
In few places, three or four vector layers were visible ( Figure 10, middle row). Apart from those locations, it was possible to compare river shape and detect potential changes. One of the differences can be seen between the second and the third cross-section profile, where the Prądnik River shifted to the west (Figure 10, bottom row). Comparing the DTM ONP from 2012 and the orthophotomap from 2017, this change occurred after 2012 and before 2017, and the shift reaches 3 m. Other significant changes are visible on cross-sections no. 62 and no. 64 ( Figure 10, top row). The river course on the orthophotomap from 2003 is similar to the BDOT10k vector layer from 2006, but the orthophotomap from 2009 shows a different course of the river in this location, so the changes occurred between 2006 and 2009. The characteristics and placement of these changes, in the middle of a farmland, suggest that they were manmade. There was a significant change in the shape of the Prądnik River, and the riverbed was displaced by 25 m. In another three areas, the Prądnik River did not change its shape; the displacements of watercourse do not exceed 1 m, and are most likely caused by differences in the water level and the orthophotomap's geolocation accuracy. There have been no other changes observed in the shape of the watercourse after 2003. The shape on the orthophotomap from 2003 is more irregular and the displacement reaches 1.5 m. In the southern part, we identified an area where the Prądnik changed its shape after 2003, it is less meandering, and these changes reach 1.5 m. In addition, one area where the Prądnik River increased its meanders after 2009 and before 2017 was identified. These changes reach 6 metres. Unfortunately, the area is not visible on the latest orthophotomap. The remaining parts of the river, visible on two or more images, do not show any significant changes. . The characteristics and placement of these changes, in the middle of a farmland, suggest that they were manmade. There was a significant change in the shape of the Prądnik River, and the riverbed was displaced by 25 m. In another three areas, the Prądnik River did not change its shape; the displacements of watercourse do not exceed 1 m, and are most likely caused by differences in the water level and the orthophotomap's geolocation accuracy. There have been no other changes observed in the shape of the watercourse after 2003. The shape on the orthophotomap from 2003 is more irregular and the displacement reaches 1.5 m. In the southern part, we identified an area where the Prądnik changed its shape after 2003, it is less meandering, and these changes reach 1.5 m. In addition, one area where the Prądnik River increased its meanders after 2009 and before 2017 was identified. These changes reach 6 metres. Unfortunately, the area is not visible on the latest orthophotomap. The remaining parts of the river, visible on two or more images, do not show any significant changes.

Vector Layers Analysis
In this study, a comparison between vector layers representing the Prądnik River was performed. Since the ONP layer was created using manual vectorisation on a high resolution DTM by the Ojcowski National Park employees, it was used as a base layer and other layers were compared to it. The ONP layer is up to date as of 2012.

Vector Layers Analysis
In this study, a comparison between vector layers representing the Prądnik River was performed. Since the ONP layer was created using manual vectorisation on a high resolution DTM by the Ojcowski National Park employees, it was used as a base layer and other layers were compared to it. The ONP layer is up to date as of 2012.

Current Data Analysis
The MPHP database vector layer was the least accurate of all analysed layers (Figure 11). It was analysed only visually, due to its low accuracy in comparison to the DTM. The reason for the low accuracy of this layer is the use of topographic maps at 1:50,000 scale as a basis for its creation. Other layers were compared in qualitative and quantitative terms. A qualitative comparison concerned determining the placement of the main course of the Prądnik River and its lateral watercourses. The results of the qualitative analysis are presented in Figure 12. The Prądnik River was divided into 12 parts: eight parts of the main watercourse and four of the lateral watercourses. The classification was done based on data recorded directly in the database (ONPw, BDOT10k). The layer had clear symbols for main and lateral watercourse, which differ in width (Topo10k92). On selected occasions only, the symbols were partially unambiguous and required a further map-based inspection (Topo10k65 and WIG25-some parts of the main watercourse were signed on the map, some prongs were marked because of the mill symbol on the map next to them). On the WIG25k map, all four lateral watercourses were marked, but two of them only schematically, due to limited amount of space on the map. The MPHP database vector layer was the least accurate of all analysed layers ( Figure 11). It was analysed only visually, due to its low accuracy in comparison to the DTM. The reason for the low accuracy of this layer is the use of topographic maps at 1:50,000 scale as a basis for its creation. Other layers were compared in qualitative and quantitative terms. A qualitative comparison concerned determining the placement of the main course of the Prądnik River and its lateral watercourses. The results of the qualitative analysis are presented in Figure 12. The Prądnik River was divided into 12 parts: eight parts of the main watercourse and four of the lateral watercourses. The classification was done based on data recorded directly in the database (ONPw, BDOT10k). The layer had clear symbols for main and lateral watercourse, which differ in width (Topo10k92). On selected occasions only, the symbols were partially unambiguous and required a further map-based inspection (Topo10k65 and WIG25--some parts of the main watercourse were signed on the map, some prongs were marked because of the mill symbol on the map next to them). On the WIG25k map, all four lateral watercourses were marked, but two of them only schematically, due to limited amount of space on the map.  Interestingly, during the field measurements, the main watercourse in part 1 (Figure 12) was completely dry (Figure 8(S 78)), and the water flowed only in the millrace 1a. The convergence of the Topo10k65 layer with the ONPw layer is significant for parts 1/1a and 3/3a. This convergence may result from the map interpretation mentioned above. Figure 12 shows the magnitude of difference between the layout of the main and lateral watercourse originating from various data sources.
A quantitative analysis was performed with reference to the ONPw layer. Table 5 contains a comparison of the vector data in terms of total sinuosity (TS). This is a parameter based on the coefficient between the length of the riverbed and the shortest distance between its beginning and end [83]. For a straight part TS equals 1.    The TS parameter was used as a measure of the detail of the Prądnik River shape. The results clearly show that only on the BDOT10k layer is the shape is close to the ONPw layer. The topographic maps have a similar level of detail, but clearly lower than the BDOT10k. The shape of the river is the least detailed in the WIG25k layer. Notable are parts 3 and 5 (Table 5) with a high value of TS. These are the fragments of the Prądnik River, where its riverbed often changed (see orthophotomaps analysis), it could also happen that the observed river meanders were intentionally exaggerated by the map editor (Topo10k65). Another qualitative analysis of the data was based on comparing the distance of the Prądnik River parts to the ONPw layer (Table 6). On every ONPw layer object, points were generated in increments of 1 meter. Then, the distances from these points to the corresponding objects on the other layers were calculated. The median and maximum were calculated for these values. Looking at the median values, we can see that, for the BDOT10k layer, for most parts it is below 1 m, which indicates a high similarity of these vectors. The maximum values for some parts sections of the BDOT10k layer may indicate some local differences with the ONPw layer. Based on the median value, only for Section 5 can one see a large discrepancy between the vectors. Table 6 shows that the representation of the Prądnik River on 1:10,000 topographic map in the PUWG 1965 projection (older ones) is more similar to the ONPw layer than the topographic map in the PL-1992 projection (more recent ones). The values from Table 6 also confirm the lowest quality of the WIG25k layer.

Old Maps Analysis
The results of the comparison are as follows: • Josephine map of Galicia-The map presents the general course of the river and its main turns.
Little meanders that are shown on the map are a type of a symbol of the river, and do not represent actual shape of the object. The largest difference between the spatial placement of the river course on the Josephine map and the ONP watercourses vector is around 170 m (Table 7).

•
Map of Ojców surroundings-The largest observed difference between the river course on the Ojców Bazaar map and the ONP watercourses vector is 100 m (Table 7). However, for most of the river course, this distance is much smaller. The shape of the river shows its main course but is more generalised and does not present smaller meanders. • Russian and Soviet "two-verst"-The largest difference between the river course on the Russian "two-verst" map and the ONP watercourses vector is 50 m ( Table 7). The differences between these two maps are particularly visible in places where the river has two courses. The second watercourse was a manmade structure created for diverging part of the water from the river to places that needed it, in this case for the use of a watermill, thus, it is classified as a millrace. • WIG25k-The largest difference between the river course on the WIG25k map and the ONP watercourses vector is around 30 m (Table 7). Overall, the WIG25k proved to be the most accurate map. The course of the Prądnik River is shown in detail, and it is possible to compare it locally with modern maps. • Godfryd Ossowski map-Due to the fact that the map range did not overlap with other maps used for this study, it was analysed only qualitatively. The general shape of the riverbed can be seen on this map. A particularly useful feature of the Ossowski map is the location of no longer existing artificial river forks associated with mills. The fragment shown in Figure 5 shows two mills and associated millraces, and one millrace without a mill.  Figure 13 shows the general change of the shape and placement of the Prądnik River course as represented on old maps. Due to the accuracies of maps (Table 7), and different methods of data acquisition during their creation, we cannot analyse smaller changes of the river course. We can, however, see that some significant changes occurred over time. This is not surprising, since it was the only source of water in an otherwise isolated part of the country, and inhabitants of Ojców used the Prądnik River as a source of energy for the mills. In addition, the Prądnik is not a deep river and with its narrow, isolated flood area it was easy to redirect it after every flood. This could also have happened, due to environmental reasons.
Remote Sens. 2020, 12, x FOR PEER REVIEW 22 of 30 its narrow, isolated flood area it was easy to redirect it after every flood. This could also have happened, due to environmental reasons. The obvious conclusion for all old maps of this area, is that rivers represented on older maps with smaller scales differ more from the ONP watercourse layer than newer and larger scale maps of the same region. Moreover, the differences between vectors from each map and the ONP watercourse layer are not systematic and differ in size. Old maps present only the main course of the river, they show curve sequence and the degree of regulation of the river, but the accuracy of those maps does not allow for comparison with the present river course. Thus, it is impossible to use them for change detection of location and shape of the Prądnik River over time. However, the Godfryd Ossowski map can be used to determine location and change of manmade millraces, which were rebuilt many times. The research confirms the fact known from the literature that old maps can be compared with contemporary data only to some degree [84].

Conclusions
The main aim of this study was to evaluate the previous and current shape of the Prądnik River-a small alluvial river surrounded by significantly higher limestone outcrops and wood The obvious conclusion for all old maps of this area, is that rivers represented on older maps with smaller scales differ more from the ONP watercourse layer than newer and larger scale maps of the same region. Moreover, the differences between vectors from each map and the ONP watercourse layer are not systematic and differ in size. Old maps present only the main course of the river, they show curve sequence and the degree of regulation of the river, but the accuracy of those maps does not allow for comparison with the present river course. Thus, it is impossible to use them for change detection of location and shape of the Prądnik River over time. However, the Godfryd Ossowski map can be used to determine location and change of manmade millraces, which were rebuilt many times. The research confirms the fact known from the literature that old maps can be compared with contemporary data only to some degree [84].

Conclusions
The main aim of this study was to evaluate the previous and current shape of the Prądnik River-a small alluvial river surrounded by significantly higher limestone outcrops and wood covered stacks. This river was used as an example since there are many objects of this kind that pose threat to the safety of people and infrastructure due to occurring floods. The monitoring of these areas is often done with the use of remote sensing techniques such as LiDAR, airborne mapping and GNSS surveys, and products derived from this techniques, such as DEMs based on point clouds, orthophotomaps, and cross-section profiles. In this study, we evaluated the quality and completeness of these data sets. We also proposed how remote sensing data could be enriched with other techniques, in order to create a full picture of the current situation of any small watercourse and its surroundings.
The results of this study show problems with current approaches to data acquisition via remote sensing methods. Here, we list them and add information on how diagnosed problems could be mitigated.
All-purpose LiDAR based point clouds do not give the full information on the shape of the watercourse, due to the problems with data acquisition in wooded areas and underneath plant life. This can be mitigated by customising the flight plan, to be performed in early spring, or by using the TLS and SfM instead, with customised flight plans.

•
LiDAR based point clouds do not show the riverbed, due to the fact that the laser beam does not go through water. This problem can be mitigated by using bathymetric scanners or providing additional terrestrial land surveys of the riverbed.

•
Orthophotomaps in mountain areas do not provide sufficient information on the shape of the riverbed due to acquisition problems. All-purpose orthophotomaps are often made throughout the year, which means that some areas will be covered with deciduous trees obscuring the river. Additionally, since there are no general guidelines for the time of the day when the flight should take place, large areas can be shaded. This problem can be mitigated with either customised flight plans, or by using UAV for providing data for maps, since this type of an aircraft needs less preparation and planning.

•
Providing data for detailed cross-section profiles, needed in a more precise analysis, cannot often be done by a GNSS based remote sensing method, due to the problems with signal, satellite availability and, at times, multipath error. This can be mitigated by using a traditional tachymetric survey or a hybrid approach, with some elements measured via GNSS (traverse points, cross-sections where the signal is available) and others via tachymetry.

•
In mountain areas, all-purpose remote sensing-based data on small and narrow objects are more susceptible to systematic and non-linear errors. They should be evaluated with regard to their quality and consistency before being used for any analysis. This should be done by performing a ground truth high accuracy survey of highly distinguishable elements around the river or cross-section profiles. These ground truth data can be used for statistical analysis.

•
In most areas, remote sensing-based approaches have only been available for the last 20 years or so. Flood monitoring processes often required longer time steps, and this lack of data can be, to some extent, filled with available archival data, such as river profiles and old maps.
Evidence of previously described conclusions can be seen in more detail while analysing the case study of the Prądnik River. The evaluation of available data is given in Table 8.   1978-1979/1983-1986 Moderate, it is possible to compare shape and order of meanders Restricted conditions on the study area, i.e., a small river within the area of the National Park, surrounded by the valley with high, forested and steep slopes limits the possible use of LiDAR based techniques and forces the use of a traditional survey measurement. Only this method allows for providing sufficient accuracy for the analysis of changes in the shape of the riverbed. In addition, none of the analysed remote sensing data sets, contain information on the depth of the river. Therefore, they cannot be effectively compared to another data sets.
Orthotomaps of the Prądnik River valley can be regarded as partially useful. They show parts of the river and its surroundings, which helped to recognise significant changes fast (Figure 10 bottom row). However, it was impossible to analyse the whole river, due to the trees and slopes that casted shadows. Furthermore, it was difficult to carry out a multitemporal analysis based on the orthophotomaps, as different parts of the river were visible on particular data sets. The orthophotomaps could be used in a detailed analysis of small portions of the river after a prior visibility assessment. They would be more useful during the study of changes based on the measurements conducted along the riverbanks, rather than cross-sectional profiles.
Among old and modern topographic maps, one can find few that are useful. Maps and databases maintained at scales of 1:10,000 to 1:25,000 seem to be useful for conducting analyses, however, modern 1:10,000 maps confirm that not every data set in these scales allows for detailed shape analyses. Old maps are useful primarily for analysing the general shape of the river. The exceptions are the high quality Ossowski map and the WIG25k map.
The analysis of data from the ALS showed that these products were created using general-purpose algorithms. Both point clouds and terrain models do not allow for the use of only these data in the high accuracy analysis of the Prądnik riverbed shape. After enriching these data with a rough course of the river (e.g., derived from the BDOT10k), they can be used for automatic detection of the riverbed.
The analysis allowed for drawing the following conclusions: • Since the 19th century, backfilling mills use has been progressing, which locally changed the course of the main riverbed of the Prądnik River.

•
In some areas, the Prądnik River changed its course by several meters. Smaller changes in the shape of the riverbed cannot be determined using available historical data.

•
Due to the accuracy of the source data, it is difficult to say with certainty whether this change was dictated by natural or anthropogenic factors. The studied area before the creation of the National Park was used as agricultural and residential areas. This suggests that some changes are caused by human activities.

•
In several places associated with millraces, the mainstream of the Prądnik River swapped with the side stream, and vice versa.
Apart from the data set created for this paper, there is no currently widely available data set, or a combination of data sets that describes the Prądnik River in enough detail to be used for proper flood mitigation planning.