Next Article in Journal
Evaluation of Synoptic Snowfall on the Antarctic Ice Sheet Based on CloudSat, In-Situ Observations and Atmospheric Reanalysis Datasets
Next Article in Special Issue
InSAR Time Series Analysis of L-Band Data for Understanding Tropical Peatland Degradation and Restoration
Previous Article in Journal
Coupling Hyperspectral Remote Sensing Data with a Crop Model to Study Winter Wheat Water Demand
Previous Article in Special Issue
On the Potential of Sentinel-1 for High Resolution Monitoring of Water Table Dynamics in Grasslands on Organic Soils
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characterizing Boreal Peatland Plant Composition and Species Diversity with Hyperspectral Remote Sensing

by
Mara Y. McPartland
1,2,*,
Michael J. Falkowski
3,
Jason R. Reinhardt
2,
Evan S. Kane
4,5,
Randy Kolka
6,
Merritt R. Turetsky
7,
Thomas A. Douglas
8,
John Anderson
9,
Jarrod D. Edwards
9,
Brian Palik
6 and
Rebecca A. Montgomery
2
1
Department of Geography, Environment & Society, University of Minnesota, Minneapolis, MN 55455, USA
2
Department of Forest Resources, University of Minnesota, Saint Paul, MN 55108, USA
3
Natural Resource Ecology Laboratory, Colorado State University, Fort Collins, CO 80523, USA
4
School of Forest Resources and Environmental Science, Michigan Technological University, Houghton, MI 49931, USA
5
Northern Research Station, USDA Forest Service, Houghton, MI 49931, USA
6
Northern Research Station, USDA Forest Service, Grand Rapids, MN 55744, USA
7
Department of Integrative Biology, University of Guelph, Guelph, ON N1G 2W1, Canada
8
Cold Regions Research and Engineering Laboratory, U.S. Army, Fort Wainwright, AK 99703, USA
9
Geospatial Research Lab, U.S. Army Corps of Engineers, Alexandria, VA 22315, USA
*
Author to whom correspondence should be addressed.
Remote Sens. 2019, 11(14), 1685; https://doi.org/10.3390/rs11141685
Submission received: 7 June 2019 / Revised: 3 July 2019 / Accepted: 5 July 2019 / Published: 16 July 2019
(This article belongs to the Special Issue Remote Sensing of Peatlands II)

Abstract

:
Peatlands, which account for approximately 15% of land surface across the arctic and boreal regions of the globe, are experiencing a range of ecological impacts as a result of climate change. Factors that include altered hydrology resulting from drought and permafrost thaw, rising temperatures, and elevated levels of atmospheric carbon dioxide have been shown to cause plant community compositional changes. Shifts in plant composition affect the productivity, species diversity, and carbon cycling of peatlands. We used hyperspectral remote sensing to characterize the response of boreal peatland plant composition and species diversity to warming, hydrologic change, and elevated CO2. Hyperspectral remote sensing techniques offer the ability to complete landscape-scale analyses of ecological responses to climate disturbance when paired with plot-level measurements that link ecosystem biophysical properties with spectral reflectance signatures. Working within two large ecosystem manipulation experiments, we examined climate controls on composition and diversity in two types of common boreal peatlands: a nutrient rich fen located at the Alaska Peatland Experiment (APEX) in central Alaska, and an ombrotrophic bog located in northern Minnesota at the Spruce and Peatland Responses Under Changing Environments (SPRUCE) experiment. We found a strong effect of plant functional cover on spectral reflectance characteristics. We also found a positive relationship between species diversity and spectral variation at the APEX field site, which is consistent with other recently published findings. Based on the results of our field study, we performed a supervised land cover classification analysis on an aerial hyperspectral dataset to map peatland plant functional types (PFTs) across an area encompassing a range of different plant communities. Our results underscore recent advances in the application of remote sensing measurements to ecological research, particularly in far northern ecosystems.

Graphical Abstract

1. Introduction

The far northern regions of the globe are warming at twice the global average, resulting in significant effects on the vegetation composition of arctic tundra and boreal ecosystems. Prominent among the changes to plant communities include the pan-arctic northward expansion of shrubs into tundra, and the spread of wetland vegetation as a result of permafrost thaw and ground subsidence [1,2,3,4]. Boreal forest ecosystems are also shifting from dominance by conifers to deciduous species as a result of increased fire prevalence [5]. The relationship between climate and plant community composition has important implications for biogeochemical cycles. For example, research indicates that elevated temperatures in arctic environments increase the rate of nitrogen mineralization, promoting the encroachment of woody plants into the tundra [6]. Increased shrub cover leads to an increase in the rate of carbon uptake, but with uncertain consequences for the overall arctic carbon cycle [4,7,8]. Shifts in arctic and boreal plant community composition and biogeochemical function can have significant impacts on wildlife [9,10] and the livelihoods of local communities [11,12]. Despite these wide-ranging impacts, resource and logistical challenges present an obstacle to the study of remote ecosystems. These challenges make the examination of factors such as plant community compositional change difficult. Indeed, given the extent and remoteness of northern ecosystems, remote sensing approaches to assessing vegetation change are necessary for understanding the ecological effects of climate change over large spatial scales.
Peatlands cover approximately 15% of land cover in the northern hemisphere and factor significantly in the global carbon cycle [13,14,15]. Peatlands are highly heterogeneous ecosystems characterized by thick organic soils and dominated by a variety of plant functional types (PFTs) depending on landscape position and hydrology [16,17]. PFTs refer to the physiological, morphological, structural, and chemical attributes that determine how plants allocate resources, take up carbon, and reproduce themselves [18,19]. In this study, we use the PFT framework as a trait-based method of vascular plant-species classification based on shared structure and function [18]. Changes in mean annual temperature and hydrologic conditions in peatlands can lead to shifts in the composition of PFTs, ultimately affecting ecosystem productivity [20,21,22,23]. Prior research in peatlands has indicated that warming and drying affect PFT composition by promoting dominance by woody shrubs [23,24,25]. Conversely, raised water tables from increased precipitation or permafrost thaw lead to an increase in graminoid species such as sedges [20,21,22]. Lowered water tables and increased shrub dominance have been associated with decreased annual gross primary productivity (GPP) relative to elevated water tables [22,26]. Studies have also demonstrated a link between the graminoid cover and methane emissions rates from wetland ecosystems [27,28]. Given their importance in regulating ecosystem function, mapping the distribution of PFTs could be an important element of predicting ecosystem responses to environmental change [29,30].
High spatial resolution hyperspectral imagery is emerging as a promising method for conducting remote assessments of PFT distribution [19,31]. Hyperspectral data leverages a far greater number of spectral bands as compared to multi-spectral data, ultimately allowing for the differentiation of vegetation characteristics with a high degree of fidelity [32]. For example, prior research has demonstrated that variation in the structural, physiological, and chemical characteristics of particular species or PFTs interact uniquely with different regions of the electromagnetic spectrum, allowing them to be distinguished from other species [33,34]. Species-specific properties that affect reflectance include unique leaf pigments such as chlorophyll and anthocyanin [35,36,37,38], leaf nitrogen content [39,40], and leaf water content [41]. Recent research has shown that phylogenetic distance both within and among species is correlated with spectral dissimilarity, a relationship that scales from leaf-level to plot-level canopy observations [42] (Schweiger et al. 2018). Leveraging this type of information could allow for species to be mapped over large areas through classification analysis of hyperspectral data, improving our ability to track changes in species’ distributions and diversity across broad spatial extents. Although mapping PFTs rather than individual species may not allow for accurate estimates of species diversity, it still affords valuable information regarding the structure and function of ecosystems as a whole [29]. The use of hyperspectral remote sensing as a tool for mapping PFTs could provide ecologists with the ability to track ecosystem responses to climate change [43]. Furthermore, these types of analyses could ultimately help inform management strategies and environmental policy approaches focused on ecosystem conservation and management in the face of climate change.
In addition to tracking compositional changes using remote sensing, a number of approaches to characterizing species and functional diversity using hyperspectral methods have also been investigated. In low diversity environments such as arctic and boreal ecosystems, the reduction of species richness may have amplified effects on ecosystem productivity due to low functional redundancy in plant types [44,45,46,47]. For example, peatland studies have found a loss of species in arctic systems in response to the warming effects of carbon cycling rates [45,48]. Thus, tracking changes in the vascular functional diversity of boreal peatlands is important for understanding the potential impacts of climate change on the function of these ecosystems [19]. Recent research using hyperspectral methods has linked spectral variation among pixels with species diversity by treating a stack of pixels as a reflectance signature that captures the optical characteristics of vegetation across hundreds of spectral bands [49,50]. Remote sensing estimation of diversity has previously been accomplished by linking spectral heterogeneity among pixels with species richness, a method based on the spectral variation hypothesis that unique species will present with strongly differentiated spectral properties [42,51,52,53,54]. Several different statistical approaches have been used in the calculation of spectral heterogeneity. A simple method of estimating spectral heterogeneity from a multi- or hyperspectral dataset is comparing the coefficient of variation of multiple measurements [55,56,57,58]. Another approach is calculating Euclidian distances among different scans or pixels derived from multivariate eigenvector methods such as principal component analysis [52]. The further a pixel lands from the centroid of the data, the greater the spectral heterogeneity of that pixel [51,52]. This approach has been correlated with alpha diversity metrics that account for both species richness and evenness, such as the Shannon index, as well as with metrics such as functional diversity or the number of different plant types present in a pixel [42,59]. Together, these studies indicate the ability of multiple different remote sensing approaches to characterize the ecological diversity of a landscape. Ultimately, remotely sensed estimates of diversity are integral to mapping essential biodiversity variables (EBVs) – structural and functional characteristics of ecosystems that predict biodiversity levels at a global scale [60,61,62].
The objective of this research was to assess the ability of remote sensing approaches for characterizing and tracking changes in the plant communities of peatland ecosystems in response to global change drivers. We collected the vegetation and hyperspectral data at two different peatland ecosystems characteristic of boreal and arctic regions across the northern hemisphere. Both sites are located within large-scale ecosystem manipulation experiments examining the effects of a suite of climate change drivers on peatland ecosystems. We analyzed our data to determine which elements of vascular plant functional composition and structure influence the reflectance spectra. We also examined relationships between plot-level species diversity and spectral variation using multiple methodological approaches. Predicated on the relationships that we demonstrated between near-earth spectral reflectance and plant functional cover, we leveraged aerial hyperspectral imagery to map the distribution of PFTs across a boreal peatland in interior Alaska, USA.

