State Key Laboratory of Remote Sensing Science

The ongoing coronavirus disease 2019 (COVID-19) outbreak expanded rapidly throughout China. Major behavioral, clinical, and state interventions were undertaken to mitigate the epidemic and prevent the persistence of the virus in human populations in China and worldwide. It remains unclear how these unprecedented interventions, including travel restrictions, affected COVID-19 spread in China. We used real-time mobility data from Wuhan and detailed case data including travel history to elucidate the role of case importation in transmission in cities across China and to ascertain the impact of control measures. Early on, the spatial distribution of COVID-19 cases in China was explained well by human mobility data. After the implementation of control measures, this correlation dropped and growth rates became negative in most locations, although shifts in the demographics of reported cases were still indicative of local chains of transmission outside of Wuhan. This study shows that the drastic control measures implemented in China substantially mitigated the spread of COVID-19.

Land surface temperature (LST) is one of the key parameters in the physics of land surface processes from local through global scales. The importance of LST is being increasingly recognized and there is a strong interest in developing methodologies to measure LST from space. However, retrieving LST is still a challenging task since the LST retrieval problem is ill-posed. This paper reviews the current status of selected remote sensing algorithms for estimating LST from thermal infrared (TIR) data. A brief theoretical background of the subject is presented along with a survey of the algorithms employed for obtaining LST from space-based TIR measurements. The discussion focuses on TIR data acquired from polar-orbiting satellites because of their widespread use, global applicability and higher spatial resolution compared to geostationary satellites. The theoretical framework and methodologies used to derive the LST from the data are reviewed followed by the methodologies for validating satellite-derived LST. Directions for future research to improve the accuracy of satellite-derived LST are then suggested.

Responding to an outbreak of a novel coronavirus [agent of coronavirus disease 2019 (COVID-19)] in December 2019, China banned travel to and from Wuhan city on 23 January 2020 and implemented a national emergency response. We investigated the spread and control of COVID-19 using a data set that included case reports, human movement, and public health interventions. The Wuhan shutdown was associated with the delayed arrival of COVID-19 in other cities by 2.91 days. Cities that implemented control measures preemptively reported fewer cases on average (13.0) in the first week of their outbreaks compared with cities that started control later (20.6). Suspending intracity public transport, closing entertainment venues, and banning public gatherings were associated with reductions in case incidence. The national emergency response appears to have delayed the growth and limited the size of the COVID-19 epidemic in China, averting hundreds of thousands of cases by 19 February (day 50).

We have produced the first 30 m resolution global land-cover maps using Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper Plus (ETM+) data. We have classified over 6600 scenes of Landsat TM data after 2006, and over 2300 scenes of Landsat TM and ETM+ data before 2006, all selected from the green season. These images cover most of the world's land surface except Antarctica and Greenland. Most of these images came from the United States Geological Survey in level L1T (orthorectified). Four classifiers that were freely available were employed, including the conventional maximum likelihood classifier (MLC), J4.8 decision tree classifier, Random Forest (RF) classifier and support vector machine (SVM) classifier. A total of 91,433 training samples were collected by traversing each scene and finding the most representative and homogeneous samples. A total of 38,664 test samples were collected at preset, fixed locations based on a globally systematic unaligned sampling strategy. Two software tools, Global Analyst and Global Mapper developed by extending the functionality of Google Earth, were used in developing the training and test sample databases by referencing the Moderate Resolution Imaging Spectroradiometer enhanced vegetation index (MODIS EVI) time series for 2010 and high resolution images from Google Earth. A unique land-cover classification system was developed that can be crosswalked to the existing United Nations Food and Agriculture Organization (FAO) land-cover classification system as well as the International Geosphere-Biosphere Programme (IGBP) system. Using the four classification algorithms, we obtained the initial set of global land-cover maps. The SVM produced the highest overall classification accuracy (OCA) of 64.9% assessed with our test samples, with RF (59.8%), J4.8 (57.9%), and MLC (53.9%) ranked from the second to the fourth. We also estimated the OCAs using a subset of our test samples (8629) each of which represented a homogeneous area greater than 500 m × 500 m. Using this subset, we found the OCA for the SVM to be 71.5%. As a consistent source for estimating the coverage of global land-cover types in the world, estimation from the test samples shows that only 6.90% of the world is planted for agricultural production. The total area of cropland is 11.51% if unplanted croplands are included. The forests, grasslands, and shrublands cover 28.35%, 13.37%, and 11.49% of the world, respectively. The impervious surface covers only 0.66% of the world. Inland waterbodies, barren lands, and snow and ice cover 3.56%, 16.51%, and 12.81% of the world, respectively.

