SEARCH · Engineering Papers

Results for “variable selection”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 631 records · Page 35

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Interfacial Chemistry Involved in Selective Separation of NMC/LMO and LCO/LMO Binary Cathode Materials by Froth Flotation Using Oleic Acid

The variability in cathode compositions within recycled lithium-ion battery (LIB) feedstocks poses a significant challenge to efficient downstream refining processes. This study demonstrates the feasibility of using froth flotation with oleic acid as a collector to selectively separate lithium nickel-manganese-cobalt oxide (NMC) and lithium cobalt oxide (LCO) from lithium manganese oxide (LMO) materials. Laboratory-scale flotation tests achieved an 80% separation efficiency in a single stage, producing a froth product with >90% purity of NMC/LCO at approximately 90% yield. Concurrently, the LMO materials were enriched in the sink product with ∼90% purity and ∼90% yield. This approach was further validated using recycled cathode materials, confirming its applicability to realistic feedstocks. The underlying mechanism governing the selective separation of NMC/LCO from LMO was investigated using ζ-potential measurements, contact angle measurements, bubble-particle attachment experiments, and X-ray photoelectron spectroscopy (XPS) analysis. Both contact angle and bubble-particle attachment results confirmed that oleic acid adsorption rendered NMC and LCO surfaces hydrophobic, thereby enhancing flotation recovery. At pH 5, oleic acid adsorbed preferentially onto NMC and LCO surfaces via electrostatic interactions, while exhibiting minimal adsorption on LMO surfaces. However, separation efficiency deteriorated at higher pH, which was attributed to the co-flotation of LMO materials caused by oleate chemisorption on MnOH + species. This work establishes froth flotation as a viable cathode/cathode separation strategy, providing a low-cost, scalable pathway to preconcentrate and enrich nickel-rich and cobalt-rich cathode active materials from incompatible cathode chemistries for direct recycling or hydrometallurgical processing. Furthermore, this study reveals, for the first time, the mechanism of oleate adsorption on the surface of different cathode materials.

PH↗

NASA NTRS1976

Toward the modeling of land use change: A spatial analysis using remote sensing and historical data

It was hypothesized that the chronological observation of land use change could be shown to follow a predictable pattern and these patterns could be correlated with other statistical data to develop transition probabilities suitable for modeling purposes. A literature review and preliminary research, however, indicated a totally stochastic approach was not practical for simulating land use change and thus a more deterministic approach was adopted. The approach used assumes the determinants of the land use conversion process are found in the market place, where land transactions among buyers and sellers occur. Only one side of the market transaction process is studied, however, namely, the purchaser's desires in securing an ideal or suitable site. The problem was to identify the ideal qualities, quantities or attributes desired in an industrial site (or housing development), and to formulate a general algorithmic statement capable of identifying potential development sites. Research procedures involved developing a list of variables previously noted in the literature to be related to site selection and streamlining the list to a set suitable for statistical testing. A sample of 157 industries which have located (or relocated) in the 16-county Knoxville metropolitan region since 1950 was selected for industrial location analysis. Using NASA color infrared photography and Tennessee Valley Authority historical aerial photography, data were collected on the spatial characteristics of each industrial location event. These data were then subjected to factor analysis to determine the interrelations of variables.

Honea, R. B.↗

NASA NTRS1989

Seasat over-land scatterometer data. II - Selection of extended area land-target sites for the calibration of spaceborne scatterometers

The post-launch performance verification for future scatterometers can use extended area land targets to calibrate antenna gain patterns and to verify and monitor deployment configurations. For the Ku-band Seasat scatterometer, a region of tropical rain forest in the Amazon basin was used as a homogeneous extended-area land target. As this region is continuously being deforested, other regions are investigated for calibrating scatterometers. The global backscatter coefficients (sigma0) are compared to classifications of natural vegetation and cultivation intensity and the variability with time during the three-month mission is studied. The statistical variability of sigma0 is compared with prior estimates resulting from the known variability of the instrument parameters and communication noise. Data from selected forested regions with relatively homogeneous sigma0 and little time-dependence are presented.

Kennett, Rosemary G.↗

DOE OSTI2022

Habitat use by female desert tortoises suggests tradeoffs between resource use and risk avoidance

Animals may select habitat to maximize the benefits of foraging on growth and reproduction, while balancing competing factors like the risk of predation or mortality from other sources. Variation in the distribution of food resources may lead animals to forage at times or in places that carry greater predation risk, with individuals in poor quality habitats expected to take greater risks while foraging. We studied Mojave desert tortoises ( Gopherus agassizii ) in habitats with variable forage availability to determine if risk aversion in their selection of habitat relative was related to abundance of forage. As a measure of risk, we examined tortoise surface activity and mortality. We also compared tortoise body size and body condition between habitats with ample forage plants and those with less forage plants. Tortoises from low forage habitats selected areas where more annual plants were nutritious herbaceous flowering plants but did not favor areas of greater perennial shrub cover that could shelter them or their burrows. In contrast, tortoises occupying high forage habitats showed no preference for forage characteristics, but used burrows associated with more abundant and larger perennial shrubs. Tortoises in high forage habitats were larger and active above ground more often but did not have better body condition. Mortality was four times higher for females occupying low forage habitat than those in high forage habitat. Our results are consistent with the idea that tortoises may minimize mortality risk where food resources are high, but may accept some tradeoff of greater mortality risk in order to forage optimally when food resources are limiting.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2018

Wildfire Mapping in Interior Alaska Using Deep Neural Networks on Imbalanced Datasets

Wildfires are the dominant disturbance impacting many regions in Alaska and is expected to intensify due to climate change. Accurate tracking and quantification of wildfires are important for climate modeling and ecological studies in this region. Remote sensing platforms (e.g., MODIS, Landsat) are an important tool for mapping wildfires in Alaska at varying spatial and temporal scales. Deep neural networks (DNN) have shown superior performance in many classification problems such as high-dimensional remote sensing data. Detection of wildfires is an imbalanced classification problem where one class contains a much smaller or larger sample size and DNNs performance can decline. We take a known weight-selection strategy during DNN training and apply them to MODIS variables (e.g., NDVI, surface reflectance) for a binary classification (i.e., wildfire or no-wildfire) across Alaska during the 2004 wildfire year, which is one of the largest on record. The method splits the input training data into sets, one for training the DNN to update weights and the other that splits into a validation set that is equally balanced between the wildfire and no-wildfire class. The performance is monitored on the validation set to select the weights based on the best validation-loss score. This approach was applied to two sampled datasets, such as where the no-wildfire class can significantly outweigh the wildfire class. The normal DNN training strategy was unable to map wildfires for the highly imbalanced dataset, however, the weight-selection strategy was able to map wildfires very accurately (0.96 recall score for 78,702 wildfire pixels (500×500 m)).

Langford, Zach↗

NASA NTRS1987

Decoupled direct method for sensitivity analysis in combustion kinetics

An efficient, decoupled direct method for calculating the first order sensitivity coefficients of homogeneous, batch combustion kinetic rate equations is presented. In this method the ordinary differential equations for the sensitivity coefficients are solved separately from , but sequentially with, those describing the combustion chemistry. The ordinary differential equations for the thermochemical variables are solved using an efficient, implicit method (LSODE) that automatically selects the steplength and order for each solution step. The solution procedure for the sensitivity coefficients maintains accuracy and stability by using exactly the same steplengths and numerical approximations. The method computes sensitivity coefficients with respect to any combination of the initial values of the thermochemical variables and the three rate constant parameters for the chemical reactions. The method is illustrated by application to several simple problems and, where possible, comparisons are made with exact solutions and those obtained by other techniques.

Radhakrishnan, Krishnan↗

NASA NTRS2012

Computational Work to Support FAP/SRW Variable-Speed Power-Turbine Development

The purpose of this report is to document the work done to enable a NASA CFD code to model the transition on a blade. The purpose of the present work is to down-select a transition model that would allow the flow simulation of a Variable-Speed Power-Turbine (VSPT) to be accurately performed. The modeling is to be ultimately performed to also account for the blade row interactions and effect on transition and therefore accurate accounting for losses. The present work is limited to steady flows. The low Reynolds number k-omega model of Wilcox and a modified version of same will be used for modeling of transition on experimentally measured blade pressure and heat transfer. It will be shown that the k-omega model and its modified variant fail to simulate the transition with any degree of accuracy. A case is therefore made for more accurate transition models. Three-equation models based on the work of Mayle on Laminar Kinetic Energy were explored and the Walters and Leylek model which was thought to be in a more mature state of development is introduced and implemented in the Glenn-HT code. Two-dimensional flat plate results and three-dimensional results for flow over turbine blades and the resulting heat transfer and its transitional behavior are reported. It is shown that the transition simulation is much improved over the baseline k-omega model.

Ameri, Ali A.↗

DOE OSTI2025

A Pentavalent HIV-1 Subtype C Vaccine Containing Computationally Selected gp120 Strains Improves the Breadth of V1V2 Region Responses

Background: HIV-1 envelope (Env) variable loops 1 and 2 (V1V2) directed non-neutralizing antibodies were a correlate of decreased transmission risk in the RV144 vaccine trial. Thus, the elicitation and breadth of antibody responses against the V1V2 of HIV-1 Env are important considerations for HIV-1 vaccine candidates. The V1V2 region’s highly variable nature and the extensive diversity of subtype C HIV-1 Envelopes (Envs) make the V1V2 response breadth a high priority for HIV-1 vaccine regimens aiming for V1V2-mediated protection in Southern Africa. Here, we determined whether the breadth of the anti-V1V2 vaccine response can be broadened by including HIV-1 Env strains computationally designed to enhance the coverage of subtype C V1V2 sequence diversity. Methods: Three subtype C Env strains were selected to maximize antibody binding coverage while complementing subtype C vaccine gp120s that were given in human clinical trials in South Africa, as well as to improve epitope accessibility. Humoral immunogenicity of a novel trivalent gp120 vaccine immunogen, a bivalent gp120 boost already in clinical trials (1086C and TV1), and a pentavalent (all five gp120s combined) were evaluated in a preclinical immunization study in guinea pigs. The pentavalent combination was further evaluated with alum versus glucopyranosyl lipid adjuvants formulated in squalene-in-water emulsion (GLA-SE) adjuvants in non-human primates. The breadth of the anti-V1V2 response was assessed using an array of cross-subtype variable loops 1&2 (V1V2) scaffold proteins and linear V2 peptides. Results: The breadth of the IgG response against V1V2 antigens of the trivalent and pentavalent groups was comparable, and both were greater than the breadth of the bivalent group. Linear epitope mapping showed that two linear epitopes in V2 were targeted by the vaccinated animals: the V2 hotspot focused at 169K that potentially correlated with decreased HIV-1 risk in RV144 and the V2.2 site (179LDV/I181) that is part of the integrin α4β7 binding site. The bivalent vaccine elicited a significantly higher magnitude of binding to the V2 hotspot compared to the trivalent vaccine whereas the trivalent vaccine elicited significantly higher binding to the V2.2 epitope compared to the bivalent vaccine, while the pentavalent recognized both regions. Conclusions: These results demonstrate that the three new computationally selected subtype C Envs successfully complemented 1086C and TV1 for broader V1V2 antibody responses, and, in concert with adjuvants that stimulate V1V2 responses, can be considered as part of a rationale immunogen design to improve V1V2 IgG coverage in future vaccine trials in South Africa.

Immunology↗

NASA NTRS1981

Calculation of sensitivity derivatives in thermal problems by finite differences

The optimum design of a structure subject to temperature constraints is considered. When mathematical optimization techniques are used, derivatives of the temperature constraints with respect to the design variables are usually required. In the case of large aerospace structures, such as the Space Shuttle, the computation of these derivatives can become prohibitively expensive. Analytical methods and a finite difference approach have been considered in studies conducted to improve the efficiency of the calculation of the derivatives. The present investigation explores two possibilities for enhancing the effectiveness of the finite difference approach. One procedure involves the simultaneous solution of temperatures and derivatives. The second procedure makes use of the optimum selection of the magnitude of the perturbations of the design variables to achieve maximum accuracy.

Haftka, R. T.↗

NASA NTRS1999

Psychological Selection of NASA Astronauts for International Space Station Missions

During the upcoming manned International Space Station (ISS) missions, astronauts will encounter the unique conditions of living and working with a multicultural crew in a confined and isolated space environment. The environmental, social, and mission-related challenges of these missions will require crewmembers to emphasize effective teamwork, leadership, group living and self-management to maintain the morale and productivity of the crew. The need for crew members to possess and display skills and behaviors needed for successful adaptability to ISS missions led us to upgrade the tools and procedures we use for astronaut selection. The upgraded tools include personality and biographical data measures. Content and construct-related validation techniques were used to link upgraded selection tools to critical skills needed for ISS missions. The results of these validation efforts showed that various personality and biographical data variables are related to expert and interview ratings of critical ISS skills. Upgraded and planned selection tools better address the critical skills, demands, and working conditions of ISS missions and facilitate the selection of astronauts who will more easily cope and adapt to ISS flights.

Galarza, Laura↗

NASA NTRS2003

Age affects severity of venous gas emboli on decompression from 14.7 to 4.3 psia

INTRODUCTION: Variables that define who we are, such as age, weight and fitness level influence the risk of decompression sickness (DCS) and venous gas emboli (VGE) from diving and aviation decompressions. We focus on age since astronauts that perform space walks are approximately 10 yr older than our test subjects. Our null hypothesis is that age is not statistically associated with the VGE outcomes from decompression to 4.3 psia. METHODS: Our data are from 7 different NASA tests where 188 men and 50 women performed light exercise at 4.3 psia for planned exposures no less than 4 h. Prebreathe (PB) time on 100% oxygen ranged from 150-270 min, including ascent time, with exercise of different intensity and length being performed during the PB in four of the seven tests with 150 min of PB. Subjects were monitored for VGE in the pulmonary artery using a Doppler ultrasound bubble detector for a 4-min period every 12 min. There were six design variables; the presence or absence of lower body adynamia and five PB variables; plus five concomitant variables on physical characteristics: age, weight height, body mass index, and gender that were available for logistic regression (LR). We used LR models for the probability of DCS and VGE, and multinomial logit (ML) models for the probability of Spencer VGE Grades 0-IV at exposure times of 61, 95, 131, 183 min, and for the entire exposure. RESULTS: Age was significantly associated with VGE in both the LR and ML models, so we reject the null hypothesis. Lower body adynamia was significant for all responses. CONCLUSIONS: Our selection of tests produced a wide range of the explanatory variables, but only age, lower body adynamia, height, and total PB time was helpful in various combinations to model the probability of DCS and VGE.

Clinical Trial↗

DOE OSTI2020

Dark Energy Survey Identification of A Low-Mass Active Galactic Nucleus at Redshift 0.823 from Optical Variability

We report the identification of a low-mass active galactic nucleus (AGN), DES J0218−0430, in a redshift z = 0.823 galaxy in the Dark Energy Survey (DES) Supernova field. We select DES J0218−0430 as an AGN candidate by characterizing its long-term optical variability alone based on DES optical broad-band light curves spanning over 6 yr. An archival optical spectrum from the fourth phase of the Sloan Digital Sky Survey shows both broad Mg ii and broad H β lines, confirming its nature as a broad-line AGN. Archival XMM–Newton X-ray observations suggest an intrinsic hard X-ray luminosity of |$L_{{\rm 2-12\, keV}}\approx 7.6\pm 0.4\times 10^{43}$| erg s^−1, which exceeds those of the most X-ray luminous starburst galaxies, in support of an AGN driving the optical variability. Based on the broad H β from SDSS spectrum, we estimate a virial black hole (BH) mass of M_• ≈ 10^6.43–10^6.72 M_⊙ (with the error denoting the systematic uncertainty from different calibrations), consistent with the estimation from OzDES, making it the lowest mass AGN with redshift > 0.4 detected in optical. We estimate the host galaxy stellar mass to be M_* ≈ 10^10.5 ± 0.3 M_⊙ based on modelling the multiwavelength spectral energy distribution. DES J0218−0430 extends the M_•–M_* relation observed in luminous AGNs at z ∼ 1 to masses lower than being probed by previous work. Our work demonstrates the feasibility of using optical variability to identify low-mass AGNs at higher redshift in deeper synoptic surveys with direct implications for the upcoming Legacy Survey of Space and Time at Vera C. Rubin Observatory.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

How will United States commercial building energy use be impacted by IPCC climate scenarios?

Climate change and anthropogenically-forced shift of weather in the future will impact energy use and resilience of both the built environment and the electric grid. The aim of this analysis is to understand how future climate scenarios will impact electricity and natural gas use of commercial buildings in the United States. Here, this study analyzes this impact for 2030, 2045, and 2100 using Representative Concentration Pathways (RCP) scenarios defined in Intergovernmental Panel on Climate Change (IPCC) Assessment Report 5. The large, gridded simulation of meteorological variables for RCPs 2.6, 4.5, 6.0, and 8.5 are selected and downscaled to make available hourly Future Meteorological Year (FMY) weather files for use and improvement in subsequent studies. High performance computing resources use these FMYs to simulate commercial prototype buildings in every American Society of Heating, Refrigeration, and Air Conditioning Engineers (ASHRAE) climate zone of the United States (US), and results are scaled to nation-wide energy use using conditioned floor area multipliers. The analysis is conducted without speculating the physical and performance traits of future buildings or the grid characteristics. This analysis quantifies the impact of climate change on source electrical and natural gas usage for commercial buildings in the United States over the next 80 years. If US commercial floorspace remained constant, total energy use by 2100 is predicted between an 1.75% decrease under the greatest emission scenario (8.5) and a 1.76% increase under the lowest emission scenario (2.6). When adjusted for anticipated urban growth by 2100, the predicted range is 65% increase (8.5) and 71% increase (2.6). Under a global temperature rise climate scenario, the warmest US climate zones will see a large increases in electricity use derived from space cooling while the coldest US climate zones will see significant decreases in natural gas use caused by the decrease in heating necessary. While climate change may ultimately require adaptations of the built environment to withstand its effects and because the United States is a country that requires more heating than cooling, from a building energy perspective, climate change (average temperature rise) is a net energy saver for the United States.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2023

Physics-based adaptivity of a spectral method for the Vlasov–Poisson equations based on the asymmetrically-weighted Hermite expansion in velocity space

We propose a spectral method for the 1D-1V Vlasov–Poisson system where the discretization in velocity space is based on asymmetrically-weighted Hermite functions, dynamically adapted via a scaling α and shifting u of the velocity variable. Specifically, at each time instant an adaptivity criterion selects new values of α and u based on the numerical solution of the discrete Vlasov–Poisson system obtained at that time step. Once the new values of the Hermite parameters α and u are fixed, the Hermite expansion is updated and the discrete system is further evolved for the next time step. The procedure is applied iteratively over the desired temporal interval. The key aspects of the adaptive algorithm are: the map between approximation spaces associated with different values of the Hermite parameters that preserves total mass, momentum and energy; and the adaptivity criterion to update α and u based on physics considerations relating the Hermite parameters to the average velocity and temperature of each plasma species. For the discretization of the spatial coordinate, we rely on Fourier functions and use the implicit midpoint rule for time stepping. The resulting numerical method possesses intrinsically the property of fluid-kinetic coupling, where the low-order terms of the expansion are akin to the fluid moments of a macroscopic description of the plasma, while kinetic physics is retained by adding more spectral terms. Moreover, the scheme features conservation of total mass, momentum and energy associated in the discrete, for periodic boundary conditions. A set of numerical experiments confirms that the adaptive method outperforms the non-adaptive one in terms of accuracy and stability of the numerical solution.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Iodine solid sorbent design: a literature review of the critical criteria for consideration

Designing sorbents for iodine capture in different conditions requires selection and optimization of a large and diverse range of variables. These variables fall into general categories (or features) of sorbent activity, sorbent stability, and the fate of the loaded material in terms of the disposal (waste form) options available. To illustrate, silver-loaded, high-porosity sorbents make for maximized iodine capture and less pressure drop in a column-based sorption system approach, however, this high porosity can lead to less mechanically stable sorbents. Additionally, waste forms containing silver must also be compliant with additional criteria for hazardous waste disposal. Thus, all these aspects must be considered simultaneously when selecting a sorbent for utilization under specific conditions. Information is given for different types of sorbent design considerations for different operating conditions and some emphasis is also given on promising alternatives for silver as the active (chemisorption-based) getter metal. Discussion is given around demonstrated options for waste forms for different metal-iodide compounds.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Muon identification using multivariate techniques in the CMS experiment in proton-proton collisions at sqrt(s) = 13 TeV

The identification of prompt and isolated muons, as well asmuons from heavy-flavour hadron decays, is an important task. Wedeveloped two multivariate techniques to provide highly efficientidentification for muons with transverse momentum greater than10 GeV. One provides a continuous variable as an alternative to acut-based identification selection and offers a betterdiscrimination power against misidentified muons. The other oneselects prompt and isolated muons by using isolation requirements toreduce the contamination from nonprompt muons arising inheavy-flavour hadron decays. Both algorithms are developed using59.7 fb$^{-1}$ of proton-proton collisions data at a centre-of-massenergy of √(s)=13 TeV collected in 2018 with the CMSexperiment at the CERN LHC.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2021

Quantifying nitrogen loss hotspots and mitigation potential for individual fields in the US Corn Belt with a metamodeling approach

The high productivity in the US Corn Belt is largely enabled by the consumption of millions of tons of manufactured fertilizer. Excessive application of nitrogen (N) fertilizer has been pervasive in this region, and the unrecovered N eventually escaped from croplands in forms of nitrous oxide (N 2 O) emission and N leaching. Mitigating these negative impacts is hindered by a lack of practical information on where to focus and how much mitigation potential to expect. At a large scale, process-based crop models are the primary tools for predicting variables required by decision making, but their applications are prohibited by expensive computational and data storage costs. To overcome these challenges, we built a series of metamodels to learn the key mechanisms regarding the carbon (C) and N cycle from a well-validated process-based biogeochemical model, ecosys. The trained metamodel captures over 98% of the variability of the ecosys simulated outputs for 99 randomly selected counties in Iowa, Illinois, and Indiana. To identify hotspots with high mitigation potential, we introduce net societal benefit (NSB) as an indicator for synthesizing the loss in yield and social benefits through emissions and pollutants avoided. Our results show that reducing N fertilizer by 10% leads to 9.8% less N 2 O emissions and 9.6% less N leaching at the cost of 4.9% more SOC depletion and 0.6% yield reduction over the study region. The estimated total annual NSB is $\$395$ M (uncertainty ranges from $\$114$ M to $\$1271$ M), including $\$334$ from social benefits (uncertainty ranges from $\$46$ M to $\$1076$ M), $\$100$ M from saving fertilizer (uncertainty ranges from $\$13$ M to $\$455$ M), and –$\$40$ M due to yield changes (uncertainty ranges from –$\$261$ M to $\$69$ M). For the median scenario, we noted that 20% of the study area accounts for nearly 50% of the NSB, and thus represent hotspot locations for targeted mitigation. Although the uncertainty range suggests that developing such a high-resolution framework is not yet settled and the scenario based estimations are not appropriate to inform the management practices for individual farmers, our efforts shed light on the new generation of analytical tools for life cycle assessment.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Viability, behavior, and color expression in the offspring of matings between common wall lizard Podarcis muralis color morphs

Abstract Color polymorphisms are widely studied to identify the mechanisms responsible for the origin and maintenance of phenotypic variability in nature. Two of the mechanisms of balancing selection currently thought to explain the long-term persistence of polymorphisms are the evolution of alternative phenotypic optima through correlational selection on suites of traits including color and heterosis. Both of these mechanisms can generate differences in offspring viability and fitness arising from different morph combinations. Here, we examined the effect of parental morph combination on fertilization success, embryonic viability, newborn quality, antipredator, and foraging behavior, as well as inter-annual survival by conducting controlled matings in a polymorphic lacertid Podarcis muralis, where color morphs are frequently assumed to reflect alternative phenotypic optima (e.g., alternative reproductive strategies). Juveniles were kept in outdoor tubs for a year in order to study inter-annual growth, survival, and morph inheritance. In agreement with a previous genome-wide association analysis, morph frequencies in the year-old juveniles matched the frequencies expected if orange and yellow expressions depended on recessive homozygosity at 2 separate loci. Our findings also agree with previous literature reporting higher reproductive output of heavy females and the higher overall viability of heavy newborn lizards, but we found no evidence for the existence of alternative breeding investment strategies in female morphs, or morph-combination effects on offspring viability and behavior. We conclude that inter-morph breeding remains entirely viable and genetic incompatibilities are of little significance for the maintenance of discrete color morphs in P. muralis from the Pyrenees.

Abalos, Javier↗