Search NASA⌕ Search

SEARCH · Search NASA

Results for “vegetation index”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Estimation of Crop Gross Primary Production (GPP): I. Impact of MODIS Observation Footprint and Impact of Vegetation BRDF Characteristics

Accurate estimation of gross primary production (GPP) is essential for carbon cycle and climate change studies. Three AmeriFlux crop sites of maize and soybean were selected for this study. Two of the sites were irrigated and the other one was rainfed. The normalized difference vegetation index (NDVI), the enhanced vegetation index (EVI), the green band chlorophyll index (CIgreen), and the green band wide dynamic range vegetation index (WDRVIgreen) were computed from the moderate resolution imaging spectroradiometer (MODIS) surface reflectance data. We examined the impacts of the MODIS observation footprint and the vegetation bidirectional reflectance distribution function (BRDF) on crop daily GPP estimation with the four spectral vegetation indices (VIs - NDVI, EVI, WDRVIgreen and CIgreen) where GPP was predicted with two linear models, with and without offset: GPP = a × VI × PAR and GPP = a × VI × PAR + b. Model performance was evaluated with coefficient of determination (R2), root mean square error (RMSE), and coefficient of variation (CV). The MODIS data were filtered into four categories and four experiments were conducted to assess the impacts. The first experiment included all observations. The second experiment only included observations with view zenith angle (VZA) ≤ 35◦ to constrain growth of the footprint size,which achieved a better grid cell match with the agricultural fields. The third experiment included only forward scatter observations with VZA ≤ 35◦. The fourth experiment included only backscatter observations with VZA ≤ 35◦. Overall, the EVI yielded the most consistently strong relationships to daily GPP under all examined conditions. The model GPP = a × VI × PAR + b had better performance than the model GPP = a × VI × PAR, and the offset was significant for most cases. Better performance was obtained for the irrigated field than its counterpart rainfed field. Comparison of experiment 2 vs. experiment 1 was used to examine the observation footprint impact whereas comparison of experiment 4 vs. experiment 3 was used to examine the BRDF impact. Changes in R2, RMSE,CV and changes in model coefficients "a" and "b" (experiment 2 vs. experiment 1; and experiment 4 vs. experiment 3) were indicators of the impacts. The second experiment produced better performance than the first experiment, increasing R2 (↑0.13) and reducing RMSE (↓0.68 g C m−2 d−1) and CV (↓9%). For each VI, the slope of GPP = a × VI × PAR in the second experiment for each crop type changed little while the slope and intercept of GPP = a × VI × PAR + b varied field by field. The CIgreen was least affected by the MODIS observation footprint in estimating crop daily GPP (R2, ↑0.08; RMSE, ↓0.42 g C m−2 d−1; and CV, ↓7%). Footprint most affected the NDVI (R2, ↑0.15; CV, ↓10%) and the EVI (RMSE, ↓0.84 g C m−2 d−1). The vegetation BRDF impact also caused variation of model performance and change of model coefficients. Significantly different slopes were obtained for forward vs. backscatter observations, especially for the CIgreen and the NDVI. Both the footprint impact and the BRDF impact varied with crop types, irrigation options, model options and VI options.

Chlorophyll↗
DOE OSTI2023

Modeling Global Vegetation Gross Primary Productivity, Transpiration and Hyperspectral Canopy Radiative Transfer Simultaneously Using a Next Generation Land Surface Model—CliMA Land

