Mapping forests in monsoon Asia with ALOS PALSAR 50-m mosaic images and MODIS imagery in 2010

An Adaptive-Parameter Pixel Unmixing Method for Mapping Evergreen Forest Fractions Based on Time-Series NDVI: A Case Study of Southern China

Spectral unmixing remains the most popular method for estimating the composition of mixed pixels. However, the spectral-based unmixing method cannot easily distinguish vegetation with similar spectral characteristics (e.g., different forest tree species). Furthermore, in large areas with significant heterogeneity, extracting a large number of pure endmember samples is challenging. Here, we implement a fractional evergreen forest cover-self-adaptive parameter (FEVC-SAP) approach to measure FEVC at the regional scale from continuous intra-year time-series normalized difference vegetation index (NDVI) values derived from moderate resolution imaging spectroradiometer (MODIS) imagery acquired over southern China, an area with a complex mixture of temperate, subtropical, and tropical climates containing evergreen and deciduous forests. Considering the cover of evergreen forest as a fraction of total forest (evergreen forest plus non-evergreen forest), the dimidiate pixel model combined with an index of evergreen forest phenological characteristics (NDVIann-min: intra-annual minimum NDVI value) was used to distinguish between evergreen and non-evergreen forests within a pixel. Due to spatial heterogeneity, the optimal model parameters differ among regions. By dividing the study area into grids, our method converts image spectral information into gray level information and uses the Otsu threshold segmentation method to simulate the appropriate parameters for each grid for adaptive acquisition of FEVC parameters. Mapping accuracy was assessed at the pixel and sub-pixel scales. At the pixel scale, a confusion matrix was constructed with higher overall accuracy (87.5%) of evergreen forest classification than existing land cover products, including GLC 30 and MOD12. At the sub-pixel scale, a strong linear correlation was found between the cover fraction predicted by our method and the reference cover fraction obtained from GF-1 images (R2 = 0.86). Compared to other methods, the FEVC-SAP had a lower estimation deviation (root mean square error = 8.6%). Moreover, the proposed method had greater estimation accuracy in densely than sparsely forested areas. Our results highlight the utility of the adaptive-parameter linear unmixing model for quantitative evaluation of the coverage of evergreen forest and other vegetation types at large scales.

Comparative Analysis of the Global Forest/Non-Forest Maps Derived from SAR and Optical Sensors. Case Studies from Brazilian Amazon and Cerrado Biomes

Global-scale forest/non-forest (FNF) maps are of crucial importance for applications like biomass estimation and deforestation monitoring. Global FNF maps based on optical remote sensing data have been produced by the wall-to-wall satellite image analyses or sampling strategies. The German Aerospace Center (DLR) and the Japan Aerospace Exploration Agency (JAXA) also made available their global FNF maps based on synthetic aperture radar (SAR) data. This paper attempted to answer the following scientific question: how comparable are the FNF products derived from optical and SAR data? As test sites we selected the Amazon (tropical rainforest) and Cerrado (tropical savanna) biomes, the two largest Brazilian biomes. Forest estimations from 2015 derived from TanDEM-X (X band; HH polarization) and ALOS-2 (L band; HV polarization) SAR data, as well as forest cover information derived from Landsat 8 optical data were compared with each other at the municipality and image sampling levels. The optical-based forest estimations considered in this study were derived from the MapBiomas project, a Brazilian multi-institutional project to map land use and land cover (LULC) classes of an entire country based on historical time series of Landsat data. In addition to the existing forest maps, a set of 1619 Landsat 8 RGB color composites was used to generate new independent comparison data composed of circular areas with 5-km diameter, which were visually interpreted after image segmentation. The Spearman rank correlation estimated the correlation among the data sets and the paired Mann–Whitney–Wilcoxon tested the hypothesis that the data sets are statistically equal. Results showed that forest maps derived from SAR and optical satellites are statistically different regardless of biome or scale of study (municipality or image sampling), except for the Cerrado´s forest estimations derived from TanDEM-X and ALOS-2. Nevertheless, the percentage of pixels classified as forest or non-forest by both SAR sensors were 90% and 80% for the Amazon and Cerrado biome, respectively, indicating an overall good agreement.

New JAXA High-Resolution Land Use/Land Cover Map for Vietnam Aiming for Natural Forest and Plantation Forest Monitoring

