Search NASA⌕ Search

SEARCH · Search NASA

Results for “ocean models”

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 127 records · Page 7

North Pacific Decadal Variability in the GEOS-5 Atmosphere-Ocean Model

This study examines the mechanisms of the Pacific decadal oscillation (PDO) in the GEOS-5 general circulation model. The model simulates a realistic PDO pattern that is resolved as the first empirical orthogonal function (EOF) of winter sea surface temperature (SST). The simulated PDO is primarily forced by Aleutian low through Ekman transport and surface fluxes, and shows a red spectrum without any preferred periodicity. This differs from the observations, which indicate a greater role of El Nino-Southern Oscillation (ENSO) forcing, and likely reflects the too short time scale of the simulated ENSO. The geostrophic transport in response to the Aleutian low is limited to the Kuroshio-Oyashio Extension, and is unlikely the main controlling factor in this model, although it reinforces the Ekman-induced SST anomalies. The delay between the Aleutian low and the PDO is relatively short (1 year) suggesting that the fast Ekman response (rather than Rossby wave propagation) sets the SST pattern immediately following an Aleutian low fluctuation. The atmospheric feedback (response to the SST) is only about 25 of the forcing and never evolves into an Aleutian low completely, instead projecting onto the North Pacific Oscillation (NPO), a meridional dipole in sea level pressure (SLP). The lack of preferred periodicity and weak atmospheric response bothindicate a coupled oscillation is an unlikely mechanism for the PDO in this model. In agreement with recent studies, the NPO is correlated with the North Pacific Gyre Oscillation (NPGO), which is another leading EOF of the North Pacific SST. A possible connection between the PDO and the NPGO is discussed.

Pacific Decadal Variability↗

Numerical integrators for Lagrangian oceanography

Abstract. A common task in Lagrangian oceanography is to calculate a large number of drifter trajectories from a velocity field precalculated with an ocean model. Mathematically, this is simply numerical integration of an ordinary differential equation (ODE), for which a wide range of different methods exist. However, the discrete nature of the modelled ocean currents requires interpolation of the velocity field in both space and time, and the choice of interpolation scheme has implications for the accuracy and efficiency of the different numerical ODE methods. We investigate trajectory calculation in modelled ocean currents with 800 m, 4 km, and 20 km horizontal resolution, in combination with linear, cubic and quintic spline interpolation. We use fixed-step Runge–Kutta integrators of orders 1–4, as well as three variable-step Runge–Kutta methods (Bogacki–Shampine 3(2), Dormand–Prince 5(4) and 8(7)). Additionally, we design and test modified special-purpose variants of the three variable-step integrators, which are better able to handle discontinuous derivatives in an interpolated velocity field. Our results show that the optimal choice of ODE integrator depends on the resolution of the ocean model, the degree of interpolation, and the desired accuracy. For cubic interpolation, the commonly used Dormand–Prince 5(4) is rarely the most efficient choice. We find that in many cases, our special-purpose integrators can improve accuracy by many orders of magnitude over their standard counterparts, with no increase in computational effort. Equivalently, the special-purpose integrators can provide the same accuracy as standard methods at a reduced computational cost. The best results are seen for coarser resolutions (4 and 20 km), thus the special-purpose integrators are particularly advantageous for research using regional to global ocean models to compute large numbers of trajectories. Our results are also applicable to trajectory computations on data from atmospheric models.

58 GEOSCIENCES↗

Skill Assessment of a Spectral Ocean-Atmosphere Radiative Model

