Next Article in Journal
DETER-R: An Operational Near-Real Time Tropical Forest Disturbance Warning System Based on Sentinel-1 Time Series Analysis
Next Article in Special Issue
Radiation Sensitivity Analysis of Ocean Wake Information Detection System Based on Visible Light Remote Sensing
Previous Article in Journal
Analyses of the Dust Storm Sources, Affected Areas, and Moving Paths in Mongolia and China in Early Spring
Previous Article in Special Issue
Bayesian Sea Ice Detection Algorithm for CFOSAT
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Characteristics and Mixing Effects of Internal Solitary Waves in the Northern Yellow Sea as Revealed by Satellite and In Situ Observations

1
School of Marine Science and Technology, Tianjin University, Tianjin 300072, China
2
School of Mechanical Engineering, Tianjin University, Tianjin 300350, China
3
College of Marine and Environmental Sciences, Tianjin University of Science and Technology, Tianjin 300457, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(15), 3660; https://doi.org/10.3390/rs14153660
Submission received: 28 June 2022 / Revised: 20 July 2022 / Accepted: 26 July 2022 / Published: 30 July 2022
(This article belongs to the Special Issue Remote Sensing Applications in Ocean Observation)

Abstract

:
This study examines the characteristics, statistics, and mixing effects of internal solitary waves (ISWs) observed in the northern Yellow Sea (YS) during the summers of 2018 and 2019. The mooring stations are located between offshore islands with rough topographic features. Throughout the observation period, the ISWs with vertical displacements of up to 10 m induced prevailing high-frequency (3–10 min period) temperature variations. Synthetic aperture radar (SAR) images showed that the observed ISWs propagate in zonal directions generated around the islands where internal-tide-generating body force is strong. The estimated ISW propagation speed ranges from 0.16 to 0.25 m s−1, which agrees with the Korteweg-de Vries (KdV) model. The ISW intensity exhibits a clear spring-neap cycle corresponding to the local tidal forcing. The constant occurrence of ISWs at low tide suggests an important generation site where the ISWs are tidally generated. The ray-tracing result indicates that this generation site appears to be located at a strait between Dahao and Xiaohao islands. A generalized KdV model successively reproduces the propagation process from the generation site to the mooring station. Following the passage of ISWs, microstructure profiling observations reveal a high turbulent kinetic energy dissipation rate (10−6 W kg−1). The prevalence of ISWs in the study area is believed to play a crucial role in regulating vertical heat and nutrient transport, thereby modulating the biogeochemical cycle.

1. Introduction

Internal solitary waves (ISWs) are ubiquitous in world coastal oceans; they are long-lived features and can travel tens/hundreds of kilometers before dissipating [1,2]. These waves can generate strong reverse flow across the pycnocline, which induces the intense turbulent vertical exchange of heat and nutrients, regulating the biogeochemical cycle [3,4,5]. Therefore, comprehending the characteristics, statistics, and mixing effects of ISWs is crucial to understanding how marine ecosystems behave.
ISWs can be produced by a variety of mechanisms, including the lee wave mechanism [6], a transformation of the internal tide [7,8], or resonance in the transcritical flow [9], depending on the local stratification, topographic features, and tidal amplitudes. The lee wave mechanism has traditionally been used to explain the generation of ISWs at an underwater sill or bank [6]. The isopycnal disturbance on the topography exhibits a large lee wave when the flow reaches a supercritical flow condition (the Froude number Fr = U/c0 > 1, where U is the tidal current and c0 is the first-mode linear internal wave phase speed). The lee wave moves upstream as the tidal flow slows, evolving into several ISWs. Therefore, locally generated ISWs are often associated with tides. Kozlov et al. found a similar distance between adjacent ISW packages in the White Sea that matched the M2 wavelength value based on synthetic aperture radar (SAR) images [10]. According to statistical findings from years of observations in the Strait of Georgia, the ISW packets tended to appear at the specific tidal phase with similar time intervals [11], which emphasized the importance of tidal forcing in generating ISWs in the ocean.
ISWs can cause significant vertical displacements across the pycnocline and exacerbate turbulence [12,13]. In the investigation of Oregon’s continental shelf, Moum et al. revealed the evolving nature of interfaces through microstructure observations. As ISWs propagate shoreward, the interfaces may become unstable and break, creating turbulent flow [14]. Lee et al. investigated the generation and longevity of nonlinear internal waves in the northern East China Sea [4], suggesting that shear instability was the generation mechanism for the observed turbulence. The turbulent mixing can modify the transport of particulates and nutrients, influencing the distribution and recruitment of various littoral larval species [15]. These dynamical and ecological consequences provide additional motivation to investigate the causes of internal turbulent mixing and their effects on ecology and aquaculture.
The Yellow Sea (YS) is a temperate shelf sea in the northwest Pacific Ocean with a mean depth of 44 m. It connects to the East China Sea in the south and the Bohai Sea in the north. Limited observations have revealed the existence of ISWs, wind-induced near-inertial waves, and semidiurnal internal tides there [12,16,17,18,19]. Hsu et al. examined the distribution of ISWs in the YS based on SAR images onboard the satellites ERS-1 and ERS-2 [16]. They found that the ISWs were widely distributed along the coast of the YS, with the islands west of the Korea Peninsula being the local hot spots for ISWs due to their topographic features. After their generation on the west coast of the Korean Peninsula, ISWs propagated into the YS. Their phase speed increased nonlinearly with the water depth, averaging a value of approximately 0.1–0.7 m/s [20]. A recent numerical study showed that the simulated internal tides in the YS have an excellent spatial consistency with the satellite-detected ISWs because both the internal tides and ISWs are of tidal origin [21]. Moreover, Liu et al. observed that ISWs with a vertical displacement of 4–5 m had induced intensified turbulence in the pycnocline by creating small-scale shear instabilities based on microstructure observations in the local shelf break of the southern YS [12].
Only limited studies have looked into the ISWs in the northern YS. Lin et al. used numerical simulations to explore that energetic internal tides were generated near the coast around the northern YS [22], with a maximum baroclinic energy flux of 45 W/m. Although the baroclinic energy flux is much lower in the northern YS due to the shallow bathymetry, there are still significant large values at the topographic changes, which indicate the local generation of ISWs and internal tides. Alpers et al. have confirmed the presence of ISW packets in the northern YS through satellite imagery [17]. These images show the characteristics of ISWs in the northern YS, which appear frequently in the summer and seldom have wavelengths longer than 1000 m due to the shallow water. However, it is still challenging to determine the variability, property, and mixing effects of the ISWs because satellite observations alone cannot shed light on the features of ISWs in the interior of oceans. In situ observations are important methods to reveal the characteristics and mixing effects of ISWs [23,24]. Although the combination of satellite observations and in situ observations has been widely used [25,26], the characteristics of ISWs and their effects on turbulent mixing are still unclear in the northern YS due to the lack of in situ observations.
In this study, we conducted long-term high-frequency temperature observation for two mooring stations in the northern YS. SAR images and mooring observations were combined to analyze the statistical characteristics of the ISWs. The effect of ISWs on turbulent mixing was also discussed. This paper is organized as follows. Section 2 describes the observation data and methods. In Section 3, the results of the in situ hydrographic observations and satellite SAR investigations are presented. The correlation between ISWs and tide, potential generation sites, and the influence of ISWs on turbulence are discussed in Section 4. Conclusions are summarized in Section 5.