Separating point clouds into ground and non-ground measurements is an essential step to generate digital terrain models (DTMs) from airborne LiDAR (light detection and ranging) data. However, most filtering algorithms need to carefully set up a number of complicated parameters to achieve high accuracy. In this paper, we present a new filtering method which only needs a few easy-to-set integer and Boolean parameters. Within the proposed approach, a LiDAR point cloud is inverted, and then a rigid cloth is used to cover the inverted surface. By analyzing the interactions between the cloth nodes and the corresponding LiDAR points, the locations of the cloth nodes can be determined to generate an approximation of the ground surface. Finally, the ground points can be extracted from the LiDAR point cloud by comparing the original LiDAR points and the generated surface. Benchmark datasets provided by ISPRS (International Society for Photogrammetry and Remote Sensing) working Group III/3 are used to validate the proposed filtering method, and the experimental results yield an average total error of 4.58%, which is comparable with most of the state-of-the-art filtering algorithms. The proposed easy-to-use filtering method may help the users without much experience to use LiDAR data and related technology in their own applications more easily.

The application of deep learning approaches to finance has received a great deal of attention from both investors and researchers. This study presents a novel deep learning framework where wavelet transforms (WT), stacked autoencoders (SAEs) and long-short term memory (LSTM) are combined for stock price forecasting. The SAEs for hierarchically extracted deep features is introduced into stock price forecasting for the first time. The deep learning framework comprises three stages. First, the stock price time series is decomposed by WT to eliminate noise. Second, SAEs is applied to generate deep high-level features for predicting the stock price. Third, high-level denoising features are fed into LSTM to forecast the next day's closing price. Six market indices and their corresponding index futures are chosen to examine the performance of the proposed model. Results show that the proposed model outperforms other similar models in both predictive accuracy and profitability performance.

Climate change refers to long-term shifts in weather conditions and patterns of extreme weather events. It may lead to changes in health threat to human beings, multiplying existing health problems. This review examines the scientific evidences on the impact of climate change on human infectious diseases. It identifies research progress and gaps on how human society may respond to, adapt to, and prepare for the related changes. Based on a survey of related publications between 1990 and 2015, the terms used for literature selection reflect three aspects--the components of infectious diseases, climate variables, and selected infectious diseases. Humans' vulnerability to the potential health impacts by climate change is evident in literature. As an active agent, human beings may control the related health effects that may be effectively controlled through adopting proactive measures, including better understanding of the climate change patterns and of the compound disease-specific health effects, and effective allocation of technologies and resources to promote healthy lifestyles and public awareness. The following adaptation measures are recommended: 1) to go beyond empirical observations of the association between climate change and infectious diseases and develop more scientific explanations, 2) to improve the prediction of spatial-temporal process of climate change and the associated shifts in infectious diseases at various spatial and temporal scales, and 3) to establish locally effective early warning systems for the health effects of predicated climate change.

Climate conditions significantly affect vegetation growth in terrestrial ecosystems. Due to the spatial heterogeneity of ecosystems, the vegetation responses to climate vary considerably with the diverse spatial patterns and the time-lag effects, which are the most important mechanism of climate-vegetation interactive effects. Extensive studies focused on large-scale vegetation-climate interactions use the simultaneous meteorological and vegetation indicators to develop models; however, the time-lag effects are less considered, which tends to increase uncertainty. In this study, we aim to quantitatively determine the time-lag effects of global vegetation responses to different climatic factors using the GIMMS3g NDVI time series and the CRU temperature, precipitation, and solar radiation datasets. First, this study analyzed the time-lag effects of global vegetation responses to different climatic factors. Then, a multiple linear regression model and partial correlation model were established to statistically analyze the roles of different climatic factors on vegetation responses, from which the primary climate-driving factors for different vegetation types were determined. The results showed that (i) both the time-lag effects of the vegetation responses and the major climate-driving factors that significantly affect vegetation growth varied significantly at the global scale, which was related to the diverse vegetation and climate characteristics; (ii) regarding the time-lag effects, the climatic factors explained 64% variation of the global vegetation growth, which was 11% relatively higher than the model ignoring the time-lag effects; (iii) for the area with a significant change trend (for the period 1982-2008) in the global GIMMS3g NDVI (P < 0.05), the primary driving factor was temperature; and (iv) at the regional scale, the variation in vegetation growth was also related to human activities and natural disturbances. Considering the time-lag effects is quite important for better predicting and evaluating the vegetation dynamics under the background of global climate change.