Highly detailed and accurate forest maps are important for various applications including forest monitoring, forestry policy, climate change, and biodiversity loss. This study demonstrates a comprehensive and geographically transferable approach to produce a 12 category high-resolution land use/land cover (LULC) map over mainland Vietnam in 2016 by remote sensing data. The map included several natural forest categories (evergreen broadleaf, deciduous (mostly deciduous broadleaf), and coniferous (mostly evergreen coniferous)) and one category representing all popular plantation forests in Vietnam such as acacia (Acacia mangium, Acacia auriculiformis, Acacia hybrid), eucalyptus (Eucalyptus globulus), rubber (Hevea brasiliensis), and others. The approach combined the advantages of various sensor data by integrating their posterior probabilities resulting from applying a probabilistic classifier (comprised of kernel density estimation and Bayesian inference) to each datum individually. By using different synthetic aperture radar (SAR) images (PALSAR-2/ScanSAR, PALSAR-2 mosaic, Sentinel-1), optical images (Sentinel-2, Landsat-8) and topography data (AW3D30), the resultant map achieved 85.6% for the overall accuracy. The major forest classes including evergreen broadleaf forests and plantation forests had a user’s accuracy and producer’s accuracy ranging from 86.0% to 95.3%. Our map identified 9.55 × 106 ha (±0.16 × 106 ha) of natural forests and 3.89 × 106 ha (±0.11 × 106 ha) of plantation forests over mainland Vietnam, which were close to the Vietnamese government’s statistics (with differences of less than 8%). This study’s result provides a reliable input/reference to support forestry policy and land sciences in Vietnam.

Spaceborne L-Band Synthetic Aperture Radar Data for Geoscientific Analyses in Coastal Land Applications: A Review

The coastal zone offers among the world’s most productive and valuable ecosystems and is experiencing increasing pressure from anthropogenic impacts: human settlements, agriculture, aquaculture, trade, industrial activities, oil and gas exploitation and tourism. Earth observation has great capability to deliver valuable data at the local, regional and global scales and can support the assessment and monitoring of land- and water-related applications in coastal zones. Compared to optical satellites, cloud-cover does not limit the timeliness of data acquisition with spaceborne Synthetic Aperture Radar (SAR) sensors, which have all-weather, day and night capabilities. Hence, active radar systems demonstrate great potential for continuous mapping and monitoring of coastal regions, particularly in cloud-prone tropical and sub-tropical climates. The canopy penetration capability with long radar wavelength enables L-band SAR data to be used for coastal terrestrial environments and has been widely applied and investigated for the following geoscientific topics: mapping and monitoring of flooded vegetation and inundated areas; the retrieval of aboveground biomass; and the estimation of soil moisture. Human activities, global population growth, urban sprawl and climate change-induced impacts are leading to increased pressure on coastal ecosystems causing land degradation, deforestation and land use change. This review presents a comprehensive overview of existing research articles that apply spaceborne L-band SAR data for geoscientific analyses that are relevant for coastal land applications.

Modelling Agriculture, Forestry and Other Land Use (AFOLU) in response to climate change scenarios for the SAARC nations

Agriculture and forestry are the two major land use classes providing sustenance to the human population. With the pace of development, these two land use classes continue to change over time. Land use change is a dynamic process under the influence of multiple drivers including climate change. Therefore, tracing the trajectory of the changes is challenging. The artificial neural network (ANN) has successfully been applied for tracing such a dynamic process to capture nonlinear responses. We test the application of the multilayer perceptron neural network (MLP-NN) to project the future Agriculture, Forestry and Other Land Use (AFOLU) for the year 2050 for the South Asian Association for Regional Cooperation (SAARC) nations which is a geopolitical union of Afghanistan, Bangladesh, Bhutan, India, Nepal, Maldives, Pakistan and Sri Lanka. The Intergovernmental Panel on Climate Change (IPCC) and Food and Agriculture Organization (FAO) use much frequently the term 'AFOLU' in their policy documents. Hence, we restricted our land use classification scheme as AFOLU for assessing the influence of climate change scenarios of the IPCC fifth assessment report (RCP 2.6, RCP 4.5, RCP 6.0 and RCP 8.5). Agricultural land would increase in all the SAARC nations, with the highest increase in Pakistan and Maldives; moderate increase in Afghanistan, India and Nepal; and the least increase in Bangladesh, Bhutan and Sri Lanka. The forestry land use will witness a decreasing trend under all scenarios in all of the SAARC nations with varying levels of changes. The study is expected to assist planners and policymakers to develop nations' specific strategy to proportionate land use classes to meet various needs on a sustainable basis.

