Exacerbated grassland degradation and desertification in Central Asia during 2000-2014.

Grassland degradation and desertification is a complex process, including both state conversion (e.g., grasslands to deserts) and gradual within-state change (e.g., greenness dynamics). Existing studies hardly separated the two components and analyzed it as a whole based on time series vegetation index data, which cannot provide a clear and comprehensive picture for grassland degradation and desertification. Here we propose an integrated assessment strategy, by considering both state conversion and within-state change of grasslands, to investigate grassland degradation and desertification process in Central Asia. First, annual maps of grasslands and sparsely vegetated land were generated to track the state conversions between them. The results showed increasing grasslands were converted to sparsely vegetated lands from 2000 to 2014, with the desertification region concentrating in the latitude range of 43-48° N. A frequency analysis of grassland vs. sparsely vegetated land classification in the last 15 yr allowed a recognition of persistent desert zone (PDZ), persistent grassland zone (PGZ), and transitional zone (TZ). The TZ was identified in southern Kazakhstan as one hotspot that was unstable and vulnerable to desertification. Furthermore, the trend analysis of Enhanced Vegetation Index during thermal growing season (EVITGS ) was investigated in individual zones using linear regression and Mann-Kendall approaches. An overall degradation across the area was found; moreover, the second desertification hotspot was identified in northern Kazakhstan with significant decreasing in EVITGS , which was located in PGZ. Finally, attribution analyses of grassland degradation and desertification were conducted by considering precipitation, temperature, and three different drought indices. We found persistent droughts were the main factor for grassland degradation and desertification in Central Asia. Considering both state conversion and gradual within-state change processes, this study provided reference information for identification of desertification hotspots to support further grassland degradation and desertification treatment, and the method could be useful to be extended to other regions.

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.

Timely monitoring of Asian Migratory locust habitats in the Amudarya delta, Uzbekistan using time series of satellite remote sensing vegetation index.

The Asian Migratory locust (Locusta migratoria migratoria L.) is a pest that continuously threatens crops in the Amudarya River delta near the Aral Sea in Uzbekistan, Central Asia. Its development coincides with the growing period of its main food plant, a tall reed grass (Phragmites australis), which represents the predominant vegetation in the delta and which cover vast areas of the former Aral Sea, which is desiccating since the 1960s. Current locust survey methods and control practices would tremendously benefit from accurate and timely spatially explicit information on the potential locust habitat distribution. To that aim, satellite observation from the MODIS Terra/Aqua satellites and in-situ observations were combined to monitor potential locust habitats according to their corresponding risk of infestations along the growing season. A Random Forest (RF) algorithm was applied for classifying time series of MODIS enhanced vegetation index (EVI) from 2003 to 2014 at an 8-day interval. Based on an independent ground truth data set, classification accuracies of reeds posing a medium or high risk of locust infestation exceeded 89% on average. For the 12-year period covered in this study, an average of 7504 km2 (28% of the observed area) was flagged as potential locust habitat and 5% represents a permanent high risk of locust infestation. Results are instrumental for predicting potential locust outbreaks and developing well-targeted management plans. The method offers positive perspectives for locust management and treatment of infested sites because it is able to deliver risk maps in near real time, with an accuracy of 80% in April-May which coincides with both locust hatching and the first control surveys. Such maps could help in rapid decision-making regarding control interventions against the initial locust congregations, and thus the efficiency of survey teams and the chemical treatments could be increased, thus potentially reducing environmental pollution while avoiding areas where treatments are most likely to cause environmental degradation.

Spatial characterization of colonies of the flying fox bat, a carrier of Nipah virus in Thailand.

BACKGROUND: A major reservoir of Nipah virus is believed to be the flying fox genus Pteropus, a fruit bat distributed across many of the world's tropical and sub-tropical areas. The emergence of the virus and its zoonotic transmission to livestock and humans have been linked to losses in the bat's habitat. Nipah has been identified in a number of indigenous flying fox populations in Thailand. While no evidence of infection in domestic pigs or people has been found to date, pig farming is an active agricultural sector in Thailand and therefore could be a potential pathway for zoonotic disease transmission from the bat reservoirs. The disease, then, represents a potential zoonotic risk. To characterize the spatial habitat of flying fox populations along Thailand's Central Plain, and to map potential contact zones between flying fox habitats, pig farms and human settlements, we conducted field observation, remote sensing, and ecological niche modeling to characterize flying fox colonies and their ecological neighborhoods. A Potential Surface Analysis was applied to map contact zones among local epizootic actors. RESULTS: Flying fox colonies are found mainly on Thailand's Central Plain, particularly in locations surrounded by bodies of water, vegetation, and safe havens such as Buddhist temples. High-risk areas for Nipah zoonosis in pigs include the agricultural ring around the Bangkok metropolitan region where the density of pig farms is high. CONCLUSIONS: Passive and active surveillance programs should be prioritized around Bangkok, particularly on farms with low biosecurity, close to water, and/or on which orchards are concomitantly grown. Integration of human and animal health surveillance should be pursued in these same areas. Such proactive planning would help conserve flying fox colonies and should help prevent zoonotic transmission of Nipah and other pathogens.

