A Comparative Study Between Regional Atmospheric Model Simulations Coupled and Uncoupled to a Regional Ocean Model of the Indian Summer Monsoon
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.
Tropical cyclones (TCs) alter upper-ocean temperature and influence ocean heat content via enhanced turbulent mixing. A better understanding of the role of TCs within the climate system requires a fully coupled modeling framework, where TC-induced ocean responses feed back to the atmosphere and subsequently to the climate mean state and variability. Here, we investigate the impacts of TC wind forcing on the global ocean and the associated feedbacks within the climate system using the fully coupled Community Earth System Model version 1.3 (CESM1.3). Using the low-resolution version of CESM1.3 (1° atmosphere and ocean grid spacing) with no intrinsic TCs, we conduct a suite of sensitivity experiments by inserting TC winds extracted from a high-resolution (0.25° atmosphere grid spacing) TC-permitting simulation into the low-resolution model. Results from the low-resolution TC experiment are compared to a low-resolution control simulation to diagnose TCs’ impact. We found that the added TC winds can increase ocean heat content by affecting ocean vertical mixing, air–sea enthalpy fluxes, and cloud amount. Here, the added TCs can influence mean SST, precipitation, ocean subsurface temperature, and ocean mixed layer depth. We found a strengthening of the wind-driven subtropical cells and a weakening of the Atlantic meridional overturning circulation due to the changes of surface buoyancy fluxes. TCs in the model cause anomalous equatorward ocean heat convergence in the deep tropics and an increase of poleward ocean heat transport out of the subtropics. Our modeling results provide new insights into the multiscale interactions between TCs and the coupled climate system.
A coupling toolkit has been developed to reduce the complexity of model component coupling, in particular among hierarchical model components. The toolkit provides the services of data registration, data matching, data filtering, and model component coupling. In addition, it can generate diagrams to reveal the "producer"-to-"consumer" relations among the components. We have tested this toolkit with the operational NASA Goddard Earth Observing System (GEOS-5), which is built on the Earth System Modeling Framework (ESMF) and consists of several-level Earth system components.
We report given the current lack of long-term disposal for spent nuclear fuel storage in the United States, the current inventory of advanced test reactor (ATR) fuel at Idaho National Laboratory (INL) is expected to be stored in a dry storage facility for an extended storage period. As part of an effort to understand the conditions during a 50-year period, a three-dimensional computational fluid dynamics (CFD) of the unsealed canisters packed with spent fuel in the storage facility was constructed, this CFD model is coupled with radiolytic chemical reactions associated with water vapor and other species. The coupled model is simulated over a range of sensitivity parameters for the fuel decay heat, relative humidity, and oxyhydroxide film thickness for undried and fully-dried fuel. The most important parameter was the initial fuel decay heat. Over the range of cases, the hydrogen and nitric acid concentrations are low enough no significant impact is expected.
An effort has been dedicated to developing a reference model for multi-physics coupled simulation of the pebble bed high-temperature gas-cooled reactor (PB-HTGR) with SAM and Griffin computer codes for safety analysis and licensing purpose. The reference problem of a prototypical reactor design serves as the foundation for the U.S. NRC (Nuclear Regulatory Commission) to verify the adequacy of computer codes and evaluation models for specific reactor types. In this work, a SAM model of the HTR-PM reactor has been developed based on publicly available design information and the multi-dimensional Pronghorn model developed by Idaho National Laboratory. The SAM HTR-PM model includes a multi-dimensional core region and 0- D/1-D fluid components. Additionally, a simplified air RCCS loop is modeled for decay heat removal. The Griffin model is based on a recent work by Idaho National Laboratory. The coupling between the models is achieved through the MOOSE MultiApp system. Both steady-state and transient scenarios were simulated to demonstrate the coupled model’s capability for multi-physics simulations.
Compound riverine and coastal flooding is usually driven by complex interactions among meteorological, hydrological, and ocean extremes. However, existing efforts to model this phenomenon often do not integrate hydrological processes across atmosphere–land–river–ocean systems, leading to substantial uncertainties that have not been fully examined. To bridge this gap, we leverage the new capabilities of the Energy Exascale Earth System Model (E3SM) that enable a multi-component framework that integrates coastally refined atmospheric, terrestrial, and oceanic components. We evaluate compound uncertainties arising from two-way land–river–ocean coupling in E3SM and track the cascading meteorological and hydrological uncertainties through ensemble simulations over the Delaware River basin and estuary during Hurricane Irene (2011). Our findings highlight the importance of two-way river–ocean coupling to compound flood modeling and demonstrate E3SM's capability in capturing compound flood extent near the coast, with a hit rate over 0.75. Our study shows the growing uncertainties that transition from atmospheric forcings to flood distribution and severity. Furthermore, an analysis based on artificial neural networks is used to assess the roles of hydrological drivers, such as infiltration and soil moisture, in the generation of compound flooding. The response of compound floods to tropical cyclones (TCs) is found to be susceptible to these often overlooked drivers. For instance, the flooded area could increase more than 2-fold (∼2.4) if Hurricane Irene were preceded by an extreme antecedent soil moisture condition (AMC). The results not only support the use of a multi-component framework for interactive flooding processes, but also underscore the necessity of broader definitions of compound flooding that encompass the simultaneous occurrence of intense precipitation, storm surge, and high AMC during TCs.
This study presents a coupled thermo-hydro-mechanical (THM) model for simulating the heating and hydration behavior of bentonite, a buffer material in deep geological repositories (DGRs). The model incorporates a new temperature-dependent soil water retention curve which captures the thermal-induced shift in water retention behavior. It also distinguishes between liquid and gas permeability, modeling intrinsic gas permeability as a function of accessible porosity to improve vapor transport and desaturation predictions. The model was validated against two large-scale tank tests, demonstrating good agreement with measured temperature, relative humidity, and water inflow data. It revealed a complex porosity evolution driven by thermal expansion, vapor movement, vapor condensation, and hydration-induced swelling during heating and hydration processes. The simulation results also suggest that the permeability of the hydration layer plays a critical role in controlling water intake. Clogging of this layer can significantly reduce the volume of water inflow during the hydration phase. Furthermore, while the model effectively captures key THM behavior, further development of the mechanical constitutive law is required to account for possible thermo-elasto-plastic volume changes and microstructural effects. Overall, the model provides a robust tool for evaluating the evolution of bentonite-based barrier material in DGRs.
We present results from the NASA GISS ModelE2.1‐G‐CC Earth System Model with coupled climate‐carbon cycle simulations that were submitted to the sixth phase of the Coupled Model Intercomparison Project (CMIP6) Coupled Climate‐Carbon Cycle MIP (C4MIP). Atmospheric CO2 concentration and carbon budgets for the land and ocean in the historical simulations were generally consistent with observations. Low simulated atmospheric CO2 concentrations during 1850‐1950 were due to excess uptake from prescribed land cover change, which erroneously replaced arid shrublands with higher biomass crops, and assumed high 2004 LAI values in vegetated lands throughout the historical simulation. At the end of the historical period, slightly higher simulated CO2 than observed resulted from the land being an insufficient net carbon sink, despite the net effect of CO2 fertilization and warming‐induced increases to leaf photosynthetic capacity. The global ocean carbon uptake agreed well with the observations with the largest discrepancies in the low latitudes. Future climate projection at 2091‐2100 agreed with CMIP5 models in the northward shift, of temperate deciduous forest climate and expansion across Eurasia along 60 °N latitude, and dramatic regional biome shifts from drying and warming in continental Europe. Carbon feedback parameters were largely similar to the CMIP5 model ensemble. For our model, the variation of land feedback parameters within the uncertainty arises from the fertilization feedback being less sensitive due to lack of increased vegetation growth, and the comparably more negative ocean carbon‐climate feedback is due to the large slowdown of the Atlantic overturning circulation.
This dataset presents a suite of downscaled hydro-climate projections over the conterminous United States (CONUS) based on multiple selected Global Climate Models (GCMs) from the Coupled Models Intercomparison Project phase 6 (CMIP6). The CMIP6 GCMs are downscaled dynamically using Regional Climate Model version 4 (RegCM4). Each ensemble member covers the 1980-2019 baseline and 2020-2059 near-future periods under the high-end (SSP585) emission scenario. This dataset is formulated to support the SECURE Water Act Section 9505 Assessment for the US Department of Energy (DOE) Water Power Technologies Office (WPTO). For further details on this dataset, please refer to Kao et al. (2022) and Rastogi et al. (2022).
This dataset presents a suite of hydropower projections for the conterminous United States (CONUS), derived from multiple downscaled and bias-corrected Global Climate Models (GCMs) from the Coupled Model Intercomparison Project Phase 6 (CMIP6). The CMIP6 GCMs are downscaled using either statistical (DBCCA) or dynamical (RegCM) approaches, based on two meteorological reference datasets (Daymet and Livneh). The resulting downscaled precipitation, temperature, and wind speed data are then used to drive two calibrated hydrologic models (VIC and PRMS), enabling simulations of projected future hydrologic responses across the CONUS. Simulated total runoff is subsequently employed to drive two hydropower models (WMP, now implemented as mosartwmpy-power, and WRES) to evaluate how climate change may affect future hydropower production for both federal and non-federal hydropower fleets. This dataset was developed to support the SECURE Water Act Section 9505 Assessment for the U.S. Department of Energy (DOE) Water Power Technologies Office (WPTO). For further details, see Broman et al. (2024), Thurber et al. (2024), Kao et al. (2022), and Zhou et al. (2023).
Near-surface mercury and ozone depletion events occur in the lowest part of the atmosphere during Arctic spring. Mercury depletion is the first step in a process that transforms long-lived elemental mercury to more reactive forms within the Arctic that are deposited to the cryosphere, ocean, and other surfaces, which can ultimately get integrated into the Arctic food web. Depletion of both mercury and ozone occur due to the presence of reactive halogen radicals that are released from snow, ice, and aerosols. In this work, we added a detailed description of the Arctic atmospheric mercury cycle to our recently published version of the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem 4.3.3) that includes Arctic bromine and chlorine chemistry and activation/recycling on snow and aerosols. The major advantage of our modelling approach is the online calculation of bromine concentrations and emission/recycling that is required to simulate the hourly and daily variability of Arctic mercury depletion. We used this model to study coupling between reactive cycling of mercury, ozone, and bromine during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) spring season in 2020 and evaluated results compared to land-based, ship-based, and remote sensing observations. The model predicts that elemental mercury oxidation is driven largely by bromine chemistry and that particulate mercury is the major form of oxidized mercury. The model predicts that the majority (74%) of oxidized mercury deposited to land-based snow is re-emitted to the atmosphere as gaseous elemental mercury, while a minor fraction (4%) of oxidized mercury that is deposited to sea ice is re-emitted during spring. Our work demonstrates that hourly differences in bromine/ozone chemistry in the atmosphere must be considered to capture the springtime Arctic mercury cycle, including its integration into the cryosphere and ocean.
The CCSD(T) coupled-cluster model with perturbative triples is considered a gold standard for computational modeling of the correlated behavior of electrons in molecular systems. A fundamental constraint is the relatively small global-memory capacity in GPUs compared to the main-memory capacity on host nodes, necessitating relatively smaller tile sizes for high-dimensional tensor contractions in NWChem's GPU-accelerated implementation of the CCSD(T) method. A coordinated redesign is described to address this limitation and associated data movement overheads, including a novel fused GPU kernel for a set of tensor contractions, along with inter-node communication optimization and data caching. The new implementation of GPU-accelerated CCSD(T) improves overall performance by 3.4x. Finally, we discuss the trade-offs in using this fused algorithm on current and future supercomputing platforms.
A major tropospheric loss of soluble species such as nitric acid results from scavenging by water droplets. Several theoretical formulations have been advanced which relate an effective time-independent loss rate for soluble species to statistical properties of precipitation such as the wet fraction and length of a precipitation cycle. In this paper, various 'effective' loss rates that have been proposed are compared with the results of detailed time-dependent model calculations carried out over a seasonal time scale. The model is a stochastic precipitation model coupled to a tropospheric photochemical model. The results of numerous time-dependent seasonal model runs are used to derive numerical values for the nitric acid residence time for several assumed sets of preciptation statistics. These values are then compared with the results obtained by utilizing theoretical 'effective' loss rates in time-independent models.
In this paper, we develop a general effective theory for two copies of the Sachdev-Ye-Kitaev (SYK) model with a time-dependent bilinear coupling. For a quantum quench problem with an initial state of the thermofield double state, we show how the evolution of the system is described by a complex reparametrization field with a classical Hamiltonian. We study correlation functions in this system and compare the large-q theory with the bulk low energy effective theory. In particular, we study the special case of a “rescued black hole”, which describes how a time-evolved thermofield double state can evolve to the ground state of a coupled SYK model by a carefully tuned time-dependent coupling. In the low energy region, there is a holographic dual interpretation, which is a geometry that crosses over from an eternal black hole to a global AdS2 vacuum. This family of geometries allow us to access the bulk region that would be the black hole interior without the rescue process. By comparing the large-q and low energy theory, we find that even in the low energy region the deviation from the low energy theory cannot be neglected if the rescue process starts late. This provides evidence that the low energy effective theory of the bulk fails near the inner horizon of the black hole. We note the possibility of a connection to a two-dimensional analog of the higher-dimensional black hole singularity.
The Sea-Viewing Wide Field-of-view Sensor (SeaWiFS) has observed 2.5 years of routine global chlorophyll observations from space. The mission was launched into a record El Nino event, which eventually gave way to one of the most intensive and longest-lasting La Nina events ever recorded. The SeaWiFS chlorophyll record captured the response of ocean phytoplankton to these significant events in the tropical Indo-Pacific basins, but also indicated significant interannual variability unrelated to the El Nino/La Nina events. This included large variability in the North Atlantic and Pacific basins, in the North Central and equatorial Atlantic, and milder patterns in the North Central Pacific. This SeaWiFS record was tracked with a coupled physical/biogeochemical/radiative model of the global oceans using near-real-time forcing data such as wind stresses, sea surface temperatures, and sea ice. This provided an opportunity to offer physically and biogeochemically meaningful explanations of the variability observed in the SeaWiFS data set, since the causal mechanisms and interrelationships of the model are completely understood. The coupled model was able to represent the seasonal distributions of chlorophyll during the SeaWiFS era, and was capable of differentiating among the widely different processes and dynamics occurring in the global oceans. The model was also reasonably successful in representing the interannual signal, especially when it was large, such as, the El Nino and La Nina events in the tropical Pacific and Indian Oceans. The model provided different phytoplankton group responses for the different events in these regions: diatoms were predominant in the tropical Pacific during the La Nina but other groups were predominant during El Nino. The opposite condition occurred in the tropical Indian Ocean. Both situations were due to the different responses of the basins to El Nino. The interannual variability in the North Atlantic, which was exhibited in SeaWiFS data as a decline in the spring/summer bloom in 1999 relative to 1998, resulted in the model from a more slowly shoaling mixed layer, allowing herbivore populations to keep pace with increasing phytoplankton populations. However, several aspects of the interannual cycle were not well-represented by the model. Explanations ranged from inherent model deficiencies, to monthly averaging of forcing fields, to biases in SeaWiFS atmospheric correction procedures.
This study examines the impact of changing the lateral diffusion coefficient ARedi on the transport of the Antarctic Circumpolar Current (ACC). The lateral diffusion coefficient ARedi is poorly constrained, with values ranging across an order of magnitude in climate models. The ACC is difficult to accurately simulate, and there is a large spread in eastward transport in the Southern Ocean (SO) in these models. This paper examines how much of that spread can be attributed to different eddy parameterization coefficients. A coarse-resolution, fully coupled model suite was run with A Redi = 400, 800, 1200, and 2400 m 2 s -1 . Additionally, two simulations were run with two-dimensional representations of the mixing coefficient based on satellite altimetry. Relative to the 400 m 2 s -1 case, the 2400 m 2 s -1 case exhibits 1) an 11% decrease in average wind stress from 50° to 65°S, 2) a 20% decrease in zonally averaged eastward transport in the SO, and 3) a 14% weaker transport through the Drake Passage. The decrease in transport is well explained by changes in the thermal current shear, largely due to increases in ocean density occurring on the northern side of the ACC. In intermediate waters these increases are associated with changes in the formation of intermediate waters in the North Pacific. We hypothesize that the deep increases are associated with changes in the wind stress curl allowing Antarctic Bottom Water to escape and flow northward.
Nitrous oxide (N 2 O) is a potent and persistent greenhouse gas, with rising atmospheric concentrations driven in part by inefficient use of synthetic nitrogen (N) fertilizers in agriculture. Predicting soil N 2 O emissions is challenging due to high spatial and temporal variability arising from complex soil biogeochemical processes. Process-based ecosystem models and standalone machine learning (ML) approaches without extensive site-specific calibration often miss high-emission episodes. Here, we show how an Ensemble Modeling System (EMS) based on outputs from an ensemble of ecosystem models coupled to an ensemble of ML models can improve predictions and understanding of N 2 O fluxes from US cropland. Trained and validated on ~12,000 N 2 O chamber measurements at 17 US Midwest sites (six crops, 35 management practices), the EMS accurately predicted daily fluxes of N 2 O at both training (R 2 = 0.84, RMSE = 16.4 g N ha −1 d −1 ) and held-out testing sites (R 2 = 0.84, RMSE = 6.2 g N ha −1 d −1 ). Analyses identified six dominant N 2 O drivers: soil organic carbon (SOC), NH 4 + , NO 3 - , water-filled pore space, temperature, and aboveground biomass production. Wet, warm soils produced large N 2 O peaks only with sufficient SOC and mineral N; in low-SOC soils, fluxes remained low. Incorporating these drivers into process-based models might significantly improve their predictive capacity. The EMS demonstrates a strong potential to predict N 2 O fluxes at unseen sites, enabling more reliable regional inventories, improved gap-filling where measurements are sparse, and enhanced understanding of mechanisms to advance targeted mitigation strategies in food, feed, and bioenergy crops.
During the last two plus decades, The Goddard Earth Observing System (GEOS) and Massachusetts Institute of Technology (MIT) modeling groups have developed, respectively, atmosphere-only and ocean-only global general circulation models. These two models (GEOS and MIT-GCM (General Circulation Model)) have demonstrated their data assimilation capabilities with the recent releases of the Modern Era Reanalysis for Research Applications, Version 2 (MERRA-2) atmospheric reanalysis and the Estimating the Circulation and Climate of the Ocean, Version 4 (ECCO-v4) ocean (and sea ice) state estimate. Independently, the two modeling groups have also produced global atmosphere-only and ocean-only simulations with km-scale grid spacing which proved invaluable for process studies and for the development of satellite and in-situ sampling strategies. Recently, a new effort has been made to couple these two models and to leverage their data-assimilation and high resolution capabilities (i.e., eddy-permitting ocean, cloud-permitting atmosphere). The focus in the model development is put on sub-seasonal to decadal time scales. In this talk, I discuss the new coupled model and present some first coupled simulation results. This will include a high-resolution coupled GEOS-MIT simulation, whereby we have coupled a cubed-sphere-720 (approximately 1/8 degrees) configuration of the GEOS atmosphere to a latitude-longitude-cap-1080 (approximately 1/12 degrees) configuration of the MIT ocean. We compare near-surface diagnostics of this fully coupled ocean-atmosphere set-up to equivalent atmosphere-only and ocean-only simulations. In the comparisons we focus in particular on the differences in air-sea interactions between sea surface temperature (SST) and wind for the coupled and uncoupled simulations.