SEARCH · Engineering Papers
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.
Final report on the operations, data collection, archival and dissemination of earthquake observations from borehole arrays in California, the Pacific Northwest, and Alaska during the project period from August 2019 through July 2022.
The new WAVEFORMS Initiative in the Ground-based Nuclear Detonation Detection (GNDD) program includes an increased emphasis on the development of Earth models and methods to predict entire seismic and acoustic waveforms more accurately. In general, this increased emphasis is predicated on the need to better characterize seismic events and provide improved model-based discrimination between event types including earthquakes and explosions. More specifically, while current moment tensor inversion methods tend to work well for larger events (M>~4) using tuned 1-D Earth models, the development of state-of-the-art 3-D models and methods is required for the prediction of shorter period waves over large areas for discrimination of smaller events. There is no standard metric for model-based waveform prediction accuracy used in the waveform modeling/inversion community. However, there are several popular waveform misfit definitions; and minimizing the corresponding objective functions is the goal of waveform inversion. Some example misfit definitions employed for adjoint waveform tomography include measures of simple travel time differences (e.g. Tape et al., 2010), cross-correlation travel time differences (e.g. Luo and Schuster, 1991), multi-taper frequency dependent methods (e.g. Lei et al., 2020), time-frequency phase misfit functions (e.g. Fichtner 2010; Rodgers et al., 2022), normalized cross-correlation methods (e.g. Tao et al., 2018), and others. In some cases, these misfit definitions also involve complicated weighting schemes and summations over multiple frequency bands making it difficult to duplicate the misfit measurement with alternative models and datasets. Although each of the misfit definitions mentioned above are useful for developing waveform models, the actual misfit values are not usually meaningful outside of a given project, model, and/or dataset. Therefore, it is difficult to understand and communicate model performance for predicting waveforms and comparing to other models and/or new model iterations with a different dataset. Therefore, there is a need for a generalized method for evaluating overall model performance that is independent from the specific misfit chosen to develop the waveform models that is also intuitive and meaningful. In this report, we describe a new metric we refer to as ‘Percent of Correlated Signal’. The following sections describe and demonstrate the metric with a case study event and a more rigorous test using a random selection of globally distributed events. While the focus here is on global tomography models, the metric is meant to applicable to regional ‘wiggle-for-wiggle’ waveform models/studies as well.
The Source Physics Experiments (SPE) were designed to improve our physics-based understanding of explosion sources for the purposes of nuclear test monitoring. Phase I consisted of 6 chemical explosions in the Climax Stock Granite of the Nevada National Security Site (NNSS), while Phase II consisted of 4 explosions in a contrasting dry alluvium geology (DAG) in Yucca Flat, providing essential data in various media and emplacement conditions to further modeling efforts. For Phase III, the Rock Valley Direct Comparison (RVDC) seeks to directly compare earthquake and explosion source types. An unusually shallow series of events in 1993 along the Rock Valley Fault Zone in the southeastern portion of the NNSS has been targeted for this direct comparison. Depth ranges for the events, previously estimated to be less than 3 km, is achievable by modern drilling techniques and accessibility to the epicentral locations would require minimal improvements to the infrastructure. The events providing this unique opportunity for direct comparison are the focus of this report.
Geophysical research on historical nuclear tests is an important aspect of future monitoring capabilities in seismic research. This research is challenging due to the limited number of digital seismic recordings during the peak of nuclear testing (1945-1992). These limited records are unique and non-reproducible data with potential high research impact. Releasing available nuclear explosion seismic records to the explosion monitoring community is thus of high value and is the motivation for this dataset release. The target of this effort was on compilation and quality control of regional seismic records of nuclear explosions recorded on Lawrence Livermore National Laboratory stations ELK, KNB, LAC, and MNV, known collectively as the Livermore National Network (LNN) (Figure 1). LNN was established in the early 1960s for the primary purpose of monitoring underground nuclear testing at the former Nevada Test Site (NTS), now known as the Nevada Nuclear Security Site (NNSS) following the signing of the Limited Test Ban Treaty (LTBT). LNN consisted initially of short-period vertical component Benioff’s recorded on film located at Mina, NV (MNV) and Kanab, Utah (KNB). LNN added two additional stations at Landers, CA (LAC) and Elko, NV (ELK) in 1967 and upgraded equipment to broadband seismometers recorded on frequency modulation (FM) tapes from 1967-1979, followed by digital recordings after 1979 (Jarpe, 1989). The digital recordings were on a variety of now obsolete media, including 9-track, Exabyte, and DAT tapes. Jarpe (1989) describes the seismic station instrumentation details over the period of deployment. LNN recorded valuable non-repeatable unique data of several hundreds of nuclear explosions at NNSS, as well as earthquakes and chemical and mining explosions (Walter, 2020). The details of these nuclear tests are provided in the Department of Energy Report NV-209 Rev 16 (DOE, 2015).
Our goal is to develop, apply and validate a holistic thermal, hydrologic, mechanical, and chemical (THMC) workflow that includes evaluation of induced seismic slip in EGS reservoirs. We will integrate experimental and modeling approaches to reduce parameter uncertainty and better predict/mitigate seismic hazard at EGS sites. Our novel approach couples 3D physics-based earthquake simulations with THMC models (THMc+E). This capability will enable improve engineering decisions at Utah-FORGE and move EGS operations toward repeatable, robust, economically viable, and socially accepted development. For example, our THMC+E models will predict circulation scenarios and related seismic hazard for a suite of flow rates and under uncertainty, thus enabling evaluation of optimal circulation strategy. Laboratory experiments will be performed to constrain key model parameters and Bayesian techniques will provide a probabilistic evaluation of parameters used in models. THMC+E simulations will enable exploration various circumstances that may hinder EGS success and develop mitigation strategies.
Seismic source characterization is an essential function of global nuclear explosion monitoring (NEM). While large events (roughly with moment magnitude, M w , greater than 5.0) can often be easily detected, located and identified with high signal-to-noise ratios at teleseismic distances (> 20°), trends in NEM research require confident source characterization at much lower magnitudes (say down to 3.0) and exploitation of sparse observations (from only a few stations) at regional distance (< 20°). Regional distance waveform inversion to characterize sources is now widely used and effective (e.g. Ford et al., 2009; Alvizuri and Tape, 2018; Alvizuri et al., 2018; Chiang et al., 2018; Ford et al., 2022). These methods obtain the magnitude, depth and seismic moment tensor, which represents the forces that excited the observed seismic waves (slip on an earthquake fault, explosion, collapse or a combination of various forces). Common to many problems in seismology, the isolation of the source 2 properties requires removal of path propagation effects that waves experience while traveling through the three-dimensional (3D) Earth (the structure exists due to different rock types, material properties, temperature and tectonic processes).
DAS records a variety of signals, including (1) earthquakes and volcanic eruptions, (2) explosions (DAG, mining, etc.), (3) acoustic waves, (4) ocean (waves, whales), (5) vehicle, ship and pedestrian traffic, (6) glaciology and landslides, and (7) oil and gas. Main questions include: (1) what is the extent of the signal content? (2) how does DAS compare to traditional sensors? (3) how do interrogators, fiber types, and fiber depth burials impact measurements? (4) special applications offshore to overcome the lack of instrumentation.
Characterizing the shallow structure of the Rock Valley region of the Nevada National Security Site is a critical component of the Rock Valley Direct Comparison project. Geophysical data of the region is needed for operational decisions, to constrain geologic models used for simulation, and to facilitate the analysis of future explosive source data. Local measurements of gravity are a key piece of geophysical information that helps to resolve the underlying geologic composition, fault structure, and density characteristics, yet, in the Rock Valley region these measurements are sparse on the scale of the testbed. In this report, we present the details of a recent gravity data acquisition survey designed to collect a dense dataset in the region of interest that complements the existing gravity work but greatly enhances our resolution. This dataset will be integrated with a complementary Los Alamos National Laboratory gravity collection and combined with the existing seismic data in a joint inversion. These measurements were conducted over two weeks with a portable gravimeter and high-resolution GPS and include repeat measurements at a USGS base station as well as reoccupation of gravity sites in the regional dataset. This collection of over 100 new dense gravity measurements will facilitate refinement of the existing Geologic Framework Model and directly complement newly acquired dense seismic data, ultimately improving the project’s ability to investigate the direct comparison of shallow earthquake and explosive sources.
The LA100 Equity Strategies project integrates community guidance with robust research, modeling, and analysis to identify strategy options that can increase equitable outcomes in Los Angeles' clean energy transition. As Los Angeles transitions toward clean energy, existing distribution grid infrastructure will need to be updated and expanded to support reliable service during routine operations, enable interconnection with distributed energy resources and electrified loads, and provide access to energy-related services during disasters. This chapter focuses on equity in distribution grid upgrades, reliability, and resilience in Los Angeles. Specifically, NREL performed grid upgrade and resilience analyses using a detailed model of the distribution grid and income-differentiated household load profiles, electric vehicle (EV) adoption patterns, distributed solar adoption, and grid reliability to explore two key questions to inform how the City of Los Angeles can ensure a resilient and reliable distribution grid for all communities during the clean energy transition: Where can distribution system upgrades can be prioritized to enable equitable access to, and adoption of, clean energy technologies and how can Los Angeles provide equitable, resilient access to electricity-related services (e.g., health care, food) during disaster events like earthquakes and flooding? The electric distribution system is the "last mile" of the grid, linking the multistate bulk power system with customers; new loads, including EVs; and distributed energy resources, such as customer and community solar and storage. This analysis focuses on the 4.8-kilovolt (kV) system, including service transformers that represent the utility-side of the grid connection for most residential customers. Chapter 17 looks at the customer-side of the grid connection with a focus on electric panel upgrade needs. The transition toward clean energy can put additional stress on the distribution system from distributed energy resources and electrification - especially EVs and increased use of electricity for heating, cooling, cooking, and hot water. This stress, measured here as the number of equipment overloads and voltage violations, correlates strongly to grid reliability and therefore is used as a proxy for understanding additional upgrades needed and to help ensure equitable access to electrification and distributed energy resources. NREL also conducted community resilience analysis to examine customer-level access to both electricity and a larger range of services, such as hospitals and grocery stores during a disaster. This analysis explicitly considers equity to understand differences in current resilience and resilience strategies to effectively improve critical services access for all Angelenos. Research was guided by input from the community engagement process, and associated equity strategies are presented in alignment with that guidance.
On April 1, 1946, an earthquake off Unimak Island, Alaska, generated a tsunami that raced across the Pacific. One of the islands in danger was Kwajalein, the logistical hub for Operation Crossroads. Commodore Ben Wyatt, the U.S. naval commander of the Marshall Islands, ordered the evacuation of “10 Red Cross workers, 19 USO women, 17 native women, 16 male [hospital] patients, four tiny babies, one small boy, and two navy doctors,” to the USS Rockbridge, which then steamed in deep water off the atoll until the danger passed. One of the Red Cross workers, Reva Hurwitz, played bridge in the ship’s wardroom until 2 a.m., a luxury since the island curfew, strictly enforced, was 10:30 pm. The tsunami bypassed Kwajalein, and everyone returned to the island less than twenty-four hours after being evacuated. Although the evacuation proved unnecessary, the women were grateful, said Hurwitz, because they had the opportunity “to use hot water and eat fresh fruit.”
As significantly more focus is placed on carbon sequestration as a means to help mitigate the progression of climate change, the need to adequately characterize an injection target for its induced seismicity potential becomes more apparent. In the Iron Mountain project in Southern Utah, a single fault is identified in the Neck of the Desert near Iron Mountain. The dip of the fault is unknown, and the geomechanical stability of the fault (i.e., distribution of Coulomb failure functions (CFFs)) is highly dependent upon the dip. For some dips, such as 15 degrees, the hazard from injection is low. For other dips, such as 75 degrees, the hazard is high. Another key consideration is the depth of penetration of the fault. Should the fault not penetrate the stress and pressure perturbed region from injection, it will not be induced to fail. Thus, understanding the dip and the depth of penetration of the fault will greatly improve our understanding of the likelihood of inducing an earthquake.
Togiak is one of hundreds of Alaskan Native, roadless communities only accessible by boat or plane. This project constructed a heat recovery (HR) system which pipes what was previously wasted heat generated at the power plant to four nearby community buildings. The power plant houses diesel fired generators which provide electricity to the community. This project designed and constructed the system that recovers a portion of the heat produced by the generators and pipes it to the end-user buildings. The heat recovery system integrates with each building’s internal heating system. The objective is to significantly reduce the heating costs in these four community facilities. By reducing the energy use and cost in these community facilities, the local economy is strengthened by keeping the funds inside the community that were formerly used to purchase imported fuel, thereby multiplying the impact of the energy cost reduction. Energy independence is enhanced and energy resilience is increased by reducing dependence on long fuel supply chains which are easily interrupted by unpredictable natural and man-made disasters such as earthquakes, climate change, road rail and sea accidents, etc.
As geologic carbon sequestration projects begin to be funded with a higher degree of frequency, microseismic monitoring will become more necessary to establish caprock integrity and induced seismicity risk. The Farnsworth, TX site, which hosts enhanced oil recovery operations, provides an ideal laboratory to test essential components of microseismic monitoring such as surface station placement, denoising, and different autopicker methods. We find that despite careful optimization, the surface station placement at the Farnsworth site was insufficient given the high level of background noise from the industrial operations at the oil field. Thus, we use only borehole geophone data in our processing. Two denoising techniques were examined: the DeepDenoiser and the continuous wavelet transform. The continuous wavelet transform was shown to be a valuable tool in converting between raw waveforms to processed denoised waveforms for microseismic monitoring. In the case of the noisy waveforms, only ten detections are found, compared with 90 in the denoised data for a two-hour window. The DeepDenoiser suffered from the fact that the training data was regional earthquake data. In addition, the PhaseNet machine-learning autopicker was applied to both the noisy and the denoised data, and this algorithm detected thousands of more arrivals in the data denoised with the continuous wavelet transform technique.
In light-water reactors (LWRs), the concrete biological shield (CBS) designates the concrete structure directly facing the reactor pressure vessel (RPV). The primary function of the CBS is to protect equipment and personnel from the neutron and gamma radiation exiting the RPV. Most CBSs in operation in the United States also provide structural support to the RPV. Based on a long-term structural analysis, the structural integrity of the CBS to transfer the in-service passive load to the foundation, dynamic loading during an earthquake, and thermal loading during a loss-of-coolant accident (LOCA) must be maintained. Based on test reactor data for accelerated conditions, the mechanical properties of irradiated concrete are affected when exposed to neutron fluence higher than ≈ 10 19 n.cm -2 (E > 0.1 MeV). At 80 years of operation, the fluence estimates at the surface of the CBS range between 1×10 19 and 7×10 19 n.cm -2 (E > 0.1 MeV): the fluence value depends mostly on the design. Although fluence attenuates through the CBS by virtue of the shielding properties of concrete, about 5% to 10% of the CBS wall depth is affected by high fluence. The main degradation mechanism is the so-called radiation-induced volumetric expansion (RIVE) caused by the amorphization of aggregate-forming minerals. RIVE causes large deformation of up to 18% in pure quartz, and it also reduces the structural properties of concrete. The structural effects of long-term exposure to fast neutron irradiation in LWRs was studied in the report entitled “Assessment of the Effect of the Irradiation-Induced Degradation on the In-Service Structural Performance of the Concrete Biological Shields” (M3LW-23OR0403013). This report illustrates that (1) irradiation-induced damage expands further than the region subject to a cumulative fluence of 10 19 n.cm -2 (E > 0.1 MeV) to reach the steel reinforcement location, and (2) the irradiation-induced damage region shows no substantial residual bearing capacity and should not be accounted for in the integrity assessment of the CBS.
Presents results from power grid simulations on synthetic grids that estimate the impact from power transformer failures due to earthquakes. In the best-case scenario, there is an 82.9% likelihood of load loss greater or equal to 51 MW. Methodologies are described as well as path forward to improving transformer resilience and, consequently, the resilience of the bulk power grid.
Ground motions recorded at one location are often extrapolated to nearby regions within a given radius. Here, we explore the appropriateness of spatial extrapolation using data from seven small aperture seismic network deployments in southern California. Six of these deployments are linear arrays of 4-13 stations, and one is a 2D array of 13 stations at Pinyon Flats Observatory. The spatial footprint array diameters are 3 km or less, and each array was operational for a year or more. From our base catalog (M2.5+ earthquakes; 4038 events; September 2010 - June 2023), automated methods remove temporally overprinted waveforms from nearby events (< 5 km) in quick succession (< 5 min) and data with nonviable waveforms. These 200 samples per second data are filtered at 0.5-25 Hz and must have signal-tonoise ratios (SNRs) of 2.5+. Peak ground acceleration (PGA) and peak ground velocity (PGV) are derived individually from the maximum absolute values of each of the 3-component waveforms (vertical, northsouth, and east-west). Five of the seven arrays traverse the San Jacinto fault, and two do not. Ground motion observations are compared with theoretical estimates from Abrahamson et al., 2014. On average, arrays deployed within and across fault zones consistently record ground motions above theoretical expectations, whereas off-fault arrays record ground motions at or slightly below theoretical expectations. We attribute these differences to site conditions because these trends prevail for the full data suites. For each network and each individual channel, the coefficient of variation indicates that the standard deviations are ~30±6% of the mean. Exploring relative ground motion contributions from all three channels (ternary plots), as expected, most data show that vertical ground motions are attenuated compared to horizontal ones. However, this is not always the case for RA array data, where vertical motions can be ~2-3 times larger than horizontal motions for select events near Cahuilla, CA. These anomalously high vertical motions are focal mechanism-related. These results suggest that ternary plots created using only a small amount of data can be used as a data quality metric and a tool to find anomalous features in three-component data.
To allow for on-site Potassium Dihydrogen Phosphate (KDP) crystal finishing capabilities, Lawrence Livermore National Laboratory (LLNL) acquired the Moore 3 diamond turning machine. Equipment such as the Moore 3 is required by LLNL to have adequate seismic bracing before permanent installation to reduce the risk of personnel injury or machine damage during an earthquake. The purpose of this project was to design a bracket that can withstand seismic forces induced by the machine without hindering function. To develop a solution, background research—including stakeholder identification, stakeholder interviews, and a literature survey—was conducted. Based on LLNL site and stakeholder requirements, several approaches were developed for bracing the Moore 3. Ideations were evaluated based on a weighted decision matrix and a final design strategy was selected. This design was an L-shaped weldment that featured provisions for four HILTI concrete anchors, three triangular gussets, and a steel construction.
Seismic traveltime anomalies of waves that traverse the uppermost 100–200 km of the outer core have been interpreted as evidence of reduced seismic velocities (relative to radial reference models) just below the core–mantle boundary (CMB). These studies typically investigate differential traveltimes of SmKS waves, which propagate as P waves through the shallowest outer core and reflect from the underside of the CMB m times. The use of SmKS and S(m- 1)KS differential traveltimes for core imaging are often assumed to suppress contributions from earthquake location errors and unknown and unmodelled seismic velocity heterogeneity in the mantle. The goal of this study is to understand the extent to which differential SmKS traveltimes are, in fact, affected by anomalous mantle structure, potentially including both velocity heterogeneity and anisotropy. Velocity variations affect not only a wave’s traveltime, but also the path of a wave, which can be observed in deviations of the wave’s incoming direction. Since radial velocity variations in the outer core will only minimally affect the wave path, in contrast to other potential effects, measuring the incoming direction of SmKS waves provides an additional diagnostic as to the origin of traveltime anomalies. Here we use arrays of seismometers to measure traveltime and direction anomalies of SmKS waves that sample the uppermost outer core. We form subarrays of EarthScope’s regional Transportable Array stations, thus measuring local variations in traveltime and direction. We observe systematic lateral variations in both traveltime and incoming wave direction, which cannot be explained by changes to the radial seismic velocity profile of the outer core. Moreover, we find a correlation between incoming wave direction and traveltime anomaly, suggesting that observed traveltime anomalies may be caused, at least in part, by changes to the wave path and not solely by perturbations in outer core velocity. Modelling of 1-D ray and 3-D wave propagation in global 3-D tomographic models of mantle velocity anomalies match the trend of the observed traveltime anomalies. Overall, we demonstrate that observed SmKS traveltime anomalies may have a significant contribution from 3-D mantle structure, and not solely from outer core structure.