Land surface phenology: What do we really 'see' from space?

Modelling asynchrony in phenology considering a dynamic representation of meteorological variables

Shifts in the timing of phenological events across many taxa and ecosystems are a result of climate change. Within a trophic network, phenological mismatches between interlinked species can have negative impacts for biodiversity, ecosystems, and the trophic network. Here we developed interaction indices that quantify the level of synchrony and asynchrony among groups of species in three interlinked trophic levels, as well as accounting for a dynamic representation of meteorology. Insect first flight, vegetation green-up and arrival of migrant birds were the phenological indicators, obtained from a combination of spatially and temporally explicit species observations from citizen science programmes and remote sensing platforms (i.e., Landsat). To determine phenological shifts in interlinked taxa we created and applied several phenological indices of synchrony-asynchrony, combining information from the phenological events and critical time windows of meteorological variables. To demonstrate our method of incorporating a meteorological component in our new interaction index, we implemented the relative sliding time window analysis, a stepwise regression model, to identify critical time windows preceding the phenological events on a yearly basis. The new indices of phenological change identified several asynchronies within trophic levels, allowing exploration of potential interactions based on synchrony among interlinked species. Our novel index of synchrony-asynchrony including a meteorological dimension could be highly informative and should open new pathways for studying synchrony among species and interaction networks.

Comparative analysis of global urban land surface phenology between the MODIS and VIIRS products and extraction methods

The reliability of land surface phenology (LSP) derived from satellite remote sensing is crucial for obtaining accurate estimates of the phenological response of vegetation to future climate change in urban ecosystems. Differences in phenological definition and extraction methodology using remote sensing can generate systemic errors in estimating the phenological temperature sensitivity to predict the biological response of vegetation. Here, we evaluated the start of the season (SOS), the end of the season (EOS), and the growing season length (GSL) between the Terra and Aqua combined Moderate Resolution Imaging Spectroradiometer (MODIS) Land Cover Dynamics (MCD12Q2) and the Suomi National Polar-Orbiting Partnership NASA Visible Infrared Imaging Radiometer Suite (VIIRS) Land Cover Dynamics (VNP22Q2) over 1470 urban clusters worldwide. We found that the availability of valid phenological estimation in urban areas was low across climates, and cities without valid phenological data were mainly distributed in tropical and dry climate zones. The phenological comparison for each urban cluster worldwide revealed significant differences between the two products. For SOS, EOS, and GSL, only a small proportion of urban clusters (6.01%, 3.54%, and 1.50%, respectively) exhibited R2 surpassing 0.5, and approximately 31.63%, 37.37%, and 58.57% showed the RMSD values exceeding two weeks. These results remained consistent and robust across various climates, latitudinal gradients, and valid pixel proportions. We also revealed discrepancies in the inter-annual variability of LSP between the two products. These disparities in phenological metrics directly lead to systematic uncertainties in the measurements of temperature sensitivity. Furthermore, urban phenology disparities primarily arise from the different phenology definitions employed by the two products. Standardizing the definition of phenology metrics extraction from remote sensing data is crucial to the comprehensive understanding of phenological responses to urban environments, which thus provides referencing resources for studying how future environmental changes affect vegetation in natural ecosystems.

Analysing long-term spatiotemporal land surface phenology patterns over the Iberian Peninsula using 250 m MODIS EVI2 data

Iberian Peninsula ecosystems are especially vulnerable to global warming, and variations in climate patterns may alter the wide variety of services they provide. Despite this, seasonal variations in Iberian ecosystems have been understudied. Thus, this study aims to characterise land surface phenology (LSP) patterns in the Iberian Peninsula over 21 years (2001-2021), considering three phenometrics: the start (SOS), the end (EOS), and the length of the growing season (LOS). These were estimated from 8-day image composites of EVI2 (Two-band Enhanced Vegetation Index), derived from the surface reflectance product MOD09Q1 at a 250-metre spatial resolution. Phenometrics and in-situ human phenological observations of plant phenophases were also compared. Pearson's correlation coefficient, p-value, and absolute differences between paired phenometrics and phenophases were calculated to quantify uncertainty between both phenological approaches. Generally, SOS and EOS dates were later in the Alpine and Atlantic biogeographic regions. SOS occurred in March-April and EOS between October-December. Natural vegetation land cover types had similar phenological dynamics, with SOS occurring between late winter and spring and EOS between autumn and early winter. However, phenometric dates were earlier in Mediterranean savannas, grasslands, and scrublands. Atlantic evergreen broadleaf vegetation showed the strongest correlations between SOS and the first (r = 0.96) and second (r = 0.68) leaf unfolding. In contrast, Mediterranean evergreen broadleaf vegetation had unclear correlations. Atlantic deciduous broadleaf vegetation showed a moderate correlation between SOS and first (r = 0.52) and second (r = 0.57) leaf unfolding. The correlation between SOS and the first (r = 0.14) and second (r = 0.13) leaf unfolding was weaker for Mediterranean deciduous broadleaf vegetation. EOS and autumn phenophases generally showed unclear consistency. Higher spatial resolution satellite data may improve the consistency between EOS and autumn phenophases in Iberian ecosystems, as well as between SOS and spring phenophases in heterogeneous Mediterranean ecosystems.

LiPheStream - A 18-month high spatiotemporal resolution point cloud time series of Boreal trees from Finland

In the present paper, we introduce a high-resolution spatiotemporal point cloud time series, acquired using a LiDAR sensor mounted 30 metres above ground on a flux observation tower monitoring a boreal forest. The dataset comprises a 18-month long (April 2020 - September 2021) time series with an average interval of 3.5 days between observations. The data acquisition, transfer, and storage systems established at Hyytiälä (Finland) are named the LiDAR Phenology station (LiPhe). The dataset consists of 103 time points of LiDAR point clouds covering a total of 458 individual trees, comprising three distinct Boreal species. Additional reference information includes the respective location, the species, and the initial height (at the first time point) of each individual tree. The processing scripts are included to outline the workflow used to generate the individual tree point clouds (LiPheKit). The presented dataset offers a comprehensive insight into inter- and intra-species variations of the individual trees regarding their growth strategies, phenological dynamics, and other functioning processes over two growth seasons.

Maize yield prediction and condition monitoring at the sub-county scale in Kenya: synthesis of remote sensing information and crop modeling

Agricultural production assessments are crucial for formulating strategies for closing yield gaps and enhancing production efficiencies. While in situ crop yield measurements can provide valuable and accurate information, such approaches are costly and lack scalability for large-scale assessments. Therefore, crop modeling and remote sensing (RS) technologies are essential for assessing crop conditions and predicting yields at larger scales. In this study, we combined RS and a crop growth model to assess phenology, evapotranspiration (ET), and yield dynamics at grid and sub-county scales in Kenya. We synthesized RS information from the Food and Agriculture Organization (FAO) Water Productivity Open-access portal (WaPOR) to retrieve sowing date information for driving the model simulations. The findings showed that grid-scale management information and progressive crop growth could be accurately derived, reducing the model output uncertainties. Performance assessment of the modeled phenology yielded satisfactory accuracies at the sub-county scale during two representative seasons. The agreement between the simulated ET and yield was improved with the combined RS-crop model approach relative to the crop model only, demonstrating the value of additional large-scale RS information. The proposed approach supports crop yield estimation in data-scarce environments and provides valuable insights for agricultural resource management enabling countermeasures, especially when shortages are perceived in advance, thus enhancing agricultural production.

Varying effects of tree cover on relationships between satellite-observed vegetation greenup date and spring temperature across Eurasian boreal forests

Earth observation satellites have facilitated the quantification of how vegetation phenology responds to climate warming on large scales. However, satellite image pixels may contain a mixture of multiple vegetation types or species with diverse phenological responses to climate variability. It is unclear how these mixed pixels affect the statistical relationships between satellite-derived vegetation phenology and climate factors. Here, we aim to investigate the impacts of percent tree cover (PTC), a measure of mixed pixel, on the statistical relationships between satellite-derived vegetation greenup date (GUD) and spring air temperature across Eurasian boreal forests at a 0.05° spatial resolution. We estimated GUD using Moderate Resolution Imaging Spectroradiometer (MODIS) time series data. The responses of GUD to interannual variation in spring temperature (April to May) during 2001-2020 were characterized by correlation coefficient (RTAM) and sensitivity (STAM). We then evaluated the local impacts of PTC on spatial variations in RTAM and STAM using partial correlation analysis through spatial moving windows. Our results indicate that, for most areas, forests with higher PTC were associated with stronger RTAM and STAM. Moreover, PTC had stronger local impacts on RTAM and STAM than mean annual temperature and temperature seasonality for 37.3% and 27.4% of the moving windows, respectively. These impacts were spatially varying and different among forest types. Specifically, deciduous broadleaf forests and deciduous needleleaf forests tend to have a higher proportion of these impacts compared to other forest types. Our findings demonstrate the nonnegligible effects of PTC on the statistical responses of GUD to temperature variability at coarse spatial resolution (0.05°) across Eurasian boreal forests.

"Residential greenness, gestational diabetes mellitus (GDM) and microbiome diversity during pregnancy"

BACKGROUND

Gestational diabetes mellitus (GDM) is associated with reduced gut microbiota richness that was also reported to differ significantly between those living in rural compared to urban environments. Therefore, our aim was to examine the associations between greenness and maternal blood glucose levels and GDM, with microbiome diversity as a possible mediator in these associations.

METHODS

Pregnant women were recruited between January 2016 and October 2017. Residential greenness was evaluated as mean Normalized Difference Vegetation Index (NDVI) within 100, 300 and 500 m buffers surrounding each maternal residential address. Maternal glucose levels were measured at 24-28 weeks of gestation and GDM was diagnosed. We estimated the associations between greenness and glucose levels and GDM using generalized linear models, adjusting for socioeconomic status and season at last menstrual period. Using causal mediation analysis, the mediation effects of four different indices of microbiome alpha diversity in first trimester stool and saliva samples were assessed.

RESULTS

Of 269 pregnant women, 27 participants (10.04%) were diagnosed with GDM. Although not statistically significant, adjusted exposure to medium tertile levels of mean NDVI at 300 m buffer had lower odds of GDM (OR = 0.45, 95% CI: 0.16, 1.26, p = 0.13) and decreased change in mean glucose levels (β = -6.28, 95% CI: 14.91, 2.24, p = 0.15) compared to the lowest tertile levels of mean NDVI. Mixed results were observed at 100 and 500 m buffers, and when comparing highest tertile levels to lowest. No mediation effect of first trimester microbiome on the association between residential greenness and GDM was observed, and a small, possibly incidental, mediation effect on glucose levels was observed.

CONCLUSION

Our study suggests possible associations between residential greenness and glucose intolerance and risk of GDM, though without sufficient evidence. Microbiome in the first trimester, while involved in GDM etiology, is not a mediator in these associations. Future studies in larger populations should further examine these associations.

Increasing temperature regulates the advance of peak photosynthesis timing in the boreal ecosystem

The shift in vegetation phenology is an essential indicator of global climate change. Numerous researches based on reflectance-based vegetation index data have explored the changes in the start (SOS) and end (EOS) of vegetation life events at long time scales, while a huge discrepancy existed between the phenological metrics of vegetation structure and function. The peak photosynthesis timing (PPT), which is crucial in regulating terrestrial ecosystem carbon balance, has not received much attention. Using two global reconstructed solar-induced chlorophyll fluorescence data (CSIF and GOSIF) directly associated with vegetation photosynthesis, the spatio-temporal dynamics in PPT as well as the key environmental controls across the boreal ecosystem during 2001-2019 were systematically explored. Multi-year mean pattern showed that PPT mainly appeared in the first half of July. Compared to the northern Eurasia, later PPT appeared in the northern North America continent for about 4-5 days. Meanwhile, spatial trend in PPT exhibited an advanced trend during the last two decades. Especially, shrubland and grassland were obvious among all biomes. Spatial partial correlation analysis revealed that preseason temperature was the dominant environmental driver of PPT trends, occupying 81.32% and 78.04% of the total pixels of PPT_CSIF and PPT_GOSIF, respectively. Attribution analysis by ridge regression again emphasized the largest contribution of temperature to PPT dynamics in the boreal ecosystem by 52.22% (PPT_CSIF) and 46.59% (PPT_GOSIF), followed by radiation (PPT_CSIF: 24.44%; PPT_GOSIF: 28.66%) and precipitation (PPT_CSIF: 23.34%; PPT_GOSIF: 24.75%). These results have significant implications for deepening our understanding between vegetation photosynthetic phenology and carbon cycling with respect to future climate change in the boreal ecosystem.

"Residential greenness and site-specific cancer: A registry based cohort of 144,427 participants with a 21-years of follow-up, Tel-Aviv district, Israel"

BACKGROUND

Few longitudinal studies evaluated the beneficial associations between cumulative residential greenness and site-specific cancer. Our objective was to evaluate the associations between cumulative residential greenness exposure and site-specific cancer incidence (lung, bladder, breast, prostate, and skin cancer) within a registry-based cohort study.

METHODS

This study was based on 144,427 participants who lived in the Tel Aviv district during 1995-2015. The residential greenness exposure was estimated for every participant, as the weighted mean residential greenness exposure, based on the mean Normalized Difference Vegetation Index (NDVI) in the residential area and the duration of the residence in this area. Cox regression models were used to evaluate the unadjusted and adjusted associations between exposure to greenness and cancer incidence during 1998-2015 (Hazard Ratios (HRs) and 95% Confidence Intervals (CIs)). Covariates included in adjusted models were selected based on prior knowledge and directed acyclic graphs. We imputed missing data and further sensitivity analyses were conducted.

RESULTS

After adjustments, beneficial associations between exposure to greenness and cancer incidence were observed. An interquartile range (IQR) increase in NDVI was associated with a lower HRs for lung cancer (HR_adj. = 0.75 95% CI: 0.66-0.85), bladder cancer (HR_adj. = 0.71, 95% CI: 0.62-0.82), breast cancer (HR_adj. = 0.81, 95% CI: 0.74-0.88), prostate cancer (HR_adj. = 0.77, 95% CI: 0.70-0.86) and skin cancer (HR_adj. = 0.78, 95% CI: 0.69-0.88). Generally, the patterns of associations were consistent between complete-case models and imputed models, when estimated for participants aged 16 years or 40 years and older at baseline, when stratified by area level socioeconomic status, when evaluated for non-movers participants and after further adjustment to social determinants of health.

