SEARCH · Engineering Papers

Results for “MODIS”

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.

SourceResearch connectionFrom yearThrough year

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

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Data and Scripts associated with a manuscript on ecosystem responses to wildfires in the Columbia River Basin

This data package is associated with the publication “Ecosystem leaf area, gross primary production, and evapotranspiration responses to wildfire in the Columbia River Basin” submitted to Biogeosciences (Shi et al., 2024; doi: 10.22541/au.171053013.30286044/v1). In this research, data products, leaf area index (LAI), gross primary production (GPP), and evapotranspiration (ET), from the Moderate Resolution Imaging Spectroradiometer (MODIS) are used to quantify the resistance and resilience of different ecosystem types in the Columbia River Basin (CRB). A machine learning algorithm, random forest (RF), was used to examine the impacts of precipitation, vapor pressure deficit (VPD), and burn severity from Monitoring Trends in Burn Severity (MTBS) on ecosystem resilience. The data package includes the processed MODIS data products, precipitation, VPD, and burn severity in 138 fire regions in CRB and the input files for RF model training. This data package includes six folders. The MODIS products are included in three MODIS_* folders with shell scripts for data clipping and *ncl files for data processing: (1) “/MODIS_LAI_CRB”; (2) “/MODIS_GPP_CRB”; and (3) “/MODIS_ET_CRB”. All the processed data for each fire event are NetCDF formatted. The MTBS burn severity data and the shell and *ncl scripts used for data processing are in the folder named (4) “MTBS_fire”. The ERA meteorological fields and the data processing scritps are in (5) “ERA_Var_CR”. All the scripts for figure development are in the format of *ncl and in the folder (6) “paper_scripts”. See the file ending in “flmd.csv” for a list of all files contained in this data package and descriptions for each. Tabular column headers and units are described in the data dictionary file ending in “dd.csv”.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

SLOPE daily and 250 m gross primary productivity (GPP) for the CONUS, 2000-2019, Carbon Monitoring System (CMS)

The SatelLite Only Photosynthesis Estimation (SLOPE) GPP product is a daily, 250 m GPP dataset covering the Contiguous United States (CONUS) from 2000 to present with 1 day latency. Gross primary productivity (GPP) quantifies the amount of carbon dioxide (CO2) fixed by plants through photosynthesis. Although as a key quantity of terrestrial ecosystems, there is a lack of high-spatial-and-temporal-resolution, real-time, and observation-based GPP products. This product has been developed to address this critical gap, leveraging a number of MODIS land and atmosphere products. There are three distinct features of the SLOPE GPP production algorithm. (1) SLOPE couples machine learning models with MODIS atmosphere and land products to accurately estimate PAR. (2) SLOPE couples highly efficient and pragmatic gap-filling and filtering algorithms with surface reflectance acquired by both Terra and Aqua MODIS satellites to derive a soil-adjusted NIRv (SANIRv) dataset. (3) SLOPE couples a temporal pattern recognition approach with a long-term Crop Data Layer (CDL) product to predict dynamic C4 crop fraction. PAR, SANIRv and C4 fraction are used to drive a parsimonious model with only two slope parameters to estimate GPP along with a quantitative uncertainty on a per-pixel and daily basis. The slope GPP product has a R2 = 0.84 and a root-mean-square error (RMSE) of 1.65 gC m-2 d-1, evaluated against from 50 AmeriFlux eddy covariance sites (332 site-years). Archived data from 2000 to 2019 are publicly available at the NASA's Oak Ridge National Laboratory Distributed Active Archive Center (ORNL DAAC). Data from 2020 are available from the authors upon request. All data are projected in the standard MODIS Land Integerized Sinusoidal tile map projection. Each processing tile is in size of 4800 pixels by 4800 pixels, representing approximately 1200 km by 1200 km land region. In addition to the GPP product, SLOPE PAR, SANIRV, and C4 fraction, along with their uncertainties, are also released.

Jiang, Chongya↗

DOE OSTI2023

Evaluation of E3SM land model snow simulations over the western United States

