Institute of Remote Sensing and Digital Earth

In response to ecosystem degradation from rapid economic development, China began investing heavily in protecting and restoring natural capital starting in 2000. We report on China's first national ecosystem assessment (2000-2010), designed to quantify and help manage change in ecosystem services, including food production, carbon sequestration, soil retention, sandstorm prevention, water retention, flood mitigation, and provision of habitat for biodiversity. Overall, ecosystem services improved from 2000 to 2010, apart from habitat provision. China's national conservation policies contributed significantly to the increases in those ecosystem services.

Aridity, which is increasing worldwide because of climate change, affects the structure and functioning of dryland ecosystems. Whether aridification leads to gradual (versus abrupt) and systemic (versus specific) ecosystem changes is largely unknown. We investigated how 20 structural and functional ecosystem attributes respond to aridity in global drylands. Aridification led to systemic and abrupt changes in multiple ecosystem attributes. These changes occurred sequentially in three phases characterized by abrupt decays in plant productivity, soil fertility, and plant cover and richness at aridity values of 0.54, 0.7, and 0.8, respectively. More than 20% of the terrestrial surface will cross one or several of these thresholds by 2100, which calls for immediate actions to minimize the negative impacts of aridification on essential ecosystem services for the more than 2 billion people living in drylands.

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.

China's terrestrial ecosystems have functioned as important carbon sinks. However, previous estimates of carbon budgets have included large uncertainties owing to the limitations of sample size, multiple data sources, and inconsistent methodologies. In this study, we conducted an intensive field campaign involving 14,371 field plots to investigate all sectors of carbon stocks in China's forests, shrublands, grasslands, and croplands to better estimate the regional and national carbon pools and to explore the biogeographical patterns and potential drivers of these pools. The total carbon pool in these four ecosystems was 79.24 ± 2.42 Pg C, of which 82.9% was stored in soil (to a depth of 1 m), 16.5% in biomass, and 0.60% in litter. Forests, shrublands, grasslands, and croplands contained 30.83 ± 1.57 Pg C, 6.69 ± 0.32 Pg C, 25.40 ± 1.49 Pg C, and 16.32 ± 0.41 Pg C, respectively. When all terrestrial ecosystems are taken into account, the country's total carbon pool is 89.27 ± 1.05 Pg C. The carbon density of the forests, shrublands, and grasslands exhibited a strong correlation with climate: it decreased with increasing temperature but increased with increasing precipitation. Our analysis also suggests a significant sequestration potential of 1.9-3.4 Pg C in forest biomass in the next 10-20 years assuming no removals, mainly because of forest growth. Our results update the estimates of carbon pools in China's terrestrial ecosystems based on direct field measurements, and these estimates are essential to the validation and parameterization of carbon models in China and globally.

Significance China has launched six key ecological restoration projects since the late 1970s, but the contribution of these projects to terrestrial C sequestration remains unknown. In this study we examined the ecosystem C sink in the project area (∼16% of the country’s land area) and evaluated the project-induced C sequestration. The total annual C sink in the project area between 2001 and 2010 was estimated to be 132 Tg C per y, over half of which (74 Tg C per y, 56%) was caused by the implementation of the six projects. This finding indicates that the implementation of the ecological restoration projects in China has significantly increased ecosystem C sequestration across the country.

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 .

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.

As a list of remotely sensed data sources is available, how to efficiently exploit useful information from multisource data for better Earth observation becomes an interesting but challenging problem. In this paper, the classification fusion of hyperspectral imagery (HSI) and data from other multiple sensors, such as light detection and ranging (LiDAR) data, is investigated with the state-of-the-art deep learning, named the two-branch convolution neural network (CNN). More specific, a two-tunnel CNN framework is first developed to extract spectral-spatial features from HSI; besides, the CNN with cascade block is designed for feature extraction from LiDAR or high-resolution visual image. In the feature fusion stage, the spatial and spectral features of HSI are first integrated in a dual-tunnel branch, and then combined with other data features extracted from a cascade network. Experimental results based on several multisource data demonstrate the proposed two-branch CNN that can achieve more excellent classification performance than some existing methods.

The inkjet technique has the capability of generating droplets in the picoliter volume range, firing thousands of times in a few seconds and printing in the noncontact manner. Since its emergence, inkjet technology has been widely utilized in the publishing industry for printing of text and pictures. As the technology developed, its applications have been expanded from two-dimensional (2D) to three-dimensional (3D) and even used to fabricate components of electronic devices. At the end of the twentieth century, researchers were aware of the potential value of this technology in life sciences and tissue engineering because its picoliter-level printing unit is suitable for depositing biological components. Currently inkjet technology has been becoming a practical tool in modern medicine serving for drug development, scaffold building, and cell depositing. In this article, we first review the history, principles and different methods of developing this technology. Next, we focus on the recent achievements of inkjet printing in the biological field. Inkjet bioprinting of generic biomaterials, biomacromolecules, DNAs, and cells and their major applications are introduced in order of increasing complexity. The current limitations/challenges and corresponding solutions of this technology are also discussed. A new concept, biopixels, is put forward with a combination of the key characteristics of inkjet printing and basic biological units to bring a comprehensive view on inkjet-based bioprinting. Finally, a roadmap of the entire 3D bioprinting is depicted at the end of this review article, clearly demonstrating the past, present, and future of 3D bioprinting and our current progress in this field.

Unmanned aerial vehicle (UAV) remote sensing has great potential for vegetation mapping in complex urban landscapes due to the ultra-high resolution imagery acquired at low altitudes. Because of payload capacity restrictions, off-the-shelf digital cameras are widely used on medium and small sized UAVs. The limitation of low spectral resolution in digital cameras for vegetation mapping can be reduced by incorporating texture features and robust classifiers. Random Forest has been widely used in satellite remote sensing applications, but its usage in UAV image classification has not been well documented. The objectives of this paper were to propose a hybrid method using Random Forest and texture analysis to accurately differentiate land covers of urban vegetated areas, and analyze how classification accuracy changes with texture window size. Six least correlated second-order texture measures were calculated at nine different window sizes and added to original Red-Green-Blue (RGB) images as ancillary data. A Random Forest classifier consisting of 200 decision trees was used for classification in the spectral-textural feature space. Results indicated the following: (1) Random Forest outperformed traditional Maximum Likelihood classifier and showed similar performance to object-based image analysis in urban vegetation classification; (2) the inclusion of texture features improved classification accuracy significantly; (3) classification accuracy followed an inverted U relationship with texture window size. The results demonstrate that UAV provides an efficient and ideal platform for urban vegetation mapping. The hybrid method proposed in this paper shows good performance in differentiating urban vegetation mapping. The drawbacks of off-the-shelf digital cameras can be reduced by adopting Random Forest and texture analysis at the same time.

) have resulted from land management failures, including forest fires and insect pests. Overall land-use change and land management have contributed about 1.45 Pg of carbon to the total carbon released from 1990 to 2010. Our results highlight the importance of improving land-use management, especially in view of the recently proposed expansion of urban areas in China.

This paper reviewed major remote sensing image classification techniques, including pixel-wise, sub-pixel-wise, and object-based image classification methods, and highlighted the importance of incorporating spatio-contextual information in remote sensing image classification. Further, this paper grouped spatio-contextual analysis techniques into three major categories, including 1) texture extraction, 2) Markov random fields (MRFs) modeling, and 3) image segmentation and object-based image analysis. Finally, this paper argued the necessity of developing geographic information analysis models for spatial-contextual classifications using two case studies.

In recent years, air pollution has become an important public health concern. The high concentration of fine particulate matter with diameter less than 2.5 µm (PM2.5) is known to be associated with lung cancer, cardiovascular disease, respiratory disease, and metabolic disease. Predicting PM2.5 concentrations can help governments warn people at high risk, thus mitigating the complications. Although attempts have been made to predict PM2.5 concentrations, the factors influencing PM2.5 prediction have not been investigated. In this work, we study feature importance for PM2.5 prediction in Tehran’s urban area, implementing random forest, extreme gradient boosting, and deep learning machine learning (ML) approaches. We use 23 features, including satellite and meteorological data, ground-measured PM2.5, and geographical data, in the modeling. The best model performance obtained was R2 = 0.81 (R = 0.9), MAE = 9.93 µg/m3, and RMSE = 13.58 µg/m3 using the XGBoost approach, incorporating elimination of unimportant features. However, all three ML methods performed similarly and R2 varied from 0.63 to 0.67, when Aerosol Optical Depth (AOD) at 3 km resolution was included, and 0.77 to 0.81, when AOD at 3 km resolution was excluded. Contrary to the PM2.5 lag data, satellite-derived AODs did not improve model performance.

Abstract Land use policies have turned southern China into one of the most intensively managed forest regions in the world, with actions maximizing forest cover on soils with marginal agricultural potential while concurrently increasing livelihoods and mitigating climate change. Based on satellite observations, here we show that diverse land use changes in southern China have increased standing aboveground carbon stocks by 0.11 ± 0.05 Pg C y −1 during 2002–2017. Most of this regional carbon sink was contributed by newly established forests (32%), while forests already existing contributed 24%. Forest growth in harvested forest areas contributed 16% and non-forest areas contributed 28% to the carbon sink, while timber harvest was tripled. Soil moisture declined significantly in 8% of the area. We demonstrate that land management in southern China has been removing an amount of carbon equivalent to 33% of regional fossil CO 2 emissions during the last 6 years, but forest growth saturation, land competition for food production and soil-water depletion challenge the longevity of this carbon sink service.

With the launch of space-borne satellites, more synthetic aperture radar (SAR) images are available than ever before, thus making dynamic ship monitoring possible. Object detectors in deep learning achieve top performance, benefitting from a free public dataset. Unfortunately, due to the lack of a large volume of labeled datasets, object detectors for SAR ship detection have developed slowly. To boost the development of object detectors in SAR images, a SAR dataset is constructed. This dataset labeled by SAR experts was created using 102 Chinese Gaofen-3 images and 108 Sentinel-1 images. It consists of 43,819 ship chips of 256 pixels in both range and azimuth. These ships mainly have distinct scales and backgrounds. Moreover, modified state-of-the-art object detectors from natural images are trained and can be used as baselines. Experimental results reveal that object detectors achieve higher mean average precision (mAP) on the test dataset and have high generalization performance on new SAR imagery without land-ocean segmentation, demonstrating the benefits of the dataset we constructed.

Polarimetry is one of the most promising types of remote sensing for improved characterization of atmospheric aerosol. Indeed, aerosol particles constitute a highly variable atmospheric component characterized by a large number of parameters describing particle sizes, morphologies (including shape and internal structure), absorption and scattering properties, amounts, horizontal and vertical distribution, etc. Reliable monitoring of all these parameters is very challenging, and therefore the aerosol effects on climate and environment are considered to be among the most uncertain factors in climate and environmental research. In this regard, observations that provide both the angular distribution of the scattered atmospheric radiation as well as its polarization state at multiple wavelengths covering the UV–SWIR spectral range carry substantial implicit information on the atmospheric composition. Therefore, high expectations in improving aerosol characterization are associated with detailed passive photopolarimetric observations. The critical need to use space-borne polarimetry for global accurate monitoring of detailed aerosol properties was first articulated in the late 1980s and early 1990s. By now, several orbital instruments have already provided polarization observations from space, and a number of advanced missions are scheduled for launch in the coming years by international and national space agencies. The first and most extensive record of polarimetric imagery was provided by POLDER-I, POLDER-II, and POLDER/PARASOL multi-angle multi-spectral polarization sensors. Polarimetric observations with the POLDER-like design intended for collecting extensive multi-angular multi-spectral measurements will be provided by several instruments, such as the MAI/TG-2, CAPI/TanSat, and DPC/GF-5 sensors recently launched by the Chinese Space Agency. Instruments such as the 3MI/MetOp-SG, MAIA, SpexOne and HARP2 on PACE, POSP, SMAC, PCF, DPC–Lidar, ScanPol and MSIP/Aerosol-UA, MAP/Copernicus CO2 Monitoring, etc. are planned to be launched by different space agencies in the coming decade. The concepts of these future instruments, their technical designs, and the accompanying algorithm development have been tested intensively and analyzed using diverse airborne prototypes. Certain polarimetric capabilities have also been implemented in such satellite sensors as GOME-2/MetOp and SGLI/GCOM-C. A number of aerosol retrieval products have been developed based on the available measurements and successfully used for different scientific applications. However, the completeness and accuracy of aerosol data operationally derived from polarimetry do not yet appear to have reached the accuracy levels implied by theoretical sensitivity studies that analyzed the potential information content of satellite polarimetry. As a result, the dataset provided by MODIS is still most frequently used by the scientific community, yet this sensor has neither polarimetric nor multi-angular capabilities. Admittedly polarimetric multi-angular observations are highly complex and have extra sensitivities to aerosol particle morphology, vertical variability of aerosol properties, polarization of surface reflectance, etc. As such, they necessitate state-of-the-art forward modeling based on first-principles physics which remains rare, and conventional retrieval approaches based on look-up tables turn out to be unsuitable to fully exploit the information implicit in the measurements. Several new-generation retrieval approaches have recently been proposed to address these challenges. These methods use improved forward modeling of atmospheric (polarized) radiances and implement a search in the continuous space of solutions using rigorous statistically optimized inversions. Such techniques provide more accurate retrievals of the main aerosol parameters such as aerosol optical thickness and yield additional parameters such as aerosol absorption. However, the operational implementation of advanced retrieval approaches generally requires a significant extra effort, and the forward-modeling part of such retrievals still needs to be substantially improved. Ground-based passive polarimetric measurements have also been evolving over the past decade. Although polarimetry helps improve aerosol characterization, especially of the fine aerosol mode, the operators of major observational networks such as AERONET remain reluctant to include polarimetric measurements as part of routine retrievals owing to their high complexity and notable increase in effort required to acquire and interpret polarization data. In addition to remote-sensing observations, polarimetric characteristics of aerosol scattering have been measured in situ as well as in the laboratory using polar nephelometers. Such measurements constitute direct observations of single scattering with no contributions from multiple scattering effects and therefore provide unique data for the validation of aerosol optical models and retrieval concepts. This article overviews the above-mentioned polarimetric observations, their history and expected developments, and the state of resulting aerosol products. It also discusses the main achievements and challenges in the exploitation of polarimetry for the improved characterization of atmospheric aerosols.

Core Ideas Heihe was the first basin‐scale integrated observatory network established in China. An intensive flux observation matrix experiment was conducted. New techniques, e.g., wireless sensor network, flux matrix, and airborne remote sensing, are used. The integrated observatory network is useful in land surface processes research. Research on land surface processes at the catchment scale has drawn much attention over the past few decades, and a number of watershed observatories have been established worldwide. The Heihe River Basin (HRB), which contains the second largest inland river in China, is an ideal natural field experimental area for investigation of land surface processes involving diverse landscapes and the coexistence of cold and arid regions. The Heihe Integrated Observatory Network was established in 2007. For long‐term observations, a hydrometeorological observatory, ecohydrological wireless sensor network, and satellite remote sensing are now in operation. In 2012, a multiscale observation experiment on evapotranspiration over heterogeneous land surfaces was conducted in the midstream region of the HRB, which included a flux observation matrix, wireless sensor network, airborne remote sensing, and synchronized ground measurements. Under an open data policy, the datasets have been publicly released following careful data processing and quality control. The outcomes highlight the integrated research on land surface processes in the HRB and include observed trends, scaling methods, high spatiotemporal resolution remote sensing products, and model–data integration in the HRB, all of which are helpful to other endorheic basins in the “Silk Road Economic Belt.” Henceforth, the goal of the Heihe Integrated Observatory Network is to develop an intelligent monitoring system that incorporates ground‐based observatory networks, unmanned aerial vehicles, and multi‐source satellites through the Internet of Things technology. Furthermore, biogeochemical processes observation will be improved, and the study of integrating ground observations, remote sensing, and large‐scale models will be promoted further.

