Geospatial Approaches to Monitoring the Spread of Invasive Species of Solidago spp.

: Global climate change inﬂuences plant invasion which spreads all over the Europe. Invasive plants are predominantly manifest negative impacts, which require increased attention not only from ecologists. The research examines the possibilities offered by geospatial technologies in mapping the spatial spread of invasive plants of the genus Solidago . Invasive plant population was investigated at two localities, Mal ý Šariš and Chminianska Nov á Ves in Slovakia, as well as the mapping of the area by multispectral imaging to determine the spectral reﬂectance curve of the monitored plant species. Using spatial analyses in the geographic information system, we evaluated changes in the plant density in the two localities. Based on the obtained results, we found that the number of individuals (ramets) in the Mal ý Šariš is signiﬁcantly increasing, while in the examined area of Chminianska Nov á Ves, there is a decrease in the number of Solidago spp. in the last monitored year. At the same time, we can state that in the areas with the highest increase in the number of ramets, the highest plant density per hectare was also recorded. We can also say that due to the spectral proximity of the surrounding vegetation, the spectral resolution in four spectral bands is insufﬁcient for the classiﬁcation of multispectral records in the case of Solidago spp. and cannot replace the advantages of high spectral resolution hyperspectral imaging, which signiﬁcantly reﬁnes the feature space for Solidago spp. and the surrounding vegetation.


Introduction
The ecological impact of invasive species has been observed in all types of ecosystems.Typically, invaders can change the niches of co-occurring species, alter the structure and function of ecosystems by degrading native communities and disrupt evolutionary processes through anthropogenic movement of species across physical and geographical barriers [1][2][3][4][5].Recent invasion research ranges from developing testable hypotheses aimed at understanding the mechanisms of invasion to providing guidelines for the control and management of invasive species.Early detection and mapping of the extent of rapidly spreading invasive populations are critical for informing management priorities, including eradication efforts [6].Conserving biodiversity, however, requires data.The ability to monitor the state of our natural resources and the impacts of human activity on these resources is fundamental to designing appropriate and optimized management strategies.The scientific community needs to be able to access global, long-term, reliable information on spatial temporal changes in the distribution of direct and indirect anthropogenic pressures to biological diversity; in the distribution, structure, composition, and functioning of ecosystems; as well as evidence of the effectiveness of various management actions [7].Biodiversity monitoring experts necessitates them to work at the interface between ecology and remote sensing.Field surveys are costly in terms of time and manpower, and georeferenced records of most species are inadequate for distribution modelling [8].The development of remote sensing technology has made it possible to map the distribution of plant species and aims to extract information about the condition and/or state of an object without requiring physical contact [9][10][11].Understanding species distributions, and being able to predict their potential shifts, relies on being able to characterize the set of conditions that drive occupancy.Doing so requires information on species' presence, and on the environmental conditions and resources available to the organisms of concern.Deriving maps detailing the bio-physical cover of the Earth's surface from satellite remote sensing (SRS) is among the oldest applications of satellite data to inform applied ecology with, for instance, the first global land cover products produced in the 1990s [12,13].New combinations of satellite remote sensing (SRS) data with methodologies such as support vector machines [14] and multisensory image fusion [15] are constantly tested for application, while investigations for developing new indices to monitor vegetation have been occurring for decades [16,17].Although SRS is increasingly popular in applied ecology, the data are still underused.Detection of specific animal and plant individuals using information collected by satellites is still not commonly achieved, although successful examples have started to accumulate [7,18].Moderate (<100 m ground sampling) resolution satellite images have for example been used to map the spatial extents of cheatgrass Bromus tectorum [19].GeoEye-1 satellite imagery (spatial resolution of 1.65 m) was suitable for the detection of medium-to large-sized wildlife as well as informing population estimation exercises in open savanna [20].Very high spatial resolution (pixel size <10 m 2 ) commercial optical sensors have provided new opportunities for habitat mapping at a finer spatial scale than previously possible [21,22].Most biodiversity researchers and conservation practitioners lack the financial resources to acquire large amounts of SRS data, especially if they need information over wide geographical extents or long periods of time [23].There are many underused products and algorithms that could benefit environmental and wildlife management (e.g., fractional vegetation cover products, radar and hyperspectral imagery, differential interferometry algorithms).Most mapping studies are performed at the local scale (100s to 1000s of hectares), while distribution modelling is typically performed at global, continental, national, or regional scales [24].Researchers and managers have embraced aerial photography [25][26][27][28][29][30][31] for mapping the distribution of invasive weeds.The potential for SRS to support natural resource management has never been greater; however, several issues limit this potential.First, although many satellite products are freely available, a significant proportion of products are not [32].Further, SRS-based data analyses can be expensive given logistical requirements (i.e., hardware, software, qualified staff and training) for the processing and analysis of large data sets.Altogether, costs can be considerable and are hampering widespread application of satellite monitoring in applied ecology and management [23,32,33].In addition to free SRS products, open source software solutions, such as R [34], QGIS [35] or GRASS GIS [36] are on the rise.
Present research considered equipment's which are more available for the investigation of the invasive plants Solidago spp.spreading in low acreage localities.There were used precise positioning of plants by RTK GNSS equipment during four growing seasons and multispectral imaging using UAV in last monitored season for remote sensing.

Monitored Plant Species
Solidago canadensis and Solidago gigantea (Solidago spp.) are dangerous plant invaders in Europe, which suppress the indigenous flora.Both species originated from North America and were spread over the Europe and Asia as beautiful ornamental flowers.They are perennial herbs.Morphology of both species is very similar.Both reach heights about 1.5 to 2 m tall, sometimes spreading by means of underground rhizomes.Solidago spp.are clonal plants.Populations with clonal growth have difficult structure and it is important to recognize the genets and ramets.Ramete is a basic clonal unit and its synonym is individual plant.It was vegetative reproduced from the clone plant and is physiologically independent.Genete is a plant colony with the same genetic, where it is not possible to recognize mother plant.It is synonymous with a bunch of plants.
It often grows from genete with no leaves at the base but numerous leaves on the stem.At the top, each shoot produces a sizable array of many small yellow flower heads, sometimes several hundred.The seeds are dispersed by wind; within established stands, clonal growth of the species can form dense stands.It is possible to find them in a wide variety of natural habitats, although it is restricted to areas with at least seasonally moist soils.

