Thermohaline Dynamics in the Northern Continental Slope of the South China Sea: A Case Study in the Qiongdongnan Slope

<jats:p>Understanding the marine hydro-thermohaline environment is essential for terrestrial meteorology and the coastal ecosystem. Here, we provide insight into the hydro-thermohaline environment at the Qiongdongnan continental slope of the northern South China Sea and the mechanism controlling it, with focus on its short-term characteristics. We employ a well-validated three-dimensional unstructured-grid-based Finite Volume Coastal Ocean Model (FVCOM) to analyze the spatial-temporal behavior of its hydro-thermohaline structures and to quantify the transport fluxes over a full tidal period. The analysis reveals a two-layer flow structure with directionally oppositely moving layers in the along-isobaths direction. Furthermore, transport patterns undergo periodic changes. During the spring tide, the downslope (along-isobaths) transport of water/heat/salt is approximately 119%/70%/120% higher (62%/62%/62% lower) than during the neap tide. From analyzing the different terms in the thermohaline balance equation, we find that the main dynamic factors controlling heat transport over a tidal period are the gravitational convention and the mean flow, while the salt transport is only dominated by the mean flow. The data of the short-term thermohaline evolution of the QDNS provided in this study may be of use for future studies of the northern SCS, including its marine ecology and marine fisheries.</jats:p>

section_content: The South China Sea (SCS) is located in the northwest of the Pacific Ocean, south of the Chinese landmass (Figure 1a ). Covering an area of 3. 5 million km 2, with an average depth of over 2000 m and a maximum depth of 5560 m, the SCS is the largest and deepest sea in that region [1]. The SCS is widely recognized for a number of thermohaline dynamic phenomena, including upwelling [2, 3], thermoclines [4, 5], and mesoscale eddies [6, 7]. A better physical understanding of them can help improve the study of marine fisheries development, marine pollutant countermeasures and marine ecosystem protection. For example, the continental slope in the SCS serves as spawning, nursery, and feeding grounds for various economically valuable fisheries [8], while temperature and salinity are important factors affecting the fish species and productivity [9] [10] [11]. Thus, the study of the thermohaline environment in the continental slope area of the SCS contributes to the ecological environment on which fish reproduction depends. In addition, the continental slope in the northern SCS is the largest reservoir of microplastics [12, 13], indicating the importance of regional hydrodynamic study in pollution prevention and ecological protection. Moreover, ocean dynamics research in the SCS plays a crucial role in understanding the regional weather and climate due to its connection with extreme weather disasters such as El Niño Southern Oscillation [14]. The archipelago and the bumpy continental slope/shelf around the deep basin lead to a complex topography in the northern SCS [15], which plays a vital role in shaping the northern SCS thermohaline structure. The continental shelf on the southeast side of Hainan Island, a broad and shallow terrace less than 200 m deep, extends 80-90 km seaward (Figure 1b ). Its isobaths are parallel to the southern coastline of Hainan Island (Figure 1b ). The Qiongdongnan continental slope (QDNS, indicated by the red arrow in Figure 1a and the lilac rectangle in Figure 1b ), with a typical slope angle of 4. 5-10 • and a correspondingly sharp increase of water depth from 200 m to 1000-1800 m (Figure 1b, c ), connects the Hainan continental shelf to the Qiongdongnan basin [16]. Slopes as steep as this are responsible for a relatively narrow cross-shore width of the surface Ekman divergence [17], significant vertical turbulent mixing and thus uplifting of the thermocline [18], and strong interior transport between the surface and bottom boundary layers [19]. In contrast, a wide continental shelf dissipates too little energy to lift the thermocline upward. In addition to steep slopes, the funnel effect caused by nearshore topography features, such as the cape, can accelerate the current, which also increases vertical mixing [20, 21]. 

Previous studies of the hydrodynamics of the South China Sea focused on large-scale and mesoscale phenomena, such as the general ocean circulation [22, 23] and mesoscale eddies [6, 24], while the dynamic features within the continental slope region were rarely discussed. Furthermore, several studies investigated interannual and seasonal averaged hydro-thermohaline variations in the SCS [25] [26] [27]. However, short-term changes, such as tidal-cycle and daily differences, have received less attention. 

This study numerically investigates the three-dimensional spatial characteristics and short-term temporal features of the hydro-thermohaline structure in the QDNS. The goal is to quantitatively analyze the spatial-temporal variation of water, heat, and salt fluxes around the QDNS, and to elucidate the hydrodynamic mechanisms controlling thermohaline transport there.
section_name: Introduction

