Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Indices”

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 325 records · Page 18

Nickel-rich, volatile depleted iron meteorites: Relationships and formation processes

Ungrouped iron meteorites Tishomingo, Willow Grove, and Chinga, and group IVB iron meteorites, are Ni-rich. Similarities include enrichments of 10–100 × CI for some refractory siderophile elements, and equivalent depletions in more volatile siderophile elements. Superimposed on the overall enrichment/depletion trend, certain siderophile elements (P, W, Fe, Mo) are depleted relative to elements of similar volatility. All three ungrouped irons derive from parent bodies formed in the early Solar System. Willow Grove and Chinga are characterized by cosmic ray exposure corrected 182 W/ 184 W consistent with metal-silicate segregation on their parent bodies within 1–3 Myr of Solar System formation, within the age range determined for segregation of magmatic iron meteorite parent bodies, including group IVB irons. Tishomingo is characterized by a younger model age 4–5 Myr subsequent to Solar System formation, reflecting either late stage melting resulting from 26 Al decay, or an impact resetting. The discovery of stishovite in Tishomingo, indicating exposure to a minimum shock pressure of 8–9 GPa, is consistent with the latter. Stishovite in Tishomingo and chromite included in troilite-daubréelite in IVB irons allows oxygen isotopic composition comparison between these meteorites. Different mass independent oxygen isotopic compositions of IVB irons and Tishomingo indicate genetically distinct parent bodies. By contrast, mass independent Mo isotopic compositions overlap within analytical uncertainties, indicating a similar, carbonaceous chondrite (CC) type genetic heritage. Molybdenum and 183 W isotopic data for Chinga and Willow Grove indicate derivation from CC type parent bodies. Willow Grove shares Mo and 183 W isotopic compositions with the Ni-rich South Byron Trio (SBT) grouplet and the Milton pallasite. These Ni-rich meteorites likely formed in the same general nebular environment as other CC planetesimals, likely the outer Solar System. Highly siderophile element (HSE) abundances of Willow Grove and Tishomingo are similar to some IVB meteorites, consistent with formation by moderate degrees of fractional crystallization from initial metallic melts with low S and P, and modestly fractionated HSE. The comparable HSE abundances of Tishomingo and Willow Grove to some IVB irons, yet substantially higher Ni concentrations, indicate formation on parent bodies with lower bulk HSE abundances or HSE concentration in proportionally smaller volumes of metal. HSE abundances in Chinga are considerably lower than in IVB irons, highly fractionated, and processes responsible for these remain elusive. For IVB irons and these ungrouped irons, high temperature condensation likely dominated the enrichment and depletion of the refractory and volatile siderophile elements, respectively. Parent body degassing may have also played a role. Relative depletion of volatile siderophile elements is not, however, a universal feature of high-Ni meteorites. The SBT and Milton pallasite are Ni-rich, but less depleted in the more volatile siderophile elements. Nickel enrichment was likely driven by oxidation of Fe metal during parent body accretion or core segregation. Oxidation of the Tishomingo and Willow Grove parent bodies may have occurred at ~ IW + 1, indicated by relative Mo and W depletions due to metal/water reaction during differentiation. Late-stage reduction, indicated by the presence of Cr-bearing sulfides in Tishomingo and IVB irons, may have resulted from exhaustion of the oxidant.

36 MATERIALS SCIENCE↗

Governance and resilience as entry points for transforming food systems in the countdown to 2030

Due to complex interactions, changes in any one area of food systems are likely to impact—and possibly depend on—changes in other areas. Here we present the first annual monitoring update of the indicator framework proposed by the Food Systems Countdown Initiative, with new qualitative analysis elucidating interactions across indicators. Since 2000, we find that 20 of 42 indicators with time series have been trending in a desirable direction, indicating modest positive change. Qualitative expert elicitation assessed governance and resilience indicators to be most connected to other indicators across themes, highlighting entry points for action—particularly governance action. Literature review and country case studies add context to the assessed interactions across diets, environment, livelihoods, governance and resilience indicators, helping different actors understand and navigate food systems towards desirable change.

Schneider, Kate R. [Johns Hopkins Univ., Washingto↗

A Climatology of Clear-Air Turbulence and Mountain Wave Turbulence Throughout the Stratosphere Using MERRA-2 Data

Forty years of MERRA-2 reanalysis data is used to calculate twenty-nine unique indices for both clear-air turbulence (CAT) and mountain wave turbulence (MWT) throughout the stratosphere. The indices are combined into a single estimate of turbulence intensity for both CAT and MWT for five latitudinal domains. These indices have previously been applied at commercial flight levels, but this work is the first time such indices have been evaluated for estimating turbulence throughout the entire stratosphere. The results indicate that turbulence throughout the stratosphere is generally around four-times more intense in the winter season and increases in intensity with altitude. The data suggests the strongest source of CAT in each hemisphere is the polar vortex in the winter stratosphere above 50 hPa. The strongest sources of MWT are the Himalayas and the Andes in the lower stratosphere. The relationship between CAT and the quasi-biennial oscillation is also examined and indicates that easterly tropical winds are associated with the most intense turbulence. The trends in the climatology of CAT over forty years were compared indicating a significant increase in turbulence intensity in the upper stratosphere, and a weaker but still significant decrease in turbulence intensity in the lower stratosphere.

54 ENVIRONMENTAL SCIENCES↗

EPCAPE-PT-LANL Measurements: Single Particle Soot Photometer

Coastal cities offer a unique environment for studying aerosol-cloud interactions and the effects of urban emissions on cloud properties. As part of the Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE), the Partitioning Thrust by Los Alamos National Laboratory (EPCAPE-PT-LANL) was conducted. Our campaign focused on measuring the optical and chemical properties of aerosols and their interactions within marine stratocumulus clouds in La Jolla, California. EPCAPE-PT-LANL enhances the primary goals of EPCAPE through innovative observations of vapor-phase transitions between aerosols and cloud droplets, the impact of black carbon on aerosol-cloud dynamics, and the effects of cloud processing on aerosol optical properties. Instrument: Single Particle Soot Photometer (Droplet Measurements Technology) Data Notes: Contact us if you want additional data products from this instrument. Reported data is the black carbon (rBC) number and mass concentration with diagnostic flags. The SP2 measures incandescence from particles that is induced with a 1064 Nd-YAG laser. Particles that absorb the laser energy and then emit radiation is assumed to contain black carbon. The amount of intensity of radiation is related to the mass of absorbing material in the particle. Here, we calibrate the incandescence intensity to size selected Regal Black (Cabot) nebulized from solution. Data provided for the EPCAPE campaign used the combined broadband high-gain + low-gain channels (BHBL). Each channel has a lower threshold of detection equivalent to 2-s of the respective channel noise. Thresholding has been applied to select between the high-gain and low-gain channels. The detection limit is 80-540 nm (Deq) or 0.36- 55 fg. It is assumed that rBC is the only aerosol type that is in significant concentration in the sampled atmosphere that absorbs the laser energy. Particle data are integrated over 10 second windows to calculate a rBC number and mass concentration. QC/QA: • Diagnostic flags that impacted measured concentration: - Sample, Sheath, and Purge Flow Rates: Despite observing fluctuations in all flows, the rBC detection and mass quantification are generally observed to be stable, although large changes in sample flow rate did impact detection efficiency. A flag was implemented for sample flow deviations >12 cm3/min from the set point averaged over 30 seconds - Laser Power: Detection efficiency decreases with laser power and deviations in laser power also affect mass quantification. Flag for laser power is set for deviations in laser current from the set value >5 mA. - Primary Detector Threshold: Primary thresholding is the signal value which determines whether a “particle” is recorded. Thresholding errors occur when the threshold value is too HIGH and real particles are ignored. • Several periods without data: - 16-18 Nov: Ultra Zero Air generator failure, flows deviated significantly from set points. rBC # conc recovered but questionable. - 20-21 Nov: Offline for calibrations for several hours each day. - 25-27 Nov: Data is missing. - 30 Nov: Power outage 3 Dec: SP2 hard-drive full. - 3 Dec: Power Outage Header: - BHBL_NumbConc[#/cc]: Refers to the number concentration of refractory black carbon (rBC) measured from combined broadband high-gain and low-gain channels, expressed in particles per cubic centimeter. - BHBL_BCMass_Conc[ug/m3]: Refers to the mass concentration of refractory black carbon (rBC) measured from combined broadband high-gain and low-gain channels, expressed in micrograms per cubic meter. - NoData_Flag[bool]: A boolean flag that indicates whether no data was recorded during a measurement. - NoBC_Flag[bool]: A boolean flag indicating whether no black carbon particles were recorded during the measurement. - Laser_Flag[bool]: A boolean flag indicating deviations in laser current during the measurement. - SampleFlow_Flag[bool]: A boolean flag indicating deviations from the set sample flow rate during the measurement. - PrimThresh_Flag[bool]: A boolean flag indicating deviations from the set primary threshold during the measurement, potentially ignoring real particles. - Manual_Flag[bool]: A boolean flag indicating manual intervention or adjustments during the measurement. - CVI_Flag[bool]: A boolean flag indicating whether the Counterflow Virtual Impactor (CVI) was active or inactive during the measurement.

54 ENVIRONMENTAL SCIENCES↗

Characterization of ferredoxins involved in electron transfer pathways for nitrogen fixation implicates differences in electronic structure in tuning 2[4Fe 4S] Fd activity

Ferredoxins (Fds) are small proteins which shuttle electrons to pathways like biological nitrogen fixation. Physical properties tune the reactivity of Fds with different pathways, but knowledge on how these properties can be manipulated to engineer new electron transfer pathways is lacking. Recently, we showed that an evolved strain of Rhodopseudomonas palustris uses a new electron transfer pathway for nitrogen fixation. This pathway involves a variant of the primary Fd of nitrogen fixation in R. palustris, Fer1, in which threonine at position 11 is substituted for isoleucine (Fer1 T11I ). To understand why this substitution in Fer1 enables more efficient electron transfer, we used in vivo and in vitro methods to characterize Fer1 and Fer1 T11I . Electrochemical characterization revealed both Fer1 and Fer1 T11I have similar redox transitions (–480 mV and – 550 mV), indicating the reduction potential was unaffected despite the proximity of T11 to an iron-sulfur (Fe—S) cluster of Fer1. Additionally, disruption of hydrogen bonding around an Fe—S cluster in Fer1 by substituting threonine with alanine (T11A) or valine (T11V) did not increase nitrogenase activity, indicating that disruption of hydrogen bonding does not explain the difference in activity observed for Fer1 T11I . Electron paramagnetic resonance spectroscopy studies revealed key differences in the electronic structure of Fer1 and Fer1 T11I , which indicate changes to the high spin states and/or spin-spin coupling between the Fe—S clusters of Fer1. Finally, our data implicates these electronic structure differences in facilitating electron flow and sets a foundation for further investigations to understand the connection between these properties and intermolecular electron transfer.

59 BASIC BIOLOGICAL SCIENCES↗

A modeling pipeline to relate municipal wastewater surveillance and regional public health data

As COVID-19 becomes endemic, public health departments require improved passive indicators, which are independent of voluntary testing data, to estimate the prevalence of COVID-19 in local communities. Quantification of SARS-CoV-2 RNA from wastewater has the potential to be a powerful passive indicator. However, connecting measured SARS-CoV-2 RNA to community prevalence is challenging. We have developed a generalized methodology to improve the predictive power of wastewater measurements and applied it to data collected from treatment plants in the Chicago area. We built and compared a set of multi-linear regression models, which incorporate pepper mild mottle virus (PMMoV) as a population biomarker, Bovine coronavirus (BCoV) as a recovery control, and wastewater system flow rate into a corrected estimate for SARS-CoV-2 RNA concentration. For our data, models with BCoV performed better than those with PMMoV, but all model estimates of prevalence significantly improved correlation compared to doing no correction. We also investigate the utility of RNA measurements in wastewater as a leading indicator of COVID-19 trends. We do this in a rolling manner for corrected wastewater data and for other prevalence indicators, and statistically compare the temporal relationship between new increases in the wastewater data and those in other prevalence indicators. We find that wastewater trends often lead other COVID-19 indicators in predicting new surges.

60 APPLIED LIFE SCIENCES↗

Imprecise global sensitivity analysis using bayesian multimodel inference and importance sampling

Global Sensitivity Analysis (GSA) aims to understand the relative importance of uncertain input variables to model response. Conventional GSA involves calculating sensitivity (Sobol’) indices for a model with known model parameter distributions. However, model parameters are affected by aleatory and epistemic uncertainty, with the latter often caused by lack of data. In this paper, we propose a new framework to quantify uncertainty in probability model-form and model parameters resulting from small datasets and integrate these uncertainties into Sobol’ index estimates. First, the process establishes, through Bayesian multimodel inference, a set of candidate probability models and their associated probabilities. Imprecise Sobol’ indices are calculated from these probability models using an importance sampling reweighting approach. This results in probabilistic Sobol’ indices, whose distribution characterizes uncertainty in the sensitivity resulting from small dataset size. The imprecise Sobol’ indices thus provide a measure of confidence in the sensitivity estimate and, moreover, can be used to inform data collection efforts targeted to minimize the impact of uncertainties. Through an example studying the parameters of a Timoshenko beam, we show that these probabilistic Sobol’ indices converge to the true/deterministic Sobol’ indices as the dataset size increases and hence, distribution-form uncertainty reduces. The approach is then applied to assess the sensitivity of the out-of-plane properties of an E-glass fiber composite material to its constituent properties. This second example illustrates the approach for an important class of materials with wide-ranging applications when data may be lacking for some input parameters.

42 ENGINEERING↗

Towards cost-competitive middle distillate fuels from ethanol within a market-flexible C2 platform-based biorefinery concept

Ethanol to middle distillates (ETMD) is a promising pathway to produce sustainable liquid fuels to decarbonize the hard-to-electrify transportation sectors due to (1) the abundant sugar/starch and lignocellulosic biomass, (2) the existing deployment scale of fuel ethanol production (similar to 29 billion gallons per year globally), and (3) emerging opportunities in C2+ alcohol synthesis from CO2. Here we report a conceptual market-responsive biorefinery centered around a new ETMD pathway based on one-step ethanol to butene-rich olefins (ETO) over a Cu-Zn-Y/Beta catalyst. Specifically, this ethanol conversion pathway comprises one-step ETO, oligomerization, and hydrotreating. This ETO is distinct from that in the conventional ethanol-to-jet process which is based on two-step ethanol to ethylene and ethylene oligomerization to butenes. Butene-rich olefins can be shifted to butadiene-rich products by simply changing the reaction atmosphere from hydrogen to inert gas over the same ETO catalyst. Leveraging the experimental results, baseline techno-economic analysis (TEA) and sensitivity analysis indicate that the ethanol conversion cost is $0.60 per gallon gasoline equivalent (GGE), with opportunities for further cost reduction via improving the liquid hydrocarbon yield and space velocities, and process optimization on balancing dewatering of ethanol feed prior to the ETO step. The minimum fuel selling price (MFSP) of liquid hydrocarbons derived from corn starch ethanol with butadiene as coproduct is $1.64 per GGE, in the range that is cost competitive with petroleum kerosene-type jet fuel. Projected MFSP for cellulosic ethanol (corn stover) derived hydrocarbons is below $3.00 per GGE and co-production of butadiene further reduces the MFSP to $1.70 per GGE. The Well-to-Wake life-cycle analysis indicates that 85% greenhouse gas emission reduction can be achieved when using corn stover compared to petroleum reference and the associated carbon credits will provide significant economic incentives to favor the cellulosic ethanol-derived hydrocarbon fuels. This study demonstrates a low-cost pathway to middle distillate fuels leveraging existing ethanol infrastructure, where catalysis innovation drives the reduction of process complexity and flexible coproduction of a value-added chemical product.

Zhang, Junyan↗

Conservation of endo -glucanase 16 (EG16) activity across highly divergent plant lineages

Plant cell walls are highly dynamic structures that are composed predominately of polysaccharides. As such, endogenous carbohydrate active enzymes (CAZymes) are central to the synthesis and subsequent modification of plant cells during morphogenesis. The endo-glucanase 16 (EG16) members constitute a distinct group of plant CAZymes, angiosperm orthologs of which were recently shown to have dual β-glucan/xyloglucan hydrolase activity. Molecular phylogeny indicates that EG16 members comprise a sister clade with a deep evolutionary relationship to the widely studied apoplastic xyloglucan endo-transglycosylases/hydrolases (XTH). A cross-genome survey indicated that EG16 members occur as a single ortholog across species and are widespread in early diverging plants, including the non-vascular bryophytes, for which functional data were previously lacking. Remarkably, enzymological characterization of an EG16 ortholog from the model moss Physcomitrella patens (PpEG16) revealed that EG16 activity and sequence/structure are highly conserved across 500 million years of plant evolution, vis-à-vis orthologs from grapevine and poplar. Ex vivo biomechanical assays demonstrated that the application of EG16 gene products caused abrupt breakage of etiolated hypocotyls rather than slow extension, thereby indicating a mode-of-action distinct from endogenous expansins and microbial endo-glucanases. As a result, the biochemical data presented here will inform future genomic, genetic, and physiological studies of EG16 enzymes.

59 BASIC BIOLOGICAL SCIENCES↗

Direction of impact for explainable risk assessment modeling

Abstract Several graphical indicators have been recently introduced to help analysts visualize the marginal effects of inputs in complex models. The insights derived from such tools may help decision‐makers and risk analysts in designing interventions. However, we know little about the adequacy and consistency of different indicators. This work investigates popular marginal effect indicators to understand whether they yield indications consistent with the properties of the quantitative model under inspection. Specifically, we examine the notions of monotonicity, Lipschitz, and concavity consistency. Surprisingly, only PD functions satisfy all these notions of consistency. However, when selecting the indicators, in addition to consistency, analysts need to consider the risk of model extrapolation. For situations where such risk is under control, we utilize individual conditional expectations together with PD plots. Two applications, on a NASA space risk assessment model and a susceptible exposed infected recovered (SEIR) model for the COVID‐19 pandemic illustrate the insights obtained from these indicators.

Mathematical Methods In Social Sciences↗

First Phase Consensus Roadmap for Development of Condition-Based Cable Reliability Assurance

The objective of this work was to develop a first phase consensus roadmap for condition-based qualification (CBQ) of electrical cables. With CBQ, qualification of Class 1E electrical cables moves from a time-based approach to a condition-based approach, which is anticipated to be safer in terms of reliability and conservatism, and more cost effective in the long run. However, due to barriers, the CBQ approach has not yet been adopted by U.S. nuclear power plants (NPPs). Based upon a review of current work evaluating CBQ, the limitation of available condition monitoring technology seems to be the largest barrier. The importance of condition monitoring, or more specifically selecting appropriate condition indicators, during CBQ cannot be understated. However, selecting appropriate condition indicators is challenged by techniques that are destructive and only evaluate cable degradation locally. Further, arguably, no one identified condition indicator fully establishes cable condition. Thus, additional work is necessary to evaluate potential condition indicators towards CBQ. In addition to the requirements of IEC/IEEE Std. 60780-323, ideal condition indicators should include a) both destructive and non-destructive approaches, b) both local and global measurements, c) real-time (i.e., online) monitoring that trends with degradation, d) enable correlation with qualified levels of degradation, and e) be established within a repository of condition indicators with applicable materials and/or components and their acceptance criteria. Additional work is needed in development of technology and methodology prior to adoption of CBQ, especially for extending qualified life of installed components. Education and early experience by the industry and regulators will be required for this change in approach as an alternative to re-analysis. A series of workshops that bring together stakeholders to identify and address gaps will be needed. The longstanding cooperative working group of cable researchers from the U.S. Department of Energy, the Electric Power Research Institute, and the Nuclear Regulatory Commission forms a valuable starting point for development of a consensus roadmap to condition-based qualification approach as a viable options for qualification of cable systems in U.S. light water reactors.

42 ENGINEERING↗

A double peak in the seasonality of California's photosynthesis as observed from space

Solar-induced chlorophyll fluorescence (SIF) has been shown to be a powerful proxy for photosynthesis and gross primary productivity (GPP). The recently launched TROPOspheric Monitoring Instrument (TROPOMI) features the required spectral resolution and signal-to-noise ratio to retrieve SIF from space. Here, we present a downscaling method to obtain 500 m spatial resolution SIF over California. We report daily values based on a 14 d window. TROPOMI SIF data show a strong correspondence with daily GPP estimates at AmeriFlux sites across multiple ecosystems in California. We find a linear relationship between SIF and GPP that is largely invariant across ecosystems with an intercept that is not significantly different from zero. Measurements of SIF from TROPOMI agree with MODerate Resolution Imaging Spectroradiometer (MODIS) vegetation indices – the normalized difference vegetation index (NDVI), enhanced vegetation index (EVI), and near-infrared reflectance of vegetation index (NIR$_v$) – at annual timescales but indicate different temporal dynamics at monthly and daily timescales. TROPOMI SIF data show a double peak in the seasonality of photosynthesis, a feature that is not present in the MODIS vegetation indices. The different seasonality in the vegetation indices may be due to a clear-sky bias in the vegetation indices, whereas previous work has shown SIF to have a low sensitivity to clouds and to detect the downregulation of photosynthesis even when plants appear green. We further decompose the spatiotemporal patterns in the SIF data based on land cover. The double peak in the seasonality of California's photosynthesis is due to two processes that are out of phase: grasses, chaparral, and oak savanna ecosystems show an April maximum, while evergreen forests peak in June. An empirical orthogonal function (EOF) analysis corroborates the phase offset and spatial patterns driving the double peak. The EOF analysis further indicates that two spatiotemporal patterns explain 84 % of the variability in the SIF data. Results shown here are promising for obtaining global GPP at sub-kilometer spatial scales and identifying the processes driving carbon uptake.

54 ENVIRONMENTAL SCIENCES↗

FUN-BioCROP model with litter decomposition parameters derived from the LIDET dataset

This repository contains the code and data necessary to run the FUN-BioCROP (Fixation and Uptake of Nitrogen-Bioenergy Carbon, Rhizosphere, Organisms, and Protection) model with litter decomposition parameters derived from a modified Monte Carlo simulation that used the Long-term Intersite Decomposition Experiment Team dataset. Related publication:Juice, S.M., Ridgeway, J.R., Hartman, M.D., Parton, W.J., Berardi, D.M., Sulman, B.N., Allen, K.E., & Brzostek, E.R. Reparameterizing litter decomposition using a simplified Monte Carlo method improves litter decay simulated by a microbial model and alters bioenergy soil carbon estimates. Description of Files: FUNBioCROP_LIDET Study.Rmd R code with FUN-BioCROP model that can be run with 10 different sets of parameters for litter decomposition (Baseline, LIDET, or eight other best parameter sets identified in the modified Monte Carlo simulation. CORPSE Functions_Bioenergy_V2.R Code with CORPSE model functions, called by FUNBioCROP_LIDET Study.Rmd Model Input Data: bulk.csv, bulk_till.csv, rhizo.csv, rhizo_till.csv, litter.csv Initial C and N (kg C or N/m2) pool values for each soil compartment, final values from spin up. All five files have the same columns: (Column - Description - Units) uFastC - Unprotected fast decomposing carbon - kg carbon/m2 uSlowC - Unprotected slow decomposing carbon - kg carbon/m2 uNecroC - Unprotected necromass carbon - kg carbon/m2 pFastC - Protected fast decomposing carbon - kg carbon/m2 pSlowC - Protected slow decomposing carbon - kg carbon/m2 pNecroC - Protected necromass carbon - kg carbon/m2 livingMicrobeC - Carbon in living microbial biomass - kg carbon/m2 uFastN - Unprotected fast decomposing nitrogen - kg nitrogen/m2 uSlowN - Unprotected slow decomposing nitrogen - kg nitrogen/m2 uNecroN - Unprotected necromass nitrogen - kg nitrogen/m2 pFastN - Protected fast decomposing nitrogen - kg nitrogen/m2 pSlowN - Protected slow decomposing nitrogen - kg nitrogen/m2 pNecroN - Protected necromass nitrogen - kg nitrogen/m2 inorganicN - Inorganic nitrogen - kg nitrogen/m2 CO2 - Carbon in carbon dioxide - kg carbon/m2 livingMicrobeN - Nitrogen in living microbial biomass - kg nitrogen/m2 Model Input Data: FluxTower_AvgSoilT.csv: Average daily soil temperature (oC) at 10 cm depth at University of Illinois Urbana-Champaign (UIUC) Energy Farm flux tower from 7/2008-3/2016. (One year of averaged data) Model Input Data: FluxTower_AvgSoilVWC.csv: Average daily soil volumetric water content (VWC) at 10 cm depth at UIUC Energy Farm flux tower from 7/2008-3/2016. (One year of averaged data) Model Input Data: input_CCS_LIDET Study.csv: This file has daily data to run FUN-BioCROP (Column - Description - Units): yr - calendar year - year doy - day of year (1 to 365) (no leap year) - day anpp - aboveground NPP (DayCent) - kg C/m2/day bnpp - belowground NPP (DayCent) - kg C/m2/day aglivc - live aboveground biomass carbon (DayCent) - kg C/m2 bglivcj - live juvenile fine root biomass carbon (DayCent) - kg C/m2 bglivcm - live mature fine root biomass carbon (DayCent) - kg C/m2 aglivn - live aboveground biomass nitrogen (DayCent) - kg N/m2 bglivnj - live juvenile fine root biomass nitrogen (DayCent) - kg N/m2 bglivnm - live mature fine root biomass nitrogen (DayCent) - kg N/m2 nyr - simulation year - year cult - indicates a cultivation event (0 or 1) crop - indicates a new crop (0 or 1) fert - indicates a fertilizer event (0 or 1) frst - indicates the first day of the growing season (0 or 1) harv - indicates a harvest event (0 or 1) last - indicates the end of the growing season (0 or 1) croptype - crop type (0=none; 1=alfalfa; 2=corn; 3=grass clover pasture; 4=soybean; 5=wheat) cropsrl - crop specific root length - mm/g root cultrhizmix - fraction of rhizosphere mixed with bulk soil during cultivation (0.0-1.0) - fraction cultlitmix - fraction of litter mixed with bulk soil during cultivation (0.0-1.0) - fraction harvremov - fraction of above ground biomass removed during harvest (0.0-1.0) - fraction fertamt - fertilization amount - g N/m2 lifehist - plant life history (0 = annual, 1 = perennial) froot_turnover_c - amount of C in fine root turnover - kg C/m2 froot_turnover_n - amount of N in fine root turnover - kg N/m2 agrd_turnover_c - amount of C in aboveground biomass turnover - kg C/m2 agrd_turnover_n - amount of N in aboveground biomass turnover - kg N/m2 leaf_litter_fastfrac - Fast decomposing fraction of leaf litter (0.0-1.0) - fraction root_litter_fastfrac - Fast decomposing fraction of root litter (0.0-1.0) - fraction root_diameter - root diameter - mm root_length - root length - mm root/m2 rhizo_frac - fraction of total soil volume that is rhizosphere (0.0 - 1.0) - fraction date - date in format YYYY-MM-DD Instructions: Save the model code ("FUN-BioCROP_LIDET Study.Rmd") and accompanything files (data streams and CORPSE function code) in the same folder. In model code "Chunk 3: Load CORPSE Data Streams" set the working directory (setwd) to the folder with the files saved in step #1. In "Chunk 5: Define LIDET parameter sets" select the litter decomposition parameter set to be used in the run, and comment out all other sets. If changing any parameter values, edit them in "Chunk 6: Load parameters." Run all chunks up to and including "Chunk 10: Prepare Data for Export." In "Chunk 11: Export Output Data" edit data frames for export and filenames, as necessary. "Chunk 12: Graph Total Soil C" makes a figure of C remaining over the model run period. Description of each model chunk (in file FUN-BioCROP_LIDET Study.Rmd): Chunk 1: Remove all functions, clear memory. Removes all functions from R environment, clears the memory. Chunk 2: Load Packages. Loads packages necessary to run the code. Chunk 3: Load CORPSE Data Streams. Sets the working directory and loads the data files necessary to run CORPSE. Chunk 4: Load CORPSE Functions. Loads the R script with CORPSE functions from the working directory, "CORPSE Functions_Bioenergy_V2.R". Chunk 5: Define LIDET parameter sets. Has ten different parameter sets for litter decomposition tested in this study: Baseline parameters, LIDET parameters, and the other 8 best performing parameter sets identified in the modified Monte Carlo. To run the model, all but one parameter set must be commented out. Chunk 6: Load Parameters. Loads all fixed parameters to run the model. Data frame with definitions of parameters is in the CORPSE function script "CORPSE Functions_Bioenergy_V2.R" Chunk 7: Prepare Data Streams. Takes data streams loaded in Chunk 3 and puts them in the format necessary to run the model. The model is coded to run at least two sites at a time, so if only one site is being run it must be run in duplicate. Individual data tables of daily values are created in this chunk from the input data file. Chunk 8: Set Initial Conditions. Creates data tables of soil C and N pools for each soil compartment (rhizo_till, rhizo, bulk_till, bulk, litter) and loads initial values into the data tables. Creates lists for each soil compartment to hold model output. Chunk 9: Load FUN Data and Set Up Matrices. Uses DayCent data to calculate FUN input data: root and leaf N demand, total N demand, plant CN, leaf N available for retranslocation, and litter production. Creates matrices for FUN model outputs. Chunk 10: Run Model. Runs the model. Chunk 11: Prepare Data for Export. Combines data from each day saved as lists into data frames for each soil compartment. Adds values from all soil compartments together to calculate total soil values, creates separate data frames for each soil C and N pool (e.g., protected slow C) for the total soil value. Adds different C and N pools together to calculate total soil C and N for all layers. Creates data frame of ratio of protected to unprotected SOC. Organizes FUN data for export. Chunk 12: Export Results. Exports CSV files of model results to the working directory. Chunk 13: Graph Total Soil C. Makes figure of C remaining over time. Related Links: Original FUN-BioCROP model: https://github.com/BrzostekEcologyLab/FUN-BioCROP LIDET dataset: https://andlter.forestry.oregonstate.edu/data/abstract.aspx?dbcode=TD023

Juice, Stephanie↗

Bias control for a memory device

Methods, systems, and devices for bias control for a memory device are described. A memory system may store indication of whether data is coherent. In some examples, the indication may be stored as metadata, where a first value indicates that the data is not coherent and a second value or a third value indicate that the data is coherent. When a processing unit or other component of the memory system processes a command to access data, the memory system may operate according to a device bias mode when the indication is the first value, and according to a host bias mode when the indication is the second value or the third value.

97 MATHEMATICS AND COMPUTING↗

Bias control for a memory device

Methods, systems, and devices for bias control for a memory device are described. A memory system may store indication of whether data is coherent. In some examples, the indication may be stored as metadata, where a first value indicates that the data is not coherent and a second value or a third value indicate that the data is coherent. When a processing unit or other component of the memory system processes a command to access data, the memory system may operate according to a device bias mode when the indication is the first value, and according to a host bias mode when the indication is the second value or the third value.

Walker, Dean↗

Verification of the sensitivity and uncertainty-based criticality safety validation techniques: ORNL’s SCALE case study

The use of sensitivity and uncertainty-based techniques has been well-developed in the neutronic community for criticality safety applications. Sensitivity techniques calculate the first-order relative variations in keff due to relative cross-section variations by isotope, reaction type, and incident neutron energy. The result is a vector of sensitivity coefficients which is folded with prior cross-section uncertainties to calculate the so-called similarity indices. Similarity indices are integral quantities which serve as useful mathematical measures for assessing the relevance of existing or new experiments to given application conditions, representing design systems for which the models are to be validated. Despite their importance for a wide range of validation activities, the peer-reviewed literature is non-existent on the development of rigorous verification analysis for their calculated values. Thus, this manuscript proposes a simple non-intrusive method to calculate similarity indices, serving as an independent verification tool for their evaluation. The manuscript derives the theoretical basis for the proposed verification analysis. Next, the verification analysis is applied to the ORNL’s SCALE code suite which employs adjoint-based methods to calculate similarity indices. Numerical experiments using a number of well-known benchmark experiments are employed to verify the calculation of the similarity indices.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Predictability and empirical dynamics of fisheries time series in the North Pacific

Previous studies have documented a strong relationship between marine ecosystems and large-scale modes of sea surface height (SSH) and sea surface temperature (SST) variability in the North Pacific such as the Pacific Decadal Oscillation and the North Pacific Gyre Oscillation. In the central and western North Pacific along the Kuroshio-Oyashio Extension (KOE), the expression of these modes in SSH and SST is linked to the propagation of long oceanic Rossby waves, which extend the predictability of the climate system to ~3 years. Using a multivariate physical-biological linear inverse model (LIM) we explore the extent to which this physical predictability leads to multi-year prediction of dominant fishery indicators inferred from three datasets (i.e., estimated biomasses, landings, and catches). We find that despite the strong autocorrelation in the fish indicators, the LIM adds dynamical forecast skill beyond persistence up to 5-6 years. By performing a sensitivity analysis of the LIM forecast model, we find that two main factors are essential for extending the dynamical predictability of the fishery indicators beyond persistence. The first is the interaction of the fishery indicators with the SST/SSH of the North and tropical Pacific. The second is the empirical relationship among the fisheries time series. This latter component reflects stock-stock interactions as well as common technological and human socioeconomic factors that may influence multiple fisheries and are captured in the training of the LIM. These results suggest that empirical dynamical models and machine learning algorithms, such as the LIM, provide an alternative and promising approach for forecasting key ecological indicators beyond the skill of persistence.

60 APPLIED LIFE SCIENCES↗