Next Article in Journal
Image Formation of Azimuth Periodically Gapped SAR Raw Data with Complex Deconvolution
Next Article in Special Issue
Mapping Spatial Variations of Structure and Function Parameters for Forest Condition Assessment of the Changbai Mountain National Nature Reserve
Previous Article in Journal
The Radiative Transfer Characteristics of the O2 Infrared Atmospheric Band in Limb-Viewing Geometry
Previous Article in Special Issue
A Multi-Temporal Object-Based Image Analysis to Detect Long-Lived Shrub Cover Changes in Drylands
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessment of Night-Time Lighting for Global Terrestrial Protected and Wilderness Areas

Key Laboratory of Geographical Processes and Ecological Security in Changbai Mountains, Ministry of Education, School of Geographical Sciences, Northeast Normal University, Changchun 130024, China
*
Authors to whom correspondence should be addressed.
Remote Sens. 2019, 11(22), 2699; https://doi.org/10.3390/rs11222699
Submission received: 1 October 2019 / Revised: 10 November 2019 / Accepted: 15 November 2019 / Published: 18 November 2019
(This article belongs to the Special Issue Remote Sensing Applications in Monitoring of Protected Areas)

Abstract

:
Protected areas (PAs) play an important role in biodiversity conservation and ecosystem integrity. However, human development has threatened and affected the function and effectiveness of PAs. The Defense Meteorological Satellite Program/Operational Linescan System (DMSP/OLS) night-time stable light (NTL) data have proven to be an effective indicator of the intensity and change of human-induced urban development over a long time span and at a larger spatial scale. We used the NTL data from 1992 to 2013 to characterize the human-induced urban development and studied the spatial and temporal variation of the NTL of global terrestrial PAs. We selected seven types of PAs defined by the International Union for Conversation of Nature (IUCN), including strict nature reserve (Ia), wilderness area (Ib), national park (II), natural monument or feature (III), habitat/species management area (IV), protected landscape/seascape (V), and protected area with sustainable use of natural resources (VI). We evaluated the NTL digital number (DN) in PAs and their surrounding buffer zones, i.e., 0–1 km, 1–5 km, 5–10 km, 10–25 km, 25–50 km, and 50–100 km. The results revealed the level, growth rate, trend, and distribution pattern of NTL in PAs. Within PAs, areas of types V and Ib had the highest and lowest NTL levels, respectively. In the surrounding 1–100 km buffer zones, type V PAs also had the highest NTL level, but type VI PAs had the lowest NTL level. The NTL level in the areas surrounding PAs was higher than that within PAs. Types Ia and III PAs showed the highest and lowest NTL growth rate from 1992 to 2013, respectively, both inside and outside of PAs. The NTL distributions surrounding the Ib and VI PAs were different from other types. The areas close to Ib and VI boundaries, i.e., in the 0–25 km buffer zones, showed lower NTL levels, for which the highest NTL level was observed within the 25–100 km buffer zone. However, other types of PAs showed the opposite NTL patterns. The NTL level was lower in the distant buffer zones, and the lowest night light was within the 1–25 km buffer zones. Globally, 6.9% of PAs are being affected by NTL. Conditions of wilderness areas, e.g., high latitude regions, Tibetan Plateau, Amazon, and Caribbean, are the least affected by NTL. The PAs in Europe, Asia, and North America are more affected by NTL than South America, Africa, and Oceania.

Graphical Abstract

1. Introduction

A protected area (PA) is defined as a geographical space that is recognized, dedicated, and managed, through legal or other effective means, to achieve the long-term conservation of nature with associated ecosystem services and cultural values in mind [1]. In general, protected areas include national parks, national forests, national seashores, all levels of natural reserves, wildlife refuges and sanctuaries, and designated areas for the conservation of native biological diversity and natural and cultural heritage and significance. Protected areas also include some of the last frontiers that have unique landscape characteristics and ecosystem functions in wilderness conditions [2].
PAs reflect the efforts to protect the world’s threatened species and their habitats. PAs are increasingly recognized as essential providers of ecosystem services and biological resources, key components in climate change mitigation strategies, and in some cases, are vehicles for protecting threatened human communities or sites of great cultural and spiritual value.
Humans have created protected areas over the past millennia for a multitude of reasons [3,4]. The establishment of Yellowstone National Park in 1872 by the US Congress ushered in the modern era of governmental protection of natural areas, which catalyzed a global movement [5,6]. However, even with the implementation of a tremendous variety of monitoring programs, as well as conservation planning efforts and achievements, species’ population decline, biodiversity loss, extinction, system degradation, pathogen spread, and state change events are occurring at unprecedented rates [7,8]. The effects are augmented by continued changes in land use, urbanization, and invasive spread alongside the direct, indirect, and interactive effects of climate change and disruption [3,4]. Protected areas become more important in serving as indicators of ecosystem conditions and functions, either by their status and/or by comparison with unprotected adjacent areas. Protected areas are highly prized by society with diversified representative characteristics. Earth’s remaining wilderness areas are becoming increasingly important buffers against the changing environment. However, they are not yet an explicit target in international policy frameworks [9].
Over the past four decades, more PAs have been established around the world. PAs play a vital role in conserving biodiversity; specifically, PAs provide a paradise for endangered and species from declining habitats and a sanctuary for over-harvested and poached species [10,11,12]. PAs have a positive impact on the local environment, such as maintaining water resources, regulating climate, and preventing forest damage [13], among other benefits [14,15].
On the other hand, ecosystems of PAs have been disturbed by human development, in particular, urbanization, such that biodiversity and habitats in those PAs have been reduced [16,17,18,19,20,21]. The effectiveness of PAs has been affected by the land use of the surrounding areas. Nightlight may impose problems to PAs [22,23]. The threat from human development to PAs is a concern. Currently, assessment studies on the impact of human development on PAs are mainly based on data, such as human population [14] and housing [24] at regional or local scales. However, accurate, reliable, and comprehensive population or housing statistics are often not available in global scale monitoring.
The Defense Meteorological Satellite Program/Operational Linescan System (DMSP/OLS) night-time stable light (NTL) dataset is a timely and effective data source for monitoring human development and mapping the dynamics of land cover on a regional or global scale [25,26,27,28]. The NTL data have a wide time span and are suitable for dynamic analysis on long-term sequences. The NTL can be measured and used as a proxy of human development. NTL imposes impacts on several taxa in terrestrial and aquatic ecosystems, including mammals, birds, reptiles, amphibians, fish, invertebrates, and plants. It is necessary to study the lighting conditions of reserves [29,30,31]. DMSP/OLS NTL is an excellent data source for analyzing the impact of human development on various ecosystems [32], vegetation [33], and biodiversity [34] over long-term sequences and over large geographical areas and in remote locations [35].
DMSP/OLS NTL data have been used as an effective indicator to evaluate the conservation of protected areas [36,37,38]. Reported studies have investigated the difference in NTL between the interior and the surrounding areas and have revealed that the NTL level of the PAs is lower than that of the surrounding areas [29]. Studies suggested that the NTL level of the boundary of PAs is particularly high [39]. In addition, some studies have combined NTL with other data to conduct research, for example, Geldmann et al. used the NTL and population data to construct a temporal human pressure index (THPI) on the time series, which quantifies the human pressure on the protected area. However, the spatial resolution of the study was 10 km2, and only attempted comparisons between 1995 and 2010 [40].
In this study, we used NTL data as an indicator of the intensity of human development to portray the characteristics of NTL in various types of PAs defined by the International Union for Conversation of Nature (IUCN). The objectives of this study were to: (1) reveal human development surrounding different types of terrestrial PAs, growth rate, and trend characteristics measured using the NTL level; and (2) make an assessment about the NTL effects on global terrestrial PAs.

2. Materials and Methods

2.1. Protected Areas and Data

We used the 2019 World Database on Protected Areas (WDPA, https://protectedplanet.net/), as the most updated data source for obtaining the coverage of global protected areas. We selected seven types of terrestrial PAs defined by the IUCN, including Ia, Ib, and II–VI, as defined in Table 1. In order to explore the spatial distribution and change of NTL within PAs and their vicinity, we examined buffer zones with distances of 0–1 km, 1–5 km, 5–10 km, 10–25 km, 25–50 km, and 50–100 km from PA boundaries. This range of buffer distances was chosen to capture different types of NTL effects. Human activities within PA boundaries and within a 1 km distance exert a direct and significant influence on protected areas, e.g., habitat loss, noise, and light pollution [41]. At further distances, artificial surfaces contribute to the landscape disturbances and effects, such as the isolation of protected areas, disruption of connectivity, and introduction and spread of invasive species. Even as far as 50 km from the protected area, human development can impose effects on PAs. PAs may be impacted even if the source of lighting lies kilometers away owing to skyglow [24].

2.2. DMSP/OLS NTL Data

The Operational Linescan System (OLS) is one of the sensors that is carried by the Defense Meteorological Satellite Program (DMSP). The DMSP/OLS NTL data came from the National Geophysical Data Center (NGDC) under the National Oceanic and Atmospheric Administration (NOAA), which eliminates accidental noise effects, such as clouds and flaring, with a 0.008333 degrees spatial resolution. The number of each pixel is not the actual light level on the ground, but rather the relative brightness level recorded as a digital number (DN) from 0 to 63. Currently, NTL data from 1992 to 2013 are available through online access. DMSP/OLS is different from Landsat, Satellite Pour l’Observation de la Terre (SPOT), and the Advanced Very High Resolution Radiometer (AVHRR) sensors that use the ground object to monitor the reflected radiation characteristics of sunlight. The DMSP/OLS sensor can work at night and capture the low intensity of urban lights and even small-scale residential areas and traffic lights, providing powerful data support for monitoring human development research on a large scale [42].

2.3. Data Processing

For origin data, discrepancies appeared between the DN values and the number of lit pixels from different satellites in the same year, and abnormal fluctuations appeared in the DN values for different years derived from the same satellite. It was necessary to calibrate the original NTL between different years and satellites. In this study, we assumed that all regions of the world that developed positive reflection in NTL would keep the DN value and the number of lit pixels either consequently increased or remained the same. The NTL were corrected using four main steps as described below.

2.3.1. Inter-Annual Correction

As the economy and population continue to develop over time, the assumption was that all PA regions would be positively developing and increasing in terms of the amount of NTL. Therefore, the DN values of each pixel in the time series would either increase or remain unchanged. With this assumption, the light pixels detected in the early image would remain in the latter image, and the DN value should not be larger than the DN value detected in their subsequent images. If the pixel that detected the light in the earlier image disappeared in the later year, the pixel with light would be considered as an unstable light pixel, and its DN value would be replaced with a value of 0. Second, the DN value of each stable light pixel was corrected to ensure that the DN value in the early image was smaller than the DN value of the corresponding position pixel in the later image. Due to the replacement of satellites between 1992 and 2013, the light sensitivity of each sensor was different, leading to the difference in the number of lit pixels and the DN values of the corresponding position pixels from different sensors. Thus, each satellite was interannually corrected using the following step [43]:
D N ( n , i ) = { 0 ,                                     D N ( n + 1   , i ) = 0   D N ( n 1   , i ) ,     D N ( n + 1   , i ) > 0       n d       D N ( n 1   , i ) > D   D N ( n , i ) D N ( n , i ) ,               o t h e r w i s e             ( n   =   1992 ,   1993 ,   ,   2013 ) ,
where DN(n−1, i), DN(n, i), and DN(n+1, i) are the DN values of the ith lit pixel of the NTL data in the (n − 1)th, the nth, and (n + 1)th years, respectively.

2.3.2. Inter-Satellite Correction

The NTL from 1992 to 2013 were derived from multiple sensors. The data were not continuous and the data between different sensors could not be directly compared. It was necessary to correct the sensors to improve the continuity and comparability of the NTL. Although the image of satellite F18 did not intersect with those of the other satellites, there was an intersection between F10 and F12, F12 and F14, F14 and F5, and F15 and F16. First, based on F10, using the data of the intersection year, 100,000 pixels were randomly selected on each continent to establish a minimum linear-square regression model between two satellites; thus, the data of F12 could be corrected, which could then be used to correct the data of F14. The same process was applied to data from other satellite pairs.
The regression formulas are as follows:
Correct F12 based on F10:
F 12 = 0.870 × F 10 0.078   ( R 2   =   0.925 )
Correct F14 based on F12:
F 14 = 0.927 × F 12 1.709   ( R 2   =   0.965 )
Correct F15 based on F14:
F 15 = 0.961 × F 14 + 1.408   ( R 2   =   0.951 )
Correct F16 based on F15:
F 16 = 1.062 × F 15 + 1.672   ( R 2   =   0.960 )

2.3.3. Intra-Annual Correction

The intra-annual composition correction used images from two satellites in the same year to remove any intra-annually unstable lit pixels. There were two images from two different satellites in 1994 and 1997–2008. We used the average DN values of the two NTL images to calculate the annual NTL data for these years. First, we examined all the lit pixels to determine whether the pixel was an unstable light pixel during the year. If only one satellite was detected, the light pixel was defined as an unstable light pixel during the year. Second, in intra-annual composites, the DN values of intra-annually unstable lit pixels were replaced with values of 0, and the DN value of each intra-annually stable lit pixel was replaced by the average DN value of two NTL images from the same year. This produced one intra-annual composite for each year [43].
D N ( n , i ) = { 0 ,                                                       D N ( n , i ) a = 0 | D N ( n , i ) b = 0 D N ( n , i ) a + D N ( n , i ) b ,                     otherwise                     ( n   =   1992 ,   1993 ,   ,   2013 ) ,
where D N ( n , i ) a and D N ( n , i ) b are the DN values of the ith lit pixel from two NTL data in the nth year, and DN(n,i) is the DN value of the ith lit pixel of the intra-annual composite in the nth year.

2.3.4. Correction for Different Sensors

We performed another first-step correction on this series of data after the first three steps, which resulted in a dataset that grew positively and with continuity.

2.3.5. Evaluation of the Calibrated NTL Time Series

A number of studies have shown that the DMSP/OLS NTL data are correlated with economic activities [44,45,46,47,48] and have used gross domestic product (GDP) and electricity consumption (EC) to evaluate the performance of intercalibration techniques [28,49,50,51,52,53]. We applied this assessment method in this study to evaluate the calibrated NTL time series. We obtained country-level GDP data (1992–2013) from the World Bank (http://data.worldbank.org/) and EC records derived from International Energy Statistics (http://www.eia.gov/beta/international/data/browser/#). We compared the raw NTL and the calibrated NTL with GDP and EC data (Figure 1) to reduce the inconsistency of the data and enhanced correlation with GDP and EC. We calculated the Pearson correlation between the raw NTL with GDP and EC, and between the calibrated NTL and the two ancillary datasets (Figure 2). The calibrated method improved the correlations between the NTL time series and GDP and EC.

2.4. Calculation of the Trend Using Sen’s Slope

The Sen’s slope method can be used to avoid the loss of time series data and the influence of data distribution on the analysis results to eliminate the interference of outliers in a time series [54]. Sen’s slope has been used in the analysis of long-term sequence data sets to detect the magnitude of the change trend. We used Sen’s slope to calculate the trend of the NTL from 1992 to 2013. The calculation formula is:
Q = m e d i a n X j X i j i                                       0 < i < j n
where Q is the Sen’s slope; and Xj and Xi are the average DN value corresponding to j and i year, respectively. If the length of the time series is n, the number of Qi is N = n(n − 1)/2; and the Q is determined by N.
                      Q = { Q [ ( N 1 ) / 2 ]                                                                                       N   i s   o d d ( Q ( N / 2 ) + Q [ ( N + 2 ) / 2 ] ) / 2                         N   i s   e v e n
where Q > 0, Q < 0 and Q = 0 indicate that there is a rising trend, a downward trend, and no obvious trend, respectively.

3. Results

3.1. Spatial Distribution and Trend of NTL for Global PAs

We superimposed the global PAs and the global NTL, except for Antarctica (Figure 3). The result revealed that the DN values for most of the PAs were lower than that of the surrounding areas. Further assessment ranked the average NTL DN values of all PAs on each continent as: Europe > Asia > North America > South America > Africa > Oceania. The DN values in Europe was much higher than those in the other continents (Figure 4). The ranking of the average change trend on all continents was consistent with the average DN value, i.e., the highest was in Europe and the lowest was in Oceania (Figure 5).

3.2. Average NTL Level in Different PAs and Buffers

Figure 6 illustrates that NTL DN values in types III and V PAs were significantly higher than that of the other five types of PAs in the interior and buffer zones. With the buffer radius increased, the average NTL DN values of all PAs increased first and then decreased (Figure 7). The highest DN values of different types of PAs appeared in the range of 1–10 km.
From the time series, the NTL level within and outside PAs increased without exception, but the NTL distributions among seven types of PAs had different characteristics (Figure 8).
Type Ia PAs represent strict nature reserves. As shown in Figure 8a, the average NTL DN value in buffer zones of type Ia PAs increased in the 0–10 km buffer zone and then decreased with over the 10–100 km buffer zone, reaching the maximum DN value in the 5–10 km buffer zone. The average DN value in the 0–100 km buffer zone of type Ia was much higher than the average DN value within the boundaries of the PAs.
Type Ib PAs represent wilderness area with large unmodified or slightly modified areas (Figure 8b). The average DN value within boundaries of type Ib was lower than that of type Ia PAs with the lowest NTL value among all types of PAs. The average DN value of the type Ib buffer increased from the 0–1 km to the 50–100 km buffer zones.
Type II PAs represent national parks (Figure 8c). The change trend of the DN value of the type II PAs buffer was similar to that of Ia (Figure 8a), suggesting a trend of increasing in the close buffer zones (0–10 km) and decreasing in distant buffer zones (10–100 km), with the maximum DN value in the 5–10 km buffer areas. The fluctuation of the NTL DN value of type II was the lowest among all types of PAs.
Type III PAs represent natural monuments or features (Figure 8d). The difference between DN values within PAs and 0–100 km buffer zones were the highest among all types of PAs. The NTL outside of type III PAs increased in the 1–5 km buffer zone and then decreased, indicating more human development in this range.
Type IV PAs represent habitat/species management areas (Figure 8e). The nearest (0–1 km) buffer zone had the lowest DN value. The area with the highest NTL level was concentrated in the 1–10 km buffer zone.
Type V PAs represent protected landscapes/seascapes (Figure 8f). The DN value of type V PAs was much higher than that of the other types of PAs (Figure 9). The average DN value and the range of the buffer zones was also the highest among all PA types (Figure 9). The brighter areas around the PAs were concentrated within 1–25 km, and the 0–1 km and 25–100 km areas were darker.
The average DN values in type VI PAs buffers were the lowest among all PAs (Figure 8g). The DN values of the 0–1 km buffer zone were the lowest among all types of PAs. As distance from the boundaries of the PAs increased, the DN value in the buffer zones increased, and reached a maximum in the 25–50 km buffer zone.
The results showed that the NTL levels within PAs were lower than that of surrounding areas, indicating that human development was limited and controlled by the boundary of the PAs. The NTL level within and outside the PAs increased from 1992 to 2013 (Figure 8). The areas with the highest NTL level in the buffer zone constantly approached the boundary of the type Ia PAs (Figure 8a).
In general, the ranking of the NTL level among the seven types of PAs was V > IV > III > VI > II > Ia > Ib (Table 2). The ranking of NTL level in the 0–100 km buffer zones outside the boundaries of types of PAs was V > III > IV > Ia > II > Ib > VI (Figure 7). The ranking of the NTL level of buffer zones and the interior was 5–10 km > 1–5 km > 10–25 km > 25–50 km > 50–100 km > 0–1 km. The NTL level within the PAs was significantly lower than that observed outside the PAs (Table 2). For most types of PAs, e.g., Ia, II, III, IV, and V, the brightest areas around the PAs were concentrated in the 1–25 km buffer zone. However, for types Ib and VI, the brightest areas around PAs were concentrated at further distances (Figure 9).

3.3. NTL Growth Rate in Different PAs and Buffers

The 1992–2013 NTL growth rates of every type of PAs and their buffer zones were doubled, except within the 0–1 km for type III PAs (Figure 10), indicating enhanced human development within and around PAs. The NTL growth rate for type Ia PAs was the highest both within and the surrounding outside among all types. Especially, NTL growth rates within the PAs and in the 0–10 km buffer zones were much higher than that of other types. The DN value of the interior Ia increased from 0.04 in 1992 to 0.19 in 2013, representing an increase of 378%. The 0–1 km buffer zone increased by 441% from 0.54 in 1992 to 2.93 in 2013. The growth rates of the 1–5 km and 5–10 km buffer zones were 372% and 274%, respectively. The NTL growth rate from 1992 to 2013 for type Ib PAs decreased from 217% for the internal area to 135% for the furthest buffer zone. The growth rate in the 0–1 km buffer zone was 215%, which was much greater than that of other buffer zones. NTL growth rates for type Ib PAs were much lower than that of Ia. The NTL growth rate for type II PAs ranged from 160% in the 5–10 km buffer zone to 226% in the outermost 50–100 km buffer zone. The growth rate within the type II PAs was 222%, ranking second. For type III PAs, the NTL growth rate was 114% within the protected boundary area. The average growth rate of the buffer zones was the lowest among all types of PAs. The growth rate increased with distance, from 90% in the 0–1 km buffer zone to 170% in the 50–100 km buffer zone. The NTL growth rate within type IV PAs (178%) was slightly higher than that of its surrounding areas (157–175%). The difference in growth rates between buffer zones was small, with a minimum of 157% in the 1–5 km buffer zone and a maximum of 175% in the 5–25 km buffer zone. The NTL growth rate within the type V PAs was 157%. The growth rate of the 0–1 km buffer zone was 164%, which was the highest among all type V PA buffers. The growth rate of the remaining buffer zones decreased with distance from the boundary, i.e., from 133% in the 1–5 km buffer zone to 151% in the 50–100 km buffer zone. From 1992 to 2013, the NTL growth rate of type VI PAs was 245%. The buffers near the boundary had a higher growth rate than that of the buffers farther away. The highest growth rate (263%) occurred in the 0–1 km buffer zone, and the minimum growth rate was 142% in the 50–100 km buffer zone.

3.4. Trends in Different PAs and Buffer Zones

NTL DN values of all PAs and their buffers showed a significant increasing trend from 1992 to 2013 (Figure 11). Except for the type V PAs, the NTL trends within the other six types of PAs were lower than that in any other buffer zone at the 0–100 km. The change trend within type Ib PAs was the lowest of all the PAs, with a value close to 0. The change trend of type V (0.13 DN/year) was the highest, which was greater than the trend in 50–100 km buffer zone. The average trend of type V PAs buffers was the lowest among all types of PAs.
The trend for type Ia PAs was significantly higher than that of type Ib, indicating that Ia was affected more by human development. The change trends of different types of PAs showed characteristics. The rank of mean values of change trends in buffer zones were V > III > Ia > IV > II > Ib > VI. The change trend of the type V PAs was the highest, e.g., 0.26/year in 1–5 km and 0.24/year in 5–10 km buffers.
The trends in type V, III, Ia, and IV PAs were concentrated in the 1–10 km buffer zone, indicating active urban development. The type VI, Ib, and II PAs had the lower change trends inside PAs, but higher change trends in the peripheral areas far from the PA boundary. The change trends in the area near the PA boundary between 0–10 km were relatively low.

4. Discussion

4.1. The NTL Distribution Pattern in Different Types of PAs

The primary objective of type Ia PAs is to conserve regionally, nationally, and globally outstanding ecosystems, species (occurrences or aggregations), and/or geodiversity features: these attributes will have been formed mostly or entirely by non-human forces [55]. However, the global NTL level inside type Ia PAs was not the lowest. The NTL value in type Ia PAs was 2.5 times higher than that of the type Ib PAs. Among seven types of PAs, the NTL growth rate inside type Ia was very high from 1992 to 2013, reaching 378%, which was far greater than the others. The NTL for 0–1 km and 1–5 km buffer zones had growth rates that were ranked highly, with values of 441% and 372%, respectively. The area with the highest NTL level around Ia consistently approached the boundary of the PAs. The brightest NTL areas between 1992 and 1999 was in the 25–50 km buffer zone, which shifted to the 5–10 km buffer zone between 2001 and 2003. According to the growth rate and trend of each buffer zone, the brightest NTL area would be encroaching in the 0–5 km area quickly. Human development in the 0–1 km outside PAs could have direct impacts [56,57]. Because the baselines of the nightlight level of different types of protected areas are very different, the actual change values corresponding to the same growth rate are very different even though the growth rates are the same. For example, the NTL growth rate of type Ia PAs exceeded the others, but the absolute increase of the NTL DN value inside and around the Ia PAs was still low in comparison to other types because its baseline of the nightlight level was far less than other types. It is noteworthy that, according to the guidelines, Ia will be degraded or destroyed when subjected to all but very limited human impact [55]. In addition, the area of category Ia was generally small, and human-induced impacts contributed more in small PAs than to stochastic processes [58].
For type Ib PAs, the average NTL DN value was 0.04, which was the lowest among all types. Also, the NTL levels in and out of the PAs boundaries was lower than type Ia PAs. The 50–100 km buffer zone was the area with the most concentrated human population around Ib PAs. The NTL distribution for type II PAs were similar to type VI and both were lower among the types of PAs. The NTL levels of types II and VI were not significantly higher than those of Ia and Ib.
The main function of type IV PAs is to protect certain species and habitats. The NTL levels inside type IV PAs and buffer zones were significantly higher than all other types except V. NTL levels within and outside type III PAs were high, second only to type V. NTL levels within and outside of type V PAs were much higher than other types.
The result suggested that the NTL levels were significantly lower within the PAs than that of the surrounding 0–100 km buffer zones. In particular, there was a big difference between inside and outside the boundary of types Ia and III PAs. For most PAs, the surrounding areas with the highest DN value were in the 1–25 km buffer zone.

4.2. Skyglow as a Biodiversity Threat

Skyglow occurs when artificial light is scattered by water, atmospheric molecules, or aerosols and returned to Earth. In this study, we used the NTL data to show the impact of urbanization on protected areas. However, we did not apply any light propagation models to integrate the skyglow effect due to indirect lights. Skyglow is of increasing concern since it is able to multiply and extend the light pollution effect and affect those areas with no direct light pollution. Because of skyglow, the biological impacts of artificial light are not limited to the vicinity of the light source and may spread over much larger extents via several mechanisms. Individual lights may be visible kilometers away from their source, and the addition of artificial skyglow can extinguish such lunar light cycles and permanently remove dark nights from a landscape [59,60]. Therefore, artificial lighting may have an effect on natural ecosystems even when the source of lighting lies kilometers away. Most organisms have evolved molecular circadian clocks controlled by natural day–night cycles. These clocks play key roles in metabolism, growth, and behavior [61]. As the world grows ever-more illuminated, many light-sensitive species will be lost, especially in or near highly illuminated urban areas [62]. Light pollution threatens biodiversity through changed night habits, including organismal movements [63,64]; foraging [65,66,67,68]; interspecific interactions [69]; communication [70,71]; reproduction [72] of insects, amphibians, fish, birds, bats, and other animals; and it can disrupt plants resulting in phenological changes by distorting their natural day–night cycle [63].
Therefore, measures are necessary to prevent or reduce the ecological impact of night-time light pollution. Maintaining and increasing natural unlit areas is likely to be the most effective option for reducing the ecological effects of lighting [73]. From the lighting differences inside and outside the PAs, it can be seen that in the process of human development, PAs have greatly reduced human interference, so PAs are undoubtedly the best place to maintain darkness. More stringent control measures should be implemented within and around PAs, such as limiting the duration of lighting, reducing the trespass of lighting, changing the intensity of lighting, and changing the spectrum of lighting [73,74].
Setting an area surrounding PAs is an ideal option to provide a buffer against the light pollution released by human development. However, it cannot be solved by using only remote sensing data. Different types of PAs are in different natural and socio-economic conditions, and as such, the appropriate buffer radius should also be different. The success of planners in reducing the ecological impacts of light pollution will ultimately depend on an assessment of the critical mechanisms and thresholds that determine those impacts in a particular environment [73]. One possible way is to combine remote sensing data with biodiversity data to explore how biodiversity in protected areas with similar natural conditions respond to different lighting patterns. This may be of great significance for solving light pollution near the PAs and may also contribute to determining how much of a buffer distance should be set around the PAs.

4.3. Limitations of Night-Time Lighting Data

The 1992–2013 NTL data used in this study were derived from six different sensors. Therefore, it is difficult to discriminate whether the small difference between the light data of different years was caused by the sensor or by changes in the field light. Moreover, we are not sure how much data error was caused by the sensor. In this study, we compared the nightlight levels between 1992 and 2013. The time span and the difference in the NTL data were sufficiently large, such that the data error caused by the different sensor could be ignored.
In addition, the light data could not detect negative light growth (i.e., reductions in light) in each PA. To reduce the error caused by the sensor when processing the light data, we assumed that all regions were developing positively, and the DN value was increasing; thus, we assigned the larger DN value of the previous year to the darker pixels corresponding to the following year. If the actual nightlight of a certain area was dark each year, the DN value would remain unchanged. On the global scale, the nightlight within and around each type of PA was increasing each year, and there was no negative growth in the brightness of the PA. However, when assessing individual PAs, there may be negative growth due to increasingly strict regulations and improved awareness of protection, e.g., lighting tools could be replaced by those with a lower brightness and lighting time could be reduced. Therefore, the lighting data may not be suitable for the study of PAs on a small spatial scale, especially in developed countries where the management of PAs is more stringent.

Author Contributions

All authors contributed to this manuscript. Specific contributions include conceptualization, L.F., J.Z., and Y.W.; methodology, J.Z. and L.F.; software, L.F.; validation, L.F. and J.Z.; formal analysis, L.F.; data curation, L.F. and J.Z.; writing—original draft preparation, L.F. and J.Z.; writing—review and editing, L.F., J.Z., Y.W., Z.R., H.Z., and X.G.; visualization, J.Z. and L.F.; supervision, J.Z.; funding acquisition, J.Z.

Funding

This research was funded by the National Natural Science Foundation of China (grant nos. 41771450, 41871330, and 41630749), the Fundamental Research Funds for the Central Universities (grant no. 2412019BJ001), the Foundation of the Education Department of Jilin Province in the 13th Five-Year project (grant no. JJKH20190282KJ), and the Science and Technology Development Project of Jilin Province (grant no. 20190802024ZG).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. UNEP-WCMC; IUCN. Protected Planet Report 2016; UNEP-WCMC: Cambridge, UK; IUCN: Gland, Switzerland, 2016; ISBN 978-92-807-3587-1. Available online: https://www.iucn.org/theme/protected-areas/publications/protected-planet-report (accessed on 12 June 2019).
  2. Wang, Y. Remote Sensing of Protected Lands; CRC Press: Boca Raton, FL, USA, 2011; ISBN 978-1-4398-4187-7. [Google Scholar]
  3. Elvidge, C.D.; Sutton, P.C.; Wagner, T.W.; Ryzner, R.; Vogelmann, J.V.; Goetz, S.J.; Smith, A.J.; Jantz, C.; Seto, K.C.; Imhoff, M.L.; et al. (Eds.) Land Change Science; Kluwer Academic Publishers: Dordrecht, The Netherlands, 2004; pp. 315–328. [Google Scholar]
  4. Crabtree, R.L.; Sheldon, J.W.; Sheldon, J.W. Monitoring and Modeling Environmental Change in Protected Areas: Integration of Focal Species Populations and Remote Sensing. In Remote Sensing of Protected Lands; CRC Press: Boca Raton, FL, USA, 2012; pp. 495–524. [Google Scholar]
  5. IUCN. Shaping a sustainable future. In The IUCN Programme 2009–2012; IUCN: Gland, Switzerland, 2008. [Google Scholar]
  6. Heinen, J.; Hite, K. Protected natural areas. In Encyclopedia of Earth; Cutler, J., Cleveland, C.J., Eds.; Environmental Information Coalition, National Council for Science and the Environment: Washington, DC, USA, 2007. [Google Scholar]
  7. Hoffmann, M.; Hilton-Taylor, C.; Angulo, A.; Böhm, M.; Brooks, T.M.; Butchart, S.H.M.; Carpenter, K.E.; Chanson, J.; Collen, B.; Cox, N.A.; et al. The Impact of Conservation on the Status of the World’s Vertebrates. Science 2010, 330, 1503. [Google Scholar] [CrossRef]
  8. Pereira, H.M.; Leadley, P.W.; Proença, V.; Alkemade, R.; Scharlemann, J.P.W.; Fernandez-Manjarrés, J.F.; Araújo, M.B.; Balvanera, P.; Biggs, R.; Cheung, W.W.L.; et al. Scenarios for Global Biodiversity in the 21st Century. Science 2010, 330, 1496. [Google Scholar] [CrossRef]
  9. Watson, J.E.M.; Venter, O.; Lee, J.; Jones, K.R.; Robinson, J.G.; Possingham, H.P.; Allan, J.R. Protect the last of the wild. Nature 2018, 563, 27. [Google Scholar] [CrossRef] [PubMed]
  10. Prendergast, J.R.; Quinn, R.M.; Lawton, J.H.; Eversham, B.C.; Gibbons, D.W. Rare species, the coincidence of diversity hotspots and conservation strategies. Nature 1993, 365, 335. [Google Scholar] [CrossRef]
  11. Finlayson, M.; Cruz, R.D.; Davidson, N.; Alder, J.; Cork, S.; Groot, R.S.D.; Leveque, C.; Milton, G.R.; Peterson, G.; Pritchard, D.; et al. Millennium Ecosystem Assessment: Ecosystems and Human Well-Being: Wetlands and Water Synthesis; World Resources Institute: Washington, DC, USA, 2005. [Google Scholar]
  12. Gaston, K.J.; Jackson, S.F.; Cantúsalazar, L.; Cruzpiñón, G. The Ecological Performance of Protected Areas. Annu. Rev. Ecol. Evol. Syst. 2008, 39, 93–113. [Google Scholar] [CrossRef]
  13. Walker, R.; Moore, N.J.; Arima, E.; Perz, S.; Simmons, C.; Caldas, M.; Vergara, D.; Bohrer, C. Protecting the Amazon with protected areas. Proc. Natl. Acad. Sci. USA 2009, 106, 10582–10586. [Google Scholar] [CrossRef] [PubMed]
  14. George, W.; Paul, E.; Bean, W.T.; Burton, A.C.O.; Brashares, J.S. Accelerated human population growth at protected area edges. Science 2008, 321, 123–126. [Google Scholar]
  15. Scherl, L.M.; Wilson, A.; Wild, R.; Blockhus, J.; Franks, P.; Mcneely, J.A.; Mcshane, T.O. Can protected areas contribute to poverty reduction? opportunities and limitations. Biodiversity 2010, 11, 5–7. [Google Scholar]
  16. Zhou, Y.; Li, X.; Asrar, G.R.; Smith, S.J.; Imhoff, M. A global record of annual urban dynamics (1992–2013) from nighttime lights. Remote Sens. Environ. 2018, 219, 206–220. [Google Scholar] [CrossRef]
  17. Li, X.; Zhou, Y. Urban mapping using DMSP/OLS stable night-time light: A review. Int. J. Remote Sens. 2017, 38, 6030–6046. [Google Scholar] [CrossRef]
  18. Rodrigues, A.S.; Andelman, S.J.; Bakarr, M.I.; Boitani, L.; Brooks, T.M.; Cowling, R.M.; Fishpool, L.D.; Da, F.G.; Gaston, K.J.; Hoffmann, M. Effectiveness of the global protected area network in representing species diversity. Nature 2004, 428, 640–643. [Google Scholar] [CrossRef] [PubMed]
  19. Bruner, A.G.; Gullison, R.E.; Rice, R.E.; Fonseca, G.A. Effectiveness of parks in protecting tropical biodiversity. Science 2001, 291, 125–128. [Google Scholar] [CrossRef] [PubMed]
  20. Andam, K.S.; Ferraro, P.J.; Pfaff, A.; Sanchez-Azofeifa, G.A.; Robalino, J.A. Measuring the Effectiveness of Protected Area Networks in Reducing Deforestation. Proc. Natl. Acad. Sci. USA 2008, 105, 16089–16094. [Google Scholar] [CrossRef] [PubMed]
  21. Brun, C.; Cook, A.R.; Lee, J.S.H.; Wich, S.A.; Lian, P.K.; Carrasco, L.R. Analysis of deforestation and protected area effectiveness in Indonesia: A comparison of Bayesian spatial models. Glob. Environ. Chang. 2015, 31, 285–295. [Google Scholar] [CrossRef]
  22. Foley, J.A.; Ruth, D.; Asner, G.P.; Carol, B.; Gordon, B.; Carpenter, S.R.; F Stuart, C.; Coe, M.T.; Daily, G.C.; Gibbs, H.K. Global consequences of land use. Science 2005, 309, 570–574. [Google Scholar] [CrossRef]
  23. Mcdonald, R.I.; Kareiva, P.; Forman, R.T.T. The implications of current and future urbanization for global protected areas and biodiversity conservation. Biol. Conserv. 2008, 141, 1695–1703. [Google Scholar] [CrossRef]
  24. Radeloff, V.C.; Stewart, S.I.; Hawbaker, T.J.; Gimmi, U.; Pidgeon, A.M.; Flather, C.H.; Hammer, R.B.; Helmers, D.P. Housing growth in and near United States protected areas limits their conservation value. Proc. Natl. Acad. Sci. USA 2010, 107, 940–945. [Google Scholar] [CrossRef]
  25. Li, X.; Zhou, Y.; Asrar, G.R.; Mao, J.; Li, X.; Li, W. Response of vegetation phenology to urbanization in the conterminous United States. Glob. Chang. Biol. 2017, 23, 2818–2830. [Google Scholar] [CrossRef]
  26. Elvidge, C.D.; Sutton, P.C.; Baugh, K.E.; Ziskin, D.C.; Anderson, S. National Trends in Satellite Observed Lighting: 1992–2009. In Proceedings of the Agu Fall Meeting, San Francisco, CA, USA, 5–9 December 2011. [Google Scholar]
  27. Hsu, F.C.; Baugh, K.E.; Ghosh, T.; Zhizhin, M.; Elvidge, C.D. DMSP-OLS Radiance Calibrated Nighttime Lights Time Series with Intercalibration. Remote Sens. 2015, 7, 1855–1876. [Google Scholar] [CrossRef]
  28. Zhang, Q.; Pandey, B.; Seto, K.C. A Robust Method to Generate a Consistent Time Series from DMSP/OLS Nighttime Light Data. IEEE Trans. Geosci. Remote Sens. 2016, 54, 5821–5831. [Google Scholar] [CrossRef]
  29. Gaston, K.J.; Bennie, J. Demographic effects of artificial nighttime lighting on animal populations. Doss. Environ. 2014, 22, 323–330. [Google Scholar] [CrossRef]
  30. Gaston, K.J.; Duffy, J.P.; Bennie, J. Quantifying the erosion of natural darkness in the global protected area system. Conserv. Biol. 2015, 29, 1132–1141. [Google Scholar] [CrossRef] [PubMed]
  31. Longcore, T.; Rich, C. Ecological Light Pollution. Front. Ecol. Environ. 2004, 2, 191–198. [Google Scholar] [CrossRef]
  32. Bennie, J.; Duffy, J.P.; Davies, T.W.; Correa-Cano, M.E.; Gaston, K.J. Global Trends in Exposure to Light Pollution in Natural Terrestrial Ecosystems. Remote Sens. 2015, 7, 2715–2730. [Google Scholar] [CrossRef]
  33. Freitas, J.R.; Bennie, J.; Mantovani, W.; Gaston, K.J. Exposure of tropical ecosystems to artificial light at night: Brazil as a case study. PLoS ONE 2017, 12, e0171655. [Google Scholar] [CrossRef]
  34. Koen, E.L.; Minnaar, C.; Roever, C.L.; Boyles, J.G. Emerging threat of the 21st century lightscape to global biodiversity. Glob. Chang. Biol. 2018, 24, 2315–2324. [Google Scholar] [CrossRef]
  35. Zhao, M.; Zhou, Y.; Li, X.; Cao, W.; He, C.; Yu, B.; Li, X.; Elvidge, C.D.; Cheng, W.; Zhou, C. Applications of satellite remote sensing of nighttime light observations: Advances, challenges, and perspectives. Remote Sens. 2019, 11, 1971. [Google Scholar] [CrossRef]
  36. Wei, J.; He, G.; Leng, W.; Long, T.; Guo, H. Characterizing Light Pollution Trends across Protected Areas in China Using Nighttime Light Remote Sensing Data. ISPRS Int. J. Geo-Inf. 2018, 7, 243. [Google Scholar]
  37. Xiang, W.; Tan, M. Changes in Light Pollution and the Causing Factors in China’s Protected Areas, 1992–2012. Remote Sens. 2017, 9, 1026. [Google Scholar] [CrossRef]
  38. Davies, T.W.; Duffy, J.P.; Bennie, J.; Gaston, K.J. Stemming the Tide of Light Pollution Encroaching into Marine Protected Areas. Conserv. Lett. 2015, 9, 164–171. [Google Scholar] [CrossRef]
  39. Raap, T.; Pinxten, R.; Eens, M. Artificial light at night causes an unexpected increase in oxalate in developing male songbirds. Conserv. Physiol. 2018, 6, coy005. [Google Scholar] [CrossRef] [PubMed]
  40. Geldmann, J.; Joppa, L.N.; Burgess, N.D. Mapping change in human pressure globally on land and within protected areas. Conserv. Biol. 2015, 28, 1604–1616. [Google Scholar] [CrossRef] [PubMed]
  41. Theobald, D.M.; Miller, J.R.; Hobbs, N.T. Estimating the cumulative effects of development on wildlife habitat. Landsc. Urban Plan. 1997, 39, 25–36. [Google Scholar] [CrossRef]
  42. He, C.; Li, J.; Chen, J.; Shi, P.; Chen, J.; Pan, Y.; Li, J.; Zhuo, L.; Toshiaki, I. The urbanization process of Bohai Rim in the 1990s by using DMSP/OLS data. J. Geogr. Sci. 2006, 16, 174–182. [Google Scholar] [CrossRef]
  43. Liu, Z.; He, C.; Zhang, Q.; Huang, Q.; Yang, Y. Extracting the dynamics of urban expansion in China using DMSP-OLS nighttime light data from 1992 to 2008. Landsc. Urban Plan. 2012, 106, 62–72. [Google Scholar] [CrossRef]
  44. Chen, X.; Nordhaus, W.D. Using luminosity data as a proxy for economic statistics. Proc. Natl. Acad. Sci. USA 2011, 108, 8589. [Google Scholar] [CrossRef] [Green Version]
  45. Pandey, B.; Joshi, P.K.; Seto, K.C. Monitoring urbanization dynamics in India using DMSP/OLS night time lights and SPOT-VGT data. Int. J. Appl. Earth Obs. Geoinf. 2013, 23, 49–61. [Google Scholar] [CrossRef]
  46. Zhang, Q.; Seto, K.C. Mapping urbanization dynamics at regional and global scales using multi-temporal DMSP/OLS nighttime light data. Remote Sens. Environ. 2011, 115, 2320–2329. [Google Scholar] [CrossRef]
  47. Elvidge, C.; Imhoff, M.; Baugh, K.; Hobson, V.; Nelson, I.; Safran, J.; Dietz, J.; Tuttle, B. Night-time lights of the world: 1994–1995. ISPRS J. Photogramm. Remote Sens. 2001, 56, 81–99. [Google Scholar] [CrossRef]
  48. Seto, K.C.; Zhang, Q. Can Night-Time Light Data Identify Typologies of Urbanization? A Global Assessment of Successes and Failures. Remote Sens. 2013, 5, 3476–3494. [Google Scholar]
  49. Lyytimäki, J.; Tapio, P.; Assmuth, T. Unawareness in environmental protection: The case of light pollution from traffic. Land Use Policy 2012, 29, 598–604. [Google Scholar] [CrossRef]
  50. de Miguel, A.S.; Zamorano, J.; Castano, J.G.; Pascual, S. Evolution of the energy consumed by street lighting in Spain estimated with DMSP-OLS data. J. Quant. Spectrosc. Radiat. Transf. 2014, 139, 109–117. [Google Scholar] [CrossRef] [Green Version]
  51. He, C.; Ma, Q.; Li, T.; Yang, Y.; Liu, Z. Spatiotemporal dynamics of electric power consumption in Chinese Mainland from 1995 to 2008 modeled using DMSP/OLS stable nighttime lights data. J. Geogr. Sci. 2012, 22, 125–136. [Google Scholar] [CrossRef]
  52. Wu, J.; He, S.; Peng, J.; Li, W.; Zhong, X. Intercalibration of DMSP-OLS night-time light data by the invariant region method. Int. J. Remote Sens. 2013, 34, 7356–7368. [Google Scholar] [CrossRef]
  53. Li, X.; Zhou, Y. A Stepwise Calibration of Global DMSP/OLS Stable Nighttime Light Data (1992–2013). Remote Sens. 2017, 9, 637. [Google Scholar]
  54. Yue, S.; Pilon, P.; Cavadias, G. Power of the Mann–Kendall and Spearman’s rho tests for detecting monotonic trends in hydrological series. J. Hydrol. 2002, 259, 254–271. [Google Scholar] [CrossRef]
  55. Dudley, N. Guidelines for Applying Protected Area Management Categories; IUCN: Gland, Switzerland, 2008; ISBN 978-2-8317-1086-0. [Google Scholar]
  56. Radeloff, V.C.; Hammer, R.B.; Stewart, S.I. Rural and Suburban Sprawl in the U.S. Midwest from 1940 to 2000 and Its Relation to Forest Fragmentation. Conserv. Biol. 2005, 19, 793–805. [Google Scholar] [CrossRef]
  57. Crooks, K.R.; Soulé, M.E. Mesopredator release and avifaunal extinctions in a fragmented system. Nature 1999, 400, 563–566. [Google Scholar] [CrossRef]
  58. Woodroffe, R. Edge Effects and the Extinction of Populations Inside Protected Areas. Science 1998, 280, 2126–2128. [Google Scholar] [CrossRef] [Green Version]
  59. Davies, T.; Bennie, J.; Inger, R.; Gaston, K. Artificial light alters regimes of natural sky brightness. Sci. Rep. 2013, 3, 1722. [Google Scholar] [CrossRef] [Green Version]
  60. Puschnig, J.; Posch, T.; Uttenthaler, S. Night sky photometry and spectroscopy performed at the Vienna University Observatory. J. Quant. Spectrosc. Radiat. Transf. 2014, 139, 64–75. [Google Scholar] [CrossRef] [Green Version]
  61. Dunlap, J. Molecular Bases for Circadian Clocks. Cell 1999, 96, 271–290. [Google Scholar] [CrossRef] [Green Version]
  62. Hölker, F.; Wolter, C.; Perkin, E. Light Pollution as a Biodiversity Threat. Trends Ecol. Evol. 2010, 25, 681–682. [Google Scholar] [CrossRef] [PubMed]
  63. Peters, A.; Verhoeven, K.J.F. Impact of Artificial Lighting on the Seaward Orientation of Hatchling Loggerhead Turtles. J. Herpetol. 1994, 28, 112–114. [Google Scholar] [CrossRef]
  64. Emma Louise, S.; Gareth, J.; Stephen, H. Street lighting disturbs commuting bats. Curr. Biol. CB 2009, 19, 1123–1127. [Google Scholar]
  65. Rydell, J. Seasonal use of illuminated areas by foraging northern bats Eptesicus nilssoni. Ecography 1991, 14, 203–207. [Google Scholar] [CrossRef]
  66. Buchanan, B.W. Effects of enhanced lighting on the behaviour of nocturnal frogs. Anim. Behav. 1993, 45, 893–899. [Google Scholar] [CrossRef] [Green Version]
  67. Bird, B.L.; Branch, L.C.; Miller, D.L. Effects of Coastal Lighting on Foraging Behaviorof Beach Mice. Conserv. Biol. 2004, 18, 1435–1439. [Google Scholar] [CrossRef]
  68. Santos, C.D.; Miranda, A.C.; Granadeiro, J.P.; Lourenco, P.M.; Saraiva, S.; Palmeirim, J.M. Effects of artificial illumination on the nocturnal foraging of waders. Acta Oecol. 2010, 36, 166–172. [Google Scholar] [CrossRef]
  69. Svensson, A.M.; Rydell, J. Mercury vapour lamps interfere with the bat defence of tympanate moths (Operophtera spp.; Geometridae). Anim. Behav. 1998, 55, 223–226. [Google Scholar] [CrossRef] [Green Version]
  70. Baker, B.J.; Richardson, J. The effect of artificial light on male breeding-season behaviour in green frogs, Rana clamitans melanota. Can. J. Zool. 2006, 84, 1528–1532. [Google Scholar] [CrossRef]
  71. Miller, M.W. Apparent Effects of Light Pollution on Singing Behavior of American Robins. Condor 2006, 108, 130. [Google Scholar] [CrossRef]
  72. Boldogh, S.; Dobrosi, D.; Samu, P. The effects of the illumination of buildings on house-dwelling bats and its conservation consequences. Acta Chiropterol. 2007, 9, 527–534. [Google Scholar] [CrossRef]
  73. Gaston, K.J.; Davies, T.W.; Bennie, J.; Hopkins, J.; Fernandez-Juricic, E. Review: Reducing the ecological consequences of night-time light pollution: Options and developments. J. Appl. Ecol. 2012, 49, 1256–1266. [Google Scholar] [CrossRef] [Green Version]
  74. Kyba, C.C.M.; Hölker, F. Do artificially illuminated skies affect biodiversity in nocturnal landscapes? Landsc. Ecol. 2013, 28, 1637–1640. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Gross domestic product (GDP; unit: trillions of dollars) and global the sum of NTL (SNTL) (1992–2013) from (a) the raw NTL and (b) the calibrated NTL. Electricity consumption (EC; unit: billion kWh) and global SNTL (1992–2013) from (c) the raw NTL and (d) the calibrated NTL. DMSP/OLS—Defense Meteorological Satellite Program/ Operational Linescan System.
Figure 1. Gross domestic product (GDP; unit: trillions of dollars) and global the sum of NTL (SNTL) (1992–2013) from (a) the raw NTL and (b) the calibrated NTL. Electricity consumption (EC; unit: billion kWh) and global SNTL (1992–2013) from (c) the raw NTL and (d) the calibrated NTL. DMSP/OLS—Defense Meteorological Satellite Program/ Operational Linescan System.
Remotesensing 11 02699 g001
Figure 2. The correlation of GDP (a) and EC (b) with the raw NTL and the calibrated NTL. The black line is the 1:1 line.
Figure 2. The correlation of GDP (a) and EC (b) with the raw NTL and the calibrated NTL. The black line is the 1:1 line.
Remotesensing 11 02699 g002
Figure 3. Spatial distribution of NTL in association with global terrestrial protected areas (PAs). The orange-colored areas represent the PAs; the dark blue-colored area (including oceans) were in darkness; the yellow-colored areas had a higher NTL digital number (DN) value than the green-colored areas; the white-colored area had no data.
Figure 3. Spatial distribution of NTL in association with global terrestrial protected areas (PAs). The orange-colored areas represent the PAs; the dark blue-colored area (including oceans) were in darkness; the yellow-colored areas had a higher NTL digital number (DN) value than the green-colored areas; the white-colored area had no data.
Remotesensing 11 02699 g003
Figure 4. (a) Average NTL DN value for every PA and (b) average NTL DN value for every continent. For (a), we used the natural breaks (Jens) method to divide the PAs into three levels according to the average DN value. The green-colored PAs were those with the average value of 0, the yellow-colored PAs had a lower than average DN value, and the red-colored PAs had a higher than average DN value. For (b), the ranking of the average NTL DN value of PAs in each continent gave the order as Europe > Asia > North America > South America > Africa > Oceania.
Figure 4. (a) Average NTL DN value for every PA and (b) average NTL DN value for every continent. For (a), we used the natural breaks (Jens) method to divide the PAs into three levels according to the average DN value. The green-colored PAs were those with the average value of 0, the yellow-colored PAs had a lower than average DN value, and the red-colored PAs had a higher than average DN value. For (b), the ranking of the average NTL DN value of PAs in each continent gave the order as Europe > Asia > North America > South America > Africa > Oceania.
Remotesensing 11 02699 g004
Figure 5. (a) Average trend of NTL of the global terrestrial PAs and (b) average trend of every continent. For (a), we used the natural breaks (Jens) method to divide the PAs into three levels according to the average trend of NTL. The blue-colored PAs had an average trend of 0, the yellow-colored PAs had a lower than average trend, and the red-colored PAs had a higher than average trend. For (b), the ranking of the average NTL trend of PAs in each continent gave the order as Europe > Asia > North America > South America > Africa > Oceania.
Figure 5. (a) Average trend of NTL of the global terrestrial PAs and (b) average trend of every continent. For (a), we used the natural breaks (Jens) method to divide the PAs into three levels according to the average trend of NTL. The blue-colored PAs had an average trend of 0, the yellow-colored PAs had a lower than average trend, and the red-colored PAs had a higher than average trend. For (b), the ranking of the average NTL trend of PAs in each continent gave the order as Europe > Asia > North America > South America > Africa > Oceania.
Remotesensing 11 02699 g005
Figure 6. NTL level on different buffers for each type of PA. Columns with different colors represent the interior of PAs and different buffer zones.
Figure 6. NTL level on different buffers for each type of PA. Columns with different colors represent the interior of PAs and different buffer zones.
Remotesensing 11 02699 g006
Figure 7. NTL level for each type of PA in buffer zones. Columns with different colors represent different types of PAs. The light blue shadowed areas represents the mean DN values of PAs and buffers.
Figure 7. NTL level for each type of PA in buffer zones. Columns with different colors represent different types of PAs. The light blue shadowed areas represents the mean DN values of PAs and buffers.
Remotesensing 11 02699 g007
Figure 8. Changes within PAs and in different buffer zones in the NTL level of different types of PAs in the time series. The numbers in the grid represent the average NTL DN values. The color of the grids is from blue to red, and the corresponding values are from small to large. (a) strict nature reserve (Ia), (b) wilderness area (Ib), (c) national park (II), (d) natural monument or feature (III), (e) habitat/species management area (IV), (f) protected landscape/seascape (V), and (g) protected area with sustainable use of natural resources (VI).
Figure 8. Changes within PAs and in different buffer zones in the NTL level of different types of PAs in the time series. The numbers in the grid represent the average NTL DN values. The color of the grids is from blue to red, and the corresponding values are from small to large. (a) strict nature reserve (Ia), (b) wilderness area (Ib), (c) national park (II), (d) natural monument or feature (III), (e) habitat/species management area (IV), (f) protected landscape/seascape (V), and (g) protected area with sustainable use of natural resources (VI).
Remotesensing 11 02699 g008aRemotesensing 11 02699 g008bRemotesensing 11 02699 g008c
Figure 9. The NTL level of different types of PAs in their interior and surrounding buffer zones. The blue shadowed area represents the fact that the 1–10 km buffer zone had the highest NTL level.
Figure 9. The NTL level of different types of PAs in their interior and surrounding buffer zones. The blue shadowed area represents the fact that the 1–10 km buffer zone had the highest NTL level.
Remotesensing 11 02699 g009
Figure 10. Growth rate of NTL within every type of PA and the buffer zones.
Figure 10. Growth rate of NTL within every type of PA and the buffer zones.
Remotesensing 11 02699 g010
Figure 11. Trends within every type of PAs and of buffers.
Figure 11. Trends within every type of PAs and of buffers.
Remotesensing 11 02699 g011
Table 1. Seven types of terrestrial protected areas defined by IUCN.
Table 1. Seven types of terrestrial protected areas defined by IUCN.
IUCN PAs TypeIUCN PAs CodeNumber of PAsArea (km2)Description
Strict Nature ReserveIa11,9213,874,328PAs that are strictly set aside to protect biodiversity and also possibly geological/geomorphological features, where human visitation, use, and impacts are strictly controlled and limited to ensure the protection of the conservation values.
Wilderness AreaIb31141,152,403PAs that are usually large unmodified or slightly modified areas, retaining their natural character and influence, without permanent or significant human habitation.
National Park55236,178,074Large natural or near natural areas set aside to protect large-scale ecological processes, along with the complement of species and ecosystems characteristic of the area, which also provide a foundation for environmentally and culturally compatible spiritual, scientific, educational, recreational, and visitor opportunities.
Natural Monument or Feature14,317434,460PAs set aside to protect specific natural monuments, which can be a landform; sea mount; submarine cavern; geological feature, such as a cave; or even a living feature, such as an ancient grove. They are generally quite small protected areas and often have a high visitor value.
Habitat/Species Management Area56,8367,250,215PAs aiming to protect particular species or habitats and management reflects this priority. Many category IV PAs will need regular, active interventions to address the requirements of particular species or to maintain habitats, but this is not a requirement of the category.
Protected Landscape/SeascapeV44,9154,351,409PAs where the interaction of people and nature over time have produced the area of distinct character with significant ecological, biological, cultural, and scenic value.
Protected Area with Sustainable use of Natural Resources717712,757,697PAs that conserve ecosystems and habitats, together with associated cultural values and traditional natural resource management systems. They are generally large, with most of the area in a natural condition, where a proportion is under sustainable natural resource management and where low-level, non-industrial use of natural resources compatible with nature conservation is seen as one of the main aims of the area.
Table 2. NTL level of interior and buffer zones for the seven types of PAs.
Table 2. NTL level of interior and buffer zones for the seven types of PAs.
PA TypeInterior0–1 km1–5 km5–10 km10–25 km25–50 km50–100 km
Ia0.0981.5722.3722.5982.4782.3441.916↑
Ib0.0390.4960.9201.3801.8972.3172.452↑
II0.1471.1291.6021.9201.9321.7021.427
III0.5944.868↑6.289↑6.0384.734↑3.272↑2.013↑
IV0.6521.3173.593↑3.3012.6312.0781.456
V3.184↑3.270↑7.151↑6.840↑5.418↑3.610↑2.357↑
VI0.1600.1691.2171.4831.7391.8111.692
Ave.0.6961.8323.3073.3662.9752.4481.902

Share and Cite

MDPI and ACS Style

Fan, L.; Zhao, J.; Wang, Y.; Ren, Z.; Zhang, H.; Guo, X. Assessment of Night-Time Lighting for Global Terrestrial Protected and Wilderness Areas. Remote Sens. 2019, 11, 2699. https://doi.org/10.3390/rs11222699

AMA Style

Fan L, Zhao J, Wang Y, Ren Z, Zhang H, Guo X. Assessment of Night-Time Lighting for Global Terrestrial Protected and Wilderness Areas. Remote Sensing. 2019; 11(22):2699. https://doi.org/10.3390/rs11222699

Chicago/Turabian Style

Fan, Liangxian, Jianjun Zhao, Yeqiao Wang, Zhoupeng Ren, Hongyan Zhang, and Xiaoyi Guo. 2019. "Assessment of Night-Time Lighting for Global Terrestrial Protected and Wilderness Areas" Remote Sensing 11, no. 22: 2699. https://doi.org/10.3390/rs11222699

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