Forest Condition Monitoring Using Very-High-Resolution Satellite Imagery in a Remote Mountain Watershed in Nepal

Satellite imagery has proven extremely useful for repetitive timeline-based data collection, because it offers a synoptic view and enables fast processing of large quantities of data. The changes in tree crown number and land cover in a very remote watershed (area 1305 ha) in Nepal were analyzed using a QuickBird image from 2006 and an IKONOS image from 2011. A geographic object-based image analysis (GEOBIA) was carried out using the region-growing technique for tree crown detection, delineation, and change assessment, and a multiresolution technique was used for land cover mapping and change analysis. The coefficient of determination for tree crown detection and delineation was 0.97 for QuickBird and 0.99 for IKONOS, calculated using a line-intercept transect method with 10 randomly selected windows (1×1 ha). The number of tree crowns decreased from 47,121 in 2006 to 41,689 in 2011, a loss of approximately 90 trees per month on average; the area of needle-leaved forest was reduced by 140 ha (23%) over the same period. Analysis of widely available very-high-resolution satellite images using GEOBIA techniques offers a cost-effective method for detecting changes in tree crown number and land cover in remote mountain valleys; the results provide the information needed to support improved local-level planning and forest management in such areas.


Introduction
The importance of forests and the problems posed by deforestation and forest degradation have been discussed widely in recent years (Hansen et al 2003;FAO 2011;Lambin and Meyfroidt 2011). In the Hindu Kush-Himalayan (HKH) region, forests play an essential role in maintaining the health of the mountain ecosystem and provide mountain people with fuel, food, construction materials, animal feed, and a source of income. They also play an essential role in the hydrological cycle and contribute to reduction of disaster risk (Smith 2013). While forest loss contributes to climate change, forests can also play an important role in mitigation and adaptation through protection and management to increase carbon sequestration (Karki et al 2014).
Despite the widespread recognition of the importance of forests, deforestation and forest degradation continue to be a challenge in the HKH region as elsewhere and are likely to be compounded in future years by climate change. Planning and management to slow or reverse forest loss requires detailed information on forest status and change over time; so does assessment of the impacts of any interventions. But there is less precise information available in the region on the extent of current forest cover or deforestation rates, particularly the detailed information required for planning and assessing locallevel activities. It is especially challenging to gather detailed information on forest extent at ground level in the remote and poorly accessible valleys of the HKH, and even more so to repeat measurements over time. A different approach is required to gather such information.
Satellite images are extremely useful for repetitive timeline-based data collection across large areas, because they offer a synoptic view and the possibility of fast processing of large quantities of data. The high correlation between the spectral bands of satellite images and vegetation parameters make this method the primary source for large-area monitoring of land cover characteristics, especially in remote and rough terrain (Baccini et al 2004;Lu 2006;Uddin et al 2015). It has major advantages over in situ measurements in enabling cost-effective and continuous monitoring of natural resources. More recently, very-high-resolution (VHR) satellite images such as Cartosat, GeoEye, IKONOS, QuickBird, SPOT-5, and WorldView have become available that can be used for detailed monitoring, for example, to extract forest structures and model individual tree parameters (Gibbs et al 2007;Leboeuf et al 2007). Using VHR satellite images is much less costly than carrying out exhaustive ground measurements and can provide extensive, accurate, and timely monitoring and measurements (Immitzer et al 2012). Geographic objectbased image analysis (GEOBIA) provides an opportunity to improve the details in land cover (LC) and land use (LU) maps and assessments of forest parameters (eg tree crown detection and delineation, tree species classification, and identification and monitoring of degraded forest land) (Erikson 2004;Erikson and Olofsson 2005;Ke and Quackenbush 2011;Immitzer et al 2012;Pu and Landry 2012).
A number of algorithms are available for satellite image classification, both pixel-and object-based. Image segmentation techniques were developed in the 1980s but were used to only a small extent in geospatial applications (Blaschke 2010). Following recent rapid development, object-or segmentation-based image classification algorithms now provide much faster and more accurate results than conventional pixel-based algorithms (Stathakis and Vasilakos 2006;Thessler et al 2008;Knorn et al 2009;Blaschke 2010;Quintano and Cuesta 2010). Various GEOBIA techniques have been used for tree crown detection and delineation, including valley following (Ke and Quackenbush 2011), region growing (Culvenor 2002;Ke and Quackenbush 2008), watershed transformation (Wang et al 2004), multiresolution segmentation (Kim et al 2009), wavelet segmentation (Coillie et al 2008), and multiscale object-specific segmentation (Hay et al 2005). Larsen et al (2011) compared 6 different algorithms for tree crown detection under a variety of forest conditions: local maxima detection, valley following, template matching, scalespace theory, Markov random fields, and marked point processes. Detection and delineation of tree crowns can be used to model tree structural variables (eg stem diameter, height, and biomass), which are useful in forest inventory and evaluation of growth performance (Culvenor 2002). Now a number of researchers are starting to use segmentation techniques integrated with active (Light Detection and Ranging [LiDAR]) and passive (multispectral) remote sensing data sets to obtain more accurate and precise results in forested areas (Erdody and Moskal 2010;Kim et al 2010;Mbaabu et al 2014;Thapa et al 2014).
Uniform or consistent data sets in a defined timeframe are essential for monitoring of natural systems that may be influenced by human activities or changes in climate. LCLU maps show the overall changes in existing land features, but VHR satellite images and algorithms are required to understand the dynamics of specific features such as trees within and outside of forested areas. Bai et al (2005) used IKONOS (2006 and WorldView-2 (2011) images for GEOBIA-based quantification and comparative analysis of urban tree species in the city of Tampa, Florida, USA. Ke and Quackenbush (2011) described 3 methods for evaluation and ensuring accuracy of detected and delineated tree crowns through satellite images: (1) how well the delineated crowns represent the position of reference crowns, (2) how well the delineated crown size represents the actual size of reference crowns, and (3) how well the results represent the forest stand properties. According to Ke and Quackenbush (2011), there is no wellestablished method for accurate assessment of detected and delineated tree crowns, but the most commonly used method is to compare the total number of detected trees against reference tree counts at plot level (Culvenor 2002;Ke and Quackenbush 2011). However, one-to-one (canopyto-canopy) assessment is very challenging for tree crown change due to both systematic and nonsystematic errors that occur during the acquisition of satellite images over different periods (Potere 2008).
Understanding the detailed dynamics of LCLU in a small area depends greatly on the spatial resolution of the satellite imagery, because spatial resolution affects image classification results. Following the advent of VHR satellite data (resolution , 5 m), it is now possible to derive detailed LCLU information at a level sufficient to support watershed-level planning. Pu et al (2011) identified significant improvements in classification accuracy using GEOBIA-compared to pixel-based methods-using IKONOS imagery. Similarly, Mathieu et al (2007) used IKONOS imagery and object-based classification to identify 20 LCLU classes for the city of Dunedin, New Zealand, with results similar to those of Pu et al (2011).
Real-time accuracy assessment based on field data is very important for identification, verification, and use of remote-sensing-derived products, although several studies have had to depend on existing field or secondary information due to limitations of resources and time (Stehman and Czaplewski 1998;Foody 2002). Various sources such as coordinates recorded through global positioning system (GPS) receivers, field photographs, topographic sheets, and interviews can be used for accuracy assessment and validation of derived products.
The main objective of the present study was to assess the feasibility of using GEOBIA techniques with VHR satellite images (QuickBird and IKONOS) to detect, delineate, and count tree crowns, prepare LC maps, and detect changes over time, in a cost-effective manner in remote and otherwise poorly accessible mountain valleys of Nepal. The Lorpa watershed in Jumla district was used as an example, with the specific objectives of detecting and delineating tree crowns size in 2006 and 2011, preparing LC maps, and using the results to assess LC change over this 5-year period.

Study area
Jumla is a mountain district located in the Mid-western Development Region of Nepal. The area is historically significant and is considered to be the center of origin of the Nepali language. It is also known for its biodiversity and spectacular scenery. It is one of the least developed districts of Nepal, ranking 68th out of 75 (National Population and Housing Census 2012). The temperature ranges from +33uC in summer to 215uC in winter; average annual rainfall is approximately 800 mm, of which about 450 mm falls between July and September. There is a marked variation in elevation, and the climate ranges from cold temperate at lower elevations (,3000 masl), through alpine at middle elevations (3000-4500 masl), to tundra at the highest elevations (. 4500 masl). The study was conducted in the Lorpa watershed in Jumla district ( Figure 1). The watershed's center point lies at 82u159360E and 29u209240N. It has a total area of 13.05 km 2 , elevation range from 2497 to 4078 masl, and range of slope from 0.5u to 77u. The watershed includes parts of the Dillichaur and Patmara village development committee areas and has 2 main villages: Mathillo Lorpa and Tallo Lorpa. Agriculture is the main source of livelihood, with some additional income from wage labor and small business. There are large areas of forest, some rangeland, and a small area of agricultural land in the flatter part of the lower catchment.

Remote sensing data and preprocessing
A QuickBird satellite image acquired on 18 November 2006 and an IKONOS satellite image acquired on 24 October 2011 were used for change detection. The QuickBird panchromatic sensor collects a single band of data in the spectral range of 405-1053 nm at 0.5 m resolution, and the QuickBird multispectral sensor collects 4 bands of data in the spectral bands of 450-545 nm (blue), 466-620 nm (green), 590-710 nm (red), and 715-918 nm (near-infrared) at 2 m resolution. The IKONOS panchromatic sensor collects a single band of data in the spectral range of 760-850 nm at 1 m resolution, and its multispectral sensor collects 4 bands of data in the spectral bands of 445-516 nm (blue), 506-595 nm (green), 632-698 nm (red), and 757-853 nm (near infrared) at 4 m spatial resolution. The shuttle radar topography mission (SRTM) digital elevation model (DEM), with add-on products such as slope and aspect, was used for topographic information. Base layers (like settlements, roads, and contours) in vector format (shapefiles) were used both as baseline information for the maps and for LC extraction. GPS coordinates and field photographs were used as ancillary information for the identification of LC features on the satellite images.
The IKONOS and QuickBird images were orthorectified individually band by band using rational polynomial coefficient files along a horizontal 90 m (3 arc second) SRTM DEM by applying the cubic convolution method in zone 44 of the UTM (universal transverse Mercator) coordinate system. WGS (world geodetic system) 84 was chosen as datum and spheroid. Although there is a large difference between the resolution of the VHR satellite images and that of the SRTM DEM, the STRM has been found to be better than any other DEM available in the public domain (Suwandana et al 2012). The accuracy of the rectified images was evaluated using GPS readings recorded on the ground.
The multispectral imagery's resolution was improved through a resolution merge process with the help of a higher-spatial-resolution panchromatic band for better spatial visualization and further processing. However, positional differences were observed after the images were overlaid. Sensor azimuth angles with 47.37u deference were observed to be 210u and 162.63u for IKONOS and QuickBird images, respectively. To overcome this positional problem, the multispectral IKONOS image was further rectified. In the 2 data sets, 26 ground control points were selected; they showed an overall root mean square error of 1.2 m for panchromatic images and 1.5 m for multispectral images. Prior to segmentation, low-pass median filters were applied to minimize oversegmentation and smooth the appearance (Platt and Schoennagel 2009). Convolution 333 low-pass filters were used to reduce local variation, remove noise, enhance tree features, and improve the quality of the analyzed satellite images.
Tree crown detection, delineation, and change analysis Tree crowns were delineated using QuickBird and IKONOS images applying the region-growing technique in eCognition Developer. Prior to segmentation, a lowpass filter was applied to smooth the image to avoid oversegmentation (Platt and Schoennagel 2009). This filter produces more homogenous image segments and may reduce the amount of convolutions in the final segmented polygons (Mora et al 2010). In this study, we developed a local maxima algorithm to identify tree tops. From the tree top, the segments were grown to identify tree crown area ( Figure 2). The first step in regiongrowing was to create image objects using multiresolution segmentation. Second, we identified the brightest pixels, which were then used as seed pixels (tree tops). Regions were "grown" from the seed pixels to the local minima, resulting in homogeneous objects based on predefined homogeneity criteria (Erikson 2003;Shih and Cheng 2005;Cui et al 2008).
The tree crown maps provided information on the occurrence of trees in both forest and nonforest categories and enabled assessment of changes over time in the overall tree count and the number of trees in different crown density categories. A 1-ha homogenous grid was created for the entire study area to compare the

Land cover classification and change analysis
For the QuickBird and IKONOS image classification, the GEOBIA technique was used following Bajracharya et al (2010) and an algorithm developed by Uddin et al (2015). The image segments thus developed were used as the basic unit of classification, and each segment was attributed to a particular LC class using a set of defined rules based on information from spectral bands, indices, shape, color, and contextual parameters. The rules were developed for each QuickBird and IKONOS scene to overcome the inherent variability that arises in reflectance values across scenes because of the different times of acquisition. A minimum of 5 reference segments were selected for each class using field information to develop the rules. The image metrics for a given class were chosen according to their relative importance and potential to delineate the class. The thresholds for each metric were fixed iteratively by visual analysis using 2-dimensional feature plots and validating with reference segments. The rules with 90% class separability (based on reference segments) were chosen for application over the entire scene. After ascertaining the class separation using the segment-based approach, classification was performed to obtain an LC classification of both scenes. LC detection using satellite data is an attempt to identify natural and human impacts on earth (Bajracharya et al 2010;Qamer et al 2012). The LC maps from 2006 and 2011 were evaluated and the area of the different classes compared. A change matrix and cross-tabulated values were used to identify changes over the 5 years. For the change matrix, the original 10 classes (Table 1) were reduced to 8, with the "built-up," "barren," and "river" classes subsumed into "other." In a change matrix, diagonal values show the stability of a class, while commission and omission values indicate a shift in area or percentage between classes (Olofsson et al 2013). Using the change matrix, "gross change, forest to nonforest," "gross change, nonforest to forest," and "net change" were extracted at watershed level. Changes in LC between 2006 and 2011 were also identified visually.

Accuracy assessment
Ten 1-ha windows were selected randomly for accuracy assessment of the detected and delineated tree crowns in the 2006 QuickBird and 2011 IKONOS images. Tree crowns were delineated manually and counted in each window using an image interpretation technique and then delineated and counted automatically in the same windows using the region-growing technique (Tiede et al 2006;Tiede et al 2007). A linear regression line was plotted between the 2 values and the coefficient of determination calculated.
Field verification was carried out at 60 randomly selected points in the watershed to assess the accuracy of the 2006 and 2011 LC maps. The geographic coordinates of the points were recorded together with the identified FIGURE 2 Image analysis to delineate tree crowns: (A) satellite image; (B) local maxima/tree tops; (C) distance from tree tops; (D) tree crowns. There are no other water bodies in this watershed.

MountainResearch
LC category. The sample points covered a range of locations and features but were insufficient to fully assess the LC accuracy. Thus, 80 additional randomly distributed points were selected in the raw IKONOS satellite image, the LC features identified by visual interpretation, and the results compared with the class on the LC map, giving a total of 140 points for the accuracy assessment.   determination between algorithm delineated and visually recognized and delineated tree crowns was 0.97 for QuickBird and 0.99 for IKONOS (Supplemental material, Figure S2: http://dx.doi.org/10.1659/MRD-JOURNAL-D-14-00074.S1). Forest was the largest LC class in the watershed ( Figure 5). The total area of needle-leaved forest decreased from 634 ha in 2006 (48.59% of the total land area) to 493 ha in 2011 (37.78% of the total land area) ( Table 2). The total area of the second largest class, grassland, increased from 347 ha in 2006 to 436 ha in 2011, an increase of 6.83% of the original cover. The proportion of the total area covered by closed shrubland decreased by 4.52% between 2006 and 2011, while that covered by open shrubland increased by 6%, an overall increase of 5%, and an 80% increase in the original shrubland area. The proportion of agricultural land decreased slightly from 79 to 72 ha. There was no change in the total area of the settlement, barren land, and river/ riverbed classes. In terms of changes, 166 ha changed from forest to nonforest and 26 ha from nonforest to forest, a net loss of forest of 140 ha ( Figure 5). These results show extensive deforestation and forest degradation in the watershed during the 5 years covered by the study, with an average loss of 90 trees per month, or more than 4% of the cover per year.
The accuracy of the LC maps was assessed by comparing the LC class in the map with the LC class identified by manual assessment at 140 randomly distributed points (60 on the ground and 80 in the IKONOS image). The accuracy assessments are shown in Tables 3 and 4.

Discussion
This study was conducted in the small Lorpa watershed in Jumla district, Nepal, to assess and monitor natural resources, especially forest. The results indicate that tree cutting remains an ongoing problem in the watershed, especially cutting of the larger, mature trees. During the late 1980s, marketing problems and poor crop returns forced many farmers to cut down their mature trees (Tiwari 2010). According to our remote sensing analysis, some of the loss is compensated by increased numbers of young trees growing to a canopy size of 25-50 m 2 and regeneration of trees with a canopy size of ,25 m. But regeneration was slightly lower in 2011 than in 2006, and the new growth did not compensate for loss of mature trees in terms of either number or area.
We used the region-growing technique for the detection and delineation of tree crowns from commercially available QuickBird 2006 and IKONOS 2011 images. A number of authors have confirmed that this technique provides higher accuracy than other algorithms. Cui et al (2008) noted that it is a robust method that can extract objects and boundaries smoothly. Erikson and Olofsson (2005) compared the application of 3 different tree crown delineation methods-template matching (Olofsson et al 2006), 2 region-growing algorithms based on fuzzy rules (Brandtberg 2002), and Brownian motion-and found that all 3 approaches showed an overall accuracy of around 80%. Ke and Quackenbush (2008) concluded that region growing was more accurate than other methods in a study comparing region growing and valley following, as did Erikson and Olofsson (2005), who compared region growing with other algorithms. Similarly, Hussin et al (2014), who compared the region-growing and valley-following approaches, concluded that the region-growing approach showed better delineation than the valley-following approach because region growing is more flexible in detecting tree crowns of different sizes, since it uses both local maxima and local minima. Culvenor (2002) also found that the total number of trees delineated in a needle-leaved forest using a region-growing technique was almost identical to the number of trees counted on the ground (Ke and Quackenbush 2011). The greater accuracy of the technique may be because the method allows different parameters (eg crown width and roundness) to be set.
Other authors have also used commercially available VHR satellite data to obtain information on such factors as tree numbers (Mathieu et al 2007;Zhou et al 2013), tree species (Erikson 2004;Immitzer et al 2012;Latif et al 2012), and forest canopy structure and health (Pouliot et al 2002;Lé vesque and King 2003;Chubey 2006) in other parts of the world for use in forest management. GEOBIA techniques have generally been found to be much more suitable for the detection and delineation of tree crowns and evaluation of LCLU than pixel-based classification algorithms (Zhou et al 2013) and were therefore the method of choice in the present study, as exemplified above. Species-level classification could be performed, but in this study, 10 optimum LC classes were mapped and compared without compromising the accuracy of each class. During digital image classification, the analyst incorporated field experience.
The change in forest cover could be seen clearly in the satellite images (Supplemental material, Figure S3: http://dx.doi.org/10.1659/MRD-JOURNAL-D-14-00074.S1). Overall in Nepal, changes in forest policy have led to immense progress in community-based conservation and restoration of forests (Agrawal and Chhatre 2006;Niraula et al 2013). As elsewhere in the HKH region, local livelihoods in the Lorpa watershed are closely linked to forests and scrub vegetation. Villagers, who need trees for fuel and for construction, have a vested interest in maintaining forest but are also dependent on it (Gilmour 1990). The study area is dominated by rangeland, with only small areas of forest  Total (  and arable land. Water, fuelwood, and other forest resources are scarce, and their availability is further reduced by the increasing demand from the 3 local communities. Agricultural production relies entirely on irrigation, and the watershed is being increasingly used at unsustainable levels, with overgrazing and overharvesting of fuelwood and timber leading to ecosystem degradation (Paudel 2007).

MountainResearch
In remote areas, there may be less awareness among policy-makers of the existence and scale of forest loss and of the long-term implications, because repeated groundbased measurements are intensive and costly and thus virtually nonexistent (Lu 2006). Alternative possibilities for meeting immediate needs for fuelwood and timber are also considerably more limited in remote areas, and even where guidelines or laws exist, possibilities for enforcement are limited (Mather 2007). The area is dominated by rangeland, with only small areas of forest and arable land. Water, fuelwood, and other forest resources are scarce, and their availability is further reduced by the increasing demand from the 3 local communities. Agricultural production relies entirely on irrigation, and the watershed is being increasingly used at unsustainable levels, with overgrazing and overharvesting of fuelwood and timber leading to ecosystem degradation (Mather 2007).
There was some indication that forest management in this area was poor. Participatory appraisal research with the local communities revealed that the current forest governance norms were inadequate, and followup by local line agencies on community-based forest operational plans was very limited (Acharya et al 2008). It appeared that timber was being "leaked" out of Jumla district, both legally and illegally. Remoteness, inadequate forest governance, and authorized and unauthorized overuse of local forests in Lorpa have resulted in continuing deforestation and forest degradation (Paudel 2007). LC maps and information derived from satellite images of the type reported here can help increase awareness and understanding of the problems, support the development of appropriate management plans, and provide a low-cost means for detailed monitoring of forest status.
To be low cost, satellite images have to be widely available. However, obtaining appropriate satellite images, free of cloud cover, snow, or other defects from the same source and month in different years for a particular area, can be problematic. In this study, the images available for the different years were from different sources and a different month (QuickBird from 18 November 2006 and IKONOS from 24 October 2011) and had a different spatial resolution (0.5 m for QuickBird and 1 m for IKONOS). However, it is unlikely that these differences had any marked effect on the results. The 3 weeks' difference in the date is unlikely to have had any noticeable effect on the analysis, because the coniferous trees retain their needles, and the visual characteristics of other classes also vary little at this time of year. The difference in resolution is also unlikely to have affected the tree count, because the great majority of trees had a tree crown diameter considerably greater than the limits of resolution (Ke and Quackenbush 2011). This assumption was further supported by the very high level of accuracy determined in the tree counts in the lowerresolution IKONOS image. Ke and Quackenbush (2011) noted that one problem often encountered with GEOBIA algorithms is off-nadir imagery. Through most of the high-resolution optical stereo images, a high-resolution DEM can be produced for orthorectification. However, these high-resolution DEMs are very costly, so the study used the SRTM DEM with integration of rational polynomial coefficient files (St-Onge et al 2008). Distortion was measured in the QuickBird image using the ground points as a reference. A deviation of 47.37u was detected in the sensor azimuth angle, which can lead to a spatial distortion in the tree comparison. To overcome this problem, the IKONOS image was further rectified geometrically, following orthorectification, using the reference points in the QuickBird image.
There are a number of effects that can cause difficulties in the interpretation of satellite images and detection and delineation of tree crowns. These include the effects of shadow from surrounding hills and mountains, the intermingling of tree crowns, topography, and GPS readings from the field (Dare 2005). In this study, some of these were addressed through manual scanning of the images, followed by geometric correction. Although our algorithm takes advantage of both spectral and spatial information to separate individual trees and to locate treetops, some potential problems deserving further research still exist (Wang et al 2004).
Tree contour detection methods using satellite images are sensitive to noise, or texture, especially in mountainous regions, so preprocessing is essential (Immitzer et al 2012). Blaschke et al (2014) used the term "trivial" in the context of classification where the object has a unique spectral "signature," but in reality it is very difficult to retain a similar spectral response for a particular feature throughout an image. Although VHR satellite images provide very refined objects through segmentation, it is still difficult to extract the range of spectral values for particular LC features . In this study, a multispectral image (4 bands) with spectral indices and DEM-based derived information were used as a basis for rule-based GEOBIA, which improved detection and delineation of the tree crowns.
Other possibilities apart from optical satellite data are also becoming available to improve detection and delineation of individual tree crowns. Integration of active (LiDAR, radar) and passive satellite data has been shown to produce more effective and efficient tree crown identification (Leckie et al 2003;Ke and Quackenbush 2011). Previously one of the major limitations of VHR satellite data was the limited number of spectral bands (Herold et al 2003), but this has now been overcome with WorldView-2 and similar products. In the future, species-level identification, monitoring, and assessment will be able to provide an effective means of support for improved forest management and associated local decision-making (Immitzer et al 2012).

Conclusion
This study showed that it is possible to use GEOBIA techniques with widely available VHR satellite images (QuickBird and IKONOS) to detect, delineate, and count tree crowns, prepare LC maps, and detect changes in tree number and LC over time, in a cost-effective manner, for the remote mountain valleys of Nepal. The GEOBIA region-growing algorithm provided effective and accurate results on the detection and delineation of tree crowns. The results showed a high level of deforestation (approximately 90 trees per month on average over 5 years), which appeared to be related to poor forest governance and the remoteness of the area. This was somewhat unexpected, as in more accessible areas of the country, deforestation and forest degradation are reported to have been reduced and even reversed in recent years as a result of the successful community forestry program.
Although there were a number of potential sources of error in the results, including orthorectification, incident angle, shadows, topography, and GPS readings, these could be addressed through knowledge of the ground realities and visual interpretation of images. The overall accuracy achieved in this study was extremely high, although monitoring at tree species level remains to be assessed. Thus the methodology provides a lowcost and feasible way to obtain the information needed to raise community awareness and support planning and decision-making for improved forest management in remote and poorly accessible mountain watersheds.

ACKNOWLEDGMENTS
We acknowledge the support of the Asian Development Bank (ADB) for the research through funding of the High Mountain Agri-business and Livelihood Improvement (HIMALI) Project (HIMALI). Thanks are due to Renate Fleiner, Govinda Joshi, and Madhav Dhakal for providing extensive field information for validation of results. Our special gratitude goes to Basanta Shrestha, director of strategic cooperation, and Birendra Bajracharya, program manager, at Mountain Environment Regional Information System (MENRIS), for their encouragement and support. We also thank our project partner Local Initiatives for Biodiversity, Research and Development (LI-BIRD) for their constructive and gracious cooperation. Our sincere thanks go to A. Beatrice Murray for language editing of the manuscript.