section_content: This study uses a computational fluid dynamics (CFD) numerical model, introduced in Section 2. 1, to simulate the hydrodynamics in the QDNS. The fluxes of water, heat, and salt transport obtained from the simulations are analyzed using the thermohaline balance equation, introduced and explained in Section 2. 2, to determine the dominant controlling factors.
section_name: Data and Methods

section_content: This study employs an unstructured grid, Finite Volume Coastal Ocean Model (FV-COM) [28], to simulate the three-dimensional, time-dependent circulation in the QDNS of the northern SCS at a relatively short temporal resolution of days to weeks. The unstructured triangular grid provides high flexibility and variability, which ensures high spatial resolution model results by giving a small-sized grid in the key study areas such as QDNS in this paper. Using a finite-volume approach, FVCOM is well-suited to maintain mass, momentum, temperature, and salinity conservation over complex topography such as QDNS, improving the calculation accuracy of the transport fluxes of water, heat and salt. The model can be run in two distinct modes using the model-split solver: an external mode which solves the vertically integrated (two-dimensional) equations of motion and an internal mode which resolves the three-dimensional structure.
section_name: Numerical Model

section_content: The model domain of about 42,315 km 2 is bounded by Hainan Island at the land side and exhibits an open boundary at the sea side, indicated by the solid black line in Figure 1b. Figure 1d depicts the unstructured triangular grid of the numerical model, which contains 24,889 nodes and 48,792 cells. The size of a horizontal grid cell is about 5000 m along the open boundary and gradually decreases to about 200 m in the QDNS and nearshore area to resolve short-term slope-induced hydrodynamics (Figure 1d ). Vertically, for the internal mode, the model employs a hybrid, terrain-following coordinate system consisting of 40 layers to match the complex bathymetry between the continental shelf area and the Qiongdongnan basin (Figure 1b ). Precisely, for sea floor locations with a depth of less than 80 m, these layers are evenly spaced, whereas for sea floor locations with a depth of more than 80 m, the near-floor and near-free-surface layers are more densely spaced (5 layers with a 2 m spacing) than those in between. 

Coastline data is obtained from the Global Self-consistent, Hierarchical, High-resolution Geography Database (GSHHG) produced by the National Geophysical Data Center (NGDC) [29]. Topography is constructed based on the NGDC's ETOPO1 Global Relief Model [30], with 1 • resolution in longitude-latitude.
section_name: Model Domain

section_content: The model run covers the time from 15 June to 1 October 2020. The first 16 days of the total simulation period are used to reach a quasi-equilibrium, and the remaining 92 days for analysis. The time step used in the external mode of the FVCOM is 1 s, which is bounded through the Courant-Friedrich Levy (CFL) stability criterion [31], while the time step of the internal mode is 5 s. Hourly output results ensure the high temporal precision needed to resolve the short-term ocean dynamics in the QDNS. 

The model is forced with initial, open boundary and free surface conditions, but coldstarted with zero tidal elevation and flow velocity. We initialized thermohaline conditions using the Global Ocean Forecasting System (GOFS 3. 1) of the HYbrid Coordinate Ocean Model (HYCOM) [32]. The database has a horizontal resolution of 1/12 • in both longitude and latitude, a vertical resolution of 41 layers, and a temporal frequency of 3 h. 

Active open boundary conditions are applied to integrate external tidal and thermohaline forces. Hourly tidal elevations along the open boundary were predicted by 13 harmonic components (K 1, O 1, P 1, Q 1, M 2, S 2, N 2, K 2, M 4, MS 4, MN 4, M f, M m ) from the TPXO 7. 2 global tide Model [33]. The external forcing of the thermohaline variables along the open boundary is also provided by HYCOM, GOFS 3. 1 [32]. 

The free surface is forced by climatological wind stress and atmospheric heat flux, calculated from the u/v-wind component (0. 2 • × 0. 2 • ) 10 m above the ground and air-sea interface thermal parameters (0. 2 • × 0. 2 • ), including latent heat flux, sensible heat flux, downward shortwave radiation flux, upward shortwave radiation flux, downward longwave radiation flux and upward longwave radiation flux at the water surface, respectively. Both wind and heat flux data are derived from the National Centers for Environmental Prediction (NCEP) Climate Forecast System Version 2 (CFSv2) [34]. The relevant information on model forcing is organized in Table 1.
section_name: Model Configurations

section_content: We calibrate the model outputs of sea level, flow velocity and direction, temperature, and salinity using the correlation coefficient (CC), the Willmott skill score (WSS), and mean relative error (MRE):

where m i represents the simulation results, and O i refers to the observation results. m(S m ) and O(S o ) are the average values (standard deviations) of the calculated and observed values, respectively. WSS and CC values close to unity mean a significant correlation between the numerical simulation and field observations [35]. Figure 2 compares the simulation results (black lines) with the reference data (red dots) for July 2020, provided by the Sanya tide gauge (109 • 30 E, 18 • 14 N) and collected by the vessel-mounted Acoustic Doppler Current Profiler (ADCP) during a hydrographic survey conducted in the summer of 2020. The survey also provides several thermohaline profiles data collected by two glider-mounted Conductivity, Temperature, Depth (CTD) sensors, which are used to calibrate the temperature and salinity field in the model (Figure 3 ). The sea level, flow, and thermohaline validate stations are indicated by the red dot, triangles and diamonds, respectively, in Figure 1b. 

dots) for July 2020, provided by the Sanya tide gauge (109°30′ E, 18°14′ N) and collected by the vessel-mounted Acoustic Doppler Current Profiler (ADCP) during a hydrographic survey conducted in the summer of 2020. The survey also provides several thermohaline profiles data collected by two glider-mounted Conductivity, Temperature, Depth (CTD) sensors, which are used to calibrate the temperature and salinity field in the model (Figure 3 ). The sea level, flow, and thermohaline validate stations are indicated by the red dot, triangles and diamonds, respectively, in Figure 1b. It can be seen that the distribution and magnitude of the simulation data are consistent with the measurements, suggesting that the model can adequately simulate the flow, sea level, and thermohaline characteristics of the study area. This is supported by CC and WSS values close to unity for the sea level and flows and small MREs for temperature and salinity (Table 2 ). Note that the flow, temperature and salt validation are based on short-term measurements, supporting that the FVCOM model can be applied to study the short-term slope dynamics in the QDNS. It can be seen that the distribution and magnitude of the simulation data are consistent with the measurements, suggesting that the model can adequately simulate the flow, sea level, and thermohaline characteristics of the study area. This is supported by CC and WSS values close to unity for the sea level and flows and small MREs for tem- perature and salinity (Table 2 ). Note that the flow, temperature and salt validation are based on short-term measurements, supporting that the FVCOM model can be applied to study the short-term slope dynamics in the QDNS.
section_name: Model Validation

section_content: The instantaneous rate of transport of salt through a unit width of a section perpendicular to the mean flow is given by Dyer [36]

where h is the depth and s is the salinity. At any depth, one can separate u and s via

where u z and s z are the observed values averaged over an hour, and u and s are their turbulent variations. u z and s z may be further separated into mean and deviation with respect to the depth averaging procedure:

where u v and s v are the depth deviations from the depth averaged values u and s. Because of tidal fluctuations, u and s will fluctuate greatly over a tidal cycle. Consequently, separating into mean and fluctuation with respect to a time-averaging procedure yields u = u + U and s = s + S, and

where U and S are the tidal cycle time fluctuations of the tidal cycle mean values u and s. 

Putting everything together yields

where the brackets denote the depth average of the products u v s v and u s, and h is the mean depth. The first term Q 1 represents the contribution to the total salt flux by the Eulerian mean flow. The second term Q 2, generally known as tidal pumping, occurs when the tide does not act as a purely standing wave. The third term Q 3 arises from the velocity and salinity variation with depth, the shear effect, and will be finite if the two are correlated. Q 3 can also be termed a gravitational mode and includes the effect of vertical diffusion. The fourth term Q 4 represents the salt flux owing to the short period of turbulence. We neglect this term because it is usually much smaller than the other terms [36]. 

The heat balance equation is similar to the salt balance equation introduced above, except that the parameter s is replaced by ρC p t, where t is the temperature, ρ is the water density, and C p = 4168 J/kg • • C is the specific heat capacity.
section_name: Thermohaline Balance

section_name: Results and Discussions

section_content: In order to study the short-term dynamics in the QDNS, we analyzed the flow field in the research area at two characteristic moments (i. e., peak flood and flood slack), with focus on the similarities and differences of the flow field between the spring tide and neap tide periods. 

At peak flood, the flow field exhibits the following characteristics: (i) a strong influx of surface water at the entire western open boundary, (ii) a strong overall influx of water from both open boundaries in the continental shelf area due to strong flows in the middle (and lower) water layers, and (iii) a strong surface flow along the isobaths in the basin area directed northeastward (Figure 4 ). These features are much more pronounced during the spring tide (Figure 4a-c ) than during the neap tide (Figure 4g-i ). At flood slack, the flow in the basin area remains, but not at the same strength, whereas the flow in the continental shelf area no longer exhibits a clear characteristic. 

A cyclonic eddy distinctly exists during the spring tide within the areas between 110 • N to 110. 8 • E, and 16. 8 • N to 17. 6 • N, which is indicated by the pink arrow in Figure 4b,e,f. The eddy moves northeastward from the spring tide (shown by the pink arrow in Figure 4b,e) to the neap tide (shown by the orange arrow in Figure 4h,k), with a substantially reduced maximum velocity (from 0. 7 m/s to 0. 3 m/s). Considering the complex topography of the QDNS, six representative sections are chosen for presenting the vertical flow structure (Figure 1c ). Section 1 in Figure 1c is approximately directed along the isobaths of the QDNS, while Sections 2-6 are perpendicular to Section 1 in the downslope direction. Contours and arrows in Figure 5 depict the water flow velocities across and along these sections, respectively, at two characteristic instances. It can be seen that the water flows upslope east of 111. 1°E and downslope west of 111. 1°E during the spring tide (Figure 5a1 ). Water flow across the isobath mainly occurs at depths of less than 200 m. The maximum value of the upslope flow velocity increases from 0. 2 m/s (Figure 5a1 ) to 0. 4 m/s (Figure 5a2 ), indicating a stronger cross-isobaths water exchange during the neap tide than the spring tide. It can also be inferred from the crosssectional flow velocity in Sections 2-6 that, along the isobaths, the surface water flows northeastward near the QDNS, while the bottom water flows southwestward in the basin area (Figure 5b-f ). Considering the complex topography of the QDNS, six representative sections are chosen for presenting the vertical flow structure (Figure 1c ). Section 1 in Figure 1c is approximately directed along the isobaths of the QDNS, while Sections 2-6 are perpendicular to Section 1 in the downslope direction. Contours and arrows in Figure 5 depict the water flow velocities across and along these sections, respectively, at two characteristic instances. It can be seen that the water flows upslope east of 111. 1 • E and downslope west of 111. 1 • E during the spring tide (Figure 5a1 ). Water flow across the isobath mainly occurs at depths of less than 200 m. The maximum value of the upslope flow velocity increases from 0. 2 m/s (Figure 5a1 ) to 0. 4 m/s (Figure 5a2 ), indicating a stronger cross-isobaths water exchange during the neap tide than the spring tide. It can also be inferred from the cross-sectional flow velocity in Sections 2-6 that, along the isobaths, the surface water flows northeastward near the QDNS, while the bottom water flows southwestward in the basin area (Figure 5b-f ).
section_name: Hydrodynamic Structure

section_content: Qiongdong upwelling pumps the cold salt water from the seafloor upwards [2], influencing the distributions of the sea surface temperature (SST) and sea surface salinity (SSS) within the study area (Figure 6a, c ). Compared with the basin area, the SST near the coast of Hainan Island, indicated by the green arrow, is lower by 5 °C while salinity is higher by 0. 3 psu. This result is consistent with reports by Jing et al. [1] and Bai et al. [37]. Outside the coastal upwelling area, the SST in the model domain is approximately 30 °C, except at peak flood of the neap tide when temperatures rise to 32 °C because of peaking air-sea heat flux (Figure 6i ). This implies that SST is temporally associated with air-sea heat flux.
section_name: Thermohaline Structure

