Feasibility of Oil Slick Detection Using BeiDou-R Coastal Simulation

Oil spills, which can cause severe immediate and long-term harm to marine ecological environments for decades after the initial accident, require rapid and accuratemonitoring. Currently, optical and radar satellite images are used tomonitor oil spills; however, remote sensing generally needs a long revisit period. Global Navigation Satellite System reflected signals (GNSS-R) can provide all-weather and all-day ocean monitoring and is therefore more suitable for oil spill monitoring. To assess the feasibility of the BeiDou Navigation Satellite System reflected signals (BeiDou-R) in detecting oil slicks, a BeiDou-R coastal simulated experiment is performed in this study on the oil slick distribution of an oil pipeline explosion accident. We set up an observation point and selected observation satellites, and a delay-Doppler map (DDM) of an oil-slicked sea surface under coastal scenarios was created by combining the mean-square slope (MSS) model for oil-slicked/clean surfaces and the Zavorotny–Voronovich (Z–V) scattering model. DDM simulation of the coastal scenarios effectively represents the scattering coefficient distribution of the presence of an oil slick.Theoretical analysis revealed that oil slicks can be detected within a radius of less than 5 km around the specular reflection point (SP) for BeiDou-R coastal simulation.


Introduction
Pollution by oil slicks in open ocean areas and coastal waters, whether accidental or deliberate, is a serious threat to marine environments, the economy, and human life.An oil spill occurs somewhere in the world almost every day of the year.The Deepwater Horizon oil spill accident in the Gulf of Mexico, which began on April 20, 2010, directly affected 176,000 km 2 of the Gulf 's surface [1].In addition, an oil pipeline explosion occurring in the Yellow Sea on November 22, 2013, heavily contaminated more than 3000 m 2 of the sea surface; the direct economic loss of this accident has been calculated as 750 million Yuan [2].To reduce the impact of oil spill disasters on the environment and human health, accurate monitoring of oil spills at sea is of particularly importance.Operational oil spill monitoring is currently accomplished using a combination of optical and radar satellite images [3,4].However, because optical satellites are easily affected by weather, they cannot detect oil slicks under adverse conditions.Radar is used for monitoring under all-weather conditions and all day; however, the observation model is often in side-scan mode, which significantly affects the accuracy of oil spill detection owing to speckle noise in Synthetic Aperture Radar (SAR) images.In contrast, Global Navigation Satellite System reflected signals (GNSS-R) provide reflectivity information of the L-band, which is capable of penetrating cloud cover.In addition, the most important difference between GNSS and traditional remote sensing applications is that GNSS remote sensing is a bistatic technique.The transmitter and the receiver are not collocated, resulting in the transmitted GNSS signal reflecting or scattering in a forward direction; traditional remote sensing applications use more conventional back-scattered signals [5].Moreover, GNSS-R payloads installed on different platforms can receive signals from the navigation satellites in large parts of the ocean and can provide long-term monitoring of oil slicks in all-weather conditions and all day.
The powerful GNSS-R was originally put forward in 1993 by Martin-Neira [6] as a new technique for ocean altimetry.Since then, reflections of GNSS signals from different types of surfaces have been widely studied.Several remote sensing applications presently using GNSS-R technology include  ocean altimetry [7][8][9][10], ocean wind [11][12][13], ocean wave sensing [14][15][16], and sea ice detection [17][18][19].In addition, the possibility of detecting an oil-slicked sea surface using GNSS-R technology simulation has been investigated by Valencia et al. [20,21].This technique is based mainly on scattering coefficient retrieval using a delay-Doppler map (DDM).Li and Huang (2013) applied the methodology of a spatial integration approach (SIA) to sea surface oil slick detection simulation using DDM obtained under a simplified scenario with elevation angle of 90 ∘ [22].This method has also been extended to general scenarios [23]; details of the DDM process simulation have been previously reported [24].
In addition to Global Positioning System (GPS), Global Navigation Satellite System (GNSS), and Galileo, other GNSS instruments currently in development include China's BeiDou constellation [25].The development of the BeiDou system provides a new signal source to monitor ocean.Previous research focuses mainly on oil spill simulation detection using GPS reflected signals (GPS-R) technology.However, GPS uses only Medium Earth Orbit (MEO) satellites.Compared with existing GNSS constellations, the Chinese BeiDou satellite navigation system consists of MEO satellites in addition to Inclined Geosynchronous Orbit (IGSO) and Geostationary Earth Orbit (GEO) satellites.Until March 2015, five GEO, five IGSO, and six MEO BeiDou navigation satellites were in orbit.Moreover, the mixed constellation structure of BeiDou allows persistent monitoring of a given ocean area from multiple angles simultaneously.Currently, very few studies are available on ocean oil slick detection simulation using BeiDou-R.To test the performance of oil slick detection by using the BeiDou-R, the present study uses coastal simulation to detect oil-slicked ocean surfaces by using BeiDou IGSO and MEO satellites.Because BeiDou GEO satellites are almost stationary relative to Earth, they cannot use relative motion to produce Doppler frequency shift of DDM and are not considered in simulation experiments.The coastal simulation in this study is based on an actual oil spill accident that caused marine oil pollution.Based on the polluted area, we set up an observation point and select the corresponding observation satellites to monitor the oil spill conditions.Simulation DDM of the oil-slicked surface is then calculated by using the mean-square slope (MSS) model for oil-slicked/clean surfaces and the Z-V scattering model.

