Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “covariates”

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 613 records · Page 34

Artificial neural networks estimate evapotranspiration for Miscanthus × giganteus as effectively as empirical model but with fewer inputs

Estimating actual evapotranspiration (ET) is particularly crucial for addressing how vegetation affects the water balance of ecosystems. ET estimation can be complex with empirical models due to their many parameters and reliance on aridity. In contrast, artificial neural networks (ANNs) could potentially estimate ET with fewer and more common meteorological parameters. In this study, we trained two ANNs, one using a feed-forward approach (FFN) and the other a nonlinear auto-regressive network (NARX), to predict ET and compared them to the commonly used empirical model Granger and Gray (GG). We trained our models on a nine-year eddy covariance (EC) dataset for Miscanthu s × giganteus ( M . × giganteus ) from Illinois (UIEF), then tested them using out-of-sample data from both UIEF and a different location in Iowa (SABR) to compare the accuracy of FFN, NARX, and GG models in estimating daily ET. A combination of air temperature (T a ) and solar radiation (R s ) was chosen as inputs due to the highest R 2 for FFN (R 2 = 0.79, 0.81, and 0.79 for training, testing, and validation, respectively) and only T a for NARX (R 2 = 0.70 for out-of-sample validation). The predictive power of the FFN model was superior to the NARX and GG models at the UIEF site (R 2 = 0.84, 0.70, and 0.83 for out-of-sample validation, respectively). Our analysis showed that ANN approaches are as accurate as empirical approaches for estimating ET but use fewer inputs.

54 ENVIRONMENTAL SCIENCES↗

Mineral-Associated Organic Matter Concentration Beneath Northern Temperate Trees Varies by Mycorrhizal Type and Leaf Habit

Mycorrhizal fungi are important drivers of soil organic matter dynamics, but it can be difficult to isolate the effects of the fungi themselves from covarying traits of their host trees. For example, many trees with an evergreen leaf habit associate with ectomycorrhizal (ECM) fungi, while many deciduous tree species associate with arbuscular mycorrhizal (AM) fungi. Because leaf habit influences the quantity and quality of organic matter inputs to soil, it is often an important factor in soil carbon and nitrogen dynamics, and thus can mask the effects of mycorrhizal fungi on soil organic matter processes. We evaluated how tree mycorrhizal associations and leaf habit separately influence the amount and composition of mineral-associated organic matter (MAOM) and particulate organic matter (POM) in forest soils in New Hampshire and Vermont, USA. We measured carbon (C) and nitrogen (N) concentrations and C/N ratios of three soil density fractions beneath six tree species that vary in mycorrhizal association and leaf habit. We found lower concentrations of MAOM C and N beneath evergreen vs. deciduous trees, but only for tree species associating with AM fungi. Further, MAOM C/N was higher beneath evergreen trees and beneath trees with ECM fungi rather than AM fungi. Furthermore, these results add to the growing body of support for mycorrhizal fungi as mediators of soil organic matter dynamics, suggesting that the MAOM fraction is more sensitive to leaf habit beneath AM-associated versus ECM-associated trees. Because MAOM decomposition is thought to be less responsive than POM decomposition to changes in soil temperature and moisture, differences in the tendency of AM- and ECM-dominated forests to support MAOM formation and persistence may lead to systematic differences in the response of these forest types to ongoing climate change.

54 ENVIRONMENTAL SCIENCES↗

Multistage distributionally robust mixed-integer programming with decision-dependent moment-based ambiguity sets

We study multistage distributionally robust mixed-integer programs under endogenous uncertainty, where the probability distribution of stage-wise uncertainty depends on the decisions made in previous stages. We first consider two ambiguity sets defined by decision-dependent bounds on the first and second moments of uncertain parameters and by mean and covariance matrix that exactly match decision-dependent empirical ones, respectively. For both sets, we show that the subproblem in each stage can be recast as a mixed-integer linear program (MILP). Moreover, we extend the general moment-based ambiguity set in to the multistage decision-dependent setting, and derive mixed-integer semidefinite programming (MISDP) reformulations of stage-wise subproblems. We develop methods for attaining lower and upper bounds of the optimal objective value of the multistage MISDPs, and approximate them using a series of MILPs. We deploy the Stochastic Dual Dynamic integer Programming (SDDiP) method for solving the problem under the three ambiguity sets with risk-neutral or risk-averse objective functions, and conduct numerical studies on multistage facility-location instances having diverse sizes under different parameter and uncertainty settings. Furthermore, our results show that the SDDiP quickly finds optimal solutions for moderate-sized instances under the first two ambiguity sets, and also finds good approximate bounds for the multistage MISDPs derived under the third ambiguity set. We also demonstrate the efficacy of incorporating decision-dependent distributional ambiguity in multistage decision-making processes.

