Yuanwei Qin

Spatiotemporal characteristics, patterns, and causes of land-use changes in China since the late 1980s

Land-use/land-cover changes (LUCCs) have links to both human and nature interactions. China’s Land-Use/cover Datasets (CLUDs) were updated regularly at 5-year intervals from the late 1980s to 2010, with standard procedures based on Landsat TMETM+ images. A land-use dynamic regionalization method was proposed to analyze major land-use conversions. The spatiotemporal characteristics, differences, and causes of land-use changes at a national scale were then examined. The main findings are summarized as follows.Land-use changes (LUCs) across China indicated a significant variation in spatial and temporal characteristics in the last 20 years (1990–2010). The area of cropland change decreased in the south and increased in the north, but the total area remained almost unchanged. The reclaimed cropland was shifted from the northeast to the northwest. The built-up lands expanded rapidly, were mainly distributed in the east, and gradually spread out to central and western China. Woodland decreased first, and then increased, but desert area was the opposite. Grassland continued decreasing. Different spatial patterns of LUC in China were found between the late 20th century and the early 21st century. The original 13 LUC zones were replaced by 15 units with changes of boundaries in some zones. The main spatial characteristics of these changes included (1) an accelerated expansion of built-up land in the Huang-Huai-Hai region, the southeastern coastal areas, the midstream area of the Yangtze River, and the Sichuan Basin; (2) shifted land reclamation in the north from northeast China and eastern Inner Mongolia to the oasis agricultural areas in northwest China; (3) continuous transformation from rain-fed farmlands in northeast China to paddy fields; and (4) effectiveness of the “Grain for Green” project in the southern agricultural-pastoral ecotones of Inner Mongolia, the Loess Plateau, and southwestern mountainous areas. In the last two decades, although climate change in the north affected the change in cropland, policy regulation and economic driving forces were still the primary causes of LUC across China. During the first decade of the 21st century, the anthropogenic factors that drove variations in land-use patterns have shifted the emphasis from one-way land development to both development and conservation.The “dynamic regionalization method” was used to analyze changes in the spatial patterns of zoning boundaries, the internal characteristics of zones, and the growth and decrease of units. The results revealed “the pattern of the change process,” namely the process of LUC and regional differences in characteristics at different stages. The growth and decrease of zones during this dynamic LUC zoning, variations in unit boundaries, and the characteristics of change intensities between the former and latter decades were examined. The patterns of alternative transformation between the “pattern” and “process” of land use and the causes for changes in different types and different regions of land use were explored.

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.

Mapping paddy rice planting area in rice-wetland coexistent areas through analysis of Landsat 8 OLI and MODIS images.

Accurate and up-to-date information on the spatial distribution of paddy rice fields is necessary for the studies of trace gas emissions, water source management, and food security. The phenology-based paddy rice mapping algorithm, which identifies the unique flooding stage of paddy rice, has been widely used. However, identification and mapping of paddy rice in rice-wetland coexistent areas is still a challenging task. In this study, we found that the flooding/transplanting periods of paddy rice and natural wetlands were different. The natural wetlands flood earlier and have a shorter duration than paddy rice in the Panjin Plain, a temperate region in China. We used this asynchronous flooding stage to extract the paddy rice planting area from the rice-wetland coexistent area. MODIS Land Surface Temperature (LST) data was used to derive the temperature-defined plant growing season. Landsat 8 OLI imagery was used to detect the flooding signal and then paddy rice was extracted using the difference in flooding stages between paddy rice and natural wetlands. The resultant paddy rice map was evaluated with in-situ ground-truth data and Google Earth images. The estimated overall accuracy and Kappa coefficient were 95% and 0.90, respectively. The spatial pattern of OLI-derived paddy rice map agrees well with the paddy rice layer from the National Land Cover Dataset from 2010 (NLCD-2010). The differences between RiceLandsat and RiceNLCD are in the range of ±20% for most 1-km grid cell. The results of this study demonstrate the potential of the phenology-based paddy rice mapping algorithm, via integrating MODIS and Landsat 8 OLI images, to map paddy rice fields in complex landscapes of paddy rice and natural wetland in the temperate region.

Mapping Deciduous Rubber Plantation Areas and Stand Ages with PALSAR and Landsat Images

Accurate and updated finer resolution maps of rubber plantations and stand ages are needed to understand and assess the impacts of rubber plantations on regional ecosystem processes. This study presented a simple method for mapping rubber plantation areas and their stand ages by integration of PALSAR 50-m mosaic images and multi-temporal Landsat TM/ETM+ images. The L-band PALSAR 50-m mosaic images were used to map forests (including both natural forests and rubber trees) and non-forests. For those PALSAR-based forest pixels, we analyzed the multi-temporal Landsat TM/ETM+ images from 2000 to 2009. We first studied phenological signatures of deciduous rubber plantations (defoliation and foliation) and natural forests through analysis of surface reflectance, Normal Difference Vegetation Index (NDVI), Enhanced Vegetation Index (EVI), and Land Surface Water Index (LSWI) and generated a map of rubber plantations in 2009. We then analyzed phenological signatures of rubber plantations with different stand ages and generated a map, in 2009, of rubber plantation stand ages ( 10 years-old) based on multi-temporal Landsat images. The resultant maps clearly illustrated how rubber plantations have expanded into the mountains in the study area over the years. The results in this study demonstrate the potential of integrating microwave (e.g., PALSAR) and optical remote sensing in the characterization of rubber plantations and their expansion over time.

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 × 106 km2. 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.

Mapping paddy rice planting area in wheat-rice double-cropped areas through integration of Landsat-8 OLI, MODIS, and PALSAR images.

As farmland systems vary over space and time (season and year), accurate and updated maps of paddy rice are needed for studies of food security and environmental problems. We selected a wheat-rice double-cropped area from fragmented landscapes along the rural–urban complex (Jiangsu Province, China) and explored the potential utility of integrating time series optical images (Landsat-8, MODIS) and radar images (PALSAR) in mapping paddy rice planting areas. We first identified several main types of non-cropland land cover and then identified paddy rice fields by selecting pixels that were inundated only during paddy rice flooding periods. These key temporal windows were determined based on MODIS Land Surface Temperature and vegetation indices. The resultant paddy rice map was evaluated using regions of interest (ROIs) drawn from multiple high-resolution images, Google Earth, and in-situ cropland photos. The estimated overall accuracy and Kappa coefficient were 89.8% and 0.79, respectively. In comparison with the National Land Cover Data (China) from 2010, the resultant map better detected changes in the paddy rice fields and revealed more details about their distribution. These results demonstrate the efficacy of using images from multiple sources to generate paddy rice maps for two-crop rotation systems.

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 distribution using multi-temporal Landsat imagery in the Sanjiang Plain, northeast China.

Information of paddy rice distribution is essential for food production and methane emission calculation. Phenology-based algorithms have been utilized in the mapping of paddy rice fields by identifying the unique flooding and seedling transplanting phases using multi-temporal moderate resolution (500 m to 1 km) images. In this study, we developed simple algorithms to identify paddy rice at a fine resolution at the regional scale using multi-temporal Landsat imagery. Sixteen Landsat images from 2010–2012 were used to generate the 30 m paddy rice map in the Sanjiang Plain, northeast China—one of the major paddy rice cultivation regions in China. Three vegetation indices, Normalized Difference Vegetation Index (NDVI), Enhanced Vegetation Index (EVI), and Land Surface Water Index (LSWI), were used to identify rice fields during the flooding/transplanting and ripening phases. The user and producer accuracies of paddy rice on the resultant Landsat-based paddy rice map were 90% and 94%, respectively. The Landsat-based paddy rice map was an improvement over the paddy rice layer on the National Land Cover Dataset, which was generated through visual interpretation and digitalization on the fine-resolution images. The agricultural census data substantially underreported paddy rice area, raising serious concern about its use for studies on food security.

Impacts of ecological restoration projects on agricultural productivity in China

The changes in cropland quantity and quality due to land use are critical concerns to national food security, particularly for China. Despite the significant ecological effects, the ecological restoration program (ERP), started from 1999, has evidently altered the spatial patterns of China’s cropland and agricultural productivity. Based on cropland dynamic data from 2000 to 2008 primarily derived from satellite images with a 30-m resolution and satellite-based net primary productivity models, we identified the impacts on agricultural productivity caused by ERP, including “Grain for Green” Program (GFGP) and “Reclaimed Cropland to Lake” (RCTL) Program. Our results indicated that the agricultural productivity lost with a rate of 132.67×104 t/a due to ERP, which accounted for 44.01% of the total loss rate caused by land use changes during 2000–2005. During 2005–2008, the loss rate due to ERP decreased to 77.18×104 t/a, which was equivalent to 58.17% of that in the first five years and 30.22% of the total loss rate caused by land use changes. The agricultural productivity loss from 2000–2008 caused by ERP was more attributed to GFGP (about 70%) than RCTL. Although ERP had a certain influence on cropland productivity during 2000–2008, its effect was still much less than that of urbanization; moreover, ERP was already converted from the project implementation phase to the consolidation phase.

A global moderate resolution dataset of gross primary production of vegetation for 2000–2016

Accurate estimation of the gross primary production (GPP) of terrestrial vegetation is vital for understanding the global carbon cycle and predicting future climate change. Multiple GPP products are currently available based on different methods, but their performances vary substantially when validated against GPP estimates from eddy covariance data. This paper provides a new GPP dataset at moderate spatial (500 m) and temporal (8-day) resolutions over the entire globe for 2000–2016. This GPP dataset is based on an improved light use efficiency theory and is driven by satellite data from MODIS and climate data from NCEP Reanalysis II. It also employs a state-of-the-art vegetation index (VI) gap-filling and smoothing algorithm and a separate treatment for C3/C4 photosynthesis pathways. All these improvements aim to solve several critical problems existing in current GPP products. With a satisfactory performance when validated against in situ GPP estimates, this dataset offers an alternative GPP estimate for regional to global carbon cycle studies.