The research, led by Piyapon Kankong of the University of Tokyo with colleagues including Toshiaki Owari, Takuya Hiroshima, and Nyo Me Htun, was conducted at the University of Tokyo Hokkaido Forest (UTHF) in Furano, Japan. This site, established in 1899, hosts an arboretum containing roughly 270 tree species and has maintained meticulous phenology records since 1930, including budburst, flowering, leaf coloration, and natural defoliation. That long observational tradition made it an ideal proving ground for testing whether satellite data could match the accuracy of trained observers standing beneath the canopy with binoculars.

The satellite at the center of the study is PlanetScope, a constellation operated by Planet Labs that delivers images at 3-meter spatial resolution with near-daily revisit times. That combination is critical. Workhorse environmental satellites such as MODIS (500 m), GCOM-C (250 m), Landsat (30 m), and even Sentinel-2 (10 m) simply cannot separate individual tree crowns, which are typically smaller than 10 meters across. PlanetScope’s 3-meter pixels, by contrast, can isolate the crowns of larger trees, and its daily coverage means researchers rarely miss the rapid transitions that mark autumn senescence. Over the three study years from 2022 to 2024, the team assembled 67, 73, and 103 cloud-free images of the site, respectively.

To validate the satellite data, the researchers selected 49 canopy trees belonging to 31 deciduous species, each with a fully exposed crown visible from above. Ground staff recorded leaf coloration as a percentage of the canopy, from 0 to 100, once or twice per week from late August through early December. These observations were classified into stages that defined three key metrics: the start of autumn phenology (SOA, roughly 20 percent coloration), the middle (MOA, about 50 percent), and the end (EOA, around 80 percent). The team also flew a DJI Mavic 3 Multispectral drone at 80 meters above the canopy to create centimeter-scale orthomosaics used to delineate each tree crown precisely; crown diameters in the sample ranged from 3.75 meters in Styrax obassia to 15.57 meters in Tilia japonica.

With crowns mapped, the researchers extracted six vegetation indices from the PlanetScope time series: NDVI, kNDVI, EVI, EVI2, GRVI, and the newer near-infrared reflectance of vegetation (NIRv). Time series were gap-filled, smoothed with a Savitzky-Golay filter, and fitted with double logistic models using the R package phenofit, with the Beck model providing the best fit. Two extraction methods were then compared: a threshold-based approach defining SOA, MOA, and EOA at 80, 50, and 20 percent of the seasonal amplitude, and a third-order derivative method that locates phenological dates from curvature changes in the fitted curves. The results were evaluated against ground observations using correlation, RMSE, and bias statistics.

The spectral story was clear. Among the raw surface reflectance bands, near-infrared performed best, correlating at r = -0.767 with the progression of leaf coloration across 2,257 observations. That makes physical sense: as leaves senesce, their internal cellular structure breaks down, reducing the scattering of near-infrared light, while chlorophyll degradation unmasks carotenoids and anthocyanins that raise reflectance in visible bands. Vegetation indices that incorporate the NIR band outperformed any single band, with EVI showing the strongest relationship (r = -0.791), followed closely by EVI2 and NIRv. GRVI, which excludes the near-infrared, tracked coloration least effectively.

When it came to extracting specific phenological dates, the threshold-based method generally outperformed the derivative approach, producing lower and more stable errors across species. The middle of autumn phenology emerged as the most reliably captured stage. For the ten best-performing species, which included Larix kaempferi, Carpinus cordata, Ginkgo biloba, and Quercus crispula, NIRv achieved the highest precision for MOA (r = 0.774 to 0.780, RMSE of 9.2 to 10.8 days), while EVI2 delivered the best overall accuracy with an RMSE of just 4.8 days. The authors describe this as a trade-off between precision and accuracy: NIRv is the most consistent index over time, but EVI2 lands closest to the true dates.

Species identity mattered enormously. Species with gradual, visually distinct coloration, such as larch, hornbeam, and maple, showed strong agreement between satellite and ground metrics, while Populus suaveolens, Prunus sargentii, and Tilia japonica, which color weakly, defoliate rapidly, or suffer understory interference, performed poorly. Average RMSE across all species ranged from 19.3 days in Larix kaempferi to 36.8 days in Populus suaveolens. Drone imagery helped explain the discrepancies: high-agreement species displayed clear intra-crown color gradients that PlanetScope pixels captured faithfully, whereas low-agreement species showed weak coloration or rapid leaf loss that blurred the satellite signal. Systematic biases also appeared, with satellite-derived start dates tending to run early, likely because indices like NIRv detect declines in photosynthesis before color becomes visible, and end dates tending to run late, likely because green understory vegetation contaminates the spectral signal after the canopy has finished turning.

Why does tree-by-tree autumn monitoring matter? Individual trees are the basic units of forest communities, and their phenological timing shapes photosynthetic seasonality, hydrological regulation, and nutrient cycling. Studies have documented large phenological variation both within and among species, and autumn phenology in particular responds to warming in complex, poorly understood ways, unlike the more predictable advance of spring budburst. Because the growing season length, set by the interval between spring and autumn events, is tightly linked to annual carbon dioxide sequestration, uncertainty about autumn timing translates directly into uncertainty about the forest carbon sink. Fine-scale satellite monitoring could help close that gap.

