SEARCH · Search NASA

Results for “global ocean”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Ice-shelf freshwater triggers for the Filchner–Ronne Ice Shelf melt tipping point in a global ocean–sea-ice model

Some ocean modeling studies have identified a potential tipping point from a low to a high basal melt regime beneath the Filchner–Ronne Ice Shelf (FRIS), Antarctica, with significant implications for subsequent Antarctic ice sheet mass loss. To date, investigation of the climate drivers and impacts of this possible event have been limited because ice-shelf cavities and ice-shelf melting are only now starting to be included in global climate models. Using a global ocean–sea-ice configuration of the Energy Exascale Earth System Model (E3SM) that represents both ocean circulations and melting within ice-shelf cavities, we explore freshwater triggers (iceberg melt and ice-shelf basal melt) of a transition to a high-melt regime at FRIS in a low-resolution (30 km in the Southern Ocean) global ocean–sea-ice model. We find that a realistic spatial distribution of iceberg melt fluxes is necessary to prevent the FRIS melt regime change from unrealistically occurring under historical-reanalysis-based atmospheric forcing. Further, improvement of the default parameterization for mesoscale eddy mixing significantly reduces a large regional fresh bias and weak Antarctic Slope Front structure, both of which precondition the model to melt regime change. Using two different stable model con figurations, we explore the sensitivity of FRIS melt regime change to regional ice-sheet freshwater fluxes. Through a series of sensitivity experiments prescribing incrementally increasing melt rates from the smaller, neighboring ice shelves in the eastern Weddell Sea, we demonstrate the potential for an ice-shelf melt “domino effect” should the upstream ice shelves experience increased melt rates. The experiments also reveal that modest ice-shelf melt biases in a model, especially at coarse ocean resolution where narrow continental shelf dynamics are not well resolved, can lead to an unrealistic melt regime change at downstream ice shelves. Thus, we find that remote connections between melt fluxes at different ice shelves are sensitive to baseline model conditions. Our results highlight both the potential and the peril of simulating prognostic Antarctic ice-shelf melt rates in a low-resolution global model.

58 GEOSCIENCES↗

DOE OSTI2024

Ice-shelf freshwater triggers for the Filchner–Ronne Ice Shelf melt tipping point in a global ocean–sea-ice model

Abstract. Some ocean modeling studies have identified a potential tipping point from a low to a high basal melt regime beneath the Filchner–Ronne Ice Shelf (FRIS), Antarctica, with significant implications for subsequent Antarctic ice sheet mass loss. To date, investigation of the climate drivers and impacts of this possible event have been limited because ice-shelf cavities and ice-shelf melting are only now starting to be included in global climate models. Using a global ocean–sea-ice configuration of the Energy Exascale Earth System Model (E3SM) that represents both ocean circulations and melting within ice-shelf cavities, we explore freshwater triggers (iceberg melt and ice-shelf basal melt) of a transition to a high-melt regime at FRIS in a low-resolution (30 km in the Southern Ocean) global ocean–sea-ice model. We find that a realistic spatial distribution of iceberg melt fluxes is necessary to prevent the FRIS melt regime change from unrealistically occurring under historical-reanalysis-based atmospheric forcing. Further, improvement of the default parameterization for mesoscale eddy mixing significantly reduces a large regional fresh bias and weak Antarctic Slope Front structure, both of which precondition the model to melt regime change. Using two different stable model configurations, we explore the sensitivity of FRIS melt regime change to regional ice-sheet freshwater fluxes. Through a series of sensitivity experiments prescribing incrementally increasing melt rates from the smaller, neighboring ice shelves in the eastern Weddell Sea, we demonstrate the potential for an ice-shelf melt “domino effect” should the upstream ice shelves experience increased melt rates. The experiments also reveal that modest ice-shelf melt biases in a model, especially at coarse ocean resolution where narrow continental shelf dynamics are not well resolved, can lead to an unrealistic melt regime change at downstream ice shelves. Thus, we find that remote connections between melt fluxes at different ice shelves are sensitive to baseline model conditions. Our results highlight both the potential and the peril of simulating prognostic Antarctic ice-shelf melt rates in a low-resolution global model.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Data-driven global ocean modeling for seasonal to decadal prediction

