Vegetation Phenology Metrics Derived from Temporally Smoothed and Gap-filled MODIS Data

Smoothed and gap-filled VI provides a good base for estimating vegetation phenology metrics. The TIMESAT software was improved by incorporating the ancillary information from MODIS products. A simple assessment of the association between retrieved greenup dates and ground observations indicates satisfactory result from improved TIMESAT software. One application example shows that mapping Nectar Flow Phenology is tractable on a continental scale using hive weight and satellite vegetation data. The phenology data product is supporting more researches in ecology, climate change fields

ENHANCING CONSERVATION WITH HIGH RESOLUTION PRODUCTIVITY DATASETS FOR THE CONTERMINOUS UNITED STATES

Human driven alteration of the earth’s terrestrial surface is accelerating through land use changes, intensification of human activity, climate change, and other anthropogenic pressures. These changes occur at broad spatio-temporal scales, challenging our ability to effectively monitor and assess the impacts and subsequent conservation strategies. While satellite remote sensing (SRS) products enable monitoring of the earth’s terrestrial surface continuously across space and time, the practical applications for conservation and management of these products are limited. Often the processes driving ecological change occur at fine spatial resolutions and are undetectable given the resolution of available datasets. Additionally, the links between SRS data and ecologically meaningful metrics are weak. Recent advances in cloud computing technology along with the growing record of high resolution SRS data enable the development of SRS products that quantify ecologically meaningful variables at relevant scales applicable for conservation and management. The focus of my dissertation is to improve the applicability of terrestrial gross and net primary productivity (GPP/NPP) datasets for the conterminous United States (CONUS). In chapter one, I develop a framework for creating high resolution datasets of vegetation dynamics. I use the entire archive of Landsat 5, 7, and 8 surface reflectance data and a novel gap filling approach to create spatially continuous 30 m, 16-day composites of the normalized difference vegetation index (NDVI) from 1986 to 2016. In chapter two, I integrate this with other high resolution datasets and the MOD17 algorithm to create the first high resolution GPP and NPP datasets for CONUS. I demonstrate the applicability of these products for conservation and management, showing the improvements beyond currently available products. In chapter three, I utilize this dataset to evaluate the relationships between land ownership and terrestrial production across the CONUS domain. The main results of this work are three publically available datasets: 1) 30 m Landsat NDVI; 2) 250 m MODIS based GPP and NPP; and 3) 30 m Landsat based GPP and NPP. My goal is that these products prove useful for the wider scientific, conservation, and land management communities as we continue to strive for better conservation and management practices

Validation of Sentinel-2, MODIS, CGLS, SAF, GLASS and C3S leaf area index products in maize crops

Altres ajuts: this research was funded by the Copernicus Global Land Service (CGLOPS-1, 199494-JRC).We proposed a direct approach to validate hectometric and kilometric resolution leaf area index (LAI) products that involved the scaling up of field-measured LAI via the validation and recalibration of the decametric Sentinel-2 LAI product. We applied it over a test study area of maize crops in northern China using continuous field measurements of LAINet along the year 2019. Sentinel-2 LAI showed an overall accuracy of 0.67 in terms of Root Mean Square Error (RMSE) and it was used, after recalibration, as a benchmark to validate six coarse resolution LAI products: MODIS, Copernicus Global Land Service 1 km Version 2 (called GEOV2) and 300 m (GEOV3), Satellite Application Facility EUMETSAT Polar System (SAF EPS) 1.1 km, Global LAnd Surface Satellite (GLASS) 500 m and Copernicus Climate Change Service (C3S) 1 km V2. GEOV2, GEOV3 and MODIS showed a good agreement with reference LAI in terms of magnitude (RMSE ≤ 0.29) and phenology. SAF EPS (RMSE = 0.68) and C3S V2 (RMSE = 0.41), on the opposite, systematically underestimated high LAI values and showed systematic differences for phenological metrics: a delay of 6 days (d), 20 d and 24 d for the start, peak and the end of growing season, respectively, for SAF EPS and an advance of −4 d, −6 d and −6 d for C3S

Land surface phenology from VEGETATION and PROBA-V data. Assessment over deciduous forests

