Search NASA⌕ Search

SEARCH · Search NASA

Results for “splitting methods”

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.

1,777 records · Page 99
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Extensions to models of the galaxy–halo connection

ABSTRACT We explore two widely used empirical models for the galaxy–halo connection, subhalo abundance matching (SHAM) and the halo occupation distribution (HOD), and compare them with the hydrodynamical simulation IllustrisTNG (TNG) for multiple statistics quantifying the galaxy distribution at $n_{\rm gal}\approx 1.3\times 10^{-3}\, ({\rm Mpc}\,h^{-1})^{-3}$. We observe that in their most straightforward implementations, both models fail to reproduce the two-point clustering measured in TNG. We find that SHAM models that use the relaxation velocity, Vrelax, and the peak velocity, Vpeak, perform best, and match the clustering reasonably well, although neither captures adequately the one-halo clustering. Splitting the total sample into sub-populations, we discover that SHAM overpredicts the clustering of high-mass, blue, star-forming, and late-forming galaxies and underpredicts that of low-mass, red, quiescent, and early-forming galaxies. We also study various baryonic effects, finding that subhaloes in the dark-matter-only simulation have consistently higher values of their SHAM-proxy properties than their full-physics counterparts. We then consider a 2D implementation of the HOD model augmented with a secondary parameter (environment, velocity anisotropy, σ2Rhalf-mass, and total potential) tuned so as to match the two-point clustering of the IllustrisTNG galaxies on large scales. We analyse these galaxy populations adopting alternative statistical tools such as galaxy–galaxy lensing, void–galaxy cross-correlations, and cumulants of the density field, finding that the hydrodynamical galaxy distribution disfavours σ2Rhalf-mass and the total potential as secondary parameters, while the environment and velocity anisotropy samples are consistent with full physics across all statistical probes examined. Our results demonstrate the power of examining multiple statistics for determining the secondary parameters that are vital for understanding the galaxy–halo connection.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The Splitting of Double-component Active Asteroid P/2016 J1 (PANSTARRS)

We present deep imaging observations, orbital dynamics, and dust-tail model analyses of the double-component asteroid P/2016 J1 (J1-A and J1-B). The observations were acquired at the Gran Telescopio Canarias (GTC) and the Canada–France–Hawaii Telescope (CFHT) from mid-March to late July of 2016. A statistical analysis of backward-in-time integrations of the orbits of a large sample of clone objects of P/2016 J1-A and J1-B shows that the minimum separation between them occurred most likely ∼2300 days prior to the current perihelion passage, i.e., during the previous orbit near perihelion. This closest approach was probably linked to a fragmentation event of their parent body. Monte Carlo dust-tail models show that those two components became active simultaneously ∼250 days before the current perihelion, with comparable maximum loss rates of ∼0.7 and ∼0.5 kg s{sup −1}, and total ejected masses of 8 × 10{sup 6} and 6 × 10{sup 6} kg for fragments J1-A and J1-B, respectively. Consequently, the fragmentation event and the present dust activity are unrelated. The simultaneous activation times of the two components and the fact that the activity lasted 6–9 months or longer, strongly indicate ice sublimation as the most likely mechanism involved in the dust emission process.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Anion exchange membrane test protocol validation

This study presents the validation of protocols for measuring ion exchange capacity (IEC) and alkaline stability of anion exchange membranes (AEMs) for low-temperature water electrolysis. While protocols are often tested within individual laboratories, their results across multiple laboratories with varying equipment, environmental conditions, and personnel qualification remain unverified. The validation involved Los Alamos National Laboratory (LANL), National Renewable Energy Laboratory (NREL), and University of Oregon (UO) using the same commercially available AEM to assess reproducibility and reliability of the protocols under diverse conditions. For the IEC protocol, results across laboratories were consistent within ±10% of the NMR-determined reference value. The alkaline stability protocol could pose greater challenges due to factors such as variations in sample collection timing, preservation methods, and analytical techniques, but consistent test results for percentage IEC loss were demonstrated across institutions. These results highlight the reliability and applicability of the protocols, emphasizing the importance of validation to ensure consistency in diverse research environments.

08 HYDROGEN↗
DOE OSTI2016