Accurate modeling of ocean dynamics is crucial for enhancing our understanding of complex ocean circulation processes, predicting climate variability, and tackling challenges posed by climate change. Although great efforts have been made to improve traditional numerical models, predicting global ocean variability over multiyear scales remains challenging. Here, we propose ORCA-DL (Oceanic Reliable foreCAst via Deep Learning), a data-driven three-dimensional ocean model for seasonal to decadal prediction of global ocean dynamics. ORCA-DL accurately simulates the three-dimensional structure of global ocean dynamics with high physical consistency and outperforms state-of-the-art numerical models in capturing extreme events, including El Niño–Southern Oscillation and upper ocean heat waves. Moreover, ORCA-DL stably emulates ocean dynamics at decadal timescales, demonstrating its potential even for skillful decadal predictions and climate projections. Our results demonstrate the high potential of data-driven models for providing efficient and accurate global ocean modeling and prediction.

Science & Technology - Other Topics↗

DOE OSTI2023

The large decline in carbonate $\delta^{238}$U from a PETM section at Tingri (South Tibet) was driven by local sea-level changes, not global oceanic anoxia

Uranium isotope compositions (δ 238 ) and Th/U in carbonates are being explored as paleoredox proxies to study global oceanic anoxia. However, the impact of changes to local depositional environments that might cause a decline in δ 238 U and an increase in Th/U has not been thoroughly investigated. Here we report a decline in δ 238 U (from +0.1 to -0.3‰) and an increase in Th/U from shallow-marine carbonates during the Paleocene-Eocene thermal maximum at Tingri (south Tibet). Changes of δ 238 U and Th/U are not fully coupled. The decline in δ 238 U coincided with a local sea-level fall, whereas the increase in Th/U occurred mainly during a sea-level rise. The decline in δ 238 U likely results from changes in redox conditions of pore waters and in primary carbonate mineralogy. The increase in Th/U is ascribed to reduced authigenic U(IV) accumulation. Here our results suggest that changes in local depositional environments can cause behavior of δ 238 U and Th/U resembling that induced by expanded global oceanic anoxia. This finding calls for caution to consider local factors before applying δ 238 U and Th/U as global paleoredox proxies.

58 GEOSCIENCES↗

DOE OSTI2024

A Scale‐Dependent Analysis of the Barotropic Vorticity Budget in a Global Ocean Simulation

Abstract The climatological mean barotropic vorticity budget is analyzed to investigate the relative importance of surface wind stress, topography, planetary vorticity advection, and nonlinear advection in dynamical balances in a global ocean simulation. In addition to a pronounced regional variability in vorticity balances, the relative magnitudes of vorticity budget terms strongly depend on the length‐scale of interest. To carry out a length‐scale dependent vorticity analysis in different ocean basins, vorticity budget terms are spatially coarse‐grained. At length‐scales greater than 1,000 km, the dynamics closely follow the Topographic‐Sverdrup balance in which bottom pressure torque, surface wind stress curl and planetary vorticity advection terms are in balance. In contrast, when including all length‐scales resolved by the model, bottom pressure torque and nonlinear advection terms dominate the vorticity budget (Topographic‐Nonlinear balance), which suggests a prominent role of oceanic eddies, which are of km in size, and the associated bottom pressure anomalies in local vorticity balances at length‐scales smaller than 1,000 km. Overall, there is a transition from the Topographic‐Nonlinear regime at scales smaller than 1,000 km to the Topographic‐Sverdrup regime at length‐scales greater than 1,000 km. These dynamical balances hold across all ocean basins; however, interpretations of the dominant vorticity balances depend on the level of spatial filtering or the effective model resolution. On the other hand, the contribution of bottom and lateral friction terms in the barotropic vorticity budget remains small and is significant only near sea‐land boundaries, where bottom stress and horizontal viscous friction generally peak.

Khatri, Hemant↗

DOE OSTI2023

Spatio–Temporal Coarse–Graining Decomposition of the Global Ocean Geostrophic Kinetic Energy

