Terrestrial remote sensing data products retrieved from radiometric measurements in the optical and thermal infrared spectrum such as vegetation spectral indices can be heavily contaminated by atmospheric conditions, including cloud and aerosol layers. This contamination results in gaps or noisy observations. The harmonic analysis of time series (HANTS) has been widely used for time series reconstruction of remote sensing imagery in recent decades. To use HANTS model, a series of parameters, such as number of frequencies (NF), fitting error tolerance (FET), degree of over-determinedness (DoD), and regularization factor (Delta), need to be defined by users. These parameters provide flexibilities, but also make it difficult for non-expert users to determine appropriate settings for specific applications. This study systematically evaluated the reconstruction performance of the model under different parameter setting scenarios by simulating pixel-wise reference and noisy NDVI time series. The results of these numerical experiments were further used to identify optimal settings and improve global NDVI reconstruction performance. The results suggested optimal settings for different areas (local optimization). If a user opts to use unique settings for global reconstruction, the setting NF = 4, FET = 0.05, DoD = 5, and Delta = 0.5 can produce the best performance across all setting scenarios (global optimization). In addition, several internal improvements, such as dynamic weighting scheme, polynomial and inter-annual harmonic components, and ancillary attributes of input data can be used to further improve the performance of reconstruction. With these results, future non-expert users can easily determine appropriate settings of HANTS for specific applications in different regions.@en

Spatiotemporal residual noise in terrestrial earth observation products, often caused by unfavorable atmospheric conditions, impedes their broad applications. Most users prefer to use gap-filled remote sensing products with time series reconstruction (TSR) algorithms. Applying currently available implementations of TSR to large-volume datasets is time-consuming and challenging for non-professional users with limited computation or storage resources. This study introduces a new open-source software package entitled ‘HANTS-GEE’ that implements a well-known and robust TSR algorithm, i.e. Harmonic ANalysis of Time Series (HANTS), on the Google Earth Engine (GEE) platform for scalable reconstruction of terrestrial earth observation data. Reconstruction tasks can be conducted on user-defined spatiotemporal extents when raw datasets are available on GEE. According to site-based and regional-based case evaluation, the new tool can effectively eliminate cloud contamination in the time series of earth observation data. Compared with traditional PC-based HANTS implementation, the HANTS-GEE provides quite consistent reconstruction results for most terrestrial vegetated sites. The HANTS-GEE can provide scalable reconstruction services with accelerated processing speed and reduced internet data transmission volume, promoting algorithm usage by much broader user communities. To our knowledge, the software package is the first tool to support full-stack TSR processing for popular open-access satellite sensors on cloud platforms.@en

The time lag between anomalies in precipitation and vegetation activity plays a critical role in early drought detection as agricultural droughts are caused by precipitation shortages. The aim of this study is to explore a new approach to estimate the time lag between a forcing (precipitation) and a response (NDVI) signal in the frequency domain by applying cross-spectral analysis. We prepared anomaly time series of image data on TRMM3B42 precipitation (accumulated over antecedent durations of 10, 60, and 150 days) and NDVI, reconstructed and interpolated MOD13A2 and MYD13A2 to daily interval using a Fourier series method to model time series affected by gaps and outliers (iHANTS) for a dry and a wet year in a drought-prone area in the northeast region of China. Then, the cross-spectral analysis was applied pixel-wise and only the phase lag of the annual component of the forcing and response signal was extracted. The 10-day antecedent precipitation was retained as the best representation of forcing. The estimated phase lag was interpreted using maps of land cover and of available soil water-holding capacity and applied to investigate the difference in phenology responses between a wet and dry year. In both the wet and dry year, we measured consistent phase lags across land cover types. In the wet year with above-average precipitation, the phase lag was rather similar for all land cover types, i.e., 7.6 days for closed to open grassland and 14.5 days for open needle-leaved deciduous or evergreen forest. In the dry year, the phase lag increased by 7.0 days on average, but with specific response signals for the different land cover types. Interpreting the phase lag against the soil water-holding capacity, we observed a slightly higher phase lag in the dry year for soils with a higher water-holding capacity. The accuracy of the estimated phase lag was assessed through Monte Carlo simulations and presented reliable estimates for the annual component.@en

The regional surface soil heat flux (G0) estimation is very important for the large-scale land surface process modeling. However, most of the regional G0 estimation methods are based on the empirical relationship between G0 and the net radiation flux. A physical model based on harmonic analysis was improved (referred to as "HM model") and applied over the Heihe River Basin northwest China with multiple remote sensing data, e.g., FY-2C, AMSR-E, and MODIS, and soil map data. The sensitivity analysis of the model was studied as well. The results show that the improved model describes the variation of G0 well. Land surface temperature (LST) and thermal inertia (Γ) are the two key input variables to the HM model. Compared with in situ G0, there are some differences, mainly due to the differences between remote-sensed LST and the in situ LST. The sensitivity analysis shows that the errors from-7 to-0.5K in LST amplitude and from-300 to 300J m-2 K-1 s-0.5 in Γ will cause about 20% errors, which are acceptable for G0 estimation.@en