Mapping paddy rice planting areas through time series analysis of MODIS land surface temperature and vegetation index data.

Knowledge of the area and spatial distribution of paddy rice is important for assessment of food security, management of water resources, and estimation of greenhouse gas (methane) emissions. Paddy rice agriculture has expanded rapidly in northeastern China in the last decade, but there are no updated maps of paddy rice fields in the region. Existing algorithms for identifying paddy rice fields are based on the unique physical features of paddy rice during the flooding and transplanting phases and use vegetation indices that are sensitive to the dynamics of the canopy and surface water content. However, the flooding phenomena in high latitude area could also be from spring snowmelt flooding. We used land surface temperature (LST) data from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor to determine the temporal window of flooding and rice transplantation over a year to improve the existing phenology-based approach. Other land cover types (e.g., evergreen vegetation, permanent water bodies, and sparse vegetation) with potential influences on paddy rice identification were removed (masked out) due to their different temporal profiles. The accuracy assessment using high-resolution images showed that the resultant MODIS-derived paddy rice map of northeastern China in 2010 had a high accuracy (producer and user accuracies of 92% and 96%, respectively). The MODIS-based map also had a comparable accuracy to the 2010 Landsat-based National Land Cover Dataset (NLCD) of China in terms of both area and spatial pattern. This study demonstrated that our improved algorithm by using both thermal and optical MODIS data, provides a robust, simple and automated approach to identify and map paddy rice fields in temperate and cold temperate zones, the northern frontier of rice planting.

Water Productivity Mapping (WPM) Using Landsat ETM+ Data for the Irrigated Croplands of the Syrdarya River Basin in Central Asia.

The overarching goal of this paper was to espouse methods and protocols for water productivity mapping (WPM) using high spatial resolution Landsat remote sensing data. In a world where land and water for agriculture are becoming increasingly scarce, growing "more crop per drop" (increasing water productivity) becomes crucial for food security of future generations. The study used time-series Landsat ETM+ data to produce WPMs of irrigated crops, with emphasis on cotton in the Galaba study area in the Syrdarya river basin of Central Asia. The WPM methods and protocols using remote sensing data consisted of: (1) crop productivity (ton/ha) maps (CPMs) involvingcrop type classification, crop yield and biophysical modeling, and extrapolating yield models to larger areas using remotely sensed data; (2) crop water use (m³/ha) maps (WUMs) (or actual seasonal evapotranspiration or actual ET) developed through Simplified Surface Energy Balance (SSEB) model; and (3) water productivity (kg/m³) maps (WPMs) produced by dividing raster layers of CPMs by WUMs. The SSEB model calculated WUMs (actual ET) by multiplying the ET fractionby reference ET. The ETfraction was determined using Landsat thermal imagery by selecting the "hot" pixels (zero ET) and "cold" pixels (maximum ET). The grass reference ET was calculated by FAO Penman-Monteith method using meteorological data. The WPMs for the Galaba study area demonstrated a wide variations (0-0.54 kg/m³) in water productivity of cotton fields with overwhelming proportion (87%) of the area having WP less than 0.30 kg/m³, 11% of the area having WP in range of 0.30-0.36 kg/m³, and only 2% of the area with WP greater than 0.36 kg/m³. These results clearly imply that there are opportunities for significant WP increases in overwhelming proportion of the existing croplands. The areas of low WP are spatially pin-pointed and can be used as focus for WP improvements through better land and water management practices.

Energy determines broad pattern of plant distribution in Western Himalaya.

