SEARCH · Engineering Papers

Results for “seismic wave”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

NASA NTRS1974

High-frequency lunar teleseismic events

A small number of seismic signals, including some of the strongest observed to date, have been identified as representing a fourth principal category of natural lunar seismic events with characteristics distinct from those produced by normal meteoroid impacts, deep moonquakes, and thermal moonquakes. These signals are much richer in high frequencies than other events observed at comparable distances, and display relatively impulsive P- and S-wave beginnings, indicating negligible seismic-wave scattering near the source. Source depths of these events may range between 0 and perhaps 300 km. These and other characteristics could represent either (1) meteoroids impacting upon outcrops of competent lunar crystal rock, (2) rare impacting objects that penetrate to competent rock below a scattering zone, or (3) shallow tectonic moonquakes.

Nakamura, Y.↗

DOE OSTI2025

Hydrogen Transportation and Distributed Energy Systems Seismic Risk Assessment for Cascadia Subduction Zone Airport Facilities

Portland International Airport in Oregon is exploring operating a fleet of 28 fuel cell electric buses to support airport operations and possibly provide backup power during outages. In this report, we evaluate the risk of hydrogen deployment at an airport considering the potential for seismic activity. We present a seismic risk assessment for pressurized piping for a generic hydrogen refueling station at the airport. Only pressurized piping was considered in the risk analysis because of its vulnerability and availability of fragility information unlike the other components. Seismic capacity of other components should be included as more information becomes available. To characterize the seismic hazard, we considered the initiating event frequency of a Cascadia Subduction Zone earthquake event based on data from the literature. Seismic stressors that translate to the site from an earthquake event are expressed as peak ground acceleration. Site disturbance is characterized as a function of soil conditions, represented by different shear wave velocity values, which affect how seismic waves propagate and impact structures. Pipe fragility curves as a function of seismic capacity were used to correlate ground acceleration and failure probability. These fragility characteristics were combined with site-specific soil conditions to calculate the probability of pipe rupture during a seismic event. It should be noted that multiple simultaneous failures due to an earthquake event as a common cause are not considered. Some epistemic uncertainties such as self-exciting shocks, aftershocks, and the time variant nature of the ground motion are not considered either.

08 HYDROGEN↗

DOE OSTI2021

Listening for the Landing: Seismic Detections of Perseverance's Arrival at Mars With InSight

The entry, descent, and landing (EDL) sequence of NASA's Mars 2020 Perseverance Rover will act as a seismic source of known temporal and spatial localization. We evaluate whether the signals produced by this event will be detectable by the InSight lander (3,452 km away), comparing expected signal amplitudes to noise levels at the instrument. Modeling is undertaken to predict the propagation of the acoustic signal (purely in the atmosphere), the seismoacoustic signal (atmosphere-to-ground coupled), and the elastodynamic seismic signal (in the ground only). Our results suggest that the acoustic and seismoacoustic signals, produced by the atmospheric shock wave from the EDL, are unlikely to be detectable due to the pattern of winds in the martian atmosphere and the weak air-to-ground coupling, respectively. However, the elastodynamic seismic signal produced by the impact of the spacecraft's cruise balance masses on the surface may be detected by InSight. The upper and lower bounds on predicted ground velocity at InSight are 2.0 × 10 -14 and 1.3 × 10 -10 m s -1 . The upper value is above the noise floor at the time of landing 40% of the time on average. The large range of possible values reflects uncertainties in the current understanding of impact-generated seismic waves and their subsequent propagation and attenuation through Mars. Uncertainty in the detectability also stems from the indeterminate instrument noise level at the time of this future event. A positive detection would be of enormous value in constraining the seismic properties of Mars, and in improving our understanding of impact-generated seismic waves.

58 GEOSCIENCES↗

DOE OSTI2022

A Methodological Overview of Seismic Analysis for Nuclear Event Detection