Abstract. Seasonal snow has crucial impacts on climate, ecosystems, and humans, but it is vulnerable to global warming. The land component (ELM) of the Energy Exascale Earth System Model (E3SM) mechanistically simulates snow processes from accumulation, canopy interception, compaction, and snow aging to melt. Although high-quality field measurements, remote sensing snow products, and data assimilation products with high spatio-temporal resolution are available, there has been no systematic evaluation of the snow properties and phenology in ELM. This study comprehensively evaluates ELM snow simulations over the western United States at 0.125∘ resolution during 2001–2019 using the Snow Telemetry (SNOTEL) in situ networks, MODIS remote sensing products (i.e., MCD43 surface albedo product), the spatially and temporally complete (STC) snow-covered area and grain size (MODSCAG) and MODIS dust and radiative forcing in snow (MODDRFS) products (STC-MODSCAG/STC-MODDRFS), and the snow property inversion from remote sensing (SPIReS) product and two data assimilation products of snow water equivalent and snow depth – i.e., University of Arizona (UA) and SNOw Data Assimilation System (SNODAS). Overall the ELM simulations are consistent with the benchmarking datasets and reproduce the spatio-temporal patterns, interannual variability, and elevation gradients for different snow properties including snow cover fraction (fsno), surface albedo (αsur) over snow cover regions, snow water equivalent (SWE), and snow depth (Dsno). However, there are large biases of fsno with dense forest cover and αsur in the Rocky Mountains and Sierra Nevada in winter, compared to the MODIS products. There are large discrepancies of snow albedo, snow grain size, and light-absorbing particle-induced snow albedo reduction between ELM and the MODIS products, attributed to uncertainties in the aerosol forcing data, snow aging processes in ELM, and remote sensing retrievals. Against UA and SNODAS, ELM has a mean bias of −20.7 mm (−35.9 %) and −20.4 mm (−35.5 %), respectively, for spring, and −13.8 mm (−27.8 %) and −10.2 mm (−22.2 %), respectively, for winter. ELM shows a relatively high correlation with SNOTEL SWE, with mean correlation coefficients of 0.69 but negative mean biases of −122.7 mm. Compared to the snow phenology of STC-MODSCAG and SPIReS, ELM shows delayed snow accumulation onset dates by 17.3 and 12.4 d, earlier snow end dates by 35.5 and 26.8 d, and shorter snow durations by 52.9 and 39.5 d, respectively. This study underscores the need for diagnosing model biases and improving ELM representations of snow properties and snow phenology in mountainous areas for more credible simulation and future projection of mountain snowpack.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Multi-scale integration of satellite remote sensing improves characterization of dry-season green-up in an Amazon tropical evergreen forest

In tropical forests, leaf phenology-particularly the pronounced dry-season green-up-strongly regulates biogeochemical cycles of carbon and water fluxes. However, uncertainties remain in the understanding of tropical forest leaf phenology at different spatial scales. Phenocams accurately characterize leaf phenology at the crown and ecosystem scales but are limited to a few sites and time spans of a few years. Time-series satellite observations might fill this gap, but the commonly used satellites (e.g. MODIS, Landsat and Sentinel-2) have resolutions too coarse to characterize single crowns. To resolve this observational challenge, we used the PlanetScope constellation with a 3m resolution and near daily nadir-view coverage. We first developed a rigorous method to cross-calibrate PlanetScope surface reflectance using daily BRDF-adjusted MODIS as the reference. We then used linear spectral unmixing of calibrated PlanetScope to obtain dry-season change in the fractional cover of green vegetation (GV) and non-photosynthetic vegetation (NPV) at the PlanetScope pixel level. We used the Central Amazon Tapajos National Forest k67 site, as all necessary data (from field to phenocam and satellite observations) was available. For this proof of concept, we chose a set of 22 dates of PlanetScope measurements in 2018 and 16 in 2019, all from the six drier months of the year to provide the highest possible cloud-free temporal resolution. Our results show that MODIS-calibrated dry-season PlanetScope data (1) accurately assessed seasonal changes in ecosystem-scale and crown-scale spectral reflectance; (2) detected an increase in ecosystem-scale GV fraction (and a decrease in NPV fraction) from June to November of both years, consistent with local phenocam observations with R 2 around 0.8; and (3) monitored large seasonal trend variability in crown-scale NPV fraction. Finally, our results highlight the potential of integrating multi-scale satellite observations to extend fine-scale leaf phenology monitoring beyond the spatial limits of phenocams.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Seasonal and long-term variations in leaf area of Congolese rainforest

It is important to understand temporal and spatial variations in the structure and photosynthetic capacity of tropical rainforests in a world of changing climate, increased disturbances and human appropriation. The equatorial rainforests of Central Africa are the second largest and least disturbed of the biodiversly-rich and highly productive rainforests on Earth. Currently, there is a dearth of knowledge about the phenological behavior and long-term changes that these forests are experiencing. Here, this study reports on leaf area seasonality and its time trend over the past two decades as assessed from multiple remotely sensed datasets. Seasonal variations of leaf area in Congolese forests derived from MODIS data co-vary with the bimodal precipitation pattern in this region, with higher values during the wet season. Independent observational evidence derived from MISR and EPIC sensors in the form of angular reflectance signatures further corroborate this seasonal behavior of leaf area. The bimodal patterns vary latitudinally within this large region. Two sub-seasonal cycles, each consisting of a dry and wet season, could be discerned clearly. These exhibit different sensitivities to changes in precipitation. Contrary to a previous published report, no widespread decline in leaf area was detected across the entire extent of the Congolese rainforests over the past two decades with the latest MODIS Collection 6 dataset. Long-term precipitation decline did occur in some localized areas, but these had minimal impacts on leaf area, as inferred from MODIS and MISR multi-angle observations.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Photosynthesis phenology, as defined by solar-induced chlorophyll fluorescence, is overestimated by vegetation indices in the extratropical Northern Hemisphere

Vegetation phenology is highly sensitive to climate change, although the data and methods used to estimate key phenological states can influence this sensitivity. Because of its direct relation to leaf photosynthetic carbon uptake, remotely sensed solar-induced chlorophyll fluorescence (SIF) can provide new insight assessing changes in vegetation phenology. In this work, we investigated the potential of using a SIF time series product named contiguous SIF (CSIF) to estimate spring, summer, and autumn phenology in the extratropical Northern Hemisphere (>30°N) and compared the results with those based on Moderate Resolution Imaging Spectroradiometer (MODIS) Normalized Difference Vegetation Index (NDVI) for the period 2001–2018. Overall, we found similar spatial patterns in phenological states. However, specific dates of key phenological events differed when using CSIF vs. MODIS NDVI data. NDVI data indicated that the growing season started earlier (by 10.1 days on average) and ended later (11.5 days on average) relative to CSIF data. This implies that actual periods for photosynthetic activity are shorter (by 21.6 days on average) than those estimated from vegetation indices more directly related to changes in canopy structure. These large differences between results from NDVI and that from CSIF suggest that vegetation indices such as NDVI seem to overestimate the period for active photosynthesis over the extratropical Northern Hemisphere. Furthermore, while phenology of the early growing season is dominated by temperature for both NDVI and CSIF data, phenology of the late growing season is mainly controlled by temperature for NDVI but by precipitation for CSIF. Our findings were further confirmed by other SIF (GOME-2 SIF) and vegetation index (MODIS EVI) datasets. Phenology modes in Earth system modelling are often parameterized using leaf unfolding and senescence from either station or satellite observations. Our results imply that canopy structure-based parameterization schemes may have overestimated photosynthesis active period, and thus productivity responses. We conclude that SIF data offers a novel and unique approach for assessing phenological change - one that is more directly tied to the carbon cycle and how it is being influenced by climate change.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

The underappreciated importance of solar radiation in constraining spring phenology of temperate ecosystems in the Northern and Eastern United States

Spring phenology of temperate ecosystems is highly sensitive to climate change, generating various impacts on many important terrestrial surface biophysical processes. Although various prognostic models relying on environmental variables of temperature and photoperiod have been developed for spring phenology, comprehensive ecosystem-scale evaluations over large landscapes and long-time periods remain lacking. Further, environmental variables other than temperature and photoperiod might also importantly constrain spring phenology modelling but remain under-investigation. To address these issues, we leveraged around 20-years datasets of environmental variables (from Daymet and GLDAS products) and the spring phenology metric (i.e., the greenup date) respectively derived from MODIS and PhenoCams across 108 sites in the Northern and Eastern United States. We firstly cross-compared MODIS-derived greenup date with official PhenoCams product with high accuracy (R 2 = 0.70). Then, we evaluated the three prognostic models (i.e., Growing Degree Date (GDD), Sequential (SEQ) and optimality-based (OPT)) with MODIS-derived spring phenology, assessed the model residuals and their associations with soil moisture, rainfall, and solar radiation, and revised the two photoperiod-relevant models (SEQ, OPT) by replacing the daylength variable with solar radiation, which was found to contribute the most to model residuals. We found that 1) all models demonstrated good capability in characterizing spring phenology, with OPT performing the best (RMSE = 8.04 ± 5.05 days), followed by SEQ (RMSE = 10.57 ± 7.77 days) and GDD (RMSE = 10.84 ± 8.42 days), 2) all models displayed high model residuals showing tight correlation with solar radiation (r = 0.45–0.75), and 3) the revised models that included solar radiation significantly performed better with an RMSE reduction by 22.08%. Such results are likely because solar radiation better constrains early growing season plant photosynthesis than photoperiod, supporting the hypothesis of spring phenology as an adaptive strategy to maximize photosynthetic carbon gain (approximated by solar radiation) while minimizing frost damage risk (captured by temperature). Collectively, our study here reveals the underappreciated importance of solar radiation in constraining spring phenology of temperate ecosystems, and suggests ways to improve spring phenology modelling and other phenology-related ecological processes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Daily, 30 m Resolution NDSI Data for the East River Watershed, CO for 2000-2020

This dataset contains daily Normalized Difference Snow Index (NDSI) values at 30 m spatial resolution for the East River watershed in Colorado, USA. The temporal range of these data includes water years 2001-2020. These data were created using the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM). This model fuses low spatial and high temporal resolution data from MODIS (500 m, daily) with high spatial and low temporal resolution data from Landsat (30 m, 16 days) to create a 30m synthetic daily snow product. This product allows for the analysis of historical snow covered area trends in the East River Watershed at fine spatiotemporal resolutions where it was not available previously. This research was performed as a part of the Department of Energy’s Subsurface Biogeochemical Research Program with the primary intent of better understanding the timing and spatial patterns of water delivery to the Critical Zone in mountain watersheds. Each .zip file contains one "water year" of data (October 1 - September 30; i.e., water year 2010 starts October 1, 2010 and ends September 30, 2011). Each zip file contains the following: STARFM daily Normalized Difference Snow Index (NDSI) fusion data files in GeoTiff format with one layer for each day between Landsat data acquisition dates (i.e., for dates of Landsat acquisition, the Landsat image is included for that date). The study area is located in an area of Landsat path overlap, so Landsat dates acquisitions are every 7-9 days. Landsat NDSI files containing the high spatial (30m), low temporal (7-9 days due to Landsat path overlap) resolution data used as input to STARFM in GeoTiff format with one layer for each day. Dates for which no Landsat data were obtained are included as NoData layers. MODIS NDSI files containing the high temporal (daily), low spatial (500m) resolution data used as input to STARFM in GeoTiff format with one layer for each day. Please note the MODIS data were resampled to 30m pixels for input into the STARFM model. The data have a scale factor of 10,000 and a no data value of -32767. The projection of all datasets is WGS 84 (EPSG: 4326), which has a latitude/longitude based degree resolution of 0.0002694946 X 0.0002694946, and approximates to the 30 m spatial resolution mentioned above. The Layer Index files in .csv format. They contain information for each layer in the above GeoTiff files regarding the corresponding date for each layer, the fraction of pixels in the image that contain valid data (missing data is due to either cloud cover or poor data quality; these values are not percent snow cover). Dates of Landsat overpass are indicated in these files. If no Landsat data were able to be obtained due to cloud cover or lack of Landsat Tier 1 data available on Google Earth Engine, this is also noted.

EARTH SCIENCE > CRYOSPHERE > SNOW/ICE↗

DOE OSTI2022

LAI, EVI, NDVI, and kNDVI in 23 pantropical forests affected by 21 cyclones

Statement of purpose: Cyclones alter the function and composition of tropical forests, making effects of intensifying cyclones on carbon-rich forests a critical topic of study. Here, we quantified cyclone-induced damage and recovery of 21 cyclone disturbances affecting 23 pantropical forest sites between 1988-2017 utilizing leaf area index (LAI), enhanced vegetation index (EVI), normalized difference vegetation index (NDVI), and transformed NDVI (kNDVI) values from Google Earth Engine. Field observations collected in a meta-analysis (Bomfim et al., 2022, in review) were used to ground-truth and test effects of soil resource availability and disturbance factors on damage and recovery. This meta-analysis also served as the basis to begin vegetation index extraction, utilizing unique site and date combinations, from tropical forests effect by cyclone disturbances. We began collecting NDVI (5km resolution) from the NOAA Climate Data Record (CDR) of AVHRR Normalized Difference Vegetation Index (NDVI), Version 5 data product (Vermote, 2019) for all case studies included, 42. Next, we began extracting Landsat data from Landsat 4, 5, and 8, courtesy of the U.S. Geological Survey, in search of higher resolution data. We selected a 3 by 3 Landsat pixel area, leading to a 90m resolution data extraction. The specific imagery used includes Landsat 4 USGS Landsat 4 TM Collection 1 Tier 1 TOA (top of atmosphere) Reflectance, Landsat 5 USGS Landsat 5 TM (thematic mapper) Collection 1 Tier 1 TOA Reflectance, and Landsat 8 USGS Landsat 8 Collection 1 Tier 1 TOA Reflectance. Within Google Earth Engine, we selected the date and location (latitude and longitude), calculated NDVI, kNDVI, and EVI utilizing Landsat bands (see metadata_NGEE-tropics_cyclones), and extracted post- and pre-cyclone values for each case study to calculate cyclone-induced change in the vegetative index. Due to limited spatial resolution of Landsat remote sensing data, MODIS products were investigated next. First, the MOD13Q1.006 Terra Vegetation Indices 16-Day Global 250m product was used to extract 250m EVI and NDVI (Didan, 2015) and then the MCD15A3H.006 MODIS Leaf Area Index/FPAR 4-Day Global 500m product product was used to extract LAI 500m (Myneni et al., 2015). Pre- and post-cyclone values, change in the vegetative index, and standard deviation for all values are included in the main csv (see case_study_data.csv) for all vegetative indices collected, including LAI 500m, EVI 250m, NDVI 250m, NDVI 90m, kNDVI 90m, EVI 90m, and NDVI 5km. Lastly, recovery values were calculated utilizing a standardization method (see metadata_NGEE-tropics_cyclones) and recovery values for MODIS (see MODIS_recovery.csv) and Landsat (Landsat_recovery.csv) data are included.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Multiscale ACI Satellite Database

The SATELLITE_EAGLES_PNNL NetCDF dataset contains a suite of satellite- and reanalysis-derived atmospheric and surface parameters on a regular latitude–longitude grid. The dataset includes core geophysical fields such as land fraction, aerosol optical depth at multiple wavelengths (465, 550, 667, and 865 nm), sea surface temperature, estimated inversion strength, and various thermodynamic and dynamic quantities (e.g., relative humidity, vertical velocity, boundary-layer height, and surface fluxes) from both MERRA and ERA reanalysis products, provided as daily-mean and instantaneous values. A major component of the dataset consists of MODIS-retrieved cloud microphysical properties, including cloud droplet number concentration, cloud effective radius, optical thickness, and liquid water path, provided for three compositing regimes (“All,” “Q06,” and “G18”). Corresponding cloud-top parameters—temperature, height, and pressure—along with total and domain-mean cloud fraction fields are also included. The file further integrates additional satellite data from AMSR-E (for cloud water, rain water, and surface precipitation retrievals) and CERES (for top-of-atmosphere radiative fluxes, cloud fractions, and albedo). This dataset is designed to evaluate aerosol–cloud interactions in warm clouds, emphasizing the use of MODIS for deriving cloud droplet number concentration and liquid water path statistics. The complementary satellite and reanalysis fields are co-located and time-matched to the same instantaneous MODIS observations, enabling consistent comparisons between cloud properties, aerosol loading, and large-scale meteorological conditions. The dataset is recently featured in Christensen et al. (2025), Machine Learning Reveals Strong Grid-Scale Dependence in the Satellite Nd–LWP Relationship, EGUsphere [preprint], https://doi.org/10.5194/egusphere-2025-3850, 2025.

