Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mapping variables”

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 91 records · Page 5

AI-Batt (Autonomous Identification of Battery Life Models) [SWR 21-36]

Autonomous Identification of Battery Life Models (AI-Batt) AI-Batt is a MATLAB code base for developing lifetime models for batteries from accelerated aging data. The code base provides many functions for processing, visualizing, and modeling battery aging data, making the data processing, exploration, and modeling workflow substantially faster. These tools are tailored for working with battery aging data sets, which usually consist of many separate time-series for each cell, with many test conditions and possible replicates at each condition, which makes it difficult to simply process or visualize the data set. Complex modeling tasks, such as cross-validation, sensitivity analysis, and uncertainty quantification have been implemented to enable thorough statistical investigation of model predictions. Additionally, several machine-learning algorithms are implemented to autonomously identify suitable models via symbolic regression. Data processing functions automatically cast data from the struct data type, which is commonly used to store experimental data, but is not an acceptable input for most algorithms, to the table data type, which can be easily used as input to any optimization algorithm. Also, the data can be separated into time-invariant and time-variant data tables, which is helpful for exploring the data set as well as developing separate models for time-variant and time-invariant aging mechanisms. For example, in aging tests with constant temperature, temperature is a time-invariant experimental condition. Visualization tools enable plotting of data, model fits, and model simulations possible with single-line function calls, empowering data exploration of complex data sets with both time-varying and time-invariant trends. Plots can be automatically generated for the whole data set, or separated by data group (groups of test replicates) or individual data series. Data points or data series can be automatically colored by the value of a variable with a variety of color maps, and model predictions can also be colored by the value of a fit statistic. Comparisons between data sets and the predictions/simulations of different models on the same data set can be easily plotted as well. Distributions of parameter values from bootstrap resampling can be plotted to visualize the reliability of parameter estimation, or determine any correlations between parameters. Modeling tools handle the complex task of creating and parsing symbolic equations for modeling battery lifetime. Equations are parsed to grab relevant data variables, parameter values, or specified sub-models for input into optimization, evaluation, or simulation functions. Models can be optimized locally (one set of parameters for each data series), bi-level (some parameters shared across the data set), or globally (single set of parameters for all data). Functions implementing symbolic regression algorithms help users to discover effective model equations, even in poorly sampled, high-dimensional data.

Smith, Kandler↗

Regional Feedstock Partnership Biomass Quality Assessment Final Report

The United States (U.S.) Department of Energy (DOE) developed the Billion-Ton Vision to enable production of one-billion tons of sustainable, reliable biomass for the bioenergy industry by 2030 (Perlack et al., 2005). The Sun Grant Regional Feedstock Partnership (RFP) was organized to fill information gaps and validate biomass yield assumptions related to the Billion-Ton Study (Owens, 2018; Owens, Karlen, and Lacey, 2016). Along with the more than 130 scientific publications generated from these studies, yield and sustainability data from the RFP field trials not only validated the Billion-Ton estimates, but were critical in developing both the U.S. Billion-Ton Update report in 2011 and the 2016 Billion-Ton Report (DOE, 2011; 2016). The intention of this biomass quality assessment report is to build on these initial successes from the RFP field trials by focusing on variability in biomass quality data necessary to evaluate conversion performance. This report contains a summary of chemical quality results from samples collected as part of the RFP field trials. This report focuses on assessment of the impact of experimental agronomic designs on biomass properties followed by analyses of the impact of environmental and production variables on biomass properties. Datasets include species and other genetic variables, fertilizer treatments, harvest information, and yield, as well as other publicly available data such as precipitation, temperature, soil properties, and drought. The key outcomes from this chemical quality focused assessment have included: • Complete evaluation of the impacts of agronomic designs, genetics, and environmental conditions on chemical properties for Miscanthus, switchgrass, sorghum, energycane, mixed perennial grasses, and shrub willow short-rotation feedstocks • Over 30 peer review publications and technical reports focused on variability in quality data • Development of spatial and temporal environmental quality prediction maps for Miscanthus and switchgrass feedstocks allowing for comprehensive evaluation of variability in feedstock chemical quality across U.S. regions and over multiple harvest years

09 BIOMASS FUELS↗

Rapid eco‐phenotypic feedback and the temperature response of biomass dynamics