Tracking Reforestation in the Loess Plateau, China after the “Grain for Green” Project through Integrating PALSAR and Landsat Imagery

An unprecedented reforestation process happened in the Loess Plateau, China due to the ecological restoration project ‘Grain for Green Project’, which has affected regional carbon and water cycles as well as brought climate feedbacks. Accurately mapping the area and spatial distribution of emerged forests in the Loess Plateau over time is essential for forest management but a very challenging task. Here we investigated the changes of forests in the Loess Plateau after the forest reconstruction project. First, we used a pixel and rule-based algorithm to identify and map the annual forests from 2007 to 2017 in the Loess Plateau by integrating 30 m Landsat data and 25 m resolution PALSAR data in this study. Then, we carried out the accuracy assessment and comparison with several existing forest products. The overall accuracy (OA) and Kappa coefficient of the resultant map, were about 91% and 0.77 in 2010, higher than those of the other forest products (FROM-GLC, GlobeLand30, GLCF-VCF, JAXA, and OU-FDL) with OA ranging from 83.57% to 87.96% and Kappa coefficients from 0.52 to 0.68. Based on the annual forest maps, we found forest area in the Loess Plateau has increased by around 15,000 km2 from 2007 to 2017. This study clearly demonstrates the advantages of data fusion between PALSAR and Landsat images for monitoring forest cover dynamics in the Loess Plateau, and the resultant forest maps with lower uncertainty would contribute to the regional forest management.

JAXA Annual Forest Cover Maps for Vietnam during 2015–2018 Using ALOS-2/PALSAR-2 and Auxiliary Data

Monitoring the temporal changes of forests is important for sustainable forest management. In this study, we investigated the potential of using multi-temporal synthetic aperture radar (SAR) images for mapping annual change in forest cover at a national scale. We assessed the robustness of using multi-temporal Phased Array L-band Synthetic Aperture Radar-2/Scanning Synthetic Aperture Radar (PALSAR-2/ScanSAR) mosaic images for forest mapping by comparison with single-temporal PALSAR-2 mosaic images for three test sites in North, Central, and Southern Vietnam. We then used a combination of multi-temporal PALSAR-2/ScanSAR images, multi-temporal Moderate Resolution Imaging Spectroradiometer (MODIS) Normalized Difference Vegetation Index (NDVI) images, and Shuttle Radar Topography Mission (SRTM) images to map annual forest cover for mainland Vietnam during 2015–2018. Average overall accuracies of our forest/non-forest (FNF) maps (86.6% ± 3.1%) were greater than recent maps of Japan Aerospace Exploration Agency (JAXA, (77.5% ± 3.2%)) and European Space Agency (ESA, (85.4% ± 1.6%)). Our estimates of mainland Vietnam’s forest area were close to that of the Vietnamese government. A comparison of the spatial distribution of forest estimated from JAXA and ESA FNF maps showed that our FNF map in 2015 agreed relatively well with the ESA map, with 77% of pixels being consistent. This study demonstrates the merit of using multi-temporal PALSAR-2/ScanSAR images for annual forest mapping at a national scale.

Capturing Coastal Dune Natural Vegetation Types Using a Phenology-Based Mapping Approach: The Potential of Sentinel-2

Coastal areas harbor the most threatened ecosystems on Earth, and cost-effective ways to monitor and protect them are urgently needed, but they represent a challenge for habitat mapping and multi-temporal observations. The availability of open access, remotely sensed data with increasing spatial and spectral resolution is promising in this context. Thus, in a sector of the Mediterranean coast (Lazio region, Italy), we tested the strength of a phenology-based vegetation mapping approach and statistically compared results with previous studies, making use of open source products across all the processing chain. We identified five accurate land cover classes in three hierarchical levels, with good values of agreement with previous studies for the first and the second hierarchical level. The implemented procedure resulted as being effective for mapping a highly fragmented coastal dune system. This is encouraging to take advantage of the earth observation through remote sensing technology in an open source perspective, even at the fine scale of highly fragmented sand dunes landscapes.

