Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ocean model”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

A USCLIVAR Project to Assess and Compare the Responses of Global Climate Models to Drought-Related SST Forcing Patterns: Overview and Results

The USCLI VAR working group on drought recently initiated a series of global climate model simulations forced with idealized SST anomaly patterns, designed to address a number of uncertainties regarding the impact of SST forcing and the role of land-atmosphere feedbacks on regional drought. Specific questions that the runs are designed to address include: What are the mechanisms that maintain drought across the seasonal cycle and from one year to the next? What is the role of the leading patterns of SST variability, and what are the physical mechanisms linking the remote SST forcing to regional drought, including the role of land-atmosphere coupling? The runs were carried out with five different atmospheric general circulation models (AGCM5), and one coupled atmosphere-ocean model in which the model was continuously nudged to the imposed SST forcing. This paper provides an overview of the experiments and some initial results focusing on the responses to the leading patterns of annual mean SST variability consisting of a Pacific El Nino/Southern Oscillation (ENSO)-like pattern, a pattern that resembles the Atlantic Multi-decadal Oscillation (AMO), and a global trend pattern. One of the key findings is that all the AGCMs produce broadly similar (though different in detail) precipitation responses to the Pacific forcing pattern, with a cold Pacific leading to reduced precipitation and a warm Pacific leading to enhanced precipitation over most of the United States. While the response to the Atlantic pattern is less robust, there is general agreement among the models that the largest precipitation response over the U.S. tends to occur when the two oceans have anomalies of opposite sign. That is, a cold Pacific and warm Atlantic tend to produce the largest precipitation reductions, whereas a warm Pacific and cold Atlantic tend to produce the greatest precipitation enhancements. Further analysis of the response over the U.S. to the Pacific forcing highlights a number of noteworthy and to some extent unexpected results. These include a seasonal dependence of the precipitation response that is characterized by signal-to-noise ratios that peak in spring, and surface temperature signal-to-noise ratios that are both lower and show less agreement among the models than those found for the precipitation response. Another interesting result concerns what appears to be a substantially different character in the surface temperature response over the U.S. to the Pacific forcing by the only model examined here that was developed for use in numerical weather prediction. The response to the positive SST trend forcing pattern is an overall surface warming over the world's land areas with substantial regional variations that are in part reproduced in runs forced with a globally uniform SST trend forcing. The precipitation response to the trend forcing is weak in all the models.

Schubert, Siegfried↗
NASA NTRS2020

CMIP6 Historical Simulations (1850‐2014) with GISS‐E2.1

Simulations of the CMIP6 historical period 1850‐‐2014, characterized by the emergence of anthropogenic climate drivers like greenhouse gases, are presented for different configurations of the NASA Goddard Institute for Space Studies (GISS) Earth System ModelE2.1. The GISS‐E2.1 ensembles are more sensitive to forcing than their CMIP5 predecessors (GISS‐E2), but warm less during recent decades due to reduced total forcing. This forcing reduction is attributed to an increase in longwave opacity in pre‐industrial simulations, resulting in an atmosphere less sensitive to further increases in opacity that result from rising greenhouse gas concentrations. This demonstrates the importance of the base climatology to forcing and forced climate trends. Most model versions match observed temperature trends since 1979 from the ocean to the stratosphere. The choice of ocean model is important to the transient climate response, as found in CMIP5 GISS‐E2: the model that more efficiently exports heat to the deep ocean shows a smaller rise in tropospheric temperature. Model sea level rise over the historical period is traced to excessive drawdown of aquifers to meet irrigation demand with a smaller contribution from thermal expansion. This shows how fully coupled models can provide indirect observational constraints upon forcing, in this case, constraining irrigation rates with observed sea level changes. The overall agreement of GISS‐E2.1 with observed trends is familiar from evaluation of its predecessors, as is the conclusion that these trends are almost entirely anthropogenic in origin.

NASA Goddard Institute for Space Studies (GISS) Ea↗
NASA NTRS2020

Computationally Efficient Modeling and Data Assimilation of Near-Surface Variability