Simulation Approach
2.1.Framework.In this paper, a realistic scenario is simulated to improve the reliability of BeiDou-R for detection oil spill, which is derived from QuickBird multispectral images.Furthermore, this paper gives comparative analysis on the performance difference between BeiDou IGSO and MEO satellites for oil slick detection.
Figure 1 shows the specific DDM simulation process for an oil slick scenario.In the proposed framework, coastal simulation experiments using two BeiDou IGSO and MEO satellites are conducted on the basis of a realistic oil slick scenario.First, according to the geographic location information of the oil spill area, we set up the coastal observation point and selected the corresponding observation satellites.The location of the specular reflection point (SP) and the scattering point needs to be implemented according to the design of the coastal scenario.Under the conditions of observing the scattering coefficient distribution of a realistic oil slick, DDM simulation of the oil-slicked surface was then calculated by combining the MSS model for oil-slicked/clean surfaces and the Z-V scattering model.Finally, DDM simulation had been used to retrieve the scattering coefficient distribution of the presence of an oil slick based on the signal delay and Doppler frequency shift.Oil slicks were detected within a radius coverage area around the SP based on the DDM.These methods are discussed in detail in the following sections.

DDM Theory.
In an actual GNSS-R operational system, GNSS signals need to be received and processed using a bistatic configuration.GNSS-R receiver measurement of the sea-scattered signal is based on scattering from the different scattering units of the sea surface.The differences in delay and Doppler frequency shift correspond to the different scattering units of the sea surface, as shown in Figure 2.
GNSS signals scattered from the ocean are affected by the reflecting surface when it presents oil slicks.In such cases, the changes induced by the oil slick surface can be described by the correlation power waveform of the scattered signal in different delay and Doppler frequency shift values.
The DDM accounts for the correlation power distribution of the sea surface signals scattered as a function of delay and Doppler frequency shift with respect to the SP by using the GNSS-R scattering model proposed by Zavorotny and Voronovich expressed as [26,27] where * indicates two-dimensional convolution and  is the ambiguity function described as where the triangle of the Λ(Δ) function is described as In the application of the GNSS-R signals, the delay and the Doppler at the SP are used as the reference point.Δ is the delay of the scattering point with respect to the SP expressed as where  is the frequency of C/A code (BeiDou  = 2.046 × 10 6 Hz, GPS  = 1.023 × 10 6 Hz),   ( ⃗ ) is the distance from the transmitter to the scattering point,   ( ⃗ ) is the distance from the receiver to the scattering point, and  is the speed of light.Δ is the Doppler frequency shift of the scattering points with respect to the SP expressed as [28] where ⃗   and ⃗   are the velocities of transmitter and receiver, respectively; m is the unit vector of the incident wave; and n is the unit vector of the scattered wave.The Doppler frequency shift of the signal at different scattering points across the surface can be accurately determined by using (4) with respect to the known positions and velocities of the transmitter, receiver, and SP.Before calculating the received power, Σ of each scattering point is determined as where   is the coherent integration time,  is the antenna radiation pattern, and ⃗  is the displacement vector of a scattering point from the SP.The SP can be calculated by using the known locations of the receiver and transmitter; details of the calculation are shown elsewhere [5].The expression of each scattering point of the scattering coefficients over the sea surface is [26 where |R| 2 is the Fresnel reflection coefficient; ⃗  is the scattering vector, which is defined as ⃗  = (n − m) = ⃗  ⊥ +   ẑ; and  is the carrier wave number. pdf (− ⃗  ⊥ /  ) represents the probability density function (PDF) of the given sea surface slope.Soulat presented the PDF of surface slopes for arbitrary wind direction [29].Here, we assume that the wind direction is along the -axis of the Earth-centered - incidence (ECXI) system.The slope PDF can be written as [26] where − ⊥, /  and − ⊥, /  represents the upwind and crosswind ocean slope, respectively;  2  is a crosswind MSS component; and  2   is an upwind MSS component.According to the Cox and Munk MSS model, the presence of oil slicks on the sea surface translates into changes in the surface MSS [30].Accordingly, the relationship among the ocean surface's MSS and the wind speed (WS) can be formulated as To model the effect of the oil slick presence, the oil slick surface's MSS is computed by using the appropriate relationship with the WS as where the second subscripts of MSS, that is,  and , represent clean and oil-slicked surfaces, respectively and  10 is the WS at the height of 10 m from the sea surface.This MSS model was originally proposed for optical wavelengths signals.To apply the model to L-band GNSS-R signals, an empirical modification was proposed by Katzberg et al. [31] in 2006.
In this paper, we adopted this empirical MSS model for the clean and slicked regions in the simulation.Oil slicks on the sea surface cause changes in the scattering coefficient distribution over the observed surface and thus change the corresponding entire DDM.By using this oil slick model, the impacts of the oil slick on the DDM can be distinguished.

Simulation Scenarios
The coastal simulation experiment is based on the oil slick distribution from an oil pipeline explosion that occurred in the district of Qingdao, China.The oil spill distribution data were used as input for the DDM simulation by using QuickBird multispectral images.According to the location of oil spill on the sea surface, we assumed two observation points, R1 and R2 receivers, for monitoring the oil spill.At the R1 observation point, we used the BeiDou IGSO C06 satellite to conduct simulation experiments; this simulation is hereafter referred to as experiment I. Figure 3 shows the Skyplot point of IGSO C06 at experiment I at 08:10 LT on November 27, 2013.At the R2 observation point, we used the BeiDou MEO C11 satellite to conduct simulation experiments; this simulation is hereafter referred to as experiment II. Figure 4 shows the Skyplot point of MEO C11 for experiment II at 23:50 LT on November 27, 2013.Figure 5 shows the analog scene monitoring of oil slicks observed on the sea surface, and the black color in the figure represents the shape of the oil slick; the inverted triangles represent the locations of the receiver as R1 and R2; and star points represent the final SPs of C06 (IGSO) and C11 (MEO) in the two experiments.According to the theory of DDM oil slick detection mentioned in Section 2.2, the relative motions of the transmitter and the receiver can produce a Doppler shift to develop DDM for the present of oil slick distinguish in our simulation.The GEO satellite is almost stationary and is therefore not suitable for coastal DDM technology.Thus, the GEO satellite was not considered in our coastal simulation.
In this work, we did not consider the simulation system noise.The influence of the Earth's curvature on the simulation experiment was also not considered because the height of the antenna is only 15 m.On the basis of nearby weather station data from the same period, we determined that the WS observed surface was 6.6 m/s and that the distribution was homogeneous over the entire observation surface.Specific descriptions of the simulation scenarios are shown in Table 1.

Simulation Results and Discussion
On the basis of the experiment I simulation parameters, we obtained the averaged scattering coefficients of both oilslicked and clean ocean surfaces over the observation surface by using ( 6), (8), and (9). Figure 6(a) shows the averaged scattering coefficients of both oil-slicked and clean ocean surfaces in the IGSO scenario.It should be noted that the scattering coefficient of the oil slick decreased more rapidly than that of the clean ocean surface.The figure clearly shows that the scattering coefficients of the two types of surfaces became identical at scattering points 6 km from the SP, which suggests that the oil slick can be easily detected within a range of about 6 km from the scattering points in our simulation conditions.
On the basis of experiment II simulation parameters, we obtained the averaged scattering coefficients of both oilslicked and clean ocean surfaces over the observation surface by using ( 6), (8), and (9). Figure 6(b) shows the averaged scattering coefficients of both oil-slicked and clean ocean surfaces in the MEO scenario.It appears that the averaged scattering coefficients of the two types of surfaces became the same as the scattering points 8 km from the SP, which suggests that oil slicks can be detected in this range.
As shown in Figure 6, the horizontal axis represents the location of each scattering point relative to the SP over the observation surface, and the vertical axis represents the averaged scattering coefficients of both oil-slicked areas and clean areas.The SP is located at the center of the coordinate system ( = 0,  = 0).
A contamination mask was derived from the QuickBird multispectral images (Figure 3) for our simulation, which was   used as input for the DDM simulation.By combining the design in Figure 6 of the two simulations of oil slick detection scenarios and observed actual oil spill distribution, we obtained the scattering coefficient distributions of each scattering point observation surface including the ocean and the oil slick surface by using ( 6), (8), and ( 9).Figures 7(a) and 7(b) show the corresponding scattering coefficient distributions of the two experiments, IGSO and MEO scenarios, respectively.The oil slick area can be clearly distinguished in Figure 7, because it presents a higher scattering coefficient than that of the clean surface.This is attributed to oil damping of the surface waves, which causes more power to be reflected in the forward direction from this observation surface, although the WS over the entire observation surface is considered uniform.In Figure 7, the horizontal (vertical) axis indicates the location of each scattering point relative to the SP over the observation surface -direction (-direction).In the figure, the red (blue) color indicates the scattering coefficient distribution of the oil slick (ocean) surface over the observation surface.The SP is ( = 0,  = 0).
Both simulation experiments I and II were employed in the DDM based on Figure 7. Figure 8 shows the simulated DDM results for simulation experiment I, which represents the scattering coefficient distribution in the presence of the oil slick shown in Figure 7(a).Figure 9 shows that for experiment II, which represents the scattering coefficient distribution in the presence of the oil slick shown in Figure 7(b).As shown in Figures 8 and 9, the area of the oil spill can be clearly distinguished from the simulated DDMs of the two scenarios, and their ranges are within the theoretical ranges of Figures 6 and 7. Hence, the results of the two simulation experiments indicate that it is feasible to use the BeiDou-R coastal setup to detect oil slicks over the observed surface.In Figures 8  and 9, the horizontal axis indicates the delay in chips and the vertical axis indicates Doppler.The different colors represent the scattered power values of the corresponding signals.
In the coastal setup, the receiver's location is fixed; thus, the satellite's height and velocity are major factors in the DDM map.A comparison of Figures 8 and 9 revealed that the Doppler shift and the output power from MEO were wider and stronger, respectively, than those from IGSO.This occurred because the velocity of MEO was faster than that of IGSO in the observed period, and the height of MEO was lower than that of IGSO.According to the DDM theory introduced in Section 2.2, if the observed surface points from MEO are more than IGSO, oil slicks can be distinguished more easily from MEO in our observation period.

Conclusion
This study used BeiDou-R simulation to prove the capability of oil slick detection under a coastal scenario based on a DDM obtained by combining the MSS model for oil-slicked/ clean surfaces and the Z-V scattering model.The simulated retrieval process was based on an oil pipeline explosion accident that occurred in 2013 in the district of Qingdao, China.We conducted two simulation experiments of oil slicks by using the BeiDou IGSO C06 and MEO C11 satellites under the same wind conditions.The simulation results showed that the oil-slicked area could be clearly distinguished by using DDM map technology from BeiDou MEO and IGSO satellites with the coastal setup.
By analysis and comparison, the oil slick was detected more easily by MEO, and oil slicks were clearly distinguishable within a coverage area radius of less than 5 km around the SP in the observation period.Because the receiver's location is fixed in the coastal setup, the satellite's height and velocity affect the DDM results and the oil slick inversion effect.Therefore, a lower satellite height and faster satellite velocity result in stronger output power and better inversion.

Figure 2 :
Figure 2: Corresponding relationship between the scattering units and the delay-Doppler units.
points on sea surfaces to the specular point (km) (a) Results of the scattering coefficients from the BeiDou Inclined Geosynchronous Orbit (IGSO) simulation experiment Oil-slicked sea surface Clean sea points on sea surfaces to the specular point (km) (b) Results of the scattering coefficients from the BeiDou Medium Earth Orbit (MEO) simulation experiment

Figure 6 :
Figure 6: Scattering coefficients of both oil-slicked areas and clean areas versus distance to the specular reflection point (SP) based on the two different simulation experiments.

Figure 7 :
Figure 7: Simulation results of scattering coefficient distribution from both scenarios.

Figure 8 :
Figure 8: Simulated DDM results for simulation experiment I for Inclined Geosynchronous Orbit (IGSO) scenario.