Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Agreements”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Comparing Experimental Phase Behavior of Ion-Doped Block Copolymers with Theoretical Predictions Based on Selective Ion Solvation

The effects of salt-doping on the morphological behavior of block copolymers are well established but remain poorly understood, partially because of the challenge of resolving electrostatics in a heterogeneous medium with low average permittivity. By employing a recently developed field theory, we analyze the phase behavior of polystyrene-b-poly(ethylene oxide) (SEO) copolymers doped with lithium bis(trifluoromethanesulfonyl)imide salts (LiTFSI). Using a single fitting parameter, the ionic solvation radius, we obtain qualitative agreement between our theory and experimental data over a range of polymer molecular weights and copolymer compositions. Such agreement supports and highlights the need of solvation free energy to accurately describe the self-assembly of ion-doped block copolymers and demonstrates that experimentally observed dependence on molecular weight, not present in neutral block copolymers, can be rationalized by solvation effects. Overall, morphological variations are stronger than those predicted by the leading, linear order theory but can be captured by the full model.

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

Solvent Isotherms and Structural Transitions in Nanoparticle Superlattice Assembly

Here, we introduce a Molecular Theory for Compressible Fluids (MOLT-CF) that enables us to compute free energies and other thermodynamic functions for nanoparticle superlattices with any solvent content, including the dry limit. Quantitative agreement is observed between MOLT-CF and united-atom molecular dynamics simulations performed to assess the reliability and precision of the theory. Among other predictions, MOLT-CF shows that the amount of solvent within the superlattice decreases approximately linearly with its vapor pressure and that in the late stages of drying, solvent-filled voids form at lattice interstitials. Applied to single-component superlattices, MOLT-CF predicts fcc-to-bcc Bain transitions for decreasing vapor pressure and for increasing ligand length, both in agreement with experimental results. We explore the stability of other single-component phases and show that the C14 Frank-Kasper phase, which has been reported in experiments, is not a global free-energy minimum. Implications for precise assembly and prediction of multicomponent nanoparticle systems are discussed.

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

Surface Oxygen Depletion of Layered Transition Metal Oxides in Li-Ion Batteries Studied by Operando Ambient Pressure X-ray Photoelectron Spectroscopy

A new operando spectro-electrochemical setup was developed to study oxygen depletion from the surface of layered transition metal oxide particles at high degrees of delithiation. An NCM111 working electrode was paired with a chemically delithiated LiFePO 4 counter electrode in a fuel cell-inspired membrane electrode assembly (MEA). A propylene carbonate-soaked Li-ion conducting ionomer served as an electrolyte, providing both good electrochemical performance and direct probing of the NCM111 particles during cycling by ambient pressure X-ray photoelectron spectroscopy. The irreversible emergence of an oxygen-depleted phase in the O 1s spectra of the layered oxide particles was observed upon the first delithiation to high state-of-charge, which is in excellent agreement with oxygen release analysis via mass spectrometry analysis of such MEAs. By comparing the metal oxide-based O 1s spectral features to the Ni 2p 3/2 intensity, we can calculate the transition metal-to-oxygen ratio of the metal oxide close to the particle surface, which shows good agreement with the formation of a spinel-like stoichiometry as an oxygen-depleted phase. This new setup enables a deeper understanding of interfacial changes of layered oxide-based cathode active materials for Li-ion batteries upon cycling.

25 ENERGY STORAGE↗
DOE OSTI2023

Quantitative Simulations of Siloxane Adsorption in Metal–Organic Frameworks

We present a transferable force field (FF) for simulating the bulk properties of linear and cyclic siloxanes and the adsorption of these species in metal–organic frameworks (MOFs). Unlike previous FFs for siloxanes, our FF accurately reproduces the vapor–liquid equilibria of each species in the bulk phase. The quality of our FF combined with the Universal Force Field using standard Lorentz–Berthelot combining rules for MOF atoms was assessed in a wide range of MOFs without open metal sites, showing good agreement with dispersion-corrected density functional theory calculations. Predictions with this FF show good agreement with the limited experimental data for siloxane adsorption in MOFs that is available. As an example of using the FF to predict adsorption properties in MOFs, we present simulations examining entropy effects in binary linear and cyclic siloxane mixture coadsorption in the large-pore MOF with structure code FOTNIN.

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

