Management practices are one of the most important factors affecting forest structure and function. Landowners in southern United States manage forests using appropriately sized areas, to meet management objectives that include economic return, sustainability, and esthetic enjoyment. Road networks spatially designate the socio-environmental elements for the forests, which represented and aggregated as forest management units. Road networks are widely used for managing forests by setting logging roads and firebreaks. We propose that common types of forest management are practiced in road-delineated units that can be determined by remote sensing satellite imagery coupled with crowd-sourced road network datasets. Satellite sensors do not always capture road-caused canopy openings, so it is difficult to delineate ecologically relevant units based only on satellite data. By integrating citizen-based road networks with the National Land Cover Database, we mapped road-delineated management units across the regional landscape and analyzed the size frequency distribution of management units. We found the road-delineated units smaller than 0.5 Ha comprised 64% of the number of units, but only 0.98% of the total forest area. We also applied a statistical similarity test (Warren’s Index) to access the equivalency of road-delineated units with forest disturbances by simulating a serious of neutral landscapes. The outputs showed that the whole southeastern U.S. has the probability of road-delineated unit of 0.44 and production forests overlapped significantly with disturbance areas with an average probability of 0.50.
Research
Mapping Regional Forest Management Units: A Road-based Framework in Southeastern Coastal Plain and Piedmont
https://doi.org/10.21203/rs.3.rs-33723/v2
This work is licensed under a CC BY 4.0 License
Journal Publication
published 25 Feb, 2021
You are reading this latest preprint version
The Southeastern United States (SEUS) forest comprises 32% of the total U.S. forestland (Oswalt et al, 2014), which combined with the productivity of the forest, places this region at the forefront of American forestry production (Fox et al, 2007). This heterogeneous landscape is composed of heavily managed forests, intensive agriculture, and multiple metropolitan areas. SEUS, although one of the most densely forested regions in the United States (Hanson, 2010), is also heavily dissected by road networks (Coffin, 2007). The diverse forest management patterns, reflecting long-term land-use legacies (Haynes et al, 2002; Josephson et al, 1989), contribute to the complex land mosaic of SEUS.
It is challenging to quantify the ecological and anthropogenic mechanisms that control the spatial structure of the forest landscape and its surrounding areas in these complex forest mosaics. Forest management is the predominant factor in forest ecology and structural patterns (Becknell et al, 2015), but little is known about how management practices are related to surrounding land-use at the regional scale. One thing that is known is that, in the SEUS, significant expansions of urban areas tend to convert forested land to urban uses and that croplands tend to transition to pine plantations (Becknell et al, 2015, Haynes, 2002; Wear and Greis, 2002 and 2013; Stanturf et al, 2003). To understand the ecological and anthropogenic influences of differently managed forests on ecosystem processes, all landscapes should be understood at multiple scales, from the local scale (forest management unit) to regional scales, with the regional scale referring to broad forest mosaics that are formed from management patches (O’Neill et al, 1996). Understanding forest management spatial patterns require defining a map-based management unit, which is the subdivision regarding effects of land use on forest ecosystems. This research seeks to further the understanding of how spatial patterns of forest management affect land use, by posing the following questions: Do roads delineate forest management units? What is the spatial distribution of road-defined forest management units in the SEUS? Moreover, how are the distributions of management units affected by different forest management approaches? And how does forest management affect nearby land use, and how does nearby land use affect forest management?
Forest management is the main driving force of forest structure in SEUS (Becknell et al, 2015) and alters forest properties and processes, which affects forest ecosystem services. (Kurz et al, 2008; Oswalt et al, 2014; Stephens et al, 2012; Kleindl et al. 2018). Forest management can be classified into four categories: production forestry, ecological forestry, passive management and preservation management (Becknell et al 2015). Production management harvests forest products and sustains the bio-productivity of the system with the sole objective of producing wood, pulp, and other forest products. In SEUS, production management based on silviculture systems, which homogenize parts of the landscape, has predominated the SEUS (Siry et al, 2002). However, forest conservation systems have evolved considerably over recent decades (Mitchener et al, 2005; Franklin et al, 2007). Ecological management uses legacies of disturbance, including intermediate stand disturbance processes such as variable density thinning and fire, and variable and appropriate recovery times to manage forests that still produce economically valuable wood products while preserving many of the values of natural forests (Franklin et al. 2007). Passive management is defined as a practice with little or no active management. We argue that all forests are managed to some degree and that doing nothing is a form of management. Preservation forestry aims to minimize the ecological footprint of society with the objectives of protecting wildlife and maintaining ecosystem services. Furthermore, certain forest management practices can suppress wildfires (Waldrop et al, 1992), prevent insect /pathogen outbreaks (Netherer et al, 2005; Faccoli et al, 2014), change water yield and hydrologic regulation (Douglass, 1983), produce wood products, provide places for hunting and recreation, and conserve habitat biodiversity.
Forestland structures, functions, and ecological processes are scale-dependent (Niemelä, 1999; Drever et al, 2006; Battaglia et al, 1998). Regionally, for the purpose of sustainable forest management, we need to develop criteria and indicators of management units. Forest management units in this study are zones or patches, which can be identified, mapped and managed according to the land-use objectives. Road networks link human activities (e.g. management practices) and surrounding physical environments (land cover). For production, preservation, and ecological forestry, in many cases, forest management units are harvest or burn units. Roads are built to create access for the managers and harvesters. For example, the preserved forest in the Ordway Swisher Biological Station (OSBS) in northcentral Florida is subdivided by road networks into management (burn) units, which is the smallest unit of land that is actively managed (OSBS User Guide, 2015). In the Joseph W. Jones Ecological Research Center in Ichauway, Georgia, and OSBS, the internal road network provides access to the research site and serves as prescribed fire breaks (Joseph W. Jones Ecological Research Center, 2010; OSBS User Guide, 2015). For passive management practices, there are currently no clear criteria defining the management units. However, in national forest systems (mostly with ecologically managed forests and multi-use production forests), existing roads and trails are used for controlling prescribed fire and wildfires (e.g. Apalachicola National Forest, Osceola National Forest, USDA, 1999). For an example of a privately-owned forest, the Red Hills region in Georgia uses roads to delineate burn units ranging approximately from five to 30 Ha (Robertson et al, 2007).
Road networks facilitate movements of humans and connect natural resources with societies and economies. As conduits for human access to nature, the physical footprint of approximately 6.6 million km of roads in the United States (U.S. Department of Transportation. 2013) has significant primary and secondary impacts on ecosystems and the distribution of species (Bennett 1991). Fifteen to twenty percent of American land is subject to the ecological effects of road networks (Forman and Alexander 1998, Forman et al. 2000). The most noticeable effects of road networks on forest structures are landscape structure changes, including reduced mean patch size, increased patch shape complexity, increased edge densities, and reduced unit connectivity. In one case, McGarigal et al. (2001) investigated the landscape structure changes of the San Juan Mountains from 1950-1993 and found that roads had a more significant ecological impact (e.g. core forest areas and patches sizes decrease) on landscape structure than logging activities.
In addition to management practices (e.g., harvesting, fertilizing), the construction of road networks divides the forested land into smaller patches, thereby increasing the potential intensity of the effects of management practices. Road networks in managed forests provide easy access for managers and harvesters to extract and regenerate resources (Demir et al, 2005). Roads may influence fire regimes by increasing fire ignition as a result of human activities (Franklin and Forman, 1987). Moreover, road networks alter the spatial configuration of management patches by functioning as firebreaks, which form new patterns in landscapes (Franklin and Forman. 1987; Nelson and Finn, 1991; Eker and Coban, 2010). By quantifying the spatial patterns of management units created by roads we may gain insight into the ecological effects of road networks on spatial forest structures within differently managed areas.
The impacts and ecological effects of roads on the landscape might be misestimated because methods measuring the road-effect zones and landscape scale effects are not yet well developed (Ries et al, 2004, Hou et al, 2013). Roads and streams may be challenging to identify, or invisible because they do not open the canopy so that many or roads are not detectable on satellite imagery or even aerial photography. The reliability of large-scale road data is also challenged due to issues of accuracy, coverage and immediacy, all of which can underestimate the extent and ecological impacts of roads on forest structures (Riitters et al, 2004). We propose that common types of forest management are practiced in road-delineated units that are detectable by remote sensing satellite images coupled with crowd-sourced road network datasets. We also describe and study the patterns of forest management units in response to land ownership and different management practices.
We focus on the Southeastern U.S. Coastal Plain and Piedmont (SEUS) region (Figure 1). The SEUS is located between Piedmont to the north and the Atlantic Ocean to the east and covers a significant portion of the southeastern United States. The SEUS is the home to the most densely production-forested region nationwide, which makes up 32% of total U.S forest cover (Oswalt et al. 2014). Based on EPA eco-region descriptions, land cover in the SEUS is a mosaic of cropland, pasture, woodland and forests (Bailey, 2004). Major silvicultural forests in SEUS are pine forests, such as slash pine (Pinus elliottii Engelm.) and loblolly pine (Pinus taeda L.) forests. European settlement and the extensive harvesting in the early 1900s removed 98% of the original longleaf pine (Pinus palustris Mill.) forests, which was one of the most dominant ecosystems in SEUS (Outcalt, 2000) and converted them to plantations of native slash pine. The SEUS forest system is a fire-dominated system with native trees adapted to short-period stand-clearing events. The primary forest management types are production and passive management due to the dominant ownership of private owners, logging companies, and investment institutions (Real Estate Investment Trusts (REITs) and Timber Investment Management Organizations (TIMOs)) (Zhang et al, 2012). This diversity of land cover types is spatially heterogenous, and patch sizes of the numerous vegetation classes vary across a wide range of scales (Figure 1).
3.1 Forest Extent
In this study, forest extent is determined by a composite of the 2006 and 2011 USGS National Land Cover Database (NLCD), which was constructed from Landsat imagery at 30-m spatial resolution (Jin et al, 2013). We aggregated 21 NLCD classes into two classes: forest (deciduous, evergreen, mixed forest and woody wetland) and non-forest (including water). Only the pixels that contain 50% or more forest area in NLCD will be considered as forested pixels. We also extracted the SEUS urban areas by using the most recent 2015 US Census Bureau’s TIGER cartographic boundary urban areas (TIGER, 2015) dataset to remove urban areas from theanalysis.
3.2 Forest Management Type
An integrated random forest classifier was built from the analysis of long-term phenological features derived from BFAST outputs and spectral entropy calculated from the Terra-MODIS enhanced vegetation index (EVI: MOD13Q data product), along with ancillary data such as land ownership, and disturbance history to classify different forest management types (Breiman, 2001; Verbesselt et al, 2010; Zaccarelli et al, 2013). The forest management type map has a spatial grain of 250 m and is a composite of phenological patterns and changes in the patterns from February 2001 through December 2016 (Figure S1). The SEUS forest management type map has an overall accuracy of 89% for a 10-fold cross-validation. The forest management raster is available for each region as georeferenced GeoTIFF rasters with a 250-meter resolution from PANGAEA (Marsik et al, 2017).
3.3 Road Networks
We selected OpenStreetMap as the primary road data and the USDA National Forest service trail and road maps (https://data.fs.usda.gov/geodata/edw/datasets.php?dsetCategory=transportation，accessed Dec 2016) as secondary in this study. OpenStreetMap is a collaborative, crowdsourced project that creates free, open, and accessible maps of road networks. OpenStreetMap is one of the most popular and well-supported Volunteered Geographic Information (VGI) datasets (Mooney et al, 2010). Community volunteers collect geographic information and submit it to the global OpenStreetMap database (Ciepluch et al, 2009). OpenStreetMap monitors road networks at near real-time and includes additional classes of roads such as private access roads and driveways in rural areas, small service roads or alleys in urban areas, and forest access roads. All those road features are critical for this research and no distinctions were drawn between the types of road, traffic volume, or other factors. OpenStreetMap shows up-to-date road networks information, which the other official road databases do not offer. The accuracy of the OpenStreetMap in our study region has been studied. For the state of Florida, Zielstra and Hochmair (2011) compared road networks dataset from different sources and concluded that OpenStreetMap was significantly better the other road databases. All OpenStreetMap data were downloaded from the website of Geofabrik:http://download.Geofabrik.de (Accessed Dec 2016). USDA National Forest Services Trails and Road Map provide the coverage of detailed transportation map in National Forests (Coghlan et al., 1997).
Road density in SEUS was measured as the total length of all roads (in kilometers) in a district divided by the total land coverage area of the district (km2) based on our developed road networks map (Figure S2).
3.4 Landscape Fire and Resource Management Planning Tools (LANDFIRE)
Disturbance data from the Landscape Fire and Resource Management Planning Tools (LANDFIRE) disturbance database were used to evaluate the management unit map. LANDFIRE is a combination of Landsat images, fire program data, and cooperator-provided field data and other ancillary databases (e.g., PAD-US)., and is a shared program between the wildland fire management programs of the U.S. Department of Agriculture Forest Services and U.S. Department of the Interior (Rollins, 2009). LANDFIRE also describes land cover/use both spatially and temporally from 1999 to 2014 and provides the existing vegetation composition map based on dominant species or group of dominant species. Spatially, LANDFIRE is a Landsat-based (30 m) database, which matches the 30-m spatial resolution of this study (https://www.landfire.gov/disturbance.php). We chose the LANDFIRE project data because the spatial scale is small enough to detect subtle changes brought about by land management practices, and large enough to reflect the characteristic variability of essential ecological processes (such as wildfire) in the appropriate spatial context. The disturbance data from LANDFIRE will be used to evaluate management units delineated by the road network in the context of intensely managed forests in SEUS (see Section 4.2, Figure S3).
3.5 Forest Ownership
Geospatial land-ownership data from federal and nongovernmental agencies were integrated for land ownership mapping (Figure S4). Forest ownership in SEUS is broadly categorized as publicly owned and privately owned according to the landowners. There are six sub-types of public ownership, which are federally protected, federal, state protected, state, military, and local. Also, there are four sub-types of private ownership: non-governmental organization, private, family, and corporate. The ownership classification implies different management objectives, as well as landowner skills, budgets and interests. Datasets from other federal and state government agencies were regrouped and classified into ten sub-types to create a comprehensive dataset that includes public land ownership and privately protected easements as well as specially designated areas and associated protection level (see Table S1). The final product is a 250-meter spatial resolution raster data depicting the forest ownership types and resampled to 30-m in this study to match the spatial resolution of the NLCD database.
4.1 Mapping Road-delineated Units
In each forest management type, the fundamental element of management practice is the management unit. In this study, we define the individual forest clusters that delineated by road-networks as “management units” and the clusters that directly derived from forest extent map as “management patches." We hereby developed two comparative methods for landscape analysis by using two sets of input data, with and without incorporating OpenStreetMap. For the method without incorporating OpenStreetMap, management patches were mapped on the forest extent map resulting from the map described in section 3.1 and the Region Group tool in ESRI ArcGIS 10.X (ESRI Inc,) to identify clusters of forest pixels that formed unique and unconnected forests.
Road-delineated forests units were mapped on the forest extent map from section 3.1 after superimposing detailed road networks with the Region Group tool to identify forest clusters as units. When superimposing road maps, all road networks were converted to one-pixel segments. After one-pixel wide road segments were derived, we converted all the forest pixels to 1 and the pixels that contained at least one road segment to non-forest pixels (30-m spatial resolution) to 0.
4.2 Geospatial Assessment
In SEUS, regional forest management activities were represented as disturbances as described by LANDFIRE data, such as clear cuts, fires, and thinning. Figure 2 shows the example view of LANDFIRE cumulative disturbance with delineated forest extent, and it can be clearly observed (with the Warren’s Index of 0.62) that the disturbed areas and delineated unit shared boundaries and show a large degree of equivalence. The spatial coincidence has been shown to facilitate the interpretation and integration of defining regional forest management units. In this study, we propose a test based on forest management unit and the geographic corresponding forest disturbances to compare the geographical similarity between the management unit and corresponding forest disturbances. The assessment of geographic image overlap is analogous to quantifying the niche similarity of two species in two dimensions. As two-dimensional rasters, both disturbance and forest extent data can be treated as homogenous and spatial-explicit datasets.
Testing the overlap between pairs of road-delineated forest extent with disturbance regions were compared using the similarity statistics of Warren’s Index (Warren et al, 2008). The values of Warren’s Index range from 0 to 1. The value of 0 means forest management units have no spatial overlap with forest disturbance areas, and 1 means all forest management units are identical to disturbance areas. The statistics of Warren’s Index assume probability distribution defined over geographic space, in which the (or ) denotes the spatial probability distribution of X (road-delineated forests compartments), or Y (probability distribution of forest disturbances) to cell . In nature, Warren’s equivalency index carries no biological assumptions concerning the parameters as being from any probability distributions. Spatially, we applied it into assessing road-delineated forest management units versus the areas of corresponding disturbances.
Firstly, the Hellinger distance was calculated to measure and compare the probability distance (Van de Vaart, 1998):
while the Warren’s statistic is:
To measure and quantify the similarity of differently managed management units with the corresponding disturbance areas, all the units were classified into four groups (ecological, passive, preservation and production management). A 10 km × 10 km grid was utilized to compute Warren’s Index over the SEUS. The inputs of this iterative Warren’s Index spatial analysis are: 1) stacked disturbance area data from LANDFIRE and 2) corresponding road-delineated forest management units. These analyses were carried out using the nivheOverlap function of dismo package in R 3.3.X (Hijmans et al., 2015).
We further tested the hypothesis of road delineated forest management units by generating neutral landscape models from two perspective as:
1) Spatially. In this study, we used the Worldwide Reference System 2 (WRS-2) row 17 path 39 (17_39 hereafter), which is one scene of Landsat Thematic Mapper (TM). As shown in Figure 3, the forests in 17_39 are heavily disturbed and fragmented, but it also shared the larges contiguous forests in SEUS (Okefenokee National Refuge). This heterogeneous landscape consists of a mixture of natural and plantation forests, urban centers, urban and rural residential areas, and commercial and small-scale agricultural operations. Instead of using LANDFIRE disturbance data, we simulated a neutral disturbance layer by mimicking the disturbance patch shape and area called “Random NLM” (Saura and Martinez-Millan, 2000; Sciaini et al, 2018). 17_39 was divided into 325 individual landscapes with the spatial resolution of 30m (10×10 km in size), which means each landscape in 17_39 contains 333×333 pixels (10,000 ha). We applied the neutral landscape generation algorithm in each landscape in 17_39 for 325 times and recalculated the Warren’s Index for each 10×10 km landscape with the simulated disturbance layer with the road networks delineated forest extent map.
2) Iteratively. We applied a range-based attribute filter to the SEUS road density map (Figure S2). The road densities range from 0 to 91.48 km/ km2, so we used the interval of 5 to randomly select a serious spatial noncorrelated 10×10 km grids with a total number of 19 landscapes in SEUS (Figure 4). For each spot, we iteratively run the Random NLM algorithm for 500 times and calculated the Warren’s Index with overlaying the simulated disturbance resulted from NLM with the road networks delineated forest extent map.
5.1 Spatial Assessment
We calculated the Warren’s Index with overlaying the road networks delineated forest extent with LANDFIRE disturbance map on a 10×10 km grid in SEUS (Figure 5). The total number of 10×10 km grids is 11,072. Figure 5 shows the probability distribution of road network delineating management unit, with an average of 0.44 and the standard deviation is 0.28. We also computed the distribution and found that most forested areas have the Warren’s Index from 0.3 to 0.5. When we extracted 17 39 (Figure 6 (a)), shown as Figure 6(b) with the mean value of 0.54 and the standard deviation of 0.17.
By replotting the spatial probability of road networks delineating management units in 17_39 (Figure 6 (c)), we recalculated the histogram of how the Warren’s Index distributed over 17_39 region after replacing the LANDFIRE data with simulated disturbance layer, as shown in Figure 6(d). The mean value of the probability distribution is 0.14 with the standard deviation of 0.082. The results show that human-derived forest disturbances can result in a non-random association with forest extent. By comparing with the LANDFIRE disturbance derived Warrans’ I map Figure 6(a), road networks make a great contribution in identifying and shaping forest patterns in the case study area 17_39, wherein SEUS reflect the management practices on the ground, and so for the road-delineated forest compartments, we hereby call them forest management units.
Road density in SEUS ranges from 0 to 91km/ km2. Table 1 lists the info of all 19 10km*10km grids with a variation of road densities, the mean values of 500 simulated landscapes. We compared the original Warren’s Index of all 19 landscapes with the mean values calculated from simulated landscapes, and found all the simulated landscape show significant difference on each landscape with the original Warren’s Index. We also plot the histogram of the landscapes to see how the simulated neutral landscape distributed (Figure 7). A linear relationship between the Warren’s Index overlaying simulated disturbances and roads with the Warren’s Index overlaying LANDFIRE disturbances with roads (Figure S5) with R2 of 0.4755. It shows the evidence that roads help to shape forest patterns by statistical test and approve the hypnosis that road delineates forest management patterns.
5.2 Forest Management Unit map
A total area of 5.24 ×107 Ha managed forest was measured in the SEUS from the forest extent map based on NLCD. When we compared the NLCD derived forest/non-forest maps between 2006 and 2011, a strong tendency of deforestation was found on our study area (1.7 × 106 Ha). The total length of the roads in SEUS is 2.26 million km, which was calculated based on the number of rasterized road network multiply the spatial resolution (30 m). The road densities in SEUS range from zero to 49.29 km/ km2 (Figure S2). We mapped the forest management units at regional scales, and the spatial distribution of those management units represented the diversity and heterogeneity of road-delineated forests clusters.
As the unit size increases, the number of units decreases in an approximately logarithmic manner (Figure 9). Overall, 177,400 units occupied the SEUS with the mean size of 29.5 Ha. We found that 35.94% of the forest management units represent more than 99% of the whole forest and ranged in area from 0.5 to 172,886 Ha (the Okefenokee National Wildlife Refuge). The remaining 64.06% of small forest compartments, which are smaller than 0.5 Ha, covers only 0.98% of the forested area.
The SEUS forests is dominated by management unit size ranging from 100 to 10,000 Ha. The forest management unit-size class map was reclassified based on the forest management unit sizes (Figure 8). From Figure 8, it can be clearly seen that riparian forests stand out as large, unbroken linear features (most of the orange colored units), which cover average unit sizes from 10,000 Ha to 100,000 Ha. A characteristic feature for the southeastern forest is that relatively small and large management units locate close to each other, surrounded with small-sized units (Figure 8 and 10). Figure 10 shows five close-ups/examples from multiple locations and landscapes. In our study area, the largest management unit is an area of the Okefenokee National Wildlife Refuge (Georgia and Florida) with the size of 172,886 Ha (Figures 8 and 10-4), followed by the Atchafalaya River basin in Louisiana at 102,367 Ha (Figure 8 and 10-3). Riparian forests are often undisturbed in SEUS without much road access because their soils are too wet for harvesting machinery and the trees are not as valuable. There is also a cluster of large forest management units (yellow: 1,000 – 10,000 Ha) along the Piedmont ecoregion, and the middle part of the southeastern plains. The Atchafalaya River basin in southern Louisiana is full of canals and river channels but no roads. There are many smaller units in Alabama and Mississippi and central Louisiana due to the relatively high road density in this region, so the landscape is broken up according to the road density map. Another example is Figure 10-5, Great Dismal Swamp National Wildlife Refuge, is the largest intact forest across southeastern Virginia and northeastern North Carolina and was established for protecting and managing the swamp’s ecosystem (USFWS, 2006)
To illustrate the contributions of road networks to our regional SEUS management unit map, a representative comparison set was done with and without incorporating road networks data. In Figure 11, we overlay the histogram of the management unit (incorporating road networks) with the histogram of the management patches (without incorporating road networks. The red histogram shows the size-based frequency distribution of patches without incorporating road networks and the purple histogram illustrates the size frequency of forest management units. When refining the management units with road networks, there are 17 times more patches (management patches) compared with the map without incorporating road networks (management units)
6.1 Management Units Under Different Management Approaches
Road density (Figure S2) has been proposed as a broad index of roads’ ecological effects in a landscape (Forman et al, 1998). The magnitude of average road densities ranges from 0.63 Km/ Km2 to 2.2 Km/ Km2 in all managed forests in SEUS. Preservation forest holds the lowest average road density of 0.63 Km/ Km2, and the passively managed forest has the highest average road density of 2.2 Km/ Km2. Ecological forest and production forest have the average road density of 1.86 Km/ Km2 and 1.29 Km/ Km2, respectively. One reason is the building of forest service roads and private logging roads, which obviously increase the road density in production forestland. Road density is a predictor of forest management intensity (Wendland et al, 2011), and the indicator of human interactions with forests (Forman et al, 2003). We compared the size-frequency distributions of management units with a map of different kinds of management (production, ecological, preservation, and passive management) derived independently. Preservation and production management had the largest patches, with means of 109.6 and 82.6 Ha, respectively. Ecological and passively managed units averaged about half as large as 73.8 and 73.0 Ha, respectively (Table 2).
We incorporated Warren’s Indexndex to assess quantitatively the geographic overlap between forest management units under different functional forest management types. For ecological management, the specific practice is designed to emulate the outcome of natural disturbance, which is to create an uneven-age stand structure to manage competition between and within multi-cohort stands. The distribution of ecological management units shows spatial heterogeneity with structurally complex stands. For passive management forest lands, as the passively managed forests mostly adopt many irregular shapes with blurred boundaries, and rupture of connectivity. For preservation management forest: mostly large government-managed land for multiple-purpose including watershed, wildlife, recreation and wilderness aspects. Accordingly, various practices may be applied to it such as harvest, cutting, retention cutting, thinning and prescribed fire. As shown in Table. 3 Warren’s Index represent the overlap between forest management units under different management types with the corresponding forest disturbance area.
The 10×10 km based spatial grid analysis of Warren’s Index is shown in Figure 5 and Figure (S6-S9). Among all four forest management types, production forestry showed the highest probability of road-delineated management units with I = 0.50, and the passive managed forests showed the lowest probability of road-delineated management units I = 0.33. As the dominant forest: production forests, the average size of plantation management units tend to be large, have a uniform composition (Figure S9), internally homogeneous and involve practice such as clearcutting and thinning. For preservation/wilderness, across the whole SEUS, the probability of road-delineated management unit is 0.44, with the standard deviation of 0.29, because the large areas of wildfires happened but with road network setting as firebreaks. Roads provide access and firebreaks, as the use of prescribed fires is widespread in SEUS and much of the prescribed fires are on private lands (Haines et al, 2001). In this study, we used a threshold of 50% of the similarity score although many useful criteria for establishing such thresholds have been proposed (Jimenez-Valverde and Lobo, 2007). In Table 3’s analysis, the criteria were just used to show binary predictions of four differently managed forests.
6.2 Ownership Representation of Forest Management Unit
Forest management activities are important links between human and environmental factors, especially at regional scale. Forest ownership patterns also explain different types of land management practices and trajectories of land cover change (Turner 1996). The aim of this part of the research is to link regional land ownership to management. We produced the SEUS ownership database (Figure S4; Table S1) by collecting the data from different sources, where the ownership was divided as Private and Public. Based on our understanding of the SEUS forest ownership, we reclassified the forest owner types into public and private with 10 sub-classes (Table.4): 1) Public: Federal protected, federal, state protected, state, military, local and NGO lands; and 2) Private: Private, family and corporate forests. Our ownership data indicate that 18.7% of the landowners are under public forest and 81.3% of private forest landowners, which covers 41.5% and 58.5% total forestland.
We can see that the type of land ownership affects units’ distribution (Table 4). The special characteristics of SEUS forest ownership patterns can result in strong contrasts in management unit distribution. The major ownership types in the region are family, corporate and state, which have different management objectives. In addition, as the road networks increased, the number of small-sized parcels also shows a substantial increase.
The mean forest management units under different ownerships range from 73.2 Ha (privately owned forests) to 115.9 Ha (state protected forests) with the standard deviation of 33.4 Ha. By comparing with the table of “Average management unit size in each forest management type”, we found that the average unit size of military land is 74.9 Ha, which is the closest to the average management unit size of ecological forestry at 73.8 Ha. The land with the ownership of state protected represents the largest average management unit at 115.9 Ha.
Quantifying forest management units under different management approaches is key step to ensure that appropriate management practices and policies are in place to maintain the array of forest ecosystem services. Regionally, forest harvesting operations are conducted by the within road-defined boundaries. Forest management practices in SEUS are specifically the activities primarily dictated by forest harvest and the needs of management for recreation and sustainability. Roads networks in SEUS are often used for setting firebreaks and timber harvest boundaries, so there should be a spatial coincidence between road-delineated management units and disturbance. In this study, we assessed the forest management from stands to regional scale, by incorporating road networks and multi-temporal disturbance remote sensing database.
A number of conclusions can be drawn from the analysis in this study:
1) Road networks play a role in delineating forests from local to regional scale. By defining the individual forest clusters delineated by road networks as “forest management units” and the clusters that directly derived from forest extent map as “management patches”, we mapped the forest extent map of both “units” and “patches” and compared them with treating “patches” map as background. There were 17 times more “units” than “patches” over the whole SEUS. And we also summarized the size distribution road delineated units, with units smaller than 0.5 Ha comprised 64% of the counts of units, these small units altogether covered only 0.98% of the total forest area.
2) We quantitatively tested the probability distribution patterns by using Warren’s Index of road-delineated management units and the corresponding forest disturbances area. The average probability of road-delineated management units is 0.44, and we also visualized the probabilities by setting a 10km×10km grid. In SEUS, the high equivalency between the road-delineated units and the corresponding areas were found at most production forests, and large-size preserved areas (e.g. Okefenokee National Wildlife Preserve and St. Marks National Wildlife Refuge).
3) The combination of remote sensing data and OpenStreetMap constitutes a useful tool to monitor, characterize and quantify land cover and management unit distributions at macrosystems scale. By using the NLCD as the forest reference data and OpenStreetMap as road networks dataset, we produced the OpenStreetMap refined management unit pattern map and analyzed the spatial size distribution of forest patterns. In addition, by incorporating the OpenStreetMap, the roads are shown to play an important role in causing fragmentation of the remnant forestlands. The size frequency distribution tells us that all of the 64% of management units are small management units (< 0.5ha) making up just 0.98% forestlands.
4) Our land ownership product indicates that 18.7% of public forest and 81.3% of private and industrial forestland owners, cover 41.5% and 58.5% area of total forestland, respectively. Management practices affected units are represented at not only the stand or local scale, but also will change the forest pattern dramatically at the regional scale. We provided substantial evidence that road networks occupy a substantial proportion of forest. From a forest management perspective, more road landscape area leads to less available land for trees at the macrosystems scale. On the other hand, logging roads and trails are an efficient way to manage forests, which we can see the relatively high road network densities in ecological and production forestry.
This study represents benefits so society in that future management decisions can be evaluated cross scales, taking account of both climate and disturbance regimes. More information on the effects of land ownership and forest management, combined with the detailed road network and a continental coverage land cover maps can aid in thwarting further forest fragmentation by promoting more reasonable road planning by land planners and decision makers.
|
GEE |
Google Earth Engine |
|
LANDFIRE |
Landscape Fire and Resource Management Planning Tools |
|
MODIS |
Moderate Resolution Imaging Spectroradiometer |
|
NLCD |
National Land Cover Database |
|
NLM |
Neutral Landscape Model |
|
OSM |
OpenStreetMap |
|
SEUS |
Southeastern U.S. Coastal Plain and Piedmont |
|
VGI |
Volunteered Geographic Information |
Acknowledgments
Acknowledgement is made of the assistance of Dr. Michael Binford and Dr. Peter Waylen, Department of Geography, University of Florida, for English writing and reviewing, for suggestions and discussion. Funding for this research was provided by the National Science Foundation Macrosystems Biology Program Grant EF #1241860.
Authors’ contributions
YD and FC conceptualized the idea for the study, FC generated the SEUS forest ownership database; YD performed data analysis and led the writing of the manuscript; FC critically reviewed the data analysis, ad contributed substantially to the writing. Both authors read and approved the final manuscript.
Funding
We acknowledge funding from the Macrosystems Biology Program Grant EF #1241860 from United States National Science Foundation (NSF).
Availability of data and materials
The dataset and code used during the current study are available at: https://doi.org/10.6084/m9.figshare.11406612.v1
Ethics approval and consent to participate
Not applicable
Consent for publication
Not applicable
Completing interests
The authors declare that they have no competing interests
- Bailey, R. G. (2004). Identifying ecoregion boundaries. Environmental management, 34(1), S14-S26.
- Battaglia M, Sands P J. Process-based forest productivity models and their application in forest management[J]. Forest Ecology and Management, 1998, 102(1): 13-32.
- Becknell J M, Desai A R, Dietze M C, et al. Assessing interactions among changing climate, management, and disturbance in forests: A macrosystems approach[J]. BioScience, 2015, 65(3): 263-274.
- Bennett A F. Roads, roadsides and wildlife conservation: a review[J]. Nature conservation 2: the role of corridors, 1991.
- Breiman, Leo. "Random forests." Machine learning 45.1 (2001): 5-32.
- Ciepluch B, Mooney P, Jacob R, et al. Using OpenStreetMap to deliver location-based environmental information in ireland[J]. SIGSPATIAL Special, 2009, 1(3): 17-22.
- CIFOR C&I Team. The CIFOR criteria and indicators generic template. No. Criteria and Indicators Toolbox Series no. 2. CIFOR, Bogor, Indonesia, 1999.
- Coffin A W. From roadkill to road ecology: a review of the ecological effects of roads[J]. Journal of transport Geography, 2007, 15(5): 396-406.
- Davis DE, Stoll MR. Southern United States: An Environmental History. Santa Barbara, Calif: ABC-CLIO; 2006.
- Demir M, Hasdemir M. Functional planning criterion of forest road network systems according to recent forestry development and suggestion in Turkey[J]. Am. J. Environ. Sci, 2005, 1(1): 22-28.
- Drever C R, Peterson G, Messier C, et al. Can forest management based on natural disturbances maintain ecological resilience? [J]. Canadian Journal of Forest Research, 2006, 36(9): 2285-2299.
- Douglass J E. The potential for water yield augmentation from forest management in the eastern United States [J]. 1983.
- Eker M, Coban H O. Impact of road network on the structure of a multifunctional forest landscape unit in southern Turkey[J]. 2010.
- FAO, 2010. Global Forest Resources Assessment 2010. FAO Forestry Paper No. 163. UN Food and Agriculture Organization, Rome
- Faccoli M, Bernardinelli I. Composition and elevation of spruce forests affect susceptibility to bark beetle attacks: Implications for forest management[J]. Forests, 2014, 5(1): 88-102.
- Forman, R. T. T., & Alexander, L. E. (1998). Roads and Their Major Ecological Effects. Annual Review of Ecology and Systematics, 29, 207–231.
- Forman, R. T. T., & Deblinger, R. D. (2000). The ecological road-effect zone of a Massachusetts (U.S.A.) suburban highway. Conservation Biology, 14(1), 36–46.
- Forman, R. T., Sperling, D., Bissonette, J. A., Clevenger, A. P., Cutshall, C. D., Dale, V. H., ... & Jones, J. (2003). Road ecology: science and solutions. Island Press.
- Fox T R, Jokela E J, Allen H L. The development of pine plantation silviculture in the southern United States[J]. Journal of Forestry, 2007, 105(7): 337-347.
- Franklin J F, Forman R T. T. Creating landscape patterns by forest cutting: ecological consequences and principles[J]. Landscape ecology, 1987, 1(1): 5-18.
- Franklin, Jerry F.; Mitchell, Robert J.; Palik, Brian J. 2007. Natural disturbance and stand development principles for ecological forestry. Gen. Tech. Rep. NRS-19. Newtown Square, PA: U.S. Department of Agriculture, Forest Service, Northern Research Station. 44 p.
- Girres J F, Touya G. Quality assessment of the French OpenStreetMap dataset[J]. Transactions in GIS, 2010, 14(4): 435-459.
- Haines, T. K., Busby, R. L., & Cleaves, D. A. (2001). Prescribed burning in the South: trends, purpose, and barriers. Southern Journal of Applied Forestry, 25(4), 149-153.
- Hanson C. Southern forests for the future[M]. 2010.
- Haynes, R.W. 2002. Forest management in the 21st century: Changing numbers, changing context. Journal of Forestry 100(2):38–43.
- Heilman G E, Strittholt J R, Slosser N C, et al. Forest Fragmentation of the Conterminous United States: Assessing Forest Intactness through Road Density and Spatial Characteristics Forest fragmentation can be measured and monitored in a powerful new way by combining remote sensing, geographic information systems, and analytical software[J]. BioScience, 2002, 52(5): 411-422
- Hijmans RJ, Phillips S, Leathwick J, Elith J, Hijmans MRJ. Package `dismo'. Circles, 9, 1. 2015.
- Hou Z, Xu Q, Nuutinen T, et al. Extraction of remote sensing-based forest management units in tropical forests[J]. Remote Sensing of Environment, 2013, 130: 1-10.
- Jin S, Yang L, Danielson P, et al. A comprehensive change detection method for updating the National Land Cover Database to circa 2011[J]. Remote Sensing of Environment, 2013, 132: 159-175.
- Joseph W. Jones Ecological Research Center. 2010. The Prescribed Fire Management Program at Ichauway 1994–2009.
- Josephson, H. (1989). A history of forestry research in the southern United States. Notes, (1462), Retrieved from http://www.srs.fs.usda.gov/pubs/3175
- Kurz, Werner A., C. C. Dymond, G. Stinson, G. J. Rampley, E. T. Neilson, A. L. Carroll, T. Ebata, and L. Safranyik. "Mountain pine beetle and forest carbon feedback to climate change." Nature 452, no. 7190 (2008): 987-990.
- Marsik, M., Staub, C., Kleindl, W., Fu, C-S., Hall, J., Yang, D., Stevens, F.R., Binford, M. PANGAEA, https://doi.org/10.1594/PANGAEA.880304 (2017).
- McGarigal, ,W. H.Romme, M.Crist, and E.Roworth. 2001. Cumulative effects of roads and logging on landscape structure in the San Juan Mountains, Colorado (USA). Landscape Ecology 16:327–349.
- Mitchener L J, Parker A J. Climate, lightning, and wildfire in the national forests of the southeastern United States: 1989-1998[J]. Physical Geography, 2005, 26(2): 147-162.
- Mooney P, Corcoran P, Winstanley A C. Towards quality metrics for OpenStreetMap[C]//Proceedings of the 18th SIGSPATIAL international conference on advances in geographic information systems. ACM, 2010: 514-517
- Coghlan, G., & Sowa, R. (1997). National forest road system and use. Draft rep. Washington, DC: US Department of Agriculture, Forest Service, Engineering staff.
- Neis P, Zielstra D. Recent developments and future trends in volunteered geographic information research: The case of OpenStreetMap[J]. Future Internet, 2014, 6(1): 76-106.
- Nelson J D, Finn S T. The influence of cut-block size and adjacency rules on harvest levels and road networks[J]. Canadian Journal of Forest Research, 1991, 21(5): 595-600.
- Netherer S, Nopp-Mayr U. Predisposition assessment systems (PAS) as supportive tools in forest management—rating of site and stand-related hazards of bark beetle infestation in the High Tatra Mountains as an example for system application and verification[J]. Forest Ecology and Management, 2005, 207(1): 99-107.
- Niemelä J. Management in relation to disturbance in the boreal forest[J]. Forest Ecology and Management, 1999, 115(2): 127-134.
- Olson D M, Dinerstein E, Wikramanayake E D, et al. Terrestrial Ecoregions of the World: A New Map of Life on Earth A new global map of terrestrial ecoregions provides an innovative tool for conserving biodiversity[J]. BioScience, 2001, 51(11): 933-938.
- O'Neill R V, Hunsaker C T, Timmins S P, et al. Scale problems in reporting landscape pattern at the regional scale[J]. Landscape ecology, 1996, 11(3): 169-180.
- OSBS User Guide V 09/2015, http://www.ordway-swisher.ufl.edu/docs/OSBS_User_Guide.pdf
- Oswalt S N, Smith W B, Miles P D, et al. Forest Resources of the United States, 2012: a technical document supporting the Forest Service 2010 update of the RPA Assessment[J]. 2014.
- Outcalt, K.W. 2000. The longleaf pine ecosystem of the South. Native Plants Journal, 1(1), pp.42-53.
- PAD-US, Protected areas database of the United States; Available from https://gapanalysis.usgs.gov/padus/data/download/
- Ries L, Fletcher Jr R J, Battin J, et al. Ecological responses to habitat edges: mechanisms, models, and variability explained[J]. Annu. Rev. Ecol. Evol. Syst., 2004, 35: 491-522.
- Riitters KH, Wickham JD, Coulston JW. 2004. Use of road maps in national assessments of forest fragmentation in the United States. Ecology and Society 9(2):13.
- Robertson, K. M., & Ostertag, T. E. (2007, November). Effects of land use on fuel characteristics and fire behavior in pinelands of southwest Georgia. In Proceedings of the 23rd Tall Timbers Fire Ecology Conference: Fire in Grassland & Shrubland Ecosystems, San Diego, California (Vol. 27, pp. 181-191).
- Rollins M G. LANDFIRE: a nationally consistent vegetation, wildland fire, and fuel assessment[J]. International Journal of Wildland Fire, 2009, 18(3): 235-249.
- Saura, S., & Martinez-Millán, J. (2000). Landscape patterns simulation with a modified random clusters method. Landscape ecology, 15(7), 661-678.
- Sciaini, M., Fritsch, M., Scherer, C., & Simpkins, C. E. (2018). NLMR and landscapetools: An integrated environment for simulating and modifying neutral landscape models in R. bioRxiv, 307306.
- Siry J P. Intensive timber management practices[J]. Southern forest resource assessment, 2002, 14: 327-340.
- Stephens S L, McIver J D, Boerner R E J, et al. The effects of forest fuel-reduction treatments in the United States[J]. BioScience, 2012, 62(6): 549-560.
- Technical report: 2013 Forest Plan Monitoring and Evaluation Report National Forests in Florida. USDA.
- TIGER Cartographic Boundary – Urban Areas, prepared by the U.S. Census Bureau, 2015; Available from https://www.census.gov/geo/maps-data/data/kml/kml_ua.html
- Turner M G, Wear D N, Flamm R O. Land Ownership and Land‐Cover Change in the Southern Appalachian Highlands and the Olympic Peninsula[J]. Ecological applications, 1996, 6(4): 1150-1172.
- United States Department of Agriculture Forest Service. (1999). Final environmental impact statement for the land and resource management plans for National Forests in Florida.
- US Fish and Wildlife Service. (2006). Great Dismal Swamp National Wildlife Refuge and Nansemond National Wildlife Refuge Final Comprehensive Conservation Plan. Available on-line: http://www. fws. gov/uploadedFiles/Region_5/NWRS/South_Zone/Great_Dismal_Swamp_Complex/Great_Dismal_Swamp/FinalCCP_GDS. pdf.
- Van der Vaart, A. W. (1998). Asymptotic statistics (Vol. 3). Cambridge university press.
- Waldrop T A, White D L, Jones S M. Fire regimes for pine-grassland communities in the southeastern United States[J]. Forest Ecology and Management, 1992, 47(1-4): 195-210.
- Warren, D. L., Glor, R. E., & Turelli, M. (2008). Environmental niche equivalency versus conservatism: quantitative approaches to niche evolution. Evolution, 62(11), 2868-2883.
- Wear D N, Greis, J. (2002). Southern forest resource assessment: summary of findings. Journal of Forestry.
- Wear D N, Greis J G. The southern forest futures project: summary report[J]. 2012.
- Wear D N, Greis J G. The southern forest futures project: Technical report[J]. 2013.
- Wendland, K. J., Lewis, D. J., Alix-Garcia, J., Ozdogan, M., Baumann, M., & Radeloff, V. C. (2011). Regional-and district-level drivers of timber harvesting in European Russia after the collapse of the Soviet Union. Global Environmental Change, 21(4), 1290-1300.
- Zaccarelli, N., Li, B. L., Petrosillo, I., & Zurlini, G. (2013). Order and disorder in ecological time-series: introducing normalized spectral entropy. Ecological Indicators, 28, 22-30.
- Zielstra D and Hochmair H H 2011 Digital street data: Free versus proprietary. GIM International 25(7): 29–33
- Zhang D, Butler B J, Nagubadi R V. Institutional timberland ownership in the US South: magnitude, location, dynamics, and management[J]. Journal of Forestry, 2012, 110(7): 355-361.
Table 1. Spatial assessment covering the different road densities.
|
ID |
Road Density |
Mean of Simulated Landscapes |
St of Simulated Landscapes |
I |
|
1 |
0 |
0.16 |
0.008 |
0.41* |
|
2 |
30.23 |
0.06 |
0.009 |
0.10* |
|
3 |
35.15 |
0.12 |
0.008 |
0.66* |
|
4 |
15.03 |
0.10 |
0.009 |
0.25* |
|
5 |
40.63982391 |
0.02 |
0.008 |
0.20* |
|
6 |
4.9943223 |
0.26 |
0.006 |
0.52* |
|
7 |
85.59053802 |
0.19 |
0.007 |
0.49* |
|
8 |
49.92194366 |
0.23 |
0.007 |
0.55* |
|
9 |
24.98642921 |
0.25 |
0.006 |
0.47* |
|
10 |
9.908494949 |
0.13 |
0.009 |
0.30* |
|
11 |
65.24546814 |
0.10 |
0.008 |
0.38* |
|
12 |
75.59457397 |
0.10 |
0.008 |
0.22* |
|
13 |
19.80970001 |
0.24 |
0.007 |
0.66* |
|
14 |
59.10257721 |
0.02 |
0.007 |
0.31* |
|
15 |
45.15394592 |
0.11 |
0.008 |
0.33* |
|
16 |
55.00973892 |
0.12 |
0.008 |
0.23* |
|
17 |
68.21031952 |
0.17 |
0.008 |
0.34* |
|
18 |
79.92221832 |
0.03 |
0.008 |
0.34* |
|
19 |
91.88208771 |
0.11 |
0.008 |
0.24* |
Note: * means significant different with a 95% confidence interval
Table. 2. Average management unit size in each forest management type
|
|
Average Management Unit Size (Ha) |
Standard Error （Size） |
Average Road Density (Km/ Km2) |
Standard Error （Density） |
|
Ecological Management |
73.8 |
0.124 |
1.86 |
0.0017 |
|
Passive Management |
73.0 |
0.037 |
2.20 |
0.0005 |
|
Preservation Management |
109.6 |
0.079 |
0.63 |
0.0012 |
|
Production Management |
82.6 |
0.015 |
1.29 |
0.0052 |
Table 3. Warren’s Index between different managed units and corresponding disturbances overlap with stacked LANDFIRE disturbance
|
Management Type |
Warrens I |
Standard Dev |
|
Ecological |
0.38 |
0.26 |
|
Passive |
0.38 |
0.28 |
|
Preservation |
0.44 |
0.29 |
|
Production |
0.50 |
0.24 |
Table. 4 Forest mean unit sizes based on different ownership types in SEUS.
|
Owner Class |
Owner Type |
Mean Unit Area (Ha) |
Stdev |
Area (Ha) |
|
Public |
Federal Protected |
95.4 |
54.0 |
1,291,925 |
|
Public |
Federal |
89.7 |
32.9 |
2,750,912 |
|
Public |
State Protected |
112.6 |
33.4 |
2,793,425 |
|
Public |
State |
96.8 |
52.1 |
1,277,169 |
|
Public |
Military |
73.4 |
40.7 |
740,831 |
|
Public |
Local |
89.4 |
45.5 |
1,635,244 |
|
Public |
NGO |
93.6 |
42.9 |
130,331 |
|
Private |
Private |
73.6 |
42.9 |
11,695,831 |
|
Private |
Family |
72.9 |
36.3 |
23,247,681 |
|
Private |
Corporate |
86.9 |
35.8 |
13,628,388 |