1. Introduction
The population of the Bathurst barren ground caribou (Rangifer tarandus groenlandicus, hereafter described as “Bathurst Caribou”) has declined 93% since 1986 [1] . Declines have also been reported for other barren ground caribou and wild reindeer herds around the globe [2] .
A combination of many natural and human factors (e.g., habitat, harvest, pre- dators, diseases, parasites, insects, extreme weather, climate change, industrial development, and pollution) may have contributed to the decline [3] [4] [5] . Among these factors, climate change and an increase in industrial development have been notable across the Arctic in recent decades [6] - [12] . The climate in the Arctic has been changing at almost twice of the global average rate during the last century [6] . Under changing climate, the area of lichen-rich mature forest over the Bathurst Caribou winter range has decreased as the burned area increases with warmer summer temperatures [7] . The accessibility to lichen, the main diet of caribou during the winter, may have also been influenced by an increase in thaw-freeze events and ice content in snow [7] . On the summer range, the high variability of forage availability during the post-calving period and fall may result in a decrease in net caribou productivity, namely, the late-winter calf: cow ratio [8] . Also, a warming climate may prolong the duration of insect activity and thus negatively impact caribou energetics, health, and growth [9] [10] .
There has been an increase in mining exploration, operation, and reclamation activities within the Bathurst Caribou range since the 1980s. Active mines include the Ekati Diamond Mine, Diavik Diamond Mine, Snap Lake Diamond Mine, and GahchoKue Diamond Mine [11] [12] . Also, there are several mining exploration projects, as well as reclamation projects for closed mines [11] . Traditional knowledge indicates that some of the potential effects of mining operations on caribou might include the shifting of caribou migration routes away from mining areas, and the deterioration of caribou health (e.g., abnormal material and smell in caribou meat, liver, and the hide linings) [13] . Boulanger et al. [14] estimated that there were fewer caribou with a zone of influence (ZOI) of about 14 km surrounding the Ekati-Diavik mining complex based on caribou aerial survey data and locations of satellite-collared cows.
However, little is known about the mechanisms that may produce a ZOI. Observations by local community members [13] suggested that the sight of, noise from, and road dust released by mining operationsmight influence caribou. In addition to road dust, other sources of airborne particulate matter may include diesel fuel combustion, living waste incineration, and the secondary chemical reactions in the atmosphere [15] . According to their sizes, particulate matter can be classified as fine particulate matter (PM2.5 with aerodynamic diameter < 2.5 μm) and coarse particulate matter (aerodynamic diameter ≥ 2.5 μm or the total suspended particulates diameter TSP-PM2.5) [15] . The concentrations of TSP and PM2.5 are usually monitored through air sampling [15] . Their deposition rates were also monitored at some locations. For example, the deposition rate of TSP (or dust deposition) was monitored at a number of locations around the Ekati Diamond Mine [16] . The depositions of both coarse and fine particulate matters might impact tundra ecosystem properties (e.g., soil pH) as well as caribou forage quality and growth. In this study, our objective is to answer the question: how have the Ekati mining operations affected the availability and quality of forage for caribou?
2. Materials and Methods
2.1. Study Area
The Bathurst caribou herd ranges over a vast area from 63˚50'N - 66˚27'N and 107˚31'W - 114˚43'W determined by locations obtained from satellite-collared cows from 1996 to 2012 [17] . Bathurst caribou calve to the west and southwest of Bathurst Inlet, winter below treeline in the Northwest Territories and the northern portion of Saskatchewan, and during the summer forage on tundra areas in between. The Ekati Diamond Mine is located in the middle the summer range of the Bathurst caribou herd (Figure 1).
Field surveys were conducted near the Ekati Diamond Mine during August 12-24, 2015, and August 15-29, 2016. The field surveys were conducted along six transects (Figure 1). Five of the transects (red lines) were along the active portion of the Misery Haul Road, indicated by km markers relative to the Main Camp and direction (3 km west, 5 km east, 8 km west, 20 km south, and 20.2 km south). Transects 20 km south and 20.2 km south are approximately 200 m apart, but their vegetation distributions were sufficiently different to be treated as two independent transects. One transect was sampled at the rarely used spur road outside the Misery Camp, and is shown in green in Figure 1.
2.2. Field Survey Methods
GPS position was measured using two Trimble NetR9 GNSS Reference Receivers coupled with Trimble antennas. These receivers are capable of recording both GPS and Glonass data. One receiver was set at a base station near a road and acquired data in static mode throughout the day. The other receiver and its antenna were placed in a backpack and carried along a transect for a total daily distance of se- veral kms. Both receivers were acquiring data simultaneously at a1-second inter- val. Base station data were also acquired at a 30-second interval for processing by NRCan Precision Point Processing service. Post processing generated tens of
Figure 1. Location of vegetation field measurement transects (red and green lines) around the Ekati diamond mine inside the Ba- thurst caribou habitat, which is located in the bordering region of NWT, Nunavut, and Saskatchewan.
thousands of 3-dimensional points (i.e., latitude, longitude, and elevation). Planimetric results are in NAD 83 CSRS, while altimetry (orthometric) is in CGDV28.
Sampling sites were located at 0.1, 5, 10, 30, 100, 200, 300, 500, 750, 1000, 1250, 1500 m along each transect. The percent cover for each vascular plant species, lichen, and moss, bare ground, and rock was visually estimated along each transect following methods described in [18] . The mean height of each vascular plant species was measured using a ruler for five representative plants per site. Down-looking digital photos with regular visible bands were then taken. These down-looking photos provide a more reliable estimate of percent cover [19] , and were used to correct any potential observer bias. Down-looking NDVI photos with the near-infrared band and visible bands were also taken at select sites. Leblanc et al. [20] have developed a method to derive surface reflectance information from digital photos, which can be used to better link field measurements with satellite remote sensing imagery and to assess if remote sensing can detect dust accumulation on leaves.
Dwarf birch (Betula nana) was the most common plant species at all these sites and was selected as an indicator species for quantifying potential effects of dust from mining operations on the quality of caribou forage. Leaves from dwarf birch were collected at each site and the amount of dust present on leaves was analyzed at the Environment Laboratory of the Ekati Diamond Mine.
The soil sample was collected from each site to determine soil pH. The top 0 to 5 cm of soil is considered the layer most affected by airborne contaminants such as dust, and thus is the layer of most interest for both human and wildlife health risk assessments [21] . We collected a 50 to 100 g soil sample from the top 5 cm of the soil layer, and near the dominant plant species at the site, which was usually dwarf birch. The sample was taken back (in a paper bag if the sample was relatively dry or a plastic Ziplock bag if the sample was wet) to the Environment Laboratory of the Ekati Diamond Mine to analyze soil pH.
2.3. Laboratory Methods
To measure the amount of dust on dwarf birch leaves in the laboratory, we weighed about 1 gram of air-dried leaves from a site and submerged them into 100 ml of distilled water. The leaves were well stirred and left in the water for 10 minutes to allow the dust to dissolve into the water. Three 10 ml samples of the water were used to measure the turbidity of the solution using an EPA compliant (http://www.coleparmer.ca/) LaMotte LTC3000we Benchtop Turbidity and Chlorine Meter. Results are in Nephelometric Turbidity Units (NTU) per 100 ml water per gram of air-dry leaves.
Soil pH was measured using a soil pH meter (HACH H135 Advanced Compact Water PR by Cole-Parmer Canada, http://www.hach.com/, or YSI Professional Plus Multiparameter Water Quality Instrument, https://www.ysi.com/proplus). The analysis followed these steps:
・ A batch of 8 to 10 beakers (300 ml) was washed using distilled water.
・ Approximately 10 - 20 g of a soil sample from a site was added to a beaker.
・ Distilled water was added to the beaker to fully cover the soil sample.
・ The solution was shaken and then allowed to settle for approximately 5 minutes.
・ A digital photo of the solution was taken, together with a site label.
・ The pH meter calibrated with standard solutions at pH = 4, 7, and 10.
・ The pH meter was then washed with distilled water after calibration and shaken dry.
・ The soil solution was transferred to a 50 ml beaker, and the pH meter was sub- merged into it and allowed to stabilize (usually <5 minutes). The pH value for the site was then recorded.
・ The above steps were repeated so that 5 soil pH measurements were conducted for each site.
2.4. Data Analysis Methods
To quantify the impacts of mining related dust on caribou forage, we had to overcome a number of challenges. First is the large natural variation in the availability of caribou forage caused by numerous factors that could mask the effects of mine related dust deposition. For example, Gould and Walker [22] reported that distribution and richness of vegetation along the Hood River might vary greatly due to influences of topography, disturbance, and soil pH. The Hood River is located at the northern extent of the Bathurst caribou range and has not been subject to industrial development. In this study, we stratified all sites into seven land cover classes: boulder field, esker, lichen-rock tundra, dwarf shrub, low-tall shrub, tussock, and graminoid wetland. Within a given land cover class; differences associated with topography, disturbance, and soil pH tend to be reduced. Therefore, we could minimize the natural variation byanalyzing the impacts of mining related dust on caribou forage within a particular land cover class. Only the dwarf shrub class had a sufficient sample size for statistical analysis.
The second challenge is the remaining within-class natural variation, particularly in the extent of % rock cover. Vascular plants and lichen that caribou feed on do not grow on the bare surface of rocks. As a result, variation in the amount of rock cover between dwarf shrub sites could significantly bias the results. We addressed this challenge by removing the rock fraction in a site from the analysis and the percent cover of a given plant species is calculated for the remaining area:
(1)
The third challenge is that coarse and fine particulate matter released from mi- ning operations may impact caribou forage differently [15] . The coarse particulate matter is mechanically produced by the break-up of larger materials. For the su- rrounding area around the Ekati Diamond Mine, dust is mainly from traffic on unpaved roads constructed using crushed granite or wind-blown over the unco- vered soil. The dominant components of granite are silicon (72% weight), aluminium (14%), potassium (4%), sodium (4%), iron (3%), calcium (2%), and magnesium (1%) [23] . Metal analysis of dust samples collected from dust deposition stations near the haul roads during 2012-2014 show similar results. The major metal components of dust were aluminum (8%), potassium (6%), sodium (2%), calcium (4%), and magnesium (8%) besides 72% silicon [16] . Silicon dioxide and aluminum oxide are not water soluble under natural environment conditions and thus neutral in pH. The pH values of other metal solutions in water (e.g., potassium hydroxide) and their carbonated compounds (e.g., potassium carbonate, potassium hydrogen carbonate) range from 8 to 12 (http://www.aqion.de/site/191). Consequently, a typical pH value of road dust from a haul road is expected to be about 9. At the Ekati Diamond Mine, road dust suppression agents have been applied on the Misery Haul Road in recently years. The pH values of road dust suppression agents are typically acidic. As a result, the pH of dust from the Misery Haul Road may vary tempo-spatially depending on the amount, timing, and types of dust suppression agents applied.
PM2.5 contains most of the emissions from fossil fuel combustion, secondary particulates, biomass burn, and the finer components of road dust [23] . The se- condary particulates are produced by the chemical reactions in the atmosphere. For example, sulfur dioxide is oxidized in the atmosphere to form sulfuric acid (H2SO4). Nitrogen dioxide (NO2) is oxidized to nitric acid (HNO3), which in turn reacts with ammonia (NH3) to form ammonium nitrate (NH4NO3). The com- positions of PM2.5 at the Edmonton McIntyre Air Monitoring Station in 2007 were measured to be organic carbon (45% weight), elemental carbon (15%), sulfur di- oxide (17%), nitrate (9%), and Ammonium (10%) [24] . The corresponding acids (i.e., sulfuric acid, nitric acid, and carbon acid) have pH values ranging from 1 to 4, comparable with the finding that the pH values of PM2.5 are usually < 3 [25] .
Because of their larger particulate sizes, road dust tends to deposit onto the land surface within a shorter distance and timeframe, relative to PM2.5. With different pH values, they can also affect soil pH differently. Regarding direct impacts on vegetation beside the change of soil pH, an acidic PM2.5 solution containing hea- vy metals may result in a decrease in lichen survival rate [26] . On the other hand, road dust may affect vegetation photosynthesis, respiration, transpiration when dust covers a leaf surface [27] . In this study, we addressed this challenge by investigating two contrasting transects, namely, those from a haul road vs. those largely dust-free from a rarely used spur road.
3. Results
3.1. Topography, Disturbance, and Land Cover Classes
The topography surrounding the Ekati Diamond Mine is a gently rolling landscape. Along the survey transects, elevation typically varied less than 20 m within 2.5 km of the mine (Figure 2).
Lichen-rock tundra class was found on the tops of hills at some sites. Hill tops tend to be drier than slopes or down in the valley (Table 1). The wettest sites are usually located at valley bottoms, where graminoid wetland, tussock, and low- high shrub classes are typically found. The dwarf shrub class occupies most of the slope sites, as well as some of the sites at the hilltop and valley bottoms. Two other classes (boulder field and esker) represent the legacy of glaciation and deglaciation processes in this region [28] . Dominated by coarse materials such as boulders, gravels, and sands, these classes are driest in terms of moisture regime (Table 1). Among all sites sampled along the six transects, dwarf shrub was the most prevalent class (81%), followed by esker (5%), and the low-high shrub and tussock classes (4% each). As a result, the dwarf shrub class was selected as an appropriate indicator for the potential effects of mine related dust on caribou forage.
3.2. The Amount of Dust on Leaves
The amount of dust per gram of air-dried dwarf birch leaves decreased with distance from the Misery Haul Road (Figure 3); however, there was substantial variation between transects. For the majority of transects from the Misery Haul Road, the highest amount of dust on leaves was recorded for sites closest to the road. Some exceptions may occur. For example, the highest amount of dust on
Figure 2. Elevation change on the east transect from the misery haul road at 5 km benchmark from the main camp, based on the precision GPS survey on August 19, 2016.
Table 1. Description of land cover classes over the area surrounding the Ekati diamond mine.
leaves collected from the 8 km west transect was found at 100 m, indicating the high special variability in the measurements of dust on leaves. In contrast, the highest dust deposition was always observed by the Ekati Diamond Mine [16] .
In addition to topography and winds, rain events may further affect the amount of dust on leaves. Our in-situ observations suggest that traffic can stir dust from an unpaved road within 12 hours after a rain event. The longer the duration after a rain event, the higher the amount of dust may be stirred up by traffic. All surveys were conducted after rain events of differing intensity and durations, and the amount of dust on leaves represents an accumulation of dust deposition after a rain event. With high variability in topography, wind direction and speed, rain
Figure 3. The amount of dust on leaves of dwarf birch plotted against the distance from the misery haul road (3 km west, 8 km west, 5 km east, 20 km south, and 20.2 km south), as well as from a rarely used end road outside of the misery camp.
intensity, and the time interval between rain events and the field survey, a single measurement of dust on leaves may be less meaningful. With this in mind, we averaged the amount of dust on leaves over a range of distances along all transects from the Misery Haul Road (Table 2).
The average amount of dust on leaves at sites within 10 m was 8.8 ± 2.2 times higher than at sites >1500 m from the Misery Haul Road (Table 2). A significantly higher amount of dust was recorded from leaves collected from >500 m to 1000 m, which was 2.6 ± 0.9 times higher than leaves collected from sites >1500 m from the Misery Haul Road. However, there was no statistically significant difference in the amount of dust from leaves collected from >1000 to 1500 m and from >1500 m, indicating the maximum distance of higher than background amount of dust on leaves was about 1000 m. Our finding is comparable with the direct dust deposition monitoring results [16] .
In contrast, there was no significant difference in the amounts of dust among sites on the transect from a rarely used spur road outside of the Misery Camp (Figure 3 and Table 2). This result confirms our finding that the amount of dust on leaves at sites within 1000 m of the Misery Haul Road was significantly higher than background level.
3.3. Impacts on Soil pH
The average pH value of dust samples taken from transects along the Misery Haul Road was 8.63 ± 0.44, ranging from 8.1 to 9.4 (Figure 4). The variation in dust pH values may be related to differences in the application of dust suppressing
Table 2. The ratio of the amount of dust on leaves averaged over sites within a given distance range to that >1500 m (a) from the Misery Haul Road, as well as (b) from a rarely used end road outside of the Misery Camp. The error values are one standard estimation error relative to the mean dust amount for sites >1500 m.
agents, construction material, and usage history. The area with the heaviest application of dust suppression agents is the airstrip. The pH value of samples collected from the airstrip was 6.1 ± 0.2, while the pH values for dust sampled from haul roads with no application of dust suppression agents ranged from 9 to 10. The variation in pH values for the dust samples collected from transects along the Misery Haul Road fell between the air strip and haul road without dust suppression, and likely related to the difference in the amount, timing, and types of dust suppression agents applied.
The pH values of dust were much higher than those of tundra soils not disturbed by natural erosion or affected by dust deposition, typically ranging from 4 to 6 [22] . The large difference in pH values between road dust and tundra soils suggests that there should be a gradient in soil pH from basic to acidic pH values as the distance from the haul road increases. A significant relationship between the soil pH and the distance from the Misery Haul Road was observed (Figure 5). Soil pH decreased logarithmically from about 9 at the roadside to about 4.5 at a distance around 1000 - 1500 m, beyond this distance, the soil pH fluctuated around 4.5.
In contrast, a decrease in soil pH with distance was not observed along the transect from a rarely used spur road outside the Misery Camp (Figure 6). In fact, we found a significantly increasing trend with distance for this transect. Simulation results from a CALPUFF air dispersion model show that the concentration of PM2.5 tends to be higher at sites close to mining operations [14] , resulting in decreased soil pH at these sites. As such, the increasing trend with distance for this transect could be explained by the influence of PM2.5, which has a pH < 3 [25] . However, we noticed that the difference in pH values between PM2.5 and natural tundra soil was significantly less than that between road dust and natural tundra soil of the dwarf shrub class. In addition, the concentration of coarse particulates (TSP-PM2.5) tends to be much higher than that of PM2.5 [16] . As a result, the impact of coarse particulates was found to be more dominant than that of PM2.5. Nevertheless, the two contrasting trends are indicative of the different impacts of
Figure 4. pH value of dust samples collected at roadside every 2 km along the misery haul road on August 20, 2015.
Figure 5. Soil pH change along the distance from the misery haul road, with the fit line given by y = −0.424ln(x) + 7.5446, R2 = 0.90, p-value < 0.01, and the sample size n = 65. The p-value was computed between measured and estimated pH values.
coarse particulates from road dust and PM2.5 on soil pH.
3.4. Relationships between Soil pH and Vegetation % Cover
We observed a significant 2nd order polynomial relationship between the percent cover of vascular plants and soil pH (Figure 7). The optimum soil pH values for vascular plants were found to range from 4.5 to 6. As soil pH increased above 6,
Figure 6. Soil pH change with distance for the transect from a rarely used end road outside of the misery camp, with the fit line given by y = 0.00015x + 4.5197, R2 = 0.26, p- value < 0.05, and the sample size n = 23.
a significant reduction in the percent cover of vascular plants was observed.
In contrast, lichen was observed to prefer acidic soil (Figure 8). As soil pH increased from 4 to 9, the percent cover of lichen decreased.
3.5. Impacts on Vegetation Percentage Cover
A significant reduction in vascular plant percent cover occurred within the first 10 m from the Misery Haul Road (Figure 9). This reduction corresponds to soil pH values > 6 at sites near the Misery Haul Road (Figure 5). In addition, the amount of dust on leaves in this zone near the road is about nine times above background levels, which can contribute to a reduction in photosynthesis rates at these nearby sites and result in poor growth near the road.
There were no changes in the percent cover of vascular plants with distance from the Misery Haul road beyond the first 10 m. The same can be said for the transect at the rarely used spur road outside of the Misery Camp. This result is expected given soil pH values beyond the first 10 m were generally in the range of 4.5 - 6, the optimum soil pH range for the growth of arctic vascular plants (Figure 7). The higher amount of dust on leaves could potentially be negative for vascular plants but appeared to be masked by the large natural variation in vascular plant cover among sites.
Lichen distribution tends to be more sensitive to dust deposition. A significant reduction in the percent cover of lichen was observed at sites near the Misery Haul Road (Figure 10). The reduction in lichen cover corresponds to increased soil pH at sites within 1000 m of the Misery Haul Road and is indicative of the potential negative association between raised soil pH and lichen cover.
Figure 7. Relationship between soil pH and vascular plant % cover for the transects from the Misery Haul Road, with the fit line given by y = −6.2513x2 + 64.172x − 73.851, R2 = 0.84, p-value < 0.01, and the sample size n = 65. The p-value was computed between measured and estimated vascular % cover. Results for the transectfrom a rarely used spur road outside the misery camp were also plotted.
Figure 8. Relationship between soil pH and lichen % cover for the transects from the Misery Haul Road, given by y = 0.0672x4 − 2.2663x3 + 28x2 − 151.03x + 302.49, R2 = 0.41, p-value < 0.01, and n = 65. The p-value was computed between measured and estimated lichen % cover. Results for the transectfrom a rarely used end road outside the misery campwere also plotted.
For the transect from the rarely used spur road outside of the Misery Camp, no significant trend in the percent cover of lichen with distance was observed (Figure 10). If soil pH were the sole factor controlling the distribution of lichen, we would have expected a higher lichen % cover at sites near the spur road than those at >2.5 km (Figure 6 and Figure 8). The fact no such an increase in lichen % cover was observed suggests that the PM2.5 deposition may have some
Figure 9. Vascular % cover changes with distance for the transects from the misery haul road. Results for the transect from a rarely used end road outside of the misery camp were also plotted. The lines show the 5-point running mean.
Figure 10. Lichen % cover changes with distance from a mining operation. The red line shows the 5-point running mean, while dashed lines for linear trends: y = 0.0137x + 1.0609, R2 = 0.5075, p-value < 0.01, and n = 48 for the first 1000 m from the Misery Haul Road; and y = −0.0005x + 7.1191, R2 = 0.02, p-value = 0.5, and n = 23 for the transect from a rarely used end road outside of the misery camp.
negative impact on lichen survival [26] , and thus partially canceled out the potential positive effect of soil pH decrease on lichen cover near the Misery Camp.
4. Discussion
4.1. Uncertainties in Field Survey Data
The analysis of dust content on leaves may be subject to several error sources. For example, when we picked the leaves from plants in the field, some of the dust may be shaken off. Alternatively, if we would cut a branch or twigs in the field with minimal disturbance, we still had to pick out leaves in the laboratory for turbidity analysis and lose some dust in the process. More dust may be lost during transport from the field to the laboratory, during the air-drying process, and when the leaf sample was weighed and used for turbidity measurement. Also, the time since last rainfall may not have been equal for all sites. Finally, it is possible that dust on leaves may not be completely dissolved in the water solution during the turbidity measurement procedure. Consequently, the amount of dust on leaves quantified by the turbidity measurement method will not give us the absolute amount of dust on leaves. Instead, they should be used as a relative indicator of the amount of dust on leaves.
Similarly, the analysis of soil pH may also subject to errors. One of the main uncertainties in soil pH results is the possible variation in the soil depth from which soil samples were collected. Soil from naturally disturbed sites may also have a high pH value [20] . If a soil sample was collected at a deeper layer, we might find a higher soil pH. We attempted to collect soil samples in the top 5 cm; however, the actual depth from which the soil sample was collected may vary from site to site, and contribute to the variability in soil pH between sites. Another potential error may occur during the measurement of soil pH. The pH meters we used were calibrated and accurate to within 0.1 pH value. Uncertainty in the pH value may occur if dust in the water solution settled down to the bottom. Maintaining the soil solution by gently stirring the mixture during the mea- surement is recommended.
The percent cover of vascular plants and lichen was visually estimated in the field and later corrected using digital site photos [17] . While visual estimation is subjective, the digital photo analysis method is more accurate and provides some verification of field results. The digital photo analysis method can be time-con- suming and can introduce additional error if there are overlaps between different vegetation types. Trends and means derived from these site measurements may be more reliable.
4.2. Comparability with Previous Studies
The significant influence of soil pH on the distribution of vascular plants and lichen was also reported by the Hood River study [22] . Despite the difference in vegetation parameters used in the two studies (species richness in their study and percent cover in our study), both suggested a negative effect on vascular plants where soil pH values are either too high or too low. Near the Ekati Diamond Mine, the optimum soil pH for vascular plants appears to be 4.5 to 6. Beyond this soil pH range, the percent cover of vascular plants decreases, especially at higher pH values nearer to the road associated with road dust. Additionally, both studies indicated that lichen increases as soil pH decreases from 9 to 4.
In terms of road dust impacts on soil pH and thus vegetation distribution, our results also agreed well with previous studies [29] [30] . For example, Myers-Smith et al. [29] reported a statistically significant logarithmic relationship between the distance from the 577-km-long Dalton Highway in Alaska and the soil pH out to 1 km from the road, similar to our finding. They also reported that plots adjacent to the road had significantly higher graminoid and cloudberry (Rubus chamaemorus) biomass, but less moss, evergreen shrub, lichen, and for b biomass. Our results showed different impacts on specific species, which is understandable given that the dominant plant functional types in between their study area and ours were somewhat different. Nevertheless, the reduction of lichen due to road dust was observed by both studies.
4.3. Interpretation of Impact Assessment
Given some of the potential uncertainties in the measurement of dust on leaves, soil pH, and percent cover of vegetation, the results are an index of the potential effects of mine related dust on caribou forage and should be interpreted with caution. Conclusions drawn from a single parameter are less reliable. Instead, the inference should rely on the convergence of results derived from several parameters.
For example, our results indicate the zone of dust on leaves for the Misery Haul Road was about 1 km. Dustfall monitoring results also show an increase in dust deposition up to 1 km [23] . A significant increase in soil pH from the Misery Haul Road occurred within 1 to 1.5 km. Lichen coverage was also reduced within 1 km of the reduction from the Misery Haul Road. Integrating all these results together, we may conclude with confidence that the zone of affected forage (ZOAF) is about 1 km from the Misery Haul Road.
Our results are comparable with several regional scale landscape studies that reporting that caribou reduced the use of areas within 5 km from infrastructure and human activity [31] . In contrast, the 14 km ZOI around the Ekati-Diavik mining complex identified by [14] appears to be much larger that the ZOAF observed in this study. On the other hand, we emphasize that this study deals only with the zone of disturbed vegetation, and that of the coarse particulate matter (i.e. TSP-PM2.5). The zones of noise, visibility of mining, and PM2.5 could potentially reach further and represent important areas to focus future research in order to further refine our understanding of their individual and combined contributions to the extent and intensity of ZOI.
5. Conclusions
On the basis of two years of field data and laboratory analyses conducted on samples collected from the Ekati Diamond Mine in the Bathurst caribou range in Canada’s Arctic, we draw the following conclusions.
・ Precision GPS measurements show the area surrounding the Ekati Diamond Mine is a gently rolling topographic landscape, with a range in elevation from valley to peak of < 40 m for all transects, with a hill or two approximately every 1 km along a transect.
・ Dwarf shrub is the predominant land cover class across the study area, accounting for 81% coverage of all sites, followed by esker (5%), and the low-high shrub and tussock classes (4% each).
・ The average amount of dust measured on dwarf birch leaves collected from sites at 0 - 10, >10 - 100, >100 - 500, >500 - 1000, and >1000 - 1500 m from the Misery Haul Road was 8.8, 6.1, 4.2, 2.6, and 1.1 times than that from >1.5 km, respectively. These results indicated that the zone of raised dust on leaves was about 1 km from the Misery Haul Road.
・ Road dust increased soil pH logarithmically, from about 4.5 at 1 - 1.5 km dis- tance to about 9 at the roadside of the Misery Haul Road.
・ The optimum soil pH for vascular plants is around 4.5 - 6. The increase in soil pH values > 6 within 10 m of the Misery Haul Road reduced vascular percent cover. An increase of dust on leaves in the zone may further contribute to this reduction.
・ Lower soil pH (to about 4) favors lichen, which corresponds to a reduction in lichen cover that is observed within ~1 km of the Misery Haul Road.
・ These distance related trends in dust deposition rate, the amount of dust on leaves, the soil pH values, and the vegetation percent covers were not observed for sites on the transect from a rarely used spur road outside the Misery Camp.
Overall, these results suggest that the ZOAF from the Misery Haul Road is about 1 km. On the other hand, no such a zone exists for the rarely used spur road outside the Misery Camp, indicating the ZOAF is subject to road use. Further investigation is needed for all other mining roads to fully quantify the ZOAF of the Ekati Diamond Mine within the Bathurst caribou summer range.
Finally, our results are comparable with the ZOI estimates of 1 - 5 km [31] , but much less than the 14 km of ZOI around the Ekati-Diavik mining complex identified by [14] . Investigation of other mechanisms (e.g., the zones of noise, visibility of mining, and PM2.5) is essential to further refine our understanding of ZOI.
Acknowledgements
The study was funded by the Northwest Territories Cumulative Impact Monitoring Program (NWT CIMP) and Natural Resources Canada. Dominion Diamond Ekati Corporation (DDEC) provided logistic supports, including travel, accommodation, laboratory facility for measuring soil pH and dust on leaves, and safety training. Wildlife technicians from DDEC (Matt Hoover, Cody Drygeese, Jeff Mantla, and Misty Sinclair) assisted us in the field surveys. Three co-op students (Charlotte Kelly of Carleton University, Anumeet Garchaand Holden Ciufo of University of Waterloo) participated the fieldwork.