Backtracking drifting objects using surface currents from high-frequency (HF) radar technology

In this work, the benefits of high-frequency (HF) radar ocean observation technology for backtracking drifting objects are analysed. The HF radar performance is evaluated by comparison of trajectories between drifter buoys versus numerical simulations using a Lagrangian trajectory model. High-resolution currents measured by a coastal HF radar network combined with atmospheric fields provided by numerical models are used to backtrack the trajectory of two dataset of surface-drifting buoys: group I (with drogue) and group II (without drogue). A methodology based on optimization methods is applied to estimate the uncertainty in the trajectory simulations and to optimize the search area of the backtracked positions. The results show that, to backtrack the trajectory of the buoys in group II, both currents and wind fields were required. However, wind fields could be practically discarded when simulating the trajectories of group I. In this case, the optimal backtracked trajectories were obtained using only HF radar currents as forcing. Based on the radar availability data, two periods ranging between 8 and 10?h were selected to backtrack the buoy trajectories. The root mean squared error (RMSE) was found to be 1.01?km for group I and 0.82?km for group II. Taking into account these values, a search area was calculated using circles of RMSE radii, obtaining 3.2 and 2.11?km² for groups I and II, respectively. These results show the positive contribution of HF radar currents for backtracking drifting objects and demonstrate that these data combined with atmospheric models are of value to perform backtracking analysis of drifting objects.     

SAR observation and model tracking of an oil spill event in coastal waters

Oil spills are a major contributor to marine pollution. The objective of this work is to simulate the oil spill trajectory of oil released from a pipeline leaking in the Gulf of Mexico with the GNOME (General NOAA Operational Modeling Environment) model. The model was developed by NOAA (National Oceanic and Atmospheric Administration) to investigate the effects of different pollutants and environmental conditions on trajectory results. Also, a Texture-Classifying Neural Network Algorithm (TCNNA) was used to delineate ocean oil slicks from synthetic aperture radar (SAR) observations. During the simulation, ocean currents from NCOM (Navy Coastal Ocean Model) outputs and surface wind data measured by an NDBC (National Data Buoy Center) buoy are used to drive the GNOME model. The results show good agreement between the simulated trajectory of the oil spill and synchronous observations from the European ENVISAT ASAR (Advanced Synthetic Aperture Radar) and the Japanese ALOS (Advanced Land Observing Satellite) PALSAR (Phased Array L-band Synthetic Aperture Radar) images. Based on experience with past marine oil spills, about 63.0% of the oil will float and 18.5% of the oil will evaporate and disperse. In addition, the effects from uncertainty of ocean currents and the diffusion coefficient on the trajectory results are also studied.     

Observation-based evaluation of surface wave effects on currents and trajectory forecasts

Knowledge of upper ocean currents is needed for trajectory forecasts and is essential for search and rescue operations and oil spill mitigation. This paper addresses effects of surface waves on ocean currents and drifter trajectories using in situ observations. The data set includes colocated measurements of directional wave spectra from a wave rider buoy, ocean currents measured by acoustic Doppler current profilers (ADCPs), as well as data from two types of tracking buoys that sample the currents at two different depths. The ADCP measures the Eulerian current at one point, as modelled by an ocean general circulation model, while the tracking buoys are advected by the Lagrangian current that includes the wave-induced Stokes drift. Based on our observations, we assess the importance of two different wave effects: (a) forcing of the ocean current by wave-induced surface fluxes and the Coriolis–Stokes force, and (b) advection of surface drifters by wave motion, that is the Stokes drift. Recent theoretical developments provide a framework for including these wave effects in ocean model systems. The order of magnitude of the Stokes drift is the same as the Eulerian current judging from the available data. The wave-induced momentum and turbulent kinetic energy fluxes are estimated and shown to be significant. Similarly, the wave-induced Coriolis–Stokes force is significant over time scales related to the inertial period. Surface drifter trajectories were analysed and could be reproduced using the observations of currents, waves and wind. Waves were found to have a significant contribution to the trajectories, and we conclude that adding wave effects in ocean model systems is likely to increase predictability of surface drifter trajectories. The relative importance of the Stokes drift was twice as large as the direct wind drag for the used surface drifter.     

