High-Resolution Observations of Upwelling and Front in Daya Bay, South China Sea

: Field observations of coastal regions are important for studying physical and biological features. Observations of high-resolution coastal phenomena were obtained by using a tow-yo instrument and a turbulence proﬁler at Daya Bay in the South China Sea in October 2015. Details of coastal phenomena, including warm water from a nuclear plant discharge, as well as an upwelling, and front, were obtained. The upwelling, with a width of 2 km, resulted in saltier and more turbid water near the bottom, with low chlorophyll-a and dissolved oxygen contents being transported upward to the surface layer and changing the local water environment. The front, with the lateral salinity variations as large as 0.7 psu across 1 km, was active at the water intersection of the South China Sea and Daya Bay. Such events commonly form during weak stratiﬁcation periods in autumn. Continuous measurements from VMP-250 proﬁler over circa 22 h revealed active fronts and an averaged dissipation rate of 8 × 10 − 8 W/kg and diffusivity of 5.8 × 10 − 5 m 2 /s (i.e., one order of magnitude larger than in the open ocean) in the thermocline. The front was accompanied by strong mixing, indicating that it had formed at the intersection of different water masses and played an important role in energy dissipation in Daya Bay, further affecting the distribution of ecological elements.


Introduction
The South China Sea (SCS) is the largest marginal sea connected to the Western Pacific and is dominated by the East Asian monsoon. In winter, northeasterly winds overlie the SCS, whereas in summer, southwesterly winds prevail across the region [1,2]. In the Northern SCS, surface currents (e.g., the Kuroshio Current, South China Sea Warm Current, and cyclonic currents) display complicated dynamics. All of these large-scale currents are generally associated with wind, wind stress curl, and the processes from the Pacific [3][4][5][6][7].
Mesoscale to small-scale structures prevail in shelf seas of the Northern SCS, including eddies, fronts, and upwellings/downwellings [8]. These features play very important roles in regulating oceanic heat and energy exchange, are vital to the transport of biochemicals, and influence interannual/decadal climate variations [9]. Daya Bay (DYB) is located in the subtropical ocean and is one of the largest and most important gulfs along the southern coast of China. Physical and biochemical processes in DYB have been the focus of many previous studies. Thermal pollution has been identified in the vicinity of a nuclear power The main data used in the study were collected by a YODA profiler and a VMP-250 turbulence profiler.
The YODA profiler (manufactured by JFE Advantech Co. Ltd., Japan; Figure 1a) is a tow-yo instrument that can be used to profile the water column with high resolution from small boats, without occupying much space [11]. The instrument is equipped with a deployment winch and sensors that measure conductivity, temperature, pressure, chlorophyll, turbidity, and DO. The brush at the top of the instrument provides a stabilizing effect on the free-fall sinking speed, which is approximately constant at 0.2 m/s. The sampling frequency is 10 Hz. All data are stored internally and downloaded to a PC through a wet connector and interface. designed for operation off small vessels with limited deck space (e.g., zodiacs) or where electrical power facilities are limited or missing (e.g., ice camps). The instrument is equipped with two shear probes, one thermistor, and a nose-mounted conductivitytemperature sensor. The VMP-250 has a constant free-fall sinking speed of approximately 0.7 m/s. The sampling frequency is 1024 Hz. All data are stored internally and downloaded to a PC through a wet connector and interface.

Observations
DYB is a large semi-enclosed bay (ca. 24 km in length and 13-25 km in width, covering an area of ca. 600 km 2 ). The bay has an irregular coastline and contains more than 50 islands. It is located along the south coast of China (114°29.7′-114°49.7′ E, 22°31.2′-22°50.0′ N), as shown in Figure 2. Measurements were taken between 20 October and 22 October 2015. On 20 October 20 2015, we obtained conductivity, temperature, pressure, chlorophyll, turbidity, and DO data, using YODA along the western boundary of DYB (S1) and a cross-section of the mouth of the bay (S2) (Figure 2). The S1 observation transect ran from the north to the south, with the survey starting at 12:20 local time and ending at 14:50, whereas S2 ran from the west to the east, with the survey starting at 16:00 and ending at 18:30. During the cruise, YODA obtained 53 and 56 profiles at S1 and S2, respectively. The average space interval between successive YODA measurements was approximately 230 m, enabling fine-scale observations of physical phenomena. The VMP-250 (manufactured by Rockland Scientific inc., Canada; Figure 1b) is a coastal-zone profiler for the measurement of microscale turbulence. The profiler is designed for operation off small vessels with limited deck space (e.g., zodiacs) or where electrical power facilities are limited or missing (e.g., ice camps). The instrument is equipped with two shear probes, one thermistor, and a nose-mounted conductivity-temperature sensor. The VMP-250 has a constant free-fall sinking speed of approximately 0.7 m/s. The sampling frequency is 1024 Hz. All data are stored internally and downloaded to a PC through a wet connector and interface.