Recent progress in satellite observations has provided unprecedented opportunities to monitor vegetation activity at global scale. However, a major challenge in fully utilizing remotely sensed data to constrain land surface models (LSMs) lies in inconsistencies between simulated and observed quantities. For example, gross primary productivity (GPP) and transpiration (T) that traditional LSMs simulate are not directly measurable from space, although they can be inferred from spaceborne observations using assumptions that are inconsistent with those LSMs. In comparison, canopy reflectance and fluorescence spectra that satellites can detect are not modeled by traditional LSMs. To bridge these quantities, we presented an overview of the next generation land model developed within the Climate Modeling Alliance (CliMA), and simulated global GPP, T, and hyperspectral canopy radiative transfer (RT; 400–2,500 nm for reflectance, 640–850 nm for fluorescence) at hourly time step and 1° spatial resolution using CliMA Land. CliMA Land predicts vegetation indices and outgoing radiances, including solar-induced chlorophyll fluorescence (SIF), normalized difference vegetation index (NDVI), enhanced vegetation index (EVI), and near infrared reflectance of vegetation (NIRv) for any given sun-sensor geometry. The spatial patterns of modeled GPP, T, SIF, NDVI, EVI, and NIRv correlate significantly with existing data-driven products (mean R 2 = 0.777 for 9 products). CliMA Land would be also useful in high temporal resolution simulations, for example, providing insights into when GPP, SIF, and NIRv diverge.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Processed Surface Imagery & Raster Imagery taken onboard TigerShark (U3) UAS

Land-surface orthomosaic imagery derived from the images collected using the Micasense Altum Imager on board the UAS TigerShark (U3) on March 11th, 2021 at an altitude of 2000 ft. The Altum takes images of the land surface across five visible bands and one long wave infrared thermal band. The photogrametry software Agisoft PhotoScan v 1.4 is used to align and stitch the images into a larger composet image using the technique of structure from motion image capture to construct a dense cloud and 3-D model of the surface, which is used to produce a digital elevation model of the terraine surveyed and orthomosaic imagery. This dataset contains the raster calculated orthomosaics of various vegetative indicies, which are used to indicate plant health and land surface characteristics, including Normalized Difference Vegetation Index (NDVI), Green Normalized Difference Vegetation Index (GNDVI), Normalized Difference Water Index (NDWI), surface temperature, and Enhanced_vegetation_index (EVI). These orthomosaics were collected west of the Starkville, MS, airport, and include natural, agricultural, and industrial land coverage. The images have a resolution of 0.2 m/pix.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Temporal covariance structure of multi-spectral phenotypes and their predictive ability for end-of-season traits in maize

Vegetation indices (VIs) derived from multi-spectral imaging (MSI) platforms can be used to study properties of crop canopy, providing non-destructive phenotypes that could be used to better understand growth curves throughout the growing season. To investigate the amount of variation present in several VIs and their relationship with important end-of-season traits, genetic and residual (co)variances for VIs, grain yield and moisture were estimated using data collected from maize hybrid trials. The VIs considered were Normalized Difference Vegetation Index (NDVI), Green NDVI, Red Edge NDVI, Soil-Adjusted Vegetation Index, Enhanced Vegetation Index and simple Ratio of Near Infrared to Red (Red) reflectance. Genetic correlations of VIs with grain yield and moisture were used to fit multi-trait models for prediction of end-of-season traits and evaluated using within site/year cross-validation. To explore alternatives to fitting multiple phenotypes from MSI, random regression models with linear splines were fit using data collected in 2016 and 2017. Heritability estimates ranging from (0.10 to 0.82) were observed, indicating that there exists considerable amount of genetic variation in these VIs. Furthermore, strong genetic and residual correlations of the VIs, NDVI and NDRE, with grain yield and moisture were found. Considerable increases in prediction accuracy were observed from the multi-trait model when using NDVI and NDRE as a secondary trait. Finally, random regression with a linear spline function shows potential to be used as an alternative to mixed models to fit VIs from multiple time points.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2008

Estimating the Effect of Gypsy Moth Defloiation Using MODIS

