Next Article in Journal
Characterizing the Macro and Micro Properties of Precipitation during the Landfall of Typhoon Lekima by Using GPM Observations
Next Article in Special Issue
Typhoon Early Warning and Monitoring Based on the Comprehensive Characteristics of Oceanic and Ionospheric Echoes from HFSWR: The Case of Typhoon Muifa
Previous Article in Journal
The Preparation Phase of the 2023 Kahramanmaraş (Turkey) Major Earthquakes from a Multidisciplinary and Comparative Perspective
Previous Article in Special Issue
Ship Formation Identification with Spatial Features and Deep Learning for HFSWR
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Surface Vector Current Retrieval by Single-Station High-Frequency Surface Wave Radar Based on Ocean Dynamics in the Taiwan Strait

1
College of Meteorology and Oceanography, National University of Defense Technology, Changsha 410073, China
2
PLA Army Chemical Defense College, Beijing 102200, China
3
School of Electronic Information, Wuhan University, Wuhan 430072, China
4
Collaborative Innovation Center of Geospatial Technology, Wuhan 430072, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Remote Sens. 2024, 16(15), 2767; https://doi.org/10.3390/rs16152767
Submission received: 21 June 2024 / Revised: 21 July 2024 / Accepted: 27 July 2024 / Published: 29 July 2024
(This article belongs to the Special Issue Innovative Applications of HF Radar (Second Edition))

Abstract

:
In order to address the issue of limited common coverage and high cost in mapping ocean surface vector current by two (or more) high-frequency surface wave radars, this paper proposes a single-station surface wave radar vector current inversion algorithm. The feasibility of this algorithm has been validated in the Taiwan Strait. Based on the ocean dynamic characteristics of the Taiwan Strait, the algorithm utilizes the radial current obtained from a high-frequency surface wave radar (HFSWR) in Fujian Province to invert the ocean surface vector current. The surface vector current can be decomposed into three primary dynamic components: tidal currents, wind-driven currents, and geostrophic currents. Firstly, tidal current forecasting models and Ekman and Stokes theories are used to calculate the tidal and wind-driven currents in the Taiwan Strait, respectively. Subsequently, the directions of geostrophic currents in the Taiwan Strait are determined with sea surface height data, and the magnitudes of the geostrophic currents are constrained using the radial current from the single HFSWR. Finally, the three components are added together to obtain the vector current. Comparative results demonstrate that the efficacy of the algorithm has been validated through field experiments (with two HFSWRs and two drifting buoys) conducted in the southwestern of the Taiwan Strait. Further research is needed on the applicability of this algorithm to other sea areas and monitoring systems.

1. Introduction

The Taiwan Strait is an important maritime passage, and observations of surface currents in this region play a significant guiding role in various aspects, such as navigational safety, maritime search and rescue, marine fisheries, pollution monitoring, and the development of marine resources [1,2,3,4]. In recent decades, high-frequency surface wave radars (HFSWRs) have been used as over-the-horizon radars for monitoring oceanic currents. Both research teams from mainland China and Taiwan have deployed several HFSWRs to monitor surface currents in the Taiwan Strait. These include the marine observation demonstration area of the Fujian Provincial Ocean Forecasting Station, which consists of two radar stations in Dongshan and Longhai, and the Taiwan Ocean Radar Observing System (TOROS), which has 19 HFSWRs deployed along the coasts of Taiwan and Penghu Islands [5,6,7,8,9]. Their observation coverage spans nearly the entire Taiwan Strait.
According to the Doppler principle, a single HFSWR can only detect the radial component of the current. This means that at least two HFSWRs are required to obtain the vector current within their common coverage area [10,11,12,13,14]. However, this poses two challenges: it increases the demands on the layout and location selection of radar arrays, leading to higher costs, and it restricts the detection range of the vector currents to jointly covered areas. Therefore, estimating the vector currents solely based on the radial currents detected by a single HFSWR is a significantly valuable and critical area of focus. Previous studies have proposed various algorithms to address these challenges and estimate vector currents based on radial currents. These algorithms can be categorized into two types: those based on algebraic-geometric relationships and those based on fluid dynamic constraints. They provide new approaches for inverting vector currents, expanding the potential applications of HFSWRs.
The major examples include the Parameter Estimation Algorithm [15], Stream Function Method [16], Optimal Interpolation (OI) method [17], Vortex Identification Method [18], and Stream Potential Function Method [19]. The Parameter Estimation Algorithm fundamentally utilizes the least squares method to estimate the magnitude and direction parameters of vector currents. It has been previously applied in shore-based and shipboard radar systems [20,21], yielding positive results. However, the algorithm’s limitations are based on the assumption of spatially uniform ocean currents, resulting in significant errors under moderate to high sea conditions. The Stream Function (SFM) is contingent upon the assumption of horizontally non-divergent ocean dynamics, while the Stream Potential Function (SPFM) incorporates horizontally divergent potential functions, making it more suitable for realistic oceanic scenarios compared to SFM. Marmain et al. achieved vector current mapping from a single HFSWR by combining non-divergent ocean currents with the Vortex Identification Method (VIM), yet lacked experimental validation [18]. Additionally, Wang et al. proposed a method based on sparse representation and unitary transformation techniques for estimating vector currents based on the stream function and verified the algorithm through simulations, but did not conduct field experiments [22]. Furthermore, Igor used radial currents as one of the inputs, significantly improving the assimilated model output of oceanic currents [23]. Han et al. developed an inversion algorithm by combining variational adjoint assimilation methods and the POM numerical model, and achieved promising results in field experiments; however, the duration of comparative validation was limited [24].
This paper presents a radar vector current retrieval algorithm based on ocean dynamic characteristics in the Taiwan Strait. Previous research suggests that the currents in the Taiwan Strait consist of three main components: tidal currents, geostrophic currents, and wind-driven currents, which can be written as below,
V = V t + V w + V g u = u t + u w + u g v = v t + v w + v g ,
where V is the total vector current, V t is the tidal current, V w is the wind-driven current and V g is the geostrophic current. Tidal currents are generated by the tide-generating forces of the sun and moon [25,26,27,28], while geostrophic currents result from the balance between the horizontal pressure gradient force and the Coriolis force, influenced by sea surface heights (SSH) [25]. Wind-driven currents are caused by friction between the sea surface wind and the sea surface. Theoretically, sea surface vector currents can be directly calculated using dynamic theory, eliminating the need for remote sensing methods. However, ideal assumptions and other factors in dynamic theory may lead to variations in actual sea surface currents. Therefore, we propose combining dynamic theory with single-station HFSWR observations to invert the surface vector currents and verify the effectiveness of the algorithm using the vector currents measured by two HFSWRs and drifters in the field experiment.
This paper is organized as follows. Section 2 presents the dataset and methodology utilized to conduct the research. Section 3 presents the results. A discussion is provided in Section 4, followed by the conclusions in Section 5.

