Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variational”

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

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

At least 91 records · Page 5

Solar Cycle Variations and Equatorial Oscillations: Modeling Study

Solar cycle activity effects (SCAE) in the lower and middle atmosphere, reported in several studies, are difficult to explain on the basis of the small changes in solar radiation that accompany the 11-year cycle, It is therefore natural to speculate that dynamical processes may come into play to produce a leverage. Such a leverage may be provided by the Quasi-Biennial Oscillation (QBO) in the zonal circulation of the stratosphere, which has been linked to solar activity variations. Driven primarily by wave mean flow interaction, the QBO period and its amplitude are variable but are also strongly influenced by the seasonal cycle in the solar radiation. This influence extends to low altitudes referred to as "downward control". Relatively small changes in solar radiative forcing can produce small changes in the period and phase of the QBO, but this in turn can produce measurable differences in the wind field. Thus, the QBO may be an amplifier of solar activity variations and a natural conduit of these variations to lower altitudes. To test this hypothesis, we conducted experiments with a 2D (two-dimensional) version of our Numerical Spectral Model that incorporates Hines' Doppler Spread Parameterization for small-scale gravity waves (GW). Solar cycle radiance variations (SCRV) are accounted for by changing the radiative heating rate on a logarithmic scale from 0.1 % at the surface to 1 % at 50 km to 10% at 100 km. With and without SCRV, but with the same GW flux, we then conduct numerical experiments to evaluate the magnitude of the SCAE in the zonal circulation. The numerical results indicate that, under certain conditions, the SCAE is significant and can extend to lower altitudes where the SCRV is inconsequential. At 20-km the differences in the modeled wind velocities are as large as 5 m/s. For a modeled QBO period of 30 months, we find that the seasonal cycle in the solar forcing (through the Semi-annual Oscillation (SAO)) acts as a strong pacemaker to lockup the phase and period of the QBO. The SCAE then shows up primarily as a distinct but relatively weak amplitude modulation. But with the QBO period between 30 and 34 (or less than 30, presumably) months, the seasonal phase lock is weak. Solar flux radiance variations in the seasonal cycle then cause variations in the QBO period and phase that amplify the SCAE to produce relatively large variations in the wind field. These variations also extend to mid latitudes.

Mayr, H. G.↗

Tactical Defenses Against Systematic Variation in Wind Tunnel Testing

This paper examines the role of unexplained systematic variation on the reproducibility of wind tunnel test results. Sample means and variances estimated in the presence of systematic variations are shown to be susceptible to bias errors that are generally non-reproducible functions of those variations. Unless certain precautions are taken to defend against the effects of systematic variation, it is shown that experimental results can be difficult to duplicate and of dubious value for predicting system response with the highest precision or accuracy that could otherwise be achieved. Results are reported from an experiment designed to estimate how frequently systematic variations are in play in a representative wind tunnel experiment. These results suggest that significant systematic variation occurs frequently enough to cast doubts on the common assumption that sample observations can be reliably assumed to be independent. The consequences of ignoring correlation among observations induced by systematic variation are considered in some detail. Experimental tactics are described that defend against systematic variation. The effectiveness of these tactics is illustrated through computational experiments and real wind tunnel experimental results. Some tutorial information describes how to analyze experimental results that have been obtained using such quality assurance tactics.

DeLoach, Richard↗

Summary of 2006 to 2010 FPMU Measurements of International Space Station Frame Potential Variations

Electric potential variations on the International Space Station (ISS) structure in low Earth orbit are dominated by contributions from interactions of the United States (US) 160 volt solar arrays with the relatively high density, low temperature plasma environment and inductive potentials generated by motion of the large vehicle across the Earth?s magnetic field. The Floating Potential Measurement Unit (FPMU) instrument suite comprising two Langmuir probes, a plasma impedance probe, and a floating potential probe was deployed in August 2006 for use in characterizing variations in ISS potential, the state of the ionosphere along the ISS orbit and its effect on ISS charging, evaluating effects of payloads and visiting vehicles, and for supporting ISS plasma hazard assessments. This presentation summarizes observations of ISS frame potential variations obtained from the FPMU from deployment in 2006 through the current time. We first describe ISS potential variations due to current collection by solar arrays in the day time sector of the orbit including eclipse exit and entry charging events, potential variations due to plasma environment variations in the equatorial anomaly, and visiting vehicles docked to the ISS structure. Next, we discuss potential variations due to inductive electric fields generated by motion of the vehicle across the geomagnetic field and the effects of external electric fields in the ionosphere. Examples of night time potential variations at high latitudes and their possible relationship to auroral charging are described and, finally, we demonstrate effects on the ISS potential due to European Space Agency and US plasma contactor devices.

Minow, Joseph I.↗

Solar UV Variations During the Decline of Cycle 23

Characterization of temporal and spectral variations in solar ultraviolet irradiance over a solar cycle is essential for understanding the forcing of Earth's atmosphere and climate. Satellite measurements of solar UV variability for solar cycles 21, 22, and 23 show consistent solar cycle irradiance changes at key wavelengths (e.g. 205 nm, 250 nm) within instrumental uncertainties. All historical data sets also show the same relative spectral dependence for both short-term (rotational) and long-term (solar cycle) variations. Empirical solar irradiance models also produce long-term solar UV variations that agree well with observational data. Recent UV irradiance data from the Solar Radiation and Climate Experiment (SORCE) Spectral Irradiance Monitor (SIM) and Solar Stellar Irradiance Comparison Experiment (SOLSTICE) instruments covering the declining phase of Cycle 23 present a different picture oflong-term solar variations from previous results. Time series of SIM and SOLSTICE spectral irradiance data between 2003 and 2007 show solar variations that greatly exceed both previous measurements and predicted irradiance changes over this period, and the spectral dependence of the SIM and SOLSTICE variations during these years do not show features expected from solar physics theory. The use of SORCE irradiance variations in atmospheric models yields substantially different middle atmosphere ozone responses in both magnitude and vertical structure. However, short-term solar variability derived from SIM and SOLSTICE UV irradiance data is consistent with concurrent solar UV measurements from other instruments, as well as previous results, suggesting no change in solar physics. Our analysis of short-term solar variability is much less sensitive to residual instrument response changes than the observations of long-term variations. The SORCE long-term UV results can be explained by under-correction of instrument response changes during the first few years of measurements, rather than requiring an unexpected change in the physical behavior of the Sun.

DeLand, Matthew, T.↗

Thermal Phase Variations of WASP-12b: Defying Predictions

We report Warm Spitzer full-orbit phase observations of WASP-12b at 3.6 and 4.5 micrometers. This extremely inflated hot Jupiter is thought to be overflowing its Roche lobe, undergoing mass loss and accretion onto its host star, and has been claimed to have a C/O ratio in excess of unity. We are able to measure the transit depths, eclipse depths, thermal and ellipsoidal phase variations at both wavelengths. The large-amplitude phase variations, combined with the planet's previously measured dayside spectral energy distribution, are indicative of non-zero Bond albedo and very poor day-night heat redistribution. The transit depths in the mid-infrared-(R(sub p)/R(sub ∗))(sup 2) = 0.0123(3) and 0.0111(3) at 3.6 and 4.5 micrometers, respectively-indicate that the atmospheric opacity is greater at 3.6 than at 4.5 micrometers, in disagreement with model predictions, irrespective of C/O ratio. The secondary eclipse depths are consistent with previous studies: F(sub day)/F(sub ∗) = 0.0038(4) and 0.0039(3) at 3.6 and 4.5 micrometers, respectively. We do not detect ellipsoidal variations at 3.6 micrometers, but our parameter uncertainties-estimated via prayer-bead Monte Carlo-keep this non-detection consistent with model predictions. At 4.5 micrometers, on the other hand, we detect ellipsoidal variations that are much stronger than predicted. If interpreted as a geometric effect due to the planet's elongated shape, these variations imply a 3:2 ratio for the planet's longest:shortest axes and a relatively bright day-night terminator. If we instead presume that the 4.5 micrometer ellipsoidal variations are due to uncorrected systematic noise and we fix the amplitude of the variations to zero, the best-fit 4.5 micrometer transit depth becomes commensurate with the 3.6 micrometer depth, within the uncertainties. The relative transit depths are then consistent with a solar composition and short scale height at the terminator. Assuming zero ellipsoidal variations also yields a much deeper 4.5 micrometer eclipse depth, consistent with a solar composition and modest temperature inversion. We suggest future observations that could distinguish between these two scenarios.

planetary↗

The Variable Crab Nebula: Evidence for a Connection Between GeV Flares and Hard X-ray Variations

In 2010, hard X-ray variations (Wilson-Hodge et al. 2011) and GeV flares (Tavani et al 2011, Abdo et al. 2011) from the Crab Nebula were discovered. Connections between these two phenomena were unclear, in part because the timescales were quite different, with yearly variations in hard X-rays and hourly to daily variations in the GeV flares. The hard X-ray flux from the Crab Nebula has again declined since 2014, much like it did in 2008-2010. During both hard X-ray decline periods, the Fermi LAT detected no GeV flares, suggesting that injection of particles from the GeV flares produces the much slower and weaker hard X-ray variations. The timescale for the particles emitting the GeV flares to lose enough energy to emit synchrotron photons in hard X-rays is consistent with the yearly variations observed in hard X-rays and with the expectation that the timescale for variations slowly increases with decreasing energy. This hypothesis also predicts even slower and weaker variations below 10 keV, consistent with the non-detection of counterparts to the GeV flares by Chandra (Weisskopf et al 2013). We will present a comparison of the observed hard X-ray variations and a simple model of the decay of particles from the GeV flares to test our hypothesis.

Variable Crab Nebula↗

Robust scalable initialization for Bayesian variational inference with multi-modal Laplace approximations

Predictive modeling typically relies on Bayesian model calibration to provide uncertainty quantification. Variational inference utilizing fully independent (“mean-field”) Gaussian distributions are often used as approximate probability density functions. This simplification is attractive since the number of variational parameters grows only linearly with the number of unknown model parameters. However, the resulting diagonal covariance structure and unimodal behavior can be too restrictive to provide useful approximations of intractable Bayesian posteriors that exhibit highly non-Gaussian behavior, including multimodality. High-fidelity surrogate posteriors for these problems can be obtained by considering the family of Gaussian mixtures. Gaussian mixtures are capable of capturing multiple modes and approximating any distribution to an arbitrary degree of accuracy, while maintaining some analytical tractability. Unfortunately, variational inference using Gaussian mixtures with full-covariance structures suffers from a quadratic growth in variational parameters with the number of model parameters. The existence of multiple local minima due to strong nonconvex trends in the loss functions often associated with variational inference present additional complications, These challenges motivate the need for robust initialization procedures to improve the performance and computational scalability of variational inference with mixture models. In this work, we propose a method for constructing an initial Gaussian mixture model approximation that can be used to warm-start the iterative solvers for variational inference. The procedure begins with a global optimization stage in model parameter space. In this step, local gradient-based optimization, globalized through multistart, is used to determine a set of local maxima, which we take to approximate the mixture component centers. Around each mode, a local Gaussian approximation is constructed via the Laplace approximation. Finally, the mixture weights are determined through constrained least squares regression. The robustness and scalability of the proposed methodology is demonstrated through application to an ensemble of synthetic tests using high-dimensional, multimodal probability density functions. Here, the practical aspects of the approach are demonstrated with inversion problems in structural dynamics.

97 MATHEMATICS AND COMPUTING↗

Variational nonlinear WKB in the Eulerian frame

Nonlinear WKB is a multiscale technique for studying locally plane-wave solutions of nonlinear partial differential equations (PDEs). Its application comprises two steps: (1) replacement of the original PDE with an extended system separating the large scales from the small and (2) reduction of the extended system to its slow manifold. In the context of variational fluid theories with particle relabeling symmetry, nonlinear WKB in the mean Eulerian frame is known to possess a variational structure. This much has been demonstrated using, for instance, the theoretical apparatus known as the generalized Lagrangian mean. On the other hand, the variational structure of nonlinear WKB in the conventional Eulerian frame remains mysterious. By exhibiting a variational principle for the extended equations from step (1) above, here we demonstrate that nonlinear WKB in the Eulerian frame is in fact variational. Remarkably, the variational principle for the extended system admits loops of relabeling transformations as a symmetry group. Noether’s theorem therefore implies that the extended Eulerian equations possess a family of circulation invariants parameterized by S 1 . As an illustrative example, we use our results to systematically deduce a variational model of high-frequency acoustic waves interacting with a larger-scale compressible isothermal flow.

97 MATHEMATICS AND COMPUTING↗

Gene-by-environment interactions influence the fitness cost of gene copy-number variation in yeast

Abstract Variation in gene copy number can alter gene expression and influence downstream phenotypes; thus copy-number variation provides a route for rapid evolution if the benefits outweigh the cost. We recently showed that genetic background significantly influences how yeast cells respond to gene overexpression, revealing that the fitness costs of copy-number variation can vary substantially with genetic background in a common-garden environment. But the interplay between copy-number variation tolerance and environment remains unexplored on a genomic scale. Here, we measured the tolerance to gene overexpression in four genetically distinct Saccharomyces cerevisiae strains grown under sodium chloride stress. Overexpressed genes that are commonly deleterious during sodium chloride stress recapitulated those commonly deleterious under standard conditions. However, sodium chloride stress uncovered novel differences in strain responses to gene overexpression. West African strain NCYC3290 and North American oak isolate YPS128 are more sensitive to sodium chloride stress than vineyard BC187 and laboratory strain BY4743. Consistently, NCYC3290 and YPS128 showed the greatest sensitivities to overexpression of specific genes. Although most genes were deleterious, hundreds were beneficial when overexpressed—remarkably, most of these effects were strain specific. Few beneficial genes were shared between the sodium chloride-sensitive isolates, implicating mechanistic differences behind their sodium chloride sensitivity. Transcriptomic analysis suggested underlying vulnerabilities and tolerances across strains, and pointed to natural copy-number variation of a sodium export pump that likely contributes to strain-specific responses to overexpression of other genes. Our results reveal extensive strain-by-environment interactions in the response to gene copy-number variation, raising important implications for the accessibility of copy-number variation-dependent evolutionary routes under times of stress.

60 APPLIED LIFE SCIENCES↗

Adaptive Variational Quantum Imaginary Time Evolution Approach for Ground State Preparation

Abstract An adaptive variational quantum imaginary time evolution (AVQITE) approach is introduced that yields efficient representations of ground states for interacting Hamiltonians on near‐term quantum computers. It is based on McLachlan's variational principle applied to imaginary time evolution of variational wave functions. The variational parameters evolve deterministically according to equations of motions that minimize the difference to the exact imaginary time evolution, which is quantified by the McLachlan distance. Rather than working with a fixed variational ansatz, where the McLachlan distance is constrained by the quality of the ansatz, the AVQITE method iteratively expands the ansatz along the dynamical path to keep the McLachlan distance below a chosen threshold. This ensures the state is able to follow the quantum imaginary time evolution path in the system Hilbert space rather than in a restricted variational manifold set by a predefined fixed ansatz. AVQITE is used to prepare ground states of H 4 , H 2 O, and BeH 2 molecules, where it yields compact variational ansätze and ground state energies within chemical accuracy. Polynomial scaling of circuit depth with system size is shown through a set of AVQITE calculations of quantum spin models. Finally, quantum Lanczos calculations are demonstrated alongside AVQITE without additional quantum resource costs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The impact of mineral reactive surface area variation on simulated mineral reactions and reaction rates

Reactive transport modeling is an essential tool to simulate complex geochemical reactions in porous media that can impact formation properties including porosity and permeability. However, simulating these reactions is challenging due to uncertainties in model parameters, particularly mineral surface areas. Imaging has emerged as a powerful means of estimating model parameters including porosity and mineral abundance, accessibility and accessible surface area. However, these parameters, particularly mineral accessible surface area, vary with image resolution. This work aims to enhance understanding of the impact of image resolution and other means of estimating mineral reactive surface area on simulated mineral reactions and reaction rates. Mineral surface areas calculated from images with resolutions of 0.34 μm and 5.71 μm were used to simulate mineral reactions in the context of geologic CO 2 sequestration in the Paluxy formation at the continuum scale. Additional simulations were carried out using BET surface areas collected from the literature and geometric surface areas. Simulations were run for 7300 days and mineral volume fractions and effluent ion concentrations tracked and compared. Variations in mineral surface areas measured from images are within 1 order of magnitude and yield similar simulation results, indicating the impact of image resolution on simulated reactions and reaction rates is minimum for the resolutions and sample considered. In comparison, surface areas obtained from BET and geometric approaches are 1–4 orders of magnitude higher than image-obtained surface areas and result in greater simulated reaction rates and extents. Minerals with high reaction rates (calcite and siderite) are most impacted by surface area values at short times where simulated mineral volume fractions at longer times agree relatively well, even for simulations with several orders of magnitude variation in surface area. Phases with lower reaction rates, such as K-feldspar and muscovite, are predominantly impacted over longer times where variations in surface areas impact reaction extents and porosity evolution. Overall, variations in surface areas due to image resolution are small and result in little variation in simulated reactions and reaction rates while there are significant variations in simulation results when other surface area estimates are used. These variations, however, depend on the reactivity of the mineral phase where surface areas of fast-reacting phases largely impact simulations at short (hours to days) timescales and surface areas of slower-reactive phases impact longer term simulations (years).

58 GEOSCIENCES↗

A systems genetics approach to deciphering the effect of dosage variation on leaf morphology in Populus

Abstract Gene copy number variation is frequent in plant genomes of various species, but the impact of such gene dosage variation on morphological traits is poorly understood. We used a large population of Populus carrying genomically characterized insertions and deletions across the genome to systematically assay the effect of gene dosage variation on a suite of leaf morphology traits. A systems genetics approach was used to integrate insertion and deletion locations, leaf morphology phenotypes, gene expression, and transcriptional network data, to provide an overview of how gene dosage influences morphology. Dosage-sensitive genomic regions were identified that influenced individual or pleiotropic morphological traits. We also identified cis-expression quantitative trait loci (QTL) within these dosage QTL regions, a subset of which modulated trans-expression QTL as well. Integration of data types within a gene co-expression framework identified co-expressed gene modules that are dosage sensitive, enriched for dosage expression QTL, and associated with morphological traits. Functional description of these modules linked dosage-sensitive morphological variation to specific cellular processes, as well as candidate regulatory genes. Together, these results show that gene dosage variation can influence morphological variation through complex changes in gene expression, and suggest that frequently occurring gene dosage variation has the potential to likewise influence quantitative traits in nature.

Biochemistry & Molecular Biology↗

Test of nuclear decay rate variation due to an antineutrino flux

Unexplained variations of the decay rate parameter for weak interaction decays such as β ± decay, electron capture, as well as strong interaction α decay have been reported. Because these variations have been presented by a number of groups, at various locations, using various types of detectors, different isotopes, and over extended periods of time, some researchers have interpreted the source of these variations as not from ambient environmental factors such as temperature, pressure, and humidity but via an unexplained fundamental interaction. The state of decay rate parameter variations experiments is reviewed, and the reported results are placed into a common comparable context by defining a phenomenological cross section. After that we make decay parameter measurements as a function of time at the level of 10 -5 in the presences of an antineutrino flux with an on-off cycle time of ≈ 30 days. This level of precision requires a detailed understanding of both systematic and statistical errors; otherwise, systematic errors in the measurement may mimic the decay events of fundamental interactions. The experiment searched for variation of the $^{54}_{25}$Mn, e - capture decay rate parameter and $^{137}_{55}$Cs, β - decay rate parameter both to a level of precision of 1 part in ≈10 5 by comparing the difference between the decay rate in the presence of an antineutrino flux ≈ 3 x 10 12 $\bar{ν}$ cm -2 s -1 and no flux measurements. The experiment was located 6.53 m from the reactor core of the High Flux Isotope Reactor at the Oak Ridge National Laboratory. Two weak interaction decays, one via electron capture and the other via β - decay, were selected because the final state and the time reverse state each contain a neutrino and antineutrino, covering arguments that the antineutrino flux may interact differently or not at all in one of the cases. The $\gamma$ spectra from both decays were collected and analyzed independently. The measured variation in the decay rate parameters are found to be δλ/λ =(0.034 ± 1.38) x10 -5 for $^{54}_{25}$Mn and δλ/λ =(0.67±1.56) x10 -5 for $^{137}_{55}$Cs. These results are consistent with no measurable decay rate parameter variation due to an antineutrino flux, yielding a 68% confidence level upper limit sensitivity for $^{54}_{25}$Mn δλ/λ ≤ 1.31 x 10 -5 or σ ≤ 1.29 x10 -25 cm 2 in cross section and for $^{137}_{55}$Cs δλ/λ ≤ 2.23 x 10 -5 or σ ≤ 5.69 x10 -27 cm 2 . The cross-section upper limit obtained in these null or no observable effect measurements are ≈ 10 4 times more sensitive than past experiments reporting positive results in 54 Mn and 137 Cs. . .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Prospects for measuring the time variation of astrophysical neutrino sources at dark matter detectors

We study the prospects for measuring the time variation of solar and atmospheric neutrino fluxes at future large-scale xenon and argon dark matter detectors. For solar neutrinos, a yearly time variation arises from the eccentricity of Earth’s orbit and, for charged current interactions, from a smaller energy-dependent day-night variation due to flavor regeneration as neutrinos travel through Earth. For a 100-ton xenon detector running for ten years with a xenon-136 fraction of ≲ 0.1 % , in the electron recoil channel a time-variation amplitude of about 0.8% is detectable with a power of 90% and the level of significance of 10%. This is sufficient to detect time variation due to eccentricity, which has amplitude of ∼ 3 % . In the nuclear recoil channel, the detectable amplitude is about 10% under current detector resolution and efficiency conditions, and this generally reduces to about 1% for improved detector resolution and efficiency, the latter of which is sufficient to detect time variation due to eccentricity. Our analysis assumes both known and unknown periods. We provide scalings to determine the sensitivity to an arbitrary time-varying amplitude as a function of detector parameters. Identifying the time variation of the neutrino fluxes will be important for distinguishing neutrinos from dark matter signals and other detector-related backgrounds and extracting properties of neutrinos that can be uniquely studied in dark matter experiments. Published by the American Physical Society 2024

Zhuang, Yi (ORCID:0000000277138724)↗

Resource selection functions based on hierarchical generalized additive models provide new insights into individual animal variation and species distributions

Habitat selection studies are designed to generate predictions of species distributions or inference regarding general habitat associations and individual variation in habitat use. Such studies frequently involve either individually indexed locations gathered across limited spatial extents and analyzed using resource selection functions (RSFs) or spatially extensive locational data without individual resolution typically analyzed using species distribution models. Both analytical methodologies have certain desirable features, but analyses that combine individual- and population-level inference with flexible non-linear functions may provide improved predictions while accounting for individual variation. Here, we describe how RSFs can be fit using hierarchical generalized additive models (HGAMs) using widely available software, providing a means to explore individual variation in habitat associations and to generate species distribution maps. We used GPS tracking data from golden eagles Aquila chrysaetos from across eastern North America with four environmental predictors to generate monthly distribution models. We considered three model structures that assumed different amounts of individual variation in the functional relationship between predictors and habitat use and used k-fold cross-validation to compare model performance. Models accounting for individual variability in shape and smoothness of functional responses performed best. Eagles exhibited the least amount of individual variation in response to land cover variables during winter months, with most individuals more closely adhering to the population-level trend. During the summer months, eagles exhibited more substantial individual variation in shape and smoothness of the functional relationships, suggesting some need to account for individual variation in eagle habitat use for both inferential and predictive purposes, during this time of year. Because they allow users to blend flexible functions with random effects structures and are well-supported by a variety of software platforms, we believe that HGAMs provide a useful addition to the suite of analyses used for modeling habitat associations or predicting species distributions.

54 ENVIRONMENTAL SCIENCES↗

Arctic Shrub Expansion, Plant Functional Trait Variation, and Effects on Belowground Carbon Cycling (Final Technical Report)

Terrestrial ecosystems are undergoing dramatic changes in response to climate warming, and these changes are expected to feedback to the atmosphere, potentially altering the trajectory of future climate change. Feedbacks from Arctic ecosystems are a major concern because the Arctic is projected to warm significantly in the 21 st century and because >50% of global belowground organic carbon is stored in permafrost and overlying soils. Warming-driven release of this carbon could drastically increase atmospheric greenhouse gas concentrations and accelerate climate warming. Plant communities are also responding to warming, as evidenced by the widely documented increase in woody-shrub growth and “greening” across much of the Arctic tundra biome. This vegetation shift may offset or amplify warming by altering carbon cycling. The direction and magnitude of shrub effects remain highly uncertain, however, due to limited understanding of the consequences of shrub expansion for belowground carbon cycling and simplification of these relationships in models. The major shrubs expanding in the Arctic (Betula, Salix, and Alnus) vary widely with respect to aboveground and belowground traits (e.g., tissue production and chemistry, rooting depth, microbial symbionts), and may also exhibit substantial intraspecific variation in these traits in response to environmental conditions. Such variation is likely to have profound implications for soil carbon cycling. The overarching goal of this project was to improve process-based understanding of the influence of shrub expansion on carbon cycling to enable improved representation of carbon dynamics in ecosystem and Earth system models. We investigated how plant functional traits vary among shrub genera, respond to environmental conditions, and affect belowground carbon and nutrient cycling by quantifying relationships among functional traits and biogeochemical cycling along edaphic gradients nested within a climate gradient in the Alaskan tundra. We found consistent differences in leaf and root traits among shrub genera and between shrubs and a widespread sedge species, indicating diverse nutrient acquisition strategies and belowground impacts among different arctic shrubs. We also found striking differences in trait values among individuals within the same species or genera within sites. Soil parameters were more important than climate parameters for predicting size and leaf trait variation, and root trait responses were less dependent on climate overall. For all but one root trait, including parameters representing aboveground traits improved the predictive ability of models. These results demonstrate that tundra shrub traits vary considerably at local scales and soil factors drive this variation, especially belowground. Furthermore, leveraging information about aboveground traits and soil conditions can improve predictions of how belowground traits will respond to climate change. Despite these differences, soil carbon and nitrogen pools in the active layer did not vary among plots dominated by different shrub or sedge genera. Instead, pool sizes generally decreased from warmer to colder sites, consistent with a productivity gradient. Patterns of isotopic N composition indicate that shrubs tighten nitrogen cycling via nitrogen resorption or immobilization of shrub litter. Overall, these results suggest that further identifying the specific shrub genera in the tundra landscape will ultimately provide better predictions of belowground dynamics across the changing arctic. We also performed simulation experiments with the Terrestrial Ecosystem Model (TEM) incorporated in the Predictive Ecosystem Analyzer (PEcAn) framework, treats model parameters as probability distributions, estimates parameters based on a synthesis of available field data, and then quantifies both model sensitivity and uncertainty to a given parameter or suite of parameters. We performed simulations across different types of tundra, including shrub tundra. One key finding was that both model sensitivity and uncertainty to a given parameter could vary within the same type of tundra, but in a different geographical location, such as over the climate gradient of shrub tundra described above. We organized a special session at the annual meeting of the Ecological Society of America in August 2019 to disseminate our results, refine recommendations for model improvement, and initiate collaborations to implement these recommendations in existing models of tundra carbon dynamics at ecosystem to Earth system scales. Our results support DOE near-term priorities by providing mechanistic insights into the role of vegetation change in the terrestrial carbon cycle in a region that is inadequately represented in Earth system models. Current models reduce the complexity of Arctic vegetation to a small number of plant functional types (PFTs). This approach implicitly assumes that each PFT represents the average ecological function of its constituent species, thus ignoring the effects of trait variation on biogeochemical cycling and potentially leading to large uncertainty in the sign and magnitude of ecosystem feedbacks to climate. By quantifying variation of plant functional traits across broad gradients of climatic and edaphic conditions and elucidating the linkages of such variation with carbon and nutrient cycling, our results illustrate the need and create a foundation for further developing trait-based modeling approaches that allow the traits of PFTs to vary as a function of environmental conditions. These approaches should improve the capacity of simulation models to offer insights into ecosystem carbon dynamics associated with novel plant communities in a rapidly changing Arctic.

54 ENVIRONMENTAL SCIENCES↗

Variational quantum solver employing the PDS energy functional

In our previous work (J. Chem. Phys. 2020, 153, 201102),we reported a new class of quantum algorithms that are based on the quantum computation of the connected moment expansion to find the ground and excited state energies. In particular, the Peeters-Devreese-Soldatov (PDS) formulation is found variational and bearing the potential for further combining with the existing variational quantum infrastructure. Following this direction, here we propose a variational quantum solver employing the PDS energy gradient. In comparison with the usual variational quantum eigensolver (VQE) and the original static PDS approach, the proposed variational quantum solver offers an effective approach to achieve high ac-curacy at finding the ground state and its energy through the rotation of the trial wave function of modest quality guided by the low order PDS energy gradients, thus improves the ac-curacy and efficiency of the quantum simulation. We demonstrate the performance of the proposed variational quantum solver for toy models, H2molecule, and strongly correlated planar H4system in some challenging situations. In all the case studies, the proposed variational quantum approach out-performs the usual VQE and static PDS calculations even at the lowest order.

Peng, Bo↗

Variation with interplanetary sector of the total magnetic field measured at the OGO 2, 4, and 6 satellites

Variations in the scalar magnetic field (delta B) from the polar orbiting OGO 2, 4, and 6 spacecraft are examined as a function of altitude for times when the interplanetary magnetic field is toward the sun and for times when the interplanetary magnetic field away from the sun. This morphology is basically the same as that found when all data, irrespective of interplanetary magnetic sector, are averaged together. Differences in delta B occur, both between sectors and between seasons, which are similar in nature to variations in the surface delta Z found by Langel (1973c). The altitude variation of delta B at sunlit local times, together with delta Z at the earth's surface, demonstrates that the delta Z and delta B which varies with sector has an ionospheric source. Langel (1973b) showed that the positive delta B region in the dark portion of the hemisphere is due to at least two sources, the westward electrojet and an unidentified non-ionospheric source(s). Comparison of magnetic variations between season/sector at the surface and at the satellite, in the dark portion of the hemisphere, indicates that these variations are caused by variations in the latitudinally narrow electrojet currents and not by variations in the non-ionospheric source of delta B.

Langel, R. A.↗