The area of North American forests affected by gypsy moth defoliation continues to expand despite efforts to slow the spread. With the increased area of infestation, ecological, environmental and economic concerns about gypsy moth disturbance remain significant, necessitating coordinated, repeatable and comprehensive monitoring of the areas affected. In this study, our primary objective was to estimate the magnitude of defoliation using Moderate Resolution Imaging Spectroradiometer (MODIS) imagery for a gypsy moth outbreak that occurred in the US central Appalachian Mountains in 2000 and 2001. We focused on determining the appropriate spectral MODIS indices and temporal compositing method to best monitor the effects of gypsy moth defoliation. We tested MODIS-based Normalized Difference Vegetation Index (NDVI), Enhanced Vegetation Index (EVI), Normalized Difference Water Index (NDWI), and two versions of the Normalized Difference Infrared index (NDIIb6 and NDIIb7, using the channels centered on 1640 nm and 2130 nm respectively) for their capacity to map defoliation as estimated by ground observations. In addition, we evaluated three temporal resolutions: daily, 8-day and 16-day data. We validated the results through quantitative comparison to Landsat based defoliation estimates and traditional sketch maps. Our MODIS based defoliation estimates based on NDIIb6 and NDIIb7 closely matched Landsat defoliation estimates derived from field data as well as sketch maps. We conclude that daily MODIS data can be used with confidence to monitor insect defoliation on an annual time scale, at least for larger patches (greater than 0.63 km2). Eight-day and 16-day MODIS composites may be of lesser use due to the ephemeral character of disturbance by the gypsy moth.

deBeurs, K. M.↗
NASA NTRS2018

A Global Spatially Contiguous Solar-Induced Fluorescence (CSIF) Dataset Using Neural Networks

Satellite-retrieved solar-induced chlorophyll fluorescence (SIF) has shown great potential to monitor the photosynthetic activity of terrestrial ecosystems. However, several issues, including low spatial and temporal resolution of the gridded datasets and high uncertainty of the individual retrievals, limit the applications of SIF. In addition, inconsistency in measurement footprints also hinders the direct comparison between gross primary production (GPP) from eddy covariance (EC) flux towers and satellite-retrieved SIF. In this study, by training a neural network (NN) with surface reflectance from the MODerate-resolution Imaging Spectroradiometer (MODIS) and SIF from Orbiting Carbon Observatory-2 (OCO-2), we generated two global spatially contiguous SIF (CSIF) datasets at moderate spatiotemporal (0.05° 4-day) resolutions during the MODIS era, one for clear-sky conditions (2000–2017) and the other one in all-sky conditions (2000–2016). The clear-sky instantaneous CSIF (CSIF(sub clear-inst)) shows high accuracy against the clear-sky OCO-2 SIF and little bias across biome types. The all-sky daily average CSIF (CSIF(sub all-daily)) dataset exhibits strong spatial, seasonal and interannual dynamics that are consistent with daily SIF from OCO-2 and the Global Ozone Monitoring Experiment-2 (GOME-2). An increasing trend (0.39 %) of annual average CSIFall-daily is also found, confirming the greening of Earth in most regions. Since the difference between satellite-observed SIF and CSIF is mostly caused by the environmental down-regulation on SIF(sub yield), the ratio between OCO-2 SIF and CSIF(sub clear-inst) can be an effective indicator of drought stress that is more sensitive than the normalized difference vegetation index and enhanced vegetation index. By comparing CSIF(sub all-daily) with GPP estimates from 40 EC flux towers across the globe, we find a large cross-site variation (c.v. = 0.36) of the GPP–SIF relationship with the highest regression slopes for evergreen needleleaf forest. However, the cross-biome variation is relatively limited (c.v. = 0.15). These two contiguous SIF datasets and the derived GPP–SIF relationship enable a better understanding of the spatial and temporal variations of the GPP across biomes and climate.

Zhang, Yao↗
NASA NTRS2022

Yellowstone Ecological Forecasting: Assessing Change in Aspen Extent in Northern Yellowstone National Park

