Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “random networks”

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 559 records · Page 31

Reliability of Open Public Electric Vehicle Direct Current Fast Chargers

The aim was to systematically evaluate the usability of all public electric vehicles (EV) direct current fast chargers (DCFC) in the San Francisco region. To achieve a rapid transition to EVs, a highly reliable and easy to use charging infrastructure is critical to building confidence among consumers. The functionality and usability of all 182 open, public DCFC charging stations with CCS connectors (combined charging system) in the 9 counties of the Bay Area were tested (655 electric vehicle service equipment (EVSE) ports). An EVSE was classified as functional if it charged an EV for 2 minutes. Overall, 73.3% of the 655 EVSEs were functional. The causes of the nonfunctioning EVSEs (23.5%) were blank or unresponsive screens or error messages; payment system failures; charge initiation failures; network failures; or broken connectors. In addition, the cable was too short to reach the EV inlet for 3.2% of the EVSEs. A random sampling of 10% of the EVSEs, approximately 8 days after the first evaluation, found no overall change in functionality. The level of functionality found with field testing conflicts with the 95–98% uptime reported by the EV service providers (EVSPs) who operate the EV charging stations. There is a need for precise and verifiable definitions of uptime, downtime, and excluded time, as applied to public EV chargers. In conclusion, the level of failure of the existing public EV DCFC charge infrastructure highlights the importance of improving the system design and maintenance to improve adoption of EVs.

33 ADVANCED PROPULSION SYSTEMS↗

Extended observability of linear time-invariant systems under recurrent loss of output data

Recurrent loss of sensor data in integrated control systems of an advanced aircraft may occur under different operating conditions that include detected frame errors and queue saturation in computer networks, and bad data suppression in signal processing. This paper presents an extension of the concept of observability based on a set of randomly selected nonconsecutive outputs in finite-dimensional, linear, time-invariant systems. Conditions for testing extended observability have been established.

Luck, Rogelio↗

Monitoring global geodynamical parameter variations with GPS and low earth satellites

A low earth satellite enhances global geodynamical parameters determination with GPS in two ways. First, it improves the GPS orbits, which in turn improve the estimates of other parameters. Secondly, a low earth satellite completes an orbit cycle in far shorter time (90 to 120 minutes) than do GPS satellites (12 hours); it observes more GPS satellites than a ground receiver does in shorter time and increases the correlation between GPS orbit errors. This reduces the error in the determination of nonrotational coordinate parameters, i.e., geocentric offset of ground tracking sites. Covariance results with the global geodynamical parameters modeled as constants, and as random-walk parameters to closer reflect the actual variations, are compared. The effects of using different GPS data quality and different ground tracking network are studied. The case of using two earth satellites in orthogonal orbital planes is also investigated.

Wu, S. C.↗

Water Across Synthetic Aperture Radar Data (WASARD): SAR Water Body Classification for the Open Data Cube

The detection of inland water bodies from Synthetic Aperture Radar (SAR) data provides a great advantage over water detection with optical data, since SAR imaging is not impeded by cloud cover. Traditional methods of detecting water from SAR data involves using thresholding methods that can be labor intensive and imprecise. This paper describes Water Across Synthetic Aperture Radar Data (WASARD): a method of water detection from SAR data which automates and simplifies the thresholding process using machine learning on training data created from Geoscience Australia’s WOFS algorithm. Of the machine learning models tested, the Linear Support Vector Machine was determined to be optimal, with the option of training using solely the VH polarization or a combination of the VH and VV polarizations. WASARD was able to identify water in the target area with a correlation of 97% with WOFS. Sentinel-1, Open Data Cube, Earth Observations, Machine Learning, Water Detection 1. INTRODUCTION Water classification is an important function of Earth imaging satellites, as accurate remote classification of land and water can assist in land use analysis, flood prediction, climate change research, as well as a variety of agricultural applications [2]. The ability to identify bodies of water remotely via satellite is immensely cheaper than contracting surveys of the areas in question, meaning that an application that can accurately use satellite data towards this function can make valuable information available to nations which would not be able to afford it otherwise. Highly reliable applications for the remote detection of water currently exist for use with optical satellite data such as that provided by LANDSAT. One such application, Geoscience Australia’s Water Observations from Space (WOFS) has already been ported for use with the Open Data Cube [6]. However, water detection using optical data from Landsat is constrained by its relatively long revisit cycle of 16 days [5], and water detection using any optical data is constrained in that it lacks the ability to make accurate classifications through cloud cover [2]. The alternative solution which solves these problems is water detection using SAR data, which images the Earth using cloud-penetrating microwaves. Because of its advantages over optical data, much research has been done into water detection using SAR data. Traditionally, this has been done using the thresholding method, which involves picking a polarization band and labeling all pixels for which this band’s value is below a certain threshold as containing water. The thresholding method works since water tends to return a much lower backscatter value to the satellite than land [1]. However, this method can be flawed since estimating the proper threshold is often imprecise, complicated, and labor intensive for the end user. Thresholding also tends to use data from only one SAR polarization, when a combination of polarizations can provide insight into whether water is present. [2] In order to alleviate these problems, this paper presents an application for the Open Data Cube to detect water from SAR data using support vector machine (SVM) classification. 2. PLATFORM WASARD is an application for the Open Data Cube, a mechanism which provides a simple yet efficient means of ingesting, storing, and retrieving remote sensing data. Data can be ingested and made analysis ready according to whatever specifications the researcher chooses, and easily resampled to artificially alter a scene’s resolution. Currently WASARD supports water detection on scenes from ESA’s Sentinel-1 and JAXA’s ALOS. When testing WASARD, Sentinel-1 was most commonly used due to its relatively high spatial resolution and its rapid 6 day revisit cycle [5]. With minor alterations to the application's code, however, it could support data from other satellites. 3. METHODOLOGY Using supervised classification, WASARD compares SAR data to a dataset pre-classified by WOFS in order to train an SVM classifier. This classifier is then used to detect water in other SAR scenes outside the training set. Accuracy was measured according to the following metrics:  Precision: a measure of what percentage of the points WASARD labels as water are truly water  Recall: a measure of what percentage of the total water cover WASARD was able to identify.  F1 Score: a harmonic average of the precision and recall scores Both precision and recall are calculated at the end of the training phase, when the trained classifier is compared to a testing dataset. Because the WOFS algorithm’s classifications are used as the truth values when training a WASARD classifier, when precision and recall are mentioned in this paper, they are always with respect to the values produced by WOFS on a similar scene of Landsat data, which themselves have a classification accuracy of 97% [6]. Visual representations of water identified by WASARD in this paper were produced using the function wasard_plot(), which is included in WASARD. 3.1 Algorithm Selection The machine learning model used by WASARD is the Linear Support Vector Machine (SVM). This model uses a supervised learning algorithm to develop a classifier, meaning it creates a vector which can be multiplied by the vector formed by the relevant data bands to determine whether a pixel in a SAR scene contains water. This classifier is trained by comparing data points from selected bands in a SAR scene to their respective labels, which in this case are “water” or “not water” as given by the WOFS algorithm. The SVM was selected over the Random Forest model, which outperformed the SVM in training speed, but had a greater classification time and lower accuracy, and the Multilayer Perceptron Artificial Neural Network, which had a slightly higher average accuracy than the SVM, but much greater training and classification times. Figure 1: Visual representation of the SVM Classifier. Each white point represents a pixel in a SAR scene. In Figure 1, the diagonal line separating pixels determined to be water from those determined not to be water represents the actual classification vector produced by the SVM. It is worth noting that once the model has been trained, classification of pixels is done in a similar manner as in the thresholding method. This is especially true if only one band was used to train the model. 3.1 Feature Selection Sentinel-1 collects data from two bands: the Vertical/Vertical polarization (VV) and the Vertical/Horizontal polarization (VH). When 100 SVM classifiers were created for each polarization individually, and for the combination of the two, the following results were achieved: Figure 2: Accuracy of classifiers trained using different polarization bands. Precision and Recall were measured with respect to the values produced by WOFS. Figure 2 demonstrates that using both the VV and VH bands trades slightly lower recall for significantly greater precision when compared with the VH band alone, and that using the VV band alone is inferior in both metrics. WASARD therefore defaults to using both the VV and VH bands, and includes the option to use solely the VH band. The VV polarization’s lower precision compared to the VH polarization is in contrast to results from previous research and may merit further analysis [4]. 3.2 Training a Classifier The steps in training a classifier with WASARD are 1. Selecting two scenes (one SAR, one optical) with the same spatial extents, and acquired close to each other in time, with a preference that the scenes are taken on the same day. 2. Using the WOFS algorithm to produce an array of the detected water in the scene of optical data, to be used as the labels during supervised learning 3. Data points from the selected bands from the SAR acquisition are bundled together into an array with the corresponding labels gathered from WOFS. A random sample with an equal number of points labeled “Water” and “Not Water” is selected to be partitioned into a training and a testing dataset 4. Using Scikit-Learn’s LinearSVC object, the training dataset is used to produce a classifier, which is then tested against the testing dataset to determine its precision and recall The result is a wasard_classifier object, which has the following attributes: 1. f1, recall, and precision: 3 metrics used to determine the classifier’s accuracy 2. Coefficient: Vector which the SVM uses to make its predictions. The classifier detects water when the dot product of the coefficient and the vector formed by the SAR bands is positive 3. Save(): allows a user to save a classifier to the disk in order to use it without retraining 4. wasard_classify(): Classifies an entire xarray of SAR data using the SVM classifier All of the above steps are performed automatically when the user creates a wasard_classifier object. 3.3 Classifying a Dataset Once the classifier has been created, it can be used to detect water in an xarray of SAR data using wasard_classify(). By taking the dot product of the classifier’s coefficients and the vector formed by the selected bands of SAR data, an array of predictions is constructed. A classifier can effectively be used on the same spatial extents as the ones where it was trained, or on any area with a similar landscape. While

Kreiser, Zachary↗

Entropy Analysis of FPGA Interconnect and Switch Matrices for Physical Unclonable Functions

Random variations in microelectronic circuit structures represent the source of entropy for physical unclonable functions (PUFs). In this paper, we investigate delay variations that occur through the routing network and switch matrices of a field-programmable gate array (FPGA). The delay variations are isolated from other components of the programmable logic, e.g., look-up tables (LUTs), flip-flops (FFs), etc., using a feature of Xilinx FPGAs called dynamic partial reconfiguration (DPR). A set of partial designs is created to fix the placement of a time-to-digital converter (TDC) and supporting infrastructure to enable the path delays through the target interconnect and switch matrices to be extracted by subtracting out common-mode delay components. Delay variations are analyzed in the different levels of routing resources available within FPGAs, i.e., local routing and across-chip routing. Data are collected from a set of Xilinx Zynq 7010 devices, and a statistical analysis of within-die variations in delay through a set of the randomly-generated and hand-crafted interconnects is presented.

97 MATHEMATICS AND COMPUTING↗

Trends in lignin modification: a comprehensive analysis of the effects of genetic manipulations/mutations on lignification and vascular integrity

A comprehensive assessment of lignin configuration in transgenic and mutant plants is long overdue. This review thus undertook the systematic analysis of trends manifested through genetic and mutational manipulations of the various steps associated with monolignol biosynthesis; this included consideration of the downstream effects on organized lignin assembly in the various cell types, on vascular function/integrity, and on plant growth and development. As previously noted for dirigent protein (homologs), distinct and sophisticated monolignol forming metabolic networks were operative in various cell types, tissues and organs, and form the cell-specific guaiacyl (G) and guaiacyl-syringyl (G-S) enriched lignin biopolymers, respectively. Regardless of cell type undergoing lignification, carbon allocation to the different monolignol pools is apparently determined by a combination of phenylalanine availability and cinnamate-4-hydroxylase/"p-coumarate-3-hydroxylase" (C4H/C3H) activities, as revealed by transcriptional and metabolic profiling. Downregulation of either phenylalanine ammonia lyase or cinnamate-4-hydroxylase thus predictably results in reduced lignin levels and impaired vascular integrity, as well as affecting related (phenylpropanoid-dependent) metabolism. Depletion of C3H activity also results in reduced lignin deposition, albeit with the latter being derived only from hydroxyphenyl (H) units, due to both the guaiacyl (G) and syringyl (S) pathways being blocked. Apparently the cells affected are unable to compensate for reduced G/S levels by increasing the amounts of H-components. The downstream metabolic networks for G-lignin enriched formation in both angiosperms and gymnosperms utilize specific cinnamoyl CoA O-methyltransferase (CCOMT), 4-coumarate:CoA ligase (4CL), cinnamoyl CoA reductase (CCR) and cinnamyl alcohol dehydrogenase (CAD) isoforms: however, these steps neither affect carbon allocation nor H/G designations, this being determined by C4H/C3H activities. Such enzymes thus fulfill subsidiary processing roles, with all (except CCOMT) apparently being bifunctional for both H and G substrates. Their severe downregulation does, however, predictably result in impaired monolignol biosynthesis, reduced lignin deposition/vascular integrity, (upstream) metabolite build-up and/or shunt pathway metabolism. There was no evidence for an alternative acid/ester O-methyltransferase (AEOMT) being involved in lignin biosynthesis.The G/S lignin pathway networks are operative in specific cell types in angiosperms and employ two additional biosynthetic steps to afford the corresponding S components, i.e. through introduction of an hydroxyl group at C-5 and its subsequent O-methylation. [These enzymes were originally classified as ferulate-5-hydroxylase (F5H) and caffeate O-methyltransferase (COMT), respectively.] As before, neither step has apparently any role in carbon allocation to the pathway; hence their individual downregulation/manipulation, respectively, gives either a G enriched lignin or formation of the well-known S-deficient bm3 "lignin" mutant, with cell walls of impaired vascular integrity. In the latter case, COMT downregulation/mutation apparently results in utilization of the isoelectronic 5-hydroxyconiferyl alcohol species albeit in an unsuccessful attempt to form G-S lignin proper. However, there is apparently no effect on overall G content, thereby indicating that deposition of both G and S moieties in the G/S lignin forming cells are kept spatially, and presumably temporally, fully separate. Downregulation/mutation of further downstream steps in the G/S network [i.e. utilizing 4CL, CCR and CAD isoforms] gives predictable effects in terms of their subsidiary processing roles: while severe downregulation of 4CL gave phenotypes with impaired vascular integrity due to reduced monolignol supply, there was no evidence in support of increased growth and/or enhanced cellulose biosynthesis. CCR and CAD downregulation/mutations also established that a depletion in monolignol supply reduced both lignin contents supply reduced both lignin contents and vascular integrity, with a concomitant shift towards (upstream) metabolite build-up and/or shunting.The extraordinary claims of involvement of surrogate monomers (2-methoxybenzaldehyde, feruloyl tyramine, vanillic acid, etc.) in lignification were fully disproven and put to rest, with the investigators themselves having largely retracted former claims. Furthermore analysis of the well-known bm1 mutation, a presumed CAD disrupted system, apparently revealed that both G and S lignin components were reduced. This seems to imply that there is no monolignol specific dehydrogenase, such as the recently described sinapyl alcohol dehydrogenase (SAD) for sinapyl alcohol formation. Nevertheless, different CAD isoforms of differing homology seem to be operative in different lignifying cell types, thereby giving the G-enriched and G/S-enriched lignin biopolymers, respectively. For the G-lignin forming network, however, the CAD isoform is apparently catalytically less efficient with all three monolignols than that additionally associated with the corresponding G/S lignin forming network(s), which can more efficiently use all three monolignols. However, since CAD does not determine either H, G, or S designation, it again serves in a subsidiary role-albeit using different isoforms for different cell wall developmental and cell wall type responses.The results from this analysis contrasts further with speculations of some early investigators, who had viewed lignin assembly as resulting from non-specific oxidative coupling of monolignols and subsequent random polymerization. At that time, though, the study of the complex biological (biochemical) process of lignin assembly had begun without any of the (bio)chemical tools to either address or answer the questions posed as to how its formation might actually occur. Today, by contrast, there is growing recognition of both sophisticated and differential control of monolignol biosynthetic networks in different cell types, which serve to underscore the fact that complexity of assembly need not be confused any further with random formation. Moreover, this analysis revealed another factor which continues to cloud interpretations of lignin downregulation/mutational analyses, namely the serious technical problems associated with all aspects of lignin characterization, whether for lignin quantification, isolation of lignin-enriched preparations and/or in determining monomeric compositions. For example, in the latter analyses, some 50-90% of the lignin components still cannot be detected using current methodologies, e.g. by thioacidolysis cleavage and nitrobenzene oxidative cleavage. This deficiency in lignin characterization thus represents one of the major hurdles remaining in delineating how lignin assembly (in distinct cell types) and their configuration actually occurs.

Review, Academic↗

Parallel-in-Time Solution of Allen-Cahn Equations by Integrating Operator Learning into the Parareal Method

While recent advances in deep learning have shown promising efficiency gains in solving time-dependent partial differential equations (PDEs), matching the accuracy of conventional numerical solvers still remains a challenge. One strategy to improve the accuracy of deep learning-based solutions for time-dependent PDEs is to use the learned model as the coarse propagator in the Parareal method and a traditional numerical method as the fine solver. However, successful integration of deep learning into the Parareal method requires consistency between the coarse and fine solvers, particularly for PDEs exhibiting rapid changes such as sharp transitions. Here, to ensure this consistency, we propose using convolutional neural networks (CNNs) to learn the fully discrete time-stepping operator defined by the same numerical scheme employed as the fine solver. We demonstrate the effectiveness of the proposed method in solving the classical and mass-conservative Allen–Cahn (AC) equations. Through iterative updates in the Parareal algorithm, our approach achieves a significant computational speedup compared to traditional fine solvers while converging to high-accuracy solutions. Our results highlight that the proposed hybrid Parareal algorithm effectively accelerates simulations, particularly when implemented on multiple GPUs, and converges to the desired accuracy in only a few iterations. Another advantage of our method is that the CNN model is trained on trajectory-based data generated from random initial conditions, such that the trained model can be used to solve the AC equations with various initial conditions without retraining. This work demonstrates the potential of integrating neural network methods into parallel-in-time frameworks for efficient and accurate simulations of time-dependent PDEs.

97 MATHEMATICS AND COMPUTING↗

Islands Within Islands: Bacterial Phylogenetic Structure and Consortia in Hawaiian Lava Caves and Fumaroles

Lava caves, tubes, and fumaroles in Hawai‘i present a range of volcanic, oligotrophic environments from different lava flows and host unexpectedly high levels of bacterial diversity. These features provide an opportunity to study the ecological drivers that structure bacterial community diversity and assemblies in volcanic ecosystems and compare the older, more stable environments of lava tubes, to the more variable and extreme conditions of younger, geothermally active caves and fumaroles. Using 16S rRNA amplicon-based sequencing methods, we investigated the phylogenetic distinctness and diversity and identified microbial interactions and consortia through co-occurrence networks in 70 samples from lava tubes, geothermal lava caves, and fumaroles on the island of Hawai‘i. Our data illustrate that lava caves and geothermal sites harbor unique microbial communities, with very little overlap between caves or sites. We also found that older lava tubes (500–800 yrs old) hosted greater phylogenetic diversity (Faith's PD) than sites that were either geothermally active or younger (<400 yrs old). Geothermally active sites had a greater number of interactions and complexity than lava tubes. Average phylogenetic distinctness, a measure of the phylogenetic relatedness of a community, was higher than would be expected if communities were structured at random. This suggests that bacterial communities of Hawaiian volcanic environments are phylogenetically over-dispersed and that competitive exclusion is the main driver in structuring these communities. This was supported by network analyses that found that taxa (Class level) co-occurred with more distantly related organisms than close relatives, particularly in geothermal sites. Network “hubs” (taxa of potentially higher ecological importance) were not the most abundant taxa in either geothermal sites or lava tubes and were identified as unknown families or genera of the phyla, Chloroflexi and Acidobacteria. These results highlight the need for further study on the ecological role of microbes in caves through targeted culturing methods, metagenomics, and long-read sequence technologies.

59 BASIC BIOLOGICAL SCIENCES↗

Genetic programming for understanding cognitive biases that generate polarization in social networks

Recent studies have applied agent-based models to infer human-interpretable explanations of individual-scale behaviors that generate macro-scale patterns in complex social systems. Genetic programming has proven to be an ideal explainable AI tool for this purpose, where primitives may be expressed in an interpretable fashion and assembled into agent rules. Evolutionary model discovery (EMD) is a tool that combines genetic programming and random forest feature importance analysis, to infer individual-scale, human-interpretable explanations from agent-based models. We deploy EMD to investigate the cognitive biases behind the emergence of ideological polarization within a population. An agent-based model is developed to simulate a social network, where agents are able to create or sever links with one another, and update an internal ideological stance based on their neighbors' stances. Agent rules govern these actions and constitute of cognitive biases. A set of 7 cognitive biases are included as genetic program primitives in the search for rules that generate hyper-polarization among the population of agents. We find that heterogeneity in cognitive biases is more likely to generate polarized social networks. Highly polarized social networks are likely to emerge when individuals with confirmation bias are exposed to those with either attentional bias, egocentric bias, or cognitive dissonance.

Gunaratne, Chathika↗

A robust synthetic data generation framework for machine learning in high-resolution transmission electron microscopy (HRTEM)

Machine learning techniques are attractive options for developing highly-accurate analysis tools for nanomaterials characterization, including high-resolution transmission electron microscopy (HRTEM). However, successfully implementing such machine learning tools can be difficult due to the challenges in procuring sufficiently large, high-quality training datasets from experiments. In this work, we introduce Construction Zone, a Python package for rapid generation of complex nanoscale atomic structures which enables fast, systematic sampling of realistic nanomaterial structures and can be used as a random structure generator for large, diverse synthetic datasets. Using Construction Zone, we develop an end-to-end machine learning workflow for training neural network models to analyze experimental atomic resolution HRTEM images on the task of nanoparticle image segmentation purely with simulated databases. Further, we study the data curation process to understand how various aspects of the curated simulated data—including simulation fidelity, the distribution of atomic structures, and the distribution of imaging conditions—affect model performance across three benchmark experimental HRTEM image datasets. Using our workflow, we are able to achieve state-of-the-art segmentation performance on these experimental benchmarks and, further, we discuss robust strategies for consistently achieving high performance with machine learning in experimental settings using purely synthetic data. Construction Zone and its documentation are available at https://github.com/lerandc/construction_zone.

36 MATERIALS SCIENCE↗

Active Learning Framework

Machine learning (ML) of interatomic potentials show great promise to accelerate scientific simulation, e.g., by emulating expensive computations at a high accuracy but much reduced computational cost. Training datasets are calculated from computationally expensive ab initio quantum mechanics methods, density functional theory (DFT). Trained on this data, an ML model can be very successful in predicting energy and forces for new atomic configurations. A critical factor is the quality and diversity of the training dataset. Thus, a highly automated approach to dataset construction based on active learning framework is designed suitable for material physics. The active learning scheme begins with fully randomized atomic configurations. Then, many Molecular Dynamics (MD) trajectories are simulated using current ML potentials, where each MD trajectory is initialized to a random disordered configuration. The temperature is varied in order to diversify the sampled configuration during these simulations. The variance of predictions for eight neural networks within an ensemble is analyzed to determine whether the model is operating as expected. This helps in determining whether collecting more data would be helpful to the model by checking the ensemble variance is greater than the threshold. In this case, the MD trajectory is terminated and the final atomic configuration is placed on a queue (SQL database) for DFT calculations and added to training dataset. Periodically, ML model is retrained to the updated training model. This Active Learning loop is iterated until the cost of MD simulations becomes prohibitively expensive. The MD simulations will hopefully be sufficiently robust to support nucleation after many active learning iterations. In this sense, active learning scheme must automatically discover the important low energy and nonequilibrium physics.

Nebgen, Benjamin↗

The Joint Assimilation of Remotely Sensed Leaf Area Index and Surface Soil Moisture into a Land Surface Model

This work tests the hypothesis that jointly assimilating satellite observations of leaf area index and surface soil moisture into a land surface model improves the estimation of land vegetation and water variables. An Ensemble Kalman Filter is used to test such hypothesis across the Contiguous United States during April 2015 – December 2018. The performance of the proposed methodology is assessed for several modeled vegetation and water variables (evapotranspiration, net ecosystem exchange, and soil moisture) in terms of random errors and anomaly correlation coefficients against a set of independent validation datasets (i.e., Global Land Evaporation Amsterdam Model, FLUXCOM, and International Soil Moisture Network). Results show that the assimilation of leaf area index mostly improves the estimation of evapotranspiration and net ecosystem exchange, whereas the assimilation of surface soil moisture alone improves surface soil moisture content, especially in the western US, in terms of both root mean squared error and anomaly correlation coefficient. The joint assimilation of vegetation and soil moisture information combines the results of individual vegetation and soil moisture assimilations and reduces errors (and increases correlations with the reference datasets) in evapotranspiration, net ecosystem exchange, and surface soil moisture simulated by the land surface model. However, because soil moisture satellite observations only provide information on the water content in the top 5 cm of the soil column, the impact of the proposed data assimilation technique on root zone soil moisture is limited. This work moves one step forward in the direction of improving our estimation and understanding of land surface interactions using a multi-variate data assimilation approach, which can be particularly useful in regions of the world where ground observations are sparse or missing altogether.

Data assimilation↗

Entanglement-fidelity limits of photonically networked atomic qubits from recoil and timing

The remote entanglement of two atomic quantum memories through photonic interactions is accompanied by atomic momentum recoil. When the interactions occur at different times, such as from the random emission over the lifetime of the atomic excited state, the difference in recoil timing can expose “which-path” information and ultimately lead to decoherence. Time-bin encoded photonic qubits can be particularly sensitive to asynchronous recoil timing. In this paper we study the limits of entanglement fidelity in atomic systems due to recoil and other timing imbalances and show how these effects can be suppressed or even eliminated through proper experimental design.

Quantum communication↗

A Bayesian model for multivariate discrete data using spatial and expert information with application to inferring building attributes

When modeling sparsely observed multivariate data, strong prior information elicited from experts can be used to bolster predictive accuracy and counteract sampling bias. Similarly, modeling autocorrelation in space can help make use of co-occurrence patterns present in many types of spatial data. To make use of both expert prior information and spatial structure, we propose a novel graphical model for a spatial Bayesian network developed specifically to address challenges in inferring the attributes of buildings from geographically sparse observational data. This model is implemented as the sum of a spatial multivariate Gaussian random field and a tabular conditional probability function in real-valued space prior to projection onto the probability simplex. This modeling form is especially suitable for the usage of prior information in the form of sets of atomic rules obtained from experts. To perform inference with missing data, we implement a Markov chain Monte Carlo scheme composed of alternating steps of Gibbs sampling of missing entries and Hamiltonian Monte Carlo for model parameters. A case study in building attribution is presented to highlight the advantages and limitations of this approach.

97 MATHEMATICS AND COMPUTING↗

Non-mercury methylating microbial taxa are integral to understanding links between mercury methylation and elemental cycles in marine and freshwater sediments

The goal of this study was to explore the role of non-mercury (Hg) methylating taxa in mercury methylation and to identify potential links between elemental cycles and Hg methylation. Statistical approaches were utilized to investigate the microbial community and biochemical functions in relation to methylmercury (MeHg) concentrations in marine and freshwater sediments. Sediments were collected from the methylation zone (top 15 cm) in four Hg-contaminated sites. Both abiotic (e.g., sulfate, sulfide, iron, salinity, total organic matter, etc.) and biotic factors (e.g., hgcA, abundances of methylating and non-methylating taxa) were quantified. Random forest and stepwise regression were performed to assess whether non-methylating taxa were significantly associated with MeHg concentration. Co-occurrence and functional network analyses were constructed to explore associations between taxa by examining microbial community structure, composition, and biochemical functions across sites. Regression analysis showed that approximately 80% of the variability in sediment MeHg concentration was predicted by total mercury concentration, the abundances of Hg methylating taxa, and the abundances of the non-Hg methylating taxa. The co-occurrence networks identified Paludibacteraceae and Syntrophorhabdaceae as keystone non Hg methylating taxa in multiple sites, indicating the potential for syntrophic interactions with Hg methylators. Strong associations were also observed between methanogens and sulfate-reducing bacteria, which were likely symbiotic associations. The functional network results suggested that non-Hg methylating taxa play important roles in sulfur respiration, nitrogen respiration, and the carbon metabolism-related functions methylotrophy, methanotrophy, and chemoheterotrophy. Interestingly, keystone functions varied by site and did not involve carbon- and sulfur-related functions only. In conclusion, our findings highlight associations between methylating and non-methylating taxa and sulfur, carbon, and nitrogen cycles in sediment methylation zones, with implications for predicting and understanding the impact of climate and land/sea use changes on Hg methylation.

54 ENVIRONMENTAL SCIENCES↗

Telecommunications media for the delivery of educational programming

The technical characteristics of various telecommunications media are examined for incorporation into educational networks. FM radio, AM radio, and VHF and UHF television are considered along with computer-aided instruction. The application of iteration networks to library systems, and microform technology are discussed. The basic principles of the communications theory are outlined, and the operation of the PLATO 4 random access system is described.

Ballard, R.↗

Dynamic Flow Management Problems in Air Transportation

In 1995, over six hundred thousand licensed pilots flew nearly thirty-five million flights into over eighteen thousand U.S. airports, logging more than 519 billion passenger miles. Since demand for air travel has increased by more than 50% in the last decade while capacity has stagnated, congestion is a problem of undeniable practical significance. In this thesis, we will develop optimization techniques that reduce the impact of congestion on the national airspace. We start by determining the optimal release times for flights into the airspace and the optimal speed adjustment while airborne taking into account the capacitated airspace. This is called the Air Traffic Flow Management Problem (TFMP). We address the complexity, showing that it is NP-hard. We build an integer programming formulation that is quite strong as some of the proposed inequalities are facet defining for the convex hull of solutions. For practical problems, the solutions of the LP relaxation of the TFMP are very often integral. In essence, we reduce the problem to efficiently solving large scale linear programming problems. Thus, the computation times are reasonably small for large scale, practical problems involving thousands of flights. Next, we address the problem of determining how to reroute aircraft in the airspace system when faced with dynamically changing weather conditions. This is called the Air Traffic Flow Management Rerouting Problem (TFMRP) We present an integrated mathematical programming approach for the TFMRP, which utilizes several methodologies, in order to minimize delay costs. In order to address the high dimensionality, we present an aggregate model, in which we formulate the TFMRP as a multicommodity, integer, dynamic network flow problem with certain side constraints. Using Lagrangian relaxation, we generate aggregate flows that are decomposed into a collection of flight paths using a randomized rounding heuristic. This collection of paths is used in a packing integer programming formulation, the solution of which generates feasible and near-optimal routes for individual flights. The algorithm, termed the Lagrangian Generation Algorithm, is used to solve practical problems in the southwestern portion of United States in which the solutions are within 1% of the corresponding lower bounds.

Patterson, Sarah Stock↗

Network-Scale Ubiquitous Volume Estimation Using Tree-Based Ensemble Learning Methods

Currently ubiquitous volume data for roadway networks remains the key missing dimension in traffic operations. Most volume data are average annual daily traffic (AADT) measures derived from the Highway Performance Monitoring System (HPMS). Although methods to factor the AADT to hourly averages for typical day of week exist, actual volume data is limited to a sparse collection of locations in which volumes are continuously recorded. This paper/poster explores the use of state-of-art machine learning techniques to estimate accurate volume measures that span the highway network providing ubiquitous coverage in space, and point-in-time measures for a specific date and time. Three tree-based ensemble learning models, random forest (RF), gradient boost machine (GBM), and extreme gradient boost (XGBoost), were tested for volume estimation by learning from combined dataset of commercial probe data provided by TomTom, the FHWA's Travel Monitoring Analysis System (TMAS) data, and other infrastructure attributes such as number of lanes, speed limit, and weather. The methods were tested on major corridors and freeways in the metropolitan area of Denver. All three machine learning methods were able to provide hourly volume estimates 24 hours a day, 7 days a week, and 365 days a year with around 18% mean absolute error to true volume and about 5% of error with respect to roadway capacity. The low error measures allow the potential application by transportation agencies.

33 ADVANCED PROPULSION SYSTEMS↗