Observing and understanding changes in Africa is a hotspot in global ecological environmental research since the early 1970s. As possible causes of environmental degradation, frequent droughts and human activities attracted wide attention. Remote sensing of nighttime light provides an effective way to map human activities and assess their intensity. To identify settlements more effectively, this study focused on nighttime light in the northern Equatorial Africa and Sahel settlements to propose a new method, namely, the patches filtering method (PFM) to identify nighttime lights related to settlements from the National Polar-orbiting Partnership Visible Infrared Imaging Radiometer Suite (NPP-VIIRS) monthly nighttime light data by separating signal components induced by biomass burning, thereby generating a new annual image in 2016. The results show that PFM is useful for improving the quality of NPP-VIIRS monthly nighttime light data. Settlement lights were effectively separated from biomass burning lights, in addition to capturing the seasonality of biomass burning. We show that the new 2016 nighttime light image can very effectively identify even small settlements, notwithstanding their fragmentation and unstable power supply. We compared the image with earlier NPP-VIIRS annual nighttime light data from the National Oceanic and Atmospheric Administration (NOAA) National Center for Environmental Information (NCEI) for 2016 and the Sentinel-2 prototype Land Cover 20 m 2016 map of Africa released by the European Space Agency (ESA-S2-AFRICA-LC20). We found that the new annual nighttime light data performed best among the three datasets in capturing settlements, with a high recognition rate of 61.8%, and absolute superiority for settlements of 2.5 square kilometers or less. This shows that the method separates biomass burning signals very effectively, while retaining the relatively stable, although dim, lights of small settlements. The new 2016 annual image demonstrates good performance in identifying human settlements in sparsely populated areas toward a better understanding of human activities.@en

Satellite-based Earth Observation systems archived a variety of vegetation products during the last 50 years, which can reveal regional to global ecosystem dynamics across diverse spatiotemporal scales. The anomaly metrics such as Vegetation Condition Index (VCI) defined by comparing the current vegetation growth condition to historical average status based on long-term EO-based vegetation products were widely used to delineate abnormal vegetation variation exerted by either climatic or anthropogenic factors (e.g., droughts, wildfires). However, currently available long-term vegetation products may differ from each other in terms of sensors (observational platform or spectral bands), bio-physical definitions (e.g., NDVI, EVI, LAI, and VOD), spatiotemporal resolution, as well as the time-spans, which results in inconsistency across these vegetation products. Taking the VCI as an example, this study evaluated the uncertainty of vegetation anomalies detected based on different vegetation products over the middle reach of the Yangtze River by explicitly considering the effect of sensors, biophysical definitions, and time-spans. The preliminary results showed that VCI derived from NDVI products from different sensors (AVHRR vs. MODIS) induced significant inconsistent anomalies over most landscapes. The differences resulting from products with different biophysical definitions (NDVI vs. EVI, LAI, and VOD) are much lower than those from different sensors but still significant over specific areas. As for the time-spans, the 20-year NDVI based VCI presented a considerable reduction in variance over the study area on average compared to VCI calculated based on 5-year NDVI. In summary, caution should be taken when applying EO-based vegetation products for vegetation anomalies mapping, especially for quantitative assessment.@en

Rainfall is a key driver of terrestrial vegetation. Clarifying the response mechanism of vegetation to rainfall can advance the understanding of expected changes in ecosystems under projected rainfall scenarios. Besides the rainfall amount over a period of time, the frequency, duration and intensity are of importance in driving ecosystem processes. The pulsating nature of rainfall forcing and of vegetation response applies particularly to arid and semi – arid regions and led to conceptualize the “pulse-reserve” paradigm, which we used to explain the approach we propose in this study. We introduce the Transfer Function Analysis (TFA) method in the frequency domain to capture the response of vegetation to rainfall at multiple temporal scales. The TFA method determines coherence, gain, and phase to characterize the existence, strength, and time-lag, respectively, of the vegetation response to rainfall at different temporal scales. Specifically, the coherence measures the existence of the response and only with a significant response are the gain and phase values significant and valuable for further analysis. The gain measures the strength of the relationship between fluctuations in vegetation growth and fluctuations in rainfall, while the phase value (i.e. the time-lag) measures how fluctuations in vegetation growth lag (or lead) fluctuations in rainfall. The TFA method was applied to the 34-years (1982–2015) NDVI3g and CHIRPS precipitation dataset in the Sahel-Sudano-Guinean region (20°W ~ 60°E, 0–25°E). The Sahelian zone was characterized by a significant vegetation response to rainfall across all inter- and intra- annual time-scales, while the Sudano-Guinean zone was dominated by significant response at annual or 6-month scales. The negative phase lag indicated that rainfall variation normally led NDVI change for most areas and across timescales. However, a positive phase observed in part of the tropical rainforest area indicated that NDVI changes led rainfall variations, which may be caused by the strong vegetation-rainfall feedback through recycling of precipitation by evapotranspiration. In summary, these results suggested that the TFA method is a powerful tool to quantify the vegetation-rainfall response regime across a range of timescales, as conceptualized by the “pulse-reserve” paradigm. Unraveling the response of fluctuations in vegetation growth to separate components of the forcing by precipitation might improve our understanding of environmental change in the past decades in the Sahel - Sudan - Guinean region.@en