CONCLUSION

Residential greenness may reduce the risk for lung, bladder, breast, prostate, and skin cancers. If our observations will be replicated, it may present a useful avenue for public-health intervention to reduce cancer burden.

LAI-Based Phenological Changes and Climate Sensitivity Analysis in the Three-River Headwaters Region

Global climate changes have a great impact on terrestrial ecosystems. Vegetation is an important component of ecosystems, and the impact of climate changes on ecosystems can be determined by studying vegetation phenology. Vegetation phenology refers to the phenomenon of periodic changes in plants, such as germination, flowering and defoliation, with the seasonal change of climate during the annual growth cycle, and it is considered to be one of the most efficient indicators to monitor climate changes. This study collected the global land surface satellite leaf area index (GLASS LAI) products, meteorological data sets and other auxiliary data in the Three-River headwaters region from 2001 to 2018; rebuilt the vegetation LAI annual growth curve by using the asymmetric Gaussian (A-G) fitting method and extracted the three vegetation phenological data (including Start of Growing Season (SOS), End of Growing Season (EOS) and Length of Growing Season (LOS)) by the maximum slope method. In addition, it also integrated Sen’s trend analysis method and the Mann-Kendall test method to explore the temporal and spatial variation trends of vegetation phenology and explored the relationship between vegetation phenology and meteorological factors through a partial correlation analysis and multiple linear regression models. The results of this study showed that: (1) the SOS of vegetation in the Three-River headwaters region is concentrated between the beginning and the end of May, with an interannual change rate of −0.14 d/a. The EOS of vegetation is concentrated between the beginning and the middle of October, with an interannual change rate of 0.02 d/a. The LOS of vegetation is concentrated between 4 and 5 months, with an interannual change rate of 0.21 d/a. (2) Through the comparison and verification with the vegetation phenological data observed at the stations, it was found that the precision of the vegetation phonology extracted by the A-G method and the maximum slope method based on GLASS LAI data is higher (MAE is 7.6 d, RMSE is 8.4 d) and slightly better than the vegetation phenological data (MAE is 9.9 d, RMSE is 10.9 d) extracted based on the moderate resolution imaging spectroradiometer normalized difference vegetation index (MODIS NDVI) product. (3) The correlation between the SOS of vegetation and the average temperature in March–May is the strongest. The SOS of vegetation is advanced by 1.97 days for every 1 °C increase in the average temperature in March–May; the correlation between the EOS of vegetation and the cumulative sunshine duration in August–October is the strongest. The EOS of vegetation is advanced by 0.07 days for every 10-h increase in the cumulative sunshine duration in August–October.

Summer temperature—but not growing season length—influences radial growth of Salix arctica in coastal Arctic tundra

Arctic climate change is leading to an advance of plant phenology (the timing of life history events) with uncertain impacts on tundra ecosystems. Although the lengthening of the growing season is thought to lead to increased plant growth, we have few studies of how plant phenology change is altering tundra plant productivity. Here, we test the correspondence between 14 years of Salix arctica phenology data and radial growth on Qikiqtaruk–Herschel Island, Yukon Territory, Canada. We analysed stems from 28 individuals using dendroecology and linear mixed-effect models to test the statistical power of growing season length and climate variables to individually predict radial growth. We found that summer temperature best explained annual variation in radial growth. We found no strong evidence that leaf emergence date, earlier leaf senescence date, or total growing season length had any direct or lagged effects on radial growth. Radial growth was also not explained by interannual variation in precipitation, MODIS surface greenness (NDVI), or sea ice concentration. Our results demonstrate that at this site, for the widely distributed species S. arctica, temperature—but not growing season length—influences radial growth. These findings challenge the assumption that advancing phenology and longer growing seasons will increase the productivity of all plant species in Arctic tundra ecosystems.

High-Throughput Remote Sensing of Vertical Green Living Walls (VGWs) in Workplaces

Vertical green living walls (VGWs)—growing plants on vertical walls inside or outside buildings—have been suggested as a nature-based solution to improve air quality and comfort in modern cities. However, as with other greenery systems (e.g., agriculture), managing VGW systems requires adequate temporal and spatial monitoring of the plants as well as the surrounding environment. Remote sensing cameras and small, low-cost sensors have become increasingly valuable for conventional vegetation monitoring; nevertheless, they have rarely been used in VGWs. In this descriptive paper, we present a first-of-its-kind remote sensing high-throughput monitoring system in a VGW workplace. The system includes low- and high-cost sensors, thermal and hyperspectral remote sensing cameras, and in situ gas-exchange measurements. In addition, air temperature, relative humidity, and carbon dioxide concentrations are constantly monitored in the operating workplace room (scientific computer lab) where the VGW is established, while data are continuously streamed online to an analytical and visualization web application. Artificial Intelligence is used to automatically monitor changes across the living wall. Preliminary results of our unique monitoring system are presented under actual working room conditions while discussing future directions and potential applications of such a high-throughput remote sensing VGW system.

Evaluation of Vegetation Indexes and Green-Up Date Extraction Methods on the Tibetan Plateau

The vegetation green-up date (GUD) of the Tibetan Plateau (TP) is highly sensitive to climate change. Accurate estimation of GUD is essential for understanding the dynamics and stability of terrestrial ecosystems and their interactions with climate. The GUD is usually determined from a time-series of vegetation indices (VIs). The adoption of different VIs and GUD extraction methods can lead to different GUDs. However, our knowledge of the uncertainty in these GUDs on TP is still limited. In this study, we evaluated the performance of different VIs and GUD extraction methods on TP from 2003 to 2020. The GUDs were determined from six Moderate Resolution Imaging Spectroradiometer (MODIS) derived VIs: normalized difference vegetation index (NDVI), enhanced vegetation index (EVI), normalized difference infrared index (NDII), phenology index (PI), normalized difference phenology index (NDPI), and normalized difference greenness index (NDGI). Four extraction methods (βmax, CCRmax, G20, and RCmax) were applied individually to each VI to determine GUD. The GUDs obtained from all VIs showed similar patterns of early green-up in the eastern and late green-up in the western plateau, and similar trend of GUD advancement in the eastern and postponement in the western plateau. The accuracy of the derived GUDs was evaluated by comparison with ground-observed GUDs from 19 agrometeorological stations. Our results show that two snow-free VIs, NDGI and NDPI, had better performance in GUD extraction than the snow-calibrated conventional VIs, NDVI and EVI. Among all the VIs, NDGI gave the highest GUD accuracy when combined with the four extraction methods. Based on NDGI, the GUD extracted by the CCRmax method was found to have the highest consistency (r = 0.62, p < 0.01, RMSE = 11 days, bias = −3.84 days) with ground observations. The NDGI also showed the highest accuracy for preseason snow-covered site-years (r = 0.71, p < 0.01, RMSE = 10.69 days, bias = −4.05 days), indicating its optimal resistance to snow cover influence. In comparison, NDII and PI hardly captured GUD. NDII was seriously affected by preseason snow cover, as indicated by the negative correlation coefficient (r = −0.34, p < 0.1), high RMSE and bias (RMSE = 50.23 days, bias = −24.25 days).