Observation-based evaluation of surface wave effects on currents and trajectory forecasts

Knowledge of upper ocean currents is needed for trajectory forecasts and is essential for search and rescue operations and oil spill mitigation. This paper addresses effects of surface waves on ocean currents and drifter trajectories using in situ observations. The data set includes colocated measurements of directional wave spectra from a wave rider buoy, ocean currents measured by acoustic Doppler current profilers (ADCPs), as well as data from two types of tracking buoys that sample the currents at two different depths. The ADCP measures the Eulerian current at one point, as modelled by an ocean general circulation model, while the tracking buoys are advected by the Lagrangian current that includes the wave-induced Stokes drift. Based on our observations, we assess the importance of two different wave effects: (a) forcing of the ocean current by wave-induced surface fluxes and the Coriolis–Stokes force, and (b) advection of surface drifters by wave motion, that is the Stokes drift. Recent theoretical developments provide a framework for including these wave effects in ocean model systems. The order of magnitude of the Stokes drift is the same as the Eulerian current judging from the available data. The wave-induced momentum and turbulent kinetic energy fluxes are estimated and shown to be significant. Similarly, the wave-induced Coriolis–Stokes force is significant over time scales related to the inertial period. Surface drifter trajectories were analysed and could be reproduced using the observations of currents, waves and wind. Waves were found to have a significant contribution to the trajectories, and we conclude that adding wave effects in ocean model systems is likely to increase predictability of surface drifter trajectories. The relative importance of the Stokes drift was twice as large as the direct wind drag for the used surface drifter.

     

Evaluation of two algorithms for a network of coastal HF radars in the Mid-Atlantic Bight

The National High Frequency (HF) Surface Current Mapping Radar Network is being developed as a backbone system within the U.S. Integrated Ocean Observing System. This paper focuses on the application of HF radar-derived surface current maps to U.S. Coast Guard Search and Rescue operations along the Mid-Atlantic coast of the USA. In that context, we evaluated two algorithms used to combine maps of radial currents into a single map of total vector currents. In situ data provided by seven drifter deployments and four bottom-mounted current meters were used to (1) evaluate the well-established unweighted least squares (UWLS) and the more recently adapted optimal interpolation (OI) algorithms and (2) quantify the sensitivity of the OI algorithm to varying decorrelation scales and error thresholds. Results with both algorithms were shown to depend on the location within the HF radar data footprint. The comparisons near the center of the HF radar coverage showed no significant difference between the two algorithms. The most significant distinction between the two was seen in the drifter trajectories. With these simulations, the weighting of radial velocities by distance in the OI implementation was very effective at reducing both the distance between the actual drifter and the cluster of simulated particles as well as the scale of the search area that encompasses them. In this study, the OI further reduced the already improved UWLS-based search areas by an additional factor of 2. The results also indicated that the OI output was relatively insensitive to the varying decorrelation scales and error thresholds tested.     

Analysis of the reliability of a statistical oil spill response model

A statistical oil spill response model is developed and validated by means of actual oil slick observations reported during the Prestige accident and trajectories of drifter buoys. The model is based on the analysis of a database of hypothetical oil spill scenarios simulated by means of a Lagrangian transport model. To carry out the simulations, a re-analysis database consisting of 44-year hindcast dataset of wind and waves and climatologic daily mean surface currents is used. The number of scenarios required to obtain statistically reliable results is investigated, finding that 200 scenarios provide an optimal balance between the accuracy of the results and the computational effort. The reliability of the model was analyzed by comparing the actual data with the numerical results. The agreement found between actual and numerical data shows that the developed statistical oil spill model is a valuable tool to support spill response planning.     

A short-term predictive system for surface currents from a rapidly deployed coastal HF radar network

In order to address the need for surface trajectory forecasts following deployment of coastal HF radar systems during emergency-response situations (e.g., search and rescue, oil spill), a short-term predictive system (STPS) based on only a few hours data background is presented. First, open-modal analysis (OMA) coefficients are fitted to 1-D surface currents from all available radar stations at each time interval. OMA has the effect of applying a spatial low-pass filter to the data, fills gaps, and can extend coverage to areas where radial vectors are available from a single radar only. Then, a set of temporal modes is fitted to the time series of OMA coefficients, typically over a short 12-h trailing period. These modes include tidal and inertial harmonics, as well as constant and linear trends. This temporal model is the STPS basis for producing up to a 12-h current vector forecast from which a trajectory forecast can be derived. We show results of this method applied to data gathered during the September 2010 rapid-response demonstration in northern Norway. Forecasted coefficients, currents, and trajectories are compared with the same measured quantities, and statistics of skill are assessed employing 16 24-h data sets. Forecasted and measured kinetic variances of the OMA coefficients typically agreed to within 10–15%. In one case where errors were larger, strong wind changes are suspected and examined as the cause. Sudden wind variability is not included properly within the STPS attack we presently employ and will be a subject for future improvement.     

Assessing the fidelity of surface currents from a coastal ocean model and HF radar using drifting buoys in the Middle Atlantic Bight

The rapid expansion of urbanization along the world’s coastal areas requires a more comprehensive and accurate understanding of the coastal ocean. Over the past several decades, numerical ocean circulation models have tried to provide such insight, based on our developing understanding of physical ocean processes. The systematic establishment of coastal ocean observation systems adopting cutting-edge technology, such as high frequency (HF) radar, satellite sensing, and gliders, has put such ocean model predictions to the test, by providing comprehensive observational datasets for the validation of numerical model forecasts. The New York Harbor Observing and Prediction System (NYHOPS) is a comprehensive system for understanding coastal ocean processes on the continental shelf waters of New York and New Jersey. To increase confidence in the system’s ocean circulation predictions in that area, a detailed validation exercise was carried out using HF radar and Lagrangian drifter-derived surface currents from three drifters obtained between March and October 2010. During that period, the root mean square (RMS) differences of both the east–west and north–south currents between NYHOPS and HF radar were approximately 15 cm s^?1. Harmonic analysis of NYHOPS and HF radar surface currents shows similar tidal ellipse parameters for the dominant M₂ tide, with a mean difference of 2.4 cm s^?1 in the semi-major axis and 1.4 cm s^?1 in the semi-minor axis and 3° in orientation and 10° in phase. Surface currents derived independently from drifters along their trajectories showed that NYHOPS and HF radar yielded similarly accurate results. RMS errors when compared to currents derived along the trajectory of the three drifters were approximately 10 cm s^?1. Overall, the analysis suggests that NYHOPS and HF radar had similar skill in estimating the currents over the continental shelf waters of the Middle Atlantic Bight during this time period. An ensemble-based set of particle tracking simulations using one drifter which was tracked for 11 days showed that the ensemble mean separation generally increases with time in a linear fashion. The separation distance is not dominated by high frequency or short spatial scale wavelengths suggesting that both the NYHOPS and HF radar currents are representing tidal and inertial time scales correctly and resolving some of the smaller scale eddies. The growing ensemble mean separation distance is dominated by errors in the mean flow causing the drifters to slowly diverge from their observed positions. The separation distance for both HF radar and NYHOPS stays below 30 km after 5 days, and the two technologies have similar tracking skill at the 95 % level. For comparison, the ensemble mean distance of a drifter from its initial release location (persistence assumption) is estimated to be greater than 70 km in 5 days.     

A high-resolution operational forecast system for oil spill response in Belfast Lough

This paper presents a high-resolution operational forecast system for providing support to oil spill response in Belfast Lough. The system comprises an operational oceanographic module coupled to an oil spill forecast module that is integrated in a user-friendly web application. The oceanographic module is based on Delft3D model which uses daily boundary conditions and meteorological forcing obtained from COPERNICUS and from the UK Meteorological Office. Downscaled currents and meteorological forecasts are used to provide short-term oil spill fate and trajectory predictions at local scales. Both components of the system are calibrated and validated with observational data, including ADCP data, sea level, temperature and salinity measurements and drifting buoys released in the study area. The transport model is calibrated using a novel methodology to obtain the model coefficients that optimize the numerical simulations. The results obtained show the good performance of the system and its capability for oil spill forecast.     

Lagrangian circulation and forbidden zone on the West Florida Shelf

This paper presents some recent results of drifters released on the West Florida Shelf during 1996–1997 and compares with the numerical model results of the wind-driven circulation. Using satellite tracked surface drifters during the one year period from February 1996 to February 1997, a drifter free region, called the “forbidden zone”, is found over the southern portion of the West Florida Shelf. This finding is consistent with historical drift bottle data and with a recent numerical model study of the West Florida Shelf circulation response to climatological wind forcing. Direct drifter simulations by numerical model during March 1996 show a good agreement with both the in situ ADCP current observation and drifter observation. Three mechanisms are proposed for the observed Lagrangian features. The primarily dynamic mechanism is the along-shore wind forcing, which induces a coastal jet that tends to leave the coast and the bottom onshore and near surface offshore transports. The second one is the convergent coastal geometry and bottom topography for the southward flow in central shelf near Tampa Bay that enforces the coastal jet and the bottom and near surface transport. The last is a kinematic one, simply due to the short along-shore Lagrangian excursion, driven by the typical synoptic weather systems. Thus near surface shelf waters over the north may not reach the southern coast of the West Florida. Implication is that surface hazard such as oil spill that may occur outside of the southern West Florida shelf may not greatly impact the southern coastal region except Florida Keys. However, the biological and chemical patches over the north that may occur in the water column such as red tides still can easily reach the southern coastal region through the subsurface and bottom waters.     

Stochastic Lagrangian trajectory model for drifting objects in the ocean

The prediction of drifting object trajectories in the ocean is a complex problem plagued with uncertainties. This problem is usually solved simulating the possible trajectories based on wind and advective numerical and/or instrumental data in real time, which are incorporated into Lagrangian trajectory models. However, both data and Lagrangian models are approximations of reality and when comparing trajectory data collected from drifter exercises with respect to Lagrangian models results, they differ considerably. This paper introduces a stochastic Lagrangian trajectory model that allows quantifying the uncertainties related to: (i) the wind and currents numerical and/or instrumental data, and (ii) the Lagrangian trajectory model. These uncertainties are accounted for within the model through random model parameters. The quantification of these uncertainties consists in an estimation problem, where the parameters of the probability distribution functions of the random variables are estimated based on drifter exercise data. Particularly, it is assumed that estimated parameters maximize the likelihood of our model to reproduce the trajectories from the exercise. Once the probability distribution parameters are estimated, they can be used to simulate different trajectories, obtaining location probability density functions at different times. The advantage of this method is that it allows: (i) site specific calibration, and (ii) comparing uncertainties related to different wind and currents predictive tools. The proposed method is applied to data collected during the DRIFTER Project (eranet AMPERA, VI Programa Marco), showing very good predictive skills.     

Temporal variation of chlorophyll a concentration in the coastal waters affected by the Hebei Spirit oil spill in the West Sea of Korea

Time series changes in chlorophyll a concentration before and after the Hebei Spirit oil spill that occurred in December 2007 were analyzed using NCEP wind and SeaWiFS/MODIS ocean color data. Prevailing southwesterly winds and northeast/southwestward tidal currents pushed the oil towards Korea’s West Sea coast of Taean. After the oil spill, daily chlorophyll a concentration decreased about 45-50% compared to the normal condition before the oil spill, and this decrease continued for about two weeks. Monthly mean chlorophyll a concentration in December 2007 was lower compared to the average value for the same month between 1998 and 2007, but, in October and November 2007 before the spill and in January-February 2008 after the spill, the concentration value was higher than average for the same period between 1998 and 2007.     

A relocatable ocean model in support of environmental emergencies

During the Costa Concordia emergency case, regional, subregional, and relocatable ocean models have been used together with the oil spill model, MEDSLIK-II, to provide ocean currents forecasts, possible oil spill scenarios, and drifters trajectories simulations. The models results together with the evaluation of their performances are presented in this paper. In particular, we focused this work on the implementation of the Interactive Relocatable Nested Ocean Model (IRENOM), based on the Harvard Ocean Prediction System (HOPS), for the Costa Concordia emergency and on its validation using drifters released in the area of the accident. It is shown that thanks to the capability of improving easily and quickly its configuration, the IRENOM results are of greater accuracy than the results achieved using regional or subregional model products. The model topography, and to the initialization procedures, and the horizontal resolution are the key model settings to be configured. Furthermore, the IRENOM currents and the MEDSLIK-II simulated trajectories showed to be sensitive to the spatial resolution of the meteorological fields used, providing higher prediction skills with higher resolution wind forcing.     

Simulating the migration of drifters deployed in the Bay of Biscay, during the Prestige crisis

A main conclusion following the oil spill from the Prestige tanker was that improvements in ocean circulation models were necessary; this was in order to predict, more accurately, the trajectories followed by the oil slicks and hence assist in fight against oil pollution operations. In this contribution, the results of the validation of a semi-empirical ocean circulation model, parameterised for the Bay of Biscay and forced with operational oceano-meteorological remote sensing observations, are shown. The model results have been validated with observations from drifting buoys, deployed in the Bay of Biscay during the crisis. The results show that the model explains a relatively large percentage of the current variability. The comparisons between the real and the estimated drifter trajectories indicate that for 3, 5 and 7 day-long trajectories, the drifter position is estimated with errors of approximately 23, 35 and 46km, respectively. The model reproduces relatively well the trajectory followed by the drifter with the shortest period (23 days).     

Numerical reconstruction of trajectory of small-size surface drifter in the Mediterranean sea

In this study, we addressed the effects of wind-induced drift on Lagrangian trajectories of surface sea objects using high-resolution ocean forecast and atmospheric data. Application of stochastic Leeway model for prediction of trajectories drift was considered for the numerical reconstruction of the Elba accident that occurred during the period 21.06.2009–22.06.2009: a person on an inflatable raft was lost in the vicinity of the Elba Island coast; from the initial position, the person on a raft drifted southwards in the open sea and later he was found on a partially deflated raft during rescue operation. For geophysical forcing, we used high-resolution currents from the Mediterranean Forecasting System and atmospheric wind from the European Centre for Medium-Range Weather Forecasts. To investigate the effect of wind on trajectory behavior, numerical simulations were performed using different categories of drifter-like particles similar to a person on an inflatable raft. An algorithm of spatial clustering was used to differentiate the most probable search areas with a high density of particles. Our results showed that the simulation scenarios using particles with characteristics of draft-limited sea drifters provided better prediction of an observed trajectory in terms of the probability density of particles.     

Surface current patterns observed by HF radar: methodology and analysis of currents to the north of the Yaeyama Islands,East China Sea

A new method was developed to compare the classifications of different kind of data maps based on the self-organizing map (SOM) analysis. The surface current maps of the northern coast of Ishigaki and Iriomote Islands, East China Sea, observed by high-frequency (HF) radar were classified by SOM. Winds, sea surface temperatures (SST), and reanalysis data were also classified by SOM. The optimum area for classification was determined objectively by the new method, which relates these patterns to HF radar current patterns. We found two typical surface patterns: the first was that northeastward flows are dominant in the observation area and the second was that a clockwise eddy was also dominant in this area. The southwestward wind pattern was strongly related to the dominant clockwise eddy pattern of the HF radar current field.     

Wind-speed inversion from HF radar first-order backscatter signal

Land-based high-frequency (HF) radars have the unique capability of continuously monitoring ocean surface environments at ranges up to 200 km off the coast. They provide reliable data on ocean surface currents and under slightly stricter conditions can also give information on ocean waves. Although extraction of wind direction is possible, estimation of wind speed poses a challenge. Existing methods estimate wind speed indirectly from the radar derived ocean wave spectrum, which is estimated from the second-order sidebands of the radar Doppler spectrum. The latter is extracted at shorter ranges compared with the first-order signal, thus limiting the method to short distances. Given this limitation, we explore the possibility of deriving wind speed from radar first-order backscatter signal. Two new methods are developed and presented that explore the relationship between wind speed and wave generation at the Bragg frequency matching that of the radar. One of the methods utilizes the absolute energy level of the radar first-order peaks while the second method uses the directional spreading of the wind generated waves at the Bragg frequency. For both methods, artificial neural network analysis is performed to derive the interdependence of the relevant parameters with wind speed. The first method is suitable for application only at single locations where in situ data are available and the network has been trained for while the second method can also be used outside of the training location on any point within the radar coverage area. Both methods require two or more radar sites and information on the radio beam direction. The methods are verified with data collected in Fedje, Norway, and the Ligurian Sea, Italy using beam forming HF WEllen RAdar (WERA) systems operated at 27.68 and 12.5 MHz, respectively. The results show that application of either method requires wind speeds above a minimum value (lower limit). This limit is radar frequency dependent and is 2.5 and 4.0 m/s for 27.68 and 12.5 MHz, respectively. In addition, an upper limit is identified which is caused by wave energy saturation at the Bragg wave frequency. Estimation of this limit took place through an evaluation of a year long database of ocean spectra generated by a numerical model (third generation WAM). It was found to be at 9.0 and 11.0 m/s for 27.68 and 12.5 MHz, respectively. Above this saturation limit, conventional second-order methods have to be applied, which at this range of wind speed no longer suffer from low signal-to-noise ratios. For use in operational systems, a hybrid of first- and second-order methods is recommended.     

Surface current variability in the Keum River Estuary (South Korea) during summer 2002 as observed by high-frequency radar and coastal monitoring buoy

High-frequency (HF) radar observations of surface currents were conducted for 3 months during summer 2002 in the Keum River estuary. A comparison between HF radar-derived currents and directly measured ones form a buoy showed that the regression slope is close to 1 and the correlation coefficient greater than 0.86, with an RMS difference less than 13 cm/s which is less than 17% of the tidal current. This fairly good agreement allows us to use HF radar observation in investigating the surface flow and circulation in this tidal-current-dominant coastal-plume area. To examine the spatial variation in tidal current characteristics, as well as currents associated with non-tidal forcing, the HF radar-derived currents were separated into tidal and sub-tidal frequency currents. The overall pattern of M₂-current ellipse distribution in the study area showed a counterclockwise rotation, with the offshore maximum current direction to the northeast. Eccentricity, the direction of maximum current, and the phase of net motion of the ellipse changed near the estuary mouth and near the gap of the Saemangeum reclamation tide dyke due to the complex coastal geometry and the out-flowing jet during the ebb period.     

A high-resolution real-time forecasting system for predicting the fate of oil spills in the Strait of Bonifacio (western Mediterranean Sea)

The Strait of Bonifacio is a long and narrow area between Corsica and Sardinia. To manage environmental emergencies related to the spill of oil from vessels, an innovative forecasting system was developed. This tool is capable of operationally predicting the dispersion of hydrocarbon spills in the coastal area of the Bonifacio Strait, either from an instantaneous or continuous spill and either in forward or backward mode. Experimental datasets, including ADCP water current measurements and the trajectories of drifter buoys released in the area, were used to evaluate the accuracy of this system. A comparison between the simulation results and experimental data revealed that both the water circulation and the surface transport processes are accurately reproduced by the model. The overall accuracy of the system in reproducing the transport of an oil spill at sea was estimated for both forward and backward prediction mode and in relation to different forecasting time lags.     

The wave climate of Liverpool Bay—observations and modelling

Directional wave measurements have been made in Liverpool Bay by means of wave buoys and acoustic instruments within the footprint of a phased-array high frequency (HF) radar system, which measures currents and waves. Several years of data have now been collected and are supplemented by an 11-year wave model hindcast. Wave parameters have been derived from the various instruments and compared: the directional waverider buoy is taken to provide the ground truth, confirming the good observations obtained from the ADCP; the HF radar wave data have a positive bias, while the model data have a negative bias. The variation of wave climate over various time-scales from seasonal and inter-annual to inter-decadal is examined. Significant wave–current interactions may occur in this area of shallow water and high tidal range and the measurements provide a good test of coupled hydrodynamic-wave models. The waves are mainly fetch-limited: largest events are due to depressions, which track across the UK from SW, generating westerly and WNW winds in the right rear quadrant. Hence, the future extreme wave events will be closely related to future North Atlantic storm tracks. Projections of 50-year return period wave heights differ between different instruments and model datasets. The future wave climate of Liverpool Bay is not expected to change much from the present day; although a slight increase in the severity of the most extreme events is projected, the frequency of extreme wind and wave events in general is slightly reduced. There is evidence for variability on decadal time-scales, with some correlation with the North Atlantic oscillation.