Next Article in Journal
Internal Waves Study on a Narrow Steep Shelf of the Black Sea Using the Spatial Antenna of Line Temperature Sensors
Previous Article in Journal
A Low-Cost Stereo Video System for Measuring Directional Wind Waves
Previous Article in Special Issue
Formation, Detection, and Modeling of Submerged Oil: A Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Application of Three-Dimensional Interpolation in Estimating Diapycnal Diffusivity in the South China Sea

1
Physical Oceanography Laboratory/CIMST, Ocean University of China, and Qingdao National Laboratory for Marine Science and Technology, Qingdao 266100, China
2
Ocean College, Zhejiang University, Zhoushan 316021, China
*
Authors to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2020, 8(11), 832; https://doi.org/10.3390/jmse8110832
Submission received: 28 September 2020 / Revised: 16 October 2020 / Accepted: 20 October 2020 / Published: 22 October 2020

Abstract

:
Diapycnal diffusivity is an important parameter to characterize oceanic turbulent mixing and vertical transport. However, due to the challenging accessibility of field observations, the observation of diapycnal diffusivity in the South China Sea (SCS) is rare. In this study, a three-dimensional field of diapycnal diffusivity in the SCS with high spatial resolution is performed by interpolating the rare field observations, which aims to provide a reference for the value of diapycnal diffusivity in ocean models. Given the anisotropy of diapycnal diffusivity and its rapid change in the magnitude in the vertical direction, several typical interpolation methods are compared in this study. Results of two cross-validation methods demonstrate that the three-dimensional (3D) thin-plate spline interpolation method yields the most reasonable and accurate results among a total of five typical methods used in this study.

1. Introduction

Turbulent mixing is a microscale process that makes the properties of water bodies more uniform, which is extremely important for regulating climate and marine productivity [1,2]. Moreover, turbulent mixing could be an important mechanism for controlling the vertical transport of nutrients, dissolved gases and particulate in water [3]. The South China Sea (SCS) is the largest marginal sea of the Northwest Pacific and the mixing in the SCS plays a key role in maintaining the strength of abyssal water transport and the Pacific circulation [3,4].
Turbulent mixing in the SCS has been a hot topic for scholars for a long time. Oceanic turbulent mixing is commonly characterized by turbulent kinetic energy dissipation rate and diapycnal diffusivity, and there is a strong correlation between these two indicators [5]. Tian et al. [3] found that the diapycnal diffusivity in the deep ocean basin of the SCS was two orders larger than that in the western Pacific, based on a fine scale parameterization and hydrographic observations obtained from the northern SCS and the western Pacific. Using the Thorpe scale method, Alford et al. [6] and Buijsman et al. [7] revealed that the diapycnal diffusivity of the Luzon Strait increases in the deep water. Lozovatsky et al. [8] studied the spatial structure and temporal variability of turbulent kinetic energy dissipation rate in limited areas of the northern SCS based on the microstructure measurements. Measurements and numerical simulations [9,10] indicated that energetic internal tides and internal waves generated near the Luzon Strait propagate into the SCS and enhance turbulent mixing. Yang et al. [4] presented the first three-dimensional (3D) patchy distribution of turbulent mixing in the SCS based on the Gregg–Henyey–Polzin (GHP) parameterization and hydrological observations from 2005 to 2012. Shang et al. [11] analyzed spatial distribution of turbulent mixing in the upper ocean of the SCS with microstructure data obtained from 26 April to 23 May 2010.
These studies have greatly improved our understanding of turbulent mixing in the SCS. However, due to the limitations of field observations, diapycnal diffusivity in the SCS calculated by the GHP parameterization method has only patchy spatial distribution. And the low spatial resolution is detrimental to the further research of turbulent mixing in the SCS. Unfortunately, turbulent mixing is a complicated process and far from fully understood. It is now often used to add a constant background diffusivity to the coefficients of vertical mixing to represent mixing processes [12], which may be a potential source of error in the numerical simulation results (diapycnal diffusivity in the SCS varies many orders of magnitude and is extremely inhomogeneous).
Due to lack of understanding of the turbulent mixing variability in the whole SCS, the products from numerical simulation based on a simple parameterization of eddy diffusivity should be treated with caution [4]. The spatial interpolation of scattered data has been focused on for many years in various fields. The advantage of the spatial interpolation method is that the 3D field of diapycnal diffusivity can be obtained with low computational cost and high spatial resolution. Against this background, it’s a method worth trying to estimate diapycnal diffusivity by interpolating the measurements. As the measured diapycnal diffusivity is sparsely distributed in the SCS, a relatively accurate interpolation method would be useful. Furthermore, the interpolated results may provide a reference for calculating vertical fluxes across the basin.
There are several methods widely used in scattered data interpolation and fitting, such as the thin-plate spline method [13], inverse distance weighted method [14], Chebyshev polynomial method [15], and trigonometric polynomial method [16]. Given the anisotropy of diapycnal diffusivity and its rapid change in the magnitude in the vertical direction, the 3D thin-plate spline (3D TPS) method, 3D inverse distance weight (3D IDW) method, 3D Chebyshev polynomial fitting (CPF) method, 3D piecewise linear (3D PL) method, and 2D thin-plate spline (2D TPS) method are used to estimate the diapycnal diffusivity in the SCS. With the benefits of the 10-fold random cross-validation method [17] and the plane cross-validation method, the optimum interpolation method can be determined. The framework in this study that uses cross-validation method to select the optimum interpolation method is universal, not just applicable to the SCS. The main aim of this study is to obtain the 3D complete spatial distribution of diapycnal diffusivity in the SCS, which would provide potential applications for improving the accuracy of numerical simulation.
This paper is organized as follows. A brief description of the data is presented in Section 2. The interpolation methods are given in Section 3. Result analysis and conclusions are presented in Section 4 and Section 5, respectively.

