High Frequency Radar Observations Research Articles

Tsunami Early Warning Using High-Frequency Ocean Radar System in the Kii Channel, Japan

Abstract The tsunami of the 2011 Tohoku earthquake was monitored by the high-frequency (HF) ocean radar system in the Kii Channel, Japan. In this study, we conducted a retrospective analysis of tsunami early warning based on HF radar observations. First, we synthesized the surface current velocity using observation data of two radar land stations: Minato and Saikazaki. After vector synthesis, the time series of current velocity showed high consistency with the results of forward simulation. Then, we used the time series starting from 08:00 (UTC), 11 March 2011, as the input for data assimilation, reconstructing the tsunami wavefield in the Kii Channel. Comparison with actual observations at Kobe showed that our approach can accurately forecast coastal tsunami waveforms at least 50 min prior to tsunami arrival. The normalized prediction error decreased over time, indicating a progressively improving forecast accuracy. This is the first study on tsunami data assimilation using HF radar for the early warning of an actual seismogenic tsunami. We also discussed the limitations of the current HF radar system in the Kii Channel and made suggestions regarding the deployment of future radar land stations.

High-frequency Radar Observations of Convergence (Downwelling) and Water Temperature Variations in Yeongil Bay

High-frequency Radar Observations of Convergence (Downwelling) and Water Temperature Variations in Yeongil Bay

On the significant impact of the 17 March 2015 St. Patrick’s Day geomagnetic storm on the ionosphere over Indian region

The response of the ionosphere in the Indian region to the extreme space weather events of 17–19 March 2015 is studied using Global Positioning System (GPS) Total electron Content (TEC) measurements, Thermosphere Ionosphere Electrodynamics-General Circulation Model (TIE-GCM) simulations and High Frequency (HF) Radar observations. The data from Thiruvananthapuram (8.5° N, 77° E geographic, 0.3° N magnetic latitude) and northern hemispheric midlatitude station Hanle (32.80° N, 79° E geographic, 23.9° N magnetic latitude) have been analyzed for the above event. The TEC over Thiruvananthapuram, a dip equatorial location exhibits a decrease during the post sunset hours on 17 March 2015 and is followed by the generation of F region irregularities at an altitude of 600 km, as observed by HF radar as well as TEC scintillations. The TEC over the mid latitude station of Hanle decreased during the storm recovery phase on 18March 2015 in response to negative storm effects. The extent of decrease in TEC as observed during the recovery phase could not be simulated by TIE-GCM over Hanle. This difference is attributed to the enhanced Joule heating taking place in the polar region during St. Patrick’s Day storm and the associated composition changes. An anomalous heating scenario was introduced and the simulated TIE-GCM results were found to match the observations. The present study clearly brings out the impact of 17 March 2015, St. Patrick’s Day geomagnetic storm on the ionization distribution over equatorial as well as mid latitude stations over the Indian region using observations as well as TIE-GCM simulations and the evolutionary features of ionospheric irregularities over equatorial location using TEC and HF radar observations. The role played by thermosphere ionosphere coupling in modulating the ionization variability over equatorial as well as mid latitude region during periods of extreme geomagnetic forcing is highlighted.

Surface Current Variations and Hydrological Characteristics of the Penghu Channel in the Southeastern Taiwan Strait

Coastal ocean dynamics application radar (CODAR) SeaSonde high-frequency (HF) radars deployed along the coast of Taiwan were used to reveal ocean surface current variations both hourly and through climatological seasons in the Penghu Channel (PHC), southeastern Taiwan Strait (TS), from December 2014 to December 2020. The ocean current in the PHC has a semidiurnal tidal cycle, and the seasonal main flow, wind direction, and wind strength significantly affect the direction and speed of the flow passing through the PHC. The speed of the tidal current in the PHC area can reach more than 1 m/s, and the monthly average flow speed in the PHC is between 0.12 (winter) and 0.24 m/s (summer). Several buoys indicated that the southward flow along the western coast of Taiwan drifted through the PHC in fall and winter. The HF radar observations confirmed the same, implying that this occurred during the strong northeastern monsoon. For a weak northerly wind or even southerly wind, the flow in the PHC can be northward. Different wind directions can affect the speed of the flow passing through the PHC and the branch flow in the northern PHC. The HF radar results are highly consistent with the spatial characteristics of satellite data regarding the sea surface temperature, sea surface salinity, and chlorophyll concentrations; however, there are significant differences from the satellite-derived ocean current.

Optimisation of Parameters in a German Bight Circulation Model by 4DVAR Assimilation of Current and Water Level Observations

Uncertain parameters in a 3D barotropic circulation model of the German Bight are estimated with a variational optimisation approach. Surface current measurements from a high frequency (HF) radar are used in combination with acoustic Doppler current profiler (ADCP) and tide gauge observations as input for a 4DVAR assimilation scheme. The required cost function gradients are estimated using an adjoint model code. The focus of the study is on systematic errors of the model with the control vector including parameters of the bathymetry, bottom roughness, open boundary forcing, meteorological forcing as well as the turbulence model. The model uses the same bathymetry, open boundary forcing, and metereological forcing as the operational model run at the Federal Maritime and Hydrographic Agency (BSH). The baroclinic BSH model is used as a reference to put the performance of the optimised model into perspective. It is shown that the optimised model has better agreement with HF radar data and tide gauge observations both within the fortnight training period and the test period 1 month later. Current profile measurements taken at two platforms indicate that both models have comparable error magnitudes at those locations. The optimised model was also compared with independent drifter data. In this case, drifter simulations based on the BSH model and the respective operational drift model including some surface wave effects were used as a reference. Again, these comparison showed very similar results overall, with some larger errors of the tuned model in very shallow areas, where no observations were used for the tuning and surface wave effects, which are only explicitly considered in the BSH model, play a more important role. The tuned model seems to be slightly more dissipative than the BSH model with more energy entering through the western boundary and less energy leaving toward the north. It also became evident that the 4DVAR cost function minimisation process is complicated by momentum advection, which leads to non-differentiable dependencies of the model with respect to the control vector. It turned out that the omission of momentum advection in the adjoint code still leads to robust estimates of descent directions.

Variational interpolation of high-frequency radar surface currents using DIVAnd

DIVAnd (Data-Interpolating Variational Analysis, in n-dimensions) is a tool to interpolate observations on a regular grid using the variational inverse method. We have extended DIVAnd to include additional dynamic constraints relevant to surface currents, including imposing a zero normal velocity at the coastline, imposing a low horizontal divergence of the surface currents, temporal coherence and simplified dynamics based on the Coriolis force, and the possibility of including a surface pressure gradient. The impact of these constraints is evaluated by cross-validation using the HF (high-frequency) radar surface current observations in the Ibiza Channel from the Balearic Islands Coastal Ocean Observing and Forecasting System (SOCIB). A small fraction of the radial current observations are set aside to validate the velocity reconstruction. The remaining radial currents from the two radar sites are combined to derive total surface currents using DIVAnd and then compared to the cross-validation dataset and to drifter observations. The benefit of the dynamic constraints is shown relative to a variational interpolation without these dynamical constraints. The best results were obtained using the Coriolis force and the surface pressure gradient as a constraint which are able to improve the reconstruction from the Open-boundary Modal Analysis, a quite commonly used method to interpolate HF radar observations, once multiple time instances are considered together.

Annual and Seasonal Surface Circulation Over the Mid‐Atlantic Bight Continental Shelf Derived From a Decade of High Frequency Radar Observations