Observations
DYB is a large semi-enclosed bay (ca. 24 km in length and 13-25 km in width, covering an area of ca. 600 km 2 ). The bay has an irregular coastline and contains more than 50 islands. It is located along the south coast of China (114 • 29.7 -114 • 49.7 E, 22 • 31.2 -22 • 50.0 N), as shown in Figure 2. Measurements were taken between 20 October and 22 October 2015. On 20 October 20 2015, we obtained conductivity, temperature, pressure, chlorophyll, turbidity, and DO data, using YODA along the western boundary of DYB (S1) and a cross-section of the mouth of the bay (S2) (Figure 2). The S1 observation transect ran from the north to the south, with the survey starting at 12:20 local time and ending at 14:50, whereas S2 ran from the west to the east, with the survey starting at 16:00 and ending at 18:30. During the cruise, YODA obtained 53 and 56 profiles at S1 and S2, respectively. The average space interval between successive YODA measurements was approximately 230 m, enabling fine-scale observations of physical phenomena.
The VMP-250 profiler was used to measure turbulence at a fixed station at the mouth of DYB (yellow triangle in Figure 2b). We obtained 74 turbulence profiles from 10:00 on 21 October to 7:30 on 22 October 2015. The time intervals between VMP-250 measurements ranged from 10 to 30 min; therefore, any rapid changes in physical phenomena could be observed. The VMP-250 profiler was used to measure turbulence at a fixed station at the mouth of DYB (yellow triangle in Figure 2b). We obtained 74 turbulence profiles from 10:00 on 21 October to 7:30 on 22 October 2015. The time intervals between VMP-250 measurements ranged from 10 to 30 min; therefore, any rapid changes in physical phenomena could be observed.
There are no major rivers discharged into DYB, but three small streams discharged into Dapeng Cove. From September to October 2015, the average rainfall in DYB is less than 30 mm. The sea surface wind was northeasterly and was less than 10 m/s during the observation period, indicating that our observation was less affected by surface runoff and sea surface wind field. Therefore, the effects of rainfall, runoff and sea surface wind are ignored in the latter analysis. The tidal current in DYB is characterized by diurnal tide and semidiurnal tide and the semidiurnal tide is more prominent (Song et al., 2016).

Calculation of Turbulence
The shear probes in the VMP-250 profiler measured high-frequency velocity fluctuations, which were used to calculate the local turbulent dissipation rate. The dissipation rate was estimated by using the following isotropic formula: where is the kinematic molecular viscosity, the overline indicates a spatial average, u is either one of the two horizontal components of velocity, z is the vertical coordinate, ( ) There are no major rivers discharged into DYB, but three small streams discharged into Dapeng Cove. From September to October 2015, the average rainfall in DYB is less than 30 mm. The sea surface wind was northeasterly and was less than 10 m/s during the observation period, indicating that our observation was less affected by surface runoff and sea surface wind field. Therefore, the effects of rainfall, runoff and sea surface wind are ignored in the latter analysis. The tidal current in DYB is characterized by diurnal tide and semidiurnal tide and the semidiurnal tide is more prominent ( [18]).

Calculation of Turbulence
The shear probes in the VMP-250 profiler measured high-frequency velocity fluctuations, which were used to calculate the local turbulent dissipation rate. The dissipation rate was estimated by using the following isotropic formula: where ν is the kinematic molecular viscosity, the overline indicates a spatial average, u is either one of the two horizontal components of velocity, z is the vertical coordinate, ψ(k) is the spectrum of vertical shear, and k is the vertical wavenumber [23]. The upper limit of the integration, k max , is variable. The k max value was calculated by using the method recommended by Shang et al. [24], which states that it should be the smallest of the following values: the lowest frequency that shows corruption of the shear signal by vibrations, a wavenumber of 150 cpm due to the spatial resolution of the shear probe, the cutoff frequency of the antialiasing filter, an estimate of the wavenumber that resolves 90% of the shear variance according to the Nasmyth spectrum, using the method of Lueck [25], or the location of the spectral minimum determined with a low-order fit to the spectrum in log-log space. Diapycnal diffusivity [23] was calculated based on the dissipation rate (ε) and stratification as follows: where Γ is the mixing efficiency, which was set to 0.2. The stratification was calculated as follows: N 2 = − g ρ ∂ρ ∂z , where g is the acceleration due to gravity. Figure 3 provides a workflow graph to recap on the logic path of the study. The workflow mainly consists of three parts, including the used instruments, the data types, and the observed costal phenomena. spectrum in log-log space. Diapycnal diffusivity [23] was calculated based on the dissipation rate (ε) and stratification as follows:

A Workflow Graph
where Γ is the mixing efficiency, which was set to 0.2. The stratification was calculated as , where g is the acceleration due to gravity. Figure 3 provides a workflow graph to recap on the logic path of the study. The workflow mainly consists of three parts, including the used instruments, the data types, and the observed costal phenomena.

High-Resolution Observations by YODA
In the section, features of coastal phenomenon, including the nuclear power station discharged warm water, upwelling and front, are revealed by high spatial resolution observations by YODA.

Water Properties
Our efforts to obtain high-resolution data for DYB by using the YODA profiler proved to be successful. Figures 4 and 5 show the observed data, which indicated that the water environment of DYB was complex, with obvious horizontal spatial and vertical differences in temperature and salinity characteristics. The temperature varied from 26.2 to 27.4 °C, and salinity varied from 32.5 to 33.2 psu. The water along transect S1 was warm and fresh to the north, cold and salty in the middle, and warm and fresh farther to the south. The water was warm and salty to the west, and cold and fresh to the east. The temperature-salinity (T-S) sections (Figures 4 and 5) and T-S diagrams ( Figure 6) indicate that there were at least five distinct water masses along transect S1 and three distinct water masses along transect S2. The different water masses might have been formed by different dynamic characteristics.

