Performance Evaluation of RTK-GNSS with Survey-Grade Receivers and Short Observation Times in Forested Areas

1. Introduction

Since their introduction in the 1980s, Global Navigation Satellite Systems (GNSS) have brought about a paradigm shift in all fields dependent on position [1]. One such area is the management of natural resources, where satellite positioning has become an invaluable tool for administrators [2]. The integration of satellite positioning in forestry operations has been somewhat slower compared to other areas, mainly due to the difficulties of obtaining accurate data under forest canopies [3], as the reduced performance of GNSS receivers in forested environments has been an issue of interest in the last few decades [4,5]. However, the continuous improvement of the technology has alleviated these concerns and led to GNSS becoming pervasive in forest management [6] and an essential component of precision forestry [7]. Of the various uses of satellite positioning in forestry, we can note: the positioning of sample plots for forest inventory to facilitate revisitation [4], mapping of ownership or management unit boundaries [8], the planning and execution of forest roads [7,9], the monitoring of harvesting operations [10,11], the integration of field measurements with remotely-sensed data [12] or the tracking and monitoring of wildlife [13]. Evidently, applications related to these operations are characterized by different requirements in terms of positioning accuracy.

Forest cover will prevent some of the radio signals from reaching the GNSS antenna on the ground and therefore lead to either a degradation of accuracy [5] or even the impossibility of obtaining location data. Therefore, it can be argued that the most important factor leading to the adoption of satellite positioning in forestry research and application is the constantly increasing number of GNSS satellites in operation. Presently, the four GNSS constellations (NAVSTAR GPS, GLONASS, Galileo and BeiDou) amount to more than 100 satellites orbiting the Earth. Other important advancements that took place in the last decades are: (1) the introduction of additional satellite signals that allow dual or triple-frequency receivers to model the ionospheric error [14], (2) the evolution of differential GNSS methods such as Real-Time Kinematic (RTK) that make use of corrections determined over known points (bases) to improve positioning accuracy [15], (3) the expansion of regional, national or even continental RTK networks [16] that allow users access to differential GNNS without the need of setting up their own bases and (4) the implementation of auxiliary systems meant to increase the accuracy and reliability of satellite positioning, such as the Wide Area Augmentation System (WAAS) in North America [17] or the European Geostationary Navigation Overlay Service (EGNOS) [18]. In addition, we must also note the continuous improvement of GNSS receivers, in terms of hardware and software capabilities.

NAVSTAR GPS, the initial GNSS system, was developed with the intention of providing position data in open spaces [19]. Nowadays, especially due to the proliferation of RTK, which is still the most accurate method and can achieve (sub)centimeter-level errors [20], we can consider the issue of accuracy to be solved in the case of open space, for all but the most highly-demanding applications (in terms of positioning accuracy). This is not the case for challenging environments, such as urban canyons, dense forest canopies or steep slopes, all of which have a significant impact on GNSS satellite signals [21] and can therefore cause severe degradations in terms of accuracy and reliability [22]. There are two main issues associated with such environments: (1) a significant reduction of the number of line-of-sight (LOS) satellites due to the presence of obstructions [23] and (2) a higher incidence of multipath, which occurs when radio signals do not reach the receiver antenna directly, but after being reflected by buildings, branches, tree trunks or the ground surface. Unlike most other error sources, the effect of multipath cannot be reduced through differential techniques [24] and is presently considered the most impactful error source [25], especially in challenging environments such as under the forest canopy [26,27].

Given the rapid changes of GNSS and GNSS-adjacent technologies, its potential has to be constantly re-evaluated [6], especially in conditions that cause signal obstructions where there is still a potential for improving accuracy. Because of this, many studies undertaken over the last few decades have focused on satellite positioning in forested terrain. The reported accuracies vary wildly, which is to be expected given the vast differences in terms of: instruments, number of satellite constellations in operation at the time of publication, use of RTK vs. static determinations, length of recording sessions and forest conditions [11]. Therefore, the following literature review will mostly focus on the qualitative findings of previous research and less so on reported numerical values of accuracy indicators. One of the first comprehensive studies of GPS accuracy under the forest canopy focused on the influence of antenna height [4]. The study reported a better signal reception (and therefore improved accuracy) at antenna heights of 4 meters, as compared to heights of 2 meters, and a linear relationship between errors and the Point Dilution of Precision (PDOP). [28] carry out an extensive survey campaign using differential GPS, collecting data over nine months for 27 sites located inside the forest and concluding that forest type is the most important factors influencing accuracy, with mean errors of 8.4-8.5m in hardwood stands and 7.8-11.3m in conifers, dependent on terrain. [5] reports a mean horizontal error of approx. 1 meter for observation times of 10 minutes, showing that canopy density, tree species and observation times influence accuracy. [29] tested the performance of GPS+GLONASS in open vs. forested areas, reporting centimeter-level accuracies but with long observation times (300-350 minutes per point) and highlighting the difficulty of GNSS signal reception under the canopy. [30] used a survey-grade receiver and differential GNSS to position tree stems in various canopy conditions, reporting more accurate positions in conifers and horizontal RMSEs in the range of 14-22 meters. [26] focused on individual tree characteristics, rather than general forest conditions, establishing that the location of nearby trees, rather than their size is the primary factor affecting horizontal accuracy. [6] also report a possible interaction between individual tree locations and accuracy, but without focusing on this matter. A certain impact of forest conditions on GNSS performance is commonly reported, but rarely with statistically significant correlations established [6,9,15,26,31,32,33]. However, in other cases no significant influence of forest composition on GNSS performance is found [3,33,34,35].

Attempts at improving accuracy of positioning in forest conditions have been made, with varying degrees of success. The most common ones involve the use of: dual-frequency receivers [1,9], Differential GNSS [5,15,25] or Precise Point Positioning (PPP) [15], signals from multiple constellations to increase the number of satellites visible [29,36], a system composed of 3 receivers mounted on a metal beam [25,37], different methods of holding the receiver [3], the addition of a metal plate under the receiver to reduce multipath [1] or the fusion of GNSS determinations with data provided by an Inertial Measurement Unit (IMU) [11].

Given the rapid evolution of satellite positioning methods and the lack of consensus with regards to the interaction of environmental variables (especially forest composition) and GNSS accuracy, this study aims to analyze the accuracy and precision of a survey-grade receiver operating in RTK-mode for fast (30 sec) point determinations in two test plots (with deciduous and coniferous vegetation cover, respectively). The short observation times used here are a consequence of focusing on practical uses of GNSS receivers in field data collection campaigns carried out for forestry applications. This also permits a larger sample size, as compared to most other previous research. Both horizontal and vertical accuracy and precision were considered, along with factors such as: forest type, ground slope, canopy density as estimated by the Normalized Difference Vegetation Index (NDVI), DOPs or the number of visible satellites. All tree stems positions were determined using traditional, high-accuracy, surveys with total station equipment, to also account for the number of and distances to trees around GNSS-determined points.

2. Materials and Methods

2.1. Study Area

This research was carried out in a forested area located in the central part of Romania (45⁰ 43’ 30’’ N, 25⁰ 39’ 19’’ E), part of the Natura 2000 site “Lempeș Citadel Hill – Hărman swamp” (Figure 1). The area is almost completely covered by dense forest vegetation and is characterized by hilly terrain, with relatively steep slopes and an elevation range of 500-700 meters. Two test plots were established in this area, one inside a Scots pine (Pinus sylvestris) forest stand and one inside a deciduous stand composed mainly of beech (Fagus sylvatica) and oak (Quercus robur). The pine plot has an area of approx. 11,000 m², while the deciduous one has an area of approx. 25,000 m². General forest characteristics for the management units in which the plots are located, as published in the latest forest management plan, are presented in Table 1.

2.2. Field Data Collection

To estimate the accuracy of GNSS positioning in forest environments, points were distributed over the test plots and marked with wooden stakes. While points were for the most part randomly chosen, care was taken not to position them too closely (< 0.5 meter) to tree trunks, which would hinder the operation of the GNSS receiver. Data was collected using a Stonex S10 survey-grade GNSS receiver, which can receive satellite signals from all four global constellations presently in use (NAVSTAR GPS, GLONASS, Galileo and BeiDou) and has a stated accuracy for RTK operation of 2.5 mm + 0.1 ppm RMS (horizontal) and 3.5 mm + 0.4 ppm RMS (vertical). Determinations were made using the RTK method, with the receiver obtaining corrections directly in the field from a Continuously Operating Reference Station (CORS) of the national RTK network in Romania (ROMPOS – Romanian Positioning System). The receiver was mounted over each point using a monopod with levelling (Figure 2A) and 30 positions (1-second epochs) were determined to obtain average positions. The receiver identified the SFGH reference station as being closest, at approx. 19 kilometers from the study area, and therefore used it for correcting obtained coordinates. A mask angle of 10⁰ was used to filter out signals from satellites close to the horizon, in order to reduce signal interference and the chance of multipath occurring [7]. The antenna was kept at a constant height of 2.50 meters. This ensures that the position of the operator next to the receiver does not influence determinations, while also reducing accidental antenna movements due to wind. It must also be noted that the purpose of the study was not to obtain the most accurate positions possible, but rather to obtain any positions and then determine their accuracy. To achieve this, a high threshold for PDOP of 8 was used, and the receiver recorded data as long as it was in connection with the reference station, even if it did not have a Fixed/Float solution. In other words, GNSS data collection was interrupted only by the occasional, short-term loss of GPRS signal which caused the receiver to enter autonomous mode. WGS84 coordinates specific to GNSS observations were transformed into planimetric coordinates (X, Y, Z), using the official grid coordinate system in Romania (Stereographic 1970) and the TransDatRO software developed by the Romanian National Center of Cartography which is integrated directly in the receiver’s operating system.

The coordinates of points were also determined by traditional ground surveys carried out with total station equipment, to serve as reference (ground truth) against which to compare GNSS positions. For the beech plot, a closed-loop traverse was surveyed, with 5 points inside the canopy and 3 points outside it. These 3 points were determined as references for traverse calculation, using RTK-GNSS and the same SFGH reference station. Each point was determined twice, once in the morning and once in the evening, to ensure significant changes of satellite geometries (Figure 2B). For each observation, a mean position was calculated based on 300 1-sec epochs (about 5 minutes of observation time). In this case, a PDOP threshold of 2 was used and only Fixed solutions were stored. The two pairs of coordinates were compared, found sufficiently close to each other for all points (Table 2), and averaged to obtain final coordinates. In a similar manner, for the pine plot another closed-loop traverse consisting of 7 points was determined, with 3 points located outside the canopy and used as reference. Traverse misclosures were in the range of 2-4 cm (horizontal) and 6-8 cm (vertical). During the topographical survey, the location of tree trunks was also determined, with 1,185 trees identified (432 in the coniferous test plot and 753 in the deciduous plot, respectively).

In total, 89 test points were determined, 65 in the deciduous and 24 in the pine plot. The relatively low number of points in the pine stand is due to the extensive presence of under-canopy vegetation, which hindered both operator movement and visibility from the total station’s location. Field data collection was collected in February-March 2024, so the deciduous stand was in leaf-off condition. Canopy conditions for the coniferous test plot are presented in Figure 3, while Figure 4 shows the conditions for the deciduous test plot.

2.3. Accuracy and Precision Estimates

Repeated GNSS observations (30 per point in this study) provide an estimate of precision, but accuracy must be established by comparison with external data as reference (in this case, coordinates obtained by traditional surveying techniques). Precision is calculated and stored directly by the GNSS receiver for the horizontal and vertical directions and is simply the standard deviation of repeated observations recorded for each point.

Regarding accuracy, this was also considered for the horizontal and vertical component and involves the calculation of horizontal (e_H, which will be also referred to as displacement) and vertical errors (e_Z):

$$e_H=\sqrt{{(X_{GNSS}-X_{Truth})}^2+{(Y_{GNSS}-Y_{Truth})}^2}$$

(1)

$$e_H=\sqrt{{(X_{GNSS}-X_{Truth})}^2+{(Y_{GNSS}-Y_{Truth})}^2}$$

(2)

$$e_Z=Z_{GNSS}-Z_{Truth}$$

(3)

where (X_GNSS, Y_GNSS, Z_GNSS) are the coordinates determined by the receiver and (X_Truth, Y_Truth, Z_Truth) are reference coordinates from the topographical survey.

Once individual errors are known, accuracy can be determined for groups of points, stratified based on forest type (i.e., test plot) or for categories of factors such as the number of Float solutions, the number of satellites and so on. In this case, accuracy is estimated by the Mean Unsigned Error (also called Bias), Mean Absolute Error (MAE) and Root Mean Square Error (RMSE):

$$Bias={\displaystyle \frac{{\sum }_{i=1}^n{e}_i}{n}}$$

(4)

$$MAE={\displaystyle \frac{{\sum }_{i=1}^n\left|e_i\right|}{n}}$$

(5)

$$RMSE=\sqrt{{\displaystyle \frac{{\sum }_{i=1}^n{e}_i^2}{n}}}$$

(6)

where e_i is the error (horizontal or vertical) corresponding to the i-th point out of the n total points included in the respective group. Note that since horizontal errors represent displacements and are always positive, Bias is only reported for vertical errors.

2.4. Potential Factors Considered as Influencing GNSS Accuracy and Precision

Several factors were considered for the analysis of GNSS accuracy in forested areas. These factors can be classified thus:

• GNSS factors: Point Dilution of Precision (PDOP), Horizontal Dilution of Precision (HDOP), Vertical Dilution of Precision (VDOP), number of visible satellites and the proportion of Float/DGPS solutions.
• Environmental factors: ground slope (calculated from a DTM interpolated at a 5-meter resolution from Z_Truth values), forest composition (which test plot the points belong to) and NDVI (calculated in Google Earth Engine, at a 20-meter resolution from cloud-free Sentinel-2 imagery collected between 2022-2024).
• Factors related to tree locations: distance to the nearest tree, number of trees within radii of 2, 4, 6, 8 and 10 meters around each point and the average distance to trees in these radii.

3. Results

3.1. GNSS Accuracy in Coniferous/Deciduous Forest Conditions

Statistical indicators of accuracy were calculated both for all 89 surveyed points and stratified by forest type (Table 3). Accuracy, especially in terms of horizontal displacement, is relatively good given the challenging conditions, with a MAE of 1.63 meters and an RMSE of 2.03 meters. The most accurate horizontal position corresponds to a beech point, which had a horizontal displacement of 0.15 m. Meanwhile, the most inaccurate position is also found in the beech/oak stand and has a displacement of 6.74 meters. It is apparent that horizontal accuracy is improved in the beech/oak stand, with RMSE increasing by 34% (from 1.84 to 2.47 m) and median error increasing by 71% (from 1.21 to 2.07 m) in the pine stand.

In terms of vertical accuracy, while 3 of the 89 points exhibit negative vertical errors (ranging from – 0.02 to – 1.462), it is evident that GNSS observations tended to overestimate elevation (with Bias values of 4.40 and 2.76 m for the beech/oak and pine points, respectively). Also, here we find a reversal of the situation highlighted for horizontal accuracy, with observations in the pine plot having better vertical accuracy than those in beech/oak, with the latter case showing an increase of 51% for RMSE (from 3.49 to 5.27 m) and 36% for the median error (from 2.65 to 3.61 m).

In terms of statistical distribution (Figure 5), while W values for the Shapiro-Wilk normality tests show that the distribution of error values is close in shape to the normal distribution, p-values indicate that the null hypothesis can be rejected, suggesting that errors significantly deviate from a normal distribution in all cases, with the sole exception of vertical errors in the pine stand (p-value = 0.22, n = 24). Given the results of the Shapiro-Wilk test, the two samples (coniferous vs. deciduous GNSS observations) were compared at a 0.05 significance level using the Wilcoxon rank-sum test, which does not assume normality [38]. This showed significant differences between the two plots, in terms of horizontal (W = 562, p-value = 0.044) and vertical accuracy (W = 1038, p-value = 0.017).

In terms of the direction of displacement, a tendency of points in the pine plot to be shifted towards the West and South is found, while points in the beech/oak plot tended to be shifted towards the East and West/South-West (Figure 6).

3.2. Variation of GNSS Conditions under the Forest Canopy

As previously discussed, several factors related to GNSS data collection and recorded by the receiver were considered: PDOP, HDOP, VDOP, the number of visible satellites and the type of GNSS solution (Table 4). Note that no Fixed solutions were obtained, with the 30 GNSS determinations for each point being split between Float and DGPS solutions.

DOPs indicate consistently appropriate satellite configurations in both test plots, with DOP values under 4 for all points (Figure 7). This is due to the high number of visible satellites, which had an average of approx. 13 and a minimum of 8 satellites, in both forest stands. The similarity of GNSS conditions for the 2 test plots is also evidenced by the Wilcoxon rank sum test, with p-values between 0.074 and 0.303 for the 5 factors considered. In terms of solutions, for the most part only DGPS solutions were achievable under the dense canopy cover, with 83% of points having 30 out of 30 DGPS solutions and only 3 with at least 10 out 30 Float solutions (1 point in pine and 2 in beech/oak, respectively).

3.3. Distribution of Trees Around Test Points

Distances between points and their nearest tree range from 0.20 to 9.44 meters, with an average of 2.76 meters. Note that, while point positions were chosen so that no tree was closer than 0.5 meters, this was estimated visually in the field, so 4 of the points (2 in each test plot) ended up closer than 0.5 meters to a tree. In terms of the number of trees located around points inside various radii (from 2 to 10 meters, in increments of 2 meters), we must note that the 2-meter radius was excluded from further analysis, as 75% of observations had 0-1 trees inside that radius (Table 5).

3.4. Importance of Factor Variation on GNSS Accuracy

External factors were tested in terms of their relative importance in explaining the variation of horizontal/vertical accuracy and precision. For brevity, only the 4 most important factors are presented (Table 6).

The factors considered in this study do a rather poor job of estimating overall accuracy (25-27% of variation explained), but this is mostly due to deciduous test plot. When only points inside the pine stand are considered, we find that the factors involved explain more than 80% of the variation in horizontal/vertical accuracy. On the other hand, GNSS precision (i.e., the standard deviation of the 30 positions determined for each point) is better correlated with the factors considered, with almost half of variation explained overall and nearly all variation (86-98%) explained for the coniferous test plot. Some additional findings are also worth highlighting. First, DOPs only influenced precision, having virtually no influence on accuracy. Second, ground slope appears to only influence vertical accuracy. Third, factors related to tree location (distance to nearest tree, number of and average distances to trees in various radii) only have somewhat of an influence in the coniferous plot, but overall do not seem to be deciding factors on accuracy or precision.

These findings are supported in part by the results of Analysis of Variance (ANOVA) tests, which identify several factors as having a significant influence on GNSS accuracy (p-values < 0.05):

• When all points are considered together: Species is significant for both horizontal and vertical accuracy, while No. of satellites and Slope are significant only for vertical accuracy.
• For the coniferous plot: No. of trees in a 4m radius and Distance to nearest tree are significant for horizontal accuracy, while no factors are identified as significant for vertical accuracy.
• For the deciduous plot: no factors are identified as significant, for both horizontal and vertical accuracy.

Regarding precision of GNSS observations, ANOVA showed that:

• When all points are considered together: PDOP and No. of satellites are significant for horizontal/vertical precision, while Species and Slope are significant only for horizontal and vertical precision, respectively.
• For the coniferous plot: No. of satellites is significant for horizontal precision, with no significant factors for vertical precision.
• For the deciduous plot: PDOP is significant for horizontal precision, while Slope is significant for vertical precision.

3.5. Level of Agreement between GNSS Accuracy and Precision

The precision of GNSS-determined positions is commonly expressed by the standard deviation of repeated observations of the same point. Therefore, standard deviations of GNSS coordinates have been compared with the horizontal/vertical accuracy (i.e., the displacement of point locations and vertical errors, respectively). Results show little to no correlation between accuracy and precision, especially for the horizontal component (Figure 8). The overall r-value for the horizontal component is 0.09, while for the vertical component an r-value of 0.35 was found. Meanwhile, in pine r-values of 0.06 (horizontal accuracy/precision) and -0.11 (vertical accuracy/precision) are found. On the other hand, accuracy is somewhat better correlated with precision in the beech/oak plot, with r-values of 0.21 (horizontal) and 0.46 (vertical).

4. Discussion

4.1. GNSS Accuracy in Coniferous/Deciduous Forest Conditions

As previously discussed, in the last few decades the term GNSS has come to cover a wide range of instruments, methods and technologies. This makes quantitative comparisons of accuracy between studies difficult and not particularly relevant, especially when also considering differences in terms of experimental design. As an example, [29] obtained centimer-level accuracy using RTK under the forest canopy, but only by allowing the receiver to resolve ambiguities, which led to much longer observation times than in our study (up to 6 hours per point). On the other hand, [30] also used RTK-GNSS in forest conditions and obtained significantly worse accuracies than those reported here (RMSE values of 16-20 meters), but this could be explained by the fact that the GNSS receiver was positioned very close to tree trunks (30 cm) in order to determine their position, which was not the case for this study.

However, by considering the qualitative (rather than quantitative) findings, our results can be put into a larger context. First, while the influence of canopy presence itself was not a focus of this study, it is well established that RTK-GNSS can easily achieve centimeter-level accuracy in open ground, so our findings (RMSE of 2.03 and MAE of 1.63 meters for horizontal accuracy) clearly indicate a significant degradation of GNSS performance in forest environments. This is in line with previous studies [3,9,34,37,39]. As far as vertical errors are concerned, results show a decrease in accuracy of approx. 2x. Similar proportions has been previously reported [29,36] and expecting a two-fold degradation of accuracy when using GNSS to determine altitude instead of location has been a sort of rule-of-thumb for some time [40]. Secondly, differences between errors depending on forest compositions were found to be statistically significant, which is in line with some [6,31], but not all previous studies [3,34]. Note that forest composition (for e.g., pine vs. beech/oak, in this case) is a necessary simplification that does not account for stem density, differences of age among trees or the presence and type of sub-canopy vegetation, among other conditions that might potentially influence GNSS performance. This could explain the seemingly contradictory findings reported previously. The fact that horizontal accuracy was better in deciduous forest is explained by the leaf-off conditions in which GNSS data was collected. However, the relatively better vertical accuracy in the pine plot is not easily explained. This behavior could potentially be attributed to the small sample size for pine (n = 24), the difference in slope (pine avg. slope: 19%; beech/oak avg. slope: 24%) and aspect (the pine plot is mostly oriented to the SW, while the beech/oak is oriented towards the NW) or most likely a of combination of these factors. Without additional data these hypotheses cannot be tested.

4.2. GNSS Solution Type under the Forest Canopy

GNSS solutions can be of three types, which in order of increasing accuracy are: DGPS (code-based positioning), Float (the receiver has locked onto satellite signals for enough time to allow for more accurate carrier-phase positioning) and Fixed (when enough time has passed so that the receiver can resolve the integer ambiguity specific to carrier-phase measurements). In this study, no Fixed solutions were achieved, with almost all epochs recorded as DGPS solutions (97% of 2,670 solutions – 30 per point). [37] also report achieving Fixed solutions only for a few percent of observations. This potentially goes towards explaining the relatively low accuracy, especially given the high number of visible satellites (13 on average) and good PDOP values (average PDOP is under 2). Time-to-fix is outside the scope of this study so was not comprehensively tested, but it should be noted that, at times, the GNSS receiver was left installed on the monopod under the canopy for as long as 1-2 hours and no Fixed solution was achieved (only Float, given sufficient time). The difficulty of carrier-phase determinations under the canopy due to signal interruption hindering integer ambiguity resolution, has been previously established [7,9,32,33]. GNSS users operating in forest ecosystems should also take into account that, even if a Fixed solution is achieved, the integer resolution could be wrong, which is more likely to happen in challenging conditions [9]. This could lead to the unexpected situation in which Fixed determinations are affected by large errors [23] and therefore less accurate than Float/DGPS ones.

4.3. Relative Importance of Factors for GNSS Accuracy/Precision