AbstractA decade (2007–2016) of hourly 6‐km‐resolution maps of the surface currents across the Mid‐Atlantic Bight (MAB) generated by a regional‐scale High Frequency Radar network are used to reveal new insights into the spatial patterns of the annual and seasonal mean surface flows. Across the 10‐year time series, temporal means and interannual and intra‐annual variability are used to quantify the variability of spatial surface current patterns. The 10‐year annual mean surface flows are weaker and mostly cross‐shelf near the coast, increasing in speed and rotating to more alongshore directions near the shelfbreak, and increasing in speed and rotating to flow off‐shelf in the southern MAB. The annual mean surface current pattern is relatively stable year to year compared to the hourly variations within a year. The 10‐year seasonal means exhibit similar current patterns, with winter and summer more cross‐shore while spring and fall transitions are more alongshore. Fall and winter mean speeds are larger and correspond to when mean winds are stronger and cross‐shore. Summer mean currents are weakest and correspond to a time when the mean wind opposes the alongshore flow. Again, intra‐annual variability is much greater than interannual, with the fall season exhibiting the most interseasonal variability in the surface current patterns. The extreme fall seasons of 2009 and 2011 are related to extremes in the wind and river discharge events caused by different persistent synoptic meteorological conditions, resulting in more or less rapid fall transitions from stratified summer to well‐mixed winter conditions.

Statistical Study of the Midlatitude Mesospheric Vertical Winds Observed by the Wuhan and Beijing MST Radars in China

AbstractBased on the long‐term observations of the Wuhan and Beijing Mesosphere‐Stratosphere‐Troposphere (MST) radars, the statistic characteristics of the mesospheric vertical winds at midlatitudes are investigated for the first time. The vertical wind data come from the Wuhan MST radar (29.51°N, 114.13°E) during 2012–2017 and from the Beijing MST radar (39.75°N, 116.96°E) during 2012–2014. The subdaily, seasonal, and annual variations of the mesospheric vertical winds at Wuhan and Beijing are presented. No obvious subdaily and seasonal variations are found and the annual variations of the mean vertical winds are downward in most of time. The mean vertical winds at Wuhan become smaller in summer and turn to upward at some heights, while those of Beijing only become smaller in summer. The difference of summer vertical winds over Wuhan and Beijing is closely related to the gravity wave activity. The Wuhan MST radar is located to the east side of the Qinghai‐Tibet plateau and also in the East Asian monsoon region, thus the mesospheric winds over Wuhan may experience more influence of gravity waves compared with those over Beijing. The possible bias sources in the vertical wind measurement of the two MST radars are carefully analyzed and the presented wind velocities are proved to be credible. The measured velocity magnitude is similar to other Very High Frequency radar observations, but much larger than that of the WACCM simulation. The Stokes drift, multicellular circulation structures, the wave influence on the radar echoes are raised as the potential sources for this discrepancy.

Evaluation of Radial Ocean Surface Currents Derived From Sentinel‐1 IW Doppler Shift Using Coastal Radar and Lagrangian Surface Drifter Observations

AbstractOcean surface radial velocities (RVLs) derived from the Sentinel‐1 A/B Interferomic Wide (IW) mode Doppler frequency shift observations are regularly acquired over the Norwegian coastal zone. These data can be used to complement existing ocean observation systems with high‐resolution (up to 1.5 × 1.5 km) spatial ocean surface current (OSC) maps. In this study, Sentinel‐1 IW Level 2 OSC retrievals were obtained from 2 months (October–November 2017) of raw Doppler shift observations acquired over the Norwegian Coastal Current (NCC). The results were evaluated using coastal high‐frequency radar (HFR) and Lagrangian ocean surface drifter observations. The analysis shows that distinct patterns of the NCC, with range directed currents reaching up to 0.7 m/s, can be detected in the SAR data. The mean bias between Synthetic Aperture Radar (SAR) and HFR observations was −0.08 m/s and the root mean square deviation (RMSD) was 0.25 m/s. In comparison, the agreement between the SAR‐derived OSC and the trajectories from Lagrangian surface drifters showed a mean bias of 0.02 m/s and an RMSD of 0.24 m/s. The accuracy of the SAR OSC retrievals rely on precise wind‐wave bias correction. Hence, the accuracy of the model wind field (speed and direction) is crucial. The sea state contribution must also be taken into account during the bias correction. A typical required accuracy of the OSC velocity is on the order of 0.1 m/s. Therefore, the comparisons demonstrate that the use of SAR for OSC retrieval is promising.

High-Frequency Radar Observations of Surface Circulation Features along the South-Western Australian Coast