2. Data

Yang et al. [4] estimated the diapycnal diffusivity in the SCS based on in situ observations from 2005 to 2012 and the GHP parameterization. Detailed information of the measuring data can be referred to Yang et al. [4].
The 3D structure of diapycnal diffusivity in the SCS given by Yang et al. [4] is shown in Figure 1. The horizontal distribution of diapycnal diffusivity ranges from 110° E to 121° E and 11° N to 21° N with a resolution of 0.5° × 0.5°, and the vertical distribution ranges from 200 m to 5000 m with a resolution of 200 m. Based on these data (a total of 2353 data), a 3D spatial interpolation analysis is performed in this paper to reconstruct the 3D spatial distribution of diapycnal diffusivity in the SCS.

3. Methodology

Several 3D spatial interpolation methods are used to estimate the diapycnal diffusivity at unknown points and the 2D TPS interpolation method is used as a comparison. The following contents are the methods that are commonly used for interpolating scattered point data.

3.1. The 3D TPS Interpolation

Radial basis functions have the advantage of being unconstrained by dimensions, and they are statistically unbiased estimates of minimum variance [18]. In this study, the radial basis function method is applied due to its good accuracy in the interpolation of scattered points [18,19]. The TPS function is a commonly used radial basis function, which guarantees the minimum bending energy of the deformation field [13]. The TPS interpolation has good numerical stability and convergence, and it is more suitable for global deformation [20]. Therefore, the TPS function is chosen as the basis function in this study. Assuming that the diffusivity is an isotropic distribution, there are n data points with measured diapycnal diffusivity K ( x i , y i , z i ) , I = 1, 2, …, n. Then, the estimated diapycnal diffusivity at other points K ( x i , y i , z i ) is calculated by 3D TPS interpolation as follow:
K ( x j , y j , z j ) = i = 1 n c i U ( r i j ) + λ 1 + λ 2 x i + λ 3 y i + λ 4 z i ,
U ( r i j ) = r i j 2 ln ( r i j ) ,
where r i j represents Euclidean distance between ( x i , y i , z i ) and ( x j , y j , z j ) , the parameters c i , λ 1 , λ 2 , λ 3 , λ 4 can be solved by
[ c λ ] = [ A     P P T 0 ] 1 [ K 0 ] ,
where
c = [ c 1 c 2 c n ] , λ = [ λ 1 λ 2 λ 3 λ 4 ] , P = [ 1     x 1     y 1     z 1 1     x 2     y 2     z 2                   1     x n     y n     z n ] , K = [ K 1 K 2 K n ] ,
A = [ 0 U ( r 12 ) U ( r 13 ) U ( r 1 n ) U ( r 21 ) 0        U ( r 23 ) U ( r 2 n )                                             U ( r n 1 ) U ( r n 2 )     0 ] .
However, the actual diapycnal diffusivity in the SCS is anisotropic, and its variability in the vertical direction is greater than that in the horizontal direction. In this study, the anisotropic of data, which has a crucial influence on the interpolation results, is taken into account in interpolation process.
As shown in Table 1, the mean absolute change rate (e.g., in vertical direction, calculate the absolute rate of change for the diapycnal diffusivity values of any two adjacent points in vertical direction, and then average all absolute change rate) of the diapycnal diffusivity values in the vertical direction is two orders of magnitude higher than that in the horizontal direction. In addition, turbulent mixing in the SCS is mainly triggered by internal waves breaking, and the horizontal feature scale of the typical movement is also two orders of magnitude larger than the vertical feature scale [2,9]. In other words, in the interpolation, diapycnal diffusivity in the horizontal direction should have a greater weight than that in the vertical direction [21]. Therefore, the weight in the vertical direction is set to 1/100 of that in the horizontal direction to achieve the anisotropy characteristics of diapycnal diffusivity, and Equation (1) is rewritten as:
K ( x j , y j , z j ) = i = 1 n c i j U ( r i j ) + λ 1 + λ 2 x i + λ 3 y i + λ 4 z i , z j = 100 z j , z i = 100 z i .
In the subsequent 3D interpolation methods used in this study, the weight in the vertical direction is always set to 1/100 of that in the horizontal direction.