The study is not without limitations, which the authors acknowledge candidly. Each 3-meter pixel can mix signals from understory plants, cloud cover reduces image availability, and radiometric inconsistencies among the many PlanetScope sensors introduce noise that demands careful filtering. The 3-meter resolution also restricts monitoring to crowns larger than 3 meters, and the analysis covered a single site over three years with deciduous species only. Future work should compare PlanetScope against drone imagery, explore the constellation’s underused red-edge band, which is sensitive to pigment change, and extend validation across multiple forest types and longer periods. Still, the demonstration that a commercial small-satellite constellation can time the middle of autumn coloration to within about five days for well-behaved species marks a genuine advance. As climate change continues to shuffle the calendar of leaf fall, scientists now have a tool that can watch it happen, one tree at a time.

Subject of Research: Satellite monitoring of autumn leaf phenology at the individual tree scale

Article Title: Monitoring autumn leaf phenology at the individual tree scale using PlanetScope satellite imagery and ground-based observations

Article References: Kankong, P., Owari, T., Hiroshima, T., & Htun, N. M. (2026). Monitoring autumn leaf phenology at the individual tree scale using PlanetScope satellite imagery and ground-based observations. Discover Ecology, 2(1), Article 1. https://doi.org/10.1007/s44396-025-00018-5

Image Credits: AI Generated

DOI: 10.1007/s44396-025-00018-5

Keywords: autumn phenology, PlanetScope, remote sensing, leaf coloration, vegetation indices, NIRv, EVI2, individual tree crowns, UAV imagery, forest carbon, climate change, Hokkaido

The core insight of the research is deceptively simple but methodologically powerful. An annual count of compound heatwave-ozone events can rise for two very different reasons: either heatwaves as a whole become more frequent, or a larger share of the heatwaves that do occur coincide with elevated ozone. These two scenarios expose vegetation to fundamentally different conditions, yet both inflate the same summary statistic. To separate them, the team constructed a statistical framework in which the total number of heatwave events was held fixed, so that the coefficient associated with compound events isolated the effect of ozone co-occurrence per se rather than general heatwave activity. This fixed-total parameterization, compared against a conventional separate-count model, revealed how much hidden information the raw compound-event count had been carrying.

The data underpinning the analysis were assembled from some of the most advanced environmental datasets available. Daily maximum temperatures came from the ERA5 atmospheric reanalysis produced by the European Centre for Medium-Range Weather Forecasts, with heatwave thresholds calibrated against the 1983 to 2012 baseline at each pixel using the 90th percentile of warm-season temperatures. A heatwave was declared when daily maxima exceeded this local threshold for at least three consecutive days. Ozone exposure was mapped using the CHAP dataset, a machine-learning reconstruction of daily ground-level ozone at one-kilometer resolution that blends ground stations, satellite retrievals, and reanalysis fields, with a cross-validation coefficient of determination of 0.87. Days were classified as ozone-stressed when the eight-hour daily maximum concentration exceeded 100 micrograms per cubic meter, a widely used operational cutoff in Chinese compound-exposure studies.

For the biological response, the researchers turned to the MODIS MCD12Q2 land-surface phenology product, which tracks the timing of seasonal canopy transitions from vegetation-index time series. Their focal variable was the end of season, or EOS, the date at which the canopy moves toward dormancy in autumn. Croplands were deliberately excluded because planting, irrigation, and harvest schedules dominate their seasonal rhythms and would obscure environmental signals. That left forests and grasslands, the two dominant natural vegetation types covering roughly a quarter and nearly thirty percent of China’s landscape respectively. Using a fixed-effects panel regression over more than 48,000 sampled pixels, the team quantified how year-to-year variation in exposure was associated with year-to-year variation in EOS within each pixel, while absorbing persistent local differences and shared annual fluctuations. Uncertainty was estimated with a spatial block bootstrap of 2,000 replicates to guard against the statistical distortion caused by neighboring pixels sharing conditions.

The mapping of compound stress itself produced striking geography. Compound heatwave-ozone events clustered most densely over the North China Plain and the Fenwei Plain, regions where dense industrial activity and vehicle fleets supply abundant nitrogen oxides and volatile organic compounds, the raw ingredients of ozone, and where persistent anticyclonic systems trap both heat and pollutants over the boundary layer. Roughly one in five land pixels fell into the upper compound-frequency classes. Heatwave-only events, by contrast, were concentrated in Southern China and the Tarim Basin, where monsoon cloudiness and stronger atmospheric ventilation suppress photochemical ozone accumulation despite frequent heat. Ozone-only stress days blanketed the mid-latitudes, with more than thirty percent of the land area experiencing over 114 such days per year. Temporally, the study identified what the authors call a misalignment-overlap pattern: ozone pollution peaks in June, heatwaves intensify in July and August, and so compound exposure compresses into a narrow midsummer window that coincides exactly with the peak growing season.

