BUILDING A BENCHMARK FOR CHANGE: Land Cover Updating in For Sustainable Future

Presentation on “Isabela Land Cover Assessment and Watershed Mapping” by Dr. J. Andres Ignacio December 07, 2020 PRESENTATION OUTLINE Structure and Key Talking Points for Discussion PRESENTATION OUTLINE • 1.) ESSC-ILCA Project Overview a.) Context of the Problem b.) Project Objectives c.) Methodologies Implemented d.) Project Scope and Data limitations

• 2.) Land Cover Updating: The Research Process a.) Data acquisition, preparation and input b.) Analysis, validation and finalization c.) Knowledge generation and method documentation d.) Knowledge sharing and collaboration

• 3.) Land Cover Assessment: Initial Results and Findings a.) Overall accuracy and validity b.) Forest and land cover statistics c.) Upland and protected areas

• 4.) Sample ILCA datasets integration: Watershed Mapping ESSC-ILCA PROJECT OVERVIEW Context of the Problem | Project Objectives | Methodologies Implemented | Project Scope and Data Limitations CONTEXT OF THE PROBLEM

Biological Diversity

NSMNP

Water Regulation Cultural Identity CONTEXT OF THE PROBLEM

• The existing spatial on the Northern Sierra Madre and the province of Isabela might be limited in terms of its timeliness and usefulness in land management initiatives.

NSMNP• More recent advances in remote sensing technology will allow to create a more holistic, realistic and robust baseline maps. PROJECT OBJECTIVES

San Pablo • ILCA’s goal is to contribute to monitoring forest and land cover change in the Sierra Madre landscape (especially to the 9 municipalities that are part of the NSMP). City • Produce an updated and detailed land cover assessment for the San Mariano province of Isabela, which encompasses the Northern Sierra Madre Natural Park. METHODOLOGIES IMPLEMENTED Recent Advances in Remote Sensing Technologies and Multi-stakeholder Engagement Approach • The classification methodology employed by ESSC, Sen2Agri, was developed by its Belgian partner Universite Catholique de Louvain and is used by the European Space Agency as its primary toolbox for land cover assessment using data from the Sentinel 2 satellite platform. Sentinel 2 allows much more frequent data collection for any point on Earth (every 5 days) and provides high accuracy and rapid land cover data.

• Specifically, it has the following strengths and advantages:

1.) High Temporal and Spatial Resolution Satellite Imagery

2.) Avoidance of fuzzy and confusing land use and/or land cover classes PROJECT SCOPE AND DATA LIMITATIONS 1.) The project focuses on creating a robust baseline and updated data for the land cover of the province. No further land cover analysis (i.e. change detection) were conducted since no existing comparable data is available. Yet, standardization is possible to make initial assumptions.

2.) Primary and secondary data used for validating identified forest formations and classifications identified is limited. Further studies on this can be explored by interested individuals or groups.

3.) Due to COVID19 restrictions, validation sessions were done remotely. Technical difficulties and recurring localized lockdowns in most of the LGUs resulted to a 73% completion rate of online consultations. PROJECT SCOPE AND DATA LIMITATIONS 4.) Data collection challenges in ground truthing activities for some of the barangays in selected municipalities/cities due to safety issues (i.e. unstable weather conditions and presence of rebel groups)

5.) Misclassification of few selected crops such as calamansi and banana identified in limited spaces, reflected by smaller pixels in the output maps.

6.) Majority of the training datasets built and gathered from the field survey were focused on accessible areas particularly near existing road networks. LAND COVER UPDATING: The Research Process Data Acquisition, Preparation, and Input | Analysis, Validation and Finalization | Knowledge Generation and Method Documentation | Knowledge Sharing and Collaboration Objective 4: Knowledge Sharing THE RESEARCH PROCESS and Collaboration Objective 3: Knowledge Generation and Method Documentation Objective 2: Analysis, Validation and Finalization Objective 1: Data Acquisition, Preparation and Input THE RESEARCH PROCESS 1.1 Inception meeting and consultation with local stakeholders and collaborators in the province

KEY ACCOMPLISHMENTS AND MILESTONES Objective 1: Consultative meetings with key local collaborators in Data Acquisition, Preparation and the province particularly Mabuwaya Foundation, ISU- Input ✓ Cabagan and CVPED Joint PAMB En Banc Resolution No. 01, series of 2019: Approving the Request of ESSC to Conduct ILCA ✓ Project in all Protected Areas in Isabela (3 DENR- PAMB’s in Isabela) Resolution No. 441, series of 2019: Signing of Memorandum of Agreement with the Provincial ✓ Government of Isabela

