Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational”

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 1,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Assessing the sensitivity and repeatability of permanganate oxidizable carbon as a soil health metric: An interlab comparison across soils

Soil organic matter is central to the soil health framework. Therefore, reliable indicators of Soil organic matter is central to the soil health framework. Therefore, reliable indicators of changes in soil organic matter are essential to inform land management decisions. Permanganate oxidizable carbon (PDXC), an emerging soil health indicator, has shown promise for being sensitive to soil management. However, strict standardization is required for widespread implementation in research and commercial contexts. Here, we used 36 soils-three from each of the 12 USDA soil orders-to determine the effects of sieve size and soil mass of analysis on PDXC results. Using replicated measurements across 12 labs in the US and the EU (n = 7951 samples), we quantified the relative importance of 1) variation between labs, 2) variation within labs, 3) effect soil mass, and 4) effect of soil sieve size on the repeatability of PDXC. We found a wide range of overall variability in PDXC values across labs (0.03 to 171.8%; mean = 13.4%), and much of this variability was attributable to within-lab variation (median = 6.5%) independently of soil mass or sieve size. Greater soil mass (2.5 g) decreased absolute PDXC values by a mean of 177 mg kg -1 soil and decreased analytical variability by 6.5%. For soils with organic carbon (SOC) >10%, greater soil mass (2.5 g) resulted in more frequent PDXC values above the limit of detection whereas the lower soil mass (0.75 g) resulted in PDXC values below the limit of detection for SOC contents <5%. A finer sieve size increased absolute values of PDXC by 124 mg kg -1 while decreasing the analytical variability by 1.8%. In general, soils with greater SOC contents had lower analytical variability. These results point to potential standardizations of the PDXC protocol that can decrease the variability of the metric. We recommend that the PDXC protocol be standardized to use 2.5 g for soils <10% SOC. Sieve size was a relatively small contributor to analytical variability and therefore we recommend that this decision be tailored to the study purpose. Tradeoffs associated with these standardizations can be mitigated, ultimately providing guidance on how to standardize PDXC for routine analysis.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Benchmarking CO₂ storage simulations: Results from the 11 th Society of Petroleum Engineers Comparative Solution Project

The 11 th Society of Petroleum Engineers Comparative Solution Project (shortened SPE11 herein) benchmarked simulation tools for geological carbon dioxide (CO 2 ) storage. A total of 45 groups from leading research institutions and industry across the globe signed up to participate, with 18 ultimately contributing valid results that were included in the comparative study reported here. This paper summarizes the SPE11 results. A comprehensive introduction and qualitative discussion of the submitted data are provided, together with an overview of online resources for accessing the full depth of data. A global metric for analyzing the relative distance between submissions is proposed and used to conduct a quantitative analysis of the submissions. This analysis attempts to statistically resolve the key aspects influencing the variability between submissions. The study shows that the major qualitative variation between the submitted results is related to thermal effects, dissolution-driven convective mixing, and resolution of facies discontinuities. Moreover, a strong dependence on grid resolution is observed across all three versions of the SPE11. However, our quantitative analysis suggests that the observed variations are predominantly influenced by factors not documented in the technical responses provided by the participants. We therefore identify that unreported variations due to human choices within the process of setting up, conducting, and reporting on the simulations underlying each SPE11 submission are at least as impactful as the computational choices reported.

Nordbotten, Jan M. [Univ. of Bergen (Norway); Norw↗
DOE OSTI2022

Development towards stable chlorine isotope measurements of astromaterials using the modified Middleton source of an accelerator mass spectrometer

The volatile element Cl can be lost during the formation and evolution of planetary bodies, leading to fractionation of its two stable isotopes 35 Cl and 37 Cl. Chlorine isotope variations (reported as δ 37 Cl in parts per thousand (‰) relative to Standard Mean Ocean Chloride, SMOC) are documented to exceed 80‰ between different lunar rock samples and have been variably interpreted as the fingerprint of degassing during accretion, magma ocean, or volcanic portions of lunar history. The large intersample and intrasample variations observed by both bulk isotope ratio mass spectrometry (IRMS) and in-situ secondary ion mass spectrometry (SIMS) methods are difficult to interpret in part because of a paucity of bulk Cl isotope measurements. This lack of high-precision bulk data is due to the relative rarity of IRMS laboratories capable of making these high precision measurements on small samples of precious planetary materials such as those returned by human or robotic exploration. Here we present a new method for performing high precision δ 37 Cl measurements using the modified Middleton ion source of an existing accelerator mass spectrometer. For samples with as little as 1 μg Cl–the equivalent of 2–4 mm 3 of a typical lunar rock sample–the average cathode accuracy is ~1‰. Finally, cathode reproducibility is typically ~1‰ (2σ) for samples with at least 10 μg of Cl, increasing to ~3–6‰ for aliquots with ~1–2 μg Cl, similar to published SIMS results and sufficient to study astromaterials from the Moon, Mars, or 4 Vesta, which have tens of ‰ observed variations.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2021

Multi-cyclone analysis and machine learning model implications of cyclone effects on forests

Past studies of cyclones (hurricanes, typhoons, tropical cyclones) disturbance showed that meteorological, topographical, and biological factors affect the patterns of forest disturbance intensity but left open the extent to which these findings were representative across different global cyclone regions. Using remote sensing data and machine learning models, we examined how these factors change over spatial scales and assessed their consistency across four major cyclones: Katrina (August 2005), Rita (September 2005), Yasi (February 2011), and María (September 2017). Our results revealed that the factors which best explained forest disturbance intensity pattern varied across these regions. Wind speed and precipitation were the dominant factors contributing to the variation in impacts of Katrina; terrain features, especially elevation, explained most of the variation in disturbance intensity of Rita; pre-disturbance vegetation condition was significant predictors of effects of Yasi; these factors played equal roles in explaining the disturbance intensity variation of María. A 40 m/s (144 km/h) wind speed threshold was proposed to split low- and high-level forest disturbance intensity. Other than wind speed, few generalizations can be made on features across multiple regions. We built several generalized hurricane impact models, which worked well with the test data from cyclones used for model development (R 2 = 0.89). However, these models did not have good predictions on other cyclones, such as Michael (October 2018) and Laura (August 2020). This study showed that each cyclone interacted with the landscape in a unique way and the challenges remained in building a generalized cyclone impact model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Spatial and temporal variability of tritium in precipitation within South Africa and it's bearing on hydrological studies

Tritium, the radioactive isotope of hydrogen, has been used to understand groundwater recharge processes for decades. The current variation of tritium in the atmosphere is largely attributed to stratospheric production and fall out rates as well as global circulation phenomena controlling the hydrological cycle. Global controls on the variability in atmospheric tritium activity are poorly suited to explain local variation and tritium activities in precipitation are often assumed to be uniform over both local and regional catchments and watersheds. This assumption can result in both over and under estimation of modern recharge within an aquifer when using tritium as the recharge proxy. In order to minimize the inherent prediction residuals associated with tritium based recharge investigations, the variability of tritium in precipitation was modelled from 127 spatial precipitation samples taken over a two year period, combined with a 76 precipitation sample group-set taken over a one year period in a single location. Precipitation events were traced backward in time, from the point of collection, using HYSPLIT modelling to ascertain the origins of moisture content as well as the altitudes of moisture origin reached along the particle track. Tritium activities, collected over a one year period in Paarl, range from 0.45 to 4.16 TU and have a mean of 1.59 TU. Spatial storm events in the Western Cape in 2017 and 2018 had a range from 0 to 2.2 and 0.37 to 3.27 TU, respectively, with mean activities of 1.18 (n = 34) and 1.25 TU (n = 32). Both storm events had similar tritium variability (σ = 0.5 n = 35 and 0.48 n = 32). Regional precipitation events had the largest range of tritium activities (0.55–12.2 TU). Although not all tritium activities can be explained by interrogating the water mass origin, this study suggests that approximately 90% of events can be completely or partially attributed to the origin of the water mass. The variability of tritium, both spatially and temporally, was higher than expected, confirming that when uniform tritium inputs are used, the groundwater system would provide inaccurate modern recharge estimates. Higher spatial resolution of tritium variation in precipitation for a particular region will improve our ability to relate tritium activities in groundwater to local precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A combined experimental and crystal plasticity study of grain size effects in magnesium alloys

This work presents a method to incorporate the micro Hall-Petch equation into the crystal plasticity finite element (CPFE) framework accounting for the microstructural features to understand the coupling between grain size, texture, and loading direction in magnesium alloys. The effect of grain size and texture is accounted for by modifying the slip resistances of individual basal and prismatic systems based on the micro Hall-Petch equation. The modification based on the micro Hall-Petch equation endows every slip system at each microstructural point with a slip system-level grain size and maximum compatibility factor, which are in turn used to modify the slip resistance. While the slip-system level grain size is a measure of the grain size, the maximum compatibility factor encodes the effect of the grain boundary on the slip system resistance modification and is computed based on the Luster-Morris factor. The model is calibrated using experimental stress-strain curves of Mg-4Al samples with three different grain sizes from which the Hall-Petch coefficients are extracted and compared with Hall-Petch coefficients predicted using original parameters from previous work. The predictability of the model is then evaluated for a Mg-4Al sample with different texture and three grain sizes subjected to loading in different directions. The calibrated parameters are then used for some parametric studies to investigate the variation of Hall-Petch slope for different degrees of simulated spread in basal texture, variation of Hall-Petch slope with loading direction relative to basal poles for a microstructure with strong basal texture, and variation of yield strength with change in grain morphology. The proposed approach to incorporate the micro Hall-Petch equation into the CPFE framework provides a foundation to quantitatively model more complicated scenarios of coupling between grain size, texture and loading direction in the plasticity of Mg alloys.

Crystal plasticity↗
DOE OSTI2023

Long–short-term memory encoder–decoder with regularized hidden dynamics for fault detection in industrial processes

The ability of recurrent neural networks (RNN) to model nonlinear dynamics of high dimensional process data has enabled data-driven RNN-based fault detection algorithms. Previous studies have focused on detecting faults by identifying the discrepancies in data distribution between the faulty and normal data, as reflected in prediction errors generated by RNN models. However, in industrial processes, variations in data distribution can also result from changes in normal control setpoints and compensatory control adjustments in response to disturbances, making it hard to differentiate between normal and faulty conditions. This paper proposes a fault detection method utilizing a long short-term memory (LSTM) encoder–decoder structure with regularized hidden dynamics and reversible instance normalization (RevIN) to compactly represent high-dimensional measurements for effective monitoring. During training, the hidden states of the model are regularized to form a low-dimensional latent space representation of the original multivariate time series data. As a result, the prediction errors of the latent states can be used to monitor the abnormal dynamic variations, while the reconstruction errors of the measured variables are used to monitor the abnormal static variations. Furthermore, the proposed indices can reflect operating conditions, even when the distribution of test data changes, which helps distinguish faults from normal adjustments and disturbances that controllers can settle. Here, data from numerical simulation and the Tennessee Eastman process are used to illustrate the effectiveness of the proposed fault detection method.

42 ENGINEERING↗
DOE OSTI2025

Apatite geochemistry as a tool for understanding the petrogenesis of layered mafic-ultramafic rocks in the Bushveld Complex, South Africa

The sources of the magmas that formed the Rustenburg Layered Suite of the Bushveld Complex in South Africa remain debated, despite decades of research. Vertical and lateral variation in bulk rock and mineral separate Sr-Nd isotopic compositions, which generally indicate enriched sources, demonstrate that the layered sequence was formed by the emplacement of multiple batches of magma, crucially resulting in episodes of PGE-Cr-V mineralisation. The Lu-Hf isotope compositions of zircon are, however, at odds with the bulk rock Sr-Nd isotopic heterogeneity as they show near homogeneous compositions throughout the layered sequence (εHf (2.06 Ga) =−8). This lack of variation in Hf isotope composition has been attributed to deep, continental lithospheric mantle-related and/or crustal contamination of plume-derived Bushveld magmas. In this study, we analysed the major, trace element and Sr-Nd isotope geochemistry of apatite in the Rustenburg Layered Suite. Apatite occurs as an intercumulus mineral in the lowermost regions and a cumulus mineral in the uppermost regions of the layered sequence and can therefore be used to test existing models for the isotopic disequilibrium between bulk rock Sr-Nd and zircon Hf isotopic compositions. Apatite is largely chlorapatite in the lowermost regions and fluorapatite in the uppermost regions of the layered sequence. The Merensky Reef is unusual in that it comprises both chlorapatite and fluorapatite. Apatite throughout the layered sequence is generally unzoned and shows no evidence of late-stage alteration. Trace element data show that apatite is enriched in L/HREE, with common negative Eu-Sr anomalies. These trace element signatures are consistent with a magmatic origin for the apatite grains, with prior, or concurrent, plagioclase crystallization from the same melt. Variability in in situ Sr and Nd isotope compositions of apatite is recorded throughout the layered sequence with εNd (2.06 Ga) compositions varying between −2.5 and − 10.2 and initial 87 Sr/ 86 Sr compositions varying between 0.7079 and 0.7103 (for the Marikana dikes only). The variability in Sr-Nd isotope compositions of apatite is consistent with the bulk rock (and mineral separate) variation in Sr-Nd isotope compositions, suggesting apatite preserves primary magmatic compositions in the Rustenburg Layered Suite.

Apatite↗
DOE OSTI2022

Closed-loop control of meltpool temperature in directed energy deposition

The objective of this work is to mitigate flaw formation in powder and laser-based directed energy deposition (DED) additive manufacturing process through close-loop control of the meltpool temperature. In this work, the meltpool temperature was controlled by modulating the laser power based on feedback signals from a coaxial two-wavelength imaging pyrometer. The utility of closed-loop control in DED is demonstrated in the context of practically inspired trapezoid-shaped stainlesssteel parts (SS 316L). We demonstrate that parts built under closed-loop control have reduced variation in porosity and uniform microstructure compared to parts built under open-loop conditions. For example, post-process characterization showed that closed-loop processed parts had a volume percent porosity ranging from 0.036% to 0.043%. In comparison, open-loop processed parts had a larger variation in volume percent porosity ranging from 0.032% to 0.068%. Further, parts built with closed-loop processing depicted consistent dendritic microstructure. By contrast, parts built with open-loop processing showed microstructure heterogeneity with the presence of both dendritic and planar grains, which in turn translated to large variation in microhardness

36 MATERIALS SCIENCE↗
DOE OSTI2026

Nondestructive ultrasonic characterization of a triple-weld-bead wire arc additively manufactured ER70S-6 S-curved wall

Understanding build-scale microstructural variation in wire arc additive manufacturing (WAAM) of low-carbon steels is essential for ensuring consistent structure–property relationships throughout large components. Conventional destructive characterization techniques, such as scanning electron microscopy (SEM) and electron backscatter diffraction (EBSD), are time-intensive and limited to localized regions, making comprehensive evaluation of large WAAM structures challenging. In this study, a nondestructive ultrasonic approach was employed to characterize a 252 mm tall ER70S-6 S-curved WAAM wall produced using a triple-bead deposition strategy. Optical and SEM analysis revealed a repeating dual-region microstructure consisting of uniform polygonal ferrite at melt pool centers and heterogeneous ferrite with coarse and fine grains near melt pool boundaries, attributed to cyclic thermal conditions. Longitudinal ultrasonic backscatter imaging was used to evaluate the continuity of this periodicity along the full build height. The ultrasonic response exhibited a consistent repeating pattern that correlated with the observed layer-wise microstructural variation. X-ray computed tomography confirmed the absence of detectable porosity, indicating that ultrasonic contrast is primarily governed by grain morphology. Overall, the results indicate that longitudinal backscatter ultrasound is a promising nondestructive characterization technique for validating microstructural variations along the s-curved WAAM wall, with significant potential for microstructure optimization and process control.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Algorithms for coordinate reconstruction in position-sensitive virtual Frisch-grid detectors

Arrays of position-sensitive virtual Frisch-grid (VFG) CdZnTe (CZT) detectors provide a cost-effective solution for integrating large-area arrays for gamma-ray imaging and spectroscopy. These detectors employ high-aspect ratio CZT crystals (bars) with thicknesses up to 50 mm and cross-sections of up to 10 × 10 mm 2 . Despite the long drift distances of charge carriers in such crystals, the detectors have demonstrated excellent performance, achieving energy resolutions better than 1 % full width at half maximum (FWHM) at 662 keV and 3D position resolutions finer than 1 mm. The high spatial resolution is a critical feature of these detectors, as it enables correction of response non-uniformities caused by crystal defects, which remain present even in the highest-quality CZT material. Dislocations and dislocation walls are the primary defects responsible for variations in charge carrier losses as they drift from the interaction points toward the charge-collecting electrodes. The mechanism by which these defects affect carrier transport is generally well understood. Dislocations and sub-grain boundaries act as sinks for carrier-trapping centers, primarily impurities and secondary phases such as tellurium inclusions and precipitates. Here, these regions exhibit significantly higher carrier-trapping rates, leading to variations in the μτ-products. Because the locations of these micron-sized regions are fixed within the detector volume, fluctuations in the total collected charge arise from the random distribution of interaction sites. This results in non-uniform detector responses and degradation of energy resolution. However, by measuring the interaction-site locations with sufficient precision, charge-loss variations can be accurately corrected, allowing recovery of nearly intrinsic energy resolution.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Sensitivity and uncertainty of the IFR-1 BISON benchmark

The fuel performance code BISON is being used to evaluate metallic fuel for a new fast-spectrum test reactor called the Versatile Test Reactor, which is being considered for adoption by the US Department of Energy. To quantify the accuracy of BISON predictions, researchers at Oak Ridge National Laboratory have been developing a series of benchmarks based on legacy metallic fuel experiments. As part of this effort, the sensitivity of BISON predictions to variations in model inputs and the uncertainties associated with BISON predictions must be established. This paper summarizes efforts to perform a comprehensive sensitivity analysis (SA) and uncertainty quantification (UQ) on a benchmark based on the IFR-1 experiment.For the SA, at least one input was chosen from every BISON model and physics module used in the benchmark. The inputs were varied individually in a series of BISON simulations. Here, the resulting variations in benchmark predictions were normalized to calculate sensitivities. These sensitivities were then used to inform input selections for the UQ.The UQ was performed using the Monte Carlo UQ method. A literature review was conducted to estimate uncertainty distributions for the selected inputs, and values were sampled randomly from each distribution in a series of BISON simulations. Variations in the benchmark predictions were used to estimate uncertainty distributions and confidence intervals. It was found that nearly 100% of the benchmark predictions matched the corresponding legacy values within the confidence intervals. However, this is at least partially because the confidence intervals associated with benchmark predictions were wide. The uncertainty contributions of assumptions in the benchmark, experimental uncertainties, and BISON models were quantified. Some analysis was performed to identify inputs that contributed to the uncertainties. Finally, recommendations are made for future benchmark and future BISON development.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Effect of 316 stainless steel powder processing conditions on microstructure

One parameter that has both economic and performance effects in the manufacturing of steel parts is the powder feedstock used. 316 L steel powder can be produced through water atomization or gas atomization, where both the powder size and the atomization process determine the cost of powder feedstock. Gas atomized powder is spherical and has a lower oxygen content than water atomized powder, making it the preferred choice for performance, but comes at a higher financial cost. Here, in this study, various 316 L powder batches are characterized to understand physical properties and microstructural variations. Density and flowability were determined using gas pycnometry and Hall flow tests. Powder morphology and porosity were qualitatively analyzed via optical microscopy and scanning electron microscopy, Finally, transmission electron microscopy was used to determine changes in phases present. Spherical gas atomized powder of 30–100 μm contained Mn$-$Si rich oxides both on the surface and within the matrix. Meanwhile, irregular water atomized powder 30 μm in diameter contained silicon oxides without manganese. All powder particles contained a nickel rich cellular structure which was shown to support the formation of additional phases. While literature has observed general microstructural features across 316 L powder particles, transmission electron microscopy (TEM) in this study has identified variations in precipitation and cell structures. Information regarding morphology, flowability, density, and phase distributions can be applied to understand variations in part consolidation, which will in turn show which powder aspects are significant for production.

Chemical segregation↗
DOE OSTI2021

Unpacking the drivers of diurnal dynamics of sun-induced chlorophyll fluorescence (SIF): Canopy structure, plant physiology, instrument configuration and retrieval methods

Sun-induced chlorophyll fluorescence (SIF) from spaceborne sensors is a promising tool for global carbon cycle monitoring, but its application is constrained by insufficient understanding of the drivers underlying diurnal SIF dynamics. SIF measurements from ground-based towers can reveal diurnal SIF dynamics across biomes and environmental conditions; however, meaningful interpretation of diurnal variations requires disentangling impacts from canopy structure, plant physiology, instrument configuration and retrieval methods, which often interact with and confound each other. This study aims to unpack these drivers using 1) concurrent ground and airborne canopy-scale and leaf-scale measurements at a corn field, 2) a mechanistic SIF model that explicitly considers the dynamics of photochemistry (via the fraction of open photosystem II reaction centers, qL) and photoprotection (via nonphotochemical quenching, NPQ) as well as their interactive dependence on the sub-canopy light environment, and 3) cross-comparison of SIF instrument configurations and retrieval methods. We found that crop row orientations and sun angles can introduce a distinctive midday dip in SIF in absence of stress, due to a midday drop of absorbed photosynthetically active radiation (APAR) when crop rows are north-south oriented. Canopy structure caused distinctive responses in both qL and NPQ at different positions within the vertical canopy that collectively influenced fluorescence quantum yield (Φ F ) at the leaf scale. Once integrated at the canopy scale, diurnal dynamics of both APAR and canopy escape probability (ε) are critical for accurately shaping diurnal SIF variations. While leaf-level qL and NPQ exhibited strong diurnal dynamics, their influence was attenuated at the canopy scale due to opposing effects on SIF at different canopy layers. Furthermore, different system configurations (i.e., bi-hemispherical vs. hemispherical-conical) and retrieval methods can bias the SIF magnitude and distort its diurnal shapes, therefore confounding the interpretation of inherent strength and dynamics of SIF emission. Our findings demonstrate the importance of crop row structures, interactive variations in canopy structure and plant physiology, instrument configuration, and retrieval method in shaping the measured dynamics of diurnal SIF. This study highlights the necessity to account for these factors to accurately interpret satellite SIF, and informs future synthesis work with different SIF instrumentation and retrieval methods across sites.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Decoupling of greenness and gross primary productivity as aridity decreases

Ecosystem primary productivity is a key ecological process influencing many ecosystem services, including carbon storage. Thus, clarifying how primary productivity in terrestrial ecosystems responds to climatic variability can reveal key mechanisms that will drive future changes in the global carbon budget. Satellite products of canopy greenness are widely used as proxies for vegetation productivity to evaluate how ecosystems respond to climate variability. However, to what degree inter-annual variations in productivity are consistent with greenness and how this relationship varies spatially remains unclear. Here we investigated the strength of the coupling between inter-annual variations in leaf area index (LAI, a measure of greenness) and ecosystem gross primary productivity (GPP) derived from eddy covariance towers, i.e., the r2 of the LAI-GPP relationship. Overall, inter-annual GPP and LAI were highly coupled (i.e., high r2) in arid grasslands, but were fully decoupled in mesic evergreen broadleaf forests, indicating that this relationship varies strongly along aridity gradients. A possible mechanism of the spatial variation in the LAI-GPP relationship is that the tradeoff between ecosystem structure (LAI) and physiology (photosynthesis per unit leaf area) becomes stronger in more humid climates. Land models overestimated the r2 of LAI-GPP correlation for most ecosystem types and failed to capture the spatial pattern along aridity gradients. We conclude that relying on greenness products for evaluating inter-annual changes in vegetation productivity may bias assessments, especially in tropical rainforest ecosystems. Our findings may also reconcile observed disparities between responses in greenness and GPP during drought in Amazon forests.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A physical method for downscaling land surface temperatures using surface energy balance theory

Fine-resolution land surface temperature (LST) derived from thermal infrared remote sensing images is a good indicator of surface water status and plays an essential role in the exchange of energy and water between land and atmosphere. A physical surface energy balance (SEB)-based LST downscaling method (DTsEB) is developed to downscale coarse remotely sensed thermal infrared LST products with fine-resolution visible and near-infrared data. Here, the DTsEB method is advantageous for its ability to mechanically interrelate surface variables contributing to the spatial variation of LST, to quantitatively weigh the contributions of each related variable within a physical framework, and to efficaciously avoid the subjective selection of scaling factors and the establishment of statistical regression relationships. The applicability of the DTsEB method was tested by downscaling 12 scenes of 990 m Moderate Resolution Imaging Spectroradiometer (MODIS) and aggregated Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) LST products to 90 m resolution at six overpass times between 2005 and 2015 over three 9.9 km by 9.9 km cropland (mixed by grass, tree, and built-up land) study areas. Three typical LST downscaling methods, namely the widely applied TsHARP, the later developed least median square regression downscaling (LMS) and the geographically weighted regression (GWR), were introduced for intercomparison. The results showed that the DTsEB method could more effectively reconstruct the subpixel spatial variations in LST within the coarse-resolution pixels and achieve a better downscaling accuracy than the TsHARP, LMS and GWR methods. The DTsEB method yielded, on average, root mean square errors (RMSEs) of 2.01 K and 1.42 K when applied to the MODIS datasets and aggregated ASTER datasets, respectively, which were lower than those obtained with the TsHARP method, with average RMSEs of 2.41 K and 1.71 K, the LMS method, with average RMSEs of 2.35 K and 1.63 K, and the GWR method, with average RMSEs of 2.38 K and 1.64 K, respectively. The contributions of the related surface variables to the subpixel spatial variation in the LST varied both spatially and temporally and were different from each other. In summary, the DTsEB method was demonstrated to outperform the TsHARP, LMS, and GWR methods and could be used as a good alternative for downscaling LST products from coarse to fine resolution with high robustness and accuracy.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Field-scale dynamics of planting dates in the US Corn Belt from 2000 to 2020

Crop planting dates are a dynamic feature of agricultural systems that respond to short- and long-term climate signals, crop and cultivar selection, and technology changes. Planting date records are essential for yield gap analyses, accurate crop modeling, and tracking farmer adaptations to weather and climate change. Although planting dates have high variation at local scales due to heterogeneity in farm resources and decision-making, available long-term data on planting dates is largely restricted to aggregated regional statistics or, at best, satellite-derived datasets with limited spatiotemporal extent and at resolutions unable to distinguish individual fields (> 250 m). Here, we generated retrospective annual field-scale (30 m) planting date maps for both maize and soybeans spanning 2000-2020 across a 12 state region in the United States Corn Belt based on Landsat satellite data and a large ground sample of over 28,000 maize and soybean fields. Using training data from 2015-2020 for model selection, we found that planting date predictions improved with harmonic regression of Landsat data and additional annual weather covariates. The preferred random forests model approximately doubled performance compared to a null model based on state median planting dates, capturing 47% of field-level variation for maize (mean absolute error, MAE = 7.4 days) and 44% for soybeans (MAE = 7.5 days) against held-out ground truth test data for 2008-2014. We also evaluated the full 2000-2020 dataset with state agricultural statistics, finding strong agreement with median planting dates for maize (R 2 = 0.76, MAE = 4.4 days) and slightly lower agreement for soybeans (R 2 = 0.65, MAE = 5.4 days) when aggregated to the state level. We then used this new dataset to analyze environmental determinants of planting dates at a finer-scale than previously possible, controlling for unobserved variation at the sub-state district level. We found that during 2000-2020, each standard deviation increase in rainfall delayed planting by ~ 2.5 days, and fields with higher soil productivity ratings tended to be planted earlier. We did not find meaningful trends over the last two decades in planting dates for maize or soybeans, in contrast to trends towards earlier planting dates late last century and predicted for this period in climate adaptation studies. We hypothesize increases in early season rainfall may have inhibited these shifts towards earlier planting. Remotely sensed planting dates will be a useful tool for yield gap analyses, crop simulation modeling, and ongoing assessment of climate adaptation.

54 ENVIRONMENTAL SCIENCES↗