Ag@ZnO core-shell nanoparticles study by first principle: The structural, magnetic and optical properties

Ag@ZnO core-shell nanoparticles of around 72 atoms have been investigated by the density functional theory, revealing proving for the first time that the core-shell structure exhibits a shrinkage phenomenon from outer shell in agreement with the other studies in literatures. Our calculations predict that the Ag@ZnO core-shell structure is a ferromagnetic spin polarized state, and the magnetism mainly stems from the spin splitting of 2p electrons of O atoms. In addition, the total and partial DOS of Ag@ZnO indicate that the nanostructure is a half-metallic nanoparticle and has the characters of the p-type semiconductor. Furthermore, the optical properties calculations show that the absorption edge of Ag@ZnO have a red shift and good photocatalysis compare to that of the bulk ZnO. These results of the Ag@ZnO core-shell structure obtain a well agreement with the experimental measurement. - Graphical abstract: Geometric structure of (a) Ag@ZnO core-shell nanostructure; (b) the core of Ag; (c) the shell of ZnO The core-shell nanoparticle Ag@ZnO contains Ag inner core of radius of 4 Å and ZnO outer shell with thickness of 2 Å. Ag@ZnO core-shell nanoparticles of around 72 atoms have been proved for the first time that the core-shell structure exhibit a shrinkage phenomenon from outer shell. Our calculations predict that the Ag@ZnO core-shell structure is a half-metallic nanoparticle and has the characters of the p-type semiconductor. The absorption edge of Ag@ZnO have a red shift and get good photo-catalysis compare to that of the bulk ZnO.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Tailoring the Surface of Metal Halide Perovskites to Enable the Atomic Layer Deposition of Metal Oxide Contacts

Replacing organic contact layers with inorganic counterparts, such as metal oxides, is one strategy for improving long-term device stability in metal halide perovskite solar cells. Often, the methods used to deposit metal oxide thin films are incompatible with metal halide perovskites, creating challenges for the fabrication of contacts above the perovskite absorber layer. In this study, we utilize a one-step, solution treatment of the top surface of Cs 0.25 FA 0.75 Pb(Br 0.20 I 0.80 ) 3 to create a thin (~1 nm) overlayer of lead sulfide (PbS) to protect the underlying perovskite during subsequent deposition. X-ray characterization of the surface region shows that the PbS overlayer limits undesirable changes to the perovskite structure and stoichiometry during atomic layer deposition (ALD) of SnO 2 . This protection enables ALD growth of SnO 2 electron contacts on top of the perovskite without an organic transport layer (e.g., C 60 ), resulting in a solar cell with a power conversion efficiency of 5.8%. This result is a marked improvement over devices with ALD SnO 2 grown directly on the perovskite without a PbS overlayer, which produce no power output. Furthermore, the interface characterization and device results in this study highlight some of the key challenges associated with ALD metal oxide growth on perovskite materials and can help inform the future design of inorganic contact layer deposition in solar photovoltaic technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Prediction of s ± -Wave Superconductivity Enhanced by Electronic Doping in Trilayer Nickelates La 4 Ni 3 O 10 under Pressure

Motivated by the recently reported signatures of superconductivity in trilayer La 4 ⁢Ni 3 ⁢O 10 under pressure, here we comprehensively study this system using ab initio and random-phase approximation techniques. Without electronic interactions, the Ni d 3z 2 –r 2 orbitals show a bonding-antibonding and nonbonding splitting behavior via the O p z orbitals inducing a “trimer” lattice in La 4⁢ Ni 3 ⁢O 10 , analogous to the dimers of La 3 ⁢Ni 2 ⁢O 7 . The Fermi surface consists of three electron sheets with mixed e g orbitals, and a hole and an electron pocket made up of the d 3⁢z 2 –r 2 orbital, suggesting a Ni two-orbital minimum model. In addition, we find that superconducting pairing is induced in the s ± -wave channel due to partial nesting between the M = (π,π) centered pockets and portions of the Fermi surface centered at the Γ = (0,0) point. With changing electronic density n, the s ± instability remains leading and its pairing strength shows a domelike behavior with a maximum around n = 4.2 ( ~6.7% electron doping). The superconducting instability disappears at the same electronic density as that in the new 1313 stacking La 3 ⁢Ni 2 ⁢O 7 , correlated with the vanishing of the hole pocket that arises from the trilayer sublattice, suggesting that the high-T c superconductivity of La 3⁢ Ni 2 ⁢O 7 does not originate from a trilayer and monolayer structure. Furthermore, we confirm the experimentally proposed spin state in La 4 ⁢Ni 3 ⁢O 10 with an in-plane (π, π) order and antiferromagnetic coupling between the top and bottom Ni layers, and spin zero in the middle layer.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Dependence of interface energetics and kinetics on catalyst loading in a photoelectrochemical system