Near-surface (< 20m) ocean exhibits high variability due to coupled interactions, for e.g., with the atmosphere, sea ice, land, etc. Here we focus on atmospheric heat and momentum (wind) forcing, which are known to cause diurnal variability within the mixed layer. Only recently with a combination of sufficiently high vertical/horizontal resolution (75L, 1/4deg) and sub-daily atmospheric forcing fields, ocean models are starting to resolve this diurnal variability. However, the computation expense of such a high vertical resolution is burdensome in the context of coupled modeling and data assimilation. An alternative approach is to parameterize this diurnal variability with a prognostic model, that is embedded into the ocean model.In the first part of this presentation, we will demonstrate results with the above two approaches, by comparing them to profiles of near-surface temperature and salinity. In the context of data assimilation and reanalysis, this modeling capability opens the door to re-examine and perhaps improve specification of background (or, ensemble) error characteristics. The second half of this talk will focus on illustrating diurnally varying errors within an ensemble DA, and possible approaches to improve localization (horizontal/vertical) to extract maximum possible observational information content from in-situ and satellite observations of sea surface temperature.

Akella, Santha↗
DOE OSTI2022

Overestimated Eddy Kinetic Energy in the Eddy‐Rich Regions Simulated by Eddy‐Resolving Global Ocean–Sea Ice Models

Abstract The performance of eddy‐resolving global ocean–sea ice models in simulating mesoscale eddies is evaluated using six eddy‐resolving experiments forced by different atmospheric reanalysis products. Interestingly, eddy‐resolving ocean general circulation models (OGCMs) tend to simulate more (less) energetic eddy‐rich (eddy‐poor) regions with a smaller (larger) spatial extent than satellite observation, which finally shows that larger (smaller) mesoscale energy intensity (EI) is simulated in the eddy‐rich (eddy‐poor) regions. Quantitatively, there is an approximately 27%–60% overestimation of EI in the eddy‐rich regions, which are mainly located in the Kuroshio–Oyashio Extension, the Gulf Stream, and the Antarctic Circumpolar Currents regions, although the global mean EI is underestimated by 25%–45%. Apparently, the eddy kinetic energy in the eddy‐poor region is underestimated. Further analyses based on coherent mesoscale eddy properties show that the overestimation in the eddy‐rich regions is mainly attributed to mesoscale eddies’ intensity and is more prominent when mesoscale eddies are in their growth stage.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2004

Research on the Natural Variability of Climate

The scheme for ocean vertical mixing presented in Canuto, Howard, Cheng & Dubovikov and Canuto, Howard, Hogan, Cheng, Dubovikov & Montenegro [GISS mixing scheme] was extended to include the latitudinal dependence reported by Gregg, Sanford & Winkel of the input to the interior ocean turbulence field of energy from internal waves. The resulting latitude dependence was implemented in our module of code for the GISS vertical mixing scheme and tested in the global NCAR-CSM ocean model with 3 degree X 3 degree resolution and 25 levels in the vertical. With the latitude dependence, the background diffusivity decreases from approx. 0.1 sq cm/sec at midlatitudes to approx. less than 0.01 sq cm/sec at the equator. A significant improvement was seen in the tropical ocean, with tropical thermoclines being sharpened in agreement with the requirements of observational studies and the needs of ENS0 modeling. At the same time global measures nf performance such as the meridional overturning and northward heat transport were not adversely affected. Results were presented at the 2004 Layered Ocean Modeling meeting and a paper has been prepared and submitted to Geophysical Research Letters.

Source record↗
NASA NTRSDate not supplied

AMOC Stability in GISS-E2-1-G: A Case of Stochastic Bifurcation