Christensen, Matthew [Pacific Northwest National L↗

DOE OSTI2025

Multiscale ACI Satellite Database

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Droplet collection efficiencies inferred from satellite retrievals constrain effective radiative forcing of aerosol–cloud interactions

Abstract. Process-oriented observational constraints for the anthropogenic effective radiative forcing due to aerosol–cloud interactions (ERFaci) are highly desirable because the uncertainty associated with ERFaci poses a significant challenge to climate prediction. The contoured frequency by optical depth diagram (CFODD) analysis supports the evaluation of model representation of cloud liquid-to-rain conversion processes because the slope of a CFODD, generated from joint MODerate Resolution Imaging Spectroradiometer (MODIS)-CloudSat cloud retrievals, provides an estimate of cloud droplet collection efficiency in single-layer warm liquid clouds. Here, we present an updated CFODD analysis as an observational constraint on the ERFaci due to warm rain processes and apply it to the U.S. Department of Energy's Energy Exascale Earth System Model version 2 (E3SMv2). A series of sensitivity experiments shows that E3SMv2 droplet collection efficiencies and ERFaci are highly sensitive to autoconversion, i.e., the rate of mass transfer from cloud liquid to rain, yielding a strong correlation between the CFODD slope and the shortwave component of ERFaci (ERFaciSW; Pearson's R=-0.91). E3SMv2's CFODD slope (0.20 ± 0.04) is in agreement with observations (0.20 ± 0.03). The strong sensitivity of ERFaciSW to the CFODD slope provides a useful constraint on highly uncertain warm rain processes, whereby ERFaciSW, constrained by MODIS-CloudSat, is estimated by calculating the intercept of the linear association between the ERFaciSW and the CFODD slopes, using the MODIS-CloudSat CFODD slope as a reference.

Beall, Charlotte M. (ORCID:0000000311370835)↗

DOE OSTI2025

Evaluation of cloud height, optical thickness, and phase retrievals from the CHROMA algorithm applied to Sentinel-3 OLCI data

We previously developed the Cloud Height Retrieval from O 2 Molecular Absorption (CHROMA) algorithm for the Ocean Color Instrument (OCI) on the new NASA Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) mission. Here, we apply CHROMA to observations from the Ocean Land Colour Instrument (OLCI) to guide expectations for PACE, as it will take some time to obtain large-scale validation data for OCI. We use cloud top height (CTH), phase, and (for liquid clouds) cloud optical thickness (COT) data from the ground-based Atmospheric Radiation Measurement (ARM) network to evaluate the OLCI retrievals. We found that OLCI and Moderate Resolution Imaging Spectroradiometer (MODIS) CTH compare similarly well to the ARM reference. OLCI has a tendency to underestimate CTH as CTH increases, and algorithm assumptions about cloud geometric thickness may contribute to this. ARM COT from multifilter shadowband radiometers (MFRSR) and Sun photometers are well-correlated with one another, albeit with a roughly 30 % offset on average; OLCI and MODIS COT agree more closely with the MFRSR data. OLCI retrieval uncertainty estimates show skill at telling low-uncertainty cases from high-uncertainty ones, although CTH uncertainties are underestimated. Additionally, we compare the OLCI data to satellite retrievals based on thermal infrared measurements from MODIS and Sea and Land Surface Temperature Radiometer (SLSTR) data. Differences are broadly consistent with physical expectations based on the A-band vs. thermal techniques, although one key challenge in such aggregated comparisons is different cloud masking sensitivities and algorithm failure rates meaning additional sampling differences are introduced. We conclude by discussing the transition to and possible enhancements for PACE OCI.

Sayer, Andrew M. [Univ. of Maryland Baltimore Coun↗

DOE OSTI2023

Snow-corrected vegetation indices for improved gross primary productivity assessment in North American evergreen forests

North American evergreen forests cover large areas and influence the global carbon cycle. Satellite remote sensing has been used to track the phenology of ecosystem photosynthesis of these forests by detecting variation in vegetation optical properties associated with physiological and structural features, and most of these methods have been closely tied to vegetation greenness. However, in evergreens, the application of satellite data to monitor photosynthetic phenology is often limited by the lack of sensitivity of greenness-based indices. In this study, we identified 47 evergreen forest flux sites in North America that had MODIS observation overlapping with the flux tower records. We then calculated four vegetation indices using MODIS MAIAC data (MCD19A1), including NDVI, CCI, NIRv, and kNDVI, for the 47 flux sites and evaluated relationships between gross primary productivity (GPP) and vegetation indices across the North American evergreen forests. Our results showed that snow had substantial effects on the performance of all vegetation indices in tracking GPP phenology, particularly in the early spring when rapid changes occurred to both GPP and snow cover. Furthermore, different vegetation indices were affected differently, indicating contradictory and confounding effects of snow on these indices. After correcting for the snow effects, both CCI and NIRv performed well in tracking GPP phenology, albeit for different reasons. CCI is sensitive to seasonal changes in the relative levels of chlorophyll and carotenoid pigments, which are closely tied to GPP phenology in evergreens. NIRv is sensitive to the absorbed photosynthetically active radiation and to the contribution of deciduous components to the overall optical properties. We also found that correlations between GPP and vegetation indices varied among ecoregions and climate classes. In general, regions with pronounced seasonal GPP patterns had stronger correlations between GPP and greenness-based indices than regions with weaker seasonal GPP patterns. These biome differences were less pronounced for CCI. The snow artifacts and complementary vegetation index effects reported here should be considered in any large-scale studies of GPP using reflectance-based indices from optical satellites.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Persistent global greening over the last four decades using novel long-term vegetation index data with enhanced temporal consistency

Advanced Very High-Resolution Radiometer (AVHRR) satellite observations have provided the longest global daily records from 1980s, but the remaining temporal inconsistency in vegetation index datasets has hindered reliable assessment of vegetation greenness trends. To tackle this, we generated novel global long-term Normalized Difference Vegetation Index (NDVI) and Near-Infrared Reflectance of vegetation (NIRv) datasets derived from AVHRR and Moderate Resolution Imaging Spectroradiometer (MODIS). We addressed residual temporal inconsistency through three-step post processing including cross-sensor calibration among AVHRR sensors, orbital drifting correction for AVHRR sensors, and machine learning-based harmonization between AVHRR and MODIS. After applying each processing step, we confirmed the enhanced temporal consistency in terms of detrended anomaly, trend and interannual variability of NDVI and NIRv at calibration sites. Our refined NDVI and NIRv datasets showed a persistent global greening trend over the last four decades (NDVI: 0.0008 yr -1 ; NIRv: 0.0003 yr -1 ), contrasting with those without the three processing steps that showed rapid greening trends before 2000 (NDVI: 0.0017 yr -1 ; NIRv: 0.0008 yr -1 ) and weakened greening trends after 2000 (NDVI: 0.0004 yr -1 ; NIRv: 0.0001 yr -1 ). These findings highlight the importance of minimizing temporal inconsistency in long-term vegetation index datasets, which can support more reliable trend analysis in global vegetation response to climate changes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Satellite-based aerosol optical depth estimates over the continental U.S. during the 2020 wildfire season: Roles of smoke and land cover

Wildfires produce smoke that can affect an area >1000 times the burn extent, with far-reaching human health, ecologic, and economic impacts. Accurately estimating aerosol load within smoke plumes is therefore crucial for understanding and mitigating these impacts. We evaluated the effectiveness of the latest Collection 6.1 MODIS Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm in estimating aerosol optical depth (AOD) across the U.S. during the historic 2020 wildfire season. We compared satellite-based MAIAC AOD to ground-based AERONET AOD measurements during no-, light-, medium-, and heavy-smoke conditions identified using the Hazard Mapping System Fire and Smoke Product. This smoke product consists of maximum extent smoke polygons digitized by analysts using visible band imagery and classified according to smoke density. We also examined the strength of the correlations between satellite- and ground-based AOD for major land cover types under various smoke density levels. MAIAC performed well in estimating AOD during smoke-affected conditions. Correlations between MAIAC and AERONET AOD were strong for medium- (r = 0.91) and heavy-smoke (r = 0.90) density, and MAIAC estimates of AOD showed little bias relative to ground-based AERONET measurements (normalized mean bias = 3 % for medium, 5 % for heavy smoke). During two high AOD, heavy smoke episodes, MAIAC underestimated ground-based AERONET AOD under mixed aerosol (i.e., smoke and dust; median bias = −0.08) and overestimated AOD under smoke-dominated (median bias = 0.02) aerosol. MAIAC most overestimated ground-based AERONET AOD over barren land (mean NMB = 48 %). Our findings indicate that MODIS MAIAC can provide robust estimates of AOD as smoke density increases in coming years. Increased frequency of mixed aerosol and expansion of developed land could affect the performance of the MAIAC algorithm in the future, however, with implications for evaluating wildfire-associated health and welfare effects and air quality standards.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Evaluation of simulations of near-surface variables using the regional climate model CCLM for the MOSAiC winter period

The ship-based experiment MOSAiC 2019/2020 was carried out during a full year in the Arctic and yielded an excellent data set to test the parameterizations of ocean/sea-ice/atmosphere interaction processes in regional climate models (RCMs). In the present paper, near-surface data during MOSAiC are used for the verification of the RCM COnsortium for Small-scale MOdel–Climate Limited area Mode (COSMO-CLM or CCLM). CCLM is used in a forecast mode (nested in ERA5) for the whole Arctic with 15 km resolution and is run with different configurations of sea ice data. These include the standard sea ice concentration taken from passive microwave data with around 6 km resolution, sea ice concentration from Moderate Resolution Imaging Spectroradiometer (MODIS) thermal infrared data and MODIS sea ice lead fraction data for the winter period. CCLM simulations show a good agreement with the measurements. Relatively large negative biases for temperature occur for November and December, which are likely associated with a too large ice thickness used by CCLM. The consideration of sea ice leads in the sub-grid parameterization in CCLM yields improved results for the near-surface temperature. ERA5 data show a large warm bias of about 2.5°C and an underestimation of the temperature variability.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Comparing Coarse-Resolution Land Surface Temperature Products over Western Australia

Satellite-derived land surface temperature (LST) has commonly been used to monitor global temperature changes. The MODIS MYD11A2 product is the most common coarse-resolution product used for this purpose. An updated MODIS product (MYD21A2) and new VIIRS (VNP21A2) product have also recently become available. This study analyses eight-day, quality-controlled, LST imagery over Western Australia (WA) for the three products for an urban and a non-urban area for the years 2013, 2016, and 2019. An analysis of the data indicates that (i) the averaged daytime LST difference between the three products for Perth city over the three years was 1.32 °C, while at night it was 0.89 °C; (ii) the averaged daytime difference over the Kimberley region was 7.02 °C with a night average difference of 2.39 °C; and (iii) both the MYD21A2 and VNP21A2 products still appear to record anomalous monthly LST values, particularly in the humid Kimberley monsoonal months. The overall objective of the National Aeronautics and Space Administration (NASA) is to ensure that the LST values of the two satellite system products are comparable, so evidence of LST value differences will require further investigation, especially if the older product is to be replaced by the newer systems.

47 OTHER INSTRUMENTATION↗