Elimination of the impact of vessels on ocean wave height inversion with X-band wave monitoring radar

Directional wave spectra and integrated wave parameters can be derived from X-band radar sea surface images.A vessel on the sea surface has a significant influence on wave parameter inversions that can be seen as intensive backscatter speckles in X-band wave monitoring radar sea surface images.A novel algorithm to eliminate the interference of vessels in ocean wave height inversions from X-band wave monitoring radar is proposed.This algorithm is based on the characteristics of the interference.The principal components(PCs) of a sea surface image sequence are extracted using empirical orthogonal function(EOF)analysis.The standard deviation of the PCs is then used to identify vessel interference within the image sequence.To mitigate the interference,a suppression method based on a frequency domain geometric model is applied.The algorithm framework has been applied to OSMAR-X,a wave monitoring system developed by Wuhan University,based on nautical X-band radar.Several sea surface images captured on vessels by OSMAR-X are processed using the method proposed in this paper.Inversion schemes are validated by comparisons with data from in situ wave buoys.The root-mean-square error between the significant wave heights(SWH) retrieved from original interference radar images and those measured by the buoy is reduced by 0.25 m.The determinations of surface gravity wave parameters,in particular SWH,confirm the applicability of the proposed method.     

A new modulation transfer function for ocean wave spectra retrieval from X-band marine radar imagery

When imaging ocean surface waves by X-band marine radar,the radar backscatter from the sea surface is modulated by the long surface gravity waves. The modulation transfer function(MTF) comprises tilt,hydrodynamic,and shadowing modulations. A conventional linear MTF was derived using HH-polarized radar observations under conditions of deep water. In this study,we propose a new quadratic polynomial MTF based on VV-polarized radar measurements taken from heterogeneous nearshore wave fields. This new MTF is obtained using a radar-observed image spectrum and in situ buoy-measured wave frequency spectrum. We validate the MTF by comparing peak and mean wave periods retrieved from X-band marine radar image sequences with those measured by the buoy. It is shown that the retrieval accuracies of peak and mean wave periods of the new MTF are better than the conventional MTF. The results also show that the bias and root mean square errors of the peak and mean wave periods of the new MTF are 0.05 and 0.88 s,and 0.32 and 0.53 s,respectively,while those of the conventional MTF are 0.61 and 0.98 s,and 1.39 and 1.48 s,respectively. Moreover,it is also shown that the retrieval results are insensitive to the coefficients in the proposed MTF.     

A joint method to retrieve directional ocean wave spectra from SAR and wave spectrometer data

This paper proposes a joint method to simultaneously retrieve wave spectra at dif ferent scales from spaceborne Synthetic Aperture Radar(SAR) and wave spectrometer data. The method combines the output from the two dif ferent sensors to overcome retrieval limitations that occur in some sea states. The wave spectrometer sensitivity coeffi cient is estimated using an ef fective signifi cant wave height(SWH), which is an average of SAR-derived and wave spectrometer-derived SWH. This averaging extends the area of the sea surface sampled by the nadir beam of the wave spectrometer to improve the accuracy of the estimated sensitivity coeffi cient in inhomogeneous sea states. Wave spectra are then retrieved from SAR data using wave spectrometer-derived spectra as fi rst guess spectra to complement the short waves lost in SAR data retrieval. In addition, the problem of 180° ambiguity in retrieved spectra is overcome using SAR imaginary cross spectra. Simulated data were used to validate the joint method. The simulations demonstrated that retrieved wave parameters, including SWH, peak wave length(PWL), and peak wave direction(PWD), agree well with reference parameters. Collocated data from ENVISAT advanced SAR(ASAR), the airborne wave spectrometer STORM, the PHAROS buoy, and the European Centre for Medium-Range Weather Forecasting(ECMWF) were then used to verify the proposed method. Wave parameters retrieved from STORM and two ASAR images were compared to buoy and ECMWF wave data. Most of the retrieved parameters were comparable to reference parameters. The results of this study show that the proposed joint retrieval method could be a valuable complement to traditional methods used to retrieve directional ocean wave spectra, particularly in inhomogeneous sea states.     

Comparison of retrieving methods of ocean wave periods from satellite altimeter with buoy measurements

