Reforestation and Deforestation in Northern Luzon, Philippines: Critical Issues as Observed from Space

: Among the richest in biodiversity globally has been the Philippine rainforest, which used to cover about 90% of the country’s land area. During the last few decades, the forest cover has been reduced to less than 10% of the original, only a fraction of which is old ‐ growth forest. The negative impacts of deforestation led to the launching of the National Greening Program (NGP) that involved the planting of more than a billion seedlings over a few million hectares of land from 2011 to 2016. To assess the success of the NGP, satellite data from Landsat and the Moderate Resolution Imaging Spectroradiometer (MODIS) were analyzed before, during, and after the NGP. Reforestation in the NGP sites was examined concurrently with observed deforestation in Luzon using forest loss data derived from Landsat for the period 2001 to 2018. The results show that losses declined from 2011 to 2015 but increased from 2016 to 2018. Because of such losses, the net effect is a balance of reforestation and deforestation or no significant gain from the NGP. Case studies were done in three sites in the Sierra Madre forest, where half of the remaining old ‐ growth forest is located, using a combination of Landsat and Very High Resolution (VHR) data. The Landsat data were classified into closed forest, open forest, and other vegetation cover types. The conversion from one vegetation cover type to another was evaluated through the use of the Sankey Diagram. While some non ‐ forest types became open or closed forests, the loss of open or closed forests is more pronounced. VHR data reveal critical issues happening within the NGP sites during the NGP period. More comprehensive data from MODIS also confirm that there was no significant increase in the forest cover in Luzon, Sierra Madre, and Cordillera from 2001 to 2018.


Introduction
The Sierra Madre and Cordillera forests in Northern Luzon were among the most extensive rainforests in the Philippines for a long time. They host highly biodiverse species of plants and animals, many of which are endemic [1][2][3][4]. Unfortunately, these forests also provided the muchcoveted timber for housing and other uses and became the site of many irresponsible logging activities. They also became sites of swidden agriculture and mining. During the 1970 to 1990 period, the percentage of remaining rainforest in the country was drastically reduced from about 70% of what it was in the 1900s to less than 10% [5]. Despite many laws and regulations aimed at protecting the remaining forest, including converting many areas to national parks, natural monuments, watershed forest reserves, marine reserves, and protected landscapes and seascapes, forest loss continued to occur up to the present [6][7][8]. The ongoing deforestation in the country has been a big issue considering the many established adverse effects, including soil erosion, landslides, flooding, loss of biodiversity, and watershed degradation. The impacts are more profound considering the critical role of forests in the sequestration of atmospheric carbon dioxide and in keeping the richness of current biodiversity intact and the ecosystem resilient to climate change. Recent studies also indicate that proper reforestation has a cooling effect of about 1 to 3 °C [9].
The national government has embarked on important programs, such as the "National Greening Program (NGP)" created by Executive Order No. 26 on 24 February 2011, in part to make amends to previous overexploitation and restore the much-depleted forest in the country to a more acceptable state [10]. The program has been regarded as the most ambitious and well thought out reforestation project ever undertaken by the Department of Environment and Natural Resources (DENR). Posted as the accomplishment at the end of the program in 2016 is the planting of 1.3 billion seedlings in around 1.7 million hectares and the creation of some four million jobs that benefitted about 558,323 individuals as hired workers [11]. It made the incentive to continue with the program that is now referred to as the Enhanced National Greening Program (ENGP) by the current administration to rehabilitate all the forestlands, estimated at 7.1 million hectares, from 2016 to 2028 [10]. Such plans are laudable, considering that the Philippines is among the four Asian countries with the largest annual loss of forested areas at 157,000 ha per year from 2000 to 2005 [12]. However, the NGP has not been without controversy since there have not been good ground confirmations of the reported successes in the areas of activities.
To have an improved assessment of the impact of NGP activities, it is necessary to have comprehensive monitoring of the forest cover changes associated with the program. Satellite data have been used in many global and regional studies of forest cover and vegetation [13][14][15][16]. The merit of satellite data is that it is comprehensive and provides continuous and consistent coverage of the region. It thus allows the assessment of seasonal, interannual, and relatively long-term changes on a large scale. Moreover, with the advent of very high-resolution data (about 1-3 m), direct confirmation of the success of the program to improve the status of the forest cover can be performed.
This study's primary goal is to use satellite observations to gain insights into the overall impacts of the NGP on the current state of the forest cover in Northern Luzon with particular emphasis on the Sierra Madre forest and, to a lesser extent, the Cordillera forest. The Sierra Madre forest is especially important to study because it is the site of about half of the surviving old-growth forest in the country [6]. The old-growth forest has been referred to as the only hope of preserving many of the remaining species in the forest [17]. In this regard, three study sites were selected, as described below, to document the continuation of losses and gains while the NGP is in progress and elucidate on the issues that hinder the program to succeed. Historical satellite data from 2001 to 2019 is used to better understand the status of forest cover in Luzon and the associated transformation of the forest due to the NGP. We also distinguish between closed forest, which is the undisturbed region of the forest where most of the old-growth forests are located, and open forest, which had been accessible to humans and may have roads and communities within the region. The NGP is an admirable program meant to provide a reversal of historical forest losses. With the use of satellite observations, we validate the program's success and, at the same time, examine the critical issues that need to be addressed.

Study Areas and NGP Sites
The state of forest cover in Northern Luzon in 2000 is depicted by the color-coded image presented in Figure 1. The image represents the percentage of tree cover derived from Landsat by Hansen et al. [14]. It depicts the distribution of forest cover in the Sierra Madre (right or eastern side of the image), which stretches from the north (top right) to the south and the Cordillera Mountains (central left) that cover much of the western region of Northern Luzon, but converges with the Sierra Madre at its southern end. The forests in the Zambales mountains are also shown in the lower left and include those in Mount Arayat and Mount Pinatubo. The locations of areas that are supposed to be off-limits to encroachment or deforestation, such as national parks, watershed reserves, and protected landscapes, are also shown. To protect the remnant of old-growth forest, a large fraction of the Sierra Madre forest has been established as national parks, such as the Fuyot-Springs National Park in the north, the Northern Sierra Madre Park in the middle, and the Aurora Memorial National Park towards the south. The Sierra Madre mountain range in the Philippines is also the longest in the country, extending about 500 km from northern Cagayan to Quezon Province [17]. It has protected inland communities from the direct impact of typhoons coming mainly from the Pacific Ocean. On the other hand, the Cordillera forest covers a broader area than the Sierra Madre but is not as long. Cordillera's forest cover is already heavily depleted and has become very patchy mainly because it has been more accessible and vulnerable to logging, swidden agriculture, and other activities.
The distribution of NGP activities is depicted as red-colored areas on the map shown in Figure  1. The map shows that most of the NGP sites are concentrated in the western part of the Cordillera Mountains, the southern part of Sierra Madre, and the west side of the Zambales mountains. It is also apparent that the focus in the Cordillera mountains is the region where the forest has been highly degraded, and similarly for the Sierra Madre and Zambales forests. Although the greening activities are less extensive in the Sierra Madre, the efforts to restore the forest in the region are significant and are focused on expanding the forest's extent on the western side. The three sites in the Sierra Madre, as mentioned earlier, are indicated by the light blue squares in the figure and were selected for case studies to illustrate in detail how the open and closed forest cover changed during the NGP period. The sites have approximately the same area and are chosen because of significant NGP activities in these sites and ample coverage of both forested and non-forested areas.

NDVI Data from Satellites
One of the most useful vegetation parameters derived from satellite data is the Normalized Difference Vegetation Index (NDVI), as first reported by Tucker [18]. The NDVI is estimated using the ratio (ρNIR − ρR)/(ρNIR + ρR), where ρNIR and ρR are the atmospherically corrected reflectance of the near-infrared (NIR) and red bands, respectively. It takes advantage of the sensitivities of red and near-infrared to the chlorophyll content, canopy density, and leaf water content of the observed surface. A forest canopy in the tropics typically has NDVI values ranging from 0.6 to 0.9 [19]. Lower NDVI values correspond to agricultural areas, grasslands, or bare soil, while negative NDVI could be either water surfaces or clouds [20]. Moreover, tropical forests do not exhibit strong seasonal patterns; hence, its corresponding NDVI values do not vary much within a year. This phenomenon is illustrated in Figure 2a   To investigate the trends and interannual changes in forest cover related to the NGP and other forest activities, a time series of annual NDVI composite maps is generated using MODIS and Landsat reflectance data. Although a MODIS tree cover product (MOD44B) is available, NDVI data (MOD13Q1) [19] are used in this study as a proxy for forest cover (NDVI ≥ 0.8) because of accuracy issues, as reported by Sexton et al. [21]. Meanwhile, the Landsat-based tree cover product is sparse and does not cover the entire period under investigation. The MODIS NDVI at 250 m and Landsat NDVI at 30 m together provide information about the forest cover canopy that can be used to assess forest loss and gain through time [14,[22][23][24][25]. Despite having a coarse resolution, these sensors can detect changes in forest and plantation areas because these sites are relatively homogenous. Furthermore, while the reflectance of seedlings is difficult to discriminate at its early stages, over time, they grow into young trees with a more expansive canopy that can be better detected by the satellites.

Data on Forest Loss
Information on annual global forest cover loss, as derived from Landsat data, has been developed and released for public use by Hansen et al. [14]. The data have been widely used because of their global coverage, reasonably good resolution (30 m), and relatively long historical record, which allows the study of interannual changes in the forest cover [8,22,23]. In this study, version 1.7 of the Hansen dataset was used to assess losses in the forested regions in the study areas of interest in Luzon for each year from 2001 up to 2018. A forest cover, as defined by Hansen et al. [14], refers to a canopy closure with vegetation taller than 5 m in height, while forest loss refers to standreplacement disturbance or a change from forest to non-forest state within a 30 m × 30 m pixel scale. These criteria pose some of the known limitations of the dataset, resulting in the inclusion of plantation areas as part of the forest and the underestimation of forest loss estimates [26]. Part of these issues was addressed by using a co-registered land classification map provided by the National Mapping and Resource Information Authority (NAMRIA) [27] to identify open and closed forests in the Hansen data. The categories of open and closed forest refer to the aggregate of broadleaf, coniferous, and mixed forest. NAMRIA defines open forest as having discontinuous tree layer formation with coverage of between 10 to 40 percent of the ground area. On the other hand, the closed forest has tree formations in the various storeys and the undergrowth covering more than 40 percent of the ground area and does not have a continuous dense grass layer. This definition of forest is consistent with that of the Food and Agriculture Organization (FAO), which describes the forest as land with an area of more than 0.5 hectares and tree crown cover of more than 10 percent. The trees should also reach a minimum height of 5 m at maturity, consistent with the definition of forest tree cover used in determining forest cover loss [14].

Forest Cover Change Analysis
To better evaluate whether the NGP has led to the desired increase in forest cover, a forest cover change analysis is done to compare the forest state before, during, and after NGP. Surface reflectance data from Landsat 5, 7, and 8 are used to generate forest cover classification maps. To obtain complete coverage of the land area, three-year composite maps are derived to characterize the pre (2008-2010), during (2014-2016), and post (2017-2019) status of the forest cover in the study sites. The data from 2011 to 2013 were not analyzed because the NGP just got started at this time, and its impact may not be as apparent as in later years. All the band values are normalized using the coefficients published by Roy et al. [28]. After this, cloud masking is applied before the image composites are generated using the statistical median.
The composite maps are then classified using four classes: closed forest, open forest, sparse vegetation, and non-vegetation. The closed and open forest classes follow the definition described earlier in Section 2.3. The sparse vegetation class constitutes agriculture areas and grasslands, while the non-vegetation category encompasses built-up areas and bare land. The classification is performed in the statistical software, R [29], using the raster package [30] to process Landsat data and the randomForest package [31] for Random Forest classification. The data from blue, green, red, NIR, short-wave infrared (SWIR) 1, and SWIR 2 bands, and NDVI, are used as input layers. Due to the lack of available ground data, the training and validation data are selected using pseudo-invariant features for the available period. These invariant features are identified by calculating the standard deviation per pixel. Calibration coefficients are first applied to minimize the effects of sensor differences [32]. NDVI is selected as the indicator of pseudo-invariant features as it is less sensitive to changes in topography, which is a common characteristic in all sites [33]. Pixels, where the standard deviation in NDVI is above 0.01, is masked out before gathering training data. Additional training data of each site are selected around its vicinity to compensate for the lack of training data for some classes. In contrast, all the validation data were gathered within the sites. Validation was performed using high-resolution images available on Google Earth. The overall accuracy of the classifier is at least 77%, with a Kappa of 0.66. From the classified maps, Sankey diagrams are generated using the googleVis package [34] to determine and visualize the changes in forest cover before, during, and post NGP.

Estimates of Forest Loss from Landsat
A color-coded map of the land cover of Luzon indicating the general locations of open and closed forests, cropland, mangroves, grassland, and shrubs in 2010, as reported by NAMRIA [27], is presented in Figure 3a. The classification map provides the means to gain insights into the forest cover condition in areas considered as open and closed forest areas. Figure 3a shows a much-depleted closed forest and a relatively more extensive open forest. In the Sierra Madre forest, the closed forest cover is confined mainly in the northern and southern portions. In the Cordillera mountains, the closed forest is not as extensive and is located primarily in the north and central regions. There are also patches of closed forest in the Zambales mountains. A large fraction of the vegetation cover shown in Figure 3a is areas of forest losses, the sizes of which are displayed at eight times the actual size for improved visibility. The data provided by Hansen et al. [14] indicate vast areas of forest losses in the Sierra Madre forest during the 2001 to 2010 period, while losses in the Cordillera forest were dominant during the 2011 to 2018 period. Interestingly, in some places, the data points of losses are scattered in different places, while some are clustered and adjacent to each other. This implies that there are deforestation activities that are done by low budget operators or farmers with relatively small-scale activities while, at the same time, there are large-scale operators that concentrate their activities in specific sections of the forest. A key concern is that some of the forest loss regions are located inside the closed forest and even in protected areas. As noted below, losses are also found even in areas where NGP sites are located. Due to the reported underestimation in the Hansen product [26], the loss reported in this study is the minimum possible, and the actual loss may even be more serious.
To better assess the magnitude of deforestation and reforestation activities, charts of yearly forest loss in Luzon for both closed (dark green) and open (light green) forests are presented in Figure  3b Such an increase in the loss of tree cover during the 2016 to 2018 period is discouraging because it occurred at the end of the first phase and the beginning of the second phase of the national greening program. Quantitatively, it is also important to note that the rate of losses is higher for the closed forest (i.e., 893 ha per year) than for the open forest (i.e., 388 ha per year). This is again a big concern because the much-valued closed forest that is closely connected to the rich biodiversity of the region is being reduced more indiscriminately than the open forest. For completeness, a similar chart is presented in Figure 3c but only in NGP sites. It appears that the unexpected phenomenon depicted in Figure 3b has also been occurring within the NGP sites. These results are disturbing since the goals of a very visible program like NGP that is intended to do reforestation were being negated concurrently by deforestation activities that are likely illegal and executed by unscrupulous entrepreneurs and other entities.

Regional and Large Scale Assessments of Greening
The availability of continuous coverage of satellite NDVI data from MODIS (which is used as a proxy for forest cover, as discussed previously) provides the opportunity to study regional and largescale forestry changes in the country. Figure 4a,b, and c show estimates of NDVI yearly averages (in black) as well as reductions in forest area annually (in red) as derived from Hansen's forest loss data and cumulative forest gain (in green) if the NGP areas are assumed to be wholly reforested in study sites 1, 2 and 3 for the period from 2011 to 2016. The NDVI plot in site 1 (Figure 4a) shows a slight increase from 2011 to 2013 but is followed by a slight decrease after that, suggesting very minimal or no gain in forest cover. The forested area is also shown to be initially 42,577 ha in 2011 but was reduced to 42,128 ha in 2016, resulting in a net loss of 449 ha, using the loss data of Hansen et al. [14]. On the other hand, if the NGP sites were planted completely by trees during the period, there would have been a gain of about 2770 ha. Overall, there should have been a net gain of about 2300 ha or a 5.4% increase. However, the NDVI data do not show a significant trend.
In Site 2, the NDVI plot in Figure 4b   If the NGP sites were entirely planted by trees as assumed in the analysis above, the NDVI data should have shown a significant increase in all sites, but this assumption is likely too optimistic since reports have shown otherwise [35]. However, it is apparent that even if the NGP was correctly implemented, the percentage change in forest cover, as indicated in the study areas, is about 20% or less. Since these results may not represent what is typical for the entire country, a similar analysis was done to the whole of Sierra Madre forest, the Cordillera forest, and the rest of Luzon Island. The results are presented in Figure 4d-f. The yearly changes in NDVI show a slight increase from 2011 to 2012, but became relatively stable after that. Quantitatively, using the Hansen et al. [14] data, the Sierra Madre forest cover shows a decrease from 1,764,000 ha to 1,747,000 ha, for a net loss of 17,000 ha. Meanwhile, a 100% implementation of reforestation in the NGP sites would cause an increase of 105,631 ha, resulting in a net gain of 88,631 ha. On the other hand, the forest cover in Cordillera forest went down from 1,751,000 ha to 1,730,000 ha, for a net loss of 21,000 ha, while the NGP project has a potential gain of 13,000 ha forest cover. For the rest of Luzon, the forest cover went down from 1,604,000 to 1,581,000 ha, for a loss of 23,000 ha, while the NGP project could provide a gain of 195,000 ha and an overall net gain of 172,000 ha. These results are not reflected in the trend of NDVI values presented in Figure 4e,f.
To better evaluate how the forest cover has changed for a more extended period, MODIS/NDVI data from April 2000 to April 2018 was used again as a proxy for forest cover, and results are presented in Figure 5. It is apparent that in the study sites 1, 2, and 3, there is significant interannual variability but, overall, trends are negligible at 0.002 per decade for sites 1 and 2, and 0.001 per decade for site 3. It is also interesting that there is a dip in 2010 before the NGP started and in 2011 during the onset of the NGP period. This means that if only data from 2010 to 2018 are considered, there would be a significant positive trend and it is only through the use of a longer data set that more reliable information about the actual trend can be derived.  The interannual variability and trend are also evaluated at a much larger scale using data from the entire Sierra Madre forest, Cordillera, and the rest of Luzon. The results are presented in Figure  5d-f, respectively. This time, the trends are slightly more positive but still negligible at 0.003 for the three regions. Interestingly, there is an even more persistent dip in 2010 and 2011. The results from the larger scale areas are consistent with those from the three study areas in the Sierra Madre forest.
The impact of the NGP is further assessed by looking at the changes in NDVI values before and after the first phase of the NGP. Figure 6a shows the composite NDVI maps from 2008 to 2010, representing the state of the forest cover in Luzon just before the NGP started, while Figure 6b shows the composite values for the period 2017-2019, representing the state of the forest cover after the first phase of the NGP. Visual inspection of the two images indicates a remarkably similar distribution of NDVI values. The changes in NDVI values from the earlier period to the later period are represented by the difference map shown in Figure 6c. The difference map indicates no change (yellow) in the generally forested areas of Sierra Madre, Cordillera, and Zambales forest. This is consistent with the results discussed previously. Some relatively minor changes are apparent both in the middle of the forest and at the edges. The more significant changes are shown to occur in the agricultural and grassland regions, but the changes are almost an even distribution of positive and negative changes. Overall, the maps indicate minimal improvement in the forest cover during the two periods.

Case Studies in Sierra Madre Forest Using Landsat
To gain insights into the forest cover changes before, during, and after phase 1 of the NGP program, the three study sites in the Sierra Madre forest are analyzed in greater detail. As indicated earlier, the variability and trends in these study sites are similar to those observed from the Sierra Madre and Cordillera forest and the rest of Luzon. In this section, classified vegetation cover from Landsat before, during, and after the first phase of the NGP program in the three sites are evaluated. Landsat data from 2008 to 2010 are used to represent the forest cover before the NGP period, while Landsat data from the 2014 to 2016 period are used to represent the forest cover during the first phase of the NGP period. As indicated earlier, data from 2011 to 2013 are not used because it is an early period of the program when changes due to NGP are only moderate. In addition, Landsat data during the 2017 to 2019 period are used to represent the post-NGP period.
The Landsat images are classified into closed forest, open forest, sparse vegetation, and nonvegetation using the technique described earlier. For each study site, a set of classified images has been generated for the three different periods, and it is apparent that transformations are going on, as part of one class is transformed into another class. Tracking down such conversions is complicated but can be simplified through the use of the Sankey Diagram, as discussed earlier. In this study, the diagram is divided into three columns representing the pre-NGP, NGP, and post-NGP periods. In each column, color-coded representation of the area of each class is presented. From one column to the next, the diagram shows how the area for each class changes in magnitude as part of the area is converted to the other classes.
The Sankey Diagram and the associated set of classified images for Site 1 are presented in Figure  7.

Monitoring NGP Sites Using VHR Images
To gain additional insight into the observed changes, a sequence of very high-resolution (VHR) images is used to assess qualitatively, but in good spatial detail, how the vegetation cover gets transformed from one period to another. The VHR images were accessed through Google Earth Pro. Enhancements were applied to make the images uniformly calibrated. The resolution of the RGB images (at about 1 to 3 m) is good enough to discriminate the barren areas and the different vegetation cover types. Figure 10  also became barren on 5 April 2016 at the end of Phase 1 of the NGP, which is consistent with the forest loss shown earlier in Figure 3. Meanwhile, several vegetation patches similar to area B were lost on 5 April 2016, as shown in Figure 10c. The area labeled C is shown to have similar vegetation cover throughout the period, but surrounding areas became barren on 5 April 2016, as indicated in Figure 10c. Burn scars are evident in some of the patches of barren land. Overall, the images from March 2014 to April 2016 show a general decline in the vegetation cover during NGP, but a recovery was observed (except in A) in July 2019 that is likely because of ENGP, but it is not much of an improvement compared to 2014. The set of VHR images shown in Figure 11, representing selected NGP areas within Site 2, illustrates an even more serious issue. The image in Figure 11a taken on 30 January 2010 depicts a relatively robust vegetation condition before the start of the NGP, but by 13 May 2014, which is during the NGP, much of the vegetation is gone, as shown in Figure 11b. This is quite contrary to what has been expected. In particular, the areas labeled A, B, and C, and other areas, became practically barren from 2010 to 2014. Some improvements are apparent by 23 February 2015 ( Figure  11c), and more significantly by 23 March 2017 (Figure 11d), especially for A and B. However, the gains in A and B, as well as in other areas, were lost again only a year later on 18 June 2018, as depicted in Figure 11e. At a closer look, it is evident that some areas surrounding the NGP sites are also used for farming. A restoration was again implemented, as observed on 17 February 2019, and shown in Figure 11f. It appears that new trees are being planted in A and B, but although there is an improvement from 2018 to 2019 in A, a similar improvement is not apparent in B. The reforestation strategy for NGP does not appear consistent, and for NGP to succeed, more effort should be spent on preventing losses during the program.  Figure 12a depicts a closed canopy cover at the start of the NGP, but patches of denuded forest started to emerge in the middle of the NGP, as shown in Figure 12b,c, which shows that these areas appear to have recovered in 2018 but only to be lost again and even expanded a year later, as shown in Figure 12d. One of the overarching issues of interest is whether there is a plan to ensure that the seedlings have a good chance to survive after planting. To illustrate what is going on in the ground, Figures 13  and 14 show VHR images from two NGP areas in Sierra Madre and Cordillera, respectively. In addition to the VHR images, geotagged photos taken by DENR of the actual plantation areas are also shown. Figure 13a is a photo that shows a densely vegetated area covering the left side and an open vegetation area on the right side with a single identifiable plant or seedling in the middle that may be accompanied by smaller plants. A series of VHR images covering the site where the geotagged photo was taken is shown in Figure 13b-d. Figure 13a, which was taken on 15 April 2012, indicates some newly planted areas that appear to have grown by 2015, as shown in Figure 13c. However, it became barren again in 2018, as depicted in Figure 13d. The photo in Figure 13a demonstrates the vulnerability of plants in open vegetation areas and the difficulty of tree seedlings to survive unless they are under proper care and attention during early stages that require the combined use of water and fertilizer, as well as pest and disease control. The corresponding VHR satellite images seem to support this finding. It is important to indicate that not all NGP sites have problems similar to those indicated previously. An example of a case where NGP was properly implemented is shown in Figure 14. The geotagged photo in Figure 14a shows relatively grown-up seedlings of pine trees in the NGP area as of 10 March 2014. A VHR image of the general location as of 22 April 2010 and shown in Figure 14b illustrates an almost barren area in the region. On 6 January 2014, about three months before the photo was taken, it is apparent that the area has much more robust vegetation, as indicated in Figure  14c. The vegetation in the entire region has also improved significantly. In this case, NGP succeeded in transforming a degraded area into a much more desirable vegetated and forested area.

Conclusions
The National Greening Program (NGP) that started in 2011 is an ambitious but a much-needed program that involves large scale restoration of forest cover in the country with the planting of more than a billion seedlings of trees. The goals of the NGP have strong scientific merit considering that the forest cover in the country has been drastically reduced to a highly deplorable state. The program has also helped the impoverished farmers and residents in the region do the planting of seedlings consisting of forest trees like Narra, Lauan, and Mahogany, and trees that provide immediate economic benefits like fruit trees, cacao, coffee, rubber, and bamboo. The reported success of the program has led to the extension of the activities to 2028.
However, the results of this study indicate that such a report is premature and may be difficult to justify. In particular, the results of this study suggest that despite the massive number of seedlings reported to have been planted, the state of the forest cover, at least in the greater Luzon region, does not show significant improvement on a large as well as regional scale. Analysis of Landsat data reveals significant losses in both open and closed forest cover during the same period when NGP was ongoing. This is alarming, especially since the extent of losses even increased significantly towards the end of phase 1 of the NGP program.
Satellite data indicate that the implementation of the NGP has not been done as effectively as desired. Time-series studies of forest cover in three study areas and the entire Sierra Madre, Cordillera, and Northern Luzon area reveal that there has not been a monotonical increase in vegetation or forest cover as expected from the highly touted NGP program. The results of trend analysis indicate practically no trend from the start of the NGP period to 2018. The use of MODIS data for more comprehensive coverage and longer record length also yielded a lack of increase in vegetation and forest cover.
A case study in the Sierra Madre forest using classified Landsat data and Sankey Diagrams shows modest but significant changes in the closed and open forest but not necessarily in the same direction. In Site 1, there was a slight increase in closed forest cover of 3% while the open forest suffered a loss of about 10%. In Site 2, there was a reversal in that the closed forest cover declined by about 9% while the open forest cover increased by about 13%. In Site 3, closed forest cover declined by 22% while the open forest increased by 31%. The net loss of open forest in Site 1 of about 1305 ha is acceptable, especially since there is a gain of 625 ha in closed forest. However, although the results in Site 2 and 3 show net increases in forest cover, the significant losses of 9% for Site 2 and 22% in site 3 in the extent of closed forest are disappointing. High-resolution data were used to confirm unambiguously that deforestation was going on while NGP was being implemented. Many seedlings were observed to be planted in open areas where it is unlikely that the plants get the kind of attention and care needed to survive.
The use of satellite data in this study is confined to optical sensors only, which are known to be vulnerable to cloud cover. This is generally acceptable because the forest cover does not change much over a relatively short period and clouds are dynamic with cloud-covered areas eventually becoming cloud-free over a reasonable period. For near real-time monitoring, the use of microwave data such as those from SAR, radar altimeter, and passive microwave radiometers would be more useful. For improved resolution in a limited area, the use of aircraft data and drones carrying similar sensors would be desirable. Detailed information about the forest cover, including canopy, tree type, and tree height, can be obtained through other types of sensors such as SAR, radar altimeter, and Lidar.
As the second phase of the NGP is being implemented, the activity must be better monitored to ensure that some of the key goals are achieved. In this regard, regular reporting of the project's success and shortcomings should be made available near real-time and more accurately assessed through the combined use of in situ and satellite data. For an effective restoration of the forest, the NGP should also be accompanied by a conscious effort to reduce, if not eliminate, deforestation. The preservation of existing forests and especially old-growth forest is more desired because of their critical role in preserving the biodiversity and ecology of the region.  main movers and supporters of the National Greening Program. With much passion, he insisted on accountability, transparency and proper application of science for the monitoring of the NGP.

Conflicts of Interest:
The authors declare no conflict of interest.