
Citation: | Zhenpeng Zhang, Chaoya Wu, Shuai Shao, Wei Liu, En-Tao Wang, Yan Li. Diversity of protease-producing bacteria in the Bohai Bay sediment and their extracellular enzymatic properties[J]. Acta Oceanologica Sinica. doi: 10.1007/s13131-020-1589-x |
Internal waves (IWs) are ubiquitous in the ocean and play important roles in energy transfer, biological processes, marine production activities, and ocean engineering. Because IWs frequently occur in the South China Sea (SCS), observing these waves in this region is highly significant (Huang et al., 2016).
In-situ observations, such as acoustic Doppler current profiles and conductivity–temperature–depth, are widely used to observe IWs. These strategies enable the observation of the amplitude, period, phase velocity, polarity, and other characteristics of IWs with high accuracy (Briscoe, 1975) but are usually limited to specific locations. Moreover, they cannot determine the characteristics of IWs at different locations.
Developments in satellite remote sensing have enabled the observation of IWs using optical and microwave satellite images. Synthetic aperture radar (SAR) can obtain high-resolution (several meters) sea surface images, rendering it suitable for observing IWs (Yang et al., 2001); additionally, it has the advantage of being able to conduct all-day, all-weather operations. Alpers (1985) presented a fundamental radar imaging theory of IWs for SAR in 1985, i.e., IW leads to sea-surface convergence or divergence, which changes the surface roughness, thereby affecting the gray level of the radar image. Many characteristics of IWs have been estimated in the literature based on this study, such as the velocity of IWs (Li et al., 2000), the characteristic half-widths of IWs (Zheng et al., 2001), the amplitudes of IWs (Jia et al., 2019), and the spatial distribution of IWs (Liu et al., 1998; Ning et al., 2020). The moderate resolution imaging spectroradiometer (MODIS) has a large swath area of several hundred kilometers, fine spatial resolution, and near-daily global coverage; hence, it is commonly used to study the IWs, in particular, to study characteristics such as the occurrence frequency (Jackson, 2007; Wang et al., 2011), propagation and sources (Sun et al., 2019; Wang et al., 2013), and spatial distribution (Zhou et al., 2016) of IWs.
X-band marine radar has high spatiotemporal resolution and can obtain sea surface images in real-time. Thus, this technology may help determine the evolution of an IW packet (Kropfli et al., 1999). The imaging mechanism of X-band marine radar is similar to that of SAR; specifically, it detects modulations in sea surface roughness associated with currents induced by IWs (Xue et al., 2013). Watson and Robinson (1990) studied the evolution of IWs in the Strait of Gibraltar using X-band marine radar images. Orr and Mignerey (2003) used X-band marine radar to track IWs in the SCS. Ramos et al. (2009) used Radon transform to calculate IW parameters, such as the propagation direction of IWs, their nonlinear velocity, the distance between solitary waves, and the number of solitary waves in each packet, from radar backscatter images. Lü et al. (2010) employed Radon transform to retrieve the propagation direction and velocity of IW packets in the northeastern SCS and then estimated the mixed layer depth by using the two-layer finite depth model. Plant et al. (2011) found at least two different types of IWs in the SCS by analyzing the images of shipborne dual-polarized coherent X-band radar deployed in the SCS. Pan and Jay (2009) used X-band radar images to track IW packets in the Colombia River and subsequently analyzed the characteristics of the detected IWs in combination with the KdV equation. Zha et al. (2012) used X-band radar images in combination with the Morison empirical formula to estimate the force exerted by IWs on a cylindrical pile. Lund et al. (2013) proposed an automated retrieval method of IWs from radar images and discussed the characteristics of IWs with amplitudes greater than 20 m. Badiey et al. (2016) estimated the speed, propagation direction, and amplitude of IWs by analyzing the data obtained from moored thermistor chains and an X-band radar.
Previous studies mostly estimated the average parameters of IWs, such as their average phase velocity and propagation direction (Ramos et al., 2009); unfortunately, these parameters cannot determine the phase velocity of different parts of the waves, which are important for studying the evolution of IWs. However, the distribution of IWs is complex as the latter may show significant changes over areas with a varied topography. Thus, obtaining the small-scale evolution of IWs using traditional methods is challenging. This study proposes a novel method to extract the phase velocity vector of each point on the crest of IWs using X-band marine radar images. The rest of this paper is organized as follows. Section 2 introduces the experimental data, Section 3 presents the method of extracting IW parameters from X-band radar images, and Section 4 verifies and discusses the method. Finally, the conclusions are presented in Section 5.
An experiment was performed in September 2017 on an oil platform in the SCS. The observed sea area lies ~100 km south of Hainan Island (Fig. 1), where the IWs appear frequently (Liang et al., 2019; Jia et al., 2021). The water depth of the region is ~90 m, and the water depth gradually decreases from southeast to northwest. The tides in the observed sea area are mainly irregular diurnal tides. IWs often appear and propagate from southeast to northwest or from east to west (Wang et al., 2013). The X-band marine radar used in the experiment is developed based on a commercial Furuno Far-2127 shipborne navigation radar. Data acquisition software is developed to convert the backscatter from the sea surface into an 8-bit gray value image (Chen et al., 2019). The main configurations of the radar are listed in Table 1. The antenna rotation speed is 42 r/min, and hence, the temporal resolution of the radar is 1.43 s. A single radar image sequence contains 32–128 radar images, and the time difference between two adjacent image sequences is 3–10 min. The radial resolution of the radar is 3.75 m, and the measurement ranges were ~6 km. As an example, Fig. 2 shows two images collected by the X-band radar when IWs pass; here, the wide adjacent bright and dark bands reflect radar echoes caused by changes in sea surface roughness due to IWs, and small wind waves are superimposed on the wide bands. The wave crests of the IW packet are regular in Fig. 2a but intersect in Fig. 2b. Thus, the variation of IWs in this area is complex.
Parameter | Value |
Polarization | HH |
Radar frequency/MHz | 9 410±30 |
Pulse width/ns | 70 |
Pulse repetition frequency/Hz | 3 000 |
Transmit power/kW | 25 |
Beamwidth (3 dB)/(°) | ~1 horizontally, ~20 vertically |
Antenna | 2.4 m slotted waveguide antenna |
Antenna rotation speed/(r·min−1) | 42 |
Sample frequency/MHz | 40 |
Given the limitations of the experimental conditions, no in-situ measurement of IWs was conducted. Instead, synchronous MODIS images of the region were selected to verify the IWs recorded by the X-band marine radar images. The MODIS sensors are onboard the National Aeronautics and Space Administration’s Earth Observing System satellites Terra and Aqua (Jackson, 2007). Terra and aqua transit at 10:30 (local time) and 13:30 (local time) every day. The MODIS has a swath width of 2 330 km and a spatial resolution of 250 m; it can be used to observe ocean IWs at a large scale. Because the optical satellite images are easily blocked by clouds, only two groups of matched X-band marine radar images and MODIS images were obtained.
This section describes the method of extracting IW parameters from X-band marine radar images, including the phase velocities on the crests of the IWs, the number of solitary waves, and their distances in the IW packet.
The radar image sequence is first averaged to reduce the influence of wind waves on the X-band marine radar images. For an X-band marine radar image sequence,
Two adjacent radar image sequences must be selected, i.e.,
The backscatter intensity of X-band radar decreases with increasing distance (Lund et al., 2013). Thus, removing the ramp of gray values with distance is necessary to extract the characteristics of IWs effectively.
The power function is used to fit the gray values of the averaged radar image for each
$$ {\bar I_{{\rm{f}}i}}\left( {r,{\varphi _0},t} \right) = a \cdot {r^b}\quad \left( {i = 1,2} \right) , $$ | (1) |
where a and b are coefficients that can be determined using the least-squares method.
The ramp of the radar images is then corrected as follows:
$$ {\bar I_{{\rm{c}}i}}\left( {r,{\varphi _0},t} \right) = {\bar I_i}\left( {r,{\varphi _0},t} \right) - {\bar I_{{\rm{f}}i}}\left( {r,{\varphi _0},t} \right)\quad \left( {i = 1,2} \right) . $$ | (2) |
Following radar backscatter ramp correction, the gray values of the radar image are adjusted to the range of 0–255 as follows:
$$ {\bar I_{{\rm{c}}i}}\left( {r,\varphi ,t} \right) = \frac{{{{\bar I}_{{\rm{c}}i}}\left( {r,\varphi ,t} \right) - \min \left( {{{\bar I}_{{\rm{c}}i}}\left( {r,\varphi ,t} \right)} \right)}}{{\max \left( {{{\bar I}_{{\rm{c}}i}}\left( {r,\varphi ,t} \right)} \right) - \min \left( {{{\bar I}_{{\rm{c}}i}}\left( {r,\varphi ,t} \right)} \right)}} \times 255 . $$ | (3) |
Two-dimensional fast Fourier transform (2D-FFT) is applied to the radar images to obtain the average propagation direction of the IWs. The 2D wavenumber spectrum of the IWs is then obtained. The average propagation direction of the IW packet can be obtained using
$$ {\varphi _{\rm{p}}} = \arctan \left( {\frac{{{k_{y\max }}}}{{{k_{x\max }}}}} \right) , $$ | (4) |
where the wavenumbers
According to Eq. (4), the propagation direction of an IW has an ambiguity of 180°. Thus, cross-spectral analysis of the two adjacent radar images is used to determine the actual propagation direction of IWs (Chen et al., 2019).
The radar radial profiles at
$$ {v_i} = \frac{{\Delta {r_i}}}{{\Delta t}}\quad \left( {i = 1,2, \ldots } \right) , $$ | (5) |
where
Owing to the small observation range of X-band marine radar, a radar image may not cover all ISWs of an IW packet. Thus, the number of ISWs and the average phase velocity of the IW packet may have to be determined using multiple X-band radar image sequences.
The optical flow method is widely used in moving-target tracking and dynamically analyzes images according to the characteristics of each pixel (Horn and Schunck, 1981). It is used herein to track IWs and extract the phase velocity vector on different locations of the IW. In this method, the gray value of the image before and after motion remains unchanged as a prerequisite,
$$ \frac{{{\rm d} \bar I}}{{{\rm d} t}} = 0 . $$ | (6) |
Using the chain rule for differentiation, we find that the equation of optical flow for a gray image observed by X-band marine radar is
$$ {\bar I_x}u + {\bar I_y}v + {\bar I_t} = 0 , $$ | (7) |
where
$$ {E_{rr}}^2 = \iint {\left[ {{{\left( {{I_x}u + {I_y}v + {I_t}} \right)}^2} + \text{λ} \left( {u_x^2 + u_y^2 + v_x^2 + v_y^2} \right)} \right]{\rm d}x{\rm d}y} , $$ | (8) |
which is usually replaced by the following equation in practice,
$$ {E_{rr}}^2 = \iint {\left[ {{{\left( {{I_x}u + {I_y}v + {I_t}} \right)}^2} + \text{λ} \left( {{{\left( {u - \bar u} \right)}^2} + {{\left( {v - \bar v} \right)}^2}} \right)} \right]{\rm d}x{\rm d}y} , $$ | (9) |
where λ is an unknown weighting factor because the value of it needs to consider the noise in the X-band radar image. Solving this extreme problem, we can determine the velocity of each point in the image as follows:
$$ \left\{\begin{array}{l}{u^{n + 1}} = {{\bar u}^n} - {{\bar I}_x}\dfrac{{{{\bar I}_x}{{\bar u}^n} + {{\bar I}_y}{{\bar v}^n} + {{\bar I}_t}}}{{{\text{λ} ^2} + \bar I_x^2 + \bar I_y^2}},\\ {v^{n + 1}} = {{\bar v}^n} - {{\bar I}_y}\dfrac{{{{\bar I}_x}{{\bar u}^n} + {{\bar I}_y}{{\bar v}^n} + {{\bar I}_t}}}{{{\text{λ} ^2} + \bar I_x^2 + \bar I_y^2}}.\\ \end{array}\right. $$ | (10) |
The average phase velocity (Eq. (5)) estimated from the radial profiles at
First, the two adjacent radar images
Second, the phase velocity of each ISW in the wave packet is estimated using the method described in Section 3.2, and the average phase velocity of the IW is used as the reference phase velocity. For convenience, each peak in the radial profile at
Third, the weighting factor is set to change from
The method to estimate IW parameters from X-band marine radar images is shown in Fig. 3.
Different types of IW images observed by X-band radar (i.e., a set of single IW packet and a set of crossed ISWs) are used to estimate the IW parameters and verify the proposed method for retrieving IW parameters from X-band radar images. Because in-situ measurements are not conducted in the experiment, the MODIS images are used to verify the distribution of IWs. The phase velocities estimated by the proposed optical flow-based method are verified using the average phase velocity of the reference point.
A MODIS image acquired on September 2, 2017, at 03:30 UTC is shown in Fig. 4; in the figure, multiple groups of IW packets may be observed. A group of IW packets located approximately 9.1 km to the west of the radar station could be observed to propagate from east to west. According to the empirical formula in (Tang, 2019), the phase velocity of IWs in the region is approximately 0.76 m/s. Thus, the IW packets would arrive at the radar at ~00:10 UTC on September 1, 2017. The X-band radar image sequence recorded on September 1, 2017, at 23:15 UTC matches the IW in the MODIS image. As shown in Fig. 2a, the propagation direction of the IW packet in the X-band radar image agrees with that determined from the MODIS image (Fig. 4), i.e., the stripes in the X-band marine radar image sequence reflect the IW packet observed by MODIS.
The X-band radar image sequence is averaged using the method described in Section 3.1.1, as shown in Fig. 5. Average processing smoothens the wind waves and swells in the image, thereby facilitating the identification of IW stripes.
The gray values are fitted with the distance determined by the radar by using Eq. (1) to perform ramp correction on the radar image. Figure 6 shows a fitting curve of this group of images. The R-squared value of the fitting results is 0.999 3, which means the selected power function model is applicable to the present case.
The 2D wavenumber spectrum is obtained via 2D FFT of the radar image in the study area (Fig. 7). The propagation direction of the IW in Fig. 7 can be calculated as 90° or 270° (Eq. (4)), and the actual propagation direction of the IW packet can be determined. Cross-spectral analysis reveals that the IW packet propagates from east to west, as shown by the red dashed line in Fig. 7.
After determining the propagation direction of the IW, the radial profiles of the two adjacent radar image sequences are selected at the azimuth of this direction, as shown in Fig. 8. The IW packet contains four ISWs, the distances between which are 219 m, 158 m, and 227 m. The phase velocities at each reference point of the IW packet are 1.01 m/s (leading wave), 1.09 m/s, 1.14 m/s, and 1.18 m/s (Eq. (5)). Typically, the leading wave in an IW packet shows the maximum amplitude; thus, the backscatter intensity caused by the leading wave is largest in the radar image while those caused by the latter ISWs gradually decrease, consistent with the changes in the gray values in Fig. 8.
Histogram equalization is then applied to the two average adjacent radar images (Fig. 9a) to reduce the attenuation of gray values with the range. The images are converted into binary images by selecting a threshold gray value (Fig. 9b), and regions with less than 50 pixels are removed from the binary images (Fig. 9c) to enhance the extraction of IW information.
The value of
Radial profiles | Optical flow-based method | X-band radar (Orr and Mignerey, 2003) | SAR observations (Liang et al., 2019) | Empirical formula (Tang, 2019) | |
Reference point A2 | 1.09 m/s | 1.11 m/s | 1.17 m/s | 0.89 m/s | 0.76 m/s |
Reference point A3 | 1.14 m/s | 1.06 m/s | 1.17 m/s | 0.89 m/s | 0.76 m/s |
Reference point A4 | 1.18 m/s | 1.11 m/s | 1.17 m/s | 0.89 m/s | 0.76 m/s |
The MODIS image obtained on September 11, 2017, at 03:24 UTC is shown in Fig. 11. Two groups of IW packets may be observed at ~13 km to the east and southeast of the radar station, and these packets propagate from east to west and from southeast to northwest, respectively. According to the empirical formula proposed by Tang (2019), the IW packet will arrive at the radar station at ~08:09 UTC on September 11, 2017. The matched X-band radar image sequence was recorded on September 11, 2017, at 07:12 UTC, as shown in Fig. 2b, and the propagation directions of the two IW packets detected agree with those shown in the MODIS image (Fig. 11). In other words, X-band marine radar correctly observed the two groups of IW packets.
The X-band radar image sequences are first preprocessed, after which 2D-FFT is used to obtain the IW wavenumber spectrum, which shows two main spectral peaks. The propagation directions of the two ISWs, which are
The radial profiles
After the radar images are processed via histogram equalization and noise reduction, λ is determined using the phase velocity at reference point B1 to be λ0 = 1.18. The phase velocities of the points on the IW crest are obtained by the optical flow-based method with the selected λ0, as shown in Fig. 14.
Figure 14 shows that the proposed method could detect the propagation directions of the points on the crest of the IWs accurately and that the distribution of phase velocities is reasonable. Some errors in the region of approximately x=0.6 and y=1–2 km may be detected because noise cannot be completely eliminated from the image (Horn and Schunck, 1981). The phase velocity estimated by the optical flow-based method for reference point B2 is 1.15 m/s, which also agrees with the phase velocity of IWs in this region determined in previous studies (Orr and Mignerey, 2003; Liang et al., 2019; Tang, 2019). The phase velocities calculated using different methods are listed in Table 3.
Radial profiles | Optical flow-based method | X-band radar (Orr and Mignerey, 2003) | SAR observations (Liang et al., 2019) | Empirical formula (Tang, 2019) | |
Reference B2 | 1.14 m/s | 1.15 m/s | 1.17 m/s | 0.89 m/s | 0.76 m/s |
Figures 10 and 14 clearly show that the proposed optical flow-based method is applicable not only to single IW packet with parallel ISWs but also to crossed IWs.
In the proposed method, the determination of λ0 is important. To study the sensitivity of the value to the method, the other reference point B2 (Fig. 12), has been chosen to determine the value of λ0. The adjacent two radial profiles along the propagation direction of the ISW B2 are shown in Fig. 15. According to the radial profiles, the average phase velocity at point B2 is estimated to be 1.14 m/s, and the value of λ0 is 1.19. Using the optical flow-based method, the phase velocities of the crossed IWs are estimated, and velocity at the reference point B1 is 0.73 m/s. Compared to the average phase velocity of B1 estimated from the adjacent radial profiles, the relative error of the phase velocity at B1 is 0.01 m/s. Therefore, the proposed method is affected little by the selection of different reference points, i.e., the optical flow-based method is stable.
Phase velocity is one of the most important parameters of IWs. Previous methods of extracting the phase velocity of ISWs mostly estimated the average phase velocity from satellite images. Although the observation range of X-band marine radar is smaller than that of satellite images, X-band marine radar can observe the sea surface in real-time and continuously, which is important for studying the evolution of IWs (Kropfli et al., 1999). Previous methods of extracting the phase velocity of ISWs from X-band radar images were mostly based on Radon transform, which can only estimate the average phase velocity of IWs and cannot obtain the phase velocity of each point on the IWs. A novel method is proposed herein to obtain the phase velocity vector of each point on IWs crest effectively, which is verified using the X-band radar images recorded on an oil platform in the SCS, and the phase velocities estimated by the optical flow-based method are close to those derived from the adjacent radial profiles, as well as the phase velocity of the IWs of the sea area measured by previous studies (Wang et al., 2013; Tang, 2019). Moreover, the proposed method is insensitive to the selection of different reference points; hence, it can reasonably estimate the phase velocities of different parts of IWs. In future research, experiments with in-situ observations of IWs will be performed, and the proposed method will be further evaluated using more data.
Parameter | Value |
Polarization | HH |
Radar frequency/MHz | 9 410±30 |
Pulse width/ns | 70 |
Pulse repetition frequency/Hz | 3 000 |
Transmit power/kW | 25 |
Beamwidth (3 dB)/(°) | ~1 horizontally, ~20 vertically |
Antenna | 2.4 m slotted waveguide antenna |
Antenna rotation speed/(r·min−1) | 42 |
Sample frequency/MHz | 40 |
Radial profiles | Optical flow-based method | X-band radar (Orr and Mignerey, 2003) | SAR observations (Liang et al., 2019) | Empirical formula (Tang, 2019) | |
Reference point A2 | 1.09 m/s | 1.11 m/s | 1.17 m/s | 0.89 m/s | 0.76 m/s |
Reference point A3 | 1.14 m/s | 1.06 m/s | 1.17 m/s | 0.89 m/s | 0.76 m/s |
Reference point A4 | 1.18 m/s | 1.11 m/s | 1.17 m/s | 0.89 m/s | 0.76 m/s |
Radial profiles | Optical flow-based method | X-band radar (Orr and Mignerey, 2003) | SAR observations (Liang et al., 2019) | Empirical formula (Tang, 2019) | |
Reference B2 | 1.14 m/s | 1.15 m/s | 1.17 m/s | 0.89 m/s | 0.76 m/s |
Parameter | Value |
Polarization | HH |
Radar frequency/MHz | 9 410±30 |
Pulse width/ns | 70 |
Pulse repetition frequency/Hz | 3 000 |
Transmit power/kW | 25 |
Beamwidth (3 dB)/(°) | ~1 horizontally, ~20 vertically |
Antenna | 2.4 m slotted waveguide antenna |
Antenna rotation speed/(r·min−1) | 42 |
Sample frequency/MHz | 40 |
Radial profiles | Optical flow-based method | X-band radar (Orr and Mignerey, 2003) | SAR observations (Liang et al., 2019) | Empirical formula (Tang, 2019) | |
Reference point A2 | 1.09 m/s | 1.11 m/s | 1.17 m/s | 0.89 m/s | 0.76 m/s |
Reference point A3 | 1.14 m/s | 1.06 m/s | 1.17 m/s | 0.89 m/s | 0.76 m/s |
Reference point A4 | 1.18 m/s | 1.11 m/s | 1.17 m/s | 0.89 m/s | 0.76 m/s |
Radial profiles | Optical flow-based method | X-band radar (Orr and Mignerey, 2003) | SAR observations (Liang et al., 2019) | Empirical formula (Tang, 2019) | |
Reference B2 | 1.14 m/s | 1.15 m/s | 1.17 m/s | 0.89 m/s | 0.76 m/s |