A 10-member ensemble simulation using the NASA GISS-E2-1-G climate model shows a clear bifurcation in the Atlantic Meridional Overturning Circulation (AMOC) strength under the SSP2-4.5 future scenario, leading to 8 stable AMOC-on and 2 stable AMOC-off climate states after 400 years of integration, despite identical forcing of each ensemble member. Ocean models have previously shown a tipping point where sufficient forcing (e.g. by freshwater inputs) can shut down the AMOC. A variety of fully coupled models have demonstrated this, either through hosing experiments or increased precipitation and warming at high latitudes due to increased concentrations of greenhouse gases. In the GISS simulations, there are no external freshwater perturbations. The bifurcation arises freely in the coupled system and is the result of stochastic variability associated with sea-ice transport and melting in the Irminger Sea following a slowing of the GHG concentration increase at the end of the twenty-first century. We believe this is the first time that a CMIP-class model has shown such a bifurcation across an initial condition ensemble.

GISS ModelE2.1↗
NASA NTRS2001

The Low-Frequency Variability of the Tropical Atlantic Ocean

Upper ocean temperature variability in the tropical Atlantic is examined from the Comprehensive Ocean Atmosphere Data Set (COADS) as well as from an ocean model simulation forced by COADS anomalies appended to a monthly climatology. Our findings are as follows: Only the sea surface temperatures (SST) in the northern tropics are driven by heat fluxes, while the southern tropical variability arises from wind driven ocean circulation changes. The subsurface temperatures in the northern and southern tropics are found to have a strong linkage to buoyancy forcing changes in the northern North Atlantic. Evidence for Kelvin-like boundary wave propagation from the high latitudes is presented from the model simulation. This extratropical influence is associated with wintertime North Atlantic Oscillation (NAO) forcing and manifests itself in the northern and southern tropical temperature anomalies of the same sign at depth of 100-200 meters as result of a Rossby wave propagation away from the eastern boundary in the wake of the boundary wave passage. The most apparent association of the southern tropical sea surface temperature anomalies (STA) arises with the anomalous cross-equatorial winds which can be related to both NAO and the remote influence from the Pacific equatorial region. These teleconnections are seasonal so that the NAO impact on the tropical SST is the largest it mid-winter but in spring and early summer the Pacific remote influence competes with NAO. However, NAO appears to have a more substantial role than the Pacific influence at low frequencies during the last 50 years. The dynamic origin of STA is indirectly confirmed from the SST-heat flux relationship using ocean model experiments which remove either anomalous wind stress forcing or atmospheric forcing anomalies contributing to heat exchange.

Haekkinen, Sirpa↗
DOE OSTI2023

Analytical population dynamics underlying harmful algal blooms triggered by prey avoidance

A modified version of the NPZD ecosystem model is used to analytically examine the effects of predation avoidance, a possible mechanism for triggering harmful algal blooms (HAB). To resolve HAB development caused by predation avoidance, an additional phytoplankton functional group is considered, one that has slower nutrient uptake and better predation avoidance characteristics than the non-harmful phytoplankton group used in traditional NPZD models. Because the two phytoplankton groups (one non-harmful and one HAB) compete for only one resource within the same system, steady state (equilibrium) conditions cannot occur without the presence of zooplankton; only the non-harmful phytoplankton group, which defeats the HAB group in the resource competition, can survive in the equilibrium. The presence of sufficient zooplankton effectively acts to replenish the nutrient pool by consuming the non-harmful phytoplankton. When this occurs, two equilibrium states are found: one with both phytoplankton groups coexisting, and one that only includes the HAB group. Further, the condition required for equilibrium is that the total nitrogen within the system should be larger than a threshold determined by model coefficients. The threshold and feasibility of the equilibrium are sensitive to the relative HAB predation avoidance coefficient. If the coefficient is larger than the ratio of net growth rates between the HAB and non-harmful phytoplankton group, the threshold becomes infinite, and an equilibrium is not feasible. The time scale for the system to reach an equilibrium state that includes a HAB group is determined asymptotically. The dependence of a threshold condition as a controlling factor may explain the regime shift of dominant species causing HABs. The ecosystem model is fully implemented into the Regional Ocean Modeling System and applied to an idealized coastal embayment (with depths and geometry taken from San Francisco Bay) to show numerically the dominance of prey avoidance dynamics in a natural shallow water environment that includes advection and diffusion. The analytical results improve strategies for HAB modeling and provide guidance for setting model coefficients necessary to resolve a HAB event.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Understanding AMOC stability: the North Atlantic Hosing Model Intercomparison Project