Mapping Annual Forest Change Due to Afforestation in Guangdong Province of China Using Active and Passive Remote Sensing Data

Accurate acquisition of spatial distribution of afforestation in a large area is of great significance to contributing to the sustainable utilization of forest resources and the evaluation of the carbon accounting. Annual forest maps (1986–2016) of Guangdong, China were generated using time series Landsat images and PALSAR data. Initially, four PALSAR-based classifiers were used to classify land cover types. Then, the optimal mapping algorithm was determined. Next, an accurate identification of forest and non-forest was carried out by combining Landsat-based phenological variables and PALSAR-based land cover classifications. Finally, the spatio-temporal distribution of forest cover change due to afforestation was created and its forest biomass dynamics changes were detected. The results indicated that the overall accuracy of forest classification of the improved model based on the PALSAR-based stochastic gradient boosting (SGB) classification and the maximum value of normalized difference vegetation index (NDVI; SGB-NDVI) were approximately 75–85% in 2005, 2010, and 2016. Compared with the Japan Aerospace Exploration Agency (JAXA) PALSAR-forest/non-forest, the SGB-NDVI-based forest product showed great improvement, while the SGB-NDVI product was the same or slightly inferior to the Global Land Cover (GLC) and vegetation tracker change (VCT)-based land cover types, respectively. Although this combination of multiple sources contained some errors, the SGB-NDVI model effectively identified the distribution of forest cover changes by afforestation events. By integrating aboveground biomass dynamics (AGB) change with forest cover, the trend in afforestation area closely corresponded with the trend in forest AGB. This technique can provide an essential data baseline for carbon assessment in the planted forests of southern China.

The NDVI-CV Method for Mapping Evergreen Trees in Complex Urban Areas Using Reconstructed Landsat 8 Time-Series Data

Evergreen trees play a significant role in urban ecological services, such as air purification, carbon and oxygen balance, and temperature and moisture regulation. Remote sensing represents an essential technology for obtaining spatiotemporal distribution data for evergreen trees in cities. However, highly developed subtropical cities, such as Nanjing, China, have serious land fragmentation problems, which greatly increase the difficulty of extracting evergreen trees information and reduce the extraction precision of remote-sensing methods. This paper introduces a normalized difference vegetation index coefficient of variation (NDVI-CV) method to extract evergreen trees from remote-sensing data by combining the annual minimum normalized difference vegetation index (NDVIann-min) with the CV of a Landsat 8 time-series NDVI. To obtain an intra-annual, high-resolution time-series dataset, Landsat 8 cloud-free and partially cloud-free images over a three-year period were collected and reconstructed for the study area. Considering that the characteristic growth of evergreen trees remained nearly unchanged during the phenology cycle, NDVIann-min is the optimal phenological node to separate this information from that of other vegetation types. Furthermore, the CV of time-series NDVI considers all of the phenologically critical phases; therefore, the NDVI-CV method had higher extraction accuracy. As such, the approach presented herein represents a more practical and promising method based on reasonable NDVIann-min and CV thresholds to obtain spatial distribution data for evergreen trees. The experimental verification results indicated a comparable performance since the extraction accuracy of the model was over 85%, which met the classification accuracy requirements. In a cross-validation comparison with other evergreen trees’ extraction methods, the NDVI-CV method showed higher sensitivity and stability.

On the Synergistic Use of Optical and SAR Time-Series Satellite Data for Small Mammal Disease Host Mapping

(1) Background: Echinococcus multilocularis (Em), a highly pathogenic parasitic tapeworm, is responsible for a significant burden of human disease. In this study, optical and time-series Synthetic Aperture Radar (SAR) data is used synergistically to model key land cover characteristics driving the spatial distributions of two small mammal intermediate host species, Ellobius tancrei and Microtus gregalis, which facilitate Em transmission in a highly endemic area of Kyrgyzstan. (2) Methods: A series of land cover maps are derived from (a) single-date Landsat Operational Land Imager (OLI) imagery, (b) time-series Sentinel-1 SAR data, and (c) Landsat OLI and time-series Sentinel-1 SAR data in combination. Small mammal distributions are analyzed in relation to the surrounding land cover class coverage using random forests, before being applied predictively over broader areas. A comparison of models derived from the three land cover maps are made, assessing their potential for use in cloud-prone areas. (3) Results: Classification accuracies demonstrated the combined OLI-SAR classification to be of highest accuracy, with the single-date OLI and time-series SAR derived classifications of equivalent quality. Random forest analysis identified statistically significant positive relationships between E. tancrei density and agricultural land, and between M. gregalis density and water and bushes. Predictive application of random forest models identified hotspots of high relative density of E. tancrei and M. gregalis across the broader study area. (4) Conclusions: This offers valuable information to improve the targeting of limited-resource disease control activities to disrupt disease transmission in this area. Time-series SAR derived land cover maps are shown to be of equivalent quality to those generated from single-date optical imagery, which enables application of these methods in cloud-affected areas where, previously, this was not possible due to the sparsity of cloud-free optical imagery.