2. Materials and Methods

2.1. Mooring Observation

The data for this study were collected at two mooring stations in the northern YS from 29 June to 18 August 2018 and 5 June to 18 August 2019 (Figure 1b). The mooring stations (Stns.) A4 (39.03°N 122.98°E) and B4 (39.04°N 122.96°E) are located between two offshore islands (Zhangzi and Haiyang islands) in water with a mean depth of about 40 m. The T-chain is used to measure temperature at the mooring stations, equipped with four temperature sensors (T, ONT18S) with a frequency of 1 Hz at Stn. A4. At Stn. B4, the T-chain is equipped with two Conductivity-Temperature-Depth instruments (CTD, RBR420), two Temperature-Depth sensors (TD, RBR duet), and two temperature sensors (T, ONT18S). Temperature measurements are valid from 1.5 m to 25 m above the bottom. Table 1 summarizes information on the corrected sensor depths and temporal resolution. Figure 1c depicts the instrument configuration diagram.
We calculated the vertical velocity (w) from vertically densely spaced temperature sensor data after converting the series T(z, t) to isotherm displacements η (z, t) [27]. This study selects a set of isotherms with constant mean spacing (0.5 °C). For each isotherm, w is calculated by dividing the displacement of the isotherms in adjacent time by the time interval
w = η t
where η is the isotherm displacements in adjacent time and t is the time interval. The velocity of each layer is acquired by linear interpolation.

2.2. VMP Observations

We deployed the Vertical Microstructure Profiler- (VMP-) at Stn. B4 on 14–15 July 2019, to investigate the role of ISWs in inducing turbulent mixing. The VMP was activated once every hour for 22 h. The VMP is outfitted with two high-resolution frequency shear probes, one high-frequency temperature probe (FP07), and one high-frequency conductivity probe (SBE7), all with sampling rates of 512 Hz. The turbulent kinetic energy dissipation rate was calculated using microstructure shear. Using the assumption of isotropic turbulence, the turbulent kinetic energy dissipation rate ε was calculated as follows:
ε = 15 2 v ( u z ) 2 ¯ = 15 2 v k 1 k 2 ϕ ( k ) d z
where v is the kinematic molecular viscosity; ( u / z ) 2 ¯ is the variance of the vertical shear; φ ( k ) is the power spectrum of velocity shear; and k1 and k2 are the lower and upper limits of the wavenumber for integration, respectively. The calculation of ε is based on the ODAS Matlab Library Manual program by Rockland Scientific International (RSI) Inc. Previous studies have extensively described the corresponding data processing theories and methods (Gregg, 1999). For details of the calculation methods, please refer to Yang et al. [28] and Xu et al. [18]. The shear spectrum is calculated over consecutive segments of 2 m with a 50% overlap. As a result, we could obtain vertical profiles of ε with a vertical resolution of 1 m.

2.3. SAR Images

The ISW-induced currents modify sea surface roughness, which can be imaged using satellite images via Bragg backscattering [29,30,31]. Synthetic aperture radar (SAR) has long been a principal sensor in the observation of ISWs, because of its all-day, all-weather imaging capability [32]. SAR maps the sea surface roughness through Bragg scattering from the capillary waves and short gravity waves. The relation is
Λ = λ 2 sin θ
where Λ and λ are the wavelength of capillary waves and radar, respectively, and θ is the local incidence angle. In satellite SAR images, the ISWs typically appear as a pair of bright and dark stripes, corresponding to the rougher (convergence) and smoother (divergence) surface zones, respectively [30,32]. In this study, we used the C-band SAR of Sentinel-1 and GF-3 satellites launched by the European Space Agency on 3 April 2014 and China Aerospace Science and Technology Corporation on 10 August 2016, respectively. Both satellites can continuously image ISWs at 10 m resolution using HH-polarization. There are 59 images taken in the study area during the observation periods in 2018 and 2019, 11 of which show clear signatures of ISW packets.

2.4. The KdV Equation

The evolution of ISWs is generally influenced by topography and horizontally variable stratification [33,34]. For shallow water approximation, the KdV equation and its generalizations are widely used for the investigation of the propagation and transformation of ISWs in the shelf and coastal oceans [10,33]. To overcome the lack of direct observations of continuous stratification in the northern YS, we used the generalized Korteweg-de Vries (gKdV) equation to simulate the propagation of ISWs in the varied topography [10,35] as follows:
A t + C 0 ( x ) A x + α ( x ) A A x + α 1 ( x ) A 2 A x + β ( x ) A x x x + ( C 0 ( x ) 2 γ ( x ) ) γ x A ζ A x x = 0
where A(x, t) is the interface amplitude; C0(x) is the linear wave speed; and parameters α ( x ) , α 1 ( x ) , β ( x ) , γ ( x ) , and ζ are the coefficients for the nonlinear, higher-order nonlinear (cubic), dispersion, transformation, and dissipation effects, respectively.
The numerical model based on the gKdV equation was discretized by the Predictor-Evaluation-Corrector-Evaluation (PECE) method and iterated by the implicit corrector scheme [36], which ensures convergence and computational efficiency of the iterative process. There is a steady summer pycnocline between the surface and intermediate water masses along the ISW propagation trajectory at depths of 5–10 m. Segur and Hammack [37] and Koop and Butler [38] found that the KdV model can be simplified to a two-layer model and predicted the solitary waves with remarkable accuracy. Thus, a two-layer system of upper quasi-homogeneous and lower layers with thicknesses h1 and h2, respectively, can be used [10,33]:
C 0 = Δ ρ g h 1 h 2 / ρ ( h 1 + h 2 )
α = [ 3 ( h 1 h 2 ) / 2 ( h 1 h 2 ) ] C 0
α 1 = [ 3 ( h 1 2 + h 2 2 + 6 h 1 h 2 ) / 8 ( h 1 h 2 ) 2 ] C 0
β = ( h 1 h 2 / 6 ) C 0
γ = { C 0 ( 0 ) C 0 1 [ h 1 1 ( 0 ) + h 2 1 ( 0 ) ] / [ h 1 1 + h 2 1 ] } 3
where g is the acceleration of gravity; ρ 1 , h1 and ρ 2 , h2 are the density and thickness of upper and lower layers, respectively.

3. Results

This section is divided into subheadings. It should provide a concise and precise description of the experimental results, their interpretation, as well as the experimental conclusions that can be drawn.

3.1. Temperature Variations

Figure 2 and Figure 3 show the temperature variation observed by the T-chain at Stns. A4 and B4, respectively. These two stations shared many characteristics due to their close distance and similar time of the year for observation. The observed temperature shows a warming trend at both stations. The temperature measured by the lowest sensors increased from 10 to 18 °C and from 8 to 20 °C at Stns. A4 and B4, respectively. The observed part of the water column is usually stratified, with the largest vertical temperature difference reaching 11.2 °C and 12.8 °C at Stns. A4 and B4, respectively. There was a significant increase in temperature throughout the water column on 13 August at Stn. B4 (Figure 3a,b). The bottom temperature changed from 15 to 19 °C, after which the temperature fluctuated and dropped gradually. The surface elevation demonstrates that semidiurnal tides dominate at both stations with a distinct spring-neap cycle (Figure 2a and Figure 3a).
One notable phenomenon observed at both stations is prevailing temperature spikes that last throughout the observation period. Close-view examples of the temperature spikes show that these spikes have frequencies ranging from 3–10 min (Figure 2c,e and Figure 3c,e). When the bottom layer is well-mixed, temperature spikes may appear only in the upper sensors (Figure 2d,f). The uppermost temperature sensor increased from 14 to 23 °C in 100 s with a period of 5 min (Figure 2c). Figure 2c,e and Figure 3c,e show that these temperature spikes are caused by the prevailing high-frequency ISWs of depression in this area. The vertical displacement of the ISWs is about 8–10 m (Figure 2c). Additionally, the leading waves have the largest vertical displacement of 13 m, corresponding to the highest temperature variation of 7 °C, as observed in Figure 3d,f. The observed prevailing ISW packets in the study area are surprising. This is because, although the previous SAR imagery studies have revealed the presence of ISW packets around offshore islands in the northern YS, this region has never been considered a hot spot for ISWs [17]. Given the prevalence of the observed ISWs, they may have significant biogeochemical implications. Next, the properties of the ISWs are examined.

3.2. Spectral Estimates

The frequency spectrum of temperature is calculated to examine the energy content as a function of frequency. The temperature in the uppermost layer is selected to calculate the frequency spectrum. Figure 4c,d depicts the frequency spectrum calculated from mooring temperature observations at Stns. A4 and B4, respectively. The frequency spectra are dominated by near-inertial/diurnal (~1 cycle per day (cpd)), semidiurnal (~2 cpd), and higher tidal harmonics at both stations. The strongest peak corresponds to the frequency of the M2 tide. The spectra generally decrease with increasing frequency, except where there seems to exist a spectral bump at ~144 to 480 cpd (a period of ~3 to 10 min). Figure 2 and Figure 3 depict that the spectral bump is associated with the high-frequency temperature variations.
The high-frequency region of the spectra was then examined to obtain a better look at the spectral bump. In order to examine the temporal variation of the high-frequency motions, we calculated the frequency spectrum for each day by dividing the temperature variation into daily segments. Figure 4a,b shows that high spectral energy appears between 3 and 10 min periods, which corresponds to the spectral bump. The spectral bump almost lasts throughout the observation period. This is consistent with the observed frequent temperature spikes, indicating that ISWs are common in the study area.

3.3. Close View of the ISWs

The SAR images and mooring observations are combined to examine the properties of ISWs. For example, we selected two typical SAR images (10 m resolution, HH-polarization) with clear ISWs acquired on 11 July and 8 August 2019. The in situ mooring observations are analyzed alongside the SAR images to provide additional information about propagation details. Figure 5d depicts the capture of two ISW packets (ISW A and ISW B) at 3.4 and 7.9 km from Stn. B4. ISW A arrives at Stn. B4 approximately 4 h after the imaging time, causing significant temperature variations. Within 3 min, the temperature measured by the uppermost sensor rises from 13 to 22 °C. Close inspection reveals that ISW A has a vertical displacement of 10 m and a period of about 10 min (Figure 5c). The ISW propagation speed is estimated to be 0.25 m s−1 based on relative packet positions and arrival times. ISW B originates on Haiyang Island and has a propagation speed of 0.16 m s−1.
The other SAR image was acquired on 8 August 2019, at 17:46, capturing at least four ISW packets (Figure 6d). Three of the four packets propagate from west to east in a similar direction. ISW C first passes through the mooring station, causing high-frequency temperature variations for 28 min. This ISW has an amplitude of approximately 9 m and a period of 7 min (Figure 6c). The propagation speed of ISW C is estimated to be 0.21 m s−1. The propagation speeds of ISWs A and B are estimated to be 0.21 and 0.17 m s−1, respectively, based on the distance revealed by SAR images and arrival time. Wavelengths appear to be variable during the propagation of the ISW. ISW A has the largest wavelength at about 500 m, and ISWs B and C have relatively small wavelengths of less than 200 m. In the same ISW packet, the wavelengths of ISWs are also different. As shown in ISW B and ISW C in Figure 6, the wavelengths are significantly larger in the southern part of the waves.
The linear phase speed at any frequency depends on the vertical density profile of the water column, based on the KdV equation for internal waves [39]. The two-layer model, consisting of the upper layer with thickness h1 and density ρ 1 , and the lower layer with thickness h2 = Hh1 and density ρ 2 > ρ 1 , is widely considered for the simulation of internal waves [40,41]. In this case, assuming the ISWs are of the KdV type, the speeds of propagation of ISWs are generally calculated from Equation (10):
c = g ( ρ 2 ρ 1 ) h 1 h 2 ρ 2 h 1 + ρ 1 h 2 + α η 0 3
where g is the acceleration due to gravity; η 0 is the vertical displacement of ISWs; ρ 1 and ρ 2 are the upper layer and lower layer density; h1 and h2 are the upper layer and lower layer thickness, respectively. The parameter α is calculated using Equation (5). At the mooring station, the phase speed calculated according to Equation (10) is about 0.25 m s−1, consistent with estimates from satellite observations.

4. Discussion

4.1. Statistics of ISWs

Mooring observations indicate that ISWs are prevailing phenomena in our study area that exist throughout the observation period in summer. We performed a statistical analysis on the time of appearance of ISWs to the tidal phases to investigate their characteristics and generation mechanisms. Wave packets were defined as ISWs with vertical displacements greater than 5 m. At Stns. A4 and B4, a total of 125 and 277 ISW packets were observed, respectively. The surface elevation shows a clear spring-neap cycle at both stations (Figure 2a and Figure 3a). The occurrence of ISWs is binned by their time relative to the nearest maximum of the fortnightly cycle of tidal magnitudes and by their time relative to the low tide. At Stn. A4, more ISWs seem to be observed during spring tide than during neap tide (Figure 7c). However, ISWs do not appear more frequently at Stn. B4 during spring tide (Figure 7f), showing no obvious correlation with tidal amplitude. The one-hour and one-day root-mean-square (RMS) vertical velocity ( σ w ) are then calculated (Figure 7b,e). The ( σ w ) shows a clear fortnightly cycle which agrees well with the variation of RMS tidal amplitude ( σ H ) .
Variations in both the occurrence frequency and intensity of the ISWs can affect the ( σ w ) . Previous research found that the likelihood of ISW packets peaked around the strongest tides in the spring-neap tide [11,17,42]. However, the likelihood of ISWs shows no apparent correlation with the spring-neap tide in the present study (Figure 7c,f). Thus, we suggest that the fortnightly variations of ( σ w ) are responsible for modulating the ISW intensity by the spring-neap tide.
The occurrence of the ISW packet was examined by comparing it with the flood-ebb cycle of the local tides. At both stations, we first calculated the average value of the surface elevation for 25 h, showing a clear dominant semidiurnal cycle (Figure 8). The occurrence of ISWs was classified according to the local tidal phases. Figure 8 shows that most ISWs are observed at both stations during low tides. This relationship appears to be stronger at Stn. B4. In other regions, the occurrence of ISWs at a particular tidal phase has been interpreted as evidence that their generation mechanisms are connected to tidal phases [11]. This strengthens the case that tides cause the observed ISWs.
Figure 9 depicts a typical example of the phase-locked ISW packets at Stn. B4. ISW packets always occurred approximately one hour later than the previous day (marked by red arrows). Although the ISW packets appear at irregular times of the day, they occur at every low-tide phase. Furthermore, the intensity of ISWs occurring during spring tides (19–22 July) tends to be greater than during neap tides (12–18 July), consistent with the statistical results. The occurrence of ISWs at the fixed low-tide phase suggests that tide-topography interactions generate ISWs.

4.2. Generation Sites of ISWs

Next, we investigated the generation and propagation of ISWs. Figure 10a depicts an ISW distribution map of the northern YS based on the SAR images collected during the observation period. ISWs propagate in various directions, implying multiple generation sites. Specifically, the observed ISWs at the mooring station are generated at the rough topography around the island rather than being locally generated. Many ISWs propagating eastwards seem to originate from Zhangzi Island. A small number of westward-propagating ISWs also appear, which seem to have been generated around Haiyang Island.
In order to investigate the generation sites of the ISWs, we examined the distribution of internal-tide-generating body force in the northern YS (Figure 10). The depth-integrated internal-tide-generating body force is calculated as follows (Baines, 1982):
F = 1 ω 1 H 2 H 0 N 2 ( z ) z d z ( Q x d H d x ) 2 + ( Q y d H dy ) 2
where ω is the M2 tidal frequency; H is the water depth; N2(z) is the squared buoyancy frequency; Qx and Qy are the zonal and meridional components of the tidal transport calculated by (Qx, Qy) = (uH, vH); u and v are the zonal and meridional components of the barotropic tidal current; and d H d x and d H d y are the zonal and meridional bottom slopes. The bottom topography is based on the 0.25′ grid general bathymetric chart of the oceans (GEBCO_2019) provided by the Intergovernmental Oceanographic Commission (IOC) and the International Hydrographic Organization (IHO). The OSU Tidal Prediction Software (OTPS) forecasts the zonal and meridional components of tidal velocity. Potential hot spots of internal wave generation sites around the world ocean have been examined using the spatial distribution of the F [43,44,45]. F > 0.25 m2 s−2 is usually regarded as the critical value for the generation of ISWs [44,46]. Figure 10a depicts that large F tends to be located around islands with striking topographic features.
That energetic internal tides were generated close to the offshore islands around the northern YS is consistent with recent numerical simulations [22]. Both SAR images and calculated F indicate that the ISWs are generated at these offshore islands and then propagate to the mooring stations.
Statistical analysis showed that the ISWs are related to tidal phases, with most ISWs occurring during low tide (Figure 8 and Figure 9). As a result, even though the SAR images show multiple ISW generation sites around the islands, there should be a major ISW generation site with tidally-generated ISWs reaching the mooring station periodically. Figure 9 depicts an example of ISWs fixed to local low tides, implying that they are generated tidally at a fixed location. The data from SAR images are then combined with in situ observations and the KdV model to provide additional insights into the generation and propagation details.
The ray-tracing method was used to find the major generation sites of ISWs in the northern YS [47,48,49]. Assuming that the ISWs propagate with a circular wave crest, the generation site of ISWs is defined as the intersection of two perpendicular lines across the wave crest. Figure 11a depicts an ISW packet captured on 11 July 2019, corresponding to the first ISW packet that appeared at low tide in Figure 9. The observed ISW packet from satellite images traces the generation site S1 to the straits between Dahao and Xiaohao islands, approximately 10.1 km from Stn. B4 (Figure 10b and Figure 11a). The generalized KdV model was then used to simulate the propagation of ISWs. The density is referred from the WOA18 database, showing that the upper and lower layer densities in the study area are approximately 1023 and 1024 kg m−3, respectively. We then simulated the propagation of ISWs with the generalized KdV model. The real topography from GEBCO was used for simulation. The density was referred to in the WOA18 database, showing that the upper and lower layer densities are about 1023 and 1024 kg m−3 in the study area. The depth of pycnocline is 5 m. After that, Equation (4) is numerically solved on a regular grid with Δ x = 10 m, Δ t = 10 s, the viscosity coefficient ζ = 0.0025 m2 s−2 [10], and a total integration time of 20 h. At x = 0 km, a train of two consecutive waves with vertical displacements of −5 and −4 m appears.
Figure 11b shows that the ISW packet undergoes several changes as they propagate seaward from S1. Their vertical displacements are significantly reduced, and they become broader and flatter. The ISW packet reaches the position corresponding to the SAR image after ~8.2 h of propagation from S1 (Figure 11a). The propagation from S1 to B4 takes about 12 h in total, with a mean propagating speed of ~0.25 m s−1, consistent with the velocity estimated from the SAR image and mooring observations (Figure 5d). Figure 11c shows the barotropic tide at S1 when ISWs are generated. At low tide, when the direction of tidal velocities is toward the northeast, ISWs are generated in the straits between the islands. Although the site between the islands was estimated to be a critical ISW generation site, ISWs coming from different directions and at different times indicate the existence of other source regions. Furthermore, although we have shown that the ISWs are tidally generated at rough topography, the exact generation mechanism of the ISWs also awaits further investigation in the future based on numerical simulations. Nevertheless, the prevalence and properties of ISWs originating from offshore islands in the northern YS are documented for the first time based on comprehensive in situ observations and SAR images.

4.3. Influences of ISWs on Internal Turbulence

Turbulent mixing in shelf seas is a critical process that controls diapycnal nutrient transport, and it plays a vital role in biological production and carbon cycles [50,51]. The one-day VMP observations enable us to investigate the impact of ISWs on internal turbulence. Figure 12b shows the time-depth evolution of the observed ε and temperature. Previous studies showed that the energy input, conversion, and radiation rates increase monotonically with ε [52]. The largest ε (~10−6 W kg−1) occurred in the well-mixed bottom boundary layer (BBL). One of the most noticeable features here is the presence of two peaks near the pycnocline. At the 4th hour after 13:00 on 14 July 2019, the ε reached the first peak with a maximum value of 10−7 W kg−1. About eleven hours later, the ε reached the other peak with a larger value of 10−6 W kg−1.
It is widely believed that turbulence in the surface boundary layer (SBL) in shelf seas can be strongly influenced by wind forcing. In order to investigate the impact of wind on the two ε peaks, we investigated the variation of wind speed at 10 m above the sea surface from the fifth generation European Centre for Medium-Range Weather Forecasts (ECMWF)’s atmospheric reanalysis of the global climate (ERA5). Figure 12a presents that the wind speed is less than 5 m s−1 during the observation period, suggesting that the wind may not be the cause of turbulence. Surface buoyancy fluxes do not dramatically change when ε is intensified, which shows that changes in surface buoyancy fluxes do not cause internal turbulence. Kay and Jay found that substantial internal turbulence occurred even during periods of weak near-bed shear [53], suggesting that internal processes must contribute significantly to the internal turbulence. Figure 12c depicts the variations in internal temperature and turbulence. The temperature variations show a semidiurnal period, consistent with the M2 tide. Around four hours (17:00 14 July 2019) and thirteen hours (02:00 15 July 2019) after the beginning of the observations, the temperature shows continuous high-frequency variations. Furthermore, we selected the depths of the σ θ = 22.5 and 20.9 kg m−3 isopycnals to represent the lower and upper boundaries of the pycnocline, respectively (as shown by the white contours in Figure 12b). The two ε peaks in the pycnocline correspond to the high-frequency temperature variations (Figure 12c). The average ε in the pycnocline also shows a similar trend as the one-day RMS vertical velocity (Figure 12d). The largest root-mean-square vertical velocity appears at 15 h, corresponding to the long-lasting strong turbulence from 15 to 17 h. Although their largest peaks do not exactly coincide with each other, we suggest that the sustaining strong turbulence is related to the high-frequency temperature variation event which has induced a long-lasting influence. We therefore suggest that ISWs cause the observed internal turbulence.
The present observations demonstrate that ISWs are quite prevalent in the study area, originating at rough topographic features near the offshore island. The one-day VMP observation shows how passing ISWs can significantly alter internal turbulence in the northern YS. Moreover, the frequent turbulent mixing induced by ISWs can increase vertical nutrient supply and influence the local biogeochemical cycle [54,55,56]. Furthermore, the areas around the Changshan Islands are crucial aquaculture regions in China [57,58]. The presence of ISWs, which cause increased vertical flux and high-frequency temperature variations during warm seasons, may have significant ecological consequences that must be investigated.

5. Conclusions

The properties, generation mechanisms, statistics, and mixing effects of ISWs are investigated using mooring observations and SAR images from the northern YS. According to the spectral analysis, previous ISWs induced numerous high-frequency temperature spikes with periods ranging from 3 to 10 min. These ISWs are mainly depression type, with vertical displacements of up to 10 m. Dramatic vertical displacements can lead to nearly 9 °C temperature variations within the pycnocline. We estimated the propagation velocity of ISWs to be 0.16 to 0.25 m s−1 by combining mooring observations and SAR images, which is consistent with the theoretical value. According to the statistical analysis, ISWs at the mooring stations do not indicate a clear preference for spring tides, but the RMS vertical velocity of ISWs (ISW intensity) exhibits a spring-neap tidal cycle and is linearly proportional to the barotropic tidal height. Furthermore, while ISWs can be found at the mooring stations at any tidal phase, they are far more prevalent during low tides.
These findings indicate that the ISWs are tidal in origin, with a major generation site. SAR images obtained during the observation period show that most ISWs in this region propagate in the zonal direction, which seems to originate from the Zhangzi and Haiyang islands. This is supported by the horizontal distribution of the depth-integrated internal-tide-generating body force, which shows consistent features of high values around the offshore islands. The ray-tracing method was used to locate the major generation site, which appeared to be located in the straits surrounding Zhangzi Island. The KdV model was used to simulate the propagation of the ISWs from the generation site to the mooring station. The ISWs take ~12 h to propagate from the generation site to the mooring station, which coincides with the observational estimates.
The ε showed two large values (~10−6 W kg−1) when there were temperatures with high-frequency variations caused by ISW packets. The similar variation trends between ε and vertical velocity confirmed that ISWs cause internal turbulence. Turbulent mixing induced by frequent ISWs can increase the vertical nutrient supply and influence the biogeochemical cycle. Considering that the regions around these offshore islands are now important bottom-aquaculture farms in China, the influence of the prevailing ISWs on aquaculture remains to be examined.

Author Contributions

Conceptualization, H.L. and W.Y.; methodology, H.L.; software, C.J.; validation, H.L., W.Y., H.W., C.L., L.Z. and C.J.; formal analysis, H.L.; investigation, H.L.; resources, W.Y.; data curation, H.L. and W.Y.; writing—original draft preparation, H.L.; writing—review and editing, W.Y., H.W., L.Z. and C.L.; visualization, H.L., W.Y., H.W., L.Z. and C.J.; supervision, W.Y.; project administration, W.Y.; funding acquisition, W.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the National Natural Science Foundation of China (41906006), the Tianjin Natural Science Foundation (21JCYBJC00500), and the National Key Research and Development Program of China (2018YFD0900704).

Data Availability Statement

