Exploring the role of biotic factors in regulating the spatial variability in land surface phenology across four temperate forest sites.

Land surface phenology (LSP), the characterization of plant phenology with satellite data, is essential for understanding the effects of climate change on ecosystem functions. Considerable LSP variation is observed within local landscapes, and the role of biotic factors in regulating such variation remains underexplored. In this study, we selected four National Ecological Observatory Network terrestrial sites with minor topographic relief to investigate how biotic factors regulate intra-site LSP variability. We utilized plant functional type (PFT) maps, functional traits, and LSP data to assess the explanatory power of biotic factors for the start and end of season (SOS and EOS) variability. Our results indicate that PFTs alone explain only 0.8-23.4% of intra-site SOS and EOS variation, whereas including functional traits significantly improves explanatory power, with cross-validation correlations ranging from 0.50 to 0.85. While functional traits exhibited diverse effects on SOS and EOS across different sites, traits related to competitive ability and productivity were important for explaining both SOS and EOS variation at these sites. These findings reveal that plants exhibit diverse phenological responses to comparable environmental conditions, and functional traits significantly contribute to intra-site LSP variability, highlighting the importance of intrinsic biotic properties in regulating plant phenology.

Urbanization imprint on land surface phenology: The urban-rural gradient analysis for Chinese cities.

Rising temperature shifts plant phenology. Chinese cities, experiencing extensive expansion and intensive warming, spanning a wide latitudinal range, might provide ideal experimental opportunities for observing and predicting phenological responses to warming temperature. Using the urban-rural gradient approach, we explored urbanization imprint on land surface phenology across the entire urbanization intensity (UI) gradient ranging from 0% to 100% in 343 Chinese cities using the VIIRS Land Surface Phenology along with MODIS Land Surface Temperature (LST) products. We found prevalent advancing and delaying trends for the start of the growing season (SOS) and the end of the growing season (EOS) with increasing UI across 343 Chinese cities, respectively. Overall, the phenology shifted earlier by 8.6 ± 0.54 days for SOS, later by 1.3 ± 0.51 days for EOS, and lengthened by 9.9 ± 0.77 days for the growing season length (GSL) in urban core areas (UI above 50%) relative to their rural counterparts (UI lower than 1%). The temperature sensitivity of SOS and EOS was 10.5 ± 0.25 days earlier and 2.9 ± 0.16 days later per 1°C LST increase in spring and autumn, respectively. Moreover, the northern cities witnessed higher temperature sensitivity for SOS and EOS than the southern ones. Both spring and autumn temperature sensitivity across these 343 cities would likely decrease with future urban warming, suggesting any projections of future phenological responses to continued warming must be approached with caution.

Comparing land surface phenology of major European crops as derived from SAR and multispectral data of Sentinel-1 and -2.

The frequent acquisitions of fine spatial resolution imagery (10 m) offered by recent multispectral satellite missions, including Sentinel-2, can resolve single agricultural fields and thus provide crop-specific phenology metrics, a crucial information for crop monitoring. However, effective phenology retrieval may still be hampered by significant cloud cover. Synthetic aperture radar (SAR) observations are not restricted by weather conditions, and Sentinel-1 thus ensures more frequent observations of the land surface. However, these data have not been systematically exploited for phenology retrieval so far. In this study, we extracted crop-specific land surface phenology (LSP) from Sentinel-1 and Sentinel-2 of major European crops (common and durum wheat, barley, maize, oats, rape and turnip rape, sugar beet, sunflower, and dry pulses) using ground-truth information from the "Copernicus module" of the Land Use/Cover Area frame statistical Survey (LUCAS) of 2018. We consistently used a single model-fit approach to retrieve LSP metrics on temporal profiles of CR (Cross Ratio, the ratio of the backscattering coefficient VH/VV from Sentinel-1) and NDVI (Normalized Difference Vegetation Index from Sentinel-2). Our analysis revealed that LSP retrievals from Sentinel-1 are comparable to those of Sentinel-2, particularly for winter crops. The start of season (SOS) timings, as derived from Sentinel-1 and -2, are significantly correlated (average r of 0.78 for winter and 0.46 for summer crops). The correlation is lower for end of season retrievals (EOS, r of 0.62 and 0.34). Agreement between LSP derived from Sentinel-1 and -2 varies among crop types, ranging from r = 0.89 and mean absolute error MAE = 10 days (SOS of dry pulses) to r = 0.15 and MAE = 53 days (EOS of sugar beet). Observed deviations revealed that Sentinel-1 and -2 LSP retrievals can be complementary; for example for winter crops we found that SAR detected the start of the spring growth while multispectral data is sensitive to the vegetative growth before and during winter. To test if our results correspond reasonably to in-situ data, we compared average crop-specific LSP for Germany to average phenology from ground phenological observations of 2018 gathered from the German Meteorological Service (DWD). Our study demonstrated that both Sentinel-1 and -2 can provide relevant and at times complementary LSP information at field- and crop-level.

