Next Article in Journal
Prediction Analysis of Sea Level Change in the China Adjacent Seas Based on Singular Spectrum Analysis and Long Short-Term Memory Network
Previous Article in Journal
Dynamic Analysis of Crane Vessel and Floating Wind Turbine during Temporary Berthing for Offshore On-Site Maintenance Operations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Method for Predicting High-Resolution 3D Variations in Temperature and Salinity Fields Using Multi-Source Ocean Data

1
College of Intelligent Systems Science and Engineering, Harbin Engineering University, Harbin 150001, China
2
Qingdao Hatran Ocean Intelligence Technology Co., Ltd., Qingdao 266400, China
3
Key Laboratory of Physical Oceanography, MOE, Institute for Advanced Ocean Study, Frontiers Science Center for Deep Ocean Multispheres and Earth System (DOMES), Ocean University of China, Qingdao 266100, China
4
The College of Ocean and Atmosphere, Ocean University of China, Qingdao 266100, China
5
Innovation and Development Center, Harbin Engineering University, Qingdao 266400, China
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2024, 12(8), 1396; https://doi.org/10.3390/jmse12081396
Submission received: 9 July 2024 / Revised: 3 August 2024 / Accepted: 12 August 2024 / Published: 14 August 2024
(This article belongs to the Section Physical Oceanography)

Abstract

:
High-resolution three-dimensional (3D) variations in ocean temperature and salinity fields are of great significance for ocean environment monitoring. Currently, AI-based 3D temperature and salinity field predictions rely on expensive 3D data, and as the prediction period increases, the stacking of high-resolution 3D data greatly increases the difficulty of model training. This paper transforms the prediction of 3D temperature and salinity into the prediction of sea surface elements and the inversion of subsurface temperature and salinity using sea surface elements, by leveraging the relationship between sea surface factors and subsurface temperature and salinity. This method comprehensively utilizes multi-source ocean data to avoid the issue of data volume caused by stacking high-resolution historical data. Specifically, the model first utilizes 1/4° low-resolution satellite remote sensing data to construct prediction models for sea surface temperature (SST) and sea level anomaly (SLA), and then uses 1/12° high-resolution temperature and salinity data as labels to build an inversion model of subsurface temperature and salinity based on SST and SLA. The prediction model and inversion model are integrated to obtain the final high-resolution 3D temperature and salinity prediction model. Experimental results show that the 20-day prediction results in the two sea areas of the coastal waters of China and the Northwest Pacific show good performance, accurately predicting ocean temperature and salinity in the vast majority of layers, and demonstrate higher resource utilization efficiency.

1. Introduction

High-resolution three-dimensional (3D) temperature and salinity fields provide refined and comprehensive data support and analysis for marine environmental monitoring [1,2,3,4,5]. Accurately predicting changes in the 3D ocean temperature and salinity fields is of great significance. Currently, there are two main methods for 3D temperature and salinity prediction. The first method is based on physical dynamics using numerical models [6,7,8], which simulate ocean movement processes through a series of complex dynamic equations. However, this approach requires discretization of the complex dynamic equations, and the iterative process consumes a significant amount of computational resources.
The second method is based on data-driven deep learning prediction, which predicts by learning the patterns of data changes. Compared with traditional numerical models, deep learning methods have significant advantages in terms of computational cost and flexibility and have been widely applied in the field of ocean temperature and salinity in recent years. For example, in temperature prediction, Zhang et al. [9] used the FC-LSTM model to estimate the time series of single-point SST, achieving higher accuracy than traditional SVM algorithms. However, single-point prediction models lack spatial correlation, and subsequent research [10,11,12,13,14,15] has improved algorithm models to further optimize the performance of sea temperature prediction algorithms.
In addition to improving model accuracy, combining physical methods has shown significant effects. For instance, Hou et al. [16] incorporated shortwave and longwave radiation into the prediction model, leveraging the physical relationship between SST and radiation to enhance SST prediction accuracy. Hu et al. [17] used deep learning algorithms to extract multi-variable physical features of the sea surface, improving the prediction accuracy of SST and other sea surface elements. Wu et al. [18] adopted a physics-guided deep learning algorithm to further enhance sea temperature prediction accuracy and improve the model’s physical consistency.
In terms of salinity, Zhang et al. [19] used DNN neural networks to retrieve sea surface salinity, demonstrating the model’s reliable performance in salinity learning. Liu et al. [20] proposed a neural network called DSE-NN, which effectively captures the spatiotemporal distribution characteristics and physical laws of salinity, achieving higher salinity reconstruction accuracy. Additionally, numerous studies [21,22,23,24,25] have shown that data-driven deep learning models can effectively learn the complex nonlinear relationships between input variables and salinity, resulting in more accurate salinity predictions.
Current research on ocean temperature and salinity prediction primarily focuses on the two-dimensional surface plane. However, for marine environmental monitoring, effective prediction in three-dimensional space is of greater significance. In recent years, Liu et al. [26] utilized TD-LSTM to predict the average seawater temperature at different depths for each observation point one month ahead. Zhang et al. [27] employed M-ConvLSTM to predict the 3D ocean temperature field. Due to the influence of thermoclines, which reduced prediction accuracy, Zuo et al. [28] proposed the SST-4D-CNN model to predict ocean temperature. Sun et al. [29] used the 3D-UNet model to predict SSbT above 400 m in the Pacific and its adjacent oceans. These methods require expensive 3D ocean data and have limitations in terms of temporal or spatial resolution of the predicted data. However, oceanic small-scale data play a critical role in detailed observations [30,31,32]. The current challenge in research is how to leverage existing ocean data to predict refined temperature and salinity information accurately.
Since many measurable sea surface elements can reflect internal oceanic changes, key internal ocean parameters can be inferred from them [33]. For example, Wu et al. [34] found that internal ocean dynamics, influenced by thermal expansion, can affect the sea surface, and SST can indirectly reflect changes in subsurface temperatures. Syst, E. [35] discovered that thermohaline expansion significantly contributes to sea level rise, which is closely related to ocean temperature. Based on this relationship, subsequent researchers [36,37,38,39,40] utilized sea surface elements to construct deep learning networks to invert subsurface temperature and salinity fields. The experimental results show that sea surface elements can effectively reconstruct 3D temperature and salinity fields.
The method proposed in this paper relies on the nonlinear relationship between SST, SLA, and subsurface 3D temperature and salinity fields, transforming 3D temperature and salinity prediction into horizontal sea surface element prediction and vertical sea surface to subsurface temperature and salinity inversion. The experimental steps are as follows: first, a 1/4° low-resolution satellite remote sensing data-based M-ViT model is used for SST and SLA prediction. The prediction results are then temporally and spatially matched with reanalysis data, and the SimVP-gsta model is used to invert the subsurface 3D temperature and salinity fields, thereby achieving a 1/12° high-resolution 3D temperature and salinity prediction model. Experimental results show that the proposed method can effectively predict changes in the subsurface 3D temperature and salinity fields within 20 days.
The rest of this paper is organized as follows: Section 2 will introduce the data used in the experiments, the research area, and the data partitioning methods. Section 3 will provide a detailed description of the model components selected for this study. Section 4 will cover the experimental parameters and the evaluation metrics used. Section 5 will validate the effectiveness of the three-dimensional forecasting model constructed in this study. Section 6 will discuss the experimental results, and finally, Section 7 will present the conclusions of the study.

2. Dataset

