Real-time positioning in logging: Effects of forest stand characteristics, topography, and line-of-sight obstructions on GNSS-RF transponder accuracy and radio signal propagation

Eloise G. Zimbelman; Robert F. Keefe

doi:10.1371/journal.pone.0191017

Abstract

Real-time positioning on mobile devices using global navigation satellite system (GNSS) technology paired with radio frequency (RF) transmission (GNSS-RF) may help to improve safety on logging operations by increasing situational awareness. However, GNSS positional accuracy for ground workers in motion may be reduced by multipath error, satellite signal obstruction, or other factors. Radio propagation of GNSS locations may also be impacted due to line-of-sight (LOS) obstruction in remote, forested areas. The objective of this study was to characterize the effects of forest stand characteristics, topography, and other LOS obstructions on the GNSS accuracy and radio signal propagation quality of multiple Raveon Atlas PT GNSS-RF transponders functioning as a network in a range of forest conditions. Because most previous research with GNSS in forestry has focused on stationary units, we chose to analyze units in motion by evaluating the time-to-signal accuracy of geofence crossings in 21 randomly-selected stands on the University of Idaho Experimental Forest. Specifically, we studied the effects of forest stand characteristics, topography, and LOS obstructions on (1) the odds of missed GNSS-RF signals, (2) the root mean squared error (RMSE) of Atlas PTs, and (3) the time-to-signal accuracy of safety geofence crossings in forested environments. Mixed-effects models used to analyze the data showed that stand characteristics, topography, and obstructions in the LOS affected the odds of missed radio signals while stand variables alone affected RMSE. Both stand characteristics and topography affected the accuracy of geofence alerts.

Citation: Zimbelman EG, Keefe RF (2018) Real-time positioning in logging: Effects of forest stand characteristics, topography, and line-of-sight obstructions on GNSS-RF transponder accuracy and radio signal propagation. PLoS ONE 13(1): e0191017. https://doi.org/10.1371/journal.pone.0191017

Editor: Budiman Minasny, The University of Sydney, AUSTRALIA

Received: August 17, 2017; Accepted: December 27, 2017; Published: January 11, 2018

Copyright: © 2018 Zimbelman, Keefe. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: All relevant data are within the paper and its Supporting Information files.