We expand on a recent determination of the first global energy spectrum of the ocean’s surface geostrophic circulation using a coarse-graining (CG) method. We compare spectra from CG to those from spherical harmonics by treating land in a manner consistent with the boundary conditions. While the two methods yield qualitatively consistent domain-averaged results, spherical harmonics spectra are too noisy at gyre-scales (> 1000 km). More importantly, spherical harmonics are inherently global and cannot provide local information connecting scales with currents geographically. CG shows that the extra-tropics mesoscales (100–500 km) have a root-mean-square (rms) velocity of ~15 cm/s, which increases to ~30–40 cm/s locally in the Gulf Stream and Kuroshio and to ~16–28 cm/s in the ACC. There is notable hemispheric asymmetry in mesoscale energy-per-area, which is higher in the north due to continental boundaries. We estimate that ≈25–50% of total geostrophic energy is at scales smaller than 100 km, and is un(der)-resolved by pre-SWOT satellite products. Spectra of the time-mean circulation show that most of its energy (up to 70%) resides in stationary eddies with characteristic scales smaller than (< 500 km). This highlights the preponderance of ‘standing’ small-scale structures in the global ocean due to the temporally coherent forcing by boundaries. By coarse-graining in space and time, we compute the first spatio-temporal global spectrum of geostrophic circulation from AVISO and NEMO. These spectra show that every length-scale evolves over a wide range of time-scales with a consistent peak at ≈200 km and ≈2–3 weeks.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Observational ozone datasets over the global oceans and polar regions (version 2024)

Studying tropospheric ozone over the remote areas of the planet, such as the open oceans and the polar regions, is crucial to understand the role of ozone as a global climate forcer and regulator of atmospheric oxidative capacity. A focus on the pristine oceanic and polar regions complements the available land-based datasets and provides insights into key photochemical and depositional loss processes that control the concentrations and spatiotemporal variability in ozone as well as the physicochemical mechanisms driving these patterns. However, an assessment of the role of ozone over the oceanic and polar regions has been hampered by a lack of comprehensive observational datasets. Here, we present the first comprehensive collection of ozone data over the oceans and the polar regions. The overall dataset consists of 77 ship cruises/buoy-based observations and 48 aircraft-based campaigns. The dataset, consisting of more than 630 000 independent ozone measurement data points covering the period from 1977 to 2022 and an altitude range from the surface to 5000 m (with a focus on the lowest 2000 m), allows systematic analyses of the spatiotemporal distribution and long-term trends over the 11 defined ocean/polar regions. The datasets from ships, buoys, and aircraft are complemented by ozonesonde data from 29 launch sites or field campaigns and by 21 non-polar and 17 polar ground-based station datasets. The datasets contain information on how long the observed air masses were isolated from land, as estimated by backward trajectories from the individual observation points. To extract observations representative of oceanic conditions, we recommend using a subset of the data with an isolation time of 72 h or longer, from the analysis with coincident radon observations. These filtered oceanic and polar data showed typically flat diurnal cycles at high latitudes, whereas daytime decreases in ozone (11 %–16 %) were observed at lower latitudes. The ship/buoy- and aircraft-based datasets presented here will supplement the land-based ones in the TOAR-II (Tropospheric Ozone Assessment Report Phase II) database to provide a fully global assessment of tropospheric ozone. The described dataset is available at https://doi.org/10.17596/0004044 (Kanaya et al., 2025).