The selected research areas for this study are the Coastal Waters of China (CWOC, 0° N–42° N, 105° E–130° E) and the Northwest Pacific Ocean (NWP, 10° S–42° N, 99° E–180° E). The sea surface forecast data used in this study come from NOAA’s OISST data (https://www.ncei.noaa.gov/data, accessed on 10 October 2023) and daily mean SLA data provided by CMEMS (https://data.marine.copernicus.eu/products, accessed on 10 October 2023), both with a resolution of 1/4°.
For the three-dimensional temperature and salinity inversion, the input data consist of OISST data interpolated using bilinear interpolation to match the spatial resolution of low-resolution satellite data and reanalysis data. The target data are high-precision reanalysis data provided by CMEMS (https://data.marine.copernicus.eu/products, accessed on 15 October 2023), with a resolution of 1/12° and comprising 50 layers. For the experiment, layers 1–48 are used as the target dataset for inversion, with each layer representing depth as detailed in Table 1.
For this forecasting and inversion experiment, data from 2006 to 2017 were selected as the training set, and data from 2018 to 2020 were used as the validation set. Data from December 2021 to November 2022, approximately one year, were used as the test set, and test results were statistically analyzed on a quarterly basis. The details of the training data partitioning are shown in the figure: Figure 1a shows the partitioning of the sea surface element prediction data, and Figure 1b shows the partitioning of the sea surface element inversion subsurface temperature and salinity data.

3. Method

Due to the characteristics of the experimental data, including large information volume, multi-source data, and a high number of grid points, the complex information carried within the data is not conducive to model learning. Simplifying the data structure for the model to learn and deriving general patterns is key to building the model in this study.

3.1. Model Construction

Currently, the encoder–decoder structure of deep learning is widely used in related fields such as computer vision, for example, the UNet model [41]. The encoding process reduces the complexity of the data features, thus decreasing computational complexity and improving the model’s learning efficiency. The decoding process converts low-resolution feature maps into high-resolution outputs to capture image details and structures, and integrates them with contextual information from the decoding stage.
The SimVP-gsta model [42] (https://github.com/chengtan9907/OpenSTL, accessed on 21 September 2023) builds upon the original encoder–decoder structure by adding the gsta module (spatio-temporal attention mechanism) to the decoding structure, which accelerates the model’s convergence speed and improves operational efficiency. This model has demonstrated excellent performance in various datasets, including video prediction. The structure of this model is shown in Figure 2, where Figure 2a illustrates the overall framework of the model, and Figure 2b depicts the gsta module in the decoding structure.
The improved model in this paper, M-ViT, has a structure similar to the SimVP-gsta model, featuring an encoder–translator–decoder framework. The translator is a lightweight Mobile–ViT module [43] (https://github.com/jaiwei98/mobile-vit-pytorch, accessed on 13 September 2023). This model performed excellently in the spatiotemporal prediction in this experiment. The specific model is shown in Figure 3, where Figure 3a depicts the overall framework of the model, and Figure 3b shows the decoder structure, which is the Mobile–ViT module.
In the experiments for constructing the intelligent forecasting and inversion models, we compared the training results of SimVP-gsta, M-ViT, and other typical models that have been used in ocean research [36,38,40,44]. The best-performing model from this comparison was selected as the foundational model for the 3D temperature and salinity forecasting model in this study.

3.2. Construction of the 3D Prediction Model

The three-dimensional temperature and salinity intelligent forecasting model is an integrated system comprising three submodels: the sea surface intelligent forecasting model, the bilinear interpolation model, and the three-dimensional temperature and salinity inversion model based on combined SST and SLA. The diagram of the three-dimensional forecasting model is shown in Figure 4. The integration method for the three-dimensional forecasting model can be implemented as follows. First, obtain the SST and SLA for the next x time steps:
S S T t + 1 , S S T t + 2 , , S S T t + x = P r e d _ M o d e l 1 ( S S T t n , S S T t n + 1 , , S S T t )
S L A t + 1 , S L A t + 2 , , S L A t + x = P r e d _ M o d e l 2 ( S L A t n , S L A t n + 1 , , S L A t )
where P r e d _ M o d e l 1 is the SST forecasting model and P r e d _ M o d e l 2 is the SLA forecasting model. Then, perform bilinear interpolation on the predicted SST and SLA for the next x time steps and match them with the reanalysis data in both the spatial and temporal dimensions:
B _ S S T t + i = B i _ M o d e l ( S S T t + i ) , i x
B _ S L A t + i = B i _ M o d e l ( S L A t + i ) , i x
where B _ S S T t + i and B _ S L A t + i represent the SST and SLA states after bilinear interpolation, respectively. The bilinear interpolation model is illustrated in Figure 5.
Assuming the data at the point ( x , y ) and its surrounding four points ( x 0 , y 0 ) , ( x 0 , y 1 ) , ( x 1 , y 0 ) , and ( x 1 , y 1 ) are f ( x 0 , y 0 ) , f ( x 0 , y 1 ) , f ( x 1 , y 0 ) , and f ( x 1 , y 1 ) , respectively, f ( x 0 , y ) = f ( x 0 , y 0 ) + y [ f ( x 0 , y 1 ) f ( x 0 , y 0 ) ] , f ( x 1 , y ) = f ( x 1 , y 0 ) + y [ f ( x 1 , y 1 ) f ( x 1 , y 0 ) ] , the bilinear interpolation is given by f ( x , y ) = f ( x 0 , y ) + x [ f ( x 1 , y ) f ( x 0 , y ) ] . The interpolated SST and SLA values for the corresponding time are then jointly input into the inversion model to obtain the future subsurface temperature and salinity information:
S T t + i , S S t + i = I n v _ M o d e l ( B _ S S T t + i , B _ S L A t + i ) , i x
where S T t + i and S S t + i represent the subsurface temperature and salinity, respectively, and I n v _ M o d e l is the inversion model.
Figure 4. Three-dimensional prediction model.
Figure 4. Three-dimensional prediction model.
Jmse 12 01396 g004
Figure 5. Bilinear interpolation model.
Figure 5. Bilinear interpolation model.
Jmse 12 01396 g005

4. Parameter Settings and Evaluation Metrics

4.1. Parameter Settings

Before training, the input and target data are first normalized using Z-Score normalization. The forecast data are organized in the shape (B, T, H, W), where B is the batch size, T is the time step, H and W are the height and width of the data, respectively. The inversion data are organized in the shape (B, L, E, H, W), where L represents the ocean depth layers, and E is the number of input features. We provide statistics for the optimal final prediction model M-ViT and the optimal inversion model SimVP-gsta. The results are as follows. The data and model parameters are shown in Table 2. Here, “dim” represents the hidden layer dimension of the SimVP-gsta model, “NT” is the number of gsta modules, and “Ns” is the number of encoder–decoder modules. “C” denotes the hidden layer dimension of the M-ViT model, “L” is the number of Transformer layers used in the Mobile–ViT model, and “P” is the size of each data block during segmentation.
In this prediction and inversion experiment, the MSE function was selected as the loss function. Due to the different distribution patterns of temperature and salinity, a dual loss function was used to constrain the model in the inversion experiment, capturing the distribution patterns of the temperature and salinity data. The specific scheme is shown in Figure 6. This experiment was conducted on a single NVIDIA A100 GPU (Sourced from NVIDIA Corporation in Santa Clara, CA, USA) for training. The initial learning rate for both the prediction and inversion model training was set to 0.001, and OneCycleLR was employed for learning rate updates, continuously adjusting the network’s learning rate to find the optimal solution. The number of iterations for training the network was set to 100, and the Adam optimizer was chosen.

4.2. Evaluation Metrics

In this experiment, Mean Absolute Error (MAE) (Equation (6)) and Root Mean Squared Error (RMSE) (Equation (7)) were used as evaluation metrics for the prediction and inversion experiment results. The formulas are as follows:
MAE = 1 N i = 1 N | y i y ^ i |
RMSE = 1 N i = 1 N ( y i y ^ i ) 2
In the above formula, N represents the total number of grid points in the sample, y i denotes the observed true values, and h a t y i represents the predicted values.

5. Experimental Results

5.1. SST Prediction

The experiments used satellite remote sensing SST data to construct a temperature forecasting model for the offshore areas of China and the Northwest Pacific region. Historical SST data for 30 days were input into the model to predict temperature changes for the next 20 days. The experiments were conducted under the same training conditions, comparing UNet, SimVP-gsta, M-ViT, Linear Regression Model (L-R), ViT Model [45], and ConvFormer Model [46]. The training results were tested using a test set, and the best-performing model among the compared models was selected as the foundational model for the three-dimensional temperature and salinity forecasting model. Before computing the metrics, the land areas were removed, and the period from December 2021 to November 2022 was used as the test set. The results of the tests for the two regions are presented below.

5.1.1. Coastal Waters of China

To better describe the model’s prediction ability for each time period, the 20-day error metrics were averaged every 5 days and tabulated as shown in Table 3.
Analysis of the error metrics in Table 3: From the overall results of the 1–20 day forecasts, the M-ViT model demonstrates the best performance in SST prediction experiments compared to other deep learning models, with the SimVP-gsta model following. Compared to the SimVP-gsta model, the M-ViT model reduces the average MAE metric by 6.7% and the RMSE metric by 5.2% for the 20-day forecast. The linear regression model exhibits smaller prediction errors in the 1–5 day time range; however, its forecast error growth rate is the largest as time increases. To more intuitively express the test results, the error metrics over the 20 days are visualized for comparison (as shown in Figure 7). The figure clearly shows that the M-ViT model has the best predictive performance.

5.1.2. Northwest Pacific

Analysis of the error metrics in Table 4: The M-ViT model demonstrates better results compared to other prediction models, exhibiting lower overall MAE and RMSE. Specifically, compared to the SimVP-gsta model, the M-ViT model reduces the average MAE by 4.6% and the RMSE by 3.2% over the 20-day forecast period. The errors over these 20 days are visualized for comparison (as shown in Figure 8), where it is evident that the M-ViT model consistently maintains the best predictive performance throughout the entire period.
From the experimental results of the above two sea regions, it can be concluded that the M-ViT model exhibits the best predictive performance in overall temperature prediction. Therefore, this paper uses the M-ViT model as the base model for temperature prediction in the 3D temperature and salinity prediction experiment.

5.2. SLA Prediction

Thirty days of historical SLA data were input into the models to predict the changes in SLA for the next 20 days. The testing procedures were the same as in the SST prediction experiments. The following are the results of the tests conducted separately for the two sea areas. The models used for training in this section were the same as those used in the SST forecasting model.

5.2.1. Coastal Waters of China

To better describe the predictive capabilities of each model in each time period, we averaged the 20-day error metrics every 5 days and present them in a table as shown in Table 5.
Analysis of the error metrics in Table 5: The M-ViT model outperformed other deep learning models in predicting SLA, demonstrating the best overall performance. Compared to the SimVP-gsta model, the M-ViT model provided slightly better results, reducing the average MAE by 1.8% and the RMSE by 1.1% over the 20-day forecast period. The errors over these 20 days are visualized for comparison (as shown in Figure 9), where it can be seen that the M-ViT model delivers superior prediction results.

5.2.2. Northwest Pacific

Analysis of the error metrics in Table 6: The M-ViT and SimVP-gsta models provide better average prediction results over the 20-day period compared to other prediction models. Specifically, the SimVP-gsta model performs better in terms of average prediction accuracy for the 1–5 day range, while the M-ViT model has an advantage in the 6–20 day range. The error metrics over the 20 days are visualized for comparison (as shown in Figure 10). Figure 10 shows that the linear regression model has lower errors in short-term predictions but performs poorly in medium- to long-term forecasts.
Based on the results of the two sea area experiments mentioned above, the M-ViT model demonstrates the best predictive performance in overall SLA predictions. Therefore, this paper will use the M-ViT model as the baseline model for SLA predictions in the 3D temperature and salinity prediction experiments.

5.3. Ocean Subsurface Temperature and Salinity Vertical Inversion

The experiment first performed bilinear interpolation on the historical SST and SLA data at a 1/4° resolution from satellite remote sensing. These interpolated data were then matched spatially with 1/12° resolution reanalysis temperature and salinity data. These combined data sources were used to construct a subsurface temperature and salinity inversion model for the Chinese coastal and Northwest Pacific regions. The interpolated historical SST data and SLA data were combined and input into the model to retrieve temperature and salinity data for 48 layers below the sea surface (approximately 5000 m deep). Under the same training conditions, models including UNet, SimVP-gsta, M-ViT, ConvFormer, Atten-UNet, and ViT were trained simultaneously. The training results were evaluated using a test set, with the best-performing model among the compared models chosen as the base model for the three-dimensional temperature and salinity forecasting model. Before calculating the metrics, the land areas were excluded, and the period from December 2021 to November 2022 was used as the test set. The following are the test results for the two regions.

5.3.1. Coastal Waters of China

Due to the different distribution characteristics of temperature and salinity at various subsurface layers, the error metrics of the three models’ predictions were averaged separately for depths shallower than 600 m and deeper than 600 m. The results were then tabulated as shown in Table 7.
Analysis of Error Metrics in Table 7: The SimVP-gsta model outperformed other deep learning models overall in the temperature and salinity vertical inversion experiments. For temperature inversion, the results of the SimVP-gsta and M-ViT models are comparable, with SimVP-gsta being slightly better. However, in the salinity inversion, the SimVP-gsta model significantly outperformed the other models. The M-ViT model exhibited smaller errors at depths shallower than 600 m, but at depths greater than 600 m, it struggled to capture salinity distribution characteristics, leading to a noticeable increase in errors.
To better illustrate the inversion results, the error metrics for the 48-layer temperature and salinity inversions are visualized for comparison (as shown in Figure 11). Figure 11 shows that under the same experimental conditions, the SimVP-gsta model achieved the best overall results in temperature and salinity inversion. While many models tend to learn the temperature distribution pattern but struggle to simultaneously learn both temperature and salinity distribution patterns, the SimVP-gsta model excels at learning both. This balanced learning capability makes the SimVP-gsta model superior in terms of overall performance.

5.3.2. Northwest Pacific

Analysis of error metrics in Table 8: In terms of overall performance, the SimVP-gsta model excels in both temperature and salinity inversion. The M-ViT model also demonstrated strong inversion performance for temperature, but it exhibited larger errors in salinity inversion, particularly at depths greater than 600 m, where it struggled to maintain a balanced learning of salinity patterns.
The error metrics for the 48-layer temperature and salinity inversions are visualized for comparison (as shown in Figure 12). Figure 12 illustrates that, under the same experimental conditions, the SimVP-gsta model outperformed other deep learning models in overall inversion results. Although the UNet model showed relatively small errors in salinity inversion, it still experienced error fluctuations. Most models exhibited an imbalance in learning salinity, indicating that they find it challenging to simultaneously capture the distribution patterns of both temperature and salinity. In contrast, the SimVP-gsta model delivered better inversion results, demonstrating superior model performance.
Based on the results of the two sea area experiments mentioned above, the SimVP-gsta model demonstrated the best performance in overall temperature and salinity vertical inversion. Therefore, this paper will use the SimVP-gsta model as the baseline model for temperature and salinity inversion in the 3D temperature and salinity prediction experiments.

5.4. Three-Dimensional Temperature and Salinity Prediction

In this section, we integrate the SST prediction model, the SLA prediction model, and the ocean subsurface vertical inversion model to develop a 3D temperature and salinity prediction model for the ocean subsurface. To optimize the prediction performance, the best models from the above experiments were chosen as the foundational models for the 3D prediction model. The basic prediction steps are as follows: first, input the 1/4° SST data and SLA data into the M-ViT model to obtain 20-day prediction data. Then, apply bilinear interpolation to the prediction data and input them into the inversion model, the SimVP-gsta model, to obtain the final high-resolution 3D temperature and salinity prediction results.
To effectively test the prediction results while maintaining seasonal continuity, the period from December 2021 to November 2022 was selected as the test dataset. The data were divided into four quarters for testing and statistics: winter (December–February), spring (March–May), summer (June–August), and autumn (September–November). Before calculating the metrics, the land areas were excluded, and the data from the entire year of 2022 were used as the test set. The following are the test results for the two specific sea areas.

5.4.1. Coastal Waters of China

The average of the 20-day prediction results for each prediction within the entire quarter was calculated, and then the 20-day prediction results were averaged every 5 days. The statistical results are shown in Table 9.
Further analysis of Table 9 reveals that, during the 1–20 day prediction period, the temperature and salinity prediction errors gradually increased over time, which is consistent with the objective laws. By comparing the prediction errors from days 1 to 5 with those from days 16 to 20, the feasibility of this 3D prediction model can be analyzed.
For the 1–5 day prediction, the average MAE for temperature at depths less than 600 m ranges between 0.4950 °C and 0.5697 °C, and the average RMSE ranges between 0.6989 °C and 0.7392 °C. At depths greater than 600 m, the average MAE ranges between 0.1043 °C and 0.1107 °C, and the average RMSE ranges between 0.1513 °C and 0.1562 °C. For the 16–20 day prediction, the average MAE for temperature at depths less than 600 m ranges between 0.5531 °C and 0.6136 °C, and the average RMSE ranges between 0.7775 °C and 0.8225 °C. At depths greater than 600 m, the average MAE ranges between 0.1051 °C and 0.1225 °C, and the average RMSE ranges between 0.1526 °C and 0.1711 °C. Compared to the 1–5 day prediction, the average MAE and RMSE for temperatures at depths less than 600 m increase by about 0.10. At depths greater than 600 m, the average MAE and RMSE increase by about 0.06.
The average MAE for salinity predictions from 1 to 5 days ranges between 0.1897 psu and 0.2485 psu for depths less than 600 m, and the average RMSE ranges between 0.2893 psu and 0.3664 psu. For depths greater than 600 m, the average MAE ranges between 0.0565 psu and 0.0622 psu, and the average RMSE ranges between 0.0765 psu and 0.0838 psu. For salinity predictions from 16 to 20 days, the average MAE for depths less than 600 m ranges between 0.1930 psu and 0.2537 psu, and the average RMSE ranges between 0.3036 psu and 0.3797 psu. For depths greater than 600 m, the average MAE ranges between 0.0559 psu and 0.0562 psu, and the average RMSE ranges between 0.0580 psu and 0.0792 psu. Compared to the 1–5 day predictions, the average MAE and RMSE for salinity at depths less than 600 m increase by about 0.02, while the salinity error at depths greater than 600 m shows a negative growth. This is because, during the backpropagation process in deep learning, the model tends to focus on areas with larger losses, and the salinity error at depths greater than 600 m is relatively small. This causes the model to emphasize the salinity distribution at depths less than 600 m, leading to a lack of consistency.
To further analyze the characteristics of the predictions in each season, the temperature and salinity error data in Table 9 were plotted in 48 layers as shown in Figure 13 and Figure 14.
From Figure 13, it can be observed that during autumn (September–November), the temperature prediction accuracy is higher above the thermocline. However, the largest errors at the thermocline occur during autumn, indicating that autumn is a complex, variable, and unstable transition period [47,48]. Figure 14 shows that in salinity predictions, the errors in autumn are consistently larger at depths less than 600 m. Beyond 600 m, the model tends to focus on areas with larger errors during training, resulting in no significant differences in salinity distribution among the four seasons within this depth range. From the error trends, it can be seen that the temperature and salinity forecast results closely match the trend of the SimVP-gsta model’s inversion results, indicating that the fundamental model determines the basic performance of the three-dimensional forecasting model.
Based on the overall error analysis, it can be concluded that the prediction method proposed in this paper has strong robustness in capturing the patterns of temperature and salinity distribution. To more intuitively display the prediction results, the 48-layer prediction results for the 1st and 20th days were selected, and the 5th, 15th, 25th, 35th, and 45th layers were visualized. Additionally, the absolute error plot for each layer was drawn as shown in Figure 15 and Figure 16. From the figures, it can be seen that the proposed method can effectively predict the 3D temperature and salinity for 20 days, demonstrating strong reliability in the prediction results.

5.4.2. Northwest Pacific

According to Table 10, the average MAE for temperature predictions from 1 to 5 days ranges between 0.4777 °C and 0.5583 °C for depths less than 600 m, with the average RMSE ranging between 0.6896 °C and 0.8370 °C. For depths greater than 600 m, the average MAE ranges between 0.0836 °C and 0.0883 °C, with the average RMSE ranging between 0.1164 °C and 0.1226 °C. For temperature predictions from 16 to 20 days, the average MAE for depths less than 600 m ranges between 0.5378 °C and 0.6385 °C, with the average RMSE ranging between 0.7910 °C and 0.9503 °C. For depths greater than 600 m, the average MAE ranges between 0.0872 °C and 0.0904 °C, with the average RMSE ranging between 0.1216 °C and 0.1261 °C. Compared to the 1–5 day predictions, the average MAE and RMSE for temperature at depths less than 600 m both increase by about 0.14, while the average MAE and RMSE at depths greater than 600 m both increase by about 0.03.
The average MAE for salinity predictions from 1 to 5 days ranges between 0.1609 psu and 0.2088 psu for depths less than 600 m, with the average RMSE ranging between 0.2381 psu and 0.2996 psu. For depths greater than 600 m, the average MAE ranges between 0.0765 psu and 0.0934 psu, with the average RMSE ranging between 0.1229 psu and 0.1416 psu. For salinity predictions from 16 to 20 days, the average MAE for depths less than 600 m ranges between 0.1636 psu and 0.2097 psu, with the average RMSE ranging between 0.2417 psu and 0.3044 psu. For depths greater than 600 m, the average MAE ranges between 0.0780 psu and 0.0850 psu, with the average RMSE ranging between 0.1233 psu and 0.1357 psu. Compared to the 1–5 day predictions, the average MAE and RMSE for salinity at depths less than 600 m both increase by about 0.01.
To further analyze the characteristics of the predictions in each season, the temperature and salinity error data in Table 10 were plotted in 48 layers as shown in Figure 17 and Figure 18.
Based on Figure 17 and Figure 18, it can be observed that autumn exhibits the highest overall errors, further indicating its complexity and variability in temperature and salinity distribution. To provide a more intuitive display of the prediction results, the 48-layer prediction results for days 1 and 20 were selected. From these, the 5th, 15th, 25th, 35th, and 45th layers were chosen for visualization. Additionally, absolute error plots were generated for each layer as depicted in Figure 19 and Figure 20.
By validating the three-dimensional temperature and salinity models in two sea areas, the growth rate of temperature and salinity errors over 20 days and the characteristics of seasonal error distribution were analyzed. The construction effectiveness of the three-dimensional temperature and salinity forecasting model was confirmed through visualization and analysis of the results.

5.4.3. Error Analysis of Prediction

Due to the complex and variable nature of oceanic regions, errors in temperature and salinity predictions are influenced. Thus, we conducted a visual analysis of the prediction errors to explore the impact of dynamic regions on model learning. Taking the Northwest Pacific Ocean as an example, the prediction errors are illustrated in Figure 21, where erro_s represents salinity absolute error, and erro_t represents temperature absolute error. From the figure, it is evident that the errors are primarily concentrated along the Chinese coastline, near the Kuroshio Current, and near the equatorial currents. The coastal region errors are mainly attributed to the irregular shape of the coastline and the variable topography, which complicate the ocean dynamic processes. Additionally, significant tidal effects in these coastal areas make it challenging for the model to accurately capture temperature information distribution [49,50].
In the vicinity of the Kuroshio Current, both temperature and salinity exhibit substantial errors. This phenomenon may be due to the formation of numerous cyclonic and anticyclonic eddies in the Kuroshio Current region. These eddies cause variations in temperature and salinity distribution [51,52,53,54,55], creating interference with the model’s learning process. In the equatorial region, salinity errors are more pronounced, possibly due to the influence of equatorial currents that lead to an imbalance in salinity distribution patterns [56,57,58,59]. Temperature, being less affected by equatorial currents, shows errors mainly in the thermocline.

5.4.4. Sensitivity Analysis

Since the predicted SST and SLA serve as intermediate variables for predicting temperature and salinity information as well as sea surface input information, we further analyzed the sensitivity of temperature and salinity predictions to the accuracy of predicted SST and SLA. We designed the following methods for testing, as detailed in Table 11. Specifically, we substituted the true values for predicted values to analyze the impact on temperature and salinity predictions, thus examining how variations in input data affect subsurface temperature and salinity predictions. Method 1 uses true SST and true SLA as model inputs. Method 2 uses true SST and predicted SLA as inputs. Method 3 uses predicted SST and true SLA as inputs. Method 4 uses predicted SST and predicted SLA as inputs.
To demonstrate the feasibility of the validation method, we conducted experimental analysis in the offshore waters of China. The RMSE errors for the 20-day temperature and salinity predictions were plotted as shown in Figure 22 and Figure 23. Figure 22 presents the 20-day temperature forecast results, while Figure 23 displays the RMSE errors for salinity over the 20 days.
From Figure 22 and Figure 23, it can be analyzed that in the upper ocean (<50 m), SST and SLA have significant impacts on the prediction accuracy of subsurface temperatures. Among these, the accuracy of SST predictions has a major influence. When the SST is predicted more accurately, it reduces the prediction error of the upper ocean temperature. This phenomenon is due to the complex interactions between the upper ocean temperature and the SST through turbulent mixing, dynamic processes, and feedback mechanisms. There is a certain nonlinear mapping relationship between them, which we attempt to learn using deep learning methods. When the SST prediction is inaccurate, the deep learning model’s ability to fit the nonlinear relationship between the ocean surface and subsurface temperatures is impaired, leading to increased temperature prediction errors as SST prediction errors grow.
In the thermocline and halocline regions, the accuracy of SLA predictions significantly affects the error in the entire layer, and it is the dominant factor. The more accurate the SLA, the smaller the errors in temperature and salinity within the layer. This is likely due to the impact of currents, heat and salt transport, and eddies, which bring water with varying temperatures and salinities, affecting the structure of the temperature and salinity layers (reference needed). Accurate SLA predictions help the model capture these dynamic processes, improving the accuracy of temperature and salinity predictions in the layer.
Additionally, from the error plots, it can be seen that when the prediction period exceeds 10 days, the prediction errors for temperature and salinity increase significantly, especially for temperature. This is primarily because the prediction model’s accuracy deteriorates for later SST and SLA forecasts, leading to a growing deviation from the true values and a gradual decline in prediction performance.

6. Discussion

This paper utilizes the nonlinear relationship between SST, SLA, and the 3D sea temperature and salinity fields in the subsurface layer to transform the prediction of 3D subsurface temperature and salinity. The approach involves first capturing the horizontal movement of the SST and SLA, and then mapping these to the vertical distribution of subsurface temperature and salinity. The advantage of this method is that it does not require training the model with four-dimensional high-resolution spatiotemporal data, which significantly reduces computational burden and makes the model easier to train. Experimental results demonstrate that this method effectively predicts the 3D temperature and salinity fields of the subsurface ocean for up to 20 days.
To analyze the impact of input data variations on temperature and salinity changes at different depths, we designed sensitivity analysis experiments, substituting true values for predicted values to observe their effects on temperature and salinity predictions. Our analysis revealed that the accuracy of the SST effectively controls changes in the upper ocean temperatures, while the accuracy of the SLA controls the temperature and salinity information in the thermocline. The results of this experiment highlighted the following limitations of our proposed method and further analysis:
(1) Error Accumulation between Prediction and Inversion Models: If there is a bias between the SST and SLA data output by the prediction model and the true data, this bias affects the accuracy of the ocean mixed layer when input into the inversion model. Since the inversion model maps sea surface information to the subsurface temperature and salinity fields using a nonlinear function, greater accuracy in the inversion model leads to more reliable subsurface temperature and salinity information. Improving the performance of both the prediction and inversion models can effectively reduce error accumulation and enhance the accuracy of temperature and salinity predictions.
(2) Information Bias between Satellite Observations and Reanalysis Data: Satellite data, obtained from satellite observations, and reanalysis data, derived from numerical models after assimilation, may have inherent information biases. This bias introduces additional noise, which hampers the deep learning model’s ability to learn the mapping from sea surface information to temperature and salinity.
(3) Prediction Errors in Coastal, Estuarine, and High-Current Areas: Prediction errors are notably higher in coastal, estuarine, or high-current regions due to the complexity of terrain and strong current dynamics. Irregular coastlines and variable topography complicate oceanographic processes, affecting accurate predictions of flow and mixing. Additionally, tidal effects in coastal areas significantly impact temperature and salinity distribution, and turbulence and mixing in strong current regions introduce uncertainties in process parameterization, increasing prediction errors.
(4) Duration Limitations: The three-dimensional temperature and salinity prediction method is primarily constrained by the accuracy of SST and SLA forecasts. Sensitivity analysis shows that prediction errors for temperature and salinity in the upper ocean and thermocline increase significantly when the prediction period exceeds 10 days. This limitation is primarily due to the accuracy constraints of SST and SLA predictions.
Based on the issues identified with the current three-dimensional prediction method, future experiments will aim to address these problems by incorporating higher-precision satellite remote sensing data and more advanced deep learning models to further improve the accuracy of three-dimensional temperature and salinity field predictions.

7. Conclusions

This paper employs a method that combines two-dimensional models with multi-source ocean data to effectively predict three-dimensional temperature and salinity fields for 20 days. This approach replaces the previous method of directly constructing prediction models from three-dimensional data. The advantage of this method is that it avoids the challenges associated with stacking high-resolution three-dimensional data over time, which increases the computational resource demands and complicates model training. Analysis of the prediction results shows that this method aligns well with physical principles and accurately predicts most temperatures and salinities. Additionally, by selecting a more advanced deep learning model, we achieved more accurate predictions and conducted a refined analysis of the model’s results using various methods.
In summary, we have expanded the understanding of current three-dimensional temperature and salinity prediction methods, offering advancements in data utilization, model architecture, and prediction result analysis.

Author Contributions

Conceptualization, X.C., C.L. and S.Z.; methodology, X.C.; formal analysis, X.C., C.L., S.Z. and F.G.; funding acquisition, C.L. and F.G.; investigation, X.C., C.L., S.Z. and F.G.; resources, C.L., S.Z. and F.G.; software, X.C.; supervision, C.L., S.Z. and F.G.; validation, C.L., S.Z. and F.G.; visualization, X.C.; writing—original draft preparation, X.C.; writing—review and editing, X.C., C.L., S.Z. and F.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Key R&D Program Projects of Heilongjiang Province (No. 2022ZX01A15), the Cultivation Project of Qingdao Science and Technology Plan Park (No. 23-1-5-yqpy-11-qy).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The daily mean SST data provided by NOAA (https://www.ncei.noaa.gov/data, accessed on 10 October 2023) and daily mean SLA data provided by CMEMS (https://data.marine.copernicus.eu/products, accessed on 10 October 2023).The high-precision reanalysis data provided by CMEMS (https://data.marine.copernicus.eu/products, accessed on 15 October 2023).

Conflicts of Interest

Author Chang Liu was employed by the company Qingdao Hatran Ocean Intelligence Technology Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationship that could be construed as a potential conflict of interest.

References

  1. Tonani, M.; Balmaseda, M.; Bertino, L.; Blockley, E.; Brassington, G.; Davidson, F.; Drillet, Y.; Hogan, P.; Kuragano, T.; Lee, T.; et al. Status and future of global and regional ocean prediction systems. J. Oper. Oceanogr. 2015, 8, s201–s220. [Google Scholar] [CrossRef]
  2. Barceló-Llull, B.; Pascual, A.; Mason, E.; Mulet, S. Comparing a multivariate global ocean state estimate with high-resolution in situ data: An anticyclonic intrathermocline eddy near the Canary Islands. Front. Mar. Sci. 2018, 5, 66. [Google Scholar] [CrossRef]
  3. Guinehut, S.; Dhomps, A.L.; Larnicol, G.; Le Traon, P.Y. High resolution 3-D temperature and salinity fields derived from in situ and satellite observations. Ocean Sci. 2012, 8, 845–857. [Google Scholar] [CrossRef]
  4. Maes, C.; Dewitte, B.; Sudre, J.; Garçon, V.; Varillon, D. Small-scale features of temperature and salinity surface fields in the Coral Sea. J. Geophys. Res. Ocean. 2013, 118, 5426–5438. [Google Scholar] [CrossRef]
  5. Edwards, M.; Beaugrand, G.; Hays, G.C.; Koslow, J.A.; Richardson, A.J. Multi-decadal oceanic ecological datasets and their application in marine policy and management. Trends Ecol. Evol. 2010, 25, 602–610. [Google Scholar] [CrossRef] [PubMed]
  6. Fan, W.; Song, J.; Li, S. A numerical study on seasonal variations of the thermocline in the South China Sea based on the ROMS. Acta Oceanol. Sin. 2014, 33, 56–64. [Google Scholar] [CrossRef]
  7. Todd, A.; Zanna, L.; Couldrey, M.; Gregory, J.; Wu, Q.; Church, J.A.; Farneti, R.; Navarro-Labastida, R.; Lyu, K.; Saenko, O.; et al. Ocean-only FAFMIP: Understanding regional patterns of ocean heat content and dynamic sea level change. J. Adv. Model. Earth Syst. 2020, 12, e2019MS002027. [Google Scholar] [CrossRef]
  8. Morioka, Y.; Taguchi, B.; Behera, S.K. Eastward propagating decadal temperature variability in the South Atlantic and Indian Oceans. J. Geophys. Res. Ocean. 2017, 122, 5611–5623. [Google Scholar] [CrossRef]
  9. Zhang, Q.; Wang, H.; Dong, J.; Zhong, G.; Sun, X. Prediction of sea surface temperature using long short-term memory. IEEE Geosci. Remote Sens. Lett. 2017, 14, 1745–1749. [Google Scholar] [CrossRef]
  10. Yang, Y.; Dong, J.; Sun, X.; Lima, E.; Mu, Q.; Wang, X. A CFCC-LSTM model for sea surface temperature prediction. IEEE Geosci. Remote Sens. Lett. 2017, 15, 207–211. [Google Scholar] [CrossRef]
  11. Qiao, B.; Wu, Z.; Tang, Z.; Wu, G. Sea surface temperature prediction approach based on 3D CNN and LSTM with attention mechanism. In Proceedings of the 2022 24th International Conference on Advanced Communication Technology (ICACT), Pyeongchang, Republic of Korea, 13–16 February 2022; pp. 342–347. [Google Scholar]
  12. Guo, X.; He, J.; Wang, B.; Wu, J. Prediction of sea surface temperature by combining interdimensional and self-attention with neural networks. Remote Sens. 2022, 14, 4737. [Google Scholar] [CrossRef]
  13. Liang, S.; Zhao, A.; Qin, M.; Hu, L.; Wu, S.; Du, Z.; Liu, R. A Graph Memory Neural Network for Sea Surface Temperature Prediction. Remote Sens. 2023, 15, 3539. [Google Scholar] [CrossRef]
  14. Shao, Q.; Li, W.; Han, G.; Hou, G.; Liu, S.; Gong, Y.; Qu, P. A deep learning model for forecasting sea surface height anomalies and temperatures in the South China Sea. J. Geophys. Res. Ocean. 2021, 126, e2021JC017515. [Google Scholar] [CrossRef]
  15. Wei, L.; Guan, L. Seven-day sea surface temperature prediction using a 3DConv-LSTM model. Front. Mar. Sci. 2022, 9, 905848. [Google Scholar] [CrossRef]
  16. Hou, S.; Li, W.; Liu, T.; Zhou, S.; Guan, J.; Qin, R.; Wang, Z. MUST: A Multi-source Spatio-Temporal data fusion Model for short-term sea surface temperature prediction. Ocean Eng. 2022, 259, 111932. [Google Scholar] [CrossRef]
  17. Hu, S.; Shao, Q.; Li, W.; Han, G.; Zheng, Q.; Wang, R.; Liu, H. Multivariate Sea Surface Prediction in the Bohai Sea Using a Data-Driven Model. J. Mar. Sci. Eng. 2023, 11, 2096. [Google Scholar] [CrossRef]
  18. Wu, S.; Zhang, X.; Bao, S.; Dong, W.; Wang, S.; Li, X. Predicting Ocean Temperature in High-Frequency Internal Wave Area with Physics-Guided Deep Learning: A Case Study from the South China Sea. J. Mar. Sci. Eng. 2023, 11, 1728. [Google Scholar] [CrossRef]
  19. Zhang, L.; Zhang, Y.; Yin, X. Aquarius sea surface salinity retrieval in coastal regions based on deep neural networks. Remote Sens. Environ. 2023, 284, 113357. [Google Scholar] [CrossRef]
  20. Liu, S.; Jia, W.; Zhang, W. DSE-NN: Discretized Spatial Encoding Neural Network for Ocean Temperature and Salinity Interpolation in the North Atlantic. J. Mar. Sci. Eng. 2024, 12, 1013. [Google Scholar] [CrossRef]
  21. Xie, H.; Xu, Q.; Cheng, Y.; Yin, X.; Fan, K. Reconstructing three-dimensional salinity field of the South China Sea from satellite observations. Front. Mar. Sci. 2023, 10, 1168486. [Google Scholar] [CrossRef]
  22. Wen, Z.; Zhang, H.; Shu, W.; Zhang, L.; Liu, L.; Lu, X.; Zhou, Y.; Ren, J.; Li, S.; Zhang, Q. The SSR Brightness Temperature Increment Model Based on a Deep Neural Network. Remote Sens. 2023, 15, 4149. [Google Scholar] [CrossRef]
  23. Tian, T.; Cheng, L.; Wang, G.; Abraham, J.; Wei, W.; Ren, S.; Zhu, J.; Song, J.; Leng, H. Reconstructing ocean subsurface salinity at high resolution using a machine learning approach. Earth Syst. Sci. Data 2022, 14, 5037–5060. [Google Scholar] [CrossRef]
  24. Jang, E.; Kim, Y.J.; Im, J.; Park, Y.G.; Sung, T. Global sea surface salinity via the synergistic use of SMAP satellite and HYCOM data based on machine learning. Remote Sens. Environ. 2022, 273, 112980. [Google Scholar] [CrossRef]
  25. Jang, E.; Kim, Y.J.; Im, J.; Park, Y.G. Improvement of SMAP sea surface salinity in river-dominated oceans using machine learning approaches. GISci. Remote Sens. 2021, 58, 138–160. [Google Scholar] [CrossRef]
  26. Liu, J.; Zhang, T.; Han, G.; Gou, Y. TD-LSTM: Temporal dependence-based LSTM networks for marine temperature prediction. Sensors 2018, 18, 3797. [Google Scholar] [CrossRef] [PubMed]
  27. Zhang, K.; Geng, X.; Yan, X.H. Prediction of 3-D ocean temperature by multilayer convolutional LSTM. IEEE Geosci. Remote Sens. Lett. 2020, 17, 1303–1307. [Google Scholar] [CrossRef]
  28. Zuo, X.; Zhou, X.; Guo, D.; Li, S.; Liu, S.; Xu, C. Ocean temperature prediction based on stereo spatial and temporal 4-D convolution model. IEEE Geosci. Remote Sens. Lett. 2021, 19, 1003405. [Google Scholar] [CrossRef]
  29. Sun, N.; Zhou, Z.; Li, Q.; Zhou, X. Spatiotemporal prediction of monthly sea subsurface temperature fields using a 3D U-Net-Based model. Remote Sens. 2022, 14, 4890. [Google Scholar] [CrossRef]
  30. Muis, S.; Apecechea, M.I.; Dullaart, J.; de Lima Rego, J.; Madsen, K.S.; Su, J.; Yan, K.; Verlaan, M. A high-resolution global dataset of extreme sea levels, tides, and storm surges, including future projections. Front. Mar. Sci. 2020, 7, 263. [Google Scholar] [CrossRef]
  31. Banesh, D.; Petersen, M.R.; Ahrens, J.; Turton, T.L.; Samsel, F.; Schoonover, J.; Hamann, B. An image-based framework for ocean feature detection and analysis. J. Geovisualization Spat. Anal. 2021, 5, 17. [Google Scholar] [CrossRef]
  32. McWilliams, J.C. The nature and consequences of oceanic eddies. Geophys. Monogr. Ser. 2008, 177, 5–15. [Google Scholar]
  33. Klemas, V.; Yan, X.H. Subsurface and deeper ocean remote sensing from satellites: An overview and new results. Prog. Oceanogr. 2014, 122, 1–9. [Google Scholar] [CrossRef]
  34. Wu, X.; Yan, X.H.; Jo, Y.H.; Liu, W.T. Estimation of subsurface temperature anomaly in the North Atlantic using a self-organizing map neural network. J. Atmos. Ocean. Technol. 2012, 29, 1675–1688. [Google Scholar] [CrossRef]
  35. Cazenave, A.; Meyssignac, B.; Ablain, M.; Balmaseda, M.; Bamber, J.; Barletta, V.; Beckley, B.; Benveniste, J.; Berthier, E.; Blazquez, A.; et al. Global sea-level budget 1993-present. Earth Syst. Sci. Data 2018, 10, 1551–1590. [Google Scholar]
  36. Mao, K.; Liu, C.; Zhang, S.; Gao, F. Reconstructing Ocean Subsurface Temperature and Salinity from Sea Surface Information Based on Dual Path Convolutional Neural Networks. J. Mar. Sci. Eng. 2023, 11, 1030. [Google Scholar] [CrossRef]
  37. Zhao, Y.; He, Z.; Jiang, Y. The Inversion of Three-Dimensional Ocean Temperature and Salinity Fields for the Assimilation of Satellite Observations. J. Mar. Sci. Eng. 2024, 12, 534. [Google Scholar] [CrossRef]
  38. Xie, H.; Xu, Q.; Cheng, Y.; Yin, X.; Jia, Y. Reconstruction of subsurface temperature field in the south China Sea from satellite observations based on an attention U-net model. IEEE Trans. Geosci. Remote Sens. 2022, 60, 4209319. [Google Scholar] [CrossRef]
  39. Zhang, J.; Ning, P.; Zhang, X.; Wang, X.; Zhang, A. Deriving Sea Subsurface Temperature Fields from Satellite Remote Sensing Data Using a Generative Adversarial Network Model. Earth Space Sci. 2022, 10, e2022EA002804. [Google Scholar] [CrossRef]
  40. Song, T.; Xu, G.; Yang, K.; Li, X.; Peng, S. Convformer: A Model for Reconstructing Ocean Subsurface Temperature and Salinity Fields Based on Multi-Source Remote Sensing Observations. Remote Sens. 2024, 16, 2422. [Google Scholar] [CrossRef]
  41. Ronneberger, O.; Fischer, P.; Brox, T. U-net: Convolutional networks for biomedical image segmentation. In Medical Image Computing and Computer-Assisted Intervention–MICCAI 2015: 18th International Conference, Munich, Germany, 5–9 October 2015; Proceedings, Part III 18; Springer: Cham, Switzerland, 2015; pp. 234–241. [Google Scholar]
  42. Gao, Z.; Tan, C.; Wu, L.; Li, S.Z. Simvp: Simpler yet better video prediction. In Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition, New Orleans, LA, USA, 18–24 June 2022; pp. 3170–3180. [Google Scholar]
  43. Mehta, S.; Rastegari, M. Mobilevit: Light-weight, general-purpose, and mobile-friendly vision transformer. arXiv 2021, arXiv:2110.02178. [Google Scholar]
  44. Aguilar-Martinez, S.; Hsieh, W.W. Forecasts of tropical Pacific sea surface temperatures by neural networks and support vector regression. Int. J. Oceanogr. 2009, 2009, 167239. [Google Scholar] [CrossRef]
  45. Arnab, A.; Dehghani, M.; Heigold, G.; Sun, C.; Lučić, M.; Schmid, C. Vivit: A video vision transformer. In Proceedings of the IEEE/CVF International Conference on Computer Vision, Montreal, QC, Canada, 10–17 October 2021; pp. 6836–6846. [Google Scholar]
  46. Li, H.; Qian, J.; Tian, Y.; Rakhlin, A.; Jadbabaie, A. Convex and non-convex optimization under generalized smoothness. Adv. Neural Inf. Process. Syst. 2024, 36. [Google Scholar]
  47. Jamshidi, S.; Abu Bakar, N. Seasonal variations in temperature, salinity and density in the southern coastal waters of the Caspian Sea. Oceanology 2012, 52, 380–396. [Google Scholar] [CrossRef]
  48. Pennington, J.T.; Chavez, F.P. Seasonal fluctuations of temperature, salinity, nitrate, chlorophyll and primary production at station H3/M1 over 1989–1996 in Monterey Bay, California. Deep Sea Res. Part II Top. Stud. Oceanogr. 2000, 47, 947–973. [Google Scholar] [CrossRef]
  49. Hu, H.; Wang, J. Modeling effects of tidal and wave mixing on circulation and thermohaline structures in the Bering Sea: Process studies. J. Geophys. Res. Ocean. 2010, 115. [Google Scholar] [CrossRef]
  50. Yndestad, H.; Turrell, W.R.; Ozhigin, V. Lunar nodal tide effects on variability of sea level, temperature, and salinity in the Faroe-Shetland Channel and the Barents Sea. Deep Sea Res. Part I Oceanogr. Res. Pap. 2008, 55, 1201–1217. [Google Scholar] [CrossRef]
  51. He, Q.; Zhan, H.; Cai, S.; He, Y.; Huang, G.; Zhan, W. A new assessment of mesoscale eddies in the South China Sea: Surface features, three-dimensional structures, and thermohaline transports. J. Geophys. Res. Ocean. 2018, 123, 4906–4929. [Google Scholar] [CrossRef]
  52. Sun, W.; An, M.; Liu, J.; Liu, J.; Yang, J.; Tan, W.; Dong, C.; Liu, Y. Comparative analysis of four types of mesoscale eddies in the Kuroshio-Oyashio extension region. Front. Mar. Sci. 2022, 9, 984244. [Google Scholar] [CrossRef]
  53. Yao, H.; Ma, C.; Jing, Z.; Zhang, Z. On the Vertical Structure of Mesoscale Eddies in the Kuroshio-Oyashio Extension. Geophys. Res. Lett. 2023, 50, e2023GL105642. [Google Scholar] [CrossRef]
  54. Seidov, D.; Mishonov, A.; Reagan, J.; Parsons, R. Eddy-resolving in situ ocean climatologies of temperature and salinity in the Northwest Atlantic Ocean. J. Geophys. Res. Ocean. 2019, 124, 41–58. [Google Scholar] [CrossRef]
  55. Liu, S.S.; Sun, L.; Wu, Q.; Yang, Y.J. The responses of cyclonic and anticyclonic eddies to typhoon forcing: The vertical temperature-salinity structure changes associated with the horizontal convergence/divergence. J. Geophys. Res. Ocean. 2017, 122, 4974–4989. [Google Scholar] [CrossRef]
  56. Chen, L.; Zhang, R.H.; Gao, C. Effects of temperature and salinity on surface currents in the equatorial pacific. J. Geophys. Res. Ocean. 2022, 127, e2021JC018175. [Google Scholar] [CrossRef]
  57. Vaughan, S.L.; Molinari, R.L. Temperature and salinity variability in the deep western boundary current. J. Phys. Oceanogr. 1997, 27, 749–761. [Google Scholar] [CrossRef]
  58. Lynn, R.J. Seasonal variation of temperature and salinity at 10 m in the California Current. Calif. Coop. Ocean. Fish. Invest. Rep 1967, 11, 31. [Google Scholar]
  59. Johnson, G.C.; McPhaden, M.J.; Rowe, G.D.; McTaggart, K.E. Upper equatorial Pacific Ocean current and salinity variability during the 1996–1998 El Nino–La Nina cycle. J. Geophys. Res. Ocean. 2000, 105, 1037–1053. [Google Scholar] [CrossRef]
Figure 1. Data partitioning for prediction and inversion. Figure (a) shows the division of forecast data, and figure (b) shows the division of inversion data.
Figure 1. Data partitioning for prediction and inversion. Figure (a) shows the division of forecast data, and figure (b) shows the division of inversion data.
Jmse 12 01396 g001
Figure 2. SimVP-gsta model. Figure (a) shows the SimVP–gsta model, and figure (b) shows the gsta module.
Figure 2. SimVP-gsta model. Figure (a) shows the SimVP–gsta model, and figure (b) shows the gsta module.
Jmse 12 01396 g002
Figure 3. M-ViT model. Figure (a) shows the M–ViT model, and figure (b) shows the Mobile–ViT module.
Figure 3. M-ViT model. Figure (a) shows the M–ViT model, and figure (b) shows the Mobile–ViT module.
Jmse 12 01396 g003
Figure 6. Loss function settings in the inversion experiment.
Figure 6. Loss function settings in the inversion experiment.
Jmse 12 01396 g006
Figure 7. Visualization of the average SST prediction error over 20 days in the Coastal Waters of China. The MAE chart is on the left, and the RMSE chart is on the right.
Figure 7. Visualization of the average SST prediction error over 20 days in the Coastal Waters of China. The MAE chart is on the left, and the RMSE chart is on the right.
Jmse 12 01396 g007
Figure 8. Visualization of the average SST prediction error over 20 days in the Northwest Pacific. The MAE chart is on the left, and the RMSE chart is on the right.
Figure 8. Visualization of the average SST prediction error over 20 days in the Northwest Pacific. The MAE chart is on the left, and the RMSE chart is on the right.
Jmse 12 01396 g008
Figure 9. Visualization of the average SLA prediction error over 20 days in the Coastal Waters of China. The MAE chart is on the left, and the RMSE chart is on the right.
Figure 9. Visualization of the average SLA prediction error over 20 days in the Coastal Waters of China. The MAE chart is on the left, and the RMSE chart is on the right.
Jmse 12 01396 g009
Figure 10. Visualization of the average SLA prediction error over 20 days in the Northwest Pacific. The MAE chart is on the left, and the RMSE chart is on the right.
Figure 10. Visualization of the average SLA prediction error over 20 days in the Northwest Pacific. The MAE chart is on the left, and the RMSE chart is on the right.
Jmse 12 01396 g010
Figure 11. Visualization of errors in 48-layer temperature and salinity inversion in the Coastal Waters of China. Figure (a) represents the temperature error, with the MAE chart on the left and the RMSE chart on the right. Figure (b) represents the salinity error, with the MAE chart on the left and the RMSE chart on the right.
Figure 11. Visualization of errors in 48-layer temperature and salinity inversion in the Coastal Waters of China. Figure (a) represents the temperature error, with the MAE chart on the left and the RMSE chart on the right. Figure (b) represents the salinity error, with the MAE chart on the left and the RMSE chart on the right.
Jmse 12 01396 g011
Figure 12. Visualization of errors in 48-layer temperature and salinity inversion in the Northwest Pacific. Figure (a) represents the temperature error, with the MAE chart on the left and the RMSE chart on the right. Figures (b) represents the salinity error, with the MAE chart on the left and the RMSE chart on the right.
Figure 12. Visualization of errors in 48-layer temperature and salinity inversion in the Northwest Pacific. Figure (a) represents the temperature error, with the MAE chart on the left and the RMSE chart on the right. Figures (b) represents the salinity error, with the MAE chart on the left and the RMSE chart on the right.
Jmse 12 01396 g012
Figure 13. Three-dimensional sea temperature prediction error curves in the Coastal Waters of China. Figures (ad) show the average MAE for temperature every 5 days from 1 to 20 days, and figures (eh) show the average RMSE for temperature.
Figure 13. Three-dimensional sea temperature prediction error curves in the Coastal Waters of China. Figures (ad) show the average MAE for temperature every 5 days from 1 to 20 days, and figures (eh) show the average RMSE for temperature.
Jmse 12 01396 g013
Figure 14. Three-dimensional salinity prediction error curves in the Coastal Waters of China. Figures (ad) show the average MAE for salinity, and figures (eh) show the average RMSE for salinity.
Figure 14. Three-dimensional salinity prediction error curves in the Coastal Waters of China. Figures (ad) show the average MAE for salinity, and figures (eh) show the average RMSE for salinity.
Jmse 12 01396 g014
Figure 15. Visualization of the 3D temperature and salinity predictions for day 1 in the Coastal Waters of China. (a) Temperature, (b) salinity.
Figure 15. Visualization of the 3D temperature and salinity predictions for day 1 in the Coastal Waters of China. (a) Temperature, (b) salinity.
Jmse 12 01396 g015
Figure 16. Visualization of the 3D temperature and salinity predictions for day 20 in the Coastal Waters of China. (a) Temperature, (b) salinity.
Figure 16. Visualization of the 3D temperature and salinity predictions for day 20 in the Coastal Waters of China. (a) Temperature, (b) salinity.
Jmse 12 01396 g016
Figure 17. Three-dimensional sea temperature prediction error curves in the Northwest Pacific. Figures (ad) show the average MAE for temperature every 5 days from 1 to 20 days, and figures (eh) show the average RMSE for temperature.
Figure 17. Three-dimensional sea temperature prediction error curves in the Northwest Pacific. Figures (ad) show the average MAE for temperature every 5 days from 1 to 20 days, and figures (eh) show the average RMSE for temperature.
Jmse 12 01396 g017
Figure 18. Three-dimensional sea salinity prediction error curves in the Northwest Pacific. Figures (ad) show the average MAE for salinity, and figures (eh) show the average RMSE for salinity.
Figure 18. Three-dimensional sea salinity prediction error curves in the Northwest Pacific. Figures (ad) show the average MAE for salinity, and figures (eh) show the average RMSE for salinity.
Jmse 12 01396 g018
Figure 19. Visualization of the 3D temperature and salinity predictions for day 1 in the Northwest Pacific. (a) Temperature, (b) salinity.
Figure 19. Visualization of the 3D temperature and salinity predictions for day 1 in the Northwest Pacific. (a) Temperature, (b) salinity.
Jmse 12 01396 g019
Figure 20. Visualization of the 3D temperature and salinity predictions for day 20 in the Northwest Pacific. (a) Temperature, (b) salinity.
Figure 20. Visualization of the 3D temperature and salinity predictions for day 20 in the Northwest Pacific. (a) Temperature, (b) salinity.
Jmse 12 01396 g020
Figure 21. Visualization of temperature and salinity prediction errors. The areas highlighted in the figure are the regions with larger forecast errors in this instance.
Figure 21. Visualization of temperature and salinity prediction errors. The areas highlighted in the figure are the regions with larger forecast errors in this instance.
Jmse 12 01396 g021
Figure 22. The RMSE curves for temperature with different methods. Figures (ad) show the errors averaged every 5 days for the 20-day temperature and salinity forecasts. Orange and red box indicate the positions where error changes are noticeably variable over time.
Figure 22. The RMSE curves for temperature with different methods. Figures (ad) show the errors averaged every 5 days for the 20-day temperature and salinity forecasts. Orange and red box indicate the positions where error changes are noticeably variable over time.
Jmse 12 01396 g022
Figure 23. The RMSE curves for salinity with different methods. Figures (ad) show the errors averaged every 5 days for the 20-day temperature and salinity forecasts. Orange and red box indicate the positions where error changes are noticeably variable over time.
Figure 23. The RMSE curves for salinity with different methods. Figures (ad) show the errors averaged every 5 days for the 20-day temperature and salinity forecasts. Orange and red box indicate the positions where error changes are noticeably variable over time.
Jmse 12 01396 g023
Table 1. Table of depth statistics for each layer of reanalysis data.
Table 1. Table of depth statistics for each layer of reanalysis data.
LayerDepthLayerDepthLayerDepthLayerDepth
11 m1325 m25185 m371452 m
22 m1429 m26222 m381684 m
33 m1534 m27216 m391941 m
45 m1640 m28318 m402225 m
56 m1747 m29380 m412533 m
67 m1855 m30453 m422865 m
79 m1965 m31541 m433220 m
811 m2077 m32643 m443597 m
913 m2192 m33763 m453992 m
1015 m22109 m34902 m464405 m
1118 m23130 m351062 m474833 m
1221 m24155 m361245 m485274 m
Table 2. Data and model parameter settings.
Table 2. Data and model parameter settings.
ParameterModelData
SimVP-gsta M-ViTBT(in, out)L(in, out)E(in, out)
dim N_T N_S CLP
NWP5126610246(7, 5)8(30, 20)(1, 48)(2, 2)
CWOC5126810246(15, 9)8(30, 20)(1, 48)(2, 2)
Table 3. Statistical results of the 20-day SST prediction error in the Coastal Waters of China. An arrow pointing downward indicates that a lower value of the indicator is better.
Table 3. Statistical results of the 20-day SST prediction error in the Coastal Waters of China. An arrow pointing downward indicates that a lower value of the indicator is better.
TimeAverage-MAE (↓)Average-RMSE (↓)
Simvp-gsta UNet M-ViT L-R ConvFormer ViT Simvp-gsta UNet M-ViT L-R ConvFormer ViT
1∼5 days0.38960.68060.32010.34570.43480.38430.53650.82990.44710.47540.57490.5292
6∼10 days0.49340.75850.46390.53860.55830.48940.67150.94960.64230.75250.74500.6683
11∼15 days0.53440.80930.52730.68400.65450.58190.72771.01820.72820.95050.85370.7720
16∼20 days0.58290.85540.55410.82230.65490.64540.79171.07640.74841.14110.87280.8511
Average0.50010.77600.46640.59770.57560.52530.68190.96850.64650.82990.76160.7052
Table 4. Statistical results of the 20-day SST prediction error in the Northwest Pacific. An arrow pointing downward indicates that a lower value of the indicator is better.
Table 4. Statistical results of the 20-day SST prediction error in the Northwest Pacific. An arrow pointing downward indicates that a lower value of the indicator is better.
TimeAverage-MAE (↓)Average-RMSE (↓)
Simvp-gsta UNet M-ViT L-R ConvFormer ViT Simvp-gsta UNet M-ViT L-R ConvFormer ViT
1∼5 days0.29720.50270.27820.29890.67300.32860.43680.65380.47930.45290.84410.4902
6∼10 days0.43470.58290.41820.49600.74630.44460.64410.78100.62790.74320.97370.6556
11∼15 days0.50990.64220.48840.62970.78990.50930.75370.86810.73140.94081.05420.7467
16∼20 days0.56980.70340.54250.76030.84380.56420.84520.94910.81561.13951.10620.8261
Average0.45290.60780.43180.54620.75580.46170.67000.81300.64860.81910.99460.6797
Table 5. Statistical results of the 20-day SLA prediction error in the Coastal Waters of China. An arrow pointing downward indicates that a lower value of the indicator is better.
Table 5. Statistical results of the 20-day SLA prediction error in the Coastal Waters of China. An arrow pointing downward indicates that a lower value of the indicator is better.
TimeAverage-MAE (↓)Average-RMSE (↓)
Simvp-gsta UNet M-ViT L-R ConvFormer ViT Simvp-gsta UNet M-ViT L-R ConvFormer ViT
1∼5 days0.02290.04290.02240.03270.03270.02810.03200.05460.03170.03850.04310.0378
6∼10 days0.03680.05230.03630.05060.04550.04410.05100.06740.05070.06870.06010.0587
11∼15 days0.04410.05890.04310.06190.05240.05480.06020.07610.05930.08160.06920.0715
16∼20 days0.04840.06470.04760.06940.05630.06140.06580.08290.06520.09030.07450.0797
Average0.03810.05470.03740.05260.04670.04710.05230.07030.05170.06980.06170.0619
Table 6. Statistical results of the 20-day SLA prediction error in the Northwest Pacific. An arrow pointing downward indicates that a lower value of the indicator is better.
Table 6. Statistical results of the 20-day SLA prediction error in the Northwest Pacific. An arrow pointing downward indicates that a lower value of the indicator is better.
TimeAverage-MAE (↓)Average-RMSE (↓)
Simvp-gsta UNet M-ViT L-R ConvFormer ViT Simvp-gsta UNet M-ViT L-R ConvFormer ViT
1∼5 days0.02130.04020.02350.02140.05070.02030.03030.05040.03260.03540.05440.0300
6∼10 days0.03270.04660.03270.03960.05070.03290.04700.06080.04640.06360.06480.0467
11∼15 days0.04160.05270.04080.05210.05680.04290.06060.07140.05870.08130.07560.0609
16∼20 days0.04760.05780.04670.06070.06110.05030.07020.08010.06800.09320.08400.0715
Average0.03510.04930.03590.04350.05480.03660.05200.06570.05140.06840.06970.0523
Table 7. Statistical results of vertical temperature and salinity inversion errors in the Coastal Waters of China. An arrow pointing downward indicates that a lower value of the indicator is better.
Table 7. Statistical results of vertical temperature and salinity inversion errors in the Coastal Waters of China. An arrow pointing downward indicates that a lower value of the indicator is better.
ElementDepthAverage-MAE (↓)Average-RMSE (↓)
Simvp-gsta UNet M-ViT Atten-UNet ConvFormer ViT Simvp-gsta UNet M-ViT Atten-UNet ConvFormer ViT
Temp<600 m0.49690.75770.50990.68070.79330.53910.69681.09160.71920.99861.03850.7881
>600 m0.10170.16050.10860.15650.15500.13680.14670.26860.15350.26780.21710.1929
Sal<600 m0.21060.41830.23930.35120.37250.32820.32610.94220.38650.88160.55160.6021
>600 m0.06900.45860.11620.22950.37140.23610.09371.00960.30630.95800.51760.4784
Table 8. Statistical results of vertical temperature and salinity inversion errors in the Northwest Pacific. An arrow pointing downward indicates that a lower value of the indicator is better.
Table 8. Statistical results of vertical temperature and salinity inversion errors in the Northwest Pacific. An arrow pointing downward indicates that a lower value of the indicator is better.
ElementDepthAverage-MAE (↓)Average-RMSE (↓)
Simvp-gsta UNet M-ViT Atten-UNet ConvFormer ViT Simvp-gsta UNet M-ViT Atten-UNet ConvFormer ViT
Temp<600 m0.51340.66840.57610.78730.82080.88230.75870.95670.84381.12411.16981.3008
>600 m0.09750.12010.11010.17850.15850.17350.13770.17860.15900.30030.25030.2826
Sal<600 m0.19350.26610.24400.40370.39600.39350.31920.36820.42440.75960.66070.6035
>600 m0.10250.13110.23810.32080.25850.33570.24430.25270.46470.70500.54420.6079
Table 9. Statistical results of the 20-day 3D temperature and salinity prediction errors in the Coastal Waters of China. An arrow pointing downward indicates that a lower value of the indicator is better.
Table 9. Statistical results of the 20-day 3D temperature and salinity prediction errors in the Coastal Waters of China. An arrow pointing downward indicates that a lower value of the indicator is better.
TimeElementDepthAverage-MAE (↓)Average-RMSE (↓)
Winter Spring Summer Autumn Winter Spring Summer Autumn
1∼5 daysTemp<600 m0.56970.54630.60020.49500.73920.74200.71350.6989
>600 m0.10880.11070.10450.10430.15290.15620.15130.1536
Sal<600 m0.20180.18970.19450.24850.33580.29980.28930.3664
>600 m0.06090.06220.05870.05650.08290.08380.07910.0765
6∼10 daysTemp<600 m0.58180.57450.52310.51600.76110.77820.74320.7269
>600 m0.11060.11040.10420.10520.15540.15610.15110.1551
Sal<600 m0.20470.19160.19670.24790.33920.30240.29580.3711
>600 m0.05840.05990.05670.05640.07900.08100.07660.0764
11∼15 daysTemp<600 m0.59260.59540.54540.53330.78360.80420.77430.7495
>600 m0.11240.11070.10440.10610.15770.15670.15150.1564
Sal<600 m0.20750.19240.20100.24960.34200.30260.30350.3754
>600 m0.05810.05920.05600.05550.07860.08020.07580.0753
16∼20 daysTemp<600 m0.61160.61360.56600.55310.81350.82250.80270.7775
>600 m0.11430.12250.10510.10620.15980.17110.15260.1565
Sal<600 m0.21040.19300.20680.25370.34610.30360.31170.3797
>600 m0.05820.05850.05590.05620.07900.07920.07580.0763
Table 10. Statistical results of the 20-day 3D temperature and salinity prediction errors in the Northwest Pacific. An arrow pointing downward indicates that a lower value of the indicator is better.
Table 10. Statistical results of the 20-day 3D temperature and salinity prediction errors in the Northwest Pacific. An arrow pointing downward indicates that a lower value of the indicator is better.
TimeElementDepthAverage-MAE (↓)Average-RMSE (↓)
Winter Spring Summer Autumn Winter Spring Summer Autumn
1∼5 daysTemp<600 m0.48340.47770.55320.55830.70020.68960.81970.8370
>600 m0.08420.08360.08500.08830.11650.11640.11740.1226
Sal<600 m0.17590.16090.18720.20880.26000.23810.26680.2996
>600 m0.08200.07650.07960.09340.12800.12520.12290.1416
6∼10 daysTemp<600 m0.49800.50420.57870.58080.72640.72410.85390.8653
>600 m0.08500.08440.08540.08870.11780.11770.11800.1231
Sal<600 m0.17460.16090.18880.20810.25900.23840.27040.3002
>600 m0.08270.07650.07920.08960.12860.12460.12540.1376
11∼15 daysTemp<600 m0.51760.53600.60830.60590.75700.77020.89970.9013
>600 m0.08620.08580.08630.08950.11990.12000.11950.1244
Sal<600 m0.17320.16180.19130.20850.25800.23920.27450.3016
>600 m0.08180.07680.07850.08670.12830.12240.12490.1338
16∼20 daysTemp<600 m0.53780.56550.63850.63100.79100.81460.95030.9429
>600 m0.08760.08720.08750.09040.12250.12290.12160.1261
Sal<600 m0.17240.16360.19410.20970.25760.24170.27860.3044
>600 m0.08060.07800.07830.08500.12700.12490.12330.1357
Table 11. Table of model sensitivity testing method. The charkmark indicates that this data is selected as the input data.
Table 11. Table of model sensitivity testing method. The charkmark indicates that this data is selected as the input data.
InputsMethod 1Method 2Method 3Method 4
SST-True
SST-Pred
SLA-True
SLA-Pred
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Cao, X.; Liu, C.; Zhang, S.; Gao, F. A Method for Predicting High-Resolution 3D Variations in Temperature and Salinity Fields Using Multi-Source Ocean Data. J. Mar. Sci. Eng. 2024, 12, 1396. https://doi.org/10.3390/jmse12081396

AMA Style

Cao X, Liu C, Zhang S, Gao F. A Method for Predicting High-Resolution 3D Variations in Temperature and Salinity Fields Using Multi-Source Ocean Data. Journal of Marine Science and Engineering. 2024; 12(8):1396. https://doi.org/10.3390/jmse12081396

Chicago/Turabian Style

Cao, Xiaohu, Chang Liu, Shaoqing Zhang, and Feng Gao. 2024. "A Method for Predicting High-Resolution 3D Variations in Temperature and Salinity Fields Using Multi-Source Ocean Data" Journal of Marine Science and Engineering 12, no. 8: 1396. https://doi.org/10.3390/jmse12081396

APA Style

Cao, X., Liu, C., Zhang, S., & Gao, F. (2024). A Method for Predicting High-Resolution 3D Variations in Temperature and Salinity Fields Using Multi-Source Ocean Data. Journal of Marine Science and Engineering, 12(8), 1396. https://doi.org/10.3390/jmse12081396

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop