Mark Broich

Time-series analysis of multi-resolution optical imagery for quantifying forest cover loss in Sumatra and Kalimantan, Indonesia

Monitoring loss of humid tropical forests via remotely sensed imagery is critical for a number of environmental monitoring objectives, including carbon accounting, biodiversity, and climate modeling science applications. Landsat imagery, provided free of charge by the U.S. Geological Survey Center for Earth Resources Observation and Science (USGS/EROS), enables consistent and timely forest cover loss updates from regional to biome scales. The Indonesian islands of Sumatra and Kalimantan are a center of significant forest cover change within the humid tropics with implications for carbon dynamics, biodiversity maintenance and local livelihoods. Sumatra and Kalimantan feature poor observational coverage compared to other centers of humid tropical forest change, such as Mato Grosso, Brazil, due to the lack of ongoing acquisitions from nearby ground stations and the persistence of cloud cover obscuring the land surface. At the same time, forest change in Indonesia is transient and does not always result in deforestation, as cleared forests are rapidly replaced by timber plantations and oil palm estates. Epochal composites, where single best observations are selected over a given time interval and used to quantify change, are one option for monitoring forest change in cloudy regions. However, the frequency of forest cover change in Indonesia confounds the ability of image composite pairs to quantify all change. Transient change occurring between composite periods is often missed and the length of time required for creating a cloud-free composite often obscures change occurring within the composite period itself. In this paper, we analyzed all Landsat 7 imagery with <50% cloud cover and data and products from the Moderate Resolution Imaging Spectroradiometer (MODIS) to quantify forest cover loss for Sumatra and Kalimantan from 2000 to 2005. We demonstrated that time-series approaches examining all good land observations are more accurate in mapping forest cover change in Indonesia than change maps based on image composites. Unlike other time-series analyses employing observations with a consistent periodicity, our study area was characterized by highly unequal observation counts and frequencies due to persistent cloud cover, scan line corrector off (SLC-off) gaps, and the absence of a complete archive. Our method accounts for this variation by generating a generic variable space. We evaluated our results against an independent probability sample-based estimate of gross forest cover loss and expert mapped gross forest cover loss at 64 sample sites. The mapped gross forest cover loss for Sumatra and Kalimantan was 2.86% of the land area, or 2.86 Mha from 2000 to 2005, with the highest concentration having occurred in Riau and Kalimantan Tengah provinces.

Reconciling forest conservation and logging in Indonesian Borneo.

Combining protected areas with natural forest timber concessions may sustain larger forest landscapes than is possible via protected areas alone. However, the role of timber concessions in maintaining natural forest remains poorly characterized. An estimated 57% (303,525 km2) of Kalimantans land area (532,100 km2) was covered by natural forest in 2000. About 14,212 km2 (4.7%) had been cleared by 2010. Forests in oil palm concessions had been reduced by 5,600 km2 (14.1%), while the figures for timber concessions are 1,336 km2 (1.5%), and for protected forests are 1,122 km2 (1.2%). These deforestation rates explain little about the relative performance of the different land use categories under equivalent conversion risks due to the confounding effects of location. An estimated 25% of lands allocated for timber harvesting in 2000 had their status changed to industrial plantation concessions in 2010. Based on a sample of 3,391 forest plots (1×1 km; 100 ha), and matching statistical analyses, 2000–2010 deforestation was on average 17.6 ha lower (95% C.I.: −22.3 ha–−12.9 ha) in timber concession plots than in oil palm concession plots. When location effects were accounted for, deforestation rates in timber concessions and protected areas were not significantly different (Mean difference: 0.35 ha; 95% C.I.: −0.002 ha–0.7 ha). Natural forest timber concessions in Kalimantan had similar ability as protected areas to maintain forest cover during 2000–2010, provided the former were not reclassified to industrial plantation concessions. Our study indicates the desirability of the Government of Indonesia designating its natural forest timber concessions as protected areas under the IUCN Protected Area Category VI to protect them from reclassification.

Remotely sensed forest cover loss shows high spatial and temporal variation across Sumatera and Kalimantan, Indonesia 2000-2008

The Indonesian islands of Sumatera and Kalimantan (the Indonesian part of the island of Borneo) are a center of significant and rapid forest cover loss in the humid tropics with implications for carbon dynamics, biodiversity conservation, and local livelihoods. The aim of our research was to analyze and interpret annual trends of forest cover loss for different sub-regions of the study area. We mapped forest cover loss for 2000?2008 using multi-resolution remote sensing data from the Landsat enhanced thematic mapper plus (ETM +) and moderate resolution imaging spectroradiometer (MODIS) sensors and analyzed annual trends per island, province, and official land allocation zone. The total forest cover loss for Sumatera and Kalimantan 2000?2008 was 5.39?Mha, which represents 5.3% of the land area and 9.2% of the year 2000 forest cover of these two islands. At least 6.5% of all mapped forest cover loss occurred in land allocation zones prohibiting clearing. An additional 13.6% of forest cover loss occurred where clearing is legally restricted. The overall trend of forest cover loss increased until 2006 and decreased thereafter. The trends for Sumatera and Kalimantan were distinctly different, driven primarily by the trends of Riau and Central Kalimantan provinces, respectively. This analysis shows that annual mapping of forest cover change yields a clearer picture than a one-time overall national estimate. Monitoring forest dynamics is important for national policy makers, especially given the commitment of Indonesia to reducing greenhouse gas emissions as part of the reducing emissions from deforestation and forest degradation in developing countries initiative (REDD +). The improved spatio-temporal detail of forest change monitoring products will make it possible to target policies and projects in meeting this commitment. Accurate, annual forest cover loss maps will be integral to many REDD + objectives, including policy formulation, definition of baselines, detection of displacement, and the evaluation of the permanence of emission reduction.