Kanaya, Yugo [Japan Agency for Marine-Earth Scienc↗

DOE OSTI2023

Effective climate sensitivity distributions from a 1D model of global ocean and land temperature trends, 1970–2021

Abstract Current theoretically based Earth system models (ESMs) produce Effective Climate Sensitivities (EffCS) that range over a factor of three, with 80% of those models producing stronger global warming trends for 1970–2021 than do observations. To make a more observationally based estimate of EffCS, a 1D time-dependent forcing-feedback model of temperature departures from energy equilibrium is used to match measured ranges of global-average surface and sub-surface land and ocean temperature trends during 1970–2021. In response to two different radiative forcing scenarios, a full range of three model free parameters are evaluated to produce fits to a range of observed surface temperature trends (± 2σ) from four different land datasets and three ocean datasets, as well as deep-ocean temperature trends and borehole-based trend retrievals over land. Land-derived EffCS are larger than over the ocean, and EffCS is lower using the newer Shared Socioeconomic Pathways (SSP245, 1.86 °C global EffCS, ± 34% range 1.48–2.15 °C) than the older Representative Concentration Pathway forcing (RCP6, 2.49 °C global average EffCS, ± 34% range 2.04–2.87 °C). The strongest dependence of the EffCS results is on the assumed radiative forcing dataset, underscoring the role of radiative forcing uncertainty in determining the sensitivity of the climate system to increasing greenhouse gas concentrations from observations alone. The results are consistent with previous observation-based studies that concluded EffCS during the observational period is on the low end of the range produced by current ESMs.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Drifter-based global ocean current energy resource assessment

Not Available

Sadoughipour, Mahsan↗

DOE OSTI2024

Global Ocean Cooling of 2.3°C During the Last Glacial Maximum

Abstract Quantitative constraints on past mean ocean temperature (MOT) critically inform our historical understanding of Earth's energy balance. A recently developed MOT proxy based on paleoatmospheric Xe, Kr, and N 2 ratios in ice core air bubbles is a promising tool rooted in the temperature dependences of gas solubilities. However, these inert gases are systematically undersaturated in the modern ocean interior, and it remains unclear how air‐sea disequilibrium may have changed in the past. Here, we carry out 30 tracer‐enabled model simulations under varying circulation, sea ice cover, and wind stress regimes to evaluate air‐sea disequilibrium in the Last Glacial Maximum (LGM) ocean. We find that undersaturation of all three gases was likely reduced, primarily due to strengthened high‐latitude winds, biasing reconstructed MOT by −0.38 ± 0.37°C (1σ). Accounting for air‐sea disequilibrium, paleoatmospheric inert gases indicate that LGM MOT was 2.27 ± 0.46°C (1σ) colder than the pre‐industrial era.

Seltzer, A. M.↗

DOE OSTI2024

West African Monsoon System’s Responses to Global Ocean–Regional Atmosphere Coupling

This study explores the added value (AV) of a regional Earth system model (ESM) compared to an atmosphere-only regional climate model (RCM) in simulating West African monsoon (WAM) rainfall. The primary goals are to foster discussions on the suitability of coupled RCMs for WAM projections and deepen our understanding of ocean–atmosphere coupling’s influence on the WAM system. The study employs results from dynamical downscaling of the ERA-Interim reanalysis and Max Plank Institute ESM, low resolution (MPI-ESM-LR), by two RCMs, atmosphere only (REMO) and REMO coupled with Max Planck Institute Ocean Model (MPIOM) (ROM), at ~25-km horizontal resolution. Results show that in regions distant from coupling domain boundaries such as West Africa (WA), constraint conditions from ERA-Interim are more beneficial than coupling effects. REMO, reliant on oceanic sea surface temperatures (SSTs) from observations and influenced by ERA-Interim, is biased under coupling conditions, although coupling offers potential advantages in representing heat and mass fluxes. Contrastingly, as intended, coupling improves SSTs and monsoon fluxes’ relationships under ESM-forced conditions. In this latter case, the coupling features a dipole-like spatial structure of AV, improving precipitation over the Guinea Coast but degrading precipitation over half of the Sahel. Our extensive examination of physical processes and mechanisms underpinning the WAM system supports the plausibility of AV. Additionally, we found that the monsoonal dynamics over the ocean respond to convective activity, with the Sahara–Sahel surface temperature gradient serving as the maintenance mechanism. While further efforts are needed to enhance the coupled RCM, we advocate for its use in the context of WAM rainfall forecasts and projections.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Subglacial Discharge Effects on Antarctic Ice‐Shelf Basal Melt and the Southern Ocean in a Global, Coupled Ocean—Sea‐Ice Model

Subglacial freshwater from beneath Antarctica enters the ocean at depth, enhancing ice-shelf melting and affecting Southern Ocean properties. To study these effects in an Antarctic-wide context, we use a continental-scale subglacial hydrology model that calculates grounding line freshwater flux for a global, coupled ocean—sea-ice model. We find that subglacial discharge impacts melt rates primarily through continental shelf temperature modification, contrasting with findings from regional studies that do not permit large-scale adjustments. The consequence is that Antarctic melt rates scale with subglacial discharge more strongly than inferred from regional studies. We also find that the addition of buoyancy at depth facilitates heat upwelling to the surface, resulting in higher sea ice volume downstream of cold ice shelves and lower sea ice volume downstream of warm ice shelves. This highlights the drawbacks of simplifications in previous global studies that deposit Antarctic meltwater at the ocean surface and find uniform ocean surface cooling and sea-ice growth. While the patterns we find are robust, we conclude that the addition of subglacial discharge at present-day rates has a small effect on basal melt rates, hydrography, and sea ice. However, stronger discharge can have significant effects and can even accelerate a shift from low to high melting for ice shelves close to such a tipping point. Finally, the importance of feedbacks between enhanced cavity overturning and continental shelf conditions poses a complication for parameterizing subglacial discharge effects on melting for ice-sheet projections that do not include a coupled ocean component.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Global seamless tidal simulation using a 3D unstructured-grid model (SCHISM v5.10.0)

We present a new 3D unstructured-grid global ocean model to study both tidal and nontidal processes, with a focus on the total water elevation. Unlike existing global ocean models, the new model resolves estuaries and rivers down to ~8 m without the need for grid nesting. The model is validated with both satellite and in situ observations for elevation, temperature, and salinity. Tidal elevation solutions have a mean complex root-mean-square error (RMSE) of 4.2 cm for M2 and 5.4 cm for all five major constituents in the deep ocean. The RMSEs for the other four constituents, S2, N2, K1, and O1, are, respectively, 2.05, 0.93, 2.08, and 1.34 cm). The nontidal residual assessed by a tide gauge dataset (GESLA) has a mean RMSE of 7 cm. For the first time ever, we demonstrate the potential for seamless simulation on a single mesh from the global ocean into several estuaries along the US West Coast. The model is able to accurately capture the total elevation, even at some upstream stations. The model can therefore potentially serve as the backbone of a global tide surge and compound flooding forecasting framework.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