Ocean phytoplankton, detrital material, and water absorb and scatter light spectrally. The Ocean- Atmosphere Spectral Irradiance Model (OASIM) is intended to provide surface irradiance over the oceans with sufficient spectral resolution to support ocean ecology, biogeochemistry, and heat exchange investigations, and of sufficient duration to support inter-annual and decadal investigations. OASIM total surface irradiance (integrated 200 nm to 4 microns) was compared to in situ data and three publicly available global data products at monthly 1-degree resolution. OASIM spectrally-integrated surface irradiance had root mean square (RMS) difference= 20.1 W/sq m (about 11%), bias=1.6 W/sq m (about 0.8%), regression slope= 1.01 and correlation coefficient= 0.89, when compared to 2322 in situ observations. OASIM had the lowest bias of any of the global data products evaluated (ISCCP-FD, NCEP, and ISLSCP 11), and the best slope (nearest to unity). It had the second best RMS, and the third best correlation coefficient. OASIM total surface irradiance compared well with ISCCP-FD (RMS= 20.7 W/sq m; bias=-11.4 W/sq m, r=0.98) and ISLSCP II (RMS =25.2 W/sq m; bias= -13.8 W/sq m; r=0.97), but less well with NCEP (RMS =43.0 W/sq m ;bias=-22.6 W/sq m; x=0.91). Comparisons of OASIM photosynthetically available radiation (PAR) with PAR derived from SeaWiFS showed low bias (-1.8 mol photons /sq m/d, or about 5%), RMS (4.25 mol photons /sq m/d ' or about 12%), near unity slope (1.03) and high correlation coefficient (0.97). Coupled with previous estimates of clear sky spectral irradiance in OASIM (6.6% RMS at 1 nm resolution), these results suggest that OASIM provides reasonable estimates of surface broadband and spectral irradiance in the oceans, and can support studies on ocean ecosystems, carbon cycling, and heat exchange.

Gregg, Watson, W.↗

A Modeling Study of Oceanic Response to Daily and Monthly Surface Forcing

The goal of this study is to investigate the effect of high-frequency surface forcing (wind stresses and heat fluxes) on upper-ocean response. We use the reduced-gravity quasi-isopycnal ocean model by Schopf and Loughe (1995) for this study. Two experiments are performed: one with daily and the other with monthly surface forcing. The two experiments are referred to as DD and MM, respectively. The daily surface wind stress is produced from the SSM/I wind data (Atlas et al. 1991) using the drag coefficient of Large and Pond (1982). The surface latent and sensible heat fluxes are estimated using the atmospheric mixed layer model by Seager et al. (1995) with the time-varying air temperature and specific humidity from the NCEP-NCAR reanalysis (Kalnay et al. 1996). The radiation is based on climatological shortwave radiation from the Earth Radiation Budget Experiment (ERBE) [Harrison et al. 1993] and the daily GEWEX SRB data. The ocean model domain is restricted to the Pacific Ocean with realistic land boundaries. At the southern boundary the model temperature and salinity are relaxed to the Levitus (1994) climatology. The time-mean SST distribution from MM is close to the observed SST climatology while the mean SST field from DD is about 1.5 C cooler. To identify the responsible processes, we examined the mean heat budgets and the heat balance during the first year (when the difference developed) in the two experiments. The analysis reveals that this is contributed by two factors. One is the difference in latent heat flux. The other is the difference in mixing processes. To further evaluate the responsible processes, we repeated the DD experiment by reducing the based vertical diffusion from 1e-4 to 0.5e-5. The resultant SST field becomes quite closer to the observed SST field. SST variability from the two experiments is generally similar, but the equatorial SST differences between the two experiments show interannual variations. We are investigating the possible mechanisms responsible for the different responses.

Sui, Chung-Hsiung↗

Streamlining Ocean Dynamics Modeling with Fourier Neural Operators: A Multiobjective Hyperparameter and Architecture Optimization Approach

Training an effective deep learning model to learn ocean processes involves careful choices of various hyperparameters. We leverage DeepHyper’s advanced search algorithms for multiobjective optimization, streamlining the development of neural networks tailored for ocean modeling. The focus is on optimizing Fourier neural operators (FNOs), a data-driven model capable of simulating complex ocean behaviors. Selecting the correct model and tuning the hyperparameters are challenging tasks, requiring much effort to ensure model accuracy. DeepHyper allows efficient exploration of hyperparameters associated with data preprocessing, FNO architecture-related hyperparameters, and various model training strategies. We aim to obtain an optimal set of hyperparameters leading to the most performant model. Moreover, on top of the commonly used mean squared error for model training, we propose adopting the negative anomaly correlation coefficient as the additional loss term to improve model performance and investigate the potential trade-off between the two terms. The numerical experiments show that the optimal set of hyperparameters enhanced model performance in single timestepping forecasting and greatly exceeded the baseline configuration in the autoregressive rollout for long-horizon forecasting up to 30 days. Utilizing DeepHyper, we demonstrate an approach to enhance the use of FNO in ocean dynamics forecasting, offering a scalable solution with improved precision.

97 MATHEMATICS AND COMPUTING↗

Effects of Atmosphere and Ocean Horizontal Model Resolution on Tropical Cyclone and Upper-Ocean Response Forecasts in Four Major Hurricanes

A coupled atmosphere–ocean model is necessary for tropical cyclone (TC) prediction to accurately characterize ocean feedback on atmospheric processes within the TC environment. Here, the ECMWF coupled global model is run at horizontal resolutions from 9 to 1.4 km in the atmosphere, as well as 25 and 8 km in the ocean, to identify how resolution impacts forecast accuracy of four observed major TCs in the Atlantic: Irma, Florence, Teddy, and Ida. Most of the resolutions used here are unprecedented for global models. GOES-16 and synthetic aperture radar (SAR) satellite images and best track data are used for atmospheric validation. Salinity and temperature observations from Air-Launched Autonomous Micro-Observer (ALAMO) floats are used to validate modeled upper-ocean response, including mixed layer deepening, sea surface cooling, and near-inertial waves in the wakes of TCs. Increasing atmospheric resolution leads to more realistic TC structure and stronger winds, significantly improving TC intensity forecasts and modestly improving track errors. Ocean resolution impacts the upper-ocean response but does not influence atmospheric forecasts for the fast-moving TCs considered here. Stronger mixing, sea surface cooling, and near-inertial oscillations are found for both higher atmosphere and ocean resolutions, provided the initial upper-ocean state is the same for the two ocean resolutions. Whether this agrees better with the ALAMO observations also depends on the realism of the initial upper-ocean state in the model, emphasizing the importance of ocean initialization for the accurate upper-ocean response. Overall, the model at all resolutions correctly predicts stronger mixing, surface cooling, and near-inertial oscillation amplitudes to the right of a TC center, as observed by ALAMO floats.

Atmosphere-ocean interaction↗

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↗