Next Article in Journal
Investigations of the Chemical Distribution in Sorbitol and Citric Acid (SorCA) Treated Wood—Development of a Quality Control Method on the Basis of Electromagnetic Radiation
Previous Article in Journal
Impacts of Climate Change on Blue Carbon Stocks and Fluxes in Mangrove Forests
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Managing Moose from Home: Determining Landscape Carrying Capacity for Alces alces Using Remote Sensing

by
David W. Kramer
1,*,
Thomas J. Prebyl
2,
Nathan P. Nibbelink
2,
Karl V. Miller
2,
Alejandro A. Royo
3 and
Jacqueline L. Frair
4
1
New York Department of Environmental Conservation, Division of Fish and Wildlife, 625 Broadway, Albany, NY 12233, USA
2
D. B. Warnell School of Forestry and Natural Resources, University of Georgia, Athens, GA 30602, USA
3
Forestry Sciences Laboratory, Northern Research Station, USDA Forest Service, Irvine, PA 16329, USA
4
College of Environmental Science and Forestry, State University of New York, 1 Forestry Drive, Syracuse, NY 13210, USA
*
Author to whom correspondence should be addressed.
Forests 2022, 13(2), 150; https://doi.org/10.3390/f13020150
Submission received: 30 November 2021 / Revised: 6 January 2022 / Accepted: 14 January 2022 / Published: 19 January 2022
(This article belongs to the Section Forest Inventory, Modeling and Remote Sensing)

Abstract

:
In temperate forests of the northeastern U.S., moose (Alces alces) populations are adapted for mixed-age heterogeneous landscapes that provide abundant herbaceous forage in warm months and coniferous forage during winter. Heterogeneity of forest stands is driven by management activities or natural disturbance, resulting in a multi-age forest at a landscape scale. Here, we present a method to estimate landscape carrying capacity of moose by combining remote sensing classification of forest cover class with literature or field-based estimates of class-specific forage abundance. We used Landsat imagery from 1991 to 2013 for the Allegheny National Forest and 2013–2018 for the Adirondack Park, and associated training polygons, to predict based on NDVI and SWI whether a forested landscape fit into one of three cover classes: mature forest, intermediate timber removal, or overstory timber removal. Our three-classes yielded a mean land cover prediction accuracy of 94.3% (Khat = 0.91) and 86.9% (Khat = 0.76) for ANFR and AP, respectively. In the AP, we applied previously calculated summer crude protein values to our predicted cover types, resulting in an estimated average carrying capacity of 760 moose (SD ± 428) across all sampling years, similar in magnitude to a density estimate of 716 moose (95% CI = 566–906) calculated during the same time. Our approach was able to accurately identify forest timber treatments across landscapes at differing spatial and temporal scales and provide an alternative method to estimate landscape-level ungulate carrying capacity. The ability to accurately identify areas of potential conflict from overbrowsing, or to highlight areas in need of land cover treatments can increase the toolset for ungulate management in managed forest landscapes.

1. Introduction

Estimating landscape carrying capacity provides a useful tool for the management of ungulates when evaluating the impacts of herbivory on regenerating forest stands [1]. Landscape carrying capacity is the number of ungulates that a landscape can nutritionally support, without individuals experiencing decreases in body condition [2,3,4]. Although managers have tried to estimate carrying capacity for forested stands in the past for moose [5,6], the estimation process is data intensive, and limited in resolution at large spatial scales. Remote sensing has been used to characterize forested landscapes and moose habitat suitability [7]; however, studies that use forage abundance data to estimate carrying capacity for the landscape have been limited [8]. Estimating landscape-level carrying capacity would allow managers to evaluate potential impacts of land management decisions on moose populations and assist in managing the complex relationship between forest regeneration and browsing populations.
Over the past twenty-five years, research has highlighted the importance of forest structure on ecosystem services and processes [9,10]. Concurrent advances in remote sensing allow researchers to map and model forest structure using satellite imagery [9,11,12,13]. Satellite-derived models have been used to estimate aboveground biomass [14,15], land cover type [16], forest structure [11], and successional stage [16,17,18], and evaluate changes in those metrics over time [19,20,21]. Satellite-derived data can similarly be utilized to monitor changes in the Normalized Difference Vegetation Index (NDVI), a modified ratio of near infrared to visible red wavelengths, to estimate “greenness” on the landscape [22]. Change in NDVI over time has been used as an index to measure the impacts of drought [23,24], assess forage availability for herbivores [25,26] and estimate agricultural production [27,28]. Most salient to our work, the difference in NDVI over time has been used to identify changes in forest plots due to harvesting activity [16,29].
Moose rely on patches of early successional habitat created by timber harvest for forage and seasonal thermoregulation [30,31]. The vegetation in early successional communities following timber harvest includes species that are preferred forage for moose and are abundant enough to support higher densities of moose compared to mature forests [32,33]. Access to abundant high-quality forage can improve body condition which is vital for species in areas of extreme winter temperatures and snow events [34,35]. However, preference for areas of timber harvest is not without nuance. Moose may avoid recently harvested areas due to continued harvest activity in neighboring stands, lack of vegetation immediately post-harvest, or compositional shifts in the early successional plant community towards non-preferred forage [36,37]. Despite these caveats, early successional habitat is sufficiently important for moose that it is often a primary influence on habitat selection [38]. Therefore, the ability to identify the distribution and quantity of early succession habitat on the landscape is essential for effective population management.
Here we present a method to estimate landscape carrying capacity by combining remote sensing classification of forest cover class with class-specific estimates of forage abundance from the literature. We initially developed the remote sensing method using the northern hardwood ecosystem of the Allegheny National Forest (ANFR) in northwestern Pennsylvania and then applied the method to the mixed temperate/boreal landscape of the Adirondack Park (AP) in northern New York to demonstrate the flexibility of the methodology. We then estimated landscape-level carrying capacity using the remote sensing predictions and data collected by Peterson [39] to demonstrate potential method applications. This relatively rapid and novel approach to estimate landscape-level carrying capacity may be useful for ungulate management where active forestry creates dynamic landscapes that can impact ungulate forage availability.

2. Materials and Methods

2.1. Study Site

2.1.1. Allegheny National Forest, Northwestern Pennsylvania

This study area was within bounds of the Allegheny National Forest Region (ANFR; Figure 1A), in northwestern Pennsylvania (Elk, Forest, McKean, and Warren counties). The ANFR covers 2077 km2 in the Allegheny Plateau ecoregion, which is dominated by secondary growth northern hardwood forest [40,41]. The region contains the Allegheny reservoir and is bordered to the west by the Allegheny River, and to the north by the state of New York. The woody overstory is predominately composed of sugar maple (Acer saccharum), red maple (A. rubrum), black cherry (Prunus serotina), and American beech (Fagus grandifolia). Areas of development are sparse and typically located at lower elevations across the landscape.
Agricultural production was limited (>7%) in the ANFR and typically occurs as small patches near towns or along roadways. The region has prominent hills and valleys with elevations ranging from 319 to 753 m [42]. The area receives an average of 107 cm of precipitation a year, including more than 150 cm of snowfall annually, and mean monthly average temperatures range from −4 °C in winter to 22 °C in summer [41,43,44]. The Allegheny region is well known for production of hardwood timber products [45]. Much of the landscape is subject to periodic timber treatments (e.g., overstory removals and shelterwood cuts; [46]). The distribution of timber treatments are dependent on ownership, stand age, species composition, and the state of neighboring stands, resulting in a dynamic heterogeneous landscape [47].