Funding: This work was funded by the U.S. Centers for Disease Control and Prevention (CDC) National Institute for Occupational Safety and Health (NIOSH) (https://www.cdc.gov/niosh/index.htm) grant number 5 U01 OH010841 at the University of Idaho (RFK). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Since the Department of Defense launched its first Navigation System with Timing and Ranging (NAVSTAR) satellite in 1978, global positioning system (GPS) technology has become an integral component of national defense, homeland security, civilian life, and scientific research [1]. Global navigation satellite system (GNSS) positioning is a more general term that encompasses all global satellite-based positioning systems such as GPS, GLONASS, Galileo, and BeiDou. Initial research evaluating the practicality of GNSS technology in forested landscapes indicated the potential use of GNSS for a range of operational and research uses in natural resources [2]. Early assessments of GNSS for forestry found it could be used to locate forest inventory plots [3,4], quickly determine timber harvest unit boundaries [4], locate forest roads [4], and track wheeled skidders [5]. Various other studies tracked harvest equipment with GNSS and used this information to assess site disturbance and to calculate productivity from time-study data [6,7]. More recent studies related to forest operations have analyzed the potential for GNSS data to quantify production efficiency [8], track log trucks [9], improve operational monitoring [10], and increase efficiency and calibrate remotely-sensed inventory data using GNSS-derived harvester head positions [11].

The accuracy of GNSS coordinate readings is dependent upon the number and geometry of satellites visible to a GNSS unit at any point in time. Positional dilution of precision (PDOP) is an index of the influence of satellite geometry on GNSS measurements [12]. In general, a lower PDOP value indicates an arrangement of satellites providing higher measurement reliability and values less than 2 are desirable [12]. PDOP values can be calculated for user-defined locations using GNSS mission planning software [13]. Another factor that could potentially affect GNSS accuracy in forested conditions is aspect. The Wide Area Augmentation System (WAAS) provides real-time GNSS data correction [14]. GNSS receivers need a clear view of a geostationary communications satellite (GEO) in order to receive WAAS correction signals [14]. GNSS receivers in the northern United States usually need an unobstructed view to the south of less than 20 degrees to receive WAAS signals because GEO satellites are low on the horizon [14]. This leads to a potential for increased error on slopes without a clear view to the south [14]. Few studies have evaluated the effect of aspect on GNSS performance, and while two studies found higher GNSS fix rates and lower location error on south aspects, these differences were not statistically significant [15,16].

GNSS use in forestry is often affected by error associated with satellite signal obstruction by the canopy or other solid objects and the reflection or diffraction of satellite signals from nearby objects or surfaces, an error known as multipathing [13,17]. Previous studies have shown that forest stand structural characteristics and terrain affect GNSS accuracy [5,18–23]. Holden et al. [21] developed a method to model GNSS precision using three canopy descriptor variables (percentage of sky obstruction, maximum canopy hole radius, and fragmentation of sky view). Lewis et al. [22] modeled the proportion of 3D GNSS fixes, PDOP, and location error using the percent canopy cover and satellite view (to represent terrain obstruction). Newer studies have taken advantage of the correlation between GNSS signal strength and forest stand characteristics by evaluating the potential to predict and map forest parameters using GNSS signals [24,25]. GNSS receiver type (survey-, mapping-, or recreation-grade) also affects the accuracy of position measurements [5,26–30]. Survey-grade receivers are capable of sub-centimeter accuracy in the open and sub-meter accuracy under mature forest conditions [30,31]. At the other end of the spectrum, recreation-grade units are the least expensive and have accuracies ranging from 2–5 m in the open [13,32–34] to 3.8–12 m in mature forests [13,30,32–34].

Advances in positioning technology for remote environments have emerged simultaneously for several uses, ranging from recreation to public safety and defense. These devices link GNSS positional information with radio frequency (RF) transmission of location coordinates (GNSS-RF) to form ad-hoc networks in which the locations of all units can be monitored on mobile phones or tablets. Like traditional GNSS devices, each GNSS-RF transponder determines its coordinates using one or more satellite-based positioning systems. However, the RF transmission is a second component that allows those coordinates to be sent to other, nearby GNSS-RF units at user-defined intervals. GNSS-RF transponders include consumer-grade units for recreational use like the Garmin Rino and Garmin Alpha 100, devices such as the Raveon Atlas PT marketed for public safety, a variety of military-grade GNSS radios designed for defense applications, and consumer-grade mobile-based solutions from goTenna and Beartooth that turn smartphones into two-way radios for voice or text communication in areas without cellular service. While traditional GNSS devices allow users to see their own positions, GNSS-RF devices enable real-time positioning through location sharing among individuals and equipment in remote locations [8,23,35–37]. Thus, either the device’s native screen or an attached tablet can display the location of other devices in the network moving in real-time.

Many GNSS-RF transponders support geofencing, in which a virtual boundary is defined around a user-defined geographic zone. Geofences are either circular or polygonal in shape, can vary widely in size depending on intended application, and can be stationary or mobile. Alert notifications are triggered as tracked mobile objects cross into or out of the geofence, and this functionality may be useful for a range of operational forestry applications such as detecting the amount of time workers spend near cable logging hazards [23], signaling when log skidders or log trucks cross harvest unit boundaries [36], and delineating tree falling hazard zones around manual fallers [37]. GNSS-RF real-time positioning and geofences have the potential to improve communication and situational awareness on logging operations, in wildland firefighting, transportation, and recreation. In this study, our focus was on logging safety, as logging is consistently ranked as one of the most hazardous occupations in the United States [38], with the highest fatal work injury rate of civilian occupations in 2015 [39]. By increasing situational awareness, the active display of real-time positioning logistics may be able to reduce hazards posed by the frequent interactions among ground workers, heavy equipment, and irregular terrain that are common on active logging operations.

GNSS-RF and related technologies pose new challenges for quantifying positional accuracy because positional error is associated with both the accuracy of GNSS locations and successful propagation of radio signals between devices. In other words, the accuracy of GNSS-RF position sharing depends not only on factors that influence GNSS accuracy as described above, but also on factors that affect radio signal propagation and attenuation. Radio signals experience diffraction, scattering and reflection as they travel through vegetation [40–42]. Leaf state and vegetation depth and density influence radio signal attenuation [42–45]. Attenuation can also be affected by wind [42,43,46], humidity [47,48], rain [46], and terrain [49,50]. It is unclear whether the same factors affect both the GNSS and RF components of emerging technologies, or if different forest stand characteristics and topographic factors affect one or the other.

The purpose of this study was to characterize the factors affecting real-time positioning on irregular, forested terrain through analysis of the effects of forest stand characteristics, topography and other line-of-sight (LOS) obstructions on the GNSS accuracy and radio signal propagation quality of multiple Raveon Atlas PT GNSS-RF transponders (Raveon Technologies, Vista, California, USA) functioning as a network. Until recently, most previous research evaluating GNSS for forestry applications has focused on stationary units. Recent studies characterizing GNSS units in motion include work by Kaartinen et al. [11] and Liu et al. [24], while others [8,23,36,37] have quantified mobile GNSS units paired with RF-based transmission. To further our understanding of location sharing among mobile GNSS units, we evaluated the time-to-signal accuracy of geofence crossing alerts in 21 randomly-selected stands on the University of Idaho Experimental Forest, under a wide spectrum of stand and topographic conditions. Time-to-signal accuracy refers to the difference between when a person or object crosses a geofence and when the alert is generated and shared with other, nearby GNSS-RF devices. Specifically, we tested three hypotheses. Our first null hypothesis was that neither forest stand characteristics, topography, nor the presence of obstructions in the line-of-sight affected the probability of successful radio signal propagation between GNSS-RF units. We tested this by analyzing the odds of missed radio signals within each stand using mixed-effects logistic regression. Our second null hypothesis was that neither stand characteristics, topography, nor physical obstructions affected the stationary positional accuracy of GNSS-RF units. To test this, we used linear mixed-effects models to determine which factors most affected root mean squared error (RMSE) of the PT. Our third null hypothesis was that neither stand characteristics, topography, nor physical obstructions affected the time-to-signal accuracy of geofence crossings. We tested this by using linear mixed-effects models to determine which factors most affected geofence intersection alert delay.

Materials and methods

Field experiment

Five Raveon Atlas PT GNSS-RF transponders collected positional data during the field experiment. As GNSS-RF units, the PTs receive their coordinates and then transmit that information to other PTs using radio frequency. The units can be attached to tablets or computers, which allows ground workers and equipment operators on logging operations to see all other positioning devices in real-time. PTs receive their coordinates from NAVSTAR GPS satellites only and have a specified 24-hour static accuracy of < 2.5 m for 50% of measurements and of < 5 m for 90% of measurements [51]. Depending on terrain, the devices can communicate up to 48 km away and position updates can be transmitted as frequently as one signal per second [51]. PTs can be used with Raveon RavTrack software, which has several options for geofencing, including different notification options.

In this study, real-time geofence alert signals were evaluated in a random sample of 21 stands on the University of Idaho Experimental Forest (UIEF) (Fig 1A). Only stands ≥ 2.02 hectares (5 acres) in size were selected. Within each stand, the timing of geofence alerts was characterized for a manual faller entering a 100 m × 300 m rectangular geofence (Fig 1B). In addition to placing one stationary PT at the geofence intersection point (Atlas PT X in Fig 1B) to record data, a compass and 100-m fiberglass tape were used to place three other PTs 100 m from the virtual boundary intersection point at angles forming the vertices of an equilateral triangle (Atlas PTs A, B, and C at triangle points A, B, and C, respectively, in Fig 1B). The first of these stationary PTs (PT A) was placed at a randomly selected bearing from the intersection point (triangle point A). The bearing was sampled from the set of whole numbers between 1 and 360, with replacement. The remaining two stationary PTs (PTs B and C) were placed 120° and 240° clockwise, respectively, from this first PT (triangle points B and C, respectively). All stationary PTs were zip-tied to wooden stakes such that the bottoms of their antennas were 1 m above the ground surface. Each stationary PT was attached to a Windows 10 Dell Venue 8 Pro tablet running Raveon RavTrack software. Finally, maps for each stand were loaded onto the tablets using 1-m resolution National Agriculture Imagery Program (NAIP) images downloaded from The National Map website [52] and all tablets were synced with the National Institute of Standards and Technology (NIST) time server [53] each day.

Download:

Fig 1. Map of the stand locations on the UIEF and illustration of the experimental setup.

(A) The 21 stands are delineated according to total basal area (m²/ha) and the UIEF unit boundaries are shown in blue. Background map is 1-m NAIP imagery. (B) Illustration of global navigation satellite system (GNSS) technology paired with radio frequency (RF) transmission (GNSS-RF). GNSS-RF transponders (Atlas PTs) receive positional information from GNSS satellites and relay this information to one another using radio frequency transmission. Atlas PT X is located at the geofence intersection point, while Atlas PTs A, B, and C are located at the triangle points A, B, and C, respectively. The manual faller carried a PT attached at the hip (Atlas PT F).

https://doi.org/10.1371/journal.pone.0191017.g001

The orientation of geofence crossing in each stand was also randomly selected from the sample of whole numbers between 1 and 360, with replacement. A rectangular geofence was established in each stand using a Suunto sight-through azimuth compass, fiberglass tape, and an Eos Arrow 100 GNSS unit (Eos Positioning Systems, Terrebonne, Quebec, Canada) with a specified accuracy of < 0.6 m [54]. One side of the geofence was centered at the geofence intersection point and was oriented perpendicular to the crossing direction. The geofence was 100 m wide at the crossing point and 300 m long.

In each stand, a manual faller carrying a PT attached at the belt crossed the geofence once by walking a 90-m route oriented perpendicular to the geofence (i.e., in the chosen geofence crossing direction), starting and ending 45 m from the intersection point. For consistency, the manual faller walked at a pace of 45 bpm, as dictated by a digital metronome. The route was established using a compass and 50-m fiberglass tape and was marked with pin flags. The observed time at which the faller crossed the geofence was recorded in the field using a custom script in R [55] running on a Windows 10 Dell Venue 8 Pro tablet synced with the NIST time server. The predicted intersection times were recorded by the tablets attached to each stationary PT (PTs A, B, C, and X). All PTs were set to collect and transmit their coordinates at a rate of once per second.

Within each stand, topographic, physical, and vegetative obstructions present along each LOS path between the geofence intersection point and the PTs located at triangle points A, B, and C were quantified during setup using a modification of the FIREMON line intercept method [56]. To mark the LOS path, a 100-m fiberglass tape was attached to two metal stakes. One stake was located at the intersection point and the second was located where the stationary PT would be placed during the experiment. The tape was secured to these two stakes and stretched taught 1 m above the ground to match the location of the stationary PT antenna height. This height was approximately equal to that of a PT when worn by a ground worker on a belt clip. The start of the tape was attached to the stake at the intersection point while the end was attached to the stake at the triangle point. For each LOS path, three 5-m sections were randomly selected from the segments shown in Fig 2 and all vegetative obstructions in these sections were classified using the key in Table 1. Only obstacles physically touching the fiberglass tape were quantified. To measure each obstacle, two meter sticks were used to hold the tape 1 m above the ground. The first meter stick was located 1 m before the start of the obstacle and the second meter stick was located 1 m past the end of the obstacle. Once the fiberglass tape was in position, the two locations at which the obstacle first and last contacted the tape were recorded to the nearest centimeter. Obstacles less than 5 cm in size (as measured along the LOS path) were not quantified. When gaps were present between nearby obstacles, the obstacles were treated as two separate obstacles only when the gap was greater than 25 cm. Obstructions less than 1 m tall were not quantified and trees that were less than 12.5-cm diameter at breast height (DBH) were classified as coniferous vegetation. To simplify analysis, vegetative obstructions within the three measured 5-m sections were summarized using Eq (1): (1) Where V_i is the measured distance of all vegetative obstructions recorded for the three 5-m sections along the ith LOS path. S_i, C_i, T_i, W_i, SC_i, SW_i, and CW_i represent the distance of vegetative obstructions defined in Table 1 as measured in the three 5-m sections along the ith LOS path. Then, to account for the fact that only 15 m of each 100-m LOS path was measured, the total distance of vegetation along each path was calculated using Eq (2): (2) Where TV_i represents the total distance of vegetation along the entire ith LOS path and V_i is the measured distance of all vegetative obstructions recorded in the three 5-m sections along the ith LOS path (Eq (1)). Lastly, all boulders, streams, and forest roads were recorded as present or absent along each LOS path, regardless of their location on the path. These were recorded because of the effects they may have as terrain changes. However, because only one boulder was measured in the LOS paths, we removed it from analysis.

Download:

Fig 2. LOS path sections.

Each LOS path was divided into 20 5-m sections and three sections were randomly selected for each LOS path. This figure shows the 20 sections and their locations along the LOS path. Sections highlighted in green represent the three randomly selected sections for which all vegetative obstructions were measured using the key in Table 1.

https://doi.org/10.1371/journal.pone.0191017.g002

Download:

Table 1. Obstacle list and key.

https://doi.org/10.1371/journal.pone.0191017.t001

Once all obstructions were quantified and recorded, each LOS path was walked carrying a Garmin Alpha 100 GNSS-RF unit (Garmin, Olathe, Kansas, USA) to record the vertical elevation profile. Using the Garmin data, each LOS path was classified in terms of the presence or absence of concavity and convexity. A LOS path was concave if the minimum elevation along the path was at least 3 m below the lower of the two path endpoints. A LOS path was convex if the maximum elevation along the path was at least 3 m above the higher of the two path endpoints. The classification criteria for concavity and convexity is illustrated in Fig 3. The percent slope and aspect were also measured at the geofence intersection point. Aspect was measured as a continuous circular variable, but was reclassified as either N (316°– 45°), E (46°– 135°), S (136°– 225°), or W (226°– 315°). Trimble’s GNSS mission planning website [57] was used to determine the predicted PDOP values for each day during the experiment. Sampling was only conducted during times with predicted PDOP values less than 4.5 to ensure consistency.

Download:

Fig 3. Illustration of slope classification as concave, convex, or both.

Blue dots represent the higher of the two LOS endpoints while red dots represent the lower of the two LOS endpoints. Green lines represent the ground surface along the LOS path.

https://doi.org/10.1371/journal.pone.0191017.g003

To quantify forest stand characteristics, a 0.03-hectare fixed-area plot was established in each stand, centered at the intersection point. The DBH, total height, and height to the base of the live crown were quantified for all trees ≥ 12.5-cm DBH within the plot. Using these measurements, total basal area (TBA), trees per hectare (TPH), mean height (Ht), and quadratic mean diameter (QMD) were calculated for each stand and used as variables representing forest stand characteristics during analysis. Quadratic mean diameter is a commonly used metric in forestry that refers to the diameter of the tree of mean basal area (stem cross-sectional surface area), as measured at breast height (1.37 m). QMD is calculated as the square root of the squared stem diameters divided by the number of stems sampled, as defined in Eq (3): (3) Where D_i is the DBH of the ith tree and n is the number of trees sampled.

To quantify missed radio signals in each stand, the number of missed position updates transmitted from the faller’s PT (PT F) to the three PTs at the triangle points (PTs A, B, and C) was calculated for the 90-s interval centered on the observed geofence intersection time. Because all units were set to transmit their coordinates at 1-s intervals, 90 position updates would have been received in this time period in the absence of missed signals.

Stationary GNSS accuracy was summarized using RMSE, which is a common measure of GNSS positional error and represents the difference between the predicted and observed coordinates of a GNSS unit. In each stand, the predicted coordinates were obtained using the stationary PT located at the geofence crossing point (denoted as Atlas PT X in Fig 1B). These coordinates were collected once per second for 5 min prior to the time of geofence crossing. The observed (true) coordinates for each of these PTs were obtained using the Eos Arrow 100 GNSS unit described above. All coordinates were converted to the Universal Transverse Mercator (UTM) projection and then the RMSE for each stationary PT was calculated using Eq (4): (4) Where RMSE_i is the RMSE value in the ith stand, x_i is the observed easting value in the ith stand (i.e., the Arrow 100 easting coordinates), is the jth predicted easting value in the ith stand (i.e., the PT easting coordinates), y_i is the observed northing value in the ith stand (i.e., the Arrow 100 northing coordinates), ŷ_ij is the jth predicted northing value in the ith stand (i.e., the PT northing coordinates), and n is the total number of PT signals received in the ith stand.

The overall geofence intersection alert delay for each stand was derived by averaging the time-to-signal delay calculated at each of the three triangle points (A, B, and C) (Eq (5)): (5) Where D_i is the overall delay for the geofence intersection in the ith stand and O_i is the observed time at which the faller crossed the geofence in the ith stand, as recorded in the field. P_ij is the predicted intersection time in the ith stand at the jth triangle point, represented by the recorded alert at triangle point A. P_ik is the predicted intersection time in the ith stand at the kth triangle point, represented by the recorded alert at triangle point B. P_il is the predicted intersection time in the ith stand at the lth triangle point, represented by the recorded alert at triangle point C. Using this formula, positive delays indicate geofence crossing alerts that were triggered after the faller crossed the geofence and negative delays indicate geofence crossing alerts that occurred before the faller intersected the geofence. We used time-to-signal delay as an integrated measure of the accuracy of mobile GNSS units sharing their locations, which differs from RMSE calculated by Kaartinen et al. [11] using known reference points along a path.

Analysis of missed radio signals

To test the null hypothesis that the probability of successful GNSS-RF signal propagation was not related to forest stand characteristics, topography, or obstructions in the line-of-sight, a binomial generalized linear mixed-effects model was used to evaluate relationships between the odds of missed signals as a function of vegetative LOS obstructions, topography, and forest stand characteristics. The model was fitted using the glmer function in the R lme4 package [58]. Variables included as fixed effects were the total distance of vegetation along each LOS path (TV_i), TBA, TPH, Ht, QMD, slope, aspect, and the presence or absence of forest roads, streams, convex slopes, and concave slopes (Table 2). These variables were included because of their potential effect on the successful propagation of radio signals. To avoid errors with model convergence, the TPH variable was multiplied by a scalar of 0.01. The stand was used as a random effect to account for unobserved variation between stands. The response was the log odds of missed position updates along each LOS path during the 90-s interval centered around the observed geofence intersection time.

Download:

Table 2. Model parameters.

https://doi.org/10.1371/journal.pone.0191017.t002

Analysis of RMSE

A linear mixed-effects model was also used to test the null hypothesis that neither forest stand characteristics, topography, nor physical obstructions affected GNSS accuracy. The model was fitted using the lmer function in the R lme4 package [58] using the day on which sampling occurred as a random effect to account for variation between days, as might occur due to changing satellite availability and geometry. Variables included as fixed effects in this model were total vegetation (TV_mean), TBA, TPH, Ht, QMD, slope, aspect, and the presence or absence of forest roads, streams, convex slopes, and concave slopes (Table 2). TV_mean for each stand was calculated by averaging the total distance of vegetation (TV_i) from the three LOS paths within each stand. Also, because variables recorded as either present or absent (forest roads, streams, convex slopes, and concave slopes) were quantified along each LOS path, these variables were also considered to be present in this stand-level analysis if they were present along any of the LOS paths. The response variable was the PT RMSE in each stand (RMSE_i) calculated using Eq (4).

Analysis of geofence intersection alert delay

To test the null hypothesis that neither forest stand characteristics, topography, nor physical obstructions affected the time-to-signal accuracy of geofence crossings, a linear mixed-effects model was used to quantify relationships between the magnitude of geofence intersection alert delay as a function of forest stand characteristics, topographic structure, and vegetative LOS obstructions. The model was fitted using the lmer function in the R lme4 package [58] using the day on which sampling occurred as a random effect to account for variation between days, as might occur due to changing satellite availability and geometry. Variables included as fixed effects were total vegetation (TV_mean), TBA, TPH, Ht, QMD, slope, aspect, and the presence or absence of forest roads, streams, convex slopes, and concave slopes (Table 2). TV_mean for each stand was calculated by averaging the total distance of vegetation (TV_i) from the three LOS paths within each stand. Also, because variables recorded as either present or absent (forest roads, streams, convex slopes, and concave slopes) were quantified along each LOS path, these variables were also considered to be present in this stand-level analysis if they were present along any of the LOS paths. The response variable was the overall intersection alert delay in each stand (D_i) calculated using Eq (5).

Model selection

For each of the three analyses, a full model was first fitted to the data using all fixed effect terms. These fixed effects were removed one at a time in order of highest p-values. The aictab function in the R AICcmodavg package [59] was used to compare all resulting models and the model with the lowest corrected Akaike Information Criterion (AICc) was selected. In terms of the geofence alert delay model, the first two models with the lowest AICc had fixed effect terms that were not significant (p-values > 0.05), in which case the model with the third lowest AICc was selected as the final model because all fixed effects had p-values ≤ 0.05. Inferences about all three final models were made using the lincon function in the R trtools package [60], which provides point estimates, standard errors, 95% confidence intervals, and p-values for each model term. In the case of the logistic regression model, the point estimates, standard errors, and confidence intervals were exponentiated to represent the effect of each variable on the odds of getting a missed signal.

Results

Analysis of missed radio signals

The proportion of missed radio signals ranged from 0/90 to 20/90, with a mean of 3.30/90. The mixed-effects logistic regression model with the lowest AICc had total distance of vegetation along the LOS path (TV_i), TPH*0.01, convex, stream, road, and aspect as fixed effects (Table 3). All fixed effects affected the odds of missed signals (p ≤ 0.05). The odds of a missed signal decreased by a factor of 0.93 per unit increase in TV_i (p = 5.65 × 10⁻⁷), while the odds of a missed signal increased by a factor of 1.10 per unit increase in TPH*0.01 (p = 1.45 × 10⁻³). The odds of a missed signal were 1.61 times higher when a slope was convex vs. not convex (p = 3.58 × 10⁻²) and 2.00 times higher in the presence of roads (p = 2.34 × 10⁻⁵). In the presence of streams, the odds of a missed signal decreased by a factor of 0.66 (p = 3.16 × 10⁻²). The odds of a missed signal were 1.05 times higher on east, 2.17 times higher on north, and 2.92 times higher on west aspects (as compared to south aspects), although this effect was only significant on north (p = 9.99 × 10⁻³) and west (p = 1.65 × 10⁻⁴) aspects.

Download:

Table 3. Summary of mixed-effects logistic regression model using stand as a random effect and the odds of missed position updates as the response.

https://doi.org/10.1371/journal.pone.0191017.t003

Analysis of RMSE

RMSE ranged from 1.81 m to 16.69 m, with a mean of 6.61 m. For the RMSE analysis, the mixed-effects model with the lowest AICc had Ht and QMD as fixed effects, both of which affected RMSE (p ≤ 0.05) (Table 4). The RMSE increased as Ht increased (p = 9.15 × 10⁻⁶) but varied indirectly with QMD (p = 3.02 × 10⁻³). Fig 4 illustrates the relationships between predicted RMSE as a function of the two explanatory variables included in the final mixed-effects model.

Download:

Fig 4. Mixed-effects model predictions for PT RMSE.

Predicted RMSE as a function of the two model variables (Ht and QMD). Predictions for each variable were made using the mean of the other predictor. 95% confidence intervals computed using the bootstrap are shown as colored bands. Points on each plot represent partial residuals.

https://doi.org/10.1371/journal.pone.0191017.g004

Download:

Table 4. Summary of mixed-effects linear regression model using day as a random effect and RMSE as the response.

https://doi.org/10.1371/journal.pone.0191017.t004

Analysis of geofence intersection alert delay

Geofence intersection alert delay ranged from −5.33 s to 66 s, with a mean of 18.62 s. The final mixed-effects model used to analyze the delay had TBA, concave, and aspect as fixed effects (Table 5). This model had the third lowest AICc, but was chosen because the two models with lower AICc values each had model terms that were not significant (p > 0.05). The geofence intersection alert delay increased as TBA increased (p = 6.49 × 10⁻⁴) and was also higher in the presence of concave slopes when compared to slopes that were not concave (p = 2.46 × 10⁻²). Finally, the delay was smaller on east, north, and south aspects (as compared to west aspects), although this effect was only significant on east aspects (p = 1.83 × 10⁻²). Fig 5 illustrates the relationships between the predicted alert delay as a function of the three explanatory variables included in the final mixed-effects model.

Download:

Fig 5. Mixed-effects model predictions for geofence intersection alert delay.

Predicted delay as a function of the three model variables (TBA, concave, and aspect). Predictions for each variable were made using the mean of the other predictors. 95% confidence intervals computed using the bootstrap are shown as colored bands. Points on each plot represent partial residuals.

https://doi.org/10.1371/journal.pone.0191017.g005

Download:

Table 5. Summary of mixed-effects linear regression model using day as a random effect and geofence intersection alert delay as the response.

https://doi.org/10.1371/journal.pone.0191017.t005

Discussion

Analysis of missed radio signals indicated that forest stand characteristics, topography, and LOS obstructions affected the odds of missed signals. The odds of missed radio signals varied directly with stand density (TPH) and varied indirectly with LOS obstructions (TV_i). Because previous work has shown that radio signal attenuation increases with greater vegetation depth and density, this result is somewhat counterintuitive [42–45]. However, the magnitude of these effects on the odds ratio was relatively small and may be a result of an interaction between stand density and vegetation in the LOS and/or external influences from wind or humidity [42,43,46–48]. Furthermore, topography influenced radio signal propagation, as convex slopes, aspect, and the presence of forest roads and streams affected the odds of missed signals. The mean number of missed signals during the observed 90-s interval was small, suggesting that relatively few position updates are missed at short distances (100 m).

Analysis of the PT RMSE indicated that only stand variables affected stationary GNSS accuracy. RMSE varied directly with Ht, but decreased with increasing QMD. Because characteristics associated with increasing forest stand density, such as canopy cover, are known to reduce GNSS accuracy, we expected that both predictors would have positive relationships with RMSE. Thus, this result is counterintuitive and may be an artefact of the data. A few sampled stands had relatively open canopies with large, mature trees at low density that may have affected the relationship between RMSE and QMD.

Both stand and topographic variables affected geofence intersection alert delay. Because of the way time delays were calculated, positive delays represent late alerts while negative delays represent early alerts. The alert delay varied directly with TBA and was higher in the presence of concave slopes. Aspect also affected alert delay, with delay being smaller on east slopes compared to west slopes. Taken together these results show that the time-to-signal accuracy of GNSS-RF geofence crossings is affected by both GNSS accuracy and radio signal propagation.

The relationships among the response and predictor variables were not strong for any of the three models. In the case of the missed radio signals, this could be due to the fact that the proportion of missed signals was generally low and the distance between radios fairly small (100 m). In terms of alert delay and RMSE, it is important to note that PTs receive coordinates from NAVSTAR GPS satellites only. Newer GNSS devices receive coordinates from multiple satellite constellations (i.e., from GPS, GLONASS, and BeiDou), which may improve accuracy and reliability in forested environments [61,62]. We collected data for a single geofence crossing in each stand, and RMSE measurements were collected for 5 min at each location. It is possible that clearer relationships would be evident if RMSE was determined using data collected over a longer period of time. Furthermore, future work should perform similar experiments using GNSS-RF transponders capable of processing multiple satellite constellations, as this may yield stronger patterns.

Our results suggest that GNSS-RF radio signal propagation is related to stand density, topography, and obstructions in the line-of-sight and that geofence alert timing is related to stand characteristics and topography. This indicates that the accuracy and successful sharing of GNSS coordinates may change depending on stand conditions and topography, both of which vary on active timber sales. Thus, real-time positioning based on consumer-grade GNSS-RF units may improve general communication and situational awareness on logging operations by allowing ground workers and equipment operators to view the relative positions of nearby workers and machines in real-time on mobile devices. However, high-resolution, mission-critical safety applications of this technology (e.g., geofencing) are not yet advisable under mature forest conditions. Future work should focus on the development of correction methods that account for the effects of forest stand characteristics on GNSS accuracy and geofence alert delay. Previous work has shown that adjustments should also be made for the angle and speed at which a tracked object approaches a geofence [36,37]. Such correction methods may improve reliability and enable a broader range of uses of GNSS-RF technology for safety applications in natural resources, especially if combined with GNSS-RF transponders that communicate with multiple satellite systems. Moreover, phone-based GNSS positioning that incorporates inertial navigation system (INS) capabilities has the potential for improving positioning accuracy. When paired with emerging Bluetooth-based augmentation technologies to create ad-hoc networks, these mobile-based solutions may be promising alternatives to dedicated GNSS-RF units for real-time positioning in remote environments.

Supporting information

S1 Table. Dataset used in analysis of missed radio signals.

https://doi.org/10.1371/journal.pone.0191017.s001

(XLSX)

S2 Table. Dataset used in analysis of RMSE.

https://doi.org/10.1371/journal.pone.0191017.s002

(XLSX)

S3 Table. Dataset used in analysis of geofence intersection alert delay.

https://doi.org/10.1371/journal.pone.0191017.s003

(XLSX)

Acknowledgments

The authors would like to thank four field technicians, Molly Rard, Rebecca Ramsey, Kevin Cannon, and Andrew Naughton, for their valuable assistance with experiment setup, stand inventory, and data collection.

References

1. nasa.gov [Internet]. NASA; c2012–2015 [cited 2017 Jul 3]. Global positioning system history; [about 1 screen]. Available from: http://www.nasa.gov/directorates/heo/scan/communications/policy/GPS_History.html.
2. Wilkie DS. Performance of a backpack GPS in a tropical rain forest. Photogramm Eng Remote Sensing. 1989;55(12): 1747–1749.
- View Article
- Google Scholar
3. Evans DL, Carraway RW, Simmons GT. Use of global positioning system (GPS) for forest plot location. South J Appl For. 1992;16(2): 67–70.
- View Article
- Google Scholar
4. Bilodeau J, Gosselin S, Lowell K, Edwards G, Gold C. Operational integration of GPS into forest management activities which use a GIS to monitor forest operations. Proceedings of the GIS’93 Symposium; 1993 Feb 15–18; Vancouver, BC. p. 196–200.
5. Veal MW, Taylor SE, McDonald TP, McLemore DK, Dunn MR. Accuracy of tracking forest machines with GPS. Trans ASAE. 2001;44(6): 1903–1911.
- View Article
- Google Scholar
6. McDonald TP, Carter EA, Taylor SE. Using the global positioning system to map disturbance patterns of forest harvesting machinery. Can J For Res. 2002;32(2): 310–319.
- View Article
- Google Scholar
7. McDonald TP, Fulton JP. Automated time study of skidders using global positioning system data. Comput Electron Agric. 2005;48(1): 19–37.
- View Article
- Google Scholar
8. Becker RM, Keefe RF, Anderson NM. Use of real-time GNSS-RF data to characterize the swing movements of forestry equipment. Forests. 2017;8(2): 44.
- View Article
- Google Scholar
9. Devlin GJ, McDonnell K. Performance accuracy of real-time GPS asset tracking systems for timber haulage trucks travelling on both internal forest road and public road networks. Int J For Eng. 2009;20(1): 45–49.
- View Article
- Google Scholar
10. Gallo R, Grigolato S, Cavalli R, Mazzetto F. GNSS-based operational monitoring devices for forest logging operation chains. J Agric Eng. 2013;44(2s): 140–144.
- View Article
- Google Scholar
11. Kaartinen H, Hyyppä J, Vastaranta M, Kukku A, Jaakkola A, Yu X, et al. Accuracy of kinematic positioning using global satellite navigation systems under forest canopies. Forests. 2015;6(9): 3218–3236.
- View Article
- Google Scholar
12. Johnson CE, Barton CC. Where in the world are my field plots? Using GPS effectively in environmental field studies. Front Ecol Environ. 2004;2(9): 475–482.
- View Article
- Google Scholar
13. Wing MG, Eklund A, Kellogg LD. Consumer-grade global positioning system (GPS) accuracy and reliability. J For. 2005;103(4): 169–173.
- View Article
- Google Scholar
14. Karsky D. Comparing four methods of correcting GPS data: DGPS, WAAS, L-band, and postprocessing. Missoula (MT): U.S. Department of Agriculture, Forest Service, Missoula Technology and Development Center; 2004 Jul. Report No.: 0471-3-2307-MTDC.
15. D’Eon RG, Serrouya R, Smith G, Kochanny CO. GPS radiotelemetry error and bias in mountainous terrain. Wildl Soc Bull. 2002;30(2): 430–439.
- View Article
- Google Scholar
16. D’Eon RG, Delparte D. Effects of radio-collar position and orientation on GPS radio-collar performance, and the implications of PDOP in data screening. J Appl Ecol. 2005;42(2): 383–388.
- View Article
- Google Scholar
17. Kos T, Markezic I, Pokrajcic J. Effects of multipath reception on GPS positioning performance. Proceedings of the ELMAR-2010 Symposium; 2010 Sep 15–17; Zadar, Croatia. IEEE; 2010.
18. Liu CJ, Brantigan R. Using differential GPS for forest traverse surveys. Can J For Res. 1995;25(11): 1795–1805.
- View Article
- Google Scholar
19. Deckert C, Bolstad PV. Forest canopy, terrain, and distance effects on global positioning system point accuracy. Photogramm Eng Remote Sensing. 1996;62(3): 317–321.
- View Article
- Google Scholar
20. Sigrist P, Coppin P, Hermy M. Impact of forest canopy on quality and accuracy of GPS measurements. Int J Remote Sens. 1999;20(18): 3595–3610.
- View Article
- Google Scholar
21. Holden NM, Martin AA, Owende PMO, Ward SM. A method for relating GPS performance to forest canopy. Int J For Eng. 2001;12(2): 51–56.
- View Article
- Google Scholar
22. Lewis JS, Rachlow JL, Garton EO, Vierling LA. Effects of habitat on GPS collar performance: Using data screening to reduce location error. J Appl Ecol. 2007;44(3): 663–671.
- View Article
- Google Scholar
23. Wempe AM, Keefe RF. Characterizing rigging crew proximity to hazards on cable logging operations using GNSS-RF: Effect of GNSS positioning error on worker safety status. Forests. 2017;8(10): 357.
- View Article
- Google Scholar
24. Liu J, Hyyppä J, Yu X, Jaakkola A, Liang X, Kaartinen H, et al. Can global navigation satellite system signals reveal the ecological attributes of forests? Int J Appl Earth Obs Geoinf. 2016;50: 74–79.
- View Article
- Google Scholar
25. Liu J, Hyyppä J, Yu X, Jaakkola A, Kukko A, Kaartinen H, et al. A novel GNSS technique for predicting boreal forest attributes at low cost. IEEE Trans Geosci Remote Sens. 2017;55(9): 4855–4867.
- View Article
- Google Scholar
26. Bolstad P, Jenks A, Berkin J, Horne K, Reading WH. A comparison of autonomous, WAAS, real-time, and post-processed global positioning systems (GPS) accuracies in northern forests. North J Appl For. 2005;22(1): 5–11.
- View Article
- Google Scholar
27. Piedallu C, Gégout J-C. Effects of forest environment and survey protocol on GPS accuracy. Photogramm Eng Remote Sensing. 2005;71(9): 1071–1078.
- View Article
- Google Scholar
28. Wing MG, Eklund A. Performance comparison of a low-cost mapping grade global positioning systems (GPS) receiver and consumer grade GPS receiver under dense forest canopy. J For. 2007;105(1): 9–14.
- View Article
- Google Scholar
29. Danskin SD, Bettinger P, Jordan TR, Cieszewski C. A comparison of GPS performance in a southern hardwood forest: Exploring low-cost solutions for forestry applications. South J Appl For. 2009;33(1): 9–16.
- View Article
- Google Scholar
30. Bettinger P, Fei S. One year’s experience with a recreation-grade GPS receiver. Math Comput For Nat Resource Sci. 2010;2(2): 153–160.
- View Article
- Google Scholar
31. Danskin S, Bettinger P, Jordan T. Multipath mitigation under forest canopies: A choke ring antenna solution. For Sci. 2009;55(2): 109–116.
- View Article
- Google Scholar
32. Wing MG. Consumer-grade global positioning systems (GPS) receiver performance. J For. 2008;106(4): 185–190.
- View Article
- Google Scholar
33. Andersen HE, Clarkin T, Winterberger K, Strunk J. An accuracy assessment of positions obtained using survey- and recreational-grade global positioning system receivers across a range of forest conditions within the Tanana Valley of interior Alaska. West J Appl For. 2009; 24(3): 128–136.
- View Article
- Google Scholar
34. Wing MG. Consumer-grade GPS receiver measurement accuracy in varying forest conditions. Res J For. 2011;5(2): 78–88.
- View Article
- Google Scholar
35. Keefe RF, Eitel JUH, Smith AMS, Tinkham WT. Applications of multi transmitter GPS-VHF in forest operations. Proceedings of the 47th International Symposium on Forestry Mechanization and 5th International Forest Engineering Conference; 2014 Sep 23–26; Gerardmer, France.
36. Grayson LM, Keefe RF, Tinkham WT, Eitel JUH, Saralecos JD, Smith AMS, et al. Accuracy of WAAS-enabled GPS-RF warning signals when crossing a terrestrial geofence. Sensors. 2016;16(6): 912. pmid:27322287
- View Article
- PubMed/NCBI
- Google Scholar
37. Zimbelman EG, Keefe RF, Strand EK, Kolden CA, Wempe AM. Hazards in motion: Development of mobile geofences for use in logging safety. Sensors. 2017;17(4): 822. pmid:28394303
- View Article
- PubMed/NCBI
- Google Scholar
38. Sygnatur EF. Logging is perilous work. Compens Work Cond. 1998;3: 3–9.
- View Article
- Google Scholar
39. U.S. Bureau of Labor Statistics, Current Population Survey, Census of Fatal Occupational Injuries, 2016. Available from: http://www.bls.gov/iif/oshcfoi1.htm.
40. Tamir T. Radio wave propagation along mixed paths in forest environments. IEEE Trans Antennas Propag. 1977;25(4): 471–477.
- View Article
- Google Scholar
41. Seker SS. Radio pulse transmission along mixed paths in a stratified forest. IEE Proceedings H (Microwaves, Antennas and Propagation). 1989;136(1): 13–18.
- View Article
- Google Scholar
42. Savage N, Ndzi D, Seville A, Vilar E, Austin J. Radio wave propagation through vegetation: Factors influencing signal attenuation. Radio Sci. 2003;38(5): 1088.
- View Article
- Google Scholar
43. Tewari RK, Swarup S, Roy MN. Radio wave propagation through rain forests of India. IEEE Trans Antennas Propag. 1990;38(4): 433–449.
- View Article
- Google Scholar
44. Al-Nuaimi MO, Stephens RBL. Measurements and prediction model optimisation for signal attenuation in vegetation media at centimetre wave frequencies. IEE Proceedings–Microwaves, Antennas and Propagation. 1998;145(3): 201–206.
- View Article
- Google Scholar
45. Azevedo JAR, Santos FES. An empirical propagation model for forest environments at tree trunk level. IEEE Trans Antennas Propag. 2011;59(6): 2357–2367.
- View Article
- Google Scholar
46. Meng YS, Lee YH, Ng BC. The effects of tropical weather on radio-wave propagation over foliage channel. IEEE Transactions on Vehicular Technology. 2009;58(8): 4023–4030.
- View Article
- Google Scholar
47. Löw K. UHF measurements of seasonal field-strength variations in forests. IEEE Transactions on Vehicular Technology. 1988;37(3): 121–124.
- View Article
- Google Scholar
48. Meng YS, Lee YH, Ng BC. Study of propagation loss prediction in forest environment. Prog Electromagn Res B Pier B. 2009;17: 117–133.
- View Article
- Google Scholar
49. Kürner T, Cichon DJ, Wiesbeck W. Concepts and results for 3D digital terrain-based wave propagation models: An overview. IEEE Journal on Selected Areas in Communications. 1993;11(7): 1002–1012.
- View Article
- Google Scholar
50. Freeland RS, Buschermohle MJ, Nichols WM. Precision agriculture: RTK base-to-tractor range limitations using RF communication. Appl Eng Agric. 2014;30(1): 97–104.
- View Article
- Google Scholar
51. Raveon Support [Internet]. Raveon Technologies; [cited 2017 Jul 14]. Available from: http://www.raveon.com/support/.
52. The National Map Download [Internet]. United States Geological Survey (USGS); [cited 2017 Jul 14]. Available from: https://viewer.nationalmap.gov/basic/.
53. NIST Internet Time Servers [Internet]. National Institute of Standards and Technology; [cited 2017 Jul 14]. Available from: http://tf.nist.gov/tf-cgi/servers.cgi.
54. Arrow 100 Subfoot GNSS [Internet]. Eos Positioning Systems; c2014–2017 [cited 2017 Jul 14]. Available from: http://www.eos-gnss.com/arrow-100/.
55. R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria, 2016. Available from: https://www.R-project.org/.
56. Lutes DC, Keane RE, Caratti JF, Key CH, Benson NC, Sutherland S, et al. FIREMON: Fire effects monitoring and inventory system. Fort Collins (CO): U.S. Department of Agriculture, Forest Service, Rocky Mountain Research Station; 2006 Jun. Report No.: RMRS-GTR-164-CD.
57. GNSS Planning Online [Internet]. Trimble; [cited 2017 Jul 5]. Available from: http://www.trimble.com/GNSSPlanningOnline/#/Settings.
58. Bates D, Maechler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67(1): 1–48.
- View Article
- Google Scholar
59. Mazerolle MJ. AICcmodavg: Model selection and multimodel inference based on (Q)AIC(c). 2016. R package version 2.1–0. Available from: https://cran.r-project.org/package=AICcmodavg.
60. Johnson T. trtools: Miscellaneous tools for teaching statistics. 2017. R package version 0.1.0. Available from: http://github.com/trobinj/trtools.
61. Næsset E, Bjerke T, Øvstedal O, Ryan LH. Contributions of differential GPS and GLONASS observations to point accuracy under forest canopies. Photogramm Eng Remote Sensing. 2000;66(4): 403–407.
- View Article
- Google Scholar
62. Li X, Zhang X, Ren X, Fritsche M, Wickert J, Schuh H. Precise positioning with current multi-constellation global navigation satellite systems: GPS, GLONASS, Galileo and BeiDou. Sci Rep. 2015;5: 1–14. pmid:25659949
- View Article
- PubMed/NCBI
- Google Scholar

[ref1] 1. nasa.gov [Internet]. NASA; c2012–2015 [cited 2017 Jul 3]. Global positioning system history; [about 1 screen]. Available from: http://www.nasa.gov/directorates/heo/scan/communications/policy/GPS_History.html.

[ref2] 2. Wilkie DS. Performance of a backpack GPS in a tropical rain forest. Photogramm Eng Remote Sensing. 1989;55(12): 1747–1749.
View Article
Google Scholar

[3] View Article

[4] Google Scholar

[ref3] 3. Evans DL, Carraway RW, Simmons GT. Use of global positioning system (GPS) for forest plot location. South J Appl For. 1992;16(2): 67–70.
View Article
Google Scholar

[6] View Article

[7] Google Scholar

[ref4] 4. Bilodeau J, Gosselin S, Lowell K, Edwards G, Gold C. Operational integration of GPS into forest management activities which use a GIS to monitor forest operations. Proceedings of the GIS’93 Symposium; 1993 Feb 15–18; Vancouver, BC. p. 196–200.

[ref5] 5. Veal MW, Taylor SE, McDonald TP, McLemore DK, Dunn MR. Accuracy of tracking forest machines with GPS. Trans ASAE. 2001;44(6): 1903–1911.
View Article
Google Scholar

[10] View Article

[11] Google Scholar

[ref6] 6. McDonald TP, Carter EA, Taylor SE. Using the global positioning system to map disturbance patterns of forest harvesting machinery. Can J For Res. 2002;32(2): 310–319.
View Article
Google Scholar

[13] View Article

[14] Google Scholar

[ref7] 7. McDonald TP, Fulton JP. Automated time study of skidders using global positioning system data. Comput Electron Agric. 2005;48(1): 19–37.
View Article
Google Scholar

[16] View Article

[17] Google Scholar

[ref8] 8. Becker RM, Keefe RF, Anderson NM. Use of real-time GNSS-RF data to characterize the swing movements of forestry equipment. Forests. 2017;8(2): 44.
View Article
Google Scholar

[19] View Article

[20] Google Scholar

[ref9] 9. Devlin GJ, McDonnell K. Performance accuracy of real-time GPS asset tracking systems for timber haulage trucks travelling on both internal forest road and public road networks. Int J For Eng. 2009;20(1): 45–49.
View Article
Google Scholar

[22] View Article

[23] Google Scholar

[ref10] 10. Gallo R, Grigolato S, Cavalli R, Mazzetto F. GNSS-based operational monitoring devices for forest logging operation chains. J Agric Eng. 2013;44(2s): 140–144.
View Article
Google Scholar

[25] View Article

[26] Google Scholar

[ref11] 11. Kaartinen H, Hyyppä J, Vastaranta M, Kukku A, Jaakkola A, Yu X, et al. Accuracy of kinematic positioning using global satellite navigation systems under forest canopies. Forests. 2015;6(9): 3218–3236.
View Article
Google Scholar

[28] View Article

[29] Google Scholar

[ref12] 12. Johnson CE, Barton CC. Where in the world are my field plots? Using GPS effectively in environmental field studies. Front Ecol Environ. 2004;2(9): 475–482.
View Article
Google Scholar

[31] View Article

[32] Google Scholar

[ref13] 13. Wing MG, Eklund A, Kellogg LD. Consumer-grade global positioning system (GPS) accuracy and reliability. J For. 2005;103(4): 169–173.
View Article
Google Scholar

[34] View Article

[35] Google Scholar

[ref14] 14. Karsky D. Comparing four methods of correcting GPS data: DGPS, WAAS, L-band, and postprocessing. Missoula (MT): U.S. Department of Agriculture, Forest Service, Missoula Technology and Development Center; 2004 Jul. Report No.: 0471-3-2307-MTDC.

[ref15] 15. D’Eon RG, Serrouya R, Smith G, Kochanny CO. GPS radiotelemetry error and bias in mountainous terrain. Wildl Soc Bull. 2002;30(2): 430–439.
View Article
Google Scholar

[38] View Article

[39] Google Scholar

[ref16] 16. D’Eon RG, Delparte D. Effects of radio-collar position and orientation on GPS radio-collar performance, and the implications of PDOP in data screening. J Appl Ecol. 2005;42(2): 383–388.
View Article
Google Scholar

[41] View Article

[42] Google Scholar

[ref17] 17. Kos T, Markezic I, Pokrajcic J. Effects of multipath reception on GPS positioning performance. Proceedings of the ELMAR-2010 Symposium; 2010 Sep 15–17; Zadar, Croatia. IEEE; 2010.

[ref18] 18. Liu CJ, Brantigan R. Using differential GPS for forest traverse surveys. Can J For Res. 1995;25(11): 1795–1805.
View Article
Google Scholar

[45] View Article

[46] Google Scholar

[ref19] 19. Deckert C, Bolstad PV. Forest canopy, terrain, and distance effects on global positioning system point accuracy. Photogramm Eng Remote Sensing. 1996;62(3): 317–321.
View Article
Google Scholar

[48] View Article

[49] Google Scholar

[ref20] 20. Sigrist P, Coppin P, Hermy M. Impact of forest canopy on quality and accuracy of GPS measurements. Int J Remote Sens. 1999;20(18): 3595–3610.
View Article
Google Scholar

[51] View Article

[52] Google Scholar

[ref21] 21. Holden NM, Martin AA, Owende PMO, Ward SM. A method for relating GPS performance to forest canopy. Int J For Eng. 2001;12(2): 51–56.
View Article
Google Scholar

[54] View Article

[55] Google Scholar

[ref22] 22. Lewis JS, Rachlow JL, Garton EO, Vierling LA. Effects of habitat on GPS collar performance: Using data screening to reduce location error. J Appl Ecol. 2007;44(3): 663–671.
View Article
Google Scholar

[57] View Article

[58] Google Scholar

[ref23] 23. Wempe AM, Keefe RF. Characterizing rigging crew proximity to hazards on cable logging operations using GNSS-RF: Effect of GNSS positioning error on worker safety status. Forests. 2017;8(10): 357.
View Article
Google Scholar

[60] View Article

[61] Google Scholar

[ref24] 24. Liu J, Hyyppä J, Yu X, Jaakkola A, Liang X, Kaartinen H, et al. Can global navigation satellite system signals reveal the ecological attributes of forests? Int J Appl Earth Obs Geoinf. 2016;50: 74–79.
View Article
Google Scholar

[63] View Article

[64] Google Scholar

[ref25] 25. Liu J, Hyyppä J, Yu X, Jaakkola A, Kukko A, Kaartinen H, et al. A novel GNSS technique for predicting boreal forest attributes at low cost. IEEE Trans Geosci Remote Sens. 2017;55(9): 4855–4867.
View Article
Google Scholar

[66] View Article

[67] Google Scholar

[ref26] 26. Bolstad P, Jenks A, Berkin J, Horne K, Reading WH. A comparison of autonomous, WAAS, real-time, and post-processed global positioning systems (GPS) accuracies in northern forests. North J Appl For. 2005;22(1): 5–11.
View Article
Google Scholar

[69] View Article

[70] Google Scholar

[ref27] 27. Piedallu C, Gégout J-C. Effects of forest environment and survey protocol on GPS accuracy. Photogramm Eng Remote Sensing. 2005;71(9): 1071–1078.
View Article
Google Scholar

[72] View Article

[73] Google Scholar

[ref28] 28. Wing MG, Eklund A. Performance comparison of a low-cost mapping grade global positioning systems (GPS) receiver and consumer grade GPS receiver under dense forest canopy. J For. 2007;105(1): 9–14.
View Article
Google Scholar

[75] View Article

[76] Google Scholar

[ref29] 29. Danskin SD, Bettinger P, Jordan TR, Cieszewski C. A comparison of GPS performance in a southern hardwood forest: Exploring low-cost solutions for forestry applications. South J Appl For. 2009;33(1): 9–16.
View Article
Google Scholar

[78] View Article

[79] Google Scholar

[ref30] 30. Bettinger P, Fei S. One year’s experience with a recreation-grade GPS receiver. Math Comput For Nat Resource Sci. 2010;2(2): 153–160.
View Article
Google Scholar

[81] View Article

[82] Google Scholar

[ref31] 31. Danskin S, Bettinger P, Jordan T. Multipath mitigation under forest canopies: A choke ring antenna solution. For Sci. 2009;55(2): 109–116.
View Article
Google Scholar

[84] View Article

[85] Google Scholar

[ref32] 32. Wing MG. Consumer-grade global positioning systems (GPS) receiver performance. J For. 2008;106(4): 185–190.
View Article
Google Scholar

[87] View Article

[88] Google Scholar

[ref33] 33. Andersen HE, Clarkin T, Winterberger K, Strunk J. An accuracy assessment of positions obtained using survey- and recreational-grade global positioning system receivers across a range of forest conditions within the Tanana Valley of interior Alaska. West J Appl For. 2009; 24(3): 128–136.
View Article
Google Scholar

[90] View Article

[91] Google Scholar

[ref34] 34. Wing MG. Consumer-grade GPS receiver measurement accuracy in varying forest conditions. Res J For. 2011;5(2): 78–88.
View Article
Google Scholar

[93] View Article

[94] Google Scholar

[ref35] 35. Keefe RF, Eitel JUH, Smith AMS, Tinkham WT. Applications of multi transmitter GPS-VHF in forest operations. Proceedings of the 47th International Symposium on Forestry Mechanization and 5th International Forest Engineering Conference; 2014 Sep 23–26; Gerardmer, France.

[ref36] 36. Grayson LM, Keefe RF, Tinkham WT, Eitel JUH, Saralecos JD, Smith AMS, et al. Accuracy of WAAS-enabled GPS-RF warning signals when crossing a terrestrial geofence. Sensors. 2016;16(6): 912. pmid:27322287
View Article
PubMed/NCBI
Google Scholar

[97] View Article

[98] PubMed/NCBI

[99] Google Scholar

[ref37] 37. Zimbelman EG, Keefe RF, Strand EK, Kolden CA, Wempe AM. Hazards in motion: Development of mobile geofences for use in logging safety. Sensors. 2017;17(4): 822. pmid:28394303
View Article
PubMed/NCBI
Google Scholar

[101] View Article

[102] PubMed/NCBI

[103] Google Scholar

[ref38] 38. Sygnatur EF. Logging is perilous work. Compens Work Cond. 1998;3: 3–9.
View Article
Google Scholar

[105] View Article

[106] Google Scholar

[ref39] 39. U.S. Bureau of Labor Statistics, Current Population Survey, Census of Fatal Occupational Injuries, 2016. Available from: http://www.bls.gov/iif/oshcfoi1.htm.

[ref40] 40. Tamir T. Radio wave propagation along mixed paths in forest environments. IEEE Trans Antennas Propag. 1977;25(4): 471–477.
View Article
Google Scholar

[109] View Article

[110] Google Scholar

[ref41] 41. Seker SS. Radio pulse transmission along mixed paths in a stratified forest. IEE Proceedings H (Microwaves, Antennas and Propagation). 1989;136(1): 13–18.
View Article
Google Scholar

[112] View Article

[113] Google Scholar

[ref42] 42. Savage N, Ndzi D, Seville A, Vilar E, Austin J. Radio wave propagation through vegetation: Factors influencing signal attenuation. Radio Sci. 2003;38(5): 1088.
View Article
Google Scholar

[115] View Article

[116] Google Scholar

[ref43] 43. Tewari RK, Swarup S, Roy MN. Radio wave propagation through rain forests of India. IEEE Trans Antennas Propag. 1990;38(4): 433–449.
View Article
Google Scholar

[118] View Article

[119] Google Scholar

[ref44] 44. Al-Nuaimi MO, Stephens RBL. Measurements and prediction model optimisation for signal attenuation in vegetation media at centimetre wave frequencies. IEE Proceedings–Microwaves, Antennas and Propagation. 1998;145(3): 201–206.
View Article
Google Scholar

[121] View Article

[122] Google Scholar

[ref45] 45. Azevedo JAR, Santos FES. An empirical propagation model for forest environments at tree trunk level. IEEE Trans Antennas Propag. 2011;59(6): 2357–2367.
View Article
Google Scholar

[124] View Article

[125] Google Scholar

[ref46] 46. Meng YS, Lee YH, Ng BC. The effects of tropical weather on radio-wave propagation over foliage channel. IEEE Transactions on Vehicular Technology. 2009;58(8): 4023–4030.
View Article
Google Scholar

[127] View Article

[128] Google Scholar

[ref47] 47. Löw K. UHF measurements of seasonal field-strength variations in forests. IEEE Transactions on Vehicular Technology. 1988;37(3): 121–124.
View Article
Google Scholar

[130] View Article

[131] Google Scholar

[ref48] 48. Meng YS, Lee YH, Ng BC. Study of propagation loss prediction in forest environment. Prog Electromagn Res B Pier B. 2009;17: 117–133.
View Article
Google Scholar

[133] View Article

[134] Google Scholar

[ref49] 49. Kürner T, Cichon DJ, Wiesbeck W. Concepts and results for 3D digital terrain-based wave propagation models: An overview. IEEE Journal on Selected Areas in Communications. 1993;11(7): 1002–1012.
View Article
Google Scholar

[136] View Article

[137] Google Scholar

[ref50] 50. Freeland RS, Buschermohle MJ, Nichols WM. Precision agriculture: RTK base-to-tractor range limitations using RF communication. Appl Eng Agric. 2014;30(1): 97–104.
View Article
Google Scholar

[139] View Article

[140] Google Scholar

[ref51] 51. Raveon Support [Internet]. Raveon Technologies; [cited 2017 Jul 14]. Available from: http://www.raveon.com/support/.

[ref52] 52. The National Map Download [Internet]. United States Geological Survey (USGS); [cited 2017 Jul 14]. Available from: https://viewer.nationalmap.gov/basic/.

[ref53] 53. NIST Internet Time Servers [Internet]. National Institute of Standards and Technology; [cited 2017 Jul 14]. Available from: http://tf.nist.gov/tf-cgi/servers.cgi.

[ref54] 54. Arrow 100 Subfoot GNSS [Internet]. Eos Positioning Systems; c2014–2017 [cited 2017 Jul 14]. Available from: http://www.eos-gnss.com/arrow-100/.

[ref55] 55. R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria, 2016. Available from: https://www.R-project.org/.

[ref56] 56. Lutes DC, Keane RE, Caratti JF, Key CH, Benson NC, Sutherland S, et al. FIREMON: Fire effects monitoring and inventory system. Fort Collins (CO): U.S. Department of Agriculture, Forest Service, Rocky Mountain Research Station; 2006 Jun. Report No.: RMRS-GTR-164-CD.

[ref57] 57. GNSS Planning Online [Internet]. Trimble; [cited 2017 Jul 5]. Available from: http://www.trimble.com/GNSSPlanningOnline/#/Settings.

[ref58] 58. Bates D, Maechler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67(1): 1–48.
View Article
Google Scholar

[149] View Article

[150] Google Scholar

[ref59] 59. Mazerolle MJ. AICcmodavg: Model selection and multimodel inference based on (Q)AIC(c). 2016. R package version 2.1–0. Available from: https://cran.r-project.org/package=AICcmodavg.

[ref60] 60. Johnson T. trtools: Miscellaneous tools for teaching statistics. 2017. R package version 0.1.0. Available from: http://github.com/trobinj/trtools.

[ref61] 61. Næsset E, Bjerke T, Øvstedal O, Ryan LH. Contributions of differential GPS and GLONASS observations to point accuracy under forest canopies. Photogramm Eng Remote Sensing. 2000;66(4): 403–407.
View Article
Google Scholar

[154] View Article

[155] Google Scholar

[ref62] 62. Li X, Zhang X, Ren X, Fritsche M, Wickert J, Schuh H. Precise positioning with current multi-constellation global navigation satellite systems: GPS, GLONASS, Galileo and BeiDou. Sci Rep. 2015;5: 1–14. pmid:25659949
View Article
PubMed/NCBI
Google Scholar

[157] View Article

[158] PubMed/NCBI

[159] Google Scholar

Figures

Abstract

Introduction

Materials and methods

Field experiment

Analysis of missed radio signals

Analysis of RMSE

Analysis of geofence intersection alert delay

Model selection

Results

Analysis of missed radio signals

Analysis of RMSE

Analysis of geofence intersection alert delay

Discussion

Supporting information

S1 Table. Dataset used in analysis of missed radio signals.

S2 Table. Dataset used in analysis of RMSE.

S3 Table. Dataset used in analysis of geofence intersection alert delay.

Acknowledgments

References