Rogier de Jong

Shifts in Global Vegetation Activity Trends

Vegetation belongs to the components of the Earth surface, which are most extensively studied using historic and present satellite records. Recently, these records exceeded a 30-year time span composed of preprocessed fortnightly observations (1981–2011). The existence of monotonic changes and trend shifts present in such records has previously been demonstrated. However, information on timing and type of such trend shifts was lacking at global scale. In this work, we detected major shifts in vegetation activity trends and their associated type (either interruptions or reversals) and timing. It appeared that the biospheric trend shifts have, over time, increased in frequency, confirming recent findings of increased turnover rates in vegetated areas. Signs of greening-to-browning reversals around the millennium transition were found in many regions (Patagonia, the Sahel, northern Kazakhstan, among others), as well as negative interruptions—“setbacks”—in greening trends (southern Africa, India, Asia Minor, among others). A minority (26%) of all significant trends appeared monotonic.

No growth stimulation of Canada’s boreal forest under half-century of combined warming and CO2 fertilization

Significance Limited knowledge about the mechanistic drivers of forest growth and responses to environmental changes creates uncertainties about the future role of circumpolar boreal forests in the global carbon cycle. Here, we use newly acquired tree-ring data from Canada’s National Forest Inventory to determine the growth response of the boreal forest to environmental changes. We find no consistent boreal-wide growth response over the past 60 y across Canada. However, some southwestern and southeastern forests experienced a growth enhancement, and some regions such as the northwestern and maritime areas experienced a growth depression. Growth–climate relationships bring evidence of an intensification of the impacts of hydroclimatic variability on growth late in the 20th century, in parallel with the rapid rise of summer temperature. Considerable evidence exists that current global temperatures are higher than at any time during the past millennium. However, the long-term impacts of rising temperatures and associated shifts in the hydrological cycle on the productivity of ecosystems remain poorly understood for mid to high northern latitudes. Here, we quantify species-specific spatiotemporal variability in terrestrial aboveground biomass stem growth across Canada’s boreal forests from 1950 to the present. We use 873 newly developed tree-ring chronologies from Canada’s National Forest Inventory, representing an unprecedented degree of sampling standardization for a large-scale dendrochronological study. We find significant regional- and species-related trends in growth, but the positive and negative trends compensate each other to yield no strong overall trend in forest growth when averaged across the Canadian boreal forest. The spatial patterns of growth trends identified in our analysis were to some extent coherent with trends estimated by remote sensing, but there are wide areas where remote-sensing information did not match the forest growth trends. Quantifications of tree growth variability as a function of climate factors and atmospheric CO2 concentration reveal strong negative temperature and positive moisture controls on spatial patterns of tree growth rates, emphasizing the ecological sensitivity to regime shifts in the hydrological cycle. An enhanced dependence of forest growth on soil moisture during the late-20th century coincides with a rapid rise in summer temperatures and occurs despite potential compensating effects from increased atmospheric CO2 concentration.

Creating Multi-Temporal Composites of Airborne Imaging Spectroscopy Data in Support of Digital Soil Mapping

An increasing demand for full spatio-temporal coverage of soil information drives the growing use of soil spectroscopy. Soil spectroscopy application performed under laboratory conditions or in-field studies in semi-arid areas have shown promising results. However, when acquiring data in temperate zones, limitations by vegetation-free coverage, variation in soil moisture and management are driving coherent spatio-temporal data collection. This study explores the use of multi-temporal imaging spectroscopy data to increase the total mapping area of bare soils in a heterogeneous agricultural landscape. Spectrally and spatially high-resolution data from the Airborne Prism Experiment (APEX) were collected in September 2013, April 2014 and April 2015. Bare soils in all acquisitions were identified. To eliminate short-term differences in soil moisture and soil surface roughness, the empirical line method was used to calibrate the reflectance values of the singular images (2013 and 2015) towards the singular image with most bare soil pixels (2014). Difference indicators show that the calibration was successful (decrease in root mean square difference and angle difference, increase in R2 and gain and offset close to one and zero). Finally, the multi-temporal composite image contained more than double the amount of bare soil pixels as compared to a singular acquisition. Summary statistics show that reflectance values of the multi-temporal composite approximate the single image data of 2014 (mean and standard deviation of 2014: 24.2 ± 8.9 vs. 24.0 ± 9.5 for the multi-temporal composite of 2013, 2014 and 2015). This indicates that global differences in soil moisture and land management have been corrected for. As a result, an improved spatial representation of soil parameters can be retrieved from the composite data. Spatial distribution of the correction factors and analysis of the spatial variability of all images, however, indicate that non-linear, short-term differences like variation in soil moisture and land management largely influence the result of the multi-temporal composite. Quantification and attribution of those factors will be required in the future to allow correcting for them.

