Access the full text.
Sign up today, get DeepDyve free for 14 days.
References for this paper are not available at this time. We will be adding them shortly, thank you for your patience.
Article Landscape Fragmentation, Ecosystem Services and Local Knowledge in the Baroro River Watershed, Northern Philippines 1, 2 2 2 Mark Anthony M. Ramirez *, Juan M. Pulhin , Josephine E. Garcia , Maricel A. Tapia , 2 2 3 4 Florencia B. Pulhin , Rex Victor O. Cruz , Catherine C. De Luna and Makoto Inoue Resources, Environment and Economics Center for Studies, Inc. (REECS), Quezon City 1109, Philippines College of Forestry and Natural Resources, University of the Philippines Los Baños, College, Laguna 4031, Philippines; email@example.com (J.M.P.); firstname.lastname@example.org (J.E.G.); email@example.com (M.A.T.); firstname.lastname@example.org (F.B.P.); email@example.com (R.V.O.C.) Interdisciplinary Studies Center for Integrated Natural Resources and Environment Management, University of the Philippines Los Baños, College, Laguna 4031, Philippines; firstname.lastname@example.org Faculty of Human Sciences, Waseda University, Tokorozawa City 359-1192, Japan; email@example.com * Correspondence: firstname.lastname@example.org; Tel.: +63-2-995-0556 Received: 30 June 2019; Accepted: 27 September 2019; Published: 1 October 2019 Abstract: Landscape fragmentation, the breaking up of land use type into smaller parcels, is damaging watersheds worldwide. Without addressing its causes, landscape fragmentation can permanently destroy habitats and compromise ecosystem services (ES) that a watershed provides. This paper aims to establish associations between watershed landscape fragmentation and ES by integrating science (satellite imageries and fragmentation analyses) and local geographic knowledge (key informant interviews and focus group discussions) at different time periods. Using the case of the Baroro River Watershed in Northern Philippines, this paper posits that local knowledge, when integrated with scientific knowledge, becomes a significant medium through which watershed landscape fragmentation and declining quality of ES can be better understood and addressed. Results also indicate that people’s experiences and knowledge on ES coincide with watershed landscape fragmentation as evidenced by satellite images and fragmentation analyses done at different time periods. This implies that people’s knowledge is well grounded on facts and complements scientific knowledge necessary in crafting more effective landscape policies that can tackle watershed fragmentation. Study results are also crucial in providing information to serve as inputs in the development of a more robust watershed management plan; particularly in implementing sustainable land uses without sacrificing the watershed’s overall integrity. Keywords: landscape fragmentation; watershed; ecosystem services; local knowledge; Philippines 1. Introduction Tropical forests are known as the richest terrestrial ecosystems in the world in terms of flora and fauna as well as structural complexity [1,2]. These forests can be found inside watersheds, where both surface and groundwater emanate and drain into a common outlet. Watersheds are landscapes composed of a mix of local ecosystems and land uses covering a wide area. Currently, most watershed landscapes are being threatened by fragmentation, mainly caused by anthropogenic activities. As species are becoming extinct due to habitat destruction, the numerous ecosystem services (ES) being provided by these watersheds are being compromised at an alarming rate, which adversely affects both socio-ecological systems and human well-being. Ecosystems are indispensable Resources 2019, 8, 164; doi:10.3390/resources8040164 www.mdpi.com/journal/resources Resources 2019, 8, 164 2 of 29 sources of various goods and services needed by humans to survive, and keep natural processes intact [3,4]. To sustain the flow of ES, there is a need to address watershed fragmentation at the landscape level. Landscape fragmentation is defined as the breaking up of contiguous land use into several patches over time which affects the effective functioning of the ecosystem [5,6]. As deforestation happens, patches of original vegetation remain, modifying both biotic and abiotic elements of the environment . This process has been aggravated by continued expansion and encroachment of agriculture, cattle ranching, and illegal logging, converting the last remnants of tropical watersheds into isolated patches [8–11]. Studies estimate that the tropics has around 50 million forest fragments and that when summed up together will produce nearly 50 million kilometers of forest edges . Watershed landscape fragmentation is found to have substantially increased carbon emissions from tropical forests by 31%, and that Asia accounts for the second largest emission due to its topography and small-scale nature of fragments . Within the tropics are watershed landscapes composed of diverse and mosaic natural and human-modified ecosystems. As the components of the watershed landscape are interconnected, a change in the flow of ES upstream will have direct or indirect impacts to midstream and downstream areas . A change in land use due to human activities usually results to fragmentation of natural habitat that affects the capacity of the watershed to provide certain ES . Thus, understanding historical land use patterns can explain current and future conditions of the ES in a given watershed . While there has been an increase in the number of research studies investigating the link between fragmentation and ES in recent years , current studies that examine the role of historical land use in shaping the current condition of the watershed are limited in the ES literature . Historical land use has an important role in influencing the current structure and function of the environment that affects its capacity to produce ES [14,15] and has a considerable effect in the way local people view and manage their resources. Historical land use can characterize the relationship between populations and the ecosystem that they inhabit in a given time period [16,17], and can determine the driving forces of land cover change . Thus, understanding the drivers of land cover change in the past can therefore explain the differences in land use practices of people, which can be established by using local and scientific knowledge. This paper therefore assesses the link between watershed landscape fragmentation and ES by comparing scientific data and local knowledge of people interacting with their environment. Examining the case of the Baroro River Watershed in Northern Philippines, the paper argues that local knowledge is a powerful medium through which watershed landscape fragmentation and the declining quality of ES can be better explained and addressed. Currently, majority of the 143 watersheds supporting the national irrigation system in the Philippines are affected by declining forest cover . Changes in the landscape are due to ballooning population, unsustainable land use practices, and lack of coordinated governance efforts resulting to loss of biodiversity and continuing poverty. Such conditions prompted the conduct of this timely and relevant study. As knowledge of local people differ based on their cultural background, geographical location, life experiences and land practices, this study selected communities to represent the entire cross- section of the Baroro River Watershed (upstream, midstream, and downstream). This method of analysis aims to capture differences in levels of local knowledge and appreciation of ES. Aside from the standard historical land use/land cover (LULC) change analysis, this study is the first scientific attempt to quantify the extent of watershed landscape fragmentation using established landscape ecology metrics at different time periods. Official LULC maps were used, except for the 1940s as satellite images were not yet available at that time. In addition, this is the only fragmentation study done in the Philippines, both at national or local levels, which relates to ES. Results of this study are expected to contribute to the buildup of knowledge relevant in designing future plans and interventions on watershed and habitat management. Resources 2019, 8, 164 3 of 29 2. Landscape Fragmentation in the Philippines: A Historical Perspective In the 1600s, Philippine forests constituted around 90% of the total land area of the country . At that time, lands were traditionally owned and managed by the indigenous people who practiced swidden farming for survival. Under the Spanish regime in the 1880s, trees were cut and used in the construction of ships for the galleon trade. This practice reduced the forest cover to around 70% at the beginning of the 20th century . Table 1 shows data on land use changes in the Philippines from 1880 to 1980. Table 1. Land use change in the Philippines in a century (1880–1980). All Totals in 1000 km 1880 1920 1950 1970 1980 2010 Total land area 300 300 300 300 300 300 Total forest area 210 193 163 138 103 77* Total net cultivated area 20 35 47 64 97 121** Population per square kilometer 20 36 68 122 160 308*** % forest area 70% 64% 54% 46% 34% 26% % net cultivated area 7% 12% 16% 21% 32% 40% Sources: Goldewijk, et al., 2004 , *Mangobay.com, 2016 , **World Bank, 2106  ***PSA, 2012 . During the Post-World War II era, Philippine forests were given importance due to the value and contribution of its timber to economic development . Large-scale logging peaked in the 1970s that resulted in high deforestation rate, said to be among the fastest in the world . With the opening up of a significant portion of the forests in the watershed areas, migrants started to come in and occupied the logged-over areas. Migration in the country was mainly driven by landlessness and lack of employment opportunities in the lowland. As a result, more mid- and upland watershed areas were planted with agricultural crops such as rice, leading to increased fragmentation . Up till now, people are living in the uplands and contribute to changes in the watershed landscape. Given its history of large-scale timber extraction, the forest cover in the Philippines is now at about 26% of its original state . This is a clear indication that forests and watersheds are the main casualties of economic development and population growth throughout the centuries. What remains of the natural tree vegetation in the country is confined in the upper elevation of the watershed landscape, with scattered distribution in the mid- and lowlands. This historical perspective on land use change from the early 17th century to the present is significant in the study of watershed landscape fragmentation in the Philippines as it explains why the conversion of watersheds into different land uses led to the breaking up of watershed landscape into numerous patches. 3. Watershed Landscape Fragmentation and Ecosystem Services Changes in the watershed landscape are brought about by natural disturbances and human activities. As the world population grows over time, supporting people’s needs is critical for survival. This eventually led to the opening up of tropical forests for extractive industries, food production, and settlements resulting to watershed landscape fragmentation. A study of Gibbs et al.  revealed that more than half of the new agricultural lands in the tropics encroached on intact forests, while an additional 28% are estimated to open up secondary forests . Hence, fragmentation and land-use intensification are integrally intertwined . Combined with climate change, landscape fragmentation may have deleterious consequences to a wide range of ES . It weakens the health of the ecosystem by increasing the edge to interior ratio and reducing the core of the forests to the point of endangering the survival of species. A recent study estimates that more than 70% of the world’s forests are now within one kilometer of an edge usually in close proximity to human modified areas . The disruption caused by fragmentation can have devastating effects to flora and fauna that cannot adjust to the level of sunlight, rain, noise, and pollution at the edge of these fragments . Some species that require large contiguous area for habitat can be permanently lost as fragmentation reduces the ability of these species to move across landscapes, interrupting its foraging and breeding Resources 2019, 8, 164 4 of 29 activities . Large predators are usually absent in fragmented landscape that can result in food- web cascades. In this case, species at the lower level of the food chain multiply out of control, contributing to pests and diseases . Thus, only those species that can thrive in the edges of the fragments can adapt to the harsh conditions, eventually replacing the original species composition . Watershed landscape fragmentation will also likely lead to local warming due to creation of non- vegetative surfaces that absorb more solar radiation, affecting the microclimate of the watershed and exacerbating the impacts of global warming . This increases evapotranspiration that alters the local water cycle, with rivers and streams becoming more prone to either floods or droughts at any month of the year . This also affects soil dwelling organisms sensitive to increasing temperature and are critical in ecosystem processes such as nutrient cycling . As there is less vegetation to hold the soil during peak flows in a fragmented landscape, soil deposition or erosion is evident when observing the quality of water in the rivers and streams. Soil erosion encompasses water pollution and siltation, loss of organic matter, and reduction of water storage capacity of the watershed . Particularly interesting in ES studies and fragmentation is the time lags between land-use change and ecosystem responses. These time lags delay the effects of fragmentation on the provision of ES supply, which means that some effects of fragmentation can be manifested not at present but can be felt in the future . Such delayed effects are more pronounced for biodiversity-based services as re-colonization and establishing connectivity are slow processes  as well as to hydrological services where groundwater turnover times are typically very slow, thereby reversing water quality degradation is difficult . 4. Methods 4.1. Study Area The Baroro River Watershed is located in the province of La Union (Northern Luzon, Philippines, see Figure 1). It is composed of the municipalities of San Gabriel, Bagulin, San Juan, San Fernando, Bacnotan, and Santol. The watershed has a total land area of 19,063 hectares, which is the main source of water for both irrigation and domestic purposes in all municipalities except for Bagulin and Santol. The Baroro River Watershed represents the case of majority of watersheds in the Philippines in terms of anthropogenic pressures and changes in land use—from a contiguous forest in the early 1900s to numerous patches of mini-forests concentrated at the upper portion of the watershed at present. Agricultural areas, production forests, and protection forests constitute the present land uses in the watershed. The study further focused on three villages, which are representatives of the cross-section of the watershed namely, 1) Barangay Lon-oy in San Gabriel for the upstream, 2) Barangay Cabaroan in San Juan for the midstream, and 3) Barangay Baroro in Bacnotan for the downstream area. Forest vegetation can still be seen in the upper and middle portions of the watershed while its downstream portion is dominated by an urban sprawl. San Juan is one of the municipalities at the receiving end of the watershed. It is also a popular tourist destination primarily because of its famous beaches for surfing. Some communities in this side of the watershed are into fishing and farming, but agricultural lands are now slowly being converted into subdivisions and business establishments due to improving economic opportunities and outcomes. Agriculture dominates the middle portion of the watershed, particularly of rice and corn varieties, interspersed with permanent and semi- permanent dwellings. Communities that inhabit the upper watershed are composed of Kankanaeys, an indigenous people and original settlers in the area. The Kankanaeys are still practicing rice terracing, particularly in sloping areas, but commercial agriculture as compared to swidden farming is the dominant practice for the rest of the watershed. The middle and lower watersheds are occupied by migrants of Ilocano and Tagalog descents. From 55,697 individuals recorded in 1990, the population in the watershed increased to 77,166 in 2015, with an annual growth rate of 2.76 % and is expected to put more pressures on the remaining natural resources. Resources 2019, 8, 164 5 of 29 Figure 1. Location and elevation maps of the Baroro River Watershed, La Union, Philippines. 4.2. Research Design The study was conducted using a two-stage research design. First, satellite images were sourced from the official land use and land cover (LULC) maps of the National Mapping and Resource Information Authority (NAMRIA). This government facility is responsible for providing mapmaking Resources 2019, 8, 164 6 of 29 services and related information to the general public. The LULC maps represent the different land use classes of Baroro River Watershed, and correspond to the years 1988, 2003, 2010, and 2015. These four time series were the same time when NAMRIA produced the official land use cover maps of the country. The maps were then cross-referenced using available sheets from Google Earth Pro (Open Geospatial Consortium, Inc., Wayland, MA, USA). The processing and analysis of reclassification of land types were done using an ArcGIS software (ESRI, New York Street, Redlands, CA, USA). Selected landscape metrics were computed using FRAGSTATS 4.2.1 due to its accuracy in calculating fragmentation metrics and straightforwardness on raster datasets [38,39]. FRAGSTATS is a spatial pattern analysis program that offers a comprehensive choice of landscape metrics and have been used to quantify landscape structure. It is used to analyze fragmentation and describe the characteristics and components of landscapes . These statistics facilitate the comparison of landscapes and the evaluation of processes. The advantage of FRAGSTATS is that the calculations are applied in a fully integrated approach in a GIS platform, which can generate a patch map. The raster datasets were then processed in ArcGIS software. For the second stage of the research design, site visits to the three selected villages were done to validate the maps and obtain additional relevant data. Land use and environmental changes in the villages were traced based on recollection of oldest inhabitants in the areas through focus group discussions (FGDs). Prior to the FGDs, key informant interviews (KIIs) with village officials were done to gather crucial information needed to craft questions for the FGDs. Information focused on reconstruction of the changes that happened at different time periods (1988, 2003, 2010, and 2015). Specifically, the timeline presented to the participants were divided into four decades to coincide with the produced maps: 1940s (pre-mapping period), 1980s (coinciding with 1988), 2000s (coinciding with 2003), and 2010s (coinciding with 2010 and 2015). This was also done to confirm whether the two types of data have overlaps. Key questions were formulated to gain a more comprehensive understanding of how local people used watershed resources, how these practices affect watershed landscape fragmentation, the ES delivery of the Baroro River Watershed, and how people in the watershed adapt to existing conditions and projected scenarios. Also, before the start of the FGDs, an introduction about the study and its objectives, definition of watershed landscape fragmentation, and presentation of fragmentation maps based on LULC were given to the participants. In terms of ES, emphasis was placed on freshwater production, soil productivity, preservation of biodiversity, provision of food, fiber and raw materials, cultural services, and maintenance of microclimate. These were the ES that the FGD participants were most familiar with and were defined in the Millennium Ecosystem Assessment (MA). An average of eight carefully selected participants joined the FGDs, with ages that range from 35 to 83 years old representing the different age groups and gender and knowledgeable of the history of their barangays. The FGDs were conducted following the protocols set by the Overseas Development Institute (ODI) to ensure that reliable information will be generated . The older participants can recall what transpired in the 1940s through oral history and life experience, which became the study’s base year. After the FRAGSTAT results, KIIs, and FGDs, interpretations regarding the impacts of watershed landscape fragmentation to ES were made. 4.3. Satellite Imagery Classification and Accuracy Assessment Geospatial and remote-sensing data are reliable sources for understanding and determining the drivers of LULC changes at any landscape . In this study, satellite images used for LULC maps from NAMRIA that represent land use classes of Baroro River Watershed from 1988 to 2015 were obtained. The land cover maps were generated based on a classification scheme consistent with international standards for global reporting and integration . For 1988 land cover data, Landsat TM data from 1982 to 1985 were processed and supported by ground measurement conducted by the then Bureau of Forest Development (BFD) with the assistance from German Technical Cooperation Agency (GTZ), and Food and Agriculture Organization (FAO). The 2010 land cover map was sourced from Landsat 7 ETM+ data from 2010 with 30-m resolution and classified based on FAO-FRA Project Field Inventory Manual with NAMRIA and Forest Management Bureau (FMB) collaboration; ALOS Resources 2019, 8, 164 7 of 29 AVNIR-2 (10 m resolution), SPOT 5 (10 m resolution) and Landsat 7 ETM+ (30 m resolution). Data were validated on the ground with accuracy assessment (calculated by province through overlaying the ground validated sample data over the preliminary map and presented the success or failure of the matches in Confusion matrix with overall accuracy of 89%). Landsat operational land imager (OLI) data from 2015 was used, with 30 m resolution with reference data from Google Earth Pro, topographic maps, ground truth data and Interferometric Synthetic Aperture Radar (IFSAR) data for coastline using OBIA eCognition software for digital classification with ground validation and accuracy assessment of 93.92% . The land cover maps were validated and presented to the Department of Environmental and Natural Resources (DENR) regional and provincial offices, with local government units (LGUs), and selected stakeholders for comments, suggestions, and approval. The 1988, 2003, 2010, and 2015 LULC maps are the officially recognized maps by the Philippine Government that can be used by the national government offices, academe and other organizations/institutions. 4.4. Image Reclassification A reclassification scheme of five classes (from 12 aggregated categories) was developed based on physiographical knowledge of the study area, supporting ancillary data, researchers’ prior local knowledge, and visual interpretation using historical function of Google Earth Pro. The five LULC classes were categorized as water bodies, grassland/barren land, mixed agriculture, forest, and built- up areas (Table 2). Table 2. Reclassified land use/land cover classes for Baroro River Watershed based on NAMRIA LULC data. LULC Class Description Water River and reservoirs Grassland/Barren land 70% grassland, open areas Mixed agriculture Annual crop, perennial crop, cereals and sugar, cropland, cultivated land with brush land, wooded land/wooded grassland, forest plantation, shrubs, open Forest forest/broadleaved Built-up area Area with residential, infrastructure development 4.5. Landscape Metrics Spatial metrics were analyzed at both class and landscape levels. This means that each patch type in the landscape mosaic of the watershed was analyzed as a whole. At the class level, eight landscape metrics were selected including class area, number of patches, edge density, total edge, median patch size, area-weighted mean shape index, mean shape index, and mean perimeter-area ratio. At the landscape level, nine indices were selected, namely: 1) number of patches, 2) mean patch size, 3) edge density, 4) area-weighted means shape index, 5) perimeter-area fractal dimension, 6) contagion, 7) aggregation index, 8) largest patch index, and 9) mean shape index. These landscape metrics were chosen because it can show changes in land use activity, as demonstrated in past studies [43,44] and as described in Table 3. Table 3. Landscape metrics used in the Baroro River Watershed. Metric Name Levels of Acronym Description Reference (Units) Application [45,46] CA Class area (ha) Total area of all patches per class Class level TE Total Edge (m) Total edge is an absolute measure of total edge Class level  length for a particular patch type (class level) or for all patch types (landscape level) [45–47] NP Number of patches Total number of patches Class and landscape level Resources 2019, 8, 164 8 of 29 Metric Name Levels of Acronym Description Reference (Units) Application MPS Mean patch size Average patch size Class and  (ha) landscape level ED Edge density Total length of edge of a certain LULC class per Class level [46,49] (m/ha) unit area (m/ha) AWMSI Area-weighted It measures the complexity of patch shape of a Class and [46,47] Mean Shape Index particular LULC class covered to a standard landscape shape (square), by weight patches according to level their size. It equals 1 when all patches are square and increases with complexity of patch shapes MSI Mean shape index MSI is equal to 1 when all patches are circular Class level [45,50] (for polygons) or square for raster (grids) and it increases with increasing patch shape irregularity PAFRAC Perimeter- area It indicates the relationship between the area Landscape [45,51] Fractal Dimension and perimeter of the patches level CONTAG Contagion (%) CONTAG is widely used in landscape ecology, Landscape [48,52,53] it measures both patch type interspersion (e.g., level the intermixing of units of different patch types) as well as patch dispersion (e.g., the spatial distribution of a patch type) at the landscape level. It indicates the aggregation of patches. AI Aggregation Index Calculates the like adjacencies of different pairs Landscape [47,54] of patch types involving the focal class level LPI Largest Patch Percentage of the landscape comprised by the Landscape  Index (%) largest patch of the corresponding LULC class level MSI Mean Shape Index MSI equals the sum of the patch perimeter (m) Landscape [45,51] divided by the square root of patch area (m2) for level each patch of the corresponding patch type, adjusted by a constant to adjust for a circular standard (vector) or square standard (raster), divided by the number of patches of the same type; in other words, MSI equals the average shape index (SHAPE) of patches of the corresponding patch type 5. Results 5.1. Land Use and Land Cover Change Dynamics The LULC classification from the five classes across four time periods indicate the conversion of forest and grassland/barren lands to mixed agriculture and built-up areas (Table 4 and Figure 2). The dominant land use/land cover that increased progressively over the study period was mixed agricultural land and built-up areas. For example, cultivated land increased by 14% from 2003 to 2010, while built-up areas increased by 100% in the same period. On the contrary, forest cover decreased by 8% from 2003 to 2010 and 5% from 2010 to 2015. Grassland/barren lands were also reduced by 94% from 1988 to 2015. In addition, grassland/barren land and forest lands were mainly converted to cultivated, built-up areas. Resources 2019, 8, 164 9 of 29 Table 4. LULC changes in the Baroro River Watershed from 1988 to 2015. 1988 2003 2010 2015 Cover Change between Periods (%) LULC Class Area (ha) Shares Area (ha) Shares Area (ha) Shares Area (ha) Shares (1988–2003) (2003–2010) (2010–2015) (1988–2015) Grassland/Barren land 117.93 0.01 0.00 0.00 3.32 0.00 7.27 0.00 −100 +100 +54 −94 Built-up area 0.00 0.00 0.00 0.00 279.26 0.01 502.03 0.03 0 +100 +44 +100 Forest area 13,596.88 0.71 14,137.35 0.74 13,034.50 0.68 12,454.41 0.65 +4 −8 −5 −8 Inland water 166.27 0.01 166.27 0.01 166.27 0.01 166.27 0.01 0 0 0 0 Mixed agriculture 5302.04 0.28 4879.50 0.25 5699.77 0.30 6053.14 0.32 −8 +14 +6 +14 Resources 2019, 8, 164 10 of 29 Figure 2. LULC maps for 1988, 2003, 2010 and 2015 in Baroro River Watershed, Northern Philippines (Source: NAMRIA). Land use change detection is an important process in monitoring and managing land resources because it provides quantitative analysis of spatial variation of land use pattern . As shown in Table 4, there has been a considerable change in land use during the 17-year period. In the upper part of Baroro River Watershed, forest was the dominant land use. However, results show that dense forest cover decreased by 8% during the study period, while mixed agriculture increased by 14%. Forest areas have been converted to cultivated and built-up areas with a remarkable increase of 100% (Figure 3). Note that the official LULC 1998 and 2003 maps from NAMRIA did not include a category of built up areas. However, this particular category was included by NAMRIA in the 2010 and 2015 maps. Secondary documents and local interviews revealed that human activities such as increasing population and proliferation of farms in the study area occurred from 1998 to 2003. From 1998 to 2003, forest cover registered a significant increase, which government and other sectors attribute to regrowth vegetation, plantations established through reforestation projects, and spontaneous tree growing activities on public and private lands . Correspondingly, the area for mixed agriculture also decreased as a direct effect of these efforts. Resources 2019, 8, 164 11 of 29 Figure 3. Cultivated and built-up areas in the upper, middle and lower part of Baroro Watershed, La Union (Source: Google Earth Pro). 5.2. Landscape Fragmentation: Spatial Metrics at Landscape Level Landscape metrics are used to determine structure and quantify changes in the configuration of a landscape . Also, landscape metrics are discrete and static quantitative indicators of landscape . Spatial metrics applied at landscape level were used in this study to measure the graininess of the landscape (e.g., the tendency of the landscape to exhibit fine versus coarse-grain texture). Metrics in Table 5 help determine the structure and fragmentation of land use/land cover of the Baroro River Watershed landscape, and all mosaic structures and related information were obtained in conjunction with GIS and FRAGSTATS. The spatial metrics at landscape level showed significant changes between four time periods (Table 5 and Figure 4). The number of patches increased from 101 in 1988 to 786 in 2015, and the mean patch size declined from 189.93 to 24.41, indicating that the landscape is dominated by small patches. This is supported by the area weighted mean shape index (AWMSI), where the value increased from 1988 to 2015, indicating an augmentation of the non-regularity of patches in the watershed. The ED increased significantly from 2003 to 2010 and may be a result of the rapid increase of NP during this period. The PAFRAC value continued to increase during the period and approached the critical status of 1.5, suggesting a decrease in the stability of the land use/land cover structure of the Baroro River Watershed landscape. The CONTAG analysis showed a series of decreasing values. The results of landscape metrics show that several small forests were removed between 1988 and 2015 resulting in reduced patches and sizes. The agricultural expansion in Baroro River Watershed and conversion of forestland into different land uses during the last two decades led to fragmentation of contiguous forests affecting the watershed landscape. Resources 2019, 8, 164 12 of 29 Table 5. Spatial metrics change at landscape level in Baroro River Watershed. Between four time periods. Metrics 1988 2003 2010 2015 Number of patches (NP) 101 120 384 786 Mean patch size (MPS; ha) 189.93 159.86 49.96 24.41 Edge density (ED; (m/ha) 11.58 15.34 50.68 75.96 Area-weighted Mean Shape Index (AWMSI) 3.46 4.95 6.96 9.85 Perimeter-area Fractal Dimension (PARFAC) 1.31 1.33 1.39 1.45 Contagion (CONTAG) 74.72 70.44 63.41 61.36 Aggregation Index (AI) 98.22 97.70 92.42 88.64 Largest Patch Index (LPI) 46.90 64.95 17.82 23.44 Mean Shape Index (MSI) 1.85 1.89 1.97 1.87 5.3. Landscape Fragmentation: Spatial Metrics at Class Level Changes in landscape metrics at class level between four time periods are presented in Table 6. There have been changes in the size, number, and spatial distribution of fragments, with different patterns for different land uses at class level. Landscape fragmentation in the upper portion of the Baroro River Watershed increased and encroached in the natural forest cover through agricultural expansion and deforestation. The forest cover is highly dynamic in patch type. Forest area is the patch type with the largest class area (CA). The total landscape area decreased from 71% (13,595.49 ha) to 65% (12,453.48 ha), while the patch number increased from 9 to 214. The mean patch size (MPS) also drastically decreased from 1510 in 1988 to 55 in 2015 for the dense forest. The increase in number of forest patches as well as decrease in mean patch size indicate that forest cover in the watershed has been fragmented. Several forest patches have also been converted into numerous small patches and isolated in recent years. AWMSI showed changes for forest type class, indicating an increase in the shape complexity and fragmentation. Table 6. Spatial metrics change at the class level in the Baroro River Watershed between four time periods. Year CA (ha) NP ED (m/ha) TE (m) MPS (ha) AWMSI MSI MPAR Grassland/Barren 117.99 1 0.31 6000 117.99 1.37 1.37 50.85 Built-up area 0.00 0 0.00 0 0.00 0.00 0.00 0.00 Forest 13,595.49 9 8.04 154,230 1510.61 2.79 1.91 260.71 Mixed agriculture 5305.23 10 9.99 191,730 452.50 3.88 2.10 500.04 Water 164.25 81 4.82 92,520 2.03 4.41 1.82 997.61 Grassland/Barren 0.00 0 0.00 0 0.00 0.00 0.00 0.00 Built-up area 0.00 0 0.00 0 0.00 0.00 0.00 0.00 Forest 14,138.01 16 13.14 252,060 1139.11 3.34 2.09 509.53 Mixed agriculture 4880.70 23 12.71 243,840 158.68 3.43 2.05 412.96 Water 164.25 81 4.82 92,520 2.03 4.41 1.82 997.61 Grassland/Barren 3.51 3 0.07 1,320 1.17 1.28 1.27 569.72 Built-up area 280.62 24 3.17 60,840 11.69 2.82 1.89 370.03 Forest 13,031.10 170 61.17 1,173,390 52.45 4.81 1.89 2290.48 Mixed agriculture 5703.48 106 32.13 616,290 28.62 3.98 1.94 970.68 Water 164.25 81 4.82 92,520 2.03 4.41 1.82 997.61 Grassland/Barren 7.29 4 0.13 2580 1.82 1.52 1.44 545.82 Built-up area 500.31 192 9.36 179,640 2.61 2.38 1.51 663.23 Forest 12453.48 214 78.54 1,506,750 55.34 5.26 2.00 463.17 Mixed agriculture 6059.34 295 59.02 1,132,260 12.09 5.69 1.76 456.43 Water 164.25 81 4.86 93,270 2.03 4.41 1.82 997.61 Resources 2019, 8, 164 13 of 29 Figure 4. Patch maps of Baroro River Watershed in 1988, 2003, 2010, and 2015. Resources 2019, 8, 164 14 of 29 Agricultural land is one of the important land use classes in the watershed and plays an important role in fragmentation of the natural forest. During the study period, agricultural land increased from 5305.23 ha (28%) in 1988 to 6059.34 ha (32%) in 2015. The number of patches (NP), edge density (ED), total edge (TE), and area-weighted mean shape index (AWMSI) increased progressively over time with the exception of the mean patch size, which showed the reverse pattern. This suggests that mixed agriculture patches have become more complex and been broken into small sizes. 5.4. Ecosystem Services at Spatial and Temporal Scales Examining the land use history of the Baroro River Watershed requires going back to past events that led to its current condition. These events may explain why watershed landscape fragmentation occurred and what ES were affected by changes in the land use and forest cover. To do this, three villages were selected to represent the cross-section of the watershed and accurately describe the connectivity of each section in the context of a watershed landscape. Barangay Lon-Oy depicts past and current situation in the upstream watershed section, Barangay Cabaroan for the midstream section, and Barangay Baroro for the downstream section. A timeline was established and divided according to decades corresponding to fragmentation maps. These were then used to gather temporal knowledge from the local people about the history and uses of the Baroro River Watershed. 5.4.1. Upstream Section or Headwater in Barangay Lon-Oy, San Gabriel The earliest inhabitant arrived in Barangay Lon-Oy from the neighboring province of Benguet in the early 1920s to search for areas suitable for kaingin or swidden farming. Stories passed orally to the next generation of Kankanaeys describe the watershed area as being covered with lush green forest, with around 90% forest cover until the 1940s. Inside the watershed were river and streams with clear water used by Kankanaeys for drinking and cooking. The FGD participants shared that the clean and steady flow of water came from the thick forest where wildlife of numerous species still existed, with some being hunted for food such as wild boars (Sus philippensis) and wild deer (Rusa marianna); see Table 7). Different kinds of fish abounded in the rivers such as eel (Anguilla marmorata), crabs (Brachyura spp.) and ipon or largesnout goby (Awaous melanocephalus), an endemic species in the Baroro River Watershed. Houses were made from bamboo (Bambusa vulgaris) and other grasses sourced from the forest and nearby areas, and there was no need to cut trees for timber. Settlers used to survive by eating root crops such as taro (Dioscorea alata) and cassava (Manihot esculenta) harvested from kaingin farms, and rice from paddies. The soil back then was described as very productive since there were limited number of farms for few families and crops were only for household consumption. Kankanaeys used caves as burial grounds based on their belief system. The microclimate was generally cool throughout the year. The forest cover started to decline in the 1980s, from the 90% in the 1940s to an estimate of 60% in the 1980s. This was attributed to the increasing number of families living inside the village. Water became polluted at this time, particularly due to use of chemical fertilizers and pesticides up to a point that it was not safe for drinking. Some farmers began to use a 60/40 ratio of organic to chemical fertilizer as the soil was not that productive. Still, food was sufficient and rice terraces were a familiar sight in the area. As population increased, so did the houses constructed using wood as primary material. Timber became for sale and trees from the forest such as narra (Pterocarpus indicus), red lauan (Shorea negrosensis), and yakal (Shorea astylosa) were cut moderately. There was an increase in hunting activities to support the growing number of families, thus the number of wildlife species began to go down. Caves were replaced by cemeteries as permanent burial place. The weather became hotter and drought was experienced by the people for the first time. In the 2000s, people in Lon-Oy observed that the forest cover in the watershed started to increase, but this was not due to a natural occurrence. The government led reforestation activities using exotic species such as gmelina (Gmelina arborea), narra (Pterocarpus indicus), and mahogany (Swietenia macrophylla). Despite this effort, the volume of water flowing from the river was declining, especially during summer. The ratio of organic to chemical fertilizer used by farmers was recorded at 50/50 Resources 2019, 8, 164 15 of 29 based on the estimates of the FGD respondents. Village officials also prohibited the cutting of dipterocarps during this period and encouraged families to plant fruit trees such as marang (Artocarpus odoratissimus), jackfruit (Artocarpus heterophyllus), guyabano (Anona muricate), and rambutan (Nephelium lappaceum). The number of wild animals waned as wild boar (Sus philippensis) and deer (Rusa marianna) were not sighted anymore. According to the respondents, the microclimate was hotter preventing them from using their beddings at night. Typhoons became stronger with increased gustiness that damaged roads and the mini-hydro plant in the village. The situation worsened during the decade following 2010. While forest cover was described by local people as improving due to additional reforestation efforts by the government and the private sector such as Holcim, a cement manufacturing company in the Philippines, water volume continued to decline and its quality worsened as algal bloom was observed in some portions of the river. Diet of the people was completely replaced by livestock as boar and other wild animals cannot be found in the forest. Weather continued to be erratic and temperature soared. 5.4.2. Midstream Section in Barangay Cabaroan, San Juan The same as with Lon-Oy, the forest cover in Barangay Cabaroan was estimated to be at 90–95% in the 1940s (see Table 8). Water in the river was used for drinking and irrigation purposes; while different species of fish thrived due to high quality and volume. Rice farming was still limited and that the soil was described as very productive. If needed, farmers used carabao manure as fertilizer. Diet was composed of fish, vegetables, and meat—some coming from animals in the wild such as frogs (Fejervarya limnocharis) and wild chicken (Gallus gallus domesticus). Houses were made of bamboos and only branches of trees were cut for fuel wood. Several varieties of timber and fruit trees abound in the area such as duhat (Syzygium cumini), guava (Psidium guajava), mabolo (Diospyros blancoi), and guyabano (Anona muricate) to name a few. People used to have a picnic near the river and freshwater can easily be dug from the soil. Climate was generally colder during this period. Reduction in the forest cover from 95% to 70% in the 1980s was observed by the participants. The river was described as in good condition and can be used for drinking, though some portions were already contaminated from the use of chemicals upstream, which coincides with what happened in Barangay Lon-Oy during this period. Farmers started to build deep wells at a depth of 10 feet. As agriculture intensified, so did the use of chemical fertilizers, especially by wealthy farmers. It was also during the later years of this period when farmers were forced to use chemical pesticides as stem borers and grasshoppers attacked crops. People were increasingly relying on the market for food as commercial livestock were available and accessible. Decreasing number of fish and absence of wild animals for food explained the shift to other food sources in the village. While some benefits were not derived from the watershed, people still go to the river for recreation. This was also the period when the climate began to become hotter. The next period saw a steep decline in forest or tree cover, as recalled by the respondents. Data showed a decline to 30% in 2000s from the previous 70% in the 1980s due to the increasing number of farms devoted to commercial agriculture. As all farmers were completely relying on chemical fertilizers and pesticides, water quality in the river further declined. In 2009, people started to buy bottled water for drinking. The period also saw an increase in output of charcoal making to support the tourism industry in the coastal areas of San Juan, resulting to more cutting of trees. There were lesser fish in the rivers and native trees were replaced by fruit trees. Picnic was not viable along the river as the quality of the surrounding areas became degraded. People started to use air-conditioning system to get them through summer. There was no change in the forest/tree cover in 2010 and the succeeding years from the previous period. However, 90% of the population relied on bottled water for drinking, costing the household more money to survive. Some farmers began shifting to organic agriculture, however, majority of the farmers were not convinced that the soil can support crop production naturally. Demand for charcoal still exists but authorities regulated its sale to prevent further cutting of trees. This decade also saw extreme typhoons passing the area with stronger wind and rains such as typhoons Ondoy, Pepeng, Yolanda, and Egay. Resources 2019, 8, 164 16 of 29 5.4.3. Downstream Section in Barangay Baroro, Bacnotan The present physical characteristic of Barangay Baroro at the downstream section of the watershed is completely different from its situation in the 1940s (see Table 9). Eighty percent (80%) of the area was said to be covered with trees as only 50 families were living in the village back then. Water quality in the river was clean and the volume was steady; making it ideal for different fish species to thrive. Open-dug wells were already existent during this period, but had shallow depths as compared now. While farmers used chemical fertilizers for rice (Oryza sativa), corn (Zea mays), tobacco (Nicotiana tabacum L.), and other crops, this was still at the 50/50 ratio. With the village near the ocean, presence of sea grass (Enhalus spp), seashells, and healthy coral reef abounded as households depended on the ocean for sustenance. Houses were made of natural materials such as bamboo (Babusa vulgaris) and pawid (Nypa fruticans). There were also several varieties of fruit-bearing trees as well as wild animals for hunting. Ipon (Awaous melanocephalus) was abundant during September to October. Respondents described the temperature to be cold and starts in the afternoon. In the 1980s, tree cover started to decline due to increasing number of houses for 200 to 300 migrants coming from different parts of Ilocos Region. Water in the river was described as clean, but quality started to decline. The river’s reach was getting wider and the depth became shallower. On the other hand, deep wells depths were recorded at 20 feet, indicating over abstraction of groundwater. Farmers were heavily relying on chemical fertilizers. As soil productivity worsened, the number of farmers in the village dropped to 14. Food was sourced from farms and market and not from the watershed. Enterprising individuals began putting fish pens in the river. The number of wildlife species for both flora and fauna, including fish, continued to decrease at an alarming rate. As a sign of declining supply, the price of ipon increased during this period. The weather started to become hotter. Storm surge was experienced in 1987. Except for some patches of trees, forest cover was almost non-existent by 2000s. Houses increased to around 400 to 500, resulting in the establishment of subdivisions. The depth of deep wells remained at 20 feet, but there was absence of groundwater during summer. This convinced the local government to put up a communal potable water supply. Drinking water was sourced through refilling stations. River reach became wider, which was particularly dangerous to people living along the riverside especially during typhoon season due to flooding. With respect to occupation, some farmers shifted to other jobs in nearby municipalities. This period also saw an increase in the number of fish pens to 16. Unregulated harvesting resulted to decreasing number of fish species in the river and ocean, particular ipon. There were limited sightings of green viper (Trimeresurus flavomaculatus), monitor lizard (Varanus indicus), dolphin (Orcaella brevirostris), and tortoise (Eretmochelys imbricata). Beginning in 2010, salt water intrusion was experienced due to overabstraction of water. This was brought about by pressure from a burgeoning population, which increased to 500-600 households composed mostly of foreigners. Due to this, water refilling stations proliferated. The water quality in the river worsened and became brownish in color. People were experiencing extreme weather events. A fish kill was recorded in 2017. 6. Discussion Results of this study show the connection between local knowledge and land use metrics. The increasing number of patches (NP), decreasing mean patch size (MPS) as well as increasing edge density (ED) at the landscape level corresponded to the declining quality of ES. This includes loss of biodiversity, declining soil productivity, and deteriorating water supply as experienced by people living inside the watershed over four time periods being studied. Conversion of previously vegetated land use to agriculture and built up areas continues up to the present, which contributed to the near critical status of the watershed; specifically approaching the PARFAC of 1.5 and stood at 1.45 in 2015. This implies that the stability of the land use/land cover of the Baroro River Watershed is continuously declining. Resources 2019, 8, 164 17 of 29 Table 7. Status of ecosystem services through time in Barangay Lon-Oy, San Gabriel, La Union (Upstream/Headwaters). Ecosystem 1940s 1980s 2000s 2010s Services • Increasing forest cover due to reforestation through the National Greening Program (NGP) and other private • 60% forest cover • Increasing forest cover due to initiatives (i.e., Holcim) • Increased number of kaingin reforestation activities of the • Through NGP, fruit trees were • 90% forest cover farms government and individual efforts planted inside the forest using Freshwater • Water was crystal clean • Water flowing from the to plant gmelina (Gmelina arborea), coffee (Coffea arabica), rambutan production • Strong but steady flow in rivers and river was not potable due to Mahogany (Swietenia macrophylla), (Nephelium lappaceum), guyabano streams pollution and narra (Pterocarpus indicus) (Anona muricate), and lanzones • Safe and potable • Cases of poisoned animals • Number of farms still increasing (Lansium domesticum) from drinking water in the • Declining water quantity flowing • Decreasing volume of water in river from the river; the river • Decreasing water quality (i.e., brownish in color) • Presence of algae in the river • Naturally growing trees for house construction, i.e., narra (Pterocarpus indicus), red lauan (Shorea • Timber harvesting was prohibited • People hunt and forage for food negrosensis), and yakal • Naturally growing trees was • Houses were made of concrete, • Bamboo (Bambusa vulgaris) and giant (Shorea astylosa) confined to upper slopes hence, limited use of timber/trees ferns (Diplazium esculentum) abound Food, fiber, • Harvesting of timber for sale • Pioneer tree species such as marang • Increase number of farmers into and used for house construction and raw • Depleting food sources (i.e., (Artocarpus odoratissimus), langka livestock raising • Root crops as food sources such as ube materials wild animals) forcing (Artocarpus heterophyllus), guyabano • Bamboo (Bambusa vulgaris) and (Dioscorea alata), carrot (Daucus carota), increased hunting activities (Anona muricat), and rambutan rattan (Calamus manan canes) and cassava (Manihot esculenta) • Declining number of (Nephelium lappaceum) occupied the started to decline in numbers • Timber harvesting was limited bamboos, giant ferns, and forest timber • Food remained sufficient Resources 2019, 8, 164 18 of 29 Ecosystem 1940s 1980s 2000s 2010s Services • Many species of fish were abundant in • Wild deer (Rusa marianna) rivers, i.e., ipon (Awaous melanocephalus was already non-existent • Practice of planting exotic tree juvenile stage), bunog (Awaous • There were monitor lizards species • No more sightings of deer (Rusa melanocephalus adult stage), eel (Anguilla (Varanus indicus) sighted • No more sightings of deer (Rusa marianna), wild boar (Sus marmorata), lobster (Nephropidae spp.), • Exotic species of trees not yet marianna), wild boar (Sus philippensis) in the forest Maintenance and crab (Brachyura spp.) planted in the forest philippensis) in the forest suggesting that they are now of biodiversity • Wild boars (Sus philippensis), deer (Rusa • Declining number of birds • Number of monitor lizard (Varanus extinct marianna), monitor lizard (Varanus were observed indicus) and birds started to decline • Continuous practice of planting indicus), monkey (Macaca fascicularis), • Harvest of ipon (Awaous • Increased number of sighting of exotic species and fruit trees and wild cat (Felis silvestris) to name a melanocephalus juvenile stage) heron (Nycticorax nycticorax) few were abundant in the forest with started to decline presence of a wild range of birds • Cemetery replaced caves as burial ground • Caves used by local people as burial • Increased number of farmers Cultural sites practicing terracing of rice • Increasing number of farmers • Farmers started to go back to tabas services • With traditional system of naming and other crops practicing terracing system of composting typhoons • Flowering of bamboo as signal that there is an impending drought • Temperature even during dawn started to get hot; (people were not using beddings and jackets at • Microclimate remained cool even in • Temperature started to rise night) Microclimate • Erratic weather is the new summertime • Case of moderate drought • Experienced strong typhoons and normal • Intensity of rain was not heavy observed every 10 years winds as compared in the past (i.e., Typhoon Ondoy, Pepeng, Perya) that damaged the community’s mini-hydro plants and roads Resources 2019, 8, 164 19 of 29 Table 8. Status of ecosystem services through time in Barangay Cabaroan, San Juan, La Union (Midstream). Ecosystem 1940s 1980s 2000s 2010s Services • 30% forest/tree cover • Some forests were converted into houses and agricultural farms • Water from the river was not • 70% forest cover potable • 30% forest cover (no • Water in the river was already • Frequent flooding change from previous • 90-95% forest cover contaminated due to brownish • Increased number of decade) Freshwater • Water from the river used for drinking cover, which people attributed households in the area • Increasing number of production • Clear and high water level in the river (average of to upstream activities • Some households people in the village, 15 feet in depth) • Water was potable constructed septic tanks to 90% of which drink • Plenty of fish in the river and its tributaries • Construction of deep-wells for prevent river pollution mineral water irrigation started (depth at 10 • Application of chlorine in • Low surface water level feet) deep wells to make the water potable; • In 2009, people started to buy drinking water; • Deep wells remained at a depth of 20 feet • Farmers who had enough capital already used chemical fertilizers as the soil cannot support • Started to use chemical production fertilizer in the 1990s • Some farmers practiced • Pests were a natural occurrence • All farmers in the barangay organic farming due to • Soil was productive without the need for inputs such as stem borer (Scirpophaga already used chemical Soil high price of produce in • Carabao manure as fertilizer incertulas), grasshopper (Caelifera fertilizers productivity the market • Limited number of kaingin spp) • Three cropping cycles for • Continuous use of • Two-cropping seasons rice rice (Oryza sativa) and corn chemical fertilizers (Oryza sativa) and corn (Zea (Zea mays) can be supported mays), peanut (Arachis hypogaea); by chemical fertilizers due to chemical inputs • Increasing number of farm lands Resources 2019, 8, 164 20 of 29 Ecosystem 1940s 1980s 2000s 2010s Services • People were dependent on fish and wild animals for food • No cutting of trees as bamboos (Bambusa vulgaris) and grasses were used for house construction Food, fiber, • Food from wild animals and fish • Cutting of trees for charcoal • Charcoal making • Vegetables were plenty such as sitaw (Vigna and raw declined making regulated by the local unguiculata ssp. sesquipedalis), talong (Solanum materials • Forest vegetation thick with • No more wild animals and government melongena), upo (Lagenaria siceraria), and kangkong bamboos (Bambusa vulgaris) fish • No more wild animals (Ipomoea aquatica) and fish • Presence of wild frogs (Fejervarya limnocharis), snakes (Boiga cynodont), and chicken (Gallus gallus domesticus); • Fish and other water species were abundant, i.e., dalag (Channa striata), hito (Clarias batrachus), carpa (Cyprinus carpio), ipon (Awaous melanocephalus • Declining number of fish in the juvenile stage), eel (Anguilla marmorata), tilapia river such as hito (Clarias (Oreochromis niloticus), sugpo (Penaeus monodon), batrachus), eel (Anguilla talangka (Synchiropus spp.), bunog (Awaous • Much lesser in number and marmorata), and ipon (Awaous melanocephalus adult stage), kampa (Rhiyacichthys kinds of fish caught in the melanocephalus juvenile stage) aspro), shells (bukkukao), and suso (Lymnaeidae and river (composition and • No more carps in the river Ampullariidae spp.) distribution of • Wild animals such as wild boar • Wild fauna: carabao (Bubalus bubalis carabanesis), macroinvertebrate Maintenance (Sus philippensis), deer (Rusa sheep (Ovis aries), chicken (Gallus domesticus), wild communities • Same condition during of biodiversity marianna), and monkeys (Macaca pig (Sus philippensis), monkey (Macaca fascicularis), • Same situation in for the previous decade fascicularis) were already non- usa (Rusa marianna), bayawak (Varanus indicus), vegetation, specifically forest existent bird, and snake (Boiga cynodon) trees: narra (Pterocarpus • Monitor lizard (Varanus indicus) • Horses were used for transportation via karitela indicus) and fruit bearing became lesser in numbers • Fruit trees abound: guava (Psidium guajava), duhat trees such as guava (Psidium • Mabolo (Diospyros blancoi), (Syzygium cumini), narra (Pterocarpus indicus), guajava) kaimito (Chrysophyllum caimito), molawin (Vitex parviflora), mabolo (Diospyros and atis (Annona squamosa) were blancoi), sineguelas (Spondias purpurea), guyabano lesser in number (Anona muricate), anonas (Annona reticulata), atis (Annona squamosa), kamachile (Pithecellobium dulce), and aratilis (Muntingia calabura) Resources 2019, 8, 164 21 of 29 Ecosystem 1940s 1980s 2000s 2010s Services • People regularly schedule Cultural picnics near the river; • People used to go to the river for family • Lesser number of people • Recreation by the river services • People used “tambuli” or gatherings/picnics doing picnics by the river was not an option due carabao horn to warn people of • They just dig the soil to get freshwater due to pollution to pollution danger such as floods • Extreme typhoons were • Much hotter climate than experienced, with Microclimate before stronger winds and • Fresh air was colder due to abundant trees • Weather was hot • Some houses started to use more intense rains (i.e., • Flooding was observed near the Barangay Hall due • Flooding was experienced in air conditioning system Typhoon Ondoy, to extraordinary intense rain 1976 • Experienced severe drought Pepeng, Yolanda, in 2009. Lando) • Hotter climate Table 9. Status of ecosystem services through time in Barangay Baroro, Bacnotan, La Union (Downstream). Ecosystem Services 1940s 1980s 2000s 2010s • Households increasing to about 400 to 500 • Foreigners started to • Households now at • Increasing number of houses due to live in the area, drawn 500-600 Freshwater the establishment of subdivisions • 80% forest/tree cover by its beaches • Saltwater intrusion production • Increased number of households to • Community was established in the 1920s • Communal Baroro has been observed around 200-300; composed of mostly • Original 50 households in the area Potable Water Supply • Increasing number migrants • Water in the river and streams was crystal clean was established to of water refilling • Deep wells were 20 feet deep to • Presence of open-dug wells provide safe drinking stations abstract water; • River reach was narrow water to households • Water from the river • River reach became wider and depth • Deep wells were 20 feet becomes brown became shallower deep to abstract water; during rainy season but no water during summer Resources 2019, 8, 164 22 of 29 Ecosystem Services 1940s 1980s 2000s 2010s • River reach continued to become wider and depth became shallower • Limited number of farmers who planted corn (Zea Soil • Heavy reliance on chemical fertilizers mays), tobacco (Nicotiana tabacum L.), rice (Oryza productivity (100%) • Farmers shifted to other sativa), and peanut (Arachis spp) • Decreasing number of farmers due to jobs as harvests • Some wealthy farmers using chemical fertilizers in land conversion continued to decline the 1960s • People became dependent on livestock. • Fish were caught farther • Presence of fish pens with milkfish than before. • Fish were plenty in rivers and ocean, i.e., ipon (Chanos chanos) and tilapia • Fish pens increased in Food, fiber and (Awaous melanocephalus juvenile stage), tilapia (Oreochromis niloticus) number raw materials (Oreochromis niloticus), and talaba (Crassostrea iredalei) • Decrease in the number of fishes in the • Same as the previous • Sea grass started to • Presence of native spinach (Talinum fruticosum) river decade. disappear • Availability of bamboo (Bambusa vulgaris) and pawid • Houses were made of concrete • Fish species same as the (Nypa fruticans) for house construction • Ipon (Awaous melanocephalus juvenile previous decade but stage) present in the river numbers declined • Sea grass can be seen in coastal areas • Healthy coral reefs • Declining number of wild • Dolphin (Orcaella • Presence of sea grass in coastal areas animal/marine species, i.e., monitor brevirostris) and tortoise • Decline in fish • Fora: lanete (Wrightia pubescens), anahaw (Saribus lizard (Varanus indicus), wild cat (Felis (Eretmochelys imbricate) harvest and number rotundifolius) camachile (Pithecellobium dulce), niyog silvestris), snakes (Boiga cynodon), disappeared of wildlife Maintenance (Cocos nucifera), acacia (Samanea saman), balete (Ficus monkeys (Macaca fascicularis), dolphin • Declining number of • Local government- of biodiversity benjamina), sampalok (Tamarindus indica), duhat (Orcaella brevirostris), and tortoise monitor lizard (Varanus initiated program to (Syzygium cumini), mushroom, madre cacao (Gliricidia (Eretmochelys imbricate) indicus), and green viper bring back the sepium), ipil (Intsia bijuga), yakal (Shorea astylosa), and • Declining catch from the ocean, i.e., (Trimeresurus pawikan narra (Pterocarpus indicus), octopus (Octupos vulgaris) flavomaculatus) (Eretmochelys • Fish species: hipon (Fenneropenaeus chinensis), ipon • Stingray (Taeniura lymma) became low • Declining number of imbricate) (Awaous melanocephalus juvenile stage), bongo in number fish catch Resources 2019, 8, 164 23 of 29 Ecosystem Services 1940s 1980s 2000s 2010s (Awaous melanocephalus adult stage), padaw (Perna viridis), Prawn (Penaeus semisulcatus), aguas (Planiliza macrolepis), umok (shell), stringray (Taeniura lymma), imkadit (coral fish), octopus (Octupos vulgaris), dugong (Dugong dugon), and dolphin (Orcaella brevirostris) • Wild fauna: labuyo (Gallus domesticus), wild cat (Felis silvestris), bayawak (Varanus indicus), sawa (Phyton reticulatus), frog (Fejervarya limnocharis), and monkeys (Macaca fascicularis) • Abundance of ipon from September to December Cultural • Floating restaurant services was established in • None identified • None identified • None identified Baroro River Watershed Microclimate • Increasing number • Local climate started to get warmer • Strong typhoons and of storm surges • Cold weather • Storm surge was experienced in 1987 flooding in 1992 • 2017 fish kill • Stronger typhoons Resources 2019, 8, 164 24 of 29 In addition, the number of agricultural patches increased and the mean patch size decreased over time. This indicates that more farms have encroached in the watershed and affected the general flow of ES. Changes in the landscape metrics over time can also be explained by the economic growth within the watershed, particularly in the middle and lower portions as it became urbanized. The relationship between local knowledge and scientific data was distinctly manifested in the upstream section of the Baroro River Watershed. The communities within are more proximate to the forest and the people benefit from the various ES that the watershed provides on a daily basis. On the other hand. people from the middle and lower portions of the watershed benefit from specific ES such as water for irrigation and soil fertility for crop production; as these are the most visible services valued [2,58]. They do not recognize other ES such as maintenance of biodiversity and regulation of microclimate unlike those living in the upstream section. This is due to the fact that the people in both midstream and downstream areas have more opportunities at their disposal. They can easily get food from the market unlike those living in the upstream who rely on forests and wildlife for sustenance. Thus, knowledge about a certain ES also depends on the location, access to resources, and socioeconomic status . Historically, the conversion of forests to agricultural farms far outweighed the benefits that the watershed can provide since the value of ES coming from the watershed are less obvious and did not readily translate into monetary terms. This became the primary driver of land use conversion in the Baroro River Watershed. Moreover, the loss of forest or tree cover in the middle section of the watershed as well as in the downstream communities coincides with the loss of knowledge about the ES from the watershed. An exception to this case are flooding and other natural disasters as people can easily relate to these adverse conditions, particularly in the upper watershed. When people experience the benefits or impacts of certain ES, they tend to put more value to a particular ES. Tapping into this inherent knowledge may translate local knowledge into actions necessary to rehabilitate the Baroro River Watershed. Particularly, interventions to reconnect the different landscape fragments to regain what was lost in terms of species richness and other indicators of watershed health can be planned. As majority of the present generation has no recollection of what transpired in the past in terms of forest cover and landscape, which led to the current condition of the watershed, providing orientation or guidance on the benefits that one can derive from the watershed is of primary importance. In this manner, the communities can learn the importance of nature and value of the ES that it provides. Nevertheless, the landscape fragmentation in the Baroro River Watershed did not only result from the breaking down of the watershed into patches, but also from a confluence of other factors such as fragmented governance, piecemeal policies, and unsustainable land management in the past. Lack of coordination, particularly among LGUs and community members on appropriate land uses in the watershed, has led to its current state. There was also lack of communication between national government agencies responsible for forestry and agriculture on how to manage the watershed landscape due to its varying and often overlapping mandates; as well as policies that were dysfunctional and not appropriate on the ground. This was evidenced by the increasing area of agriculture over time, particularly on inorganic farming and resulting to encroachment on forestlands. Also, as majority of the population in the watershed are migrants, land practices of indigenous groups are slowly becoming obsolete. In addition, local knowledge about the ES in the watershed was evident from the results of the study’s FGDs, KIIs, and fragmentation analyses. Due to increased awareness based on the study’s results, there was a collective effort from the three LGUs to halt landscape fragmentation by initiating action to declare the Baroro River Watershed as a protected area. At this point, reconnection of fragments will start at the institutional level. Lastly, this study showed that local knowledge on the history of land use in the watershed and the use of land use metrics provide complementary information in understanding the causes of fragmentation and its impacts to certain ES. Local people identified not only changes in the environment but also the time frame when these changes happened and its effect on livelihood. Land use metrics are therefore important inputs that can be used by decision makers to justify initiatives Resources 2019, 8, 164 25 of 29 or plans, while local knowledge on ES can provide entry points for various initiatives on both watershed landscape rehabilitation and biodiversity conservation. 7. Conclusions The current research is a preliminary attempt to establish a link between watershed landscape fragmentation and ES by integrating science and local knowledge as methods of analysis. Results provide scientific basis that the increasing fragmentation of the Baroro River Watershed highly compromised the condition of certain ES due to loss of ecosystem functions. This adds to the understanding of the relationship between watershed landscape fragmentation and its effects on the delivery of ES in the watershed, particularly in the tropics where similar studies are scant. From 1940s to 2015, the forests of the Baroro River Watershed were lost to the expansion of agriculture and built up areas, and mostly evident in the middle and lower sections of the watershed based on land use metrics. Local knowledge from people living within the watershed confirmed this finding as they experienced the adverse impacts of land use conversion and watershed fragmentation through declining water supply, loss of biodiversity, incidence of pests and diseases, cases of floods and droughts, and worsening micro-climate. Reversing these conditions, however, requires more detailed studies on ecosystem functions as time delays govern the changes in ES . Given the increasing number of fragments in the watershed landscape, it is seemingly a daunting task to attain a contiguous and intact forest of 50,000 hectares in the near future; which is said to be the minimum requirement to ensure that watershed landscape and ecosystems are properly functioning [59,60]. For one, this will entail the rehabilitation of the entire Baroro River Watershed and adjacent watersheds to satisfy the area requirement. Understanding the dynamics of small patches is therefore of critical importance as it provides ES such as water quality regulation and carbon storage. At the watershed level, “reconnecting” these fragments would allow for an increase in area where biodiversity can thrive and re-populate, creating biodiversity corridors. While local governments, farmers, and the private sector have already initiated reforestation initiatives in Baroro River Watershed, these efforts should be deliberately targeted to connect fragments and reduce edge effects. One possible solution is to have an integrated watershed landscape management plan anchored on the management of these fragments and with primary focus on sustainable land uses that address both the environment and social requirements towards a healthy watershed. To this end, the study proponents suggest the use of landscape metrics and historical land use analysis as framework for watershed management planning to help adequately explain the historical land use pattern and anticipate future impacts brought about by time lags in ES supply. Fragmentation does not only happen at the biophysical level but also at the governance, policy, and management levels. It is not only a result of increasing population but rather of other factors such as uncoordinated watershed management between local governments and national government agencies, unsustainable land use policies and lack of awareness of stakeholders vis-à-vis the situation of the Baroro River Watershed that happened in the past. Similar with other watersheds in the country, the governance in the study area started when problems on increasing fragmentation and deforestation manifested itself through decreasing flow of ES readily felt by people. While actions had been taken by the national government, particularly the local DENR office and the communities upstream to plant trees, efforts to manage the development of midstream and downstream areas were neglected due to lack of coordinated planning. In most cases and not just in the Philippines but also in other tropical countries, forests are the main casualties of development. The challenge now is how to avert the further fragmentation of the watershed to ensure that ES will remain intact for future generations. Results from this study posit that reconnection can happen by first addressing these institutional issues in the watershed. This is beginning to take shape in the Baroro River Watershed as different actors have realized that past actions contributed to the fragmentation of the watershed, and that they can manage their natural resources to arrest the declining levels of ES. Resources 2019, 8, 164 26 of 29 In addition, this proposed framework would also create social and institutional corridors that will tap existing resources coming from the different stakeholders in the watershed landscape. In this context, local knowledge on watershed resources and ES can be used as entry points for the process as well as in effectively crafting incentives and disincentives through policies at the local level. Determining the current appreciation of the various stakeholders of the watershed functions would help in crafting effective education campaign, especially among the younger generation, especially those who are seemingly detached from the natural environment. Author Contributions: All authors participated in the collection, processing and analysis of data for the study as well as contributed in writing the article. Funding: This research was funded by the Department of Agriculture–Bureau of Agricultural Research (DA- BAR) as part of the project entitled, “Community-Based Adaptation: Enhancing Climate Change Resilience of Communities and Ecosystems through Participatory Watershed Management”. Acknowledgments: We would like to thank the Provincial Government of La Union and the Municipal Governments of San Gabriel, San Juan and Bacnotan, La Union, Philippines for allowing us to conduct the study in Baroro Watershed and in committing to continue the research initiative in the area. We would also like to extend our gratitude to the Provincial Environment and Natural Resources Office–Department of Environment and Natural Resources of Region 1 for helping us to access their data and maps used in this article. Conflicts of Interest: The authors declare no conflict of interest. References 1. Whitmore, T.C. An Introduction to Tropical Rain Forests; Clarendon Press: Oxford, UK, 1990. 2. Silvano, R.; Udvardy, S.; Ceroni, M.; Farley, J. An Ecological Integrity Assessment of a Brazilian Atlantic Forest Watershed Based on Surveys of Stream Health and Local Farmers’ Perceptions: Implications for Management. Ecol. Econ. 2005, 53, 369–385. 3. Renó, V.; Novo, E.; Escada, M. Forest Fragmentation in the Lower Amazon Floodplain: Implications for Biodiversity and Ecosystem Service Provision to Riverine Populations. Remote Sens. 2016, 8, 886. 4. Mitchell, M.; Bennett, E.; Gonzalez, A. Strong and Nonlinear Effects of Fragmentation on Ecosystem Service Provision at Multiple Scales. Environ. Res. Lett. 2015, 10, 094014. 5. Davidson, C. Issues in Measuring Landscape Fragmentation. Wildl. Soc. Bull. 1998, 26, 32–37. 6. Maryland-National Capital Park and Planning Commission (MNCPPC). Forest Fragmentation in Prince George’s County: Measuring Forest Cores and Edges to Determine Fragmentation Trends; MNCPPC: Riverdale, MD, USA, 2016. 7. Bélanger, L.; Grenier, M. Agriculture Intensification and Forest Fragmentation in the St. Lawrence Valley, Quebec, Canada. Landsc Ecol. 2002, 17, 495–507. 8. FAO (Food and Agriculture Organization). Manual on Deforestation, Degradation, and Fragmentation Using Remote Sensing and GIS; FAO: Rome, Italy, 2007. 9. Gibbs, H.K.; Ruesch, A.S.; Achard, F.; Clayton, M.K.; Holmgren, P.; Ramankutty, N.; Foley, J.A. Tropical Forests were the Primary Sources of New Agricultural Land in the 1980s and 1990s. Proc. Natl. Acad. Sci. USA 2010, 107, 16732–16737. 10. Da Ponte, E.; Kuenzer, C.; Parker, A.; Rodas, O.; Oppelt, N.; Fleckenstein, M. Forest Cover Loss in Paraguay and Perception of Ecosystem Services: A Case Study of the Upper Parana Forest. Ecosyst. Serv. 2017, 24, 200–212. 11. Sunderland, T.; Abdoula, R.; Ahammad, R.; Asaha, S.; Baudron, F.; Deakin, E.; Duriaux, J.; Eddy, I.; Foli, S.; Gumbo, D.; et al. A Methodological Approach for Assessing Cross-Site Landscape Change: Understanding Socio-Ecological Systems. For. Policy Econ. 2017, 84, 83–91. 12. Brinck, K.; Fischer, R.; Groeneveld, J.; Lehmann, S.; De Paula, M.D.; Pu ̈tz, S.; Sexton, J.O.; Song, D.; Huth, A. High-Resolution Analysis of Tropical Forest Fragmentation and its Impact on the Global Carbon Cycle. Nat. Commun. 2017, 8, 14855. 13. Azevedo, J.C.; Pinto, M.A.; Perera, A.H. Forest Landscape Ecology and Global Change: An Introduction. In Forest Landscapes and Global Change; Springer: New York, NY, USA, 2014. 14. Ziter, C.; Graves, R.; Turner, M. How do land-use legacies affect ecosystem services in United States cultural landscapes? Landsc. Ecol. 2017, 32, 2205, doi:10.1007/s10980-017-0545-4. Resources 2019, 8, 164 27 of 29 15. Bürgi, M.; Östlund, L.; Mladenoff, D. Legacy Effects of Human Land Use: Ecosystems as Time-lagged Systems. Ecosystems 2017, 20, 94–103. 16. Castillo, A.; Magaña, A.; Pujadas, A.; Martínez, L.; and Godínez, C. Understanding the interaction of rural people with ecosystems: A case study in a tropical dry forest of Mexico. Ecosystems 2005, 8, 630643. 17. Caballero-Serrano, V.; Alday, J.; Amigo, J.; Caballero, D.; Carrasco, J.C.; McLaren, B.; Onaindia, M. Social Perceptions of Biodiversity and Ecosystem Services in the Ecuadorian Amazon. Hum Ecol. 2017, 45, 475– 486, doi:10.1007/s10745-017-9921-6. 18. Villafuerte, L.R. House Bill 1795. In Proceedings of the House of Representatives, Seventeenth Congress, Quezon City, Philippines, 20 July 2016. 19. Verburg, P.H.; Veldkamp, A. Projecting land use transitions at forest fringes in the Philippines at two spatial scales. Landsc. Ecol. 2004, 19, 77–98. 20. International Partnership for the Satoyama Initiative (IPSI). Use and Management of “Muyong” in Ifugao Province, Northern Luzon Island in the Philippines. 2016. Available online: http://satoyama- initiative.org/en/use-and-management-of-muyong-in-ifugao-province-northern-luzon-island/ (accessed on 15 January 2016). 21. Goldewijk, K.K.; Ramakutty, N. Land Use Changes During the Past 300 Years. In Natural Resources Policy and Management; Verheye, W., Ed.; Encyclopedia of Life Support Systems (EOLSS), Eolss Publishers: Oxford, UK, 2004. Available online: http://www.eolss.net (accessed on 6 February 2016). 22. Mangobay.com. Philippine Forest and Information Data. 2016. Available online: http://rainforests.mongabay.com/deforestation/2000/Philippines.htm (accessed on 6 February 2016). 23. World Bank. Agricultural Land Area of the Philippines. 2016. Available online: http://www.tradingeconomics.com/philippines/agricultural-land-percent-of-land-area-wb-data.html (accessed on 6 February 2016). 24. PSA (Philippine Statistics Authority). 2010 Philippine Population Density; PSA: Quezon City, Philippines, 25. Borlagdan, S.B.; Guiang, E.S.; Pulhin, J.M. Preliminary Assessment of Community-Based Forest Management in the Philippines; Institute of Philippine Culture, Ateneo de Manila University, 2001. Available online: http://dlc.dlib.indiana.edu/dlc/bitstream/handle/10535/7541/CBFM%20Preliminary%20Assessment.pdf?se quence=1 (accessed on 12 June 2016). 26. Espaldon, V.O.; Smit, B. Community reforestation in the Philippines: An Evaluation of Community Contracts. Knowl. Policy Int. J. Knowl. Transf. Util. 1998, 10, 34–42. 27. FAO (Food and Agriculture Organization). State of the World’s Forests 2011; FAO: Rome, Italy, 2011. 28. Didham, R.K. The ecological consequences of habitat fragmentation. Encycl. Life Sci. A21904 2010, doi:10.1002/9780470015902.a0021904. 29. Leadley, P.; Pereira, H.M.; Alkemade, R.; Fernandez-Manjarrés, J.F.; Proença, V.; Scharlemann, J.P.W.; Walpole, M.J. Biodiversity scenarios: Projections of 21st century change in biodiversity and associated ecosystem services. In Secretariat of the Convention on Biological Diversity, Montreal. Technical Series No. 50; Secretariat of the Convention on Biological Diversity: Montreal, QC, Canada, 2010. 30. Haddad, N.M.; Brudvig, L.A.; Clobert, J.; Davies, K.F.; Gonzalez, A.; Holt, R.D.; Lovejoy, T.E.; Sexton, J.O.; Austin, M.P.; Collins, C.D.; et al. Habitat Fragmentation and its Lasting Impact on Earth’s Ecosystems. Sci. Adv. 2015, 1, e15000052. 31. Fischer, J.; Lindenmayer, D.B. Landscape modification and habitat fragmentation: A synthesis. Glob. Ecol. Biogeogr. 2007, 16, 265–280. 32. Mitchell, M.; Suarez Castro, A.; Martinez-Harms, M.; Maron, M.; McAlpine, C.; Gaston, K.; Johansen, K.; Rhodes, J. Reframing Landscape Fragmentation’s Effects on Ecosystem Services; Opinion Cell Press: Cambridge, MA, USA, 2015. 33. Ellis-Cockcroft, I.; Cotter, J. Tropical Forest Fragmentation; Implications for Ecosystem Function. In Greenpeace Research Laboratories Technical Report (Review) 02-2014; Greenpeace: London, UK, 2014; p. 15. 34. Frelich, E. Forest and Terrestrial Ecosystem Impacts of Mining; The University of Minnesota Center Forest Ecology: Minneapolis, MN, USA, 2014. 35. Mitchell, M.; Bennett, E.; Gonzalez, A. Forest Fragments Modulate the Provision of Multiple Ecosystem Services. J. Appl. Ecol. 2014, 51, 909–918, doi:10.1111/1365-2664.12241. Resources 2019, 8, 164 28 of 29 36. Decocq, G.; Andrieu, E.; Brunet, J.; Chabrerie, O.; De Frenne, P.; De Smedt, P.; Deconchat, M.; Diekmann, M.; Ehrmann, S.; Giffard, B.; et al. Ecosystem Services from Small Forest Patches in Agricultural Landscapes. Curr. For. Rep. 2016, 2, 30–44, doi 10.1007/s40725-016-0028-x. 37. Aguirre-Gutie’rrez, J.; Biesmeijer, J.C.; van Loon, E.E.; Reemer, M.; WallisDeVries, M.F.; Carvalheiro, L.G. Susceptibility of Pollinators to ongoing Landscape Changes Depends on Landscape History. Divers. Distrib. 2015, 21, 1129–1140. 38. McGarigal, K.; Cushman, S.A.; Ene, E. FRAGSTATS v42.1: Spatial Pattern Analysis Program for Categorical and Continuous Maps, Computer Software Program Produced by the Authors at the University of Massachusetts, Amherst. 2012. Available online: http://www.umass.edu/ landeco/research/fragstats/fragstats.html (accessed on 20 July 2018). 39. MacLean, M.G.; Congalton, R.G.A. Comparison of Landscape Fragmentation Analysis Programs for Identifying Possible Invasive Plant Species Locations in Forest Edge. Landsc. Ecol. 2015, 30, 1241–1256. 40. Omar, D. Focus Group Discussion in Built Environment Qualitative Research Practice. IOP Conf. Ser.: Earth Environ. Sci. 2018, 117, 012050, doi:10.1088/1755-1315/117/1/012050. 41. Hansen, M.C.; DeFries, R.S.; Townshend, J.R.G.; Sohlberg, R. Global Land Cover Classification at 1 km Spatial Resolution Using a Classification Tree Approach. Int. J. Remote. Sens. 2000, 21, 1331–1364. 42. Santos, R.N. National Mapping Efforts: The Philippines Land Cover/Land Use Changes (LU/LUC)) and its Impacts on Environment in South/Southeast Asia; International Regional Science Meeting: Quezon City, Philippines, 43. Kamusoko, C.; Aniya, M. Land Use/Cover Change and Landscape Fragmentation Analysis in the Bindura District, Zimbabwe. Land Degrad. Dev. 2007, 18, 221–233. 44. Tolessa, T.; Senbeta, F.; Kidane, M. Landscape Composition and Configuration in the Central Highlands of Ethiopia. Ecol. Evol. 2016, 6, 7409–7421, doi:10.1002/ece3.2477. 45. McGarigal, K.; Marks, B.J. Fragstats: Spatial Pattern Analysis Program for Quantifying Landscape Structure. Reference Manual. For. Sci. Dep. Or. State Univ. Corvallis Or. 1995, 59, 1–122. 46. Karatassiou, M.; Galidaki, G.N.; Sklavou, P. Landscape Pattern Changes in Response to Transhumance Abandonment on Mountain Vermio (Nort Greece). Sustainability 2015, 7, 15652–15673. 47. Smiraglia, D.; Ceccarelli, T.; Bajocco, S.; Perini, L.; Salvati, L. Unraveling Landscape Complexity: Land Use/Land Cover Changes and Landscape Pattern Dynamics (1954–2008) in Contrasting Peri-urban and Agro-forest Regions of Northern Italy. Environ. Manag. 2015, 56, 916–932. 48. Leitão, A.B.; Miller, J.; Ahern, J.; McGarigal, K. Measuring Landscapes—A Planner’s Handbook; Island Press: Washington, DC, USA, 2006. 49. Syrbe, R.; Walz, U.; 2012. Spatial Indicators for the Assessment of Ecosystem Services: Providing, Benefiting and Connecting Areas and Landscape Metrics. Ecol. Indic. 2012, 21, 80–88, doi:10.1016/j.ecolind.2012.02.013. 50. Gökyer, E. Understanding Landscape Structure Using Landscape Metrics. In Advances in Landscape Architecture; IntechOpen: London, UK, 2013. 51. Wang, X.; Blanchet, G.F.; Koper, N. Measuring Habitat Fragmentation: An Evaluation of Landscape Pattern Metrics. Methods Ecol. Evol. 2014, 5, 634–646. 52. Turner, M.G. Landscape ecology: The effect of pattern on process. Ann. Rev. Ecol. Syst. 1989, 20, 171–197. 53. O’Neill, R.V.; Hunsaker, C.; Timmins, S.; Jackson, B.; Jones, K.; Ritters, K.; Wickham, J. Scale Problems in Reporting Landscape Pattern at the Regional Scale. Landsc. Ecol. 1996, 11, 169–180. 54. Ghosh, A.; Munshi, M.; Areendran, G.; Joshi, P.K. Pattern Space Analysis of Landscape Metrics for Detecting Changes in Forests of Himalayan Foothills. Asian J. Geoinformatics 2012, 12, 1. 55. Zhang, Z.; Gao, J. Linking Landscape Structures and Ecosystem Service Value Using Multivariate Regression Analysis: A Case Study of the Chaohu Lake Basin, China. Environ. Earth Sci. 2016, 75, 3. 56. Chokkalingam, U.; Carandang, A.P.; Pulhin, J.M.; Lasco, R.D.; Peras, R.J.J.; Toma, T. One Century of Forest Rehabilitation in the Philippines: Approaches, Outcomes and Lessons. In Country Case Studies on Review of Forest Rehabilitation Initiatives: Lessons from the Past; Center for International Forestry Research (CIFOR): Situ Gede, Sindang Barang Bogor Barat, Indonesia, 2016. 57. McGarigal, K.; Cushman, S.A. Comparative Evaluation of Experimental Approaches to the Study of Habitat Fragmentation Effects. Ecol. Appl. 2002, 12, 335–345. 58. Lewan, L.; Soderqvist, T. Knowledge and Recognition of Ecosystem Services Among the General Public in a Drainage Basin in Scania, Southern Sweden. Ecol. Econ. 2002, 42, 459–467. Resources 2019, 8, 164 29 of 29 59. Gibson, L.; Lee, T.M.; Koh, L.P.; Brook, B.W.; Gardner, T.A.; Barlow, J.; Peres, C.A.; Bradshaw, C.J.A.; Laurance, W.F.; Lovejoy, T.E.; et al. Primary Forests are Irreplaceable for Sustaining Tropical Biodiversity. Nature 2011, 478, 378–381. 60. Laurance, W.F.; Camargo, J.L.C.; Luizão, R.C.C.; Laurance, Susan, G.; Pimm, S.L.; Bruna, E.M.; Stouffer, P.C.; Bruce Williamson, G.; Benítez-Malvido, J.; Vasconcelos; Heraldo, L. The Fate of Amazonian Forest Fragments: A 32-year Investigation. Biol. Conserv. 2011, 144, 56–67. © 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
Resources – Multidisciplinary Digital Publishing Institute
Published: Oct 1, 2019
Access the full text.
Sign up today, get DeepDyve free for 14 days.