A new merged high-frequency radar (HFR) data set collected using SeaSonde and WERA (WEllen RAdar) systems was used to examine the ocean surface circulation at diurnal, seasonal and inter-annual time scales along the south-west coast of Australia (SWWA), between 29°–32° S. Merging was performed after resampling WERA data on the coarser SeaSonde HFR grid and averaging data from the two HFR systems in the area of common overlap. Direct comparisons between WERA and SeaSonde vectors in their overlapping areas provided scalar and vector correlation values in the range Ru = [0.24, 0.76]; Rv = [0.39, 0.83]; ρ = [0.44, 0.75], with mean bias between velocity components in the range [−0.02, 0.28] ms−1, [−0.16, 0.16] ms−1 for the U, V components, respectively. The lower agreement between vectors was obtained in general at the boundaries of the HFR domains, where the combined effects of the bearing errors, geometrical constraints, and the limited angular field of view were predominant. The combined data set allowed for a novel characterization of the dominant features in the region, such as the warmer poleward-flowing Leeuwin Current (LC), the colder Capes Current (CC) and its northward extensions, the presence of sub-mesoscale to mesoscale eddies and their generation and aggregation areas, along with the extent offshore of the inertial-diurnal signal. The contribution of tides was weak within the entire HFR domain (<10% total variance), whilst signatures of significant inertial- and diurnal-period currents were present due to diurnal–inertial resonance. A clear discontinuity in energy and variance distribution occurred at the shelf break, which separates the continental shelf and deeper offshore regions, and defined the core of the LC. Confined between the LC and the coastline, the narrower and colder CC current was a feature during the summer months. Persistent (lifespan greater than 1 day) sub-mesoscale eddies (Rossby number O (1)) were observed at two main regions, north and south of 31.5° S, offshore of the 200 m depth contour. The majority of these eddies had diameters in the range 10–20 km with 50% more counter clockwise rotating (CCW) eddies compared to clockwise (CW) rotating eddies. The northern region was dominated by CCW eddies that were present throughout the year whilst CW eddies were prevalent in the south with lower numbers during the summer months.

Annual and Seasonal Surface Circulation Over the Mid-Atlantic Bight Continental Shelf Derived From a Decade of High Frequency Radar Observations

Annual and Seasonal Surface Circulation Over the Mid-Atlantic Bight Continental Shelf Derived From a Decade of High Frequency Radar Observations

Tide-current-eddy interaction: A seasonal study using high frequency radar observations along the western Bay of Bengal near 16°N

The seasonally reversing boundary currents along the western Bay of Bengal (BoB) provide a background for studying the rich and complex relationship between tidally-driven currents, mean current structure and mesoscale eddies. A year-long high-resolution surface current data set from high-frequency radars (HFR) during 2015 along the Andhra Pradesh coast (near 82°E, 16°N) has been utilized to investigate the interaction of tidal currents with bathymetry and eddies during three seasons: January, March–April and November–December. The seasonal variability shows the reversibility of boundary currents with speeds of ~1.58 (1.24) m s−1 during poleward (equatorward) flow associated with an offshore anticyclonic (inshore cyclonic) eddy in this region. The amplitudes of the tidal constituents from the HFR currents are highly correlated (R > 0.80) with those from the sea level observations from a nearby tide gauge. Among the semi-diurnal and diurnal tidal constituents, M2 (5.25–6.51 cm s−1) and K1 (2.36–4.32 cm s−1) are the dominant two respectively, in all three seasons. The tidal variance of 40–70% is observed throughout the domain, with highest contribution from the M2 tidal currents. The M2 tidal ellipse patterns match well with the mean circulation pattern, i.e., cyclonic (anticyclonic) during November–December (March–April). In the nearshore regions, the tidal currents amplify mainly through shallow-water tidal constituents (M4 and MS4) interacting with bathymetry; whereas in the offshore regions, they are dominated by M2, which co-vary with the mean seasonal circulation pattern. The highest amplitude of M2 tidal currents during November–December is predominantly due to high stratification along the western BoB. In the presence of stratification, the loss of tidal kinetic energy into turbulent processes is reduced, resulting in the amplification of M2 tidal currents. The shallow water constituents are also observed with significant amplitudes of 0.50–3.50 cm s−1 in the nearshore region. The amplitudes of overtides (M3, M4 and S4) associated with an offshore anti-cyclonic eddy are higher during March–April, whereas those for the compound tides (MK3, SK3, MS4, 2SK5, 2MN6 and 2SM6) are associated with an inshore cyclonic eddy dominate during November–December.