Monitored Localities
Localities were selected based on the following criteria: (1) locality which is freely accessible; (2) locality without management; (3) locality with appearance of Solidago spp. with potential to spread.There were four selected localities in the sub-urban area of Prešov (Slovakia), but after two years of observations, two of them were excluded based on unexpected management provided on them.The two remaining localities (Malý Šariš and Chminianska Nová Ves) are localized in the Western part of Prešov (Figure 1).
America and were spread over the Europe and Asia as beautiful ornamental flowers.Th are perennial herbs.Morphology of both species is very similar.Both reach heights abo 1.5 to 2 m tall, sometimes spreading by means of underground rhizomes.Solidago spp. a clonal plants.Populations with clonal growth have difficult structure and it is importa to recognize the genets and ramets.Ramete is a basic clonal unit and its synonym is in vidual plant.It was vegetative reproduced from the clone plant and is physiologica independent.Genete is a plant colony with the same genetic, where it is not possible recognize mother plant.It is synonymous with a bunch of plants.
It often grows from genete with no leaves at the base but numerous leaves on t stem.At the top, each shoot produces a sizable array of many small yellow flower hea sometimes several hundred.The seeds are dispersed by wind; within established stan clonal growth of the species can form dense stands.It is possible to find them in a wi variety of natural habitats, although it is restricted to areas with at least seasonally mo soils.

Monitored Localities
Localities were selected based on the following criteria: (1) locality which is fre accessible; (2) locality without management; (3) locality with appearance of Solidago sp with potential to spread.There were four selected localities in the sub-urban area Prešov (Slovakia), but after two years of observations, two of them were excluded bas on unexpected management provided on them.The two remaining localities (Malý Ša and Chminianska Nová Ves) are localized in the Western part of Prešov (Figure 1).

Locality of Malý Šariš
The monitored locality of Malý Šariš is situated 5 km west of Prešov and 289 m abo sea level.The locality has a triangle shape with an area of 0.81 ha.Two sites (South-W and South-East) are bordered by field, while a third site (North) neighbors them with t road number E18.The research locality is oriented to the north, predominantly flat wit

Locality of Malý Šariš
The monitored locality of Malý Šariš is situated 5 km west of Prešov and 289 m above sea level.The locality has a triangle shape with an area of 0.81 ha.Two sites (South-West and South-East) are bordered by field, while a third site (North) neighbors them with the road number E18.The research locality is oriented to the north, predominantly flat with a slight slope of up to 2 • .The locality consists of diluvial sediments as well as claystone's and sandstones.There are fluvisols and planosols occur in the research area.

Locality Chminianska Nová Ves
Monitored locality Chminianska Nová Ves is situated 11.7 km west of the Prešov in the 362 m a.s.l.The area of 2.84 ha has a western slope of 11 • .The site has an irregular longitudinal shape in a west-east direction.It is bordered by the highway D1 in the south and dense banks with the vegetation of the Svinka river in the north.The locality is divided in the middle by a small section of dense shrub vegetation.The locality consists of mediumand coarse-grained sandstones but also a flysch, a sequence of sedimentary rock layers of claystone's and sandstones.It is possible to find mainly cambisols and dystrict leptosols and, to a lesser extent, fluvisols around the Svinka watercourse.

Monitoring of the Solidago spp. Spreading by Using GNSS
Empirical research consists of recording the location of Solidago's genets.Direct positions of the Solidago genets were recorded during the four-year period (2016-2019).Monitoring was provided in the period of plant flowering, when the yellow flowers of Solidago spp.Are easy to recognize (Table 1.).There was recorded position of the genets as well as the number of the shoots in the genets.Plant location recording was performed via the Sygic mobile application, version 11.2.6, in the first monitoring year in 2016.In the next three years (2017-2019) we used RTK GNSS (Real-Time Kinematic Global Navigation Satellite System-a type of differential GNSS) instrument to locate the genets of Solidago spp., followed by using position corrections from the Slovak Spatial Observation Service (SKPOS).We were able to record the position of ramets and shoots with an accuracy of less than 2 cm.The coordinates of plants were recorded in a text document from where we imported them and processed as a point vector layer in the geographic information system QGIS 3.Each point of the vector layer that represented the genet also contained information on the number of shoots.There was plotted the occurrence and number of plants cartographically, in all mappings.The obtained data were divided into 6 categories using the cartogram method and each genet was marked with a different point size and color depending on the number of ramets.

Analysis of the Changes in Solidago spp. Spreading by Using Geospatial Analysis in GIS
We used spatial analysis of GIS to evaluate changes in the distribution of Solidago spp.during the three years-2017, 2018 and 2019-at both localities.For the first year (2016), analysis was omitted due to the fact that in that year was the position of individuals recorded by using the Sygic mobile application, which caused the measurement inaccuracy to reach several meters.

Kernel Analysis of Plant Density
The first step in the analysis was to calculate the standard distance for each layer with the function spatial point pattern analysis by tools of SAGA GIS software (System for Automated Geoscientific Analyses).Subsequently, we used the formula: where (h) is used to calculate the radius based on the standard distance (σ) and the number of the noticed points (n) [37].Reached radius was then used in the function Kernel density estimation.Plant density was then calculated based on focal map algebra with the using focal average, while the result was the density of plant shoots at squared meter.Using a raster calculator, which work on the principle of local map algebra, the value of each cell was converted to the number of plants per hectare.In this way, we determined the plant density for all three monitored years in both mentioned localities.

Evaluation of Spectral Reflectance of Selected Invasive Plants Using Multispectral Drone Imaging
We used methods and remote sensing tools for the monitoring at the locality of Malý Šariš within the vegetation season 2019.The aim of the approach was to investigate the possibilities of the mapping of invasive Solidago spp.Spreading by using multispectral image collection following by supervised classification.Due to the small size of the monitoring locality as well as the low penetration of Solidago spp.Spread, scanning with a sufficiently high spatial resolution of the resulting multispectral images was necessary.For this purpose, we used the configuration of the system of unmanned aerial vehicle DJI Matrice 200, which carried a Parrot Sequoia multispectral camera via a gyroscopic gymbal.This minimized sensing errors due to vibrations that occur during drone flight.The Parrot Sequoia camera is primarily intended for monitoring vegetation by scanning in four spectral bands (green-550 nm, red-660 nm, red edge-735 nm and near-infrared-790 nm), which are most suitable for the analysis of the reflective properties of vegetation.The camera also has a standard RGB camera, which allows for image collection in the visible spectrum of electromagnetic radiation and subsequent creation of orthophotos.The camera captured the locality of Malý Šariš according to a predefined image flight from a height of 30 m, which in the end resulted in final images with a spatial resolution of 3 cm.Taking pictures also required a series of ground measurements aimed at positioning artificially signaled ground control points via RTK GNSS device of the Topcon brand Hiper HR.Position information on plants was also used as part of the ground measurements of Solidago spp.from 2019.During the image flight, 350 images were created for each spectral band as well as images in the visible spectrum.These images were subsequently processed via Agisoft Photoscan software.This software uses the method Structure from Motion to derive a cloud of 3D points using central raw projection images and their mutual overlap defined when planning the parameters of the image flight.This point cloud software was used to interpolate the digital elevation model which used images of the area in the final orthorectification.In the process, multispectral images are radiometrically calibrated by means of a calibration scan of the target before the actual image collection, to eliminate the interfering effects of reflections from the surroundings of the scanned area, which distort the real reflectivity of the area of interest.Finally, multispectral images have a radiometric resolution of 16 bits, which allows us to split the values of reflected electromagnetic radiation up to 65536 levels.Scanning monitoring of the locality of Malý Šariš was performed in two time periods, 14 June 2019, i.e., at the beginning of the Solidago spp.growing season, and 23 August 2019, i.e., on its blooming phases.The above-mentioned time periods was chosen for the comparison of reflection properties of Solidago spp. in different stages of its vegetation cycle.To determine the spectral reflectance, a network of points was used to which the reflectance level values were taken electromagnetic radiation via the Add raster values to points functioning the environment of the geographic information system QGIS.In the case of the Solidago spp., it was an irregular network of points from location data aimed through RTK GNSS.In the case of other vegetation, we covered the area of the monitoring site with a regular network point with a spacing of 0.5 m.These points were divided into tree or shrub class vegetation, spontaneous grass vegetation and mown grass vegetation.We monitored 1258 shoots of Solidago spp. in 2016 in Malý Šariš.After calculating, it was determined that there were 1553 shoots per hectare.Individual plants, called ramets, numbered 37 and the rest were unevenly placed in 221 genets, while the lowest number of shoots in genet was two and the highest number was 32 (Figure 2a).In the next year, 2017, we evaluated the increased number of Solidago spp.Shoots, which was 1865.From this number, there were 161 ramets and 265 genets.The highest number of shoots in genet was 48 and the lowest number was two.After calculating, there were 2331 shoots per hectare (Figure 2b) [38].In the third monitoring year, 2018, we evaluated 2523 shoots.A total of 180 plants were noted as ramets, and 2343 plants were a part of 412 genets.The shoots are in genets distributed unevenly, while the lowest number of shoots in genet was two and the highest number was 50 (Figure 2c).In the year 2019, we evaluated 3934 shoots of Solidago spp., which presented 4856 shoots per hectare.We noticed 324 ramets and 661 genets with a different number of shoots.The minimum number of shoots in genet was two and the maximum was 42 (Figure 2d).After four-year monitoring of Solidago spp., we can notice the trifold number in the locality of Malý Šariš.This is very dangerous for the native flora.In the figures below, the monitored area is bordered by the yellow color, ramets are marked with the red triangle and genets based on the number of the shoots by different color.The arrow indicates the expected direction of distribution of the Solidago spp.population from the source site [39].

Mapping the
Remote Sens. 2021, 13, x FOR PEER REVIEW 6 of 20 electromagnetic radiation via the Add raster values to points functioning the environment of the geographic information system QGIS.In the case of the Solidago spp., it was an irregular network of points from location data aimed through RTK GNSS.In the case of other vegetation, we covered the area of the monitoring site with a regular network point with a spacing of 0.5 m.These points were divided into tree or shrub class vegetation, spontaneous grass vegetation and mown grass vegetation.We monitored 1258 shoots of Solidago spp. in 2016 in Malý Šariš.After calculating, it was determined that there were 1553 shoots per hectare.Individual plants, called ramets, numbered 37 and the rest were unevenly placed in 221 genets, while the lowest number of shoots in genet was two and the highest number was 32 (Figure 2a).In the next year, 2017, we evaluated the increased number of Solidago spp.Shoots, which was 1865.From this number, there were 161 ramets and 265 genets.The highest number of shoots in genet was 48 and the lowest number was two.After calculating, there were 2331 shoots per hectare (Figure 2b) [38].In the third monitoring year, 2018, we evaluated 2523 shoots.A total of 180 plants were noted as ramets, and 2343 plants were a part of 412 genets.The shoots are in genets distributed unevenly, while the lowest number of shoots in genet was two and the highest number was 50 (Figure 2c).In the year 2019, we evaluated 3934 shoots of Solidago spp., which presented 4856 shoots per hectare.We noticed 324 ramets and 661 genets with a different number of shoots.The minimum number of shoots in genet was two and the maximum was 42 (Figure 2d).After four-year monitoring of Solidago spp., we can notice the trifold number in the locality of Malý Šariš.This is very dangerous for the native flora.In the figures below, the monitored area is bordered by the yellow color, ramets are marked with the red triangle and genets based on the number of the shoots by different color.The arrow indicates the expected direction of distribution of the Solidago spp.population from the source site [39].

