The Effect of Earthquake Kinematics on Tsunami Propagation
Explore the source record for details and available documents.
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
This Web-based tsunami simulation system allows users to remotely run a model on JPL s supercomputers for a given undersea earthquake. At the time of this reporting, predicting tsunamis on the Internet has never happened before. This new code directly couples the earthquake model and the ocean model on parallel computers and improves simulation speed. Seismometers can only detect information from earthquakes; they cannot detect whether or not a tsunami may occur as a result of the earthquake. When earthquake-tsunami models are coupled with the improved computational speed of modern, high-performance computers and constrained by remotely sensed data, they are able to provide early warnings for those coastal regions at risk. The software is capable of testing NASA s satellite observations of tsunamis. It has been successfully tested for several historical tsunamis, has passed all alpha and beta testing, and is well documented for users.
A collaborative effort between Pacific Northwest National Laboratory (PNNL) and Oak Ridge National Laboratory (ORNL) is underway to provide a technical basis and methodology for the criticality safety community to use the sum of fractions (SoF) method for generating limits for mixtures of “selected actinide nuclides” included in the ANSI/ANS 8.15 standard. The PNNL scope in this project is to define a range of mixtures of 233 U, 235 U, and 239 Pu moderated with either light water or polyethylene and to examine the critical masses for these mixtures. The ORNL scope is primarily to provide validation support for these studies. More complete discussion of the project and its validation aspects will be presented at the upcoming International Conference on Nuclear Criticality Safety (ICNC) this Fall in Sendai, Japan. Clear differences in benchmark similarity to application systems as assessed by the integral parameter ck are noted in the validation studies performed as part of this project as a function of moderator. The c k value is a correlation coefficient that represents that amount of shared uncertainty in k eff due to cross sections between two systems. Individual nuclide-reaction contributions between the two systems can be simply summed to arrive at the total c k value. Specifically, the c k values for light-water–moderated solution experiments are higher for a water-moderated application than for a polyethylene-moderated application. This result is neither totally unexpected nor surprising, but the magnitude of the difference was difficult to anticipate. The TSUNAMI sequence, in the SCALE 6.2.4 code package developed by ORNL, was used to generate eigenvalues and reactivity effects with perturbation-theory based approach through sensitivity coefficients for all nuclides in the system with all reactions and energy groups. The TSUNAMI-Indices and Parameters (IP) sequence then uses the sensitivity data generated through TSUNAMI to generate relational parameters (i.e., c k ) to determine the degree of similarity between systems. One detail of the SCALE material and data implementation must be discussed at this point. Several thermal scattering laws (TSLs) are available for 1 H. SCALE uses a different nuclide ID number for each TSL; essentially, each version of 1 H is treated as a unique nuclide. For example, 1 H bound in water ( 1 H-H2O) is assigned the nuclide ID 1001, whereas 1 H bound in polyethylene (h-poly) is assigned the nuclide ID 9001001. The same cross section data are used for all reactions in 1 H, regardless of TSL, except for scattering below the TSL cutoff energy. TSUNAMI-IP treats different nuclide IDs as different nuclides; thus, no uncertainty is shared between 1 H-H2O and h-poly, despite much of the same data, including covariance data, being used for both nuclides. This presents a question: how much of the difference in assessed similarity between water- and polyethylene-moderated systems is due to the moderators, and how much is caused by the treatment of 1 H-H2O and h poly with cross section and covariance data. The extended edits generated by TSUANMI-IP allow for an investigation of this issue specifically, as well as a demonstration of the general techniques available within TSUNAMI to understand the results of the similarity assessment. This paper presents and analyzes the similarity assessment of both water- and polyethylene-moderated systems for a single benchmark: PST-002-001.
The fragmentation of a small asteroid in the atmosphere greatly increases its cross sections for aerodynamic braking and energy dissipation. At a typical impact velocity of 22 km/s, the atmosphere absorbs more than half the kinetic energy of stony meteoroids with diameters, D(sub m), less than 220 m and iron meteoroids with D(sub m) less than 80 m. The corresponding diameter for comets with impact velocity 50 km/s is D(sub m) less than 1600 m. Most of the atmospheric energy dissipation occurs in a fraction of a scale height, so large meteors appear to 'explode' or 'flare' at the end of their visible paths. This dissipation of energy in the atmosphere protects the earth from direct impact damage (e.g., craters), but it produces a blast wave that can do considerable damage. The area of destruction around the impact point in which the over-pressure in the blast wave exceeds 4 lb/sq in = 2.8 x 10(exp 5) dynes/cu cm, which is enough to knock over trees and destroy buildings, increases rapidly from zero for chondritic meteoroids less than 56 m in diameter (15 megatons) to about 200 sq km for those 80 m in diameter (48 megatons); the probable diameter of the tunguska impactor of 1908 is about 80 m. Crater formation and earthquakes are not significant in land impacts by stony asteroids less than about 200 m in diameter because of the air protection. A tsunami is probably the most devastating type of damage for asteroids 200 m to 1 km in diameter. An impact by an asteroid this size anywhere in the Atlantic would devastate coastal areas on both sides of the ocean. An asteroid a few kilometers across would produce a tsunami that would reach the foothills of the Appalachian Mountains in the upper half of the East Coast of the United States. Most of Florida is protected from a tsunami by the gradual slope of the ocean off its coast, which causes most of the tsunami energy to be reflected back into the Atlantic. The atmosphere plume produced by asteroids with diameters exceeding about 120 m cannot be contained by the atmosphere, so this bubble of high-temperature gas forms a new layer on top of the atmosphere. The dust entrapped in this hot gas is likely to have optical depths exceeding tau = 10 for asteroids with diameters exceeding about 0.5 to 1 km. The optical flux from asteroids 60 m or more in diameter is enough to ignite pine forests. However, the blast wave from an impacting asteroid goes beyond the radius in which the fire starts. The blast wave tends to blow out the fire, so it is likely that the impact will char the forest, as at Tunguska, but the impact will not produce a sustained fire. Because comets dissipate their energy much higher in the atmosphere than asteroids, they illuminate a much larger region and their blast wave is weaker. So they are much more effective in producing large fires. This suggests that the KT impactor was a comet rather than an asteroid.
This guide was prepared for the field trip to the KT elastic sequence of northeastern Mexico, 5-8 February 1994, in conjunction with the Conference on New Developments Regarding the KT Event and Other Catastrophes in Earth History, held in Houston, Texas. The four-day excursion offers an invaluable opportunity to visit three key outcrops: Arroyo El Mimbral, La Lajilla, and El Pinon. These and other outcrops of this sequence have recently been interpreted as tsunami deposits produced by the meteorite impact event that produced the 200 to 300-km Chicxulub basin in Yucatan, and distributed ejecta around the world approximately 65 m.y. ago that today is recorded as a thin clay layer found at the K/T boundary. The impact tsunami interpretation for these rocks has not gone unchallenged, and others examining the outcrops arrive at quite different conclusions: not tsunami deposits but turbidites; not KT at all but 'upper Cretaceous.' Indeed, it is in hopes of resolving this debate through field discussion, outcrop evaluation, and sampling that led the organizers of the conference to sanction this field trip. This field guide provides participants with background information on the KT clastic sequence outcrops and is divided into two sections. The first section provides regional and logistical context for the outcrops and a description of the clastic sequence. The second section presents three representative interpretations of the outcrops by their advocates. There is clearly no way that these models can be reconciled and so two, if not all three, must be fundamentally wrong. Readers of this guide should keep in mind that many basic outcrop observations that these models are based upon remain unresolved. While great measures were taken to ensure that the information in the description section was as objective as possible, many observations are rooted in interpretations and the emphasis placed on certain observations depends to some degree upon the perspective of the author.
The TSUNAMI suite within the SCALE code package includes several methods for generating sensitivity data, including multigroup (MG) and continuous-energy (CE) capabilities. For generating sensitivities with CE data, three methods are available in SCALE 6.3.0: (1) the iterated fission probability (IFP) method with the KENO Monte Carlo transport solver, (2) IFP with the Shift Monte Carlo transport solver, and (3) the Contributon-Linked eigenvalue sensitivity/Uncertainty estimation via Tracklength importance Characterization (CLUTCH) with the KENO Monte Carlo transport solver. Currently, it is difficult to generate accurate sensitivities with large reflectors when using the CLUTCH method, specifically with fissionable and hydrogenous materials. To address this issue, the work presented herein examines a methodology to calculate the adjoint flux externally with the 3D deterministic SN transport code DENOVO in SCALE; the result is then read directly into the CLUTCH-TSUNAMI sequence. This hybridization method replaces the Monte Carlo F*(r) calculation in CLUTCH while still utilizing the forward calculation. The critical benchmark HEU-MET-FAST-028-001 is used to generate sensitivities based on the inability of CLUTCH to generate accurate sensitivities. Results from the hybrid method appear to generate sensitivity values that are in excellent agreement with direct perturbations. Although further testing is needed, the method provides promising results for the development and utility of a hybrid method for use in TSUNAMI.
This study addresses the need for comprehensive investigations into TRi-structural ISOtropic (TRISO) fuel pebble transportation validation. In this work, an exploratory model, the pebble tanker(PT), was developed with the aim of facilitating the validation of nuclear criticality safety calculations in the context of industrial-scale transportation of TRISO fuel. The PT model was designed to investigate the availability and applicability of critical benchmark experiments crucial for assessing the transportation of these pebbles. This work incorporated sensitivity/uncertainty (S/U) similarity studies to quantify the applicability of critical benchmark experiments and to address nuclear data uncertainties in the context of TRISO transportation. Two container models were investigated: one for the Hermes-type pebble and one for the Pebble Bed Modular Reactor (PBMR)–type pebble. The models were simplified, considering fuel, containment, and either water or air, to enable a focus on the underlying physics of applications involving TRISO fuel pebbles using the PT model. A crucial aspect under consideration was the capacity of the transport package to hold pebbles while ensuring subcriticality in the flooded state. An approach in the criticality validation process involves assessing the similarity between systems through an integral index parameter evaluation. This involves calculating a correlation coefficient (referred to as c k ) based on shared nuclear data–induced uncertainty between a benchmark experiment and the application of the PT model. To facilitate this analysis, the SCALE tools, particularly the CSAS6-Shift, TSUNAMI-3D-Shift, and TSUNAMI-IP sequences, were employed for comprehensive studies in neutronics and S/U analysis. Our findings showed that there are sufficient critical experimental benchmarks to perform this validation of the PT model in the most reactive state, i.e. when the tanker is flooded. This paper provides valuable insights into validating a transport package for Generation IV TRISO fuel pebbles.
This presentation is on the Oak Ridge National Labatory (ORNL) S/U capabilities for critical experiment design. It covers a brief intro to TSUNAMI analysis for NCS validation. Additionally, it shows similarity indices and historical use of TSUNAMI for critical experiment design. It concludes with recent designs: IER-441, IER-304, and IER-554 (on-going).
Rapid and reliable estimation of large earthquake magnitude (above 8) is key to mitigating the risks associated with strong shaking and tsunamis. Standard early warning systems based on seismic waves fail to rapidly estimate the size of such large earthquakes. Geodesy-based approaches provide better estimations, but are also subject to large uncertainties and latency associated with the slowness of seismic waves. Recently discovered speed-of-light prompt elastogravity signals (PEGS) have raised hopes that these limitations may be overcome, but have not been tested for operational early warning. Here we show that PEGS can be used in real time to track earthquake growth instantaneously after the event reaches a certain magnitude. We develop a deep learning model that leverages the information carried by PEGS recorded by regional broadband seismometers in Japan before the arrival of seismic waves. After training on a database of synthetic waveforms augmented with empirical noise, we show that the algorithm can instantaneously track an earthquake source time function on real data. Our model unlocks ‘true real-time’ access to the rupture evolution of large earthquakes using a portion of seismograms that is routinely treated as noise, and can be immediately transformative for tsunami early warning.
The SCALE code package offers multiple nuclear data libraries and sensitivity and uncertainty (S/U) methods supporting and derived from Monte Carlo (MC) transport. The CSAS and TSUNAMI-3D sequences use KENO MC, utilizing either continuous-energy (CE) cross sections or multigroup (MG) cross section libraries. TSUNAMI-3D has two CE calculational methods: the iterated fission probability (IFP) method, and the Contribution-Linked eigenvalue sensitivity/Uncertainty estimation via Tracklength importance CHaracterization (CLUTCH) method. Previous work has shown poor agreement between CLUTCH and confirmatory direct perturbation calculations in specific applications (e.g., fissionable and polyethylene reflectors). The HEU-MET-FAST-084 (HMF-084) International Criticality Safety Benchmark Evaluation Project evaluation consists of 27 cylindrical highly enriched uranium metal cores with 14 unique reflector materials of 0.5 and 1 in. thicknesses. Included in this list of reflector materials are natural uranium and polyethylene. This work utilized SCALE 6.2.4 models of the HMF-084 evaluation, with additional non-physical configurations to test both the MG bias and CLUTCH functionality across a variety of reflector material thicknesses. The evaluation’s use of concentric cylinders allowed for examination of several MG self-shielding methods: infinite homogenous, cylindrical, and spherical. The results indicate that the use of polyethylene reflectors with CLUTCH is not fundamentally impossible but sensitive to geometry. The poor performance of CLUTCH with fissionable reflectors was reaffirmed. The 2 in. and greater polyethylene-reflected calculations demonstrate the necessity of using the 302-group library for fast systems. The nickel MG bias was substantial, as discussed in a companion paper, as were cobalt and iron.
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.”
A joint effort by PNNL and ORNL examined mixtures of 233 U, 235 U, and 239 Pu to determine subcritical mass limits based on similarities to critical benchmark experiments through the Sum-of-Fractions method. Results of the validation efforts by ORNL were used to determine areas of applicability based on EALF and moderator/reflector type. TSUNAMI-1D, TSUNAMI-IP, and VADER sequences used for generating sensitivities, similarity assessments, and statistical testing. Validation efforts provided here using the Sum-of-Fractions method results in k eff values below established calculational margins and provides additional flexibility to limits provided in ANSI/ANS-8.15.
This report documents a study performed to investigate the requirements for criticality safety benchmark experiments for high-assay, low-enriched uranium (HALEU) fuel in transportation applications. In this work, an exploratory application model, the “Pebble Tanker,” was developed to represent TRISO fuel in a transportation scenario for an analysis of the validation basis in industrial quantities. An aspect of the criticality validation process involves assessing the “similarity” between application and experimental benchmark systems through an integral index parameter evaluation. Here, this includes propagating nuclear data uncertainties and calculating a correlation coefficient (hereinafter referred to as “c k ”) to evaluate the similarity of benchmark experiments compared with the application Pebble Tanker model. Finding sufficient critical benchmark experiments allows for the evaluation of bias and bias uncertainty, thus determining the upper subcritical limit (USL) of the transportation package. A target k eff of ~0.94 was used in this work to establish appropriate modeling conditions, reflecting a reasonable estimate for a USL. Two container models were investigated: one with the Hermes-type pebble and one with the Pebble Bed Modular Reactor (PBMR)–type pebble. The models were simplified, considering only fuel, containment structure, and either water or air. This allows a focus on the underlying physics of applications involving TRISO fuel pebbles using the Pebble Tanker model. A crucial consideration is the transport package's ability to safely hold pebbles while flooded, maintaining subcritical conditions. Tools available in the SCALE 6.3.1 suite—the CSAS6-Shift, TSUNAMI-3D-Shift, and TSUNAMI-IP sequences—were employed for neutronics and sensitivity and uncertainty (S/U) analysis of the Pebble Tanker. Findings demonstrated sufficient available critical experiment benchmarks to perform a validation of the Pebble Tanker in the most reactive state, i.e., when the Tanker is flooded.
This report describes an assessment of sodium fast reactor assemblies loaded within an AREVA MAP package. The assessment was performed to determine the state of the validation basis for sodium fast reactor fuel within a transportation package originally intended for fresh, light-water reactor fuel assemblies. A similar study is being simultaneously released (Cumberland, 2025), and substantial parts of the explanatory text are identical to that parallel work, which follows the same workflow. Both studies employed the TSUNAMI-3D and TSUNAMI-IP sequences from the SCALE code system to determine correlation coefficients between various configurations and databases of validation assessments.
This report describes an assessment of sodium fast reactor assemblies loaded inside a Transnuclear B1 package. The assessment was performed to determine the state of the validation basis for sodium fast reactor fuel within a transportation package originally intended for fresh, light-water reactor fuel assemblies. The study employed the TSUNAMI-3D and TSUNAMI-IP sequences from the SCALE code system to determine correlation coefficients between various configurations and databases of validation assessments.