Several factors describe the broad pattern of diversity in plant species distribution. We explore these determinants of species richness in Western Himalayas using high-resolution species data available for the area to energy, water, physiography and anthropogenic disturbance. The floral data involves 1279 species from 1178 spatial locations and 738 sample plots of a national database. We evaluated their correlation with 8-environmental variables, selected on the basis of correlation coefficients and principal component loadings, using both linear (structural equation model) and nonlinear (generalised additive model) techniques. There were 645 genera and 176 families including 815 herbs, 213 shrubs, 190 trees, and 61 lianas. The nonlinear model explained the maximum deviance of 67.4% and showed the dominant contribution of climate on species richness with a 59% share. Energy variables (potential evapotranspiration and temperature seasonality) explained the deviance better than did water variables (aridity index and precipitation of the driest quarter). Temperature seasonality had the maximum impact on the species richness. The structural equation model confirmed the results of the nonlinear model but less efficiently. The mutual influences of the climatic variables were found to affect the predictions of the model significantly. To our knowledge, the 67.4% deviance found in the species richness pattern is one of the highest values reported in mountain studies. Broadly, climate described by water-energy dynamics provides the best explanation for the species richness pattern. Both modeling approaches supported the same conclusion that energy is the best predictor of species richness. The dry and cold conditions of the region account for the dominant contribution of energy on species richness.

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.

Priority regions for research on dryland cereals and legumes.

Dryland cereals and legumes  are important crops in farming systems across the world.  Yet they are frequently neglected among the priorities for international agricultural research and development, often due to lack of information on their magnitude and extent. Given what we know about the global distribution of dryland cereals and legumes, what regions should be high priority for research and development to improve livelihoods and food security? This research evaluated the geographic dimensions of these crops and the farming systems where they are found worldwide. The study employed geographic information science and data to assess the key farming systems and regions for these crops. Dryland cereal and legume crops should be given high priority in 18 farming systems worldwide, where their cultivated area comprises more than 160 million ha. These regions include the dryer areas of South Asia, West and East Africa, the Middle East and North Africa, Central America and other parts of Asia. These regions are prone to drought and heat stress, have limiting soil constraints, make up half of the global population and account for 60 percent of the global poor and malnourished. The dryland cereal and legume crops and farming systems merit more research and development attention to improve productivity and address development problems. This project developed an open access dataset and information resource that provides the basis for future analysis of the geographic dimensions of dryland cereals and legumes.

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

Extensive forest changes have occurred in monsoon Asia, substantially affecting climate, carbon cycle and biodiversity. Accurate forest cover maps at fine spatial resolutions are required to qualify and quantify these effects. In this study, an algorithm was developed to map forests in 2010, with the use of structure and biomass information from the Advanced Land Observation System (ALOS) Phased Array L-band Synthetic Aperture Radar (PALSAR) mosaic dataset and the phenological information from MODerate Resolution Imaging Spectroradiometer (MOD13Q1 and MOD09A1) products. Our forest map (PALSARMOD50 m F/NF) was assessed through randomly selected ground truth samples from high spatial resolution images and had an overall accuracy of 95%. Total area of forests in monsoon Asia in 2010 was estimated to be ~6.3 × 10(6 )km(2). The distribution of evergreen and deciduous forests agreed reasonably well with the median Normalized Difference Vegetation Index (NDVI) in winter. PALSARMOD50 m F/NF map showed good spatial and areal agreements with selected forest maps generated by the Japan Aerospace Exploration Agency (JAXA F/NF), European Space Agency (ESA F/NF), Boston University (MCD12Q1 F/NF), Food and Agricultural Organization (FAO FRA), and University of Maryland (Landsat forests), but relatively large differences and uncertainties in tropical forests and evergreen and deciduous forests.

A 50-m forest cover map in Southeast Asia from ALOS/PALSAR and its application on forest fragmentation assessment.

Southeast Asia experienced higher rates of deforestation than other continents in the 1990s and still was a hotspot of forest change in the 2000s. Biodiversity conservation planning and accurate estimation of forest carbon fluxes and pools need more accurate information about forest area, spatial distribution and fragmentation. However, the recent forest maps of Southeast Asia were generated from optical images at spatial resolutions of several hundreds of meters, and they do not capture well the exceptionally complex and dynamic environments in Southeast Asia. The forest area estimates from those maps vary substantially, ranging from 1.73×10(6) km(2) (GlobCover) to 2.69×10(6) km(2) (MCD12Q1) in 2009; and their uncertainty is constrained by frequent cloud cover and coarse spatial resolution. Recently, cloud-free imagery from the Phased Array Type L-band Synthetic Aperture Radar (PALSAR) onboard the Advanced Land Observing Satellite (ALOS) became available. We used the PALSAR 50-m orthorectified mosaic imagery in 2009 to generate a forest cover map of Southeast Asia at 50-m spatial resolution. The validation, using ground-reference data collected from the Geo-Referenced Field Photo Library and high-resolution images in Google Earth, showed that our forest map has a reasonably high accuracy (producer's accuracy 86% and user's accuracy 93%). The PALSAR-based forest area estimates in 2009 are significantly correlated with those from GlobCover and MCD12Q1 at national and subnational scales but differ in some regions at the pixel scale due to different spatial resolutions, forest definitions, and algorithms. The resultant 50-m forest map was used to quantify forest fragmentation and it revealed substantial details of forest fragmentation. This new 50-m map of tropical forests could serve as a baseline map for forest resource inventory, deforestation monitoring, reducing emissions from deforestation and forest degradation (REDD+) implementation, and biodiversity.

