Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Limb”

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 19 records

CRCNS21 Computational Models of Multisensory Integration by Upper Limb in Humanoids and Amputees

This international collaborative research project between Johns Hopkins University (JHU) and the Technical University of Munich (TUM) investigated how the human brain processes and integrates multiple types of sensory information, such as touch and force, with the goal of improving prosthetic limbs for amputees and advancing sensory capabilities in humanoid robots. The research advanced our understanding of how the brain responds to sensory feedback in upper-limb amputees. Through experiments in which amputees received electrical stimulation while performing phantom hand movements, we demonstrated that sensory feedback activates the cortical sensorimotor and multisensory regions, and that these regions communicate dynamically during stimulation. Experiments with intact-limb participants explored the integration of visual, haptic, and force feedback, as well as in virtual reality motor training, further showing how the brain processes multimodal sensory information. In addition, this research inspired work on examining the reliability of where amputees perceive sensations over time, which contributed to a successful doctoral fellowship for continued investigation. Our collaborators at TUM improved multimodal sensor technology combining tactile and thermal feedback for humanoid robots, demonstrating the feasibility of integrating multiple sensor types into a unified system for detecting and responding to environmental stimuli. The experimental methods and analysis techniques developed across both teams, including functional network analysis and multimodal sensor integration, provide a foundation for future research in prosthetics and robotics. This research benefits the public by generating knowledge about how amputees process restored sensory information. Advances in humanoid sensing contribute to safer human-robot interaction. The project also fostered international collaboration and cross-disciplinary training: one TUM doctoral student spent a summer at JHU working on multimodal sensor integration, while two JHU students traveled to TUM to host workshops on neuromorphic sensory encoding and sensory integration.

42 ENGINEERING

Ion Acceleration in Fermi-LAT Behind-the-Limb Solar Flares: The Role of Coronal Shock Waves

Here, we investigate the relationship between the gamma-ray emission measured with Fermi-LAT and radio signatures of coronal shock waves in four behind-the-limb (BTL) solar flares. All events were associated with metric type II radio burst. Both start and end times of the radio bursts were synchronized with the gamma-ray emission. The type II bursts associated with the BTL gamma-ray flares had higher speeds and lower formation heights than those of an average sample. These findings support the notion that the highly relativistic ions that produce the gamma-rays in BTL flares are accelerated at CME-driven propagating coronal shock waves rather than in large-scale coronal loops.

Corona

Range extender mediates long-distance enhancer activity

Although most mammalian transcriptional enhancers regulate their cognate promoters over distances of tens of kilobases, some enhancers act over distances in the megabase range1. The sequence features that enable such long-distance enhancer–promoter interactions remain unclear. Here we used in vivo enhancer-replacement experiments at the mouse Shh locus to show that short- and medium-range limb enhancers cannot initiate gene expression at long-distance range. We identify a cis-acting element, range extender (REX), that confers long-distance regulatory activity and is located next to a long-range limb enhancer of Sall1. The REX element has no endogenous enhancer activity. However, addition of the REX to other short- and mid-range limb enhancers substantially increases their genomic interaction range. In the most extreme example observed, addition of REX increased the range of an enhancer by an order of magnitude from its native 73 kb to 848 kb. The REX element contains highly conserved [C/T]AATTA homeodomain motifs that are critical for its activity. These motifs are enriched in long-range limb enhancers genome-wide, including the ZRS (zone of polarizing activity (ZPA) regulatory sequence), a benchmark long-range limb enhancer of Shh2. The ZRS enhancer with mutated [C/T]AATTA motifs maintains limb activity at short range, but loses its long-range activity, resulting in severe limb reduction in knock-in mice. In summary, we identify a sequence signature associated with long-range enhancer–promoter interactions and describe a prototypical REX element that is necessary and sufficient to confer long-distance activation by remote enhancers.

Bower, Grace

Arm and shoulder muscle segmentation in axial MRI with UNet deep learning model