Biodiversity promotes primary productivity and growing season lengthening at the landscape scale.

Experiments have shown positive biodiversity-ecosystem functioning (BEF) relationships in small plots with model communities established from species pools typically comprising few dozen species. Whether patterns found can be extrapolated to complex, nonexperimental, real-world landscapes that provide ecosystem services to humans remains unclear. Here, we combine species inventories from a large-scale network of 447 1-km2 plots with remotely sensed indices of primary productivity (years 2000-2015). We show that landscape-scale productivity and its temporal stability increase with the diversity of plants and other taxa. Effects of biodiversity indicators on productivity were comparable in size to effects of other important drivers related to climate, topography, and land cover. These effects occurred in plots that integrated different ecosystem types (i.e., metaecosystems) and were consistent over vast environmental and altitudinal gradients. The BEF relations we report are as strong or even exceed the ones found in small-scale experiments, despite different community assembly processes and a species pool comprising nearly 2,000 vascular plant species. Growing season length increased progressively over the observation period, and this shift was accelerated in more diverse plots, suggesting that a large species pool is important for adaption to climate change. Our study further implies that abiotic global-change drivers may mediate ecosystem functioning through biodiversity changes.

Examining land surface phenology in the tropical moist forest eco-zone of South America.

Using leaf area index (LAI) data from 1981 to 2014 in the tropical moist forest eco-zone of South America, we extracted start (SOS) and end (EOS) dates of the active growing season in forest and savanna at each pixel. Then, we detected spatiotemporal characteristics of SOS and EOS in the two vegetation types. Moreover, we analyzed relationships between interannual variations of SOS/EOS and climatic factors, and simulated SOS/EOS time series based on preceding mean air temperature and accumulated rainfall. Results show that mean SOS and EOS ranged from 260 to 330 day of year (DOY) and from 150 to 260 DOY across the study region, respectively. From 1981 to 2014, SOS advancement is more extensive than SOS delay, while EOS advancement and delay are similarly extensive. For most pixels of forest and savanna in tropical moist forest eco-zone, preceding rainfall correlates predominantly negatively with SOS but positively with EOS, while the relationship between preceding temperature and phenophases is location-specific. In addition, preceding rainfall is more extensive than preceding temperature in simulating SOS, while both preceding rainfall and temperature play an important role for simulating EOS. This study highlights the reliability of using LAI data for long-term phenological analysis in the tropical moist forest eco-zone.

Monitoring Spatio-Temporal Changes of Terrestrial Ecosystem Soil Water Use Efficiency in Northeast China Using Time Series Remote Sensing Data.

Soil water use efficiency (SWUE) was proposed as an effective proxy of ecosystem water use efficiency (WUE), which reflects the coupling of the carbon⁻water cycle and function of terrestrial ecosystems. The changes of ecosystem SWUE at the regional scale and their relationships with the environmental and biotic factors are yet to be adequately understood. Here, we aim to estimate SWUE over northeast China using time-series Moderate Resolution Imaging Spectroradiometer (MODIS) gross primary productivity data and European Space Agency climate change initiative (ESA CCI) soil moisture product during 2007⁻2015. The spatio-temporal variations in SWUE and their linkages to multiple factors, especially the phenological metrics, were investigated using trend and correlation analysis. The results showed that the spatial heterogeneity of ecosystem SWUE in northeast China was obvious. SWUE distribution varied among vegetation types, soil types, and elevation. Forests might produce higher photosynthetic productivity by utilizing unit soil moisture. The seasonal variations of SWUE were consistent with the vegetation growth cycle. Changes in normalized difference vegetation index (NDVI), land surface temperature, and precipitation exerted positive effects on SWUE variations. The earlier start (SOS) and later end (EOS) of the growing season would contribute to the increase in SWUE. Our results help complement the knowledge of SWUE variations and their driving forces.

Deciduous forest responses to temperature, precipitation, and drought imply complex climate change impacts.

Changes in spring and autumn phenology of temperate plants in recent decades have become iconic bio-indicators of rapid climate change. These changes have substantial ecological and economic impacts. However, autumn phenology remains surprisingly little studied. Although the effects of unfavorable environmental conditions (e.g., frost, heat, wetness, and drought) on autumn phenology have been observed for over 60 y, how these factors interact to influence autumn phenological events remain poorly understood. Using remotely sensed phenology data from 2001 to 2012, this study identified and quantified significant effects of a suite of environmental factors on the timing of fall dormancy of deciduous forest communities in New England, United States. Cold, frost, and wet conditions, and high heat-stress tended to induce earlier dormancy of deciduous forests, whereas moderate heat- and drought-stress delayed dormancy. Deciduous forests in two eco-regions showed contrasting, nonlinear responses to variation in these explanatory factors. Based on future climate projection over two periods (2041-2050 and 2090-2099), later dormancy dates were predicted in northern areas. However, in coastal areas earlier dormancy dates were predicted. Our models suggest that besides warming in climate change, changes in frost and moisture conditions as well as extreme weather events (e.g., drought- and heat-stress, and flooding), should also be considered in future predictions of autumn phenology in temperate deciduous forests. This study improves our understanding of how multiple environmental variables interact to affect autumn phenology in temperate deciduous forest ecosystems, and points the way to building more mechanistic and predictive models.

Land surface phenology of Northeast China during 2000-2015: temporal changes and relationships with climate changes.

As an important crop growing area, Northeast China (NEC) plays a vital role in China's food security, which has been severely affected by climate change in recent years. Vegetation phenology in this region is sensitive to climate change, and currently, the relationship between the phenology of NEC and climate change remains unclear. In this study, we used a satellite-derived normalized difference vegetation index (NDVI) to obtain the temporal patterns of the land surface phenology in NEC from 2000 to 2015 and validated the results using ground phenology observations. We then explored the relationships among land surface phenology, temperature, precipitation, and sunshine hours for relevant periods. Our results showed that the NEC experienced great phenological changes in terms of spatial heterogeneity during 2000-2015. The spatial patterns of land surface phenology mainly changed with altitude and land cover type. In most regions of NEC, the start date of land surface phenology had advanced by approximately 1.0 days year-1, and the length of land surface phenology had been prolonged by approximately 1.0 days year-1 except for the needle-leaf and cropland areas, due to the warm conditions. We found that a distinct inter-annual variation in land surface phenology related to climate variables, even if some areas presented non-significant trends. Land surface phenology was coupled with climate variables and distinct responses at different combinations of temperature, precipitation, sunshine hours, altitude, and anthropogenic influence. These findings suggest that remote sensing and our phenology extracting methods hold great potential for helping to understand how land surface phenology is sensitive to global climate change.

Variability and evolution of global land surface phenology over the past three decades (1982-2012).

Monitoring land surface phenology (LSP) is important for understanding both the responses and feedbacks of ecosystems to the climate system, and for representing these accurately in terrestrial biosphere models. Moreover, by shedding light on phenological trends at a variety of scales, LSP provides the potential to fill the gap between traditional phenological (field) observations and the large-scale view of global models. In this study, we review and evaluate the variability and evolution of satellite-derived growing season length (GSL) globally and over the past three decades. We used the longest continuous record of Normalized Difference Vegetation Index data available to date at global scale to derive LSP metrics consistently over all vegetated land areas and for the period 1982-2012. We tested GSL, start- and end-of-season metrics (SOS and EOS, respectively) for linear trends as well as for significant trend shifts over the study period. We evaluated trends using global environmental stratification information in place of commonly used land cover maps to avoid circular findings. Our results confirmed an average lengthening of the growing season globally during 1982-2012 - averaging 0.22-0.34 days yr(-1), but with spatially heterogeneous trends. About 13-19% of global land areas displayed significant GSL change, and over 30% of trends occurred in the boreal/alpine biome of the Northern Hemisphere, which showed diverging GSL evolution over the past three decades. Within this biome, the 'Cold and Mesic' environmental zone appeared as an LSP change hotspot. We also examined the relative contribution of SOS and EOS to the overall changes, finding that EOS trends were generally stronger and more prevalent than SOS trends. These findings constitute a step towards the identification of large-scale phenological drivers of vegetated land surfaces, necessary for improving phenological representation in terrestrial biosphere models.

Comparison of the driving forces of spring phenology among savanna landscapes by including combined spatial and temporal heterogeneity.

Understanding spatial and temporal dynamics of land surface phenology (LSP) and its driving forces are critical for providing information relevant to short- and long-term decision making, particularly as it relates to climate response planning. With the third generation Global Inventory Monitoring and Modeling System (GIMMS3g) Normalized Difference Vegetation Index (NDVI) data and environmental data from multiple sources, we investigated the spatio-temporal changes in the start of the growing season (SOS) in southern African savannas from 1982 through 2010 and determined its linkage to environmental factors using spatial panel data models. Overall, the SOS occurs earlier in the north compared to the south. This relates in part to the differences in ecosystems, with northern areas representing high rainfall and dense tree cover (mainly tree savannas), whereas the south has lower rainfall and sparse tree cover (mainly bush and grass savannas). From 1982 to 2010, an advanced trend was observed predominantly in the tree savanna areas of the north, whereas a delayed trend was chiefly found in the floodplain of the north and bush/grass savannas of the south. Different environmental drivers were detected within tree- and grass-dominated savannas, with a critical division being represented by the 800 mm isohyet. Our results supported the importance of water as a driver in this water-limited system, specifically preseason soil moisture, in determining the SOS in these water-limited, grass-dominated savannas. In addition, the research pointed to other, often overlooked, effects of preseason maximum and minimum temperatures on the SOS across the entire region. Higher preseason maximum temperatures led to an advance of the SOS, whereas the opposite effects of preseason minimum temperature were observed. With the rapid increase in global change research, this work will prove helpful for managing savanna landscapes and key to predicting how projected climate changes will affect regional vegetation phenology and productivity.

Strong contribution of autumn phenology to changes in satellite-derived growing season length estimates across Europe (1982-2011).

Land Surface Phenology (LSP) is the most direct representation of intra-annual dynamics of vegetated land surfaces as observed from satellite imagery. LSP plays a key role in characterizing land-surface fluxes, and is central to accurately parameterizing terrestrial biosphere-atmosphere interactions, as well as climate models. In this article, we present an evaluation of Pan-European LSP and its changes over the past 30 years, using the longest continuous record of Normalized Difference Vegetation Index (NDVI) available to date in combination with a landscape-based aggregation scheme. We used indicators of Start-Of-Season, End-Of-Season and Growing Season Length (SOS, EOS and GSL, respectively) for the period 1982-2011 to test for temporal trends in activity of terrestrial vegetation and their spatial distribution. We aggregated pixels into ecologically representative spatial units using the European Landscape Classification (LANMAP) and assessed the relative contribution of spring and autumn phenology. GSL increased significantly by 18-24 days decade(-1) over 18-30% of the land area of Europe, depending on methodology. This trend varied extensively within and between climatic zones and landscape classes. The areas of greatest growing-season lengthening were the Continental and Boreal zones, with hotspots concentrated in southern Fennoscandia, Western Russia and pockets of continental Europe. For the Atlantic and Steppic zones, we found an average shortening of the growing season with hotspots in Western France, the Po valley, and around the Caspian Sea. In many zones, changes in the NDVI-derived end-of-season contributed more to the GSL trend than changes in spring green-up, resulting in asymmetric trends. This underlines the importance of investigating senescence and its underlying processes more closely as a driver of LSP and global change.

A satellite-derived dataset on vegetation phenology across Central Asia from 2001 to 2023.

Satellite-observed land surface phenology (LSP) data have helped us better understand terrestrial ecosystem dynamics at large scales. However, uncertainties remain in comprehending LSP variations in Central Asian drylands. In this article, an LSP dataset covering Central Asia (45-100°E, 33-57°N) is introduced. This LSP dataset was produced based on Moderate Resolution Imaging Spectroradiometer (MODIS) 0.05-degree daily reflectance and land cover data. The phenological dynamics of drylands were tracked using the seasonal profiles of near-infrared reflectance of vegetation (NIRv). NIRv time series processing involved the following steps: identifying low-quality observations, smoothing the NIRv time series, and retrieving LSP metrics. In the smoothing step, a median filter was first applied to reduce spikes, after which the stationary wavelet transform (SWT) was used to smooth the NIRv time series. The SWT was performed using the Biorthogonal 1.1 wavelet at a decomposition level of 5. Seven LSP metrics were provided in this dataset, and they were categorized into the following three groups: (1) timing of key phenological events, (2) NIRv values essential for the detection of the phenological events throughout the growing season, and (3) NIRv value linked to vegetation growth state during the growing season. This LSP dataset is useful for investigating dryland ecosystem dynamics in response to climate variations and human activities across Central Asia.

Comparison of change-based and shape-based data fusion methods in fine-resolution land surface phenology monitoring with Landsat and Sentinel-2 data.

Fine-resolution land surface phenology (LSP) is urgently required for applications on agriculture management and vegetation-climate interaction, especially over heterogeneous areas, such as agricultural lands and fragmented forests. The critical challenge of fine-resolution LSP monitoring is how to reconstruct the spatiotemporal continuous vegetation index time series. To solve this problem, various data fusion methods have been devised; however, the comprehensive inter-comparison is lacking across different spatial heterogeneity, data quality, and vegetation types. We divide these methods into two main categories: the change-based methods fusing satellite observations with different spatiotemporal resolutions, and the shape-based methods fusing prior knowledge of shape models and satellite observations. We selected four methods to rebuilt two-band enhanced vegetation index (EVI2) series based on the harmonized Landsat and Sentinel-2 (HLS) data, including two change-based methods, namely the Spatial and temporal Adaptive Reflectance Fusion Model (STARFM), the Flexible Spatiotemporal DAta Fusion (FSDAF), and two shape-based methods, namely the Multiple-year Weighting Shape-Matching (MWSM), and the Spatiotemporal Shape-Matching Model (SSMM). Four phenological transition dates were extracted, evaluated with PhenoCam observations and the 500 m Visible Infrared Imaging Radiometer Suite (VIIRS) phenology product. The 30 m transition dates show more spatial details and reveal more apparent intra-class and inter-class phenology variation compared with 500 m product. The four transition dates of SSMM and FSDAF (R2>0.74, MAD<15 days) show better agreement with PhenoCam-derived dates. The performance difference between fusion methods over various application scenarios are then analyzed. Fusion results are more robust when temporal frequency is higher than 15 observations per year. The shape-based methods are less sensitive to temporal sampling irregularity than change-based methods. Both change-based methods and shape-based methods cannot perform well when the region is heterogeneous. Among different vegetation types, SSMM-like methods have the highest overall accuracy. The findings in this paper can provide references for regional and global fine-resolution phenology monitoring.

Revegetation affects the response of land surface phenology to climate in Loess Plateau, China.

Land surface phenology (LSP), defined as the plant's growth rhythm retrieved from satellite sensing products, is proven to shift with climate change and affect the carbon cycles of terrestrial ecosystems. Global afforested area is largely increasing and consequently affecting local and global climate. However, how and to what extent revegetation affects LSP remains relatively unexplored. Here we investigated the difference in four LSPs (i.e., greenup, maturity, senescence, and dormancy) and the response of LSP to climate between restored and native vegetation on Loess Plateau, China, where a remarkable process of vegetation restoration happened during 1982-2015. Most study regions showed a longer growing season (LOS) over time, specifically, with a slight delay in greenup but a relatively large delay in senescence. We found that air temperature was the dominant factor affecting greenup and maturity, while precipitation mostly controlled the senescence and dormancy in the study area. Under similar climate conditions, the LSP of restored vegetation (i.e., restored forest and grassland) showed a significant difference (p < 0.05) from native ones during 1999-2015. Compared to the native forest, restored forest from cropland and grassland showed a delayed greenup date by 0.3 and 3.6 days (p < 0.05) and an advanced dormancy date of 6.6 and 9.0 days (p < 0.05), respectively. Furthermore, the restored vegetation became less sensitive to air temperature than native vegetation, while the restored forest was more sensitive to precipitation, and its growth was affected by the water limitation to a larger extent in the study area. Our study highlights the necessity of considering land use management and its effect on the LSP change to better understand the effect of afforestation on global climate and carbon cycles.

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.

Multi-Season Phenology Mapping of Nile Delta Croplands Using Time Series of Sentinel-2 and Landsat 8 Green LAI.

Space-based cropland phenology monitoring substantially assists agricultural managing practices and plays an important role in crop yield predictions. Multitemporal satellite observations allow analyzing vegetation seasonal dynamics over large areas by using vegetation indices or by deriving biophysical variables. The Nile Delta represents about half of all agricultural lands of Egypt. In this region, intensifying farming systems are predominant and multi-cropping rotations schemes are increasing, requiring a high temporal and spatial resolution monitoring for capturing successive crop growth cycles. This study presents a workflow for cropland phenology characterization and mapping based on time series of green Leaf Area Index (LAI) generated from NASA's Harmonized Landsat 8 (L8) and Sentinel-2 (S2) surface reflectance dataset from 2016 to 2019. LAI time series were processed for each satellite dataset, which were used separately and combined to identify seasonal dynamics for a selection of crop types (wheat, clover, maize and rice). For the combination of L8 with S2 LAI products, we proposed two time series smoothing and fitting methods: (1) the Savitzky-Golay (SG) filter and (2) the Gaussian Processes Regression (GPR) fitting function. Single-sensor and L8-S2 combined LAI time series were used for the calculation of key crop Land Surface Phenology (LSP) metrics (start of season, end of season, length of season), whereby the detection of cropland growing seasons was based on two established threshold methods, i.e., a seasonal or a relative amplitude value. Overall, the developed phenology extraction scheme enabled identifying up to two successive crop cycles within a year, with a superior performance observed for the seasonal than for the relative threshold method, in terms of consistency and cropland season detection capability. Differences between the time series collections were analyzed by comparing the phenology metrics per crop type and year. Results suggest that L8-S2 combined LAI data streams with GPR led to a more precise detection of the start and end of growing seasons for most crop types, reaching an overall detection of 74% over the total planted crops versus 69% with S2 and 63% with L8 alone. Finally, the phenology mapping allowed us to evaluate the spatial and temporal evolution of the croplands over the agroecosystem in the Nile Delta.

Monitoring Cropland Phenology on Google Earth Engine Using Gaussian Process Regression.

Monitoring cropland phenology from optical satellite data remains a challenging task due to the influence of clouds and atmospheric artifacts. Therefore, measures need to be taken to overcome these challenges and gain better knowledge of crop dynamics. The arrival of cloud computing platforms such as Google Earth Engine (GEE) has enabled us to propose a Sentinel-2 (S2) phenology end-to-end processing chain. To achieve this, the following pipeline was implemented: (1) the building of hybrid Gaussian Process Regression (GPR) retrieval models of crop traits optimized with active learning, (2) implementation of these models on GEE (3) generation of spatiotemporally continuous maps and time series of these crop traits with the use of gap-filling through GPR fitting, and finally, (4) calculation of land surface phenology (LSP) metrics such as the start of season (SOS) or end of season (EOS). Overall, from good to high performance was achieved, in particular for the estimation of canopy-level traits such as leaf area index (LAI) and canopy chlorophyll content, with normalized root mean square errors (NRMSE) of 9% and 10%, respectively. By means of the GPR gap-filling time series of S2, entire tiles were reconstructed, and resulting maps were demonstrated over an agricultural area in Castile and Leon, Spain, where crop calendar data were available to assess the validity of LSP metrics derived from crop traits. In addition, phenology derived from the normalized difference vegetation index (NDVI) was used as reference. NDVI not only proved to be a robust indicator for the calculation of LSP metrics, but also served to demonstrate the good phenology quality of the quantitative trait products. Thanks to the GEE framework, the proposed workflow can be realized anywhere in the world and for any time window, thus representing a shift in the satellite data processing paradigm. We anticipate that the produced LSP metrics can provide meaningful insights into crop seasonal patterns in a changing environment that demands adaptive agricultural production.

Land surface phenology and greenness in Alpine grasslands driven by seasonal snow and meteorological factors.

Snow accumulation and melt have multiple impacts on Land Surface Phenology (LSP) and greenness in Alpine grasslands. Our understanding of these impacts and their interactions with meteorological factors are still limited. In this study, we investigate this topic by analyzing LSP dynamics together with potential drivers, using satellite imagery and other data sources. LSP (start and end of season) and greenness metrics were extracted from time series of vegetation and leaf area index. As explanatory variables we used snow accumulation, snow cover melt date and meteorological factors. We tested for inter-annual co-variation of LSP and greenness metrics with seasonal snow and meteorological metrics across elevations and for four sub-regions of natural grasslands in the Swiss Alps over the period 2003-2014. We found strong positive correlations of snow cover melt date and snow accumulation with the start of season, especially at higher elevation. Autumn temperature was found to be important at the end of season below 2000 m above sea level (m asl), while autumn precipitation was relevant above 2000 m asl, indicating climatic growth limiting factors to be elevation dependent. The effects of snow and meteorological factors on greenness revealed that this metric tends to be influenced by temperatures at high elevations, and by snow melt date at low elevations. Given the high sensitivity of alpine grassland ecosystems, these results suggest that alpine grasslands may be particularly affected by future changes in seasonal snow, to varying degree depending on elevation.

Detecting distant sources of airborne pollen for Poland: Integrating back-trajectory and dispersion modelling with a satellite-based phenology.

Airborne pollen might be transported over thousands of kilometres, which has important ecological, evolutionary and clinical consequences. The long-distance transport (LDT) of birch (Betula sp.) pollen has been described in detail for northern Europe. However, a comprehensive analysis of this transport from other European regions is lacking. This study focused on the post-seasonal LDT of birch pollen to Poland (central Europe), with special attention paid to determining potential source areas of pollen and describing the causal mechanism favouring LDT episodes. Pollen monitoring (1997-2016) was conducted in Poznan and Rzeszów (500 km away from each other) using volumetric traps. The LDT episodes were characterized by analysing the (1) bi-hourly backward air mass trajectories using the Hybrid Single Particle Lagrangian Integrated Trajectory model (HYSPLIT); (2) sea level pressure (SLP) and 500 hPa geopotential height (z500) anomalies; and (3) patterns of the Enhanced Vegetation Index to determine the birch flowering time along the moving air mass trajectories. The potential locations of birch populations within broadleaved forests were estimated with GLOBCOVER data. Finally, the movement of pollen emitted from potential source areas was simulated using the HYSPLIT dispersion model. LDT episodes were mainly recorded in the first fortnight of May. The main source areas of pollen to Poland were western Russia, Belarus and to a lesser extent the eastern Baltic republics and the Scandinavian Peninsula. In most cases, a high-pressure centre located over Scandinavia and an elevated z500 over Germany-Denmark-Sweden favoured pollen transport. On average, the post-seasonal LDT episodes of birch pollen to Poland occur almost every year (Poznan) or every second year (Rzeszów). The episodes are highly variable in time; thus, the pollen concentration may unexpectedly cause allergy symptoms in sensitized patients. In some cases, these episodes may be extremely severe, thereby prolonging and strengthening the exposure to birch pollen allergens.

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

Land surface phenology (LSP) provides bio-indication of ongoing climate change. It uses space-borne greenness proxies to monitor plant phenology at the landscape level from the regional to global scale. However, several unconsidered methodological and observational -related limitations may lead to misinterpretation of the satellite-derived signals. For instance, changes in species composition within a pixel could result in a change in the time series of the greenness proxy, due to the distinct phenology of the plant species involved. The change in the signal would then be misinterpreted as a phenological change while it is actually related to changes in species composition within the pixel. Other limitations include the selection of the smoothing technique and the method used to extract the LSP metrics. These not only may affect the timing of the LSP metrics but also the sign of the observed LSP change. Another and much less known limitation is related to the mixed signal from multi-canopy layers. Satellites may detect changes that corresponds to the understorey layer in complex vertical vegetation systems while the 'real' contribution of this layer (in terms of ecosystem functioning and dynamics) might be small compared to the undetected overstorey layer in cases of a late overstorey development. Here, some of the LSP basics are reviewed with emphasis on these (and other) potential sources of misinterpretation. Several aids to overcome these limitations, which include suggestions for multi methods analysis and the integration of information from satellite and ground-based sensors are provided alongside some prospective future LSP research directions.