Numerical Terradynamic Simulation Group 11-2016 Mapping Arctic Tundra Vegetation Communities Using Field Spectroscopy and Multispectral Satellite Data in North Alaska , USA

The Arctic is currently undergoing intense changes in climate; vegetation composition and productivity are expected to respond to such changes. To understand the impacts of climate change on the function of Arctic tundra ecosystems within the global carbon cycle, it is crucial to improve the understanding of vegetation distribution and heterogeneity at multiple scales. Information detailing the fine-scale spatial distribution of tundra communities provided by high resolution vegetation mapping, is needed to understand the relative contributions of and relationships between single vegetation community measurements of greenhouse gas fluxes (e.g., ~1 m chamber flux) and those encompassing multiple vegetation communities (e.g., ~300 m eddy covariance measurements). The objectives of this study were: (1) to determine whether dominant Arctic tundra vegetation communities found in different locations are spectrally distinct and distinguishable using field spectroscopy methods; and (2) to test which combination of raw reflectance and vegetation indices retrieved from field and satellite data resulted in accurate vegetation maps and whether these were transferable across locations to develop a systematic method to map dominant vegetation communities within larger eddy covariance tower footprints distributed along a 300 km transect in northern Alaska. We showed vegetation community separability primarily in the 450–510 nm, 630–690 nm and 705–745 nm regions of the spectrum with the field spectroscopy data. This is line with the different traits of these arctic tundra communities, with the drier, often non-vascular plant dominated communities having much higher reflectance in the 450–510 nm and 630–690 nm regions due to the lack of photosynthetic material, whereas the low reflectance values of the vascular plant dominated communities highlight the strong light absorption found here. High classification accuracies of 92% to 96% were achieved using linear discriminant analysis with raw and rescaled spectroscopy reflectance data and derived vegetation indices. However, lower classification accuracies (~70%) resulted when using the coarser 2.0 m WorldView-2 data inputs. The results from this study suggest that tundra vegetation communities are separable using plot-level spectroscopy with hand-held sensors. These results also show that tundra vegetation mapping can be scaled from the plot level (<1 m) to patch level (<500 m) using spectroscopy data rescaled to match the wavebands of the multispectral satellite remote sensing. We find that developing a consistent method for classification of vegetation communities across the flux tower sites is a challenging process, given the Remote Sens. 2016, 8, 978; doi:10.3390/rs8120978 www.mdpi.com/journal/remotesensing Remote Sens. 2016, 8, 978 2 of 24 spatial variability in vegetation communities and the need for detailed vegetation survey data for training and validating classification algorithms. This study highlights the benefits of using fine-scale field spectroscopy measurements to obtain tundra vegetation classifications for landscape analyses and use in carbon flux scaling studies. Improved understanding of tundra vegetation distributions will also provide necessary insight into the ecological processes driving plant community assemblages in Arctic environments.