Land surface phenology has been widely retrieved although no consensus exists on the optimal satellite dataset and the method to extract phenology metrics. This study is the first comprehensive comparison of vegetation variables and methods to retrieve land surface phenology for 1999-2017 time series of Copernicus Global Land products derived from SPOT-VEGETATION and PROBA-V data. We investigated the sensitivity of phenology to (I) the input vegetation variable: normalized difference vegetation index (NDVI), leaf area index (LAI), fraction of absorbed photosynthetically active radiation (FAPAR), and fraction of vegetation cover (FCOVER); (II) the smoothing and gap filling method for deriving seasonal trajectories; and (III) the method to extract phenological metrics: thresholds based on a percentile of the annual amplitude of the vegetation variable, autoregressive moving averages, logistic function fitting, and first derivative methods. We validated the derived satellite phenological metrics (start of the season (SoS) and end of the season (EoS)) using available ground observations of Betula pendula, B. alleghaniensis, Acer rubrum, Fagus grandifolia, and Quercus rubra in Europe (Pan-European PEP725 network) and the USA (National Phenology Network, USA-NPN). The threshold-based method applied to the smoothed and gap-filled LAI V2 time series agreed best with the ground phenology, with root mean square errors of ˜10 d and ˜25 d for the timing of SoS and EoS respectively. This research is expected to contribute for the operational retrieval of land surface phenology within the Copernicus Global Land Servic

SATELLITE MICROWAVE MEASUREMENT OF LAND SURFACE PHENOLOGY: CLARIFYING VEGETATION PHENOLOGY RESPONSE TO CLIMATIC DRIVERS AND EXTREME EVENTS

The seasonality of terrestrial vegetation controls feedbacks to the climate system including land-atmosphere water, energy and carbon (CO2) exchanges with cascading effects on regional-to-global weather and circulation patterns. Proper characterization of vegetation phenology is necessary to understand and quantify changes in the earthÆs ecosystems and biogeochemical cycles and is a key component in tracking ecological species response to climate change. The response of both functional and structural vegetation phenology to climatic drivers on a global scale is still poorly understood however, which has hindered the development of robust vegetation phenology models. In this dissertation I use satellite microwave vegetation optical depth (VOD) in conjunction with an array of satellite measures, Global Positioning System (GPS) reflectometry, field observations and flux tower data to 1) clarify vegetation phenology response to water, temperature and solar irradiance constraints, 2) demonstrate the asynchrony between changes in vegetation water content and biomass and changes in greenness and leaf area in relation to land cover type and climate constraints, 3) provide enhanced assessment of seasonal recovery of vegetation biomass following wildfire and 4) present a method to more accurately model tropical vegetation phenology. This research will establish VOD as a useful and informative parameter for regional-to-global vegetation phenology modeling, more accurately define the drivers of both structural and functional vegetation phenology, and help minimize errors in phenology simulations within earth system models. This dissertation also includes the development of Gross Primary Productivity (GPP) and Net Primary Productivity (NPP) vegetation health climate indicators as part of a NASA funded project entitled Development and Testing of Potential Indicators for the National Climate Assessment; Translating EOS datasets into National Ecosystem Biophysical Indicators

Field-scale mapping of evaporative stress indicators of crop yield: An application over Mead, NE, USA

The Evaporative Stress Index (ESI) quantifies temporal anomalies in a normalized evapotranspiration (ET) metric describing the ratio of actual-to-reference ET (fRET) as derived from satellite remote sensing. At regional scales (3–10 km pixel resolution), the ESI has demonstrated the capacity to capture developing crop stress and impacts on regional yield variability in water-limited agricultural regions. However, its performance in some regions where the vegetation cycle is intensively managed appears to be degraded due to spatial and temporal limitations in the standard ESI products. In this study, we investigated potential improvements to ESI by generating maps of ET, fRET, and fRET anomalies at high spatiotemporal resolution (30-m pixels, daily time steps) using a multi-sensor data fusion method, enabling separation of landcover types with different phenologies and resilience to drought. The study was conducted for the period 2010–2014 covering a region around Mead, Nebraska that includes both rainfed and irrigated crops. Correlations between ESI and measurements of maize yield were investigated at both the field and county level to assess the potential of ESI as a yield forecasting tool. To examine the role of crop phenology in yield-ESI correlations, annual input fRET time series were aligned by both calendar day and by biophysically relevant dates (e.g. days since planting or emergence). At the resolution of the operational U.S. ESI product (4 km), adjusting fRET alignment to a regionally reported emergence date prior to anomaly computation improves r2 correlations with county-level yield estimates from 0.28 to 0.80. At 30-m resolution, where pure maize pixels can be isolated from other crops and landcover types, county-level yield correlations improved from 0.47 to 0.93 when aligning fRET by emergence date rather than calendar date. Peak correlations occurred 68 days after emergence, corresponding to the silking stage for maize when grain development is particularly sensitive to soil moisture deficiencies. The results of this study demonstrate the utility of remotely sensed ET in conveying spatially and temporally explicit water stress information to yield prediction and crop simulation models

