Search NASA⌕ Search

SEARCH · Search NASA

Results for “test vector”

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 451 records · Page 25

A measurement of $H_0$ from DESI DR1 using energy densities

We present a new measurement of the Hubble constant, independent of standard rulers and robust to pre-recombination modifications such as Early Dark Energy (EDE), obtained by calibrating the total energy density of the Universe. We start using the present-day photon density as an anchor, and use the baryon-to-photon ratio from Big Bang Nucleosynthesis based measurements and the baryon-to-matter ratio from the baryons' imprint on galaxy clustering to translate to a physical matter density at present day. We then compare this to measurements of the ratio of the matter density to the critical density ($Ω_{\mathrm{m}}$), calculated using the relative positions of the baryon acoustic oscillations, to measure the critical density of the universe and hence $H_0$. The important measurements of the evolution of the energy density all happen at low redshift, so we consider this a low-redshift measurement. We validate our method both on a suite of $N$-body mocks and on noiseless theory vectors generated across a wide range of Hubble parameters in both $Λ$CDM and EDE cosmologies. Using DESI DR1 data combined with the angular CMB acoustic scale and the latest BBN constraints, we find $H_0 = 69.0 \pm 2.5$ km s$^{-1}$ Mpc$^{-1}$, consistent with existing early and late-time determinations of the Hubble constant. We consider the impact of non-standard dark energy evolution on our measurement. Future data, including that from further iterations of DESI and from Euclid, will add to these results providing a powerful test of the Hubble tension.

Krolewski, Alex [U. Waterloo (main); Waterloo U., ↗

First M87 Event Horizon Telescope Results. V. Physical Origin of the Asymmetric Ring

The Event Horizon Telescope (EHT) has mapped the central compact radio source of the elliptical galaxy M87 at 1.3 mm with unprecedented angular resolution. Here we consider the physical implications of the asymmetric ring seen in the 2017 EHT data. To this end, we construct a large library of models based on general relativistic magnetohydrodynamic (GRMHD) simulations and synthetic images produced by general relativistic ray tracing. We compare the observed visibilities with this library and confirm that the asymmetric ring is consistent with earlier predictions of strong gravitational lensing of synchrotron emission from a hot plasma orbiting near the black hole event horizon. The ring radius and ring asymmetry depend on black hole mass and spin, respectively, and both are therefore expected to be stable when observed in future EHT campaigns. Overall, the observed image is consistent with expectations for the shadow of a spinning Kerr black hole as predicted by general relativity. If the black hole spin and M87's large scale jet are aligned, then the black hole spin vector is pointed away from Earth. Models in our library of non-spinning black holes are inconsistent with the observations as they do not produce sufficiently powerful jets. At the same time, in those models that produce a sufficiently powerful jet, the latter is powered by extraction of black hole spin energy through mechanisms akin to the Blandford-Znajek process. We briefly consider alternatives to a black hole for the central compact object. Analysis of existing EHT polarization data and data taken simultaneously at other wavelengths will soon enable new tests of the GRMHD models, as will future EHT campaigns at 230 and 345 GHz.

79 ASTRONOMY AND ASTROPHYSICS↗

ModuleOT: A Hardware Security Module for Operational Technology: Preprint

With increasing penetration levels of distributed energy resources (DERs) on the distribution grid, as well as new technological advancements in the cyber space, new cyberattack vectors are being introduced, and the available attack surface is constantly increasing. Despite this increasing risk, the standard IEEE 1547-2018 does not yet recommend cybersecurity measures for DERs. To address this and to better protect data on the distribution grid - from the standpoints information security as well as operational security - ModuleOT has been developed. The module aims to significantly reduce cyberattack vectors by improving data privacy for user applications. This is accomplished by performing the core functions of encryption, authentication, authorization, certificate management, and user access control. The module integrates a custom security application with hardware cryptographic acceleration. The application secures all communications using Transmission Control Protocol over Internet Protocol (TCP/IP). These include the three most commonly used communications protocols for power systems information exchange: Modbus, Distributed Network Protocol 3 (DNP3), and Smart Energy Profile 2.0 (SEP2.0). These three protocols are also supported by IEEE 1547-2018 for all DER devices. This paper tests the data encryption/decryption feature on a physical networking test bed with emulated Modbus devices reporting grid data and presents the results.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Glancing-incidence focussed ion beam milling: A coherent X-ray diffraction study of 3D nano-scale lattice strains and crystal defects

Here, this study presents a detailed examination of the lattice distortions introduced by glancing incidence Focussed Ion Beam (FIB) milling. Using non-destructive multi-reflection Bragg coherent X-ray diffraction we probe damage formation in an initially pristine gold micro-crystal following several stages of FIB milling. These experiments allow access to the full lattice strain tensor in the micro-crystal with ~25 nm 3D spatial resolution, enabling a nano-scale analysis of residual lattice strains and defects formed. Our results show that 30 keV glancing incidence milling produces fewer large defects than normal incidence milling at the same energy. However the resulting residual lattice strains have similar magnitude and extend up to ~50 nm into the sample. At the edges of the milled surface, where the ion-beam tails impact the sample at near-normal incidence, large dislocation loops with a range of Burgers vectors are formed. Further glancing incidence FIB polishing with 5 keV ion energy removes these dislocation loops and reduces the lattice strains caused by higher energy FIB milling. However, even at the lower ion energy, damage-induced lattice strains are present within a ~20 nm thick surface layer. These results highlight the need for careful consideration and management of FIB damage. They also show that low-energy FIB-milling is an effective tool for removing FIB-milling induced lattice strains. This is important for the preparation of micro-mechanical test specimens and strain microscopy samples.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A compute-bound formulation of Galerkin model reduction for linear time-invariant dynamical systems

This work aims to advance computational methods for projection-based reduced-order models (ROMs) of linear time-invariant (LTI) dynamical systems. For such systems, current practice relies on ROM formulations expressing the state as a rank-1 tensor (i.e., a vector), leading to computational kernels that are memory bandwidth bound and, therefore, ill-suited for scalable performance on modern architectures. This weakness can be particularly limiting when tackling many-query studies, where one needs to run a large number of simulations. This work introduces a reformulation, called rank-2 Galerkin, of the Galerkin ROM for LTI dynamical systems which converts the nature of the ROM problem from memory bandwidth to compute bound. We present the details of the formulation and its implementation, and demonstrate its utility through numerical experiments using, as a test case, the simulation of elastic seismic shear waves in an axisymmetric domain. We quantify and analyze performance and scaling results for varying numbers of threads and problem sizes. In conclusion, we present an end-to-end demonstration of using the rank-2 Galerkin ROM for a Monte Carlo sampling study. We show that the rank-2 Galerkin ROM is one order of magnitude more efficient than the rank-1 Galerkin ROM (the current practice) and about 970 times more efficient than the full-order model, while maintaining accuracy in both the mean and statistics of the field.

97 MATHEMATICS AND COMPUTING↗

Clostridioides difficile Toxin A Remodels Membranes and Mediates DNA Entry Into Cells to Activate Toll-Like Receptor 9 Signaling

Background & Aims: Clostridioides difficile toxin A (TcdA) activates the innate immune response. TcdA co-purifies with DNA. Toll-like receptor 9 (TLR9) recognizes bacterial DNA to initiate inflammation. We investigated whether DNA bound to TcdA activates an inflammatory response in murine models of C difficile infection via activation of TLR9. Methods: We performed studies with human colonocytes and monocytes and macrophages from wild-type and TLR9 knockout mice incubated with TcdA or its antagonist (ODN TTAGGG) or transduced with vectors encoding TLR9 or small-interfering RNAs. Cytokine production was measured with enzyme-linked immunosorbent assay. We studied a transduction domain of TcdA (TcdA 57-80 ), which was predicted by machine learning to have cell-penetrating activity and confirmed by synchrotron small-angle X-ray scattering. Intestines of CD1 mice, C57BL6J mice, and mice that express a form of TLR9 that is not activated by CpG DNA were injected with TcdA, TLR9 antagonist, or both. Enterotoxicity was estimated based on loop weight to length ratios. A TLR9 antagonist was tested in mice infected with C difficile. We incubated human colon explants with an antagonist of TLR9 and measured TcdA-induced production of cytokines. Results: The TcdA 57-80 protein transduction domain had membrane remodeling activity that allowed TcdA to enter endosomes. TcdA-bound DNA entered human colonocytes. TLR9 was required for production of cytokines by cultured cells and in human colon explants incubated with TcdA. TLR9 was required in TcdA-induced mice intestinal secretions and in the survival of mice infected by C difficile. Even in a protease-rich environment, in which only fragments of TcdA exist, the TcdA 57-80 domain organized DNA into a geometrically ordered structure that activated TLR9. Conclusions: TcdA from C difficile can bind and organize bacterial DNA to activate TLR9. TcdA and TcdA fragments remodel membranes, which allows them to access endosomes and present bacterial DNA to and activate TLR9. Rather than inactivating the ability of DNA to bind TLR9, TcdA appears to chaperone and organize DNA into an inflammatory, spatially periodic structure.

59 BASIC BIOLOGICAL SCIENCES↗

Accurate Machine Learning for Predicting the Viscosities of Deep Eutectic Solvents

Deep eutectic solvents (DESs) are emerging as environmentally friendly designer solvents for mass transport and heat transfer processes in industrial applications; however, the lack of accurate tools to predict and thus control their viscosities under both a range of environmental factors and formulations hinders their general application. While DESs may serve as designer solvents, with nearly unlimited combinations, this unfortunately makes it experimentally infeasible to comprehensively measure the viscosities of all DESs of potential industrial interest. To assist in the design of DESs, we have developed several new machine learning (ML) models that accurately and rapidly predict the viscosities of a diverse group of DESs at different temperatures and molar ratios using, to date, one of the most comprehensive data sets containing the properties of over 670 DESs over a wide range of temperatures (278.15–385.25 K). Three ML models, including support vector regression (SVR), feed forward neural networks (FFNNs), and categorical boosting (CatBoost), were developed to predict DES viscosity as a function of temperature and molar ratio and contrasted with multilinear and two-factor polynomial regression baselines. Further, quantum chemistry-based, COSMO-RS-derived sigma profile (σ-profile) features were used as inputs for the ML models. The CatBoost model is excellent at externally predicting DES viscosity, as indicated by high R 2 (0.99) and low root-mean-square-error (RMSE) and average absolute relative deviations (AARD) (5.22%) values for the testing data sets, and 98% of the data points lie within the 15% of AARD deviations. Furthermore, SHapley additive explanation (SHAP) analysis was employed to interpret the ML results and rationalize the viscosity predictions. The result is an ML approach that accurately predicts viscosity and will aid in accelerating the design of appropriate DESs for industrial applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Improving the Performance of Integrated Power-Hardware-in-the-Loop and Quasi-Static Time-Series Simulations

This paper introduces a novel multirate cosimulation architecture that overcomes previous challenges integrating faster real-time, microseconds-scale, power-hardware-in-the-loop (PHIL) with larger scale but slower, real-time, seconds-scale, quasi-static time-series (QSTS) simulation. Specifically, an intermediate reduced-equivalent electromagnetic transient (EMT) model duplicates the key power system topology to capture high-speed dynamics and converge hardware-power system interactions between QSTS updates. Exchanging multidimensional parameter vectors (loads, control status, etc.) between the QSTS and EMT models enables capturing the interactions of the rich, high-node-count (thousands of electrical nodes) QSTS model with the hardware. This architecture offers full spatial resolution from QSTS, including individual load dynamics, actual distribution management system controls, and nodal voltages. Simultaneously, the intermediate EMT model provides more detailed high-speed transients of key distributed energy resource hardware interactions. In addition, we use careful PHIL interface design and the exchange of complex power data, rather than current, to further improve performance. This cosimulation architecture is demonstrated by testing a simulated distribution system with an interconnected 500-kVA advanced photovoltaic (PV) inverter in PHIL. In conclusion, the architecture successfully captured PV local volt-volt ampere reactive (volt-VAR) control interactions with the larger network under time-varying electrical and weather conditions.

42 ENGINEERING↗

Vehicle Position Detection Based on Machine Learning Algorithms in Dynamic Wireless Charging

Dynamic wireless charging (DWC) has emerged as a viable approach to mitigate range anxiety by ensuring continuous and uninterrupted charging for electric vehicles in motion. DWC systems rely on the length of the transmitter, which can be categorized into long-track transmitters and segmented coil arrays. The segmented coil array, favored for its heightened efficiency and reduced electromagnetic interference, stands out as the preferred option. However, in such DWC systems, the need arises to detect the vehicle’s position, specifically to activate the transmitter coils aligned with the receiver pad and de-energize uncoupled transmitter coils. This paper introduces various machine learning algorithms for precise vehicle position determination, accommodating diverse ground clearances of electric vehicles and various speeds. Through testing eight different machine learning algorithms and comparing the results, the random forest algorithm emerged as superior, displaying the lowest error in predicting the actual position.

47 OTHER INSTRUMENTATION↗

Naturally light Dirac and pseudo-Dirac neutrinos from left-right symmetry

We develop a class of left-right symmetric theories based on the gauge group SU(3) c x SU(2) L x SU(2) R x U(1) with a generalized seesaw mechanism for generating the charged fermion masses. Neutrinos are naturally Dirac particles in this setup with their small masses arising from two-loop quantum corrections. We evaluate these two-loop diagrams exactly and analyze the flavor structure of the lepton sector. We find excellent fits to neutrino oscillation data, independent of the right-handed gauge symmetry breaking scale. We also explore the possibility that neutrinos are pseudo-Dirac particles in this framework, with the tiny mass splittings between active and sterile neutrinos arising from Planck-induced corrections and find possible realizations. These models can be tested in the near future with precision cosmological measurements of ΔN eff in CMB which is predicted to be ≃ 0.14. This class of models allows for a solution to the strong CP problem via parity symmetry without the need for an axion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

KiDS-1000 cosmology: Combined second- and third-order shear statistics

Aims.In this work, we perform the first cosmological parameter analysis of the fourth release of Kilo Degree Survey (KiDS-1000) data with second- and third-order shear statistics. This paper builds on a series of studies aimed at describing the roadmap to third-order shear statistics. Methods.We derived and tested a combined model of the second-order shear statistic, namely, the COSEBIs and the third-order aperture mass statistics 〈ℳ ap 3 〉 in a tomographic set-up. We validated our pipeline withN-body mock simulations of the KiDS-1000 data release. To model the second- and third-order statistics, we used the latest version of HMCODE2020 for the power spectrum and BIHALOFITfor the bispectrum. Furthermore, we used an analytic description to model intrinsic alignments and hydro-dynamical simulations to model the effect of baryonic feedback processes. Lastly, we decreased the dimension of the data vector significantly by considering only equal smoothing radii for the 〈ℳ ap 3 〉 part of the data vector. This makes it possible to carry out a data analysis of the KiDS-1000 data release using a combined analysis of COSEBIs and third-order shear statistics. Results.We first validated the accuracy of our modelling by analysing a noise-free mock data vector, assuming the KiDS-1000 error budget, finding a shift in the maximum of the posterior distribution of the matter density parameter, ΔΩ m < 0.02 σ Ω m , and of the structure growth parameter, ΔS 8 < 0.05 σ S 8 . Lastly, we performed the first KiDS-1000 cosmological analysis using a combined analysis of second- and third-order shear statistics, where we constrained Ω m = 0.248 −0.055 +0.062 andS 8 = σ 8 √(Ω m /0.3 )= 0.772 ± 0.022. The geometric average on the errors of Ω m andS 8 of the combined statistics decreases, compared to the second-order statistic, by a factor of 2.2.

Astronomy & Astrophysics↗

Stacked Low-Inertia Converter or Solid-State Transformer: Modeling and Model Predictive Priority-Shifting Control for Voltage Balance

This paper presents control challenges of stacked low-inertia converter (SLIC) or cascaded reduced dc-link solid-state transformer (SST) and proposes a novel model predictive priority-shifting (MPPS) control with implicit modulator and a discrete-time large-signal model for voltage balancing and dc-link regulation. Low-inertia converters, featuring small electrolytic capacitor-less dc links, dramatically reduce cost, size, and weight compared to conventional solutions. However, without a large dc-link buffer, the input and output are tightly coupled, leading to significant control challenges. The control becomes even more challenging with these converters stacked input-series output-parallel (ISOP) for medium-voltage (MV) grid, which causes coupling between the modules besides the coupling within each module. This paper analyzes the multi-objective, multi-degree of freedom control problem, using the modular soft-switching solid-state transformer (M-S4T) as an example of the SLIC. First, distribution of control efforts under controller saturation is critical because multiple control objectives can be conflicting, especially when the module voltages are unbalanced and are being restored. The MPPS can shift the priorities to address this issue. Second, due to the low inertia and high dc-link ripple, classic space vector pulse-width modulation (SVPWM), average model with small-ripple assumption, and control design based on small-signal model cannot accurately modulate, model, and control the nonlinear reduced dc link. Therefore, a discrete-time large-signal model of the M-S4T is established to derive the predictive control in the MPPS. The MPPS and the PI control are compared in MV simulations to show the issue of applying the PI to the SLIC and the effectiveness of the MPPS for voltage balancing and dc-link regulation in a deadbeat manner. Finally, the proposed control is tested on a 5 kV ISOP SiC SST prototype to verify priority shifting to address controller saturation issue and fast and robust voltage balancing.

14 SOLAR ENERGY↗

Calibration and Localization of Optically Pumped Magnetometers Using Electromagnetic Coils

In this paper, we propose a method to estimate the position, orientation, and gain of a magnetic field sensor using a set of (large) electromagnetic coils. We apply the method for calibrating an array of optically pumped magnetometers (OPMs) for magnetoencephalography (MEG). We first measure the magnetic fields of the coils at multiple known positions using a well-calibrated triaxial magnetometer, and model these discreetly sampled fields using vector spherical harmonics (VSH) functions. We then localize and calibrate an OPM by minimizing the sum of squared errors between the model signals and the OPM responses to the coil fields. We show that by using homogeneous and first-order gradient fields, the OPM sensor parameters (gain, position, and orientation) can be obtained from a set of linear equations with pseudo-inverses of two matrices. The currents that should be applied to the coils for approximating these low-order field components can be determined based on the VSH models. Computationally simple initial estimates of the OPM sensor parameters follow. As a first test of the method, we placed a fluxgate magnetometer at multiple positions and estimated the RMS position, orientation, and gain errors of the method to be 1.0 mm, 0.2°, and 0.8%, respectively. Lastly, we calibrated a 48-channel OPM array. The accuracy of the OPM calibration was tested by using the OPM array to localize magnetic dipoles in a phantom, which resulted in an average dipole position error of 3.3 mm. The results demonstrate the feasibility of using electromagnetic coils to calibrate and localize OPMs for MEG.

47 OTHER INSTRUMENTATION↗

Phylogenetic and ecological drivers of the avian lung mycobiome and its potentially pathogenic component

Vertebrate lungs contain diverse microbial communities, but little is known about the drivers of community composition or consequences for health. Microbiome assembly by processes such as dispersal, coevolution, and host-switching can be probed with comparative surveys; however, few studies exist for lung microbiomes, particularly for the fungal component, the mycobiome. Distinguishing among fungal taxa that are generalist or specialist symbionts, potential pathogens, or incidentally inhaled spores is urgent because of potential for emerging diseases. Here, we characterize the avian lung mycobiome and test the relative influences of environment, phylogeny, and functional traits. We used metabarcoding and culturing from 195 lung samples representing 32 bird species across 20 families. We identified 526 fungal taxa as estimated by distinct sequence types (zOTUs) including many opportunistic pathogens. These were predominantly from the phylum Ascomycota (79%) followed by Basidiomycota (16%) and Mucoromycota (5%). Yeast and yeast-like taxa (Malassezia, Filobasidium, Saccharomyces, Meyerozyma, and Aureobasidium) and filamentous fungi (Cladosporium, Alternaria, Neurospora, Fusarium, and Aspergillus) were abundant. Lung mycobiomes were strongly shaped by environmental exposure, and further modulated by host identity, traits, and phylogenetic affinities. Our results implicate migratory bird species as potential vectors for long-distance dispersal of opportunistically pathogenic fungi.

59 BASIC BIOLOGICAL SCIENCES↗

J-PLUS: Support vector regression to measure stellar parameters

Stellar parameters are among the most important characteristics in studies of stars which, in traditional methods, are based on atmosphere models. However, time, cost, and brightness limits restrain the efficiency of spectral observations. The Javalambre Photometric Local Universe Survey (J-PLUS) is an observational campaign that aims to obtain photometry in 12 bands. Owing to its characteristics, J-PLUS data have become a valuable resource for studies of stars. Machine learning provides powerful tools for efficiently analyzing large data sets, such as the one from J-PLUS, and enables us to expand the research domain to stellar parameters. The main goal of this study is to construct a support vector regression (SVR) algorithm to estimate stellar parameters of the stars in the first data release of the J-PLUS observational campaign. The training data for the parameter's regressions are featured with 12-waveband photometry from J-PLUS and are crossidentified with spectrum-based catalogs. These catalogs are from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope, the Apache Point Observatory Galactic Evolution Experiment, and the Sloan Extension for Galactic Understanding and Exploration. We then label them with the stellar effective temperature, the surface gravity, and the metallicity. Ten percent of the sample is held out to apply a blind test. We develop a new method, a multi-model approach, in order to fully take into account, the uncertainties of both the magnitudes and the stellar parameters. The method utilizes more than 200 models to apply the uncertainty analysis. We present a catalog of 2 493 424 stars with the root mean square error of 160 K in the effective temperature regression, 0.35 in the surface gravity regression, and 0.25 in the metallicity regression. We also discuss the advantages of this multi-model approach and compare it to other machine-learning methods.

79 ASTRONOMY AND ASTROPHYSICS↗

Cyber-Attack Detection for Photovoltaic Farms Based on Power-Electronics-Enabled Harmonic State Space Modeling

Here in this paper, a physics-data-based detection method is proposed to detect a variety of cyber-attacks in Photovoltaic (PV) farms using the power electronics-enabled harmonic state space (HSS) models, which, to our knowledge, is original. At the device level, HSS-based detection is developed to monitor harmonic vectors of individual PV converter with minimum sensor measurements, thus improving accuracy and robustness compared to Kalman Filter-based detection. At the system level that involves multiple PV converters, a clustering approach is developed to investigate attack propagation and accurately locate attack sources within a PV farm. The proposed approach is one of the first attempts to address PV security through interaction between the device and system, maximizing the accuracy and robustness at different levels. To verify the feasibility, a comprehensive attacks model is built, including single attack, coordinated attacks, and replay attacks. Besides, the impacts of irradiance changes are taken into consideration in the test scenarios. With the real-time data acquisition and hardware-in-the-loop testbed, comprehensive test results are provided to verify the feasibility of the proposed detection methodology.

42 ENGINEERING↗

Using the Endogenous CRISPR-Cas System of Heliobacterium modesticaldum To Delete the Photochemical Reaction Center Core Subunit Gene

In Heliobacterium modesticaldum , as in many Firmicutes , deleting genes by homologous recombination using standard techniques has been extremely difficult. The cells tend to integrate the introduced plasmid into the chromosome by a single recombination event rather than perform the double recombination required to replace the targeted locus. Transformation with a vector containing only a homologous recombination template for replacement of the photochemical reaction center gene pshA produced colonies with multiple genotypes, rather than a clean gene replacement. To address this issue, we required an additional means of selection to force a clean gene replacement. In this study, we report the genetic structure of the type I-A and I-E CRISPR-Cas systems from H. modesticaldum , as well as methods to leverage the type I-A system for genome editing. In silico analysis of the CRISPR spacers revealed a potential consensus protospacer adjacent motif (PAM) required for Cas3 recognition, which was then tested using an in vivo interference assay. Introduction of a homologous recombination plasmid that carried a miniature CRISPR array targeting sequences in pshA (downstream of a naturally occurring PAM sequence) produced nonphototrophic transformants with clean replacements of the pshA gene with ~80% efficiency. Mutants were confirmed by PCR, sequencing, optical spectroscopy, and growth characteristics. This methodology should be applicable to any genetic locus in the H. modesticaldum genome. The heliobacteria are the only phototrophic members of the largely Gram-positive phylum Firmicutes , which contains medically and industrially important members, such as Clostridium difficile and Clostridium acetobutylicum . Heliobacteria are of interest in the study of photosynthesis because their photosynthetic system is unique and the simplest known. Since their discovery in the early 1980s, work on the heliobacteria has been hindered by the lack of a genetic transformation system. The problem of introducing foreign DNA into these bacteria has been recently rectified by our group; however, issues still remained for efficient genome editing. The significance of this work is that we have characterized the endogenous type I CRISPR-Cas system in the heliobacteria and leveraged it to assist in genome editing. Using the CRISPR-Cas system allowed us to isolate transformants with precise replacement of the pshA gene encoding the main subunit of the photochemical reaction center.

59 BASIC BIOLOGICAL SCIENCES↗

Huge ensembles – Part 1: Design of ensemble weather forecasts using spherical Fourier neural operators

Abstract. Simulating low-likelihood high-impact extreme weather events in a warming world is a significant and challenging task for current ensemble forecasting systems. While these systems presently use up to 100 members, larger ensembles could enrich the sampling of internal variability. They may capture the long tails associated with climate hazards better than traditional ensemble sizes. Due to computational constraints, it is infeasible to generate huge ensembles (comprised of 1000–10 000 members) with traditional, physics-based numerical models. In this two-part paper, we replace traditional numerical simulations with machine learning (ML) to generate hindcasts of huge ensembles. In Part 1, we construct an ensemble weather forecasting system based on spherical Fourier neural operators (SFNOs), and we discuss important design decisions for constructing such an ensemble. The ensemble represents model uncertainty through perturbed-parameter techniques, and it represents initial condition uncertainty through bred vectors, which sample the fastest-growing modes of the forecast. Using the European Centre for Medium-Range Weather Forecasts Integrated Forecasting System (IFS) as a baseline, we develop an evaluation pipeline composed of mean, spectral, and extreme diagnostics. With large-scale, distributed SFNOs with 1.1 billion learned parameters, we achieve calibrated probabilistic forecasts. As the trajectories of the individual members diverge, the ML ensemble mean spectra degrade with lead time, consistent with physical expectations. However, the individual ensemble members' spectra stay constant with lead time. Therefore, these members simulate realistic weather states during the rollout, and the ML ensemble passes a crucial spectral test in the literature. The IFS and ML ensembles have similar extreme forecast indices, and we show that the ML extreme weather forecasts are reliable and discriminating. These diagnostics ensure that the ensemble can reliably simulate the time evolution of the atmosphere, including low-likelihood high-impact extremes. In Part 2, we generate a huge ensemble initialized each day in summer 2023, and we characterize the simulations of extremes.

Mahesh, Ankur↗