The GFDL‐CM4X Climate Model Hierarchy, Part I: Model Description and Thermal Properties

We present the GFDL‐CM4X (Geophysical Fluid Dynamics Laboratory Climate Model version 4X) coupled climate model hierarchy. The primary application for CM4X is to investigate ocean and sea ice physics as part of a realistic coupled Earth climate model. CM4X utilizes an updated MOM6 (Modular Ocean Model version 6) ocean physics package relative to CM4.0, and there are two members of the hierarchy: one that uses a horizontal grid spacing of 0.25° (referred to as CM4X‐p25) and the other that uses a grid 0.125° (CM4X‐p125). CM4X also refines its atmospheric grid from the nominally 100 km (cubed sphere C96) of CM4.0–50 km (C192). Finally, CM4X simplifies the land model to allow for a more focused study of the role of ocean changes to global mean climate. CM4X‐p125 reaches a global ocean area mean heat flux imbalance of -0.02 W m -2 within $\mathcal{O}$ (150) years in a pre‐industrial simulation, and retains that thermally equilibrated state over the subsequent centuries. This 1850 thermal equilibrium is characterized by roughly 400 ZJ less ocean heat than present‐day, which corresponds to estimates for anthropogenic ocean heat uptake between 1870 and present‐day. CM4X‐p25 approaches its thermal equilibrium only after more than 1000 years, at which time its ocean has roughly 1100 ZJ more heat than its early 21st century ocean initial state. Furthermore, the root‐mean‐square sea surface temperature bias for historical simulations is roughly 20% smaller in CM4X‐p125 relative to CM4X‐p25 (and CM4.0). We offer the mesoscale dominance hypothesis for why CM4X‐p125 shows such favorable thermal equilibration properties.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Simulations of ocean deoxygenation in the historical era: insights from forced and coupled models