For validating the results of retrieved mean wave period, four empirical algorithms established previously are introduced. Based on the data of over five years derived from TOPEX satellite altimeter for the entire East China Sea, ocean wave periods were calculated and statistical comparison among them was performed. The retrieved mean wave period obtained with our new distribution parameters showed better agreement with the wave period TB measured by buoy than that calculated by other three algorithms. The difference between the mean values of and that of TB is 0.16 s and the RMSE (root mean square error) of is the lowest value (0.48).     

The study of single station inverting the sea surface current by HF ground wave radar based on adjoint assimilation technology

This paper introduces the assimilation technology in an ocean dynamics model and discusses the feasibility of inverting the sea surface current in the detection zone by assimilating the sea current radial velocity detected by single station HF ground wave radar in ocean dynamics model. Based on the adjoint assimilation and POM model, the paper successfully inverts the sea surface current through single station HF ground wave radar in the Zhoushan sea area. The single station HF radar inversion results are also compared with the bistatic HF radar composite results and the fixed point measured results by Annderaa current meter. The error analysis shows that acquisition of flow velocity and flow direction data from the single station HF radar based on adjoint assimilation and POM model is viable and the data obtained have a high correlation and consistency with the flow field observed by HF radar.     

Study on numerical simulation and dynamic mechanism of winter-time circulation in the eastern China seas

An MOM2 based 3-dimentional prognostic baroclinic Z-ordinate model was established to study the circulation in eastern China seas, considering the topography, inflow and outflow on the open boundary, wind stress, temperature and salinity exchange on the sea surface. The results were consistent with observation and showed that the Kuroshio intrudes in large scale into the East China Sea continental shelf East China, during which its water is exchanged ceaselessly with outer sea water along Ryukyu Island. The Tsushima Warm Current is derived from several sources, a branch of the Kuroshio, part of the Taiwan Warm Current, and Yellow Sea mixed water coming from the west of Cheju Island. The water from the west of Cheju Island contributes approximately 13% of the Isushima Warm Current total transport through the Korea Strait. The circulation in the Bohai Sea and Yellow Sea is basically cyclonic circulation, and is comprised of coastal currents and the Yellow Sea Warm Current. Besides simulation of the real circulation, numerical experiments were conducted to study the dynamic mechanism. The numerical experiments indicated that wind directly drives the East China Sea and Yellow Sea Coastal Currents, and strengthens the Korea Coastal Current and Yellow Sea Warm Current. In the no wind case, the kinetic energy of the coastal current area and main YSWC area is only 1% of that of the wind case.Numerical experiments also showed that the Tsushima Warm Current is of great importance to the formation of the Korea Coastal Current and Yellow Sea Warm Current.     

A new method for the estimation of oceanic mixed-layer depth using shipboard X-band radar images

Shipboard X-band radar images acquired on 24 June 2009 are used to study nonlinear internal wave characteristics in the northeastern South China Sea. The studied images show three nonlinear internal waves in a packet. A method based on the Radon Transform technique is introduced to calculate internal wave parameters such as the direction of propagation and internal wave velocity from backscatter images. Assuming that the ocean is a two-layer finite depth system, we can derive the mixed-layer depth by applying the internal wave velocity to the mixed-layer depth formula. Results show reasonably good agreement with in-situ thermistor chain and conductivity-temperature-depth data sets.     

Impact of tropical cyclone Matmo on mixed zone of the Yellow and Bohai Seas

The Bohai Sea is a low-lying semi-enclosed sea area that is linked to the Yellow Sea via the Bohai straits(mixed zone). Its of fshore seabed is shallow, which makes it vulnerable to serious marine meteorological disasters associated with the northward passage of Pacific tropical cyclones. Analyses on data of remote sensing and buoy of the mixed zone of the Yellow and Bohai seas indicate that all the wind speed, significant wave height, and salinity(SAL) increased, sea surface temperature decreased, and wind energy density changed considerably during the passage of tropical cyclone Matmo on July 25, 2014. It was found that the SAL inversion layer in the mixed zone of the Yellow and Bohai Seas was caused by the tropical cyclone. Furthermore, it was found that the tropical cyclone transported the northern Yellow Sea cold water mass(NYSCWM) into the mixed zone of the Yellow and Bohai Seas. The NYSCWM has direct influence on both the aquaculture and the ecological environment of the region. Therefore, further research is needed to establish the mechanism behind the formation of the SAL inversion layer in the mixed zone, and to determine the influence of tropical cyclones on the NYSCWM.     

Contrast validation test for retrieval method of high frequency ground wave radar

1 INTRODUCTION Ocean wave and sea wind, are important oce-anic dynamic phenomena having great influence on the development of marine economy, exploitation of marine resources, and location selection, planning and designing, construction and operation of marine projects; so study on measuring methods of ocean wave and sea wind is important. High frequency ground wave radar (HF radar) was a technique developed in the last decades for the detection of oceanic environment. Long radio wave (mu…     

A subsurface intensity index of the cold eddy in the East China Sea

The Cold Eddy in the East China Sea (CEECS) is located about 150 km southwest of Cheju Island. This region experiences a complex current system under the influences of the Yellow Sea Warm Current (YSWC), the Yellow Sea Coastal Current (YSCC), and the northward Kuroshio residual. To identify the strength of the CEECS, a simple subsurface intensity index is developed. Because the CEECS can be traced by temperature minimum, the intensity index is determined by the difference in sea temperature averaged across cores within and surrounding the region. Based on SODA, the CEECS subsurface intensity index time series can be calculated, with which the inter-annual variation of the CEECS is analyzed.     

A New Multigrid 3D-VAR Optimization Method for Bottom Friction Using HF Radar Current Observation

This paper proposes a new method for data assimilation of the surface radial current observed by High Frequency ground wave radar and optimization of the bottom friction coefficient. In this method, the shallow water wave equation is introduced into the cost function of the multigrid three-dimensional variation data assimilation method as the weak constraint term, the surface current and the bottom friction coefficient are defined as the analytical variables, and the high spatiotemporal resolution surface radial flow observed by the high-frequency ground wave radar is used to optimize the surface current and bottom friction coefficient. This method can effectively consider the spatiotemporal correlation of radar data and extract multiscale information from surface radial flow data from long waves to short waves. Introducing the shallow water wave equation into the cost function as a weak constraint condition can adjust both the momentum and mass fields simultaneously to obtain more reasonable analysis information. The optimized bottom friction coefficient is introduced into the regional ocean numerical model to carry out numerical experiments. The test results show that the bottom friction coefficient obtained by this method can effectively improve the accuracy of the numerical simulation of sea surface height in the offshore area and reduce the simulation error.     

Wave steepness retrieved from scatterometer data in a genetic algorithm

Wave steepness is an important characteristic of a high sea state, and is widely applied on wave propagations at ports, ships, offshore platforms, and CO2 circulation in the ocean. Obtaining wave steepness is a difficult task that depends heavily on theoretical research on wavelength distribution and direct observations. Development of remote-sensing techniques provides new opportunities to study wave steepness. At present, two formulas are proposed to estimate wave steepness from QuikSCAT and ERS-1/2 scatterometer data. We found that wave steepness retrieving is not affected by radar band, and polarization method, and that relationship of wave steepness with radar backscattering cross section is similar to that with wind. Therefore, we adopted and modified a genetic algorithm for relating wave steepness with radar backscattering cross section. Results show that the root-mean-square error of the wave steepness retrieved is 0.005 in two cases from ERS-1/2 scatterometer data and from QuikSCAT scatterometer data.     

Impacts of SIS and CICE as Sea Ice Components in BCC_CSM on the Simulation of the Arctic Climate

Sea ice is a sensitive indicator of climate change and an important component of climate system models. The Los Alamos Sea Ice Model 5.0(CICE5.0) was introduced to the Beijing Climate Center Climate System Model(BCC_CSM) as a new alternative to the Sea Ice Simulator(SIS). The principal purpose of this paper is to analyze the impacts of these two sea ice components on simulations of basic Arctic sea ice, atmosphere, and ocean states. Two sets of experiments were conducted with the same configurations except for the sea ice component used, i.e., SIS and CICE. The distributions of sea ice concentration and thickness reproduced by the CICE simulations in both March and September were closer to actual observations than those reproduced by SIS simulations, which presented a very thin sea ice cover in September. Changes in sea ice conditions also brought about corresponding modifications to the atmosphere and ocean circulation. CICE simulations showed higher agreement with the reference datasets than did SIS simulations for surface air temperature, sea level pressure, and sea surface temperature in most parts of the Arctic Ocean. More importantly, compared with simulations with SIS, BCC_CSM with CICE revealed stronger Atlantic meridional overturning circulation(AMOC), which is more consistent with actual observations. Thus, CICE shows better performance than SIS in BCC_ CSM. However, both components demonstrate a number of common weaknesses, such as overestimation of the sea ice cover in winter, especially in the Nordic Sea and the Sea of Okhotsk. Additional studies and improvements are necessary to develop these components further.     

Numerical simulation and validation of ocean waves measured by an Along-Track Interferometric Synthetic Aperture Radar

A new nonlinear integral transform of ocean wave spectra into Along-Track Interferometric Synthetic Aperture Radar （ATI-SAR） image spectra is described. ATI-SAR phase image spectra are calculated for various sea states and radar configurations based on the nonlinear integral transform. The numerical simulations show that the slant range to velocity ratio （R/V）, significant wave height to ocean wavelength ratio （Hi2）, the baseline （2B） and incident angle （0） affect ATI-SAR imaging. The ATI-SAR imaging theory is validated by means of Two X-band, HH-polarized ATI-SAR phase images of ocean waves and eight C-band, HH-polarized ATI-SAR phase image spectra of ocean waves. It is shown that ATI-SAR phase image spectra are in agreement with those calculated by forward mapping in situ directional wave spectra collected simultaneously with available ATI-SAR observations. ATI-SAR spectral correlation coefficients between observed and simulated are greater than 0.6 and are not sensitive to the degree of nonlinearity. However, the ATI-SARoPhase image spectral turns towards the range direction, even if the real ocean wave direction is 30. It is also shown that the ATI-SAR imaging mechanism is significantly affected by the degree of velocity bunching nonlinearity, especially for high values of R/V and H/2.     

Dynamics of anchor last deployment of submersible buoy system

Submersible buoy systems are widely used for oceanographic research,ocean engineering and coastal defense.Severe sea environment has obvious effects on the dynamics of submersible buoy systems.Huge tension can occur and may cause the snap of cables,especially during the deployment period.This paper studies the deployment dynamics of submersible buoy systems with numerical and experimental methods.By applying the lumped mass approach,a three-dimensional multi-body model of submersible buoy system is developed considering the hydrodynamic force,tension force and impact force between components of submersible buoy system and seabed.Numerical integration method is used to solve the differential equations.The simulation output includes tension force,trajectory,profile and dropping location and impact force of submersible buoys.In addition,the deployment experiment of a simplified submersible buoy model was carried out.The profile and different nodes' velocities of the submersible buoy are obtained.By comparing the results of the two methods,it is found that the numerical model well simulates the actual process and conditions of the experiment.The simulation results agree well with the results of the experiment such as gravity anchor's location and velocities of different nodes of the submersible buoy.The study results will help to understand the conditions of submersible buoy's deployment,operation and recovery,and can be used to guide the design and optimization of the system.     

Underwater topography detection of Taiwan Shoal with SAR images

Under suitable conditions of tidal current and wind, underwater topography can be detected by synthetic aperture radar (SAR) indirectly. Underwater topography SAR imaging includes three physical processes: radar ocean surface backscattering, the modulation of sea surface short wave spectrum by the variations in sea surface currents, and the modulation of sea surface currents by the underwater topography. The first process is described usually by Bragg scattering theory because the incident angle of SAR is always between 20°–70°. The second process is described by the action balance equation. The third process is described by an ocean hydrodynamic model. Based on the SAR imaging mechanism for underwater topography, an underwater topography SAR detection model and a simplified method for its calculation are introduced. In the detection model, a two-dimensional hydrodynamic model — the shallow water model is used to describe the motion of tidal current. Due to the difficulty of determining the expression of SAR backscattering cross section in which some terms can not be determined, the backscattering cross section of SAR image used in the underwater topography SAR detection is pro-processed by the simulated SAR image of the coarse-grid water depth to simplify the calculation. Taiwan Shoal, located at the southwest outlet of Taiwan Strait, is selected as an evaluation area for this technique due to the occurrence of hundreds of sand waves. The underwater topography of Taiwan Shoal was detected by two scenes of ERS-2 SAR images which were acquired on 9 January 2000 and 6 June 2004. The detection results are compared with in situ measured water depths for three profiles. The average absolute and relative errors of the best detection result are 2.23 m and 7.5 %, respectively. These show that the detection model and the simplified method introduced in the paper is feasible.     

Tidal current observation in the southern Yellow Sea in the summers of 2001 and 2003

Direct current observations in the Yellow Sea interior are very scarce due to intense fishing and trawling activities. Most previous studies on tides in the area were based on coastal measurements or satellite altimeter sea levels and have not been rigorously compared with direct measurements. In this paper, tidal currents are studied with current profiles from three bottom-moored Sontek Acoustic Doppler Profilers （ADPs） deployed in the southern Yellow Sea in summer of 2001 and 2003. The measured current series were dominated by tidal currents. Maximum velocities are between 40-80 cm/s at the mooring stations. M2 current is the most dominant primary tidal constituent, while MS4 and M4 are the most significant shallow water tides with much smaller amplitudes than the primary tides.     

Barotropic current fluctuations coupled with sea level drawdown in Yellow and Bohai Seas

Sub-tidal barotropic current variations coupled with residual sea level fluctuation in the Bohai and Yellow Seas during wintertime are addressed in this study.The temporal evolution and spatial distribution of current fluctuation are investigated using moored acoustic Doppler current profiler data in a three-dimensional numerical model.It is found that a southward current followed by a northward current occurred in the northern Yellow Sea during the fluctuation,concurrent with a significant outflow followed by inflow through the Bohai Strait.The process is consistent from surface to bottom and is coupled with remarkable residual sea level fluctuation.This quasi three-day fluctuation with amplitude 0.2-0.3 m/s leads to 1 m/1.2 m drawdown in the northern Yellow and Bohai Seas,respectively,strongly influencing water exchange between those seas.Because this a prominent feature in the seas,it is necessary to evaluate its effect on fluctuation during winter in future studies,in particular,the northward current during the recovery phase of sea level in the Bohai and Yellow Seas regarding seasonal variation.     

Wind wave characteristics and engineering environment of the South China Sea

Wave simulation was conducted for the period 1976 to 2005 in the South China Sea (SCS) using the wave model, WAVEWATCH-III. Wave characteristics and engineering environment were studied in the region. The wind input data are from the objective reanalysis wind datasets, which assimilate meteorological data from several sources. Comparisons of significant wave heights between simulation and TOPEX/Poseidon altimeter and buoy data show a good agreement in general. By statistical analysis, the wave characteristics, such as significant wave heights, dominant wave directions, and their seasonal variations, were discussed. The largest significant wave heights are found in winter and the smallest in spring. The annual mean dominant wave direction is northeast (NE) along the southwest (SW)-NE axis, east northeast in the northwest (NW) part of SCS, and north northeast in the southeast (SE) part of SCS. The joint distributions of wave heights and wave periods (directions) were studied. The results show a single peak pattern for joint significant wave heights and periods, and a double peak pattern for joint significant wave heights and mean directions. Furthermore, the main wave extreme parameters and directional extreme values, particularly for the 100-year return period, were also investigated. The main extreme values of significant wave heights are larger in the northern part of SCS than in the southern part, with the maximum value occurring to the southeast of Hainan Island. The direction of large directional extreme H _s values is focus in E in the northern and middle sea areas of SCS, while the direction of those is focus in N in the southeast sea areas of SCS.     

Tidal Effects on the Bottom Thermal Front of North Yellow Sea Cold Water Mass near Zhangzi Island in Summer 2009

Three seabed-mounted TD/CTD chains and two upward-looking acoustic Doppler current profilers(ADCPs) in the southwest of Zhangzi Island are used and a simultaneous cruise observation in the northern North Yellow Sea(NYS) is conducted to study temperature variation in the bottom thermal front zone of the NYS Cold Water Mass(NYSCWM) during the summer of 2009. In the flood-ebb tidal cycles, the bottom temperature decreases(increases) during flood(ebb) tides, which are dominated by the tidal-current induced horizontal advection. The ebb tide-induced temperature increase is larger than the flood tide-induced temperature decrease due to seasonal warming. In the spring-neap tidal cycles, the temperature and the vertical temperature structure show notable fortnightly variation from 16 July to 25 August. The bottom temperature increases from neap to spring tides and decreases from spring to neap. The Richardson number demonstrates strengthened vertical mixing during spring tides but enhanced stratification during neap tides. The spring-neap variation in vertical shear caused by tidal current is the dominant factor that induces the fortnightly variation in vertical mixing and thus bottom temperature.