3.2. The 3D IDW Interpolation

The 3D IDW interpolation method is based on the hypothesis that each measurement point has a local influence, which is inversely proportional to the distance between the point to be interpolated and the measurement point. All distances in this study are calculated in Cartesian coordinate system. The diapycnal diffusivity K ( x j , y j , z j ) estimated by the 3D IDW interpolation method is performed as follows:
K ( x j , y j , z j ) = i = 1 n K i ( x i , y i , z i ) d i j α i = 1 n 1 d i j α ,
where d i j = ( x j x i ) 2 + ( y j y i ) 2 + ( z j z i ) 2 , α = 2 . The distances in this article are all calculated under the Cartesian coordinate system, and the spherical coordinate system is not used.

3.3. The 3D CPF

In this paper, the diapycnal diffusivity estimated by the 3D CPF method is performed as follows:
{ K ~ C P F = k = 0 K max s = 0 S max l = 0 L max β k , s , l T k ( x ) P s ( y ) M l ( z ) x = 2 x x max x min x max x min , y = 2 y y max y min y max y min , z = 2 z z max z min z max z min T k ( x ) = cos ( k arccos ( x ) ) , P s ( y ) = cos ( s . arccos ( y ) ) , M l ( z ) = cos ( l arccos ( z ) ) ,
where K max , S max , and L max represent the highest order of polynomials in zonal, meridional and vertical direction, respectively, which are limited to 10 in this study. β k , s , l are coefficients that can be solved by the least squares method.

3.4. The 2D TPS Interpolation

The 3D diapycnal diffusivity data can be decomposed into multiple 2D data at different water depths and then 2D interpolation schemes can be adopted. To make a comparison with the 3D TPS method, the 2D TPS method is used to interpolate the diapycnal diffusivity at each depth. The 2D TPS method is introduced as follows: assuming that there are n data points with measured diapycnal diffusivity K ( x i , y i ) at a certain depth, diapycnal diffusivity at other points K ( x j , y j ) can be estimated in Section 3.1, and the only difference is that the z dimension is deleted.

3.5. Cross-Validation

In this study, the 10-fold random cross-validation [17] is employed to assess the interpolation results. The measured diapycnal diffusivity dataset is randomly divided into ten subsets, nine of which are selected as set A and one as set B. The data in set A are used to interpolate, while the accuracy of estimation is evaluated with the data in set B. The random cross-validation is repeated 50 times. By calculating the averaged error of interpolated results in the 50 experiments, the effectiveness and accuracy of the aforementioned interpolation methods can be assessed.
In order to verify the effectiveness of interpolation in where an entire section is removed, as a supplement to the 10-fold random cross validation scheme, another cross-validation scheme (plan cross-validation) is also implemented in this study. In this scheme, the measured diapycnal diffusivity on a whole certain horizontal plan is retained to evaluate the accuracy of each interpolation method, meanwhile the remaining data is used to interpolate. In order to ensure that the estimate is interpolated and there are enough measure points in the evaluated, water depth planes ranging from 400 m to 3800 m are selected for evaluation and the plane cross-validation is repeated 18 times.

3.6. Statistical Analysis

Two statistical parameters, namely the mean absolute error (MAE) and the mean relative error (MRE), are used to quantify the interpolation performance. The MAE and MRE are calculated as follows:
M A E = i = 1 m | K ~ i K i m | ,
M R E = 1 m i = 1 m | K i ~ K i K i | ,
where K i and K ~ i are measured and estimated diapycnal diffusivity, respectively, and m is the number of validation data in set B.

4. Result Analysis and Discussion

Given the rapid change of diapycnal diffusivity in the magnitude in the vertical direction (the measured diapycnal diffusivity varying between 10−5 m2/s and 10−1 m2/s), the interpolation of log10 (K) may produce better results than the interpolation of K. In addition, the prediction error could be used to determine whether the data should be transformed [22]. The 3D TPS method is taken as an example to compare the interpolation results of K and log10 (K). Equations (1) and (3) are rewritten as (11) and (12) when the interpolation is carried out on log10 (K).
log 10 ( K ( x j , y j , z j ) ) = i = 1 n c i j U ( r i j ) + λ 1 + λ 2 x i + λ 3 y i + λ 4 z i ,
[ c λ ] = [ A     P P T 0 ] 1 [ log 10 ( K )      0 ] .
The 10-fold random cross-validation results of the 3D TPS method are shown in Figure 2 which clearly shows that the interpolation errors for the logarithm of K are much smaller than those for K, suggesting the reliability of interpolating the logarithm of K. Noting that the interpolated log10 (K) values are converted back to K values before calculating the interpolation errors. Therefore, the logarithmic pretreatment of diapycnal diffusivity is taken in the following experiments.
All the combinations of K max = 0 , 1 , , 10 , S max = 0 , 1 , , 10 , and L max = 0 , 1 , , 10 are tested when we fit the diapycnal diffusivity by using the 3D CPF method. The M A E ¯ ( M A E ¯ = ( i = 1 N M A E ( i ) / N , N is the total number of random cross-validation experiments) and the M R E ¯ ( M R E ¯ = ( i = 1 N M R E ( i ) ) / N , N is the total number of random cross-validation experiments) for each combination are calculated during the 10-fold random cross-validation process [16]. In order to facilitate the comparison with the 3D TPS method, M A E ¯ and M R E ¯ with the minimum are regarded as fitting errors and shown in Figure 3. The degrees corresponding to the minimum values of M A E ¯ and M R E ¯ in the 10-fold random cross-validation process are given in Table 2.
As shown in Figure 3, the errors of the 3D CPF method ( M A E ¯ = 0.0030   m 2 / s , M R E ¯ = 85.06 % ) are unreliable and unstable during the random cross-validation process, which suggests that the 3D CPF method may not be a good approach in this study. In order to expand the measured diapycnal diffusivity into a spatial complete field, other 3D interpolation methods were applied. Among them, the 3D piecewise linear (3D PL) interpolation method was implemented by MATLAB griddata function (linear interpolation as an interpolation method) and used as a comparison. The 10-fold random cross-validation results of these 3D interpolation methods are shown in Figure 4. It is clearly shown that the 3D TPS method performs best in the 10-fold random cross-validation with M A E ¯ = 0.0013   m 2 / s and M R E ¯ = 32.65 % , followed by the 3D PL method with M A E ¯ = 0.0016   m 2 / s and M R E ¯ = 45.67 % .
The M A E ¯ s and M R E ¯ s of the 2D TPS method at different depths are illustrated in Figure 5. The MAE and MRE of all points to be estimated are 0.0043 m2/s and 172.51%, respectively, which are much larger than the errors of the 3D TPS method. In addition, it is obvious that the 2D interpolation method cannot accurately describe the true distribution of geographic phenomena in the 3D space, and the interpolation effect is easily affected by the number of measured points on the plane. Therefore, the 2D TPS method is not suitable for 3D reconstruction of the diapycnal diffusivity field.
Although the background diffusivity is 10−5 m2/s or even 10−4 m2/s, considering that the average value of the measured diapycnal diffusivity is 0.0052 m2/s and the M R E ¯ of the 3D TPS method is 32.65%, the M A E ¯ of the 3D TPS method is reasonable. The first random cross-validation experiment is taken as an example: Figure 6 shows that the diapycnal diffusivity obtained by the 3D TPS method has a good agreement with the measured diapycnal diffusivity in set B and the relative errors of most points to be inspected fluctuate in [−50%, 50%].
The statistical results of the overestimation and underestimation of the 3D interpolation methods are shown in Table 3. According to the statistical results in Table 3, the interpolated results are more likely to be significantly overestimated (the average proportion of relative error greater than 50%) for this diapycnal diffusivity field. Among the 3D TPS, PL, and IDW methods, the 3D TPS method yields the smallest proportions for both overestimation and underestimation, although the low error proportion of the 3D TPS method is only a slight improvement on that of the 3D PL method.
Results of plane cross-validation experiments are shown in Figure 7. All cross-validation processes suggest that the 3D TPS method produces better results (the MAE and MRE of all points to be estimated are 0.0011 m2/s and 27.70%, respectively) than the 3D PL and IDW methods. The MAE and MRE of all points with 3D PL method are 0.0015 m2/s and 39.08% respectively. Considering that the 3D PL method requires the least amount of calculation among the 3D interpolation methods, the error level of the 3D PL method is encouraging.
It is interesting to find that the errors of the 3D IDW are much larger than those of the 3D PL and TPS methods. It can be seen from Figure 8 that the diapycnal diffusivity is obviously overestimated in the weak mixing area, while it is underestimated in the high mixing area. This phenomenon can be seen more clearly from random cross-validation (Figure 9). The main reason is that turbulent mixing generally increases with depth in the SCS, which leads to the fact that the measured points in the deep sea have a greater impact on the interpolation points than points closer to the interpolation points (inconsistent with the assumption of the 3D IDW method). Moreover, the interpolated values with the 3D IDW method within the data set are bounded by min ( K i ) < K j < max ( K i ) , where K i and K j are measured and estimated diapycnal diffusivity, respectively. In other words, the 3D IDW method is essentially composed of smoothing procedures [14], which cause the diapycnal diffusivity to be easily underestimated in the high mixing area and overestimated in the weak mixing area.
The 3D TPS and PL methods are selected to reconstruct spatial fields of diapycnal diffusivity in the plane cross-validation. Figure 10 shows the depth-averaged diffusivity reconstructed by interpolation in the upper (400–600 m), intermediate (600–1600 m), and deep layers (>1600 m). Since the 3D TPS method ensures that the first derivative of the space field is continuous [23], the maps obtained by the 3D TPS method are smoother than those obtained by the 3D PL method, and the gradient changes are smaller. In addition, both cross-validation methods show that the errors of the 3D TPS method are smaller than those of the 3D PL method. Moreover, since the derivative of the diapycnal diffusivity appears in the equations of motions, the application effect of the 3D TPS method may be better than that of the 3D PL method. In a word, results calculated by the 3D TPS are a slight improvement over those calculated by the 3D PL method.
Combining all the above analysis, the best interpolation method of diapycnal diffusivity in the SCS is the 3D TPS method, and both cross-validation methods provide evidence for the rationality of this method. The interpolated results of the 3D TPS method with horizontal resolution of 0.25° × 0.25° and vertical resolution of 100 m are show in Figure 11. The interpolated field of diapycnal diffusivity can be obtained with high spatial resolution, but it does not include the information about the diapycnal diffusivity changing with time. Therefore, the cross-validation results can only prove that the interpolation accuracy is satisfactory. The real assessment will be whether the interpolation results of the 3D TPS method improve the model’s representation of reality.

5. Conclusions

In order to expand the 3D scatter data of diapycnal diffusivity into a spatial complete field, the 3D TPS method which accounts for anisotropy (vertical scales smaller than the horizontal ones) in the data is carried out to estimate the diapycnal diffusivity in the SCS. Given the rapid change of diapycnal diffusivity in the magnitude, better results are obtained by logarithmic pretreatment of diapycnal diffusivity. The selection of spatial interpolation method for estimating the diapycnal diffusivity is also discussed in this paper. Results of 10-fold random cross-validation experiments illustrate that the 3D TPS method has the smallest errors ( M A E ¯ = 0.0013   m 2 / s , M R E ¯ = 32.65 % ), followed by the 3D PL method ( M A E ¯ = 0.0016   m 2 / s , M R E ¯ = 45.67 % ). Results of plane cross-validation experiments also indicate that the 3D TPS method can reconstruct the diapycnal diffusivity field relatively accurately. Moreover, the interpolated results may provide a reference for calculating vertical fluxes across the basin, and we will study what the error induced in total flux is by having transects in the subsequent work.
Both cross-validation methods show that the diapycnal diffusivity obtained by the 3D interpolation has a good agreement with the measured diapycnal diffusivity. The diapycnal diffusivity field reconstructed by the 3D TPS method is differentiable, with lower errors than those with the 3D PL method. However, note that for the diapycnal diffusivity field, the 3D TPS method is a slight improvement over the 3D PL method. The 3D CPF method is considered to be inappropriate in estimating the diapycnal diffusivity due to its excessive and unstable errors. Meanwhile, the diapycnal diffusivity map reconstructed by the 3D IDW method is extremely inconsistent with the measured data, and this shortcoming is related to the nature of the IDW function itself.
The 3D TPS method applied in this study, which accounts for anisotropy in the data, may provide potential applications for improving the mixing configuration in the numerical model used for simulating a realistic ocean state in the SCS. Moreover, the 3D interpolation methods used in this study might be also applicable to data analysis in any basin (adjust the ratio of horizontal weight and vertical weight appropriately), but their interpolation accuracy needs to be evaluated. A forthcoming validation is planned using the interpolation results of the 3D TPS method to test an ocean model of the SCS against standard model configurations for the diapycnal diffusivity.

Author Contributions

Conceptualization, X.L.; Formal analysis, S.L.; Funding acquisition, X.L.; Investigation, J.G.; Methodology, J.G.; Software, J.G.; Supervision, X.L.; Writing—original draft, J.G.; Writing—review & editing, Y.N. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the National Natural Science Foundation of China (Grant No. 41876003 and U1806214) and the National Key Research and Development Plan of China (Grant No. 2017YFA0604101).

Acknowledgments

We are grateful to Jiwei Tian for providing the measured data of diapycnal diffusivity. We would like to thank Ruichen Cao and Anzhou Cao for their friendly help. We deeply thank the reviewers and editor for their constructive suggestions on the earlier version of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gregg, M.C.; Sanford, T.B.; Winkel, D.P. Reduced mixing from the breaking of internal waves in equatorial waters. Nature 2003, 422, 513–515. [Google Scholar] [CrossRef] [PubMed]
  2. Garrett, C. Mixing with latitude. Nature 2003, 422, 477–478. [Google Scholar] [CrossRef] [PubMed]
  3. Tian, J.; Yang, Q.; Zhao, W. Enhanced diapycal mixing in the South China Sea. J. Phys. Oceanogr. 2009, 39, 3191–3203. [Google Scholar] [CrossRef]
  4. Yang, Q.; Zhao, W.; Liang, X.; Tian, J. Three-dimensional distribution of turbulent mixing in the South China Sea. J. Phys. Oceanogr. 2016, 46, 769–788. [Google Scholar] [CrossRef] [Green Version]
  5. Osborn, T.R. Estimates of the Local Rate of Vertical Diffusion from Dissipation Measurements. J. Phys. Oceanogr. 1980, 10, 83–89. [Google Scholar] [CrossRef] [Green Version]
  6. Alford, M.H.; MacKinnon, J.A.; Nash, J.D.; Simmons, H.; Pickering, A.; Klymak, J.M.; Pinkel, R.; Sun, O.; Rainville, L.; Musgrave, R.; et al. Energy flux and dissipation in luzon strait: Two tales of two ridges. J. Phys. Oceanogr. 2011, 41, 2211–2222. [Google Scholar]
  7. Buijsman, M.C.; Legg, S.; Klymak, J. Double-Ridge Internal tide interference and its effect on dissipation in Luzon strait. J. Phys. Oceanogr. 2012, 42, 1337–1356. [Google Scholar] [CrossRef]
  8. Lozovatsky, I.; Liu, Z.; Fernando, H.J.; Hu, J.; Wei, H. The TKE dissipation rate in the northern South China Sea. Ocean Dyn. 2013, 63, 1189–1201. [Google Scholar] [CrossRef]
  9. Alford, M.H.; Coauthors. The formation and fate of internal waves in the South China Sea. Nature 2015, 521, 65–69. [Google Scholar] [CrossRef]
  10. Chang, M.H.; Lien, R.C.; Tang, T.Y.; Dasaro, E.A.; Yang, Y.J. Energy flux of nonlinear internal waves in northern South China Sea. Geophys. Res. Lett. 2006, 33, 155–170. [Google Scholar] [CrossRef] [Green Version]
  11. Shang, X.D.; Liang, C.R.; Chen, G.Y. Spatial distribution of turbulent mixing in the upper ocean of the South China Sea. Ocean Sci. 2017, 13, 503–519. [Google Scholar] [CrossRef] [Green Version]
  12. Ezer, T. On the seasonal mixed layer simulated by a basin-scale ocean model and the Mellor-Yamada turbulence scheme. J. Geophys. Res. 2000, 105, 16843–16855. [Google Scholar] [CrossRef]
  13. Bookstein, F.L. Principal Warps: Thin-Plate Splines and the Decomposition of Deformations. IEEE Trans. Pattern Anal. Mach. Intell. 1989, 11, 567–585. [Google Scholar] [CrossRef] [Green Version]
  14. Lam, N.S.N. Spatial interpolation methods: A review. Am. Cartogr. 1983, 10, 129–150. [Google Scholar] [CrossRef]
  15. Li, B.; Liu, Y.; Wang, X.; Fu, Q.; Lv, X. Application of the orthogonal polynomial fitting method in estimating PM2.5 concentrations in central and southern regions of China. Int. J. Environ. Res. Public Health 2019, 16, 1418. [Google Scholar]
  16. Nie, Y.; Wang, Y.; Lv, X. Acquiring the Arctic-Scale Spatial Distribution of Snow Depth Based on AMSR-E Snow Depth Product. J. Atmos. Ocean. Technol. 2019, 36, 1957–1965. [Google Scholar] [CrossRef]
  17. Geisser, S. The Predictive Sample Reuse Method with Applications. J. Am. Stat. Assoc. 1975, 70, 320–328. [Google Scholar] [CrossRef]
  18. Franke, R. Scattered Data Interpolation: Tests of Some Method. Math. Comput. 1982, 38, 181–200. [Google Scholar]
  19. Wu, Z.; Schaback, R. Local error estimates for radial basis function interpolation of scattered data. IMA J. Numer. Anal. 1993, 13, 13–27. [Google Scholar] [CrossRef]
  20. Bartoli, A.; Perriollat, M.; Chambon, S. Generalized thin-plate spline warps. Int. J. Comput. Vis. 2010, 88, 85–110. [Google Scholar] [CrossRef]
  21. Merwade, V.M.; Maidment, D.R.; Maidment, J.A. Anisotropic considerations while interpolating river channel bathymetry. J. Hydrol. 2006, 331, 731–741. [Google Scholar] [CrossRef]
  22. Li, J.; Heap, A.D. Spatial interpolation methods applied in the environmental sciences: A review. Environ. Model. Softw. 2014, 53, 173–189. [Google Scholar] [CrossRef]
  23. Franke, R. Smooth interpolation of scattered data by local thin plate spline. Comput. Math. Appl. 1982, 8, 273–281. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The 3D structure of diapycnal diffusivity (m2/s; logarithmic scale); the meridional sections (upper panel) and the zonal sections (lower panel).
Figure 1. The 3D structure of diapycnal diffusivity (m2/s; logarithmic scale); the meridional sections (upper panel) and the zonal sections (lower panel).
Jmse 08 00832 g001
Figure 2. The 10-fold random cross-validation results of the interpolation of K and log10 (K).
Figure 2. The 10-fold random cross-validation results of the interpolation of K and log10 (K).
Jmse 08 00832 g002
Figure 3. MAEs and MREs of the 3D CPF method for 50 groups of experiments.
Figure 3. MAEs and MREs of the 3D CPF method for 50 groups of experiments.
Jmse 08 00832 g003
Figure 4. MAEs (a) and MREs (b) of the 3D interpolation methods for 50 groups of experiments.
Figure 4. MAEs (a) and MREs (b) of the 3D interpolation methods for 50 groups of experiments.
Jmse 08 00832 g004
Figure 5. M A E ¯ s and M R E ¯ s at different water depths of the 2D TPS method.
Figure 5. M A E ¯ s and M R E ¯ s at different water depths of the 2D TPS method.
Jmse 08 00832 g005
Figure 6. Measured diapycnal diffusivity (m2/s; logarithmic scale) in set B versus diapycnal diffusivity obtained by the 3D TPS interpolation method for the first random cross-validation experiment. The red line and blue lines are 1:1 line and ±50% relative error lines, respectively.
Figure 6. Measured diapycnal diffusivity (m2/s; logarithmic scale) in set B versus diapycnal diffusivity obtained by the 3D TPS interpolation method for the first random cross-validation experiment. The red line and blue lines are 1:1 line and ±50% relative error lines, respectively.
Jmse 08 00832 g006
Figure 7. MAEs (a) and MREs (b) of the 3D interpolation methods at different levels in the plane cross-validation experiments.
Figure 7. MAEs (a) and MREs (b) of the 3D interpolation methods at different levels in the plane cross-validation experiments.
Jmse 08 00832 g007
Figure 8. (a) Map of the measured diapycnal diffusivity (m2/s; logarithmic scale) at a depth of 600 m. (b) Map of the diapycnal diffusivity reconstructed by the 3D IDW method at a depth of 600 m.
Figure 8. (a) Map of the measured diapycnal diffusivity (m2/s; logarithmic scale) at a depth of 600 m. (b) Map of the diapycnal diffusivity reconstructed by the 3D IDW method at a depth of 600 m.
Jmse 08 00832 g008
Figure 9. Measured diapycnal diffusivity (m2/s; logarithmic scale) in set B versus diapycnal diffusivity obtained by the 3D IDW interpolation method (the first random cross-validation experiment is taken as an example). The red line is 1:1 line.
Figure 9. Measured diapycnal diffusivity (m2/s; logarithmic scale) in set B versus diapycnal diffusivity obtained by the 3D IDW interpolation method (the first random cross-validation experiment is taken as an example). The red line is 1:1 line.
Jmse 08 00832 g009
Figure 10. (ac) Maps of depth-averaged measured diapycnal diffusivity (m2/s; logarithmic scale). (df) Maps of depth-averaged diapycnal diffusivity reconstructed by the 3D TPS method. (gi) Maps of depth-averaged diapycnal diffusivity reconstructed by the 3D PL method. The order from top to bottom are the upper layer (400–600 m), intermediate layer (600–1600 m), and deep layer (>1600 m) respectively.
Figure 10. (ac) Maps of depth-averaged measured diapycnal diffusivity (m2/s; logarithmic scale). (df) Maps of depth-averaged diapycnal diffusivity reconstructed by the 3D TPS method. (gi) Maps of depth-averaged diapycnal diffusivity reconstructed by the 3D PL method. The order from top to bottom are the upper layer (400–600 m), intermediate layer (600–1600 m), and deep layer (>1600 m) respectively.
Jmse 08 00832 g010
Figure 11. The meridional sections of the 3D structure of diapycnal diffusivity (m2/s; logarithmic scale) reconstructed by the 3D TPS method.
Figure 11. The meridional sections of the 3D structure of diapycnal diffusivity (m2/s; logarithmic scale) reconstructed by the 3D TPS method.
Jmse 08 00832 g011
Table 1. The mean absolute rate of change of the diffusivity values in meridional, zonal and vertical directions.
Table 1. The mean absolute rate of change of the diffusivity values in meridional, zonal and vertical directions.
DirectionsThe Mean Absolute Rate of Change (m2s−1km−1)
Zonal direction 6.5 × 10 5
Meridional direction 6.7 × 10 5
Vertical direction 9.1 × 10 3
Table 2. The degrees corresponding to M A E ¯ and M R E ¯ with the minimum of the 3D CPF method during the 10-fold cross-validation process.
Table 2. The degrees corresponding to M A E ¯ and M R E ¯ with the minimum of the 3D CPF method during the 10-fold cross-validation process.
ErrorsKmaxSmaxLmax
M A E ¯ 493
M R E ¯ 593
Table 3. The average proportion of relative errors of the 3D interpolation methods for 50 groups of experiments.
Table 3. The average proportion of relative errors of the 3D interpolation methods for 50 groups of experiments.
3D Interpolation MethodsThe Average Proportion of Relative Error Greater than 50%The Average Proportion of Relative Error Less than −50%
3D TPS13.20%4.29%
3D PL19.87%6.34%
3D IDW36.36%26.95%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Guo, J.; Nie, Y.; Li, S.; Lv, X. Application of Three-Dimensional Interpolation in Estimating Diapycnal Diffusivity in the South China Sea. J. Mar. Sci. Eng. 2020, 8, 832. https://doi.org/10.3390/jmse8110832

AMA Style

Guo J, Nie Y, Li S, Lv X. Application of Three-Dimensional Interpolation in Estimating Diapycnal Diffusivity in the South China Sea. Journal of Marine Science and Engineering. 2020; 8(11):832. https://doi.org/10.3390/jmse8110832

Chicago/Turabian Style

Guo, Junting, Yafei Nie, Shuang Li, and Xianqing Lv. 2020. "Application of Three-Dimensional Interpolation in Estimating Diapycnal Diffusivity in the South China Sea" Journal of Marine Science and Engineering 8, no. 11: 832. https://doi.org/10.3390/jmse8110832

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