Current and Past Projects

Mapping Chlorophyll and Carotenoid pigment dynamics across North America using MODIS MAIAC data

This goal of this project is to examine whether satellite observations of photoprotective pigments can improve monitoring of drought stress and mortality risk in western US forests. We are developing and validating a new dataset of the Chlorophyll-Carotenoid Index (CCI) from MODIS satellite data and evaluating its ability to track forest photosynthetic responses to water stress. The work investigates how CCI compares to existing vegetation indices for detecting drought-induced physiological changes and whether it can provide early warning signals for insect-induced tree mortality in evergreen conifer forests. A gridded dataset over North America (Canada, United States, Mexico) and accompanying manuscript are in preparation.

Map of North America shaded green to show the Chlorophyll-Carotenoid Index (CCI) derived from MODIS MAIAC data, with the CCI formula and a color scale from -0.3 to 0.2.
Photo of a montane conifer forest above a chart of seasonal SIF anomalies for three years against a pre-drought baseline, showing a sharp mid-summer decline in Year 2.

Impacts of wildfire and bark beetle mortality on forest carbon uptake

This work investigated whether satellite observations of solar-induced chlorophyll fluorescence (SIF), a remote sensing proxy for plant photosynthesis, can effectively detect and characterize tree mortality from wildfires and bark beetle outbreaks in western US forests. By comparing SIF responses across different disturbance types and severity levels, we found significant negative SIF anomalies up to 2-3 years prior to the onset of visible bark beetle mortality. We also found a strong agreement between the magnitude of annual SIF decline and the severity of vegetation loss due to wildfire.https://doi.org/10.1016/j.rse.2026.115550

Connecting satellite measurements of Solar-induced Chlorophyll Fluorescence (SIF) to photosynthesis in montane conifer forests

We looked at whether high-resolution satellite observations of solar-induced fluorescence (SIF) can effectively track photosynthesis patterns across complex mountain terrain in California's Sierra Nevada. We evaluated SIF against gross primary productivity (GPP) estimates from eddy covariance towers, machine-learning models, and land surface models across a 3000-meter elevation gradient. We showed that SIF can serve as a useful tool for carbon monitoring in montane evergreen conifer forests where traditional vegetation indices may not accurately capture photosynthetic activity. This work was published in Environmental Research Letters. https://doi.org/10.1088/1748-9326/ad07b4

Diagram of elevation- and season-dependent growth limits (cold and light limited, moisture limited, ideal growing conditions) above a chart comparing seasonal SIF to GPP from eddy covariance and FLUXCOM.
Map of the Salt Lake City monitoring network colored by simulated CO2 flux error, with six site locations marked, above a chart of atmospheric sensitivity versus hours before particle release.

Bayesian inverse modeling of urban CO2 emissions

We developed and tested a Bayesian inverse modeling framework to estimate urban CO2 emissions via an Observing System Simulation Experiment (OSSE) over Salt Lake City, Utah. This study explored how different ways of representing spatiotemporal uncertainties in the a priori emissions and in the observation-transport representation affect the model's ability to recover known emissions patterns. We also examined how individual monitoring sites within an urban network contribute to constraining emissions, including how nearby point sources can introduce biases into domain-wide estimates. This study was published in Elementa: Science of the Anthropocene. https://doi.org/10.1525/elementa.375

Bayesian inverse modeling tutorial & code in R

I worked with folks at NOAA's Global Monitoring Division to develop instructional materials and R scripts for running Bayesian inversions of atmospheric tracers (CO2, CH4, N2O, etc). See the code here: https://github.com/lkunik/bayesian-inversion-R

Publications