skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Learning ocean circulation models with reservoir computing
Two elementary models of ocean circulation, the well-known double-gyre stream function model and a single-layer quasi-geostrophic (QG) basin model, are used to generate flow data that sample a range of possible dynamical behavior for particular flow parameters. A reservoir computing (RC) machine learning algorithm then learns these models from the stream function time series. In the case of the QG model, a system of partial differential equations with three physically relevant dimensionless parameters is solved, including Munk- and Stommel-type solutions. The effectiveness of a RC approach to learning these ocean circulation models is evident from its ability to capture the characteristics of these ocean circulation models with limited data including predictive forecasts. Further assessment of the accuracy and usefulness of the RC approach is conducted by evaluating the role of both physical and numerical parameters and by comparison with particle trajectories and with well-established quantitative assessments, including finite-time Lyapunov exponents and proper orthogonal decomposition. The results show the capability of the methods outlined in this article to be applied to key research problems on ocean transport, such as predictive modeling or control.  more » « less
Award ID(s):
2121919 2121923
PAR ID:
10391030
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Physics of Fluids
Volume:
34
Issue:
11
ISSN:
1070-6631
Page Range / eLocation ID:
116604
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract A simplified quasigeostrophic (QG) analytical model together with an idealized numerical model are used to study the effect of uneven ice–ocean stress on the temporal evolution of the geostrophic current under sea ice. The tendency of the geostrophic velocity in the QG model is given as a function of the lateral gradient of vertical velocity and is further related to the ice–ocean stress with consideration of a surface boundary layer. Combining the analytical and numerical solutions, we demonstrate that the uneven stress between the ice and an initially surface-intensified, laterally sheared geostrophic current can drive an overturning circulation to trigger the displacement of isopycnals and modify the vertical structure of the geostrophic velocity. When the near-surface isopycnals become tilted in the opposite direction to the deeper ones, a subsurface velocity core is generated (via geostrophic setup). This mechanism should help understand the formation of subsurface currents in the edge of Chukchi and Beaufort Seas seen in observations. Furthermore, our solutions reveal a reversed flow extending from the bottom to the middepth, suggesting that the ice-induced overturning circulation potentially influences the currents in the deep layers of the Arctic Ocean, such as the Atlantic Water boundary current. 
    more » « less
  2. Abstract Parameterization of mesoscale eddies in coarse resolution ocean models is necessary to include the effect of eddies on the large‐scale oceanic circulation. We propose to use a multiple‐scale Quasi‐Geostrophic (MSQG) model to capture the eddy dynamics that develop in response to a prescribed large‐scale flow. The MSQG model consists in extending the traditional quasi geostrophic (QG) dynamics to include the effects of a variable Coriolis parameter and variable background stratification. Solutions to this MSQG equation are computed numerically and compared to a full primitive equation model. The large‐scale flow field permits baroclinically unstable QG waves to grow. These instabilities saturate due to non‐linearities and a filtering method is applied to remove large‐scale structures that develop due to the upscale cascade. The resulting eddy field represents a dynamically consistent response to the prescribed background flow, and can be used to rectify the large‐scale dynamics. Comparisons between Gent‐McWilliams eddy parameterization and the present solutions show large regions of agreement, while also indicating areas where the eddies feed back onto the large scale in a manner that the Gent‐McWilliams parameterization cannot capture. Also of interest is the time variability of the eddy feedback which can be used to build stochastic eddy parameterizations. 
    more » « less
  3. null (Ed.)
    In this paper, we propose Test-Time Training, a general approach for improving the performance of predictive models when training and test data come from different distributions. We turn a single unlabeled test sample into a self-supervised learning problem, on which we update the model parameters before making a prediction. This also extends naturally to data in an online stream. Our simple approach leads to improvements on diverse image classification benchmarks aimed at evaluating robustness to distribution shifts. 
    more » « less
  4. Existing physics-based modeling approaches do not have a good compromise between performance and computational efficiency in predicting the seismic response of reinforced concrete (RC) frames, where high-fidelity models (e.g., fiber-based modeling method) have reasonable predictive performance but are computationally demanding, while more simplified models (e.g., shear building model) are the opposite. This paper proposes a novel artificial intelligence (AI)-enhanced computational method for seismic response prediction of RC frames which can remedy these problems. The proposed AI-enhanced method incorporates an AI technique with a shear building model, where the AI technique can directly utilize the real-world experimental data of RC columns to determine the lateral stiffness of each column in the target RC frame while the structural stiffness matrix is efficiently formulated via the shear building model. Therefore, this scheme can enhance prediction accuracy due to the use of real-world data while maintaining high computational efficiency due to the incorporation of the shear building model. Two data-driven seismic response solvers are developed to implement the proposed approach based on a database including 272 RC column specimens. Numerical results demonstrate that compared to the experimental data, the proposed method outperforms the fiber-based modeling approach in both prediction capability and computational efficiency and is a promising tool for accurate and efficient seismic response prediction of structural systems. 
    more » « less
  5. Abstract Surface runoff and infiltrated water en route to the stream interact with dynamic landscape properties, ranging from vegetation and microbial activities to soil and geological attributes. Stream solute concentrations are highly variable and interconnected due to these interactions, flow paths, and residence times, and often exhibit hysteresis with flow. Significant unknowns remain about how point measurements of stream solute chemistry reflect interdependent hydrobiogeochemical and physical processes, and how signatures are encapsulated as nonlinear dynamical relationships between variables. We take a Machine Learning (ML) approach to understand and capture these dynamical relationships and improve predictions of solutes at short and long time scales. We introduce a physical process‐based “flow‐gate” into an Long Short‐Term Memory (LSTM) model, which enables the model to learn hysteresis behaviors if they exist. Further, we use information‐theoretic metrics to detect how solutes are interdependent and iteratively select source solutes that best predict a given target solute concentration. The “flow‐gate LSTM” model improves model predictions (1%–32% decreases in RMSE) relative to the standard LSTM model for all nine solutes included in the study. The predictive improvements from the flow‐gate LSTM model highlight the importance of lagged concentration and discharge relationships for certain solutes. It also indicates a potential limitation in the traditional LSTM model approach since flow rates are always provided as input sources, but this information is not fully utilized. This work provides a starting point for a predictive understanding of geochemical interdependencies using machine‐learning approaches and highlights potential improvements in model architecture. 
    more » « less