Improving Surface Current Resolution Using Direction Finding Algorithms for Multiantenna High-Frequency Radars

AbstractWhile land-based high-frequency (HF) radars are the only instruments capable of resolving both the temporal and spatial variability of surface currents in the coastal ocean, recent high-resolution views suggest that the coastal ocean is more complex than presently deployed radar systems are able to reveal. This work uses a hybrid system, having elements of both phased arrays and direction finding radars, to improve the azimuthal resolution of HF radars. Data from two radars deployed along the U.S. East Coast and configured as 8-antenna grid arrays were used to evaluate potential direction finding and signal, or emitter, detection methods. Direction finding methods such as maximum likelihood estimation generally performed better than the well-known multiple signal classification (MUSIC) method given identical emitter detection methods. However, accurately estimating the number of emitters present in HF radar observations is a challenge. As MUSIC’s direction-of-arrival (DOA) function permits simple empirical tests that dramatically aid the detection process, MUSIC was found to be the superior method in this study. The 8-antenna arrays were able to provide more accurate estimates of MUSIC’s noise subspace than typical 3-antenna systems, eliminating the need for a series of empirical parameters to control MUSIC’s performance. Code developed for this research has been made available in an online repository.

The impact of remote sensing observations on cross-shelf transport estimates from 4D-Var analyses of the Mid-Atlantic Bight

This paper explores the impact of the individual components of a coastal ocean observing system on estimates of the circulation derived from a state-of-the-art analysis and forecast system for the Mid-Atlantic Bight and Gulf of Maine. The foundation of these activities is the Regional Ocean Modeling System 4-dimensional variational (4D-Var) data assimilation platform, which is run in support of the Mid-Atlantic Regional Association Coastal Ocean Observing System as part of the U.S. Integrated Ocean Observing System. The specific focus of this study is on the impact of remote sensing observations from both space- and land-based platforms on estimates of cross-shelf transport in the vicinity of the National Science Foundation Ocean Observatories Initiative Pioneer array. Sea surface temperature (SST) and sea surface height (SSH) were found to have, on average, a similar impact on the transport estimates. However, during a typical 3-day 4D-Var assimilation cycle, approximately two orders of magnitude more observations of SST than SSH are used in the model, and closer analysis shows that each altimeter measurement has approximately 50 times more impact on the transport estimates than an individual SST observation. This highlights the value of altimetry data for ocean state estimation, and the significance of expanding the altimeter constellation. The observations that are most impactful of all are in situ measurements of temperature and salinity, which have typically 3–4 times more impact than an individual SSH datum. A robust geographical distribution of the observation impacts emerges across a range of transport metrics which results from the combined influence of space-time dynamical interpolation and error covariance information within the 4D-Var system. The observation impact calculations suggest that High Frequency (HF) radar estimates of surface currents have relatively little direct influence on cross-shelf transport estimates. However, quantification of the sensitivity of these same estimates to changes in the observing system indicate that HF radar observations indirectly provide important information. This is understood in the current system by appealing to the idea of borrowing strength from the field of statistics in which some observations (satellite remote sensing in the case considered here) can borrow strength from other, seemingly less important observations.

Three-dimensional exchange flows in a semi-enclosed bay: Numerical simulations and high frequency radar observations

In this work we present the exchange flows between a semi-enclosed Bay and the adjacent ocean as measured by a real-time coastal observing system in Baja California, Mexico. Several sources of information as numerical modeling, currents profile data, high frequency radar observations and drifters tracking, all together show an interesting three-dimensional structure to the exchange flows. The northern entrance (relatively shallow, ∼40 m at its deepest portion) showed an estuarine-like exchange pattern, with a net outflow at the shallow parts and inflow at the deepest, suggesting that the frictional effect is more important than the earth's rotation, as typically observed in relatively small bodies of water like estuaries and coastal lagoons. On the other hand, the southern entrance (with a ∼400 m deep canyon) is more complex, with recurrent anticyclonic exchange pattern with inflows at the western limit and outflows at the east. We suggest that the vertical structure of this circulation may be sustained by the deep canyon which lowers the thermocline, but further investigation is needed.