Developing Forest Cover Composites through a Combination of Landsat-8 Optical and Sentinel-1 SAR Data for the Visualization and Extraction of Forested Areas

Mapping the distribution of forested areas and monitoring their spatio-temporal changes are necessary for the conservation and management of forests. This paper presents two new image composites for the visualization and extraction of forest cover. By exploiting the Landsat-8 satellite-based multi-temporal and multi-spectral reflectance datasets, the Forest Cover Composite (FCC) was designed in this research. The FCC is an RGB (red, green, blue) color composite made up of short-wave infrared reflectance and green reflectance, specially selected from the day when the Normalized Difference Vegetation Index (NDVI) is at a maximum, as the red and blue bands, respectively. The annual mean NDVI values are used as the green band. The FCC is designed in such a way that the forested areas appear greener than other vegetation types, such as grasses and shrubs. On the other hand, the croplands and barren lands are usually seen as red and water/snow is seen as blue. However, forests may not necessarily be greener than other perennial vegetation. To cope with this problem, an Enhanced Forest Cover Composite (EFCC) was designed by combining the annual median backscattering intensity of the VH (vertical transmit, horizontal receive) polarization data from the Sentinel-1 satellite with the green term of the FCC to suppress the green component (mean NDVI values) of the FCC over the non-forested vegetative areas. The performances of the FCC and EFCC were evaluated for the discrimination and classification of forested areas all over Japan with the support of reference data. The FCC and EFCC provided promising results, and the high-resolution forest map newly produced in the research provided better accuracy than the extant MODIS (Moderate Resolution Imaging Spectroradiometer) Land Cover Type product (MCD12Q1) in Japan. The composite images proposed in the research are expected to improve forest monitoring activities in other regions as well.

Synergizing community-based forest monitoring with remote sensing: a path to an effective REDD+ MRV system

BACKGROUND

The reliable monitoring, reporting and verification (MRV) of carbon emissions and removals from the forest sector is an important part of the efforts on reducing emissions from deforestation and forest degradation (REDD+). Forest-dependent local communities are engaged to contribute to MRV through community-based monitoring systems. The efficiency of such monitoring systems could be improved through the rational integration of the studies at permanent plots with the geospatial technologies. This article presents a case study of integrating community-based measurements at permanent plots at the foothills of central Nepal and biomass maps that were developed using GeoEye-1 and IKONS satellite images.

RESULTS

The use of very-high-resolution satellite-based tree cover parameters, including crown projected area (CPA), crown density and crown size classes improves salience, reliability and legitimacy of the community-based survey of 0.04% intensity at the lower cost than increasing intensity of the community-based survey to 0.14% level (2.5 USD/ha vs. 7.5 USD/ha).

CONCLUSION

The proposed REDD+ MRV complementary system is the first of its kind and demonstrates the enhancement of information content, accuracy of reporting and reduction in cost. It also allows assessment of the efficacy of community-based forest management and extension to national scale.

Forest cover mapping in post-Soviet Central Asia using multi-resolution remote sensing imagery

