Russell Doughty

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.

Temporal Consistency Between Gross Primary Production and Solar-Induced Chlorophyll Fluorescence in the Ten Most Populous Megacity Areas over Years

The gross primary production (GPP) of vegetation in urban areas plays an important role in the study of urban ecology. It is difficult however, to accurately estimate GPP in urban areas, mostly due to the complexity of impervious land surfaces, buildings, vegetation, and management. Recently, we used the Vegetation Photosynthesis Model (VPM), climate data, and satellite images to estimate the GPP of terrestrial ecosystems including urban areas. Here, we report VPM-based GPP (GPPvpm) estimates for the world’s ten most populous megacities during 2000–2014. The seasonal dynamics of GPPvpm during 2007–2014 in the ten megacities track well that of the solar-induced chlorophyll fluorescence (SIF) data from GOME-2 at 0.5° × 0.5° resolution. Annual GPPvpm during 2000–2014 also shows substantial variation among the ten megacities, and year-to-year trends show increases, no change, and decreases. Urban expansion and vegetation collectively impact GPP variations in these megacities. The results of this study demonstrate the potential of a satellite-based vegetation photosynthesis model for diagnostic studies of GPP and the terrestrial carbon cycle in urban areas.

Divergent trends of open-surface water body area in the contiguous United States from 1984 to 2016

Significance Strong variations in open-surface water body areas have impacted United States agriculture, economy, society, and ecosystems. This study presents the uneven water-resource distribution across the contiguous United States with the western half of the United States having less water body area but stronger interannual variability compared with the eastern half. Divergent trends of open-surface water body area in the last three decades, mainly driven by climate, indicated that the water-poor regions of the Southwest and Northwest United States were getting poorer, while the water-rich regions of Southeast and far north Great Plains were getting richer. Surface water body shrinkage in drought years led to massive groundwater mining and the rapid decrease of land water storage in California and the southern Great Plains. The contiguous United States (CONUS), especially the West, faces challenges of increasing water stress and uncertain impacts of climate change. The historical information of surface water body distribution, variation, and multidecadal trends documented in remote-sensing images can aid in water-resource planning and management, yet is not well explored. Here, we detected open-surface water bodies in all Landsat 5, 7, and 8 images (∼370,000 images, >200 TB) of the CONUS and generated 30-meter annual water body frequency maps for 1984–2016. We analyzed the interannual variations and trends of year-long water body area, examined the impacts of climatic and anthropogenic drivers on water body area dynamics, and explored the relationships between water body area and land water storage (LWS). Generally, the western half of the United States is prone to water stress, with small water body area and large interannual variability. During 1984–2016, water-poor regions of the Southwest and Northwest had decreasing trends in water body area, while water-rich regions of the Southeast and far north Great Plains had increasing trends. These divergent trends, mainly driven by climate, enlarged water-resource gaps and are likely to continue according to climate projections. Water body area change is a good indicator of LWS dynamics in 58% of the CONUS. Following the 2012 prolonged drought, LWS in California and the southern Great Plains had a larger decrease than surface water body area, likely caused by massive groundwater withdrawals. Our findings provide valuable information for surface water-resource planning and management across the CONUS.

Continued decrease of open surface water body area in Oklahoma during 1984–2015

Oklahoma contains the largest number of manmade lakes and reservoirs in the United States. Despite the importance of these open surface water bodies to public water supply, agriculture, thermoelectric power, tourism and recreation, it is unclear how these water bodies have responded to climate change and anthropogenic water exploitation in past decades. In this study, we used all available Landsat 5 and 7 images (16,000 scenes) from 1984 through 2015 and a water index- and pixel-based approach to analyze the spatial-temporal variability of open surface water bodies and its relationship with climate and water exploitation. Specifically, the areas and numbers of four water body extents (the maximum, year-long, seasonal, and average extents) were analyzed to capture variations in water body area and number. Statistically significant downward trends were found in the maximum, year-long, and annual average water body areas from 1984 through 2015. Furthermore, these decreases were mainly attributed to the continued shrinking of large water bodies (>1km2). There were also significant decreases in maximum and year-long water body numbers, which suggested that some of the water bodies were vanishing year by year. However, remarkable inter-annual variations of water body area and number were also found. Both water body area and number were positively related to precipitation, and negatively related to temperature. Surface water withdrawals mainly influenced the year-long water bodies. The smaller water bodies have a higher risk of drying under a drier climate, which suggests that small water bodies are more vulnerable under climate-warming senarios.

Spatial-temporal consistency between gross primary productivity and solar-induced chlorophyll fluorescence of vegetation in China during 2007–2014

Accurately estimating spatial-temporal patterns of gross primary production (GPP) is important for the global carbon cycle. Satellite-based light use efficiency (LUE) models are regarded as an efficient tool in simulating spatial-temporal dynamics of GPP. However, the accuracy assessment of GPP simulations from LUE models at both spatial and temporal scales remains a challenge. In this study, we simulated GPP of vegetation in China during 2007-2014 using a LUE model (Vegetation Photosynthesis Model, VPM) based on MODIS (moderate-resolution imaging spectroradiometer) images with 8-day temporal and 500-m spatial resolutions and NCEP (National Center for Environmental Prediction) climate data. Global Ozone Monitoring Instrument 2 (GOME-2) solar-induced chlorophyll fluorescence (SIF) data were used to compare with VPM simulated GPP (GPPVPM) temporally and spatially using linear correlation analysis. Significant positive linear correlations exist between monthly GPPVPM and SIF data over a single year (2010) and multiple years (2007-2014) in most areas of China. GPPVPM is also significantly positive correlated with GOME-2 SIF (R2 > 0.43) spatially for seasonal scales. However, poor consistency was detected between GPPVPM and SIF data at yearly scale. GPP dynamic trends have high spatial-temporal variation in China during 2007-2014. Temperature, leaf area index (LAI), and precipitation are the most important factors influence GPPVPM in the regions of East Qinghai-Tibet Plateau, Loss Plateau, and Southwestern China, respectively. The results of this study indicate that GPPVPM is temporally and spatially in line with GOME-2 SIF data, and space-borne SIF data have great potential for evaluating LUE-based GPP models.

Identifying Establishment Year and Pre-Conversion Land Cover of Rubber Plantations on Hainan Island, China Using Landsat Data during 1987–2015

Knowing the stand age of rubber tree (Hevea brasiliensis) plantations is vitally important for best management practices, estimations of rubber latex yields, and carbon cycle studies (e.g., biomass, carbon pools, and fluxes). However, the stand age (as estimated from the establishment year of rubber plantation) is not available across large regions. In this study, we analyzed Landsat time series images from 1987–2015 and developed algorithms to identify (1) the establishment year of rubber plantations; and (2) the pre-conversion land cover types, such as old rubber plantations, evergreen forests, and cropland. Exposed soil during plantation establishment and linear increases in canopy closure during non-production periods (rubber seedling to mature plantation) were used to identify the establishment year of rubber plantations. Based on the rubber plantation map for 2015 (overall accuracy = 97%), and 1981 Landsat images since 1987, we mapped the establishment year of rubber plantations on Hainan Island (R2 = 0.85/0.99, and RMSE = 2.34/0.54 years at pixel/plantation scale). The results show that: (1) significant conversion of croplands and old rubber plantations to new rubber plantations has occurred substantially in the northwest and northern regions of Hainan Island since 2000, while old rubber plantations were mainly distributed in the southeastern inland strip; (2) the pattern of rubber plantation expansion since 1987 consisted of fragmented plantations from smallholders, and there was no tendency to expand towards a higher altitude and steep slope regions; (3) the largest land source for new rubber plantations since 1988 was old rubber plantations (1.26 × 105 ha), followed by cropland (0.95 × 105 ha), and evergreen forests (0.68 × 105 ha). The resultant algorithms and maps of establishment year and pre-conversion land cover types are likely to be useful in plantation management, and ecological assessments of rubber plantation expansion in China. Remote Sens. 2018, 10, 1240; doi:10.3390/rs10081240 www.mdpi.com/journal/remotesensing Remote Sens. 2018, 10, 1240 2 of 23

Quantifying spatial-temporal changes of tea plantations in complex landscapes through integrative analyses of optical and microwave imagery

Abstract High demand for tea has driven the expansion of tea plantations in the tropical and subtropical regions over the past few decades. Tea plant cultivation promotes economic development and creates job opportunities, but tea plantation expansion has significant impacts on biodiversity, carbon and water cycles, and ecosystem services. Mapping the spatial distribution and extent of tea plantations in a timely fashion is crucial for land use management and policy making. In this study, we mapped tea plantation expansion in Menghai County, Yunnan Province, China. We analyzed the structure and features of major land cover types in this tropical and subtropical region using (1) the HH and HV gamma-naught imagery from the Advanced Land Observation Satellite (ALOS) Phased Array L-band Synthetic Aperture Radar (PALSAR) and (2) time series Landsat TM/ETM+/OLI imagery. Tea plantation maps for 2010 and 2015 were generated using the pixel-based support vector machine (SVM) approach at 30 m resolution, which had high user/producer accuracies of 83.58%/91.67% and 87.50%/90.83%, respectively. The resultant maps show that tea plantation area increased by 33.56% (∼9335 ha), from ∼27,817 ha in 2010 to ∼37,152 ha in 2015. The additional tea plantation area was mainly converted from forest (32.50%) and cropland (67.50%). The results showed that the combination of PALSAR and optical data performed better in tea plantation mapping than using optical data only. This study provides a promising new approach to identify and map tea plantations in complex tropical landscapes at high spatial resolution.

Enhanced gross primary production and evapotranspiration in juniper encroached grasslands

Woody plant encroachment (WPE) into grasslands has been occurring globally and may be accelerated by climate change in the future. This land cover change is expected to alter the carbon and water cycles, but it remains uncertain how and to what extent the carbon and water cycles may change with WPE into grasslands under current climate. In this study, we examined the difference of vegetation indices (VIs), evapotranspiration (ET), gross primary production (GPP), and solar-induced chlorophyll fluorescence (SIF) during 2000-2010 between grasslands and juniper-encroached grasslands. We also quantitatively assessed the changes of GPP and ET for grasslands with different proportions of juniper encroachment (JWPE). Our results suggested that JWPE increased the GPP, ET, greenness-related VIs, and SIF of grasslands. Mean annual GPP and ET were, respectively, ~55% and ~45% higher when grasslands were completely converted into juniper forests under contemporary climate during 2000-2010. The enhancement of annual GPP and ET for grasslands with JWPE varied over years ranging from about +20% GPP (~+30% for ET) in the wettest year (2007) to about twice as much GPP (~+55% for ET) in the severe drought year (2006) relative to grasslands without encroachment. Additionally, the differences in GPP and ET showed significant seasonal dynamics. During the peak growing season (May-August), GPP and ET for grasslands with JWPE were ~30% and ~40% higher on average. This analysis provided insights into how and to what degree carbon and water cycles were impacted by JWPE, which is vital to understanding how JWPE and ecological succession will affect the regional and global carbon and water budgets in the future.

Expansion dynamics of deciduous rubber plantations in Xishuangbanna, China during 2000–2010

Monoculture rubber plantations have been replacing tropical rain forests substantially in Southern China and Southeast Asia over the past several decades, which have affected human wellbeing and ecosystem services. However, to the best of our knowledge on the extent of rubber plantation expansion and their stand ages is limited. We tracked the spatiotemporal dynamics of deciduous rubber plantations in Xishuangbanna, the second largest natural rubber production region in China, from 2000 to 2010 using time-series data from the Phased Array type L-band Synthetic Aperture Radar (PALSAR), Landsat, and Moderate Resolution Imaging Spectroradiometer (MODIS). We found that rubber plantations have been expanding across a gradient from the low-elevation plains to the high elevation mountains. The areas of deciduous rubber plantations with stand ages ≤5, 6–10, and ≥11-year old were ~1.2 × 105 ha, ~0.8 × 105 ha, and ~2.9 × 105 ha, respectively. Older rubber plantations were mainly located in low-elevation and species-rich regions (500–900 m) and younger rubber trees were distributed in areas of relative high-elevation with fragile ecosystems. Economic and market factors have driven the expansion of rubber plantations, which is not only a threat to biodiversity and environmental sustainability, but also a trigger for climatic disasters. This study illustrates that the integration of microwave, optical, and thermal data is an effective method for mapping deciduous rubber plantations in tropical mountainous regions and determining their stand ages. Our results demonstrate the spatiotemporal pattern of rubber expansions over the first decade of this century.