2. Methods

2.1. Study Sites

We collected data at two sites representing typical northern boreal and arctic peatland ecosystems. The first site was at the Alaska Peatland Experiment (APEX) at the Bonanza Creek Long Term Ecological Research Station, located 30 miles west of Fairbanks, Alaska (64.82°N, 147.87°W) (Figure 1). The site is located in a sedge- and shrub-dominated rich fen on the floodplain of the Tanana River. APEX was initiated in 2005 as a long-term manipulation to study the impact of altered hydrology on peatland ecology and biogeochemistry. APEX consists of three water table manipulation treatments: a raised, a lowered, and a control treatment, each 120 m2 in size (Figure 2, left panel). The lowered treatment plot is drained by a trench that borders the plot. Water from the trench is then pumped into the raised water table treatment. The water table in the control plot is not manipulated. We performed field data collection within the APEX water table manipulation treatment plots.
The second study site was at the Spruce and Peatland Response Under Changing Environments (SPRUCE) project located at the Marcell Experimental Forest in northern Minnesota (47.30°N, 93.29°W) (Figure 1). SPRUCE was initiated in 2014 as a climate change experiment to explore the effects of increased temperatures and elevated CO2 levels on ecological and biogeochemical processes in peatlands at the southern edge of the boreal region (Figure 2, right panel). The experiment is located within an ombrotophic bog dominated by Sphagnum mosses and black spruce (Picea mariana) with an understory of ericaceous shrubs, graminoids, and forbs. The SPRUCE experiment is a regression-based model in which ecological changes are measured in response to a broad range of temperatures [63]. Temperatures and CO2 levels are manipulated in 10 chambered treatment plots (Figure 2, right panel). An additional three unenclosed control plots were also tracked in this experiment. The chambers consist of octagonal enclosures 12.8 m in diameter and 7 m tall. The enclosures penetrate completely through the peat soils down to mineral material below. Peat soils average 2.5 m in depth across the site. Ambient temperatures and CO2 levels are maintained in one control chamber while another control chamber has ambient temperatures and elevated CO2. The remaining eight treatment chambers are maintained at +2.25, +4.5 +6.75, and +9 °C above ambient temperatures to simulate a range of possible warmer climate scenarios. Carbon dioxide levels are maintained at ambient levels in half of the temperature treatment chambers and elevated to between 800 and 900 parts per million in the other half of the chambers. Each chamber, as well as the unchambered control enclosures, contain three 2-m2 vegetation plots, where sampling was conducted for this experiment.

2.2. Data Collection

2.2.1. Vegetation Cover Sampling

Species and PFT cover were measured differently at the two study sites owing to differences in their respective data collection protocols (Table 1). In particular, non-vascular plants at APEX (which accounted for less than 10% of total cover) were recorded where present, whereas, at SPRUCE, mosses carpeted the soil surface and were therefore considered to have 100% cover across the site and were not included in the sampling protocol. At APEX, species cover was recorded using a point bar vegetation survey method [64]. Species were recorded along two perpendicular transects positioned across a series of replicate plots in each treatment. A metal frame was positioned above the canopy 1 m above the surface of the peat. A laser pointer was inserted into 10 evenly spaced holes in the metal frame along the transect. Species and height above the ground were recorded for all interceptions of the laser beam. This was repeated for each “hit” from the canopy to the surface of the peat. If the target was not photosynthetic it was recorded as either “standing dead”, indicating that it was non-photosynthetic biomass but still standing, or as “litter”, indicating it was part of the thatch forming at the surface of the peat. For the lowered and raised treatments, four replicate plots were sampled, and in the control treatment five replicates were sampled. Each species was categorized by PFT, and percent cover for each type was calculated based on the percentage of hits of each PFT relative to the entire sample (see Table 2 for species list). In addition, we calculated the leaf area index (LAI) for each point along each transect as the number of living layers of vegetation present between the canopy and the soil surface [25].
Vegetation cover at SPRUCE was sampled in mid-July across the range of temperature and carbon dioxide treatments. Three 2-m2 plots were sampled in each of the 13 chambers for a total of 33 plots in the final analysis. Each plot was sampled using a 1- x 2-m frame that was set onto PVC pipes installed in the experiment to facilitate repeated data collection in the exact same area. The frame was divided into 50 cells 20 x 20 cm in size. The presence of all species was recorded within each cell of the grid, and species were later grouped into PFTs (Table 2). If a species or PFT was present in half of the cells it was considered to occupy 50% of the plot. Our sampling method at SPRUCE did not allow for the calculation of the LAI.

2.2.2. Spectral Sampling

To link spectral reflectance with functional composition and species diversity, sampling of spectral reflectance was performed at APEX and SPRUCE in the same plots where vegetation cover data were collected. Spectral reflectance measurements were taken at APEX using an Analytical Spectral Devices Fieldspec Pro that measured reflectance in 1-nm bandwidths between 300 and 2500 nm (Figure 3). Data were collected during peak growing season on 29 June 2016. Scans were performed during a 1-h window on each side of solar noon, which occurred at approximately 12:30 p.m. local time. The sun at the time of data collection was at an azimuth angle of approximately 151° east of north. Three scans were performed at each plot at the same location to acquire a local average. The instrument was positioned 1 m above the peat surface, and captured a field of view of 44 cm in diameter (Equation (1)) [65].
G F O V = 2 * h e i g h t * t a n ( 25 ° ) 2
Data were processed such that water and atmospheric absorption bands that had a relative reflectance of greater than one or less than zero were automatically excluded, as well as areas that exhibited noise from the instrument. This resulted in three regions of the electromagnetic spectrum being included in the analysis (Figure 3). Thirteen plots were sampled for vegetation and spectral reflectance among the three APEX water table manipulation treatments: five plots in the control treatment and four plots each of the lowered and raised plots.
Spectral data were collected at SPRUCE on 22 September 2016 under clear sky conditions beginning at 1:00 p.m. local time. The sun at the time of data collection was at an azimuth angle of 178° east. Although it was past peak growing season, senescence had only just begun at the site, and experimentally warmed plots were still near peak growing conditions. Data collection was performed using a PP Systems UniSpec-DC spectroradiometer that detected incoming and reflected solar radiation in 3–4 nm bandwidths between 310 and 1100 nm (Figure 3). Three scans were performed above each vegetation plot in the same sampling location to acquire a local average. Details on data collection and processing are provided in McPartland et al. 2018 [25]. Reflectance was calculated for the average of the three scans following the methods of Harris et al. (2014) [16] and Wang et al. (2016 & 2018) [49,50]. Bands at the beginning and end of the spectrum were eliminated due to noise, resulting in a data range of 400–1000 nm used in the analysis (Figure 3). Data for SPRUCE chamber number 4 (+4.5 °C, elevated CO2) were also excluded because reflectance data for that plot presented as a significant outlier, with values that were far outside of the typical reflectance distribution collected at the other plots. This was likely due to changing sky conditions between the time that the white reference scan was taken, and the time of measurement.

2.2.3. Aerial Hyperspectral Data Collection

We used a hyperspectral dataset collected in July of 2014 at the Bonanza Creek Long Term Ecological Research Station, where the APEX field site is located. (Figure 2, right panel). The aerial dataset included the APEX study site and surrounding peatland ecosystems (Figure 4). The aerial image encompassed a variety of peatland cover types including black spruce and tamarack bogs, sedge meadows, tussock grasslands, and birch and willow shrub communities. Airborne hyperspectral data were collected by the remote sensing company SpecTir. The instrument used in the data collection was a ProSpecTir VNIR SWIR dual sensor push broom imaging spectrometer that collected data at 360 channels from approximately 400 to 2450 nm at a spectral resolution of approximately 2.9 nm in the visible and near-infrared, and 8.5 nm in the short-wave infrared portions of the electromagnetic spectrum. The imagery was collected approximately 200 m above ground surface, resulting in a spatial resolution of 1 m2 [66]. The study area is approximately 2 km2 in size. Data were processed using the Spectral Angle Mapper algorithm in the Python-based Spectral Classification Plug-in [67]. Ground control was performed by SpecTir and participants in this study. More detailed information on the geolocation and post-processing of the aerial hyperspectral data is provided in Anderson et al. (in review) [68].

2.3. Data Analysis

2.3.1. Analysis of Vegetation Composition and Species Diversity

For both sites, we calculated the percent cover of each PFT as a fraction of the total observations made at the plot level. The total number of plots included in our analysis was 13 across the three treatments at APEX, and 33 plots across 13 chambers at SPRUCE. At APEX, the effect of water table treatment on PFT cover was analyzed using a one-way ANOVA with a post-hoc pairwise means comparison to determine whether significant differences existed among treatments. At SPRUCE, we used a linear mixed effects modeling framework to evaluate the response of PFT cover to warming and elevated CO2. Temperature was added to the model as a continuous variable, and CO2 level was included as a factor. We included treatment chamber in the model as a random effect and evaluated model fit using the Akaike information criterion, as well as both marginal and conditional R2, which describe the variance explained by the fixed effects and full model, respectively [69]. Linear mixed effects modeling was done using the lme4 package in R.
We calculated species diversity using the Shannon diversity index which incorporates both species richness and evenness into a single measure of alpha diversity:
H =   i = 1 s p i log ( p i )
In this equation, Shannon diversity (H’) is calculated at the plot-level in which pi represents the proportion of the population represented by species i [70]. Shannon diversity was calculated for each plot using the R Vegan package for community ecology [71]. We examined whether significant differences in species diversity among plots existed using a one-way ANOVA with a post-hoc pairwise means comparison at APEX, and linear mixed effects modeling at SPRUCE. We assessed significance based on an alpha level of 0.05 in the ANOVA, and evaluated model coefficients in the linear mixed effects model.

2.3.2. Spectral Data Analysis

For both sites, we performed principal component analyses (PCA) on the full reflectance spectra to determine which PFTs contributed most strongly to variation among spectral scans (Figures S1 and S2) [72,73]. This analysis reduced the dimensionality of the full dataset while preserving variation among plots within each experiment [74,75]. We were then able to determine the strength of the relationship between spectral variation and PFT cover. We performed the analyses using the R vegan package [71], and used permutation tests to determine the significance of different cover types on the spectral PCA results. This modeled the relationship between the PCA axes of the spectral data and PFT cover as a linear relationship. By this method, significance values are assigned by automatically comparing the sample data with a randomly permuted dataset. We considered different PFTs to be significant drivers of spectral variation at a p-value of less than or equal to 0.05.
We also examined the relationship between Shannon diversity and spectral variation through several measures of community and spectral heterogeneity following similar methods to those used by Wang et al. (2016) [49], who found a positive relationship between plot-level plant species diversity and the coefficient of variation of the associated spectral scans. Based on their findings, we predicted that greater coefficients of variation calculated for our average per-plot spectra would correlate with higher Shannon diversity index values [49]. We also performed a similar analysis based on research done by Rocchini et al. (2010) [53] that theorized that greater community heterogeneity would be correlated with greater spectral heterogeneity. In this analysis, we used measured heterogeneity as centroid distances derived from a principal component analysis where heterogeneity was measured as the distance from the center point of the PCA [53,76]. Basic linear regression was used to assess the relationship between diversity and spectral coefficients of variation, as well as between spectral heterogeneity and community heterogeneity measured using Euclidian distances from the global mean. Histograms were used to determine whether the data were normally distributed, and we used Anderson-Darling tests to reject the null hypothesis of uniformity in the distribution of our data (Figure S3).

2.3.3. Hyperspectral Image Analysis

For the APEX field site, we used the aerial imagery to perform a supervised classification analysis using a random forests (RF) model to classify PFTs across the landscape [77]. We designated four different classes representative of dominant site PFTs. The classes were coniferous forests, low shrub, tussock meadow, and graminoid fen (Figure S4). Before performing the analysis, we manually identified approximately one hundred training points within each class using a high-resolution aerial photograph collected in July of 2016 using an eBee senseFly unmanned aerial system (UAS) mounted with a Canon S110 camera (Figure 4, left panel). Although identifying training points in the field would have been ideal, we are confident that we were able to select training points that accurately reflected the ascribed vegetation cover class based on the image and our knowledge of the site. The UAS was flown at a height of approximately 92 m above the ground. The resulting imagery had a 3-cm ground resolution and was mosaicked together and georeferenced using a series of high-resolution GPS points that were collected via total station at the time of data collection. The UAS image was mosaicked using the Agisoft Photscan software platform. RF analyses were completed using the randomForest R package [78] and data were processed using the rgdal R package [79]. Standard model parameters were used in the RF model. Five hundred trees were created, and the number of variables considered at each node was the square root of the number of parameters. We used a model selection approach, as implemented in the rfUtilities R package [80], to find the most parsimonious and explanatory set of predictors. We used the overall out-of-bag accuracy and Cohen’s kappa coefficient to evaluate model performance [81]. Variable importance was assessed using the mean decrease in accuracy and mean decrease in impurity upon iterated variable permutation, as implemented in the randomForest R package; variables that are more important will tend to cause higher mean decreases in both accuracy and node impurity when they permuted [77]. All associated mapmaking was performed using Quantum GIS [82].

3. Results

3.1. Response of Plant Functional Composition and Species Diversity to Experimental Maniuplation

Analysis of the vegetation cover data indicated that plant communities at both sites had diverged as a result of the experimental manipulation of water tables, temperatures, and CO2 levels. At APEX, the most striking trend was the overall higher percent cover of litter and reduced cover of nearly all other PFTs in the lowered water table treatment as compared with the control and raised water table treatments (Table 3, Figure 5). Sedges, equisetum, forbs, and shrubs all showed trends towards higher abundance in the control treatment, and mosses were significantly more abundant in the control than other treatment (Table 3). At SPRUCE, we found a trend towards greater shrub cover with warming, coinciding with a decrease in forb cover (Figure 6, Table 4). Graminoid species remained fairly consistent across the treatments (Figure 6).
Our analysis of diversity also showed a trend towards greater alpha diversity in the control treatment as compared to the raised or lowered treatments at APEX (ANOVADiversity F2,10 = 3.52, p = 0.07) (Figure 7). We did not find any trends in plot-level diversity with warming or CO2 level at the SPRUCE site (Table 3).

3.2. Relationship between Community Composition and Spectral Response

The results of the principal component analyses (PCAs) of our plot-level spectral datasets indicated that variation in plot-level vegetation characteristics was related to differentiation in spectral reflectance at both sites (Figure 8). At APEX, results from our permutation tests indicated that spectral reflectance varied by treatment (r2 = 0.45, p = 0.005), with the most variation in reflectance explained by plot-level species diversity (Table 5). Variation in the spectra at APEX was closely associated with the percent cover of mosses and forbs. Mosses and forbs together accounted for less than 10% on average of the cover at APEX, but had a strong effect on spectral reflectance, particularly in the control and raised water table treatment plots. The lowered water table treatment was associated with a greater litter cover, coinciding with lower plot-level species diversity (Figure 8, left panel). Overall, the first principal component explained 79% of the variation in the dataset, driven primarily by greater diversity in the control and raised water table treatments, and higher litter cover in the lowered water table treatment (Figure 8, left panel, Table 5).
PCA results from the reflectance data at SPRUCE showed significant spectral differences across temperature (r2 = 0.28, p = 0.01) and CO2 treatments (r2 = 0.19, p = 0.01) (Figure 8, right panel). Like APEX, reflectance spectra responded to differences in forb cover (Table 6). Unlike APEX, shrub cover also emerged as a significant vector in the SPRUCE data. Overall, the first principal component explained 96% of total variance among spectral scans. Since we did not find any strong differences in diversity among treatments, diversity did not emerge as a significant vector in the analysis. These results suggest that the relative cover of shrubs versus forbs associated with temperature and CO2 level is the main driver of spectral differences among treatments at SPRUCE.

3.3. Relationship between Species Diversity and Spectral Variation

We found that plot-level species diversity was a significant predictor of spectral variation among plots at the APEX site, but not at the SPRUCE site (Table 5 and Table 6). Based on this result we examined the relationship between species diversity and spectral diversity at the APEX site. The relationship between Shannon diversity and spectral CV showed a positive relationship (R2adj = 0.36, F2,11 = 4.23, p = 0.06, RMSE = 0.30) (Figure 9, left panel). The relationship between community and spectral heterogeneity showed a positive relationship (R2adj = 0.21, F2.11 = 4.11, p = 0.07, RMSE = 0.78) (Figure 9, right panel). It should be noted that differences in heterogeneity are relative to a global mean value. Points closest to zero along both axes are less heterogeneous than those further from the global mean.

3.4. Hyperspectral Image Analysis—Mapping of PFTs

The results of the mapping project indicate that the random forests supervised classification analysis was highly successful at distinguishing among cover types (Figure 10). The random forests model selection algorithm was set to separate cover into four classes based on the entire hyperspectral dataset of 360 spectral bands. The overall model had an out-of-bag error rate of 9.85% and a Cohen’s kappa coefficient of 0.869. Forest and shrub classes were classified with the greatest accuracy, and graminoid fen was the least accurate (Table 7). We found several spectral bands to be consistently important in making these classifications, as estimated by both the mean decrease in accuracy as well as the mean decrease in node impurity due to variable permutation (Table 8). Nearly all these bands fell within the 1800–2500 nm range representing the short-wave infrared region of the electromagnetic spectrum.

4. Discussion

4.1. Hyperspectral Remote Sensing of Peatland Response to Climate Change

By leveraging two large ecosystem manipulation experiments we have demonstrated that hyperspectral remote sensing can capture the response of peatland functional composition and species diversity to hydrologic change and increasing temperatures. Prior research has already demonstrated that hydrology, temperature, and CO2 affect peatland ecosystem structure and function [22,23,26,83,84,85]. Remote sensing studies in peatlands have also demonstrated the efficacy of hyperspectral methods in mapping peatland flora [17,86,87]. However, we have not encountered any study that has used hyperspectral remote sensing to track peatland ecological response to the diverse suite of global change drivers explored here. Our field results from the APEX site indicate that long-term hydrologic change in rich fen peatlands leads to shifts in species diversity and ecosystem productivity. The dominant source of variation within the principal component analysis of the APEX spectral data appeared largely driven by the greater diversity and overall cover associated with the control and raised plots, and a higher percent cover of litter associated with the lowered plots, which describes a gradient of ecosystem GPP [47]. At the SPRUCE experiment, we found that species composition varied among treatments. The results of our vegetation cover sampling indicated an increasing trend in shrub cover, and an associated decreasing trend in forb cover with warming. These results align with previous research linking warmer air temperatures and increases in soil nitrogen availability with an increased shrub growth in far northern ecosystems [4,88,89]. In our analysis of the SPRUCE spectral data, shrub and forb cover both emerged as statistically significant vectors within the multi-variate analysis of reflectance. Our results support the feasibility of tracking changes in PFT cover, in particular of shrubs, using hyperspectral methods.
We did not track non-vascular moss cover at SPRUCE because they were ubiquitous across the site. However, prior research indicates that the moisture status of peat-forming mosses can strongly influence reflectance spectra [85,90]. Further research at SPRUCE would characterize the role of moss moisture status on hyperspectral reflectance properties. However, focusing solely on vascular plant cover, our results highlight the potential to use hyperspectral data to identify and track a range of different climate change-induced impacts to northern boreal peatlands.

4.2. Remote Sensing of Boreal Peatland Species Diversity

In low-diversity ecosystems such as peatlands, the loss of species diversity may have cascading effects on ecosystem function and productivity [44]. The low functional redundancy of peatlands is such that the loss of a species may also represent the loss of an entire plant type, thereby negatively effecting ecological productivity [47]. We demonstrate that drought in peatlands negatively affects diversity, particularly of mosses and herbaceous species. Specifically, variation and heterogeneity in spectral reflectance increased with increasing Shannon diversity among plots. Increasing spectral heterogeneity has been previously used as a proxy for diversity in studies that have mapped peatland ecological diversity using remote sensing [91,92,93]. This type of study is based on the species–area relationship in which species richness increases with scale across a heterogeneous landscape [94]. These studies have extrapolated species richness by calculating indices such as the normalized difference vegetation index (NDVI), then calculating heterogeneity for pixel aggregates of varying spatial scales [53,55]. Where these approaches have modeled diversity as a function of area, our approach captures diversity directly by comparing the spectral heterogeneity among plots with varying levels of diversity [53]. Alternative approaches have attempted to map diversity through the identification of unique species [54,95,96,97]. However, this approach is challenging in non-forested systems in which individual organisms do not occupy an entire pixel. Our results indicate remote sensing may be used to estimate relative levels of diversity directly by comparing the spectral heterogeneity of measurements performed over a variety of different canopy types. We suggest that these results offer the ability to track diversity in systems with mixed or heterogeneous community assemblages that include trees, shrubs, graminiods, or other plant types [98].

4.3. Hyperspectral Characterization and Mapping of Plant Functional Types

Our results are in conversation with a growing body of research leveraging remotely sensed data to track vegetation change in peatland ecosystems [17,86,87]. The results of our analyses of both our plot-level and aerial datasets demonstrate the efficacy of using hyperspectral measurements to characterize the distribution of peatland PFTs. Through our close-range hyperspectral data collection we saw that the cover of PFTs was a driver of spectral variation. Despite relatively weak trends in our vegetation cover data at both sites, the reflectance spectra were highly responsive to subtle changes in the vegetation cover resulting from treatment. In particular, at SPRUCE, differences in the relative cover of shrubs versus forbs had a strong effect on spectral reflectance. This is likely due to the differences in stature, leaf structure, and foliar chemistry between high-growing woody shrubs compared with ground-layer herbaceous PFTs such as forbs [36,99]. We also predict that collecting data in the early fall improved our ability to detect variation among treatments because carotenoids and anthocyanins present in leaf tissues around senescence may have increased optical diversity, allowing for greater distinction among PFTs. Future research could examine the role of seasonality in improving the efficacy of hyperspectral remote sensing in making species-level distinctions.
Our plot-level findings prompted us to apply hyperspectral remote sensing to map peatland plant functional types over a larger spatial scale. The cover types represented within the aerial hyperspectral dataset included the PFTs that were present in our field measurements. We found the land-cover classification analysis was successful at distinguishing among different PFTs. This was particularly true for forest and shrub cover, which had the lowest error rate in the classification analysis. The percent cover of shrubs was also one of the strongest sources of variation in the field-based measurements at SPRUCE. We predict that the forested and shrub-covered regions of our study area were the most successfully mapped because their stature and leaf structure are clearly distinguishable from other PFTs, whereas graminiod fen and tussock grasses were less optically distinct from each other.
Our results indicate that reflectance in the short-wave infrared regions of the electromagnetic spectrum are particularly important in mapping peatland PFTs. Variable importance indicators such as model accuracy indicated that several SWIR bands between 1800 and 2500 nm were instrumental in creating the final classification analysis. These results support prior research in peatlands that has shown SWIR reflectance to be particularly sensitive to changes in peatland plant cover and moisture status [85]. Given the importance of changing shrub cover within the greater arctic and boreal landscape, these results support the application of hyperspectral measurements to track changes in the distribution of shrubs [4].
In arctic and boreal ecosystems where the environment is undergoing rapid change, the methods described here could prove invaluable to developing assessments of shifts in species composition over time. Due to the large spatial extent and low population density of the far north, examining changes to the landscape across broad spatial scales and over time must largely be achieved using remote sensing techniques. The use of hyperspectral imagery is integral to this approach because of the ability to select spectral bands that best capture variation across community types [90]. The high spatial resolution of aerial over satellite data is also key to this approach because it allows for fine-scale distinctions among cover classes. There are multiple research programs that collect hyperspectral data and make it publicly available to scientists. For example, NASA’s Airborne Visible/Infrared Imaging Spectrometer collected aerial hyperspectral imagery over the APEX field site in 2017 and 2018 as part of the Arctic and Boreal Vulnerability Experiment (ABoVE) [100]. With the launch of NASA’s Hyperspectral InfraRed Imager (HySPIRI) planned for 2021, global satellite hyperspectral imagery will become available to the remote sensing research community within the next several years [101]. The National Ecological Observatory Network also recently launched its hyperspectral Aerial Observation Platform with the mission of bridging spatial scales from organisms to landscapes [102]. A 2017 survey published by the National Academies of Science, Engineering and Medicine regarding the future of earth observation placed a high priority on hyperspectral data acquisition for the purposes of tracking diversity and changes in PFT cover (National Academies Press 2017). Our results establish statistical relationships between remotely sensed hyperspectral measurements, species diversity, and plant functional cover that could be used in future studies to identify and track community- and ecosystem-level changes in arctic and boreal regions.

5. Conclusions

The results from our study demonstrate that warming and hydrologic change lead to detectable changes in species composition and ecosystem productivity in boreal peatlands. We have further demonstrated the utility and efficacy of applying remotely sensed data toward characterizing and tracking changes in plant functional cover over time. Remote sensing is emerging as an effective tool to support ecological research, given its growing ability to characterize key elements of ecosystem structure and function such as leaf area, foliar chemistry, species diversity, and, more recently, the cover of different species and PFTs. We have shown that remote sensing can be used to detect the effects of a range of global change drivers on the productivity and functional composition of two types of peatland ecosystems representative of extensive areas of the northern boreal environment. We have also shown the feasibility of using hyperspectral remote sensing to map the distribution of PFTs over large spatial extents. In light of rapid global change, this approach could provide valuable insight into the changing ecology of the far north.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-4292/11/14/1685/s1, Figure S1: Mean treatment reflectance spectra for the APEX site. Five plots are averaged in the control, and four are averaged in the lowered and raised plots. Each plot-level average is comprised of three scans taken in succession, Figure S2: Mean treatment reflectance spectra for the SPRUCE site. Each spectra represents six plots across two chambers. Each plot-level average is comprised of three scans taken in succession, Figure S3: Distribution of data included in the linear modeling of the relationship between hyperspectral and plot-level characteristics. We also performed Anderson-Darling tests to support our use of these data in our models. We were able to reject the null hypothesis that the data were uniformly distributed in all cases, Figure S4: Spectral response curves for the cover types represented in the aerial hyperspectral dataset to demonstrate differentiation among cover types.

Author Contributions

Conceptualization, M.Y.M., M.J.F., E.S.K., M.R.T., T.A.D. and R.A.M.; Data curation, M.Y.M., J.A. and J.D.E.; Formal analysis, M.Y.M. and J.R.R.; Funding acquisition, M.J.F., E.S.K., R.K., M.R.T., T.A.D., B.P. and R.A.M.; Investigation, M.Y.M., M.J.F., E.S.K., J.A., J.D.E. and R.A.M.; Methodology, M.Y.M, J.R.R and R.A.M.; Project administration, M.Y.M; Resources, M.J.F, R.K., M.R.T., T.A.D., B.P. and R.A.M.; Supervision, M.J.F., E.S.K. and R.A.M.; Visualization, M.Y.M.; Writing – original draft, M.Y.M.; Writing – review & editing, M.J.F., J.R.R., E.S.K., R.K., M.R.T., T.A.D., J.A., J.D.E., B.P. and R.A.M.

Funding

This research was supported by the National Science Foundation grant #DEB LTREB 1354370 to (to E.S.K. and M.R.T.) and NASA grant #NNX14AF96G to (M.J.F. and E.S.K.), as well as cooperative agreements between the USDA Forest Service Northern Research Station and the University of Minnesota and Michigan Tech University. R.A.M. was supported by the Minnesota Agricultural Experiment Station project MIN-42-060. APEX is supported by funding for the Bonanza Creek LTER (funded jointly by National Science Foundation grant DEB-0423442 and USDA Forest Service, Pacific Northwest Research grant PNW01-JV11261952-231). SPRUCE is sponsored by the US Department of Energy Office of Science, by the Office of Biological and Environmental Research, and by the Oak Ridge National Laboratory. Generous funding for M.Y.M. was made possible by the University of Minnesota Department of Forest Resources. T.A.D. acknowledges the Engineer Research and Development Center Army Direct (6.1) Research Program and the Strategic Environmental Research and Development Program (projects RC-2110 and 18-1170).

Acknowledgments

We acknowledge the research groups at both the APEX and SPRUCE sites, as well as the efforts of numerous contributors including Danielle Rupp, Lucas Albano, Natalie Zwanenburg, Jamie Hollingsworth, Lynette Potvin, Joseph A. Kendrick, Samantha Radermacher, Anna Schweiger, Brian Benscoter, Doug Kastendick, Deacon Kyllander, and William Robert Nettles. We would also like to thank the contributions of four anonymous commentators for improving the quality of this manuscript during the review processes. The authors claim no conflicts of interest.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tape, K.; Sturm, M.; Racine, C. The evidence for shrub expansion in Northern Alaska and the Pan-Arctic. Glob. Chang. Biol. 2006, 12, 686–702. [Google Scholar] [CrossRef]
  2. Myers-Smith, I.H.; Harden, J.W.; Wilmking, M.; Fuller, C.C.; McGuire, A.D.; Chapin, F.S.I. Wetland Succession in a Permafrost Collapse: Interations between Fire and Thermokarst. Biogeosciences 2008, 5, 1273–1286. [Google Scholar] [CrossRef]
  3. Myers-Smith, I.H.; Elmendorf, S.C.; Beck, P.S.A.; Wilmking, M.; Hallinger, M.; Blok, D.; Tape, K.D.; Rayback, S.A.; Macias-Fauria, M.; Forbes, B.C.; et al. Climate sensitivity of shrub growth across the tundra biome. Nat. Clim. Chang. 2015, 5, 887–891. [Google Scholar] [CrossRef]
  4. Myers-Smith, I.H.; Forbes, B.C.; Wilmking, M.; Hallinger, M.; Lantz, T.; Blok, D.; Tape, K.D.; Macias-Fauria, M.; Sass-Klaassen, U.; Lévesque, E.; et al. Shrub Expansion in Tundra Ecosystems: Dynamics, Impacts and Research Priorities. Environ. Res. Lett. 2011, 6, 045509. [Google Scholar] [CrossRef]
  5. Johnstone, J.F.; Hollingsworth, T.N.; Mack, M.C.; Turetsky, M.; Chapin, F.S.; Romanovsky, V. Fire, climate change, and forest resilience in interior Alaska This article is one of a selection of papers from The Dynamics of Change in Alaska’s Boreal Forests: Resilience and Vulnerability in Response to Climate Warming. Can. J. For. Res. 2010, 40, 1302–1312. [Google Scholar] [CrossRef]
  6. Weintraub, M.N.; Schimel, J.P. Nitrogen Cycling and the Spread of Shrubs Control Changes in the Carbon Balance of Arctic Tundra Ecosystems. BioScience 2005, 55, 408–415. [Google Scholar] [CrossRef]
  7. Mack, M.C.; Schuur, E.A.G.; Bret-Harte, M.S.; Shaver, G.R.; Chapin, F.S. Ecosystem carbon storage in arctic tundra reduced by long-term nutrient fertilization. Nature 2004, 431, 440–443. [Google Scholar] [CrossRef] [PubMed]
  8. Douglas, T.A.; Jones, M.C.; Hiemstra, C.A.; Arnold, J.R. Sources and sinks of carbon in boreal ecosystems of interior Alaska: A review. Elem. Sci. Anth. 2014, 2. [Google Scholar] [CrossRef]
  9. Van Hemert, C.; Flint, P.L.; Udevitz, M.S.; Koch, J.C.; Atwood, T.C.; Oakley, K.L.; Pearce, J.M. Forecasting Wildlife Response to Rapid Warming in the Alaskan Arctic. BioScience 2015, 65, 718–728. [Google Scholar] [CrossRef] [Green Version]
  10. Berteaux, D.; Gauthier, G.; Domine, F.; Ims, R.A.; Lamoureux, S.F.; Lévesque, E.; Yoccoz, N. Effects of Changing Permafrost and Snow Conditions on Tundra Wildlife: Critical Places and Times. Arct. Sci. 2016, 3, 65–90. [Google Scholar] [CrossRef]
  11. Chapin, F.S., III; Sommerkorn, M.; Robards, M.D.; Hillmer-Pegram, K. Ecosystem stewardship: A resilience framework for arctic conservation. Glob. Environ. Chang. 2015, 34, 207–217. [Google Scholar] [CrossRef]
  12. Horstkotte, T.; Utsi, T.A.; Larsson-Blind, Å.; Burgess, P.; Johansen, B.; Käyhkö, J.; Oksanen, L.; Forbes, B.C. Human–Animal Agency in Reindeer Management: Sámi Herders’ Perspectives on Vegetation Dynamics under Climate Change. Ecosphere 2017, 8, e01931. [Google Scholar] [CrossRef]
  13. Gorham, E. Northern Peatlands: Role in the Carbon Cycle and Probable Responses to Climatic Warming. Ecol. Appl. 1991, 1, 182–195. [Google Scholar] [CrossRef] [PubMed]
  14. Yu, Z.; Loisel, J.; Brosseau, D.P.; Beilman, D.W.; Hunt, S.J. Global peatland dynamics since the Last Glacial Maximum. Geophys. Res. Lett. 2010, 37. [Google Scholar] [CrossRef]
  15. Treat, C.C.; Kleinen, T.; Broothaerts, N.; Dalton, A.S.; Dommain, R.; Douglas, T.A.; Drexler, J.Z.; Finkelstein, S.A.; Grosse, G.; Hope, G.; et al. Widespread global peatland establishment and persistence over the last 130,000 y. Proc. Natl. Acad. Sci. USA 2019, 116, 4822–4827. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Harris, A.; Gamon, J.A.; Pastorello, G.; Wong, C. Retrieval of the photochemical reflectance index for assessing xanthophyll cycle activity: A comparison of near-surface optical sensors. Biogeosciences 2014, 11, 6277–6292. [Google Scholar] [CrossRef]
  17. Erudel, T.; Fabre, S.; Houet, T.; Mazier, F.; Briottet, X. Criteria Comparison for Classifying Peatland Vegetation Types Using In Situ Hyperspectral Measurements. Remote Sens. 2017, 9, 748. [Google Scholar] [CrossRef]
  18. Kattge, J.; Díaz, S.; Lavorel, S.; Prentice, I.C.; Leadley, P.; Bonisch, G.; Garnier, E.; Westoby, M.; Reich, P.B.; Wright, I.J.; et al. TRY—A global database of plant traits. Glob. Chang. Biol. 2011, 17, 2905–2935. [Google Scholar] [CrossRef]
  19. Jetz, W.; Cavender-Bares, J.; Pavlick, R.; Schimel, D.; Davis, F.W.; Asner, G.P.; Guralnick, R.; Kattge, J.; Latimer, A.M.; Moorcroft, P.; et al. Monitoring plant functional diversity from space. Nat. Plants 2016, 2, 16024. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Weltzin, J.F.; Bridgham, S.D.; Pastor, J.; Chen, J.; Harth, C. Potential effects of warming and drying on peatland plant community composition. Glob. Chang. Biol. 2003, 9, 141–151. [Google Scholar] [CrossRef]
  21. Weltzin, J.F.; Pastor, J.; Harth, C.; Bridgham, S.D.; Updegraff, K.; Chapin, C.T. Response of bog and fen plant communities to warming and water-table manipulations. Ecology 2000, 81, 3464–3478. [Google Scholar] [CrossRef]
  22. Churchill, A.C.; Turetsky, M.R.; McGuire, A.D.; Hollingsworth, T.N. Response of Plant Community Structure and Primary Productivity to Experimental Drought and Flooding in an Alaskan Fen. Can. J. For. Res. 2014, 45, 185–193. [Google Scholar] [CrossRef]
  23. Dieleman, C.M.; Branfireun, B.A.; McLaughlin, J.W.; Lindo, Z. Climate Change Drives a Shift in Peatland Ecosystem Plant Community: Implications for Ecosystem Function and Stability. Glob. Chang. Biol. 2015, 21, 388–395. [Google Scholar] [CrossRef] [PubMed]
  24. Radu, D.D.; Duval, T.P. Precipitation frequency alters peatland ecosystem structure and CO2 exchange: Contrasting effects on moss, sedge, and shrub communities. Glob. Chang. Biol. 2018, 24, 2051–2065. [Google Scholar] [CrossRef] [PubMed]
  25. McPartland, M.Y.; Kane, E.S.; Falkowski, M.J.; Kolka, R.; Turetsky, M.R.; Palik, B.; Montgomery, R.A. The Response of Boreal Peatland Community Composition and NDVI to Hydrologic Change, Warming and Elevated Carbon Dioxide. Glob. Chang. Biol. 2019, 25, 93–107. [Google Scholar] [CrossRef]
  26. Olefeldt, D.; Euskirchen, E.S.; Harden, J.; Kane, E.; McGuire, A.D.; Waldrop, M.P.; Turetsky, M.R. A decade of boreal rich fen greenhouse gas fluxes in response to natural and experimental water table variability. Glob. Chang. Biol. 2017, 23, 2428–2440. [Google Scholar] [CrossRef] [PubMed]
  27. Christensen, T.R.; Johansson, T.; Åkerman, H.J.; Mastepanov, M.; Svensson, B.H.; Malmer, N.; Friborg, T.; Crill, P. Thawing sub-arctic permafrost: Effects on vegetation and methane emissions. Geophys. Res. Lett. 2004, 31. [Google Scholar] [CrossRef]
  28. Goud, E.M.; Moore, T.R.; Roulet, N.T. Predicting Peatland Carbon Fluxes from Non-Destructive Plant Traits. Funct. Ecol. 2018, 31, 1824–1833. [Google Scholar] [CrossRef]
  29. Chapin, F.S.; Hobbie, S.E.; Zhong, H.; Bret-Harte, M.S. Plant functional types as predictors of transient responses of arctic vegetation to global change. J. Veg. Sci. 1996, 7, 347–358. [Google Scholar] [CrossRef]
  30. Pielke, R.; Baldocchi, D.; Hobbie, S.E.; Roulet, N.; Eugster, W.; Kasischke, E.; Rastetter, E.; Zimov, S.A.; Running, S.W.; Chapin, F.S.; et al. Arctic and boreal ecosystems of western North America as components of the climate system. Glob. Chang. Biol. 2000, 6, 211–223. [Google Scholar] [CrossRef] [Green Version]
  31. Turner, W. Sensing Biodiversity. Science 2014, 346, 301–302. [Google Scholar] [CrossRef] [PubMed]
  32. Goetz, A.F.; Vane, G.; Solomon, J.E.; Rock, B.N. Imaging Spectrometry for Earth Remote Sensing. Science 1985, 228, 1147–1153. [Google Scholar] [CrossRef] [PubMed]
  33. Ustin, S.L.; Gamon, J.A. Remote sensing of plant functional types. New Phytol. 2010, 186, 795–816. [Google Scholar] [CrossRef] [PubMed]
  34. Ustin, S.L.; Roberts, D.A.; Gamon, J.A.; Asner, G.P.; Green, R.O. Using Imaging Spectroscopy to Study Ecosystem Processes and Properties. BioScience 2004, 54, 523–534. [Google Scholar] [CrossRef]
  35. Lichtenthaler, H.K.; Gitelson, A.; Lang, M. Non-Destructive Determination of Chlorophyll Content of Leaves of a Green and an Aurea Mutant of Tobacco by Reflectance Measurements. J. Plant Physiol. 1996, 148, 483–493. [Google Scholar] [CrossRef]
  36. Sims, D.A.; Gamon, J.A. Relationships between leaf pigment content and spectral reflectance across a wide range of species, leaf structures and developmental stages. Remote Sens. Environ. 2002, 81, 337–354. [Google Scholar] [CrossRef]
  37. Blackburn, G.A. Hyperspectral Remote Sensing of Plant Pigments. J. Exp. Bot. 2007, 58, 855–867. [Google Scholar] [CrossRef]
  38. Serbin, S.P.; Singh, A.; McNeil, B.E.; Kingdon, C.C.; Townsend, P.A. Spectroscopic determination of leaf morphological and biochemical traits for northern temperate and boreal tree species. Ecol. Appl. 2014, 24, 1651–1669. [Google Scholar] [CrossRef] [Green Version]
  39. Martin, M.E.; Aber, J.D. High spectral resolution remote sensing of forest canopy lignin, nitrogen, and ecosystem processes. Ecol. Appl. 1997, 7, 431–443. [Google Scholar] [CrossRef]
  40. Knyazikhin, Y.; Schull, M.A.; Stenberg, P.; Mõttus, M.; Rautiainen, M.; Yang, Y.; Marshak, A.; Carmona, P.L.; Kaufmann, R.K.; Lewis, P.; et al. Hyperspectral Remote Sensing of Foliar Nitrogen Content. Proc. Natl. Acad. Sci. USA 2013, 110, 811–812. [Google Scholar] [CrossRef]
  41. Peñuelas, J.; Pinol, J.; Ogaya, R.; Filella, I. Estimation of plant water concentration by the reflectance Water Index WI (R900/R970). Int. J. Remote Sens. 1997, 18, 2869–2875. [Google Scholar] [CrossRef]
  42. Schweiger, A.K.; Cavender-Bares, J.; Townsend, P.A.; Hobbie, S.E.; Madritch, M.D.; Wang, R.; Tilman, D.; Gamon, J.A. Plant spectral diversity integrates functional and phylogenetic components of biodiversity and predicts ecosystem function. Nat. Ecol. Evol. 2018, 2, 976–982. [Google Scholar] [CrossRef] [PubMed]
  43. Turner, W.; Spector, S.; Gardiner, N.; Fladeland, M.; Sterling, E.; Steininger, M. Remote sensing for biodiversity science and conservation. Trends Ecol. Evol. 2003, 18, 306–314. [Google Scholar] [CrossRef]
  44. Tilman, D. The Influence of Functional Diversity and Composition on Ecosystem Processes. Science 1997, 277, 1300–1302. [Google Scholar] [CrossRef] [Green Version]
  45. Walker, M.D.; Wahren, C.H.; Hollister, R.D.; Henry, G.H.R.; Ahlquist, L.E.; Alatalo, J.M.; Bret-Harte, M.S.; Calef, M.P.; Callaghan, T.V.; Carroll, A.B.; et al. Plant community responses to experimental warming across the tundra biome. Proc. Natl. Acad. Sci. USA 2006, 103, 1342–1346. [Google Scholar] [CrossRef] [Green Version]
  46. Pajunen, A.; Oksanen, J.; Virtanen, R. Impact of shrub canopies on understorey vegetation in western Eurasian tundra. J. Veg. Sci. 2011, 22, 837–846. [Google Scholar] [CrossRef]
  47. Reich, P.B.; Tilman, D.; Isbell, F.; Mueller, K.; Hobbie, S.E.; Flynn, D.F.B.; Eisenhauer, N. Impacts of Biodiversity Loss Escalate Through Time as Redundancy Fades. Science 2012, 336, 589–592. [Google Scholar] [CrossRef] [Green Version]
  48. Ward, S.E.; Ostle, N.J.; Oakley, S.; Quirk, H.; Henrys, P.A.; Bardgett, R.D. Warming effects on greenhouse gas fluxes in peatlands are modulated by vegetation composition. Ecol. Lett. 2013, 16, 1285–1293. [Google Scholar] [CrossRef] [Green Version]
  49. Wang, R.; Gamon, J.A.; Emmerton, C.A.; Li, H.; Nestola, E.; Pastorello, G.Z.; Menzer, O. Integrated Analysis of Productivity and Biodiversity in a Southern Alberta Prairie. Remote Sens. 2016, 8, 214. [Google Scholar] [CrossRef]
  50. Wang, R.; Gamon, J.A.; Schweiger, A.K.; Cavender-Bares, J.; Townsend, P.A.; Zygielbaum, A.I.; Kothari, S. Influence of species richness, evenness, and composition on optical diversity: A simulation study. Remote Sens. Environ. 2018, 211, 218–228. [Google Scholar] [CrossRef]
  51. Palmer, M.W.; Earls, P.G.; Hoagland, B.W.; White, P.S.; Wohlgemuth, T. Quantitative tools for perfecting species lists. Environmetrics 2002, 13, 121–137. [Google Scholar] [CrossRef]
  52. Rocchini, D. Effects of spatial and spectral resolution in estimating ecosystem α-diversity by satellite imagery. Remote Sens. Environ. 2007, 111, 423–434. [Google Scholar] [CrossRef]
  53. Rocchini, D.; Balkenhol, N.; Carter, G.A.; Foody, G.M.; Gillespie, T.W.; He, K.S.; Kark, S.; Levin, N.; Lucas, K.; Luoto, M.; et al. Remotely sensed spectral heterogeneity as a proxy of species diversity: Recent advances and open challenges. Ecol. Inform. 2010, 5, 318–329. [Google Scholar] [CrossRef]
  54. Baldeck, C.A.; Asner, G.P. Estimating Vegetation Beta Diversity from Airborne Imaging Spectroscopy and Unsupervised Clustering. Remote Sens. 2013, 5, 2057–2071. [Google Scholar] [CrossRef] [Green Version]
  55. Gould, W. Remote sensing of vegetation, plant species richness, and regional biodiversity hotspots. Ecol. Appl. 2000, 10, 1861–1870. [Google Scholar] [CrossRef]
  56. Carter, G.A.; Knapp, A.K.; Anderson, J.E.; Hoch, G.A.; Smith, M.D. Indicators of plant species richness in AVIRIS spectra of a mesic grassland. Remote Sens. Environ. 2005, 98, 304–316. [Google Scholar] [CrossRef]
  57. Levin, N.; Shmida, A.; Levanoni, O.; Tamari, H.; Kark, S. Predicting mountain plant richness and rarity from space using satellite-derived vegetation indices. Divers. Distrib. 2007, 13, 692–703. [Google Scholar] [CrossRef]
  58. Lucas, K.; Carter, G. The use of hyperspectral remote sensing to assess vascular plant species richness on Horn Island, Mississippi. Remote Sens. Environ. 2008, 112, 3908–3915. [Google Scholar] [CrossRef]
  59. Oldeland, J.; Wesuls, D.; Rocchini, D.; Schmidt, M.; Jürgens, N. Does using species abundance data improve estimates of species diversity from remotely sensed spectral heterogeneity? Ecol. Indic. 2010, 10, 390–396. [Google Scholar] [CrossRef]
  60. Pereira, H.M.; Ferrier, S.; Walters, M.; Geller, G.N.; Jongman, R.H.G.; Scholes, R.J.; Bruford, M.W.; Brummitt, N.; Butchart, S.H.M.; Cardoso, A.C.; et al. Essential Biodiversity Variables. Science 2013, 339, 277–278. [Google Scholar] [CrossRef] [Green Version]
  61. Proença, V.; Martin, L.J.; Pereira, H.M.; Fernandez, M.; McRae, L.; Belnap, J.; Böhm, M.; Brummitt, N.; García-Moreno, J.; Gregory, R.D.; et al. Global biodiversity monitoring: From data sources to Essential Biodiversity Variables. Biol. Conserv. 2017, 213, 256–263. [Google Scholar] [CrossRef] [Green Version]
  62. Vihervaara, P.; Auvinen, A.-P.; Mononen, L.; Törmä, M.; Ahlroth, P.; Anttila, S.; Böttcher, K.; Forsius, M.; Heino, J.; Heliölä, J.; et al. How Essential Biodiversity Variables and Remote Sensing Can Help National Biodiversity Monitoring. Glob. Ecol. Conserv. 2017, 10, 43–59. [Google Scholar] [CrossRef]
  63. Hanson, P.J.; Riggs, J.S.; Nettles, W.R.; Phillips, J.R.; Krassovski, M.B.; Hook, L.A.; Gu, L.; Richardson, A.D.; Aubrecht, D.M.; Ricciuto, D.M.; et al. Attaining whole-ecosystem warming using air and deep-soil heating methods with an elevated CO2 atmosphere. Biogeosciences 2017, 14, 861–883. [Google Scholar] [CrossRef]
  64. Jonasson, S. Evaluation of the Point Intercept Method for the Estimation of Plant Biomass. Oikos 1988, 52, 101–106. [Google Scholar] [CrossRef]
  65. ASD. Fieldspec 3 User Manual; ASD Inc.: Boulder, CO, USA, 2007. [Google Scholar]
  66. Herweg, J.A.; Kerekes, J.P.; Weatherbee, O.; Messinger, D.; Van Aardt, J.; Ientilucci, E.; Ninkov, Z.; Faulring, J.; Raqueno, N.; Meola, J. SpecTIR hyperspectral airborne Rochester experiment data collection campaign. SPIE Def. Secur. Sens. 2012, 8390, 839028. [Google Scholar] [CrossRef]
  67. Kruse, F.; Lefkoff, A.; Boardman, J.; Heidebrecht, K.; Shapiro, A.; Barloon, P.; Goetz, A. The spectral image processing system (SIPS)—Interactive visualization and analysis of imaging spectrometer data. Remote Sens. Environ. 1993, 44, 145–163. [Google Scholar] [CrossRef]
  68. Anderson, J.E.; Douglas, T.A.; Barbato, R.A.; Saari, S.; Jones, R.M. Vegetation Mapping and Seasonal Thaw Estimates in Interior Alaska Permafrost. Remote Sens. Environ. In review.
  69. Nakagawa, S.; Schielzeth, H. A General and Simple Method for Obtaining R2 from Generalized Linear Mixed-Effects Models. Methods Ecol. Evol. 2013, 4, 133–142. [Google Scholar] [CrossRef]
  70. Gotelli, N.J.; Ellison, A.M. A Primer of Ecological Statistics, 2nd ed.; Sinauer Associates Is an Imprint of Oxford University Press: Sunderland, MA, USA, 2012; ISBN 978-1-60535-064-6. [Google Scholar]
  71. Oksanen, J.; Blanchet, F.G.; Friendly, M.; Kindt, R.; Legendre, P.; McGlinn, D.; Minchin, P.R.; O’Hara, R.B.; Simpson, G.L.; Solymos, P.; et al. vegan: Community Ecology Package. R package version 2.4-4. 2017. Available online: https://CRAN.R-project.org/package=vegan (accessed on 16 July 2019).
  72. Zuur, A.F.; Ieno, E.N.; Smith, G.M. (Eds.) Principal component analysis and redundancy analysis. In Analysing Ecological Data; Statistics for Biology and Health; Springer: New York, NY, USA, 2007; pp. 193–224. [Google Scholar]
  73. Legendre, P.; Legendre, L.F.J. Numerical Ecology; Elsevier: Oxford, UK, 2012; ISBN 978-0-444-53869-7. [Google Scholar]
  74. Harsanyi, J.C.; Chang, C.-I. Hyperspectral image classification and dimensionality reduction: An orthogonal subspace projection approach. IEEE Trans. Geosci. Remote. Sens. 1994, 32, 779–785. [Google Scholar] [CrossRef]
  75. Li, W.; Prasad, S.; Fowler, J.E.; Bruce, L.M. Locality-Preserving Dimensionality Reduction and Classification for Hyperspectral Image Analysis. IEEE Trans. Geosci. Remote Sens. 2012, 50, 1185–1198. [Google Scholar] [CrossRef]
  76. Li, H.; Reynolds, J.F. On Definition and Quantification of Heterogeneity. Oikos 1995, 73, 280–284. [Google Scholar] [CrossRef]
  77. Breiman, L. Random forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef]
  78. Liaw, A.; Wiener, M. Classification and Regression by RandomForest. R News 2003, 2, 18–22. [Google Scholar]
  79. Bivand, R.; Keitt, T.; Rowlingson, B. rgdal: Bindings for the ‘Geospatial’ Data Abstraction Library. R package version 1.2-16. 2017. Available online: https://CRAN.R-project.org/package=rgdal (accessed on 16 July 2019).
  80. Evans, J.S.; Murphy, M.A.; Holden, Z.A.; Cushman, S.A. Modeling Species Distribution and Change Using Random Forest. In Predictive Species and Habitat Modeling in Landscape Ecology: Concepts and Applications; Springer: Drew, CA, USA, 2011; ISBN 978-1-4419-7390-0. [Google Scholar]
  81. Mitchell, M.W. Bias of the Random Forest Out-of-Bag (OOB) Error for Certain Input Parameters. Open J. Stat. 2011, 1, 205–211. [Google Scholar] [CrossRef] [Green Version]
  82. QGIS Development Team. QGIS Geographic Information System. Open Source Geospatial Foundation Project. 2018. Available online: http://qgis.osgeo.org (accessed on 16 July 2019).
  83. Boelman, N.T.; Stieglitz, M.; Rueth, H.M.; Sommerkorn, M.; Griffin, K.L.; Shaver, G.R.; Gamon, J.A. Response of NDVI, biomass, and ecosystem gas exchange to long-term warming and fertilization in wet sedge tundra. Oecologia 2003, 135, 414–421. [Google Scholar] [CrossRef] [PubMed]
  84. Chivers, M.R.; Turetsky, M.R.; Waddington, J.M.; Harden, J.W.; McGuire, A.D. Effects of Experimental Water Table and Temperature Manipulations on Ecosystem CO2 Fluxes in an Alaskan Rich Fen. Ecosystems 2009, 12, 1329–1342. [Google Scholar] [CrossRef]
  85. Meingast, K.M.; Falkowski, M.J.; Kane, E.S.; Potvin, L.R.; Benscoter, B.W.; Smith, A.M.; Bourgeau-Chavez, L.L.; Miller, M.E. Spectral detection of near-surface moisture content and water-table position in northern peatland ecosystems. Remote Sens. Environ. 2014, 152, 536–546. [Google Scholar] [CrossRef]
  86. Schmidtlein, S.; Feilhauer, H.; Bruelheide, H. Mapping Plant Strategy Types Using Remote Sensing. J. Veg. Sci. 2012, 23, 395–405. [Google Scholar] [CrossRef]
  87. Harris, A.; Charnock, R.; Lucas, R. Hyperspectral remote sensing of peatland floristic gradients. Remote Sens. Environ. 2015, 162, 99–111. [Google Scholar] [CrossRef] [Green Version]
  88. Chapin, F.S., III; Shaver, G.R.; Giblin, A.E.; Nadelhoffer, K.J.; Laundre, J.A. Responses of Arctic Tundra to Experimental and Observed Changes in Climate. Ecology 1995, 76, 694–711. [Google Scholar] [CrossRef]
  89. Shaver, G.R.; Bret-Harte, M.S.; Chapin, F.S. Primary and secondary stem growth in arctic shrubs: Implications for community response to environmental change. J. Ecol. 2002, 90, 251–267. [Google Scholar] [CrossRef]
  90. Bratsch, S.N.; Epstein, H.E.; Buchhorn, M.; Walker, D.A. Differentiating among Four Arctic Tundra Plant Communities at Ivotuk, Alaska Using Field Spectroscopy. Remote Sens. 2016, 8, 51. [Google Scholar] [CrossRef]
  91. Poulin, M.F.; Careau, D.; Rochefort, L.; DesRochers, A. From Satellite Imagery to Peatland Vegetation Diversity: How Reliable Are Habitat Maps? Conserv. Ecol. 2002, 6. [Google Scholar] [CrossRef] [Green Version]
  92. Cabezas, J.; Galleguillos, M.; Perez-Quezada, J.F. Predicting Vascular Plant Richness in a Heterogeneous Wetland Using Spectral and Textural Features and a Random Forest Algorithm. IEEE Geosci. Remote Sens. Lett. 2016, 13, 1–5. [Google Scholar] [CrossRef]
  93. Lopatin, J.; Kattenborn, T.; Galleguillos, M.; Perez-Quezada, J.F.; Schmidtlein, S. Using aboveground vegetation attributes as proxies for mapping peatland belowground carbon stocks. Remote Sens. Environ. 2019, 231, 111217. [Google Scholar] [CrossRef]
  94. Whittaker, R.H. Evolution and Measurement of Species Diversity. Taxon 1972, 21, 213–251. [Google Scholar] [CrossRef]
  95. Baldeck, C.A.; Colgan, M.S.; Féret, J.-B.; Levick, S.R.; Martin, R.E.; Asner, G.P.; Asner, G. Landscape-scale variation in plant community composition of an African savanna from airborne species mapping. Ecol. Appl. 2014, 24, 84–93. [Google Scholar] [CrossRef]
  96. Carlson, K.M.; Asner, G.P.; Hughes, R.F.; Ostertag, R.; Martin, R.E. Hyperspectral Remote Sensing of Canopy Biodiversity in Hawaiian Lowland Rainforests. Ecosystems 2007, 10, 536–549. [Google Scholar] [CrossRef]
  97. Féret, J.-B.; Asner, G.P. Mapping Tropical Forest Canopy Diversity Using High-Fidelity Imaging Spectroscopy. Ecological Applications 2014, 24(6), 1289–1296. [Google Scholar] [CrossRef]
  98. Feilhauer, H.; Faude, U.; Schmidtlein, S. Combining Isomap ordination and imaging spectroscopy to map continuous floristic gradients in a heterogeneous landscape. Remote Sens. Environ. 2011, 115, 2513–2524. [Google Scholar] [CrossRef]
  99. Kattenborn, T.; Fassnacht, F.E.; Schmidtlein, S. Differentiating Plant Functional Types Using Reflectance: Which Traits Make the Difference? Remote Sens. Ecol. Conserv. 2019, 5, 5–19. [Google Scholar] [CrossRef]
  100. Miller, C.E.; Green, R.O.; Thompson, D.R.; Thorpe, A.K.; Eastwood, M.; Mccubbin, I.B.; Olson-Duvall, W.; Bernas, M.; Sarture, C.M.; Nolte, S.; et al. ABoVE: Hyperspectral Imagery from AVIRIS-NG, Alaskan and Canadian Arctic, 2017–2018. ORNL DAAC 2018. [Google Scholar] [CrossRef]
  101. Lee, C.M.; Cable, M.L.; Hook, S.J.; Green, R.O.; Ustin, S.L.; Mandl, D.J.; Middleton, E.M. An Introduction to the NASA Hyperspectral InfraRed Imager (HyspIRI) Mission and Preparatory Activities. Remote Sens. Environ. 2015, 167, 6–19. [Google Scholar] [CrossRef]
  102. Kampe, T.U.; Johnson, B.R.; Kuester, M.A.; Keller, M. NEON: The first continental-scale ecological observatory with airborne remote sensing of vegetation canopy biochemistry and structure. J. Appl. Remote Sens. 2010, 4, 043510. [Google Scholar] [CrossRef]
Figure 1. Location of the Alaska Peatland Experiment (APEX) and Spruce and Peatland Responses Under Changing Environments (SPRUCE) peatland field sites. APEX is located in central Alaska. SPRUCE is located in northern Minnesota.
Figure 1. Location of the Alaska Peatland Experiment (APEX) and Spruce and Peatland Responses Under Changing Environments (SPRUCE) peatland field sites. APEX is located in central Alaska. SPRUCE is located in northern Minnesota.
Remotesensing 11 01685 g001
Figure 2. (Left) Aerial image of APEX water table manipulation treatments and the surrounding area (photo credit: Evan Schjins). (Right) Aerial image of SPRUCE experimental chambers (photo credit: US Department of Energy, Oak Ridge National Laboratory).
Figure 2. (Left) Aerial image of APEX water table manipulation treatments and the surrounding area (photo credit: Evan Schjins). (Right) Aerial image of SPRUCE experimental chambers (photo credit: US Department of Energy, Oak Ridge National Laboratory).
Remotesensing 11 01685 g002
Figure 3. Bands included in SPRUCE and APEX spectral sampling. Dashed line indicates SPRUCE data and dotted lines indicate APEX data. Data represent sample spectral scans collected at each site.
Figure 3. Bands included in SPRUCE and APEX spectral sampling. Dashed line indicates SPRUCE data and dotted lines indicate APEX data. Data represent sample spectral scans collected at each site.
Remotesensing 11 01685 g003
Figure 4. (Left) Unmanned aerial image rendered in three visible bands, including a near-infrared band. Image has been clipped to the extent of the hyperspectral aerial image. The UAS image was used to select supervised classification training points. (Right) Extent of hyperspectral aerial dataset, rendered in true-color RGB bands at 457, 570, and 650 nm. The APEX water table manipulation areas can be seen as the white squares to the north of the images.
Figure 4. (Left) Unmanned aerial image rendered in three visible bands, including a near-infrared band. Image has been clipped to the extent of the hyperspectral aerial image. The UAS image was used to select supervised classification training points. (Right) Extent of hyperspectral aerial dataset, rendered in true-color RGB bands at 457, 570, and 650 nm. The APEX water table manipulation areas can be seen as the white squares to the north of the images.
Remotesensing 11 01685 g004
Figure 5. Results of APEX vegetation cover sampling. Bars represent the treatment mean. Error bars represent standard deviation. Percent cover was calculated for each PFT as a fraction of the total number of observations. Data were analyzed for statistical significance using a one-way ANOVA with a post-hoc pairwise means comparison for significant difference among treatments. Bars with the same symbol are significantly different from each other (ANOVA, p < 0.05).
Figure 5. Results of APEX vegetation cover sampling. Bars represent the treatment mean. Error bars represent standard deviation. Percent cover was calculated for each PFT as a fraction of the total number of observations. Data were analyzed for statistical significance using a one-way ANOVA with a post-hoc pairwise means comparison for significant difference among treatments. Bars with the same symbol are significantly different from each other (ANOVA, p < 0.05).
Remotesensing 11 01685 g005
Figure 6. Percent PFT cover in response to warming at the SPRUCE project. Linear mixed effects models were fit to the PFT data with temperature as a continuous variable, CO2 level as a factor, and chamber included as random effect.
Figure 6. Percent PFT cover in response to warming at the SPRUCE project. Linear mixed effects models were fit to the PFT data with temperature as a continuous variable, CO2 level as a factor, and chamber included as random effect.
Remotesensing 11 01685 g006
Figure 7. Shannon alpha diversity in response to treatment at APEX. Boxes represent the interquartile range, the central black bar indicates the treatment mean, and error bars represent standard deviation. Number of plots included per treatment: control = 5, lowered = 4, raised = 4.
Figure 7. Shannon alpha diversity in response to treatment at APEX. Boxes represent the interquartile range, the central black bar indicates the treatment mean, and error bars represent standard deviation. Number of plots included per treatment: control = 5, lowered = 4, raised = 4.
Remotesensing 11 01685 g007
Figure 8. (Left) Principal component analysis for APEX spectral data demonstrating distribution of treatment plots with community variables mapped. Arrow length represents the strength of the gradient. (Right) Principal component analysis of spectral data at SPRUCE, demonstrating the distribution of temperature treatments (shown in °C) with community variables mapped on. The length of the arrow represents the strength of the gradient.
Figure 8. (Left) Principal component analysis for APEX spectral data demonstrating distribution of treatment plots with community variables mapped. Arrow length represents the strength of the gradient. (Right) Principal component analysis of spectral data at SPRUCE, demonstrating the distribution of temperature treatments (shown in °C) with community variables mapped on. The length of the arrow represents the strength of the gradient.
Remotesensing 11 01685 g008
Figure 9. Linear relationship between spectral and community variation. (Left) Shannon diversity and spectral coefficient of variation. (Right) Spectral heterogeneity and community heterogeneity measured as Euclidian distance for a global centroid calculated using a principal component analysis.
Figure 9. Linear relationship between spectral and community variation. (Left) Shannon diversity and spectral coefficient of variation. (Right) Spectral heterogeneity and community heterogeneity measured as Euclidian distance for a global centroid calculated using a principal component analysis.
Remotesensing 11 01685 g009
Figure 10. Results of supervised classification analysis of PFT distribution at APEX and the surrounding area made using a random forests modeling approach.
Figure 10. Results of supervised classification analysis of PFT distribution at APEX and the surrounding area made using a random forests modeling approach.
Remotesensing 11 01685 g010
Table 1. Description of two study sites, data types, and differing data collection methods between sites.
Table 1. Description of two study sites, data types, and differing data collection methods between sites.
APEXSPRUCE
Peatland typeRich fenOmbrotrophic bog
LocationAlaska, USAMinnesota, USA
Experimental designWater table manipulation with 120 m2 control, lowered and raised treatmentsRegression-based factorial between increasing temperature and CO2 level
Vegetation samplingPoint bar laser survey method2 m2 sampling frame method
Spectral reflectanceASD Fieldspec ProUnispec DC
Table 2. List of vascular plant species present at both study sites and plant functional type groupings used in analysis.
Table 2. List of vascular plant species present at both study sites and plant functional type groupings used in analysis.
SPRUCEAPEX
PFTSpeciesPFTSpecies
ForbDrosera rotundifoliaEquisetumEquisetum fluviatile
ForbMaianthemum trifoliumForbGalium trifidum
GraminoidCarex magellanicaForbPotamogeton gramineus
GraminoidCarex oligospermaGraminoidCalamagrostis canadensis
GraminoidCarex trispermaGraminoidCarex loliacea
GraminoidEriophorum vaginatumGraminoidCarex utriculata
GraminoidEriophorum virginicumMossSphagnum spp.
ShrubAndromeda polifoliaShrubPotentilla palustris
ShrubChamaedaphne calyculata
ShrubKalmia polifolia
ShrubRhododendron groenlandicum
ShrubVaccinium angustifolium
ShrubVaccinium oxycoccos
Table 3. APEX ANOVA results. Italicized values indicate statistical significance at an alpha level of 0.05.
Table 3. APEX ANOVA results. Italicized values indicate statistical significance at an alpha level of 0.05.
PFTSum SqMean SqF2, 10P
Forb19.249.621.4780.274
Sedge182.991.462.0940.174
Shrub13.366.732.7240.114
Equisetum68.7734.392.030.182
Grass402.7201.41.9970.186
Moss182.791.378.7790.006
Litter980.9490.54.1350.049
Diversity0.6500.333.5160.07
Table 4. SPRUCE linear mixed effects model results for the response of cover to warming and CO2 level, with chamber added as an error term.
Table 4. SPRUCE linear mixed effects model results for the response of cover to warming and CO2 level, with chamber added as an error term.
PFTMarginal R2Conditional R2AICβTEMPβCO2βTEMP:CO2
Forb0.2780.605342−4.444−6.1431.170
Sedge0.0170.1534120.75815.178−3.484
Shrub0.1550.251407−9.932−28.63312.935
Diversity0.0230.285−37−0.0030.024−0.005
Table 5. Results from the APEX permutation test of the relationship between variance within spectral principal component analysis and the percent cover of different plant functional types. Asterisk indicates statistical significance at a 95% confidence level.
Table 5. Results from the APEX permutation test of the relationship between variance within spectral principal component analysis and the percent cover of different plant functional types. Asterisk indicates statistical significance at a 95% confidence level.
VariablePC1PC2r2p
Diversity–0.997–0.0740.4950.038 *
Equisetum−0.9900.1420.1850.370
Forbs−0.839−0.5450.1510.458
Sedges−0.405−0.9140.1380.491
Grasses0.2240.9750.0160.910
Litter0.8980.4400.3050.169
Moss−0.837−0.5470.3640.108
Shrubs−0.922−0.3860.1630.403
LAI−0.313−0.9500.1260.486
Table 6. Results from the SPRUCE permutation test of the relationships between the spectral principal component analysis, percent cover of different plant functional types, and plot-level species diversity. One asterisk indicates statistical significance at a 95% confidence level, and two asterisks indicate significance at a 99% confidence interval. LAI refers to leaf area index.
Table 6. Results from the SPRUCE permutation test of the relationships between the spectral principal component analysis, percent cover of different plant functional types, and plot-level species diversity. One asterisk indicates statistical significance at a 95% confidence level, and two asterisks indicate significance at a 99% confidence interval. LAI refers to leaf area index.
VariablePC1PC2r2p
Forbs0.9260.3770.2440.011 *
Sedges0.2710.9630.0390.520
Shrubs0.347–0.9380.2620.007 **
Trees0.906−0.4230.0330.581
Diversity0.7710.6370.0950.204
Table 7. Results of supervised random forests classification analysis.
Table 7. Results of supervised random forests classification analysis.
ForestGraminoidShrubTussockClass Error
Forest945000.051
Graminoid fen7802100.192
Shrub139500.040
Tussock Grass191880.111
Table 8. Random forests model output of variable importance for spectral bands used to perform the land cover classification analysis.
Table 8. Random forests model output of variable importance for spectral bands used to perform the land cover classification analysis.
Model AccuracyGini/Impurity
Band numberMean DecreaseBand numberMean decrease
12107.5312614.68
12336.9912104.32
12616.7715463.78
11936.7612553.71
24416.7112333.60
12046.6615523.49
24476.6611873.47
24306.6112443.40
24126.5612673.35
4746.5212153.29