section_content: Qiongdong upwelling pumps the cold salt water from the seafloor upwards [2], influencing the distributions of the sea surface temperature (SST) and sea surface salinity (SSS) within the study area (Figure 6a, c ). Compared with the basin area, the SST near the coast of Hainan Island, indicated by the green arrow, is lower by 5 • C while salinity is higher by 0. 3 psu. This result is consistent with reports by Jing et al. [1] and Bai et al. [37]. Outside the coastal upwelling area, the SST in the model domain is approximately 30 • C, except at peak flood of the neap tide when temperatures rise to 32 • C because of peaking air-sea heat flux (Figure 6i ). This implies that SST is temporally associated with air-sea heat flux. The SSS distributions in Figure 6c, g, k, o are strongly spatially correlated with the surface flow field in Figure 4a, d, g, j, suggesting that advection contributes to salt transport. As the lilac arrows in Figures 4a and 6c show, the flow comes from the southwest open boundary and the Qiongdong upwelling area brings fresh and salt water, respectively. Thus, the SSS on the left side of Section 1 in Figure 1c differs from that on the right side (Figure 7a2 ). Vertically, the halocline appears at 100~250 m. In the along-isobaths direction, the halocline east of 110. 8° E is thicker and saltier than that of the west (Figure 7a2 ). In contrast, Sections 2-6 presented a weaker spatial distinction of the halocline in the along-section direction. The SSS distributions in Figure 6c, g, k, o are strongly spatially correlated with the surface flow field in Figure 4a, d, g, j, suggesting that advection contributes to salt transport. As the lilac arrows in Figures 4a and 6c show, the flow comes from the southwest open boundary and the Qiongdong upwelling area brings fresh and salt water, respectively. Thus, the SSS on the left side of Section 1 in Figure 1c differs from that on the right side (Figure 7a2 ). Vertically, the halocline appears at 100~250 m. In the along-isobaths direction, the halocline east of 110. 8 • E is thicker and saltier than that of the west (Figure 7a2 ). In contrast, Sections 2-6 presented a weaker spatial distinction of the halocline in the alongsection direction.
section_name: Thermohaline Structure

section_content: The water ( f W ), heat ( f H ), and salt ( f S ) fluxes can be calculated as
section_name: Water, Heat, and Salt Fluxes

section_content: The water (W f ), heat (H f ), and salt (S f ) fluxes can be calculated as [38]

where u ⊥ represents the cross-sectional flow velocity component, h is the water depth, x the section width, ρ the density of water, and C p = 4168 J/kg • • C the specific heat capacity. Table 3 shows the net water/heat/salt transport across Sections 1-6 of Figure 1c, calculated by numerical integration of the water/heat/salt flux curves shown in Figure 8. 

In Section 1, the water flux peaks at -3. 59 Sv at UTC 2020-07-19 07:00 (i. e., the peak flood of the spring tide) when the flow velocity reached its maximum value. The negative values of the water volume transport across Section 1 during the spring tide (-2. 86 × 10 11 m 3 ) and neap tide (-1. 31 × 10 11 m 3 ) indicate a residual current in the downslope direction (Table 3 ). The net water fluxes are negative in Sections 2-6 (except for Section 6 during the spring tide), demonstrating that water flows southwestward in the along-isobaths direction of the QDNS (Figure 8b-f ). Considering our previous conclusion (Section 3. 1) that the surface water flows northeastward near the QDNS, while the bottom water flows in the direction opposite to it in the basin area, we can now deduce that the net water transport in the along-isobaths direction of the QDNS is dominated by the southwestward bottom water in the basin area. Similar behavior can be seen in the hydro-thermohaline curves (Figure 8b-f ), but the water, heat, and salt transports during the neap tide gradually increase from Section 2 to Section 6 (Table 3 ) due to increasing averaged water depth. 

The water, heat, and salt transports across Section 6 during the spring tide are 1. 08 × 10 12 m 3, 4. 44 × 10 14 W, and 3. 81 × 10 13 kg, respectively. This anomaly can be understood when taking a look at the flow variations of Points A and B in Figure 1c (yellow stars). In fact, during the neap tide, the depth-averaged flow velocity there is negative (Figure 9c ) due to the dominant southwestward flow (Figure 4i ). During the spring tide, the eddy discussed in Section 3. 1 causes a constantly changing flow direction at Point B (red lines in Figure 9b, c ). In contrast, the flow at Point A (black lines in Figure 9b, c ), which is not covered by the eddy, is dominated by the surface current, resulting in a positive depth-averaged velocity. As discussed, the northeastward water, heat, and salt transports during the spring tide in Section 6 can be attributable to the eddy. 

The quantitative calculation of the net transport fluxes leads to the following conclusions. Spatially, hydro-thermohaline transport in the downslope direction (Section 1) is typically smaller than those along the isobaths (Sections 2-6). Temporally, the spring tide makes a 119%/70%/120% higher contribution to water/heat/salt transport than the neap tide in Section 1. In Sections 2-6, the net water/heat/salt transported during the spring tide was generally smaller (i. e., 38%/38%/38% on average) than that during the neap tide (Figure 10 ). 