Mapping the Spread of Solidago spp. via GNSS at the Chminianska Nová Ves Locality
There were monitored 116 shoots of Solidago spp. in 2016 at the locality Chminianska Nová Ves.After calculation, it was determined that there were 41 shoots per hectare.We evaluated eight genets and no ramets (Figure 3a).The minimum number of shoots was two and maximum was 61.During the second monitoring year in 2017, we evaluated increased number of shoots 156.We noticed three ramets and 13 genets with minimum two shoots and maximum 85 shoots.This presents 54.9 shoots per hectare (Figure 3b).In August 2018, we evaluated 244 shoots of Solidago spp.; none were evaluated as ramet.Plants were in 17 genets, in which the minimum number of shoots was 2 and the maximum number of shoots was 100.The number of shoots that occurred in this locality after calculating increased by 85 per hectare (Figure 3c).In the last year of mapping, 2019, Solidago spp.spread was evaluated to be 209 shoots, from which, seven were ramets.The rest of the plants were distributed in the 24 genets with a maximum of 41 shoots (Figure 3d).In the figure below, the monitored area is bordered by the yellow color, while ramets and genets are marked with different colors depending on quantity.
Remote Sens. 2021, 13, x FOR PEER REVIEW 7 of 20 There were monitored 116 shoots of Solidago spp. in 2016 at the locality Chminianska Nová Ves.After calculation, it was determined that there were 41 shoots per hectare.We evaluated eight genets and no ramets (Figure 3a).The minimum number of shoots was two and maximum was 61.During the second monitoring year in 2017, we evaluated increased number of shoots 156.We noticed three ramets and 13 genets with minimum two shoots and maximum 85 shoots.This presents 54.9 shoots per hectare (Figure 3b).In August 2018, we evaluated 244 shoots of Solidago spp.; none were evaluated as ramet.Plants were in 17 genets, in which the minimum number of shoots was 2 and the maximum number of shoots was 100.The number of shoots that occurred in this locality after calculating increased by 85 per hectare (Figure 3c).In the last year of mapping, 2019, Solidago spp.spread was evaluated to be 209 shoots, from which, seven were ramets.The rest of the plants were distributed in the 24 genets with a maximum of 41 shoots (Figure 3d).In the figure below, the monitored area is bordered by the yellow color, while ramets and genets are marked with different colors depending on quantity.Within the monitored localities, we learned several facts by comparison of data obtained by mapping during the years 2016 to 2019.The number of shoots increased by 48.3% at the monitored locality of Malý Šariš, between 2016 and 2017.As can be seen in Table 2, a rapid quadruple transmission was recorded in the numbers of ramets, but in comparison with the numbers of genets, an increase of only 19% was achieved [36].A comparison of data from 2018 and 2019 showed an increase in the number of shoots by 80% from 180 to 324, but also an increase in the number of genets by 60%.During the four monitored years, we can observe that the number of all recorded shoots increases up to threefold.Within the Chminianska Nová Ves locality, a lower number of Solidago spp.shoots were recorded compared to the locality of Malý Šariš.During the first two monitored years, the total number of shoots increased by 34.5%.The highest plant growth was Within the monitored localities, we learned several facts by comparison of data obtained by mapping during the years 2016 to 2019.The number of shoots increased by 48.3% at the monitored locality of Malý Šariš, between 2016 and 2017.As can be seen in Table 2, a rapid quadruple transmission was recorded in the numbers of ramets, but in comparison with the numbers of genets, an increase of only 19% was achieved [36].A comparison of data from 2018 and 2019 showed an increase in the number of shoots by 80% from 180 to 324, but also an increase in the number of genets by 60%.During the four monitored years, we can observe that the number of all recorded shoots increases up to threefold.Within the Chminianska Nová Ves locality, a lower number of Solidago spp.shoots were recorded compared to the locality of Malý Šariš.During the first two monitored years, the total number of shoots increased by 34.5%.The highest plant growth was recorded in 2018, when the number of shoots increased by 56.4%.On the other hand, in 2019, we found that the number of shoots decreased by 14.3%.The highest density of the investigated plant species was in the range from 8500 to 9000 plants per hectare at the locality of Malý Šariš.In the year 2017, we evaluated the density of up to 500 plants per hectare.As seen in Figure 4a, the highest density was noticed in the Northern site close to trees vegetation.The following year, we evaluated increased plant density up to 1500 plants per hectare.In Figure 4b are visible areas with orange and red color, which represents the density of more than 3500 plants per hectare.The most visible changes in plant density were noticed in 2019 (Figure 4c) when we evaluated a decrease in the areas with lighter color and an increase in dark colors.The areas with plant density higher than 4500 plants per hectare increased, which are highlighted in red.We noticed areas in the South-West of the investigated locality with increased occurrence of the Solidago spp. in 2019.Before this year, the locality was affected by human activity and no invasive plants were spread.Plant density decreased in the North-East and in central area of the investigated locality, where the tree and shrub vegetation is located.We, therefore, assume that higher vegetation may represent a certain barrier to the transmission and spread of seeds by wind from the source site, as well as over dense of different plant species could be not impenetrable for the genets spreading.In this part of the investigated area, there is a greater capture of seeds, which are transmitted by cars passing on the adjacent road.Another area of increased density of Solidago spp. is located along the eastern boundary of the site, where there is a field road separating our research area from the so-called seed bank, i.e., the source site, which means that the released seeds do not have to overcome any significant obstacle to their spread.

Analysis of Solidago spp. Changes in Propagation through Spatial Analyzes by GIS in Chminianska Nová Ves
Due to the lower number of recorded shoots, a lower plant density was found in the Chminianska Nová Ves locality than in Malý Šariš.Plant density ranged between 1 and 305 plants per hectare.It is possible to see that the density shown mainly in light yellow color, from which we can conclude that the density of plants in 2017 was mostly in the range of up to 30 individuals per hectare, while the genets were distributed almost over the entire monitored area (Figure 5a).A year later, in 2018, there were only very slight changes in the density of individuals (Figure 5b).The most significant changes in the density of Solidago spp. in this locality occurred in 2019.As we can see in Figure 5c, the light-colored areas with low plant density decreased significantly and, on the contrary, the orange-to red-colored parts of the area were added, where the density ranges from 100 plants per hectare and more.This year, we recorded a lower number of shoots at the locality compared to previous years, which was also reflected in the density of plants.The number of genets decreased, which were in the peripheral parts of the territory and thus also at a greater distance from the area where a higher concentration of genets was recorded.There is no source (seed bank) locality identified near this area from which Solidago spp.seeds would spread and, at the same time, the whole locality is bordered by tree and shrub vegetation, which prevents the spread of seeds to a large extent.The investigated locality is not maintained, and the herbaceous vegetation reaches dimensions of more than half a meter and there is no space for the successful growth of new ramets, which may represent a barrier to the spread of invasive species.Therefore, the number of shoots increases only in small parts, where the most shoots were recorded during all monitored years, i.e., it is the area with the highest plant density per hectare.Although the locality is bordered by the motorway, we can say that it does not significantly contribute to the spreading of Solidago spp., as was the case in the locality of Malý Šariš.

Analysis of Solidago spp. Changes in Propagation through Spatial Analyzes by GIS in Chminianska Nová Ves
Due to the lower number of recorded shoots, a lower plant density was found in the Chminianska Nová Ves locality than in Malý Šariš.Plant density ranged between 1 and 305 plants per hectare.It is possible to see that the density shown mainly in light yellow color, from which we can conclude that the density of plants in 2017 was mostly in the range of up to 30 individuals per hectare, while the genets were distributed almost over the entire monitored area (Figure 5a).A year later, in 2018, there were only very slight changes in the density of individuals (Figure 5b).The most significant changes in the density of Solidago spp. in this locality occurred in 2019.As we can see in Figure 5c, the lightcolored areas with low plant density decreased significantly and, on the contrary, the orange-to red-colored parts of the area were added, where the density ranges from 100 plants per hectare and more.This year, we recorded a lower number of shoots at the locality compared to previous years, which was also reflected in the density of plants.The number of genets decreased, which were in the peripheral parts of the territory and thus also at a greater distance from the area where a higher concentration of genets was recorded.There is no source (seed bank) locality identified near this area from which Solidago In addition to the calculation and representation of plant density for all three years, we also determined the difference in plant density between 2019 and 2017 for both selected localities.The difference in density is also visualized in figures by a color scale from white, which represents the lowest differences in density, to a dark color, which in turn shows the highest differences in density.Furthermore, by classifying the rasters, we calculated their acreage, which we were able to quantify the growth levels of Solidago spp. in these localities, the largest differences in plant density were found in parts near trees and shrubs in the locality Malý Šariš, which we also see in Figure 6.As mentioned above, we believe that higher plant stands prevent the transfer of seeds and ramets, which thus remain in place in front of such a barrier.At the same time, these parts are near the source site, from where we assume that the invasive species is spreading.We also observed that the areas with the highest density difference are identical to the areas where the largest increase in the number of shoots was recorded during the monitored years.The largest area (5327.32m 2 ) is occupied by areas where the plant density ranges from 1000 shoots per hectare.Conversely, areas with a density of over 6000 plants per hectare, which are shown in Figure 6 a distinctive blue color, occupy only 13.251 m 2 (Table 3).In addition to the calculation and representation of plant density for all three years, we also determined the difference in plant density between 2019 and 2017 for both selected localities.The difference in density is also visualized in figures by a color scale from white, which represents the lowest differences in density, to a dark color, which in turn shows the highest differences in density.Furthermore, by classifying the rasters, we calculated their acreage, which we were able to quantify the growth levels of Solidago spp. in these localities, the largest differences in plant density were found in parts near trees and shrubs in the locality Malý Šariš, which we also see in Figure 6.As mentioned above, we believe that higher plant stands prevent the transfer of seeds and ramets, which thus remain in place in front of such a barrier.At the same time, these parts are near the source site, from where we assume that the invasive species is spreading.We also observed that the areas with the highest density difference are identical to the areas where the largest increase in the number of shoots was recorded during the monitored years.The largest area (5327.32m 2 ) is occupied by areas where the plant density ranges from 1000 shoots per hectare.Conversely, areas with a density of over 6000 plants per hectare, which are shown in Figure 6 a distinctive blue color, occupy only 13.251 m 2 (Table 3).The highest difference in plant density was recorded in areas where we located several genets in the Chminianska Nová Ves locality, from which the seeds do not spread over long distances but remain near these clumps (Figure 7).The reason for spreading over such short distances may also be the fact that the original vegetation has created a  The highest difference in plant density was recorded in areas where we located several genets in the Chminianska Nová Ves locality, from which the seeds not spread over long distances but remain near these clumps (Figure 7).The reason for spreading over such short distances may also be the fact that the original vegetation has created a continuous dense stand, through which the seeds of Solidago spp.cannot penetrate the soil.We can, therefore, assume that during the following vegetation periods, the number of individuals will increase mainly in the vicinity of these genets and not in the rest of the territory.Areas with a density of up to 50 shoots per hectare occupy the largest part of the territory, specifically, 8144.4 m 2 (Table 3).Interestingly, however, in the lowest area, the density ranges from 150 to 200 shoots per hectare, while the areas with the highest plant density (more than 2000 shoots /ha) occupy an area of 293.53 m 2 .