The removal and reintroduction of the gray wolf (Canis lupus) in Yellowstone National Park have shaped the ecological composition of this distinct landscape, representing a textbook example of trophic dynamics. With particular importance to conservation science, researchers have studied the trophic cascades between wolves and species such as elk (Cervus canadensis) and quaking aspen (Populus tremuloides). In conjunction with the National Park Service, Yellowstone National Park, Utah State University, and the University of Wisconsin–Stevens Point, this project utilized satellite remote sensing to investigate the long-term trends in aspen extent. Through random forest modeling and phenological approaches, Landsat 5 Thematic Mapper (TM; years 1986–2011) and Sentinel-2 Multispectral Instrument (MSI; years 2017–2019) datasets were used to derive color composites, Normalized Difference Vegetation Index (NDVI), Enhanced Vegetation Index (EVI), and Tasseled Cap Indices (Brightness, Greenness, Wetness). The International Space System (ISS) Global Ecosystem Dynamics Investigation (GEDI) provided canopy height data. The team consolidated results into maps and time-series which provide an in-depth depiction of aspen stand extent. The National Park Service will use these end products to assist in its management practices and inform wildlife restoration decisions within and beyond Yellowstone National Park.

Kyle Steen↗
NASA NTRS1987

Comparison of North and South American biomes from AVHRR observations

Previous analysis of the North American continent with AVHRR-derived vegetation index measurements showed a strong relation between known patterns of vegetation seasonality, productivity and the spectral vegetation index measurements. This study extends that analysis to South America to evaluate the degree to which these findings extend to tropical regions. The results show that the spectral vegetation index measurements provide a general indicator of vegetation activity across the major biomes of the Western Hemisphere of the earth, including tropical regions. The satellite-observed patterns are strongly related to the known climatology of the continents and may offer a means to improve understanding of global bioclimatology. For example, South America is shown to have a longer growing season with much earlier spring green-up than North America. The time integral of the measurements, computed from 12 composited monthly values, produces a value that is related to published net primary productivity data. However, limited net primary production data does not allow complete evaluation of satellite-observed contrasts between North and South American biomes. These results suggest that satellite-derived spectral vegetation index measurements are of great potential value in improving knowledge of the earth's biosphere.

Goward, Samuel N.↗
NASA NTRS2021

Colorado River Basin Water Resources: Pairing Phenology Citizen Science Observations with NASA Earth Observations to Monitor Riparian Vegetation

Native plant species in the Colorado River Basin are currently threatened by non-native species that are changing resource competition dynamics. Monitoring native species is essential for understanding changes in plant health throughout the year as well as how species’ distributions change inter-annually. However, current monitoring of species is limited to local efforts and cannot be efficiently conducted across more extensive regions of the basin. Project partners, including the USA National Phenology Network, United States Geological Survey, and Educating Children Outdoors, provided data and on-the-ground expertise to assist the project. NASA Earth observations, including Landsat 8 Operational Land Imager (OLI), Aqua and Terra Moderate Resolution Imaging Spectroradiometer (MODIS), and Sentinel-2 MultiSpectral Instrument (MSI) from the European Space Agency, provided remote sensing data that covers the large spatial extent of the basin. This study combined both data sources over the study period of 2013 – 2019 and explored the relationship between on-the-ground phenology observations and remotely sensed phenology signals. The team found that pairing the datasets required spatial coordinates to be accurately logged and that it was necessary to have consistent observations throughout the year by citizen scientists. The team found promising relationships between phenology observations and the two vegetation indices, normalized difference vegetation index (NDVI) and enhanced vegetation index (EVI), when study locations were in an area with a high proportion of vegetation within a satellite pixel. The results and methods outlined in this study create new opportunities to leverage remote sensing alongside citizen science data for widespread monitoring of non-native species.

Tyler Pantle↗
NASA NTRS1985

The relationship of global green leaf biomass to atmospheric CO2 concentrations

Advanced very high resolution radiometer data from NOAA's polar orbiting meteorological satellite have been obtained globally for a 21 month period, processed to produce a green leaf biomass spectral vegetative index for the entire terrestrial surface by month, zonally aggregated by latitude, and compared to atmospheric CO2 concentrations from observing stations. A strong inverse association was found between the monthly Pt. Barrow CO2 concentrations and the vegetation index measurements from 50 deg N to 80 deg N, between the monthly Mauna Loa CO2 concentrations and the vegetation index measurements from 10 deg N to 30 deg N, 10 deg N to 80 deg N, and the global total, and between the globally averaged CO2 concentrations and the globally averaged vegetation index. No relationships between atmospheric CO2 concentrations and the vegetative index measurements from any latitude zone or combinations of zones were found for the South Pole station.

Tucker, C. J.↗
DOE OSTI2021

Satellite solar-induced chlorophyll fluorescence and near-infrared reflectance capture complementary aspects of dryland vegetation productivity dynamics

Mounting evidence indicates dryland ecosystems play an important role in driving the interannual variability and trend of the terrestrial carbon sink. Nevertheless, our understanding of the seasonal dynamics of dryland ecosystem carbon uptake through photosynthesis [gross primary productivity (GPP)] remains relatively limited due in part to the limited availability of long-term data and unique challenges associated with satellite remote sensing across dryland ecosystems. Here, we comprehensively evaluated longstanding and emerging satellite vegetation proxies in their ability to capture seasonal dryland GPP dynamics. Specifically, we evaluated: 1) reflectance-based proxies normalized difference vegetation index (NDVI), soil adjusted vegetation index (SAVI), near infrared reflectance index (NIR v ), and kernel NDVI (kNDVI) from the MODerate resolution Imaging Spectroradiometer (MODIS); and 2) newly available physiologically-based proxy solar-induced chlorophyll fluorescence (SIF) from the TROPOspheric Monitoring Instrument (TROPOMI). As a performance benchmark, we used GPP estimates from a robust network of 21 western United States eddy covariance tower sites that span representative gradients in dryland ecosystem climate and functional composition. We found that NIR v and SIF were the best performing GPP proxies and captured complementary aspects of seasonal GPP dynamics across dryland ecosystem types. NIR v offered better performance than the other proxies across relatively low-productivity, sparsely non-evergreen vegetated sites (R 2 = 0.59 ± 0.13); whereas SIF best captured seasonal dynamics across relatively high-productivity sites, including evergreen-dominated sites (R 2 = 0.74 ± 0.07). Notably, across grass-dominated sites, all reflectance-based proxies (NDVI, SAVI, NIRv and kNDVI) showed significant seasonal bias (hysteresis) that strengthened with the total fraction of woody vegetation cover, likely due to seasonal patterns in woody vegetation reflectance that are unrelated to or decoupled from GPP. In conclusion, future efforts to fully integrate the complementary strengths of NIR v and SIF could significantly improve our understanding and representation of dryland GPP dynamics in satellite-based models.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2014

Scientific Impact of MODIS C5 Calibration Degradation and C6+ Improvements

The Collection 6 (C6) MODIS (Moderate Resolution Imaging Spectroradiometer) land and atmosphere data sets are scheduled for release in 2014. C6 contains significant revisions of the calibration approach to account for sensor aging. This analysis documents the presence of systematic temporal trends in the visible and near-infrared (500 m) bands of the Collection 5 (C5) MODIS Terra and, to lesser extent, in MODIS Aqua geophysical data sets. Sensor degradation is largest in the blue band (B3) of the MODIS sensor on Terra and decreases with wavelength. Calibration degradation causes negative global trends in multiple MODIS C5 products including the dark target algorithm's aerosol optical depth over land and Ångstrom exponent over the ocean, global liquid water and ice cloud optical thickness, as well as surface reflectance and vegetation indices, including the normalized difference vegetation index (NDVI) and enhanced vegetation index (EVI). As the C5 production will be maintained for another year in parallel with C6, one objective of this paper is to raise awareness of the calibration-related trends for the broad MODIS user community. The new C6 calibration approach removes major calibrations trends in the Level 1B (L1B) data. This paper also introduces an enhanced C6C calibration of the MODIS data set which includes an additional polarization correction (PC) to compensate for the increased polarization sensitivity of MODIS Terra since about 2007, as well as detrending and Terra- Aqua cross-calibration over quasi-stable desert calibration sites. The PC algorithm, developed by the MODIS ocean biology processing group (OBPG), removes residual scan angle, mirror side and seasonal biases from aerosol and surface reflectance (SR) records along with spectral distortions of SR. Using the multiangle implementation of atmospheric correction (MAIAC) algorithm over deserts, we have also developed a detrending and cross-calibration method which removes residual decadal trends on the order of several tenths of 1% of the top-of-atmosphere (TOA) reflectance in the visible and near-infrared MODIS bands B1-B4, and provides a good consistency between the two MODIS sensors. MAIAC analysis over the southern USA shows that the C6C approach removed an additional negative decadal trend of Terra (Delta)NDVI approx.0.01 as compared to Aqua data. This change is particularly important for analysis of vegetation dynamics and trends in the tropics, e.g., Amazon rainforest, where the morning orbit of Terra provides considerably more cloud-free observations compared to the afternoon Aqua measurements.

Terra and Aqua↗
NASA NTRS1977

Use of near infrared/red radiance ratios for estimating vegetation biomass and physiological status

The application of photographic infrared/red (ir/red) reflectance or radiance ratios for the estimation of vegetation biomass and physiological status were investigated by analyzing in situ spectral reflectance data from experimental grass plots. Canopy biological samples were taken for total wet biomass, total dry biomass, leaf water content, dry green biomass, dry brown biomass, and total chlorophyll content at each sampling date. Integrated red and photographic infrared radiances were regressed against the various canopy or plot variables to determine the relative significance between the red, photographic infrared, and the ir/red ratio and the canopy variables. The ir/red ratio is sensitive to the photosynthetically active or green biomass, the rate of primary production, and actually measures the interaction between the green biomass and the rate of primary production within a given species type. The ir/red ratio resulted in improved regression significance over the red or the ir/radiances taken separately. Only slight differences were found between ir/red ratio, the ir-red difference, the vegetation index, and the transformed vegetation index. The asymptotic spectral radiance properties of the ir, red, ir/red ratio, and the various transformations were evaluated.

Tucker, C. J.↗
NASA NTRS1978

Red and photographic infrared linear combinations for monitoring vegetation

In situ collected spectrometer data were used to evaluate and quantify the relationships between various linear combinations of red and photographic infrared radiances and experimental plot biomass, leaf water content, and chlorophyll content. The radiance variables evaluated included the red and photographic infrared (IR) radiance and the linear combinations of the IR/red ratio, the square root of th IR/red difference, the vegetation index, and the transformed vegetation index. In addition, the corresponding green and red linear combinations were evaluated for comparative purposes. Three data sets were used from June, September, and October sampling periods. Regression analysis showed the increase utility of the IR and red linear combinations vis-a-vis the same green and red linear combinations. The red and IR linear combinations had 7% and 14% greater regression significance than the green and red linear combinations for the June and September sampling periods, respectively. The VI, TVI, and square root of the IR/red ration were the most significant followed closely by the IR/red ratio. Less than 6% difference separated the highest and lowest of these four IR and red linear combinations. The use of these linear combinations was shown to be sensitive primarily to the green leaf area or green leaf biomass.

Tucker, C. J.↗
NASA NTRS1985

Analysis of Terrestrial Conditions and Dynamics

An ecological model is developed to estimate annual net primary productivity of vegetation in twelve major North American biomes. Three models are adapted and combined, each addressing a different factor known to govern primary productivity, i.e., photosynthesis, respiration, and moisture availability. Measures of intercepted photosynthetically active radiation (1PAR) for input to the photosynthesis model are derived from spectral vegetation index data. Normalized Difference Vegetation Index (NDVI) data are produced from NOAA-7 Advanced Very High Resolution Radiometer (AVHRR) observations for April 1982 through March 1983. NDVI values are sampled from within the biomes at locations for which climatological data are available. Monthly estimates of Net Primary Productivity (NPP) for each sample location are generated and summed over the twelve month period. These monthly estimates are averaged to produce a single annual estimated NPP value for each biomes. Comparison of estimated NPP values with figures reported in the literature produces a correlation coefficient of 85.

Goward, S. N.↗
NASA NTRS1979

Red and photographic infrared linear combinations for monitoring vegetation

The relationships between various linear combinations of red and photographic infrared radiances and vegetation parameters are investigated. In situ spectrometers are used to measure the relationships between linear combinations of red and IR radiances, their ratios and square roots, and biomass, leaf water content and chlorophyll content of a grass canopy in June, September and October. Regression analysis shows red-IR combinations to be more significant than green-red combinations. The IR/red ratio, the square root of the IR/red ratio, the vegetation index (IR-red difference divided by their sum) and the transformed vegetation index (the square root of the vegetation index + 0.5) are found to be sensitive to the amount of photosynthetically active vegetation. The accumulation of dead vegetation over the year is found to have a linearizing effect on the various vegetation measures.

Tucker, C. J.↗
DOE OSTI2026

Data for Yield from Iowa’s first commercial miscanthus fields: implications of spatial variability for productivity and sustainability beyond research plots

This dataset contains biomass yield measurements and associated vegetation index data collected from commercial Miscanthus × giganteus fields in eastern Iowa during the 2022–2023 growing seasons. The data support the analyses presented in the article: “Yield From Iowa's First Commercial Miscanthus Fields: Implications of Spatial Variability for Productivity and Sustainability Beyond Research Plots.” We collected 105 ground-truth biomass samples from four mature commercial fields (>4 years old) covering 92.81 ha. Samples were taken from 3 m² quadrats that were hand-harvested in alignment with commercial harvest timing. Stem biomass (excluding leaves) was weighed, moisture-corrected, and converted to dry-matter yield expressed in Mg DM ha⁻¹. Sampling locations were selected to capture spatial variability visible in aerial imagery and were recorded using RTK GPS. Each biomass observation was paired with vegetation indices derived from high-resolution PlanetScope satellite imagery (3 m resolution). Images were acquired throughout the growing season, and indices were calculated to evaluate their ability to predict end-of-season biomass yield. Statistical and machine learning approaches were used to identify key predictors, and a linear regression model based on end-of-July Green Normalized Difference Vegetation Index (GNDVI) was developed and evaluated. This repository includes the data used in that modeling workflow. Management practices, economic data, full imagery time series, and additional methodological details are described in the associated publication and are not included here. The dataset consists of three comma-separated value (CSV) files: 1. Combine_Groundtruth_Yield_VI_22_23.csv This file contains ground-truth biomass yield measurements and associated key vegetation index values collected during the 2022 and 2023 growing seasons. Rows: 105 observations Columns: Year — Year of observation (2022 or 2023) Field — Field location identifier Sample_number — Unique sample identifier GNDVI_End_Jul — Green Normalized Difference Vegetation Index calculated at end of July GNDVI_End_Aug — Green Normalized Difference Vegetation Index calculated at end of August NDRE_End_Aug — Normalized Difference Red Edge index calculated at end of August Biomass_Stem_Yield_MgDM/ha — Measured stem biomass yield (megagrams dry matter per hectare) 2. trainData_GNDVI.csv This file contains the subset of observations used to train the predictive relationship between July GNDVI and biomass yield. Rows: 76 observations Columns: Unnamed: 0 — Row index retained from the original data processing workflow GNDVI_End_Jul — GNDVI at end of July Stem_Yield_MgDM/ha — Observed stem biomass yield (Mg DM ha⁻¹) 3. testData_GNDVI.csv This file contains the test dataset used to evaluate model performance. Rows: 29 observations Columns: Unnamed: 0 — Row index retained from the original data processing workflow GNDVI_End_Jul — GNDVI at end of July Predicted_Yield_MgDM/ha — Model-predicted stem biomass yield (Mg DM ha⁻¹) Observed_Yield_MgDM/ha — Measured stem biomass yield (Mg DM ha⁻¹)

Potential yield, yield gap, in-field management, y↗