Underground explosions generate potentially detectable signatures, including energy waves that travel through the Earth’s subsurface (i.e., seismic waves), low-frequency sound waves (i.e., infrasound and hydroacoustic waves), and radioactive gases and/or particles that might leak from the test cavity (if the event was nuclear). There can also be intelligence indicators of a test, such as observations of modified patterns of life and activity at a suspected test site. If all of these detectable signatures and intelligence indicators are present and self-consistent, then analysts have high confidence in classifying a signature generating event as an explosion. However, because only partial information about an event is likely to be available, determining whether an event was natural (e.g., an earthquake or landslide) or manmade (e.g., a chemical or nuclear explosion) is much more challenging. This primer describes how one category of event signatures—seismic signatures—can augment event analyses. While universities and government organizations have generated detailed technical descriptions of seismic analytic techniques, we seek to translate seismic event analysis for a broad, non-technical audience. When the geologic conditions near an event are well-characterized, seismic data can be used to calculate critical information, such as event location and depth, with relatively high accuracy. Moreover, specific features within seismic datasets can help determine whether an event was an explosion. However, a key challenge in seismic analysis is that geologic site conditions are often poorly characterized, complicating the ability to discern the true nature of the event. To overcome this challenge, geologists answer a series of questions (discussed in section 1) to guide seismic event analysis and determine the most probable nature of an event. As more information is gathered during each analytic step, confidence grows regarding the nature of the event. Section 2 addresses uncertainties in seismic analysis and the vital nature of high-fidelity geologic data for accurate seismic event analysis.

58 GEOSCIENCES↗

DOE OSTI2024

Applicability of 1D site response analysis to shallow sedimentary basins: A critical evaluation through physics-based 3D ground motion simulations

One-dimensional site response analysis (1D SRA) remains the standard practice in considering the effect of local soil deposits and predicting site-specific ground motions, although its range of applicability to realistic seismic wavefields is still in question. In this 1D approach, horizontal and vertical ground shaking are assumed to be induced by vertically propagating shear and compressional waves, respectively. A recent study based on analytical two-dimensional (2D) plane waves and simple point source earthquake simulations has shown two mechanistic limitations in this 1D modelling technique for general inclined seismic waves, that is, systematic over-prediction of the vertical motion and wave trapping in the 1D soil column. In this article, we evaluate in detail the applicability of this 1D modelling approach to realistic three-dimensional (3D) simulated seismic wavefields in shallow sedimentary basins. Linear-viscoelastic 1D SRA predictions using two types of input motions that are commonly used in practice—rock outcrop and in-column motions, are compared with the reference true site response results from 3D earthquake simulations in terms of various measures in the frequency and time domain. It is shown that the horizontal motion in the 3D seismic wavefield exhibits dominant shear wave propagation phenomenon, while the vertical motion is a combined effect of compressional and shear waves and can be over-predicted by the 1D approach when the incident seismic waves are inclined. Direct evidence of the wave refraction process that leads to the vertical motion over-prediction is provided. 1D SRA with in-column inputs can yield motions that have significantly longer duration compared to the true 3D site response solution due to trapped waves, casting in doubt the frequent need for increased soil damping in existing site studies to compensate for wave attenuation due to scattering alone. Sensitivity investigation on the increase of soil profile damping by a multiplier D mul shows D mul values compatible with those found in the literature for both horizontal and vertical motions. It is shown that the level of D mul optimized for a best match of the spectral acceleration is dependent on the characteristic of the input motion and a larger D mul is typically required for the vertical component. In contrast, 1D SRA with outcrop motions predicts motions with shorter significant duration due to its inability to capture the basin-edge generated surface waves. A suite of ground motion simulations was performed to assess the sensitivity of the observations to the basin geologic structure including the velocity gradient, rock-basin impedance contrast and basin depth. The analysis results show that the accuracy of the simplified 1D procedure is dependent on the wavefield composition of both the input motions and the true 3D site response solution. While the horizontal motions in shallow sedimentary basins can, to the first order, be reasonably captured by the simplified 1D approach, 1D SRA for the vertical component is in general not reliable and contributions from inclined shear waves should be accounted for in site-specific evaluation of the vertical design ground motion.

1D assumption↗

DOE OSTI2023

High-resolution lithospheric structure of continental China from joint inversion of surface wave and gravity data