Spectral Analysis of Solidago spp. Stands in the Locality of Malý Šariš
The increase in reflectance is higher from 1.5% for visible spectra (green and red) to 6% for border red and near-infrared radiation in the case of Solidago spp.(Figure 8).Significant differences in the reflectivity of Solidago spp.beyond the visible spectrum are identical with the general course of the reflectivity of vegetation, which is manifested by a higher reflectivity in the infrared part of the spectrum.For this reason, this spectral region is used in the analysis of vegetation properties.On the other hand, a larger range in reflectivity represents a complication in the automated classification of multispectral images, which manifests itself in the form of a higher probability when the Solidago spp.spectral space overlaps with other species of surrounding vegetation.

Spectral Analysis of Solidago spp. Stands in the Locality of Malý Šariš
The increase in reflectance is higher from 1.5% for visible spectra (green and red) to 6% for border red and near-infrared radiation in the case of Solidago spp.(Figure 8).Significant differences in the reflectivity of Solidago spp.beyond the visible spectrum are identical with the general course of the reflectivity of vegetation, which is manifested by a higher reflectivity in the infrared part of the spectrum.For this reason, this spectral region is used in the analysis of vegetation properties.On the other hand, a larger range in reflectivity represents a complication in the automated classification of multispectral images, which manifests itself in the form of a higher probability when the Solidago spp.spectral space overlaps with other species of surrounding vegetation.