Ionic Associations and Hydration in the Electrical Double Layer of Water-in-Salt Electrolytes

Water-in-Salt-Electrolytes (WiSEs) are an exciting class of concentrated electrolytes finding applications in energy storage devices because of their expanded electrochemical stability window, good conductivity and cation transference number, and fire-extinguishing properties. These distinct properties are thought to originate from the presence of an anion-dominated ionic network and interpenetrating water channels for cation transport, which indicates that associations in WiSEs are crucial to understanding their properties. Currently, associations have mainly been investigated in the bulk, while little attention has been given to the electrolyte structure near electrified interfaces. Here, we develop a theory for the electrical double layer (EDL) of WiSEs, where we consistently account for the thermoreversible associations of species into Cayley tree aggregates. The theory predicts an asymmetric structure of the EDL. At negative voltages, hydrated Li + dominates, and cluster aggregation is initially slightly enhanced before disintegration at larger voltages. At positive voltages, when compared to the bulk, clusters are strictly diminished. Performing atomistic molecular dynamics (MD) simulations of the EDL of WiSE provides EDL data for validation and bulk data for parametrization of our theory. Validating the predictions of our theory against MD showed good qualitative agreement. Furthermore, we performed electrochemical impedance measurements to determine the differential capacitance of the studied LiTFSI WiSE and also found reasonable agreement with our theory. Overall, the developed approach can be used to investigate ionic aggregation and solvation effects in the EDL, which, among other properties, can be used to understand the precursors for solid-electrolyte interphase formation.

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

Ni and Co Incorporation in Forsterite: A Density Functional Theory Study with Hubbard Correction

Ni and Co are critical elements for the world economy and modern technologies. Mafic and ultramafic deposits represent low-grade yet abundant alternatives to traditional Ni and Co ores. Here, in this work, density functional theory (DFT) with the Hubbard U correction (DFT+U) was used to simulate the incorporation of Ni and Co in forsterite (Mg 2 SiO 4 ), the Mg endmember of olivine, a common mineral in mafic and ultramafic rocks. Hubbard U terms for Ni and Co were parametrized using a series of oxide, hydroxide, carbonate, silicate, and sulfide minerals relevant to extraction and recovery of Ni and Co from mafic and ultramafic deposits. Electronic, energetic, magnetic, and structural properties were considered in the parametrization. For each of Ni and Co, an effective Hubbard correction (U eff ) value that optimized agreement with either experimental data or a hybrid exchange-correlation functional for all of the minerals considered is reported. DFT+U ab initio molecular dynamics (AIMD) simulations of Ni and Co incorporated into the M1 and M2 octahedral sites of forsterite were then performed. Ni and Co substitution in the M1 site was more energetically favorable than substitution in the M2 site, in agreement with published partition coefficients. AIMD trajectories were used to compute extended X-ray absorption fine structure (EXAFS) spectra of Ni in the M1 and M2 sites for direct fitting to a published experimental spectrum of Ni in a natural San Carlos olivine sample. The results of the fit indicated that ordering of Ni in the M1 site was not as strong at the low Ni concentrations relevant to mafic and ultramafic silicate minerals as that at the higher concentrations of the Ni-Mg olivine solid solutions studied to date.

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

Validated Reactive Force Field Quantifies MXene Interfacial Properties, Mechanics, and Thermal Transport

MXenes combine rich surface chemistry, mechanical strength, and high conductivity for a multitude of emerging applications. Predictive modeling supports accelerated materials designs and has been limited by the absence of validated and transferable force fields. Here, we introduce an interpretable, reactive INTERFACE force field (IFF and IFF-R) for Ti 3 C 2 T x MXenes that is trained based on chemical knowledge and achieves quantitative agreement with experiments across lattice parameters (<0.5%), density (<0.2%), liquid contact angles, Raman spectra, and the in-plane elastic modulus (∼320 GPa). The models cover surface terminations from hydroxyl (−OH) to fluorine (−F) groups and are extensible to other chemistries. We introduce pH-resolved surface chemistry and identify dopamine adsorption mechanisms at MXene–aqueous interfaces supported by QCM-D and UV–Vis experiments. The data reveal coplanar and perpendicular binding modes and concentration-dependent multilayer assembly. We predict previously inaccessible properties, including termination-dependent cleavage energies, interlayer shear moduli and dynamic shear failure, nanoindentation and brittle fracture, anisotropic in-plane and out-of-plane thermal conductivities, including the role of defects. Agreement with available experimental data is consistently close and exceeds DFT accuracy across the benchmark properties examined. The IFF/IFF-R model is compatible with CHARMM, AMBER, OPLS, and CVFF force fields for simulations of MXenes with diverse surface terminations, electrolyte interfaces, biointerfaces, and polymer composites without additional parameters. Parameter sets, 3D models, and analysis scripts are provided for community use. The validated, reactive, and transferable IFF framework facilitates predictive design of MXene-based films, membranes, sensing interfaces, and composites.

MXene↗
DOE OSTI2020

Quantitative Determination of Biomass-derived Renewable Carbon in Fuels from Coprocessing of Bio-oils in Refinery Using a Stable Carbon Isotopic Approach

Increasing renewable carbon incorporation into conventional fuels through coprocessing with vacuum gas oil (VGO, a petroleum refining feedstock) is a critical step in biofuels development, scaling-up, adoption and associated GHG reduction. Optimization of the co-processing parameters maximizes incorporation of the renewable carbon in the fuel products. Quantitative determination of the renewable carbon content in the co-processed products provides direct evaluation of the parameters. The co-processing bio-oil with VGO through hydrocracking (HC) or fluid catalytic cracking (FCC) system resulted in carbon isotopic fractionation that prevented the direct use of the isotope mixing model for quantifying the renewable carbon. Here, we report an algorithm of using a stable carbon isotope approach to quantify the renewable carbon content in co-processing biofuel products through high-precision ?13C analysis. A controlled experiment carried out by blending a fossil diesel (-29.013‰) with a bio-diesel (-30.099‰) at various blending levels up to 98.0/2.0 wt% is presented and has demonstrated the applicability of this approach. The carbon isotope fractionation factors for the bio-oil co-processing were obtained by using a 14C-derived isotope-mixing model. The ?13C method was tested by co-processing 13C-labeled bio-crude and natural woody biomass-derived fast pyrolysis (FP) and catalytic fast pyrolysis (CFP) bio-oils with VGO. The results were verified by 14C accelerator mass spectrometry (AMS) method (ASTM-D6866) and compared with the yield mass balance (YMB) method. Strong agreement between d13C and 14C AMS methods demonstrated the applicability of the ?13C method to quantify renewable carbon content in co-processing fuel products and guide the co-processing optimization

Li, Zhenghua↗
DOE OSTI2025

Molecular-Level Insights into the NMR Relaxivity of Gadobutrol Using Quantum and Classical Molecular Simulations

MRI is an indispensable diagnostic tool in modern medicine; however, understanding the molecular-level processes governing NMR relaxation of water in the presence of MRI contrast agents remains a challenge, hindering the molecular-guided development of more effective contrast agents. By using quantum-based polarizable force fields, the first-of-its-kind molecular dynamics (MD) simulations of Gadobutrol are reported where the 1 H NMR longitudinal relaxivity r 1 of the aqueous phase is determined without any adjustable parameters. The MD simulations of r 1 dispersion (i.e., frequency dependence) show good agreement with measurements at frequencies of interest in clinical MRI. Importantly, the simulations reveal key insights into the molecular level processes leading to r 1 dispersion by decomposing the NMR dipole–dipole autocorrelation function G(t) into a discrete set of molecular modes, analogous to the eigenmodes of a quantum harmonic oscillator. The molecular modes reveal important aspects of the underlying mechanisms governing r 1 , such as its multiexponential nature and the importance of the second eigenmodal decay. By simply analyzing the MD trajectories on a parameter-free approach, the Gadobutrol simulations show that the outer-shell water contributes ∼50% of the total relaxivity r 1 compared to the inner-shell water, in contrast to simulations of (nonchelated) gadolinium-aqua where the outer shell contributes only ∼15% of r 1 . The deviation between simulations and measurements of r 1 below clinical MRI frequencies is used to determine the low-frequency electron-spin relaxation time for Gadobutrol, in good agreement with independent studies.

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

Data Assimilative Optimization of WSA Source Surface and Interface Radii using Particle Filtering

TheWang-Sheeley-Arge (WSA) model estimates solar wind speed and interplanetary magnetic field polarity in the inner heliosphere using global photospheric magnetic field maps. WSA employs the Potential Field Source Surface (PFSS) and Schatten Current Sheet (SCS) models to determine the Sun's global coronal magnetic field. The PFSS and SCS models are connected through two radial parameters, the source surface and interface radii, which specify the overlap region between the inner SCS and outer PFSS models. Though both radii values are adjustable, they have typically been fixed to 2.5 solar radii. Our work highlights how solar wind predictions improve when the radii are allowed to vary over time. Data assimilation using particle filtering (sequential Monte Carlo) is used to infer optimal values over a fixed time window. Solar wind model predictions and satellite observations are compared with a newly developed quality-of-agreement prediction metric. The agreement metric between the model and observations is assumed to correspond to the probability of the two key WSA model parameters, the source surface and interface radii, where the highest metric value implies the optimal radii.We find that the optimal particle filter values of solar radii can perform twice as well as standard values for an exploratory period during Carrington Rotation 1901, with these values also reducing nonphysical kinking effects seen in solar magnetic field lines. Data assimilation choices of input realization and time frame have implications for variation in the solar wind over time. We present this work's theoretical context and practical applications for prediction accuracy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Predictive Inverse Model for Advective Heat Transfer in a Short–Circuited Fracture: Dimensional Analysis, Machine Learning, and Field Demonstration

Identifying fluid flow maldistribution in planar geometries is a well–established problem in subsurface science/engineering. Of particular importance to the thermal performance of enhanced (or “engineered”) geothermal systems is identifying the existence of nonuniform (i.e., heterogeneous) permeability and subsequently predicting advective heat transfer. Here, machine learning via a genetic algorithm (GA) identifies the spatial distribution of an unknown permeability field in a two–dimensional Hele–Shaw geometry (i.e., parallel plates). The inverse problem is solved by minimizing the L2 norm between simulated residence time distribution (RTD) and measurements of an inert tracer breakthrough curve (BTC) (C–Dot nanoparticle). Principal component analysis (PCA) of spatially correlated permeability fields enabled reduction of the parameter space by more than a factor of 10 and restricted the inverse search to reservoir–scale permeability variations. Thermal experiments and tracer tests conducted at the mesoscale Altona Field Laboratory (AFL) demonstrate that the method accurately predicts the effects of extreme flow channeling on heat transfer in a single bedding–plane rock fracture. However, this is only true when the permeability distributions provide adequate matches to both tracer RTD and frictional pressure loss. Without good agreement to frictional pressure loss, it is still possible to match a simulated RTD to measurements, but subsequent predictions of heat transfer are grossly inaccurate. Here, the results of this study suggest that it is possible to anticipate the thermal effects of flow maldistribution, but only if both simulated RTDs and frictional pressure loss between fluid inlets and outlets are in good agreement with measurements.

42 ENGINEERING↗
DOE OSTI2021

Linear Response Function Reveals the Most Effective Remote Forcing in Causing September Arctic Sea Ice Melting in CESM

We apply the linear response function method to investigate the most excitable mode of the September Arctic sea ice and its corresponding remote oceanic forcing in climate models. The method is useful in identifying the fundamental, dynamical and causal relationship between the remote forcing and sea ice response. We find that the most excitable mode of the September Arctic sea ice preferentially takes place over the Pacic side of the Arctic and its remote forcing corresponds to a dipole pattern of precipitation anomaly in the tropics with an increase of precipitation over the western and central tropical and subtropical Pacic ocean while a decrease over the Maritime Continent. The tropical precipitation anomaly likely drives a Rossby wave train propagating toward higher latitudes and leads to a ridge anomaly over the Pacic side of the Arctic, resulting in poleward atmospheric heat transport, enhanced downward longwave radiation and thus melting of the sea ice. In addition, a good agreement is found with the leading tropical Arctic teleconnection mode detected in a pre-industrial control simulation, which supports the usefulness and robustness of the linear response function method and the importance of tropical forcing in shaping the sea ice variability.

Wu, Yutian↗
DOE OSTI2022

Are Terrestrial Biosphere Models Fit for Simulating the Global Land Carbon Sink?

The Global Carbon Project estimates that the terrestrial biosphere has absorbed about one-third of anthropogenic CO 2 emissions during the 1959–2019 period. This sink-estimate is produced by an ensemble of terrestrial biosphere models and is consistent with the land uptake inferred from the residual of emissions and ocean uptake. The purpose of our study is to understand how well terrestrial biosphere models reproduce the processes that drive the terrestrial carbon sink. One challenge is to decide what level of agreement between model output and observation-based reference data is adequate considering that reference data are prone to uncertainties. To define such a level of agreement, we compute benchmark scores that quantify the similarity between independently derived reference data sets using multiple statistical metrics. Models are considered to perform well if their model scores reach benchmark scores. Our results show that reference data can differ considerably, causing benchmark scores to be low. Model scores are often of similar magnitude as benchmark scores, implying that model performance is reasonable given how different reference data are. While model performance is encouraging, ample potential for improvements remains, including a reduction in a positive leaf area index bias, improved representations of processes that govern soil organic carbon in high latitudes, and an assessment of causes that drive the inter-model spread of gross primary productivity in boreal regions and humid tropics. The success of future model development will increasingly depend on our capacity to reduce and account for observational uncertainties.

58 GEOSCIENCES↗
DOE OSTI2023

Cloud Condensation Nuclei Closure Study Using Airborne Measurements Over the Southern Great Plains

Abstract Airborne measurements of non‐refractory bulk aerosol chemical composition, aerosol size distributions, and cloud condensation nuclei (CCN) were conducted onboard a research aircraft during the Holistic Interactions of Shallow Clouds, Aerosols and Land Ecosystems field campaign in the spring and summer of 2016. A CCN closure study for the entire campaign period was performed where measured CCN concentrations at 0.24% and 0.46% supersaturation were compared with the predicted CCN concentrations calculated using κ ‐Köhler theory, three different assumptions of aerosol mixing state, and various assumptions about hygroscopicity, density, and the insoluble fraction of organic particles. We found that the Closure Ratio (CR), calculated as the ratio of predicted to measured CCN concentrations, was equal to one under two different aerosol mixing state assumptions: (a) all particles are composed of 100% organic particles, and (b) particles are externally mixed and composed of pure sulfates, nitrates, and organic particles assuming hygroscopicity values for organic particles ( ) between 0.04 and 0.17. The use of internal mixing state assumption often led to overprediction of CCN concentrations but the agreement within ±20% with the measured CCN concentrations was observed under certain closure permutations. A similar agreement, that is, within ±20%, was also observed using permuted parameters concerning density (1 and 1.5 g cm −3 ) and an insoluble fraction (0% and 20%) of organic particles. These findings may provide constraints on to predict CCN concentrations at a remote continental Southern Great Plains site.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Methane and Carbon Dioxide Fluxes in a Temperate Tidal Salt Marsh: Comparisons Between Plot and Ecosystem Measurements

Tidal wetlands are comprised of complex interdependent pathways where measurements of carbon exchange are often scale dependent. Common data collection methods (i.e., chambers and eddy covariance) are inherently constrained to different spatial and temporal scales which could generate biased information for applications of carbon accounting, identifying functional relationships and predicting future responses to climate change. Consequently, it is needed to systematically evaluate measurements derived from multiple approaches to identify differences and how techniques complement each other to reconcile interpretations. For our study, to accomplish this, we tested ecosystem-scale eddy covariance with plot-scale chamber measurements within a temperate salt marsh. We found good agreement (R 2 = 0.71–0.95) when comparing measurements of CH 4 emissions and CO 2 exchange but this agreement was dependent upon canopy phenology with discrepancies mainly arising during senescence and dormancy phenophases. The environmental drivers for CH 4 and CO 2 fluxes were mostly preserved across different measurement techniques, but the number of drivers increases while their individual strength decreases at the ecosystem scale. Empirical upscaling models parameterized with chamber measurements overestimated annual net ecosystem exchange (NEE; 108%) and gross primary production (GPP; 12%) while underestimating ecosystem respiration (Reco; 14%) and CH 4 emissions (69%) compared to eddy covariance measurements. Our results suggest that the environmental complexity of CH 4 and CO 2 fluxes in salt marshes may be underestimated by chamber-based measurements, and highlights how different techniques are complementary while considering limitations at each level of measurement.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Evaluation and Optimization of Snow Albedo Scheme in Noah-MP Land Surface Model Using In Situ Spectral Observations in the Colorado Rockies

The Biosphere-Atmosphere Transfer Scheme (BATS) ground snow albedo algorithm is commonly used in land-surface models (LSM), weather forecasting and research applications. This study addresses key uncertainties in BATS simulated ground snow albedo within the Noah-MP LSM framework through evaluation and optimization of the Noah-MP BATS ground snow albedo formulation using 2-band (visible and near-infrared (NIR)) in situ albedo observations at Rocky Mountain field stations. The Noah-MP BATS ground snow albedo scheme is extremely sensitive to its input parameters. Namely, an ensemble generated by varying BATS input parameters within potentially plausible ranges provides an average daily range (maximum ensemble member minus minimum ensemble member) of ground snow albedo exceeding 0.45 in visible and NIR bands. Parameter optimization improves agreement between simulated and in situ observed ground snow albedo in visible, NIR and broadband spectrums. Importantly, optimized parameters result in reduced biases relative to observed fresh-snow albedo and better agreement with observed albedo decay. Our analysis across different sites supports that the optimized BATS ground snow albedo parameters are appropriate to transfer in space and time, at least within the region studied (the central-southern Rocky Mountains). The primary error source remaining after parameter optimization is that observed fresh-snow albedo is highly variable, particularly in the NIR spectrum, whereas BATS fresh-snow albedo is constant, an issue which requires further investigation. This study shows significant correlations between observed fresh-snow albedo and surface meteorological conditions (e.g., downward shortwave radiation and temperature) which can support future model development that attempts to include a time-varying formulation for fresh-snow albedo.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Statistical Study of the Nonthermal Continuum Radiation Beaming Angle Measured by the High Frequency Receiver on Van Allen Probes-A

The nonthermal continuum (NTC) radiation beaming angle is computed over the entire Van Allen Probes-A mission when the spacecraft was in the dawn sector. The conditions in the dawn sector are favorable for the wave vector to lie near/in the spacecraft's spin plane allowing a favorable estimate of the beaming angle, and the dawn sector is also advantageous in that previous studies show NTC occurrence to peak in this sector. We found that scatter plots, over the entire mission, of beaming angle versus magnetic latitude form a distinct inverted-V pattern, with the apex at/near the magnetic equator. This pattern was sharpest for frequencies ( f ) ≲ 100 kHz. Using the NTC beaming formula from the linear mode conversion theory (LMCT), we show that such an inverted-V pattern is expected due to the large variation in the plasmapause location over the entire mission. The theoretical derived pattern qualitatively reproduces the observed pattern but not quantitatively. The lack of quantitative agreement is discussed and is attributed to several factors, one factor is off-centered emissions from the radio window. The qualitative agreement strongly supports LMCT as being the dominant mechanism generating NTC for f ≲ 100 kHz. For f ≳ 100 kHz, the inverted-V pattern becomes less distinct, and strong near-equatorial beaming is observed. Finally, after considering contamination of our selections by left-handed polarized AKR, our study suggests that besides LMCT another unidentified NTC generation mechanism becomes important for f ≳ 100 kHz.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Cascadia Subduction Zone Fault Heterogeneities From Newly Detected Small Magnitude Earthquakes

The Cascadia subduction zone (CSZ) is known to host M9 megathrust ruptures; however, no such event has occurred in historical observation. The distribution and characteristics of small- to moderate-sized earthquakes can be used to determine the behavior of the megathrust fault but are notably absent offshore the CSZ due to the distance from onshore seismometers. We use automated subspace detection coupled with an onshore-offshore seismic deployment to find small-magnitude earthquakes in the offshore seismogenic zone and analyze their locations in the context of interseismic locking and seismogenic zone extent. Here we detected and located 5,282 earthquakes, 4,096 of which had been previously undetected. We find that the downdip extent of the seismogenic zone as defined by interplate seismicity agrees with the 20% locking contour of the Schmalzle et al. (2014, https://doi.org/10.1002/2013GC005172) geodetic model and extends deeper than predicted by previous thermal models. We cannot determine the updip extent of the seismogenic zone; this may be due to a lack of templates for detection in the updip source area, stress shadows updip of asperity loading, and/or strong locking to the trench. We present a map of possible asperities determined by the small earthquakes in this study. Our asperity locations and extents show some, but not complete, agreement with the asperities modeled from the 1700 M9 rupture and geodetic locking models, and good agreement with the paleo-rupture extents determined from offshore turbidites and forearc basin-based asperity estimates. This highlights the need of continued offshore observations over time, and to elucidate fine-scale variation in locking.

58 GEOSCIENCES↗