Ocean deoxygenation due to anthropogenic warming represents a major threat to marine ecosystems and fisheries. Challenges remain in simulating the modern observed changes in the dissolved oxygen (O 2 ). Here, we present an analysis of upper ocean (0-700m) deoxygenation in recent decades from a suite of the Coupled Model Intercomparison Project phase 6 (CMIP6) ocean biogeochemical simulations. The physics and biogeochemical simulations include both ocean-only (the Ocean Model Intercomparison Project Phase 1 and 2, OMIP1 and OMIP2) and coupled Earth system (CMIP6 Historical) configurations. We examine simulated changes in the O 2 inventory and ocean heat content (OHC) over the past 5 decades across models. The models simulate spatially divergent evolution of O 2 trends over the past 5 decades. The trend (multi-model mean and spread) for upper ocean global O 2 inventory for each of the MIP simulations over the past 5 decades is 0.03 ± 0.39×10 14 [mol/decade] for OMIP1, -0.37 ± 0.15×10 14 [mol/decade] for OMIP2, and -1.06 ± 0.68×10 14 [mol/decade] for CMIP6 Historical, respectively. The trend in the upper ocean global O2 inventory for the latest observations based on the World Ocean Database 2018 is -0.98×10 14 [mol/decade], in line with the CMIP6 Historical multi-model mean, though this recent observations-based trend estimate is weaker than previously reported trends. A comparison across ocean-only simulations from OMIP1 and OMIP2 suggests that differences in atmospheric forcing such as surface wind explain the simulated divergence across configurations in O 2 inventory changes. Additionally, a comparison of coupled model simulations from the CMIP6 Historical configuration indicates that differences in background mean states due to differences in spin-up duration and equilibrium states result in substantial differences in the climate change response of O 2 . Finally, we discuss gaps and uncertainties in both ocean biogeochemical simulations and observations and explore possible future coordinated ocean biogeochemistry simulations to fill in gaps and unravel the mechanisms controlling the O 2 changes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

A nonhydrostatic formulation for MPAS-Ocean

The Model for Prediction Across Scales-Ocean (MPAS-Ocean) is an open-source, global ocean model and is one component of a family of climate models within the MPAS framework, including atmosphere, sea-ice, and land-ice models. Here, in this work, a new formulation for the ocean model is presented that solves the nonhydrostatic, incompressible Boussinesq equations on an unstructured, staggered, z-level grid. The introduction of this nonhydrostatic capability is necessary for the resolution of internal wave dynamics and large eddy simulations. Compared to the standard, hydrostatic formulation, a nonhydrostatic pressure solver and a vertical momentum equation are added, where the PETSc (Portable Extensible Toolkit for Scientific Computation) library is used for the inversion of a large sparse system for the nonhydrostatic pressure. Numerical results on a stratified seiche, internal solitary wave, overflow and lock-exchange test cases are presented, and the parallel efficiency of the code is evaluated using up to 1024 processors.

3D Poisson equation↗

DOE OSTI2025

A nonhydrostatic formulation for MPAS-Ocean

3D Poisson equation↗

DOE OSTI2025

Development and evaluation of a new 4DEnVar-based weakly coupled ocean data assimilation system in E3SMv2

The development, implementation, and evaluation of a new weakly coupled ocean data assimilation (WCODA) system for the fully coupled Energy Exascale Earth System Model version 2 (E3SMv2) utilizing the four-dimensional ensemble variational (4DEnVar) method are presented in this study. The 4DEnVar method, based on the dimension-reduced projection four-dimensional variational (DRP-4DVar) approach, replaces the adjoint model with the ensemble technique, thereby reducing computational demands. Monthly mean ocean temperature and salinity data from the EN4.2.1 reanalysis are integrated into the ocean component of E3SMv2 from 1950 to 2021 with the goal of providing realistic initial conditions for decadal predictions and predictability studies. The performance of the WCODA system is assessed using various metrics, including the reduction rate of the cost function, root mean square error (RMSE) differences, correlation differences, and model biases. Results indicate that the WCODA system effectively assimilates the reanalysis data into the climate model, consistently achieving negative reduction rates of the cost function and notable improvements in RMSE and correlation across various ocean layers and regions. Significant enhancements are observed in the upper ocean layers across the majority of global ocean regions, particularly in the north Atlantic, north Pacific, and Indian Ocean. Model biases in sea surface temperature and salinity are also substantially reduced. For sea surface temperature, cold biases in the north Pacific and north Atlantic are diminished by about 1–2 °C, and warm biases in the Southern Ocean are corrected by approximately 1.5–2.5 °C. In terms of salinity, improvements are observed with bias reductions of about 0.5–1 psu in the north Atlantic and north Pacific and up to 1.5 psu in parts of the Southern Ocean. The ultimate goal of the WCODA system is to advance the predictive capabilities of E3SM for subseasonal to decadal climate predictions, thereby supporting research on strategic energy-sector policies and planning.

54 ENVIRONMENTAL SCIENCES↗