THE RESEARCH PROCESS 1.2 Researched (desktop and on-site) and conducted ground truthing activities to gather, acquire or purchase existing paper and digital maps, studies, reports, tables and data in other formats Objective 1: Data Acquisition, KEY ACCOMPLISHMENTS AND MILESTONES Preparation and Input 67 relevant studies and reports were acquired from local stakeholders, partners and collaborators in the ✓ field (i.e. Mabuwaya Foundation) 302 additional geospatial documents were collected from both national line agencies and local ✓ government units Collected 13,088 geotagged photos from 36 out of 37 municipalities (97%) and 914 out of 1,055 barangays ✓ (89%)

THE RESEARCH PROCESS 2.1 Employed a classification procedure which takes the most updated training and validation polygons Objective 2: into account (i.e. April 2018 – February 2020) to Analysis, Validation produce the best classification output and Finalization

KEY ACCOMPLISHMENTS AND MILESTONES 589 satellite images were pre-processed but only 350 images were further used for feature extraction and ✓ analyzed; 10 iterations conducted Around 12,208 training polygons were created and ✓ prepared prior and during the project 19 land cover features were distinguished; labels for ✓ major lakes and other identifiable rivers were added THE RESEARCH PROCESS 2.2 Improvement and enhancement of the training datasets were made through independent accuracy Objective 2: assessments and stakeholders’ validation Analysis, Validation and Finalization KEY ACCOMPLISHMENTS AND MILESTONES Online training course was designed to involve interested local collaborators in the enhancement of ✓ built training datasets, resulting to an additional 748 training polygons based on the capstone project outputs In lieu of the multi-stakeholder validation sessions, ILCA Virtual ‘Mapathon’ was held from August to ✓ November. A total of 27 out of 37 (73%) participated in individual online consultation

THE RESEARCH PROCESS Objective 3: 3.1 Created a simple interactive web maps Knowledge Generation and with uploaded geospatial datasets for easy Method Documentation access of interested local stakeholders, partners and collaborators.

KEY ACCOMPLISHMENTS AND MILESTONES Uploaded output maps and shapefiles along with metadata and process documentation, in a separate ✓ web page dedicated for ILCA datasets

THE RESEARCH PROCESS Objective 4: 4.1 Local stakeholders were given technical Knowledge Sharing training that focuses on integrating data and Collaboration outputs from the project. Also, to maximize the potential of the datasets local stakeholders and collaborators are encouraged to continually update and validate the GIS database

KEY ACCOMPLISHMENTS AND MILESTONES Maintained strategic partnerships with local stakeholders and collaborators by strengthening ✓ ‘online’ community of practice (i.e ILCA-GIS group) Initial commitments and interests were forwarded by identified local partners in using the ILCA datasets ✓ produced for further studies, analysis and verification (i.e. Mabuwaya Foundation and ISU-CVPED) LAND COVER ASSESSMENT: Initial Results and Findings Overall Accuracy and Validity | Land Cover Statistics | Upland and protected areas INITIAL RESULTS AND FINDINGS: Overall Accuracy and Validity

Overall Accuracy and Validity of ILCA Training Datasets 1

0.9 0.95 0.9 0.8 0.83 0.79 0.7

0.6 NT 0.5 1st OT 2nd OT 0.4 VM 0.3

0.2

0.1

0 ANNUAL OA ANNUAL KAPPA PERENNIAL OA PERENNIAL KAPPA NT – No training datasets 1st OT – 1st online training data enhancement Source: Isabela Land Cover Assessment 2020, ESSC-ILCA 2nd OT – 2nd online training data enhancement VM – Virtual ‘mapathon’ series INITIAL RESULTS AND FINDINGS: Overall Accuracy and Validity • Participation of local stakeholders in the validation of initial post- processed output maps increases the overall accuracy and precision of the training datasets

• Involvement of local stakeholders in such process foster transparency in the process of delivering outputs. This also adds credibility and shared ownership of the final outcomes

• Collaboration with local stakeholders and collaborators also promotes data acceptability which will facilitate the use of the final outputs of the ILCA project INITIAL RESULTS AND FINDINGS: Isabela Land Cover Statistics 2020 ISABELA LAND COVER STATISTICS Forest Areas Annual Cropland Other Land Land with Tree Cover Perennial Cropland

4% 10%

14% 41%

31%

Source: Isabela Land Cover Assessment 2019, ESSC-ILCA INITIAL RESULTS AND FINDINGS: Forests Land Area Statistics in Isabela, 2019

% Share to the Forest Land Area of Isabela

