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 1,297 records · Page 72

High Resolution Radio Imaging of the Merging Galaxies NGC3256 and NGC4194

We present new 6cm and 4cm radio continuum images of the central regions of the merging galaxy systems NGC3256 and NGC4194. NGC3256 is imaged with a resolution of approx. 1 in. or approx. 190pc; NGC4194 is imaged with a resolution of approx. 0.3 in. or approx. 50pc. In both systems, we detect numerous compact radio sources embedded in more diffuse radio emission. We detect 65 compact sources in NGC3256 at 6cm and we detect 46 compact sources in NGC4194, both to a limiting luminosity of approx. 5 x 10(exp 18) W/ Hz or approx. 5 times the luminosity of Cas A. Most of the compact radio sources are loosely associated with active star forming regions but not with specific optical emission sources. Several compact radio sources in NGC3256 are near positions of compact X-ray sources detected by Lira et al.. In both NGC3256 and NGC4194, we are able to measure reliable spectral indices for the stronger sources. We find in NGC3256 approx. 20% have nominally flat radio spectral indices (indicating they are dominated by thermal radio emission from HII regions) while approx. 80% have nominally steep spectral indices (indicating they are dominated by nonthermal emission from supernova remnants). In NGC4194, half the compact radio sources have flat spectral indices and half have steep indices. For the flat-spectrum sources, we estimate the number of young massive stars and the associated ionized gas masses. For the steep-spectrum sources, we estimate supernova rates. We compare these results with those from other well-studied merging galaxy systems. We gratefully acknowledge use of the NRAO Very Large Array (VLA) and the VLA Archive. NRAO is a facility of the National Science Foundation, operated under cooperative agreement by Associated Universities, Inc.

Neff, S. G.↗

Results for the Subjective Habitability Acceptability Questionnaire From HERA C5 and C6

BACKGROUND The impact of habitability on psychological and performance outcomes is recognized as an important factor in habitat and vehicle design—especially as missions increase in length. Formal evaluation of habitability’s impact is possible using the Spaceflight Habitability Acceptability Questionnaire (SHAQ), which was psychometrically validated by Roma et al. in 2022 [1]. The SHAQ survey systematically quantifies the relationships between habitat areas (e.g., bedroom, kitchen, hygiene area, office, etc.) and behavioral health and performance (BHP) outcomes (i.e., individual and team performance, mood, stress, sleep, social interaction), and considers the moderating effects of key habitability aspects (i.e., privacy, control, convenience, efficiency, comfort, social density) on that relationship. METHOD Respondents were HERA participants from Campaign 5 and 6 (n = 32). Data from three environments were used in our SHAQ analyses: SHAQ in HERA (in-mission), SHAQ at Home (pre-mission), and SHAQ in Hotel (pre- and post-mission for C6, just pre-mission for C5). For each habitat area, respondents are asked how that area layout impacts six BHP outcomes on a 201-point visual analog scale with the anchors: impairs, has no effect, and facilitates. They also had the option to respond “Not Applicable”. Follow up questions asked why six key habitability aspects (e.g., privacy) influenced the way the habitat area affected a given outcome. Participants indicated whether a habitat area’s habitability aspects were inadequate (low score) to adequate (midpoint) to ideal (high score), using a 201-point visual analog scale as well but were binned to a 7-point scale for analyses. RESULTS AND CONCLUSION When comparing across environments, ratings of BHP outcomes were worse in HERA compared to home and pre-mission hotel, but rebounded for the post-mission hotel such that post-mission hotel rating were often higher than all previous ratings. Within HERA, ratings of how a habitat area impacted an outcome skewed towards facilitation with communal areas (i.e., galley and office) rated most facilitating. The follow up questions about why six key aspects related to habitability moderated outcomes indicated that most habitat areas were rated as having “adequate” habitability characteristics. Privacy was the only moderator with average scores below adequate for all habitat areas. The office had the most moderators scoring above adequate for all habitat areas. Ratings indicated small to moderate effects of layouts on BHP outcomes: with facilitating layouts averaging between 4 and 5, and impairing layouts averaging between 2 and 4. The sleeping module layout facilitated sleep, with ratings indicating adequate comfort but inadequate privacy. The galley layout facilitated social interaction and mood, with ratings indicating adequate efficiency and convenience. The office layout facilitated social interaction, team and individual performance, with ratings similarly indicating adequate efficiency and convenience. The hygiene module worsened individual performance, mood and stress, with ratings indicating less than adequate privacy and social density. Our findings suggest there is value in having dedicated team and individual spaces. REFERENCES [1] Roma, P. G., Landon, L.B., Spencer, C. A., Whitmire, A. M., & Williams, T. J. (2022). SHAQ: Development and Validation. Human Factors: Special Issue on Spaceflight, 65(6), 1074-1104. ACKNOWLEDGEMENTS This project was funded by the NASA Human Research Program directed task Human Factors and Behavioral Performance in HERA (PI = S. Bell) to the BHP Lab. KBR Inc. and JES Tech authors were supported by KBR’s Human Health and Performance Contract NNJ15HK11B with NASA.