Solar hydrogen production by the photoelectrochemical method promises a means to store solar energy. While it is generally understood that the process is highly sensitive to the nature of the interface between the semiconductor and the electrolyte, a detailed understanding of this interface is still missing. For instance, few prior studies have established a clear relationship between the interface energetics and the catalyst loading amount. Here we aim to study this relationship on a prototypical Si-based photoelectrochemical system. Two types of interfaces were examined, one with GaN nanowires as a protection layer and one without. It was found that when GaN was present, higher Pt loading (> 0.1 μg/cm 2 ) led to not only better water reduction (and, hence, hydrogen evolution) kinetics but also more favorable interface energetics for greater photovoltages. In the absence of the protection layer, by stark contrast, increased Pt loading exhibited no measurable influence on the interface energetics, and the main difference was observed only in the hydrogen evolution kinetics. Furthermore, the study sheds new light on the importance of interface engineering for further improvement of photoelectrochemical systems, especially concerning the role of catalysts and protection layers.

14 SOLAR ENERGY↗
DOE OSTI2016

A Methodology for Loading the Advanced Test Reactor Driver Core for Experiment Analysis

The Advanced Test Reactor (ATR) at the Idaho National Laboratory (INL) is one of the premier irradiation facilities in the world. The ATR hosts experiments from a variety of laboratories, companies, and disciplines. Because of its flexibility in providing customized irradiation conditions to experiments, such as coolant chemistry, temperature, and reactor power, analysis for experiments, both from a programmatic and safety perspective, can be a challenging endeavor. Because of these experiment variables, there is always an effort to improve upon the existing analysis methodology in order to provide the sponsor of any given experiment with the most accurate predictions possible. The paper discusses the effort to improve upon the methodology for loading the ATR driver core for experiment analysis. Experiment analysis at INL is generally broken up into three separate disciplines: neutronics, thermal-hydraulics, and structural. The focus of this paper is neutronics analysis. After the experiment has a preliminary design, the neutronics analyst will create a model in the Monte Carlo N-Particle (MCNP) transport code to calculate fuel and material heating rates from MCNP tallies and provide those to the thermal-hydraulics analyst, and calculate reactivity worths, perturbation on the axial fission profile of the ATR driver core by the experiment, void worths, and experiment backup worths. Also, there are usually several programmatic parameters that the experiment sponsor requests that the neutronics analyst calculates. These can include heat generation rates, burnup and displacements per atom (DPA). In order to accurately analyze these parameters, an important input into the physics analyst's MCNP model is the loading of the ATR driver core. The ATR driver core consists of 40, 19-plate highly-enriched uranium (HEU) fuel elements, arranged in a serpentine shape, as pictured in Figure 1. The ATR driver core is responsible for the vast majority of the thermal power of the reactor and the nuclear interactions with the experiments being irradiated in ATR. A unique feature of ATR is the 'power tilt' across the reactor. A power tilt means that the Northwest, Northeast, Southwest, and Southeast Lobes of ATR frequently operate at different powers, providing different irradiation conditions for the experiments in and around the lobes of the reactor. In order to achieve lobe power splits, different fuel elements for each of the 40 positions are chosen by ATR Reactor Engineering so that lobe power splits can be maintained during an operating cycle. A typical fuel loading for an ATR cycle consists of 14 to 18 new fuel elements and 2 to 26 depleted fuel elements. These loadings are typically chosen from three different types of fuel elements. All three ATR fuel element types share a common design, spelled out in Chapter 4 of the ATR Safety Analysis Report (SAR). The fuel element meat is HEU UAlx dispersion fuel sandwiched between aluminum cladding with a length of 121.92 cm (4 ft), with 19 plates per fuel element. The three variations of this common design are designated by the first two letters in their serial numbers: XA, YA, and NB. XA elements are most commonly used. These elements have nominally 1075 g of uranium-235 and 0.66 g of boron-10. The NB elements have the same uranium-235 loading as the XA elements, but no boron in any plates. These elements are used when excess reactivity is needed. The YA elements have the same loading as the XA elements, except that the outermost plate (Plate 19) has no boron or uranium. These elements are used to mitigate aged reflector damage. In order to accurately model the ATR driver core, all three types of fuel elements must be accounted for in the MCNP model. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

Recent Measurement of Flavor Asymmetry of Antiquarks in the Proton by Drell–Yan Experiment SeaQuest at Fermilab

A measurement of the flavor asymmetry of the antiquarks ($\bar{d}$ and $\bar{u}$) in the proton is described in this thesis. The proton consists of three valence quarks, sea quarks, and gluons. Antiquarks in the proton are sea quarks. They are generated from the gluon splitting: g → q + $\bar{q}$. According to QCD (Quantum Chromodynamics), the gluon splitting is independent of quark flavor. It suggests that the amounts of $\bar{d}$ and $\bar{u}$ should be the same in the proton. However, the NMC experiment at CERN found that the amount of $\bar{d}$ is larger than that of $\bar{u}$ in the proton using the deep inelastic scattering in 1991. This result is obtained for $\bar{d}$ and $\bar{u}$ integrated over Bjorken x. Bjorken x is the fraction of the momentum of the parton to that of the proton. The NA51 experiment (x ~ 0.2) at CERN and E866/NuSea experiment (0.015 < x < 0.35) at Fermilab measured the flavor asymmetry of the antiquarks ($\bar{d}$/$\bar{u}$) in the proton as a function of x using Drell–Yan process. The experiments reported that the flavor symmetry is broken over all measured x values. Understanding the flavor asymmetry of the antiquarks in the proton is a challenge of the QCD. The theo- retical investigation from the first principle of QCD such as lattice QCD calculation is important. In addition, the QCD effective models and hadron models such as the meson cloud model can also be tested with the flavor asymmetry of antiquarks. From the experimental side, it is important to measure with higher accuracy and in a wider x range. The SeaQuest (E906) experiment measures $\bar{d}$/$\bar{u}$ at large x (0.15 < x < 0.45) accurately to understand its behavior. The SeaQuest experiment is a Drell–Yan experiment at Fermi National Accelerator Laboratory (Fermilab). In the Drell–Yan process of proton-proton reaction, an antiquark in a proton and a quark in another proton annihilate and create a virtual photon, which then decays into a muon pair (q$\bar{q}$ → γ* → µ + µ - ). The SeaQuest experiment uses a 120 GeV proton beam extracted from Fermilab’s Main Injector. The proton beam interacts with hydrogen and deuterium targets. The SeaQuest spectrometer detects the muon pairs from the Drell–Yan process. The $\bar{d}$/$\bar{u}$ ratio at 0.1 < x < 0.58 is extracted from the number of detected Drell–Yan muon pairs. After the detector construction, commissioning run and detector upgrade, the SeaQuest experiment started the physics data acquisition from 2013. We finished so far three periods of physics data acquisition. The fourth period is in progress. The detector construction, detector performance evaluation, data taking and data analysis for the flavor asymmetry of the antiquarks $\bar{d}$/$\bar{u}$ in the proton are my contribution to SeaQuest. The cross section ratio of Drell–Yan process in p- p and p-d reactions is obtained from dimuon yields. In the experiment with high beam intensity, it is important to control the tracking efficiency of charged particles through the magnetic spectrometer. The tracking efficiency depends on the chamber occupancy, and the appropriate method for the correction is important. The chamber occupancy is the number of hits in drift chambers. A new method of the correction for the tracking efficiency is developed based on the occupancy, and applied to the data. This method reflects the real response of the drift chambers. Therefore, the systematic error is well controlled by this method. The flavor asymmetry of antiquarks is obtained at 0.1 < x < 0.58. At 0.1 < x < 0.45, the result is $\bar{d}$/$\bar{u}$ > 1. The result at 0.1 < x < 0.24 agrees with the E866 result. The result at x > 0.24, however, disagrees with the E866 result. The result at 0.45 < x < 0 the statistical errors. u¯ results extracted from experiments are used to investigate the validity of the theoretical models. The present experimental result provides the data points in wide x region. It is useful for understanding the proton structure in the light of QCD and effective hadron models. The present result has a practical application as well. Antiquark distributions are important as inputs to simulations of hadron reactions such as W± production in various experiments. The new knowledge on antiquark distributions helps to improve the precision of the simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Structural and photoelectrochemical properties in the thin film system Cu–Fe–V–O and its ternary subsystems Fe–V–O and Cu–V–O

Thin-film material libraries in the ternary and quaternary metal oxide systems Fe–V–O, Cu–V–O, and Cu–Fe–V–O were synthesized using combinatorial reactive co-sputtering with subsequent annealing in air. Their compositional, structural, and functional properties were assessed using high-throughput characterization methods. Prior to the investigation of the quaternary system Cu–Fe–V–O, the compositions (Fe 61 V 39 )O x and (Cu 52 V 48 )O x with promising photoactivity were identified from their ternary subsystems Fe–V–O and Cu–V–O, respectively. Two Cu–Fe–V–O material libraries with (Cu 29-72 Fe 4-27 V 22-57 )O x and (Cu 11-55 Fe 27-73 V 12-34 )O x composition spread were investigated. Seven mixed ternary and quaternary phase regions were identified: I (α-Cu 3 FeV 6 O 26 /FeVO 4 ), II (Cu 5 V 2 O 10 /FeVO 4 /α-Cu 3 Fe 4 V 6 O 26 ), III (Cu 5 V 2 O 10 ), IV (Cu 5 V 2 O 10 /FeVO 4 , V (FeVO 4 /γ-Cu 2 V 2 O 7 /α-Cu 3 Fe 4 V 6 O 26 ), VI (β-Cu 2 V 2 O 7 /α-Cu 3 Fe 4 V 6 O 26 /FeVO 4 ), and VII (β-Cu 3 Fe 4 V 6 O 26 /FeVO 4 ). Furthermore, in the investigated composition range, two photoactive regions, (Cu 53 Fe 7 V 40 )O x and (Cu 45 Fe 21 V 34 )O x , were identified, exhibiting 103 μ A/cm 2 and 108 μ A/cm 2 photocurrent density for the oxygen evolution reaction at 1.63 V vs reversible hydrogen electrode, respectively. The highest photoactive region (Cu 45 Fe 21 V 34 )O x comprises the dominant α-Cu 3 Fe 4 V 6 O 24 phase and minor FeVO 4 phase. This photoactive region corresponds to having an indirect bandgap of 1.87 eV and a direct bandgap of 2.58 eV with an incident photon-to-current efficiency of 30% at a wavelength of 310 nm.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Small Scale Field Test Demonstrating CO 2 Sequestration In Arbuckle Saline Aquifer And By CO 2 -Eor At Wellington Field, Sumner County, Kansas

The objectives of this project are to understand the processes that occur when a maximum of 70,000 metric tonnes of CO 2 are injected into two different formations to evaluate the response in different lithofacies and depositional environments. The evaluation will be accomplished through the use of both in situ and indirect MVA (monitoring, verification, and accounting) technologies. The project will optimize for carbon storage accounting for 99% of the CO 2 using lab and field testing and comprehensive characterization and modeling techniques. Site characterization and CO 2 injection should demonstrate state-of-the-art MVA tools and techniques to monitor and visualize the injected CO 2 plume and to refine geomodels developed using nearly continuous core, exhaustive wireline logs, and well tests and a multi-component 3-D seismic survey. Reservoir simulation studies will map the injected CO 2 plume and estimate tonnage of CO 2 stored in solution, as residual gas, and by mineralization and integrate MVA results and reservoir models shall be used to evaluate CO 2 leakage. A rapid-response mitigation plan was developed to minimize CO 2 leakage and provide a comprehensive risk management strategy. The CO 2 was intended to be supplied from a reliable facility and have an adequate delivery and quality of CO 2 . However, several unforeseen circumstances complicated this plan: (1) the initially negotiated CO 2 supply facility went offline and contracts associated with CO 2 supply had to be renegotiated, (2) a UIC Class VI permit proved to be difficult to obtain due to the experimental nature of the project. Both subjects are detailed in separate deliverables attached to this report. The CO 2 enhanced oil recovery (EOR) and geologic storage in Mississippian carbonate reservoir was sucessully deployed. Approximately 20,000 metric tons of CO 2 was injected in the upper part of the Mississippian reservoir to verify CO 2 EOR viability in carbonate reservoirs and evaluate a potential of transitioning to geologic CO 2 storage through EOR. A total of 1,101 truckloads, 19,803 metric tons—an average of 120 tonnes per day—were delivered over the course of injection that lasted from January 9 to June 21, 2016. After cessation of CO 2 injection, the KGS 2-32 well was converted to water injector and continues to operate. CO 2 EOR progression in the field was monitored weekly with fluid level, temperature, and production recording and formation fluid composition sampling. It is important to note that normally, CO 2 EOR pilots are less efficient than commercial operations due to lack of directional and precise well control, lack of surface facilities for CO 2 recycling, and other factors. As a result of this pilot CO 2 injection, the observed incremental average oil production increase was ~68% with only ~18% of injected CO 2 produced back. Decline curve analysis forecasts of additional cumulative oil produced were 32.44M STB to the end of 2027. Wellington Mississippian pilot efficiency by the end of forecast calculations is 11 MCF per barrel of produced oil. Using 32M STB oil production and $1,964,063 cost of CO 2 , CO 2 EOR cost per barrel of oil production is ~$60. Simple but robust monitoring technologies proved to be very efficient in detecting and locating CO 2 . High CO 2 reservoir retentions with low yields within an actively producing field could help to estimate real-world risks of CO 2 geological storage for future projects. The Wellington Field CO 2 EOR was executed in a controlled environment with high efficiency. This case study proves that CO 2 EOR could be successfully applied in Kansas carbonate reservoirs if CO 2 sources and associated infrastructure are available. Recent developments in unconventional resources development in Mid-Continent USA and associated large volume disposal of backflow water and the resulting seismic activity have brought more focus and attention to the Arbuckle Group in southern Kansas. Despite the commercial interest, limited essential information about reservoir properties and structural elements has impeded the management and regulation of disposal, an issue brought to the forefront by recent seismicity in and near areas of large volumes and rates of brine disposal. The Kansas Geological Survey (KGS) collected, compiled, and analyzed available data, including well logs, core data, step rate tests, drill stem tests, 2-D and 3-D seismic data, water level measurements, and others types of data. Several exploratory wells were drilled and core was collected and modern suites of logs were analyzed. Reservoir properties were populated into several site-specific geological models. The geological models illustrate the highly heterogeneous nature of the Arbuckle Group. Vertical and horizontal variability results in several distinct hydro-stratigraphic units that are the result of both depositional and diagenetic processes. During the course of this project, it has been demonstrated that advanced seismic interpretation methods can be used successfully for characterization of the Mississippian reservoir and Arbuckle saline aquifer. Analysis of post-stack 3-D seismic data at the Mississippian reservoir showed the response of a gradational velocity transition. Pre-stack gather analysis showed that porosity zones of the Mississippian and Arbuckle reservoirs exhibit characteristic amplitude versus offset (AVO) response. Simultaneous AVO inversion estimated P- and S-impedances. The 3-D survey gather azimuthal anisotropy analysis (AVAZ) provided information about the fault and fracture network and showed good agreement to the regional stress field and well data. Mississippian reservoir porosity and fracture predictions agreed well with the observed mobility of injected CO 2 in KGS well 2-32. Fluid substitution modeling predicted acoustic impedance reduction in the Mississippian carbonate reservoir introduced by the presence of CO 2 . Seismicity in the United States midcontinent has increased by orders of magnitude over the past decade. Spatiotemporal correlations of seismicity to wastewater injection operations have suggested that injection-related pore fluid pressure increases are inducing the earthquakes. In this investigation, we examine earthquake occurrence in southern Kansas and northern Oklahoma and its relation to the change in pore pressure. The main source of data comes from the Wellington Array in the Wellington oil field, in Sumner County, Kansas, which has monitored for earthquakes in central Sumner County, Kansas, since early 2015. The seismometer array was established to monitor CO 2 injection operations at Wellington Field. Although no seismicity was detected in association with the spring 2016 Mississippian CO 2 injection, the array has recorded more than 2,500 earthquakes in the region and is providing valuable understanding to induced seismicity. A catalog of earthquakes was built from this data and was analyzed for spatial and temporal changes, stress information, and anisotropy information. The region of seismic concern has been shown to be expanding through use of the Wellington earthquake catalog, which has revealed a northward progression of earthquake activity reaching the metropolitan area of Wichita. The stress orientation was also calculated from this earthquake catalog through focal mechanism inversion. The calculated stress orientation was confirmed through comparison to other stress measurements from well data and previous earthquake studies in the region. With this knowledge of the stress orientation, the anisotropy in the basement could be understood. This allowed for the anisotropy measurements to be correlated to pore pressure increases. The increase in pore pressure was monitored through time-lapse shear-wave anisotropy analysis. Since the onset of the observation period in 2010, the orientation of the fast shear wave has rotated 90°, indicating a change associated with critical pore pressure build up. The time delay between fast and slow shear wave arrivals has increased, indicating a corresponding increase in anisotropy induced by pore pressure rise. In-situ near-basement fluid pressure measurements corroborate the continuous pore pressure increase revealed by the shear-wave anisotropy analysis over the earthquake monitoring period. This research is the first to identify a change in pore fluid pressure in the basement using seismological data and it was recently published in the AAAS journal Science Advances (Nolte et al., 2017). The shear-wave splitting analysis is a novel application of the technique, which can be used in other regions to identify an increase in pore pressure. This increasing pore fluid pressure has become more regionally extensive as earthquakes are occurring in southern Kansas, where they previously were absent. These monitoring techniques and analyses provide new insight into mitigating induced seismicity’s impact on society.

02 PETROLEUM↗
DOE OSTI2020

Planck intermediate results

We present the NPIPE processing pipeline, which produces calibrated frequency maps in temperature and polarization from data from the Planck Low Frequency Instrument (LFI) and High Frequency Instrument (HFI) using high-performance computers. NPIPE represents a natural evolution of previous Planck analysis efforts, and combines some of the most powerful features of the separate LFI and HFI analysis pipelines. For example, following the LFI 2018 processing procedure, NPIPE uses foreground polarization priors during the calibration stage in order to break scanning-induced degeneracies. Similarly, NPIPE employs the HFI 2018 time-domain processing methodology to correct for bandpass mismatch at all frequencies. In addition, NPIPE introduces several improvements, including, but not limited to: inclusion of the 8% of data collected during repointing manoeuvres; smoothing of the LFI reference load data streams; in-flight estimation of detector polarization parameters; and construction of maximally independent detector-set split maps. For component-separation purposes, important improvements include: maps that retain the CMB Solar dipole, allowing for high-precision relative calibration in higher-level analyses; well-defined single-detector maps, allowing for robust CO extraction; and HFI temperature maps between 217 and 857 GHz that are binned into 0'.9 pixels (Nside = 4096), ensuring that the full angular information in the data is represented in the maps even at the highest Planck resolutions. The net effect of these improvements is lower levels of noise and systematics in both frequency and component maps at essentially all angular scales, as well as notably improved internal consistency between the various frequency channels. Based on the NPIPE maps, we present the first estimate of the Solar dipole determined through component separation across all nine Planck frequencies. The amplitude is (3366.6 ± 2.7) μK, consistent with, albeit slightly higher than, earlier estimates. From the large-scale polarization data, we derive an updated estimate of the optical depth of reionization of τ = 0.051 ± 0.006, which appears robust with respect to data and sky cuts. There are 600 complete signal, noise and systematics simulations of the full-frequency and detector-set maps. As a Planck first, these simulations include full time-domain processing of the beam-convolved CMB anisotropies. The release of NPIPE maps and simulations is accompanied with a complete suite of raw and processed time-ordered data and the software, scripts, auxiliary data, and parameter files needed to improve further on the analysis and to run matching simulations.

79 ASTRONOMY AND ASTROPHYSICS↗