MODIS vegetation products as proxies of photosynthetic potential along a gradient of meteorologically and biologically driven ecosystem productivity

© 2016 Author(s). A direct relationship between gross ecosystem productivity (GEP) estimated by the eddy covariance (EC) method and Moderate Resolution Imaging Spectroradiometer (MODIS) vegetation indices (VIs) has been observed in many temperate and tropical ecosystems. However, in Australian evergreen forests, and particularly sclerophyll and temperate woodlands, MODIS VIs do not capture seasonality of GEP. In this study, we re-evaluate the connection between satellite and flux tower data at four contrasting Australian ecosystems, through comparisons of GEP and four measures of photosynthetic potential, derived via parameterization of the light response curve: ecosystem light use efficiency (LUE), photosynthetic capacity (Pc), GEP at saturation (GEPsat), and quantum yield (α) with MODIS vegetation satellite products, including VIs, gross primary productivity (GPPMOD) leaf area index (LAIMOD), and fraction of photosynthetic active radiation (fPARMOD). We found that satellite-derived biophysical products constitute a measurement of ecosystem structure (e.g. leaf area index-quantity of leaves) and function (e.g. leaf level photosynthetic assimilation capacity-quality of leaves), rather than GEP. Our results show that in primarily meteorological-driven (e.g. photosynthetic active radiation, air temperature, and/or precipitation) and relatively aseasonal ecosystems (e.g. evergreen wet sclerophyll forests), there were no statistically significant relationships between GEP and satellite-derived measures of greenness. In contrast, for phenology-driven ecosystems (e.g. tropical savannas), changes in the vegetation status drove GEP, and tower-based measurements of photosynthetic activity were best represented by VIs. We observed the highest correlations between MODIS products and GEP in locations where key meteorological variables and vegetation phenology were synchronous (e.g. semi-arid Acacia woodlands) and low correlation at locations where they were asynchronous (e.g. Mediterranean ecosystems). However, we found a statistical significant relationship between the seasonal measures of photosynthetic potential (Pc and LUE) and VIs, where each ecosystem aligns along a continuum; we emphasize here that knowledge of the conditions in which flux tower measurements and VIs or other remote sensing products converge greatly advances our understanding of the mechanisms driving the carbon cycle (phenology and climate drivers) and provides an ecological basis for interpretation of satellite-derived measures of greenness

Remote sensing environmental change in southern African savannahs : a case study of Namibia

