Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational inference”

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

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

At least 1,009 records · Page 56

Causes of the F region winter anomaly

Measurements made by the Atmosphere Explorer satellites are used to determine the causes of the F-region winter anomaly, i.e., seasonal variations which result in peak electron densities which are larger in winter than in summer. It is found that for the data examined there are three causes of the winter anomaly: (1) seasonal changes in neutral composition; (2) an increase in the vibrational temperature of N2 and hence in the rate coefficient for the reaction O(+)(4S) + N2 yields NO(+) + N; (3) an increase in winter in the production of O(+)(4S) due to enhanced electron quenching of O(+)(2D). The presence of vibrational excitation is inferred from an increase in the rate coefficient for the reaction O(+) + N2 yields NO(+) + N from winter to summer at heights above 260 km. It is demonstrated quantitatively that the initial increase in electron density produced in winter by causes (1) and (2) results in an increase in the quenching of O(+)(4S), thereby enhancing the source of O(+)(4S) ions. It is estimated that causes (1) to (3) account for approximately 55%, 20% and 25% respectively of the seasonal change in O(+)(4S) from summer to winter at solar minimum.

Torr, D. G.↗

Wind transport near the poles of Mars: Timescales of changes in deposition and erosion

Movement of sediment into and out of polar deposits is closely linked to the polar volatile budget and to changes in wind systems over the course of astronomically induced climate cycles. The present observations of the morphology of polar layered deposits, mantling sediments, dune fields, and variable surface features are the basis of inferences on the efficacy of polar sediment transport mechanisms. The time scales of formation of these features vary from days to perhaps 10(exp 6) yr, and latitudinal banding of dune fields near the poles may have been formed on timescales of 10(exp 7) yr. Orientations of intracrater dunes, dune crests, and wind streaks have been measured for latitudes -45 to -90 to compare features of likely different timescales of formation with models of wind flow from the south polar region. There is a complex variation with latitude of the indicated wind directions and of the efficacy of the resultant winds in orienting dune fields that suggests influence of frost cover on the ability of winds to move sediment in the spring and fall. Because of changes in the relative effectiveness of spring and fall winds expected with progression of the season of perihelion, the latitudinal variation in transport efficiency may mean that sediments at different latitudes dominantly respond to wind erosion and transport at different times during the perihelion cycle.

Thomas, Peter C.↗

High sensitivity dynamic spectral search for flare star radio

We observed ten well-known flare stars with the Arcibo radio telescope at 1.4 GHz and 5 GHz, using a special observing technique to discriminate between real flares and radio freqeuncy interference. With a high sensitivity of 5.5 K/Jy at 1.4 GHz when averaged over a 50 MHz band, we are able to recognize flux enhancements as weak as approximately 6 mJy above the sky background variations. In about 85 hours of observation, about a dozen bursts were detected, only from AD Leo. All had flux densities lower than 70 mJy, which probably explains their lack of fine structures (except for the strongest one), such as were reported in the literature for stronger flares. Half of the bursts that we recorded are 100% circularly polarized, and half are not circularly polarized. Our results are a first attempt of reliable statistics on dMe flare rates at 1.4 GHz. The high brightness temperatures we infer for the observed bursts are interpreted in terms of coherent emission processes, either the cyclotron maser instability or plasma radiation. Efficiencies are comparable to those of solar or planetary radio emissions in the case of the cyclotron maser, and higher than the solar efficiency in the case of plasma radiation, with the caveat that there are great uncertainties in the coronal model and the source size.

Abada-Simon, M.↗

The Distribution of Methane in Titan's Stratosphere from Cassini/CIRS Observations

Cassini/CIRS spectra in the far- and mid-infrared region are used to determine the abundance of methanein Titans lower stratosphere and investigate its distribution with latitude. The CIRS spectra include emissionfrom both the CH4 v4 band at 7.7 micron and pure rotational lines longwards of 50 micron, which show differential sensitivities to thermal profile and methane mole fraction. We analyze nadir and limb data taken over the first part of the Cassini mission (August 2005 to June 2010), including a selection of 12 latitudes that provides a reasonably complete and regular sampling of both hemispheres. Unexpectedly, but in a consistent manner for limb and nadir geometries, large variations of the methane mole fraction near 15 mbar (approx. 85 km) are found, with values ranging from approx. 1.0% (at low latitudes and near +/-50-55 deg) to approx. 1.5% (at +/-30-35 deg and polar latitudes). Error bars on the retrieved methane mole fraction are 0.07-0.12% at low latitudes in the Southern hemisphere and 0.14-0.21% northward of 40 deg N. A 1.0% methanemole fraction at low latitudes permits us to reconcile the HASI-measured temperatures below 147 kmaltitude (2.7 mbar) with inferences from CIRS. The roughly hemispherically-symmetric distribution ofmethane gas is reminiscent of that observed or predicted for the tropospheric methane clouds, whichon a yearly-averaged basis, show preferential occurrences at tropical and polar latitudes. We speculatethat convective events at these latitudes result into local stratospheric methane enrichment, whichmay persist year-round due to dynamical mixing times in the lower stratosphere only moderately shorterthan a Titan year.

Methane↗

Comparing and Integrating Constraint Programming and Temporal Planning for Quantum Circuit Compilation

Recently, the makespan-minimization problem of compiling a general class of quantum algorithms into near-term quantum processors has been introduced to the AI community. The research demonstrated that temporal planning is a strong solution approach for the studied class of quantum circuit compilation (QCC) problems. In this paper, we explore the use of methods from operations research, specifically constraint programming (CP), as an alternative and complementary approach to temporal planning. We also extend previous work by introducing two new problem variations that incorporate important characteristics identified by the quantum computing community. We apply temporal planning and CP to the baseline and extended QCC problems as both stand-alone and hybrid approaches. The hybrid method uses solutions found by temporal planning to warm-start CP, leveraging the ability of temporal planning to find satisficing solutions to problems with a high degree of task optionality, an area that CP typically struggles with. These solutions are then used to seed the CP formulation which significantly benefits from inferred bounds on planning horizon and task counts provided by the warm-start. Our extensive empirical evaluation indicates that while stand-alone CP is not competitive with temporal planning, except for the smallest problems, CP in a hybrid setting is beneficial for all temporal planners in all problem classes.

quantum circuit compilation↗

Structure of axisymmetric mantle plumes

The structure of axisymmetric subsolidus thermal plumes in the earth's lower mantle is inferred from calculations of axisymmetric thermal plumes in an infinite Prandtl number fluid with thermally activated viscosity. The velocity and temperature distribution is determined for axisymmetric convection above a heated disk in an incompressible fluid cylinder 2,400 km in height and 1,200 km in diameter. Several calculations of plumes with heat transport in the range 100-400 GW, similar to the advective heat transport at the Hawaiian hotspot, are presented. Hotspot formation by plumes originating at the base of the mantle requires both large viscosity variations and a minimum heat transport.

Olson, Peter↗

The multifrequency spectral evolution of blazar 3C 345

The blazar 3C 345 underwent a 2.5 mag optical outburst between 1990 November and 1991 May. We have obtained 10 nearly simultaneous multifrequency spectra during the course of the outburst in order to study the multifrequency spectral variations of 3C 345 as a function of time. Although our observations were not sampled frequently enough to completely resolve the variations in every frequency band, the general rise and decline of the outburst were seen in the UV through radio with differing rise times. Simulations of an electron distribution injected into a tangled magnetic field show a relationship between frequency and characteristic timescale that was also observed in the radio variations of 3C 345. The two X-ray observations made during the monitoring period showed no evidence of variability. The multifrequency spectrum was modeled with two major components: a relativistic jet and a relativistic thermal accretion disk. Models calculated for each spectrum indicate that the outburst can be explained in terms of these models by varying the high-energy cutoff of the injected electron distribution in the jet model, while also varying the mass accretion rate in the disk model. There is marginal evidence that the inferred accretion rate varies with the jet luminosity.

Webb, J. R.↗

A Review of Evidence for High Life Coefficients on Propeller and Rotor Blades Under Static Thrust Conditions with Some New Experimental Results

Interest has increased recently in the thrust-producing capability of rotors at very high collective pitch angles. An early reference noted this behaviour in rotors and offered alternative models for section lift characteristics to explain it. The same phenomenon was coincidentally noted and used in a propeller code, resulting in very good correlation with static thrust data. The proposed paper will present experimental data demonstrating the pronounced persistence of thrust for propellers at increasing collective pitch angles. Comparisons with blade element/momentum theory will be made. These results are expected to point to the need to define (ultimately to explain) aerodynamic lift and drag behaviour in a rotating environment. Experimental measurements made by the U.S. Army Aeroflightdynamics Directorate at the Ames Research Center have shown that locally measured normal force coefficients along the span of a highly twisted rotor blade continue to increase at high values of collective pitch. In some cases these coefficients exceed expected values for the same type of airfoil tested under two dimensional conditions. To date no one to the authors' knowledge has defined the variation of C(n) with pitch for very high angles (to 45 deg) in a rotating environment and for a blade of reasonably high aspect ratio; however, total propeller thrust measurements support the idea that stalling does not occur in the same way as on a wing. This paper will present experimental data in the form of surface pressure distributions as well as flow visualization (microtufts) to explore the aerodynamic behavior of the rotating airfoil at high values of blade incidence. This paper also reviews experimental evidence and infers some high lift coefficient behavior from it. Comparisons between predicted thrust, utilizing modified airfoil characteristics and a blade element model, and measured thrust for both rotors and propellers that cover the extremes of collective pitch are shown and discussed.

Talbot, Peter D.↗

Phase Function Determination in Support of Orbital Debris Size Estimation

To recover the size of a space debris object from photometric measurements, it is necessary to determine its albedo and basic shape: if the albedo is known, the reflective area can be calculated; and if the shape is known, the shape and area taken together can be used to estimate a characteristic dimension. Albedo is typically determined by inferring the object s material type from filter photometry or spectroscopy and is not the subject of the present study. Object shape, on the other hand, can be revealed from a time-history of the object s brightness response. The most data-rich presentation is a continuous light-curve that records the object s brightness for an entire sensor pass, which could last for tens of minutes to several hours: from this one can see both short-term periodic behavior as well as brightness variations with phase angle. Light-curve interpretation, however, is more art than science and does not lend itself easily to automation; and the collection method, which requires single-object telescope dedication for long periods of time, is not well suited to debris survey conditions. So one is led to investigate how easily an object s brightness phase function, which can be constructed from the more survey-friendly point photometry, can be used to recover object shape. Such a recovery is usually attempted by comparing a phase-function curve constructed from an object s empirical brightness measurements to analytically-derived curves for basic shapes or shape combinations. There are two ways to accomplish this: a simple averaged brightness-versus phase curve assembled from the empirical data, or a more elaborate approach in which one is essentially calculating a brightness PDF for each phase angle bin (a technique explored in unpublished AFRL/RV research and in Ojakangas 2011); in each case the empirical curve is compared to analytical results for shapes of interest. The latter technique promises more discrimination power but requires more data; the former can be assembled in its essentials from fewer measurements but will be less definitive in its assignments. The goal of the present study is to evaluate both techniques under debris survey conditions to determine their relative performance and, additionally, to learn precisely how a survey should be conducted in order to maximize their performance. Because the distendedness of objects has more of an effect than their precise shape in calculating a characteristic dimension, one is interested in the techniques discrimination ability to distinguish between an elongated rectangular prism and a short rectangular prism or cube, or an elongated cylinder from a squat cylinder or sphere. Sensitivity studies using simulated data will be conducted to determine discrimination power for both techniques as a function of amount of data collected and range (and specific region) of phase angles sampled. Empirical GEODSS photometry data for distended objects (dead payloads with solar panels, rocket bodies) and compact objects (cubesats, calibration spheres, squat payloads) will also be used to test this discrimination ability. The result will be a recommended technique and data collection paradigm for debris surveys in order to maximize this type of discrimination.

Hejduk, M. D.↗

Variability of Basal Melt Beneath the Pine Island Glacier Ice Shelf, West Antarctica

Observations from satellite and airborne platforms are combined with model calculations to infer the nature and efficiency of basal melting of the Pine Island Glacier ice shelf, West Antarctica, by ocean waters. Satellite imagery shows surface features that suggest ice-shelf-wide changes to the ocean s influence on the ice shelf as the grounding line retreated. Longitudinal profiles of ice surface and bottom elevations are analyzed to reveal a spatially dependent pattern of basal melt with an annual melt flux of 40.5 Gt/a. One profile captures a persistent set of surface waves that correlates with quasi-annual variations of atmospheric forcing of Amundsen Sea circulation patterns, establishing a direct connection between atmospheric variability and sub-ice-shelf melting. Ice surface troughs are hydrostatically compensated by ice-bottom voids up to 150m deep. Voids form dynamically at the grounding line, triggered by enhanced melting when warmer-than-average water arrives. Subsequent enlargement of the voids is thermally inefficient (4% or less) compared with an overall melting efficiency beneath the ice shelf of 22%. Residual warm water is believed to cause three persistent polynyas at the ice-shelf front seen in Landsat imagery. Landsat thermal imagery confirms the occurrence of warm water at the same locations.

Bindschadler, Robert↗

General-Purpose Bayesian Tensor Learning With Automatic Rank Determination and Uncertainty Quantification

A major challenge in many machine learning tasks is that the model expressive power depends on model size. Low-rank tensor methods are an efficient tool for handling the curse of dimensionality in many large-scale machine learning models. The major challenges in training a tensor learning model include how to process the high-volume data, how to determine the tensor rank automatically, and how to estimate the uncertainty of the results. While existing tensor learning focuses on a specific task, this paper proposes a generic Bayesian framework that can be employed to solve a broad class of tensor learning problems such as tensor completion, tensor regression, and tensorized neural networks. We develop a low-rank tensor prior for automatic rank determination in nonlinear problems. Our method is implemented with both stochastic gradient Hamiltonian Monte Carlo (SGHMC) and Stein Variational Gradient Descent (SVGD). We compare the automatic rank determination and uncertainty quantification of these two solvers. We demonstrate that our proposed method can determine the tensor rank automatically and can quantify the uncertainty of the obtained results. We validate our framework on tensor completion tasks and tensorized neural network training tasks.

Bayesian inference↗

Spatial Variability in Column CO2 Inferred from High Resolution GEOS-5 Global Model Simulations: Implications for Remote Sensing and Inversions

Column CO2 observations from current and future remote sensing missions represent a major advancement in our understanding of the carbon cycle and are expected to help constrain source and sink distributions. However, data assimilation and inversion methods are challenged by the difference in scale of models and observations. OCO-2 footprints represent an area of several square kilometers while NASA s future ASCENDS lidar mission is likely to have an even smaller footprint. In contrast, the resolution of models used in global inversions are typically hundreds of kilometers wide and often cover areas that include combinations of land, ocean and coastal areas and areas of significant topographic, land cover, and population density variations. To improve understanding of scales of atmospheric CO2 variability and representativeness of satellite observations, we will present results from a global, 10-km simulation of meteorology and atmospheric CO2 distributions performed using NASA s GEOS-5 general circulation model. This resolution, typical of mesoscale atmospheric models, represents an order of magnitude increase in resolution over typical global simulations of atmospheric composition allowing new insight into small scale CO2 variations across a wide range of surface flux and meteorological conditions. The simulation includes high resolution flux datasets provided by NASA s Carbon Monitoring System Flux Pilot Project at half degree resolution that have been down-scaled to 10-km using remote sensing datasets. Probability distribution functions are calculated over larger areas more typical of global models (100-400 km) to characterize subgrid-scale variability in these models. Particular emphasis is placed on coastal regions and regions containing megacities and fires to evaluate the ability of coarse resolution models to represent these small scale features. Additionally, model output are sampled using averaging kernels characteristic of OCO-2 and ASCENDS measurement concepts to create realistic pseudo-datasets. Pseudo-data are averaged over coarse model grid cell areas to better understand the ability of measurements to characterize CO2 distributions and spatial gradients on both short (daily to weekly) and long (monthly to seasonal) time scales

Ott, L.↗

Population structure limits the use of genomic data for predicting phenotypes and managing genetic resources in forest trees

There is overwhelming evidence that forest trees are locally adapted to climate. Thus, genecological models based on population phenotypes have been used to measure local adaptation, infer genetic maladaptation to climate, and guide assisted migration. However, instead of phenotypes, there is increasing interest in using genomic data for gene resource management. We used whole-genome resequencing and common-garden experiments to understand the genetic architecture of adaptive traits in black cottonwood. We studied the potential of using genome-wide association studies (GWAS) and genomic prediction to detect causal loci, identify climate-adapted phenotypes, and inform gene resource management. We analyzed population structure by partitioning phenotypic and genomic (single-nucleotide polymorphism) variation among 840 genotypes collected from 91 stands along 16 rivers. Most phenotypic variation (60 to 81%) occurred among populations and was strongly associated with climate. Population phenotypes were predicted well using genomic data (e.g., predictive abilityr> 0.9) but almost as well using climate or geography (r> 0.8). In contrast, genomic prediction within populations was poor (r< 0.2). We identified many GWAS associations among populations, but most appeared to be spurious based on pooled within-population analyses. Hierarchical partitioning of linkage disequilibrium and haplotype sharing suggested that within-population genomic prediction and GWAS were poor because allele frequencies of causal loci and linked markers differed among populations. Given the urgent need to conserve natural populations and ecosystems, our results suggest that climate variables alone can be used to predict population phenotypes, delineate seed zones and deployment zones, and guide assisted migration.

Science & Technology - Other Topics↗

Suprathermal electron energy distribution within the dayside Venus ionosphere

The suprathermal electron energy distribution for the dayside ionosphere has been derived from data returned by the Pioneer-Venus orbiter retarding potential analyzer. The shape and magnitude of the spectrum are consistent with the assumption that solar EUV radiation is the only significant source. The magnitude of the spectrum and its variation with altitude suggest that significant vertical transport occurs, with the electrons being lost through the ionopause. In turn, significant vertical transport suggests that the effective vertical electron heat conductivity may be comparable to the field-free value. The heat input to the thermal electron gas from the measured suprathermal electron flux is too small by a factor of at least five to maintain the observed electron temperature profile if the electron thermal conductivity is assumed to be close to the field-free value. It is thus inferred that most of the heat is supplied by the solar wind.

Knudsen, W. C.↗

Theoretical properties of Omega-loops in the convective zone of the Sun. 1: Emerging bipolar magnetic regions

It is proposed that the observed 500 G intensity of the magnetic fields emerging through the surface of the Sun can be understood from the Bernoulli effect in the upwelling Omega-loops of magnetic field. It is also proposed that the inferred 10(exp 5) G azimuthal flux bundles below the base of the convective zone can be understood as a consequence of the large-scale buoyancy associated with the upwelling fluid in and around the rising Omega-loop. The process fits in naturally with the Babcock-Leighton form of the solar alpha-omega-dynamo. The emerging Omega-loop implies the coherence of the upwelling all the way from the bottom of the convective zone, enhancing the convective heat transport to account for the observed variation of the solar irradiance by about two parts in 10(exp 3).

Parker, E. N.↗

Beyond carbon flux partitioning: Carbon allocation and nonstructural carbon dynamics inferred from continuous fluxes

Carbon (C) allocation and nonstructural carbon (NSC) dynamics play essential roles in plant growth and survival under stress and disturbance. However, quantitative understanding of these processes remains limited. Here, in this work, we propose a framework where we connect commonly measured carbon cycle components (eddy covariance fluxes of canopy CO 2 exchange, soil CO 2 efflux, and allometry-based biomass and net primary production) by a simple mass balance model to derive ecosystem-level NSC dynamics (NSC i ), C translocation (dC i ), and the biomass production efficiency (BPE i ) in above- and belowground plant (i = agp and bgp) compartments. We applied this framework to two long-term monitored loblolly pine (Pinus taeda) plantations of different ages in North Carolina and characterized the variations of NSC and allocation in years under normal and drought conditions. The results indicated that the young stand did not have net NSC flux at the annual scale, whereas the mature stand stored a near-constant proportion of new assimilates as NSC every year under normal conditions, which was comparable in magnitude to new structural growth. Roots consumed NSC in drought and stored a significant amount of NSC post drought. The above- and belowground dC i and BPE i varied more from year to year in the young stand and approached a relatively stable pattern in the mature stand. The belowground BPE bgp differed the most between the young and mature stands and was most responsive to drought. With the internal C dynamics quantified, this framework may also improve biomass production estimation, which reveals the variations resulting from droughts. Overall, these quantified ecosystem-scale dynamics were consistent with existing evidence from tree-based manipulative experiments and measurements and demonstrated that combining the continuous fluxes as proposed here can provide additional information about plant internal C dynamics. Given that it is based on broadly available flux data, the proposed framework is promising to improve the allocation algorithms in ecosystem C cycle models and offers new insights into observed variability in soil–plant–climate interactions.

54 ENVIRONMENTAL SCIENCES↗

Ecosystem Water‐Saving Timescale Varies Spatially With Typical Drydown Length

Abstract Stomatal optimization theory is a commonly used framework for modeling how plants regulate transpiration in response to the environment. Most stomatal optimization models assume that plants instantaneously optimize a reward function such as carbon gain. However, plants are expected to optimize over longer timescales given the rapid environmental variability they encounter. There are currently no observational constraints on these timescales. Here, a new stomatal model is developed and is used to analyze the timescales over which stomatal closure is optimized. The proposed model assumes plants maximize carbon gain subject to the constraint that they cannot draw down soil moisture below a critical value. The reward is integrated over time, after being weighted by a discount factor that represents the timescale ( τ ) that a plant considers when optimizing stomatal conductance to save water. The model is simple enough to be analytically solvable, which allows the value of τ to be inferred from observations of stomatal behavior under known environmental conditions. The model is fitted to eddy covariance data in a range of ecosystems, finding the value of τ that best predicts the dynamics of evapotranspiration at each site. Across 82 sites, the climate metrics with the strongest correlation to τ are measures of the average number of dry days between rainfall events. Values of τ are similar in magnitude to the longest such dry period encountered in an average year. The results here shed light on which climate characteristics shape spatial variations in ecosystem‐level water use strategy.

54 ENVIRONMENTAL SCIENCES↗

Using Machine Learning to Generate a GISS ModelE Calibrated Physics Ensemble (CPE)

A neural network (NN) surrogate of the NASA GISS ModelE atmosphere (version E3) is trained on a perturbed parameter ensemble (PPE) spanning 45 physics parameters and 36 outputs. The NN is leveraged in a Markov Chain Monte Carlo (MCMC) Bayesian parameter inference framework to generate a second posterior constrained ensemble coined a “calibrated physics ensemble,” or CPE. The CPE members are characterized by diverse parameter combinations and are, by definition, close to top-of-atmosphere radiative balance, and must broadly agree with numerous hydrologic, energy cycle and radiative forcing metrics simultaneously. Global observations of numerous cloud, environment, and radiation properties (provided by global satellite products) are crucial for CPE generation. The inference framework explicitly accounts for discrepancies (or biases) in satellite products during CPE generation. We demonstrate that product discrepancies strongly impact calibration of important model parameter settings (e.g., convective plume entrainment rates; fall speed for cloud ice). Structural improvements new to E3 are retained across CPE members (e.g., stratocumulus simulation). Notably, the framework improved the simulation of shallow cumulus and Amazon rainfall while not degrading radiation fields, an upgrade that neither default parameters nor Latin Hypercube parameter searching achieved. Analyses of the initial PPE suggested several parameters were unimportant for output variation. However, many “unimportant” parameters were needed for CPE generation, a result that brings to the forefront how parameter importance should be determined in PPEs. From the CPE, two diverse 45-dimensional parameter configurations are retained to generate radiatively-balanced, auto-tuned atmospheres that were used in two E3 submissions to CMIP6.

54 ENVIRONMENTAL SCIENCES↗