Abstract. The Atlantic meridional overturning circulation (AMOC) is an important part of our climate system. The AMOC is predicted to weaken under climate change; however, theories suggest that it may have a tipping point beyond which recovery is difficult, hence showing quasi-irreversibility (hysteresis). Although hysteresis has been seen in simple models, it has been difficult to demonstrate in comprehensive global climate models. Here, we outline a set of experiments designed to explore AMOC hysteresis and sensitivity to additional freshwater input as part of the North Atlantic Hosing Model Intercomparison Project (NAHosMIP). These experiments include adding additional freshwater (hosing) for a fixed length of time to examine the rate and mechanisms of AMOC weakening and whether the AMOC subsequently recovers once hosing stops. Initial results are shown from eight climate models participating in the Sixth Coupled Model Intercomparison Project (CMIP6). The AMOC weakens in all models as a result of the freshening, but once the freshening ceases, the AMOC recovers in half of the models, and in the other half it stays in a weakened state. The difference in model behaviour cannot be explained by the ocean model resolution or type nor by details of subgrid-scale parameterisations. Likewise, it cannot be explained by previously proposed properties of the mean climate state such as the strength of the salinity advection feedback. Instead, the AMOC recovery is determined by the climate state reached when hosing stops, with those experiments where the AMOC is weakest not experiencing a recovery.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2020

Decay of the Snow Cover Over Arctic Sea Ice From ICESat‐2 Acquisitions During Summer Melt in 2019

From the onset of melt in early June, corresponding declines in Ice, Cloud, and Land Elevation Satellite‐2 (ICESat‐2) freeboard and surface albedo can be seen over the entire Arctic sea ice cover. In the 2019 summer, area‐averaged freeboard decreased from 34 cm prior to melt to a minimum of 12 cm in August while the area‐averaged albedo decreased from ~0.7 to 0.38 for the same period. Calculations using ICESat‐2 freeboards and modeled ice thickness from Pan‐Arctic Ice Ocean Modeling and Assimilation System (PIOMAS) give area‐averaged snow depths ranging from 17 cm prior to melt to 3 cm in August over seasonal ice and from 34 to 4 cm over multiyear ice. Mean rates of snow ablation (including evaporation) in mid‐June were as high as 2 cm/day, comparable to field records from other years. Increases in freeboard after mid‐August in the high latitude (>80°N) multiyear ice cover, north of the Greenland coast, are likely due to earlier freeze‐up and snow accumulation in these regions with shorter melt seasons.

R Kwok↗
NASA NTRS1999

Arctic Summer Ice Processes

The primary objective of this study is to estimate the flux of heat and freshwater resulting from sea ice melt in the polar seas. The approach taken is to examine the decay of sea ice in the summer months primarily through the use of spaceborne Synthetic Aperture Radar (SAR) imagery. The improved understanding of the dynamics of the melt process can be usefully combined with ice thermodynamic and upper ocean models to form more complete models of ice melt. Models indicate that more heat is absorbed in the upper ocean when the ice cover is composed of smaller rather than larger floes and when there is more open water. Over the course of the summer, floes disintegrate by physical forcing and heating, melting into smaller and smaller sizes. By measuring the change in distribution of floes together with open water over a summer period, we can make estimates of the amount of heating by region and time. In a climatic sense, these studies are intended to improve the understanding of the Arctic heat budget which can then be eventually incorporated into improved global climate models. This work has two focus areas. The first is examining the detailed effect of storms on floe size and open water. A strong Arctic low pressure storm has been shown to loosen up the pack ice, increase the open water concentration well into the pack ice, and change the distribution of floes toward fewer and smaller floes. This suggests episodic melting and the increased importance of horizontal (lateral) melt during storms. The second focus area is related to an extensive ship-based experiment that recently took place in the Arctic called Surface Heat Budget of the Arctic (SHEBA). An icebreaker was placed purposely into the older pack ice north of Alaska in September 1997. The ship served as the base for experimenters who deployed extensive instrumentation to measure the atmosphere, ocean, and ice during a one-year period. My experiment will be to derive similar measurements (floe size, open water, temporal change) using spaceborne SAR data obtained during the summer of 1998, and compare these results with an ocean and ice model of summer melt. Additional information is contained in the original.

Holt, Benjamin↗
NASA NTRS1999

Sensitivity of the Tropical Pacific Ocean to Precipitation Induced Freshwater Flux

We have performed a series of experiments using an ocean model to study the sensitivity of tropical Pacific Ocean to variations in precipitation induced freshwater fluxes. Variations in these fluxes arise from natural causes on all time scales. In addition, estimates of these fluxes are uncertain because of differences among measurement techniques. The model used is a quasi-isopycnal model, covering the Pacific from 40 S to 40 N. The surface forcing is constructed from observed wind stress, evaporation, precipitation, and surface temperature (SST) fields. The heat flux is produced with an iterative technique so as to maintain the model close to the observed climatology, but with only a weak damping to that climatology. Climatological estimates of evaporation are combined with various estimates of precipitation to determine the net surface freshwater flux. Results indicate that increased freshwater input decreases salinity as expected, but increases temperatures in the upper ocean. Using the freshwater flux estimated from the Microwave Sounding Unit leads to a warming of up to 0.6 C in the western Pacific over a case with zero net freshwater flux. SST is sensitive to the discrepancies among different precipitation observations, with root-mean-square differences in SST on the order of 0.2-0.3 C. The change in SST is more pronounced in the eastern Pacific, with differences of over 1 C found among the various precipitation products. Interannual variation in precipitation during El Nino events leads to increased warming. During the winter of 1982-83, freshwater flux accounts for about 0.4 C (approximately 10-15% of the maximum warming) of the surface warming in the central-eastern Pacific. Thus, the error of SST caused by the discrepancies in precipitation products is more than half of the SST anomaly produced by the interannual variability of observed precipitation. Further experiments, in which freshwater flux anomalies are imposed in the western, central, and eastern Pacific, show that the influence of net freshwater flux is also spatially dependent. The imposition of freshwater flux in the far western Pacific leads to a trapping of salinity anomaly to the surface layers near the equator. An identical flux imposed in the central Pacific produces deeper and off-equatorial salinity anomalies. The contrast between these two simulations is consistent with other simulations of the western Pacific barrier layer information.

Yang, Song↗
DOE OSTI2024

Estuarine hurricane wind can intensify surge-dominated extreme water level in shallow and converging coastal systems

Based on the projected increase in hurricane landfall frequency on the middle to lower US east coast, we examined the crucial role of the estuarine wind field in exacerbating coastal flooding. A regionally refined atmospheric and two high-resolution hydrology and ocean models are integrated to provide plausible and physically consistent ensembles of hurricane events and the associated flooding inside the Delaware Bay and River, a US mid-Atlantic estuary. Model results show that the hurricane propagation direction, estuarine geometry, remote surge from the open ocean, and direct nearshore upwind stress could magnify the flood magnitude. More specifically, inland-bound tracks that make landfall before reaching the mid-Atlantic coast produce a more significant surge within Delaware Bay than the shore-parallel tracks, where the estuarine wind direction plays the primary role in surge amplification. Ultimately, this study emphasized the need for integrated models to capture the nonlinear dynamics and interactions in flood hazard modeling.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1979

A large-scale numerical model of sea ice

The described large-scale sea ice model, which is capable of coupling with atmospheric and oceanic models of comparable resolution, simulates the yearly cycle of ice in both the Northern and Southern Hemispheres. Model results for the yearly cycle of sea ice thickness and extent in both the Arctic and the Antarctic are presented. Horizontally the model resolution is approximately 200 km, while vertically four layers - ice, snow, ocean, and atmosphere - are considered. Thermodynamic processes based on energy balances at the various interfaces and dynamic processes based on wind stress, water stress, Coriolis force, internal ice resistance, and the stress from the tilt of the sea surface are incorporated. It is assumed that the ice within a given grid square is of uniform thickness, although each square has a variable percentage of its area ice free.