Savannah biomes cover a fifth of Earth’s surface, harbour many of the world’s most iconic species and most of its livestock and rangeland, while sustaining the livelihoods of an important proportion of its human population. They provide essential ecosystem services and functions, ranging from forest, grazing and water resources, to global climate regulation and carbon sequestration. However, savannahs are highly sensitive to human activities and climate change. Across sub-Saharan Africa, climatic shifts, destructive wars and increasing anthropogenic disturbances in the form of agricultural intensification and urbanization, have resulted in widespread land degradation and loss of ecosystem services. Yet, these threatened ecosystems are some of the least studied or protected, and hence should be given high conservation priority. Importantly, the scale of land degradation has not been fully explored, thereby comprising an important knowledge gap in our understanding of ecosystem services and processes, and effectively impeding conservation and management of these biodiversity hotspots. The primary drivers of land degradation include deforestation, triggered by the increasing need for urban and arable land, and concurrently, shrub encroachment, a process in which the herbaceous layer, a defining characteristic of savannahs, is replaced with hardy shrubs. These processes have significant repercussions on ecosystem service provision, both locally and globally, although the extents, drivers and impacts of either remain poorly quantified and understood. Additionally, regional aridification anticipated under climate change, will lead to important shifts in vegetation composition, amplified warming and reduced carbon sequestration. Together with a growing human population, these processes are expected to compound the risk of land degradation, thus further impacting key ecosystem services. Namibia is undergoing significant environmental and socio-economic changes. The most pervasive change processes affecting its savannahs are deforestation, degradation and shrub encroachment. Yet, the extent and drivers of such change processes are not comprehensively quantified, nor are the implications for rural livelihoods, sustainable land management, the carbon cycle, climate and conservation fully explored. This is partly due to the complexities of mapping vegetation changes with satellite data in savannahs. They are naturally spatially and temporally variable owing to erratic rainfall, divergent plant functional type phenologies and extensive anthropogenic impacts such as fire and grazing. Accordingly, this thesis aims to (i) quantify distinct vegetation change processes across Namibia, and (ii) develop methodologies to overcome limitations inherent in savannah mapping. Multi-sensor satellite data spanning a range of spatial, temporal and spectral resolutions are integrated with field datasets to achieve these aims, which are addressed in four journal articles. Chapters 1 and 2 are introductory. Chapter 3 exploits the Landsat archive to track changes in land cover classes over five decades throughout the Namibian Kalahari woodlands. The approach addresses issues implicit in change detection of savannahs by capturing the distinct phenological phases of woody vegetation and integrating multi-sensor, multi-source data. Vegetation extent was found to have decreased due to urbanization and small-scale arable farming. An assessment of the limitations leads to Chapter 4, which elaborates on the previous chapter by quantifying aboveground biomass changes associated with deforestation and shrub encroachment. The approach centres on fusing multiple satellite datasets, each acting as a proxy for distinct vegetation properties, with calibration/validation data consisting of concurrent field and LiDAR measurements. Biomass losses predominate, demonstrating the contribution of land management to ecosystem carbon changes. To identify whether biomass is declining across the country, Chapter 5 focuses on regional, moderate spatial resolution time-series analyses. Phenological metrics extracted from MODIS data are used to model observed fractional woody vegetation cover, a proxy for biomass. Trends in modelled fractional woody cover are then evaluated in relation to the predominant land-uses and precipitation. Negative trends slightly outweighed positive trends, with decreases arising largely in protected, urban and communal areas. Since precipitation is a fundamental control on vegetation, Chapter 6 investigates its relation to NDVI, by assessing to what extent observed trends in vegetation cover are driven by rainfall. NDVI is modelled as a function of precipitation, with residuals assumed to describe the fraction of NDVI not explained by rainfall. Mean annual rainfall and rainfall amplitude show a positive trend, although extensive “greening” is unrelated to rainfall. NDVI amplitude, used as a proxy for vegetation density, indicates a widespread shift to a denser condition. In Chapter 7, trend analysis is applied to a Landsat time-series to overcome spatial and temporal limitations characteristic of the previous approaches. Results, together with those of the previous chapters, are synthesized and a synopsis of the main findings is presented. Vegetation loss is predominantly caused by demand for urban and arable land. Greening trends are attributed to shrub encroachment and to a lesser extent conservation laws, agroforestry and rangeland management, with precipitation presenting little influence. Despite prevalent greening, degradation processes associated with shrub encroachment, including soil erosion, are likely to be widespread. Deforestation occurs locally while shrub encroachment occurs regionally. This thesis successfully integrates multi-source data to map, measure and monitor distinct change processes across scales

Estimating aboveground woody biomass change in Kalahari woodland: combining field, radar, and optical data sets

Maps that accurately quantify aboveground vegetation biomass (AGB) are essential for ecosystem monitoring and conservation. Throughout Namibia, four vegetation change processes are widespread, namely, deforestation, woodland degradation, the encroachment of the herbaceous and grassy layers by woody strata (woody thickening), and woodland regrowth. All of these vegetation change processes affect a range of key ecosystem services, yet their spatial and temporal dynamics and contributions to AGB change remain poorly understood. This study quantifies AGB associated with the different vegetation change processes over an 8-year period, for a region of Kalahari woodland savannah in northern Namibia. Using data from 101 forest inventory plots collected during two field campaigns (2014–2015), we model AGB as a function of the Advanced Land Observing Satellite Phased Array L-band synthetic aperture radar (PALSAR and PALSAR-2) and dry season Landsat vegetation index composites, for two periods (2007 and 2015). Differences in AGB between 2007 and 2015 were assessed and validated using independent data, and changes in AGB for the main vegetation processes are quantified for the whole study area (75,501 km2). We find that woodland degradation and woody thickening contributed a change in AGB of −14.3 and 2.5 Tg over 14% and 3.5% of the study area, respectively. Deforestation and regrowth contributed a smaller portion of AGB change, i.e. −1.9 and 0.2 Tg over 1.3% and 0.2% of the study area, respectively

A new global fAPAR and LAI dataset derived from optimal albedo estimates: comparison with MODIS products