Responses of Land Surface Phenology to Wildfire Disturbances in the Western United States Forests

Land surface phenology (LSP) characterizes the seasonal dynamics in the vegetation communities observed for a satellite pixel and it has been widely associated with global climate change. However, LSP and its long-term trend can be influenced by land disturbance events, which could greatly interrupt the LSP responses to climate change. Wildfire is one of the main disturbance agents in the western United States (US) forests, but its impacts on LSP have not been investigated yet. To gain a comprehensive understanding of the LSP responses to wildfires in the western US forests, this dissertation focused on three research objectives: (1) to perform a case study of wildfire impacts on LSP and its trend by comparing the burned and a reference area, (2) to investigate the distribution of wildfire impacts on LSP and identify control factors by analyzing all the wildfires across the western US forests, and (3) to quantify the contributions of land cover composition and other environmental factors to the spatial and interannual variations of LSP in a recently burned landscape. The results reveal that wildfires play a significant role in influencing spatial and interannual variations in LSP across the western US forests. First, the case study showed that the Hayman Fire significantly advanced the start of growing season (SOS) and caused an advancing SOS trend comparing with a delaying trend in the reference area. Second, summarizing \u3e800 wildfires found that the shifts in LSP timing were divergent depending on individual wildfire events and burn severity. Moreover, wildfires showed a stronger impact on the end of growing season (EOS) than SOS. Last, LSP trends were interrupted by wildfires with the degree of impact largely dependent on the wildfire occurrence year. Third, LSP modeling showed that land cover composition, climate, and topography co-determine the LSP variations. Specifically, land cover composition and climate dominate the LSP spatial and interannual variations, respectively. Overall, this research improves the understanding of wildfire impacts on LSP and the underlying mechanism of various factors driving LSP. This research also provides a prototype that can be extended to investigate the impacts on LSP from other disturbances

Temporal Requirements for Future Landsat Systems for Agricultural Monitoring

Agricultural monitoring is an important application of earth-observing satellite systems, which may be used for stress and disease detection, growth stage monitoring, and yield prediction in crops at a fraction of the time and cost it would take to survey fields manually. Satellites within the Landsat program are frequently used for agricultural monitoring, but they do not always collect imagery often enough to capture rapid changes in vegetation health. To address this limitation, an increase in revisit rate is being considered for future Landsat systems. This research aims to determine the necessary overpass frequency for a future Landsat sensor for agricultural growth stage monitoring and yield prediction. Two experiments were conducted to study the effects of temporal resolution on the accuracy of these tasks. The first experiment investigated the impact of imaging frequency on growth stage monitoring. Image-derived plot-average Normalized Difference Vegetation Index (NDVI) time-series data collected over a small corn field were used to estimate phenological transition dates. Images were then removed from the original time-series, and dates were recalculated from the resampled data. Using PlanetScope surface reflectance imagery, the average range of estimated dates increased by a day and the average absolute deviation between estimated dates increased by 1/3 of a day for every day of increase in average revisit interval. Using the higher-quality PlanetScope L3H surface reflectance product, these rates of increase were approximately halved. Higher imaging frequency and higher radiometric quality both led to greater precision in estimates. The second experiment investigated the impacts of imaging frequency and time-series end date on yield correlation accuracy. Plot-average Green Normalized Difference Vegetation Index (GNDVI) time-series data collected over a small corn field during two different growing seasons were resampled to different revisit intervals, gap-filled and smoothed using two different methods, and correlated with plot-average yield at each day of the growing season. These experiments were then repeated with images removed from the end of the time-series. All methods tested performed well on time-series ending 65-72 days or more after green-up, and performed poorly for time-series ending prior to the day of peak GNDVI. Mean R-squared values for GNDVI-yield correlations decreased with increasing revisit intervals. This effect was stronger for the more typical 2019 data, as well as for time-series ending earlier in the growing season. The findings of this study, along with cloud contamination statistics, were used to recommend an overpass frequency of 1-4 days for future yield-monitoring satellite systems. The optimal frequency within this range depends on the specific task being attempted