Characterising the Land Surface Phenology of Middle Eastern Countries Using Moderate Resolution Landsat Data

Global change impacts including climate change, increased CO2 and nitrogen deposition can be determined through a more precise characterisation of Land Surface Phenology (LSP) parameters. In addition, accurate estimation of LSP dates is being increasingly used in applications such as mapping vegetation types, yield forecasting, and irrigation management. However, there has not been any attempt to characterise Middle East vegetation phenology at the fine spatial resolution appropriate for such applications. Remote-sensing based approaches have proved to be a useful tool in such regions since access is restricted in some areas due to security issues and their inter-annual vegetation phenology parameters vary considerably because of high uncertainty in rainfall. This study aims to establish for the first time a comprehensive characterisation of the vegetation phenological characteristics of the major vegetation types in the Middle East at a fine spatial resolution of 30 m using Landsat Normalized Difference Vegetation Index (NDVI) time series data over a temporal range of 20 years (2000–2020). Overall, a progressive pattern in phenophases was observed from low to high latitude. The earliest start of the season was concentrated in the central and east of the region associated mainly with grassland and cultivated land, while the significantly delayed end of the season was mainly distributed in northern Turkey and Iran corresponding to the forest, resulting in the prolonged length of the season in the study area. There was a significant positive correlation between LSP parameters and latitude, which indicates a delay in the start of the season of 4.83 days (R2 = 0.86, p < 0.001) and a delay in the end of the season of 6.54 days (R2 = 0.83, p < 0.001) per degree of latitude increase. In addition, we have discussed the advantages of fine resolution LSP parameters over the available coarse datasets and showed how such outputs can improve many applications in the region. This study shows the potential of Landsat data to quantify the LSP of major land cover types in heterogeneous landscapes of the Middle East which enhances our understanding of the spatial-temporal dynamics of vegetation dynamics in arid and semi-arid settings in the world.

A global increase in tree cover extends the growing season length as observed from satellite records

Plant phenology provides information on the seasonal dynamics of plants, and changes herein are important for understanding the impact of climate change and human management on the biosphere. Land surface phenology is the study of plant phenology across large spatial scales estimated by satellite observations. However, satellite observations (pixels) are often composed of a mixture of vegetation types, like woody vegetation and herbaceous vegetation, having different phenological characteristics. Therefore, any changes in tree cover presumably impact land surface phenology, as trees usually have a different seasonal cycle compared to herbaceous vegetation. On the other hand, changes in land surface phenology are often interpreted as a result of climate change-induced impacts on the photosynthetic activity of vegetation. Therefore, it is important to better understand the role of changes in vegetation cover (here, the proportion between tree and short vegetation cover) in satellite-derived land surface phenology analysis. We studied the impact of changes in tree cover on satellite observed land surface phenology at a global scale over the past three decades. We found an extension of the growing season length in 36.6% of the areas where tree cover increased, whereas only 20.1% of the areas where tree cover decreased showed an increase in growing season length. Furthermore, the ratio between tree cover and short vegetation cover was found to affect changes in the length of the growing season, with the denser tree cover showing a more pronounced extension of the growing season length (especially in boreal forests). These results highlight the importance of changes in tree cover when analyzing the impact of climate change on vegetation phenology. Our study thereby addresses a critical knowledge gap for an improved understanding of changes in land surface phenology during recent decades in the context of climate and human-induced global land cover change.

Macrophenology: insights into the broad-scale patterns, drivers, and consequences of phenology

Plant phenology research has surged in recent decades, in part due to interest in phenological sensitivity to climate change and the vital role phenology plays in ecology. Many local-scale studies have generated important findings regarding the physiology, responses, and risks associated with shifts in plant phenology. By comparison, our understanding of regional- and global-scale phenology has been largely limited to remote sensing of green-up without the ability to differentiate among plant species. However, a new generation of analytical tools and data sources-including enhanced remote sensing products, digitized herbarium specimen data, and public participation in science-now permits investigating patterns and drivers of phenology across extensive taxonomic, temporal, and spatial scales, in an emerging field that we call macrophenology. Recent studies have highlighted how phenology affects dynamics at broad scales, including species interactions and ranges, carbon fluxes, and climate. At the cusp of this developing field of study, we review the theoretical and practical advances in four primary areas of plant macrophenology: (1) global patterns and shifts in plant phenology, (2) within-species changes in phenology as they mediate species' range limits and invasions at the regional scale, (3) broad-scale variation in phenology among species leading to ecological mismatches, and (4) interactions between phenology and global ecosystem processes. To stimulate future research, we describe opportunities for macrophenology to address grand challenges in each of these research areas, as well as recently available data sources that enhance and enable macrophenology research.

Estimating and Monitoring Land Surface Phenology in Rangelands: A Review of Progress and Challenges

Land surface phenology (LSP) has been extensively explored from global archives of satellite observations to track and monitor the seasonality of rangeland ecosystems in response to climate change. Long term monitoring of LSP provides large potential for the evaluation of interactions and feedbacks between climate and vegetation. With a special focus on the rangeland ecosystems, the paper reviews the progress, challenges and emerging opportunities in LSP while identifying possible gaps that could be explored in future. Specifically, the paper traces the evolution of satellite sensors and interrogates their properties as well as the associated indices and algorithms in estimating and monitoring LSP in productive rangelands. Findings from the literature revealed that the spectral characteristics of the early satellite sensors such as Landsat, AVHRR and MODIS played a critical role in the development of spectral vegetation indices that have been widely used in LSP applications. The normalized difference vegetation index (NDVI) pioneered LSP investigations, and most other spectral vegetation indices were primarily developed to address the weaknesses and shortcomings of the NDVI. New indices continue to be developed based on recent sensors such as Sentinel-2 that are characterized by unique spectral signatures and fine spatial resolutions, and their successful usage is catalyzed with the development of cutting-edge algorithms for modeling the LSP profiles. In this regard, the paper has documented several LSP algorithms that are designed to provide data smoothing, gap filling and LSP metrics retrieval methods in a single environment. In the future, the development of machine learning algorithms that can effectively model and characterize the phenological cycles of vegetation would help to unlock the value of LSP information in the rangeland monitoring and management process. Precisely, deep learning presents an opportunity to further develop robust software packages such as the decomposition and analysis of time series (DATimeS) with the abundance of data processing tools and techniques that can be used to better characterize the phenological cycles of vegetation in rangeland ecosystems.

Forecasting fire risk with machine learning and dynamic information derived from satellite vegetation index time-series

Fire risk mapping - mapping the probability of fire occurrence and spread - is essential for pre-fire management as well as for efficient firefighting efforts. Most fire risk maps are generated using static information on variables such as topography, vegetation density, and fuel instantaneous wetness. Satellites are often used to provide such information. However, long-term vegetation dynamics and the cumulative dryness status of the woody vegetation, which may affect fire occurrence and spread, are rarely considered in fire risk mapping. Here, we investigate the impact of two satellite-derived metrics that represent long-term vegetation status and dynamics on fire risk mapping - the long-term mean normalized difference vegetation index (NDVI) of the woody vegetation (NDVI_W) and its trend (NDVI_T). NDVI_W represents the mean woody density at the grid cell, while NDVI_T is the 5-year trend of the woody NDVI representing the long-term dryness status of the vegetation. To produce these metrics, we decompose time-series of satellite-derived NDVI following a method adjusted for Mediterranean woodlands and forests. We tested whether these metrics improve fire risk mapping using three machine learning (ML) algorithms (Logistic Regression, Random Forest, and XGBoost). We chose the 2007 wildfires in Greece for the analysis. Our results indicate that XGBoost, which accounts for variable interactions and non-linear effects, was the ML model that produced the best results. NDVI_W improved the model performance, while NDVI_T was significant only when NDVI_W was high. This NDVI_W-NDVI_T interaction means that the long-term dryness effect is meaningful only in places of dense woody vegetation. The proposed method can produce more accurate fire risk maps than conventional methods and can supply important dynamic information that may be used in fire behavior models.

Developing a New Method to Identify Flowering Dynamics of Rapeseed Using Landsat 8 and Sentinel-1/2

Identifying the rapeseed (Brassica napus L.) flowering dates are important for planting area estimation, growth monitoring, and yield estimation. However, there is currently a lack of data on rapeseed flowering dates at the parcel scale. In this study, a new spectral index (Normalized Rapeseed Flowering Index, NRFI) is proposed to detect rapeseed flowering dates from time series data generated from Landsat 8 OLI and Sentinel-2 sensors. This study also analyzed the feasibility of using the backscattering coefficients (VV, VH, and VV/VH) of Sentinel-1 to detect the flowering dates of rapeseed at the parcel scale. Based on the spectral and polarization characteristics of 718 rapeseed parcels collected in 2018, we developed a method to automatically identify peak flowering dates by the local maximum of NRFI series and the local minimum of VH and VV, along with the maximum of VV/VH. The results show that most of the peak flowering dates derived from Sentinel-1 and Sentinel-2 can be confirmed by the in-situ phenological observations at the Deutscher Wetterdienst (DWD) stations in Germany. The NRFI outperforms the Normalized Difference Yellow Index (NDYI) in identifying the peak flowering dates from Landsat 8. The derived medians of peak flowering dates by NRFI, NDYI (Sentinel-2), and VH are similar, while a systematic delay is observed by NDYI (Landsat 8). The method with the spectrum and backscattering coefficients will be a potential tool to identify crop flowering dynamics and map crop planting area.

Strengthening the evidence base for temperature-mediated phenological asynchrony and its impacts

Climate warming has caused the seasonal timing of many components of ecological food chains to advance. In the context of trophic interactions, the match-mismatch hypothesis postulates that differential shifts can lead to phenological asynchrony with negative impacts for consumers. However, at present there has been no consistent analysis of the links between temperature change, phenological asynchrony and individual-to-population-level impacts across taxa, trophic levels and biomes at a global scale. Here, we propose five criteria that all need to be met to demonstrate that temperature-mediated trophic asynchrony poses a growing risk to consumers. We conduct a literature review of 109 papers studying 129 taxa, and find that all five criteria are assessed for only two taxa, with the majority of taxa only having one or two criteria assessed. Crucially, nearly every study was conducted in Europe or North America, and most studies were on terrestrial secondary consumers. We thus lack a robust evidence base from which to draw general conclusions about the risk that climate-mediated trophic asynchrony may pose to populations worldwide.

A Novel Approach to Modelling Mangrove Phenology from Satellite Images: A Case Study from Northern Australia

Around the world, the effects of changing plant phenology are evident in many ways: from earlier and longer growing seasons to altering the relationships between plants and their natural pollinators. Plant phenology is often monitored using satellite images and parametric methods. Parametric methods assume that ecosystems have unimodal phenologies and that the phenology model is invariant through space and time. In evergreen ecosystems such as mangrove forests, these assumptions may not hold true. Here we present a novel, data-driven approach to extract plant phenology from Landsat imagery using Generalized Additive Models (GAMs). Using GAMs, we created models for six different mangrove forests across Australia. In contrast to parametric methods, GAMs let the data define the shape of the phenological curve, hence showing the unique characteristics of each study site. We found that the Enhanced Vegetation Index (EVI) model is related to leaf production rate (from in situ data), leaf gain and net leaf production (from the published literature). We also found that EVI does not respond immediately to leaf gain in most cases, but has a two- to three-month lag. We also identified the start of season and peak growing season dates at our field site. The former occurs between September and October and the latter May and July. The GAMs allowed us to identify dual phenology events in our study sites, indicated by two instances of high EVI and two instances of low EVI values throughout the year. We contribute to a better understanding of mangrove phenology by presenting a data-driven method that allows us to link physical changes of mangrove forests with satellite imagery. In the future, we will use GAMs to (1) relate phenology to environmental variables (e.g., temperature and rainfall) and (2) predict phenological changes.

Improved Estimates of Arctic Land Surface Phenology Using Sentinel-2 Time Series

The high spatial resolution and revisit time of Sentinel-2A/B tandem satellites allow a potentially improved retrieval of land surface phenology (LSP). The biome and regional characteristics, however, greatly constrain the design of the LSP algorithms. In the Arctic, such biome-specific characteristics include prolonged periods of snow cover, persistent cloud cover, and shortness of the growing season. Here, we evaluate the feasibility of Sentinel-2 for deriving high-resolution LSP maps of the Arctic. We extracted the timing of the start and end of season (SoS and EoS, respectively) for the years 2019 and 2020 with a simple implementation of the threshold method in Google Earth Engine (GEE). We found a high level of similarity between Sentinel-2 and PhenoCam metrics; the best results were observed with Sentinel-2 enhanced vegetation index (EVI) (root mean squared error (RMSE) and mean error (ME) of 3.0 d and −0.3 d for the SoS, and 6.5 d and −3.8 d for the EoS, respectively), although other vegetation indices presented similar performances. The phenological maps of Sentinel-2 EVI compared well with the same maps extracted from the Moderate Resolution Imaging Spectroradiometer (MODIS) in homogeneous landscapes (RMSE and ME of 9.2 d and 2.9 d for the SoS, and 6.4 and −0.9 d for the EoS, respectively). Unreliable LSP estimates were filtered and a quality flag indicator was activated when the Sentinel-2 time series presented a long period (>40 d) of missing data; discontinuities were lower in spring and early summer (9.2%) than in late summer and autumn (39.4%). The Sentinel-2 high-resolution LSP maps and the GEE phenological extraction method will support vegetation monitoring and contribute to improving the representation of Artic vegetation phenology in land surface models.

Applications of the Google Earth Engine and Phenology-Based Threshold Classification Method for Mapping Forest Cover and Carbon Stock Changes in Siem Reap Province, Cambodia

Digital and scalable technologies are increasingly important for rapid and large-scale assessment and monitoring of land cover change. Until recently, little research has existed on how these technologies can be specifically applied to the monitoring of Reducing Emissions from Deforestation and Forest Degradation (REDD+) activities. Using the Google Earth Engine (GEE) cloud computing platform, we applied the recently developed phenology-based threshold classification method (PBTC) for detecting and mapping forest cover and carbon stock changes in Siem Reap province, Cambodia, between 1990 and 2018. The obtained PBTC maps were validated using Google Earth high resolution historical imagery and reference land cover maps by creating 3771 systematic 5 × 5 km spatial accuracy points. The overall cumulative accuracy of this study was 92.1% and its cumulative Kappa was 0.9, which are sufficiently high to apply the PBTC method to detect forest land cover change. Accordingly, we estimated the carbon stock changes over a 28-year period in accordance with the Good Practice Guidelines of the Intergovernmental Panel on Climate Change. We found that 322,694 ha of forest cover was lost in Siem Reap, representing an annual deforestation rate of 1.3% between 1990 and 2018. This loss of forest cover was responsible for carbon emissions of 143,729,440 MgCO2 over the same period. If REDD+ activities are implemented during the implementation period of the Paris Climate Agreement between 2020 and 2030, about 8,256,746 MgCO2 of carbon emissions could be reduced, equivalent to about USD 6-115 million annually depending on chosen carbon prices. Our case study demonstrates that the GEE and PBTC method can be used to detect and monitor forest cover change and carbon stock changes in the tropics with high accuracy.

Status of Phenological Research Using Sentinel-2 Data: A Review

Remote sensing of plant phenology as an indicator of climate change and for mapping land cover has received significant scientific interest in the past two decades. The advancing of spring events, the lengthening of the growing season, the shifting of tree lines, the decreasing sensitivity to warming and the uniformity of spring across elevations are a few of the important indicators of trends in phenology. The Sentinel-2 satellite sensors launched in June 2015 (A) and March 2017 (B), with their high temporal frequency and spatial resolution for improved land mapping missions, have contributed significantly to knowledge on vegetation over the last three years. However, despite the additional red-edge and short wave infra-red (SWIR) bands available on the Sentinel-2 multispectral instruments, with improved vegetation species detection capabilities, there has been very little research on their efficacy to track vegetation cover and its phenology. For example, out of approximately every four papers that analyse normalised difference vegetation index (NDVI) or enhanced vegetation index (EVI) derived from Sentinel-2 imagery, only one mentions either SWIR or the red-edge bands. Despite the short duration that the Sentinel-2 platforms have been operational, they have proved their potential in a wide range of phenological studies of crops, forests, natural grasslands, and other vegetated areas, and in particular through fusion of the data with those from other sensors, e.g., Sentinel-1, Landsat and MODIS. This review paper discusses the current state of vegetation phenology studies based on the first five years of Sentinel-2, their advantages, limitations, and the scope for future developments.

Using Landsat satellites to assess the impact of check dams built across erosive gullies on vegetation rehabilitation

Gully erosion, a process of soil removal due to water accumulation and runoff, is a worldwide problem affecting agricultural lands. Building check dams perpendicular to the flow direction is one of the suggested control practices to stabilize this process. Though there are many studies on the effect of erosive controls on land stabilization, few examine its effect on the rehabilitation of vegetation. Here we use information from the satellites Landsat-7 (1999-2018) and Landsat-8 (2013-2018) to assess the effect of soil check dams built during 2012 across three gullies with distinct structures in a dryland area on vegetative cover and water status. We use a time series analysis technique to decompose Landsat-derived soil adjusted vegetation index (SAVI) into woody (SAVI_W) and herbaceous (iSAVI_H) contributions. The integral over the seasonal signal of the normalized difference water index (iNDWI) was used to assess changes in water status in the gully. We used herbaceous biomass collected in the field in 2014-2017 to validate iSAVI_H as a proxy of herbaceous biomass. Our results show that following the construction of the check dams, the change in woody vegetation cover is best described by a sigmoid model with an increase of ~57% (95% CI: 39%-76%; p < 0.0001), while the herbaceous vegetation increases linearly at a rate of ~71% per year (95% CI: 48%-93% y^-1; p < 0.0001). The correlation between iSAVI_H and herbaceous biomass (R² = 0.56; n = 16; p < 0.001) corroborates this increase. We found higher herbaceous productivity in the deeper gully compared to the shallower gullies but not statistically different increase rates. An increase in iNDWI of ~68% (95% CI: 43%-95%; p < 0.0001) likely implies an improved water infiltration rate that favored the vegetation expansion. Our satellite-based approach can be used to assess the impact of erosive control practices on vegetation rehabilitation in heterogeneous gullies.

Self-adapting extraction of cropland phenological transitions of rotation agroecosystems using dynamically fused NDVI images

Remote sensing can be used to monitor cropland phenological characteristics; however, tradeoffs between the spatial and temporal resolutions of cloudless satellite images limit the accuracy of their retrieval. In this study, an improved enhanced spatial and temporal adaptive reflectance fusion model (ESTARFM) was applied to human-dominated Xiong'an New Area to develop a self-adapting algorithm automating the extraction of main phenological transition points (greenup, maturity, senescence, and dormancy). The analyses of cropland phenological characteristics were performed utilizing the Softmax classification method. By examining three different phases of fusion images, it was found that the improved ESTARFM was more accurate than the original ESTARFM (correlation coefficient > 0.76; relative root mean square error < 0.25; structural similarity index > 0.79). The reconstructed normalized difference vegetation indexes were consistent with that acquired by the Moderate Resolution Imaging Spectroradiometer (average discrepancy: 0.1136, median absolute deviation: 0.0110). The greenup, maturity, senescence, and dormancy points were monitored in 5-day resolution and 50-day length on a 30-m grid scale, and their average day of year (DOY) were 67, 119, 127, and 166 for wheat; 173, 224, 232, and 283 for single-season corn; and 189, 227, 232, and 285 for rotation corn, respectively. The corresponding median absolute deviations were 2, 3, 2, and 2 days for wheat; 2, 5, 3, and 4 days for single-season corn; and 2, 5, 2, and 2 days for rotation corn, respectively, while all coefficients of variation did not exceed 6%. The proposed self-adapting approach can be used for identifying the planting modes at grid level in rotation agroecosystems and cropland phenological dynamics on a global or regional scale.

Detecting flowering phenology in oil seed rape parcels with Sentinel-1 and -2 time series

A novel methodology is proposed to robustly map oil seed rape (OSR) flowering phenology from time series generated from the Copernicus Sentinel-1 (S1) and Sentinel-2 (S2) sensors. The time series are averaged at parcel level, initially for a set of 229 reference parcels for which multiple phenological observations on OSR flowering have been collected from April 21 to May 19, 2018. The set of OSR parcels is extended to a regional sample of 32,355 OSR parcels derived from a regional S2 classification. The study area comprises the northern Brandenburg and Mecklenburg-Vorpommern (N) and the southern Bavaria (S) regions in Germany. A method was developed to automatically compute peak flowering at parcel level from the S2 time signature of the Normalized Difference Yellow Index (NDYI) and from the local minimum in S1 VV polarized backscattering coefficients. Peak flowering was determined at a temporal accuracy of 1 to 4 days. A systematic flowering delay of 1 day was observed in the S1 detection compared to S2. Peak flowering differed by 12 days between the N and S. Considerable local variation was observed in the N-S parcel-level flowering gradient. Additional in-situ phenology observations at 70 Deutscher Wetterdienst (DWD) stations confirm the spatial and temporal consistency between S1 and S2 signatures and flowering phenology across both regions. Conditions during flowering strongly determine OSR yield, therefore, the capacity to continuously characterize spatially the timing of key flowering dates across large areas is key. To illustrate this, expected flowering dates were simulated assuming a single OSR variety with a 425 growing degree days (GDD) requirement to reach flowering. This GDD requirement was calculated based on parcel-level peak flowering dates and temperatures accumulated from 25-km gridded meteorological data. The correlation between simulated and S2 observed peak flowering dates still equaled 0.84 and 0.54 for the N and S respectively. These Sentinel-based parcel-level flowering parameters can be combined with weather data to support in-season predictions of OSR yield, area, and production. Our approach identified the unique temporal signatures of S1 and S2 associated with OSR flowering and can now be applied to monitor OSR phenology for parcels across the globe.

Agricultural productivity in relation to climate and cropland management in West Africa

The climate of West Africa is expected to become more arid due to increased temperature and uncertain rainfall regimes, while its population is expected to grow faster than the rest of the world. As such, increased demand for food will likely coincide with declines in agricultural production in a region where severe undernutrition already occurs. Here, we attempt to discriminate between the impacts of climate and other factors (e.g. land management/degradation) on crop production across West Africa using satellite remote sensing. We identify trends in the land surface phenology and climate of West African croplands between 2000 and 2018. Using the combination of a an attribution framework and residual trend anlaysis, we discriminate between climate and other impacts on crop productivity. The combined effect of rainfall, land surface temperature and solar radiation explains approximately 40% of the variation in cropland productivity over West Africa at the 95% significance level. The largest proportions of croplands with greening trends were observed in Mali, Niger and Burkina Faso, and the largest proportions with browning trends were in Nigeria, The Gambia and Benin. Climate was responsible for 52% of the greening trends and 25% of the browning trends. Within the other driving factors, changes in phenology explained 18% of the greening and 37% of the browning trends across the region, the use of inputs and irrigation explained 30% of the greening trends and land degradation 38% of the browning trends. These findings have implications for adaptation policies as we map out areas in need of improved land management practices and those where it has proven to be successful.

The polar regions in a 2°C warmer world

Over the past decade, the Arctic has warmed by 0.75°C, far outpacing the global average, while Antarctic temperatures have remained comparatively stable. As Earth approaches 2°C warming, the Arctic and Antarctic may reach 4°C and 2°C mean annual warming, and 7°C and 3°C winter warming, respectively. Expected consequences of increased Arctic warming include ongoing loss of land and sea ice, threats to wildlife and traditional human livelihoods, increased methane emissions, and extreme weather at lower latitudes. With low biodiversity, Antarctic ecosystems may be vulnerable to state shifts and species invasions. Land ice loss in both regions will contribute substantially to global sea level rise, with up to 3 m rise possible if certain thresholds are crossed. Mitigation efforts can slow or reduce warming, but without them northern high latitude warming may accelerate in the next two to four decades. International cooperation will be crucial to foreseeing and adapting to expected changes.

Text Mining in Remotely Sensed Phenology Studies: A Review on Research Development, Main Topics, and Emerging Issues

As an interdisciplinary field of research, phenology is developing rapidly, and the contents of phenological research have become increasingly abundant. In addition, the potentiality of remote sensing technologies has largely contributed to the growth and complexity of this discipline, in terms of the scale of analysis, techniques of data processing, and a variety of topics. As a consequence, it is increasingly difficult for scientists to get a clear picture of remotely sensed phenology (rs+pheno) research. Bibliometric analysis is increasingly used for the study of a discipline and its conceptual dynamics. This review analyzed the last 40 years (1979–2018) of publications in the rs+pheno field retrieved from the Scopus database; such publications were investigated by means of a text mining approach, both in terms of bibliographic and text data. Results demonstrated that rs+pheno research is exponentially growing through time; however, it is primarily considered a subset of remote sensing science rather than a branch of phenology. In this framework, in the last decade, agriculture is becoming more and more a standalone science in rs+pheno research, independently from other related topics, e.g., classification. On the contrary, forestry struggles to gain its thematic role in rs+pheno studies and remains strictly connected with climate change issues. Classification and mapping represent the major rs+pheno topic, together with the extraction and the analysis of phenological metrics, like the start of the growing season. To the contrary, forest ecophysiology, in terms of ecosystem respiration and net ecosystem exchange, results as the most relevant new topic, together with the use of the red edge band and SAR (Synthetic Aperture Radar) data in rs+pheno agricultural studies. Some niche emerging rs+pheno topics may be recognized in the ocean and arctic investigations linked to phytoplankton blooming and ice cover dynamics. The findings of this study might be applicable for planning and managing remotely sensed phenology research; scientists involved in such discipline might use this study as a reference to consider their research domain in a broader dynamical network.

Remotely Sensed Water Limitation in Vegetation: Insights from an Experiment with Unmanned Aerial Vehicles (UAVs)

Unmanned aerial vehicles (UAVs) equipped with multispectral sensors present an opportunity to monitor vegetation with on-demand high spatial and temporal resolution. In this study we use multispectral imagery from quadcopter UAVs to monitor the progression of a water manipulation experiment on a common shrub, Baccharis pilularis (coyote brush) at the Blue Oak Ranch Reserve (BORR) ~20 km east of San Jose, California. We recorded multispectral imagery at several altitudes with nearly hourly intervals to explore the relationship between two common spectral indices, NDVI (normalized difference vegetation index) and NDRE (normalized difference red edge index), leaf water content and water potential as physiological metrics of plant water status, across a gradient of water deficit. An examination of the spatial and temporal thresholds at which water limitations were most detectable revealed that the best separation between levels of water deficit were at higher resolution (lower flying height), and in the morning (NDVI) and early morning (NDRE). We found that both measures were able to identify moisture deficit across treatments; however, NDVI was better able to distinguish between treatments than NDRE and was more positively correlated with field measurements of leaf water content. Finally, we explored how relationships between spectral indices and water status changed when the imagery was scaled to courser resolutions provided by satellite-based imagery (PlanetScope).We found that PlanetScope data was able to capture the overall trend in treatments but unable to capture subtle changes in water content. These kinds of experiments that evaluate the relationship between direct field measurements and UAV camera sensitivity are needed to enable translation of field-based physiology measurements to landscape or regional scales.

Determination of Vegetation Thresholds for Assessing Land Use and Land Use Changes in Cambodia using the Google Earth Engine Cloud-Computing Platform

As more data and technologies become available, it is important that a simple method is developed for the assessment of land use changes because of the global need to understand the potential climate mitigation that could result from a reduction in deforestation and forest degradation in the tropics. Here, we determined the threshold values of vegetation types to classify land use categories in Cambodia through the analysis of phenological behaviors and the development of a robust phenology-based threshold classification (PBTC) method for the mapping and long-term monitoring of land cover changes. We accessed 2199 Landsat collections using Google Earth Engine (GEE) and applied the Enhanced Vegetation Index (EVI) and harmonic regression methods to identify phenological behaviors of land cover categories during the leaf-shedding phenology (LSP) and leaf-flushing phenology (LFS) seasons. We then generated 722 mean phenology EVI profiles for 12 major land cover categories and determined the threshold values for selected land cover categories in the mid-LSP season. The PBTC pixel-based classified map was validated using very high-resolution (VHR) imagery. We obtained a cumulative overall accuracy of more than 88% and a cumulative overall accuracy of the referenced forest cover of almost 85%. These high accuracy values suggest that the very first PBTC map can be useful for estimating the activity data, which are critically needed to assess land use changes and related carbon emissions under the Reducing Emissions from Deforestation and forest Degradation (REDD+) scheme. We found that GEE cloud-computing is an appropriate tool to use to access remote sensing big data at scale and at no cost.

Using Hidden Markov Models for Land Surface Phenology: An Evaluation Across a Range of Land Cover Types in Southeast Spain

Land Surface Phenology (LSP) metrics are increasingly being used as indicators of climate change impacts in ecosystems. For this purpose, it is necessary to use methods that can be applied to large areas with different types of vegetation, including vulnerable semiarid ecosystems that exhibit high spatial variability and low signal-to-noise ratio in seasonality. In this work, we evaluated the use of hidden Markov models (HMM) to extract phenological parameters from Moderate Resolution Imaging Spectroradiometer (MODIS) derived Normalized Difference Vegetation Index (NDVI). We analyzed NDVI time-series data for the period 2000–2018 across a range of land cover types in Southeast Spain, including rice croplands, shrublands, mixed pine forests, and semiarid steppes. Start of Season (SOS) and End of Season (EOS) metrics derived from HMM were compared with those obtained using well-established smoothing methods. When a clear and consistent seasonal variation was present, as was the case in the rice croplands, and when adjusting average curves, the smoothing methods performed as well as expected, with HMM providing consistent results. When spatial variability was high and seasonality was less clearly defined, as in the semiarid shrublands and steppe, the performance of the smoothing methods degraded. In these cases, the results from HMM were also less consistent, yet they were able to provide pixel-wise estimations of the metrics even when comparison methods did not.

Seasonal patterns of canopy photosynthesis captured by remotely sensed sun-induced fluorescence and vegetation indexes in mid-to-high latitude forests: A cross-platform comparison

Characterized by the noticeable seasonal patterns of canopy photosynthesis, mid-to-high latitude forests are sensitive to climate change and crucial for understanding the global carbon cycle. To monitor the seasonal cycle of the canopy photosynthesis from space, several remotely sensed indexes, such as normalized difference vegetation index (NDVI), enhanced vegetation index (EVI), and leaf area index (LAI) have been implemented within the past decades. Recently, satellite-derived sun-induced fluorescence (SIF) has shown great potential of providing retrievals that are more related to photosynthesis process. However, the potentials of different canopy measurements have not been thoroughly assessed in the context of recent advances of new satellites and proposals of improved indexes. At 15 forested sites, we present a cross-platform intercomparison of one emerging remote sensing based index of phenology index (PI) and two SIF datasets against the conventional indexes such as NDVI, EVI, and LAI to capture the seasonal cycles of canopy photosynthesis. NDVI, EVI, LAI, and PI were calculated from Moderate Resolution Imaging Spectroradiometer (MODIS) measurements, while SIF were evaluated from Global Ozone Monitoring Experiment-2 (GOME-2) and Orbiting Carbon Observatory-2 (OCO-2) observations. Results indicated that GOME-2 SIF was highly correlated with gross primary production (GPP) and absorbed photosynthetically active radiation during the growing seasons. The SIF-GPP relationship can generally be considered linear at the 16-day scale. Key phenological metrics such as start of the seasons and end of the seasons captured by SIF from GOME-2 and OCO-2 matched closely with photosynthesis phenology as inferred by GPP. However, the applications of OCO-2 SIF for phenological studies may be limited only for a small range of sites (at site-level) due to a limited spatial sampling. Among the MODIS estimations, PI and NDVI provided most reliable predictions of start of growing seasons, while no indexes accurately captured the end of growing seasons.

Using Time Series of High-Resolution Planet Satellite Images to Monitor Grapevine Stem Water Potential in Commercial Vineyards

Spectral-based vegetation indices (VI) have been shown to be good proxies of grapevine stem water potential (Ψstem), assisting in irrigation decision-making for commercial vineyards. However, VI-Ψstem correlations are mostly reported at the leaf or canopy scales, using proximal canopy-based sensors or very-high-spatial resolution images derived from sensors mounted on small airplanes or drones. Here, for the first time, we take advantage of high-spatial resolution (3-m) near-daily images acquired from Planet’s nano-satellite constellation to derive VI-Ψstem correlations at the vineyard scale. Weekly Ψstem was measured along the growing season of 2017 in six vines each in 81 commercial vineyards and in 60 pairs of grapevines in a 2.4 ha experimental vineyard in Israel. The Clip application programming interface (API), provided by Planet, and the Google Earth Engine platform were used to derive spatially continuous time series of four VIs—GNDVI, NDVI, EVI and SAVI—in the 82 vineyards. Results show that per-week multivariable linear models using variables extracted from VI time series successfully tracked spatial variations in Ψstem across the experimental vineyard (Pearson’s-r = 0.45–0.84; N = 60). A simple linear regression model enabled monitoring seasonal changes in Ψstem along the growing season in the vineyard (r = 0.80–0.82). Planet VIs and seasonal Ψstem data from the 82 vineyards were used to derive a ‘global’ model for in-season monitoring of Ψstem at the vineyard-level (r = 0.78; RMSE = 18.5%; N = 970). The ‘global’ model, which requires only a few VI variables extracted from Planet images, may be used for real-time weekly assessment of Ψstem in Mediterranean vineyards, substantially improving the efficiency of conventional in-field monitoring efforts.

Economic Assessment of Fire Damage to Urban Forest in the Wildland–Urban Interface Using Planet Satellites Constellation Images

The wildland-urban interface (WUI)—the area where wildland vegetation and urban buildings intermix—is at a greater risk of fire occurrence because of extensive human activity in that area. Although satellite remote sensing has become a major tool for assessing fire damage in wildlands, it is unsuitable for WUI fire monitoring due to the low spatial resolution of the images from satellites that provide frequent information which is relevant for timely fire monitoring in WUI. Here, we take advantage of frequent (i.e., ca. daily), high-spatial-resolution (3 m) imagery acquired from a constellation of nano-satellites operated by Planet Labs (“Planet”) to assess fire damage to urban trees in the WUI of a Mediterranean city in Israel (Haifa). The fire occurred at the end of 2016, consuming ca. 17,000 of the trees (152 trees ha−1) within the near-by wildland and urban parts of the city. Three vegetation indices (GNDVI, NDVI and GCC) from Planet satellite images were used to derive a burn severity map for the WUI area after applying a subpixel discrimination method to distinguish between woody and herbaceous vegetation. The produced burn severity map was successfully validated with information acquired from an extensive field survey in the WUI burnt area (overall accuracy and kappa: 87% and 0.75%, respectively). Planet’s vegetation indices were calibrated using in-field tree measurements to obtain high spatial resolution maps of burned trees and consumed woody biomass in the WUI. These were used in conjunction with an ecosystem services valuation model (i-Tree) to estimate spatially-distributed and total economic loss due to damage to urban trees caused by the fire. Results show that nearly half of the urban trees were moderately and severely burned (26% and 22%, respectively). The total damage to the urban forest was estimated at ca. 41 ± 10 M USD. We conclude that using the method developed in this study with high-spatial-resolution Planet images has a great potential for WUI fire economic assessment.