Abstract. The land and ocean absorb on average just over half of the anthropogenic emissions of carbon dioxide (CO2) every year. These CO2 "sinks" are modulated by climate change and variability. Here we use a suite of nine dynamic global vegetation models (DGVMs) and four ocean biogeochemical general circulation models (OBGCMs) to estimate trends driven by global and regional climate and atmospheric CO2 in land and oceanic CO2 exchanges with the atmosphere over the period 1990–2009, to attribute these trends to underlying processes in the models, and to quantify the uncertainty and level of inter-model agreement. The models were forced with reconstructed climate fields and observed global atmospheric CO2; land use and land cover changes are not included for the DGVMs. Over the period 1990–2009, the DGVMs simulate a mean global land carbon sink of −2.4 ± 0.7 Pg C yr−1 with a small significant trend of −0.06 ± 0.03 Pg C yr−2 (increasing sink). Over the more limited period 1990–2004, the ocean models simulate a mean ocean sink of −2.2 ± 0.2 Pg C yr−1 with a trend in the net C uptake that is indistinguishable from zero (−0.01 ± 0.02 Pg C yr−2). The two ocean models that extended the simulations until 2009 suggest a slightly stronger, but still small, trend of −0.02 ± 0.01 Pg C yr−2. Trends from land and ocean models compare favourably to the land greenness trends from remote sensing, atmospheric inversion results, and the residual land sink required to close the global carbon budget. Trends in the land sink are driven by increasing net primary production (NPP), whose statistically significant trend of 0.22 ± 0.08 Pg C yr−2 exceeds a significant trend in heterotrophic respiration of 0.16 ± 0.05 Pg C yr−2 – primarily as a consequence of widespread CO2 fertilisation of plant production. Most of the land-based trend in simulated net carbon uptake originates from natural ecosystems in the tropics (−0.04 ± 0.01 Pg C yr−2), with almost no trend over the northern land region, where recent warming and reduced rainfall offsets the positive impact of elevated atmospheric CO2 and changes in growing season length on carbon storage. The small uptake trend in the ocean models emerges because climate variability and change, and in particular increasing sea surface temperatures, tend to counter\\-act the trend in ocean uptake driven by the increase in atmospheric CO2. Large uncertainty remains in the magnitude and sign of modelled carbon trends in several regions, as well as regarding the influence of land use and land cover changes on regional trends.

A major research plan entitled “Integrated research on the ecohydrological process of the Heihe River Basin” was launched by the National Natural Science Foundation of China in 2010. One of the key aims of this research plan is to establish a research platform that integrates observation, data management, and model simulation to foster twenty-first-century watershed science in China. Based on the diverse needs of interdisciplinary studies within this research plan, a program called the Heihe Watershed Allied Telemetry Experimental Research (HiWATER) was implemented. The overall objective of HiWATER is to improve the observability of hydrological and ecological processes, to build a world-class watershed observing system, and to enhance the applicability of remote sensing in integrated ecohydrological studies and water resource management at the basin scale. This paper introduces the background, scientific objectives, and experimental design of HiWATER. The instrumental setting and airborne mission plans are also outlined. The highlights are the use of a flux observing matrix and an eco-hydrological wireless sensor network to capture multiscale heterogeneities and to address complex problems, such as heterogeneity, scaling, uncertainty, and closing water cycle at the watershed scale. HiWATER was formally initialized in May 2012 and will last four years until 2015. Data will be made available to the scientific community via the Environmental and Ecological Science Data Center for West China. International scientists are welcome to participate in the field campaign and use the data in their analyses.

It is more and more acknowledged that land-use/cover dynamic change has become a key subject urgently to be dealt with in the study of global environmental change. Supported by the Landsat TM digital images, spatial patterns and temporal variation of land-use change during 1995 –2000 are studied in the paper. According to the land-use dynamic degree model, supported by the 1km GRID data of land-use change and the comprehensive characters of physical, economic and social features, a dynamic regionalization of land-use change is designed to disclose the spatial pattern of land-use change processes. Generally speaking, in the traditional agricultural zones, e.g., Huang-Huai-Hai Plains, Yangtze River Delta and Sichuan Basin, the built-up and residential areas occupy a great proportion of arable land, and in the interlock area of farming and pasturing of northern China and the oases agricultural zones, the reclamation of arable land is conspicuously driven by changes of production conditions, economic benefits and climatic conditions. The implementation of “returning arable land into woodland or grassland” policies has won initial success in some areas, but it is too early to say that the trend of deforestation has been effectively reversed across China. In this paper, the division of dynamic regionalization of land-use change is designed, for the sake of revealing the temporal and spatial features of land-use change and laying the foundation for the study of regional scale land-use changes. Moreover, an integrated study, including studies of spatial pattern and temporal process of land-use change, is carried out in this paper, which is an interesting try on the comparative studies of spatial pattern on change process and the change process of spatial pattern of land-use change.

ABSTRACT Land cover changes ( LCCs ) play an important role in the climate system. Research over recent decades highlights the impacts of these changes on atmospheric temperature, humidity, cloud cover, circulation, and precipitation. These impacts range from the local‐ and regional‐scale to sub‐continental and global‐scale. It has been found that the impacts of regional‐scale LCC in one area may also be manifested in other parts of the world as a climatic teleconnection. In light of these findings, this article provides an overview and synthesis of some of the most notable types of LCC and their impacts on climate. These LCC types include agriculture, deforestation and afforestation, desertification, and urbanization. In addition, this article provides a discussion on challenges to, and future research directions in, assessing the climatic impacts of LCC .

Leaf area index (LAI) products at regional and global scales are being routinely generated from individual instrument data acquired at a specific time. As a result of cloud contamination and other factors, these LAI products are spatially and temporally discontinuous and are also inaccurate for some vegetation types in many areas. A better strategy is to use multi-temporal data. In this paper, a method was developed to estimate LAI from time-series remote sensing data using general regression neural networks (GRNNs). A database was generated from Moderate-Resolution Imaging Spectroradiometer (MODIS) and CYCLOPES LAI products as well as MODIS reflectance products of the BELMANIP sites during the period from 2001-2003. The effective CYCLOPES LAI was first converted to true LAI, which was then combined with the MODIS LAI according to their uncertainties determined from the ground-measured true LAI. The MODIS reflectance was reprocessed to remove remaining effects. GRNNs were then trained over the fused LAI and reprocessed MODIS reflectance for each biome type to retrieve LAI from time-series remote sensing data. The reprocessed MODIS reflectance data from an entire year were inputted into the GRNNs to estimate the 1-year LAI profiles. Extensive validations for all biome types were carried out, and it was demonstrated that the method is able to estimate temporally continuous LAI profiles with much improved accuracy compared with that of the current MODIS and CYCLOPES LAI products. This new method is being used to produce the Global Land Surface Satellite LAI products in China.

Abstract Urban boundaries, an essential property of cities, are widely used in many urban studies. However, extracting urban boundaries from satellite images is still a great challenge, especially at a global scale and a fine resolution. In this study, we developed an automatic delineation framework to generate a multi-temporal dataset of global urban boundaries (GUB) using 30 m global artificial impervious area (GAIA) data. First, we delineated an initial urban boundary by filling inner non-urban areas of each city. A kernel density estimation approach and cellular-automata based urban growth modeling were jointly used in this step. Second, we improved the initial urban boundaries around urban fringe areas, using a morphological approach by dilating and eroding the derived urban extent. We implemented this delineation on the Google Earth Engine platform and generated a 30 m resolution global urban boundary dataset in seven representative years (i.e. 1990, 1995, 2000, 2005, 2010, 2015, and 2018). Our extracted urban boundaries show a good agreement with results derived from nighttime light data and human interpretation, and they can well delineate the urban extent of cities when compared with high-resolution Google Earth images. The total area of 65 582 GUBs, each of which exceeds 1 km 2 , is 809 664 km 2 in 2018. The impervious surface areas account for approximately 60% of the total. From 1990 to 2018, the proportion of impervious areas in delineated boundaries increased from 53% to 60%, suggesting a compact urban growth over the past decades. We found that the United States has the highest per capita urban area (i.e. more than 900 m 2 ) among the top 10 most urbanized nations in 2018. This dataset provides a physical boundary of urban areas that can be used to study the impact of urbanization on food security, biodiversity, climate change, and urban health. The GUB dataset can be accessed from http://data.ess.tsinghua.edu.cn .

Air pollution over China has attracted wide interest from public and academic community. PM2.5 is the primary air pollutant across China. Quantifying interactions between meteorological conditions and PM2.5 concentrations are essential to understand the variability of PM2.5 and seek methods to control PM2.5. Since 2013, the measurement of PM2.5 has been widely made at 1436 stations across the country and more than 300 papers focusing on PM2.5-meteorology interactions have been published. This article is a comprehensive review on the meteorological impact on PM2.5 concentrations. We start with an introduction of general meteorological conditions and PM2.5 concentrations across China, and then seasonal and spatial variations of meteorological influences on PM2.5 concentrations. Next, major methods used to quantify meteorological influences on PM2.5 concentrations are checked and compared. We find that causality analysis methods are more suitable for extracting the influence of individual meteorological factors whilst statistical models are good at quantifying the overall effect of multiple meteorological factors on PM2.5 concentrations. Chemical Transport Models (CTMs) have the potential to provide dynamic estimation of PM2.5 concentrations by considering anthropogenic emissions and the transport and evolution of pollutants. We then comprehensively examine the mechanisms how major meteorological factors may impact the PM2.5 concentrations, including the dispersion, growth, chemical production, photolysis, and deposition of PM2.5. The feedback effects of PM2.5 concentrations on meteorological factors are also carefully examined. Based on this review, suggestions on future research and major meteorological approaches for mitigating PM2.5 pollution are made finally.

The tasselled cap transformation (TCT) is a useful tool for compressing spectral data into a few bands associated with physical scene characteristics with minimal information loss. TCT was originally evolved from the Landsat multi-spectral scanner (MSS) launched in 1972 and is widely adapted to modern sensors. In this study, we derived the TCT coefficients for the newly launched (2013) operational land imager (OLI) sensor on-board Landsat 8 for at-satellite reflectance. A newly developed standardized mechanism was used to transform the principal component analysis (PCA)-based rotated axes through Procrustes rotation (PR) conformation according to the Landsat thematic mapper (TM)-based tasselled cap space. Firstly, OLI data were transformed into TM TCT space directly and considered as a dummy target. Then, PCA was applied on the original scene. Finally, PR was applied to get the transformation results in the best conformation to the target image. New coefficients were analysed in detail to confirm Landsat 8-based TCT as a continuity of the original tasselled cap idea. Results show that newly derived set of coefficients for Landsat OLI is in continuation of its predecessors and hence provide data continuity through TCT since 1972 for remote sensing of surface features such as vegetation, albedo and water. The newly derived TCT for OLI will also be very useful for studying biomass estimation and primary production for future studies.

Abstract. Fine particulate matter with aerodynamic diameters ≤2.5 µm (PM2.5) has adverse effects on human health and the atmospheric environment. The estimation of surface PM2.5 concentrations has made intensive use of satellite-derived aerosol products. However, it has been a great challenge to obtain high-quality and high-resolution PM2.5 data from both ground and satellite observations, which is essential to monitor air pollution over small-scale areas such as metropolitan regions. Here, the space–time extremely randomized trees (STET) model was enhanced by integrating updated spatiotemporal information and additional auxiliary data to improve the spatial resolution and overall accuracy of PM2.5 estimates across China. To this end, the newly released Moderate Resolution Imaging Spectroradiometer Multi-Angle Implementation of Atmospheric Correction AOD product, along with meteorological, topographical and land-use data and pollution emissions, was input to the STET model, and daily 1 km PM2.5 maps for 2018 covering mainland China were produced. The STET model performed well, with a high out-of-sample (out-of-station) cross-validation coefficient of determination (R2) of 0.89 (0.88), a low root-mean-square error of 10.33 (10.93) µg m−3, a small mean absolute error of 6.69 (7.15) µg m−3 and a small mean relative error of 21.28 % (23.69 %). In particular, the model captured well the PM2.5 concentrations at both regional and individual site scales. The North China Plain, the Sichuan Basin and Xinjiang Province always featured high PM2.5 pollution levels, especially in winter. The STET model outperformed most models presented in previous related studies, with a strong predictive power (e.g., monthly R2=0.80), which can be used to estimate historical PM2.5 records. More importantly, this study provides a new approach for obtaining high-resolution and high-quality PM2.5 dataset across mainland China (i.e., ChinaHighPM2.5), important for air pollution studies focused on urban areas.

A near-real time drought monitoring approach is developed using Terra–Moderate Resolution Imaging Spectoradiometer (MODIS) Normalized Difference Vegetation Index (NDVI) and Land Surface Temperature (LST) products. The approach is called Vegetation Temperature Condition Index (VTCI), which integrates land surface reflectance and thermal properties. VTCI is defined as the ratio of LST differences among pixels with a specific NDVI value in a sufficiently large study area; the numerator is the difference between maximum LST of the pixels and LST of one pixel; and the denominator is the difference between maximum and minimum LSTs of the pixels. VTCI is lower for drought and higher for wet conditions. The ground-measured precipitation data from a study area covering parts of the states of Texas and Oklahoma in the southern Great Plains, USA, are used to validate the drought monitoring approach. Linear correlation analysis between VTCI, and total monthly precipitation and departure from normal monthly precipitation shows that VTCI is not only closely related to recent rainfall events but also related to past rainfall amounts, and indicates that VTCI might be a better and a near-real time drought monitoring approach.

As one of the key tools for regulating human-ecosystem relations, environmental conservation policies can promote ecological rehabilitation across a variety of spatiotemporal scales. However, quantifying the ecological effects of such policies at the regional level is difficult. A case study was conducted at the regional level in the ecologically vulnerable region of the Loess Plateau, China, through the use of several methods including the Universal Soil Loss Equation (USLE), hydrological modeling and multivariate analysis. An assessment of the changes over the period of 2000-2008 in four key ecosystem services was undertaken to determine the effects of the Chinese government's ecological rehabilitation initiatives implemented in 1999. These ecosystem services included water regulation, soil conservation, carbon sequestration and grain production. Significant conversions of farmland to woodland and grassland were found to have resulted in enhanced soil conservation and carbon sequestration, but decreased regional water yield under a warming and drying climate trend. The total grain production increased in spite of a significant decline in farmland acreage. These trends have been attributed to the strong socioeconomic incentives embedded in the ecological rehabilitation policy. Although some positive policy results have been achieved over the last decade, large uncertainty remains regarding long-term policy effects on the sustainability of ecological rehabilitation performance and ecosystem service enhancement. To reduce such uncertainty, this study calls for an adaptive management approach to regional ecological rehabilitation policy to be adopted, with a focus on the dynamic interactions between people and their environments in a changing world.

Recently, five Global LAnd Surface Satellite (GLASS) products have been released: leaf area index (LAI), shortwave broadband albedo, longwave broadband emissivity, incident short radiation, and photosynthetically active radiation (PAR). The first three products cover the years 1982–2012 (LAI) and 1981–2010 (albedo and emissivity) at 1–5 km and 8-day resolutions, and the last two radiation products span the period 2008–2010 at 5 km and 3-h resolutions. These products have been evaluated and validated, and the preliminary results indicate that they are of higher quality and accuracy than the existing products. In particular, the first three products have much longer time series, and are therefore highly suitable for various environmental studies. This paper outlines the algorithms, product characteristics, preliminary validation results, potential applications and some examples of initial analysis of these products.