Spectral Analysis of Solidago spp. Stands in the Locality of Malý Šariš
The increase in reflectance is higher from 1.5% for visible spectra (green and red) to 6% for border red and near-infrared radiation in the case of Solidago spp.(Figure 8).Significant differences in the reflectivity of Solidago spp.beyond the visible spectrum are identical with the general course of the reflectivity of vegetation, which is manifested by a higher reflectivity in the infrared part of the spectrum.For this reason, this spectral region is used in the analysis of vegetation properties.On the other hand, a larger range in reflectivity represents a complication in the automated classification of multispectral images, which manifests itself in the form of a higher probability when the Solidago spp.spectral space overlaps with other species of surrounding vegetation.By comparing the reflectivity of our classified vegetation classes in the monitored locality of Malý Šariš, we can state their relatively high spectral proximity (Figure 9).A more significant difference occurs only in the case of tree and shrub (or dense) vegetation.
Remote Sens. 2021, 13, x FOR PEER REVIEW 13 of 20 By comparing the reflectivity of our classified vegetation classes in the monitored locality of Malý Šariš, we can state their relatively high spectral proximity (Figure 9).A more significant difference occurs only in the case of tree and shrub (or dense) vegetation.This spectral proximity indicates complications in the automated process of classifying multispectral records, which was confirmed by using the method of supervised classification of the image from summer period.For these purposes, we used the Semi-Automatic Classification Plugin (SCP) extension in the QGIS geographic information system environment (https://fromgistors.blogspot.com/p/semi-automaticclassification-plugin.html,accessed on 23 November 2021).This is a robust extension with a wide range of tools for preprocessing and classifying multispectral records through guided classification.Defining the signatures of individual vegetation classes in the Malý Šariš locality was aimed at distinguishing Solidago spp.from other species of surrounding vegetation.For the Solidago spp.class, data from field measurements were imported as training areas.Other classes (tree and shrub vegetation, mowed and spontaneous grass vegetation) were defined manually in the SCP extension environment.As already indicated by the comparison of the spectral curves of the reflectivity of these classes, their spectral proximity caused complications in the unambiguous spectral definition of This spectral proximity indicates complications in the automated process of classifying multispectral records, which was confirmed by using the method of supervised classification of the image from summer period.For these purposes, we used the Semi-Automatic Classification Plugin (SCP) extension in the QGIS geographic information system environment (https://fromgistors.blogspot.com/p/semi-automatic-classification-plugin.html,accessed on 23 November 2021).This is a robust extension with a wide range of tools for preprocessing and classifying multispectral records through guided classification.Defining the signatures of individual vegetation classes in the Malý Šariš locality was aimed at distinguishing Solidago spp.from other species of surrounding vegetation.For the Solidago spp.class, data from field measurements were imported as training areas.Other classes (tree and shrub vegetation, mowed and spontaneous grass vegetation) were defined manually in the SCP extension environment.As already indicated by the comparison of the spectral curves of the reflectivity of these classes, their spectral proximity caused complications in the unambiguous spectral definition of individual vegetation classes.The spectral ranges of these classes overlapped considerably over the entire spectral resolution of the image, and this shortcoming could not be significantly minimized even by the offered classification algorithms.The resulting classification of the image record (Figure 10) shows considerable inaccuracies when the primarily observed Solidago spp.class is classified to a greater extent than it actually is.Additionally, the class of tree and shrub vegetation is classified in places with significantly dense grass vegetation.On the other hand, sparse grass vegetation is often included in the class of mown grass.The overall accuracy of the classification is very low, reaching less than 28% (Table 4).

Discussion
Management of invasive plant species remains a major challenge globally.Some invasive species have the potential to establish and spread rapidly in new environments, and to alter native plant communities by outcompeting native taxa, thus reducing local biodiversity [40].
There is a need for cost-effective, easy to use, inclusive and repeatable approaches to map invasive alien plants [41].A high-spatial-resolution multispectral image provides rich spectral and textural information [42][43][44].Multispectral and hyperspectral images do not contain three-dimensional structural information, such as canopy height and vertical structures.LiDAR data, a collection of points, are a three-dimensional representation of an object.The features provided by different types of data contain a large amount of key information [11].
In the present study, the four-year monitoring of Solidago spp. in two selected localities was provided.Modeling of the potential influence of climate change predicts increased frequency and severity of heat stress and drought for Central Europe [45,46] The existing research focused on autumn phenological response to summer drought and heat, resulting in indications that the meteorological drought conditioned by lower total precipitation and higher evapotranspiration shifted phenological stages, causing them to occur earlier [47,48].This may have resulted in a shift in timing of the Solidago spp.flowering by about one month during the four-year period of the present research.While, in 2016, we noticed flowering of Solidago spp. in late September (25,27,28), in 2017, flowering occurred in early September (8), and in 2018 and 2019, flowering was noticed in late August (23 and 27).
In connection with previous studies, an increasing number of shoots in the population was also noticed.Plant density in the monitored localities increased minimum three-fold over the four-year period of study.An eight-year study of the relative species Solidago altissima evaluated changes in demographic parameters of shoots over time [49].The proportion of rhizomes producting shoots rose (from c. 50 % to more than 90 %) such that annual shoot density increased almost five-fold by the end of the study, when it had almost reached equilibrium.By the final year of the study, a density of more than 100 shoots per m 2 was reached.Another study focused on the populations of S. gigantea.Comparable field survey between origin country (USA) and introduced country (Europe) was provided.The shoot density of the entire population was estimated from the mean of three randomly placed 1m 2 plots; in very dense populations (more than 100 shoots per m 2 ), the plot size was reduced to 0.25 m.Population varied considerably in their size and density.Population (patch) size varied from 2 to 50,000 m 2 in Europe and from 1.5 to 15,000 m 2 in North America.The mean shoot density was significantly higher in European populations, namely 78.5 ± 4.6 ramets per m 2 compared to 35.6 ± 4.0 ramets per m 2 in American populations [40].
In the present research, state-of-the-art geospatial technologies were used to monitor the spreading of invasive Solidago spp.The aim was to continue monitoring invasive plants population dispersion in selected localities, but mainly to approach the possibilities of using a geographic information system and remote sensing of the Earth in the study of this issue.Sentinel satellite time-series or unmanned aerial vehicle (UAV) data offer the opportunity to characterize vegetation classes at unprecedented spatio-temporal scales.Subsequently, combining data from different types of sensors, such as optical and synthetic-aperture radar (SAR), improves classification accuracy [50].Multispectral sensors on UAVs can map invasive species and their data can be an effective alternative to field data for fitting and validating OCCs (one-class classifiers) based on satellite imagery.However, UAVs can currently cover only a few ha (1 ha equals 10,000 square meters), which reduces their usefulness for local monitoring [51][52][53].
Based on our findings, the multispectral record in four spectral bands (green-550 nm, red-660 nm, red edge-735 nm and near-infrared-790 nm) is insufficient.Vegetation as a whole shows higher values of reflectivity in the summer.On the other hand, a comparison of the spectral curve of Solidago spp. with other herbaceous vegetation shows similar spectral behavior.A transdisciplinary approach combining various remote sensing methods with field trips to provide an accurate and up-to-date understanding of the occurrence and density of invasive alien trees in a cape in South Africa at a 20 m resolution was provided [41].The rich spectral information that they acquired demonstrated that the red edge and shortwave infrared parts of the spectrum were critical for alien tree discrimination, as the key traits distinguishing these alien trees from indigenous shrubland vegetation were differences in biomass and water usage.Discrimination between different types of alien trees and discriminating alien trees from native vegetation with high water use (e.g., wetlands and mountain forests) remains a challenge for the given spatial and spectral resolutions [41].It is more appropriate to use hyperspectral imaging, where the spectral resolution ranges in the order of tens to hundreds of narrowly defined spectral bands, which would significantly refine the feature space Solidago spp., but also the surrounding vegetation, and allow for a clearer distinction of these feature spaces from each other [54].On the other hand, different publications suggest that the use of remote sensing methods with a suitable system configuration (hyperspectral camera) represents an effective way to map the spread of invasive plants.Comparative studies have found that remote sensing variables can perform as well as field measurements [55].It is expected that variations in spectral signatures (shape and the depth of the shape) should be found across environmental gradients or taxonomic lines [56,57].One of the relatives, an invasive goldenrod, Solidago altissima, is abundant in urban environments of Japan and a threat to native biodiversity.There was a study conducted that investigated the ability to upscale the in situ hyperspectral reflectance signature obtained at crown level for landscape scale detection.The detection of S. altissima distribution at landscape scale by direct upscaling crown scale hyperspectral signature was possible using high resolution satellite image with availability of green (~480 nm) and yellow (~600 nm) spectrum bands at fine resolution (~5 m).The detection, however, was influenced by the phenology state of the flowering stages, community size and adjacent plant that possesses similar carotenoid characteristics [58].Earlier study of relative species Solidago altissima was assessed by aerial hyperspectral imagery with high spatial (1.5 m) and spectral (8.9 nm) resolutions for detecting and mapping the early invasion of understory vegetation in moist tall grassland.Generalized linear models (GLMs) were constructed to predict S. altissima occurrence using 1.5 m pixels from hyperspectral data collected during the spring when understory vegetation was directly observable from above.The results suggest that the aerial hyperspectral images obtained during spring before the seasonal development of the grass canopy are useful for the early detection and mapping of S. altissima invading moist tall grassland [59].For the world's worst invasive alien species, Japanese knotweed (Fallopia japonica), a remote sensing method was proposed to detect local occurrences by low-cost digital orthophotos taken in early spring and summer by concurrently exploring its temporal, spectral, and spatial characteristics.Temporal characteristics of the species were quantified by a band ratio calculated from the green and red spectral channels of both images.The combination of all signatures resulted in a producer accuracy of 90.3% and a user accuracy of 98.1% for F. japonica when validation was based on independent regions of interest.A producer accuracy of 61.4% was obtained for F. japonica when comparing the classification result with all occurrences of F. japonica identified during a field validation campaign.The possibility to use orthophotos, which are commonly available and easily accessible in most municipalities, facilitates an immediate implementation of the approach in situations where intervention is urgent [60].The invasive species Lepidium latifolium was monitored by the remote sensing for the more efficient management of this weed in San Francisco (USA).Advanced remote sensing datasets were shown to be sufficient for species distribution modelling.Remote sensing offers powerful tools that deserve wider use in ecological research and management [24].We used the late fall Landsat TM near-infrared band to map the branches of a Eurasian riparian tree Tamarix ramosissima in Colorado, USA when the plants were leafless [61].Despite these successes, in most cases, observing alien plants requires data collected from sensors pushing the limits of at least one type of resolution (spatial, temporal, or spectral resolution) since the profiles of these species may be quite similar to those of native plants from a remote sensing perspective [62].Everitt and his colleagues conducted monitoring studies in the USA using aerial photographs taken during the flowering seasons of Eurasian Euphorbia esula and Asian Tamarix chinensis [63,64].They found that the visible-wavelength (400-700 nm) reflectance of infested settings was significantly higher due to the bright-colored inflorescence.Color infrared (CIR) aerial photographs are also commonly employed, with resolutions ranging from a few centimeters (aerial videography) to ~2 m [65].In the previous study, time-series panchromatic and CIR aerial photographs were used that were acquired during the onset of Heracleum mantegazzianum invasion (the largest central European forb native to the western Caucasus) in Czech Republic [30].The brightness of H. mantegazzianum was much higher than the neighboring vegetation, not only during the flowering, but in the early fruiting season (June-August) due to its distinct structure of fruiting umbels.The timing of data acquisition is crucial for aerial photograph analysis [66].Reviews suggest that alien plant invasions can be studied using remote sensing when the invader presents a novel structure, phenology, or biochemistry relative to neighboring native vegetation.A specific set of remotely sensed data and techniques can be utilized to study each type of invasion characteristic [67].Aerial remote sensing approach would be useful for effective monitoring and management of the species.However, it is challenging for remote sensing techniques to detect a specific single plant species and population size in highly heterogeneous urban and suburban environments [58].

Conclusions
The four-year monitoring observed the density of Solidago spp.and its spreading in two low acreage localities in East Slovakia.The highest density and the highest number of ramets in the northern part of the locality of Malý Šariš were found near tree and shrub vegetation.This may represent a barrier to wind transmission or, worse, the penetration of growing roots, which subsequently accumulates plants in the part before the barrier.In the southern part, which has been affected by human activity, the incidence of Solidago spp.Has already increased.We can, therefore, assume that if there are no significant changes in this location in the coming years, the dispersion of Solidago spp.will increase.In any case, human activity of landscaping of permanent management of the areas (mowing) is needed, which will weaken the growth and vigor of the species under study of the part of the source site in recent months, which may lead to a reduction in seed transfer.However, the studied area of Chminianska Nová Ves has no source locality in its vicinity from which Solidago spp.would expand to a large extent, which is a signicant positive.During the last monitoring year, we noticed a decrease in the number of Solidago spp.The whole locality is covered with dense vegetation, thanks to which the seeds are not able to enter the soil and thus do not have sufficiently suitable conditions for their growth.We could, therefore, say that the original plant species in this case displace invasive plants, which have not found suitable conditions for their intensive spread.During the four-year monitoring, we noticed that the flowering of Solidago spp. is influenced by climate change, which was reflected in the earlier beginning of flowering each year.While, in 2016, flowered plants were recorded in mid-September, in 2019, this flowering was recorded in the second half of August.Through this research and subsequent analysis of the obtained data by state-of-the-art tools, we have proved that geography, with its methods and modern technologies, has enormous possibilities for use in other scientific disciplines.Based on this, we can say that in the coming years, the geographic information system will be an important and sought-after source data.
Spread of Solidago spp.via GNSS at the Research Localities 3.1.1.Mapping the Spread of Solidago spp.via GNSS at the Locality of Malý Šariš Spread of Solidago spp.via GNSS at the research localities 3.1.1.Mapping the Spread of Solidago spp.via GNSS at the Locality of Malý Šariš

Figure 4 .
Figure 4. Solidago spp.density in locality of Malý Šariš in 2017 (a), 2018 (b) and 2019 (c), red arrow shows the direction of the seeds spreading from seed bank.

Figure 4 .
Figure 4. Solidago spp.density in locality of Malý Šariš in 2017 (a), 2018 (b) and 2019 (c), red arrow shows the direction of the seeds spreading from seed bank.

Figure 8 .
Figure 8.The spectral reflectance curve of Solidago spp. at the beginning and culmination of growing season.

Figure 8 .
Figure 8.The spectral reflectance curve of Solidago spp. at the beginning and culmination of growing season.

Figure 9 .
Figure 9.The spectral reflectance curves of classified class of vegetation in the Malý Šariš locality.

Figure 9 .
Figure 9.The spectral reflectance curves of classified class of vegetation in the Malý Šariš locality.

Table 1 .
Dates of monitoring of Solidago spp.In the 4-year period.

Table 2 .
Comparison of Solidago spp.dispersion at a selected localities during mapped period of four years.Analysis of Solidago spp.Changes in Propagation through Spatial Analyzes Using GIS in Research Localities 3.2.1.Analysis of Solidago spp.Changes in Propagation through Spatial Analyzes Using GIS in Malý Šariš

Table 4 .
A confusion matrix of classified multispectral image of the locality of Malý Šariš.

Table 4 .
A confusion matrix of classified multispectral image of the locality of Malý Šariš.