Assessing representation of remote sensing derived forest structure and land cover across a network of protected areas.

Protected areas (PA) are an effective means of conserving biodiversity and protecting suites of valuable ecosystem services. Currently, many nations and international governments use proportional area protected as a critical metric for assessing progress towards biodiversity conservation. However, the areal and other common metrics do not assess the effectiveness of PA networks, nor do they assess how representative PA are of the ecosystems they aim to protect. Topography, stand structure, and land cover are all key drivers of biodiversity within forest environments, and are well-suited as indicators to assess the representation of PA. Here, we examine the PA network in British Columbia, Canada, through drivers derived from freely-available data and remote sensing products across the provincial biogeoclimatic ecosystem classification system. We examine biases in the PA network by elevation, forest disturbances, and forest structural attributes, including height, cover, and biomass by comparing a random sample of protected and unprotected pixels. Results indicate that PA are commonly biased towards high-elevation and alpine land covers, and that forest structural attributes of the park network are often significantly different in protected versus unprotected areas (426 out of 496 forest structural attributes found to be different; p < 0.01). Analysis of forest structural attributes suggests that establishing additional PA could ensure representation of various forest structure regimes across British Columbia's ecosystems. We conclude that these approaches using free and open remote sensing data are highly transferable and can be accomplished using consistent datasets to assess PA representations globally.