2.1.2. Adirondack Park, Northern New York

This study area was within the bounds of the Adirondack Park (24,281 km2; Figure 1B), in northern New York. Elevation in the Adirondack Park (AP) ranges from 100 m in low-lying lake shores to over 1600 m in the higher ranges. The park consists of large glacial valleys that gradually rise in elevation to the High Peaks region in the east-central part of the park. Average monthly temperatures range from −9 °C in winter to 18 °C in summer. The region receives an average of 100 cm of rainfall precipitation a year, with an additional 290 cm of snowfall annually [48].
The AP includes both publicly (61%) and privately owned land (39%), where all publicly owned lands are protected by Article XIV of the New York State Constitution as ‘forever wild forest’ which precludes resource extraction or development. The AP is comprised of a patchwork of the northern boreal ecosystem interspersed with temperate deciduous forests and large peatland complexes. Lower elevations with fertile soils support a diverse array of tree species dominated by American beech, yellow birch (Betula allegheniensis), paper birch (B. papyrifera), sugar maple and red maple. Higher elevations are more coniferous, dominated by species such as red spruce (Picea rubens), balsam fir (Abies balsamea), white pine (Pinus strobus), and eastern hemlock (Tsuga canadensis; [33,48]). Because public land acquisition continues within the AP, portions of public land could have been exposed to resource extraction or development immediately prior to acquisition by New York State.
Approximately 25% (13% of all AP land) of the forested private lands within the AP are participants with the New York State Conservation Easement Program. Private properties enrolled in an easement participate in a structured forest management program, which allows for timber harvest and other associated activities. Lands that are subjected to harvest are predominately composed of marketable timber species, such as sugar maple, red maple, red oak (Quercus rubra), white ash (Fraxinus americana), black cherry (Prunus serotina), and white pine [33]. Timber harvest methods included a combination of shelterwood removal, overstory removal, single tree selection, and salvage thinning.

2.2. Imagery Selection

2.2.1. Allegheny National Forest, Northwestern Pennsylvania

Landsat 5 and 8 satellite imagery (30 m resolution) was acquired from the USGS Global Visualization Viewer for imagery from March to October 1991–2013. We chose one cloudless image (<5% cloud cover) from the acquired Landsat imagery for each year during the focal period, if available. If multiple cloudless images were available within the same year, we selected an image that occurred during peak growing season (July-September) to reduce the seasonal variation associated with spring budding and fall leaf-off to maximize green reflectance. We did not select imagery from consecutive years due to an inability to acquire cloud-free imagery for older imagery at regular intervals. The resulting search yielded images for seven years during 1991–2013: 1991, 1993, 1996, 2000, 2002, 2006, 2009, and 2013 (Appendix A).

2.2.2. Adirondack Park, Northern New York

Landsat 8 satellite imagery (30 m resolution) was acquired annually from the USGS Global Visualization Viewer for imagery from late May to early October 2014–2018. Due to the size of the AP, we used 4 separate Landsat scenes per year (Figure 2). We chose one cloudless image (<5% cloud cover) per scene from the acquired Landsat imagery for each year during the focal period. If multiple cloudless images were available within the same year, we selected an image that occurred during peak growing season (July–September). The resulting data selection process yielded images for each of the five years for two of the four Landsat scenes. For the other two scenes, we only found cloudless images from 2015 to 2018 for one scene and images for 2013 and 2015–2018 for the second. Due to baseline reflectance variation across the four scenes within each year, we analyzed each of the four panels separately through time rather than creating one annual mosaicked landscape to increase prediction accuracy.

2.3. Remote Sensing Classification

We used random forests models [49,50,51] to classify each of the images across the time series into predictions of forest structure (similar to [16]) for both study sites. Covariates used in the random forests classifiers included NDVI [52,53], two shortwave infrared (SWIR) bands, and the previous predicted map (after 1991) (Table 1). NDVI provides an index of green biomass across a landscape. The use of NDVI helps discern changes in canopy greenness associated with timber harvest.
NDVI = ( near infrared ) ( red ) ( near infrared ) + ( red )
We also used the two shortwave infrared (SWIR) bands (approximately 1.55–1.75 mm and 2.09–2.35 mm) and took the ratio of each SWIR band divided by the near infrared (NIR) band to reduce the noise from shadows and aerosols. SWIR reflectance can aide in the detection of soil differences following forest harvest activities [54]. We additionally calculated the differences in NDVI and the two SWIR ratios between a given year and the previous time step. Lastly, we included the predicted cover class raster from the previous prediction as a predictor variable in each subsequent model (e.g., cover prediction from 2013 was used in the 2014 model prediction for the AP), using the first panel to create a baseline that would allow for future scenario planning rather than historical evaluations.
The training data for each time period consisted of training polygons derived from landowner data that was provided by two regional timber companies (Lyme Timber Co. and Kane Hardwoods). The data consisted of all timber treatments conducted during and immediately prior to the study. Data for shelterwood and overstory removals were ground-truthed using NAIP Imagery and classified into three forest cover classes: mature (<30% canopy openness, typically no removal), intermediate removal (~30–60% removal), and overstory removal (>60% removal). To ensure adequate representation of the variation within each classification type we used a minimum of 75 polygons representing each of the 3 forest classes. This resulted in 245–270 (x = 260) training polygons per panel for the ANFR study and 408–1556 (x = 990) polygons per panel for the AP Study). A mean value was calculated for each training polygon for each predictor raster (NDVI, two SWIR bands, the 3 difference bands, and previous year’s prediction, if applicable). The polygon means were aggregated within the respective training polygon classification to create a range of values per classification. For each year, we used a random forest model with 5000 trees, a node size = 2, and proportional sampling of training data where each tree was trained with 70% of the available data for each cover class and the remaining 30% was withheld for out-of-bag testing. After each model was fit, we generated predicted cover classes back out to the landscape as raster values. We measured the accuracy of the model by comparing the value of the training polygon, from each respective sampling year, to the predicted majority class within each polygon and calculated user’s (probability that a polygon is included to an incorrect classification type), producer’s (probability that a polygon is excluded from the correct classification type), and overall accuracy (correctly classified polygons/total number of polygons) using a confusion matrix [55]. Lastly, we calculate the Kappa value to evaluate the degree of agreement between our prediction and on the ground training polygons [56].

2.4. Estimation of Landscape-Level Carrying Capacity for Moose