Fernerkundung der Vegetationsphänologie über MODIS NDVI Daten - Herausforderungen bei der Datenverarbeitung und -validierung mittels Bodenbeobachtungen zahlreicher Arten und LiDAR

Phenology, the cyclic events in living organisms is triggered by climatic conditions and indicators of climate change. They are important factors influencing species interactions and ecosystem functioning. This thesis deals with the estimation of phenological metrics (Land Surface Phenology or LSP) from MODIS based time series NDVI data. Results of data analysis emphasises the role of ground observations, topography and LiDAR characteristics of forest stand in describing the variability in LSP.Phänologie, die zyklischen Stadien von lebenden Organismen werden über klimatische Verhältnisse gesteuert und dienen als Indikatoren des Klimawandels. Diese Faktoren beeinflussen maßgeblich die Interaktionen zwischen Arten und sind für das Funktionieren von Ökosystemen ausschlaggebend. Diese Arbeit behandelt die Bestimmung von phänologischen Metriken (Phänologie der Landoberfläche oder LSP) unter Verwendung von MODIS basierten NDVI Zeitreihen. Die Ergebnisse der Datenanalyse hebt die Wichtigkeit von Bodenbeobachtungen, Topographie und LiDAR Merkmalen von Waldbeständen bei der Beschreibung der LSP Variabilität hervor

Global Analysis of Bioclimatic Controls on Ecosystem Productivity Using Satellite Observations of Solar-Induced Chlorophyll Fluorescence

Ecosystem productivity models rely on regional climatic information to estimate temperature and moisture constraints influencing plant growth. However, the productivity response to these environmental factors is uncertain at the global scale and has largely been defined using limited observations from sparse monitoring sites, including carbon flux towers. Recent studies have shown that satellite observations of Solar-Induced chlorophyll Fluorescence (SIF) are highly correlated with ecosystem Gross Primary Productivity (GPP). Here, we use a relatively long-term global SIF observational record from the Global Ozone Monitoring Experiment-2 (GOME-2) sensors to investigate the relationships between SIF, used as a proxy for GPP, and selected bio-climatic factors constraining plant growth at the global scale. We compared the satellite SIF retrievals with collocated GPP observations from a global network of tower carbon flux monitoring sites and surface meteorological data from model reanalysis, including soil moisture, Vapor Pressure Deficit (VPD), and minimum daily air temperature (Tmin). We found strong correspondence (R2 \u3e 80%) between SIF and GPP monthly climatologies for tower sites characterized by mixed, deciduous broadleaf, evergreen needleleaf forests, and croplands. For other land cover types including savanna, shrubland, and evergreen broadleaf forest, SIF showed significant but higher variability in correlations between sites. In order to analyze temperature and moisture related effects on ecosystem productivity, we divided SIF by photosynthetically active radiation (SIFp) and examined partial correlations between SIFp and the climatic factors across a global range of flux tower sites, and over broader regional and global extents. We found that productivity in arid ecosystems is more strongly controlled by soil water content to an extent that soil moisture explains a higher proportion of the seasonal cycle in productivity than VPD. At the global scale, ecosystem productivity is affected by joint climatic constraint factors so that VPD, Tmin, and soil moisture were significant (p \u3c 0.05) controls over 60, 59, and 35 percent of the global domain, respectively. Our study identifies and confirms dominant climate control factors influencing productivity at the global scale indicated from satellite SIF observations. The results are generally consistent with climate response characteristics indicated from sparse global tower observations, while providing more extensive coverage for verifying and refining global carbon and climate model assumptions and predictions