A Production Efficiency Model-Based Method for Satellite Estimates of Corn and Soybean Yields in the Midwestern US

Remote sensing techniques that provide synoptic and repetitive observations over large geographic areas have become increasingly important in studying the role of agriculture in global carbon cycles. However, it is still challenging to model crop yields based on remotely sensed data due to the variation in radiation use efficiency (RUE) across crop types and the effects of spatial heterogeneity. In this paper, we propose a production efficiency model-based method to estimate corn and soybean yields with MODerate Resolution Imaging Spectroradiometer (MODIS) data by explicitly handling the following two issues: (1) field-measured RUE values for corn and soybean are applied to relatively pure pixels instead of the biome-wide RUE value prescribed in the MODIS vegetation productivity product (MOD17); and (2) contributions to productivity from vegetation other than crops in mixed pixels are deducted at the level of MODIS resolution. Our estimated yields statistically correlate with the national survey data for rainfed counties in the Midwestern US with low errors for both corn (R2 = 0.77; RMSE = 0.89 MT/ha) and soybeans (R2 = 0.66; RMSE = 0.38 MT/ha). Because the proposed algorithm does not require any retrospective analysis that constructs empirical relationships between the reported yields and remotely sensed data, it could monitor crop yields over large areas.

Surface water network structure, landscape resistance to movement and flooding vital for maintaining ecological connectivity across Australia’s largest river basin

ContextLandscape-scale research quantifying ecological connectivity is required to maintain the viability of populations in dynamic environments increasingly impacted by anthropogenic modification and environmental change.ObjectiveTo evaluate how surface water network structure, landscape resistance to movement, and flooding affect the connectivity of amphibian habitats within the Murray–Darling Basin (MDB), a highly modified but ecologically significant region of south-eastern Australia.MethodsWe evaluated potential connectivity network graphs based on circuit theory, Euclidean and least-cost path distances for two amphibian species with different dispersal abilities, and used graph theory metrics to compare regional- and patch-scale connectivity across a range of flooding scenarios.ResultsCircuit theory graphs were more connected than Euclidean and least-cost equivalents in floodplain environments, and less connected in highly modified or semi-arid regions. Habitat networks were highly fragmented for both species, with flooding playing a crucial role in facilitating landscape-scale connectivity. Both formally and informally protected habitats were more likely to form important connectivity “hubs” or “stepping stones” compared to non-protected habitats, and increased in importance with flooding.ConclusionsSurface water network structure and the quality of the intervening landscape matrix combine to affect the connectivity of MDB amphibian habitats in ways which vary spatially and in response to flooding. Our findings highlight the importance of utilising organism-relevant connectivity models which incorporate landscape resistance to movement, and accounting for dynamic landscape-scale processes such as flooding when quantifying connectivity to inform the conservation of dynamic and highly modified environments.

Spatiotemporal dynamics of surface water networks across a global biodiversity hotspot—implications for conservation

The concept of habitat networks represents an important tool for landscape conservation and management at regional scales. Previous studies simulated degradation of temporally fixed networks but few quantified the change in network connectivity from disintegration of key features that undergo naturally occurring spatiotemporal dynamics. This is particularly of concern for aquatic systems, which typically show high natural spatiotemporal variability. Here we focused on the Swan Coastal Plain, a bioregion that encompasses a global biodiversity hotspot in Australia with over 1500 water bodies of high biodiversity. Using graph theory, we conducted a temporal analysis of water body connectivity over 13 years of variable climate. We derived large networks of surface water bodies using Landsat data (1999–2011). We generated an ensemble of 278 potential networks at three dispersal distances approximating the maximum dispersal distance of different water dependent organisms. We assessed network connectivity through several network topology metrics and quantified the resilience of the network topology during wet and dry phases. We identified ‘stepping stone’ water bodies across time and compared our networks with theoretical network models with known properties. Results showed a highly dynamic seasonal pattern of variability in network topology metrics. A decline in connectivity over the 13 years was noted with potential negative consequences for species with limited dispersal capacity. The networks described here resemble theoretical scale-free models, also known as ‘rich get richer’ algorithm. The ‘stepping stone’ water bodies are located in the area around the Peel-Harvey Estuary, a Ramsar listed site, and some are located in a national park. Our results describe a powerful approach that can be implemented when assessing the connectivity for a particular organism with known dispersal distance. The approach of identifying the surface water bodies that act as ‘stepping stone’ over time may help prioritize surface water bodies that are essential for maintaining regional scale connectivity.