Mangrove 0.14

Ultramafic Forest 13.4

Limestone Forest 4.75

Secondary Forest 41.77

Primary Forest 28.22

Mossy/Montane Forest 11.69

0 5 10 15 20 25 30 35 40 45

Source: Isabela Land Cover Assessment 2019, ESSC-ILCA INITIAL RESULTS AND FINDINGS: Forests Land Area Statistics in Isabela, 2019

% Share to the Forest Land Area of Isabela

Mangrove 0.14599.92 (in ha.)

Ultramafic Forest 56,295.9513.4 (in ha.)

Limestone Forest 19,4.75 953.91 (in ha.)

Secondary Forest 175, 447.2241.77 (in ha.)

Primary Forest 118, 530.4928.22 (in ha.)

Mossy/Montane Forest 49, 109.2411.69 (in ha.)

0 5 10 15 20 25 30 35 40 45

Source: Isabela Land Cover Assessment 2019, ESSC-ILCA INITIAL RESULTS AND FINDINGS: Annual Cropland Area Statistics in Isabela, 2019

% Share to the Annual Cropland Area of Isabela

Sugarcane 5.49

Cassava 0.15

Rice 50.45

Corn 43.91

0 10 20 30 40 50 60

Source: Isabela Land Cover Assessment 2019, ESSC-ILCA INITIAL RESULTS AND FINDINGS: Annual Cropland Area Statistics in Isabela, 2019

% Share to the Annual Cropland Area of Isabela

5.49 Sugarcane 17, 479.92 (in ha.)

Cassava 0.15462.92 (in ha.)

Rice 160,50.45 737.83 (in ha.)

Corn 139,43.91 927.36 (in ha.)

0 10 20 30 40 50 60

Source: Isabela Land Cover Assessment 2019, ESSC-ILCA INITIAL RESULTS AND FINDINGS: Other Land Area Statistics in Isabela, 2019

% Share to the Other Land Area of Isabela

Water 12.14

Built-up Surfaces 24.44

Grassland 60.08

Bare Soil/Sediment 3.32

0 10 20 30 40 50 60 70

Source: Isabela Land Cover Assessment 2019, ESSC-ILCA INITIAL RESULTS AND FINDINGS: Other Land Area Statistics in Isabela, 2019

% Share to the Other Land Area of Isabela

Water 17,12.14 428.82 (in ha.)

Built-up Surfaces 35,24.44 090.86 (in ha.)

Grassland 86, 60.08256.94 (in ha.)

Bare Soil/Sediment 3.324,775.62 (in ha.)

0 10 20 30 40 50 60 70

Source: Isabela Land Cover Assessment 2019, ESSC-ILCA INITIAL RESULTS AND FINDINGS: Perennial Cropland Area Statistics in Isabela, 2019

% Share to the Perennial Cropland Area of Isabela

Coconut 29.17

Banana 69.24

Calamansi 1.58

0 10 20 30 40 50 60 70 80

Source: Isabela Land Cover Assessment 2019, ESSC-ILCA INITIAL RESULTS AND FINDINGS: Perennial Cropland Area Statistics in Isabela, 2019

% Share to the Perennial Cropland Area of Isabela

Coconut 11, 513.6429.17 (in ha.)

Banana 27, 327.0969.24 (in ha.)

Calamansi 1.58625.34 (in ha.)

0 10 20 30 40 50 60 70 80

Source: Isabela Land Cover Assessment 2019, ESSC-ILCA INITIAL RESULTS AND FINDINGS: Land and Forest Cover in Isabela

• 19 land cover features extracted based from the analysis

• Verified location of forest types using secondary data provided by NSMNP, Mabuwaya Foundation and LGUs

• Forest and annual cropland areas dominate the land cover of the province INITIAL RESULTS AND FINDINGS: Upland and protected areas INITIAL RESULTS AND FINDINGS: Upland and protected areas Key Discussions on the Results and Findings: Land Cover, Upland and Protected Areas in Isabela • The entire Isabela including the NSMNP landscape is in a critical state because of existing and newly emerging threats to the rich biodiversity and forest areas need to be examined for its future impacts

• The extensive reach of the NSMNP calls for an integrated approach in managing and monitoring its natural resources not only by the concerned local agencies, but also by the cities/municipalities situated at the buffer zones

• Economic and ecological trade-offs of forest usage happen as a rural livelihood mechanism among affected communities Sample ILCA Datasets Integration: Watershed Mapping Example of ILCA datasets integration and analysis SAMPLE ILCA DATASETS INTEGRATION: Watershed Mapping in Isabela End of Presentation