When the phenological associations were computed, forests and grasslands parted company in an unexpected way. In forests, holding the total number of heatwave events constant, each additional ozone-coincident heatwave was associated with an earlier end of season, by about 0.293 days per event, with a confidence interval comfortably excluding zero. Grasslands showed the opposite tendency, a positive coefficient of about 0.114 days per event, though it was statistically uncertain. The formal grassland-minus-forest contrast was 0.407 days per event, a clear and significant ecosystem difference. Perhaps most revealing was what happened when the conventional parameterization was compared with the fixed-total one: in grasslands, the positive compound-event association largely evaporated once overall heatwave frequency was accounted for, indicating that most of the apparent grassland signal had been borrowed from the general heatwave regime. In forests, by contrast, the negative association actually strengthened, suggesting that ozone co-occurrence carries genuine information about autumn canopy decline beyond the sheer number of heatwaves.

The physiological logic behind this forest response is plausible, though not yet fully demonstrated. Ozone enters leaves almost exclusively through stomata, the microscopic pores that plants open to absorb carbon dioxide, and once inside it triggers oxidative damage that accelerates senescence. Stomatal behavior, however, is not a passive conduit. It responds strongly to temperature, vapor pressure deficit, and the plant’s water status, and hot, dry conditions can partially close stomata, limiting ozone uptake even as ambient concentrations climb. Forests and grasslands differ markedly in rooting depth, canopy structure, and water-use regulation, which is one reason the same atmospheric exposure might translate into different canopy outcomes. The authors caution that resolving the mechanism will require direct measurements of stomatal conductance, plant water status, and stomatal ozone flux, since event-based counts of threshold exceedances compress away the intensity and duration of exposure that determine the physiological dose plants actually receive.

One of the study’s most vivid findings concerns elevation. In forests, the negative association between ozone co-occurrence and end-of-season timing was strongest at low elevations and progressively weakened upward, at a rate of roughly 0.242 days per event per thousand meters. Across the forest elevation range sampled, from about 155 to 3,258 meters, the conditional coefficient swung from minus 0.518 to plus 0.234 days per event, although the high-elevation estimate carried wide uncertainty. Topography reshapes surface energy balance, water redistribution, and evaporative demand, and previous satellite analyses have shown that elevational patterns of autumn phenology cannot be explained by temperature alone. The pattern suggests that forests are most phenologically sensitive to ozone-coincident heatwaves when those events occur under warmer, more water-demanding lowland conditions, a finding consistent with evidence from the 2022 European heat and drought extremes, where forest carbon uptake and autumn canopy activity depended heavily on the seasonal hydrothermal background.

The study’s authors are candid about its limits. The 2014 to 2023 window is short for characterizing rare compound extremes, and annual event counts inevitably compress differences in event timing, duration, intensity, and antecedent conditions, all of which are known to modulate vegetation responses. Satellite-derived end-of-season dates, moreover, capture canopy greenness rather than the exact cessation of photosynthesis, and different remote-sensing proxies can disagree about the timetable of senescence. Integrating solar-induced chlorophyll fluorescence, eddy-covariance carbon flux measurements, and flux-based ozone metrics such as phytotoxic ozone dose would provide a more physiologically grounded picture of how compound exposure propagates from leaf chemistry to landscape-scale phenology.

Even with those caveats, the implications are substantial. Climate projections indicate that heatwaves will grow more frequent and severe, and that ozone episodes will intensify in many rapidly developing regions, making compound exposure an increasingly routine feature of the growing season rather than a statistical curiosity. The finding that the ecological meaning of a compound-event count depends on whether it reflects more heatwaves or more ozone within those heatwaves offers modelers a sharper instrument: vegetation projections that treat compound events as a single undifferentiated category may misattribute the drivers of phenological change. For China’s forests, particularly the densely exposed lowland forests of the north, the message is that air-quality management is climate-adaptation policy. curbing ozone precursors during heatwave episodes could ease an invisible but measurable pressure that pulls the growing season to an early close, with consequences for carbon storage, autumn land-surface energy exchange, and the ecosystems that depend on the length of the green season.

Subject of Research: Associations between compound heatwave-ozone exposure and autumn vegetation phenology across Chinese forests and grasslands

Article Title: Distinguishing heatwave occurrence from ozone co-occurrence reveals contrasting autumn phenological associations across Chinese ecosystems

Article References: Hou, W., Hu, Z., Wang, L., Wang, M., & Liu, X. (2026). Distinguishing heatwave occurrence from ozone co-occurrence reveals contrasting autumn phenological associations across Chinese ecosystems. Environmental and Sustainability Indicators, 32, Article 101528. https://doi.org/10.1016/j.indic.2026.101528

Image Credits: AI Generated

DOI: 10.1016/j.indic.2026.101528

Keywords: heatwaves, ground-level ozone, autumn phenology, end of season, forests, grasslands, China, compound events, satellite remote sensing, fixed-effects regression, elevation gradient, stomatal ozone uptake