Assessing SARAL/AltiKa Data in the Coastal Zone: Comparisons with HF Radar Observations

We present an initial assessment of SARAL/AltiKa data in the coastal band. The study focuses on the Ibiza Channel where the north-south water exchanges play a key role in controlling the circulation variability in the western Mediterranean. In this area, the track 16 of SARAL/AltiKa intercepts the domain covered by a coastal high-frequency (HF) radar system, which provides surface currents with a range up to 60 km. We evaluate the performance of the SARAL/AltiKa Ssalto/Duacs delayed-time along-track products compared to the HF radar surface velocity fields. SARAL/AltiKa data are retrieved at a distance of only 7 km from the coast, putting in evidence the emerging capabilities of the new altimeter. The derived velocities resolved the general features of the seasonal mesoscale variability with reasonable agreement with HF radar fields (significant correlations of 0.54). However, some discrepancies appear, which might be caused by instrumental hardware radar errors, ageostrophic velocities as well as inaccurate corrections and editing in the altimeter data. Root mean square (rms) differences between the estimated SARAL/AltiKa and the HF radar velocities are about 13 cm/s. These results are consistent with recent studies in other parts of the ocean applying similar approaches to Topex/Poseidon and Jason-1 missions and using coastal altimeter corrections.

Surface Current and Wave Measurement during Cyclone Phailin by High Frequency Radars along the Indian Coast

Cyclone Phailin originated in the east central Bay of Bengal (BoB) and crossed into the Indian mainland after traversing through the BoB. High frequency radar (HFR) operated by the National Institute of Ocean Technology could track the surface currents and high wave activity within its measuring limits. The radar data provide valuable information on the surface dynamics during the cyclone period. The HFR observations compare well with those of wave rider buoy. This opens up opportunities for observing the wave conditions during the cyclonic period over longer distances from the shore. This method is relatively more robust as HFR is less likely to be disrupted due to the passage of cyclones, unlike moored systems.

HF radar observation of tidal current over the Beibu Gulf, South China Sea

In this paper, a new fully-digital high frequency (HF) radar system OSMAR-SD (ocean state measuring and analyzing radar, S-smaller, D-digital) is employed for observation of the tidal current over the Beibu Gulf, South China Sea. This system can operate at two frequencies to achieve wider range and better resolution of measurement. Experimental results show that the radar measurements agree with the tide level data well. The environmental noise level is also analyzed to demonstrate its detrimental effect upon the measurement, the period of this effect is about 24 h. Moreover, it’s the first time that the characteristic of the irregular diurnal tide over Beibu Gulf is clearly measured by HF radar.

Estimation of currents using SARAL/AltiKa in the coastal regions of India

Abstract. The present study explores the possibility of deriving the across track currents along the Indian coastal region from SARAL/AltiKa mission. The across track surface geostrophic currents obtained from along track SARAL altimeter data are directly compared (qualitatively) with high frequency (HF) radar observations of surface currents in the Bay of Bengal. The velocity component from HF radar which is perpendicular to the altimeter tracks is considered. Since the ageostrophic velocity contribution is ignored for the moment, the surface geostrophic currents SARAL compare well only under low wind conditions. Due to high along track resolution of SARAL there are large variations in velocity which are not captured by the HF radar observations. In general, the magnitude and variations in surface currents derived from SARAL altimeter are comparable with HF radar observations.

High-frequency radar observations of the June 2013 US East Coast meteotsunami

We report here on the observation and offline detection of the meteotsunami off the New Jersey coast on June 13, 2013, using coastal radar systems and tide gauges. This work extends the previous observations of tsunamis originating in Japan and Indonesia. The radars observed the meteotsunami 23 km offshore, 47 min before it arrived at the coast. Subsequent observations showed it moving onshore. The neighboring tide gauge height reading provides confirmation of the radar observations near the shore.