In terms of external factors considered, our main finding is that, as far as accuracy is concerned, environmental conditions (such as tree stand species, NDVI or ground slope) tend to explain more in terms of GNSS performance than GNSS-specific indicators such as the number of visible satellites or DOPs. This is to be expected, as the accuracy of satellite positioning is a function of the environment in which determinations are carried out [37]. [19] also reported that forest structure parameters have a higher influence on satellite positioning than factors related to GNSS, while [33] showed that the main driver of both accuracy and precision of GNSS are the canopy characteristics. Overall, the factors considered here explain most of the errors found in coniferous conditions (over 85% of variance explained, for both horizontal and vertical accuracy) but this does not apply to the deciduous plot, where only 17-29% of error variance is explainable. As far as precision is concerned, GNSS-specific factors (especially PDOP, HDOP and VDOP) explain significantly more of the variation in performance for both test plots. The only case in which precision is poorly explained is for vertical precision in the deciduous plot (29% of variance explained). This is another indicator that the unexpectedly large vertical errors for this plot (significantly higher than the ones in the coniferous plot) are due to unconsidered factors, warranting further study.

Factors related to the distribution of trees around the GNSS receiver rarely showed predictive power, especially for the case of horizontal accuracy/precision. This would indicate that, rather than the tree trunks, factors related to the tree crowns (size, shape and so on) which were not considered for this study would likely have more of an impact on GNSS performance. [26] report some predictive power for the location of nearby trees but establish that the size of trunks has no correlation to errors, which we believe supports this hypothesis.

4.4. On Accuracy Versus Precision of GNSS-Determined Positions

The issue of accuracy/precision is an important one, as in practice users would rarely have access to reference, ground-truth data to estimate the accuracy of GNSS positioning. Rather, precision has to be relied upon as a substitute for accuracy. While it is expected that precision underestimates accuracy, given that accuracy is a summation of precision and systematic and gross errors [41], it is desirable that the two are correlated. This would at least imply that precision estimates (given for example by the standard deviation of coordinate values for repeated determinations of the same points) allow comparisons of the relative accuracy achieved. However, our results indicate that precision estimates have very little in common with accuracy, especially for latitude/longitude determinations. Standard deviations of computed coordinates have been previously reported as insufficient for estimating GNSS accuracy [31,42,43]. On a related note, Dilution of Precision (DOP) indicators were found to be better correlated with precision, rather than accuracy. As an example, r increases from 0.09 (correlation of PDOP with horizontal errors) to 0.45 (correlation of PDOP with horizontal precision). Previous studies also found PDOP to have little or no predictive power for accuracy [7,42]. With the increase of GNSS satellites in orbit, improper satellite geometries would only be found in the most challenging environments, meaning that nowadays the utility of DOPs in explaining GNSS performance is questionable.

5. Conclusions

GNSS benefits from constant technological advances and breakthroughs, therefore testing its performance is an ongoing concern [44]. The aim of this study was to evaluate the performance of GNSS positioning using a survey-grade receiver, for the challenging case of relatively steep terrain covered by dense forest vegetation. The experimental design (short observation times of 30 seconds/point using RTK with corrections provided by a CORS of Romania’s national system of GNSS reference stations) aimed to mimic, at least in part, practical usage of GNSS by operators in the forestry sector. Results show that, using this methodology, meter or even sub-meter accuracy is achievable, even with code-based determinations due to the satellite signal interference due to the presence of vegetation. Horizontal/vertical errors show statistically significant differences between coniferous and deciduous stands. However, more forest structural parameters need to be considered to better explain the variation of GNSS performance. Our analysis highlighted that, under dense forest canopies, environmental conditions might affect accuracy to a larger extent than GNSS-specific factors such as the number of satellites and PDOP values. Presently, the most critical issue for satellite positioning in forest areas would be the mitigation of multipath, which is an active area of research [24].

Overall, using survey-grade receivers and RTK-GNSS with corrections from national/regional networks allows rapid positioning under the canopy at adequate accuracies (especially if altitude data is not required), depending of course on application requirements. Further advancements of GNSS are expected, especially in terms of antenna design, instruments or software processing, which imply that a continuous reassessment of the technology’s potential is warranted [3]. The added benefit of increasing observation time or the potential of modern, RTK-enabled handheld receivers (such as those integrated into smartphones or even smartwatches) to achieve a similar performance in forest conditions to that of more cumbersome survey-grade receivers, are issues to be addressed in further studies.