Comparison of vegetation phenological metrics extracted from GIMMS NDVIg and MERIS MTCI data sets over China

Evaluating vegetation phenology is crucial for a better understanding of the effects of climate change on the terrestrial ecosystem. The scientific community has used various vegetation index data sets from different sensors to quantify vegetation phenology from regional to global scales. The normalized difference vegetation index (NDVI) related to photosynthetic activities is the most widely used index. Recently, a number of published articles have used the Medium Resolution Imaging Spectrometer (MERIS) terrestrial chlorophyll index (MTCI) to measure vegetation phenology. MTCI can closely represent the red-edge position (REP). Unlike NDVI, MTCI is more sensitive to high values of chlorophyll content. However, the consistency of vegetation phenological metrics derived from MTCI and NDVI needs to be further explored. This study compared two phenological metrics, i.e. onset of greenness (OG) and end of senescence (ES), extracted from MERIS MTCI data and Advanced Very High Resolution Radiometer (AVHRR) Global Inventory Modeling and Mapping Studies (GIMMS) first generation NDVI (NDVIg) data, which has the longest time records, at nine regions in China from 2003 to 2006. The results showed that the differences of OG and ES vary between different vegetation types, regions, and years, although both NDVI and MTCI time series capture the growth patterns well for most vegetation types. Compared to ES, the OG estimates are more consistent. NDVI yields in general later ES estimates than MTCI.

Spatial Differentiation of Arable Land and Permanent Grassland to Improve a Land Management Model for Nutrient Balancing

Agroecosystems play an important role in providing economic and ecosystem services, which directly impact society. Inappropriate land use and unsustainable agricultural management with associated nutrient cycles can jeopardize important soil functions such as food production, livestock feeding, and conservation of biodiversity. The objective of this study was to integrate remotely sensed land cover information into a regional land management model (LMM) to improve the assessment of spatially explicit nutrient balances for agroecosystems. Remotely sensed data and an optimized parameter set contributed to an improved LMM output, allowing for a better land allocation within the model. The best input parameter combination was based on two different land cover classifications with overall accuracies of 98%, improving the land allocation performance compared with using nonspatially explicit input. We conclude that the combined use of remote sensing data and the LMM has the potential to provide valuable guidance for farm practices. It further helps to generate a spatial description of farm-level nutrient balance, a crucial ability when choosing policy options related to sustainable management of agricultural soils.

Barest Pixel Composite for Agricultural Areas Using Landsat Time Series

Many soil remote sensing applications rely on narrow-band observations to exploit molecular absorption features. However, broadband sensors are invaluable for soil surveying, agriculture, land management and mineral exploration, amongst others. These sensors provide denser time series compared to high-resolution airborne imaging spectrometers and hold the potential of increasing the observable bare-soil area at the cost of spectral detail. The wealth of data coming along with these applications can be handled using cloud-based processing platforms such as Earth Engine. We present a method for identifying the least-vegetated observation, or so called barest pixel, in a dense time series between January 1985 and March 2017, based on Landsat 5, 7 and 8 observations. We derived a Barest Pixel Composite and Bare Soil Composite for the agricultural area of the Swiss Plateau. We analysed the available data over time and concluded that about five years of Landsat data were needed for a full-coverage composite (90% of the maximum bare soil area). Using the Swiss harmonised soil data, we derived soil properties (sand, silt, clay, and soil organic matter percentages) and discuss the contribution of these soil property maps to existing conventional and digital soil maps. Both products demonstrate the substantial potential of Landsat time series for digital soil mapping, as well as for land management applications and policy making.

Altitude‐dependent influence of snow cover on alpine land surface phenology

Snow cover impacts alpine land surface phenology in various ways but our knowledge about the effect of snow cover on alpine land surface phenology is still limited. We studied this relationship in the European Alps using satellite-derived metrics of Snow Cover Phenology (SCP), namely First Snow Fall, Last Snow Day and Snow Cover Duration (FSF, LSD and SCD, respectively), in combination with Land Surface Phenology (LSP), namely Start Of Season, End Of Season and Length Of Season (SOS, EOS and LOS, respectively) for the period of 2003–2014. We tested the dependency of inter-annual differences (Δ) of SCP and LSP metrics with altitude (up to 3000 meter above sea level (m a.s.l.)) for seven natural vegetation types, four main climatic subregions and four terrain expositions. We found that 25.3% of all pixels showed significant (p < 0.05) correlation between ΔSCD and ΔSOS and 15.3% between ΔSCD and ΔLOS across the entire study area. Correlations between ΔSCD and ΔSOS as well as ΔSCD and ΔLOS are more pronounced in the northern subregions of the Alps, at high altitudes, and on north- and west-facing terrain – or more generally, in regions with longer SCD. We conclude that snow cover has a greater effect on alpine phenology at higher than at lower altitudes, which may be attributed to the coupled influence of snow cover with underground conditions and air temperature. Alpine ecosystems may therefore be particularly sensitive to future change of snow cover at high altitudes under climate warming scenarios.

Predicting missing values in spatio-temporal satellite data

Continuous, consistent, and long time-series from remote sensing are essential to monitoring changes on Earth’s surface. However, analyzing such data sets is often challenging due to missing values introduced by cloud cover, missing orbits, sensor geometry artifacts, and so on. We propose a new and accurate spatio-temporal prediction method to replace missing values in remote sensing data sets. The method exploits the spatial coherence and temporal seasonal regularity that are inherent in many data sets. The key parts of the method are: 1) the adaptively chosen spatio-temporal subsets around missing values; 2) the ranking of images within the subsets based on a scoring algorithm; 3) the estimation of empirical quantiles characterizing the missing values; and 4) the prediction of missing values through quantile regression. One advantage of quantile regression is the robustness to outliers, which enables more accurate parameter retrieval in the analysis of remote sensing data sets. In addition, we provide bootstrap-based quantification of prediction uncertainties. The proposed prediction method was applied to a Normalized Difference Vegetation Index data set from the Moderate Resolution Imaging Spectroradiometer and assessed with realistic test data sets featuring between 20% and 50% missing values. Validation against established methods showed that the proposed method has a good performance in terms of the root-mean-squared prediction error and significantly outperforms its competitors. This paper is accompanied by the open-source R package gapfill, which provides a flexible, fast, and ready-to-use implementation of the method.

Relative Influence of Timing and Accumulation of Snow on Alpine Land Surface Phenology

Timing and accumulation of snow are among the most important phenomena influencing land surface phenology in mountainous ecosystems. However, our knowledge on their influence on alpine land surface phenology is still limited, and much remains unclear as to which snow metrics are most relevant for studying this interaction. In this study, we analyzed five snow and phenology metrics, namely, timing (snow cover duration (SCD) and last snow day), accumulation of snow (mean snow water equivalent, SWEm), and mountain land surface phenology (start of season and length of season) in the Swiss Alps during the period 2003–2014. We examined elevational and regional variations in the relationships between snow and alpine land surface phenology metrics using multiple linear regression and relative weight analyses and subsequently identified the snow metrics that showed strongest associations with variations in alpine land surface phenology of natural vegetation types.We found that the relationships between snow and phenology metrics were pronounced in high-elevational regions and alpine natural grassland and sparsely vegetated areas. Start of season was influenced primarily by SCD, secondarily by SWEm, while length of season was equally affected by SCD and SWEm across different elevational bands. We conclude that SCD plays the most significant role compared to other snow metrics. Future variations of snow cover and accumulation are likely to influence alpine ecosystems, for instance, their species composition due to changes in the potential growing season. Also, their spatial distribution may change as a response to the new environmental conditions if these prove persistent.

Comparative study of three satellite image time-series decomposition methods for vegetation change detection

ABSTRACT Satellite image time-series (SITS) methods have contributed notably to detection of global change over the last decades, for instance by tracking vegetation changes. Compared with multi-temporal change detection methods, temporally highly resolved SITS methods provide more information in a single analysis, for instance on the type and consistency of change. In particular, SITS decomposition methods show a great potential in extracting various components from non-stationary time series, which allows for an improved interpretation of the temporal variability. Even though many case studies have applied SITS decomposition methods, a systematic comparison of common algorithms is still missing. In this study, the seasonal trend loess (STL), breaks for additive season and trend (BFAST) and multi-resolution analysis-wavelet transform (MRA-WT) were explored in order to evaluate their performance in modelling, monitoring and detecting land-cover changes with pronounced seasonal variations from simulated normal difference vegetation index time series. The selected methods have all proven their ability to characterize the non-stationary vegetation dynamics along with different physical processes driving the vegetation dynamics. Our results indicated that BFAST is the most accurate method for the examined simulated dataset in terms of RMSE, whereas MRA-WT showed a great potential for the extraction of multi-level vegetation dynamics. Considering the computational efficiency, both STL and MRA-WT outperformed BFAST.