Biomass dynamics capture information on population dynamics and ecosystem-level processes (e.g., changes in production over time). Understanding how rising temperatures associated with global climate change influence biomass dynamics is thus a pressing issue in ecology. The total biomass of a species depends on its density and its average mass. Consequently, disentangling how biomass dynamics responds to increasingly warm and variable temperatures ultimately depends on understanding how temperature influences both density and mass dynamics. Here, we address this issue by keeping track of experimental microbial populations growing to carrying capacity for 15 days at two different temperatures, and in the presence and absence of temperature variability. We develop a simple mathematical expression to partition the contribution of changes in density and mass to changes in biomass and assess how temperature responses in either one influence biomass shifts. Moreover, we use time-series analysis (Convergent Cross Mapping) to address how temperature and temperature variability influence reciprocal effects of density on mass and vice versa. We show that temperature influences biomass through its effects on density and mass dynamics, which have opposite effects on biomass and can offset each other. We also show that temperature variability influences biomass, but that effect is independent of any effects on density or mass dynamics. Last, we show that reciprocal effects of density and mass shift significantly across temperature regimes, suggesting that rapid and environment-dependent eco-phenotypic dynamics underlie biomass responses. Overall, our results connect temperature effects on population and phenotypic dynamics to explain how biomass responds to temperature regimes, thus shedding light on processes at play in cosmopolitan and abundant microbes as the world experiences increasingly warm and variable temperatures.

54 ENVIRONMENTAL SCIENCES↗

Parallel physics-informed neural networks via domain decomposition

Here we develop a distributed framework for the physics-informed neural networks (PINNs) based on two recent extensions, namely conservative PINNs (cPINNs) and extended PINNs (XPINNs), which employ domain decomposition in space and in time-space, respectively. This domain decomposition endows cPINNs and XPINNs with several advantages over the vanilla PINNs, such as parallelization capacity, large representation capacity, efficient hyperparameter tuning, and is particularly effective for multi-scale and multi-physics problems. Here, we present a parallel algorithm for cPINNs and XPINNs constructed with a hybrid programming model described by MPI + X, where X ∈ {CPUs, GPUs}. The main advantage of cPINN and XPINN over the more classical data and model parallel approaches is the flexibility of optimizing all hyperparameters of each neural network separately in each subdomain. We compare the performance of distributed cPINNs and XPINNs for various forward problems, using both weak and strong scalings. Our results indicate that for space domain decomposition, cPINNs are more efficient in terms of communication cost but XPINNs provide greater flexibility as they can also handle time-domain decomposition for any differential equations, and can deal with any arbitrarily shaped complex subdomains. To this end, we also present an application of the parallel XPINN method for solving an inverse diffusion problem with variable conductivity on the United States map, using ten regions as subdomains.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Estimating Submicron Aerosol Mixing State at the Global Scale With Machine Learning and Earth System Modeling

Abstract This study integrates machine learning and particle‐resolved aerosol simulations to develop emulators that predict submicron aerosol mixing state indices from the Earth system model (ESM) simulations. The emulators predict aerosol mixing state using only quantities that are predicted by the ESM, including bulk aerosol species concentrations, which do not by themselves carry mixing state information. We used PartMC‐MOSAIC as the particle‐resolved model and NCAR's CESM as the ESM. We trained emulators for three different mixing state indices for submicron aerosol in terms of chemical species abundance ( χ a ), the mixing of optically absorbing and nonabsorbing species ( χ o ), and the mixing of hygroscopic and nonhygroscopic species ( χ h ). Our global mixing state maps show considerable spatial and seasonal variability unique to each mixing state index. Seasonal averages varied spatially between 13% and 94% for χ a , between 38% and 94% for χ o , and between 20% and 87% for χ h with global annual averages of 67%, 68%, and 56%, respectively. High values in one index can be consistent with low values in another index depending on the grouping of species and their relative abundance, meaning that each mixing state index captures different aspects of the population mixing state. Although a direct validation with observational data has not been possible yet, our results are consistent with mixing state index values derived from ambient observations. This work is a prototypical example of using machine learning emulators to add information to ESM simulations.

54 ENVIRONMENTAL SCIENCES↗

Forward Modeling of 3-D Ion Properties in Jupiter’s Magnetosphere Using Juno/JADE-I Data

The Jovian Auroral Distributions Experiment Ion sensor (JADE-I) on NASA’s Juno mission provides in-situ measurements of ions from 0.1 to 46.2 keV/q inside Jupiter’s magnetosphere. JADE-I is used to study the plasma with two types of datasets from the same measurement: Time-of-flight (TOF) and SPECIES. The TOF dataset provides mass-per-charge measurements with a range of 1–64 amu/q but oversamples particles over 6π steradian viewing per spacecraft spin and has little directional information. On the other hand, the SPECIES dataset can provide a good measurement of the flow direction but does not provide mass-per-charge information due to the telemetry limit. In this study, we developed a 2-step forward modeling method that combines the advantages and avoids the disadvantages of TOF and SPECIES data to derive the 3-D properties of heavy ions. Assuming that the ion velocity distribution can be described with the kappa distribution, we first perform the forward model fit of the TOF data to calculate the relative abundance of heavy ion species. Then we fix the relative abundance and perform the second forward model fit on the SPECIES data. Here, using this method, we obtain the densities of different heavy ions, the shared temperature and kappa value, and the 3-D flow velocity vector. Some data examples of the equatorial plasma disk before Perijove 24 are included to demonstrate the method. Plasma properties can then be mapped to explore spatial and temporal variabilities in Jupiter’s magnetosphere.

79 ASTRONOMY AND ASTROPHYSICS↗

Area-Preserving Mapping of 3D Carotid Ultrasound Images Using Density-Equalizing Reference Map

Carotid atherosclerosis is a focal disease at the bifurcations of the carotid artery. To quantitatively monitor the local changes in the vessel-wall-plus-plaque thickness (VWT) and compare the VWT distributions for different patients or for the same patients at different ultrasound scanning sessions, a mapping technique is required to adjust for the geometric variability of different carotid artery models. In this work, we propose a novel method called density-equalizing reference map (DERM) for mapping 3D carotid surfaces to a standardized 2D carotid template, with an emphasis on preserving the local geometry of the carotid surface by minimizing the local area distortion. The initial map was generated by a previously described arc-length scaling (ALS) mapping method, which projects a 3D carotid surface onto a 2D non-convex L-shaped domain. A smooth and area-preserving flattened map was subsequently constructed by deforming the ALS map using the proposed algorithm that combines the density-equalizing map and the reference map techniques. This combination allows, for the first time, one-to-one mapping from a 3D surface to a standardized non-convex planar domain in an area-preserving manner. Evaluations using 20 carotid surface models show that the proposed method reduced the area distortion of the flattening maps by over 80% as compared to the ALS mapping method.

60 APPLIED LIFE SCIENCES↗

Bayesian reduced-order deep learning surrogate model for dynamic systems described by partial differential equations

We propose a reduced-order deep-learning surrogate model for dynamic systems described by time-dependent partial differential equations. This method employs space–time Karhunen–Loève expansions (KLEs) of the state variables and space-dependent KLEs of space-varying parameters to identify the reduced (latent) dimensions. Subsequently, a deep neural network (DNN) is used to map the parameter latent space to the state variable latent space. An approximate Bayesian method is developed for uncertainty quantification (UQ) in the proposed KL-DNN surrogate model. The KL-DNN method is tested for the linear advection–diffusion and nonlinear diffusion equations, and the Bayesian approach for UQ is compared with the deep ensembling (DE) approach, commonly used for quantifying uncertainty in DNN models. It was found that the approximate Bayesian method provides a more informative distribution of the PDE solutions in terms of the coverage of the reference PDE solutions (the percentage of nodes where the reference solution is within the confidence interval predicted by the UQ methods) and log predictive probability. The DE method is found to underestimate uncertainty and introduce bias. For the nonlinear diffusion equation, we compare the KL-DNN method with the Fourier Neural Operator (FNO) method and find that KL-DNN is 10% more accurate and needs less training time than the FNO method.

97 MATHEMATICS AND COMPUTING↗

An Automated Method for Mapping Giant Kelp Canopy Dynamics from UAV

Satellite and aerial imagery have been used extensively for mapping the abundance and distribution of giant kelp ( Macrocystis pyrifera ) in southern California. There is now great potential for using unoccupied aerial vehicles (UAVs) to map kelp canopy at very high resolutions. However, tides and currents have been shown to affect the amount of floating kelp canopy on the water surface, and the impacts of these processes on remotely sensed kelp estimates in this region have not been fully quantified. UAVs were used to map fine-scale changes in canopy area due to tidal height and current speed at kelp forests off the coast of Palos Verdes, CA and Santa Barbara, CA. An automated method for detecting kelp canopy was developed that was 67% accurate using red-green-blue (RGB) UAV imagery and 93% accurate using multispectral UAV imagery across a range of weather, ocean, and illumination conditions. Increases in tidal height of 1 m reduced the amount of floating kelp canopy by 15% in Santa Barbara and by over 30% in Palos Verdes. The effect of current speed on visible kelp canopy was inconclusive, but there was a trend towards lower canopy area with increased current speed. Therefore, while tidal height and current speed can introduce significant variability to estimates of kelp abundance, the magnitude of this variability is site specific. Still, UAVs are a valuable tool for mapping of kelp canopy and can provide greater spatial resolution and temporal coverage than is possible from many satellite sensors. This data can provide insight into the patterns and drivers of high frequency fluctuations in kelp abundance.

54 ENVIRONMENTAL SCIENCES↗

Machine learning for the identification of phase transitions in interacting agent-based systems: A Desai-Zwanzig example

Deriving closed-form analytical expressions for reduced-order models, and judiciously choosing the closures leading to them, has long been the strategy of choice for studying phase- and noise-induced transitions for agent-based models (ABMs). In this paper, we propose a data-driven framework that pinpoints phase transitions for an ABM—the Desai-Zwanzig model—in its mean-field limit, using a smaller number of variables than traditional closed-form models. To this end, we use the manifold learning algorithm Diffusion Maps to identify a parsimonious set of data-driven latent variables, and we show that they are in one-to-one correspondence with the expected theoretical order parameter of the ABM. We then utilize a deep learning framework to obtain a conformal reparametrization of the data-driven coordinates that facilitates, in our example, the identification of a single parameter-dependent ordinary differential equation (ODE) in these coordinates. Additionally, we identify this ODE through a residual neural network inspired by a numerical integration scheme (forward Euler). We then use the identified ODE—enabled through an odd symmetry transformation—to construct the bifurcation diagram exhibiting the phase transition.

97 MATHEMATICS AND COMPUTING↗

Poisson’s ratio of eTPU molded bead foams in compression via in situ synchrotron X-ray microtomography

In situ synchrotron X-ray microtomography was used to characterize the bulk deformation behavior by computing the Poisson's ratio of expanded thermoplastic polyurethane (eTPU) molded bead foams used in footwear midsole during compression. Quantitative data on morphological characteristics were obtained using an iterative image processing workflow. Image correlation on the 4D datasets using DVC was performed to calculate the volumetric and axial strain to estimate the Poisson ratio. Strain maps from DVC showed the influence of variability in ligament thickness distribution on the global mechanical behavior exhibited which dominated the response seen in these bead foams. Lastly, our results showed a strong correlation between Poisson ratio and distribution of ligament thickness in foams.

36 MATERIALS SCIENCE↗

Revealing causes of macroscale heterogeneity in lithium ion pouch cells via synchrotron X-ray diffraction

Heterogeneous battery performance is a critical issue for maximization of cell lifetime capacity and safety. Using high energy synchrotron X-ray diffraction, the influence of charge rate, voltage limit, uneven stack pressure, and gas generation on the lithium transport properties was quantified in single-layer graphite/LiNi 0.5 Mn 0.3 Co 0.2 O 2 pouch cells. A freshly formatted cell tracked in operando during initial fast charge cycles indicated variable position-dependent performances, while lateral mapping showed a significant fast charge (6C) heterogeneity compared to slow charge (C/2). Pressure effects were non-dominant compared to charge rate. Maps of previously aged and rested cells indicate that lateral heterogeneity slowly equilibrates at rest, but regenerates upon further cycling at fast charge rate. Furthermore, an unformatted cell was mapped at charge and discharge during its first formation cycle to analyze the effect of byproduct gases on the heterogeneous lithium transport. Gas was observed as randomly interspersed “bubbles” which locally hindered lithium intercalation and caused significant heterogeneity. Electrode architectures and charging protocols that promote homogeneous intercalation are critical for predictable high-performance and long-life batteries.

25 ENERGY STORAGE↗

Assessing Uncertainties and Approximations in Solar Heating of the Climate System

In calculating solar radiation, climate models make many simplifications, in part to reduce computational cost and enable climate modeling, and in part from lack of understanding of critical atmospheric information. Whether known errors or unknown errors, the community's concern is how these could impact the modeled climate. The simplifications are well known and most have published studies evaluating them, but with individual studies it is difficult to compare. Here, we collect a wide range of such simplifications in either radiative transfer modeling or atmospheric conditions and assess potential errors within a consistent framework on climate-relevant scales. We build benchmarking capability around a solar heating code (Solar-J) that doubles as a photolysis code for chemistry and can be readily adapted to consider other errors and uncertainties. The broad classes here include: use of broad wavelength bands to integrate over spectral features; scattering approximations that alter phase function and optical depths for clouds and gases; uncertainty in ice-cloud optics; treatment of fractional cloud cover including overlap; and variability of ocean surface albedo. We geographically map the errors in W m -2 using a full climate re-creation for January 2015 from a weather forecasting model. For many approximations assessed here, mean errors are ~2 W m -2 with greater latitudinal biases and are likely to affect a model’s ability to match the current climate state. Combining this work with previous studies, we make priority recommendations for fixing these simplifications based on both the magnitude of error and the ease or computational cost of the fix.

58 GEOSCIENCES↗

Renormalized site density functional theory for models of ion hydration

Development of accurate statistical mechanics models of molecular liquid systems is a problem of great practical and fundamental importance. Site-density functional theory (SDFT) is one of the promising directions in this area, but its success hinges upon the ability to efficiently reconcile co-existence of two distinct intra- and inter-molecular interaction regimes in a molecular liquid. Renormalized formulation of site-density functional theory (RSDFT), that we have recently developed, resolves this problem by introducing an additional potential field variable that decouples two interaction scales and maps the molecular liquid problem onto the effective simple liquid mixture. This work provides a critical assessment of RSDFT for the hydrated ion system - a problem, which historically has always been one of the most difficult cases for SDFT applications. Using two-site model of water, we perform a comprehensive analysis of hydrated alkali metal and halogen ions, including both structural and free energy based characteristics. The results indicate that RSDFT provides a significant improvement over conventional 3D-RISM implementations, and may prove useful in coarse grained simulations based on two-site solvent models.

Chuev, Gennady N.↗

Phase-resolving spin-wave microscopy using infrared strobe light

The need for sensitively and reliably probing magnetization dynamics has been increasing in various contexts such as studying novel hybrid magnonic systems, in which the spin dynamics strongly and coherently couple to other excitations, including microwave photons, light photons, or phonons. Recent advances in quantum magnonics also highlight the need for employing the magnon phase as quantum state variable, which is to be detected and mapped out with high precision in on-chip micro- and nanoscale magnonic devices. Here, in this study, we demonstrate a facile optical technique that can directly perform concurrent spectroscopic and imaging functionalities with spatial and phase resolutions, using infrared strobe light operating at 1550-nm wavelength. To showcase the methodology, we spectroscopically studied the phaseresolved spin dynamics in a bilayer of Permalloy and yttrium iron garnet Y 3 Fe 5 O 12 (YIG), and spatially imaged the backward-volume spin-wave modes of YIG in the dipolar spin-wave regime. Using the strobe light probe, the detected precessional phase contrast can be directly used to construct the map of the spin wave's wave front, in the continuous-wave regime of spin-wave propagation and in the stationary state, without needing any optical reference path. By selecting the applied field, frequency, and detection phase, the spin-wave images can be made sensitive to the precession amplitude and phase. Our results demonstrate that infrared optical strobe light can serve as a versatile platform for magneto-optical probing of magnetization dynamics, with potential implications in investigating hybrid magnonic systems.

Xiong, Yuzan [University of North Carolina, Chapel↗

Developing socially and economically generative, resilient PV-energy systems for low- and moderate-income communities: Applications for Puerto Rico (Final Technical Report)

The nexus of energy and poverty is one of the most significant long-term challenges facing the energy and electricity sectors, with negative feedback loops between energy-related costs, risks, and insecurities, on the one hand, and diverse other social and economic insecurities, on the other, helping to undermine health, wellbeing, and resilience in low- and moderate-income communities around the US and the world (Biswas et al. 2022). The transition of energy systems to clean energy alternatives offers a potential opportunity to redesign future energy arrangements in ways that, instead, create beneficial social value for communities, reverses the energy-poverty nexus, and enhances community wellbeing (Biswas et al. 2020; Miller et al. 2018). The promise of clean energy transitions to address the energy-poverty nexus is especially significant when designed using principles of just energy transitions. To leverage just clean energy transitions to accomplish a reversal of the energy-poverty nexus requires improvement in three capabilities: (1) Understanding and mapping the community-scale dynamics and variability of the energy-poverty nexus, which vary considerably from community to community; (2) Evaluating the technical, social, and economic potential for different clean energy system designs to deliver social value to specific communities; and (3) Working collaboratively with communities to imagine and design clean energy systems solutions that meet their needs, fit their capabilities and contexts, and deliver multiple forms of community-desired social value/benefits that help reverse key facets of the energy-poverty nexus. This project developed a portfolio of novel tools, methodological approaches, and collaborative partnerships with four communities in Puerto Rico that further the objective of building these three capabilities.

14 SOLAR ENERGY↗

McMillan map and nonlinear Twiss parameters

In this article we consider two dynamical systems: the McMillan sextupole and octupole integrable mappings originally introduced by Edwin McMillan; the second one is also known as canonical McMillan map. Both of them are simplest symmetric McMillan maps with only one intrinsic parameter, the trace of the Jacobian at the fixed point. While these dynamical systems have numerous of applications and are used in many areas of math and physics, some of their dynamical properties have not been described yet. We fulfill the gap and provide complete description of all stable trajectories including parametrization of invariant curves, Pioncaré rotation numbers and canonical action-angle variables. In the second part we relate these maps with general chaotic map in McMillan-Turaev form. We show that McMillan sextupole and octupole mappings are first order approximations of dynamics around the fixed point, in a similar way as linear map and quadratic invariant (Courant-Snyder invariant in accelerator physics) is the zeroth order approximation (known as linearization). Finally we suggest the new formalism of nonlinear Twiss parameters which incorporate dependence of rotation number as a function of amplitude, in contrast to e.g. betatron phase advance used in accelerator physics which is independent of amplitude. Specifically in application to accelerator physics this new formalism is capable of predicting dynamical aperture around 1-st, 2-nd, 3-rd and 4-th order resonances for flat beams, which is critical for beam injection/extraction.

43 PARTICLE ACCELERATORS↗

Adaptation Strategies Strongly Reduce the Future Impacts of Climate Change on Simulated Crop Yields

Abstract Simulations of crop yield due to climate change vary widely between models, locations, species, management strategies, and Representative Concentration Pathways (RCPs). To understand how climate and adaptation affects yield change, we developed a meta‐model based on 8703 site‐level process‐model simulations of yield with different future adaptation strategies and climate scenarios for maize, rice, wheat and soybean. We tested 10 statistical models, including some machine learning models, to predict the percentage change in projected future yield relative to the baseline period (2000–2010) as a function of explanatory variables related to adaptation strategy and climate change. We used the best model to produce global maps of yield change for the RCP4.5 scenario and identify the most influential variables affecting yield change using Shapley additive explanations. For most locations, adaptation was the most influential factor determining the projected yield change for maize, rice and wheat. Without adaptation under RCP4.5, all crops are expected to experience average global yield losses of 6%–21%. Adaptation alleviates this average projected loss by 1–13 percentage points. Maize was most responsive to adaptive practices with a projected mean yield loss of −21% [range across locations: −63%, +3.7%] without adaptation and −7.5% [range: −46%, +13%] with adaptation. For maize and rice, irrigation method and cultivar choice were the adaptation types predicted to most prevent large yield losses, respectively. When adaptation practices are applied, some areas are predicted to experience yield gains, especially at northern high latitudes. These results reveal the critical importance of implementing adequate adaptation strategies to mitigate the impact of climate change on crop yields.

54 ENVIRONMENTAL SCIENCES↗