High-Resolution Observations by YODA
In the section, features of coastal phenomenon, including the nuclear power station discharged warm water, upwelling and front, are revealed by high spatial resolution observations by YODA.

Water Properties
Our efforts to obtain high-resolution data for DYB by using the YODA profiler proved to be successful. Figures 4 and 5 show the observed data, which indicated that the water environment of DYB was complex, with obvious horizontal spatial and vertical differences in temperature and salinity characteristics. The temperature varied from 26.2 to 27.4 • C, and salinity varied from 32.5 to 33.2 psu. The water along transect S1 was warm and fresh to the north, cold and salty in the middle, and warm and fresh farther to the south. The water was warm and salty to the west, and cold and fresh to the east. The temperature-salinity (T-S) sections (Figures 4 and 5) and T-S diagrams ( Figure 6) indicate that there were at least five distinct water masses along transect S1 and three distinct water masses along transect S2. The different water masses might have been formed by different dynamic characteristics.       The different water masses along transect S1 could be separated into the following five distinct areas. Area 1 had a latitude greater than 22.56 • N, with warmer and fresher water found at the surface. Water in this area was affected by thermal discharges from the nuclear power station. Area 2 was found between the latitudes 22.525 • N and 22.535 • N, with saltier water present throughout the whole depth of the water column and no obvious changes in temperature with depth. Area 3 was found between the latitudes 22.505 • N and 22.525 • N, with warmer water present near the surface and no obvious changes in salinity with depth. Area 4 had a latitude lower than 22.505 • N, with cooler water present throughout the whole depth of the water column. Area 5 was found between the latitudes 22.535 • N and 22.56 • N and represented a transition region. The different water masses along transect S1 could be separated into the following five distinct areas. Area 1 had a latitude greater than 22.56° N, with warmer and fresher water found at the surface. Water in this area was affected by thermal discharges from the nuclear power station. Area 2 was found between the latitudes 22.525° N and 22.535° N, with saltier water present throughout the whole depth of the water column and no obvious changes in temperature with depth. Area 3 was found between the latitudes 22.505°N and 22.525° N, with warmer water present near the surface and no obvious changes in salinity with depth. Area 4 had a latitude lower than 22.505° N, with cooler water present throughout the whole depth of the water column. Area 5 was found between the latitudes 22.535° N and 22.56° N and represented a transition region.
The different water masses along transect S2 could be separated into the following three distinct areas. Area 1 had a longitude greater than 114.7° E, with cold and fresher water in the upper layer. Area 2 was found between the longitudes 114.655° E and 114.7° E, with saltier water present near the bottom and no obvious changes in temperature with depth. Area 3 had a longitude less than 114.655° E, with warmer and fresher water present near the surface. Generally, water on a shelf is cooler and fresher than in the adjoining ocean, leading to the formation of a shelf break front [26]. Therefore, it was reasonable to assume that the cold water at longitudes greater than 114.7° E originated from the bay. The area between 114.655° E and 114.7° E was considered a transition region.

Influence of the Nuclear Power Station
Although the observation area was small, the water characteristics of DYB exhibited obvious spatial variability over a distance of approximately 12.5 km.
A significant feature of water temperature was that it was clearly affected by the thermal discharge from the nuclear power station in the north, and the temperature was significantly higher in the north than in other places ( Figure 4a). The surface water temperature was as high as 27.4 °C in the northernmost observation position, which was approximately 2 km from the nuclear power station. This was 1 °C higher than the water temperature in the south. The warm water from the thermal discharge was mainly concentrated in the upper 10 m and extended approximately 5 km to the south of the release orifice, which corresponded to 22.56° N (Figure 4a). This range at which water discharge from the nuclear power station influenced water temperature was consistent with the satellite remote-sensing study of [10]. It has been reported that during the winter months, the thermal plume formed by water discharge from the nuclear power station is localized to an area within a few kilometers of the power plant. During the summer months, there is a larger thermal plume that extends 8-10 km south along the coast [10]. The different water masses along transect S2 could be separated into the following three distinct areas. Area 1 had a longitude greater than 114.7 • E, with cold and fresher water in the upper layer. Area 2 was found between the longitudes 114.655 • E and 114.7 • E, with saltier water present near the bottom and no obvious changes in temperature with depth. Area 3 had a longitude less than 114.655 • E, with warmer and fresher water present near the surface. Generally, water on a shelf is cooler and fresher than in the adjoining ocean, leading to the formation of a shelf break front [26]. Therefore, it was reasonable to assume that the cold water at longitudes greater than 114.7 • E originated from the bay. The area between 114.655 • E and 114.7 • E was considered a transition region.