2. Data and Methods

2.1. HFSWR Experimental

Filed observations were conducted using two HFSWRs deployed along the coast of Dongshan and Longhai in Fujian Province from 22 August 2021 to 27 August 2021 (UTC+8). The HFSWRs operated at the 8.75 MHz frequency, with an output of the ocean surface currents at a spatial resolution of 5 km every 10 min. The radial currents derived in this study represent the projected vector currents toward the Dongshan station. Specific parameters of HFSWRs can be found in the literature [29].
Within the combined coverage area of these two radars, drifters were deployed from vessels to validate the effectiveness of the proposed algorithm in this paper. The drifters used were LTWDB01 long endurance surface drifters (as shown in Figure 1b), capable of floating on the ocean surface for extended periods to measure the surface ocean current speed and direction. These buoys weigh no more than 5 kg and utilize the BeiDou II satellite positioning system, providing an accuracy of 10 m and a data update rate of 5 min. Figure 1c depicts the trajectories of the two drifters. It is evident that the deployment positions of the two drifters differed, and both buoys followed trajectories oriented from the southwest to the northeast. The timestamps for the drifter data can be found in Table 1 (UTC+8). Additionally, Figure 1a shows the wind rose during the experimental period.
There is no moored buoy in this experiment, and there is no ADCP data available for verification in this sea area. Our proposed algorithm needs to be verified by comparing it with HFSWR current products and drifter currents. In this study, Dongshan radar station measurements were used to constrain and adjust the speed of the geostrophic current.
Prior to employing the radial current to constrain and adjust the geostrophic current, an assessment is conducted on the HF radar observation currents each time, assigning a credibility value C. We utilize a scale of 0 to 1 to denote the credibility (C) of the ocean currents at each spatial point at each time, where 0 signifies low confidence, and 1 signifies high confidence.
The HFSWR utilizes the first-order spectrum in the echo range Doppler spectrum to retrieve ocean surface currents. The signal-to-noise ratio (SNR) in the first-order spectrum serves as a critical indicator of radar signal quality. A high SNR implies that the signal strength significantly exceeds the background noise, directly impacting the extraction of velocity information using Doppler frequency shifts. Moreover, in scenarios with elevated noise levels, the signal may be engulfed by noise, rendering the effective detection of Doppler frequency shifts unfeasible.
The performance of the Multiple Signal Classification (MUSIC) algorithm in estimating the Direction-of-Arrival (DOA) is crucial. By applying spatial filtering and processing to the signals, the MUSIC algorithm can more precisely determine the location of ocean surface currents, thereby facilitating the extraction of directional information of the ocean currents.
To synthesize vector currents in the common coverage area, at least two radar stations positioned at different locations are required. The Geometric Dilution of Precision (GDOP) is also a vital parameter, as it directly influences the accuracy of radar measurements. A lower GDOP value indicates that the radar system can provide more precise information regarding the current speed and direction.
With reference to the method outlined in [30,31], the credibility of the ocean currents, denoted as (C), can be represented by the following formula:
d v 1 = 1 v v 1 v 2 c o s θ 1 θ 2 s i n 2 θ 1 θ 2 Δ v 1 d v 2 = 1 v v 1 v 2 c o s θ 1 θ 2 s i n 2 θ 1 θ 2 Δ v 1 d θ 1 = 1 v v 1 v 2 s i n θ 1 θ 2 v 1 2 + v 2 2 2 v 1 v 2 c o s θ 1 θ 2 s i n 3 θ 1 θ 2 c o s θ 1 θ 2 Δ θ 1 d θ 2 = 1 v v 1 v 2 s i n θ 1 θ 2 + v 1 2 + v 2 2 2 v 1 v 2 c o s θ 1 θ 2 s i n 3 θ 1 θ 2 c o s θ 1 θ 2 Δ θ 2 C = d v 1 + d v 2 + d θ 1 + d θ 2
where, d v 1 and d v 2 are the errors of radial current velocities measured by two radar sites respectively, and d θ 1 and d θ 2 are the errors of DOA. Radial current errors d v 1 and d v 2 are characterized by the first-order spectral SNR corresponding to the single radar site radial current.
In Figure 1c, the color overlay illustrates the spatial distribution of the mean credibility during the experimental period. It can be observed that the central region of the radar footprint exhibits higher credibility, exceeding 0.7, while the edge areas display lower credibility, dropping below 0.4.

2.2. Tidal Current

Seawater experiences a regular horizontal motion caused by tidal forces, known as tidal currents. The current creates an elliptical pattern when connecting the endpoints of the tidal current vector over one tidal cycle. The pattern is used to represent the rotational direction and magnitude of the tidal current. In coastal areas, tidal currents are constrained by the coastline, resulting in very small semi-minor magnitudes of the tidal ellipses. However, in more distant regions from the coastline, surface currents are influenced by the Coriolis force and predominantly exhibit rotary currents. Surface tidal currents dominate most of the region of the Taiwan Strait [25]. During high tide, offshore tidal currents flow into the central Taiwan Strait, with a northeastward current in the seas off the southern coast of Fujian and a northward deviation within the central and southern parts of the strait. Conversely, during low tide, the tidal currents in the Taiwan Strait exhibit an outflow pattern, with a southwestward current in the seas off the southern coast of Fujian and a southward deviation within the central and southern parts of the strait.
It is generally believed that the tidal current is the superposition of cosine oscillations of various tidal components. Therefore, by expanding the tidal force according to frequency, various tidal components can be obtained. Based on the periodicity of these component tides during their propagation and coupled with observational data (such as tidal gauge stations, buoys, and satellite data), tidal currents can be obtained through numerical simulation methods that consider the impact of actual topography. The tidal inversion and prediction system OTIS (The Oregon State University Tidal Inversion System), developed by Oregon State University, assimilates coastal tidal gauge data along with satellite remote sensing data (TPXO8.0-atlas), enabling the output of tides and tidal currents in the northern South China Sea and East China Sea [32] with a spatial resolution of 1/30° in the Taiwan Strait region. The output includes the tidal ellipse parameters of eight major tidal constituents, such as M2, S2, O1, K1, N2, K2, P1, Q1. OTIS has been validated through multiple comparisons, including satellite remote sensing, high-frequency radar, and moored buoys, all of which have produced positive results since its introduction in 2004 [33,34,35].
In this study, OTIS is adopted to obtain the tidal currents in the Taiwan Strait and support the synthesis of surface vector currents. It is noted that the authors performed spatial interpolation on OTIS to maintain consistency between its spatial latitude and longitude grid points and the radial currents detected by the DS high-frequency radar.

2.3. Wind-Driven Current

The annual average wind speed in the Taiwan Strait ranges from 8.0 to 8.5 m/s. Higher wind speeds are generally experienced in the central and southwestern regions, while the northeastern sea areas have slightly weaker winds [36]. From mid-May to August, the prevailing wind is the southwesterly monsoon, with wind directions mainly from the south to the southwest across most sea areas. The distribution of average wind speeds is greater in the north and south and smaller in the central area, ranging from 6.0 to 9.0 m/s. The average wind speed in the central part of the Taiwan Strait, Kinmen, and Penghu ranges from 8.0 to 9.0 m/s. During the period from October to April of the following year, the prevailing wind in the Taiwan Strait is the northeasterly monsoon, with wind directions mainly from the north to the northeast. This is accompanied by strong winds that reach their peak from December to January of the following year. In January, the frequency of northeasterly winds in the Taiwan Strait can reach up to 60% to 78.8%, and the average wind speed during the northeasterly monsoon period is 9.0 to 13.0 m/s [37].
The wind-driven current V w consists of Ekman and Stokes components, which can be computed based on Ekman theory [38,39] and Stokes theory [40], respectively. Ekman components U E can be represented as:
z A z U E z = i f U E U E z = h = 0 A z U E z z = 0 = τ ρ w ,
where, z represents the water depth, f is the Coriolis parameter, A z is the eddy viscosity [41]. The following two formulas are the boundary conditions. Where h is the total water depth, ρ w is the density of seawater. τ = ρ a C D U 10 2 represents wind stress. ρ a is the air density. C D is the drag coefficient. U 10 is wind speed at a height of 10 m above ocean surface. According to Wu et al. [42,43], C D = 0.8 + 0.065 U 10 × 10 3 is suitable for both high and low wind speeds. According to boundary conditions, it is possible to solve U E for different water depths:
U E = τ e i π / 4 ρ w f A z sinh [ 1 + i h z f / 2 A z ] cosh [ 1 + i h f / 2 A z ]
Due to the nonlinearity of surface gravity waves, Stokes drift U S z occurs in the wave propagation direction. U S z can be calculated from the non-directional ocean wave spectrum F ω [44,45]:
U S z = 2 g 0 ω 3 e 2 k ^ z F ω d ω
where, g is the acceleration due to gravity, ω is the angular frequency. k ^ is the unit vector in the direction of wave propagation θ . In this paper, JONSWAP wave spectrum is used to calculate U S z [46]:
F ( k ) J o n s w a p = 0.0081 2 k 4 exp [ 5 4 ( k m k ) 2 ] γ exp [ 1 2 σ 2 k m k k m ) 2 , k m = 22 g 8 / 3 U 19.5 L 2 / 3 ,   σ = 0.07 ,   k k m 0.09 ,   k > k m ,   λ = 3.3
where k m is the peak wave number, U 19.5 is the sea surface wind speed at 19.5 m height, L = 100   km is the fetch length, σ is a parameter that controls the shape of the non-directional ocean wave spectrum.
The surface wind data required for calculating wind-driven currents are the sea level 10 m wind U/V components from the ECMWF reanalysis data ERA5, with a temporal resolution of 1 h and a spatial resolution of 0.25°. The sea surface wind speed U 19.5 at the height of 19.5 m required in Equation (6) can be calculated based on the relationship between wind and height [47]. The wind-driven currents are interpolated to match the spatial and temporal characteristics of the radial current at the DS Station. Distance-weighted interpolation is employed for spatial alignment, while linear interpolation is used for temporal alignment.

2.4. Geostrophic Current