Table 3 shows the net water/heat/salt transport across Sections 1-6 of Figure 1c, calculated by numerical integration of the water/heat/salt flux curves shown in Figure 8. In Section 1, the water flux peaks at -3. 59 Sv at UTC 2020-07-19 07:00 (i. e., the peak flood of the spring tide) when the flow velocity reached its maximum value. The negative values of the water volume transport across Section 1 during the spring tide (-2. 86 × 10 11 m 3 ) and neap tide (-1. 31 × 10 11 m 3 ) indicate a residual current in the downslope direction (Table 3 ). The net water fluxes are negative in Sections 2-6 (except for Section 6 during the spring tide), demonstrating that water flows southwestward in the along-isobaths direction of the QDNS (Figure 8b-f ). Considering our previous conclusion (Section 3. 1) that the surface water flows northeastward near the QDNS, while the bottom water flows in the direction opposite to it in the basin area, we can now deduce that the net water transport in the along-isobaths direction of the QDNS is dominated by the southwestward bottom water in the basin area. Similar behavior can be seen in the hydro-thermohaline curves (Figure 8b-f ), but the water, heat, and salt transports during the neap tide gradually increase from Section 2 to Section 6 (Table 3 ) due to increasing averaged water depth. 

The water, heat, and salt transports across Section 6 during the spring tide are 1. 08 × 10 12 m 3, 4. 44 × 10 14 W, and 3. 81 × 10 13 kg, respectively. This anomaly can be understood when taking a look at the flow variations of Points A and B in Figure 1c (yellow stars). In fact, during the neap tide, the depth-averaged flow velocity there is negative (Figure 9c ) due to the dominant southwestward flow (Figure 4i ). During the spring tide, the eddy discussed in Section 3. 1 causes a constantly changing flow direction at Point B (red lines in Figure 9b, c ). In contrast, the flow at Point A (black lines in Figure 9b, c ), which is not covered by the eddy, is dominated by the surface current, resulting in a positive depthaveraged velocity. As discussed, the northeastward water, heat, and salt transports during the spring tide in Section 6 can be attributable to the eddy. The quantitative calculation of the net transport fluxes leads to the following conclusions. Spatially, hydro-thermohaline transport in the downslope direction (Section 1) is typically smaller than those along the isobaths (Sections 2-6). Temporally, the spring tide makes a 119%/70%/120% higher contribution to water/heat/salt transport than the neap tide in Section 1. In Sections 2-6, the net water/heat/salt transported during the spring tide was generally smaller (i. e., 38%/38%/38% on average) than that during the neap tide (Figure 10 ).
section_name: Water, Heat, and Salt Fluxes

section_content: As discussed above, we have established a strong correlation between the spatial distributions and cross-sectional transport fluxes of water and salt. In this section, we discuss the mechanisms by which the water flow affects the salinity as well as the temperature. Based on the thermohaline balance equation introduced in Section 2. 2, Figures 11 and 12 show the hydrodynamic terms contributing to salt and heat fluxes, respectively, where 16-days-averaged simulation results from 14 July to 29 July 2020 have been used for the calculation.
section_name: Thermohaline Balance

section_content: As discussed above, we have established a strong correlation between the spatial distributions and cross-sectional transport fluxes of water and salt. In this section, we discuss the mechanisms by which the water flow affects the salinity as well as the temperature. Based on the thermohaline balance equation introduced in Section 2. 2, Figures 11 and 12 show the hydrodynamic terms contributing to salt and heat fluxes, respectively, where 16-daysaveraged simulation results from 14 July to 29 July 2020 have been used for the calculation. 

was generally smaller (i. e., 38%/38%/38% on average) than that during the neap tide (Figure 10 ).
section_name: Thermohaline Balance

