Research Visualizations

Visualizations from dissertation and collaborative research on rice paddy carbon flux, phenology, and remote sensing in Arkansas. Combines eddy covariance ground truth, Landsat/Sentinel/MODIS satellite imagery, machine learning, and process-based modeling to estimate GPP, planting and harvest dates, and CH₄ flux at field-to-state scale.

Methods represented:
Remote Sensing Random Forest / XGBoost Harmonic Regression GAM / Non-linear Regression Time Series Analysis Eddy Covariance Vegetation Photosynthesis Model (VPM) GIS & Spatial Analysis

Methodology Process Flow Diagram

Project Workflow & Predictive Modeling Framework — Spatial GPP, Growing Season Length, and CH₄

Four data sources — eddy covariance ground truth, gridded environmental data, satellite reflectance, and rice land cover — integrated across two phases: Phase 1 site-scale CH₄/GPP estimation and validation; Phase 2 state-scale predictive modeling using Sentinel, Landsat, and MODIS vegetation indices alongside soil and agronomic datasets.

Remote Sensing Eddy Covariance Machine Learning GIS
Conceptual Framework Site to State Scale

Conceptual Framework for Ecosystem Flux Scaling

Upscaling from site-level, high-frequency Eddy Covariance measurements of CH₄ and NEE to regional, state-scale GPP and growing season length maps. Predictors span climatological, plant status (EVI, NDVI), water status (LSWI), and agronomic variables.

Eddy Covariance Spatial Scaling
Phenological Stages

Vegetation Phenological Stages

Maps key phenological terminology — SOS, GUD, POS, EOS — onto the plant growth curve from emergence through senescence. Applies to single cropping systems.

Phenology

Planting Harvest Date Violin Boxplot

Planting & Harvest Date Distribution

Violin-boxplot composite showing the spread, median, and IQR of planting dates (mean DOY 120, N=978) and harvest dates (mean DOY 257, N=673) across Arkansas rice fields. Derived from USDA survey records.

Phenology Descriptive Statistics
kNDVI max to PD HD Histograms

Lag from kNDVI Peak to Planting & Harvest Dates

Histograms of the lag (days) between the kNDVI maximum and observed planting date (mean 88.6 d, SD 16.7) and harvest date (mean 51.5 d, SD 14.4), demonstrating the predictive signal embedded in satellite vegetation indices.

Remote Sensing Phenology Statistics
Lag vs Growing Season Length

Lag between SOS and Planting Date vs. Growing Season Length

Three-panel scatterplot comparing growing season length (HD–PD) against the lag between planting date and three SOS definitions (derivative, threshold, UD methods). Pearson r = 0.43–0.60, indicating that earlier-detected green-up relative to planting is associated with longer growing seasons.

Phenology Regression Remote Sensing
RMSE vs Climate Anomaly

Climate Drivers of Planting & Harvest Date Model Error

Six-panel figure (leave-2-years-out cross-validation) showing how interannual climate anomalies (min air temp, min soil temp, VPD) relate to model RMSE for PD (panel A, r = –0.14 to –0.42) and HD (panel B, r = –0.56 to –0.72, p < 0.05). Warmer and drier years yield more accurate harvest date predictions.

Machine Learning Cross-validation Phenology

VPM Yield Area

Relationship between Yield and GPP

Spatial relationship between rice yield and VPM-modeled GPP across Arkansas production areas.

VPMSpatial
Temporal GPP VPM

Time Series of GPP — Three Methods

Temporal comparison of GPP estimates derived from three independent approaches, showing consistency and divergence across the growing season.

VPMEddy Covariance
Interannual Yield GPP

Interannual Trends of Yield & GPP

Multi-year trajectory of rice yield alongside GPP, highlighting co-variation patterns and anomalous seasons.

Time Series
Region Boxplot

GPP Comparison — Rice Production Regions of Arkansas

Distribution of GPP across distinct rice-growing regions, revealing spatial variability in productivity potential.

RegionalStatistics
GPP Variety Plot

Factors Affecting GPP by Cultivar

How cultivar variety and biophysical factors modulate GPP, compared against eddy covariance ground truth.

Eddy Covariance
GPP Multiple Sources

Temporal Dynamics of GPP — Multiple Models

Side-by-side temporal profiles of GPP predicted by different model sources, enabling cross-model comparison across the growing season.

Multi-modelTime Series
Rice LSWI

LSWI Influence on GPP across Clay Content Levels

Three-panel figure: (1) histogram of clay content in Arkansas rice fields, (2) LSWI–GPP relationship across clay strata, and (3) growing frequency vs. LSWI residuals across clay content levels.

LSWIRegressionSoil

LUE Biophysical

Light Use Efficiency vs. Eight Biophysical Factors

Eight-panel scatterplot matrix (colored by days after planting) showing LUE responses to VPD (r = –0.32), relative humidity (r = 0.54), DBSI (r = –0.73), AWEInsh (r = –0.65), IAVI (r = 0.82), air temperature (r = 0.54), evapotranspiration (r = 0.42), and cumulative GDD (r = 0.67). GAM smooths overlaid. IAVI and VPD emerge as dominant controls.

GAM Regression Vegetation Indices LUE
fAPAR distributions

fAPAR Distributions — EVI, NDVI, and LAI Parameterizations

Three-panel histogram comparing fAPAR computed from EVI (panel A), NDVI (panel B), and LAI (panel C). Highlights the proportion of observations where fAPAR > 1 (gold, physically implausible) and LUE > 1 (pink), informing parameterization choices for the VPM.

fAPAR VPM Distribution
fAPAR vs LAI

fAPAR–LAI Relationship with GPP and Phenological Timing

Scatter plot of fAPAR against Leaf Area Index (m²/m²) fitted with the Beer-Lambert exponential model (fAPAR = 1 – e–K×LAI). Points colored by days after planting (yellow→teal) and sized by EC-measured GPP, showing canopy light interception dynamics across the growing season.

fAPAR · LAI GPP EC Non-linear Fit
LUE over Time

Light Use Efficiency Across Time

Seasonal and interannual pattern of LUE — a key driver of GPP in the VPM framework — relative to phenological stages.

LUETime Series

GPP VI Dual Axis

Relationship between GPP and Vegetation Indices

Dual-axis visualization linking GPP dynamics to satellite-derived vegetation indices (NDVI, EVI, kNDVI), with regression equations for each relationship.

NDVI · EVI · kNDVIRegression
Harmonic Fit Field 30

Harmonic Fit of kNDVI — Field 30

Harmonic regression applied to the kNDVI time series of a single rice field, extracting phenological transition dates from a smooth seasonal curve.

kNDVIHarmonic Regression

Combined violin lag hist

Planting & Harvest Date — Violin, Lag, and Histogram

Multi-panel composite integrating violin plots, lag distributions, and histograms to characterize the full temporal spread of planting and harvest dates across the Arkansas rice dataset.

PhenologyDistribution
Correlation Plot

Pairwise Correlation Plot

Triangular correlation matrix among planting/harvest dates, meteorological indicators, and vegetation indices. Strong positive correlations (yellow) between SOS, EOS, and thermal variables (GDD, soil temperature); soil organic carbon shows limited co-variation. Built with ggplot2 in R.

CorrelationPhenology