Geostrophic currents, which result from the balance between the horizontal pressure gradient force and Coriolis force, are present not only in open ocean regions but also in nearshore areas, where they contribute to local circulation along with tidal currents. These currents play a significant role in upper ocean mixing and vertical transport [48,49]. A study by Li and Liao found a correlation between sea surface height and winter monsoons, concluding that the persistent influence of monsoons causes sea surface height to rise in the northern Taiwan Bank, similar to how strong winds lift water masses [50,51].
In the southwest region of the Taiwan Strait, offshore and alongshore currents flow eastward or northeastward when the pressure gradient force is balanced with the Coriolis force. According to Oey’s work, the geostrophic balance causes an eastward alongshore flow during the monsoon relaxation period, which then turns north after crossing the Penghu Channel [52]. Wang utilized vector currents obtained from high-frequency radars to verify the timing and extent of stable geostrophic currents. Wang concluded that the geostrophic balance dominates the residual current in the across-strait direction, while in the alongshore direction, both the geostrophic balance and wind stress explain the residual current in this region [53].
Daily datasets were extracted from the Global Ocean Physics Reanalysis Products GLOBAL-REANALYSIS-PHY-001-030 (GLORYS12V1) provided by the Copernicus Marine Environment Monitoring Service (CMEMS, http://copernicus.eu, accessed on 10 September 2023). The GLORYS12V1 product is a CMEMS global ocean eddy-resolving (1/12° horizontal resolution and 50 vertical levels) reanalysis covering the altimetry era of 1993–2024. The calculation of geostrophic currents based on the sea surface height gradient is a well-established theoretical method, demonstrated by Equation (7),
f v = g ζ x   f u = g ζ y ,
where ζ represents the SSH, x, and y are the axes of the Cartesian coordinate system. u and v represent the geostrophic currents calculated using Equation (7) in the equatorial and meridional directions. These values establish the directions of the geostrophic currents within the radar footprint. The spatial distribution of geostrophic currents calculated from sea surface height at UTC 12:00 22 August 2021 is shown in Figure 2, where the majority of the geostrophic currents flow eastward or northeastward. Particularly, the geostrophic current exceeding 30 cm/s south of the DS radar station flows eastward along the southern edge of the Taiwan Bank. Upon nearing the Penghu Islands, it shifts its direction to the northeast and intensifies to 60 cm/s, resembling the formation of a weak cyclonic feature.
Based on the results of the ‘Fujian Coastal Zone Survey’ conducted by the Oceanographic Division of Fujian Provincial Science Committee from 1960 to 1964 and the ‘Comprehensive Survey of the Western Taiwan Strait’ carried out from 1984 to 1985 [48], a persistent northeastward geostrophic current has been identified in the Taiwan Strait [49,50,51]. However, it is important to note that the magnitude of this geostrophic current is not constant, a phenomenon that has also been confirmed by numerous scholars. Assuming that the dynamic theoretical calculations of tidal currents and wind-driven currents are accurate, the determination of the magnitude of the geostrophic current is the only requirement for inferring sea surface vector currents.
In this study, the directions of geostrophic currents existing in the footprint of the HFSWRs are retrieved from SSH. Then, we set the initial geostrophic current as 20 cm/s and introduce a scaling factor η . We determine the magnitude of the geostrophic current using an iterative method. The value of η needs to be determined based on the radial current V E _ s s w . If the value of η is correct, the projection of V t + V w + η V g in the radial direction at the DS Station, denoted as V E _ p r o j , will have the minimum error with respect to V E _ s s w . In this study, the η value is determined using radial currents with credibility greater than 0.5 within the detection range. The iterative process is illustrated in Figure 3. Here, χ 2 represents the cost function:
χ 2 = 1 C ( [ V E _ s s w V E _ p r o j ] 2 )
where V E _ s s w represents the radial current at the DS radar Station, V E _ p r o j denotes the projection of the sea current vector in the radial current direction, and c denotes the total number of grid points where the credibility is equal to or greater than 0.5. We apply this iterative algorithm to determine the value of η at every moment, and the value η is constant for all radial current positions at one moment. Setting η as 3 initially, the range for the geostrophic current search is between 0~60 cm/s. By employing the iterative algorithm, we calculate the η value that minimizes the error χ 2 , which represents the specific η at that moment.
The temporal and spatial resolutions of the three current components are shown in Table 2, as well as those of the HF radar radial currents. To match the radial currents, the three current components are interpolated to match the spatial characteristics of the radial current at the DS radar station. Distance-weighted interpolation is employed for spatial alignment, while linear interpolation is used for temporal alignment.

3. Results

Based on the methods and data introduced in Section 2, we conducted an inversion of the surface vector currents in the southwestern Taiwan Strait. The calculated wind-driven currents, tidal currents, and geostrophic currents at 14:20 on 22 August 2021 (UTC+8) are shown in Figure 4. Wind-driven currents flowed toward the northeast in most regions, and large wind-driven currents exceeding 60 cm/s were found in the northern region, while small ones varied in the central and southern regions. The tidal currents in the Taiwan Bank were relatively strong, with the maximum velocities exceeding 100 cm/s. Following this are the tidal currents in the Xiamen-Penghu Depression located north of the Taiwan Bank, where the tidal velocities were approximately 40 cm/s. The weakest tidal currents are found south of the Taiwan Bank, with speeds of less than 10 cm/s.
The spatial distribution of the ocean surface current credibility at this moment is depicted in Figure 5. For clarity, Figure 5 displays only regions with credibility equal to or greater than 0.5, which covers most of the radar footprint. Areas with credibility below 0.5 primarily correspond to ocean regions farther away from the radar station (e.g., east of the Taiwan Bank).
Given this, we obtain the total retrieval currents by adding these three kinds of currents. In Figure 6, the blue arrows represent the vector currents (HF) measured by HFSWRs, while the red arrows depict the retrieval results (RE). The mean root mean square errors were calculated in the northward and eastward directions according to the method in Hatayama’s study [54], defined as
R M S E u = 1 m 1 m U R E m U H F m 2 R M S E v = 1 m 1 m V R E m V H F m 2
to evaluate the differences. Here U and V are the eastward and northward components of the vector current, and m is the total number of effective currents within the radar footprint.
In the central sea areas covered by radar where HF generally could demonstrate higher accuracy [8], RE shows good agreement with HF, with root mean square errors (RMSE) of 11.5 cm/s in the eastward direction (U component) and 18.2 cm/s in the northward direction (V component). The relative errors of the U and V components are 19.05% and 21.17%, respectively. However, in the marginal zones, the agreement between RE and HF is not as strong as that in the central areas.
Another example of three components of currents is shown in Figure 7. The direction of the wind-driven currents predominantly aligns with the northwesterly direction, with magnitudes that are generally modest, not exceeding a velocity of 30 cm/s. In comparison to the previous example, the tidal currents near the Taiwan Bank and the Penghu Islands are less intense, with maximum values reaching up to 70 cm/s. Within the Xiamen-Penghu Depression, tidal current velocities range between 40 and 50 cm/s.
Similar to Figure 5, Figure 8 also displays only the regions with credibility greater than or equal to 0.5. This area is slightly smaller than that shown in Figure 5, particularly near the DS radar station. The variability in credibility is influenced by numerous factors, including SNR at different times, radio interference, and the operational status of the radar equipment. Subsequently, we compare the ocean currents of the RE and HF within this area.
Figure 9 displays the current pattern at 09:00 on 26 August 2021, showing a southwest-to-northeast direction in a large portion of the sea area. Similar to the previous case, the agreement between RE and HF in the marginal zones was relatively poor. Within the common radar coverage areas, the RMSE of the U component is 11.3 cm/s, and the RMSE of the V component is 17.3 cm/s. The relative errors of the U and V components are 15.62% and 18.62%, respectively.
In order to validate the algorithm, we compared RE and HF with vector current measurements obtained from drifters. Before comparison, spatiotemporal matching was conducted due to the mobility of the drifters. The nearest neighbor principle was used to determine the spatial and temporal closest ocean current data for comparison. The comparisons of U and V components are illustrated in Figure 10. It is evident that both the U and V components exhibit strong periodicity, confirming the dominance of tidal currents in the Taiwan Strait, a region known for strong tidal influence. Furthermore, in Figure 10b,d, both the HF and RE exhibit a closer alignment with the drifters in the V direction. However, the alignment appears less pronounced in the U direction. This could be attributed to significant disturbances in the U direction during the experimental period, as indicated by abrupt changes in the U component, as shown in Figure 10c. Table 3 presents the RMSE of the HF and RE in comparison with the two drifters. The RMSE values listed before and after “/” correspond to drifters 1 and 2, respectively. From Table 3, it is observed that the RMSE of the V component exceeds that of the U component, and the RMSE obtained from the comparison with Drifter2 is higher than that of Drifter1. The accuracy of HF slightly surpasses that of RE, with a maximum difference of not more than 4 cm/s.
Figure 11 presents scatter plots comparing the HF, RE, and drifters in the U and V components, with the color of the scatter points indicating the quantity of samples. Figure 11a–f depict the comparison results with Drifter1, while Figure 11g–l depict the comparison results with Drifter2. The solid black line in each plot represents y = x , with the correlation coefficient (cc) also indicated in each subplot. Their p-values are extremely small, approaching zero, suggesting that the positive correlation between the two is not a random occurrence but should be considered statistically significant. From Figure 11, it is evident that the results for the V component are more concentrated along the black line, resulting in higher correlation coefficients for the V component compared to the U component. When comparing with the drifters, it is observed that, except for the RE correlation coefficients being higher than HF in Figure 11a,b, the remaining plots show higher correlation coefficients between HF and the drifter results. Additionally, the correlation coefficients between RE and HF are higher than those between RE and drifters, indicating that the variations in RE and HF are more closely aligned, given that RE is obtained through the inversion of radial flow results.

4. Discussion

The dynamic theory of tides is mature, and OTIS is widely utilized and substantiated by scholars; thus, we do not alter the tidal current output through iteration. The contribution of wind-driven currents is much smaller than that of tidal currents in the Taiwan Strait during the experimental period, and the theory for wind-driven currents is also mature. Wind datasets exhibit good spatial continuity and high temporal resolution. Consequently, the calculated results for wind-driven currents are relatively accurate. As a result, we also do not alter the wind-driven currents through iteration. Although the dynamic theory of geostrophic current is relatively mature, high-time-resolution sea surface height data are difficult to obtain, and the contribution of geostrophic currents is larger than that of wind-driven currents sometimes in the Taiwan Strait [37]. Therefore, using the radial current as a constraint, we determine the magnitude of the geostrophic currents through an iterative process. In a specific area of the Taiwan Strait, this method of determining geostrophic currents is effective, as most of the geostrophic currents in this area flow east or northeast all year and are not very complex. However, in areas under complex dynamic conditions, more effective methods for determining geostrophic currents are necessary.
The comparison reveals a high correlation between HF, RE, and drifter currents, particularly for the V component, where the correlation between HF and RE exceeds 0.9. Additionally, the lower precision of RE compared to HF is foreseeable, as the precision of RE depends on the detection accuracy of the radial current. Nevertheless, the difference in precision between RE and HF is not significant, with variations of up to 4 cm/s and 1 cm/s for maximum and minimum values, respectively, indicating the excellent performance of the proposed algorithm. This algorithm not only effectively inverses the temporal and spatial variations of surface vector currents but also achieves high accuracy.
Both the HF and RE exhibit clearer periodic characteristics and higher correlations in the V component. However, the U component generally shows periodicity but with a relatively lower correlation. This could be due to the performance differences between the two HFSWRs, as the DS radar station detects the radial component in the U direction for most sea areas. Drifters mainly drift on the ocean surface without the use of drogues; they are predominantly influenced by wind (Stokes drift) and tides. This factor may have contributed to the relatively large errors observed in areas with strong geostrophic currents, such as the eastern region of the DS radar station. Overall, the algorithm performs well in the Taiwan Strait, but its applicability to other ocean areas requires validation using experimental data from these areas. Additionally, the algorithm, developed based on ocean dynamics, may produce significant errors during extreme phenomena, such as typhoons, when classic drift theory-based wind-driven currents may become inaccurate and the contribution of wind-driven currents will significantly increase.
For the same sea area, not only HFSWRs but also X-band radars and satellite-based sensors [55] are capable of detecting the radial components of surface vector currents. Exploring the feasibility of this algorithm in various detection systems holds great potential for future research.
The characteristics of surface currents in the Taiwan Strait exhibit seasonal variations; however, the three primary constituent components remain unchanged. We believe that this method is also applicable to winter monsoon conditions. The fundamental principle of this method lies in reconstructing vector ocean currents through the utilization of ocean dynamics theory in combination with observations. In contrast to the prevailing direction of summer monsoons, winter monsoons are typically dominated by northeasterly winds. However, other influencing factors on the final inversion of ocean currents remain similar. We believe that both the Ekman theory and Stokes theory remain applicable under winter monsoon conditions, rendering this method suitable for such conditions. Furthermore, as indicated in the Introduction, the calculation methods for these components take into account real-time variations in tides, winds, and radial currents, resulting in the algorithm’s applicability across different seasons without the need for additional adjustments.

5. Conclusions

The Taiwan Strait is characterized by stable ocean currents and hydrodynamic features, including seasonal strong monsoons, powerful tidal waves, perennial northeastward geostrophic currents, and distinctive shallow bank topography. Based on an understanding of the oceanographic conditions in the Taiwan Strait, this study proposes an algorithm for computing the two principal components of surface vector currents: tidal currents and wind-driven currents. The algorithm utilizes radial current observations from a single-station HFSWR to determine the magnitude of the geostrophic currents, thus retrieving the surface vector current.
This algorithm not only reduces costs but also provides accurate surface vector currents comparable to those obtained by dual-site high-frequency radar observations, demonstrating its practical applicability. Although its efficacy has been validated in the context of the Taiwan Strait, additional research is imperative to evaluate its suitability across diverse oceanic regions.

Author Contributions

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

Funding

This work was supported in part by a Project Funded by China Postdoctoral Science Foundation (2023M734315 and 2023M741774), and in part by the National Natural Science Foundation of China, grant number (42305150 and 61771352).

Data Availability Statement

We thank the CMEMS SSH data obtained from https://doi.org/10.48670/moi-00021, accessed on 10 September 2023. The sub-region could be from 115°E to 125°E, from 21°N to 25°N, time span is from 2 April to 6 April 2017, the variable is sea surface height above geoid. The surface wind from the European Centre for Medium-Range Weather Forecasts (ECMWF) hourly reanalysis single level (ERA5) is available at https://doi.org/10.24381/cds.adbb2d47, accessed on 10 September 2023. The sub-region could be from 115°E to 125°E, from 21°N to 25°N, time span is from 2 April to 6 April 2017, the variables are 10 m u and v component of wind. The bathymetry data are from GEBCO_2022, which can be downloaded from https://download.gebco.net, accessed on 10 September 2023. The sub-region could be from 115°E to 125°E, from 21°N to 25°N. All the data used in this study including radial currents, wind, SSH and bathymetry, as well as codes are available and can be accessed via the REDOC catalog at https://github.com/MengyanF/REDOC, accessed on 10 September 2023.

Acknowledgments

We thank Xiaoyan Li for processing the HFSWR data during field experiments.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Barrick, D.E. A coastal radar system for tsunami warning. Remote Sens. Environ. 1979, 8, 353–358. [Google Scholar] [CrossRef]
  2. Kondo, M. Oceanographic investigations of fishing grounds in the East China Sea and the Yellow Sea, 1: Characteristics of the mean temperature and salinity distributions measured at 50 m and near the bottom. Bull. Seikai Reg. Fish. Res. Lab. 1985, 62, 19–66. [Google Scholar]
  3. Qiu, Y.; Li, L.; Chen, C.-T.A.; Guo, X.; Jing, C. Currents in the Taiwan Strait as observed by surface drifters. J. Oceanogr. 2011, 67, 395–404. [Google Scholar] [CrossRef]
  4. Hu, Z.; Qi, Y.; He, X.; Wang, Y.-H.; Wang, D.-P.; Cheng, X.; Liu, X.; Wang, T. Characterizing surface circulation in the Taiwan Strait during NE monsoon from Geostationary Ocean Color Imager. Remote Sens. Environ. 2019, 221, 687–694. [Google Scholar] [CrossRef]
  5. Huang, W.; Wu, S.; Gill, E.; Wen, B.; Hou, J. HF radar wave and wind measurement over the Eastern China Sea. IEEE Trans. Geosci. Remote Sens. 2002, 40, 1950–1955. [Google Scholar] [CrossRef]
  6. Gurgel, K.W.; Antonischki, G.; Essen, H.H.; Schlick, T. Wellen Radar (WERA): A new ground-wave HF radar for ocean remote sensing. Coast. Eng. 1999, 37, 219–234. [Google Scholar] [CrossRef]
  7. Evans, M.W.; Georges, T.M. Coastal Ocean Dynamics Radar (CODAR): NOAA’s Surface Current Mapping System. In Proceedings of the OCEANS’79, San Diego, CA, USA, 17–19 September 1979. [Google Scholar]
  8. Li, M.; Wu, X.; Zhang, L.; Yue, X.; Li, C.; Liu, J. A New Algorithm for Surface Currents Inversion with High-Frequency Over-the-Horizon Radar. IEEE Geosci. Remote Sens. Lett. 2017, 14, 1303–1307. [Google Scholar] [CrossRef]
  9. Li, M.; Zhang, L.; Wu, X.; Yue, X.; Emery, W.J.; Yi, X.; Liu, J.; Yang, G. Ocean Surface Current Extraction Scheme with High-Frequency Distributed Hybrid Sky-Surface Wave Radar System. IEEE Trans. Geosci. Remote Sens. 2018, 56, 4678–4690. [Google Scholar] [CrossRef]
  10. Feng, M.; Ai, W.; Lu, W.; Shan, C.; Ma, S.; Chen, G. Sea surface temperature retrieval based on simulated microwave polarimetric measurements of a one-dimensional synthetic aperture microwave radiometer. Acta Oceanol. Sin. 2021, 40, 122–133. [Google Scholar] [CrossRef]
  11. Qazi, W.A.; Emery, W.J.; Hansen, M.W. Extracting Ocean Surface Currents from SAR: MCC and Doppler Centroid Tracking. In AGU Fall Meeting Abstracts; AGU: Washington DC, USA, 2012. [Google Scholar]
  12. Martin, A.C.H.; Gommenginger, C. Towards wide-swath high-resolution mapping of total ocean surface current vectors from space: Airborne proof-of-concept and validation. Remote Sens. Environ. 2017, 197, 58–71. [Google Scholar] [CrossRef]
  13. Ren, Y.; Li, X.M.; Gao, G.; Busche, T.E. Derivation of Sea Surface Tidal Current from Spaceborne SAR Constellation Data. IEEE Trans. Geosci. Remote Sens. 2017, 55, 3236–3247. [Google Scholar] [CrossRef]
  14. Qazi, W.A.; Emery, W.J.; Fox-Kemper, B. Computing Ocean Surface Currents Over the Coastal California Current System Using 30-Min-Lag Sequential SAR Images. IEEE Trans. Geosci. Remote Sens. 2014, 52, 7559–7580. [Google Scholar] [CrossRef]
  15. Lipa, B.; Barrick, D. Least-squares methods for the extraction of surface currents from CODAR crossed-loop data: Application at ARSLOE. IEEE J. Ocean. Eng. 1983, 8, 226–253. [Google Scholar] [CrossRef]
  16. Liu, L.; Wu, X.; Cheng, F.; Yang, S.; Ke, H. Algorithm for HF radar vector current measurements. J. Oceanogr. 2007, 63, 47–66. [Google Scholar] [CrossRef]
  17. Kim, S.Y.; Terrill, E.J.; Cornuelle, B.D. Mapping surface currents from HF radar radial velocity measurements using optimal interpolation. J. Geophys. Res. Ocean. 2008, 113, C10. [Google Scholar] [CrossRef]
  18. Marmain, J.; Forget, P.; Molcard, A. Characterization of ocean surface current properties from single site HF/VHF radar. Ocean Dyn. 2011, 61, 1967–1979. [Google Scholar] [CrossRef]
  19. Shi, J.; Chen, X.; Li, X.; Chen, P.; Shao, Q. Construction of sea surface current vectors using a single HF radar: A theoretical study. iScience 2022, 25, 105722. [Google Scholar] [CrossRef]
  20. Chang, G.; Li, M.; Zhang, L.; Ji, Y.; Xie, J. Measurements of ocean surface currents using shipborne High-Frequency radar. In Proceedings of the 2014 IEEE Radar Conference, Cincinnati, OH, USA, 19–23 May 2014; pp. 1067–1070. [Google Scholar]
  21. Yongtan, Z.Z.L. Computational Method for the Extraction of Suface-Current by HF Ground-Wave-Radar %J Marine Science Bulletin. Mar. Sci. Bull. 1997, 16, 79–85. [Google Scholar]
  22. Wang, Z.; Xie, J.; Ji, Z.; Quan, T. Remote sensing of surface currents with single shipborne high-frequency surface wave radar. Ocean. Dyn. 2016, 66, 27–39. [Google Scholar] [CrossRef]
  23. Shulman, I.; Paduan, J.D. Assimilation of HF radar-derived radials and total currents in the Monterey Bay area. Deep. Sea Res. Part II 2009, 56, 149–160. [Google Scholar] [CrossRef]
  24. Han, S.; Yang, H.; Xue, W.; Wang, X. The study of single station inverting the sea surface current by HF ground wave radar based on adjoint assimilation technology. J. Ocean. Univ. China 2017, 16, 383–388. [Google Scholar] [CrossRef]
  25. Wang, L.; Pawlowicz, R.; Wu, X.; Yue, X. Wintertime Variability of Currents in the Southwestern Taiwan Strait. J. Geophys. Res. Ocean. 2021, 126, e2020JC016586. [Google Scholar] [CrossRef]
  26. Jan, S.; Chern, C.-S.; Wang, J.; Chao, S.-Y. The anomalous amplification of M2 tide in the Taiwan Strait. Geophys. Res. Lett. 2004, 31. [Google Scholar] [CrossRef]
  27. Zhang, W.-Z.; Shi, F.; Hong, H.-S.; Shang, S.-P.; Kirby, J.T. Tide-surge Interaction Intensified by the Taiwan Strait. J. Geophys. Res. Ocean. 2010, 115. [Google Scholar] [CrossRef]
  28. Lin, S.F.; Tang, T.Y.; Jan, S.; Chen, C.J. Taiwan strait current in winter. Cont. Shelf Res. 2005, 25, 1023–1042. [Google Scholar] [CrossRef]
  29. Wang, M.; Wu, X.; Zhang, L.; Yue, X.; Yi, X.; Xie, X.; Yu, L. Measurement and Analysis of Antenna Pattern for MIMO HF Surface Wave Radar. IEEE Antennas Wirel. Propag. Lett. 2023, 22, 1788–1792. [Google Scholar] [CrossRef]
  30. Wang, L.; Wu, X.; Ai, W. A Scheme for Credibility of Surface Currents Derived from High Frequency Radars. J. Phys. Conf. Ser. 2024, 2718, 012009. [Google Scholar] [CrossRef]
  31. Gong, Z.; Ke, H.; Hou, J. Error analyze of current measurements by hf surface wave radar. J. Wuhan Univ. Nat. Sci. Ed. 2007, 53, 356. [Google Scholar]
  32. Egbert, G.D.; Erofeeva, S.Y. Efficient Inverse Modeling of Barotropic Ocean Tides. J. Atmos. Ocean. Technol. 2002, 19, 183–204. (In English) [Google Scholar] [CrossRef]
  33. Zhao, Z.; Alford, M.H. Source and propagation of internal solitary waves in the northeastern South China Sea. J. Geophys. Res. Ocean. 2006, 111, C11. [Google Scholar] [CrossRef]
  34. Padman, L.; Erofeeva, S. A barotropic inverse tidal model for the Arctic Ocean. Geophys. Res. Lett. 2004, 31. [Google Scholar] [CrossRef]
  35. Jiang, L.; Wang, M. Diurnal Currents in the Bohai Sea Derived from the Korean Geostationary Ocean Color Imager. IEEE Trans. Geosci. Remote Sens. 2017, 55, 1437–1450. [Google Scholar] [CrossRef]
  36. Cheng, K.-S.; Ho, C.-Y.; Teng, J.-H. Wind Characteristics in the Taiwan Strait: A Case Study of the First Offshore Wind Farm in Taiwan. Energies 2020, 13, 6492. [Google Scholar] [CrossRef]
  37. Dong, S.; Gong, Y.; Wang, Z. Long-term variations of wind and wave conditions in the Taiwan Strait. Reg. Stud. Mar. Sci. 2020, 36, 101256. [Google Scholar] [CrossRef]
  38. Price, J.F.; Weller, R.A.; Schudlich, R.R. Wind-Driven Ocean Currents and Ekman Transport. Science 1987, 238, 1534–1538. [Google Scholar] [CrossRef] [PubMed]
  39. Lenn, Y.-D.; Chereskin, T.K. Observations of Ekman Currents in the Southern Ocean. J. Phys. Oceanogr. 2009, 39, 768–779. (In English) [Google Scholar] [CrossRef]
  40. Stokes, G.G. On the Theory of Oscillatory Waves. In Mathematical and Physical Papers; Stokes, G.G., Ed.; Cambridge Library Collection—Mathematics; Cambridge University Press: Cambridge, UK, 2009; Volume 1, pp. 197–229. [Google Scholar]
  41. Chen, Z.; Yan, X.-H.; Jiang, Y. Coastal cape and canyon effects on wind-driven upwelling in northern Taiwan Strait. J. Geophys. Res. Ocean. 2014, 119, 4605–4625. [Google Scholar] [CrossRef]
  42. Wu, J. Wind-Stress coefficients over Sea surface near Neutral Conditions—A Revisit. J. Phys. Oceanogr. 1980, 10, 727–740. (In English) [Google Scholar] [CrossRef]
  43. Wu, J. Wind-stress coefficients over sea surface from breeze to hurricane. J. Geophys. Res. Ocean. 1982, 87, 9704–9706. [Google Scholar] [CrossRef]
  44. Kenyon, K.E. Stokes drift for random gravity waves. J. Geophys. Res. (1896–1977) 1969, 74, 6991–6994. [Google Scholar] [CrossRef]
  45. Breivik, Ø.; Bidlot, J.-R.; Janssen, P.A.E.M. A Stokes drift approximation based on the Phillips spectrum. Ocean. Model. 2016, 100, 49–56. [Google Scholar] [CrossRef]
  46. Hasselmann, D.E.; Dunckel, M.; Ewing, J.A. Directional Wave Spectra Observed during JONSWAP 1973. J. Phys. Oceanogr. 1980, 10, 1264–1280. (In English) [Google Scholar] [CrossRef]
  47. Wu, J. Wind stress and surface roughness at air-sea interface. J. Geophys. Res. (1896–1977) 1969, 74, 444–455. [Google Scholar] [CrossRef]
  48. Wang, S. Characteristics of tidal and residual currents in western Taiwan strait. J. Oceanogr. Taiwan Strait 1989, 8, 207–210. [Google Scholar]
  49. Zhu, D.; Li, L.; Li, Y.; Guo, X. Seasonal variation of surface currents in the southwestern Taiwan Strait observed with HF radar. Chin. Sci. Bull. 2008, 53, 2385–2391. [Google Scholar] [CrossRef]
  50. Li, L.; Guo, X.; Liao, E.; Jiang, Y. Subtidal variability in the Taiwan Strait induced by combined forcing of winter monsoon and topography. Sci. China Earth Sci. 2018, 61, 483–493. [Google Scholar] [CrossRef]
  51. Liao, E.; Oey, L.Y.; Yan, X.-H.; Li, L.; Jiang, Y. The Deflection of the China Coastal Current over the Taiwan Bank in Winter. J. Phys. Oceanogr. 2018, 48, 1433–1450. [Google Scholar] [CrossRef]
  52. Guan, B.; Fang, G. Winter counter-wind currents off the southeastern China coast: A review. J. Oceanogr. 2006, 62, 1–24. [Google Scholar] [CrossRef]
  53. Lin, X.; Yan, X.-H.; Jiang, Y.; Zhang, Z. Performance assessment for an operational ocean model of the Taiwan Strait. Ocean. Model. 2016, 102, 27–44. [Google Scholar] [CrossRef]
  54. Hatayama, T.; Awaji, T.; Akitomo, K. Tidal currents in the indonesian seas and their effect on transport and mixing. J. Geophys. Res. 1996, 101, 12353–12373. [Google Scholar] [CrossRef]
  55. Tan, Z.; Wang, J.; Guo, J.; Liu, C.; Zhang, M.; Ma, S. Improving satellite-retrieved cloud base height with ground-based cloud radar measurements. Adv. Atmos. Sci. 2024, 41, 1890–1909. [Google Scholar] [CrossRef]
Figure 1. (a) Wind rose during the radar experiment. (b) LTWDB−01 drifter. (c) Mean credibility distribution and drifter trajectory of the buoy in the southwestern Taiwan Strait. Dots indicate the start of the trajectory. DS and LH represent Dongshan and Longhai radar stations, respectively.
Figure 1. (a) Wind rose during the radar experiment. (b) LTWDB−01 drifter. (c) Mean credibility distribution and drifter trajectory of the buoy in the southwestern Taiwan Strait. Dots indicate the start of the trajectory. DS and LH represent Dongshan and Longhai radar stations, respectively.
Remotesensing 16 02767 g001
Figure 2. Sea surface height and geostrophic current pattern at UTC 12:00 22 August 2021.
Figure 2. Sea surface height and geostrophic current pattern at UTC 12:00 22 August 2021.
Remotesensing 16 02767 g002
Figure 3. Flowchart of the iterative process to determine η .
Figure 3. Flowchart of the iterative process to determine η .
Remotesensing 16 02767 g003
Figure 4. Retrieval current patterns of wind-driven currents (a), tidal currents (b), and geostrophic currents (c) at 14:20 on 22 August 2021 (UTC+8). The units are cm/s, and the magnitude of speed is shown by color. For clarity, the arrows are not on the same scale.
Figure 4. Retrieval current patterns of wind-driven currents (a), tidal currents (b), and geostrophic currents (c) at 14:20 on 22 August 2021 (UTC+8). The units are cm/s, and the magnitude of speed is shown by color. For clarity, the arrows are not on the same scale.
Remotesensing 16 02767 g004
Figure 5. The spatial distribution of credibility exceeding 0.5 is illustrated at 14:20 on 22 August 2021 (UTC+8). Contour lines representing credibility equal to 0.9 are also marked in black curves.
Figure 5. The spatial distribution of credibility exceeding 0.5 is illustrated at 14:20 on 22 August 2021 (UTC+8). Contour lines representing credibility equal to 0.9 are also marked in black curves.
Remotesensing 16 02767 g005
Figure 6. Comparison between HF and RE currents at 14:20 on 22 August 2021 (UTC+8). The blue arrows represent the vector currents (HF) measured by high-frequency radars; the red arrows are the retrieval results (RE).
Figure 6. Comparison between HF and RE currents at 14:20 on 22 August 2021 (UTC+8). The blue arrows represent the vector currents (HF) measured by high-frequency radars; the red arrows are the retrieval results (RE).
Remotesensing 16 02767 g006
Figure 7. Retrieval current patterns of wind-driven currents (a), tidal currents (b), and geostrophic currents (c) at 09:00 on 26 August 2021 (UTC+8). The units are cm/s, and the magnitude of speed is shown by color. For clarity, the arrows are not on the same scale.
Figure 7. Retrieval current patterns of wind-driven currents (a), tidal currents (b), and geostrophic currents (c) at 09:00 on 26 August 2021 (UTC+8). The units are cm/s, and the magnitude of speed is shown by color. For clarity, the arrows are not on the same scale.
Remotesensing 16 02767 g007
Figure 8. The spatial distribution of credibility exceeding 0.5 is illustrated at 09:00 on 26 August 2021 (UTC+8). Contour lines representing credibility equal to 0.9 are also marked in black curves.
Figure 8. The spatial distribution of credibility exceeding 0.5 is illustrated at 09:00 on 26 August 2021 (UTC+8). Contour lines representing credibility equal to 0.9 are also marked in black curves.
Remotesensing 16 02767 g008
Figure 9. Comparison between HF and RE currents at 09:00 on 26 August 2021 (UTC+8). The blue arrows represent the vector currents (HF) measured by high-frequency radars; the red arrows are the retrieval results (RE).
Figure 9. Comparison between HF and RE currents at 09:00 on 26 August 2021 (UTC+8). The blue arrows represent the vector currents (HF) measured by high-frequency radars; the red arrows are the retrieval results (RE).
Remotesensing 16 02767 g009
Figure 10. Comparison between HF, Re, and drifter currents. (a) U component for Drifter1, (b) V component for Drifter1, (c) U component for Drifter2, and (d) V component for Drifter2.
Figure 10. Comparison between HF, Re, and drifter currents. (a) U component for Drifter1, (b) V component for Drifter1, (c) U component for Drifter2, and (d) V component for Drifter2.
Remotesensing 16 02767 g010
Figure 11. Scatters of the HF, RE, and drifter currents with correlation coefficients are shown in each subplot. (ac) U component for Drifter1, (df) V component for Drifter1, (gi) U component for Drifter2, and (jl) V component for Drifter2.
Figure 11. Scatters of the HF, RE, and drifter currents with correlation coefficients are shown in each subplot. (ac) U component for Drifter1, (df) V component for Drifter1, (gi) U component for Drifter2, and (jl) V component for Drifter2.
Remotesensing 16 02767 g011
Table 1. Information of drifters.
Table 1. Information of drifters.
DriftersStart TimeEnd TimePositioning
Accuracy
117:16, 22 August 202116:25, 26 August 2021BeiDou II, <10 m
210:57, 22 August 202102:23, 27 August 2021BeiDou II, <10 m
Table 2. The temporal and spatial resolutions of the data used in this study.
Table 2. The temporal and spatial resolutions of the data used in this study.
Current ComponentTemporal
Resolution
Spatial
Resolution
Data
Source
Tidal currents1 h1/30°OTIS (TPXO8.0-atlas)
Wind-driven currents1 h1/4°ECMWF (ERA5)
Geostrophic currents24 h1/12°CMEMS (GLORYS12V1)
Radial currents10 min5 kmHF radar
Table 3. RMSE of the comparison, which listed before and after “/” correspond to Drifter 1 and 2.
Table 3. RMSE of the comparison, which listed before and after “/” correspond to Drifter 1 and 2.
RMSE_U (cm/s)RMSE_V (cm/s)
HF11.9/16.513.1/22.2
RE12.8/18.513.9/24.1
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, L.; Feng, M.; Ai, W.; Wu, X.; Zhao, X.; Hu, S. Surface Vector Current Retrieval by Single-Station High-Frequency Surface Wave Radar Based on Ocean Dynamics in the Taiwan Strait. Remote Sens. 2024, 16, 2767. https://doi.org/10.3390/rs16152767

AMA Style

Wang L, Feng M, Ai W, Wu X, Zhao X, Hu S. Surface Vector Current Retrieval by Single-Station High-Frequency Surface Wave Radar Based on Ocean Dynamics in the Taiwan Strait. Remote Sensing. 2024; 16(15):2767. https://doi.org/10.3390/rs16152767

Chicago/Turabian Style

Wang, Li, Mengyan Feng, Weihua Ai, Xiongbin Wu, Xianbin Zhao, and Shensen Hu. 2024. "Surface Vector Current Retrieval by Single-Station High-Frequency Surface Wave Radar Based on Ocean Dynamics in the Taiwan Strait" Remote Sensing 16, no. 15: 2767. https://doi.org/10.3390/rs16152767

APA Style

Wang, L., Feng, M., Ai, W., Wu, X., Zhao, X., & Hu, S. (2024). Surface Vector Current Retrieval by Single-Station High-Frequency Surface Wave Radar Based on Ocean Dynamics in the Taiwan Strait. Remote Sensing, 16(15), 2767. https://doi.org/10.3390/rs16152767

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