section_content: As discussed above, we have established a strong correlation between the spatial distributions and cross-sectional transport fluxes of water and salt. In this section, we discuss the mechanisms by which the water flow affects the salinity as well as the temperature. Based on the thermohaline balance equation introduced in Section 2. 2, Figures 11 and 12 show the hydrodynamic terms contributing to salt and heat fluxes, respectively, where 16-days-averaged simulation results from 14 July to 29 July 2020 have been used for the calculation. The Eulerian mean flow term (Q 1 ) contributes the majority of the overall salt flux (~±5 m 2 /s) and heat flux (~±5 × 10 9 Jm -1 s -1 ). In contrast, the magnitude of the tidal pumping term (Q 2 ) was only one-thousandth of Q 1, indicating a negligible effect on heat and salt transports in the study domain. The Gravitational Convection term (Q 3 ) is two orders of magnitude smaller than Q 1 for salt transport. However, for heat transport, Q 3 has the same magnitude as Q 1 (Figure 12 ), explaining why the correlation between the water and salt fluxes is much stronger than between the water and heat fluxes. on the heat flux but the opposite effect on the salt flux. East of 110. 9° E, 3 Q contributes a maximum of 8. 12 × 10 9 Jm -1 s -1 to the heat flux, more than three times of the maximum 1 Q

. 

In Sections 2-6, 3 Q contributes significantly to the northeastward heat flux, with a maximum value up to 10 10 Jm -1 s -1 (Figure 12 ). 3 Q has the same effect as 1 Q on the salt flux but the opposite effect on the heat flux.
section_name: Thermohaline Balance

section_content: The hydro-thermohaline environment in the QDNS of the SCS has been investigated using the well-established FVCOM model. This study reports the short-term spatial-temporal behavior of the hydro-thermohaline characteristics and transport fluxes of water, heat and salt in the QDNS during tide cycles. The mechanisms responsible for them are studied by comparing different terms in the thermohaline balance equation. 

The surface water near the QDNS moves northeastward along the isobaths, while the bottom water in the basin area, which dominates the net flux of the water transport, moves southwestward. In the downslope direction, water, heat, and salt transports during the spring tide are substantially larger (119%, 70% and 120% higher, respectively) than during In Section 1 of Figure 1c, Q 1 constitutes a negative (downslope) contribution to the thermohaline flux west of 110. 9 • E and a positive (upslope) one east of it (Figures 11a and 12a ), similar to the flow velocity distribution in Figure 5a1. Q 3 has the same effect as Q 1 on the heat flux but the opposite effect on the salt flux. East of 110. 9 • E, Q 3 contributes a maximum of 8. 12 × 10 9 Jm -1 s -1 to the heat flux, more than three times of the maximum Q 1. 

In Sections 2-6, Q 3 contributes significantly to the northeastward heat flux, with a maximum value up to 10 10 Jm -1 s -1 (Figure 12 ). Q 3 has the same effect as Q 1 on the salt flux but the opposite effect on the heat flux.
section_name: Conclusions

section_content: The hydro-thermohaline environment in the QDNS of the SCS has been investigated using the well-established FVCOM model. This study reports the short-term spatialtemporal behavior of the hydro-thermohaline characteristics and transport fluxes of water, heat and salt in the QDNS during tide cycles. The mechanisms responsible for them are studied by comparing different terms in the thermohaline balance equation. 

The surface water near the QDNS moves northeastward along the isobaths, while the bottom water in the basin area, which dominates the net flux of the water transport, moves southwestward. In the downslope direction, water, heat, and salt transports during the spring tide are substantially larger (119%, 70% and 120% higher, respectively) than during the neap tide, while the opposite is true (62%, 62% and 62% lower, respectively) for transports along the isobaths. A stronger cross-isobaths water exchange during the neap tide is responsible for the former, while the latter is caused by dominating southwestward flow and the motion of a subsurface cyclonic eddy. Moreover, salt transport generally strongly correlates with water transport since the mean flow dominates the salt balance equation. In addition, for heat transport, gravitational convection also plays a substantial role in the heat balance, weakening the correlation between transports of heat and water. 

The data of the short-term thermohaline evolution of the QDNS provided in this study may be of use for future studies of the northern SCS, including its marine ecology environment protection, marine pollution prevention, and ocean ranching development.
section_name: Conclusions

section_name: Acknowledgments:

section_content: We appreciate Samuel Ukpong Okon 's comments on the manuscript and the formatting of the figures.
section_name: Acknowledgments:

section_content: Funding: This research was funded by the National Key Research and Development Program of China : No. 2021YFE0206200 ; National Key Research and Development Program of China : No. 2017YFC0305905 ; Key Research and Development Project of Hainan Province : No. ZDYF2022SHFZ045. Institutional Review Board Statement: Not applicable. Informed Consent Statement: Not applicable.

section_content: Data Availability Statement: Not applicable.

section_name: Conflicts of Interest:

section_content: The authors declare no conflict of interest.
section_name: Conflicts of Interest: