SEARCH · Search NASA
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.
This paper introduces Omega, the Ocean Model for E3SM Global Applications. Omega is a new ocean model designed to run efficiently on high performance computing (HPC) platforms, including exascale heterogeneous architectures with accelerators, such as Graphics Processing Units (GPUs). Omega is written in C and uses the Kokkos performance portability library. These were chosen because they are well-supported and will help future-proof Omega for upcoming HPC architectures. Omega will eventually replace the Model for Prediction Across Scales-Ocean (MPAS-Ocean) in the US Department of Energy's (DOE's) Energy Exascale Earth System Model (E3SM). Omega runs on unstructured horizontal meshes with variable-resolution capability and implements the same horizontal discretization as MPAS-Ocean. This work documents the design and performance of Omega Version 0.1.0 (Omega-V0), which solves the shallow water equations with passive tracers and is the first step towards the full primitive equation ocean model. On Central Processing Units (CPUs), Omega-V0 is 1.4 times faster than MPAS-Ocean with the same configuration. Omega-V0 is more efficient on GPUs than CPUs on a per-watt basis – by a factor of 5.3 on Frontier and 3.6 on Aurora, two of the world's fastest exascale computers.
Monitoring Sahelian rainfall variability is increasingly critical as climate extremes intensify across the region. Here, we develop the Sahelian Monsoon Ocean-Pressure Index (SMOPI), a novel global synthetic indicator constructed from five dynamically coherent sea-level pressure regions statistically linked to June-September Sahel monsoon rainfall. SMOPI captures intra-seasonal and interannual variability, and crucially, reflects the influence of both regional processes and large-scale teleconnections on monsoon dynamics. It aligns with the dominant rainfall variability mode in reanalyses and 29 CMIP6 models. Strong/positive SMOPI phases coincide with wet years and are associated with enhanced convergence, favorable jet configurations, and robust Pacific, Atlantic, and Indian Ocean teleconnections. Conversely, weak/negative SMOPI phases correspond to drought conditions and divergent moisture fluxes. SMOPI exposes model failures in reproducing historical droughts and offers new physical insights into rainfall-driving mechanisms. It stands out as a scalable, potentially transferable diagnostic tool for monitoring/forecasting and evaluating Sahelian monsoon rainfall under global warming.
Abstract Symmetric instability (SI) extracts kinetic energy from fronts in the surface mixed layer (SML), potentially affecting the SML structure and dynamics. Here, a global submesoscale-permitting ocean model named MITgcm LLC4320 simulation is used to examine the Stone linear prediction of the maximum SI scale to estimate grid spacings needed to begin resolving SI. Furthermore, potential effects of SI on the usable wind work are estimated roughly: this estimate of SI “activity” is useful for assessing if these modes should be resolved or parameterized. The maximum SI scale varies by latitude with median values from 568 to 23 m. Strong seasonality is observed in the SI scale and activity. The median scale in winter is 188 m globally, 2.5 times of that of summer (75 m). SI is more active in winter: 15% of the time compared with 6% in summer. The strongest SI activity is found in the western Pacific, western Atlantic, and Southern Oceans. The required grid spacings for a global model to begin resolving SI eddies in the SML are 24 m (50% of regions resolved) and 7.9 m (90%) in winter, decreasing to 9.4 m (50%) and 3.6 m (90%) in summer. It is also estimated that SI may reduce usable wind work by an upper bound of 0.83 mW m −2 globally, or 5% of the global magnitude. The sensitivity of these estimates to empirical thresholds is provided in the text.
For coastal adaptation purposes, it is important to estimate the climate related changes in extreme sea levels due to storm surges and ocean waves, in addition to mean sea level rise. In this study, we provide the first consistent and continuous estimation of projected changes in global storm surges and ocean waves from the past to the warmer future, based on an extremely high resolution global climate model. The spatial pattern in the trend of annual maximum sea surface heights and wave heights is predominantly driven by changes in tropical cyclone (TC) frequency. In the western North Pacific, future TC frequencies are projected to decrease, and the annual maximum sea surface heights and wave heights show decreasing trends (-20 cm/century and -200 cm/century). Although highly intense tropical cyclones are enhanced in the warmer climate, highly extreme storm surges and wave heights do not necessarily increase due to the large natural variability.
Global and regional water cycle includes precipitation, water vapor divergence, and change of column water vapor in the atmosphere, and land surface evapotranspiration, terrestrial water storage change, and river discharge, which is linked to ocean salinity near the river mouth. The water cycle is a crucial component of the Earth system, and numerous studies have addressed its individual components (e.g., precipitation). Here we assess, for the first time, if remote sensing and reanalysis datasets can accurately and self-consistently portray the Amazon water cycle. This is further assisted with satellite ocean salinity measurements near the mouth of the Amazon River. The widely used practice of taking the mean of an ensemble of datasets to represent water cycle components (e.g., precipitation) can produce large biases in water cycle closure. Additionally, closure is achieved with only a small subset of data combinations (e.g., ERA5 precipitation and evapotranspiration plus GRACE satellite terrestrial water storage), which rules out the lower precipitation and higher evapotranspiration estimates, providing valuable constraints on assessments of precipitation, evapotranspiration, and their ratio. The common approach of using the Óbidos stream gauge (located hundreds of kilometers from the river mouth) multiplied by a constant (1.25) to represent the entire Amazon discharge is found to misrepresent the seasonal cycle, and this can affect the apparent influence of Amazon discharge on tropical Atlantic salinity.
Abstract The global tropical cyclone (TC) number has historically been relatively constant from year‐to‐year, however, the reason remains unknown. Furthermore, climate projections are inconclusive regarding future global TC frequency changes. Here, we investigated years in which observed global TC activity deviated from the mean and potential links to ocean drivers from 1980 to 2021. We found that the annual global number of named storm days and accumulated cyclone energy (ACE) were significantly linked with El Niño–Southern Oscillation (ENSO) and the Atlantic Meridional Mode (AMM). La Niña and positive AMM are associated with the bottom percentiles of both TC metrics, and vice versa for El Niño and negative AMM. The ENSO Longitude Index explains variability in annual global named storm days and ACE as well as the Niño 3.4 index. This research reveals that reliable future projections of ENSO are necessary, but not sufficient, to understand future changes in global TC frequency.
The US Department of Energy (DOE) conducts climate simulations on some of the world’s largest supercomputers. These exascale machines use heterogeneous architectures with both CPUs and GPUs, and scientific codes must adapt to make full use of this computing power. Los Alamos National Lab is developing Omega: The Ocean Model for E3SM Global Applications, which is specifically designed for modern exascale computers. It uses external libraries that have been optimized for a variety of architectures to run on different supercomputers. Omega is an unstructured-mesh ocean model based on TRiSK numerical methods. It will be the new ocean component of the DOE’s Energy Exascale Earth System Model (E3SM). The algorithms in Omega follow those of the current ocean component, MPAS-Ocean, but it will be written in C++ rather than Fortran to take advantage of the Kokkos performance portability library. Omega spatial operators are written as Kokkos kernels to run efficiently on both CPUs and GPUs. Work on Omega began in 2023 with a new C++ framework for unstructured mesh partitioning, halo exchanges, parallel IO, and Kokkos interfaces. The current version, Omega-0, is being developed to solve the shallow water equations and at present includes all of the tendency terms but not time stepping. Here we share the results of Omega-0 verification and performance testing. Verification includes unit tests implemented with CTest as well as convergence tests in Polaris, an in-house python package with a large suite of test problems. Performance tests compare simulations conducted on CPUs versus GPUs and across different architectures: tests are run on Frontier, which has AMD “Optimized 3rd Gen EPYC” CPUs and AMD MI250X GPUs, as well as Perlmutter, which is composed of AMD EPYC 7763 CPUs and NVIDIA A100 GPUs.
Omega version 0, the first step toward a new ocean model, has been developed to take full advantage of next-generation supercomputers, including those powered by GPUs. Omega is rebuilt from the ground up, redesigning computational aspects of ocean models to run faster and more efficiently across diverse high-performance computing systems.
Abstract Dissolved organic carbon (DOC) in coastal waters is integral to biogeochemical cycling, but global and regional drivers of DOC are still uncertain. In this study we explored spatial and temporal differences in DOC concentrations and stocks across the global coastal ocean, and how these relate to temperature and salinity. We estimated a global median coastal DOC stock of 3.15 Pg C (interquartile range (IQR) = 0.85 Pg C), with median DOC concentrations being 2.2 times higher than in open ocean surface waters. Globally and seasonally, salinity was the main driver of DOC with concentrations correlated negatively with salinity, without a clear relationship to temperature. DOC concentrations and stocks varied with region and season and this pattern is likely driven by riverine inputs of DOC and nutrients that stimulate coastal phytoplankton production. Temporally, high DOC concentrations occurred mainly in months with high freshwater input, with some exceptions such as in Eastern Boundary Current margins where peaks are related to primary production stimulated by nutrients upwelled from the adjacent ocean. No spatial trend between DOC and temperature was apparent, but many regions (19 out of 25) had aligned peaks of seasonal temperature and DOC, related to increased phytoplankton production and vertical stratification at high temperatures. Links of coastal DOC with salinity and temperature highlight the potential for anthropogenic impacts to alter coastal DOC concentration and composition, and thereby ecosystem status.
An important technique used by climate modelers to isolate the impacts of increasing greenhouse gasses on Earth System processes is to simulate the impact of an abrupt increase in carbon dioxide. The spatial pattern of change provides a “fingerprint” that is generally much larger than natural variability. Insofar as the response to radiative forcing is linear (the impact of quadrupling CO 2 is twice the impact of doubling CO 2 ) this fingerprint can then be used to estimate the impact of historical greenhouse gas forcing. However, the degree to which biogeochemical cycles respond linearly to radiative forcing has rarely been tested. In this paper, we evaluate which ocean biogeochemical fields are likely to respond linearly to changing radiative forcing, which ones do not, and where linearity breaks down. We also demonstrate that the representation of lateral mixing by mesoscale eddies, which varies significantly across climate models, plays an important role in modulating the breakdown of linearity. Globally integrated surface rates of biogeochemical cycling (primary productivity, particulate export) respond in a relatively linear fashion and are only moderately sensitive to mixing. By contrast, the habitability of the interior ocean (as determined by hypoxia and calcite supersaturation) behaves non-linearly and is very sensitive to mixing. This is because the deep ocean, as well as certain regions in the surface ocean, are very sensitive to the magnitude of deep wintertime convection. The cessation of convection under global warming is strongly modulated by the representation of eddy mixing.
The Joint Task Force, Science Monitoring And Reliable Telecommunications (JTF SMART) Subsea Cables, is working to integrate environmental sensors for ocean bottom temperature, pressure, and seismic acceleration into submarine telecommunications cables. The purpose of SMART Cables is to support climate and ocean observation, sea level monitoring, observations of Earth structure, and tsunami and earthquake early warning and disaster risk reduction, including hazard quantification. Recent advances include regional SMART pilot systems that are the first steps to trans-ocean and global implementation. Examples of pilots include: InSEA wet demonstration project off Sicily at the European Multidisciplinary Seafloor and water column Observatory Western Ionian Facility; New Caledonia and Vanuatu; French Polynesia Natitua South system connecting Tahiti to Tubaui to the south; Indonesia starting with short pilot systems working toward systems for the Sumatra-Java megathrust zone; and the CAM-2 ring system connecting Lisbon, Azores, and Madeira. This paper describes observing system simulations for these and other regions. Funding reflects a blend of government, development bank, philanthropic foundation, and commercial contributions. In addition to notable scientific and societal benefits, the telecommunications enterprise’s mission of global connectivity will benefit directly, as environmental awareness improves both the integrity of individual cable systems as well as the resilience of the overall global communications network. SMART cables support the outcomes of a predicted, safe, and transparent ocean as envisioned by the UN Decade of Ocean Science for Sustainable Development and the Blue Economy. As a continuation of the OceanObs’19 conference and community white paper (Howe et al., 2019, 10.3389/fmars.2019.00424), an overview of the SMART programme and a description of the status of ongoing projects are given.
Not Available
In response to greenhouse gas forcing, climate models predict that poleward ocean heat transport (OHT) weakens in the Southern Ocean but increases in the Arctic. The role of ocean circulation changes in this OHT response has been evaluated by decomposing OHT anomalies into a dynamic component (holding ocean temperature fixed while circulation evolves) and a thermodynamic component (holding ocean circulation fixed while temperature evolves). However, ocean temperature changes are themselves shaped by circulation changes through redistribution of the existing heat reservoir and subsequent effects on air–sea heat fluxes. The thermodynamic component can therefore be influenced by circulation changes, making the standard thermodynamic–dynamic decomposition incomplete for isolating the role of circulation changes in OHT anomalies. To address this issue, we use a passive–active decomposition to assess the relative contributions of ocean circulation and passive ocean temperature changes to polar OHT anomalies in a fully coupled climate model. Passive temperature changes are defined as those thermally forced by the atmosphere in the absence of circulation changes. In this passive–active decomposition, an advective term involving both circulation and passive temperature changes remains ambiguous—classifying it as active implies circulation changes dominate Southern Ocean OHT anomalies, whereas classifying it as passive implies temperature changes dominate. However, both interpretations imply that ocean circulation changes have a much weaker effect on polar OHT anomalies than inferred from the standard decomposition. In conclusion, these results help reconcile conclusions from studies using the standard decomposition with those using passive tracer methods to assess the role of circulation changes in polar OHT anomalies.
In a warming climate, the Global Meridional Overturning Circulation (GMOC) is expected to change significantly with a risk of disrupting the global redistribution of ocean properties that sustains marine ecosystems, carbon cycle, and others. Here we make a novel attempt to utilize a diagnostic ocean & sea-ice model to estimate the GMOC and its interdecadal changes since the mid-1950s that are consistent with historical hydrographic observations. We find that significant changes in the GMOC have already occurred, most notably in the upper and lower overturning cells in the Southern Ocean. The former has expanded poleward and into denser water and strengthened by 3–4 Sv since the mid-1970s, while the latter has contracted and weakened by a similar rate during the same period. These changes are driven by the increasing Southern Hemisphere (SH) Ferrel cell and associated increases in the westerlies and the surface buoyancy loss over its sinking branch, and the increasing Antarctic meltwater discharge, in response to ozone depletion in the SH stratosphere and increasing atmospheric CO 2 . A large-scale readjustment of the GMOC seems to be underway in the South Atlantic and Indo-Pacific Oceans since the mid-2000s in response to the Southern Ocean changes.
Anthropogenic climate change will drive extensive mass loss across both the Antarctic (AIS) and Greenland Ice Sheets (GrIS), with the potential for global climate system feedbacks, especially in polar regions. Historically, the high-latitude North Atlantic and Southern Ocean have been critical regions for anthropogenic heat and carbon uptake, but our understanding of how this uptake will be altered by future freshwater discharge is incomplete. We assess each ice sheet's impact on global ocean anthropogenic heat and carbon storage for a high-emission scenario over the 21 st -century using a coupled Earth system model. We explore the impact of contemporaneous mass loss from both ice sheets on anthropogenic heat and carbon storage and quantify their linear and nonlinear contributions. Notably, added freshwater reduces ocean heat and carbon storage by 2,100, and the sum of individual freshwater effects differ from those induced by simultaneous freshwater discharge from both ice sheets. Combined AIS and GrIS freshwater engenders distinct anthropogenic storage anomalies—particularly in the high-latitude Southern Ocean and North Atlantic. From 2080 to 2100, GrIS freshwater exerts primary control on the temporal evolution of global ocean heat storage, while global ocean carbon storage is modulated by the linear AIS and GrIS freshwater impacts. Nonlinear impacts of simultaneous ice sheet discharge have a non-negligible contribution to the evolution of global ocean heat storage. Further, anthropogenic heat changes are realized more quickly in response to ice sheet discharge than anthropogenic carbon. Our results highlight the need to incorporate both ice sheets actively in climate models to accurately project future global climate.
Anomalously warm ocean waters have been implicated as the cause of accelerated ablation along the margin of the Greenland ice sheet, as well as the cause of increasing basal melt rate, mass loss, and grounding line retreat of many Antarctic ice shelves. The overarching science objective of this project is to investigate how the interplay of regional processes and decadal changes in local and remote forcing impact the delivery and end member composition of waters over the continental shelves of Greenland and Antarctica. To address this goal, an existing global eddy-active ocean general circulation model (OGCM), coupled to a thermodynamic/dynamic ice model and forced with atmospheric reanalysis fluxes from 1948-2009, was analyzed to understand these processes around the continental margins of Greenland and Antarctica. This simulation did not contain a representation of ice sheet melt or resolve mesoscale eddies over these continental margins, so two new coupled OGCM/sea-ice models were constructed and analyzed to address the roles of small scale features and processes (< 5km) and freshwater flux surpluses.
While wind is the primary driver of the oceanic general circulation, it kills the ocean’s most energetic motions.