The SAR images from Sentinel-1 and GF-3 satellites were obtained online (https://scihub.copernicus.eu/dhus/#/home (accessed on 1 December 2020) and https://osdds.nsoas.org.cn/#/ (accessed on 1 December 2020)). The ERA5 data were obtained online (https://cds.climate.copernicus.eu/cdsapp#!/home, accessed on 1 May 2021). Mooring and VMP data analyzed in this study are available at https://zenodo.org/record/6609695 (accessed on 30 May 2022).

Acknowledgments

The acquisition of mooring data was supported by Laboratory of Ocean Dynamics and Ecology, LODE. We are grateful to the entire crew for their expertise and hard work. The authors thank the editor and three anonymous reviewers for their constructive suggestions to improve the quality of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Helfrich, K.R.; Melville, W.K. Long nonlinear internal waves. Annu. Rev. Fluid Mech. 2006, 38, 395–425. [Google Scholar] [CrossRef]
  2. Jackson, C. Internal wave detection using the Moderate Resolution Imaging Spectroradiometer (MODIS). J. Geophys. Res. Atmos. 2007, 112, C11012. [Google Scholar] [CrossRef] [Green Version]
  3. Duda, T.F.; Lynch, J.F.; Irish, J.D.; Beardsley, R.C.; Ramp, S.R.; Chiu, C.S.; Tang, T.Y.; Yang, Y.J. Internal tide and nonlinear internal wave behavior at the continental slope in the northern south China Sea. IEEE J. Ocean. Eng. 2005, 29, 1105–1130. [Google Scholar] [CrossRef] [Green Version]
  4. Lee, J.H.; Lozovatsky, I.; Jang, S.T.; Chan, J.J.; Fernando, H. Episodes of nonlinear internal wave in the northern East China Sea. Geophys. Res. Lett. 2006, 33. [Google Scholar] [CrossRef]
  5. Dong, J.; Zhao, W.; Chen, H.; Meng, Z.; Shi, X.; Tian, J. Asymmetry of internal waves and its effects on the ecological environment observed in the northern South China Sea. Deep-Sea Res. Part I 2015, 98, 94–101. [Google Scholar] [CrossRef]
  6. Maxworthy, T. A note on the internal solitary waves produced by tidal flow over a three-dimensional ridge. J. Geophys. Res. Oceans 1979, 84, 338–346. [Google Scholar] [CrossRef]
  7. Holloway, P.E. Internal hydraulic jumps and solitons at a shelf break region on the Australian North West Shelf. J. Geophys. Res. Atmos. 1987, 92, 5405–5416. [Google Scholar] [CrossRef]
  8. Sandstrom, H.; Elliot, J.A.; Cchrane, N.A. Observing groups of solitary internal waves and turbulence with BATFISH and echo-sounder. J. Phys. Oceanogr. 1989, 19, 987–997. [Google Scholar] [CrossRef] [Green Version]
  9. Melville, W.K.; Helfrich, K.R. Transcritical two-layer flow over topography. J. Fluid Mech. 1987, 178, 31–52. [Google Scholar] [CrossRef] [Green Version]
  10. Kozlov, I.; Romanenkov, D.; Zimin, A.; Chapron, B. SAR observing large-scale nonlinear internal waves in the White Sea. Remote Sens. Environ. 2014, 147, 99–107. [Google Scholar] [CrossRef] [Green Version]
  11. Li, L.; Pawlowicz, R.; Wang, C. Seasonal variability and generation mechanisms of nonlinear internal waves in the strait of Georgia. J. Geophys. Res. Oceans 2018, 123, 5706–5726. [Google Scholar] [CrossRef]
  12. Liu, Z.; Wei, H.; Lozovatsky, I.D.; Fernando, H.J.S. Late summer stratification, internal waves, and turbulence in the Yellow Sea. J. Mar. Syst. 2009, 77, 459–472. [Google Scholar] [CrossRef]
  13. Bourgault, D.; Janes, D.C.; Galbraith, P.S. Observations of a large-amplitude internal wave train and its reflection off a steep slope. J. Phys. Oceanogr. 2011, 41, 586–600. [Google Scholar] [CrossRef]
  14. Moum, J.N.; Farmer, D.M.; Smyth, W.D.; Armi, L.; Vagle, S. Structure and generation of turbulence at interfaces strained by internal solitary waves propagating shoreward over the continental shelf. J. Phys. Oceanogr. 2003, 33, 2093–2112. [Google Scholar] [CrossRef]
  15. Scotti, A.; Pineda, J. Obervation of very large and steep waves of elevation near the Massachusetts coast. Geophys. Res. Lett. 2004, 31. [Google Scholar] [CrossRef] [Green Version]
  16. Hsu, M.; Liu, A.K.; Liu, C. A study of internal waves in the China Seas and Yellow Sea using SAR. Cont. Shelf Res. 2000, 20, 389–410. [Google Scholar] [CrossRef]
  17. Alpers, W.; He, M.X.; Zeng, K.; Guo, L.F.; Li, X.M. The distribution of internal waves in the East China Sea and the Yellow Sea studied by multi-sensor satellite images. In Proceedings of the IEEE International Geoscience & Remote Sensing Symposium, Seoul, Korea, 29 July 2005. [Google Scholar] [CrossRef]
  18. Xu, P.; Yang, W.; Zhu, B.; Wei, H.; Zhao, L.; Nie, H. Turbulent mixing and vertical nitrate flux induced by the semidiurnal internal tides in the southern Yellow Sea. Cont. Shelf Res. 2020, 208, 104240. [Google Scholar] [CrossRef]
  19. Yang, W.; Wei, H.; Liu, Z.; Li, G. Intermittent intense thermocline shear associated with wind-forced near-inertial internal waves in a summer stratified temperate shelf sea. J. Geophys. Res. Oceans 2021, 126, e2021JC017576. [Google Scholar] [CrossRef]
  20. Kim, H.; Son, Y.B.; Jeong, J.; Jo, Y. Comparison of internal waves in various ocean fields around the Korean Peninsula. J. Coast. Res. 2018, 85, 466–470. [Google Scholar] [CrossRef]
  21. Liu, K.; Sun, J.; Guo, C.; Yang, Y.; Yu, W.; Wei, Z. Seasonal and spatial variations of the M2 internal tide in the Yellow Sea. J. Geophys. Res. Oceans 2019, 124, 1115–1138. [Google Scholar] [CrossRef]
  22. Lin, F.; Asplin, L.; Wei, H. Summertime M2 internal tides in the Northern Yellow Sea. Front. Mar. Sci. 2021, 8, 798504. [Google Scholar] [CrossRef]
  23. Lien, C.R. Energy of nonlinear internal waves in the South China Sea. Geophys. Res. Lett. 2005, 32, 215–236. [Google Scholar] [CrossRef] [Green Version]
  24. McSweeney, J.M.; Lerczak, J.A.; Barth, J.A.; Becherer, J.; Colosi, J.A.; MacKinnon, J.A.; MacMahan, J.H.; Moum, J.N.; Pierce, S.D.; Waterhouse, A.F. Observations of shoaling nonlinear internal bores across the Central California inner shelf. J. Phys. Oceanogr. 2020, 50, 111–132. [Google Scholar] [CrossRef]
  25. Lai, Z.; Jin, G.; Huang, Y.; Chen, H.; Shang, X.; Xiong, X. The generation of nonlinear internal waves in the South China Sea: A three-dimensional, nonhydrostatic numerical study. J. Geophys. Res. Oceans 2019, 124, 8949–8968. [Google Scholar] [CrossRef]
  26. Bai, X.; Liu, Z.; Zheng, Q.; Hu, J.; Lamb, K.G.; Cai, S. Fission of shoaling internal waves on the Northeastern Shelf of the South China Sea. J. Geophys. Res. Oceans 2019, 124, 4529–4545. [Google Scholar] [CrossRef]
  27. Pinkel, R. Observations of the near-surface internal wavefield. J. Phys. Oceanogr. 1981, 11, 1248–1257. [Google Scholar] [CrossRef] [Green Version]
  28. Yang, W.; Wei, H.; Zhao, L.; Zhang, J. Turbulence and vertical nitrate flux adjacent to the Changjiang Estuary during fall. J. Mar. Syst. 2020, 212, 103427. [Google Scholar] [CrossRef]
  29. Valenzuela, G.R. Theories for the interaction of electromagnetic and oceanic waves—A review. Bound.-Lay. Meteorol. 1978, 13, 61–85. [Google Scholar] [CrossRef]
  30. Alpers, W. Theory of radar imaging of internal waves. Nature 1985, 314, 245–247. [Google Scholar] [CrossRef]
  31. Chang, M.H.; Lien, R.C.; Yang, Y.J.; Tang, T.Y.; Wang, J. A composite view of surface signatures and interior properties of nonlinear internal waves: Observations and applications. J. Atmos. Ocean. Technol. 2008, 25, 523–531. [Google Scholar] [CrossRef]
  32. Fu, L.L.; Holt, B. Internal waves in the Gulf of California—Observations from a spaceborne radar. J. Geophys. Res. Oceans 1984, 89, 2053. [Google Scholar] [CrossRef]
  33. Liu, A.K.; Chang, Y.S.; Hsu, M.K.; Liang, N.K. Evolution of nonlinear internal waves in the East and South China Seas. J. Geophys. Res. Oceans 1998, 103, 7995–8008. [Google Scholar] [CrossRef]
  34. Zhao, Z. Satellite observation of internal solitary waves converting polarity. Geophys. Res. Lett. 2003, 30. [Google Scholar] [CrossRef]
  35. Grimshaw, R.; Pelinovsky, E.; Talipova, T. Modelling internal solitary waves in the Coastal Ocean. Surv. Geophys. 2007, 28, 273–298. [Google Scholar] [CrossRef]
  36. Djidjeli, K.; Price, W.G.; Twizell, E.H.; Wang, Y. Numerical methods for the solution of the third- and fifth-order dispersive Korteweg-de Vries equations. J. Comput. Appl. Math. 1995, 58, 307–336. [Google Scholar] [CrossRef] [Green Version]
  37. Segur, H.; Hammack, J.L. Soliton models of long internal waves. J. Fluid Mech. 1982, 118, 285–304. [Google Scholar] [CrossRef]
  38. Koop, C.G.; Butler, G. An investigation of internal solitary waves in a two-fluid system. J. Fluid Mech. 1981, 112, 225–251. [Google Scholar] [CrossRef]
  39. Apel, J.R.; Ostrovsky, L.A.; Stepanyants, Y.A.; Lynch, J.F. Internal solitons in the ocean and their effect on underwater sound. J. Acoust. Soc. Am. 2007, 121, 695–722. [Google Scholar] [CrossRef]
  40. Liu, B.; Yang, H.; Zhao, Z.; Li, X. Internal solitary wave propagation observed by tandem satellites. Geophys. Res. Lett. 2014, 41, 2077–2085. [Google Scholar] [CrossRef]
  41. Shen, H.; Perrie, W.; Johnson, C.L. Predicting internal solitary waves in the gulf of maine. J. Geophys. Res. Oceans 2020, 125, e2019JC015941. [Google Scholar] [CrossRef] [Green Version]
  42. Colosi, J.A.; Kumar, N.; Suanda, S.H.; Freismuth, T.M.; MacMahan, J.H. Statistics of internal tide bores and internal solitary waves observed on the inner continental shelf off Point Sal, California. J. Phys. Oceanogr. 2018, 48, 123–143. [Google Scholar] [CrossRef] [Green Version]
  43. Sherwin, T.J.; Vlasenko, V.I.; Stashchuk, N.; Jeans, D.R.G.; Jones, B. Along-slope generation as an explanation for some unusually large internal tides. Deep-Sea Res. Part I: Oceanogr. Res. Pap. 2002, 49, 1787–1799. [Google Scholar] [CrossRef]
  44. Li, X.; Zhao, Z.; Han, Z.; Xu, L. Internal solitary waves in the East China Sea. Acta Oceanol. Sin. 2008, 27, 51–59. [Google Scholar] [CrossRef]
  45. Bai, X.; Liu, Z.; Li, X.; Hu, J. Generation sites of internal solitary waves in the southern Taiwan Strait revealed by MODIS true-colour image observations. Int. J. Remote Sens. 2014, 35, 4086–4098. [Google Scholar] [CrossRef]
  46. Lozovatsky, I.; Liu, Z.; Fernando, H.; Armengol, J.; Roget, E. Shallow water tidal currents in close proximity to the seafloor and boundary-induced turbulence. Ocean Dyn. 2012, 62, 177–191. [Google Scholar] [CrossRef]
  47. Da Silva, J.C.B.; Buijsman, M.C.; Magalhaes, J.M. Internal waves on the upstream side of a large sill of the Mascarene Ridge: A comprehensive view of their generation mechanisms and evolution. Deep Sea Res. Part I Oceanogr. Res. Pap. 2015, 99, 87–104. [Google Scholar] [CrossRef] [Green Version]
  48. Sun, L.; Zhang, J.; Meng, J. A study of the spatial-temporal distribution and propagation characteristics of internal waves in the Andaman Sea using MODIS. Acta Oceanol. Sin. 2019, 38, 121–128. [Google Scholar] [CrossRef]
  49. Zhang, X.; Li, X.; Zhang, T. Characteristics and generations of internal wave in the Sulu Sea inferred from optical satellite images. J. Oceanol. Limnol. 2020, 38, 1435–1444. [Google Scholar] [CrossRef]
  50. Rippeth, T.P.; Palmer, M.R.; Simpson, J.H.; Fisher, N.R.; Sharples, J. Thermocline mixing in summer stratified continental shelf sea. Geophys. Res. Lett. 2005, 32, 159–171. [Google Scholar] [CrossRef]
  51. Sharples, J.; Moore, M.C.; Hickman, A.E.; Holligan, P.M.; Simpson, J.H. Internal tidal mixing as a control on continental margin ecosystems. Geophys. Res. Lett. 2009, 36, L23603. [Google Scholar] [CrossRef] [Green Version]
  52. Chen, Z.; Xie, J.; Xu, J.; Zhan, J.; Cai, S. Energetics of nonlinear internal waves generated by tidal flow over topography. Ocean Model. 2013, 68, 1–8. [Google Scholar] [CrossRef]
  53. Kay, D.J. Interfacial mixing in a highly stratified estuary 1. Characteristics of mixing. J. Geophys. Res. Oceans 2003, 108. [Google Scholar] [CrossRef] [Green Version]
  54. Williams, C.; Sharples, J.; Mahaffey, C.; Rippeth, T. Wind-driven nutrient pulses to the subsurface chlorophyll maximum in seasonally stratified shelf seas. Geophys. Res. Lett. 2013, 40, 5467–5472. [Google Scholar] [CrossRef] [Green Version]
  55. Shroyer, E.L.; Benoit-Bird, K.J.; Nash, J.D.; Moum, J.N. Stratification and mixing regimes in biological thin layers over the Mid-Atlantic Bight. Limnol. Oceanogr. 2014, 59, 1349–1363. [Google Scholar] [CrossRef] [Green Version]
  56. Sharples, J.; Zeldis, J.R. Variability of internal tide energy, mixing and nitrate fluxes in response to changes in stratification on the northeast New Zealand continental shelf. N. Z. J. Mar. Freshw. 2019, 55, 51–58. [Google Scholar] [CrossRef] [Green Version]
  57. Zhao, Y.; Zhang, J.; Lin, F.; Ren, J.S.; Sun, K.; Liu, Y.; Wu, W.; Wang, W. An ecosystem model for estimating shellfish production carrying capacity in bottom culture systems. Ecol. Model. 2019, 393, 1–11. [Google Scholar] [CrossRef]
  58. Nan, X.; Wei, H.; Zhang, H.; Nie, H. Factors influencing the interannual variation in biomass of bottom-cultured yesso scallop (Patinopecten yessoensis) in the Changhai Sea Area, China. Front. Mar. Sci. 2022, 8, 798359. [Google Scholar] [CrossRef]
Figure 1. (a) Bathymetry of the Yellow Sea. (b) Inset of the outlined square in (a), indicating the location of the mooring station (red triangle). (c) Instrument configuration schematic of the T-chain at Stns. A4 and B4.
Figure 1. (a) Bathymetry of the Yellow Sea. (b) Inset of the outlined square in (a), indicating the location of the mooring station (red triangle). (c) Instrument configuration schematic of the T-chain at Stns. A4 and B4.
Remotesensing 14 03660 g001
Figure 2. (a) Time-depth variations of observed temperature at Stn. A4. The black curve represents sea surface elevation. (b) Temporal variation of temperature as monitored by four temperature sensors. (cf) Close views of the high-frequency variation of the temperature.
Figure 2. (a) Time-depth variations of observed temperature at Stn. A4. The black curve represents sea surface elevation. (b) Temporal variation of temperature as monitored by four temperature sensors. (cf) Close views of the high-frequency variation of the temperature.
Remotesensing 14 03660 g002
Figure 3. The same as in Figure 2 but for Stn. B4.
Figure 3. The same as in Figure 2 but for Stn. B4.
Remotesensing 14 03660 g003
Figure 4. Temperature-frequency spectra of daily temperature variations at (a) Stn. A4 and (b) Stn. B4. The frequency spectra during the two-month-long observations at (c) Stn. A4 and (d) Stn. B4. For reference, the frequencies of the essential components are indicated by horizontal dashed lines.
Figure 4. Temperature-frequency spectra of daily temperature variations at (a) Stn. A4 and (b) Stn. B4. The frequency spectra during the two-month-long observations at (c) Stn. A4 and (d) Stn. B4. For reference, the frequencies of the essential components are indicated by horizontal dashed lines.
Remotesensing 14 03660 g004
Figure 5. (a) Time-depth variations of temperature from 11 July 2019 at 7:51 to 12 July 2019 at 3:51 (UTC + 8). (b) The corresponding temporal variations of temperature. (c) Inset of the shaded area in (a) showing the typical ISWs marked with a black line. (d) SAR image was acquired over the northern YS on 11 July 2019 (05:51 UTC + 8). Red stars denote mooring stations. Notably, the x-axis has been changed to hours after this SAR image time.
Figure 5. (a) Time-depth variations of temperature from 11 July 2019 at 7:51 to 12 July 2019 at 3:51 (UTC + 8). (b) The corresponding temporal variations of temperature. (c) Inset of the shaded area in (a) showing the typical ISWs marked with a black line. (d) SAR image was acquired over the northern YS on 11 July 2019 (05:51 UTC + 8). Red stars denote mooring stations. Notably, the x-axis has been changed to hours after this SAR image time.
Remotesensing 14 03660 g005
Figure 6. Same as in Figure 5 but for another typical SAR image. The SAR image (ad) was acquired on 8 August 2019 (17:46 UTC + 8).
Figure 6. Same as in Figure 5 but for another typical SAR image. The SAR image (ad) was acquired on 8 August 2019 (17:46 UTC + 8).
Remotesensing 14 03660 g006
Figure 7. Statistical characteristics of ISWs. (a,d) The temporal variation of surface elevation (black line) and occurrence of ISWs (red points) relative to the time of low tide. (b,e) RMS vertical velocity in one-hour (grey curves) and one-day intervals (black curves). The red curve represents the 1.8 × RMS variations of the sea surface. (c,f) The number of observed ISWs during a fortnightly cycle, binned by days before/after the time of maximum tidal range.
Figure 7. Statistical characteristics of ISWs. (a,d) The temporal variation of surface elevation (black line) and occurrence of ISWs (red points) relative to the time of low tide. (b,e) RMS vertical velocity in one-hour (grey curves) and one-day intervals (black curves). The red curve represents the 1.8 × RMS variations of the sea surface. (c,f) The number of observed ISWs during a fortnightly cycle, binned by days before/after the time of maximum tidal range.
Remotesensing 14 03660 g007
Figure 8. The occurrence of ISWs relative to the tidal phase at (a) Stn. A4 and (b) Stn. B4. Black curves represent the time series of the tidal elevation simulated using the Ohio State University Tidal Prediction Software (OTPS). The number of ISWs is indicated by vertical bars.
Figure 8. The occurrence of ISWs relative to the tidal phase at (a) Stn. A4 and (b) Stn. B4. Black curves represent the time series of the tidal elevation simulated using the Ohio State University Tidal Prediction Software (OTPS). The number of ISWs is indicated by vertical bars.
Remotesensing 14 03660 g008
Figure 9. (a) Stack plots of the temperature at Stn. B4 from 11–22 July 2019 (UTC + 8). Red arrows indicate the arrival times of the ISW packets, and the white line indicates the surface elevation. (b) The corresponding time series of the barotropic tide at Stn. B4.
Figure 9. (a) Stack plots of the temperature at Stn. B4 from 11–22 July 2019 (UTC + 8). Red arrows indicate the arrival times of the ISW packets, and the white line indicates the surface elevation. (b) The corresponding time series of the barotropic tide at Stn. B4.
Remotesensing 14 03660 g009
Figure 10. (a) The spatial distribution of the depth-integrated internal-tide-generating body force. The superposed curves represent the ISWs detected by SAR images during the observation periods of 2018 and 2019. ISWs propagating eastward and westward are shown in red and black, respectively. (b) The three-dimensional topography shows the surrounding study area. The ISW generation site, as traced by the ray-tracing method, is represented by S1. The locations of mooring stations are also indicated.
Figure 10. (a) The spatial distribution of the depth-integrated internal-tide-generating body force. The superposed curves represent the ISWs detected by SAR images during the observation periods of 2018 and 2019. ISWs propagating eastward and westward are shown in red and black, respectively. (b) The three-dimensional topography shows the surrounding study area. The ISW generation site, as traced by the ray-tracing method, is represented by S1. The locations of mooring stations are also indicated.
Remotesensing 14 03660 g010
Figure 11. (a) SAR image was taken on 11 July 2019 (05:51 UTC + 8). The red circle indicates the estimated generation site (S1) traced by the ray-tracing method. The red star denotes location of Stn. B4. (b) The propagation process of ISWs simulated by the gKdV model. The dashed line indicates the topography. (c) The corresponding tidal current at S1 at the moment of ISW generation. Colors represent the surface elevation, and the arrow indicates the direction of tidal variations.
Figure 11. (a) SAR image was taken on 11 July 2019 (05:51 UTC + 8). The red circle indicates the estimated generation site (S1) traced by the ray-tracing method. The red star denotes location of Stn. B4. (b) The propagation process of ISWs simulated by the gKdV model. The dashed line indicates the topography. (c) The corresponding tidal current at S1 at the moment of ISW generation. Colors represent the surface elevation, and the arrow indicates the direction of tidal variations.
Remotesensing 14 03660 g011
Figure 12. (a) Temporal variation of the wind speed (purple line) and buoyancy flux (green line). (b) Time-depth variation of ε with the isotherms (solid lines) superposed. White lines represent pycnocline boundaries. (c) The replotted ε around the pycnocline (black bars) and temperature. (d) Temporal variation of the averaged ε within the pycnocline (green line) and 1 h RMS vertical velocity (purple line).
Figure 12. (a) Temporal variation of the wind speed (purple line) and buoyancy flux (green line). (b) Time-depth variation of ε with the isotherms (solid lines) superposed. White lines represent pycnocline boundaries. (c) The replotted ε around the pycnocline (black bars) and temperature. (d) Temporal variation of the averaged ε within the pycnocline (green line) and 1 h RMS vertical velocity (purple line).
Remotesensing 14 03660 g012
Table 1. Stations, sensors, frequency, and sensor height of the mooring observations.
Table 1. Stations, sensors, frequency, and sensor height of the mooring observations.
StationsDepth (m)SensorsSampling Period (s)Sensor Height (m)
A439T11.5
T19
T117
T124.5
B438.8CTD3001.5
TD108
T1014
T1018
CTD1021
TD1025
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, H.; Yang, W.; Wei, H.; Jiang, C.; Liu, C.; Zhao, L. On Characteristics and Mixing Effects of Internal Solitary Waves in the Northern Yellow Sea as Revealed by Satellite and In Situ Observations. Remote Sens. 2022, 14, 3660. https://doi.org/10.3390/rs14153660

AMA Style

Liu H, Yang W, Wei H, Jiang C, Liu C, Zhao L. On Characteristics and Mixing Effects of Internal Solitary Waves in the Northern Yellow Sea as Revealed by Satellite and In Situ Observations. Remote Sensing. 2022; 14(15):3660. https://doi.org/10.3390/rs14153660

Chicago/Turabian Style

Liu, Heping, Wei Yang, Hao Wei, Chengfei Jiang, Changgen Liu, and Liang Zhao. 2022. "On Characteristics and Mixing Effects of Internal Solitary Waves in the Northern Yellow Sea as Revealed by Satellite and In Situ Observations" Remote Sensing 14, no. 15: 3660. https://doi.org/10.3390/rs14153660

APA Style

Liu, H., Yang, W., Wei, H., Jiang, C., Liu, C., & Zhao, L. (2022). On Characteristics and Mixing Effects of Internal Solitary Waves in the Northern Yellow Sea as Revealed by Satellite and In Situ Observations. Remote Sensing, 14(15), 3660. https://doi.org/10.3390/rs14153660

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