In just the past five years, the field of Earth observation has progressed beyond the offerings of conventional space agency based platforms to include a plethora of sensing opportunities afforded by CubeSats, Unmanned Aerial Vehicles (UAVs), and smartphone technologies that are being embraced by both for-profit companies and individual researchers. Over the previous decades, space agency efforts have brought forth well-known and immensely useful satellites such as the Landsat series and the Gravity Research and Climate Experiment (GRACE) system, with costs typically on the order of one billion dollars per satellite and with concept-to-launch timelines on the order of two decades (for new missions). More recently, the proliferation of smartphones has helped to miniaturise sensors and energy requirements, facilitating advances in the use of CubeSats that can be launched by the dozens, while providing ultra-high (3-5 m) resolution sensing of the Earth on a daily basis. Start-up companies that did not exist five years ago now operate more satellites in orbit than any space agency, and at costs that are a mere fraction of the cost of traditional satellite missions. With these advances come new space-borne measurements, such as real-time high-definition video for tracking air pollution, storm-cell development, flood propagation, precipitation monitoring, or even for constructing digital surfaces using structure-from-motion techniques. Closer to the surface, measurements from small unmanned drones and tethered balloons have mapped snow depths, floods, and estimated evaporation at sub-meter resolutions, pushing back on spatio-temporal constraints and delivering new process insights. At ground level, precipitation has been measured using signal attenuation between antennae mounted on cell phone towers, while the proliferation of mobile devices has enabled citizen-scientists to catalogue photos of environmental conditions, estimate daily average temperatures from battery state, and sense other hydrologically important variables such as channel depths using commercially available wireless devices. Global internet access is being pursued via high altitude balloons, solar planes, and hundreds of planned satellite launches, providing a means to exploit the Internet of Things as an entirely new measurement domain. Such global access will enable real-time collection of data from billions of smartphones or from remote research platforms. This future will produce petabytes of data that can only be accessed via cloud storage and will require new analytical approaches to interpret. The extent to which today's hydrologic models can usefully ingest such massive data volumes is unclear. Nor is it clear whether this deluge of data will be usefully exploited, either because the measurements are superfluous, inconsistent, not accurate enough, or simply because we lack the capacity to process and analyse them. What is apparent is that the tools and techniques afforded by this array of novel and game-changing sensing platforms present our community with a unique opportunity to develop new insights that advance fundamental aspects of the hydrological sciences. To accomplish this will require more than just an application of the technology: in some cases, it will demand a radical rethink on how we utilise and exploit these new observing systems to enhance our understanding of the Earth and its linked processes.

A general framework for processing high and very-high resolution imagery in support of a Global Human Settlement Layer (GHSL) is presented together with a discussion on the results of the first operational test of the production workflow. The test involved the mapping of 24.3 million km <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">2</sup> of the Earth surface spread in four continents, corresponding to an estimated population of 1.3 billion people in 2010. The resolution of the input image data ranges from 0.5 to 10 meters, collected by a heterogeneous set of platforms including satellite SPOT (2 and 5), CBERS 2B, RapidEye (2 and 4), WorldView (1 and 2), GeoEye 1, QuickBird 2, Ikonos 2, and airborne sensors. Several imaging modes were tested including panchromatic, multispectral and pan-sharpened images. A new fully automatic image information extraction, generalization and mosaic workflow is presented that is based on multiscale textural and morphological image features extraction. New image feature compression and optimization are introduced, together with new learning and classification techniques allowing for the processing of HR/VHR image data using low-resolution thematic layers as reference. A new systematic approach for quality control and validation allowing global spatial and thematic consistency checking is proposed and applied. The quality of the results are discussed by sensor, band, resolution, and eco-regions. Critical points, lessons learned and next steps are highlighted.

Restoring disturbed and over-exploited ecosystems is important to mitigate human pressures on natural ecosystems. China has launched an ambitious national ecosystem restoration program called Grain to Green Program (GTGP) over the last decade. By using remote sensing techniques and ecosystem modelling, we quantitatively evaluated the changes in ecosystem carbon sequestration since China's GTGP program during period of 2000-2008. It was found the NPP and NEP in this region had steadily increased after the initiative of the GTGP program, and a total of 96.1 Tg of additional carbon had been sequestered during that period. Changes in soil carbon storage were lagged behind and thus insignificant over the period, but was expected to follow in the coming decades. As a result, the Loess Plateau ecosystem had shifted from a net carbon source in 2000 to a net carbon sink in 2008. The carbon sequestration efficiency was constrained by precipitation, and appropriate choices of restoration types (trees, shrubs, and grasses) in accordance to local climate are critical for achieving the best benefit/cost efficiency.