The monthly nighttime light data observed by the Visible Infrared Imaging Radiometer Suite (VIIRS) sensor onboard the National Polar-orbiting Partnership (NPP) satellite are useful to assess human activities. However, noise in the VIIRS data and unstable power supply in many areas often limit its application in detecting anthropic activities. This study aims to develop a new method (i.e., Patches Denoising Method, PDM) to remove noise in the VIIRS monthly nighttime light data towards better information on anthropic activities. The new synthesized data were examined in the Sahel region where the intensity of nighttime light is often affected by unstable power supply. The results showed that the new data could capture more settlements correctly compared with the annual data released by the National Centers for Environmental Information (NCEI) of the National Oceanic and Atmospheric Administration (NOAA). This suggests that the Patches De-noising Method is a useful way to remove noise in the NPP-VIIRS monthly nighttime light data.@en

This paper evaluated the accuracy of multiple satellite-based precipitation products including the tropical rainfall measuring mission multisatellite precipitation analysis (TMPA) (TMPA 3B42RT and TMPA 3B42 version 7) and the Climate Prediction Center MORPHing technique (CMORPH) (CMORPH RAW and CMORPH BLD version 1.0) datasets and investigated the impact of the accuracy and temporal coverage of these data products on the reliability of the standardized precipitation index (SPI) estimates. The satellite-based SPI was compared with the SPI estimate using in situ precipitation observations from 2221 meteorological observation sites across China from 1998 to 2014. The SPI values calculated from the products calibrated with rain gauge measurements (TMPA 3B42 and CMORPH BLD) are generally more consistent with the SPI obtained with in situ measurements than those obtained using noncalibrated products (TMPA 3B42RT and CMORPH RAW products). The short data record of satellite precipitation data products is not the primary source of large errors in the SPI estimates, suggesting that the SPI estimate using satellite precipitation data products can be applied to drought assessment and monitoring. Satellite-based SPI estimates are more accurate in eastern China than in western China because of larger uncertainties in precipitation retrievals in western China, characterized by arid and semiarid climate conditions and complex landscapes. The satellite-based SPI can capture typical drought events throughout China, with the limitation that it is based on precipitation only and that different durations of antecedent precipitation are only suitable for specific drought conditions.@en

This paper evaluated the accuracy of multiple satellite-based precipitation products including the tropical rainfall measuring mission multisatellite precipitation analysis (TMPA) (TMPA 3B42RT and TMPA 3B42 version 7) and the Climate Prediction Center MORPHing technique (CMORPH) (CMORPH RAW and CMORPH BLD version 1.0) datasets and investigated the impact of the accuracy and temporal coverage of these data products on the reliability of the standardized precipitation index (SPI) estimates. The satellite-based SPI was compared with the SPI estimate using in situ precipitation observations from 2221 meteorological observation sites across China from 1998 to 2014. The SPI values calculated from the products calibrated with rain gauge measurements (TMPA 3B42 and CMORPH BLD) are generally more consistent with the SPI obtained with in situ measurements than those obtained using noncalibrated products (TMPA 3B42RT and CMORPH RAW products). The short data record of satellite precipitation data products is not the primary source of large errors in the SPI estimates, suggesting that the SPI estimate using satellite precipitation data products can be applied to drought assessment and monitoring. Satellite-based SPI estimates are more accurate in eastern China than in western China because of larger uncertainties in precipitation retrievals in western China, characterized by arid and semiarid climate conditions and complex landscapes. The satellite-based SPI can capture typical drought events throughout China, with the limitation that it is based on precipitation only and that different durations of antecedent precipitation are only suitable for specific drought conditions.@en

The Pan-Eurasian and African Continents are characterized by large ranges of climates varying from humid, semi-humid, semi-arid and arid regions, and great challenges exist in water allocation for different sectors that related to water resource and food security, which depends strongly on the water use information. Quantitative information on water use is also important to understand the effectiveness of water allocation and further to prevent from water stress resulted by drought in water-scarce regions. Explosive development of satellite remote sensing observations provide great chance to provide useful spatiotemporal information for quantifying the water use at regional to global scales. In this paper, a process-based model ETMonitor was used in combination with biophysical and hydrological parameters retrieved from earth observations to estimate the actual evapotranspiration, i.e. the agricultural and ecological water use. The total water use is also partitioned into beneficial part, e.g. plant transpiration, and non-beneficial part, e.g. soil evaporation and canopy rainfall interception, according to the water accounting framework. The estimated water use show good agreements with the ground observation, indicating the ability of ETMonitor for global and continental scale water use estimation. The spatial and temporal patterns of the water use in the Pan-Eurasian and African Continents were further analysed, while large spatial variation of water use was convinced. Current study also highlights the great capability of satellite observations in studying the regional water resource and continental water cycle.@en