Despite rapid advances and large-scale initiatives in forest mapping, reliable cross-border information about the status of forest resources in Central Asian countries is lacking. We produced consistent Central Asia forest cover (CAFC) maps based on a cost-efficient approach using multi-resolution satellite imagery from Landsat and MODIS during 2009–2011. The spectral-temporal metrics derived from 2009–2011 Landsat imagery (overall accuracy of 0.83) was used to predict sub-pixel forest cover on the MODIS scale for 2010. Accuracy assessment confirmed the validity of MODIS-based forest cover map with a normalized root-mean-square error of 0.63. A general paucity of forest resources in post-Soviet Central Asia was indicated, with 1.24% of the region covered by forest. In comparison to the CAFC map, a regional map derived from MODIS Vegetation Continuous Fields tended to underestimate forest cover, while the Global Forest Change product matched well. The Global Forest Resources Assessments, based on individual country reports, overestimated forest cover by 1.5 to 147 times, particularly in the more arid countries of Turkmenistan and Uzbekistan. Multi-resolution imagery contributes to regionalized assessment of forest cover in the world’s drylands while developed CAFC maps (available at https://data.zef.de/) aim to facilitate decisions on biodiversity conservation and reforestation programs in Central Asia.

Spatiotemporal patterns of paddy rice croplands in China and India from 2000 to 2015

Due to rapid population growth and urbanization, paddy rice agriculture is experiencing substantial changes in the spatiotemporal pattern of planting areas in the two most populous countries-China and India-where food security is always the primary concern. However, there is no spatially explicit and continuous rice-planting information in either country. This knowledge gap clearly hinders our ability to understand the effects of spatial paddy rice area dynamics on the environment, such as food and water security, climate change, and zoonotic infectious disease transmission. To resolve this problem, we first generated annual maps of paddy rice planting areas for both countries from 2000 to 2015, which are derived from time series Moderate Resolution Imaging Spectroradiometer (MODIS) data and the phenology- and pixel-based rice mapping platform (RICE-MODIS), and analyzed the spatiotemporal pattern of paddy rice dynamics in the two countries. We found that China experienced a general decrease in paddy rice planting area with a rate of 0.72 million (m) ha/yr from 2000 to 2015, while a significant increase at a rate of 0.27mha/yr for the same time period happened in India. The spatial pattern of paddy rice agriculture in China shifted northeastward significantly, due to simultaneous expansions in paddy rice planting areas in northeastern China and contractions in southern China. India showed an expansion of paddy rice areas across the entire country, particularly in the northwestern region of the Indo-Gangetic Plain located in north India and the central and south plateau of India. In general, there has been a northwesterly shift in the spatial pattern of paddy rice agriculture in India. These changes in the spatiotemporal patterns of paddy rice planting area have raised new concerns on how the shift may affect national food security and environmental issues relevant to water, climate, and biodiversity.

Mapping paddy rice planting area in northeastern Asia with Landsat 8 images, phenology-based algorithm and Google Earth Engine

Area and spatial distribution information of paddy rice are important for understanding of food security, water use, greenhouse gas emission, and disease transmission. Due to climatic warming and increasing food demand, paddy rice has been expanding rapidly in high latitude areas in the last decade, particularly in northeastern (NE) Asia. Current knowledge about paddy rice fields in these cold regions is limited. The phenology- and pixel-based paddy rice mapping (PPPM) algorithm, which identifies the flooding signals in the rice transplanting phase, has been effectively applied in tropical areas, but has not been tested at large scale of cold regions yet. Despite the effects from more snow/ice, paddy rice mapping in high latitude areas is assumed to be more encouraging due to less clouds, lower cropping intensity, and more observations from Landsat sidelaps. Moreover, the enhanced temporal and geographic coverage from Landsat 8 provides an opportunity to acquire phenology information and map paddy rice. This study evaluated the potential of Landsat 8 images on annual paddy rice mapping in NE Asia which was dominated by single cropping system, including Japan, North Korea, South Korea, and NE China. The cloud computing approach was used to process all the available Landsat 8 imagery in 2014 (143 path/rows, ~3290 scenes) with the Google Earth Engine (GEE) platform. The results indicated that the Landsat 8, GEE, and improved PPPM algorithm can effectively support the yearly mapping of paddy rice in NE Asia. The resultant paddy rice map has a high accuracy with the producer (user) accuracy of 73% (92%), based on the validation using very high resolution images and intensive field photos. Geographic characteristics of paddy rice distribution were analyzed from aspects of country, elevation, latitude, and climate. The resultant 30-m paddy rice map is expected to provide unprecedented details about the area, spatial distribution, and landscape pattern of paddy rice fields in NE Asia, which will contribute to food security assessment, water resource management, estimation of greenhouse gas emissions, and disease control.