97 MATHEMATICS AND COMPUTING↗

Estimating and Evaluating Roughness Length and Displacement Height in Heterogeneous Urban Environments

The roughness length (z 0 ) and displacement height (z d ) are essential surface-layer parameters in numerical models (e.g., weather, climate, wall-modeled LES, etc.). This work evaluates the consistency of z 0 and z d estimates from morphometric and anemometric methods using data from two eddy-covariance flux towers (AmeriFlux US-INg and US-INc) in Indianapolis, IN. Results show inconsistencies in estimated z 0 and z d values depending on the chosen method. The two evaluated anemometric methods estimate non-physical values of z d when compared to roughness elements surrounding both towers. Additionally, predictions of mean wind speed using surface-layer similarity theory with morphometric estimates exhibit a bias during near-neutral and stable conditions relative to observations. The overestimation of mean wind speed by surface layer similarity theory is consistent with previous observational and modeling studies in urban areas, suggesting that the application of similarity theories to urban environments may have limitations. Differentiation of vegetation from built structures appears to impact morphometric z 0 and z d estimates, particularly where vegetation is abundant; however, it has little impact on correcting biases in the similarity theory. Specifically, we find that existing similarity theories using morphometric estimates underestimate integral velocity and length scales, and the degree of underestimation depends on the stability conditions. Accounting for the degree of anisotropy in surface-layer turbulence helps reduce the biases between similarity theories and observations during unstable conditions, but not in near-neutral cases. Future work is needed to identify the cause of such biases for near-neutral conditions.

Aerodynamic roughness length↗

Realization of the Werner-Holevo and Landau-Streater Quantum Channels for Qutrits on Quantum Computers.

We realize the Landau-Streater (LS) and Werner-Holevo (WH) quantum channels for qutrits on IBM quantum computers. These channels correspond to the interaction between a qutrit and its environment that results in the globally unitarily covariant qutrit transformation violating the multiplicativity of the maximal p-norm. Our realization of the LS and WH channels is based on embedding the qutrit states into states of two qubits and using the single-qubit and two-qubit CNOT gates to implement the specific interaction. We employ the standard quantum gates, hence the developed algorithm suits any quantum computer. We run our algorithm on a 5-qubit computer and a 20-qubit computer, as well as on a simulator. We quantify the quality of the implemented channels comparing their action on different input states with theoretical predictions. The overall efficiency is quantified by the fidelity between the theoretical and experimental Choi states implemented on the 20-qubit computer.

completely-positive trace-preserving map↗

Remote quantification of Cm(III) and HNO 3 by fluorescence spectroscopy and chemometrics

A unique approach to remotely quantify Cm(III) (0–100 µg mL −1 ) in HNO 3 (1–12 M) using steady-state laser fluorescence spectroscopy and multivariate regression models was developed. Photoluminescence is amenable to remote measurements using fiber-optic cables and is sensitive to numerous lanthanide and actinide species. In-line measurements can provide feedback to support complex processing in harsh environments (e.g., hot cells) to help guide and optimize radiochemical separations. In this work, Cm(III) spectra were acquired remotely in a glove box as a function of HNO 3 concentration to better understand spectral characteristics and evaluate the utility of multivariate regression models in this system. Furthermore, the Cm(III) fluorescence peak shape, width, position, and intensity changed significantly as a function of HNO 3 concentration, likely because of the displacement of emission quenching inner-sphere water molecules and complexation with nitrate ions. Despite significant covariance and nonlinearity in the data, a D-optimal design strategy successfully minimized training set sample size and was used to build effective partial least squares regression models for Cm(III) and HNO 3 concentrations without a priori knowledge of solution conditions. Chemometrics for modeling complex fluorescence spectra are promising and may find widespread applicability for online analysis in numerous chemical systems found in the nuclear field.

Actinide↗

Estimation of soil classes and their relationship to grapevine vigor in a Bordeaux vineyard: advancing the practical joint use of electromagnetic induction (EMI) and NDVI datasets for precision viticulture

Working within a vineyard in the Pessac Léognan Appellation of Bordeaux, France, this study documents the potential of using simple statistical methods with spatially-resolved and increasingly available electromagnetic induction (EMI) geophysical and normalized difference vegetation index (NDVI) datasets to accurately estimate Bordeaux vineyard soil classes and to quantitatively explore the relationship between vineyard soil types and grapevine vigor. First, co-located electrical tomographic tomography (ERT) and EMI datasets were compared to gain confidence about how the EMI method averaged soil properties over the grapevine rooting depth. Then, EMI data were used with core soil texture and soil-pit based interpretations of Bordeaux soil types (Brunisol, Redoxisol, Colluviosol and Calcosol) to estimate the spatial distribution of geophysically-identified Bordeaux soil classes. A strong relationship (r = 0.75, p < 0.01) was revealed between the geophysically-identified Bordeaux soil classes and NDVI (both 2 m resolution), showing that the highest grapevine vigor was associated with the Bordeaux soil classes having the largest clay fraction. The results suggest that within-block variability of grapevine vigor was largely controlled by variability in soil classes, and that carefully collected EMI and NDVI datasets can be exceedingly helpful for providing quantitative estimates of vineyard soil and vigor variability, as well as their covariation. The method is expected to be transferable to other viticultural regions, providing an approach to use easy-to-acquire, high resolution datasets to guide viticultural practices, including routine management and replanting.

54 ENVIRONMENTAL SCIENCES↗

Scalable computations for nonstationary Gaussian processes

Nonstationary Gaussian process models can capture complex spatially varying dependence structures in spatial datasets. However, the large number of observations in modern datasets makes fitting such models computationally intractable with conventional dense linear algebra. In addition, derivative-free or even first-order optimization methods can be very slow to converge when estimating many spatially varying parameters. In this paper, we present a computational framework which couples an algebraic block diagonal plus low-rank covariance matrix approximation with stochastic trace estimation to facilitate the efficient use of second-order solvers for maximum likelihood estimation of Gaussian process models with many parameters. We demonstrate the effectiveness of these methods by simultaneously fitting 192 parameters in the popular nonstationary model of Paciorek and Schervish using 107,600 sea surface temperature anomaly measurements.

97 MATHEMATICS AND COMPUTING↗

Practical Guide to Measuring Wetland Carbon Pools and Fluxes

Abstract Wetlands cover a small portion of the world, but have disproportionate influence on global carbon (C) sequestration, carbon dioxide and methane emissions, and aquatic C fluxes. However, the underlying biogeochemical processes that affect wetland C pools and fluxes are complex and dynamic, making measurements of wetland C challenging. Over decades of research, many observational, experimental, and analytical approaches have been developed to understand and quantify pools and fluxes of wetland C. Sampling approaches range in their representation of wetland C from short to long timeframes and local to landscape spatial scales. This review summarizes common and cutting-edge methodological approaches for quantifying wetland C pools and fluxes. We first define each of the major C pools and fluxes and provide rationale for their importance to wetland C dynamics. For each approach, we clarify what component of wetland C is measured and its spatial and temporal representativeness and constraints. We describe practical considerations for each approach, such as where and when an approach is typically used, who can conduct the measurements (expertise, training requirements), and how approaches are conducted, including considerations on equipment complexity and costs. Finally, we review key covariates and ancillary measurements that enhance the interpretation of findings and facilitate model development. The protocols that we describe to measure soil, water, vegetation, and gases are also relevant for related disciplines such as ecology. Improved quality and consistency of data collection and reporting across studies will help reduce global uncertainties and develop management strategies to use wetlands as nature-based climate solutions.

Bansal, Sheel (ORCID:0000000312331707)↗

Using ground-based lidar data to investigate the water–vapor budget in the daytime atmospheric boundary layer

The moisture advection term in the water–vapor budget equation is investigated with a combination of a vertically-staring water–vapor lidar and Doppler lidar systems. These instruments make it possible to get the mean profile of moisture tendency and the latent heat flux (LHF) divergence. We use data of the Land–Atmosphere Feedback Experiment (LAFE) at the Atmospheric Radiation Measurement (ARM) Program’s Southern Great Plains (SGP) site, Oklahoma, USA, collected on 30 August 2017 between 15 and 24 UTC, which corresponds to 09 to 18 LT. The lidars provide turbulence resolving profiles of moisture and vertical wind fluctuations. The LHF profile is derived from the covariance of these moisture and vertical wind fluctuations. The mean boundary layer height z i is determined from the peak of the moisture variance. The results demonstrate that the combination of two remote sensing instruments can be applied for determining the dominant water–vapor budget terms, namely moisture tendency, latent heat flux divergence and moisture advection.

Advection↗

Using the red chromatic coordinate to characterize the phenology of forest canopy photosynthesis

Vegetation phenology has received increasing attention in climate change research. Near-surface sensing using digital repeat photography has proven to be useful for ecosystem-scale monitoring of vegetation phenology. However, our understanding of the link between phenological metrics derived from digital repeat photography and the phenology of forest canopy photosynthesis is still incomplete, especially for evergreen plant species. Using 49 site-years of digital images from the PhenoCam Network from eight evergreen needleleaf forest (ENF) and six deciduous broadleaf forest (DBF) sites in North America, we explored the potential of the green chromatic (GCC) and red chromatic coordinates (RCC) in tracking forest canopy photosynthesis by comparing camera-derived start- and end-of-growing season (SOS and EOS, respectively) with corresponding estimates derived from eddy covariance-derived daily gross primary productivity (GPP). We found that for DBF sites, both GCC and RCC performed comparable in capturing SOS and EOS. However, similar to earlier studies, GCC had limited potential in capturing GPP-based SOS or EOS for ENF sites. In contrast, we found RCC was a better predictor of both GPP-based SOS and EOS for ENF sites. Environmental and ecological explanations were both provided that phenological transitions derived from RCC were highly correlated with spring and autumn meteorological conditions, as well as having overall higher correlations with phenology based on LAI, a critical variable for describing canopy development. Our results demonstrate that RCC is an underappreciated metric for tracking vegetation phenology, especially for ENF sites where GCC failed to provide reliable estimates for GPP-based SOS or EOS. Our results improve confidence in using digital repeat photography to characterize the phenology of canopy photosynthesis across forest types.

54 ENVIRONMENTAL SCIENCES↗

Energy balance in the DSSAT-CSM-CROPGRO model

One potential way to improve crop growth models is for the models to predict energy balance and evapotranspiration (ET) from first principles, thus serving as a check on “engineered” ET methodology. In this paper, we present new implementations and the results of an energy balance model (EBL) developed by Jagtap and Jones (1989) and then implemented in DSSAT’s CROPGRO (CG-EBL) model by Pickering et al. (1995) as a linked energy balance-photosynthesis model that has not been field-tested until now. The energy balance code computes evapotranspiration and other energy balance components, as well as a canopy air temperature, based on three sources (sunlit leaves, shaded leaves, soil surface). Model performance was evaluated with measured biomass and energy fluxes from two sites in Nebraska, namely, the US-Ne2 irrigated maize-soybean rotation field and the US-Ne3 rainfed maize-soybean rotation field, which are part of the Ameriflux eddy covariance network (htt ps://ameriflux.lbl.gov/sites). After implementing new aerodynamic resistances and the stomatal conductance model of the Ball–Berry–Leuning, crop growth, evapotranspiration and soil temperature were simulated well by the EBL model. The EBL improved ET predictions slightly over the often-used FAO56 method [Penman–Monteith (Allen et al., 1998)] for 4 of the 5 years evaluated for both irrigated and rainfed conditions. Further, a significant improvement was achieved using EBL for the simulation of soil temperature at the various depths compared to STEMP, the original subroutine in DSSAT for simulating soil temperature. Compared to the other available DSSAT methods, the EBL explicitly simulates the impacts of crop morphology, physiology and management on the crop’s environment and energy and mass exchange, which in turn directly affect the water use and irrigation requirements, phenology, photosynthesis, growth, sterility, and yield of the crop.

54 ENVIRONMENTAL SCIENCES↗

Solar position confounds the relationship between ecosystem function and vegetation indices derived from solar and photosynthetically active radiation fluxes

Vegetation indices derived from solar and photosynthetically active radiation (PAR) sensors (i.e. radiation derived) have been under-utilized in inferring ecosystem function, despite measurement capability at hundreds of sites. This under-utilization may be attributed to reported mismatches among the seasonality of radiation- and satellite-derived vegetation indices and canopy photosynthesis; herein referred to as measurement biases. In this work, biases in radiation derived reflectance and vegetation indices were assessed using a decadal record of satellite and ground based spectroradiometer data, ecosystem phenology and CO 2 fluxes, and radiation derived vegetation indices (i.e. the Normalized Difference Vegetation Index [NDVI], the two band Enhanced Vegetation Index [EVI2]) from a high latitude tundra site (i.e. Imnaviat). At Imnaviat, we found poor correspondence between the three types of reflectance and vegetation indices, especially during the latter part of the growing season. Radiation derived vegetation indices resulted in incorrect estimates of phenological timing of up to a month and poor relationships with canopy photosynthesis (i.e. Gross Ecosystem Exchange (GEE)). These mismatches were attributed to solar position (i.e. solar zenith and azimuth angle) and a method, based on the diel visible and near-infrared albedo variation, was developed to improve the performance of the vegetation indices. The ability of radiation derived vegetation indices to infer GEE and phenological dates drastically improved once radiation derived vegetation indices were corrected for solar position associated biases at Imnaviat. Moreover, radiation derived vegetation indices became better aligned with MODerate resolution Imaging Spectroradiometer (MODIS) satellite estimates after solar position associated biases were corrected at Imnaviat and at 25 Fluxnet sites (~90 site years) across North America. Corrections developed here provide a way forward in understanding daily ecosystem function or filling large gaps in eddy covariance data at a significant number of Fluxnet sites.

54 ENVIRONMENTAL SCIENCES↗

Quantifying carbon budget, crop yields and their responses to environmental variability using the ecosys model for U.S. Midwestern agroecosystems

As one of the major agricultural production areas in the world, the United States (U.S.) Midwest plays a vital role in the global food supply and agricultural ecosystem services. Although significant efforts have been made in modeling the carbon cycle dynamics over this area, large uncertainty still exists in the previous simulations in terms of reproducing individual components of the carbon cycle and their responses to environmental variability. Here we evaluated the performance of an advanced agroecosystem model, ecosys, in simulating carbon budgets over the U.S. Midwest, considering both the magnitude of carbon flux/yield and its response to environmental (climate and soil) variability. We conducted model simulations and evaluations at 7 cropland eddy-covariance sites as well as over 293 counties of Illinois, Indiana, and Iowa in the U.S. Midwest. The site-level simulations showed that ecosys captured both the magnitude and seasonal patterns of carbon fluxes (i.e., net ecosystem carbon exchange (NEE), ecosystem gross primary production (GPP), and ecosystem respiration (Reco)), leaf area index (LAI), and dynamic plant carbon allocation processes, with R2 equal to 0.92, 0.87, 0.87, and 0.78 for GPP, NEE, Reco, and LAI, respectively across all the sites compared with the observations. For regional scale simulations, ecosys reproduced the spatial distribution and interannual variability of corn and soybean yields with the constraints of observed yields and a new remotely sensed GPP product, with R2 of multi-year averaged simulated and observed yield equal 0.83 and 0.80 for corn and soybean, respectively. The simulated responses of carbon cycle dynamics to environmental variability were consistent with that from the empirical observations at both site and regional scales. Furthermore, our results demonstrated the applicability of ecosys in simulating the carbon cycle and soil carbon dynamics of the U.S. Midwestern agroecosystems under different climate and soil conditions.

54 ENVIRONMENTAL SCIENCES↗

Aboveground and belowground contributions to ecosystem respiration in a temperate deciduous forest

In this study, we developed a three-way carbon dioxide (CO 2 ) flux-partitioning algorithm that separates net ecosystem exchange (NEE) into aboveground plant respiration (R above ), belowground root and soil respiration (R below ), and gross primary production (GPP). We applied this algorithm to a coupled dataset of continuous chamber-measured soil respiration and eddy covariance (EC)-measured NEE of CO 2 in an oak-hickory (Quercus-Carya) deciduous broadleaf forest from 2006 to 2015. We found that on annual time scale, R below dominated over R above with the former accounting for 66.9–86.4% and the latter 13.6–33.1%, of the total ecosystem respiration (R eco ). The ratio of R below to R above varied seasonally, ranging from 1.77 to 7.25 in growing season, and 1.02 to 4.57 in non-growing season. The temperature sensitivity (E 0 ) of R below was significantly higher than that of R above , and E 0 of R eco responded differently to air and soil temperature. Over the whole study period, annual mean R above , R below , and GPP were 243, 806, and 1170 g C m –2 , respectively, with annual R eco accounting for 89.6% of GPP, of which 68.8% was lost as R below and 20.8% lost as R above , and leaving only 10% of the carbon fixation in ecosystems. Furthermore, these estimates, however, did not consider potential light inhibition of leaf respiration. If we accept the presence of light inhibition, then the daytime three-way partitioning method would underestimate annual R above by 20.4% whereas the nighttime method would overestimate R above by 23.9% and GPP by 4.7%, compared with estimates accounting for light inhibition in leaves.

54 ENVIRONMENTAL SCIENCES↗

Drought-induced decoupling between carbon uptake and tree growth impacts forest carbon turnover time

The ability of forests to withstand, and recover from, acute drought stress is a critical uncertainty regarding the impacts of climate change on the terrestrial carbon (C) cycle, but it is unclear how drought responses scale from individual trees to whole forests. Here, we assembled a dataset of tree-ring chronologies co-located within the footprint of eddy covariance towers across North America and Europe, with the aim of quantifying the sensitivity of tree radial growth versus gross primary productivity (GPP) during and following drought. We found that drought induced a large decoupling across C cycle processes, whereby GPP was relatively resistant to water stress despite large reductions in tree-ring widths. This decoupling also occurred in the year following drought (i.e., a ‘drought legacy effect’), and was similar in magnitude in response to both summer and winter droughts. By modeling whole-forest C turnover time, we show that a radial growth-GPP decoupling has important ramifications for the forest C cycle, especially if the C not used to support radial growth is instead allocated towards pools with short residence times. Furthermore, our results demonstrate that quantifications of drought impacts that rely solely on C uptake are missing this fundamental pathway through which drought alters the forest C cycle and the resulting feedbacks to the climate system.

54 ENVIRONMENTAL SCIENCES↗

Dryness controls temperature-optimized gross primary productivity across vegetation types

Temperature response of gross primary productivity (GPP) is a well-known property of ecosystem, but GPP at the optimum temperature (GPP_T opt ) has not been fully discussed. Our understanding of how GPP_T opt responds to warming and water availability is highly limited. Here, in this study, we analyzed data at 326 globally distributed eddy covariance sites (79°N-37°S), to identify controlling factors of GPP_T opt . Although GPP_T opt was significantly influenced by soil moisture, global solar radiation, mean annual temperature, and vapor pressure deficit in a non-linear pattern (R 2 = 0.47), the direction and magnitude of these climate variables’ effects on GPP_T opt depend on the dryness index (DI), a ratio of potential evapotranspiration to precipitation. The spatial pattern showed that soil moisture did not affect GPP_T opt across energy-limited sites with DI < 1 while dominated GPP_T opt across water-limited sites with DI >1. The temporal pattern showed that GPP_T opt was lowered by warming or low precipitation in water-limited sites while energy-limited sites tended to maintain a stable GPP_T opt regardless of changes in air temperature. Vegetation types in humid climates tended to have higher GPP_T opt and were more likely to benefit from a warmer climate since it was not restricted by water conditions. This study highlights that the response of GPP_T opt to global warming depends on the dryness conditions, which explains the nonlinear control of water and temperature over GPP_T opt . Our finding is essential to realistic prediction of terrestrial carbon uptake under future climate and vegetation conditions.

54 ENVIRONMENTAL SCIENCES↗

Fine resolution remote sensing spectra improves estimates of gross primary production of croplands

Gross primary production (GPP) is a fundamental measure of the terrestrial carbon cycle critical to our understanding of ecosystem function under the changing climate and land use. Remote sensing enables access to continuous spatial coverage, but remains challenged in heterogeneous croplands. Coarse resolution products, like MOD17A (500 m), may aggregate fragmented land cover types commonly found in heavily managed landscapes and misrepresent their respective contribution to carbon production. Consequently, this study demonstrates the capability of fine-resolution imagery (20-30 m) and available red-edge vegetation indices to characterize GPP across seven Midwest cropping systems. Four sites were established on a 22-year-old USDA Conservation Reserve Program (CRP); and the other three on land conventionally farmed with corn-soybean-wheat rotation (AGR). We compare in situ GPP estimates from eddy-covariance towers with ten satellite models: eight variants of the vegetation photosynthesis models (VPM), of which five include a red-edge vegetation index, as well as conventional products Landsat CONUS GPP (30 m) and MOD17A2H V6 (500 m). Daily and cumulative fine-resolution imagery integrated within VPM generally agreed with tower-based GPP in heterogeneous landscapes more than those from MODIS 500 m VPM or conventional GPP products from MOD17AH V6 or Landsat 8 CONUS. Here, replacing EVI2 with red-edge indices NDRE2, NDRE1, and MTCI in Sentinel 2 VPMs notably improved explanation of variance and estimation of cumulative GPP. While existing methods using MODIS- and Landsat-derived GPP are important baselines for regional and global studies, future research may benefit from the higher spatial, temporal, and radiometric resolution.

54 ENVIRONMENTAL SCIENCES↗