We present the first comparison between new fAPAR and LAI products derived from the GlobAlbedo dataset and the widely-used MODIS fAPAR and LAI and products. The GlobAlbedo derived products are produced using a 1D two-stream radiative transfer (RT) scheme designed explicitly for global parameter retrieval from albedo, with consistency between RT model assumptions and observations, as well as with typical large-scale land surface model RT schemes. The approach does not require biome-specific structural assumptions (e.g. cover, clumping, understory), unlike more detailed 3D RT model approaches. GlobAlbedo-derived values of fAPAR and LAI are compared with MODIS values over 2002-2011 at multiple flux tower sites within selected biomes, over 1200 × 1200 km regions and globally. GlobAlbedo-derived fAPAR and LAI values are temporally more stable than the MODIS values due to (1) the smoothness of the underlying albedo, derived via optimal estimation (assimilation) using an a priori estimate of albedo derived from an albedo ‘climatology’ (composited multi-year albedo observations) and (2) space-time invariant prior information in the inversion of the two-stream RT scheme. Parameters agree closely in timing but with GlobAlbedo values consistently lower than MODIS, particularly for LAI. Larger differences occur in winter (when values are lower) and in the Southern hemisphere. Globally, we find that: GlobAlbedo-derived fAPAR is ~0.9- 1.01 × MODIS fAPAR with an offset of ~0.03; GlobAlbedo-derived LAI is ~0.6 × MODIS LAI with an offset of ~0.2. Differences arise due to the RT model assumptions underlying the products, meaning care is required in interpreting either set of values, particularly when comparing to finescale ground-based estimates. We present global calibrations between GlobAlbedo-derived and MODIS products.JRC.H.5-Land Resources Managemen

Seasonal adaptation of the thermal‐based two‐source energy balance model for estimating evapotranspiration in a semiarid tree‐grass ecosystem

© 2020 by the authors.The thermal-based two-source energy balance (TSEB) model has accurately simulated energy fluxes in a wide range of landscapes with both remote and proximal sensing data. However, tree-grass ecosystems (TGE) have notably complex heterogeneous vegetation mixtures and dynamic phenological characteristics presenting clear challenges to earth observation and modeling methods. Particularly, the TSEB modeling structure assumes a single vegetation source, making it difficult to represent the multiple vegetation layers present in TGEs (i.e., trees and grasses) which have different phenological and structural characteristics. This study evaluates the implementation of TSEB in a TGE located in central Spain and proposes a new strategy to consider the spatial and temporal complexities observed. This was based on sensitivity analyses (SA) conducted on both primary remote sensing inputs (local SA) and model parameters (global SA). The model was subsequently modified considering phenological dynamics in semi-arid TGEs and assuming a dominant vegetation structure and cover (i.e., either grassland or broadleaved trees) for different seasons (TSEB-2S). The adaptation was compared against the default model and evaluated against eddy covariance (EC) flux measurements and lysimeters over the experimental site. TSEB-2S vastly improved over the default TSEB performance decreasing the mean bias and root-mean-square-deviation (RMSD) of latent heat (LE) from 40 and 82 W m−2 to −4 and 59 W m−2, respectively during 2015. TSEB-2S was further validated for two other EC towers and for different years (2015, 2016 and 2017) obtaining similar error statistics with RMSD of LE ranging between 57 and 63 W m−2. The results presented here demonstrate a relatively simple strategy to improve water and energy flux monitoring over a complex and vulnerable landscape, which are often poorly represented through remote sensing models.The research received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No 721995. It was also funded by Ministerio de Economía y Competitividad through FLUXPEC CGL2012-34383 and SynerTGE CGL2015-G9095-R (MINECO/FEDER, UE) projects. The research infrastructure at the measurement site in Majadas de Tiétar was partly funded through the Alexander von Humboldt Foundation, ELEMENTAL (CGL 2017-83538-C3-3-R, MINECO-FEDER) and IMAGINA (PROMETEU 2019; Generalitat Valenciana).Peer reviewe

Remote sensing phenology at European northern latitudes - From ground spectral towers to satellites

Plant phenology exerts major influences on carbon, water, and energy exchanges between atmosphere and ecosystems, provides feedbacks to climate, and affects ecosystem functioning and services. Great efforts have been spent in studying plant phenology over the past decades, but there are still large uncertainties and disputations in phenology estimation, trends, and its climate sensitivities. This thesis aims to reduce these uncertainties through analyzing ground spectral sampling, developing methods for in situ light sensor calibration, and exploring a new spectral index for reliable retrieval of remote sensing phenology and climate sensitivity estimation at European northern latitudes. The ground spectral towers use light sensors of either nadir or off-nadir viewing to measure reflected radiation, yet how plants in the sensor view contribute differently to the measured signals, and necessary in situ calibrations are often overlooked, leading to great uncertainties in ground spectral sampling of vegetation. It was found that the ground sampling points in the sensor view follow a Cauchy distribution, which is further modulated by the sensor directional response function. We proposed in situ light sensor calibration methods and showed that the user in situ calibration is more reliable than manufacturer’s lab calibration when our proposed calibration procedures are followed. By taking the full advantages of more reliable and standardized reflectance, we proposed a plant phenology vegetation index (PPI), which is derived from a radiative transfer equation and uses red and near infrared reflectance. PPI shows good linearity with canopy green leaf area index, and is correlated with gross primary productivity, better than other vegetation indices in our test. With suppressed snow influences, PPI shows great potentials for retrieving phenology over coniferous-dominated boreal forests. PPI was used to retrieve plant phenology from MODIS nadir BRDF-adjusted reflectance at European northern latitudes for the period 2000-2014. We estimated the trend of start of growing season (SOS), end of growing season (EOS), length of growing season (LOS), and the PPI integral for the time span, and found significant changes in most part of the region, with an average rate of -0.39 days·year-1 in SOS, 0.48 days·year-1 in EOS, 0.87 days·year-1 in LOS, and 0.79%·year-1 in the PPI integral over the past 15 years. We found that the plant phenology was significantly affected by climate in most part of the region, with an average sensitivity to temperature: SOS at -3.43 days·°C-1, EOS at 1.27 days·°C-1, LOS at 3.16 days·°C-1, and PPI integral at 2.29 %·°C-1, and to precipitation: SOS at 0.28 days∙cm-1, EOS at 0.05 days∙cm-1, LOS at 0.04 days∙cm-1, and PPI integral at -0.07%∙cm-1. These phenology variations were significantly related to decadal variations of atmospheric circulations, including the North Atlantic Oscillation and the Arctic Oscillation. The methods developed in this thesis can help to improve the reliability of long-term field spectral measurements and to reduce uncertainties in remote sensing phenology retrieval and climate sensitivity estimation

Spatiotemporal Patterns and Phenology of Tropical Vegetation Solar-Induced Chlorophyll Fluorescence across Brazilian Biomes Using Satellite Observations

Solar-induced fluorescence (SIF) has been empirically linked to gross primary productivity (GPP) in multiple ecosystems and is thus a promising tool to address the current uncertainties in carbon fluxes at ecosystem to continental scales. However, studies utilizing satellite-measured SIF in South America have concentrated on the Amazonian tropical forest, while SIF in other regions and vegetation classes remain uninvestigated. We examined three years of Orbiting Carbon Observatory-2 (OCO-2) SIF data for vegetation classes within and across the six Brazilian biomes (Amazon, Atlantic Forest, Caatinga, Cerrado, Pampa, and Pantanal) to answer the following: (1) how does satellite-measured SIF differ? (2) What is the relationship (strength and direction) of satellite-measured SIF with canopy temperature (T can), air temperature (T air), and vapor pressure deficit (VPD)? (3) How does the phenology of satellite-measured SIF (duration and amplitude of seasonal integrated SIF) compare? Our analysis shows that OCO-2 captures a significantly higher mean SIF with lower variability in the Amazon and lower mean SIF with higher variability in the Caatinga compared to other biomes. OCO-2 also distinguishes the mean SIF of vegetation types within biomes, showing that evergreen broadleaf (EBF) mean SIF is significantly higher than other vegetation classes (deciduous broadleaf (DBF), grassland (GRA), savannas (SAV), and woody savannas (WSAV)) in all biomes. We show that the strengths and directions of correlations of OCO-2 mean SIF to T can , T air , and VPD largely cluster by biome: negative in the Caatinga and Cerrado, positive in the Pampa, and no correlations were found in the Pantanal, while results were mixed for the Amazon and Atlantic Forest. We found mean SIF most strongly correlated with VPD in most vegetation classes in most biomes, followed by T can. Seasonality from time series analysis reveals that OCO-2 SIF measurements capture important differences in the seasonal timing of SIF for different classes, details masked when only examining mean SIF differences. We found that OCO-2 captured the highest base integrated SIF and lowest seasonal pulse integrated SIF in the Amazon for all vegetation classes, indicating continuous photosynthetic activity in the Amazon exceeds other biomes, but with small seasonal increases. Surprisingly, Pantanal EBF SIF had the highest total integrated SIF of all classes in all biomes due to a large seasonal pulse. Additionally, the length of seasons only accounts for about 30% of variability in total integrated SIF; thus, integrated SIF is likely captures differences in photosynthetic activity separate from structural differences. Our results show that satellite measurements of SIF can distinguish important functioning and phenological differences in vegetation classes and thus has the potential to improve our understanding of productivity and seasonality in the tropics

The development of a temporal-BRDF model-based approach to change detection, an application to the identification and delineation of fire affected areas.

Although large quantities of southern Africa burn every year, minimal information is available relating to the fire regimes of this area. This study develops a new, generic approach to change detection, applicable to the identification of land cover change from high temporal and moderate spatial resolution satellite data. Traditional change detection techniques have several key limitations which are identified and addressed in this work. In particular these approaches fail to account for directional effects in the remote sensing signal introduced by variations in the solar and sensing geometry, and are sensitive to underlying phenological changes in the surface as well as noise in the data due to cloud or atmospheric contamination. This research develops a bi-directional, model-based change detection algorithm. An empirical temporal component is incorporated into a semi-empirical linear BRDF model. This may be fitted to a long time series of reflectance with less sensitivity to the presence of underlying phenological change. Outliers are identified based on an estimation of noise in the data and the calculation of uncertainty in the model parameters and are removed from the sequence. A "step function kernel" is incorporated into the formulation in order to detect explicitly sudden step-like changes in the surface reflectance induced by burning. The change detection model is applied to the problem of locating and mapping fire affected areas from daily moderate spatial resolution satellite data, and an indicator of burn severity is introduced. Monthly burned area datasets for a 2400km by 1200km area of southern Africa detailing the day and severity of burning are created for a five year period (2000-2004). These data are analysed and the fire regimes of southern African ecosystems during this time are characterised. The results highlight the extent of the burning which is taking place within southern Africa, with between 27-32% of the study area burning during each of the five years of observation. Higher fire frequencies are exhibited by savanna and grassland ecosystems, while more dense vegetation types such as shrublands and deciduous broadleaf forests burn less frequently. In addition the areas which burn more frequently do so with a greater severity, with a positive relationship identified between the frequency and the severity of burning

Sun-angle effects on remote-sensing phenology observed and modelled using himawari-8

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. Satellite remote sensing of vegetation at regional to global scales is undertaken at considerable variations in solar zenith angle (SZA) across space and time, yet the extent to which these SZA variations matter for the retrieval of phenology remains largely unknown. Here we examined the effect of seasonal and spatial variations in SZA on retrieving vegetation phenology from time series of the Normalized Difference Vegetation Index (NDVI) and Enhanced Vegetation Index (EVI) across a study area in southeastern Australia encompassing forest, woodland, and grassland sites. The vegetation indices (VI) data span two years and are from the Advanced Himawari Imager (AHI), which is onboard the Japanese Himawari-8 geostationary satellite. The semi-empirical RossThick-LiSparse-Reciprocal (RTLSR) bidirectional reflectance distribution function (BRDF) model was inverted for each spectral band on a daily basis using 10-minute reflectances acquired by H-8 AHI at different sun-view geometries for each site. The inverted RTLSR model was then used to forward calculate surface reflectance at three constant SZAs (20°, 40°, 60°) and one seasonally varying SZA (local solar noon), all normalised to nadir view. Time series of NDVI and EVI adjusted to different SZAs at nadir view were then computed, from which phenological metrics such as start and end of growing season were retrieved. Results showed that NDVI sensitivity to SZA was on average nearly five times greater than EVI sensitivity. VI sensitivity to SZA also varied among sites (biome types) and phenological stages, with NDVI sensitivity being higher during the minimum greenness period than during the peak greenness period. Seasonal SZA variations altered the temporal profiles of both NDVI and EVI, with more pronounced differences in magnitude among NDVI time series normalised to different SZAs. When using VI time series that allowed SZA to vary at local solar noon, the uncertainties in estimating start, peak, end, and length of growing season introduced by local solar noon varying SZA VI time series, were 7.5, 3.7, 6.5, and 11.3 days for NDVI, and 10.4, 11.9, 6.5, and 8.4 days for EVI respectively, compared to VI time series normalised to a constant SZA. Furthermore, the stronger SZA dependency of NDVI compared with EVI, resulted in up to two times higher uncertainty in estimating annual integrated VI, a commonly used remote-sensing proxy for vegetation productivity. Since commonly used satellite products are not generally normalised to a constant sun-angle across space and time, future studies to assess the sun-angle effects on satellite applications in agriculture, ecology, environment, and carbon science are urgently needed. Measurements taken by new-generation geostationary (GEO) satellites offer an important opportunity to refine this assessment at finer temporal scales. In addition, studies are needed to evaluate the suitability of different BRDF models for normalising sun-angle across a broad spectrum of vegetation structure, phenological stages and geographic locations. Only through continuous investigations on how sun-angle variations affect spatiotemporal vegetation dynamics and what is the best strategy to deal with it, can we achieve a more quantitative remote sensing of true signals of vegetation change across the entire globe and through time