Rivers and flooded areas identified by medium-resolution remote sensing improve risk prediction of the highly pathogenic avian influenza H5N1 in Thailand.

Thailand experienced several epidemic waves of the highly pathogenic avian influenza (HPAI) H5N1 between 2004 and 2005. This study investigated the role of water in the landscape, which has not been previously assessed because of a lack of high-resolution information on the distribution of flooded land at the time of the epidemic. Nine Landsat 7 - Enhanced Thematic Mapper Plus scenes covering 174,610 km(2) were processed using k-means unsupervised classification to map the distribution of flooded areas as well as permanent lakes and reservoirs at the time of the main epidemic HPAI H5N1 wave of October 2004. These variables, together with other factors previously identified as significantly associated with risk, were entered into an autologistic regression model in order to quantify the gain in risk explanation over previously published models. We found that, in addition to other factors previously identified as associated with risk, the proportion of land covered by flooding along with expansion of rivers and streams, derived from an existing, sub-district level (administrative level no. 3) geographical information system database, was a highly significant risk factor in this 2004 HPAI epidemic. These results suggest that water-borne transmission could have partly contributed to the spread of HPAI H5N1 during the epidemic. Future work stemming from these results should involve studies where the actual distribution of small canals, rivers, ponds, rice paddy fields and farms are mapped and tested against farm-level data with respect to HPAI H5N1.

Improving risk models for avian influenza: the role of intensive poultry farming and flooded land during the 2004 Thailand epidemic.

Since 1996 when Highly Pathogenic Avian Influenza type H5N1 first emerged in southern China, numerous studies sought risk factors and produced risk maps based on environmental and anthropogenic predictors. However little attention has been paid to the link between the level of intensification of poultry production and the risk of outbreak. This study revised H5N1 risk mapping in Central and Western Thailand during the second wave of the 2004 epidemic. Production structure was quantified using a disaggregation methodology based on the number of poultry per holding. Population densities of extensively- and intensively-raised ducks and chickens were derived both at the sub-district and at the village levels. LandSat images were used to derive another previously neglected potential predictor of HPAI H5N1 risk: the proportion of water in the landscape resulting from floods. We used Monte Carlo simulation of Boosted Regression Trees models of predictor variables to characterize the risk of HPAI H5N1. Maps of mean risk and uncertainty were derived both at the sub-district and the village levels. The overall accuracy of Boosted Regression Trees models was comparable to that of logistic regression approaches. The proportion of area flooded made the highest contribution to predicting the risk of outbreak, followed by the densities of intensively-raised ducks, extensively-raised ducks and human population. Our results showed that as little as 15% of flooded land in villages is sufficient to reach the maximum level of risk associated with this variable. The spatial pattern of predicted risk is similar to previous work: areas at risk are mainly located along the flood plain of the Chao Phraya river and to the south-east of Bangkok. Using high-resolution village-level poultry census data, rather than sub-district data, the spatial accuracy of predictions was enhanced to highlight local variations in risk. Such maps provide useful information to guide intervention.

Flying over an infected landscape: distribution of highly pathogenic avian influenza H5N1 risk in South Asia and satellite tracking of wild waterfowl.

Highly pathogenic avian influenza (HPAI) H5N1 virus persists in Asia, posing a threat to poultry, wild birds, and humans. Previous work in Southeast Asia demonstrated that HPAI H5N1 risk is related to domestic ducks and people. Other studies discussed the role of migratory birds in the long distance spread of HPAI H5N1. However, the interplay between local persistence and long-distance dispersal has never been studied. We expand previous geospatial risk analysis to include South and Southeast Asia, and integrate the analysis with migration data of satellite-tracked wild waterfowl along the Central Asia flyway. We find that the population of domestic duck is the main factor delineating areas at risk of HPAI H5N1 spread in domestic poultry in South Asia, and that other risk factors, such as human population and chicken density, are associated with HPAI H5N1 risk within those areas. We also find that satellite tracked birds (Ruddy Shelduck and two Bar-headed Geese) reveal a direct spatio-temporal link between the HPAI H5N1 hot-spots identified in India and Bangladesh through our risk model, and the wild bird outbreaks in May-June-July 2009 in China (Qinghai Lake), Mongolia, and Russia. This suggests that the continental-scale dynamics of HPAI H5N1 are structured as a number of persistence areas delineated by domestic ducks, connected by rare transmission through migratory waterfowl.