Quantifying individual upper-limb muscle volumes from MRI provides key insight into muscle-specific strength, deficits, and adaptations. Manual delineation is the gold standard but time‑intensive, and the performance of current deep learning approaches, particularly for small or anatomically complex muscles, remains incompletely characterized. We evaluated a state‑of‑the‑art deep learning framework across the entire upper limb and analyzed factors governing segmentation performance, with attention to the forearm. Three previously published MRI datasets (1.5 T, 3D GRE T1‑weighted; total n = 39) spanning young, middle‑aged, and older adults were curated and quality‑checked, including expert manual segmentations for 31 muscles. Following multiclass mask reconstruction, we trained three 3D nnU‑Net multiclass models matched to the muscle subsets present across datasets, using five‑fold cross‑validation and a composite Dice Similarity Coefficient (DSC) + cross entropy loss. Segmentation accuracy was assessed with DSC. Performance varied across muscles (mean DSC = 0.806 ± 0.098), ranging from 0.920 (Deltoid) to 0.461 (Extensor pollicis brevis). In uncertainty‑weighted regressions, muscle volume was positively associated with DSC (R2 = 0.36, p < 0.001), whereas training segmentation count and muscle orientation showed negligible associations (R2 ≤ 0.06). A weighted mixed‑effects model identified volume as the strongest evaluated predictor, explaining 23.9% of variance in DSC; orientation and training count each contributed <1%, leaving 61.5% unexplained. These results indicate that deep learning–based segmentation can accurately quantify muscle volume for many upper‑limb muscles but remains constrained for small, low‑contrast forearm muscles.

Gillespie, Samuel

Intercomparison of flood inundation models across land use types and hydrological flood stages

Flood Inundation Mapping (FIM) model selection is a key operational decision because accurate, rapid mapping underpins early warning and resource allocation. FIM performance is context-dependent and can vary with hydrograph phase, land-use/land-cover (LULC), and the evaluation benchmark. Intercomparison studies typically assess a single near-peak snapshot against one reference dataset. Here, we provide a context-stratified intercomparison across (i) multiple hydrograph phases, (ii) LULC classes, and (iii) benchmark types, for five FIM approaches spanning a wide range of physical complexity and operational cost (TRITON, LISFLOOD-FP, HEC-RAS 2D, ARC-Curve2Flood, and OWP HAND-FIM). We use the Hurricane Matthew flood (2016) in the Neuse River Basin, North Carolina, USA, as a case study. Using high-resolution remote sensing-derived flood inundation maps, hand-labeled points, and building footprints, we assess model skill across two rising and two falling hydrograph limbs and across major LULC types. Results show that model rankings shift systematically across contexts: LISFLOOD-FP ranks highest in three of four flood phases, while TRITON leads during one rising limb phase; LISFLOOD-FP performs best in vegetated areas, whereas HEC-RAS improves relative performance in agricultural and urban areas; and benchmark choice influences conclusions, with LISFLOOD-FP performing best for flooded-building detection in the late falling limb, while TRITON ranks highest against hand-labeled points. We also report representative wall-clock runtimes for each workflow to provide use-case context for operational feasibility. Together, these results offer transferable guidance for model selection and for designing large-scale, benchmark-aware FIM intercomparison studies.

Nikrou, Parvaneh [University of Alabama]

A gene desert required for regulatory control of pleiotropic Shox2 expression and embryonic survival

Approximately a quarter of the human genome consists of gene deserts, large regions devoid of genes often located adjacent to developmental genes and thought to contribute to their regulation. However, defining the regulatory functions embedded within these deserts is challenging due to their large size. Here, we explore the cis-regulatory architecture of a gene desert flanking the Shox2 gene, which encodes a transcription factor indispensable for proximal limb, craniofacial, and cardiac pacemaker development. We identify the gene desert as a regulatory hub containing more than 15 distinct enhancers recapitulating anatomical subdomains of Shox2 expression. Ablation of the gene desert leads to embryonic lethality due to Shox2 depletion in the cardiac sinus venosus, caused in part by the loss of a specific distal enhancer. The gene desert is also required for stylopod morphogenesis, mediated via distributed proximal limb enhancers. In summary, our study establishes a multi-layered role of the Shox2 gene desert in orchestrating pleiotropic developmental expression through modular arrangement and coordinated dynamics of tissue-specific enhancers.

59 BASIC BIOLOGICAL SCIENCES

The Role of Snowmelt and Subsurface Heterogeneity in Headwater Hydrology of a Mountainous Catchment in Colorado: A Model‐Data Integration Approach

Mountainous headwater streams are sustained by both snowmelt‐driven streamflow and groundwater discharge in the Upper Colorado River Basin. However, predicting headwater stream discharge magnitude and peak flow timing is challenging in mountainous terrains, where snowmelt rates vary with vegetation type and elevation, and heterogeneous subsurface physical properties influence groundwater storage and its release. We used a model‐data integration approach to investigate the roles of snowmelt and subsurface structure in stream discharge and groundwater level. We ran an ensemble of 100 integrated surface‐subsurface hydrologic models for a mountainous headwater catchment near Crested Butte, Colorado, USA. We also evaluated and calibrated these models against observed data sets, including snow depth measurements using distributed temperature probes, stream discharge, and groundwater levels. Calibration with multiple data sources using neural density estimators has further constrained uncertainty in subsurface properties and snowmelt rates. Results indicated that observed slower snowmelt rates in evergreen forests delayed the peak flow and baseflow onset. In upstream areas with lower subsurface permeability, water was stored within the subsurface but was not released as interflow or shallow groundwater flow, and thereby not contributing to downstream streamflow during recession limb periods. Double peaks in groundwater occurred in areas with spatial subsurface heterogeneity, in our case due to the contrast between granodiorite and Mancos shale. These process‐based insights into groundwater and snowmelt dynamics in mountainous headwaters will help improve predictions of headwater hydrology.

Wang, Lijing [University of Connecticut, Storrs, C

In vivo mapping of mutagenesis sensitivity of human enhancers

Distant-acting enhancers are central to human development1. However, our limited understanding of their functional sequence features prevents the interpretation of enhancer mutations in disease2. Here we determined the functional sensitivity to mutagenesis of human developmental enhancers in vivo. Focusing on seven enhancers that are active in the developing brain, heart, limb and face, we created over 1,700 transgenic mice for over 260 mutagenized enhancer alleles. Systematic mutation of 12-base-pair blocks collectively altered each sequence feature in each enhancer at least once. We show that 69% of all blocks are required for normal in vivo activity, with mutations more commonly resulting in loss (60%) than in gain (9%) of function. Using predictive modelling, we annotated critical nucleotides at the base-pair resolution. The vast majority of motifs predicted by these machine learning models (88%) coincided with changes in in vivo function, and the models showed considerable sensitivity, identifying 59% of all functional blocks. Taken together, our results reveal that human enhancers contain a high density of sequence features that are required for their normal in vivo function and provide a rich resource for further exploration of human enhancer logic.

Kosicki, Michael

Altered post-fracture systemic bone loss in a mouse model of osteocyte dysfunction

Femur fracture leads to loss of bone at uninjured skeletal sites, which may increase risk of subsequent fracture. Osteocytes, the most abundant bone cells, can directly resorb bone matrix and regulate osteoclast and osteoblast activity, but their role in systemic bone loss after fracture remains poorly understood. In this study we used a transgenic (TG+) mouse model that overexpresses human B-cell lymphoma 2 (BCL-2) in osteoblasts and osteocytes. This causes enhanced osteoblast proliferation, followed by disruption in lacunar-canalicular connectivity and massive osteocyte death by 10 wk of age. We hypothesized that reduced viable osteocyte density would decrease the magnitude of systemic bone loss after femur fracture, reduce perilacunar remodeling, and alter callus formation. Bone remodeling was assessed using serum biomarkers of bone formation and resorption at 5 d post-fracture. We used micro-computed tomography, high resolution x-ray microscopy, mechanical testing, and Raman spectroscopy to quantify the magnitude of systemic bone loss, as well as changes in osteocyte lacunar volume, bone strength, and bone composition 2 wk post-fracture. Fracture was associated with a reduction in circulating markers of bone resorption in non-transgenic (TG-) animals. TG+ mice exhibited high bone mass in the limbs, greater cortical elastic modulus and reduced post-yield displacement. After fracture, TG+ mice lost less trabecular bone than TG- mice, but conversely TG+ mice exhibited trends toward a lower yield point and reduced femoral cortical thickness after fracture, though these were not statistically significant. Lacunar density was greater in TG+ mice, but fracture did not alter lacunar volume in TG+ or TG- mice. These findings suggest that osteocytes potentially play a significant role in the post-traumatic systemic response to fracture, though the effects differ between trabecular and cortical bone.

60 APPLIED LIFE SCIENCES

slewpy

slewpy is a Python package that allows the simulation of the science operations of an astrophysics space satellite mission. slewpy allows users to specify an astronomical target list with observing priorities and a satellite configuration (i.e., orbit and various satellite parameters). Taking these inputs, slewpy can be used to run a time-resolved simulation of an astrophysics mission and outputs simulated target observations given constraints such as satellite slewing rates between targets, observing time on a given target, and Sun-, Earth-, and Moon-limb observing constraints. slewpy provides the tools to test astrophysics space satellite mission designs against observing requirements for a given science case.

Geringer-Sameth, Alex [Lawrence Livermore National

Cellular signaling within aged skeletal muscle reveals a dysregulated stress-induced remodeling response following volumetric muscle loss in female mice

Severe muscle trauma disrupts endogenous repair mechanisms, producing chronic functional deficits that are incompletely characterized in aged populations. This study investigated inflammatory, molecular, and physiological responses to volumetric muscle loss in young adult and aged female mice. Serum cytokine profiling revealed elevated baseline inflammation in aged animals and a blunted response to injury, with cytokines such as IL-6 increasing 4.4-fold in young versus 1.8-fold in aged mice at day 3 compared to baseline. By day 28 post-injury, histological analyses showed comparable reductions in muscle size and increases in fibrosis across ages. Despite these similar tissue-level outcomes, age-dependent differences emerged in downstream functional and molecular responses. Muscle functional testing demonstrated persistent force deficits independent of age but altered muscle relaxation kinetics in aged mouse muscles (p < 0.001), suggesting dysregulated excitation-contraction coupling. Additionally, mice displayed age-associated differences in post-injury limb loading. Global proteomic analyses further confirmed age-associated enrichment of complement and antigen-processing pathways alongside metabolic dysfunction (p < 0.05). Phosphoproteomic profiling revealed reduced basal kinase activity in the muscles of aged mice yet exaggerated injury-induced phosphorylation of Mapk1-associated phosphosites, indicating a dysregulated stress response. Collectively, these findings suggest that aged murine muscles function within a heightened inflammatory and perturbed kinase-signaling environment that may hinder the coordination of regenerative programs, underscoring the need for regenerative strategies that address age-specific molecular contexts to improve functional recovery across the lifespan.

Habing, Krista M.

Groundwater elevation data for monitoring wells within the East and Taylor River basins, Colorado (USA)

This dataset is comprised of temporal variations in groundwater elevation data for the 24 monitoring wells located throughout the East River watershed. Seasonal to annual variations in groundwater elevations are a critical property of mountainous watersheds needed to understand both hydrological and below ground biogeochemical processes. Such data serve as a critical constraint for numerical models describing coupled groundwater-surface water behavior within the watershed. Additionally, the offset between the maximum and minimum groundwater elevations defines the extent of the bedrock weathering zone, with annual excursions in the groundwater hydrographic (i.e., the rising and falling hydrographic limbs) imposing primary controls on bedrock saturation state and redox conditions that govern biogeochemical reactions impacting nitrogen, carbon, and metals cycling. Manufacturer-specific software is used to download pressure data from each transducer, with broadly available spreadsheet software (e.g. Microsoft Excel) used to convert temporal variations in water pressure to elevations in units of meters above mean sea level. As additional monitoring wells are installed within the East River watershed and new groundwater monitoring wells are installed in the Taylor River watershed, temporal groundwater elevation data will be included as a part of this master dataset. Details regarding the metadata associated with each monitoring well location, including well depths, screened intervals, well location coordinates, and bedrock type, are included, as is a standard operating procedure for generating groundwater elevation data from water pressure values recorded by the pressure transducers. This dataset includes: (1) a zip file (East_River_Watershed_Compiled_Groundwater_Elevation_Data_Plots.zip), containing (a) PNG of groundwater hydrographs, (b) a CSV file with groundwater elevation data, and (c) CSV file containing metadata organized by location; (2) an Excel file (East_River_Watershed_Compiled_Groundwater_Elevation_Data_Plots.xlsx) with the groundwater elevation data, groundwater hydrographs, and metadata organized by location; (3) a Word file (Groundwater_elevation_data_protocols.docx) and a PDF file version (Groundwater_elevation_data_protocols.pdf) containing field protocols and methods; (4) a location metadata (locations.csv) file; (5) a file level metadata (flmd.csv); and (6) data dictionary (dd.csv) file.This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

54 ENVIRONMENTAL SCIENCES

Z-Target Radiography Postprocessing With A Deep Convolution Neural Network

Analyzing X-ray radiographs is crucial for understanding target behavior in Inertial Confinement Fusion (ICF) and High Energy Density (HED) platforms. However, the density of Magneto Raleigh Taylor (MRT) bands and limitations of target materials often obscure relevant spike growth and density information. To address this issue, machine learning postprocessing techniques can be applied to remove darkened regions in radiography images. In this study, a novel method is presented for removing MRT darkened regions from z-target radiographs using a convolutional neural network (CNN). The CNN, consisting of six layers, treats the darkened regions as noise and employs a mixed loss function and end-to-end frameworks to suppress them while preserving sharpness. The six-layer architecture is designed to effectively learn features when provided with a larger volume of learning space. Each layer is optimized using a mixed loss function that combines a standard loss pixel approach with a multi-scaled structural similarity index loss, which considers luminance, contrast, and structure in local neighborhoods. This approach is particularly beneficial for capturing the stochastic structure of MRT limbs. Due to the limited availability of experimental data, training is conducted using synthetic target radiography from 3D Alegra simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Diagenesis is key to unlocking outcrop fracture data suitable for quantitative extrapolation to geothermal targets

Exceptionally large, well-exposed sandstone outcrops in New York provide insights into folds, deformation bands, and fractures that could influence permeability, heat exchange, and stimulation outcomes of geothermal reservoir targets. Cambrian Potsdam Sandstone with <5% porosity contains decimeter-scale open, angular-limbed monoclines <0.5 km apart with associated low-porosity mm-wide cataclastic deformation bands. Crossing and abutting relationships among sub-vertical opening-mode fractures show four chronological Sets A–D, striking NNW, NE, NW, and ENE, respectively. Fracture lengths and heights range from millimeters to tens of meters. Sets A and C macro-fractures, and possibly B and D, contain quartz deposits. All sets have abundant associated quartz cemented microfractures that also record set orientations and crosscutting relations. Quartz cement deposits—evidence of diagenesis—are the key to identifying attributes of outcrop fractures suitable for extrapolation to geothermal targets in sandstones because they show which fractures formed in the subsurface. Set A fluid inclusion homogenization temperatures (120°C–129°C) are compatible with fracture at >3 km depth. Fractures are stiff and those ≥0.05 mm (Set C) and ≥0.1 mm (Set A) are open and potentially conducive to flow. Sets A and D are abundant in outcrops with close fracture spacing—0.18 m and 0.68 m, respectively—and define a rectangular connectivity network dominated by crossing and abutting X and Y nodes. Set A aperture distributions follow a power law with slope –0.8 up to 0.15 mm; other sets have lognormal distributions. Set A and D microfractures are weakly clustered, while macro-fractures commonly have 1D anticlustered (regular or periodic) arrangements at shorter length scales (<0.2 m). Sub-horizontal fractures are barren and may have formed near the surface. Fracture heights, lengths, and spatial arrangements show good trace connectivity but low open connectivity. For geothermal applications, outcrop results predict low initial well-test permeabilities owing to quartz disconnecting open fractures, but stimulation of closely spaced microfractures and partly open macro-fractures could yield high surface area for heat exchange. Quantitative extrapolation of key fracture attributes like abundance, orientation, spatial arrangement, length, and open fracture connectivity is possible from outcrops to fractured reservoirs if differing thermal histories and diagenesis are accounted for.

02 PETROLEUM

Hunga Tonga–Hunga Ha′apai Volcano Impact Model Observation Comparison (HTHH-MOC) project: experiment protocol and model descriptions

The 2022 Hunga volcanic eruption injected a significant amount of water vapor and a moderate amount of sulfur dioxide into the stratosphere, causing observable responses in the climate system. We have developed a model–observation comparison project to investigate the evolution of volcanic water and aerosols and their impacts on atmospheric dynamics, chemistry, and climate, using several state-of-the-art chemistry climate models. The project goals are (1) to evaluate the current chemistry–climate models to quantify their performance in comparison to observations and (2) to understand atmospheric responses in the Earth system after this exceptional event and investigate the potential impacts in the projected future. To achieve these goals, we designed specific experiments for direct comparisons to observations, for example from balloons and the Microwave Limb Sounder satellite instrument. Experiment 1 consists of two sets of free-running ensemble experiments from 2022 to 2031: one with fixed sea-surface temperatures and sea ice and one with coupled ocean. These experiments will help to understand the long-term evolution of water vapor and aerosols; quantify HTHH effects on stratospheric and mesospheric temperatures, dynamics, and transport; understand the impact of dynamic changes on ozone chemistry; quantify the net radiative forcings; and evaluate any surface climate impact. Experiment 2 is a nudged-run experiment from 2022 to 2023 using observed meteorology. To allow participation of more climate models with varying complexities of aerosol simulation, we include two sets of simulations in Experiment 2: Experiment 2a is designed for models with internally generated aerosol, while Experiment 2b is designed for models using prescribed aerosol surface area density. This experiment will help to analyze H 2 O and aerosol evolution, quantify the net radiative forcings, understand the impacts on mid-latitude and polar O 3 chemistry, and allow close comparisons with observations.

Zhu, Yunqian [Univ. of Colorado, Boulder, CO (Unit

Empowering Lineworkers: The Case for Active Exoskeletons in Utility Work

Exoskeletons have evolved from early medical prototypes to advanced systems capable of addressing physical demands in various industries. This report explores the potential of active exoskeleton technology within the utility sector, focusing on its application for linemen who face significant risks of work-related musculoskeletal disorders (WMSDs). By analyzing existing literature on exoskeletons across industries such as construction, manufacturing, and military, the study identifies a gap in utility-specific applications. Task-specific design features like gravity compensation, limb support, and advanced safety measures, improve exoskeletons’ potential to alleviate physical strain, reduce workplace injuries, and enhance productivity. This review emphasizes the need for targeted research and development to optimize exoskeleton designs for the utility sector to provide benefits for workers, companies, and the broader community.

60 APPLIED LIFE SCIENCES

Cultivating clarity: understanding the impact of land cost assumptions on biofuel viability

The transition to sustainable energy has increased interest in biofuel production to reduce greenhouse gas emissions, decrease reliance on imported oil, and ensure energy resilience. Here, this study examines the often-overlooked impact of land cost assumptions on the economic viability of biofuel production. Using a discounted cash flow techno-economic framework, we evaluated three land cost scenarios—no land costs, land rental costs, and land purchase costs—across six bioenergy feedstocks (corn, soybeans, switchgrass, miscanthus, poplar, and microalgae) and three biofuel products (corn ethanol, soybean biodiesel, and sustainable aviation fuel) at the county level for the contiguous United States. The analysis reveals substantial variation in minimum fuel selling prices due to these scenarios. High-yield crops like algae showed low sensitivity to land costs, while low-yield crops such as soybeans were highly sensitive. Geographical differences were significant, with minimum fuel selling price increases most pronounced in high-value land regions like the Corn Belt. Case studies further illustrate the influence of local productivity and land costs on economic outcomes across the United States. These findings emphasize the importance of maintaining consistent land cost assumptions in biofuel economic assessments. By quantifying the interplay between land value, crop productivity, and economic feasibility, this study provides essential insights for policymakers and stakeholders to advance sustainable energy solutions.

09 BIOMASS FUELS

Multi-objective optimization of sustainable aviation fuel production pathways in the U.S. Corn Belt

As a potential source of low-carbon transportation energy, biofuels offer certain advantages over vehicle electrification (e.g., lower societal vulnerability to grid failures, and improved range of sustainable aviation), but also several challenges, including cost, carbon intensity, and land usage. There are also well-founded concerns that biofuel supply chains could be disrupted if extreme weather events impact feedstock yields. In this paper, we explore the use of multi-objective optimization to identify biofuel production pathways that balance cost, greenhouse gas emissions, and supply vulnerability to extreme weather. We compare the use of three different many-objective evolutionary algorithms and linear programming in optimizing biomass cultivation decisions in the U.S. Corn Belt under weather uncertainty using historical, modeled, and synthetic yield data. We consider four feedstock choices (corn, soy, switchgrass, and algae) with two land types (agricultural and marginal lands) and evaluate decisions using three alternative spatial resolutions (ranging from the USDA agricultural district level to the state level). Results show that feedstock choice is the primary driver of objective performance (i.e., the position and shape of 3D, approximate Pareto frontiers). Spatial diversification is a less effective tool in reducing exposure to weather-caused drops in crop yield.

09 BIOMASS FUELS