Introduction
The Arctic is currently undergoing dramatic changes in climate [1][2][3].As such, changes in vegetation composition and productivity are predicted to occur due to a lengthening of the growing season [4][5][6], increases in air and soil temperature and active layer depth [7,8], increased shrub cover [9,10], and acceleration of below ground microbial activity and nutrient cycling [3,11,12].These changes influence the terrestrial carbon cycle and thus alter the relationship between arctic ecosystems and global climate.For example, an increasing growing season length can enhance carbon uptake through photosynthesis and thereby reduce atmospheric carbon dioxide (CO 2 ) concentrations [13].Conversely, increasing temperatures may heighten methane (CH 4 ) and CO 2 release from thawing permafrost, enhancing plant-mediated transport of CH 4 [14,15] and contributing to rising levels of greenhouse gases in the atmosphere.
Carbon cycle models used to inform global climate predictions typically run at a ~1 km to 0.5 • grid cell resolution, yet ecosystem carbon fluxes are most frequently measured at the plot scale (~1 m).Many of these models assume uniform vegetation composition within a cell footprint, yet in reality Arctic landscapes have high spatial heterogeneity [16].A better understanding of the role of vegetation community composition on greenhouse gas exchange is needed to guide scaling from plot levels to coarser resolution footprints [17][18][19].For example, Shaver et al. [20] note that changes in vegetation community composition could lead to higher ecosystem productivity, alongside increased long-term sequestration of carbon [1].Furthermore, Ström et al. [21] note that a high proportion of carbon assimilated through photosynthesis is allocated belowground in Arctic environments, which can be impacted by different vegetation community composition.Climate warming and a lengthening growing season length may reduce cold temperature constraints on stored belowground carbon, substantially increasing ecosystem CO 2 respiration or loss as CH 4 [22].Vegetation can be important in regulating CH 4 emissions as plants provide labile substrates for methanogenic production.Tundra plants have structures, such as aerenchyma, which may also provide conduits for transport of CH 4 from the soil to the atmosphere [14].Although the distribution of Arctic tundra vegetation is relatively well established at coarser scales (e.g., 1:7,500,000) and documented in databases such as the Circumpolar Arctic Vegetation Map ((CAVM team, 2003; [23]) and others [24,25], knowledge of local tundra community distributions remains ambiguous at many locations.Understanding the fine-scale spatial distribution of tundra communities, combined with improved high resolution vegetation mapping, is needed to understand relative contributions of and relationships between single vegetation community measurements of greenhouse gas fluxes (e.g., ~1 m chamber flux measurements; [15]) and those encompassing multiple vegetation communities (e.g., ~100 m eddy covariance measurements; [26]).
Remote sensing allows for the assessment of vegetation across the Arctic at a variety of spatial and temporal scales under a rapidly changing climate.Moderate resolution satellite-based multispectral sensors such as the 1 km resolution Advanced Very High Resolution Radiometer (AVHRR), the 250 m to 1 km Moderate Resolution Imaging Spectroradiometer (MODIS) and the 30 m Landsat series imagers Remote Sens. 2016, 8, 978 3 of 24 observe tundra ecosystems at daily to bi-weekly repeat cycles.These remote sensing data have been used successfully to map and monitor regional vegetation change across Arctic tundra systems [27,28], including decadal increases in biomass and greenness [29][30][31].Although coarse spatial resolution imagery (>100 m) may be suitable for determining broad-scale ecosystem distributions and regional change, it is insufficient to capture the high spatial heterogeneity occurring across many Arctic tundra landscapes [16,32,33].High resolution satellite data (<5 m) such as WorldView-2 may be able to capture this spatial heterogeneity, as tundra vegetation communities low stature may make it difficult to identify and quantify using coarser products [19].However, there can be a trade-off between fine and coarser spatial resolution satellite data through temporal sampling frequency.Satellite data are often collected bi-weekly to monthly in comparison to daily repeat cycles using moderate resolution (e.g., MODIS) products.Finally, hyperspectral data have also been used successfully to map tundra vegetation, and has worked across different geographies [28,[34][35][36][37], but it is only available through hand held or airborne sensors.
In complement to satellite remote sensing data, field spectroscopy provides a useful tool to characterize and map vegetation communities [29,38,39], species composition [40,41], biophysical/ ecological properties [3], and plant traits [42,43] from local to landscape scales.Field spectroscopy data are not widely collected and used in Arctic tundra vegetation studies [44] due to the often difficult logistics and weather related constraints for field sampling in remote high latitude environments [45].None-the-less, recent studies using field spectroscopy report improved spectral differentiation of vegetation communities at the plot scale [28,39,46].For example, Bratsch et al. [39] successfully used field spectroscopy to discriminate between four tussock tundra communities at Ivotuk, Alaska with overall classification accuracies of 84%-94%.Langford et al. [47] used LiDAR data in combination with multispectral satellite imagery to map plant functional types (PFTs) at the Barrow-BEO site (also featured in this paper).However, given the tundra variability across the Alaskan landscape it is important to attempt such mapping exercises to multiple sites and combining field spectral measurements with satellite data.Further, and because of the variable conditions of the Arctic, it is important that surveys are carried-out by a consistent team, and that the timings of the surveys and imagery acquisition are well matched.Through this combination, we can test to which extent field spectroscopy can be used to map Arctic tundra communities across large landscape areas.
The goal of this research is to assess whether we could successfully map tundra vegetation communities at plot and landscape scales using detailed vegetation community analysis in conjunction with field spectroscopy and multispectral remote sensing data, across a variety of tundra types.The first objective is to first determine whether dominant tundra vegetation communities found in different locations across the Arctic Coastal Plain are spectrally distinct and distinguishable using field spectroscopy methods.Secondly, to test which combination of raw reflectance and vegetation indices retrieved from field and satellite data allows accurately mapping vegetation and across locations.To test these objectives, we developed a systematic method to map dominant vegetation communities within larger (~300 m) eddy covariance tower footprints distributed along a 300 km transect in Northern Alaska.The ability to map Arctic tundra vegetation at finer spatial scales is necessary to better understand the contribution of tundra communities to the global carbon cycle and scale related dependencies of ecological processes in various plant assemblages.

Study Area
This study was undertaken at four tower eddy covariance field sites across the North Slope of Alaska (Figure 1).Two tower sites were located near Barrow.Two additional sites were located further south towards the Brooks Range, in Atqasuk and Ivotuk.
The Barrow Environmental Observatory (Barrow-BEO; 71 • 16 51.61"N,156 • 36 44.44"W) tower site consists of an ice-wedge polygonal tundra landscape [48], with well-developed high-, low-and flat-center polygons, whereas the other Barrow tower (hereafter Barrow-BES; 71 • 16 51.17"N,156 • 35 47.28"W) is situated in a vegetated drained lake basin [49] approximately 300 m to the east of Barrow-BEO.Mean annual air temperature and mean summer precipitation at the Barrow sites for the 1948-2013 period was −11.3 • C and 72 mm, respectively [50].The vegetation communities at Barrow-BEO and Barrow-BES are dominated by sedges, grasses and various mosses [23].The third site, Atqasuk (70 • 28 40"N, 157 • 25 05"W), is located approximately 100 km south of the Barrow sites.This site includes primarily low centered, well developed polygonal tundra with well-drained high edges [24,51] and has a mean annual air temperature of −10.8 • C and a mean summer precipitation of 100 mm for the 1999-2006 period.The vegetation communities found at this site are dominated by sedges, grasses, mosses and dwarf shrubs (<40 cm) [23].The fourth site, Ivotuk (68.49• N, 155.74 • W), is located in the foothills of the Brooks Range Mountains, approximately 300 km south of Barrow.There is no substantial polygon formation located here, with the site consisting of a gentle northwest facing slope and a wet meadow on the margins of a stream.The mean annual air temperature and mean summer precipitation at this site are −8.9 • C and 210 mm respectively.The vegetation communities at Ivotuk are dominated by tussock sedge, dwarf shrubs and mosses [23].All four tower sites are located on areas of continuous permafrost with an active layer thaw depth of approximately 37 cm [52] for Barrow-BEO/Barrow-BES, 43 cm for Atqasuk [53] and 25 cm for Ivotuk [50].
Remote Sens. 2016, 8, 978 4 of 24 site, Atqasuk (70°28′40″N, 157°25′05″W), is located approximately 100 km south of the Barrow sites.This site includes primarily low centered, well developed polygonal tundra with well-drained high edges [24,51] and has a mean annual air temperature of −10.8 °C and a mean summer precipitation of 100 mm for the 1999-2006 period.The vegetation communities found at this site are dominated by sedges, grasses, mosses and dwarf shrubs (<40 cm) [23].The fourth site, Ivotuk (68.49°N, 155.74°W), is located in the foothills of the Brooks Range Mountains, approximately 300 km south of Barrow.
There is no substantial polygon formation located here, with the site consisting of a gentle northwest facing slope and a wet meadow on the margins of a stream.The mean annual air temperature and mean summer precipitation at this site are −8.9 °C and 210 mm respectively.The vegetation communities at Ivotuk are dominated by tussock sedge, dwarf shrubs and mosses [23].All four tower sites are located on areas of continuous permafrost with an active layer thaw depth of approximately 37 cm [52] for Barrow-BEO/Barrow-BES, 43 cm for Atqasuk [53] and 25 cm for Ivotuk [50].

Vegetation Data
Vegetation surveys were conducted at the four tower sites between 18 and 31 July 2014.The main vegetation types and micro-topographic features (e.g., polygon troughs, rims and centers) were identified by an intensive walkover survey and use of previous vegetation maps of the areas [24,54] within a 300 m tower footprint, chosen as the "fetch", i.e., the distance from the tower that carbon fluxes are expected to originate from, calculated as 100 m multiplied by the height of the tower (in this case, maximum tower height 3.5 m) [50,55].A precautionary approach was taken during the vegetation survey by placing ten 1 m 2 quadrats for which each vegetation type (wet sedge meadow for example) and feature was recorded.In total, this resulted in 5 vegetation types/50 quadrats at Barrow-BEO, 1 vegetation type/10 quadrats at Barrow-BES (which was only accessible from a central boardwalk), 2 vegetation types/30 quadrats (an additional 10 quadrats were placed in one type to capture the variability) at Atqasuk and 3 vegetation types over 30 quadrats at Ivotuk (Figure S1).
The percentage cover of all live vascular and non-vascular plant species and other materials including standing dead, standing water and bare ground was categorized as 0.1 (present), 1 (occasional, with 1 < x < 5 individuals), or 3 (containing more than 10 individuals).The nearest 5% of all live vascular, non-vascular and other materials thereafter were recorded.All vascular and nonvascular plant species identifications were made during the field campaign according to Hultén [56] and Vitt et al. [57] respectively.

Vegetation Data
Vegetation surveys were conducted at the four tower sites between 18 and 31 July 2014.The main vegetation types and micro-topographic features (e.g., polygon troughs, rims and centers) were identified by an intensive walkover survey and use of previous vegetation maps of the areas [24,54] within a 300 m tower footprint, chosen as the "fetch", i.e., the distance from the tower that carbon fluxes are expected to originate from, calculated as 100 m multiplied by the height of the tower (in this case, maximum tower height 3.5 m) [50,55].A precautionary approach was taken during the vegetation survey by placing ten 1 m 2 quadrats for which each vegetation type (wet sedge meadow for example) and feature was recorded.In total, this resulted in 5 vegetation types/50 quadrats at Barrow-BEO, 1 vegetation type/10 quadrats at Barrow-BES (which was only accessible from a central boardwalk), 2 vegetation types/30 quadrats (an additional 10 quadrats were placed in one type to capture the variability) at Atqasuk and 3 vegetation types over 30 quadrats at Ivotuk (Figure S1).
The percentage cover of all live vascular and non-vascular plant species and other materials including standing dead, standing water and bare ground was categorized as 0.1 (present), 1 (occasional, with 1 < x < 5 individuals), or 3 (containing more than 10 individuals).The nearest 5% of all live vascular, non-vascular and other materials thereafter were recorded.All vascular and nonvascular plant species identifications were made during the field campaign according to Hultén [56] and Vitt et al. [57] respectively.
For mapping purposes, dominant plant communities were defined using a combination of cluster analyses and ordination of species composition data.Groups of similar samples were identified using a one-way hierarchical cluster analyses performed using PC-ORD software version 6 (MjM Software, Gleneden Beach, OR, USA).Ordination was then performed using non-metric multi-dimensional scaling (NMDS) as this approach does not assume unimodal species response curves and is well suited to ecological datasets of this type [58].The specifications used in the NMDS were selected according to the characteristics of the dataset and were: (i) Bray-Curtis distance measure; (ii) number of dimensions equals 2 (following examination of the stress plots of three runs); and (iii) data auto-transformation (Wisconsin double and square root transformations).The NMDS was performed in R version 3.1.0(R Core Team 2014, [59]) using the VEGAN community ecology package [60].The resulting plots for each site, on which samples are arranged primarily according to their floristic dissimilarity, were examined visually and compared with the groupings from cluster analysis.
For this analysis, classifications were conducted on the plant community level following the approach of Buchhorn et al. [28] and Bratsch et al. [39].Although PFTs are commonly used in Earth System and other carbon flux models to prescribe vegetation parameters and system response to changing environmental conditions [61], the spatial heterogeneity of the tundra ecosystem is evident even at sub-meter scales, making single PFT classifications not good representatives of vegetation composition, distribution and properties.Furthermore, classifications based on landscape-level PFTs may result in higher classification errors at the tower footprint scale (~300-500 m) relative to classifications based on vegetation community type [28].The vegetation community, chosen as the basis for our classification scheme, also better matches upscaling schemes used in local and regional carbon assessments, where land cover maps are used in conjunction with small area (e.g., 1 m chamber) flux measurements to obtain landscape level estimates of carbon exchange and greenhouse gas emissions [15].Although the classifications resulting from this study do not directly reflect PFT types, our vegetation communities can be re-classified to PFT types for input in carbon flux models (in Section 3.2, a comparison of vegetation communities and PFTs is provided).

Field Spectroscopy Data Collection and Pre-Processing
Field spectroscopy measurements were collected using a UniSpec DC Spectrometer Analysis System (PP systems, Amesbury, MA, USA) having 256 discrete bands spanning from 450 nm to 1040 nm at the four field sites on a single occasion in July 2014 (Table 1).All measurements were collected under near to clear sky conditions and as close to solar noon as possible.For each vegetation type (Section 2.2), 30 spectral points were collected at quadrant locations matching the vegetation survey, within a 300 m radius footprint of the eddy covariance towers.In total, 150 samples were collected at Barrow-BEO (11 July 2014), 30 samples were collected at Barrow-BES (11 July 2014), 60 samples were collected at Atqasuk (29 July 2014), and 90 samples were collected at Ivotuk (16 July 2014) (Figure S1).Vegetation canopy reflectance was measured using downward facing fibre optic (UNI684 straight fibre optic, 2M, HCS LOH, SMA-custom ferrule with 100 mm tip) Hypo Tube (UNI688), which was used to control the field of view (FOV).Spectra were recorded from approximately 1.5 m above the ground with a FOV of 9.5 • , corresponding to a circular surface with an approximate 50 cm diameter.A white standard (UNI240, PP systems, Amesbury, MA, USA) was used for white reference normalization.At each location, a photograph was taken to illustrate the vegetation community and a corresponding differential GPS measurement (dGPS) was acquired (Trimble R7, Trimble Navigation Limited, Sunnyvale, CA, USA).
As the field spectrometer used (UniSpec DC) did not account for incoming radiation, the raw spectral data was corrected using the mean white reference standard normalization measurement taken during each measurement period.One mean reflectance spectrum for each vegetation community (calculated from between 30 and 60 spectral measurements dependent on vegetation community characteristics and sampling quality) was calculated per location.Wavebands that fell outside the range of 450-1040 nm in the raw reflectance data (hereby referred to as UniSpec) were removed during quality control processing.The UniSpec data were then rescaled to match the 7 bands of the WorldView-2 satellite data (hereby referred to as UniSpecWV2) to test whether downgrading narrowband data to multispectral data would still produce accurate tundra vegetation maps.We degraded the UniSpec data to the 7 WorldView-2 bands by averaging band reflectance as follows: Blue: 450-510 nm; Green: 510-580 nm; Yellow; 585-625 nm; Red: 630-690 nm; Red Edge: 705-745 nm; NIR1: 770-895 nm; and NIR2: 860-1040 nm.

Satellite Data Collection and Pre-Processing
The WorldView-2 multispectral data (2 m spatial resolution) acquired for Barrow-BEO and Barrow-BES, Atqasuk and Ivotuk were also used to classify tundra vegetation (Table 1).The image acquisition dates were selected for the time period between vegetation sampling and field spectral measurement campaigns.All sites had imagery for 2014 with the exception of Ivotuk, where imagery from 2013 was used as there was no suitable, cloud-free satellite imagery available during the following year.The orthorectified WorldView-2 data were georeferenced and calibrated to Top-of-Atmosphere Reflectance following the procedures outlined by Digital Globe (Digital Globe Corporate.Longmont, CO, USA).To assess whether we could relate the TOA reflectance to the spectrometer data we used a dark object subtraction and then compared reflectance values at known locations (see Supplementary Material for a detailed description).We extracted the reflectance data from the WorldView-2 satellite imagery using dGPS locations of the field spectrometry measurements (extracted data is hereby referred to as WorldView-2).The spatial accuracy of georeferenced dGPS points was within 30 cm of "truth" using the satellite imagery.

Vegetation Indices
We used the UniSpec, UniSpecWV2 and WorldView-2 reflectance data to calculate three vegetation indices (VIs), namely the normalized difference vegetation index (NDVI) [62], the normalized difference water index (NDWI) [63] and the enhanced vegetation index (EVI) [64].We chose these indices as NDVI is commonly used to describe the quantity of green vegetation at plot and landscape levels [25,29,65].The EVI was included to improve on areas where NDVI is limited, due to canopy structure variation [64,66], i.e., non-homogeneous Arctic tundra vegetation can range from <1 cm to >5 m in height therefore affecting the reflectance signal, and sensitivity to background reflectance from soils, etc. [67].The NDWI estimates plant water stress and associated impacts on plant productivity, phenology, and land cover classifications.This is especially useful in Arctic tundra ecosystems where water table depth and soil moisture can be intrinsically linked to vegetation cover [37,68].We used the field spectroscopy bands 665 nm and 845 nm [69] and the Red and NIR1 from the UniSpecWV2 and WorldView-2 data to calculate NDVI, NDWI and EVI.The advantage of using vegetation indices in combination with raw reflectance data is that the indices themselves relate parts of the spectra to reflect biochemical properties and other properties of the spectra that cannot be accounted for when analyzing each individual band.

Spectral Separability and Image Classification
Our objective was to map the following vegetation communities: mesic-sedge-grass-herb meadow, dry lichen heath, wet sedge meadow, tussock tundra, and mixed shrub-sedge tussock tundra.These were the vegetation community classes resulting from the vegetation clustering analysis (described in Section 3.1).We first tested whether there was adequate spectral separability among tundra vegetation communities to inform a classification algorithm to produce a vegetation map.To do this, we applied a Principal Component Analysis (PCA) to UniSpec data to assess whether the vegetation communities were separable.The PCA identifies whether each vegetation community is in a different area of the reflectance data space, and the loadings of the PCA axes indicate the contribution of reflectance at each of the spectral bands to the final PCA plot [69].Secondly, we used a one way analysis of variance (ANOVA) and a pair-wise post hoc Tukey HSD test (95% confidence level) to assess which wavelengths differed significantly among the vegetation communities.Once the region(s) were defined, a second PCA was used on the selected wavelengths.
Finally, we used a linear discriminant analysis (LDA) to discriminate amongst the vegetation communities.We chose this method because it allows for the identification of significant differences among groups based on relatively large numbers of variables [70], through variable reduction to linear combinations of those variables that explain the majority of variation in the data [71].It has been successfully used to classify vegetation communities and species in several studies [41,72,73], including tundra [39].We classified the communities using input UniSpec data from the significantly different wavelengths, as well as the UniSpecWV2 reflectance and WorldView-2 image reflectance data.We also trained a LDA using a combination of each type of reflectance data and the three VIs.Initially, one classifier was developed across-sites, but this was deemed unsuitable due to low classification accuracies (see Supplementary Material, Figure S3), therefore site-specific classifiers were developed.Two LDA functions were obtained from each analysis, producing three classes each for Barrow-BEO/Barrow-BES and Ivotuk and two classes for Atqasuk.
Classification accuracy was calculated as the number of ground truth points correctly classified in a given vegetation community divided by the number actually in that reference vegetation community.The Kappa coefficient metric [74] was used as a statistical measure summarizing the classification accuracy.Kappa coefficient values range from −1 to 1, with ≤0 indicating no agreement, and 0.01-0.20 as none to slight, 0.21-0.40 as fair, 0.41-0.60 as moderate, 0.61-0.80 as substantial and 0.81-1 as almost perfect agreement [74].
We applied the resulting LDA discriminant functions (UniSpecWV2 and UniSpecWV2 plus VI) to the seven WorldView-2 bands to produce the final vegetation community map for each field site.For each image we ran all the corresponding LDA functions (as many as number of vegetation communities at the site), and images were classified with 4-6 classes.As only two of the many communities found at Atqasuk [75] were found in the tower footprint, one resulting LDA function was produced.Finally, the LDA maps were combined by layering each raster together and the final vegetation classes were mapped.Pair-wise error matrices (and corresponding Kappa coefficients) for each analysis type (LDA including or excluding VIs) and site were calculated based on errors of omission and commission between the classified maps and the ground-truth data.Analysis was undertaken using ENVI v5.2 (Exelis Visual Information Solutions, Boulder, CO, USA) and ArcMap v.10.2.1 (ESRI 2011.ArcGIS Desktop, Environmental Systems Research Institute, Redlands, CA USA).All statistical analyses were carried out using the packages CRAN and MASS in R version 3.1.0(R Core Team 2014) [59].

Vegetation Communities
Hierarchical one-way clustering and NMDS analysis showed a total of eight close groupings, including mesic-sedge-grass-herb meadow, dry lichen heath, wet sedge meadow, tussock tundra (two separate communities at Atqasuk (sandy substrates) and Ivotuk (non-sandy substrates), and mixed shrub-sedge tussock tundra communities (Figure 2).At Barrow-BEO/Barrow-BES, the dry lichen heath community was dominated by Polytrichum moss and various lichen species, with little vascular plants present, while the mesic-sedge-grass herb meadow community was dominated by a mixture of graminoid species (Eriophorum russeolum and Poa arctica) with various mosses, lichens and liverworts present.The wet sedge community here was dominated by Carex aquatilis (sedge) with an underlying Sphagnum spp.and Drepanocladus spp.moss carpet below the canopy.At Atqasuk, the tussock sedge (sandy substrates) community was dominated by Eriophorum vaginatum, Aulocomnion turgidum, Rubus chamaemorus and various evergreen shrubs.The wet sedge community here does not contain a moss layer, but consisted of Eriophorum angustifolium and Eriophoporum russeolum.At Ivotuk, the tussock sedge (non-sandy substrates) also contained Eriophorum vaginatum tussocks interspersed with Sphagnum spp., as well as various evergreen and deciduous shrubs.The wet sedge community here consisted of tall Carex aquatilis, with low-growing deciduous shrubs such as Salix pulchra and a moss carpet containing Sphagnum spp.The mixed shrub-sedge tussock community was also dominated by large Salix pulchra and Betula nana, with Eriophorum vaginatum tussock.Details on species frequency and abundance are provided in Supplementary Material (Table S1) and details on PFT frequency and abundance within all eight communities are provided in Table 2.At Barrow-BEO/Barrow-BES, the dry lichen heath community was dominated by Polytrichum moss and various lichen species, with little vascular plants present, while the mesic-sedge-grass herb meadow community was dominated by a mixture of graminoid species (Eriophorum russeolum and Poa arctica) with various mosses, lichens and liverworts present.The wet sedge community here was dominated by Carex aquatilis (sedge) with an underlying Sphagnum spp.and Drepanocladus spp.moss carpet below the canopy.At Atqasuk, the tussock sedge (sandy substrates) community was dominated by Eriophorum vaginatum, Aulocomnion turgidum, Rubus chamaemorus and various evergreen shrubs.The wet sedge community here does not contain a moss layer, but consisted of Eriophorum angustifolium and Eriophoporum russeolum.At Ivotuk, the tussock sedge (non-sandy substrates) also contained Eriophorum vaginatum tussocks interspersed with Sphagnum spp., as well as various evergreen and deciduous shrubs.The wet sedge community here consisted of tall Carex aquatilis, with low-growing deciduous shrubs such as Salix pulchra and a moss carpet containing Sphagnum spp.The mixed shrub-sedge tussock community was also dominated by large Salix pulchra and Betula nana, with Eriophorum vaginatum tussock.Details on species frequency and abundance are provided in Supplementary Material (Table S1) and details on PFT frequency and abundance within all eight communities are provided in Table 2.

Spectral Separability in Tundra Vegetation Communities
Average spectral profiles for all vegetation communities are shown in Figure 3. Vegetation communities differed in the shape of their spectral signature in the visible region of the spectrum, mostly due to the mesic sedge-grass-herb meadow and dry lichen heath communities showing higher reflectance in the green (about 10%, at 600 nm) and in the red (about 12%, at 680 nm) regions relative to the other vegetation communities.Communities differed in the NIR region, with the wet sedge meadow community (Barrow-BEO/Barrow-BES) showing much lower reflectance (around 25%), followed by the mesic sedge-grass-herb meadow and wet sedge meadow (Atqasuk) communities (between 24% and 28%, respectively).The other communities showed similar reflectance values in the NIR region.Mesic sedge-grass-herb meadow and wet sedge tundra has a characteristic water absorption feature at ~970 nm, likely from the large amount of standing water (between 30% and 80% coverage; Table 2).

Principal Components Analysis
The PCA using UniSpec data at all four sites showed limited separability between vegetation communities, with 48.9% of the variation in the data being explained by the first component alone, and 70.2% explained by the first two components (Figure S2).At the individual site level and using the three significant regions (blue, red and red edge), the PCA showed much higher community separability (Figure 5).Communities at Atqasuk had clearer separation than at the other sites (Figure 5b).For Barrow-BEO/Barrow-BES, 93.3% of the variation in the data was explained by the first component alone, and 99.7% was explained by the first two components.For Ivotuk, 67.9% of the variation was explained by the first component and 86.5% was explained by the first two components.For Atqasuk, 74.9% and 94.8% of the variation was explained by the first component and the first two components, respectively.
The regions of the spectrum that most contributed to PCA1 at Barrow-BEO/Barrow-BES and Atqasuk were blue, red and red edge; red edge (450-510, 630-690 and 705-745 nm, respectively) contributed most to PCA2 at both Barrow-BEO/Barrow-BES and Atqasuk.At Ivotuk, it was only the blue and red, which contributed most to PCA1 and red edge contributed most to PCA2.   4).4).

Linear Discriminant Analysis
Due to the high variability across the field sites, site-specific classifiers were required as we found one classifier to be inadequate for characterizing all classes across the tower sites (see Supplementary Material; Figure S4).
The individual LDA results for the UniSpec and UniSpecWV2 data across all sites showed minimal overlap in the vegetation communities (Figure 6a,b,d,e,g,h,j,k).This discrimination was achieved with accuracies between 90% and 100% (Kappa coefficients ranging between 0.86 and 1; Table 4, Figure 7) for the UniSpecWV2 data.Adding NDVI, NDWI and EVI further improved the results, resulting in classification accuracies between 90% and 100% (Kappa coefficient values ranging between 0.88 and 1; Table 4, Figure 7).
At Barrow-BEO/Barrow-BES and Ivotuk, the WorldView-2 data resulted in low classification accuracies, ranging between 50% and 64% (Kappa coefficient values of between 0.24 and 0.37; Table 4).Adding the vegetation indices resulted in an equally poor classification, with accuracies ranging between 39% and 67% (Kappa coefficient values between 0.26 and 0.4).Atqasuk, however, had results similar to the UniSpecWV2 data, with classification accuracies ranging between 86% and 92% (Kappa coefficient values between 0.71 and 0.74; Table 4).

Linear Discriminant Analysis
Due to the high variability across the field sites, site-specific classifiers were required as we found one classifier to be inadequate for characterizing all classes across the tower sites (see Supplementary Material; Figure S4).
The individual LDA results for the UniSpec and UniSpecWV2 data across all sites showed minimal overlap in the vegetation communities (Figure 6a,b,d,e,g,h,j,k).This discrimination was achieved with accuracies between 90% and 100% (Kappa coefficients ranging between 0.86 and 1; Table 4, Figure 7) for the UniSpecWV2 data.Adding NDVI, NDWI and EVI further improved the results, resulting in classification accuracies between 90% and 100% (Kappa coefficient values ranging between 0.88 and 1; Table 4, Figure 7).At Barrow-BEO/Barrow-BES and Ivotuk, the WorldView-2 data resulted in low classification accuracies, ranging between 50% and 64% (Kappa coefficient values of between 0.24 and 0.37; Table 4).
Adding the vegetation indices resulted in an equally poor classification, with accuracies ranging between 39% and 67% (Kappa coefficient values between 0.26 and 0.4).Atqasuk, however, had results similar to the UniSpecWV2 data, with classification accuracies ranging between 86% and 92% (Kappa coefficient values between 0.71 and 0.74; Table 4).

Vegetation Map Validation
Higher classification accuracy was achieved when using the UniSpecWV2 reflectance data to classify dry lichen heath, wet sedge meadow (Barrow-BEO/Barrow-BES) and tussock tundra (sandy substrates).Vegetation community classification accuracies between each classifier varied quite considerably.At Barrow-BEO and Barrow-BES, using the LDA function from the UniSpecWV2 reflectance data (including three VIs) improved classification accuracy for both the mesic sedge-

Vegetation Map Validation
Higher classification accuracy was achieved when using the UniSpecWV2 reflectance data to classify dry lichen heath, wet sedge meadow (Barrow-BEO/Barrow-BES) and tussock tundra (sandy substrates).
Vegetation community classification accuracies between each classifier varied quite considerably.At Barrow-BEO and Barrow-BES, using the LDA function from the UniSpecWV2 reflectance data (including three VIs) improved classification accuracy for both the mesic sedge-grass-herb meadow and dry lichen heath, yet decreased the accuracy for the wet sedge Remote Sens. 2016, 8, 978 14 of 24 community (Table 5).The overall classification accuracy for these two sites was 65% (Kappa coefficient 0.43) without VIs and 64% (Kappa coefficient 0.43) including VIs.
At Atqasuk, including VIs in the LDA decreased accuracy for the tussock tundra (sandy substrates) community, but increased it for the wet sedge community (Table 5).The overall classification accuracies for both techniques were 88% (Kappa coefficient 0.71) and 80% (Kappa coefficient 0.56), respectively.
At Ivotuk, an improvement in the classification results occurred when including VIs for both the wet sedge meadow and the mixed shrub-sedge tussock communities, but a decrease in accuracy occurred for the tussock tundra (non-sandy substrates) (Table 5).The overall classification accuracies for both techniques were 67% (Kappa coefficient 0.39) and 73% (Kappa coefficient 0.52), respectively.
The LDA results using the UniSpecWV2 reflectance data showed little variation in vegetation community distribution, with Barrow-BEO and Barrow-BES showing a 92% agreement between with or without inclusion of VIs.Atqasuk showed a little more variation, with 78% agreement using the rescaled data, between either including or excluding VIs.Finally, Ivotuk had similar accuracy using the UniSpecWV2 reflectance data, with 77% agreement in vegetation community distribution.Vegetation community distributions (percentage cover of the tower footprint area) using both analyses techniques are shown in Table 6 and finalized maps are shown in Figure 7.
grass-herb meadow and dry lichen heath, yet decreased the accuracy for the wet sedge community (Table 5).The overall classification accuracy for these two sites was 65% (Kappa coefficient 0.43) without VIs and 64% (Kappa coefficient 0.43) including VIs.
At Atqasuk, including VIs in the LDA decreased accuracy for the tussock tundra (sandy substrates) community, but increased it for the wet sedge community (Table 5).The overall classification accuracies for both techniques were 88% (Kappa coefficient 0.71) and 80% (Kappa coefficient 0.56), respectively.
At Ivotuk, an improvement in the classification results occurred when including VIs for both the wet sedge meadow and the mixed shrub-sedge tussock communities, but a decrease in accuracy occurred for the tussock tundra (non-sandy substrates) (Table 5).The overall classification accuracies for both techniques were 67% (Kappa coefficient 0.39) and 73% (Kappa coefficient 0.52), respectively.
The LDA results using the UniSpecWV2 reflectance data showed little variation in vegetation community distribution, with Barrow-BEO and Barrow-BES showing a 92% agreement between with or without inclusion of VIs.Atqasuk showed a little more variation, with 78% agreement using the rescaled data, between either including or excluding VIs.Finally, Ivotuk had similar accuracy using the UniSpecWV2 reflectance data, with 77% agreement in vegetation community distribution.Vegetation community distributions (percentage cover of the tower footprint area) using both analyses techniques are shown in Table 6 and finalized maps are shown in Figure 7.

Discussion
This study demonstrates that it is possible to successfully distinguish between the dominant vegetation communities at the four Alaska tundra tower field sites using a combination of the blue, red and red edge regions of the electromagnetic spectrum obtained from UniSpec and UniSpecWV2 data.Furthermore, using UniSpec bands within a LDA, successful discrimination between tundra vegetation communities was achieved with classification accuracies between 83% and 98%.These results were used to map tundra vegetation communities within the flux tower footprints using 2 m resolution WorldView-2 imagery.However, we were not able to produce high accuracy (classification accuracies over 80%) vegetation maps for 300 m tower footprints using the same classifier at all four sites due to the sub-pixel heterogeneity of the land cover [76].Finally, scaling from plot (e.g., UniSpec) to landscape (WorldView-2) was challenging because of the difference in spatial resolution between the field spectroscopy narrow band reflectance data and the multispectral satellite imagery.These results indicate that multispectral satellite imagery having a spatial resolution of 2 m is not spatially adequate to represent the innate sub-meter vegetation heterogeneity found in these ecosystems as described by Muster et al. [76].Furthermore, spectral aggregation could also be a factor.For example, even if a pixel contained a homogenous vegetation community, there will inherently be a loss of spectral information when aggregating spectral data across a wider bandwidth.
Using the UniSpec data, the Arctic tundra vegetation communities were separable from a spectral perspective, in particular within the blue, red and red edge regions of the spectrum, which is in agreement with other studies having similar vegetation communities [4,28,39].Although spectral profiles of the primary vegetation community types (wet sedge meadow, tussock tundra, etc.) showed certain diagnostic reflectance and absorption features, it can be challenging to distinguish between community spectral signatures [28,77].In this study, we were not able to distinguish between community spectral signatures using all bands, but obtained successful separation between community spectral signatures in the blue, red and red edge regions.The blue and red regions can be linked to light absorption by plants, while the red edge region (680-740 nm) can be a useful metric of chlorophyll content [78].These regions are a useful metric for the vitality of vegetation, while the short-wave infra-red region can be especially useful in tundra ecosystems, as the reflectance of vegetation such as the dry lichen heath communities are influenced by absorption of cellulose and lignin in this region, as well as reduced absorption by water and chlorophyll [79].This creates a spectral signature that is distinguishable in comparison to green, vascular vegetation, which is shown in our results.Previous studies have also shown high NIR backscatter in the reflectance spectra of shrub tundra communities with decreased slope between beginning and end of NIR reflectance plateau [28].Additionally, both wet sedge and tussock sedge communities in this study show NIR reflectance plateaus with stronger slopes between beginning and end of plateaus.
The mesic sedge-grass-herb meadow and dry lichen heath had a high proportion of low-chlorophyll-rich species, high non-vascular plant cover and higher proportion of dead plant material, which corresponds to lower reflectance in the blue and green regions and higher reflectance in the red region.This is highlighted by Buchhorn et al. [28], who show that areas with higher dead plant material also had limited development of a green reflectance peak due to reduced leaf pigment absorption.Arguably the "greenest" and most productive vegetation, the mixed shrub-sedge tussock community at Ivotuk, had the deepest chlorophyll absorption.Buchhorn et al. [28] also highlight that spectral behavior in the NIR region is intrinsically linked with soil moisture and mean vegetation height.At our study sites, the community with lowest reflectance within the NIR region was the wet sedge meadow at Barrow-BEO/Barrow-BES which had high surface moisture content as well as having tall sedges present.
The classification results in this study may have been affected by mixed signals caused by standing water, different species, soil surface features and moisture content of surrounding vegetation and soil [79].Within certain communities (for example, the wet sedge meadow at Barrow-BEO/Barrow-BES) there was a substantial amount of open water when the measurements were taken, therefore leading to "noise" in the spectral signal.Furthermore, the structure of the tussock sedge communities can also impact the reflectance, as their erectophile structure, along with high density of dead plant material may cause additional backscattering in the NIR region [80].Schaepman-Strub et al. [77] highlight there can be limitations when trying to obtain fractional cover of graminoid and shrub-dominated communities within peatlands.The fine-scale differences between the plant traits, community dynamics and spectral properties create the unique spectral signatures that allow each community to be differentiated.This can be especially useful as many areas of the Arctic are undergoing an increase in shrub-cover [13], allowing for potential mapping of this expansion, as well as partitioning vegetation community contribution to landscape level carbon budgets.
Adding VIs improved classification accuracy as they reduce spectral noise and can be directly related to plant physiology [45], and therefore linked to detailed arctic tundra vegetation properties [81].VIs have shown that the simple transformation of band reflectance can be less sensitive to external variables [45], subsequently producing potentially more accurate classifications for heterogeneous landscapes over large areas.The improved performance when including VIs in this study is especially useful in regards to monitoring plant productivity and ecosystem functioning.It has been shown previously [82,83] that NDVI can be linked heavily with net ecosystem exchange (NEE) [84,85] and CH 4 fluxes [25] for example.NDVI can be highly correlated with vegetation leaf area index (LAI), but also is invariant at higher LAI ranges and can be very sensitive to background variations in vegetation communities such as background soil and water [8].The inclusion of the EVI and NDWI in our analysis improved classification accuracies across the majority of sites, as they allow for more detailed analysis of canopy structure variation [64] and plant water stress/soil water content [37], respectively.
Higher accuracy classifications (between 92% and 96%) occurred from the UniSpecWV2 data.These results suggest that it is possible for remotely-sensed tundra vegetation to be scaled from plot (<1 m) to the patch (~500 m) scale due to the reduced noise in the signal found in the extracted data, allowing for better separation between communities.The UniSpecWV2 reflectance data was used to produce the maps of the Arctic tundra vegetation with and without inclusion of VIs (classification accuracies for Barrow-BEO/Barrow-BES, 65% and 64%; Atqasuk, 88% and 80%; and Ivotuk, 67% and 73%, respectively).Although the classifier used in the analysis (with input UniSpecWV2 data) was obtained using narrow band field spectroscopy data rescaled to match the broad bands of the satellite imagery, they still maintain a level of accuracy that is found when using the narrow bands (>80% accuracy).However, at the individual field sites, classification accuracies were lower, possibly due to mixed signals caused by the heterogeneous nature of the vegetation communities, with spectral characteristics in each reflecting a combination of different vegetation communities within the 2 m resolution of each pixel.Where individual community separation was less successful, this could be due to the heterogeneous nature of the land cover.For example, at Ivotuk, the wet sedge meadow contains some areas having substantial shrub cover, causing a mixing of the spectral signature [78], with both communities showing similar spectral signatures, especially in the red edge region.Kushida et al. [86] showed no significant difference between sedge and shrub plots located at an Alaskan study site, even though both communities have very different structural properties.At first glance, many of the tundra communities evaluated in this study appeared similar, containing similar species.However, at the spectral level they were significantly different.For example, the mesic sedge-grass-herb meadow and dry lichen heath tundra communities at Barrow-BEO/Barrow-BES and the tussock tundra communities at both Atqasuk and Ivotuk could be considered physiognomically similar respectively, due to the communities containing similar species/abundance, yet they differed significantly when evaluating the spectral profiles.This could be due to the proportions of non-vascular plant species in the dry lichen heath community in comparison to the mesic-sedge-grass-herb meadow, which have an impact on the spectral signal [3].The spectral reflectance of non-vascular communities is influenced by absorption of cellulose and lignin (and lack of absorption by chlorophyll as found in "green" vascular plant communities [79].
The results of the LDA using the WorldView-2 data indicate a loss in classification accuracy.This again highlights the "ecologically complex" nature of land cover across these field sites, with the 2 m satellite imagery resolution being insufficient to capture the fine-scale heterogeneity in vegetation cover.However, Atqasuk showed similar classification accuracies between using the UniSpecWV2 and WorldView-2 data.This may be due to only having two dominant communities at this site, therefore the landscape is much less heterogeneous, with each community being spectrally different from one another.Vegetation across these landscapes ranged from short stature communities in the high Arctic, to communities containing larger shrubs and tussock sedges as one moves further south and with varying levels of moisture, all linked to further ecosystem processes.This may have resulted in the inability to produce an all-site classifier, due to the lack of representation of all vegetation communities at each site.However, the vegetation communities described in this study are representative of large areas across the Alaskan Arctic.The spatial heterogeneity of vegetation communities across these landscapes means they may contribute differently to the landscape-scale carbon exchange [77,87], highlighting the necessity to be able to map vegetation communities at a fine-scale.Therefore, the classifiers created are valuable with regards to their ability to map fine-scale vegetation communities.
One limitation to the research presented here is there is no temporal aspect to the vegetation maps, as they are created using data collected at one time period during 2014 (satellite imagery for Ivotuk was obtained in 2013 due to no suitable imagery available for 2014).The microclimatic differences in climate (from year to year), hydrology and topography can influence vegetation community composition and biomass/greenness [12,88,89].Therefore, data should be collected within a suitable time period (ideally within the same year/month) during the growing season so it is as representative as it can be of the community present [89].However, changes in Arctic vegetation community composition occurs at a slower rate than elsewhere due to low temperatures and short growing seasons, and the vegetation communities found in these locations could be considered "stable", with the distribution of each community unlikely to change rapidly enough to make the findings of this study non-applicable.Even if the communities are not stable [88,90], the growth rate is so slow that it still makes the study applicable across a short time span like the difference in data collection dates at Ivotuk.Potential expansion of wet sedge meadow communities could occur with further permafrost thaw with areas becoming wetter [91,92] over a decadal time scale [93], but we believe the vegetation maps created within this study were likely not majorly affected by the temporal variation across these sites, making them applicable for using within carbon cycle prediction models.Given the substantial cloud cover of these arctic ecosystems, previous studies (including Langford et al. [47] have used imagery and ground truthing data collected in different years.Ideally the use of airborne hyperspectral imagery [19] or unmanned aerial vehicles (UAV) with very high spatial resolution (sub-centimeter for UAV), collected at the same time of ground data will improve the accuracy of vegetation mapping, especially at the patch scale [94].In addition, airborne LiDAR [47] might also be used to improve classification accuracies by providing data of vegetation structure and biomass profiles, in addition to terrain features that may influence landscape temperature and moisture conditions and consequently vegetation characteristics.

Conclusions
We show that field spectroscopy can be used to successfully discriminate and map dominant tundra vegetation communities for four sites in northern Alaska.The vegetation communities found at Ivotuk represents the dominant tundra vegetation type in Alaska, while communities at Barrow and Atqasuk represent approximately 60% of all Arctic wetlands [23].In this study, these vegetation communities were classified with patch scale (~500 m) accuracies ranging from over 50% to 80%.This highlights the applicability and representativeness of this methodology for mapping these heterogeneous environments at this scale.There are however, limitations to scaling up and mapping tundra vegetation communities because of the requirements for high spatial resolution satellite or airborne data across the landscape.This study highlights the benefits of using fine-scale spectroscopy over coarse areas (300 m eddy covariance tower footprints) for tundra vegetation classification and demonstrates the ability to map vegetation communities using field spectroscopy in combination with high spatial resolution multispectral satellite imagery.Although coarse scale vegetation mapping across Arctic tundra ecosystems is suitable for broad-scale vegetation distribution monitoring, there is still a necessity to include fine-scale vegetation community mapping in further carbon cycle models in order to accurately assess their influence, thus facilitating the monitoring of changes occurring as a result of a warming Arctic and understanding large scale ecological processes in Arctic plant community assemblage.

Figure 1 .
Figure 1.Map showing locations of the tower field sites included in this study (upper left) and corresponding 300 m diameter eddy covariance tower footprints (lower left) for: (a) Barrow-BEO; (b) Barrow-BES; (c) Atqasuk; and (d) Ivotuk site images (right) showing land cover types and topography.The lower left panel shows overlying WorldView-2 satellite imagery for each site obtained on 25 July 2014 (Barrow-BEO/Barrow-BES), 9 July 2014 (Atqasuk) and 21 June 2013 (Ivotuk).

Figure 1 .
Figure 1.Map showing locations of the tower field sites included in this study (upper left) and corresponding 300 m diameter eddy covariance tower footprints (lower left) for: (a) Barrow-BEO; (b) Barrow-BES; (c) Atqasuk; and (d) Ivotuk site images (right) showing land cover types and topography.The lower left panel shows overlying WorldView-2 satellite imagery for each site obtained on 25 July 2014 (Barrow-BEO/Barrow-BES), 9 July 2014 (Atqasuk) and 21 June 2013 (Ivotuk).

Figure 3 .
Figure 3. Spectral profiles for all vegetation communities.Inset photographs show the dominant vegetation community type.The blank lines show the mean value of the spectral signatures (individual spectral samples are shown in grey).The sample locations at Barrow include: (a) Dry lichen heath; (b) Mesic sedge-grass-herb meadow; and (c) wet sedge tundra.The sample locations at Atqasuk include: (d) tussock tundra (sandy substrates); and (e) wet sedge meadow.Samples at Ivotuk include: (f) tussock tundra (non-sandy substrates); (g) mixed shrub-sedge tussock; and (h) wet sedge meadow.

Figure 3 .Figure 3 .
Figure 3. Spectral profiles for all vegetation communities.Inset photographs show the dominant vegetation community type.The blank lines show the mean value of the spectral signatures (individual spectral samples are shown in grey).The sample locations at Barrow include: (a) Dry lichen heath; (b) Mesic sedge-grass-herb meadow; and (c) wet sedge tundra.The sample locations at Atqasuk include: (d) tussock tundra (sandy substrates); and (e) wet sedge meadow.Samples at Ivotuk include: (f) tussock tundra (non-sandy substrates); (g) mixed shrub-sedge tussock; and (h) wet sedge meadow.

Table 1 .
Acquisition dates for both the UniSpec DC field spectroscopy measurements and the corresponding WorldView-2 satellite imagery.Due to lack of suitable cloud-free imagery available for Ivotuk during 2014, the closest cloud-free acquisition date from 2013 was used instead.

Table 2 .
Frequency of occurrence and abundance for each plant functional type (PFT) or tundra vegetation category across the four Alaska flux tower sites.

Table 3 .
Analysis of Variance between each vegetation community and WorldView-2 multispectral imagery bands.

Table 4 .
Classification accuracies and Kappa coefficients of the results from the Linear Discriminant Analysis (LDA) for the UniSpec, UniSpecWV2 and WorldView-2 data when including or excluding NDVI, NDWI and EVI.

Table 4 .
Classification accuracies and Kappa coefficients of the results from the Linear Discriminant Analysis (LDA) for the UniSpec, UniSpecWV2 and WorldView-2 data when including or excluding NDVI, NDWI and EVI.

Table 5 .
Error matrix and associated accuracy for LDA classifiers for all four flux tower sites.The left hand tables are the LDA classifier without inclusion of VIs and the right hand tables are the LDA classifier with VIs.VIs are vegetation indices; normalized difference vegetation index (NDVI), normalized difference water index (NDWI) and enhanced vegetation index (EVI).

Table 6 .
Vegetation community distributions derived from the Linear Discriminant Analysis and the Linear Discriminant Analysis including NDVI, NDWI and EVI vegetation indices.