Share and Cite

MDPI and ACS Style

McPartland, M.Y.; Falkowski, M.J.; Reinhardt, J.R.; Kane, E.S.; Kolka, R.; Turetsky, M.R.; Douglas, T.A.; Anderson, J.; Edwards, J.D.; Palik, B.; et al. Characterizing Boreal Peatland Plant Composition and Species Diversity with Hyperspectral Remote Sensing. Remote Sens. 2019, 11, 1685. https://doi.org/10.3390/rs11141685

AMA Style

McPartland MY, Falkowski MJ, Reinhardt JR, Kane ES, Kolka R, Turetsky MR, Douglas TA, Anderson J, Edwards JD, Palik B, et al. Characterizing Boreal Peatland Plant Composition and Species Diversity with Hyperspectral Remote Sensing. Remote Sensing. 2019; 11(14):1685. https://doi.org/10.3390/rs11141685

Chicago/Turabian Style

McPartland, Mara Y., Michael J. Falkowski, Jason R. Reinhardt, Evan S. Kane, Randy Kolka, Merritt R. Turetsky, Thomas A. Douglas, John Anderson, Jarrod D. Edwards, Brian Palik, and et al. 2019. "Characterizing Boreal Peatland Plant Composition and Species Diversity with Hyperspectral Remote Sensing" Remote Sensing 11, no. 14: 1685. https://doi.org/10.3390/rs11141685

APA Style

McPartland, M. Y., Falkowski, M. J., Reinhardt, J. R., Kane, E. S., Kolka, R., Turetsky, M. R., Douglas, T. A., Anderson, J., Edwards, J. D., Palik, B., & Montgomery, R. A. (2019). Characterizing Boreal Peatland Plant Composition and Species Diversity with Hyperspectral Remote Sensing. Remote Sensing, 11(14), 1685. https://doi.org/10.3390/rs11141685

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop