Search NASA⌕ Search

SEARCH · Search NASA

Results for “variable threshold”

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

Mineral Dust Cycle in the Multiscale Online Nonhydrostatic AtmospheRe CHemistry Model (Monarch) Version 2.0

We present the dust module in the Multiscale Online Non-hydrostatic AtmospheRe CHemistry model (MONARCH) version 2.0, a chemical weather prediction system that can be used for regional and global modeling at a range of resolutions. The representations of dust processes in MONARCH were upgraded with a focus on dust emission (emission parameterizations, entrainment thresholds, considerations of soil moisture and surface cover), lower boundary conditions (roughness, potential dust sources), and dust–radiation interactions. MONARCH now allows modeling of global and regional mineral dust cycles using fundamentally different paradigms, ranging from strongly simplified to physics-based parameterizations. We present a detailed description of these updates along with four global benchmark simulations, which use conceptually different dust emission parameterizations, and we evaluate the simulations against observations of dust optical depth. We determine key dust parameters, such as global annual emission/deposition flux, dust loading, dust optical depth, mass-extinction efficiency, single-scattering albedo, and direct radiative effects. For dust-particle diameters up to 20 µm, the total annual dust emission and deposition fluxes obtained with our four experiments range between about 3500 and 6000 Tg, which largely depend upon differences in the emitted size distribution. Considering ellipsoidal particle shapes and dust refractive indices that account for size-resolved mineralogy, we estimate the global total (longwave and shortwave) dust direct radiative effect (DRE) at the surface to range between about −0.90 and −0.63 W m−2 and at the top of the atmosphere between −0.20 and −0.28 W m−2. Our evaluation demonstrates that MONARCH is able to reproduce key features of the spatiotemporal variability of the global dust cycle with important and insightful differences between the different configurations.

Dust cycle↗

High Solids Loading Aqueous Slurry Formation of Corn Stover Before Pretreatment in a Fed-batch Bioreactor

With the increase on population the world will depend on renewable sources to meet the increasingly energy needs. Use of lignocellulosic biomass as a renewable source has been proven efficient for conversion to cellulosic ethanol and capable of contributing to achieve the threshold on energy demand while reducing greenhouse gases in 90% when compared with fossil fuels (Wang et al., 2007, 2012). However, the processing of biomass encounters limitations in feeding and flow within biorefineries due to the system plugging by compaction and slurry high yield stress, preventing transport of biomass materials and in some cases unexpected plant shutdowns that results in high operational costs (dos Santos et al., 2021; Ximenes et al., 2021). Different solutions for biomass handling and slurry formation from densified materials like pellets have been studied. Cellulosic biomass residues are typically processed in pretreatment reactors to which acid or base is added (Humbird et al., 2002). Other pretreatments such as liquid hot water (LHW) and steam explosion that act without the addition of chemicals use pressures above the saturation vapor of water, to disrupt biomass structure (Ruiz et al., 2021). Although these methods are widely known, the use of pretreatments like acid digestion and LHW present challenges like waste disposal and could be energy inefficient (Mosier et al., 2005). Alternatively other approaches like enzymatic liquefaction have been applied for biomass transformation into slurries with promising results and without the use of pretreatments (dos Santos et al., 2021). Nevertheless, enzymatic liquefaction implementation faces difficulties due to the recalcitrant properties of biomass, its variability, and the release of enzyme inhibitors. While it has been proven that chemical composition properties in biomass have effects on enzyme inhibition (Huang et al., 2022; Kim et al., 2011; Zhai et al., 2018), hindering the efficiency of the liquefaction, not ample research has been conducted in understanding the physical properties as particle size, porosity, water adsorption retention and their effect in the ability to form a slurry at high solids concentration (Yan et al., 2020). In order to improve lignocellulosic biomass handling and formation of biomass slurries, enzyme assisted liquefaction for slurry creation from corn stover at solids loadings up to 30% is reported in this work. Two different kinds of biomass (pelleted corn stover and cobs) were liquefied in a fed-batch process using commercial enzymes Celluclast 1.5L or Ctec-2 at 1FPU or 3 FPU per gram of dry solids in 10 mM sodium citrate buffer solution (pH 4.8). Pellets were fed into a 1 L stirred bioreactor according to a predefined fed-batch protocol over the first 5 hours until reaching 30% of solids loading. After 6, 24 and 96 hours, samples were taken and characterized with respect to their sugar composition, rheology and water absorption. Successful slurry creation with dramatically reduced yield stress was achieved for corn stover for both assessed enzymes. Yield stresses of 178±7 Pa (3 FPU, Celluclast 1.5L) and 79±6 Pa (3 FPU, Ctec-2) were measured for corn stover at 24 hours, compared to 6,000 Pa for samples without enzyme. Yield stress was 155± 29 Pa (3FPU, Ctec-2) and 257 ± 72 Pa (1 FPU, Celluclast 1.5L) for corn cobs at 24 hours. Yield stress decreased when residence time increased with an enhanced fluidity noted for higher enzyme concentrations. A profile for 6, 24 and 96h of yield stress measurements is presented.

Gutierrez, Diana↗

Artificial Neural Networks to Predict Cognitive Impairment of Rodents Subjected to Space Radiation

INTRODUCTION We use artificial neural networks (ANNs) as an example machine learning (ML) tool to predict the cognitive performance impairment of rats induced by irradiation. The experimental data in the analyses is attentional set-shifting (ATSET) test scores from a rodent model exposed to ≤15 cGy of individual galactic cosmic radiation (GCR) ions: 4He, 28Si, or 56Fe, expected for a Lunar or Mars mission [1]. This work investigates rats at a subject-based level and uses applied dose and performance scores taken before irradiation to predict whether a rat will be impaired when irradiated. The results of this study are significant to crewed space missions as they support the potential of predicting an astronaut’s impairment in a specific task before spaceflight through the implementation of appropriately trained ML tools. METHODS Data used in this work are scores from the ATSET, a multi-stage constrained cognitive flexibility test [2]. Our computational model utilizes the number of attempts to reach the criterion to pass a stage as a behavioral performance measure for rats. We use the post-irradiation scores, generate thresholds from cumulative distribution plots of non-irradiated rats, and calculate the percent of irradiated rats whose scores fall below the threshold to infer how each radiation type/dose affects a population. Rats scoring above the threshold are labeled impaired while the others are non-impaired. We then employ ANNs as a typical ML technique, and use each subject’s individual scores taken before radiation along with the applied dose, to predict their personal susceptibility to cognitive impairment due to space radiation exposure. RESULTS AND CONCLUSION A significant finding is the exhibition of a dose-dependent increasing probability of impairment for 1 to 10 cGy of 28Si or 56Fe in the simple discrimination (SD) stage of the ATSET, and for 1 to 10 cGy of 56Fe in the compound discrimination (CD) stage. On a subject-based level, implementing ML classifiers such as ANNs identifies rats that have a higher tendency for impairment after GCR exposure [1]. The receiver operating characteristic (ROC) and the precision-recall (PR) curves of the ML models show a better prediction of impairment when 56Fe is the ion in question in both SD (Figure 1) and CD stages. They, however, do not depict impairment due to 4He in SD (Figure 1) and 28Si in CD, suggesting no dose-dependent impairment response in these cases. In this work, “good” prediction pertains to “better-than-random-chance”, due to the limited sample size and the high inter- and intra-individual variabilities in response to brain stimulation paradigms, as applicable to both animals and humans. More behavioral tests and biomarkers should be investigated on the same subjects, to be fed to the ML models to capture the agents responsible for performance alterations of some individuals versus others.

machine learning↗

Machine Learning Models to Predict Cognitive Impairment of Rodents Subjected to Space Radiation

INTRODUCTION We use artificial neural networks (ANNs) as an example machine learning (ML) tool to predict the cognitive performance impairment of rats induced by irradiation. The experimental data in the analyses is attentional set-shifting (ATSET) test scores from a rodent model exposed to ≤15 cGy of individual galactic cosmic radiation (GCR) ions: 4He, 28Si, or 56Fe, expected for a Lunar or Mars mission [1]. This work investigates rats at a subject-based level and uses applied dose and performance scores taken before irradiation to predict whether a rat will be impaired when irradiated. The results of this study are significant to crewed space missions as they support the potential of predicting an astronaut’s impairment in a specific task before spaceflight through the implementation of appropriately trained ML tools. METHODS Data used in this work are scores from the ATSET, a multi-stage constrained cognitive flexibility test [2]. Our computational model utilizes the number of attempts to reach the criterion to pass a stage as a behavioral performance measure for rats. We use the post-irradiation scores, generate thresholds from cumulative distribution plots of non-irradiated rats, and calculate the percent of irradiated rats whose scores fall below the threshold to infer how each radiation type/dose affects a population. Rats scoring above the threshold are labeled impaired while the others are non-impaired. We then employ ANNs as a typical ML technique, and use each subject’s individual scores taken before radiation along with the applied dose, to predict their personal susceptibility to cognitive impairment due to space radiation exposure. RESULTS AND CONCLUSION A significant finding is the exhibition of a dose-dependent increasing probability of impairment for 1 to 10 cGy of 28Si or 56Fe in the simple discrimination (SD) stage of the ATSET, and for 1 to 10 cGy of 56Fe in the compound discrimination (CD) stage. On a subject-based level, implementing ML classifiers such as ANNs identifies rats that have a higher tendency for impairment after GCR exposure [1]. The receiver operating characteristic (ROC) and the precision-recall (PR) curves of the ML models show a better prediction of impairment when 56Fe is the ion in question in both SD (Figure 1) and CD stages. They, however, do not depict impairment due to 4He in SD (Figure 1) and 28Si in CD, suggesting no dose-dependent impairment response in these cases. In this work, “good” prediction pertains to “better-than-random-chance”, due to the limited sample size and the high inter- and intra-individual variabilities in response to brain stimulation paradigms, as applicable to both animals and humans. More behavioral tests and biomarkers should be investigated on the same subjects, to be fed to the ML models to capture the agents responsible for performance alterations of some individuals versus others.

machine learning↗

A New Approach to Monitoring Coastal Marshes for Persistent Flooding

Many areas in coastal Louisiana are below sea level and protected from flooding by a system of natural and man-made levees. Flooding is common when the levees are overtopped by storm surge or rising rivers. Many levees in this region are further stressed by erosion and subsidence. The floodwaters can become constricted by levees and trapped, causing prolonged inundation. Vegetative communities in coastal regions, from fresh swamp forest to saline marsh, can be negatively affected by inundation and changes in salinity. As saltwater persists, it can have a toxic effect upon marsh vegetation causing die off and conversion to open water types, destroying valuable species habitats. The length of time the water persists and the average annual salinity are important variables in modeling habitat switching (cover type change). Marsh type habitat switching affects fish, shellfish, and wildlife inhabitants, and can affect the regional ecosystem and economy. There are numerous restoration and revitalization projects underway in the coastal region, and their effects on the entire ecosystem need to be understood. For these reasons, monitoring persistent saltwater intrusion and inundation is important. For this study, persistent flooding in Louisiana coastal marshes was mapped using MODIS (Moderate Resolution Imaging Spectroradiometer) time series of a Normalized Difference Water Index (NDWI). The time series data were derived for 2000 through 2009, including flooding due to Hurricane Rita in 2005 and Hurricane Ike in 2008. Using the NDWI, duration and extent of flooding can be inferred. The Time Series Product Tool (TSPT), developed at NASA SSC, is a suite of software developed in MATLAB(R) that enables improved-quality time series images to be computed using advanced temporal processing techniques. This software has been used to compute time series for monitoring temporal changes in environmental phenomena, (e.g. NDVI times series from MODIS), and was modified and used to compute the NDWI indices and also the Normalized Difference Soil Index (NDSI). Coastwide Reference Monitoring System (CRMS) water levels from various hydrologic monitoring stations and aerial photography were used to optimize thresholds for MODIS-derived time series of NDWI and to validate resulting flood maps. In most of the profiles produced for post-hurricane assessment, the increase in the NDWI index (from storm surge) is accompanied by a decrease in the vegetation index (NDVI) and then a period of declining water. The NDSI index represents non-green or dead vegetation and increases after the hurricane s destruction of the marsh vegetation. Behavior of these indices over time is indicative of which areas remain flooded, which areas recover to their former levels of vegetative vigor, and which areas are stressed or in transition. Tracking these indices over time shows the recovery rate of vegetation and the relative behavior to inundation persistence. The results from this study demonstrated that identification of persistent marsh flooding, utilizing the tools developed in this study, provided an approximate 70-80 percent accuracy rate when compared to the actual days flooded at the CRMS stations.

Kalcic, M. T.↗

Comparison of ISCCP and other cloud amounts

A new 8-year global cloud climatology has been produced by the International Satellite Cloud Climatology Project (ISCCP) that provides information every 3 h at 280-km spatial resolution covering the period from July 1983 through June 1991. If cloud detection errors and differences in area sampling are neglected, individual ISCCP cloud amounts agree with individual surface observations to within 15% rms with biases of only a few percent. When measurements of small-scale, broken clouds are isolated in the comparison, the rms differences between satellite and surface cloud amounts are about 25%, similar to the rms difference between ISCCP and Landsat determinations of cloud amount. For broken clouds, the average ISCCP cloud amounts are about 5% smaller than estimated by surface observers (difference between earth cover and sky cover), but about 5% larger than estimated from very high spatial resolution satellite observations (overestimate due to low spatial resolution offset by underestimate due to finite radiance thresholds). Detection errors caused by errors in the ISCCP clear-sky radiances or incorrect radiance threshold magnitudes, are the dominant source of error in monthly average cloud amounts. The ISCCP cloud amounts appear to be too low over land by about 10%, somewhat less in summer and somewhat more in winter, and about right (maybe slightly low) over oceans. In polar regions, ISCCP cloud amounts are probably too low by about 15%-25% in summer and 5%-10% in winter. Comparison of the ISCCP climatology to three other cloud climatologies shows excellent agreement in the geographic distribution and seasonal variation of cloud amounts; there is little agreement about day/night contrasts in cloud amount. Notable results from ISCCP are that the gobal annual mean cloud amount is about 63%, being about 23% higher over oceans than over land, that it varies by less than 1% rms from month to month, and that it has varied by about 4% on a time scale approximately equal to 2-4 years. The magnitude of interannual variations of local (280-km scale) monthly mean cloud amounts is about 7%-9%. Longitudinal contrasts in cloud amount are just as large as latitudinal contrasts. The largest seasonal variation of cloud amount occurs in the tropics, being larger in summer than in winter; the seasonal variation in the middle latitudes has the opposite phase. Polar regions may have little seasonable variability in cloud amount. The ISCCP results show slightly more nightime than daytime cloud amount over oceans and more daytime than nightime cloud amount over land.

Rossow, William B.↗

Finite boundary effects on the spherical Rayleigh–Taylor instability between viscous fluids

For the Rayleigh–Taylor unstable arrangement of a viscous fluid sphere embedded in a finite viscous fluid spherical shell with a rigid boundary and a radially directed acceleration, a dispersion relation is developed from a linear stability analysis using the method of normal modes. [Formula: see text] is the radially directed acceleration at the interface. ρ i denotes the density, μ i is the viscosity, and R i is the radius, where i = 1 is the inner sphere and i = 2 is the outer sphere. The dispersion relation is a function of the following dimensionless variables: viscosity ratio [Formula: see text], density ratio [Formula: see text], spherical harmonic mode n, [Formula: see text], [Formula: see text], and the dimensionless growth rate [Formula: see text], where σ is the exponential growth rate. We show that the boundedness provided by the outer spherical shell has a strong influence on the instability behavior, which is reflected not only in the modulation of the growth rate but also in the selection of the most unstable modes that are physically possible. This outer boundary effect is quantified by the relative magnitude of the radius ratio H. We find that when H is close to unity, lower order harmonics are excluded from becoming the most unstable within a vast region of the parameter space. In other words, the effect of H has precedence over the other controlling parameters d, B, and a wide range of s in establishing what the lowest most unstable mode can be. When H ~ 1, low order harmonics can become the most unstable only for s >> 1. However, in the limit when s → ∞, we show that the most unstable mode is n = 1 and derive the dispersion relation in this limit. The exclusion of most unstable low order harmonics caused by a finite outer boundary is not realized when the outer boundary extends beyond a certain threshold length-scale in which case all modes are equally possible depending on the value of B.

74 ATOMIC AND MOLECULAR PHYSICS↗

Custom-trained Machine-learning Interatomic Potentials: ZnCl2 Aqueous Solution

This dataset was generated using an iterative active-learning strategy implemented in the ArcaNN software package (https://github.com/arcann-chem/arcann_training) to train machine-learning interatomic potentials for aqueous ZnCl2 solutions. Each active-learning cycle consisted of three stages: training, exploration, and labeling. The initial training set combined configurations generated in this work from enhanced-sampling ab initio molecular dynamics simulations with configurations from a previously reported neural-network-potential study of aqueous ZnCl2. The enhanced-sampling ab initio molecular dynamics simulations involved Zn–Cl separation and the chloride coordination number around Zn²? as collective variables. These configurations served as the seed dataset. Subsequent active-learning cycles expanded the training set by identifying and labeling configurations that were poorly represented by the current models, thereby improving coverage of ion-association states and changes in local coordination and charge-state environments relevant to the solution free-energy landscape. For all selected configurations, single-point calculations of the total energies and atomic forces were performed within density functional theory using the CP2K Quickstep module. Reference calculations employed the revPBE-D3 and r2SCAN exchange-correlation functionals. Motivated by recent work on aqueous Zn²?, the main revPBE calculations omitted D3 dispersion contributions involving Zn²?, while retaining the D3 correction for water and chloride. For comparison, fully dispersion-corrected revPBE-D3 reference calculations were also performed, with D3 applied to all species, including Zn²?. Valence electrons were treated explicitly, while core electrons were represented using norm-conserving Goedecker–Teter–Hutter pseudopotentials. The wave functions were expanded using the mixed Gaussian-and-plane-wave scheme with TZV2P-MOLOPT basis sets for all elements and a 600 Ry auxiliary plane-wave cutoff for the electron density. Self-consistent-field convergence was accelerated using the orbital-transformation and Direct Inversion in the Iterative Subspace algorithms, with a convergence threshold of 10?6. All single-point calculations were performed in periodic orthorhombic cells. The CELL_REF keyword in CP2K was used to define a fixed reference cell with a box length of 25 Å. This treatment ensured a consistent reference for configurations extracted from NpT trajectories with fluctuating cell dimensions. The resulting DFT energies and atomic forces constitute the ground-truth labels used to train the MLIPs. The resulting MLIP was trained for aqueous ZnCl2 solutions spanning concentrations from 0 to 30 molal and a broad pH range, from strongly acidic to strongly basic conditions. Representative examples of configurations included in the MLIP training dataset are provided below. These include 1) Representative configurations from the dataset labeled at the revPBE-D3 level, with D3 dispersion interactions involving Zn2+ excluded (revPBE-wo-D3). 2) Representative configurations from the dataset labeled at the fully dispersion-corrected revPBE-D3 level, with D3 interactions applied to all species, including Zn2+ (revPBE-D3). 3) Representative configurations from the dataset labeled at the r2SCAN level of theory (r2SCAN).

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗

Cost of Fish Exclusion and Passage Technologies for Hydropower

Hydropower represents a reliable source of renewable energy and accounts for approximately 7% of the total electrical generation in the United States. Future expansion of hydropower is likely to be in the form of either smaller new stream development projects or powering existing non-powered dams. For these new projects to be successful, careful analysis of risks, costs, and uncertainty to offset reduced power production as well as ensuring the protection and safe passage of migratory fish to gain public support, will be required. Exclusion and passage are two common approaches to protect fish from entrainment and impingement at hydropower facilities. The thresholds for entrainment risk and requirements for exclusion and passage often differ depending on the species involved, the characteristics of the facility, and the goals of stakeholders. While the costs associated with environmental mitigations represent a large proportion of the total costs required for the licensing of hydropower facilities, little quantitative information is present within the literature regarding the specific costs of fish exclusion and passage. Working with FOA awardee Natel Energy, scientists at Oak Ridge National Laboratory were tasked with assessing the capital construction costs for downstream fish exclusion and passage infrastructure. This report used keyword searches of an existing environmental mitigation cost data set and manual extraction of additional cost data associated with protection, mitigation, and enhancement (PM&E) measures related to positive barrier screening and passage from regulatory licensing documents available in the Federal Energy Regulatory Commission (FERC) eLibrary. This approach yielded a total of 50 PM&E mitigation measures with estimated capital construction costs pertaining to positive barrier screens, 142 pertaining to passage studies, and 26 pertaining to passage-related studies. PM&E measures associated with positive barrier screens represented <10% of the 171 total FERC project dockets available in the data set. These data were highly skewed toward conventional relicensing projects, as <7% were associated with new stream development (NSD) projects. Results from these data indicate highly variable costs associated with fish screening, with flow-normalized costs one to two orders of magnitude higher for screening with the highest exclusion capability (≤0.09 in. spacing) compared with coarser screening (1 to 2 in.). Furthermore, estimated capital costs of passage infrastructure were positively related to the scale of the project based on installed capacity for some, but not all, types of passage. These data provide an initial baseline for estimating exclusion and passage costs for hydropower development and may help developers consider options for more fish-friendly generation technologies, though gaps remain relating to a lack of data, particularly for NSD projects. More data may still be available within the FERC eLibrary, but significant effort will be required to manually identify and extract the data for future analyses.

13 HYDRO ENERGY↗

The role of interannual ENSO events in decadal timescale transitions of the Interdecadal Pacific Oscillation

The build-up of decadal timescale upper ocean heat content in the off-equatorial western tropical Pacific can provide necessary conditions for interannual El Niño/Southern Oscillation (ENSO) events to contribute to decadal timescale transitions of tropical Pacific SSTs to the opposite phase of the Interdecadal Pacific Oscillation (IPO). This can be viewed as a corollary to subseasonal westerly wind burst events contributing to El Niño interannual timescale transitions. In this work, a long pre-industrial control run with CESM1 is analyzed to show that there is a greater chance of ENSO activity to contribute to an IPO transition when off-equatorial western Pacific Ocean heat content reaches either a maximum (for El Niño to contribute to a transition to positive IPO) or minimum (for La Niña to contribute to a transition to negative IPO) as seen in observations. If above a necessary ocean heat content threshold, the convergence associated with westerly anomaly near-equatorial surface winds associated with El Niño activity can draw that heat content equatorward to sustain anomalously warm western and central Pacific SSTs. These are associated with positive precipitation and convective heating anomalies, a Gill-type response and wind stress curl anomalies that continue to feed warm water into the near-equatorial western Pacific. These conditions then sustain the decadal-timescale transition to positive IPO (with the opposite sign for transition to negative IPO). Associated central equatorial Pacific convective heating anomalies produce SLP and wind stress anomalies in the North and South Pacific that can excite westward-propagating off-equatorial oceanic Rossby waves to contribute to the western Pacific thermocline depth and consequent heat content anomalies.

54 ENVIRONMENTAL SCIENCES↗

Long-Lag, Wide-pulse Gamma-Ray Bursts

Currently, the best available probe of the early phase of gamma-ray burst (GRB) jet attributes is the prompt gamma-ray emission, in which several intrinsic and extrinsic variables determine GRB pulse evolution. Bright, usually complex bursts have many narrow pulses that are difficult to model due to overlap. However, the relatively simple, long spectral lag, wide-pulse bursts may have simpler physics and are easier to model. In this work we analyze the temporal and spectral behavior of wide pulses in 24 long-lag bursts, using a pulse model with two shape parameters - width and asymmetry - and the Band spectral model with three shape parameters. We find that pulses in long-lag bursts are distinguished both temporally and spectrally from those in bright bursts: the pulses in long spectral lag bursts are few in number, and approximately 100 times wider (10s of seconds), have systematically lower peaks in vF(v), harder low-energy spectra and softer high-energy spectra. We find that these five pulse descriptors are essentially uncorrelated for our long-lag sample, suggesting that at least approximately 5 parameters are needed to model burst temporal and spectral behavior. However, pulse width is strongly correlated with spectral lag; hence these two parameters may be viewed as mutual surrogates. We infer that accurate formulations for estimating GRB luminosity and total energy will depend on several gamma-ray attributes, at least for long-lag bursts. The prevalence of long-lag bursts near the BATSE trigger threshold, their predominantly low vF(v) spectral peaks, and relatively steep upper power-law spectral indices indicate that Swift will detect many such bursts.

Norris, J. P.↗

Fault-Tolerant Operation of Bosonic Qubits with Discrete-Variable Ancillae

Fault-tolerant quantum computation with bosonic qubits often necessitates the use of noisy discrete-variable ancillae. In this work, we establish a comprehensive and practical fault-tolerance framework for such a hybrid system and synthesize it with fault-tolerant protocols by combining bosonic quantum error correction (QEC) and advanced quantum control techniques. We introduce essential building blocks of error-corrected gadgets by leveraging ancilla-assisted bosonic operations using a generalized variant of path-independent quantum control. Using these building blocks, we construct a universal set of error-corrected gadgets that tolerate a single-photon loss and an arbitrary ancilla fault for four-legged cat qubits. Notably, our construction requires only dispersive coupling between bosonic modes and ancillae, as well as beam-splitter coupling between bosonic modes, both of which have been experimentally demonstrated with strong strengths and high accuracy. Moreover, each error-corrected bosonic qubit is comprised of only a single bosonic mode and a three-level ancilla, featuring the hardware efficiency of bosonic QEC in the full fault-tolerant setting. We numerically demonstrate the feasibility of our schemes using current experimental parameters in the circuit-QED platform. Finally, we present a hardware-efficient architecture for fault-tolerant quantum computing by concatenating the four-legged cat qubits with an outer qubit code utilizing only beam-splitter couplings. Our estimates suggest that the overall noise threshold can be reached using existing hardware. These developed fault-tolerant schemes extend beyond their applicability to four-legged cat qubits and can be adapted for other rotation-symmetrical codes, offering a promising avenue toward scalable and robust quantum computation with bosonic qubits. Published by the American Physical Society 2024

Physics↗

Assessing the Viability of Using GEOS-Forecast Product for Landslides Forecasting: A Step Toward Early Warning System

Landslides across the globe are mostly triggered by extreme rainfall events affecting infrastructure, transportation and livelihoods. The risks are rarely quantified due to lack of data, analytical skills and limited modeling techniques. Knowledge of local to global scale landslide risks provides communities and national agencies the ability to adapt disaster management practices to mitigate and recover from these hazards. In order to minimize the risks and improve characterization of community resilience to landslides, it is vital to have reliable information about the factors triggering landslides such as rainfall, well ahead in time. Forecasting potential landslide activity and impacts can be achieved through reliable precipitation forecast models. However, it is challenging because of the temporal and spatial variability of precipitation, an important factor in triggering landslides. Evaluation of the precipitation field, associated errors, and sampling uncertainties is integral for development of efficient and reliable landslide forecasting and early warning system. This study develops a methodology to assess the viability of using a precipitation field provided by a global model and its potential integration in the landslide forecasting system. The study focuses on the comparison between the IMERG (Integrated Multi-satellitE Retrievals for Global Precipitation Mission) and GEOS (NASA Goddard Earth Observing System)-Forecast product over contiguous United States (CONUS). GEOS model assimilates new observations every 6 hours, at 00, 06, 12, and 18 UTC. The framework is tested on the GEOS-Forecast Model initialized at 00 UTC using daily IMERG early product as reference using both categorical and continuous statistics. The categorical statistics includes the probability of detection (POD), success ratio (SR), critical success index (CSI), and the hit bias. Continuous statistics such as correlation, normalized standard deviation, and root-mean-square error are also evaluated. Overall, GEOS-Forecast precipitation field over the analysis period (~1 year) show underestimation with respect to IMERG early for the daily accumulated rainfall. However, the probability distribution function and cumulative distribution function of both show similar patterns. In terms of correlations, POD, SR, CSI, hit bias, the performance varies with respect to the rainfall threshold used.

Sana Khan↗

Long-Lag, Wide-pulse Gamma-Ray Bursts

The best available probe of the early phase of gamma-ray burst (GRB) jet attributes is the prompt gamma-ray emission, in which several intrinsic and extrinsic variables determine observed GRB pulse evolution, including at least: jet opening angle, profiles of Lorentz factor and matter/field density, distance of emission region from central source, and viewing angle. Bright, usually complex bursts have many narrow pulses that are difficult to model due to overlap. However, the relatively simple, long spectral lag, wide-pulse bursts may have simpler physics and are easier to model. We have analyzed the temporal and spectral behavior of wide pulses in 24 long-lag bursts from the BATSE sample, using a pulse model with two shape parameters - width and asymmetry - and the Band spectral model with three shape parameters. We find that pulses in long-lag bursts are distinguished both temporally and spectrally from those in bright bursts: the pulses in long spectral lag bursts are few in number, and approximately 100 times wider (10s of seconds), have systemtically lower peaks in nu*F(nu), harder low-energy spectra and softer high-energy spectra. These five pulse descriptors are essentially uncorrelated for our long-lag sample, suggesting that at least approximately 5 parameters are needed to model burst temporal and spectral behavior, roughly commensurate with the theoretical phase space. However, we do find that pulse width is strongly correlated with spectral lag; hence these two parameters may be viewed as mutual surrogates. The prevalence of long-lag bursts near the BATSE trigger threshold, their predominantly low nu*F(nu) spectral peaks, and relatively steep upper power-law spectral indices indicate that Swiift will detect many such bursts.

Norris, J. P.↗

Lower soil moisture and deep soil temperatures in thermokarst features increase old soil carbon loss after 10 years of experimental permafrost warming

Almost half of the global terrestrial soil carbon (C) is stored in the northern circumpolar permafrost region, where air temperatures are increasing two times faster than the global average. As climate warms, permafrost thaws and soil organic matter becomes vulnerable to greater microbial decomposition. Long-term soil warming of ice-rich permafrost can result in thermokarst formation that creates variability in environmental conditions. Consequently, plant and microbial proportional contributions to ecosystem respiration may change in response to long-term soil warming. Natural abundance δ 13 C and Δ 14 C of aboveground and belowground plant material, and of young and old soil respiration were used to inform a mixing model to partition the contribution of each source to ecosystem respiration fluxes. We employed a hierarchical Bayesian approach that incorporated gross primary productivity and environmental drivers to constrain source contributions. We found that long-term experimental permafrost warming introduced a soil hydrology component that interacted with temperature to affect old soil C respiration. Old soil C loss was suppressed in plots with warmer deep soil temperatures because they tended to be wetter. When soil volumetric water content significantly decreased in 2018 relative to 2016 and 2017, the dominant respiration sources shifted from plant aboveground and young soil respiration to old soil respiration. The proportion of ecosystem respiration from old soil C accounted for up to 39% of ecosystem respiration and represented a 30-fold increase compared to the wet-year average. Our findings show that thermokarst formation may act to moderate microbial decomposition of old soil C when soil is highly saturated. However, when soil moisture decreases, a higher proportion of old soil C is vulnerable to decomposition and can become a large flux to the atmosphere. Furthermore, as permafrost systems continue to change with climate, we must understand the thresholds that may propel these systems from a C sink to a source.

permafrost↗

Spiking Neurons for Analysis of Patterns

Artificial neural networks comprising spiking neurons of a novel type have been conceived as improved pattern-analysis and pattern-recognition computational systems. These neurons are represented by a mathematical model denoted the state-variable model (SVM), which among other things, exploits a computational parallelism inherent in spiking-neuron geometry. Networks of SVM neurons offer advantages of speed and computational efficiency, relative to traditional artificial neural networks. The SVM also overcomes some of the limitations of prior spiking-neuron models. There are numerous potential pattern-recognition, tracking, and data-reduction (data preprocessing) applications for these SVM neural networks on Earth and in exploration of remote planets. Spiking neurons imitate biological neurons more closely than do the neurons of traditional artificial neural networks. A spiking neuron includes a central cell body (soma) surrounded by a tree-like interconnection network (dendrites). Spiking neurons are so named because they generate trains of output pulses (spikes) in response to inputs received from sensors or from other neurons. They gain their speed advantage over traditional neural networks by using the timing of individual spikes for computation, whereas traditional artificial neurons use averages of activity levels over time. Moreover, spiking neurons use the delays inherent in dendritic processing in order to efficiently encode the information content of incoming signals. Because traditional artificial neurons fail to capture this encoding, they have less processing capability, and so it is necessary to use more gates when implementing traditional artificial neurons in electronic circuitry. Such higher-order functions as dynamic tasking are effected by use of pools (collections) of spiking neurons interconnected by spike-transmitting fibers. The SVM includes adaptive thresholds and submodels of transport of ions (in imitation of such transport in biological neurons). These features enable the neurons to adapt their responses to high-rate inputs from sensors, and to adapt their firing thresholds to mitigate noise or effects of potential sensor failure. The mathematical derivation of the SVM starts from a prior model, known in the art as the point soma model, which captures all of the salient properties of neuronal response while keeping the computational cost low. The point-soma latency time is modified to be an exponentially decaying function of the strength of the applied potential. Choosing computational efficiency over biological fidelity, the dendrites surrounding a neuron are represented by simplified compartmental submodels and there are no dendritic spines. Updates to the dendritic potential, calcium-ion concentrations and conductances, and potassium-ion conductances are done by use of equations similar to those of the point soma. Diffusion processes in dendrites are modeled by averaging among nearest-neighbor compartments. Inputs to each of the dendritic compartments come from sensors. Alternatively or in addition, when an affected neuron is part of a pool, inputs can come from other spiking neurons. At present, SVM neural networks are implemented by computational simulation, using algorithms that encode the SVM and its submodels. However, it should be possible to implement these neural networks in hardware: The differential equations for the dendritic and cellular processes in the SVM model of spiking neurons map to equivalent circuits that can be implemented directly in analog very-large-scale integrated (VLSI) circuits.

Huntsberger, Terrance↗

Modifications to Solar Titan-130 Combustion Systems for Efficient, High Turndown Operation

The project team of Southwest Research Institute® (SwRI®), Solar Turbines Incorporated (Solar), the Electric Power Research Institute (EPRI), the University of California, Irvine (UCI), and the Georgia Institute of Technology (Georgia Tech) investigated methods to allow higher efficiency part-load operation of a Solar Titan 130 gas turbine. The objective was to develop a low-emission combustion system capable of sustaining combustion and avoiding lean blowout during high turndown operation, which would allow the gas turbine to operate as efficiently as possible at part load. Currently, electric utility markets are beginning to experience substantial increases in renewable energy generation. Some of these renewable energy sources have highly variable output in an uncontrolled manner. In order to maintain grid stability, there is a need for power plants to ramp up power to the grid rapidly to make up for drops in renewable generation. This is often termed spinning reserve, but the size of this reserve may need to increase as renewable penetration into the electric utility market increases. Small combined heat and power (CHP) power plants provide a promising option for meeting this spinning reserve requirement. In order to operate in spinning reserve while still meeting the heat requirements for the CHP, the gas turbine needs to operate efficiently at very low loads. Efficient, high turndown operations in this engine are limited by the lean flammability limit of the premixed combustion system. This project sought enhance the lean operability range of the Titan 130 combustor. First, the project team participated in a brainstorming activity and ultimately selected two concepts to explore: fuel augmentation with hydrogen (H2) to improve the stability at lean operating conditions and modifications to the fuel nozzle to improve the emissions performance at lean operating conditions. Analytical and laboratory investigations were accomplished by UCI to investigate the efficacy of H2 addition at improving lean blow out (LBO) limits and the resulting emissions. These investigations used a variety of chemical reactor network (CRN) and CFD models, validated against laboratory data, to model the impact of H 2 and inform the experimental efforts accomplished by SwRI and Solar. Ultimately, both the CRN and CFD models yielded generally good agreement with the experimental data below a particular temperature threshold. Atmospheric tests of a full-scale T130 annular combustor were performed at SwRI facilities in San Antonio, Texas, to investigate the use of H 2 addition. For these tests, the T130 combustion system remained largely unchanged; minor modifications were performed to the fuel ducting to allow for the safe use of H 2 . The test ultimately demonstrated that the addition of H 2 to the fuel mixture significantly increased the AFR ratio at which the combustor could operate. This improvement to the LBO limit should allow for less use of compressor bleed and less throttling needed by the inlet guide vanes (IGV). This in turn could result in more efficient operation of the gas turbine at lower load points. The second modification explored in this work was a direct modification to the T130 injector. The project team hypothesized that modifications to the pilot of the T130 injector could provide lower emissions at high turn-down operations. These modifications were manufactured and explored by the team at Solar. High pressure rig tests, originally slated to occur at SwRI, were ultimately accomplished by Solar to maintain overall project budget and mitigate cost growth attributable to supply chain issues and inflation. The pressurized rig tests ultimately showed that the SwRI Project No. 18.24153 - DE-EE0008415 Page 2 Final Technical Report January 24, 2024 modifications did not significantly alter the performance of the combustion system at the high turn-down conditions; both the modified injectors and the baseline configuration exhibited elevated emissions comparted to the full-load operating condition. A final set of studies performed by EPRI investigated the benefit-cost of flexible CHP as well as a grid interconnection study for the California Independent System Operator (CAISO) grid. These studies considered: traditional CHP with no spinning reserve available for on-demand grid support, 50% flexible CHP where 50% of the machine’s capacity is consumed by on-site baseload operations while providing an additional 50% capacity for on-demand grid support, and 70% flexible CHP where 70% of capacity is consumed on-site by baseload operations and 30% is available for on-demand grid support. In all cases, the analyses showed a benefit-to-cost ratio greater than unity implying a positive net present value for all configurations. However, the traditional CHP showed the most economic benefit. These results are sensitive to several factors, many of which are not fully known and may vary over time. Thus site owners must be convinced that taking up the increased costs and risks from flexible CHP would be worth implementing. As the grid in California and across the country transition to incorporate larger renewable energy generation, flexible CHP can provide much needed operating reserves and dispatchability. Alternative fuel options, such as hydrogen blending and biofuels, may also lower carbon intensities of CHP. Flexible CHP should be examined in the evolving market to understand innovative business models, changes market rules and services, and new technologies.

20 FOSSIL-FUELED POWER PLANTS↗

Studies of Antarctic Sea Ice Concentrations from Satellite Data and Their Applications

Large changes in the sea ice cover have been observed recently. Because of the relevance of such changes to climate change studies it is important that key ice concentration data sets used for evaluating such changes are interpreted properly. High and medium resolution visible and infrared satellite data are used in conjunction with passive microwave data to study the true characteristics of the Antarctic sea ice cover, assess errors in currently available ice concentration products, and evaluate the applications and limitations of the latter in polar process studies. Cloud-free high resolution data provide valuable information about the natural distribution, stage of formation, and composition of the ice cover that enables interpretation of the large spatial and temporal variability of the microwave emissivity of Antarctic sea ice. Comparative analyses of co-registered visible, infrared and microwave data were used to evaluate ice concentrations derived from standard ice algorithms (i.e., Bootstrap and Team) and investigate the 10 to 35% difference in derived values from large areas within the ice pack, especially in the Weddell Sea, Amundsen Sea, and Ross Sea regions. Landsat and OLS data show a predominance of thick consolidated ice in these areas and show good agreement with the Bootstrap Algorithm. While direct measurements were not possible, the lower values from the Team Algorithm results are likely due to layering within the ice and snow and/or surface flooding, which are known to affect the polarization ratio. In predominantly new ice regions, the derived ice concentration from passive microwave data is usually lower than the true percentage because the emissivity of new ice changes with age and thickness and is lower than that of thick ice. However, the product provides a more realistic characterization of the sea ice cover, and are more useful in polar process studies since it allows for the identification of areas of significant divergence and polynya activities. Also, heat and salinity fluxes are proportionately increased in these areas compared to those from the thicker ice areas. A slight positive trend in ice extent and area from 1978 through 2000 is observed consistent with slight continental cooling during the period. However, the confidence in this result is only moderate because the overlap period for key instruments is just one month and the sensitivity to changes in sensor characteristics, calibration and threshold for the ice edge is quite high.

Comiso, Josefino C.↗