The ability to monitor and predict sea temperature is crucial for determining the likelihood that ocean-related events will occur. However, most studies have focused on predicting sea surface temperature, and less attention has been paid to predicting sea subsurface temperature (SSbT), which can reflect the thermal state of the entire ocean. In this study, we use a 3D U-Net model to predict the SSbT in the upper 400 m of the Pacific Ocean and its adjacent oceans for lead times of 12 months. Two reconstructed SSbT products are added to the training set to solve the problem of insufficient observation data. Experimental results indicate that this method can predict the ocean temperature more accurately than previous methods in most depth layers. The root mean square error and mean absolute error of the predicted SSbT fields for all lead times are within 0.5–0.7 °C and 0.3–0.45 °C, respectively, while the average correlation coefficient scores of the predicted SSbT profiles are above 0.96 for almost all lead times. In addition, a case study qualitatively demonstrates that the 3D U-Net model can predict realistic SSbT variations in the study area and, thus, facilitate understanding of future changes in the thermal state of the subsurface ocean.

Sea surface temperature (SST) is an important factor in the marine environment and has significant impacts on climate, ecology, and maritime activities. Most existing SST prediction methods consider the ocean as a uniform field and use a uniform grid to predict SST. However, the marine environment is a complex system, and factors such as solar radiation, differences in land and sea thermal properties, and ocean circulation lead to uneven spatial distributions of SSTs. We propose a non-uniform grid construction method based on an SST spatial gradient to encode SST data, as well as a Non-uniform Grid Graph Convolutional Network (NGGCN) model. The NGGCN consists of two spatiotemporal modules, each of which extracts spatial features from the GCN module, captures temporal correlations through the GRU module, and performs feature restoration and output results through the fully connected module. We selected data from the Yellow Sea and Bohai Sea to validate the effectiveness of the NGGCN in predicting SST at different time scales and prediction steps. The results indicate that our model shows a significant improvement in prediction performance compared to other models.

Directional spreading function of the gravity-capillary wave spectrum can provide the high-wavenumber wave energy distribution among different directions on the sea surface. The existing directional spreading functions have been mainly developed for the low-wavenumber gravity wave with buoy data. In this paper, we use radar observations to derive the directional spreading function of the gravity-capillary wave spectrum, which is expressed as the second-order Fourier series expansion. So far the standard form of the second-order harmonic coefficient has not been proposed to correctly unify the gravity and gravity-capillary wave. Our strategy is to introduce a correcting term to replace the inaccurate gravity-capillary spectral component in Elfouhaily’s directional spreading function. The second-order harmonic coefficient at L, C and Ku band calculated by the radar observation is used to fit the correcting term to obtain one at the full gravity-capillary wave region. According to our proposed the directional spreading function, there is a spectral region between the gravity and gravity-capillary range where it signifies the negative upwind–crosswind asymmetry at low and moderate speed range. And this is not reflected by the previous models, but has been confirmed by radar observations. The Root Mean Square Difference (RMSD) of the proposed second-order harmonic coefficient versus the radar-observed one at L, C band Ku band is 0.0438, 0.0263 and 0.0382, respectively. The overall bias and RMSD are −0.0029 and 0.0433 for the whole second-order harmonic coefficient range, respectively. The result verifies the accuracy of the proposed directional spreading function at L, C band Ku band.