Influence of the Nuclear Power Station
Although the observation area was small, the water characteristics of DYB exhibited obvious spatial variability over a distance of approximately 12.5 km.
A significant feature of water temperature was that it was clearly affected by the thermal discharge from the nuclear power station in the north, and the temperature was significantly higher in the north than in other places ( Figure 4a). The surface water temperature was as high as 27.4 • C in the northernmost observation position, which was approximately 2 km from the nuclear power station. This was 1 • C higher than the water temperature in the south. The warm water from the thermal discharge was mainly concentrated in the upper 10 m and extended approximately 5 km to the south of the release orifice, which corresponded to 22.56 • N (Figure 4a). This range at which water discharge from the nuclear power station influenced water temperature was consistent with the satellite remote-sensing study of [10]. It has been reported that during the winter months, the thermal plume formed by water discharge from the nuclear power station is localized to an area within a few kilometers of the power plant. During the summer months, there is a larger thermal plume that extends 8-10 km south along the coast [10]. The strong seasonal difference in the thermal plume is related to vertical mixing in the water column in winter and stratification in summer. In addition to the warm water being discharged from the nuclear power station, a cloud of warm water with a radius of 2 km at the surface between 22.5 • N and 22.52 • N was also observed. This area of water was 0.3 • C warmer than the surrounding water. The temperature difference between the surface and bottom was less than 1 • C, indicating weak stratification of the water column during the observation period.
The salinity and chlorophyll-a content were also affected by the power plant. In the area influenced by discharge water, the salinity ( Figure 4b) and chlorophyll-a content ( Figure 7c) were lower than in the south. There was no significant difference in DO (Figure 7d).
water column in winter and stratification in summer. In addition to the warm water being discharged from the nuclear power station, a cloud of warm water with a radius of 2 km at the surface between 22.5° N and 22.52° N was also observed. This area of water was 0.3 °C warmer than the surrounding water. The temperature difference between the surface and bottom was less than 1°C, indicating weak stratification of the water column during the observation period.
The salinity and chlorophyll-a content were also affected by the power plant. In the area influenced by discharge water, the salinity ( Figure 4b) and chlorophyll-a content ( Figure 7c) were lower than in the south. There was no significant difference in DO ( Figure  7d).

The Front
Water in DYB is cooler and fresher than in the adjoining SCS, resulting in a strong salinity front that is the most notable hydrographic feature in this region ( Figure 4). There was a lateral density gradient identified along both transects S1 and S2. For transect S1, YODA measurements revealed lateral salinity variations as large as 0.7 psu across 1 km between 22.53 • N and 22.54 • N. Transect S2 also had a similar variation in salinity between 114.65 • E and 114.66 • E.
The width of the front was approximately 1 km, and the water depth, h, was 20 m. The Rossby deformation radius R = g h/ f was approximately 6 km. The width of the front was much less than the Rossby deformation radius, indicating that it was not maintained by geostrophic balance. Here, g is the reduced gravity, which is given by g = g(r 2 − r 1 )/r, where r 1 and r 2 are the densities of the upper and lower layers, respectively, r is the mean density, g = 9.81 m/s 2 , and f is the Coriolis parameter.
In the front area along transect S1, chlorophyll-a increased from 3 to 8 g/L, and the turbidity decreased from 4 to 1 FTU. We assumed that the strong front in S1 was formed at the intersection of different water masses and was strengthened by an upwelling system, in which cooler and saltier water at the bottom was dragged to the upper layer (Section 3.1.4). The presence of the front also prevented water from the adjoining SCS from entering the inner area of DYB.

The Upwelling
A cold salty water mass dominated the area between 22.52 • N and 22.54 • N. The turbidity was relatively high from the bottom to the sea surface, whereas the DO and chlorophyll-a concentrations were relatively low. The distributions of the parameters indicated a pumping process. The width of the upwelling was approximately 2 km. The most direct evidence of upwelling was the turbidity signal. Figure 7e shows that the turbidity was much higher at depth than in the upper layer. It was induced by friction between water layers, making the water in the bottom boundary layer turbid. In the area between 22.52 • N and 22.54 • N, the whole water column was more turbid than the surrounding water. This could only be explained by the turbidity of the entire water body being derived from the turbid water body at the bottom.
The occurrence of an upwelling is closely related to wind conditions. In the Northern Hemisphere, when the wind blows parallel to the coast for a long time, under the action of a geostrophic deflection force, the wind drift formed by the wind causes the surface water to leave the coast. This in turn causes the shallower water near the shore to rise and form an upwelling. In the North SCS, the most notable upwellings are the Eastern Guangdong upwelling and the Eastern Hainan upwelling, which occur in summer months [27][28][29]. In some case, the wind in particular direction acting on the vorticity field can generate vertical motions. For example, Chang et al. [30] reported that the seaward increasement of cyclonic vorticity near the Kuroshio's western edge favors a stronger Ekman transport away from the jet, producing upwelling at the shelfbreak under a northeasterly wind. The passage of high-energy atmospheric systems also may stir up the upwelling [31]. In this study, however, the upwelling mechanism described above cannot explain the upwelling observed here in October. The upwelling areas caused by large-scale wind fields are more extensive than the upwelling area observed in this study, which had a width of approximately 2 km. The sea surface wind field far from the high-energy atmospheric system during the observation period is also not a cause of the vertical motions of water.
When the surface water is blocked by the coast or an island, sea water will gather and diverge in the horizontal direction, and an upwelling may occur in the vertical direction. In the study region, the terrain is flat in the area where the upwelling was observed (Figure 1), which indicates that the upwelling was induced by other dynamic mechanisms.
We assumed that the upwelling was formed by the interaction between tidal currents and the topography. Upwelling can bring nutrients from the bottom layer to the surface layer, make plankton grow in large quantities, and provide food for fish. Therefore, upwelling areas are often important fishing grounds, such as the Peru fishing grounds, which benefit from the cold Peru Current (upwelling compensation current). In our study area, the upwelling may also have an important impact on the local environment. As shown in Figures 7 and 8, the upwelling changed the distribution of chlorophyll-a and DO in the local area and carried turbid material from the bottom up to the surface. and the topography. Upwelling can bring nutrients from the bottom layer to the surface layer, make plankton grow in large quantities, and provide food for fish. Therefore, upwelling areas are often important fishing grounds, such as the Peru fishing grounds, which benefit from the cold Peru Current (upwelling compensation current). In our study area, the upwelling may also have an important impact on the local environment. As shown in Figures 7 and 8, the upwelling changed the distribution of chlorophyll-a and DO in the local area and carried turbid material from the bottom up to the surface. Figure 8. (a) Temperature, salinity, (b) chlorophyll-a, dissolved oxygen, and turbidity at a depth of 3 m along transect S1. The filled backgrounds indicate the upwelling (gray), front (orange), and warm discharge water from the nuclear reactor (pink).

High-Resolution Observations by the VMP-250 Profiler
In the section, feature of coastal phenomenon of front and its accompanied turbulence are revealed by high temporal resolution observations by YODA.

Water Properties
We used the VMP-250 profiler to observe turbulence for 22 h at a fixed location and obtained 74 profiles. The time intervals varied between 10 and 30 min (Figure 9c). A nosemounted conductivity-temperature sensor was installed on the VMP-250 profiler, enabling the analysis of temperature and salinity characteristics observed at continuous monitoring stations.
The measurements indicated that both the temperature and salinity changed rapidly with time. Using density data, we calculated the depth of the surface mixed layer and the thickness of the bottom mixed layer (Figure 9c). The depth of the mixed layer was defined as the point at which the density varied by less than 0.05 kg/m 3 . The depth of the upper mixed layer was 8 m, whereas the thickness of the bottom mixed layer extended to 8 m (the water depth at the observation station was ca. 20 m). Sometimes, a bottom mixed layer could not be discerned. Activities in the upper water extended to the sea bottom at times (e.g., at 20:00 on 21 October 2015). The bottom mixed layer was mainly affected by tides, with the mixed layer depth increasing at high tide (Figure 9c,f). The largest buoyancy fluxes were in the stratified region just above the bottom mixed layer ( Figure  9d).
The T-S sections and T-S diagrams ( Figure 10) indicate that the water masses at various stations comprised distinct water masses. For example, between 21:00 on 21

High-Resolution Observations by the VMP-250 Profiler
In the section, feature of coastal phenomenon of front and its accompanied turbulence are revealed by high temporal resolution observations by YODA.

Water Properties
We used the VMP-250 profiler to observe turbulence for 22 h at a fixed location and obtained 74 profiles. The time intervals varied between 10 and 30 min (Figure 9c). A nose-mounted conductivity-temperature sensor was installed on the VMP-250 profiler, enabling the analysis of temperature and salinity characteristics observed at continuous monitoring stations.
The measurements indicated that both the temperature and salinity changed rapidly with time. Using density data, we calculated the depth of the surface mixed layer and the thickness of the bottom mixed layer (Figure 9c). The depth of the mixed layer was defined as the point at which the density varied by less than 0.05 kg/m 3 . The depth of the upper mixed layer was 8 m, whereas the thickness of the bottom mixed layer extended to 8 m (the water depth at the observation station was ca. 20 m). Sometimes, a bottom mixed layer could not be discerned. Activities in the upper water extended to the sea bottom at times (e.g., at 20:00 on 21 October 2015). The bottom mixed layer was mainly affected by tides, with the mixed layer depth increasing at high tide (Figure 9c,f). The largest buoyancy fluxes were in the stratified region just above the bottom mixed layer (Figure 9d).
The T-S sections and T-S diagrams ( Figure 10) indicate that the water masses at various stations comprised distinct water masses. For example, between 21:00 on 21 October 2015 and 3:00 on 22 October 2015, the water temperature was low and salinity was high, especially in the bottom mixed layer. This was a period of high tide, indicating that SCS water was dominant during this period. From 6:00 to 8:00 on 22 October 2015, the upper layer of water exhibited a low temperature and low salinity. This was the ebb tide period, indicating that DYB water was the dominant water mass. Compared with the Figure 6, the characteristics of the seawater was confirmed. The observation station was located at the confluence of the SCS and DYB, which changed constantly under the action of tides. With the interactions at the confluence, an abundance of thermohaline structures was formed.
was high, especially in the bottom mixed layer. This was a period of high tide, indicating that SCS water was dominant during this period. From 6:00 to 8:00 on 22 October 2015, the upper layer of water exhibited a low temperature and low salinity. This was the ebb tide period, indicating that DYB water was the dominant water mass. Compared with the Figure 6, the characteristics of the seawater was confirmed. The observation station was located at the confluence of the SCS and DYB, which changed constantly under the action of tides. With the interactions at the confluence, an abundance of thermohaline structures was formed.  Isopycnal contours (solid lines) in (a-e) represent density. The front is marked with an inverted triangle in (a). In (c-e), the red line indicates the surface mixed layer, whereas the yellow line indicates the bottom mixed layer. Vertical dotted lines in (f) indicate the time points of VMP-250 observations. The water-level data in (f,g) were derived at the tidal elevation station T1.

The Front
In addition to YODA measurements, lateral density gradients were also identified in the study area from the fixed station observations (Figure 9). These gradients were controlled by variations in both salinity and temperature. Stratification persisted throughout the upper water column, with considerable variability along isopycnal. VMP-250 measurements indicated lateral temperature variations as large as 0.4 °C and lateral salinity variations as large as 0.3 psu over 30-min periods at 11:30 and 20:30 on 21 October 2015. It is noted that the conventional definition of density front (dρ/dx) is related to the variation of horizontal distance of density. Here, we assume that the density front does not change in a certain period of time when it moves in the flow direction, dρ/dt = Udρ/dx. Thus, the larger the variation in dρ/dt, the stronger the density front. There were two strong fronts at these two moments, with density variations of approximately 0.2 kg/m 3 over a 30-min period. At 20:30 on 21 October 2015, the front extended to the bottom. The width of the front was approximately 1 km (Section 3.1.3), and it took approximately 30 min for the front to cross the observation position. Therefore, we inferred that the speed of the front was approximately 0.5 m/s, which is similar to the tidal current velocity [18].

Turbulence
Examples of dissipation spectra at specified depth segments from the VMP-250 measurement at 12:15 on 21 October 2015 are shown in Figure 11. There were two shear probes installed at the head of profiler, and data from both probes are shown. The spectra were approximately consistent with the Nasmyth universal spectrum, indicating the credibility of the measurement. The triangles indicate the upper limit of integration. The raw microscale velocity shear was weak between 4 and 8 m; therefore, the estimated dissipation rate was small (ε~10 −7 W/kg). The microscale velocity shear strengthened between 8 and 12 m, as did the dissipation (ε~10 −6 W/kg).

The Front
In addition to YODA measurements, lateral density gradients were also identified in the study area from the fixed station observations (Figure 9). These gradients were controlled by variations in both salinity and temperature. Stratification persisted throughout the upper water column, with considerable variability along isopycnal. VMP-250 measurements indicated lateral temperature variations as large as 0.4 • C and lateral salinity variations as large as 0.3 psu over 30-min periods at 11:30 and 20:30 on 21 October 2015. It is noted that the conventional definition of density front (dρ/dx) is related to the variation of horizontal distance of density. Here, we assume that the density front does not change in a certain period of time when it moves in the flow direction, dρ/dt = Udρ/dx. Thus, the larger the variation in dρ/dt, the stronger the density front. There were two strong fronts at these two moments, with density variations of approximately 0.2 kg/m 3 over a 30-min period. At 20:30 on 21 October 2015, the front extended to the bottom. The width of the front was approximately 1 km (Section 3.1.3), and it took approximately 30 min for the front to cross the observation position. Therefore, we inferred that the speed of the front was approximately 0.5 m/s, which is similar to the tidal current velocity [18].

Turbulence
Examples of dissipation spectra at specified depth segments from the VMP-250 measurement at 12:15 on 21 October 2015 are shown in Figure 11. There were two shear probes installed at the head of profiler, and data from both probes are shown. The spectra were approximately consistent with the Nasmyth universal spectrum, indicating the credibility of the measurement. The triangles indicate the upper limit of integration. The raw microscale velocity shear was weak between 4 and 8 m; therefore, the estimated dissipation rate was small (ε~10 −7 W/kg). The microscale velocity shear strengthened between 8 and 12 m, as did the dissipation (ε~10 −6 W/kg). It has been reported that dissipation and mixing on the shelf are mainly controlled by near-bottom friction and shear from low-frequency internal waves and solibores [26]. However, according to our observations, the front that formed at the intersection of the SCS and DYB also plays an important role in local energy dissipation. Figure 9c shows the dissipation rates and associated diapycnal diffusivities for each profile. The dissipation rates of the thermocline and bottom boundary layer were of the same order of magnitude. The thickness of the bottom mixing layer has the characteristics of periodic variation, and its increase is consistent with the period of high tide. The dissipation rate at the bottom also increase during the high tide periods (Figure 9e). Taken together with Figure 9c,e,f, it was concluded that the turbulence in the bottom mixed layer was driven by tidal currents. In the upper mixing layer, ocean mixing is mainly controlled by buoyancy flux and wind momentum flux [32]. In this study, we focused on the dissipation in the thermocline.
The dissipation rate of the thermocline was mainly between 1 × 10 −9 and 1 × 10 −6 W/kg, with an average value of 8 × 10 −8 W/kg. The maximum value was mainly recorded at 11:30, 16:00, and 20:30 on 21 October 2015, as well as 3:00 on 22 October 2015, at which point the density profile changed substantially. Peak dissipation at 5:00 on 22 October 2015 was induced in a patch near the bottom. Average diffusivity profiles were calculated from the average dissipation rate and stratification. Because diffusivity was not well defined in the well-mixed boundary layers, average diffusivity profiles are only shown for data excluding the boundary layer (Figure 12a). The average diffusivity associated with midcolumn, stratified turbulence was 5.8 × 10 −5 m 2 /s for the whole profile, which is one order of magnitude greater than the level in an open-ocean thermocline [33,34]. It has been reported that dissipation and mixing on the shelf are mainly controlled by near-bottom friction and shear from low-frequency internal waves and solibores [26]. However, according to our observations, the front that formed at the intersection of the SCS and DYB also plays an important role in local energy dissipation. Figure 9c shows the dissipation rates and associated diapycnal diffusivities for each profile. The dissipation rates of the thermocline and bottom boundary layer were of the same order of magnitude. The thickness of the bottom mixing layer has the characteristics of periodic variation, and its increase is consistent with the period of high tide. The dissipation rate at the bottom also increase during the high tide periods (Figure 9e). Taken together with Figure 9c,e,f, it was concluded that the turbulence in the bottom mixed layer was driven by tidal currents. In the upper mixing layer, ocean mixing is mainly controlled by buoyancy flux and wind momentum flux [32]. In this study, we focused on the dissipation in the thermocline.
The dissipation rate of the thermocline was mainly between 1 × 10 −9 and 1 × 10 −6 W/kg, with an average value of 8 × 10 −8 W/kg. The maximum value was mainly recorded at 11:30, 16:00, and 20:30 on 21 October 2015, as well as 3:00 on 22 October 2015, at which point the density profile changed substantially. Peak dissipation at 5:00 on 22 October 2015 was induced in a patch near the bottom. Average diffusivity profiles were calculated from the average dissipation rate and stratification. Because diffusivity was not well defined in the well-mixed boundary layers, average diffusivity profiles are only shown for data excluding the boundary layer (Figure 12a). The average diffusivity associated with mid-column, stratified turbulence was 5.8 × 10 −5 m 2 /s for the whole profile, which is one order of magnitude greater than the level in an open-ocean thermocline [33,34]. To analyze the relationship between front and turbulent mixing in the thermocline, the density of water was averaged, excluding the boundary layer. Changes in average density with time (dρ/dt) were calculated. Active fronts existed on 21 November but there were fewer fronts on 22 November (Figure 12b). The fronts detected between 15:00 and 18:00 on 21 October 2015 were affected by the variation in density near the surface mixed layer. It can be seen from Figure 10 that the dissipation rate and variation in dρ/dt were positively correlated, i.e., a larger dρ/dt tended to produce a larger dissipation rate. This indicates that a density front in the study area is important for locating the sources of energy dissipation.
The turbulent mixing process in the ocean determines the redistribution of seawater characteristics. This plays a key role in the redistribution of nutrients, DO, and other ecological elements and substances in the ocean, which then affects marine productivity [35][36][37][38]. In this study, we did not make observations of ecological elements when VMP-250 measurements were recorded. However, it was apparent that the front regulated the distribution of ecological elements (Figure 7), while turbulent mixing also played an important role.

Comparison with Observations in July 2019
Daya Bay is located in the tropical monsoon region. In this area, the northeast monsoon is strong in winter, whereas the weaker southwest monsoon is dominant in summer [39]. During winter, the water in DYB is vertically mixed under the influence of the northeast monsoon [40]. Water stratification in DYB varies seasonally. The development of physical phenomena in the ocean is closely related to stratification, and seasonal variation in the formation of a front and its effect on energy dissipation in the observation area was obvious.
To confirm the seasonal variation in the front in DYB, we analyzed observations of the area (114.68° E, 22.58° N) from July 2019 (Figure 1). At the observation site, 138 profiles were obtained by YODA within 48 h, with a time interval of 20 min. The temperature difference between the surface and bottom was more than 6 °C, and the salinity difference was more than 2 psu (Figure 13), which were much larger than the values measured in October ( Figure 10). The fresh water input during the summer rainy season diminished the surface salinity, and water with a lower temperature and high salinity intruded into the bay along the bottom from the SCS. Unlike the observations made in October, there was a stable bottom mixed layer in July that was more than 10 m thick (Figure 13d). Above the bottom mixing layer was a thermocline reflecting strong stratification. In the July 2019 To analyze the relationship between front and turbulent mixing in the thermocline, the density of water was averaged, excluding the boundary layer. Changes in average density with time (dρ/dt) were calculated. Active fronts existed on 21 November but there were fewer fronts on 22 November (Figure 12b). The fronts detected between 15:00 and 18:00 on 21 October 2015 were affected by the variation in density near the surface mixed layer. It can be seen from Figure 10 that the dissipation rate and variation in dρ/dt were positively correlated, i.e., a larger dρ/dt tended to produce a larger dissipation rate. This indicates that a density front in the study area is important for locating the sources of energy dissipation.
The turbulent mixing process in the ocean determines the redistribution of seawater characteristics. This plays a key role in the redistribution of nutrients, DO, and other ecological elements and substances in the ocean, which then affects marine productivity [35][36][37][38]. In this study, we did not make observations of ecological elements when VMP-250 measurements were recorded. However, it was apparent that the front regulated the distribution of ecological elements (Figure 7), while turbulent mixing also played an important role.

Comparison with Observations in July 2019
Daya Bay is located in the tropical monsoon region. In this area, the northeast monsoon is strong in winter, whereas the weaker southwest monsoon is dominant in summer [39]. During winter, the water in DYB is vertically mixed under the influence of the northeast monsoon [40]. Water stratification in DYB varies seasonally. The development of physical phenomena in the ocean is closely related to stratification, and seasonal variation in the formation of a front and its effect on energy dissipation in the observation area was obvious.
To confirm the seasonal variation in the front in DYB, we analyzed observations of the area (114.68 • E, 22.58 • N) from July 2019 (Figure 1). At the observation site, 138 profiles were obtained by YODA within 48 h, with a time interval of 20 min. The temperature difference between the surface and bottom was more than 6 • C, and the salinity difference was more than 2 psu (Figure 13), which were much larger than the values measured in October ( Figure 10). The fresh water input during the summer rainy season diminished the surface salinity, and water with a lower temperature and high salinity intruded into the bay along the bottom from the SCS. Unlike the observations made in October, there was a stable bottom mixed layer in July that was more than 10 m thick (Figure 13d). Above the bottom mixing layer was a thermocline reflecting strong stratification. In the July 2019 observations, there was no obvious density front, but a high-frequency internal wave at the thermocline was observed. As shown in Figure 13c, the density (in the region of 1022.5 kg/m 2 ) oscillated within a thickness of 4 m in the thermocline. According to observations from October 2015, vertical fluctuation of the isotherm exceeded 10 m, and the density changed substantially at different times and at certain depths. observations, there was no obvious density front, but a high-frequency internal wave at the thermocline was observed. As shown in Figure 13c, the density (in the region of 1022.5 kg/m 2 ) oscillated within a thickness of 4 m in the thermocline. According to observations from October 2015, vertical fluctuation of the isotherm exceeded 10 m, and the density changed substantially at different times and at certain depths. In autumn, the upwelling along the coast of Guangdong disappeared due to the transition of the monsoon. Low-temperature and high-salinity water could not intrude into the bay along the bottom from the SCS, leading to the temperature and salinity of the bottom water being warmer and lower, respectively, than in summer (Figures 9 and 13). The larger temperature and salinity differences between the surface and bottom in summer led to stronger stratification in July 2019 compared to October 2015. As shown in Figure 14, the buoyancy frequency in summer in the thermocline was more than one order of magnitude higher. The strong stratification prevented the upper dynamic waves from crossing the thermocline, and the front was absent in summer. Thus, it was concluded that weak stratification in autumn is important for the formation of a front. In the final analysis, the transition of the monsoon was the primary factor causing changes in stratification. In autumn, the upwelling along the coast of Guangdong disappeared due to the transition of the monsoon. Low-temperature and high-salinity water could not intrude into the bay along the bottom from the SCS, leading to the temperature and salinity of the bottom water being warmer and lower, respectively, than in summer (Figures 9 and 13). The larger temperature and salinity differences between the surface and bottom in summer led to stronger stratification in July 2019 compared to October 2015. As shown in Figure 14, the buoyancy frequency in summer in the thermocline was more than one order of magnitude higher. The strong stratification prevented the upper dynamic waves from crossing the thermocline, and the front was absent in summer. Thus, it was concluded that weak stratification in autumn is important for the formation of a front. In the final analysis, the transition of the monsoon was the primary factor causing changes in stratification.

Conclusions
We used a tow-yo CTD instrument, the YODA profiler, to record high-resolution physical and biological data in DYB from a small vessel in October 2015. Observations in the DYB revealed complex coastal phenomena, which included warm-water discharge from a nuclear plant, as well as an upwelling, and strong density front. Those coastal phenomena changed the distribution of ecological elements. The upwelling resulted in saltier and more turbid water near the bottom, with low chlorophyll-a and DO contents, being transported upward to the surface layer and changing the local water environment. The front had a width of 1 km and formed at the intersection of the SCS and DYB. It was strengthened by the upwelling system.
Turbulence over a 22-h period at the mouth of DYB was also observed by using a microstructure turbulence profiler. The water column was divided into a surface mixed layer, thermocline, and bottom mixed layer. The turbulence was focused within the thermocline. The dissipation rate and diapycnal diffusivity in the thermocline displayed clear variability. The average (in time and space) dissipation rate and diapycnal diffusivity in the thermocline were 8 × 10 −8 W/kg and 5.8 × 10 −5 m 2 /s, respectively, which were one order of magnitude larger than the respective levels in an open-ocean thermocline. Active density fronts were observed in the thermocline and were associated with strong turbulent mixing. The active front was induced by weak stratification during the observation period. The results indicate that, rather than internal waves, it was the front that had formed at the intersection of different water bodies that played the primary role in energy dissipation at the mouth of DYB.

Conclusions
We used a tow-yo CTD instrument, the YODA profiler, to record high-resolution physical and biological data in DYB from a small vessel in October 2015. Observations in the DYB revealed complex coastal phenomena, which included warm-water discharge from a nuclear plant, as well as an upwelling, and strong density front. Those coastal phenomena changed the distribution of ecological elements. The upwelling resulted in saltier and more turbid water near the bottom, with low chlorophyll-a and DO contents, being transported upward to the surface layer and changing the local water environment. The front had a width of 1 km and formed at the intersection of the SCS and DYB. It was strengthened by the upwelling system.
Turbulence over a 22-h period at the mouth of DYB was also observed by using a microstructure turbulence profiler. The water column was divided into a surface mixed layer, thermocline, and bottom mixed layer. The turbulence was focused within the thermocline. The dissipation rate and diapycnal diffusivity in the thermocline displayed clear variability. The average (in time and space) dissipation rate and diapycnal diffusivity in the thermocline were 8 × 10 −8 W/kg and 5.8 × 10 −5 m 2 /s, respectively, which were one order of magnitude larger than the respective levels in an open-ocean thermocline. Active density fronts were observed in the thermocline and were associated with strong turbulent mixing. The active front was induced by weak stratification during the observation period.
The results indicate that, rather than internal waves, it was the front that had formed at the intersection of different water bodies that played the primary role in energy dissipation at the mouth of DYB.