behavioral health↗

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↗

Estimation of photosynthetic capacity using polarization

The authors examine the hypothesis that some part of the ecosystem-dependent variability of the vegetation indices is attributable to the effects of specular reflection of sunlight by leaves. A new class of vegetation indices, or 'minus specular' vegetation indices, is defined to account for the effects of specularly reflected light. Results show that the 'minus specular' indices, when compared to the traditional vegetation indices, potentially provide better estimates of the photosynthetic activity of a canopy than the traditional vegetation indices, particularly as a function of sun and view angles.

Rondeaux, G. M.↗

A discussion of plausible solar irradiance variations, 1700-1992

From satellite observations the solar total irradiance is known to vary. Sunspot blocking, facular emission, and network emission are three identified causes for the variations. In this paper we examine several different solar indices measured over the past century that are potential proxy measures for the Sun's irradiance. These indices are (1) the equatorial solar rotation rate, (2) the sunspot structure, the decay rate of individual sunspots, and the number of sunspots without umbrae, and (3) the length and decay rate of the sunspot cycle. Each index can be used to develop a model for the Sun's total irradiance as seen at the Earth. Three solar indices allow the irradiance to be modeled back to the mid-1700s. The indices are (1) the length of the solar cycle, (2) the normalized decay rate of the solar cycle, and (3) the mean level of solar activity. All the indices are well correlated, and one possible explanation for their nearly simultaneous variations is changes in the Sun's convective energy transport. Although changes in the Sun's convective energy transport are outside the realm of normal stellar structure theory (e.g., mixing length theory), one can imagine variations arising from even the simplest view of sunspots as vertical tubes of magnetic flux, which would serve as rigid pillas affecting the energy flow patterns by ensuring larger-scale eddies. A composite solar irradiance model, based upon these proxies, is compared to the northern hemisphere temperature depatures for 1700-1992. Approximately 71% of the decadal variance in the last century can be modeled with these solar indices, although this analysis does not include anthropogenic or other variations which would affect the results. Over the entire three centuries, approx. 50% of the variance is modeled. Both this analysis and previous similar analyses have correlations of model solar irradiances and measured Earth surface temperatures that are significant at better than the 95% confidence level. To understand our present climate variations, we must place the anthropogenic variations in the context of natural variability from solar, volcanic, oceanic, and other sources.

Hoyt, Douglas V.↗

Low-Cost Evaluation of EO-1 Hyperion and ALI for Detection and Biophysical Characterization of Forest Logging in Amazonia (NCC5-481)

Major uncertainties exist regarding the rate and intensity of logging in tropical forests worldwide: these uncertainties severely limit economic, ecological, and biogeochemical analyses of these regions. Recent sawmill surveys in the Amazon region of Brazil show that the area logged is nearly equal to total area deforested annually, but conversion of survey data to forest area, forest structural damage, and biomass estimates requires multiple assumptions about logging practices. Remote sensing could provide an independent means to monitor logging activity and to estimate the biophysical consequences of this land use. Previous studies have demonstrated that the detection of logging in Amazon forests is difficult and no studies have developed either the quantitative physical basis or remote sensing approaches needed to estimate the effects of various logging regimes on forest structure. A major reason for these limitations has been a lack of sufficient, well-calibrated optical satellite data, which in turn, has impeded the development and use of physically-based, quantitative approaches for detection and structural characterization of forest logging regimes. We propose to use data from the EO-1 Hyperion imaging spectrometer to greatly increase our ability to estimate the presence and structural attributes of selective logging in the Amazon Basin. Our approach is based on four "biogeophysical indicators" not yet derived simultaneously from any satellite sensor: 1) green canopy leaf area index; 2) degree of shadowing; 3) presence of exposed soil and; 4) non-photosynthetic vegetation material. Airborne, field and modeling studies have shown that the optical reflectance continuum (400-2500 nm) contains sufficient information to derive estimates of each of these indicators. Our ongoing studies in the eastern Amazon basin also suggest that these four indicators are sensitive to logging intensity. Satellite-based estimates of these indicators should provide a means to quantify both the presence and degree of structural disturbance caused by various logging regimes. Our quantitative assessment of Hyperion hyperspectral and ALI multi-spectral data for the detection and structural characterization of selective logging in Amazonia will benefit from data collected through an ongoing project run by the Tropical Forest Foundation, within which we have developed a study of the canopy and landscape biophysics of conventional and reduced-impact logging. We will add to our base of forest structural information in concert with an EO-1 overpass. Using a photon transport model inversion technique that accounts for non-linear mixing of the four biogeophysical indicators, we will estimate these parameters across a gradient of selective logging intensity provided by conventional and reduced impact logging sites. We will also compare our physical ly-based approach to both conventional (e.g., NDVI) and novel (e.g., SWIR-channel) vegetation indices as well as to linear mixture modeling methods. We will cross-compare these approaches using Hyperion and ALI imagers to determine the strengths and limitations of these two sensors for applications of forest biophysics. This effort will yield the first physical ly-based, quantitative analysis of the detection and intensity of selective logging in Amazonia, comparing hyperspectral and improved multi-spectral approaches as well as inverse modeling, linear mixture modeling, and vegetation index techniques.

Asner, Gregory P.↗

Geology and geochemistry of paleosols developed on the Hekpoort Basalt, Pretoria Group, South Africa

The Hekpoort paleosols comprise a regional paleoweathering horizon developed on 2.224 +/- 0.021 Ga basaltic andesite lavas at the top of the Hekpoort Formation of the Pretoria Group, Transvaal Supergroup, South Africa. In five separate profiles, from outcrops along road cuts near Waterval Onder and the Daspoort Tunnel and in three drill cores from the Bank Break Area (BB3, BB8, and BB14), the top of the paleosol is a sericite-rich zone. The sericite zone grades downward into a chlorite-rich zone. In core BB8 and in the road cut at the Daspoort Tunnel, we sampled the underlying or parent basaltic andesite into which the chlorite zone grades. We did not obtain samples of the parent material at Waterval Onder and in cores BB3 and BB14, but chemical analyses indicate that the chlorite and sericite zones in these profiles derive from underlying lavas similar to the ones we sampled in core BB8 and at the Daspoort Tunnel. The presence of apparent rip-up clasts of the paleosol in the overlying ironstones of the Strubenkop Formation in the cores from Bank Break makes it very unlikely that most of the alteration was a result of interactions with hydrothermal fluids. Desiccation cracks at the top of the paleosol that were filled with sand during the deposition of the overlying sediments at Waterval Onder point to a subaerial weathering origin. Very little, if any, Al, Ti, Zr, V, or Cr moved a discernible distance during weathering of any of the five profiles. The vertical distribution of Fe, Mg, Mn, Ni, and Co indicates that these elements were largely removed from the top of the soil during weathering. The overall abundance of these elements in each of the profiles indicates that a significant fraction of the complement lost from the top subsequently reprecipitated in the lower portion of the soil as constituents of an Fe2(+) -rich smectite. The loss of Fe from the top of the soil during weathering of the Hekpoort paleosols indicates that atmospheric PO2 was less than 8 x 10(-4) atm about 2.22 Ga. Fe2(+) -rich smectite should only precipitate during soil formation if atmospheric PCO2 is less than or equal to 2 x 10(-2) atm (Rye, Kuo, and Holland, 1995). Ca and Na were largely lost during weathering. Some Na was apparently added to the sericite zone in cores BB3, BB8, and BB14 after weathering. All five profiles are enriched in K and Rb, and most are enriched in Ba. The distribution of these elements indicates that they all were added during post-weathering hydrothermal metasomatism. Rb-Sr analysis of the paleosol at the Daspoort Tunnel indicates that metasomatism last affected that profile 1.925 +/- 0.032 Ga (Macfarlane and Holland, 1991).

Non-NASA Center↗

Continued Research into Characterizing the Preturbulence Environment for Sensor Development, New Hazard Algorithms and Experimental Flight Planning

The purpose of the research was to develop and test improved hazard algorithms that could result in the development of sensors that are better able to anticipate potentially severe atmospheric turbulence, which affects aircraft safety. The research focused on employing numerical simulation models to develop improved algorithms for the prediction of aviation turbulence. This involved producing both research simulations and real-time simulations of environments predisposed to moderate and severe aviation turbulence. The research resulted in the following fundamental advancements toward the aforementioned goal: 1) very high resolution simulations of turbulent environments indicated how predictive hazard indices could be improved resulting in a candidate hazard index that indicated the potential for improvement over existing operational indices, 2) a real-time turbulence hazard numerical modeling system was improved by correcting deficiencies in its simulation of moist convection and 3) the same real-time predictive system was tested by running the code twice daily and the hazard prediction indices updated and improved. Additionally, a simple validation study was undertaken to determine how well a real time hazard predictive index performed when compared to commercial pilot observations of aviation turbulence. Simple statistical analyses were performed in this validation study indicating potential skill in employing the hazard prediction index to predict regions of varying intensities of aviation turbulence. Data sets from a research numerical model where provided to NASA for use in a large eddy simulation numerical model. A NASA contractor report and several refereed journal articles where prepared and submitted for publication during the course of this research.

Kaplan, Michael L.↗