Parkinson, C. L.↗
DOE OSTI2020

River Plume Liftoff Dynamics and Surface Expressions

The water surface expression of liftoff and its dependence on discharge are examined using numerical simulations with the Regional Ocean Modeling System (ROMS). Liftoff is the process by which buoyant river water separates from the bed and flows over denser saltwater. During low-discharge conditions liftoff occurs in the river and is accompanied by a change in the surface slope. During high-discharge conditions liftoff occurs outside the mouth and generates a ridge on the water surface. The location and height of the ridge can be described by analytical equations in terms of discharge, shelf slope, and river mouth aspect ratio. The offshore distance and height of the ridge are proportional to the river discharge and vary inversely with river mouth aspect ratio. For steep shelf slopes liftoff occurs close to the river mouth and generates a large ridge. The ridge is modified, but not eliminated, by the presence of tides. The water surface slope change at the ridge peak is large enough to be detected by the upcoming Surface Water and Ocean Topography (SWOT) altimeter and can be used to identify the liftoff location during high discharge. However, during low discharge the water surface slope change at the liftoff location is too small to be detected by SWOT. These results indicate that remote measurements of the presence or absence of the ridge may be useful to distinguish between low and high flows, and remote measurements of the ridge location or height could be used to estimate freshwater discharge.

42 ENGINEERING↗
NASA NTRS2000

A Coupled Ocean General Circulation, Biogeochemical, and Radiative Model of the Global Oceans: Seasonal Distributions of Ocean Chlorophyll and Nutrients

A coupled ocean general circulation, biogeochemical, and radiative model was constructed to evaluate and understand the nature of seasonal variability of chlorophyll and nutrients in the global oceans. Biogeochemical processes in the model are determined from the influences of circulation and turbulence dynamics, irradiance availability. and the interactions among three functional phytoplankton groups (diatoms. chlorophytes, and picoplankton) and three nutrients (nitrate, ammonium, and silicate). Basin scale (greater than 1000 km) model chlorophyll results are in overall agreement with CZCS pigments in many global regions. Seasonal variability observed in the CZCS is also represented in the model. Synoptic scale (100-1000 km) comparisons of imagery are generally in conformance although occasional departures are apparent. Model nitrate distributions agree with in situ data, including seasonal dynamics, except for the equatorial Atlantic. The overall agreement of the model with satellite and in situ data sources indicates that the model dynamics offer a reasonably realistic simulation of phytoplankton and nutrient dynamics on synoptic scales. This is especially true given that initial conditions are homogenous chlorophyll fields. The success of the model in producing a reasonable representation of chlorophyll and nutrient distributions and seasonal variability in the global oceans is attributed to the application of a generalized, processes-driven approach as opposed to regional parameterization and the existence of multiple phytoplankton groups with different physiological and physical properties. These factors enable the model to simultaneously represent many aspects of the great diversity of physical, biological, chemical, and radiative environments encountered in the global oceans.

Gregg, Watson W.↗
NASA NTRSDate not supplied

How strong is ocean mixing in the Pacific Ocean cold tongue, and does it matter?

Microstructure observations within the Pacific Ocean cold tongue reveal that strong surface boundary layer turbulence often penetrates into the thermocline producing hundreds of W/m^2 of downward heat transport during nighttime and early morning. Hence, it has been suggested that ocean turbulent mixing contributes significantly to the sea-surface temperature budget and modulates air-sea interaction on a wide range of timescales with consequences for global climate and weather. However, most observations of ocean microstructure in the cold tongue are from 0 N,140 W. In this presentation, I will attempt to characterize ocean mixing at and beyond 0 N,140 W using models. Questions to be addressed include: 1) can we model ocean mixing at 0 N,140 W? And 2) how does ocean mixing vary in space and time beyond 0 N,140 W? Finally, 3) how does ocean mixing compare to the air-sea heat flux in the sea-surface temperature budget? Results will be derived from existing observations at 0 N, 140 W and a new hierarchy of regional ocean process simulations including a 20-year-long submesoscale-permitting regional model and month-long turbulence-permitting large eddy simulations embedded in a regional model.