Indices of Vegetation Activity

In this chapter we explain satellite-based vegetation indices (VIs) as dynamic spectral measures of vegetation activity. VIs are among the most widely used satellite products in monitoring ecosystems and agriculture, resource management, and estimations of many biophysical canopy properties. A theoretical basis for their formulation is presented and we describe how VIs are processed and composited from satellite imagery. Recent trends in their validation and quality assessment using in situ tower measurements are also discussed. Finally, a cross section of major findings involving the use of satellite VIs in ecological and climate science is presented and we conclude with research challenges and environmental issues that will drive future uses of satellite VIs.

Modeling multidecadal surface water inundation dynamics and key drivers on large river basin scale using multiple time series of Earth‐observation and river flow data

Periodically inundated floodplain areas are hot spots of biodiversity and provide a broad range of ecosystem services but have suffered alarming declines in recent history. Despite their importance, their long-term surface water (SW) dynamics and hydroclimatic drivers remain poorly quantified on continental scales. In this study, we used a 26 year time series of Landsat-derived SW maps in combination with river flow data from 68 gauges and spatial time series of rainfall, evapotranspiration and soil moisture to statistically model SW dynamics as a function of key drivers across Australias Murray-Darling Basin (∼1 million km2). We fitted generalized additive models for 18,521 individual modeling units made up of 10 × 10 km grid cells, each split into floodplain, floodplain-lake, and nonfloodplain area. Average goodness of fit of models was high across floodplains and floodplain-lakes (r2 > 0.65), which were primarily driven by river flow, and was lower for nonfloodplain areas (r2 > 0.24), which were primarily driven by rainfall. Local climate conditions were more relevant for SW dynamics in the northern compared to the southern basin and had the highest influence in the least regulated and most extended floodplains. We further applied the models of two contrasting floodplain areas to predict SW extents of cloud-affected time steps in the Landsat series during the large 2010 floods with high validated accuracy (r2 > 0.97). Our framework is applicable to other complex river basins across the world and enables a more detailed quantification of large floods and drivers of SW dynamics compared to existing methods.

Patterns of tree-cover loss along the Indonesia–Malaysia border on Borneo

Borneo Island is experiencing rapid tree-cover loss. This loss has been quantified for the Indonesian part of the island at Landsat spatial resolution, but no recent study exists that extends across the border into Malaysia. This research focused on quantifying patterns of tree-cover loss in the Indonesia–Malaysia border zone on Borneo. The methods used for quantifying 2000–2010 tree-cover loss within 20 km on either side of the border are an internally consistent mapping algorithm used on Landsat imagery and a local indicator of spatial autocorrelation to quantify the concentration of loss. Within the 20 km zone on either side of the border, tree-cover loss rates in lowlands were high in the two countries (19.8% and 14.4%, in Indonesia and Malaysia, respectively), but rates in the Malaysian uplands were an order of magnitude higher than in the Indonesian uplands (2.95% and 0.25%, respectively). Clusters of tree-cover loss in the Malaysian uplands were considerably larger than in the Indonesian uplands.

Evaluating static and dynamic landscape connectivity modelling using a 25-year remote sensing time series

ContextDespite calls for landscape connectivity research to account for spatiotemporal dynamics, studies have overwhelmingly evaluated the importance of habitats for connectivity at single or limited moments in time. Remote sensing time series represent a promising resource for studying connectivity within dynamic ecosystems. However, there is a critical need to assess how static and dynamic landscape connectivity modelling approaches compare for prioritising habitats for conservation within dynamic environments.ObjectivesTo assess whether static landscape connectivity analyses can identify similar important areas for connectivity as analyses based on dynamic remotely sensed time series data.MethodsWe compared degree and betweenness centrality graph theory metric distributions from four static scenarios against equivalent results from a dynamic 25-year remotely sensed surface-water time series. Metrics were compared at multiple spatial aggregation scales across south-eastern Australia’s 1 million km2 semi-arid Murray–Darling Basin and three sub-regions with varying levels of hydroclimatic variability and development.ResultsWe revealed large differences between static and dynamic connectivity metric distributions that varied by static scenario, region, spatial scale and hydroclimatic conditions. Static and dynamic metrics showed particularly low overlap within unregulated and spatiotemporally variable regions, although similarities increased at coarse aggregation scales.ConclusionsIn regions that exhibit high spatiotemporal variability, static connectivity modelling approaches are unlikely to serve as effective surrogates for more data intensive approaches based on dynamic, remotely sensed data. Although this limitation may be moderated by spatially aggregating static connectivity outputs, our results highlight the value of remotely sensed time series for assessing connectivity in dynamic landscapes.