WPS8569 Policy Research Working Paper 8569 Potential Forest Loss and Biodiversity Risks from Road Improvement in Lao PDR Stephen Danyo Susmita Dasgupta David Wheeler Development Economics Development Research Group & Environment and Natural Resources Global Practice August 2018 Policy Research Working Paper 8569 Abstract This paper develops and applies a spatial econometric model through shorter transport times to market and lower vehi- that links road upgrading to forest clearing and biodiversity cle maintenance costs. The estimated response parameters loss in the Lao People’s Democratic Republic. The paper and a composite biodiversity indicator are used to compute uses 500-meter cells to estimate the relationship between an index of expected biodiversity loss from road upgrading the rate of forest clearing in a cell and its distance to the in each 500meter cell. The results identify areas in the Lao closest point on the nearest road link, the quality of that People’s Democratic Republic where high expected bio- link, the cell’s legal protection status, transport cost to the diversity loss may warrant additional protection as road nearest urban center, the agricultural opportunity value of upgrading continues. This analysis will help policy makers the land, and terrain elevation. The parameter estimates are in the country to weigh context-specific trade-offs between all robust, with the expected signs and very high statistical development and conservation objectives associated with significance. The paper highlights the results that measure road improvement. the impact of improved road quality on forest clearing This paper is a product of the Development Research Group, Development Economics and the Environment and Natural Resources Global Practice. It is part of a larger effort by the World Bank to provide open access to its research and make a contribution to development policy discussions around the world. Policy Research Working Papers are also posted on the Web at http://www.worldbank.org/research. The authors may be contacted at sdasgupta@worldbank.org. The Policy Research Working Paper Series disseminates the findings of work in progress to encourage the exchange of ideas about development issues. An objective of the series is to get the findings out quickly, even if the presentations are less than fully polished. The papers carry the names of the authors and should be cited accordingly. The findings, interpretations, and conclusions expressed in this paper are entirely those of the authors. They do not necessarily represent the views of the International Bank for Reconstruction and Development/World Bank and its affiliated organizations, or those of the Executive Directors of the World Bank or the governments they represent. Produced by the Research Support Team Potential Forest Loss and Biodiversity Risks from Road Improvement in Lao PDR Stephen Danyo* Senior Environmental Specialist, World Bank Susmita Dasgupta Lead Environmental Economist, Development Research Group, World Bank David Wheeler Senior Fellow, World Resources Institute Keywords: Deforestation; Biodiversity Loss; Infrastructure; Road; Green Growth, Forest, Protected Areas, Lao PDR JEL Classification: Q23; R49; Z3 * Authors’ names in alphabetical order. Corresponding Author: Susmita Dasgupta. Email: sdasgupta@worldbank.org; Telephone: +1-202-473-2679 We are thankful to Qinghua Zhao, Senior Information Officer, Development Research Group, World Bank, for development of the new World Bank biodiversity database used in this paper. The findings, interpretations, and conclusions expressed in this paper are entirely those of the authors. They do not necessarily represent the views of the International Bank for Reconstruction and Development/World Bank and its affiliated organizations, or those of the Executive Directors of the World Bank or the governments they represent. 1. Introduction Protected-area strategies seek to conserve biodiversity in tropical forests by restricting infrastructure development in legally-demarcated zones. However, conflicts can arise when forested areas have significant agricultural potential. Within conservation areas, regulatory enforcement may be weakened by dominant economic interests. And some areas with high biodiversity value may be left unprotected so that development can proceed. By implication, both development and conservation may be hindered by a dichotomous policy regime that treats large areas as either completely protected or completely open for development. In reality, valuable economic and ecological resources have non-uniform, overlapping spatial distributions. Policy makers in the Lao People’s Democratic Republic should therefore be equipped to weigh context- specific trade-offs between development and conservation objectives. This paper seeks to contribute by focusing on the trade-offs associated with road improvement. We develop and apply a high-resolution spatial model of road improvement impacts that includes both ecological risks and the economics of forest clearing. The remainder of the paper is organized as follows. Section 2 reviews prior empirical research on the economics of forest clearing globally. In Section 3, we motivate our exercise with a theoretical economic model of road improvement and deforestation. Sections 2 and 3 largely replicate the discussion in Dasgupta and Wheeler (2016). Section 4 describes our spatially- formatted database for Lao PDR, while Section 5 specifies and estimates a deforestation model that incorporates the impact of road improvement. In Section 6, we prepare for the assessment of ecological impacts by developing a composite high-resolution ecological risk indicator for Lao PDR. Section 7 uses the results of Sections 5 and 6 to develop high-resolution estimates of potential ecological risks from a national road improvement program in Lao PDR. Section 8 summarizes and concludes the paper. 2 2. Prior Research Empirical research has provided many useful insights about the determinants of forest clearing. The results are generally consistent with an economic model in which the conversion of forested land varies with potential profitability. Nelson and Chomitz (2009) and Rudel et al. (2009) have studied this relationship across countries over multi-year intervals. Within countries, numerous econometric studies have estimated the impact of economic, social and geographic drivers on deforestation during multi-year intervals. Some studies have used aggregate data for states, provinces or sub-provinces (e.g., studies for Brazilian municipios by Pfaff (1997) and Igliori (2006), and Mexican states by Barbier and Burgess (1996)). Many studies have also used GIS-based techniques to obtain multi-year estimates at a higher level of spatial disaggregation (e.g., Cropper et. al. (1999, 2001) for Thailand; Agwaral et al. (2002) for Madagascar; Deininger and Minton (1999, 2002), Chowdhury (2006) and Vance and Geoghegan (2002) for Mexico; Kaimowitz et al. (2002) for Bolivia; and De Pinto and Nelson (2009) for Panama). In rarer cases, studies have used annual national or regional aggregate time series for extended periods (e.g. Zikri (2009) for Indonesia; Ewers et al. (2008) for Brazil). While econometric work on long-run deforestation drivers is well-advanced, previous data problems limited treatments of economic dynamics to theoretical work and simulation. Arcanda et al. (2008) and others studied the theoretical relationships between macroeconomic drivers and forest clearing. Notable simulation exercises include Cattaneo (2001) for Brazil and San et al. (2000) for Indonesia. Recently, the advent of monthly and annual remote sensing databases has led to the first spatial estimation exercises that explicitly incorporate economic dynamics (Wheeler et al., 2011; Dasgupta et al., 2014). Direct impact studies in Latin America using the new databases have 3 included high-resolution work on new road construction and deforestation in Brazil, where satellite monitoring has been available for a longer period (Laurance et al., 2009) and Bolivia, Panama, Paraguay and Peru (Reymondin et al., 2013). More recently, Li et al. (2015) have investigated the potential impact of decreased travel time from improvement of road links in the Democratic Republic of Congo, but their cross-sectional analysis employs land use information (JRC 2003) that predates the new high-resolution satellite data. Damania et al. (2018) and Dasgupta and Wheeler (2016) have advanced the state of the art by estimating the potential impact of road upgrading on forest clearing and biodiversity in Bolivia, Cameroon, Myanmar and eight Congo Basin countries. Their deforestation models incorporate the latest high-resolution satellite data on forest cover changes, as well as quality estimates for specific links in the regional road system. They also create a biodiversity threat index that combines and synthesizes several measures of potential biodiversity loss. This paper builds on the work of Damania et al., (2018) and Dasgupta and Wheeler (2016), while extending it with a significantly-broader index of potential biodiversity loss and explicit adjustment of biodiversity risk assessment to incorporate the latest information about remaining forested areas. 3. Modeling the Economics of Road Improvement and Deforestation From a formal analytical perspective, this research treats road improvement as a problem of competitive selection among corridors that traverse the same region. The corridors differ in length, construction cost conditions, biodiversity value and potential agricultural income. The optimum corridor choice reflects maximization of a social utility function that values both income and biodiversity, subject to a fixed budget constraint, feasible road quality improvement in each corridor (reflecting the budget constraint), and the corridor-specific impacts of road quality 4 improvement on potential biodiversity loss, expected income growth in the corridor, and expected income growth from increased trade between areas connected by the corridor. We motivate the following econometric analysis with a model that specifies the basic determinants of road improvement impacts and explores their economic and ecological dimensions.1 For expositional clarity and simplicity, we use the Cobb-Douglas (constant- elasticity) specification for all profit and cost functions. This specification is linear in the logs of equation variables and is a first-order approximation to more general specifications (CES, translog, etc.). To simplify the theoretical analysis of alternatives, we assume a constant road improvement budget. Applied assessments of road improvement alternatives might also consider cases where road quality is held constant and improvement expenditures differ because corridors differ in length. In such cases, estimated improvement costs would be subtracted from estimated income gains to provide the appropriate comparison. 3.1 Interregional Trade We specify the economic value of trade between areas connected by a road corridor with a standard gravity model, augmented by a measure of road quality. For prior applications of this approach to overland trade expansion, see Damania et al., (2018), Dasgupta and Wheeler (2016), and Buys, Deichmann and Wheeler (2010). ∝ ∝ ∝ ∝ 1 ∝ where α0, α1, α2, α3, α4 > 0 Tij = Value of trade between areas i and j Ei, Ej = Economic scale of areas i and j dij = Road distance between areas i and j qij = Quality of the road linking areas i and j (0≤q≤1, where 0 → forest track and 1 → paved 1st class road) 1 The theoretical model in this section was first developed by one of the co-authors (Wheeler) in collaborative work with Richard Damania on road upgrading and forest clearing in the Congo Basin countries (Damania et al., 2018). Many thanks to Richard for his insights and suggestions. 5 An improvement in road quality (q) lowers transport costs and increases the potential profitability of agricultural production in the road corridor. The budget for maintaining a particular road quality level is given by: 2 (γ1 > 1) For a fixed road quality budget B: 3 In (3), given a fixed road quality budget, road quality declines with road length. The rate of decline depends on the cost elasticity of road quality. Substituting (3) into (1) and simplifying, an increase in the road improvement budget (B1 → B2) has the following impact on trade between areas i and j: ∝ 4 The interpretation of (4) is straightforward: Trade between areas i and j will expand more or less proportionately with the road improvement budget, depending upon whether the road quality elasticity of trade is greater than or less than the cost elasticity of road quality. 3.2 Local Profitability and Forest Clearing An improvement in road quality lowers transport costs and increases the expected profitability of agricultural production in a road corridor. The expected profitability of agricultural production increases with the size of a cleared parcel, which is determined by the distance of the forest clearing limit from the road. With both factors taken into account, the expected profit function for each unit of road frontage is: 6 5 where β0, β1, β2 > 0, 0 < β3 ≤ 1 π = Potential profit from land clearing v = Potential agricultural value per hectare h = Distance of forest-clearing limit from the road The proprietor of each road-front parcel also confronts incremental forest clearing costs that are constant or increasing with distance from the road and, for agricultural production, incremental commodity transport costs that increase with distance from the road. The proprietor’s composite cost function is: 6 (δ1 > 1) With (5) and (6) specified as present values, the optimal clearing distance from the road is determined by the equality of marginal expected profits and clearing costs. Taking appropriate derivatives, substituting and simplifying, the optimal clearing distance is given by: ∗ 7 where / , δ1 - β3 > 0 (from (5) and (6)) and θ is a composite of equation constants. In (7), the composite elasticities of B, v and d are positive, positive and negative, respectively. Responsiveness increases in each case as the difference narrows between the cost and profit elasticities of forest clearing (δ1 and β3). The responsiveness of forest clearing to greater spending on road quality and greater road distance increases with the ratio of the profit and cost elasticities of road quality. Appropriate substitution and simplification yield the following expression for total profit in the road corridor: ∗ 8 П 2 where P is a composite of equation constants and (> 1). 7 In (8), total expected profits rise less than proportionally with road distance, since expected profits for each roadside land parcel fall as road distance increases. The sign of the distance elasticity in (8) is ambiguous, since ρμ can be greater or less than one. In the former case, total profits will actually fall as road distance increases. From (8), an increase in the road improvement budget (B1 → B2) has the following impact on total profits in a road corridor: П 9 П In (9), total road corridor profits increase more or less proportionately with the road quality budget, depending upon whether ρμ is greater or less than one. The overall interpretation of (9) is as follows: For a road corridor, the responsiveness of total profits to an increase in the road quality budget is positively related to the ratio of profit and cost elasticities of road quality, positively related to the profit elasticity of forest clearing, and negatively related to the difference between the cost and profit elasticities of forest clearing. 3.3 Biodiversity Loss Substituting from (7), total biodiversity loss in a road corridor is given by: ∗ 10 2 2 where Ɛ = measured biodiversity intensity. In (10), the relationship between total biodiversity loss and road distance depends on whether the ratio of the profit and cost elasticities of road quality is greater than the difference between the cost and profit elasticities of forest clearing. If this is the case (given a fixed budget for road quality), total biodiversity loss will decline with road distance. If the converse is true, increasing road distance will increase biodiversity loss. From (10), the relationship between a road budget increase and total biodiversity loss is given by: 8 11 . As above, biodiversity loss will increase more than proportionately with a road budget increase if the ratio of profit and cost elasticities of road quality is greater than the difference between the cost and profit elasticities of forest clearing. For a fixed road quality budget, the relative profit and biodiversity loss equations for two competitive road corridors (A and B) are as follows: П 12 П 13 Ceteris paribus (i.e., for identical road budgets and biodiversity intensities), the corridor with greater potential agricultural value per hectare will have both higher total profits and greater total biodiversity loss. The effects of road distance are more complex for two reasons. First: 1 1 14 0 Given (14), equations (12) and (13) incorporate three possibilities. Let η1 = ρμ and η2 = μ/(δ1-β3). If η1 and η2 are both less than one, the corridor with the longer road will have greater total profits and greater biodiversity loss. If η1 and η2 are both greater than one, the corridor with the longer road will have lower profits and lower biodiversity loss. In the third possible case -- η1 > 1 and η2 < 1 -- the corridor with the shorter road will be favored unambiguously because it will have higher profits and lower biodiversity loss. It is important to note that this is a ceteris paribus result, since the corridor with the shorter road could have sufficiently greater biodiversity intensity to reverse the result. A complete calculation would also have to include consideration of trade- related regional income, which would favor the shorter road. Given the relatively complex set of 9 factors that affect these results, the impacts of road improvement on income gains, forest clearing and biodiversity loss seem likely to vary significantly across countries, and across regions within countries. In this paper, we focus on the forest clearing and biodiversity components of the problem. 4. Data 4.1 Road Networks Table 1 and Figure 1 summarize 2015 road network information for Lao PDR in digital maps (shapefiles) provided by Delorme, Inc. As Table 1 shows, the maps identify 2,055 road network links. Network links are graded by status: primary roads (purple in Figure 1), secondary roads (red), and tertiary roads (orange). Table 1: Delorme road links by network status: Lao PDR Road Type Links Primary 370 Secondary 309 Tertiary 1,376 Total 2,055 10 Figure 1: Lao PDR road networks 11 4.2 Forest Clearing Hansen et al. (2013) publish annual high-resolution estimates of global forest clearing. The data are currently available at 30 m spatial resolution for 2001-2016. We use the Hansen estimates to compute forest clearing rates through 2016 in 500 m cells.2,3 Figures 2 and 3 display our results for Lao PDR in 2000 and 2016, aggregated to 0.1-degree cells for ease of interpretation. Each cell is color coded by total percent cleared in the relevant year. Figure 2 (2000) reveals an initial pattern of forest clearing in road corridors that becomes much more pronounced in Figure 3 (2016), particularly in areas of high road network density. 5. Estimation of the Forest Clearing Model The estimation exercise in this paper draws on the insights of previous research and the theoretical model in Section 3 to incorporate five critical determinants of forest clearing in road corridors: distance from the road, road quality, legal protection status, transport time to the nearest market center and terrain elevation. Following Damania et al. (2018) and Dasgupta and Wheeler (2016), we extend previous work on road transport cost and deforestation by introducing explicit estimation of road quality effects. 2 We use the 500 m approximation for expositional convenience. Our cell sides are .005 decimal degrees in length. This translates to 556.6 meters on a side, yielding cell areas of 309,804 square meters or 30.98 hectares. We produce Hansen estimates at .005-degree resolution by resampling Hansen 30 m raster values (1 if clearing occurred from 2001 to 2016; 0 otherwise). 3 Our estimate also incorporates the Hansen 30 m estimates of tree cover in 2000. 12 Figure 2: Road network and forest clearing in Lao PDR, 2000 13 Figure 3: Road network and forest clearing in Lao PDR, 2016 \ 14 5.1 Model Specification We specify our estimation model as follows. (15) a h ε Expected signs: β1, β4, β5, β7 <0 ; β3< β2<0; β6 >0 where ci = Cumulative percent cleared in 2016, cell i di = Distance of cell i from nearest road segment qk = Quality dummy variable for nearest road segment [q2: secondary road; q3: tertiary road] (q1 excluded (effect absorbed in the regression constant) to prevent total collinearity) pi = Legal protection status of cell i [1 if protected; 0 otherwise] ti = Transport time from cell i to the nearest urban market center ai = Agricultural opportunity value of land in cell i hi = Elevation of cell i Ɛi = Random error term A priori, we expect the impact on forest clearing to be negative for distance from the nearest road (β1< 0); decreasing as road quality declines (β3< β2<0); negative for legally‐protected areas (β4< 0); negative for transport distance to the nearest urban center (β5 < 0), positive for the agricultural opportunity value of the land (β6>0); and negative for more elevated terrain (β7 < 0). 5.2 Variable Measures Our data and sources are as follows. Cumulative percent cleared, 2016: this incorporates two Hansen measures: percent cleared prior to 2001 and annual clearing from 2001 to 2016 (Hansen et al. (2013)). We produce estimates at .005-degree resolution by resampling Hansen 30 m raster values. Distance from road segment: distance from the centroid of each cell to the nearest road segment, calculated in ArcGIS 10. 15 Road quality: proxied by primary, secondary or tertiary status assigned by Delorme (Table 1). Legal protection status: 1 if the cell is in a protected area identified by the World Database on Protected Areas (WDPA); 0 otherwise. The WDPA shapefile has been downloaded from http://www.protectedplanet.net/. Travel time: time from the cell centroid to the nearest urban center with a population of 50,000 or greater, as estimated by Uchida and Nelson (2009). Raster resolution: .0083 decimal degrees. Agricultural opportunity value: mean value for a cell, calculated from the high-resolution global grid developed by Deveny et al. (2009). Raster resolution: .0025 decimal degrees. Elevation: Lao PDR data extracted from the CGIAR-SRTM data set (3 arc-seconds (approximately 90 m) resolution), downloaded from http://srtm.csi.cgiar.org. 5.3 Estimation For this exercise, we utilize an 11% random sample to minimize potential estimation bias and inconsistency produced by spatial autocorrelation among neighboring cells.4 The dependent variable is continuously distributed in the range [0,1], so the appropriate estimator is fractional logit.5 Table 2 presents our results, which are uniformly robust: All parameter signs match our prior expectations, and all have high statistical significance. For 500 m Hansen cells, forest clearing declines with distance from the nearest road segment, transport time to the nearest urban center and elevation. We have excluded Delorme primary status, so the results for road quality dummy variables should be interpreted as deviations from the impact of primary road segments: Ceteris paribus, secondary roads have lower impact on forest clearing than primary roads, and the impact of tertiary roads is lower still. As expected, clearing is greater on land with higher agricultural opportunity value. Our result for protected status is striking and holds independent 4 The full data set contains 795,072 observations for Lao PDR, georeferenced by centroids of Hansen-resampled .005- degree cells. We randomly select 100,000 observations from the data set. After elimination of rows with missing values for at least one regression variable, the regression data set contains 90,330 observations -- 11.4% of the full data set. 5 We employ the Stata estimator glm [link(logit); family(binomial); vce(robust)]. 16 interest: Forest clearing is much lower in protected areas than the other model variables would predict, and the effect has high statistical significance. Table 2: Regression results: Forest clearing in Lao PDR Dependent variable: Cumulative percent cleared, 2016 Regression variables: Log distance from nearest road segment -0.208 (63.30)** Secondary road segment -0.151 (10.67)** Tertiary road segment -0.224 (18.65)** Legally protected -0.720 (70.17)** Transport time -0.001 (36.26)** Agricultural opportunity value 0.003 (52.72)** Elevation -0.001 (84.07)** Constant 0.967 (68.68)** Observations 90,330 Robust t statistics in parentheses * significant at 5%; ** significant at 1% 5.4 Implications for Forest Clearing To assess the implications of our results, we upgrade all secondary and tertiary road links to primary status and predict the resulting changes in clearing for all Hansen cells in our database. Figures 4 and 5 display predicted cumulative clearing and the change from 2016, respectively. Predicted impacts vary widely, with little or no change in areas far from the Delorme road network6 6 The econometric results for log (distance from road) in Table 2 quantify a continuous relationship, so upgrading secondary and tertiary roads will have small predicted effects in areas that are relatively far from those roads. 17 and areas dominated by roads that already have primary status. Figure 5 shows that the largest predicted changes in forest clearing are concentrated in northern Lao PDR, with the largest clusters visible in the northwest, southwest and eastern parts of the northern region. Figure 4: Cumulative forest clearing after road upgrading 18 Figure 5: Forest clearing attributable to road upgrading 19 6. Incorporating Biodiversity Drawing on a new World Bank database, we develop a multidimensional biodiversity indicator for Lao PDR to identify concentrations of diverse species and their vulnerability to encroachment and habitat loss. The database incorporates information from thousands of digital range maps (shapefiles) provided by the International Union for the Conservation of Nature (IUCN) and BirdLife International. We also incorporate ecoregion maps provided by WWF International. We use GIS overlays to compute our biodiversity indicator for each Hansen cell in our forest clearing database. Our approach follows the methodology developed in Dasgupta and Wheeler (2016), but draws on the new World Bank database to incorporate many more dimensions of vulnerability and risk. Our approach incorporates separate measures for amphibians, mammals, reptiles, birds and WWF ecoregions. 6.1 Species Density We measure species density as the species count for each 500 m Hansen cell. 6.2 Species Vulnerability Species density provides critical information for assessing ecological risks, but at least two other elements are needed: (1) Geographic vulnerability, which can be proxied by endemicity: the proportion of each species’ range that lies within each cell. Species that reside in very few cells may be particularly vulnerable to habitat encroachment. For this exercise, endemicity treats all species equally at the national level, since each species has a total count of 1. Total endemicity for each cell -- the sum of its species endemicity measures -- assigns higher values to cells inhabited by species whose ranges are relatively limited. By implication, forest clearing in higher-value cells may be particularly destructive for remaining critical habitat. 20 (2) A measure of extinction risk that adds the insights of the international scientific community. We convert IUCN Red List status codes to extinction probabilities using the methodology of Mooers et al. (2008).7 For species indicator construction, we normalize these probabilities so that a weight of 1.0 is assigned to species in the highest category (Critically Endangered). Table 3 tabulates conversions from Red List codes to normalized species weights, using four probability assignments. Three employ IUCN estimates to derive measures of extinction probability over the next 50, 100 and 500 years. The fourth draws on Isaac et al. (2007), who combine a direct extinction risk measure with a measure of each species' isolation on a phylogenetic tree.8 Table 3: Normalized species aggregation weightsa Normalized Extinction Probabilities IUCN: Future Years IUCN Isaacb 50 100 500 Code Status CR Critically Endangered 1.00000 1.00000 1.00000 1.00000 EN Endangered 0.50000 0.43299 0.66770 0.99600 VU Vulnerable 0.25000 0.05155 0.10010 0.39000 NT Near Threatened 0.12500 0.00412 0.01000 0.02000 LC Least Concern 0.06250 0.00005 0.00010 0.00050 Rounded Weight Ratios CR:EN 2 2 1 1 CR:VU 4 19 10 3 CR:NT 8 243 100 50 CR:LC 16 20,000 10,000 2,000 a Data source: Mooers et al. (2008). b From calculations by Mooers et al., based on Isaac et al. (2007). 7 The IUCN's current classification categories are Critically Endangered, Endangered, Vulnerable, Near Threatened and Least Concern. 8 A phylogenetic tree is a branching tree diagram that traces the evolutionary descent of different species from a common ancestor. Species in sparse (isolated) branches of a phylogenetic tree are relatively unique, since they share common descent patterns with fewer other species. 21 Table 3 shows that Isaac’s inclusion of the phylogenetic isolation factor changes the weight ratios substantially, particularly for species in the lowest threat category (Least Concern). We explore the implications for hypothetical areas A and B in Table 4. A is populated by only 2 species, both rated as Critically Endangered. B is populated by 20,000 species, but all are rated as of Least Concern. Our extinction risk indicator for each area is the sum of normalized extinction probabilities for resident species. Assignment of weights for Mooers' IUCN-derived 50-year extinction probabilities yields a total risk indicator of 2 for A -- twice the total for B, because each Critically Endangered species is weight-equivalent to 20,000 Least Concern species. In contrast, assignment of the Isaac weights yields an overall risk rating for B (1,250) that is 625 times greater than the rating for A (2), because each Critically Endangered Species is weight-equivalent to 16 species of Least Concern. The other two cases are intermediate, but far closer to the 50-year IUCN case. Table 4: Implications of alternative weighting schemes Total Scores IUCN Extinction Probabilities: Future Years Status Area Species (Uniform Count Within Areas) Isaac 50 100 500 A 2 CR 2 2 2 2 B 20,000 LC 1,250 1 2 10 6.3 Biome Vulnerability Measures based on animal species alone provide an incomplete accounting of biodiversity. A more complete measure would incorporate plants and insects, using indices similar to those we have developed for animals. Although no such indices exist at the requisite geographic scale, WWF has provided a first approximation by segmenting the world into 825 terrestrial ecoregions 22 (Figure 6). WWF defines an ecoregion as “a large unit of land or water containing a geographically distinct assemblage of species, natural communities, and environmental conditions.”9 Accordingly, we adopt the ecoregion as a general proxy for distinctive plant and inspect species, as well as animal species that are not represented in the range maps provided by IUCN and BirdLife International. Our method for incorporating WWF ecoregions resembles our treatment of species endemicity. For this exercise, we identify all moist and dry forest ecoregions in Lao PDR. We compute the percent of total forest area in Lao PDR accounted for by each ecoregion. Then we compute its vulnerability index as the inverse of its area share and assign the appropriate index value to each cell. This accounting assigns high values to cells in smaller ecoregions, where clearing single cells may pose more significant threats to biome integrity. 6.4 A Composite Biodiversity Indicator The example in Table 4 suggests that alternative vulnerability indicators may yield significantly different risk metrics. It is important for a vulnerability indicator methodology to accommodate different risk-weighting schemes in a consistent and plausible way. From the World Bank database, we obtain 50 spatially-formatted indicators: Overall species indicators (8): total cell counts and total endemicity measures for four vertebrate classes: amphibians, mammals, reptiles and birds10; Species indicators by IUCN status (40): cell counts and endemicity measures for species within each vertebrate class in five IUCN risk status categories (per Table 3): CR (critically endangered), EN (endangered), VU (vulnerable), NT (near-threatened); and LC (least concern). 9 Complete information about the WWF terrestrial ecoregions is available online at http://wwf.panda.org/about_our_earth/ecoregions/. 10 The total indicators include species in IUCN threat status class DD (data deficient). 23 Ecoregion indicators (2): Cell counts and endemicity measures for WWF ecoregions. 24 Figure 6: WWF Terrestrial Ecoregions 25 We recognize the potential diversity of stakeholder views about the relative importance of species-, ecoregion- and risk-related factors in this context. To accommodate diverse concerns, we adopt a conservative strategy for indicator construction. First, we divide each of our 50 indicators by its maximum value and multiply by 100 to create an index in the range 0-100. This ensures comparability in measurement. Then, for each 500 m Hansen cell in our database, we select the maximum value among 50 indices as our risk indicator. This ensures significant representation for the concerns of stakeholders whose focus may range from protection of critically-endangered amphibians to conservation of small Laotian biomes. Figure 7 displays the geographic distribution of the resulting indicator, aggregated to 0.1 degree grid cells for ease of interpretation.11 We have overlaid a province map of Lao PDR as a geographic reference. Figure 7 reveals a dominant wedge-shaped “axis” that narrows from the extreme northwest (where the highest indicator values are concentrated in Louang Namtha and Bokeo) to east-central Lao PDR (with the highest indicator values in Bolikhamxai and Khammouan). Large areas with high indicator values are also evident in Phôngsali, Oudômxai and Louangphrabang. Figure 7 also reveals a much smaller “axis” in southeastern Lao PDR, with the highest indicator values concentrated in the eastern parts of Xékong and Attapu. Pockets of high indicator value are also visible in the southern frontier areas of Champasak and Attapu. These may well reflect small Laotian portions of ranges for species and/or ecoregions that extend into Cambodia. Our methodology focuses on Laotian national considerations, so the small Laotian portions of these regions/ecoregions would be assigned high endemicity value. 11 Our aggregative indicator is the maximum indicator value within each 0.1-degree grid cell. 26 6.5 Accounting for Forest Clearing The composite biodiversity indicator in Figure 7 incorporates range and ecoregion maps that have not been explicitly adjusted for forest clearing. To provide a more accurate view of current biodiversity status, we combine the information sets reflected in Figures 3 and 7. Within each 500 m Hansen cell, we compute the product of the biodiversity indicator value and the estimated percent that remains forested in 2016. We aggregate to 0.1-degree grid cells by summing products within each aggregate cell and normalizing results so that the maximum cell value is 100. Comparison of the result in Figure 8 with Figure 7 shows that incorporation of forest clearing has a very significant effect: The previous, simple “wedge” geometry is replaced by a complex pattern with a broad distribution of the highest values. The highest-value areas in Figure 7 retain importance, but they are joined in Figure 8 by high-value clusters in Houaphan, Vientiane and Xaisômboun. In summary, the effect of incorporating forest clearing is wide dispersal of the remaining high-priority biodiversity conservation areas. In the clearing-adjusted map displayed by Figure 8, every Laotian province has clusters that rank among the highest-priority areas in Lao PDR. 27 Figure 7: Composite biodiversity indicator, Lao PDR 28 Figure 8: Composite biodiversity indicator adjusted for forest clearing, 2016 29 7. Road Upgrading and Ecological Risk We calculate expected biodiversity loss for a cell as the product of its previous-clearing- adjusted biodiversity indicator value (Figure 8) and the change in forest clearing induced by road upgrading (Figure 5). Figure 9 displays the geographic distribution of predicted biodiversity indicator values after road upgrading, while Figure 10 displays the change in indicator value that is attributable to upgrading. The map in Figure 10 has been created by subtracting indicator values in Figure 8 from those in Figure 9, and re-normalizing the result to the range [0,100]. Once again, the results indicate a broad geographic distribution of the largest biodiversity impacts. Significant impact clusters are visible in southern Phôngsali, northwestern Louang Namtha, western Oudômxai, east-central Louangphrabang, northern and southeastern Houaphan, northeastern Xaisômboun, southwest Vientiane, northwest Khammouan, southeast Savannakhét, northeast Champasak and southern Xékong. Clusters of lower but substantial impact are also visible in widely-scattered locations. 8. Summary and Conclusions In this paper, we have developed, estimated and applied a spatially-explicit model that links road upgrading to forest clearing and potential biodiversity loss in Lao PDR. Our approach combines high-resolution forest-clearing information from Hansen et al. (2013); information on biodiversity and protected areas from IUCN and WWF International; digital road maps from Delorme, Inc.; digital elevation maps from CGIAR; and information on agricultural opportunity values from Deveny et al. (2009). Using the 30 m Hansen data, we estimate cumulative forest clearing rates for all 500 m cells in Lao PDR. For each 500 m cell, we measure the distance to the closest point on the nearest road; the quality of the road link; legal protection status; travel time to the nearest urban market center; the mean agricultural opportunity value of the land; and elevation. 30 Figure 9: Previous-clearing-adjusted biodiversity indicator after road upgrading 31 Figure 10: Change in previous-clearing-adjusted biodiversity indicator after ipgrading (Normalized to [0,100]) 32 We use an appropriately-specified econometric model to estimate the effects of these variables on the cumulative percent of forest cleared in 2016. Our estimation exercise uses a small random sample to minimize potential spatial autocorrelation problems. The resulting parameter estimates are extremely robust, with the expected signs and very high statistical significance in all cases. Our result for legal protection status is particularly strong, indicating that forest conservation in protected areas has been robust in Lao PDR. We focus on distance from the nearest road link and its Delorme primary, secondary or tertiary status, which we use as a proxy for road quality. We find an important role for road quality, which reflects the effects of average vehicle speeds and differential maintenance costs in transport to market.12 Econometric measurement of road quality effects enables us to estimate the deforestation impact of programs that upgrade secondary and tertiary roads to primary status. To illustrate our methodology, we explore the impact of upgrading all secondary and tertiary roads to primary roads in Lao PDR. We find a highly non-uniform spatial distribution, with several large impact clusters in northern Lao PDR and a few smaller clusters in southern Lao PDR. With our model-based predictions in hand, we measure biodiversity risk in each 500 m cell using ecoregion information from WWF International and thousands of species range maps from IUCN and Birdlife International. This information is summarized in a new World Bank database that provides georeferenced indicators for ecoregions and species intensity, extinction risk and 12 We cannot completely discount simultaneity bias in this estimation exercise, since anticipated returns from future forest clearing might have provided one motivation for past road construction in Lao PDR. However, the relative paucity of forest clearing beyond established urban centers and principal road corridors as recently as 2000 (Figure 2) strongly suggests that other factors were the primary drivers of historical road construction in Lao PDR. Any remaining upward bias from simultaneity will produce some tendency for our methodology to overestimate the potential impact of road improvement on forest clearing. In any case, since “extinctions are forever”, we believe that invocation of the precautionary principle in conservation is more than sufficient to accommodate some potential over- estimation. 33 endemicity for amphibians, mammals, reptiles and birds. We develop 50 measures that incorporate biome status, species density, endemicity, and extinction risk. We normalize each measure to the range [0,100] and develop a composite biodiversity index for each 500 m cell by selecting the maximum normalized value among the 50 measures for that cell. This approach ensures representation for diverse stakeholder concerns, which may range from preservation of critically-endangered amphibians to conservation of small ecoregions that are endemic to Lao PDR. Available species range maps may not reflect the latest information on forest clearing, which is best represented by the global Hansen data set. Accordingly, we adjust our biodiversity indicator for each 500 m cell by multiplying its biodiversity indicator value by the Hansen estimate of percent forest cover remaining as of 2016.13 The result is our summary estimate of current biodiversity risk. To assess the potential impact of road improvement on biodiversity, we compute an expected loss measure for each cell by multiplying its forest-clearing-adjusted biodiversity index by the predicted increase in forest clearing percent from secondary and tertiary road upgrading. We map the results to identify priority areas with high expected losses. We find highly diverse patterns of expected forest clearing and biodiversity loss that would be difficult to anticipate without this kind of analysis. We believe that these results highlight several important dimensions of ecological risk assessment and road investment policy in Lao PDR. First, newly-available data on forest clearing, species risks and ecoregions provide an unprecedented opportunity to assess ecological impacts at high spatial resolution. Second, assessment can be founded on a rigorous methodology that uses 13 As we explain in the paper, this requires combining the Hansen estimate of forest cover in 2000 with the estimate of forest clearing from 2001 to 2016. 34 Laotian data to derive econometric estimates of responsiveness to road upgrading while controlling for a host of economically-significant forest clearing determinants. Third, a comprehensive biodiversity assessment can incorporate many dimensions of biodiversity risk into a comprehensive index that ensures representation for broad stakeholder concerns. In this paper, we have developed and illustrated an assessment that incorporates all of these factors. We should also note that our approach is not limited to national-scale assessments. Since our methodology uses high-resolution spatial databases, it can also be “zoomed” for assessment of forest-clearing, ecological risks, and road improvement impacts in specific road corridors or small geographic regions. It can also be extended to other road network databases (e.g. open street maps), and to the analysis of potential impacts for proposed roads that do not currently exist. In summary, since road upgrading will inevitably accompany rural development programs, we believe that explicit identification of ecologically-vulnerable areas and transport investment impacts can provide at least two valuable types of information. First, with limited budgets, it can help steer road upgrading programs toward corridors where expected biodiversity losses will be minimized. Second, it can inform the adoption of appropriate protection measures in vulnerable road corridors and neighboring areas. 35 References Arcanda, Jean-Louis, Patrick Guillaumont, Sylviane Guillaumont Jeanneney. 2008. Deforestation and the real exchange rate. Journal of Development Economics, 86(2): 242-262. Barbier, Edward and Joanne C. Burgess. 1996. Economic analysis of deforestation in Mexico. Environment and Development Economics, 1:203-239. Buys, Piet, Uwe Deichmann and David Wheeler. 2010. Road Network Upgrading and Overland Trade Expansion in Sub-Saharan Africa. Journal of African Economies 19(3): 399-432. June. Cattaneo, Andrea. 2001. Deforestation in the Brazilian Amazon: Comparing the Impacts of Macroeconomic Shocks, Land Tenure, and Technological Change. Land Economics 77(2): 219- 240. Chowdhury, Rinku Roy. 2006. Driving forces of tropical deforestation: The role of remote sensing and spatial models. Singapore Journal of Tropical Geography, 27:82–101. Cropper, Maureen, Jyotsna Puri, and Charles Griffiths. 2001. Predicting the Location of Deforestation: The Role of Roads and Protected Areas in North Thailand. Land Economics, 77(2):172-186. Cropper, Maureen, Charles Griffiths and Muthukumara Mani. 1999. Roads, Population Pressures, and Deforestation in Thailand, 1976-1989. Land Economics, 75(1): 58-73. Damania, Richard, David Wheeler, Jason Russ and Alvaro Federico Barra. 2018. “The Road to Growth: Measuring the Tradeoffs between Economic Growth and Ecological Destruction”. World Development. 101(January): 351-376 Dasgupta, Susmita and David Wheeler. 2016. Minimizing ecological damage from road improvement in tropical forests. World Bank Policy Research Working Paper No. 7826. Dasgupta, Susmita, Daniel Hammer, Robin Kraft and David Wheeler. 2014. Vyaghranomics in space and time: Estimating habitat threats for Bengal, Indochinese, Malayan and Sumatran tigers. Journal of Policy Modeling (forthcoming). Deininger, Klaus and Bart Minten. 1999. Poverty, Policies, and Deforestation: The Case of Mexico. Economic Development and Cultural Change, 47:313–344. Deininger, Klaus and Bart Minten. 2002. Determinants of Deforestation and the Economics of Protection: An Application to Mexico. American Journal of Agricultural Economics, 84(4): 943- 960. De Pinto, Alessandro and Gerald C. Nelson. 2009. Land Use Change with Spatially Explicit Data: A Dynamic Approach. Environmental and Resource Economics, 43(2): 209-229. 36 Deveny, A., J. Nackoney, N. Purvis, R. Kopp, M. Gusti, G. Kindermann,, M. Macauley, E. Myers Madeira, M. Obersteiner and A. Stevenson. 2009. Forest Carbon Index: The geography of forests in climate solutions. Joint report by Resources for the Future and Climate Advisers. December Ewers, Robert, William Laurance and Carlos Souza. 2008. Temporal fluctuations in Amazonian deforestation rates. Environmental Conservation, 35:303-310. Hansen, M.C., P. V. Potapov, R. Moore, M. Hancher, S. A. Turubanova, A. Tyukavina, D. Thau, S. V. Stehman, S. J. Goetz, T. R. Loveland, A. Kommareddy, A. Egorov, L. Chini, C. O. Justice and J. R. G. Townshend. 2013. High-Resolution Global Maps of 21st-Century Forest Cover Change. Science 342 (6160): 850-853. Igliori, Danilo Camargo. 2006. Deforestation, Growth and Agglomeration Effects: Evidence from Agriculture in the Brazilian Amazon. Department of Land Economy, University of Cambridge. Isaac N., S. Turvey, B. Collen, C. Waterman and J. Baillie. 2007. Mammals on the EDGE: conservation priorities based on threat and phylogeny. PLoS One 2: e29. Joint Research Center (JRC). 2003. European Commission Global Land Cover 2000 Project. http://bioval.jrc.ec.europa.eu/products/glc2000/glc2000.php. Kaimowitz, D., P. Mendez, A. Puntodewo and J. Vanclay, 2002. Spatial regression analysis of deforestation in Santa Cruz, Bolivia. In C.H. Wood and R. Porro (eds) Land Use and Deforestation in the Amazon. University Press of Florida, p. 41-65. Laurance, William F. 2012. As Roads Spread in Rainforests,The Environmental Toll Grows. Yale Environment 360. 19 January. http://e360.yale.edu/feature/as_roads_spread_in_tropical_rain_forests_environmental_toll_grows /2485/ Laurance, William F., Miriam Goosem and Susan G.W. Laurance. 2009. Impacts of roads and linear clearings on tropical forests. Trends in Ecology & Evolution 24(12): 659-669. 11 September Laurance, William F., Mark A. Cochrane, Scott Bergen, Philip M. Fearnside, Patricia Delamônica, Christopher Barber, Sammya D'Angelo and Tito Fernandes. 2001. The Future of the Brazilian Amazon. Science 291(5503): 438-439. 19 January. Li, Man, Alessandro De Pinto, John M. Ulimwengu, Liangzhi You and Richard D. Robertson. 2015. Impacts of Road Expansion on Deforestation and Biological Carbon Loss in the Democratic Republic of Congo. Environmental and Resource Economics 60: 433-469 37 Mooers, A., D. Faith and W. Maddison. 2008. Converting Endangered Species Categories to Probabilities of Extinction for Phylogenetic Conservation Prioritization. PLoS ONE 3(11): e3700. doi:10.1371/journal.pone.0003700. Nelson, Andrew and Kenneth M. Chomitz. 2009. Protected Area Effectiveness in Reducing Tropical Deforestation: A Global Analysis of the Impact of Protection Status. Evaluation Brief 7, Independent Evaluation Group, World Bank, October. Pfaff, Alexander. 1997. What Drives Deforestation in the Brazilian Amazon? Evidence from Satellite and Socioeconomic Data. Policy Research Working Paper No. 1772, World Bank. May. Reymondin, Louis, Karolina Argote, Andy Jarvis, Carolina Navarrete. Alejandro Coca, Denny Grossman, Alberto Villalba, Paul Suding. 2013. Road Impact Assessment Using Remote Sensing Methodology for Monitoring Land-Use Change in Latin America: Results of Five Case Studies. Inter-American Development Bank Environmental Safeguards Unit, Technical Note No. IDB - TN - 561. July. RFF (Resources for the Future and Climate Advisors). 2011. The Geography of Forests in Climate Solutions: Price Geography. http://www.forestcarbonindex.org/maps.html Rudel, Thomas, Laura Schneider, Maria Uriarte, B. L. Turner II, Ruth DeFries, Deborah Lawrence, Jacqueline Geoghegan, Susanna Hecht, Amy Ickowitz, Eric F. Lambin, Trevor Birkenholtz, Sandra Baptista, and Ricardo Grau. 2009. Agricultural intensification and changes in cultivated areas, 1970–2005. Proceedings of the National Academy of Science. 106 (49): 20675–20680. San, Nu Nu, Hans Löfgren and Sherman Robinson. 2000. Structural Adjustment, Agriculture, and Deforestation in the Sumatera Regional Economy. TMD Discussion Paper No. 52, International Food Policy Research Institute, March. Stifel, D., and B. Minten. 2008. Isolation and Agricultural Productivity. Agricultural Economics 39 (1): 1–15. Stifel, David, Bart Minten and Bethlehem Koro. 2012. Economic Benefits and Returns to Rural Feeder Roads: Evidence from a Quasi-Experimental Setting in Ethiopia. International Food Policy Research Institute and Ethiopian Development Research Institute. ESSP II Working Paper No. 40. May Uchida, H. and A. Nelson. 2009. Agglomeration Index: Towards a New Measure of Urban Concentration. Background paper for World Development Report 2009. Washington, DC: World Bank. University of Maryland. 2014. Global Forest Change 2000–2012: Data Download. Department of Geographical Sciences. Available online at http://www.earthenginepartners.appspot.com/science-2013-global-forest/download.html 38 Vance, Colin and Jacqueline Geoghegan. 2002. Temporal and spatial modelling of tropical deforestation: a survival analysis linking satellite and household survey data. Agricultural Economics, 27(3): 317-332. Wheeler, David, Daniel Hammer, Robin Kraft, Susmita Dasgupta, Brian Blankespoor. 2011. Economic Dynamics and Forest Clearing: A Spatial Econometric Analysis for Indonesia. Center for Global Development Working Paper No. 280. http://www.cgdev.org/content/publications/detail/1425820 Zikri, Muhammad. 2009. An Econometric Model for Deforestation in Indonesia. Working Paper in Economics and Development Studies, No. 200903. Center for Economics and Development Studies, Department of Economics, Padjadjaran University. July. 39