Following the creation of the remote sensed AP map, we incorporated additional land cover layers for non-predicted cover types to create a more accurate depiction of the landscape. First, we removed areas of state-owned public lands acquired ≥ 20 years ago from the prediction raster since those lands were not subject to recent timber treatments. All public forested stands were assumed to be mature given that any lands cut prior to public acquisition have since matured. Additionally, this reduced the amount of model misclassification, since the majority of the high elevation ‘high peaks’ regions are older public lands. We used NLCD 2016 data for agriculture, developed, conifer, grass, and water classes [57]. We reclassified each of the NLCD datasets into one of seven land cover types to simplify land cover at a large spatial scale (Table 2). In addition to the seven land cover types used in the ANFR, we used a wetland layer developed by the Adirondack Park Agency to ensure that wetlands were adequately identified, given their significance to moose [58].
We merged the simplified NLCD layer with the AP wetlands layer. The AP wetlands data were higher resolution and derived with more focused methodology, therefore we prioritized the AP wetlands layer; retaining AP wetlands over NLCD data when both were present. We combined the AP wetlands/NLCD layer with four of the non-public lands prediction years (2015–2018) separately into their own annual maps, prioritizing the timber harvest prediction layers over the AP wetlands/NLCD layer in cases where NLCD land cover was classified as deciduous forest, mixed forest, conifer forest and scrubland. Following the merging and classification of the map layers, we then applied known values of animal use days (AUD) for lactating female moose to each cover class [39]. We only derived AUD using summer crude protein nutritional estimate per cover type, as it was determined to be the limiting nutritional value for moose in the Adirondack region of New York (Table 3; [39]). To match the landscape classes of Peterson [39], we used a digital elevation model to differentiate between lowland and upland deciduous and mixed forest at an elevation of 497 m.
The carrying capacity values for the land cover classes of agriculture, developed/barren, and grassland were not calculated in Peterson [39], and the three cover classes combined accounted for less than five percent of the AP. Additionally, the limited agriculture production in the region focused on hay and corn (1.2% of the landscape), and given other previous works on moose browsing and forage availability [59], we deemed it appropriate to assign those carrying capacity values equal to zero. Lastly, we combined intermediate removal and overstory removal forest cover types to align with the classification of ‘Regenerating Forest’ used in Peterson [39]. We estimated the number of AUD for the Adirondacks for each of the four prediction years by applying the cover capacity estimates to the proportion of each cover type on the landscape and calculated the average AUD value across years. We did not apply a landscape-level moving window because a park-wide estimate of capacity was required, rather than site specific values of capacity. Lastly, we compared our estimates for AUD to moose density estimates calculated using aerial distance sampling during the same time period in the Adirondacks [60].

3. Results

3.1. Allegheny National Forest, Northwestern Pennsylvania

Our three-class prediction (mature forest, intermediate removal and overstory removal) for each of the selected years resulted in overall prediction accuracy of 94.3% across all years (Table 4; Appendix B). The average class producer’s accuracy and user’s accuracy were high across all years (70.6–100%), with intermediate removal cover type having the lowest average user’s accuracy (x = 81.4) and overstory removal had the lowest producer’s accuracy (x = 88.6). The majority of the classification confusion occurred between the overstory removal and intermediate removal classes, while mature forest classification had a user’s accuracy of 100% in all but one (2009) of the seven predictive land covers. The average Cohen’s kappa value across all years and panels suggest that the predictions were in agreement (Khat = 0.91) with the reality on the ground. The variable importance plots indicated that NDVI was the most important predictor for land cover in a majority of the predictions (n = 5). Following the merging with NLCD spatial layers, the average land cover proportion for mature forest was 70% across all years (Table 5A). The two timber harvest prediction layers, intermediate removal and overstory removal, made up 2.9% and 6.2% of the landscape.

3.2. Adirondack Park, Northern New York

Our predictive land cover for mature forest, intermediate removal and overstory removal in the AP yielded a mean overall accuracy of 86.9% across all years/panels. The average class producer’s accuracy was high across all years and panels, ranging from 80.7 to 91.7%. However, the average user’s accuracy for the three predictive classes of mature forest, intermediate removal and overstory removal were 92.5%, 84.7%, and 62.8% respectively. The average Cohen’s kappa value across all years and panels suggest that the predictions were in ‘substantial’ agreement (Khat = 0.76) with the reality on the ground [56]. The variable importance plots indicated that SWIR bands were the most important predictors for land cover in a majority of the predictions (n = 15).
After combining the four spatial layers (prediction, DEC lands, NLCD, and AP Wetlands) and reclassifying to match Peterson [39] cover classes, regenerating forest and upland deciduous/mixed forest accounted for an average 10.6% and 31.9% of the landscape (Table 5B), respectively. The average AUD value across the four years was 759.7 (±427.8 SD; Figure 3). The AUD estimate was lowest in 2015 (AUD = 679.4) and highest for the 2017 prediction (AUD = 868.6). Additionally, the proportion of predicted regenerating forest was the lowest in 2015 (8.9%) and greatest in 2017 (12.9%). Our estimated moose population in the AP was of similar magnitude as the estimate from aerial distance sampling of 716 (95% CI = 566–906) moose within the Adirondack Park [60].

4. Discussion

For ungulates that derive the bulk of their necessary forage from only a few land cover types, silvicultural treatments that alter and create those types can inordinately affect forage availability on the landscape. However, assessment of those changes is not possible from general land cover classifications (e.g., NLCD) due to the lack of specific land cover predictions within the national land cover layers [57]. Additionally, compiling silvicultural treatment data at a landscape scale can be difficult due to data variation among multiple stakeholders and the complex public-private nature of timber harvest networks. Here, we show that it is possible to create and update land cover maps that accurately predict the spatial distribution in land cover types in ways that meaningfully relate to moose carrying capacity at large spatial scales. This work highlights the utility of combining multiple forms of remote sensing data and on the ground data collection to (1) refine our understanding of habitat condition by better classifying forest cover into biologically relevant sub-categories; and (2) by incorporating known estimates of forest cover type carrying capacity to estimate landscape-scale carrying capacity. Our method provided a useful tool that can aid in the management of ungulates in dynamic forested landscapes by giving managers the ability to evaluate the current land cover capacity for ungulate populations and to assess long-term impacts by identifying patches of landscape that may transition from one cover class to another (i.e., harvest of a mature stand, significantly increasing capacity, or transitioning from early successional stands into mature forest, significantly decreasing capacity).
Our work emphasizes the role that timber treatments play in influencing landscape-scale carrying capacity for ungulates, and that active engagement with forest management companies and private landowners can produce high-resolution land cover identification. By engaging corporate and private landowners and pooling data across sources, more accurate landscape maps can be developed to help aid in wildlife management and research [37]. Providing managers with the ability to accurately evaluate current land cover composition can help to provide insight on the future trajectory of ungulate populations, thereby helping to determine appropriate actions to meet management objectives. While our work focuses on anthropogenic-driven forest change, the methodology could be applied to blow down, fire or other large-scale forest landscape change and determine how those impacts may influence ungulate populations. Lastly, our method bridges the gap between landowner-driven localized forest management and large-scale population management objectives developed by state agencies.
Our timber harvest modeling approach was applied to two study sites that varied in spatial extent and temporal scale, but yielded predictions that were fairly accurate despite those differences. There was greater prediction variation within the AP, likely because of the AP’s land area is ten times greater than the ANFR and contains greater variation in forest types and elevation. The temporal differences between the studies, ranging from annual predictions to up three-year intervals, did not reduce our ability to identify anthropogenic forest change on the landscape or negatively influence prediction accuracy. Despite the greater variation at the larger spatial scale, the demonstrated predictive ability at multiple spatial scales makes our classification approach and associated carrying capacity estimates a viable tool for the management of large ungulates, as large ungulates are managed at multiple spatial scales, from localized impacts to larger scale management units [61]. Additionally, the temporal flexibility allows prediction of the impacts of recent timber harvest on an annual basis or to periodically evaluate landscape changes when setting long-term management objectives [32]. Recent approaches have utilized LiDAR to create more comprehensive forest cover predictions that can identify more subtle variations in stand structure [51,62,63], but LiDAR can be cost prohibitive when predicting forest cover at large spatial scales or at multiple points in time [12].

5. Conclusions

In both of our case studies, only two or three broadly classified land cover types contributed to the carrying capacity for moose at larger spatial scales. While localized forage abundance may impact individual survival and fecundity [64,65], fine-scale variations may be nearly undetectable at the large spatial scales at which large ungulates are managed. It is even possible that a single, potentially rare land cover class (e.g., regenerating forest) can influence carrying capacity on a massive landscape as we identified in the case of moose in the Adirondacks. This can result in substantial shifts in landscape-level carrying capacity in very short-order in dynamic predominately hardwood forests that is controlled by anthropogenic change (i.e., timber harvest regimes). The role of one or two cover classes highlights the need for managers to regularly evaluate the current landscape-level capacity and determine how future cover change can impact large ungulates. Many state agencies tend to regulate ungulates by setting and altering harvest tag allocation, often based on current ungulate browsing impacts. By using remote sensing to predict current forest condition and evaluate potential changes in ungulate carrying capacity, agencies can increase or decrease harvest to reflect anticipated carrying capacity, rather than managing for impacts post-hoc, when landscape change has already been significant and/or animal population size or health have already responded to change.

Author Contributions

Conceptualization, D.W.K., K.V.M., A.A.R. and N.P.N.; methodology, D.W.K. and T.J.P.; writing—original draft preparation, D.W.K.; writing—review and editing, D.W.K., A.A.R., N.P.N., K.V.M. and J.L.F.; funding acquisition, A.A.R., N.P.N., K.V.M. and J.L.F. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the USDA Forest Service Northern Research Station, the University of Georgia Warnell School of Forestry and Natural Resources (USDA-AFRI Award #12-IA-11242302-093), SUNY-ESF, and NYS-DEC (Federal Aid in Wildlife Restoration Grant W-173-G).

Institutional Review Board Statement

No institutional review was required to complete this research.

Informed Consent Statement

This research did not require informed consent.

Data Availability Statement

Landsat data can be found at https://earthexplorer.usgs.gov/. Data used for training polygons are private property of The Forestland Group, Generations Forestry, Hancock Forest Management, Landvest, Lyme Adirondacks, Kane Hardwoods and Molpus Woodlands and not publically available.

Acknowledgments

This project would not have been possible without the cooperation of numerous commercial and public foresters including Allegheny National Forest, the Pennsylvania Bureau of Forestry, Bradford Water Authority, Forest Investment Associates, The Forestland Group, Generations Forestry, Hancock Forest Management, Landvest, Lyme Adirondacks, Kane Hardwoods and Molpus Woodlands. We thank countless other technicians, consultants, and volunteers who conducted annual pellet surveys and moose browse surveys that were essential for the completion of this research. We are grateful for the review and comments provided by J.E. Hurst. Landsat data are available from the U.S. Geological Survey. Questions concerning the use or redistribution of USGS data should be directed to: [email protected] or 1-888-ASK-USGS (1-888-275-8747).

Conflicts of Interest

The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A

List of Landsat 5 and 8 scene images used in a time series regression to estimate forest classes across the Allegheny National Forest, Pennsylvania and the Adirondack Park, New York.
PanelYearDateScene
Case Study 1: Allegheny National Forest, Pennsylvania
01703119917-AprLT50170311991097XXX03
017031199330-MayLT50170311993150PAC03
017031199626-AugLT50170311996239XXX01
01703120006-SepLT50170312000250XXX03
017031200212-SepLT50170312002255LGS01
01703120069-OctLT50170312006282GNC01
017031200923-MarLT50170312009082GNC01
017031201326-SepLC80170312013269
Case Study 2: Adirondack Park, New York
01402920148-SepLC08_L1TP_014029_20140908_20170303_01_T1
014029201527-SepLC08_L1TP_014029_20150927_20170225_01_T1
014029201613-SepLC08_L1TP_014029_20160913_20180130_01_T1
014029201730-JulLC08_L1TP_014029_20170730_20170811_01_T1
01402920185-OctLC08_L1TP_014029_20181005_20181010_01_T1
01403020146-JulLC08_L1TP_014030_20140706_20170304_01_T1
014030201527-SepLC08_L1TP_014030_20150927_20170225_01_T1
014030201613-SepLC08_L1TP_014030_20160913_20180130_01_T1
014030201730-JulLC08_L1TP_014030_20170730_20170811_01_T1
01403020185-OctLC08_L1TP_014030_20181005_20181010_01_T1
015029201516-JulLC08_L1TP_015029_20150716_20170226_01_T1
01502920164-SepLC08_L1TP_015029_20160904_20170221_01_T1
015029201723-SepLC08_L1TP_015029_20170923_20171013_01_T1
015029201822-JunLC08_L1TP_015029_20180622_20180703_01_T1
015030201328-SepLC08_L1TP_015030_20130928_20170308_01_T1
015030201516-JulLC08_L1TP_015030_20150716_20170226_01_T1
015030201620-SepLC08_L1TP_015030_20160920_20170221_01_T1
015030201723-SepLC08_L1TP_015030_20170923_20171013_01_T1
015030201821-MayLC08_L1TP_015030_20180521_20180605_01_T1

Appendix B

Classification error matrix for the three-class random forest prediction model for each of the sampling periods in northwestern Pennsylvania from 1993–2013 (A) and for each of the four panels and sampling periods in northern New York from 2013–2018 (B).
A.Pennsylvania
Reference
1993MatureOverstory RemovalIntermediate RemovalUser’s Accuracy
ClassificationMature107001.00
Overstory Removal011380.93
Intermediate Removal30190.86
Producer’s Accuracy0.971.000.71
Overall Accuracy (%)95.5
Kappa Index0.92
Reference
1996MatureOverstory RemovalIntermediate RemovalUser’s Accuracy
ClassificationMature100001.00
Overstory Removal011380.93
Intermediate Removal53200.71
Producer’s Accuracy0.950.970.71
Overall Accuracy (%)93.3
Kappa Index0.89
Reference
2000MatureOverstory RemovalIntermediate RemovalUser’s Accuracy
ClassificationMature98001.00
Overstory Removal09950.95
Intermediate Removal20460.97
Producer’s Accuracy0.981.000.91
Overall Accuracy (%)97.4
Kappa Index0.96
Reference
2002MatureOverstory RemovalIntermediate RemovalUser’s Accuracy
ClassificationMature89001.00
Overstory Removal289140.85
Intermediate Removal00571.00
Producer’s Accuracy0.981.000.80
Overall Accuracy (%)93.7
Kappa Index0.9
Reference
2006MatureOverstory RemovalIntermediate RemovalUser’s Accuracy
ClassificationMature94001.00
Overstory Removal08320.98
Intermediate Removal84590.83
Producer’s Accuracy0.920.950.97
Overall Accuracy (%)94.5
Kappa Index0.92
Reference
2009MatureOverstory RemovalIntermediate RemovalUser’s Accuracy
ClassificationMature121010.99
Overstory Removal149150.75
Intermediate Removal31580.93
Producer’s Accuracy0.970.970.78
Overall Accuracy (%)91.2
Kappa Index0.86
Reference
2013MatureOverstory RemovalIntermediate RemovalUser’s Accuracy
ClassificationMature109001.00
Overstory Removal056130.81
Intermediate Removal13680.94
Producer’s Accuracy0.990.950.84
Overall Accuracy (%)93.1
Kappa Index0.89
B.New York
Reference
Panel 014029—2014MatureOverstory RemovalIntermediate RemovalUser’s Accuracy
ClassificationMature2672130.95
Overstory Removal2341380.84
Intermediate Removal12142510.69
Producer’s Accuracy0.880.960.68
Overall Accuracy (%)83.9
Kappa Index0.75
Reference
Panel 014029—2015MatureOverstory RemovalIntermediate RemovalUser’s Accuracy
ClassificationMature3220130.96
Overstory Removal4844220.89
Intermediate Removal057380.6
Producer’s Accuracy0.870.930.89
Overall Accuracy (%)88.6
Kappa Index0.8
Reference
Panel 014029—2016MatureOverstory RemovalIntermediate RemovalUser’s Accuracy
ClassificationMature3120230.93
Overstory Removal3634350.92
Intermediate Removal060350.63
Producer’s Accuracy0.90.950.88
Overall Accuracy (%)89.3
Kappa Index0.81
Reference
Panel 014029—2017MatureOverstory RemovalIntermediate RemovalUser’s Accuracy
ClassificationMature3013310.9
Overstory Removal3064380.92
Intermediate Removal065300.68
Producer’s Accuracy0.910.880.88
Overall Accuracy (%)88.9
Kappa Index0.80
Reference
Panel 014029—2018MatureOverstory RemovalIntermediate RemovalUser’s Accuracy
ClassificationMature3363260.92
Overstory Removal3433970.91
Intermediate Removal573270.7
Producer’s Accuracy0.90.920.88
Overall Accuracy (%)89.2
Kappa Index0.81
Reference
Panel 014030—2014MatureOverstory RemovalIntermediate RemovalUser’s Accuracy
ClassificationMature3350220.94
Overstory Removal5246250.92
Intermediate Removal0105700.6
Producer’s Accuracy0.870.960.87
Overall Accuracy (%)87.8
Kappa Index0.78
Reference
Panel 014030—2015MatureOverstory RemovalIntermediate RemovalUser’s Accuracy
ClassificationMature5540320.82
Overstory Removal7753630.95
Intermediate Removal074750.5
Producer’s Accuracy0.770.940.88
Overall Accuracy (%)84.0
Kappa Index0.72
Reference
Panel 014030—2016MatureOverstory RemovalIntermediate RemovalUser’s Accuracy
ClassificationMature5970330.95
Overstory Removal8333360.8
Intermediate Removal051770.4
Producer’s Accuracy0.880.940.75
Overall Accuracy (%)83.4
Kappa Index0.70
Reference
Panel 014030—2017MatureOverstory RemovalIntermediate RemovalUser’s Accuracy
ClassificationMature5860440.93
Overstory Removal7923410.81
Intermediate Removal049790.38
Producer’s Accuracy0.880.960.73
Overall Accuracy (%)82.7
Kappa Index0.69
Reference
Panel 014030—2018MatureOverstory RemovalIntermediate RemovalUser’s Accuracy
ClassificationMature9262120.99
Overstory Removal7317570.39
Intermediate Removal562260.67
Producer’s Accuracy0.920.770.6
Overall Accuracy (%)88.6
Kappa Index0.63
Reference
Panel 015029—2015MatureOverstory RemovalIntermediate RemovalUser’s Accuracy
ClassificationMature6840370.95
Overstory Removal63144670.86
Intermediate Removal1216740.74
Producer’s Accuracy0.910.960.81
Overall Accuracy (%)87.9
Kappa Index0.80
Reference
Panel 015029—2016MatureOverstory RemovalIntermediate RemovalUser’s Accuracy
ClassificationMature6800410.94
Overstory Removal49134820.89
Intermediate Removal1216740.74
Producer’s Accuracy0.930.940.81
Overall Accuracy (%)88.6
Kappa Index0.82
Reference
Panel 015029—2017MatureOverstory RemovalIntermediate RemovalUser’s Accuracy
ClassificationMature7540470.94
Overstory Removal61263790.81
Intermediate Removal5208760.72
Producer’s Accuracy0.920.890.75
Overall Accuracy (%)86.2
Kappa Index0.77
Reference
Panel 015029—2018MatureOverstory RemovalIntermediate RemovalUser’s Accuracy
ClassificationMature8964520.94
Overstory Removal68232520.73
Intermediate Removal15190560.73
Producer’s Accuracy0.920.880.7
Overall Accuracy (%)86.0
Kappa Index0.74
Reference
Panel 015030—2013MatureOverstory RemovalIntermediate RemovalUser’s Accuracy
ClassificationMature1110270.8
Overstory Removal1941920.89
Intermediate Removal014510.25
Producer’s Accuracy0.850.780.74
Overall Accuracy (%)75.8
Kappa Index0.57
Reference
Panel 015030—2015MatureOverstory RemovalIntermediate RemovalUser’s Accuracy
ClassificationMature1530180.89
Overstory Removal1222650.95
Intermediate Removal029110.72
Producer’s Accuracy0.930.940.9
Overall Accuracy (%)91.2
Kappa Index0.84
Reference
Panel 015030—2016MatureOverstory RemovalIntermediate RemovalUser’s Accuracy
ClassificationMature2440130.95
Overstory Removal2041690.88
Intermediate Removal029110.72
Producer’s Accuracy0.920.880.88
Overall Accuracy (%)90.2
Kappa Index0.82
Reference
Panel 015030—2017MatureOverstory RemovalIntermediate RemovalUser’s Accuracy
ClassificationMature2150160.93
Overstory Removal1241950.92
Intermediate Removal034140.71
Producer’s Accuracy0.950.890.87
Overall Accuracy (%)90.6
Kappa Index0.84
Reference
Panel 015030—2018MatureOverstory RemovalIntermediate RemovalUser’s Accuracy
ClassificationMature2380130.95
Overstory Removal2381410.82
Intermediate Removal150160.75
Producer’s Accuracy0.910.860.83
Overall Accuracy (%)87.6
Kappa Index0.79

References

  1. De Calesta, D.S.; Stout, S.L. Relative deer density and sustainability: A conceptual framework for integrating deer management with ecosystem management. Wildl. Soc. Bull. 1997, 25, 252–285. [Google Scholar]
  2. McCullough, D.R. The George Reserve Deer Herd: Population Ecology of a K-Selected Species; University of Michigan Press: Ann Arbor, MI, USA, 1979; p. 271. [Google Scholar]
  3. Macnab, J. Carrying capacity and related slippery shibboleths. Wildl. Soc. Bull. 1985, 13, 403–410. [Google Scholar]
  4. Hobbs, N.T.; Swift, D.M. Estimates of habitat carrying capacity incorporating explicit nutritional constraints. J. Wildl. Manag. 1985, 49, 814–822. [Google Scholar] [CrossRef]
  5. Crête, M. Approximation of K carrying capacity for moose in eastern Quebec. Can. J. Zool. 1989, 67, 373–380. [Google Scholar] [CrossRef]
  6. Wam, K.H.; Hjeljord, O.; Soldberg, E.J. Differential forage use makes carrying capacity equivocal on ranges of Scandanavian moose (Alces alces). Can. J. Zool. 2010, 88, 1179–1191. [Google Scholar] [CrossRef]
  7. Razenkova, E.; Radeloff, V.; Dubinin, M.; Bragina, E.; Allen, A.; Clayton, M.; Pidgeon, A.; Baskin, L.; Coops, N.; Hobi, M. Vegetation productivity summarized by the Dynamic Habitat Indices explains broad-scale patterns of moose abundance across Russia. Sci. Rep. 2020, 10, 836. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Doan, T.; Guo, X. Understanding bison carrying capacity estimation in Northern Great Plains using remote sensing and GIS. Can. J. Remote Sens. 2019, 45, 139–162. [Google Scholar] [CrossRef]
  9. O’Hara, L.O.; Latham, P.A.; Hessburg, P.; Smith, B.G. A structural classification for inland Northwest forest vegetation. West. J. Appl. For. 1996, 11, 97–102. [Google Scholar] [CrossRef] [Green Version]
  10. Brockerhoff, E.G.; Barbaro, L.; Castagneyrol, B.; Forrester, D.; Gardiner, B.; González-Olabarria, J.; Lyver, P.; Meurisse, N.; Oxbrough, A.; Taki, H.; et al. Forest biodiversity, ecosystem functioning and the provision of ecosystem services. Biodivers. Conserv. 2017, 26, 3005–3035. [Google Scholar] [CrossRef] [Green Version]
  11. Hummel, S.; Hudak, A.; Uebler, E.; Falkowski, M.; Megown, K. A comparison of accuracy and cost of LiDAR versus stand exam data for landscape management on the Malhuer National Forest. J. For. 2011, 109, 267–273. [Google Scholar]
  12. Beland, M.; Parker, G.; Sparrow, B.; Harding, D.; Chasmer, L.; Phinn, S.; Antonarakis, A.; Strahler, A. On promoting the use of LiDAR systems in forest ecosystem research. For. Ecol. Manag. 2019, 450, 117484. [Google Scholar] [CrossRef]
  13. LaRue, E.A.; Wagner, F.; Fei, S.; Atkins, J.; Fahey, R.; Gough, C.; Hardiman, B.B. Compatibility of aerial and terrestrial LiDAR for quantifying forest structural diversity. Remote Sens. 2020, 12, 1407. [Google Scholar] [CrossRef]
  14. Coops, N.C.; Hilker, T.; Wulder, M.A.; St-Onge, B.; Newnham, G.; Siggins, A.; Trofymow. J.A. Estimating canopy structure of Douglas-fir forests from discrete-return LiDAR. Trees 2007, 21, 295–310. [Google Scholar] [CrossRef] [Green Version]
  15. Clark, D.B.; Olivas, P.C.; Oberbauer, S.F.; Clark, D.A.; Ryan, M.G. First direct landscape-scale measurement of tropical rain forest leaf area index, a key driver of global primary productivity. Ecol. Lett. 2008, 11, 163–172. [Google Scholar] [CrossRef]
  16. Falkowski, M.J.; Evans, J.S.; Martinuzzi, S.; Gessler, P.E.; Hudak, A.T. Characterizing forest succession with LIDAR data: An evaluation for Inland Northwest, USA. Remote Sens. Environ. 2009, 113, 946–956. [Google Scholar] [CrossRef] [Green Version]
  17. Bergen, K.M.; Dronova, I. Observing succession on aspen-dominated landscape using remote sensing-ecosystem approach. Landsc. Ecol. 2007, 22, 1395–1410. [Google Scholar] [CrossRef]
  18. Hao, Z.; Xhang, J.; Song, B.; Ye, J.; Li, B. Vertical structure and spatial associations of dominant tree species in an old-growth temperature forest. For. Ecol. Manag. 2007, 252, 1–11. [Google Scholar] [CrossRef]
  19. Baker, W.L. A review of models of landscape change. Landsc. Ecol. 1989, 1, 111–133. [Google Scholar] [CrossRef]
  20. Shugart, H.H. The importance of structure in the longer-term dynamics of ecosystems. J. Geophys. Res–Atmos. 2000, 105, 20065–20075. [Google Scholar] [CrossRef]
  21. Homer, C.G.; Dewitz, J.; Yang, L.; Jin, S.; Danielson, P.; Xian, G.; Coulston, J.; Herold, N.; Wickham, J.; Megown, K. Completion of the 2011 National Land Cover Database for the conterminous United States-Representing a decade of land cover change information. Photogramm. Eng. Remote Sens. 2015, 81, 345–354. [Google Scholar]
  22. Madden, M. (Ed.) Manual of Geographic Information Systems; American Society for Photogrammetry and Remote Sensing: Bethesda, MD USA, 2009; p. 1330. [Google Scholar]
  23. Peters, A.J.; Walter-Shea, E.A.; Ji, L.; Vina, A.; Hayes, M.; Svoboda, M.D. Drought monitoring with NDVI-based standardized vegetation index. Photogramm. Eng. Remove Sens. 2002, 68, 71–75. [Google Scholar]
  24. Gu, Y.; Brown, J.F.; Verdin, J.P.; Wardlow, B. A five-year analysis of MODIS NDVI and NDWI for grassland drought assessment over the central Great Plains of the United States. Geophys. Res. Lett. 2007, 34, 6. [Google Scholar] [CrossRef] [Green Version]
  25. Hebblewhite, M.; Merrill, E.; McDermid, G. A multi-scale test of the forage maturation hypothesis in a partially migratory ungulate population. Ecol. Monogr. 2008, 78, 141–166. [Google Scholar] [CrossRef] [Green Version]
  26. Hansen, B.B.; Aanes, R.; Herfindal, I.; Saether, B.; Henriksen, S. Winter habitat–space use in a large arctic herbivore facing contrasting forage abundance. Polar Biol. 2009, 32, 971–984. [Google Scholar] [CrossRef]
  27. Quarmby, N.A.; Milnes, M.; Hindle, T.L.; Silleos, N. The use of multi-temporal NDVI measurements from AVHRR data for crop yield estimation and prediction. Int. J. Remote Sens. 1993, 14, 199–210. [Google Scholar] [CrossRef]
  28. Hayes, M.J.; Decker, W.L. Using NOAA AVHRR data to estimate maize production in the United States Corn Belt. Int. J. Remote Sens. 1996, 17, 3189–3200. [Google Scholar] [CrossRef]
  29. Wilson, E.H.; Sader, S.A. Detection of forest harvest type using multiple dates of Landsat TM imagery. Remote Sens. Environ. 2002, 80, 385–396. [Google Scholar] [CrossRef]
  30. Renecker, L.A.; Hudson, R.J. Seasonal energy expenditures and thermoregulatory responses of moose. Can. J. Zool. 1986, 64, 322–327. [Google Scholar] [CrossRef]
  31. Fisher, J.T.; Wilkinson, L. The response of mammals to forest fire and timber harvest in the North America boreal forest. Mammal Rev. 2005, 35, 51–81. [Google Scholar] [CrossRef]
  32. Schrempp, T.V.; Rachlow, J.; Johnson, R.; Shipley, L.; Long, R.; Aycrigg, J.; Hurley, M. Linking forest management to moose population trends: The role of the nutritional landscape. PLoS ONE 2019, 14, e0219128. [Google Scholar] [CrossRef] [Green Version]
  33. Peterson, S.; Kramer, D.; Hurst, J.; Frair, J. Browse selection by moose in the Adirondack Park, New York. Alces 2020, 56, 107–126. [Google Scholar]
  34. Regelin, W.L.; Schwartz, C.; Franzmann, A. Effects of forest succession on nutritional dynamics of moose forage. Swed. Wildl. Res. Suppl. 1987, 1, 247–264. [Google Scholar]
  35. Saether, B.E.; Andersen, R. Resource limitation in a generalist herbivore, the moose, (Alces alces): Ecological constraints on behavioral decisions. Can. J. Zool. 1990, 68, 993–999. [Google Scholar] [CrossRef]
  36. Milligan, H.R.; Koricheva, J. Effects of tree species richness and composition on moose winter browsing damage and foraging selectivity: An experimental study. J. Anim. Ecol. 2013, 82, 739–748. [Google Scholar] [CrossRef]
  37. Mumma, M.A.; Gillingham, M.; Marshall, S.; Proctor, C.; Bevington, A.; Scheideman, M. Regional moose (Alces alces) responses to forestry cutblocks are driven by landscape-scale patterns of vegetation composition and regrowth. For. Ecol. Manag. 2021, 481, 118763. [Google Scholar] [CrossRef]
  38. Thompson, I.D.; Stewart, R.W. Management of Moose Habitat. In Ecology and Management of the North American Moose, 2nd ed.; Schwartz, C., Franzmann, A., McCabe, R., Eds.; Smithsonian Institution Press: Washington, DC, USA, 1997; pp. 377–401. [Google Scholar]
  39. Peterson, S. Browse Selection and Constraints for Moose (Alces alces) in the Adirondack Park, New York, USA. Master’s Thesis, State University of New York, College of Environmental Science and Forestry, Syracuse, New York, NY, USA, 2018. [Google Scholar]
  40. Hough, A.; Forbes, R. The ecology and silvics of forests in the high plateau of Pennsylvania. Ecol. Monogr. 1943, 13, 299–320. [Google Scholar] [CrossRef]
  41. Anacker, B.L.; Kirschbaum, C.D. Vascular flora of the Kinzua Quality Deer Cooperative, northwestern Pennsylvania, U.S.A. Bartonia 2006, 63, 11–28. [Google Scholar]
  42. Redding, J. History of Deer Population Trends and Forest Cutting on the Allegheny National Forest; General Technical Report NE-197; U.S. Department of Agriculture, Forest Service, Northeastern Forest Experiment Station: Warren, PA, USA, 1995.
  43. Bjorkbom, J.C.; Larson, R.G. The Tionesta Scenic and Research Natural Areas; General Technical Report NE-31; USDA Forest Service, Northeastern Forest Experiment Station: Warren, PA, USA, 1977.
  44. Menne, M.J.; Williams, C.N.; Vose, R.S. The United States historical climatology network monthly temperature data-version 2. Bull. Am. Meteorol. Soc. 2009, 90, 993–1107. [Google Scholar] [CrossRef] [Green Version]
  45. Lumber Heritage Region of Pennsylvania. Management Action Plan, May 2001; Lumber Heritage Region: Emporium, PA, USA, 2001; p. 86. [Google Scholar]
  46. Nyland, R.D. Silviculture: Concepts and Applications; McGraw-Hill Higher Education: Burr Ridge, IL, USA, 2002. [Google Scholar]
  47. Marquis, D.; Ernst, R.; Stout, S. Prescribing Silvicultural Treatments in Hardwood Stands of the Alleghenies (Revised); General Technical Report NE-96; USDA Forest Service, Northeastern Forest Experiment Station: Radnor, PA, USA, 1992.
  48. Jenkins, J.; Keal, A. The Adirondack Atlas: A Geographic Portrait of the Adirondack Park; Syracuse University Press: Syracuse, NY, USA, 2004. [Google Scholar]
  49. Breiman, L. Random forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef] [Green Version]
  50. Liaw, A.; Wiener, M. Classification and regression by Random Forest. R News 2002, 2, 18–22. [Google Scholar]
  51. Shanley, C.S.; Eacker, D.; Reynolds, C.; Bennetsen, B.; Gilbert, S. Using LiDAR and Random Forest to improve deer habitat models in a managed forest landscape. For. Ecol. Manag. 2021, 499, 119580. [Google Scholar] [CrossRef]
  52. Chavez, P.S., Jr.; MacKinnon, D.J. Automatic detection of vegetation changes in the Southwestern United States using remotely sensed images. Photogramm. Eng. Remote Sens. 1994, 60, 571–583. [Google Scholar]
  53. Eidenshink, J.C.; Faudeen, J.L. The 1-km AVHRR global land data set: First stages in implementation. Int. J. Remote Sens. 1994, 15, 3443–3462. [Google Scholar] [CrossRef]
  54. Drake, N.A.; Mackin, S.; Settle, J.J. Mapping vegetation, soils, and geology in semiarid shrublands using spectral matching and mixture modeling of SWIR AVIRIS imagery. Remove Sens. Environ. 1999, 68, 12–25. [Google Scholar] [CrossRef]
  55. Lewis, H.G.; Brown, M. A generalized confusion matrix for assessing area estimate from remotely sensed data. Int. J. Remote Sens. 2001, 22, 3223–3235. [Google Scholar] [CrossRef]
  56. Landis, J.R.; Koch, G. An application of hierarchical kappa-type statistics in the assessment of majority agreement among multiple observers. Biometrics 1977, 33, 363–374. [Google Scholar] [CrossRef] [PubMed]
  57. Homer, C.G.; Dewitz, J.; Jin, S.; Xian, G.; Costello, C.; Danielson, P.; Gass, L.; Funk, M.; Wickham, J.; Stehman, S.; et al. Conterminous United States land cover change patterns 2001–2016 from the 2016 National Land Cover Database. ISPRS J. Photogramm. Remote Sens. 2020, 162, 184–199. [Google Scholar] [CrossRef]
  58. Adirondack Park Agency. Wetlands Effects Database and GIS for the Adirondack Park. 2004. Available online: https://apa.ny.gov/Research/ParkwideGISFinalReport.pdf (accessed on 5 January 2018).
  59. Renecker, L.A.; Schwartz, C.C. Food Habits and Feeding Behavior. In Ecology and Management of the North American Moose, 2nd ed.; Schwartz, C., Franzmann, A., McCabe, R., Eds.; Smithsonian Institution Press: Washington, DC, USA, 1997; pp. 403–440. [Google Scholar]
  60. Hinton, J.H.; Wheat, R.E.; Schuette, P.; Hurst, J.; Kramer, D.; Frair, J. Challenges and opportunities for robust population monitoring of moose along their southern range in eastern North America. J. Wildl. Manag. Accepted.
  61. Hanberry, B.B. Addressing regional relationships between white-tailed deer densities and land classes. Ecol. Evol. 2021, 11, 13570–13578. [Google Scholar] [CrossRef]
  62. Hill, R.A.; Thomson, A.G. Mapping woodland species composition and structure using airborne spectral and LiDAR data. Int. J. Remote Sens. 2005, 26, 3763–3779. [Google Scholar] [CrossRef]
  63. Blouin, J.; DeBow, J.; Rosenblatt, E.; Hines, J.; Alexander, C.; Gieder, K.; Fortin, N.; Murdoch, J.; Donovan, T. Moose habitat selection and fitness consequences during two critical winter tick life stages in Vermont, United States. Front. Ecol. Evol. 2021, 9, 642276. [Google Scholar] [CrossRef]
  64. Peek, J.M. Habitat Relationships. In Ecology and Management of the North American Moose, 2nd ed.; Schwartz, C., Franzmann, A., McCabe, R., Eds.; Smithsonian Institution Press: Washington, DC, USA, 1997; pp. 351–376. [Google Scholar]
  65. Stewart, K.M.; Bowyer, R.T.; Wiesburg, P.J. Biology and Management of White-Tailed Deer, 1st ed.; CRC Press: Boca Raton, FL, USA, 2011; pp. 181–218. [Google Scholar]
Figure 1. Study site locations of (A) Allegheny National Forest in northwestern Pennsylvania from 1993 to 2013 and (B) Adirondack Park in northern New York from 2015 to 2018.
Figure 1. Study site locations of (A) Allegheny National Forest in northwestern Pennsylvania from 1993 to 2013 and (B) Adirondack Park in northern New York from 2015 to 2018.
Forests 13 00150 g001
Figure 2. Landsat panels that were required to complete coverage of the Adirondack Park in northern New York from 2015 to 2018.
Figure 2. Landsat panels that were required to complete coverage of the Adirondack Park in northern New York from 2015 to 2018.
Forests 13 00150 g002
Figure 3. Predicted estimates (mean ± standard deviation) derived from remote sensing imagery for the total number of moose that can be supported in the Adirondack Park, New York.
Figure 3. Predicted estimates (mean ± standard deviation) derived from remote sensing imagery for the total number of moose that can be supported in the Adirondack Park, New York.
Forests 13 00150 g003
Table 1. List of raster variables derived from Landsat satellite imagery used to classify three forest treatment types (i.e., overstory removal, intermediate removal, and mature forest) in northwestern Pennsylvania from 1993 to 2013.
Table 1. List of raster variables derived from Landsat satellite imagery used to classify three forest treatment types (i.e., overstory removal, intermediate removal, and mature forest) in northwestern Pennsylvania from 1993 to 2013.
Predictors
NDVI (target scene year)
NDVI (previous scene year)
Difference NDVI (previous target)
Difference in 1.55–1.75 mm Band (previous target)
Difference in 2.09–2.35 mm Band (previous target)
Previous panel prediction (e.g., 2015 landscape for 2016 target scene)
Table 2. Reclassification of NLCD land cover into one of the seven bins for the Adirondack Park, NY, from 2015 to 2018, used to estimate associated carrying capacity.
Table 2. Reclassification of NLCD land cover into one of the seven bins for the Adirondack Park, NY, from 2015 to 2018, used to estimate associated carrying capacity.
NLCD ClassReclassified
Open waterWater
Developed, open spaceDeveloped
Developed, low intensityDeveloped
Developed, medium intensityDeveloped
Developed, high intensityDeveloped
Rock/clay/sandDeveloped
Deciduous forestMature Forest
Evergreen forestConifer
Mixed forestMature Forest
ScrublandGrass/Scrub
GrasslandGrass/Scrub
Pasture/hayAgriculture
Cultivated cropsAgriculture
Woody wetlandsWetlands
Herbaceous wetlandsGrass
Table 3. Land cover class and the estimated number of animal use days (AUD) per kilometer for moose in six different cover types [39] within the Adirondack Park, New York, from 2015 to 2018.
Table 3. Land cover class and the estimated number of animal use days (AUD) per kilometer for moose in six different cover types [39] within the Adirondack Park, New York, from 2015 to 2018.
Land Cover TypeMoose/km2
Conifer forest0
Upland deciduous forest/mixed forest0.0028
Lowland deciduous forest/mixed forest0
Wooded wetland0
Open wetland0
Regenerating forest0.0195
Table 4. The overall prediction accuracy of forest timber treatments (mature forest, intermediate removal or overstory removal) using Landsat 5 and Landsat 8 data for 1993–2013 in northwestern Pennsylvania (A) and for four different Landsat scenes for 2013–2018 in northern New York (B).
Table 4. The overall prediction accuracy of forest timber treatments (mature forest, intermediate removal or overstory removal) using Landsat 5 and Landsat 8 data for 1993–2013 in northwestern Pennsylvania (A) and for four different Landsat scenes for 2013–2018 in northern New York (B).
A. Study 1: Allegheny National Forest, PA
YearOverall Accuracy (%)Cohen’s Kappa
199395.50.92
199693.30.89
200097.40.96
200293.70.9
200694.50.92
200991.20.86
201393.10.89
B. Study 2: Adirondack Park, New York
YearOverall Accuracy (%)Cohen’s Kappa
Panel 014029
210483.90.75
201588.60.8
201689.30.81
201788.90.8
201889.20.81
Panel 014030
210487.80.78
2015840.72
201683.40.7
201782.70.69
201888.60.63
Panel 015029
201587.90.8
201688.60.82
201786.20.77
2018860.74
Panel 015030
201375.80.57
201591.20.84
201690.20.82
201790.60.84
201887.60.79
Table 5. Post remote sensing prediction land cover class landscape proportions after incorporating NLCD data for Allegheny National Forest in northwestern Pennsylvania and after incorporating NLCD data and the Adirondack Park Agency wetlands data for Adirondack Park in northern New York.
Table 5. Post remote sensing prediction land cover class landscape proportions after incorporating NLCD data for Allegheny National Forest in northwestern Pennsylvania and after incorporating NLCD data and the Adirondack Park Agency wetlands data for Adirondack Park in northern New York.
Land Cover ClassYear
Allegheny National Forest1993199620002002200620092013AVG.STDEV
Mature Forest75.5279.2471.5672.7257.0668.1563.2869.657.54
Overstory Removal1.12.01.00.813.21.01.12.94.6
Intermediate Removal7.32.73.82.86.77.812.26.23.4
Conifer7.37.36.16.16.06.56.06.40.6
Grass/Scrubland0.50.57.07.06.66.17.15.03.1
Water0.60.60.80.80.80.80.70.70.1
Developed1.01.03.23.23.23.23.22.61.1
Agriculture6.86.86.66.66.66.66.56.60.1
Adirondack Park2015201620172018AVG.STDEV
Conifer Forest16.215.515.015.915.70.5
Upland Decid/Mixed Forest32.332.132.132.031.90.5
Lowland Decid/Mixed Forest20.319.918.819.519.60.6
Wooded Wetland12.412.412.412.412.40.0
Open Wetland5.35.35.35.35.30.0
Regenerating Forest8.910.212.910.410.61.7
Ag/Developed/Grass/Scrub4.64.64.44.54.50.1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kramer, D.W.; Prebyl, T.J.; Nibbelink, N.P.; Miller, K.V.; Royo, A.A.; Frair, J.L. Managing Moose from Home: Determining Landscape Carrying Capacity for Alces alces Using Remote Sensing. Forests 2022, 13, 150. https://doi.org/10.3390/f13020150

AMA Style

Kramer DW, Prebyl TJ, Nibbelink NP, Miller KV, Royo AA, Frair JL. Managing Moose from Home: Determining Landscape Carrying Capacity for Alces alces Using Remote Sensing. Forests. 2022; 13(2):150. https://doi.org/10.3390/f13020150

Chicago/Turabian Style

Kramer, David W., Thomas J. Prebyl, Nathan P. Nibbelink, Karl V. Miller, Alejandro A. Royo, and Jacqueline L. Frair. 2022. "Managing Moose from Home: Determining Landscape Carrying Capacity for Alces alces Using Remote Sensing" Forests 13, no. 2: 150. https://doi.org/10.3390/f13020150

APA Style

Kramer, D. W., Prebyl, T. J., Nibbelink, N. P., Miller, K. V., Royo, A. A., & Frair, J. L. (2022). Managing Moose from Home: Determining Landscape Carrying Capacity for Alces alces Using Remote Sensing. Forests, 13(2), 150. https://doi.org/10.3390/f13020150

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