Strong↗
DOE OSTI2024

Assessing Effects of Climate Change on Legacy Waste at the Enewetak Atoll

The Republic of the Marshall Islands (RMI) is in the central Pacific Ocean ~4,500 km west of Hawaii. The Enewetak Atoll, located in the northwest part of the RMI, was the site for 43 nuclear weapon tests between 1948 and 1958. Fallout and deposition from the tests contaminated the island surfaces, lagoon waters and sediment, and nearby ocean waters at the atoll. In the 1970s, a cleanup effort collected radioactive waste and placed it in the Cactus Crater on Runit Island (also called the Runit Dome). In December 2021, Congress directed the U.S. Department of Energy to study the impacts of climate change on the Runit Dome nuclear waste disposal site. Pacific Northwest National Laboratory (PNNL) assembled a multidisciplinary team of climate scientists, ocean modelers, environmental scientists, and health physicists to assess the likely effects of remaining radionuclides at the Enewetak Atoll. PNNL’s approach focused on effects of tropical cyclones that were postulated to mobilize and transport contaminated lagoon sediments and result in human and biota exposure. PNNL’s study estimated (1) the radionuclide source term, (2) the effects of climate change on severe storms, (3) mobilization and transport of radionuclides, and (4) radiation dose to humans and biota. Radionuclides in the lagoon and/or ocean waters of the Enewetak Atoll were characterized by the U.S. Atomic Energy Commission (AEC) in 1972, Woods Hole Oceanographic Institution in 2015, and Lawrence Livermore National Laboratory in 2018. The RMI Nationwide Radiological Study was conducted in the early 1990s for radionuclides remaining in island soils. The 1972 AEC survey remains the most comprehensive source of radionuclide data on lagoon sediments. Climate change modeling at a regional scale in the central Pacific Ocean is limited. PNNL climate scientists simulated severe historical storms postulated to occur both in a recent climate (2015) and in the future (2090) using the Advanced Research Weather Research and Forecasting (WRF-ARW) model, employing a pseudo-global-warming technique. A postulated complete, future failure of the Runit Dome was also considered. PNNL developed a high-resolution regional ocean hydrodynamics model covering the entire RMI extended economic zone using the Finite Volume Coastal Ocean Model (FVCOM). The FVCOM model was run using global reanalysis data for current climate and WRF-ARW simulation for the future climate. PNNL also developed a radionuclide fate and transport model using the FVCOM Integrated Compartment Model (FVCOM-ICM) to simulate the current and future mobilization and transport of radionuclides sorbed to lagoon sediments and the exchange of radionuclides between the water and sediment. FVCOM-ICM-predicted radionuclide concentrations were then used to estimate radiation dose to humans and biota at all islands of the Enewetak Atoll. Under current climate conditions, annual radiation exposures for the southern islands including Enewetak (Fred) and Medren (Elmer) were below the current U.S. standards. Radiation doses were somewhat elevated starting at Runit Island northward and westward to Enjebi Island (Janet). The islands in the northwest quadrant, particularly Bokoluo (Alice) and Bokombako (Belle), remain relatively contaminated. The islands in the southwestern quadrant have low contamination. The highest contribution to radiation doses comes from consumption of locally grown foods. Two radionuclides, 90Sr and 137Cs, contributed the greatest fraction for most terrestrial foods. In current climate conditions, the storms temporarily increased radionuclide concentrations in the lagoon waters, increasing the radiation dose slightly. In future conditions, doses are expected to be smaller, primarily because of the radioactive decay of the shorter-lived radioisotopes of 90Sr and 137Cs. This could make all islands in the far northwest of the atoll – except Bokombako (Belle) and perhaps Bokoluo (Alice) – suitable for residency. For the f

Prasad, Rajiv↗