The lithospheric structure of continental China has been previously determined by seismic travel time tomography, surface wave tomography, and joint inversion of body wave and surface wave data. However, due to the inherent limitations of seismic data, the lithospheric structure of continental China is still not well resolved in the shallow part and in some regions where the station coverage is relatively sparse. In this study, we aim at improving the lithospheric structure by joint inversion of seismic surface wave data and satellite gravity data to take advantage of the uniform distribution and complementary strength of the gravity data. Here, the empirical relationship between velocity and density is used as a bridge for joint inversion of surface wave and gravity data. The joint inversion shear-wave velocity Vs (density) model, named as USTClitho1.0g, can fit both surface wave and gravity data well. This high-resolution Vs model can better fit the active airgun source seismic arrival times and better delineate some features in continental China, such as the velocity contrast across the north–south gravity lineament (NSGL), the lithosphere thinning in eastern China, the crustal footprint of the Hainan mantle plume, the likely magma chamber beneath volcanos in northeast China, the middle-lower crust low velocity layer beneath the Tibetan plateau, and the tearing of subducted Indian mantle lithosphere. Our joint inversion Vs model can provide a reference model for geosciences in continental China and surrounding areas.

58 GEOSCIENCES↗

DOE OSTI2021

Determining fast-S and slow-S propagation directions with SV-P data produced by buried explosives and recorded with vertical geophones

We have evaluated the concept of practicing S-wave reflection seismology with legacy 3D seismic data generated by a P-wave source and recorded with only vertical geophones. This type of S-wave imaging is based on the principle that seismic P-wave sources not only produce a downgoing illuminating P wavefield, but they also simultaneously produce a downgoing illuminating SV wavefield that, in almost all cases, is suitable for S-wave reflection imaging. The S-mode used in this study is the SV-P, or converted-P, mode. This mode involves a downgoing illuminating SV wavefield and an upgoing reflected P-mode that is recorded by vertical geophones. In flat-layered stratigraphy, the lengths of the SV and P raypaths in SV-P imaging are identical to the lengths of the SV and P raypaths in P-SV imaging with P-sources and 3C geophones. P-SV imaging of deep rocks has been practiced for more than two decades; SV-P imaging is a new concept. SV-P data should provide the same options for investigating deep rocks as do P-SV data. We have determined one of the equivalences between SV-P data extracted from vertical-geophone data and P-SV data extracted from horizontal geophones: that both modes react to azimuth-dependent variations in the S velocity in anisotropic rocks. Azimuthal variations in the SV-P traveltime can be used to define the polarization direction of the fast-S-wave mode, which is also the azimuth of the maximum horizontal stress (SHmax). Our investigation demonstrates a noninvasive method for monitoring changes in the SHmax azimuth across a CO 2 storage reservoir, or any targeted porous rock, as fluids are cycled into, and then out of, that rock’s pore space.

Geochemistry & Geophysics↗

NASA NTRS1978

Lunar Seismology: the Internal Structure of the Moon

The direct P and S wave arrival times are the primary data set that can be measured on the seismograms of natural lunar seismic events. Polarization filtering techniques allow the enhancement of secondary body wave arrivals and record curves to identify the secondary phases and deduce structural information. Finally, shear wave amplitude vs. distance curves yield information on the location and magnitude of seismic velocity gradients in the interior. The results of these analyses show that the moon appears to have a two-layer crust at all four seismic stations: a 20 km upper crust that seems to be constant at all sites and a lower crust that is 40 km thick at stations 12 and 14 (mare), 55 + or - 10 km at station 16 (highland), and tentatively either 40 km or 70 km at station 15. The lower mantle extends from 480 km to at least 1100 km depth which is the maximum depth of penetration of all but a few seismic waves used as data. No definitive evidence for or against a lunar core exists.

Goins, N. R.↗

DOE OSTI2021

Compressional wave velocity and effective stress in unsaturated soil: Potential application for monitoring moisture conditions in vadose zone sediments

Compressional seismic (P-wave) velocities were measured experimentally at different water potentials and confining pressure states using a variety of sediments collected from the US Department of Energy Hanford Site in southeastern Washington State. P-wave velocity was measured in variably saturated sediments using ultrasonic piezoelectric transducers and water potential was measured with a heat dissipation sensor. We propose a model to relate changes in P-wave velocity to soil stress conditions using an effective stress formulation that includes capillary stress as well as adsorptive stress. We show that compressional wave velocities can be related to water potential and confining stress using a model with five fitted parameters. A single set of parameters provides an excellent fit to data acquired from all of the sediment samples. Under very dry conditions, corresponding to low water potentials, water flow occurs in adsorptive thin films associated with very low unsaturated hydraulic conductivities. The results presented here suggest the utility of seismic methods for establishing low water potential and flow conditions within variably saturated Hanford formation sediments and providing feedback on the performance of mitigation efforts such as soil desiccation and surface infiltration barriers. Seismic methods can provide minimally invasive, three-dimensional information on relevant flow conditions in vadose zone sediments as well as changes that occur over time.

Linneman, Dorothy C.↗

DOE OSTI2023

Source Characterization of the Declared North Korean Nuclear Tests From Regional Distance Coda Wave Spectral Ratios

Abstract Seismic observations of underground nuclear explosions provide crucial data on source yield and depth that cannot easily be estimated from other geophysical methods. However, it is difficult to obtain reliable yield estimates for test sites for which we do not have direct seismic calibration experiments. To obtain source information from uncalibrated sites and paths, local and regional seismic records of six, proximal, declared underground nuclear explosions in North Korea are used to compute spectral ratios of narrow‐band waveform envelopes of body‐wave coda that remove path and site effects to reveal precise, relative source moment. The yields of these explosions are obtained from the observed source ratios by simultaneously fitting the classical source model of Mueller and Murphy (1971), https://doi.org/10.1785/bssa0610061675 to all event pairs. The source model provides an impressive fit to the observations considering that the P phase coda derived source spectral ratios did not require a prior knowledge of the source or regionally calibrated corrections to be applied to the data. However, the observed corner frequencies from S wave coda spectral ratios are lower than the source model predictions, but are well fit by models calculated using the corner frequency consistent with the Fisk conjecture. The results presented here provide novel constraints on the spectral distributions of the energy radiated by the sources of the DPRK test series, and allows for an independent evaluation of existing estimated source model parameters.

58 GEOSCIENCES↗

DOE OSTI2021

Rock Valley Accelerated Weight Drop Seismic Data Processing and Picking of P-wave and S-wave Arrival Times

Rock Valley, in the southern end of the Nevada National Security Site, hosts a fault system that was responsible for a shallow (< 3 km below surface ) magnitude 3.7 earthquake in May 1993. In order to better understand this system, seismic properties of the shallow subsurface need to be better constrained. In April and May of 2021, accelerated weight drop (AWD) active-source seismic data were recorded in order to measure P- and S-wave travel-times for the area. This report describes the processing and phase picking of the recorded seismic waveforms. In total, we picked 7,982 P-wave arrivals at offsets up to ~2500 m, and 4,369 S-wave arrivals at offsets up to ~2200 m. These travel-time picks can be inverted for shallow P-wave and S-wave velocity structure in future studies.

58 GEOSCIENCES↗

DOE OSTI2026

Machine learning for seismic low-frequency extrapolation

The cycle-skipping problem that plagues full waveform inversion (FWI) can be at least partially mitigated if low frequencies (which encode the kinematics of wave propagation in seismic data) are recorded. However, seismic sources and receivers are band-limited, so seismic data does not generally include signals down to 0 Hz. To improve our ability to solve the seismic inverse problem, one can synthesize this missing low-frequency (LF) content from the recorded high-frequency (HF) data using machine learning (ML) models. Deep learning models such as convolutional neural networks (CNNs) demonstrate impressive ability to perform low frequency extrapolation. However, such models require powerful hardware (GPU machines) and careful training. We assess the extrapolation capabilities of three different ML models that do not require GPU machines, namely, random forest, Gaussian process regression and gradient boosting, on both synthetic and real data. Experimental results on two synthetic data sets (generated from a low velocity lens embedded in a homogeneous medium, and the Marmousi model) demonstrate that FWI applied to the extrapolated data consistently improves inversion accuracy relative to FWI applied to the original data sets that do not contain low frequencies. Application of low-frequency extrapolation to real data from the Northwest Shelf of Australia demonstrates that tree-based ML models such as gradient boosting can outperform CNNs in terms of both accuracy and computational cost on non-GPU architectures.

58 GEOSCIENCES↗

DOE OSTI2022

Data for Monitoring Fracture Saturation with Internal Transportable Seismic Sources and Twin Neural Networks

Seismic coda-wave analysis is a well-developed method for detecting subtle physical changes in complex media by measuring arrival times in the late-arriving energy from multiply-scattered or reflected waves. However, a challenge arises when multiply-scattered waves are not sufficiently separated in time from the direct arrivals to provide a clear coda wave train. Additional complications for monitoring changes in fracture systems arise when the signals originate from unsynchronized internal sources, such as natural induced seismicity, from acoustic emission, or from transportable intra-fracture sources (chattering dust), that generate uncontrolled signals that vary in arrival time, amplitude and frequency content. Here, we use a twin neural network (also known as a Siamese neural network) for dimensionality reduction to analyze signals from chattering dust to classify the fluid saturation state of a synthetic fracture system. The twin neural network with shared weights generates a low-dimensional representation of the data input by minimizing contrastive loss, serving as the input to a multiclass classifier that accurately classifies whether multiple fractures in a fracture system are fully saturated or partially saturated, or whether a change in saturation has occurred in different fractures in the system. These results show that information buried in unresolved codas from uncontrolled sources can be extracted using machine learning to monitor the evolution of fracture systems caused by physical and chemical processes even when the scattered and direct wave fields overlap.

58 GEOSCIENCES↗

DOE OSTI2023

Examining Infrasound Propagation at High Spatial Resolution Using a Nodal Seismic Array

Infrasound—acoustic waves in the atmosphere below 20 Hz—is a useful monitoring tool. Topography and atmospheric structure strongly control infrasound propagation, and at common source–receiver distances neither of these effects can be ignored when quantitative source constraints are sought. Detailed spatial measurements of the infrasound wavefield would inform propagation models and improve source estimates. However, the “large-N” deployment strategy now well-known in seismology has not yet been realized for infrasound studies. Here, we use the 900-node seismic array from the 2014 Imaging Magma Under St. Helens (iMUSH) experiment as a proxy for a large-N infrasound network, by leveraging acoustic–seismic coupled arrivals. The active-source component of iMUSH consisted of 23 shallowly buried explosions around Mount Saint Helens volcano; these explosions produced epicentral infrasound recorded on the nodes. We find that the bulk presence of ground-coupled infrasound on the nodes is controlled by wind noise and source–receiver distance, with observed arrivals for eight explosions. Explosions with the most extensive coupling produce complex spatial waveform patterns across the array. These patterns are related to both topographic and atmospheric propagation effects, as well as spatially variable site (coupling) effects. We compare our observations to simple topographic diffraction and high-resolution wind advection models, and full-wave numerical simulations. We find strong spatial correlations between (a) coupled arrival strength and modeled topographic obstruction and (b) coupled arrival time and along-path winds. Our seismoacoustic analyses and results are applicable to other existing and future nodal seismic data sets and can expand the utility of such deployments.

58 GEOSCIENCES↗

DOE OSTI2022

Monitoring Fracture Saturation With Internal Seismic Sources and Twin Neural Networks

Seismic coda-wave analysis is a well-developed method for detecting subtle physical changes in complex media by measuring arrival times in the late-arriving energy from multiply scattered or reflected waves. However, a challenge arises when multiply scattered waves are not sufficiently separated in time from the direct arrivals to provide a clear coda wave train. Additional complications for monitoring changes in fracture systems arise when the signals originate from unsynchronized internal sources, such as natural or induced seismicity, from acoustic emission, or from transportable intra-fracture sources (chattering dust), that generate uncontrolled signals that vary in arrival time, amplitude and frequency content. Here, we use a twin neural network (TNN also known as a Siamese neural network) for dimensionality reduction to analyze signals from chattering dust to classify the fluid saturation state of a synthetic fracture system. The TNN with shared weights generates a low-dimensional representation of the data input by minimizing contrastive loss, serving as the input to a multiclass classifier that accurately classifies whether multiple fractures in a fracture system are fully saturated or partially saturated, or whether a change in saturation has occurred in different fractures in the system. Furthermore, these results show that information buried in unresolved codas from uncontrolled sources can be extracted using machine learning to monitor the evolution of fracture systems caused by physical and chemical processes even when the scattered and direct wave fields overlap.

58 GEOSCIENCES↗

DOE OSTI2021

The First Detection of an Earthquake From a Balloon Using Its Acoustic Signature

Extreme temperature and pressure conditions on the surface of Venus present formidable technological challenges against performing ground-based seismology. Efficient coupling between the Venusian atmosphere and the solid planet theoretically allows the study of seismically generated acoustic waves using balloons in the upper atmosphere, where conditions are far more clement. However, earthquake detection from a balloon has never been demonstrated. We present the first detection of an earthquake from a balloon-borne microbarometer near Ridgecrest, CA in July 2019 and include a detailed analysis of the dependence of seismic infrasound, as measured from a balloon on earthquake source parameters, topography, and crustal and atmospheric structure. Our comprehensive analysis of seismo-acoustic phenomenology demonstrates that seismic activity is detectable from a high-altitude platform on Earth, and that Rayleigh wave-induced infrasound can be used to constrain subsurface velocities, paving the way for the detection and characterization of such signals on Venus.

58 GEOSCIENCES↗

DOE OSTI2026

Coupling geophysical, geological, geochemical and mineralogical assessments to examine preferential contaminant transport pathways in interbedded fractured bedrock

This study shows that a multi-faceted approach, combining borehole geophysical logging and surface seismic P-wave first-arrival tomography with confirmatory coring, well installation, and chemical and mineralogical analysis, is effective for identifying difficult-to-locate preferential contaminant transport pathways in deeper fractured bedrock. Seismic tomography detected porous 10–20 m wide elongated fractured conduits that allow acidic groundwater contaminated with uranium (U) and nitrate (NO 3 - ) to migrate within interbedded shale-limestone bedrock over 1000 m from a former disposal facility (S-3 Ponds site) located at the DOE Y-12 National Security Complex in Tennessee (USA). Conventional drilling techniques would easily miss these conduits because they are oriented parallel with fractured bedding planes. Synchrotron analysis of aquifer solids revealed that > 95 % of the U is hexavalent (U VI ). This uranyl (UO 2 2+ ) species is coordinated with carbonate, iron oxide, silicate and phosphate minerals within cemented white to yellow precipitates, which contain U concentrations as high as ∼21.6 % by weight fraction. Identifying the presence of these mineral phases, enables a further understanding of the potential effectiveness of remediation actions. The combination of methodologies presented here can also be applied to other explorations, such as the detection of water supply.

Fate and transport↗

DOE OSTI2020

Beirut Explosion Yield and Mushroom Cloud Height - Effects of the Source Environment

I use crater dimensions to estimate the yield of the August 4th, 2020 Beirut explosion to be equivalent to approximately 1.4 kilotons of TNT with a lower bound of about 0.7 kilotons. Based on the amount of ammonium nitrate reported to have been stored at the Beirut harbor, I assume an upper bound for the yield of 2.75 kilotons. However, it is highly likely that the yield was less than 2.75 kilotons, since reported values for TNT equivalence of ammonium nitrate are typically much less than one hundred percent. The crater-size based yield estimates are based on crater radius estimates from satellite imagery and empirical curves and data for scaled crater radius from past chemical and nuclear explosions. I present evidence that suggests that the relatively large crater radius is due to a high degree of coupling of shock wave energy to the surrounding medium and a reduction of the effective stress because of a high level of saturation of the geologic media beneath the explosion. I provide yield estimates based on seismic body-wave magnitude and crater depth as corroborating evidence. I compare preliminary estimates for the maximum debris cloud height, based on cell phone videos/images, with predicted maximum heights for this yield range from empirical formulas and numerical cloud-rise models. Based on a preliminary analysis of cell phone footage, the observed maximum cloud height appears to be approximately 1600 m. This is much lower than that predicted using standard empirical formulas and buoyant cloud rise models. I present results from a modified buoyant cloud rise model that more accurately predicts the maximum cloud height by allowing for the inclusion of a fixed amount of air and/or water into the fireball at the start of cloud rise. The amount of mass that needs to be added at the start, to reproduce the observed maximum cloud height, is relatively small compared to the total mass entrained during cloud rise. A much greater amount of dry air or debris is required, relative to water, for an equivalent reduction in maximum cloud height. The ammonium nitrate is one possible source for water in the fireball since it was being stored in a very humid environment and ammonium nitrate is known to be hygroscopic. The ground beneath the explosion, especially if it were saturated, and the nearby harbor could also have been sources for water or debris in the fireball.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